Dynamics of dilute water in carbon tetrachloride

A dilute solution of water in a hydrophobic solvent, such as carbon tetrachloride (CCl4), presents an opportunity to study the rotational properties of water without the complicating effects of hydrogen bonds. We report here the results of theoretical, experimental, and semiempirical studies of a 0.03 mole percent solution of water in CCl4. It is shown that for this solution there are negligible water-water interactions or water-CCl4 interactions; theoretical and experimental values for proton NMR chemical shifts (deltaH) are used to confirm the minimal interactions between water and the CCl4. Calculated ab initio values and semiempirical values for oxygen-17 and deuterium quadrupole coupling constants (chi) of water/CCl4 clusters are reported. Experimental values for the 17O, 2H, and 1H NMR spin-lattice relaxation times, T1, of 0.03 mole percent water in dilute CCl4 solution at 291 K are 94+/-3 ms, 7.0+/-0.2 s, and 12.6+/-0.4 s, respectively. These T1 values for bulk water are also referenced. "Experimental" values for the quadrupole coupling constants and relaxation times are used to obtain accurate, experimental values for the rotational correlation times for two orthogonal vectors in the water molecule. The average correlation time, tauc, for the position vector of 17O (orthogonal to the plane of the molecule) in monomer water, H2(17)O, is 91 fs. The average value for the deuterium correlation time for the deuterium vector in 2H2O is 104 fs; this vector is along the OD bond. These values indicate that the motion of monomer water in CCl4 is anisotropic. At 291 K, the oxygen rotational correlation time in bulk 2H2(17)O is 2.4 ps, the deuterium rotational correlation time in the same molecule is 3.25 ps. (Ropp, J.; Lawrence, C.; Farrar, T. C.; Skinner, J. L. J. Am. Chem. Soc. 2001, 123, 8047.) These values are a factor of about 20 longer than the tauc value for dilute monomer water in CCl4.     

Rotational motion in liquid water is anisotropic: a nuclear magnetic resonance and molecular dynamics simulation study

Experimental NMR measurements of the deuterium and (17)O T(1) relaxation times in deuterium-enriched liquid water have been performed from 275 to 350 K. These relaxation times can yield rotational correlation times of appropriate molecule-fixed unit vectors if the quadrupole coupling constants and asymmetry parameters are known. We determine the latter from ab initio studies of water clusters and experimental chemical shift measurements. We find that the rotational correlation time for the OD bond vector in D(2)(16)O varies from 5.8 ps at 275 K to 0.86 ps at 350 K, and that the rotational correlation time for the out-of-plane vector of dilute D(2)(17)O in D(2)(16)O varies from 4.4 ps at 275 K to 0.64 ps at 350 K. These results indicate that the rotational motion of water is anisotropic. Molecular dynamics simulations of liquid water are in good agreement with these experiments at the higher temperatures, but the simulation results are considerably faster than experiment at the lower temperatures.     

Dynamics and thermodynamics of water in PAMAM dendrimers at subnanosecond time scales

