Gadolinium (III) ion in liquid water: Structure, dynamics, and magnetic interactions from first principles

Effect of hydration equilibria on the relaxometric properties of Gd(III) complexes: new insights into old systems

We present a detailed relaxometric and computational investigation of three Gd(III) complexes that exist in solution as an equilibrium of two species with a different number of coordinated water molecules: [Gd(H2O)q]3+ (q = 8, 9), [Gd(EDTA)(H2O)q]- and [Gd(CDTA)(H2O)q]- (q = 2, 3). 1H nuclear magnetic relaxation dispersion (NMRD) data were recorded from aqueous solutions of these complexes using a wide Larmor frequency range (0.01-500 MHz). These data were complemented with 17O transverse relaxation rates and chemical shifts recorded at different temperatures. The simultaneous fit of the NMRD and 17O NMR data was guided by computational studies performed at the DFT and CASSCF/NEVPT2 levels, which provided information on Gd⋯H distances, 17O hyperfine coupling constants and the zero-field splitting (ZFS) energy, which affects electronic relaxation. The hydration equilibrium did not have a very important effect in the fits of the experimental data for [Gd(H2O)q]3+ and [Gd(CDTA)(H2O)q]-, as the hydration equilibrium is largely shifted to the species with the lowest hydration number (q = 8 and 2, respectively). The quality of the analysis improves however considerably for [Gd(EDTA)(H2O)q]- upon considering the effect of the hydration equilibrium. As a result, this study provides for the first time an analysis of the relaxation properties of this important model system, as well as accurate parameters for [Gd(H2O)q]3+ and [Gd(CDTA)(H2O)q]-.

Ionic Contraction across the Lanthanide Series Decreases the Temperature-Induced Disorder of the Water Coordination Sphere

In liquid, temperature affects the structures of lanthanide complexes in multiple ways that depend upon complex interactions between ligands, anions, and solvent molecules. The relative simplicity of lanthanide aqua ions (Ln³⁺) make them well suited to determine how temperature induces structural changes in lanthanide complexes. We performed a combination of ab initio molecular dynamics (AIMD) simulations and extended X-ray absorption fine structure (EXAFS) measurements, both at 25 and 90 °C, to determine how temperature affects the first- and second-coordination spheres of three Ln³⁺ (Ce³⁺, Sm³⁺, and Lu³⁺) aqua ions. AIMD simulations show first lanthanide coordination spheres that are similar at 25 and 90 °C, more so for the Lu³⁺ ion that remains as eight-coordinate than for the Ce³⁺ and Sm³⁺ ions that change their preferred coordination number from nine (at 25 °C) to eight (at 90 °C). The measured EXAFS spectra are very similar at 25 and 90 °C, for the Ce³⁺, Sm³⁺, and Lu³⁺ ions, suggesting that the dynamical disorder of the Ln³⁺ ions in liquid water is sufficient such that temperature-induced changes do not clearly manifest changes in the structure of the three ions. Both AIMD simulations and EXAFS measurements show very similar structures of the first coordination sphere of the Lu³⁺ ion at 25 and 90 °C.

Predicting 1H NMR relaxation in Gd3+-aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict ¹H NMR T₁ relaxation of water from paramagnetic Gd³⁺ ions in solution at 25 °C. Simulations of the T₁ relaxivity dispersion function r₁ computed from the Gd³⁺-¹H dipole-dipole autocorrelation function agree within ≃8% of measurements in the range f₀ ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. The simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f₀ ≲ 5 MHz, the simulation overestimates r₁ compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r₁ at high frequencies in chelated Gd³⁺ contrast-agents used for clinical MRI.

Dynamic Quantum Sensing of Paramagnetic Species Using Nitrogen-Vacancy Centers in Diamond

