Open access peer-reviewed chapter

Basics of Chemical Exchange Saturation Transfer (CEST) Magnetic Resonance Imaging

Written By

Kenya Murase

Submitted: June 10th, 2017 Reviewed: October 12th, 2017 Published: December 20th, 2017

DOI: 10.5772/intechopen.71645

Chapter metrics overview

1,977 Chapter Downloads

View Full Metrics

Abstract

Chemical exchange saturation transfer (CEST) is one of the contrast mechanisms in magnetic resonance imaging (MRI) and has been used to detect dilute proteins through the interaction between bulk water and labile solute protons. Amide proton transfer (APT) MRI has been developed for imaging diseases such as acute stroke. Moreover, various CEST agents have been explored to enhance the CEST effect. The contrast mechanism of CEST or APT MRI, however, is complex and depends not only on the concentration of amide protons or CEST agents and exchange properties, but also varies with imaging parameters such as radiofrequency (RF) power and magnetic field strength. When there are multiple exchangeable pools within a single CEST system, the contrast mechanism of CEST becomes even more complex. Numerical simulations are useful and effective for analyzing the complex contrast mechanism of CEST and for investigating the optimal imaging parameter values. In this chapter, we present the basics of CEST or APT MRI and a simple and fast numerical method for solving the time-dependent Bloch-McConnell equations for analyzing the behavior of magnetization and/or contrast mechanism in CEST or APT MRI. We also present a method for analyzing the behavior of magnetization in spin-locking CEST MRI.

Keywords

  • Bloch-McConnell equations
  • numerical solution
  • chemical exchange saturation transfer (CEST) MRI
  • amide proton transfer (APT) MRI
  • spin-locking

1. Introduction

Chemical exchange saturation transfer (CEST) is one of the contrast mechanisms in magnetic resonance imaging (MRI) [1] and has been increasingly used to detect dilute proteins through the interaction between bulk water protons and labile solute protons [2, 3, 4]. Amide proton transfer (APT) MRI has been developed for imaging diseases such as acute stroke and cancer, and is now under intensive evaluation for clinical translation [5, 6]. APT MRI is a particular type of CEST MRI that specifically probes labile amide protons of endogenous mobile proteins and peptides in tissue [5, 6]. In addition to APT MRI [5, 6], useful CEST MRI contrast for clinical imaging can be generated from amine protons [7], hydroxyl protons [8], glycosaminoglycans [9], and glutamate [10], as well as from changes in creatine and lactate concentrations [11]. Glucose and iopamidol have been used as exogenous CEST agents that have been administered to patients [12, 13]. Moreover, various CEST agents have been energetically developed to detect the parameters that reflect tissue molecular environment such as hydrogen ion exponent (pH) and/or to enhance the CEST effect [14].

In CEST or APT MRI, the exchangeable proton spins are saturated, and the saturation is transferred upon chemical exchange to the bulk water pool [1, 15]. As a result, a large contrast enhancement in bulk water can be achieved. The contrast mechanism of CEST or APT MRI, however, is complex and depends not only on the concentration of amide protons or CEST agents, relaxation, and exchange properties but also varies with imaging parameters such as radiofrequency (RF) power and magnetic field strength [15]. When there are multiple exchangeable pools within a single CEST system, the contrast mechanism of CEST becomes all the more complex [16]. Numerical simulations are useful and effective for analyzing the complex CEST contrast mechanism and for investigating the optimal imaging parameter values [17, 18]. In order to perform extensive numerical simulations for CEST or APT MRI, it requires the development of a simple and fast numerical method for obtaining the solutions to the time-dependent Bloch-McConnell equations.

In this chapter, we present the basics of CEST or APT MRI and a simple and fast numerical method for solving the time-dependent Bloch-McConnell equations for analyzing the behavior of magnetization and/or contrast mechanism in CEST or APT MRI. We also present it in SL CEST MRI.

Advertisement

2. Bloch-McConnell equations in the presence of CEST

2.1. Two-pool chemical exchange model

A two-pool chemical exchange model is illustrated in Figure 1 . A and B in Figure 1 represent the pools of bulk water protons and labile solute protons, respectively. The time-dependent Bloch-McConnell equations for the two-pool chemical exchange model in CEST or APT MRI are expressed as [17, 18].

dMxatdt=R2aMxatkabMxat+ωaMyatω1yMzat+kbaMxbtdMyatdt=ωaMxatR2aMyatkabMyat+ω1xMzat+kbaMybtdMzatdt=ω1yMxatω1xMyatR1aMzatkabMzat+kbaMzbt+R1aM0adMxbtdt=kabMxatkbaMxatR2bMxbt+ωbMybtω1yMzbtdMybtdt=kabMyatωbMxbtR2bMybtkbaMybt+ω1xMzbtdMzbtdt=kabMzat+ω1yMybtω1xMybtR1bMzbtkbaMzbt+R1bM0bE1

Figure 1.

Two-pool chemical exchange model.

where superscripts a and b show the parameters in pool A and pool B, respectively. For example, Mxat, Myat, and Mzatare the x, y, and z magnetization components in pool A at time t, respectively. R1aand R2aare the reciprocals of the longitudinal (T1a) and transverse relaxation times (T2a), that is, the longitudinal and transverse relaxation rates in pool A, respectively. kaband kbadenote the exchange rate from spins in pool A to those in pool B and that from spins in pool B to those in pool A, respectively ( Figure 1 ). M0aand M0bare the thermal equilibrium z magnetization components in pool A and pool B, respectively. Δωa = ωa − ωand Δωb = ωb − ω, where ωa, ωb, and ωdenote the Larmor frequencies in pool A and pool B, and the frequency of the RF-pulse irradiation, respectively. ω1xand ω1yare the x and y components of the amplitude of the RF-pulse irradiation (ω1), respectively. Note that ω1 = γB1, where γand B1 are the gyromagnetic ratio (γ/2π = 42.58 MHz/T) and RF power, respectively. When the RF pulse is applied along an angle φfrom the x-axis of the rotating frame as illustrated in Figure 2 , ω1xand ω1yare represented by ω1x = ω1 cos ϕand ω1y = ω1 sin ϕ, respectively. When the RF pulse is applied along the x-axis of the rotating frame, ω1xand ω1ybecome ω1 and 0, respectively.

Figure 2.

Parameters for analyzing the behavior of magnetization in the rotating frame.

The differential equations given by Eq. (1) can be combined into one vector equation (homogeneous linear differential equation) [18]:

dMtdt=Aωω1ϕMt,E2

where

Mt=MxatMyatMzatMxbtMybtMzbt1TE3

and

Aωω1ϕ=R2a+kabωaω1ykba000ωaR2a+kabω1x0kba00ω1yω1xR1a+kab00kbaR1aM0akab00R2b+kbaωbω1y00kab0ωbR2b+kbaω1x000kabω1yω1xR1b+kbaR1bM0b0000000.E4

Tin Eq. (3) denotes the matrix transpose.

For simplicity, we assume that the RF pulse is applied along the x-axis of the rotating frame, that is, ϕ = 0. According to Koss et al. [19], the matrix A(ω, ω1, 0) can be given by.

Aωω10=EC00,E5

where Eis the evolution matrix and Cis the constant-term matrix. Furthermore, Eis given by.

E=R+K.E6

In the case of Agiven by Eq. (4), Ris reduced to.

R=Ra00Rb,E7

where

Ra=R2aωa0ωaR2aω10ω1R1a,E8

and

Rb=R2bωb0ωbR2bω10ω1R1b,E9

Kin Eq. (6) is given by

K=kabkbakabkbaI,E10

where Iis a 3-by-3 identity matrix and ⊗ denotes the Kronecker tensor product. Cin Eq. (5) is given by.

C=R1aM0aR1bM0bT001T.E11

The solution of Eq. (2) with ϕbeing 0 can be given by [18].

Mt=eAωω10tM0,E12

where trepresents the so-called saturation time and M(0) is the matrix of initial values at t = 0. eA(ω, ω1, 0)tis the matrix exponential.

It should be noted that mass balance imposes the following relationship between the exchange rates (kaband kba) of pool A and pool B [17]:

kab=kab+kbaM0bM0a+M0bE13

and

kba=kab+kbaM0aM0a+M0bE14

2.2. Three-pool chemical exchange model

Figure 3 illustrates a three-pool chemical exchange model in which pool a represents the bulk water pool. In this case, Rand Kare given by [19].

R=Ra000Rb000RcE15

Figure 3.

Three pool chemical exchange model.

and

K=kabkackbakcakabkba0kac0kacI,E16

respectively. Rc in Eq. (15) is given by Eq. (8) in which the subscript aand superscript aare replaced by c. Cis given by.

C=R1aM0aR1bM0bR1cM0cT001T.E17

The solutions of other multi-pool chemical exchange models such as an hour-pool chemical exchange model are described in Ref. [20].

2.3. Calculation of Z-spectrum, MTRasym, and PTR

The CEST effect has usually been analyzed using the so-called Z-spectrum [18]. The Z-spectrum is given by the following equation:

Zspectrum=MzaωoffM0a,E18

where MzaΔωoffis the z magnetization component of bulk water protons (pool A) at Δωoff. Note that Δωoff =  − Δωa.

The magnetization transfer asymmetry (MTRasym) analysis has been performed using the following equation [18]:

MTRasym=MzaωoffMzaωoffM0a.E19

Instead of MTRasym, the following equation for proton transfer ratio (PTR) has also been used for analyzing the CEST effect [18]:

PTR=MzaωoffMzaωoffMzaωoff,E20

where MzaΔωoffdenotes the z magnetization component of pool A at the opposite side of the water resonance (Δωoff).

Figure 4(a) shows Z-spectra as a function of offset frequency (Δωoff) for various saturation times (0.5, 1, 2, 5, and 10 s) in the two-pool chemical exchange model ( Figure 1 ). Figure 4(b) shows Z-spectra as a function of Δωofffor various ω1 values (25, 50, 100, 150, and 200 Hz). It should be noted that because B1 = ω1/γ, ω1 values of 25, 50, 100, 150, and 200 Hz correspond to B1 values of 0.59, 1.17, 2.35, 3.52, and 4.70 μT, respectively. Figure 4(c) shows Z-spectra as a function of Δωofffor various M0b/M0avalues (1/500, 1/250, 1/125, 1/100, and 1/50).

Figure 4.

Z-spectra as a function ofΔωofffor various values of saturation timea,ω1b,andM0b/M0a(c) in the two-pool chemical exchange model.

In the above simulations, we assumed that T1aand T2awere 3 s and 100 ms, respectively, and T1b=1sand T2b=15ms[16]. The chemical shift of protons in pool B was set to be 4 ppm. It should be noted that the chemical shift of 4 ppm corresponds to Δωoffof 1192.8 Hz for the magnetic field strength of 7 T. Unless otherwise indicated, kab + kbawas assumed to be 100 Hz. M0aand M0bwere assumed to be 1 and 1/250, respectively. The saturation time and ω1 were taken as 2 s and 100 Hz, respectively. The matrix exponential and Kronecker tensor product were calculated using the MATLAB® functions “expm” and “kron,” respectively.

The peaks at 0 Hz (0 ppm) and 1192.8 Hz (4 ppm) in Figure 4 derived from pool A and pool B, respectively. As shown in Figure 4(a) and Figure 4(b) , Z-spectra changed largely depending on the saturation time and ω1, that is, Z-spectra became broad and tended to saturate with increasing saturation time and ω1. As shown in Figure 4(c) , the peaks at 1192.8 Hz increased with increasing M0b/M0avalue.

Figure 5 shows cases for the three-pool chemical exchange model ( Figure 3 ) consisting of bulk water (pool A) and two labile proton pools (pool B and pool C). In these cases, we assumed that T1a=3s, T2a=100ms, T1b=T1c=1s, and T2b=T2c=15ms[16]. The chemical shifts of two labile proton pools were set to be 4 ppm (Δωoff = 1192.8 Hz for the magnetic field strength of 7 T) and 5 ppm (Δωoff = 1491.0 Hz for 7 T). Unless otherwise indicated, kab + kba, kac + kca, and kbc + kcbwere assumed to be 100 Hz, 300 Hz, and 100 Hz, respectively. M0a, M0b, and M0cwere assumed to be 1, 1/250, and 1/500, respectively. The saturation time and ω1 were taken as 5 s and 50 Hz, respectively.

Figure 5.

Z-spectra as a function ofΔωofffor various values of saturation timea,ω1b,andM0c/M0a(c) in the three-pool chemical exchange model.

Figure 5(a) shows Z-spectra as a function of Δωofffor various saturation times (0.5, 1, 2, 5, and 10 s). The peaks at 0 Hz (0 ppm), 1192.8 Hz (4 ppm), and 1491.0 Hz (5 ppm) derive from pool A, pool B, and pool C, respectively. As shown in Figure 5(a) , Z-spectra changed largely depending on the saturation time, that is, Z-spectra became broad and tended to saturate with increasing saturation time. Figure 5(b) shows Z-spectra as a function of Δωofffor various ω1 values (25, 50, 100, 150, and 200 Hz). As in Figure 4(b) , Z-spectra became broad with increasing ω1 value. Figure 5(c) shows Z-spectra as a function of Δωofffor various M0c/M0avalues (1/500, 1/250, 1/125, 1/100, and 1/50). The peaks at 1491.0 Hz increased with increasing M0c/M0avalue.

Figure 6(a) shows the MTRasym values given by Eq. (19) as a function of ω1 for various saturation times (0.5, 1, 2, 5, and 10 s) in the two-pool chemical exchange model ( Figure 1 ), whereas Figure 6(b) shows those as a function of saturation time for various ω1 values (25, 50, 100, 150, and 200 Hz). As shown in Figure 6(a) , when ω1 was small, MTRasym tended to increase with increasing ω1 and saturation time. However, when ω1 was large, MTRasym tended to saturate or decrease with increasing ω1 value, depending on the saturation time. As shown in Figure 6(b) , MTRasym tended to saturate with increasing saturation time for all ω1 values.

Figure 6.

(a) MTRasym values as a function ofω1 for various saturation times in the two − pool chemical exchange model. (b) MTRasym values as a function of saturation time for various ω1 values.

Figure 7(a) shows the PTR values given by Eq. (20) as a function of ω1 for various saturation times (0.5, 1, 2, 5, and 10 s) in the two-pool chemical exchange model ( Figure 1 ), whereas Figure 7(b) shows those as a function of saturation time for various ω1 values (25, 50, 100, 150, and 200 Hz). As shown in Figure 7 , although PTR showed almost the same tendency with MTRasym ( Figure 6 ), the change in the PTR value depending on the saturation time or ω1 was larger than that in the MTRasym value.

Figure 7.

(a) PTR values as a function ofω1 for various saturation times in the two − pool chemical exchange model. (b) PTR values as a function of saturation time for various ω1 values.

In this study, we presented a simple equation for solving the time-dependent Bloch-McConnell equations, in which our previous method [18] and the approach presented by Koss et al. [19] were combined. Our method can be easily expanded to multi-pool chemical exchange models by modifying the matrix Ain Eq. (2). We previously reported that the solutions obtained by our method agreed with the analytical solutions given by Mulkern and Williams, [21] and the numerical solutions obtained using a fourth/fifth-order Runge–Kutta-Fehlberg (RKF) algorithm [18], indicating the validity of our method. In addition, our method considerably reduced the computation time as compared with the RKF algorithm [18]. These results suggest that our method will be useful in calculating the parameters such as the exchange rate of CEST agents using the non-linear least-squares fitting method [17].

As previously described, the so-called Z-spectrum has usually been used to analyze the CEST effect [18]. The Z-spectrum is obtained by plotting the z magnetization component of bulk water protons (Mza) in the form of Mzaversus Δωoff[Eq. (18)]. Figure 4(a) and Figure 5(a) showed that the saturation time affected the Z-spectra, and the CEST effect increased and saturated with increasing saturation time. The fact that the CEST effect saturates with increasing saturation time is more clearly confirmed by the relationship between MTRasym or PTR, and the saturation time shown in Figure 6(b) or Figure 7(b) . As shown in Figure 4(b) and Figure 5(b) , ω1 also affected the Z-spectra. Although the CEST effect increased with increasing ω1 value, the separation among peaks in the Z-spectrum plots degraded with increasing ω1 value. The influence of ω1 on the CEST effect is also clearly demonstrated by the relationship between MTRasym and PTR, and ω1 shown in Figure 6(a) or Figure 7(a) . The use of large ω1 may directly saturate bulk water protons, causing the so-called spillover effect [18]. The results shown in Figures 4 7 suggest that the values of imaging parameters in CEST MRI such as the saturation time and ω1 must be determined in consideration of both the CEST effect and spillover effect. Our method is useful for determining the optimal values of imaging parameters in CEST MRI.

2.4. Calculation of R1ρand R2ρ

The longitudinal relaxation rate in the rotating frame (R1ρ) can be obtained from the negative of the largest (least negative) real eigenvalue (λ1) of the matrix Ain Eq. (2), that is, R1ρ =  − λ1 [19, 22].

The transverse relaxation rate in the rotating frame (R2ρ) can be obtained from the absolute value of the largest real part of the complex eigenvalue (λ2) of the matrix Ain Eq. (2), that is, R2ρ = |Re(λ2)| [22], where Re denotes the real part of a complex number.

Figure 8 shows the common logarithm of R1ρ(a) and R2ρ(b) as a function of Δωofffor saturation times of 0.5, 1, 2, 5, and 10 s in the two-pool chemical exchange model ( Figure 1 ). The peaks at 0 Hz (0 ppm) and 1192.8 Hz (4 ppm) derive from pool A and pool B, respectively. As shown in Figure 8 , R1ρand R2ρwere not affected by the saturation time.

Figure 8.

(a) Common logarithm ofR1ρand (b) R2ρ values as a function of Δωofffor various saturation times in the two-pool chemical exchange model.

Figure 9 shows the common logarithm of R1ρ(a) and R2ρ(b) as a function of Δωofffor ω1 values of 25, 50, 100, 150, and 200 Hz in the two-pool chemical exchange model ( Figure 1 ). As shown in Figure 9 , both parameters became broad with increasing ω1 value.

Figure 9.

(a) Common logarithm ofR1ρand (b) R2ρ values as a function of Δωoff for various ω1 values in the two-pool chemical exchange model.

As described above, R1ρand R2ρcan be obtained from the negative of the largest (least negative) real eigenvalue and the absolute value of the largest real part of the complex eigenvalue of the matrix Ain Eq. (2), respectively. We previously reported that there was good agreement between the R1ρand R2ρvalues thus obtained and those obtained numerically [22]. These results appear to indicate the validity of these procedures.

As shown in Figure 8 , R1ρand R2ρwere not affected by the saturation time, because the matrix Ain Eq. (2) is independent of the saturation time. When ω1 was varied, the influence of ω1 on R1ρand R2ρincreased with increasing ω1 value ( Figure 9 ). Especially, the separation between peaks in the R1ρplots degraded with increasing ω1 value [ Figure 9(a) ]. This also appears to be due to the spillover effect.

Advertisement

3. Spin-locking CEST MRI

3.1. Principle of spin-locking

Longitudinal relaxation time in the rotating frame (T1ρ) has been demonstrated to be effective for probing the slow-motion interactions between motion-restricted water molecules and their local macromolecular environment [23] and provides novel image contrast that is not available from conventional MRI techniques. The imaging of biologic tissue based on T1ρis currently being investigated for various tissues, including articular cartilage, breast, and head and neck [24, 25, 26]. In T1ρ-weighted MRI of tissues, the image is sensitive to molecular processes that occur over a range of frequencies determined by the amplitude of an applied SL pulse [23].

As pointed out by Jin et al. [27], the SL approach is useful for improving the signal-to-noise ratio (SNR) in CEST MRI. Furthermore, Kogan et al. [28] demonstrated that a combination of the CEST and SL approaches is useful for detecting proton exchange in the slow-to intermediate-exchange regimes.

As earlier described, the Bloch-McConnell equations for the two-pool chemical exchange model ( Figure 1 ) in the rotating frame with the same frequency as that of the RF-pulse irradiation is given by Eq. (2) [18, 29]. The solution of Eq. (2) can be given by [18]

Mt=eAωω1ϕtM0.E21

Figure 10 illustrates the image of the pulse sequence with SL. We assume that the SL pulse (frequency: ω, amplitude: ω1, and frequency offset: Ω) is applied on the x-axis ( Figure 2 ). The effective magnetic field (B1eff)anditsangle with respect to thezaxis(θ) are given by B1eff=ω12+Ω2/γand θ = tan−1(ω1/Ω), respectively ( Figure 2 ). To achieve SL, the magnetization is first flipped by the θ-degree RF pulse (frequency: ωand amplitude: ω1θ) to the x-z plane, then spin locked by B1effforaduration oftSL, and then flipped back to the z-axis for imaging ( Figure 10 ). The θ-degree RF pulse for flipping is applied on the –y axis, that is, ϕ =  − π/2, whereas the θ-degree RF pulse for flipping back is applied on the y axis, that is, ϕ = π/2. The θ-degree rotation matrix for flipping [R(θ)] is given by [30].

Rθ=eAωaω1θπ/2tθ,E22

Figure 10.

Diagram of spin-locking pulse sequence.

where ω1θand tθdenote the amplitude and the duration of the θ-degree RF-pulse irradiation, respectively ( Figure 10 ), and ω1θ×tθ=θ. Thus, we obtain the magnetization vector immediately after SL for a duration of tSL[M(tSL)] as.

MtSL=eAωω10tSLRθM0.E23

The θ-degree rotation matrix for flipping back to the z-axis [R(−θ)] is given by.

Rθ=eAωaω1θπ/2tθ,E24

Thus, the magnetization vector after flipping back to the z-axis [M+(tSL)] is given by.

M+tSL=RθMtSL=RθeAωω10tSLRθM0.E25

Note that Ω and θare taken to be 0 and π/2, respectively, for an on-resonance SL sequence, whereas the saturation pulse is applied without flipping the magnetization in the sequence without SL such as the conventional CEST sequence [15]. Therefore, the magnetization vector after the saturation pulse [M(tSAT)] in the conventional CEST MRI is simply expressed as.

MtSAT=eAωω10tSATM0,E26

where tSATdenotes the duration of saturation.

3.2. Calculation of T

T1ρcan be obtained numerically by fitting the z component of magnetization for tSL[M+(tSL) given by Eq. (25)] in pool A [MzatSL] to the following equation [30]:

MzatSL=M0aMzssaetSL/T1ρ+Mzssa,E27

where Mzssadenotes the steady-state z component of magnetization in pool A. In this study, we used the Simplex method [31] to calculate T1ρfrom Eq. (27).

The approximate solution for T1ρhas been derived by Trott and Palmer [29]:

T1ρ1R1cos2θ+R2+Rexsin2θ,E28

where θ = tan−1(ω1/Ω), Rex=PaPbΔω2kex/ωae2ωbe2/ωe2+kex2, ωae=ω12+Δωa2, ωbe=ω12+Δωb2, ωe=ω12+Ω2, Ω=ω¯ω, ω¯=Paωa+Pbωb, Δω = Δωb − Δωa = ωb − ωa, and kex = kab + kba. Paand Pbare the fractional sizes of pool A and pool B, and are given by Pa=M0a/M0a+M0band Pb=M0b/M0a+M0b, respectively. R1 and R2 are the population-averaged relaxation rates, and are given by R1=PaR1a+PbR1band R2=PaR2a+PbR2b, respectively. It should be noted that Ω is the population-averaged offset frequency in this case. Mzssain Eq. (27) is approximated by [27].

MzssaM0aR1cos2θR1ρ.E29

Figure 11 shows an example of the three-dimensional plots of the magnetization vector in pool A in the two-pool chemical exchange model ( Figure 1 ). Figure 11(a) and 11(b) show cases without and with SL, respectively. In these cases, the relaxation time constants were assumed to be T1a=1.5s, T2a=0.06s, T1b=0.77s, and T2b=0.033s[32]. tθin Eq. (22) and (24) was taken as 200 μs [27]. ω1θin Eq. (22) and (24) was calculated from ω1θ=θ/tθ. Unless specifically stated, Δω(= ωb − ωa) and ω1 were assumed to be 2400 and 1000 Hz, respectively. Ω was assumed to be 2000 Hz. Thus, θwas tan−1(ω1/Ω) = tan−1(1000/2000) ≈ 26.6 degrees. kex(= kab + kba) was assumed to be 1500 Hz, and kabwas assumed to be given by kab=M0b/M0akba[18]. M0b/M0awas assumed to be 0.03. As shown in Figure 11(a) , when the SL pulse was not applied, the magnetization vector rotated largely around B1eff. On the other hand, when the SL pulse was applied [ Figure 11(b) ], the magnetization vector moved along B1eff, and the rotation around B1effwas suppressed.

Figure 11.

Three-dimensional plots of the magnetization vector in pool A in the two-pool chemical exchange model. (a) and (b) show cases without and with spin-locking pulse, respectively.

When M0b/M0awas 0.003, there was good agreement between the T1ρvalues calculated from Eq. (27) and Trott and Palmer’s solutions given by Eq. (28) (data not shown). When M0b/M0awas 0.03, some difference was observed between them in the off-resonance case. When M0b/M0awas 0.3, large differences were observed between them in both the on- and off-resonance cases [30].

In this study, we developed a simple and fast method for calculating the magnetization vector in SL CEST MRI, in which a simple matrix equation was derived for solving the time-dependent Bloch-McConnell equations in SL MRI [Eq. (25)] and the θ-degree rotation matrix [Eq. (22)] was introduced for considering the effect of the θ-degree RF pulse for flipping the magnetization. As shown in Figure 11 , the trajectory of the magnetization vector in the sequence with SL could be visualized by calculating M(tSL) using Eq. (23), whereas that in the sequence without SL could be visualized by calculating M(tSAT) using Eq. (26). Although Figure 11 shows the three-dimensional plots observed from one direction, we can observe the trajectory of the magnetization vector from various directions by rotating the plot. If we compared the three-dimensional plots with and without SL ( Figure 11 ), then the effect of SL is well understood. Therefore, our method is helpful for visually understanding the effect of SL. In addition, as our method allows us to simply and quickly calculate the time evolution of the magnetization vector under various study conditions in SL CEST MRI, our method can also be useful for optimizing the study conditions in SL CEST MRI.

As previously described, when M0b/M0awas small, that is, when the population of two pools was highly asymmetric, the T1ρvalues calculated from Eq. (27) agreed with the solutions given by Eq. (28). However, the difference between them increased with increasing M0b/M0a[30]. This finding appears to be due to the fact that Trott and Palmer’s solution [Eq. (28)] was derived by approximating the parameters such as relaxation rates using their population-averaged values, and thus the validity of this approximation decreases with decreasing asymmetry in the populations of the two pools.

Although we treated the two-pool chemical exchange model ( Figure 1 ) for analyzing T1ρor R1ρin SL CEST MRI, recent investigations have shown the importance of improved theoretical approaches for describing multi-site chemical exchange phenomena [33, 34]. Thus, Trott and Palmer [33] have tried to generalize their approach for T1ρor R1ρ[29]. For such purposes, it is necessary to expand the Bloch-McConnell equations to those based on multi-pool chemical exchange models. Our method can be easily expanded to multi-pool chemical exchange models by modifying the matrix Agiven by Eq. (4) [20] as previously described, and it is helpful for testing the validity of newly developed approaches for analyzing multi-site chemical exchange phenomena.

Advertisement

4. Correction of B0 and B1

As previously described, the CEST effect has usually been analyzed using MTRasym [Eq. (19)] or PTR [Eq. (20)]. However, these parameters are susceptible to the B0 inhomogeneity of the static magnetic field. When there exists the B0 inhomogeneity, the spillover effect is no longer symmetric. Furthermore, the B1 inhomogeneity of the RF pulse may also cause spatial variation in labeling efficiency and spillover factor [35]. Apart from the efforts in improving magnetic field inhomogeneities using hardware-based methods, such as parallel transmit technologies [36], post-processing algorithms have been developed for field inhomogeneity correction [37, 38].

Kim et al. [37] showed that direct water saturation imaging allows measurement of the absolute water frequency in each voxel, allowing proper centering of Z-spectra on a voxel-by-voxel basis independent of spatial B0 field variations, and that the B0 inhomogeneity in CEST MRI can be corrected on a voxel-by-voxel basis through the centering of Z-spectra. This method is called “water saturation shift referencing (WASSR)” approach. This method, however, would require acquisition of saturation images at 20–40 frequencies [38]. Since the SNR of CEST MRI is low, multiple acquisitions for each frequency offset of complete Z-spectra would be needed, which is not practical in the clinical setting. Zhou et al. demonstrated that a practical six-offset multi-acquisition method combined with a single reference Z-spectrum to acquire high-SNR CEST MRI can accomplish improved CEST MRI with B0 inhomogeneity correction within an acceptable scanning time [38].

A B1-correction of CEST contrasts is crucial for the evaluation of data obtained in clinical studies at high field strengths with strong B1-inhomogeneities. To correct for the B1 inhomogeneity, a B1 map is acquired for correction of Z-spectra using either a calibration [39] or an interpolation approach [40]. Singh et al. [39] developed an approach for B1 inhomogeneity correction based on acquiring calibration data at a coarsely sampled B1 values in conjunction with the measured B1 maps, whereas Windschuh et al. [40] developed an approach based on Lorentzian line fits.

The comprehensive methods like simultaneous mapping of B0 and B1 fields [35, 41], and model-based correction algorithm, [42] have also been developed to improve the accuracy of MTRasym or PTR.

Advertisement

Acknowledgments

This work was supported, in part, by a Grant-in-Aid for Challenging Exploratory Research (Grant No. 25670532) from the Japan Society for the Promotion of Science.

