Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides.

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.