Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Magn Reson Med. Author manuscript; available in PMC 2012 October 1.
Published in final edited form as:
PMCID: PMC3135736

Simulation and optimization of pulsed radio frequency (RF) irradiation scheme for chemical exchange saturation transfer (CEST) MRI – demonstration of pH-weighted pulsed-amide proton CEST MRI in an animal model of acute cerebral ischemia


Chemical exchange saturation transfer (CEST) MRI is capable of measuring dilute labile protons and microenvironment properties; however, the CEST contrast is also dependent upon experimental conditions, particularly, the RF irradiation scheme. Although continuous-wave (CW) RF irradiation has been conventionally utilized, the RF pulse duration or duty cycle are limited on most clinical systems, for which pulsed RF irradiation must be chosen. Here, conventional numerical simulation was extended to describe pulsed-CEST MRI contrast as a function of RF pulse parameters (i.e., RF pulse duration and flip angle) and labile proton properties (i.e., exchange rate and chemical shift). For diamagnetic CEST agents undergoing slow/intermediate chemical exchange, our simulation showed a linear regression relationship between the optimal mean RF power for pulsed-CEST MRI and that of CW-CEST MRI. Worth noting, the optimized pulsed-CEST contrast was approximately equal to that of CW-CEST MRI for exchange rates below 50 s−1, as confirmed experimentally using a multi-compartment pH phantom. Moreover, acute stroke animals were imaged with both pulsed- and CW- amide protons CEST MRI, which showed similar contrast. In summary, our study elucidated the RF irradiation dependence of pulsed-CEST MRI contrast, providing useful insights to guide its experimental optimization and quantification.

Keywords: amide proton transfer (APT), chemical exchange saturation transfer (CEST), pH, radio frequency (RF), stroke


Chemical exchange saturation transfer (CEST) imaging is an emerging MRI technique capable of detecting dilute labile proton groups and microenvironment properties such as pH and temperature (13). For exchangeable proton groups undergoing slow and intermediate chemical exchange with bulk water protons, CEST contrast is nearly proportional to the labile proton concentration (f) and exchange rate (ksw), and therefore provides sizeable sensitivity enhancement for detecting dilute CEST agents (4,5). For instance, it has been shown that certain polypeptides and metabolites may be detected via CEST MRI (69). In addition, CEST MRI has been increasingly explored for in vivo applications (1012). Importantly, amide proton transfer (APT) imaging, a specific form of CEST MRI that probes the amide proton chemical exchange, is sensitive to local pH (13,14). As such, endogenous APT MRI has been increasingly used for imaging ischemic acidosis, particularly useful for characterizing impaired tissue metabolism during acute stroke, complementary to the commonly used perfusion and diffusion MRI (1518).

However, CEST MRI contrast is complex. In addition to labile proton concentration and exchange rate, the measureable CEST contrast also depends on experimental parameters such as magnetic field strength, RF irradiation power/amplitude, duration and scheme (19). More specifically, RF irradiation may not only saturate labile protons but also directly attenuate the bulk water signal (spillover), reducing the sensitivity and specificity of CEST MRI. Conventionally, CEST imaging has been implemented using fast image readout with a long continuous-wave (CW) RF irradiation after sufficient relaxation recovery, such that the steady state CEST contrast is obtained. In fact, mathematical solutions previously developed were mainly for CW-CEST MRI (1923). On the other hand, CW RF irradiation may not be applicable on systems where the maximal RF pulse duration and duty cycle are limited, particularly on clinical scanners, for which, the pulsed RF irradiation scheme must be applied instead (2427). Consequently, it is crucial to quantify and optimize pulsed-CEST MRI for in vivo applications.

For CW-CEST MRI, it has been shown that experimentally obtainable CEST contrast peaks at an optimal RF amplitude, under which the saturation efficiency and RF spillover factor are balanced (19). The goal of our study is to establish a quantitative algorithm to describe pulsed-CEST MRI. It is important to point out that whereas there is a single parameter to optimize for CW-CEST MRI (i.e., RF amplitude), multiple variables must be evaluated for pulsed-CEST MRI, including RF flip angle and pulse duration (Fig. 1). Because the transient solution is tedious and complex, our study modified the numerical solution, previously developed for CW-CEST MRI, to describe the pulsed-CEST MRI contrast (28,29). Specifically, we simulated pulsed-CEST MRI as a function of RF pulse train parameters (RF flip angle and duration) and CEST agent properties (exchange rate and chemical shift), assuming a commonly used Gaussian-shaped RF pulse. In addition, we compared the optimal equivalent RF power for pulsed-CEST MRI with that of CW-CEST MRI, and the simulation was validated with a multi-compartment pH phantom. Furthermore, pH-weighted endogenous APT MRI was obtained with both pulsed- and CW-CEST MRI sequence, which showed similar contrast, consistent with simulation and phantom results.

