The dominant interaction between peptide and urea is electrostatic in nature: a molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea-peptide compared to water-peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.