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 Geant4-DNA low energy extension of the Geant4 Monte Carlo (MC) toolkit is a continuously evolving MC simulation code permitting mechanistic studies of cellular radiobiological effects. Geant4-DNA considers the physical, chemical, and biological stages of the action of ionizing radiation (in the form of x- and γ-ray photons, electrons and β^±-rays, hadrons, α-particles, and a set of heavier ions) in living cells towards a variety of applications ranging from predicting radiotherapy outcomes to radiation protection both on earth and in space. In this work, we provide a brief, yet concise, overview of the progress that has been achieved so far concerning the different physical, physicochemical, chemical, and biological models implemented into Geant4-DNA, highlighting the latest developments. Specifically, the “dnadamage1” and “molecularDNA” applications which enable, for the first time within an open-source platform, quantitative predictions of early DNA damage in terms of single-strand-breaks (SSBs), double-strand-breaks (DSBs), and more complex clustered lesions for different DNA structures ranging from the nucleotide level to the entire genome. These developments are critically presented and discussed along with key benchmarking results. The Geant4-DNA toolkit, through its different set of models and functionalities, offers unique capabilities for elucidating the problem of radiation quality or the relative biological effectiveness (RBE) of different ionizing radiations which underlines nearly the whole spectrum of radiotherapeutic modalities, from external high-energy hadron beams to internal low-energy gamma and beta emitters that are used in brachytherapy sources and radiopharmaceuticals, respectively.

Accurately modeling the radiobiological mechanisms responsible for the induction of DNA damage remains a major scientific challenge, particularly for understanding the effects of low doses of ionizing radiation on living beings, such as the induction of carcinogenesis. A computational approach based on the Monte Carlo technique to simulate track structures in a biological medium is currently the most reliable method for calculating the early effects induced by ionizing radiation on DNA, the primary cellular target of such effects. The Geant4-DNA Monte Carlo toolkit can simulate not only the physical, but also the physico-chemical and chemical stages of water radiolysis. These stages can be combined with simplified geometric models of biological targets, such as DNA, to assess direct and indirect early DNA damage. In this study, DNA damage induced in a human fibroblast cell was evaluated using Geant4-DNA as a function of incident particle type (gammas, protons, and alphas) and energy. The resulting double-strand break yields as a function of linear energy transfer closely reproduced recent experimental data. Other quantities, such as fragment length distribution, scavengeable damage fraction, and time evolution of damage within an analytical repair model also supported the plausibility of predicting DNA damage using Geant4-DNA.The complete simulation chain application “molecularDNA”, an example for users of Geant4-DNA, will soon be distributed through Geant4.

The ROSIRIS project aims to study the radiobiology of integrated systems for medical treatment optimisation using ionising radiations and evaluate the associated risk. In the framework of this project, one research focus is the interpretation of the initial radio-induced damage in DNA created by ionising radiation (and detected by γ-H2AX foci analysis) from the track structure of the incident particles. In order to calculate the track structure of ionising particles at a nanometric level, the Geant4 Monte Carlo toolkit was used. Geant4 (Object Oriented Programming Architecture in C++) offers a common platform, available free to all users and relatively easy to use. Nevertheless, the current low-energy threshold for electromagnetic processes in GEANT4 is set to 1 keV (250 eV using the Livermore processes), which is an unsuitable value for nanometric applications. To lower this energy threshold, the necessary interaction processes and models were identified, and the corresponding available cross sections collected from the literature. They are mostly based on the plane-wave Born approximation (first Born approximation, or FBA) for inelastic interactions and on semi-empirical models for energies where the FBA fails (at low energies). In this paper, the extensions that have been introduced into the 9.3 release of the Geant4 toolkit are described, the so-called Geant4-DNA extension, including a set of processes and models adapted in this study and permitting the simulation of electron (8 eV-1 MeV), proton (100 eV-100 MeV) and alpha particle (1 keV-10 MeV) interactions in liquid water.

High-Z metallic nanoparticles (NPs) are new players in the therapeutic arsenal against cancer, especially radioresistant cells. Indeed, the presence of these NPs inside malignant cells is believed to enhance the effect of ionizing radiation by locally increasing the dose deposition. In this context, the potential of platinum nanoparticles (PtNPs) as radiosensitizers was investigated in two breast cancer cell lines, T47D and MDA-MB-231, showing a different radiation sensitivity. PtNPs were internalized in the two cell lines and localized in lysosomes and multivesicular bodies. Analyses of cell responses in terms of clonogenicity, survival, mortality, cell-cycle distribution, oxidative stress, and DNA double-strand breaks did not reveal any significant enhancement effect when cells were pre-exposed to PtNPs before being irradiated, as compared to radiation alone. This result is different from that reported in a previous study performed, under the same conditions, on cervical cancer HeLa cells. This shows that the efficacy of radio-enhancement is strongly cell-type-dependent. Simulation of the early stage ionization processes, taking into account the irradiation characteristics and realistic physical parameters in the biological sample, indicated that PtNPs could weakly increase the dose deposition (by 3%) in the immediate vicinity of the nanoparticles. Some features that are potentially responsible for the biological effect could not be taken into account in the simulation. Thus, chemical and biological effects could explain this discrepancy. For instance, we showed that, in these breast cancer cell lines, PtNPs exhibited ambivalent redox properties, with an antioxidant potential which could counteract the radio-enhancement effect. This work shows that the efficacy of PtNPs for enhancing radiation effects is strongly cell-dependent and that no effect is observed in the case of the breast cancer cell lines T47D and MDA-MB-231. Thus, more extensive experiments using other relevant biological models are needed in order to evaluate such combined strategies, since several clinical trials have already demonstrated the success of combining nanoagents with radiotherapy in the treatment of a range of tumor types.

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 (10⁹) 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.