|Home | About | Journals | Submit | Contact Us | Français|
Magnetic Resonance Imaging (MRI) has been immensely valuable in diagnostic clinical imaging over the last few decades owing to its exceptional spatial and anatomical resolution. The signal in MRI is generated by relaxation of the transverse component of the net magnetization of protons present in the body, predominantly from bulk water. Thus, any agent or process that affects the net magnetization of the water protons in body tissues will also influence image contrast. Gd3+-based contrast agents shorten both the longitudinal and transverse relaxation times (T1 and T2) of water protons to approximately the same extent, in essence by relaxing all nearby proton spins. This effect is detected as increased signal intensity in T1 weighted MRI images when the appropriate pulse sequence is applied. Over the past 25 years Gd3+-complexes have been spectacularly successful as extracellular or blood pool T1 agents but their relative insensitivity to changes in environment coupled with the fact that they are never completely silent limits their applicability in the design of responsive MRI agents. A conceptually different approach to contrast enhancement is based on chemical exchange saturation transfer (CEST). This technique relies on dynamic chemical exchange processes inherent in biological tissues to transfer saturated 1H spins into the bulk water proton pool, which leads to a decrease of net magnetization and is detected as a negative contrast (darkening of the image) in MRI. Originally exchangeable -NH and -OH protons of various biomolecules were used to generate CEST contrast (DIACEST). However, these agents suffer from a few drawbacks, particularly in association with the small, usually less than 6 ppm, chemical shift difference between the two exchanging pools. The great benefit of using paramagnetic hyperfine shifting lanthanide complexes as CEST agents (PARACEST) is that the chemical shift difference between the two exchanging pools can potentially be much larger, up to several hundred ppm, facilitating easy saturation of one of the exchangeable spin pools without partial saturation of the bulk water pool. Another advantage of PARACEST is that the exchangeable sites are not limited to -NH or -OH protons but sites with faster exchange rates such as a Ln3+-bound H2O molecule, in particular, can also be considered. Since the water exchange rate on lanthanide complexes is extremely sensitive to the chemical environment, this has created unprecedented opportunities in the design of responsive PARACEST agents. In addition, multi-frequency MRI imaging is inherent to PARACEST: multiple agents present in the body can be imaged in one experiment by selectively turning on and off each agent by applying the appropriate saturation frequency.
This article on alternatives to Gd3+-based MRI contrast agents focuses mainly on paramagnetic CEST (PARACEST) agents with special emphasis on their tremendous advantages as responsive agents. Since the conception of CEST for contrast enhancement in MRI,1 numerous reviews have dealt with different aspects of CEST-based contrast agents, both diamagnetic and paramagnetic.2–12 The goal of this review is to offer a comprehensive discussion of the theoretical background of chemical exchange as applied to CEST and in turn examine the fundamental aspects of the different parameters of paramagnetic lanthanide complexes that aid in the design of novel PARACEST agents. The first section on paramagnetic properties of the lanthanide ions provides a basic understanding of the advantages of Ln3+ complexes in general as MRI contrast agents. The historical and theoretical background of chemical exchange, the initial use of CEST for contrast enhancement in diamagnetic molecules and the various factors influencing water exchange processes in lanthanide ions facilitate a better understanding of the rational design of novel PARACEST agents that are discussed in extensive detail in the following sections. An entire section is devoted to the evaluation of the responsive behavior of these agents reported up to the publication date of this review, targeting ten different classes of responsive targets. Similar to Gd3+-based contrast agents, small molecule PARACEST agents also face problems associated with detection sensitivity and retention time in vivo. Several examples of the efforts of various groups to tackle this problem and the success attained so far are discussed. We also delve into the important in vivo characteristics of these compounds by discussing the kinetic inertness, thermodynamic stability and biodistribution of selected lanthanide complexes before exploring the problems associated with the in vivo imaging of these agents. Finally, a few other paramagnetic lanthanide complexes based on Tm3+, Dy3+ and Eu2+, that could potentially be used as MRI contrast agents and not based on CEST are discussed briefly in a separate section.
All Ln3+ ions, with the exception of La3+ and Lu3+, have unpaired electrons and are thus paramagnetic. These unpaired electrons are located in the 4f orbitals, which are lower in energy than the unfilled 5d orbitals and do not overlap significantly with ligand orbitals. For this reason, the 4f electrons delocalize only to a small extent onto ligand orbitals, rendering the bonding in Ln3+ complexes predominantly ionic. Consequently, the coordination geometry is determined largely by steric rather than electronic factors. As the f-orbitals are largely shielded from the ligand, the crystal field effects are insignificant and the electronic and magnetic properties of Ln3+-ions are not influenced appreciably by the coordinating ligand.13
Unlike their chemical behavior, the magnetic properties of the lanthanides differ dramatically along the series. Each of the lanthanide elements has a different magnetic moment and magnetic susceptibility that is determined by the electron configuration (Figure 1). For example, trivalent lanthanide ions that have half-odd integer values viz Nd3+, Sm3+, Gd3+, Dy3+, Er3+ and Yb3+ (Kramers’ ions) are associated with a permanent magnetic moment in the absence of a magnetic field, while the rest of the lanthanide ions are not (non-Kramers’ ions).14–15 Magnetic susceptibility is a 2nd rank tensor that governs the magnitude and directionality of the interaction between the magnetic moment and the external magnetic field. The Ln3+-ion affects nuclei in its vicinity mainly through paramagnetic relaxation enhancements (PRE) and lanthanide induced shifts (LIS), both of which depend on the magnetic susceptibility of the lanthanide ion.16–19 PRE is generated as a result of the interaction between fluctuating magnetic fields created by the unpaired electrons with the nuclear spin and is manifested as the shortening of the nuclear relaxation time 1/Ti,M, where i = 1, 2. The relaxation enhancing effect can be expressed as the sum of three components: contact (scalar), dipolar and Curie relaxation enhancement (Eq 1).17
The contact relaxation enhancement is a through-bond effect and can be quite significant for directly coordinated atoms but its strength rapidly decreases as the number of bonds increases. For nuclei that are not directly coordinated to the paramagnetic center, the other two mechanisms (dipolar and Curie relaxation enhancement) depend on the electronic relaxation time T1e of the paramagnetic metal ion. Qualitatively, the dipolar relaxation enhancement (Solomon effect) is most efficient when the electronic relaxation rate of the metal matches the Larmor frequency of the nucleus to be relaxed. Therefore, the dipole-dipole relaxation enhancement is dominant only for metal ions with relatively long electronic relaxation times. Gd3+ is unique among the lanthanides because its symmetric seven electron ground state (8S7/2) results in an electronic relaxation rate that is six orders of magnitude slower than the other Ln3+ ions.17–20 The unpaired electrons in Ln3+ ions other than Gd3+ by necessity populate the f orbitals anisotropically which gives rise to low energy excited states. This in turn results in a strong magnetic anisotropy and fast electronic relaxation with very short T1e, on the order of 10−13 s.21–22 The relaxation enhancement of these metal ions can have significant contribution from the Curie-spin relaxation mechanism, which arises from the interaction of the nuclear spin with the thermal average of the electron spin.17,19,23 While the dipolar relaxation enhancement is proportional to the isotropic magnitude of the magnetic susceptibility tensor (χ) of the Ln3+-ion, the Curie-spin relaxation is proportional to the square of χ. The Curie spin relaxation effect is significant at lower temperatures and higher magnetic fields and for ions whose magnetic moment is large. Generally, Curie spin relaxation affects the transverse relaxation more than the longitudinal relaxation. As a result, paramagnetic Ln3+-ions other than Gd3+ are less efficient T1 relaxation agents; but the strongly paramagnetic ions (Tb3+, Dy3+, Ho3+, Er3+) can have a significant effect on NMR linewidths. Such line broadening often makes it difficult or even impossible to detect nuclei that are located within a certain distance of the lanthanide ion. For example, the radius of this “blind sphere” in which heteronuclear 1H-13C cross peaks cannot be detected due to excessive line broadening of the 1H signal (>100 kHz) ranges from less than 3 Å for the weakly paramagnetic Eu3+ to about 12–14 Å for the strongly paramagnetic Dy3+. However, because paramagnetic relaxation effects scale with r−6 while pseudocontact shifts scale with r−3, lanthanide induced shifts are easily detectable at a distance where the relaxation effect is already negligible.17,19,24 The dipolar relaxation enhancement (Solomon effect) and Curie spin relaxation enhancement has been reviewed in detail elsewhere.17,19–20,25–26
The diamagnetic shift is caused by inductive and electric field effects and its value can be estimated using a diamagnetic lanthanide ion (La3+, Lu3+) or the pseudo-lanthanide Y3+. The origin of contact shift is the through bond interactions of the unpaired 4f electron spins with the observed nucleus. Contact shift can be quite large for nuclei directly coordinated to the lanthanide ion but its magnitude decreases rapidly as the number of bonds between the lanthanide and the observed nucleus increases. The strength of the contact interaction is dependent on the hyperfine coupling constant A/ђ and the spin expectation value <SZ> which is characteristic for each lanthanide (Eq 3).28–29
The magnitude of the Fj factor is determined by the nature of the nucleus i and independent of the lanthanide j and the structure of the ligand. The isotropic hyperfine coupling constant A/ħ is included in Fi (Eq 4).19
where μB is the Bohr magneton, k the Boltzmann constant, γi the gyromagnetic ratio of nucleus i and T the absolute temperature.
where r, θ and ϕ are the polar coordinates of nucleus i in the reference system formed by the principal axes of the magnetic susceptibility tensor. D1j and D2j are magnetic anisotropy constants proportional to Bleaney’s constant, Cj, and and are the second order crystal field coefficients, which reflect the strength of interaction between the metal and donor atoms. For a series of isostructural complexes the G factor and the crystal field coefficients are independent of the lanthanide and the pseudocontact shift produced by a Ln3+-ion is proportional to its Cj value. The lanthanide independent parameters are often combined to give Eq 8 (McConnell-Robertson). 24,27
The lanthanide dependent <Sz> and Cj values have been determined for each lanthanide and are listed in Table 1.
According to Bleaney’s theory,31 the second order crystal field parameters are proportional to the anisotropic (axial and rhombic) components (χax and χrh) of the magnetic susceptibility tensor and so the magnitude of the pseudocontact shift can be expressed as in Eq 9.
Since the magnitude of the tensor anisotropy Δχ is quite different for each lanthanide ion, the magnitude of the pseudocontact shift at each nucleus in a complex will vary significantly depending upon its distance from the lanthanide ion, ranging from only those nuclei in the immediate vicinity of the complex (Sm3+) to over 40 Å (Dy3+). When θ = 54.74° and ϕ = 45, the dipolar field vanishes and pseudocontact shift is zero (Figure 2). The Ln3+-induced dipolar field thus has a spatial directionality and consequently the pseudocontact shifts experienced by each nuclei of the ligand provides three-dimensional structural information about the complex. Hence, Eq 9 forms the basis of the application of lanthanide shift reagents in the determination of three dimensional structures of proteins and other molecules.17–18,21,32
For axially symmetric complexes or when conformational averaging of the Ln3+-shift reagent-substrate adduct occurs, the second term is zero. In this case the dipolar field is determined by 3(cos2 θ − 1) and the paramagnetic shift is given by Eq 10.
For a series of Ln3+ complexes, the lanthanide induced shift data should give a linear relationship when plotted according to Eq. 11 or 12 provided that the complexes are isostructural. In this case Fi and are invariant along the lanthanide series and their values can be determined by fitting either equation, thereby allowing separation of the contact and pseudocontact shifts.33–35 However, the assumption that the complexes are isostructural along the lanthanide series does not always hold true and gradual changes in the coordination geometry can give rise to variations of the crystal field coefficient and the hyperfine coupling constant across the series. Hence, plots of versus or versus do not always give a linear correlation but a break is often observed around the middle of the series near gadolinium. Crystal field independent methods that treat this deviation from linearity have been developed. 36–39 Another general method of LIS analysis that requires no assumptions about the orientation of principle magnetic axis system uses a linear least squares fitting of the measured LIS data to equations involving components of the magnetic susceptibility tensor. The best fit of the LIS data is determined through permutations of the LIS values over a number of nuclei. A computational method (SHIFT ANALYSIS program) based on this approach has also been developed.40–41 The SHIFT ANALYSIS program combined with molecular mechanics (MM+) calculations is a powerful method for the solution structure determination of various Ln3+-complexes and its effectiveness has been demonstrated with several lanthanide complexes including Ln3+-DOTA-tetraamide type chelates.41–42
Finally, it is worth noting that a simple method to elucidate the number of inner sphere water molecules in lanthanide complexes under fast exchange conditions is based on the measurement of the bulk water 17O NMR shifts.43–44
Since clinical MRI mainly involves the imaging of bulk water protons because of their high abundance in tissues, our focus here will be on contrast mechanisms that involve bulk water in chemical exchange with the labile protons (amines, amides, guanidinium protons, alcohols), either endogenous or exogenous. Chemical exchange in such cases can involve the conventional magnetization transfer (MT) techniques, which involves exchange of magnetization between a semi-solid macromolecular phase and bulk water or the more recent CEST techniques that involve exchange between protons of solutes and bulk water.45–46 Conventional magnetization transfer experiments have found numerous applications since first demonstrated for in vivo imaging by Balaban.47 One of the major applications is magnetic resonance angiography (MRA) and characterization of white matter disease in the brain (mainly multiple sclerosis) among others.45 Although the exact mechanism underlying magnetization transfer has been under debate, it is generally assumed that both chemical exchange and dipolar interactions play an important role.12,48 There have been numerous reviews focused on magnetization transfer techniques and their applications;45,49–50 here we discuss saturation transfer techniques as applied to exchange between protons of solutes (both endogenous and exogenous) and bulk water, while briefly sketching the effects of inherent MT effects in vivo on CEST. These techniques have garnered significant attention since the first report by Balaban and Ward demonstrated that saturation of the labile protons of a low concentration solute (contrast agent) can result in a collective buildup of spin saturation in bulk water protons as a result of chemical exchange.1,51–52 Before discussing CEST spectroscopy and imaging, a brief background on the theory of chemical exchange as well as a description of the modified Bloch equations for two-site chemical exchange is given.
The use of NMR to study chemical exchange has been around since the early 1950s, discovered around the same time as scalar spin-spin coupling.53–56 Gutowsky and co-workers were among the first to cite examples of intramolecular exchange as indicated by the 1H or 19F multiplets in compounds such as PH3, PF3, F2PO(OH) and BrF5.54 Since then different types of chemical exchange processes have been studied by NMR including intramolecular processes involving motions of side chains of proteins, conformational equillibria and unfolding of proteins or intermolecular processes involving enzyme catalyzed reactions, binding of small molecules to macromolecules, protonation/deprotonation reactions and exchange of labile protons on a macromolecule with the solvent. NMR chemical exchange refers to any process whereby a nucleus exchanges between two or more environments with different NMR parameters (chemical shift, relaxation or scalar coupling).57 Figure 3 shows a classic example of chemical exchange between two equally populated spin states. At low temperatures, the two methyl groups are magnetically non-equivalent and exhibit two separate signals suggesting slow exchange. As the temperature is raised, the two peaks begin to coalesce and eventually sharpen into a single peak at higher temperatures. The coalescence at higher temperatures is a result of the reaction rate matching the difference in Larmor frequencies between the two sites (magnetic environment).
Many NMR techniques have been used to quantify exchange processes, including lineshape analysis, two-dimensional NMR spectroscopy and double resonance techniques among others.54,58–59 These techniques have been reviewed in numerous other reports and textbooks.53,60–62 An alternate practice to study exchange was first performed by Forsen and Hoffman in 1963, who used saturation transfer or magnetization transfer to measure proton transfer rates between salicylaldehyde and water.7,63 Saturation transfer experiments usually involve perturbation of one of the pools by an external radiofrequency field; in other words, the longitudinal relaxation pertaining to this signal is selectively perturbed. Upon exchange, saturation transfer can be monitored by observing the change in intensity of the other signal. This technique has since been exploited in numerous investigations.47,51,64–66
CEST can occur in any slow-to-intermediate chemical exchanging system involving nuclei in two or more magnetically nonequivalent environments or “pools”67 and typically requires separate resonances for each of the pools under consideration for magnetization transfer (Figure 3, bottom spectrum). The condition for slow-to-intermediate exchange on the NMR timescale is classically defined by Eq 13: the rate of exchange (kex) must be no greater than the frequency difference between the two pools, Δω (ω is expressed in units of radians.sec−1).
Consider two pools of protons, A and B with pool A representing the bulk water protons and pool B representing the pool of labile protons associated with a contrast agent. In this special case, Δω in Eq 13 can be considered to represent the frequency difference between the exchange site on the contrast agent and bulk water. Eq 13 can be re-written as
Here, τm represents the lifetime of the labile proton on the contrast agent and is just the reciprocal of the exchange rate constant, kex. As briefly described above, the basic CEST experiment first involves the application of a selective RF pulse at the Larmor frequency of pool B for a sufficiently long period of time so as to reach a new level of steady state magnetization. At that point, the magnetization along all three rotating frames are time-independent. This results in the distribution of spins in pool B moving closer to saturation; in other words, the number of spins aligned with and against the field in the Boltzmann distribution of nuclear spins becomes almost equal, resulting in a decrease in the bulk magnetization of pool B. If the exchange condition (Eqs 13 and 14) is satisfied, there will be a subsequent decrease in the bulk magnetization of pool A, thereby translating into a decrease in the signal intensity of the bulk water protons.4 In most cases involving imaging, the bulk water protons (pool A) are always present at much higher concentration in tissues (~80 M) than pool B protons, typically no more than millimolar. Thus, there is constant accumulation of non-saturated protons in pool B from A accompanied by replacement of a large number saturated protons in pool A from B and, since one always monitors bulk water, even a 10% saturation of the pool B protons can result in a significant decrease in pool A intensity.
The frequency difference between the two pools of protons (Δω) plays an important role in determining the efficiency of saturation transfer. The larger the frequency difference between the bulk water protons (pool A) and the protons of the contrast agent (pool B), the faster the exchange rate can be without approaching the fast exchange condition (Eqs 13 and 14). A large frequency difference also implies greater specificity of the selective radiofrequency pulse used to saturate pool B, without directly saturating the bulk water protons. Also, as will be discussed later, large Δω values results in a significant reduction of the macromolecular-water background (tissue MT effect).1 In other words, if the Δω values are sufficiently large, then magnetization transfer due to saturation of macromolecules and their subsequent exchange with bulk water has a smaller contribution to the CEST spectrum.
One can use the well-known NMR Bloch equations modified for chemical exchange to predict the NMR behavior of the two exchanging pools of proton spins. The Bloch equations, traditionally used to derive nuclear magnetizations along the x, y, and z axes,68–69 were first modified by McConnell and Thompson in 1957 to include the terms for chemical exchange.70 The simplest case of a two-site exchange was described by Woessner to obtain the variation of x, y and z magnetizations with time for each of the two pools (Eqs 15–20).12,67
T1a and T2a are the spin-lattice and transverse relaxation times of the bulk water pool A in the absence of exchange, ω is the Larmor frequency at which the selective radiofrequency irradiation pulse is applied, ωa is the Larmor frequency of the protons of poolA and ωb is the Larmor frequency of pool B. ω1 is the nutation rate of the selective RF irradiation , Ca is the rate at which spins leave pool A, and Cb is rate at which spins leave pool B. By definition then, Cb is the inverse of τb, the lifetime of protons in pool B (also often referred to as τm). The Z magnetization at thermal equilibrium in the two pools, and , are proportional to the number of protons in each pool as determined by the chemical composition of the exchanging system.4 Since chemical exchange is restricted by mass balance, Eq 23 holds true.
In a CEST experiment, when a selective RF pulse of frequency ω is applied at the Lar mor frequency of the pool B protons until steady state is achieved, the time derivatives of all nuclear magnetizations are equal to zero (left hand side of Eqs 15–20). The resulting set of equations, solved algebraically, give the net magnetization of pool A protons (bulk water) along the z axis, for a selected range of irradiation frequencies. It is clear from these equations that the saturation transfer efficiency of a contrast agent depends on numerous factors: (a) the time constant for proton exchange τb, (b) the saturation power B1, (c) the Larmor frequency difference between the two exchanging sites Δω, (d) the number of exchanging sites per agent ( , where c is the concentration of the contrast agent and q is the number of exchanging sites in the agent), and (e) the relaxation rates of both proton pools.
Owing to the complexity of these equations, and the factors that affect chemical exchange, numerous approximations have been made.71–75 The earliest approximation and most widely used approach made by Forsen and Hoffman in their landmark paper on chemical exchange involves the assumption that complete saturation of pool B protons can be achieved by using a frequency selective, weak saturation pulse such that there is negligible direct saturation of pool A protons.63,76 With these assumptions, the last two terms of Eq 19 can be omitted: can be dropped since we assume complete saturation of the pool B protons ( ) and can be neglected when there is negligible direct saturation of the pool A (ω1 ~ 0, Eq 25). Under steady state conditions, the time derivative equals zero and a simplified solution is obtained (Eqs 25–28).2
Here, c is the concentration of the contrast agent and q is the number of exchangeable proton sites, with 55.5 representing the molar concentration of bulk water. In other words, the protons that leave pool A with the z-magnetization will be replaced by protons from pool B with no magnetization; consequently decreasing the bulk magnetization of pool A at the rate, . At the same time, the spin-lattice relaxation will simultaneously return the z-magnetization of pool A to its equilibrium value at a rate , so at steady state, these rates become equal (Eq 29).
While Eq 29 provides a simple solution, it only makes predictions for limiting cases. In particular, it should be noted that the pool B protons can reach complete saturation only using a strong RF field, and under ideal situations and in most cases this system never actually reaches complete saturation.12
In CEST experiments, the effect on bulk water is usually displayed as a function of saturation frequency (ω) relative to the magnitude of the remaining water proton signal (Ms/M0), where Ms is the water signal intensity with saturation and M0 is the signal intensity without saturation. Grad and Bryant first coined the term Z-spectrum and used it to describe the change of intensity of one resonance as a function of the frequency of a saturating pulse applied selectively to another resonance in chemical exchange with the observed resonance.77 To obtain a Z-spectrum, one simply applies a selective radiofrequency pulse for a favorable duration followed immediately by a traditional spin-echo observation pulse to attain the net signal intensity of the bulk water. Repeating this experiment over a range of frequencies gives a complete Z-spectrum (equivalent to a CEST spectrum) for the contrast agent under study. An example of a simulated Z-spectrum is displayed in Figure 4 (simulations done using MATLAB 7, Mathworks Inc., Natick, MA).67 The CEST effect is also more commonly represented as a percent decrease in total bulk water intensity (Eq 30).
As first demonstrated by Wolff and Balaban, and explained later in detail in Section 22.214.171.124, imaging proton chemical exchange processes requires collection of two images; the first collected after selective saturation at the exchanging site of interest (Ms), while the second image is collected at the opposite side of the water resonance from the site of interest (M0).51,64 Pixel-by-pixel division (Ms/M0) or subtraction of the images (M0−Ms) yields a CEST image which reflects the contrast enhancement as a result of the agent.1 Taking the simulated spectrum in Figure 4 as an example, a selective radiofrequency pulse is first applied at +25 ppm corresponding to the decrease in magnetization of pool A bulk water protons due to exchange with pool B. This results in decreased bulk water signal intensity (MS) and consequently a darker image (negative contrast). Following this, another radiofrequency pulse is applied at −25 ppm corresponding to the opposite end of the Z-spectrum where no CEST effect is seen. Since there are no exchanging sites at this frequency, there is no decrease in magnetization of bulk water and hence a brighter image is obtained (M0). Subtraction of these two images, (M0 − Ms), yields a CEST image which reflects only the effects of the contrast agent.
The use of CEST techniques in the design of contrast agents for MRI was first demonstrated by Balaban and Ward in early 2000 (Figure 5).1 They evaluated different types of exchangeable proton sites in small molecules such as sugars, amino acids, nucleosides and other diamagnetic molecules. The authors found that the nitrogen bound protons of amines, amides and guanidinium groups were most effective in reducing the signal intensity of bulk water by saturation transfer. Table 2 lists the findings of some representative molecules reported in their work. It was shown that the CEST properties of these small molecules are dependent on numerous factors including pH, temperature and, of course, concentration of contrast agent. Some general conclusions can be drawn from this small library of CEST compounds:12 (a) these agents are inherently pH and temperature sensitive, (b) a single compound with multiple types of exchange sites can display different CEST properties depending on the frequency of saturation, (c) as the exchange rate of the compound increases, the CEST peak becomes broad and shifts towards bulk water, (d) the CEST efficiency of the agent increases linearly with its concentration, but the line shape of the CEST peak is maintained, and (e) one of the most attractive features of these contrast agents over traditional T1 or T2-based contrast agents is the ability to turn the CEST contrast effect on and off at will, either by turning off the presaturation pulse or by changing its frequency. Following this initial report on diamagnetic small molecules (DIACEST agents) as contrast agents for MRI, numerous reports on different techniques for improving the CEST efficiency of these molecules were reported.52,66,75,78–83 Balaban and co-workers demonstrated that proton exchange processes between endogenous metabolites, like urea and ammonia, and water could be sensitively detected ex vivo and in vivo.65–66 In an interesting account, Aime and co-workers reported the use of a common X-ray contrast agent, iopamidol, as a potential CEST agent.84 Iopamidol has two types of labile protons (amide –NH and alcoholic –OH) capable of exchange with bulk water. The -OH protons exhibit much faster water exchange rates and fall in the fast exchange regime, , and are hence not suitable for CEST. The amide protons however exhibit much slower water exchange kinetics ( estimated at 2560 s−1) and can act as CEST exchange sites. Initial experiments measuring the T1 and T2 relaxation capabilities of iopamidol revealed that while not useful as a T1 agent, it did show promise as a T2-shortening contrast agent. An MR image with a multi-echo type experiment (TE = 120 ms) confirmed that a negative contrast of ~25 % could be detected in the presence of 13 mM agent. Additionally, a phantom image of the same agent at different concentrations revealed that iopamidol can also function as an efficient CEST agent. A comparison of the two phantom images is shown in Figure 6.
Goffeney and coworkers subsequently demonstrated the use of macromolecular cationic polymeric molecules with hundreds of exchanging sites per molecule as a mechanism to enhance the molecular sensitivity of DIACEST agents.78 Cationic polymers have gained significant attention as non-viral DNA delivery systems85 and complexes of cationic polymers with DNA plasmids (anionic polymers), called polyplexes, have been effective in gene therapy. Owing to their importance, non-invasive in vivo monitoring of gene delivery could provide crucial information regarding distribution, local concentration and kinetics. Taking advantage of the large number of amide protons present in macromolecular cationic polymers, Goffeney et al. showed that CEST sensitivity enhancements as large as 500,000 can be obtained compared to monomeric DIACEST agents. For example, poly-L-lysine (PLL, molecular weight of 500 kD) with roughly 2000 amide protons and twice as many amine protons, can result in nearly a 40% change in the signal intensity of bulk water at concentrations as low as 100 μM, a significant improvement over the initial DIACEST molecules proposed by Ward and Balaban. Since this report, other macromolecular DIACEST agents have been explored for use as contrast agents in MRI and have been reviewed elsewhere.12,79,86–87
Glycogen is the readily metabolizable storage form of excess glucose in animals and humans, analogous to starch in plants. It is a large, branched polysaccharide made up of glucose residues linked by α-1,4-glycosidic bonds that can be broken down to yield glucose molecules when needed. Glycogen metabolism in liver plays a critical role in maintaining the blood glucose levels. Excess glucose is converted to glycogen predominantly in liver and muscle cells.88 In the muscle, glycogen provides fuel for muscle contraction whereas in the liver, glycogen is broken down and exported from liver into blood whenever glucose is needed by other organs. While inherited disorders of glycogen metabolism occur in rare cases, glycogen content can also be abnormal in conditions such as obesity, insulin resistance and Type 2 Diabetes (T2D).89–90 Consequently, a non-invasive method to quantify the glycogen content in vivo could have a major impact in understanding the pathophysiology of these common diseases.91 13C MRS can be used as a non-invasive technique in directly detecting the concentration of glycogen in vivo by utilizing the natural abundance 13C 1-C of glycogen, but the requirement of specialized coils and amplifiers on clinical scanners ultimately limits this tool to research sites.
van Zijl and co-workers hypothesized that it should be possible to detect chemical exchange between the hydroxyl groups of glycogen and water using CEST principles and use it as a tool for non-invasive imaging of glycogen in vivo using a standard clinical MRI scanner.91 Conventional 1H NMR spectra of glycogen in water indicated that two potential exchange sites for glycogen, one at 1.2 ppm reflecting the secondary -OH protons and another at 0.7 ppm reflecting the CH2OH group, could potentially be used for CEST (spectrum referenced to water at 0 ppm). The 1.2 ppm -OH resonance is fairly sharp at 4°C but gradually broadens as the temperature is raised to 37°C, indicative of relatively slow exchange at the lower temperature and faster exchange at physiological temperature. Interestingly, the CEST spectrum of glycogen shows only a small CEST exchange peak at 1.2 ppm at 4°C but the effect increases in magnitude at higher temperatures up to 37°C. The CEST spectra recorded at 4.7 T and 9.4 T clearly illustrate the advantage of working in higher magnetic fields owing not to better adherence to the slow exchange condition but to better separation of the resonances from bulk water (larger Δω), thus greatly improving selectivity (Figure 7).
To evaluate the utility of glycogen detection in tissue using CEST principles (glycoCEST), the breakdown of glycogen (glycogenolysis) stimulated by the hormone glucagon was followed in livers isolated from fed mice and perfused with a normal Krebs-Henseleit buffer. The liver of a fed mouse has excess glycogen that is easily degraded and exported for use by other tissues. When the blood glucose levels are low, the hormone glucagon (produced by pancreas), initiates glycogenolysis and the export of glucose. This process can be monitored in perfused livers by following the CEST properties of glycogen before and after injection of glucagon. A series of colorized glycoCEST images of the ex vivo perfused liver as a result of presaturation at ±1 ppm are shown in Figure 8.
The first grayscale image marks the t = 0 perfusion with glucagon. The orange pixels in the colorized images indicate regions of glycogen concentration, suggesting either a non-uniform distribution of glycogen in the liver or variations in proton exchange rates depending on the size or location of the glycogen particle.7 As glucagon perfusion continues over a period of time, the regions of orange pixels substantially reduce and then become constant. These glycoCEST images acquired at regular intervals after administering glucagon clearly illustrate that the magnitude of CEST is reduced following the depletion of glycogen. This technique thus presents an attractive alternative to 13C MRS for following glycogen metabolism in tissues using standard clinical MRI scanners.
The DIACEST examples described above all involve transfer of proton spins but, as we shall see (Section 4.4), exchange of water molecules between multiple pools offers an alternative mechanism for initiating CEST contrast. The resonance frequencies of these pools obviously have to be different and paramagnetic hyperfine shifting lanthanide complexes in which the metal bound water corresponds to one of the pools could potentially be very advantageous because the chemical shift of the lanthanide bound water protons can be very large. The successful application of these complexes for CEST, however, requires that the water exchange rates satisfy the slow to intermediate exchange condition (Eq. 14). For this reason, it is important to briefly review what is known about water molecule exchange in lanthanide complexes. Given the fundamental importance of water exchange in determining the relaxivity of Gd3+-based contrast agents, water exchange in various lanthanide complexes, particularly the Gd3+ chelates, has been studied in great detail.92–94 Variable temperature and variable pressure 17O NMR has proven to be an efficient way to determine the exchange rate and mechanism respectively on various metal complexes. The volume of activation (ΔV‡), defined as the difference between the partial molar volumes of the transition state and the combined partial volumes of the reactants), can be determined from the pressure dependence of the kex by 17O transverse relaxation measurements. ΔV‡ involves the effect of change in bond lengths due to formation and breaking of bonds in the transition state as well as the effect of solvation. The latter is negligible when the ground and transition states have the same charge. In this case, ΔV‡ is a direct measure of the bond formation and breaking occurring in the transition state and provides valuable information about the mechanism of the exchange process. A large positive activation volume is indicative of an expanded transition state of a limiting dissociative mechanism (D) while a large negative activation volume is characteristic for a limiting associative process (A) with a shrunken transition state. Between these two extremes there is a continuous spectrum of interchange (I) mechanisms with smaller positive, zero or negative activation volumes (Id, I, Ia) depending on the extent of contribution of the incoming water molecule. 92–94
The coordination number of lanthanide aqua ions has been the focus of active research for many years. It is now widely accepted that, due to the gradual decrease in ionic radius along the lanthanide series (lanthanide contraction), the lighter (larger) Ln3+- aqua ions (La3+ to Nd3+) are nine coordinate with a tricapped trigonal prism (TTP) geometry while the heavier (smaller) ones (Gd3+ to Lu3+) are eight coordinate with a square antiprism (SAP) geometry.95–98 The Sm3+ and Eu3+ aqua ions exist as an equilibrium mixture of eight and nine coordinate species, presumably due to the particular size requirements of these two species and the ionic radii of these two ions. This change in water coordination number that occurs near Sm3+/Eu3+ influences the mechanism of water exchange. The water exchange rate, kex, for the eight coordinate, heavier lanthanide ions decreases from Gd3+ to Yb3+ and the kinetic data (activation entropy and volume) obtained by 17O NMR relaxation measurements are in accordance with an interchange associative Ia mechanism, which proceeds through a nine coordinate TTP transition state.43,99–103 Unfortunately, for the lighter lanthanides the kinetic effect is not large enough to allow the accurate determination of water exchange rate by 17O NMR measurements, but interchange rate constants with SO42− and the lower limit of kex determined by 17O NMR at 14.7 T for Pr(H2O)9 and Nd(H2O)9 104 suggest that water exchange is the fastest for Eu3+ and gradually decreases towards La3+ (Figure 9).104–105
Water exchange in the lighter lanthanides likely occurs by an interchange dissociative (Id) mechanism with an eight coordinate SAP transition state. The rate of water exchange is determined by the free energy difference between the eight coordinate SAP and the nine coordinate TTP species for a given Ln3+ ion. For the heavier lanthanides, [Ln(H2O)8]3+ has the lowest energy and the energy difference between the eight coordinate ground and nine coordinate transition state increases with decreasing ion size going from Gd3+ to Lu3+ because it becomes more and more difficult to accommodate nine inner sphere water molecules with decreasing size. For the lighter lanthanide ions, the nine coordinate [Ln(H2O)9]3+ is the ground state and the exchange is assumed to occur by a dissociatively activated pathway through an eight coordinate transition state. The water exchange rate is expected to increase with decreasing ionic size from La3+ to Nd3+ as the energy difference between the eight coordinate and nine coordinate geometry increases. In the crossover region (Sm3+ − Eu3+), the energy difference between the eight and nine coordinate species is very low (they exist in equilibrium) so exchange is fastest at these ions.92–94 It is worth noting that a gradual change in the kex along the lanthanide series would be expected if the hydration state of all the Ln3+ ions were the same. However, the observed trend of water exchange rate is a classic example of the so called gadolinium break, which describes an abrupt change in the properties (partial molar volume, heat capacity, molar entropy and extraction behavior) of lanthanide compounds at or near Gd3+.36,106–110 According to a recent extended X-ray absorption fine structure (EXAFS) spectroscopy study,111 all the lanthanide ions exist in aqueous solution as tricapped trigonal prismatic (TTP) species and, contrary to the gadolinium break hypothesis, there is no sudden change in coordination geometry. The water molecules in the TTP capped positions in the lighter lanthanide aqua ions (La3+ − Nd3+) are equidistant but distortions of the regular TTP geometry start to become more and more pronounced from Sm3+ to Dy3+. Gradually increasing water deficiency in the capping positions for the heaviest lanthanides will eventually lead to only one strongly bound capping water molecule for Lu3+. Thus, the gadolinium break is due to increasingly distorted coordination geometry and increasing differences in the Ln3+-O (capping) – Ln3+-O (prismatic) distances.
According to the Solomon-Bloembergen theory, water exchange rate (kex) in the various chelated forms of Gd3+ is one of the determining factors of relaxivity (the others are rotational correlation time τr and electronic relaxation time τe). Theory predicts a kex value of around 108 s−1 for optimal relaxivity. Originally, the water exchange rate of Gd3+-based MRI contrast agents was assumed to be even faster than that of the Gd3+ aqua ion, mainly because multidentate spectator ligands in transition metal complexes speed up the exchange of coordinating solvent molecules. Later, direct water exchange rate measurements performed in Merbach’s laboratory by 17O NMR revealed that water exchange in lanthanide complexes is actually significantly slower than that in the corresponding aqua ion. Current Gd3+-based contrast agents have water exchange rates that are about one order of magnitude slower than the optimal value predicted by the Solomon-Bloembergen theory.20,93 This disparity prompted an extensive worldwide research effort to develop “high relaxivity” agents. As a result, a large number of new Gd3+-complexes have been synthesized and their water exchange kinetics measured. Typically, most ligands for Gd3+ were based on a poly-aza poly-carboxylate framework and, in virtually all cases, water exchange was found to be slower than that required for optimal relaxivity; so a major effort was directed at accelerating water exchange in these complexes toward the optimal value. With the advent of PARACEST agents, the focus of research attention shifted to the design of lanthanide complexes that exhibit relatively slow water exchange. Thus, a considerable amount of knowledge has been gained about the structures of various lanthanide complexes and the rate at which water exchanges in those complexes as well the various factors that affect the water exchange rate. These will be discussed in detail in the following sections.
Lanthanide chelates with octadentate poly-amino poly-carboxylate ligands (including all commercial Gd3+-based contrast agents) generally have one inner sphere water molecule. Activation volume measurements on Gd 3+-1 (1 is commonly referred to as DTPA: diethylenetriamine pentaacetic acid) and Gd3+-2 (2 is commonly referred to as DOTA: 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) indicate a limiting dissociative (D) mechanism in which the incoming water molecule does not participate in the exchange.112 This results in an increase in activation energy and consequently slower exchange rates compared to the Gd3+-aqua ion, where the water exchange proceeds by an interchange mechanism (I). This is not surprising because the cage of donor atoms surrounding the Gd3+-ion is expected to prevent the entering water molecule from assisting the leaving one (ten coordinate transition state is not possible). If the coordination number is lower than 9, then the incoming water molecule can assist the outgoing molecule resulting in a faster exchange rate. This effect is nicely illustrated by a comparison of the exchange rates for two related complexes, the nine coordinate Gd3+-3 (τm = 70 ns) with two bound water molecules and the eight coordinate Gd3+-4 (τm = 6 ns).113 The latter has only one metal bound water molecule resulting in a water exchange rate two orders of magnitude faster than the nine-coordinate complex. It is assumed that the two complexes have different exchange mechanisms: water exchange on the eight-coordinate complex proceeds by an associative mechanism because the transition state can easily accommodate another H2O molecule, which is not possible for the nine-coordinate species.
As expected, for a dissociative water exchange mechanism (D), increased steric crowding around the water binding site will accelerate water exchange because the coordinated water molecule must leave before the incoming molecule can bind. Insertion of an extra methylene group into the polyamine backbone or the pendant arm of the ligand is a simple way to increase the steric compression around the binding site.114–117 For example, addition of a methylene group to the macrocyclic amine backbone of DOTA gives rise to the ligand 5 whose Gd3+ complex has a water exchange rate almost two orders of magnitude faster than Gd3+-2.114–115 Interestingly, the Gd3+-O (bound water) distances are quite similar in these complexes (0.248 nm in Gd3+-5 and 0.245 nm in Gd3+-2) so the faster exchange in Gd3+-5 is not due to an increased Gd3+-O bond distance but rather has been traced to the extra methylene group which shifts the plane defined by the four carboxylate oxygens closer to the bound water, thereby hastening dissociation.114 The small negative value of activation volume for the exchange process indicates a strong participation of the incoming water molecule (interchange mechanism).115 A further increase in steric compression around the water coordination site by adding more methylene groups results in complexes which do not contain an inner sphere water molecule. An example of this is given by Gd3+-6.112,118 Introducing a methylene group into DTPA has a similar effect: both Gd3+-7 and Gd3+-8 have much faster water exchange rates than Gd3+-1, while Gd3+-9 does not contain a inner sphere water molecule.115–117 A comparison of these complexes suggests that backbone modification should have a stronger effect by about an order of magnitude. Substituting the bulkier phosphonate groups for acetate is another approach to increase the steric crowding around the water binding site. Thus, replacing one acetate sidearm of DOTA with a phosphonate group as in Gd3+-10 results in a significant increase in the rate of water exchange.119–120 Replacing more acetates with phosponates again sterically blocks the inner sphere water coordination position such that lanthanide complexes of 11 and 12 do not have a metal bound water.120–121
While these examples demonstrate that it is possible to alter the water exchange rate of DOTA and DTPA-like complexes by adding a methylene group to the ligand or replacing the acetates with other functional groups, the resulting complexes usually have decreased thermodynamic stability and/or kinetic inertness compared to the corresponding parent ligand.115,122 However, lanthanide complexes of DOTA and its derivatives exhibit a very unique type of coordination isomerism which allows more subtle control of the steric compression around the water coordination site. The relative orientation of the acetate sidearms and the macrocyclic ring results in two diastereomeric coordination isomers that differ in the twist angle between the basal N4 and capped O4 squares. A N4/O4 twist angle of around 39° defines a square antiprismatic (SAP) geometry, while −29° defines the twisted square antriprismatic (TSAP) coordination geometry (Figure 10). Complexes of DOTA and DOTA-like ligands have two independent components of chirality originating from the conformation of the five-membered N-C-C-N chelate ring and the helicity of the side-arms. Depending on the sign of the N-C-C-N torsion angle, the conformation of each ethylene group in the macrocyclic ring is left-handed, designated as λ (negative N-C-C-N torsion angle) or right-handed, δ (positive N-C-C-N torsion angle), and so the macrocyclic ring can have either (λλλλ) or (δδδδ) conformation. Similarly, the orientation of the pendant arms can be clockwise Δ (positive N-C-C-O torsion angle) or counter-clockwise Λ (negative N-C-C-O torsion angle). Hence, there are two enantiomeric pairs of diastereoisomers: the Δ (λλλλ) and Λ (δδδδ) enantiomeric pair adopts the SAP geometry while the Λ(λλλλ) and Δ(δδδδ) enantiomeric pair adopts the TSAP geometry.123–124 These coordination isomers can interconvert by arm rotation (Λ ↔ Δ) or ring inversion [(λλλλ) ↔ (δδδδ)]. The SAP/TSAP ratio can easily be estimated from the 1H NMR spectra because the lanthanide induced shifts observed for the ligand protons (and the bound water protons, when observable) are substantially larger in the SAP isomer. Interestingly, the water exchange rates in the TSAP isomers are about two orders of magnitude faster than in the SAP isomers because the latter has a more compact structure with less steric crowding around the bound water. 125–128 In fact, the Gd3+-bound water lifetime of the TSAP isomer of Gd3+-2 is close to the optimal value predicted by Solomon Bloembergen Morgen theory. However, the population of the two species is strongly influenced by steric factors including the size of the lanthanide ion. As expected, the lighter (larger) ions prefer the TSAP geometry while the heavier (smaller) lanthanide ions favor the SAP configuration. 129 To control the water exchange rate across the lanthanide series, one must be able to select one of the coordination isomers independently from the size of the lanthanide. This was elegantly achieved with DOTA-like ligands by freezing out the molecular motions through which the inter-conversion occurs. In lanthanide complexes of ligand 13, chiral substitution of the alpha carbon of all four acetate sidearms with methyl groups inhibits the arm rotation, while the nitrobenzyl group on the macrocyclic backbone eliminates the ring inversion. Since the orientation of the arms and the conformation of the macrocyclic ring is determined by the chirality of the asymmetric carbons present in the ligand (the alpha carbons of sidearm on the carbon bearing the nitrobenzyl substituent), appropriate selection of the chirality will control the coordination geometry, and thus the water exchange rate.130 Lanthanide complexes of the bridged ligand 14 represent an example of locking the conformation of the ligand without the use of chiral carbon. The original intent of this work was to slow water exchange by hindering the departure and approach of water molecules to the lanthanide ion, yet it was found that the bridge was so rigid that water was entirely excluded from the Ln3+ coordination sphere. The crystal structure of the Yb3+ complex revealed that it exists exclusively as the TSAP isomer. This unexpected geometry for Yb3+ is due to the presence of the bridging unit, which forces the same orientation on the conformation of the macrocyclic ring and the helicity of the pendent arms affording only the two TSAP enantiomeric pairs, Λ(λλλλ) and Δ(δδδδ). The “conformation locking” effect of the bridge can be understood by considering the two ethylene bridges that are shared by the macrocycles (the 12- and 21-membered) of the bicyclic system, which will achieve the lowest energy when the two rings adopt the same conformation (λ or δ). In addition, two of the amide side-arms are part of the 21-membered ring and so their helicity will be determined by the conformation of the macrocyclic backbone. Interestingly, the interconversion between the two enantiomers [Λ(λλλλ) ↔ Δ(δδδδ))] must occur by sequential arm rotation and ring inversion through the SAP isomer.
The water exchange rates of lanthanide complexes of both DTPA and DOTA amides show a strong dependence on the lanthanide ionic radii.131 The exchange rates in the lanthanide complexes of 15 is nearly constant for the lighter lanthanides (Ln = Nd3+, Eu3+), but increases sharply for the heavier ones from Gd3+ to Ho3+ (Figure 11a). This trend was explained by a change in the mechanism near the center of the lanthanide series; activation volume of exchange (Figure 11b) was consistent with an interchange mechanism (I) for the Nd3+ complex, while water exchange occurs by a dissociative mechanism (D) for the remaining lanthanide complexes. Water exchange in the LnDOTA-tetra-glycinate complexes (ligand 16) shows an even more surprising trend, opposite that observed for the aqua ions (Figure 12).132 Based on the measured bound water lifetimes (τM), the complexes can be divided into three distinct groups. The lifetime increases from Pr3+ to Eu3+, decreases dramatically between Eu3+ and Tb3+, then decreases very slightly between Tb3+ and Yb3+. The high resolution NMR spectra report that the entire series of Ln3+-16 complexes exist in solution largely as the SAP isomers, so the observed trend of kex is not due to a switch in coordination geometry but rather is more likely due to a change (or changes) in exchange mechanism that occurs along this series. The activation volumes of exchange measured for Eu3+-17 and Eu3+-18 (+ 4.9 cm3/mol and +6.9 cm3/mol, respectively) indicates that water exchange in these two DOTA-tetraamide systems occurs via a dissociative activation mechanism.128 Assuming the same exchange mechanism is operative for the remaining heavier lanthanide complexes of these DOTA-tetraamide ligands, one would predict an increase in the Ln3+-O (bound water) distance with decreasing ionic radii (from Eu3+ to Lu3+) and hence faster water exchange as observed for the complexes with ligand 16 (Figure 12). For the lighter Ln3+ ions water exchange slows from Sm3+-16 to Eu3+-16 and this likely reflects at least a partial contribution from an associative mechanism whereby the incoming water molecule assists the leaving molecule. The fact that water exchange is so finely tuned with a minimum at Eu3+ has important implications for the design of paramagnetic magnetization transfer agents. While the hyperfine shifts induced by Eu3+ are smaller than the ones produced by some other lanthanides (Tb3+, Dy3+ and Tm3+), the slowest water exchange rate of Eu3+ combined with its low paramagnetic relaxation enhancement (T1 and T2 shortening effect) makes Eu3+ the paramagnetic ion of choice in the design of PARACEST agents involving presaturation of the bound water molecule.
A recent study of the positively charged complex, Gd3+-17, with different counter ions X (X = OTf−, Cl−, OAc−, Br− and I−) demonstrated that the counter ion has a surprisingly big effect on water exchange (Figure 13).133 It was proposed that anions modulate water exchange by altering the structural ordering or the extent of the second sphere hydration - the more ordered the second sphere hydration, the slower the water exchange rate. This was supported by the X-ray crystal structures of the complexes where the degree of structural order and local hydration decreased from Cl− to I−, consistent with faster water exchange. Interestingly, no significant anion effect was found for the more hydrophobic (SSSS)-19, which was explained by the overriding effect of the hydrophobic side-arms.
Ultimately, the water exchange rate in lanthanide complexes is determined by the electron deficiency around the central lanthanide ion. Qualitatively, greater electron deficiency on the lanthanide ion increases the electrostatic attraction between the metal ion and the bound water resulting in a shorter Ln3+-O (bound water) bond distance and slower water exchange. If Gd3+-2 is considered as the reference complex, ligating groups that render the central metal ion more electron-deficient compared to carboxylates are expected to decelerate water exchange. Accordingly, neutral coordinating groups that are less capable of electron donation to the central lanthanide ion compared to the negatively charged carboxylate should decelerate water exchange. The oxygen of the amide functionality is a nice example of such a donor atom. The weak donor strength of amides compared to carboxylates correlates well with the basicity of these functional groups. The pKa of the conjugate acid of amides is around −0.5 while the pKa of carboxylic acids is roughly 4.134 These and other properties of the amide group (restricted rotation about the C-N bond, short C-N bond length, stability toward hydrolysis and protonation at the oxygen rather than nitrogen) can adequately be explained by the structure shown in Figure 14, in which the contribution of the C=N resonance form is around 30% even though there remains considerable ongoing debate about the validity of this description.135 The decelerating effect of a coordinating amide to water exchange in DOTA and DTPA-amide complexes is surprisingly strong. This is nicely illustrated by a series of both DOTA- and DTPA-based Gd3+-complexes (Gd3+-20, Gd3+-21, Gd3+-22, and Gd3+-23-Gd3+-34) in which progressive substitution of an acetate side-arm by an amide side-arm slows water exchange substantially with each substitution. As a general rule of thumb, substitution of an amide for a carboxylate results in about 3 to 4 fold decrease in kex.5,136 In the case of PCTA derivatives, which have two bound water molecules, the effect of amide substitution is significantly weaker, as the water exchange rate of Gd3+-35 is only about half of that for Gd3+-3.137
High resolution 1H NMR studies indicate that lanthanide complexes of DOTA-tetraamide ligands usually adopt the SAP geometry in solution when the amide side-arms are not sterically demanding. However, increasing steric requirements of the side-arms favor the TSAP isomer with bulky primary amines or secondary amines showing a marked preference for the TSAP geometry. Alpha substitution on the acetate sidearms also increases the steric crowding around the central lanthanide ion and DOTA-tetraamide complexes such as 36 exist in solution largely as the TSAP isomer.125,128,138 The crystal structure of various LnDOTA-tetraamide complexes indicates a preference for SAP geometry in the solid state.125,139–142 The crystal structures of the lanthanide complexes with ligand 19 show that all complexes adopt the SAP geometry in solid state.139 As a consequence of increased steric crowding around the Ln3+-ion with decreasing ionic radius, the size of the O4 square and O-Ln3+-O angle decreases while the Ln3+-O (bound water) distance gradually increases the along the series. Although the distance between the N4-O4 squares seems to be relatively invariant across the series, the heavier, smaller Ln3+ ions are positioned closer to the N4 plane. NMR studies confirmed that the Eu3+ complex adopts the same SAP geometry in solution yet, interestingly, the observed 1H NMR shift values for the Yb3+ complex show that it adopts the TSAP geometry in solution with one Yb3+-bound water molecule.143 The SAP isomer could not be detected in solution but clearly it must be present in very low concentration because the complex crystallized in this form. The preferential crystallization of the minor isomer in this case clearly demonstrates that the crystal structure is not always representative of the solution structure.
The crystal structure of Eu3+-19 triflate indicates that this complex also adopts the typical SAP geometry but the Eu3+-O (bound water) distance is slightly shorter than usual.139 Interestingly, the four pendant arms are not equivalent as one of them is positioned perpendicular to the Ln3+-O (bound water) bond while the other three occupy an equatorial position. The solution structure of Yb3+-16 was studied by hyperfine shift analysis using the SHIFT ANALYSIS method and a comparison of the experimental and predicted LIS values revealed that in acetonitrile the side-arms are positioned axially, while in water, which is more polar than acetonitrile, an equatorial arrangement is preferred. The dependence of the average position of the side-arms on the solvent may explain why the exchange rate is slower in acetonitrile.42
Recent studies of EuDOTA-tetraamide derivatives revealed that the water exchange rate can be further modulated by electronic effects and the degree of hydrophobicity, polarity and steric bulk of the substituent attached to the amide nitrogen atoms of the ligand. The available experimental data strongly suggest that a complex interplay of each of these factors will significantly affect water exchange and should be considered in the design of new PARACEST agents. Modulation of water exchange through electronic effects was first demonstrated with the Gd3+-complexes of p-nitrobenzyl- and p-aminobenzyl-DO3A (37). In the nitro derivative, the electron deficiency around the Gd 3+-ion is greater than in the p-amino derivative and this stabilizes the nine-coordinate ground state of the complex resulting in slower water exchange.144 Later it was shown that inductive and mesomeric effects of remote substituents can also affect water exchange rates in a series of EuDOTA-tetraamide (38) complexes provided the substituent can electronically communicate with the donor atoms.145 Although the largest difference in the water exchange rates in this series of complexes was only about two fold, the kex values showed a good correlation with the electron donating/withdrawing effects of the substituent. These data suggest that modulation of the water exchange rates by electronic effects is a feasible approach to the design of responsive PARACEST agents.145 Intuitively, one would expect that bulky, hydrophobic substituents on the amide would slow down water exchange by restricting access of water to the metal ion. Indeed, a good correlation between the kex and the calculated solvent accessible surface area was found for a series of EuDOTA-tetraamides with different amide substituents (Figure 15). This effect was attributed to altered organization of the second hydration sphere, which increases the free energy of the transition state for the exchange.146 However, recent data on other DOTA-tetraamide derivatives with varying hydrophobicity and polarity (43–48) strongly suggest that “smaller the accessible surface area, slower the water exchange rate” rule is not always valid and other factors are also involved. In addition to hydrophobic effects, another important factor is polarity. In fact, it was shown that introduction of polar substituents capable of participating in hydrogen-bonding interactions with second sphere water molecules slows down water exchange, stabilizing the second coordination hydration sphere which, in turn, stabilizes the single inner-sphere water molecule.138
Water molecule exchange by definition involves the exchange of oxygen atoms and therefore the water exchange rate represents the lower limit of the prototropic exchange rate. Obviously, the exchange of protons can occur independently from the coordinated water but this requires an acid or base catalyst. So, it is normally assumed that the prototropic exchange rate is equal to the water exchange rate over a wide range of pH, about 3–8.32,147 However, certain functional groups on a ligand (as in 22 and 49) can act as a catalyst to speed up proton exchange.148–149 For example, the residence lifetime of the bound water in Gd3+-22 (26 μs as determined by 17O NMR) is not pH dependent but the water 1H relaxivity of the complex shows strong pH dependence between pH 5 – 8. The pH responsive behavior of this complex is due to catalysis of prototropic exchange by the extended phosphonate groups that undergo protonation in this pH range.
Although numerous experiments involving diamagnetic agents have been performed using CEST techniques, they are faced with one major drawback. A CEST experiment is bound by the condition that the exchangeable proton needs to be selectively irradiated while maintaining the intermediate-to-slow exchange condition (kex ≤ Δω). Accordingly, the frequency difference (Δω) between the two exchanging pools needs to be sufficiently large. Since Δω values are typically no more than ~6 ppm, very slow exchange rates are required to satisfy the intermediate-to-slow exchange condition. While the use of higher magnetic fields is obviously advantageous as discussed earlier (Figure 7a), the practical feasibility of moving to a high enough field to meet this requirement may be problematic. Another intrinsic problem associated with using DIACEST agents in vivo is that it may be difficult to distinguish the CEST effects of the agent from the inherent tissue MT effect. Therefore, a large Δω value is clearly advantageous for CEST not only because it would permit the use of more rapidly exchanging species, but it may allow for easy differentiation of CEST from tissue MT.
While it has been known for more than 40 years that some of the paramagnetic lanthanides are able to induce extremely large (up to several hundred ppm) hyperfine shifts in nearby nuclei, it was only relatively recently realized that the exchange rate of a Ln3+-bound water molecule in certain lanthanide complexes may be slow enough to satisfy the kex ≤ Δω condition for CEST. The first generation of such compounds consists of lanthanide complexes with a variety of DOTA-tetraamide ligands. Such systems typically consist of one Ln3+-bound water molecule and –NH groups. The Ln3+-bound water molecule is in a capping position and aligned with the principle magnetic axis of the complex and, for this reason, experiences the largest hyperfine shifts, while the amide protons are situated somewhat further away from the Ln3+ ion and are oriented in the direction where the dipolar field of the Ln3+ ion nearly vanishes (Figure 16). Consequently, the amide proton shifts are roughly an order of magnitude smaller than those observed for the Ln3+-bound water molecule (Refer Section 3, Figure 2 for a more detailed explanation on lanthanide-induced shifts). In breakthrough experiments that set the stage for PARACEST, two independent research groups concurrently discovered that the water molecule in Eu(III) complexes of certain DOTA-tetraamide ligands (16, 17, and 18) could be observed directly using 1H and 17O NMR.42,126,150
The separate resonance for bound water cannot be observed in the NMR spectra of Ln3+-DOTA or Ln3+-DTPA because water exchange is too fast in these poly-carboxylate systems to be seen on the NMR timescale. Lanthanide complexes of the DOTA-tetraamide ligands, on the other hand, often display slow enough water exchange to allow direct detection of the bound water molecule by 1H and 17O NMR.42,125–126,128,143,150–152 As one would expect, the conditions necessary for detection depend on the exchange rate and the paramagnetic line broadening properties of the lanthanide ion itself. Often, the NMR spectra are recorded in wet acetonitrile at temperatures ranging from 233 K to 333 K,126,150,152 but, in some cases, when the exchange is slow enough, the bound water can be observed in pure water at room temperature.42 Based on the 1H NMR data of several Ln3+-DOTA-tetraamide complexes, an interesting and quite useful empirical relationship has been established between the 1H chemical shift of the bound water and the axial H4 macrocyclic proton; the shift of the bound water proton resonance is about 2-fold larger than that of the H4 proton.152 This relationship is very helpful in estimating the chemical shift of the bound water protons even when they cannot be directly observed. The kinetic parameters characterizing the exchange process (kex, activation volume activation enthalpy and entropy) in such systems can be determined from variable temperature and variable pressure 1H and 17O NMR measurements.126,128,150 Separate Eu3+-bound water resonances were detected for the SAP and TSAP coordination isomers in Eu3+-17 and as anticipated, the hyperfine shifts of the two bound water resonances were quite different with the SAP isomer displaying a much larger hyperfine shift than the TSAP isomer. In acetonitrile and at low temperatures (232 K), the bound water 1H and 17O NMR resonances were identified at 1H δ = 84.1 ppm, 17O δ ≈ −1150 ppm (SAP isomer) and 1H δ = 19 ppm, 17O δ ≈ −950 ppm (TSAP isomer), respectively.126,150 Variable temperature and pressure 1H and 17O NMR studies of the two isomers not only afforded the kex values of each isomer but also provided direct insight into the correlation between the water exchange and SAP-TSAP interconversion. Simultaneous fitting of the NMR data for water exchange and the SAP/TSAP interconversion revealed that both processes proceed through a common intermediate that does not have an inner sphere water molecule, with interconversion of isomers occurring predominantly though arm rotation.126 Subsequent research after the first reports of direct NMR observation of the metal bound water soon led to the application of slowly exchanging lanthanide complexes in CEST experiments.
A high resolution 1H NMR spectrum of Eu3+-16 in water at room temperature is illustrated in Figure 17. The relatively broad peak near 50 ppm that integrates as two protons and disappears upon addition of a drop of D2O to this sample can be assigned to a single, inner-sphere Eu3+-bound water molecule. RF saturation of this highly shifted peak prior to collection of the NMR spectrum resulted in a 20 % reduction in signal intensity at concentrations of agent as low as 10 mM.1,151 If one performs this presaturation experiment over a wide range of saturation frequencies, one obtains the CEST spectrum shown in Figure 18. While the Eu3+ complexes of various DOTA-tetraamide complexes typically display a Eu3+-bound water resonance around 50±10 ppm depending upon temperature, other lanthanide ions show even larger hyperfine shifted water resonances, ranging from +500 ppm to −720 ppm (Table 3). The CEST spectrum of Dy3+-17 shows a broad CEST peak at −720 ppm from the bound water protons (Figure 19) owing to the fact that Dy3+ induces about 25-fold larger pseudocontact shifts than Eu3+ (see Table 1).4 These examples illustrate that exchange sites covering a wide range of Δω values are potentially available for CEST activation. This expands the number of possible exchanging species that can potentially meet the CEST requirement to include functional groups that normally would exchange too fast in a diamagnetic molecule. As discussed in detail below, the rate at which a Ln3+-bound water molecule exchanges with bulk solvent is much more sensitive to subtle changes in the structure of the complex and the chemical environment than amide protons, making water exchange a versatile platform for the design of responsive PARACEST agents.
It is interesting to note that CEST activation via an exchanging water molecule has only been detected in lanthanide complexes (mostly Eu3+) of DOTA-type ligands. As indicated earlier (Section 126.96.36.199), Ln3+ complexes of DTPA pentaamide derivatives like 30 and 31 display water exchange that is too fast for CEST. The reason for this, likely, lies in the unique properties of DOTA and its derivatives. The physico-chemical properties of the complex (thermodynamic stability, kinetic behavior, magnitude of the induced shifts and spectroscopic properties) are determined by its molecular structure. The macrocyclic DOTA ligand is more rigid, preorganized and symmetric than the acyclic DTPA and it is quite plausible that the higher flexibility of DTPA-amide derivatives may allow faster water exchange by a different mechanism. The effect of rigidity has not been studied in detail for DOTA-tetraamide derivatives, although there is some indication that rigidifying the macrocyclic ring may actually accelerate the exchange since the water exchange on the more rigid, backbone substituted complex Gd3+-50 is about 20 times faster than on Gd3+-22.153
Although not apparent in the CEST spectrum of Eu 3+-16, the protons of the coordinated water molecule are not the only exchanging site in this complex. The amide protons of the complex also show significant reduction in water intensity, as can be seen from the CEST spectrum of the analogous acid Eu3+-51 (Figure 20). While it is difficult to observe the NH-based CEST effect clearly in Eu 3+-based PARACEST agents, other lanthanide complexes induce much larger shifts. From a comprehensive report by Terreno and co-workers of the amide-based PARACEST effects across the lanthanide series for the ligand 51 along (Table 3), it is evident that most of the lanthanide complexes in the series satisfy the exchange condition required to observe CEST resulting from amide proton exchange.154 Another classic example of amide CEST can be illustrated by considering the Yb3+ complex of 17.155 The 1H NMR spectrum of the complex is dominated by the presence of two non-equivalent amide protons of each -NH2 group at −14.5 ppm and −17.7 ppm, with both protons being highly sensitive to pH (Figure 21).
ka, k b and kc represent the rate constants due to acid-catalyzed exchange (ka), base-catalyzed exchange (kb) or other possible contributions (kc). The rate determining step for amide protons is deprotonation under base-catalyzed conditions and protonation under acid-catalyzed conditions. Non-catalyzed proton exchange with water is thought to be negligible for the amide protons.158–159 The pH dependence of these NH-based exchange mechanisms has been studied in detail for many years as a means to understand protein structures in solution – these have been reviewed elsewhere in detail.156–157,160–161 It has been shown that amide proton exchange is largely base-catalyzed at pH > 5 so Eq 31 can be reduced to Eq 32.
Consequently, while NH-based PARACEST effects are not as strongly influenced by the lanthanide induced shifts when compared to bound-water PARACEST effects, they may offer some advantages owing to their innate pH-sensitivity and their slower exchange kinetics.
Woods and co-workers were able to identify another potential exchange group in lanthanide complexes of DOTA-type ligands when they replaced the coordinating amide side groups with hydroxyethyl arms (ligands 52 and 53).162 It was shown earlier that the hydroxyethyl arms retain their protons even upon coordination of the –OH groups to a lanthanide ion leaving them open for exchange with bulk water.163–165 However, high-resolution 1H NMR spectra also show that the coordinated –OH protons are not highly shifted in this complex and consequently –OH proton exchange is fast compared to Δω. However, with similar complexes (Eu-53) dissolved in slightly wet acetonitrile, CEST exchange peaks were detected for all three hydroxyl groups plus the Eu3+-bound water molecule (Figure 22). As more water was added to this system, the rate of proton exchange with the water pool gradually increased and the CEST effects completely disappeared in pure water. While structurally significant and an interesting example of –OH based CEST, these ligands need to be further modified to provide stronger lanthanide induced shifts and slower water exchange rates for practical applications.
One of the most important advantages of PARACEST agents over their diamagnetic counterparts lies in the fact that faster exchanging systems can be used to generate CEST effects due to the large Δω values induced by the paramagnetic lanthanide ions. While this is true, it is clear from the modified Bloch equations for chemical exchange (Section 4.1.2) that other factors such as relaxation rates, applied power and duration of the presaturation pulse also play an important role. For PARACEST agents to be useful in a clinical setting, one of the critical considerations is the specific absorption rate (SAR) power limitations (B1) of any future human application. There is general agreement among users that a 3 to 5 % change in the water signal intensity is sufficient to see contrast in an MR image. In order to understand the parameters that need to be optimized to maximize CEST effects for clinical applications, Woods et al. used a previously published two-site exchange model to make several important predictions.4 Eq 33 (a rearranged form of Eq 27) predicts that a maximum CEST effect (largest decrease in bulk water intensity smallest Mza/Moa) will be observed for a PARACEST agent whenever the water residence lifetime of a pool of exchangeable protons (τb) is as small as possible and the concentration of the contrast agent ( ) and bulk water relaxation (T1a) remain as large as possible. The observed CEST effect will also depend on the extent of saturation of the PARACEST agent (pool B) - incomplete saturation of the pool B protons, determined mainly by the power and duration of the presaturation pulse B1, will result in less than optimal CEST.
Although CEST increases with higher B1 values, this value needs to be chosen to achieve the right balance between optimum CEST and minimum indirect irradiation of the bulk water pool of protons (pool A) while still staying within the limitations imposed by SAR for irradiation power.67 Thus, the chemical parameters that can be optimized for optimal CEST include the residence lifetime of the protons on the contrast agent, τb, the chemical shift of the exchangeable protons, Δω, and the bulk paramagnetic relaxation properties of the contrast agent itself. A numerical analysis presented by Woessner et al. predicts quite simply that the optimal τb depends on the applied B1 (Eq 34).67 Simulated two-site exchange CEST spectra show that for a given value of B1, the residual magnetization of the pool A protons passes through a minimum (maximum CEST effect) when Eq 34 holds true.
This important relationship, shown graphically in Figure 23 allows us to predict an optimal value of τb for any a given value of B1. Thus, from a group of PARACEST agents with known water exchange kinetics, we should be able to select the optimal agent for any given B1 restriction (determined by the experimental setup, i.e., instrument, coil, SAR limitations). Numerous chemical modifications to modulate water exchange rates have already been discussed in detail earlier in this review. As an example, a recent work by Mani et al. elegantly demonstrated this prediction that simple structural modifications to the ligand framework allows us to modulate the water exchange rates for a series of Eu(III)-based PARACEST agents.138
One of the primary advantages of CEST as a contrast mechanism is that the effect is detectable only when a selective RF irradiation pulse is applied at the specific frequency of the exchangeable protons. Consequently, if the resonance frequencies of the exchangeable protons are well separated, multiple agents can be used in one imaging experiment and selectively turned on and off by applying the appropriate frequency-selective RF pulse. PARACEST agents based on Ln3+ ions are particularly attractive in this respect because a variety of lanthanide ions across the series can be used to generate markedly diverse complexes with exchangeable protons covering a wide range of frequencies. The proof of principle for multi-frequency MRI was demonstrated by Aime and co-workers in phantoms and living cells with Eu3+ and Tb3+-based PARACEST agents.166 Phantom images of the two complexes, Eu3+-51 and Tb3+ -51, both individually and as a cocktail mixture clearly presented the selectivity of the PARACEST agents; while the Eu3+ complex was irradiated only at the resonance frequency of its bound water proton at 50 ppm, the Tb3+ complex was selectively irradiated at −600 ppm (Figure 24).
“Multi-Frequency” MRI using PARACEST agents would be extremely useful for imaging multiple biomarkers simultaneously, since any parameter (pH, temperature, redox, metabolite concentration) that alters the water exchange kinetics will in turn alter image contrast.145,167–170 To date, however, with a few exceptions, only Eu3+ complexes have displayed sufficiently slow water exchange kinetics to satisfy the kex ≤ Δω requirement in pure water. Complexes of other lanthanide ions that induce large paramagnetic shifts [Tm 3+, Tb3+, Pr3+ and Yb3+] have faster than optimal water-exchange kinetics and usually do not show CEST from bound water. This has limited the development of multi- frequency PARACEST agents based on water exchange. An alternate approach to overcome this problem was shown by Aime et al. using LIPOCEST agents (Section 4.4.8).171–172 They demonstrated that encapsulating different neutral polynuclear Tm3+-based shift reagents in the inner core of liposomes produces significant changes in the frequencies of the water protons encapsulated within the core.
In a recent account, Viswanathan et al. demonstrated that the resonance frequency of the Eu3+-bound water is sensitive to relatively minor structural changes in the periphery of the ligand (54–56).173 All three Eu3+ complexes used in this study exhibited three distinct, well separated resonance frequencies for the Eu3+-bound water protons: 45 ppm for Eu3+-54, 54 ppm for Eu3+-55 and 64 ppm for Eu3+-56. In a mixture containing all three complexes, each agent could selectively be activated without interference from the others by applying the appropriate presaturation frequency to produce multi-frequency PARACEST images (Figure 25).
Several methods, virtually all NMR based, have been used for the measurement of water exchange rates of lanthanide complexes. Merbach and co-workers pioneered the use of 17O variable temperature and pressure experiments for the measurement of kinetic and activation parameters that describe the water exchange processes in lanthanide macrocyclic complexes.93,112 In relatively fast exchanging systems, the bound water signal cannot be observed in the 17O spectra. In such cases, the water exchange rates for the coordinated water molecule can be directly obtained from the variable temperature measurements of the transverse 17O relaxation rates of the bulk water using the Swift Connick approach, as reviewed elsewhere.92–94,102,115,126,131,174–175 This technique is applicable to the determination of kex values over a wide range except for extremely fast or slow systems. It is especially useful for Ln3+ complexes where the relaxation enhancement of water is relatively large (Gd3+, Tb3+, and Dy3+). If the complex exists as two interconverting isomers, the bulk water is obviously affected by both isomers present in the solution, the kex thus obtained contains contributions from both isomers and represents a weighted average of the individual kex values (Eq 35):125,127
where P is the molar fraction of each isomer.126
Under slow exchange conditions when the Ln3+-bound water 1H and 17O NMR signals are observable, the exchange rate can be determined by the temperature dependent line shape analysis of either or both of the 1H or 17O resonances using a simple two-site chemical exchange model with unequal populations.42,126,128 This method has been used to measure the exchange rate of several LnDOTA-tetraamide complexes in acetonitrile at low temperatures and in some cases, when the bound water peaks due to both the SAP and TSAP isomers are observable, the kex value of each isomer can be determined separately.
Other techniques based on NMRD profiles, variable temperature 1H NMR relaxation time and luminescence measurements have also been used to measure water exchange rates in lanthanide complexes but the scope of these methods are more limited. 176–177 In this section, we concentrate on two methods in particular, that have wide applications for PARACEST agents - numerical solutions to the modified Bloch equations and Omega plots.
Most PARACEST agents typically have more than one type of exchanging site, protons from the amide group and bound water being the most common sites encountered. In order to measure the exchange rates for such a three pool system, Woessner and co-workers extended the modified Bloch equations described earlier to a three-pool case to include exchange for a three-pool system between pool A (bulk water) and pool B (bound water) and between pool A (bulk water) and pool C (-NH). 67 It is generally assumed that there is no exchange between the pools B and C (Eqs 36–44).67
Cb and Cc represent the transition rates of the protons leaving pools B and C and entering pool A. Two new terms (Cab and Cac) have been introduced in the three pool system when compared to the two pool system described earlier (Section 4.1.2). Cab is the transition rate of the pool A protons entering pool B and Cac is the transition rate of the pool A protons entering pool C. All the other definitions remain similar to the two-pool case. Similar to the previous descriptions for the two pool case, the thermal equilibrium Z magnetizations ( and ) are directly proportional to the number of protons in each pool (Eqs 45–47).
At steady state, with the time derivatives of all nuclear magnetizations at zero, solutions to the resulting set of nine equations gives . Nevertheless, solutions to these equations without algebraic simplification for a three pool system or even more complicated systems are painstaking. To obtain a simplified solution without approximations, Woessner et al. applied Cramer’s rule to numerically solve the modified Bloch equations for the steady state case. Cramer’s rule simply involves writing the augmented matrix that contains the coefficients of all the nuclear magnetizations and constant terms. The numerical X, Y, and Z magnetizations for all the pools are included in the matrix shown in Eq 48. Using commercial software such as MATLAB or PSI-PLOT, these matrix equations can be re-written as a linear equation of the form dY/dt = A * Y + b and solved by using the matrix left division (\) and matrix power (expm 1) functions of MATLAB. One of the most important advantages of this method lies in the fact that it does not require an assumption that the system is at steady-state. Instead, the matrix solution to the modified Bloch equations accounts for the length of the presaturation pulse by the use of an exponential operator that expresses the time dependence of the equations.
Fitting experimental CEST spectra to the modified Bloch equations enable one to obtain water exchange rates for any number of complicated exchange processes involving multiple exchanging pools for a wide variety of situations, without approximations.82–83
In a recent report,178 Dixon and co-workers detailed a new concentration-independent method to determine the water exchange rates of CEST contrast agents. Consider first the direct saturation of a single pool of bulk water protons (pool A) without chemical exchange, as discussed by Woessner et al. (Eq 49).67
If this pool encounters chemical exchange with another pool during this process, an additional term, τa (the lifetime of the protons in pool A) is introduced into Eq 49.
Now consider the case of PARACEST agents, where bound water pool protons are in chemical exchange with protons of the bulk water. After prolonged saturation of the bound water protons, such that they achieve steady state, the net magnetization remaining in the bulk water pool is determined by the loss of magnetization of the pool B protons and the rate of recovery of the bulk water magnetization due to T1a. Thus,
Rearranging the above equation into a linear form of the form y = mx + c,
Thus, a plot of vs should give a straight line with a slope of 55.5Cb/c and a Y-axis intercept of 55.5/cCbT1a, while the X-axis intercept provides a direct readout of the exchange rate . Once the exchange rate is determined, the concentration can be determined from the slope or the Y-intercept. The Omega plot has the advantage in that it provides a concentration independent approach to determine the water exchange rates of PARACEST agents. This may prove to be especially important in vivo.
The current focus of many imaging techniques involves the design of sensors to specifically measure certain physiological parameters such as pH and temperature or to signal the presence of metabolites and biomarkers in tissues. Many techniques including electrochemical,179–180 fluorescence,181–182 magnetic resonance,183 and PARACEST have been employed to selectively target nitric oxide,184 glucose,169,185 lactate,170 and metal ions186 among others.187 In applying these techniques to a biological target, there are a number of obstacles that must be resolved before they could be considered useful clinically. There are several properties a good sensor must possess: (a) the agent must be kinetically and thermodynamically stable in biological media, (b) it must be nontoxic, (c) should not disturb the biological system of interest at a molecular level, and (d) should be minimally invasive for patient comfort.148,188 Direct detection of the analyte is generally considered best, although this principle may not always hold true. Ideally, interactions of the analyte with the sensor should be reversible, particularly in cases where the analyte plays a key role in biological cascades needed for homeostasis, and the binding affinity should approximately match the concentration of the analyte typical of normal biology.
Several biological sensors based on PARACEST agents have been reported.4 The ability of a PARACEST agent to be selectively turned on and off by varying the frequency of the selective presaturation pulse is a substantial positive for this platform because this offers the possibility of detecting multiple biological responses in a single experiment. It is thus reasonable to suggest that further discoveries will result in useable sensors that respond to early biomarkers of disease, report on the effectiveness of disease treatment, and provide insights into a variety of metabolic processes.
pH plays a key role in biological homeostasis, with small changes serving as indicators of various abnormalities. For example, healthy tissues have an extracellular pH around 7.4 while cancerous tissues are more acidic, typically having an extracellular pH of 6.8–6.9.189 Therefore, it is not surprising that pH mapping was one of the earliest and most important targets for contrast agents. The earliest studies of diamagnetic CEST molecules already recognized that the magnitude of CEST response was pH dependent.1,52 Consequently, in theory, any molecule that shows a pH dependent CEST response could potentially be used as a pH imaging agent.
A number of Gd3+-based complexes display a change in relaxivity with changes in pH and can therefore be used as pH imaging agents.189–190 However, the contrast enhancement observed can either be a result of change in local concentration of the agent or a pH effect or both. One approach to overcome this problem is to use a dual contrast agent method involving the sequential injection of two different agents with identical pharmacokinetics. The pH sensitive agent is injected first followed by a second agent that mimics the tissue concentration distribution of the first agent, but is pH insensitive.6 While pH mapping studies have been done using this protocol in animals, it is impractical for clinical applications.. Unlike Gd3+-based agents, two different PARACEST agents can be administered and imaged simultaneously negating the need to know the exact tissue concentration of the sensor. As an alternative, single agent ratiometric methods can also be used, provided that the agent has two CEST sites whose pH response is different. In this case, two images can be collected and the ratio of those two CEST responses should directly report tissue pH. Two ratiometric methods have been proposed for monitoring pH in tissues based upon using two exchanging pools either on two separate molecules or within the same molecule.
Aime and co-workers were the first to propose the use of the ratiometric method involving two PARACEST agents, Yb3+-51 and Eu3+-51, to measure pH in vitro.167 They investigated a series of Ln3+-51 complexes and found that the Yb3+ complex showed the most promising CEST properties for the amide protons, based on their exchange rates and pH dependence, vide infra. While Yb3+-51 displays an amide-based CEST response at −16 ppm, no Yb3+-OH2 peak was observed, presumably because water exchange in this complex is too fast to be observed on the NMR timescale. It is assumed that the pH response of Yb3+-51 is base-catalyzed as demonstrated for a variety of other -NH protons (Section 188.8.131.52). A linear relationship between CEST and pH was observed over the pH range from 5.5 to 8.1 for Yb3+-51. This first generation pH sensor showed nice differences CEST contrast as shown in Figure 26. A maximum CEST contrast effect was seen at pH 8.1 but is almost completely diminished at pH ~ 6. The Eu3+-51 contrast agent, needed purely for ratiometric analysis in this case, shows a bound water CEST response that is virtually independent of changes in pH at ~50 ppm downfield from bulk water. Figure 27 illustrates the CEST spectrum resulting from a combination of 16 mM Eu3+-51 and 20 mM Yb3+-51 at pH 8.1. The water exchange resonance from Eu3+-51 is visible as a broad shoulder at ~ 50 ppm while the –NH exchange component from Yb3+-51 is less well-resolved on upfield shoulder of the bulk water resonance. Nevertheless, this study set the groundwork for other PARACEST pH sensors by establishing the ratiometric method. Additionally, in vivo applications of this dual sensor system are theoretically possible because it is relatively safe to assume identical bio-distributions for these two complexes because they have identical charges and hydrophilic/lipophilic balances.
An alternative method which would not require the use of two imaging agents in vivo is possible if the two CEST responsive pools were in the same molecule. The Yb3+ complex of 17 proved to have two observable CEST effects for the amide protons.155 These protons are magnetically non-equivalent and result in two broad resonances at −14.5 and −17.7 ppm. As shown in the inset in Figure 21, these protons are quite sensitive to changes in pH. Surprisingly, Zhang et al. observed that τM varies only slightly in the range of pH 6.46 – 7.84, indicating that these resonances should serve as efficient antennae for transferring magnetization to bulk water at pH value relevant to biological conditions.155 In fact, a 30 mM solution of Yb3+-17 showed 52 and 42 % decrease in bulk water signal when the –NH protons were irradiated separately and a 70% decrease upon simultaneous radiation. Such promising results led to the use of this complex in dendrimer studies (Section 4.4.9). Nevertheless, this is a perfect example of using one molecule with two exchangeable sites for the ratiometric detection of pH.
Another promising example of a PARACEST agent with two detectable exchanging sites is the Pr3+ complex of 51. This complex shows a measurable CEST response from the metal bound water molecule as well as the amide protons, making it a dual pool sensor.154 As shown in Figure 28, the bound water CEST effect is independent of pH while that of the amides is highly pH dependent in the range of 5.0 – 8.0. Here, water exchange is independent of pH while base-catalyzed exchange of the amide protons in the same molecule has their usual pH dependency. As the two CEST exchange peaks in this complex are well separated in terms of frequency, the use of this complex to monitor pH in vivo would be exciting. However, the CEST response of a 30 mM sample at pH 7.4 is only about 40% likely making this agent too insensitive for in vivo studies. Nevertheless, if concentration issues can be solved, a sensor of this type would be optimal for in vivo applications.
Understanding the relationship between variations in tissue temperatures and metabolism, physiology and disease states is a topic of great interest.168 This is, in part, due to recent studies that show elevated temperatures in cancer as well as in some heart conditions.191–193 In vivo tissue measurements include invasive thermocouples that are highly sensitive but are limited in that they measure only local tissue temperature and not simultaneous organ temperature gradients.194–195 A number of MRI techniques have been proposed as non-invasive methods for imaging temperature. The most widely used, proton resonance frequency (PRF) shifts, has shown promise in this arena with measurements that show an accuracy of ±1–3 °C. This technique is based on the observation of Hindman,195 who noted that the frequency of water proton resonance has a linear relationship with temperature over the range of 37 – 43 °C. This range of accuracy would be appropriate for tumor ablation but would not be widely applicable for imaging corticoid atherosclerotic plaques which only differ in temperature about 0.2–0.3 °C compared to normal tissues. Other techniques, such as hyperfine shift agents or paramagnetic relaxation agents, are beleaguered by poor spatial resolution and high concentrations need for imaging.196 PARACEST thermometry imaging, however, uses the bulk water, present at 40–50 M in vivo, for detection so this improves sensitivity considerably.
There are two factors that response to changes in temperature in these complexes, the rate of water/proton exchange and the magnitude of the lanthanide-induced hyperfine shift. Thus, both the magnitude and the shape of the CEST exchange peaks are extremely sensitive to temperature variations. As all PARACEST agents response to changes in temperature, the following examples are intended to illustrate the basic principles of using these systems in vivo. Terreno et al. took advantage of the two pools of exchangeable protons in Pr3+-51 and utilized a ratiometric approach to measure temperature gradients over a range of 295 – 330 K.154 An inverse linear relationship between temperature and CEST from water exchange was found in this complex (Figure 29). The amide protons, however, showed a more complicated behavior with changes in temperature. Here, the CEST response increased with temperature to a maximum at around 315K and then decreased with further increases in temperature. By taking the ratio of these two pools, -NH and Pr3+-OH2, a concentration-independent measure of temperature was described. Using a somewhat different approach, Zhang, et al. demonstrated that the frequency of the CEST water exchange peak in Eu3+-51 could be used as a temperature sensor.168 Given that the frequency of this exchange peak is independent of concentration (even though its magnitude is concentration dependent), one can extract an exact temperature measurement by simply scanning the frequency of the presaturation pulse over a narrow range of frequencies to locate the peak maximum. To demonstrate this, temperature maps were generated on phantom samples containing 10 mM Eu3+-51, pH 7 at different controlled temperatures. The images in Figure 30 show that even small temperature gradients are easily detected, an excellent property for biological thermometry. This rather simple first generation PARACEST thermometer showed a ~690 fold increase in temperature sensitivity over the more widely used PRF technique.
Li and co-workers set out to design a temperature sensitive agent capable of producing a maximum CEST response at physiological temperature. To achieve this, an additional amino acid phenylalanine was added to 51 to produce 57.194 This resulted in a slowing of water exchange and a larger maximal CEST effect at temperatures near physiological range (see Table 11). The temperature dependence of the chemical shift of Eu3+-57 was shown to be linear but not quite as sensitive to temperature changes compared to Eu3+-51 (Figure 31). This new complex provided greatest sensitivity near 38°C with an ~75% increase in CEST effect compared to Eu3+-51 under the same experimental conditions, a substantial advantage for biological applications. It is also obvious from the figure that above 40°C, rate of exchange between the bound water and the bulk water becomes fast and results in a decrease in intensity of the CEST signal. This trend can be explained by referring back to Figure 23(b) where a maximum CEST effect is seen when the rate of water exchange is equal to 2πB1 (Section 4.4.4), with any other data point showing less than optimum CEST. Though the bound water in Eu3+-51 is more sensitive to changes in temperature, greater B1 power is needed to realize the same CEST effect as that found with Eu3+-57 (assuming equivalent agent concentrations). Standard spin-echo imaging with phantoms showed that this complex was accurate to ±0.5 °C of the known value and the standard deviation of the pixel temperature was ±0.2 °C (Figure 32). These studies show extensive progress in the area of thermometric PARACEST agents and hold great promise for in vivo applications.
The central importance of glucose production, storage, transport, and utilization in mammalian tissues provides motivation to develop non-invasive methods for the detection of glucose in vivo.169,197 Additionally, how these pathways respond to treatment of disease states such as diabetes and cancer would invariably aid in optimizing modes of treatment. Several glucose detection techniques are currently available. Glucose sensors have been developed using ultraviolet, fluorescence, circular dichroism (CD), and electrochemical methods but are invasive and not feasible for in vivo applications.198–205 Noninvasive techniques such as 1H magnetic resonance spectroscopy (MRS) has also shown potential for biological applications but is complicated by multiple overlapping metabolite signals with similar chemical shifts.206 The use of PARACEST imaging agents would provide a non-invasive solution with potentially higher spatial resolution. Unlike the pH and temperature sensors described above, an imaging agent for sensing glucose was designed around the concept that water exchange would be altered upon binding of glucose to specific receptors on the agent surface. Trokowski, et al. first illustrated this design by ligand 58 which contains two phenylboronate groups for binding reversibly with sugar cis-diol groups.185,197 The absence of a Eu3+-OH2 peak in the 1H NMR spectrum of Eu3+-58 indicates that water exchange is too fast to be detected on the NMR time scale in the absence of glucose but, interestingly, the water exchange rate slows about 2 fold upon addition of glucose. This was attributed to glucose forming a capped complex with the two phenyboronate groups, thereby hindering Eu3+-OH2 exchange as depicted in Figure 33. This model was supported by the appearance of multiplets due to 2J-coupling between glucose and the boronate arms and by circular dichoism studies confirming that a single molecule of sugar binds to the two arms of Eu3+-58.169 While phenylboronates have been shown to bind preferentially to fructose, competitive sugar binding studies showed preferential binding with Eu3+-58 falls in the order, glucose > fructose > galactose. These studies indicate that the macrocyclic ligand helps to position the phenylboronate arms in such a manner that is advantageous for glucose binding. To further support this model, Trokowski et al. showed that a ligand with two phenyboronates cis to one in the macrocycle exhibits preferential sugar binding in the order of fructose > galactose > glucose. This technology was successfully tested recently in ex vivo perfused mouse livers.197,207 In this experiment, hyperglycemia was modeled by perfusing the isolated liver from a fed mouse with 10 mM glucose and 10 mM sensor. A control liver was isolated from a 24 hr fasted mouse and perfused with the same amount of sensor, but no glucose. As shown in Figure 34(a), the CEST spectroscopy of the liver perfused with glucose and agent showed a distinct CEST peak near 42 ppm, while the second liver perfused only with the agent showed no clear water exchange peak. The appearance of a water exchange peak in the CEST spectrum after glucose addition demonstrates qualitatively that water exchange is slower when glucose is bound to the complex. Quantitatively, the change in CEST intensity induced by glucose was a 17 % decrease in water signal. Further analysis with the proposed sensor demonstrated that at 10 mM concentration of the sensor, changes in glucose concentration over 0–20 mM can be detected (binding constant for the Eu3+-60-glucose system is 2275 ± 266 M−1 at pH 10.2 and 339 ± 29 M−1 at pH 7). This is a reasonable range for biological detection where a typical blood glucose level in the diabetic state is 7–10 mM. This work helped to establish the principle that designing complexes based on chemical principles of reactivity is an elegant way to design “smart” imaging agents.
Redox active species including metal ions and oxygen species are tightly regulated in vivo under most circumstances. However, the amount of energy available in biological redox systems can have incredibly detrimental effects if left unregulated.207–208 More specifically, disruption of redox homeostasis is suspected to play a key role in Alzheimer’s disease, Parkinson’s disease, cancer, stroke, and atherosclerosis.209 Understanding the relationship between redox active species and their function and/or malfunction would undoubtedly advance treatment of many diseases. Therefore, the development of non-invasive imaging techniques to monitor tissue redox has been the subject of many investigations in recent years. Optical techniques are quite sensitive for redox active species such as NADH at sub-cellular levels but are plagued by attenuation with increasing tissue depth limiting their use in vivo.210–213 Paramagnetic redox-sensitive nitroxides have been used as contrast agents in MRI but are limited by low relaxivity and relatively short lifetimes in vivo.214–215 Ratnakar et al. recently demonstrated that the p-nitro moiety of Eu3+-38 may be reduced to the p-amino derivative and that the resulting complexes have remarkable different CEST properties (Figure 35).145 While this reaction is not biocompatible, this experiment demonstrates the potential of developing redox-sensitive PARACEST agents based on this platform for use in vivo.
L-Lactate, the end-product of anaerobic glycolysis, is known to be elevated in stroke, cancer, cysts, during brain activation, and a few other metabolic disorders.216 Elegant methods have been published for detecting lactate by 1H spectroscopy and imaging but overlap of the lactate proton resonances with those of fats and other tissue metabolites makes quantitative lactate imaging quite challenging. Lactate is easier to uniquely detect by 13C NMR but even with enrichment, the sensitivity of 13C for detection of lactate in vivo is difficult except in the most favorable circumstances. 217 Recent advances in hyperpolarized 13C for detection of lactate in tumors has shown considerable promise.218 The addition of a lactate-specific imaging agent as an alternative method could be beneficial in some cases. Aime et al.170 recently reported a PARACEST detection system for lactate. The amide protons of the lactate-free complex, Yb3+-59 that appears as a broad resonance in the 1H NMR spectra at −28.5 ppm is transformed into a narrow set of doublets at −14 ppm and −20 ppm upon binding of L-lactate to the Yb3+ coordination sphere (lactate acts as a bidentate ligand in this case). The separation of ~10 ppm between Yb3+-59 and the ternary adduct Yb3+-59-lactate indicates that this agent can detect both the bound and free forms of the complex. Figure 36 shows the CEST spectra of the free and bound forms of Yb3+-59 at pH 7. However, unlike the significant difference observed in the high resolution spectra between the free and bound forms, the CEST spectra show partial overlap in the CEST peaks. Nevertheless, pre-irradiation at −29.1 ppm showed significant changes in the CEST effect as a function of lactate concentration (from 60 % to 0 %) while irradiation at −15.5 ppm shows only modest changes. As later observed by Woods et al, this system of irradiation at two different frequencies could also be used to give a concentration independent ratiometric method for the determination of lactate concentration in solution.4
Nitric oxide (NO) plays a key role as a signal transmitter not only in the vascular endothelium, but also in the central and peripheral neurons. A deficiency in NO has been shown to occur in atherosclerosis, subarachnoid hemorrhage, and ischemia-reperfusion injury. Detrimental affects also appear if NO levels are too high. Excess NO is believed to induce epileptic seizures and even cause neuronal injury after a stroke.219 The NO radical is produced in μM quantities by the enzyme nitric oxide synthase (NOS) in the presence of oxygen.220 In addition to its low concentration, the extremely short lifetime of this biological messenger molecule makes it very difficult to detect in vivo, especially by MRI. One way to detect NO using MRI was successfully demonstrated by Fujii et al. by using spin-trapping reagents, in a technique they called “MRI spin-trapping”.221 NO has been shown to form reasonably stable complexes with various Fe-based compounds. In fact, it is known that hemoglobin and deoxyhemoglobin act as natural spin-traps for NO in the body.222 Using this idea, Fujii et al. used an Fe(II) chelate N-methyl-D-glucamine to form a spin-trap complex with NO, (MGD)2-Fe(II)-NO, and successfully visualized and mapped the site of NO generation in rats.
Another approach to detect the short-lived NO radical would be to use responsive contrast agents that undergo irreversible reactions with NO leading to products that can be detected easily using MRI, as shown recently by Liu et al. using a PARACEST agent.184 They report the use of the Yb3+-DO3A-orthoaminoanilide derivative, Yb3+-60, for NO detection. The orthoaminoanilide derivative, with two different types of exchangeable protons, shows CEST exchange peaks at −11 ppm and +8 ppm due to the amide and amine functional groups, respectively. Upon reaction with nitric oxide in the presence of oxygen, the compound undergoes irreversible conversion to a triazene derivative, leading to permanent loss of both the CEST peaks (Scheme 1).
Figure 37 shows the phantom images obtained before and after treatment of the complexes with NO and O2. In an approach similar to Aime’s work on pH sensitive agents,167 Pagel also used a Tm3+-51 complex that is unresponsive to changes in NO levels to track the pharmacokinetics and biodistribution of these agents. It can be clearly observed from the phantom images that significant loss of signal is observed for both the amine and the amide peaks of the Yb3+ complex after treatment with NO while the Tm3+ complex remained unaffected. While loss of the amine signal at +8 ppm can be easily explained by the conversion of the aromatic amine to a triazene molecule, the loss of amide peak at −11 ppm is attributed to conformational changes that result in the amide peak shifting further away from the lanthanide ion. Thus, irreversible smart agents such as this have the added benefit that they can detect short-lived species like the NO radical by detecting an irreversible product that accumulates with time.
Divalent zinc, ranking as the 2nd most abundant transition metal ion in the body (after iron), plays ubiquitous biological roles.223 The average human is thought to have about 2.3 grams of total zinc.208 Zn2+ participates in apoptosis, regulation of synaptic transmission and cell death.189 An understanding of its tissue biodistribution would facilitate a better understanding the biological regulation of this ion. Although a number of Gd3+-based Zn2+ sensors have been reported, they suffer from the same problem as most other small molecule Gd3+ agents in that the change in T1 relaxivity is typically quite small. To test the potential advantages of PARACEST for Zn2+ detection, Trokowski et al. prepared Eu3+-61 as a sensor for Zn2+. The agent, like earlier Gd3+-based systems, contained two dipyridyl amine moieties that show strong affinity for zinc ions.186 It was anticipated that binding of Zn2+ to the two dipyridyl amine (dpa) arms above the Eu3+-inner sphere water coordination site would result in a significant change in the rate of water exchange and thereby CEST. The CEST spectrum of a 20 mM buffered solution of Eu3+-61 shows a characteristic CEST effect due to bound water at 50 ppm in the absence of zinc. Upon addition of equimolar quantities of zinc ions to the sample at pH 7.1 resulted in significant broadening of the CEST exchange peak at 50 ppm (Figure 38a). At pH 8, this same CEST peak was so broad that it nearly disappears from the spectrum, consistent with an increase in water exchange (Figure 38b). One possible reason for dramatic broadening of the bound water peak could be attributed to changes in the coordination geometry of the Eu 3+ complex upon Zn2+ binding. It is well known that only the SAP isomer of DOTA-tetraamide lanthanide complexes exhibit slow enough water exchange for CEST so it is possible that binding of Zn2+ at its binding site induces a change in structure in the Eu3+ coordination sphere from the SAP to TSAP isomer, thereby increasing water exchange. However, high-resolution NMR spectra indicate that this was not the case and the Eu3+-retains the SAP geometry upon binding of Zn2+. Further studies indicated that, while most ligand systems containing the dpa arm exhibit a 2:1 coordination with Zn2+, Eu3+-61 forms only a mononuclear complex. In fact, potentiometric studies indicate that above pH 6, a Zn2+-OH species is formed and it is this species that is responsible the increase in water exchange that occurs upon binding Zn2+ binding through base-catalyzed exchange of protons between t he Eu3+-bound water and bulk solvent.
Most PARACEST agents discussed so far require analyte concentration on the order of 1–10 mM for detection, in most cases limited by the binding affinity between the analyte and the sensor. This places serious limits on potential applications in vivo for many molecular targets that exist in lower concentrations. A different approach in sensor design was developed by Yoo et al. using PARACEST active complexes.224–225 They set out to exploit the use of enzymatic catalysis to alter the structure of the sensor molecule, thereby changing its spectral properties. Their general strategy involved linking an enzyme substrate to a lanthanide complex. In the presence of the enzyme, the substrate would be consumed while the sensor is converted into a new structure with different CEST properties. This approach allows indirect detection of relatively low concentrations of enzyme activity. Given that proton exchange between a primary amine and water is generally faster than between an amide and water, the basic concept was to construct the DOTA-peptide conjugate that could be cleaved by an enzyme leaving behind a primary amino group. The system chosen by Yoo et al. (shown in Scheme 2) consists of the TmDOTA complex with one arm of the four attached to a tetrapeptide, Asp-Glu-Val-Asp (DEVD),225 that is cleaved by caspase-3, an “executioner” of the metabolic death cascade during cell apoptosis. Imaging the levels of this enzyme could be important in evaluating patient response to apoptosis-promoting antitumor therapies. The intact non-cleaved Tm3+ complex shows a CEST exchange peak near −51 ppm corresponding to the amide proximal to the lanthanide ion. Incubation of 48 nM caspase-3 with the sensor (25 mM) for 1 h (37°, pH 7.4) showed a significant decrease of the PARACEST effect at −51 ppm and a new response at +8 ppm. The new CEST signal was consistent with the conversion of the TmDOTA-DEVD amide 62 to the TmDOTA-amine derivative shown in Scheme 2. As shown in Figure 39, a 14.5 % decrease in water signal intensity in observed in the parametric map obtained after irradiation at −51 ppm in the absence of the enzyme, while no significant change in the water signal intensity was observed after the enzyme reaction. In this preliminary design, the targeted substrate ligand arm involved conversion of an amide to an amine, which would accelerate the rate of chemical exchange from ~300 to ~3000 s−1 as well as change the frequency of the exchange. Unfortunately, the reduction in CEST signal upon reaction with caspase-3 is opposite what one would prefer in a biological sensor so this aspect will likely limit the application of this sensor in vivo. However, the principle behind this PARACEST design is quite unique. In a later report, Yoo et al. extended this work to a two sensor system. The second sensor, Yb3+-62 shows an amide based CEST response at −16 ppm.224 Upon addition of enzyme, (1) the CEST response at −51 ppm dissipates, (2) a new CEST response at +8 ppm is observed and (3) the magnitude of the CEST response at −16 ppm remains unchanged. Monitoring the CEST effect of two complexes simultaneously would allow a better visualization of enzymatic activity. In fact, preliminary PARACEST MRI experiments on phantom samples showed that the activity of caspase-3 and caspase-8 could be monitored using this two sensor system.
A similar method was employed by Chauvin et al. with Yb3+-63.226 The approach couples an enzyme-specific substrate, in this case αBz-βGal, to a lanthanide chelating moiety via a self-immolative spacer. This Yb3+ derivative shows no CEST response initially, but upon introduction of the β-galactosidase enzyme, the spacer is cleaved and eliminated via the enzymatic action producing amide based CEST effects at −16.7 and −20.5 ppm. The enzymatic cleavage reaction is relatively fast, t1/2 = 68 min. As shown in Scheme 3, the subsequent electron cascade of the enzyme converts the amide to an amine, with new spectral properties. CEST images obtained from 20 mM samples of Yb3+-63 before and after enzymatic reaction are shown in Figure 40. Analysis of these images indicated that no significant change in water signal intensity is observed for Yb3+-63 but a 29% CEST effect is observed after reaction with enzyme. One main advantage of this system is that the presence of enzyme activity induces a new signal so that one is not simply monitoring loss of signal of the starting material. Although it still remains to be seen whether these agents can be used at the concentrations that are more compatible for in vivo applications, this methodology offers tremendous advantages in terms of the fact that a whole library of enzyme specific sensors could quickly be built and deployed giving a plethora of data from a single experiment.
A popular approach to blood pool Gd3+-based T1 agents for vascular imaging traditionally involves reversible non-covalent conjugation of the contrast agent with high molecular weight macromolecules like albumin.227–229 Serum albumin is one of the most popular macromolecular binding targets in vivo owing to its abundance in blood and its reversible binding characteristics to a wide variety of biomolecules. Non-covalent conjugation of a Gd3+ complex to albumin increases the “effective” concentration of the agent and results in extended “half-life” of the compound in the blood, enabling imaging of multiple vascular regions after a single injection of agent.230 One could envision an analogous improvement in CEST sensitivity and significant changes in CEST by reversible binding of a PARACEST agent with albumin. To test this hypothesis, Ali et al. used the o-benzyl functionality well-known to impart high binding affinity for albumin,228,231 as the design basis for an albumin binding PARACEST agent.232 They designed two DOTA-based tetraamide derivatives (64 and 65), either di- or tetra-substituted with the o-benzyl functionality.
The Gd3+ complexes of both the ligands show considerable affinity for site II of HSA (Human Serum Albumin) with binding constants of 3.14 × 104 M−1 for Gd3+-64 and 1.04 × 104 M−1 for Gd3+-65 and very poor affinity for site 1.232 They exhibited a modest ~two-fold increase in relaxivity upon binding to HSA (smaller than most other HSA binding T1 agents), consistent with relatively slow water exchange in typical GdDOTA-tetraamide complexes. Working with the more soluble di-substituted Eu3+-64 complex, the authors found no change in the CEST properties of the complex upon binding with HSA (Figure 41). τm values determined from fitting the CEST spectra to Bloch equations indicated a two-fold increase in water exchange kinetics upon albumin binding. While there is still a need to address the sensitivity of these albumin binding PARACEST agents, the fact that there is very little effect on the CEST properties of these Eu3+ complexes upon albumin binding makes them potentially useful as vascular imaging agents.230
Peptide-based targeting vectors are increasingly becoming popular for therapeutic and diagnostic applications and a large number of DOTA-peptide conjugates have been synthesized and successfully used for targeted MRI and nuclear medicine applications.233–234 In addition to the development of molecular and enzyme specific responsive agents, an entirely new class of PARACEST imaging agents is beginning to emerge. These new PARACEST agents composed of oligopeptide-based systems offer an attractive alternative to the design of biocompatible contrast agents. An interesting array of oligopeptide-derivatized DOTA-tetraamide ligands (57, 66–75) was recently reported by Bartha et al. to identify new PARACEST complexes of different Ln3+ ions as targeting vectors.235–237 While Eu3+-based PARACEST agents have been most widely studied, Wojciechowski et al. ventured into other Ln3+ ions with an intention to discovering new Ln3+-based PARACEST agents. They reported that the different Ln3+ complexes of the oligopeptide derivatives show different CEST effects owing to their different water exchange kinetics.235
In order to understand the effects of varying the amino acid side chains on the CEST properties of these oligopeptide derivatized ligands, Suchy et al. conceptualized a new series of ligands by changing one amino acid from their original series of dipeptide-based ligands. By using ligand 57 as the basis for comparison, owing to the large CEST effect of its Eu3+-bound water protons and favorable temperature dependent CEST properties,194 they synthesized ligands 70 and 71 by replacing amino acid phenylalanine with tyrosine and tryptophan respectively. Interestingly, while the tyrosine-based Eu3+ derivative retained a similar CEST effect as Eu3+-57 (Figure 42a), the corresponding tryptophan derivative did not (Figure 42b).236
In another report, Suchy et al. synthesized a new ligand to alter the in vivo biodistribution of their targeted oligopeptide derivatives.237 They conjugated ligand 57 to a cell-penetrating peptide, a cystamine derivative, with the expectation that the disulfide bond when cleaved in a reductive intracellular environment may result in trapping the agent within the cell. The CEST properties of the neutral complex Eu3+-75 compares favorably with its anionic precursor, Eu3+-57, thus providing a new methodology for advanced PARACEST agents with favorable in vivo properties.
Like other MR imaging agents, an intrinsic limitation of small molecule PARACEST agents is sensitivity. The lowest detection limit reported for a low molecular weight PARACEST agent gave a 5% decrease in water signal intensity using a concentration of 0.5 mM agent.238 The usual injected dose of clinical Gd3+-agents is 0.1 mM/kg, but for MR angiography, doses up to 0.4 mM/kg are used. This translates to peak blood Gd3+ concentrations in the low mM range.239–240 Preliminary data suggest that the concentration requirements for PARACEST agents for in vivo experiments are similar to those of the Gd3+-based T1 contrast agents currently in use. In vitro experiments of several Eu3+-DOTA tetraamide PARACEST agents in plasma were shown to have easily detectable CEST effect in this concentration range (for example 8.9 % decrease in bulk water intensity was measured for complex Eu3+-47 in 1 mM solution at 20°C, with 5 s irradiation time and B1 = 14.1 μT).138 Theoretical considerations also suggest that the detection limit of a single PARACEST exchanging species with optimal water exchange rate, chemical shift and relaxation properties should be comparable to a single Gd3+-based T1 agent.67,241 Nonetheless, molecular imaging often requires detection of biomarkers present at extremely low concentrations; so there is clearly a need for the development of PARACEST contrast agents with enhanced sensitivity and improved targeting abilities.
One approach to increase the sensitivity of CEST agents, demonstrated by Goffeney et al. for DIACEST agents, involves increasing the number of exchangeable proton sites per CEST molecule (Section 4.2.2).78 Using cationic polymers (50–800 kD MW), with a large number of amide and amine protons sites, they established that a CEST effect of almost 50% could be attainable with micromolar concentrations of the agent. In an attempt to improve on the small Δω values of diamagnetic CEST macromolecules, Aime et al. exploited the formation of supramolecular adducts between the cationic polymers used by van Zijl et al. and a shift reagent.238 Lanthanides have a long history as frequency shift reagents for a variety of NMR applications.242–244 Hence, including a shift reagent capable of forming an adduct with a macromolecular target having multiple exchange sites could potentially increase Δω substantially without eliminating the advantages offered by polymerization (higher number of exchanging sites). It is also expected that exchange kinetics of the shift reagent need not meet the slow exchange condition for CEST. Aime et al. reported the use of Tm3+-12, a widely-used shift reagent for aqueous cations, and poly-L-arginine (PLR) as the macromolecular target.238 PLR has two potentially different CEST-active sites, the amide protons of the backbone and the guanidine protons of the side-chains. In the absence of Tm3+-12, PLR has a CEST effect from the amide backbone protons (depending upon pH) but none from the guanidine protons. This is a classic example of a situation where kex > Δω for the guanidine protons. Addition of Tm3+-12 to PLR results in the formation of a tightly associated ion-pair and the appearance of a CEST peak at 30 ppm (Figure 43). In addition to the dipolar interaction of the complex with the polymer leading to an increase in Δω, it is probable that the exchange rate of the guanidyl protons with the solvent is also affected by the tight association with the complex. Regardless of mechanism, the CEST requirement, Δω > kex is then met and CEST is observed. Using this binary system, it was found that only micromolar concentrations of polymer (1.7 μM) and shift reagent (30 μM) were required to produce a 5% decrease in water signal intensity.
Greater improvements in sensitivity can also be made by maximizing the number of PARACEST exchanging species at the targeted site. Low molecular weight PARACEST agents typically contain between one and four exchange sites per molecule, thereby requiring a high concentration of the agent to produce significant contrast. By increasing the number of metal centers, significant enhancement in contrast is expected. Multiple approaches to resolve this strategy have been identified, involving conjugation of the PARACEST agents to high-molecular weight scaffolds like polymers or dendrimers,241,245 or encapsulation of a large number of metal centers in the inner core of high-molecular weight systems like liposomes, 246 or apoferritin.247
Aime and co-workers pioneered the effort in this field with their work on liposome-based PARACEST agents, which they termed as LIPOCEST agents.9,246 Liposomes (Figure 44), discovered by Bangham et al. in 1961,248 are tiny nano-sized bilayer vesicles with an aqueous core surrounded by phospholipid bilayers. They have attracted widespread attention mainly as drug carriers, particles for targeting, and as contrast agents for MRI.9,249 One of the advantages of using liposomes as carriers for contrast agents comes from the fact that thousands of imaging reporter molecules can be encapsulated within the aqueous core of a single liposome, significantly increasing the local concentration of these agents. These biocompatible nano-vesicles have a large number of water molecules in their inner core that are in relatively slow exchange with bulk water. Aime and co-workers first reported the use of a new class of contrast agent composed of liposomes loaded with a paramagnetic lanthanide shift reagent (SR).246 Upon entrapment of a SR within the liposome, two signals are observed in the 1H NMR spectrum, one from the exterior bulk water and another from the water within the liposome core. The resonance of the water protons inside the liposomes is shifted from that of the external bulk water to an extent that depends mainly on the shift properties of the trapped reagent. In this case, the rate of water exchange between a coordinated water molecule on the SR and the entrapped water molecules must be in the intermediate-to-fast exchange regime so that all water molecules inside the liposome experience the paramagnetic effects of the SR as quickly as possible. Therefore, the CEST requirement (Δω kex) in this case must be met by the exchange rate of water or protons across the lipid bilayer.
Liposomes containing 0.1 M Tm3+-76 results in a shift of the entrapped water resonance to 3.2 ppm (Δintralipo), which is well separated from the bulk resonance (Figure 45).246 Phantom images of the LIPOCEST particles at different concentrations showed that 90 pM (liposome concentration) could be detected by CEST imaging, a tremendous improvement over traditional small molecule PARACEST agents. Thus, it was proposed that by using different types of paramagnetic shift agents, one should be able to prepare LIPOCEST agents with different CEST activation frequencies, giving rise to highly sensitive, “multi-frequency” CEST agents. 166,246 Nevertheless, the relatively small chemical shift of the intraliposomal water resonance (similar in magnitude to typical DIACEST compounds) remains a significant disadvantage for the in vivo use of these first generation LIPOCEST agents. The chemical shift of the intraliposomal water resonance results from the sum of two contributions (Eq 54).
For fast exchanging systems, the dipolar contribution ( ) is an average of the chemical shifts of the SR-bound water and bulk water; so it is limited by the amount of SR that can be encapsulated within the liposomes. Under optimal conditions, the magnitude of contribution amounts to no more than ± 4 ppm. The bulk magnetic susceptibility (BMS) contribution that arises from the alignment of the magnetic moment of the compartmentalized SR to the external static field, negligible for a spherical liposome, can be significantly increased by changing the alignment of the magnetic moment of the SR center with the external magnetic field8 by entrapping it in a non-spherical compartment.172 Non-spherical liposomes are most easily prepared by shrinking spherical liposomes under hyper-osmotic conditions.250–251 In addition to the above mentioned two factors, Δintralipo is also governed by both the magnetic moment (μeff) of the Ln3+ and orientation of the SR within the liposome, which in turn can be modulated by incorporating amphiphilic paramagnetic complexes within the liposome membrane. A significant Δintralipo shift of the bulk water resonance was observed at ~18 ppm for a non-spherical liposome that included Tm3+-77 as part of the lipid preparation (becomes part of the membrane) with Tm3+-78 in the inner core.172 An alternate non-spherical liposome with Dy3+-77 incorporated into the membrane and Dy3+-78 encapsulated in the inner core resulted in a Δintralipo shift of −45 ppm (Figure 46). This result provides a good insight into the dependence of the magnitude and sign of Δintralipo on the magnetic moment of different lanthanide ions.
Deli Castelli et al. also demonstrated that the sign and magnitude of Δintralipo are not only controlled by the magnetic moment of the paramagnetic metal ion, but it also depends on the orientation of the liposomes with respect to the B1. This can be modulated by using suitable ligands where the conformation of the ligand pendant arms and their mode of coordination change the orientation of the vesicles.252 When the non-spherical liposome is parallel to the axis of the magnetic field, the Δintralipo shifts downfield and vice-versa as depicted in Figure 47. It must be noted that while the three lanthanide complexes, Tm3+-78, Dy3+-78 and Gd3+-78, show positive, negative and null shifts respectively when encapsulated in spherical liposomes, non-spherical liposome preparations with the same complexes show random trends as a result of the orientation of the liposomes with respect to B1. A series of amphiphiles based on the Tm3+, Dy3+ and Gd3+ complexes with ligands 79–83 incorporated in the liposome membrane with the corresponding Ln3+-78 entrapped in the inner core, reveals that the orientation of the non-spherical liposome depends on the type of coordination cage provided by the ligand (open-chain or macrocyclic) and the identity of the hydrophobic side-chains (Table 4).
A recent comparison of mono-, di- and tri-nuclear Tm3+ complexes revealed that the induced shift of the inner core water protons is due not only to a BMS contribution plus the magnetic properties of the encapsulated SR, but it also depends on the concentration of the intraliposomal shift reagents. 171 The degree of the observed shift resulting from the binuclear complex Tm3+-84 is twice that of its mononuclear counterpart Tm3+-78 since the effective concentration of the metal in the binuclear complex is more that its mononuclear counterpart. However, the tri-nuclear Tm3+-85 complex did not follow this same trend because of its different coordination environment. This effect is due to the fact that the trinuclear complex, with the presence of the linear chelating groups between the macrocyclic cages, has a more flexible coordinating geometry than the macrocyclic DOTA type rigid coordination environment. Thus, the coordinating environment of the macrocyclic cage aligns the paramagnetic contact shift field of the lanthanide in an optimal manner to shift the water resonance more effectively than the linear coordination environment (Figure 48).
When the osmotic pressure of a suspension of liposomes made of Tm3+-80 and encapsulating Tm3+-78 is increased, the liposome loses its spherical nature and shifts Δintralipo away from the bulk water until an isotonic point is reached (Figure 49).253 Further, varying the concentration of the incorporated Tm3+-80 complex, while keeping the concentration of Tm 3+-78 constant results in a linear relationship with Δintralipo and saturation transfer ratio. On the other hand, when the concentration of encapsulated Tm3+-78 was changed with respect to the Tm3+-80 there was an opposite trend in both the shift of the intraliposomal water resonance and saturation transfer ratio. The concentration of Tm3+-78 was found to regulate the osmotic pressure of the solution which, in turn, governs the shape of the liposome. As the concentration of Tm3+-78 is increased, the osmotic stress felt by the liposomes gradually decreases and they tend to become more spherical.253
In another paper, Terreno et al. reported the first ex vivo detection of two LIPOCEST agents differing in their Δintralipo when co-injected in the bovine muscle.254 The spherical liposome with entrapped Tm3+-78 complex shows Δintralipo at +3 ppm, while the non-spherical Tm3+-80 agent shows a corresponding peak at +18 ppm. After injection of these liposomal preparations into bovine muscle, each LIPOCEST agent could be imaged independent of the other agent by applying the appropriate frequency selective presaturation pulse identified for agent. An analytical model for the size dependence of the liposomes on the magnetization transfer rate (MTR) for spherical liposomes containing Tm3+-2 was reported by Zhao et al.255 Three sizes of spherical liposomes containing the same concentration of Tm3+-2 showed a difference in the MTRasym (Figure 50) with the Δintralipo increasing as the liposome size decreases.
Langereis et al. reported the first temperature sensitive non-spherical LIPOCEST agent and 19F contrast agent as a dual contrast agent and potential carrier system for MRI guided drug delivery.256 At room temperature, the dual contrast agent displays the expected CEST effect due to the entrapped shift reagent Tm3+-78. Interestingly, two different CEST peaks are observed in the CEST spectrum of the agent at +11 ppm and −17 ppm, presumably due to two different orientations of the non-spherical liposomes with respect to the magnetic field. Upon increasing the temperature to 315 K, the liposomal melting temperature, the shift reagent was released from the inner core and the CEST signal consequently disappeared (Figure 51). As anticipated, the encapsulated 19F reporter molecule (NH4PF6) showed the opposite behavior. At room temperature with the paramagnetic shift reagent and the 19F reporter molecule entrapped in the inner core of the liposome, the 19F NMR signal was too broad to be detected. Upon warming the liposome to its melting temperature near 315 K, a sharp increase in 19F intensity was observed with concomitant loss of the CEST signal. Phantom MRI images of the temperature variation of the CEST effect and 19F fluorescence intensity are shown in Figure 52. This example illustrates the basic principles behind image-guided drug delivery.
Apo-ferritin, a naturally occurring nano-sized carrier, offers similar advantages as liposomes in the transport of a large number of paramagnetic metal centers trapped in a compartment along with limited number of water molecules.9,247 An earlier report demonstrated that a Gd3+-based T1 agent entrapped within the inner core of apo-ferritin has a substantially higher water relaxivity (80 mM−1s−1) compared to the agent in water.257 It was hypothesized that such a dramatic increase in relaxivity must reflect catalytic prototropic exchange of water protons by residues on the inner surface of apo-ferritin. This catalytic effect on water exchange could prove extremely beneficial to PARACEST agents whose CEST efficiency is highly dependent on the rate of water exchange. To test this hypothesis, Vasalatiy, et al. encapsulated Eu3+-16, a slow water-exchange system, within the negatively charged inner core of apo-ferritin and found that the CEST signal from Eu3+-16 completely disappeared upon encapsulation (Figure 53).247 Further studies showed that the water and/or proton exchange rate was significantly enhanced when the complex is trapped within the protein. Indeed these studies provide sufficient proof that (a) the water exchange process across the ferritin pores is extremely fast and (b) the exchange process between the Eu3+-bound water molecule and the water within the core is catalyzed by the residues on the inner surface of the protein core. While this catalytic increase in water exchange prove beneficial to the Gd3+-based T1 agents,9 the need for slow exchange kinetics make these systems unlikely carriers for PARACEST agents.
As an alternative to nano-sized PARACEST agents involving liposomes and apo-ferritin, Winter and co-workers successfully demonstrated the use of paramagnetic CEST agents conjugated to perfluorocarbon encapsulated nanoparticles for the detection of fibrin. 258 They used a methoxy-benzyl substituted DOTA-tetraamide derivative with a lipophilic tail for incorporation on a perfluorocarbon nanoparticle, Eu3+-86. The PARACEST nanoparticle emulsion, with a diameter of 294 nm and a Eu3+ concentration of 2.5 mM, shows significant CEST at +52 ppm, while control nanoparticles without the attached PARACEST agent showed no CEST. Two-chambered phantom images of the PARACEST vs control nanoparticles provided further proof of the saturation transfer properties of the PARACEST agent on the exterior surface of the nanoparticles, with over 10% change in signal intensity observed for the bulk water. Winter et al. also demonstrated that these lipid-encapsulated perfluorocarbon nanoparticles, similar to their earlier report with Gd3+-based contrast agents,259 are capable of targeting fibrin. Clot images treated with both control and PARACEST nanoparticles reveal that while the clots themselves display a uniform hypo-intensity with respect to the surrounding buffer, the surface of the clots could be clearly distinguished in the PARACEST treated images, thereby providing the first example of a targeted PARACEST agent for detection of fibrin (Figure 54).
While the use of physical aggregates like liposomes and other naturally occurring nano-sized carriers is very attractive, shear force and other environmental effects such as high-dilution, temperature and pressure make them vulnerable in vivo.260 An alternative approach to enhance CEST sensitivity would be to conjugate PARACEST agents by covalent modification to scaffolds such as dendrimers or polymers. Numerous studies of Gd3+-based T1 agents conjugated to dendrimers have revealed that they are excellent molecular substrates to attach multiple metal centers with their symmetric, nano-sized, well-defined and highly branched architectures.229 Taking advantage of this fact, Pikkemaat et al. designed a Yb3+ complex with exchangeable amide protons (CEST effect observed at −15 to −20 ppm) conjugated to a symmetrical poly(propylene imine) dendrimer to serve as PARACEST agents.245 In addition, the local pH of these agents can be modulated by the tertiary amines of the dendritic core of the poly(propylene imine) dendrimers, making them attractive as pH-responsive contrast agents. As expected, the G3 (3rd generation) dendrimer with 16 end-terminated Yb 3+-17-modified groups (structure 89) is 16 times more sensitive than its mono-nuclear counterpart Yb 3+-87, while the G1 dendrimer (1st generation) is 4 times more sensitive (structure 88). All three compounds are pH sensitive with the maximum changes in CEST observed over the physiological pH range, 6.5 to 7.5 (Figure 55a). The mono-nuclear compound Yb3+-87 shows a maximum CEST effect at a pH of 7.5, while the G3 dendrimer (structure 89) shows a maximum at pH 6.5; this difference can be attributed to the fact that the presence of multiple tertiary amino groups create a more basic micro-environment which in turn could result in larger exchange rates and also shifts the pH profile (Figure 55b). While these poly-nuclear dendrimer-based PARACEST agents offer only limited increase in sensitivity over the supramolecular adduct designed by Aime (millimolar for small dendrimers versus low micromolar range for LIPOCEST), fine-tuning of the dendritic core and increasing the size of the framework remains an attractive possibility.
Like dendrimers, polymer-based imaging and therapeutic probes have seen significant advances in the last decade for the diagnosis and treatment of various diseases.260–262 The combination of well-defined chemical compositions, increasing sophistication in polymer chemistry and ease of chemical modification has resulted in an increased interest in polymeric imaging probes. In an effort to use this technology to improve the efficiency and detection limit of PARACEST agents, Wu et al. designed Eu3+-based PARACEST agents with one polymerizable side-chain moiety, 90 and 91.241 Free radical polymerization of the monomer with different percentages of the initiator yielded polymers, 92 and 93, of differing molecular weight. It was shown that while the CEST efficiency on a per Eu3+ basis remains essentially the same in the polymer as that of the corresponding monomer, the molecular sensitivity was of course significantly enhanced. With a 5% change in water signal intensity set as a lower limit for CEST detection, it was found that the longer polymers exhibit detection limits in the low micromolar range, approaching the levels required for targeted molecular imaging applications (Figure 56).
While all the examples discussed in this section clearly demonstrate the significant improvements to sensitivity of PARACEST agents, further optimization of the chemistry and other parameters is needed to fully realize the promise of CEST-based agents for in vivo studies.
Given that the vast majority of all PARACEST agents reported to date have been derived from DOTA-tetraamide ligands, a brief review of the synthetic methods to prepare these ligands is in order. Tetraamide derivatives of DOTA are traditionally obtained by the alkylation of cyclen (1,4,7,10-tetraazacyclododecane) with readily available bromo- or chloroacetyl amides in a dipolar aprotic solvent such as acetonitrile or DMF in the presence of a base such as K2CO3 or iPr2NEt (Scheme 4).42,138,140,145,263–265 This method is quite general, and gives reasonable yields for both primary and secondary amines with a wide variety of functional groups. In addition, it can be extended to the synthesis of DOTA amides with mixed sidearms using selective protection-alkyation-deprotection-alkylation steps starting from cyclen.266–267 This route was also adapted for the synthesis of peptide-based DOTA tetraamides. Utilizing this strategy, a large library of ligands and their subsequent lanthanide complexes were produced in excellent yields.235–237
Another approach to amide derivatives of DOTA involves the activation and conversion of the carboxylate groups of DOTA into amides using a coupling agent such as BOP or HBTU. This method is particularly useful in cases where haloacetyl amides are hard to synthesize or in cases where quarternization of tertiary amines is a significant problem. (Scheme 5).173,268
A new solid-phase route to DOTA-peptide based ligands has been developed that utilizes an interesting DOTA derivative functionalized with a NH2-group on the alpha carbon of one of the acetate sidearms. The synthesis starts with the DO3A-derivative attached to a Wang resin through one of the carboxylate arms.269 A peptide chain was then grown from the amino group of the DOTA unit using standard SPPS chemistry followed by cleavage from the resin using TFA. Finally, Tm3+-complexation was performed in water at pH 5 to 7.5 (Scheme 6). The structure of the ligand was specifically designed so that the amide group on the acetate sidearm would be positioned in the proximity of the central Tm3+-ion in the final complex in order to achieve greater Ln3+-induced shift of the amide proton. This derivative was then used to detect enzyme activity as discussed in Section 184.108.40.206.
The Ln3+-complexes of neutral DOTA tetraamides are usually prepared by reaction of the ligand with equivalent amount of Ln3+-salts (anhydrous triflates are commonly used) in acetonitrile or methanol under anhydrous conditions.145 Due to slow formation rates, the complexation is often run for several days at 40 to 50°C. Ln3+-complexes of anionic tetraamides such as 51 are usually synthesized in water while maintaining the pH between 5–7.42,125,140,145,173,263–265
The synthesis and CEST properties of a bifunctional chelator for PARACEST agents has also been reported.270 Here, the ligand has an aromatic nitro group that can be reduced and converted to an active isothiocyanato group for later bioconjugation (structure 94). The Eu3+ complex of ligand 94 adopts predominantly (90%) the SAP geometry as reported by 1H NMR. Interestingly, CEST spectra of this nitro-derivative showed that alpha substitution on one side arm also impacts water exchange in the resulting Eu3+ complex. Water exchange was found to be slower in the bifunctional derivative when compared to the parent complex, 16.
Since PARACEST agents are ultimately intended for in vivo applications, their tolerance, biodistribution, toxicity and metal ion release characteristics must be critically evaluated. The stability of a metal complex is usually described by two constants, a thermodynamic stability constant and a kinetic rate constant that characterizes dissociation of the complex into a “free” ligand and metal ion. The latter is referred to as kinetic inertness. Since the amount of metal ion released from the complex is largely determined by its kinetic inertness, for in vivo applications this is the more important factor and the in vivo safety of Ln3+-complexes should be judged by their resistance towards dissociation. 271 The ligand protonation constants and thermodynamic stabilities of several Ln3+-DOTA tetraamide complexes have been studied in detail by pH-potentiometry.264,272–276 In the parent ligand DOTA, the negatively charged acetate sidearms stabilize the protonated macrocyclic nitrogen atoms through hydrogen bonding interactions and thus the first and second pKa’s of the ligand are relatively high.277 The neutral amide functionality is not capable of forming analogous hydrogen bonds with the protonated nitrogen atoms and therefore, the first two pKa values of all DOTA tetraamide ligands is about 3 logK units lower than those of DOTA (Table 5).275
The total basicity (the sum of pKa values) of the tetraamides is also significantly lower than that of DOTA and, for this reason, the thermodynamic stability of Ln3+-DOTA tetraamides is about ten orders of magnitude less than that of the corresponding DOTA complexes. There is only a slight variation in the stability constants with decreasing ionic size of the lanthanide ion along the series. The stabilities increase for the lighter lanthanides and then remain almost constant for the heavier ones with a small decrease seen at the very end of the series (Table 6). For each ligand, the Eu3+ complex appears to be the most thermodynamically stable of all lanthanide complexes, likely reflecting the best match between the size of coordination cavity of the ligand and this metal ion. It is worth noting that unlike DOTA, the stability of the Cu2+ and Zn2+ DOTA-tetraamide complexes is somewhat higher than that of the corresponding Ln3+-complexes. This may be of concern for in vivo applications since in this case the displacement of the Ln3+-ion by these endogenous metal ions is thermodynamically favorable. However, it has been shown that endogenous metal ions do not play a role in dissociation of Ln3+-chelates of DOTA and other DOTA-like ligands based on 12-membered macrocyclic tetra-amines (cyclen and pyclen) in vivo. In fact, the dominant pathway for release of metal ions in vivo for macrocyclic Ln3+-chelates involves the acid-catalyzed dissociation and the acid independent spontaneous dissociation of the complex. 278
The kinetic behavior, both formation and dissociation, of Ln3+-DOTA-tetraamides has been studied by UV-Vis spectrophotometry. Formation of the Ln3+-51 complexes is about 10-fold slower than that of the corresponding DOTA complexes. The formation reaction follows the same mechanism as the formation of DOTA complexes involving base-catalyzed deprotonation and rearrangement of a protonated intermediate in the rate determining step.275 The structure and stability of this intermediate has been studied in detail for DOTA complexes and it is now well established that the metal ion is coordinated to four oxygens of the acetate sidearms, two N-atoms of the macrocyclic ring and three water molecules.279–281 Although the structure of the intermediate complex has not been determined for Ln3+-51, the carboxylate groups in this case are located further away for the coordination cage, consistent with a lower stability of the intermediate. This in turn results in a slower rate of complex formation even though the basicity of the macrocyclic N atoms in 51 is significantly lower than those for 2.275
Interestingly, formation of Ln3+-complexes with neutral DOTA-tetraamides occurs considerably slower (about 2 to 3 order of magnitude) than that of the DOTA complexes and occurs through the direct reaction between the fully deprotonated ligand and the lanthanide ion (Table 7). The participation of mono and diprotonated species is negligible even at pH values where the concentration of the protonated ligand species is several orders of magnitude higher than that of the fully deprotonated ligand.276 The lack of a protonated intermediate is extremely surprising because the diprotonated complex Gd3+-H217 has been prepared and its structure in the solid state was determined.282 This X-ray structure shows that the Gd3+ ion is coordinated by four water molecules and four amide oxygens while two of the macrocyclic N atoms remain protonated. This may seem contradictory to the results of the kinetic measurements but one must realize that the existence or the crystallization of such a diprotonated intermediate does not necessarily imply that the formation of the final product occurs through this intermediate, rather, it can be a “dead-end” complex, which is not involved in the pathway leading to the final product.277
As already pointed out, the in vivo fate of a metal complex is largely dependent on its kinetic inertness. Considering the potential application of DOTA-tetraamide complexes as MRI agents, evaluation of their kinetic inertness is essential, especially in light of their low thermodynamic stability. Since Ln3+-chelates of DOTA-like ligands dissociate in vivo only by acid-catalyzed dissociation and acid-independent spontaneous dissociation, the endogenous metal ions (Cu2+ and Zn2+) do not affect the dissociation rates.278,280 The acid-catalyzed dissociation kinetics of some LnDOTA-tetraamides have been studied in strongly acid solutions, below pH 1 (Table 8).264,274–275 In this pH range, the complexes are thermodynamically unstable and the measurements could be performed in a reasonable amount of time. The kinetic data were collected in the presence of a large excess of acid and, under these pseudo first order conditions, the dissociation rates could be described by Eq 55.
The rates of dissociation were found to be linearly proportional to the acid concentration as described by Eq 56.
where k0 represents an acid-independent dissociation rate constant (spontaneous dissociation) and k1 represents acid-assisted dissociation of the protonated complex. This kinetic behavior is similar to that reported for the LnDOTA chelates and is consistent with a rapid formation of a monoprotonated intermediate with protonation of an amide oxygen as the initial step.280 Rearrangement of the protonated species transfers the proton from the acetate arm to a ring nitrogen atom and electrostatic repulsion between the protonated nitrogen and the metal ion results in the formation of an “out-of-basket” complex which then undergoes dissociation. A comparison of the k1 values reveals that Eu3+-51 dissociates about 15-fold faster than Eu3+-18, likely due to the fact that at this pH all four glycine carboxylates are protonated and proton transfer from the protonated glycinate carboxylates occurs relatively easily (Table 8). Even though the spontaneous dissociation of LnDOTA-tetraamides appears to be somewhat faster than dissociation of LnDOTA, acid-catalyzed dissociation of these complexes is significantly slower (up to 30 fold, depending on the structure) than dissociation of the analogous LnDOTA chelates. The slower acid-catalyzed dissociation of the LnDOTA-tetraamide complexes in general reflects the more difficult protonation of the neutral amide ligating arms which in turn improves kinetic inertness.
As judged by their favorable in vitro kinetic inertness, one would predict that the LnDOTA-tetraamides should be safe for in vivo use provided they are eliminated from the body as rapidly as the LnDOTA chelates. Biodistribution studies with radioactive 177 Lu3+-51, 153 Gd3+-22 and 177 Lu3+-95 demonstrated that negatively charged (anionic) complexes are excreted rapidly by renal filtration and are well-tolerated even at the relatively high concentrations required for MRI (>0.1 mmol/kg).283 However, 153 Gd3+-22 showed a marked preference for bone (33 % injected dose/organ), likely due to the bone targeting effect of the phosphonates. On the other hand, 177 Lu3+-95 in which phosphonate monoester groups were substituted for the phosphonate did not localize on bone. Cationic complexes of the ligands 16, 49, and 96 were also primarily excreted through the renal system and well tolerated at low doses (1–35 μmol). They showed slightly higher uptake in bone liver and spleen than 177 Lu3+-51. Since these organs are associated with free Ln3+ uptake as a consequence of in vivo dissociation, it is more likely that the complexes remain intact but, due to their 3+ charge, they mimic the biodistribution of a free Ln3+ ion. However, at higher doses (0.1 mmol/kg) the cationic complexes were acutely toxic, causing death within a minute of injection by interfering with heart function. Interestingly, the cationic complex, Lu3+-97 did not show any acute toxic effects when administered at a typical MRI dose.264 The exact reasons for this are not yet known but it was suggested that the decreased toxicity is due to the presence of the 12 peripheral alcoholic OH groups. This peculiar behavior may have significant implications for the design of PARACEST agents for in vivo applications because modification of a complex with a large number of peripheral alcoholic OH groups may lower its undesirable interaction with cell surfaces.
Free Ln ions are toxic because they act as a surrogate for Ca2+ ions and subsequently accumulate in bone and teeth. They can also form colloids with serum proteins that are eventually eliminated by the reticuloendothelial system.284 Free Gd3+ due to the dissociation of the contrast agent has also been implicated in Nephrogenic Systemic Fibrosis (NSF) although it has also been suggested that the intact complexes may play a role in triggering NSF. NSF is a systemic fibrosing disorder which characteristically exhibits fibrosis of the skin and the connective tissues, but it also affects other organs.285–287 It is important to emphasize that NSF is exceedingly rare and only affects a small percentage of patients with severely compromised renal function. Choice of contrast agent appears to have a great effect on the occurrence of NSF: there is growing evidence that NSF is associated with MRI contrast agents that exhibit low kinetic inertness, specifically Gd3+-complexes of open chain ligands derived from DTPA. So far, NSF has not been confirmed in patients who received macrocyclic contrast agents.287–288 The exceptionally high kinetic inertness of Gd3+ complexes of DOTA and DOTA-like ligands in general is the result of the preorganized, conformationally rigid structure of the macrocyclic cyclen backbone. Since Ln-complexes are eliminated by the kidney, renal failure leads to extended circulation times. The elimination half-life of Gd3+ chelates in healthy humans is around 1.5 hours, whereas it can be as long as 30 hours or more in patients with compromised renal function.239,289–291 If the kinetic inertness of the complex is not sufficiently high, the longer circulation times caused by kidney failure allow the metal-ligand system to approach equilibrium as defined by the conditional thermodynamic stability constant at a given pH. The released free Gd3+ then causes NSF.271,292–293 While it is erroneously accepted that high thermodynamic stability renders a complex stable in vivo, experimental evidence strongly suggests that the in vitro acid-catalyzed dissociation rates more accurately predict the fate of a complex in vivo. It was first reported 18 years ago that the extent of long term Gd3+-deposition in mice correlated well with the in vitro acid-catalyzed dissociation rates of the complexes. The deposition of Gd3+ was the lowest for Gd3+-2, which is in agreement with its extremely slow proton assisted dissociation rate.294
Since lanthanides other than Gd3+ can also trigger NSF, the in vivo stability and toxicity of Eu3+-51 has also been studied with regard to NSF.271 The pharmacokinetics of Eu3+-51 doped with tracer levels of 177 Lu3+-51 was followed in rats at 0.7 μmol/kg and 0.1 mmol/kg doses.271 The complex had an elimination half-life of 20 minutes and Eu3+-51 is completely eliminated from blood within 2 hours post-IV dosing. An analysis of residual 177 Lu3+-51 in tissue showed that the complex was found mainly in blood, kidney, and urine samples at 30 minutes, while at 2 hours the residual amount in these organs had decreased significantly. Such findings suggest rapid excretion of the complex by renal filtration. This pharmacokinetic behavior is very similar to that of other macrocyclic complexes such as Gd 3+-2. The chronic toxicity of Eu3+-51 was also evaluated in rats at 0.1 mmol/kg, 0.5 mmol/kg and 1.0 mmol/kg single doses, which are estimated to be required for in vivo applications. The compound was well tolerated in all cases and doses up to 1.0 mmol/kg were not associated with any toxicity. Even though the thermodynamic stability constant of Eu3+-51 is approximately 10 orders of magnitude lower than that of Gd3+ -2, Eu3+-51 dissociates about 25 times slower than Gd3+-2 (at pH = 1). While high thermodynamic stability combined with high kinetic inertness is certainly the best case scenario, thermodynamic stability is not the sole contributor to the safe application of these complexes in biological systems. Even with a low thermodynamic stability, the complex may dissociate so slowly that the amount of released metal ion is negligible before the complex is completely eliminated. Thus, kinetic inertness is a better predictor than thermodynamic stability for in vivo safety of such lanthanide complexes.271,292–293
In conclusion, PARACEST agents based on DOTA tetraamides do not exhibit overt toxicity at doses that would be required for their in vivo application (0.1–1.0 mmol/kg) and it should be emphasized that the biggest obstacle to their in vivo applications is not sensitivity or toxicity issues, but rather problems associated with inherent tissue MT effects (Section 220.127.116.11). Once these are solved, either by designing new complexes that can shift the bound water peak away from the MT window or by using improved imaging pulse sequences, all the in vitro studies of responsive agents presented in Section 4.4.7 should work in vivo.
In vivo MR imaging techniques based on magnetization transfer techniques have been in use for many years47 largely based on the presence of a semi-solid macromolecular component in tissue that gives a broad, underlying 1H NMR signal from water associated with the solid-like component that spans a radiofrequency range of ~100 kHz (Figure 57). It has been demonstrated that when a weak RF saturation pulse is applied at a particular frequency of this macromolecular phase, transfer of this magnetization to the bulk water signal is thought to occur by two mechanisms; (1) through dipolar coupling between the protons that have restricted motion due to entrapment of water or protons in macromolecular cavities and/or (2) through the chemical exchange between the water protons and macromolecular side chains containing –NH, -OH exchanging groups.45 It is difficult to distinguish between these two mechanisms. While the conventional MT effect of the semi-solid macromolecular phase in the tissues usually occurs over a large frequency range of ~100 kHz (determined by dipolar interaction, magnetic susceptibility and the chemical anisotropy of the macromolecules in the tissues), the CEST effects of DIACEST and PARACEST agents usually occur only over a very narrow frequency range. As recently observed by Zhou and van Zijl,12 one of the most significant differences between the inherent MT effects of biological tissues and CEST is the frequency specificity of a CEST agent as compared to the mainly symmetric appearance of the conventional MT proton signal.
One of the challenges of imaging the CEST effect from a diamagnetic CEST agent in vivo is the small Δω values of the exchanging sites in these contrast agents. As a result, there is a possibility that the MT effect from tissues can largely obscure CEST from an agent.12 Initial efforts to avoid this problem for in vivo imaging involved the collection of two MR images with a presaturation offset pulse applied at equal frequencies on either side of the bulk water resonance (Figure 57).65,81 While the MT signal from tissue is reasonably symmetrical about the free bulk water resonance, the contribution from a CEST agent should result in asymmetry in the CEST spectrum. Thus, by subtracting images collected at a positive offset frequency from another collected at a negative offset frequency, a true CEST image devoid of MT effects should be obtained. Using this approach, Guivel-Scharen et al. demonstrated that tissue metabolites can be detected and imaged via water protons using the signal amplification properties of saturation transfer in the presence of in vivo MT effects.65 Images of the medulla of the kidney were obtained using a difference imaging technique between a control irradiation (− 1.74 ppm) and on-resonance irradiation (+ 1.74 ppm) using the protons of urea that are in chemical exchange with the water proton in the tissues (Figure 58).
It is generally expected that the intrinsic problem of MT effects might be offset to some extent by using PARACEST agents owing to the large Δω values for their exchanging sites. In one of the earliest reports of imaging PARACEST agents in tissues, Li et al. described the temperature dependence of Eu 3+-57 in brain tissue.194 Their use of a two-dimensional Fast Low Angle Shot (FLASH) pulse sequence preceded by a continuous presaturation pulse showed consistency in the average temperature of the tissue despite the inherent MT effects of the tissues. While the phantom image and the corresponding CEST spectrum in aqueous solution shows the expected result, images of the agent mixed with brain tissue show significantly overlap with the inherent MT effects in tissue. Even with this overlap, they were able distinguish asymmetry in the CEST spectrum near 50 ppm corresponding to exchange of the bound water molecule on Eu3+-57 (Figure 59).
To further validate these exchange phenomena in biological systems, Li et al. described a four-pool model that basically combined the three-pool model described earlier by Woessner et al. with an extra pool to describe exchange between the macromolecular protons and bulk water, the MT pool. 67,295 Simulations of the CEST spectra in the presence and absence of macromolecular protons indicate that while the agent showed little variation in the CEST effect as a function of the exchangeable proton chemical shift Δω (Figure 60a), the net CEST effect decreased considerably for chemical shifts within ±100 ppm of the bulk water as a result of the interference from the macromolecular MT (Figure 60b). On the other hand, for Δω values greater than ± 100 ppm, the effects due to tissue MT effects showed minimal interference with CEST. As discussed previously (Table 3), some paramagnetic lanthanide complexes display extremely large shifts for exchangeable water/proton resonances, especially the bound water protons in some Tb3+, Tm3+, Dy3+ and Yb3+ complexes. Thus, while Eu3+ complexes have been most extensively studied so far owing to their slow exchange kinetics, it would be extremely beneficial to use lanthanide ions that induce larger shifts to improve CEST efficiency in vivo by avoiding the MT effects from endogenous macromolecules.
Li et al. analyzed the experimental CEST results obtained for their contrast agent in vivo using the four pool model and fit the spectra using Lorentzian, Gaussian and super-Lorentzian line shapes (Figure 61). The calculated values of bound water chemical shift using the super Lorentizian fitting were in reasonably good agreement with the experimental results. The average T2 of the protons associated with the endogenous macromolecules are also in excellent agreement with the values obtained in earlier measurements, giving validity to the use of the four-pool modeling system to quantify PARACEST spectra.296 Although the CEST effects of the bound water protons of the Eu3+ complex are still considerably quenched by the endogenous macromolecular MT effects, these results represent a foray into in vivo imaging using PARACEST agents and offer sufficient hope for further improvements.
Another method to get good temporal resolution in spite of the MT effects is the use of either a fast or slow selective presaturation time using interleaved repetitive short saturation pulses and multiple echo acquisition, as reported by Liu et al.297 The use of a short echo time and long saturation time is an ideal condition for the system to attain steady-state magnetization and to give the maximum contrast in vivo. The choice of imaging sequence is governed by the T1 of the tissue and the contrast agent. This report compares the effective improvement in PARACEST contrast between two imaging pulse sequences: one involving a short presaturation time with a long echo time using the RARE (Rapid Acquisition with Relaxation Enhancement) sequence suitable for environments with long T1, and the other involving a short presaturation time with a short echo time using the FLASH (Fast Low Angle Shot) pulse sequence that lengthens the total acquisition time, but offsets the loss of PARACEST contrast due to short T1.298–300 These modified pulse sequences provide major inroads into improving the temporal resolution and scan times of in vivo images.
While significant knowledge has been gained in these studies about in vivo imaging of PARACEST agents, there is much to be done before these agents will be used in a clinical setting. The inherent MT effect of tissues still remains as the main barrier towards the use of these agents for in vivo studies. Overcoming this problem will greatly facilitate the translation of these agents into in vivo applications, especially with respect to the responsive agents.
So far, we have discussed in great detail the use of paramagnetic lanthanide complexes as contrast agents using a new contrast enhancement technique called chemical exchange saturation transfer. While Gd3+ complexes with faster water exchange kinetics are ideal as T1 agents, other Ln3+ complexes with slow exchange kinetics have shown great promise as PARACEST agents. In this section we briefly discuss a few different types of lanthanide complexes that are viable substitutes to Gd3+-based T1 agents and PARACEST agents.
Many of the PARACEST agents discussed so far satisfy the slow-to-intermediate water exchange condition and are activated by saturation of a highly shifted Ln3+-bound water resonance. Using the modified Bloch equations, it can be predicted that the efficiency of saturation for CEST agents with faster water exchange can be significantly improved if a strong RF pulse is used for off-resonance saturation. This places obvious limitations on the use of conventional off-resonance saturation techniques for these fast exchanging systems owing to the fact that the power required to achieve saturation of the exchanging pool would likely exceed the acceptable specific absorption rate (SAR) limit for humans. While we have discussed the use of off-resonance saturation transfer as one of the techniques to study chemical exchange effects on slow exchanging systems as a major part of this review, many other techniques have been used to study slow exchanging processes.58,82,301–302 An alternative way to detect chemical exchange in fast exchanging paramagnetic lanthanide complex was recently demonstrated by Vinogradov et al. using a method called On Resonance Paramagnetic Agent Chemical Exchange Effect (OPARACHEE).303 The basic hypothesis is that it is possible to use exchange and relaxation-sensitive pulses directly on the bulk water resonance to observe contrast due to chemical exchange. 304–307 The principle behind OPARACHEE can be qualitatively understood as follows: application of a long, low power 360° pulse (or equivalent) on the bulk water resonance should result in a full 360 rotation of the net magnetization in the absence of exchange. But when chemical exchange occurs, some of the exchanging protons “leak away” and do not experience the full 360 degree rotation. This effect will result in a decrease of the Z-magnetization and negative contrast enhancement, quite similar to conventional CEST. In practice, this is accomplished by using a modified WALTZ-16 pulse train (WALTZ-16*), where the pulse length is fixed at 2.5 ms, which does not necessarily correspond to calibration of the 90 pulse as in the original WALTZ-16 sequence.308–309 The effectiveness of this method is dependent on several parameters, including the longitudinal and transverse relaxation times, water residence times of both the bound and bulk water pools, the chemical shift difference between the bound and bulk water resonances (Δω), and concentration of the paramagnetic agent. Simulations of the magnitude of the magnetization ratio as a function of the bulk water residence lifetime and Δω suggest that maximal effect is attained when the bulk water residence time is ~0.5−1 s and Δω is as large as possible. In practice this means that Tm 3+ and Dy3+ complexes of DOTA should in principle perform best and this was indeed demonstrated in spectroscopy and imaging experiments.
An aqueous solution of Tm3+-17 at a concentration of 12 μM showed that a 4% reduction in bulk water can be achieved using a relatively low power B1 pulse of only 200 Hz. These results were also confirmed in an imaging experiment with three tubes of phantom containing different concentrations of Tm3+-17 (Figure 62). The lowest concentration (12.5 μM) that could be seen in the spectroscopy experiment could not be differentiated with respect to bulk water alone in the phantoms. This can be attributed to the expected contribution of B0 inhomogeneity in the imaging experiments. The same experiments performed with Dy3+-17 resulted in a 12 % decrease in bulk water signal intensity using 125 μM complex.
To evaluate the practical advantages of OPERACHEE over the “off resonance” presaturation CEST imaging techniques, mice were imaged after injection of Tm3+-51 using a WALTZ-16* pulse for excitation of the water signal followed by a standard spin echo acquisition. The kidney tissues appeared bright in such images due to the short T2 of this tissue; further delineation of different tissue types within the kidney could be achieved due to the difference in the magnetization transfer effect and the relaxation times. These results show that OPERACHEE could be used as an alternate imaging technique using other types of lanthanide complexes such Tm3+ and Dy3+ that have a fast bound water exchange and short relaxation times (Figure 63). Some advantages of OPARACHEE over CEST for imaging the exchange effects of such agents can be summarized as follows: 1) the low power WALTZ-16* pulse train causes much less RF power deposition than conventional CEST imaging, 2) the method expands the range of usable Ln3+-complexes to include fast exchanging systems that cannot be pre-saturated in the more typical CEST experiment, 3) the detection limit appears to be about an order of magnitude lower and 4) the need to know the exact frequency of Ln3+-bound water molecule is removed. A major disadvantage, however, is the high sensitivity of WALTZ-16* to short relaxing species, which removes some tissue signals from the image even in the absence of the paramagnetic complex. These experiments reveal that the OPERACHEE imaging sequence for fast exchanging systems presents a feasible alternative to the more traditional CEST method. The use of low B1 power levels of 100 to 300 Hz with agent concentrations as low as 30–200 μM greatly increase the viability of these systems when compared to other off resonance techniques.
Europium is the only lanthanide that exhibits extensive chemistry in the divalent state in aqueous solution. Eu2+ and Sr2+ have almost exactly the same ionic radius and the coordination chemistry of these two ions is quite similar. Eu2+ is paramagnetic and has the same 4f7 electron configuration as Gd3+ (isoelectronic). Given the more favorable water exchange kinetics of Eu2+ compared to Gd3+, this divalent ion may present an attractive alternative to Gd3+ for certain T1 contrast agents applications. Eu2+ does however have a shorter electronic relaxation time than Gd3+; so the relaxivity of Eu2+ complexes are expected to be somewhat lower than the corresponding Gd3+ complexes. The Eu2+ aqua ion has an average coordination number of 7.2, suggesting that 7 and 8 coordinate species exist in equilibrium. This means that the free energy difference between the ground and transition state is very low resulting in an extremely fast exchange rate, in fact, the fastest ever measured by 17O NMR.310 Activation volume measurements have shown that water exchange proceeds by a limiting associative mechanism (A). Water exchange in Eu2+-polyamino polycarboxylate complexes is also generally much faster than the corresponding Gd3+-complexes due to the larger ionic radius and smaller ionic charge of Eu2+. From 17O NMR measurements, the water exchange in Eu2+-2 follows an interchange mechanism (Ia) with an activation volume of +0.1 ± 1.0 cm3mol−1. The extremely fast water exchange rate of Eu2+-2 is likely due to the longer than usual Eu2+-O (bound water) distance in this complex, which adopts a TSAP geometry.
Practical applications of Eu2+ complexes as T1 agents are severely hampered by the extreme sensitivity of this ion toward oxidation. This chemical feature makes Eu2+ an attractive platform for creating T1-based redox sensors. This would be especially advantageous because the oxidized product, Eu3+, has negligible relaxation effects. Generally, the sensitivity of Eu2+ complexes of poly-amino poly-carboxylate ligands (DTPA, DOTA) to oxidation is even greater than the Eu2+ aqua ion. Ligands based on the azacrown ether 1,4,10,13-tetraoxa-7,16-diazacyclooctadecane are promising but the complexes still have somewhat more negative redox potential than the aqua ion. Water exchange in Eu2+-98 proceeds by an associative mechanism (A) because there is enough room in this complex for another water molecule to coordinate. So far, the least sensitive complex is the cryptate, Eu2+-99. The variable temperature 17O NMR data fittings revealed that this complex has the slowest water exchange rate when compared with the earlier reported complexes.311 The activation volume indicates that the water exchange occurs via an interchange mechanism (Ia) where the rate determining step involves the incoming water molecule (Table 9). While biomedical applications of Eu2+ complexes may not ultimately prove feasible, the data collected in these studies could improve our understanding of how electronic relaxation mechanism and water exchange rate contributes to proton relaxivity beyond that already learned in studies of Gd3+ complexes.
As explained in our discussion of the paramagnetic properties of lanthanide ions in Section 3, conventional Gd3+-based contrast agents have longer electronic relaxation times due to their isotropic electronic ground state and as a result proton relaxation is reasonably well-described by Solomon Bloembergen Morgan theory. These agents affect both T1 and T2. Most of the remaining paramagnetic lanthanide ions, especially Dy3+, Pr3+, Sm3+, Ho3+, Er3+ and Yb3+ that are characterized by shorter electronic relaxation times, relax protons via a Curie mechanism that primarily affects T2.312 The contribution of Curie relaxation mechanism increases substantially with the external magnetic field B0 and is proportional to the square of the magnetic moment of the lanthanide ion. 23,313 Among the trivalent lanthanide ions, Dy3+ has the shortest electronic relaxation time (~0.5 ps) and highest magnetic moment (10.6 Bohr magnetons) which in combination results in very efficient T2 relaxation. Vander Elst et al.314 reported a series of Dy3+ complexes based on 1 (101–104) and showed that at low field that the transverse relaxivities of the complexes are similar to their longitudinal relaxivities, while at high field there is a marked dependence of r2 on both the magnetic field and τM. According to theoretical considerations, the low field transverse relaxivity of a Dy3+-complex is determined by its electronic relaxation time and should therefore be small and similar to its longitudinal relaxivity. At higher magnetic fields, however, the Curie contribution becomes more significant and the transverse relaxivity increases with the square of the external magnetic field and attains an optimal value when τM is between 0.1 and 1 μsec. Longer or shorter bound water residence lifetimes are detrimental to r2 (Figure 64). This was clearly demonstrated by field dependent proton r2 relaxivity measurements on Dy3+ complexes of 1, DTPA bisamides and 16. At low field, all Dy3+ complexes showed small and comparable r2 values but at higher fields, the r2 of Dy3+-1 and Dy3+-16 increased only slightly because of their faster and slower than optimal exchange rate respectively while the r2 of Dy3+-101 showed a marked increase with the square of B0 (Table 10). Imaging experiments confirmed the much more efficient T2 shortening effect of Dy3+-101 compared to Dy3+-1: signal attenuation in the sample containing 10 mM of Dy3+-101 could already be observed on the first echo image and on the fifth echo the signal could no longer be detected (Figure 65).314
The Curie spin contribution to r2 is also dependent on the rotational correlation time τR of the complex. This was demonstrated by Caravan et al. who measured the r2 of Dy3+-105 in the presence and absence of albumin.315 As reported earlier, Gd3+-106 has a high affinity for HSA. Dy3+-106 showed a 3–8 times increase r2 in the presence of HSA due to the increase in τR upon binding. Field dependent relaxivity measurements on Dy3+-106 indicated that the r1 of Dy3+-106 is largely field independent, while r2 increases with increasing field because the contribution of the Curie spin effect becomes larger as theoretically predicted. On the other hand, the r1 and r2 of Dy3+-105, a related complex with no metal bound water, did not change with the magnetic field. This indicates, as expected, the Curie relaxation effect of the Dy3+ is transferred to the bulk water through the exchange of the metal bound water molecule.
Dy3+-1 was also proposed for use in MRI as a T2* negative contrast agent to differentiate between healthy and damaged tissues. The complex produces a susceptibility effect due to its localization in extracellular compartments of healthy tissues. In damaged tissues where cell membranes are disrupted, the complex is distributed homogenously. As a result of compartmentalization, a local field gradient is generated because of the large magnetic susceptibility of the metal ion, which in turn results in a decrease in T2* and increasing negative contrast.316–317 However, since small chelates such as Dy3+-1 extravasate from the vascular space quite rapidly, a PAMAM dendrimer based Dy3+ complex with slower clearance rate has also been proposed as blood pool susceptibility agent.318
Although Gd3+ complexes have been extremely successful in altering contrast in T1-weighted images in both basic research and in clinical applications, the new era of molecular imaging makes more demands on MRI that cannot be satisfied with current generation Gd3+ agents. Over the past several years, chemical exchange saturation transfer (CEST) agents that make use of slowly exchanging –NH, -OH or bound water protons have evolved as a potentially new class of imaging agents and a promising alternative to Gd3+-based contrast agents. Even though chemical exchange principles have been widely described in the NMR literature, it was the pioneering work of Ward and Balaban on the application of small molecule diamagnetic molecules for CEST that stimulated new interest in developing such agents for molecular imaging of biological processes. The discovery of paramagnetic lanthanide DOTA-tetraamide complexes with extremely slow water exchange kinetics has further stimulated new ideas about CEST-based contrast agents. PARACEST agents have significant advantages over their diamagnetic counterparts, mainly due to the large lanthanide induced frequency shifts of the proton exchange sites in these molecules. One of the main reasons for the explosive growth of PARACEST agents can be attributed to the fact that so much has already been studied in great detail about water exchange mechanisms for Gd3+ and its complexes. This knowledge has proven extremely useful to translate most of the known concepts involving Gd3+ to PARACEST. In addition, the fact that virtually any paramagnetic lanthanide ion (except Gd3+) can be used for CEST imaging makes these agents extremely attractive as multi-frequency contrast agents to report multiple biological parameters such as pH, redox, temperature or enzyme activity in a single imaging experiment. It would be difficult to conceive such experiments with conventional T1 or T2 agents. Finally, arguably the most attractive feature of PARACEST is the ability to turn image contrast on or off with an RF pulse, a feature that is impossible with the never silent Gd3+ agents.
As water exchange is extremely sensitive to many types of environmental or structural changes, these agents have created new opportunities for developing responsive PARACEST agents that signal changes in biological or physiological events by MRI. Several successful in vitro experiments using novel responsive PARACEST agents provide ample evidence that these agents show great promise for reporting physiological phenomena using existing imaging technologies. Initial in vivo bio-distribution experiments suggest the kinetic inertness of these complexes will result in low toxicity for these agents, similar to their Gd3+ counterparts. While considerably more must be done to realize the full potential of these contrast agents in clinical settings, it can be envisioned that PARACEST agents are paving the way for MRI to become more than a tool for anatomical imaging but a true contender for molecular imaging as well.
The authors acknowledge financial support from the National Institutes of Health (CA-115531, CA-126608, RR-02584 and EB-004582) and the Robert A. Welch Foundation (AT-584) throughout the writing of this review.
Subha Viswanathan earned her B.Sc. from Stella Maris College, Madras University, Chennai in 2000. In 2007, she earned a Ph.D in Inorganic Chemistry from Syracuse University under the guidance of Dr. Ana de Bettencourt-Dias. While there she obtained a doctoral prize for her work on the synthesis and characterization of nitro- and oxazoline-derivatized ligands for lanthanide complexes. She was also honored with the William D. Johnson award for Outstanding Teaching Assistant in 2006. In 2007 she came to the Advanced Imaging Research Center at UT Southwestern Medical Center to work as a postdoctoral researcher for Dr. Dean Sherry. Her primary research interests include the synthesis and characterization of peptide-based macrocyclic lanthanide complexes for PARACEST.
Zoltán Kovács was born in Kisvárda, Hungary in 1963. He received his Ph.D from Lajos Kossuth University, Debrecen, Hungary (now University of Debrecen) in 1992 under the direction of Dr. Béla Győri, working on cyano- and carboxy-borane complexes. Immediately thereafter he moved to the lab of Dr. Dean Sherry at the Department of Chemistry, University of Texas at Dallas as a post-doctoral researcher and later as Research Scientist to conduct research on the synthesis of bifunctional ligands for lanthanide complexes. He also worked as a Senior Research Scientist for Macrocyclics, Dallas, Texas. In 2006, he moved to the Advanced Imaging Research Center at UT Southwestern Medical Center as Assistant Professor where his primary research interests include the synthesis and application of hyperpolarized yttrium-89 complexes and carbon-13 labeled metabolic intermediates and the design and development of new ligand systems for PARACEST, T1 and radiopharmaceutical imaging.
Kayla N. Green received her B.S. degree from Tarleton State University in 2003. In 2007, she obtained her Ph.D. in Inorganic Chemistry from Texas A&M University under the direction of Marcetta Y. Darensbourg where she studied immobilized biomimetic complexes of metalloenzymes. While there she was she was honored with the U.S. Senator Phil Gramm Doctoral Fellowship Award (2007) and a Chemistry Biology Interface - NIH Training Grant (2004). In 2008 she began her postdoctoral work at UT Southwestern Medical Center. Her research currently focuses on the synthesis of lanthanide complexes for molecular imaging using novel NMR techniques and the development of computational methods for the evaluation of these complexes.
S. James Ratnakar earned his Bachelors in Science (BS) in Chemistry (April 1995) from Madras University, India and his Masters in Science (MS) in Chemistry (May 1997) Bharathidasan University, India. He received his Ph.D in July 2004 from Madras University, India for his research on the synthesis and relaxivity studies of gadolinium (III) complexes of polyazapolycarboxylate macrocycles, under the guidance of Dr. V Alexander at Loyola College Chennai, India. After a brief postdoctoral work with Dr. V. Alexander, in 2006 he joined Dr. Dean Sherry at the Advanced Imaging Research Center at UT Southwestern Medical School, Dallas, Texas. The broad area of his research interest is design and synthesis of macrocyclic ligands and their lanthanide complexes for use as imaging agents mainly for MRI.
A. Dean Sherry earned a PhD in Inorganic Chemistry at Kansas State University and won a National Institutes of Health Postdoctoral Fellowship to study metal-protein interactions before joining the chemistry faculty at UT-Dallas in 1972. He served as Chairman of the UT Dallas chemistry department from 1979–1990. In 1990, he accepted an appointment as Professor of Radiology at UT-Southwestern Medical Center where, for the next 15 years, he spent 50% of his time help build a metabolism MR program and 50% time teaching chemistry and doing research at UT Dallas. He was recognized for his outstanding achievements in chemistry by winning the Doherty Award from the DFW Section of the American Chemical Society in 1990 and was honored by receiving the Chancellor’s Outstanding Teaching Award at UT-Dallas in 1994. He currently holds a Cecil & Ida Green Distinguished Chair in Systems Biology at UT Dallas and serves as Director of the Advanced Imaging Research Center (AIRC) on the campus of UT Southwestern Medical Center. He and his collaborators have published ~320 peer-reviewed research articles and his research has been funded by the National Institutes of Health and a variety of private companies and foundations. He is the Scientific Founder of Macrocyclics, a local biomedical technology & development company that supplies bifunctional chelates and other key research intermediates to the biomedical imaging and nuclear medicine imaging & therapy communities.
Dedicated to the memory of Dr. Donald E Woessner who guided us in much of our science.