Towards an analytic solution for pulsed CEST

Rapid and quantitative CEST-MRI sequence using water presaturation

PURPOSE

Despite the significant potential for in vivo metabolic imaging in preclinical and clinical applications, CEST MRI suffers from long scan time and inaccurate quantification. This study aims to suppress the contaminations among signals under different frequencies, which could shorten the TR and thereby facilitate CEST imaging acceleration and quantification.

METHODS

A novel sequence is proposed by applying a water-presaturation (WPS) module at the beginning of each TR. WPS CEST quickly knocks down the residual signal from previous TRs so that the magnetization of all TRs recovers from zero, which aligns well with the formula of quasi-steady-state theorem and enables accurate quantification within shorter TR. WPS CEST was assessed by simulations, creatine phantom, and healthy human brain scans at 3 T.

RESULTS

In simulation and phantom experiment, WPS CEST allows accurate estimation of exchange rate (ksw) using omega plot and using shorter delay time (Td) and saturation time (Ts) (e.g., 1 s/1 s) compared with the conventional CEST. Simulations further showed that WPS CEST could obtain consistent spin-lock relaxation (R1ρ) values over varied Tds and Tss. Six human scans indicated that R1ρ collected from conventional sequences showed significant differences between two groups with Td and Ts of (1 s/1 s) and (2 s/2 s) (amide: 1.721 ± 0.051 s-1 vs. 1.622 ± 0.050 s-1, p = 0.001; nuclear Overhauser enhancement: 1.792 ± 0.046 s-1 vs. 1.687 ± 0.053 s-1, p = 0.004), whereas WPS CEST scans using these 2 Td/Ts values obtained the same mean R1ρ (amide: 1.616 ± 0.053 s-1 vs. 1.616 ± 0.048 s-1, p = 0.862; nuclear Overhauser enhancement: 1.688 ± 0.064 s-1 vs. 1.684 ± 0.054 s-1, p = 0.544).

CONCLUSION

WPS CEST demonstrated accurate quantitation within shorter TR compared with conventional sequences, and thereby may allow rapid quantitative CEST scans in various situations.

Adjustment of rotation and saturation effects (AROSE) for CEST imaging

PURPOSE

Endogenous CEST signal usually has low specificity due to contaminations from the magnetization transfer contrast (MTC) and other labile protons with overlapping or close Larmor frequencies. We propose to improve CEST signal specificity with adjustment of rotation and saturation effects (AROSE).

METHODS

The AROSE approach measures the difference between CEST signals acquired with the same average irradiation power but largely different duty cycles, for example, a continuous wave or a high duty cycle pulse train versus a low duty cycle pulse train with a flip angle φ. Simulation, phantom, and in vivo rodent studies were performed to evaluate the characteristics of the AROSEφ signal.

RESULTS

Simulation and experimental results show that AROSE2π is a low-pass filter that can suppress fast exchanging processes (e.g., >3000 s-1 ), whereas AROSEπ is a band-pass filter suppressing both fast and slow exchange (e.g., <30 s-1 ) rates. For other φ angles, the sensitivity and the exchange-rate filtering effect of AROSEφ falls between AROSEπ and AROSE2π . AROSE can also minimize MTC and improve the Larmor frequency selectivity of the CEST signal. The linewidth of the AROSE1.5π spectrum is about 60% to 65% when compared to the CEST spectrum measured by continuous wave. Depending on the needs of an application, the sensitivity, exchange-rate filtering, and Larmor frequency selectivity can be adjusted by varying the flip angle, duty cycle, and average irradiation power.

CONCLUSION

Compared to conventional CEST signals, AROSE can minimize MTC and improve exchange rate filtering and Larmor frequency specificity.

Individually tailored spatial-spectral pulsed CEST MRI for ratiometric mapping of myocardial energetic species at 3T

PURPOSE

CEST MRI has been used to probe changes in cardiac metabolism via assessment of CEST contrast from Cr. However, B1 variation across the myocardium leads to spatially variable Cr CEST contrast in healthy myocardium.

METHODS

