Abstract
Metabotropic glutamate receptor subtype 1 (mGluR1) is a crucial molecular target in the central nervous system disorders. 4-18F-fluoro-N-[4-[6-(isopropylamino)pyrimidin-4-yl]-1,3-thiazol-2-yl]-N-methylbenzamide (18F-FITM) has been recently developed as a useful PET ligand for mGluR1 imaging in our laboratory. In this study, we aimed to measure the affinity and density of mGluR1 using PET with 18F-FITM in rat brain under the in vivo conditions. Methods: Binding potentials (BPND) and amounts of specific binding (bound ligand concentration) at equilibrium state in brain regions were noninvasively estimated using the equilibrium analysis combined with the receptor-blocked approach (EA RBA) for kinetic analysis of 18F-FITM PET results in place of reference tissue methods. Using BPND and specific binding values of rats treated with multidose ligand, we performed Scatchard analyses for in vivo measurements of mGluR1 density (maximum number of binding sites, or Bmax) and ligand affinity (dissociation constant, or Kd) in brain regions, respectively. Results: The pretreatment of rats with unlabeled FITM (1 mg/kg) occupied an mGluR1 binding site of 18F-FITM by more than 99% and did not affect the input function. Hence, we used the tissue time–activity curve for receptor-blocked rats as representative of the nondisplaceable (free and nonspecific binding of radioligand) compartment. The BPND based on EA RBA showed a high correlation with the BPND based on invasive Logan plot graphical analysis in the thalamus, hippocampus, striatum, and cingulate cortex. The Kd (nM) and Bmax (pmol/mL) obtained by the Scatchard analyses with the multidose ligand assays were 2.1 and 36.3, respectively, for the thalamus; 2.1 and 27.5, respectively, for the hippocampus; 1.5 and 22.2, respectively, for the striatum; and 1.5 and 20.5, respectively, for the cingulate cortex with a high confidence. Conclusion: Our study is the first to our knowledge to measure the in vivo affinity (Kd and binding potential) of 18F-FITM and mGluR1 density (Bmax) with a high correlation to in vitro values in rat brain regions. This measurement using PET with 18F-FITM would be a useful index for research about mGluR1 functions in central nervous system disorders and development of new pharmaceuticals.
The excitatory amino acid neurotransmitter glutamate is known to act via ionotropic receptors (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid, N-methyl-d-aspartic acid, and kainite) and metabotropic receptors (metabotropic glutamate receptor subtypes 1–8 [mGluR1–8]) (1). MGluRs are members of the group C family of G-protein–coupled receptors and are divided into 3 different groups based on receptor sequence homology and physiologic activity. The mGluR1, along with mGluR5, is classified as group I (2). The activation of mGluR1 by endogenous glutamate stimulates intraneuronal cascades, which release Ca2+ from intracellular stores and activate protein kinase C (2). Thus, inhibition of mGluR1 activation has been indicated as potentially useful for neuroprotection against various injuries to the central nervous system, such as stroke (3–8), Parkinson disease (9–12), and Huntington disease (13).
Several radioligands for PET have been synthesized to investigate the in vivo localization and biologic characteristics of mGluR1 (14–17). Recently, we developed 4-18F-fluoro-N-[4-[6-(isopropylamino)pyrimidin-4-yl]-1,3-thiazol-2-yl]-N-methylbenzamide (18F-FITM) (Supplemental Fig. 1; supplemental materials are available online only at http://jnm.snmjournals.org) (18) as a useful PET ligand for mGluR1 imaging in the brain.
PET is an advanced molecular imaging modality that enables the study of the living human brain and, in particular, specific proteins involved in pathophysiology or as targets for therapeutic interventions. In addition to its diagnostic applications, PET has been used to assess new therapeutic agents for targeting neuropsychiatric and neurodegenerative disorders in the brain (19). Tracking the dynamic biodistribution of radiolabeled compounds on maps of living animal brains has been made possible by the development of high-resolution PET scanners (20,21). Thus, in vivo pharmacologic indices, such as the occupancy of the target binding sites by tested compounds and receptor density, can be commonly determined by PET with radioligands (22,23). Moreover, multiple PET scans using living rodents allow high-accuracy data to be easily acquired, providing obvious advantages over ex vivo assays. Although mGluR1 is a crucial target molecule for neuropharmacology, in vivo measurement studies of mGluR1 in the brain using PET with a specific radioligand have never been performed.
In the present study, we measured the binding potential (BP), density (maximum number of binding sites, or Bmax), and dissociation constant (Kd) of mGluR1 in the rat brain using 18F-FITM PET studies.
MATERIALS AND METHODS
General
18F-FITM (Supplemental Fig. 1) was synthesized by 18F-fluorination of a nitro precursor with 18F-KF/Kryptofix 222 (Merck) as described previously (18). Unlabeled FITM was synthesized in our laboratory according to previously reported procedures (24). The compound was dissolved in saline containing 20% ethanol and 10% polyoxyethylene sorbitan monooleate (polysorbate 80) and used for in vivo studies.
Sprague–Dawley rats were purchased from Japan SLC and kept in a temperature-controlled environment with a 12-h light–dark cycle; the rats were fed a standard diet (MB-1; Funabashi Farm). Animal experiments were performed according to the recommendations of the Committee for the Care and Use of Laboratory Animals, National Institute of Radiological Sciences.
Small-Animal PET Experiments
PET Procedure
A Sprague–Dawley rat was secured in a custom-designed chamber and placed in a small-animal PET scanner (Inveon; Siemens Medical Solutions). Body temperature was maintained using a 40°C water circulation system (T/Pump TP401; Gaymar Industries). A 24-gauge intravenous catheter (Terumo Medical Products) was placed in the tail vein of the rat for a bolus injection.
A bolus of 18F-FITM (17–18 MBq, 30–40 pmol, 0.1 mL) was injected via catheter in the tail vein of the rat, and dynamic emission scans in 3-dimensional list mode were obtained for 90 min (10 s × 6 scans, 15 s × 4 scans, 1 min × 5 scans, 2 min × 4 scans, and 5 min × 15 scans).
PET images were reconstructed using ASIPro VM (Analysis Tools and System Setup/Diagnostics Tool; Siemens Medical Solutions). Volumes of interest referencing the MRI template (25) were drawn on the cerebellum, thalamus, hippocampus, striatum, cingulate cortex, and pons. The sizes (cm2) in each region of interest were 1.5 for the cerebellum, 1.2 for the thalamus, 0.3 for the hippocampus, 0.7 for the striatum, 0.2 for the cingulate cortex, and 0.5 for the pons. Each PET image was overlaid on the MRI template, and the time–activity curve for each region was characterized. The radioactivity was decay-corrected to the injection time and is expressed as the percentage of injected dose per milliliter (%ID/mL).
Arterial Blood Sampling
Before PET assessment, a Sprague–Dawley rat (n = 4; male; age, 8 wk; weight, 270–290 g) had a polyethylene catheter (FR2; Imamura) inserted into the left femoral artery for blood sampling. Arterial blood (0.5–1 mL) was manually sampled at 10, 20, 30, 40, and 50 s and 1, 2, 3, 4, 5, 10, 15, 30, 60, and 90 min after the PET scan started. The radioactivity in whole blood and plasma was counted using a 1480 Wizard Autogamma scintillation counter (Perkin-Elmer). Radioactivity was corrected for decay. Plasma treatment and metabolite analysis were performed according to a previous report (26). Similarly, the blocking study using rats (Sprague–Dawley; n = 4; male; age, 8 wk; weight, 275–285 g) pretreated with unlabeled FITM (1 mg/kg) was also performed as described above.
Estimation of BPND Using PET Kinetics
Invasive Method
PET kinetic analyses were performed by Logan plot graphical analysis (GA) (27). The distribution volume (VT) of brain regions in the baseline and blocking studies was obtained using a plasma time–activity curve and tissue time–activity curve. All kinetic analyses were performed using PMOD software (version 3.2; PMOD Technologies).
To estimate nondisplaceable distribution volume (VND), the Lassen plot GA was performed as described in the supplemental file.
The distribution volume of specific binding (VS) was obtained by subtracting the VND based on the Lassen plot GA from the VT based on the Logan plot GA. The BPND, which equals the distribution volume ratio minus 1, was determined according to the following equation:
Noninvasive Method
The BPND was noninvasively estimated to perform the Scatchard analysis with the multidose ligand assays. The equilibrium analysis (EA), in which saturated specific binding of the radioligand at equilibrium state is calculated from the specific binding time curve (Cb) obtained by subtracting the tissue time–activity curve of the reference region (CND) from the tissue time–activity curve in a region of interest (CT), has been previously introduced (28). In this analysis, the equilibrium state was defined as occurring at the time point when dCb/dt for specific binding was 0. The equilibrium time was defined by fitting to a 3-exponential equation. On the other hand, the receptor-blocked approach (RBA), in which the nondisplaceable compartment could be obtained by saturating a receptor with unlabeled ligand instead of using a reference tissue, has been previously reported (29). When a receptor binding site was significantly blocked by excess unlabeled ligand treatment, binding on-rate (k3) of the radioligand would be close to zero, indicating that its tissue time–activity curve was parallel to the nondisplaceable compartment. In this study, we combined EA with RBA because there was no reference region using PET with this radioligand. According to this approach, the Cb in each brain region were determined by subtracting the tissue time–activity curve of the blocked rats as CND from the tissue time–activity curve of the baseline rats as CT, obtaining an equilibrium time using a 3-exponential fitting. The BPND based on EA RBA was calculated according to the following equations:
To validate BPND based on the EA RBA, we generated a correlational scatterplot of BPND based on the EA RBA against BPND based on the Logan plot GA in each brain region.
In Vivo Measurement of Bmax and Kd by 18F-FITM PET with Multidose Ligand Treatment
The PET procedure was performed as described in the “Small-Animal PET Experiments” section. Sprague–Dawley rats (n = 3; male; age, 7–10 wk; weight, 220–310 g) were treated with different doses of unlabeled FITM (0, 1, 5, or 30 μg/kg or 1 mg/kg) just before a bolus injection of 18F-FITM (17–18 MBq, 30–40 pmol, 0.1 mL). Estimations of equilibrium state and BPND were acquired in accordance with the approach as described in the “Estimation of BPND Using PET Kinetics” section.
Scatchard analyses were performed to estimate the Bmax and Kd of mGluR1 with 18F-FITM. For the analysis, the amount of specific binding was plotted against BPND instead of the ratio of bound ligand concentration to free ligand concentration under the in vitro condition. Values were calculated according to the following equations:
Linear regression for the estimation of Bmax and Kd values for mGluR1 with 18F-FITM in the thalamus, hippocampus, striatum, and cingulate cortex was performed using GraphPad Prism 5 (GraphPad Software). In this plot, Kd is calculated as the multiplicative inverse of the slope of the regression line, and Bmax is determined from the x-axis intercept.
RESULTS
Estimation of BPND Using PET Kinetics
Figure 1 shows the input functions (plasma time–activity curve) of 18F-FITM in baseline and blocked rat brains.
Time–activity curves of 18F-FITM in plasma of baseline and blocked rats. Blocked rats were pretreated with unlabeled FITM (1 mg/kg). Plasma input function was corrected by injection time and metabolites of 18F-FITM. Mean ± SE (n = 4 in each group).
Plasma time–activity curves in both baseline and blocked rats showed initial uptake by 7%–8 %ID/mL just after a bolus injection and rapid clearance by roughly 0.1 %ID/mL at 15 min. The pretreatment of rats with unlabeled FITM (1 mg/kg) barely influenced the input function. In the baseline tissue time–activity curves, remarkable accumulation of radioactivity was detected in the cerebellum; moderate uptake was detected in the thalamus, striatum, hippocampus, and cingulate cortex; and little uptake was found in the pons (Supplemental Fig. 2A). The pretreatment of rats with unlabeled FITM decreased uptake significantly and made it difficult to distinguish each region of the brain (Supplemental Fig. 2B).
The VT based on Logan plot GA in baseline and blocked rats is shown in Table 1. In the baseline rats, the rank order of VS and BPND using the VND value based on the Lassen plot GA (Supplemental Fig. 3) among brain regions was cerebellum >> thalamus > striatum > hippocampus > cingulate cortex >> pons (Table 1), corresponding to the reported mGluR1 densities in rat brain (31).
Results of Logan Plot GA in Brain Regions of Baseline and Blocked Rats
Figure 2 shows the results of the EA RBA in the thalamus, hippocampus, striatum, and cingulate cortex.
Identification of time point at which equilibrium occurs for mGluR1 with 18F-FITM in thalamus (A), hippocampus (B), striatum (C), and cingulate cortex (D). Measured values in baseline experiments were subtracted from measured values in blocked experiments. Difference was defined as specific binding and fitted to set of 3 exponentials. Each equilibrium time is indicated by black arrow. Points are expressed as mean ± SE (n = 4).
The equilibrium time was 90 min for the thalamus, 76 min for the hippocampus, 74 min for the striatum, and 65 min for the cingulate cortex. BPND according to Equation 3 was 18.4 ± 1.9 for the thalamus, 13.6 ± 1.5 for the hippocampus, 14.5 ± 1.5 for the striatum, and 13.0 ± 1.4 for the cingulate cortex.
Figure 3 shows the results of linear regression using plots of BPND based on the EA RBA against BPND based on the Logan plot GA. The regression line had good confidence (R2 = 0.93), and the slope was 0.90. Although the BPND based on EA RBA showed a result of roughly 10% underestimations, compared with that based on Logan plot GA, the good linear regression result supported the use of the EA RBA instead of reference tissue methods. The EA RBA could not be used in the cerebellum and pons because of the poor regression in those regions.
Correlational scatterplots between BPND values based on EA RBA and those based on Logan plot GA. Plots show average datum points, with SEs in x- and y-axes (n = 4). ○ = thalamus; □ = hippocampus; △ = striatum; ⋄ = cingulate cortex.
In Vivo Measurement of Bmax and Kd by 18F-FITM PET with Multidose Ligand Treatment
Figure 4 shows representative PET images of the thalamus/hippocampus, striatum, and cingulate cortex on the coronal slices in baseline rats and rats pretreated with different doses of unlabeled FITM (1, 5, and 30 μg/kg and 1 mg/kg, Figs. 4B, 4C, 4D, and 4E, respectively). With increasing doses of unlabeled FITM, the radioactivity showed a stepwise attenuation in all regions.
PET of 18F-FITM distribution in rat brains after intravenous injection of unlabeled FITM. Rats were pretreated with different doses of unlabeled FITM: baseline (A), 1 μg/kg (B), 5 μg/kg (C), 30 μg/kg (D), and 1 mg/kg (E) administered). Coronal PET images were generated by averaging dynamic data at 0–90 min after bolus injection of 18F-FITM. Volumes of interest were placed on hippocampus, thalamus, striatum, and cingulate cortex. Hi = hippocampus; Th = thalamus; St = striatum; Ci = cingulate cortex.
The dose-dependent reduction of 18F-FITM binding to mGluR1 in the multidose ligand assays was shown by alterations to the time curve of specific binding in the thalamus, hippocampus, striatum, and cingulate cortex (Supplemental Fig. 4). The transient equilibrium of specific binding (dCb/dt = 0) based on EA RBA was reached at 70–90 min in the thalamus, 55–75 min in the hippocampus, 60–80 min in the striatum, and 50–70 min in the cingulate cortex (Supplemental Fig. 4).
Figure 5 shows the Scatchard plots for the thalamus, hippocampus, striatum, and cingulate cortex. Table 2 shows Bmax, Kd, the coefficient of determination (R2), and BP (Bmax/Kd) values in each brain region. The Kd values in all brain regions were in the range of 1.5–2.1, with good confidence of linear fit (R2 = 0.87–0.96).
Scatchard analyses for determination of mGluR1 Bmax and Kd of 18F-FITM in thalamus (A), hippocampus (B), striatum (C), and cingulate cortex (D). Regression lines are indicated for each region. The plots were expressed by average datum points, with standard errors in x- and y-axis (n = 3). [B] = bound ligand concentration.
In Vivo Measurements of Kd, Bmax, and BP by Scatchard Analyses in Rat Brain
DISCUSSION
To our knowledge, the present 18F-FITM PET studies provided in vivo measurements of Bmax and Kd in rat brain for the first time.
We have developed 18F-FITM as a novel radioligand of mGluR1 imaging and demonstrated its potential in in vitro and in vivo studies of mGluR1 in rodent and nonhuman primate brains (18,32). In those studies, we confirmed the pharmacologic selectivity of 18F-FITM on the basis of competition with a mGluR1-selective ligand such as JNJ16259685. In a preliminary experiment, we examined brain uptake of 18F-FITM in wild-type and mGluR1 knockout mice. Brain uptake of radioactivity was barely detected in the PET image of the mGluR1 knockout mouse (Supplemental Fig. 5). These results directly demonstrate the high specificity of 18F-FITM for mGluR1, indicating that this radioligand is a promising tracer suitable for the accurate in vivo measurement of pharmacologic indices for mGluR1 in the brain.
BP, defined as the ratio of bound ligand concentration to free ligand concentration, was introduced in the earliest days of PET application to the study of neuroreceptors in the brain (33). Scatchard analysis has often been performed for the in vivo determination of Bmax and Kd using BP values and the amount of specific binding of the ligand (34–36). In general, the BP, under typical assumptions, can be directly calculated only from brain data using a variety of reference tissue methods such as BPND in vivo.
The brain stem is known to have negligible expression of mGluR1 by immunohistochemistry (31) and autoradiograms with several radioligands including 18F-FITM (14,15,17,18,37). Therefore, we assumed that the pons, a part of the brain stem, could be used as a reference region. To verify the aptitude of the pons, we measured the VT of rat brain regions by Logan plot GA in baseline and blocking studies. Unfortunately, apparent specific binding in the pons was detected; the VT values in baseline and blocking studies were 6.2 and 2.7, respectively (Table 1). Thus, the reference region methods were not adaptable to in vivo measurement for pharmacologic indices of this radioligand because of an underestimation of BPND. In fact, the BPND based on the simplified reference tissue model (38) using the tissue time–activity curve of the pons was underestimated by roughly 40%–50% in all brain regions, compared with the BPND (distribution volume ratio − 1) based on a 2-tissue-compartment model.
To noninvasively measure BPND, we used the RBA, which was able to estimate the distribution volume of the free and nonspecific binding compartment, instead of using the reference tissue (29). However, with this approach, there is a potential concern that the input function may be changed by treatment with excess unlabeled ligand and that the injected unlabeled ligand dose may not fully occupy a binding site of the receptor. In the present study, the plasma time–activity curve of 18F-FITM with blockade of specific binding did not show an acute increase because of overflow from the peripherals at the initial phase (Fig. 1). Although the initial peak in the plasma time–activity curve in the blocked rats occurred earlier than that in the baseline rats, this potential concern is unnecessary because all tissue time–activity curves peaked more than 45 min after the 18F-FITM injection (Supplemental Fig. 2A). Another concern is the ligand occupancy of the targeted receptor. As shown in the tissue time–activity curves of the blocked rats (Supplemental Fig. 2B), the uptake was significantly lower than that in baseline tissue time–activity curves in each brain region measured. The ligand occupancy and VND were calculated using Lassen plot GA with linear regression (Supplemental Fig. 3). It was shown that pretreatment with 1 mg of unlabeled FITM per kilogram occupied mGluR1 binding sites at more than 99% in each brain region. This regression analysis showed that the tissue time–activity curve for each brain region in the blocked rats was almost equal to the nondisplaceable compartment, which can be used instead of reference tissue.
To perform the Scatchard analysis with the multidose ligand assays, we estimated BPND based on the EA RBA in rat brain regions. Specific binding time curves were obtained by EA (Fig. 2) (28). For validation of the EA RBA, we generated a correlational scatterplot between BPND based on the EA RBA and that based on the Logan plot GA (Fig. 3). The regression line showed good correlation in the thalamus, hippocampus, striatum, and cingulate cortex. However, plots of the cerebellum, a region with abundant mGluR1, and the pons, a region with negligible mGluR1, were removed from the regression line because of their nonadaptability to the EA RBA. In the cerebellum, the equilibrium state of 18F-FITM has been previously shown to be over 180 min (32). The in vivo measurements of Bmax and Kd of mGluR1 with 18F-FITM in the cerebellum were not suited to the EA RBA, because the EA could be applied in vivo only when a specific binding reached a peak within the time span of a PET experiment (28). More importantly, PET with this radioligand was suited to quantification of mGluR1 in moderate- or low-density regions, although the in vivo quantification of mGluR1 with 18F-FITM in the cerebellum would be inferior to that with 18F-MK-1312 (16), a promising PET ligand for cerebellar mGluR1. 18F-MK-1312 PET analysis with a 2-tissue-compartment model in a monkey showed good kinetic parameters with high confidence in the cerebellum, but the directly estimated BPND (k3/k4) was under 1.0. In addition, significant specific binding of 18F-MK-1312 was detected in the cerebellum but not in the thalamus and striatum.
On the other hand, the pons, a region with negligible mGluR1, is located near the cerebellum, the region with the highest accumulation of radioactivity. Therefore, nonsaturated radioactive signals in the cerebellum would spill into the pons. This profile might be the cause of the nonadaptability of reference tissue models to the pons. However, the rat brain is considerably smaller than the human brain, and partial-volume errors cause spill-in of radioactivity from adjacent areas. The larger human brain may have regions with no (or limited) displaceable uptake that could hypothetically be used as a reference region. Although 18F-FITM with EA RBA was not useful in the cerebellum and pons, which in the moderate- or low-mGluR1 brain regions—such as the thalamus, hippocampus, striatum, and cingulate cortex—demonstrated the ability to noninvasively quantify mGluR1.
Interestingly, the expression levels of mGluR1 in animal models of central nervous system disorders, such as Parkinson disease and Alzheimer disease, are altered in the thalamus, striatum, and cerebral cortex but not in the cerebellum (9,39,40). To measure Bmax and Kd values in vivo in the thalamus, hippocampus, striatum, and cingulate cortex, we performed Scatchard analyses with multidose ligand assays using the BPND and the amount (pmol/mL) of specific binding based on the EA RBA (Fig. 5). As shown in Table 2, Kd and Bmax values in each region were calculated from a regression line with high correlation (R2 = 0.87–0.96).
We previously measured the in vitro Kd value (1.9 nM) of 18F-FITM using a rat brain homogenate (37). Theoretically, the Kd value is a fundamentally impartial value in both in vitro and in vivo conditions. However, in vivo measurements are sometimes different from in vitro measurements because of the metabolism of the radioligand, high amounts of nonspecific binding, the influence of transporters, and a variety of other reasons. In the present study, the in vivo Kd values in the thalamus, hippocampus, striatum, and cingulate cortex were all in the range of 1.5–2.1, approximating the in vitro measured value. These results suggested that PET with 18F-FITM is not affected by in vivo factors. The EA RBA with 18F-FITM demonstrated relatively accurate measurement of mGluR1 density (pmol/mL) in the diencephalon and telencephalon, such as 36.3 for the thalamus, 27.5 for the hippocampus, 22.2 for the striatum, and 20.5 for the cingulate cortex.
CONCLUSION
We have determined the Bmax and Kd of mGluR1 in the thalamus, hippocampus, striatum, and cingulate cortex using 18F-FITM PET combined with the EA RBA using multidose ligand assays. 18F-FITM PET is a promising tool for the further understanding of the mGluR1 function under pathophysiologic conditions. The present results will aid the development of pharmaceuticals targeting central nervous system disorders, for which it is important to measure indices such as mGluR1 occupancy and to determine effective dose in a target brain region.
DISCLOSURE STATEMENT
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
We thank Tatsumi Hirata (National Institute of Genetics, Shizuoka, Japan) for development of the mGluR1 knockout mouse. We are also grateful to the staff of the National Institute of Radiological Sciences for assistance with the cyclotron operation, radiosynthesis, and animal experiments. This work was partly supported by a Grant-in-Aid for Basic Research (C) (22591379) from the Ministry of Education, Culture, Sports, Science and Technology of the Japanese Government. No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Aug. 22, 2012.
- © 2012 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication March 13, 2012.
- Accepted for publication May 22, 2012.