"dsbandrepair" - An updated Geant4-DNA simulation tool for evaluating the radiation-induced DNA damage and its repair.

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.

Geant4-DNA simulation of human cancer cells irradiation with helium ion beams.

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.

A framework for prediction of personalized pediatric nuclear medical dosimetry based on machine learning and Monte Carlo techniques.

Objective:A methodology is introduced for the development of an internal dosimetry prediction toolkit for nuclear medical pediatric applications. The proposed study exploits Artificial Intelligence techniques using Monte Carlo simulations as ground truth for accurate prediction of absorbed doses per organ prior to the imaging acquisition considering only personalized anatomical characteristics of any new pediatric patient.Approach:GATE Monte Carlo simulations were performed using a population of computational pediatric models to calculate the specific absorbed dose rates (SADRs) in several organs. A simulated dosimetry database was developed for 28 pediatric phantoms (age range 2-17 years old, both genders) and 5 different radiopharmaceuticals. Machine Learning regression models were trained on the produced simulated dataset, with leave one out cross validation for the prediction model evaluation. Hyperparameter optimization and ensemble learning techniques for a variation of input features were applied for achieving the best predictive power, leading to the development of a SADR prediction toolkit for any new pediatric patient for the studied organs and radiopharmaceuticals.Main results. SADR values for 30 organs of interest were calculated via Monte Carlo simulations for 28 pediatric phantoms for the cases of five radiopharmaceuticals. The relative percentage uncertainty in the extracted dose values per organ was lower than 2.7%. An internal dosimetry prediction toolkit which can accurately predict SADRs in 30 organs for five different radiopharmaceuticals, with mean absolute percentage error on the level of 8% was developed, with specific focus on pediatric patients, by using Machine Learning regression algorithms, Single or Multiple organ training and Artificial Intelligence ensemble techniques. Significance: A large simulated dosimetry database was developed and utilized for the training of Machine Learning models. The developed predictive models provide very fast results (<2 s) with an accuracy >90% with respect to the ground truth of Monte Carlo, considering personalized anatomical characteristics and the biodistribution of each radiopharmaceutical. The proposed method is applicable to other medical dosimetry applications in different patients' populations.

Standardization and Validation of Brachytherapy Seeds' Modelling Using GATE and GGEMS Monte Carlo Toolkits.

This study aims to validate GATE and GGEMS simulation toolkits for brachytherapy applications and to provide accurate models for six commercial brachytherapy seeds, which will be freely available for research purposes. The AAPM TG-43 guidelines were used for the validation of two Low Dose Rate (LDR), three High Dose Rate (HDR), and one Pulsed Dose Rate (PDR) brachytherapy seeds. Each seed was represented as a 3D model and then simulated in GATE to produce one single Phase-Space (PHSP) per seed. To test the validity of the simulations' outcome, referenced data (provided by the TG-43) was compared with GATE results. Next, validation of the GGEMS toolkit was achieved by comparing its outcome with the GATE MC simulations, incorporating clinical data. The simulation outcomes on the radial dose function (RDF), anisotropy function (AF), and dose rate constant (DRC) for the six commercial seeds were compared with TG-43 values. The statistical uncertainty was limited to 1% for RDF, to 6% (maximum) for AF, and to 2.7% (maximum) for the DRC. GGEMS provided a good agreement with GATE when compared in different situations: (a) Homogeneous water sphere, (b) heterogeneous CT phantom, and (c) a realistic clinical case. In addition, GGEMS has the advantage of very fast simulations. For the clinical case, where TG-186 guidelines were considered, GATE required 1 h for the simulation while GGEMS needed 162 s to reach the same statistical uncertainty. This study produced accurate models and simulations of their emitted spectrum of commonly used commercial brachytherapy seeds which are freely available to the scientific community. Furthermore, GGEMS was validated as an MC GPU based tool for brachytherapy. More research is deemed necessary for the expansion of brachytherapy seed modeling.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE.

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

