Synopsis
Accurate arterial input function (AIF) measurement in Dynamic Contrast
Enhanced MRI (DCE-MRI) remains challenging. This hinders DCE-MRI’s wider
adoption. Since the contrast reagent
(CR) is detected indirectly through water proton R1 relaxation rate
constant change, DCE-MRI intrinsically works as a dual-probe (CR and water)
method. In this study, we demonstrate
that while the common pharmacokinetic parameters associated with CR
extravasation are highly sensitive to AIF accuracy, the transcytolemmal water
exchange parameter is not. With the
recent correlation of water exchange kinetics and cellular metabolic activity,
this current work demonstrates the practicability of high-resolution metabolic imaging of the
prostate. Purpose
Accurate arterial input function
(AIF) measurement in prostate Dynamic Contrast Enhanced Magnetic Resonance
Imaging (DCE-MRI) remains challenging. This hinders DCE-MRI’s wider clinical adoption.
Since Gadolinium based contrast reagent (CR) is detected indirectly through
water proton R
1 relaxation rate constant change, DCE-MRI is
intrinsically a dual-molecular probe (CR and water) technique. In this study, we demonstrate that while the
most common pharmacokinetic model parameters like K
trans (CR
transfer constant) and v
e (extravascular, extracellular, volume
fraction) are highly AIF-dependent for prostate DCE‑MRI, the transcytolemmal water
exchange parameter is not. With the recent correlation of water exchange kinetics and cellular
metabolic activity, this work demonstrates the practicability of a
high-resolution metabolic imaging approach for the prostate.
Methods
Prostate DCE-MRI data were acquired on
13 subjects with a Siemens TIM Trio (3T) system under an IRB-approved
protocol. RF receiving used a combination
of Spine Matrix and flexible Body Matrix coil arrays. The DCE-MRI acquisition employed a 3D
TurboFLASH pulse sequence with a 256*144*16 matrix size and a 360*203 mm
2
FOV; slice thickness 3 or 3.2 mm. Imaging
parameters are: TR/TE/FA, 5.0 ms/1.57 ms/15º; and image frame sampling interval
is 6.3 s. Other details are reported
previously.
1 All subjects subsequently underwent standard
ultrasound-guided systematic ten-core
prostate biopsies. Malignancy was found
in 5 subjects. One region of
interest (ROI) was selected for each subject, resulting in 5 malignant and 8 benign
ROI time-courses.
1
The water exchange pharmacokinetic model
assuming a single compartment of water signal
2 extracts K
trans,
v
e, and τ
i (mean intracellular water lifetime) parameters. Monte Carlo simulations were used to
investigate parameter changes associated with AIF magnitude variation. The “true” AIF (AIF
0) was that determined
directly from the data
1. A
noiseless tissue time-course was first generated based on data acquisition
detail and pharmacokinetic equations. Simulations then fitted pharmacokinetic
parameters to the time-course with various AIF amplitudes scaled from 0.8 to
2.0 of AIF
0 (step-size 0.1). For
each amplitude-adjusted AIF (e.g., AIF
1
= 130%∙AIF
0), four hundred fittings were
performed with random Gaussian noise
added to the simulated time-course and a different initial guess parameter
value set each time.
Results
The
Figure 1 bar graph shows the K
trans,
v
e, and τ
i accuracy (mean)
and precision (SD) changes responding to a 30% AIF amplitude increase (AIF
1). K
trans
= 0.4 min
-1, ve = 0.25, and τ
i = 0.45 s base values were used to generate the simulated tissue DCE
time-course. K
trans and v
e
numerical values decreased nearly 30% with a 30% AIF magnitude increase,
while the change for τ
i remains negligibly small (< 3%). The τ
i reciprocal (τ
i-1)
is the unidirectional rate constant for equilibrium
cellular water efflux (k
io).
4
Figure 2a shows a portion of an axial T
2-weighted
pelvic image (anterior up) of a 66 year old subject with left prostate
malignancy.
Figures 2b – 2d
show zoomed parametric maps of a tumor-encompassing ROI in the left prostate
indicated by the
2a yellow rectangle. In
Figure 3 lesion ROI-averaged
in
vivo k
io values are plotted against subsequently obtained
ex vivo Gleason scores, for the five
malignant cases. A clear trend of
increasing k
io with Gleason score (Spearman’s correlation
coefficient, r = 0.908; p = 0.0333) is observed (the blue dashed curve is
intended only to guide the eye). The
error bars show k
io standard errors and estimated Gleason score
errors.
3 Discussion
Although the AIF is hard to measure accurately from the CR bolus
first-pass, an AIF error of less than 30% can be readily managed from the
washout phase, taking into account the total blood volume and CR elimination
rate. Furthermore, if k
io increases are dominated by increases in the
mean Na+,K+-ATPase [NKA] turnover for the cells in the
voxel,
4 this work suggests significant impact and unprecedented
insight into prostate tumor progression in
vivo. Active trans-membrane water
cycling accompanies the continuous, homeostatic active trans-membrane ion and
osmolyte cycling enzymatic activity, driven by the membrane P-type ATPase ion
pump, constitutively characteristic of a living cell.
4 Fig. 3 indicates a greater Gleason score
correlates with greater mean NKA activity.
Taken together, the k
io map shows
in vivo high-resolution metabolic imaging that is insensitive to
AIF uncertainty (thus easy to quantify in that sense). It is important to point out that precision
in the k
io map is greatest when CR extravasation is extensive (large
K
trans),
4 as always seen in prostate.
1 Acknowledgements
Grant
Support:
NIH:
RO1-NS040801, UO1-CA154602. R44 CA180425. References
1.
Li, et al. Magn Reson Med., 69:171-178,
2013. 2. Yankeelov, et
al. Magn Reson Med ., 50:1151-1169,
2003. 3. Goodman et al.
The Prostate 72:1389-1398, 2012. 4. Springer, et al.
NMR Biomed. 27:760-73, 2014.