Synergistic approach to improve "alchemical" free energy calculation in rugged energy surface.

The authors present an integrated approach to "alchemical" free energy simulation, which permits efficient calculation of the free energy difference on rugged energy surface. The method is designed to obtain efficient canonical sampling for rapid free energy convergence. The proposal is motivated by the insight that both the exchange efficiency in the presently designed dual-topology alchemical Hamiltonian replica exchange method (HREM), and the confidence of the free energy determination using the overlap histogramming method, depend on the same criterion, viz., the overlaps of the energy difference histograms between all pairs of neighboring states. Hence, integrating these two techniques can produce a joint solution to the problems of the free energy convergence and conformational sampling in the free energy simulations, in which lambda parameter plays two roles to simultaneously facilitate the conformational sampling and improve the phase space overlap for the free energy determination. Specifically, in contrast with other alchemical HREM based free energy simulation methods, the dual-topology approach can ensure robust conformational sampling. Due to these features (a synergistic solution to the free energy convergence and canonical sampling, and the improvement of the sampling efficiency with the dual-topology treatment), the present approach, as demonstrated in the model studies of the authors, is highly efficient in obtaining accurate free energy differences, especially for the systems with rough energy landscapes.

[1]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[2]  L. Scriven,et al.  Efficient molecular simulation of chemical potentials , 1989 .

[3]  B. Berne,et al.  Novel methods of sampling phase space in the simulation of biological systems. , 1997, Current opinion in structural biology.

[4]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .

[5]  K. Tai Conformational sampling for the impatient. , 2004, Biophysical chemistry.

[6]  Stefan Boresch,et al.  The Role of Bonded Terms in Free Energy Simulations. 2. Calculation of Their Influence on Free Energy Differences of Solvation , 1999 .

[7]  Ulrich H E Hansmann,et al.  Efficient sampling of protein structures by model hopping. , 2005, Physical review letters.

[8]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[9]  Wei Yang,et al.  Generalized ensembles serve to improve the convergence of free energy simulations , 2003 .

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

[11]  Thomas Simonson,et al.  Free energy simulations come of age: protein-ligand recognition. , 2002, Accounts of chemical research.

[12]  Jonathan W. Essex,et al.  The development of replica-exchange-based free-energy methods , 2003 .

[13]  David A. Pearlman,et al.  A Comparison of Alternative Approaches to Free Energy Calculations , 1994 .

[14]  D. Zuckerman,et al.  Single-ensemble nonequilibrium path-sampling estimates of free energy differences. , 2004, The Journal of chemical physics.

[15]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[16]  Yuko Okamoto,et al.  Generalized-ensemble algorithms: enhanced sampling techniques for Monte Carlo and molecular dynamics simulations. , 2003, Journal of molecular graphics & modelling.

[17]  Thomas Simonson,et al.  Free energy of particle insertion , 1993 .

[18]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

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

[20]  T. Schlick,et al.  Biomolecular free energy profiles by a shooting/umbrella sampling protocol, "BOLAS". , 2004, The Journal of chemical physics.

[21]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[22]  Stefan Boresch,et al.  THE ROLE OF BONDED TERMS IN FREE ENERGY SIMULATIONS : 1. THEORETICAL ANALYSIS , 1999 .

[23]  Wei Yang,et al.  Conformational sampling via a self-regulating effective energy surface. , 2006, The Journal of chemical physics.

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

[25]  Régis Pomès,et al.  Absolute free energy calculations by thermodynamic integration in four spatial dimensions. , 2005, The Journal of chemical physics.

[26]  Michael R. Shirts,et al.  Comparison of efficiency and bias of free energies computed by exponential averaging, the Bennett acceptance ratio, and thermodynamic integration. , 2005, The Journal of chemical physics.

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

[28]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[29]  Régis Pomès,et al.  Enhancing the accuracy, the efficiency and the scope of free energy simulations. , 2005, Current opinion in structural biology.

[30]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[31]  Alan E. Mark,et al.  Calculation of Relative Free-Energy Via Indirect Pathways , 1991 .

[32]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[33]  K. Shing,et al.  The chemical potential in dense fluids and fluid mixtures via computer simulation , 1982 .

[34]  Thomas B Woolf,et al.  Using overlap and funnel sampling to obtain accurate free energies from nonequilibrium work measurements. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Martin Karplus,et al.  Optimal estimates of free energies from multistate nonequilibrium work data. , 2006, Physical review letters.

[36]  Peter T. Cummings,et al.  Quantitative comparison and optimization of methods for evaluating the chemical potential by molecular simulation , 1997 .

[37]  Charles L. Brooks,et al.  Efficient Sampling of Ligand Orientations and Conformations in Free Energy Calculations Using the λ-Dynamics Method , 2000 .

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

[39]  Anthony K. Felts,et al.  Temperature weighted histogram analysis method, replica exchange, and transition paths. , 2005, The journal of physical chemistry. B.

[40]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[41]  Stefan Boresch,et al.  Alchemical free energy calculations and multiple conformational substates. , 2005, The Journal of chemical physics.

[42]  Sidney Yip,et al.  Optimized free-energy evaluation using a single reversible-scaling simulation , 1999 .

[43]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

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

[45]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[46]  Seokmin Shin,et al.  Replica-exchange method using the generalized effective potential. , 2003, Physical review letters.

[47]  Ron Elber,et al.  Long-timescale simulation methods. , 2005, Current opinion in structural biology.

[48]  Chris Oostenbrink,et al.  Single‐step perturbations to calculate free energy differences from unphysical reference states: Limits on size, flexibility, and character , 2003, J. Comput. Chem..

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

[50]  Eric Darve,et al.  Assessing the efficiency of free energy calculation methods. , 2004, The Journal of chemical physics.

[51]  G. Bester,et al.  Interlayer contraction in MgB2 upon replacement of Mg by Al: Effect of the covalent bond energy , 2005 .

[52]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[53]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[54]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[55]  I. Andricioaei,et al.  On Monte Carlo and molecular dynamics methods inspired by Tsallis statistics: Methodology, optimization, and application to atomic clusters , 1997 .

[56]  H. Scheraga,et al.  Global optimization of clusters, crystals, and biomolecules. , 1999, Science.

[57]  T. Woolf,et al.  Using non-Gaussian density functional fits to improve relative free energy calculations. , 2005, The Journal of chemical physics.

[58]  G. Crooks Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.