Fig. 1
Illustration of CW- and pulsed-CEST MRI sequences, which includes a long RF irradiation module followed by fast image readout. In comparison with the CW-CEST MRI, the pulsed-CEST MRI sequence has multiple variable parameters to optimize including the ...


Numerical Simulation

We simulated pulsed- and CW- CEST MRI in Matlab (MathWorks, Natick, MA) with modified Bloch-McConnell solution (20,21). Specifically, pulsed-CEST MRI was simulated piecewise for each repetition of RF pulse and interpulse delay. Each RF pulse was divided into 64 steps and the spin evolution within each step (Δτ=τp/64) was modeled assuming a constant B1 amplitude, B1(i)i=1:N=ϕγ·e(i·τp/Nτp/2)2/c2i=1:Ne(i·τp/Nτp/2)2/c2·τpN, where γ is the gyromagnetic ratio, ϕ is RF flip angle, and c=τp2ln(tl) with tl being the truncation level, chosen to be 1%. In addition, transverse magnetization was set to zero immediately before inter-pulse delay to represent the dephasing caused by crusher gradients, while its longitudinal magnetization relaxes towards its equilibrium state. The magnetization at the end of inter-pulse delay was then used as the initial state for iterated simulation till the end of RF pulse train. Moreover, we used representative relaxation parameters, with T1 and T2 of 3 s and 100 ms, respectively, for bulk water, and 1 s and 15 ms for the labile proton group. The labile proton concentration with respect to bulk water protons was 1:500. Finally, we assumed a total RF saturation time (TS) of 3 s and duty cycle of 50%, for a magnetic field strength of 4.7T.

RF pulse flip angle (Fig. 2)

Fig. 2
Evaluation of pulsed-CEST MRI contrast as a function of RF flip angle. a) Four simulated Z-spectra with RF flip angles of being 90°, 180°, 360° and 540°, respectively. We assumed Ksw=50 s−1, Δωs ...

Z-spectra were simulation for representative flip angles of 90°, 180°, 360° and 450°, assuming a commonly used Gaussian shaped pulse and a typical pulse duration of 15 ms. The chemical shift was varied from −5 to 5 ppm with intervals of 0.125 ppm. In addition, CEST contrast was solved as a function of RF flip angle, for three representative exchange rates: 25, 50 and 100 s−1. Moreover, the RF flip angle (ϕ) was varied from 150° to 250° with intervals of 10°, assuming a representative pulse duration of 15 ms and the optimal flip angle was solved as a function of exchange rate, for three typical chemical shifts. Finally, the optimal flip angle was evaluated as a function of chemical shift, from 1 to 5 ppm per 0.25 ppm, for three exchange rates, 25, 50 and 100 s−1.

RF pulse duration (Fig. 3)

Fig. 3
Evaluation of pulsed-CEST MRI contrast as a function of RF irradiation pulse duration. a) Three simulated Z-spectra with RF duration of 5, 15 and 30 ms. We assumed Ksw=50 s−1, Δωs=1.9 ppm and Φ=180°. b) Simulated ...

Z-spectra were simulation for pulse duration of 5, 15 and 30 ms, assuming a commonly used Gaussian shaped pulse with a flip angle of 180°, and the chemical shift was varied from −5 to 5 ppm with intervals of 0.125 ppm. In addition, CEST contrast was solved as a function of RF pulse duration, for three representative exchange rates of 25, 50 and 100 s−1. Moreover, the RF pulse duration was increased stepwise from 5 to 30 ms, with intervals of 1 ms and the optimal pulse duration was solved as a function of exchange rate for three typical chemical shift, 1.9, 3.5 and 5 ppm. Finally, the optimal pulse duration was evaluated as a function of chemical shift, from 1 to 5 ppm per 0.25 ppm, for three representative exchange rates, 25, 50 and 100 s−1.

