[go: up one dir, main page]
Skip to main content
Springer Nature Link
Log in

Alchemical prediction of hydration free energies for SAMPL

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

Abstract

Hydration free energy calculations have become important tests of force fields. Alchemical free energy calculations based on molecular dynamics simulations provide a rigorous way to calculate these free energies for a particular force field, given sufficient sampling. Here, we report results of alchemical hydration free energy calculations for the set of small molecules comprising the 2011 Statistical Assessment of Modeling of Proteins and Ligands challenge. Our calculations are largely based on the Generalized Amber Force Field with several different charge models, and we achieved RMS errors in the 1.4–2.2 kcal/mol range depending on charge model, marginally higher than what we typically observed in previous studies (Mobley et al. in J Phys Chem B 111(9):2242–2254, 2007, J Chem Theory Comput 5(2):350–358, 2009, J Phys Chem B 115:1329–1332, 2011; Nicholls et al. in J Med Chem 51:769–779, 2008; Klimovich and Mobley in J Comput Aided Mol Design 24(4):307–316, 2010). The test set consists of ethane, biphenyl, and a dibenzyl dioxin, as well as a series of chlorinated derivatives of each. We found that, for this set, using high-quality partial charges from MP2/cc-PVTZ SCRF RESP fits provided marginally improved agreement with experiment over using AM1-BCC partial charges as we have more typically done, in keeping with our recent findings (Mobley et al. in J Phys Chem B 115:1329–1332, 2011). Switching to OPLS Lennard–Jones parameters with AM1-BCC charges also improves agreement with experiment. We also find a number of chemical trends within each molecular series which we can explain, but there are also some surprises, including some that are captured by the calculations and some that are not.

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

Similar content being viewed by others

Notes

  1. These are perhaps not true predictions, as the data is typically available in the literature, but often in relatively obscure places, requiring individual attention from a skilled practitioner to extract hydration free energies and evaluate the accuracy of the data.

  2. the FEP branch of git, commit 9e5241b9dfc8503f887d54640183bc8e397af261

  3. A full table of all computed hydration free energies is available in the Supporting Information (Table 2).

