Adaptively biased molecular dynamics for free energy calculations.

We present an adaptively biased molecular dynamics (ABMD) method for the computation of the free energy surface of a reaction coordinate using nonequilibrium dynamics. The ABMD method belongs to the general category of umbrella sampling methods with an evolving biasing potential and is inspired by the metadynamics method. The ABMD method has several useful features, including a small number of control parameters and an O(t) numerical cost with molecular dynamics time t. The ABMD method naturally allows for extensions based on multiple walkers and replica exchange, where different replicas can have different temperatures and/or collective variables. This is beneficial not only in terms of the speed and accuracy of a calculation, but also in terms of the amount of useful information that may be obtained from a given simulation. The workings of the ABMD method are illustrated via a study of the folding of the Ace-GGPGGG-Nme peptide in a gaseous and solvated environment.

[1]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[2]  A. Laio,et al.  Equilibrium free energies from nonequilibrium metadynamics. , 2006, Physical Review Letters.

[3]  Young K. Truong,et al.  Polynomial splines and their tensor products in extended linear modeling: 1994 Wald memorial lecture , 1997 .

[4]  Deprotonation of solvated formic acid: Car-Parrinello and metadynamics simulations. , 2006, The journal of physical chemistry. B.

[5]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[7]  Alessandro Laio,et al.  Azulene-to-naphthalene rearrangement: the Car-Parrinello metadynamics method explores various reaction mechanisms. , 2004, Chemphyschem : a European journal of chemical physics and physical chemistry.

[8]  A. Laio,et al.  Flexible docking in solution using metadynamics. , 2005, Journal of the American Chemical Society.

[9]  B. Silverman.  Density estimation for statistics and data analysis , 1986 .

[10]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[11]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[12]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[13]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[14]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[15]  Gabriel Stoltz,et al.  Computation of free energy profiles with parallel adaptive dynamics. , 2007, The Journal of chemical physics.

[16]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[17]  A. Laio,et al.  Assessing the accuracy of metadynamics. , 2005, The journal of physical chemistry. B.

[18]  M. Parrinello,et al.  Proton transfer in heterocycle crystals. , 2004, Physical review letters.

[19]  Andrew E. Torda,et al.  Local elevation: A method for improving the searching properties of molecular dynamics simulation , 1994, J. Comput. Aided Mol. Des..

[20]  A. Laio,et al.  A bias-exchange approach to protein folding. , 2007, The journal of physical chemistry. B.

[21]  Volodymyr Babin,et al.  The free energy landscape of small peptides as obtained from metadynamics with umbrella sampling corrections. , 2006, The Journal of chemical physics.

[22]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[23]  T. Ikeda,et al.  Hydration structure of Y3+ and La3+ compared: an application of metadynamics. , 2005, The Journal of chemical physics.

[24]  Alessandro Laio,et al.  Efficient exploration of reactive potential energy surfaces using Car-Parrinello molecular dynamics. , 2003, Physical review letters.

[25]  M. Parrinello,et al.  Ab Initio Study of Dehydroxylation−Carbonation Reaction on Brucite Surface , 2004 .

[26]  Exploring intramolecular reactions in complex systems with metadynamics: the case of the malonate anions. , 2005, The journal of physical chemistry. A.

[27]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[28]  A. Laio,et al.  Efficient reconstruction of complex free energy landscapes by multiple walkers metadynamics. , 2006, The journal of physical chemistry. B.

[29]  A. Preusser,et al.  Algorithm 671: FARB-E-2D: fill area with bicubics on rectangles—a contour plot program , 1989, TOMS.

[30]  Alessandro Laio,et al.  A minimum free energy reaction path for the E2 reaction between fluoro ethane and a fluoride ion. , 2004, Journal of the American Chemical Society.

[31]  Alessandro Laio,et al.  Microscopic Mechanism of Antibiotics Translocation through a Porin. , 2004, Biophysical journal.