Pre- and post-treatment image-based dosimetry in90Y-microsphere radioembolization using the TOPAS Monte Carlo toolkit

Advances and Emerging Techniques in Y-90 Radioembolization for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the fourth leading cause of cancer deaths worldwide. Despite the high incidence of HCC, mortality remains high, with an estimated 5-year survival rate of less than 20%. Surgical resection represents a potential curative treatment for HCC; however, less than 20% of patients with HCC are candidates for surgical resection. In patients with unresectable HCC, Yttrium-90 (Y90) transarterial radioembolization (TARE) has emerged as an innovative treatment option. This locoregional therapy delivers high doses of radiation directly to liver tumors via intra-arterial injection, allowing for the targeted destruction of malignant cells while sparing surrounding healthy tissue. In this review, we will explore the latest advances in Y90 TARE for the treatment of HCC, focusing on key developments such as the following: (1) improvements in radiation lobectomy and segmentectomy techniques, (2) the introduction of personalized dosimetry, (3) the integration of combination therapies, (4) the use of imageable microspheres, (5) pressure-enabled Y90 delivery systems, and (6) the application of Y90 surrogates.

SPECT/CT Dosimetry of Bronchial Artery 99mTc Macroaggregated Albumin Injection in Pulmonary Malignancies: Feasibility Evaluation of Bronchial Artery 90Y Radioembolization

Background External beam radiation therapy for primary and secondary pulmonary malignancies has limited utility for treating ultracentral tumors (ie, adjacent to the proximal bronchial tree or heart) or multiple metastases due to either radiation to central organs at risk (OARs) or extensive lung tissue exposure. Bronchial artery yttrium 90 (90Y) radioembolization may be a therapeutic option for these patients. Purpose To evaluate the feasibility of bronchial artery 90Y radioembolization using technetium 99m (99mTc) macroaggregated albumin (MAA) injection as a surrogate for 90Y microspheres and to use SPECT/CT dosimetry to assess 99mTc-MAA distribution and calculated anticipated 90Y doses to tumors and OARs. Materials and Methods In this prospective clinical trial, study participants with either primary lung cancer or pulmonary metastases were enrolled between August 2020 and October 2023 from a single academic medical center. All participants underwent bronchial artery embolization for malignancy-induced hemoptysis prophylaxis. 99mTc-MAA was injected via bronchial arteries, followed by bland embolization. SPECT/CT imaging and Monte Carlo simulations were used to evaluate 99mTc-MAA distribution and predict 90Y doses. Predicted 90Y doses to tumors and OARs are reported as means ± SDs. Results Eight participants (mean age, 63.0 years ± 13.58; six [75%] male participants) were included. All had ultracentral tumors, and four had four or more tumors. SPECT/CT revealed a concentrated 99mTc-MAA accumulation in tumors, with a mean tumor-to-normal tissue ratio of 22.71 ± 20.17. Simulations indicated that a 90Y biologically effective dose (α/β ratio of 10 Gy) of 175.7-3173.6 Gy (mean, 778.8 Gy ± 981.9) could be delivered to all tumors while remaining under OAR toxicity thresholds. Conclusion SPECT/CT dosimetric analysis of 99mTc-MAA injected via the bronchial artery for pulmonary malignancy suggests that bronchial artery 90Y radioembolization is feasible and could be an alternative treatment for patients unable to receive external beam radiation therapy. ClinicalTrials.gov Identifier: NCT04105283 © RSNA, 2025 Supplemental material is available for this article.

Monte Carlo methods for medical imaging research

In radiation-based medical imaging research, computational modeling methods are used to design and validate imaging systems and post-processing algorithms. Monte Carlo methods are widely used for the computational modeling as they can model the systems accurately and intuitively by sampling interactions between particles and imaging subject with known probability distributions. This article reviews the physics behind Monte Carlo methods, their applications in medical imaging, and available MC codes for medical imaging research. Additionally, potential research areas related to Monte Carlo for medical imaging are discussed.