Atomistic molecular dynamics simulations are used to study generation 5 polyamidoamine (PAMAM) dendrimers immersed in a bath of water. We interpret the results in terms of three classes of water: buried water well inside of the dendrimer surface, surface water associated with the dendrimer-water interface, and bulk water well outside of the dendrimer. We studied the dynamic and thermodynamic properties of the water at three pH values: high pH with none of the primary or tertiary amines protonated, intermediate pH with only the primary amines protonated, and low pH with all amines protonated. For all pH values we find that both buried and surface water exhibit two relaxation times: a fast relaxation ( approximately 1 ps) corresponding to the libration motion of the water and a slow ( approximately 20 ps) diffusional component related to the escaping of water from one domain to another. In contrast for bulk water the fast relaxation is approximately 0.4 ps while the slow relaxation is approximately 14 ps. These results are similar to those found in biological systems, where the fast relaxation is found to be approximately 1 ps while the slow relaxation ranges from 20 to 1000 ps. We used the 2PT MD method to extract the vibrational (power) spectrum and found substantial differences for the three classes of water. The translational diffusion coefficient for buried water is 11-33% (depending on pH) of the bulk value while the surface water is about 80%. The change in rotational diffusion is quite similar: 21-45% of the bulk value for buried water and 80% for surface water. This shows that translational and rotational dynamics of water are affected by the PAMAM-water interactions as well as due to the confinement in the interior of the dendrimer. We find that the reduction of translational or rotational diffusion is accompanied by a blue shift of the corresponding libration motions ( approximately 10 cm(-1) for translation, approximately 35 cm(-1) for rotation), indicating higher local force constants for these motions. These effects are most pronounced for the lowest pH, probably because of the increased rigidity caused by the internal charges. From the vibrational density of states we also calculate the enthalpies and entropies of the various waters. We find that water molecules are enthalpically favored near the PAMAM dendrimer: energy for surface water is approximately 0.1 kcal/mol lower to that in the bulk, and approximately 0.5-0.9 kcal/mol lower for buried water. In contrast, we find that both the buried and surface water are entropically unfavored: buried water is 0.9-2.2 kcal/mol lower than the bulk while the surface water is 0.1-0.2 kcal/mol lower. The net result is a thermodynamically unfavored state of the water surrounding the PAMAM dendrimer: 0.4-1.3 kcal/mol higher for buried water and 0.1-0.2 kcal/mol for surface water. This excess free energy of the surface and buried waters is released when the PAMAM dendrimer binds to DNA or metal ions, providing an extra driving force.     

The combined effect of cations and anions on the dynamics of water

We studied the orientational relaxation of the OD-stretch vibration of HDO molecules in concentrated solutions of alkali-halide salts (NaCl, NaI, CsCl and KI) in isotopically diluted water (4% D(2)O in H(2)O), using polarization-resolved femtosecond infrared pump-probe spectroscopy (fs-IR). We were able to distinguish the orientational dynamics of the water molecules solvating the halide ions from the dynamics of the bulk water and the water solvating the cations. We found that the reorientation of the halide-bound molecules shows two strongly different components (2.0 ± 0.3 ps and 9 ± 1 ps), related to a wiggling motion of the OD group hydrogen-bonded to the anion, and rotational diffusion of the molecule over the charged anion surface, respectively. The relative amplitudes of the two components are dependent on the nature of both the anion and cation, and on the concentration. These results show that cations can have a profound effect on the solvation shell dynamics of their counter-ions.     

Spectral diffusion in a fluctuating charge model of water

We apply the combined electronic structure/molecular dynamics approach of Corcelli, Lawrence, and Skinner [J. Chem. Phys. 120, 8107 (2004)] to the fluctuating charge (SPC-FQ) model of liquid water developed by Rick, Stuart, and Berne [J. Chem. Phys. 101, 6141 (1994)]. For HOD in H(2)O the time scale for the long-time decay of the OD stretch frequency time-correlation function, which corresponds to the time scale for hydrogen-bond rearrangement in the liquid, is about 1.5 ps. This result is significantly longer than the 0.9 ps decay previously calculated for the nonpolarizable SPC/E water model. Our results for the SPC-FQ model are in better agreement with recent vibrational echo experiments.     

13C NMR, infrared, solvation and theoretical investigation of the conformational isomerism in 1-haloacetones (X = Cl, Br and I)

The solvent dependence of the 13C NMR spectra of chloroacetone (CA), bromoacetone (BA) and iodoacetone (IA) are reported and the 3J(CH) couplings analysed using ab initio calculations and solvation theory. In CA the energy difference (E(cis) - E(gauche)) between the cis (Cl-C-C=O 0 degrees) and gauche (Cl-C-C=O 155 degrees) conformers is 1.7 kcal mol(-1) in the vapour, decreasing to 0.8 kcal mol(-1) in CCl4 solution and to -1.0 kcal mol(-1) in the pure liquid. The conformational equilibrium, in BA, is between the more polar cis (Br-C-C=O 0 degrees) and gauche (Br-C-C=O 132 degrees) conformations. The energy difference (E(cis) - E(gauche)) is 1.8 kcal mol(-1) in the vapour, decreasing to 0.9 kcal mol(-1) in CCl4 solution and to -0.4 kcal mol(-1) in the pure liquid. The energy difference (E(cis) - E(gauche)), in IA, between the cis (I-C-C=O 0 degrees) and gauche (I-C-C=O 104 degrees) conformers is 1.1 kcal mol(-1) in the vapour phase, decreasing to 0.5 kcal mol(-1) in CCl4 solution and to -0.5 kcal mol(-1) in the pure liquid. The vapour state energy difference for BA [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)] and for IA [1.6 kcal mol(-1) at B3LYP/6-311++G(d,p)/LANL2DZ)] are in very good agreement with the above values. For CA the agreement is also satisfactory [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)].     

