Search

Xuan Wei,1,2,* Yu Zhong,3,* Xiaofei Yi,1 Tingting Li,1 Zhougui Ling,2 Moyu Ming,2 Shuang Zhang,1 Zhiyi He1 1Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Guangxi Medical University, Nanning, People’s Republic of China; 2Department of Respiratory and Critical Care Medicine, The Fourth Affiliated Hospital of Guangxi Medical University, Liuzhou, People’s Republic of China; 3Department of Emergency Medicine, The Fourth Affiliated Hospital of Guangxi Medical University, Liuzhou, People’s Republic of China*These authors contributed equally to this workCorrespondence: Zhiyi He, Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Guangxi Medical University, NO. 6 Shuangyong Road, Qingxiu District, Nanning City, Guangxi Zhuang Autonomous Region, People’s Republic of China, 530000, Tel +86 187 7801 7698, Email zhiyi-river@163.comObjective: Chronic obstructive pulmonary disease (COPD) is a chronic respiratory disease with high prevalence, morbidity, and mortality. Chuankezhi (CKZ) injection, a Chinese patent medicine, has been commonly used for treating COPD. This study evaluated the clinical efficacy of CKZ injections in COPD patients and explored potential underlying mechanisms by integrating meta-analysis and network pharmacology.Research Methods: Randomized controlled trials (RCTs) were search in database by Web of Science, Cochrane Library and PubMed as of November 2022 for literature collection, and the Review Manager 5.4 was used to analyze the data. Through the network pharmacology method, the chemical components and their targets, as well as the disease targets were further analyzed.Results: A total of 15 RCTs including 1212 patients were included. The results of meta-analysis showed that CKZ injection can significantly improve the clinical effective rate (RR = 1.25, 95% CI: 1.14 to 1.36), and the clinical advantage was that it can significantly reduced acute exacerbation rate (RR = 0.29, 95% CI: 0.12 to 0.70) and COPD assessment test (CAT) scores (MD =− 4.62, 95% CI:-8.966 to-0.28). A total of 31 chemical compounds and 178 potential targets for CKZ injection were obtained from the online databases. Molecular docking revealed that most key components and targets could form stable structure.Conclusion: This systematic review with meta-analysis and network pharmacology demonstrates that CKZ could effectively improve the clinical efficacy and safety in the treatment of COPD. Such efficacy may be related to an anti-inflammatory effect and immunoregulation of CKZ via multiple components, multiple targets and multiple pathways. Keywords: meta-analysis, network pharmacology, chronic obstructive pulmonary disease, Chuankezhi injection

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Purpose: We previously demonstrated that 12 months of aromatase inhibitor (AI) treatment was not associated with a difference in body compositi on or other markers of cardiometabolic health when compared to controls. Here we report on the pre-planned extension of the study. The pre-specified primary hypothesis was that AI therapy for 24 months would lead to increased visceral adipose tissue (VAT) area when compared to controls. Methods: We completed a 12-month extension to our prospective 12-month cohort study of 52 women commencing AI treatment (median age 64.5 years) and 52 women with breast pathology not requiring endocrine therapy (63.5 years). Our primary outcome of interest was VAT area. Secondary and exploratory outcomes included other measures of body composition, hepatic steatosis, measures of atherosclerosis and vascular reactivity. Using mixed models and the addition of a fourth time point, we increased the number of study observations by 79 and were able to rigorously determine the treatment effect. Results: Among study completers (AI = 39, controls = 40), VAT area was comparable between groups over 24 months, the mean-adjusted difference was −1.54 cm2 (95% CI: −14.9; 11.9, P = 0.79). Both groups demonstrated parallel and continuous increa ses in VAT area over the observation period that did not diverge or change between groups. No statistically significant difference in our secondary and expl oratory outcomes was observed between groups. Conclusions: While these findings provide reassurance that short-to-medium-term exposure to AI therapy is not associated with metabolically adverse changes when compared to controls, risk evolution should be less focussed on the AI-asso ciated effect and more on the general development of cardiovascular risk over time.

