Efficient Solvation Free Energy Calculations of Amino Acid Analogs by Expanded Ensemble Molecular Simulation.

We present an efficient, automated expanded ensemble method to calculate the residual chemical potential or solvation free energy by molecular dynamics simulation. The methodology is validated by computing the residual chemical potential of 13 amino acid analogs in water at 300 K and 1 bar and comparing to reference simulation data. Overall agreement is good, with the methodology of the present study reaching limiting precisions of less than 0.1 kBT in half of the total simulation time of the reference simulation study which utilized Bennett's acceptance ratio method. The apparent difference in the efficiencies is a result of the inherent advantages of the expanded ensemble method, which creates an improved decorrelation of simulation data and improves the sampling of the important regions of the configurational phase space of each subensemble. The present adaptation utilizes histograms of proposed transition energies collected throughout the entire simulation, to make extremely precise calculations of the relative free energy between neighboring subensembles.

[1]  Sven P. Jacobsson,et al.  Solubility of Organic Compounds in Water/Octanol Systems. A Expanded Ensemble Molecular Dynamics Simulation Study of log P Parameters , 2001 .

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

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

[4]  E. Maginn,et al.  Monte Carlo simulations of gas solubility in the ionic liquid 1-n-butyl-3-methylimidazolium hexafluorophosphate. , 2005, The journal of physical chemistry. B.

[5]  P M Cullis,et al.  Affinities of amino acid side chains for solvent water. , 1981, Biochemistry.

[6]  David L Mobley,et al.  Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent. , 2007, The journal of physical chemistry. B.

[7]  William C. Swope,et al.  The role of long ranged forces in determining the structure and properties of liquid water , 1983 .

[8]  F. Escobedo,et al.  Expanded ensemble and replica exchange methods for simulation of protein-like systems , 2003 .

[9]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[10]  Andrew S. Paluch,et al.  Comparing the Use of Gibbs Ensemble and Grand-Canonical Transition-Matrix Monte Carlo Methods to Determine Phase Equilibria , 2008 .

[11]  Alexander P. Lyubartsev,et al.  Solvation free energies of methane and alkali halide ion pairs: An expanded ensemble molecular dynamics simulation study , 1998 .

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

[13]  Fernando A Escobedo,et al.  Variance minimization of free energy estimates from optimized expanded ensembles. , 2008, The journal of physical chemistry. B.

[14]  Density-of-states Monte Carlo method for simulation of fluids , 2002, cond-mat/0201470.

[15]  Jeffrey R Errington.  Evaluating surface tension using grand-canonical transition-matrix Monte Carlo simulation and finite-size scaling. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Alexander P. Lyubartsev,et al.  Determination of Free Energy from Chemical Potentials: Application of the Expanded Ensemble Method , 1996 .

[17]  Jaeeon Chang,et al.  The calculation of chemical potential of organic solutes in dense liquid phases by using expanded ensemble Monte Carlo simulations. , 2009, The Journal of chemical physics.

[18]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[19]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[20]  Jeffrey R. Errington,et al.  Direct calculation of liquid–vapor phase equilibria from transition matrix Monte Carlo simulation , 2003 .

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

[22]  F. Escobedo,et al.  On the use of Bennett's acceptance ratio method in multi-canonical-type simulations. , 2004, The Journal of chemical physics.

[23]  Alexander P. Lyubartsev,et al.  Free energy calculations for Lennard-Jones systems and water using the expanded ensemble method A Monte Carlo and molecular dynamics simulation study , 1994 .

[24]  D. Landau,et al.  Determining the density of states for classical statistical models: a random walk algorithm to produce a flat histogram. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  J. Errington.  Prewetting transitions for a model argon on solid carbon dioxide system. , 2004, Langmuir : the ACS journal of surfaces and colloids.

[26]  Di Wu,et al.  Phase-space overlap measures. II. Design and implementation of staging methods for free-energy calculations. , 2005, The Journal of chemical physics.

[27]  Q. Yan,et al.  Fast calculation of the density of states of a fluid by Monte Carlo simulations. , 2003, Physical review letters.

[28]  B. Hess,et al.  Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models. , 2006, The journal of physical chemistry. B.

[29]  B. Efron.  Nonparametric estimates of standard error: The jackknife, the bootstrap and other methods , 1981 .

[30]  Di Wu,et al.  Phase-space overlap measures. I. Fail-safe bias detection in free energies calculated by molecular simulation. , 2005, The Journal of chemical physics.

[31]  C. Abreu,et al.  On the use of transition matrix methods with extended ensembles. , 2006, The Journal of chemical physics.

[32]  D. N. Card,et al.  Monte Carlo Estimation of the Free Energy by Multistage Sampling , 1972 .

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

[34]  K. Dill,et al.  The protein folding problem. , 1993, Annual review of biophysics.

[35]  Abraham M Lenhoff,et al.  Solvation free energy of amino acids and side-chain analogues. , 2007, The journal of physical chemistry. B.

[36]  Aatto Laaksonen,et al.  Determination of solvation free energies by adaptive expanded ensemble molecular dynamics. , 2004, The Journal of chemical physics.

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

[38]  C. Anfinsen.  Principles that govern the folding of protein chains. , 1973, Science.

[39]  J. Errington,et al.  Determination of Henry's law constants through transition matrix Monte Carlo simulation , 2005 .

[40]  Michael R. Shirts,et al.  Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. , 2005, The Journal of chemical physics.

[41]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[42]  Michael R. Shirts,et al.  Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. , 2003, Physical review letters.

[43]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[44]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[45]  V. K. Shen,et al.  Metastability and instability in the Lennard-Jones fluid investigated by transition-matrix Monte Carlo , 2004 .

[46]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[47]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

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

[49]  Peter T. Cummings,et al.  Precision and accuracy of staged free-energy perturbation methods for computing the chemical potential by molecular simulation , 1998 .

[50]  David L Mobley,et al.  Nonlinear scaling schemes for Lennard-Jones interactions in free energy calculations. , 2007, The Journal of chemical physics.

[51]  P. G. Debenedetti,et al.  An improved Monte Carlo method for direct calculation of the density of states , 2003 .

[52]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[53]  Alexander P. Lyubartsev,et al.  M.DynaMix – a scalable portable parallel MD simulation package for arbitrary molecular mixtures , 2000 .

[54]  Andrew S. Paluch,et al.  A method for computing the solubility limit of solids: application to sodium chloride in water and alcohols. , 2010, The Journal of chemical physics.

[55]  Athanassios Z Panagiotopoulos,et al.  Generalization of the Wang-Landau method for off-lattice simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[56]  R. R. Picard,et al.  Canonical transition probabilities for adaptive Metropolis simulation , 1999 .

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

[58]  R. R. Picard,et al.  Monte Carlo Transition Dynamics and Variance Reduction , 2000 .

[59]  D. Kofke,et al.  Optimal intermediates in staged free energy calculations , 1999 .

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

[61]  David A. Kofke,et al.  Appropriate methods to combine forward and reverse free-energy perturbation averages , 2003 .

[62]  Charles L. Brooks,et al.  λ‐dynamics: A new approach to free energy calculations , 1996 .

[63]  B. Efron.  Computers and the Theory of Statistics: Thinking the Unthinkable , 1979 .

[64]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[65]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

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

[67]  David L Mobley,et al.  Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. , 2009, Journal of chemical theory and computation.