scholarly journals How and why plants and human N-glycans are different: Insight from molecular dynamics into the “glycoblocks” architecture of complex carbohydrates

2020 ◽  
Vol 16 ◽  
pp. 2046-2056 ◽  
Author(s):  
Carl A Fogarty ◽  
Aoife M Harbison ◽  
Amy R Dugdale ◽  
Elisa Fadda
Keyword(s):  
Molecular Dynamics ◽  
3D Structure ◽  
Md Simulations ◽  
Structural Disorder ◽  
Dynamic Features ◽  
Structural Dynamic ◽  
Alternative Description ◽  
Spatial Environment ◽  
3D Architecture ◽  
Complex Carbohydrates

The N-glycosylation is one of the most abundant and diverse post-translational modifications of proteins, implicated in protein folding and structural stability, and mediating interactions with receptors and with the environment. All N-glycans share a common core from which linear or branched arms stem from, with functionalization specific to different species and to the cells’ health and disease state. This diversity generates a rich collection of structures, all diversely able to trigger molecular cascades and to activate pathways, which also include adverse immunogenic responses. These events are inherently linked to the N-glycans’ 3D architecture and dynamics, which remain for the large part unresolved and undetected because of their intrinsic structural disorder. In this work we use molecular dynamics (MD) simulations to provide insight into N-glycans’ 3D structure by analysing the effects of a set of very specific modifications found in plants and invertebrate N-glycans, which are immunogenic in humans. We also compare these structural motifs and combine them with mammalian N-glycan motifs to devise strategies for the control of the N-glycan 3D structure through sequence. Our results suggest that the N-glycans’ architecture can be described in terms of the local spatial environment of groups of monosaccharides. We define these “glycoblocks” as self-contained 3D units, uniquely identified by the nature of the residues they comprise, their linkages and structural/dynamic features. This alternative description of glycans’ 3D architecture can potentially lead to an easier prediction of sequence-to-structure relationships in complex carbohydrates, with important implications in glycoengineering design.

scholarly journals How and why plants and human N-glycans are different: Insight from molecular dynamics into the “glycoblocks” architecture of complex carbohydrates

2020 ◽  
Author(s):  
CA Fogarty ◽  
AM Harbison ◽  
AR Dugdale ◽  
E Fadda
Keyword(s):  
Molecular Dynamics ◽  
3D Structure ◽  
Md Simulations ◽  
Structural Disorder ◽  
Dynamic Features ◽  
Structural Dynamic ◽  
Alternative Description ◽  
Spatial Environment ◽  
3D Architecture ◽  
Complex Carbohydrates

AbstractN-glycosylation is one of the most abundant and diverse post-translational modifications of proteins, implicated in protein folding and structural stability, and mediating interactions with receptors and with the environment. All N-glycans share a common core from which linear or branched arms stem from, with functionalization specific to different species and to the cells’ health and disease state. This diversity generates a rich collection of structures, all diversely able to trigger molecular cascades and to activate pathways, which also include adverse immunogenic responses. These events are inherently linked to the N-glycans 3D architecture and dynamics, which remain for the large part unresolved and undetected because of their intrinsic structural disorder. In this work we use molecular dynamics (MD) simulations to provide insight into N-glycans 3D structure by analysing the effects of a set of very specific modifications found in plants and invertebrate N-glycans, which are immunogenic in humans. We also compare these structural motifs and combine them with mammalian N-glycans motifs to devise strategies for the control of the N-glycan 3D structure through sequence. Our results suggest that the N-glycans architecture can be described in terms of the local spatial environment of groups of monosaccharides. We define these “glycoblocks” as self-contained 3D units, uniquely identified by the nature of the residues they comprise, their linkages and structural/dynamic features. This alternative description of glycans 3D architecture can potentially lead to an easier prediction of sequence-to-structure relationships in complex carbohydrates, with important implications in glycoengineering design.

