Purpose:
This study aimed to develop a computational environment for the accurate simulation of human cancer cell irradiation using Geant4-DNA. New cell geometrical models were developed and irradiated by alpha particle beams to induce DNA damage. The proposed approach may help further investigation of the benefits of external alpha irradiation therapy.

Methods:
The Geant4-DNA Monte Carlo (MC) toolkit allows the simulation of cancer cell geometries that can be combined with accurate modelling of physical, physicochemical and chemical stages of liquid water irradiation, including radiolytic processes. Geant4-DNA is used to calculate direct and non-direct DNA damage yields, such as single and double strand breaks, produced by the deposition of energy or by the interaction of DNA with free radicals.

Results:
In this study, the “molecularDNA” example application of Geant4-DNA was used to quantify early DNA damage in human cancer cells upon irradiation with alpha particle beams, as a function of linear energy transfer (LET). The MC simulation results are compared to experimental data, as well as previously published simulation data. The simulation results agree well with the experimental data on DSB yields in the lower LET range, while the experimental data on DSB yields are lower than the results obtained with the “molecularDNA” example in the higher LET range.

Conclusion:
This study explored and demonstrated the possibilities of the Geant4-DNA toolkit together with the “molecularDNA” example to simulate the helium beam irradiation of cancer cell lines, to quantify the early DNA damage, or even the following DNA damage response.

Purpose
The scientific community shows great interest in the study of DNA damage induction, DNA damage repair, and the biological effects on cells and cellular systems after exposure to ionizing radiation. Several in silico methods have been proposed so far to study these mechanisms using Monte Carlo simulations. This study outlines a Geant4-DNA example application, named “molecularDNA”, publicly released in the 11.1 version of Geant4 (December 2022).

Methods
It was developed for novice Geant4 users and requires only a basic understanding of scripting languages to get started. The example includes two different DNA-scale geometries of biological targets, namely “cylinders” and “human cell”. This public version is based on a previous prototype and includes new features, such as: the adoption of a new approach for the modeling of the chemical stage, the use of the standard DNA damage format to describe radiation-induced DNA damage, and upgraded computational tools to estimate DNA damage response.

Results
Simulation data in terms of single-strand break and double-strand break yields were produced using each of the available geometries. The results were compared with the literature, to validate the example, producing less than 5% difference in all cases. Conclusion: “molecularDNA” is a prototype tool that can be applied in a wide variety of radiobiology studies, providing the scientific community with an open-access base for DNA damage quantification calculations. New DNA and cell geometries for the “molecularDNA” example will be included in future versions of Geant4-DNA.

The advancement of multidisciplinary research fields dealing with ionising radiation induced biological damage – radiobiology, radiation physics, radiation protection and, in particular, medical physics – requires a clear mechanistic understanding of how cellular damage is induced by ionising radiation. Monte Carlo (MC) simulations provide a promising approach for the mechanistic simulation of radiation transport and radiation chemistry, towards the in silico simulation of early biological damage. We have recently developed a fully integrated MC simulation that calculates early single strand breaks (SSBs) and double strand breaks (DSBs) in a fractal chromatin based human cell nucleus model. The results of this simulation are almost equivalent to past MC simulations when considering direct/indirect strand break fraction, DSB yields and fragment distribution. The simulation results agree with experimental data on DSB yields within 13.6% on average and fragment distributions agree within an average of 34.8%.

We extended a generic Geant4 application for mechanistic DNA damage simulations to an Escherichia coli cell geometry, finding electron damage yields and proton damage yields largely in line with experimental results. Depending on the simulation of radical scavenging, electrons double strand breaks (DSBs) yields range from 0.004 to 0.010 DSB Gy⁻¹ Mbp⁻¹, while protons have yields ranging from 0.004 DSB Gy⁻¹ Mbp⁻¹ at low LETs and with strict assumptions concerning scavenging, up to 0.020 DSB Gy⁻¹ Mbp⁻¹ at high LETs and when scavenging is weakest. Mechanistic DNA damage simulations can provide important limits on the extent to which physical processes can impact biology in low background experiments. We demonstrate the utility of these studies for low dose radiation biology calculating that in E. coli, the median rate at which the radiation background induces double strand breaks is 2.8 × 10⁻⁸ DSB day⁻¹, significantly less than the mutation rate per generation measured in E. coli, which is on the order of 10⁻³.

Mechanistic modelling of DNA damage in Monte Carlo simulations is highly sensitive to the parameters that define DNA damage. In this work, we use a simple testing geometry to investigate how different choices of physics models and damage model parameters can change the estimation of DNA damage in a mechanistic DNA damage simulation built in Geant4-DNA. The choice of physics model can lead to variations by up to a factor of two in the yield of physically induced strand breaks, and the parameters that determine scavenging, and physical and chemical single strand break induction can have even larger consequences. Using low energy electrons as primary particles, a variety of parameters are tested in this geometry in order to arrive at a parameter set consistent with past simulation studies. We find that the modelling of scavenging can play an important role in determining results, and speculate that high-scavenging regimes, where only chemical radicals within 1 nm of DNA are simulated, could provide a good means of testing mechanistic DNA simulations.

A novel treatment planning framework, the Relative Biological Effective Dose (RBED), for high Z nanoparticle (NP)-enhanced photon radiotherapy is developed and tested in silico for the medical exemplar of neoadjuvant (preoperative) breast cancer MV photon radiotherapy. Two different treatment scenarios, conventional and high Z NP enhanced, were explored with a custom Geant4 application that was developed to emulate the administration of a single 2 Gy fraction as part of a 50 Gy radiotherapy treatment plan. It was illustrated that there was less than a 1% difference in the dose deposition throughout the standard and high Z NP-doped adult female phantom. Application of the RBED framework found that the extent of possible biological response with high Z NP doping was great than expected via the dose deposition alone. It is anticipated that this framework will assist the scientific community in future high Z NP-enhanced in-silico, pre-clinical and clinical trials.