Multidimensional free‐energy calculations using the weighted histogram analysis method

The recently formulated weighted histogram analysis method (WHAM)1 is an extension of Ferrenberg and Swendsen's multiple histogram technique for free‐energy and potential of mean force calculations. As an illustration of the method, we have calculated the two‐dimensional potential of mean force surface of the dihedrals gamma and chi in deoxyadenosine with Monte Carlo simulations using the all‐atom and united‐atom representation of the AMBER force fields. This also demonstrates one of the major advantages of WHAM over umbrella sampling techniques. The method also provides an analysis of the statistical accuracy of the potential of mean force as well as a guide to the most efficient use of additional simulations to minimize errors. © 1995 John Wiley & Sons, Inc.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  A. Cooper Dynamics of Proteins and Nucleic Acids , 1988 .

[3]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[4]  Wolfram Saenger,et al.  Principles of Nucleic Acid Structure , 1983 .

[5]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[6]  Peter A. Kollman,et al.  Evaluating the assumptions underlying force field development and application using free energy conformational maps for nucleosides , 1991 .

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

[8]  J. Valleau,et al.  A Monte Carlo method for obtaining the interionic potential of mean force in ionic solution , 1975 .

[9]  J Hermans,et al.  Microfolding: Conformational probability map for the alanine dipeptide in water from molecular dynamics simulations , 1988, Proteins.

[10]  C. Brooks,et al.  Constant-temperature free energy surfaces for physical and chemical processes , 1993 .

[11]  Mihaly Mezei,et al.  Monte Carlo determination of the free energy and internal energy of hydration for the Ala dipeptide at 25.degree.C , 1985 .

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

[13]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[14]  Charles L. Brooks,et al.  Molecular dynamics with internal coordinate constraints , 1988 .

[15]  T. Ackermann C. L. Brooks III, M. Karplus, B. M. Pettitt. Proteins: A Theoretical Perspective of Dynamics, Structure and Thermodynamics, Volume LXXI, in: Advances in Chemical Physics, John Wiley & Sons, New York 1988. 259 Seiten, Preis: US $ 65.25 , 1990 .

[16]  Kumar,et al.  Efficient Monte Carlo methods for the computer simulation of biological molecules. , 1992, Physical review. A, Atomic, molecular, and optical physics.