Molecular Dynamics Simulations of HLA-CW4-Β2M-KIR2DL1 Protein and Homology Modeling of a Complex Associated with Psoriasis Disease (HLA-CW6-Β 2M-KIR2DS1)

2021 ◽  
Author(s):  
Mansour H Almatarneh ◽  
Ahmad M Alqaisi ◽  
Enas K Ibrahim ◽  
Ghada G Kayed ◽  
Joshua W Hollett
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Complex Structure ◽  
3D Structure ◽  
Md Simulations ◽  
Salt Bridges ◽  
Small Peptide ◽  
Β2 Microglobulin ◽  
Dynamics Simulations ◽  
The Right

Molecular dynamics (MD) simulation was used to study the interactions of two immune proteins of HLA-Cw4-β2m-KIR2DL1 complex with small peptide QYDDAVYKL (nine amino acids) in an aqueous solution. This study aims to gain a detailed information about the conformational changes and the dynamics of the complex. The right parameters and force field for performing the MD simulations that was needed to calibrate the complex structure were determined. The non-bonded interactions (Electrostatic and van der Waals contributions), H-bond formation, and salt bridges between the ligand HLA-Cw4 and the receptor KIR2DL1 were estimated using the obtained MD trajectories. The buried surface area due to binding was calculated to get insight into the causes of specificity of receptor to ligand and explains mutations experiment. The study concluded that β2-microglobulin, one part of the complex, is not directly interacting with the peptide at the groove; therefore, it could be neglected from simulation. Our results showed that β2-microglobulin does not have any significant effect on the dynamics of the 3D-structure of the complex. This project will help in understanding to optimize candidate drug design, a small peptide that disrupts the interaction, for the optimal biological effect.

scholarly journals Molecular Insights into Human Transmembrane Protease Serine-2 (TMPS2) Inhibitors against SARS-CoV2: Homology Modelling, Molecular Dynamics, and Docking Studies

Molecules ◽  
2020 ◽  
Vol 25 (21) ◽  
pp. 5007
Author(s):  
Safaa M. Kishk ◽  
Rania M. Kishk ◽  
Asmaa S. A. Yassen ◽  
Mohamed S. Nafie ◽  
Nader A. Nemr ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Serine Protease ◽  
Tertiary Structure ◽  
Virus Disease ◽  
Homology Modelling ◽  
3D Structure ◽  
Md Simulations ◽  
Homology Model ◽  
Docking Studies ◽  
In Silico Studies

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV2), which caused novel corona virus disease-2019 (COVID-19) pandemic, necessitated a global demand for studies related to genes and enzymes of SARS-CoV2. SARS-CoV2 infection depends on the host cell Angiotensin-Converting Enzyme-2 (ACE2) and Transmembrane Serine Protease-2 (TMPRSS2), where the virus uses ACE2 for entry and TMPRSS2 for S protein priming. The TMPRSS2 gene encodes a Transmembrane Protease Serine-2 protein (TMPS2) that belongs to the serine protease family. There is no crystal structure available for TMPS2, therefore, a homology model was required to establish a putative 3D structure for the enzyme. A homology model was constructed using SWISS-MODEL and evaluations were performed through Ramachandran plots, Verify 3D and Protein Statistical Analysis (ProSA). Molecular dynamics simulations were employed to investigate the stability of the constructed model. Docking of TMPS2 inhibitors, camostat, nafamostat, gabexate, and sivelestat, using Molecular Operating Environment (MOE) software, into the constructed model was performed and the protein-ligand complexes were subjected to MD simulations and computational binding affinity calculations. These in silico studies determined the tertiary structure of TMPS2 amino acid sequence and predicted how ligands bind to the model, which is important for drug development for the prevention and treatment of COVID-19.

scholarly journals Fine-tuning of the AMBER RNA Force Field with a New Term Adjusting Interactions of Terminal Nucleotides

