Molecular dynamics simulations of the protein chymotrypsin inhibitor 2 in 8 M urea at 60°C were undertaken to investigate the molecular basis of chemical denaturation. The protein unfolded rapidly under these conditions, but it retained its native structure in a control simulation in water at the same temperature. The overall process of unfolding in urea was similar to that observed in thermal denaturation simulations above the protein's T _m of 75°C. The first step in unfolding was expansion of the hydrophobic core. Then, the core was solvated by water and later by urea. The denatured structures in both urea and at high temperature contained residual native helical structure, whereas the β-structure was completely disrupted. The average residence time for urea around hydrophilic groups was six times greater than around hydrophobic residues and in all cases greater than the corresponding water residence times. Water self-diffusion was reduced 40% in 8 M urea. Urea altered water structure and dynamics, thereby diminishing the hydrophobic effect and encouraging solvation of hydrophobic groups. In addition, through urea's weakening of water structure, water became free to compete with intraprotein interactions. Urea also interacted directly with polar residues and the peptide backbone, thereby stabilizing nonnative conformations. These simulations suggest that urea denatures proteins via both direct and indirect mechanisms.

The ability to diagnose and treat prion diseases is limited by our current understanding of the conversion process of the protein from healthy to harmful isoform. Whereas the monomeric, benign species is well characterized, the misfolded conformations responsible for infectivity and neurodegeneration remain elusive. There is mounting evidence that fibrillization intermediates, or protofibrils, but not mature fibrils or plaques, are the pathogenic species in amyloid diseases. Here, we use molecular dynamics to simulate the conversion of the prion protein. Molecular dynamics simulation produces a scrapie prion protein-like conformation enriched in β-structure that is in good agreement with available experimental data. The converted conformation was then used to model a protofibril by means of the docking of hydrophobic patches of the template structure to form hydrogen-bonded sheets spanning adjacent subunits. The resulting protofibril model provides a non-branching aggregate with a 3 ₁ axis of symmetry that is in good agreement with a wide variety of experimental data; importantly, it was derived from realistic simulation of the conversion process.

Proteins are very sensitive to their solvent environments. Urea is a common chemical denaturant of proteins, yet some animals contain high concentrations of urea. These animals have evolved an interesting mechanism to counteract the effects of urea by using trimethylamine N -oxide (TMAO). The molecular basis for the ability of TMAO to act as a chemical chaperone remains unknown. Here, we describe molecular dynamics simulations of a small globular protein, chymotrypsin inhibitor 2, in 8 M urea and 4 M TMAO/8 M urea solutions, in addition to other control simulations, to investigate this effect at the atomic level. In 8 M urea, the protein unfolds, and urea acts in both a direct and indirect manner to achieve this effect. In contrast, introduction of 4 M TMAO counteracts the effect of urea and the protein remains well structured. TMAO makes few direct interactions with the protein. Instead, it prevents unfolding of the protein by structuring the solvent. In particular, TMAO orders the solvent and discourages it from competing with intraprotein H bonds and breaking up the hydrophobic core of the protein.

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions.