Skip to main content

Main menu

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI

User menu

  • Subscribe
  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Nuclear Medicine
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI
  • Subscribe
  • My alerts
  • Log in
  • My Cart
Journal of Nuclear Medicine

Advanced Search

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • View or Listen to JNM Podcast
  • Visit JNM on Facebook
  • Join JNM on LinkedIn
  • Follow JNM on Twitter
  • Subscribe to our RSS feeds
Research ArticleClinical Investigation

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

Deni Hardiansyah, Elham Yousefzadeh-Nowshahr, Felix Kind, Ambros J. Beer, Juri Ruf, Gerhard Glatting and Michael Mix
Journal of Nuclear Medicine April 2024, 65 (4) 566-572; DOI: https://doi.org/10.2967/jnumed.123.266268
Deni Hardiansyah
1Medical Physics and Biophysics, Physics Department, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok, Indonesia;
2Medical Radiation Physics, Department of Nuclear Medicine, Ulm University, Ulm, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elham Yousefzadeh-Nowshahr
2Medical Radiation Physics, Department of Nuclear Medicine, Ulm University, Ulm, Germany;
3Department of Nuclear Medicine, Ulm University, Ulm, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Felix Kind
4Department of Nuclear Medicine, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ambros J. Beer
3Department of Nuclear Medicine, Ulm University, Ulm, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Juri Ruf
4Department of Nuclear Medicine, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gerhard Glatting
2Medical Radiation Physics, Department of Nuclear Medicine, Ulm University, Ulm, Germany;
3Department of Nuclear Medicine, Ulm University, Ulm, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael Mix
4Department of Nuclear Medicine, Medical Center – University of Freiburg, Faculty of Medicine, University of Freiburg, Freiburg, Germany; and
5Nuclear Medicine Division, Department of Medical Imaging and Clinical Oncology, Faculty of Medicine and Health Science, Stellenbosch University, Cape Town, South Africa
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF
Loading

Visual Abstract

Figure
  • Download figure
  • Open in new tab
  • Download powerpoint

Abstract

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 Embedded Image 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.

  • STP dosimetry
  • Akaike weight
  • model selection
  • NLME modeling
  • PSMA

Molecular radiotherapy targeting the prostate-specific membrane antigen (PSMA) is a promising palliative treatment option in patients with metastatic castration-resistant prostate cancer (1), using the overexpression of the PSMA in prostate cancer cells (2). [177Lu]Lu-PSMA-617 has shown high treatment efficacy and tolerability in both phase II (3–5) and phase III (6) trials. Although [177Lu]Lu-PSMA-617 therapy is generally well tolerated, several organs, such as kidneys, salivary glands, and bone marrow, are considered organs at risk (7). It has been reported that renal and hematopoietic toxicity are the main side effects of [177Lu]Lu-PSMA therapy and peptide receptor radionuclide therapy (8). Furthermore, the established kidney tolerance dose of 23 Gy based on data from external radiation therapy (5) constitutes the main dose-limiting factor for repeated treatment over multiple cycles. Thus, various protocols for kidney protection using amino acids (9) and dosimetry during [177Lu]Lu-PSMA therapy have been developed.

Calculation of the absorbed dose in molecular radiotherapy dosimetry is based on the time-integrated activity (TIA) and the dose factor (10). Estimating individual absorbed doses in molecular radiotherapy is desirable for therapy optimization and kidney protection. However, the individual absorbed dose is often not calculated in routine clinical practice, as it needs repeated imaging to determine the TIA (11,12). This repeated imaging corresponds to a high workload and additional patient burden, requiring several patient visits. Therefore, a method with fewer imaging sessions is highly desirable.

The feasibility of using a reduced number of imaging data per patient for TIA determination has been investigated for [177Lu]Lu-PSMA therapy (13,14), peptide receptor radionuclide therapy (15–24), radioiodine therapy (25,26), and radioimmunotherapy (27). One of the favorable approaches for dosimetry based on single-time-point (STP) imaging data is the nonlinear mixed-effects (NLME) method (28). For STP dosimetry, it has been shown (17) that the application of NLME modeling resulted in a low number of TIAs with an error exceeding 10% compared with the all-time-point (ATP) TIAs. The percentage of TIAs with a 10% error was approximately 3 times lower than that using the dose mapping methods introduced by Hänscheid et al. (18) and the curve-fitting methods introduced by Madsen et al. (23). Recently, the implementation of NLME and a physiologically based pharmacokinetic model in STP dosimetry could lead to an accurate determination of TIAs in various organs, including kidneys and tumors (11).

Although NLME modeling in STP dosimetry has shown promising results, it has thus far been implemented only in populations with relatively few patients (i.e., 10 (17) or fewer (11) patients). Implementation of the technique in larger patient populations is needed to validate the concept and increase the statistical significance of the accuracy of STP dosimetry with NLME modeling. In addition, identifying the optimal fit function through a population-based model selection (PBMS) (29) with large patient populations is highly desirable for increasing reproducibility. Therefore, in this study, we implemented a PBMS method in the NLME framework to identify the best model describing the [177Lu]Lu-PSMA-617 biokinetic data in kidneys and performed STP dosimetry using NLME modeling in a larger population of 63 patients.

