The effect of aqueous solutions of trimethylamine-N-oxide on pressure induced modifications of hydrophobic interactions

To understand the mechanism of protein protection by the osmolyte trimethylamine-N-oxide (TMAO) at high pressure, using molecular dynamics (MD) simulations, solvation of hydrophobic group is probed in aqueous solutions of TMAO over a wide range of pressures relevant to protein denaturation. The hydrophobic solute considered in this study is neopentane which is a considerably large molecule. The concentrations of TMAO range from 0 to 4 M and for each TMAO concentration, simulations are performed at five different pressures ranging from 1 atm to 8000 atm. Potentials of mean force are calculated and the relative stability of solvent-separated state over the associated state of hydrophobic solute are estimated. Results suggest that high pressure reduces association of hydrophobic solutes. From computations of site-site radial distribution function followed by analysis of coordination number, it is found that water molecules are tightly packed around the nonpolar particle at high pressure and the hydration number increases with increasing pressure. On the other hand, neopentane interacts preferentially with TMAO over water and although hydration of neopentane reduces in presence of this osmolyte, TMAO does not show any tendency to prevent the pressure-induced dispersion of neopentane moieties. It is also observed that TMAO molecules prefer a side-on orientation near the neopentane surface, allowing its oxygen atom to form favorable hydrogen bonds with water while maintaining some hydrophobic contacts with neopentane. Analysis of hydrogen-bond properties and solvation characteristics of TMAO reveals that TMAO can form hydrogen bonds with water and it reduces the identical nearest neighbor water molecules caused by high hydrostatic pressures. Moreover, TMAO enhances life-time of water-water hydrogen bonds and makes these hydrogen bonds more attractive. Implication of these results for counteracting effect of TMAO against protein denaturation at high pressures are discussed.     

甘油水溶液氢键特性的分子动力学模拟

为了研究低温保护剂溶液的结构和物理化学特性, 以甘油为保护剂, 采用分子动力学方法, 对不同浓度的甘油和水的二元体系进行了模拟. 得到了不同浓度的甘油水溶液在2 ns内的分子动力学运动轨迹, 通过对后1 ns内运动轨迹的分析, 得到了各个原子对的径向分布函数和甘油分子的构型分布. 根据氢键的图形定义, 分析了氢键的结构和动力学特性. 计算了不同浓度下体系中平均每个原子(O和H)和分子(甘油和水)参与氢键个数的百分比分布及其平均值. 同时还计算了所有氢键、水分子之间的氢键以及甘油与水分子之间的氢键的生存周期.     

Anomalous Density and Permittivity Effects on the Structure of Water

The relationship between the structural peculiarities of the hydrogen bond net and the anomalous behavior of density and dielectric permittivity of water is investigated. The degree of ordering in the network of hydrogen bonds is described in terms of the structural functions, among which the most important ones are tetrahedricity and the mean number of hydrogen bonds per molecule. The temperature dependence of the number of hydrogen bonds per molecule is discussed in terms of the analysis of the experimental temperature dependences of density and permittivity of water on the saturation line. The estimated density of hexagonal ice at the melting point is reproduced in terms of the concepts suggested for analysis; the estimate is rather close to the experimental value. Possible applications of the new approach are discussed.     

Structural studies on the hydration of L-glutamic acid in solution

A combination of neutron diffraction augmented with isotopic substitution and computer modeling using empirical potential structure refinement has been used to extract detailed structural information for L-glutamic acid dissolved in 2 M NaOH solution. This work shows that the tetrahedral hydrogen bonding network in water is severely disrupted by the addition of glutamic acid and NaOH, with the number of water-water hydrogen bonds being reduced from 1.8 bonds per water molecule in pure water to 1.4 bonds per water molecule in the present solution. In the glutamic acid molecule, each carboxylate oxygen atom forms an average of three hydrogen bonds with the surrounding water solvent with one of these hydrogens being shared between the two oxygen atoms on each carboxylate group, while each amine hydrogen forms a single hydrogen bond with the surrounding water solvent. Additionally, the average conformation of the glutamic acid molecules in these solutions is extracted.     

