Patient-specific dosimetry in peptide receptor radionuclide therapy: a clinical review

Kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy based on multiple small VOIs

PURPOSE

The aim was to investigate the use of multiple small VOIs for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

METHOD

The study was based on patient and simulated SPECT images in anthropomorphic geometries. Images were reconstructed using two reconstruction programs (local LundaDose and commercial Hermia) using OS-EM with and without resolution recovery (RR). Five small VOIs were placed to determine the average activity concentration (AC) in each kidney. The study consisted of three steps: (i) determination of the number of iterations for AC convergence based on simulated images; (ii) determination of recovery-coefficients (RCs) for 2 mL VOIs using a separate set of simulated images; (iii) assessment of operator variability in AC estimates for simulated and patient images. Five operators placed the VOIs, using for guidance: a) SPECT/CT with RR, b) SPECT/CT without RR, and c) CT only. For simulated images, time-integrated ACs (TIACs) were evaluated. For patient images, estimated ACs were compared with results of a previous method based on whole-kidney VOIs.

RESULTS

Eight iterations and ten subsets were sufficient for both programs and reconstruction settings. Mean RCs (mean ± SD) with RR were 1.03 ± 0.02 (LundaDose) and 1.10 ± 0.03 (Hermia), and without RR 0.91 ± 0.03 (LundaDose) and 0.94 ± 0.03 (Hermia). Most stable and accurate estimates of the AC were obtained using five 2-mL VOIs guided by SPECT/CT with RR, applying them to images without RR, and including an explicit RC for recovery correction.

CONCLUSION

The small VOI method based on five 2-mL VOIs was found efficient and sufficiently accurate for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

Prediction of [177Lu]Lu-DOTA-TATE therapy response using the absorbed dose estimated from [177Lu]Lu-DOTA-TATE SPECT/CT in patients with metastatic neuroendocrine tumour

BACKGROUND

Peptide receptor radionuclide therapy (PRRT) with [177Lu]Lu-DOTA-TATE has shown efficacy in patients with metastatic neuroendocrine tumours (NETs). Personalised dosimetry is crucial to optimise treatment outcomes and minimise adverse events. In this study, we investigated the correlation between the tumour-absorbed dose (TAD) estimated from [177Lu]Lu-DOTA-TATE SPECT/CT and the therapeutic response.

METHOD

A retrospective analysis was conducted on patients with advanced well-differentiated NETs grades 1-3 who underwent PRRT and exhibited greater uptake than liver on pre-therapeutic [68Ga]Ga-DOTA-TOC PET/CT. Target lesions were selected based on the RECIST 1.1 and PERCIST 1.0 criteria using [177Lu]Lu-DOTA-TATE SPECT/CT and pre-therapeutic contrast-enhanced CT scans. For anatomical image analysis, the sum of the longest diameter (SLD) of the target lesions was measured using the RECIST 1.1 criteria for patient-based analysis and the longest diameter (LD) of the target lesion using the RECIST-L criteria for lesion-based analysis. Standardised uptake values (SUVs) were measured on SPECT/CT images, and TADs were calculated based on the SUVs. Dosimetry was performed using a single SPECT/CT imaging time point at day 4-5 post-therapy. Statistical analyses were conducted to investigate correlations and determine the target lesion responses.

RESULTS

Twenty patients with primary tumour sites and hepatic metastases were included. Fifty-five target lesions, predominantly located in the pancreas and liver, were analysed. The cumulative TAD (lesion-based analysis: r = 0.299-0.301, p = 0.025-0.027), but not the cycle 1 SUV (lesion-based analysis: r = 0.198-0.206, p = 0.131-0.147) or cycle 1 TAD (lesion-based analysis: r = 0.209-0.217, p = 0.112-0.126), exhibited a significant correlation with the change in LD of the target lesion. Binary logistic regression analysis identified the significance of the cumulative TAD in predicting disease control according to the RECIST-L criteria (odds ratio = 1.031-1.051, p = 0.024-0.026).

CONCLUSIONS

The cumulative TAD estimated from [177Lu]Lu-DOTA-TATE SPECT/CT revealed a significant correlation with change in LD, which was significantly higher for the cumulative TAD than for the cycle 1 SUV or TAD. A higher cumulative TAD was associated with disease control in the target lesion. However, considering the limitations inherent to a confined sample size, careful interpretation of these findings is required. Estimation of the cumulative TAD of [177Lu]Lu-DOTA-TATE therapy could guide the platform towards personalised therapy.

