Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: description of the folding pathway.

Residual structure in the denatured state of a protein may contain clues about the early events in folding. We have simulated by molecular dynamics the denatured state of barnase, which has been studied by NMR spectroscopy. An ensemble of 10(4) structures was generated after 2 ns of unfolding and following for a further 2 ns. The ensemble was heterogeneous, but there was nonrandom, residual structure with persistent interactions. Helical structure in the C-terminal portion of helix alpha1 (residues 13-17) and in helix alpha2 as well as a turn and nonnative hydrophobic clustering between beta3 and beta4 were observed, consistent with NMR data. In addition, there were tertiary contacts between residues in alpha1 and the C-terminal portion of the beta-sheet. The simulated structures allow the rudimentary NMR data to be fleshed out. The consistency between simulation and experiment inspires confidence in the methods. A description of the folding pathway of barnase from the denatured to the native state can be constructed by combining the simulation with experimental data from phi value analysis and NMR.