Abstract
The aggregation of monomeric amyloid β protein (Aβ) peptide into oligomers and amyloid fibrils in the mammalian brain is associated with Alzheimer's disease. Insight into the thermodynamic stability of the Aβ peptide in different polymeric states is fundamental to defining and predicting the aggregation process. Experimental determination of Aβ thermodynamic behavior is challenging due to the transient nature of Aβ oligomers and the low peptide solubility. Furthermore, quantitative calculation of a thermodynamic phase diagram for a specific peptide requires extremely long computational times. Here, using a coarse-grained protein model, molecular dynamics (MD) simulations are performed to determine an equilibrium concentration and temperature phase diagram for the amyloidogenic peptide fragment Aβ16-22 Our results reveal that the only thermodynamically stable phases are the solution phase and the macroscopic fibrillar phase, and that there also exists a hierarchy of metastable phases. The boundary line between the solution phase and fibril phase is found by calculating the temperature-dependent solubility of a macroscopic Aβ16-22 fibril consisting of an infinite number of β-sheet layers. This in silico determination of an equilibrium (solubility) phase diagram for a real amyloid-forming peptide, Aβ16-22, over the temperature range of 277-330 K agrees well with fibrillation experiments and transmission electron microscopy (TEM) measurements of the fibril morphologies formed. This in silico approach of predicting peptide solubility is also potentially useful for optimizing biopharmaceutical production and manufacturing nanofiber scaffolds for tissue engineering.