Monte Carlo simulations of the solution structure of simple alcohols in water-acetonitrile mixtures

Monte Carlo simulations have been performed to explore the solution structure of ethyl, isopropyl, isobutyl, and tertiary butyl alcohols in pure water, pure acetonitrile, and different mixtures of the two solvents. The explicit solvent studies in NpT ensembles at T = 298 K illustrate that the solute "discriminates" the solvent's components and that the composition of the first solvation shell differs from that of the bulk solution. Since the polarizable continuum dielectric method (PCM) does not presently model the solvation of molecules with both polar and apolar sites in mixed protic solvents, we suggest a direction for further program development wherein a continuum dielectric method would accept more than one solvent and the solute sites would be solvated by user-defined solvent components. The prevailing solvation model will be determined upon the lowest free energy calculated for a particular solvation pattern of the solute having a specific conformational/tautomeric state. Characterization of equilibrium hydrogen-bond formation becomes a complicated problem that depends on the chemical properties of the solute and its conformation, as well as upon the varying nature of the first solvation shell. For example, while the number of hydrogen bonds to secondary and tertiary alcohol solutes are nearly constant in pure water and in water-acetonitrile mixtures with at least 50% water content, the number of hydrogen bonds to primary alcohols gradually decreases for most of their conformations when acetonitrile content is increased. Nonetheless, the calculations indicate that O-H...O(water) hydrogen bonds are still possible in a small fraction of the arrangements for the solution models with water content of 30% or less. The isopentene solute does not form any observable hydrogen bonds, despite having an electron-rich, double-bond site.     

Microwave dielectric permittivity and relaxation of aqueous potassium trimethylacetate solutions

The MW dielectric properties of aqueous potassium trimethylacetate (pivalate) solutions have been measured at six frequencies (10–25 GHz) at 288, 298, and 308 K. The static dielectric constants and dielectric relaxation times and activation parameters have been calculated. Trimethylacetate ion leads to the decrease in the mobility of water molecules and strengthening of their hydrogen bonds. These changes of water are similar to those in solutions of other carboxylates with a large number of nonpolar groups. The hydrophobic hydration of the trimethylacetate ion is maximal in this series.     

Rate and thermodynamic studies of the interaction between water and iso-butyl cellosolve by ultrasonic methods

Ultrasonic absorption and velocity measurements in aqueous solution of iso-butyl cellosolve (ethylene glycol iso-butyl ether) as a function of the concentration are reported. The two relaxational absorptions have been attributed to the perturbation of the equilibria expressed by AB?A+B and Aα(1/n)A_n where A is the solute, B is the solvent, AB is the complex and A _n is the solute aggregate. The rate constants for each step have been determined. From the concentration dependence of the maximum excess absorption per wave length, the enthalpy change and the volume change for the reaction between the solute and the solvent have been determined for aqueous solutions of butyl cellosolve (ethylene glycol n-butyl ether), iso-butyl cellosolve and propyl cellosolve (ethylene glycol n-propyl ether). The results are consistent with a hydrogen bonding reaction. The effect of the ethers on water structure are considered and it is clear that the fraction of water molecules which can hydrogen bond to the solute decreases with the increasing hydrophobicity of the solute.     

Influence of concentration and anion size on hydration of H+ ions and water structure

Neutron diffraction experiments with hydrogen isotope substitution on aqueous solutions of HCl and HBr have been performed at concentrations ranging from 1:17 to 1:83 solute per water molecules, at ambient conditions. Data are analyzed using the empirical potential structure refinement technique in order to extract information on both the ion hydration shells and the microscopic structure of the solvent. It is found that the influence of these solutes on the water structure is less concentration dependent than that of salts or hydroxides. Moreover protons readily form a strong H-bond with a water molecule upon solvation, at all proportions. The majority of them is also bonded via a longer bond to another water molecule, giving a prepeak in the g(OwOw). At high solute concentration, the second water molecule may be substituted by the counterion. In particular at solute concentrations of the order of 1:17 or higher, all protons have an anion within a distance of 4.5 A.     

