Effect of Conformational Equilibrium on Solvation Properties of 1,2-DCE in Water: A Solvation Thermodynamics and 3D-RISM Study

The contributions of the enthalpy and entropy of solvation for the study of chemical and biological systems are important in the prediction, interpretation, and manipulation of these processes. The relation between solvation Gibbs energies, enthalpies, and entropies of solvation, and their rigorous relation with the conformational equilibrium, are derived for the first time and applied with a computational method, in accordance with the Solvation Thermodynamics previous results, to 1,2-dichloroethane solvation in water. The rigid conformer calculations in solution were performed by using PC+/3D-RISM approach, with the conformational averaged results for enthalpy and solvation Gibbs energy reproducing the experimental results quite successfully. A qualitative agreement in the entropy of solvation predictions was also observed.

Correlating Antiagglomerant Performance with Gas Hydrate Cohesion

Although inhibiting hydrate formation in hydrocarbon-water systems is paramount in preventing pipe blockage in hydrocarbon transport systems, the molecular mechanisms responsible for antiagglomerant (AA) performance are not completely understood. To better understand why macroscopic performance is affected by apparently small changes in the AA molecular structure, we perform molecular dynamics simulations. We quantify the cohesion energy between two gas hydrate nanoparticles dispersed in liquid hydrocarbons in the presence of different AAs, and we achieve excellent agreement against experimental data obtained at high pressure using the micromechanical force apparatus. This suggests that the proposed simulation approach could provide a screening method for predicting, in silico, the performance of new molecules designed to manage hydrates in flow assurance. Our results suggest that entropy and free energy of solvation of AAs, combined in some cases with the molecular orientation at hydrate-oil interfaces, are descriptors that could be used to predict performance, should the results presented here be reproduced for other systems as well. These insights could help speed up the design of new AAs and guide future experiments.

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 Shape and the Hydrophobic Effect

This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.

Temperature-dependent solvation modulates the dimensions of disordered proteins

For disordered proteins, the dimensions of the chain are an important property that is sensitive to environmental conditions. We have used single-molecule Förster resonance energy transfer to probe the temperature-induced chain collapse of five unfolded or intrinsically disordered proteins. Because this behavior is sensitive to the details of intrachain and chain-solvent interactions, the collapse allows us to probe the physical interactions governing the dimensions of disordered proteins. We find that each of the proteins undergoes a collapse with increasing temperature, with the most hydrophobic one, λ-repressor, undergoing a reexpansion at the highest temperatures. Although such a collapse might be expected due to the temperature dependence of the classical "hydrophobic effect," remarkably we find that the largest collapse occurs for the most hydrophilic, charged sequences. Using a combination of theory and simulation, we show that this result can be rationalized in terms of the temperature-dependent solvation free energies of the constituent amino acids, with the solvation properties of the most hydrophilic residues playing a large part in determining the collapse.

Sensitivity of local hydration behaviour and conformational preferences of peptides to choice of water model

Secondary structural preferences of the beta-hairpin of the 2GB1 protein in the folded and unfolded ensembles are shown to be sensitive to the choice of water model.

Different aggregation dynamics of benzene–water mixtures

The differences between the molecular aggregations in benzene–water mixtures are identified using all-atom molecular dynamics simulations.

Soft collective fluctuations governing hydrophobic association

The interaction between two associating hydrophobic particles has traditionally been explained in terms of the release of entropically frustrated hydration shell water molecules. However, this picture cannot account for the kinetics of hydrophobic association and is therefore not capable of providing a microscopic description of the hydrophobic interaction (HI). Here, Monte Carlo simulations of a pair of molecular-scale apolar solutes in aqueous solution reveal the critical role of collective fluctuations in the hydrogen bond (HB) network for the microscopic picture of the HI. The main contribution to the HI is the relaxation of solute-water translational correlations. The existence of a heat capacity maximum at the desolvation barrier is shown to arise from softening of non-HB water fluctuations and the relaxation of many-body correlations in the labile HB network. The microscopic event governing the kinetics of hydrophobic association has turned out to be a relatively large critical collective fluctuation in hydration water displacing a substantial fraction of HB clusters from the inner to the outer region of the first hydration shell.

Physical origin underlying the entropy loss upon hydrophobic hydration

The hydrophobic effect (HE) is commonly associated with the demixing of oil and water at ambient conditions and plays the leading role in determining the structure and stability of biomolecular assembly in aqueous solutions. On the molecular scale HE has an entropic origin. It is believed that hydrophobic particles induce order in the surrounding water by reducing the volume of configuration space available for hydrogen bonding. Here we show with computer simulation results that this traditional picture, based on average structural features of hydration water, configurational properties of single water molecules, and up to pairwise correlations, is not correct. Analyzing collective fluctuations in water clusters we are able to provide a fundamentally new picture of HE based on pronounced many-body correlations affecting the switching of hydrogen bonds (HBs) between molecules. These correlations emerge as a nonlocal compensation of reduced fluctuations of local electrostatic fields in the presence of an apolar solute. We propose an alternative view which may also be formulated as a maximization principle: The electrostatic noise acting on water molecules is maximized under the constraint that each water molecule on average maintains as many HBs as possible. In the presence of the solute the maximized electrostatic noise is a result of nonlocal fluctuations in the labile HB network giving rise to strong correlations among at least up to four water molecules.

