Processing math: 100%

0994

Quantitative modeling of exchange in bSSFPX
Shu Zhang1, Robert E Lenkinski1,2, and Elena Vinogradov1,2

1Department of Radiology, UT Southwestern Medical Center, Dallas, TX, United States, 2Advanced Imaging Research Center, UT Southwestern Medical Center, Dallas, TX, United States

Synopsis

bSSFP was shown to be sensitive to exchange and is explored as an alternative way for CEST/T experiments (the bSSFPX method). In this abstract, an analytical solution is derived for the magnetization behavior in the bSSFPX. The solution describes the transient signal fluctuations for short saturation times. The solution is verified by comparing it to full, step-wise computation of Bloch-McConnell Equations. Overall, this solution is in good agreement with Bloch-McConnell Equations and follows the transient signal oscillations well. Work is in progress to use this solution to quantify exchange rates experimentally.

Purpose

Balanced steady state free precession (bSSFP) was shown to be sensitive to exchange and is explored as an alternative way for chemical exchange saturation transfer (CEST) and T experiments (the bSSFPX method 1,2). In this method, the data are acquired during saturation, generating exchange-dependent images at different saturation times in a “single-shot”. 2 This allows a QUEST-like quantification. 3 In this abstract, an analytical solution is derived for the magnetization behavior in the bSSFPX. The solution takes T into consideration (often ignored in the existing models for exchange 4). Therefore, it can describe the transient signal fluctuations for short saturation times or agents with long T2’s. The solution is verified by comparing it to full, step-wise computation of Bloch-McConnell Equations (BME). It lays the foundation for exchange rate quantification using bSSFPX. Although this solution is primarily derived for bSSFPX, it can be expanded in a straightforward manner to the modeling of pulsed CEST experiments.

Theory

Consider only 1 pool first. In a bSSFP sequence, 5 if the relaxation is ignored, Mn+1=Rz(θ)Rx(α)Mn=Re(Φ)Mn=eHeffTRMn    (1)where Mn+1 and Mn are the magnetization vectors at the end of n+1 and n cycle. Rz(θ) and Rx(α) are the standard z- and x-rotation matrices respectively, where θ=2πΔTR is the precession angle during TR and Δ is the off-resonance shift. According to the Euler theorem, Rz(θ)Rx(α) can be represented by a single rotation Re(Φ), where e is the directionality vector and Φ is the rotation angle about e. The effective field (Fig.1) is in the same direction as e and with the magnitude ωbSSFPeff: ωbSSFPeff=Φ/TR    (2)The eHeffTR is the matrix exponential form of Re(Φ).

When relaxation is included, Eq.1 becomes: Mn+1=e(Heff˜R)TRMn[Ie(Heff˜R)TR](Heff˜R)1R1M0   (3)where ˜R=diag(R2,R2,R1) is the standard relaxation matrix and M0=[0,0,M0]T. From Eq.3, the steady-state magnetization Mss can be obtained: Mss=(Heff˜R)1R1M0   (4) Heff can be expressed in its eigenbasis as Heff=DΛeffD1 where D is the diagonalization matrix and Λeff is the eigenvalue diagonal matrix. ˜R can also be diagonalized in the eigenbasis of Heff, which yields ˜R=DΛRD1 where ΛR=diag(R1ρ,R2ρ,R2ρ), R1ρ=cos2(Θ)R1+sin2(Θ)R2, R2ρ=sin2(Θ)R1+cos2(Θ)R2, and Θ is the angle between the effective field and the Z-axis (Fig.1). Using the eigensystem of Heff, the Mss becomes: Mss=DΛ1D1R1M0   (5)where Λ=ΛeffΛR.

Furthermore, the time-dependent magnetization can be derived from Eqs.3 and 5: Mn=Mss+DeΛnTRD1(M0Mss)   (6)

Now consider 2 pools, water and solute, coupled by exchange. When exchange exists, an additional term Rex is introduced to the water via R1ρ and R2ρ: R1ρ,w=cos2(Θw)R1w+sin2(Θw)(R2w+Rex) and R2ρ,w=sin2(Θw)R1w+cos2(Θw)(R2w+Rex). Here we follow Trott and Palmer’s work 6 for Rex (asymmetric population limit): Rex=pwpsΔ2wsk(wbSSFPeff,s/2π)2+k2=pwpsΔ2wskΦ2s/(2πTR)2+k2   (7)where pw=M0w/(M0w+M0s), ps=M0s/(M0w+M0s), Δws is the chemical shift between water and solute, and k=ksw+kws, ksw and kws are exchange rates from solute to the water pool and vice-versa.

