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.

Changes in molecular conformation and electronic structure of DNA under 12 C ions based on first-principles calculations

Investigation of DNA Damage Induced by ^7li and ^12C Ions

A molecular dynamics simulation of DNA damage induction by ionizing radiation

Ion Type Dependence of DNA Electronic Excitation in Water under Proton, alpha-particle, and Carbon Ion Irradiation: A First-Principles Simulation Study

Ion irradiation induced direct damage to DNA

Progress in the Mechanism Study on DNA Strand Breaks in Vitro Induced by Ionizing Radiation

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

Formation of clustered DNA damage in vivo upon irradiation with ionizing radiation: Visualization and analysis with atomic force microscopy

Ultrafast Nonadiabatic Dynamics of Λ-Dna Upon Low Energy Proton Irradiation

Radiation and DNA Origami Nanotechnology: Probing Structural Integrity at the Nanoscale

Quantum mechanical calculations related to ionization and charge transfer in DNA

Stretching effects on non-adiabatic electron dynamic behavior in poly(dG)-poly(dC) DNA upon the proton irradiation

Research on Monte Carlo Simulation of DNA Strand Break Induced by Low Energy Photon

Radiation‐Induced Molecular Processes in DNA: A Perspective on Gas‐Phase Interaction Studies

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

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

DFT calculation on the energy thresholds of DNA damages under irradiation conditions

Kinetics of Ionizing Irradiation-induced DNA Strand Breaks.

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

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

Quantitative analysis of dose dependent DNA fragmentation in dry pBR322 plasmid using long read sequencing and Monte Carlo simulations