Theoretical calculations of homoconjugation equilibrium constants in systems modeling acid-base interactions in side chains of biomolecules using the potential of mean force

The potentials of mean force (PMFs) were determined for systems forming cationic and anionic homocomplexes composed of acetic acid, phenol, isopropylamine, n-butylamine, imidazole, and 4(5)-methylimidazole, and their conjugated bases or acids, respectively, in three solvents with different polarity and hydrogen-bonding propensity: acetonitrile (AN), dimethyl sulfoxide (DMSO), and water (H(2)O). For each pair and each solvent a series of umbrella-sampling molecular dynamics simulations with the AMBER force field, explicit solvent, and counterions added to maintain a zero net charge of a system were carried out and the PMF was calculated by using the Weighted Histogram Analysis Method (WHAM). Subsequently, homoconjugation-equilibrium constants were calculated by numerical integration of the respective PMF profiles. In all cases but imidazole stable homocomplexes were found to form in solution, which was manifested as the presence of contact minima corresponding to hydrogen-bonded species in the PMF curves. The calculated homoconjugation constants were found to be greater for complexes with the OHO bridge (acetic acid and phenol) than with the NHN bridge and they were found to decrease with increasing polarity and hydrogen-bonding propensity of the solvent (i.e., in the series AN > DMSO > H(2)O), both facts being in agreement with the available experimental data. It was also found that interactions with counterions are manifested as the broadening of the contact minimum or appearance of additional minima in the PMF profiles of the acetic acid-acetate, phenol/phenolate system in acetonitrile, and the 4(5)-methylimidazole/4(5)-methylimidzole cation conjugated base system in dimethyl sulfoxide.     

Vibrational and orientational dynamics of water in aqueous hydroxide solutions

We report the vibrational and orientational dynamics of water molecules in isotopically diluted NaOH and NaOD solutions using polarization-resolved femtosecond vibrational spectroscopy and terahertz time-domain dielectric relaxation measurements. We observe a speed-up of the vibrational relaxation of the O-D stretching vibration of HDO molecules outside the first hydration shell of OH(-) from 1.7 ± 0.2 ps for neat water to 1.0 ± 0.2 ps for a solution of 5 M NaOH in HDO:H(2)O. For the O-H vibration of HDO molecules outside the first hydration shell of OD(-), we observe a similar speed-up from 750 ± 50 fs to 600 ± 50 fs for a solution of 6 M NaOD in HDO:D(2)O. The acceleration of the decay is assigned to fluctuations in the energy levels of the HDO molecules due to charge transfer events and charge fluctuations. The reorientation dynamics of water molecules outside the first hydration shell are observed to show the same time constant of 2.5 ± 0.2 ps as in bulk liquid water, indicating that there is no long range effect of the hydroxide ion on the hydrogen-bond structure of liquid water. The terahertz dielectric relaxation experiments show that the transfer of the hydroxide ion through liquid water involves the simultaneous motion of ~7 surrounding water molecules, considerably less than previously reported for the proton.     

Rotational structure of water in a hydrophobic environment: carbon tetrachloride

Infrared spectroscopy has been used to probe the interaction between water and the hydrophobic solvent, carbon tetrachloride. At room temperature, water exists as monomers in carbon tetrachloride, presenting a system for studying the rotational properties of water free of strong hydrogen-bonding. The rotational structure suggests a very anisotropic motion consisting of essentially free rotation about the symmetry axis and highly hindered rotation about the two perpendicular axes of the asymmetric water molecule. The rotational lifetime is significantly shortened relative to gas-phase water. An upper limit of 0.93 ps is deduced from the spectrum. Interaction with carbon tetrachloride also slightly enhances the intensity of the symmetric stretch. The results are compared with results of interactions between water and the cations Li+, Na+, K+, and Cs+. It is concluded that the attractive interaction is between the oxygen of water and the electropositive carbon of carbon tetrachloride.     

