We evaluate the molecular response of hydration water as a function of temperature and proximity to the surface of the peptide N-acetyl-leucine-methylamide (NALMA) when in the presence of the kosmotrope co-solvent glycerol or the chaotrope co-solvent dimethyl sulfoxide (DMSO), using molecular dynamics simulation with a polarizable force field. These detailed microscopic studies complement established thermodynamic analysis on the role of co-solvents in shifting the equilibrium for proteins away from or towards the native folded state. We find that the structure of the water at the peptide interfaces reflects an increase in hydration number in the glycerol solution and a decrease in hydration numbers in the DMSO solution. While the water dynamics around NALMA in the presence of both co-solvents is slower than that observed with the water solvent alone, in the DMSO mixture we no longer measure a separation in water motion time scales at low temperatures as is seen in the pure water solvent, but rather one single relaxation time. In the glycerol, however, we do observe a separation of time scales at low temperatures, supporting the hypothesis that hydration water near a hydrophobic solute evolves on a separate time scale than the extensive hydrogen-bonding network of more bulk-like water. Our simulation studies highlight the differences in the two co-solvent solutions due to the relative frequency of water contacts with the hydrophobic vs. hydrophilic peptide surface, and direct water interactions with the co-solvents.

We report quasi-elastic neutron scattering experiments to contrast the water dynamics as a function of temperature for hydrophilic and amphiphilic peptides under the same level of confinement, as models for understanding hydration dynamics near chemically heterogeneous protein surfaces. We find that the hydrophilic peptide shows only a single non-Arrhenius translational process with no evidence of spatial heterogeneity unlike the amphiphilic peptide solution that exhibits two translational relaxations with an Arrhenius and non-Arrhenius dependence on temperature. Together these results provide experimental proof that heterogeneous dynamical signatures near protein surfaces arise in part from chemical heterogeneity (energy disorder) as opposed to mere topological roughness of the protein surface.

We perform neutron diffraction and quasi-elastic neutron scattering (QENS) to probe hydration water structure, and dynamics down to supercooled temperatures, of a concentrated amphiphilic peptide system with the co-solvents glycerol and dimethyl sulfoxide. We find that the kosmotropic co-solvent glycerol preserves the hydration structure near the peptide that is observed in the water solvent alone, that in turn preserves the dynamical temperature trends of two water relaxation processes--one corresponding to a localized relaxation process of the peptide bound surface water and a second relaxation process of the outer hydration layers. By contrast the chaotropic co-solvent, by disrupting the hydration layer near the peptide surface, eliminates the inner hydration layer relaxation process induced by the peptide, to show a single time scale for translational water dynamics.

We performed classical molecular dynamics simulations using both fixed-charge and polarizable water and protein force fields to contrast the hydration dynamics near hydrophilic and amphiphilic peptides as a function of temperature. The high peptide concentrations we use serve as a model for the surface of folded proteins where hydration layers around each residue overlap significantly. Through simulation we determine that there are notable differences in the water dynamics analyzed from the outer and inner hydration layer regions of the amphiphilic peptide solution that explains the experimentally observed presence of two translational relaxations, while the hydrophilic peptide solution shows only a single non-Arrhenius translational process with no distinction between hydration layers. Given that water dynamics for the amphiphilic peptide system reproduces all known rotational and translational hydration dynamical anomalies exhibited by hydration water near protein surfaces, our analysis provides strong evidence that dynamical signatures near biological interfaces arises because of frustration in the hydration dynamics induced by chemical heterogeneity, as opposed to just topological roughness, of the protein surface.