Modeling of receptor ligand data in PET and SPECT imaging: a review of major approaches

A role for foregut tyrosine metabolism in glucose tolerance

Objective

We hypothesized that DA and L-DOPA derived from nutritional tyrosine and the resultant observed postprandial plasma excursions of L-DOPA and DA might affect glucose tolerance via their ability to be taken-up by beta cells and inhibit glucose-stimulated β-cell insulin secretion.

Methods

To investigate a possible circuit between meal-stimulated 3,4-dihydroxy-L-phenylalanine (L-DOPA) and dopamine (DA) production in the GI tract and pancreatic β-cells, we: 1) mapped GI mucosal expression of tyrosine hydroxylase (TH) and aromatic amino acid decarboxylase (AADC); 2) measured L-DOPA and DA content of GI mucosal tissues following meal challenges with different L-tyrosine (TYR) content, 3) determined whether meal TYR content impacts plasma insulin and glucose excursions; and 4) characterized postprandial plasma excursions of L-DOPA and DA in response to meal tyrosine content in rodents and a population of bariatric surgery patients. Next, we characterized: 1) the metabolic transformation of TYR and L-DOPA into DA in vitro using purified islet tissue; 2) the metabolic transformation of orally administrated stable isotope labeled TYR into pancreatic DA, and 3) using a nuclear medicine technique, we studied endocrine beta cells in situ release and binding of DA in response to a glucose challenge.

Results

We demonstrate in rodents that intestinal content and circulatory concentrations L-DOPA and DA, plasma glucose and insulin are responsive to the tyrosine (TYR) content of a test meal. Intestinal expression of two enzymes, Tyrosine hydroxylase (TH) and Aromatic Amino acid Decarboxylase (AADC), essential to the transformation of TYR to DA was mapped and the metabolism of metabolism of TYR to DA was traced in human islets and a rodent beta cell line in vitro and from gut to the pancreas in vivo. Lastly, we show that β cells secrete and bind DA in situ in response to glucose stimulation.

Conclusions

We provide proof-of-principle evidence for the existence of a novel postprandial circuit of glucose homeostasis dependent on nutritional tyrosine. DA and L-DOPA derived from nutritional tyrosine may serve to defend against hypoglycemia via inhibition of glucose-stimulated β-cell insulin secretion as proposed by the anti-incretin hypothesis.

•

Nutritional tyrosine is metabolized to L DOPA and DA in the foregut.

•

Postprandial L-DOPA and DA plasma concentrations rise in response to tyrosine.

•

Oral stable isotope labeled tyrosine is found postprandially in the pancreas as DA.

•

L-DOPA and DA are inhibitors of beta cell glucose-stimulated insulin secretion.

•

Postprandial L-DOPA and DA excursions are muted in certain bariatric surgery patients.

Recent advances in parametric neuroreceptor mapping with dynamic PET: basic concepts and graphical analyses

Tracer kinetic modeling in dynamic positron emission tomography (PET) has been widely used to investigate the characteristic distribution patterns or dysfunctions of neuroreceptors in brain diseases. Its practical goal has progressed from regional data quantification to parametric mapping that produces images of kinetic-model parameters by fully exploiting the spatiotemporal information in dynamic PET data. Graphical analysis (GA) is a major parametric mapping technique that is independent on any compartmental model configuration, robust to noise, and computationally efficient. In this paper, we provide an overview of recent advances in the parametric mapping of neuroreceptor binding based on GA methods. The associated basic concepts in tracer kinetic modeling are presented, including commonly-used compartment models and major parameters of interest. Technical details of GA approaches for reversible and irreversible radioligands are described, considering both plasma input and reference tissue input models. Their statistical properties are discussed in view of parametric imaging.

Evaluation of safety and tolerability, pharmacokinetics, and pharmacodynamics of BMS-820836 in healthy subjects: a placebo-controlled, ascending single-dose study

RATIONALE

BMS-820836, a novel triple monoamine reuptake inhibitor, is an experimental monotherapy for sufferers of major depressive disorder who have had an inadequate response to an existing antidepressant treatment.

OBJECTIVES

This study was conducted to evaluate the safety and tolerability, pharmacokinetics (PK), and serotonin transporter (SERT) and dopamine transporter (DAT) occupancy for single doses of BMS-820836 in healthy subjects.

METHODS