Dosimetry of 177Lu-DOTATOC first circle treatment in patients with advanced metastatic neuroendocrine tumors: A pilot study in China

First cycle dosimetry calculation of ¹⁷⁷Lu-DOTATOC (single activity:1.59-3.49 GBq) was carried out in eight patients with advanced neuroendocrine tumors (NETs) who underwent whole-body planar (0.5, 24, 48, 72 h) and SPECT/CT scans (24 h). Focal uptake of ¹⁷⁷Lu-DOTATOC was found in primary and metastatic tumors. Organs with the highest absorbed doses per injected activity were tumors (1.293 ± 0.862 mGy/MBq) and spleen (0.461 ± 0.408 mGy/MBq), while low absorbed doses were observed in kidneys (0.384 ± 0.112 mGy/MBq) and bone marrow (0.0297 ± 0.0123 mGy/MBq). ¹⁷⁷Lu-DOTATOC is safe, well-tolerated and appropriate in Chinese NETs patients for PRRT.

Modelling of energy-dependent spectral resolution for SPECT Monte Carlo simulations using SIMIND

PURPOSE

Monte Carlo (MC) modelling techniques have been used extensively in Nuclear Medicine (NM). The theoretical energy resolution relationship ( ∝ 1 / E ), does not accurately predict the gamma camera detector response across all energies. This study aimed to validate the accuracy of an energy resolution model for the SIMIND MC simulation code emulating the Siemens Symbia T16 dual-head gamma camera.

METHODS

Measured intrinsic energy resolution data (full width half maximum (FWHM) values), for Ba-133, Lu-177, Am-241, Ga-67, Tc-99m, I-123, I-131 and F-18 sources in air, were used to create a fitted model of the energy response of the gamma camera. Both the fitted and theoretical models were used to simulate intrinsic and extrinsic energy spectra using three different scenarios (source in air; source in simple scatter phantom and a clinical voxel-based digital patient phantom).

RESULTS

The results showed the theoretical model underestimated the FWHM values at energies above 160.0 keV up to 23.5 keV. In contrast, the fitted model better predicted the measured FWHM values with differences less than 3.3 keV. The I-131 in-scatter energy spectrum simulated with the fitted model better matched the measured energy spectrum. Higher energy photopeaks, (I-123: 528.9 keV and I-131: 636.9 keV) simulated with the fitted model, more accurately resembled the measured photopeaks. The voxel-based digital patient phantom energy spectra, simulated with the fitted and theoretical models, showed the potential impact of an incorrect energy resolution model when simulating isotopes with multiple photopeaks.

CONCLUSION

Modelling of energy resolution with the proposed fitted model enables the SIMIND user to accurately simulate NM images. A great improvement was seen for high-energy photon emitting isotopes (e.g. I-131), as well as isotopes with multiple photopeaks (e.g. Lu-177, I-131 and Ga-67) in comparison to the theoretical model. This will result in accurate evaluation of radioactivity quantification, which is vital for dosimetric purposes.

Auger electrons for cancer therapy - a review

Background

Auger electrons (AEs) are very low energy electrons that are emitted by radionuclides that decay by electron capture (e.g. ¹¹¹In, ⁶⁷Ga, ^99mTc, ^195mPt, ¹²⁵I and ¹²³I). This energy is deposited over nanometre-micrometre distances, resulting in high linear energy transfer (LET) that is potent for causing lethal damage in cancer cells. Thus, AE-emitting radiotherapeutic agents have great potential for treatment of cancer. In this review, we describe the radiobiological properties of AEs, their radiation dosimetry, radiolabelling methods, and preclinical and clinical studies that have been performed to investigate AEs for cancer treatment.

Results

