Linear Boltzmann equation solver for voxel-level dosimetry in radiopharmaceutical therapy: Comparison with Monte Carlo and kernel convolution.

BACKGROUND: With recent interest in patient-specific dosimetry for radiopharmaceutical therapy (RPT) and selective internal radiation therapy (SIRT), an increasing number of voxel-based algorithms are being evaluated. Monte Carlo (MC) radiation transport, generally considered to be the most accurate among different methods for voxel-level absorbed dose estimation, can be computationally inefficient for routine clinical use. PURPOSE: This work demonstrates a recently implemented grid-based linear Boltzmann transport equation (LBTE) solver for fast and accurate voxel-based dosimetry in RPT and SIRT and benchmarks it against MC. METHODS: A deterministic LBTE solver (Acuros MRT) was implemented within a commercial RPT dosimetry package (Velocity 4.1). The LBTE is directly discretized using an adaptive mesh refined grid and then the coupled photon-electron radiation transport is iteratively solved inside specified volumes to estimate radiation doses from both photons and charged particles in heterogeneous media. To evaluate the performance of the LBTE solver for RPT and SIRT applications, 177 Lu SPECT/CT, 90 Y PET/CT, and 131 I SPECT/CT images of phantoms and patients were used. Multiple lesions (2-1052 mL) and normal organs were delineated for each study. Voxel dosimetry was performed with the LBTE solver, dose voxel kernel (DVK) convolution with density correction, and a validated in-house MC code using the same time-integrated activity and density maps as input to the different dose engines. The resulting dose maps, difference maps, and dose-volume-histogram (DVH) metrics were compared, to assess the voxel-level agreement. Evaluation of mean absorbed dose included comparison with structure-level estimates from OLINDA. RESULTS: In the phantom inserts/compartments, the LBTE solver versus MC and DVK convolution demonstrated good agreement with mean absorbed dose and DVH metrics agreeing to within 5% except for the D90 and D70 metrics of a very low activity concentration insert of 90 Y where the agreement was within 15%. In the patient studies (five patients imaged after 177 Lu DOTATATE RPT, five after 90 Y SIRT, and two after 131 I radioimmunotherapy), in general, there was better agreement between the LBTE solver and MC than between LBTE solver and DVK convolution for mean absorbed dose and voxel-level evaluations. Across all patients for all three radionuclides, for soft tissue structures (kidney, liver, lesions), the mean absorbed dose estimates from the LBTE solver were in good agreement with those from MC (median difference < 1%, maximum 9%) and those from DVK (median difference < 5%, maximum 9%). The LBTE and OLINDA estimates for mean absorbed dose in kidneys and liver agreed to within 10%, but differences for lesions were larger with a maximum 14% for 177 Lu, 23% for 90 Y, and 26% for 131 I. For bone regions, the agreement in mean absorbed doses between LBTE and both MC and DVK were similar (median < 11%, max 11%) while for lung the agreement between LBTE and MC (median < 1%, max 8%) was substantially better than between LBTE and DVK (median < 16%, max 33%). Voxel level estimates for soft tissue structures also showed good agreement between the LBTE solver and both MC and DVK with a median difference < 5% (maximum < 13%) for the DVH metrics with all three radionuclides. The largest difference in DVH metrics was for the D90 and D70 metric in lung and bone where the uptake was low. Here, the difference between LBTE and MC had a median value < 14% (maximum 23%) for bone and < 4% (maximum 37%) for lung, while the corresponding differences between LBTE and DVK were < 23% (maximum 31%) and < 67% (maximum 313%), respectively. For a typical patient with a matrix size of 166 × 166 × 129 (voxel size 3 × 3 × 3 mm3 ), voxel dosimetry using the LBTE solver was as fast as ∼2 min on a desktop computer. CONCLUSION: Having established good agreement between the LBTE solver and MC for RPT and SIRT applications, the LBTE solver is a viable option for voxel dosimetry that can be faster than MC. Further analysis is being performed to encompass the broad range of radionuclides and conditions encountered clinically.

Lesion Dosimetry for [177Lu]Lu-PSMA-617 Radiopharmaceutical Therapy Combined with Stereotactic Body Radiotherapy in Patients with Oligometastatic Castration-Sensitive Prostate Cancer.

