Bennett's acceptance ratio and histogram analysis methods enhanced by umbrella sampling along a reaction coordinate in configurational space.

Free energy perturbation, a method for computing the free energy difference between two states, is often combined with non-Boltzmann biased sampling techniques in order to accelerate the convergence of free energy calculations. Here we present a new extension of the Bennett acceptance ratio (BAR) method by combining it with umbrella sampling (US) along a reaction coordinate in configurational space. In this approach, which we call Bennett acceptance ratio with umbrella sampling (BAR-US), the conditional histogram of energy difference (a mapping of the 3N-dimensional configurational space via a reaction coordinate onto 1D energy difference space) is weighted for marginalization with the associated population density along a reaction coordinate computed by US. This procedure produces marginal histograms of energy difference, from forward and backward simulations, with higher overlap in energy difference space, rendering free energy difference estimations using BAR statistically more reliable. In addition to BAR-US, two histogram analysis methods, termed Bennett overlapping histograms with US (BOH-US) and Bennett-Hummer (linear) least square with US (BHLS-US), are employed as consistency and convergence checks for free energy difference estimation by BAR-US. The proposed methods (BAR-US, BOH-US, and BHLS-US) are applied to a 1-dimensional asymmetric model potential, as has been used previously to test free energy calculations from non-equilibrium processes. We then consider the more stringent test of a 1-dimensional strongly (but linearly) shifted harmonic oscillator, which exhibits no overlap between two states when sampled using unbiased Brownian dynamics. We find that the efficiency of the proposed methods is enhanced over the original Bennett's methods (BAR, BOH, and BHLS) through fast uniform sampling of energy difference space via US in configurational space. We apply the proposed methods to the calculation of the electrostatic contribution to the absolute solvation free energy (excess chemical potential) of water. We then address the controversial issue of ion selectivity in the K(+) ion channel, KcsA. We have calculated the relative binding affinity of K(+) over Na(+) within a binding site of the KcsA channel for which different, though adjacent, K(+) and Na(+) configurations exist, ideally suited to these US-enhanced methods. Our studies demonstrate that the significant improvements in free energy calculations obtained using the proposed methods can have serious consequences for elucidating biological mechanisms and for the interpretation of experimental data.

[1]  K. Dill,et al.  Binding of small-molecule ligands to proteins: "what you see" is not always "what you get". , 2009, Structure.

[2]  G. Ciccotti,et al.  Mapping the network of pathways of CO diffusion in myoglobin. , 2010, Journal of the American Chemical Society.

[3]  Roderick MacKinnon,et al.  Energetic optimization of ion conduction rate by the K+ selectivity filter , 2001, Nature.

[4]  P. Dixit,et al.  Ion selectivity in the KcsA potassium channel from the perspective of the ion binding site. , 2009, Biophysical journal.

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

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

[7]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[8]  P. Bolhuis.  Transition-path sampling of β-hairpin folding , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[9]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

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

[11]  Emilio Gallicchio,et al.  Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[12]  P. A. Bash,et al.  Calculation of the relative change in binding free energy of a protein-inhibitor complex. , 1987, Science.

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

[14]  D. Beglov,et al.  Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations , 1994 .

[15]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

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

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

[18]  David M. Ferguson,et al.  On the use of acceptance ratio methods in free energy calculations , 1993 .

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

[20]  Andrew S. Paluch,et al.  Efficient Solvation Free Energy Calculations of Amino Acid Analogs by Expanded Ensemble Molecular Simulation. , 2011, Journal of chemical theory and computation.

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

[22]  G. Crooks.  Nonequilibrium Measurements of Free Energy Differences for Microscopically Reversible Markovian Systems , 1998 .

[23]  Stefan Bruckner,et al.  Unorthodox uses of Bennett's acceptance ratio method , 2009, J. Comput. Chem..

[24]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[25]  Ilsoo Kim,et al.  On the selective ion binding hypothesis for potassium channels , 2011, Proceedings of the National Academy of Sciences.

[26]  Alistair P. Rendell,et al.  The potassium channel: Structure, selectivity and diffusion , 2000 .

[27]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

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

