Metadynamics with Adaptive Gaussians.

Metadynamics is an established sampling method aimed at reconstructing the free-energy surface relative to a set of appropriately chosen collective variables. In standard metadynamics, the free-energy surface is filled by the addition of Gaussian potentials of preassigned and typically diagonal covariance. Asymptotically the free-energy surface is proportional to the bias deposited. Here, we consider the possibility of using Gaussians whose variance is adjusted on the fly to the local properties of the free-energy surface. We suggest two different prescriptions: one is based on the local diffusivity and the other on the local geometrical properties. We further examine the problem of extracting the free-energy surface when using adaptive Gaussians. We show that the standard relation between the bias and the free energy does not hold. In the limit of narrow Gaussians an explicit correction can be evaluated. In the general case, we propose to use instead a relation between bias and free energy borrowed from umbrella sampling. This relation holds for all kinds of incrementally deposited bias. We illustrate on the case of alanine dipeptide the advantage of using adaptive Gaussians in conjunction with the new free-energy estimator both in terms of accuracy and speed of convergence.

[1]  Juan J. de Pablo,et al.  Flux Tempered Metadynamics , 2011 .

[2]  Michele Parrinello,et al.  A self-learning algorithm for biased molecular dynamics , 2010, Proceedings of the National Academy of Sciences.

[3]  Gerhard Hummer,et al.  Coordinate-dependent diffusion in protein folding , 2009, Proceedings of the National Academy of Sciences.

[4]  Massimiliano Bonomi,et al.  Reconstructing the equilibrium Boltzmann distribution from well‐tempered metadynamics , 2009, J. Comput. Chem..

[5]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[6]  M. Parrinello,et al.  From four- to six-coordinated silica: Transformation pathways from metadynamics , 2007 .

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

[8]  Carsten Hartmann,et al.  Comment on two distinct notions of free energy , 2007 .

[9]  Francesco Luigi Gervasio,et al.  From A to B in free energy space. , 2007, The Journal of chemical physics.

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

[11]  M. Parrinello,et al.  Crystal structure transformations in SiO2 from classical and ab initio metadynamics , 2006, Nature materials.

[12]  G. Ciccotti,et al.  String method in collective variables: minimum free energy paths and isocommittor surfaces. , 2006, The Journal of chemical physics.

[13]  Simone Marsili,et al.  Self-healing umbrella sampling: a non-equilibrium approach for quantitative free energy calculations. , 2006, The journal of physical chemistry. B.

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

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

[16]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[17]  Eric Vanden-Eijnden,et al.  Transition state theory: variational formulation, dynamical corrections, and error estimates. , 2005, The Journal of chemical physics.

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

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

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

[21]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

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

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

[24]  Paul E. Smith,et al.  THE ALANINE DIPEPTIDE FREE ENERGY SURFACE IN SOLUTION , 1999 .

[25]  Amedeo Caflisch,et al.  Calculation of conformational transitions and barriers in solvated systems: Application to the alanine dipeptide in water , 1999 .

[26]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[27]  Christian Bartels,et al.  Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations , 1997, J. Comput. Chem..

[28]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[29]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

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

[31]  Charles L. Brooks,et al.  Conformational equilibrium in the alanine dipeptide in the gas phase and aqueous solution : a comparison of theoretical results , 1992 .

[32]  Charles L. Brooks,et al.  Reaction paths and free energy profiles for conformational transitions: An internal coordinate approach , 1991 .

[33]  S. Kearsley.  On the orthogonal transformation used for structural comparisons , 1989 .

[34]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[35]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[36]  G. N. Ramachandran,et al.  Stereochemistry of polypeptide chain configurations. , 1963, Journal of molecular biology.