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

Background. Radiation-induced DNA damages such as Single Strand Break (SSB), Double Strand Break (DSB) and Complex DSB (cDSB) are critical aspects of radiobiology with implications in radiotherapy and radiation protection applications.Materials and Methods. This study presents a thorough investigation into the effects of protons (0.1-100 MeV/u), helium ions (0.13-100 MeV/u) and carbon ions (0.5-480 MeV/u) on DNA of human fibroblast cells using Geant4-DNA track structure code coupled with DBSCAN algorithm and Monte Carlo Damage Simulations (MCDS) code. Geant4-DNA-based simulations consider 1μm × 1μm × 0.5μm water box as the target to calculate energy deposition on event-by-event basis and the three-dimensional coordinates of the interaction location, and then DBSCAN algorithm is used to calculate yields of SSB, DSB and cDSB in human fibroblast cell. The study investigated the influence of Linear Energy Transfer (LET) of protons, helium ions and carbon ions on the yields of DNA damages. Influence of cellular oxygenation on DNA damage patterns is investigated using MCDS code.Results. The study shows that DSB and SSB yields are influenced by the LET of the particles, with distinct trends observed for different particles. The cellular oxygenation is a key factor, with anoxic cells exhibiting reduced SSB and DSB yields, underscoring the intricate relationship between cellular oxygen levels and DNA damage. The study introduced DSB/SSB ratio as an informative metric for evaluating the severity of radiation-induced DNA damage, particularly in higher LET regions.Conclusions. The study highlights the importance of considering particle type, LET, and cellular oxygenation in assessing the biological effects of ionizing radiation.

GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species

Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.

A matter of space: how the spatial heterogeneity in energy deposition determines the biological outcome of radiation exposure

The outcome of the exposure of living organisms to ionizing radiation is determined by the distribution of the associated energy deposition at different spatial scales. Radiation proceeds through ionizations and excitations of hit molecules with an ~ nm spacing. Approaches such as nanodosimetry/microdosimetry and Monte Carlo track-structure simulations have been successfully adopted to investigate radiation quality effects: they allow to explore correlations between the spatial clustering of such energy depositions at the scales of DNA or chromosome domains and their biological consequences at the cellular level. Physical features alone, however, are not enough to assess the entity and complexity of radiation-induced DNA damage: this latter is the result of an interplay between radiation track structure and the spatial architecture of chromatin, and further depends on the chromatin dynamic response, affecting the activation and efficiency of the repair machinery. The heterogeneity of radiation energy depositions at the single-cell level affects the trade-off between cell inactivation and induction of viable mutations and hence influences radiation-induced carcinogenesis. In radiation therapy, where the goal is cancer cell inactivation, the delivery of a homogenous dose to the tumour has been the traditional approach in clinical practice. However, evidence is accumulating that introducing heterogeneity with spatially fractionated beams (mini- and microbeam therapy) can lead to significant advantages, particularly in sparing normal tissues. Such findings cannot be explained in merely physical terms, and their interpretation requires considering the scales at play in the underlying biological mechanisms, suggesting a systemic response to radiation.

Biosensor for deconvolution of individual cell fate in response to ion beam irradiation

Clonogenic survival assay constitutes the gold standard method for quantifying radiobiological effects. However, it neglects cellular radiation response variability and heterogeneous energy deposition by ion beams on the microscopic scale. We introduce "Cell-Fit-HD4D" a biosensor that enables a deconvolution of individual cell fate in response to the microscopic energy deposition as visualized by optical microscopy. Cell-Fit-HD4D enables single-cell dosimetry in clinically relevant complex radiation fields by correlating microscopic beam parameters with biological endpoints. Decrypting the ion beam's energy deposition and molecular effects at the single-cell level has the potential to improve our understanding of radiobiological dose concepts as well as radiobiological study approaches in general.

Towards the characterization of neutron carcinogenesis through direct action simulations of clustered DNA damage

Neutron exposure poses a unique radiation protection concern because neutrons have a large, energy-dependent relative biological effectiveness (RBE) for stochastic effects. Recent computational studies on the microdosimetric properties of neutron dose deposition have implicated clustered DNA damage as a likely contributor to this marked energy dependence. So far, publications have focused solely on neutron RBE for inducing clusters of DNA damage containing two or more DNA double strand breaks (DSBs). In this study, we have conducted a novel assessment of neutron RBE for inducing all types of clustered DNA damage that contain two or more lesions, stratified by whether the clusters contain DSBs (complex DSB clusters) or not (non-DSB clusters). This assessment was conducted for eighteen initial neutron energies between 1 eV and 10 MeV as well as a reference radiation of 250 keV x-rays. We also examined the energy dependence of cluster length and cluster complexity because these factors are believed to impact the DNA repair process. To carry out our investigation, we developed a user-friendly TOPAS-nBio application that includes a custom nuclear DNA model and a novel algorithm for recording clustered DNA damage. We found that neutron RBE for inducing complex DSB clusters exhibited similar energy dependence to the canonical neutron RBE for stochastic radiobiological effects, at multiple depths in human tissue. Qualitatively similar results were obtained for non-DSB clusters, although the quantitative agreement was lower. Additionally we identified a significant neutron energy dependence in the average length and complexity of clustered lesions. These results support the idea that many types of clustered DNA damage contribute to the energy dependence of neutron RBE for stochastic radiobiological effects and imply that the size and constituent lesions of individual clusters should be taken into account when modeling DNA repair. Our results were qualitatively consistent for (i) multiple radiation doses (including a low-dose 0.1 Gy irradiation), (ii) variations in the maximal lesion separation distance used to define a cluster, and (iii) two distinct collections of physics models used to govern particle transport. Our complete TOPAS-nBio application has been released under an open-source license to enable others to independently validate our work and to expand upon it.

