|
|
||||||||
Clinical Investigations |
1 Department of Nuclear Medicine, University of Mainz, Mainz, Germany
2 Division of Nuclear Medicine, Department of Radiology, Johns Hopkins Medical Institutions, Baltimore, Maryland
3 Department of Psychiatry, University of Mainz, Mainz, Germany
4 Institute for Nuclear Chemistry, University of Mainz, Mainz, Germany
5 PET Center and Center for Functionally Integrative Neuroscience, Aarhus, Denmark
6 Department of Psychiatry and Psychotherapy, RWTH Aachen University, Aachen, Germany
| ABSTRACT |
|---|
|
|
|---|
Key Words: fallypride raclopride desmethoxyfallypride quantification binding potential parametric
| INTRODUCTION |
|---|
|
|
|---|
Diverse methods for the quantification of 11C-RAC binding in living brain have been reported (5,1115). These methods include relatively invasive techniques requiring a metabolite-corrected arterial input function and more convenient techniques using a nonbinding reference tissue or simple binding ratios. The choice of method can influence the trade-off between precision and accuracy in the estimation of binding potential (BP), the parameter of interest, which is proportional to the density of dopamine receptors, reduced to some extent by competition from endogenous dopamine. Parametric mapping of receptor BP is useful for statistical comparisons of different populations and correlation to external covariates such as clinical rating, without a priori hypotheses about the anatomic locus of such correlations (4,16,17). Several of the methods developed for 11C-RAC analysis are suitable for the voxelwise calculation of BP maps, allowing for the undirected search for altered patterns of D2/3 receptor availability (7).
In contrast to 11C-RAC, 18F-DMFP and 18F-FP are relatively new radiopharmaceuticals with distinct pharmacokinetic properties in vivo. There has been no systematic comparison of results of the several possible methods for the voxelwise quantification of the binding of these 3 tracers in human brain. The purpose of this study was therefore to compare the results of different methods suitable for voxelwise parametric mapping of 11C-RAC, 18F-DMFP, and 18F-FP to provide the basis for selection of standard methods for planned clinical investigations. Therefore, a conventional standard 2-tissue-compartment model with metabolite-corrected arterial input was used to estimate the distribution volume (Vd) and BP of the 3 ligands in groups of healthy volunteers (11,18). Several voxelwise mapping procedures were tested, including reference tissue and arterial input methods. Finally, results of the voxelwise mapping were compared with volume-of-interest (VOI) analysis to test for bias.
| MATERIALS AND METHODS |
|---|
|
|
|---|
18F-DMFP
18F-DMFP was synthesized as described in detail previously (9). PET studies were performed on 16 healthy male volunteers with an age range of 2444 y (mean ± SD, 37 ± 9 y). Subjects were positioned in a Siemens ECAT EXACT PET scanner so that the transaxial slices were parallel to the cantomeatal line, and the head was partially immobilized with elastic polystyrene straps. A 20-min-long transmission scan with an external 68Ge ring source was acquired for attenuation correction of the emission recording. A 124-min-long dynamic emission recording in 3-dimensional (3D) mode consisting of 30 frames increasing in duration from 20 to 600 s was then initiated upon intravenous bolus injection of 209 ± 50 MBq 18F-DMFP (specific activity, 134834 GBq/µmol). Images were reconstructed in 47 planes with an interplane spacing of 3.375 mm, a matrix of 128 x 128 voxels, and a final resolution (full width at half maximum [FWHM]) of 5.4 mm. Arterial radioactivity was recorded online during the first 10 min (ALLOGG AB blood sampler; Allogg Mariefred). Thereafter, samples were collected manually at 1-min intervals up to 20 min and then at 10-min intervals until the end of the recordings. After centrifugation, the plasma radioactivity was measured in a well counter cross-calibrated to the tomograph. Portions of plasma collected at 2, 5, 10, 20, 30, 40, 60, 90, and 120 min were analyzed for the fraction of untransformed 18F-DMFP using thin-layer chromatography (TLC) as described previously (9). The continuous unmetabolized radiotracer input function was calculated by exponential fitting of the measured fractions (19).
18F-FP
Chemistry.
18F-FP was synthesized by a novel high-yield modification of the method for the synthesis of 18F-DMFP (9). In brief, the tosylated precursor ((S)-N-[(1-allyl)-2-pyrrolidinyl)methyl]-5-(3-toluenesulfonyloxy-propyl)-2-methoxybenzamide (5 mg, 10 µmol) was dissolved in 1 mL acetonitrile, treated for 5 min at 65°C with potassium carbonate (5 mg, 36 µmol), and subsequently transferred into a 5-mL vial containing 18F-fluoride using the method of Hamacher et al. (20). The reaction mixture was heated for 20 min at 85°C, diluted with 1 mL phosphoric acid (10%), and separated using high-performance liquid chromatography (HPLC) (250 x 10 mm, RP8; CH3CN: 0.25 mol/L ammonium acetate buffer + 5 mL acetic acid/L, 30:70; 5 mL/min). The fraction containing 18F-FP was isolated, diluted with 0.15 mol/L disodium hydrogen phosphate buffer, and adsorbed on a C18 cartridge to remove the HPLC solvent. The column was washed with 2 mL water and the product was eluted with 1 mL ethanol. The eluant was diluted with 9 mL of an isotonic sodium chloride solution and sterilized by filtration (0.22 µm). Quality control before injection included determination of the chemical and radiochemical purity, specific activity, pH, and absence of pyrogens.
Scanning.
Ten healthy male volunteers with an age range of 2341 y (mean ± SD, 32 ± 9.5 y) were recruited for the study of 18F-FP kinetics. After a 20-min transmission scan, a 240-min-long dynamic emission recording in 3D mode consisting of 45 frames increasing in duration from 20 to 600 s was initiated with the ECAT EXACT tomograph upon intravenous bolus injection of 184 ± 43 MBq 18F-FP (specific activity, 42258 GBq/µmol). Arterial blood sampling was as for 18F-DMFP, with additional blood samples collected at 30-min intervals between 120 and 240 min. Metabolite correction was performed as for 18F-DMFP.
11C-RAC
Six healthy volunteers (4 males, 2 females) with an age range of 1951 y (mean ± SD, 33 ± 11 y) were recruited for the study. A 90-min-long dynamic emission recording consisting of 20 frames increasing in duration from 30 to 300 s was recorded with the Scanditronics 4096+ tomograph in 2-dimensional (2D) mode after intravenous bolus injection of 745 ± 18 MBq 11C-RAC prepared by conventional methods (specific activity, 35169 GBq/µmol). After reconstruction, the final images had a spatial resolution of 6 x 6 x 6 mm FWHM with a matrix of 128 x 128 voxels and 15 slices. A metabolite-corrected arterial input function was obtained as described previously (21).
Image Processing
All dynamic PET recordings were transferred to a Sun Workstation (Sun Microsystems) for further data analysis. Although great care was taken to minimize head movement, displacement (mainly in z-direction) was noted in some cases. Correction for this movement was accomplished as follows: Assuming no head movement during the first 10 min, the summed emission recording in this interval was manually resliced in the anterior commissureposterior commissure orientation. All subsequent emission frames were realigned to this summed recording to optimize matching of the cerebral contours. To check for the quantity of misregistration, we applied a relatively large VOI in the striatum and the quality of the final alignment was assessed framewise by visual inspection. Furthermore, striatal timeradioactivity curves (timeactivity curve [TAC]) were extracted and evaluated for movement artifacts; if present, realignment was performed iteratively until no further improvement could be obtained. Final alignment was obtained by semiautomatic realignment of 30-min blocks of the emission sequence to the initial 10-min summed image. This was done semiautomatically, with manual adjustment where necessary, using commercially available software (MPI-Tool; ATV GmbH). The final dynamic emission recordings were summed and then automatically coregistered to the subjects individual MRI scans.
VOIs drawn on the coregistered individual MR images were used to extract dynamic TACs from the caudate (circular), putamen (irregular), thalamus (elliptic), and cerebellum (irregular). These VOIs were defined in 3 adjacent transaxial slices passing through the regions of greatest extension of the 4 structures.
The magnitudes of Vd and BP in the VOIs were calculated by fitting several models to the timeradioactivity curves and, except where specified, parametric maps of Vd and BP were calculated from the entire dynamic emission recordings. Finally, the mean magnitudes of Vd and BP in the parametric maps were calculated for each anatomic VOI.
Kinetic Analyses
One- and 2-Tissue-Compartment Model.
A 1-tissue-compartment, 3 parameter model with the plasma volume fixed at 0.05 mL g1 (and using the whole blood radioactivity to calculate the contribution of vascular volume) was first fitted to the cerebellum TAC to determine the magnitude of Vd. Next, a 2-tissue-compartment, 5-parameter model, with the K1/k2 ratio constrained to the cerebellum observation Vd, was fitted to the timeradioactivity curves from the striatum and thalamus VOIs (9,18), and the BP was calculated as the ratio k3/k4that is the ratio of ligand association and dissociation rate constants.
Logan Plot with Arterial Input Function.
The tracer Vd values in the VOIs were calculated by the linear graphical method of Logan et al. (13), and the magnitude of BP was calculated as the distribution volume ratio (DVR) relative to cerebellum, minus one. Voxelwise maps of BP were thus calculated using the commercially available software package PMOD (PMOD Technologies Ltd.) (22). Data recorded during the first 20 min were excluded from the linearizations.
Logan Plot with Reference Region.
The magnitude of the DVR in the VOIs was calculated using the reference tissue linearization of Logan et al. (23), and the magnitude of BP was calculated as the DVR, minus one. Data recorded during the first 30 min were excluded from the linearizations. Voxelwise maps of BP were calculated using PMOD software.
Simplified Reference Region Compartmental Method (SRTM).
The magnitude of BP was calculated in the VOIs using the simplified reference region (SRTM) compartmental method (11), and parametric maps of BP were made using the method of Gunn et al. (12), with the cerebellum serving as the reference region.
Linear Regression with Spatial Constraint Modified SRTM (LRSC).
We have developed a modification of the SRTM for optimal calculation of parametric maps of BP by linear regression with spatial constraint (LRSC) (14). Applications of this method are generally based on an analytic solution of the SRTM equation with parameters estimated by multilinear regression. In this method, based on the same assumptions used to derive the SRTM, a new set of operational equations of integral form is derived with parameters directly estimated by conventional weighted linear regression combined with general ridge regression with spatial constraints (24). The LRSC algorithm is implemented for parametric imaging to reduce the effects of high noise levels in voxelwise timeactivity curves that are typical of dynamic PET data.
Simple Ratio.
The simple ratio of specifically bound radioactivity in a receptor-containing (CRC) and a receptor-free (CRF) region, (CRC CRF)/CRF, approaches the magnitude of BP at equilibrium. Using the VOIs, we calculated this ratio in the thalamus and striatum using the summed frames recorded during the interval 3060 min in the case of 11C-RAC, 6090 min for 18F-DMFP, and 150180 min for 18F-FP.
Statistics
To compare results of the different methods, Pearson correlation coefficients were calculated between the individual magnitudes for each ligand and brain region, using the commercial software package SPSS (SPSS Inc.). Significance was calculated after Fischer transformation.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
The results of correlation analysis of the several methods are presented in Table 3. In general, the several estimates of 11C-RAC BPs correlated less well with each other than did those for 18F-DMFP and 18F-FP. The standard 2-tissue-compartment arterial input model BPs correlated significantly with results from the Logan arterial input method only for the case of 18F-DMFP. However, the 2-tissue-compartment BPs correlated significantly with SRTM BPs for all 3 ligands. Highest correlations for all 3 ligands were between the SRTM and Logan cerebellum input methods.
|
| DISCUSSION |
|---|
|
|
|---|
Whereas the apparent in vivo affinity of 18F-DMFP for D2/3 receptors (15 nmol/L) is close to that of 11C-RAC (26 nmol/L) (27), the relatively novel compound 18F-FP has substantially higher affinity (33 pmol/L) than the other 2 benzamides in the present study (2830). This high affinity provides the basis for the particular sensitivity of 18F-FP for the detection of extrastriatal D2/3 receptors (Fig. 3). However, this same high affinity of 18F-FP imposes a requirement for lengthy emission recordings to obtain an equilibrium binding estimate in the striatum.
In addition to their differing affinities for the dopamine D2/3 receptors, the present 3 benzamides are distinct with respect to their nonspecific binding in the cerebellum. Whereas 11C-RAC has a relatively low Vd in the cerebellum, that of 18F-FP was 2-fold higher and that of 18F-DMFP was 10-fold higher. Given the structural similarities and similar log P values (30) of these benzamides, we cannot readily account for the 10-fold difference in the Vd but suggest that additional factors, such as the (unmeasured) plasma protein binding, may account for the high Vd of 18F-DMFP. The high variance in the calculation of the Vd for 18F-DMFP and 18F-FP in the cerebellum and striatum might conceivably be attributed to imprecision in the calculation of the metabolite-corrected arterial input, especially at late times when the concentration of untransformed tracer in plasma is vanishingly small. However, using TLC, we were able to obtain adequate metabolite correction throughout the emission recordings for all 3 tracers.
In general, high nonspecific binding reduces the signal-to-noise ratio in the calculation of specific binding. Thus, the superior precision of the estimates of 11C-RAC BP in the thalamus and striatum may be attributed in part to its low nonspecific binding. This phenomenon may also contribute to the poor sensitivity of 18F-DMFP for the detection of dopamine receptors in the thalamus, where they are not abundant. However, the pharmacologic specificity of the 11C-RAC binding detected in the thalamus has not been demonstrated in previous displacement studies or pharmacologic challenge studies, although we have reported a partial displacement of N-11C-methylspiperone binding in pig thalamus by methylenedioxymethamphetamine challenge (7). The ratio of BP in the striatum to that in the thalamus is much higher for 18F-FP (12:1) than for 11C-RAC (6:1) irrespective of the method of analysis, suggesting that thalamic binding may be overestimated using the latter tracer.
In the striatum, where dopamine receptors are most abundant, the trade-off between precision and accuracy in the estimation of BP is different than that in the thalamus; the coefficient of variation of the BP in the striatum was somewhat less favorable for 18F-FP than for 18F-DMFP, in spite of the former compounds lower nonspecific binding. Here, the very slow approach to equilibrium binding for 18F-FP may compromise the precision of the estimate of its BP in the striatum, even given the present 4-h-long emission recordings.
Irrespective of the method used for parametric mapping, the BP for 11C-RAC was in the range 34 in normal human striatum, as reported previously by others in VOIs using bolus infusion (11,31,32) or continuous infusion methods (26). The SRTM approach for quantifying 11C-RAC binding was validated by Gunn et al. (12). In the present study, we likewise find substantial agreement between the estimates obtained in parametric maps and in striatal VOIs for the case of SRTM and other models. This suggests that the present use of voxelwise mapping does not bias the estimation of the magnitude of the 11C-RAC BP in the striatum. However, the 2-tissue-compartment model gave higher BP estimates than did the other models, albeit with a penalty in precision in the cases of 11C-RAC and 18F-DMFP. Although a relatively high BP for 11C-RAC based on the 2-compartment model has been noted previously (7), overspecification of this model may compromise the precision of the estimate of BP. Thus, on the basis of the present results, we conclude that SRTM presents optimal properties for the generation of parametric maps with 11C-RAC, based on the high correlation between SRTM results and the gold standard of the 2-tissue-compartment model. It has earlier been argued that the spatial constraint approach yields lower image noise in BP maps and is, therefore, preferable to SRTM for accurate stereotactic normalization (14). In the present study, the spatial constraint approach gave slightly (5%) higher BPs for 11C-RAC. The performance of SRTM is dependent on the noise level of kinetics; as noise levels decrease, results calculated by the LRSC approach those obtained by SRTM and the Logan plot with reference region input. SRTM and the Logan plot can have higher noise-induced bias (usually underestimation), but this bias is minimized by LRSC or BP estimates (14).
In the present study, the cerebral kinetics of 18F-DMFP binding have been evaluated in a large number of healthy volunteers (n = 16). As in an earlier report on only 10 subjects (9), there was considerable agreement between several arterial input and reference tissue methods for the calculation of BP. However, as discussed, the precision of the estimate in the striatum was lower than that for the other 2 ligands tested in the present studya difference that we attribute to the high nonspecific binding of 18F-DMFP. Given the similar in vitro affinities reported for 11C-RAC and 18F-DMFP, it is unclear why 18F-DMFP had a lower BP than 11C-RAC did in the present study of living human brain, confirming earlier studies on rhesus monkeys (8). The BP of 18F-DMFP may be underestimated because of its high nonspecific binding relative to the other tracers, which, we speculate, could be related to difference in plasma protein binding since the log P values are quite similar (30). We have found 18F-DMFP and the SPECT ligand iodobenzamide to be equally suitable for the discrimination of idiopathic and atypical parkinsonian syndromes using a simple ratio method (33). Though, on the basis of the present results, 18F-DMFP may not be superior to 11C-RAC for the estimation of availability of dopamine D2/3 receptors in the striatum, we have argued that it presents a potential practical advantage due to the greater transportability of 18F. Furthermore, 18F-DMFP may present an advantage over 11C-RAC in the single bolus or continuous infusion competition paradigms (34,35). The occurrence of a prolonged equilibrium for 18F-DMFP binding in the striatum compared with 11C-RAC predicts that the effects of pharmacologic challenge on dopamine release in the striatum might sensitively be detected in the course of an 18F-DMFP recording of several hours duration. This advantage may outweigh the higher radiation burden imposed by 18F-DMFP (effective dose equivalent: 11C-RAC, 0.01 mSv/MBq; 18F-DMFP, 0.03 mSv/MBq; 18F-FP, 0.017 mSv/MBq). However, the present findings suggest that 18F-DMFP is particularly unsuited for the detection of extrastriatal dopamine receptors.
Depending on the method used for analysis, in the present study we found 18F-FP to have a BP of 2132 in human striatum, as reported previously by Mukherjee et al. (30). The 2-tissue-compartment model BPs correlated significantly with the Logan reference tissue results and with the 2 SRTM results. Therefore, we conclude that arterial inputs are not necessary for good quantitation of the 18F-FP BP in human striatum. In contrast to 11C-RAC and 18F-DMFP (9,11), it might be that not all criteria for estimation of BP by reference tissue methods are fulfilled by 18F-FP (e.g., k2 > k4). Nevertheless, in our group of subjects, the Logan plots were always linear during the time interval used for analysis, which tends to validate the use of these equilibrium methods. Furthermore, the reference tissue Logan plot has been used for 18F-FP quantitation in several previous publications (e.g., 30).
The present estimates of the BP of 18F-FP in the thalamus (2.02.8) were somewhat lower than those in an earlier report (30). Diverse methodologic differences, especially related to the method for drawing the VOI, or in the selection of the frames to which the linearization was applied may be responsible for this difference.
Simple ratio techniques for receptor quantitation are suited to the requirements of clinical investigations, where patients cannot remain in the tomograph for extended periods. Therefore, we tested the correlation between the simple ratio relative to cerebellum measured during a 30-min static recording and several BP estimates. This correlation was not significant for 11C-RAC in the striatum but was significant in the case of several estimates obtained by voxelwise methods for 18F-FP and 18F-DMFP. This phenomenon may be related to the superior signal available from 18F. However, the ratio method overestimated the BP of 18F-DMFP and 11C-RAC in the thalamus and, to a lesser extent, also in the striatum. Since the time required to reach equilibrium is dependent on receptor density, the simple ratio measured at an optimal time point for striatal binding need not perfectly reflect the BP in regions with lower D2 receptor density (36). Thus, the performance of the simple ratio method is dependent on the chosen time frame. Furthermore, differences in the magnitude of K1 (the unidirectional clearance from vascular space to brain tissue) on the estimation of BP may also contribute to the bias from this simple method (34). Thus, the present ratio approach seems inadequate for assessing dopamine D2/3 receptor availability in brain regions with different receptor densities.
The practical necessity of using different subject groups in the present comparison of 3 different ligands imposes some limitations on the interpretations to be drawn from the study. However, since there was no statistically significant difference in age (P > 0.05) between the groups, a confound due to age can be excluded, especially since the declines in BP with increasing age are, in any event, relatively small (8% per decade) (37,38). Nevertheless, some of the present intergroup differences could reflect individual differences in competition from endogenous dopamine at the time of scanning.
Intrinsic differences in binding parameters between study groups maybe contribute to the higher variability in 18F-DMFP and 18F-FP binding between subjects compared with 11C-RAC. Intrasubject variability (i.e., testretest) has been reported to be 5% for 11C-RAC (37) and 10% for 18F-FP (30), suggesting that the binding of the latter ligand may be intrinsically more variable in a normal population. Alternately, steady-state 18F-FP binding estimates may be more difficult to measure due to the slow approach to its equilibrium. To our knowledge, the present study constitutes the first report on intersubject variability in 18F-DMFP binding, indicating a somewhat higher intrasubject variability than is typical for 11C-RAC. Intrinsic variability may be a factor in selecting a probe for the sensitive detection of differences in clinical populations and under conditions of pharmacologic activation studies.
One limitation of the present study is that 11C-RAC data were acquired in 2D mode using a different tomograph than that used for testing the new ligands, 18F-DMFP and 18F-FP. Nonetheless, the present estimates of 11C-RAC binding are very similar to previous estimates in healthy subjects (15) obtained using a tomograph (ECAT 931) with a performance similar to that of the present instrument used for 18F-DMFP and 18F-FP acquisitions. Since the resolution of the PC4096+ (39) is comparable to that of the ECAT EXACT (40), we conclude that the potential bias introduced by instrumentation is not a major factor in the present comparison of the 3 ligands.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Thomas Siessmeier MD, Department of Nuclear Medicine, University of Mainz, Langenbeckstrasse 1, 55131 Mainz, Germany.
E-mail: siessmeier{at}nuklear.klinik.uni-mainz.de
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
G. Grunder, C. Fellows, H. Janouschek, T. Veselinovic, C. Boy, A. Brocheler, K. M. Kirschbaum, S. Hellmann, K. M. Spreckelmeyer, C. Hiemke, et al. Brain and Plasma Pharmacokinetics of Aripiprazole in Patients With Schizophrenia: An [18F]Fallypride PET Study Am J Psychiatry, August 1, 2008; 165(8): 988 - 995. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Fehr, I. Yakushev, N. Hohmann, H.-G. Buchholz, C. Landvogt, H. Deckers, A. Eberhardt, M. Klager, M. N. Smolka, A. Scheurich, et al. Association of Low Striatal Dopamine D2 Receptor Availability With Nicotine Dependence Similar to That Seen With Other Drugs of Abuse Am J Psychiatry, April 1, 2008; 165(4): 507 - 514. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |