Abstract
We have performed 20ns of fully atomistic molecular dynamics simulations of Hen Egg-White Lysozyme in 0, 10, 20, 30 and 100% by weight of glycerol in water to better understand the microscopic physics behind the bioprotection offered by glycerol to naturally occuring biological systems. The sovlent exposure of protein surface residues changes when glycerol is introduced. The dynamic behavior of the protein, as quantified by the Incoherent Intermediate Scattering Function, shows a non-monotonic dependence on glycerol content. The fluctuations of the protein residues with respect to each other were found to be similar in all water– containing solvents; but different from the pure glycerol case. The increase in the number of protein–glycerol hydrogen bonds in glycerol–water binary mixtures explains the slowing down of protein dynamics as the glycerol content increases. We also explored the dynamic behavior of the hydration layer. We show that the short-length scale dynamics of this layer are insenstive to glycerol concentration. However, the long-length scale behavior shows a significant dependence on glycerol content. We also provide insights into the behavior of bound and mobile water molecules.