2020 ◽  
Author(s):  
Vojtěch Mlýnský ◽  
Petra Kührová ◽  
Tomáš Kühr ◽  
Michal Otyepka ◽  
Giovanni Bussi ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
State Of The Art ◽  
Dynamic Properties ◽  
Md Simulations ◽  
Experimental Methods ◽  
Fine Tuning ◽  
Enhanced Sampling ◽  
Structural Dynamic ◽  
Nmr Data

ABSTRACTDetermination of RNA structural-dynamic properties is challenging for experimental methods. Thus atomistic molecular dynamics (MD) simulations represent a helpful technique complementary to experiments. However, contemporary MD methods still suffer from limitations of force fields (ffs), including imbalances in the non-bonded ff terms. We have recently demonstrated that some improvement of state-of-the-art AMBER RNA ff can be achieved by adding a new term for H-bonding called gHBfix, which increases tuning flexibility and reduces the risk of side-effects. Still, the first gHBfix version did not fully correct simulations of short RNA tetranucleotides (TNs). TNs are key benchmark systems due to availability of unique NMR data, although giving too much weight on improving TN simulations can easily lead to over-fitting to A-form RNA. Here we combine the gHBfix version with another term called tHBfix, which separately treats H-bond interactions formed by terminal nucleotides. This allows to refine simulations of RNA TNs without affecting simulations of other RNAs. The approach is in line with adopted strategy of current RNA ffs, where the terminal nucleotides possess different parameters for the terminal atoms than the internal nucleotides. The combination of gHBfix with tHBfix significantly improves the behavior of RNA TNs during well-converged enhanced-sampling simulations. TNs mostly populate canonical A-form like states while spurious intercalated structures are largely suppressed. Still, simulations of r(AAAA) and r(UUUU) TNs show some residual discrepancies with the primary NMR data which suggests that future tuning of some other ff terms might be useful.

scholarly journals Investigating the conformational dynamics of Zika virus NS4B protein

2021 ◽  
Author(s):  
Taniya Bhardwaj ◽  
Prateek Kumar ◽  
Rajanish Giri
Keyword(s):  
Molecular Dynamics ◽  
Zika Virus ◽  
Conformational Dynamics ◽  
3D Structure ◽  
Md Simulations ◽  
Cd Spectroscopy ◽  
Host Immune System ◽  
Genome Replication ◽  
Intrinsically Disordered ◽  
Ns4b Protein

Zika virus (ZIKV) NS4B protein is a membranotropic protein having multifunctional roles such as evasion of host-immune system, and induction of host membrane rearrangement for viral genome replication and processing of polyprotein. Despite its versatile functioning, its topology and dynamics are not entirely understood. Presently, there is no NMR or X-ray crystallography-based structure available for any flaviviral NS4B protein. Therefore, in this study, we have investigated the structural dynamics of Zika Virus NS4B protein through 3D structure models using molecular dynamics (MD) simulations approach and experiments. Subsequently, we employed a reductionist approach to understand the dynamics of ZIKV NS4B protein. For this, we studied its N-terminal (residues 1-38), C-terminal (residues 194-251), and cytosolic (residues 131-169) regions in isolation. Further, we have performed experiments to prove the maximum dynamics in its cytosolic region. Using a combination of computational tools and circular dichroism (CD) spectroscopy, we validate the cytosolic region as an intrinsically disordered protein region (IDR). The microsecond-long all atoms molecular dynamics and replica-exchange MD simulations complement the experimental observations. We have also analysed its behaviour under the influence of differently charged liposomes and macromolecular crowding agents which may have significance on its overall dynamics. Lastly, we have proposed a ZIKV NS4B protein model illustrating its putative topology consisting of various membrane-spanning and non-membranous regions.

FLEXIBILITY AND MOBILITY IN MESOPHILIC AND THERMOPHILIC HOMOLOGOUS PROTEINS FROM MOLECULAR DYNAMICS AND FOLDUNFOLD METHOD

2010 ◽  
Vol 08 (03) ◽  
pp. 377-394 ◽  
Author(s):  
TATYANA B. MAMONOVA ◽  
ANNA V. GLYAKINA ◽  
MARIA G. KURNIKOVA ◽  
OXANA V. GALZITSKAYA
Keyword(s):  
Molecular Dynamics ◽  
3D Structure ◽  
Three Dimensional ◽  
High Mobility ◽  
Md Simulations ◽  
Homologous Proteins ◽  
Knowledge Based ◽  
Theoretical Approaches ◽  
High Flexibility ◽  
A Site

To function properly protein molecules require both flexibility and rigidity, therefore fast and accurate prediction of protein rigidity/flexibility is one of the important problems in protein science. In this work we used two theoretical approaches to determine flexible regions in four homologous pairs of proteins from thermophilic and mesophilic organisms. Protein pairs chosen in this study were selected to represent four typical folding classes. Our first approach, FoldUnfold, uses amino acid sequence and statistical information on the density of contacts of amino acids in tertiary structures of known globular proteins. The main advantages of such knowledge-based methodology are its computational speed and ability to make predictions in the absence of three-dimensional (3D) structure of a protein. The second approach uses a graph theory-based rigid cluster decomposition termed FIRST, applied together with Molecular Dynamics (MD) simulations of proteins with known structure. While MD simulations are time-consuming, they are the most direct way of studying physical properties of proteins, including their rigidity/flexibility. Flexible regions predicted by both methods in this work were in good agreement with each other. We also showed that high mobility of a site is not necessarily indicative of its high flexibility and vice versa. In our simulations thermophile proteins were less flexible than their mesophilic homologues. Longer flexible loops were found in mesophilic proteins of all classes.

scholarly journals Human hydroxymethylbilane synthase: Molecular dynamics of the pyrrole chain elongation identifies step-specific residues that cause AIP

2018 ◽  
Vol 115 (17) ◽  
pp. E4071-E4080 ◽  
Author(s):  
Navneet Bung ◽  
Arijit Roy ◽  
Brenden Chen ◽  
Dibyajyoti Das ◽  
Meenakshi Pradhan ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Active Site ◽  
3D Structure ◽  
Acute Intermittent Porphyria ◽  
Md Simulations ◽  
Chain Elongation ◽  
Active Site Residues ◽  
Autosomal Dominant Disorder ◽  
Random Acceleration

Hydroxymethylbilane synthase (HMBS), the third enzyme in the heme biosynthetic pathway, catalyzes the head-to-tail condensation of four molecules of porphobilinogen (PBG) to form the linear tetrapyrrole 1-hydroxymethylbilane (HMB). Mutations in human HMBS (hHMBS) cause acute intermittent porphyria (AIP), an autosomal-dominant disorder characterized by life-threatening neurovisceral attacks. Although the 3D structure of hHMBS has been reported, the mechanism of the stepwise polymerization of four PBG molecules to form HMB remains unknown. Moreover, the specific roles of each of the critical active-site residues in the stepwise enzymatic mechanism and the dynamic behavior of hHMBS during catalysis have not been investigated. Here, we report atomistic studies of HMB stepwise synthesis by using molecular dynamics (MD) simulations, mutagenesis, and in vitro expression analyses. These studies revealed that the hHMBS active-site loop movement and cofactor turn created space for the elongating pyrrole chain. Twenty-seven residues around the active site and water molecules interacted to stabilize the large, negatively charged, elongating polypyrrole. Mutagenesis of these active-site residues altered the binding site, hindered cofactor binding, decreased catalysis, impaired ligand exit, and/or destabilized the enzyme. Based on intermediate stages of chain elongation, R26 and R167 were the strongest candidates for proton transfer to deaminate the incoming PBG molecules. Unbiased random acceleration MD simulations identified R167 as a gatekeeper and facilitator of HMB egress through the space between the enzyme’s domains and the active-site loop. These studies identified the specific active-site residues involved in each step of pyrrole elongation, thereby providing the molecular bases of the active-site mutations causing AIP.

