Uncertainty Analysis of Time-Integrated Activity Coefficient in Single-Time-Point Dosimetry Using Bayesian Fitting Method

Purpose

Calculation of the uncertainty of the individual time-integrated activity coefficient (TIACs) is desirable in molecular radiotherapy. However, the calculation of TIAC's uncertainty in single-time-point (STP) method has never been reported in the literature. This study presents a method based on the Bayesian fitting (BF) to calculate the standard deviation (SD) of individual TIACs in the STP dosimetry.

Methods

Biokinetic data of 177Lu-DOTATATE in kidneys were obtained from PMID33443063. BF methods with extended objective function, which optimize the fitting using prior knowledge of the function's parameters, were used. Reference TIACs (rTIACs) were calculated by fitting a mono-exponential function to the all-time-point data. The goodness of fit was checked based on the visual inspection and the coefficient of variations (CV) of the fitted parameters < 0.5. BF with relative (BFr) and absolute-based (BFa) variance methods were used to obtain the calculated TIACs (cTIACs) from the STP dosimetry. Performance of the STP method was obtained by calculating the relative deviation (RD) between cTIACs and rTIACs.

Results

Visual inspection showed a good fit for all patients with CV of fitted parameters less than 50%. The mean ± SD of cTIAC's %RD were 7.0 ± 25.2 for BFr and 2.6 ± 8.9 for BFa. The range of %CV of the individual cTIAC's SD for BFr and BFa methods was 36-78% and 22-33%, respectively, while the %CV of the rTIAC SD was 0.8-49%.

Conclusion

We introduce the BF method to calculate the SD of individual TIACs in STP dosimetry. The presented method might be used as an alternative method for uncertainty analysis in STP dosimetry.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13139-024-00851-8.

Biodistribution and radiation dosimetry of [99mTc]Tc-N4-BTG in patients with biochemical recurrence of prostate cancer

BACKGROUND

In patients with prostate cancer (PCa), imaging with gastrin-releasing peptide receptor (GRPR) ligands is an alternative to PSMA-targeted tracers, particularly if PSMA expression is low or absent. [99mTc]Tc-N4-BTG is a newly developed GRPR-directed probe for conventional scintigraphy and single photon emission computed tomography (SPECT) imaging. The current study aims to investigate the safety, biodistribution and dosimetry of [99mTc]Tc-N4-BTG in patients with biochemical recurrence (BCR) of PCa.

RESULTS

No adverse pharmacologic effects were observed. Injection of [99mTc]Tc-N4-BTG resulted in an effective dose of 0.0027 ± 0.0002 mSv/MBq. The urinary bladder was the critical organ with the highest mean absorbed dose of 0.028 ± 0.001 mGy/MBq, followed by the pancreas with 0.0043 ± 0.0015 mGy/MBq, osteogenic cells with 0.0039 ± 0.0005 mGy/MBq, the kidneys with 0.0034 ± 0.0003 mGy/MBq, and the liver with 0.0019 ± 0.0004 mGy/MBq, respectively. No focal tracer uptake suggestive of PCa recurrence could be revealed for any of the patients.

CONCLUSION

[99mTc]Tc-N4-BTG appears to be a safe diagnostic agent. Compared to GRPR-targeted PET tracers, this 99mTc-labelled SPECT agent could contribute to a broader application and better availability of this novel approach. Further research to assess its clinical value is warranted.

EANM guidance document: dosimetry for first-in-human studies and early phase clinical trials

The numbers of diagnostic and therapeutic nuclear medicine agents under investigation are rapidly increasing. Both novel emitters and novel carrier molecules require careful selection of measurement procedures. This document provides guidance relevant to dosimetry for first-in human and early phase clinical trials of such novel agents. The guideline includes a short introduction to different emitters and carrier molecules, followed by recommendations on the methods for activity measurement, pharmacokinetic analyses, as well as absorbed dose calculations and uncertainty analyses. The optimal use of preclinical information and studies involving diagnostic analogues is discussed. Good practice reporting is emphasised, and relevant dosimetry parameters and method descriptions to be included are listed. Three examples of first-in-human dosimetry studies, both for diagnostic tracers and radionuclide therapies, are given.

Single-Time-Point Renal Dosimetry Using Nonlinear Mixed-Effects Modeling and Population-Based Model Selection in [177Lu]Lu-PSMA-617 Therapy

The aim of this study was to investigate the accuracy of single-time-point (STP) renal dosimetry imaging using SPECT/CT data, a nonlinear mixed-effects (NLME) model, and a population-based model selection (PBMS) in a large population for 177Lu-labeled prostate-specific membrane antigen therapy. Methods: Biokinetic data (mean ± SD) of [177Lu]Lu-PSMA-617 in kidneys at time points 1 (1.8 ± 0.8 h), 2 (18.7 ± 0.9 h), 3 (42.6 ± 1.0 h), 4 (66.3 ± 0.9 h), and 5 (160.3 ± 24.2 h) after injection were obtained from 63 patients with metastatic castration-resistant prostate cancer using SPECT/CT. Thirteen functions were derived from various parameterizations of 1- to 5-exponential functions. The function's parameters were fitted in the NLME framework to the all-time-point (ATP) data. The PBMS NLME method was performed using the goodness-of-fit test and Akaike weight to select the best function fitting the data. The best function from ATP fitting was used to calculate the reference time-integrated activity and absorbed doses. In STP dosimetry, the parameters of a particular patient with STP data were fitted simultaneously to the STP data at different time points of that patient with ATP data of all other patients. The parameters from STP fitting were used to calculate the STP time-integrated activity and absorbed doses. Relative deviations (RDs) and root-mean-square errors (RMSEs) were used to analyze the accuracy of the calculated STP absorbed dose compared with the reference absorbed dose obtained from the best-fit ATP function. The performance of STP dosimetry using PBMS NLME modeling was compared with the Hänscheid and Madsen methods. Results: The function [Formula: see text] was selected as the best-fit ATP function, with an Akaike weight of 100%. For STP dosimetry, the STP measurement by SPECT/CT at time point 3 (42.6 ± 1.0 h) showed a relatively low mean RD of -4.4% ± 9.4% and median RD of -0.7%. Time point 3 had the lowest RMSE value compared with those at the other 4 time points. The RMSEs of the absorbed dose RDs for time points 1-5 were 23%, 16%, 10%, 20%, and 53%, respectively. The STP dosimetry using the PBMS NLME method outperformed the Hänscheid and Madsen methods for all investigated time points. Conclusion: Our results show that a single measurement of SPECT/CT at 2 d after injection might be used to calculate accurate kidney-absorbed doses using the NLME method and PBMS.

A PBPK model for PRRT with [177Lu]Lu-DOTA-TATE: Comparison of model implementations in SAAM II and MATLAB/SimBiology

Physiologically based pharmacokinetic (PBPK) models offer the ability to simulate and predict the biodistribution of radiopharmaceuticals and have the potential to enable individualised treatment planning in molecular radiotherapy. The objective of this study was to develop and implement a whole-body compartmental PBPK model for peptide receptor radionuclide therapy (PRRT) with [177Lu]Lu-DOTA-TATE in SimBiology to allow for more complex analyses. The correctness of the model implementation was ensured by comparing its outputs, such as the time-integrated activity (TIA), with those of a PBPK model implemented in SAAM II software.

METHODS

A combined PBPK model for [68Ga]Ga-DOTA-TATE and [177Lu]Lu-DOTA-TATE was developed and implemented in both SAAM II and SimBiology. A retrospective analysis of 12 patients with metastatic neuroendocrine tumours (NETs) was conducted. First, time-activity curves (TACs) and TIAs from the two software were calculated and compared for identical parameter values. Second, pharmacokinetic parameters were fitted to activity concentrations, analysed and compared.

RESULTS

The PBPK model implemented in SimBiology produced TIA results comparable to those generated by the model implemented in SAAM II, with a relative deviation of less than 0.5% when using the same input parameters. The relative deviation of the fitted TIAs was less than 5% when model parameter values were fitted to the measured activity concentrations.

CONCLUSION

The proposed PBPK model implemented in SimBiology can be used for dosimetry in radioligand therapy and TIA prediction. Its outputs are similar to those generated by the PBPK model implemented in SAAM II, confirming the correctness of the model implementation in SimBiology.

Fit of biokinetic data in molecular radiotherapy: a machine learning approach

BACKGROUND

In literature are reported different analytical methods (AM) to choose the proper fit model and to fit data of the time-activity curve (TAC). On the other hand, Machine Learning algorithms (ML) are increasingly used for both classification and regression tasks. The aim of this work was to investigate the possibility of employing ML both to classify the most appropriate fit model and to predict the area under the curve (τ).

METHODS

Two different ML systems have been developed for classifying the fit model and to predict the biokinetic parameters. The two systems were trained and tested with synthetic TACs simulating a whole-body Fraction Injected Activity for patients affected by metastatic Differentiated Thyroid Carcinoma, administered with [131I]I-NaI. Test performances, defined as classification accuracy (CA) and percentage difference between the actual and the estimated area under the curve (Δτ), were compared with those obtained using AM varying the number of points (N) of the TACs. A comparison between AM and ML were performed using data of 20 real patients.

RESULTS

As N varies, CA remains constant for ML (about 98%), while it improves for F-test (from 62 to 92%) and AICc (from 50 to 92%), as N increases. With AM, [Formula: see text] can reach down to - 67%, while using ML [Formula: see text] ranges within ± 25%. Using real TACs, there is a good agreement between τ obtained with ML system and AM.

CONCLUSIONS

The employing of ML systems may be feasible, having both a better classification and a better estimation of biokinetic parameters.

MIB Guides: Preclinical Radiopharmaceutical Dosimetry

Preclinical dosimetry is essential for guiding the design of animal radiopharmaceutical biodistribution, imaging, and therapy experiments, evaluating efficacy and/or toxicities in such experiments, ensuring compliance with ethical standards for animal research, and, perhaps most importantly, providing reasonable initial estimates of normal-organ doses in humans, required for clinical translation of new radiopharmaceuticals. This MIB Guide provides a basic protocol for obtaining preclinical dosimetry estimates with organ-level dosimetry software.

On the Use of Artificial Intelligence for Dosimetry of Radiopharmaceutical Therapies

Routine clinical dosimetry along with radiopharmaceutical therapies is key for future treatment personalization. However, dosimetry is considered complex and time-consuming with various challenges amongst the required steps within the dosimetry workflow. The general workflow for image-based dosimetry consists of quantitative imaging, the segmentation of organs and tumors, fitting of the time-activity-curves, and the conversion to absorbed dose. This work reviews the potential and advantages of the use of artificial intelligence to improve speed and accuracy of every single step of the dosimetry workflow.

