Monte Carlo simulations of energy deposition and DNA damage using TOPAS-nBio

Estimation of the relative biological effectiveness for double strand break induction of clinical kilovoltage beams using Monte Carlo simulations

BACKGROUND

The Relative Biological Effectiveness (RBE) of kilovoltage photon beams has been previously investigated in vitro and in silico using analytical methods. The estimated values range from 1.03 to 1.82 depending on the methodology and beam energies examined.

PURPOSE

The focus of this work was to independently estimate RBE values for a range of clinically used kilovoltage beams (70-200 kVp) while investigating the suitability of using TOPAS-nBio for this task.

METHODS

Previously validated spectra of clinical beams were used to generate secondary electron spectra at several depths in a water tank phantom via TOPAS Monte Carlo (MC) simulations. Cell geometry was irradiated with the secondary electrons in TOPAS-nBio MC simulations. The deposited dose and the calculated number of DNA strand breaks were used to estimate RBE values.

RESULTS

Monoenergetic secondary electron simulations revealed the highest direct and indirect double strand break yield at approximately 20 keV. The average RBE value for the kilovoltage beams was calculated to be 1.14.

CONCLUSIONS

TOPAS-nBio was successfully used to estimate the RBE values for a range of clinical radiotherapy beams. The calculated value was in agreement with previous estimates, providing confidence in its clinical use in the future.

Estimation of relative biological effectiveness of 225Ac compared to 177Lu during [225Ac]Ac-PSMA and [177Lu]Lu-PSMA radiopharmaceutical therapy using TOPAS/TOPAS-nBio/MEDRAS

AIM

Over recent years, [225Ac]Ac-PSMA and [177Lu]Lu-PSMA radiopharmaceutical therapy have evolved as a promising treatment option for advanced prostate cancer. Especially for alpha particle emitter treatments, there is still a need for improving dosimetry, which requires accurate values of relative biological effectiveness (RBE). To achieve that, consideration of DNA damages in the cell nucleus and knowledge of the energy deposition in the location of the DNA at the nanometer scale are required. Monte Carlo particle track structure simulations provide access to interactions at this level. The aim of this study was to estimate the RBE of 225Ac compared to 177Lu. The initial damage distribution after radionuclide decay and the residual damage after DNA repair were considered.

METHODS

This study employed the TOol for PArtcile Simulation (TOPAS) based on the Geant4 simulation toolkit. Simulation of the nuclear DNA and damage scoring were performed using the TOPAS-nBio extension of TOPAS. DNA repair was modeled utilizing the Python-based program MEDRAS (Mechanistic DNA Repair and Survival). Five different cell geometries of equal volume and two radionuclide internalization assumptions as well as two cell arrangement scenarios were investigated. The radionuclide activity (number of source points) was adopted based on SPECT images of patients undergoing the above-mentioned therapies.

RESULTS

Based on the simulated dose-effect curves, the RBE of 225Ac compared to 177Lu was determined in a wide range of absorbed doses to the nucleus. In the case of spherical geometry, 3D cell arrangement and full radionuclide internalization, the RBE based on the initial damage had a constant value of approximately 2.14. Accounting for damage repair resulted in RBE values ranging between 9.38 and 1.46 for 225Ac absorbed doses to the nucleus between 0 and 50 Gy, respectively.

CONCLUSION

In this work, the consideration of DNA repair of the damage from [225Ac]Ac-PSMA and [177Lu]Lu-PSMA revealed a dose dependency of the RBE. Hence, this work suggested that DNA repair is an important aspect to understand response to different radiation qualities.

Biophysical characterization of collimated and uncollimated fields in pencil beam scanning proton therapy

Objective. The lateral dose fall-off in proton pencil beam scanning (PBS) technique remains the preferred choice for sparing adjacent organs at risk as opposed to the distal edge due to the proton range uncertainties and potentially high relative biological effectiveness. However, because of the substantial spot size along with the scattering in the air and in the patient, the lateral penumbra in PBS can be degraded. Combining PBS with an aperture can result in a sharper dose fall-off, particularly for shallow targets.Approach. The aim of this work was to characterize the radiation fields produced by collimated and uncollimated 100 and 140 MeV proton beams, using Monte Carlo simulations and measurements with a MiniPIX-Timepix detector. The dose and the linear energy transfer (LET) were then coupled with publishedin silicobiophysical models to elucidate the potential biological effects of collimated and uncollimated fields.Main results. Combining an aperture with PBS reduced the absorbed dose in the lateral fall-off and out-of-field by 60%. However, the results also showed that the absolute frequency-averaged LET (LET_F) values increased by a maximum of 3.5 keVμm^-1in collimated relative to uncollimated fields, while the dose-averaged LET (LET_D) increased by a maximum of 7 keVμm^-1. Despite the higher LET values produced by collimated fields, the predicted DNA damage yields remained lower, owing to the large dose reduction.Significance. This work demonstrated the dosimetric advantages of combining an aperture with PBS coupled with lower DNA damage induction. A methodology for calculating dose in water derived from measurements with a silicon-based detector was also presented. This work is the first to demonstrate experimentally the increase in LET caused by combining PBS with aperture, and to assess the potential DNA damage which is the initial step in the cascade of events leading to the majority of radiation-induced biological effects.

Monte Carlo simulation of physical dose enhancement in core-shell magnetic gold nanoparticles with TOPAS

The application of metal nanoparticles (MNPs) as sensitization materials is a common strategy that is used to study dose enhancement in radiotherapy. Recent in vitro tests have revealed that magnetic gold nanoparticles (NPs) can be used in cancer therapy under a magnetic field to enhance the synergistic efficiency in radiotherapy and photothermal therapy. However, magnetic gold NPs have rarely been studied as sensitization materials. In this study, we obtained further results of the sensitization properties of the magnetic gold NPs (Fe3O4@AuNPs) with or without magnetic field using the TOPAS-nBio Monte Carlo (MC) toolkit. We analyzed the properties of Fe3O4@AuNP in a single NP model and in a cell model under monoenergetic photons and brachytherapy, and we investigated whether the magnetic field contributes to the physical sensitization process. Our results revealed that the dose enhancement factor (DEF) of Fe3O4@AuNPs was lower than that of gold nanoparticles (AuNPs) in a single NP and in a cell irradiated by monoenergetic photons. But it’s worth mentioning that under a magnetic field, the DEF of targeted Fe3O4@AuNPs in a cell model with a clinical brachytherapy source was 22.17% (cytoplasm) and 6.89% (nucleus) higher than those of AuNPs (50 mg/mL). The Fe3O4@AuNPs were proved as an effective sensitization materials when combined with the magnetic field in MC simulation for the first time, which contributes to the research on in vitro tests on radiosensitization as well as clinical research in future.

Contemporary biomedical engineering perspective on volitional evolution for human radiotolerance enhancement beyond low-earth orbit

A primary objective of the National Aeronautics and Space Administration (NASA) is expansion of humankind's presence outside low-Earth orbit, culminating in permanent interplanetary travel and habitation. Having no inherent means of physiological detection or protection against ionizing radiation, humans incur capricious risk when journeying beyond low-Earth orbit for long periods. NASA has made large investments to analyze pathologies from space radiation exposure, emphasizing the importance of characterizing radiation's physiological effects. Because natural evolution would require many generations to confer resistance against space radiation, immediately pragmatic approaches should be considered. Volitional evolution, defined as humans steering their own heredity, may inevitably retrofit the genome to mitigate resultant pathologies from space radiation exposure. Recently, uniquely radioprotective genes have been identified, conferring local or systemic radiotolerance when overexpressed in vitro and in vivo. Aiding in this process, the CRISPR/Cas9 technique is an inexpensive and reproducible instrument capable of making limited additions and deletions to the genome. Although cohorts can be identified and engineered to protect against radiation, alternative and supplemental strategies should be seriously considered. Advanced propulsion and mild synthetic torpor are perhaps the most likely to be integrated. Interfacing artificial intelligence with genetic engineering using predefined boundary conditions may enable the computational modeling of otherwise overly complex biological networks. The ethical context and boundaries of introducing genetically pioneered humans are considered.

The status of medical physics in radiotherapy in China

PURPOSE

To present an overview of the status of medical physics in radiotherapy in China, including facilities and devices, occupation, education, research, etc. MATERIALS AND METHODS: The information about medical physics in clinics was obtained from the 9-th nationwide survey conducted by the China Society for Radiation Oncology in 2019. The data of medical physics in education and research was collected from the publications of the official and professional organizations.

RESULTS

By 2019, there were 1463 hospitals or institutes registered to practice radiotherapy and the number of accelerators per million population was 1.5. There were 4172 medical physicists working in clinics of radiation oncology. The ratio between the numbers of radiation oncologists and medical physicists is 3.51. Approximately, 95% of medical physicists have an undergraduate or graduate degrees in nuclear physics and biomedical engineering. 86% of medical physicists have certificates issued by the Chinese Society of Medical Physics. There has been a fast growth of publications by authors from mainland of China in the top international medical physics and radiotherapy journals since 2018.

CONCLUSIONS

Demand for medical physicists in radiotherapy increased quickly in the past decade. The distribution of radiotherapy facilities in China became more balanced. High quality continuing education and training programs for medical physicists are deficient in most areas. The role of medical physicists in the clinic has not been clearly defined and their contributions have not been fully recognized by the community.

Monte Carlo study of TG-43 dosimetry parameters of GammaMed Plus high dose rate 192 Ir brachytherapy source using TOPAS

PURPOSE

To develop a simulation model for GammaMed Plus high dose rate ¹⁹² Ir brachytherapy source in TOPAS Monte Carlo software and validate it by calculating the TG-43 dosimetry parameters and comparing them with published data.

METHODS

We built a model for GammaMed Plus high dose rate brachytherapy source in TOPAS. The TG-43 dosimetry parameters including air-kerma strength S_K , dose-rate constant Λ, radial dose function g_L (r), and 2D anisotropy function F(r,θ) were calculated using Monte Carlo simulation with Geant4 physics models and NNDC ¹⁹² Ir spectrum. Calculations using an old ¹⁹² Ir spectrum were also carried out to evaluate the impact of incident spectrum and cross sections. The results were compared with published data.

RESULTS

For calculations using the NNDC spectrum, the air-kerma strength per unit source activity S_K /A and Λ were 1.0139 × 10^-7 U/Bq and 1.1101 cGy.h^-1 .U^-1 , which were 3.56% higher and 0.62% lower than the reference values, respectively. The g_L (r) agreed with reference values within 1% for radial distances from 2 mm to 20 cm. For radial distances of 1, 3, 5, and 10 cm, the agreements between F(r,θ) from this work and the reference data were within 1.5% for 15° < θ < 165°, and within 4% for all θ values. The discrepancies were attributed to the updated source spectrum and cross sections. They caused deviations of the S_K /A of 2.90% and 0.64%, respectively. As for g_L (r), they caused average deviations of -0.22% and 0.48%, respectively. Their impact on F(r,θ) was not quantified for the relatively high statistical uncertainties, but basically they did not result in significant discrepancies.

CONCLUSION

A model for GammaMed Plus high dose rate ¹⁹² Ir brachytherapy source was developed in TOPAS and validated following TG-43 protocols, which can be used for future studies. The impact of updated incident spectrum and cross sections on the dosimetry parameters was quantified.