Thermodynamic phase diagram of amyloid-β (16-22) peptide

The aggregation of monomeric 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 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. To our knowledge, this is the first in silico determination of an equilibrium (solubility) thermodynamic phase diagram for a real amyloid-forming peptide. Furthermore, the in silico prediction of Aβ16-22 solubilities over the temperature range of 277-330K agrees well with fibrillation experiments and transmission electron microscopy 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