Outline
Comptes Rendus

Medical sciences / Sciences médicales
A method to quantify at late imaging a release rate of 18F-FDG in tissues
Comptes Rendus. Biologies, Volume 328 (2005) no. 8, pp. 767-772.

Abstracts

This theoretical work shows that the rate constant for the 18F-FDG release in tissues can be assessed without needing any arterial blood sampling. The method requires that the clearance of 18F-FDG from plasma has occurred, whereas 18F-FDG is still present in the tissue. This condition can be met dating from 3 h after 18F-FDG injection, when hydration and/or phlorizin injection are applied after the routine static acquisition. The release rate constant can be obtained from a graphical analysis performed at the later decreasing phase of the tissue tracer activity. A two-compartment and a three-compartment model are developed, both in accordance with one another.

Le principe d'une quantification du taux de libération du 18F-FDG par un tissu ne nécessitant pas de prélèvement de sang artériel est présenté. La méthode est valide lorsque la disparition du traceur du plasma est intervenue, alors que le traceur est toujours présent dans le tissu. Cette condition peut être remplie 3 h après l'injection du traceur, avec une hydratation et/ou l'injection de phlorizine pratiquées après l'acquisition statique de routine. Le taux de libération peut être alors obtenu par une analyse graphique de la phase tardive de la courbe d'activité du 18F-FDG dans le tissu. Deux modèles sont proposés, à deux puis à trois compartiments, qui sont en accord.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crvi.2005.06.004
Keywords: PET, 18F-FDG release, Kinetic modelling, k4
Mots-clés : TEP, Relaxation du 18F-FDG, Modèle cinétique, $ {k}_{4}$

Éric Laffon 1, 2; Michèle Allard 1; Roger Marthan 2; Dominique Ducassou 1

