Culation of partial molar volume and its components for molecular dynamics models of dilute solutions

This paper is a review of our recent computational studies of volumetric characteristics using computer models of dilute solutions. Partial molar volume (PMV) and its components are calculated for simple and complex molecules in water (methane, noble gases, surfactants, polypeptides). Advantages and disadvantages of various computational methods are discussed. It is proposed to use the Voronoi-Delaunay technique to determine the reasonable boundary between a solute molecule and solvent molecules and to identify the PMV components related to the molecule, the boundary layer, and the solvent. It is noted that the observed increase in PMV with temperature for large molecules is due to an increase in the volume of voids in the boundary layer, i.e., due to the “thermal volume.” In this case, the solvent gives a negative contribution to the PMV. In contrast, for simple molecules (methane), the contribution from the solvent is positive and is the main factor in the increase in the PMV, which is associated with a specific change in water structure around a spherical hydrophobic particle outside the boundary layer. For surfactant molecules, the contribution from the solvent changes sign (from negative to positive) with increasing temperature.     

Volumetric and compressibility studies for (L-arginine + D-maltose monohydrate + water) system in the temperature range of (298.15 to 308.15) K

The density and speed of sound of L-arginine (0.025–0.2 mol kg^?1) in aqueous + D-maltose (0–6 mass% of maltose in water) were obtained at temperatures of (298.15, 303.15 and 308.15) K. The apparent molar volume, limiting apparent molar volume, transfer volume, as well as apparent molar compressibility, limiting apparent molar compressibility, transfer compressibility, pair and triple interaction coefficients, partial molar expansibilities, coefficient of thermal expansion and also the hydration number, were calculated using the experimental density and speed of sound values. The results have been discussed in terms of solute–solute and solute–solvent interactions in these systems. Solute–solvent (hydrophilic–ionic group and hydrophilic–hydrophilic group) interactions were found to be dominating over solute–solute (hydrophobic–hydrophilic group) interactions in the solution, which increases with increase in maltose concentration.     

Theoretical study of the cosolvent effect on the partial molar volume change of staphylococcal nuclease associated with pressure denaturation

We explain the molecular mechanism of the effect of urea and glycerol cosolvents on the partial molar volume (PMV) change associated with the pressure denaturation of staphylococcal nuclease (SNase) protein recently observed in experiments. Native and denatured conformations of SNase are produced by using molecular dynamics simulations in water, and the PMV is obtained from the integral equation theory of molecular liquids called 3D-RISM, which is based on statistical mechanics. The PMV of the native SNase in water predicted by 3D-RISM theory is in good agreement with experiment. The PMV changes associated with pressure denaturation in water and in water-urea and water-glycerol mixtures are qualitatively reproduced. By analyzing the results obtained, we found two interesting cosolvent effects on the PMV: (1) both urea and glycerol cosolvents increase the PMVs of both native and denatured SNase compared to those in water and (2) both urea and glycerol cosolvents increase the PMV of denatured SNase more than that of native SNase. We also showed that these two observations can be explained in terms of the thermal volume, which is related to the packing effect of solvent molecules.     

Different Views on the Stability of Protein Conformations and Hydrophobic Effects

Apparently contradictory statements about the thermodynamics of aqueous protein solutions and of hydrophobic effect are quoted and discussed. Some credibility is found in the divergent points of view and it is pointed out that they focus attention on different aspects of the complicated conditions in aqueous solutions, some of which are more important than others for the stability of protein conformations.The importance of characteristics of solvent water is emphasized, in particular (1) the strong mutual cohesion of water molecules, and (2) structural changes of water induced by (nonpolar) solute molecules. It is stressed that consideration of only one of these effects and an inexpedient choice of standard states are origins of confusion in the literature about aqueous systems. A simple approach to hydrophobic effects considering both of the above mentioned effects, is proposed.     

Solvent reorganization contributions in solute transfer thermodynamics: inferences from the solvent equation of state

Solvation enthalpies of simple solutes contain contributions from (1) solute-solvent interactions and (2) solute-induced modifications of solvent-solvent interactions (solvent reorganization). It has recently been suggested in the literature that these contributions can, under certain conditions, be estimated with additional experimental data on thermodynamic response functions of the pure solvent (coefficient of thermal expansion, isothermal compressibility) and the solute solvation volume. We analyze and discuss these conditions based on computer simulations of a series of polar and nonpolar solutes in a polar and nonpolar liquid solvent.     

Solvation of hydrophobes in water and simple liquids