Optimal RF flip angle and pulse duration (Fig. 4)

Fig. 4
Comparison of pulsed- and CW-CEST MRI contrast and their optimal experimental conditions. a) Although the optimized pulsed-CEST MRI contrast is approximately equal to that of CW-CEST MRI for slow exchange rates, it decreases when the exchange rate becomes ...

Both CW- and pulsed-CEST MRI were simulated with the exchange rate varied from 10 to 100 s−1, at intervals of 5 s−1, for three representative chemical shifts of 1.9, 3.5 and 5 ppm. For each exchange rate and chemical shift, we used 2-dimensional variants of the RF flip angle and pulse duration: 150° to 250° with an interval of 10°, and 5 to 25 ms with an interval of 1 ms, respectively.


We prepared a pH phantom using creatine and low gelling point (LGP) agarose, as reported previously (14,27). In short, we added 1% agarose into phosphate-buffered saline (PBS) solution, doped with a trace amount of CuSO4 · 5H2O (Sigma Aldrich, St Louis, MO). The mixture was microwave heated, and then immersed into a water bath at 50°C (Cole-Parmer, Vernon Hills, IL). When the temperature of the gel solution stabilized, we added creatine and mixed well until a final concentration of 50 mM. We serially titrated the pH of the creatine-gel solution to 5.5, 6, 6.25, 6.5, 6.75 and 7 at 50°C (EuTech Instrument, Singapore), and separately transferred the solution into multiple centrifuge tubes. The tubes were immediately sealed and inserted into a tight-fitting phantom holder, filled with 1% agarose PBS solution to minimize susceptibility mismatch. It is important to note that whereas the phantom pH at room temperature may be slightly different from values titrated at 50°C, the pH contrast remains among pH compartments and is suitable for evaluating pH-weighted CEST contrast (14).


Fig 1 illustrates a typical CEST MRI sequence, which includes a long RF saturation period followed by image acquisition. Under a long RF irradiation, the CW-CEST MRI contrast is at its steady state and the only variable is the RF amplitude/power (ω1). In comparison, pulsed RF irradiation scheme is complex, with multiple parameters to optimize, including pulse shape, flip angle (ϕ), duration (τp) and interpulse delay (τd). In addition, crusher gradients are often applied between pulses to suppress residual transverse magnetization. As Gaussian pulse is widely used, and has been previously explored for CEST MRI (13,18), our current study focused on Gaussian pulsed-CEST MRI. In addition, we studied a representative RF duty cycle of 50%, namely τp= τd (26).

All images were acquired at 4.7 T (Bruker Biospec, Billerica, MA), which is capable of performing both CW and pulsed RF irradiation. We used single-shot single-slice echo planner image (EPI) as image readout (slice thickness = 10 mm, field of view (FOV) = 76 × 76 mm, bandwidth=200 kHz). The image matrix was 64 × 64, and the repetition time (TR) and echo time (TE) were 10 s and 28 ms, respectively, and the RF saturation time (TS) of 5 s (NA=2). We probed creatine amine labile protons at 6.6 ppm (1.9 ppm from bulk water resonance) for CEST MRI. For CW-CEST MRI, the RF amplitude was varied from 0.75, 1, 1.25, 1.5 and 2 μT. For pulsed-CEST MRI, the RF pulse duration was varied from 10, 15, 20, 25 to 30 ms, and for each pulse duration, the flip angle was varied from 90° to 540°, with an interval of 45°. The CEST contrast was calculated as CEST asymmetry CESTRasym=(Iref−Ilabel)/I0, where Iref and Ilabel are reference and saturation scans, and I0 is the control scan without CEST RF irradiation. In addition, B1 field was calibrated with a 1 ms RF pre-pulse, whose flip angle (θ) was varied from 10° to 180° with an interval of 10°, and the image intensity was solved by fitting I(θ) = I0 · |cosγ(ηB1 + ΔB1)τ|, where η and ΔB1 are the scaling factor and shift of B1 field, respectively. We found η=0.95 and ΔB1≈0.

