Dataset Information

Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization.

ABSTRACT: The molecular simulation of chemical reaction equilibrium (CRE) is a challenging and important problem of broad applicability in chemistry and chemical engineering. The primary molecular-based approach for solving this problem has been the reaction ensemble Monte Carlo (REMC) algorithm [Turner et al. Molec. Simulation2008, 34, (2), 119-146], based on classical force-field methodology. In spite of the vast improvements in computer hardware and software since its original development almost 25 years ago, its more widespread application is impeded by its computational inefficiency. A fundamental problem is that its MC basis inhibits the implementation of significant parallelization, and its successful implementation often requires system-specific tailoring and the incorporation of special MC approaches such as replica exchange, expanded ensemble, umbrella sampling, configurational bias, and continuous fractional component methodologies. We describe herein a novel CRE algorithm (reaction ensemble molecular dynamics, ReMD) that exploits modern computer hardware and software capabilities, and which can be straightforwardly implemented for systems of arbitrary size and complexity by exploiting the parallel computing methodology incorporated within many MD software packages (herein, we use GROMACS for illustrative purposes). The ReMD algorithm utilizes these features in the context of a macroscopically inspired and generally applicable free energy minimization approach based on the iterative approximation of the system Gibbs free energy function by a mathematically simple convex ideal solution model using the composition at each iteration as a reference state. Finally, we additionally describe a simple and computationally efficient a posteriori method to estimate the equilibrium concentrations of species present in very small amounts relative to others in the primary calculation. To demonstrate the algorithm, we show its application to two classic example systems considered previously in the literature: the N₂-O₂-NO system and the ammonia synthesis system.

SUBMITTER: Smith WR

PROVIDER: S-EPMC6161046 | biostudies-literature | 2018 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization.

Smith William R WR Qi Weikai W

ACS central science 20180823 9

The molecular simulation of chemical reaction equilibrium (CRE) is a challenging and important problem of broad applicability in chemistry and chemical engineering. The primary molecular-based approach for solving this problem has been the reaction ensemble Monte Carlo (REMC) algorithm [Turner et al. Molec. Simulation2008, 34, (2), 119-146], based on classical force-field methodology. In spite of the vast improvements in computer hardware and software since its original development almost 25 yea ...[more]

PMID: 30276252

Similar Datasets

Project description:RNA hairpin loops are one of the most abundant secondary structure elements and participate in RNA folding and protein-RNA recognition. To characterize the free energy surface of RNA hairpin folding at an atomic level, we calculated the potential of mean force (PMF) as a function of the end-to-end distance, by using umbrella sampling simulations in explicit solvent. Two RNA hairpins containing tetraloop cUUCGg and cUUUUg are studied with AMBER ff99 and CHARMM27 force fields. Experimentally, the UUCG hairpin is known to be significantly more stable than UUUU. In this study, the calculations using AMBER force field give a qualitatively correct description for the folding of two RNA hairpins, as the calculated PMF confirms the global stability of the folded structures and the resulting relative folding free energy is in quantitative agreement with the experimental result. The hairpin stabilities are also correctly differentiated by the more rapid molecular mechanics-Poisson Boltzmann-surface area approach, but the relative free energy estimated from this method is overestimated. The free energy profile shows that the native state basin and the unfolded state plateau are separated by a wide shoulder region, which samples a variety of native-like structures with frayed terminal basepair. The calculated PMF lacks major barriers that are expected near the transition regions, and this is attributed to the limitation of the 1-D reaction coordinate. The PMF results are compared with other studies of small RNA hairpins using kinetics method and coarse grained models. The two RNA hairpins described by CHARMM27 are significantly more deformable than those represented by AMBER. Compared with the AMBER results, the CHARMM27 calculated DeltaG(fold) for the UUUU tetraloop is in better agreement with the experimental results. However, the CHARMM27 calculation does not confirm the global stability of the experimental UUCG structure; instead, the extended conformations are predicted to be thermodynamically stable in solution. This finding is further supported by separate unrestrained CHARMM27 simulations, in which the UUCG hairpin unfolds spontaneously within 10 ns. The instability of the UUCG hairpin originates from the loop region, and propagates to the stem. The results of this study provide a molecular picture of RNA hairpin unfolding and reveal problems in the force field descriptions for the conformational energy of certain RNA hairpin.

Project description:Antimalarial drug resistance has thrown a spanner in the works of malaria elimination. New drugs are required for ancillary support of existing malaria control efforts. Plasmodium falciparum requires host glucose for survival and proliferation. On this basis, P. falciparum hexose transporter 1 (PfHT1) protein involved in hexose permeation is considered a potential drug target. In this study, we tested the antimalarial activity of some compounds against PfHT1 using computational techniques. We performed high throughput virtual screening of 21,352 small-molecule compounds against PfHT1. The stability of the lead compound complexes was evaluated via molecular dynamics (MD) simulation for 100 nanoseconds. We also investigated the pharmacodynamic, pharmacokinetic and physiological characteristics of the compounds in accordance with Lipinksi rules for drug-likeness to bind and inhibit PfHT1. Molecular docking and free binding energy analyses were carried out using Molecular Mechanics with Generalized Born and Surface Area (MMGBSA) solvation to determine the selectivity of the hit compounds for PfHT1 over the human glucose transporter (hGLUT1) orthologue. Five important PfHT1 inhibitors were identified: Hyperoside (CID5281643); avicularin (CID5490064); sylibin (CID5213); harpagoside (CID5481542) and quercetagetin (CID5281680). The compounds formed intermolecular interaction with the binding pocket of the PfHT1 target via conserved amino acid residues (Val314, Gly183, Thr49, Asn52, Gly183, Ser315, Ser317, and Asn48). The MMGBSA analysis of the complexes yielded high free binding energies. Four (CID5281643, CID5490064, CID5213, and CID5481542) of the identified compounds were found to be stable within the PfHT1 binding pocket throughout the 100 nanoseconds simulation run time. The four compounds demonstrated higher affinity for PfHT1 than the human major glucose transporter (hGLUT1). This investigation demonstrates the inhibition potential of sylibin, hyperoside, harpagoside, and avicularin against PfHT1 receptor. Robust preclinical investigations are required to validate the chemotherapeutic properties of the identified compounds.

Dataset Information

Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization.

Publications

Molecular Simulation of Chemical Reaction Equilibrium by Computationally Efficient Free Energy Minimization.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets