|
|
||||||||
Basic Science Investigations |
Department of Radiology, University of Washington, Seattle, Washington
| ABSTRACT |
|---|
|
|
|---|
Key Words: 3'-deoxy-3'-fluorothymidine thymidine kinase 1 compartmental modeling kinetic analysis
| INTRODUCTION |
|---|
|
|
|---|
FLT is transported between blood and cells by the same nucleoside transporters as thymidine (11), along with a small amount of facilitated diffusion. Once intracellular, FLT tracks the salvage pathway of DNA synthesis and undergoes the same early fate as thymidine, which is phosphorylation by TK1. FLT directly competes with thymidine for TK1 phosphorylation. TK1 activity is directly linked to cellular biosynthesis, where it has been shown to increase 10-fold during the DNA synthetic (S phase) part of the cell replication cycle (1214). The phosphorylated product of TK1 metabolism, FLT-monophosphate (FLTMP), is retained in tissues (15) and may continue to be phosphorylated to FLT-diphosphate (FLTDP) and FLT-triphosphate (FLTTP) (Fig. 1). At this point in the pathway to DNA, FLT differs from thymidine. FLTTP is not incorporated into DNA (11) because of the 3' substituent and its metabolism goes no further in the exogenous pathway. Because FLTTP is not incorporated into DNA, the rate-limiting step for FLT retention is the initial phosphorylation by TK1 (15). Phosphorylated FLT can be dephosphorylated back through the pathway to FLT by 5'-deoxynucleotidase enzymes. However, this occurs at a slow rate relative to TK1 activity (11,16,17). Thus, the retention of FLT nucleotides is a measure of the TK1 activity.
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
FLT is not a substrate for thymidine phosphorylase, the enzyme that breaks the bond between the pyrimidine and deoxyribose. However, FLT is glucoronidated in the liver and then delivered to the blood. FLT-glucuronide is restricted to the vascular space and is cleared by the kidneys. A sizeable fraction of the total plasma activity appears in the form of FLT-glucoronide late in the imaging sequence, so that the total plasma activity during the imaging study becomes a combination of FLT and FLT-glucoronide. No other labeled metabolites of FLT have been observed in human plasma (8,11,20). The proportion of FLT present as the glucuronide in plasma can be assessed by chromatographic analysis. The information is used to derive the metabolite-corrected time course of FLT activity in plasma. The whole blood and plasma activity concentrations have been reported to be statistically identical, showing no partitioning (8,10). The FLT plasma clearance curve provides the input to the compartmental model to characterize FLT tissue uptake. The total blood activity curve is also used in the modeling process to account for the contamination of vascular activity within each PET image region.
The essential features of the model are founded on the biochemical behavior of FLT described here. The application of our model to the estimation of regional cellular proliferation requires some additional assumptions:
An essential premise of the imaging experiment is that FLTMP and other FLT-phosphates are trapped in the tissue for the duration of the imaging study. Grierson et al. (17) observed in a lung carcinoma cell line (A549) that the phosphorylated products are essentially retained in the cell with negligible activity measured in the extracellular fluid. Basic cell studies suggest that phosphorylated intermediates, such as nucleotides, would not be observed in blood (17). Our group has also measured the slow rate of dephosphorylation of FLTMP by dNT (17). Others have characterized the relative rate of phosphorylation by TK1 for TdR and FLT, providing biochemical data for modeling FLT (14,21,22).
Anytime after injection of FLT, the total concentration of FLT in tissue, Ct, is comprised of the sum of the amounts of free FLT in the precursor pool, Qe, and in phosphorylated products of FLT, Qm:
![]() | (Eq. 1) |
is the tissue-specific gravity (g/mL).
The rate of change of the amount of free FLT in the tissue, dQe/dt, is equal to the difference between the rate of accumulation of FLT by transport into tissue from plasma or from dephosphorylation of FLTMP, the rate of loss of FLT by transfer of intracellular FLT back to the plasma, or phosphorylation by TK1 to provide FLTMP. The following equation applies:
![]() | (Eq. 2) |
The initial transfer constant in the model, K1, reflects the net transport of FLT from blood through various barriers into cells. There is negligible protein binding of FLT or TdR in the blood (23). CpFLT is the average capillary concentration of FLT and is assumed to be the FLT concentration in peripheral arterial plasma. A similar assumption that arterial plasma thymidine concentration reflects the capillary concentration was used in the determination of the synthetic rate of DNA from TdR uptake (18). K1 is comprised of blood flow, diffusion into the interstitial space, and nucleoside membrane transport (2426). In somatic tissues, transport of thymidine into the cell is rapid (27) and K1 reflects primarily blood flow (18,28). Transfer of FLT is similar to the transport of thymidine and involves facilitated diffusion through both equilibrative and concentrative transporters (26,29). Equilibrative transporters, predominantly found in human tumor cell lines, are thought to be the primary process mediating the facilitated influx and efflux of nucleosides (24,29).
FLT can also leave the tissue by way of a similar transport process in the reverse direction, k2 in the model. However, once FLT is in the tissue, it can also be phosphorylated to FLTMP by TK1. This process is represented in the model by transfer from Qe to the compartment Qm. The transfer constant for this phosphorylation reaction is given by k3 and characterizes the rate-limiting step in the uptake and retention of FLT (11,15,30). Qm includes the additional trapped intermediates FLTDP and FLTTP (Fig. 1), but these are not rate-limiting steps. These intermediates may also be dephosphorylated back through the pathway to FLTMP. The rate constant k4 describes the reverse reaction catalyzed by dNT to yield FLT, which then has the possibility of leaving the PET ROI mostly through the same equilibrative transporter or passive diffusion process. The mathematic description of the rate of change of FLT nucleotides as a function of time, dQm/dt, is given by the following formula:
![]() | (Eq. 3) |
Equations 13 provide the basis for the determination of the rate constants K1, k2, k3, and k4. They describe the time course of the total amount of radioactivity in the tissue as a function of time. The total concentration of FLT predicted by the model (Ct) is fitted to the measured concentration of FLT for each time point from the dynamic PET sequence, Ci. A nonlinear least-squares optimization process based on the Levenberg-Marquart approach computes the best-fitting values of the kinetic parameters. This method minimizes the weighted sum of the square error (WSSE) using a weighting function (
) derived from the total counts of a PET image ROI in a time frame and the duration of the time frame (31):
![]() | (Eq. 4) |
The FLT flux estimated by the model is given by the product of the transport rates of FLT from blood to tissue and phosphorylation of FLT by TK1 divided by the sum of the rates of efflux from the exchangeable tissue compartment. This FLT model composition is similar to the FDG model proposed by Sokoloff et al. (32) and the TdR model proposed by Mankoff et al. (18). By analogy, the flux constant for the proposed FLT model is calculated using Equation 5:
![]() | (Eq. 5) |
The model of FLT kinetics described here assumes no local metabolism of FLT other than phosphorylation. Thus, the model would be inappropriate for the liver, where metabolism produces FLT-glucoronide that is able to leave the tissue region and enter the blood. Under these conditions, the assumption that retention is correlated with TK1 activity would not hold, and an alternative model would have to be proposed and validated.
Parameter Ranges and Typical Values
In any compartmental model, it is necessary to establish the range of realistic parameter values expected in typical imaging conditions (Table 2). Of the parameters used to calculate KFLT in Equation 5, the expected values can be estimated from previously published rates. The transport of FLT, K1, can be estimated from the transfer constant of thymidine from blood through several barriers into the tissue. We previously estimated a typical K1 for thymidine in proliferating tissue of 0.20 mL/min/g, with a range of 0.050.30 mL/min/g (18). Transport of FLT has been estimated in cell experiments to be approximately 0.08 mL/min/g (15), and in vivo in dogs it ranged from 0.14 to 0.22 mL/min/g (9). Because the 3 tissues examined in this study potentially have great differences in FLT uptake, the parameter ranges are expected to be wider than previous model implementations using thymidine. Therefore, we can expect the range of K1 to be wider for FLT than for TdR, but the initial model parameter should be less than TdR. We therefore use a typical value for K1 of 0.1 mL/min/g as the initial starting parameter, with a range of 0.050.5 mL/min/g.
|
The transfer constant k3 represents the rate-limiting step of FLT retention via the salvage pathway: FLTMP production from FLT by TK1. This is different from the thymidine model, in which the incorporation of thymidine triphosphate (TTP) into DNA is the rate-limiting step along the pathway to DNA synthesis. Based on the in vitro relationship of the Michaelis-Menten constants (Km, Vmax) between TdR and FLT with respect to TK1 activity (21,22), the phosphorylation ratio (PRTK1) of TK1 can be measured and used as a starting point for estimating the relative rate of metabolism of FLT by TK1 (33,34). The PR for TK1 can be calculated from earlier reports as 0.08 (21). From this relationship, the starting value of k3 was set at 0.1 min1 (range, 0.0011.0 min1).
The activity in the metabolic compartment, Qm, includes all 3 FLT nucleotides. FLTDP and FLTTP are not available for extracellular transfer but they add to the activity of the compartment. FLTMP is the product of TK1 activity on FLT in the first few minutes, after which there is progressive production of FLTTP (17). A component of this metabolism, the dephosphorylation of FLTMP by dNT, is defined in our model as k4. It is important since FLT resulting from dephosphorylation can be lost from the tissue region through a metabolic pathway independent of the formation of FLTMP. This loss has been estimated in vitro by monitoring the egress of activity from cells labeled with FLT nucleotides over time (17). Thus, based on the in vitro data, k4 was given the initial value of 0.02 min1 (range, 0.0010.200 min1). Typical values and ranges of K1, Vd, k3, and k4 that were used as the starting condition in the model optimization process are summarized in Table 2.
Model Analysis
Model Input Function.
The FLT plasma clearance curve serves as the input to the model. For use in model analysis and simulations, we estimated a typical FLT plasma clearance curve from patient results. Nine patients with arterial catheters each had 30 arterial blood samples assessed for radioactivity throughout the 120 min of imaging in the following sequence: eight 15-s, two 30-s, three 1-min, three 3-min, five 5- min, and nine 10-min samples. For each blood sample, 0.1 mL plasma was assayed for radioactivity. As reported in our companion paper, there were no statistical differences between blood and plasma for the 120-min imaging period, but the plasma curve showed less sampling variability. The 9 arterial plasma activity curves were peak normalized and time shifted to align the peak of the curves. No curve was shifted more than ±0.5 min. The average activity at each time was generated by interpolating the activity of the nearest time points. The average total blood time-activity curve was then scaled to a 185-MBq injection of FLT used in further model simulations.
Because the catabolism of FLT in the liver produces a metabolite that contributes to blood activity, the total measured activity in blood was not entirely labeled FLT during the imaging period. To assess the proportion of FLT and labeled metabolites in blood over time, an aliquot of plasma from several samples was assayed for FLT content. As previously described (10), the individual curves for FLT fraction versus time after injection were empirically fitted to a monoexponential function. An average curve for FLT blood clearance, accounting for metabolites, was generated to provide the FLT input function (CpFLT) used for further modeling analysis.
The proposed model was evaluated to determine the extent to which the information obtained from a typical imaging study is sufficient to produce a unique solution with identifiable parameters. The model was characterized with respect to parameter sensitivity, identifiability, and susceptibility to noise to establish the most reliable approach for parameter estimation.
Sensitivity Analysis.
Sensitivity analysis measures the degree to which a change in an individual input parameter results in a change in the output of a model. Sensitivity functions for each parameter in the model were generated by determining the change in model output resulting from a 1% deviation in the parameter using the following formula (35):
![]() | (Eq. 6) |
TAC(t) is the difference in the model output of the time-activity curves resulting from a ±1% change in the parameter k (
k
/k
= 0.02) at time t, TAC(t) is the activity in the tissue (time-activity curve) predicted by the model at time t. The starting parameter values given in Table 2 are based on our initial analysis of human studies (10). The parameter values used for sensitivity analysis produced a simulated curve representative of an average uptake profile for 3 different representative tissues: marrow, tumor, and muscle. Transfer of FLT into the tissue, washout of unmetabolized FLT, phosphorylation by TK1 to FLTMP, and dephosphorylation of FLTMP to FLT all influence the curve and each has its maximal influence on the time-activity curve at different times. Some changes are opposite in sign to others, suggesting some capacity for the model to estimate each parameter. The required length of the experimental data collection and analysis may be indicated from the sensitivity functions. When parameter sensitivity to the model solution is not changing for a long duration near the end of the scan, the imaging time may be shortened without sacrificing information.
If each of the sensitivity functions for the parameters are sufficiently different over time, then the parameters are identifiable. To test the ability of the model to estimate parameters independently, we used identifiability analysis of the sensitivity curves to generate a sensitivity matrix as previously described (18,36,37). Each element in the sensitivity matrix represents the covariance of any 2 parameters to produce a change in the output of the model. The sensitivity matrix was inverted and normalized to give the correlation matrix with values that range from 1 to 1. Pairs of parameters within the correlation matrix elements close to ±1 will have estimates that covary considerably and are therefore difficult to estimate independently.
Monte Carlo Error Analysis.
The expected SDs of the rate constants and of the macroparameters (Vd, KFLT) were derived using a Monte Carlo technique (38) and simulated tissue time-activity curves. Poisson noise was added to each simulated time-activity curve based on a level derived from the total counts and the duration of each time frame (31). Briefly, each tissue curve was re-sampled based on the noise characteristics at each time frame, and the optimization process was repeated 200 times. The SD of the tabulated parameter estimations from each run was then determined. The bias and coefficient of variation (COV) of each parameter were computed for each tissue type, which gives an approximation of the parameter error.
Model Accuracy.
To investigate the ability to estimate key parameters of FLT metabolism accurately across the entire expected clinical range, we generated a set of noisy simulated tissue curves. The simulated PET time-activity curves were created from randomly selected parameter values inside the expected ranges typical of a tumor time-activity curve from our lung cancer patients. Each parameter was randomly sampled over a range of values observed in a study of 17 patients (10). For this analysis, 200 simulated tissue curves at high temporal resolution (0.02-min sampling) with known parameters were generated and rebinned to yield activity curves similar to our PET time frames (eight 15-s, four 30-s, six 60-s, two 5-min, and ten 10-min bins). This resulted in simulated curves with flux constants ranging between 0.005 and 0.200 mL/min/g, a range similar to that seen in our initial analysis of patient studies. Poisson noise was then added to the simulated PET curves at levels typical for clinical studies for a 185-MBq injection and a 70-kg subject. The error distribution was based on an estimation of variance in the counts using the activity and duration of each time frame in the simulated PET sequence (31,37). This approach produced a 1%3% COV for a simulated 10-min imaging frame at 60 min, which was representative of our PET measurements in patients with lung tumors. The curves were optimized using the initial starting parameters in Table 2. The true (simulated input) and estimated parameter values were compared by examining the correlation coefficient for estimated versus true parameter, the percentage bias (the difference in the measurements for estimated vs. true parameter), and the precision (SEE for estimated vs. true parameter divided by the mean true parameter, SEE/mean), which is an approximation of the associated error for each parameter estimate. Several model conditions were examined, including the 4P model with 120 min of simulated data and a reduced dataset of 60 min. An earlier report implemented a model of FLT metabolism based on a 3-parameter model (3P) excluding k4 (k4 = 0) in the optimization process with 60 min of data (8). We examined these alternative model conditions for bias, precision, and associated error (SEE/mean).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The model shows good sensitivity and precision for estimating K1 and KFLT but has difficulty reliably estimating k3 and Vd (k2) independently, similar to what we found for 11C-thymidine (18,37). Monte Carlo error analysis for marrow shows a large COV for parameters other than K1 and KFLT, including k4.
Simulations using a noise level typical of untreated lung cancer in human PET FLT time-activity curves showed accurate estimates of K1 and KFLT with little bias and SEs of <15% over a realistic range of values. In a prior report assessing the change in cellular proliferation, treatment responders had a 40%50% decrease in the analogous parameter for 11C-thymidine, KTDR (3). This suggests that applying KFLT as a measure of proliferation has the potential for sufficient precision to be useful in measuring therapy response, where anticipated changes in KFLT will likely exceed the errors in the parameter estimate.
The model had difficulty estimating loss of labeled nucleotides (k4), even with 120 min of data, which is similar to FDG (39). Although estimates of k4 are not as accurate as other model parameters, it accounts for the effect of label lost from the tissue for the best possible determination of KFLT. The ability to estimate small k4 values with 60 min of data is limited, with errors of >80%.
Estimation of k4 in the 4P model may be important for interpreting FLT metabolism (KFLT)especially in tumors, where therapy may have an effect on nucleotide reutilization. The simulations suggested that longer imaging times may be helpful to reduce the bias and associated error (SEE/mean) in estimation of KFLT. Eliminating k4 in lung tumor simulations leads to a bias in the KFLT estimate (underestimated by approximately 30%). This bias could mask a biologic response to therapy, as estimated by a 3P model or a nonmodeling approach, emphasizing the need to assess possible changes in the k4 estimate over the course of treatment. For the 60-min dataset, the k4 estimate had a large error; however, the resulting KFLT estimate had errors of approximately 10% with little bias and was half of the KFLT error using a 3P model.
Marrow showed less bias in the KFLT estimate than in tumor simulations. It is important to note that simulations based on data from untreated lung cancer patients may not represent other pathologic disease, such as cancer of the breast or colon, where the uptake behavior may be closer to marrow than lung cancer. Examination of full-length datasets in a variety of cancers will be necessary to guide future studies using FLT and determine the trade-off between the improved accuracy in parameter estimation with longer imaging times versus the practical need to minimize the time of clinical imaging studies. Our results, using model parameters appropriate for lung cancer, suggest that shortening the imaging time from 90 to 60 min may result in increasing the error in the flux estimate from approximately 5% to approximately 10%, and this would decrease the potential for FLT to measure response to cytostatic therapy.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Mark Muzi, MS, University of Washington Medical Center, Box 356465, 1959 N.E. Pacific St., Seattle, WA 98195-6465.
E-mail: muzi{at}u.washington.edu
| REFERENCES |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
Y. Menda, L. L. Boles Ponto, K. J. Dornfeld, T. J. Tewson, G. L. Watkins, M. K. Schultz, J. J. Sunderland, M. M. Graham, and J. M. Buatti Kinetic Analysis of 3'-Deoxy-3'-18F-Fluorothymidine (18F-FLT) in Head and Neck Cancer Patients Before and Early After Initiation of Chemoradiation Therapy J. Nucl. Med., July 1, 2009; 50(7): 1028 - 1035. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P.S. Dunphy and J. S. Lewis Radiopharmaceuticals in Preclinical and Clinical Development for Monitoring of Therapy with PET J. Nucl. Med., May 1, 2009; 50(Suppl_1): 106S - 121S. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. J. Kim, J. S. Lee, K. C. Im, S.-Y. Kim, S.-A. Park, S. J. Lee, S. J. Oh, D. S. Lee, and D. H. Moon Kinetic Modeling of 3'-Deoxy-3'-18F-Fluorothymidine for Quantitative Cell Proliferation Imaging in Subcutaneous Tumor Models in Mice J. Nucl. Med., December 1, 2008; 49(12): 2057 - 2066. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Mankoff and J. F. Eary Proliferation Imaging to Measure Early Cancer Response to Targeted Therapy Clin. Cancer Res., November 15, 2008; 14(22): 7159 - 7160. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Herrmann, F. Eckel, S. Schmidt, K. Scheidhauer, B. J. Krause, J. Kleeff, T. Schuster, H.-J. Wester, H. Friess, R. M. Schmid, et al. In Vivo Characterization of Proliferation for Discriminating Cancer from Pancreatic Pseudotumors J. Nucl. Med., September 1, 2008; 49(9): 1437 - 1444. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. F. Shields, J. M. Lawhorn-Crews, D. A. Briston, S. Zalzala, S. Gadgeel, K. A. Douglas, T. J. Mangner, L. K. Heilbrun, and O. Muzik Analysis and Reproducibility of 3'-Deoxy-3'-[18F]Fluorothymidine Positron Emission Tomography Imaging in Patients with Non-Small Cell Lung Cancer Clin. Cancer Res., July 15, 2008; 14(14): 4463 - 4468. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Bading and A. F. Shields Imaging of Cell Proliferation: Status and Prospects J. Nucl. Med., June 1, 2008; 49(Suppl_2): 64S - 80S. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Buck, K. Herrmann, C. M. z. Buschenfelde, M. E. Juweid, M. Bischoff, G. Glatting, G. Weirich, P. Moller, H.-J. Wester, K. Scheidhauer, et al. Imaging Bone and Soft Tissue Tumors with the Proliferation Marker [18F]Fluorodeoxythymidine Clin. Cancer Res., May 15, 2008; 14(10): 2970 - 2977. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Bradbury, D. Hambardzumyan, P. B. Zanzonico, J. Schwartz, S. Cai, E. M. Burnazi, V. Longo, S. M. Larson, and E. C. Holland Dynamic Small-Animal PET Imaging of Tumor Proliferation with 3'-Deoxy-3'-18F-Fluorothymidine in a Genetically Engineered Mouse Model of High-Grade Gliomas J. Nucl. Med., March 1, 2008; 49(3): 422 - 429. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Mankoff, J. F. Eary, J. M. Link, M. Muzi, J. G. Rajendran, A. M. Spence, and K. A. Krohn Tumor-Specific Positron Emission Tomography Imaging in Patients: [18F] Fluorodeoxyglucose and Beyond Clin. Cancer Res., June 15, 2007; 13(12): 3460 - 3469. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Muzi, A. M. Spence, F. O'Sullivan, D. A. Mankoff, J. M. Wells, J. R. Grierson, J. M. Link, and K. A. Krohn Kinetic Analysis of 3'-Deoxy-3'-18F-Fluorothymidine in Patients with Gliomas J. Nucl. Med., October 1, 2006; 47(10): 1612 - 1621. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. J. Kelloff, K. A. Krohn, S. M. Larson, R. Weissleder, D. A. Mankoff, J. M. Hoffman, J. M. Link, K. Z. Guyton, W. C. Eckelman, H. I. Scher, et al. The Progress and Promise of Molecular Imaging Probes in Oncologic Drug Development Clin. Cancer Res., November 15, 2005; 11(22): 7967 - 7985. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Tseng, M. Dandekar, M. Subbarayan, Z. Cheng, J. M. Park, S. Louie, and S. S. Gambhir Reproducibility of 3'-Deoxy-3'-18F-Fluorothymidine MicroPET Studies in Tumor Xenografts in Mice J. Nucl. Med., November 1, 2005; 46(11): 1851 - 1857. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Muzi, H. Vesselle, J. R. Grierson, D. A. Mankoff, R. A. Schmidt, L. Peterson, J. M. Wells, and K. A. Krohn Kinetic Analysis of 3'-Deoxy-3'-Fluorothymidine PET Studies: Validation Studies in Patients with Lung Cancer J. Nucl. Med., February 1, 2005; 46(2): 274 - 282. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |