A Model of Chromosome Aberration Induction: Applications to Space Research

A method to predict space radiation biological effectiveness for non-cancer effects following intense Solar Particle Events

In addition to the continuous exposure to cosmic rays, astronauts in space are occasionally exposed to Solar Particle Events (SPE), which involve less energetic particles but can deliver much higher doses. The latter can exceed several Gy in a few hours for the most intense SPEs, for which non-stochastic effects are thus a major concern. To identify adequate shielding conditions that would allow respecting the dose limits established by the various space agencies, the absorbed dose in the considered organ/tissue must be multiplied by the corresponding Relative Biological Effectiveness (RBE), which is a complex quantity depending on several factors including particle type and energy, considered biological effect, level of effect (and thus absorbed dose), etc. While in several studies only the particle-type dependence of RBE is taken into account, in this work we developed and applied a new approach where, thanks to an interface between the FLUKA Monte Carlo transport code and the BIANCA biophysical model, the RBE dependence on particle energy and absorbed dose was also considered. Furthermore, we included in the considered SPE spectra primary particles heavier than protons, which in many studies are neglected. This approach was then applied to the October 2003 SPE (the most intense SPE of solar cycle 23, also known as "Halloween event") and the January 2005 event, which was characterized by a lower fluence but a harder spectrum, i.e., with higher-energy particles. The calculation outcomes were then discussed and compared with the current dose limits established for skin and blood forming organs in case of 30-days missions. This work showed that the BIANCA model, if interfaced to a radiation transport code, can be used to calculate the RBE values associated to Solar Particle Events. More generally, this work emphasizes the importance of taking into account the RBE dependence on particle energy and dose when calculating equivalent doses.

Radiation Damage in Biomolecules and Cells 2.0

It is well known that ionizing radiation, when it hits living cells, causes a plethora of different damage types at different levels [...].

A Mission to Mars: Prediction of GCR Doses and Comparison with Astronaut Dose Limits

Long-term human space missions such as a future journey to Mars could be characterized by several hazards, among which radiation is one the highest-priority problems for astronaut health. In this work, exploiting a pre-existing interface between the BIANCA biophysical model and the FLUKA Monte Carlo transport code, a study was performed to calculate astronaut absorbed doses and equivalent doses following GCR exposure under different shielding conditions. More specifically, the interface with BIANCA allowed us to calculate both the RBE for cell survival, which is related to non-cancer effects, and that for chromosome aberrations, related to the induction of stochastic effects, including cancer. The results were then compared with cancer and non-cancer astronaut dose limits. Concerning the stochastic effects, the equivalent doses calculated by multiplying the absorbed dose by the RBE for chromosome aberrations ("high-dose method") were similar to those calculated using the Q-values recommended by ICRP. For a 650-day mission at solar minimum (representative of a possible Mars mission scenario), the obtained values are always lower than the career limit recommended by ICRP (1 Sv), but higher than the limit of 600 mSv recently adopted by NASA. The comparison with the JAXA limits is more complex, since they are age and sex dependent. Concerning the deterministic limits, even for a 650-day mission at solar minimum, the values obtained by multiplying the absorbed dose by the RBE for cell survival are largely below the limits established by the various space agencies. Following this work, BIANCA, interfaced with an MC transport code such as FLUKA, can now predict RBE values for cell death and chromosome aberrations following GCR exposure. More generally, both at solar minimum and at solar maximum, shielding of 10 g/cm² Al seems to be a better choice than 20 g/cm² for astronaut protection against GCR.

Radiobiological damage by space radiation: extension of the BIANCA model to heavy ions up to iron, and pilot application to cosmic ray exposure

