Abstract:We present a multi-scale simulation of the early stage of DNA damages by the indirect action of hydroxyl ((•)OH) free radicals generated by electrons and protons. The computational method comprises of interfacing the Geant4-DNA Monte Carlo with ReaxFF molecular dynamics software. A clustering method was employed to map the coordinates of (•)OH-radicals extracted from the ionization-track-structures onto nano-meter simulation voxels filled with DNA and water molecules. The molecular dynamics simulation provides the time-evolution and chemical reactions in individual simulation voxels as well as the energy-landscape accounted for the DNA-(•)OH chemical reaction that is essential for the first-principle enumeration of hydrogen abstractions, chemical bond breaks, and DNA-lesions induced by collection of ions in clusters less than the critical dimension which is approximately 2-3 Å. We show that the formation of broken bonds leads to DNA-base and backbone damages that collectively propagate to DNA single and double-strand breaks. For illustration of the methodology, we focused on particles with an initial energy of 1 MeV. Our studies reveal a qualitative difference in DNA damage induced by low energy electrons and protons. Electrons mainly generate small pockets of (•)OH-radicals, randomly dispersed in the cell volume. In contrast, protons generate larger clusters along a straight-line parallel to the direction of the particle. The ratio of the total DNA double-strand breaks induced by a single proton and electron track is determined to be ≈4 in the linear scaling limit. In summary, we have developed a multi-scale computational model based on first-principles to study the interaction of ionizing radiation with DNA molecules. The main advantage of our hybrid Monte Carlo approach using Geant4-DNA and ReaxFF is the multi-scale simulation of the cascade of both physical and chemical events which result in the formation of biological damage. The tool developed in this work can be used in the future to investigate the relative biological effectiveness of light and heavy ions that are used in radiotherapy.

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

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

Monte Carlo damage models of different complexity levels predict similar trends in radiation induced DNA damage

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

Modeling of DNA Damage Repair and Cell Response in Relation to p53 System Exposed to Ionizing Radiation

A molecular dynamics simulation of DNA damage induction by ionizing radiation

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

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

Monte Carlo simulations of cell survival in proton SOBP

Renormalization of radiobiological response functions by energy loss fluctuations and complexities in chromosome aberration induction: deactivation theory for proton therapy from cells to tumor control

A mechanistic model of high dose irradiation damage

On the biophysical interpretation of lethal DNA lesions induced by ionizing radiation

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

Radiation-induced DNA damage by proton, helium and carbon ions in human fibroblast cell: Geant4-DNA and MCDS-based study

RBE with Non-Poisson Distribution of Radiation Induced Strand Breaks

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

Effect of overdispersion of lethal lesions on cell survival curves

Modelling of Cellular Survival Following Radiation-Induced DNA Double-Strand Breaks

Detailed analysis of the cell-inactivation mechanism by accelerated protons and light ions

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

The complexity of DNA damage by radiation follows a Gamma distribution: insights from the Microdosimetric Gamma Model