AEs are most lethal to cancer cells when emitted near the cell nucleus and especially when incorporated into DNA (e.g. ¹²⁵I-IUdR). AEs cause DNA damage both directly and indirectly via water radiolysis. AEs can also kill targeted cancer cells by damaging the cell membrane, and kill non-targeted cells through a cross-dose or bystander effect. The radiation dosimetry of AEs considers both organ doses and cellular doses. The Medical Internal Radiation Dose (MIRD) schema may be applied. Radiolabelling methods for complexing AE-emitters to biomolecules (antibodies and peptides) and nanoparticles include radioiodination (¹²⁵I and ¹²³I) or radiometal chelation (¹¹¹In, ⁶⁷Ga, ^99mTc). Cancer cells exposed in vitro to AE-emitting radiotherapeutic agents exhibit decreased clonogenic survival correlated at least in part with unrepaired DNA double-strand breaks (DSBs) detected by immunofluorescence for γH2AX, and chromosomal aberrations. Preclinical studies of AE-emitting radiotherapeutic agents have shown strong tumour growth inhibition in vivo in tumour xenograft mouse models. Minimal normal tissue toxicity was found due to the restricted toxicity of AEs mostly on tumour cells targeted by the radiotherapeutic agents. Clinical studies of AEs for cancer treatment have been limited but some encouraging results were obtained in early studies using ¹¹¹In-DTPA-octreotide and ¹²⁵I-IUdR, in which tumour remissions were achieved in several patients at administered amounts that caused low normal tissue toxicity, as well as promising improvements in the survival of glioblastoma patients with ¹²⁵I-mAb 425, with minimal normal tissue toxicity.

Conclusions

Proof-of-principle for AE radiotherapy of cancer has been shown preclinically, and clinically in a limited number of studies. The recent introduction of many biologically-targeted therapies for cancer creates new opportunities to design novel AE-emitting agents for cancer treatment. Pierre Auger did not conceive of the application of AEs for targeted cancer treatment, but this is a tremendously exciting future that we and many other scientists in this field envision.

Tumor Drug Penetration Measurements Could Be the Neglected Piece of the Personalized Cancer Treatment Puzzle

Precision medicine aims to use patient genomic, epigenomic, specific drug dose, and other data to define disease patterns that may potentially lead to an improved treatment outcome. Personalized dosing regimens based on tumor drug penetration can play a critical role in this approach. State-of-the-art techniques to measure tumor drug penetration focus on systemic exposure, tissue penetration, cellular or molecular engagement, and expression of pharmacological activity. Using in silico methods, this information can be integrated to bridge the gap between the therapeutic regimen and the pharmacological link with clinical outcome. These methodologies are described, and challenges ahead are discussed. Supported by many examples, this review shows how the combination of these techniques provides enhanced patient-specific information on drug accessibility at the tumor tissue level, target binding, and downstream pharmacology. Our vision of how to apply tumor drug penetration measurements offers a roadmap for the clinical implementation of precision dosing.

Side Effects of CXC-Chemokine Receptor 4-Directed Endoradiotherapy with Pentixather Before Hematopoietic Stem Cell Transplantation

The chemokine receptor CXC-chemokine receptor 4 (CXCR4) is a transmembrane receptor involved in survival, proliferation, and dissemination of different cancers, including hematopoietic malignancies. Relapsed or refractory hematopoietic cancers are frequently resistant to conventional therapy, and novel highly active strategies are urgently needed. CXCR4-directed endoradiotherapy constitutes a highly promising targeted therapeutic concept. Here, we investigated the adverse effects of this novel treatment approach. Methods: Twenty-two patients with heavily pretreated lymphoproliferative or myeloid malignancies were treated with ¹⁷⁷Lu- or ⁹⁰Y-pentixather-a CXCR4-directed therapeutic radioligand-before conventional conditioning therapy followed by autologous or allogeneic hematopoietic stem cell transplantation. Twenty-five CXCR4-directed endoradiotherapies were administered to those patients. Adverse events occurring between endoradiotherapy and the start of conventional conditioning therapy were retrospectively analyzed and graded for the estimation of the safety profile. Results: CXCR4-directed endoradiotherapy with pentixather showed a favorable toxicity profile. As expected, the hematopoietic system was most affected, with all subjects developing cytopenias. Except for 1 acute kidney failure, grade 3, due to tumor lysis syndrome, overall nephro- and hepatotoxicity was low. Other higher-grade adverse events were either transient and resolved or easily manageable. Conclusion: Therapy with radiolabeled pentixather appears to be well tolerated and easily applicable when preceding conventional conditioning regimens for hematopoietic stem cell transplantation.

Radionuclide therapy: current status and prospects for internal dosimetry in individualized therapeutic planning

The efficacy and toxicity of radionuclide therapy are believed to be directly related to the radiation doses received by target tissues; however, nuclear medicine therapy continues to be based primarily on the administration of empirical activities to patients and less frequently on the use of internal dosimetry for individual therapeutic planning. This review aimed to critically describe the techniques and clinical evidence of dosimetry as a tool for therapeutic planning and the main limitations to its implementation in clinical practice. The present article is a nonsystematic review of voxel-based dosimetry. Clinical evidence pointing to a correlation between the radiation dose and therapeutic response in various diseases, such as thyroid carcinoma, neuroendocrine tumors and prostate cancer, is reviewed. Its limitations include technical aspects related to image acquisition and processing and the lack of randomized clinical trials demonstrating the impact of dosimetry on patient therapy. A more widespread use of dosimetry in therapeutic planning involves the development of user-friendly dosimetric protocols and confirmation that dose estimation implies good efficacy and low treatment-related toxicity.