Healthy subjects were assigned to seven BMS-820836 dose panels (0.025, 0.1, 0.5, 1, 2, 3, and 5 mg; n = 8 each), in which subjects were randomly allocated 3:1 to a single BMS-820836 dose or matched placebo. Serial blood samples were collected on Days 1, 2, 3, 4, 7, and 14 to characterize the PK of BMS-820836. Following evaluation of the maximum tolerated dose, SERT occupancy was determined by applying [(11)C]DASB positron emission tomography (PET) after single-dose BMS-820836 (0.5 or 3 mg; n = 3 each) and DAT occupancy by applying [(11)C]PE2I PET after single-dose BMS-820836 (3 mg; n = 6).

RESULTS

Single oral doses of BMS-820836 (0.025-3 mg) were generally safe and well tolerated. BMS-820836 had a median T max of 5.0-7.2 h and a mean apparent terminal T 1/2 of 34-57 h. Mean striatal SERT occupancies were 19 ± 9 % and 82 ± 8 % after single doses of 0.5 and 3 mg BMS-820836, respectively. The mean striatal DAT occupancy was 19 ± 9 % after a single 3 mg BMS-820836 dose.

CONCLUSIONS

Single doses of BMS-820836 have meaningful SERT and DAT occupancy and demonstrate an acceptable safety and tolerability profile in healthy control subjects.

An activity-subspace approach for estimating the integrated input function and relative distribution volume in PET parametric imaging

Dynamic positron emission tomography (PET) imaging technique enables the measurement of neuroreceptor distributions corresponding to anatomic structures, and thus, allows image-wide quantification of physiological and biochemical parameters. Accurate quantification of the concentration of neuroreceptor has been the objective of many research efforts. Compartment modeling is the most widely used approach for receptor binding studies. However, current compartment-model-based methods often either require intrusive collection of accurate arterial blood measurements as the input function, or assume the existence of a reference region. To obviate the need for the input function or a reference region, in this paper, we propose to estimate the input function. We propose a novel concept of activity subspace, and estimate the input function by the analysis of the intersection of the activity subspaces. Then, the input function and the distribution volume (DV) parameter are refined and estimated iteratively. Thus, the underlying parametric image of the total DV is obtained. The proposed method is compared with a blind estimation method, iterative quadratic maximum-likelihood (IQML) via simulation, and the proposed method outperforms IQML. The proposed method is also evaluated in a brain PET dataset.

A Factor-Image Framework to Quantification of Brain Receptor Dynamic PET Studies

The positron emission tomography (PET) imaging technique enables the measurement of receptor distribution or neurotransmitter release in the living brain and the changes of the distribution with time and thus allows quantification of binding sites as well as the affinity of a radioligand. However, quantification of receptor binding studies obtained with PET is complicated by tissue heterogeneity in the sampling image elements (i.e., voxels, pixels). This effect is caused by a limited spatial resolution of the PET scanner. Spatial heterogeneity is often essential in understanding the underlying receptor binding process. Tracer kinetic modeling also often requires an intrusive collection of arterial blood samples. In this paper, we propose a likelihood-based framework in the voxel domain for quantitative imaging with or without the blood sampling of the input function. Radioligand kinetic parameters are estimated together with the input function. The parameters are initialized by a subspace-based algorithm and further refined by an iterative likelihood-based estimation procedure. The performance of the proposed scheme is examined by simulations. The results show that the proposed scheme provides reliable estimation of factor time-activity curves (TACs) and the underlying parametric images. A good match is noted between the result of the proposed approach and that of the Logan plot. Real brain PET data are also examined, and good performance is observed in determining the TACs and the underlying factor images.

The runner's high: opioidergic mechanisms in the human brain

The runner's high describes a euphoric state resulting from long-distance running. The cerebral neurochemical correlates of exercise-induced mood changes have been barely investigated so far. We aimed to unravel the opioidergic mechanisms of the runner's high in the human brain and to identify the relationship to perceived euphoria. We performed a positron emission tomography "ligand activation" study with the nonselective opioidergic ligand 6-O-(2-[(18)F]fluoroethyl)-6-O-desmethyldiprenorphine ([(18)F]FDPN). Ten athletes were scanned at 2 separate occasions in random order, at rest and after 2 h of endurance running (21.5 +/- 4.7 km). Binding kinetics of [(18)F]FDPN were quantified by basis pursuit denoising (DEPICT software). Statistical parametric mapping (SPM2) was used for voxelwise analyses to determine relative changes in ligand binding after running and correlations of opioid binding with euphoria ratings. Reductions in opioid receptor availability were identified preferentially in prefrontal and limbic/paralimbic brain structures. The level of euphoria was significantly increased after running and was inversely correlated with opioid binding in prefrontal/orbitofrontal cortices, the anterior cingulate cortex, bilateral insula, parainsular cortex, and temporoparietal regions. These findings support the "opioid theory" of the runner's high and suggest region-specific effects in frontolimbic brain areas that are involved in the processing of affective states and mood.

