Effects of all-atom force fields on amyloid oligomerization: replica exchange molecular dynamics simulations of the Aβ(16-22) dimer and trimer.

The aim of this work is to investigate the effects of molecular mechanics force fields on amyloid peptide assembly. To this end, we performed extensive replica exchange molecular dynamics (REMD) simulations on the monomer, dimer and trimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ(16-22), using the AMBER99, GROMOS96 and OPLS force fields. We compared the force fields by analysing the resulting global and local structures as well as the free energy landscapes at 300 K. We show that AMBER99 strongly favors helical structures for the monomer and does not predict any β-sheet structure for the dimer and trimer. In contrast, the dimer and trimer modeled by GROMOS96 form antiparallel β-sheet structures, while OPLS predicts diverse structures. Overall, the free energy landscapes obtained by three force fields are very different, and we also note a weak structural dependence of our results on temperature. The implications of this computational study on amyloid oligomerization, fibril growth and inhibition are also discussed.

[1]  Normand Mousseau,et al.  Coarse-grained protein molecular dynamics simulations. , 2007, The Journal of chemical physics.

[2]  Normand Mousseau,et al.  Thermodynamics and dynamics of amyloid peptide oligomerization are sequence dependent , 2009, Proteins.

[3]  John M. Hancock,et al.  K -Means Clustering. , 2010 .

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

[5]  R. Leapman,et al.  Amyloid fibril formation by A beta 16-22, a seven-residue fragment of the Alzheimer's beta-amyloid peptide, and structural characterization by solid state NMR. , 2000, Biochemistry.

[6]  Gerhard Hummer,et al.  Molecular dynamics simulations of Alzheimer's beta-amyloid protofilaments. , 2005, Journal of molecular biology.

[7]  R. Leapman,et al.  Amyloid Fibril Formation by Aβ16-22, a Seven-Residue Fragment of the Alzheimer's β-Amyloid Peptide, and Structural Characterization by Solid State NMR† , 2000 .

[8]  J Tirado-Rives,et al.  Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water. , 1991, Biochemistry.

[9]  A. Garcia,et al.  Helix‐coil transition of alanine peptides in water: Force field dependence on the folded and unfolded structures , 2005, Proteins.

[10]  Y. Sugita,et al.  Comparisons of force fields for proteins by generalized-ensemble simulations , 2004 .

[11]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[12]  Indrajit Ghosh,et al.  Discrepancies between conformational distributions of a polyalanine peptide in solution obtained from molecular dynamics force fields and amide I' band profiles. , 2010, The journal of physical chemistry. B.

[13]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[14]  Philippe Derreumaux,et al.  Effects of G33A and G33I mutations on the structures of monomer and dimer of the amyloid-β fragment 29-42 by replica exchange molecular dynamics simulations. , 2011, The journal of physical chemistry. B.

[15]  P. Lansbury,et al.  Amyloid fibrillogenesis: themes and variations. , 2000, Current opinion in structural biology.

[16]  Mai Suan Li,et al.  Relationship between population of the fibril-prone conformation in the monomeric state and oligomer formation times of peptides: insights from all-atom simulations. , 2010, The Journal of chemical physics.

[17]  Bert L de Groot,et al.  Secondary structure propensities in peptide folding simulations: a systematic comparison of molecular mechanics interaction schemes. , 2009, Biophysical journal.

[18]  J. Straub,et al.  Influence of preformed Asp23-Lys28 salt bridge on the conformational fluctuations of monomers and dimers of Abeta peptides with implications for rates of fibril formation. , 2009, The journal of physical chemistry. B.

[19]  S. Santini,et al.  Pathway Complexity of Alzheimer's β-Amyloid Aβ16-22 Peptide Assembly , 2004 .

[20]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[21]  Alessandro Laio,et al.  Stability and structure of oligomers of the Alzheimer peptide Abeta16-22: from the dimer to the 32-mer. , 2006, Biophysical journal.

[22]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[23]  Joseph A. Bank,et al.  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 , 2022 .

[24]  Yilin Yan,et al.  The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study. , 2007, Journal of molecular biology.

[25]  D. Klimov,et al.  Dissociation of Abeta(16-22) amyloid fibrils probed by molecular dynamics. , 2007, Journal of molecular biology.

[26]  C. Ross,et al.  Protein aggregation and neurodegenerative disease , 2004, Nature Medicine.

[27]  R. Leapman,et al.  Multiple quantum solid-state NMR indicates a parallel, not antiparallel, organization of β-sheets in Alzheimer's β-amyloid fibrils , 2000 .

[28]  N. Mousseau,et al.  Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation. , 2009, The journal of physical chemistry. B.

[29]  B Urbanc,et al.  Elucidation of amyloid beta-protein oligomerization mechanisms: discrete molecular dynamics study. , 2010, Journal of the American Chemical Society.

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