Naturally occurring paramagnetic species (PS), such as free radicals and paramagnetic metalloproteins, play an essential role in a multitude of critical physiological processes including metabolism, cell signaling, and immune response. These highly dynamic species can also act as intrinsic biomarkers for a variety of disease states, while synthetic paramagnetic probes targeted to specific sites on biomolecules enable the study of functional information such as tissue oxygenation and redox status in living systems. The work presented herein describes a new sensing method that exploits the spin-dependent emission of photoluminescence (PL) from an ensemble of nitrogen-vacancy centers in diamond for rapid, nondestructive detection of PS in living systems. Uniquely this approach involves simple measurement protocols that assess PL contrast with and without the application of microwaves. The method is demonstrated to detect concentrations of paramagnetic salts in solution and the widely used magnetic resonance imaging contrast agent gadobutrol with a limit of detection of less than 10 attomol over a 100 μm × 100 μm field of view. Real-time monitoring of changes in the concentration of paramagnetic salts is demonstrated with image exposure times of 20 ms. Further, dynamic tracking of chemical reactions is demonstrated via the conversion of low-spin cyanide-coordinated Fe3+ to hexaaqua Fe3+ under acidic conditions. Finally, the capability to map paramagnetic species in model cells with subcellular resolution is demonstrated using lipid membranes containing gadolinium-labeled phospholipids under ambient conditions in the order of minutes. Overall, this work introduces a new sensing approach for the realization of fast, sensitive imaging of PS in a widefield format that is readily deployable in biomedical settings. Ultimately, this new approach to nitrogen vacancy-based quantum sensing paves the way toward minimally invasive real-time mapping and observation of free radicals in in vitro cellular environments.

Ion Association in Lanthanide Chloride Solutions

A better understanding of the solution chemistry of the lanthanide (Ln) salts in water would have wide ranging implications in materials processing, waste management, element tracing, medicine and many more fields. This is particularly true for minerals processing, given governmental concerns about lanthanide security of supply and the drive to identify environmentally sustainable processing routes. Despite much effort, even in simple systems, the mechanisms and thermodynamics of LnIII association with small anions remain unclear. In the present study, molecular dynamics (MD), using a newly developed force field, provide new insights into LnCl3 (aq) solutions. The force field accurately reproduces the structure and dynamics of Nd3+ , Gd3+ and Er3+ in water when compared to calculations using density functional theory (DFT). Adaptive-bias MD simulations show that the mechanisms for ion pairing change from dissociative to associative exchange depending upon cation size. Thermodynamics of association reveal that whereas ion pairing is favourable, the equilibrium distribution of species at low concentration is dominated by weakly bound solvent-shared and solvent-separated ion pairs, rather than contact ion pairs, reconciling a number of contrasting observations of LnIII -Cl association in the literature. In addition, we show that the thermodynamic stabilities of a range of inner sphere and outer sphere LnCl x ( 3 - x ) + coordination complexes are comparable and that the kinetics of anion binding to cations may control solution speciation distributions beyond ion pairs. The techniques adopted in this work provide a framework with which to investigate more complex solution chemistries of cations in water.

Extraction of Gd3+ and UO22+ Ions Using Polystyrene Grafted Dibenzo Crown Ether (DB18C6) with Octanol and Nitrobenzene: A Molecular Dynamics Study

Atomistic molecular dynamics (MD) simulations are performed in order to derive thermodynamic properties important to understand the extraction of gadolinium (Gd³⁺) and uranium dioxide (UO₂) with dibenzo crown ether (DBCE) in nitrobenzene (NB) and octanol (OCT) solvents. The effect of polystyrene graft length, on DBCE, on the binding behavior of Gd³⁺ and UO₂²⁺ is investigated for the first time. Our simulation results demonstrate that the binding of Gd³⁺ and UO₂²⁺ onto the oxygens of crown ethers is favorable for polystyrene grafted crown ether in the organic solvents OCT and NB. The metal ion binding free energy (ΔG_Binding) in different solvent environments is calculated using the thermodynamic integration (TI) method. ΔG_Binding becomes more favorable in both solvents, NB and OCT, with an increase in the polystyrene monomer length. The metal ion transferability from an aqueous phase to an organic phase is estimated by calculating transfer free-energy calculations (ΔG_Transfer). ΔG_Transfer is significantly favorable for both Gd³⁺ and UO₂²⁺ for the transfer from the aqueous phase to the organic phase (i.e., NB and OCT) via ion-complexation to DBCE with an increase in polystyrene length. The partition coefficient (log P) values for Gd³⁺ and UO₂²⁺ show a 5-fold increase in separation capacity with polystyrene grafted DBCE. We corroborate the observed behavior by further analyzing the structural and dynamical properties of the ions in different phases.