MATERIALS AND METHODS

Biokinetic Data

In brief, biokinetic data of [177Lu]Lu-PSMA-617 in kidneys from 63 patients with metastasized castration-resistant prostate cancer in the first treatment cycle were used in this study (8,30). A [177Lu]Lu-PSMA-617 activity with a median of 6 GBq/cycle was administered intravenously to the patients, followed by SPECT/CT imaging measurements (mean ± SD) at time points 1 (1.8 ± 0.8 h), 2 (18.7 ± 0.9 h), 3 (42.6 ± 1.0 h), 4 (66.2 ± 0.9 h), and 5 (160.3 ± 24.2 h) after injection. Acquisitions were performed on a SPECT/CT scanner (BrightView XCT; Philips Healthcare) using a medium-energy general-purpose collimator. An energy window of ±10% around the 208-keV peak was used for the measurements. Cone-beam CT scans with 30 mAs at 120 kV were used for attenuation correction. Scatter was corrected during reconstruction with the vendor implementation of the effective-scatter source estimation method (31). SPECT was reconstructed using the ordered-subsets expectation maximization algorithm without resolution recovery (8). The calibration factor for 177Lu was measured using the National Electrical Manufacturers Association image-quality phantom. The calculated calibration factor for the imaging protocol used was 9.9 ± 0.4 cps/MBq. Kidney segmentation and analysis were done using the software package Rover ABX (ABX Advanced Biochemical Compounds GmbH) (32). The kidneys were delineated on CT images. At each time point, the activity was measured in both the left and the right kidneys. To compensate for spill-out effects caused by the limited resolution of the SPECT system, all segmented kidney volumes of interest were extended by 2 SPECT voxels in all directions. The enlargement factor was determined by a phantom study. Finally, the fraction of injected activity in both kidneys was used in this study. The institutional review board (vote no. 326/18) approved the retrospective study, and all subjects gave written informed consent.

Fit Function Set: Sums-of-Exponential (SOE) Functions

The following SOE functions were used to describe the biokinetics of [177Lu]Lu-PSMA-617 in kidneys (Eqs. 1–13).Embedded Image Eq. 1Embedded Image Eq. 2Embedded Image Eq. 3Embedded Image Eq. 4Embedded Image Eq. 5Embedded Image Eq. 6Embedded Image Eq. 7Embedded Image Eq. 8Embedded Image Eq. 9Embedded Image Eq. 10Embedded Image Eq. 11Embedded Image Eq. 12Embedded Image Eq. 13where Embedded Image is a fit function, withEmbedded Image being the total number of estimated parameters; Embedded Image are the coefficients of the respective exponential terms; Embedded Image are the biologic clearance or uptake rates of the radiopharmaceutical in the kidneys;Embedded Image is the rate of blood circulation set to Embedded Image; and Embedded Image is the physical decay constant of 177Lu (Embedded Image, with the 177Lu half-life of 6.6443 d (33)). The Embedded Image and Embedded Image parameters were constrained to have values greater than 0.

Prior knowledge of radiopharmaceutical kinetics was used as the foundation of the chosen functions and constraints in this study, as suggested by Burnham et al. (34). Therefore, the following were implemented: only SOE functions were considered (35,36); a constraint of Embedded Image was used, except for Embedded Image and Embedded Image; a rapid increase in activity in the kidneys with a half-life of 1 min (ln(2)/Embedded Image) was added, which is caused by the blood circulation time in humans (29,37); and the radioactive decay was considered as an exponential factor (38–40). The functions Embedded Image and Embedded Image were added for comparison, as these mono- and biexponential functions, respectively, are often used in the literature (18,23,35).

Study Workflow