Space research seems to be object of a renewed interest, also considering that human missions to the Moon, and possibly Mars, are being planned. Among the risks affecting such missions, astronauts' exposure to space radiation is a major concern. In this work, the question of the evaluation of biological damage by Galactic Cosmic Rays (GCR) was addressed by a biophysical model called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), which simulates the induction of cell death and chromosome aberrations by different ions. While previously BIANCA has been validated for calculating cell death along hadrontherapy beams up to oxygen, herein the approach was extended up to Fe ions. Specifically, experimental survival curves available in literature for V79 cells irradiated by Si-, Ne-, Ar- and Fe-ions were reproduced, and a reference radiobiological database describing V79 cell survival as a function of ion type (1 ⩽Z⩽ 26), energy and dose was constructed. Analogous databases were generated for Chinese hamster ovary hamster cells and human skin fibroblasts, finding good agreement between simulations and data. Concerning chromosome aberrations, which are regarded as radiation risk biomarkers, dicentric data in human lymphocytes irradiated by heavy ions up to iron were reproduced, and a radiobiological database allowing calculation of lymphocyte dicentric yields as a function of dose, ion type (1 ⩽Z⩽ 26) and energy was constructed. Following interface between BIANCA and the FLUKA Monte Carlo transport code, a feasibility study was performed to calculate the relative biological effectiveness (RBE) of different GCR spectrum components, for both dicentrics and cell death. Fe-ions, although representing only 10% of the total absorbed dose, were found to be responsible for about 35%-40% of the RBE-weighted dose. Interestingly, the RBE for dicentrics was higher than that for cell survival. More generally, this work shows that BIANCA can calculate RBE values for cell death and lymphocyte dicentrics not only for ion therapy, but also for space radiation.

Healthy Tissue Damage Following Cancer Ion Therapy: A Radiobiological Database Predicting Lymphocyte Chromosome Aberrations Based on the BIANCA Biophysical Model

Chromosome aberrations are widely considered among the best biomarkers of radiation health risk due to their relationship with late cancer incidence. In particular, aberrations in peripheral blood lymphocytes (PBL) can be regarded as indicators of hematologic toxicity, which is a major limiting factor of radiotherapy total dose. In this framework, a radiobiological database describing the induction of PBL dicentrics as a function of ion type and energy was developed by means of the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model, which has been previously applied to predict the effectiveness of therapeutic-like ion beams at killing tumour cells. This database was then read by the FLUKA Monte Carlo transport code, thus allowing us to calculate the Relative Biological Effectiveness (RBE) for dicentric induction along therapeutic C-ion beams. A comparison with previous results showed that, while in the higher-dose regions (e.g., the Spread-Out Bragg Peak, SOBP), the RBE for dicentrics was lower than that for cell survival. In the lower-dose regions (e.g., the fragmentation tail), the opposite trend was observed. This work suggests that, at least for some irradiation scenarios, calculating the biological effectiveness of a hadrontherapy beam solely based on the RBE for cell survival may lead to an underestimation of the risk of (late) damage to healthy tissues. More generally, following this work, BIANCA has gained the capability of providing RBE predictions not only for cell killing, but also for healthy tissue damage.

First benchmarking of the BIANCA model for cell survival prediction in a clinical hadron therapy scenario

In the framework of RBE modelling for hadron therapy, the BIANCA biophysical model was extended to O-ions and was used to construct a radiobiological database describing the survival of V79 cells as a function of ion type (1  ⩽  Z  ⩽  8) and energy. This database allowed performing RBE predictions in very good agreement with experimental data. A method was then developed to construct analogous databases for different cell lines, starting from the V79 database as a reference. Following interface to the FLUKA Monte Carlo radiation transport code, BIANCA was then applied for the first time to predict cell survival in a typical patient treatment scenario, consisting of two opposing fields of range-equivalent protons or C-ions. The model predictions were found to be in good agreement with CHO cell survival data obtained at the Heidelberg ion-beam therapy (HIT) centre, as well as predictions performed by the local effect model (version LEM IV). This work shows that BIANCA can be used to predict cell survival and RBE not only for V79 and AG01522 cells, as shown previously, but also, in principle, for any cell line of interest. Furthermore, following interface to a transport code like FLUKA, BIANCA can provide predictions of 3D biological dose distributions for hadron therapy treatments, thus laying the foundations for future applications in clinics.

Analysis of Radiation-Induced Chromosomal Aberrations on a Cell-by-Cell Basis after Alpha-Particle Microbeam Irradiation: Experimental Data and Simulations

There is a continued need for further clarification of various aspects of radiation-induced chromosomal aberration, including its correlation with radiation track structure. As part of the EMRP joint research project, Biologically Weighted Quantities in Radiotherapy (BioQuaRT), we performed experimental and theoretical analyses on chromosomal aberrations in Chinese hamster ovary cells (CHO-K1) exposed to α particles with final energies of 5.5 and 17.8 MeV (absorbed doses: ∼2.3 Gy and ∼1.9 Gy, respectively), which were generated by the microbeam at the Physikalisch-Technische Bundesanstalt (PTB) in Braunschweig, Germany. In line with the differences in linear energy transfer (approximately 85 keV/μm for 5.5 MeV and 36 keV/μm for 17.8 MeV α particles), the 5.5 MeV α particles were more effective than the 17.8 MeV α particles, both in terms of the percentage of aberrant cells (57% vs. 33%) and aberration frequency. The yield of total aberrations increased by a factor of ∼2, although the increase in dicentrics plus centric rings was less pronounced than in acentric fragments. The experimental data were compared with Monte Carlo simulations based on the BIophysical ANalysis of Cell death and chromosomal Aberrations model (BIANCA). This comparison allowed interpretation of the results in terms of critical DNA damage [cluster lesions (CLs)]. More specifically, the higher aberration yields observed for the 5.5 MeV α particles were explained by taking into account that, although the nucleus was traversed by fewer particles (nominally, 11 vs. 25), each particle was much more effective (by a factor of ∼3) at inducing CLs. This led to an increased yield of CLs per cell (by a factor of ∼1.4), consistent with the increased yield of total aberrations observed in the experiments.

BIANCA, a biophysical model of cell survival and chromosome damage by protons, C-ions and He-ions at energies and doses used in hadrontherapy

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^-1), 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.

Proximity effects in chromosome aberration induction: Dependence on radiation quality, cell type and dose

It is widely accepted that, in chromosome-aberration induction, the (mis-)rejoining probability of two chromosome fragments depends on their initial distance, r. However, several aspects of these "proximity effects" need to be clarified, also considering that they can vary with radiation quality, cell type and dose. A previous work performed by the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model has suggested that, in human lymphocytes and fibroblasts exposed to low-LET radiation, an exponential function of the form exp(-r/r₀), which is consistent with free-end (confined) diffusion, describes proximity effects better than a Gaussian function. Herein, the investigation was extended to intermediate- and high-LET. Since the r₀ values (0.8 μm for lymphocytes and 0.7 μm for fibroblasts) were taken from the low-LET study, the results were obtained by adjusting only one model parameter, i.e. the yield of "Cluster Lesions" (CLs), where a CL was defined as a critical DNA damage producing two independent chromosome fragments. In lymphocytes, the exponential model allowed reproducing both dose-response curves for different aberrations (dicentrics, centric rings and excess acentrics), and values of F-ratio (dicentrics to centric rings) and G-ratio (interstitial deletions to centric rings). In fibroblasts, a good correspondence was found with the dose-response curves, whereas the G-ratio (and, to a lesser extent, the F-ratio) was underestimated. With increasing LET, F decreased and G increased in both cell types, supporting their role as "fingerprints" of high-LET exposure. A dose-dependence was also found at high LET, where F increased with dose and G decreased, possibly due to inter-track effects. We therefore conclude that, independent of radiation quality, in lymphocytes an exponential function can describe proximity effects at both inter- and intra-chromosomal level; on the contrary, in fibroblasts further studies (experimental and theoretical) are needed to explain the strong bias for intra-arm relative to inter-arm exchanges.

Proximity effects in chromosome aberration induction by low-LET ionizing radiation

Although chromosome aberrations are known to derive from distance-dependent mis-rejoining of chromosome fragments, evaluating whether a certain model describes such "proximity effects" better than another one is complicated by the fact that different approaches have often been tested under different conditions. Herein, a biophysical model ("BIANCA", i.e. BIophysical ANalysis of Cell death and chromosome Aberrations) was upgraded, implementing explicit chromosome-arm domains and two new models for the dependence of the rejoining probability on the fragment initial distance, r. Such probability was described either by an exponential function like exp(-r/r₀), or by a Gaussian function like exp(-r²/2σ²), where r₀ and σ were adjustable parameters. The second, and last, parameters was the yield of "Cluster Lesions" (CL), where "Cluster Lesion" defines a critical DNA damage producing two independent chromosome fragments. The model was applied to low-LET-irradiated lymphocytes (doses: 1-4Gy) and fibroblasts (1-6.1Gy). Good agreement with experimental yields of dicentrics and centric rings, and thus their ratio ("F-ratio"), was found by both the exponential model (with r₀=0.8μm for lymphocytes and 0.7μm for fibroblasts) and the Gaussian model (with σ=1.1μm for lymphocytes and 1.3μm for fibroblasts). While the former also allowed reproducing dose-responses for excess acentric fragments, the latter substantially underestimated the experimental curves. Both models provided G-ratios (ratio of acentric to centric rings) higher than those expected from randomness, although the values calculated by the Gaussian model were lower than those calculated by the exponential one. For lymphocytes the calculated G-ratios were in good agreement with the experimental ones, whereas for fibroblasts both models substantially underestimated the experimental results, which deserves further investigation. This work suggested that, although both models performed better than a step model (which previously allowed reproducing the F-ratio but underestimated the G-ratio), an exponential function describes proximity effects better than a Gaussian one.

