Aqueous sodium oleate–sodium dehydrocholate mixtures at low concentration

The aqueous mixed system sodium dehydrocholate (NaDHC)–sodium oleate (NaOL) was studied by several methods to determine the influence of the hydrophobic structure of both surfactants in the mixed micellization and the formation of the mixed monolayer adsorbed at the air–water interface. The molecular area at the critical micelle concentration in pure surfactant solutions suggests that the adsorbed oleate chain was folded to allow the double bond in the middle of the molecule to remain in contact with water, and that the NaDHC molecule was situated with its plane laying parallel to the water surface, allowing the three carbonyl groups in the hydrocarbon backbone to form hydrogen bonds with water. The interaction was repulsive at the surface, and in the mixed monolayer some molecules must move away the less hydrophilic groups from water (double bond of NaOL, carbonyl groups of NaDHC). The interaction in mixed micelles was strongly attractive, showing a preferential composition roughly equimolar. The hydrolysis in mixed micelles was augmented in comparison with pure surfactants systems, which could be explained by assuming the existence of a more hydrophobic mixed micelle core. The mixed micelle degree of ionization was below that of the pure micelles, thus indicating a high surface charge density.     

Solvation and hydrophobic hydration of different types of alkylamines inN,N-dimethylformamide and water mixtures

Enthalpies of solution of twelve amines of different type have been determined at 25°C in mixtures of N,N-dimethylformamide and water over the whole composition range. The enthalpies of transfer from water to the mixtures deviate substantially from a linear dependence on the mole fraction of water. These deviations appear to contain additive contributions of the different alkyl groups. By application of a simple hydration model the enthalpic effect of hydrophobic hydration has been calculated for each amine. For alkylamines this is determined by the number and size of the alkyl groups present in the molecule. The contribution of each alkyl group is the same in primary, secondary and tertiary amines. Results for the different alkyl groups show a close relationship with values for alcohols obtained previously. Differences between alcohols and amines can be attributed to differences in the hydrophobic hydration of the parts of the solute molecules which are adjacent to the polar group. The influence of the polar group does not seem to extend beyond the second carbon atom.     

WATsite: Hydration site prediction program with PyMOL interface

Water molecules that mediate protein–ligand interactions or are released from the binding site on ligand binding can contribute both enthalpically and entropically to the free energy of ligand binding. To elucidate the thermodynamic profile of individual water molecules and their potential contribution to ligand binding, a hydration site analysis program WATsite was developed together with an easy‐to‐use graphical user interface based on PyMOL. WATsite identifies hydration sites from a molecular dynamics simulation trajectory with explicit water molecules. The free energy profile of each hydration site is estimated by computing the enthalpy and entropy of the water molecule occupying a hydration site throughout the simulation. The results of the hydration site analysis can be displayed in PyMOL. A key feature of WATsite is that it is able to estimate the protein desolvation free energy for any user specified ligand. The WATsite program and its PyMOL plugin are available free of charge from http://people.pnhs.purdue.edu/~mlill/software . © 2014 Wiley Periodicals, Inc.     

Lignite humic acids aggregates studied by high resolution ultrasonic spectroscopy

The thermodynamic stability of lignite humic acids (sodium salt) aggregates was studied by high resolution ultrasonic spectroscopy within the temperature interval from 5 to 90°C. The changes in differential ultrasonic velocity (U12) showed strong differences among humic solutions within the concentration range from 0.005 to 10 g L⁻¹. Measurement revealed several transitions which were attributed to the weakening of humic secondary structure. Concentration around 1 g L⁻¹ seemed to be a limit under which the change of the prevalence and importance of hydration occurred. Above this concentration the difference in U12 decreased following the temperature increase which was explained as a dominance of hydrophilic hydration. In contrast, below this concentration, the temperature dependence of U12 resulted in increasing tendency which was attributed to the prevalence of hydrophobic hydration, i.e. uncovering of apolar groups towards surrounding water. Additional experiments in which the humic sample was modified by hydrochloric acid resulted in a slight structural stabilization which lead to the conclusion that humic micelle-like subaggregates form an open-layer assemblies easily accessible for interaction with an extraneous molecule. That was partly verified by addition of propionic acid which brought about even larger reconformation of humic aggregates and exhibition of polar groups towards hydration water. The reversible changes in humate solutions induced by elevated temperatures provided the evidence about the existence of significant physical interactions among humic molecules resulting in formation of various kinds of aggregates. The nature of aggregates, mainly the stability and conformation, strongly depends on the concentration. Evidently, the changes observed in this work cannot be simply explained as expansions or conformational changes of macromolecular coils.     