The COVID-19 pandemic, driven by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spurred an urgent need for effective therapeutic interventions. The spike glycoprotein of the SARS-CoV-2 is crucial for infiltrating host cells, rendering it a key candidate for drug development. By interacting with the human angiotensin-converting enzyme 2 (ACE2) receptor, the spike initiates the infection of SARS-CoV-2. Linoleate is known to bind the spike glycoprotein, subsequently reducing its interaction with ACE2. However, the detailed mechanisms underlying the protein-ligand interaction remain unclear. In this study, we characterized the pathways of ligand dissociation and the conformational changes associated with the spike glycoprotein by using ligand Gaussian accelerated molecular dynamics (LiGaMD). Our simulations resulted in eight complete ligand dissociation trajectories, unveiling two distinct ligand unbinding pathways. The preference between these two pathways depends on the gate distance between two α-helices in the receptor binding domain (RBD) and the position of the N-linked glycan at N343. Our study also highlights the essential contributions of K417, N121 glycan, and N165 glycan in ligand unbinding, which are equally crucial in enhancing spike-ACE2 binding. We suggest that the presence of the ligand influences the motions of these residues and glycans, consequently reducing accessibility for spike-ACE2 binding. These findings enhance our understanding of ligand dissociation from the spike glycoprotein and offer significant implications for drug design strategies in the battle against COVID-19.

To maintain the lipid asymmetry of the cell envelope in Gram-negative bacteria, the MlaC protein serves as a lipid transfer factor and delivers phospholipids from the outer to the inner membrane. A strategy of antibiotic discovery is to design a proper compound that can tightly bind to the MlaC protein and inhibit the MlaC function. In this study, we performed virtual screening on multiple MlaC structures obtained from molecular dynamics simulations to identify potential MlaC binders. Our results suggested that clorobiocin is a compound that could bind to the MlaC protein. Through the comparison of the bound geometry between clorobiocin and novobiocin, we pointed out that the methyl-pyrrole group of the noviose sugar in clorobiocin forms hydrophobic interactions with amino acids in the phospholipid binding pocket, which allows the compound to bind deep in the active site. This also explains why clorobiocin shows a tighter binding affinity than novobiocin. Our study highlights a practical path of antibiotic development against Gram-negative bacteria.

The authors find that the insulin receptor forms dynamic clusters during insulin signaling and that these clusters become dysfunctional in insulin resistance. This dysfunction is partially rescued by metformin, a drug used to treat type 2 diabetes.

Through adding a harmonic boost potential to smooth the system potential energy surface, Gaussian accelerated molecular dynamics (GaMD) provides enhanced sampling and free energy calculation of biomolecules without the need of predefined reaction coordinates. This work continues to improve the acceleration power and energy reweighting of the GaMD by combining the GaMD with replica exchange algorithms. Two versions of replica exchange GaMD (rex-GaMD) are presented: force constant rex-GaMD and threshold energy rex-GaMD. During simulations of force constant rex-GaMD, the boost potential can be exchanged between replicas of different harmonic force constants with fixed threshold energy. However, the algorithm of threshold energy rex-GaMD tends to switch the threshold energy between lower and upper bounds for generating different levels of boost potential. Testing simulations on three model systems, including the alanine dipeptide, chignolin, and HIV protease, demonstrate that through continuous exchanges of the boost potential, the rex-GaMD simulations not only enhance the conformational transitions of the systems but also narrow down the distribution width of the applied boost potential for accurate energetic reweighting to recover biomolecular free energy profiles.

