Structural Order of Water Molecules around Hydrophobic Solutes: Length-Scale Dependence and Solute-Solvent Coupling

It has been suggested that the structure and thermodynamics of the water molecules in the hydration layer of simple hydrophobic solutes undergo an order-disorder transition around a nanometer length-scale of the solute size. Using extensive atomistic molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulation studies, we have probed this order-disorder transition around model hydrophobic solutes of varying size and shape (spherical, planar, and linear), as well as flexible hydrophobic homopolymer chains (n-alkanes), where the conformational fluctuations are likely to create both spatial and temporal heterogeneity on the solvent accessible surface. We have explored the structural response of the water molecules in the hydration shell due to the local variations of the length-scale (or curvature) upon hydrophobic collapse and/or local conformational changes of these polymers. We have shown that the tetrahedral order of the water molecules in the hydration shell is practically independent of the polymer size in the extended state of the polymer due to the availability of a subnanometer cross-sectional length-scale, allowing the water molecules to form hydrogen bonds around the polymer chain. Beyond a certain length of the polymer chains, the collapsed states (associated with larger solute length-scale) start to induce disorder in the surface water molecules. We demonstrate that the local structure (both local number density and tetrahedral order) of the hydration layer is dynamically coupled to the local topology of the polymer. Thus, we envisage that in a flexible (bio)polymer, the hydration shell properties will be sensitive to the local conformational state of the molecule (both spatially and temporally), and the overall observed water structure and dynamics will be dependent on the topological/chemical heterogeneity, and the time-scale of fluctuations in the local curvature (length-scale) of the solvent accessible surface. Moreover, we have demonstrated the direct coupling between the local density fluctuations of water and the local hydrophobic collapse of the polymer. For the extended state of the polymer, the local solvent density fluctuation is practically independent of the solute coordinate (length-scale), and the hydrophobic collapse of the polymer is prompted by a "local dewetting" process induced by these fluctuations.