Investigating coordination flexibility of glycerophosphodiesterase (GpdQ) through interactions with mono-, di-, and triphosphoester (NPP, BNPP, GPE, and paraoxon) substrates

2019 ◽  
Vol 21 (10) ◽  
pp. 5499-5509 ◽  
Author(s):  
Gaurav Sharma ◽  
Qiaoyu Hu ◽  
Vindi M. Jayasinghe-Arachchige ◽  
Thomas J. Paul ◽  
Gerhard Schenk ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Md Simulations ◽  
Catalytically Active

Interactions of the catalytically active binuclear form of glycerophosphodiesterase (GpdQ) with chemically diverse substrates, i.e. phosphomono-, phosphodi-, and phosphotriester have been investigated using molecular dynamics (MD) simulations.

scholarly journals Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules

2021 ◽  
Vol 8 ◽  
Author(s):  
Łukasz Nierzwicki ◽  
Giulia Palermo
Keyword(s):  
Molecular Dynamics ◽  
Md Simulations ◽  
Atomic Level ◽  
Conformational States ◽  
Variable Regions ◽  
Biologically Relevant ◽  
Active Structure ◽  
Conformational Variability ◽  
Challenges And Opportunities ◽  
Catalytically Active

Single-particle cryogenic electron microscopy (cryo-EM) has revolutionized the field of the structural biology, providing an access to the atomic resolution structures of large biomolecular complexes in their near-native environment. Today’s cryo-EM maps can frequently reach the atomic-level resolution, while often containing a range of resolutions, with conformationally variable regions obtained at 6 Å or worse. Low resolution density maps obtained for protein flexible domains, as well as the ensemble of coexisting conformational states arising from cryo-EM, poses new challenges and opportunities for Molecular Dynamics (MD) simulations. With the ability to describe the biomolecular dynamics at the atomic level, MD can extend the capabilities of cryo-EM, capturing the conformational variability and predicting biologically relevant short-lived conformational states. Here, we report about the state-of-the-art MD procedures that are currently used to refine, reconstruct and interpret cryo-EM maps. We show the capability of MD to predict short-lived conformational states, finding remarkable confirmation by cryo-EM structures subsequently solved. This has been the case of the CRISPR-Cas9 genome editing machinery, whose catalytically active structure has been predicted through both long-time scale MD and enhanced sampling techniques 2 years earlier than cryo-EM. In summary, this contribution remarks the ability of MD to complement cryo-EM, describing conformational landscapes and relating structural transitions to function, ultimately discerning relevant short-lived conformational states and providing mechanistic knowledge of biological function.

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.

scholarly journals On the Use of the Discrete Constant pH Molecular Dynamics to Describe the Conformational Space of Peptides

Polymers ◽  
2020 ◽  
Vol 13 (1) ◽  
pp. 99
Author(s):  
Cristian Privat ◽  
Sergio Madurga ◽  
Francesc Mas ◽  
Jaime Rubio-Martínez
Keyword(s):  
Molecular Dynamics ◽  
Amino Acids ◽  
Md Simulations ◽  
Point Of View ◽  
Conformational Space ◽  
Conformational Sampling ◽  
Protonation State ◽  
Constant Ph ◽  
Biochemical Systems ◽  
Constant Ph Molecular Dynamics

Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.

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 Understanding Covalent Grafting of Nanotubes onto Polymer Nanocomposites: Molecular Dynamics Simulation Study

Sensors ◽  
2021 ◽  
Vol 21 (8) ◽  
pp. 2621
Author(s):  
Seunghwa Yang
Keyword(s):  
Molecular Dynamics ◽  
Load Transfer ◽  
Cell Model ◽  
Md Simulations ◽  
Dynamics Simulation ◽  
Covalent Grafting ◽  
Modelling Strategies ◽  
Multi Walled Carbon Nanotubes ◽  
Walled Carbon Nanotubes ◽  
External Loads

Here, we systematically interrogate the effects of grafting single-walled (SWNT) and multi-walled carbon nanotubes (MWNT) to polymer matrices by using molecular dynamics (MD) simulations. We specifically investigate key material properties that include interfacial load transfer, alteration of nanotube properties, and dispersion of nanotubes in the polymer matrix. Simulations are conducted on a periodic unit cell model of the nanocomposite with a straight carbon nanotube and an amorphous polyethylene terephthalate (PET) matrix. For each type of nanotube, either 0%, 1.55%, or 3.1% of the carbon atoms in the outermost nanotubes are covalently grafted onto the carbon atoms of the PET matrix. Stress-strain curves and the elastic moduli of nanotubes and nanocomposites are determined based on the density of covalent grafting. Covalent grafting promotes two rivalling effects with respect to altering nanotube properties, and improvements in interfacial load transfer in the nanocomposites are clearly observed. The enhanced interface enables external loads applied to the nanocomposites to be efficiently transferred to the grafted nanotubes. Covalent functionalization of the nanotube surface with PET molecules can alter the solubility of nanotubes and improve dispersibility. Finally, we discuss the current limitations and challenges in using molecular modelling strategies to accurately predict properties on the nanotube and polymers systems studied here.

scholarly journals Ligand-induced structural changes analysis of ribose-binding protein as studied by molecular dynamics simulations

2021 ◽  
pp. 1-12
Author(s):  
Haiyan Li ◽  
Zanxia Cao ◽  
Guodong Hu ◽  
Liling Zhao ◽  
Chunling Wang ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Changes ◽  
Binding Protein ◽  
Network Models ◽  
Md Simulations ◽  
Protein Domain ◽  
Opening Angle ◽  
Conformation Changes ◽  
Wide Range ◽  
Ribose Binding Protein

BACKGROUND: The ribose-binding protein (RBP) from Escherichia coli is one of the representative structures of periplasmic binding proteins. Binding of ribose at the cleft between two domains causes a conformational change corresponding to a closure of two domains around the ligand. The RBP has been crystallized in the open and closed conformations. OBJECTIVE: With the complex trajectory as a control, our goal was to study the conformation changes induced by the detachment of the ligand, and the results have been revealed from two computational tools, MD simulations and elastic network models. METHODS: Molecular dynamics (MD) simulations were performed to study the conformation changes of RBP starting from the open-apo, closed-holo and closed-apo conformations. RESULTS: The evolution of the domain opening angle θ clearly indicates large structural changes. The simulations indicate that the closed states in the absence of ribose are inclined to transition to the open states and that ribose-free RBP exists in a wide range of conformations. The first three dominant principal motions derived from the closed-apo trajectories, consisting of rotating, bending and twisting motions, account for the major rearrangement of the domains from the closed to the open conformation. CONCLUSIONS: The motions showed a strong one-to-one correspondence with the slowest modes from our previous study of RBP with the anisotropic network model (ANM). The results obtained for RBP contribute to the generalization of robustness for protein domain motion studies using either the ANM or PCA for trajectories obtained from MD.

Sign in / Sign up

Export Citation Format

Share Document