All animal experiments were carried out in accordance with guidelines approved by the Subcommittee on Research Animal Care of the Massachusetts General Hospital (SRAC, MGH). Standard filament middle cerebral artery occlusion (MCAO) was induced in anesthetized adult male Wistar rats (n=3, Charles River, Wilmington, MA), with both heart rate and blood SpO2 monitored online (Nonin Pulse Oximeter 8600, Plymouth, MN). In addition, body temperature was maintained within the normal physiological range with a circulating warm water jacket positioned above the torso. Multislice MRI (5 slices, slice thickness/gap = 1.8/0.2 mm, field of view = 25×25 mm2, acquisition matrix = 64×64) was obtained with single-shot echo-planar imaging (EPI) (receiver bandwidth = 200 kHz). Endogenous APT MRI was obtained using unevenly segmented RF irradiation scheme (Tr/TE/TS1/TS2=5000/15/4500/500ms), with both CW- and pulsed- CEST MRI (30). For APT MRI with CW-CEST acquisition, we had B1=0.75 μT (17), while for the pulsed-CEST MRI, we used τ=15 ms and ϕ=180°. In addition, RF irradiation was applied at ±3.5 ppm (700 Hz at 4.7T) and CEST asymmetry was calculated. Moreover, isotropic diffusion weighted MRI was obtained with two b values (250 and 1,000 mm2/s, TR/TE=3250/54ms, NA=16), and apparent diffusion coefficient (ADC) was solved by least square fitting.


We evaluated pulsed-CEST MRI contrast as a function of RF pulse flip angle. Fig. 2a shows four simulated Z-spectra, with the RF flip angle being 90°, 180°, 360° and 540°, assuming an RF pulse duration of 15 ms. These results show that whereas prominent CEST contrast can be observed with moderate flip angles (90° and 180°), the contrast decreases significantly at larger flip angles such as 360° and 540°, likely due to concomitant RF spillover effect. It is important to note that the water signal on resonance with RF saturation deviated from zero at the 360°flip angle because the RF pulse restores the bulk water magnetization along the Z-axis before the crusher gradient. Fig. 2b shows that for three representative exchange rates of 25, 50 and 100 s−1. The pulsed-CEST MRI contrast peaks when the RF flip angle is about 180–220°, but it bottoms out when the flip angle is beyond 360°. We solved the optimal RF flip angle for a typical range of exchange rate from 10 to 100 s−1. Fig. 2c shows that the optimal flip angle increases slightly with the exchange rate, and similarly, with chemical shift (Fig. 2d).

We further evaluated how the pulsed-CEST MRI contrast varies with the RF pulse duration, for a representative RF flip angle of 180°. Fig. 3a shows three Z-spectra for RF pulse durations of 5, 15 and 30 ms. It is worth pointing out that the RF saturation bandwidth is significantly broadened at short pulse duration, which results in greatly attenuated signal intensity around the bulk water resonance. Fig. 3b shows pulsed-CEST contrast as a function of RF pulse duration for three representative exchange rates of 25, 50 and 100 s−1. The optimal RF duration is prolonged at slower exchange rate, suggesting lower RF amplitude, similar to findings from CW-CEST MRI studies (19). We also evaluated how the optimal RF duration varies with exchange rate and chemical shift. Our results showed that the optimal RF pulse duration shortens with both exchange rate (Fig. 3c) and chemical shift (Fig. 3d). This is expected given that for the same RF flip angle (i.e., 180°), the equivalent RF power/amplitude is higher at shorter RF duration, and hence, is suitable for imaging CEST agents of higher exchange rate and larger chemical shift. Given that RF bandwidth is associated with the pulse duration, our results suggest that the optimal RF bandwidth also varies with exchange rate and chemical shift. Specifically, the bandwidth is 1610 Hz for a 1 ms Gaussian inversion pulse (RF Shape tool, Bruker Biospin, Billerica MA), which allows conversion of the optimal RF duration into RF bandwidth.