It is important to determine the binding pathways and mechanisms of ligand molecules to target proteins to effectively design therapeutic drugs. Molecular dynamics (MD) is a promising computational tool that allows us to simulate protein-drug binding at an atomistic level. However, the gap between the time scales of current simulations and those of many drug binding processes has limited the usage of conventional MD, which has been reflected in studies of the HIV protease. Here, we have applied a robust enhanced simulation method, Gaussian accelerated molecular dynamics (GaMD), to sample binding pathways of the XK263 ligand and associated protein conformational changes in the HIV protease. During two of 10 independent GaMD simulations performed over 500-2500 ns, the ligand was observed to successfully bind to the protein active site. Although GaMD-derived free energy profiles were not fully converged because of insufficient sampling of the complex system, the simulations still allowed us to identify relatively low-energy intermediate conformational states during binding of the ligand to the HIV protease. Relative to the X-ray crystal structure, the XK263 ligand reached a minimum root-mean-square deviation (RMSD) of 2.26 Å during 2.5 μs of GaMD simulation. In comparison, the ligand RMSD reached a minimum of only ∼5.73 Å during an earlier 14 μs conventional MD simulation. This work highlights the enhanced sampling power of the GaMD approach and demonstrates its wide applicability to studies of drug-receptor interactions for the HIV protease and by extension many other target proteins.

Malate dehydrogenase (MDH) and citrate synthase (CS) are two pacemaking enzymes involved in the tricarboxylic acid (TCA) cycle. Oxaloacetate (OAA) molecules are the intermediate substrates that are transferred from the MDH to CS to carry out sequential catalysis. It is known that, to achieve a high flux of intermediate transport and reduce the probability of substrate leaking, a MDH-CS metabolon forms to enhance the OAA substrate channeling. In this study, we aim to understand the OAA channeling within possible MDH-CS metabolons that have different structural orientations in their complexes. Three MDH-CS metabolons from native bovine, wild-type porcine, and recombinant sources, published in recent work, were selected to calculate OAA transfer efficiency by Brownian dynamics (BD) simulations and to study, through electrostatic potential calculations, a possible role of charges that drive the substrate channeling. Our results show that an electrostatic channel is formed in the metabolons of native bovine and recombinant porcine enzymes, which guides the oppositely charged OAA molecules passing through the channel and enhances the transfer efficiency. However, the channeling probability in a suggested wild-type porcine metabolon conformation is reduced due to an extended diffusion length between the MDH and CS active sites, implying that the corresponding arrangements of MDH and CS result in the decrease of electrostatic steering between substrates and protein surface and then reduce the substrate transfer efficiency from one active site to another.

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.

The composition of the outer membrane in Gram-negative bacteria is asymmetric, with the lipopolysaccharides found in the outer leaflet and phospholipids in the inner leaflet. The MlaC protein transfers phospholipids from the outer to inner membrane to maintain such lipid asymmetry in the Mla pathway. In this work, we have performed molecular dynamics simulations on apo and phospholipid-bound systems to study the dynamical properties of MlaC. Our simulations show that the phospholipid forms hydrophobic interactions with the protein. Residues surrounding the entrance of the binding site exhibit correlated motions to control the site opening and closing. Lipid binding leads to increase of the binding pocket volume and precludes entry of the water molecules. However, in the absence of the phospholipid, water molecules can freely move in and out of the binding site when the pocket is open. Dehydration occurs when the pocket closes. This study provides dynamic information of the MlaC protein and may facilitate the design of antibiotics against the Mla pathway of Gram-negative bacteria.

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations.

Alpha/beta hydrolase domain-containing 5 (ABHD5), also termed CGI-58, is the key upstream activator of adipose triglyceride lipase (ATGL), which plays an essential role in lipid metabolism and energy storage. Mutations in ABHD5 disrupt lipolysis and are known to cause the Chanarin-Dorfman syndrome. Despite its importance, the structure of ABHD5 remains unknown. In this work, we combine computational and experimental methods to build a 3D structure of ABHD5. Multiple comparative and machine learning-based homology modeling methods are used to obtain possible models of ABHD5. The results from Gaussian accelerated molecular dynamics and experimental data of the apo models and their mutants are used to select the most likely model. Moreover, ensemble docking is performed on representative conformations of ABHD5 to reveal the binding mechanism of ABHD5 and a series of synthetic ligands. Our study suggests that the ABHD5 models created by deep learning-based methods are the best candidate structures for the ABHD5 protein. The mutations of E41, R116, and G328 disturb the hydrogen bonding network with nearby residues and suppress membrane targeting or ATGL activation. The simulations also reveal that the hydrophobic interactions are responsible for binding sulfonyl piperazine ligands to ABHD5. Our work provides fundamental insight into the structure of ABHD5 and its ligand-binding mode, which can be further applied to develop ABHD5 as a therapeutic target for metabolic disease and cancer.

