Hydration Thermodynamics of Non-Polar Aromatic Hydrocarbons: Comparison of Implicit and Explicit Solvation Models

Thermodynamics of π–π Interactions of Benzene and Phenol in Water

The π–π interaction is a major driving force that stabilizes protein assemblies during protein folding. Recent studies have additionally demonstrated its involvement in the liquid–liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs). As the participating residues in IDPs are exposed to water, π–π interactions for LLPS must be modeled in water, as opposed to the interactions that are often established at the hydrophobic domains of folded proteins. Thus, we investigated the association of free energies of benzene and phenol dimers in water by integrating van der Waals (vdW)-corrected density functional theory (DFT) and DFT in classical explicit solvents (DFT-CES). By comparing the vdW-corrected DFT and DFT-CES results with high-level wavefunction calculations and experimental solvation free energies, respectively, we established the quantitative credibility of these approaches, enabling a reliable prediction of the benzene and phenol dimer association free energies in water. We discovered that solvation influences dimer association free energies, but not significantly when no direct hydrogen-bond-type interaction exists between two monomeric units, which can be explained by the enthalpy–entropy compensation. Our comprehensive computational study of the solvation effect on π–π interactions in water could help us understand the molecular-level driving mechanism underlying the IDP phase behaviors.

Structure, Dynamics, and Wettability of Water at Metal Interfaces

The water/metal interface often governs important chemophysical processes in various technologies. Therefore, from scientific and engineering perspectives, the detailed molecular-level elucidation of the water/metal interface is of high priority, but the related research is limited. In experiments, the surface-science techniques, which can provide full structural details of the surface, are not easy to directly apply to the interfacial systems under ambient conditions, and the well-defined facets cannot be entirely free from contamination at the contact with water. To answer long-standing debates regarding the wettability, structure, and dynamics of water at metal interfaces, we here develop reliable first-principles-based multiscale simulations. Using the state-of-the-art simulations, we find that the clean metal surfaces are actually superhydrophilic and yield zero contact angles. Furthermore, we disclose an inadequacy of widespread ice-like bilayer model of the water adlayers on metal surfaces from both averaged structural and dynamic points of view. Our findings on the nature of water on metal surfaces provide new molecular level perspectives on the tuning and design of water/metal interfaces that are at the heart of many energy applications.

Enthalpy-Entropy Interplay in π-Stacking Interaction of Benzene Dimer in Water

Aromatic groups can engage in an interesting class of noncovalent interactions termed π-π interactions, which play a pivotal role in stabilizing a variety of molecular architectures, including nucleic acids, proteins, and supramolecular assemblies. When the aromatic compounds interact with each other in an aqueous environment, their association is facilitated by the hydrophobic effect-the trend of nonpolar solutes to aggregate in a polar solution. To develop an in-depth understanding of hydrophobic association, we investigate in the present work π-π interactions in water, employing as a paradigm the benzene dimer. Using DFT-CES, a mean-field QM/MM method recently developed by our group, we describe the benzene solute at a quantum-mechanical level. Full consideration of detailed solute-electron density enables an optimal description of the solute-solvent interactions, leading to an accurate prediction of hydration free energies. In π-stacking of benzene, we find an entropic stabilization associated with the shrinkage of the solvent-excluded volume, which agrees with the theory of hydrophobic effect at subnanoscales. However, at the equilibrium binding distance of the benzene dimer, we find that the entropic stabilization nearly cancels out due to the enthalpic cost required for dewetting the internal space. Such an enthalpy-entropy compensation leads the association free energy to be predominantly dictated by the solute-solute interaction enthalpy. The present work offers new insight into the mechanistic role of water and the primary thermodynamic driving force of hydrophobic association.