Analysis of the Electron Density of a Water Molecule Encapsulated by Two Cholic Acid Residues

Cholic acid is a trihydroxy bile acid with a nice peculiarity: the average distance between the oxygen atoms (O7 and O12) of the hydroxy groups located at C7 and C12 carbon atoms is 4.5 Å, a value which perfectly matches with the O/O tetrahedral edge distance in Ih ice. In the solid phase, they are involved in the formation of hydrogen bonds with other cholic acid units and solvents. This fact was satisfactorily used for designing a cholic dimer which encapsulates one single water molecule between two cholic residues, its oxygen atom (Ow) being exactly located at the centroid of a distorted tetrahedron formed by the four steroid hydroxy groups. The water molecule participates in four hydrogen bonds, with the water simultaneously being an acceptor from the 2 O12 (hydrogen lengths are 2.177 Å and 2.114 Å) and a donor towards the 2 O7 (hydrogen bond lengths are 1.866 Å and 1.920 Å). These facts suggest that this system can be a nice model for the theoretical study of the formation of ice-like structures. These are frequently proposed to describe the water structure found in a plethora of systems (water interfaces, metal complexes, solubilized hydrophobic species, proteins, and confined carbon nanotubes). The above tetrahedral structure is proposed as a reference model for those systems, and the results obtained from the application of the atoms in molecules theory are presented here. Furthermore, the structure of the whole system allows a division into two interesting subsystems in which water is the acceptor of one hydrogen bond and the donor of another. The analysis of the calculated electron density is performed through its gradient vector and the Laplacian. The calculation of the complexation energy used correction of the basis set superposition error (BSSE) with the counterpoise method. As expected, four critical points located in the H…O bond paths were identified. All calculated parameters obey the proposed criteria for hydrogen bonds. The total energy for the interaction in the tetrahedral structure is 54.29 kJ/mol, while the summation obtained of the two independent subsystems and the one between the alkyl rings without water is only 2.5 kJ/mol higher. This concordance, together with the calculated values for the electron density, the Laplacian of the electron density, and the lengths of the oxygen atom and the hydrogen atom (involved in the formation of each hydrogen bond) to the hydrogen bond critical point, suggests that each pair of hydrogen bonds can be considered independent of each other.

Investigating water/oil interfaces with opto-thermophoresis

Charging of interfaces between water and hydrophobic media is a mysterious feature whose nature and origin have been under debate. Here, we investigate the fundamentals of the interfacial behaviors of water by employing opto-thermophoretic tweezers to study temperature-gradient-induced perturbation of dipole arrangement at water/oil interfaces. With surfactant-free perfluoropentane-in-water emulsions as a model interface, additional polar organic solvents are introduced to systematically modify the structural aspects of the interface. Through our experimental measurements on the thermophoretic behaviors of oil droplets under a light-generated temperature gradient, in combination with theoretical analysis, we propose that water molecules and mobile negative charges are present at the water/oil interfaces with specific dipole arrangement to hydrate oil droplets, and that this arrangement is highly susceptible to the thermal perturbation due to the mobility of the negative charges. These findings suggest a potential of opto-thermophoresis in probing aqueous interfaces and could enrich understanding of the interfacial behaviors of water.

Water as the reaction medium in organic chemistry: from our worst enemy to our best friend

A review presenting water as the logical reaction medium for the future of organic chemistry. A discussion is offered that covers both the "on water" and "in water" phenomena, and how water is playing unique roles in each, specifically with regard to its use in organic synthesis.

Water Structure and Properties at Hydrophilic and Hydrophobic Surfaces

The properties of water on both molecular and macroscopic surfaces critically influence a wide range of physical behaviors, with applications spanning from membrane science to catalysis to protein engineering. Yet, our current understanding of water interfacing molecular and material surfaces is incomplete, in part because measurement of water structure and molecular-scale properties challenges even the most advanced experimental characterization techniques and computational approaches. This review highlights progress in the ongoing development of tools working to answer fundamental questions on the principles that govern the interactions between water and surfaces. One outstanding and critical question is what universal molecular signatures capture the hydrophobicity of different surfaces in an operationally meaningful way, since traditional macroscopic hydrophobicity measures like contact angles fail to capture even basic properties of molecular or extended surfaces with any heterogeneity at the nanometer length scale. Resolving this grand challenge will require close interactions between state-of-the-art experiments, simulations, and theory, spanning research groups and using agreed-upon model systems, to synthesize an integrated knowledge of solvation water structure, dynamics, and thermodynamics.

Forces between zinc sulphide surfaces; amplification of the hydrophobic attraction by surface charge