1 Service de médecine nucléaire, hôpital du Haut-Lévêque, centre hospitalo-universitaire de Bordeaux, av. de Magellan, 33600 Pessac, France
2 Laboratoire de physiologie cellulaire respiratoire, Inserm E356, université Bordeaux-2, 33076 Bordeaux, France
@article{CRBIOL_2005__328_8_767_0,
     author = {\'Eric Laffon and Mich\`ele Allard and Roger Marthan and Dominique Ducassou},
     title = {A method to quantify at late imaging a release rate of {\protect\textsuperscript{18}F-FDG} in tissues},
     journal = {Comptes Rendus. Biologies},
     pages = {767--772},
     publisher = {Elsevier},
     volume = {328},
     number = {8},
     year = {2005},
     doi = {10.1016/j.crvi.2005.06.004},
     language = {en},
}
TY  - JOUR
AU  - Éric Laffon
AU  - Michèle Allard
AU  - Roger Marthan
AU  - Dominique Ducassou
TI  - A method to quantify at late imaging a release rate of 18F-FDG in tissues
JO  - Comptes Rendus. Biologies
PY  - 2005
SP  - 767
EP  - 772
VL  - 328
IS  - 8
PB  - Elsevier
DO  - 10.1016/j.crvi.2005.06.004
LA  - en
ID  - CRBIOL_2005__328_8_767_0
ER  - 
%0 Journal Article
%A Éric Laffon
%A Michèle Allard
%A Roger Marthan
%A Dominique Ducassou
%T A method to quantify at late imaging a release rate of 18F-FDG in tissues
%J Comptes Rendus. Biologies
%D 2005
%P 767-772
%V 328
%N 8
%I Elsevier
%R 10.1016/j.crvi.2005.06.004
%G en
%F CRBIOL_2005__328_8_767_0
Éric Laffon; Michèle Allard; Roger Marthan; Dominique Ducassou. A method to quantify at late imaging a release rate of 18F-FDG in tissues. Comptes Rendus. Biologies, Volume 328 (2005) no. 8, pp. 767-772. doi : 10.1016/j.crvi.2005.06.004. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2005.06.004/

Version originale du texte intégral

1 Introduction

Assuming that the 2-[18F]fluoro-2-deoxy-d-glucose (18F-FDG) is trapped in cancer cells in an irreversible manner, PET-scan allows quantification of the metabolic rate of glucose in tumours, by several methods that aimed at quantifying the 18F-FDG uptake rate [1–5]. Is the irreversibility of the 18F-FDG uptake an oversimplification? Usually, most biological reactions do not occur in a completely irreversible manner, which has been clearly established in neuroreceptor quantification models [6]. Moreover, in the particular case of 18F-FDG PET imaging, the tracer metabolism by hepatic tumours challenges this notion of irreversibility [7]. A first quantification method of the 18F-FDG release rate constant has been proposed by Phelps et al. [8], as an extension of Sokoloff's model [1]. Then, Patlak has proposed a graphical method for reversible tracers [2], as an extension of his previous method for irreversible ones. These latter methods both require a dynamic acquisition, which is achieved during the first hour after the tracer injection, and an invasive blood tracer sampling.

The aim of this theoretical work is to propose a graphical method allowing for the assessment of the 18F-FDG release in tissues, on certain conditions, but without needing any blood tracer sampling. It is an extension of a recently published method that aimed at quantifying the rate constant for an irreversible tracer uptake [5]. It requires the acquisition of a dynamic PET scan at the later decreasing phase of the tissue tracer activity. Both a two-compartment and a three-compartment model are developed, and simulated data are provided.

2 Two-compartment model

Let us consider the 18F-FDG uptake from blood to an intracellular compartment (Fig. 1a). The 18F-FDG uptake rate constant is K, and the 18F-FDG release rate constant is kB. At the steady state, the rate of trapped tracer per tissue volume unit dCT/dt is:

dCT/dt=KCP(t)kBCT(t)λCT(t)(1)
where CP(t) represents the tracer concentration in the plasma, and λ is the 18F physical decay constant (λ=ln2/110min−1).

Fig. 1

(a) The two-compartment model roughly considers compartments (A) and (B), which involve blood free 18F-FDG and 18F-FDG trapped in intracellular cytoplasm, respectively. The 18F-FDG uptake and release rate constants, K and kB, respectively, are the transport rate constants between the two compartments. (b) In comparison with the two-compartment model, the three-compartment model introduces an intermediate reversible compartment (R), which considers free 18F-FDG in interstitial fluid and intracellular cytoplasm. The rate constants k1 and k2 account for forward and reversed transport between compartments (A) and (R), respectively. The rate constants k3 and k4 account for forward and reversed transport between compartments (R) and (B), i.e. for 18F-FDG phosphorylation by hexokinase and 18F-FDG-6-PO4 dephosphorylation by G-6-phosphatase, respectively.

In patients with a normal renal and hepatic function, the time duration required to reach the steady state is estimated to last about 5 min [3], and the amount of the tracer that accumulates in the tissue during this time duration is assumed to be negligible in comparison with that which accumulates during the whole examination. Moreover, for simplicity, let us assume that the 18F-FDG blood time-activity curve (TAC) exponentially decays with a mean constant α, after the peak (see simulated data for a multiexponential decay). In a nearly similar manner to that used in a previous work on 18F-FDG (irreversible) uptake [5], when t=0 is set at the blood tracer peak, the solution of Eq. (1) is:

CT(t)=KCP(t=0)[e(λ+kB)teαt]/(αλkB)(2)

Eq. (2) is very similar to that of a radioactive chain equation, but with a whole decay constant equal to λ+kB. When e(λ+kB)t is much greater than eαt, i.e. when the clearance of 18F-FDG from plasma has occurred whereas 18F-FDG is still present in the tissue, the later phase of the trapped tracer TAC can be approximated by:

λCT(t)λKCP(t=0)e(λ+kB)t/(αλkB)(3)

Eq. (3) shows that a single exponential fitting to the later phase of the 18F-FDG TAC can provide the release rate constant kB, since λ is known, without needing any arterial blood sampling. (If an arterial blood sampling is carried out, a graphical analysis can also provide the 18F-FDG uptake rate constant K [5].) Let us note that, at the later phase of the tracer activity in the tissue: (i) neglecting the part of the tracer amount accumulated in the tissue during the time duration to reach the steady state is justified, since both 18F physical decay and 18F-FDG release have occurred, and (ii) λCT(t) actually corresponds to the measured trapped 18F-FDG activity, since the blood 18F-FDG is negligible (see simulated data). Furthermore, Eq. (2) allows us to predict that the time Tm corresponding to the maximal 18F-FDG activity in the tissue is lower than that estimated when no 18F-FDG release occurs (Tm) [5]:

Tm=ln[α/(λ+kB)]/[α(λ+kB)](4)
Tm=ln[α/λ]/[αλ](5)

3 Three-compartment model

In a three-compartment model (Fig. 1b), the 18F-FDG transport rate constants between the compartments are different from those of a two-compartment model (Fig. 1a). However, relationships can be established [2]:

K=k1k3/(k2+k3)(6)
kB=k4k2/(k2+k3)(7)

Following Patlak [2], let us suppose that k4 is much lower than the other rate constants of the system. Furthermore, let us consider measurements are started up at the later phase of the 18F-FDG activity in the tissue, when the clearance of 18F-FDG from plasma has occurred whereas 18F-FDG is still present in the tissue, in other words, when the tracer uptake from the blood into the tissue is much lower than the tracer release from the tissue. Thereby, when the tracer uptake from the blood into the tissue can be ignored, the rate of change of trapped tracer per tissue volume unit dCT/dt is:

dCT(t)/dt=k3Crev(t)k4CT(t)λCT(t)(8)
where Crev is the free 18F-FDG in the tissue. The rate of change of the free tracer per tissue volume unit dCrev(t)/dt is:
dCrev(t)/dt=k4CT(t)k3Crev(t)k2Crev(t)λCrev(t)(9)
The whole solutions of Eqs. (8) and (9) are, respectively:
CT(t)=CT0e(λ+k4)t+k3e(λ+k4)t]0tCrev(τ)e(λ+k4)τdτ(10)
Crev(t)=CR0e(λ+k2+k3)t+k4e(λ+k2+k3)t]0tCT(τ)e(λ+k2+k3)τdτ(11)
where CT0 and CR0 are constants representing trapped and free 18F-FDG concentration in the tissue at the beginning of the measurement, respectively. This measurement start, i.e. t=0 in the three-compartment model, occurs at the steady state of the later decreasing phase of the tissue tracer TAC, i.e. when the tracer uptake from the blood into the tissue is much lower than the tracer release from the tissue, hence when Crev mainly depends on the release of trapped tracer. Consequently, (a) in Eq. (10) the second term of the right-hand side is much lower than the first one (also see below), (b) in Eq. (11) (representing the whole solutions of Eq. (9)), the first term of the right-hand side is not relevant and can be cancelled. Thus, introducing Eq. (10) in Eq. (11) yields the following solution for free 18F-FDG in the tissue, when the second order term is neglected:
Crev(t)k4e(λ+k2+k3)t0tCT0e(λ+k4)τ+(λ+k2+k3)τdτ(12)
Then, integrating Eq. (12) leads to:
Crev(t)[k4/(k2+k3k4)]CT0[e(λ+k4)te(λ+k2+k3)t](13)
hence, since k4k2+k3:
Crev(t)[k4/(k2+k3)]CT0e(λ+k4)t(14)

