Convergence of Folding Free Energy Landscapes Via Application of Enhanced Sampling Methods in a Distributed Computing Environment

Overview

Journal J Chem Phys

Publisher American Institute of Physics

Specialties Biophysics
Chemistry

Date 2008 Jun 3

PMID 18513049

Citations 28

Authors

Xuhui Huang

Gregory R Bowman

Vijay S Pande

Affiliations

Soon will be listed here.

Abstract

We have implemented the serial replica exchange method (SREM) and simulated tempering (ST) enhanced sampling algorithms in a global distributed computing environment. Here we examine the helix-coil transition of a 21 residue alpha-helical peptide in explicit solvent. For ST, we demonstrate the efficacy of a new method for determining initial weights allowing the system to perform a random walk in temperature space based on short trial simulations. These weights are updated throughout the production simulation by an adaptive weighting method. We give a detailed comparison of SREM, ST, as well as standard MD and find that SREM and ST give equivalent results in reasonable agreement with experimental data. In addition, we find that both enhanced sampling methods are much more efficient than standard MD simulations. The melting temperature of the Fs peptide with the AMBER99phi potential was calculated to be about 310 K, which is in reasonable agreement with the experimental value of 334 K. We also discuss other temperature dependent properties of the helix-coil transition. Although ST has certain advantages over SREM, both SREM and ST are shown to be powerful methods via distributed computing and will be applied extensively in future studies of complex bimolecular systems.

Citing Articles

Efficient Sampling of Cavity Hydration in Proteins with Nonequilibrium Grand Canonical Monte Carlo and Polarizable Force Fields.

Deng J, Cui Q J Chem Theory Comput. 2024; 20(5):1897-1911.

PMID: 38417108 PMC: 11663258. DOI: 10.1021/acs.jctc.4c00013.

Continuous millisecond conformational cycle of a DEAH box helicase reveals control of domain motions by atomic-scale transitions.

Becker R, Hub J Commun Biol. 2023; 6(1):379.

PMID: 37029280 PMC: 10082070. DOI: 10.1038/s42003-023-04751-z.

Directed kinetic transition network model.

Zhou H, Wang F, Bennett D, Tao P J Chem Phys. 2019; 151(14):144112.

PMID: 31615261 PMC: 6800283. DOI: 10.1063/1.5110896.

Massive-Scale Binding Free Energy Simulations of HIV Integrase Complexes Using Asynchronous Replica Exchange Framework Implemented on the IBM WCG Distributed Network.

Xia J, Flynn W, Gallicchio E, Uplinger K, Armstrong J, Forli S J Chem Inf Model. 2019; 59(4):1382-1397.

PMID: 30758197 PMC: 6496938. DOI: 10.1021/acs.jcim.8b00817.

Predicting the Conformational Variability of Abl Tyrosine Kinase using Molecular Dynamics Simulations and Markov State Models.

Meng Y, Gao C, Clawson D, Atwell S, Russell M, Vieth M J Chem Theory Comput. 2018; 14(5):2721-2732.

PMID: 29474075 PMC: 6317529. DOI: 10.1021/acs.jctc.7b01170.

References

Snow C, Nguyen H, Pande V, Gruebele M . Absolute comparison of simulated and experimental protein-folding dynamics. Nature. 2002; 420(6911):102-6. DOI: 10.1038/nature01160. View

Williams S, Causgrove T, Gilmanshin R, Fang K, Callender R, Woodruff W . Fast events in protein folding: helix melting and formation in a small peptide. Biochemistry. 1996; 35(3):691-7. DOI: 10.1021/bi952217p. View

Schug A, Herges T, Wenzel W . All-atom folding of the three-helix HIV accessory protein with an adaptive parallel tempering method. Proteins. 2004; 57(4):792-8. DOI: 10.1002/prot.20290. View

Garcia A, Sanbonmatsu K . Alpha-helical stabilization by side chain shielding of backbone hydrogen bonds. Proc Natl Acad Sci U S A. 2002; 99(5):2782-7. PMC: 122425. DOI: 10.1073/pnas.042496899. View

Zhou R . Trp-cage: folding free energy landscape in explicit water. Proc Natl Acad Sci U S A. 2003; 100(23):13280-5. PMC: 263783. DOI: 10.1073/pnas.2233312100. View

Hansmann U, Okamoto Y . New Monte Carlo algorithms for protein folding. Curr Opin Struct Biol. 1999; 9(2):177-83. DOI: 10.1016/S0959-440X(99)80025-6. View

Chodera J, Swope W, Pitera J, Seok C, Dill K . Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. J Chem Theory Comput. 2015; 3(1):26-41. DOI: 10.1021/ct0502864. View

Swendsen , Wang . Replica Monte Carlo simulation of spin glasses. Phys Rev Lett. 1986; 57(21):2607-2609. DOI: 10.1103/PhysRevLett.57.2607. View

Gallicchio E, Levy R, Parashar M . Asynchronous replica exchange for molecular simulations. J Comput Chem. 2007; 29(5):788-94. PMC: 2977925. DOI: 10.1002/jcc.20839. View

10.

Wang F, Landau D . Efficient, multiple-range random walk algorithm to calculate the density of states. Phys Rev Lett. 2001; 86(10):2050-3. DOI: 10.1103/PhysRevLett.86.2050. View

11.

Zhou R, Berne B, Germain R . The free energy landscape for beta hairpin folding in explicit water. Proc Natl Acad Sci U S A. 2001; 98(26):14931-6. PMC: 64961. DOI: 10.1073/pnas.201543998. View

12.

Shirts M, Pande V . COMPUTING: Screen Savers of the World Unite!. Science. 2007; 290(5498):1903-4. DOI: 10.1126/science.290.5498.1903. View

13.

Jayachandran G, Vishal V, Garcia A, Pande V . Local structure formation in simulations of two small proteins. J Struct Biol. 2006; 157(3):491-9. PMC: 1931413. DOI: 10.1016/j.jsb.2006.10.001. View

14.

Hoover . Canonical dynamics: Equilibrium phase-space distributions. Phys Rev A Gen Phys. 1985; 31(3):1695-1697. DOI: 10.1103/physreva.31.1695. View

15.

Nymeyer H, Garcia A . Simulation of the folding equilibrium of alpha-helical peptides: a comparison of the generalized Born approximation with explicit solvent. Proc Natl Acad Sci U S A. 2003; 100(24):13934-9. PMC: 283524. DOI: 10.1073/pnas.2232868100. View

16.

Sorin E, Pande V . Nanotube confinement denatures protein helices. J Am Chem Soc. 2006; 128(19):6316-7. DOI: 10.1021/ja060917j. View

17.

Zheng W, Andrec M, Gallicchio E, Levy R . Simulating replica exchange simulations of protein folding with a kinetic network model. Proc Natl Acad Sci U S A. 2007; 104(39):15340-5. PMC: 2000486. DOI: 10.1073/pnas.0704418104. View

18.

Jayachandran G, Vishal V, Pande V . Using massively parallel simulation and Markovian models to study protein folding: examining the dynamics of the villin headpiece. J Chem Phys. 2006; 124(16):164902. DOI: 10.1063/1.2186317. View

19.

Lin C, Hu C, Hansmann U . Parallel tempering simulations of HP-36. Proteins. 2003; 52(3):436-45. DOI: 10.1002/prot.10351. View

20.

Li H, Fajer M, Yang W . Simulated scaling method for localized enhanced sampling and simultaneous "alchemical" free energy simulations: a general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations. J Chem Phys. 2007; 126(2):024106. DOI: 10.1063/1.2424700. View