Computational studies of liquid water and diluted water in carbon tetrachloride

Molecular dynamics simulations were carried out to study solvent effects on the energetic and dynamical properties of water molecules in liquid water and in carbon tetrachloride (CCl4). In these studies, the free-energy profiles or potentials of mean force (PMF) for water dimers in both solvents were computed. The computed PMF results showed a stable minimum near 3 A for the O-O separation, with a minimum free energy of about -2.8 kcal/mol in CCl4, as compared to a value of -0.5 kcal/mol in liquid water. The difference in free energy in water as compared to that in CCl4 was expected and is the result of competition from surrounding water molecules that are capable of forming hydrogen bonds in the liquid water. This capability is absent in the diluted water found in CCl4. We found that the rotational motions of H2O/D2O were nonisotropic, with the out-of-plane vector correlation times in H2O/D2O varying from 5.6/5.8 ps at 250 K to 0.57/0.56 ps at 350 K and the corresponding OH/OD bond vectors varying from 6.5/7.7 ps to 0.75/0.75 ps. The results compare reasonably well to the available NMR experimental and computer simulation data on the same system (Farrar; Skinner; et al. J. Am. Chem. Soc. 2001, 123, 8047). For diluted water in CCl4, we found the computed rotational correlation times also were nonisotropic and much longer than the corresponding NMR experimental values at the same concentration (Farrar; et al. J. Phys. Chem. A 2007, 111, 6146). Upon analyzing the water hydrogen-bonding patterns as a function of water concentration, we conclude that the differences in the rotational correlation times mainly result from the formation of water hydrogen-bonding networks as the water concentration is increased in liquid CCl4. In addition, we found the rotational correlation times to be substantially faster in liquid CCl4 than in liquid water.     

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.     

Physicochemical characterization of the dimeric lanthanide complexes [en{Ln(DO3A)(H2O)}2] and [pi{Ln(DTTA)(H2O)}2]2-: a variable-temperature 17O NMR study

The Gd(III) complexes of the two dimeric ligands [en(DO3A)2] {N,N'-bis[1,4,7-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecan-10-yl-methylcarbonyl]-N,N'-ethylenediamine} and [pi(DTTA)2]8- [bisdiethylenetriaminepentaacetic acid (trans-1,2-cyclohexanediamine)] were synthesized and characterized. The 17O NMR chemical shift of H2O induced by [en{Dy(DO3A)}2] and [pi{Dy(DTTA)}2]2- at pH 6.80 proved the presence of 2.1 and 2.2 inner-sphere water molecules, respectively. Water proton spin-lattice relaxation rates for [en{Gd(DO3A)(H2O)}2] and [pi{Gd(DTTA)(H2O)}2]2- at 37.0 +/- 0.1 degrees C and 20 MHz are 3.60 +/- 0.05 and 5.25 +/- 0.05 mM(-1) s(-1) per Gd, respectively. The EPR transverse electronic relaxation rate and 17O NMR transverse relaxation time for the exchange lifetime of the coordinated H2O molecule and the 2H NMR longitudinal relaxation rate of the deuterated diamagnetic lanthanum complex for the rotational correlation time were thoroughly investigated, and the results were compared with those reported previously for other lanthanide(III) complexes. The exchange lifetimes for [en{Gd(DO3A)(H2O)}2] (769 +/- 10 ns) and [pi{Gd(DTTA)(H2O)}2]2- (910 +/- 10 ns) are significantly higher than those of [Gd(DOTA)(H2O)]- (243 ns) and [Gd(DTPA)(H2O)]2- (303 ns) complexes. The rotational correlation times for [en{Gd(DO3A)(H2O)}2] (150 +/- 11 ps) and [pi{Gd(DTTA)(H2O)}2]2- (130 +/- 12 ps) are slightly greater than those of [Gd(DOTA)(H2O)]- (77 ps) and [Gd(DTPA)(H2O)]2- (58 ps) complexes. The marked increase in relaxivity (r1) of [en{Gd(DO3A)(H2O)}2] and [pi{Gd(DTTA)(H2O)}2]2- result mainly from their longer rotational correlation time and higher molecular weight.     

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.     