[99mTc]Tc-PentixaTec: development, extensive pre-clinical evaluation, and first human experience

PURPOSE

The clinical success non-invasive imaging of CXCR4 expression using [68 Ga]Ga-PentixaFor-PET warrants an expansion of the targeting concept towards conventional scintigraphy/SPECT with their lower cost and general availability. To this aim, we developed and comparatively evaluated a series of 99mTc-labeled cyclic pentapeptides based on the PentixaFor scaffold.

METHODS

Six mas3-conjugated CPCR4 analogs with different 4-aminobenzoic acid (Abz)-D-Ala-D-Arg-aa3 linkers (L1-L6) as well as the corresponding HYNIC- and N4-analogs of L6-CPCR4 were synthesized via standard SPPS. Competitive binding studies (IC50 and IC50inv) were carried out using Jurkat T cell lymphoma cells and [125I]FC-131 as radioligand. Internalization kinetics were investigated using hCXCR4-overexpressing Chem-1 cells. Biodistribution studies and small animal SPECT/CT imaging (1 h p.i.) were carried out using Jurkat xenograft bearing CB17/SCID mice. Based on the preclinical results, [99mTc]Tc-N4-L6-CPCR4 ([99mTc]Tc-PentixaTec) was selected for an early translation to the human setting. Five patients with hematologic malignancies underwent [99mTc]Tc-N4-L6-CPCR4 SPECT/planar imaging with individual dosimetry.

RESULTS

Of the six mas3-conjugated peptides, mas3-L6-CPCR4 (mas3-dap-r-a-Abz-CPCR4) showed the highest CXCR4 affinity (IC50 = 5.0 ± 1.3 nM). Conjugation with N4 (N4-L6-CPCR4) further improved hCXCR4 affinity to 0.6 ± 0.1 nM. [99mTc]Tc-N4-L6-CPCR4 also showed the most efficient internalization (97% of total cellular activity at 2 h) and the highest tumor accumulation (8.6 ± 1.3% iD/g, 1 h p.i.) of the compounds investigated. Therefore, [99mTc]Tc-N4-L6-CPCR4 (termed [99mTc]Tc-PentixaTec) was selected for first-in-human application. [99mTc]Tc-PentixaTec was well tolerated, exhibits a favorable biodistribution and dosimetry profile (2.1-3.4 mSv per 500 MBq) and excellent tumor/background ratios in SPECT and planar imaging.

CONCLUSION

The successive optimization of the amino acid composition of the linker structure and the N-terminal 99mTc-labeling strategies (mas3 vs HYNIC vs N4) has provided [99mTc]Tc-PentixaTec as a novel, highly promising CXCR4-targeted SPECT agent for clinical application. With its excellent CXCR4 affinity, efficient internalization, high uptake in CXCR4-expressing tissues, suitable clearance/biodistribution characteristics, and favorable human dosimetry, it holds great potential for further clinical use.

Tandem Isotope Therapy with 225Ac- and 177Lu-PSMA-617 in a Murine Model of Prostate Cancer

Radionuclide therapy targeting prostate-specific membrane antigen (PSMA) is a promising option for metastatic castration-resistant prostate cancer. Clinical experience using 177Lu or 225Ac has demonstrated encouraging treatment responses; however, responses are not durable. Dual-isotope combinations, or "tandem" approaches, may improve tolerability while retaining a high tumor dose. In this study, we directly compared α- versus β-particle treatment, as well as a combination thereof, at different stages of disease in a murine model of disseminated prostate cancer. Methods: First, to determine comparable injected activities from 177Lu- and 225Ac-PSMA-617, ex vivo biodistribution studies were performed at 5 time points after treatment of C4-2 subcutaneous tumor-bearing NSG mice. To establish a more representative model of metastatic prostate cancer, NSG mice were inoculated with luciferase-expressing C4-2 cells in the left ventricle, leading to disseminated visceral and bone lesions. At either 3 or 5 wk after inoculation, the mice were treated with equivalent tumor dose-depositing activities of 177Lu- or 225Ac-PSMA-617 alone or in combination (35 MBq of 177Lu, 40 kBq of 225Ac, or 17 MBq of 177Lu + 20 kBq 225Ac; 10/group). Disease burden was assessed by weekly bioluminescence imaging. Treatment efficacy was evaluated using whole-body tumor burden and overall survival. Results: The ex vivo biodistribution studies revealed that 35 MBq of 177Lu and 40 kBq of 225Ac yield equivalent absorbed tumor doses in a subcutaneous C4-2 model. The disease burden of mice treated at 3 wk after inoculation (microscopic disease) with 177Lu was not significantly different from that of untreated mice. However, 225Ac-PSMA-617 both as a single agent and in combination with 177Lu (17 MBq of 177Lu + 20 kBq of 225Ac) were associated with significant whole-body tumor growth retardation and survival benefit (overall survival, 8.3 wk for nontreatment, 9.4 wk for 177Lu, 15.3 wk for 225Ac alone, and 14.1 wk for tandem therapy). When treated at 5 wk after inoculation (macroscopic disease), all treatment groups showed retarded tumor growth and improved survival, with no significant differences between 225Ac alone and administration of half the 225Ac activity in tandem with 177Lu (overall survival, 7.9 wk for nontreatment, 10.3 wk for 177Lu, 14.6 wk for 225Ac alone, and 13.2 wk for tandem therapy). Conclusion: Treatment of a disseminated model of prostate cancer with simultaneous 225Ac- and 177Lu-PSMA-617 results in significantly decreased tumor growth compared with 177Lu, which was ineffective as a single agent against microscopic lesions. Mice treated later in the disease progression and bearing macroscopic, millimeter-sized lesions experienced significant tumor growth retardation and survival benefit in both monoisotopic and tandem regimens of 177Lu and 225Ac. Although the greatest benefits were observed with the single agent 225Ac, the tandem arm experienced no significant difference in disease burden or survival benefit, suggesting that the reduced activity of 225Ac was adequately compensated in the tandem arm. The superior therapeutic efficacy of 225Ac in this model suggests a preference for α-emitters alone, or possibly in combination, in the microscopic disease setting.

MIB Guides: Preclinical radiopharmaceutical dosimetry

Preclinical dosimetry is essential for guiding the design of animal radiopharmaceutical biodistribution, imaging, and therapy experiments, evaluating efficacy and/or toxicities in such experiments, ensuring compliance with ethical standards for animal research, and providing reasonable initial estimates of normal-organ doses in humans, required for clinical translation of new radiopharmaceuticals. This MIB guide provides a basic protocol for obtaining preclinical dosimetry estimates with organ-level dosimetry software.

Metabolic kinetic modeling of [11C]methionine based on total-body PET in multiple myeloma

PURPOSE

Multiple myeloma (MM) is a malignant disease characterized by the secretion of monoclonal immunoglobulins and has a high demand for amino acids. [11C]methionine total-body PET is capable of noninvasive dynamic monitoring of radiotracer in vivo, thus providing a way to reveal the dynamic changes of myeloma metabolism. This study aims to analyze the metabolic process of [11C]methionine based on kinetic modeling, and to preliminary reveal its application value in MM.

METHODS

Dynamic total-body [11C]methionine PET/CT was conducted with uEXPLORER in 12 subjects (9 MM patients and 3 controls). The tissue time activity curves (TACs) of organs and bone marrows were extracted. Model fitting of TACs was operated using PMOD Kinetic Modeling. After validation by Goodness of fit (GOF), the reversible two-tissue compartment model (2T4k) was used to further analysis. R software was used to analyze the correlation between kinetic parameters and clinical indicators.

RESULTS

The 2T4k has passed the criterion of GOF and was used to fit the data of 0-20 minutes. The [11C]methionine net uptake rate (Ki) was significantly higher in the MM lesions than in the non-myeloma controls (control: 0.040±0.007 mL/g/min, MM: 0.171±0.108 mL/g/min, p=0.009). The Ki values were found to be correlated with M protein levels in MM patients. MM patients with t(4;14) translocations had an elevated k4 value compared with t(4;14) negative patients.

CONCLUSION

MM lesions have a propensity for uptake of [11C]methionine. The serum levels of M protein are correlated with [11C]methionine uptake rate in myeloma. Metabolic classification based on the k4 value may be a promising strategy for risk stratification in MM.

Population-based model selection for an accurate estimation of time-integrated activity using non-linear mixed-effects modelling

PURPOSE

Personalized treatment planning in Molecular Radiotherapy (MRT) with accurately determining the absorbed dose is highly desirable. The absorbed dose is calculated based on the Time-Integrated Activity (TIA) and the dose conversion factor. A crucial unresolved issue in MRT dosimetry is which fit function to use for the TIA calculation. A data-driven population-based fitting function selection could help solve this problem. Therefore, this project aims to develop and evaluate a method for accurately determining TIAs in MRT, which performs a Population-Based Model Selection within the framework of the Non-Linear Mixed-Effects (NLME-PBMS) model.

METHODS

Biokinetic data of a radioligand for the Prostate-Specific Membrane Antigen (PSMA) for cancer treatment were used. Eleven fit functions were derived from various parameterisations of mono-, bi-, and tri-exponential functions. The functions' fixed and random effects parameters were fitted (in the NLME framework) to the biokinetic data of all patients. The goodness of fit was assumed acceptable based on the visual inspection of the fitted curves and the coefficients of variation of the fitted fixed effects. The Akaike weight, the probability that the model is the best among the whole set of considered models, was used to select the fit function most supported by the data from the set of functions with acceptable goodness of fit. NLME-PBMS Model Averaging (MA) was performed with all functions having acceptable goodness of fit. The Root-Mean-Square Error (RMSE) of the calculated TIAs from individual-based model selection (IBMS), a shared-parameter population-based model selection (SP-PBMS) reported in the literature, and the functions from NLME-PBMS method to the TIAs from MA were calculated and analysed. The NLME-PBMS (MA) model was used as the reference as this model considers all relevant functions with corresponding Akaike weights.

RESULTS

The function [Formula: see text] was selected as the function most supported by the data with an Akaike weight of (54 ± 11) %. Visual inspection of the fitted graphs and the RMSE values show that the NLME model selection method has a relatively better or equivalent performance than the IBMS or SP-PBMS methods. The RMSEs of the IBMS, SP-PBMS, and NLME-PBMS (f3a) methods are 7.4%, 8.8%, and 2.4%, respectively.

CONCLUSION

A procedure including fitting function selection in a population-based method was developed to determine the best fit function for calculating TIAs in MRT for a given radiopharmaceutical, organ and set of biokinetic data. The technique combines standard practice approaches in pharmacokinetics, i.e. an Akaike-weight-based model selection and the NLME model framework.

