Non-productive Binding Modes as a Prominent Feature of Aβ1-40 Fiber Elongation: Insights from Molecular Dynamics Simulation

Abstract

AbstractAmyloid formation has been implicated in a number of neurodegenerative diseases. The elongation of amyloid fibers is thermodynamically strongly favorable but kinetic traps exist where the incoming monomer binds in an incompatible conformation that blocks further elongation. Unfortunately, this process is difficult to follow experimentally at the atomic level. It is also too complex to simulate in full detail and thus so far has been explored either through coarse-grained simulations, which may miss many important interactions, or full atomic simulations in which the incoming peptide is constrained to be near the ideal fiber geometry. Here we use an alternate approach starting from a docked complex in which the monomer is from an experimental NMR structure of one of the major conformations in the unbound ensemble, a largely unstructured peptide with the central hydrophobic region in a 310 helix. A 1000 ns full atomic simulation in explicit solvent shows the formation of a metastable intermediate by sequential, concerted movements of both the fiber and monomer. A Markov state model shows the unfolded monomer is trapped at the end of the fiber in a set of interconverting anti-parallel {beta}-hairpin conformations. The simulation here may serve as a model for the binding of other non-{beta}-sheet conformations to amyloid fibers.