Water and other tetrahedral liquids: order, anomalies and solvation

In order to understand the common features of tetrahedral liquids with water-like anomalies, the relationship between local order and anomalies has been studied using molecular dynamics simulations for three categories of such liquids: (a) atomistic rigid-body models for water (TIP4P, TIP4P/2005, mTIP3P, SPC/E), (b) ionic melts, BeF(2) (TRIM model) and SiO(2) (BKS potential) and (c) Stillinger-Weber liquids parametrized to model water (mW) and silicon. Rigid-body, atomistic models for water and the Stillinger-Weber liquids show a strong correlation between tetrahedral and pair correlation order and the temperature for the onset of the density anomaly is close to the melting temperature. In contrast, the ionic melts show weaker and more variable degrees of correlation between tetrahedral and pair correlation metrics, and the onset temperature for the density anomaly is more than twice the melting temperature. In the case of water, the relationship between water-like anomalies and solvation is studied by examining the hydration of spherical solutes (Na(+), Cl(-), Ar) in water models with different temperature regimes of anomalies (SPC/E, TIP4P and mTIP3P). For both ionic and nonpolar solutes, the local structure and energy of water molecules is essentially the same as in bulk water beyond the second-neighbour shell. The local order and binding energy of water molecules are not perturbed by the presence of a hydrophobic solute. In the case of ionic solutes, the perturbation is largely localized within the first hydration shell. The binding energies for the ions are strongly dependent on the water models and clearly indicate that the geometry of the partial charge distributions, and the associated multipole moments, play an important role. However the anomalous behaviour of the water network has been found to be unimportant for polar solvation.

Predicting cavity formation free energy: how far is the Gaussian approximation valid?

We examine the range of validity of the Gaussian model for various water-like liquids whose intermolecular potentials differ from SPC/E water, to provide insight into the temperature dependence of the hydrophobic effect for small hard sphere solutes. We find that low compressibility liquids that have more close-packed network structures show much larger deviations from Gaussian fluctuations for low or zero occupancies relative to more compressible fluids with more open networks. Water appears to be a unique molecular fluid in possessing equilibrium density fluctuations that are faithfully described by the Gaussian theory. We ascribe this success to the fact, shown here, that the orientational correlations near a small hard sphere solute involve remarkably little reorganization from the bulk, which is a consequence of water's low solvent reorganization enthalpy and entropy.

Comparison of Tetrahedral Order, Liquid State Anomalies, and Hydration Behavior of mTIP3P and TIP4P Water Models

The relationship between local tetrahedral order, tagged particle potential energy, and coordination number is studied for mTIP3P and TIP4P models of water in the bulk as well as in the neighborhood of a small peptide. The tendency of water molecules with different binding or tagged particle potential energies to occupy environments with different degrees of disorder can be effectively illustrated by constructing tetrahedral order distributions and corresponding entropy metrics conditional on restricted ranges of local binding energy. At the state point corresponding to the onset of the density anomaly, the correlation between tetrahedral entropy versus tagged potential energy is strong and virtually identical for mTIP3P and TIP4P. In TIP4P, this correlation is retained up to temperatures as high as 300 K, while it is lost by 250 K in mTIP3P. In the 250-300 K regime that is important for biomolecular simulations, mTIP3P behaves essentially as a simple liquid while TIP4P shows the density and related anomalies characteristic of water. We also study the number of water molecules, the tetrahedral order, and the tagged molecule potential energies for water molecules as a function of the distance from the peptide for the 16-residue β-hairpin fragment of 2GB1 in mTIP3P and TIP4P solvents. The hydration shell coordination profiles (n(r)) of the number of water molecules are almost identical in the two solvents, but the radial variation in the local energies and local order show significant differences. The residue-wise variation in the tagged potential energy of water molecules within the first hydration shell is qualitatively similar in the two models. A comparison of the tetrahedral order distributions of water molecules lying at different distances from the biomolecular solute shows that the perturbation in the local tetrahedral order distributions of the bulk solvent due to the presence of the solute is marginal. Thus, in the 250-300 K regime, the mTIP3P and TIP4P water models show qualitatively different behavior in terms of the relationship between tetrahedral order and local energy, but as solvents in the neighborhood of a biomolecular solute, the differences between the two models are only quantitative and not qualitative.