Assembly and Kinetic Folding Pathways of a Tetrameric β-Sheet Complex: Molecular Dynamics Simulations on Simplified Off-Lattice Protein Models

We have performed discontinuous molecular dynamic simulations of the assembly and folding kinetics of a tetrameric β-sheet complex that contains four identical four-stranded antiparallel β-sheet peptides. The potential used in the simulation is a hybrid Go-type potential characterized by the bias gap parameter g, an artificial measure of a model protein's preference for its native state, and the intermolecular contact parameter η, which measures the ratio of intermolecular to intramolecular native attractions. The formation of the β-sheet complex and its equilibrium properties strongly depend on the size of the intermolecular contact parameter η. The ordered β-sheet complex in the folded state and nonaligned β-sheets or tangled chains in the misfolded state are distinguished by measuring the squared radius of gyration and the fraction of native contacts Q. The folding yield for the folded state is high at intermediate values of η, but is low at both small and large values of η. The folded state at small η is liquid-like, but is solid-like at both intermediate and large η. The misfolded state at small η contains nonaligned β-sheets and tangled chains with poor secondary structure at large η. Various folding pathways via dimeric and trimeric intermediates are observed, depending on η. Comparison with experimental results on protein aggregation indicates that intermediate η values are most appropriate for modeling fibril formation and small η values are most appropriate for modeling the formation of amorphous aggregates.