Construction of the Free Energy Landscape of Peptide Aggregation from Molecular Dynamics Simulations
journal contributionposted on 10.04.2012 by Laura Riccardi, Phuong H. Nguyen, Gerhard Stock
Any type of content formally published in an academic journal, usually following a peer-review process.
To describe the structure and dynamics of oligomers during peptide aggregation, a method is proposed that considers both the intramolecular and intermolecular structures of the multimolecule system and correctly accounts for its degeneracy. The approach is based on the “by-parts” strategy, which partitions a complex molecular system into parts, determines the metastable conformational states of each part, and describes the overall conformational state of the system in terms of a product basis of the states of the parts. Starting from a molecular dynamics simulation of n molecules, the method consists of three steps: (i) characterization of the intramolecular structure, that is, of the conformational states of a single molecule in the presence of the other molecules (e.g., β-strand or random coil); (ii) characterization of the intermolecular structure through the identification of all occurring aggregate states of the peptides (dimers, trimers, etc.); and (iii) construction of the overall conformational states of the system in terms of a product basis of the n “single-molecule” states and the aggregate states. Considering the Alzheimer β-amyloid peptide fragment Aβ16–22 as a first application, about 700 overall conformational states of the trimer (Aβ16–22)3 were constructed from all-atom molecular dynamics simulation in explicit water. Based on these states, a transition network reflecting the free energy landscape of the aggregation process can be constructed that facilitates the identification of the aggregation pathways.