Estimation of in vivo nonspecific binding in positron emission tomography studies without requiring a reference region

Exact parameter identification in PET pharmacokinetic modeling using the irreversible two tissue compartment model. Phys Med Biol 2024;69:165008. [PMID: 38830366 PMCID: PMC11288174 DOI: 10.1088/1361-6560/ad539e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Objective.In quantitative dynamic positron emission tomography (PET), time series of images, reflecting the tissue response to the arterial tracer supply, are reconstructed. This response is described by kinetic parameters, which are commonly determined on basis of the tracer concentration in tissue and the arterial input function. In clinical routine the latter is estimated by arterial blood sampling and analysis, which is a challenging process and thus, attempted to be derived directly from reconstructed PET images. However, a mathematical analysis about the necessity of measurements of the common arterial whole blood activity concentration, and the concentration of free non-metabolized tracer in the arterial plasma, for a successful kinetic parameter identification does not exist. Here we aim to address this problem mathematically.Approach.We consider the identification problem in simultaneous pharmacokinetic modeling of multiple regions of interests of dynamic PET data using the irreversible two-tissue compartment model analytically. In addition to this consideration, the situation of noisy measurements is addressed using Tikhonov regularization. Furthermore, numerical simulations with a regularization approach are carried out to illustrate the analytical results in a synthetic application example.Main results.We provide mathematical proofs showing that, under reasonable assumptions, all metabolic tissue parameters can be uniquely identified without requiring additional blood samples to measure the arterial input function. A connection to noisy measurement data is made via a consistency result, showing that exact reconstruction of the ground-truth tissue parameters is stably maintained in the vanishing noise limit. Furthermore, our numerical experiments suggest that an approximate reconstruction of kinetic parameters according to our analytic results is also possible in practice for moderate noise levels.Significance.The analytical result, which holds in the idealized, noiseless scenario, suggests that for irreversible tracers, fully quantitative dynamic PET imaging is in principle possible without costly arterial blood sampling and metabolite analysis.

Development of a minimally invasive simultaneous estimation method for quantifying translocator protein binding with [18F]FEPPA positron emission tomography

BACKGROUND

The purpose of this study was to assess the feasibility of using a minimally invasive simultaneous estimation method (SIME) to quantify the binding of the 18-kDa translocator protein tracer [¹⁸F]FEPPA. Arterial sampling was avoided by extracting an image-derived input function (IDIF) that was metabolite-corrected using venous blood samples. The possibility of reducing scan duration to 90 min from the recommended 2-3 h was investigated by assuming a uniform non-displaceable distribution volume (V_ND) to simplify the SIME fitting.

RESULTS

SIME was applied to retrospective data from healthy volunteers and was comprised of both high-affinity binders (HABs) and mixed-affinity binders (MABs). Estimates of global V_ND and regional total distribution volume (V_T) from SIME were not significantly different from values obtained using a two-tissue compartment model (2CTM). Regional V_T estimates were greater for HABs compared to MABs for both the 2TCM and SIME, while the SIME estimates had lower inter-subject variability (41 ± 17% reduction). Binding potential (BP_ND) values calculated from regional V_T and brain-wide V_ND estimates were also greater for HABs, and reducing the scan time from 120 to 90 min had no significant effect on BP_ND. The feasibility of using venous metabolite correction was evaluated in a large animal model involving a simultaneous collection of arterial and venous samples. Strong linear correlations were found between venous and arterial measurements of the blood-to-plasma ratio and the remaining [¹⁸F]FEPPA fraction. Lastly, estimates of BP_ND and the specific distribution volume (i.e., V_S = V_T - V_ND) from a separate group of healthy volunteers (90 min scan time, venous-scaled IDIFs) agreed with estimates from the retrospective data for both genotypes.

CONCLUSIONS

The results of this study demonstrate that accurate estimates of regional V_T, BP_ND and V_S can be obtained by applying SIME to [¹⁸F]FEPPA data. Furthermore, the application of SIME enabled the scan time to be reduced to 90 min, and the approach worked well with IDIFs that were scaled and metabolite-corrected using venous blood samples.

Simultaneous multifactor Bayesian analysis (SiMBA) of PET time activity curve data

Positron emission tomography (PET) is an in vivo imaging method essential for studying the neurochemical pathophysiology of psychiatric and neurological disease. However, its high cost and exposure of participants to radiation make it unfeasible to employ large sample sizes. The major shortcoming of PET imaging is therefore its lack of power for studying clinically-relevant research questions. Here, we introduce a new method for performing PET quantification and analysis called SiMBA, which helps to alleviate these issues by improving the efficiency of PET analysis by exploiting similarities between both individuals and regions within individuals. In simulated [11C]WAY100635 data, SiMBA greatly improves both statistical power and the consistency of effect size estimation without affecting the false positive rate. This approach makes use of hierarchical, multifactor, multivariate Bayesian modelling to effectively borrow strength across the whole dataset to improve stability and robustness to measurement error. In so doing, parameter identifiability and estimation are improved, without sacrificing model interpretability. This comes at the cost of increased computational overhead, however this is practically negligible relative to the time taken to collect PET data. This method has the potential to make it possible to test clinically-relevant hypotheses which could never be studied before given the practical constraints. Furthermore, because this method does not require any additional information over and above that required for traditional analysis, it makes it possible to re-examine data which has already previously been collected at great expense. In the absence of dramatic advancements in PET image data quality, radiotracer development, or data sharing, PET imaging has been fundamentally limited in the scope of research hypotheses which could be studied. This method, especially combined with the recent steps taken by the PET imaging community to embrace data sharing, will make it possible to greatly improve the research possibilities and clinical relevance of PET neuroimaging.

Kinetic modeling and parameter estimation of TSPO PET imaging in the human brain

PURPOSE

Translocator protein 18-kDa (TSPO) imaging with positron emission tomography (PET) is widely used in research studies of brain diseases that have a neuro-immune component. Quantification of TSPO PET images, however, is associated with several challenges, such as the lack of a reference region, a genetic polymorphism affecting the affinity of the ligand for TSPO, and a strong TSPO signal in the endothelium of the brain vessels. These challenges have created an ongoing debate in the field about which type of quantification is most useful and whether there is an appropriate simplified model.

METHODS

This review focuses on the quantification of TSPO radioligands in the human brain. The various methods of quantification are summarized, including the gold standard of compartmental modeling with metabolite-corrected input function as well as various alternative models and non-invasive approaches. Their advantages and drawbacks are critically assessed.

RESULTS AND CONCLUSIONS

Researchers employing quantification methods for TSPO should understand the advantages and limitations associated with each method. Suggestions are given to help researchers choose between these viable alternative methods.

Specific and non-specific binding of a tracer for the translocator-specific protein in schizophrenia: an [11C]-PBR28 blocking study

OBJECTIVE

The mitochondrial 18-kDa translocator protein (TSPO) is expressed by activated microglia and positron emission tomography enables the measurement of TSPO levels in the brain. Findings in schizophrenia have shown to vary depending on the outcome measure used and this discrepancy in TSPO results could be explained by lower non-displaceable binding (VND) in schizophrenia, which could obscure increases in specific binding. In this study, we have used the TSPO ligand XBD173 to block the TSPO radioligand [11C]-PBR28 and used an occupancy plot to quantify VND in patients with schizophrenia.

METHODS

A total of 7 patients with a diagnosis of schizophrenia were recruited for this study. Each patient received two separate PET scans with [11C]PBR28, one at baseline and one after the administration of the TSPO ligand XBD173. All patients were high-affinity binders (HABs) for the TSPO gene. We used an occupancy plot to quantify the non-displaceable component (VND) using 2TCM kinetic estimates with and without vascular correction. Finally we computed the VND at a single subject level using the SIME method.

RESULTS

All patients showed a global and generalized reduction in [11C]PBR28 uptake after the administration of XBD173. Constraining the VND to be equal for all patients, the population VND was estimated to be 1.99 mL/cm3 (95% CI 1.90 to 2.08). When we used vascular correction, the fractional TSPO occupancy remained similar.

CONCLUSIONS

In schizophrenia patients, a substantial component of the [11C]PBR28 signal represents specific binding to TSPO. Furthermore, the VND in patients with schizophrenia is similar to that previously reported in healthy controls. These results suggest that changes in non-specific binding between schizophrenia patients and healthy controls do not account for discrepant PET findings in this disorder.

Test-retest variability and reference region-based quantification of 18F-BCPP-EF for imaging mitochondrial complex I in the human brain

Mitochondrial complex I (MC-I) is an essential regulator of brain bioenergetics and can be quantified in the brain using PET radioligand ¹⁸F-BCPP-EF. Here we evaluate the test-retest reproducibility of ¹⁸F-BCPP-EF in humans, and assess the use of a non-invasive quantification method (standardised uptake value ratio - SUVR). Thirty healthy volunteers had a 90-min dynamic ¹⁸F-BCPP-EF scan with arterial blood sampling, five of which received a second scan to be included in the test-retest analysis. Time-activity curves (TAC) were analysed using multilinear analysis 1 (MA1) and the two-tissue compartment model (2TC) to estimate volumes of distribution (V_T). Regional SUVR-1 values were calculated from the 70 to 90-min TAC data using the centrum semiovale as a pseudo reference region, and compared to kinetic analysis-derived outcome measures. The mean absolute test-retest variability of V_T ranged from 12% to 18% across regions. Both DVR-1and SUVR-1 had improved test-retest variability in the range 2%-7%. SUVR-1 was highly correlated with DVR-1 (r² = 0.97, n = 30). In conclusion, ¹⁸F-BCPP-EF has suitable test-retest reproducibility and can be used to quantify MC-I in clinical studies.

Nondisplaceable Binding Is a Potential Confounding Factor in 11C-PBR28 Translocator Protein PET Studies

The PET ligand ¹¹C-PBR28 (N-((2-(methoxy-¹¹C)-phenyl)methyl)-N-(6-phenoxy-3-pyridinyl)acetamide) binds to the 18-kDa translocator protein (TSPO), a biomarker of glia. In clinical studies of TSPO, the ligand total distribution volume, V_T, is frequently the reported outcome measure. Since V_T is the sum of the ligand-specific distribution volume (V_S) and the nondisplaceable-binding distribution volume (V_ND), differences in V_ND across subjects and groups will have an impact on V_T Methods: Here, we used a recently developed method for simultaneous estimation of V_ND (SIME) to disentangle contributions from V_ND and V_S Data from 4 previously published ¹¹C-PBR28 PET studies were included: before and after a lipopolysaccharide challenge (8 subjects), in alcohol use disorder (14 patients, 15 controls), in first-episode psychosis (16 patients, 16 controls), and in Parkinson disease (16 patients, 16 controls). In each dataset, regional V_T estimates were obtained with a standard 2-tissue-compartment model, and brain-wide V_ND was estimated with SIME. V_S was then calculated as V_T - V_ND V_ND and V_S were then compared across groups, within each dataset. Results: A lower V_ND was found for individuals with alcohol-use disorder (34%, P = 0.00084) and Parkinson disease (34%, P = 0.0032) than in their corresponding controls. We found no difference in V_ND between first-episode psychosis patients and their controls, and the administration of lipopolysaccharide did not change V_ND Conclusion: Our findings suggest that in TSPO PET studies, nondisplaceable binding can differ between patient groups and conditions and should therefore be considered.

Quantification of Positron Emission Tomography Data Using Simultaneous Estimation of the Input Function: Validation with Venous Blood and Replication of Clinical Studies

PURPOSE

To determine if one venous blood sample can substitute full arterial sampling in quantitative modeling for multiple positron emission tomography (PET) radiotracers using simultaneous estimation of the input function (SIME).

PROCEDURES

Participants underwent PET imaging with [¹¹C]ABP688, [¹¹C]CUMI-101, and [¹¹C]DASB. Full arterial sampling and additional venous blood draws were performed for quantification with the arterial input function (AIF) and SIME using one arterial or venous (vSIME) sample.

RESULTS

Venous and arterial metabolite-corrected plasma activities were within 6 % of each other at varying time points. vSIME- and AIF-derived outcome measures were in good agreement, with optimal sampling times of 12 min ([¹¹C]ABP688), 90 min ([¹¹C]CUMI-101), and 100 min ([¹¹C]DASB). Simulation-based power analyses revealed that SIME required fewer subjects than the AIF method to achieve statistical power, with significant reductions for [¹¹C]CUMI-101 and [¹¹C]DASB with vSIME. Replication of previous findings and test-retest analyses bolstered the simulation analyses.

CONCLUSIONS

We demonstrate the feasibility of AIF recovery using SIME with one venous sample for [¹¹C]ABP688, [¹¹C]CUMI-101, and [¹¹C]DASB. This method simplifies PET acquisition while allowing for fully quantitative modeling, although some variability and bias are present with respect to AIF-based quantification, which may depend on the accuracy of the single venous blood measurement.

Quantification of [11C]PBR28 data after systemic lipopolysaccharide challenge

BACKGROUND

Lipopolysaccharide (LPS) is a classic immune stimulus. LPS combined with positron emission tomography (PET) 18 kDa translocator protein (TSPO) brain imaging provides a robust human laboratory model to study neuroimmune signaling. To evaluate optimal analysis of these data, this work compared the sensitivity of six quantification approaches.

METHODS

[11C]PBR28 data from healthy volunteers (N = 8) were collected before and 3 h after LPS challenge (1.0 ng/kg IV). Quantification approaches included total volume of distribution estimated with two tissue, and two tissue plus irreversible uptake in whole blood, compartment models (2TCM and 2TCM-1k, respectively) and multilinear analysis-1 (MA-1); binding potential estimated with simultaneous estimation (SIME); standardized uptake values (SUV); and SUV ratio (SUVR).

RESULTS

The 2TCM, 2TCM-1k, MA-1, and SIME approaches each yielded substantive effect sizes for LPS effects (partial η2 = 0.56-0.89, ps <. 05), whereas SUV and SUVR did not.

CONCLUSION

These findings highlight the importance of incorporating AIF measurements to quantify LPS-TSPO studies.

Brain serotonin transporter binding, plasma arachidonic acid and depression severity: A positron emission tomography study of major depression

BACKGROUND

Serotonin transporter (5-HTT) binding and polyunsaturated fatty acids (PUFAs) are implicated in major depressive disorder (MDD). Links between the two systems in animal models have not been investigated in humans.

METHODS

Using positron emission tomography (PET) and [¹¹C]DASB, we studied relationships between 5-HTT binding potential and plasma levels of PUFAs docosahexaenoic acid (DHA), eicosapentaenoic acid (EPA), and arachidonic acid (AA) in medication-free MDD patients (n = 21). PUFAs were quantified using transesterification and gas chromatography. Binding potential BP_P, and alternative outcome measures BP_F and BP_ND, were determined for [¹¹C]DASB in six a priori brain regions of interest (ROIs) using likelihood estimation in graphical analysis (LEGA) to calculate radioligand total distribution volume (V_T), and a validated hybrid deconvolution approach (HYDECA) that estimates radioligand non-displaceable distribution volume (V_ND) without a reference region. Linear mixed models used PUFA levels as predictors and binding potential measures as outcomes across the specified ROIs; age and sex as fixed effects; and subject as random effect to account for across-region binding correlations. As nonlinear relationships were observed, a quadratic term was added to final models.

RESULTS

AA predicted both 5-HTT BP_P and depression severity nonlinearly, described by an inverted U-shaped curve. 5-HTT binding potential mediated the relationship between AA and depression severity.

LIMITATIONS

Given the small sample and multiple comparisons, results require replication.

CONCLUSIONS

Our findings suggest that AA status may impact depression pathophysiology through effects on serotonin transport. Future studies should examine whether these relationships explain therapeutic effects of PUFAs in the treatment of MDD.

Accuracy and reliability of [11C]PBR28 specific binding estimated without the use of a reference region

[¹¹C]PBR28 is a positron emission tomography radioligand used to examine the expression of the 18 kDa translocator protein (TSPO). TSPO is located in glial cells and can function as a marker for immune activation. Since TSPO is expressed throughout the brain, no true reference region exists. For this reason, an arterial input function is required for accurate quantification of [¹¹C]PBR28 binding and the most common outcome measure is the total distribution volume (V_T). Notably, V_T reflects both specific binding and non-displaceable binding. Therefore, estimates of specific binding, such as binding potential (e.g. BP_ND) and specific distribution volume (V_S) should theoretically be more sensitive to underlying differences in TSPO expression. It is unknown, however, if unbiased and accurate estimates of these outcome measures are obtainable for [¹¹C]PBR28. The Simultaneous Estimation (SIME) method uses time-activity-curves from multiple brain regions with the aim to obtain a brain-wide estimate of the non-displaceable distribution volume (V_ND), which can subsequently be used to improve the estimation of BP_ND and V_S. In this study we evaluated the accuracy of SIME-derived V_ND, and the reliability of resulting estimates of specific binding for [¹¹C]PBR28, using a combination of simulation experiments and in vivo studies in healthy humans. The simulation experiments, based on data from 54 unique [¹¹C]PBR28 examinations, showed that V_ND values estimated using SIME were both precise and accurate. Data from a pharmacological competition challenge (n = 5) showed that SIME provided V_ND values that were on average 19% lower than those obtained using the Lassen plot, but similar to values obtained using the Likelihood-Estimation of Occupancy technique. Test-retest data (n = 11) showed that SIME-derived V_S values exhibited good reliability and precision, while larger variability was observed in SIME-derived BP_ND values. The results support the use of SIME for quantifying specific binding of [¹¹C]PBR28, and suggest that V_S can be used in complement to the conventional outcome measure V_T. Additional studies in patient cohorts are warranted.

Molecular imaging of serotonin degeneration in mild cognitive impairment

Neuropathological and neuroimaging studies have consistently demonstrated degeneration of monoamine systems, especially the serotonin system, in normal aging and Alzheimer's disease. The evidence for degeneration of the serotonin system in mild cognitive impairment is limited. Thus, the goal of the present study was to measure the serotonin transporter in vivo in mild cognitive impairment and healthy controls. The serotonin transporter is a selective marker of serotonin terminals and of the integrity of serotonin projections to cortical, subcortical and limbic regions and is found in high concentrations in the serotonergic cell bodies of origin of these projections (raphe nuclei). Twenty-eight participants with mild cognitive impairment (age 66.6±6.9, 16 males) and 28 healthy, cognitively normal, demographically matched controls (age 66.2±7.1, 15 males) underwent magnetic resonance imaging for measurement of grey matter volumes and high-resolution positron emission tomography with well-established radiotracers for the serotonin transporter and regional cerebral blood flow. Beta-amyloid imaging was performed to evaluate, in combination with the neuropsychological testing, the likelihood of subsequent cognitive decline in the participants with mild cognitive impairment. The following hypotheses were tested: 1) the serotonin transporter would be lower in mild cognitive impairment compared to controls in cortical and limbic regions, 2) in mild cognitive impairment relative to controls, the serotonin transporter would be lower to a greater extent and observed in a more widespread pattern than lower grey matter volumes or lower regional cerebral blood flow and 3) lower cortical and limbic serotonin transporters would be correlated with greater deficits in auditory-verbal and visual-spatial memory in mild cognitive impairment, not in controls. Reduced serotonin transporter availability was observed in mild cognitive impairment compared to controls in cortical and limbic areas typically affected by Alzheimer's disease pathology, as well as in sensory and motor areas, striatum and thalamus that are relatively spared in Alzheimer's disease. The reduction of the serotonin transporter in mild cognitive impairment was greater than grey matter atrophy or reductions in regional cerebral blood flow compared to controls. Lower cortical serotonin transporters were associated with worse performance on tests of auditory-verbal and visual-spatial memory in mild cognitive impairment, not in controls. The serotonin system may represent an important target for prevention and treatment of MCI, particularly the post-synaptic receptors (5-HT4 and 5-HT6), which may not be as severely affected as presynaptic aspects of the serotonin system, as indicated by the observation of lower serotonin transporters in MCI relative to healthy controls.

A hybrid deconvolution approach for estimation of in vivo non-displaceable binding for brain PET targets without a reference region

BACKGROUND AND AIM

Estimation of a PET tracer's non-displaceable distribution volume (VND) is required for quantification of specific binding to its target of interest. VND is generally assumed to be comparable brain-wide and is determined either from a reference region devoid of the target, often not available for many tracers and targets, or by imaging each subject before and after blocking the target with another molecule that has high affinity for the target, which is cumbersome and involves additional radiation exposure. Here we propose, and validate for the tracers [11C]DASB and [11C]CUMI-101, a new data-driven hybrid deconvolution approach (HYDECA) that determines VND at the individual level without requiring either a reference region or a blocking study.

METHODS

HYDECA requires the tracer metabolite-corrected concentration curve in blood plasma and uses a singular value decomposition to estimate the impulse response function across several brain regions from measured time activity curves. HYDECA decomposes each region's impulse response function into the sum of a parametric non-displaceable component, which is a function of VND, assumed common across regions, and a nonparametric specific component. These two components differentially contribute to each impulse response function. Different regions show different contributions of the two components, and HYDECA examines data across regions to find a suitable common VND. HYDECA implementation requires determination of two tuning parameters, and we propose two strategies for objectively selecting these parameters for a given tracer: using data from blocking studies, and realistic simulations of the tracer. Using available test-retest data, we compare HYDECA estimates of VND and binding potentials to those obtained based on VND estimated using a purported reference region.

RESULTS

For [11C]DASB and [11C]CUMI-101, we find that regardless of the strategy used to optimize the tuning parameters, HYDECA provides considerably less biased estimates of VND than those obtained, as is commonly done, using a non-ideal reference region. HYDECA test-retest reproducibility is comparable to that obtained using a VND determined from a non-ideal reference region, when considering the binding potentials BPP and BPND.

CONCLUSIONS

HYDECA can provide subject-specific estimates of VND without requiring a blocking study for tracers and targets for which a valid reference region does not exist.

Positron emission tomography quantification of serotonin transporter binding in medication-free bipolar disorder

OBJECTIVES

Bipolar disorder (BD) is associated with abnormalities in the serotonin transporter (5-HTT), but specific in vivo findings have been discrepant. Using positron emission tomography (PET) and [(11)C]DASB, we compared 5-HTT binding between unmedicated depressed BD subjects and healthy volunteers (HVs).

EXPERIMENTAL DESIGN

5-HTT binding in six brain regions was compared between 17 depressed, unmedicated BD subjects and 31 HVs, using the outcome measure of VT/fP (proportional to the total number of available transporters). Alternative outcome measures were examined as well. 47% of BD were BP I; and 65% reported a prior suicide attempt.

PRINCIPAL OBSERVATIONS

5-HTT binding (VT/fP ) did not differ between BD and HV groups considering six brain regions of interest simultaneously (P = 0.24). In contrast, alternative outcome measures (BPF*, BPP*, and BPND*) indicated lower binding in BD compared with HV across these six regions of interest (BPF*: P = 0.047; BPP*: P = 0.032; BPND*: P = 0.031). 5-HTT binding was unrelated to suicide attempt history, depression severity, bipolar subtype, or history of past substance use disorder.

CONCLUSIONS

Choice of outcome measure strongly affects comparisons of serotonin transporter binding using PET with [(11)C]DASB. We do not find evidence of abnormal 5-HTT binding in bipolar depression using our primary outcome measure, VT /fP . However, we did observe lower 5-HTT binding in BD with alternative outcome measures that are frequently used with [(11)C]DASB. Relative merits and assumptions of different outcome measures are discussed. Evaluation in larger samples and during different mood states, including remission, is warranted.