Abstract
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, β and α 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 β 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 α and β at values of ψ 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 β/α varies between 1.5 and 2.2, the easy pass lies at ψ 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 ϕ = π and ψ = π, whereas the amber force field gives two additional densely sampled states near ϕ = ±100° and ψ = 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.
Keywords
MeSH Terms
Affiliated Institutions
Related Publications
A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations
Abstract Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third‐generation point...
Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations
Abstract Computational studies of proteins based on empirical force fields represent a powerful tool to obtain structure–function relationships at an atomic level, and are centr...
Energy landscape of a peptide consisting of α‐helix, 3<sub>10</sub>‐helix, β‐turn, β‐hairpin, and other disordered conformations
Abstract The energy landscape of a peptide [Ace‐Lys‐Gln‐Cys‐Arg‐Glu‐Arg‐Ala‐Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER...
Peptide free-energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95
The C96 and AMBER95 force fields were compared with small model peptides Ac-(Ala)n-NMe (Ac = CH3CO, NMe = NHCH3, n=2 and 3) in vacuo and in TIP3P water by computing the free-ene...
ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB
Molecular mechanics is powerful for its speed in atomistic simulations, but an accurate force field is required. The Amber ff99SB force field improved protein secondary structur...
Publication Info
- Year
- 2003
- Type
- article
- Volume
- 50
- Issue
- 3
- Pages
- 451-463
- Citations
- 263
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1002/prot.10279
- PMID
- 12557187