Single-time-point dosimetry using model selection and nonlinear mixed-effects modelling: a proof of concept

PURPOSE

This project aims to develop and evaluate a method for accurately determining time-integrated activities (TIAs) in single-time-point (STP) dosimetry for molecular radiotherapy. It performs a model selection (MS) within the framework of the nonlinear mixed-effects (NLME) model (MS-NLME).

METHODS

Biokinetic data of [¹¹¹In]In-DOTATATE activity in kidneys at T1 = (2.9 ± 0.6) h, T2 = (4.6 ± 0.4) h, T3 = (22.8 ± 1.6) h, T4 = (46.7 ± 1.7) h, and T5 = (70.9 ± 1.0) h post injection were obtained from eight patients using planar imaging. Eleven functions were derived from various parameterisations of mono-, bi-, and tri-exponential functions. The functions' fixed and random effects parameters were fitted simultaneously (in the NLME framework) to the biokinetic data of all patients. The Akaike weights were used to select the fit function most supported by the data. The relative deviations (RD) and the root-mean-square error (RMSE) of the calculated TIAs for the STP dosimetry at T3 = (22.8 ± 1.6) h and T4 = (46.7 ± 1.7) h p.i. were determined for all functions passing the goodness-of-fit test.

RESULTS

The function [Formula: see text] with four adjustable parameters and [Formula: see text] was selected as the function most supported by the data with an Akaike weight of (45 ± 6) %. RD and RMSE values show that the MS-NLME method performs better than functions with three or five adjustable parameters. The RMSEs of TIA_NLME-PBMS and TIA_3-parameters were 7.8% and 10.9% (for STP at T3), and 4.9% and 10.7% (for STP at T4), respectively.

CONCLUSION

An MS-NLME method was developed to determine the best fit function for calculating TIAs in STP dosimetry for a given radiopharmaceutical, organ, and patient population. The proof of concept was demonstrated for biokinetic ¹¹¹In-DOTATATE data, showing that four-parameter functions perform better than three- and five-parameter functions.

Investigation of image-based lesion and kidney dosimetry protocols for 177Lu-PSMA-I&T therapy with and without a late SPECT/CT acquisition

BACKGROUND

¹⁷⁷Lu-PSMA therapy has been successfully used to prolong the survival of patients with metastatic castration-resistant prostate cancer. Patient-specific dosimetry based on serial quantitative SPECT/CT imaging can support the understanding of dose-effect relationships. However, multiple SPECT/CT measurements can be challenging for patients, which motivates the investigation of efficient sampling schedules and their impact on dosimetry. In this study, different time samplings with respect to the number and timing of SPECT/CT acquisitions with and without a late measurement were investigated.

MATERIALS AND METHODS

In total, 43 lesions and 10 kidneys of 5 patients receiving ¹⁷⁷Lu-PSMA-I&T therapy were investigated. Whole-body SPECT/CT measurements were performed at 1, 2, 3 and 7 days post-injection. For both lesions (isocontour-based segmentation) and kidneys (CT-based segmentation), a reference model was employed including all four time points. To identify the best-matching fit function out of a pre-defined set of models, visual inspection, coefficients of variation and sum of squared errors were considered as goodness-of-fit criteria. Biologically effective doses (BEDs) calculated with different time samplings (days 1, 2, 3/1, 2, 7/1, 3, 7/2, 3, 7 and 1, 2/1, 3/1, 7) were compared to the reference.

RESULTS

The best-fit function was found to be a mono-exponential model for lesions and a bi-exponential model with a population-based parameter and two free parameters for kidneys. The BEDs calculated with the time sampling 1, 3, 7 days showed the lowest deviations from the reference for lesions with 4 ± 5%. Without day 7, still 86% of all lesions showed deviations from the reference < 10%. The outlier deviations showed a positive correlation with the effective half-life of the respective lesions. For kidneys, including days 1, 2, 3 achieved the best results with 0 ± 1%. Generally, deviations for kidneys were found to be small for all time samplings (max. 13%).

CONCLUSIONS

For combined optimization of the SPECT/CT time sampling for kidney and lesion dosimetry during ¹⁷⁷Lu-PSMA-I&T therapy, the sampling with days 1, 3, 7 showed the smallest deviation from the reference. Without a late acquisition, using the schedule with days 1, 2, 3 is likewise feasible.

Single-time-point estimation of absorbed doses in PRRT using a non-linear mixed-effects model

INTRODUCTION

Estimation of accurate time-integrated activity coefficients (TIACs) and radiation absorbed doses (ADs) is desirable for treatment planning in peptide-receptor radionuclide therapy (PRRT). This study aimed to investigate the accuracy of a simplified dosimetry using a physiologically-based pharmacokinetic (PBPK) model, a nonlinear mixed effect (NLME) model, and single-time-point imaging to calculate the TIACs and ADs of ⁹⁰Y-DOTATATE in various organs of dosimetric interest and tumors.

MATERIALS & METHODS

Biokinetic data of ¹¹¹In-DOTATATE in tumors, kidneys, liver, spleen, and whole body were obtained from eight patients using planar scintigraphic imaging at T1 = (2.9 ± 0.6), T2 = (4.6 ± 0.4), T3 = (22.8 ± 1.6), T4 = (46.7 ± 1.7) and T5 = (70.9 ± 1.0) h post injection. Serum activity concentration was measured at 5 and 15 min; 0.5, 1, 2, and 4 h; and 1, 2, and 3 d p.i.. A published PBPK model for PRRT, NLME, and a single-time-point imaging datum at different time points were used to calculate TIACs in tumors, kidneys, liver, spleen, whole body, and serum. Relative deviations (RDs) (median [min, max]) between the calculated TIACs from single-time-point imaging were compared to the TIACs calculated from the all-time-points fit. The root mean square error (RMSE) of the difference between the computed ADs from the single-time-point imaging and reference ADs from the all-time point fittings were analyzed. A joint root mean square error RMSE_joint of the ADs was calculated with the RSME from both the tumor and kidneys to sort the time points concerning accurate results for the kidneys and tumor dosimetry. The calculations of TIACs and ADs from the single-time-point dosimetry were repeated using the sum of exponentials (SOE) approach introduced in the literature. The RDs and the RSME of the PBPK approach in our study were compared to the SOE approach.

RESULTS

Using the PBPK and NLME models and the biokinetic measurements resulted in a good fit based on visual inspection of the fitted curves and the coefficient of variation CV of the fitted parameters (<50%). T4 was identified being the time point with a relatively low median and range of TIACs RDs, i.e., 5 [1, 21]% and 2 [-15, 21]% for kidneys and tumors, respectively. T4 was found to be the time point with the lowest joint root mean square error RMSE_joint of the ADs. Based on the RD and RMSE, our results show a similar performance as the SOE and NLME model approach.

SUMMARY

In this study, we introduced a simplified calculation of TIACs/ADs using a PBPK model, an NLME model, and a single-time-point measurement. Our results suggest a single measurement might be used to calculate TIACs/ADs in the kidneys and tumors during PRRT.

SPECT/CT-based dosimetry of salivary glands and iodine-avid lesions following 131I therapy

Objective

The purpose was to provide uptake and radiation dose estimates to salivary glands (SG) and pathologic lesions following radioiodine therapy (RIT) of differentiated thyroid cancer patients (DTC).

Methods

A group of DTC patients (n = 25) undergoing ¹³¹I therapy joined this study with varying amounts of therapeutic activity. Sequential SPECT/CT scans were acquired at 4 ± 2, 24 ± 2, and 168 ± 3 h following administration of 3497-9250 MBq 131I. An earlier experiment with Acrylic glass body phantom (PET Phantom NEMA 2012 / IEC 2008) was conducted for system calibration including scatter, partial volume effect and count loss correction. Dose calculation was made via IDAC-Dose 2.1 code.

Results

The absorbed dose to parotid glands was 0.04-0.97 Gy/GBq (median: 0.26 Gy/GBq). The median absorbed dose to submandibular glands was 0.14 Gy/GBq (0.05 to 0.56 Gy/GBq). The absorbed dose to thyroid residues was from 0.55 to 399.5 Gy/GBq (median: 21.8 Gy/GBq), and that to distal lesions ranged from 0.78 to 28.0 Gy/GBq (median: 3.12 Gy/GBq). 41% of the thyroid residues received dose > 80 Gy, 18% between 70-80 Gy, 18% between 40-70 Gy, and 23% has dose < 40 Gy. In contrast, 18% of the metastases exhibited a dose > 80 Gy, 9% between 40-60 Gy, and the dose to the vast majority of lesions (64%) was < 40 Gy.

Conclusion

It was inferred that dose estimation after RIT with SPECT/CT is feasible to apply, together with good agreement with published ¹²⁴I PET/CT dose estimates. A broad and sub-effective dose range was estimated for thyroid residues and distal lesions. Moreover, the current methodology might be useful for establishing a dose-effect relationship and radiation-induced salivary glands damage after RIT.

Personalized dosimetry for a deeper understanding of metastatic response to high activity 131I-mIBG therapy in high risk relapsed refractory neuroblastoma

Background

Dosimetry in molecular radiotherapy for personalized treatment is assuming a central role in clinical management of aggressive/relapsed tumors. Relapsed/refractory metastatic high-risk neuroblastoma (rrmHR-NBL) has a poor prognosis and high-activity ¹³¹I-mIBG therapy could represent a promising strategy. The primary aim of this case series study was to report the absorbed doses to whole-body (D_WB ), red marrow (D_RM ) and lesions (D_Lesion ). A secondary aim was to correlate D_Lesion values to clinical outcome.

Methods

Fourteen patients affected by rrmHR-NBL were treated with high-activity ¹³¹I-mIBG therapy (two administrations separated by 15 days). The first administration was weight-based whereas the second one was dosimetry-based (achieving D_WB equals to 4 Gy). In all patients D_WB and D_RM was assessed; 9/14 patients were selected for D_Lesion evaluation using planar dosimetric approach (13 lesions evaluated). Treatment response was classified as progressive and stable disease (PD and SD), partial and complete response (PR and CR) according to the International Neuroblastoma Response Criteria. Patients were divided into two groups: Responder (CR, PR, SD) and Non-Responder (PD), correlating treatment response to D_Lesion value.

Results

The cumulative D_WB , D_RM and D_Lesion ranged from (1.5; 4.5), (1.0; 2.6) and (44.2; 585.8) Gy. A linear correlation between D_WB and D_RM and a power law correlation between the absorbed dose to WB normalized for administered activity and the mass of the patient were observed. After treatment 3, 2, 4 and 5 patients showed CR, PR, SD and PD respectively, showing a correlation between D_Lesion and the two response group.

Conclusions

Our experience demonstrated feasibility of high activity therapy of ¹³¹I-mIBG in rrmHR-NBL children as two administration intensive strategy. Dosimetric approach allowed a tailored high dose treatment maximizing the benefits of radionuclide therapy for pediatric patients while maintaining a safety profile. The assesment of D_Lesion contributed to have a deeper understaning of metabolic treatment effects.

A population-based method to determine the time-integrated activity in molecular radiotherapy

BACKGROUND

The calculation of time-integrated activities (TIAs) for tumours and organs is required for dosimetry in molecular radiotherapy. The accuracy of the calculated TIAs is highly dependent on the chosen fit function. Selection of an adequate function is therefore of high importance. However, model (i.e. function) selection works more accurately when more biokinetic data are available than are usually obtained in a single patient. In this retrospective analysis, we therefore developed a method for population-based model selection that can be used for the determination of individual time-integrated activities (TIAs). The method is demonstrated at an example of [¹⁷⁷Lu]Lu-PSMA-I&T kidneys biokinetics. It is based on population fitting and is specifically advantageous for cases with a low number of available biokinetic data per patient.

METHODS

Renal biokinetics of [¹⁷⁷Lu]Lu-PSMA-I&T from thirteen patients with metastatic castration-resistant prostate cancer acquired by planar imaging were used. Twenty exponential functions were derived from various parameterizations of mono- and bi-exponential functions. The parameters of the functions were fitted (with different combinations of shared and individual parameters) to the biokinetic data of all patients. The goodness of fits were assumed as acceptable based on visual inspection of the fitted curves and coefficients of variation CVs < 50%. The Akaike weight (based on the corrected Akaike Information Criterion) was used to select the fit function most supported by the data from the set of functions with acceptable goodness of fit.

RESULTS

The function [Formula: see text] with shared parameter [Formula: see text] was selected as the function most supported by the data with an Akaike weight of 97%. Parameters [Formula: see text] and [Formula: see text] were fitted individually for every patient while parameter [Formula: see text] was fitted as a shared parameter in the population yielding a value of 0.9632 ± 0.0037.

CONCLUSIONS

The presented population-based model selection allows for a higher number of parameters of investigated fit functions which leads to better fits. It also reduces the uncertainty of the obtained Akaike weights and the selected best fit function based on them. The use of the population-determined shared parameter for future patients allows the fitting of more appropriate functions also for patients for whom only a low number of individual data are available.

Quantitative SPECT/CT for dosimetry of peptide receptor radionuclide therapy

Neuroendocrine tumors (NETs) are uncommon malignancies of increasing incidence and prevalence. As these slow growing tumors usually overexpress somatostatin receptors (SSTRs), the use of ⁶⁸Ga-DOTA-peptides (gallium-68 chelated with dodecane tetra-acetic acid to somatostatin), which bind to the SSTRs, allows for PET based imaging and selection of patients for peptide receptor radionuclide therapy (PRRT). PRRT with radiolabeled somatostatin analogues such as ¹⁷⁷Lu-DOTATATE (lutetium-177-[DOTA,Tyr³]-octreotate), is mainly used for the treatment of metastatic or inoperable NETs. However, PRRT is generally administered at a fixed injected activity in order not to exceed dose limits in critical organs, which is suboptimal given the variability in radiopharmaceutical uptake among patients. Advances in SPECT (single photon emission computed tomography) imaging enable the absolute quantitative measure of the true radiopharmaceutical distribution providing for PRRT dosimetry in each patient. Personalized PRRT based on patient-specific dosimetry could improve therapeutic efficacy by optimizing effective tumor absorbed dose while limiting treatment related radiotoxicity.

A Physiologically Based Pharmacokinetic Model for In Vivo Alpha Particle Generators Targeting Neuroendocrine Tumors in Mice

In vivo alpha particle generators have great potential for the treatment of neuroendocrine tumors in alpha-emitter-based peptide receptor radionuclide therapy (α-PRRT). Quantitative pharmacokinetic analyses of the in vivo alpha particle generator and its radioactive decay products are required to address concerns about the efficacy and safety of α-PRRT. A murine whole-body physiologically based pharmacokinetic (PBPK) model was developed for ²¹²Pb-labeled somatostatin analogs (²¹²Pb-SSTA). The model describes pharmacokinetics of ²¹²Pb-SSTA and its decay products, including specific and non-specific glomerular and tubular uptake. Absorbed dose coefficients (ADC) were calculated for bound and unbound radiolabeled SSTA and its decay products. Kidneys received the highest ADC (134 Gy/MBq) among non-target tissues. The alpha-emitting ²¹²Po contributes more than 50% to absorbed doses in most tissues. Using this model, it is demonstrated that α-PRRT based on ²¹²Pb-SSTA results in lower absorbed doses in non-target tissue than α-PRRT based on ²¹²Bi-SSTA for a given kidneys absorbed dose. In both approaches, the energies released in the glomeruli and proximal tubules account for 54% and 46%, respectively, of the total energy absorbed in kidneys. The ²¹²Pb-SSTA-PBPK model accelerates the translation from bench to bedside by enabling better experimental design and by improving the understanding of the underlying mechanisms.

Role of Artificial Intelligence in Theranostics:: Toward Routine Personalized Radiopharmaceutical Therapies

We highlight emerging uses of artificial intelligence (AI) in the field of theranostics, focusing on its significant potential to enable routine and reliable personalization of radiopharmaceutical therapies (RPTs). Personalized RPTs require patient-specific dosimetry calculations accompanying therapy. Additionally we discuss the potential to exploit biological information from diagnostic and therapeutic molecular images to derive biomarkers for absorbed dose and outcome prediction; toward personalization of therapies. We try to motivate the nuclear medicine community to expand and align efforts into making routine and reliable personalization of RPTs a reality.

A Whole-Body Physiologically Based Pharmacokinetic Model for Alpha Particle Emitting Bismuth in Rats

Introduction/Aim: α particle emitting bismuth (²¹²Bi) as decay product of ²¹²Pb-labeled pharmaceuticals has been effective in targeted α particle therapy (TAT). Estimating the contribution of ²¹²Bi released from its chelator to the absorbed doses in nontarget tissues is challenging in TAT. Physiologically based pharmacokinetic (PBPK) modeling can help overcome this limitation. Therefore, a whole-body ²¹²Bi-PBPK model was developed to describe the pharmacokinetics (PKs) of ²¹²Bi in rats. Materials and Methods: The rat ²¹²Bi-PBPK model was implemented using the modeling software SAAM II with data and parameter values from the literature. Besides other mechanisms, ²¹²Bi interactions with red blood cells, high molecular weight plasma protein, and intracellular biological thiols are described. Important PK parameters were fitted to time-activity data. Absorbed dose coefficients (ADCs) were calculated for injecting 0.774 fmol of ²¹²Bi. Results: ²¹²Bi uptake rates of liver, bone, small intestine, bone marrow, skin, and muscle were (0.86 ± 0.13), (3.85 ± 0.63), (0.27 ± 0.05), (1.44 ± 0.29), (0.04 ± 0.01), and (0.007 ± 0.007) per min with corresponding ADCs of 0.09, 0.03, 0.03, 0.07, 0.01, and 0.003 mGy/kBq, respectively. An ADC of 0.70 mGy/kBq was determined for kidneys. Conclusion: Kidneys are the dose-limiting organs in ²¹²Bi-based TAT. The ²¹²Bi-PBPK model is an effective tool to investigate the ²¹²Bi biodistribution in murine models. Integrating the ²¹²Bi-PBPK model into other murine and human PBPK models of α particle generators can help study the efficacy and safety of TAT.

Important pharmacokinetic parameters for individualization of 177 Lu-PSMA therapy: A global sensitivity analysis for a physiologically-based pharmacokinetic model

PURPOSE

The knowledge of the contribution of anatomical and physiological parameters to interindividual pharmacokinetic differences could potentially be used to improve individualized treatment planning for radionuclide therapy. The aim of this study was therefore to identify the physiologically based pharmacokinetic (PBPK) model parameters that determine the interindividual variability of absorbed doses (ADs) to kidneys and tumor lesions in therapy with ¹⁷⁷ Lu-labeled PSMA-targeting radioligands.

METHODS

A global sensitivity analysis (GSA) with the extended Fourier Amplitude Sensitivity Test (eFAST) algorithm was performed. The whole-body PBPK model for PSMA-targeting radioligand therapy from our previous studies was used in this study. The model parameters of interest (input of the GSA) were the organ receptor densities [R₀ ], the organ blood flows f, and the organ release rates λ. These parameters were systematically sampled NE times according to their distribution in the patient population. The corresponding pharmacokinetics were simulated and the ADs (model output) to kidneys and tumor lesions were collected. The main effect S i and total effect S Ti were calculated using the eFAST algorithm based on the variability of the model output: The main effect S i of input parameter i represents the reduction in variance of the output if the "true" value of parameter i would be known. The total effect S Ti of an input parameter i represents the proportion of variance remaining if the "true" values of all other input parameters except for i are known. The numbers of samples NE were increased up to 8193 to check the stability (i.e., convergence) of the calculated main effects S i and total effects S Ti .

RESULTS

From the simulations, the relative interindividual variability of ADs in the kidneys (coefficient of variation CV = 31%) was lower than that of ADs in the tumors (CV up to 59%). Based on the GSA, the most important parameters that determine the ADs to the kidneys were kidneys flow ( S i  = 0.36, S Ti  = 0.43) and kidneys receptor density ( S i  = 0.25, S Ti  = 0.30). Tumor receptor density was identified as the most important parameter determining the ADs to tumors ( S i and S Ti up to 0.72).

CONCLUSIONS

The results suggest that an accurate measurement of receptor density and flow before therapy could be a promising approach for developing an individualized treatment with ¹⁷⁷ Lu-labeled PSMA-targeting radioligands.

Comparison of commercial dosimetric software platforms in patients treated with 177 Lu-DOTATATE for peptide receptor radionuclide therapy

Purpose

The aim of this study was to quantitatively compare five commercial dosimetric software platforms based on the analysis of clinical datasets of patients who benefited from peptide receptor radionuclide therapy (PRRT) with ¹⁷⁷Lu‐DOTATATE (LUTATHERA^®).

Methods

The dosimetric analysis was performed on two patients during two cycles of PRRT with ¹⁷⁷Lu. Single photon emission computed tomography/computed tomography images were acquired at 4, 24, 72, and 192 h post injection. Reconstructed images were generated using Dosimetry Toolkit^® (DTK) from Xeleris™ and HybridRecon‐Oncology version_1.3_Dicom (HROD) from HERMES. Reconstructed images using DTK were analyzed using the same software to calculate time‐integrated activity coefficients (TIAC), and mean absorbed doses were estimated using OLINDA/EXM V1.0 with mass correction. Reconstructed images from HROD were uploaded into PLANET® OncoDose from DOSIsoft, STRATOS from Phillips, Hybrid Dosimetry Module™ from HERMES, and SurePlan™ MRT from MIM. Organ masses, TIACs, and mean absorbed doses were calculated from each application using their recommendations.

Results

The majority of organ mass estimates varied by <9.5% between all platforms. The highest variability for TIAC results between platforms was seen for the kidneys (28.2%) for the two patients and the two treatment cycles. Relative standard deviations in mean absorbed doses were slightly higher compared with those observed for TIAC, but remained of the same order of magnitude between all platforms.

Conclusions

When applying a similar processing approach, results obtained were of the same order of magnitude regardless of the platforms used. However, the comparison of the performances of currently available platforms is still difficult as they do not all address the same parts of the dosimetric analysis workflow. In addition, the way in which data are handled in each part of the chain from data acquisition to absorbed doses may be different, which complicates the comparison exercise. Therefore, the dissemination of commercial solutions for absorbed dose calculation calls for the development of tools and standards allowing for the comparison of the performances between dosimetric software platforms.

Influence of sampling schedules on [177Lu]Lu-PSMA dosimetry

BACKGROUND

Individualized dosimetry is recommended for [¹⁷⁷Lu]Lu-PSMA radioligand therapy (RLT) which is resource-intensive and protocols are often not optimized. Therefore, a simulation study was performed focusing on the determination of efficient optimal sampling schedules (OSS) for renal and tumour dosimetry by investigating different numbers of time points (TPs).

METHODS

Sampling schedules with 1-4 TPs were investigated. Time-activity curves of the kidneys and two tumour lesions were generated based on a physiologically based pharmacokinetic (PBPK) model and biokinetic data of 13 patients who have undergone [¹⁷⁷Lu]Lu-PSMA I&T therapy. Systematic and stochastic noise of different ratios was considered when modelling time-activity data sets. Time-integrated activity coefficients (TIACs) were estimated by simulating the hybrid planar/SPECT method for schedules comprising at least two TPs. TIACs based on one single SPECT/CT measurement were estimated using an approximation for reducing the number of fitted parameters. For each sampling schedule, the root-mean-squared error (RMSE) of the deviations of the simulated TIACs from the ground truths for 1000 replications was used as a measure for accuracy and precision.

RESULTS

All determined OSS included a late measurement at 192 h p.i., which was necessary for accurate and precise tumour TIACs. OSS with three TPs were identified to be 3-4, 96-100 and 192 h with an additional SPECT/CT measurement at the penultimate TP. Kidney and tumour RMSE of 6.4 to 7.7% and 6.3 to 7.8% were obtained, respectively. Shortening the total time for dosimetry to e.g. 96 h resulted in kidney and tumour RMSE of 6.8 to 8.3% and 9.1 to 11%, respectively. OSS with four TPs showed similar results as with three TPs. Planar images at 4 and 68 h and a SPECT/CT shortly after the 68 h measurement led to kidney and tumour RMSE of 8.4 to 12% and 12 to 16%, respectively. One single SPECT/CT measurement at 52 h yielded good approximations for the kidney TIACs (RMSE of 7.0%), but led to biased tumour TIACs.

CONCLUSION

OSS allow improvements in accuracy and precision of renal and tumour dosimetry for [¹⁷⁷Lu]Lu-PSMA therapy with potentially less effort. A late TP is important regarding accurate tumour TIACs.

Technical Note: Rapid multiexponential curve fitting algorithm for voxel-based targeted radionuclide dosimetry

BACKGROUND

Dosimetry in nuclear medicine often relies on estimating pharmacokinetics based on sparse temporal data. As analysis methods move toward image-based three-dimensional computation, it becomes important to interpolate and extrapolate these data without requiring manual intervention; that is, in a manner that is highly efficient and reproducible. Iterative least-squares solvers are poorly suited to this task because of the computational overhead and potential to optimize to local minima without applying tight constraints at the outset.

METHODOLOGY

This work describes a fully analytical method for solving three-phase exponential time-activity curves based on three measured time points in a manner that may be readily employed by image-based dosimetry tools. The methodology uses a series of conditional statements and a piecewise approach for solving exponential slope directly through measured values in most instances. The proposed algorithm is tested against a purpose-designed iterative fitting technique and linear piecewise method followed by single exponential in a cohort of ten patients receiving ¹⁷⁷ Lu-DOTA-Octreotate therapy.

RESULTS

Tri-exponential time-integrated values are shown to be comparable to previously published methods with an average difference between organs when computed at the voxel level of 9.8 ± 14.2% and -3.6 ± 10.4% compared to iterative and interpolated methods, respectively. Of the three methods, the proposed tri-exponential algorithm was most consistent when regional time-integrated activity was evaluated at both voxel- and whole-organ levels. For whole-body SPECT imaging, it is possible to compute 3D time-integrated activity maps in <5 min processing time. Furthermore, the technique is able to predictably and reproducibly handle artefactual measurements due to noise or spatial misalignment over multiple image times.

CONCLUSIONS

An efficient, analytical algorithm for solving multiphase exponential pharmacokinetics is reported. The method may be readily incorporated into voxel-dose routines by combining with widely available image registration and radiation transport tools.

ICRP Publication 140: Radiological Protection in Therapy with Radiopharmaceuticals

Radiopharmaceuticals are increasingly used for the treatment of various cancers with novel radionuclides, compounds, tracer molecules, and administration techniques. The goal of radiation therapy, including therapy with radiopharmaceuticals, is to optimise the relationship between tumour control probability and potential complications in normal organs and tissues. Essential to this optimisation is the ability to quantify the radiation doses delivered to both tumours and normal tissues. This publication provides an overview of therapeutic procedures and a framework for calculating radiation doses for various treatment approaches. In radiopharmaceutical therapy, the absorbed dose to an organ or tissue is governed by radiopharmaceutical uptake, retention in and clearance from the various organs and tissues of the body, together with radionuclide physical half-life. Biokinetic parameters are determined by direct measurements made using techniques that vary in complexity. For treatment planning, absorbed dose calculations are usually performed prior to therapy using a trace-labelled diagnostic administration, or retrospective dosimetry may be performed on the basis of the activity already administered following each therapeutic administration. Uncertainty analyses provide additional information about sources of bias and random variation and their magnitudes; these analyses show the reliability and quality of absorbed dose calculations. Effective dose can provide an approximate measure of lifetime risk of detriment attributable to the stochastic effects of radiation exposure, principally cancer, but effective dose does not predict future cancer incidence for an individual and does not apply to short-term deterministic effects associated with radiopharmaceutical therapy. Accident prevention in radiation therapy should be an integral part of the design of facilities, equipment, and administration procedures. Minimisation of staff exposures includes consideration of equipment design, proper shielding and handling of sources, and personal protective equipment and tools, as well as education and training to promote awareness and engagement in radiological protection. The decision to hold or release a patient after radiopharmaceutical therapy should account for potential radiation dose to members of the public and carers that may result from residual radioactivity in the patient. In these situations, specific radiological protection guidance should be provided to patients and carers.

A simulation-based method to determine optimal sampling schedules for dosimetry in radioligand therapy

AIM

For dosimetry in radioligand therapy, the time-integrated activity coefficients (TIACs) for organs at risk and for tumour lesions have to be determined. The used sampling scheme affects the TIACs and therefore the calculated absorbed doses. The aim of this work was to develop a general and flexible method, which analyses numerous clinically applicable sampling schedules using true time-activity curves (TACs) of virtual patients.

METHOD

Nine virtual patients with true TACs of the tumours were created using a physiologically-based pharmacokinetic (PBPK) model and individual biokinetic data of five patients with neuroendocrine tumours and four with meningioma. ¹¹¹In-DOTATATE was used for pre-therapeutic dosimetry. In total, 15,120 sampling schemes, each consisting of 4 time points, were investigated. Gaussian noise of different levels was added to the corresponding true time-activity points. A bi-exponential function was used to fit the simulated time-activity data. For each investigated sampling schedule, 1000 replications were performed. Patient-specific and population-specific optimal sampling schedules were determined using the relative root-mean-square error (rRMSE). Furthermore, the fractions of TIACs a˜ deviating >5% (f_Δa˜>5%) and >10% (f_Δa˜>10%) from the true TIAC a˜_true were used for additional evaluations e.g. to investigate the effect of varying single time points.

RESULTS

Almost all patient-specific and all population-specific optimal sampling schedules have t₄≥96h for all noise levels. Changing the latest time point from the population-specific optimal time to e.g. 48h leads to a median increase of f_Δa˜>10% from 0.1% to 88% for the lowest investigated noise level. Using the determined population-specific optimal schedules, results in more accurate and precise results than established schedules from the literature.

CONCLUSION

A method of determining the optimal sampling schedule for dosimetry, which considers clinical working hours and measurement uncertainties, has been developed and applied. The simulation study shows that optimised sampling schedules result in high accuracy and precision of the determined TIACs.

In Vivo Biokinetics of 177Lu-OPS201 in Mice and Pigs as a Model for Predicting Human Dosimetry

Introduction

¹⁷⁷Lu-OPS201 is a high-affinity somatostatin receptor subtype 2 antagonist for PRRT in patients with neuroendocrine tumors. The aim is to find the optimal scaling for dosimetry and to compare the biokinetics of ¹⁷⁷Lu-OPS201 in animals and humans.

Methods

Data on biokinetics of ¹⁷⁷Lu-OPS201 were analyzed in athymic nude Foxn1^nu mice (28 F, weight: 26 ± 1 g), Danish Landrace pigs (3 F-1 M, weight: 28 ± 2 kg), and patients (3 F-1 M, weight: 61 ± 17 kg) with administered activities of 0.19–0.27 MBq (mice), 97–113 MBq (pigs), and 850–1086 MBq (patients). After euthanizing mice (up to 168 h), the organ-specific activity contents (including blood) were measured. Multiple planar and SPECT/CT scans were performed until 250 h (pigs) and 72 h (patients) to quantify the uptake in the kidneys and liver. Blood samples were taken up to 23 h (patients) and 300 h (pigs). In pigs and patients, kidney protection was applied. Time-dependent uptake data sets were created for each species and organ/tissue. Biexponential fits were applied to compare the biokinetics in the kidneys, liver, and blood of each species. The time-integrated activity coefficients (TIACs) were calculated by using NUKFIT. To determine the optimal scaling, several methods (relative mass scaling, time scaling, combined mass and time scaling, and allometric scaling) were compared.

Results

A fast blood clearance of the compound was observed in the first phase (<56 h) for all species. In comparison with patients, pigs showed higher liver retention. Based on the direct comparison of the TIACs, an underestimation in mice (liver and kidneys) and an overestimation in pigs' kidneys compared to the patient data (kidney TIAC: mice = 1.4 h, pigs = 7.7 h, and patients = 5.8 h; liver TIAC: mice = 0.7 h, pigs = 4.1 h, and patients = 5.3 h) were observed. Most similar TIACs were obtained by applying time scaling (mice) and combined scaling (pigs) (kidney TIAC: mice = 3.9 h, pigs = 4.8 h, and patients = 5.8 h; liver TIAC: mice = 0.9 h, pigs = 4.7 h, and patients = 5.3 h).

Conclusion

If the organ mass ratios between the species are high, the combined mass and time scaling method is optimal to minimize the interspecies differences. The analysis of the fit functions and the TIACs shows that pigs are better mimicking human biokinetics.

EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations

A framework is proposed for modelling the uncertainty in the measurement processes constituting the dosimetry chain that are involved in internal absorbed dose calculations. The starting point is the basic model for absorbed dose in a site of interest as the product of the cumulated activity and a dose factor. In turn, the cumulated activity is given by the area under a time-activity curve derived from a time sequence of activity values. Each activity value is obtained in terms of a count rate, a calibration factor and a recovery coefficient (a correction for partial volume effects). The method to determine the recovery coefficient and the dose factor, both of which are dependent on the size of the volume of interest (VOI), are described. Consideration is given to propagating estimates of the quantities concerned and their associated uncertainties through the dosimetry chain to obtain an estimate of mean absorbed dose in the VOI and its associated uncertainty. This approach is demonstrated in a clinical example.

ESTIMATION OF THE ORGAN ABSORBED DOSES AND EFFECTIVE DOSE FROM 68Ga-PSMA-11 PET SCAN

68Ga-PSMA-11 PET/CT has been proven to have high clinical value for imaging of prostate cancer and rapidly gained popularity. In this study, we aimed to investigate absorbed doses of 68Ga-PSMA-11. Seven patients (mean age = 66.9 ± 6.6 years, range: 57-79 years) were enrolled in the study. Whole body PET images were acquired with multiple time points. MIRD method, NUKFIT and OLINDA/EXM software were used for dosimetry calculations. Kidneys, bladder wall, salivary and lacrimal glands received the highest absorbed dose. Estimated absorbed doses to these organs after injection of 150 MBq 68Ga-PSMA-11 were 37.0, 12.6, 14.4 and 6.3 mSv, respectively. Effective dose from PET scanning with 150 MBq injected 68Ga-PSMA-11 was 2.5 mSv. In conclusion, 68Ga-PSMA-11 has a favorable dosimetry profile similar to the 68Ga labeled octreotide analogs, which are used safely in routine clinical practices for many years. No adverse effects were reported. The kidneys were the dose-limiting organs.

Activity estimation and biokinetic analysis of 99mTc-DMSA in renal infant patients using a gamma camera

Biokinetic data from the administration of radiopharmaceuticals is essential in nuclear medicine dosimetry. It has particular significance in children, as their metabolism is very different from adults. Biokinetic models for paediatric patients could therefore need to be adapted to better reflect their absorption, retention and excretion functions, when compared to adults. Obtaining quality in vivo infant or paediatric biokinetic data is then essential to improve the available reference models, which in turn can lead to the optimization of paediatric procedures and protocols in clinical practice. This study analyses the biokinetic behaviour of ^99mTc-dimercaptosuccinic acid (DMSA), in 8 infants aged 4 months to 2 years old, through an imaging study using a gamma camera, and compares the obtained values with those obtained with the reference ICRP biokinetic model. The in vivo data was treated using an adapted methodology from the MIRD 16 pamphlet. Activity curves for the liver, the kidney and the whole body, were built, and new effective absorption, retention and excretion half-lives were estimated, and compared with the reference biokinetic parameters of ICRP 128. The obtained residence time in the kidneys of 2.56 h, has a deviation of 30.8% to the ICRP 128 value of 3.70 h. The obtained maximum uptake in the kidneys was of 0.22/A₀, which compares to the value of 0.31/A₀ for ICRP. The obtained biokinetic parameters were used to estimate the absorbed dose. The obtained dose values are smaller than the reference ICRP 128 ones by 32.1% in the kidneys, and 18.4% in the liver.

Evaluation of the dosimetry approaches in ablation treatment of thyroid cancer

In this study, we aimed to evaluate dosimetric approaches in ablation treatment of Differentiated Thyroid Carcinoma (DTC) without interrupting the clinical routine. Prior to therapy, 10.7 MBq ¹³¹I in average was orally given to 24 patients suffering from DTC. MIRD formalism was used for dosimetric calculations. For blood and bone marrow dosimetry, blood samples and whole‐body counts were collected at 2, 24, 72, and 120 h after I‐131 administration. For remnant tissue dosimetry, uptake measurements were performed at the same time intervals. To estimate the remnant volume, anterior and lateral planar gamma camera images were acquired with a reference source within the field of view at 24 h after I‐131 administration. Ultrasound imaging was also performed. Treatment activities determined with the fixed activity method were administered to the patients. Secondary cancer risk relative to applied therapy was evaluated for dosimetric approaches. The average dose to blood and bone marrow were determined as 0.15 ± 0.04 and 0.11 ± 0.04 Gy/GBq, respectively. The average remnant tissue dose was 0.58 ± 0.52 Gy/MBq and the corresponding required activity to ablate the remnant was approximately 1.3 GBq of ¹³¹I. A strong correlation between 24th‐hour uptake and time‐integrated activity coefficient values was obtained. Compared to fixed activity method, approximately five times higher secondary cancer risk was determined in bone marrow dosimetry, while the risk was about three times lower in lesion‐based dosimetry.

Microscale radiosynthesis, preclinical imaging and dosimetry study of [18F]AMBF3-TATE: A potential PET tracer for clinical imaging of somatostatin receptors

BACKGROUND

Peptides labeled with positron-emitting isotopes are emerging as a versatile class of compounds for the development of highly specific, targeted imaging agents for diagnostic imaging via positron-emission tomography (PET) and for precision medicine via theranostic applications. Despite the success of peptides labeled with gallium-68 (for imaging) or lutetium-177 (for therapy) in the clinical management of patients with neuroendocrine tumors or prostate cancer, there are significant advantages of using fluorine-18 for imaging. Recent developments have greatly simplified such labeling: in particular, labeling of organotrifluoroborates via isotopic exchange can readily be performed in a single-step under aqueous conditions and without the need for HPLC purification. Though an automated synthesis has not yet been explored, microfluidic approaches have emerged for ¹⁸F-labeling with high speed, minimal reagents, and high molar activity compared to conventional approaches. As a proof-of-concept, we performed microfluidic labeling of an octreotate analog ([¹⁸F]AMBF₃-TATE), a promising ¹⁸F-labeled analog that could compete with [⁶⁸Ga]Ga-DOTATATE with the advantage of providing a greater number of patient doses per batch produced.

METHODS

Both [¹⁸F]AMBF₃-TATE and [⁶⁸Ga]Ga-DOTATATE were labeled, the former by microscale methods adapted from manual labeling, and were imaged in mice bearing human SSTR2-overexpressing, rat SSTR2 wildtype, and SSTR2-negative xenografts. Furthermore, a dosimetry analysis was performed for [¹⁸F]AMBF₃-TATE.

RESULTS

The micro-synthesis exhibited highly-repeatable performance with radiochemical conversion of 50 ± 6% (n = 15), overall decay-corrected radiochemical yield of 16 ± 1% (n = 5) in ~40 min, radiochemical purity >99%, and high molar activity. Preclinical imaging with [¹⁸F]AMBF₃-TATE in SSTR2 tumor models correlated well with [⁶⁸Ga]Ga-DOTATATE. The favorable biodistribution, with the highest tracer accumulation in the bladder followed distantly by gastrointestinal tissues, resulted in 1.26 × 10^-2 mSv/MBq maximal estimated effective dose in human, a value lower than that reported for current clinical ¹⁸F- and ⁶⁸Ga-labeled compounds.

CONCLUSIONS

The combination of novel chemical approaches to ¹⁸F-labeling and microdroplet radiochemistry have the potential to serve as a platform for greatly simplified development and production of ¹⁸F-labeled peptide tracers. Favorable preclinical imaging and dosimetry of [¹⁸F]AMBF₃-TATE, combined with a convenient synthesis, validate this assertion and suggest strong potential for clinical translation.

The Effect of Total Tumor Volume on the Biologically Effective Dose to Tumor and Kidneys for 177Lu-Labeled PSMA Peptides

The aim of this work was to simulate the effect of prostate-specific membrane antigen (PSMA)-positive total tumor volume (TTV) on the biologically effective doses (BEDs) to tumors and organs at risk in patients with metastatic castration-resistant prostate cancer who are undergoing ¹⁷⁷Lu-PSMA radioligand therapy. Methods: A physiologically based pharmacokinetic model was fitted to the data of 13 patients treated with ¹⁷⁷Lu-PSMA I&T (a PSMA inhibitor for imaging and therapy). The tumor, kidney, and salivary gland BEDs were simulated for TTVs of 0.1-10 L. The activity and peptide amounts leading to an optimal tumor-to-kidneys BED ratio were also investigated. Results: When the TTV was increased from 0.3 to 3 L, the simulated BEDs to tumors, kidneys, parotid glands, and submandibular glands decreased from 22 ± 15 to 11.0 ± 6.0 Gy_1.49, 6.5 ± 2.3 to 3.7 ± 1.4 Gy_2.5, 11.0 ± 2.7 to 6.4 ± 1.9 Gy_4.5, and 10.9 ± 2.7 to 6.3 ± 1.9 Gy_4.5, respectively (where the subscripts denote that an α/β of 1.49, 2.5, or 4.5 Gy was used to calculate the BED). The BED to the red marrow increased from 0.17 ± 0.05 to 0.32 ± 0.11 Gy₁₅ For patients with a TTV of more than 0.3 L, the optimal amount of peptide was 273 ± 136 nmol and the optimal activity was 10.4 ± 4.4 GBq. Conclusion: This simulation study suggests that in patients with large PSMA-positive tumor volumes, higher activities and peptide amounts can be safely administered to maximize tumor BEDs without exceeding the tolerable BED to the organs at risk.

