Optimal estimates of free energies from multistate nonequilibrium work data.

We derive the optimal estimates of the free energies of an arbitrary number of thermodynamic states from nonequilibrium work measurements; the work data are collected from forward and reverse switching processes and obey a fluctuation theorem. The maximum likelihood formulation properly reweights all pathways contributing to a free energy difference and is directly applicable to simulations and experiments. We demonstrate dramatic gains in efficiency by combining the analysis with parallel tempering simulations for alchemical mutations of model amino acids.

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

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

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

[4]  Debra J. Searles,et al.  The Fluctuation Theorem , 2002 .

[5]  Christian Bartels,et al.  Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations , 1997 .

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

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

[8]  Carlos Bustamante,et al.  Direct Observation of the Three-State Folding of a Single Protein Molecule , 2005, Science.

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

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

[11]  J. Anderson.  Separate sample logistic discrimination , 1972 .

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

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

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

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

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

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

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

[19]  T. Straatsma,et al.  Multiconfiguration thermodynamic integration , 1991 .

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