Calculation of free‐energy differences from computer simulations of initial and final states

A class of simple expressions of increasing accuracy for the free‐energy difference between two states is derived based on numerical thermodynamic integration. The implementation of these formulas requires simulations of the initial and final (and possibly a few intermediate) states. They involve higher free‐energy derivatives at these states which are related to the moments of the probability distribution of the perturbation. Given a specified number of such derivatives, these integration formulas are optimal in the sense that they are exact to the highest possible order of free‐energy perturbation theory. The utility of this approach is illustrated for the hydration free energy of water. This problem provides a quite stringent test because the free energy is a highly nonlinear function of the charge so that even fourth order perturbation theory gives a very poor estimate of the free‐energy change. Our results should prove most useful for complex, computationally demanding problems where free‐energy differences arise primarily from changes in the electrostatic interactions (e.g., electron transfer,charging of ions, protonation of amino acids in proteins).

[1]  Kenneth M. Merz,et al.  Determination of pKas of ionizable groups in proteins: The pKa of Glu 7 and 35 in hen egg white lysozyme and Glu 106 in human carbonic anhydrase II , 1991 .

[2]  D. Chandler,et al.  Diabatic surfaces and the pathway for primary electron transfer in a photosynthetic reaction center , 1993 .

[3]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

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

[5]  Michele Migliore,et al.  Monte Carlo study of free energy of hydration for Li+, Na+, K+, F−, and Cl− with ab initio potentials , 1988 .

[6]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

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

[8]  D. D. Stancu,et al.  Quadrature formulas with simple Gaussian nodes and multiple fixed nodes , 1963 .

[9]  K. Sharp,et al.  On the calculation of pKas in proteins , 1993, Proteins.

[10]  M. Karplus,et al.  pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model. , 1990, Biochemistry.

[11]  Wilfred F. van Gunsteren,et al.  Predictions of free energy differences from a single simulation of the initial state , 1994 .

[12]  R. Marcus,et al.  Electron transfers in chemistry and biology , 1985 .

[13]  T. Straatsma,et al.  Free energy of ionic hydration: Analysis of a thermodynamic integration technique to evaluate free energy differences by molecular dynamics simulations , 1988 .

[14]  H. Berendsen,et al.  THERMODYNAMICS OF CAVITY FORMATION IN WATER - A MOLECULAR-DYNAMICS STUDY , 1982 .

[15]  Jiali Gao,et al.  Absolute free energy of solvation from Monte Carlo simulations using combined quantum and molecular mechanical potentials , 1992 .

[16]  R. A. Kuharski,et al.  Molecular model for aqueous ferrous–ferric electron transfer , 1988 .

[17]  G. King,et al.  Calculation of electrostatic free energy differences with a time-saving approximate method , 1993 .

[18]  A. Warshel,et al.  Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process. , 1986, Biochemistry.

[19]  Bhyravabhotla Jayaram,et al.  Free energy calculations of ion hydration: an analysis of the Born model in terms of microscopic simulations , 1989 .

[20]  A. Warshel,et al.  Computer simulations of electron-transfer reactions in solution and in photosynthetic reaction centers. , 1991, Annual review of physical chemistry.

[21]  K. Merz,et al.  Transferability of ion models , 1993 .

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

[23]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[24]  G. Hummer,et al.  Ion pair potentials-of-mean-force in water , 1994, chem-ph/9404001.

[25]  P. Beroza,et al.  Electrostatic calculations of amino acid titration and electron transfer, Q-AQB-->QAQ-B, in the reaction center. , 1995, Biophysical journal.

[26]  A. Warshel,et al.  Microscopic examination of free-energy relationships for electron transfer in polar solvents , 1987 .

[27]  Steven W. Rick,et al.  The Aqueous Solvation of Water: A Comparison of Continuum Methods with Molecular Dynamics , 1994 .

[28]  Ronald M. Levy,et al.  Gaussian fluctuation formula for electrostatic free‐energy changes in solution , 1991 .

[29]  Jaroslav Kautsky,et al.  Algorithm 655: IQPACK: FORTRAN subroutines for the weights of interpolatory quadratures , 1987, TOMS.

[30]  A. Warshel,et al.  Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603. , 1992, Protein engineering.

[31]  Huan‐Xiang Zhou,et al.  Microscopic formulation of Marcus’ theory of electron transfer , 1995 .

[32]  Gerhard Hummer,et al.  The hydration free energy of water , 1995 .

[33]  M Karplus,et al.  Chromophore-protein interactions and the function of the photosynthetic reaction center: a molecular dynamics study. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[35]  Minoru Saito,et al.  Molecular dynamics/free energy study of a protein in solution with all degrees of freedom and long-range Coulomb interactions , 1995 .

[36]  R. Levy,et al.  Intrinsic pKas of ionizable residues in proteins: An explicit solvent calculation for lysozyme , 1994, Proteins.