Evaluating force field accuracy with long-time simulations of a β-hairpin tryptophan zipper peptide
Abstract
We have combined graphics processing unit-accelerated all-atom molecular dynamics with parallel tempering to explore the folding properties of small peptides in implicit solvent on the time scale of microseconds. We applied this methodology to the synthetic β-hairpin, trpzip2, and one of its sequence variants, W2W9. Each simulation consisted of over 8 μs of aggregated virtual time. Several measures of folding behavior showed good convergence, allowing comparison with experimental equilibrium properties. Our simulations suggest that the intramolecular interactions of tryptophan side chains are responsible for much of the stability of the native fold. We conclude that the ff99 force field combined with ff96 ϕ and ψ dihedral energies and an implicit solvent can reproduce plausible folding behavior in both trpzip2 and W2W9.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- January 2011
- DOI:
- 10.1063/1.3532931
- arXiv:
- arXiv:1009.0303
- Bibcode:
- 2011JChPh.134c5103H
- Keywords:
-
- molecular biophysics;
- molecular configurations;
- molecular dynamics method;
- proteins;
- 87.15.Cc;
- 87.15.B-;
- 87.15.K-;
- 87.15.hm;
- 87.15.ap;
- Folding and sequence analysis;
- Structure of biomolecules;
- Molecular interactions;
- membrane-protein interactions;
- Folding dynamics;
- Molecular dynamics simulation;
- Quantitative Biology - Biomolecules;
- Physics - Biological Physics
- E-Print:
- 10 pages, 11 figures, submitted to the Journal of Chemical Physics on June 28, 2010