90Y-DOTATOC Dosimetry-Based Personalized Peptide Receptor Radionuclide Therapy

Pretherapy PET with ⁸⁶Y-DOTATOC is considered the ideal dosimetry protocol for ⁹⁰Y-DOTATOC therapy; however, its cost, limited availability, and need for infusion of amino acids to mimic the therapy administration limit its use in the clinical setting. The goal of this study was to develop a dosimetric method for ⁹⁰Y-DOTATOC using ⁹⁰Y-DOTATOC PET/CT and bremsstrahlung SPECT/CT and to determine whether dosimetry-based administered activities differ significantly from standard administered activities. Methods: This was a prospective phase 2 trial of ⁹⁰Y-DOTATOC therapy in patients with somatostatin receptor-positive tumors. ⁹⁰Y-DOTATOC was given in 3 cycles 6-8 wk apart. In the first cycle of therapy, adults received 4.4 GBq and children received 1.85 GBq/m²; the subsequent administered activities were adjusted according to the dosimetry of the preceding cycle so as not to exceed a total kidney dose of 23 Gy and bone marrow dose of 2 Gy. The radiation dose to the kidneys was determined from serial imaging sessions consisting of time-of-flight ⁹⁰Y-DOTATOC PET/CT at 5 h after therapy and ⁹⁰Y-DOTATOC bremsstrahlung SPECT/CT at 6, 24, 48, and 72 h. The PET/CT data were used to measure the absolute concentration of ⁹⁰Y-DOTATOC and to calibrate the bremsstrahlung SPECT kidney clearance data. The radiation dose to the kidneys was determined by multiplying the time-integrated activity (from the fitted biexponential curve of renal clearance of ⁹⁰Y-DOTATOC) with the energy emitted per decay, divided by the mass of the kidneys. Results: The radiation dose to the kidneys per cycle of ⁹⁰Y-DOTATOC therapy was highly variable among patients, ranging from 0.32 to 3.0 mGy/MBq. In 17 (85%) of the 20 adult patients who received the second and the third treatment cycles of ⁹⁰Y-DOTATOC, the administered activity was modified by at least 20% from the starting administered activity. Conclusion: Renal dosimetry of ⁹⁰Y-DOTATOC is feasible using ⁹⁰Y-DOTATOC time-of-flight PET/CT and bremsstrahlung SPECT/CT and has a significant impact on the administered activity in treatment cycles.

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.

Parametric Net Influx Rate Images of 68Ga-DOTATOC and 68Ga-DOTATATE: Quantitative Accuracy and Improved Image Contrast

⁶⁸Ga-DOTATOC and ⁶⁸Ga-DOTATATE are radiolabeled somatostatin analogs used for the diagnosis of somatostatin receptor-expressing neuroendocrine tumors (NETs), and SUV measurements are suggested for treatment monitoring. However, changes in net influx rate (K_i) may better reflect treatment effects than those of the SUV, and accordingly there is a need to compute parametric images showing K_i at the voxel level. The aim of this study was to evaluate parametric methods for computation of parametric K_i images by comparison to volume of interest (VOI)-based methods and to assess image contrast in terms of tumor-to-liver ratio. Methods: Ten patients with metastatic NETs underwent a 45-min dynamic PET examination followed by whole-body PET/CT at 1 h after injection of ⁶⁸Ga-DOTATOC and ⁶⁸Ga-DOTATATE on consecutive days. Parametric K_i images were computed using a basis function method (BFM) implementation of the 2-tissue-irreversible-compartment model and the Patlak method using a descending aorta image-derived input function, and mean tumor K_i values were determined for 50% isocontour VOIs and compared with K_i values based on nonlinear regression (NLR) of the whole-VOI time-activity curve. A subsample of healthy liver was delineated in the whole-body and K_i images, and tumor-to-liver ratios were calculated to evaluate image contrast. Correlation (R²) and agreement between VOI-based and parametric K_i values were assessed using regression and Bland-Altman analysis. Results: The R² between NLR-based and parametric image-based (BFM) tumor K_i values was 0.98 (slope, 0.81) and 0.97 (slope, 0.88) for ⁶⁸Ga-DOTATOC and ⁶⁸Ga-DOTATATE, respectively. For Patlak analysis, the R² between NLR-based and parametric-based (Patlak) tumor K_i was 0.95 (slope, 0.71) and 0.92 (slope, 0.74) for ⁶⁸Ga-DOTATOC and ⁶⁸Ga-DOTATATE, respectively. There was no bias between NLR and parametric-based K_i values. Tumor-to-liver contrast was 1.6 and 2.0 times higher in the parametric BFM K_i images and 2.3 and 3.0 times in the Patlak images than in the whole-body images for ⁶⁸Ga-DOTATOC and ⁶⁸Ga-DOTATATE, respectively. Conclusion: A high R² and agreement between NLR- and parametric-based K_i values was found, showing that K_i images are quantitatively accurate. In addition, tumor-to-liver contrast was superior in the parametric K_i images compared with whole-body images for both ⁶⁸Ga-DOTATOC and ⁶⁸Ga DOTATATE.

