Purpose: This study aims to develop a comprehensive simulation framework to connect radiation effects from the microscopic to the nanoscopic scale. Method: The process begins with a Geant4-DNA simulation based on the example ”molecularDNA”, producing a dataset of twelve different types of early DNA damages within an Escherichia coli (E. coli) bacterium, generated by proton irradiation at different kinetic energies, giving a nano-scale view of the particle–matter interaction. Then we pass to the micro-scale with a Geant4 simulation, based on the example ”radiobiology”, providing a microscopic view of proton interactions with matter through the Linear Energy Transfer (LET). Then GANDALF (Generative ANsatz for DNA damage evALuation and Forecast) Machine Learning (ML) toolkit, a Neural Network (NN)-based regression system, is employed to correlate the micro-scale LET data with the nano-scale occurrences of DNA damages in the E. coli bacterium. Results: The trained ML algorithm provides a practical tool to convert LET curves versus depth in a water phantom into DNA damage curves for twelve distinct types of DNA damage. To assess the performance, we evaluated the choice and optimization of the regression system based on its interpolation and extrapolation capabilities, ensuring the model could reliably predict DNA damage under various conditions. Conclusions: Through the synergistic integration of Geant4, Geant4-DNA and ML, the study provides a tool to easily convert the results at the micro-scale of Geant4 to those at the nano-scale of Geant4-DNA without having to deal with the high CPU time requirements of the latter.

The early DNA damage on the surface of the Moon due to GCR protons and alpha particles were assessed using a multiscale approach in Geant4. This consisted of three simulation stages. A periodic boundary conditions approach was used to obtain the radiation field on the surface and inside a proposed lunar habitat. The radiation field on the cellular scale was obtained in the organs of male and female astronauts using the ICRP145 tetrahedral mesh phantoms. This was subsequently simulated using a full human cell model in Geant4-DNA to obtain the early DNA damage. Geant4-DNA track structure ionisation models upper energy limits were extended to be able to model the interactions of the GCR at sub-cellular level, covering an energy range from a few eV up to 1 TeV. Hadronic interactions and the modelling of induced radiochemical species were also implemented. The early DNA damage was assessed using the Geant4-DNA molecularDNA example. A greater yield of DNA damage was observed on the lunar surface compared with the habitat, and indirect damage due to induced hydroxyl radicals constituted most of the damage. This study demonstrates a complete simulation pipeline for the assessment of early DNA damage in astronauts in the space radiation environment.

Purpose: Based on considerable interest to enlarge the experimental database of radioresistant cells after their irradiation with helium ions, HTB140, MCF-7 and HTB177 human malignant cells are exposed to helium ion beams having different linear energy transfer (LET). Materials and methods: The cells are irradiated along the widened 62 MeV/u helium ion Bragg peak, providing LET of 4.9, 9.8, 23.4 and 36.8 keV/µm. Numerical simulations with the Geant4 toolkit are used for the experimental design. Cell survival is evaluated and compared with reference γ-rays. DNA double strand breaks are assessed via γ-H2AX foci. Results: With the increase of LET, surviving fractions at 2 Gy decrease, while RBE (2 Gy, γ) gradually increase. For HTB140 cells, above the dose of 4 Gy, a slight saturation of survival is observed while the increase of RBE (2 Gy, γ) remains unaffected. With the increase of LET the increase of γ-H2AX foci is revealed at 0.5 h after irradiation. There is no significant difference in the number of foci between the cell lines for the same LET. From 0.5 to 24 h, the number of foci drops reaching its residual level. For each time point, there are small differences in DNA DSB among the three cell lines. Conclusion: Analyses of data acquired for the three cell lines irradiated by helium ions, having different LET, reveal high elimination capacity and creation of a large number of DNA DSB with respect to γ-rays, and are between those reported for protons and carbon ions.

