Hydrophobic interactions play a central role in bioinspired strategies for molecular self-assembly in water, yet how these interactions are encoded by chemically heterogeneous interfaces is poorly understood. We report an experimental investigation of the influence of immobilized polar groups (amine) and cations (ammonium and guanidinium) on enthalpic and entropic contributions to hydrophobic interactions mediated by methyl-terminated surfaces at temperatures ranging from 298 K to 328 K and pH values between 3.5 to 10.5. We use our measurements to calculate the change in free energy (and enthalpic and entropic components) that accompanies transfer of each surface from aqueous TEA containing 60 vol% methanol into aqueous TEA ( i.e. , transfer free energy that characterizes hydrophobicity). We find the thermodynamic signature of the pure methyl surface (positive transfer enthalpy and entropy) to be altered qualitatively by incorporation of amine or guanidinium groups into the surface (negative transfer enthalpy and near zero transfer entropy). In contrast, ammonium groups immobilized on a methyl surface do not change the thermodynamic signature of the hydrophobic interaction. Compensation of entropy and enthalpy is clearly evident in our results, but the overall trends in the transfer free energies are dominated by enthalpic effects. This observation and others lead us to hypothesize that the dominant effect of the immobilized charged or polar groups in our experiments is to influence the number or strength of hydrogen bonds formed by interfacial water molecules adjacent to the nonpolar domains. Overall, these results provide insight into entropy–enthalpy compensation at chemically heterogeneous surfaces, and generate hypotheses and a rich experimental dataset for further exploration via simulation.
more »
« less
Structural features of interfacial water predict the hydrophobicity of chemically heterogeneous surfaces
The hydrophobicity of an interface determines the magnitude of hydrophobic interactions that drive numerous biological and industrial processes. Chemically heterogeneous interfaces are abundant in these contexts; examples include the surfaces of proteins, functionalized nanomaterials, and polymeric materials. While the hydrophobicity of nonpolar solutes can be predicted and related to the structure of interfacial water molecules, predicting the hydrophobicity of chemically heterogeneous interfaces remains a challenge because of the complex, non-additive contributions to hydrophobicity that depend on the chemical identity and nanoscale spatial arrangements of polar and nonpolar groups. In this work, we utilize atomistic molecular dynamics simulations in conjunction with enhanced sampling and data-centric analysis techniques to quantitatively relate changes in interfacial water structure to the hydration free energy (a thermodynamically well-defined descriptor of hydrophobicity) of chemically heterogeneous interfaces. We analyze a large data set of 58 self-assembled monolayers (SAMs) composed of ligands with nonpolar and polar end groups of different chemical identity (amine, amide, and hydroxyl) in five mole fractions, two spatial patterns, and with scaled partial charges. We find that only five features of interfacial water structure are required to accurately predict hydration free energies. Examination of these features reveals mechanistic insights into the interfacial hydrogen bonding behaviors that distinguish different surface compositions and patterns. This analysis also identifies the probability of highly coordinated water structures as a unique signature of hydrophobicity. These insights provide a physical basis to understand the hydrophobicity of chemically heterogeneous interfaces and connect hydrophobicity to experimentally accessible perturbations of interfacial water structure.
more »
« less
- PAR ID:
- 10398083
- Date Published:
- Journal Name:
- Chemical Science
- Volume:
- 14
- Issue:
- 5
- ISSN:
- 2041-6520
- Page Range / eLocation ID:
- 1308 to 1319
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Triacylglycerols (TG) are the primary neutral lipids in lipid droplets (LDs), organelles responsible for lipid storage, metabolism, and signaling. Molecular dynamics (MD) simulations have provided valuable insight into LD structure, but fixed-charge force fields struggle to capture TG behavior across both hydrophobic cores and polar interfaces. Here, we develop and evaluate a polarizable TG model using the Drude2023 lipid force field and benchmark its performance against experimental measurements of bulk density, TG−water interfacial tension, core hydration, and monolayer expansion. The Drude model accurately reproduces the experimental properties and captures key monolayer features such as surface-oriented TGs (SURF-TGs) and chemically distinct membrane packing defects. Compared to fixed-charge models such as C36-standard and C36-cutoff, the Drude polarizable model is the only force field able to capture the dual nature of TG at polar−nonpolar interfaces like the LD monolayer and more homogeneous hydrophobic environments, like the LD core. However, C36-standard is consistent with the Drude results for the LD monolayer, while C36-cutoff is consistent with the decreased hydration in the LD core. Even with large applied surface tensions, C36-cutoff does not produce Drude-like LD monolayer properties. These results highlight the importance of dynamic polarizability and establish Drude2023 as a more reliable framework for simulating TG in heterogeneous systems like LDs.more » « less
-
Hydrophobic interactions drive the binding of nonpolar ligands to the oily pockets of proteins and supramolecular species in aqueous solution. As such, the wetting of host pockets is expected to play a critical role in determining the thermodynamics of guest binding. Here we use molecular simulations to examine the impact of pressure on the wetting and dewetting of the nonpolar pockets of a series of deep-cavity cavitands in water. The portals to the cavitand pockets are functionalized with both nonpolar (methyl) and polar (hydroxyl) groups oriented pointing either upward or inward toward the pocket. We find wetting of the pocket is favored by the hydroxyl groups and dewetting is favored by the methyl groups. The distribution of waters in the pocket is found to exhibit a two-state-like equilibrium between wet and dry states with a free energy barrier between the two states. Moreover, we demonstrate that the pocket hydration of the cavitands can be collapsed onto a unified adsorption isotherm by assuming the effective pressures within each cavitand pocket differ by a shift pressure that depends on the chemical identity and number of functional groups placed about the portal. These observations support the development of a twostate capillary evaporation model that accurately describes the equilibrium between states and naturally gives rise to the effective shift pressures observed from simulation. This work demonstrates that the hydration of host pockets can be tuned following simple design rules that in turn are expected to impact the thermodynamics of guest complexation.more » « less
-
Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.more » « less
-
The aversion of hydrophobic solutes for water drives diverse interactions and assemblies across materials science, biology, and beyond. Here, we review the theoretical, computational, and experimental developments that underpin a contemporary understanding of hydrophobic effects. We discuss how an understanding of density fluctuations in bulk water can shed light on the fundamental differences in the hydration of molecular and macroscopic solutes; these differences, in turn, explain why hydrophobic interactions become stronger upon increasing temperature. We also illustrate the sensitive dependence of surface hydrophobicity on the chemical and topographical patterns the surface displays, which makes the use of approximate approaches for estimating hydrophobicity particularly challenging. Importantly, the hydrophobicity of complex surfaces, such as those of proteins, which display nanoscale heterogeneity, can nevertheless be characterized using interfacial water density fluctuations; such a characterization also informs protein regions that mediate their interactions. Finally, we build upon an understanding of hydrophobic hydration and the ability to characterize hydrophobicity to inform the context-dependent thermodynamic forces that drive hydrophobic interactions and the desolvation barriers that impede them.more » « less
An official website of the United States government