Web of Science, Scopus, and Google Scholar citation rates: a case study of medical physics and biomedical engineering: what gets cited and what doesn't?

There are often differences in a publication's citation count, depending on the database accessed. Here, aspects of citation counts for medical physics and biomedical engineering papers are studied using papers published in the journal Australasian physical and engineering sciences in medicine. Comparison is made between the Web of Science, Scopus, and Google Scholar. Papers are categorised into subject matter, and citation trends are examined. It is shown that review papers as a group tend to receive more citations on average; however the highest cited individual papers are more likely to be research papers.

The role of patient-based treatment planning in peptide receptor radionuclide therapy

BACKGROUND

Accurate treatment planning is recommended in peptide-receptor radionuclide therapy (PRRT) to minimize the toxicity to organs at risk while maximizing tumor cell sterilization. The aim of this study was to quantify the effect of different degrees of individualization on the prediction accuracy of individual therapeutic biodistributions in patients with neuroendocrine tumors (NETs).

METHODS

A recently developed physiologically based pharmacokinetic (PBPK) model was fitted to the biokinetic data of 15 patients with NETs after pre-therapeutic injection of (111)In-DTPAOC. Mathematical phantom patients (MPP) were defined using the assumed true (true MPP), mean (MPP 1A) and median (MPP 1B) parameter values of the patient group. Alterations of the degree of individualization were introduced to both mean and median patients by including patient-specific information as a priori knowledge: physical parameters and hematocrit (MPP 2A/2B). Successively, measurable individual biokinetic parameters were added: tumor volume V tu (MPP 3A/3B), glomerular filtration rate GFR (MPP 4A/4B), and tumor perfusion f tu (MPP 5A/5B). Furthermore, parameters of MPP 5A/5B and a simulated (68)Ga-DOTATATE PET measurement 60 min p.i. were used together with the population values used as Bayesian parameters (MPP 6A/6B). Therapeutic biodistributions were simulated assuming an infusion of (90)Y-DOTATATE (3.3 GBq) over 30 min to all MPPs. Time-integrated activity coefficients were predicted for all MPPs and compared to the true MPPs for each patient in tumor, kidneys, spleen, liver, remainder, and whole body to obtain the relative differences RD.

RESULTS

The large RD values of MPP 1A [RDtumor = (625 ± 1266)%, RDkidneys = (11 ± 38)%], and MPP 1B [RDtumor = (197 ± 505)%, RDkidneys = (11 ± 39)%] demonstrate that individual treatment planning is needed due to large physiological differences between patients. Although addition of individual patient parameters reduced the deviations considerably [MPP 5A: RDtumor = (-2 ± 27)% and RDkidneys = (16 ± 43)%; MPP 5B: RDtumor = (2 ± 28)% and RDkidneys = (7 ± 40)%] errors were still large. For the kidneys, prediction accuracy was considerably improved by including the PET measurement [MPP 6A/MPP 6B: RDtumor = (-2 ± 22)% and RDkidneys = (-0.1 ± 0.5)%].

CONCLUSION

Individualized treatment planning is needed in the investigated patient group. The use of a PBPK model and the inclusion of patient specific data, e.g., weight, tumor volume, and glomerular filtration rate, do not suffice to predict the therapeutic biodistribution. Integrating all available a priori information in the PBPK model and using additionally PET data measured at one time point for tumor, kidneys, spleen, and liver could possibly be sufficient to perform an individualized treatment planning.