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.

Theoretical description of proton and helium ion elastic scattering in biological media

Configuration interaction calculations of cross sections for proton stopping reactions

Simulating the nanometric track-structure of carbon ion beams in liquid water at energies relevant for hadrontherapy

Electromagnetic interaction models for Monte Carlo simulation of protons and alpha particles

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

Inelastic scattering of electrons in water from first principles: cross sections and inelastic mean free path for use in Monte Carlo track-structure simulations of biological damage

Elastic, charge transfer, and related transport cross sections for proton impact of atomic hydrogen for astrophysical and laboratory plasma modeling

Non-perturbative scaling behavior for net ionization of biologically relevant molecules by multiply-charged heavy-ion impact

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

Calculation of lineal energies for water and DNA bases using the Rudd model cross sections integrated within the Geant4-DNA processes

A New Approach to Modeling the Microdosimetry of Proton Therapy Beams

Electron transport in DNA bases: An extension of the Geant4-DNA Monte Carlo toolkit

The neutron cross section of the hydrogen liquids: substantial improvements and perspectives

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

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

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

Experimental characterisation of a proton kernel model for pencil beam scanning techniques

Monte Carlo Calculation for Production Cross-Sections of Projectile’s Isotopes from Therapeutic Carbon and Helium Ion Beams in Different Materials

Ionization cross sections of water clusters as an alternative to liquid phase water: A Geant4-DNA simulation study

Ion-induced electron production in tissue-like media and DNA damage mechanisms

A molecular dynamics simulation of DNA damage induction by ionizing radiation