Hydration of phosphatidyl serine multilayers and its modulation by conformational change induced by correlated electrostatic interaction.

Hydration of the various residues of phospholipids was inferred from the shift in the wave number of their vibration bands, obtained from the amplitudes of their positive and negative peaks in the difference spectra between those of the hydrated and the dry phospholipid multibilayers. The effect of aligned phospholipid layers on the orientation of their hydrating water molecules was inferred from the dichroic ratio of the OH stretching band, measured by polarized attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) with a germanium prism, as a function of the water-to-lipid ratio in the surface film. The results indicate that about seven water molecules are oriented by one phosphatidyl serine molecule in the surface film. About 8 to 11 additional water molecules contribute to the hydration of the polar residues as revealed by the effect on the difference spectra. The hydration appears to be cooperative. A water molecule that initiates hydration of a site facilitates access of additional water molecules, until the hydration of the whole site composed of many different interacting polar residues is completed.     

Potential of mean force of hydrophobic association: dependence on solute size

The potentials of mean force (PMFs) were determined for systems involving formation of nonpolar dimers composed of methane, ethane, propane, isobutane, and neopentane, respectively, in water, using the TIP3P water model, and in vacuo. A series of umbrella-sampling molecular dynamics simulations with the AMBER force field was carried out for each pair in either water or in vacuo. The PMFs were calculated by using the weighted histogram analysis method (WHAM). The shape of the PMFs for dimers of all five nonpolar molecules is characteristic of hydrophobic interactions with contact and solvent-separated minima and desolvation maxima. The positions of all these minima and maxima change with the size of the nonpolar molecule, that is, for larger molecules they shift toward larger distances. The PMF of the neopentane dimer is similar to those of other small nonpolar molecules studied in this work, and hence the neopentane dimer is too small to be treated as a nanoscale hydrophobic object. The solvent contribution to the PMF was also computed by subtracting the PMF determined in vacuo from the PMF in explicit solvent. The molecular surface area model correctly describes the solvent contribution to the PMF together with the changes of the height and positions of the desolvation barrier for all dimers investigated. The water molecules in the first solvation sphere of the dimer are more ordered compared to bulk water, with their dipole moments pointing away from the surface of the dimer. The average number of hydrogen bonds per water molecule in this first hydration shell is smaller compared to that in bulk water, which can be explained by coordination of water molecules to the hydrocarbon surface. In the second hydration shell, the average number of hydrogen bonds is greater compared to bulk water, which can be explained by increased ordering of water from the first hydration shell; the net effect is more efficient hydrogen bonding between the water molecules in the first and second hydration shells.     

Models of hydration and isomeric transitions of glucose molecules in aqueous solutions

The kinematic viscosity of aqueous glucose solutions is studied. It is found that the hydrodynamic radius of monosaccharide molecule in an aqueous solution depends on temperature in the range of 290–355 K. Using a bimodal model of the energy states of the volume in which the glucose molecule is located and local equilibrium is established, it is shown that the above-mentioned dependence can be attributed to disturbances in the equilibrium of isomeric transitions, induced by variations in temperature. The parameters of isomeric transitions for a glucose molecule in an aqueous solvent, the probability of “chair” and “boat” configurations occurring for glucose molecules, and the number of water molecules in the hydration shells of these configurations are calculated; the strain of the chemical bonds in the chair configuration of a glucose molecule is estimated.     

The influence of microhydration on the ionization energy thresholds of uracil and thymine