We also compared the pulsed-CEST contrast with the more commonly used CW-CEST MRI. Fig. 4a shows that for a representative chemical exchange rate range from 10 to 100 s−1, the CW-CEST contrast is approximately equal to or higher than the pulsed-CEST contrast. We plotted the difference between pulsed-CEST and CW-CEST contrast as a function of exchange rate, which showed that the pulsed-CEST MRI became less effective than CW-CEST MRI, particularly at higher exchange rate. In addition, the loss of contrast is greater at higher chemical shift (Fig. 4b). Fig. 4c shows the optimal RF amplitude for the CW-CEST MRI increases with exchange rate for three representative chemical shifts, consistent with the notion that larger RF power is needed to efficiently saturate faster exchange. In addition, the equivalent RF amplitude of pulsed-CEST MRI is calculated as B1 = [var phi]/(γτp). Fig. 4d compares the equivalent RF amplitude of pulsed-CEST MRI with that of CW-CEST MRI. Importantly, the RF amplitude for pulsed-CEST MRI displayed an excellent polynomial relationship with that of CW-CEST MRI, found to be B1_pulse_opt0.28B1_cw_opt2+1.04B1_cw_opt. This linear regression relationship allows estimation of the optimal pulsed-CEST MRI parameter based on the better-understood CW-CEST MRI.

We validated the simulation using a creatine-gel phantom with serially titrated pH, whose chemical exchange rate covers a representative range of diamagnetic CEST agents (Fig. 5a). CEST MRI indeed detected the difference in pH among the compartments, as per the CW-CEST and pulsed-CEST MRI shown in Figs. 5b and 5 c, respectively. For CW-CEST MRI, the RF amplitude is 1.25 μT, while for the pulsed-CEST MRI, we used an RF flip angle of 257° (after RF field calibration) and pulse duration of 20 ms. It is important to note that pulsed-CEST MRI showed noticeably lower contrast than the CW-CEST MRI, particularly for higher exchange rates/pH. Fig. 5d shows that pulsed-CEST contrast initially increased with RF flip angle, but decreased after peaking around an intermediate flip angle. Similarly, pulsed-CEST MRI contrast increased with pulse duration, but decreased when pulse duration became too long (Fig. 5e). Fig. 5f shows a 2-D histogram of optimal RF flip angle and pulse duration for all six pH compartments, and optimal flip angle and pulse duration were found to be approximately 257°(after RF field calibration) and 20 ms, respectively. The optimal flip angle was slightly larger than that estimated from simulation(Fig 2c), likely caused by the relatively coarse intervals of pulse flip angle and duration. In addition, the chemical exchange rate can be estimated using a previously calibrated formula for base-catalyzed creatine amine proton exchange (Fig. 5g) (14). Fig. 5h showed that CW-CEST contrast increases approximately linearly with exchange rate, suggesting excellent pH response. In comparison, whereas the pulsed-CEST MRI contrast was approximately equal to that obtained with CW-CEST MRI for exchange rates below 50 s−1, it became significantly less than the CW-CEST contrast at higher exchange rates/pH (Fig. 5i). The difference in CEST contrast between pulsed-CEST and CW-CEST MRI was evaluated as ΔCESTR= CESTRpulse − CESTRcw, which showed sizeable loss of contrast for exchange rates above 50 s−1, as predicted in Fig. 4.

