Abstract:Amyloid peptides are known to self-assemble into larger aggregates that are linked to the pathogenesis of many neurodegenerative disorders. In contrast to primary nucleation, recent experimental and theoretical studies have shown that many toxic oligomeric species are generated through secondary processes on a pre-existing fibrillar surface. Nucleation, for example, can also occur along the surface of a pre-existing fibril-secondary nucleation-as opposed to the primary one. However, explicit pathways are still not clear. In this study, we use molecular dynamics simulation to explore the free energy landscape of a free Abeta monomer binding to an existing fibrillar surface. We specifically look into several potential Abeta structural precursors that might precede some secondary events, including elongation and secondary nucleation. We find that the overall process of surface-dependent events can be described at least by the following three stages: 1. Free diffusion 2. Downhill guiding 3. Dock and lock. And we show that the outcome of adding a new monomer onto a pre-existing fibril is pathway-dependent, which leads to different secondary processes. To understand structural details, we have identified several monomeric amyloid precursors over the fibrillar surfaces and characterize their heterogeneity using a probability contact map analysis. Using the frustration analysis (a bioinformatics tool), we show that surface heterogeneity correlates with the energy frustration of specific local residues that form binding sites on the fibrillar structure. We further investigate the helical twisting of protofilaments of different sizes and observe a length dependence on the filament twisting. This work presents a comprehensive survey over the properties of fibril growth using a combination of several openMM-based platforms, including the GPU-enabled openAWSEM package for coarse-grained modeling, MDTraj for trajectory analysis, and pyEMMA for free energy calculation. This combined approach makes long-timescale simulation for aggregation systems as well as all-in-one analysis feasible. We show that this protocol allows us to explore fibril stability, surface binding affinity/heterogeneity, as well as fibrillar twisting. All these properties are important for understanding the molecular mechanism of surface-catalyzed secondary processes of fibril growth.

Self-assembly of amyloid-forming peptides by molecular dynamics simulations.

Energy Landscapes of the Monomer and Dimer of the Alzheimer’s Peptide Aβ(1–28)

Self-Assembling ABETA(30-36) Peptides: A Combined All-Atom and Coarse-Grained Simulation Study

Exploring the Early Steps of Aggregation of Amyloid-Forming Peptide Kffe

Molecular Simulations of Amyloid Structures, Toxicity, and Inhibition

Thermodynamics and Dynamics of Amyloid Peptide Oligomerization Are Sequence Dependent

Structural, Thermodynamical, and Dynamical Properties of Oligomers Formed by the Amyloid NNQQ Peptide: Insights from Coarse-Grained Simulations

Fibril Surface-Dependent Amyloid Precursors Revealed by Coarse-Grained Molecular Dynamics Simulation

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Molecular Mechanism of the Early Stage of Amyloidogenic Hexapeptides (NFGAIL) Aggregation

Sampling the Self-Assembly Pathways of KFFE Hexamers

Peptide Self-Assembly into Amyloid Fibrils: Unbiased All-Atom Simulations

Self-Assembly Pathways of Β-Sheet-rich Amyloid-β(1-40) Dimers: Markov State Model Analysis on Millisecond Hybrid-Resolution Simulations.

Computational simulations of the early steps of protein aggregation.

Molecular Dynamics Simulations on the Oligomer-Formation Process of the GNNQQNY Peptide from Yeast Prion Protein Sup35.

Computational Simulations of the Early Steps of Protein Aggregation

Pathway Complexity of Alzheimer's β-Amyloid Aβ16-22 Peptide Assembly

Molecular Dynamics Simulations of Alzheimer Abeta40 Elongation and Lateral Association.

Oligomerization of amyloid Abeta peptides using hydrogen bonds and hydrophobicity forces

Molecular Dynamics Simulations of Alanine Rich Beta-Sheet Oligomers: Insight into Amyloid Formation.

Computer Simulation Studies of Aβ37–42 Aggregation Thermodynamics and Kinetics in Water and Salt Solution