We developed a spatial-spectral (SPSP) saturation pulsed CEST protocol to compensate for B1 variation. Flip angle maps were used to individually tailor SPSP pulses comprised of a train of one-dimensional spatially selective subpulses selective along the principal B1 gradient dimension. Complete Z-spectra in the hearts of (n = 10) healthy individuals were acquired using conventional Gaussian saturation and SPSP schemes and supported by phantom studies.

RESULTS

In simulations, the use of SPSP pulses reduced the average SD of the effective saturation B1 values within the myocardium (n = 10) from 0.12 ± 0.02 μT to 0.05 ± 0.01 μT (p < 0.01) and reduced the average SD of Cr CEST contrast in vivo from 10.0 ± 4.3% to 6.1 ± 3.5% (p < 0.05). Results from the hearts of human subjects showed a significant reduction of CEST contrast distribution at 2 ppm, as well as amplitude, when using SPSP saturation. Corresponding phantom experiments revealed PCr-specific contrast generation at body temperature when SPSP saturation was used but combined PCr and Cr contrast generation when Gaussian saturation was used.

CONCLUSION

The use of SPSP saturation pulsed CEST resulted in PCr-specific contrast generation and enabled ratiometric mapping of PCr to total Cr CEST contrast in the human heart at 3T.

Numerical simulation-based assessment of pH-sensitive chemical exchange saturation transfer MRI quantification accuracy across field strengths

Chemical exchange saturation transfer (CEST) MRI detects dilute labile protons via their exchange with bulk water, conferring pH sensitivity. Based on published exchange and relaxation properties, a 19-pool simulation was used to model the brain pH-dependent CEST effect and assess the accuracy of quantitative CEST (qCEST) analysis across magnetic field strengths under typical scan conditions. First, the optimal B1 amplitude was determined by maximizing pH-sensitive amide proton transfer (APT) contrast under the equilibrium condition. Apparent and quasi-steady-state (QUASS) CEST effects were then derived under the optimal B1 amplitude as functions of pH, RF saturation duration, relaxation delay, Ernst flip angle, and field strength. Finally, CEST effects, particularly the APT signal, were isolated with spinlock model-based Z-spectral fitting to evaluate the accuracy and consistency of CEST quantification. Our data showed that QUASS reconstruction significantly improved the consistency between simulated and equilibrium Z-spectra. The residual difference between QUASS and equilibrium CEST Z-spectra was, on average, 30 times less than that of the apparent CEST Z-spectra across field strengths, saturation, and repetition times. Also, the spinlock fitting of the QUASS CEST effect significantly reduced the residual errors 9-fold. Furthermore, the isolated APT amplitude from QUASS reconstruction was consistent and higher than the apparent CEST analysis under nonequilibrium conditions. To summarize, this study confirmed that QUASS reconstruction facilitates accurate determination of the CEST system under different scan protocols across field strengths, with the potential to help standardize CEST quantification.

Detection of tissue pH with quantitative chemical exchange saturation transfer magnetic resonance imaging

Chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) has emerged as a novel means for sensitive detection of dilute labile protons and chemical exchange rates. By sensitizing to pH-dependent chemical exchange, CEST MRI has shown promising results in monitoring tissue statuses such as pH changes in disorders like acute stroke, tumor, and acute kidney injury. This article briefly reviews the basic principles for CEST imaging and quantitative measures, from the simplistic asymmetry analysis to multipool Lorentzian decoupling and quasi-steady-state reconstruction. In particular, the advantages and limitations of commonly used quantitative approaches for CEST applications are discussed.

Theory of chemical exchange saturation transfer MRI in the context of different magnetic fields

Chemical exchange saturation transfer (CEST) magnetic resonance imaging (MRI) is a versatile MRI method that provides contrast based on the level of molecular and metabolic activity. This contrast arises from indirect measurement of protons in low concentration molecules that are exchanging with the abundant water proton pool. The indirect measurement is based on magnetization transfer of radio frequency (rf)-prepared magnetization from the small pool to the water pool. The signal can be modeled by the Bloch-McConnell equations combining standard magnetization dynamics and chemical exchange processes. In this article, we review analytical solutions of the Bloch-McConnell equations and especially the derived CEST signal equations and their implications. The analytical solutions give direct insight into the dependency of measurable CEST effects on underlying parameters such as the exchange rate and concentration of the solute pools, but also on the system parameters such as the rf irradiation field B1 , as well as the static magnetic field B0 . These theoretical field-strength dependencies and their influence on sequence design are highlighted herein. In vivo results of different groups making use of these field-strength benefits/dependencies are reviewed and discussed.

