Abstract
Chemical exchange saturation transfer (CEST) MRI has been increasingly applied to detect dilute solutes and physicochemical properties, with promising in vivo applications. Whereas CEST imaging has been implemented with continuous wave (CW) radio-frequency irradiation on preclinical scanners, pulse-train irradiation is often chosen on clinical systems. Therefore, it is necessary to optimize pulse-train CEST imaging, particularly important for translational studies. Because conventional Bloch-McConnell formulas are not in the form of homogeneous differential equations, the routine simulation approach simulates the evolving magnetization step by step, which is time consuming. Herein we developed a computationally efficient numerical solution using matrix iterative analysis of homogeneous Bloch-McConnell equations. The proposed algorithm requires simulation of pulse-train CEST MRI magnetization within one irradiation repeat, with 99% computation time reduction from that of conventional approach under typical experimental conditions. The proposed solution enables determination of labile proton ratio and exchange rate from pulse-train CEST MRI experiment, within 5% from those determined from quantitative CW-CEST MRI. In addition, the structural similarity index analysis shows that the dependence of CEST contrast on saturation pulse flip angle and duration between simulation and experiment was 0.98 ± 0.01, indicating that the proposed simulation algorithm permits fast optimization and quantification of pulse-train CEST MRI.