Efficient Treatment of Fixed and Induced Dipolar Interactions

2000 ◽  
Vol 653 ◽  
Author(s):  
Celeste Sagui ◽  
Thoma Darden
Keyword(s):  
Molecular Dynamics ◽  
Md Simulations ◽  
Lagrangian Formalism ◽  
Dipolar Interactions ◽  
Time Step ◽  
Efficient Treatment ◽  
Particle Mesh Ewald ◽  
Fixed Charges ◽  
Self Consistency ◽  
Induced Dipoles

AbstractFixed and induced point dipoles have been implemented in the Ewald and Particle-Mesh Ewald (PME) formalisms. During molecular dynamics (MD) the induced dipoles can be propagated along with the atomic positions either by interation to self-consistency at each time step, or by a Car-Parrinello (CP) technique using an extended Lagrangian formalism. The use of PME for electrostatics of fixed charges and induced dipoles together with a CP treatment of dipole propagation in MD simulations leads to a cost overhead of only 33% above that of MD simulations using standard PME with fixed charges, allowing the study of polarizability in largemacromolecular systems.

Split the Charge Difference in Two! a Rule of Thumb for Adding Proper Amounts of Ions in MD Simulations

2020 ◽  
Author(s):  
Matías R. Machado ◽  
Sergio Pantano
Keyword(s):  
Molecular Dynamics ◽  
Ionic Strength ◽  
Molecular Simulations ◽  
Md Simulations ◽  
Good Practices ◽  
Rule Of Thumb ◽  
Ionic Concentrations

<p> Despite the relevance of properly setting ionic concentrations in Molecular Dynamics (MD) simulations, methods or practical rules to set ionic strength are scarce and rarely documented. Based on a recently proposed thermodynamics method we provide an accurate rule of thumb to define the electrolytic content in simulation boxes. Extending the use of good practices in setting up MD systems is promptly needed to ensure reproducibility and consistency in molecular simulations.</p>

Assessing the Antimalarial Potentials of Phytochemicals: Virtual Screening, Molecular Dynamics and In-Vitro Investigations

2019 ◽  
Vol 16 (3) ◽  
pp. 291-300
Author(s):  
Saumya K. Patel ◽  
Mohd Athar ◽  
Prakash C. Jha ◽  
Vijay M. Khedkar ◽  
Yogesh Jasrai ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Integrity ◽  
Md Simulations ◽  
Tylophora Indica ◽  
Multidrug Resistance Protein 1 ◽  
Receptor Complexes ◽  
Resistance Protein ◽  
Dynamics Simulations ◽  
Stable Trajectory

Background: Combined in-silico and in-vitro approaches were adopted to investigate the antiplasmodial activity of Catharanthus roseus and Tylophora indica plant extracts as well as their isolated components (vinblastine, vincristine and tylophorine). </P><P> Methods: We employed molecular docking to prioritize phytochemicals from a library of 26 compounds against Plasmodium falciparum multidrug-resistance protein 1 (PfMDR1). Furthermore, Molecular Dynamics (MD) simulations were performed for a duration of 10 ns to estimate the dynamical structural integrity of ligand-receptor complexes. </P><P> Results: The retrieved bioactive compounds viz. tylophorine, vinblastin and vincristine were found to exhibit significant interacting behaviour; as validated by in-vitro studies on chloroquine sensitive (3D7) as well as chloroquine resistant (RKL9) strain. Moreover, they also displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations. </P><P> Conclusion: We anticipate that the retrieved phytochemicals can serve as the potential hits and presented findings would be helpful for the designing of malarial therapeutics.

Sign in / Sign up

Export Citation Format

Share Document