β-hairpin folding simulations in atomistic detail using an implicit solvent model1

Abstract We have used distributed computing techniques and a supercluster of thousands of computer processors to study folding of the C-terminal β-hairpin from protein G in atomistic detail using the GB/SA implicit solvent model at 300 K. We have simulated a total of nearly 38 μs of folding time and obtained eight complete and independent folding trajectories. Starting from an extended state, we observe relaxation to an unfolded state characterized by non-specific, temporary hydrogen bonding. This is followed by the appearance of interactions between hydrophobic residues that stabilize a bent intermediate. Final formation of the complete hydrophobic core occurs cooperatively at the same time that the final hydrogen bonding pattern appears. The folded hairpin structures we observe all contain a closely packed hydrophobic core and proper β-sheet backbone dihedral angles, but they differ in backbone hydrogen bonding pattern. We show that this is consistent with the existing experimental data on the hairpin alone in solution. Our analysis also reveals short-lived semi-helical intermediates which define a thermodynamic trap. Our results are consistent with a three-state mechanism with a single rate-limiting step in which a varying final hydrogen bond pattern is apparent, and semi-helical off-pathway intermediates may appear early in the folding process. We include details of the ensemble dynamics methodology and a discussion of our achievements using this new computational device for studying dynamics at the atomic level.

[1]  Samuel Karlin,et al.  A First Course on Stochastic Processes , 1968 .

[2]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[3]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[4]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[5]  Michael L. McKee,et al.  Theoretical investigation of the thermal ring opening of bicyclobutane to butadiene. Evidence for a nonsynchronous process , 1988 .

[6]  T. Creighton,et al.  Protein Folding , 1992 .

[7]  A. Gronenborn,et al.  A novel, highly stable fold of the immunoglobulin binding domain of streptococcal protein G. , 1993, Science.

[8]  Axel T. Brunger,et al.  X-PLOR Version 3.1: A System for X-ray Crystallography and NMR , 1992 .

[9]  M. Sutcliffe,et al.  Determination of the solution structures of domains II and III of protein G from Streptococcus by 1H nuclear magnetic resonance. , 1992, Journal of Molecular Biology.

[10]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[11]  D. Wigley,et al.  The third IgG-binding domain from streptococcal protein G. An analysis by X-ray crystallography of the structure alone and in a complex with Fab. , 1994, Journal of molecular biology.

[12]  G L Gilliland,et al.  Two crystal structures of the B1 immunoglobulin-binding domain of streptococcal protein G and comparison with NMR. , 1994, Biochemistry.

[13]  L. Serrano,et al.  A short linear peptide that folds into a native stable β-hairpin in aqueous solution , 1994, Nature Structural Biology.

[14]  P. Argos,et al.  Knowledge‐based protein secondary structure assignment , 1995, Proteins.

[15]  L Serrano,et al.  Folding of protein G B1 domain studied by the conformational characterization of fragments comprising its secondary structure elements. , 1995, European journal of biochemistry.

[16]  Gerhard Hummer,et al.  Multi-basin dynamics of a protein in a crystal environment , 1997 .

[17]  C. Brooks,et al.  A molecular dynamics simulation study of segment B1 of protein G , 1997, Proteins.

[18]  W. C. Still,et al.  The GB/SA Continuum Model for Solvation. A Fast Analytical Method for the Calculation of Approximate Born Radii , 1997 .

[19]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[20]  V. Muñoz,et al.  Folding dynamics and mechanism of β-hairpin formation , 1997, Nature.

[21]  V. Muñoz,et al.  A statistical mechanical model for β-hairpin kinetics , 1998 .

[22]  M. Karplus,et al.  Protein Folding: A Perspective from Theory and Experiment. , 1998, Angewandte Chemie.

[23]  C L Brooks,et al.  Calculations on folding of segment B1 of streptococcal protein G. , 1998, Journal of molecular biology.

[24]  Charles L. Brooks,et al.  Molecular picture of folding of a small α/β protein , 1998 .

[25]  A. Voter.  Parallel replica method for dynamics of infrequent events , 1998 .

[26]  V S Pande,et al.  Molecular dynamics simulations of unfolding and refolding of a beta-hairpin fragment of protein G. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Andrea Amadei,et al.  A molecular dynamics study of the 41‐56 β‐hairpin from B1 domain of protein G , 1999, Protein science : a publication of the Protein Society.

[28]  C. Brooks,et al.  Folding Free Energy Surface of a Three-Stranded β-Sheet Protein , 1999 .

[29]  I. Kuntz,et al.  Inclusion of Solvation in Ligand Binding Free Energy Calculations Using the Generalized-Born Model , 1999 .

[30]  K. Hall,et al.  Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model. , 1999, Biophysical journal.

[31]  A Kolinski,et al.  Dynamics and thermodynamics of beta-hairpin assembly: insights from various simulation techniques. , 1999, Biophysical journal.

[32]  M. Karplus,et al.  Understanding beta-hairpin formation. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[33]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[34]  Alex Kentsis,et al.  D/H amide kinetic isotope effects reveal when hydrogen bonds form during protein folding , 2000, Nature Structural Biology.

[35]  R Nussinov,et al.  Molecular dynamics simulations of a beta-hairpin fragment of protein G: balance between side-chain and backbone forces. , 2000, Journal of molecular biology.

[36]  K. Hall,et al.  Experimental and computational studies of the G[UUCG]C RNA tetraloop. , 2000, Journal of molecular biology.

[37]  Vijay S. Pande,et al.  Mechanical Unfolding of a β-Hairpin Using Molecular Dynamics , 2000 .

[38]  S. Honda,et al.  Role of side-chains in the cooperative beta-hairpin folding of the short C-terminal fragment derived from streptococcal protein G. , 2000, Biochemistry.

[39]  A. Caflisch,et al.  Folding simulations of a three-stranded antiparallel β-sheet peptide , 2000 .

[40]  C. Brooks,et al.  Comparative Study of the Folding Free Energy Landscape of a Three-Stranded β-Sheet Protein with Explicit and Implicit Solvent Models , 2000 .

[41]  D. Baker,et al.  Critical role of β-hairpin formation in protein G folding , 2000, Nature Structural Biology.

[42]  Lisa J. Lapidus,et al.  Fast kinetics and mechanisms in protein folding. , 2000, Annual review of biophysics and biomolecular structure.

[43]  Ray Luo,et al.  Comparison of generalized born and poisson models: Energetics and dynamics of HIV protease , 2000 .

[44]  S. Honda,et al.  Thermodynamics of a beta-hairpin structure: evidence for cooperative formation of folding nucleus. , 2000, Journal of molecular biology.

[45]  K. Sanbonmatsu,et al.  Exploring the energy landscape of a β hairpin in explicit solvent , 2001 .

[46]  Michael R. Shirts,et al.  Mathematical analysis of coupled parallel simulations. , 2001, Physical review letters.

[47]  Sebastian Doniach,et al.  Simulation of protein folding by reaction path annealing , 2001 .