IDDRRA: A novel platform, based on Geant4-DNA to quantify DNA damage by ionizing radiation.

PURPOSE: This study proposes a novel computational platform that we refer to as IDDRRA (DNA Damage Response to Ionizing RAdiation), which uses Monte Carlo (MC) simulations to score radiation induced DNA damage. MC simulations provide results of high accuracy on the interaction of radiation with matter while scoring the energy deposition based on state-of-the-art physics and chemistry models and probabilistic methods. METHODS: The IDDRRA software is based on the Geant4-DNA toolkit together with new tools that were developed for the purpose of this study, including a new algorithm that was developed in Python for the design of the DNA molecules. New classes were developed in C++ to integrate the GUI and produce the simulation's output in text format. An algorithm was also developed to analyze the simulation's output in terms of energy deposition, Single Strand Breaks (SSB), Double Strand Breaks (DSB) and Cluster Damage Sites (CDS). Finally, a new tool was developed to implement probabilistic SSB and DSB repair models using MC techniques. RESULTS: This article provides the first benchmarks that the user of the IDDRRA tool can use to validate the functionality of the software as well as to provide a starting point to produce different types of DNA simulations. These benchmarks incorporate different kind of particles (e-, e+, protons, electron spectrum) and DNA molecules. CONCLUSION: We have developed the IDDRRA tool and demonstrated its use to study various aspects of the modeling and simulation of a DNA irradiation experiment. The tool is expandable and can be expanded by other users with new benchmarks and applications based on the user's needs and experience. New functionality will be added over time, including the quantification of the indirect damage.

Ionizing Radiation and Complex DNA Damage: Quantifying the Radiobiological Damage Using Monte Carlo Simulations.

Ionizing radiation is a common tool in medical procedures. Monte Carlo (MC) techniques are widely used when dosimetry is the matter of investigation. The scientific community has invested, over the last 20 years, a lot of effort into improving the knowledge of radiation biology. The present article aims to summarize the understanding of the field of DNA damage response (DDR) to ionizing radiation by providing an overview on MC simulation studies that try to explain several aspects of radiation biology. The need for accurate techniques for the quantification of DNA damage is crucial, as it becomes a clinical need to evaluate the outcome of various applications including both low- and high-energy radiation medical procedures. Understanding DNA repair processes would improve radiation therapy procedures. Monte Carlo simulations are a promising tool in radiobiology studies, as there are clear prospects for more advanced tools that could be used in multidisciplinary studies, in the fields of physics, medicine, biology and chemistry. Still, lot of effort is needed to evolve MC simulation tools and apply them in multiscale studies starting from small DNA segments and reaching a population of cells.

Quantification of DNA double-strand breaks using Geant4-DNA.

PURPOSE: This study aims to standardize the simulation procedure in measuring DNA double-strand breaks (DSBs), by using advanced Monte Carlo toolkits, and newly introduced experimental methods for DNA DSB measurement. METHODS: For the experimental quantification of DNA DSB, an innovative DNA dosimeter was used to produce experimental data. GATE in combination with Geant4-DNA toolkit were exploited to simulate the experimental environment. The PDB4DNA example of Geant4-DNA was upgraded and investigated. Parameters of the simulation such energy threshold (ET) for a strand break and base pair threshold (BPT) for a DSB were evaluated, depending on the dose. RESULTS: Simulations resulted to minimum differentiation in comparison to experimental data for ET = 19 ± 1 eV and BPT = 10 bp, and high differentiation for ET<17.5 eV or ET>22.5 eV and BPT = 10 bp. There was also small differentiation for ET = 17.5 eV and BPT = 6 bp. Uncertainty has been kept lower than 3%. CONCLUSIONS: This study includes first results on the quantification of DNA double-strand breaks. The energy spectrum of a LINAC was simulated and used for the first time to irradiate DNA molecules. Simulation outcome was validated on experimental data that were produced by a prototype DNA dosimeter.