Single- and Double-Strand Breaks of Dry DNA Exposed to Protons at Bragg-Peak Energies

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

Exposure to ionizing radiation can induce genetic aberrations via unrepaired DNA strand breaks. To investigate quantitatively the dose-effect relationship at the molecular level, we irradiated dry pBR322 plasmid DNA with 3 MeV protons and assessed fragmentation yields at different radiation doses using long-read sequencing from Oxford Nanopore Technologies. This technology applied to a reference DNA model revealed dose-dependent fragmentation, as evidenced by read length distributions, showing no discernible radiation sensitivity in specific genetic sequences. In addition, we propose a method for directly measuring the single-strand break (SSB) yield. Furthermore, through a comparative study with a collection of previous works on dry DNA irradiation, we show that the irradiation protocol leads to biases in the definition of ionizing sources. We support this scenario by discussing the size distributions of nanopore sequencing reads in the light of Geant4 and Geant4-DNA simulation toolkit predictions. We show that integrating long-read sequencing technologies with advanced Monte Carlo simulations paves a promising path toward advancing our comprehension and prediction of radiation-induced DNA fragmentation.

Effects of Differing Underlying Assumptions in In Silico Models on Predictions of DNA Damage and Repair

The induction and repair of DNA double-strand breaks (DSBs) are critical factors in the treatment of cancer by radiotherapy. To investigate the relationship between incident radiation and cell death through DSB induction many in silico models have been developed. These models produce and use custom formats of data, specific to the investigative aims of the researchers, and often focus on particular pairings of damage and repair models. In this work we use a standard format for reporting DNA damage to evaluate combinations of different, independently developed, models. We demonstrate the capacity of such inter-comparison to determine the sensitivity of models to both known and implicit assumptions. Specifically, we report on the impact of differences in assumptions regarding patterns of DNA damage induction on predicted initial DSB yield, and the subsequent effects this has on derived DNA repair models. The observed differences highlight the importance of considering initial DNA damage on the scale of nanometres rather than micrometres. We show that the differences in DNA damage models result in subsequent repair models assuming significantly different rates of random DSB end diffusion to compensate. This in turn leads to disagreement on the mechanisms responsible for different biological endpoints, particularly when different damage and repair models are combined, demonstrating the importance of inter-model comparisons to explore underlying model assumptions.

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

The electronic excitations caused by DNA when exposed to ion radiation is essential to DNA damage. In this paper, we investigated the energy deposition and electron excitation process of DNA with reasonable stretching range upon proton irradiation based on time-dependent density functional theory. Stretching changes the strength of hydrogen bonding between the DNA base pairs, which in turn affects the Coulomb interaction between the projectile and DNA. As a semi-flexible molecule, the way of energy deposition is weakly sensitive to the stretching rate of DNA. However, the increase of stretching rate causes the increase of charge density along the trajectory channel, sequentially resulting in an increase in proton resistance along the intruding channel. The Mulliken charge analysis indicates that the guanine base and guanine ribose are ionized, meanwhile the cytosine base and cytosine ribose are reduced at all stretching rates. In a few femtoseconds, there exists an electron flow passing through the guanine ribose, guanine, cytosine base and the cytosine ribose in turn. This electron flow increases electron transfer and DNA ionization, promoting the side chain damage of the DNA upon ion irradiation. Our results provide a theoretical insight for deciphering the physical mechanism of the early stage of the irradiation process, and are also of great significance for the study of particle beam cancer therapy in different biological tissues.

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

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keV μm−1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Applications of nanodosimetry in particle therapy planning and beyond

This topical review summarizes underlying concepts of nanodosimetry. It describes the development and current status of nanodosimetric detector technology. It also gives an overview of Monte Carlo track structure simulations that can provide nanodosimetric parameters for treatment planning of proton and ion therapy. Classical and modern radiobiological assays that can be used to demonstrate the relationship between the frequency and complexity of DNA lesion clusters and nanodosimetric parameters are reviewed. At the end of the review, existing approaches of treatment planning based on relative biological effectiveness (RBE) models or dose-averaged linear energy transfer are contrasted with an RBE-independent approach based on nandosimetric parameters. Beyond treatment planning, nanodosimetry is also expected to have applications and give new insights into radiation protection dosimetry.

Ionization clustering on charged particle tracks as a seed for biologically relevant radiation effects