References

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

    Article  CAS  Google Scholar 

  2. Nicholls A, Mobley DL, Guthrie JP, Chodera JD, Bayly CI, Cooper MD, Pande VS (2008) Predicting small-molecule solvation free energies: an informal blind test for computational chemistry. J Med Chem 51:769–779

    Article  CAS  Google Scholar 

  3. Mobley DL, Bayly CI, Cooper MD, Shirts MR, Dill KA (2009) 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 

  4. Klimovich P, Mobley DL (2010) Predicting hydration free energies using all-atom molecular dynamics simulations and multiple starting conformations. J Comput Aided Mol Design 24(4):307–316

    Article  CAS  Google Scholar 

  5. Mobley DL, Dumont É, Chodera JD, Dill KA (2011) Comparison of charge models for fixed-charge force fields: small molecule hydration free energies in explicit solvent. J Phys Chem B 115:1329–1332

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  7. 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 Comput Chem 24:1930–1935

    Article  CAS  Google Scholar 

  8. Shirts MR, Pitera JW, Swope WC, Pande VS (2003) Extremely precise free energy calculations of amino acid side chain analogs: comparison of common molecular mechanics force fields for proteins. J Chem Phys 119(11):5740

    Article  CAS  Google Scholar 

  9. Deng Y, Roux B (2004) 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:16567–16576

    Article  CAS  Google Scholar 

  10. 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 

  11. 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 

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

    Article  CAS  Google Scholar 

  13. Shivakumar D, Deng Y, Roux B (2009) Computations of absolute solvation free energies of small molecules using explicit and implicit solvent model. J Chem Theory Comput 5(4):919–930

    Article  CAS  Google Scholar 

  14. Shivakumar D, Williams J, Wu Y, Damm W, Shelley J, Sherman W (2010) Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J Chem Theory Comput 6:1509–1519

    Article  CAS  Google Scholar 

  15. Knight JL, Brooks CL III (2011) Surveying implicit solvent models for estimating small molecule absolute hydration free energies. J Comput Chem 32:2909–2923

    Article  CAS  Google Scholar 

  16. Deng Y, Roux B (2006) Calculation of standard binding free energies: aromatic molecules in the T4 lysozyme L99A mutant. J Chem Theory Comput 2:1255–1273

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  18. 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 

  19. Gallicchio E, Lapelosa M, Levy RM (2010) Binding energy distribution analysis method (BEDAM) for estimation of Protein−Ligand binding affinities. J Chem Theory Comput 6:2961–2977

    Article  CAS  Google Scholar 

  20. Chodera JD, Mobley DL, Shirts MR, Dixon RW, Branson K, Pande VS (2011) Alchemical free energy methods for drug discovery: progress and challenges. Curr Opin Struct Biol 21:150–160

    Article  CAS  Google Scholar 

  21. Li J, Zhu T, Hawkins GD, Winget P, Liotard DA, Truhlar DG (1999, June) Extension of the platform of applicability of SM5.42R universal solvation model

  22. 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 

  23. Geballe MT, Skillman AG, Nicholls A, Guthrie JP, Taylor PJ (2010) The SAMPL2 blind prediction challenge: introduction and overview. J Comput Aided Mol Design 24:259–279

    Article  CAS  Google Scholar 

  24. Wang J, Wolf R, Caldwell J, Kollman P, Case D (2004) Development and testing of a general amber force field. J Comput Chem 25(9):1157–1174

    Article  CAS  Google Scholar 

  25. Wang J, Wang W, Kollman P, Case DA (2006) Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model 25:247–260

    Google Scholar 

  26. Jakalian A, Bush B, Jack D, Bayly CI (2000) Fast, efficient generation of high-quality atomic charges. AM1-BCC model: I. Method. J Comput Chem 21(2):132–146

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  28. Bayly CI, Cieplak P, Cornell WD, Kollman P (1993) A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. J Phys Chem 97(40):10269–10280

    Article  CAS  Google Scholar 

  29. Dixon R (2011, March) OEAntechamber: assign and generate AMBER atom types and structural parameters. http://SimTk.org—OEAntechamber: assign and generate AMBER atom types and structural parameters

  30. Jorgensen WL, Maxwell DS, Tirado-Rives J (1996) Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J Am Chem Soc 118:11225–11236

    Article  CAS  Google Scholar 

  31. Shirts MR, Chodera JD (2008) Statistically optimal analysis of samples from multiple equilibrium states. J Chem Phys 129:124105. http://dx.doi.org/10.1063/1.2978177 See code at: http://simtk.org/home/pymbar

    Google Scholar 

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

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  34. Hess B, Bekker H, Berendsen HJC, Fraaije JGEM (1997) LINCS: A linear constraint solver for molecular simulations. J Comput Chem 18:1463–1472

    Article  CAS  Google Scholar 

  35. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA Jr, Vreven T, Kuden KN, Burant JC, Milliam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomberts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng Cy, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA (2004, September) Gaussian03, c.02 edn. Gaussian, Inc., Wallingford

  36. Cancès E, Mennucci B, Tomasi J (1997) A new integral equation formalism for the polarizable continuum model: theoretical background and applications to isotropic and anisotropic dielectrics. J Chem Phys 107(8):3032

    Article  Google Scholar 

  37. Mennucci B, Cancès E, Tomasi J (1997) Evaluation of solvent effects in isotropic and anisotropic dielectrics and in ionic solutions with a unified integral equation method: theoretical bases, computational implementation, and numerical applications. J Phys Chem B 101:10506–10517

    Article  CAS  Google Scholar 

  38. 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 

  39. Karamertzanis PG, Raiteri P, Galindo A (2010) The use of anisotropic potentials in modeling water and free energies of hydration. J Chem Theory Comput 6:1590–1607

    Article  CAS  Google Scholar 

  40. Kemp DD, Gordon MS (2008) An interpretation of the enhancement of the water dipole moment due to the presence of other water molecules. J Phys Chem A 112:4885–4894

    Article  CAS  Google Scholar 

  41. Wang J, Tingjun H (2011) Application of molecular dynamics simulations in molecular property prediction I: density and heat of vaporization. J Chem Theory Comput 7:2151–2165

    Article  CAS  Google Scholar 

  42. Kehoe C (2011, July) Personal communication. Discussion of computed hydration free energies via Skype and e-mail

  43. Kehoe CW, Fennell CJ, Dill KA (2011) Testing the semi-explicit assembly solvation model in the SAMPL3 community blind test. J Comput Aided Mol Design. doi:10.1007/s10822-011-9536-8

  44. Geballe M (2011, August) Overview of SAMPL hydration challenge. In: SAMPL 2011, Stanford

  45. Case DA, Cheatham TE, Darden T, Gohlke H, Luo R, Merz KM, Onufriev A, Simmerling C, Wang B, Woods RJ (2005) The amber biomolecular simulation programs. J Comput Chem 26:1668–1688

    Article  CAS  Google Scholar 

  46. Guthrie JP (2011, August) Overview of experimental data for hydration. In: SAMPL 2011, Stanford

  47. Cerutti DS, Baker NA, McCammon JA (2007) Solvent reaction field potential inside an uncharged globular protein: a bridge between implicit and explicit solvent models?. J Chem Phys 127(15):155101

    Article  Google Scholar 

Download references

Acknowledgments

DLM and SL acknowledge the Louisiana Board of Regents Research Competitiveness and Research Enhancement Subprograms as well as the Louisiana Optical Network Initiative (supported by the Louisiana Board of Regents Post-Katrina Support Fund Initiative grant LEQSF(2007-12)- ENH-PKSFI-PRS-01), and the National Science Foundation under NSF EPSCoR Cooperative Agreement No. EPS-1003897 with additional support from the Louisiana Board of Regents. DSC acknowledges the help of David A. Case (Rutgers) in developing the solvent reaction field QM calculations, and the support of NIH grants RR12255 and RR05969 (to Dr. Case).

Author information

Authors and Affiliations

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

    David L. Mobley & Shaui Liu

  2. Department of Chemistry and Chemical Biology, BioMaPS Institute, Rutgers, The State University of New Jersey, 610 Taylor Road, Piscataway, NJ, 08854-8066, USA

    David S. Cerutti

  3. IBM Almaden Research Center, 650 Harry Road, San Jose, CA, 95120, USA

    William C. Swope & Julia E. Rice

Authors
  1. David L. Mobley
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Shaui Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. David S. Cerutti
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. William C. Swope
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Julia E. Rice
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to David L. Mobley.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary material 1 (PDF 181 KB)

Supplementary material 2 (GZ 4089 KB)

Rights and permissions

Reprints and permissions

About this article

Cite this article

Mobley, D.L., Liu, S., Cerutti, D.S. et al. Alchemical prediction of hydration free energies for SAMPL. J Comput Aided Mol Des 26, 551–562 (2012). https://doi.org/10.1007/s10822-011-9528-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10822-011-9528-8

Keywords

Search

Navigation