Peptide folding simulations

Current Opinion in Structural Biology

Volumn 13, Issue 2, 2003, Pages 168-174

  Gnanakaran, S ^a   Nymeyer, Hugh ^a   Portman, John ^a   Sanbonmatsu, Kevin Y ^a   García, Angel E ^a  

^a LOS ALAMOS NATIONAL LABORATORY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

PEPTIDE;

ALGORITHM; CHEMICAL MODEL; CLINICAL PROTOCOL; DNA HELIX; KINETICS; MOLECULAR DYNAMICS; NUCLEAR OVERHAUSER EFFECT; PEPTIDE SYNTHESIS; PRIORITY JOURNAL; PROBLEM SOLVING; PROTEIN CONFORMATION; PROTEIN FOLDING; REVIEW; SIMULATION; THERMODYNAMICS; VALIDATION PROCESS;

EID: 0037397987     PISSN: 0959440X     EISSN: None     Source Type: Journal    
DOI: 10.1016/S0959-440X(03)00040-X     Document Type: Review

References (65)

1. de Alba E., Antoro J., Rico M.A., Jimenez M.A. De novo design of a monomeric three-stranded anti-parallel beta-sheet. Protein Sci. 8:1999;854-865.
2. Ramirez-Alvarado M., Blanco F.J., Serrano L. Elongation of the BH8 beta-hairpin peptide: electrostatic interactions in beta-hairpin formation and stability. Protein Sci. 10:2001;1381-1392.
3. Pastor M.T., de la Paz M.L., Lacroix E., Serrano L., Perez-Paya E. Combinatorial approaches: a new tool to search for highly structured beta-hairpin peptides. Proc. Natl. Acad. Sci. U.S.A. 99:2002;614-619.
4. Neidigh J.W., Fesinmeyer R.M., Andersen N.H. Designing a 20-residue protein. Nat. Struct. Biol. 9:2002;425-430.
5. Williams S., Causgrove T.P., Gilmanshin R., Fang K.S., Callender R.H., Woodruff W.H., Dyer R.B. Fast events in protein folding: helix melting and formation in a small peptide. Biochemistry. 35:1996;691-697.
6. Munoz V., Thompson P.A., Hofrichter J.A., Eaton W.A. Folding dynamics and mechanism of beta-hairpin formation. Nature. 390:1997;196-199.
7. Crane J.C., Koepf E.K., Kelly J.W., Gruebele M. Mapping the transition state of the WW domain beta-sheet. J. Mol. Biol. 298:2000;283-292.
8. Lapidus L.J., Eaton W.A., Hofrichter J. Measuring the rate of intramolecular contact formation in polypeptides. Proc. Natl. Acad. Sci. U.S.A. 97:2000;7220-7225.
9. Sporlein S., Carstens H., Satzger H., Renner C., Behrendt R., Moroder L., Tavan P., Zinth W., Wachtveitl J. Ultrafast spectroscopy reveals subnanosecond peptide conformation dynamics and validates molecular dynamics simulation. Proc. Natl. Acad. Sci. U.S.A. 99:2002;7998-8002 A careful comparison between molecular dynamics simulations of DMSO-solvated linear and cyclic APB peptides, and femtosecond time-resolved spectroscopy attempted to validate the simulations. Absorbance spectra were compared for several 1 ns simulations in a spectral region that possessed no ground state absorption. The simulations were in good agreement with the experiments, which occurred on identical timescales, leading the authors to associate ballistic cis/trans isomerization motions with the observed spectral features.
10. Snow C.D., Nguyen N., Pande V.S., Gruebele M. Absolute comparison of simulated and experimental protein-folding dynamics. Nature. 42:2002;102-106 The authors carried out a detailed comparison of the measured and calculated folding rates for the mini-protein BBA5. The calculation consists of thousands of 10-50 ns (700 μs total simulation time) simulations using a GB/SA implicit solvent model. There is a very good agreement between measured and calculated folding rates. This study might represent the first example of a truly synergistic collaboration between theoretical and experimental work.
11. Qiu L.L., Pabit S.A., Roitberg A.E., Hagen S.J. Smaller and faster: the 20-residue Trp-cage protein folds in 4 micros. J. Am. Chem. Soc. 124:2002;12952-12953.
12. Werner J.H., Dyer R.B., Fesinmeyer R.M., Andersen N.H. Dynamics of the primary processes of protein folding: helix nucleation. J. Phys. Chem. B. 106:2002;487-494.
13. Sugita Y., Okamoto Y. Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape. Chem. Phys. Lett. 329:2000;261-270.
14. Shirts M.R., Pande V.S. Mathematical analysis of coupled parallel simulations. Phys. Rev. Lett. 86:2001;4983-4987.
15. 5] using an implicit solvent model. The forcefield is good enough to fold both β and α peptides. Folding and unfolding rates show the non-Arrhenius temperature dependence seen in simple lattice models. Folding is more cooperative for the β peptide than the α peptide, and is slower for the β peptide than the α peptide.
16. Ferrara P., Caflisch A. Folding simulations of a three-stranded antiparallel beta-sheet peptide. Proc. Natl. Acad. Sci. U.S.A. 20:2000;10780-10785.
17. Ferrara P., Caflisch A. Native topology or specific interactions: what is more important for protein folding? J. Mol. Biol. 306:2001;837-850.
18. Cavalli A., Ferrara P., Caflisch A. Weak temperature dependence of the free energy surface and folding pathways of structured peptides. Proteins. 47:2002;305-314.
19. 6. J. Phys. Chem. B. 104:2000;10080-10086.
20. Shirts M., Pande V.S. Computing: screen savers of the world unite! Science. 290:2000;1903-1904.
21. 2 peptides using Amber and Charm forcefields. Whereas the potentials of mean force from the two potentials were markedly different, the end-to-end contact times were found to be quite similar. The authors found that the measured quenching rates were consistent with the simulations [40].
22. Voter A.F. Parallel replica method for dynamics of infrequent events. Phys. Rev. B-Cond. Matt. 57:1998;13985-13988.
23. Zagrovic B., Sorin E.J., Pande V.S. β-hairpin folding simulations in atomistic detail using an implicit solvent model. J. Mol. Biol. 313:2002;151-169 In this work, the folding of a peptide has been simulated for the longest time that has ever been reported for a peptide of this size. This simulation has allowed the authors to directly observe folding events that occur on the microsecond timescale.
24. Fersht A.R. On the simulation of protein folding by short time scale molecular dynamics and distributed computing. Proc. Natl. Acad. Sci. U.S.A. 99:2002;14122-14125 In this paper, Fersht criticizes the use of multiple short simulations to study folding kinetics. He argues that protein folding is not an elementary kinetic step, because it involves a series of early conformational steps that lead to lag phases at the beginning of the kinetics. The presence of these lag phases can bias short simulations toward selecting minor pathways that have fewer or faster lag steps, and so the major folding pathways could be missed.
25. Boczko E.M., Brooks C.L. First-principles calculation of the folding free-energy of a three-helix bundle protein. Science. 269:1995;393-396.
26. Bursulaya B.D., Brooks C.L. Comparative study of the folding free energy landscape of a three-stranded beta-sheet protein with explicit and implicit solvent models. J. Phys. Chem. B. 104:2000;12378-12383 An excellent illustration of the use of umbrella sampling to sample the free energy surface of a peptide, betanova. The free energy surfaces obtained with explicit and implicit solvent models are very similar.
27. Brooks C.L. Protein and peptide folding explored with molecular simulations. Acc. Chem. Res. 35:2002;447-454.
28. Shea J.E., Brooks C.L. From folding theories to folding proteins: a review and assessment of simulation studies of protein folding and unfolding. Annu. Rev. Phys. Chem. 52:2001;499-535.
29. Berg B.A., Neuhaus T. Multicanonical algorithms for 1st order phase-transitions. Phys. Lett. B. 267:1991;249-253.
30. Hansmann U.H.E., Okamoto Y. Prediction of peptide conformation by multicanonical algorithm: new approach to the multiple-minima problem. J. Comput. Chem. 14:1993;1333-1338.
31. Alves N.A., Hansmann U.H.E. Helix formation and folding in an artificial peptide. J. Chem. Phys. 117:2002;2337-2343.
32. Kamiya N., Higo J., Nakamura H. Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water. Protein Sci. 11:2002;2297-2307.
33. Mitsutake A., Sugita Y., Okamoto Y. Generalized-ensemble algorithms for molecular simulations of biopolymers. Biopolymers. 60:2001;96-123.
34. Garcia A.E., Sanbonmatsu K.Y. Exploring the energy landscape of a beta hairpin in explicit solvent. Proteins. 42:2001;345-354 An REMD simulation of the GB1 peptide in explicit solvent shows that the peptide samples multiple free energy minima with different secondary structure.
35. Zaloj V., Elber R. Parallel computations of molecular dynamics trajectories using the stochastic path approach. Comp. Phys. Comm. 128:2000;118-127.
36. Eastman P., Gronbech-Jensen N., Doniach S. Simulation of protein folding by reaction path annealing. J. Chem. Phys. 114:2001;3823-3841.
37. Hansmann U.H.E. Generalized ensemble techniques and protein folding simulations. Comp. Phys. Comm. 147:2002;604-607.
38. Hansmann U.H.E., Wille L.T. Global optimization by energy landscape paving. Phys. Rev. Lett. 88:2002;068105.
39. Krivov S.V., Chekmarev S.F., Karplus M. Potential energy surfaces and conformational transitions in biomolecules: a successive confinement approach applied to a solvated tetrapeptide. Phys. Rev. Lett. 88:2002;8101-8104.
40. Lapidus L.J., Steinbach P.J., Eaton W.A., Szabo A., Hofrichter J. Effects of chain stiffness on the dynamics of loop formation in polypeptides. Testing a 1-dimensional diffusion model of peptide dynamics. J. Phys. Chem. B. 106:2002;11628-11640.
41. Bieri O., Wirz J., Hellrung B., Schutkowski M., Drewello M., Kiefhaber T. The speed limit for protein folding measured by triplet-triplet energy transfer. Proc. Natl. Acad. Sci. U.S.A. 96:1999;9597-9601.
42. Hudgins R.R., Huang F., Gramlich G., Nau W.M. A fluorescence-based method for direct measurement of submicrosecond intramolecular contact formation in biopolymers: an exploratory study with polypeptides. J. Am. Chem. Soc. 124:2002;556-564.
43. Lapidus L.J., Eaton W.A., Hofrichter J. Dynamics of intramolecular contact formation in polypeptides: distance dependence of quenching rates in a room-temperature glass. Phys. Rev. Lett. 87:2001;258101.
44. 3 were analyzed using bounds on the mean contact times for a general Gaussian chain. Although correlation functions of the end-to-end vector were found to be non-Gaussian, the simulated diffusion-limited contact time agrees with the Gaussian bounds if the diffusion coefficient is reduced to match the simulated end-to-end distance relaxation.
45. Garcia A.E., Sanbonmatsu K.Y. α-Helical stabilization by side chain shielding of backbone hydrogen bonds. Proc. Natl. Acad. Sci. U.S.A. 99:2002;2782-2787 REMD simulations in explicit solvent examined the mechanism of the helix-coil transition and addressed the issues of adequate sampling and forcefield accuracy. It was found that Parm94 exaggerates helix formation, whereas Parm96, designed to address this deficiency, produces predominately β-hairpin conformations for the Ala-21 system. A simple modification of the Parm94 forcefield was shown to produce σ values with the correct order of magnitude (0.004) and reasonable values for the α-helix melting temperature.
46. Vila J.A., Ripoll D.R., Scheraga H.A. Physical reasons for the unusual alpha-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides. Proc. Natl. Acad. Sci. U.S.A. 97:2002;13075-13079.
47. Smith L.J., Daura X., van Gunsteren W.F. Assessing equilibrium and convergence in biomolecular simulations. Proteins. 48:2002;487-496 The authors attempt to establish criteria for assessing the quality of sampling achieved during a simulation. The intermolecular energy, rmsd, number of clusters as determined from a cluster analysis, number of intermolecular hydrogen bonds and number of torsion angle transitions were used as measures of equilibration. These quantities were compared for various molecular dynamics simulations of α helices and β hairpins at several different temperatures. This study revealed that different measures have different convergence times. Whereas the rmsd and intermolecular interaction energy stabilize quickly (∼3 ns), the number of hydrogen bonds and number of clusters require significantly more time to equilibrate (30-50 ns).
48. Favrin G., Irback A., Wallin S. Folding of a small helical protein using hydrogen bonds and hydrophobicity forces. Proteins. 47:2002;99-105.
49. Peng Y., Hansmann U.H.E. Solvation model dependency of helix-coil transition in polyalanine. Biophys. J. 82:2002;3269-3276.
50. Jang S., Shin S., Pak Y. Molecular dynamics study of peptides in implicit water: ab initio folding of beta-hairpin, beta-sheet, and beta-beta-alpha-motif. J. Am. Chem. Soc. 124:2002;4976-4979.
51. Ma B., Nussinov R. Molecular dynamics simulations of a beta-hairpin fragment of protein G: balance between side-chain and backbone forces. J. Mol. Biol. 296:2000;1091-1104.
52. Zhou R.H., Berne B.J., Germain R. The free energy landscape for beta hairpin folding in explicit water. Proc. Natl. Acad. Sci. U.S.A. 98:2001;14931-14936.
53. Zhou R.H., Berne B.J. Can a continuum solvent model reproduce the free energy landscape of a beta-hairpin folding in water? Proc. Natl. Acad. Sci. U.S.A. 99:2002;12777-12782 REMD calculations of the free energy surface of the GB1 β hairpin in explicit and implicit (GB/SA) solvent models are in great disagreement. The authors argue that GB/SA and implicit solvent models need to improve on the energetics of ion pair formation.
54. Wu H.W., Wang S.M., Brooks B.R. Direct observation of the folding and unfolding of a beta-hairpin in explicit water through computer simulation. J. Am. Chem. Soc. 124:2002;5282-5283.
55. Roccatano D., Colombo G., Fioroni M., Mark A.E. Mechanism by which 2,2,2-trifluoroethanol/water mixtures stabilize secondary-structure formation in peptides: a molecular dynamics study. Proc. Natl. Acad. Sci. U.S.A. 99:2002;12179-12184.
56. Lee J., Shin S. Two-dimensional correlation analysis of peptide unfolding: molecular dynamics simulations of beta hairpins. J. Phys. Chem. B. 106:2002;8796-8802.
57. Plotkin S.S., Onuchic J.N. Understanding protein folding with energy landscape theory: part I: basic concepts. Q Rev. Biophys. 35:2002;111-167.
58. Still W.C., Tempczyk A., Hawley R.C., Hendrickson T. Semianalytical treatment of solvation for molecular mechanics and dynamics. J. Am. Chem. Soc. 112:1990;6127-6129.
59. Bashford D., Case D.A. Generalized Born models of macromolecular solvation effects. Annu. Rev. Phys. Chem. 51:2000;129-152.
60. Ghosh A., Rapp C.S., Friesner R.A. Generalized Born model based on a surface integral formulation. J. Phys. Chem. 102:1998;10983-10990.
61. Shen M.Y., Freed K.F. Long time dynamics of met-enkephalin: comparison of explicit and implicit solvent models. Biophys. J. 82:2002;1791-1808.
62. Daura X., Gademann K., Schafer H., Jaun B., Seebach D., van Gunsteren W.F. The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation. J. Am. Chem. Soc. 123:2001;2393-2404 A detailed conformational study on a β-hexapeptide in methanol at two different temperatures, 298K and 340K. Reversible folding to the hairpin conformation was observed at both temperatures. It was found that the ensemble of structures from the simulation at 340K reproduced the room temperature NMR-derived data, such as NOEs, better than those obtained from the simulation at 298K. It was suggested that the disagreement in the temperature could be a result of improper parameterization.
63. Peter C., Daura X., van Gunsteren W.F. Calculation of NMR-relaxation parameters for flexible molecules from molecular dynamics simulations. J. Biomol. NMR. 20:2001;297-310.
64. 6> for clustered conformations and, subsequently, spectral densities were obtained from a Lorentz function using these averages. In the third approach, the average distances are assumed to be proportional to intensities. It was found that, regardless of the methodology chosen, the correlations between experimental and theoretical intensities were identical.
65. Garcia A.E. Large-amplitude non-linear motions in proteins. Phys. Rev. Lett. 68:1992;2696-2700.