Improving individualised dosimetry in radioiodine therapy for hyperthyroidism using population biokinetic modelling

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.

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.

Dosimetry during adjuvant 131I therapy in patients with differentiated thyroid cancer-clinical implications

The activity of radioiodine (¹³¹I) used in adjuvant therapy for thyroid cancer ranges between 30 mCi (1.1 GBq) and 150 mCi (5.5 GBq). Dosimetry based on Marinelli's formula, taking into consideration the absorbed dose in the postoperative tumour bed (D) should systematise the determination of ¹³¹I activity. Retrospective analysis of 57 patients with differentiated thyroid cancer (DTC) after thyreidectomy and adjuvant ¹³¹I therapy with the fixed activity of 3.7 GBq. In order to calculate D from Marinelli's formula, the authors took into account, among other things, repeated dosimetry measurements (after 6, 24, and 72 h) made during scintigraphy and after administration of the therapeutic activity or radioiodine. In 75% of the patients, the values of D were > 300 Gy (i.e. above the value recommended by current guidelines). In just 16% of the patients, the obtained values fell between 250 and 300 Gy, whereas in 9% of the patients, the value of D was < 250 Gy. The therapy was successful for all the patients (stimulated Tg < 1 ng/ml and ¹³¹I uptake < 0.1% in the thyroid bed in follow-up examination). Dosimetry during adjuvant ¹³¹I therapy makes it possible to diversify the therapeutic activities of ¹³¹I in order to obtain a uniform value of D.