This result justifies that in Eq. (10) the second term of the right-hand side can be considered much lower than the first one. Furthermore, introducing Eq. (14) into Eq. (10) and rearranging yields the following expression for the trapped tracer activity:

λCT(t)λCT0e(λ+k4)t(1+kRt)(15)
kR=k4k3/(k2+k3)(16)
where kR can be considered a 18F-FDG ‘retrapping’ rate constant, by comparison with Eq. (6). Eq. (15) also approximates the measured tissue activity, since the part of the free tracer activity can be considered to be negligible in comparison with that of the trapped tracer.

Furthermore, Eq. (17) hereafter shows that the transport rate constants available from Eqs. (15) and (16), k4 and kR in a three-compartment model, are in agreement with that available from Eq. (3), kB in a two-compartment model:

e(λ+kB)t=e(λ+k4kR)te(λ+k4)t[1+kRt](17)

4 Simulated data

The proposed two-compartment model assumes, for simplicity, that the 18F-FDG blood TAC decays with a mean constant α. However, it has been shown that the 18F-FDG plasma clearance was multi-exponential, i.e. that the input function (IF) was three-exponential [9]. When such an IF is taken into account, Eq. (2) becomes:

CT(t)=i=13KAi(e(λ+kB)teαit)/(αiλkB)(18)

