An understanding of the hydrophobicity of complex heterogeneous molecular assemblies is crucial to characterize and predict interactions between biomolecules. As such, uncovering the subtleties of assembly processes hinges on an accurate classificati
on of the relevant interfaces involved, and much effort has been spent on developing so-called hydrophobicity maps. In this work, we introduce a novel electrostatics-based mapping of aqueous interfaces that focuses on the collective, long-wavelength electrostatic response of water to the presence of nearby surfaces. In addition to distinguishing between hydrophobic and hydrophilic regions of heterogeneous surfaces, this electrostatic mapping can also differentiate between hydrophilic regions that polarize nearby waters in opposing directions. We therefore expect this approach to find use in predicting the location of possible water-mediated hydrophilic interactions, in addition to the more commonly emphasized hydrophobic interactions that can also be of significant importance.
Liquid water can become metastable with respect to its vapor in hydrophobic confinement. The resulting dewetting transitions are often impeded by large kinetic barriers. According to macroscopic theory, such barriers arise from the free energy requir
ed to nucleate a critical vapor tube that spans the region between two hydrophobic surfaces - tubes with smaller radii collapse, whereas larger ones grow to dry the entire confined region. Using extensive molecular simulations of water between two nanoscopic hydrophobic surfaces, in conjunction with advanced sampling techniques, here we show that for inter-surface separations that thermodynamically favor dewetting, the barrier to dewetting does not correspond to the formation of a (classical) critical vapor tube. Instead, it corresponds to an abrupt transition from an isolated cavity adjacent to one of the confining surfaces to a gap-spanning vapor tube that is already larger than the critical vapor tube anticipated by macroscopic theory. Correspondingly, the barrier to dewetting is also smaller than the classical expectation. We show that the peculiar nature of water density fluctuations adjacent to extended hydrophobic surfaces - namely, the enhanced likelihood of observing low-density fluctuations relative to Gaussian statistics - facilitates this non-classical behavior. By stabilizing isolated cavities relative to vapor tubes, enhanced water density fluctuations thus stabilize novel pathways, which circumvent the classical barriers and offer diminished resistance to dewetting. Our results thus suggest a key role for fluctuations in speeding up the kinetics of numerous phenomena ranging from Cassie-Wenzel transitions on superhydrophobic surfaces, to hydrophobically-driven biomolecular folding and assembly.
We use appropriately defined short ranged reference models of liquid water to clarify the different roles local hydrogen bonding, van der Waals attractions, and long ranged electrostatic interactions play in the solvation and association of apolar so
lutes in water. While local hydrogen bonding in- teractions dominate hydrophobic effects involving small solutes, longer ranged electrostatic and dis- persion interactions are found to be increasingly important in the description of interfacial structure around large solutes. The hydrogen bond network sets the solute length scale at which a crossover in solvation behavior between these small and large length scale regimes is observed. Unbalanced long ranged forces acting on interfacial water molecules are also important in hydrophobic association, illustrated here by analysis of the association of model methane and buckminsterfullerene solutes.
An understanding of density fluctuations in bulk water has made significant contributions to our understanding of the hydration and interactions of idealized, purely repulsive hydrophobic solutes. To similarly inform the hydration of realistic hydrop
hobic solutes that have dispersive interactions with water, here we characterize water density fluctuations in the presence of attractive fields that correspond to solute-water attractions. We find that when the attractive field acts only in the solute hydration shell, but not in the solute core, it does not significantly alter water density fluctuations in the solute core region. We further find that for a wide range of solute sizes and attraction strengths, the free energetics of turning on the attractive fields in bulk water are accurately captured by linear response theory. Our results also suggest strategies for more efficiently estimating hydration free energies of realistic solutes in bulk water and at interfaces.