Fig. 5
Experimental evaluation of pulsed-CEST MRI. a) Illustration of multi-compartment creatine-gel pH phantom. b) pH-weighted CW-CEST MRI of the phantom (B1=1.25 μT). c) pH-weighted pulsed-CEST MRI of the phantom (ϕ=257°(after RF field ...

Finally, we evaluated pH-weighted endogenous APT MRI with both pulsed- and CW-CEST MRI sequences, using acute ischemic stroke animals. Ischemic lesion can be appreciated as hypointensity in ADC (Fig. 6a), CW- (Fig. 6b, B1=0.75 μT) and pulsed-APT (Fig. 6c, τ=15 ms and ϕ=180°) images, from a representative coronal slice. A region of interest (ROI) was manually selected in the striatum region based on ADC decrease, being 0.47 ± 0.08 μm2/ms. Its APT contrast was found to be −6.4 ± 1.4% (mean ± S.D.) and −5.9 ± 1.2%, for pulsed- and CW-CEST acquisitions, respectively. Whereas the pulsed-APT is slightly lower than that of CW-APT MRI, it provides reasonably clear depiction of ischemic lesion, and therefore serves as a viable alternative for MRI systems unable to operate CW RF irradiation.

Fig. 6
in vivo evaluation of pH-weighted APT MRI. MCAO induced ischemic tissue damage can be observed in ADC map (a), which also displays hypointensity in pH-weighted APT maps acquired with CW-CEST (b, B1=0.75 μT) and pulsed-CEST (c, τ=15 ms, ...


Our study developed numerical simulation to elucidate pulsed-CEST MRI. We showed that whereas the optimal experimental conditions for RF flip angle and duration vary with chemical shift and exchange rate, they can be reasonably estimated based on the better-understood CW-CEST MRI. Specifically, for typical chemical shift (< 5ppm) and slow exchange rate, the optimal RF flip angle is approximately 180°, and the RF pulse duration can be derived from the optimal mean RF amplitude for pulsed-CEST MRI based on a linear regression relationship with that of CW-CEST MRI. Moreover, our study simulated CEST MRI for a range of chemical shifts (i.e., 1–5 ppm) assuming a typical field strength of 4.7T; the results we obtained can be easily extended to other field strengths such as 3T. For instance, the creatine amine proton chemical shift (1.9 ppm) was 380 Hz at 4.7T, similar to the composite amide proton chemical shift (3.5 ppm) of 450 Hz at 3T. It is necessary to point out that optimization of pulsed-CEST MRI with numerical simulation is very time-consuming given that the RF saturation scheme must be simulated piecewise; the results shown in Fig. 4 took over 10 hours on a Dell Precision T7400 Desktop due to multi-dimensional nature of the simulation (involving variable exchange rate, chemical shift, RF flip angle and duration). Therefore, by elucidating the relationship between pulsed- and CW-CEST MRI, our study has enabled reasonable estimation of the optimal pulsed RF saturation scheme, which may effectively guide experimental optimization.

It is very interesting to observe the linear regression relationship between the mean RF amplitude of pulsed-CEST MRI and that of CW-CEST MRI. For a representative CW RF amplitude of 1 μT, the equivalent B1 amplitude for optimal pulsed-CEST MRI will be 1.34 μT. Given a typical RF duty cycle of 50%, the ratio between the flip angle of pulsed-CEST MRI and CW-CEST MRI within a repetition of RF pulse irradiation and interpulse delay was 67%. However, the ratio between pulsed- and CW-CEST MRI RF irradiation power can be estimated using p=t=0τp+τdB12(t)dt, which by our calculation yielded 211%. As such, the optimal pulsed-CEST experimental condition seems to compromise between RF amplitude and power, which indicates that the contrast mechanism of pulsed-CEST MRI is complex and before an exact analytical solution can be obtained, it is necessary to optimize it numerically. Whereas CEST contrast varies strongly with bulk water T1, the optimal RF power showed very little T1 dependence (data not shown). This is not surprising since the optimal RF power depends on the competition of labeling coefficient and spillover factor. In addition, although bulk water T2 may influence the RF spillover effect, our simulation found that it has relatively small effect upon the optimal RF power, for moderate chemical shifts. Nevertheless, our results showed that the optimal pulsed-CEST MRI condition can be reasonably estimated from that of better understood CW-CEST MRI, which may help guide experimental optimization.

Our study showed that although the pulsed-CEST MRI contrast is persistently lower than that of optimized CW-CEST MRI due to the less-efficient saturation efficiency and loss of CEST contrast during the interpulse delay, pulsed-CEST MRI provided comparable contrast for an exchange rate less than 50 s−1. This observation is encouraging for clinical translation of pH-weighted APT MRI. Specifically, the endogenous amide proton exchange rate has been estimated to be about 30 s−1, and therefore, is suitable for pulsed-APT MRI (13). While on the other hand, it is important to point out that pulsed-CEST MRI may not be optimal for imaging CEST agents of higher exchange rate such as certain paramagnetic CEST (PARACEST) agents (31,32). To address this limitation, we need to either implement CW-RF irradiation or design novel RF irradiation scheme. Moreover, although our study evaluated Gaussian shaped RF pulse, which has been commonly used due to its simplicity and smooth frequency profile, the algorithm we developed can be readily extended to investigate RF pulse of a given shape for pulsed-CEST MRI (13,18). Indeed, we also evaluated block and sinc RF pulses, both of which are less effective than Gaussian shaped RF pulses. Specifically, block RF pulse induces a periodic oscillation pattern in Z-spectrum due to its sidelobes in the frequency domain. For the sinc RF pulse we simulated, its peak power is significantly higher than that of Gaussian or block pulses, thereby, more susceptible to RF spillover effects. Therefore, our results suggest that the frequency spectrum of the irradiation pulse should be reasonably smooth yet narrow for pulsed-CEST MRI. Whereas we can not be exhaustive in evaluating all commonly used RF pulse shapes, our study showed that for imaging of diamagnetic CEST agents undergoing slow/intermediate chemical exchange, pulsed-CEST MRI with optimized Gaussian RF irradiation provides similar contrast as CW-CEST MRI.


Our study developed a numerical method and elucidated the pulsed-CEST MRI contrast. We showed that the optimal RF irradiation scheme for pulsed-CEST MRI can be reasonably estimated from the better-understood CW-CEST MRI, which greatly simplifies numerical optimization and experimental validation of pulsed-CEST MRI. In addition, we showed that for diamagnetic CEST agents undergoing slow chemical exchange, the optimized pulsed-CEST MRI contrast is comparable to that of CW-CEST MRI. The simulation has been further validated with both phantom and in vivo applications. In summary, our results showed that pulsed-CEST MRI is feasible with optimized RF irradiation, and remains promising for clinical translation.


This study was supported in part by grants from AHA/SDG 0835384N, NIH/NIBIB 1K01EB009771, NIH/NINDS R21NS061119 and NIH/NCRR P41RR14075. This work has been presented at the Gordon Research Conference of in vivo MR, July 2010.


1. Ward KM, Aletras AH, Balaban RS. A new class of contrast agents for MRI based on proton chemical exchange dependent saturation transfer (CEST) J Magn Reson. 2000;143:79–87. [PubMed]
2. Ward KM, Balaban RS. Determination of pH using water protons and chemical exchange dependent saturation transfer (CEST) Magn Reson Med. 2000;44:799–802. [PubMed]
3. Zhang S, Malloy CR, Sherry AD. MRI Thermometry Based on PARACEST Agents. J Am Chem Soc. 2005;127(50):17572–17573. [PMC free article] [PubMed]
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. Magn Reson Med. 2003;49:998–1005. [PubMed]
5. Terreno E, Castelli D, Aime S. Encoding the frequency dependence in MRI contrast media: the emerging class of CEST agents. Contrast Media Mol Imaging. 2010;5(2):78–98. [PubMed]
6. Zhang S, Trokowski R, Sherry AD. A Paramagnetic CEST Agent for Imaging Glucose by MRI. J Am Chem Soc. 2003;125(50):15288–15289. [PubMed]
7. Aime S, Delli Castelli D, Fedeli F, Terreno E. A paramagnetic MRI-CEST agent responsive to lactate concentration. J Am Chem Soc. 2002;124(32):9364–9365. [PubMed]
8. Gilad AA, McMahon MT, Walczak P, Winnard PT, Jr, Raman V, van Laarhoven HW, Skoglund CM, Bulte JW, van Zijl PC. Artificial reporter gene providing MRI contrast based on proton exchange. Nat Biotechnol. 2007;25(2):217–219. [PubMed]
9. McMahon MT, Gilad AA, DeLiso MA, Berman S, Bulte JWM, van Zijl PCM. New “multicolor” polypeptide diamagnetic chemical exchange saturation transfer (DIACEST) contrast agents for MRI. Magn Reson Med. 2008;60(4):803–812. [PMC free article] [PubMed]
10. Zhou J, Lal B, Wilson DA, Laterra J, van Zijl PCM. Amide proton transfer (APT) contrast for imaging of brain tumors. Magn Reson Med. 2003;50:1120–1126. [PubMed]
11. Vinogradov E, He H, Lubag A, Balschi JA, Sherry AD, Lenkinski RE. MRI detection of paramagnetic chemical exchange effects in mice kidneys in vivo. Magn Reson Med. 2007;58(4):650–655. [PubMed]
12. van Zijl PCM, Jones CK, Ren J, Malloy CR, Sherry AD. MRI detection of glycogen in vivo by using chemical exchange saturation transfer imaging (glycoCEST) Proc Natl Acad Sci. 2007;104(11):4359–4364. [PubMed]
13. Zhou J, Payen JF, Wilson DA, Traystman RJ, van Zijl PC. Using the amide proton signals of intracellular proteins and peptides to detect pH effects in MRI. Nat Med. 2003;9(8):1085–1090. [PubMed]
14. Sun PZ, Sorensen AG. Imaging pH using the Chemical Exchange Saturation Transfer (CEST) MRI: Correction of Concomitant RF Irradiation Effects to Quantify CEST MRI for Chemical Exchange Rate and pH. Magn Reson Med. 2008;60(2):390–397. [PubMed]
15. Hossmann KA. Viability thresholds and the penumbra of focal ischemia. Ann Neurol. 1994;36(4):557–565. [PubMed]
16. Sun PZ, Zhou J, Sun W, Huang J, van Zijl PC. Detection of the ischemic penumbra using pH-weighted MRI. J Cereb Blood Flow Metab. 2007;27(6):1129–1136. [PubMed]
17. Sun PZ, Zhou J, Huang J, van Zijl P. Simplified Quantitative Description of Amide Proton Transfer (APT) Imaging During Acute Ischemia. Magn Reson Med. 2007;57(2):405–410. [PubMed]
18. Jokivarsi KT, Gröhn HI, Gröhn OH, Kauppinen RA. Proton transfer ratio, lactate, and intracellular pH in acute cerebral ischemia. Magn Reson Med. 2007;57(4):647–653. [PubMed]
19. Sun PZ, van Zijl PCM, Zhou J. Optimization of the irradiation power in chemical exchange dependent saturation transfer experiments. J Magn Reson. 2005;175(2):193–200. [PubMed]
20. 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. Magn Reson Med. 2005;53(4):790–799. [PubMed]
21. Sun PZ. Simplified and scalable numerical solution for describing multi-pool chemical exchange saturation transfer (CEST) MRI contrast. J Magn Reson. 2010;205(2):235–241. [PMC free article] [PubMed]
22. Li AX, Hudson RHE, Barrett JW, Johns CK, Pasternak SH, Bartha R. Four-pool modeling of proton exchange processes in biological systems in the presence of MRI-paramagnetic chemical exchange saturation transfer (PARACEST) agents. Magn Reson Med. 2008;60(5):1197–1206. [PubMed]
23. McMahon M, Gilad A, Zhou J, Sun PZ, Bulte J, van Zijl PC. 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. [PMC free article] [PubMed]
24. Gochberg DF, Kennan RP, Roberson PL, Gore JC. Quantitative imaging magnetization transfer multiple selective pulses. Magn Reson Med. 1999;41:1065–1072. [PubMed]
25. Smith SA, Farrell JAD, Jones CK, Reich DS, Calabresi PA, van Zijl PCM. Pulsed magnetization transfer imaging with body coil transmission at 3 Tesla: Feasibility and application. Magn Reson Med. 2006;56(4):866–875. [PubMed]
26. Sun PZ, Benner T, Kumar A, Sorensen AG. Investigation of optimizing and translating pH-sensitive pulsed-chemical exchange saturation transfer (CEST) imaging to a 3T clinical scanner. Magn Reson Med. 2008;60(4):834–841. [PMC free article] [PubMed]
27. Mougin OE, Coxon RC, Pitiot A, Gowland PA. Magnetization transfer phenomenon in the human brain at 7 T. NeuroImage. 2010;49(1):272–281. [PubMed]
28. Helms G, Dathe H, Hagberg GE. Pulsed saturation of the standard two-pool model for magnetization transfer. part II: the transition to steady state. Concepts Magn Reson A. 2004;21:50–62.
29. Murase K, Tanki N. Numerical solutions to the time-dependent Bloch equations revisited. Magn Reson Imaging. 2011;29(1):126–131. [PubMed]
30. Sun PZ, Cheung JS, Wang E, Benner T, Sorensen AG. Fast multi-slice pH-weighted chemical exchange saturation transfer (CEST) MRI with unevenly segmented RF irradiation. Magn Reson Med. 2010 in press. [PMC free article] [PubMed]
31. Aime S, Calabi L, Biondi L, Miranda MD, Ghelli S, Paleari L, Rebaudengo C, Terreno E. Iopamidol: Exploring the potential use of a well-established x-ray contrast agent for MRI. Magn Reson Med. 2005;53(4):830–834. [PubMed]
32. Dixon TW, Ren J, Lubag A, JM, Ratnakar J, Vinogradov E, Hancu I, Lenkinski RE, Sherry AD. A concentration-independent method to measure exchange rates in PARACEST agents. Magn Reson Med. 2010;63(3):625–632. [PMC free article] [PubMed]