Abstract
It has been demonstrated in vitro that the dopamine D2 receptor has 2 interconvertible affinity states for endogenous dopamine, referred to as the high- and the low-affinity states. 11C-(R)-2-CH3O-N-n-propylnorapomorphine (11C-MNPA) is a new agonist radioligand for in vivo imaging of the high-affinity state of dopamine D2 receptors using PET. In the present study, the kinetics of 11C-MNPA were examined for the first time, to our knowledge, in the human brain and analyzed using quantitative approaches with or without an arterial input function. Methods: A 90-min dynamic PET scan was obtained for 10 healthy men after an intravenous injection of 11C-MNPA. The binding potential (BPND) was calculated using the indirect kinetic method, a kinetic compartment analysis with a metabolite-corrected arterial input function. BPND was also calculated by the simplified reference tissue model (SRTM) and transient equilibrium methods, both with the cerebellum as the reference brain region. The results of the quantitative methods were compared in a cross-validation approach. Results: The highest regional radioactivity was observed in the putamen. BPND values obtained by kinetic analysis were 0.82 ± 0.09, 0.59 ± 0.11, and 0.28 ± 0.06, respectively, in the putamen, caudate, and thalamus. BPND values obtained by the SRTM and transient equilibrium methods were in good agreement with those obtained by the indirect kinetic method (r = 0.98 and r = 0.93, respectively). For all quantification methods, the BPND values based on data acquired from 0 to 60 min were in good agreement with those based on data acquired from 0 to 90 min (r = 0.90–0.99). Conclusion: The regional distribution of 11C-MNPA binding was in good agreement with previous PET studies of dopamine D2 receptors in the human brain using antagonist radioligands. The results support routine use of the SRTM and transient equilibrium methods, that is, methods that do not require an arterial input function and need a scan time of only about 60 min. 11C-MNPA should thus be useful for clinical research on the pathophysiology of neuropsychiatric disorders and estimation of dopamine D2 receptor occupancy by dopaminergic drugs.
The dopaminergic system in the brain plays a significant role in the physiologic regulation of motor functions, cognition, emotion, and personality (1). Alterations of dopaminergic neurotransmission have been implicated in several pathologic conditions, such as schizophrenia (2), Parkinson disease (3), and addiction (4). Dopamine receptors are classified into 5 subtypes, D1, D2, D3, D4, and D5 (5). The dopamine D2 receptor is a main therapeutic target for currently used antipsychotic drugs and has long been suggested to be involved in the pathophysiology of schizophrenia (2,6–11).
Several observations in vitro indicate that the dopamine D2 receptor exists in 2 interconvertible affinity states, distinguished by the affinity of endogenous dopamine and referred to as the high- and low-affinity states. It has further been suggested that the high-affinity state is the functionally active form of the receptor and is thus more relevant for clinical studies (12). Studies of dopamine D2 receptors using PET have almost exclusively been performed with antagonist radioligands, such as 11C-raclopride and 11C-FLB457 (13–15). Antagonist radioligands do, however, bind with equal affinity to the 2 conformational states of the dopamine D2 receptor. To examine the high-affinity state of the dopamine D2 receptor in vitro and in vivo, agonist radioligands such as (−)-N-11C-propyl-norapomorphine (11C-NPA) and 11C-(+)-4-propyl-3,4,4a,5,6,10b-hexahydro-2H-naphtho[1,2-b][1,4]oxazin-9-ol (11C-PHNO) have recently been developed (16,17).
11C-(R)-2-CH3O-N-n-propylnorapomorphine (11C-MNPA) is another new agonist PET radioligand with high affinity and selectivity for the dopamine D2 receptor (inhibitory concentration of 50%, 1.02 nM; inhibition constant 0.17 nM, respectively) (18,19) and was recently characterized in nonhuman primates (20). PET measurements in cynomolgus monkeys showed high uptake in the striatum, with a striatum-to-cerebellum ratio of 2:2. The striatal uptake of 11C-MNPA could be inhibited by the injection of unlabeled raclopride, confirming that the striatal binding is reversible and specific for dopamine D2 receptors. Subsequent applied studies using amphetamine-induced dopamine release showed that 11C-MNPA was more sensitive than 11C-raclopride, thus supporting the indication that 11C-MNPA is a promising radioligand for PET of the high-affinity state of the dopamine D2 receptor in vivo (21). The aim of the present study was to examine the regional distribution and kinetics of 11C-MNPA in the human brain. Ten control subjects were included, and data were analyzed using kinetic compartment analyses with a metabolite-corrected input function and 2 quantitative methods with the cerebellum as a reference brain region.
MATERIALS AND METHODS
Subjects
Ten healthy men (age range, 22–35 y; mean ± SD, 27.7 ± 5.4y) participated in this study. On the basis of their medical history and MRI of the brain, all subjects were free of any somatic, neurologic, or psychiatric disorders, and they had no history of current or previous drug abuse. Written informed consent was obtained from all subjects after the study was completely described. The study was approved by the Ethics and Radiation Safety Committee of the National Institute of Radiological Sciences, Chiba, Japan.
PET Procedure
11C-MNPA was synthesized as described in detail previously (20). In brief, 11C-MNPA was synthesized by the methylation of (R)-(−)-2,10,11-trihydroxy-N-n-propylnoraporphine-acetonide with 11C-methyl triflate and subsequent cleavage of the acetonide protecting group with the addition of hydrochloric acid. An ECAT EXACT HR+ PET system (CTI-Siemens) was used for all measurements. A head-fixation device was used to minimize head movements during data acquisition. A transmission scan for attenuation correction was obtained using a 68Ge–68Ga source. Dynamic PET scans were obtained after a 1-min intravenous slow bolus injection of 11C-MNPA (204.3–232.1 MBq; mean ± SD, 219.3 ± 8.2 MBq). The specific radioactivity of 11C-MNPA was 197.5–335.0 GBq/μmol (261.1 ± 43.3 GBq/μmol) at the time of injection. Brain radioactivity was measured from 0 to 90 min (20 s × 9, 1 min × 5, 2 min × 4, 4 min × 11, and 5 min × 6).
MR images of the brain were acquired with a 1.5-T MRI scanner (Gyroscan NT; Philips). T1-weighted images were obtained at 1-mm slices acquired in 3 dimensions.
Arterial Blood Sampling and Metabolite Analysis
To obtain the arterial input function, a series of arterial blood samples were taken manually from a catheter 32 times (2.5 mL × 22 times for the measurement of radioactivity concentration in whole blood and plasma; 5.0 mL × 10 times for the determination of the percentage of unchanged 11C-MNPA in plasma) during the 90-min PET scan. Each blood sample was centrifuged to obtain plasma and blood cell fractions, and the concentrations of radioactivity in whole blood and plasma were measured.
The percentage of unchanged 11C-MNPA in plasma was determined by high-performance liquid chromatography (HPLC) in 10 of the blood samples. Acetonitrile was added to each plasma sample, and the samples were then centrifuged. The supernatant was subjected to radio-HPLC analysis (column, XBridge Prep C18; Waters) (mobile phase, 48:52 90% acetonitrile:50 mM ammonium acetate). The arterial plasma input function was defined as the radioactivity of plasma multiplied by the percentage of unchanged radioligand.
Regions of Interest (ROIs)
All MR images were coregistered to the PET images using the statistical parametric mapping system (SPM2; Wellcome Trust Centre for Neuroimaging) (22). ROIs were drawn manually on summated PET images with reference to the coregistered MR images. ROIs were defined for the cerebellar cortex, putamen, caudate, and thalamus. Regional radioactivity was calculated for each frame, corrected for decay, and plotted versus time.
Kinetic Compartment Analysis of 11C-MNPA Binding
To describe the kinetics of 11C-MNPA in the brain, the 2-tissue-compartment model with 4 rate constants, K1, k2, k3, and k4, was used. The 3 compartments include CP, the radioactivity concentration of unchanged radioligand in plasma (arterial input function); CND, the radioactivity concentration of nondisplaceable radioligand in the brain (including nonspecifically bound and free radioligand); and CS, the radioactivity concentration of radioligand specifically bound to receptors. The rate constants K1 and k2 represent the influx and efflux rates for radioligand diffusion across the blood–brain barrier. The rate constants k3 and k4 represent radioligand transfer between the compartments for nondisplaceable and specifically bound radioligand.
Calculation of the 11C-MNPA Binding Potential
11C-MNPA binding was expressed by the indirect kinetic, SRTM, and transient equilibrium methods. In these methods, 11C-MNPA bindings were expressed as binding relative to nondisplaceable binding (BPND) (23). The BPND of the radioligand is proportional to the product of the receptor density (Bmax) and reciprocal affinity (1/Kd), and the ratio of Bmax to Kd corresponds to the ratio of k3 to k4, as expressed by the following equation:Eq. 1where fND is the free fraction of radioligand in the nondisplaceable compartment. The cerebellum (cerebellar cortex) was used as a reference brain region because it is a structure with negligible D2 dopamine receptor density (24). The software package PMOD (PMOD Technologies) was used for the indirect kinetic, SRTM, and transient equilibrium quantitative methods.
Indirect Kinetic Method.
In the present cross-validation approach, the indirect kinetic method was used as the standard method (25). BPND was defined as the ratio of k3 to k4 as calculated using the 2-tissue-compartment model. Because the ratio of k3 to k4 is sensitive for noise in the PET data, BPND was calculated using the indirect kinetic method. With the cerebellum as a reference region, BPND can be expressed as:Eq. 2where VT(region) is the total distribution volume (= (K1/k2)(k3/k4 + 1)) in a target region, and VT(cerebellum) is the total distribution volume (= K1/k2) in the cerebellum. The rate constants K1, k2, k3, and k4 in the putamen, caudate, and thalamus were determined by nonlinear curve fitting in a least-squares sense to the regional time–activity curves as described in the literature (15). K1 and k2 values in the cerebellum were also determined by nonlinear curve fitting in a least-squares sense but by using the 1-tissue-compartment model, assuming that the cerebellum has negligible D2 dopamine receptor density (24). To improve the stability of the curve fitting in the nonlinear curve-fitting procedure, the ratio of K1 to k2 was fixed for each subject to the value obtained in the cerebellum by the kinetic analysis with the 1-tissue-compartment model (range of K1/k2, 5.0–7.9 mL/cm3; mean ± SD, 6.5 ± 0.75 mL/cm3). In this analysis, blood volume, which depends on the first-pass extraction fraction of the tracer, was estimated using the radioactivity of whole blood to diminish the influence of tracer remaining in the blood (26).
SRTM Method.
The SRTM method, assuming that both target and reference regions have the same level of nondisplaceable binding, can be used to interpret time–activity curves in the target region as follows (27):Eq. 3where CT(t) is the total radioactivity concentration in a brain region measured by PET, R1 is the ratio of K1 to K1′ (K1, influx rate constant for the brain region; K1′, influx rate constant for the reference region), CR(t) is the radioactivity concentration in the reference region (cerebellum), and * denotes the convolution integral. The parameters R1, k2, and BPND in this model were estimated by the nonlinear curve-fitting procedure.
Simulation Study
A simulation study was performed to estimate errors of BPND calculated by the SRTM and transient equilibrium methods. The assumed values and intervals examined were acquired from the results of the kinetic approach. Regional tissue time–activity curves (0–90 min) were generated according to the 2-tissue-compartment model. We assumed that the value of K1/k2 equaled 6.6 mL/cm3 and that of k4 equaled 0.18 min−1. The tissue time–activity curves were generated with K1 values between 0.20 and 0.60 mL/cm3/min in 6 steps and with k3 values between 0.02 and 0.20 min−1 in 10 steps. A tissue time–activity curve for the cerebellum was generated according to the 1-tissue-compartment model with 2 rate constants, using 0.44 mL/cm3/min for K1 and 0.067 min−1 for k2. The average arterial input function (n = 10), corrected for labeled metabolites, was used to generate the tissue time–activity curves. BPND was then calculated by applying the SRTM and transient equilibrium methods to the generated tissue time–activity curves. The estimated BPND values were compared with the BPND values calculated by the indirect kinetic method.
Variations in K1 values between a brain region and the cerebellum due to differences in cerebral blood flow (CBF) between the 2 regions might affect the BPND calculated by the SRTM and transient equilibrium methods. Furthermore, changes in K1 due to changes in CBF might be caused by neurologic or psychiatric diseases. The K1 value for 11C-MNPA was about 0.44 mL/cm3/min in gray matter. When the CBF value in gray matter is assumed to be 0.50 mL/cm3/min (29), the first-pass extraction fraction of 11C-MNPA is 88%. The capillary permeability–surface area product (PS) value, using this extraction fraction and a K1 value of 0.44 mL/cm3/min, was calculated (30,31). With the PS value of 1.06 mL/cm3/min, the K1 range of 0.20–0.60 mL/cm3/min corresponds to the CBF range of 0.20–0.85 mL/cm3/min (28).
RESULTS
All 10 subjects participated in the study according to the protocol. Representative summated PET images (60–90 min) and T1-weighted MR images are shown in Figure 1, and the corresponding regional time–activity curves are shown in Figure 2. Regional radioactivity was highest in the putamen and lower in the caudate and thalamus. The average percentage of unchanged 11C-MNPA in plasma was 95.1% ± 2.1% at 3 min, decreasing to 25.1% ± 12.0% at 90 min (i.e., at the end of PET data acquisition) (Fig. 3). Other than MNPA, there were no more lipophilic-labeled metabolites in the plasma.
Representative summated PET images 60–90 min after intravenous injection of 11C-MNPA (221 MBq) and corresponding T1-weighted MR images in control human subject. Upper panel shows horizontal section through striatum, and lower panel shows section through cerebellum.
Representative regional time–activity curves after intravenous injection of 11C-MNPA (232 MBq) binding in control human subject. Fitted curves using 2-tissue- and 1-tissue-compartment models are also shown for target regions and cerebellum, respectively.
Average (n = 10) percentage of unchanged 11C-MNPA in plasma vs. time. Bars indicate SD.
After an intravenous injection of 11C-MNPA, total radioactivity in the brain peaked at 6.7 ± 1.2 min (range, 4.5–9.0 min), and the fraction of uptake in the brain was 6.0% ± 1.0% (range, 4.3%−7.3%) of the injected radioactivity.
The blood volume and rate constants for each brain region obtained by conventional nonlinear least-squares fit of the 2-compartment model are shown in Table 1. The BPND values of the putamen, caudate, and thalamus calculated by the 3 different methods are shown in Table 2. Specific binding, as defined by the transient equilibrium method, reached a peak within 60 min in the putamen, caudate, and thalamus (Table 2).
Rate Constants Obtained by Conventional Nonlinear Least-Squares Fit of 2-Tissue-Compartment Model
BPND Values Obtained by Different Methods and Scan Times
BPND values determined by the SRTM method on the basis of data acquired for 90 and 60 min and those determined by the transient equilibrium method were compared with values calculated by the indirect kinetic method. BPND values obtained by the SRTM method were in good agreement with those obtained by the indirect kinetic method with data obtained for 90 and 60 min (Fig. 4), and BPND values obtained by the transient equilibrium method were in good agreement with those obtained by the indirect kinetic method with data for 90 min (Fig. 5). The highest coefficient of correlation was observed between the SRTM and the indirect kinetic methods with data acquired for 90 min (r = 0.98, P < 0.001).
Comparison of BPND values in 3 regions (putamen, caudate, and thalamus) of 10 control subjects calculated by indirect kinetic and SRTM methods on the basis of data acquired over 90 (A) and 60 min (B).
Comparison of BPND values in 3 regions (putamen, caudate, and thalamus) in 10 control subjects calculated by indirect kinetic and transient equilibrium methods.
When BPND values determined by the indirect kinetic and SRTM methods with 60-min data were compared with values determined by the same 2 methods with 90-min data, good agreement was observed (r = 0.99, r = 0.92, P < 0.001) (Fig. 6).
Comparison of BPND values in 3 regions (putamen, caudate, and thalamus) in 10 control subjects calculated by indirect kinetic (A) and SRTM (B) methods on the basis of data acquired over 90 and 60 min.
To estimate the sensitivity of the SRTM and transient equilibrium methods for rate constants (indirect blood flow) over an interval with values lower and higher than average, a simulation study was performed. BPND values determined by the indirect kinetic method with data acquired for 60 and 90 min were compared with BPND values determined by the SRTM and transient equilibrium methods with data acquired for 60 and 90 min from simulated time–activity curves. The error in BPND calculated by the SRTM method with data acquired for 90 min was smallest (−24.8% to 1.5%; mean, −4.3%), and the difference in K1 between the brain region and cerebellum had only a minor effect on BPND. The error in BPND calculated by the transient equilibrium method was smallest when the K1 value was 0.44, but BPND was overestimated when the K1 value was lower than 0.36 and was underestimated when it was higher than 0.52 (Fig. 7).
Comparison of simulated BPND values calculated by indirect kinetic and SRTM methods on the basis of data acquired over 90 (A) and 60 min (B). (C) Comparison of simulated BPND values calculated by indirect kinetic and transient equilibrium methods.
DISCUSSION
Studies using agonist radioligands such as 11C-PHNO to examine the high-affinity state of the dopamine D2 receptor in the human brain have been reported previously (32). Our study describes the first, to our knowledge, PET examination using the agonist radioligand 11C-MNPA to visualize binding to G-protein–coupled receptors in the human brain. After the intravenous injection of 11C-MNPA, radioactivity appeared rapidly in the brain and was washed out in a fashion similar to that previously reported in nonhuman primates (20). Radioactivity was highest in the putamen and slightly lower in the caudate, moderate in the thalamus, and lowest in the cerebellum. This regional distribution is similar to that shown in nonhuman primates with 11C-MNPA (20,21) and is in accordance with the known distribution of dopamine D2 receptors, as demonstrated with antagonist radioligands such as 11C-raclopride in the human brain (33). Finnema also reported blocking data with a dopamine D2 antagonist in nonhuman primates (20). The pretreatment with raclopride, compared with the baseline condition, demonstrated high specific binding of the dopamine D2 receptor by reducing the striatum-to-cerebellum ratio. The striatal BPND values of 11C-MNPA were about one third of those in previous studies with the antagonist radioligand 11C-raclopride (33). The Kd value of 11C-raclopride in the human brain in vivo has been reported to be 9.1 nM (34), and the Kd value of 11C-MNPA in the monkey brain in vivo has been reported to be 11.6 nM (35). Because the Kd values of 11C-MNPA and 11C-raclopride are similar in vitro, the difference in striatal BPND between 11C-MNPA and 11C-raclopride may reflect a difference in the density of available receptors (Bmax) of the 2 radioligands. This interpretation is in line with the view that an agonist radioligand labels only the receptors in the high-affinity state, whereas an antagonist radioligand labels both high- and low-affinity-state dopamine D2 receptors (21,35).
Ginovart et al. reported that 11C-PHNO and 11C-NPA in the cat were more sensitive to amphetamine-induced dopamine release than was 11C-raclopride (36). The observation that 11C-MNPA in nonhuman primates is also more sensitive to amphetamine-induced dopamine release than is 11C-raclopride (21) has been taken as evidence for selective labeling of D2 receptors in the high-affinity state. The relatively low BPND in the present study corroborates this view.
In this study, the indirect kinetic method with arterial blood sampling was used as the gold standard (25). Because arterial blood sampling is invasive, we examined the accuracy of the SRTM and transient equilibrium methods for quantifying 11C-MNPA binding using the cerebellum as the reference brain region. The SRTM and transient equilibrium methods had previously been validated for antagonist radioligands such as 11C-raclopride and 11C-FLB457 (15,25,28). In the present study, BPND of 11C-MNPA obtained by the SRTM method was in good agreement with that obtained by the indirect kinetic method with data acquired for 60 and 90 min. The BPND value obtained by the transient equilibrium method was also in good agreement with the value obtained by the indirect kinetic method with data acquired for 90 min. Thus, it should be possible to use simplified protocols with no arterial blood sampling in applied clinical studies in humans.
In the simulation study, BPND calculated by the SRTM method was in good agreement with that calculated by the indirect kinetic method, although BPND was slightly overestimated. These results demonstrate the validity of the SRTM method for quantitating 11C-MNPA binding also when blood flow and rate constants might be deviant. The present observation is in line with an 11C-FLB457 study showing that the BPND value calculated by the SRTM method was not greatly affected by differences in K1 between the brain regions and the cerebellum (25). Thus, the SRTM method is suitable for quantifying 11C-MNPA binding when using a reference brain region without arterial blood sampling.
BPND calculated by the transient equilibrium method was not in good agreement with that calculated by the indirect kinetic method in the simulation study when the K1 value in the brain region was small. The errors in BPND calculated by the transient equilibrium method were within the range of −15% to +15% when the K1 value was 0.44 and 0.52 mL/cm3/min, corresponding to 0.5–0.65 mL/cm3/min of CBF. Although the transient equilibrium method might not be suitable for determining BPND in patients with low CBF, it is still a useful method for determining BPND without arterial blood sampling.
For clinical research, a short scanning time is preferred. In the present study, the BPND values calculated by the SRTM method with data acquired for 90 min were in good agreement with those obtained with data acquired for 60 min. In the simulation study, the BPND values obtained by the SRTM method with data acquired for 60 min were in good agreement with BPND values obtained by the indirect kinetic method, except with extremely low K1. These results suggest that the SRTM method with data acquired for 60 min is valid for clinical studies in patients with neuropsychiatric disorders such as schizophrenia and depression.
CONCLUSION
The regional distribution of 11C-MNPA was in good agreement with previous PET studies of dopamine D2 receptors in the human brain using antagonist radioligands such as 11C-raclopride and 11C-FLB457. The BPND values measured by the indirect kinetic model were in good agreement with those measured by the SRTM method with data acquired for 60 and 90 min. The BPND values measured by the transient equilibrium method also corresponded well with those measured by the indirect kinetic model with data acquired for 90 min. Simulation studies showed that errors in BPND measured by the SRTM method were small. The SRTM method with data acquired for 60 and 90 min is suitable for estimation of dopamine D2 receptor bindings using 11C-MNPA.
Acknowledgments
We thank Katsuyuki Tanimoto, Takahiro Shiraishi, and Akira Ando for their assistance in performing the PET experiments at the National Institute of Radiological Sciences. We also thank Yoshiko Fukushima of the National Institute of Radiological Sciences for her help as clinical research coordinator. This study was supported by a consignment expense for the Molecular Imaging Program on “Research Base for PET Diagnosis” from the Ministry of Education, Culture, Sports, Science and Technology, Japanese government.
Footnotes
-
COPYRIGHT © 2009 by the Society of Nuclear Medicine, Inc.
References
- Received for publication September 25, 2008.
- Accepted for publication January 21, 2009.