Fig. 2 shows simulated data from Eq. (18) and references [9,10], for different values of kB and α3. At the later phase of the 18F-FDG TAC, when exp[(λ+kB)t] is much greater than exp[αit], a single exponential fitting can provide the 18F-FDG release rate constant kB, since λ is known, without needing any arterial blood sampling. One hour after injection, in comparison with λ+kB=0.0113min−1 [10], the two first IF exponential constants satisfy this condition, α1=1.568min−1 and α2=0.1438min−1 [9], but not α3=0.0109min−1 [9].

Fig. 2

Simulation of the tissue 18F-FDG TAC in a two-compartment model (Eq. (18)), for different values of α3 and kB [9–11]: (1) α3=0.0109, kB=0min−1; (2) α3=0.0109, kB=0.0050min−1; (3) 2.25×α3=0.0245, kB=0.0050min−1. The axis origin is set at the blood tracer TAC peak (not shown).

The greater the constant α3, the faster the tracer plasma clearance 1 h after injection, the earlier the implementation of the method. This feature highlights the major part of hydration and/or phlorizin injection [11,12], which should be made in practice from the end of the routine static acquisition to the beginning of the dynamic acquisition, in order to increase α3. It has been shown indeed that phlorizin injection resulted in an increase by a factor of 2.25 in the plasma clearance rate of 18F-FDG in rats, and that urinary excretion of 18F-FDG in phlorizin-treated rats was very near that of hydrated rats [11]. Consequently, assuming that at least the condition:

exp[(λ+kB)t]=10×exp[α3t](19)
is required, when α3 is increased to a value of 2.25×0.0109=0.0245min−1, the time landmark tl to apply the proposed method is about 3 h after 18F-FDG injection.

From data involved in curve 3 of Fig. 2 (kB=0.0050 and α3=0.0245min−1, i.e. with hydration and/or phlorizin injection), single exponential fittings to the 180–240, 240–300, 300–360, and 360–420-min time ranges lead to estimates of kB equal to 0.0042, 0.0046, 0.0048 and 0.0049 min−1, respectively, which should be compared with the value of 0.0050 min−1. Note that curve 2 in Fig. 2 (kB=0.0050 and α3=0.0109min−1, i.e. without hydration and/or phlorizin injection), leads in the 360–420-min time range to an (under)estimate of 0.0023 min−1.

5 Discussion

Without needing any blood tracer sampling, when hydration and/or phlorizin injection are applied after the routine static PET acquisition, a two-compartment model analysis allows us to quantify the 18F-FDG release rate constant kB from a dynamic acquisition occurring at about 3 h after 18F-FDG injection. The later the dynamic acquisition, the more accurate the kB estimate, but, the greater the 18F-FDG release from the tissue, the more underestimated the value of kB. Moreover, this theoretical work presents a three-compartment model analysis, which supports the two-compartment one, and allows us to differentiate between the rate constants for 18F-FDG dephosphorylation and rephosphorylation, k4 and kR, respectively. The principle of the present model analysis could be applied to other 18F-labelled tracers than the 18F-FDG, at the later decreasing phase of the tissue TAC, when the clearance of the tracer from plasma has occurred whereas it is still present in the tissue.

From Eq. (18) and literature data [9–11], Fig. 2 shows quite different 18F-FDG tissue TACs, related to different values of kB and of the late IF exponential constant α3. To the very best of our knowledge, the later decreasing phase of 18F-FDG tissue TACs and the significant parts of kB and α3 have not been studied so far. It is suggested that, given the limited understanding of tumour metabolism, in addition to the quantification of the 18F-FDG uptake, the quantification of the 18F-FDG release could be investigated to characterize tumours, to differentiate malignant processes from infectious or inflammatory processes [13], and to assess the effects of treatment [7]. Furthermore, Patlak's graphical analysis for assessing 18F-FDG release rate constant requires an invasive arterial blood sampling and the acquisition of a dynamic PET scan that are performed during the first hour after the tracer injection [2–4]. In contrast, a similar assessment that can be obtained from the proposed graphical analysis appears to be easier to implement in a busy clinical environment, since it is non-invasive and requires a dynamic acquisition that occurs far from the static acquisition.