Monte Carlo dosimetric analyses on the use of90Y-IsoPet intratumoral therapy in canine subjects

Objective.To investigate different dosimetric aspects of90Y-IsoPet™ intratumoral therapy in canine soft tissue sarcomas, model the spatial spread of the gel post-injection, evaluate absorbed dose to clinical target volumes, and assess dose distributions and treatment efficacy.Approach.Six canine cases treated with90Y-IsoPet™ for soft tissue sarcoma at the Veterinary Health Center, University of Missouri are analyzed in this retrospective study. The dogs received intratumoral IsoPet™ injections, following a grid pattern to achieve a near-uniform dose distribution in the clinical target volume. Two dosimetry methods were performed retrospectively using the Monte Carlo toolkit OpenTOPAS: imaging-based dosimetry obtained from post-injection PET/CT scans, and stylized phantom-based dosimetry modeled from the planned injection points to the gross tumor volume. For the latter, a Gaussian parameter with variable sigma was introduced to reflect the spatial spread of IsoPet™. The two methods were compared using dose-volume histograms (DVHs) and dose homogeneity, allowing an approximation of the closest sigma for the spatial spread of the gel post-injection. In addition, we compared Monte Carlo-based dosimetry with voxel S-value (VSV)-based dosimetry to investigate the dosimetric differences.Main results.Imaging-based dosimetry showed differences between Monte Carlo and VSV calculations in tumor high-density areas with higher self-absorption. Stylized phantom-based dosimetry indicated a more homogeneous target dose with increasing sigma. The sigma approximation of the90Y-IsoPet™ post-injection gel spread resulted in a median sigma of approximately 0.44 mm across all cases to reproduce the dose heterogeneity observed in Monte Carlo calculations.Significance.The results indicate that dose modeling based on planned injection points can serve as a first-order approximation for the delivered dose in90Y-IsoPet™ therapy for canine soft tissue sarcomas. The dosimetry evaluation highlights the non-uniformity of absorbed doses despite the gel spread, emphasizing the importance of considering tumor dose heterogeneity in treatment evaluation. Our findings suggest that using Monte Carlo for dose calculation seems more suitable for this type of tumor where high-density areas might play an important role in dosimetry.

Characterizing the voxel-based approaches in radioembolization dosimetry with reDoseMC

BACKGROUND

Yttrium-90 (90 Y $^{90}{\rm {Y}}$ ) represents the primary radioisotope used in radioembolization procedures, while holmium-166 (166 Ho $^{166}{\rm {Ho}}$ ) is hypothesized to serve as a viable substitute for90 Y $^{90}{\rm {Y}}$ due to its comparable therapeutic potential and improved quantitative imaging. Voxel-based dosimetry for these radioisotopes relies on activity images obtained through PET or SPECT and dosimetry methods, including the voxel S-value (VSV) and the local deposition method (LDM). However, the evaluation of the accuracy of absorbed dose calculations has been limited by the use of non-ideal reference standards and investigations restricted to the liver. The objective of this study was to expand upon these dosimetry characterizations by investigating the impact of image resolutions, voxel sizes, target volumes, and tissue materials on the accuracy of90 Y $^{90}{\rm {Y}}$ and166 Ho $^{166}{\rm {Ho}}$ dosimetry techniques.

METHODS

A specialized radiopharmaceutical dosimetry software called reDoseMC was developed using the Geant4 Monte Carlo toolkit and validated by benchmarking the generated90 Y $^{90}{\rm {Y}}$ kernels with published data. The decay spectra of both90 Y $^{90}{\rm {Y}}$ and166 Ho $^{166}{\rm {Ho}}$ were also compared. Multiple VSV kernels were generated for the liver, lungs, soft tissue, and bone for isotropic voxel sizes of 1 mm, 2 mm, and 4 mm. Three theoretical phantom setups were created with 20 or 40 mm activity and mass density inserts for the same three voxel sizes. To replicate the limited spatial resolutions present in PET and SPECT images, image resolutions were modeled using a 3D Gaussian kernel with a Full Width at Half Maximum (FWHM) ranging from 0 to 16 mm and with no added noise. The VSV and LDM dosimetry methods were evaluated by characterizing their respective kernels and analyzing their absorbed dose estimates calculated on theoretical phantoms. The ground truth for these estimations was calculated using reDoseMC.