Smooth Zinc Sulphide (ZnS) surfaces were prepared by magnetron sputtering and the interaction forces were measured between them as a function of pH. At the isoelectric point (iep) of pH 7.1 the attractive force was well described by the van der Waals interaction calculated using Lifshitz theory for a layered system. Away from the iep, the forces were fitted using DLVO theory extended to account for surface roughness. At pH 9.8 the surfaces acquire a negative charge and an electrostatic repulsion is evident. Below the iep the surfaces acquire a positive charge leading to electrostatic repulsion. The forces in the range 3.8 < pH < 4.8 show an additional attraction on approach and much greater adhesion than at other pH values. This is attributed to the hydrophobic attraction being amplified by a small degree of charge on the surface as has previously been reported for adhesion measurements. The range of the measured forces is attributed to the long-range orientational order of water (>5 nm).

Decoding signatures of structure, bulk thermodynamics, and solvation in three-body angle distributions of rigid water models

A tetrahedral structure resulting from hydrogen bonding is a hallmark of liquid water and plays a significant role in determining its unique thermophysical properties. This water feature has helped understand anomalous properties and physically interpret and model hydrophobic solvation thermodynamics. Tetrahedrality is well described by the geometric relationship of any central water molecule with two of its nearest neighbors in the first coordination shell, as defined by the corresponding "three-body" angle. While order parameters and even full water models have been developed using specific or average features of the three-body angle distribution, here we examine the distribution holistically, tracking its response to changes in temperature, density, and the presence of model solutes. Surprisingly, we find that the three-body distribution responds by varying primarily along a single degree of freedom, suggesting a remarkably simplified view of water structure. We characterize three-body angle distributions across temperature and density space and identify principal components of the variations with state conditions. We show that these principal components embed physical significance and trace out transitions between tetrahedral and simple-fluid-like behavior. Moreover, we find that the ways three-body angles vary within the hydration shells of model colloids of different types and sizes are nearly identical to the variations seen in bulk water across density and temperature. Importantly, through the principal directions of these variations, we find that perturbations to the hydration-water distributions well predict the thermodynamics associated with colloid solvation, in particular, the relative entropy of this process that captures indirect, solvent-mediated contributions to the hydration free energy.

Molecular dynamics simulation of proteins under high pressure: Structure, function and thermodynamics

BACKGROUND

Molecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.

SCOPE OF REVIEW

First, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.

MAJOR CONCLUSIONS

Recent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.

GENERAL SIGNIFICANCE

MD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.

Pressure response of the THz spectrum of bulk liquid water revealed by intermolecular instantaneous normal mode analysis

The radial distribution functions of liquid water are known to change significantly their shape upon hydrostatic compression from ambient conditions deep into the kbar pressure regime. It has been shown that despite their eye-catching changes, the fundamental locally tetrahedral fourfold H-bonding pattern that characterizes ambient water is preserved up to about 10 kbar (1 GPa), which is the stability limit of liquid water at 300 K. The observed increase in coordination number comes from pushing water molecules into the first coordination sphere without establishing an H-bond, resulting in roughly two such additional interstitial molecules at 10 kbar. THz spectroscopy has been firmly established as a powerful experimental technique to analyze H-bonding in aqueous solutions given that it directly probes the far-infrared lineshape and thus the prominent H-bond network mode around 180 cm^-1. We, therefore, set out to assess pressure effects on the THz response of liquid water at 10 kbar in comparison to the 1 bar (0.1 MPa) reference, both at 300 K, with the aim to trace back the related lineshape changes to the structural level. To this end, we employ the instantaneous normal mode approximation to rigorously separate the H-bonding peak from the large background arising from the pronounced librational tail. By exactly decomposing the total molecular dynamics into hindered translations, hindered rotations, and intramolecular vibrations, we find that the H-bonding peak arises from translation-translation and translation-rotation correlations, which are successively decomposed down to the level of distinct local H-bond environments. Our utmost detailed analysis based on molecular pair classifications unveils that H-bonded double-donor water pairs contribute most to the THz response around 180 cm^-1, whereas interstitial waters are negligible. Moreover, short double-donor H-bonds have their peak maximum significantly shifted toward higher frequencies with respect to such long H-bonds. In conjunction with an increasing relative population of these short H-bonds versus the long ones (while the population of other water pair classes is essentially pressure insensitive), this explains not only the blue-shift of the H-bonding peak by about 20-30 cm^-1 in total from 1 bar to 10 kbar but also the filling of the shallow local minimum of the THz lineshape located in between the network peak and the red-wing of the librational band at 1 bar. Based on the changing populations as a function of pressure, we are also able to roughly estimate the pressure-dependence of the H-bond network mode and find that its pressure response and thus the blue-shifting are most pronounced at low kbar pressures.