Effects of ion concentration on the hydrogen bonded structure of water in the vicinity of ions in aqueous NaCl solutions

Molecular dynamics simulations of dilute and concentrated aqueous NaCl solutions are carried out to investigate the changes of the hydrogen bonded structures in the vicinity of ions for different ion concentrations. An analysis of the hydrogen bond population in the first and second solvation shells of the ions and in the bulk water is done. Although essentially no effect of ions on the hydrogen bonding is observed beyond the first solvation shell of the ions for the dilute solutions, for the concentrated solutions a noticeable change in the average number of water-water hydrogen bonds is observed in the second solvation shells of the ions and even beyond. However, the changes in the average number of hydrogen bonds are found to be relatively less when both water-water and ion-water hydrogen bonds are counted. Thus, the changes in the total number of hydrogen bonds per water are not very dramatic beyond the first solvation shell even for concentrated solutions.     

Aquaperoxovanadium(V) complexes: Quantum-chemical study

The addition of a water molecule to mono-, di-, and triperoxovanadium(V) complexes has been studied at the density functional theory (B3LYP/6-31G**) and Møller-Plesset perturbation theory (MP2/6-31G**) levels. It has been demonstrated that the H₂O…V donor-acceptor interaction cannot compete with hydrogen bonds and becomes weaker with an increase in the number of peroxo groups in the complex. In the most stable isomers of aquaperoxo complexes, water is mainly held by intermediate hydrogen bonds. The energy of addition of a water molecule to peroxovanadate is lower than or close to the heat of evaporation of water; i.e., the formation of stable aquaperoxovanadium complexes in aqueous solutions is improbable. This conclusion is consistent with the mass spectra of aqueous solutions of peroxovanadates, which show that the concentration of water-free peroxo complexes considerably exceeds the concentration of complexes with a coordinated water molecule. The coordination of the water molecule through the V…OH₂ donor-acceptor interaction is prevented by the cis effect of the peroxo group, which has the HOMO orbital of symmetry suitable for interaction with the LUMO orbital of the VO group.     

The role of electrolyte concentration in the decomposition of water and generation of electrons under conditions of anodic microdischarges

The dependence of the yield of water decomposition (from the measured yield of formation of hydrogen peroxide and yield of scavenging of OH radicals by isopropanol), as well as the yield of electron scavenging by a solute (nitrobenzene), on the concentration of aqueous electrolyte (sodium tetraborate and ammonium tetraborate) solutions under the action of anodic microdischarges was determined. Nonmonotonic functions were found, which exhibited maximums at 0.03 and 0.1 M concentrations for water decomposition and electrons, respectively. The results were explained by an increase in the fraction of anionic current through the discharge plasma-liquid cathode interface as the electrolyte concentration was increased.Translated from Khimiya Vysokikh Energii, Vol. 39, No. 2, 2005, pp. 140–142.Original Russian Text Copyright © 2005 by Polyakov, Badalyan, Bakhturova.This revised version was published online in April 2005 with a corrected cover date.     

Dependence of the number of hydrogen bonds per water molecule on its distance to a hydrophobic surface and a thereupon-based model for hydrophobic attraction