Radiation track, DNA damage and response-a review

The purpose of this paper has been to review the current status and progress of the field of radiation biophysics, and draw attention to the fact that physics, in general, and radiation physics in particular, with the aid of mathematical modeling, can help elucidate biological mechanisms and cancer therapies. We hypothesize that concepts of condensed-matter physics along with the new genomic knowledge and technologies and mechanistic mathematical modeling in conjunction with advances in experimental DNA (Deoxyrinonucleic acid molecule) repair and cell signaling have now provided us with unprecedented opportunities in radiation biophysics to address problems in targeted cancer therapy, and genetic risk estimation in humans. Obviously, one is not dealing with 'low-hanging fruit', but it will be a major scientific achievement if it becomes possible to state, in another decade or so, that we can link mechanistically the stages between the initial radiation-induced DNA damage; in particular, at doses of radiation less than 2 Gy and with structural changes in genomic DNA as a precursor to cell inactivation and/or mutations leading to genetic diseases. The paper presents recent development in the physics of radiation track structure contained in the computer code system KURBUC, in particular for low-energy electrons in the condensed phase of water for which we provide a comprehensive discussion of the dielectric response function approach. The state-of-the-art in the simulation of proton and carbon ion tracks in the Bragg peak region is also presented. The paper presents a critical discussion of the models used for elastic scattering, and the validity of the trajectory approach in low-electron transport. Brief discussions of mechanistic and quantitative aspects of microdosimetry, DNA damage and DNA repair are also included as developed by the authors' work.

Calculating Variations in Biological Effectiveness for a 62 MeV Proton Beam

A biophysical model of radiation-induced cell death and chromosome aberrations [called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA)] was further developed and applied to therapeutic protons. The model assumes a pivotal role of DNA cluster damage, which can lead to clonogenic cell death following three main steps: (i) a DNA “cluster lesion” (CL) produces two independent chromosome fragments; (ii) fragment mis-rejoining within a threshold distance d gives rise to chromosome aberrations; (iii) certain aberration types (dicentrics, rings, and large deletions) lead to clonogenic inactivation. The yield of CLs and the probability, f, that a chromosome fragment remains un-rejoined even if other fragment(s) are present within d, were adjustable parameters. The model, implemented as a MC code providing simulated dose–responses directly comparable with experimental data, was applied to pristine and modulated Bragg peaks of the proton beam used to treat eye melanoma at INFN-LNS in Catania, Italy. Experimental survival curves for AG01522 cells exposed to the Catania beam were reproduced, supporting the model assumptions. Furthermore, cell death and chromosome aberrations at different depths along a spread-out Bragg peak (SOBP) dose profile were predicted. Both endpoints showed an increase along the plateau, and high levels of damage were found also beyond the distal dose fall-off, due to low-energy protons. Cell death and chromosome aberrations were also predicted for V79 cells, in the same irradiation scenario as that used for AG01522 cells. In line with other studies, this work indicated that assuming a constant relative biological effectiveness (RBE) along a proton SOBP may be sub-optimal. Furthermore, it provided qualitative and quantitative evaluations of the dependence of the beam effectiveness on the considered endpoint and dose. More generally, this work represents an example of therapeutic beam characterization avoiding the use of experimental RBE values, which can be source of uncertainties.

Modeling radiation-induced cell death: role of different levels of DNA damage clustering

Some open questions on the mechanisms underlying radiation-induced cell death were addressed by a biophysical model, focusing on DNA damage clustering and its consequences. DNA "cluster lesions" (CLs) were assumed to produce independent chromosome fragments that, if created within a micrometer-scale threshold distance (d), can lead to chromosome aberrations following mis-rejoining; in turn, certain aberrations (dicentrics, rings and large deletions) were assumed to lead to clonogenic cell death. The CL yield and d were the only adjustable parameters. The model, implemented as a Monte Carlo code called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), provided simulated survival curves that were directly compared with experimental data on human and hamster cells exposed to photons, protons, α-particles and heavier ions including carbon and iron. d = 5 μm, independent of radiation quality, and CL yields in the range ~2-20 CLs Gy(-1) cell(-1), depending on particle type and energy, led to good agreement between simulations and data. This supports the hypothesis of a pivotal role of DNA cluster damage at sub-micrometric scale, modulated by chromosome fragment mis-rejoining at micrometric scale. To investigate the features of such critical damage, the CL yields were compared with experimental or theoretical yields of DNA fragments of different sizes, focusing on the base-pair scale (related to the so-called local clustering), the kbp scale ("regional clustering") and the Mbp scale, corresponding to chromatin loops. Interestingly, the CL yields showed better agreement with kbp fragments rather than bp fragments or Mbp fragments; this suggests that also regional clustering, in addition to other clustering levels, may play an important role, possibly due to its relationship with nucleosome organization in the chromatin fiber.

The BIANCA model/code of radiation-induced cell death: application to human cells exposed to different radiation types

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.

A model of radiation-induced cell killing: insights into mechanisms and applications for hadron therapy

A mechanism-based, two-parameter biophysical model of cell killing was developed with the aim of elucidating the mechanisms underlying radiation-induced cell death and predicting cell killing by different radiation types, including protons and carbon ions at energies and doses of interest for cancer therapy. The model assumed that certain chromosome aberrations (dicentrics, rings and large deletions, called "lethal aberrations") lead to clonogenic inactivation, and that aberrations derive from μm-scale misrejoining of chromatin fragments, which in turn are produced by "dirty" double-strand breaks called "cluster lesions" (CLs). The average numbers of CLs per Gy per cell were left as a semi-free parameter and the threshold distance for chromatin-fragment rejoining was defined the second parameter. The model was "translated" into Monte Carlo code and provided simulated survival curves, which were compared with survival data on V79 cells exposed to protons, carbon ions and X rays. The agreement was good between simulations and survival data and supported the assumptions of the model at least for doses up to a few Gy. Dicentrics, rings and large deletions were found to be lethal not only for AG1522 cells exposed to X rays, as already reported by others, but also for V79 cells exposed to protons and carbon ions of different energies. Furthermore, the derived CL yields suggest that the critical DNA lesions leading to clonogenic inactivation are more complex than "clean" DSBs. After initial validation, the model was applied to characterize the particle and LET dependence of proton and carbon cell killing. Consistent with the proton data, the predicted fraction of inactivated cells after 2 Gy protons was 40-50% below 7.7 keV/μm, increased by a factor ∼1.6 between 7.7-30.5 keV/μm, and decreased by a factor ∼1.1 between 30.5-34.6 keV/μm. These LET values correspond to proton energies below a few MeV, which are always present in the distal region of hadron therapy spread-out Bragg peaks (SOBP). Consistent with the carbon data, the predicted fraction of inactivated cells after 2 Gy carbon was 40-50% between 13.7-32.4 keV/μm, it increased by a factor ∼1.7 between 32.4-153.5 keV/μm, and decreased by a factor ∼1.1 between 153.5-339.1 keV/μm. Finally, we applied the model to predict cell death at different depths along a carbon SOBP used for preclinical experiments at HIMAC in Chiba, Japan. The predicted fraction of inactivated cells was found to be roughly constant (less than 10%) along the SOBP, suggesting that this approach may be applied to predict cell killing of therapeutic carbon beams and that, more generally, dicentrics, rings and deletions at the first mitosis may be regarded as a biological dose for these beams. This study advanced our understanding of the mechanisms of radiation-induced cell death and characterized the particle and LET dependence of proton and carbon cell killing along a carbon SOBP. The model does not use RBE values, which can be a source of uncertainty. More generally, this model is a mechanism-based tool that in minutes can predict cell inactivation by protons or carbon ions of a given energy and dose, based on an experimental photon curve and in principle, a single (experimental) survival point for the considered ion type and energy.

