Abstract
11C-N,N-Dimethyl-2-(2-amino-4-methylphenylthio)benzylamine (11C-MADAM) is a newly synthesized radioligand with high selectivity and specificity for the serotonin transporter (5-HTT) in a monkey model. The purpose of this study in humans was to examine the suitability and potential of 11C-MADAM for quantitative PET studies of 5-HTT in applied clinical studies on the pathophysiology and treatment of neuropsychiatric disorders. Methods: PET examination was performed on each of 9 male subjects after intravenous injection of 11C-MADAM with high specific radioactivity. Radioactive metabolites in plasma were determined with high-performance liquid chromatography. A metabolite-corrected arterial input function was used in kinetic 2- and 3-compartment analyses. Cerebellum was used as the reference region in a cross-validation of 6 reference tissue approaches. Results: The highest radioactivity concentration was detected in the raphe nuclei, followed consecutively by the striatum, hippocampal complex, cingulate cortex, neocortex, and cerebellum. The time–activity curve for the fraction of unchanged 11C-MADAM in plasma was best described by a sigmoid function. After 50 min, the fraction was 40%. The labeled metabolites were more polar than the mother compound. The compartment model approaches converged, and could describe the time–activity curves in all regions. The total volume of distribution (Vt) was similar to the regional distribution volumes obtained by the linear graphic analysis. The binding potentials (BPs) for 6 different approaches yielded similar values in all regions but the raphe nuclei, where the 2 equilibrium methods provided lower values. Conclusion: The regional binding distribution of 11C-MADAM is consistent with postmortem data acquired with 3H-MADAM as well as with that of other reference ligands in vitro. The time–activity curves were well described by current major quantitative approaches. The suitability of the cerebellum as a reference region for nonspecific 11C-MADAM binding could be confirmed, thus paving the way for experimentally less demanding approaches, such as the simplified reference tissue model, for applied clinical studies.
The serotonin transporter (5-HTT) is of central interest in the pathophysiology and treatment of psychiatric disorders such as anxiety and depression. This interest has been driven by early findings of decreased 5-HTT densities in suicide victims and in patients with major depression postmortem. Furthermore, the 5-HTT is the primary target for specific serotonin reuptake inhibitors (SSRIs), a group of drugs widely prescribed for the treatment of mood and anxiety disorders (1).
Studies on the pathophysiology and pharmacology of the 5-HTT in humans in vivo has for several years been hampered by the lack of suitable radioligands for PET and SPECT. Of radioligands developed earlier, the diphenyl sulfide derivative 3-amino-4-(2-dimethylaminomethylphenylsulfanyl)-benzonitrile (11C-DASB) (2) has shown the most promising characteristics (3,4).
We have recently reported on N,N-dimethyl-2-(2-amino-4-methylphenylthio)benzylamine (MADAM), a new radioligand with high selectivity and specificity for the 5-HTT (5,6). Initial characterization of 11C-MADAM in monkey brain has shown high specific binding in the brainstem, thalamus, and striatum and lower binding in cortical regions (7). The rank order of specific binding was in agreement with that of human postmortem data (8–11).
The aim of the present PET study on 9 healthy subjects was to examine the potential of 11C-MADAM for quantitative studies of 5-HTT in the human brain. The standard 2- and 3-compartment models (2CM, 3CM) and 6 derived approaches were used and compared in an attempt to select the most appropriate method for applied quantitative studies. The prospects for quantification of 5-HTT binding in the hippocampal complex and the raphe nuclei were given particular attention.
MATERIALS AND METHODS
Subjects and Design
The study was approved by the Ethics and the Radiation Safety Committees of the Karolinska Hospital. Nine male subjects, 22–55 y old, participated after giving informed consent. All subjects were healthy according to history, psychiatric interview, physical examination, blood and urine analysis, and MRI of the brain. They did not use any medication and they were all nonsmokers. Each subject was examined once with PET and 11C-MADAM.
MRI and Head Fixation System
The MRI system used for all subjects was a Signa, 1.5 T (GE Healthcare). T2-weighted and proton density images (156 × 1.0 mm slices) were obtained. To allow the same head positioning in the 2 imaging modalities, a head fixation system with an individual plaster helmet was used both in the PET and the MRI measurements (12).
Radiochemistry
11C-MADAM (for chemical structure, see (13)) was obtained by methylation of 2-((2-((dimethylamino)methyl)phenyl)thio)-5-iodophenylamine (ADAM) using 11C-methyl iodide, as described previously (13). Between 282 and 318 MBq were injected intravenously. The specific radioactivity of the radioligand injected varied between 5.4 and 78 GBq/mmol, corresponding to an injected mass of 0.74–11.9 μg.
PET Experimental Procedure
The PET system used was an ECAT Exact HR 47 (Siemens), which was run in the 3-dimensional mode (14). The in-plane and axial resolution is about 3.8 and 4.0 mm, respectively, full width at half maximum. The reconstructed volume was displayed as 47 sections with a center-to-center distance of 3.125 mm.
In each PET measurement, the subject was placed recumbent with his head in the PET system. One cannula was inserted into the left brachial artery and another was inserted into the right antecubital vein. Sterile physiologic phosphate buffer solution (pH 7.4) containing 11C-MADAM was injected as a bolus over 2 s into the right antecubital vein. The cannula was then immediately flushed with 10 mL saline.
Brain radioactivity was measured in a series of consecutive time frames for 93 min. The frame sequence consisted of three 1-min frames, four 3-min frames, and thirteen 6-min frames. However, for one subject, the measurement lasted for 87 min, resulting in a reduction of one 6-min frame compared with the original frame sequence.
Arterial Blood Sampling
To obtain the arterial input function, an automated blood sampling system was used during the first 5 min of each PET measurement (15). After the first 5 min, arterial blood samples (2 mL) were taken manually at the midpoint of each frame until the end of the measurement (16).
Determination of Radioactive Metabolites in Plasma
The fractions of plasma radioactivity corresponding to unchanged 11C-MADAM and metabolites were determined as has been described previously for other radioligands (17). In brief, arterial plasma samples (2 mL taken at 4, 10, 20, 30, 40, and 50 min) were deproteinized with acetonitrile and analyzed by gradient high-performance liquid chromatography high-performance liquid chromatography (HPLC). Unlabeled MADAM was used as a reference to provide retention times for possible labeled metabolites in plasma.
Regions of Interest (ROIs)
ROIs were defined according to anatomic boundaries for frontal cortex, cingulate cortex, hippocampal complex, putamen, raphe nuclei, and cerebellum. All ROIs except for the raphe nuclei were delineated in 3 consecutive sections on the MR images and transferred to the corresponding reconstructed PET images. On MR images the raphe nuclei cannot be differentiated from surrounding tissue. Therefore, these ROIs were delineated directly on the PET images in 3–5 sections, as a standardized circular ROI of 6 mm was applied.
Time–Activity Curves
To obtain the radioactivity concentration for each region, data for all sectional ROIs of that region were pooled. Regional radioactivity was calculated for each frame, corrected for decay, and plotted versus time, thus providing regional time–activity curves.
Before entering quantitative analysis, the time–activity curves were corrected for the effect of radioactivity in the cerebral blood volume (CBV) using the following equation: Eq. 1 where Cpet(t) represents the radioactivity measured by the PET system, a represents the CBV, and Ca(t) is the time curve for arterial whole blood. a was set to 0.04, the average blood volume in gray matter (16,18,19).
Quantitative Analyses
The 2CM, 3CM, and 6 derived approaches were compared to select the most appropriate method for applied quantitative studies. The binding potential (BP) was the parameter compared.
Kinetic Compartment Analyses.
The data were analyzed using compartment models. The general configuration is a conventional 4-compartment model with 6 rate constants (K1, k2, k3, k4, k5, k6; Fig. 1A). K1 (mL·mL−1·min−1) and k2 (min−1) correspond to the influx and efflux rates for radioligand diffusion through the blood–brain barrier, respectively. The rate constants k3 (min−1) and k4 (min−1) correspond to the rates for radioligand transfer between unbound radioligand in brain (CF, radioactivity concentration of free (unbound) radioligand in brain) and radioligand specifically bound to receptors (CB, radioactivity concentration of radioligand specifically bound to receptors). Radioligand transfer between CF and nonspecifically bound radioligand in brain (CNS, radioactivity concentration of nonspecifically bound radioligand) is described by k5 and k6 (both with the dimension min−1). All 4 compartments were assumed homogeneous in concentration and all concentrations have the dimension 37 Bq/mL (nCi/mL).
A common assumption is that the 2 compartments CF and CNS equilibrate rapidly, thus forming one effective compartment (20), which corresponds to nondisplaceable radioligand in brain (CN). Assuming that CF and CNS equilibrate rapidly, the model can be simplified into 3 compartments and 4 rate constants (Fig. 1A). This model was used to describe the time–activity curves for 11C-MADAM. The radioactivity concentration in plasma was not corrected for plasma protein binding.
On the basis of this model, the following differential equations can be expressed: Eq. 2 Eq. 3 Eq. 4 where CT(t) is the radioactivity concentration in brain, corrected for CBV.
The 2CM (Fig. 1B) is a simplification of the 3CM based on the assumption that all of the compartments, CF, CNS, and CB, equilibrate rapidly to form one effective compartment, CT. The 2CM was used as an alternative approach to describe the time–activity curves for 11C-MADAM (19). Here K1 corresponds to the influx rate of radioligand diffusion through the blood–brain barrier. The rate constant k′2 corresponds to the efflux rate, and its relation to k2, k3, and k4 of the 3CM is given by Equation 5: Eq. 5
In the 3CM analysis the 4 rate constants were determined by curve fitting in a nonlinear least-squares minimization procedure using the Simplex algorithm (21) with constraints for K1 being restricted between 0 and 99.9 and for k2, k3, and k4 between 0 and 9.9. The initial values for K1, k2, k3, and k4 were 0.5, 0.5, 2.0, and 0.5, respectively, and the local minimum of the sum of the squared residuals was determined by an iterative procedure.
Calculation of BP.
Classical receptor binding parameters, such as receptor density (Bmax) and affinity (Kd), cannot be differentiated on the basis of one PET measurement with high specific radioactivity (22). The ratio Bmax/Kd is often referred to as the BP. It corresponds to the ratio k3/k4 in the kinetic analysis. In the present study, BPs calculated directly from k3 and k4 of the 3-compartment analysis are referred to as BPDirect(1).
A variant of a 3-compartment analysis, 3CM(2), was applied with the assumption that CN is identical among all ROIs. In this approach, the K1/k2 ratio was fixed to the ratio for the reference region obtained with the 2CM. The BPs calculated from the thereby acquired k3 and k4 values are referred to as BPDirect(2).
11C-MADAM binding was also expressed using the concept of the total volume of distribution, which for the 3CM is defined as: Eq. 6
The ratio of Vt in the ROI (VtROI) and the Vt in the reference region (VtREF) was entered into the following equation to calculate the BPIndirect: Eq. 7
The BP was also calculated by means of the simplified reference tissue model (BPSRTM) and by the linear graphic analysis for reversible ligand binding to receptors developed by Logan et al., as described in the literature (23,24). The linear graphic analysis provides the regional distribution volume (DV) from which the BP can be calculated (BPLogan).
The ratio CB/CN is equal to the BP when dCB/dt = 0. In this moment—that is, the time of transient equilibrium—the number of molecules associating with the receptors is equal to the number dissociating (16,25,26). This is the basis for the transient equilibrium approach (BPEqA and BPEqB), which also was used to determine the BP as described in the literature. The BP is referred to as BPEqA. A flat shape of CB(t) may make it difficult to define time of peak equilibrium with high reliability. As an alternative approach, the area under the time–activity curve for CB(t) and CN(t) during the time interval 57–93 min was used to calculate BPEqB.
Cerebellum as Reference Region
The cerebellum is a region with very low density of 5-HTT in the human brain in vitro (8,9,11). In a pretreatment experiment with citalopram, 5 mg/kg, in monkey, radioactivity in all examined brain regions was reduced to the level of the cerebellum (7,27). Thus, there is in vitro and in vivo support for a role of the cerebellum as the reference region in calculations of BP.
RESULTS
After injection of 11C-MADAM with high specific activity, radioactivity appeared rapidly in the brain (Fig. 2). The highest radioactivity concentration was found in the raphe nuclei, followed by the putamen, hippocampal complex, cingulate cortex, and frontal cortex. Radioactivity was lowest in the cerebellum (Fig. 3).
The fraction of radioactivity in plasma representing unchanged 11C-MADAM was measured by HPLC. There was a decrease with time and, at 50 min after 11C-MADAM injection, the fraction was about 40% (Fig. 4). At samples drawn after 50 min, the signal-to-noise ratio was not sufficiently high to accurately determine the fractions. The labeled metabolites revealed by the HPLC were more polar than the mother compound.
The plasma radioactivity curves were corrected for radioligand metabolism and plotted versus time (Fig. 3). After an immediate increase, there was a rapid decrease of radioactivity during the first 5 min after radioligand injection. After this time, the decrease continued throughout the experiment but at a slower rate.
In the kinetic analysis, the 2CM and 3CM converged and could describe all regional time–activity curves (Fig. 5). The estimated rate constants are given in Table 1. In 16 of the 63 ROIs studied, the 3CM(1) was statistically superior in describing the time–activity curve, in terms of F statistics as well as the Akaike and Schwarz criteria. The 16 cases in which the 3CM(1) was statistically superior were distributed over all regions including the cerebellum, but not the raphe. A detailed presentation of data obtained for the cerebellum is given in Table 2.
In the second 3-compartment approach, 3CM(2), the ratio of K1/k2 was constrained and defined by the cerebellum. This approach was statistically superior to the 2CM in 16 of the 54 regions (not applicable for the cerebellum).
The total volume of distribution, Vt, obtained by the 3CM(1) and the 3CM(2) showed a significant correlation within all regions (Pearson correlation coefficients, 0.988–0.999; P < 0.001, 2-tailed; Table 3). In the linear graphic analysis, a linear phase was observed for all regions (Fig. 6). The regional distribution volumes (DV) obtained by the linear graphic analysis were correlated with the Vt values obtained by the compartment analyses (Pearson correlation coefficients, 0.960–0.999 within all regions; P < 0.001; 2-tailed; Table 3).
The BP was calculated for the 7 approaches (Table 4). In all analyzed regions, the BPDirect(1) gave 2–6 times higher values than the other methods. For several regions, the SD exceeded the mean value. All other methods yielded similar regional BP values, with the exception of the raphe region, where the equilibrium approaches provided low values when compared with other methods (Table 4).
DISCUSSION
After intravenous injection of 11C-MADAM, there was a marked uptake in all neocortical regions, the striatum, and the raphe nuclei (Figs. 2 and 3). The PET system provided a good signal-to-noise ratio in all regions examined, with the highest ratio in the raphe, where it was higher than 3, and the lowest ratio in the frontal cortex, where it was about 2. The rank order is in good agreement with postmortem studies of radioligand binding to 5-HTT in the human brain (8–11) (Table 5), including a recent study using 3H-MADAM (6). In concurrence with previous characterization of 11C-MADAM in a primate model (7,27), this regional distribution provides support that 11C-MADAM binds specifically to 5-HTT in the human brain in vivo.
The time–activity curve for the fraction of unchanged 11C-MADAM in plasma was best described by a sigmoid function (Fig. 4). Similar observations have been made previously for other radioligands with affinity for 5-HTT and may be explained by high initial uptake of radioligand in lung tissue (31). Of particular interest is that affinity for lung tissue has been described in rodents for the structurally akin radioligand 123I-ADAM (32). A consequence of high initial uptake in lung tissue is a low exposure of radioligand to metabolic enzymes, which thus may explain the low initial metabolic rate. The uptake in lung may also explain the relatively slow distribution of 11C-MADAM to brain (Fig. 3), when compared with other radioligands, such as 11C-raclopride, 11C-SCH 23390, and 11C-WAY 100635 (33,34).
Reversibility of 11C-MADAM binding has been demonstrated in primates by a displacement experiment using citalopram (27). In the present study in humans, reversibility is supported by k4 values above zero and by the linear graphic analysis since the latter part of the curves can be described by a linear function and the slopes do not approach infinity (Fig. 6).
Cerebellum was used as the reference region because this region has been demonstrated to have negligible 5-HTT density in vitro (Table 5). The validity of CN for cerebellum as an index for CN in other regions has been confirmed by experiments in nonhuman primates showing that the time–activity curves for all regions with specific 11C-MADAM binding decreased to the level of the cerebellum after pretreatment with citalopram (5 mg/kg) (7,27).
For the cerebellum, the time–activity curves could be described by the 3CM(1) and, in this respect, there was a high degree of consistency between individuals (Fig. 5; Table 2). The 3CM(1) was statistically superior to the 2CM in only 2 of the 9 subjects (Table 2). This observation may support the model assumption of a rapid equilibrium between the CNS and CF in the cerebellum—that is, for a reference region devoid of specific binding, a collapse of 3 compartments into 2 (20). Clearly, the compartment model analysis does not give support for the existence of kinetically distinguishable nonspecific binding, as has been shown for other radioligands, such as 11C-WAY 100635 (34).
Importantly, despite the observations regarding the validity of the cerebellum as a reference region, calculation of regional CN based on the 3CM(1) data resulted in a curve clearly under the level of the cerebellum time–activity curve (Fig. 7A). This suggests that, although the 3CM(1) does converge, it does not deliver reliable rate constants. As given by the model approach, 3CM(2) provided rate constants that describe CN(t) in accordance with the cerebellum curve (Fig. 7B). This approach gave similar BP values as the 5 other methods, which depends on the cerebellum as the reference tissue.
Macrosystem parameters, such as the distribution volume (K1/k2[1 + k3/k4]), have been acknowledged as being more reliable than individual parameters, such as single rate constants (35). Notably, the distribution volumes based on 3CM(1) and 3CM(2) are very similar to the ones calculated using the linear graphic method (Table 3). The previously mentioned initial accumulation in lung tissue and slow redistribution to brain may yield an unfavorable input function and be a reason for the failure to obtain reliable individual rate constants with 3CM(1). This suggestion is supported by the literature reporting a similar problem in work with other radioligands also belonging to the diphenyl sulfide derivative class, which, accordingly, might share affinity for lung tissue (2,36–38).
The coefficient of variation was highest for the raphe nuclei and, in some approaches, the hippocampal complex (Table 4). These are small regions that also had the highest level of noise (defined by the squared differences between the obtained data and the 3CM(2) model output). Partial-volume effect (PVE) correction was not applied in this study, partly due to the relative homogeneity of the material examined and partly due to the particular attention given to the raphe nuclei, and because it is not possible to define on the MR images, it is not available to PVE correction. Still, because the raphe nuclei is of key interest for the mechanism of action of SSRIs and because the hippocampal complex may have a role in mood disorder, the relatively poor reliability of 5-HTT quantification in these regions has to be taken in account when comparing individual data or small samples.
An analysis based on 3CM(2) with an arterial input function provides a full set of parameters. Nevertheless, in applied studies a simple protocol independent of arterial blood sampling offers practical advantages. Interestingly, several approaches independent of arterial blood sampling gave BP values in accordance with 3CM(2) (Table 4). Among these, the SRTM may be the most advantageous, as it is the only reference method that considers the entire time–activity curve. This is a feature of special importance when analyzing regions with a late equilibrium, such as the raphe nuclei. Thus, the SRTM may be advantageous for clinical studies on regional 5-HTT binding in the human brain.
11C-MADAM belongs to the group of diphenyl sulfide derivatives, as does 11C-DASB, a ligand already applied for studies on 5-HTT occupancy (2,3). This class of molecules has clearly marked an opening in the quest for putative radioligands for quantification of 5-HTT in human brain using PET. One example of this advancement is the signal-to-noise-ratio that has been reported by Frankle et al. to be higher for 11C-DASB than for the previously rather widely used 5-HTT radioligand 11C-McN 5652 (39). Interestingly, a comparison of data from the present study with those of Frankle et al. (39) indicates even higher signal-to-noise-ratios for 11C-MADAM: BPIndirect (based on mean values) is 2–3 times higher in most regions when compared with 11C-DASB. The total volume of distribution in cerebellum, 5.9 ± 1.2 (mean ± SD, n = 9), is about half of that for 11C-DASB (10.1 ± 2.0, n = 6), suggesting low nonspecific binding as a likely explanation for favorably high BP. 11C-MADAM seems to be the only member of this class for which convergence can be reached when applying a standard 3CM to PET data (2,36–38). Reaching convergence is a prerequisite for comparing the distribution volumes of the 3CM(1) and 3CM(2), thus providing support for the BPDirect(2) results. This validation has so far been presented only for 11C-MADAM.
Rapid brain uptake is generally seen as advantageous for a radioligand and particularly advantageous for a peak or ratio analysis. 11C-MADAM and 11C-DASB both satisfy criteria for peak equilibrium—however, at a rather late time point after injection—when compared with established radioligands, such as 11C-raclopride or 11C-WAY 100635 (16,34). Of the 2 radioligands, 11C-MADAM seems to reach equilibrium somewhat earlier than 11C-DASB (Table 6).
CONCLUSION
The regional binding distribution of 11C-MADAM is consistent with postmortem data acquired with 3H-MADAM (6) as well as the distribution demonstrated with other reference ligands in vitro (Table 5). The time–activity curves could be described by current major quantitative approaches, and the suitability of the cerebellum as a reference region for nonspecific 11C-MADAM binding was confirmed, thus paving the way for simplified approaches advantageous for applied clinical studies, such as the SRTM.
Acknowledgments
The excellent technical assistance of Kjerstin Lind is gratefully acknowledged, as is the entire PET group at the Karolinska Institute. This work was supported by the Swedish Research Council (grant 09114).
Footnotes
Received Mar. 8, 2005; revision accepted May 24, 2005.
For correspondence or reprints contact: Johan Lundberg, MD, Psykiatricentrum Karolinska, Karolinska University Hospital Solna, S-171 76 Stockholm, Sweden.
E-mail: johan.lundberg{at}cns.ki.se