A water molecule in the vicinity of a hydrophobic surface forms fewer hydrogen bonds than a bulk molecule because the surface restricts the space available for other water molecules necessary for its hydrogen-bonding. In this vicinity, the number of hydrogen bonds per water molecule depends on its distance to the surface. Considering the number of hydrogen bonds per bulk water molecule (available experimentally) as the only reference quantity, we propose an improved probabilistic approach to water hydrogen-bonding that allows one to obtain an analytic expression for this dependence. (The original version of this approach [Y. S. Djikaev and E. Ruckenstein, J. Chem. Phys. 130, 124713 (2009)] provides the number of hydrogen bonds per water molecule in the vicinity of a hydrophobic surface as an average over all possible locations and orientations of the molecule.) This function (the number of hydrogen bonds per water molecule versus its distance to a hydrophobic surface) can be used to develop analytic models for the effect of hydrogen-bonding on the hydration of hydrophobic particles and their solvent-mediated interaction. Presenting a model for the latter, we also examine the temperature effect on the solvent-mediated interaction of two parallel hydrophobic plates.     

Dielectric spectroscopy in the GHz region on fully hydrated zwitterionic amino acids

The complex dielectric permittivity of eight different amino acids in water solutions was determined in the frequency range from 0.2 to 20 GHz at room temperature, trying to span the whole range of solubility in each case. Two relaxations were observed at room temperature in this frequency range, which can be mainly assigned to the rotation of amino acids in the aqueous environment, and the reorientational motion of water molecules, respectively. Although the amino acids have a charged (zwitterionic) nature with huge dipole moments, the tendency towards dipolar alignment seems to be very weak, over the investigated concentration ranges. For these small bio-molecules, water screens solute-solute interactions and amino acids remain typically as isolated hydrated monomers. The dielectric results were used to estimate the number of water molecules restrained by each solute molecule. Finally, the comparison between the amino acid relaxation times made it possible to discuss the relationship between rotational dynamics and the structure and hydrodynamic coupling of the amino acid studied.     

New Method for Determining H-Bond Energy: Fluorescence from Neutral Red in Water Compared to Fluorescence from Neutral Red in Benzene

A new method is presented for determining the H-bond energy between a fluorescent solute and water. In the current example the hydrogen-bond energy between neutral red (NR) and water was measured. The H-bond energy was obtained from the temperature dependence of the ratio of the 625 nm fluorescence intensity from NR in water, which increases with increasing temperature, to that of the 530 nm intensity from NR in benzene, which decreases with increasing temperature. The H-bond energy so determined is (10050 ± 60) J-mol⁻¹, which compares favorably with an H-bond E_HB value for pure water of 10170 J-mol⁻¹. Intermolecular hydrogen-bonding between water and NR increases the nonradiative deactivation process of the excited molecule due to fast energy dissipation through the vibrations associated with the hydrogen bonds. Some of the hydrogen bonds are broken with an increase in temperature, resulting in an increase in the number of free NR molecules, thus increasing the fluorescence quantum yield.     

水与乙醇分子在金管内吸附作用的分子动力学研究

使用分子动力学研究了乙醇与水分子在纳米金管内按照不同比例混合时的吸附现象,并利用径向密度分布函数及水和乙醇分子所形成的平均氢键数来探讨纳米限制效应.结果表明,径向密度分布函数和氢键数目受纳米金管影响较大.另外,水与金管之间的作用力比乙醇与金管之间的大,导致水分子形成的平均氢键数不同于乙醇分子的.     

Solvation free energies and hydration structure of N-methyl-p-nitroaniline

Solvation Gibbs energies of N-methyl-p-nitroaniline (MNA) in water and 1-octanol are calculated using the expanded ensemble molecular dynamics method with a force field taken from the literature. The accuracy of the free energy calculations is verified with the experimental Gibbs free energy data and found to reproduce the experimental 1-octanol∕water partition coefficient to within ±0.1 in log unit. To investigate the hydration structure around N-methyl-p-nitroaniline, an independent NVT molecular dynamics simulation was performed at ambient conditions. The local organization of water molecules around the solute MNA molecule was investigated using the radial distribution function (RDF), the coordination number, and the extent of hydrogen bonding. The spatial distribution functions (SDFs) show that the water molecules are distributed above and below the nitrogen atoms parallel to the plane of aromatic ring for both the methylamino and nitro functional groups. It is found that these groups have a significant effect on the hydration of MNA with water molecules forming two weak hydrogen bonds with both the methylamino and nitro groups. The hydration structures around the functional groups in MNA in water are different from those that have been found for methylamine, nitrobenzene, and benzene in aqueous solutions, and these differences together with weak hydrogen bonds explain the lower solubility of MNA in water. The RDFs together with SDFs provide a tool for the understanding the hydration of MNA (and other molecules) and therefore their solubility.     