Rotational dynamics in supercooled water from nuclear spin relaxation and molecular simulations

Structural dynamics in liquid water slow down dramatically in the supercooled regime. To shed further light on the origin of this super-Arrhenius temperature dependence, we report high-precision (17)O and (2)H NMR relaxation data for H(2)O and D(2)O, respectively, down to 37 K below the equilibrium freezing point. With the aid of molecular dynamics (MD) simulations, we provide a detailed analysis of the rotational motions probed by the NMR experiments. The NMR-derived rotational correlation time τ(R) is the integral of a time correlation function (TCF) that, after a subpicosecond librational decay, can be described as a sum of two exponentials. Using a coarse-graining algorithm to map the MD trajectory on a continuous-time random walk (CTRW) in angular space, we show that the slowest TCF component can be attributed to large-angle molecular jumps. The mean jump angle is ～48° at all temperatures and the waiting time distribution is non-exponential, implying dynamical heterogeneity. We have previously used an analogous CTRW model to analyze quasielastic neutron scattering data from supercooled water. Although the translational and rotational waiting times are of similar magnitude, most translational jumps are not synchronized with a rotational jump of the same molecule. The rotational waiting time has a stronger temperature dependence than the translation one, consistent with the strong increase of the experimentally derived product τ(R)?D(T) at low temperatures. The present CTRW jump model is related to, but differs in essential ways from the extended jump model proposed by Laage and co-workers. Our analysis traces the super-Arrhenius temperature dependence of τ(R) to the rotational waiting time. We present arguments against interpreting this temperature dependence in terms of mode-coupling theory or in terms of mixture models of water structure.     

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.     

Unexpected aggregation of neutral, xylene-cored dinuclear GdIII chelates in aqueous solution

We have synthesized ditopic ligands L(1), L(2), and L(3) that contain two DO3A(3-) metal-chelating units with a xylene core as a noncoordinating linker (DO3A(3-) = 1,4,7,10-tetraazacyclododecane-1,4,7-triacetate; L(1) = 1,4-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzene; L(2) = 1,3-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzene; L(3) = 3,5-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzoic acid). Aqueous solutions of the dinuclear Gd(III) complexes formed with the three ligands have been investigated in a variable-temperature, multiple-field (17)O NMR and (1)H relaxivity study. The (17)O longitudinal relaxation rates measured for the [Gd(2)L(1-3)(H2O)(2)] complexes show strong field dependence (2.35-9.4 T), which unambiguously proves the presence of slowly tumbling entities in solution. The proton relaxivities of the complexes, which are unexpectedly high for their molecular weight, and in particular the relaxivity peaks observed at 40-50 MHz also constitute experimental evidences of slow rotational motion. This was explained in terms of self-aggregation related to hydrophobic interactions, pi stacking between the aromatic linkers, or possible hydrogen bonding between the chelates. The longitudinal (17)O relaxation rates of the [Gd(2)L(1-3)(H2O)(2)] complexes have been analysed with the Lipari-Szabo approach, leading to local rotational correlation times tau(1)(298) of 150-250 ps and global rotational correlation times tau(g)(298) of 1.6-3.4 ns (c(Gd): 20-50 mM), where tau(1)(298) is attributed to local motions of the Gd segments, while tau(g)(298) describes the overall motion of the aggregates. The aggregates can be partially disrupted by phosphate addition; however, at high concentrations phosphate interferes in the first coordination sphere by replacing the coordinated water. In contrast to the parent [Gd(DO3A)(H2O)(1.9)], which presents a hydration equilibrium between mono- and dihydrated species, a hydration number of q = 1 was established for the [Ln(2)L(1-3)(H2O)(2)] chelates by (17)O chemical shift measurements on Ln = Gd and UV/Vis spectrophotometry for Ln = Eu. The exchange rate of the coordinated water is higher for [Gd(2)L(1-3)(H2O)(2)] complexes k(ex)(298) = 7.5-12.0 x 10(6) s(-1)) than for [Gd(DOTA)(H2O)](-). The proton relaxivity of the [Gd(2)L(1-3)(H2O)(2)] complexes strongly decreases with increasing pH. This is related to the deprotonation of the inner-sphere water, which has also been characterized by pH potentiometry. The protonation constants determined for this process are logK(OH) = 9.50 and 10.37 for [Gd(2)L(1)(H2O)(2)] and [Gd(2)L(3)(H2O)(2)], respectively.     