Can metal ion complexes be used as polarizing agents for solution DNP? A theoretical discussion

Dynamic nuclear polarization (DNP) can be used to dramatically increase the NMR signal intensities in solutions and solids. DNP is usually performed using nitroxide radicals as polarizing agents, characterized by sharp EPR lines, fast rotation, fast diffusion, and favorable distribution of the unpaired electron. These features make the nitroxide radicals ideally suited for solution DNP. Here, we report some theoretical considerations on the different behavior of some inorganic compounds with respect to nitroxide radicals. The relaxation profiles of slow relaxing paramagnetic metal aqua ions [copper(II), manganese(II), gadolinium(III) and oxovanadium(IV)] and complexes have been re-analyzed according to the standard theory for dipolar and contact relaxation, in order to estimate the coupling factor responsible for the maximum DNP enhancement that can be achieved in solution and its dependence on field, temperature and relative importance of outer-sphere versus inner-sphere relaxation.

Hyperfine coupling constants on inner-sphere water molecules of a triazacyclononane-based Mn(II) complex and related systems relevant as MRI contrast agents

We report the synthesis of the ligand H2MeNO2A (1,4-bis(carboxymethyl)-7-methyl-1,4,7-triazacyclononane) and a detailed experimental and computational study of the hyperfine coupling constants (HFCCs) on the inner-sphere water molecules of [Mn(MeNO2A)] and related Mn(2+) complexes relevant as potential contrast agents in magnetic resonance imaging (MRI). Nuclear magnetic relaxation dispersion (NMRD) profiles, (17)O NMR chemical shifts, and transverse relaxation rates of aqueous solutions of [Mn(MeNO2A)] were recorded to determine the parameters governing the relaxivity in this complex and the (17)O and (1)H HFCCs. DFT calculations (TPSSh model) performed in aqueous solution (PCM model) on the [Mn(MeNO2A)(H2O)]·xH2O and [Mn(EDTA)(H2O)](2-)·xH2O (x = 0-4) systems were used to determine theoretically the (17)O and (1)H HFCCs responsible for the (17)O NMR chemical shifts and the scalar contributions to (17)O and (1)H NMR relaxation rates. The use of a mixed cluster/continuum approach with the explicit inclusion of a few second-sphere water molecules is critical for an accurate calculation of HFCCs of coordinated water molecules. The impact of complex dynamics on the calculated HFCCs was evaluated with the use of molecular dynamics simulations within the atom-centered density matrix propagation (ADMP) approach. The (17)O and (1)H HFCCs calculated for these complexes and related systems show an excellent agreement with the experimental data. Both the (1)H and (17)O HFCCs (A(iso) values) are dominated by the spin delocalization mechanism. The A(iso) values are significantly affected by the distance between the oxygen atom of the coordinated water molecule and the Mn(2+) ion, as well as by the orientation of the water molecule plane with respect to the Mn-O vector.

1H and 17O NMR relaxometric and computational study on macrocyclic Mn(II) complexes