Chirality transfer through hydrogen-bonding: experimental and ab initio analyses of vibrational circular dichroism spectra of methyl lactate in water

The infrared vibrational absorption (VA) and vibrational circular dichroism (VCD) spectra of methyl lactate were measured in the 1000-1800 cm(-1) region in the CCl(4) and H(2)O solvents, respectively. In particular, the chirality transfer effect, i.e. the H-O-H bending bands of the achiral water subunits that are hydrogen-bonded to the methyl lactate molecule exhibit substantial VCD strength, was detected experimentally. A series of density functional theory calculations using B3PW91 and B3LYP functionals with 6-311++G(d,p) and aug-cc-pVTZ basis sets were carried out to simulate the VA and VCD spectra of the methyl lactate monomer and the methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3. The population weighted VA and VCD spectra of the methyl lactate monomer are in good agreement with the experimental spectra in CCl(4). Implicit polarizable continuum model was found to be inadequate to account for the hydrogen-bonding effect in the observed VA and VCD spectra in H(2)O. The methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3 were used to model the explicit hydrogen-bonding. The population weighted VA and VCD spectra of the methyl lactate-H(2)O binary complex are shown to capture the main spectral features in the observed spectra in aqueous solution. The theoretical modeling shows that the extent of chirality transfer depends sensitively on the specific binding sites taken by the achiral water molecules. The observation of chirality transfer effect opens a new spectral window to detect and to model the hydrogen-bonding solvent effect on VCD spectra of chiral molecules.     

Anomalous dielectric relaxation of water molecules at the surface of an aqueous micelle

Dielectric relaxation of aqueous solutions of micelles, proteins, and many complex systems shows an anomalous dispersion at frequencies intermediate between those corresponding to the rotational motion of bulk water and that of the organized assembly or macromolecule. The precise origin of this anomalous dispersion is not well-understood. In this work we employ large scale atomistic molecular dynamics simulations to investigate the dielectric relaxation (DR) of water molecules in an aqueous micellar solution of cesium pentadecafluorooctanoate. The simulations clearly show the presence of a slow component in the moment-moment time correlation function [PhiMW(t)] of water molecules, with a time constant of about 40 ps, in contrast to only 9 ps for bulk water. Interestingly, the orientational time correlation function [Cmu(t)] of individual water molecules at the surface exhibits a component with a time constant of about 19 ps. We show that these two time constants can be related by the well-known micro-macrorelations of statistical mechanics. In addition, the reorientation of surface water molecules exhibits a very slow component that decays with a time constant of about 500 ps. An analysis of hydrogen bond lifetime and of the rotational relaxation in the coordinate frame fixed on the micellar body seems to suggest that the 500 ps component owes its origin to the existence of an extended hydrogen bond network of water molecules at the surface. However, this ultraslow component is not found in the total moment-moment time correlation function of water molecules in the solution. The slow DR of hydration water is found to be well correlated with the slow solvation dynamics of cesium ions at the water-micelle interface.     

17O nuclear quadrupole coupling constants of water bound to a metal ion: a gadolinium(III) case study

Rotational correlation times of metal ion aqua complexes can be determined from 17O NMR relaxation rates if the quadrupole coupling constant of the bound water oxygen-17 nucleus is known. The rotational correlation time is an important parameter for the efficiency of Gd3+ complexes as magnetic resonance imaging contrast agents. Using a combination of density functional theory with classical and Car-Parrinello molecular dynamics simulations we performed a computational study of the 17O quadrupole coupling constants in model aqua ions and the [Gd(DOTA)(H2O)]- complex used in clinical diagnostics. For the inner sphere water molecule in the [Gd(DOTA)(H2O)]- complex the determined quadrupole coupling parameter chi square root of (1 + eta2/3) of 8.7 MHz is very similar to that of the liquid water (9.0 MHz). Very close values were also predicted for the the homoleptic aqua ions of Gd3+ and Ca2+. We conclude that the 17O quadrupole coupling parameters of water molecules coordinated to closed shell and lanthanide metal ions are similar to water molecules in the liquid state.     

