Thermodynamic and Structural Evidence for Reduced Hydrogen Bonding among Water Molecules near Small Hydrophobic Solutes

Microscopic insights on clathrate hydrate growth from non-equilibrium molecular dynamics simulations

Clathrate hydrates form and grow at interfaces. Understanding the relevant molecular processes is crucial for developing hydrate-based technologies. Many computational studies focus on hydrate growth within the aqueous phase using the 'direct coexistence method', which is limited in its ability to investigate hydrate film growth at hydrocarbon-water interfaces. To overcome this shortcoming, a new simulation setup is presented here, which allows us to study the growth of a methane hydrate nucleus in a system where oil-water, hydrate-water, and hydrate-oil interfaces are all simultaneously present, thereby mimicking experimental setups. Using this setup, hydrate growth is studied here under the influence of two additives, a polyvinylcaprolactam oligomer and sodium dodecyl sulfate, at varying concentrations. Our results confirm that hydrate films grow along the oil-water interface, in general agreement with visual experimental observations; growth, albeit slower, also occurs at the hydrate-water interface, the interface most often interrogated via simulations. The results obtained demonstrate that the additives present within curved interfaces control the solubility of methane in the aqueous phase, which correlates with hydrate growth rate. Building on our simulation insights, we suggest that by combining data for the potential of mean force profile for methane transport across the oil-water interface and for the average free energy required to perturb a flat interface, it is possible to predict the performance of additives used to control hydrate growth. These insights could be helpful to achieve optimal methane storage in hydrates, one of many applications which are attracting significant fundamental and applied interests.

New insights on the catalytic center of proteins from peptidylprolyl isomerase group based on the FOD-M model

Generating the structure of the hydrophobic core is based on the orientation of hydrophobic residues towards the central part of the protein molecule with the simultaneous exposure of polar residues. Such a course of the protein folding process takes place with the active participation of the polar water environment. While the self-assembly process leading to the formation of micelles concerns freely moving bi-polar molecules, bipolar amino acids in polypeptide chain have limited mobility due to the covalent bonds. Therefore, proteins form a more or less perfect micelle-like structure. The criterion is the hydrophobicity distribution, which to a greater or lesser extent reproduces the distribution expressed by the 3D Gaussian function on the protein body. The vast majority of proteins must ensure solubility, so a certain part of it-as it is expected-should reproduce the structuring of micelles. The biological activity of proteins is encoded in the part that does not reproduce the micelle-like system. The location and quantitative assessment of the contribution of orderliness to disorder is of critical importance for the determination of biological activity. The form of maladjustment to the 3D Gauss function may be varied-hence the obtained high diversity of specific interactions with strictly defined molecules: ligands or substrates. The correctness of this interpretation was verified on the basis of the group of enzymes Peptidylprolyl isomerase-E.C.5.2.1.8. In proteins representing this class of enzymes, zones responsible for solubility-micelle-like hydrophobicity system-the location and specificity of the incompatible part in which the specific activity of the enzyme is located and coded were identified. The present study showed that the enzymes of the discussed group show two different schemes of the structure of catalytic center (taking into account the status as defined by the fuzzy oil drop model).

Crustwater: Modeling Hydrophobic Solvation

We describe Crustwater, a statistical mechanical model of nonpolar solvation in water. It treats bulk water using the Cage Water model and introduces a crust, i.e., a solvation shell of coordinated partially structured waters. Crustwater is analytical and fast to compute. We compute here solvation vs temperature over the liquid range, and vs pressure and solute size. Its thermal predictions are as accurate as much more costly explicit models such as TIP4P/2005. This modeling gives new insights into the hydrophobic effect: (1) that oil-water insolubility in cold water is due to solute-water (SW) translational entropy and not water-water (WW) orientations, even while hot water is dominated by WW cage breaking, and (2) that a size transition at the Angstrom scale, not the nanometer scale, takes place as previously predicted.

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.

Diffusivity and hydrophobic hydration of hydrocarbons in supercritical CO2 and aqueous brine

CO₂ injection (EOR and sequestration technique) creates the amalgamation of hydrocarbons, CO₂, and aqueous brine in the subsurface. In this study, molecular dynamics (MD) simulations were used to investigate the diffusivity of hydrocarbon molecules in a realistic scenario of supercritical CO₂ (SC-CO₂) injection in the subsurface over a wide range of pressures (50 < P < 300 bar) and aqueous brine concentrations (0, 2, and 5% brine). To overcome existing challenges in traditional diffusivity calculation approaches, we took advantage of fundamental molecular-based methods, along with further verification of results by previously published experimental data. In this regard, computational methods and MD simulations were employed to compute diffusion coefficients of hydrocarbons (benzene and pentane). It was found that the presence of water and salt affects the thermodynamic properties of molecules where the intermolecular interactions caused the hydrophobic hydration of hydrocarbons coupled with ionic hydration due to hydrogen bond and ion-dipole interactions. Based on these results, it is demonstrated that the formation of water clusters in the SC-CO₂ solvent is a major contributor to the diffusion of hydrophobic molecules. The outcome at different pressure conditions showed that hydrocarbons always would diffuse less in the presence of water. The slopes of linearly fitted MSD of benzene and pentane infinitely diluted in SC-CO₂ is around 13 to 20 times larger than the slope with water molecules (4 wt%). When pressure increases (100–300 bar), the diffusion coefficients (D) of benzene and pentane decreases (around 1.2 × 10⁻⁹ to 0.4 × 10⁻⁹ and 2 × 10⁻⁹ to 1 × 10⁻⁹ m² s⁻¹, respectively). Furthermore, brine concentration generally plays a negative role in reducing the diffusivity of hydrocarbons due to the formation of water clusters as a result of hydrophobic and ionic hydration. Under the SC-CO₂ rich (injection) system in the shale reservoir, the diffusion of hydrocarbon is correlated to the efficiency of hydrocarbon flow/recovery. Ultimately, this study will guide us to better understand the phenomena that would occur in nanopores of shale that undergo EOR or are becoming a target of CO₂ sequestration.

