scholarly journals Comparing state-of-the-art approaches to back-calculate SAXS spectra from atomistic molecular dynamics simulations

2021 ◽  
Vol 94 (9) ◽  
Author(s):  
Mattia Bernetti ◽  
Giovanni Bussi
Keyword(s):  
Molecular Dynamics ◽  
Solvent Effects ◽  
Conformational Dynamics ◽  
Md Simulations ◽  
Implicit Solvent ◽  
Remarkable Effect ◽  
Rna Molecules ◽  
Reaction Coordinates ◽  
Effective Instrument ◽  
Dynamics Simulations

Abstract Small-angle X-ray scattering (SAXS) experiments are arising as an effective instrument in the structural characterization of biomolecules in solution. However, they suffer from limited resolution, and complementing them with molecular dynamics (MD) simulations can be a successful strategy to obtain information at a finer scale. To this end, tools that allow computing SAXS spectra from MD-sampled structures have been designed over the years, mainly differing in how the solvent contribution is accounted for. In this context, RNA molecules represent a particularly challenging case, as they can have a remarkable effect on the surrounding solvent. Herein, we provide a comparison of SAXS spectra computed through different available software packages for a prototypical RNA system. RNA conformational dynamics is intentionally neglected so as to focus on solvent effects. The results highlight that solvent effects are important also at relatively low scattering vector, suggesting that approaches explicitly modeling solvent contribution are advisable when comparing with experimental data, while more efficient implicit-solvent methods can be a better choice as reaction coordinates to improve MD sampling on-the-fly. Graphic abstract

scholarly journals Molecular dynamics simulations disclose early stages of the photo-activation of cryptochrome 4

2018 ◽  
Author(s):  
D. R. Kattnig ◽  
C. Nielsen ◽  
I. A. Solov’yov
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Large Scale ◽  
Conformational Dynamics ◽  
Md Simulations ◽  
Quantum Processes ◽  
Helical Segment ◽  
High Betweenness ◽  
Dynamics Simulations ◽  
Fad Binding

AbstractBirds appear to be equipped with a light-dependent, radical-pair-based magnetic compass that relies on truly quantum processes. While the identity of the sensory protein has remained speculative, cryptochrome 4 has recently been identified as the most auspicious candidate. Here, we report on allatom molecular dynamics (MD) simulations addressing the structural reorganisations that accompany the photoreduction of the flavin cofactor in a model of the European robin cryptochrome 4 (ErCry4). Extensive MD simulations reveal that the photo-activation of ErCry4 induces large-scale conformational changes on short (hundreds of nanoseconds) timescales. Specifically, the photo-reduction is accompanied with the release of the C-terminal tail, structural rearrangements in the vicinity of the FAD-binding site, and the noteworthy formation of an α-helical segment at the N-terminal part. Some of these rearrangements appear to expose potential phosphorylation sites. We describe the conformational dynamics of the protein using a graph-based approach that is informed by the adjacency of residues and the correlation of their local motions. This approach reveals densely coupled reorganisation communities, which facilitate an efficient signal transduction due to a high density of hubs. These communities are interconnected by a small number of highly important residues characterized by high betweenness centrality. The network approach clearly identifies the sites restructuring upon photoactivation, which appear as protrusions or delicate bridges in the reorganisation network. We also find that, unlike in the homologous cryptochrome from D. melanogaster, the release of the C-terminal domain does not appear to be correlated with the transposition of a histidine residue close to the FAD cofactor.

scholarly journals Insights from Molecular Dynamics Simulations for GPCR Drug Discovery

2019 ◽  
Author(s):  
Ye Zou ◽  
John Ewalt ◽  
Ho-Leung Ng
Keyword(s):  
Molecular Dynamics ◽  
Drug Discovery ◽  
G Protein ◽  
Drug Targets ◽  
Conformational Dynamics ◽  
Md Simulations ◽  
Protein Activation ◽  
Downstream Signaling ◽  
Dynamics Simulations ◽  
G Protein Coupled

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. There is evidence that some ligands that bind to the GPCRs activate different downstream signaling pathways. G protein activation or -arrestin biased signaling involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. Molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions between receptors and ligands. Here, we highlight recent simulation studies and methods used to study biased G protein-coupled receptor signaling and their conformational dynamics. We also highlight applications of MD simulations to drug discovery.

scholarly journals The Assembly Switch Mechanism of FtsZ Filament Revealed by All-Atom Molecular Dynamics Simulations and Coarse-Grained Models