The dielectric response of interfacial water—from the ordered structures to the single hydrated shell

Water is the universal solvent in nature. Does this imply, however, that its interaction with its environment is also a universal feature? While this question maybe too fundamental to be answered by one method only, we present evidence that the broadening of the dielectric spectra of water presents universal features of dipolar interactions with different types of matrixes. If in aqueous solutions the starting point of water’s state can be considered as bulk, with only partial interactions with the solute, then the state of water adsorbed in heterogeneous materials is determined by various hydration centers of the inhomogeneous material (the matrix) and it is significantly different from the bulk. In both cases, the dielectric spectrum of water is symmetrical and can be described by the Cole–Cole (CC) function. The phenomenological model that describes a physical mechanism of the dipole–matrix interaction in complex systems underlying the CC behavior has been applied to water adsorbed in porous glasses. It was then extended to analyses of the dynamic and structural behavior of water in nonionic and ionic aqueous solutions. The same model is then used to analyze the CC relaxation processes observed in clays, aqueous solutions of nucleotides, and amino acids.     

Structural changes of the molecular complexes and consolute phenomena in binary liquid mixtures

Consolute phenomena in the aqueous solutions of the polymers are considered In view of the temperature induced structural changes of the hydrogen bonds between water and functional groups of polymer. The lower and upper critical consolute points are attributed to the appearance of the “critical” concentration of the complexes with one hydrogen bond between single water molecule and functional group of polymer. Namely such kind of the hydrogen bonds are responsible for the formation of the strongly associated water clusters, that may be followed by phase separation. Experimentally observed dependences of the critical consolute temperatures for the aqueous solutions of polyethylene glycol on the molecular weight of polymer and adding of salts are well reproduced in the framework of the proposed model.     

Molecular dynamics simulation of cellobiose in water

The conformational behavior of cellobiose was studied by molecular dynamics simulation in a periodic box of waters. Several different initial conformations were used and the results compared with equivalent vacuum simulations. The average positions and rms fluctuations within single torsional conformations of cellobiose were affected only slightly by the solvent. However, water damped local torsional librations and transitions. The conformational energies of the solute and their fluctuations were also sensitive to the presence of solvent. Intramolecular hydrogen bonding was weakened relative to that observed in vacuo due to competition with solvating waters. All cellobiose hydroxyl groups participated in intermolecular hydrogen bonds with water, with approximately eight hydrogen bonds formed per glucose ring. The hydrogen bonding was predominantly between water hydrogens and solute hydroxyl oxygens. Intermolecular hydrogen bonding to ring and bridge oxygens was seldom present. The diffusion coefficients of both water and solute agree closely with experimental values. Water interchanged rapidly between the solvating first shell and the bulk on the picosecond time scale. © 1993 John Wiley & Sons, Inc.     

甘油对冰晶生长抑制作用的分子动力学模拟

采用分子模拟方法研究了正交晶系冰晶(020)生长面在不同浓度甘油水溶液中的生长情况. 通过统计分析氢键数、 密度分布函数、 均方根偏差和原子间径向分布函数研究了水分子和甘油分子的动态行为. 结果表明, 甘油分子在水溶液中可与水分子形成大量氢键, 这使水分子间的氢键作用受到抑制, 降低了水分子的扩散性, 致使冰晶不易成核和生长; 另外, 一些甘油分子可代替水分子吸附在晶面上, 甚至占据晶格位点, 这种行为打破了冰晶的对称性并且降低了冰晶的生长速率. 因此, 甘油可同时在晶面和液相中抑制冰晶的生长.