Abstract
In this study, we have investigated the sampling of π-helical conformations, with i, i + 5 hydrogen bonding, using empirical force fields. From replica exchange molecular dynamics simulations of the helical peptide acetyl-(AAQAA)3-amide using the CHARMM22 force field and implicit solvent, we find rapid conversion from an initial α-helix to π-helical conformations. While this confirms similar studies of different peptides in explicit solvent, it does not agree with experimental data where π-helices are rarely observed. The sampling of π-helices is significantly diminished in favor of canonical α-helices when a newly extended CHARMM22/CMAP force field is used, where the backbone dihedral φ/ψ potential map in a vacuum is matched exactly to high-level quantum mechanical data for an alanine dipeptide model system. Energetic analysis shows that the difference between α- and π-helical conformations of the alanine dipeptide in a vacuum is 2.6 kcal/mol according to quantum mechanical calculations while both conformations are energetically equivalent in the unmodified CHARMM22 potential. Similar trends are also found in a number of other empirical force fields, suggesting that the observation of π-helical conformations in molecular dynamics simulations is mostly a force field artifact.
Keywords
Affiliated Institutions
Related Publications
α-Helical stabilization by side chain shielding of backbone hydrogen bonds
We study atomic models of the thermodynamics of the structural transition of peptides that form α-helices. The effect of sequence variation on α-helix formation for alanine-rich...
Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone ϕ, ψ and Side-Chain χ<sub>1</sub> and χ<sub>2</sub> Dihedral Angles
While the quality of the current CHARMM22/CMAP additive force field for proteins has been demonstrated in a large number of applications, limitations in the model with respect t...
CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model
Abstract A fluctuating charge (FQ) force field is applied to molecular dynamics simulations for six small proteins in explicit polarizable solvent represented by the TIP4P‐FQ po...
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...
Improved side‐chain torsion potentials for the Amber ff99SB protein force field
Abstract Recent advances in hardware and software have enabled increasingly long molecular dynamics (MD) simulations of biomolecules, exposing certain limitations in the accurac...
Publication Info
- Year
- 2003
- Type
- article
- Volume
- 107
- Issue
- 12
- Pages
- 2831-2836
- Citations
- 205
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1021/jp027293y