Massively parallel molecular dynamics simulations of lysozyme unfolding

We have performed molecular dynamics simulations for a total duration of more than 10 µs (with most molecular trajectories being 1 µs in duration) to study, the effect of a single mutation on hen lysozyme protein stability and denaturing, using an IBM Blue Gene/L™ supercomputer. One goal of this study was to assess the use of certain force fields to reproduce experimental results of protein unfolding using thermal denaturing techniques. A second and more important goal was to gain microscopic insights into the mechanism of protein misfolding using both thermal and chemical denaturing techniques. We found that the thermal denaturing results were robust and reproducible with various force fields. The chemical denaturing results explained why the single amino-acid mutation on residue Trp62 causes the disruption of long-range interactions in the tertiary structure. Simulation results revealed that the Trp62 residue was the key to a cooperative long-range interaction within the wild-type protein. Specifically, Trp62 acts as a bridge between two neighboring basic residues through a π-type H-bond or π-cation interaction to form an Arg-Trp-Arg "sandwich-like" structure. Our findings support the general conclusions of the experiment and provide an interesting molecular depiction of the disruption of the long-range interactions.

[1]  H. Scheraga,et al.  Contribution of unusual Arginine-Arginine short-range interactions to stabilization and recognition in proteins , 1994, Journal of protein chemistry.

[2]  D. A. Dougherty,et al.  Cation-π interactions in structural biology , 1999 .

[3]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

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

[5]  S. Buldyrev,et al.  Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model. , 2004, Biophysical journal.

[6]  V. Daggett,et al.  Long timescale simulations. , 2000, Current opinion in structural biology.

[7]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[8]  Robert S. Germain,et al.  Blue Matter, an application framework for molecular simulation on Blue Gene , 2003, J. Parallel Distributed Comput..

[9]  Christopher M. Dobson,et al.  A camelid antibody fragment inhibits the formation of amyloid fibrils by human lysozyme , 2003, Nature.

[10]  A Caflisch,et al.  Structural details of urea binding to barnase: a molecular dynamics analysis. , 1999, Structure.

[11]  Tadashi Ueda,et al.  A particular hydrophobic cluster in the residual structure of reduced lysozyme drastically affects the amyloid fibrils formation. , 2007, Biochemical and biophysical research communications.

[12]  J. Hermans,et al.  Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine “dipeptides” (Ace‐Ala‐Nme and Ace‐Gly‐Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution , 2003, Proteins.

[13]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[14]  D. A. Dougherty,et al.  Cation-pi interactions in structural biology. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[15]  B. Berne,et al.  The free energy landscape for beta hairpin folding in explicit water. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[16]  B. Kagan,et al.  Amyloidosis and Protein Folding , 2005, Science.

[17]  B. Berne,et al.  The free energy landscape for β hairpin folding in explicit water , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[19]  Ajay K. Royyuru,et al.  Blue Gene: A vision for protein science using a petaflop supercomputer , 2001, IBM Syst. J..

[20]  C. Dobson,et al.  Rationalization of the effects of mutations on peptide andprotein aggregation rates , 2003, Nature.

[21]  J. Goodfellow,et al.  Simulations of human lysozyme: probing the conformations triggering amyloidosis. , 2003, Biophysical journal.

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

[23]  C Venclovas,et al.  Some measures of comparative performance in the three CASPs , 1999, Proteins.

[24]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[25]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[26]  A. Fersht,et al.  Protein Folding and Unfolding at Atomic Resolution , 2002, Cell.

[27]  Ruhong Zhou,et al.  Thermal denaturing of mutant lysozyme with both the OPLSAA and the CHARMM force fields. , 2006, Journal of the American Chemical Society.

[28]  Lorna J. Smith,et al.  Long-Range Interactions Within a Nonnative Protein , 2002, Science.

[29]  D. Yee,et al.  Principles of protein folding — A perspective from simple exact models , 1995, Protein science : a publication of the Protein Society.

[30]  Ruhong Zhou,et al.  Destruction of long-range interactions by a single mutation in lysozyme , 2007, Proceedings of the National Academy of Sciences.

[31]  Ruhong Zhou,et al.  Observation of a dewetting transition in the collapse of the melittin tetramer , 2005, Nature.

[32]  J A McCammon,et al.  Molecular dynamics simulations of the hyperthermophilic protein sac7d from Sulfolobus acidocaldarius: contribution of salt bridges to thermostability. , 1999, Journal of molecular biology.

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

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

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

[36]  Ruhong Zhou,et al.  Hydrophobic Collapse in Multidomain Protein Folding , 2004, Science.

[37]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[38]  H Yamada,et al.  Effect of chemical modifications of tryptophan residues on the folding of reduced hen egg-white lysozyme. , 1990, Journal of biochemistry.