Comparing a simple theoretical model for protein folding with all-atom molecular dynamics simulations

Advances in computing have enabled microsecond all-atom molecular dynamics trajectories of protein folding that can be used to compare with and test critical assumptions of theoretical models. We show that recent simulations by the Shaw group (10, 11, 14, 15) are consistent with a key assumption of an Ising-like theoretical model that native structure grows in only a few regions of the amino acid sequence as folding progresses. The distribution of mechanisms predicted by simulating the master equation of this native-centric model for the benchmark villin subdomain, with only two adjustable thermodynamic parameters and one temperature-dependent kinetic parameter, is remarkably similar to the distribution in the molecular dynamics trajectories.