From radiation-induced chromosome damage to cell death: modelling basic mechanisms and applications to boron neutron capture therapy

Cell death is a crucial endpoint in radiation-induced biological damage: on one side, cell death is a reference endpoint to characterise the action of radiation in biological targets; on the other side, any cancer therapy aims to kill tumour cells. Starting from Lea's target theory, many models have been proposed to interpret radiation-induced cell killing; after briefly discussing some of these models, in this paper, a mechanistic approach based on an experimentally observed link between chromosome aberrations and cell death was presented. More specifically, a model and a Monte Carlo code originally developed for chromosome aberrations were extended to simulate radiation-induced cell death applying an experimentally observed one-to-one relationship between the average number of 'lethal aberrations' (dicentrics, rings and deletions) per cell and -ln S, S being the fraction of surviving cells. Although such observation was related to X rays, in the present work, the approach was also applied to protons and alpha particles. A good agreement between simulation outcomes and literature data provided a model validation for different radiation types. The same approach was then successfully applied to simulate the survival of cells enriched with boron and irradiated with thermal neutrons at the Triga Mark II reactor in Pavia, to mimic a typical treatment for boron neutron capture therapy.

In vitro and in vivo studies of boron neutron capture therapy: boron uptake/washout and cell death

Boron neutron capture therapy (BNCT) is a binary radiotherapy based on thermal-neutron irradiation of cells enriched with (10)B, which produces α particles and (7)Li ions of short range and high biological effectiveness. The selective uptake of boron by tumor cells is a crucial issue for BNCT, and studies of boron uptake and washout associated with cell survival studies can be of great help in developing clinical applications. In this work, boron uptake and washout were characterized both in vitro for the DHDK12TRb (DHD) rat colon carcinoma cell line and in vivo using rats bearing liver metastases from DHD cells. Despite a remarkable uptake, a large boron release was observed after removal of the boron-enriched medium from in vitro cell cultures. However, analysis of boron washout after rat liver perfusion in vivo did not show a significant boron release, suggesting that organ perfusion does not limit the therapeutic effectiveness of the treatment. The survival of boron-loaded cells exposed to thermal neutrons was also assessed; the results indicated that the removal of extracellular boron does not limit treatment effectiveness if adequate amounts of boron are delivered and if the cells are kept at low temperature. Cell survival was also investigated theoretically using a mechanistic model/Monte Carlo code originally developed for radiation-induced chromosome aberrations and extended here to cell death; good agreement between simulation outcomes and experimental data was obtained.

From DNA radiation damage to cell death: theoretical approaches

Some representative models of radiation-induced cell death, which is a crucial endpoint in radiobiology, were reviewed. The basic assumptions were identified, their consequences on predicted cell survival were analyzed, and the advantages and drawbacks of each approach were outlined. In addition to “historical” approaches such as the Target Theory, the Linear-Quadratic model, the Theory of Dual Radiation Action and Katz' model, the more recent Local Effect Model was discussed, focusing on its application in Carbon-ion hadrontherapy. Furthermore, a mechanistic model developed at the University of Pavia and based on the relationship between cell inactivation and chromosome aberrations was presented, together with recent results; the good agreement between model predictions and literature experimental data on different radiation types (photons, protons, alpha particles, and Carbon ions) supported the idea that asymmetric chromosome aberrations like dicentrics and rings play a fundamental role for cell death. Basing on these results, a reinterpretation of the TDRA was also proposed, identifying the TDRA “sublesions” and “lesions” as clustered DNA double-strand breaks and (lethal) chromosome aberrations, respectively.

First steps towards systems radiation biology studies concerned with DNA and chromosome structure within living cells

For the understanding of radiation action on biological systems like cellular macromolecules (e.g., DNA in its higher structures) a synergistic approach of experiments and quantitative modelling of working hypotheses is necessary. Further on, the influence on calculated results of certain assumptions in such working hypotheses must critically be evaluated. In the present work, this issue is highlighted in two aspects for the case of DNA damage in single cells. First, yields of double-strand breaks and frequency distributions of DNA fragment lengths after ion irradiation were calculated using different assumptions on the DNA target model. Compared to a former target model now a moderate effect due to the inclusion of a spherical chromatin domain model has been found. Second, the influence of assumptions on particular geometric chromosome models on calculated chromosome aberration data is illustrated with two target-modelling approaches for this end point.