Conformational landscape of (R,R)-pterocarpans with biological activity in vacuo and in aqueous solution (PCM and/or water clusters)

All possible combinations of stable dihedral values have been considered in vacuo at the B3LYP/6-31G level for 3,9-dihydroxy-4,8-diprenylpterocarpan (erybraedin C), whose hydroxy out-out conformation had been examined earlier together with the conformational preferences of 3,9-dimethoxy-4-prenylpterocarpan (bitucarpin A) at the same level (Phys. Chem. Chem. Phys. 2004, 6, 2849). The structure with O5 trans with respect to H6a (O(t)) is about 2 kcal/mol less stable in vacuo than that with one of the H6 trans to it (H(t)); in aqueous solution its energy gap is nearly conserved. The in-in arrangement of the hydroxyl groups of erybraedin turns out to be preferred in vacuo (even considering zero point and thermal effects), where pseudo H-bonds are formed between hydroxy hydrogens and pi electron distributions of prenyl groups. The continuum solvent effect (water) at the IEF-PCM/B3LYP/6-31G level on the relative stability of the various rotamers is very limited both on bitucarpin and erybraedin. Considering the dihydrated derivatives, significant differences in the solvation energy are found between the distinct hydration sites, increasing in the order: methoxy O, ring O, hydroxy O, and hydroxy H. In hydroxy-water interactions, in fact, water prefers to behave as an H-bond acceptor unless nearby bulky groups prevent its approach. Interestingly enough, a bridging water molecule between the hydroxy H of erybraedin and the prenyl group can be found. The inclusion of BSSE corrections in hydroxy-water interactions decidedly favors out-out hydrated arrangements, followed by out-in and in-out ones. Bulk solvent effects with IEF-PCM about the dihydrated systems almost invert the stability order found in vacuo. When a four-water cluster is considered using QM methods, waters gather in H-bonded pairs around the solute OH groups. MD simulations, carried out on a pterocarpan solute (J. Phys. Chem. B 2005, 109, 16918), supply water adducts consistent with a liquid state that have also been embedded in the continuum solvent.     

How does internal motion influence the relaxation of the water protons in Ln(III)DOTA-like complexes?

Taking advantage of the Curie contribution to the relaxation of the protons in the Tb(III) complex, and the quadrupolar relaxation of the 17O and 2H nuclei on the Eu(III) complex, the effect of the internal motion of the water molecule bound to [Ln(DOTAM)(H2O)]3+ complexes was quantified. The determination of the quadrupolar coupling constant of the bound water oxygen chi(Omicron)(1 + eta(Omicron)2/3)1/2 = 5.2 +/- 0.5 MHz allows a new analysis of the 17O and 1H NMR data of the [Gd(DOTA)(H2O)]- complex with different rotational correlation times for the Gd(III)-O(water) and Gd(III)-H(water) vectors. The ratio of the rotational correlation times for the Ln(III)-H(water) vector and the overall rotational correlation time is calculated tau(RH)/tau(RO) = 0.65 +/- 0.2. This could have negative consequences on the water proton relaxivity, which we discuss in particular for macromolecular systems. It appears that the final effect is actually attenuated and should be around 10% for such large systems undergoing local motion of the chelating groups.     

A coupled reference interaction site model/molecular dynamics study of the potential of mean force curve of the SN2 Cl- + CH3Cl reaction in water

An application of the coupled reference interaction site model (RISM)/simulation methodology to the calculation of the potential of mean force (PMF) curve in aqueous solution for the identity nucleophilic substitution reaction Cl(-) + CH(3)Cl is performed. The free energy of activation is calculated to be 27.1 kcal/mol which compares very well with the experimentally determined barrier height of 26.6 kcal/mol. Furthermore, the calculated PMF is almost superimposed with that previously calculated using the computationally rigorous Monte Carlo with importance sampling method (Chandrasekhar, J.; Smith, S. F.; Jorgensen, W. L. J. Am. Chem. Soc. 1985, 107, 154). Using the calculated PMF, a crude estimate of the solvated kinetic transmission coefficient also compares well with that of previous more accurate simulations. These results indicate that the coupled RISM/simulation method provides a cost-effective methodology for studying reactions in solution.     

