Abstract
N-methyl d-aspartate (NMDA) ion channels play a key role in a wide range of physiologic (e.g., memory and learning tasks) and pathologic processes (e.g., excitotoxicity). To date, suitable PET markers of NMDA ion channel activity have not been available. 18F-GE-179 is a novel radioligand that selectively binds to the open/active state of the NMDA receptor ion channel, displacing the binding of 3H-tenocyclidine from the intrachannel binding site with an affinity of 2.4 nM. No significant binding was observed with 10 nM GE-179 at 60 other neuroreceptors, channels, or transporters. We describe the kinetic behavior of the radioligand in vivo in humans. Methods: Nine healthy participants (6 men, 3 women; median age, 37 y) each underwent a 90-min PET scan after an intravenous injection of 18F-GE-179. Continuous arterial blood sampling over the first 15 min was followed by discrete blood sampling over the duration of the scan. Brain radioactivity (KBq/mL) was measured in summation images created from the attenuation- and motion-corrected dynamic images. Metabolite-corrected parent plasma input functions were generated. We assessed the abilities of 1-, 2-, and 3-compartment models to kinetically describe cerebral time–activity curves using 6 bilateral regions of interest. Parametric volume-of-distribution (VT) images were generated by voxelwise rank-shaping regularization of exponential spectral analysis (RS-ESA). Results: A 2-brain-compartment, 4-rate-constant model best described the radioligand’s kinetics in normal gray matter of subjects at rest. At 30 min after injection, 37% of plasma radioactivity represented unmetabolized 18F-GE-179. The highest mean levels of gray matter radioactivity were seen in the putamina and peaked at 7.5 min. A significant positive correlation was observed between K1 and VT (Spearman ρ = 0.398; P = 0.003). Between-subject coefficients of variation of VT ranged between 12% and 16%. Voxelwise RS-ESA yielded similar VTs and coefficients of variation. Conclusion: 18F-GE-179 exhibits high and rapid brain extraction, with a relatively homogeneous distribution in gray matter and acceptable between-subject variability. Despite its rapid peripheral metabolism, quantification of 18F-GE-179 VT is feasible both within regions of interest and at the voxel level. The specificity of 18F-GE-179 binding, however, requires further characterization with in vivo studies using activation and disease models.
N-methyl-d-aspartate (NMDA) receptors for l-glutamate, the major excitatory neurotransmitter in the central nervous system (CNS), are linked to ligand- and voltage-gated ion channels (1). Excitatory synaptic transmission via these receptors mediates neuroplasticity and is necessary for learning and memory. In pathologic circumstances, however, NMDA receptor overactivation may be associated with cell death (2).
Excessive NMDA receptor activation mediates excitotoxic neuronal injury after acute cerebral insults (3) and is thought to contribute to disorders of neuronal hyperexcitability (e.g., epilepsy) and chronic neurodegenerative (e.g., Alzheimer, Huntington) (4) and psychotic (5) disorders. Hence, there is interest in the development of radioligands to allow assessment of NMDA receptor function in humans in vivo.
Imaging NMDA ion channel function in vivo in humans has been challenging because of low brain uptake, low affinity for the NMDA receptor, high rates of dissociation, rapid metabolism, and high nonspecific binding of candidate tracers. To date, only 6 of more than 60 radioligands designed for in vivo imaging of the NMDA receptor system (6) have progressed to PET or SPECT studies in humans (Table 1). Of these, 5 act at the phencyclidine recognition site (7), which lies within the NMDA ion channel pore, and hence these radioligands require receptor activation to allow binding. The diarylguanidine CNS 5161, a noncompetitive antagonist, is one such use-dependent ligand that demonstrated high uptake in the putamen and thalamus in human studies in vivo (8–10). Quantification of 11C-CNS 5161 proved to be difficult, its radiochemical yield low, and its metabolism fast; thus, it was not pursued as a PET radioligand.
In vitro evaluation of a 18F-labeled analog of CNS 5161 (18F-GE-179 (11)) showed that 18F-GE-179 displaces 3H-TCP (-3H-N-(1-[thienyl] cyclohexyl)piperidine) from the phencyclidine site with a Ki (inhibition constant) of 2.4 nM (similar to that seen with 3H-MK-801 (3H-(+)-5-methyl-10,11-dihydro-5H-dibenzo[a,d]cyclo-hepten-5,10-imine maleate) and CNS 5161 (12)) and GE-179 inhibits glutamate-mediated influx of Ca2+ in teratocarcinoma 2 (NT2) cells differentiated to express human NMDA receptors in a dose-dependent manner (supplemental data; supplemental materials are available at http://jnm.snmjournals.org). GE-179 has a distribution coefficient (log D7.4) of 2.49 ± 0.1 (11), similar to that of CNS 5161 (log P = 1.92 ± 0.26 (13)) and indicative of moderate lipophilicity. In a binding study to characterize GE-179 binding to CNS targets other than the NMDA phencyclidine site, 10 nM GE-179 did not significantly bind to any of the 60 receptors, channels, and transporters assessed. Binding of 3H-ifenprodil to the rat σ 2 receptor was only inhibited by 8%; all other targets were either inhibited by 3% or less at 10 nM or had not shown significant binding (±20% inhibition) at the higher dose of 1 μM (supplemental data).
Here, to our knowledge, we report the first-in-human use of 18F-GE-179 PET in healthy human subjects. Our primary objective was to evaluate the kinetic behavior of the radioligand in vivo in normal human brain.
MATERIALS AND METHODS
This study was approved by the Research Ethics Committees of the Royal Marsden Hospital, Imperial College Healthcare NHS Trust, and University College London Hospitals NHS Foundation Trust. Permission to administer 18F-GE-179 was obtained from the Administration of Radioactive Substances Advisory Committee, U.K. All participants provided written, informed consent before participation.
Nine healthy participants without history of neurologic or psychiatric illness were scanned. The demographic and descriptive data for the subjects are presented in Table 1.
18F-GE-179 was synthesized by Hammersmith Imanet Limited as previously described (11). Details of the injectate are presented in Table 1.
PET Image Acquisition and Preprocessing
Images were acquired using an ECAT EXACT3D 962 HR+ PET camera (Siemens/CTI) with a 15.5-cm axial and 58.0-cm transaxial field of view. After a transmission scan, each subject underwent a 90-min dynamic emission scan following an intravenous bolus injection of median 186 MBq of 18F-GE-179. Thirty-four frames of increasing duration were acquired. Data were reconstructed using a Fourier rebinning algorithm (FORE; (14)) and 2-dimensioanl filtered backprojection (ramp, kernel 2.0 mm in full width at half maximum). The voxel size of the reconstructed images was 2.092 × 2.092 × 2.42 mm.
The images were acquired in a quiet, dimly lit room without background noise. Head position was monitored throughout with the camera’s laser light. Continuous arterial blood sampling at a rate of 5 mL per minute was performed from 0 to 15 min. Discrete arterial blood samples were taken at baseline and a further 8 times during the next 90.5 min. High-performance liquid chromatography was used to quantify radiolabeled metabolites in the plasma.
The fraction of plasma radioactivity attributable to the parent 18F-GE-179 was fitted to a sigmoid function normalized to unity at 0 min using CLICKFIT version1.7 (Hinz R, Cunningham VJ, Imaging Research Solutions Limited, London, U.K.) running in MATLAB 6.5 (The MathWorks Inc.) to generate a metabolite model for each participant.
The area under the metabolite model curve (AUCmetabs) was used as a measure of the rate of metabolism for each individual. Correlations between the AUCmetabs and age, weight, and body mass index (BMI) were examined by Spearman rank (ρ) correlation coefficient using SPSS (version 16.0 [SPSS Inc., IBM Corp.] for Windows [Microsoft]). The threshold of statistical significance (P < 0.05) was corrected for multiple comparisons using the Bonferroni method.
Continuous decay- and metabolite-corrected parent plasma input functions were generated for all subjects using CLICKFIT, as described previously (15). Briefly, the final input functions were generated by multiplication of the calibrated time course of radioactivity in the blood with the fits to the plasma-to-blood ratios and metabolite fractions. Rodent in vivo metabolism data showed that the likelihood of ingress of radiolabeled metabolites in to the CNS was low; this, therefore, was not accounted for (see supplemental data). The areas under the decay- and metabolite-corrected parent plasma input function curves (AUCIF) were calculated for each subject by summation of the counts in the parent plasma input function file, after correction for 18F radiodecay.
Correlations between the AUCIF and age, weight, BMI, and injected dose were examined by Spearman rank (ρ) correlation coefficient. We hypothesized a priori that the AUCIF would be negatively correlated with weight or BMI and positively correlated with injected dose; these correlations were not, therefore, corrected for multiple comparisons.
Attenuation-corrected dynamic PET images were corrected for head motion using a post hoc frame-to-frame realignment method running in MATLAB 7.4, with the corresponding nonattenuated image used as a reference (Piwave8.0, internal software (16)). Summation images (KBq/mL) that were weighted relative to frame duration were created for 5–30, 30–60, and 60–90 min from the attenuation- and motion-corrected dynamic image file using Receptor Parametric Mapping (version 6 (17)) software running in MATLAB 5.0.
Global and Regional Distribution of Radioactivity
Global (whole-brain) radioactivity concentration (i.e., uptake [KBq/mL]) in the summation images was computed using a whole-brain mask, which was created for each of the 9 participants from their 5- to 30-min summed images using the extract region tool implemented in Analyze8.1 (AnalyzeDirect, Inc.); the mean over the entire matrix within this mask was then computed.
Radioactivity concentration in the gray matter of the summation images was computed as follows. First, an individualized gray matter mask was created for each of the 9 participants by segmentation of their T1-weighted MR image according to the Unified Segmentation method (18) implemented in SPM8 (Statistical Parametric Mapping 8; Wellcome Trust Centre for Neuroimaging, UCL, London) and by thresholding of the resulting gray matter probability map at 0.5 using Analyze 8.1. The 9 individual participants’ MR images and gray matter masks were coregistered to their 5- to 30-min summation images, before calculation of the mean radioactivity concentration. White matter radioactivity concentrations were computed with white matter masks that were created simultaneously.
An 83-region gray matter–only region-of-interest (ROI) map was produced for each of the 9 participants by transformation of 30 manually created atlases (19) to the participant’s native MR imaging space using multiatlas propagation with enhanced registration and decision fusion (20). The resulting atlas was then multiplied by the thresholded gray matter component of the MR image and used to sample the radioactivity concentration (KBq/mL) in each of the 83 regions.
The median radioactivity concentrations were calculated for the participants, surveying all 83 regions and combining left and right homologs where possible.
The median asymmetry indices (AIs) were calculated to compare radioactivity concentration (RC) in the ROIs highlighted by the survey, according to the following formula:
Compartmental Modeling of 18F-GE-179 Cerebral Tissue Kinetics
Six gray matter–only ROIs were defined a priori for modeling purposes according to the known distribution of the NMDA receptor (21) and sampled with left and right sides combined. In order of decreasing approximate receptor density, the ROIs were the putamina, thalami, superior frontal gyri, parahippocampal gyri, occipital lobes, and cerebellum. The entire brain gray matter and white matter ROIs were also sampled. Each participant’s attenuation- and motion-corrected dynamic PET image was sampled over all 34 time frames. Time–activity curves for each ROI were produced using CLICKFIT from the scanner information file and a weights file that had been generated for each ROI from the sampled dynamic PET data.
Analyses of 18F-GE-179 cerebral tissue kinetics were performed using the following standard compartmental models: 1-brain-compartment, 2-rate constant (1c2kbv, reversible binding); 2-brain-compartment, 3-rate-constant (2c3kbv, irreversible binding); 2-brain-compartment, 4-rate-constant (2c4kbv, reversible binding); and 3-brain-compartment, 6-rate-constant (3c6kbv, reversible binding). All models incorporated a variable blood volume component and were used to calculate the delay between plasma and tissue time–activity curves.
The model fit was assessed using the Akaike Information Criterion (AICw) and the intraregion and between-subject coefficient of variation (CV) of volume-of-distribution (VT) estimates.
Strengths of regional VT and K1 correlations were assessed by Spearman rank (ρ) correlation coefficients, pooling all ROI data for the 9 participants.
The correlations between gray matter VT and age, weight, BMI, injected dose, specific radioactivity at the time of injection, and AUCIF were similarly interrogated with the Spearman rank (ρ) statistic, with correction for multiple comparisons using Bonferroni method.
Generation of Parametric VT Images
Parametric VT images were generated by voxelwise rank-shaping regularization of exponential spectral analysis (RS-ESA) (22). RS-ESA is a development of the model-free exponential spectral analysis estimation method (23) that has been optimized for noisy datasets by incorporating a singular value decomposition of the exponential base. Time constants were specified as 5 and 5,100 s. The noise fraction was specified as 0.15.
The median VT was computed within the ROIs and compared with the VTs derived via compartmental modeling. The correlations between regional VTs derived using each method were assessed by Spearman rank (ρ) correlation coefficient, pooling all ROI data for the 9 participants.
RESULTS
Metabolism
Unmetabolized 18F-GE-179 accounted for a mean of 50% of the plasma radioactivity at approximately 16 min and 25% at 44 min (Fig. 1).
There was a significant negative correlation between the AUCmetabs and BMI (r = −0.86, P = 0.002). There were no significant correlations with other variables.
Decay- and Metabolite-Corrected Parent Plasma Input Functions
The parent plasma input functions peaked at a median of 71.5 s (interquartile range [IQR], 65.6–81.7 s) and decreased to less than 10% of the peak activity concentration within 5 min (Fig. 2). The radioactivity concentrations in plasma at the peak of the parent plasma input functions had a median of 47.3 kBq/mL (IQR, 43.4–65.1 kBq/mL), equivalent to 0.02% of the injected dose/mL.
The median AUCIF was 4782.0 kBq⋅cm−3 (IQR, 4,165.4–5,754.8 kBq/mL). There was a significant negative correlation between the AUCIF and weight and AUCIF and BMI (both r = −0.721, P = 0.019). There were no significant correlations with other variables.
Regional Distribution of Radioactivity
The median gray matter radioactivity concentration was 7.3 kBq/mL (IQR, 5.9–9.2 kBq/mL) in the 5- to 30-min summation images. The highest median radioactivity concentration was identified in the putamina, followed by (in descending order) the cunei, lingual gyri, thalami, and cerebellum (Table 2). The between-subject CVs for the gray matter radioactivity concentration ranged from 25.4% to 29.2%.
The lowest median radioactivity concentration in gray matter was observed in the anterior medial temporal lobes, followed by (in ascending order) the parahippocampal gyri, hippocampi, and superior frontal gyri (Table 2).
Left–Right Asymmetry of 18F-GE-179 Radioactivity Concentration
The range of median asymmetry indices in the regions reported above was from –0.03 (IQR, –0.05–0.01) in the thalami to 0.03 (IQR, 0.00–0.05) in the lingual gyri (Table 2).
Time–Activity Curves
The mean gray matter radioactivity concentration peaked at 8.9 kBq/mL, equivalent to a mean of 0.0048 percentage injected dose (%ID) at 7.5 min. The mean white matter radioactivity peaked at 5.1 kBq/mL (0.0028 %ID) at 13 min. The radioactivity declined more quickly in gray than in the white matter.
The highest mean gray matter radioactivity concentration was seen in the putamina, which peaked at 10.4 kBq/mL (0.0056 %ID) at 7.5 min (Fig. 3), followed by the thalami, occipital lobes, cerebellum, superior frontal gyri, and parahippocampal gyri. The radioactivity concentration in the putamina decreased to approximately 58% of the peak at 91.5 min, to 63% in the thalami, 53% in the occipital lobes, 51% in the cerebellum, and 64% in the superior frontal gyri. The radioactivity concentration peaked at 13 min in the parahippocampal gyri, later than in other ROIs, but also declined more slowly than in the other gray matter regions (68% of the peak at 91.5 min).
Compartmental Modeling of 18F-GE-179 Cerebral Tissue Kinetics
The 2c4kbv model best described the radioligand’s kinetics in gray matter (AICw = −207; Table 3; Fig. 4), with a slight bias overestimation of peak radioactivity concentration. In contrast, the 1c2kbv and 2c3kbv models were characterized by consistently higher AICw, for example, −114, in gray matter, and the 3c6kbv model was characterized by a large number of outlier values within each ROI.
The median estimates of VT in each ROI, derived by application of the 2c4kbv models to individual time–activity curves, are shown in Table 3. The between-subject CVs for VT in gray matter regions ranged from 12.9% to 16.0%.
A positive correlation was detected between VT and K1 (Spearman ρ = 0.398; P = 0.003; Fig. 5).
There was a positive correlation between gray matter VT and age (Spearman ρ = 0.803; P = 0.009). There were no significant correlations with other variables.
Voxelwise Rank-Shaping (RS-ESA)
A representative example of an 18F-GE-179 VT image is provided in Figure 6. RS-ESA yielded VTs (Table 4) that correlated with those derived from regional compartmental modeling (Spearman ρ = 0.901; P < 0.001; Fig. 5). The between-subject CVs for VT in gray matter regions ranged from 11.6% to 14.8%. A bias toward underestimation of VT, relative to that derived from the 2c4kbv model, was observed (Fig. 7).
DISCUSSION
To our knowledge, we describe the first use of 18F-GE-179 as a PET radioligand in humans.
Brain penetration of 18F-GE-179 was high, with a mean peak radioactivity concentration of 8.9 kBq/mL in gray matter, which is greater than that achieved when using the NMDA radioligands 18F-memantine (24), 11C-AcL703 (25), and 11C-CNS 5161 (8), normalized to injected volume. The mean peak global and gray matter radioactivity concentrations have not been reported with other radioligands targeted at the NMDA receptor (26–30).
The distribution of radioactivity in gray matter after injection of 18F-GE-179 was relatively homogeneous, with the highest concentration (in the absence of correction for partial-volume effects) observed in the putamina, a region of moderate NMDA receptor density (31). Peak and summed (5–30 min) radioactivity concentrations in the white matter, often considered a reference region for nonspecific binding (32), were just over half that of the putamina. High radioactivity concentration was observed in the cerebellar cortex, in keeping with the documented existence of NMDA receptors in this region (33).
The channel maximal open probability of the NR2A subunit is approximately 2–5 times higher than that of the NR2B subunit (34), 50 times higher than that of the NR2C (35), and 9 times higher than that of NR2D subunits (36). NMDA receptors are particularly concentrated in the hippocampi (31), with NR2B subunits present in higher concentrations than NR2A (37); the lower channel opening probabilities of this subunit might underlie the modest radioactivity concentration observed in medial temporal regions, which was nearly equivalent to that seen in white matter. We also observed low radioactivity concentration in the brain stem, which is consistent with the lack of NR2A and NR2B expression in this region (37), suggesting specificity of gray matter binding. However, our finding of a weak-to-moderate positive correlation between K1 and VT suggests partial cerebral blood flow-dependency.
Ideally, the specificity of 18F-GE-179 binding in vivo would be quantified by coadministration of unlabeled GE-179 or another unlabeled agent with high affinity and selectivity for the phencyclidine-binding site. However, this is not feasible in humans because of poor tolerability of blocking doses of NMDA antagonists. A dose-escalation safety study with CNS 5161 was abandoned after a sustained systolic blood pressure increase in 1 of 2 participants administered 750 μg (38). The performance of a cognitive task, such as an episodic memory task, during radioligand injection might be expected to increase specific binding in certain brain areas, that is, in the medial temporal lobe. However, such a paradigm would have to be performed for approximately 20 min or more because of the rate of uptake of the radioligand, a second baseline scan at rest would be required, and the increased local cerebral blood flow would confound the analyses. Focal epileptic activity could also be hypothesized to increase (specific) binding in the medial temporal lobe regions. 18F-GE-179 will be investigated in participants with epilepsy in a subsequent study.
Additional limitations of this first-in-human study include the small sample size and lack of test–retest data. A further study is required to demonstrate acceptable within-subject (i.e., test– retest) variability for the quantification of 18F-GE-179 VT. However, in view of the low NMDA–channel opening probability at baseline, test–retest data at rest without control over cognitive and emotional state might not be helpful.
The rate of metabolism of 18F-GE-179 was less rapid than that of 11C-CNS 5161 (8) and 123I-CNS 1261 (29) but not 18F-memantine (24) (data not reported for other NMDA-receptor radioligands). Our analyses assumed radiolabeled metabolites did not penetrate the brain, although 1 more polar compound was identified in rodent brains at 30 min after injection (GE Healthcare, unpublished data, on file, 2009). Parent 18F-GE-179 accounted for 84% of the cerebral radioactivity at this point. In humans, the slower plasma metabolism of the radioligand suggests it should account for an even greater proportion at this time point.
18F-GE-179 exhibited faster gray matter tissue uptake and washout than 11C-CNS 5161 (8) (and also 123I-MK-801 (28,30), 18F-memantine (24), and 11C-AcL703 (25)). Model fits indicated a better fit for the reversible model. Our data indicate that quantification of regional 18F-GE-179 VT is achievable with a 2c4kbv with variable blood component. Voxelwise quantification of 18F-GE-179 VT by RS-ESA yielded comparable estimates and parametric maps of acceptable visual quality; this will also allow for comparison of VT between populations without a priori delineation of ROI whereas voxelwise compartmental modeling is hampered by noise. When ROI quantification is preferred, we suggest that the 2c4kbv model is used to avoid the negative bias associated with RS-ESA. However, the choice of quantification method will also be informed by test–retest variability, when such data become available.
18F-GE-179 uptake and VT have acceptable between-subject variability in healthy participants, as evidenced by regional between-subject CVs of less than 16%. A positive correlation between VT and age was observed and so future studies, especially if involving poorly matched subgroups, should model this as a covariate in the analyses. This finding appears at odds with the age-related declines in cerebral blood flow (39) and NMDA receptor function (40). However, it is possible that the relationship between 18F-GE-179 VT and age is nonlinear. In such a case, the positive correlation we have identified might reflect the limited age spread (range, 25–62 y) of our cohort. As expected, participants of higher weight or BMI had a lower parent plasma input than the remaining controls. Moreover, regional VT was seen to correlate with the first rate constant (i.e., linked to cerebral blood flood).
CONCLUSION
Our first-in-human evaluation of 18F-GE-179 has demonstrated several properties that are desirable in PET radioligands, such as high brain uptake and acceptable between-subject variability. Quantification of VT appears feasible within ROIs and at the voxel level. The specificity of 18F-GE-179 binding requires further characterization and could potentially limit its application. Further in vivo evaluation of 18F-GE-179 is warranted.
DISCLOSURE
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. This study was supported by the Medical Research Council (MRC) Clinical Sciences Centre; GE Healthcare plc; the Department of Health NIHR Biomedical Research Centers’ funding scheme; a MRC Doctoral Training Account (3 + 1) studentship that was awarded by Imperial College London; Epilepsy Society, UCL, UCL Hospitals, and UCLH/UCL Biomedical Research Centre; an MRC Clinician Scientist Fellowship (G108/585); and the Neurodis Foundation. Sajinder K. Luthra, Paul A Jones, William Trigg, and David J. Brooks are employees of GE Healthcare plc, which they report as a conflict of interest. Colm J. McGinnity, John S. Duncan, and Matthias J. Koepp have received fees from GE Healthcare plc but are not current or former employees of the organization. John S. Duncan has received fees for organizing symposia and lecturing for UCB Pharma, Eisai, GSK, and GE Healthcare. Alexander Hammers receives license fees for atlas variants that are not used in this study from Imperial Innovations. No other potential conflict of interest relevant to this article was reported.
Acknowledgments
We thank the staff of Hammersmith Imanet Limited and the Epilepsy Society MRI Unit for their assistance with data acquisition and preparation.
Footnotes
Published online Feb. 13, 2014.
- © 2014 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication August 27, 2013.
- Accepted for publication October 4, 2013.