In conclusion, the principle of a graphical analysis that aimed at quantifying the rate constant for 18F-FDG release in tissues is presented, which is an extension of a recently published method for assessing the tracer uptake rate [5]. This whole method emphasizes the interest of the later decreasing phase of the 18F-FDG activity in tissues.

Acknowledgement

The authors wish to thank Schering-CISbiointernational for helpful discussion, and H. Dupouy and I. Hessling for technical assistance. They also thank an anonymous reviewer for his comments on the manuscript.


References

[1] L. Sokoloff; M. Reivich; C. Kennedy; M.H. Des Rosiers; C.S. Patlak; K.D. Pettigrew; O. Sakurada; M. Shinohara The |14C]deoxyglucose method for the measurement of local cerebral glucose utilization: theory, procedure, and normal values in the conscious and anesthetized albino rat, J. Neurochem., Volume 28 (1977), pp. 897-916

[2] C.S. Patlak; R.G. Blasberg Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data: generalizations, J. Cereb. Blood Flow Metab., Volume 5 (1985), pp. 584-590

[3] N. Sadato; T. Tsuchida; S. Nakamura; A. Waki; H. Uematsu; N. Takahashi; N. Hayashi; Y. Yonekura; Y. Ishii Non-invasive estimation of the net influx constant using the standardized uptake value for quantification of FDG uptake of tumours, Eur. J. Nucl. Med., Volume 25 (1998), pp. 59-564

[4] N.M. Freedman; S.K. Sundaram; K. Kurdziel; J.A. Carrasquillo; M. Whatley; J.M. Carson; D. Sellers; S.K. Libutti; J.C. Yang; S.L. Bacharach Comparison of SUV and Patlak slope for monitoring of cancer therapy using serial PET scans, Eur. J. Nucl. Med., Volume 30 (2003), pp. 46-53

[5] E. Laffon; M. Allard; R. Marthan; D. Ducassou A method to quantify the uptake rate of 2-[18F]fluoro-2-deoxy-d-glucose in tissues, Nucl. Med. Commun., Volume 25 (2004), pp. 851-854

[6] M. Ichise; J.-H. Meyer; Y. Yonekura An introduction to PET and SPECT neuroreceptor quantification models, J. Nucl. Med., Volume 42 (2001), pp. 755-763

[7] S. Okazumi; K. Isono; K. Enomoto; T. Kituchi; M. Osaki; H. Yamamoto; H. Hayashi; T. Asano; M. Ryu Evaluation of liver tumors using fluorine-18-fluorodeoxyglucose PET: characterization of tumor and assessment of effect of treatment, J. Nucl. Med., Volume 33 (1992), pp. 333-339

[8] M.E. Phelps; S.C. Huang; E.J. Hoffman; C. Selin; L. Sokoloff; D.E. Kuhl Tomographic measurement of local cerebral glucose metabolic rate in humans with (F-18)2-fluoro-2-deoxy-d-glucose: validation of method, Ann. Neurol., Volume 6 (1979), pp. 371-388

[9] S. Eberl; A.R. Anayat; R.R. Fulton; P.K. Hooper; M.J. Fulham Evaluation of two population-based input functions for quantitative neurological FDG PET studies, Eur. J. Nucl. Med., Volume 24 (1997), pp. 299-304

[10] L.G. Strauss; A. Dimitrakopoulou-Strauss; D. Koczan; L. Bernd; U. Haberkorn; V. Ewerbeck; H.-J. Thiesen 18F-FDG kinetics and gene expression in giant cell tumors, J. Nucl. Med., Volume 45 (2004), pp. 1528-1535

[11] J.K. Moran; H.B. Lee; M.D. Blaufox Optimization of urinary FDG excretion during PET imaging, J. Nucl. Med., Volume 40 (1999), pp. 1352-1357

[12] J.R.L. Ehrenkranz; N.G. Lewis; C.R. Kahn; J. Roth Phlorizin: a review, Diabetes Metab. Res. Rev., Volume 21 (2005), pp. 31-38

[13] A. Matthies; M. Hickeson; A. Cuchiara; A. Alavi Dual time point 18F-FDG PET for the evaluation of pulmonary nodules, J. Nucl. Med., Volume 43 (2002), pp. 871-875


Comments - Policy