A single-institution prospective pilot clinical trial was performed to demonstrate the feasibility of combining [177Lu]Lu-PSMA-617 radiopharmaceutical therapy (RPT) with stereotactic body radiotherapy (SBRT) for the treatment of oligometastatic castration-sensitive prostate cancer. Methods: Six patients with 9 prostate-specific membrane antigen (PSMA)-positive oligometastases received 2 cycles of [177Lu]Lu-PSMA-617 RPT followed by SBRT. After the first intravenous infusion of [177Lu]Lu-PSMA-617 (7.46 ± 0.15 GBq), patients underwent SPECT/CT at 3.2 ± 0.5, 23.9 ± 0.4, and 87.4 ± 12.0 h. Voxel-based dosimetry was performed with calibration factors (11.7 counts per second/MBq) and recovery coefficients derived from in-house phantom experiments. Lesions were segmented on baseline PSMA PET/CT (50% SUVmax). After a second cycle of [177Lu]Lu-PSMA-617 (44 ± 3 d; 7.50 ± 0.10 GBq) and an interim PSMA PET/CT scan, SBRT (27 Gy in 3 fractions) was delivered to all PSMA-avid oligometastatic sites, followed by post-PSMA PET/CT. RPT and SBRT voxelwise dose maps were scaled (α/ß = 3 Gy; repair half-time, 1.5 h) to calculate the biologically effective dose (BED). Results: All patients completed the combination therapy without complications. No grade 3+ toxicities were noted. The median of the lesion SUVmax as measured on PSMA PET was 16.8 (interquartile range [IQR], 11.6) (baseline), 6.2 (IQR, 2.7) (interim), and 2.9 (IQR, 1.4) (post). PET-derived lesion volumes were 0.4-1.7 cm3 The median lesion-absorbed dose (AD) from the first cycle of [177Lu]Lu-PSMA-617 RPT (ADRPT) was 27.7 Gy (range, 8.3-58.2 Gy; corresponding to 3.7 Gy/GBq, range, 1.1-7.7 Gy/GBq), whereas the median lesion AD from SBRT was 28.1 Gy (range, 26.7-28.8 Gy). Spearman rank correlation, ρ, was 0.90 between the baseline lesion PET SUVmax and SPECT SUVmax (P = 0.005), 0.74 (P = 0.046) between the baseline PET SUVmax and the lesion ADRPT, and -0.81 (P = 0.022) between the lesion ADRPT and the percent change in PET SUVmax (baseline to interim). The median for the lesion BED from RPT and SBRT was 159 Gy (range, 124-219 Gy). ρ between the BED from RPT and SBRT and the percent change in PET SUVmax (baseline to post) was -0.88 (P = 0.007). Two cycles of [177Lu]Lu-PSMA-617 RPT contributed approximately 40% to the maximum BED from RPT and SBRT. Conclusion: Lesional dosimetry in patients with oligometastatic castration-sensitive prostate cancer undergoing [177Lu]Lu-PSMA-617 RPT followed by SBRT is feasible. Combined RPT and SBRT may provide an efficient method to maximize the delivery of meaningful doses to oligometastatic disease while addressing potential microscopic disease reservoirs and limiting the dose exposure to normal tissues.

Personalizing Radiotherapy Prescription Dose Using Genomic Markers of Radiosensitivity and Normal Tissue Toxicity in NSCLC.

INTRODUCTION: Cancer sequencing efforts have revealed that cancer is the most complex and heterogeneous disease that affects humans. However, radiation therapy (RT), one of the most common cancer treatments, is prescribed on the basis of an empirical one-size-fits-all approach. We propose that the field of radiation oncology is operating under an outdated null hypothesis: that all patients are biologically similar and should uniformly respond to the same dose of radiation. METHODS: We have previously developed the genomic-adjusted radiation dose, a method that accounts for biological heterogeneity and can be used to predict optimal RT dose for an individual patient. In this article, we use genomic-adjusted radiation dose to characterize the biological imprecision of one-size-fits-all RT dosing schemes that result in both over- and under-dosing for most patients treated with RT. To elucidate this inefficiency, and therefore the opportunity for improvement using a personalized dosing scheme, we develop a patient-specific competing hazards style mathematical model combining the canonical equations for tumor control probability and normal tissue complication probability. This model simultaneously optimizes tumor control and toxicity by personalizing RT dose using patient-specific genomics. RESULTS: Using data from two prospectively collected cohorts of patients with NSCLC, we validate the competing hazards model by revealing that it predicts the results of RTOG 0617. We report how the failure of RTOG 0617 can be explained by the biological imprecision of empirical uniform dose escalation which results in 80% of patients being overexposed to normal tissue toxicity without potential tumor control benefit. CONCLUSIONS: Our data reveal a tapestry of radiosensitivity heterogeneity, provide a biological framework that explains the failure of empirical RT dose escalation, and quantify the opportunity to improve clinical outcomes in lung cancer by incorporating genomics into RT.

SPECT/CT image-based dosimetry for Yttrium-90 radionuclide therapy: Application to treatment response.

This work demonstrates the efficacy of voxel-based 90 Y microsphere dosimetry utilizing post-therapy SPECT/CT imaging and applies it to the prediction of treatment response for the management of patients with hepatocellular carcinoma (HCC). A 90 Y microsphere dosimetry navigator (RapidSphere) within a commercial platform (Velocity, Varian Medical Systems) was demonstrated for three microsphere cases that were imaged using optimized bremsstrahlung SPECT/CT. For each case, the 90 Y SPECT/CT was registered to follow-up diagnostic MR/CT using deformable image registration. The voxel-based dose distribution was computed using the local deposition method with known injected activity. The system allowed the visualization of the isodose distributions on any of the registered image datasets and the calculation of dose-volume histograms (DVHs). The dosimetric analysis illustrated high local doses that are characteristic of blood-flow directed brachytherapy. In the first case, the HCC mass demonstrated a complete response to treatment indicated by a necrotic region in follow-up MR imaging. This result was dosimetrically predicted since the gross tumor volume (GTV) was well covered by the prescription isodose volume (V150 Gy = 85%). The second case illustrated a partial response to treatment which was characterized by incomplete necrosis of an HCC mass and a remaining area of solid enhancement in follow-up MR imaging. This result was predicted by dosimetric analysis because the GTV demonstrated incomplete coverage by the prescription isodose volume (V470 Gy = 18%). The third case demonstrated extrahepatic activity. The dosimetry indicated that the prescription (125 Gy) isodose region extended outside of the liver into the duodenum (178 Gy maximum dose). This was predictive of toxicity as the patient later developed a duodenal ulcer. The ability to predict outcomes and complications using deformable image registration, calculated isodose distributions, and DVHs, points to the clinical utility of patient-specific dose calculations for 90 Y radioembolization treatment planning.