bSSFP was shown to be sensitive to exchange and is explored as an alternative way for CEST/T1ρ 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 T1ρ 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 T2ρ 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=R→e(Φ)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 R→e(Φ), 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 R→e(Φ).
When relaxation is included, Eq.1 becomes: Mn+1=e(Heff−˜R)TRMn−[I−e(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ΛeffD−1 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ΛRD−1 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Λ−1D−1R1M0 (5)where Λ=Λeff−ΛR.
Furthermore, the time-dependent magnetization can be derived from Eqs.3 and 5: Mn=Mss+DeΛ⋅nTRD−1(M0−Mss) (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.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 T2ρ and T1ρ 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.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.