Treatment planning in PRRT based on simulated PET data and a PBPK model

Aim: To investigate the accuracy of treatment planning in peptide-receptor radionuclide therapy (PRRT) based on simulated PET data (using a PET noise model) and a physiologically based pharmacokinetic (PBPK) model. Methods: The parameters of a PBPK model were fitted to the biokinetic data of 15 patients. True mathematical phantoms of patients (MPPs) were the PBPK model with the fitted parameters. PET measurements after bolus injection of 150 MBq 68Ga-DOTATATE were simulated for the true MPPs. PET noise with typical noise levels was added to the data (i.e. c=0.3 [low], 3, 30 and 300 [high]). Organ activity data in the kidneys, tumour, liver and spleen were simulated at 0.5, 1 and 4 h p.i. PBPK model parameters were fitted to the simulated noisy PET data to derive the PET-predicted MPPs. Therapy was simulated assuming an infusion of 3.3 GBq of 90Y-DOTATATE over 30 min. Time-integrated activity coefficients (TIACs) of simulated therapy in tumour, kidneys, liver, spleen and remainder were calculated from both, true MPPs (true TIACs) and predicted MPPs (predicted TIACs). Variability v between true TIACs and predicted TIACs were calculated and analysed. Variability< 10 % was considered to be an accurate prediction. Results: For all noise level, variabilities for the kidneys, liver, and spleen showed an accurate prediction for TIACs, e.g. c=300: vkidney=(5 ± 2)%, vliver=(5 ± 2)%, vspleen=(4 ± 2)%. However, tumour TIAC predictions were not accurate for all noise levels, e.g. c=0.3: vtumour=(8 ± 5)%. Conclusion: PET based treatment planning with kidneys as the dose limiting organ seems possible for all reported noise levels using an adequate PBPK model and previous knowledge about the individual patient.

Voxel-based multimodel fitting method for modeling time activity curves in SPECT images

PURPOSE

Estimating the biodistribution and the pharmacokinetics from time-sequence SPECT images on a per-voxel basis is useful for studying activity nonuniformity or computing absorbed dose distributions by convolution of voxel kernels or Monte-Carlo radiation transport. Current approaches are either region-based, thus assuming uniform activity within the region, or voxel-based but using the same fitting model for all voxels.

METHODS

We propose a voxel-based multimodel fitting method (VoMM) that estimates a fitting function for each voxel by automatically selecting the most appropriate model among a predetermined set with Akaike criteria. This approach can be used to compute the time integrated activity (TIA) for all voxels in the image. To control fitting optimization that may fail due to excessive image noise, an approximated version based on trapezoid integration, named restricted method, is also studied. From this comparison, the number of failed fittings within images was estimated and analyzed. Numerical experiments were used to quantify uncertainties and feasibility was demonstrated with real patient data.

RESULTS

Regarding numerical experiments, root mean square errors of TIA obtained with VoMM were similar to those obtained with bi-exponential fitting functions, and were lower (< 5% vs.  > 10%) than with single model approaches that consider the same fitting function for all voxels. Failure rates were lower with VoMM and restricted approaches than with single-model methods. On real clinical data, VoMM was able to fit 90% of the voxels and led to less failed fits than single-model approaches. On regions of interest (ROI) analysis, the difference between ROI-based and voxel-based TIA estimations was low, less than 4%. However, the computation of the mean residence time exhibited larger differences, up to 25%.

CONCLUSIONS

The proposed voxel-based multimodel fitting method, VoMM, is feasible on patient data. VoMM leads organ-based TIA estimations similar to conventional ROI-based method. However, for pharmacokinetics analysis, studies of spatial heterogeneity or voxel-based absorbed dose assessment, VoMM could be used preferentially as it prevents model overfitting.

Individualized Dosimetry for Theranostics: Necessary, Nice to Have, or Counterproductive?

In 2005, the term theragnostics (theranostics) was introduced for describing the use of imaging for therapy planning in radiation oncology. In nuclear medicine, this expression describes the use of tracers for predicting the absorbed doses in molecular radiotherapy and, thus, the safety and efficacy of a treatment. At present, the most successful groups of isotopes for this purpose are ¹²³I/¹²⁴I/¹³¹I, ⁶⁸Ga/¹⁷⁷Lu, and ¹¹¹In/⁸⁶Y/⁹⁰Y. The purpose of this review is to summarize available data on the dosimetry and dose-response relationships of several theranostic compounds, with a special focus on radioiodine therapy for differentiated thyroid cancer and peptide receptor radionuclide therapy. These are treatment modalities for which dose-response relationships for healthy tissues and tumors have been demonstrated. In addition, available data demonstrate that posttherapeutic dosimetry after a first treatment cycle predicts the absorbed doses in further cycles. Both examples show the applicability of the concept of theranostics in molecular radiotherapies. Nevertheless, unanswered questions need to be addressed in clinical trials incorporating dosimetry-related concepts for determining the amount of therapeutic activity to be administered.

A single-source photon source model of a linear accelerator for Monte Carlo dose calculation

PURPOSE

To introduce a new method of deriving a virtual source model (VSM) of a linear accelerator photon beam from a phase space file (PSF) for Monte Carlo (MC) dose calculation.

MATERIALS AND METHODS

A PSF of a 6 MV photon beam was generated by simulating the interactions of primary electrons with the relevant geometries of a Synergy linear accelerator (Elekta AB, Stockholm, Sweden) and recording the particles that reach a plane 16 cm downstream the electron source. Probability distribution functions (PDFs) for particle positions and energies were derived from the analysis of the PSF. These PDFs were implemented in the VSM using inverse transform sampling. To model particle directions, the phase space plane was divided into a regular square grid. Each element of the grid corresponds to an area of 1 mm2 in the phase space plane. The average direction cosines, Pearson correlation coefficient (PCC) between photon energies and their direction cosines, as well as the PCC between the direction cosines were calculated for each grid element. Weighted polynomial surfaces were then fitted to these 2D data. The weights are used to correct for heteroscedasticity across the phase space bins. The directions of the particles created by the VSM were calculated from these fitted functions. The VSM was validated against the PSF by comparing the doses calculated by the two methods for different square field sizes. The comparisons were performed with profile and gamma analyses.

RESULTS

The doses calculated with the PSF and VSM agree to within 3% /1 mm (>95% pixel pass rate) for the evaluated fields.

CONCLUSION

A new method of deriving a virtual photon source model of a linear accelerator from a PSF file for MC dose calculation was developed. Validation results show that the doses calculated with the VSM and the PSF agree to within 3% /1 mm.

Prediction of time-integrated activity coefficients in PRRT using simulated dynamic PET and a pharmacokinetic model

PURPOSE

To investigate the accuracy of predicted time-integrated activity coefficients (TIACs) in peptide-receptor radionuclide therapy (PRRT) using simulated dynamic PET data and a physiologically based pharmacokinetic (PBPK) model.

METHODS

PBPK parameters were estimated using biokinetic data of 15 patients after injection of (152±15)MBq of ¹¹¹In-DTPAOC (total peptide amount (5.78±0.25)nmol). True mathematical phantoms of patients (MPPs) were the PBPK model with the estimated parameters. Dynamic PET measurements were simulated as being done after bolus injection of 150MBq ⁶⁸Ga-DOTATATE using the true MPPs. Dynamic PET scans around 35min p.i. (P₁), 4h p.i. (P₂) and the combination of P₁ and P₂ (P₃) were simulated. Each measurement was simulated with four frames of 5min each and 2 bed positions. PBPK parameters were fitted to the PET data to derive the PET-predicted MPPs. Therapy was simulated assuming an infusion of 5.1GBq of ⁹⁰Y-DOTATATE over 30min in both true and PET-predicted MPPs. TIACs of simulated therapy were calculated, true MPPs (true TIACs) and predicted MPPs (predicted TIACs) followed by the calculation of variabilities v.

RESULTS

For P₁ and P₂ the population variabilities of kidneys, liver and spleen were acceptable (v<10%). For the tumours and the remainders, the values were large (up to 25%). For P₃, population variabilities for all organs including the remainder further improved, except that of the tumour (v>10%).

CONCLUSION

Treatment planning of PRRT based on dynamic PET data seems possible for the kidneys, liver and spleen using a PBPK model and patient specific information.

Time-integrated activity coefficient estimation for radionuclide therapy using PET and a pharmacokinetic model: A simulation study on the effect of sampling schedule and noise

PURPOSE

The aim of this study was to investigate the accuracy of PET-based treatment planning for predicting the time-integrated activity coefficients (TIACs).

METHODS

The parameters of a physiologically based pharmacokinetic (PBPK) model were fitted to the biokinetic data of 15 patients to derive assumed true parameters and were used to construct true mathematical patient phantoms (MPPs). Biokinetics of 150 MBq (68)Ga-DOTATATE-PET was simulated with different noise levels [fractional standard deviation (FSD) 10%, 1%, 0.1%, and 0.01%], and seven combinations of measurements at 30 min, 1 h, and 4 h p.i. PBPK model parameters were fitted to the simulated noisy PET data using population-based Bayesian parameters to construct predicted MPPs. Therapy simulations were performed as 30 min infusion of (90)Y-DOTATATE of 3.3 GBq in both true and predicted MPPs. Prediction accuracy was then calculated as relative variability vorgan between TIACs from both MPPs.

RESULTS

Large variability values of one time-point protocols [e.g., FSD = 1%, 240 min p.i., vkidneys = (9 ± 6)%, and vtumor = (27 ± 26)%] show inaccurate prediction. Accurate TIAC prediction of the kidneys was obtained for the case of two measurements (1 and 4 h p.i.), e.g., FSD = 1%, vkidneys = (7 ± 3)%, and vtumor = (22 ± 10)%, or three measurements, e.g., FSD = 1%, vkidneys = (7 ± 3)%, and vtumor = (22 ± 9)%.

CONCLUSIONS

(68)Ga-DOTATATE-PET measurements could possibly be used to predict the TIACs of (90)Y-DOTATATE when using a PBPK model and population-based Bayesian parameters. The two time-point measurement at 1 and 4 h p.i. with a noise up to FSD = 1% allows an accurate prediction of the TIACs in kidneys.

Synthesis and preclinical evaluation of an Al18F radiofluorinated GLU-UREA-LYS(AHX)-HBED-CC PSMA ligand

PURPOSE

The aim of this study was to synthesize and preclinically evaluate an 18F-PSMA positron emission tomography (PET) tracer. Prostate-specific membrane antigen (PSMA) specificity, biodistribution, and dosimetry in healthy and tumor-bearing mice were determined.