Prospects for Radiopharmaceuticals as Effective and Safe Therapeutics in Oncology and Challenges of Tumor Resistance to Radiotherapy

The rapid advances in nuclear medicine have resulted in significant advantages for the field of oncology. The focus is on the application of radiopharmaceuticals as therapeuticals. In addition, the latest developments in cell biology (the understanding of the cell structure, function, metabolism, genetics, signaling, transformation) have given a strong scientific boost to radiation oncology. In this regard, the article discusses what is soon going to be a new jump in radiation oncology based on the already accumulated considerable knowledge at the cellular level about the mechanisms of cell transformation and tumor progression, cell response to radiation, cell resistance to apoptosis and radiation and cell radio-sensitivity. The mechanisms of resistance of tumor cells to radiation and the genetically determined individual sensitivity to radiation in patients (which creates the risk of radiation-induced acute and late side effects) are the 2 major challenges to overcome in modern nuclear medicine. The paper focuses on these problems and makes a detailed summary of the significance of the differences in the ionizing properties of radiopharmaceuticals and the principle of their application in radiation oncology that will shed additional light on how to make the anti-cancer radiotherapies more efficient and safe, giving some ideas for optimizations.

PROPOSAL FOR A EUROPEAN METROLOGY NETWORK ON BIOLOGICAL IONISING RADIATION EFFECTS

Progress in the field of ionising radiation (IR) metrology achieved in the BioQuaRT project raised the question to what extent radiobiological investigations would benefit from metrological support of the applied methodologies. A panel of experts from the medical field, fundamental research and radiation protection attended a workshop at Physikalisch-Technische Bundesanstalt to consult on metrology needs related to biological radiation effects. The panel identified a number of metrological needs including the further development of experimental and computational techniques for micro- and nanodosimetry, together with the determination of related fundamental material properties and the establishment of rigorous uncertainty budgets. In addition to this, a call to develop a metrology support for assisting quality assurance of radiobiology experiments was expressed. Conclusions from the workshop were presented at several international conferences for further discussion with the scientific community and stakeholder groups that led to an initiative within the metrology community to establish a European Metrology Network on biological effects of IR.

Assessment of Radio-Induced Damage in Endothelial Cells Irradiated with 40 kVp, 220 kVp, and 4 MV X-rays by Means of Micro and Nanodosimetric Calculations

The objective of this work was to study the differences in terms of early biological effects that might exist between different X-rays energies by using a mechanistic approach. To this end, radiobiological experiments exposing cell monolayers to three X-ray energies were performed in order to assess the yields of early DNA damage, in particular of double-strand breaks (DSBs). The simulation of these irradiations was set in order to understand the differences in the obtained experimental results. Hence, simulated results in terms of microdosimetric spectra and early DSB induction were analyzed and compared to the experimental data. Human umbilical vein endothelial cells (HUVECs) were irradiated with 40, 220 kVp, and 4 MV X-rays. The Geant4 Monte Carlo simulation toolkit and its extension Geant4-DNA were used for the simulations. Microdosimetric calculations aiming to determine possible differences in the variability of the energy absorbed by the irradiated cell population for those photon spectra were performed on 10,000 endothelial cell nuclei representing a cell monolayer. Nanodosimetric simulations were also carried out using a computation chain that allowed the simulation of physical, physico-chemical, and chemical stages on a single realistic endothelial cell nucleus model including both heterochromatin and euchromatin. DNA damage was scored in terms of yields of prompt DSBs per Gray (Gy) and per giga (10⁹) base pair (Gbp) and DSB complexity was derived in order to be compared to experimental data expressed as numbers of histone variant H2AX (γ-H2AX) foci per cell. The calculated microdosimetric spread in the irradiated cell population was similar when comparing between 40 and 220 kVp X-rays and higher when comparing with 4 MV X-rays. Simulated yields of induced DSB/Gy/Gbp were found to be equivalent to those for 40 and 220 kVp but larger than those for 4 MV, resulting in a relative biological effectiveness (RBE) of 1.3. Additionally, DSB complexity was similar between the considered photon spectra. Simulated results were in good agreement with experimental data obtained by IRSN (Institut de radioprotection et de sûreté nucléaire) radiobiologists. Despite differences in photon energy, few differences were observed when comparing between 40 and 220 kVp X-rays in microdosimetric and nanodosimetric calculations. Nevertheless, variations were observed when comparing between 40/220 kVp and 4 MV X-rays. Thanks to the simulation results, these variations were able to be explained by the differences in the production of secondary electrons with energies below 10 keV.