Ecstasy-induced reduction of the availability of the brain serotonin transporter as revealed by [11C](+)McN5652-PET and the multi-linear reference tissue model: loss of transporters or artifact of tracer kinetic modelling?

In a previous positron emission tomography (PET) study with the serotonin transporter (SERT) ligand [(11)C](+)McN5652, we found protracted reduction of the availability of the brain SERT in users of the drug ecstasy. However, the multi-linear reference tissue method for the quantification of SERT availability used in this study is prone to effects of altered levels of statistical noise that could simulate reduction of SERT. The aim of the present study was to take into account this confound by re-evaluation of the data now using a modelling approach that is rather insensitive to alterations in the level of statistical noise. A total of 116 subjects (30 current, 29 former ecstasy users, 29 drug-naive, 28 polydrug controls) in whom [(11)C](+)McN5652-PET had been performed previously were re-evaluated. The equilibrium specific-to-non-specific partition coefficient V"( 3) was obtained voxel-wise by application of the simplified reference tissue method (SRTM), which provides quite unbiased results up to rather large noise levels. Voxel-based comparisons between the groups were performed using statistical parametric mapping. V"(3) was reduced in the striatum and in the thalamus in current ecstasy users. This was confirmed by volume-of-interest-based analysis. This result suggests that the ecstasy-induced reduction of SERT availability in SERT-rich brain regions reported previously indicates reduced SERT binding potential rather than being an artifact of tracer kinetic modelling. SRTM analysis did not confirm previous findings in neocortical brain areas.

The tyrosine kinase inhibitor PD153035: implication of labeling position on radiometabolites formed in vitro

INTRODUCTION

The epidermal growth factor receptor is highly expressed in several types of cancers. Molecules with high affinity to its intracellular tyrosine kinase domain are being developed as in vivo imaging probes. The 4-anilinoquinazoline PD153035 has promising in vitro and in vivo properties for development as a reversible radioligand. Labeling it with carbon-11 in either of its two methoxy positions can potentially give rise to different radiometabolites and, consequently, different imaging capabilities. An evaluation of the radiotracers' metabolism was needed to determine the potential significance of the labeling position.

METHODS

PD153035 was labeled in the 6- and 7-O-methoxy positions by reacting the corresponding O-desmethyl precursors with [(11)C]methyl iodide. The two radiolabeled compounds were each incubated for 1 h with human and rat liver microsomes. At five time points, the radiolabeled metabolites were examined using radio-liquid chromatography. One metabolite was isolated and subjected to mass spectroscopic analysis.

RESULTS

A major polar metabolite was obtained in all incubations. Its molecular weight was consistent with an addition of oxygen, and its fragmentation was consistent with an N-oxidation rather than an aromatic hydroxylation. Regioselective 7-O-dealkylation was also observed, albeit in substantial amounts only in the assay using human microsomes.

CONCLUSIONS

Radiolabeling in the 7-O-methoxy position is advocated, since the labeled metabolites produced in the 7-O-demethylation are polar and probably rapidly cleared. The differences observed in the incubations with rat and human microsomes suggest that in vivo positron emission tomography studies with (11)C-labeled PD153035 in rodents may not be directly predictive for studies in humans.

Quantification of adenosine A2A receptors in the human brain using [11C]TMSX and positron emission tomography

PURPOSE

[7-methyl-(11)C]-(E)-8-(3,4,5-trimethoxystyryl)-1,3,7-trimethylxanthine ([(11)C]TMSX) is a positron-emitting adenosine A(2A) receptor (A2AR) antagonist for visualisation of A2AR distribution by positron emission tomography (PET). The aims of this paper were to use a kinetic model to analyse the behaviour of [(11)C]TMSX in the brain and to examine the applicability of the Logan plot. We also studied the applicability of a simplified Logan plot by omitting metabolite correction and arterial blood sampling.

METHODS

