PACE, a hybrid force field which couples united-atom protein models with coarse-grained (CG) solvent, has been further optimized, aiming to improve itse ciency for folding simulations. Backbone hydration parameters have been re-optimized based on hydration free energies of polyalanyl peptides through atomistic simulations. Also, atomistic partial charges from all-atom force fields were combined with PACE in order to provide a more realistic description of interactions between charged groups. Using replica exchange molecular dynamics (REMD), ab initio folding using the new PACE has been achieved for seven small proteins (16 – 23 residues) with different structural motifs. Experimental data about folded states, such as their stability at room temperature, melting point and NMR NOE constraints, were also well reproduced. Moreover, a systematic comparison of folding kinetics at room temperature has been made with experiments, through standard MD simulations, showing that the new PACE may speed up the actual folding kinetics 5-10 times. Together with the computational speedup benefited from coarse-graining, the force field provides opportunities to study folding mechanisms. In particular, we used the new PACE to fold a 73-residue protein, 3D, in multiple 10 – 30 μs simulations, to its native states (Cα RMSD ~ 0.34 nm). Our results suggest the potential applicability of the new PACE for the study of folding and dynamics of proteins.

A critical step of β-amyloid fibril formation is fibril elongation in which amyloid-β monomers undergo structural transitions to fibrillar structures upon their binding to fibril tips. The atomic detail of the structural transitions remains poorly understood. Computational characterization of the structural transitions is limited so far to short Aβ segments (5–10 aa) owing to the long time scale of Aβ fibril elongation. To overcome the computational time scale limit, we combined a hybrid-resolution model with umbrella sampling and replica exchange molecular dynamics and performed altogether ∼1.3 ms of molecular dynamics simulations of fibril elongation for Aβ17–42. Kinetic network analysis of biased simulations resulted in a kinetic model that encompasses all Aβ segments essential for fibril formation. The model not only reproduces key properties of fibril elongation measured in experiments, including Aβ binding affinity, activation enthalpy of Aβ structural transitions and a large time scale gap (τlock/τdock = 103–104) between Aβ binding and its structural transitions, but also reveals detailed pathways involving structural transitions not seen before, namely, fibril formation both in hydrophobic regions L17-A21 and G37-A42 preceding fibril formation in hydrophilic region E22-A30. Moreover, the model identifies as important kinetic intermediates strand–loop–strand (SLS) structures of Aβ monomers, long suspected to be related to fibril elongation. The kinetic model suggests further that fibril elongation arises faster at the fibril tip with exposed L17-A21, rather than at the other tip, explaining thereby unidirectional fibril growth observed previously in experiments.

A further parametrization of a united-atom protein model coupled with coarse-grained water has been carried out to cover all amino acids (AAs). The local conformational features of each AA have been fitted on the basis of restricted coil-library statistics of high-resolution X-ray crystal structures of proteins. Potential functions were developed on the basis of combined backbone and side chain rotamer conformational preferences, or rotamer Ramachandran plots (ϕ, Ψ, χ1). Side chain-side chain and side chain-backbone interaction potentials were parametrized to fit the potential mean forces of corresponding all-atom simulations. The force field has been applied in molecular dynamics simulations of several proteins of 56-108 AA residues whose X-ray crystal and/or NMR structures are available. Starting from the crystal structures, each protein was simulated for about 100 ns. The Cα RMSDs of the calculated structures are 2.4-4.2 Å with respect to the crystal and/or NMR structures, which are still larger than but close to those of all-atom simulations (1.1-3.6 Å). Starting from the PDB structure of malate synthase G of 723 AA residues, the wall-clock time of a 30 ns simulation is about three days on a 2.65 GHz dual-core CPU. The RMSD to the experimental structure is about 4.3 Å. These results implicate the applicability of the force field in the study of protein structures.