2.1 and 1.8 Å Average Cα RMSD Structure Predictions on Two Small Proteins, HP-36 and S15

Abstract

On two different small proteins, the 36-mer villin headpiece domain (HP-36) and the 65-mer structured region of ribosomal protein (S15), several model predictions from the ab initio approach Rosetta were subjected to molecular dynamics simulations for refinement. After clustering the resulting trajectories into conformational families, the average molecular mechanics−Poisson Boltzmann/surface area (MM−PBSA) free energies and alpha carbon (C_α) RMSDs were then calculated for each family. Those conformational families with the lowest average free energies also contained the best C_α RMSD structures (1.4 Å for S15 and HP-36 core) and the lowest average C_α RMSDs (1.8 Å for S15, 2.1 Å for HP-36 core). For comparison, control simulations starting with the two experimental structures were very stable, each consisting of a single conformational family, with an average C_α RMSD of 1.3 Å for S15 and 1.2 Å for HP-36 core (1.9 Å over all residues). In addition, the average free energies' ranks (Spearman rank, r_s) correlate well with the average C_α RMSDs (r_s = 0.77 for HP-36, r_s = 0.83 for S15). Molecular dynamics simulations combined with the MM−PBSA free energy function provide a potentially powerful tool for the protein structure prediction community in allowing for both high-resolution structural refinement and accurate ranking of model predictions. With all of the information that genomics is now providing, this methodology may allow for advances in going from sequence to structure.