Abstract
Molecular dynamics (MD) simulations in conjunction with thermodynamic perturbation approach was used to calculate relative solvation free energies of five pairs of small molecules, namely; (1) methanol to ethane, (2) acetone to acetamide, (3) phenol to benzene, (4) 1,1,1 trichloroethane to ethane, and (5) phenylalanine to isoleucine. Two studies were performed to evaluate the dependence of the convergence of these calculations on MD simulation length and starting configuration. In the first study, each transformation started from the same well-equilibrated configuration and the simulation length was varied from 230 to 2,540 ps. The results indicated that for transformations involving small structural changes, a simulation length of 860 ps is sufficient to obtain satisfactory convergence. In contrast, transformations involving relatively large structural changes, such as phenylalanine to isoleucine, require a significantly longer simulation length (>2,540 ps) to obtain satisfactory convergence. In the second study, the transformation was completed starting from three different configurations and using in each case 860 ps of MD simulation. The results from this study suggest that performing one long simulation may be better than averaging results from three different simulations using a shorter simulation length and three different starting configurations.
Similar content being viewed by others
References
Zwanzig RJ (1954) J Chem Phys 22:1420–1426
Tembe BL, McCammon JA (1984) Comput Chem 8:281–286
Pearlman DA (2001) In: Reddy MR, Erion MD (eds) Free energy calculations in rational drug design. Plenum Press, New York, pp 9–35
Jorgensen WL, Ravimohan C (1985) J Chem Phys 83:3050–3054
Rao BG, Singh UC (1990) J Am Chem Soc 112:3803–3810
Agarwal A, Brown FB, Reddy MR (2001) In: Reddy MR, Erion MD (eds) Free energy calculations in rational drug design. Plenum Press, New York, pp 97–117
Rizzo RC, Tirado-Rives J, Jorgensen WL (2001) J Med Chem 44:145
Rao BG, Tilton RF, Singh UC (1992) J Am Chem Soc 114:4447
Reddy MR, Erion MD, Agarwal A (2000) In: Lipkowitz KB, Boyd DB (eds) Reviews in computational chemistry, vol 16. Wiley-VCH Inc., New York, pp 217–304
Jorgensen WL (2004) Science 303:1813–1818
Simonson T, Archontis G, Karplus M (2002) Acc Chem Res 35:430–437
Reddy MR, Erion MD (eds) (2001) Free energy calculations in rational drug design. Plenum Press, New York
Reddy MR, Erion MD (2001) J Am Chem Soc 123:6246–6252
Alvarej J, Shoichest B (eds) (2005) Virtual screening in drug discovery. CRC Press, Boca Raton, pp 1–453
Doman TN, McGovern SL, Witherbee BJ, Kasten TP, Kurumbail R, Stallings WC, Connolly DT, Shoichet BK (2002) J Med Chem 45:2213–2221
Bowen JP, Allinger NL (1991) In: Lipkowitz KB, Boyd DB (eds) Reviews in computational chemistry, vol 2. Wiley-VCH, New York, pp 81–97
Todebush PM, Bowen JP (2001) In: Reddy MR, Erion MD (eds) Free energy calculations in rational drug design. Plenum Press, New York, pp 37–59
Warshel A, Levitt M (1976) J Mol Biol 103:227–249
Aqvist J, Warshel A (1993) Chem Rev 93:2523–2544
Singh UC, Kollman PA (1986) J Comput Chem 7:718–730
Harrison MJ, Burton NA, Hillier IH (1997) J Am Chem Soc 119:12285–12291
Rosta E, Klahn M, Warshel A (2006) J Phys Chem B 110:2934–2941
Mo Y, Gao J (2006) J Phy Chem B 110:2976–2980
Devi-Kesavan LS, Gao J (2003) J Am Chem Soc 125:1532–1540
Li G, Cui Q (2003) J Phys Chem B 107:14521–14528
Riccardi D, Schaefer P, Yang Y, Yu H, Ghosh N, Prat-Resina X, Konig P, Li G, Xu D, Guo H, Elstner M, Cui Q (2006) J Phys Chem B 110:6458
Reddy MR, Singh UC, Erion MD (2004) J Am Chem Soc 126:6224–6225
Reddy MR, Singh UC, Erion MD (2007) J Comput Chem 28:491–494
Reddy MR, Erion MD (2007) J Am Chem Soc 129:7296–7297
Bash PA, Singh UC, Brown FK, Langridge R, Kollman PA (1987) Science 235:574–575
Reddy MR, Viswanadhan VN, Weinstein JN (1991) Proc Natl Acad Sci USA 88:10287–10291
Beveridge DL, DiCapua FM (1989) Annu Rev Biophys Biophys Chem 18:431–492
Rao BG, Singh UC (1989) J Am Chem Soc 111:3125–3131
Mitchell MJ, McCammon JA (1991) J Comput Chem 12:271–275
Pearlman DA, Kollman PA (1991) J Chem Phys 94:4532–4545
Guarnieri F, Still WC (1994) J Comput Chem 15:1302–1310
Pearlman DA (1994) J Comput Chem 15:105–123
Chipot C, Millot C, Maigret B, Kollman PA (1994) J Phys Chem 98:11362
Chipot C, Kollman PA, Pearlman DA (1996) J Comput Chem 17:1112–1131
Reddy MR, Erion MD (1999) J Comput Chem 20:1018–1022
Singh UC, Benkovic SJ (1988) Proc Natl Acad Sci USA 85:9519–9523
Singh UC (1988) Proc Natl Acad Sci USA 85:4280–4284
Erion MD, van Poelje PD, Reddy MR (2000) J Am Chem Soc 122:6114–6115
Erion MD, Dang Q, Reddy MR, Kasibhatla S, Jingwei Huang SJ, Lipscomb WN, van Poelje PD (2007) 133:7296–7297
Singh UC, Kollman PA (1984) J Comput Chem 5:129–145
Singh UC, Kollman PA (1986) QUEST (Version 1.1). University of California, San Francisco
Berendsen HJC, Grigera JR, Straatsma TP (1987) J Phys Chem 91:6269–6271
Reddy MR, Berkowitz M (1989) Chem Phys Lett 155:173–176
Martin ME, Aguilar MA, Chalmet S, Ruiz-Lopez MF (2002) Chem Phys 284:607–614
Verlet L (1967) Phys Rev 159:98–103
Ryckaert JP, Ciccotti G, Berendsen HJC (1997) J Comput Phys 23:327–341
Hine J, Mookerjee PK (1975) J Org Chem 40:292–298
Cabani S, Gianni P, Mollica V, Lepori L (1981) J Solut Chem 10:563–595
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the link to the electronic supplementary material.
10822_2009_9300_MOESM1_ESM.doc
Partial atomic charges and optimized geometries of all the solute molecules are submitted in the supplementary material. (DOC 40 kb)
Rights and permissions
About this article
Cite this article
Reddy, M.R., Erion, M.D. Relative solvation free energies calculated using an ab initio QM/MM-based free energy perturbation method: dependence of results on simulation length. J Comput Aided Mol Des 23, 837–843 (2009). https://doi.org/10.1007/s10822-009-9300-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10822-009-9300-5