[31]  P. Lansbury,et al.  Atomic force microscopic imaging of seeded fibril formation and fibril branching by the Alzheimer's disease amyloid-beta protein. , 1997, Chemistry & biology.

[32]  C. Dobson,et al.  In the Footsteps of Alchemists , 2004, Science.

[33]  P. Gupta,et al.  Effect of denaturant and protein concentrations upon protein refolding and aggregation: A simple lattice model , 1998, Protein science : a publication of the Protein Society.

[34]  Gérard Vergoten,et al.  A vibrational molecular force field of model compounds with biological interest. I. Harmonic dynamics of crystalline urea at 123 K , 1990 .

[35]  Joan-Emma Shea,et al.  Effect of beta-sheet propensity on peptide aggregation. , 2009, The Journal of chemical physics.

[36]  D. Thirumalai,et al.  Dissecting the assembly of Abeta16-22 amyloid peptides into antiparallel beta sheets. , 2003, Structure.

[37]  Fabrizio Chiti,et al.  Studies of the aggregation of mutant proteins in vitro provide insights into the genetics of amyloid diseases , 2002, Proceedings of the National Academy of Sciences of the United States of America.

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

[39]  D. Otzen,et al.  Designed protein tetramer zipped together with a hydrophobic Alzheimer homology: a structural clue to amyloid assembly. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[40]  G. Favrin,et al.  Oligomerization of amyloid Abeta16-22 peptides using hydrogen bonds and hydrophobicity forces. , 2004, Biophysical journal.

[41]  Angel E. Garcia,et al.  Characterization of non-alpha helical conformations in Ala peptides , 2004 .

[42]  Amedeo Caflisch,et al.  Interpreting the aggregation kinetics of amyloid peptides. , 2006, Journal of molecular biology.

[43]  D. Thirumalai,et al.  Dissecting the assembly of A β 16-22 amyloid peptides into antiparallel β-sheets , 2002 .

[44]  R. Leapman,et al.  Multiple quantum solid-state NMR indicates a parallel, not antiparallel, organization of beta-sheets in Alzheimer's beta-amyloid fibrils. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[45]  C. Hall,et al.  Molecular dynamics simulations of spontaneous fibril formation by random-coil peptides. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[46]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

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

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

[49]  J. Shea,et al.  New insights into the mechanism of Alzheimer amyloid-beta fibrillogenesis inhibition by N-methylated peptides. , 2007, Biophysical journal.

[50]  Sarah A. Petty,et al.  Intersheet rearrangement of polypeptides during nucleation of β-sheet aggregates , 2005 .

[51]  R. Best,et al.  Balance between alpha and beta structures in ab initio protein folding. , 2010, The journal of physical chemistry. B.

[52]  P. Nguyen,et al.  Energy landscape of a small peptide revealed by dihedral angle principal component analysis , 2004, Proteins.

[53]  P. Derreumaux,et al.  Targeting the early steps of Aβ16–22 protofibril disassembly by N‐methylated inhibitors: A numerical study , 2009, Proteins.

[54]  D. van der Spoel,et al.  A temperature predictor for parallel tempering simulations. , 2008, Physical chemistry chemical physics : PCCP.

[55]  C. Dobson,et al.  Protein misfolding, functional amyloid, and human disease. , 2006, Annual review of biochemistry.

[56]  M. Hecht,et al.  De novo amyloid proteins from designed combinatorial libraries. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[57]  F. Cohen,et al.  Conformational propagation with prion‐like characteristics in a simple model of protein folding , 2001, Protein science : a publication of the Protein Society.

[58]  Martin T Zanni,et al.  Automated 2D IR spectroscopy using a mid-IR pulse shaper and application of this technology to the human islet amyloid polypeptide , 2007, Proceedings of the National Academy of Sciences.

[59]  Y. Kallberg,et al.  Prediction of Amyloid Fibril-forming Proteins* , 2001, The Journal of Biological Chemistry.

[60]  Christos Boutsidis,et al.  Atomic-level characterization of the ensemble of the Aβ(1-42) monomer in water using unbiased molecular dynamics simulations and spectral algorithms. , 2011, Journal of molecular biology.

[61]  A. Caflisch,et al.  The role of side-chain interactions in the early steps of aggregation: Molecular dynamics simulations of an amyloid-forming peptide from the yeast prion Sup35 , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[62]  Ruth Nussinov,et al.  Simulations as analytical tools to understand protein aggregation and predict amyloid conformation. , 2006, Current opinion in chemical biology.

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

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

[65]  C. Blake,et al.  The structure of amyloid fibrils by electron microscopy and X-ray diffraction. , 1997, Advances in protein chemistry.

[66]  D. Selkoe.  Folding proteins in fatal ways , 2003, Nature.

[67]  S. Santini,et al.  In Silico Assembly of Alzheimer's Aβ16-22 Peptide into β-Sheets , 2004 .

[68]  Haruki Nakamura,et al.  Peptide free‐energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95 , 2000 .