The solvation of nonpolar molecules in water and that in simple liquids are compared and contrasted. First, solvation thermodynamics is reviewed in a way that focuses on how the enthalpy and entropy of solvation depend on the choice of microscopic volume change v in the solvation process--including special choices v being zero (fixed-volume condition) and v being the partial molecular volume of a solute molecule (fixed-pressure condition)--and how the solvation quantities are related with temperature derivatives of the solvation free energy. Second, the solvation free energy and the solvation enthalpy of a Lennard-Jones (LJ) atom in model water are calculated in the parameter space representing the solute size and the strength of the solute-solvent interaction, and the results are compared with those for an LJ atom in the LJ solvent. The solvation diagrams showing domains of different types of solvation in the parameter space are obtained both for the constant-volume condition and for the constant-pressure condition. Similarities between water and the simple liquid are found when the constant-volume solvation is considered while a significant difference manifests itself in the fixed-pressure solvation. The domain of solvation of hydrophobic character in the parameter space is large in the constant-volume solvation both for water and for the simple liquid. When switched to the constant-pressure condition accompanying a microscopic volume change, the hydrophobic domain remains large in water but it becomes significantly small in the simple liquid. The contrasting results are due to the smallness of the thermal pressure coefficient of water at low temperatures.     

Water's hydrogen bonds in the hydrophobic effect: a simple model

We propose a simple analytical model to account for water's hydrogen bonds in the hydrophobic effect. It is based on computing a mean-field partition function for a water molecule in the first solvation shell around a solute molecule. The model treats the orientational restrictions from hydrogen bonding, and utilizes quantities that can be obtained from bulk water simulations. We illustrate the principles in a 2-dimensional Mercedes-Benz-like model. Our model gives good predictions for the heat capacity of hydrophobic solvation, reproduces the solvation energies and entropies at different temperatures with only one fitting parameter, and accounts for the solute size dependence of the hydrophobic effect. Our model supports the view that water's hydrogen bonding propensity determines the temperature dependence of the hydrophobic effect. It explains the puzzling experimental observation that dissolving a nonpolar solute in hot water has positive entropy.     

Lysozyme pattern formation in evaporating drops

Liquid droplets containing suspended particles deposited on a solid, flat surface generally form ringlike structures due to the redistribution of solute during evaporation (the "coffee ring effect"). The forms of the deposited patterns depend on interactions between solute(s), solvent, and substrate. In this study, deposition patterns from droplets of a simplified model biological fluid (DI water + lysozyme) are examined by scanning probe and optical microscopy. The overall lysozyme residue morphology is complex (with both a perimeter "rim" and undulating interior) but varies little with concentration. However, the final packing of lysozyme molecules is strongly dependent on initial concentration.     

Structural characteristics of a 0.23 mole fraction aqueous solution of tetrahydrofuran at 20 degrees C

Hydrogen/deuterium isotopic substitution neutron diffraction techniques were used to measure the structural correlation functions in a 0.23 mole fraction solution of tetrahydrofuran in water at room temperature. Empirical potential structure refinement (EPSR) was used to build a three-dimensional model of the liquid structure that is consistent with the experimental data. Detailed analysis shows a preference for nonpolar interactions between the cyclic ether molecules plus polar interactions between the ether and solvent water and hydrophobic hydration of the nonpolar regions of the solute. The increase in the number of hydrogen-bond-acceptor sites relative to the number of hydrogen-bond-donor sites in this system, compared to the balanced situation that would be found in pure water, has a marked compressive effect on the structure of the solvent. Despite the small size of the solvent water molecules, the 0.23 mole fraction aqueous solution is still found to contain small voids akin to those in pure liquid tetrahydrofuran. In contrast to the positive surface charge of the voids in the pure system, the average void in this aqueous solution is found to have a net negative charge. This is due to contributions from the water oxygen atoms that are negatively polarized by their intramolecular bonding.     

Molecular origin of anticooperativity in hydrophobic association

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.     

Effect of pressure on rotational isomers in solution: infrared and Raman spectroscopy

The effect of pressure on the conformational equilibria between rotational isomers in solution has been studied by the infrared and Raman spectroscopy. The experimental volume changes for the transformation are explained theoretically by the three effects contributing to these volume changes, namely the overlap of the atoms with the internal rotation angle, the change of local packing of the solute and solvent, and the change of the electrostatic interaction with the solvent. The purpose of this paper is to review what is known about the volume changes for small molecules with internal rotation angles.     

A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. II. Noncollisional mechanism for the relaxation of a polar solute in supercritical water

