An automated voxelized dosimetry tool for radionuclide therapy based on serial quantitative SPECT/CT imaging

Clinical Trial Protocol for VIOLET: A Single-Center, Phase I/II Trial Evaluation of Radioligand Treatment in Patients with Metastatic Castration-Resistant Prostate Cancer with [161Tb]Tb-PSMA-I&T

[177Lu]Lu-PSMA is an effective class of therapy for patients with metastatic castration-resistant prostate cancer (mCRPC); however, progression is inevitable. The limited durability of response may be partially explained by the presence of micrometastatic deposits, which are energy-sheltered and receive low absorbed radiation with 177Lu due to the approximately 0.7-mm mean pathlength. 161Tb has abundant emission of Auger and conversion electrons that deposit a higher concentration of radiation over a shorter path, particularly to single tumor cells and micrometastases. 161Tb has shown in vitro and in vivo efficacy superior to that of 177Lu. We aim to demonstrate that [161Tb]Tb-PSMA-I&T will deliver effective radiation to sites of metastatic prostate cancer with an acceptable safety profile. Methods: This single-center, single-arm, phase I/II trial will recruit 30 patients with mCRPC. Key eligibility criteria include a diagnosis of mCRPC with progression after at least one line of taxane chemotherapy (unless medically unsuitable) and androgen receptor pathway inhibitor; prostate-specific membrane antigen-positive disease on [68Ga]Ga-PSMA-11 or [18F]DCFPyL PET/CT (SUVmax ≥ 20); no sites of discordance on [18F]FDG PET/CT; adequate bone marrow, hepatic, and renal function; an Eastern Cooperative Oncology Group performance status of no more than 2, and no prior treatment with another radioisotope. The dose escalation is a 3 + 3 design to establish the safety of 3 prespecified activities of [161Tb]Tb-PSMA-I&T (4.4, 5.5, and 7.4 GBq). The maximum tolerated dose will be defined as the highest activity level at which a dose-limiting toxicity occurs in fewer than 2 of 6 participants. The dose expansion will include 24 participants at the maximum tolerated dose. Up to 6 cycles of [161Tb]Tb-PSMA-I&T will be administered intravenously every 6 wk, with each subsequent activity reduced by 0.4 GBq. The coprimary objectives are to establish the maximum tolerated dose and safety profile (Common Terminology Criteria for Adverse Events version 5.0) of [161Tb]Tb-PSMA-I&T. Secondary objectives include measuring absorbed radiation dose (Gy), evaluating antitumor activity (prostate-specific antigen 50% response rate, radiographic and prostate-specific antigen progression-free survival, overall survival, objective response rate), and evaluating pain (Brief Pain Inventory-Short Form) and health-related quality of life (Functional Assessment of Cancer Therapy-Prostate and Functional Assessment of Cancer Therapy-Radionuclide Therapy). Conclusion: Enrollment was completed in February 2024. Patients are still receiving [161Tb]Tb-PSMA-I&T.

Gastrin Releasing Peptide Receptors-targeted PET Diagnostics and Radionuclide Therapy for Prostate Cancer Management: Preclinical and Clinical Developments of the Past 5 Years

Each tumor has its own distinctive molecular identity. Treatment, therefore, should be tailored to this unique cancer phenotype. Theragnostics uses the same compound for targeted imaging and treatment, radiolabeled to an appropriate radionuclide, respectively. Gastrin-releasing peptide receptors (GRPRs) are overexpressed in prostate cancer, and radiolabeled GRPR antagonists have shown high diagnostic performance at staging and biochemical recurrence. Several GRPR-targeting theragnostic compounds have been developed preclinically. Their translation into clinics is underway with 4 clinical trials recruiting participants with GRPR-expressing tumors.

Impact of different models based on blood samples and images for bone marrow dosimetry after 177Lu-labeled somatostatin-receptor therapy

BACKGROUND

Peptide receptor radionuclide therapy with 177Lu-DOTATATE is a recognized option for treating neuroendocrine tumors and has few toxicities, except for the kidneys and bone marrow. The bone marrow dose is generally derived from a SPECT/CT image-based method with four timepoints or from a blood-based method with up to 9 timepoints, but there is still no reference method. This retrospective single-center study on the same cohort of patients compared the calculated bone marrow dose administered with both methods using mono, bi- or tri-exponential models. For the image-based method, the dose was estimated using Planetdose© software. Pearson correlation coefficients were calculated. We also studied the impact of late timepoints for both methods.

RESULTS

The bone marrow dose was calculated for 131 treatments with the blood-based method and for 17 with the image-based method. In the former, the median absorbed dose was 15.3, 20.5 and 28.3 mGy/GBq with the mono-, bi- and tri-exponential model, respectively. With the image-based method, the median absorbed dose was 63.9, 41.9 and 60.8 with the mono-, bi- and tri-exponential model, respectively. Blood samples after 24h post-injection did not evidence any change in the absorbed bone marrow dose with the bi-exponential model. On the contrary, the 6-day post-injection timepoint was more informative with the image-based model.

CONCLUSION

This study confirms that the estimated bone marrow dose is significantly lower with the blood-based method than with the image-based method. The blood-based method with a bi-exponential model proved particularly useful, without the need for blood samples after 24h post-injection. Nevertheless, this blood-based method is based on an assumption that needs to be more validated. The important difference between the two methods does not allow to determine the optimal one to estimate the true absorbed dose and further studies are necessary to compare with biological effects.

Administering [177Lu]Lu-PSMA-617 Prior to Radical Prostatectomy in Men with High-risk Localised Prostate Cancer (LuTectomy): A Single-centre, Single-arm, Phase 1/2 Study

BACKGROUND

High-risk localised prostate cancer (HRCaP) has high rates of biochemical recurrence; [177Lu]Lu-PSMA-617 is effective in men with advanced prostate cancer.

OBJECTIVE

To investigate the dosimetry, safety, and efficacy of upfront [177Lu]Lu-PSMA-617 in men with HRCaP prior to robotic radical prostatectomy (RP).

DESIGN, SETTING, AND PARTICIPANTS

In this single-arm, phase I/II trial, we recruited men with HRCaP (any of prostate-specific antigen [PSA] >20 ng/ml, International Society of Urological Pathology (ISUP) grade group [GG] 3-5, and ≥cT2c), with high tumour uptake on [68Ga]Ga-PSMA-11 positron emission tomography/computed tomography (PSMA PET/CT), and scheduled for RP.

INTERVENTION

Cohort A (n = 10) received one cycle and cohort B (n = 10) received two cycles of [177Lu]Lu-PSMA-617 (5 GBq) followed by surgery 6 weeks later.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

The primary endpoint was tumour radiation absorbed dose. Adverse events (AEs; Common Terminology Criteria for Adverse Events (CTCAE) version 5.0), surgical safety (Clavien-Dindo), imaging, and biochemical responses were evaluated (ClinicalTrials.gov: NCT04430192).

RESULTS AND LIMITATIONS

Between May 29, 2020 and April 28, 2022, 20 patients were enrolled. The median PSA was 18 ng/ml (interquartile range [IQR] 11-35), Eighteen (90%) had GG ≥3, and six (30%) had N1 disease. The median (IQR) highest tumour radiation absorbed dose after cycle 1 for all lesions was 35.5 Gy (19.5-50.1), with 19.6 Gy (11.3-48.4) delivered to the prostate. Five patients received radiation to lymph nodes. Nine (45%) patients achieved >50% PSA decline. The most common AEs related to [177Lu]Lu-PSMA-617 were grade 1 fatigue in eight (40%), nausea in seven (35%), dry mouth in six (30%), and thrombocytopenia in four (20%) patients. No grade 3/4 toxicities or Clavien 3-5 complications occurred. Limitations include small a sample size.

CONCLUSIONS

In men with HRCaP and high prostate-specific membrane antigen (PSMA) expression, [177Lu]Lu-PSMA-617 delivered high levels of targeted radiation doses with few toxicities and without compromising surgical safety. Further studies of [177Lu]Lu-PSMA-617 in this population are worthwhile to determine whether meaningful long-term oncological benefits can be demonstrated.

PATIENT SUMMARY

In this study, we demonstrate that up to two cycles of [177Lu]Lu-PSMA-617 given prior to radical prostatectomy in patients with high-risk localised prostate cancer are safe and deliver targeted doses of radiation to tumour-affected tissues. It is tolerated well with minimal treatment-related adverse events, and surgery is safe with a low rate of complications. Activity measured through PSA reduction, repeat PSMA PET/CT, and histological response is promising.

Quantitative calibration of Tb-161 SPECT/CT in view of personalised dosimetry assessment studies

BACKGROUND

Terbium-161 (161Tb)-based radionuclide therapy poses an alternative to current Lutetium-177 (177Lu) approaches with the additional benefit of secondary Auger and conversion electron emissions capable of delivering high doses of localised damage to micro-metastases including single cells. Quantitative single-photon emission computed tomography, paired with computed tomography (SPECT/CT), enables quantitative measurement from post-therapy imaging. In view of dosimetry extrapolations, a Tb-161 sensitivity SPECT/CT camera calibration was performed using a method previously validated for 177Lu.

