Supporting Online Material Materials and Methods Figs. S1 to S10 Table S1 References Movies S1 to S3 Atomic-level Characterization of the Structural Dynamics of Proteins

Following Folding Fast Many protein functions involve conformational changes that occur on time-scales between tens of microseconds and milliseconds. This has limited the usefulness of all-atom molecular dynamics simulations, which are performed over shorter time-scales. Shaw et al. (p. 341) now report millisecond-scale, all-atom molecular dynamics simulations in an explicitly represented solvent environment. Simulation of the folding of a WW domain showed a well-defined folding pathway and simulation of the dynamics of bovine pancreatic trypsin inhibitor showed interconversion between distinct conformational states. Millisecond-scale simulations capture biologically relevant structural transitions during protein folding. Molecular dynamics (MD) simulations are widely used to study protein motions at an atomic level of detail, but they have been limited to time scales shorter than those of many biologically critical conformational changes. We examined two fundamental processes in protein dynamics—protein folding and conformational change within the folded state—by means of extremely long all-atom MD simulations conducted on a special-purpose machine. Equilibrium simulations of a WW protein domain captured multiple folding and unfolding events that consistently follow a well-defined folding pathway; separate simulations of the protein’s constituent substructures shed light on possible determinants of this pathway. A 1-millisecond simulation of the folded protein BPTI reveals a small number of structurally distinct conformational states whose reversible interconversion is slower than local relaxations within those states by a factor of more than 1000.

[1]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[2]  F. Young Biochemistry , 1955, The Indian Medical Gazette.

[3]  D. Henning Metabolism , 1972, Introduction to a Phenomenology of Life.

[4]  M. Levitt,et al.  Computer simulation of protein folding , 1975, Nature.

[5]  K. Wüthrich,et al.  NMR investigations of the dynamics of the aromatic amino acid residues in the basic pancreatic trypsin inhibitor , 1975, FEBS letters.

[6]  M. Karplus,et al.  Dynamics of folded proteins , 1977, Nature.

[7]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[8]  James T. Hynes,et al.  The stable states picture of chemical reactions. II. Rate constants for condensed and gas phase reaction models , 1980 .

[9]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[10]  Robert Huber,et al.  Structure of bovine pancreatic trypsin inhibitor , 1984 .

[11]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[12]  R. Huber,et al.  Structure of bovine pancreatic trypsin inhibitor. Results of joint neutron and X-ray refinement of crystal form II. , 1984, Journal of molecular biology.

[13]  K. Wüthrich,et al.  Reinvestigation of the aromatic side-chains in the basic pancreatic trypsin inhibitor by heteronuclear two-dimensional nuclear magnetic resonance. , 1987, Journal of molecular biology.

[14]  A. Fersht,et al.  Mapping the transition state and pathway of protein folding by protein engineering , 1989, Nature.

[15]  K Wüthrich,et al.  Protein hydration in aqueous solution. , 1991, Science.

[16]  W. Tiller The science of crystallization : macroscopic phenomena and defect generation , 1992 .

[17]  P. Wolynes,et al.  The energy landscapes and motions of proteins. , 1991, Science.

[18]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[19]  AC Tose Cell , 1993, Cell.

[20]  K. Wüthrich,et al.  Disulfide bond isomerization in BPTI and BPTI(G36S): an NMR study of correlated mobility in proteins. , 1993, Biochemistry.

[21]  M. Klein,et al.  Constant pressure molecular dynamics algorithms , 1994 .

[22]  M Karplus,et al.  Protein folding dynamics: The diffusion‐collision model and experimental data , 1994, Protein science : a publication of the Protein Society.

[23]  Jan Pieter Abrahams,et al.  Structure at 2.8 Â resolution of F1-ATPase from bovine heart mitochondria , 1994, Nature.

[24]  A. Fersht,et al.  Conformational pathway of the polypeptide chain of chymotrypsin inhibitor-2 growing from its N terminus in vitro. Parallels with the protein folding pathway. , 1995, Journal of molecular biology.