Herein we report a detailed 1H and 17O relaxometric investigation of Mn(II) complexes with cyclen-based ligands such as 2-(1,4,7,10-tetraazacyclododecan-1-yl)acetic acid (DO1A), 2,2'-(1,4,7,10-tetraazacyclododecane-1,4-diyl)diacetic acid (1,4-DO2A), 2,2'-(1,4,7,10-tetraazacyclododecane-1,7-diyl)diacetic acid (1,7-DO2A), and 2,2',2"-(1,4,7,10-tetraazacyclododecane-1,4,7-triyl)triacetic acid (DO3A). The Mn(II) complex with the heptadentate ligand DO3A does not have inner sphere water molecules (q = 0), and therefore, the metal ion is most likely seven-coordinate. The hexadentate DO2A ligand has two isomeric forms: 1,7-DO2A and 1,4-DO2A. The Mn(II) complex with 1,7-DO2A is predominantly six-coordinate (q = 0). In aqueous solutions of [Mn(1,4-DO2A)], a species with one coordinated water molecule (q = 1) prevails largely, whereas a q = 0 form represents only about 10% of the overall population. The Mn(II) complex of the pentadentate ligand DO1A also contains a coordinated water molecule. DFT calculations (B3LYP model) are used to obtain information about the structure of this family of closely related complexes in solution, as well as to determine theoretically the 17O and 1H hyperfine coupling constants responsible for the scalar contribution to 17O and 1H NMR relaxation rates and 17O NMR chemical shifts. These calculations provide 17O A/ħ values of ca. 40 × 10(6) rad s(-1), in good agreement with experimental data. The [Mn(1,4-DO2A)(H2O)] complex is endowed with a relatively fast water exchange rate (k(ex)298 = 11.3 × 10(8) s(-1)) in comparison to the [Mn(EDTA)(H2O)]2- analogue (k(ex)298 = 4.7 × 10(8) s(-1)), but about 5 times lower than that of the [Mn(DO1A)(H2O)]+ complex (k(ex)298 = 60 × 10(8) s(-1)). The water exchange rate measured for the latter complex represents the highest water exchange rate ever measured for a Mn(II) complex.

Lanthanoids(III) and actinoids(III) in water: Diffusion coefficients and hydration enthalpies from polarizable molecular dynamics simulations

By using polarizable molecular dynamics (MD) simulations of lanthanoid(III) and actinoid(III) ions in water, we obtained ionic diffusion coefficients and hydration enthalpies for both series. These values are in good agreement with experiments. Simulations thus allow us to relate them to microscopic structure. In particular, across the series the diffusion coefficients decrease, reflecting the increase of ion–water interaction. Hydration enthalpies also show that interactions increase from light to heavy ions in agreement with experiment. The apparent contradictory result of the decrease of the diffusion coefficient with decreasing ionic radius is tentatively explained in terms of dielectric friction predominance on Stokes’ diffusive regime.

MRI stem cell tracking for therapy in experimental cerebral ischemia

Magnetic resonance has an established role in investigations on the evolution of stroke and the assessment of therapeutic strategies in experimental animals. Here we show that the technique has also an important place for the study of stem cell-mediated regenerative therapies after stroke. We review the literature by bridging from the methodological aspects of stem cell labeling via grafting and monitoring of cell dynamics after implantation into the brain all the way to MRI's role in analyzing the stem cell-mediated functional improvement. Thus, we have aimed at a view combining the focus on the monitoring of the cell activities with the aspect of lesion evolution while including also the essence of a potential functional improvement by the implantation of stem cells following stroke.

What first principles molecular dynamics can tell us about EXAFS spectroscopy of radioactive heavy metal cations in water

In this paper we show how molecular dynamics simulation can improve comprehension of structure and dynamics of water solvent around heavy cations. In particular, metal-water radial distribution function obtained from molecular dynamics can be used into EXAFS equation to improve the experimental signal fitting. Here we show results on structure and dynamics of Co2+, that is a radiocontaminant cation in its isotopic form 60Co, and lanthanoids(III) that are the chemical analogues of actinides(III) in aqueous solution.