Predicting DNA damage foci and their experimental readout with 2D microscopy: a unified approach applied to photon and neutron exposures

The consideration of how a given technique affects results of experimental measurements is a must to achieve correct data interpretation. This might be challenging when it comes to measurements on biological systems, where it is unrealistic to have full control (e.g. through a software replica) of all steps in the measurement chain. In this work we address how the effectiveness of different radiation qualities in inducing biological damage can be assessed measuring DNA damage foci yields, only provided that artefacts related to the scoring technique are adequately considered. To this aim, we developed a unified stochastic modelling approach that, starting from radiation tracks, predicts both the induction, spatial distribution and complexity of DNA damage, and the experimental readout of foci when immunocytochemistry coupled to 2D fluorescence microscopy is used. The approach is used to interpret γ-H2AX data for photon and neutron exposures. When foci are reconstructed in the whole cell nucleus, we obtain information on damage characteristics "behind" experimental observations, as the average damage content of a focus. We reproduce how the detection technique affects experimental findings, e.g. contributing to the saturation of foci yields scored at 30 minutes after exposure with increasing dose and to the lack of dose dependence for yields at 24 hours.

PROSPECTS FOR METROLOGY RELATED TO BIOLOGICAL RADIATION EFFECTS OF ION BEAMS

In recent years, several approaches have been proposed to provide an understanding of the enhanced relative biological effectiveness of ion beams based on multi-scale models of their radiation effects. Among these, the BioQuaRT project was the only one which focused on developing metrology for a multi-scale characterization of particle track structure. The progress made within the BioQuaRT project has motivated the formation of a department 'Radiation Effects' at PTB dedicated to metrological research on ionizing radiation effects. This paper gives an overview of the department's present research directions and shortly discusses ideas for the future development of metrology related to biological effects of ion beams that are based on a stakeholder consultation.

WHAT ROLES FOR TRACK-STRUCTURE AND MICRODOSIMETRY IN THE ERA OF -omics AND SYSTEMS BIOLOGY?

Ionizing radiation is a peculiar perturbation when it comes to damage to biological systems: it proceeds through discrete energy depositions, over a short temporal scale and a spatial scale critical for subcellular targets as DNA, whose damage complexity determines the outcome of the exposure. This lies at the basis of the success of track structure (and nanodosimetry) and microdosimetry in radiation biology. However, such reductionist approaches cannot account for the complex network of interactions regulating the overall response of the system to radiation, particularly when effects are manifest at the supracellular level and involve long times. Systems radiation biology is increasingly gaining ground, but the gap between reductionist and holistic approaches is becoming larger. This paper presents considerations on what roles track structure and microdosimetry can have in the attempt to fill this gap, and on how they can be further exploited to interpret radiobiological data and inform systemic approaches.

From Energy Deposition of Ionizing Radiation to Cell Damage Signaling: Benchmarking Simulations by Measured Yields of Initial DNA Damage after Ion Microbeam Irradiation

Advances in accelerator technology, which have enabled conforming radiotherapy with charged hadronic species, have brought benefits as well as potential new risks to patients. To better understand the effects of ionizing radiation on tumor and surrounding tissue, it is important to investigate and quantify the relationship between energy deposition at the nanometric scale and the initial biological events. Monte Carlo track structure simulation codes provide a powerful tool for investigating this relationship; however, their success and reliability are dependent on their improvement and development accordingly to the dedicated biological data to which they are challenged. For this aim, a microbeam facility that allows for fluence control, down to one ion per cell nucleus, was used to evaluate relative frequencies of DNA damage after interaction between the incoming ion and DNA according to radiation quality. Primary human cells were exposed to alpha particles of three different energies with respective linear energy transfers (LETs) of approximately 36, 85 or 170 keV·µm^-1 at the cells' center position, or to protons (19 keV·µm^-1). Statistical evaluation of nuclear foci formation (53BP1/γ-H2AX), observed using immunofluorescence and related to a particle traversal, was undertaken in a large population of cell nuclei. The biological results were adjusted to consider the factors that drive the experimental uncertainties, then challenged with results using Geant4-DNA code modeling of the ionizing particle interactions on a virtual phantom of the cell nucleus with the same mean geometry and DNA density as the cells used in our experiments. Both results showed an increase of relative frequencies of foci (or simulated DNA damage) in cell nuclei as a function of increasing LET of the traversing particles, reaching a quasi-plateau when the LET exceeded 80-90 keV·µm^-1. For the LET of an alpha particle ranging from 80-90 to 170 keV·µm^-1, 10-30% of the particle hits did not lead to DNA damage inducing 53BP1 or γ-H2AX foci formation.

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.