Skip to main content

Advertisement

Springer Nature Link
Log in

Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations

  • Published:
Journal of Computer-Aided Molecular Design Aims and scope Submit manuscript

An Erratum to this article was published on 23 January 2013

Abstract

Molecular dynamics simulations in explicit solvent were applied to predict the hydration free energies for 23 small organic molecules in blind SAMPL2 test. We found good agreement with experimental results, with an RMS error of 2.82 kcal/mol over the whole set and 1.86 kcal/mol over all the molecules except several hydroxyl-rich compounds where we find evidence for a systematic error in the force field. We tested two different solvent models, TIP3P and TIP4P-Ew, and obtained very similar hydration free energies for these two models; the RMS difference was 0.64 kcal/mol. We found that preferred conformation of the carboxylic acids in water differs from that in vacuum. Surprisingly, this conformational change is not adequately sampled on simulation timescales, so we apply an umbrella sampling technique to include free energies associated with the conformational change. Overall, the results of this test reveal that the force field parameters for some groups of molecules (such as hydroxyl-rich compounds) still need to be improved, but for most compounds, accuracy was consistent with that seen in our previous tests.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7

Similar content being viewed by others

References

  1. Bennett CH (1976) Efficient estimation of free energy differences from Monte Carlo data. J Comp Phys 22:245–268

    Article  Google Scholar 

  2. Beutler TC, Mark AE, van Schaik RC, Gerber PR, van Gunsteren WF (1994) Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations. Chem Phys Lett 222:529–539

    Article  CAS  Google Scholar 

  3. Boyce SE, Mobley DL, Rocklin GJ, Graves AP, Dill KA, Shoichet BK (2009) Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site. J Mol Biol 394:747–763

    Article  CAS  Google Scholar 

  4. Budavari S (eds) (1996) The Merck index: an encyclopedia of chemicals, drugs and biologicals, 12th edn. Merck Research Laboratories

  5. Buemi G (2009) DFT study of the hydrogen bond strength and ir spectra of formic, oxalic, glyoxylic and pyruvic acids in vacuum, acetone and water solution. J Phys Org Chem 22:933–947

    Article  CAS  Google Scholar 

  6. Guthrie JP (2009) A blind challenge for computational solvation free energies: introduction and overview. J Phys Chem B 113(14):4501–4507

    Article  CAS  Google Scholar 

  7. Hess B, van der Vegt NFA (2006) Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models. J Phys Chem B 110:17616–17626

    Article  CAS  Google Scholar 

  8. Horn HW, Swope WC, Pitera JW, Madura JD, Dick TJ, Hura GL, Head-Gordon T (2004) Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J Chem Phys 120(20):9665

    Article  CAS  Google Scholar 

  9. Jakalian A, Bush BL, Jack DB, Bayly CI (2000) Fast, efficient generation of high-quality atomic charges. AM1-BCC model: I. method. J Comput Chem 21:132–146

    Article  CAS  Google Scholar 

  10. Jakalian A, Jack DB, Bayly CI (2002) Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. parameterization and validation. J Comput Chem 23:1623–1641

    Article  CAS  Google Scholar 

  11. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML (1983) Comparison of simple potential functions for simulating liquid water. J Chem Phys 79:926–935

    Article  CAS  Google Scholar 

  12. Klamt A, Eckert F, Diedenhofen M (2009) Prediction of the free energy of hydration of a challenging set of pesticide-like compounds. J Phys Chem B 113:4508–4510

    Article  CAS  Google Scholar 

  13. Klotz IM, Askounis T (1947) Absorption spectra and tautomerism of cyanuric acid, melamine and some related compounds. J Am Chem Soc 69:801–803

    Article  CAS  Google Scholar 

  14. Llinàs A, Glen RC, Goodman JM (2008) Solubility challenge: Can you predict solubilities of 32 molecules using a database of 100 reliable measurements?. J Chem Inf Model 48:1289–1303

    Article  Google Scholar 

  15. Maccallum JL, Tieleman DP (2003) Calculation of the water-cyclohexane transfer free energies of neutral amino acid side-chain analogs using the OPLS all-atom force field. J Comp Chem 24(15):1930–1935

    Article  CAS  Google Scholar 

  16. Marenich AV, Cramer CJ, Truhlar DG (2009) Performance of SM6, SM8, and SMD on the SAMPL-1 test set for the prediction of small-molecule solvation free energies. J Phys Chem B 113:4438–4543

    Google Scholar 

  17. Mobley DL, Chodera JD, Dill KA (2006) On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. J Chem Phys 125:084902

    Article  Google Scholar 

  18. Mobley DL, Chodera JD, Dill KA (2007a) Confine-and-release method: Obtaining correct binding free energies in the presence of protein conformational change. J Chem Theory Comput 3:1231–1235

    Article  CAS  Google Scholar 

  19. Mobley DL, Dumont È, Chodera JD, Dill KA (2007b) Comparison of charge models for fixed-charge force fields: small-molecule hydration free energies in explicit solvent. J Phys Chem B 111:2242–2254

    Article  CAS  Google Scholar 

  20. Mobley DL, Graves AP, Chodera JD, McReynolds A, Shoichet BK, Dill KA (2007c) Predicting absolute ligand binding free energies to a simple model site. J Mol Biol 371:1118–1134

    Article  CAS  Google Scholar 

  21. Mobley DL, Barber II AE, Fennell CJ, Dill KA (2008a) Charge asymmetries in hydration of polar solutes. J Phys Chem B 112:2405–2414

    Article  CAS  Google Scholar 

  22. Mobley DL, Chodera JD, Dill KA (2008b) Entropy and conformational change in implicit solvent simulations of small molecules. J Phys Chem B 112:938–946

    Article  CAS  Google Scholar 

  23. Mobley DL, Bayly CI, Cooper MD, Dill KA (2009a) Predictions of hydration free energies from all-atom molecular dynamics simulations. J Phys Chem B 113:4533–4537

    Article  CAS  Google Scholar 

  24. Mobley DL, Bayly CI, Cooper MD, Shirts MR, Dill KA (2009b) Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. J Chem Theory Comput 5(2):350–358

    Article  CAS  Google Scholar 

  25. Nicholls A, Mobley DL, Guthrie JP, Chodera JD, Pande VS (2008) Predicting small-molecule solvation free energies: a blind challenge test for computational chemistry. J Med Chem 51:769–778

    Article  CAS  Google Scholar 

  26. Nicholls A, Wlodek S, Grant JA (2009) The sample solvation challenge: Further lessons regarding the pitfalls of parametrization. J Phys Chem B 113:4521–4532

    Article  CAS  Google Scholar 

  27. Sancier KM, Brady AP, Lee WW (1964) Absorption spectra of solutions of cyanuric acid and its chlorinated derivatives. Spectrochimica Acta 20:397–403

    Article  CAS  Google Scholar 

  28. Shirts MR, Chodera JD (2008) Statistically optimal analysis of samples from multiple equilibrium states. J Chem Phys 129:124105

    Article  Google Scholar 

  29. Shirts MR, Pande VS (2005) Solvation free energies of amino acid side chain analogs for common molecular mechanics water models. J Chem Phys 122:134508

    Article  Google Scholar 

  30. Shirts MR, Mobley DL, Chodera JD (2007a) Alchemical free energy calculations: ready for prime time? Ann Rep Comp Chem 3:41–69

    Article  CAS  Google Scholar 

  31. Shirts MR, Mobley DL, Chodera JD, Pande VS (2007b) Accurate and efficient corrections for missing dispersion interactions in molecular simulations. J Phys Chem B 111:13052–13063

    Article  CAS  Google Scholar 

  32. Shirts MR, Mobley DL, Brown SP (2010) Free energy calculations in structure-based drug design. In: Kenneth M, Merz DR, Reynolds CH (eds) Structure based drug design, Cambridge University Press, New York (forthcoming)

  33. van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC (2005) Gromacs: fast, flexible, and free. J Comp Chem 26:1701–1718

    Article  Google Scholar 

  34. Steinbrecher T, Mobley DL, Case DA (2008) Nonlinear scaling schemes for Lennard-Jones interactions in free energy calculations. J Chem Phys 127:214108

    Article  Google Scholar 

  35. Sulea T, Wanapun D, Dennis S, Purisima EO (2009) Prediction of SAMPL-1 hydration free energies using a continuum electrostatics-dispersion model. J Phys Chem B 113:4511–4520

    Article  CAS  Google Scholar 

  36. Tan C, Tan YH, Luo R (2007) Implicit nonpolar solvent models. J Phys Chem B 111:12263–12274

    Article  CAS  Google Scholar 

  37. Tetko IV, Tanchuk VY, Villa AE (2001) Prediction of n-octanol water partition coefficients from PHYSPROP database using artificial neural networks and e-state indices. J Chem Inform Comput Sci 41(5):1407–1421

    CAS  Google Scholar 

  38. Villa A, Mark AE (2002) Calculation of the free energy of solvation for neutral analogs of amino acid side chains. J Comp Chem 23(5):548–553

    Article  CAS  Google Scholar 

  39. Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA (2004) Development and testing of a general AMBER force field. J Comput Chem 25:1157–1174

    Article  CAS  Google Scholar 

  40. Wang J, Wang W, Kollman PA, Case DA (2006) Automatic atom type and bond type perception in molecular mechanical calculations. J of Mol Graphics Modell 26:247–260

    Article  Google Scholar 

