An algorithm to find minimum free-energy paths using umbrella integration.

The calculation of free-energy barriers by umbrella sampling and many other methods is hampered by the necessity for an a priori choice of the reaction coordinate along which to sample. We avoid this problem by providing a method to search for saddle points on the free-energy surface in many coordinates. The necessary gradients and Hessians of the free energy are obtained by multidimensional umbrella integration. We construct the minimum free-energy path by following the gradient down to minima on the free-energy surface. The change of free energy along the path is obtained by integrating out all coordinates orthogonal to the path. While we expect the method to be applicable to large systems, we test it on the alanine dipeptide in vacuum. The minima, transition states, and free-energy barriers agree well with those obtained previously with other methods.

[1]  Jerry M. Parks,et al.  Quantum mechanics/molecular mechanics minimum free-energy path for accurate reaction energetics in solution and enzymes: sequential sampling and optimization on the potential of mean force surface. , 2008, The Journal of chemical physics.

[2]  Weitao Yang,et al.  Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface , 2000 .

[3]  P. Tavan,et al.  Ligand Binding: Molecular Mechanics Calculation of the Streptavidin-Biotin Rupture Force , 1996, Science.

[4]  M. Tuckerman,et al.  IN CLASSICAL AND QUANTUM DYNAMICS IN CONDENSED PHASE SIMULATIONS , 1998 .

[5]  A. Bell,et al.  Efficient methods for finding transition states in chemical reactions: comparison of improved dimer method and partitioned rational function optimization method. , 2005, The Journal of chemical physics.

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

[7]  I. Todorov,et al.  The DL_POLY molecular dynamics package , 2005 .

[8]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[9]  G. Henkelman,et al.  A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives , 1999 .

[10]  H. Jónsson,et al.  Nudged elastic band method for finding minimum energy paths of transitions , 1998 .

[11]  K. Fukui.  Variational principles in a chemical reaction , 2009 .

[12]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[13]  Johannes Kästner.  Umbrella integration with higher-order correction terms. , 2012, The Journal of chemical physics.

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

[15]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[16]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[17]  Karunesh Arora,et al.  Harmonic Fourier beads method for studying rare events on rugged energy surfaces. , 2006, The Journal of chemical physics.

[18]  K Schulten,et al.  Reconstructing potential energy functions from simulated force-induced unbinding processes. , 1997, Biophysical journal.

[19]  Weitao Yang,et al.  QM/MM Minimum Free Energy Path: Methodology and Application to Triosephosphate Isomerase. , 2007, Journal of chemical theory and computation.

[20]  Jan Hermans,et al.  Simple analysis of noise and hysteresis in (slow-growth) free energy simulations , 1991 .

[21]  T. Straatsma,et al.  Free energy of hydrophobic hydration: A molecular dynamics study of noble gases in water , 1986 .

[22]  E. Vanden-Eijnden,et al.  On-the-fly string method for minimum free energy paths calculation , 2007 .

[23]  Walter Thiel,et al.  Analysis of the statistical error in umbrella sampling simulations by umbrella integration. , 2006, The Journal of chemical physics.

[24]  Paul Sherwood,et al.  Superlinearly converging dimer method for transition state search. , 2008, The Journal of chemical physics.

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

[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]  G. Henkelman,et al.  Comparison of methods for finding saddle points without knowledge of the final states. , 2004, The Journal of chemical physics.

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

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

[30]  G. Hummer.  Fast-growth thermodynamic integration: Error and efficiency analysis , 2001 .

[31]  E. Evans,et al.  Dynamic strength of molecular adhesion bonds. , 1997, Biophysical journal.

[32]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[33]  K. Fukui.  The path of chemical reactions - the IRC approach , 1981 .

[34]  K. Schulten,et al.  Molecular dynamics study of unbinding of the avidin-biotin complex. , 1997, Biophysical journal.

[35]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[36]  Thomas B Woolf,et al.  Theory of a systematic computational error in free energy differences. , 2002, Physical review letters.

[37]  S. C. Rogers,et al.  QUASI: A general purpose implementation of the QM/MM approach and its application to problems in catalysis , 2003 .

[38]  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.

[39]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[40]  P. Fleurat‐Lessard,et al.  Tracing the minimum-energy path on the free-energy surface. , 2005, The Journal of chemical physics.

[41]  Walter Thiel,et al.  Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: "Umbrella integration". , 2005, The Journal of chemical physics.

[42]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[43]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[44]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[45]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

[46]  Johannes Kästner,et al.  Umbrella integration in two or more reaction coordinates. , 2009, The Journal of chemical physics.