A cell-based dosimetry model for radium-223 dichloride therapy using bone micro-CT images and GATE simulations

A GATE simulation study for dosimetry in cancer cell and micrometastasis from the 225Ac decay chain

BACKGROUND

Radiopharmaceutical therapy (RPT) with alpha-emitting radionuclides has shown great promise in treating metastatic cancers. The successive emission of four alpha particles in the ²²⁵Ac decay chain leads to highly targeted and effective cancer cell death. Quantifying cellular dosimetry for ²²⁵Ac RPT is essential for predicting cell survival and therapeutic success. However, the leading assumption that all ²²⁵Ac progeny remain localized at their target sites likely overestimates the absorbed dose to cancer cells. To address limitations in existing semi-analytic approaches, this work evaluates S-values for ²²⁵Ac's progeny radionuclides with GATE Monte Carlo simulations.

METHODS

The cellular geometries considered were an individual cell (10 µm diameter with a nucleus of 8 µm diameter) and a cluster of cells (micrometastasis) with radionuclides localized in four subcellular regions: cell membrane, cytoplasm, nucleus, or whole cell. The absorbed dose to the cell nucleus was scored, and self- and cross-dose S-values were derived. We also evaluated the total absorbed dose with various degrees of radiopharmaceutical internalization and retention of the progeny radionuclides ²²¹Fr (t_1/2 = 4.80 m) and ²¹³Bi (t_1/2 = 45.6 m).

RESULTS

For the cumulative ²²⁵Ac decay chain, our self- and cross-dose nuclear S-values were both in good agreement with S-values published by MIRDcell, with per cent differences ranging from - 2.7 to - 8.7% for the various radionuclide source locations. Source location had greater effects on self-dose S-values than the intercellular cross-dose S-values. Cumulative ²²⁵Ac decay chain self-dose S-values increased from 0.167 to 0.364 GyBq^-1 s^-1 with radionuclide internalization from the cell surface into the cell. When progeny migration from the target site was modelled, the cumulative self-dose S-values to the cell nucleus decreased by up to 71% and 21% for ²²¹Fr and ²¹³Bi retention, respectively.

CONCLUSIONS

Our GATE Monte Carlo simulations resulted in cellular S-values in agreement with existing MIRD S-values for the alpha-emitting radionuclides in the ²²⁵Ac decay chain. To obtain accurate absorbed dose estimates in ²²⁵Ac studies, accurate understanding of daughter migration is critical for optimized injected activities. Future work will investigate other novel preclinical alpha-emitting radionuclides to evaluate therapeutic potency and explore realistic cellular geometries corresponding to targeted cancer cell lines.

Heterogeneity of dose distribution in normal tissues in case of radiopharmaceutical therapy with alpha-emitting radionuclides

Heterogeneity of dose distribution has been shown at different spatial scales in diagnostic nuclear medicine. In cancer treatment using new radiopharmaceuticals with alpha-particle emitters, it has shown an extensive degree of dose heterogeneity affecting both tumour control and toxicity of organs at risk. This review aims to provide an overview of generalized internal dosimetry in nuclear medicine and highlight the need of consideration of the dose heterogeneity within organs at risk. The current methods used for patient dosimetry in radiopharmaceutical therapy are summarized. Bio-distribution and dose heterogeneities of alpha-particle emitting pharmaceutical 223Ra (Xofigo) within bone tissues are presented as an example. In line with the strategical research agendas of the Multidisciplinary European Low Dose Initiative (MELODI) and the European Radiation Dosimetry Group (EURADOS), future research direction of pharmacokinetic modelling and dosimetry in patient radiopharmaceutical therapy are recommended.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Modeling bystander effects that cause growth delay of breast cancer xenografts in bone marrow of mice treated with radium-223

RATIONALE

The role of radiation-induced bystander effects in cancer therapy with alpha-particle emitting radiopharmaceuticals remains unclear. With renewed interest in using alpha-particle emitters to sterilize disseminated tumor cells, micrometastases, and tumors, a better understanding of the direct effects of alpha particles and the contribution of the bystander responses they induce is needed to refine dosimetric models that help predict clinical benefit. Accordingly, this work models and quantifies the relative importance of direct effects (DE) and bystander effects (BE) in the growth delay of human breast cancer xenografts observed previously in the tibiae of mice treated with 223RaCl2.

METHODS

A computational model of MDA-MB-231 and MCF-7 human breast cancer xenografts in the tibial bone marrow of mice administered 223RaCl2 was created. A Monte Carlo radiation transport simulation was performed to assess individual cell absorbed doses. The responses of the breast cancer cells to direct alpha particle irradiation and gamma irradiation were needed as input data for the model and were determined experimentally using a colony-forming assay and compared to the responses of preosteoblast MC3T3-E1 and osteocyte-like MLO-Y4 bone cells. Using these data, a scheme was devised to simulate the dynamic proliferation of the tumors in vivo, including DE and BE propagated from the irradiated cells. The parameters of the scheme were estimated semi-empirically to fit experimental tumor growth.

RESULTS

A robust BE component, in addition to a much smaller DE component, was required to simulate the in vivo tumor proliferation. We also found that the relative biological effectiveness (RBE) for cell killing by alpha particle radiation was greater for the bone cells than the tumor cells.

CONCLUSION

This modeling study demonstrates that DE of radiation alone cannot explain experimental observations of 223RaCl2-induced growth delay of human breast cancer xenografts. Furthermore, while the mechanisms underlying BE remain unclear, the addition of a BE component to the model is necessary to provide an accurate prediction of the growth delay. More complex models are needed to further comprehend the extent and complexity of 223RaCl2-induced BE.

Effect of microdistribution of alpha and beta-emitters in targeted radionuclide therapies on delivered absorbed dose in a GATE model of bone marrow

Acute hematologic toxicity is a frequent adverse effect of beta-emitter targeted radionuclide therapies (TRTs). Alpha emitters have the potential of delivering high linear energy transfer (LET) radiation to the tumor attributed to its shorter range. Antibody-based TRTs have increased blood-pool half-lives, and therefore increased marrow toxicity, which is a particular concern with alpha emitters. Accurate 3D absorbed dose calculations focusing on the interface region of blood vessels and bone can elucidate energy deposition patterns. Firstly, a cylindrical geometry model with a central blood vessel embedded in the trabecular tissue was modeled. Monte Carlo simulations in GATE were performed considering beta (177Lu, 90Y) and alpha emitters (211At, 225Ac) as sources restricted to the blood pool. Subsequently, the radioactive sources were added in the trabecular bone compartment in order to model bone marrow metastases infiltration (BMMI). Radial profiles, dose-volume histograms and voxel relative differences were used to evaluate the absorbed dose results. We demonstrated that alpha emitters have a higher localized energy deposition compared to beta emitters. In the cylindrical geometry model, when the sources are confined to the blood pool, the dose to the trabecular bone is greater for beta emitting radionuclides, as alpha emitters deposit the majority of their energy within 70 μm of the vessel wall. In the BMMI model, alpha emitters have a lower dose to untargeted trabecular bone. Our results suggest that when alpha emitters are restricted to the blood pool, as when labeled to antibodies, hematologic toxicities may be lower than expected due to differences in the microdistribution of delivered absorbed dose.