Mixed quantum-classical molecular dynamics method has been applied to vibrational relaxation of a hydrophilic model NO in supercritical water at various densities along an isotherm above the critical temperature. The relaxation rate was determined based on Fermi's golden rule at each state point and showed an inverse S-shaped curve as a function of bulk density. The hydration number was also calculated as a function of bulk density based on the calculated radial distribution function, which showed a good correlation with the relaxation rate. Change of the survival probability of the solute vibrational state was analyzed as a function of time together with the trajectory of the solvent water and the interaction with it. We will show that the solvent molecule resides near the solute molecule for a while and the solvent contributes to the relaxation by the random-noiselike Coulombic interaction only when it stays near the solute. After the solvent leaves the solute, it shows no contribution to the relaxation. The relaxation mechanism for this system is significantly different from the collisional one found for a nonpolar solute in nonpolar solvent in Paper I. Then, the relaxation rate is determined, on average, by the hydration number or local density of the solvent. Thus, the density dependence of the relaxation rate for the polar solute in supercritical water is apparently similar to that found for the nonpolar solute in nonpolar solvent, although the molecular process is quite different from each other.     

Thermal signature of hydrophobic hydration dynamics

Hydrophobic hydration, the perturbation of the aqueous solvent near an apolar solute or interface, is a fundamental ingredient in many chemical and biological processes. Both bulk water and aqueous solutions of apolar solutes behave anomalously at low temperatures for reasons that are not fully understood. Here, we use (2)H NMR relaxation to characterize the rotational dynamics in hydrophobic hydration shells over a wide temperature range, extending down to 243 K. We examine four partly hydrophobic solutes: the peptides N-acetyl-glycine-N'-methylamide and N-acetyl-leucine-N'-methylamide, and the osmolytes trimethylamine N-oxide and tetramethylurea. For all four solutes, we find that water rotates with lower activation energy in the hydration shell than in bulk water below 255 +/- 2 K. At still lower temperatures, water rotation is predicted to be faster in the shell than in bulk. We rationalize this behavior in terms of the geometric constraints imposed by the solute. These findings reverse the classical "iceberg" view of hydrophobic hydration by indicating that hydrophobic hydration water is less ice-like than bulk water. Our results also challenge the "structural temperature" concept. The two investigated osmolytes have opposite effects on protein stability but have virtually the same effect on water dynamics, suggesting that they do not act indirectly via solvent perturbations. The NMR-derived picture of hydrophobic hydration dynamics differs substantially from views emerging from recent quasielastic neutron scattering and pump-probe infrared spectroscopy studies of the same solutes. We discuss the possible reasons for these discrepancies.     

Universal model based on the mobile order and disorder theory for predicting lipophilicity and partition coefficients in all mutually immiscible two-phase liquid systems

The quantitative thermodynamic development of the mobile order and disorder theory in H-bonded liquids has been extended in order to predict partition coefficients. The model enables "a priori" estimation of the partition coefficient (log P) of neutral solutes, not only in the conventional 1-octanol/water reference but also in all mutually saturated two-phase systems made up of largely immiscible solvents. The model is obtained from the thermodynamic treatment of the various physicochemical free energy processes encoded in the overall distribution process and accordingly provides a useful tool for better understanding both the origin and the factors, such as the solute molar volume, that determine the partition coefficient of nonelectrolytes in a given system. From the comparison of the relative magnitude of the processes contributing to the log P value, a lot of information can also be gained regarding the variation in log P of the same substance partitioned between different solvent systems. As a demonstration, the model has been successfully applied to predict the log P of a great number of chemicals of varying structure, size, and chemical nature partitioned in a large set of essentially immiscible solvent pairs, differing either by their nonpolar or by their polar phase. In the systems involving water as the polar phase, the hydrophobic effect is always the driving force that governs the distribution process irrespective of the interacting or noninteracting nature of the substances studied. In the other two-phase systems, the partitioning of complexing solutes in particular appears to be ruled rather by their hydrogen-bonding capabilities than by their hydrophobicities.     

Pressure and temperature dependence of hydrophobic hydration: volumetric, compressibility, and thermodynamic signatures

The combined effect of pressure and temperature on hydrophobic hydration of a nonpolar methanelike solute is investigated by extensive simulations in the TIP4P model of water. Using test-particle insertion techniques, free energies of hydration under a range of pressures from 1 to 3000 atm are computed at eight temperatures ranging from 278.15 to 368.15 K. Corresponding enthalpy, entropy, and heat capacity accompanying the hydration process are estimated from the temperature dependence of the free energies. Partial molar and excess volumes calculated using pressure derivatives of the simulated free energies are consistent with those determined by direct volume simulations; but direct volume determination offers more reliable estimates for compressibility. At 298.15 K, partial molar and excess isothermal compressibilities of methane are negative at 1 atm. Partial molar and excess adiabatic (isentropic) compressibilities are estimated to be also negative under the same conditions. But partial molar and excess isothermal compressibilities are positive at high pressures, with a crossover from negative to positive compressibility at approximately 100-1000 atm. This trend is consistent with experiments on aliphatic amino acids and pressure-unfolded states of proteins. For the range of pressures simulated, hydration heat capacity exhibits little pressure dependence, also in apparent agreement with experiment. When pressure is raised at constant room temperature, hydration free energy increases while its entropic component remains essentially constant. Thus, the increasing unfavorability of hydration under raised pressure is seen as largely an enthalpic effect. Ramifications of the findings of the authors for biopolymer conformational transitions are discussed.     