Simple physics-based analytical formulas for the potentials of mean force for the interaction of amino acid side chains in water. IV. Pairs of different hydrophobic side chains

The potentials of mean force of 21 heterodimers of the molecules modeling hydrophobic amino acid side chains: ethane (for alanine), propane (for proline), isobutane (for valine), isopentane (for leucine and isoleucine), ethylbenzene (for phenylalanine), methyl propyl sulfide (for methionine), and indole (for tryptophane) were determined by umbrella-sampling molecular dynamics simulations in explicit water as functions of distance and orientation. Analytical expressions consisting of the Gay-Berne term to represent effective van der Waals interactions and the cavity term proposed in our earlier work were fitted to the potentials of mean force. The positions and depths of the contact minima and the positions and heights of the desolvation maxima, including their dependence on the orientation of the molecules, are well represented by the analytical expressions for all systems; large deviations between the MD-determined PMF and the analytical approximations are observed for pairs involving the least spheroidal solutes: ethylbenzene, indole, and methyl propyl sulfide at short distances at which the PMF is high and, consequently, these regions are rarely visited. When data from the PMF within only 10 kcal/mol above the global minimum are considered, the standard deviation between the MD-determined and the fitted PMF is from 0.25 to 0.55 kcal/mol (the relative standard deviation being from 4% to 8%); it is larger for pairs involving nonspherical solute molecules. The free energies of contact computed from the PMF surfaces are well correlated with those determined from protein-crystal data with a slope close to that relating the free energies of transfer of amino acids (from water to n-octanol) to the average contact free energies determined from protein-crystal data. These observations justify future use of the determined potentials in coarse-grained protein-folding simulations.     

Hydrogen bonding and cation coordination effects in primary and secondary amines dissolved in carbon tetrachloride

Raman and infrared spectroscopy were used to investigate hydrogen-bonding interactions and cation coordination effects in solutions of lithium triflate (LiCF3SO3) dissolved in two primary amines, hexylamine (HEXA) and N,N-dimethylethylenediamine (DMEDA), and in a secondary amine, dipropylamine (DPA). Strong intermolecular hydrogen-bonding interactions and weaker intramolecular hydrogen-bonding interactions that occur only in DMEDA were spectroscopically distinguished in a comparison of pure HEXA, pure DMEDA, and the dilute solutions of these amines in CCl4. The spectroscopic shifts in intensity and frequency in the NH stretching region of DPA and DPA diluted in CCl4 were similar to those of HEXA. Dilute electrolyte solutions in carbon tetrachloride were prepared to analyze specifically the cation coordination effect. In these solutions, limited intermolecular hydrogen-bonding interactions are present, and the observed spectral shifts correspond primarily to the cation-induced shifts. The symmetric SO3 stretching region of the triflate anion was investigated to probe further the coordination of the cation. The local structures of the triflate ions and the amine groups in the electrolyte solutions dissolved in CCl4 are similar to the local structures in the corresponding amine-salt crystals previously reported by us.     

Molecular reorientation of liquid water studied with femtosecond midinfrared spectroscopy

The molecular reorientation of liquid water is key to the hydration and stabilization of molecules and ions in aqueous solution. A powerful technique to study this reorientation is to measure the time-dependent anisotropy of the excitation of the O-H/O-D stretch vibration of HDO dissolved in D2O/H2O using femtosecond midinfrared laser pulses. In this paper, we present and discuss experiments in which this technique is used to study the correlation between the molecular reorientation of the water molecules and the strength of the hydrogen-bond interactions. On short time scales (<200 fs), it was found that the anisotropy shows a partial decay due to librational motions of the water molecules that keep the hydrogen bond intact. On longer time scale (>200 fs), the anisotropy shows a complete decay with an average time constant of 2.5 ps. From the frequency dependence of the anisotropy dynamics, it follows that a subensemble of the water molecules shows a fast reorientation that is accompanied by a large change of the vibrational frequency. This finding agrees with the molecular jumping mechanism for the reorientation of liquid water that has recently been proposed by Laage and Hynes.