CO₂ injection (EOR and sequestration technique) creates the amalgamation of hydrocarbons, CO₂, and aqueous brine in the subsurface.

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.

Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure

Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.

Contrasting the hydration thermodynamics of methane and methanol

The hydration thermodynamics of methane and methanol depend on the cavity creation work and energy of van der Waals and H-bonding attractions.

Interfacial structural crossover and hydration thermodynamics of charged C60 in water

Classical molecular dynamics simulations of the hydration thermodynamics, structure, and dynamics of water in hydration shells of charged buckminsterfullerenes are presented in this study. Charging of fullerenes leads to a structural transition in the hydration shell, accompanied by creation of a significant population of dangling O-H bonds pointing toward the solute. In contrast to the well accepted structure-function paradigm, this interfacial structural transition causes nearly no effect on either the dynamics of hydration water or on the solvation thermodynamics. Linear response to the solute charge is maintained despite significant structural changes in the hydration shell, and solvation thermodynamic potentials are nearly insensitive to the altering structure. Only solvation heat capacities, which are higher thermodynamic derivatives of the solvation free energy, indicate some sensitivity to the local hydration structure. We have separated the solvation thermodynamic potentials into direct solute-solvent interactions and restructuring of the hydration shell and analyzed the relative contributions of electrostatic and nonpolar interactions to the solvation thermodynamics.

Hydrophobic Hydration and the Effect of NaCl Salt in the Adsorption of Hydrocarbons and Surfactants on Clathrate Hydrates

Adsorption of functional molecules on the surface of hydrates is key in the understanding of hydrate inhibitors. We investigate the adsorption of a hydrocarbon chain, nonionic and ionic surfactants, and ions at the hydrate-aqueous interface. Our results suggest a strong connection between the water ordering around solutes in bulk and the affinity for the hydrates surface. We distinguish two types of water ordering around solutes: (i) hydrophobic hydration where water molecules form a hydrogen bond network similar to clathrate hydrates, and (ii) ionic hydration where water molecules align according to the polarity of an ionic group. The nonionic surfactant and the hydrocarbon chain induce hydrophobic hydration and are favorably adsorbed on the hydrate surface. Adsorption of ions and the ionic headgroups on the hydrate surface is not favorable because ionic hydration and the hydrogen bond structure of hydrates are incompatible. The nonionic surfactant is adsorbed by the headgroup and tail while adsorption of the ionic surfactants is not favorable through the head. Water ordering is analyzed using the hydrogen bond and tetrahedral density profiles as a function of the distance to the chemical groups. The adsorption of solutes is studied through the free energy profiles as a function of the distance to the hydrate surface. Salt lowers the melting temperature of hydrates, disrupts hydrophobic hydration, reduces the solubility of solutes in the aqueous solution, and increases the propensity of solutes to be adsorbed on hydrate surfaces. Our studies are performed by the unbiased and steered molecular dynamics simulations. The results are in line with experiments on the effect of salt and alkanes in hydrate antiagglomeration.

Water structure around hydrophobic amino acid side chain analogs using different water models

The water structure around hydrophobic groups governs various biochemical processes. There is an ongoing debate on whether water molecules near hydrophobic groups are more ordered with greater participation in water-water hydrogen bonding with respect to water in the pure bulk state. The water structure around six different hydrophobic amino acid side chain analog molecules has been studied in pure water using molecular dynamics simulations. The analysis of water tetrahedral order parameter and the number of hydrogen bonds formed by the individual water molecules in the first hydration shell of the hydrophobic analogs provide evidence that both ordering and hydrogen bonds involving water molecules are to some extent reduced in the hydrophobic hydration shell. It is revealed that the water tetrahedrality in the outer part of the first hydrophobic hydration shell is equivalent to bulk water for all the water models except for the TIP4P-2005 model which shows marginally higher tetrahedrality. However, irrespective of the model employed, water tetrahedrality has always been found to be reduced in the inner part of the first hydration shell, which eventually makes the overall water tetrahedrality in the first hydrophobic hydration shell marginally lower than that observed for pure bulk water. Importantly, it is noticed that the decrease in water structuring exhibits solute size dependencies. Around a small solute like methane, the water tetrahedral ordering or hydrogen bonding propensity is quite similar to that of the bulk state. The effect, reduction in water structuring, is however more pronounced for relatively larger solutes.