The dynamical properties of the aromatic hydrogen bond in NH4(C6H5)4B from quasielastic neutron scattering

NH(4)(C(6)H(5))(4)B represents a prototypical system for understanding aromatic H bonds. In NH(4)(C(6)H(5))(4)B an ammonium cation is trapped in an aromatic cage of four phenyl rings and each phenyl ring serves as a hydrogen bond acceptor for the ammonium ion as donor. Here the dynamical properties of the aromatic hydrogen bond in NH(4)(C(6)H(5))(4)B were studied by quasielastic incoherent neutron scattering in a broad temperature range (20< or =T< or =350 K). We show that in the temperature range from 67 to 350 K the ammonium ions perform rotational jumps around C(3) axes. The correlation time for this motion is the lifetime of the "transient" H bonds. It varies from 1.5 ps at T=350 K to 150 ps at T=67 K. The activation energy was found to be 3.14 kJ mol, which means only 1.05 kJ mol per single H bond for reorientations around the C(3) symmetry axis of the ammonium group. This result shows that the ammonium ions have to overcome an exceptionally low barrier to rotate and thereby break their H bonds. In addition, at temperatures above 200 K local diffusive reorientational motions of the phenyl rings, probably caused by interaction with ammonium-group reorientations, were found within the experimental observation time window. At room temperature a reorientation angle of 8.4 degrees +/-2 degrees and a correlation time of 22+/-8 ps were determined for the latter. The aromatic H bonds are extremely short lived due to the low potential barriers allowing for molecular motions with a reorientational character of the donors. The alternating rupture and formation of H bonds causes very strong damping of the librational motion of the acceptors, making the transient H bond appear rather flexible.     

Synthesis and characterization of the novel monoamide derivatives of Gd-TTDA