2021 ◽  
Vol 12 ◽  
Author(s):  
Dashuai Lv ◽  
Jingyuan Li ◽  
Sheng Ye
Keyword(s):  
Molecular Dynamics ◽  
Bound States ◽  
Conformational Dynamics ◽  
Md Simulations ◽  
Coarse Grained ◽  
Multi Scale ◽  
Switch Mechanism ◽  
Dynamics Simulations ◽  
Cooperative Assembly ◽  
Computational Strategy

Bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, and assembles into dynamic filaments that are essential for cell division. Here, we used a multi-scale computational strategy that combined all-atom molecular dynamics (MD) simulations and coarse-grained models to reveal the conformational dynamics of assembled FtsZ. We found that the top end of a filament is highly dynamic and can undergo T-to-R transitions in both GTP- and GDP-bound states. We observed several subcategories of nucleation related dimer species, which leading to a feasible nucleation pathway. In addition, we observed that FtsZ filament exhibits noticeable amounts of twisting, indicating a substantial helicity of the FtsZ filament. These results agree with the previously models and experimental data. Anisotropy network model (ANM) analysis revealed a polymerization enhanced assembly cooperativity, and indicated that the cooperative motions in FtsZ are encoded in the structure. Taken together, our study provides a molecular-level understanding of the diversity of the structural states of FtsZ and the relationships among polymerization, hydrolysis, and cooperative assembly, which should shed new light on the molecular basis of FtsZ’s cooperativity.

Molecular Basis of Glucose Transport Mechanism in Plants

2019 ◽  
Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla
Keyword(s):  
Molecular Dynamics ◽  
Free Energy ◽  
Molecular Dynamics Simulations ◽  
Conformational Changes ◽  
Transport Mechanism ◽  
Conformational Dynamics ◽  
Site Directed Mutagenesis ◽  
Free Energy Barrier ◽  
Transport Cycle ◽  
Dynamics Simulations

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</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.

Schistosomal Sulfotransferase Interaction with Oxamniquine Involves Hybrid Mechanism of Induced-fit and Conformational Selection

2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu
Keyword(s):  
Molecular Dynamics ◽  
Amino Acid ◽  
Binding Energy ◽  
Binding Site ◽  
Conformational Dynamics ◽  
Side Chain ◽  
Amino Acid Residues ◽  
Conformational Selection ◽  
Hybrid Mechanism ◽  
Dynamics Simulations

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.

scholarly journals Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau
Keyword(s):  
Carbon Dioxide ◽  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Hydrodynamic Radius ◽  
Md Simulations ◽  
Molecular Models ◽  
Experimental Values ◽  
Carbon Dioxide Co2 ◽  
Dynamics Simulations ◽  
Apparent Disagreement

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.

Mechanism of stacking fault annihilation in 3C-SiC epitaxially grown on Si(001) by molecular dynamics simulations

CrystEngComm ◽  
2021 ◽  
Author(s):  
Andrey Sarikov ◽  
Anna Marzegalli ◽  
Luca Barbisan ◽  
Massimo Zimbone ◽  
Corrado Bongiorno ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Stacking Faults ◽  
Md Simulations ◽  
Stacking Fault ◽  
Dynamics Simulations

In this work, annihilation mechanism of stacking faults (SFs) in epitaxial 3C-SiC layers grown on Si(001) substrates is studied by molecular dynamics (MD) simulations. The evolution of SFs located in...

scholarly journals Molecular Dynamics Simulations of Crystal Nucleation near Interfaces in Incompatible Polymer Blends

Polymers ◽  
2021 ◽  
Vol 13 (3) ◽  
pp. 347
Author(s):  
Wenlin Zhang ◽  
Lingyi Zou
Keyword(s):  
Molecular Dynamics ◽  
Polymer Blends ◽  
Md Simulations ◽  
Nucleation And Growth ◽  
Crystal Nucleation ◽  
Crystalline Order ◽  
Mobile Species ◽  
Deep Supercooling ◽  
Crystallizable Polymer ◽  
Dynamics Simulations

We apply molecular dynamics (MD) simulations to investigate crystal nucleation in incompatible polymer blends under deep supercooling conditions. Simulations of isothermal nucleation are performed for phase-separated blends with different degrees of incompatibility. In weakly segregated blends, slow and incompatible chains in crystallizable polymer domains can significantly hinder the crystal nucleation and growth. When a crystallizable polymer is blended with a more mobile species in interfacial regions, enhanced molecular mobility leads to the fast growth of crystalline order. However, the incubation time remains the same as that in pure samples. By inducing anisotropic alignment near the interfaces of strongly segregated blends, phase separation also promotes crystalline order to grow near interfaces between different polymer domains.

Sign in / Sign up

Export Citation Format

Share Document