Water structuring above solutes with planar hydrophobic surfaces

Understanding the physics of hydrophobic solvation

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations. The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones solvent, and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid-vapor coexistence in the macroscopic planar limit, i.e., as the solute radius R_s → ∞. Focusing on the radial density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), our binding potential analysis provides explicit predictions for the manner in which the key features of ρ(r) and χ(r) scale with R_s, the strength of solute-water attraction ɛ_sf, and the deviation from liquid-vapor coexistence of the chemical potential, δμ. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such, our theory provides a firm basis for understanding the physics of hydrophobic solvation.

A Form of Non-Volatile Solid-like Hexadecane Found in Micron-Scale Silica Microtubule

Anomalous solid-like liquids at the solid-liquid interface have been recently reported. The mechanistic factors contributing to these anomalous liquids and whether they can stably exist at high vacuum are interesting, yet unexplored, questions. In this paper, thin slices of silica tubes soaked in hexadecane were observed under a transmission electron microscope at room temperature. The H-spectrum of hexadecane in the microtubules was measured by nuclear magnetic resonance. On the interior surface of these silica tubes, 0.2-30 μm in inside diameter (ID), a layer (12-400 nm) of a type of non-volatile hexadecane was found with thickness inversely correlated with the tube ID. A sample of this anomalous hexadecane in microtubules 0.4 μm in ID was found to be formable by an ion beam. Compared with the nuclear magnetic resonance H-spectroscopy of conventional hexadecane, the characteristic peaks of this abnormal hexadecane were shifted to the high field with a broader characteristic peak, nuclear magnetic resonance hydrogen spectroscopy spectral features typical of that of solids. The surface density of these abnormal hexadecanes was found to be positively correlated with the silanol groups found on the interior silica microtubular surface. This positive correlation indicates that the high-density aggregation of silanol is an essential factor for forming the abnormal hexadecane reported in this paper.

Density Depletion and Enhanced Fluctuations in Water near Hydrophobic Solutes: Identifying the Underlying Physics

We investigate the origin of the density depletion and enhanced density fluctuations that occur in water in the vicinity of an extended hydrophobic solute. We argue that both phenomena are remnants of the critical drying surface phase transition that occurs at liquid-vapor coexistence in the macroscopic planar limit, i.e., as the solute radius R_{s}→∞. Focusing on the density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), we develop a scaling theory which expresses the extent of the density depletion and enhancement in compressibility in terms of R_{s}, the strength of solute-water attraction ϵ_{s}, and the deviation from liquid-vapor coexistence δμ. Testing the predictions against results of classical density functional theory for a simple solvent and grand canonical Monte Carlo simulations of a popular water model, we find that the theory provides a firm physical basis for understanding how water behaves at a hydrophobe.

Effects of the Chemical and Structural Properties of Silane Monolayers on the Organization of Water Molecules and Ions at Interfaces, from Molecular Dynamics Simulations

Understanding the organization of the hydration layer at functionalized silica surfaces is relevant for a large range of biosensing applications or surface phenomena such as biomolecule adsorption. Silane monolayers are widely used to functionalize silica surfaces. Using molecular dynamics simulations, we have investigated the role of silane molecule head-group charge, alkyl chain length, and surface coverage in the structural organization and dynamic properties of Na⁺ ions, Cl^- ions, and water molecules at the interface. The silane molecules studied are 3-aminopropyldimethylethoxysilane, n-propyldimethylmethoxysilane, octadecyldimethylmethoxysilane, and (dimethylamino)dimethylsilylundecanoate. Our results suggest that the distribution of interfacial ions is sensitive to the 2D dispersion of the silane-charged head groups. Also, while charged silane monolayers show a strong orientation of interfacial water molecules, which leads to a rupture in the hydrogen bond network and disturbs their tetrahedral organization, the arrangement of water molecules at the interface with uncharged silane monolayers seems to be related to the surface roughness and to alkyl chain length. In line with these results, the diffusion of ions and water molecules is higher at the CH₃ long monolayer interface than at the CH₃ short monolayer interface and at the charged monolayer interfaces. Also, whatever the silane molecules studied, bulk properties are recovered around 0.7 nm above the interface. The interfacial water organization is known to impact biomolecule adsorption. Therefore, these results could further help in optimizing the functionalization layers to capture analytes.

Insignificant Effect of Temperature on the Structure and Angular Jumps of Water near a Hydrophobic Cation

The ambiguity in the behavior of water molecules around hydrophobic solutes is a matter of interest for many studies. Motivated by the earlier results on the dynamics of water molecules around tetramethylammonium (TMA) cation, we present the effect of temperature on the structure and angular jumps of water due to hydrophobicity using first principles molecular dynamics simulations. The average intermolecular distance between the central oxygen and four nearest neighbors is found to be the highest for water molecules in the solvation shell of TMA at 400 K, followed by the same at 330 K. The hydrogen bond (HB) donor-acceptor count, HB per water molecule, and tetrahedral order parameter suggests the loss of tetrahedrality in the solvation shell. Elevated temperature affects the tetrahedral parameter in local regions. The HB jump mechanism is studied for methyl hydrogen and water molecules in the solvation shell. Observations hint at the presence of dangling water molecules in the vicinity of the hydrophobic cation, and no evidence is found for the enhanced structural ordering of nearby water molecules.

Physico-chemical properties of aqueous drug solutions: From the basic thermodynamics to the advanced experimental and simulation results

The physical chemical properties of aqueous solutions of model compounds are illustrated in relation to hydration and solubility issues by using three perspectives: thermodynamic, spectroscopic and molecular dynamics simulations. The thermodynamic survey of the fundamental backgrounds of concentration dependence and experimental solubility results show some peculiar behavior of aqueous solutions with several types of similar solutes. Secondly, the use of a variety of experimental spectroscopic devices, operating under different experimental conditions of dimension and frequency, has produced a large amount of structural and dynamic data on aqueous solutions showing the richness of the information produced, depending on where and how the experiment is carried out. Finally, the use of molecular dynamics computational work is presented to highlight how the different types of solute functional groups and surface topologies organize adjacent water molecules differently. The highly valuable contribution of computer simulation studies in providing molecular explanations for experimental deductions, either of a thermodynamic or spectroscopic nature, is shown to have changed the current knowledge of many aqueous solution processes. While this paper is intended to provide a collective view on the latest literature results, still the presentation aims at a tutorial explanation of the potentials of the three methodologies in the field of aqueous solutions of pharmaceutical molecules.