Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods

A challenge in free energy calculation for complex molecular systems by computer simulation is to obtain a reliable estimate within feasible computational time. In this study, we suggest an answer to this challenge by exploring a simple method, overlap sampling (OS), for producing reliable free‐energy results in an efficient way. The formalism of the OS method is based on ensuring sampling of important overlapping phase space during perturbation calculations. This technique samples both forward and reverse free energy perturbation (FEP) to improve the free‐energy calculation. It considers the asymmetry of the FEP calculation and features an ability to optimize both the precision and the accuracy of the measurement without affecting the simulation process itself. The OS method is tested at two optimization levels: no optimization (simple OS), and full optimization (equivalent to Bennett's method), and compared to conventional FEP techniques, including the widely used direct FEP averaging method, on three alchemical mutation systems: (a) an anion transformation in water solution, (b) mutation between methanol and ethane, and (c) alchemical change of an adenosine molecule. It is consistently shown that the reliability of free‐energy estimates can be greatly improved using the OS techniques at both optimization levels, while the performance of Bennett's method is particularly striking. In addition, the efficiency of a calculation can be significantly improved because the method is able to (a) converge to the right answer quickly, and (b) work for large perturbations. The basic two‐stage OS method can be extended to admit additional stages, if needed. We suggest that the OS method can be used as a general perturbation technique for computing free energy differences in molecular simulations. © 2003 Wiley Periodicals, Inc. J Comput Chem 1: 28–39, 2004

[1]  T. Woolf,et al.  Staging Is More Important than Perturbation Method for Computation of Enthalpy and Entropy Changes in Complex Systems , 2003 .

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

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

[4]  J. G. Powles,et al.  The computation of the chemical potential of a fluid , 1982 .

[5]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

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

[7]  David A. Kofke,et al.  Accuracy of free-energy perturbation calculations in molecular simulation. II. Heuristics , 2001 .

[8]  J. Valverde.  Molecular Modelling: Principles and Applications , 2001 .

[9]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[10]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

[11]  W. L. Jorgensen,et al.  Monte Carlo simulation of differences in free energies of hydration , 1985 .

[12]  Exploration of the phase space of molecular systems: Assessment of established and new methods , 1992 .

[13]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

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

[15]  P. Kollman,et al.  Water–water and water–ion potential functions including terms for many body effects , 1985 .

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

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

[18]  Robert H. Wood,et al.  Systematic errors in free energy perturbation calculations due to a finite sample of configuration space: sample-size hysteresis , 1991 .

[19]  Peter A. Kollman,et al.  A new method for carrying out free energy perturbation calculations: Dynamically modified windows , 1989 .

[20]  P. Schleyer.  Encyclopedia of computational chemistry , 1998 .

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

[22]  W. L. Jorgensen,et al.  Origin of the strong binding of adenine to a molecular tweezer , 1990 .

[23]  M. Williams.  Adenosine and Adenosine Receptors , 1990, The Receptors.

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

[25]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[26]  Peter A. Kollman,et al.  Free energy calculation methods: A theoretical and empirical comparison of numerical errors and a new method qualitative estimates of free energy changes , 1997, J. Comput. Chem..

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

[28]  William H. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[29]  C. Chipot,et al.  Proton Transfer in the Mono- and the Dihydrated Complexes of HF and HCl: An MP2/6-31+G ab Initio Study in the Self-Consistent Reaction Field Model of Solvation , 1994 .

[30]  David A. Kofke,et al.  Accuracy of free-energy perturbation calculations in molecular simulation. I. Modeling , 2001 .

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

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

[33]  Christopher A. Reynolds,et al.  Free energy calculations in molecular biophysics , 1992 .

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

[35]  Frederick E. Petry,et al.  Principles and Applications , 1997 .

[36]  David A. Pearlman,et al.  Free energy derivatives: A new method for probing the convergence problem in free energy calculations , 1994, J. Comput. Chem..

[37]  Daniel M. Zuckerman,et al.  Overcoming finite-sampling errors in fast-switching free-energy estimates: extrapolative analysis of a molecular system , 2002 .

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

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

[40]  C. Brooks.  Computer simulation of liquids , 1989 .

[41]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

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

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

[44]  Free energies of hydration using restrained electrostatic potential derived charges via free energy perturbations and linear response , 1999 .

[45]  K. Binder,et al.  Monte Carlo and molecular dynamics of condensed matter systems : euroconference on computer simulation in condensed matter physics and chemistry, Como, 3-28 July 1995 , 1996 .

[46]  P. Kollman,et al.  Encyclopedia of computational chemistry , 1998 .