scholarly journals Sensitivity of Molecular Dynamics Simulations to Equilibration Scheme: A Case Study of Bromodomain Protein BRD4-Ligand Complex System

2021 ◽  
Author(s):  
Yunhui Ge ◽  
Gaetano Calabro ◽  
Christophe Bayly ◽  
David Mobley
Keyword(s):  
Complex System ◽  
Molecular Dynamics ◽  
Binding Site ◽  
Md Simulations ◽  
Careful Consideration ◽  
Ligand Complex ◽  
Thermodynamic Conditions ◽  
Temperature And Pressure ◽  
Dynamics Simulations

In molecular dynamics (MD) simulations, an equilibration phase is used to bring the system to the desired conditions (e.g. temperature and pressure) before collecting data in more extensive production phases. Ideally, the equilibration phase would bring the system to the correct equilibrium ensemble appropriate for the target thermodynamic conditions. Many studies give relatively little attention to details of the equilibration protocol as long as the system eventually reaches stabilization at the correct temperature, pressure, volume, etc. However, in a previous study we found a surprising instability of two ligands in a binding site. That led us to study the origin of this instability more, and in the present work we have traced it to details of the equilibration protocol. We found for the studied system, different equilibration schemes caused dramatic differences in sampled configurations of the system in production runs, highlighting the importance of careful consideration of equilibration schemes in MD simulations.

scholarly journals Accelerated molecular dynamics simulations of ligand binding to a muscarinic G-protein-coupled receptor

2015 ◽  
Vol 48 (4) ◽  
pp. 479-487 ◽  
Author(s):  
Kalli Kappel ◽  
Yinglong Miao ◽  
J. Andrew McCammon
Keyword(s):  
Molecular Dynamics ◽  
Ligand Binding ◽  
Binding Site ◽  
G Protein ◽  
Partial Agonist ◽  
Md Simulations ◽  
G Protein Coupled Receptor ◽  
Accelerated Molecular Dynamics ◽  
Dynamics Simulations ◽  
G Protein Coupled

AbstractElucidating the detailed process of ligand binding to a receptor is pharmaceutically important for identifying druggable binding sites. With the ability to provide atomistic detail, computational methods are well poised to study these processes. Here, accelerated molecular dynamics (aMD) is proposed to simulate processes of ligand binding to a G-protein-coupled receptor (GPCR), in this case the M3 muscarinic receptor, which is a target for treating many human diseases, including cancer, diabetes and obesity. Long-timescale aMD simulations were performed to observe the binding of three chemically diverse ligand molecules: antagonist tiotropium (TTP), partial agonist arecoline (ARc) and full agonist acetylcholine (ACh). In comparison with earlier microsecond-timescale conventional MD simulations, aMD greatly accelerated the binding of ACh to the receptor orthosteric ligand-binding site and the binding of TTP to an extracellular vestibule. Further aMD simulations also captured binding of ARc to the receptor orthosteric site. Additionally, all three ligands were observed to bind in the extracellular vestibule during their binding pathways, suggesting that it is a metastable binding site. This study demonstrates the applicability of aMD to protein–ligand binding, especially the drug recognition of GPCRs.

scholarly journals Pharmacophore Models Derived from Molecular Dynamics Simulations of Protein-Ligand Complexes: A Case Study

2016 ◽  
Vol 11 (10) ◽  
pp. 1934578X1601101 ◽  
Author(s):  
Marcus Wieder ◽  
Ugo Perricone ◽  
Thomas Seidel ◽  
Thierry Langer
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulation ◽  
Pharmacophore Model ◽  
Dynamics Simulation ◽  
Ligand Complex ◽  
Naturally Occurring ◽  
Pharmacophore Hypothesis ◽  
Dynamics Simulations ◽  
Pharmacophore Models