Four new X-ray structures of tryptophan synthase (TS) crystallized with varying numbers of the amphipathic N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F6) molecule are presented. These structures show one of the F6 ligands threaded into the tunnel from the β-site and reveal a distinct hydrophobic region. Over this expanse, the interactions between F6 and the tunnel are primarily nonpolar, while the F6 phosphoryl group fits into a polar pocket of the β-subunit active site. Further examination of TS structures reveals that one portion of the tunnel (T1) binds clusters of water molecules, whereas waters are not observed in the nonpolar F6 binding region of the tunnel (T2). MD simulation of another TS structure with an unobstructed tunnel also indicates the T2 region of the tunnel excludes water, consistent with a dewetted state that presents a significant barrier to the transfer of water into the closed β-site. We conclude that hydrophobic molecules can freely diffuse between the α- and β-sites via the tunnel, while water does not. We propose that exclusion of water serves to inhibit reaction of water with the α-aminoacrylate intermediate to form ammonium ion and pyruvate, a deleterious side reaction in the αβ-catalytic cycle. Finally, while most TS structures show βPhe280 partially blocking the tunnel between the α- and β-sites, new structures show an open tunnel, suggesting the flexibility of the βPhe280 side chain. Flexible docking studies and MD simulations confirm that the dynamic behavior of βPhe280 allows unhindered transfer of indole through the tunnel, therefore excluding a gating role for this residue.

The importance of protonation states and proton transfer in pyridoxal 5'-phosphate (PLP)-chemistry can hardly be overstated. Although experimental approaches to investigate pKa values can provide general guidance for assigning proton locations, only static pictures of the chemical species are available. To obtain the overall protein dynamics for the interpretation of detailed enzyme catalysis in this study, guided by information from solid-state NMR, we performed molecular dynamics (MD) simulations for the PLP-dependent enzyme tryptophan synthase (TRPS), whose catalytic mechanism features multiple quasi-stable intermediates. The primary objective of this work is to elucidate how the position of a single proton on the reacting substrate affects local and global protein dynamics during the catalytic cycle. In general, proteins create a chemical environment and an ensemble of conformational motions to recognize different substrates with different protonations. The study of these interactions in TRPS shows that functional groups on the reacting substrate, such as the phosphoryl group, pyridine nitrogen, phenolic oxygen and carboxyl group, of each PLP-bound intermediate play a crucial role in constructing an appropriate molecular interface with TRPS. In particular, the protonation states of the ionizable groups on the PLP cofactor may enhance or weaken the attractions between the enzyme and substrate. In addition, remodulation of the charge distribution for the intermediates may help generate a suitable environment for chemical reactions. The results of our study enhance knowledge of protonation states for several PLP intermediates and help to elucidate their effects on protein dynamics in the function of TRPS and other PLP-dependent enzymes.

The exocyst complex regulates the last steps of exocytosis, which is essential to organisms across kingdoms. In humans, its dysfunction is correlated with several significant diseases, such as diabetes and cancer progression. Investigation of the dynamic regulation of the evolutionarily conserved exocyst-related processes using mutants in genetically tractable organisms such as Arabidopsis thaliana is limited by the lethality or the severity of phenotypes. We discovered that the small molecule Endosidin2 (ES2) binds to the EXO70 (exocyst component of 70 kDa) subunit of the exocyst complex, resulting in inhibition of exocytosis and endosomal recycling in both plant and human cells and enhancement of plant vacuolar trafficking. An EXO70 protein with a C-terminal truncation results in dominant ES2 resistance, uncovering possible distinct regulatory roles for the N terminus of the protein. This study not only provides a valuable tool in studying exocytosis regulation but also offers a potentially new target for drugs aimed at addressing human disease.

