Hydrophobicity of proteins and nanostructured solutes is governed by topographical and chemical context.

Hydrophobic interactions drive many important biomolecular self-assembly phenomena. However, characterizing hydrophobicity at the nanoscale has remained a challenge due to its nontrivial dependence on the chemistry and topography of biomolecular surfaces. Here we use molecular simulations coupled with enhanced sampling methods to systematically displace water molecules from the hydration shells of nanostructured solutes and calculate the free energetics of interfacial water density fluctuations, which quantify the extent of solute-water adhesion, and therefore solute hydrophobicity. In particular, we characterize the hydrophobicity of curved graphene sheets, self-assembled monolayers (SAMs) with chemical patterns, and mutants of the protein hydrophobin-II. We find that water density fluctuations are enhanced near concave nonpolar surfaces compared with those near flat or convex ones, suggesting that concave surfaces are more hydrophobic. We also find that patterned SAMs and protein mutants, having the same number of nonpolar and polar sites but different geometrical arrangements, can display significantly different strengths of adhesion with water. Specifically, hydroxyl groups reduce the hydrophobicity of methyl-terminated SAMs most effectively not when they are clustered together but when they are separated by one methyl group. Hydrophobin-II mutants show that a charged amino acid reduces the hydrophobicity of a large nonpolar patch when placed at its center, rather than at its edge. Our results highlight the power of water density fluctuations-based measures to characterize the hydrophobicity of nanoscale surfaces and caution against the use of additive approximations, such as the commonly used surface area models or hydropathy scales for characterizing biomolecular hydrophobicity and the associated driving forces of assembly.

Water-mediated ion-ion interactions are enhanced at the water vapor-liquid interface.

There is overwhelming evidence that ions are present near the vapor-liquid interface of aqueous salt solutions. Charged groups can also be driven to interfaces by attaching them to hydrophobic moieties. Despite their importance in many self-assembly phenomena, how ion-ion interactions are affected by interfaces is not understood. We use molecular simulations to show that the effective forces between small ions change character dramatically near the water vapor-liquid interface. Specifically, the water-mediated attraction between oppositely charged ions is enhanced relative to that in bulk water. Further, the repulsion between like-charged ions is weaker than that expected from a continuum dielectric description and can even become attractive as the ions are drawn to the vapor side. We show that thermodynamics of ion association are governed by a delicate balance of ion hydration, interfacial tension, and restriction of capillary fluctuations at the interface, leading to nonintuitive phenomena, such as water-mediated like charge attraction. "Sticky" electrostatic interactions may have important consequences on biomolecular structure, assembly, and aggregation at soft liquid interfaces. We demonstrate this by studying an interfacially active model peptide that changes its structure from α-helical to a hairpin-turn-like one in response to charging of its ends.

Structure and dynamics of single hydrophobic/ionic heteropolymers at the vapor-liquid interface of water.

We focus on the conformational stability, structure, and dynamics of hydrophobic/charged homopolymers and heteropolymers at the vapor-liquid interface of water using extensive molecular dynamics simulations. Hydrophobic polymers collapse into globular structures in bulk water but unfold and sample a broad range of conformations at the vapor-liquid interface of water. We show that adding a pair of charges to a hydrophobic polymer at the interface can dramatically change its conformations, stabilizing hairpinlike structures, with molecular details depending on the location of the charged pair in the sequence. The translational dynamics of homopolymers and heteropolymers are also different, whereas the homopolymers skate on the interface with low drag, the tendency of charged groups to remain hydrated pulls the heteropolymers toward the liquid side of the interface, thus pinning them, increasing drag, and slowing the translational dynamics. The conformational dynamics of heteropolymers are also slower than that of the homopolymer and depend on the location of the charged groups in the sequence. Conformational dynamics are most restricted for the end-charged heteropolymer and speed up as the charge pair is moved toward the center of the sequence. We rationalize these trends using the fundamental understanding of the effects of the interface on primitive pair-level interactions between two hydrophobic groups and between oppositely charged ions in its vicinity.

Trimethylamine N-oxide (TMAO) and tert-butyl alcohol (TBA) at hydrophobic interfaces: insights from molecular dynamics simulations.

TMAO, a potent osmolyte, and TBA, a denaturant, have similar molecular architecture but somewhat different chemistry. We employ extensive molecular dynamics simulations to quantify their behavior at vapor-water and octane-water interfaces. We show that interfacial structure-density and orientation-and their dependence on solution concentration are markedly different for the two molecules. TMAO molecules are moderately surface active and adopt orientations with their N-O vector approximately parallel to the aqueous interface. That is, not all methyl groups of TMAO at the interface point away from the water phase. In contrast, TBA molecules act as molecular amphiphiles, are highly surface active, and, at low concentrations, adopt orientations with their methyl groups pointing away and the C-O vector pointing directly into water. The behavior of TMAO at aqueous interfaces is only weakly dependent on its solution concentration, whereas that of TBA depends strongly on concentration. We show that this concentration dependence arises from their different hydrogen bonding capabilities-TMAO can only accept hydrogen bonds from water, whereas TBA can accept (donate) hydrogen bonds from (to) water or other TBA molecules. The ability to self-associate, particularly visible in TBA molecules in the interfacial layer, allows them to sample a broad range of orientations at higher concentrations. In light of the role of TMAO and TBA in biomolecular stability, our results provide a reference with which to compare their behavior near biological interfaces. Also, given the ubiquity of aqueous interfaces in biology, chemistry, and technology, our results may be useful in the design of interfacially active small molecules with the aim to control their orientations and interactions.