The authors regret having incorrectly written in equation (7) that the fraction of surviving cells (SF) is obtained using the natural logarithm of the negative value of the number of lethal lesions Lf. The natural logarithm must be replaced by the exponential function as follows: [Formula presented] Although the equation is incorrectly written in the paper, its correct expression has been used in the results of the analysis shown in Figure 4. This error has no impact on the other results and conclusions of the paper. The authors would like to apologise for any inconvenience caused.

Background: This study aimed to develop a novel human cell geometry for the Geant4-DNA simulation toolkit that explicitly incorporates all 23 chromosome pairs of the human cell. This approach contrasts with the existing, default human cell, geometrical model, which utilizes a continuous Hilbert curve. Methods: A Python-based tool named “complexDNA” was developed to facilitate the design of both simple and complex DNA geometries. This tool was employed to construct a human cell geometry with individual pairs of chromosomes. Subsequently, the performance of this chromosomal model was compared to the standard human cell model provided in the “molecularDNA” Geant4-DNA example. Results: Simulations using the new chromosomal model revealed minimal discrepancies in DNA damage yield and fragment size distribution compared to the default human cell model. Notably, the chromosomal model demonstrated significant computational efficiency, requiring approximately three times less simulation time to achieve equivalent results. Conclusions: This work highlights the importance of incorporating chromosomal structure into human cell models for radiation biology research. The “complexDNA” tool offers a valuable resource for creating intricate DNA structures for future studies. Further refinements, such as implementing smaller voxels for euchromatin regions, are proposed to enhance the model's accuracy.

Purpose: Interdisciplinary scientific communities have shown large interest to achieve a mechanistic description of radiation-induced biological damage, aiming to predict biological results produced by different radiation quality exposures. Monte Carlo track-structure simulations are suitable and reliable for the study of early DNA damage induction used as input for assessing DNA damage. This study presents the most recent improvements of a Geant4-DNA simulation tool named “dsbandrepair”. Methods: “dsbandrepair” is a Monte Carlo simulation tool based on a previous code (FullSim) that estimates the induction of early DNA single-strand breaks (SSBs) and double-strand breaks (DSBs). It uses DNA geometries generated by the DNAFabric computational tool for simulating the induction of early single-strand breaks (SSBs) and double-strand breaks (DSBs). Moreover, the new tool includes some published radiobiological models for survival fraction and un-rejoined DSB. Its application for a human fibroblast cell and human umbilical vein endothelial cell containing both heterochromatin and euchromatin was conducted. In addition, this new version offers the possibility of using the new IRT-syn method for computing the chemical stage. Results: The direct and indirect strand breaks, SSBs, DSBs, and damage complexity obtained in this work are equivalent to those obtained with the previously published simulation tool when using the same configuration in the physical and chemical stages. Simulation results on survival fraction and un-rejoined DSB are in reasonable agreement with experimental data. Conclusions: “dsbandrepair” is a tool for simulating DNA damage and repair, benchmarked against experimental data. It has been released as an advanced example in Geant4.11.2.

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.

Purpose: Track structure Monte Carlo (MC) codes have achieved successful outcomes in the quantitative investigation of radiation-induced initial DNA damage. The aim of the present study is to extend a Geant4-DNA radiobiological application by incorporating a feature allowing for the prediction of DNA rejoining kinetics and corresponding cell surviving fraction along time after irradiation, for a Chinese hamster V79 cell line, which is one of the most popular and widely investigated cell lines in radiobiology. Methods: We implemented the Two-Lesion Kinetics (TLK) model, originally proposed by Stewart, which allows for simulations to calculate residual DNA damage and surviving fraction along time via the number of initial DNA damage and its complexity as inputs. Results: By optimizing the model parameters of the TLK model in accordance to the experimental data on V79, we were able to predict both DNA rejoining kinetics at low linear energy transfers (LET) and cell surviving fraction. Conclusion: This is the first study to demonstrate the implementation of both the cell surviving fraction and the DNA rejoining kinetics with the estimated initial DNA damage, in a realistic cell geometrical model simulated by full track structure MC simulations at DNA level and for various LET. These simulation and model make the link between mechanistic physical/chemical damage processes and these two specific biological endpoints.

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.

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.