|Home | About | Journals | Submit | Contact Us | Français|
Exchange of water molecules between the frequency-shifted inner-sphere of a paramagnetic lanthanide ion and aqueous solvent can shorten the T2 of bulk water protons. The magnitude of the line-broadening T2 exchange (T2exch) is determined by the lanthanide concentration, the chemical shift of the exchanging water molecule, and the rate of water exchange between the two pools. A large T2exch contribution to the water linewidth was initially observed in experiments involving Eu31-based paramagnetic chemical exchange saturation transfer agents in vivo at 9.4 T. Further in vitro and in vivo experiments using six different Eu31 complexes having water exchange rates ranging from zero (no exchange) to 5 × 106 s–1 (fast exchange) were performed. The results showed that the exchange relaxivity (r2exch) is small for complexes having either very fast or very slow exchange, but reaches a well-defined maximum for complexes with intermediate water exchange rates. These experimental results were verified by Bloch simulations for two site exchange. This new class of T2exch agent could prove useful in the design of responsive MRI contrast agents for molecular imaging of biological processes.
Diamagnetic chemical exchange saturation transfer (CEST) agents and paramagnetic CEST (PARACEST) agents use proton spin saturation and water exchange to produce exogenous-based contrast in MR images (1,2). Water protons bound to the inner sphere of a paramagnetic lanthanide ion (Ln ≠ La, Gd, or Lu) are typically shifted well away from the resonance frequency of bulk water protons (3). Radiofrequency saturation at the bound water frequency can produce indirect partial saturation of the bulk water signal through chemical exchange, resulting in negative contrast (i.e., darkening) in regions where the PARACEST agent is present. The percent reduction in bulk water signal is a function of agent concentration, radiofrequency saturation power and duration, and water exchange rate. Paramagnetic complexes can act as PARACEST agents if the water molecule exchange rate is slow when compared with the frequency difference between the two water pools (kex << Δω) (1,4). One potential advantage of PARACEST agents over Gd-based T1 agents is that image contrast can be turned on and off by inserting a frequency-selective radiofrequency saturation pulse before the imaging sequence. This frequency selection creates the potential of imaging multiple lanthanide-based agents simultaneously, each responding to a different physiological parameter or biochemical event. Such agents hold great potential to further extend the functional and molecular imaging capabilities of MR (5). Some published applications of PARACEST include measures of tissue pH (6,7), Zn2+ ion concentration (8), beta-cell function (9,10), and the tissue distribution of glucose (11–13).
During initial in vivo studies of Eu3+-based PARACEST agents in mice, we observed a significant loss (>50%) of MR signal in the kidneys (after intravenous injection) and tumors grown from human cancer cell xenografts (after intratumoral injection). This signal loss was present even before a frequency-selective radiofrequency saturation pulse was applied and appeared to be caused by a local decrease in water T2 because of the presence of the PARACEST agent. Similar phenomena in mouse kidneys were recently observed using a Tm3+-based PARACEST agent and the OPARACHEE (on-resonance paramagnetic agent chemical exchange effect) pulse sequence (14,15). We also observed that this effect was more pronounced when using polymerized versions of the Eu3+-based agents that had higher local concentrations of Eu3+ ions and thus more water exchange sites per molecule (16). As the Eu3+ ion in these complexes has a relatively small paramagnetic relaxation enhancement when compared with other lanthanide ions (e.g., Gd3+, Tb3+, Dy3+ or Tm3+) (17), it was clear that the T2 line-broadening effect must originate with some other chemical feature of these complexes. This led to the hypothesis that the same slow-to-intermediate water exchange characteristics of these complexes that enable the PARACEST effect can also cause a local decrease in the bulk water signal though T2 exchange (T2exch).
The effects of chemical exchange with diamagnetic molecules containing exchangeable —NH or —OH protons on transverse relaxation times have been known for over 40 years (18,19). Many early applications focused on using such molecules to suppress the strong water proton signal in high-resolution proton NMR experiments (20–24). Selective reduction of bulk water T2 by chemical exchange was first proposed as a method for MRI contrast by Aime et al. (25) and has recently been studied in vitro using the —NH-based proton exchange of iopamidol (26). In this article, we have extended this concept to weakly paramagnetic Eu3+ agents having variable water exchange rates. Both in vitro relaxation data and in vivo imaging data are presented along with simulations of T2exch using the Bloch equations modified for chemical exchange. It is shown that the transverse relaxivity (r2) of these agents is a function of water exchange kinetics and that r2 has a well-defined maximum at a specific exchange rate.
The water exchange properties of the six different Eu3+ complexes used in this study are summarized in Table 1. The exchange rates were previously determined by fitting CEST-spectra (Z-spectra) of these complexes, taken at several different applied amplitude of (excitation) radiofrequency fields (B1), to the Bloch equations modified for chemical exchange (17). Five different concentrations of each complex were prepared (e.g., 1.25, 2.5, 5, 10, and 20 mM). The pH and osmolarity for each sample were adjusted to approximately 7.0 and 300 mOsm, respectively. The T2 of each sample was measured at 0.54 T and 37°C using a Maran Ultra 23 MHz spectrometer, the Carr–Purcell–Meiboom–Gill sequence, and a nonlinear least-squares fitting procedure. The T2 of each sample was also measured at 9.4 T and 37°C using a Varian 400 MHz vertical bore spectrometer, the Carr–Purcell–Meiboom–Gill sequence, and a nonlinear least-squares fitting procedure.
Simulations of the T2exch contribution to T2 shortening of the water resonance using the Bloch-McConnell equations were written in MATLAB (MathWorks, Natick, MA). Details of the method were presented by Woessner et al. (27). The Mxy components of the magnetization, as predicted by the simulation, were fit to a single exponential decay to yield a T2 estimate as a function of chemical shift and exchange rate.
MRI data were acquired on a Varian (Santa Clara, CA) 9.4 T (400 MHz) small animal system using a 38-mm diameter quadrature volume coil. The mice were maintained at 37°C by a flow of warmed air and their temperature was monitored using a rectal thermocouple from Small Animal Instruments (Stony Brook, NY). Images were acquired using a fast spin-echo sequence (pulse repetition time/echo time = 2500/9.12 ms, echo train = 8, average = 4). Each coronal image of the kidneys (128 × 128 × 1 pixel matrix, 38 × 38 × 2 mm field-of-view) had an in-plane resolution of 297 μm and took 3 min to acquire. Three preinjection baseline images were acquired before administration of a 0.25 mmol kg–1 dose of agent to the healthy female black six mice via tail vein injection. Images were then continuously acquired (one every 3 min) for 1 hour. A region-of-interest drawn over each kidney was used to measure the mean pixel intensity in each image. Previous PET/CT studies of untargeted Eu3+ PARACEST agents (labeled with 64Cu) show that on tail vein injection of similar doses, the kidneys immediately begin to filter the agent from the blood and into the bladder with signal remaining in the kidneys for more than 1 hour (28).
Figure 1a shows a plot of transverse relaxation rate (R2) versus Eu3+ concentration (mM) for the six compounds listed in Table 1 collected at 0.54 T (23 MHz) and 37°C. Linear fits were used to calculate the transverse relaxivity (r2) for each compound. It can be seen that r2 is small (<0.005 s–1 mM–1; Table 1) for compounds with very fast or no water exchange (EuDOTA– and EuTETA–, respectively) but larger (0.020–0.050 s–1 mM–1) for compounds with intermediate water exchange. Figure 1b shows a similar plot for all six compounds collected at 9.4 T (400 MHz) and 37°C. It can be seen that r2 is again small (<0.055 s–1 mM–1) for compounds with very fast or no exchange and larger (0.216–0.527 s–1 mM–1) for compounds with intermediate exchange.
There are two notable differences between the data collected at these two B0 fields. First, the r2 values at 9.4 T are larger than those at 0.54 T, as expected. The second less obvious, yet highly informative difference is that the order of r2 for the four intermediate exchange rate compounds is different at the two fields. The latter observation provides new insights into the mechanism and potential uses of such agents for molecular imaging. The reasons for these r2 differences are thus explained.
The total water proton T2 (T2tot) of an aqueous solution containing a paramagnetic complex undergoing water chemical exchange is given by:
where T2water is the transverse relaxation time of pure water and the latter two terms are the paramagnetic (T2para) and exchange (T2exch) contributions due to the paramagnetic ions. Given that the Eu3+ ion is only weakly paramagnetic, the 1/T2para term can be neglected except for those situations where the concentration of Eu3+ is very high. Equation 1 can then be simplified and rewritten as:
where R2B is the transverse relaxation rate of the bound water, τB is the bound water lifetime, Δω is the chemical shift difference (in rad s–1) between the two water exchange sites, and PB is the mole fraction of the paramagnetic ions (e.g., 1 mM/55.5 M). For the limiting case of Δω >> R2B, Eq. 3 reduces to:
so that T2exch is now only a function of τB, Δω, and PB (19).
Equation 4 can be used to illustrate the dependence of R2 on the bound water lifetime (τB) for a given molar fraction of paramagnetic ions. By using a 1-mM concentration of paramagnetic ions Eq. 4 can be expressed in terms of exchange transverse relaxivity (r2exch) as a function of τB:
A plot of r2exch versus τB (hereafter called a “Swift–Connick” plot based on the assumptions leading to Eq. 5) for Eu3+ at 0.54 T (Δω = 42.5 ppm or 6.14 × 103 rad s–1) is shown in Fig. 2a. The peak in the Swift–Connick plot can be used to delineate fast, intermediate, and slow water exchange regions. It can be seen that for fast and slow exchange rates (τB < 10–6 s and τB > 10–2 s, respectively) r2exch approaches zero, while for intermediate exchange rates (10–6 s < τB < 10–2 s) r2exch reaches a well-defined maximum (0.055 s–1 mM–1) at τB Δω–1 (163 μs). The previously measured bound water lifetimes (Table 1) and transverse relaxivities (Fig. 1a) were used to plot markers for the six Eu3+ complexes in Fig. 2a. The measured data are in good agreement with the Swift–Connick plot showing that r2exch is highly dependent on the water exchange rates. Figure 2a predicts that a Eu3+ complex with a bound water lifetime of 163 μs at 37°C would have the largest r2exch at 0.54 T.
A plot of r2exch versus τB for Eu3 at 9.4 T (Δω = 42.5 ppm or 1.07 × 105 rad s–1) is shown in Fig. 2b (with the data from 0.54 T superimposed). Figure 2b illustrates that when moving from 0.54 to 9.4 T the maximum r2exch value increases proportionally (from 0.055 to 0.962 s–1 mM–1) with the increase in B0 and that the bound water lifetime for the maximum r2exch moves toward faster water exchange (from 163 to 9.4 μs). The reduction in peak τB at 9.4 T causes the r2 relations between the four intermediate compounds to be inverted when compared with the 0.54 T data. The water exchange rates for the four intermediate compounds place them on the rising (or faster) side of the Swift–Connick peak at 0.54 T (Fig. 2a), while the same water exchange rates place them on the falling (or slower) side of the Swift–Connick peak at 9.4 T (Fig. 2b). This is interesting because it means that while has the smallest r2 value of the four intermediate exchange complexes at 0.54 T, it has the largest r2 value at 9.4 T. Clearly, the experimental data are in good agreement with Swift–Connick theory.
To confirm that the observed variations in transverse relaxivity were due to differences in water exchange rates, r2 was also measured at different temperatures for two of the six Eu3+ compounds studied. Figure 3a shows a plot of transverse relaxation rate (R2) versus Eu3+ concentration (mM) for at three different temperatures collected at 9.4 T. It can be seen that r2 is proportional to temperature, consistent with being on the falling or “slow” side of the Swift–Connick peak at 9.4 T in Fig. 2b. Figure 3b shows a similar plot for EuDOTA–, where r2 is now inversely proportional to temperature, consistent with EuDOTA– being on the rising or “fast” side of the Swift–Connick peak at 9.4 T. The bound water lifetimes (τB) predicted by fitting the measured r2 values to the Swift–Connick peak in Fig. 2b are in agreement with recent exchange rate data for similar EuDOTA-based compounds (30).
The Bloch-McConnell simulation results were in exact agreement with predictions made by the Swift–Connick equations describing T2 in exchanging systems, emphasizing the robustness of the Swift–Connick model.
Figure 4 shows coronal slices through the kidneys of three different mice before and 30 min after the injection of (a) EuTETA–, (b) , and (c) EuDOTA– (one compound per mouse). During renal clearance of either EuTETA– or EuDOTA– (slow and fast exchange, respectively), there was only a slight darkening of the renal pelvis. However, dramatic darkening of the renal pelvis, medulla, and cortex occurred after injection of , the intermediate water exchange compound. These contrast effects are readily explained by the differences in water exchange rates and thus r2exch of each compound. These results also explain our prior observations where severe darkening of a tumor occurred after direct injection of a polymeric form of in vivo (16,28). These data provide a benchmark for water exchange conditions that will be required in creating PARACEST agents that can be activated in vivo without producing a significant T2exch effect.
The excellent agreement between the experimental in vitro data at both 0.54 T and 9.4 T and the Swift–Connick predictions supports our hypothesis that the negative contrast observed in vivo is largely due to T2 exchange (T2exch) for Eu3+ complexes. As the T2exch effect is concentration dependent, the kidneys are a likely tissue in which to observe this phenomenon as the Eu3+ concentration within the renal nephrons is much higher than that of normal vascular space. Also, as the effect depends on B0, the observed contrast at 9.4 T for a given Eu3+ concentration will be significantly reduced at 3.0 T or 1.5 T. This suggests that PARACEST agents with intermediate water exchange rates might be best observed at lower magnetic fields where the T2exch effects are less prominent.
The Swift–Connick equation, giving a simple relationship between r2exch and τB as illustrated in Eq. 5, can be used to evaluate water exchange rates for new complexes where this parameter is unknown. The r2 of a new compound can be obtained by measuring R2 as a function of agent concentration (Fig. 1) and the water exchange rate (or τB) estimated by finding the intersection of this r2 on the Swift–Connick plot. The 2-fold degeneracy (i.e., fast or slow side of the Swift–Connick peak) can be resolved by comparison with other similar chemical structures with known bound water lifetimes, measuring the r2 at two distinct B0 fields (e.g., 0.54 T and 9.4 T), or by measuring the temperature dependence of r2 over a small range of temperatures (e.g., ±10°C) (31).
The dependence of r2exch on Δω can be predicted by plotting the Swift–Connick equation for other lanthanide ions as each ion will have a distinct Δω for the exchanging water molecule (3). Figure 5 shows that as Δω increases, the maximum r2exch increases and moves toward faster exchange, with Dy3+ giving the largest maximum r2exch (16.3 s–1 mM–1) for complexes with fast water exchange (τB = 0.8 μs). Chelates of Dy3+ (e.g., Dy-DTPA-BMA or sprodiamide) (32) have been proposed as magnetic susceptibility agents in both ex vivo and in vivo studies (32–37), where the concept of T2exch was suggested to explain some of the observed differences when compared with chelates of Gd3+ (38,39). Although the theoretical maximum r2exch of a Dy3+-based compound is an order of magnitude less than that of super paramagnetic ion oxide nanoparticles (40), the contrast mechanisms between the two are quite different. Also, by using a paramagnetic rather than a diamagnetic T2exch agent, the same level of contrast can be achieved at much lower concentrations (26).
It should be noted that Fig. 5 illustrates the effect due to r2exch alone. Strongly paramagnetic ions like Ho3+, Tb3+, and Dy3+ (where T2para can not be set to zero as in Eq. 2) will also have a paramagnetic transverse relaxivity term (r2para) added to the total r2. Including the r2para term, which is independent of the exchange rate, would simply shift the Swift–Connick plots for these ions in the positive vertical direction by an amount proportional to their paramagnetic relaxation enhancement (17). Magnetic field inhomogeneity also contributes an additional line-broadening term (effective transverse relaxation time, ) to T2 in vivo but this contribution is small when compared with the r2exch contribution for these intermediate-to-slow water exchange Eu3+ complexes as illustrated by the in vivo data of Fig. 4. This may not be the case for other lanthanide-based agents such as Dy3+, where may play a more significant role.
We have presented the first in vitro and in vivo data from weakly paramagnetic Eu3+ compounds acting as T2exch contrast agents. We have shown that T2exch contrast is proportional to the water exchange rate and that there is an optimal exchange rate for each paramagnetic lanthanide ion and magnetic field (both determined by Δω). Although the T2exch contrast mechanism is relatively weak when compared with paramagnetic gadolinium T1 agents or iron oxide T2* agents, this discovery could nonetheless be useful in the design of new types of molecular imaging agents for MRI. One example would be to create a low molecular weight polymer of a DyDOTA-based agent with an intermediate water exchange rate (r2exch approximately 10–16 s–1 mM–1 at 9.4 T; Fig. 5) to amplify the T2exch contribution (16). A modest 20-fold degree of polymerization would result in a total molecular r2exch of 200–320 s–1 mM–1 that would add to the inherent paramagnetic T2 relaxation effects of Dy3+. By using such an approach, one might be able to create a highly sensitive T2 contrast agent with an r2 as comparable with (or even greater than) super paramagnetic ion oxide nanoparticles while retaining the benefits of using a small molecule. We believe this platform could prove useful for detecting cancer-related enzyme activities in vivo.
The authors thank Dr. James Ratanakar and Dr. Subha Viswanathan for their assistance with the T2 measurements and Tomoyasu Mani for synthesizing one of the compounds used in this work.
Grant sponsor: National Institutes of Health; Grant number: RR-02584 and CA-115531; Grant sponsor: Robert A. Welch Foundation; Grant number: AT-584.