Abstract
Quantitative free energy computation involves both using a model that is sufficiently faithful to the experimental system under study (accuracy) and establishing statistically meaningful measures of the uncertainties resulting from finite sampling (precision). We use large-scale distributed computing to access sufficient computational resources to extensively sample molecular systems and thus reduce statistical uncertainty of measured free energies. In order to examine the accuracy of a range of common models used for protein simulation, we calculate the free energy of hydration of 15 amino acid side chain analogs derived from recent versions of the OPLS-AA, CHARMM, and AMBER parameter sets in TIP3P water using thermodynamic integration. We achieve a high degree of statistical precision in our simulations, obtaining uncertainties for the free energy of hydration of 0.02–0.05 kcal/mol, which are in general an order of magnitude smaller than those found in other studies. Notably, this level of precision is comparable to that obtained in experimental hydration free energy measurements of the same molecules. Root mean square differences from experiment over the set of molecules examined using AMBER-, CHARMM-, and OPLS-AA-derived parameters were 1.35 kcal/mol, 1.31 kcal/mol, and 0.85 kcal/mol, respectively. Under the simulation conditions used, these force fields tend to uniformly underestimate solubility of all the side chain analogs. The relative free energies of hydration between amino acid side chain analogs were closer to experiment but still exhibited significant deviations. Although extensive computational resources may be needed for large numbers of molecules, sufficient computational resources to calculate precise free energy calculations for small molecules are accessible to most researchers.
References
107
Referenced
609
10.1080/08927020211974
/ Mol. Simul. (2002)10.1006/jmbi.1995.0252
/ J. Mol. Biol. (1995)10.1063/1.457251
/ J. Chem. Phys. (1989)10.1002/(SICI)1096-987X(19960715)17:9<1112::AID-JCC4>3.0.CO;2-V
/ J. Comput. Chem. (1996){'key': '2024021107520059900_r5'}
10.1021/ja00124a002
/ J. Am. Chem. Soc. (1995)10.1021/jp973084f
/ J. Phys. Chem. B (1998){'key': '2024021107520059900_r8'}
10.1021/ja9621760
/ J. Am. Chem. Soc. (1996)10.1002/jcc.10083
/ J. Comput. Chem. (2002)10.1038/nsb0902-646
/ Nat. Struct. Biol. (2002)10.1016/S0959-440X(02)00308-1
/ Curr. Opin. Struct. Biol. (2002)10.1021/jp990822m
/ J. Phys. Chem. B (1999)10.1021/jp003882x
/ J. Phys. Chem. B (2001)10.1021/jp001044x
/ J. Phys. Chem. B (2000)10.1063/1.476429
/ J. Chem. Phys. (1998)10.1021/jp990988n
/ J. Phys. Chem. B (1999)10.1021/bi00507a030
/ Biochemistry (1981)10.1063/1.447824
/ J. Chem. Phys. (1984)10.1007/BF00646936
/ J. Solution Chem. (1981)10.1021/jo00891a006
/ J. Org. Chem. (1975){'key': '2024021107520059900_r22', 'first-page': '291', 'volume': '2', 'year': '1990', 'journal-title': 'J. Chem. Soc., Perkin Trans. 2'}
/ J. Chem. Soc., Perkin Trans. 2 (1990)10.1021/ci980154m
/ J. Chem. Inf. Comput. Sci. (1999){'key': '2024021107520059900_r24'}
10.1021/ja00124a015
/ J. Am. Chem. Soc. (1995)10.1021/jp9717655
/ J. Phys. Chem. B (1998)10.1080/00268979400101331
/ Mol. Phys. (1994)10.1021/jp972440r
/ J. Phys. Chem. B (1997){'key': '2024021107520059900_r29'}
10.1016/0022-2836(76)90191-1
/ J. Mol. Biol. (1976)10.1021/bi00405a042
/ Biochemistry (1988){'key': '2024021107520059900_r32', 'first-page': '6424', 'volume': '96', 'year': '1992', 'journal-title': 'J. Phys. Chem.'}
/ J. Phys. Chem. (1992)10.1002/jcc.540141013
/ J. Comput. Chem. (1993)10.1002/jcc.10052
/ J. Comput. Chem. (2002)10.1021/j100100a043
/ J. Phys. Chem. (1994)10.1002/jcc.540140207
/ J. Comput. Chem. (1993)10.1021/ja9738539
/ J. Am. Chem. Soc. (1998)10.1021/ja001120+
/ J. Am. Chem. Soc. (2000){'key': '2024021107520059900_r39'}
10.1126/science.290.5498.1903
/ Science (2000){'key': '2024021107520059900_r41'}
10.1038/nature01160
/ Nature (London) (2002)10.1006/jmbi.2001.5033
/ J. Mol. Biol. (2001)10.1002/bip.10219
/ Biopolymers (2003)10.1016/S0022-2836(02)00888-4
/ J. Mol. Biol. (2002){'key': '2024021107520059900_r46'}
10.1110/ps.0203902
/ Protein Sci. (2002)10.1063/1.439486
/ J. Chem. Phys. (1980)10.1016/0009-2614(94)00397-1
/ Chem. Phys. Lett. (1994){'key': '2024021107520059900_r50'}
10.1021/jp003919d
/ J. Phys. Chem. B (2001)10.1063/1.476437
/ J. Chem. Phys. (1998)10.1063/1.455539
/ J. Chem. Phys. (1988)- See EPAPS Document No. E-JCPSA6-119-521329 for parameter files and a brief description. A direct link to this document may be found in the online article’s HTML reference section. The document may also be reached via the EPAPS homepage (http://www.aip.org/pubservs/epaps.html) or from ftp.aip.org in the directory /epaps/. See the EPAPS homepage for more information.
10.1063/1.445869
/ J. Chem. Phys. (1983)10.1063/1.481505
/ J. Chem. Phys. (2000)10.1063/1.442716
/ J. Chem. Phys. (1982)10.1016/0021-9991(83)90014-1
/ J. Comput. Phys. (1983)10.1080/00268978600101851
/ Mol. Phys. (1986)10.1063/1.444728
/ J. Chem. Phys. (1983){'key': '2024021107520059900_r61', 'first-page': '345', 'volume': '4', 'year': '1986', 'journal-title': 'Comput. Phys. Rep.'}
/ Comput. Phys. Rep. (1986){'key': '2024021107520059900_r62'}
{'key': '2024021107520059900_r63'}
10.1002/jcc.540150702
/ J. Comput. Chem. (1994){'key': '2024021107520059900_r65'}
{'key': '2024021107520059900_r66'}
10.1002/(SICI)1096-987X(19980730)19:10<1179::AID-JCC6>3.0.CO;2-J
/ J. Comput. Chem. (1998){'key': '2024021107520059900_r68'}
10.1063/1.461148
/ J. Chem. Phys. (1991)10.1021/jp952332+
/ J. Phys. Chem. (1996)10.1080/08927020211973
/ Mol. Simul. (2002)10.1006/jmbi.1994.1430
/ J. Mol. Biol. (1994)10.1063/1.1431588
/ J. Chem. Phys. (2002){'key': '2024021107520059900_r74'}
10.1063/1.1514572
/ J. Chem. Phys. (2002)10.1002/jcc.10117
/ J. Comput. Chem. (2002){'key': '2024021107520059900_r77'}
10.1063/1.465902
/ J. Chem. Phys. (1993)10.1063/1.466975
/ J. Chem. Phys. (1994)10.1021/ja972457n
/ J. Am. Chem. Soc. (1998)10.1021/ma000172g
/ Macromolecules (2000)10.1063/1.1328069
/ J. Chem. Phys. (2000)10.1063/1.1542598
/ J. Chem. Phys. (2003)10.1021/ma981745g
/ Macromolecules (1999)10.1063/1.1542611
/ J. Chem. Phys. (2003)10.1002/(SICI)1096-987X(199802)19:3<334::AID-JCC7>3.0.CO;2-U
/ J. Comput. Chem. (1998){'key': '2024021107520059900_r87'}
{'key': '2024021107520059900_r88'}
10.1016/S0016-7037(99)00330-0
/ Geochim. Cosmochim. Acta (2000)10.1021/ja00334a030
/ J. Am. Chem. Soc. (1984)10.1021/jp961169v
/ J. Phys. Chem. (1996)10.1063/1.1308516
/ J. Chem. Phys. (2000)10.1002/1097-0282(2001)60:2<96::AID-BIP1007>3.0.CO;2-F
/ Biopolymers (2001)10.1016/0021-9991(77)90121-8
/ J. Comput. Phys. (1977)10.1063/1.480139
/ J. Chem. Phys. (1999)10.1021/ja962310g
/ J. Am. Chem. Soc. (1997)10.1002/(SICI)1096-987X(199905)20:7<730::AID-JCC8>3.0.CO;2-T
/ J. Comput. Chem. (1999)10.1080/00268979809482364
/ Mol. Phys. (1998){'key': '2024021107520059900_r99'}
{'key': '2024021107520059900_r100'}
{'key': '2024021107520059900_r101'}
{'key': '2024021107520059900_r102'}
10.1063/1.441276
/ J. Chem. Phys. (1981)10.1063/1.441277
/ J. Chem. Phys. (1981)10.1016/0010-4655(90)90090-N
/ Comput. Phys. Commun. (1990){'key': '2024021107520059900_r106'}
{'key': '2024021107520059900_r107'}
Dates
Type | When |
---|---|
Created | 21 years, 11 months ago (Aug. 29, 2003, 10:15 a.m.) |
Deposited | 1 year, 6 months ago (Feb. 11, 2024, 2:52 a.m.) |
Indexed | 3 weeks, 3 days ago (July 30, 2025, 8:23 p.m.) |
Issued | 21 years, 11 months ago (Sept. 15, 2003) |
Published | 21 years, 11 months ago (Sept. 15, 2003) |
Published Print | 21 years, 11 months ago (Sept. 15, 2003) |
@article{Shirts_2003, title={Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins}, volume={119}, ISSN={1089-7690}, url={http://dx.doi.org/10.1063/1.1587119}, DOI={10.1063/1.1587119}, number={11}, journal={The Journal of Chemical Physics}, publisher={AIP Publishing}, author={Shirts, Michael R. and Pitera, Jed W. and Swope, William C. and Pande, Vijay S.}, year={2003}, month=sep, pages={5740–5761} }