Average saturation efficiency filter (ASEF) for CEST imaging

PURPOSE

Endogenous CEST signal usually has low specificity due to contamination from the magnetization transfer effect and from fast exchanging labile protons with close Larmor frequencies. We propose to improve CEST signal specificity with an average saturation efficiency filter (ASEF).

METHODS

ASEF measures the difference between CEST signals acquired with similar average irradiation power but largely different duty cycles (DC), for example, a continuous wave or a high DC pulse train versus a low DC one. Simulation and Cr phantom studies were performed to evaluate the characteristics of ASEF for CEST.

RESULTS

Theoretical and simulation studies show that ASEF can suppress fast exchanging processes, with only a small loss of chemical exchange contrast for slow-to-intermediate exchange rates if the difference in DC is large. In the RF offset range of 2 to 5 ppm with an averaged saturation power of 0.8 and 1.6 microteslas, there is a mismatch of ∼0.1% to 2% in the magnetization transfer signal between saturation by continuous wave and a pulse train with DC = 15% and pulse duration of 24 ms, respectively. This mismatch can be minimized by careful selection of saturation power, pulse duration, and DC differences or by applying a small fudge factor between the 2 irradiation powers. Phantom studies of Cr confirmed that ASEF can minimize the magnetization transfer effect and reduce sensitivity to fast exchange processes.

CONCLUSION

ASEF can improve the specificity of slow-to-intermediate exchanging CEST signal with a relatively small loss of sensitivity.

Mapping pH using stimulated echoes formed via chemical exchange

PURPOSE

RACETE (refocused acquisition of chemical exchange transferred excitations) is a recently developed approach to imaging solute exchange with water. However, it lacks biophysical specificity, as it is sensitive to exchange rates, relaxation rates, solute concentration, and macromolecular content. We modified this sequence and developed a protocol and corresponding metric with specific sensitivity to the solute exchange rate and hence a means for mapping pH.

THEORY AND METHODS

RACETE splits the two gradients traditionally used in a stimulated-echo sequence into one applied after exciting solutes and one applied after exciting water, hence requiring exchange for echo formation. In this work, we leverage the dependence of the stimulated-echo signal on the exchange process. By preserving the total irradiation power and using a ratio metric, the other signal dependencies cancel, leaving a specific measure of exchange rate. Additionally, artifacts due to off-resonance excitation of water are addressed using a phase cancelling approach; and a gradient-echo imaging sequence with a variable flip angle excitation is tailored for a fast read-out of RECETE prepared signals. This method is validated using numerical simulations and salicylic acid phantom experiments at 9.4 T.

RESULTS

Numerical simulations and phantom experiments demonstrate that the ratio-metric is a single-variable function of exchange rate with extremely low dependence on confounding factors. Additionally, artifacts due to direct water excitation are removed and robustness to B0 and B1 inhomogeneities is demonstrated.

CONCLUSION

The proposed method can be used for fast pH mapping with robustness against the confounding effects that widely exist in other methods.

An end-to-end AI-based framework for automated discovery of rapid CEST/MT MRI acquisition protocols and molecular parameter quantification (AutoCEST)

PURPOSE

To develop an automated machine-learning-based method for the discovery of rapid and quantitative chemical exchange saturation transfer (CEST) MR fingerprinting acquisition and reconstruction protocols.

METHODS

An MR physics-governed AI system was trained to generate optimized acquisition schedules and the corresponding quantitative reconstruction neural network. The system (termed AutoCEST) is composed of a CEST saturation block, a spin dynamics module, and a deep reconstruction network, all differentiable and jointly connected. The method was validated using a variety of chemical exchange phantoms and in vivo mouse brains at 9.4T.

RESULTS