Download references

Acknowledgments

We are grateful to Huafeng Xu (D. E. Shaw Research and Development) for providing his DESMOND data, and to Christopher J. Fenell (University of California, San Francisco) for bringing the tautomerism of cyanuric acid to our attention. We are grateful to the Louisiana Optical Network Initiative Institute, supported by the Louisiana Board of Regents Post-Katrina Support Fund Initiative grant LEQSF(2007-12)-ENH-PKSFI-PRS-01 for an allocation of computer time (to DLM) and a graduate fellowship (to PVK). We thank the University of New Orleans for startup funding and for a SCoRE grant.

Author information

Authors and Affiliations

  1. Department of Chemistry, University of New Orleans, 2000 Lakeshore Drive, 70148, New Orleans, LA, USA

    Pavel V. Klimovich & David L. Mobley

Authors
  1. Pavel V. Klimovich
    View author publications

    You can also search for this author inPubMed Google Scholar

  2. David L. Mobley
    View author publications

    You can also search for this author inPubMed Google Scholar

Corresponding author

Correspondence to David L. Mobley.

Additional information

An erratum to this article can be found online at http://dx.doi.org/10.1007/s10822-012-9622-6.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Klimovich, P.V., Mobley, D.L. Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations. J Comput Aided Mol Des 24, 307–316 (2010). https://doi.org/10.1007/s10822-010-9343-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10822-010-9343-7

Keywords