Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine ldquodipeptidesrdquo (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution Export

@article{citeulike:1007761, abstract = {We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, beta and alpha separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of beta of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between alpha and beta at values of psi near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio beta/alpha varies between 1.5 and 2.2, the easy pass lies at psi near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around \&phis; = pi and psi = pi, whereas the amber force field gives two additional densely sampled states near \&phis; = ±100° and psi = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in high-resolution protein structures than is the case for any of the MM simulations. Proteins 2003;50:451-463. {\copyright} 2003 Wiley-Liss, Inc.}, address = {Department of Biochemistry and Biophysics, School of Medicine, University of North Carolina, Chapel Hill, North Carolina; Department of Theoretical Physics, University of Paderborn, Paderborn, Germany}, author = {Hu, Hao and Elstner, Marcus and Hermans, Jan}, citeulike-article-id = {1007761}, citeulike-linkout-0 = {http://dx.doi.org/10.1002/prot.10279}, citeulike-linkout-1 = {http://view.ncbi.nlm.nih.gov/pubmed/12557187}, citeulike-linkout-2 = {http://www.hubmed.org/display.cgi?uids=12557187}, citeulike-linkout-3 = {http://www3.interscience.wiley.com/cgi-bin/abstract/102523900/ABSTRACT}, day = {15}, doi = {10.1002/prot.10279}, issn = {1097-0134}, journal = {Proteins: Structure, Function, and Bioinformatics}, keywords = {amino, dialanine, force-fields}, month = {February}, number = {3}, pages = {451--463}, posted-at = {2009-07-26 16:12:40}, priority = {2}, title = {Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine ldquodipeptidesrdquo (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution}, url = {http://dx.doi.org/10.1002/prot.10279}, volume = {50}, year = {2003} }

TY - JOUR ID - citeulike:1007761 L3 - citeulike-article-id:1007761 N2 - We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, beta and alpha separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of beta of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between alpha and beta at values of psi near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio beta/alpha varies between 1.5 and 2.2, the easy pass lies at psi near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around &phis; = pi and psi = pi, whereas the amber force field gives two additional densely sampled states near &phis; = ±100° and psi = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in high-resolution protein structures than is the case for any of the MM simulations. Proteins 2003;50:451-463. © 2003 Wiley-Liss, Inc. IS - 3 JF - Proteins: Structure, Function, and Bioinformatics SN - 1097-0134 AD - Department of Biochemistry and Biophysics, School of Medicine, University of North Carolina, Chapel Hill, North Carolina; Department of Theoretical Physics, University of Paderborn, Paderborn, Germany EP - 463 TI - Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine ldquodipeptidesrdquo (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution VL - 50 SP - 451 KW - amino KW - dialanine KW - force-fields AU - Hu, Hao AU - Elstner, Marcus AU - Hermans, Jan PY - 2003/02/15/ UR - http://dx.doi.org/10.1002/prot.10279 ER -