METHODS

Serial imaging of a NEMA/IEC body phantom with Tb-161 was performed on SPECT/CT with low-energy high-resolution collimators employing a photopeak of 75 keV with a 20% width. Quantitative stability and recovery coefficients were investigated over a sequence of 19 scans with buffered 161Tb solution at total phantom activity ranging from 70 to 4990 MBq.

RESULTS

Sphere recovery coefficients were 0.60 ± 0.05, 0.52 ± 0.07, 0.45 ± 0.07, 0.39 ± 0.07, 0.28 ± 0.08, and 0.20 ± 0.08 for spheres 37, 28, 22, 17, 13, and 10mm, respectively, when considered across all activity and scan durations with dual-energy window scatter correction. Whole-field reconstructed sensitivity was calculated as 1.42E-5 counts per decay. Qualitatively, images exhibited no visual artefacts and were comparable to 177Lu SPECT/CT.

CONCLUSIONS

Quantitative SPECT/CT of 161Tb is feasible over a range of activities enabling dosimetry analogous to 177Lu whilst also producing suitable imaging for clinical review. This has been incorporated into a prospective trial of 161Tb-PSMA for men with metastatic prostate cancer.

The relationship between tumour dosimetry, response, and overall survival in patients with unresectable Neuroendocrine Neoplasms (NEN) treated with 177Lu DOTATATE (LuTate)

Peptide Receptor Radionuclide Therapy (PRRT) delivers targeted radiation to Somatostatin Receptor (SSR) expressing Neuroendocrine Neoplasms (NEN). We sought to assess the predictive and prognostic implications of tumour dosimetry with respect to response by 68 Ga DOTATATE (GaTate) PET/CT molecular imaging tumour volume of SSR (MITVSSR) change and RECIST 1.1, and overall survival (OS).

METHODS

Patients with gastro-entero-pancreatic (GEP) NEN who received LuTate followed by quantitative SPECT/CT (Q-SPECT/CT) the next day (Jul 2010 to Jan 2019) were retrospectively reviewed. Single time-point (STP) lesional dosimetry was performed for each cycle using population-based pharmacokinetic modelling. MITVSSR and RECIST 1.1 were measured at 3-months post PRRT.

RESULTS

Median of 4 PRRT cycles were administered to 90 patients (range 2-5 cycles; mean 27.4 GBq cumulative activity; mean 7.6 GBq per cycle). 68% received at least one cycle with radiosensitising chemotherapy (RSC). RECIST 1.1 partial response was 24%, with 70% stable and 7% progressive disease. Cycle 1 radiation dose in measurable lesions was associated with local response (odds ratio 1.5 per 50 Gy [95% CI: 1.1-2.0], p = 0.002) when adjusted by tumour grade and RSC. Median change in MITVSSR was -63% (interquartile range -84 to -29), with no correlation with radiation dose to the most avid lesion on univariable or multivariant analyses (5.6 per 10 Gy [95% CI: -1.6, 12.8], p = 0.133). OS at 5-years was 68% (95% CI: 56-78%). Neither baseline MITVSSR (hazard ratio 1.1 [95% CI: 1.0, 1.2], p = 0.128) nor change in baseline MITVSSR (hazard ratio 1.0 [95% CI: 1.0, 1.1], p = 0.223) were associated with OS when adjusted by tumour grade and RSC but RSC was (95% CI: 0.2, 0.8, p = 0.012).

CONCLUSION

Radiation dose to tumour during PRRT was predictive of radiologic response but not survival. Survival outcomes may relate to other biological factors. There was no evidence that MITVSSR change was associated with OS, but a larger study is needed.

Safety of Lutetium-177 prostate-specific membrane antigen-617 (PSMA-617) radioligand therapy in the setting of severe renal impairment: a case report and literature review

Reported here is a case of rapidly progressive metastatic castration-resistant prostate cancer treated with [177Lu]Lu-PSMA-617 in the setting of severe renal impairment and impending ureteric obstruction. PSMA is expressed on renal tubular cells, raising the possibility of radiation-induced nephrotoxicity, and this level of renal impairment would typically exclude the patient from [177Lu]Lu-PSMA-617 therapy. Multidisciplinary input, individualized dosimetry, and patient-specific dose reduction were used to ensure the cumulative dose to the kidneys remained within acceptable limits. He was initially planned for treatment with six cycles of [177Lu]Lu-PSMA-617. However, he had an excellent response to therapy following four cycles of treatment and the last two cycles were omitted. He has been followed for 1-year posttherapy without evidence of disease recurrence. No acute or chronic nephrotoxicity was observed. This case report highlights the utility of [177Lu]Lu-PSMA-617 therapy in severe renal impairment and provides evidence of relative safety in patients who would otherwise not be considered candidates for therapy.

Safety and Therapeutic Optimization of Lutetium-177 Based Radiopharmaceuticals

Peptide receptor radionuclide therapy (PRRT) using Lutetium-177 (¹⁷⁷Lu) based radiopharmaceuticals has emerged as a therapeutic area in the field of nuclear medicine and oncology, allowing for personalized medicine. Since the first market authorization in 2018 of [¹⁷⁷Lu]Lu-DOTATATE (Lutathera®) targeting somatostatin receptor type 2 in the treatment of gastroenteropancreatic neuroendocrine tumors, intensive research has led to transfer innovative ¹⁷⁷Lu containing pharmaceuticals to the clinic. Recently, a second market authorization in the field was obtained for [¹⁷⁷Lu]Lu-PSMA-617 (Pluvicto®) in the treatment of prostate cancer. The efficacy of ¹⁷⁷Lu radiopharmaceuticals are now quite well-reported and data on the safety and management of patients are needed. This review will focus on several clinically tested and reported tailored approaches to enhance the risk-benefit trade-off of radioligand therapy. The aim is to help clinicians and nuclear medicine staff set up safe and optimized procedures using the approved ¹⁷⁷Lu based radiopharmaceuticals.

Patient-specific dosimetry adapted to variable number of SPECT/CT time-points per cycle for [Formula: see text]Lu-DOTATATE therapy

BACKGROUND

The number of SPECT/CT time-points is important for accurate patient dose estimation in peptide receptor radionuclide therapy. However, it may be limited by the patient's health and logistical reasons. Here,  an image-based dosimetric workflow adapted to the number of SPECT/CT acquisitions available throughout the treatment cycles was proposed, taking into account patient-specific pharmacokinetics and usable in clinic for all organs at risk.

METHODS

Thirteen patients with neuroendocrine tumors were treated with four injections of 7.4 GBq of [Formula: see text]Lu-DOTATATE. Three SPECT/CT images were acquired during the first cycle (1H, 24H and 96H or 144H post-injection) and a single acquisition (24H) for following cycles. Absorbed doses were estimated for kidneys (LK and RK), liver (L), spleen (S), and three surrogates of bone marrow (L2 to L4, L1 to L5 and T9 to L5) that were compared. 3D dose rate distributions were computed with Monte Carlo simulations. Voxel dose rates were averaged at the organ level. The obtained Time Dose-Rate Curves (TDRC) were fitted with a tri-exponential model and time-integrated. This method modeled patient-specific uptake and clearance phases observed at cycle 1. Obtained fitting parameters were reused for the following cycles, scaled to the measure organ dose rate at 24H. An alternative methodology was proposed when some acquisitions were missing based on population average TDRC (named STP-Inter). Seven other patients with three SPECT/CT acquisitions at cycles 1 and 4 were included to estimate the uncertainty of the proposed methods.

RESULTS

