|
|
|||||||||
Basic Science Investigations |
1 Department of Clinical Neuroscience, Section of Psychiatry, Karolinska Institutet, Stockholm, Sweden
2 Department of Sensory and Integrative Medicine, Division of Functional Imaging, Niigata University Graduate School of Medicine and Dental Sciences, Asahimachi-dori Niigata, Japan
| ABSTRACT |
|---|
|
|
|---|
Key Words: brain humans serotonin transporter serotonin PET
| INTRODUCTION |
|---|
|
|
|---|
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 (811).
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 |
|---|
|
|
|---|
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 x 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.7411.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 35 sections, as a standardized circular ROI of 6 mm was applied.
TimeActivity 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 timeactivity curves.
Before entering quantitative analysis, the timeactivity curves were corrected for the effect of radioactivity in the cerebral blood volume (CBV) using the following equation:
![]() | (Eq. 1) |
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·mL1·min1) and k2 (min1) correspond to the influx and efflux rates for radioligand diffusion through the bloodbrain barrier, respectively. The rate constants k3 (min1) and k4 (min1) 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 min1). All 4 compartments were assumed homogeneous in concentration and all concentrations have the dimension 37 Bq/mL (nCi/mL).
|
On the basis of this model, the following differential equations can be expressed:
![]() | (Eq. 2) |
![]() | (Eq. 3) |
![]() | (Eq. 4) |
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 timeactivity curves for 11C-MADAM (19). Here K1 corresponds to the influx rate of radioligand diffusion through the bloodbrain 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 momentthat is, the time of transient equilibriumthe 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 timeactivity curve for CB(t) and CN(t) during the time interval 5793 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.
Statistics
Three statistical methods were used to compare the 2 models: the Akaike information criterion (AIC) (28), the Schwarz criterion (29), and F statistics on the squared differences between obtained data and model output (30).
| RESULTS |
|---|
|
|
|---|
|
|
|
In the kinetic analysis, the 2CM and 3CM converged and could describe all regional timeactivity 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 timeactivity 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.
|
|
|
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.9880.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.9600.999 within all regions; P < 0.001; 2-tailed; Table 3).
|
|
|
| DISCUSSION |
|---|
|
|
|---|
|
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 timeactivity 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 timeactivity 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 cerebellumthat 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 timeactivity 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.
|
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 timeactivity 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 23 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,3638). 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 equilibriumhowever, at a rather late time point after injectionwhen 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 |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
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
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. Kalbitzer, C. Svarer, V. G. Frokjaer, D. Erritzoe, W. F.C. Baare, J. Madsen, S. G. Hasselbalch, and G. M. Knudsen A Probabilistic Approach to Delineating Functional Brain Regions J. Nucl. Med. Technol., June 1, 2009; 37(2): 91 - 95. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Nye, J. R. Votaw, N. Jarkas, D. Purselle, V. Camp, J. D. Bremner, C. D. Kilts, C. B. Nemeroff, and M. M. Goodman Compartmental Modeling of 11C-HOMADAM Binding to the Serotonin Transporter in the Healthy Human Brain J. Nucl. Med., December 1, 2008; 49(12): 2018 - 2025. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Okumura, R. Arakawa, H. Ito, C. Seki, H. Takahashi, H. Takano, E. Haneda, R. Nakao, H. Suzuki, K. Suzuki, et al. Quantitative Analysis of NK1 Receptor in the Human Brain Using PET with 18F-FE-SPA-RQ J. Nucl. Med., November 1, 2008; 49(11): 1749 - 1755. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Arakawa, M. Okumura, H. Ito, C. Seki, H. Takahashi, H. Takano, R. Nakao, K. Suzuki, Y. Okubo, C. Halldin, et al. Quantitative Analysis of Norepinephrine Transporter in the Human Brain Using PET with (S,S)-18F-FMeNER-D2 J. Nucl. Med., August 1, 2008; 49(8): 1270 - 1276. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. V. Parsey, A. Ojha, R. T. Ogden, K. Erlandsson, D. Kumar, M. Landgrebe, R. Van Heertum, and J. J. Mann Metabolite Considerations in the In Vivo Quantification of Serotonin Transporters Using 11C-DASB and PET in Humans J. Nucl. Med., November 1, 2006; 47(11): 1796 - 1802. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | RSS | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |