Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called “indirect” damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our under-standing of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

An upgraded version of the BIANCA II biophysical model, which describes more realistically interphase chromosome organization and the link between chromosome aberrations and cell death, was applied to V79 and AG01522 cells exposed to protons, C-ions and He-ions over a wide LET interval (0.6-502 keV µm), as well as proton-irradiated U87 cells. The model assumes that (i) ionizing radiation induces DNA 'cluster lesions' (CLs), where by definition each CL produces two independent chromosome fragments; (ii) fragment (distance-dependent) mis-rejoining, or un-rejoining, produces chromosome aberrations; (iii) some aberrations lead to cell death. The CL yield, which mainly depends on radiation quality but is also modulated by the target cell, is an adjustable parameter. The fragment un-rejoining probability, f, is the second, and last, parameter. The value of f, which is assumed to depend on the cell type but not on radiation quality, was taken from previous studies, and only the CL yield was adjusted in the present work. Good agreement between simulations and experimental data was obtained, suggesting that BIANCA II is suitable for calculating the biological effectiveness of hadrontherapy beams. For both V79 and AG01522 cells, the mean number of CLs per micrometer was found to increase with LET in a linear-quadratic fashion before the over-killing region, where a less rapid increase, with a tendency to saturation, was observed. Although the over-killing region deserves further investigation, the possibility of fitting the CL yields is an important feature for hadrontherapy, because it allows performing predictions also at LET values where experimental data are not available. Finally, an approach was proposed to predict the ion-response of the cell line(s) of interest from the ion-response of a reference cell line and the photon response of both. A pilot study on proton-irradiated AG01522 and U87 cells, taking V79 cells as a reference, showed encouraging results.

This paper presents a biophysical model of radiation-induced cell death, implemented as a Monte Carlo code called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), based on the assumption that some chromosome aberrations (dicentrics, rings, and large deletions, called ‘‘lethal aberrations’’) lead to clonogenic inactivation. In turn, chromosome aberrations are assumed to derive from clustered, and thus severe, DNA lesions (called ‘‘cluster lesions,’’ or CL) interacting at the micrometer scale; the CL yield and the threshold distance governing CL interaction are the only model parameters. After a pilot study on V79 hamster cells exposed to protons and carbon ions, in the present work the model was extended and applied to AG1522 human cells exposed to photons, He ions, and heavier ions including carbon and neon. The agreement with experimental survival data taken from the literature supported the assumptions. In particular, the inactivation of AG1522 cells was explained by lethal aberrations not only for X-rays, as already reported by others, but also for the aforementioned radiation types. Furthermore, the results are consistent with the hypothesis that the critical initial lesions leading to cell death are DNA cluster lesions having yields in the order of *2 CL Gy-1 cell-1 at low LET and*20 CL Gy-1 cell-1 at high LET, and that the processing of these lesions is modulated by proximity effects at the micrometer scale related to interphase chromatin organization. The model was then applied to calculate the fraction of inactivated cells, as well as the yields of lethal aberrations and cluster lesions, as a function of LET; the results showed a maximum around 130 keV/lm, and such maximum was much higher for cluster lesions and lethal aberrations than for cell inactivation.