A single, merged pharmacophore hypothesis is derived combining 2000 pharmacophore models obtained during a 20 ns molecular dynamics simulation of a protein-ligand complex with one pharmacophore model derived from the initial PDB structure. This merged pharmacophore model contains all features that are present during the simulation and statistical information about the dynamics of the pharmacophore features. Based on the dynamics of the pharmacophore features we derive two distinctive feature patterns resulting in two different pharmacophore models for the analyzed system – the first model consists of features that are obtained from the PDB structure and the second uses two features that can only be derived from the molecular dynamics simulation. Both models can distinguish between active and decoy molecules in virtual screening. Our approach represents an objective way to add/remove features in pharmacophore models and can be of interest for the investigation of any naturally occurring system that relies on ligand-receptor interactions for its biological activity.

scholarly journals Deciphering collaborative sidechain motions in proteins during molecular dynamics simulations

2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Bruck Taddese ◽  
Antoine Garnier ◽  
Hervé Abdi ◽  
Daniel Henrion ◽  
Marie Chabbert
Keyword(s):  
Molecular Dynamics ◽  
Large Scale ◽  
Conformational Transition ◽  
Transmembrane Helix ◽  
Principal Component ◽  
Dynamic Structure ◽  
Md Simulations ◽  
Outward Motion ◽  
Dynamics Simulations

Abstract The dynamic structure of proteins is essential for their functions and may include large conformational transitions which can be studied by molecular dynamics (MD) simulations. However, details of these transitions are difficult to automatically track. To facilitate their analysis, we developed two scores of correlation between sidechain dihedral angles. The CIRCULAR and OMES scores are computed from, respectively, dihedral angle values and rotamer distributions. As a case study, we applied our methods to an activation-like transition of the chemokine receptor CXCR4, observed during accelerated MD simulations. The principal component analysis of the correlation matrices was consistent with the networking structure of the top ranking pairs. Both scores identify a set of residues whose “collaborative” sidechain rotamerization immediately preceded or accompanied the conformational transition of CXCR4. Detailed analysis of the sequential order of these rotamerizations suggests that an allosteric mechanism, involving the outward motion of an asparagine residue in transmembrane helix 3, might be a prerequisite to the large scale conformational transition of CXCR4. This case study provides the proof-of-concept that the correlation methods developed here are valuable exploratory techniques to help decipher complex reactional pathways.

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.

Design of Inhibitors for Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH) Enzyme of Leishmania mexicana

2020 ◽  
Vol 16 (6) ◽  
pp. 784-795
Author(s):  
Krisnna M.A. Alves ◽  
Fábio José Bonfim Cardoso ◽  
Kathia M. Honorio ◽  
Fábio A. de Molfetta
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulations ◽  
Binding Site ◽  
Point Of View ◽  
New Drugs ◽  
Leishmania Mexicana ◽  
Receptor Complexes ◽  
Starting Point ◽  
Dynamics Simulations ◽  
Selection Of

Background:: Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease. Objective:: The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana). Methods: A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations. Results:: Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity. Conclusion:: he use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.

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.

scholarly journals Molecular dynamics simulations and CD spectroscopy reveal hydration-induced unfolding of the intrinsically disordered LEA proteins COR15A and COR15B from Arabidopsis thaliana

2016 ◽  
Vol 18 (37) ◽  
pp. 25806-25816 ◽  
Author(s):  
Carlos Navarro-Retamal ◽  
Anne Bremer ◽  
Jans Alzate-Morales ◽  
Julio Caballero ◽  
Dirk K. Hincha ◽  
...  
Keyword(s):  
Experimental Data ◽  
Molecular Dynamics ◽  
Arabidopsis Thaliana ◽  
Molecular Dynamics Simulations ◽  
Md Simulations ◽  
Cd Spectroscopy ◽  
Lea Proteins ◽  
Intrinsically Disordered ◽  
Dynamics Simulations ◽  
Crowded Environment

Unfolding of intrinsically unstructured full-length LEA proteins in a differentially crowded environment can be modeled by 30 ns MD simulations in accordance with experimental data.

Sign in / Sign up

Export Citation Format

Share Document