Signaling in cells often involves co-localization of the signaling molecules. Most experimental evidence has shown that intracellular compartmentalization restricts the range of action of the second messenger, 3'-5'-cyclic adenosine monophosphate (cAMP), which is degraded by phosphodiesterases (PDEs). The objective of this study is to understand the details of molecular encounter that may play a role in efficient operation of the cAMP signaling apparatus. The results from electrostatic potential calculations and Brownian dynamics simulations suggest that positive potential of the active site from PDE enhances capture of diffusing cAMP molecules. This electrostatic steering between cAMP and the active site of a PDE plays a major role in the enzyme-substrate encounter, an effect that may be of significance in sequestering cAMP released from a nearby binding site or in attracting more freely diffusing cAMP molecules.

The acid-base chemistry that drives catalysis in pyridoxal-5'-phosphate (PLP)-dependent enzymes has been the subject of intense interest and investigation since the initial identification of PLP's role as a coenzyme in this extensive class of enzymes. It was first proposed over 50 years ago that the initial step in the catalytic cycle is facilitated by a protonated Schiff base form of the holoenzyme in which the linking lysine ε-imine nitrogen, which covalently binds the coenzyme, is protonated. Here we provide the first (15)N NMR chemical shift measurements of such a Schiff base linkage in the resting holoenzyme form, the internal aldimine state of tryptophan synthase. Double-resonance experiments confirm the assignment of the Schiff base nitrogen, and additional (13)C, (15)N, and (31)P chemical shift measurements of sites on the PLP coenzyme allow a detailed model of coenzyme protonation states to be established.

Revealing the processes of ligand-protein associations deepens our understanding of molecular recognition and binding kinetics. Hydrogen bonds (H-bonds) play a crucial role in optimizing ligand-protein interactions and ligand specificity. In addition to the formation of stable H-bonds in the final bound state, the formation of transient H-bonds during binding processes contributes binding kinetics that define a ligand as a fast or slow binder, which also affects drug action. However, the effect of forming the transient H-bonds on the kinetic properties is little understood. Guided by results from coarse-grained Brownian dynamics simulations, we used classical molecular dynamics simulations in an implicit solvent model and accelerated molecular dynamics simulations in explicit waters to show that the position and distribution of the H-bond donor or acceptor of a drug result in switching intermolecular and intramolecular H-bond pairs during ligand recognition processes. We studied two major types of HIV-1 protease ligands: a fast binder, xk263, and a slow binder, ritonavir. The slow association rate in ritonavir can be attributed to increased flexibility of ritonavir, which yields multistep transitions and stepwise entering patterns and the formation and breaking of complex H-bond pairs during the binding process. This model suggests the importance of conversions of spatiotemporal H-bonds during the association of ligands and proteins, which helps in designing inhibitors with preferred binding kinetics.

Protein adsorption on nanoparticles is closely associated with the physicochemical properties of particles, in particular, their surface properties. We synthesized two batches of polyacrylic acid-coated nanoparticles under almost identical conditions except for the heating duration and found differences in the head-group structure of the polyacrylic acid. The structure change was confirmed by NMR and MS. The two batches of particles had varied binding affinities to a selected group of proteins. Computational work confirmed that the head group of the polymer on the surface of a nanoparticle could directly interact with a protein, and small structural changes in the head group were sufficient to result in a significant difference in the free energy of binding. Our results demonstrate that protein adsorption is so sensitive to the surface properties of particles that it can reveal even small variations in the structure of a nanoparticle surface ligand, and should be useful for quick assessment of nanoparticle properties.

The regulation of a series of cellular events requires specific protein-protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.