Volumetric,Acoustic and Transport Properties of Metformin Hydrochloride Drug in Aqueous <Emphasis Type="SmallCaps">d</Emphasis>-Glucose Solutions at <Emphasis Type="Italic">T</Emphasis> = (298.15–318.15) K

Apparent molar volumes, apparent molar adiabatic compressibilities and viscosity B-coefficients for metformin hydrochloride in aqueous d-glucose solutions were determined from solution densities, sound velocities and viscosities measured at T = (298.15–318.15) K and at pressure p = 101 kPa as a function of the metformin hydrochloride concentrations. The standard partial molar volumes (\( \phi_{V}^{0} \)) and slopes (\( S_{V}^{*} \)) obtained from the Masson equation were interpreted in terms of solute–solvent and solute–solute interactions, respectively. Solution viscosities were analyzed using the Jones–Dole equation and the viscosity A and B coefficients discussed in terms of solute–solute and solute–solvent interactions, respectively. Adiabatic compressibility (\( \beta_{s} \)) and apparent molar adiabatic compressibility (\( \phi_{\kappa }^{{}} \)), limiting apparent molar adiabatic compressibility (\( \phi_{\kappa }^{0} \)) and experimental slopes (\( S_{\kappa }^{*} \)) were determined from sound velocity data. The standard volume of transfer (\( \Delta_{t} \phi_{V}^{0} \)), viscosity B-coefficients of transfer (\( \Delta_{t} B \)) and limiting apparent molar adiabatic compressibility of transfer (\( \Delta_{t} \phi_{\kappa }^{0} \)) of metformin hydrochloride from water to aqueous glucose solutions were derived to understand various interactions in the ternary solutions. The activation parameters of viscous flow for the studied solutions were calculated using transition state theory. Hepler’s coefficient \( (d\phi /dT)_{p} \) indicated the structure making ability of metformin hydrochloride in the ternary solutions.     

Modeling simple amphiphilic solutes in a Jagla solvent

Methanol is an amphiphilic solute whose aqueous solutions exhibit distinctive physical properties. The volume change upon mixing, for example, is negative across the entire composition range, indicating strong association. We explore the corresponding behavior of a Jagla solvent, which has been previously shown to exhibit many of the anomalous properties of water. We consider two models of an amphiphilic solute: (i) a "dimer" model, which consists of one hydrophobic hard sphere linked to a Jagla particle with a permanent bond, and (ii) a "monomer" model, which is a limiting case of the dimer, formed by concentrically overlapping a hard sphere and a Jagla particle. Using discrete molecular dynamics, we calculate the thermodynamic properties of the resulting solutions. We systematically vary the set of parameters of the dimer and monomer models and find that one can readily reproduce the experimental behavior of the excess volume of the methanol-water system as a function of methanol volume fraction. We compare the pressure and temperature dependence of the excess volume and the excess enthalpy of both models with experimental data on methanol-water solutions and find qualitative agreement in most cases. We also investigate the solute effect on the temperature of maximum density and find that the effect of concentration is orders of magnitude stronger than measured experimentally.     

Effects of nonpolar solutes on the thermodynamic response functions of aqueous mixtures

We investigate the effect of adding nonpolar solutes at atmospheric pressure on water's temperature of maximum density, isothermal compressibility, and isobaric heat capacity, using a statistical mechanical model of water solutions [H. S. Ashbaugh, T. M. Truskett, and P. G. Debenedetti, J. Chem. Phys. 116, 2907 (2002)]. We find that the temperature of maximum density increases with solute hydrophobicity, as characterized by its size, and decreases with its van der Waals attractive parameter a, in agreement with experiment. We predict similar trends for the addition of solutes on the isothermal compressibility and isobaric heat capacity: solute hydrophobicity causes an upward shift in water's anomalies, whereas dispersive interactions as measured by the solute's van der Waals attractive parameter shift the anomalies to lower temperatures. The locus along which the competing contributions of solute size sigma and interaction strength a to the shift in water's response functions balance each other obeys the scaling relationship sigma6 approximately a.