The acquisition times for AutoCEST optimized schedules ranged from 35 to 71 s, with a quantitative image reconstruction time of only 29 ms. The resulting exchangeable proton concentration maps for the phantoms were in good agreement with the known solute concentrations for AutoCEST sequences (mean absolute error = 2.42 mM; Pearson's r=0.992 , p<0.0001 ), but not for an unoptimized sequence (mean absolute error = 65.19 mM; Pearson's r=-0.161 , p=0.522 ). Similarly, improved exchange rate agreement was observed between AutoCEST and quantification of exchange using saturation power (QUESP) methods (mean absolute error: 35.8 Hz, Pearson's r=0.971 , p<0.0001 ) compared to an unoptimized schedule and QUESP (mean absolute error = 58.2 Hz; Pearson's r=0.959 , p<0.0001 ). The AutoCEST in vivo mouse brain semi-solid proton volume fractions were lower in the cortex (12.77% ± 0.75%) compared to the white matter (19.80% ± 0.50%), as expected.

CONCLUSION

AutoCEST can automatically generate optimized CEST/MT acquisition protocols that can be rapidly reconstructed into quantitative exchange parameter maps.

An end-to-end AI-based framework for automated discovery of rapid CEST/MT MRI acquisition protocols and molecular parameter quantification (AutoCEST)

PURPOSE

To develop an automated machine-learning-based method for the discovery of rapid and quantitative chemical exchange saturation transfer (CEST) MR fingerprinting acquisition and reconstruction protocols.

METHODS

An MR physics-governed AI system was trained to generate optimized acquisition schedules and the corresponding quantitative reconstruction neural network. The system (termed AutoCEST) is composed of a CEST saturation block, a spin dynamics module, and a deep reconstruction network, all differentiable and jointly connected. The method was validated using a variety of chemical exchange phantoms and in vivo mouse brains at 9.4T.

RESULTS

The acquisition times for AutoCEST optimized schedules ranged from 35 to 71 s, with a quantitative image reconstruction time of only 29 ms. The resulting exchangeable proton concentration maps for the phantoms were in good agreement with the known solute concentrations for AutoCEST sequences (mean absolute error = 2.42 mM; Pearson's r = 0.992 , p < 0.0001 ), but not for an unoptimized sequence (mean absolute error = 65.19 mM; Pearson's r = - 0.161 , p = 0.522 ). Similarly, improved exchange rate agreement was observed between AutoCEST and quantification of exchange using saturation power (QUESP) methods (mean absolute error: 35.8 Hz, Pearson's r = 0.971 , p < 0.0001 ) compared to an unoptimized schedule and QUESP (mean absolute error = 58.2 Hz; Pearson's r = 0.959 , p < 0.0001 ). The AutoCEST in vivo mouse brain semi-solid proton volume fractions were lower in the cortex (12.77% ± 0.75%) compared to the white matter (19.80% ± 0.50%), as expected.

CONCLUSION

AutoCEST can automatically generate optimized CEST/MT acquisition protocols that can be rapidly reconstructed into quantitative exchange parameter maps.

A hybrid numeric-analytic solution for pulsed CEST

Chemical exchange saturation transfer (CEST) methods measure the effect of magnetization exchange between solutes and water. While CEST methods are often implemented using a train of off-resonant shaped RF pulses, they are typically analyzed as if the irradiation were continuous. This approximation does not account for exchange of rotated magnetization, unique to pulsed irradiation and exploited by chemical exchange rotation transfer methods. In this work, we derive and test an analytic solution for the steady-state water signal under pulsed irradiation by extending a previous work to include the effects of pulse shape. The solution is largely accurate at all offsets, but this accuracy diminishes at higher exchange rates and when applying pulse shapes with large root-mean-squared to mean ratios (such as multi-lobe sinc pulses).

Quasi-steady-state CEST (QUASS CEST) solution improves the accuracy of CEST quantification: QUASS CEST MRI-based omega plot analysis

PURPOSE

CEST MRI omega plot quantifies the labile proton fraction ratio (f_r ) and exchange rate (k_sw ), yet it assumes long RF saturation time (Ts) and relaxation delay (Td). Our study aimed to test if a quasi-steady-state (QUASS) CEST analysis that accounts for the effect of finite Ts and Td could improve the accuracy of CEST MRI quantification.