We construct a statistical framework for investigating the physical origins of radiation effects on biological materials and report the fit of an analytical statistical model to rates of simple lesions in DNA. Modeling primary ionization damage as trails of electron vacancies left on the trajectories of fast charged particles, we derive the dependence of rates of spatial clustering on ionization density [linear energy transfer (LET)]; a clustering scale parameter, r_{0}; and number per cluster. Published experimental results on rates of single strand breaks and base lesions in dry DNA over a range of LET are fitted with the derived functions, assuming clusters of 1 or ≥1. The fits yield reasonable goodness of fit and values of r_{0} that are consistent with expectations. Limitations of the model and future developments are discussed. This framework may ultimately contribute to an improved understanding of the physical origins of biological radiation effects.

First-Principles Modeling of Electronic Stopping in Complex Matter under Ion Irradiation

Electronic stopping refers to the dynamical energy-transfer process to electrons in matter from highly energetic charged particles such as high-velocity protons. We discuss recent progress in theoretical studies of electronic stopping in condensed matter under ion irradiation, focusing on modern electronic structure theory's role in enabling the study of electronic excitation dynamics that result from the energy transfer. In the last few decades, first-principles simulation approaches based on real-time time-dependent density functional theory have greatly advanced the field. While linear response theory is widely used to study electronic stopping processes, especially for simple solids, novel first-principles dynamics approaches now allow us to study chemically complex systems and also yield detailed descriptions of electronic excitations at the molecular scale. Outstanding challenges for further advancement of electronic stopping modeling are also discussed from the viewpoint of electronic structure theory.

Role of the pre-plasma on electron beam currents from a biased laser-plasma

We are investigating laser-plasmas produced in the interaction of a 1 J 9 ns Nd:YAG laser with a solid metal target as a source of electrons. An electron beam pulsed at the laser repetition rate is produced by biasing the target and making the plasma expand in an electric field. In this paper, we focus on the measured beam currents and charge surface distribution of the beam. The peak beam currents are much higher than what is given by a simplified toy model based on the Child-Langmuir limit in a vacuum and the charge surface distributions are inhomogeneous. Both these observations are explained by the presence of a positive preplasma expanding ahead of the laser-plasma front edge.

Clinically relevant nanodosimetric simulation of DNA damage complexity from photons and protons

Relative Biological Effectiveness (RBE), the ratio of doses between radiation modalities to produce the same biological endpoint, is a controversial and important topic in proton therapy. A number of phenomenological models incorporate variable RBE as a function of Linear Energy Transfer (LET), though a lack of mechanistic description limits their applicability. In this work we take a different approach, using a track structure model employing fundamental physics and chemistry to make predictions of proton and photon induced DNA damage, the first step in the mechanism of radiation-induced cell death. We apply this model to a proton therapy clinical case showing, for the first time, predictions of DNA damage on a patient treatment plan. Our model predictions are for an idealised cell and are applied to an ependymoma case, at this stage without any cell specific parameters. By comparing to similar predictions for photons, we present a voxel-wise RBE of DNA damage complexity. This RBE of damage complexity shows similar trends to the expected RBE for cell kill, implying that damage complexity is an important factor in DNA repair and therefore biological effect.

Relative Biological Effectiveness (RBE) is a controversial and important topic in proton therapy. This work uses Monte Carlo simulations of DNA damage for protons and photons to probe this phenomenon, providing a plausible mechanistic understanding.

Mechanistic DNA damage simulations in Geant4-DNA Part 2: Electron and proton damage in a bacterial cell

We extended a generic Geant4 application for mechanistic DNA damage simulations to an Escherichia coli cell geometry, finding electron damage yields and proton damage yields largely in line with experimental results. Depending on the simulation of radical scavenging, electrons double strand breaks (DSBs) yields range from 0.004 to 0.010 DSB Gy-1 Mbp-1, while protons have yields ranging from 0.004 DSB Gy-1 Mbp-1 at low LETs and with strict assumptions concerning scavenging, up to 0.020 DSB Gy-1 Mbp-1 at high LETs and when scavenging is weakest. Mechanistic DNA damage simulations can provide important limits on the extent to which physical processes can impact biology in low background experiments. We demonstrate the utility of these studies for low dose radiation biology calculating that in E. coli, the median rate at which the radiation background induces double strand breaks is 2.8 × 10-8 DSB day-1, significantly less than the mutation rate per generation measured in E. coli, which is on the order of 10-3.