This Special Report presents a description of Geant4-DNA user applications dedicated to the simulation of track structures (TS) in liquid water and associated physical quantities (e.g., range, stopping power, mean free path…). These example applications are included in the Geant4 Monte Carlo toolkit and are available in open access. Each application is described and comparisons to recent international recommendations are shown (e.g., ICRU, MIRD), when available. The influence of physics models available in Geant4-DNA for the simulation of electron interactions in liquid water is discussed. Thanks to these applications, the authors show that the most recent sets of physics models available in Geant4-DNA (the so-called "option4" and "option 6" sets) enable more accurate simulation of stopping powers, dose point kernels, and W-values in liquid water, than the default set of models ("option 2") initially provided in Geant4-DNA. They also serve as reference applications for Geant4-DNA users interested in TS simulations.

In order to improve the understanding of the mechanisms involved in the generation of early DNA damage, a new calculation chain based on the Geant4-DNA toolkit was developed. This work presents for the first time the simulation of the physical, physicochemical and chemical stages of early radiation damage at the scale of an entire human genome (fibroblast, male) and using Geant4-DNA models. The DnaFabric software was extended to generate and export this nucleus model to a text file with a specific format that can be read by Geant4 user applications. This calculation chain was used to simulate the irradiation of the nucleus by primary protons of different energies (0,5; 0,7; 0,8; 1; 1,5; 2; 3; 4; 5; 10; 20 MeV) and the results, in terms of DNA double strand breaks, agree with experimental data found in the literature (pulsed field electrophoresis technique). These results show that the simulation is consistent and that its parameters are well balanced. Among the different parameters that can be adjusted, our results demonstrate that the criterion used to select direct strand break appears to have a very significant role on the final number of simulated double strand breaks.

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%.

Purpose In this work, we present simulated double‐strand breaks (DSBs) obtained for two human cell nucleus geometries. The first cell nucleus represents fibroblasts, filled with DNA molecules in different compaction forms: heterochromatin or euchromatin only. The second one represents an endothelial cell nucleus, either filled with heterochromatin only or with a uniform distribution of 48% of heterochromatin and 52% of euchromatin, obtained from measurements carried out at IRSN. Protons and alpha particles of different energies were used as projectiles. Each cell nucleus model includes a multi‐scale description of the DNA target from the molecular level to the whole human genome representation. Methods The cell nucleus models were generated using an extended version of the DnaFabric software in which a new model of euchromatin was implemented in addition to the existing model of heterochromatin. Thus, each nucleus model contains the complete human genome (a total of 6 Gbp) in the G0/G1 phase of the cycle, filled with a continuous chromatin fiber per chromosome that can take into account the heterochromatin and the euchromatin compaction. These geometries were then exported to a simulation chain using the Monte Carlo toolkit Geant4‐DNA to perform computations of the physical, physicochemical, and chemical stages, in order to evaluate the influence of chromatin compaction on DSB induction and the contribution of direct and indirect damage, as well as DSB complexity. Results More direct damage and less indirect damage were observed in the heterochromatin than in the euchromatin. Nevertheless, no difference in terms of DSB complexity was observed between those formed in the heterochromatin or the euchromatin models. Yields of DSB/Gy/Gbp show an increase when both heterochromatin and euchromatin models are taken into account, compared to when only heterochromatin is considered. Conclusions The results presented indicate that the chromatin compaction decreases DNA damage generated by ionizing radiation and thus, DNA compaction should be considered for the simulation of DNA repair and other cellular outcomes.

The objective of this work was to study the differences in terms of early biological effects that might exist between different X-rays energies by using a mechanistic approach. To this end, radiobiological experiments exposing cell monolayers to three X-ray energies were performed in order to assess the yields of early DNA damage, in particular of double-strand breaks (DSBs). The simulation of these irradiations was set in order to understand the differences in the obtained experimental results. Hence, simulated results in terms of microdosimetric spectra and early DSB induction were analyzed and compared to the experimental data. Human umbilical vein endothelial cells (HUVECs) were irradiated with 40, 220 kVp, and 4 MV X-rays. The Geant4 Monte Carlo simulation toolkit and its extension Geant4-DNA were used for the simulations. Microdosimetric calculations aiming to determine possible differences in the variability of the energy absorbed by the irradiated cell population for those photon spectra were performed on 10,000 endothelial cell nuclei representing a cell monolayer. Nanodosimetric simulations were also carried out using a computation chain that allowed the simulation of physical, physico-chemical, and chemical stages on a single realistic endothelial cell nucleus model including both heterochromatin and euchromatin. DNA damage was scored in terms of yields of prompt DSBs per Gray (Gy) and per giga (109) base pair (Gbp) and DSB complexity was derived in order to be compared to experimental data expressed as numbers of histone variant H2AX (γ-H2AX) foci per cell. The calculated microdosimetric spread in the irradiated cell population was similar when comparing between 40 and 220 kVp X-rays and higher when comparing with 4 MV X-rays. Simulated yields of induced DSB/Gy/Gbp were found to be equivalent to those for 40 and 220 kVp but larger than those for 4 MV, resulting in a relative biological effectiveness (RBE) of 1.3. Additionally, DSB complexity was similar between the considered photon spectra. Simulated results were in good agreement with experimental data obtained by IRSN (Institut de radioprotection et de sûreté nucléaire) radiobiologists. Despite differences in photon energy, few differences were observed when comparing between 40 and 220 kVp X-rays in microdosimetric and nanodosimetric calculations. Nevertheless, variations were observed when comparing between 40/220 kVp and 4 MV X-rays. Thanks to the simulation results, these variations were able to be explained by the differences in the production of secondary electrons with energies below 10 keV.