[25]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[26]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[27]  Kazuhiko Kinosita,et al.  Direct observation of the rotation of F1-ATPase , 1997, Nature.

[28]  V. Pande,et al.  On the transition coordinate for protein folding , 1998 .

[29]  A. Fersht,et al.  Acquisition of native-like interactions in C-terminal fragments of barnase. , 1999, Journal of molecular biology.

[30]  J. Ashby References and Notes , 1999 .

[31]  M. Karplus,et al.  Native proteins are surface-molten solids: application of the Lindemann criterion for the solid versus liquid state. , 1999, Journal of molecular biology.

[32]  J. Clarke,et al.  Two proteins with the same structure respond very differently to mutation: the role of plasticity in protein stability. , 2000, Journal of molecular biology.

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

[34]  J. Onuchic,et al.  Landscape approaches for determining the ensemble of folding transition states: success and failure hinge on the degree of frustration. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[35]  M. Karplus,et al.  Three key residues form a critical contact network in a protein folding transition state , 2001, Nature.

[36]  S. Pongor,et al.  A normalized root‐mean‐spuare distance for comparing protein three‐dimensional structures , 2001, Protein science : a publication of the Protein Society.

[37]  Jason C. Crane,et al.  The folding mechanism of a -sheet: the WW domain1 , 2001 .

[38]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[39]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

[40]  Valerie Daggett,et al.  Unifying features in protein-folding mechanisms , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[41]  A. Palmer,et al.  Disulfide bond isomerization in basic pancreatic trypsin inhibitor: multisite chemical exchange quantified by CPMG relaxation dispersion and chemical shift modeling. , 2003, Journal of the American Chemical Society.

[42]  Kresten Lindorff-Larsen,et al.  Calculation of mutational free energy changes in transition states for protein folding. , 2003, Biophysical journal.

[43]  K. Wüthrich,et al.  Dynamics of the aromatic amino acid residues in the globular conformation of the basic pancreatic trypsin inhibitor (BPTI) , 1976, Biophysics of structure and mechanism.

[44]  Greg L. Hura,et al.  Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. , 2004, The Journal of chemical physics.

[45]  Philip E. Dawson,et al.  Context-dependent contributions of backbone hydrogen bonding to β-sheet folding energetics , 2004, Nature.

[46]  A. Fersht,et al.  Phi-value analysis and the nature of protein-folding transition states. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[47]  G. Hummer From transition paths to transition states and rate coefficients. , 2004, The Journal of chemical physics.

[48]  J. Hofrichter,et al.  The protein folding 'speed limit'. , 2004, Current opinion in structural biology.

[49]  Michele Vendruscolo,et al.  Determination of the folding transition states of barnase by using PhiI-value-restrained simulations validated by double mutant PhiIJ-values. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[50]  Amedeo Caflisch,et al.  Phi-value analysis by molecular dynamics simulations of reversible folding. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[51]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]/Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

[52]  A. Auerbach,et al.  A speed limit for conformational change of an allosteric membrane protein. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[53]  Physikalische Gesellschaft Position-dependent diffusion coefficients and free energies from Bayesian analysis of equilibrium and replica molecular dynamics simulations , 2005 .

[54]  M. DePristo,et al.  Simultaneous determination of protein structure and dynamics , 2005, Nature.

[55]  R. Dror,et al.  Gaussian split Ewald: A fast Ewald mesh method for molecular simulation. , 2005, The Journal of chemical physics.

[56]  G. Hummer,et al.  Reaction coordinates and rates from transition paths. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[57]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[58]  J. Hofrichter,et al.  Sub-microsecond protein folding. , 2006, Journal of molecular biology.

[59]  Yan Zhang,et al.  Structure-function-folding relationship in a WW domain. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[60]  Amanda L. Jonsson,et al.  Φ-Analysis at the Experimental Limits: Mechanism of β-Hairpin Formation , 2006 .