For this study, the N'-monoamide derivatives of TTDA (3,6,10-tri(carboxymethyl)-3,6,10-triazadodecanedioic acid), N'-methylamide (TTDA-MA), N'-benzylamide (TTDA-BA), and N'-2-methoxybenzylamide (TTDA-MOBA), were synthesized. Their protonation constants and stability constants (log K(ML)'s) formed with Ca(2+), Zn(2+), Cu(2+), and Gd(3+) were determined by potentiometric titration in 0.10 M Me(4)NCl at 25.0 +/- 0.1 degrees C. The relaxivity values of [Gd(TTDA-MA)](-), [Gd(TTDA-BA)](-), and [Gd(TTDA-MOBA)](-) remained constant with respect to pH changes over the range 4.5-12.0. The (17)O NMR chemical shift of H(2)O induced by [Dy(TTDA-MA)(H(2)O)](-) at pH 6.80 showed 0.9 inner-sphere water molecules. Water proton relaxivity values for [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) at 37.0 +/- 0.1 degrees C and 20 MHz are 3.89, 4.21, and 4.25, respectively. The water-exchange lifetime (tau(M)) and rotational correlation time (tau(R)) of [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) are obtained from reduced the (17)O relaxation rate and chemical shifts of H(2)(17)O. The (2)H NMR longitudinal relaxation rates of the deuterated diamagnetic lanthanum complexes for the rotational correlation time were also thoroughly investigated. The water-exchange rates (K(298)(ex) for [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) are lower than that of [Gd(TTDA)(H(2)O)](2)(-) but significantly higher than those of [Gd(DTPA)(H(2)O)](2)(-) and [Gd(DTPA-BMA)(H(2)O)]. The rotational correlation times for [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) are significantly longer than those of [Gd(TTDA)(H(2)O)](2)(-) and [Gd(DTPA)(H(2)O)](2)(-) complexes. The marked increase of the relaxivity of [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) results mainly from their longer rotational correlation time. The noncovalent interaction between human serum albumin (HSA) and [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) complexes containing a hydrophobic substituent was investigated by measuring the water proton relaxation rate of the aqueous solutions. The binding association constant (K(A)) values are 1.0 +/- 0.2 x 10(3) and 1.3 +/- 0.2 x 10(3) M(-1) for [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-), which indicates a stronger interaction of [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) with HSA.     

A benzene-core trinuclear GdIII complex: towards the optimization of relaxivity for MRI contrast agent applications at high magnetic field

A novel ligand, H(12)L, based on a trimethylbenzene core bearing three methylenediethylenetriamine-N,N,N',N'-tetraacetate moieties (-CH(2)DTTA(4-)) for Gd(3+) chelation has been synthesized, and its trinuclear Gd(3+) complex [Gd(3)L(H(2)O)(6)](3-) investigated with respect to MRI contrast agent applications. A multiple-field, variable-temperature (17)O NMR and proton relaxivity study on [Gd(3)L(H(2)O)(6)](3-) yielded the parameters characterizing water exchange and rotational dynamics. On the basis of the (17)O chemical shifts, bishydration of Gd(3+) could be evidenced. The water exchange rate, k(ex)(298)=9.0+/-3.0 s(-1) is around twice as high as k(ex)(298) of the commercial [Gd(DTPA)(H(2)O)](2-) and comparable to those on analogous Gd(3+)-DTTA chelates. Despite the relatively small size of the complex, the rotational dynamics had to be described with the Lipari-Szabo approach, by separating global and local motions. The difference between the local and global rotational correlation times, tau(lO)(298)=170+/-10 ps and tau(gO)(298)=540+/-100 ps respectively, shows that [Gd(3)L(H(2)O)(6)](3-) is not fully rigid; its flexibility originates from the CH(2) linker between the benzene core and the poly(amino carboxylate) moiety. As a consequence of the two inner-sphere water molecules per Gd(3+), their close to optimal exchange rate and the appropriate size and limited flexibility of the molecule, [Gd(3)L(H(2)O)(6)](3-) has remarkable proton relaxivities when compared with commercial contrast agents, particularly at high magnetic fields (r(1)=21.6, 17.0 and 10.7 mM(-1)s(-1) at 60, 200 and 400 MHz respectively, at 25 degrees C; r(1) is the paramagnetic enhancement of the longitudinal water proton relaxation rate, referred to 1 mM concentration of Gd(3+)).     

Water Exchange and Rotational Dynamics of the Dimeric Gadolinium(III) Complex [BO{Gd(DO3A)(H(2)O)}(2)]: A Variable-Temperature and -Pressure (17)O NMR Study(1)

Rapid water exchange and slow rotation are essential for high relaxivity MRI contrast agents. A variable-temperature and -pressure (17)O NMR study at 14.1, 9.4, and 1.4 T has been performed on the dimeric BO(DO3A)(2), 2,11-dihydroxy-4,9-dioxa-1,12-bis[1,4,7,10-tetraaza-4,7,10-tris(carboxymethyl)cyclododecyl]dodecane, complex of Gd(III). This complex is of relevance to MRI as an attempt to gain higher (1)H relaxivity by slowing down the rotation of the molecule compared to monomeric Gd(III) complexes used as contrast agents. From the (17)O NMR longitudinal and transverse relaxation rates and chemical shifts we determined the parameters characterizing water exchange kinetics and the rotational motion of the complex, both of which influence (1)H relaxivity. The rate constant and the activation enthalpy for the water exchange, k(ex) and DeltaH(), are (1.0 +/- 0.1) x 10(6) s(-)(1)and (30.0 +/- 0.2) kJ mol(-)(1), respectively, and the activation volume, DeltaV(), of the process is (+0.5 +/- 0.2) cm(3) mol(-)(1), indicating an interchange mechanism. The rotational correlation time becomes about three times longer compared to monomeric Gd(III) polyamino-polyacetate complexes studied so far: tau(R) = (250 +/- 5) ps, which results in an enhanced proton relaxivity by raising the correlation time for the paramagnetic interaction.     

Molecular dynamics study of translational and rotational diffusion in liquid ortho-terphenyl

NVT molecular dynamics simulations were performed on liquid o-terphenyl as a function of temperature in the range 320-480 K. Computed translational diffusion coefficients displayed the non-Arrhenius behavior expected of a fragile glass-forming liquid and were in good, semiquantitative agreement with experimental results. Rotational correlation functions calculated for various vectors within the molecule exhibited a very short time (0-1 ps) initial decay, followed by a reversal, which corresponds to free reorientation within the "solvent" cage prior to collision with a wall. Rotational correlation times of three orthogonal vectors fixed on the central benzene were close to equal at all temperatures, indicating nearly isotropic overall molecular reorientation. The average correlation times exhibited a non-Arrhenius temperature dependence and were in very good agreement with experimental values derived from 2D and 1H NMR relaxation times. Correlation times of vectors located on the lateral phenyl rings were used to calculate the "spinning" internal rotation diffusion coefficients, which were approximately twice as great as the overall rotational diffusion constants, indicating rapid internal rotation of the phenyl side groups over wide ranges of angle in the liquid.     

Dielectric behavior of some small ketones as ideal polar molecules

The dielectric behaviors of some small symmetric ketone molecules, including acetone, 3-pentanone, cyclopentanone, 4-heptanone, and cyclohexanone, were investigated as a function of temperature (T) over a wide frequency range from 50 MHz (3.14 × 10(8) s(-1), in angular frequency) to 3 THz (1.88 × 10(13) s(-1)). The temperature dependencies of the rotational diffusion times (τ(r)) determined using (17)O NMR spin-lattice relaxation time (T(1)) measurements and viscosities of the ketones were also examined. The obtained temperature dependencies of the parameters for the ketones were compared with those of ideal polar molecules, which obey the Stokes-Einstein-Debye (SED) relationship without the formation of intermolecular dimeric associations and without orientational correlations between dipoles (molecular axes), that is, free rotation. Kirkwood correlation factors (g(K)) of only acetone and 3-pentanone were close to unity over a wide temperature range, whereas those of other ketones were obviously less than unity. These results revealed that no correlations exist between the rotational motions of dipoles in acetone and 3-pentanone, as expected in ideal polar molecules. However, other ketones exhibited orientational correlations in their dipoles because of dipole-dipole interactions via antiparallel configurations. Furthermore, because acetone and 3-pentanone satisfied the SED relationship and because their microscopic dielectric relaxation times (τ(μ)), which were calculated from the determined dielectric relaxation times (τ(D)) via the relationship τ(μ) = τ(D)g(K)(-1), were identical to 3τ(r) and were proportional to Vη(k(B)T)(-1) over the wide temperature range examined, where V, k(B), and η represent the effective molecular volume, Boltzmann's constant, and the viscosity of the liquid molecules, respectively, these two ketone molecules behave as ideal polar molecules. In addition, other ketones not significantly larger than acetone and 3-pentanone in molecular size likely form dimeric intermolecular associations with antiparallel cyclic configurations, which demonstrates the g(K) values less than unity.     

PICOSECOND RESOLUTION OF INDOLE ANISOTROPY DECAYS AND SPECTRAL RELAXATION BY 2 GHz FREQUENCY-DOMAIN FLUOROMETRY

Abstract— We used frequency-domain fluorescence spectroscopy to measure rotational diffusion and time-resolved emission spectra of indole in methanol on the picosecond timescale. The indole emission was quenched by acrylamide to allow measurements to the instrumental limit of 2 GHz and to eliminate emission from the longer-lived indole molecules, which can no longer provide information on the picosecond (ps) processes. The resolution was adequate to measure rotational correlation times as short as 8 ps at 80†C, and spectral relaxation times as short as 16 ps at 5†C.     

Rotational dynamics account for pH-dependent relaxivities of PAMAM dendrimeric, Gd-based potential MRI contrast agents

The EPTPA5) chelate, which ensures fast water exchange in GdIII complexes, has been coupled to three different generations (5, 7, and 9) of polyamidoamine (PAMAM) dendrimers through benzylthiourea linkages (H5EPTPA = ethylenepropylenetriamine-N,N,N',N',N'-pentaacetic acid). The proton relaxivities measured at pH 7.4 for the dendrimer complexes G5-(GdEPTPA)111, G7-(GdEPTPA)253 and G9-(GdEPTPA)1157 decrease with increasing temperature, indicating that, for the first time for dendrimers, slow water exchange does not limit relaxivity. At a given field and temperature, the relaxivity increases from G5 to G7, and then slightly decreases for G9 (r1 = 20.5, 28.3 and 27.9 mM(-1) s(-1), respectively, at 37 degrees C, 30 MHz). The relaxivities show a strong and reversible pH dependency for all three dendrimer complexes. This originates from the pH-dependent rotational dynamics of the dendrimer skeleton, which was evidenced by a combined variable-temperature and multiple-field 17O NMR and 1H relaxivity study performed at pH 6.0 and 9.9 on G5-(GdEPTPA)111. The longitudinal 17O and 1H relaxation rates of the dendrimeric complex are strongly pH-dependent, whereas they are not for the [Gd(EPTPA)(H2O)]2- monomer chelate. The longitudinal 17O and 1H relaxation rates have been analysed by the Lipari-Szabo spectral density functions and correlation times have been calculated for the global motion of the entire macromolecule (tau(gO)) and the local motion of the GdIII chelates on the surface (tau(lO)), correlated by means of an order parameter S2. The dendrimer complex G5-(GdEPTPA)111 has a considerably higher tau(gO) under acidic than under basic conditions (tau(298)gO = 4040 ps and 2950 ps, respectively), while local motions are less influenced by pH (tau(298)lO = 150 and 125 ps). The order parameter, characterizing the rigidity of the macromolecule, is also higher at pH 6.0 than at pH 9.9 (S2 = 0.43 vs 0.36, respectively). The pH dependence of the global correlation time can be related to the protonation of the tertiary amine groups in the PAMAM skeleton, which leads to an expanded and more rigid dendrimeric structure at lower pH. The increase of tau(gO) with decreasing pH is responsible for the pH dependent proton relaxivities. The water exchange rate on G5-(GdEPTPA)111(k(298)ex = 150 x 10(6) s(-1)) shows no significant pH dependency and is similar to the one measured for the monomer [Gd(EPTPA)(H2O)]2-. The proton relaxivity of G5-(GdEPTPA)111 is mainly limited by the important flexibility of the dendrimer structure, and to a small extent, by a faster than optimal water exchange rate.     

Thermodynamic, structural, and dynamic properties of supercooled water confined in mesoporous MCM-41 studied with calorimetric, neutron diffraction, and neutron spin echo measurements

Thermodynamic, structural, and dynamic properties of heavy water (D(2)O) confined in mesoporous silica glass MCM-41 C10, C12, and C14 were investigated by differential scanning calorimetry, neutron diffraction, and neutron spin echo (NSE) measurements, respectively. The DSC data showed that no crystallization of D(2)O confined in C10 occurs in a temperature range between 298 and 180 K, and that crystalline ice is formed at 204 and 221 K for C12 and C14, respectively. For C10, the neutron radial distribution functions of confined D(2)O suggested a structural change in the supercooled state between 223 and 173 K. For C10 sample, it has been found that the tetrahedral-like water structure is partially enhanced in the central part of pores at 173 K. For all the samples, the intermediate scattering functions from the NSE measurements are fitted by the Kohlrausch-Williams-Watts stretched exponential function which implies that confined supercooled D(2)O exhibits a wide distribution of relaxation times. For C10, C12, and C14 samples, between 298 and 240 K, the relaxation times of supercooled D(2)O follow remarkably well the Vogel-Fulcher-Tamman equation; for C10 sample, below 240 K, the relaxation times of nonfreezing D(2)O show an Arrhenius type behavior. From the present experimental results on calorimetric, structural, and dynamic properties, it has been concluded that supercooled D(2)O confined in MCM-41 C10 experiences a transition from high-density to low-density hydrogen-bonded structure at around 229 K.     

Relaxation dynamics in (HF)x(H2O)1-x solutions

The high-frequency dynamics of (HF)(x)(H(2)O)(1-x) solutions has been investigated by inelastic x-ray scattering. The measurements have been performed as a function of the concentration in the range x = 0.20-0.73 at fixed temperature T = 283 K. The results have been compared with similar data in pure water (x = 0) and pure hydrogen fluoride (x = 1). A viscoelastic analysis of the data highlights the presence of a relaxation process characterized by a relaxation time and a strength directly related to the presence of a hydrogen-bond network in the system. The comparison with the data on water and hydrogen fluoride shows that the structural relaxation time continuously decreases at increasing concentration of hydrogen fluoride passing from the value for water to the one for hydrogen fluoride tau(alphaHF), which is three times smaller. This is the consequence of a gradual decreasing number of constraints of the hydrogen-bond networks in passing from one liquid to the other.     

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.     

Towards Binuclear Polyaminocarboxylate MRI Contrast Agents? Spectroscopic and MD Study of the Peculiar Aqueous Behavior of the LnIII Chelates of OHEC (Ln=Eu,Gd, and Tb): Implications for Relaxivity

We report the study of binuclear Ln(III) chelates of OHEC (OHEC=octaazacyclohexacosane-1,4,7,10,14,17,20,23-octaacetate). The interconversion between two isomeric forms, which occurs in aqueous solution, has been studied by NMR, UV/Vis, EPR, and luminescence spectroscopy, as well as by classical molecular dynamics (MD) simulations. For the first time we have characterized an isomerization equilibrium for a Ln(III) polyaminocarboxylate complex (Ln(III)=Y, Eu, Gd and Tb) in which the metal centre changes its coordination number from nine to eight, such that: [Ln(2)(ohec)(H(2)O)(2)](2-) r<==>[Ln(2)(ohec)](2-)+2 H(2)O. The variable temperature and pressure NMR measurements conducted on this isomerization reaction give the following thermodynamic parameters for Eu(III): K(298)=0.42+/-0.01, DeltaH(0)=+4.0+/-0.2 kJ mol(-1), DeltaS(0)=+6.1+/-0.5 J K(-1) mol(-1) and DeltaV(0)=+3.2+/-0.2 cm(3) mol(-1). The isomerization is slow and the corresponding kinetic parameters obtained by NMR spectroscopy are: k(298)(is)=73.0+/-0.5 s(-1), DeltaH++(is)=75.3+/-1.9 kJ mol(-1), DeltaS++(is)= +43.1+/-5.8 J K(-1) mol(-1) and DeltaV++(is)=+7.9+/-0.7 cm(3) mol(-1). Variable temperature and pressure (17)O NMR studies have shown that water exchange in [Gd(2)(ohec)(H(2)O)(2)](2-) is slow, k(298)(ex)=(0.40+/-0.02)x10(6) s(-1), and that it proceeds through a dissociative interchange I(d) mechanism, DeltaV( not equal )=+7.3+/-0.3 cm(3) mol(-1). The anisotropy of this oblong binuclear complex has been highlighted by MD simulation calculations of different rotational correlation times. The rotational correlation time directed on the Gd-Gd axis is 24 % longer than those based on the axes orthogonal to the Gd-Gd axis. The relaxivity of this binuclear complex has been found to be low, since 1) only [Gd(2)(ohec)(H(2)O)(2)](2-), which constitutes 70 % of the binuclear complex, contributes to the inner-sphere relaxivity and 2) the anisotropy of the complex prevents water molecules from having complete access to both Gd(III) cages; this decreases the outer-sphere relaxivity. Moreover, EPR measurements for the Gd(III) and for the mixed Gd(III)/Y(III) binuclear complexes have clearly shown that the two Gd(III) centres interact intramolecularly; this enhances the electronic relaxation of the Gd(III) electron spins.