Discrete molecular dynamics simulations of peptide aggregation

We study the aggregation of peptides using the discrete molecular dynamics simulations. Specifically, at temperatures above the α-helix melting temperature of a single peptide, the model peptides aggregate into a multilayer parallel β-sheet structure. This structure has an interstrand distance of 4.8 Å and an intersheet distance of 10 Å, which agree with experimental observations. Our model explains these results as follows: hydrogen-bond interactions give rise to the interstrand spacing in β sheets, while Gō interactions between side chains make β strands parallel to each other and allow β sheets to pack into layers. An important feature of our results is that the aggregates contain free edges, which may allow for further aggregation of model peptides to form elongated fibrils.