RESULTS

The decay spectra obtained through reDoseMC showed less than a 1% difference when compared to previously published experimental data for energies below 1.9 MeV in the case of90 Y $^{90}{\rm {Y}}$ and less than 1% for energies below 1.5 MeV for166 Ho $^{166}{\rm {Ho}}$ . Additionally, the validation kernels for90 Y $^{90}{\rm {Y}}$ VSV exhibited results similar to those found in published Monte Carlo codes, with source dose depositions having less than a 3% error margin. Resolution thresholds (FWHM thresh s ${\rm {FWHM}}_\mathrm{thresh}{\rm {s}}$ ), defined as resolutions that resulted in similar dose estimates between the LDM and VSV methods, were observed for90 Y $^{90}{\rm {Y}}$ . They were 1.5 mm for bone, 2.5 mm for soft tissue and liver, and 8.5 mm for lungs. For166 Ho $^{166}{\rm {Ho}}$ , the accuracy of absorbed dose deposition was found to be dependent on the contributions of absorbed dose from photons. Volume errors due to variations in voxel size impacted the final dose estimates. Larger target volumes yielded more accurate mean doses than smaller volumes. For both radioisotopes, the radial dose profiles for the VSV and LDM approximated but never matched the reference standard.

CONCLUSIONS

reDoseMC was developed and validated for radiopharmaceutical dosimetry. The accuracy of voxel-based dosimetry was found to vary widely with changes in image resolutions, voxel sizes, chosen target volumes, and tissue material; hence, the standardization of dosimetry protocols was found to be of great importance for comparable dosimetry analysis.

MIDOS: a novel stochastic model towards a treatment planning system for microsphere dosimetry in liver tumors

PURPOSE

Transarterial radioembolization (TARE) procedures treat liver tumors by injecting radioactive microspheres into the hepatic artery. Currently, there is a critical need to optimize TARE towards a personalized dosimetry approach. To this aim, we present a novel microsphere dosimetry (MIDOS) stochastic model to estimate the activity delivered to the tumor(s), normal liver, and lung.

METHODS

MIDOS incorporates adult male/female liver computational phantoms with the hepatic arterial, hepatic portal venous, and hepatic venous vascular trees. Tumors can be placed in both models at user discretion. The perfusion of microspheres follows cluster patterns, and a Markov chain approach was applied to microsphere navigation, with the terminal location of microspheres determined to be in either normal hepatic parenchyma, hepatic tumor, or lung. A tumor uptake model was implemented to determine if microspheres get lodged in the tumor, and a probability was included in determining the shunt of microspheres to the lung. A sensitivity analysis of the model parameters was performed, and radiation segmentectomy/lobectomy procedures were simulated over a wide range of activity perfused. Then, the impact of using different microspheres, i.e., SIR-Sphere®, TheraSphere®, and QuiremSphere®, on the tumor-to-normal ratio (TNR), lung shunt fraction (LSF), and mean absorbed dose was analyzed.

RESULTS

Highly vascularized tumors translated into increased TNR. Treatment results (TNR and LSF) were significantly more variable for microspheres with high particle load. In our scenarios with 1.5 GBq perfusion, TNR was maximum for TheraSphere® at calibration time in segmentectomy/lobar technique, for SIR-Sphere® at 1-3 days post-calibration, and regarding QuiremSphere® at 3 days post-calibration.

CONCLUSION

This novel approach is a decisive step towards developing a personalized dosimetry framework for TARE. MIDOS assists in making clinical decisions in TARE treatment planning by assessing various delivery parameters and simulating different tumor uptakes. MIDOS offers evaluation of treatment outcomes, such as TNR and LSF, and quantitative scenario-specific decisions.