Nonfibrillar soluble oligomers, which are intermediates in the transition from monomers to amyloid fibrils, may be the toxic species in Alzheimer's disease. To monitor the early events that direct assembly of amyloidogenic peptides we probe the dynamics of formation of ( A β _16–22 ) _n by adding a monomer to a preformed ( A β _16–22 ) _{n −1} ( n = 4–6) oligomer in which the peptides are arranged in an antiparallel β-sheet conformation. All atom molecular dynamics simulations in water and multiple long trajectories, for a cumulative time of 6.9 μs, show that the oligomer grows by a two-stage dock–lock mechanism. The largest conformational change in the added disordered monomer occurs during the rapid (≈50 ns) first dock stage in which the β-strand content of the monomer increases substantially from a low initial value. In the second slow-lock phase, the monomer rearranges to form in register antiparallel structures. Surprisingly, the mobile structured oligomers undergo large conformational changes in order to accommodate the added monomer. The time needed to incorporate the monomer into the fluid-like oligomer grows even when n = 6, which suggests that the critical nucleus size must exceed six. Stable antiparallel structure formation exceeds hundreds of nanoseconds even though frequent interpeptide collisions occur at elevated monomer concentrations used in the simulations. The dock–lock mechanism should be a generic mechanism for growth of oligomers of amyloidogenic peptides.

The aim of this work is to investigate the effects of molecular mechanics force fields on amyloid peptide assembly. To this end, we performed extensive replica exchange molecular dynamics (REMD) simulations on the monomer, dimer and trimer of the seven-residue fragment of the Alzheimer's amyloid-β peptide, Aβ_{16 –22}, using the AMBER99, GROMOS96 and OPLS force fields. We compared the force fields by analysing the resulting global and local structures as well as the free energy landscapes at 300 K. We show that AMBER99 strongly favors helical structures for the monomer and does not predict any β-sheet structure for the dimer and trimer. In contrast, the dimer and trimer modeled by GROMOS96 form antiparallel β-sheet structures, while OPLS predicts diverse structures. Overall, the free energy landscapes obtained by three force fields are very different, and we also note a weak structural dependence of our results on temperature. The implications of this computational study on amyloid oligomerization, fibril growth and inhibition are also discussed.