Abstract Background Phosphopeptide-binding domains mediate many vital cellular processes such as signal transduction and protein recognition. We studied three well-known domains important for signal transduction: BRCT repeats, WW domain and forkhead-associated (FHA) domain. The first two recognize both phosphothreonine (pThr) and phosphoserine (pSer) residues, but FHA has high specificity for pThr residues. Here we used molecular dynamics (MD) simulations to reveal how FHA exclusively chooses pThr and how BRCT and WW recognize both pThr/pSer. The work also investigated the energies and thermodynamic information of intermolecular interactions. Results Simulations carried out included wide-type and mutated systems. Through analysis of MD simulations, we found that the conserved His residue defines dual loops feature of the FHA domain, which creates a small cavity reserved for only the methyl group of pThr. These well-organized loop interactions directly response to the pThr binding selectivity, while single loop (the 2nd phosphobinding site of FHA) or in combination with α-helix (BRCT repeats) or β-sheet (WW domain) fail to differentiate pThr/pSer. Conclusions Understanding the domain pre-organizations constructed by conserved residues and the driving force of domain-phosphopeptide recognition provides structural insight into pThr specific binding, which also helps in engineering proteins and designing peptide inhibitors.

This research is a preliminary study on the preparation of wood plastic composites (WPC) from corn husk fiber (CHF) and recycled polystyrene foam (rPS). The effects of fiber content and alkaline treatment on tensile, thermal, water absorption and morphological properties of the composites were investigated. The rPS/CHF composites were prepared using melt compounding and compression moulding processes. The results showed that an increase of fiber content increased the tensile strength, modulus, and thermal stability (Td50%) of composites. However, the water absorption of composites increased vastly as the fiber content increases. The addition of more fiber also caused an earlier thermal degradation to composites. Alkaline treatment has improved the tensile strength, modulus, thermal stability (Td10%), and also reduced water absorption of rPS/CHF composites. The WPS prepared from rPS and treated CHF shows better tensile and thermal properties with lower water absorption.

Natural fibre reinforced composites (NFRC) possess the tendency of being highly flammable in nature. The current paper investigates the effectiveness of expandable graphite (EG) on the enhancement of the flame retardancy of palm empty fruit bunch (EFB) fibre filled diglycidyl ether of bisphenol A (DGEBA) thermoset epoxy composites. The EFB was treated with 5 wt% sodium hydroxide (NaOH) to improve fibre surface adhesion with the epoxy matrix. In this research, two approaches were employed to introduce EG with concentration ranging from 1 to 5 wt% into the epoxy composite filled with 20 wtn% of EFB. The first approach involved infusion of a pre-mixed EG filled epoxy resin into EFB mat during composite fabrication while the second approach considered was the application of EG filled epoxy resin as a surface coating onto post cured EFB epoxy composite. The flammability, thermal and mechanical properties were studied through a series of experiments, such as the Bunsen burner test, dynamic mechanical analysis (DMA), tensile, flexural and impact tests. The results showed that increasing EG concentration capable of enhancing the flame retardancy and mechanical behaviour in both approaches, with 5 wt% EG demonstrated superior performance in fire resistant and mechanical properties amongst other formulations.

The eco-composites were prepared from corn stalk (CS) and low density polyethylene (LDPE) using Brabender internal mixer. An eco-degradant was used as degradable additive in LDPE/CS eco-composites. This study was focused on the effect of filler loading and eco-degradant on tensile, thermal, morphological and water absorption properties of LDPE/CS eco-composites. The incorporation of CS deteriorated the tensile strength and elongation at break of LDPE/CS ecocomposites, but increased the Young’s modulus of eco-composites. The water absorption of LDPE/CS eco-composites also increased with increases of CS content. However, the addition of CS has contributed to increase of the crystallinity of LDPE matrix which due to nucleating effect of filler. The addition of eco-degradant has increased the tensile strength and elongation at break in average of 17.4% and 78.6%, respectively. However, the Young’s modulus was decreased. The present of eco-degradant also improved the water resistivity and crystallinity of eco-composites. The increase of crystallinity in LDPE/CS ecocomposites causes the increase in melting temperature. The filler dispersion and filler-matrix adhesion of LDPE/CS eco-composites were improved due to addition of eco-degradant which evidenced by scanning electron microscope.