Figure 1 shows the workflow of the study. The parameters of the SOE functions (Eqs. 1–13) were fitted to the ATP biokinetic data of [177Lu]Lu-PSMA-617 in kidneys (Biokinetic Data section) using the NLME method (Supplemental Eqs. 1–3) (41,42). In this study, parameters of the SOE functions in Equations 1–13 were modeled as a combination of fixed and random effects. NLME model fittings and simulations were performed in MATLAB software, version R2020a (https://www.mathworks.com/help/stats/nlmefitsa.html). The PBMS NLME approach was used to select the best-fit SOE function describing the data (43). The fittings were randomly repeated 1,000 times to find the best-fit starting values. The SOE function passing the goodness-of-fit test (Supplemental Eqs. 4–6) (35,44,45,47) with the highest Akaike weight was selected as the fit function most supported by the data (35,43,48).

FIGURE 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1.

Workflow of study. PBMS NLME approach was used to select best fit function based on goodness of fit and Akaike weight. Reference absorbed doses were based on ATP fittings using best-fit function. STP absorbed doses were calculated analogously using only STP data for investigated patient. RDs and RMSEs were used to analyze accuracy of STP absorbed doses against reference absorbed doses. AICc = corrected Akaike information criterion; i = time point; N = total number of patients; NLMEM = PBMS NLME; rAD = reference absorbed dose; rTIA = reference TIA; sAD = STP absorbed dose; sTIA = STP TIA.

The best-fit SOE function from the model selection and ATP fitting was then used to calculate the reference TIA (Fig. 1). Calculation of the TIAs was based on the analytic solution of the SOE functions. In STP dosimetry, the parameters of a specific patient with STP data were fitted simultaneously to the STP data at different time points of that patient with all data points of all other patients.

The STP dosimetry methods introduced by Hänscheid et al. and Madsen et al., which have been widely adopted in STP studies (20,22,49–51), were implemented in this study for comparison. The methods were initially used to describe the biokinetics of [177Lu]Lu-DOTATATE/TOC and [90Y]Y-DOTATOC (18,23). In this study, the methods were assumed to be applicable for [177Lu]Lu-PSMA-617. The corresponding TIAs for the Hänscheid and Madsen methods were calculated using the following equations:Embedded Image Eq. 14Embedded Image Eq. 15where Embedded Image is the activity measured at time point i, Embedded Image is the time of measurement i, and Embedded Image is the population effective clearance rate constant calculated using the leave-one-out method.

The S value of 2.93 × 10−5 Gy × min−1 × MBq−1 (46,52) for male patients was used to calculate the absorbed dose using NLME ATP fitting (reference absorbed doses), NLME STP fitting (STP absorbed doses), and the STP Hänscheid and Madsen methods.

Relative deviations (RDs) and root-mean-square errors (RMSEs) were used to analyze the accuracy of the calculated STP absorbed doses, STP Hänscheid method, and STP Madsen method with the reference absorbed doses as the reference. The RDs and the RMSEs were calculated according toEmbedded Image Eq. 16Embedded Image Eq. 17where sAD is STP absorbed dose, rAD is reference absorbed dose, Embedded Image is the RD of the STP method at time point k, Embedded Image is the root-mean square of Embedded Image over all patients, and Embedded Image is the SD of Embedded Image. The reference TIA was computed through population fitting, using a notably high parameter-to-data ratio within our dataset, specifically 315 data points and 13 parameters (Table 1). In contrast, individual patient fitting using a monoexponential function as in the Hänscheid and Madsen methods uses a significantly lower ratio, such as 3 data points and 2 parameters. This results in a higher uncertainty of the resulting TIA.

View this table:
  • View inline
  • View popup
TABLE 1.

PBMS with NLME Model Results for ATP fitting

RESULTS

The NLME fittings for functions f4b, f5a, Embedded Image, Embedded Image,Embedded Image, Embedded Image, Embedded Image, and Embedded Image did not pass the goodness-of-fit test. On the basis of the PBMS NLME method, the function Embedded Image was selected as the best-fit function, with an Akaike weight of 100% (Table 1; Supplemental Fig. 1; supplemental materials are available at http://jnm.snmjournals.org). The estimated fixed-effect values are presented in Table 2. The estimated random-effect values ranged from 0.08 to 0.52. The intraindividual variability (a in Supplemental Eq. 3) was 0.079. Figure 2 shows the percentage RD of the TIAs from ATP fitting with functions f2, f3, f4a, and f4c to the TIAs from ATP fitting with function f6a. The mean percentage RDs in Figure 2 are 1.3% ± 14.4%, 1.1% ± 14.2%, 0.1% ± 1.9%, and 0.2% ± 2.0% for functions f2, f3, f4a, and f4c, respectively.

View this table:
  • View inline
  • View popup
TABLE 2.

Parameters Estimated from ATP Fitting Obtained from Best Function Derived Using PBMS with NLME Modeling (i.e., Embedded Image)

FIGURE 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2.

RD values of TIAs obtained from ATP fitting using functions f2, f3, f4a, and f4c to TIAs obtained from ATP fitting using best function according to PBMS NLME method, that is, f6a. Box plots show median, 25th and 75th percentiles, whiskers at 5th and 95th percentiles, and outliers. Percentage RDs for functions f4a or f4c are relatively small because of marginal contribution of uptake phase. Function f3 was added for comparison in Figure 2, as this function is often used in the literature (17,35,39). For broad range of covered sampling times after injection in population investigated here, single biologic elimination phase (as is also case for f2) is insufficient.

Figure 3 shows the population RD of the absorbed dose for STP dosimetry at different time points. The mean of the population RD for STP dosimetry is shown in Supplemental Table 1. The medians of the population RD for STP dosimetry at time points 1–5 were −2% (range, −55% to 59%), −4% (range, −53% to 10%), −0.7% (range, −39% to 20%), −0.1% (range, −32% to 106%), and 0.4% (range, −40% to 348%), respectively. The RMSE values of the absorbed dose RDs for time points 1–5 were 23%, 16%, 10%, 20%, and 53%, respectively. Thus, the STP dosimetry at time point 3 has the lowest mean RD and RMSE values.

FIGURE 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3.

RD values of TIAs obtained using PBMS NLME (NLMEM) function Embedded Image (Eq. 8), Hänscheid method (Eq. 14), and Madsen method (Eq. 15) for STP dosimetry at time points 1–5. TIAs from ATP fittings calculated using PBMS NLME function f6a were used as reference values. Box plots show median, 25th and 75th percentiles, whiskers at 2.5th and 97.5th percentiles, and outliers. TP = time point.

Use of STP dosimetry with the PBMS NLME method outperformed the Hänscheid and Madsen methods for all time points (Fig. 3; Supplemental Table 1). The STP NLME model at time point 3 could not accurately calculate the TIA in a few outlier patients (outside the mean ± 2 SDs of the RD), that is, patients 12, 36, 41, and 44, with absorbed dose RDs of 26%, 39%, 23% and 29%, respectively (Fig. 4).

FIGURE 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 4.

Time–activity curves from ATP (green lines) and STP (black lines) fittings at time point 3 simulated with function f6a for 4 outlier patients. Estimated individual parameters indicate that these outliers had lowest biologic clearance Embedded Image (4.8 × 10−3/h to 6.7 × 10−3/h) in population and were relatively lower than fixed effect (9.1 × 10−3/h) (Table 2), a well-known behavior in STP dosimetry (17). P = patient.

DISCUSSION

In this study, 315 biokinetic data were used and the SOE function Embedded Image was selected from the PBMS NLME method based on the goodness-of-fit test and Akaike weight of 100% (Table 1). The SOE function Embedded Image shows 2 uptake and clearance phases of the radiopharmaceutical in kidneys. This function shows a higher parameter-to-data ratio in our patient population than the frequently used monoexponential function in the literature (Fig. 2) (18,23). The total number of exponential terms and parameters in this study for [177Lu]Lu-PSMA-617 (4 exponentials and 6 parameters) is higher than found in our recent study for [177Lu]Lu-PSMA-I&T (Embedded Image [2 exponentials and 3 parameters]) (43). These differences can be attributed to the total number of data and the schedule and range of the measurements between the studies, that is, 315 biokinetic data of [177Lu]Lu-PSMA-617 measured between 0.7 h after injection and 236 h after injection (this study) versus 46 biokinetic data of [177Lu]Lu-PSMA-I&T measured between 0.5 h after injection and 168 h after injection (43). This shows that adding more data covering a broader time range is beneficial as it allows the identification of a function with a larger number of parameters. This function should in turn be able to better reflect the biokinetics of the radiopharmaceutical.

It has been shown that STP data at 48 h after injection would lead to a relatively accurate prediction of kidney TIAs for [177Lu]Lu-PSMA (20). For [177Lu]Lu-PSMA-I&T therapy, simulation results by Rinscheid et al. showed STP data at around 52 h after injection as the optimal time point for calculating kidney TIAs (14). This work identified time points 2 and 3, having a relatively low mean and median RD of absorbed doses (Fig. 3) in our [177Lu]Lu-PSMA-617 dataset. Time point 3 (42.6 ± 1.0 h) was identified as having the lowest RMSE value. Time point 3 in our study for kidneys is consistent with the optimal time point for STP dosimetry reported in the literature (14,20).

The percentages of the total number of patients with absolute RDs higher than 10% for STP dosimetry using time points 1–5 were 59% (n = 37), 40% (n = 25), 24% (n = 15), 35% (n = 22), and 38% (n = 24), respectively (Supplemental Table 1). Time point 3 was selected as optimal for STP dosimetry as it has the lowest RMSE and mean RD value. Therefore, PBMS NLME and STP imaging of [177Lu]Lu-PSMA-617 kinetics might be used to determine absorbed doses in kidneys accurately, with only 4 outliers (RDs outside the threshold of mean ± 2 SDs) at time point 3 (Fig. 4).

The fractions of percentage RDs lower than 10% for STP dosimetry at time point 3 using the methods of Hänscheid et al. and Madsen et al. were 0.44 and 0.41, respectively, whereas the best model according to the PBMS NLME approach demonstrated that the fraction of percentage RDs lower than 10% for STP dosimetry at time point 3 was 0.76. This finding aligns with the literature, suggesting that NLME modeling can reduce the number of patients with a percentage RD lower than 10% calculated using the Hänscheid and Madsen methods (17). One contributing factor to these results is the use of monoexponential functions (f2) (18,23). The monoexponential function had a lower quality in describing the biokinetic data in our patient population than the best model Embedded Image derived from PBMS NLME (Fig. 2; Supplemental Table 1). The superiority of our method in this study may also be attributed to the different radiopharmaceuticals (including the different functions used) and the different study populations.

In this study, calculation of the absorbed dose was based on the self-dose. The kidney non–self-dose has been reported to contribute up to 10% to the kidney dose for [177Lu]Lu-DOTATATE (53). Therefore, further study is needed to investigate how to add the non–self-dose using STP dosimetry with the PBMS NLME method.

NLME modeling uses a population pharmacokinetic approach, wherein instead of individually fitting the data from each individual, all data within the population are simultaneously fitted. The benefit of using NLME models is that they do not require the data to be comprehensive or densely sampled, as is often the case in nonlinear regression analyses with individual data points (41). Thus, the reference TIA was computed through population fitting, using a notably higher parameter-to-data ratio within our dataset, specifically 315 data points and 13 parameters (Table 1). In contrast, individual patient fitting using a monoexponential function uses a significantly lower ratio, such as 3 data points and 2 parameters. This results in a higher uncertainty of the resulting TIA.

In a clinical context, for a new radiopharmaceutical or a new patient population, the STP method introduced here plays a crucial role in the following sequence of steps. The first step is to collect biokinetic activity data during therapy from the patient population of interest. The second is to determine the most suitable fit function based on the data using the PBMS NLME method presented here. This will yield the corresponding fixed-effect and random-effect parameters (Table 2). The third is to conduct STP dosimetry for new patients by applying NLME modeling, using the identified optimal function along with the fixed-effect and random-effect parameters determined in step 2.

Thus, for a known radiopharmaceutical and patient population, step 3 is sufficient, as once the data are obtained in step 2, these can be used for STP dosimetries in every center. Therefore, our results can be used for other studies if the radiopharmaceutical and patient populations are equivalent and appropriate image quantification is ensured.

A prerequisite of NLME modeling is the need for a historical cohort with multiple-time-point [177Lu]Lu-PSMA-617 SPECT/CT data (step 1). Because of the reimbursement policies, the availability of posttreatment SPECT/CT imaging is increasing (17). Additionally, promoting data exchange among different health care centers, exemplified by the information offered in this publication, can streamline the model development process.

CONCLUSION

In this study, we developed and implemented an STP dosimetry using NLME modeling and PBMS to calculate kidney-absorbed doses of [177Lu]Lu-PSMA-617. Our results show that a single measurement of SPECT/CT at 2 d after injection might be used to calculate accurate kidney-absorbed doses.

KEY POINTS

QUESTION: Can a nonlinear fixed-effects model with PBMS using many biokinetic data covering a broad postinjection time range provide accurate STP-based kidney-absorbed dose estimates for patients treated with [177Lu]Lu-PSMA-617?

PERTINENT FINDINGS: A single measurement of SPECT/CT at 2 d after injection may be used to calculate accurate kidney-absorbed doses using NLME modeling. The mean RD for STP at 2 d was −4% ± 9%, with only 24% of the total number of patients having RDs higher than 10% (15/63 patients) and only 10% of the total number of patients having RDs higher than 20% (6/63 patients).

IMPLICATIONS FOR PATIENT CARE: The STP method presented in this study holds great potential to be integrated into routine clinical practice, thereby increasing the acceptance and use of individualized dosimetry, which will benefit patients by predicting related risks.

DISCLOSURE

This work was supported by a PUTI Q1 grant 2023 from Universitas Indonesia (grant NKB-467/UN2.RST/HKP.05.00/2023). No other potential conflict of interest relevant to this article was reported.

Footnotes

  • ↵* Contributed equally to this work.

  • Published online Feb. 29, 2024.

  • © 2024 by the Society of Nuclear Medicine and Molecular Imaging.

REFERENCES

  1. 1.↵
    1. Afshar-Oromieh A,
    2. Babich JW,
    3. Kratochwil C,
    4. et al
    . The rise of PSMA ligands for diagnosis and therapy of prostate cancer. J Nucl Med. 2016;57(suppl 3):79S–89S.
    OpenUrl
  2. 2.↵
    1. Ghosh A,
    2. Heston WD
    . Tumor target prostate specific membrane antigen (PSMA) and its regulation in prostate cancer. J Cell Biochem. 2004;91:528–539.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Violet J,
    2. Sandhu S,
    3. Iravani A,
    4. et al
    . Long-term follow-up and outcomes of retreatment in an expanded 50-patient single-center phase II prospective trial of 177Lu-PSMA-617 theranostics in metastatic castration-resistant prostate cancer. J Nucl Med. 2020;61:857–865.
    OpenUrlAbstract/FREE Full Text
  4. 4.
    1. Rahbar K,
    2. Ahmadzadehfar H,
    3. Kratochwil C,
    4. et al
    . German multicenter study investigating 177Lu-PSMA-617 radioligand therapy in advanced prostate cancer patients. J Nucl Med. 2017;58:85–90.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Hofman MS,
    2. Emmett L,
    3. Sandhu S,
    4. et al
    . [177Lu]Lu-PSMA-617 versus cabazitaxel in patients with metastatic castration-resistant prostate cancer (TheraP): a randomised, open-label, phase 2 trial. Lancet. 2021;397:797–804.
    OpenUrlCrossRef
  6. 6.↵
    1. Sartor O,
    2. de Bono J,
    3. Chi KN,
    4. et al
    . 177Lutetium–PSMA-617 for metastatic castration-resistant prostate cancer. N Engl J Med. 2021;385:1091–1103.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Kratochwil C,
    2. Fendler WP,
    3. Eiber M,
    4. et al
    . EANM procedure guidelines for radionuclide therapy with 177Lu-labelled PSMA-ligands (177Lu-PSMA-RLT). Eur J Nucl Med Mol Imaging. 2019;46:2536–2544.
    OpenUrl
  8. 8.↵
    1. Mix M,
    2. Renaud T,
    3. Kind F,
    4. et al
    . Kidney doses in 177Lu-based radioligand therapy in prostate cancer: is dose estimation based on reduced dosimetry measurements feasible? J Nucl Med. 2022;63:253–258.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Bodei L,
    2. Mueller-Brand J,
    3. Baum RP,
    4. et al
    . The joint IAEA, EANM, and SNMMI practical guidance on peptide receptor radionuclide therapy (PRRNT) in neuroendocrine tumours. Eur J Nucl Med Mol Imaging. 2013;40:800–816.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Bolch WE,
    2. Eckerman KF,
    3. Sgouros G,
    4. Thomas SR
    . MIRD pamphlet no. 21: a generalized schema for radiopharmaceutical dosimetry—standardization of nomenclature. J Nucl Med. 2009;50:477–484.
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. Hardiansyah D,
    2. Riana A,
    3. Beer A,
    4. Glatting G
    . Single-time-point estimation of absorbed doses in PRRT using a non-linear mixed-effects model. Z Med Phys. 2023;33:70–81.
    OpenUrl
  12. 12.↵
    1. Zvereva A,
    2. Kamp F,
    3. Schlattl H,
    4. Zankl M,
    5. Parodi K
    . Impact of interpatient variability on organ dose estimates according to MIRD schema: uncertainty and variance-based sensitivity analysis. Med Phys. 2018;45:3391–3403.
    OpenUrl
  13. 13.↵
    1. Jackson PA,
    2. Hofman MS,
    3. Hicks RJ,
    4. Scalzo M,
    5. Violet JA
    . Radiation dosimetry in 177Lu-PSMA-617 therapy using a single posttreatment SPECT/CT scan: a novel methodology to generate time- and tissue-specific dose factors. J Nucl Med. 2020;61:1030–1036.
  14. 14.↵
    1. Rinscheid A,
    2. Kletting P,
    3. Eiber M,
    4. Beer AJ,
    5. Glatting G
    . Influence of sampling schedules on 177Lu-PSMA dosimetry. EJNMMI Phys. 2020;7:41.
    OpenUrl
  15. 15.↵
    1. Chicheportiche A,
    2. Sason M,
    3. Zidan M,
    4. et al
    . Impact of single-time-point estimates of 177Lu-PRRT absorbed doses on patient management: validation of a trained multiple-linear-regression model in 159 patients and 477 therapy cycles. J Nucl Med. 2023;64:1610–1616.
    OpenUrlAbstract/FREE Full Text
  16. 16.
    1. Wang C,
    2. Peterson AB,
    3. Wong KK,
    4. Roseland ME,
    5. Schipper MJ,
    6. Dewaraja YK
    . Single-time-point imaging for dosimetry after [177Lu]Lu-DOTATATE: accuracy of existing methods and novel data-driven models for reducing sensitivity to time-point selection. J Nucl Med. 2023;64:1463–1470.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Devasia TP,
    2. Dewaraja YK,
    3. Frey KA,
    4. Wong KK,
    5. Schipper MJ
    . A novel time-activity information sharing approach using nonlinear mixed models for patient-specific dosimetry with reduced imaging time points: application in SPECT/CT imaging post 177Lu-DOTATATE. J Nucl Med. 2021;62:1118–1125.
    OpenUrl
  18. 18.↵
    1. Hänscheid H,
    2. Lapa C,
    3. Buck AK,
    4. Lassmann M,
    5. Werner RA
    . Dose mapping after endoradiotherapy with 177Lu-DOTATATE/DOTATOC by a single measurement after 4 days. J Nucl Med. 2018;59:75–81.
    OpenUrlAbstract/FREE Full Text
  19. 19.
    1. Hardiansyah D,
    2. Maass C,
    3. Attarwala AA,
    4. et al
    . The role of patient-based treatment planning in peptide receptor radionuclide therapy. Eur J Nucl Med Mol Imaging. 2016;43:871–880.
    OpenUrl
  20. 20.↵
    1. Hou X,
    2. Brosch J,
    3. Uribe C,
    4. et al
    . Feasibility of single-time-point dosimetry for radiopharmaceutical therapies. J Nucl Med. 2021;62:1006–1011.
    OpenUrlAbstract/FREE Full Text
  21. 21.
    1. Maaß C,
    2. Sachs JP,
    3. Hardiansyah D,
    4. Mottaghy FM,
    5. Kletting P,
    6. Glatting G
    . Dependence of treatment planning accuracy in peptide receptor radionuclide therapy on the sampling schedule. EJNMMI Res. 2016;6:30.
    OpenUrl
  22. 22.↵
    1. Sundlöv A,
    2. Gustafsson J,
    3. Brolin G,
    4. et al
    . Feasibility of simplifying renal dosimetry in 177Lu peptide receptor radionuclide therapy. EJNMMI Phys. 2018;5:12.
    OpenUrl
  23. 23.↵
    1. Madsen MT,
    2. Menda Y,
    3. O’Dorisio TM,
    4. O’Dorisio MS
    . Technical note: single time point dose estimate for exponential clearance. Med Phys. 2018;45:2318–2324.
    OpenUrl
  24. 24.↵
    1. Willowson KP,
    2. Eslick E,
    3. Ryu H,
    4. Poon A,
    5. Bernard EJ,
    6. Bailey DL
    . Feasibility and accuracy of single time point imaging for renal dosimetry following 177Lu-DOTATATE (‘Lutate’) therapy. EJNMMI Phys. 2018;5:33.
    OpenUrl
  25. 25.↵
    1. Merrill S,
    2. Horowitz J,
    3. Traino AC,
    4. Chipkin SR,
    5. Hollot CV,
    6. Chait Y
    . Accuracy and optimal timing of activity measurements in estimating the absorbed dose of radioiodine in the treatment of Graves’ disease. Phys Med Biol. 2011;56:557–571.
    OpenUrlPubMed
  26. 26.↵
    1. Melgar Pérez J,
    2. Orellana Salas A,
    3. Santaella Guardiola Y,
    4. Antoranz Callejo JC
    . Improving individualised dosimetry in radioiodine therapy for hyperthyroidism using population biokinetic modelling. Phys Med. 2019;62:33–40.
    OpenUrl
  27. 27.↵
    1. Glatting G,
    2. Kull T,
    3. Blumstein NM,
    4. et al
    . Dosimetry with 188Re-labelled monoclonal anti-CD66 antibodies. A simplified approach based on a single measurement 3 h p.i. Nuklearmedizin. 2006;45:134–138.
    OpenUrlPubMed
  28. 28.↵
    1. Siebinga H,
    2. de Wit-van der Veen BJ,
    3. Stokkel MDM,
    4. Huitema ADR,
    5. Hendrikx J
    . Current use and future potential of (physiologically based) pharmacokinetic modelling of radiopharmaceuticals: a review. Theranostics. 2022;12:7804–7820.
    OpenUrl
  29. 29.↵
    1. Hardiansyah D,
    2. Riana A,
    3. Beer AJ,
    4. Glatting G
    . Single-time-point dosimetry using model selection and nonlinear mixed-effects modelling: a proof of concept. EJNMMI Phys. 2023;10:12.
    OpenUrl
  30. 30.↵
    1. Harzmann S,
    2. Garcia-Aviles CA,
    3. Nemer U,
    4. Ruf J,
    5. Meyer P,
    6. Mix M
    . Comparison of dead time behaviour of two SPECT/CT systems for Lu-177. Nuklearmedizin. 2016;55:A68.
    OpenUrl
  31. 31.↵
    1. Kadrmas DJ,
    2. Frey EC,
    3. Karimi SS,
    4. Tsui BM
    . Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction. Phys Med Biol. 1998;43:857–873.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Hofheinz F,
    2. Pötzsch C,
    3. Oehme L,
    4. et al
    . Automatic volume delineation in oncological PET. Evaluation of a dedicated software tool and comparison with manual delineation in clinical data sets. Nuklearmedizin. 2012;51:9–16.
    OpenUrlPubMed
  33. 33.↵
    1. Kondev FG
    . Nuclear data sheets for A=177. Nucl Data Sheets. 2019;159:1–412.
    OpenUrl
  34. 34.↵
    1. Burnham KP,
    2. Anderson DR
    . Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. Springer-Verlag; 2002:1–48.
  35. 35.↵
    1. Kletting P,
    2. Schimmel S,
    3. Kestler HA,
    4. et al
    . Molecular radiotherapy: the NUKFIT software for calculating the time-integrated activity coefficient. Med Phys. 2013;40:102504.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Strand SE,
    2. Zanzonico P,
    3. Johnson TK
    . Pharmacokinetic modeling. Med Phys. 1993;20:515–527.
    OpenUrlCrossRefPubMed
  37. 37.↵
    1. Ivashchenko OV,
    2. O’Doherty J,
    3. Hardiansyah D,
    4. et al
    . Time-activity data fitting in molecular radiotherapy: methodology and pitfalls. Phys Med. 2024:117:103192.
  38. 38.↵
    1. Glatting G,
    2. Reske SN
    . Treatment of radioactive decay in pharmacokinetic modeling: influence on parameter estimation in cardiac 13N-PET. Med Phys. 1999;26:616–621.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Hardiansyah D,
    2. Riana A,
    3. Kletting P,
    4. et al
    . A population-based method to determine the time-integrated activity in molecular radiotherapy. EJNMMI Phys. 2021;8:82.
    OpenUrl
  40. 40.↵
    1. Williams LE,
    2. Odom-Maryon TL,
    3. Liu A,
    4. et al
    . On the correction for radioactive decay in pharmacokinetic modeling. Med Phys. 1995;22:1619–1626.
    OpenUrlPubMed
  41. 41.↵
    1. Bonate PL
    . Pharmacokinetic-Pharmacodynamic Modeling and Simulation. Springer; 2011:205–265.
  42. 42.↵
    1. Mould DR,
    2. Upton RN
    . Basic concepts in population modeling, simulation, and model-based drug development: part 2—introduction to pharmacokinetic modeling methods. CPT Pharmacometrics Syst Pharmacol. 2013;2:e38.
    OpenUrl
  43. 43.↵
    1. Hardiansyah D,
    2. Riana A,
    3. Eiber M,
    4. Beer A,
    5. Glatting G
    . Population-based model selection for an accurate estimation of time-integrated activity using non-linear mixed-effects modelling. Z Med Phys. February 20, 2023 [Epub ahead of print].
  44. 44.↵
    1. Anscombe FJ
    . Graphs in statistical analysis. Am Stat. 1973;27:17–21.
    OpenUrlCrossRef
  45. 45.↵
    1. Jiménez-Franco LD,
    2. Kletting P,
    3. Beer AJ,
    4. Glatting G
    . Treatment planning algorithm for peptide receptor radionuclide therapy considering multiple tumor lesions and organs at risk. Med Phys. 2018;45:3516–3523.
    OpenUrl
  46. 46.↵
    1. Kletting P,
    2. Kull T,
    3. Maass C,
    4. et al
    . Optimized peptide amount and activity for 90Y-labeled DOTATATE therapy. J Nucl Med. 2016;57:503–508.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Kletting P,
    2. Muller B,
    3. Erentok B,
    4. et al
    . Differences in predicted and actually absorbed doses in peptide receptor radionuclide therapy. Med Phys. 2012;39:5708–5717.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Kletting P,
    2. Glatting G
    . Model selection for time-activity curves: the corrected Akaike information criterion and the F-test. Z Med Phys. 2009;19:200–206.
    OpenUrlPubMed
  49. 49.↵
    1. Brosch-Lenz J,
    2. Delker A,
    3. Völter F,
    4. et al
    . Toward single-time-point image-based dosimetry of 177Lu-PSMA-617 therapy. J Nucl Med. 2023;64:767–774.
    OpenUrlAbstract/FREE Full Text
  50. 50.
    1. Sandström M,
    2. Freedman N,
    3. Fröss-Baron K,
    4. Kahn T,
    5. Sundin A
    . Kidney dosimetry in 777 patients during 177Lu-DOTATATE therapy: aspects on extrapolations and measurement time points. EJNMMI Phys. 2020;7:73.
    OpenUrl
  51. 51.↵
    1. Peterson AB,
    2. Mirando DM,
    3. Dewaraja YK
    . Accuracy and uncertainty analysis of reduced time point imaging effect on time-integrated activity for 177Lu-DOTATATE PRRT in patients and clinically realistic simulations. EJNMMI Res. 2023;13:57.
    OpenUrl
  52. 52.↵
    1. Stabin MG,
    2. Sparks RB,
    3. Crowe E
    . OLINDA/EXM: the second-generation personal computer software for internal dose assessment in nuclear medicine. J Nucl Med. 2005;46:1023–1027.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Sandström M,
    2. Garske-Román U,
    3. Johansson S,
    4. Granberg D,
    5. Sundin A,
    6. Freedman N
    . Kidney dosimetry during 177Lu-DOTATATE therapy in patients with neuroendocrine tumors: aspects on calculation and tolerance. Acta Oncol. 2018;57:516–521.
    OpenUrl
  • Received for publication July 5, 2023.
  • Revision received January 23, 2024.
PreviousNext
Back to top

In this issue

Journal of Nuclear Medicine: 65 (4)
Journal of Nuclear Medicine
Vol. 65, Issue 4
April 1, 2024
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Complete Issue (PDF)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Journal of Nuclear Medicine.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Single-Time-Point Renal Dosimetry Using Nonlinear Mixed-Effects Modeling and Population-Based Model Selection in [177Lu]Lu-PSMA-617 Therapy
(Your Name) has sent you a message from Journal of Nuclear Medicine
(Your Name) thought you would like to see the Journal of Nuclear Medicine web site.
Citation Tools
Single-Time-Point Renal Dosimetry Using Nonlinear Mixed-Effects Modeling and Population-Based Model Selection in [177Lu]Lu-PSMA-617 Therapy
Deni Hardiansyah, Elham Yousefzadeh-Nowshahr, Felix Kind, Ambros J. Beer, Juri Ruf, Gerhard Glatting, Michael Mix
Journal of Nuclear Medicine Apr 2024, 65 (4) 566-572; DOI: 10.2967/jnumed.123.266268

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Single-Time-Point Renal Dosimetry Using Nonlinear Mixed-Effects Modeling and Population-Based Model Selection in [177Lu]Lu-PSMA-617 Therapy
Deni Hardiansyah, Elham Yousefzadeh-Nowshahr, Felix Kind, Ambros J. Beer, Juri Ruf, Gerhard Glatting, Michael Mix
Journal of Nuclear Medicine Apr 2024, 65 (4) 566-572; DOI: 10.2967/jnumed.123.266268
Twitter logo Facebook logo LinkedIn logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One
Bookmark this article

Jump to section

  • Article
    • Visual Abstract
    • Abstract
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • CONCLUSION
    • DISCLOSURE
    • Footnotes
    • REFERENCES
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF

Related Articles

  • PubMed
  • Google Scholar

Cited By...

  • An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 4: Impact of Fitting Functions in Estimated Absorbed Doses
  • Google Scholar

More in this TOC Section

  • First-in-Human Study of 18F-Labeled PET Tracer for Glutamate AMPA Receptor [18F]K-40: A Derivative of [11C]K-2
  • Detection of HER2-Low Lesions Using HER2-Targeted PET Imaging in Patients with Metastatic Breast Cancer: A Paired HER2 PET and Tumor Biopsy Analysis
  • [11C]Carfentanil PET Whole-Body Imaging of μ-Opioid Receptors: A First in-Human Study
Show more Clinical Investigation

Similar Articles

Keywords

  • STP dosimetry
  • Akaike weight
  • model selection
  • NLME modeling
  • PSMA
SNMMI

© 2025 SNMMI

Powered by HighWire