The centrum semiovale was used as a reference region on the basis of a post-mortem study showing that it has a negligibly low density of A2ARs. Compartmental analysis was performed in five normal subjects. Parametric images of A2AR binding potential (BP) were also generated using a Logan plot with or without metabolite correction and with or without arterial blood sampling. To omit arterial blood sampling, we applied a method to extract the plasma-related information using independent component analysis (EPICA).

RESULTS

The estimated K (1)/k (2) was confirmed to be common in the centrum semiovale and main cortices. The three-compartment model was well fitted to the other regions using the fixed value of K (1)/k (2) estimated from the centrum semiovale. The estimated BPs using the Logan plot matched those derived from compartment analysis. Without the metabolite correction, the estimate of BP underestimated the true value by 5%. The estimated BPs agreed regardless of arterial blood sampling.

CONCLUSION

A three-compartment model with a reference region, the centrum semiovale, describes the kinetic behaviour of [(11)C]TMSX PET images. A2ARs in the human brain can be visualised as a BP image using [(11)C]TMSX PET without arterial blood sampling.

Serotonin transporter imaging with [123I]beta-CIT SPECT before and after one year of citalopram treatment of obsessive-compulsive disorder

BACKGROUND

Two thirds of patients with obsessive-compulsive disorder (OCD) respond to treatment with selective serotonin reuptake inhibitors (SSRIs). The neurobiological mechanisms of SSRI action and failure to respond to SSRI treatment remain to be elucidated.

OBJECTIVES

The aim of this pilot study was to quantify changes in the availability of serotonin transporter (SERT) in the course of SSRI treatment and to relate these changes to improvements of clinical symptoms.

METHODS

Ten patients with OCD were investigated at baseline and 5 of them after 1 year of SSRI treatment with citalopram 60 mg per day using brain single photon emission computed tomography and [123I]beta-CIT. Specific-to-nondisplaceable [123I]beta-CIT binding ratios (V3'') were calculated in SERT-rich brainstem, midbrain and thalamus using a magnetic resonance imaging-based region of interest (ROI) analysis. Symptom severity was evaluated with the Yale-Brown Obsessive-Compulsive Scale (Y-BOCS).

RESULTS

V3'' differed significantly between pretreatment and follow-up scans in all three brain regions, thalamus, midbrain as well as in brainstem. In thalamic ROI, differ ences in SERT availability and Y-BOCS ratings correlated. In midbrain, a trend toward a significant association was found. In brainstem, no relationship was revealed.

CONCLUSIONS

Higher occupancy of SERT by citalopram seems to be associated with better clinical response after 1 year of SSRI treatment of patients with OCD.

Quantitative analysis of adenosine A1 receptors in human brain using positron emission tomography and [1-methyl-11C]8-dicyclopropylmethyl-1-methyl-3-propylxanthine

Fully quantitative analysis of the adenosine A(1) receptor (A1R) in the brain with (11)C-MPDX and positron emission tomography is reported. The kinetics is described using a two-tissue three-compartment model, and estimated binding potentials correspond well with the estimates made by Logan plot. The image of the binding potential of the MPDX is physiologically reasonable. We conclude that MPDX is applicable to the visualization of the A1Rs in the brain with Logan plot.

Quantification of [18F]diprenorphine kinetics in the human brain with compartmental and non-compartmental modeling approaches

6-O-(2-[(18)F]fluoroethyl)-6-O-desmethyldiprenorphine ([(18)F]FDPN) is a nonselective opiate ligand that binds to postsynaptic micro, kappa and delta opiate receptors. Due to the longer half-life of F-18, compared to C-11, labeling DPN with F-18 allows for alternative experimental protocols and potentially the evaluation of endogenous opioid release. The applicability of this compound to assorted experimental protocols motivated the evaluation of [(18)F]FDPN kinetics with compartmental and non-compartmental models. The results indicate that a two-tissue compartmental model best characterizes the data obtained following a bolus injection of [(18)F]FDPN (120-min scanning protocol). Estimates of distribution volume (DV) were robust, being highly correlated for the one-tissue compartmental model as well as the invasive Logan model and the basis function method. Furthermore, the DV estimates were also stable under a shortened protocol of 60 min, showing a significant correlation with the full protocol. The binding potential (BP) values showed more variability between methods and in some cases were more sensitive to protocol length. In conclusion, this evaluation of [(18)F]FDPN kinetics illustrates that DV values can be estimated robustly using compartmental modeling, the basis function method or the invasive Logan modeling approach on a volume of interest level. BP values were also found to correlate with DV values; however, these results should be interpreted with the understanding that specific binding in the reference region (occipital region) may exist.