References

  1. 1. Ward K, Aletras A, Balaban R. A new class of contrast agents for MRI based on proton chemical exchange dependent saturation transfer (CEST). Journal of Magnetic Resonance. 2000;143:79-87
  2. 2. Goffeney N, Bulte JWM, Duyn J, Bryant LH, van Zijl PCM. Sensitive NMR detection of cationic-polymer-based gene delivery systems using saturation transfer via proton exchange. Journal of the American Chemical Society. 2001;123:8628-8629
  3. 3. Aime S, Barge A, Delli Castelli D, Fedeli F, Mortillaro A, Nielsen FU, Terreno E. Paramagnetic lanthanide (III) complexes as pH-sensitive chemical exchange saturation transfer (CEST) contrast agents for MRI applications. Magnetic Resonance in Medicine. 2002;47:639-648
  4. 4. Snoussi K, Bulte JWM, Gueron M, van Zijl PCM. Sensitive CEST agents based on nucleic acid imino proton exchange: Detection of poly(rU) and of a dendrimer-poly(rU) model for nucleic acid delivery and pharmacology. Magnetic Resonance in Medicine. 2003;49: 998-1005
  5. 5. Zhou J, Lal B, Wilson DA, Laterra J, van Zijl PCM. Amide proton transfer (APT) contrast for imaging of brain tumors. Magnetic Resonance in Medicine. 2003;50:1120-1126
  6. 6. Sun PZ, Murata Y, Lu J, Wang X, Lo EH, Sorensen AG. Relaxation-compensated fast multislice amide proton transfer (APT) imaging of acute ischemic stroke. Magnetic Resonance in Medicine. 2008;59:1175-1182
  7. 7. Harris RJ, Cloughesy TF, Liau LM, Prins RM, Antonios JP, Li D, Yong WH, Pope WB, Lai A, Nghiemphu PL, Ellingson BM. pH-weighted molecular imaging of gliomas using amine chemical exchange saturation transfer MRI. Neuro Oncol. 2015;17:1514-1524
  8. 8. Van Zijl PCM, Jones CK, Ren J, Malloy CR, Sherry AD. MRI detection of glycogenin vivoby using chemical exchange saturation transfer imaging (glycoCEST). Proceedings of the National Academy of Sciences of the United States of America. 2007;104:4359-4364
  9. 9. Ling W, Regatte RR, Navon G, Jerschow A. Assessment of glycosaminoglycan concentration in vivo by chemical exchange-dependent saturation transfer (gagCEST). Proceedings of the National Academy of Sciences of the United States of America. 2008;105:2266-2270
  10. 10. Cai K, Haris M, Singh A, Kogan F, Greenberg JH, Hariharan H, Detre JA, Reddy R. Magnetic resonance imaging of glutamate. Nature Medicine. 2012;18:302-306
  11. 11. Jones KM, Pollard AC, Pagel MD. Clinical applications of chemical exchange saturation transfer (CEST) MRI. Journal of Magnetic Resonance Imaging. 2017 (in press)
  12. 12. Walker-Samuel S, Ramasawmy R, Torrealdea F, Rega M, Rajkumar V, Johnson SP, Richardson S, Gonçalves M, Parkes HG, Arstad E, Thomas DL, Pedley RB, Lythgoe MF, Golay X. In vivo imaging of glucose uptake and metabolism in tumors. Nature Medicine. 2013;19:1067-1072
  13. 13. Chen LQ, Howison CM, Jeffery JJ, Robey IF, Kuo PH, Pagel MD. Evaluations of extracellular pH within in vivo tumors using acidoCEST MRI. Magnetic Resonance in Medicine. 2014;72:1408-1417
  14. 14. Maruyama S, Ueda J, Kimura A, Murase K. Development and characterization of novel LipoCEST agents based on thermosensitive liposomes. Magnetic Resonance in Medical Sciences. 2016;15:324-334
  15. 15. Sun PZ. Simultaneous determination of labile proton concentration and exchange rate utilizing optimal RF power: Radio frequency power (RFP) dependence of chemical exchange saturation transfer (CEST) MRI. Journal of Magnetic Resonance. 2010;202:155-161
  16. 16. Sun PZ. Simplified and scalable numerical solution for describing multi-pool chemical exchange saturation transfer (CEST) MRI contrast. Journal of Magnetic Resonance. 2010;205:235-241
  17. 17. Woessner DE, Zhang S, Merritt ME, Sherry AD. Numerical solution of the Bloch equations provides insights into the optimum design of PARACEST agents for MRI. Magnetic Resonance in Medicine. 2005;53:790-799
  18. 18. Murase K, Tanki N. Numerical solutions to the time-dependent Bloch equations revisited. Magnetic Resonance Imaging. 2011;29:126-131
  19. 19. Koss H, Rance M, Palmer AG. General expression forR1ρ relaxation forN-site chemical exchange and the special case of linear chains. Journal of Magnetic Resonance. 2017;274:36-45
  20. 20. Murase K. Numerical analysis of the magnetization behavior in magnetic resonance imaging in the presence of multiple chemical exchange pools. Open Journal of Applied Sciences. 2017;7:1-14
  21. 21. Mulkern RV, Williams ML. The general solution to the Bloch equation with constant rf and relaxation terms: Application to saturation and slice selection. Medical Physics. 1993;20:5-13
  22. 22. Murase K. A theoretical and numerical consideration of the longitudinal and transverse relaxations in the rotating frame. Magnetic Resonance Imaging. 2013;31:1544-1558
  23. 23. Li X, Han ET, Busse RF, Majumdar S.In vivo T1ρ mapping in cartilage using 3D magnetization prepared angle-modulated partitionedk-space spoiled gradient echo snapshots (3D MAPSS). Magnetic Resonance in Medicine. 2008;59:298-307
  24. 24. Regatte RR, Akella SV, Borthakur A, Kneeland JB, Reddy R.In vivoproton MR three-dimensionalT1ρ mapping of human articular cartilage: Initial experience. Radiology. 2003;229:269-274
  25. 25. Santyr GE, Henkelman RM, Bronskill MJ. Spin locking for magnetic resonance imaging with application to human breast. Magnetic Resonance in Medicine. 1989;12:25-37
  26. 26. Markkola AT, Aronen HJ, Paavonen T, Hopsu E, Sipila LM, Tanttu JI, Sepponen RE. Spin lock and magnetization transfer imaging of head and neck tumors. Radiology. 1996;200:369-375
  27. 27. Jin T, Autio J, Obata T, Kim S-G. Spin-locking versus chemical exchange saturation transfer MRI for investigating chemical exchange process between water and labile metabolite protons. Magnetic Resonance in Medicine. 2011;65:1448-1460
  28. 28. Kogan F, Singh A, Cai K, Haris M, Hariharan H, Reddy R. Investigation of chemical exchange at intermediate exchange rates using a combination of chemical exchange saturation transfer (CEST) and spin-locking methods (CESTrho). Magnetic Resonance in Medicine. 2012;68:107-119
  29. 29. Trott O, Palmer AG.R1ρrelaxation outside of the fast-exchange limit. Journal of Magnetic Resonance. 2002;154:157-160
  30. 30. Murase K. Behavior of the magnetization in spin-locking magnetic resonance imaging using numerical solutions to the time-dependent Bloch equations. Physics in Medicine and Biology. 2012;57:N481-N492
  31. 31. Nelder JA, Mead R. A simplex method for function minimization. The Computer Journal. 1965;7:308-313
  32. 32. Sun PZ, Zhou J, Huang J, van Zijl P. Simplified quantitative description of amide proton transfer (APT) imaging during acute ischemia. Magnetic Resonance in Medicine. 2007;57:405-410
  33. 33. Trott O, Palmer AG. Theoretical study ofR rotating-frame andR2 free-precession relaxation in the presence of n-site chemical exchange. Journal of Magnetic Resonance. 2004;170:104-112
  34. 34. Grey MJ, Wang C, Palmer AG. Disulfide bond isomerization in basic pancreatic trypsin inhibitor: Multisite chemical exchange quantified by CPMG relaxation dispersion and chemical shift modeling. Journal of the American Chemical Society. 2003;125:14324-14335
  35. 35. Ji Y, Zhou IY, Qiu B, Sun PZ. Progress toward quantitative in vivo chemical exchange saturation transfer (CEST) MRI. Israel Journal of Chemistry. 2017;57:1-17
  36. 36. Padormoa F, Beqiria A, Hajnalab JV, Malika SJ. Parallel transmission for ultrahigh-field imaging. NMR in Biomedicine. 2016;29:1145-1161
  37. 37. Kim M, Gillen J, Landman BA, Zhou J, Van Zijl PC. Water saturation shift referencing (WASSR) for chemical exchange saturation transfer (CEST) experiments. Magnetic Resonance in Medicine. 2009;61:1441-1450
  38. 38. Zhou J, Blakeley JO, Hua J, Kim M, Laterra J, Pomper MG, Van Zijl PC. Practical data acquisition method for human brain tumor amide proton transfer (APT) imaging. Magnetic Resonance in Medicine. 2008;60:842-849
  39. 39. Singh A, Cai K, Haris M, Hariharan H, Reddy R. OnB1 Inhomogeneity correction ofin vivohuman brain glutamate chemical exchange saturation transfer contrast at 7T. Magnetic Resonance in Medicine. 2013;69:818-824
  40. 40. Windschuh J, Zaiss M, Meissner J-E, Paech D, Radbruch A, Ladd ME, Bachert P. Correction of B1-inhomogeneities for relaxation-compensated CEST imaging at 7T. NMR in Biomedicine. 2015;28:529-537
  41. 41. Schuenke P, Windschuh J, Roeloffs V, Ladd ME, Bachert P, Zaiss M. Simultaneous mapping of water shift and B1 (WASABI)-application to field-inhomogeneity correction of CEST MRI data. Magnetic Resonance in Medicine. 2017;77:571-580
  42. 42. Sun PZ, Farrar CT, Sorensen AG. Correction for artifacts induced by B0 and B1 field inhomogeneities in pH-sensitive chemical exchange saturation transfer (CEST) imaging. Magnetic Resonance in Medicine. 2007;58:1207-1215

Written By

Kenya Murase

Submitted: June 10th, 2017 Reviewed: October 12th, 2017 Published: December 20th, 2017