METHODS

Several conditions for the labeling of 18F-PSMA-11 via 18F-AlF-complexation were screened to study the influence of reaction temperature, peptide amount, ethanol volume, and reaction time. After synthesis optimization, biodistribution and dosimetry studies were performed in C57BL6 mice. For proof of PSMA-specificity, mice were implanted with PSMA-negative (PC3) and PSMA-positive (LNCaP) tumors in contralateral flanks. Static and dynamic microPET/computed tomography (CT) imaging was performed.

RESULTS

Quantitative labeling yields could be achieved with >97 % radiochemical purity. The 18F-PSMA-11 uptake was more than 24-fold higher in PSMA-high LNCaP than in PSMA-low PC3 tumors (18.4 ± 3.3 %ID/g and 0.795 ± 0.260 %ID/g, respectively; p < 4.2e-5). Results were confirmed by ex vivo gamma counter analysis of tissues after the last imaging time point. The highest absorbed dose was reported for the kidneys. The maximum effective dose for an administered activity of 200 MBq was 1.72 mSv.

CONCLUSION

18F-PSMA-11 using direct labeling of chelate-attached peptide with aluminum-fluoride detected PSMA-expressing tumors with high tumor-to-liver ratios. The kidneys were the dose-limiting organs. Even by applying the most stringent dosimetric calculations, injected activities of up to 0.56 GBq are feasible.

Novel Approach to Repeated Arterial Blood Sampling in Small Animal PET: Application in a Test-Retest Study with the Adenosine A1 Receptor Ligand [(11)C]MPDX

PURPOSE

Small animal positron emission tomography (PET) can be used to detect small changes in neuroreceptor availability. This often requires rapid arterial blood sampling. However, current catheterization procedures do not allow repeated blood sampling. We have developed a procedure which allows arterial sampling on repeated occasions in the same animal.

PROCEDURES

Eleven male Wistar rats were two times catheterized via a superficial branch of a femoral artery and scanned with [(11)C]MPDX and blood sampling. PET images were co-registered to a magnetic resonance imaging (MRI) template. Regional tracer distribution volumes (V T) in the brain were calculated by the Logan analysis. The procedure was repeated after 1 week.

RESULTS

Surgery was successful in 90 % of the cases, and discomfort was minor. The V T data showed small differences between test and retest, low between subject variability, and a strong agreement between and within subjects.

CONCLUSION

Repeated quantitative imaging with a high reproducibility is possible with this approach.

Investigating the Effect of Ligand Amount and Injected Therapeutic Activity: A Simulation Study for 177Lu-Labeled PSMA-Targeting Peptides

In molecular radiotherapy with ¹⁷⁷Lu-labeled prostate specific membrane antigen (PSMA) peptides, kidney and/or salivary glands doses limit the activity which can be administered. The aim of this work was to investigate the effect of the ligand amount and injected activity on the tumor-to-normal tissue biologically effective dose (BED) ratio for ¹⁷⁷Lu-labeled PSMA peptides. For this retrospective study, a recently developed physiologically based pharmacokinetic model was adapted for PSMA targeting peptides. General physiological parameters were taken from the literature. Individual parameters were fitted to planar gamma camera measurements (¹⁷⁷Lu-PSMA I&T) of five patients with metastasizing prostate cancer. Based on the estimated parameters, the pharmacokinetics of tumor, salivary glands, kidneys, total body and red marrow was simulated and time-integrated activity coefficients were calculated for different peptide amounts. Based on these simulations, the absorbed doses and BEDs for normal tissue and tumor were calculated for all activities leading to a maximal tolerable kidney BED of 10 Gy_2.5/cycle, a maximal salivary gland absorbed dose of 7.5 Gy/cycle and a maximal red marrow BED of 0.25 Gy₁₅/cycle. The fits yielded coefficients of determination > 0.85, acceptable relative standard errors and low parameter correlations. All estimated parameters were in a physiologically reasonable range. The amounts (for 2⁵−2⁹ nmol) and pertaining activities leading to a maximal tumor dose, considering the defined maximal tolerable doses to organs of risk, were calculated to be 272±253 nmol (452±420 μg) and 7.3±5.1 GBq. Using the actually injected amount (235±155 μg) and the same maximal tolerable doses, the potential improvement for the tumor BED was 1–3 fold. The results suggest that currently given amounts for therapy are in the appropriate order of magnitude for many lesions. However, for lesions with high binding site density or lower perfusion, optimizing the peptide amount and activity might improve the tumor-to-kidney and tumor-to-salivary glands BED ratio considerably.

Sensitivity Analysis of a Physiologically Based Pharmacokinetic Model Used for Treatment Planning in Peptide Receptor Radionuclide Therapy

The aim of this work was to evaluate the sensitivity of time-integrated activity coefficients (TIACs) on the erroneously chosen prior knowledge in a physiologically based pharmacokinetic (PBPK) model used for treatment planning in peptide receptor radionuclide therapy (PRRT). Parameters of the PBPK model were fitted to the biokinetic data of 15 patients after the injection of (111)In-DTPAOC. The fittings were performed using fixed parameter values taken from literature as prior knowledge (reference case, Ref). The fixed parameters were gender, physical information (e.g., body weight), dissociation rate koff, dissociation constant KD, fraction of blood flow, and spleen and liver volumes. The fittings were repeated with changed fixed parameters (Changed). The relative deviations (RDs) of TIACs calculated from Changed and Ref were analyzed for kidneys, tumor, liver, spleen, remainder, whole body, and serum. A changed koff has the largest effect on RD, the largest RD values were found for changed koff = 0.001 L/min: RDkidneys = (3 ± 3)%, RDtumor = (0.5 ± 4)%, RDliver = (6 ± 9)%, RDspleen = (5 ± 5)%, RDremainder = (2 ± 31)%, RDserum = (-4 ± 25)%, and RDwholebody = (3 ± 16)%. For other changed parameters, the maximum RDs were <1%. The calculation of organ TIACs in PRRT using the PBPK model was little affected by assigning wrong prior knowledge to the evaluated patients. The calculation of bone marrow-absorbed doses could be affected by the inaccurate TIACs of serum and remainder in the case of an inadequate koff.

The Influence of Early Measurements Onto the Estimated Kidney Dose in [(177)Lu][DOTA(0),Tyr(3)]Octreotate Peptide Receptor Radiotherapy of Neuroendocrine Tumors

PURPOSE

Multiple measurements have been required to estimate the radiation dose to the kidneys resulting from [(177)Lu]DOTATATE therapy for neuroendocrine tumors. The aim of this study was to investigate the influence of early time-point measurement in the renal dose calculation.

PROCEDURES

Anterior/posterior whole-body planar scintigraphy images were acquired at approx. 1, 24, 48, and 72 h after administration of [(177)Lu]DOTATATE. Furthermore, we acquired planar 1-bed dynamic recordings in 12 frames (5 min each) during the first hour. We assessed kidney exposure with a three-phase model consisting of a linear increase to the maximum within the initial minutes p.i., followed a bi-exponential decline. This three-phase-model served as reference for evaluating accuracy of dose estimates in 105 kidneys calculated by conventional mono-exponential fitting of the final three and four whole-body images.

RESULTS

Mean effective half-life times for the reference model were 25.8 ± 12.0 min and 63.9 ± 17.6 h, predicting a mean renal dose of 5.7 ± 2.1 Gy. The effective half-life time was 46.3 ± 15.4 h for the last four and 63.3 ± 17.0 h for the last three data points. The mean start of the first whole-body measurement was 1.2 ± 0.1 h p.i. The ratio of fast to slow phases was 28.1 ± 23.9% at this time point, which caused a mean absolute percentage dose deviation of 12.4% for four data points, compared to 3.1% for three data points. At a mean time of 2.4 h p.i. (max 5.1 h), the ratio of fast to slow phase declined below 5%.

CONCLUSIONS

Kinetic analysis of renal uptake using dynamic planar scans from the first hour after injection revealed a fast and a slow washout phase. Although the fast phase did not contribute substantially to the estimated renal dose, it could influence planar measurements performed within the first hours. We found that the presence of two clearance phases can hamper accurate dose estimation based on a single-phase model, resulting in approximately 12.4% dose underestimation, thus potentially resulting in overtreatment. In the absence of dynamic initial recordings, the first dosimetry measurements should therefore be obtained later than 3-5 h after [(177)Lu]DOTATATE injection. Omitting the early whole-body image reduced the dose estimation error to 3.1%.

177Lu-OPS201 targeting somatostatin receptors: in vivo biodistribution and dosimetry in a pig model

Background

¹⁷⁷Lu is used in peptide receptor radionuclide therapies for the treatment of neuroendocrine tumors. Based on the recent literature, SST2 antagonists are superior to agonists in tumor uptake. The compound OPS201 is the novel somatostatin antagonist showing the highest SST2 affinity. The aim of this study was to measure the in vivo biodistribution and dosimetry of ¹⁷⁷Lu-OPS201 in five anesthetized Danish Landrace pigs as an appropriate substitute for humans to quantitatively assess the absorbed doses for future clinical applications.

Results

¹⁷⁷Lu-OPS201 was obtained with a specific activity ranging from 10 to 17 MBq/μg. Prior to administration, the radiochemical purity was measured as s > 99.7 % in all cases. After injection, fast clearance of the compound from the blood stream was observed. Less than 5 % of the injected activity was presented in blood 10 min after injection. A series of SPECT/CT and whole-body scans conducted until 10 days after intravenous injection showed uptake mostly in the liver, spine, and kidneys. There was no visible uptake in the spleen. Blood samples were taken to determine the time-activity curve in the blood. Time-activity curves and time-integrated activity coefficients were calculated for the organs showing visible uptake. Based on these data, the absorbed organ dose coefficients for a 70-kg patient were calculated with OLINDA/EXM. For humans after an injection of 5 GBq ¹⁷⁷Lu-OPS201, the highest predicted absorbed doses are obtained for the kidneys (13.7 Gy), the osteogenic cells (3.9 Gy), the urinary bladder wall (1.8 Gy), and the liver (1.0 Gy). No metabolites of ¹⁷⁷Lu-OPS201 were found by radio HPLC analysis. None of the absorbed doses calculated will exceed organ toxicity levels.

Conclusions

The ¹⁷⁷Lu-OPS201 was well tolerated and caused no abnormal physiological or behavioral signs. In vivo distributions and absorbed doses of pigs are comparable to those observed in other publications. According to the biodistribution data in pigs, presented in this work, the expected radiation exposure in humans will be within the acceptable range.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-016-0204-9) contains supplementary material, which is available to authorized users.