Protein hydration elucidated by molecular dynamics simulation

Abstract

Molecular dynamics (MD) simulation covering a wide range of hydration indicate that myoglobin is fully hydrated by 350 water molecules, in agreement with experiment. These waters, originally placed uniformly about the protein, form clusters that hydrate every charged group throughout the entire simulation. Some atoms in charged groups are hydrated by two water layers while 37% of the protein surface remains uncovered. The locations of the 350 waters are consistent with those of crystallographic waters resolved by x-ray and neutron diffraction. Hydration by 350 waters at 300 K stabilizes the conformation of carboxymyoglobin measured by x-ray diffraction throughout the entire protein, halves the rate of torsional transitions, and promotes alternative conformations for surface atoms. The glass transition observed experimentally in hydrated myoglobin near 220 K is also seen in the simulations and correlates with an increase in the number of dihedral angles undergoing transitions. The anharmonic protein motion above 220 K is enhanced by protein hydration.