In the present study the ionization energy thresholds (IET's) of uracil and thymine have been calculated (with the B3LYP, PMP2, and P3 levels of theory using the standard 6-31++G(d,p) basis set) with one to three water molecules placed in the first hydration shell. Then (B3LYP) polarizable continuum model (PCM) calculations were performed with one to three waters of the hydration shell included. Calculations show there is a distinct effect of microhydration on uracil and thymine. For uracil, one added water results in a decrease in the IET of about 0.15 eV. The second and third water molecules cause a further decrease by about 0.07 eV each. For thymine, the first water molecule is seen to decrease the IET by about 0.1 eV, while the second and third water molecules cause a further decrease of less than 0.1 eV each. The changes in IET calculated here for thymine with one to three waters of hydration are smaller than the experimental values determined by Kim et al. (Kim, S. K.; Lee, W.; Herschbach, D. R. J. Phys. Chem. 1996, 100, 7933). Preliminary results presented here indicate that the experimental results may involve keto-enol tautomers of thymine. The results of placing the microhydrated structures of uracil and thymine in a PCM cavity was seen to make very little difference in the IET when compared to the IET of ordinary uracil or thymine in a PCM cavity. The implications are that accurate calculations of the IET's of uracil and thymine can be obtained by simply considering long-range solvation effects.     

Molecular dynamic study of hydrophilicity of the NaCl molecule

The NaCl molecule is represented as a physical dipole in which the distance between the ions equals that in the crystal lattice of NaCl. The water molecules, initially uniformly distributed in the sphere around the NaCl molecule, in 1 or 2 ps form a hydration shell around the positively charged sodium ion, leaving the negatively charged chloride ion naked. When the number of water molecules around the sodium ion reaches 14, the NaCl(H₂O)₁₄ cluster becomes sensitive only to thermal perturbations that are due to the librational motion of water molecules. It is shown that the dependence of the cluster volume on the number of water molecules in the cluster does not admit simple linear approximation. Institute of Thermal Physics, Ural Branch, Russian Academy of Sciences. Translated fromZhurnal Strukturnoi Khimii, Vol. 37, No. 2, pp. 299–309, March–April, 1996. Translated by L. Smolina     

Microscopic wetting of self-assembled monolayers with different surfaces: a combined molecular dynamics and quantum mechanics study

Molecular dynamics simulations are used to study the micronature of the organization of water molecules on the flat surface of well-ordered self-assembled monolayers (SAMs) of 18-carbon alkanethiolate chains bound to a silicon (111) substrate. Six different headgroups (-CH(3), -C═C, -OCH(3), -CN, -NH(2), -COOH) are used to tune the character of the surface from hydrophobic to hydrophilic, while the level of hydration is consistent on all six SAM surfaces. Quantum mechanics calculations are employed to optimize each alkyl chain of the different SAMs with one water molecule and to investigate changes in the configuration of each headgroup under hydration. We report the changes of the structure of the six SAMs with different surfaces in the presence of water, and the area of the wetted surface of each SAM, depending on the terminal group. Our results suggest that a corrugated and hydrophobic surface will be formed if the headgroups of SAM surface are not able to form H-bonds either with water molecules or between adjacent groups. In contrast, the formation of hydrogen bonds not only among polar heads but also between polar heads and water may enhance the SAM surface hydrophilicity and corrugation. We explicitly discuss the micromechanisms for the hydration of three hydrophilic SAM (CN-, NH(2)- and COOH-terminated) surfaces, which is helpful to superhydrophilic surface design of SAM in biomimetic materials.     

Gas-Phase Reactions of hydrated alkaline earth metal ions, M2+(H2O) n (M = Mg, Ca, Sr, Ba and n = 4–7), with benzene

Gas-phase reactions of hydrated divalent alkaline earth metal ions and benzene were investigated by electrospray ionization Fourier-transform mass spectrometry. Rate constants for solvent-exchange reactions were determined as a function of hydration extent for Mg2+, Ca2+, Sr2+, and Ba2+ clusters containing four to seven water molecules each. All of the strontium and barium clusters react quickly with benzene. Barium reacts slightly faster than the corresponding strontium cluster with the same number of water molecules attached. For calcium, clusters with four and five water molecules react quickly, whereas those with six and seven water molecules do not. Magnesium with four water molecules reacts quickly, but not when five through seven water molecules are attached. The slow reactivity observed for some of these clusters indicates that the cation-pi interaction between the metal ion and benzene is partially screened by the surrounding water molecules. The reactivity of magnesium with seven water molecules is intermediate that of the hexa- and pentahydrate and the tetrahydrate. This result is consistent with the seventh water molecule being in the outer shell and much more weakly bound. The unusual trend in reactivity observed for magnesium may be due to the presence of mixed shell structures observed previously. These results are the first to provide information about the relative importance of cation-pi interactions in divalent metal ions as a function of metal hydration extent. Such studies should also provide a model and some insight into the relative binding affinities of divalent metal ions to aromatic residues on peptides and proteins.     

Dielectric spectroscopy study on ionic liquid microemulsion composed of water, TX-100, and BmimPF6

We report here a broadband dielectric spectroscopy study on an ionic liquid microemulsion (ILM) composed of water, Triton X-100 (TX-100), and 1-butyl-3-methylimidazolium hexafluorophosphate (bmimPF(6)). It is found that the phase behavior of this ILM can be easily identified by its dielectric response. The dielectric behavior of the ILM in the GHz range is consistent with that of TX-100∕water mixtures with comparable water-to-TX-100 weight ratio. It consists of the relaxations due to ethylene oxide (EO) unit relaxation, hydration water dynamics, and∕or free water dynamics. The water content dependence of the EO unit relaxation suggests that this relaxation involves dynamics of hydration water molecules. In the IL-in-water microemulsion phase, it is found that bmimPF(6) molecules are preferentially dissolved in water when their concentration in water is lower than the solubility. An additional dielectric relaxation that is absent in the TX-100∕water mixtures is observed in the frequency range of 10(7)-10(8) Hz for this ILM. This low-frequency relaxation is found closely related to the bmimPF(6) molecule and could be attributed to the hopping of its cations∕anions between the anionic∕cationic sites.     

Monte Carlo study of the aqueous hydration of dimethylphosphate conformations

Monte Carlo computer simulations were performed on dilute aqueous solutions of the dimethylphosphate anion and the sodium dimethylphosphate ion pair, with the two phosphodiester torsional angles in the gauche–gauche, gauche–trans, and trans–trans conformations. The structural and energetic aspects of the aqueous hydration of each molecule were analyzed in terms of quasi component distribution functions based on the proximity criterion and partitioned into ionic, hydrophilic, and hydrophobic contributions to facilitate an understanding of the hydration pattern and conformational trends in these multifunctional solutes. Special attention was also paid to methodological issues affecting hydration, such as statistical uncertainty in the determined hydration indices, choice of partial atomic charges for the solute atoms, and solute–water interaction potentials adopted in the simulations. The results showed that gauche–trans and gauche–gauche forms are equally favorable for the dimethylphosphate anion with the trans extended form destabilized by hydration. The sodium dimethylphosphate ion pair hydration energetically favors the trans–trans conformation.     

Desorption of Water from CD/DRUG Inclusion Complexes

Thermogravimetric and differential scanning calorimetric traces were recorded for several crystalline hydrated cyclodextrin complexes containing drug substances (salts of diclofenac, meclofenamate sodium, (L)-menthol) as guests. It was possible to reconcile observed thermal events for complex dehydration with three dimensional hydration patterns deduced from available X-ray crystal structures of the complexes. For complexes containing drug salts, strong retention of water molecules is attributed to their coordination to metal ions.     

Monolayer of aerosol-OT surfactants adsorbed at the air/water interface: an atomistic computer simulation study

An atomistic molecular dynamics (MD) simulation has been carried out to investigate the structural and dynamical properties of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature and at a surface coverage corresponding to that at its critical micelle concentration (78 A(2)/molecule). The estimated thickness of the adsorbed layer is in good agreement with neutron reflection data. The study shows that the surfactants exhibit diffusive motion in the plane of the interface. It is observed that the surfactant monolayer has a strong influence in restricting both the translational and reorientational motions of the water molecules close to the interface. A drastic difference in the dipolar reorientational motion of water molecules in the aqueous layer is observed with a small variation of the distance from the surfactant headgroups. It has been observed that the water molecules in the first hydration layer (region 1) form strong hydrogen bonds with surfactant headgoups. This results in the slower structural relaxation of water-water hydrogen bonds in the first hydration layer compared to that in the pure bulk water. Most interestingly, we notice that the water molecules present in the layer immediately after the first hydration layer form weaker hydrogen bonds and thus relax faster than even pure bulk water.     

Two-particle entropy and structural ordering in liquid water

Entropies of simple point charge (SPC) water were calculated over the temperature range 278-363 K using the two-particle correlation function approximation. Then, the total two-particle contribution to the entropy of the system was divided into three parts, which we call translational, configurational, and orientational. The configurational term describes the contribution to entropy, which originates from spatial distribution of surrounding water molecules (treated as points, represented by the center of mass) around the central one. It has been shown that this term can serve as the metric of the overall orientational ordering in liquid water. Analyzing each of these three terms as a function of intermolecular distance, r, we also find a rational definition of the hydration shell around the water molecule; the estimated radii of the first and second hydration shells are 0.35 nm and 0.58 nm, respectively. We find, moreover, that the first hydration shell around the water molecule participates roughly in 70% of the total orientational entropy of water, and this rate is roughly temperature independent.     

Hydration of Bromine Ions in Water

A mass spectrographic method electrospraying electrolyte solutions in vacuum was used to determine the enthalpy of hydration of bromine ions by different numbers of water molecules. The mass spectra of negative ions from aqueous solutions of NaBr were measured over the temperature range of 15–50°C and showed that the most intense peaks correspond to the hydrated bromine ions with one and two water molecules. It was found that the Van't Hoff equation holds for bromine hydrates containing up to three water molecules, whereas addition of the fourth water molecule leads to a violation of the Van't Hoff equation. The enthalpy of bromine ion hydration by the first and the second water molecules was found to be –(28.5 ± 1.7) and –(21 ± 4.2) kJ/mol, respectively.     

Extracting hydration sites around proteins from explicit water simulations

Two new methods are assessed for determining the location of hydration sites around proteins from computer simulation. Current methods extract hydration sites from peaks in the water density constructed in the protein frame. However, the dynamic nature of the water molecules, the nearby protein residues, and the protein reference frame as a whole tend to smear out the water density, making it more difficult to resolve sites. Two techniques are introduced to better resolve the water density. The first is to construct the water density from the time-averaged position of each water molecule in the protein frame while the water remains within a given distance of this averaged position. The second technique is to construct the water density from the time-averaged position of each water in the reference frame only of the nearby residues. Criteria for determining hydration sites from the water density are examined. Both techniques are found to significantly improve the detail in the water density and the number of hydration sites detected.     

Adsorption of water molecules on selected charged sodium-chloride clusters

The adsorption of water molecules (H(2)O) on sodium chloride cluster cations and anions was studied at 298 K over a mass range of 100-1200 amu using a custom-built laser desorption ionization reactor and mass spectrometer. Under the conditions used, the cations Na(3)Cl(2)(+) and Na(4)Cl(3)(+) bind up to three water molecules, whereas the larger cations, Na(5)Cl(4)(+) to Na(19)Cl(18)(+), formed hydrates with one or two only. The overall trend is a decrease in hydration with increasing cluster size, with an abrupt drop occurring at the closed-shell Na(14)Cl(13)(+). As compared to the cluster cations, the cluster anions showed almost no adsorption. Among smaller clusters, a weak adsorption of one water molecule was observed for the cluster anions Na(6)Cl(7)(-) and Na(7)Cl(8)(-). In the higher mass region, a substantial adsorption of one water molecule was observed for Na(14)Cl(15)(-). Density functional theory (DFT) computations were carried out for the adsorption of one molecule of H(2)O on the cations Na(n)Cl(n-1)(+), for n = 2-8, and the anions Na(n)Cl(n+1)(-), for n = 1-7. For each ion, the structure of the hydrate, the hydration energy, and the standard-state enthalpy, entropy, and Gibbs energy of hydration at 298 K were computed. In addition, it was useful to compute the distortion energy, defined as the electronic energy lost due to weakening of the Na-Cl bonds upon adsorption of H(2)O. The results show that strong adsorption of a H(2)O molecule occurs for the linear cations only at an end Na ion and for the nonlinear cations only at a corner Na ion bonded to two Cl ions. An unexpected result of the theoretical investigation for the anions is that certain low-energy isomers of Na(6)Cl(7)(-) and Na(7)Cl(8)(-) bind H(2)O strongly enough to produce the observed weak adsorption. The possible implications of these results for the initial hydration of extended NaCl surfaces are discussed.