METHODS

We modeled the MRI signal evolution using a typical CEST EPI sequence. The signal relaxes toward its thermal equilibrium following the bulk water relaxation rate during Td, and then toward its CEST steady state following the spin-lock relaxation rate during Ts from which the QUASS CEST effect is derived. Both f_r and k_sw were solved from simulated conventional apparent CEST and QUASS CEST MRI. We also performed MRI experiments from a Cr-gel phantom under serially varied Ts and Td times from 1.5 to 7.5 s.

RESULTS

Simulation showed that, although k_sw could be slightly overestimated (3%-15%) for the range of Ts and Td, f_r could be substantially underestimated by as much as 67%. In contrast, the QUASS solution provided accurate k_sw and f_r determination within 2%. The CEST MRI experiments confirmed that the QUASS solution enabled robust quantification of k_sw and f_r , superior over the omega plot analysis based on the conventional apparent CEST MRI measurements.

CONCLUSIONS

The QUASS CEST MRI algorithm corrects the effect of finite Ts and Td times on CEST measurements, thereby allowing robust and accurate CEST quantification.

Pulse sequences for measuring exchange rates between proton species: From unlocalised NMR spectroscopy to chemical exchange saturation transfer imaging

Within the field of NMR spectroscopy, the study of chemical exchange processes through saturation transfer techniques has a long history. In the context of MRI, chemical exchange techniques have been adapted to increase the sensitivity of imaging to small fractions of exchangeable protons, including the labile protons of amines, amides and hydroxyls. The MR contrast is generated by frequency-selective irradiation of the labile protons, which results in a reduction of the water signal associated with transfer of the labile protons' saturated magnetization to the protons of the surrounding free water. The signal intensity depends on the rate of chemical exchange and the concentration of labile protons as well as on the properties of the irradiation field. This methodology is referred to as CEST (chemical exchange saturation transfer) imaging. Applications of CEST include imaging of molecules with short transverse relaxation times and mapping of physiological parameters such as pH, temperature, buffer concentration and chemical composition due to the dependency of this chemical exchange effect on all these parameters. This article aims to describe these effects both theoretically and experimentally. In depth analysis and mathematical modelling are provided for all pulse sequences designed to date to measure the chemical exchange rate. Importantly, it has become clear that the background signal from semi-solid protons and the presence of the Nuclear Overhauser Effect (NOE), either through direct dipole-dipole mechanisms or through exchange-relayed signals, complicates the analysis of CEST effects. Therefore, advanced methods to suppress these confounding factors have been developed, and these are also reviewed. Finally, the experimental work conducted both in vitro and in vivo is discussed and the progress of CEST imaging towards clinical practice is presented.

Assessment of a clinically feasible Bayesian fitting algorithm using a simplified description of Chemical Exchange Saturation Transfer (CEST) imaging

Fitting a model based on the Bloch-McConnell (BM) equations to Chemical Exchange Saturation Transfer (CEST) spectra allows for the quantification of metabolite concentration and exchange rate as well as simultaneous correction of field inhomogeneity, direct water saturation and magnetization transfer. Employing a Bayesian fitting approach permits the integration of prior information into the analysis to incorporate expected parameter distributions and to prevent over-fitting. However, the analysis can be time consuming if a general numerical solution of the BM equations is applied. In this study, we combined a Bayesian fitting algorithm with approximate analytical solutions of the BM equations to achieve feasible computational times. To evaluate the accuracy and speed of the suggested approach, phantoms including Iodipamide, Taurine and Creatine were tested in addition to simulated data with continuous-wave (CW) and pulsed saturation with Gaussian pulses. A significant reduction of computational time was achieved when fitting CW data (about 50-fold) and pulsed saturation data (more than 100-fold) with the analytical model while the estimated parameters were largely consistent with the parameters from the general numerical solution. The increased speed of the algorithm facilitates the Bayesian analysis of CEST data within clinically feasible processing times. Other analytical models valid for different parameter regimes may be employed to extend the applicability to a wider range of CEST agents.