Chromosome aberrations and cell death by ionizing radiation: Evolution of a biophysical model

The manuscript summarizes and discusses the various versions of a radiation damage biophysical model, implemented as a Monte Carlo simulation code, originally developed for chromosome aberrations and subsequently extended to cell death. This extended version has been called BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations). According to the basic assumptions, complex double-strand breaks (called “Cluster Lesions”, or CLs) produce independent chromosome free-ends, mis-rejoining within a threshold distance d (or un-rejoining) leads to chromosome aberrations, and “lethal aberrations” (i.e., dicentrics plus rings plus large deletions) lead to clonogenic cell death. The mean number of CLs per Gy and per cell is an adjustable parameter. While in BIANCA the threshold distance d was the second parameter, in a subsequent version, called BIANCA II, d has been fixed as the mean distance between two adjacent interphase chromosome territories, and a new parameter, f, has been introduced to represent the chromosome free-end un-rejoining probability. Simulated dose-response curves for chromosome aberrations and cell survival obtained by the various model versions were compared with literature experimental data. Such comparisons provided indications on some open questions, including the role of energy deposition clustering at the nm and the µm level, the probability for a chromosome free-end to remain un-rejoined, and the relationship between chromosome aberrations and cell death. Although both BIANCA and BIANCA II provided cell survival curves in general agreement with human and hamster fibroblast survival data, BIANCA II allowed for a better reproduction of dicentrics, rings and deletions considered separately. Furthermore, the approach adopted in BIANCA II for d is more consistent with estimates reported in the literature. After testing against aberration and survival data, BIANCA II was applied to investigate the depth-dependence of the radiation effectiveness for a proton SOBP used to treat eye melanoma in Catania, Italy. The survival of AG01522 cells at different depths was reproduced, and the survival of V79 cells was predicted. For both cell lines, the simulations also predicted yields of chromosome aberrations, some of which can be regarded as indicators of the risk to normal tissues.