Abstract
PET studies of the serotonin (5-hydroxytryptamine, or 5-HT) transporter are increasingly using 11C-3-amino-4-(2-dimethylaminomethylphenylsulfanyl)benzonitrile (DASB). We noted that the percentage of unmetabolized 11C-DASB is often lower at 2 min after injection than at 12 min. We hypothesized that this is due to initial “trapping” of the unmetabolized 11C-DASB compound in the lung, a major 5-HT transporter site and dose-limiting organ. To determine whether binding to an extracranial pool of 5-HT transporters contributes to the lower initial level of unmetabolized 11C-DASB, we examined the effects of sertraline. Methods: Eleven healthy volunteers had 2 11C-DASB PET scans on the same day, and 6 of the 11 had a third scan after sertraline administration. The unmetabolized 11C-DASB fraction was measured in arterial plasma as a function of time and was fit with 2 exponentials with no damping, power function damping, or no damping with the first point removed. Results: Power function damping best fit the data as assessed by visual inspection and residuals and resulted in greater distribution volumes than did no damping with the first point removed. Test–retest reproducibility improved when power function damping was used, as compared with no damping with the first point removed. Oral sertraline raised the 2-min unmetabolized 11C-DASB percentage. Conclusion: Measurement and fitting of early metabolism time points improves curve fitting, significantly affects volume-of-distribution determination, and improves test–retest reproducibility. Saturation of lung 5-HT transporters by sertraline prevents the initial trapping of 11C-DASB. Initial trapping of high-affinity radioligands may be important in the quantification of the binding of other ligands with a high concentration of binding sites in the lungs.
The serotonin (5-hydroxytryptamine, or 5-HT) transporter is the primary target for selective serotonin reuptake inhibitors (SSRIs). SSRIs are an effective treatment for major depressive disorder, obsessive compulsive disorder, posttraumatic stress disorder, generalized anxiety disorder, and other mood and anxiety disorders (1). SSRIs bind to the transporter, blocking the reuptake of 5-HT, and stimulate a chain of intracellular events including desensitization of 5-HT1A autoreceptors, increases in neuronal firing rate, and downstream trophic effects that are temporally associated with clinical response (2,3). In vivo imaging studies report less 5-HT transporter binding in the midbrain (4,5) and amygdala in major depression (5). PET studies have reported on the relationship of 5-HT transporter occupancy in vivo to oral doses and plasma levels of SSRIs (6–8).
11C-3-Amino-4-(2-dimethylaminomethylphenylsulfanyl)benzonitrile (DASB) binds with high affinity (1.1 nmol/L) and selectivity to the 5-HT transporter and offers advantages over other PET 5-HT transporter radioligands (6,9–12). During the course of our determination of the optimal region of interest and voxel-based modeling methods for quantifying 5-HT transporter binding using PET and 11C-DASB (13), we observed that the initial metabolite point (the percentage of unmetabolized 11C-DASB) was lower at the 2-min time point than at the 12-min time point, a phenomenon we have not observed with other radioligands (Fig. 1A) (14,15). Also, a 1-tissue-compartment model did not adequately describe the rising portion of the time–activity curve (Fig. 2A). Because the lung is the critical organ for 11C-DASB (16) and the tracer is injected intravenously, making its first pass through the pulmonary vasculature, one explanation of this phenomenon is that the unmetabolized 11C-DASB is “trapped,” briefly, in the lung.
(A) Average unmetabolized 11C-DASB fraction vs. time for baseline studies. Unlike other ligands, fraction is lower at first time point (2 min) than at second time point (12 min). (B) This lower fraction is caused by lower level of unmetabolized parent compound at 2 min than at 12 min (solid lines) in all but 5 studies (dashed lines).
Time–activity curves for amygdala using ND − 1 (A) and power-function–damped (B) metabolite fits. Time–activity curve fitting using ND − 1metabolite curve fitting is clearly inadequate for early time points. Data are truncated at 60 min; full datasets are shown in insets.
In dogs, the lungs take up at least 90% of plasma 5-HT and are considered a site of 5-HT reuptake from the blood (17). In the human lung, the maximal concentration of 5-HT transporter is 7.50 pmol/mg of protein (18), a concentration that is at least an order of magnitude greater than the maximal concentration of 5-HT transporters in the dorsal raphe nucleus, 329 fmol/mg of tissue (19). In addition, SSRIs such as fluoxetine and venlafaxine have been reported to accumulate in the lungs of mice (20) and humans (21).
In this study, we examined 2 modeling methods to fit this unique unmetabolized fraction curve, and we assessed their impact on the quantification and reproducibility of 11C-DASB binding in the brain. We also tested the effect of prior administration of an SSRI, sertraline, to block the lung sites and determine whether there is an effect on 11C-DASB metabolism.
MATERIALS AND METHODS
Subjects
Healthy volunteers were recruited for the study. Inclusion criteria were assessed by the following: history, Structured Clinical Interview for DSM-IV (22), review of systems, physical examination, routine blood tests, pregnancy test, urine toxicology, and electrocardiography. Inclusion criteria included an age between 18 and 65 y, absence of any axis I DSM-V psychiatric diagnoses, absence of any medications for at least 2 wk, absence of a lifetime history of alcohol or substance abuse or dependence, absence of a lifetime exposure to 3,4-methylenedioxymethamphetamine (“ecstasy”), absence of significant medical conditions, absence of pregnancy, and capacity to provide informed consent. The Institutional Review Boards of Columbia University Medical Center and the New York State Psychiatric Institute approved the protocol. Subjects gave written informed consent after receiving an explanation of the study. Eleven subjects were enrolled in the study (4 men and 7 women; mean age ± SD, 27 ± 6 y; range, 18−36 y). All subjects underwent 2 identical scans, test–retest, on the same day. Six subjects had occupancy scans taken after 2 d of receiving 50 mg/d and then after 4 d of receiving 100 mg/d of sertraline. These subjects were part of a larger study of 5-HT transporter occupancy (8) and model determination (13).
Radiochemistry
11C-DASB was prepared as previously described (23). The chemical purity of 11C-DASB was greater than 99%, and the radiochemical purity was greater than 90%. The mean specific activity at injection was 55.5 MBq/nmol, and the mean injected dose was 570 ± 104 MBq.
PET Protocol
A venous catheter was used for injection of the radiotracer, and an arterial catheter was used to obtain samples for the input function. A polyurethane head-holder system (Soule Medical) was molded around the subject's head for immobilization. PET was performed with an ECAT HR+ (Siemens/CTI). A 10-min transmission scan was obtained before radiotracer injection. At the end of the transmission scan, between 185 and 740 MBq of 11C-DASB were administered intravenously as a bolus over 30 s. Emission data were collected in 3-dimensional mode for 120 min using 21 frames of increasing duration: 3 at 20 s, 3 at 1 min, 3 at 2 min, 2 at 5 min, and 10 at 10 min. Images were reconstructed to a 128 × 128 matrix (pixel size, 2.5 × 2.5 mm). Reconstruction was performed with attenuation correction using the transmission data, and scatter was corrected using a model-based method (24). The reconstruction filter and estimated image filter were Shepp 0.5 (2.5 mm in full width at half maximum), the Z filter was all-pass 0.4 (2.0 mm in full width at half maximum), and the zoom factor was 4.0, leading to a final image resolution of 5.1 mm in full width at half maximum at the center of the field of view (25).
Input Function Measurement
Arterial samples were collected by an automated sampling system every 5 s for the first 2 min and manually thereafter for a total of 31 samples. Centrifuged plasma samples were collected in 200-μL aliquots, and the radioactivity was measured in a well counter. A high-pressure liquid chromatography assay of 6 samples (collected at 2, 12, 20, 50, 80, and 100 min) provided unmetabolized parent and metabolite compound levels. The supernatant liquid obtained after centrifugation of the blood sample at 2,000 rpm for 1 min was transferred (0.5 mL) into a tube and mixed with acetonitrile (0.7 mL). The resulting mixture was stirred with a vortex mixer for 10 s and centrifuged at 14,000 rpm for 4 min. Approximately 0.2 mL of supernatant was removed, and radioactivity was measured in a well counter; 0.8 mL was injected onto the high-pressure liquid chromatography column (Phenomenex Prodigy C18 ODS [10 × 250 mm, 10 μm]; mobile phase, 50:50 acetonitrile:water 0.1 mol/L ammonium formate; flow rate, 5 mL/min; retention time of 11C-DASB, ∼7 min). The metabolite and unmetabolized parent compound were collected in separate vials, and radioactivity was quantified using a Bioscan γ-detector after correction for background radioactivity. To ensure the retention time of the parent fraction, we injected corresponding quality control and standard samples at the beginning and end of the study. We calculated unmetabolized parent fraction as the activity in the parent vial divided by the parent plus metabolite vials. In addition, we measured total parent radioactivity and, to facilitate comparison across subjects, normalized to the total radioactivity injected in the high-pressure liquid chromatography column (Fig. 1B). The 6 unmetabolized parent fraction levels were fit with equations described in the following section. The input function was calculated as the product of the interpolated parent fraction and the total plasma counts. After fitting a straight line between time zero and the peak of the input function, subsequent time points of the measured input function values were fitted to a sum of 3 exponentials, and the fitted values were used as input to the analyses.
Metabolite Curve Fitting
For the first 2 models, the unmetabolized fraction data are fit with a sum of 2 exponentials with coefficients constrained so that the function is 1 at time zero. The smaller of the 2 time constants of the 2 exponentials is constrained to equal the difference between λcerebellum, the terminal rate of washout of cerebellar activity, and λtotal, the smallest elimination rate constant of the total plasma (26):where λcal is the calculated time constant (λcal = λcerebellum − λtotal) and C and λ1 are estimated from the data. This baseline model is referred to as the no-damping (ND) model. The model with the first metabolite point removed is referred to as ND − 1.
Because this model dictates that the unmetabolized portion begins at 1 and decreases monotonically thereafter, it does not adequately account for early trapping of the compound in the lung. We propose to account for this effect by multiplying the baseline model by a term that will damp the early times. We used a power function for this purpose and dropped the constraint that the sum of the coefficients of the exponentials is 1:For small values of t, the power function dominates and f(t) is near zero. As t increases, the exponential portions dominate the power function. Thus, f(t) decays toward zero as t increases, as we would predict. This increases the number of model parameters by 2. We call this the power-function–damped model.
Metabolite models are fit using standard nonlinear least-squares methodology with weights derived from SE values returned from the γ-counter used in metabolite analysis.
MRI
MR images were acquired on a 1.5-T Signa Advantage system (GE Healthcare). A sagittal scout image (localizer) was obtained to identify the anterior commissure–posterior commissure plane (1 min), followed by a transaxial T1-weighted sequence with a 1.5-mm slice thickness in a coronal plane orthogonal to the anterior commissure–posterior commissure plane over the entire brain. The parameters were as follows: 3-dimensional spoiled gradient-recalled acquisition in the steady state; repetition time, 34 ms; echo time, 5 ms; flip angle, 45°; slice thickness, 1.5 mm with no gap; 124 slices; field of view, 22 × 16 cm; matrix, 256 × 192 reformatted to 256 × 256, yielding a voxel size of 1.5 × 0.9 × 0.9 mm; and time of acquisition, 11 min.
Image Analysis
Images were analyzed using MEDX software (Sensor Systems, Inc.). The last 13 frames of an individual study were coregistered to the eighth frame using a tool available from the Oxford Centre for Functional Magnetic Resonance Imaging of the Brain (FMRIB)—the FMRIB linear image registration tool (FLIRT), version 5.0 (27)—to correct for patient motion during the scan. A mean motion-corrected PET image was created and coregistered to its corresponding MR image using FLIRT. The resulting transformation was applied to all motion-corrected frames. Regions of interest were traced on the basis of brain atlases (28,29) and published reports (30,31) and included the anterior cingulate, cingulate cortex, amygdala, hippocampus, insular cortex, temporal cortex, dorsal caudate, midbrain, ventral striatum, thalamus, dorsal putamen, dorsolateral prefrontal cortex, occipital cortex, orbital prefrontal cortex, entorhinal cortex, cerebellar gray matter, medial prefrontal cortex, parietal lobe, parahippocampal gyrus, and posterior parahippocampal gyrus. MR images were cropped using the Exbrain utility (version 2; UCL Institute of Neurology (32)) and then segmented using the FMRIB automated segmentation tool (33). Within the cortex, only gray matter voxels were used to measure the PET activity distribution. Regions of interest were drawn by 3 different data analysts who were unaware of the condition. Test–retest variability in the volume of distribution (VT) between the 3 raters over all regions of interest was less than 3%.
Outcome Measures
Several methods for the quantification of 11C-DASB binding have been described (12,34,35). In this article, 11C-DASB regional VT values were derived using a 1-tissue-compartment model and likelihood estimation in the graphical approach, which reduces the noise-dependent bias inherent in the graphical approach (36–38). The 1-tissue-compartment configuration included the arterial plasma compartment (Ca) and 1 tissue compartment (CT) that included the intracerebral free, nonspecifically bound, and specifically bound compartments. Brain activity was corrected for the contribution of plasma activity assuming a 5% blood volume in the regions of interest (39).
The total regional VT (mL of plasma/g of tissue) was defined as the ratio of the tracer concentration in this region to the metabolite-corrected plasma concentration at equilibrium:Kinetic derivations of total VT were obtained as
where K1 (mL·g−1·min−1) is the unidirectional fractional rate constants for the transfer between Ca and CT, and k2′ (min−1) is the fractional rate constant for the transfer from CT to Ca (14).
Statistics
VT values estimated using the different metabolite models were compared using a linear mixed-effects model on log(VT) with subject and scan (nested within subject) as random effects and region and metabolite model as fixed effects. The log transform was used to approximately equalize the variance among regions.
RESULTS
Generally, when a tracer is administered to a subject as a bolus, the entirety of the compound initially injected is the original, or parent, compound, and the compound is metabolized over time. With 11C-DASB, we observed in 77% of the studies that the first metabolite point, at 2 min, indicated less unmetabolized parent fraction (Fig. 1A) than did the next time point, at 12 min. This finding was due to a lower level of total unmetabolized parent at the 2-min time point in all but 5 studies (Fig. 1B). Fitting the change in percentage of unmetabolized parent compound as a function of time with the sum of 2 exponentials (ND), as is routinely done with this ligand (12) and other 5-HT transporter ligands (14), insufficiently describes this metabolism time course (Fig. 3A). Instead we used 2 alternative methods. In all cases, the terminal portion of the metabolite curve is described by the sum of 2 exponentials where the time constant of the second exponential is constrained. In the alternative methods, we either remove the first time point and fit the remaining points with ND (ND − 1, Fig. 3B) or describe the rising portion of the curve with power function damping (Fig. 3C). Power function damping has the lowest weighted residuals at each time point relative to the ND − 1 (Fig. 4).
Unmetabolized 11C-DASB fraction vs. time data (circles) fit with ND (A), ND with first metabolite point removed (B), and power-function–damped (C) metabolite fits (lines) for study. ND is clearly inadequate for fitting and will not be used in subsequent analyses.
Graph of weighted residuals vs. time for each metabolite fit, jittered on time axis. Metabolite data were fit with ND − 1 (A) and power-function–damped (B) models.
We then explored whether the improved fitting of the metabolite data by the power-function–damped method made any difference in outcome measures. We used a 1-tissue-compartment model to determine the VT and found that without a damping function (ND − 1), VT is underestimated relative to power function damping (F1,838 = 545, P < 0.0001; Fig. 5A) but in very close agreement with other groups (r2 = 0.99 compared with the findings of Frankle et al. (12)). A similar underestimation was observed if we used likelihood estimation in the graphical approach to model the data (F1,838 = 801, P < 0.0001; Fig. 5B). The mean test–retest reproducibility (measured as the difference between the 2 observations divided by the mean across all regions) is also improved by using power function damping (5.9% ± 2.5%) instead of ND − 1 (10.9% ± 9.1%). Similar results were obtained using likelihood estimation in the graphical approach (5.7% ± 1.2% [power function damping] vs. 9.5% ± 10.4% [ND − 1]).
Comparison of VT values with power-function–damped (PFD) metabolite fit against those determined with ND − 1 using 1-tissue-compartment model (A) and likelihood estimation in graphical approach (B). All regions of interest for all studies are plotted. Line is identity. An outlier is excluded from both graphs.
If the large number of 5-HT transporters in the lungs are acting as a trap for the initial intravenous bolus injection of 11C-DASB, then pretreating subjects with high doses of an SSRI such as sertraline could ameliorate this initial trapping. Metabolite curves for subjects in our previous study, who had plasma sertraline levels greater than or equal to 13 ng/mL, are shown in Figure 6. In these subjects, the 2-min time point is lower than the 12-min time point on baseline scan days (Fig. 6A). In scans performed when these subjects were treated with sertraline, the unmetabolized parent fraction at 2 min is higher than the 12-min time point (Fig. 6B; P = 0.015).
Plots of unmetabolized 11C-DASB vs. time for same subjects on baseline scan days (A) and occupancy scan days (B). Subjects had plasma sertraline level ≥ 13 ng/mL on day of occupancy study. On baseline days, most subjects had lower unmetabolized 11C-DASB at 2 min than at 12 min. After subjects received saturating doses of sertraline, pattern was reversed. Subjects are identified with unique markers.
DISCUSSION
We have shown that a power-function–damped 2-exponential fit describes 11C-DASB metabolite data, resulting in both higher VT determinations and improved test–retest reproducibility. An 11C-DASB injection in the presence of high plasma levels of sertraline raises the unmetabolized fraction at the 2-min time point. This finding supports the hypothesis that the large number of 5-HT transporter binding sites in the lung transiently traps the 5-HT transporter radioligand 11C-DASB after it is intravenously injected as a bolus. Conversely, the 11C-DASB metabolite lacks affinity for 5-HT transporter and is not delayed in passage through the lungs. As can be seen in Figure 2 of Lu et al. (16), over 50% of the injected dose of 11C-DASB in humans is in the lungs at initial time points, compared with just slightly over 5% in the organ with the next highest uptake, the heart. Uptake in the brain does not peak until approximately 15 min after injection.
Only one other 11C-DASB study measured metabolites at 2 min (12). Although that study found a larger SD at the 2-min time point than at the other time points, the mean was greater than that at the 12-min time point. That initial publication involved only 6 subjects, and most subjects since then have had lower 2-min metabolite points (Gordon Frankle, oral communication, 2006). All other groups quantifying 5-HT transporter with 11C-DASB either have not had arterial input data or have measured metabolites starting at 5 min. We found that fitting the metabolite data with all data points with ND or fitting the data without the 2-min data point (ND − 1) underestimated VT compared with power function damping (the slope was 0.89 for ND − 1 VT vs. power-function–damped VT and 0.83 for ND VT vs. power-function–damped VT). Although several alternative models for describing the metabolite curve were considered, we could see no room for improvement in the fitting by power function damping and do not present the data. This metabolite pattern has not been observed with other PET ligands in humans, such as 11C-labeled hexahydro-6-[4-(methylthio)phenyl]pyrrolo[2,1-a]isoquinoline (McN5652) (14), 11C-labeled 2-(2-((dimethylamino)methyl)phenylthio)-5-methylaniline (MADAM) (40), and 123I-labeled 2-(2-((dimethylamino)methyl)phenylthio)-5-iodoaniline (ADAM) (41), possibly because of slower initial metabolism or slight differences in affinity.
Although analytic high-pressure liquid chromatography columns are widely used for the analysis of metabolites, we have used a semipreparative column for 11C-DASB as has been done in the past with other compounds (42–44). With an analytic column (Phenomenex Prodigy C18 ODS), a significant difference in the retention time of the parent compound with respect to standard 11C-DASB was observed. With a semipreparative column, the retention time of the parent agreed with that of standard 11C-DASB (7–8 min).
One consequence of using power function damping to fit the metabolite data is that the initial rise of the time–activity curve is better fit by the model. Effectively, this improvement is due to higher K1 values; our range was between 0.25 and 0.95 mL/min/g, compared with 0.62−0.81 mL/min/g as described by Ginovart et al. (34). Because K1 is the product of extraction and blood flow, K1 can never be greater than blood flow. Assuming an extraction of 1, none of our K1 values across any region of interest was above a generally acceptable blood flow value of 1.07 mL/min/g (45).
If correct, our hypothesis regarding initial lung trapping of a radioligand could affect other radioligands for binding sites that have a high concentration of their target in the lungs, such as the 5-HT2A system; our dosimetry studies in baboons suggest that, like the 5-HT transporter, this system has the lung as the critical organ as well (data not shown).
CONCLUSION
Although we have suggested one possible cause for the atypical metabolism of 11C-DASB—selective trapping of the parent compound by the large number of 5-HT transporters in the lungs—we realize other explanations are possible. Regardless of the biologic cause, careful examination and fitting of the rate of metabolism of radioligands can improve the fitting of time–activity curves and the reproducibility of outcome measures and should be applied to previously acquired data as well as to future studies using this ligand.
Acknowledgments
We thank the employees of the Brain Imaging Core of the NIMH-funded Conte Translational Neuroscience Center, the Kreitchman PET Center, and the Columbia University Radioligand Laboratory for expert help. This work was supported in part by an investigator-initiated grant from Pfizer Inc.
Footnotes
-
COPYRIGHT © 2006 by the Society of Nuclear Medicine, Inc.
References
- Received for publication April 5, 2006.
- Accepted for publication August 14, 2006.