Abstract:PurposeCalculations of deoxyribonucleic acid (DNA) damages involve many parameters in the computation process. As these parameters are often subject to uncertainties, it is of central importance to comprehensively quantify their impacts on DNA single‐strand break (SSB) and double‐strand break (DSB) yields. This has been a challenging task due to the required large number of simulations and the relatively low computational efficiency using CPU‐based MC packages. In this study, we present comprehensive evaluations on sensitivities and uncertainties of DNA SSB and DSB yields on 12 parameters using our GPU‐based MC tool, gMicroMC.MethodsWe sampled one electron at a time in a water sphere containing a human lymphocyte nucleus and transport the electrons and generated radicals until 2 Gy dose was accumulated in the nucleus. We computed DNA damages caused by electron energy deposition events in the physical stage and the hydroxyl radicals at the end of the chemical stage. We repeated the computations by varying 12 parameters: (a) physics cross section, (b) cutoff energy for electron transport, (c)–(e) three branching ratios of hydroxyl radicals in the de‐excitation of excited water molecules, (f) temporal length of the chemical stage, (g)–(h) reaction radii for direct and indirect damages, (i) threshold energy defining the threshold damage model to generate a physics damage, (j)–(k) minimum and maximum energy values defining the linear‐probability damage model to generate a physics damage, and (l) probability to generate a damage by a radical. We quantified sensitivity of SSB and DSB yields with respect to these parameters for cases with 1.0 and 4.5 keV electrons. We further estimated uncertainty of SSB and DSB yields caused by uncertainties of these parameters.ResultsUsing a threshold of 10% uncertainty as a criterion, threshold energy in the threshold damage model, maximum energy in the linear‐probability damage model, and probability for a radical to generate a damage were found to cause large uncertainties in both SSB and DSB yields. The scaling factor of the cross section, cutoff energy, physics reaction radius, and minimum energy in the linear‐probability damage model were found to generate large uncertainties in DSB yields.ConclusionsWe identified parameters that can generate large uncertainties in the calculations of SSB and DSB yields. Our study could serve as a guidance to reduce uncertainties of parameters and hence uncertainties of the simulation results.

A molecular dynamics simulation of DNA damage induction by ionizing radiation

Reactive Molecular Dynamics study on the first steps of DNA-damage by free hydroxyl radicals

A simulation approach for determining the spectrum of DNA damage induced by protons

A Geant4-DNA Evaluation of Radiation-Induced DNA Damage on a Human Fibroblast

Investigate the Radiation-Induced Damage on an Atomistic DNA Model by Using Geant4-DNA Toolkit

Simulation of DNA damage using Geant4-DNA: an overview of the "molecularDNA" example application

Monte Carlo Simulation of Ionizing Radiation Induced DNA Strand Breaks Utilizing Coarse Grained High-Order Chromatin Structures

Accessing radiation damage to biomolecules on the nanoscale by particle-scattering simulations

Electronic Excitation Response of DNA to High-energy Proton Radiation in Water

Monte Carlo Simulations of DNA Damage and Cellular Response to Hadron Irradiation

A mechanistic simulation of induced DNA damage in a bacterial cell by X- and gamma rays: a parameter study

Cellular Response to Proton Irradiation: A Simulation Study with TOPAS-nBio

Simulation and Experimental Verification of DNA Damage Due to X-Rays Interaction

Multiscale QM/MM Molecular Dynamics Study on the First Steps of Guanine-Damage by Free Hydroxyl Radicals in Solution

Simulation of Radiation-Induced DNA Damage and Protection by Histones Using the Code RITRACKS

Modeling the effect of ion-induced shock waves and DNA breakage with the reactive CHARMM force field

Application of a simple DNA damage model developed for electrons to proton irradiation

A closed parameterization of DNA-damage by charged particles as a function of energy - A geometrical approach

A Parameter Sensitivity Study For Simulating Dna Damage After Proton Irradiation Using Topas-Nbio

DNA Strand Breaks Induced by Electrons Simulated with Nanodosimetry Monte Carlo Simulation Code: NASIC

A New Open-Source Gpu-Based Microscopic Monte Carlo Simulation Tool for the Calculations of Dna Damages Caused by Ionizing Radiation - Part Ii: Sensitivity and Uncertainty Analysis