[69]  R. Riek,et al.  3D structure of Alzheimer's amyloid-β(1–42) fibrils , 2005 .

[70]  Alanine polypeptide structural fingerprints at room temperature: what can be gained from non-harmonic Car-Parrinello molecular dynamics simulations. , 2008, The journal of physical chemistry. A.

[71]  D Thirumalai,et al.  Factors governing fibrillogenesis of polypeptide chains revealed by lattice models. , 2010, Physical review letters.

[72]  Nicolas L. Fawzi,et al.  Protofibril assemblies of the arctic, Dutch, and Flemish mutants of the Alzheimer's Abeta1-40 peptide. , 2008, Biophysical journal.

[73]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[74]  Ehud Gazit,et al.  A possible role for π‐stacking in the self‐assembly of amyloid fibrils , 2002, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[75]  L. Serpell,et al.  The protofilament substructure of amyloid fibrils. , 2000, Journal of molecular biology.

[76]  Xin Jin,et al.  K-Means Clustering , 2010, Encyclopedia of Machine Learning.

[77]  Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment. , 2008, The journal of physical chemistry. B.

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

[79]  M. Emmerling,et al.  Morphology and Toxicity of Aβ-(1-42) Dimer Derived from Neuritic and Vascular Amyloid Deposits of Alzheimer's Disease* , 1996, The Journal of Biological Chemistry.

[80]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[81]  Robert A. Grothe,et al.  Structure of the cross-beta spine of amyloid-like fibrils. , 2005, Nature.

[82]  J R Ghilardi,et al.  1H NMR of A beta amyloid peptide congeners in water solution. Conformational changes correlate with plaque competence. , 1995, Biochemistry.

[83]  P. Derreumaux,et al.  Impact of the mutation A21G (Flemish variant) on Alzheimer's beta-amyloid dimers by molecular dynamics simulations. , 2006, Biophysical journal.

[84]  N. Mousseau,et al.  Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent. , 2009, The journal of physical chemistry. B.

[85]  Gerhard Stock,et al.  Conformational dynamics of trialanine in water. 2. Comparison of AMBER, CHARMM, GROMOS, and OPLS force fields to NMR and infrared experiments , 2003 .

[86]  D Thirumalai,et al.  Monomer adds to preformed structured oligomers of Aβ-peptides by a two-stage dock–lock mechanism , 2007, Proceedings of the National Academy of Sciences.

[87]  Warner L. Peticolas,et al.  Conformational studies of neuroactive ligands. 1. Force field and vibrational spectra of crystalline acetylcholine , 1989 .

[88]  D. Thirumalai,et al.  Exploring protein aggregation and self‐propagation using lattice models: Phase diagram and kinetics , 2002, Protein science : a publication of the Protein Society.

[89]  S. Lipton,et al.  Molecular pathways to neurodegeneration , 2004, Nature Medicine.

[90]  D Thirumalai,et al.  Probing the mechanisms of fibril formation using lattice models. , 2008, The Journal of chemical physics.

[91]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[92]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[93]  J. A. Hartigan,et al.  A k-means clustering algorithm , 1979 .

[94]  Ruth Nussinov,et al.  Atomic-Level Description of Amyloid β-Dimer Formation , 2006 .

[95]  Gerhard Hummer,et al.  Molecular dynamics simulations of Alzheimer's β-amyloid protofilaments , 2005 .

[96]  N. Mousseau,et al.  Energy landscapes of the monomer and dimer of the Alzheimer's peptide Abeta(1-28). , 2008, The Journal of chemical physics.

[97]  Brigita Urbanc,et al.  In silico study of amyloid β-protein folding and oligomerization , 2004 .

[98]  S. Decatur,et al.  Intersheet rearrangement of polypeptides during nucleation of {beta}-sheet aggregates. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[99]  Joan-Emma Shea,et al.  Diversity of kinetic pathways in amyloid fibril formation. , 2009, The Journal of chemical physics.

[100]  Y. Okamoto,et al.  Folding simulations of three proteins having all α-helix, all β-strand and α/β-structures , 2010 .

[101]  Fumio Hirata,et al.  The effects of solvent on the conformation and the collective motions of protein: normal mode analysis and molecular dynamics simulations of melittin in water and in vacuum , 1991 .

[102]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[103]  S. Gnanakaran,et al.  Validation of an all-atom protein force field: From dipeptides to larger peptides , 2003 .

[104]  R. Leapman,et al.  A structural model for Alzheimer's β-amyloid fibrils based on experimental constraints from solid state NMR , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[105]  Jean-Philip Piquemal,et al.  Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. , 2010, Journal of Chemical Theory and Computation.

[106]  N. Mousseau,et al.  Role of the Region 23-28 in Aβ Fibril Formation: Insights from Simulations of the Monomers and Dimers of Alzheimers Peptides Aβ40 and Aβ42 , 2008 .