Enhancing the accuracy, the efficiency and the scope of free energy simulations

Current Opinion in Structural Biology

Volumn 15, Issue 2, 2005, Pages 164-170

  Rodinger, Tomas ^a   Pomès, Régis ^a  

^a HOSPITAL FOR SICK CHILDREN   (Canada)

Author keywords

[No Author keywords available]

Indexed keywords

ACCURACY; ALGORITHM; CALCULATION; ENERGY; ENTROPY; MATHEMATICAL COMPUTING; METHODOLOGY; PRIORITY JOURNAL; REVIEW; SIMULATION; THERMODYNAMICS;

EID: 17044382128     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.sbi.2005.03.001     Document Type: Review

References (77)

1. D.L. Beveridge, and F.M. Dicapua Free-energy via molecular simulation - applications to chemical and biomolecular systems Annu Rev Biophys Biophys Chem 18 1989 431 492
2. T.P. Straatsma, and J.A. McCammon Computational alchemy Annu Rev Phys Chem 43 1992 407 435
3. C. Chipot, and D.A. Pearlman Free energy calculations. The long and winding gilded road Mol Simul 28 2002 1 12
4. W.F. van Gunsteren, X. Daura, and A.E. Mark Computation of free energy Helv Chim Acta 85 2002 3113 3129
5. D. Frenkel, and B. Smit Understanding Molecular Simulation: From Algorithms To Applications edn 2 2002 Academic Press San Diego
6. T. Simonson, G. Archontis, and M. Karplus Free energy simulations come of age: protein-ligand recognition Acc Chem Res 35 2002 430 437
7. R.W. Zwanzig High-temperature equation of state by a perturbation method 1. Nonpolar gases J Chem Phys 22 1954 1420 1426
8. C. Jarzynski Equilibrium free-energy differences from nonequilibrium measurements: a master-equation approach Phys Rev E 56 1997 5018 5035
9. C. Jarzynski Nonequilibrium equality for free energy differences Phys Rev Lett 78 1997 2690 2693
10. G.E. Crooks Nonequilibrium measurements of free energy differences for microscopically reversible Markovian systems J Stat Phys 90 1998 1481 1487
11. G. Hummer, and A. Szabo Free energy reconstruction from nonequilibrium single-molecule pulling experiments Proc Natl Acad Sci USA 98 2001 3658 3661
12. J. Liphardt, S. Dumont, S.B. Smith, I. Tinoco, and C. Bustamante Equilibrium information from nonequilibrium measurements in an experimental test of Jarzynski's equality Science 296 2002 1832 1835
13. J. Gore, F. Ritort, and C. Bustamante Bias and error in estimates of equilibrium free-energy differences from nonequilibrium measurements Proc Natl Acad Sci USA 100 2003 12564 12569
14. G. Hummer Fast-growth thermodynamic integration: error and efficiency analysis J Chem Phys 114 2001 7330 7337
15. S. Park, and K. Schulten Calculating potentials of mean force from steered molecular dynamics simulations J Chem Phys 120 2004 5946 5961 This paper reviews and extends the theoretical foundations of non-equilibrium methods. SMD is used to irreversibly transform the system from the initial to the final state over a predefined configurational pathway. Results for the unfolding of a deca-alanine peptide by pulling with stiff springs clearly indicate that CPU resources are used most efficiently when fewer non-equilibrium trajectories are used, and when each is transformed slower and thus closer to thermodynamic equilibrium.
16. D.M. Zuckerman, and T.B. Woolf Systematic finite-sampling inaccuracy in free energy differences and other nonlinear quantities J Stat Phys 114 2004 1303 1323
17. D.M. Zuckerman, and T.B. Woolf Theory of a systematic computational error in free energy differences Phys Rev Lett 89 2002 180602
18. D.M. Zuckerman, and T.B. Woolf Overcoming finite-sampling errors in fast-switching free-energy estimates: extrapolative analysis of a molecular system Chem Phys Lett 351 2002 445 453
19. N. Lu, J. Adhikari, and D.A. Kofke Variational formula for the free energy based on incomplete sampling in a molecular simulation Phys Rev E 68 2003 026122
20. N. Lu, J.K. Singh, and D.A. Kofke Appropriate methods to combine forward and reverse free-energy perturbation averages J Chem Phys 118 2003 2977 2984 Lu et al. develop overlap sampling, a general method for eliminating the systematic sampling bias arising from exponential averaging in FEP calculations. In its simplest form, the approach is easy to apply and is formally equivalent to introducing in the perturbation a virtual intermediate state consisting of a subset of configurations common to both the reference and target states; however, it only requires two (forward and reverse) simulations between these. Bennett's acceptance ratio is shown to correspond to an optimal version of overlap sampling for averaging forward and reverse simulations, in terms of both accuracy and precision (statistical error).
21. ••]).
22. N. Lu, D. Wu, T.B. Woolf, and D.A. Kofke Using overlap and funnel sampling to obtain accurate free energies from nonequilibrium work measurements Phys Rev E 69 2004 05772 This paper extends and generalizes the overlap sampling method developed for FEP simulations to the more general context of NEW methods.
23. D. Rodriguez-Gomez, E. Darve, and A. Pohorille Assessing the efficiency of free energy calculation methods J Chem Phys 120 2004 3563 3578 The authors compare Jarzynski's identity method with the adaptive biasing force method of Darve and Pohorille [39], and find that the latter works better when the system is easily driven away from equilibrium.
24. W.L. Jorgensen, and C. Ravimohan Monte Carlo simulation of differences in free energies of hydration J Chem Phys 83 1985 3050 3054
25. C.H. Bennett Efficient estimation of free-energy differences from Monte-Carlo data J Comp Phys 22 1976 245 268
26. F.M. Ytreberg, and D.M. Zuckerman Single-ensemble nonequilibrium path-sampling estimates of free energy differences J Chem Phys 120 2004 10876 10879 In this paper, the NEW method for a transformation of interest is optimized by iteratively refining reaction pathways using transition path sampling. Applications considered here and in [27] show a very substantial improvement in sampling efficiency.
27. F.M. Ytreberg, and D.M. Zuckerman Efficient use of nonequilibrium measurement to estimate free energy differences for molecular systems J Comput Chem 25 2004 1749 1759
28. L.R. Pratt A statistical-method for identifying transition-states in high dimensional problems J Chem Phys 85 1986 5045 5048
29. Dellago C, Bolhuis PG, Geissler PL: Transition path sampling. In Advances in Chemical Physics, vol 123. John Wiley & Sons; 2002:1-78.
30. R. Radhakrishnan, and T. Schlick Biomolecular free energy profiles by a shooting/umbrella sampling protocol "BOLAS" J Chem Phys 121 2004 2436 2444 A novel approach combining transition path sampling with umbrella sampling is presented and applied to a sugar puckering transition.
31. S. Park, F. Khalili-Araghi, E. Tajkhorshid, and K. Schulten Free energy calculation from steered molecular dynamics simulations using Jarzynski's equality J Chem Phys 119 2003 3559 3566
32. B. Isralewitz, M. Gao, and K. Schulten Steered molecular dynamics and mechanical functions of proteins Curr Opin Struct Biol 11 2001 224 230
33. M.O. Jensen, S. Park, E. Tajkhorshid, and K. Schulten Energetics of glycerol conduction through aquaglyceroporin GlpF Proc Natl Acad Sci USA 99 2002 6731 6736
34. R. Amaro, E. Tajkhorshid, and Z. Luthey-Schulten Developing an energy landscape for the novel function of a (beta/alpha)(8) barrel: ammonia conduction through HisF Proc Natl Acad Sci USA 100 2003 7599 7604
35. J.G. Kirkwood Statistical mechanics of fluid mixtures J Chem Phys 3 1935 300 313
36. M.J. Mitchell, and J.A. McCammon Free-energy difference calculations by thermodynamic integration - difficulties in obtaining a precise value J Comput Chem 12 1991 271 275
37. B.L. Tembe, and J.A. McCammon Ligand receptor Interactions Comput Chem 8 1984 281 283
38. G.M. Torrie, and J.P. Valleau Non-physical sampling distributions in Monte-Carlo free-energy estimation - umbrella sampling J Comput Phys 23 1977 187 199
39. E. Darve, and A. Pohorille Calculating free energies using average force J Chem Phys 115 2001 9169 9183
40. R. Bitetti-Putzer, W. Yang, and M. Karplus Generalized ensembles serve to improve the convergence of free energy simulations Chem Phys Lett 377 2003 633 641 To improve the statistical convergence of free energy calculations, a generalized ensemble is introduced in which the integration parameter is treated as a dynamic variable. Free energy differences are obtained from TI. The availability of extra degrees of freedom, together with the artificial temperature scaling resulting from conjugate momenta of the integration variable, helps overcome kinetic traps. Reduced energy barriers lead to enhanced conformational sampling of physical space. The rate of convergence of the hydration free energy of histidine improves by an order of magnitude compared to standard TI.
41. J. Hénin, and C. Chipot Overcoming free energy barriers using unconstrained molecular dynamics simulations J Chem Phys 121 2004 2904 2914 The authors combine unconstrained simulations with adaptive biasing forces for homogeneous sampling of phase space along a spatial reaction coordinate. They show how the method can be used effectively in calculating PMFs for the unfolding of a deca-alanine peptide in vacuo, the association of guanidinium and acetate, the association of two methane molecules in water and the migration of methane across the water-vapor interface.
42. X.J. Kong, and C.L. Brooks Lambda-dynamics: a new approach to free energy calculations J Chem Phys 105 1996 2414 2423
43. Z.Y. Guo, and C.L. Brooks Rapid screening of binding affinities: application of the lambda-dynamics method to a trypsin-inhibitor system J Am Chem Soc 120 1998 1920 1921
44. R. Pomès, E. Eisenmesser, C.B. Post, and B. Roux Calculating excess chemical potentials using dynamic simulations in the fourth dimension J Chem Phys 111 1999 3387 3395
45. H. Nymeyer, S. Gnanakaran, and A.E. Garcia Atomic simulations of protein folding, using the replica exchange algorithm Methods Enzymol 383 2004 119 149
46. S. Gnanakaran, H. Nymeyer, J. Portman, K.Y. Sanbonmatsu, and A.E. Garcia Peptide folding simulations Curr Opin Struct Biol 13 2003 168 174
47. A.E. Garcia, and J.N. Onuchic Folding a protein in a computer: an atomic description of the folding/unfolding of protein A Proc Natl Acad Sci USA 100 2003 13898 13903 REMD simulations provide a thorough description of the folding-unfolding mechanism of fragment B of protein A, a three-helix bundle of 46 amino acid residues. This is the largest protein studied to date with explicit solvent. Free energy, entropy and enthalpy surfaces computed over a wide temperature range characterize the thermal dependence of the folding-unfolding equilibrium and provide a quantitative picture of the enthalpy-entropy compensation modulating the equilibrium.
48. D. Paschek, and A.E. Garcia Reversible temperature and pressure denaturation of a protein fragment: a replica exchange molecular dynamics simulation study Phys Rev Lett 93 2004 238105
49. S.V. Krivov, and M. Karplus Hidden complexity of free energy surfaces for peptide (protein) folding Proc Natl Acad Sci USA 101 2004 14766 14770
50. C.J. Woods, J.W. Essex, and M.A. King The development of replica-exchange-based free-energy methods J Phys Chem B 107 2003 13703 13710 This work combines the high sampling efficiency afforded by RE methods with finite-difference TI and FEP calculations. The method is tested and compared to more conventional approaches. It is argued that the method works better than adaptive methods because the latter are driven (i.e. non-equilibrium).
51. C.J. Woods, J.W. Essex, and M.A. King Enhanced configurational sampling in binding free-energy calculations J Phys Chem B 107 2003 13711 13718
52. W. Yang, R. Bitetti-Putzer, and M. Karplus Free energy simulations: Use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence J Chem Phys 120 2004 2618 2628 This paper presents a rigorous and straightforward procedure for removing non-equilibrium contributions from ensemble averages in equilibrium methods.
53. M.R. Shirts, J.W. Pitera, W.C. Swope, and V.S. Pande Extremely precise free energy calculations of amino acid side chain analogs: comparison of common molecular mechanics force fields for proteins J Chem Phys 119 2003 5740 5761 The authors exploit the massive computing power provided by distributed computing to calculate the hydration free energy of amino acid sidechains in several commonly used biomolecular force fields.
54. Y.Q. Deng, and B. Roux Hydration of amino acid side chains: nonpolar and electrostatic contributions calculated from staged molecular dynamics free energy simulations with explicit water molecules J Phys Chem B 108 2004 16567 16576
55. J.L. MacCallum, and D.P. Tieleman Calculation of the water-cyclohexane transfer free energies of neutral amino acid side chain analogs using the OPLS all-atom force field J Comput Chem 24 2003 1930 1935
56. H. Nymeyer, and A.E. García Simulation of the folding equilibrium of a-helical peptides: a comparison of the generalized Born approximation with explicit solvent Proc Natl Acad Sci USA 100 2003 13934 13939
57. C.M. Stultz An assessment of potential of mean force calculations with implicit solvent models J Phys Chem B 108 2004 16525 16532
58. M. Feig, and C.L. Brooks Recent advances in the development and application of implicit solvent models in biomolecule simulations Curr Opin Struct Biol 14 2004 217 224
59. C. Oostenbrink, and W.F. van Gunsteren Single-step perturbations to calculate free energy differences from unphysical reference states: limits on size, flexibility, and character J Comput Chem 24 2003 1730 1739
60. C. Oostenbrink, and W.F. van Gunsteren Free energies of binding of polychlorinated biphenyls to the estrogen receptor from a single simulation Proteins 54 2004 237 246 This paper reports the calculation of relative binding affinities of a family of polychlorinated biphenyl compounds to the estrogen receptor from a single FEP simulation with multiple ligands. The judicious choice of an unphysical reference state representative of the 17 different ligands leads to relatively small errors for all but the largest ligands.
61. J. Hermans, and S. Shankar The free-energy of xenon binding to myoglobin from molecular-dynamics simulation Isr J Chem 27 1986 225 227
62. J. Hermans, and L. Wang Inclusion of the loss of translational and rotational freedom in theoretical estimates of free energies of binding. Application to a complex of benzene and mutant T4 lysozyme J Am Chem Soc 119 1997 2707 2714
63. M.K. Gilson, J.A. Given, B.L. Bush, and J.A. McCammon The statistical-thermodynamic basis for computation of binding affinities: a critical review Biophys J 72 1997 1047 1069
64. J.A. McCammon Theory of biomolecular recognition Curr Opin Struct Biol 8 1998 245 249
65. B. Roux, M. Nina, R. Pomès, and J.C. Smith Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study Biophys J 71 1996 670 681
66. S. Boresch, F. Tettinger, M. Leitgeb, and M. Karplus Absolute binding free energies: a quantitative approach for their calculation J Phys Chem B 107 2003 9535 9551 This work proposes a general approach for imposing biases that restrict the translational and rotational motion of bound ligands in the calculation of absolute binding free energies from alchemical molecular dynamics simulations. A simple method to account for these restraints is shown to address the standard-state dependence of binding equilibria. Test cases include the binding of benzene to the L99A mutant of T4 lysozyme.
67. J.M. Swanson, R.H. Henchman, and J.A. McCammon Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculations of the association free energy Biophys J 86 2004 67 74 This paper examines the conformational free energy of both ligand and receptor molecules upon binding from a long molecular dynamics simulation. This quantity is hard to evaluate because it requires knowledge of changes in the configurational entropy of the system.
68. Z. Li, and T. Lazaridis Thermodynamic contributions of the ordered water molecule in HIV-1 protease J Am Chem Soc 125 2003 6636 6637
69. P.M. Petrone, and A.E. García MHC-peptide binding is assisted by bound water molecules J Mol Biol 338 2004 419 435
70. L.R. Olano, and S.W. Rick Hydration free energies and entropies for water in protein interiors J Am Chem Soc 126 2004 7991 8000
71. H.J. Woo, A.R. Dinner, and B. Roux Grand canonical Monte Carlo simulations of water in protein environments J Chem Phys 121 2004 6392 6400 The chemical potential of water in protein cavities in equilibrium with a bulk water reservoir is calculated in the grand canonical ensemble and illustrated for several biomolecular systems.
72. D. Hamelberg, and J.A. McCammon Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method J Am Chem Soc 126 2004 7683 7689 The authors present a rigorous method for computing the standard free energies of water molecules in the binding sites of enzyme-ligand complexes. Implications for drug design are discussed.
73. H. Schafer, L.J. Smith, A.E. Mark, and W.F. van Gunsteren Entropy calculations on the molten globule state of a protein: side-chain entropies of alpha-lactalbumin Proteins 46 2002 215 224
74. H. Gohlke, and D.A. Case Converging free energy estimates: MM-PB(GB)SA studies on the protein-protein complex Ras-Raf J Comput Chem 25 2004 238 250
75. W. Chen, C.E. Chang, and M.K. Gilson Calculation of cyclodextrin binding affinities: Energy, entropy, and implications for drug design Biophys J 87 2004 3035 3049 This paper shows that changes in host-guest conformational entropy oppose binding by as much as 24 kcal/mol in calculations of cyclodextrin binding affinities for a range of ligands. Interestingly, these changes result from a narrowing of energy wells upon binding rather than from a decrease in the number of conformers.
76. N. Lu, D.A. Kofke, and T.B. Woolf Staging is more important than perturbation method for computation of enthalpy and entropy changes in complex systems J Phys Chem B 107 2003 5598 5611 The authors examine various protocols to compute entropy and enthalpy changes from single-step and multistage FEP simulations of simple model systems. They analyze the source of errors and argue that perturbations should be performed towards the state with lower entropy.
77. C. Peter, C. Oostenbrink, A. van Dorp, and W.F. van Gunsteren Estimating entropies from molecular dynamics simulations J Chem Phys 120 2004 2652 2661 The authors examine various approaches for entropy calculations from molecular dynamics simulations. They conclude that a promising approach to obtaining entropy in alchemical transformations is to determine the temperature dependence of the free energy difference by TI.