The COVID-19 pandemic, driven by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has spurred an urgent need for effective therapeutic interventions. The spike glycoprotein of the SARS-CoV-2 is crucial for infiltrating host cells, rendering it a key candidate for drug development. By interacting with the human angiotensin-converting enzyme 2 (ACE2) receptor, the spike initiates the infection of SARS-CoV-2. Linoleate is known to bind the spike glycoprotein, subsequently reducing its interaction with ACE2. However, the detailed mechanisms underlying the protein-ligand interaction remain unclear. In this study, we characterized the pathways of ligand dissociation and the conformational changes associated with the spike glycoprotein by using ligand Gaussian accelerated molecular dynamics (LiGaMD). Our simulations resulted in eight complete ligand dissociation trajectories, unveiling two distinct ligand unbinding pathways. The preference between these two pathways depends on the gate distance between two α-helices in the receptor binding domain (RBD) and the position of the N-linked glycan at N343. Our study also highlights the essential contributions of K417, N121 glycan, and N165 glycan in ligand unbinding, which are equally crucial in enhancing spike-ACE2 binding. We suggest that the presence of the ligand influences the motions of these residues and glycans, consequently reducing accessibility for spike-ACE2 binding. These findings enhance our understanding of ligand dissociation from the spike glycoprotein and offer significant implications for drug design strategies in the battle against COVID-19. Author summary: Our work explores the intricate process of ligand dissociation from the spike glycoprotein of SARS-CoV-2, the virus responsible for the ongoing COVID-19 pandemic. The spike glycoprotein plays a crucial role in the virus's ability to infect human cells by interacting with the ACE2 receptor. Our study investigates the molecular interactions involving the spike glycoprotein and linoleate, a small molecule with the potential to hinder viral entry by blocking ACE2 interaction. Using advanced computational simulations, we uncover the dynamic mechanisms governing linoleate dissociation from the spike protein, unveiling two distinct unbinding pathways. Our findings illuminate the structural changes of the spike glycoprotein during the dissociation process. Notably, we highlight the crucial roles of a gate formed by two α-helices and an N-linked glycan at position N343, which control the dissociation process and determine the unbinding pathway. These insights have significant implications for the development of novel antiviral therapies and improved COVID-19 vaccines. By unraveling the molecular details of this interaction, our work contributes to the global efforts aimed at combatting the COVID-19 pandemic. [ABSTRACT FROM AUTHOR]

Inhibition of the protein-protein interaction (PPI) mediated by breast-cancer-gene 1 C-terminal (BRCT) is an attractive strategy to sensitize breast and ovarian cancers to chemotherapeutic agents that induce DNA damage. Such inhibitors could also be used for studies to understand the role of this PPI in DNA damage response. However, design of BRCT inhibitors is challenging because of the inherent flexibility associated with this domain. Several studies identified short phosphopeptides as tight BRCT binders. Here we investigated the thermodynamic properties of 18 phosphopeptides or peptide with phosphate mimic and three compounds with phosphate groups binding to BRCT to understand promiscuous molecular recognition and guide inhibitor design. We performed molecular dynamics (MD) simulations to investigate the interactions between inhibitors and BRCT and their dynamic behavior in the free and bound states. MD simulations revealed the key role of loops in altering the shape and size of the binding site to fit various ligands. The mining minima (M2) method was used for calculating binding free energy to explore the driving forces and the fine balance between configuration entropy loss and enthalpy gain. We designed a rigidified ligand, which showed unfavorable experimental binding affinity due to weakened enthalpy. This was because it lacked the ability to rearrange itself upon binding. Investigation of another phosphate group containing compound, C1, suggested that the entropy loss can be reduced by preventing significant narrowing of the energy well and introducing multiple new compound conformations in the bound states. From our computations, we designed an analog of C1 that introduced new intermolecular interactions to strengthen attractions while maintaining small entropic penalty. This study shows that flexible compounds do not always encounter larger entropy penalty, compared with other more rigid binders, and highlights a new strategy for inhibitor design.