Absorbed doses (in Gy) per cycle available were: 3.1 ± 1.1 (LK), 3.4 ± 1.5 (RK), 4.5 ± 2.8 (L), 4.6 ± 1.8 (S), 0.3 ± 0.2 (bone marrow). There was a significant difference between bone marrow surrogates (L2 to L4 and L1 to L5, Wilcoxon's test: p value < 0.05), and while depicting very doses, all three surrogates were significantly different than dose in background (p value < 0.01). At cycle 1, if the acquisition at 24H is missing and approximated, medians of percentages of dose difference (PDD) compared to the initial tri-exponential function were inferior to 3.3% for all organs. For cycles with one acquisition, the median errors were smaller with a late time-point. For STP-Inter, medians of PDD were inferior to 7.7% for all volumes, but it was shown to depend on the homogeneity of TDRC.

CONCLUSION

The proposed workflow allows the estimation of organ doses, including bone marrow, from a variable number of time-points acquisitions for patients treated with [Formula: see text]Lu-DOTATATE.

EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands

The purpose of the EANM Dosimetry Committee is to provide recommendations and guidance to scientists and clinicians on patient-specific dosimetry. Radiopharmaceuticals labelled with lutetium-177 (177Lu) are increasingly used for therapeutic applications, in particular for the treatment of metastatic neuroendocrine tumours using ligands for somatostatin receptors and prostate adenocarcinoma with small-molecule PSMA-targeting ligands. This paper provides an overview of reported dosimetry data for these therapies and summarises current knowledge about radiation-induced side effects on normal tissues and dose-effect relationships for tumours. Dosimetry methods and data are summarised for kidneys, bone marrow, salivary glands, lacrimal glands, pituitary glands, tumours, and the skin in case of radiopharmaceutical extravasation. Where applicable, taking into account the present status of the field and recent evidence in the literature, guidance is provided. The purpose of these recommendations is to encourage the practice of patient-specific dosimetry in therapy with 177Lu-labelled compounds. The proposed methods should be within the scope of centres offering therapy with 177Lu-labelled ligands for somatostatin receptors or small-molecule PSMA.

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.

DblurDoseNet: A deep residual learning network for voxel radionuclide dosimetry compensating for single-photon emission computerized tomography imaging resolution

PURPOSE

Current methods for patient-specific voxel-level dosimetry in radionuclide therapy suffer from a trade-off between accuracy and computational efficiency. Monte Carlo (MC) radiation transport algorithms are considered the gold standard for voxel-level dosimetry but can be computationally expensive, whereas faster dose voxel kernel (DVK) convolution can be suboptimal in the presence of tissue heterogeneities. Furthermore, the accuracies of both these methods are limited by the spatial resolution of the reconstructed emission image. To overcome these limitations, this paper considers a single deep convolutional neural network (CNN) with residual learning (named DblurDoseNet) that learns to produce dose-rate maps while compensating for the limited resolution of SPECT images.

METHODS

We trained our CNN using MC-generated dose-rate maps that directly corresponded to the true activity maps in virtual patient phantoms. Residual learning was applied such that our CNN learned only the difference between the true dose-rate map and DVK dose-rate map with density scaling. Our CNN consists of a 3D depth feature extractor followed by a 2D U-Net, where the input was 11 slices (3.3 cm) of a given Lu-177 SPECT/CT image and density map, and the output was the dose-rate map corresponding to the center slice. The CNN was trained with nine virtual patient phantoms and tested on five different phantoms plus 42 SPECT/CT scans of patients who underwent Lu-177 DOTATATE therapy.

RESULTS

When testing on virtual patient phantoms, the lesion/organ mean dose-rate error and the normalized root mean square error (NRMSE) relative to the ground truth of the CNN method was consistently lower than DVK and MC, when applied to SPECT images. Compared to DVK/MC, the average improvement for the CNN in mean dose-rate error was 55%/53% and 66%/56%; and in NRMSE was 18%/17% and 10%/11% for lesion and kidney regions, respectively. Line profiles and dose-volume histograms demonstrated compensation for SPECT resolution effects in the CNN-generated dose-rate maps. The ensemble noise standard deviation, determined from multiple Poisson realizations, was improved by 21%/27% compared to DVK/MC. In patients, potential improvements from CNN dose-rate maps compared to DVK/MC were illustrated qualitatively, due to the absence of ground truth. The trained residual CNN took about 30 s on a single GPU (Tesla V100) to generate a 512 × 512 × 130 dose-rate map for a patient.

CONCLUSION

The proposed residual CNN, trained using phantoms generated from patient images, has potential for real-time patient-specific dosimetry in clinical treatment planning due to its demonstrated improvement in accuracy, resolution, noise, and speed over the DVK/MC approaches.

Radiation Dosimetry in 177Lu-PSMA-617 Therapy

Radionuclide therapy using the small molecule PSMA bound to the beta-emitting radionuclide, Lutetium-177 (¹⁷⁷Lu-PSMA) has demonstrated efficacy and survival benefit castrate resistant metastatic disease and represents a novel new line of therapy. Whilst dosimetry was critical for early development, it was not incorporated into either the TheraP or VISION randomized studies, highlighting the difficulty of adopting dosimetry in routine clinical practice. Accumulated clinical experience has also shown that the common (and generally low grade) toxicities such as nausea, xerostomia, and cytopenias are not readily predicted on the basis of dosimetry estimates. The majority of dosimetry and clinical literature deals with the radiopharmaceutical ¹⁷⁷Lu-PSMA-617 which displays relatively consistent patterns of retention among normal tissues and high specificity for metastatic prostate cancer phenotypes. Population dosimetry incorporating estimates to the kidneys, salivary glands, and bone marrow have been widely reported the typical range of doses is becoming well established. There is growing interest on tumor dosimetry in ¹⁷⁷Lu-PSMA-617 therapy as an overall modest side-effect profile from primary organ retention has been observed. A focus away from normal organ dosimetry to whole body tumor dosimetry may enable early prediction of treatment failure. Given the safety of ¹⁷⁷Lu-PSMA there is also potential to escalate administered radioactivity to further improve outcomes. Importantly, the variability of uptake between individuals, both to tumor and normal organs, has also been highlighted which provides some rationale for the utility of personalized radiation analysis to optimize treatment based on potential toxicity thresholds or tumor control. Methods to perform dosimetry using serial post treatment imaging may incorporate planar, 3D SPECT, or hybrid datasets. Reliable measurements may be obtained through either method, however, continued developments in computational analysis are better suited to fully 3D imaging; particularly in conjunction with volumetric CT to assist with alignment and contouring. Dose analysis over sequential treatment cycles is vital to understand the radiobiology of these treatments which is unique compared to external beam therapy due to dose rate, fractionation scheme, and potential for intratumoral nonuniformity.

Dosimetry in Clinical Radiopharmaceutical Therapy of Cancer: Practicality Versus Perfection in Current Practice

The use of radiopharmaceutical therapies (RPTs) in the treatment of cancers is growing rapidly, with more agents becoming available for clinical use in last few years and many new RPTs being in development. Dosimetry assessment is critical for personalized RPT, insofar as administered activity should be assessed and optimized in order to maximize tumor-absorbed dose while keeping normal organs within defined safe dosages. However, many current clinical RPTs do not require patient-specific dosimetry based on current Food and Drug Administration-labeled approvals, and overall, dosimetry for RPT in clinical practice and trials is highly varied and underutilized. Several factors impede rigorous use of dosimetry, as compared with the more convenient and less resource-intensive practice of empiric dosing. We review various approaches to applying dosimetry for the assessment of activity in RPT and key clinical trials, the extent of dosimetry use, the relative pros and cons of dosimetry-based versus fixed activity, and practical limiting factors pertaining to current clinical practice.

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.

Managing radioiodine refractory thyroid cancer: the role of dosimetry and redifferentiation on subsequent I-131 therapy

Poor responses to iodine-131 (I-131) therapy can relate to either low iodine uptake and retention in thyroid cancer cells or to increased radioresistance. Both mechanisms are currently termed radioactive iodine (RAI)-refractory (RAI-R) thyroid cancer but the first reflects unsuitability for I-131 therapy that can be evaluated in advance of treatment, whereas the other can only be identified post hoc. Management of both represents a considerable challenge in clinical practice as failure of I-131 therapy, the most effective treatment of metastatic thyroid cancer, is associated with a poor overall prognosis. The development of targeted therapies has shown substantial promise in the treatment of RAI-R thyroid cancer in progressive patients. Recent studies show that selective tyrosine kinase inhibitors (TKIs) targeting B-type rapidly accelerated fibrosarcoma kinase (BRAF) and mitogen-activated protein kinase (MEK) can be used as redifferentiation agents to re-induce RAI uptake, thereby (re)enabling I-131 therapy. The use of dosimetry prior- and post-TKI treatment can assist in quantifying RAI uptake and improve identification of patients that will benefit from I-131 therapy. It also potentially offers the prospect of calculating individualized therapeutic administered activities to enhance efficacy and limit toxicity. In this review, we present an overview of the regulation of RAI uptake and clinically investigated redifferentiation agents, both reimbursed and in experimental setting, that induce renewed RAI uptake. We describe the role of dosimetry in redifferentiation and subsequent I-131 therapy in RAI-R thyroid cancer, explain different dosimetry approaches and discuss limitations and considerations in the field.

Therapeutic Efficacy of a Bivalent Inhibitor of Prostate-Specific Membrane Antigen Labeled with 67Cu

Radionuclide therapy targeting prostate-specific membrane antigen (PSMA) is promising for prostate cancer. We previously reported a ligand, ⁶⁴Cu-CuSarbisPSMA, featuring 2 lysine-ureido-glutamate groups. Here, we report the therapeutic potential of ⁶⁷Cu-CuSarbisPSMA. Methods: Growth of PSMA-positive xenografts was evaluated after treatment with ⁶⁷Cu-CuSarbisPSMA or ¹⁷⁷Lu-LuPSMA imaging and therapy (I&T). Results: At 13 d after injection, tumor growth was similarly inhibited by the 2 tracers in a dose-dependent manner. Survival was comparable after single (30 MBq) or fractionated (2 × 15 MBq, 2 wk apart) administrations. Conclusion: ⁶⁷Cu-CuSarbisPSMA is efficacious in a PSMA-expressing model of prostate cancer.

Comparison of different methods for post-therapeutic dosimetry in [177Lu]Lu-PSMA-617 radioligand therapy

Background

Dosimetry is of high importance for optimization of patient-individual PSMA-targeted radioligand therapy (PSMA-RLT). The aim of our study was to evaluate and compare the feasibility of different approaches of image-based absorbed dose estimation in terms of accuracy and effort in clinical routine.

Methods

Whole-body planar images and SPECT/CT images were acquired from 24 patients and 65 cycles at 24h, 48h, and ≥96h after administration of a mean activity of 6.4 GBq [¹⁷⁷Lu]Lu-PSMA-617 (range 3–10.9 GBq). Dosimetry was performed by use of the following approaches: 2D planar-based dosimetry, 3D SPECT/CT-based dosimetry, and hybrid dosimetry combining 2D and 3D data. Absorbed doses were calculated according to IDAC 2.1 for the kidneys, the liver, the salivary glands, and bone metastases.

Results

Mean absorbed doses estimated by 3D dosimetry (the reference method) were 0.54 ± 0.28 Gy/GBq for the kidneys, 0.10 ± 0.05 Gy/GBq for the liver, 0.81 ± 0.34 Gy/GBq for the parotid gland, 0.72 ± 0.39 Gy/GBq for the submandibular gland, and 1.68 ± 1.32 Gy/GBq for bone metastases. Absorbed doses of normal organs estimated by hybrid dosimetry showed small, non-significant differences (median up to 4.0%) to the results of 3D dosimetry. Using 2D dosimetry, in contrast, significant differences (median up to 10.9%) were observed. Regarding bone metastases, small, but significant differences (median up to 7.0%) of absorbed dose were found for both, 2D dosimetry and hybrid dosimetry. Bland-Altman analysis revealed high agreement between hybrid dosimetry and 3D dosimetry for normal organs and bone metastases, but substantial differences between 2D dosimetry and 3D dosimetry.

Conclusion

Hybrid dosimetry provides high accuracy in estimation of absorbed dose in comparison to 3D dosimetry for all important organs and is therefore feasible for use in individualized PSMA-RLT.

Description of the methodology for dosimetric quantification in treatments with 177Lu-DOTATATE

Implementation of dosimetry calculations in the daily practice of Nuclear Medicine Departments is, at this time, a controversial issue, partly due to the lack of a standardized methodology that is accepted by all interested parties (patients, nuclear medicine physicians and medical physicists). However, since the publication of RD 601/2019 there is a legal obligation to implement it, despite the fact that it is a complex and high resource consumption procedure. The aim of this article is to review the theoretical bases of in vivo dosimetry in treatments with 177Lu-DOTATATE. The exposed methodology is the one proposed by the MIRD Committee (Medical Internal Radiation Dose) of the SNMMI (Society of Nuclear Medicine & Molecular Imaging). According to this method, the absorbed dose is obtained as the product of 2factors: the time-integrated activity of the radiopharmaceutical present in a source region and a geometrical factor S. This approach, which a priori seems simple, in practice requires several SPECT/CT acquisitions, several measurements of the whole body activity and taking several blood samples, as well as hours of image processing and computation. The systematic implementation of these calculations, in all the patients we treat, will allow us to obtain homogeneous data to correlate the absorbed doses in the lesions with the biological effect of the treatment. The final purpose of the dosimetry calculations is to be able to maximize the therapeutic effect in the lesions, controlling the radiotoxicity in the organs at risk.

Dose voxel kernel prediction with neural networks for radiation dose estimation

BACKGROUND

Currently there is an ever increasing interest in Lu-177 targeted radionuclide therapies, which target neuro-endocrine and prostate tumours. For a patient-specific treatment, an individual dosimetry based on SPECT/CT imaging is necessary. The aim of this study is to introduce a dosimetry method, where dose voxel kernels (DVK) are predicted by a neural network.

METHODS

Kidneys are considered one of the most critical organs in any radionuclide therapy. Hence we chose kidneys of 26 patients, who underwent Lu-177-DOTATOC or PSMA therapy, as target organs for our dosimetric method. First of all, density kernels with a size of 9×9×9 voxels were considered, and the corresponding DVKs were calculated by Monte Carlo simulations. These kernels were used to train a neural network (NN), which received a density kernel as input and predicted a DVK at the output. To predict the dose distribution of an entire kidney, the organ had to be partitioned into a large number of density kernels. Afterwards the DVKs were predicted by a trained NN, and employed to reconstruct the whole-organ dose distribution by convolution with the activity distribution. This method was compared to the standard method where the activity distribution is convolved with a DVK based on a homogeneous soft tissue kernel.

RESULTS

The number of training kernels amounted to 52,274 density kernels with corresponding MC-derived DVKs. The results serve as proof of principle of the newly proposed method and showed that the NN approach yielded superior results compared to the standard method with no additional computational effort.

CONCLUSION

The NN approach is an accurate and highly competitive dosimetric method to precisely estimate absorbed radiation dose in critical organs like kidneys in clinical routine. To further improve the results, a larger number of DVKs needs to be computed by Monte Carlo simulations. An extension of the method to other organs is easily conceivable.

Dosimetry for Optimized, Personalized Radiopharmaceutical Therapy

Radiopharmaceutical therapy (RPT) has grown rapidly over the last decade for treatment of numerous cancer types. Dosimetric guidance, as with other radiotherapy modalities, has benefitted patients by reducing the incidence of side effects and improving overall survival in populations treated under this paradigm. Development of tools and techniques for dosimetry-guided therapy is ongoing, with numerous the Food and Drug Administration-cleared products reaching the U.S. market in 2019. Safe use of commercial dosimetry platforms requires a deep understanding of the underlying physical principles and thoroughly vetted input data. Likewise, interpretation of dosimetry results relies on an understanding of radiobiological principles, and the principles of uncertainty propagation. In this article, we review strategies commonly employed for dosimetry-guided RPT - including quantitative imaging, dose calculation methods, and modeling of dose across time-points. Additionally, we review recent literature evidence (2013-2020) demonstrating the efficacy of personalized RPT.

Meeting report from the Prostate Cancer Foundation PSMA theranostics state of the science meeting

INTRODUCTION

The Prostate Cancer Foundation (PCF) convened a PCF prostate-specific membrane antigen (PSMA) Theranostics State of the Science Meeting on 18 November 2019, at Weill Cornell Medicine, New York, NY.

METHODS

The meeting was attended by 22 basic, translational, and clinical researchers from around the globe, with expertise in PSMA biology, development and use of PSMA theranostics agents, and clinical trials. The goal of this meeting was to discuss the current state of knowledge, the most important biological and clinical questions, and critical next steps for the clinical development of PSMA positron emission tomography (PET) imaging agents and PSMA-targeted radionuclide agents for patients with prostate cancer.

RESULTS

Several major topic areas were discussed including the biology of PSMA, the role of PSMA-targeted PET imaging in prostate cancer, the physics and performance of different PSMA-targeted PET imaging agents, the current state of clinical development of PSMA-targeted radionuclide therapy (RNT) agents, the role of dosimetry in PSMA RNT treatment planning, barriers and challenges in PSMA RNT clinical development, optimization of patient selection for PSMA RNT trials, and promising combination treatment approaches with PSMA RNT.

DISCUSSION

This article summarizes the presentations from the meeting for the purpose of globally disseminating this knowledge to advance the use of PSMA-targeted theranostic agents for imaging and treatment of patients with prostate cancer.

Comparison of different calculation techniques for absorbed dose assessment in patient specific peptide receptor radionuclide therapy

Aim

The present work concerns the comparison of the performances of three systems for dosimetry in RPT that use different techniques for absorbed dose calculation (organ-level dosimetry, voxel-level dose kernel convolution and Monte Carlo simulations). The aim was to assess the importance of the choice of the most adequate calculation modality, providing recommendations about the choice of the computation tool.

Methods

The performances were evaluated both on phantoms and patients in a multi-level approach. Different phantoms filled with a ¹⁷⁷Lu-radioactive solution were used: a homogeneous cylindrical phantom, a phantom with organ-shaped inserts and two cylindrical phantoms with inserts different for shape and volume. A total of 70 patients with NETs treated by PRRT with ¹⁷⁷Lu-DOTATOC were retrospectively analysed.

Results

The comparisons were performed mainly between the mean values of the absorbed dose in the regions of interest. A general better agreement was obtained between Dose kernel convolution and Monte Carlo simulations results rather than between either of these two and organ-level dosimetry, both for phantoms and patients. Phantoms measurements also showed the discrepancies mainly depend on the geometry of the inserts (e.g. shape and volume). For patients, differences were more pronounced than phantoms and higher inter/intra patient variability was observed.

Conclusion

This study suggests that voxel-level techniques for dosimetry calculation are potentially more accurate and personalized than organ-level methods. In particular, a voxel-convolution method provides good results in a short time of calculation, while Monte Carlo based computation should be conducted with very fast calculation systems for a possible use in clinics, despite its intrinsic higher accuracy. Attention to the calculation modality is recommended in case of clinical regions of interest with irregular shape and far from spherical geometry, in which Monte Carlo seems to be more accurate than voxel-convolution methods.

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.

Peptide Receptor Radionuclide Therapy with 67Cu-CuSarTATE Is Highly Efficacious Against a Somatostatin-Positive Neuroendocrine Tumor Model

Peptide receptor radionuclide therapy (PRRT) using radiolabeled octreotate is an effective treatment for somatostatin receptor 2-expressing neuroendocrine tumors. The diagnostic and therapeutic potential of ⁶⁴Cu and ⁶⁷Cu, respectively, offers the possibility of using a single somatostatin receptor-targeted peptide conjugate as a theranostic agent. A sarcophagine cage amine ligand, MeCOSar (5-(8-methyl-3,6,10,13,16,19-hexaaza-bicyclo[6.6.6]icosan-1-ylamino)-5-oxopentanoic acid), conjugated to (Tyr³)-octreotate, called ⁶⁴Cu-CuSarTATE, was demonstrated to be an imaging agent and potential prospective dosimetry tool in 10 patients with neuroendocrine tumors. This study aimed to explore the antitumor efficacy of ⁶⁷Cu-CuSarTATE in a preclinical model of neuroendocrine tumors and compare it with the standard PRRT agent, ¹⁷⁷Lu-LuDOTA-Tyr³-octreotate (¹⁷⁷Lu-LuTATE). Methods: The antitumor efficacy of various doses of ⁶⁷Cu-CuSarTATE in AR42J (rat pancreatic exocrine) tumor-bearing mice was compared with ¹⁷⁷Lu-LuTATE. Results: Seven days after a single administration of ⁶⁷Cu-CuSarTATE (5 MBq), tumor growth was inhibited by 75% compared with vehicle control. Administration of ¹⁷⁷Lu-LuTATE (5 MBq) inhibited tumor growth by 89%. Survival was extended from 12 d in the control group to 21 d after treatment with both ⁶⁷Cu-CuSarTATE and ¹⁷⁷Lu-LuTATE. In a second study, the efficacy of fractionated delivery of PRRT was assessed, comparing the efficacy of 30 MBq of ⁶⁷Cu-CuSarTATE or ¹⁷⁷Lu-LuTATE, either as a single intravenous injection or as two 15-MBq fractions 2 wk apart. Treatment of tumors with 2 fractions significantly improved survival over delivery as a single fraction (⁶⁷Cu-CuSarTATE: 47 vs. 36 d [P = 0.036]; ¹⁷⁷Lu-LuTATE: 46 vs. 29 d [P = 0.040]). Conclusion: This study demonstrates that ⁶⁷Cu-CuSarTATE is well tolerated in BALB/c nude mice and highly efficacious against AR42J tumors in vivo. Administration of ⁶⁷Cu-CuSarTATE and ¹⁷⁷Lu-LuTATE divided into 2 fractions over 2 wk was more efficacious than administration of a single fraction. The antitumor activity of ⁶⁷Cu-CuSarTATE in the AR42J tumor model demonstrated the suitability of this novel agent for clinical assessment in the treatment of somatostatin receptor 2-expressing neuroendocrine tumors.

SCreg: a registration-based platform to compare unicondylar knee arthroplasty SPECT/CT scans

Background

A combination of conventional computed tomography and single photon emitted computed tomography (SPECT/CT) provides simultaneous data on the intensity and location of osteoblastic activity. Currently, since SPECT/CT scans are not spatially aligned, scans following knee arthroplasty are compared by extracting average and maximal values of osteoblastic activity intensity from large subregions of the structure of interest, which leads to a loss of resolution, and hence, information. Therefore, this paper describes the SPECT/CT registration platform (SCreg) based on the principle of image registration to spatially align SPECT/CT scans following unicondylar knee arthroplasty (UKA) and allow full resolution intra-subject and inter-subject comparisons.

Methods

SPECT-CT scans of 20 patients were acquired before and 1 year after UKA. Firstly, scans were pre-processed to account for differences in voxel sizes and divided in volumes of interest. This was followed by optimization of registration parameters according to their volumetric agreement, and alignment using a combination of rigid, affine and non-rigid registration. Finally, radiotracer uptakes were normalized, and differences between pre-operative and post-operative activity were computed for each voxel. Wilcoxon signed rank sum test was performed to compare Dice similarity coefficients pre- and post-registration.

Results

Qualitative and quantitative validation of the platform assessing the correct alignment of SPECT/CT scans resulted in Dice similarity coefficient values over 80% and distances between predefined anatomical landmarks below the fixed threshold of (2;2;0) voxels. Locations of increased and decreased osteoblastic activity obtained during comparisons of osteoblastic activity before and after UKA were mainly consistent with literature.

Conclusions

Thus, a full resolution comparison performed on the platform could assist surgeons and engineers in optimizing surgical parameters in view of bone remodeling, thereby improving UKA survivorship.

BIGDOSE: software for 3D personalized targeted radionuclide therapy dosimetry

Background

Advance 3D quantitative radionuclide imaging techniques boost the accuracy of targeted radionuclide therapy (TRT) dosimetry to voxel level. The goal of this work is to develop a comprehensive 3D dosimetric software, BIGDOSE, with new features of image registration and virtual CT for patient-specific dosimetry.

Methods

BIGDOSE includes a portable graphical user interface written in Python, integrating (I) input of sequential ECT/CT images; (II) segmentation; (III) non-rigid image registration; (IV) curve fitting and voxel-based integration; (V) dose conversion and (VI) 3D dose analysis. The accuracy of the software was evaluated using a simulation study with 9 XCAT phantoms. We simulated SPECT/CT acquisitions at 1, 12, 24, 72 and 144-hrs post In-111 Zevalin injection with inter-scans misalignments using an analytical projector for medium energy general purpose (MEGP) collimator, modeling attenuation, scatter and collimator-detector response. The SPECT data were reconstructed using quantitative OS-EM method. A CT organ-based registration was performed before the dose calculation. Organ absorbed doses for the corresponding Y-90 therapeutic agent were calculated on target organs and compared with those obtained from OLINDA/EXM, using dose measured from GATE as the gold standard. One patient with In-111 DTPAOC injection as well as two patients with Y-90 microsphere embolization were used to demonstrate the clinical effectiveness of our software.

Results

In the simulation, the organ dose errors of BIGDOSE were -9.59%±9.06%, -8.36±5.82%, -23.41%±6.67%, -6.05%±2.06% for liver, spleen, kidneys and lungs, while they were -25.72%±12.52%, -14.93%±10.91%, -28.63%±12.97% and -45.30%±5.84% for OLINDA/EXM. Cumulative dose volume histograms, dose maps and iso-dose contours provided 3D dose distribution information on the simulated and patient data.

Conclusions

BIGDOSE provides a one-stop platform for voxel-based dose estimation with enhanced functions. It is a promising tool to streamline the current clinical TRT dosimetric practice with high accuracy, incorporating 3D personalized imaging information for improved treatment outcome.

Particle filter de-noising of voxel-specific time-activity-curves in personalized 177Lu therapy

BACKGROUND

Currently, there is a high interest in ¹⁷⁷Lu targeted radionuclide therapies, which could be attributed to favorable results obtained from ¹⁷⁷Lu compounds targeting neuro-endocrine and prostate tumors. SPECT based dosimetry could be used for deriving dose values for individual voxels, as is the standard in external-beam radiation-therapy (EBRT). For this a time-activity-curve (TAC) at voxel resolution and also a voxel-wise modeling of radiation energy deposition are necessary. But a voxel-wise determination of TACs is problematic, since several confounding factors exist, such as e.g. poor count-statistics or registration inaccuracies, which add noise to the observed activity states. A particle filter (PF) is a class of methods which applies regularization based on a model of the temporal evolution of activity states. The aim of this study is to introduce the application of PFs for de-noising of per-voxel time-activity curves.

METHODS

We applied a PF for de-noising the TACs of 26 patients, who underwent ¹⁷⁷Lu-DOTATOC or -PSMA therapy. The TACs were obtained from fully-quantitative, serial SPECT(/CT) data, acquired at 4h, 24h, 48h, 72h p.i. The model used in the PF was a mono-exponential decay and its free parameters were determined based on objective criteria. The time-integrated activities (TIA) resulting from the PF (PFF) were compared to the results of a mono-exponential fit (SF) of individual voxels in several volumes of interest (kidneys, spleen, tumors). Additionally, an organ-averaged TIA was derived from whole-organ VOIs and subsequent curve-fitting. This whole-organ TIA was also compared to the whole-organ TIAs obtained from summation of the voxel-wise TIAs from PFF and SF.

RESULTS

The number of particles was set to 1000. Optimal values for noise of observations and noise of the model were 0.25 and 0.5, respectively. The deviation of whole-organ TIAs from conventional organ-based dosimetry and the summation of the voxel-wise TIAs was substantial for SF (kidneys -22.3%, spleen -49.6%, tumor -60.0%), as well as for PFF (kidneys -37.1%, spleen -57.9%, tumor -70.9%). The distribution of voxel-wise half-lives resulting from the PFF method was considerably closer to the organ-averaged value, and the number of implausibly long half-lives (>physical HL) was reduced.

CONCLUSION

The PFF leads to voxel-wise half-lives, which are more plausible than those resulting from SF. However, one has to admit that voxel-wise fitting generally leads to considerable deviations from the organ-averaged TIA as obtained by conventional whole-organ evaluation. Unfortunately, we did not have ground-truth TIA of our patient data and proper ground-truth could even be impossible to obtain. Nevertheless, there are strong indicators that particle filtering can be used for reducing voxel-wise TAC noise.

What You See Is Not What You Get: On the Accuracy of Voxel-Based Dosimetry in Molecular Radiotherapy

Improvements in quantitative SPECT/CT have aroused growing interest in voxel-based dosimetry for radionuclide therapies, because it promises visualization of absorbed doses at a voxel level. In this work, SPECT/CT-based voxel-level dosimetry of a 3-dimensional (3D) printed 2-compartment kidney phantom was performed, and the resulting absorbed dose distributions were examined. Additionally, the potential of the PETPVC partial-volume correction tool was investigated. Methods: Both kidney compartments (70% cortex, 30% medulla) were filled with different activity concentrations, and SPECT/CT imaging was performed. The images were reconstructed using varying settings (iterations, subsets, and postfiltering). On the basis of these activity concentration maps, absorbed dose distributions were calculated with precalculated ¹⁷⁷Lu voxel S values and an empiric kidney half-life. An additional set of absorbed doses was calculated after applying PETPVC for partial-volume correction of the SPECT reconstructions. Results: SPECT/CT imaging blurs the 2 discrete suborgan absorbed dose values into a continuous distribution. Although this effect is slightly improved by applying more iterations, it is enhanced by additional postfiltering. By applying PETPVC, the absorbed dose values are separated into 2 peaks. Although this leads to a better agreement between SPECT/CT-based and nominal values, considerable discrepancies remain. In contrast to the calculated nominal absorbed doses of 7.8 and 1.6 Gy (in the cortex and medulla, respectively), SPECT/CT-based voxel-level dosimetry resulted in mean absorbed doses of 3.0-6.6 Gy (cortex) and 2.7-5.1 Gy (medulla). PETPVC led to improved ranges of 6.1-8.9 Gy (cortex) and 2.1-5.4 Gy (medulla). Conclusion: Our study showed that ¹⁷⁷Lu quantitative SPECT/CT imaging leads to voxel-based dose distributions largely differing from the real organ distribution. SPECT/CT imaging and reconstruction deficiencies might directly translate into unrealistic absorbed dose distributions, thus questioning the reliability of SPECT-based voxel-level dosimetry. Therefore, SPECT/CT reconstructions should be adapted to ensure an accurate quantification of the underlying activity and, therefore, absorbed dose in a volume of interest of the expected object size (e.g., organs, organ substructures, lesions, or voxels). As an example, PETPVC largely improves the match between SPECT/CT-based and nominal dose distributions. In conclusion, the concept of voxel-based dosimetry should be treated with caution. Specifically, one should remember that the absorbed dose distribution is mainly a convolved version of the underlying SPECT reconstruction.

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

Calculation of radiation dosimetry in targeted nuclear medicine therapies is traditionally resource-intensive, requiring multiple posttherapy SPECT acquisitions. An alternative approach is to take advantage of existing pharmacokinetic data from these smaller cohorts to enable dose computation from a single posttreatment scan in a manner that may be applied to a much broader patient population. Methods: In this work, a technical description of simplified dose estimation is presented and applied to the assessment of ¹⁷⁷Lu-prostate-specific membrane antigen (PSMA)-617 therapy for metastatic prostate cancer. By normalizing existing time-activity curves to a single measurement time, it is possible to calculate a mean and range of time-integrated activity values that relate to absorbed radiation dose. To assist with accurate pharmacokinetic modeling of the training cohort, a method for contour-guided image registration was developed. Results: Tissue-specific dose conversion factors for common posttreatment imaging times are reported along with a characterization of added uncertainty in comparison to a traditional serial imaging protocol. Single-time-point dose factors for tumor were determined to be 11.0, 12.1, 13.6, and 15.2 Gy per MBq/mL at image times of 24, 48, 72, and 96 h, respectively. For normal tissues, parotid gland factors were 6.7, 9.4, 13.3, and 19.3 Gy per MBq/mL at those times, and kidneys were 7.1, 10.3, 15.0, and 22.0 Gy per MBq/mL. Tumor dose estimates were most accurate using delayed scanning at times beyond 72 h. Dose to healthy tissues is best characterized by scanning patients in the first 2 d of treatment because of the larger degree of tracer clearance in this early phase. Conclusion: This work demonstrates a means for efficient dose estimation in ¹⁷⁷Lu-PSMA-617 therapy. By providing methods to simplify and potentially automate radiation dosimetry, we hope to accelerate the understanding of radiobiology and development of dose-response models in this unique therapeutic context.

Bone Marrow Absorbed Doses and Correlations with Hematologic Response During 177Lu-DOTATATE Treatments Are Influenced by Image-Based Dosimetry Method and Presence of Skeletal Metastases

This study aimed to compare different image-based methods for bone marrow dosimetry and study the dose–response relationship during treatment with ¹⁷⁷Lu-DOTATATE in patients with and without skeletal metastases. Methods: This study included 46 patients with advanced neuroendocrine tumors treated with at least 2 fractions of ¹⁷⁷Lu-DOTATATE at Sahlgrenska University Hospital. High- and low-uptake compartments were automatically outlined in planar images collected at 2, 24, 48, and 168 h after injection. The bone marrow absorbed doses were calculated from the cross doses of the high- and low-uptake compartments and the self-dose, using the time–activity concentration curve for the low-uptake compartment. This time–activity concentration curve was adjusted using a fixed constant of 1.8 for the planar dosimetry method and using the activity concentrations in vertebral bodies in SPECT images at 24 h after injection of ¹⁷⁷Lu-DOTATATE in 4 hybrid methods: L4-SPECT used the activity concentration in the L4 vertebra, whereas V-SPECT, L-SPECT, and T-SPECT used the median activity concentration in all visible vertebrae, lumbar vertebrae, and thoracic vertebrae, respectively. Results: Using the planar method, L4-SPECT, V-SPECT, L-SPECT, and T-SPECT, the estimated median bone marrow absorbed doses were 0.19, 0.36, 0.40, 0.39, and 0.46 Gy/7.4 GBq, respectively, with respective ranges of 0.12–0.33, 0.15–1.44, 0.19–1.71, 0.21–1.60, and 0.18–2.12 Gy/7.4 GBq. For all methods, the bone marrow absorbed dose significantly correlated with decreased platelet counts. This correlation increased after treatment fraction 2: the Spearman correlation (r_s) were −0.49 for the planar method, −0.61 for L4-SPECT, −0.63 for V-SPECT, −0.63 for L-SPECT, and −0.57 for T-SPECT. A separate analysis revealed an increased correlation for patients without skeletal metastases using the planar method (r_s = −0.67). In contrast, hybrid methods had poor correlations for patients without metastases and stronger correlations for patients with skeletal metastases (r_s = −0.61 to −0.74). The mean bone marrow absorbed doses were 3%–69% higher for patients with skeletal metastases than for patients without. Conclusion: The estimated bone marrow absorbed doses by image-based techniques and the correlation with platelets are influenced by the choice of measured vertebrae and the presence of skeletal metastases.

A dosimetry procedure for organs-at-risk in 177Lu peptide receptor radionuclide therapy of patients with neuroendocrine tumours

PURPOSE

Peptide receptor radionuclide therapy with ¹⁷⁷Lu-DOTATATE has become a standard treatment modality in neuroendocrine tumours (NETs). No consensus has yet been reached however regarding the absorbed dose threshold for lesion response, the absorbed dose limit to organs-at-risk, and the optimal fractionation and activity to be administered. This is partly due to a lack of uniform and comparable dosimetry protocols. The present article details the development of an organ-at-risk dosimetry procedure, which could be implemented and used routinely in a clinical context.

METHODS

Forty-seven patients with NETs underwent ¹⁷⁷Lu-DOTATATE therapy. Three SPECT/CT images were acquired at 4, 24 and 144-192 h post-injection. Three blood samples were obtained together with the SPECT/CT acquisitions and 2 additional samples were obtained around 30 min and 1 h post-injection. A bi-exponential fit was used to compute the source organ time-integrated activity coefficients. Coefficients were introduced into OLINDA/EXM software to compute organ-at-risk absorbed doses. Median values for all patients were computed for absorbed dose coefficient D/A₀ and for late effective half-life T_1/2eff for kidneys, spleen and red marrow.

RESULTS

Dosimetry resulted in a median[interquartile range] of 0.78[0.35], 1.07[0.58] and 0.028[0.010] Gy/GBq for D/A₀ and of 55[9], 71[9] and 52[18] h for T_1/2eff for kidneys, spleen and red marrow respectively.

CONCLUSIONS

A dosimetry procedure for organs-at-risk in ¹⁷⁷Lu-DOTATATE therapy based on serial SPECT/CT images and blood samples can be implemented routinely in a clinical context with limited patient burden. The results obtained were in accordance with those of other centres.

64Cu-SARTATE PET Imaging of Patients with Neuroendocrine Tumors Demonstrates High Tumor Uptake and Retention, Potentially Allowing Prospective Dosimetry for Peptide Receptor Radionuclide Therapy

Imaging of somatostatin receptor expression is an established technique for staging of neuroendocrine neoplasia and determining the suitability of patients for peptide receptor radionuclide therapy. PET/CT using ⁶⁸Ga-labeled somatostatin analogs is superior to earlier agents, but the rapid physical decay of the radionuclide poses logistic and regulatory challenges. ⁶⁴Cu has attractive physical characteristics for imaging and provides a diagnostic partner for the therapeutic radionuclide ⁶⁷Cu. Based on promising preclinical studies, we have performed a first-time-in-humans trial of ⁶⁴Cu-MeCOSar-Tyr³-octreotate (⁶⁴Cu-SARTATE) to assess its safety and ability to localize disease at early and late imaging time-points. Methods: In a prospective trial, 10 patients with known neuroendocrine neoplasia and positive for uptake on ⁶⁸Ga-DOTA-octreotate (⁶⁸Ga-DOTATATE) PET/CT underwent serial PET/CT imaging at 30 min, 1 h, 4 h, and 24 h after injection of ⁶⁴Cu-SARTATE. Adverse reactions were recorded, and laboratory testing was performed during infusion and at 1 and 7 d after imaging. Images were analyzed for lesion and normal-organ uptake and clearance to assess lesion contrast and perform dosimetry estimates. Results: ⁶⁴Cu-SARTATE was well tolerated during infusion and throughout the study, with 3 patients experiencing mild infusion-related events. High lesion uptake and retention were observed at all imaging time-points. There was progressive hepatic clearance over time, providing the highest lesion-to-liver contrast at 24 h. Image quality remained high at this time. Comparison of ⁶⁴Cu-SARTATE PET/CT obtained at 4 h to ⁶⁸Ga-DOTATATE PET/CT obtained at 1 h indicated comparable or superior lesion detection in all patients, especially in the liver. As expected, the highest early physiologic organ uptake was in the kidneys, liver, and spleen. Conclusion: ⁶⁴Cu-SARTATE is safe and has excellent imaging characteristics. High late-retention in tumor and clearance from the liver suggest suitability for diagnostic studies and for prospective dosimetry for ⁶⁷Cu-SARTATE peptide receptor radionuclide therapy, and the half-life of ⁶⁴Cu would also facilitate good-manufacturing-practice production and distribution to sites without access to ⁶⁸Ga.

Dosimetry of 177Lu-PSMA-617 in Metastatic Castration-Resistant Prostate Cancer: Correlations Between Pretherapeutic Imaging and Whole-Body Tumor Dosimetry with Treatment Outcomes

¹⁷⁷Lu-prostate-specific membrane antigen (PSMA)-617 enables targeted delivery of β-particle radiation to prostate cancer. We determined its radiation dosimetry and relationships to pretherapeutic imaging and outcomes. Methods: Thirty patients with prostate cancer receiving ¹⁷⁷Lu-PSMA-617 within a prospective clinical trial (ACTRN12615000912583) were studied. Screening ⁶⁸Ga-PSMA-11 PET/CT demonstrated high PSMA expression in all patients. After therapy, patients underwent quantitative SPECT/CT at 4, 24, and 96 h. Pharmacokinetic uptake and clearance at a voxel level were calculated and translated into absorbed dose using voxel S values. Volumes of interest were drawn on normal tissues and tumor to assess radiation dose, and a whole-body tumor dose was defined. Correlations between PSMA PET/CT parameters, dosimetry, and biochemical and therapeutic response were analyzed to identify relationships between absorbed dose, tumor burden, and patient physiology. Results: Mean absorbed dose to kidneys, submandibular and parotid glands, liver, spleen, and bone marrow was 0.39, 0.44, 0.58, 0.1, 0.06, and 0.11 Gy/MBq, respectively. Median whole-body tumor-absorbed dose was 11.55 Gy and correlated with prostate-specific antigen (PSA) response at 12 wk. A median dose of 14.1 Gy was observed in patients achieving a PSA decline of at least 50%, versus 9.6 Gy for those achieving a PSA decline of less than 50% (P < 0.01). Of 11 patients receiving a tumor dose of less than 10 Gy, only one achieved a PSA response of at least 50%. On screening PSMA PET, whole-body tumor SUV_mean correlated with mean absorbed dose (r = 0.62), and SUV_max of the parotids correlated with absorbed dose (r = 0.67). There was an inverse correlation between tumor volume and mean dose to the parotids (r = -0.41) and kidneys (r = -0.43). The mean parotid dose was also reduced with increasing body mass (r = -0.41) and body surface area (r = -0.37). Conclusion: ¹⁷⁷Lu-PSMA-617 delivers high absorbed doses to tumor, with a significant correlation between whole-body tumor dose and PSA response. Patients receiving less than 10 Gy were unlikely to achieve a fall in PSA of at least 50%. Significant correlations between aspects of screening ⁶⁸Ga-PET/CT and tumor and normal tissue dose were observed, providing a rationale for patient-specific dosing. Reduced salivary and kidney doses were observed in patients with a higher tumor burden. The parotid dose also reduced with increasing body mass and body surface area.

A method for tumor dosimetry based on hybrid planar-SPECT/CT images and semiautomatic segmentation

PURPOSE

A hybrid planar-SPECT/CT method for tumor dosimetry in 177 Lu-DOTATATE therapy, applicable to datasets consisting of multiple conjugate-view images and one SPECT/CT, is developed and evaluated.

METHODS

The imaging protocol includes conjugate-view imaging at 1, 24, 96, and 168 h post infusion (p.i.) and a SPECT/CT acquisition 24 h p.i. The dosimetry method uses the planar images to estimate the shape of the time-activity concentration curve, which is then rescaled to absolute units using the SPECT-derived activity concentration. The resulting time-integrated activity concentration coefficient (TIACC) is used to calculate the tumor-absorbed dose. Semiautomatic segmentation techniques are applied for tumor delineation in both planar and SPECT images, where the planar image segmentation is accomplished using an active-rays-based technique. The selection of tumors is done by visual inspection of planar and SPECT images and applying a set of criteria concerning the tumor visibility and possible interference from superimposed activity uptakes in the planar images. Five different strategies for determining values from planar regions of interest (ROIs), based on entire or partial ROIs, and with and without background correction, are evaluated. Evaluation is performed against a SPECT/CT-based method on data from six patients where sequential conjugate-view and SPECT/CT imaging have been performed in parallel and against ground truths in Monte Carlo simulated images. The patient data are also used to evaluate the interoperator variability and to assess the validity of the developed criteria for tumor selection.

RESULTS

For patient images, the hybrid method produces TIACCs that are on average 6% below those of the SPECT/CT only method, with standard deviations for the relative TIACC differences of 8%-11%. Simulations show that the hybrid and SPECT-based methods estimate the TIACCs to within approximately 10% for tumors larger than around 10 ml, while for smaller tumors, all methods underestimate the TIACCs due to underestimations of the activity concentrations in the SPECT images. The planar image segmentation has a low operator dependence, with a median Dice similarity coefficient of 0.97 between operators. The adopted criteria for tumor selection manage to discriminate the tumors for which the absorbed-dose deviations between the hybrid and SPECT methods are the highest.

CONCLUSIONS

The hybrid method is found suitable for studies of tumor-absorbed doses in radionuclide therapy, provided that selection criteria regarding the visibility and overlapping activities in the planar images are applied.

Dosimetry methods and clinical applications in peptide receptor radionuclide therapy for neuroendocrine tumours: a literature review

Background

The main challenge for systemic radiation therapy using radiopharmaceuticals (SRT) is to optimise the dose delivered to the tumour, while minimising normal tissue irradiation. Dosimetry could help to increase therapy response and decrease toxicity after SRT by individual treatment planning. Peptide receptor radionuclide therapy (PRRT) is an accepted SRT treatment option for irresectable and metastatic neuroendocrine tumours (NET). However, dosimetry in PRRT is not routinely performed, mainly due to the lack of evidence in literature and clinical implementation difficulties. The goal of this review is to provide insight in dosimetry methods and requirements and to present an overview of clinical aspects of dosimetry in PRRT for NET.

Methods

A PubMed query including the search criteria dosimetry, radiation dose, peptide receptor radionuclide therapy, and radionuclide therapy was performed. Articles were selected based on title and abstract, and description of dosimetric approach.

Results

A total of 288 original articles were included. The most important dosimetry methods, their main advantages and limitations, and implications in the clinical setting are discussed. An overview of dosimetry in clinical studies regarding PRRT treatment for NET is provided.

Conclusion

Clinical dosimetry in PRRT is feasible and can result in improved treatment outcomes. Current clinical dosimetry studies focus on safety and apply non-voxel-based dosimetry methods. Personalised treatment using sophisticated dosimetry methods to assess tumour and normal tissue uptake in clinical trials is the next step towards routine dosimetry in PRRT for NET.

Electronic supplementary material

The online version of this article (10.1186/s13550-018-0443-z) contains supplementary material, which is available to authorized users.

A Comparison of 2D and 3D Kidney Absorbed Dose Measures in Patients Receiving 177Lu-DOTATATE

Objective(s):

To investigate and compare quantitative accuracy of kidney absorbed dose measures made from both 2D and 3D imaging in patients receiving ¹⁷⁷Lu-DOTATATE (Lutate) for treatment of neuroendocrine tumours (NETs).

Methods:

Patients receiving Lutate therapy underwent both whole body planar imaging and SPECT/CT imaging over the kidneys at time points 0.5, 4, 24, and 96-120 hours after injection. Planar data were corrected for attenuation using transmission data, and were converted to units of absolute activity via two methods, using either a calibration standard in the field of view or relative to pre-voiding image total counts. Hand drawn regions of interest were used to generate time activity curves and kidney absorbed dose estimates in OLINDA-EXM. Fully quantitative SPECT data were generated using CT-derived corrections for both scatter and attenuation, before correction for dead time and application of a camera specific sensitivity factor to convert data to units of absolute activity. Volumes of interest were defined for kidney using the co-registered x-ray CT, before time activity curves and absorbed dose measures were generated in OLINDA-EXM, both with and without corrections made to the model for patient specific kidney volumes. Quantitative SPECT data were also used to derive dose maps through dose kernel convolution (DKC), which was treated as the gold standard.

Results:

A total of 50 studies were analysed, corresponding to various cycles of treatment from 21 patients. Planar absorbed dose estimates were consistently higher than SPECT derived estimates by, on average, a factor of 3.

Conclusion:

Quantitative SPECT is considered the gold standard approach for organ specific dosimetry however often relies on in house software. As such planar methods for estimating absorbed dose are much more widely available, and in particular, are often the only source of reference in previously published data. For the case of Lutate dosimetry, planar measures may lead to a three-fold increase in measures of kidney absorbed dose.

Deep Learning Renal Segmentation for Fully Automated Radiation Dose Estimation in Unsealed Source Therapy

BACKGROUND

Convolutional neural networks (CNNs) have been shown to be powerful tools to assist with object detection and-like a human observer-may be trained based on a relatively small cohort of reference subjects. Rapid, accurate organ recognition in medical imaging permits a variety of new quantitative diagnostic techniques. In the case of therapy with targeted radionuclides, it may permit comprehensive radiation dose analysis in a manner that would often be prohibitively time-consuming using conventional methods.

METHODS

An automated image segmentation tool was developed based on three-dimensional CNNs to detect right and left kidney contours on non-contrast CT images. Model was trained based on 89 manually contoured cases and tested on a cohort of patients receiving therapy with 177Lu-prostate-specific membrane antigen-617 for metastatic prostate cancer. Automatically generated contours were compared with those drawn by an expert and assessed for similarity based on dice score, mean distance-to-agreement, and total segmented volume. Further, the contours were applied to voxel dose maps computed from post-treatment quantitative SPECT imaging to estimate renal radiation dose from therapy.

RESULTS

Neural network segmentation was able to identify right and left kidneys in all patients with a high degree of accuracy. The system was integrated into the hospital image database, returning contours for a selected study in approximately 90 s. Mean dice score was 0.91 and 0.86 for right and left kidneys, respectively. Poor performance was observed in three patients with cystic kidneys of which only few were included in the training data. No significant difference in mean radiation absorbed dose was observed between the manual and automated algorithms.

CONCLUSION

Automated contouring using CNNs shows promise in providing quantitative assessment of functional SPECT and possibly PET images; in this case demonstrating comparable accuracy for radiation dose interpretation in unsealed source therapy relative to a human observer.

Development and Validation of RAPID: A Patient-Specific Monte Carlo Three-Dimensional Internal Dosimetry Platform

This work describes the development and validation of a patient-specific Monte Carlo internal dosimetry platform called RAPID (Radiopharmaceutical Assessment Platform for Internal Dosimetry). RAPID utilizes serial PET/CT or SPECT/CT images to calculate voxelized three-dimensional (3D) internal dose distributions with the Monte Carlo code Geant4. RAPID's dosimetry calculations were benchmarked against previously published S-values and specific absorbed fractions (SAFs) calculated for monoenergetic photon and electron sources within the Zubal phantom and for S-values calculated for a variety of radionuclides within spherical tumor phantoms with sizes ranging from 1 to 1000 g. The majority of the S-values and SAFs calculated in the Zubal Phantom were within 5% of the previously published values with the exception of a few 10 keV photon SAFs that agreed within 10%, and one value within 16%. The S-values calculated in the spherical tumor phantoms agreed within 2% for ¹⁷⁷Lu, ¹³¹I, ¹²⁵I, ¹⁸F, and ⁶⁴Cu, within 3.5% for ²¹¹At and ²¹³Bi, within 6.5% for ¹⁵³Sm, ¹¹¹In, ⁸⁹Zr, and ²²³Ra, and within 9% for ⁹⁰Y, ⁶⁸Ga, and ¹²⁴I. In conclusion, RAPID is capable of calculating accurate internal dosimetry at the voxel-level for a wide variety of radionuclides and could be a useful tool for calculating patient-specific 3D dose distributions.

SPECT/CT: an update on technological developments and clinical applications

Functional nuclear medicine imaging with single-photon emission CT (SPECT) in combination with anatomical CT has been commercially available since the beginning of this century. The combination of the two modalities has improved both the sensitivity and specificity of many clinical applications and CT in conjunction with SPECT that allows for spatial overlay of the SPECT data on good anatomy images. Introduction of diagnostic CT units as part of the SPECT/CT system has also potentially allowed for a more cost-efficient use of the equipment. Most of the SPECT systems available are based on the well-known Anger camera principle with NaI(Tl) as a scintillation material, parallel-hole collimators and multiple photomultiplier tubes, which, from the centroid of the scintillation light, determine the position of an event. Recently, solid-state detectors using cadmium-zinc-telluride became available and clinical SPECT cameras employing multiple pinhole collimators have been developed and introduced in the market. However, even if new systems become available with better hardware, the SPECT reconstruction will still be affected by photon attenuation and scatter and collimator response. Compensation for these effects is needed even for qualitative studies to avoid artefacts leading to false positives. This review highlights the recent progress for both new SPECT cameras systems as well as for various data-processing and compensation methods.

Evaluation of sequential SPECT and CT for targeted radionuclide therapy dosimetry

PURPOSE

In targeted radionuclide therapy (TRT), a prior knowledge of the absorbed dose biodistribution is essential for pre-therapy treatment planning. Previously, we showed that non-rigid organ-by-organ registration in sequential quantitative SPECT images improved dose estimation. This study aims to investigate if sequential CT can further improve TRT dosimetric accuracy.

METHODS

We simulated SPECT/CT acquisitions at 1, 12, 24, 72 and 144 h In-111 Zevalin post-injection using an analytical MEGP projector, modeling attenuation, scatter and collimator-detector response. We later recruited a patient injected with 222 MBq In-111 DTPAOC imaged at 3 SPECT/CT sessions for clinical evaluations. Four registration schemes were evaluated: whole-body-based registration performed on sequential (1) SPECT (WB-SPECT) or (2) CT (WB-CT) images; organ-based registration applied on organs individually segmented from sequential (3) SPECT (O-SPECT) or (4) CT (O-CT) images. Voxel-by-voxel integration was performed followed by Y-90 voxel-S-kernel convolution. Organ-absorbed doses, iso-dose curves, dose-volume histograms (DVHs) were generated for targeted organs for analysis.

RESULTS

In simulation study, organ-absorbed dose errors were (- 8.66 ± 2.83)%, (- 2.51 ± 3.69)%, (- 9.23 ± 3.28)%, (- 7.17 ± 2.53)% for liver, (- 14.81 ± 4.91)%, (- 3.60 ± 4.37)%, (- 18.13 ± 4.44)%, (- 11.34 ± 4.22)% for spleen, for O-SPECT, O-CT, WB-SPECT and WB-CT registrations, respectively. For all organs, O-CT showed superior results. Results of iso-dose contour, DVHs were in accordance with the organ-absorbed doses. In clinical studies, the results were also consistent which showed O-CT method deviated the most from the result with no registration.

CONCLUSIONS

We conclude that if both sequential SPECT/CT scans are available, CT organ-based registration method can more effectively improve the 3D dose estimation. Sequential low-dose CT scans might be considered to be included in the standard TRT protocol.