[61]  Douglas L. Theobald,et al.  THESEUS: maximum likelihood superpositioning and analysis of macromolecular structures , 2006, Bioinform..

[62]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[63]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[64]  Lewis E. Kay,et al.  New Tools Provide New Insights in NMR Studies of Protein Dynamics , 2006, Science.

[65]  Michael P Eastwood,et al.  A common, avoidable source of error in molecular dynamics integrators. , 2007, The Journal of chemical physics.

[66]  John L. Klepeis,et al.  Anton, a special-purpose machine for molecular dynamics simulation , 2007, ISCA '07.

[67]  Francesco Luigi Gervasio,et al.  From A to B in free energy space. , 2007, The Journal of chemical physics.

[68]  Hongxing Lei,et al.  Improved sampling methods for molecular simulation. , 2007, Current opinion in structural biology.

[69]  John L. Klepeis,et al.  A scalable parallel framework for analyzing terascale molecular dynamics simulation trajectories , 2008, 2008 SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

[70]  Characterisation of transition state structures for protein folding using 'high', 'medium' and 'low' {Phi}-values. , 2008, Protein engineering, design & selection : PEDS.

[71]  Amelia A. Fuller,et al.  An experimental survey of the transition between two-state and downhill protein folding scenarios , 2008, Proceedings of the National Academy of Sciences.

[72]  R. Dror,et al.  Microsecond molecular dynamics simulation shows effect of slow loop dynamics on backbone amide order parameters of proteins. , 2008, The journal of physical chemistry. B.

[73]  Gerhard Hummer,et al.  Peptide folding kinetics from replica exchange molecular dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[74]  B. Halle,et al.  Nanosecond to microsecond protein dynamics probed by magnetic relaxation dispersion of buried water molecules. , 2008, Journal of the American Chemical Society.

[75]  D. Case,et al.  Evaluating rotational diffusion from protein MD simulations. , 2008, The journal of physical chemistry. B.

[76]  G. Hummer,et al.  Are current molecular dynamics force fields too helical? , 2008, Biophysical journal.

[77]  Peter L. Freddolino,et al.  Ten-microsecond molecular dynamics simulation of a fast-folding WW domain. , 2008, Biophysical journal.

[78]  Feng Liu,et al.  The transition state transit time of WW domain folding is controlled by energy landscape roughness. , 2009, The Journal of chemical physics.

[79]  R. Dror,et al.  Long-timescale molecular dynamics simulations of protein structure and function. , 2009, Current opinion in structural biology.

[80]  J. Kelly,et al.  Sequence determinants of thermodynamic stability in a WW domain—An all‐β‐sheet protein , 2009, Protein science : a publication of the Protein Society.

[81]  A. Cavalli,et al.  Protein conformational transitions: the closure mechanism of a kinase explored by atomistic simulations. , 2009, Journal of the American Chemical Society.

[82]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

[83]  Peter L. Freddolino,et al.  Force field bias in protein folding simulations. , 2009, Biophysical journal.

[84]  R. Rosenfeld Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[85]  John L. Klepeis,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[86]  William A Eaton,et al.  Experimental determination of upper bound for transition path times in protein folding from single-molecule photon-by-photon trajectories , 2009, Proceedings of the National Academy of Sciences.

[87]  Vijay S Pande,et al.  The Fip35 WW domain folds with structural and mechanistic heterogeneity in molecular dynamics simulations. , 2009, Biophysical journal.

[88]  Wladimir Labeikovsky,et al.  Transient Non-native Hydrogen Bonds Promote Activation of a Signaling Protein , 2009, Cell.

[89]  G. Hummer,et al.  Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides. , 2009, The journal of physical chemistry. B.

[90]  Stefano Piana,et al.  Automated Event Detection and Activity Monitoring in Long Molecular Dynamics Simulations. , 2009, Journal of chemical theory and computation.

[91]  P. Bolhuis,et al.  (Un)Folding mechanisms of the FBP28 WW domain in explicit solvent revealed by multiple rare event simulation methods. , 2010, Biophysical journal.

[92]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.