[29]  David D L Minh,et al.  Optimized free energies from bidirectional single-molecule force spectroscopy. , 2008, Physical review letters.

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

[31]  Eric Vanden-Eijnden,et al.  Free energy and kinetics of conformational transitions from Voronoi tessellated milestoning with restraining potentials. , 2009, Journal of chemical theory and computation.

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

[33]  Arieh Warshel,et al.  Enzyme millisecond conformational dynamics do not catalyze the chemical step , 2009, Proceedings of the National Academy of Sciences.

[34]  David L Mobley,et al.  The Confine-and-Release Method: Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change. , 2007, Journal of chemical theory and computation.

[35]  Michael K. Gilson,et al.  Theory of free energy and entropy in noncovalent binding. , 2009, Chemical reviews.

[36]  B. Montgomery Pettitt,et al.  Electrostatic solvation free energy of amino acid side chain analogs: Implications for the validity of electrostatic linear response in water , 2011, J. Comput. Chem..

[37]  M. Saito,et al.  Hydration free energy calculations by the acceptance ratio method , 1990 .

[38]  C. Jarzynski.  Rare events and the convergence of exponentially averaged work values. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

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

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

[42]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[43]  Chris Oostenbrink,et al.  Free energy calculations of protein-ligand interactions. , 2011, Current opinion in chemical biology.

[44]  William L. Jorgensen,et al.  Free energy of TIP4P water and the free energies of hydration of CH4 and Cl- from statistical perturbation theory , 1989 .

[45]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[46]  Stefan Bruckner,et al.  Avoiding the van der Waals endpoint problem using serial atomic insertion , 2011, J. Comput. Chem..

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

[48]  S. Varma,et al.  Tuning ion coordination architectures to enable selective partitioning. , 2006, Biophysical journal.

[49]  Benoît Roux,et al.  Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations. , 2010, Journal of chemical theory and computation.

[50]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[51]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

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

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

[54]  M. Mezei.  Adaptive umbrella sampling: Self-consistent determination of the non-Boltzmann bias , 1987 .

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

[56]  David L Bostick,et al.  Selectivity in K+ channels is due to topological control of the permeant ion's coordinated state , 2007, Proceedings of the National Academy of Sciences.

[57]  Accurate calculations of free-energy differences by the distribution method. , 2008, The Journal of chemical physics.

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

[59]  B. Roux,et al.  Energetics of ion conduction through the K + channel , 2022 .

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

[61]  Benoît Roux,et al.  Hydration of Amino Acid Side Chains: Nonpolar and Electrostatic Contributions Calculated from Staged Molecular Dynamics Free Energy Simulations with Explicit Water Molecules , 2004 .

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

[63]  Benoit Roux,et al.  Importance of Hydration and Dynamics on the Selectivity of the KcsA and NaK Channels , 2007, The Journal of general physiology.

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

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

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

[67]  Ilsoo Kim,et al.  Mechanism of potassium-channel selectivity revealed by Na+ and Li+ binding sites within the KcsA pore , 2009, Nature Structural &Molecular Biology.

[68]  B. Roux,et al.  Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands , 2004, Nature.

[69]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[70]  Peter C. Jordan,et al.  Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study. , 2008, Biophysical journal.

[71]  Emilio Gallicchio,et al.  Advances in all atom sampling methods for modeling protein-ligand binding affinities. , 2011, Current opinion in structural biology.

[72]  C. Jarzynski,et al.  Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies , 2005, Nature.

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

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

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

[76]  A. Warshel,et al.  Toward computer-aided site-directed mutagenesis of enzymes. , 1986, Proceedings of the National Academy of Sciences of the United States of America.

[77]  Christophe Chipot,et al.  Good practices in free-energy calculations. , 2010, The journal of physical chemistry. B.

[78]  Ian R. McDonald,et al.  Machine Calculation of Thermodynamic Properties of a Simple Fluid at Supercritical Temperatures , 1967 .

[79]  G. Crooks.  Path-ensemble averages in systems driven far from equilibrium , 1999, cond-mat/9908420.

[80]  K. Shing,et al.  The chemical potential in non-ideal liquid mixtures , 1983 .

[81]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

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

[83]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

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