Methods

Simulations were performed using two models: a 1-pool model with M0w=1, T1w/T2w=2/0.05s and a 2-pool model with M0w=1, T1w/T2w=2/0.05s, M0s=0.005, T1s/T2s=1/1s, ∆ws=1ppm, ksw=100Hz. For all the simulations, rectangular pulses with duration 25μs, α=30o and TR=2ms were used. The magnetizations were calculated for each TR for total saturation time tsat=10s (5 times of T1w or 5T1w). To show the time evolution of the Z- and XY-spectrum, the spectra with tsat=3T2w, 5T2w, 0.25T1w, 0.5T1w and 3T1w were plotted.

Results and Discussion

Fig.2 compares the steady-state spectra obtained using the analytical solution with the exact BME simulations. Fig.3 shows the time evolution of the Z- and XY-spectra. The analytical solution is in excellent agreement with BME for both Z- and XY-spectrum for the 1-pool model.

For the exchanging 2-pools, the solution is in good qualitative agreement with BME. At the steady-state, the solution follows BME for the most part but deviates slightly downward at the star labeled frequencies close to water (Fig.2b,d). In Fig.3, the overall behavior and the agreement at both short and long tsat is excellent, indicating that the solution correctly depicts transient oscillations with T and T relaxations. However, large deviations were observed close to the saturation bands for short tsat (tsat<3T1w). These results indicate that the exchange contribution (Eq.7) needs to be further improved.

Conclusion

In the abstract, an analytical solution to the BME is derived for bSSFPX based on the effective field. Overall, this solution is in good agreement with BME and follows the transient signal oscillations well. Deviations observed in the regions close to saturation bands at short tsat and at star labelled frequencies (Fig.2b,d *) in the 2-pool model indicate that the exchange contribution needs to be further improved. Work is in progress to use this solution to quantify the exchange rate experimentally.

Acknowledgements

No acknowledgement found.

References

1. Zhang S, Zheng L, Lenkinski RE, et al. in: Proc Intl Soc Mag Reson Med. 2015;23:0785.

2. Zhang S, Lenkinski RE, Vinogradov E. in: Proc Intl Soc Mag Reson Med. 2016;0304.

3. McMahon MT, Gilad AA, Zhou J, et al. Quantifying exchange rates in chemical exchange saturation transfer agents using the saturation time and saturation power dependencies of the magnetization transfer effect on the magnetic resonance imaging signal (QUEST and QUESP): pH calibration for poly-L-lysine and a starburst dendrimer. Magn Reson Med. 2006;55(4):836-847.

4. Zaiss M, Bachert P. Chemical exchange saturation transfer (CEST) and MR Z-spectroscopy in vivo: a review of theoretical approaches and methods. Phys Med Biol. 2013;58:R221-R269.

5. Freeman R, Hill H. Phase and intensity anomalies in Fourier transform NMR. J Magn Reson (1969). 1971;4(3):366-383.

6. Trott O, Palmer III AG. R1ρ Relaxation outside of the Fast-Exchange Limit. J Magn Reson. 2002;154:157-160.

Figures

Figure 1. The effective RF field of bSSFP. Θ is the angle between the effective field and the Z-axis. θ=2πΔTR is the precession angle during TR, where Δ is the off-resonance shift.

Figure 2. The simulated steady-state Z- (a,b) and XY-spectra (c,d) for 1-pool model (a,c) and 2-pool model (b,d). Circle: Bloch-McConnell equations. Solid line: the derived analytical solution. *: solution deviates slightly downward from the Bloch-McConnell equations.

Figure 3. The simulated time-dependent Z- (a,b) and XY-spectra (c,d) for 1-pool model (a,c) and 2-pool model (b,d). Circle: Bloch-McConnell equations. Solid line: the derived analytical solution. Different colors from light to dark correspond to different saturation times: 0.15, 0.25, 0.5, 1, and 6s.

Proc. Intl. Soc. Mag. Reson. Med. 25 (2017)
0994