[85]  Vijay S Pande,et al.  Bayesian update method for adaptive weighted sampling. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[86]  J. D. de Pablo,et al.  Potential of mean force between a spherical particle suspended in a nematic liquid crystal and a substrate: sphere size effects. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[87]  David L Mobley,et al.  Small Molecule Solvation Free Energy: Enhanced Conformational Sampling Using Expanded Ensemble Molecular Dynamics Simulation. , 2011, Journal of chemical theory and computation.

[88]  Wei Yang,et al.  Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. , 2004, The Journal of chemical physics.

[89]  Overcoming barriers in macromolecular simulations: non-Boltzmann thermodynamic integration , 1998 .

[90]  Gerhard König,et al.  Non‐Boltzmann sampling and Bennett's acceptance ratio method: How to profit from bending the rules , 2011, J. Comput. Chem..

[91]  R. MacKinnon,et al.  Chemistry of ion coordination and hydration revealed by a K+ channel–Fab complex at 2.0 Å resolution , 2001, Nature.

[92]  P. Kollman,et al.  Exhaustive mutagenesis in silico: Multicoordinate free energy calculations on proteins and peptides , 2000, Proteins.

[93]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[94]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[95]  Benoît Roux,et al.  Ion selectivity of the KcsA channel: a perspective from multi-ion free energy landscapes. , 2010, Journal of molecular biology.

[96]  H. Then,et al.  Measuring the convergence of Monte Carlo free-energy calculations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[97]  N. Matubayasi,et al.  Distribution-function approach to free energy computation. , 2011, The Journal of chemical physics.

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

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

[100]  Wei Yang,et al.  Simulated scaling method for localized enhanced sampling and simultaneous "alchemical" free energy simulations: a general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations. , 2007, The Journal of chemical physics.

[101]  G. Hummer,et al.  Water conduction through the hydrophobic channel of a carbon nanotube , 2001, Nature.

[102]  Berend Smit,et al.  Molecular Dynamics Simulations , 2002 .

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

[104]  E. Vanden-Eijnden,et al.  Solvent coarse-graining and the string method applied to the hydrophobic collapse of a hydrated chain , 2007, Proceedings of the National Academy of Sciences.

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

[106]  Martin Karplus,et al.  Probability Distributions for Complex Systems: Adaptive Umbrella Sampling of the Potential Energy , 1998 .

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

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

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

[110]  Wilfred F. van Gunsteren,et al.  Comparison of three enveloping distribution sampling Hamiltonians for the estimation of multiple free energy differences from a single simulation , 2009, J. Comput. Chem..

[111]  G. Hummer,et al.  Free energy reconstruction from nonequilibrium single-molecule pulling experiments , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[113]  G. Pezzotti,et al.  Analysis of microstresses in cement paste by fluorescence piezospectroscopy. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

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

[116]  Gerhard Hummer,et al.  Error and efficiency of replica exchange molecular dynamics simulations. , 2009, The Journal of chemical physics.

[117]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[118]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

[119]  J. Hermans,et al.  Excess free energy of liquids from molecular dynamics simulations. Application to water models. , 1988, Journal of the American Chemical Society.

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

[121]  C. Jarzynski,et al.  Escorted free energy simulations. , 2011, The Journal of chemical physics.

[122]  V. Luzhkov.  On relation between the free-energy perturbation and Bennett's acceptance ratio methods: Tracing the influence of the energy gap. , 2010, The Journal of chemical physics.

[123]  Gerhard Hummer,et al.  Error and efficiency of simulated tempering simulations. , 2010, The Journal of chemical physics.

[124]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[125]  I. Shrivastava,et al.  K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study. , 2002, Biophysical journal.

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

[127]  Alexander D. MacKerell,et al.  An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids , 2000 .

[128]  Klaus Schulten,et al.  Free-energy cost for translocon-assisted insertion of membrane proteins , 2011, Proceedings of the National Academy of Sciences.

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

[130]  V. Luzhkov,et al.  Ion permeation mechanism of the potassium channel , 2000, Nature.

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

[132]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[133]  Arieh Warshel,et al.  Exploring the origin of the ion selectivity of the KcsA potassium channel , 2003, Proteins.

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