Modelling chromosomal aberration induction by ionising radiation: the influence of interphase chromosome architecture

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.

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.

Biological effectiveness of He-3 and He-4 ion beams for cancer hadrontherapy: a study based on the BIANCA biophysical model

While cancer therapy with protons and C-ions is continuously spreading, in the near future patients will be also treated with He-ions which, in comparison to photons, combine the higher precision of protons with the higher relative biological effectiveness (RBE) of C-ions. Similarly to C-ions, also for He-ions the RBE variation along the beam must be known as precisely as possible, especially for active beam delivery systems. In this framework the BIANCA biophysical model, which has already been applied to calculate the RBE along proton and C-ion beams, was extended to⁴He-ions and, following interface with the FLUKA code, was benchmarked against cell survival data on CHO normal cells and Renca tumour cells irradiated at different positions along therapeutic-like⁴He-ion beams at the Heidelberg Ion-beam Therapy centre, where the first He-ion patient will be treated soon. Very good agreement between simulations and data was obtained, showing that BIANCA can now be used to predict RBE following irradiation with all ion types that are currently used, or will be used soon, for hadrontherapy. Thanks to the development of a reference simulation database describing V79 cell survival for ion and photon irradiation, these predictions can be cell-type specific because analogous databases can be produced, in principle, for any cell line. Furthermore, survival data on CHO cells irradiated by a He-3 beam were reproduced to compare the biophysical properties of He-4 and He-3 beams, which is currently an open question. This comparison showed that, at the same depth, He-4 beams tend to have a higher RBE with respect to He-3 beams, and that this difference is also modulated by the considered physical dose, as well as the cell radiosensitivity. However, at least for the considered cases, no significant difference was found for the ratio between the RBE-weighted dose in the SOBP and that in the entrance plateau.

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.

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.

Chromatin structure and radiation-induced DNA damage: from structural biology to radiobiology

Genomic DNA in eukaryotic cells is basically divided into chromosomes, each consisting of a single huge nucleosomal fiber. It is now clear that chromatin structure and dynamics play a critical role in all processes involved in DNA metabolism, e.g. replication, transcription, repair and recombination. Radiation is a useful tool to study the biological effects of chromatin alterations. Conversely, radiotherapy and radiodiagnosis raise questions about the influence of chromatin integrity on clinical features and secondary effects. This review focuses on the link between DNA damage and chromatin structure at different scales, showing how a comprehensive multiscale vision is required to understand better the effect of radiations on DNA. Clinical aspects related to high- and low-dose of radiation and chromosomal instability will be discussed. At the same time, we will show that the analysis of the radiation-induced DNA damage distribution provides good insight on chromatin structure. Hence, we argue that chromatin "structuralists" and radiobiological "clinicians" would each benefit from more collaboration with the other. We hope that this focused review will help in this regard.

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.

Complex interchanges as a complex function of chromosome organisation

A Monte Carlo technique was developed for biophysical modelling of structural organisation of 46 chromosomes within human lymphocyte interphase nucleus. The technique takes into account: different levels of chromatin organisation, non-random localisation of particular chromosomes and chromosome loci dynamics. All chromosomes in a nucleus were modelled as polymer globules. A dynamic pattern of intra/interchromosomal contacts was simulated to predict radiation-induced chromosomal exchange aberrations (CA). Distance dependence of interaction probability was calculated directly with taking DNA break repair into account. Dose-response for simple and complex CA was calculated to analyse mFISH data for human lymphocytes. Calculated simple CA frequencies fitted the experimental data well. Unexpectedly, complex aberrations were underestimated, despite the dense packaging of chromosome territories within a nucleus. To study sensitivity of dose-response to uncertainty of chromosome organisation, CA were recalculated for the alternative concept of nucleus organisation, the SCD model. The simulation showed the underestimation of complex CA yield as well. The movement of damaged loci from different chromosomes to common repair factories was proposed as an additional mechanism of complex CA formation.

A Model of Chromosome Aberration Induction: Applications to Space Research

A mechanistic model and Monte Carlo code simulating chromosome aberration induction in human lymphocytes is presented. The model is based on the assumption that aberrations arise from clustered DNA lesions and that only the free ends of clustered lesions created in neighboring chromosome territories or in the same territory can join and produce exchanges. The lesions are distributed in the cell nucleus according to the radiation track structure. Interphase chromosome territories are modeled as compact intranuclear regions with volumes proportional to the chromosome DNA contents. Both Giemsa staining and FISH painting can be simulated, and background aberrations can be taken into account. The good agreement with in vitro data provides validation of the model in terms of both the assumptions adopted and the simulation techniques. As an application in the field of space research, the model predictions were compared with aberration yields measured among crew members of long-term missions on board Mir and ISS, assuming an average radiation quality factor of 2.4. The agreement obtained also validated the model for in vivo exposure scenarios and suggested possible applications to the prediction of other relevant aberrations, typically translocations.

Role of DNA organisation and environmental scavenging capacity in the evolution of radiobiological damage: models and simulations

BACKGROUND AND PURPOSE

Theoretical models and Monte Carlo simulations were developed, aimed to investigate the role played by the organisation of interphase DNA and the environmental scavenging capacity conditions in the induction of radiobiological damage.

METHODS

The induction of single- and double-strand breaks by gamma rays impinging on different DNA structures (e.g. linear DNA, SV40 minichromosome and cellular DNA) was simulated as a function of the environment scavenging capacity. Furthermore, yields of chromosome aberrations (CA) induced by gamma rays and light ions were simulated with a purposely developed MC code that explicitly takes into account the DNA higher-order organisation as chromosome territories.

RESULTS AND CONCLUSIONS

Simulations performed with the PARTRAC code allowed quantification of the dependence of dsb and ssb both on the target structure, and on the scavenging capacity. The results relative to CA showed the importance of DNA damage complexity (nanometre scale) and interphase chromosome domains (micrometre scale) in the process of aberration formation. Very good agreement was found between the model predictions on ssb, dsb and CA and available experimental data.

A model of chromosome aberration induction and chronic myeloid leukaemia incidence at low doses

Some chromosome aberration types, generally translocations, are correlated with specific cancers. An example is provided by chronic myeloid leukemia (CML) cells, most of which carry a translocation involving the ABL gene on chromosome 9 and the BCR gene on chromosome 22. The hypothesis of a causal relationship between CML and the chimeric protein product of the BCR-ABL translocation has recently received strong support. In this framework, a mechanistic model and Monte-Carlo code simulating radiation-induced chromosome aberrations in human lymphocytes will be presented. The current version of the model can predict dose-response curves for the main aberration types following acute irradiation with gamma rays and light ions of different energies. The model is based on the assumption that only clustered DNA lesions can lead to aberrations and that only lesion free ends in neighbouring chromosome territories can join and form exchanges. Such lesions are distributed within the cell nucleus according to the radiation track structure, i.e. randomly for low-LET radiation and along straight lines for high-LET light ions. Interphase chromosome territories are explicitly simulated and background aberrations are taken into account. Very good agreement was found with experimental data taken from the literature that provided a further validation of the model. As an application, yields of BCR-ABL translocations were calculated. Preliminary results led to a CML induction dose-response that is approximately quadratic below 0.1 Gy and essentially linear at higher doses up to 1 Gy. The numerical values obtained for the probability of CML induction are consistent with values obtained by other groups with different approaches.

Comparing DNA Damage-Processing Pathways by Computer Analysis of Chromosome Painting Data

Chromosome aberrations are large-scale illegitimate rearrangements of the genome. They are indicative of DNA damage and informative about damage processing pathways. Despite extensive investigations over many years, the mechanisms underlying aberration formation remain controversial. New experimental assays such as multiplex fluorescent in situ hybridyzation (mFISH) allow combinatorial "painting" of chromosomes and are promising for elucidating aberration formation mechanisms. Recently observed mFISH aberration patterns are so complex that computer and graph-theoretical methods are needed for their full analysis. An important part of the analysis is decomposing a chromosome rearrangement process into "cycles." A cycle of order n, characterized formally by the cyclic graph with 2n vertices, indicates that n chromatin breaks take part in a single irreducible reaction. We here describe algorithms for computing cycle structures from experimentally observed or computer-simulated mFISH aberration patterns. We show that analyzing cycles quantitatively can distinguish between different aberration formation mechanisms. In particular, we show that homology-based mechanisms do not generate the large number of complex aberrations, involving higher-order cycles, observed in irradiated human lymphocytes.

Models of chromosome aberration induction: an example based on radiation track structure

A few examples of models of chromosome aberration induction are summarised and discussed on the basis of the three main theories of aberration formation, that is "breakage-and-reunion", "exchange" and "one-hit". A model and code developed at the Universities of Milan and Pavia is then presented in detail. The model provides dose-response curves for different aberration types (dicentrics, translocations, rings, complex exchanges and deletions) induced in human lymphocytes by gamma rays, protons and alpha particles of different energies, both as monochromatic fields and as mixed fields. The main assumptions are that only clustered - and thus severe - DNA breaks ("Complex Lesions", CL) can participate in the production of aberrations, and that only break free ends in neighbouring chromosome territories can interact and form exchanges. The yields of CLs induced by the various radiation types of interest are taken from a previous modelling work. These lesions are distributed within a sphere representing the cell nucleus according to the radiation track structure, e.g. randomly for gamma rays and along straight lines for light ions. Interphase chromosome territories are explicitly simulated and configurations are obtained in which each chromosome occupies an intranuclear domain with volume proportional to its DNA content. In order to allow direct comparisons with experimental data, small fragments can be neglected since usually they cannot be detected in experiments. The presence of a background level of aberrations is also taken into account. The results of the simulations are in good agreement with experimental dose-response curves available in the literature, that provides a validation of the model both in terms of the adopted assumptions and in terms of the simulation techniques. To address the question of "true" incompleteness, simulations were also run in which all fragments were assumed to be visible.

Quantitative analysis of radiation-induced chromosome aberrations

We review chromosome aberration modeling and its applications, especially to biodosimetry and to characterizing chromosome geometry. Standard results on aberration formation pathways, randomness, dose-response, proximity effects, transmissibility, kinetics, and relations to other radiobiological endpoints are summarized. We also outline recent work on graph-theoretical descriptions of aberrations, Monte-Carlo computer simulations of aberration spectra, software for quantifying aberration complexity, and systematic links of apparently incomplete with complete or truly incomplete aberrations.

Role of shielding in modulating the effects of solar particle events: Monte Carlo calculation of absorbed dose and DNA complex lesions in different organs

Distributions of absorbed dose and DNA clustered damage yields in various organs and tissues following the October 1989 solar particle event (SPE) were calculated by coupling the FLUKA Monte Carlo transport code with two anthropomorphic phantoms (a mathematical model and a voxel model), with the main aim of quantifying the role of the shielding features in modulating organ doses. The phantoms, which were assumed to be in deep space, were inserted into a shielding box of variable thickness and material and were irradiated with the proton spectra of the October 1989 event. Average numbers of DNA lesions per cell in different organs were calculated by adopting a technique already tested in previous works, consisting of integrating into "condensed-history" Monte Carlo transport codes--such as FLUKA--yields of radiobiological damage, either calculated with "event-by-event" track structure simulations, or taken from experimental works available in the literature. More specifically, the yields of "Complex Lesions" (or "CL", defined and calculated as a clustered DNA damage in a previous work) per unit dose and DNA mass (CL Gy-1 Da-1) due to the various beam components, including those derived from nuclear interactions with the shielding and the human body, were integrated in FLUKA. This provided spatial distributions of CL/cell yields in different organs, as well as distributions of absorbed doses. The contributions of primary protons and secondary hadrons were calculated separately, and the simulations were repeated for values of Al shielding thickness ranging between 1 and 20 g/cm2. Slight differences were found between the two phantom types. Skin and eye lenses were found to receive larger doses with respect to internal organs; however, shielding was more effective for skin and lenses. Secondary particles arising from nuclear interactions were found to have a minor role, although their relative contribution was found to be larger for the Complex Lesions than for the absorbed dose, due to their higher LET and thus higher biological effectiveness.

Chromosome aberrations as biomarkers of radiation exposure: modelling basic mechanisms

The space radiation environment is a mixed field consisting of different particles having different energies, including high charge and energy (HZE) ions. Conventional measurements of absorbed doses may not be sufficient to completely characterise the radiation field and perform reliable estimates of health risks. Biological dosimetry, based on the observation of specific radiation-induced endpoints (typically chromosome aberrations), can be a helpful approach in case of monitored exposure to space radiation or other mixed fields, as well as in case of accidental exposure. Furthermore, various ratios of aberrations (e.g. dicentric chromosomes to centric rings and complex exchanges to simple exchanges) have been suggested as possible fingerprints of radiation quality, although all of them have been subjected to some criticisms. In this context a mechanistic model and a Monte Carlo code for the simulation of chromosome aberration induction were developed. The model, able to provide dose-responses for different aberrations (e.g. dicentrics, rings, fragments, translocations, insertions and other complex exchanges), was further developed to assess the dependence of various ratios of aberrations on radiation quality. The predictions of the model were compared with available data, whose experimental conditions were faithfully reproduced. Particular attention was devoted to the scoring criteria adopted in different laboratories and to possible biases introduced by interphase death and mitotic delay. This latter aspect was investigated by taking into account both metaphase data and data obtained with Premature Chromosome Condensation (PCC).

Modelling the influence of shielding on physical and biological organ doses

Distributions of "physical" and "biological" dose in different organs were calculated by coupling the FLUKA MC transport code with a geometrical human phantom inserted into a shielding box of variable shape, thickness and material. While the expression "physical dose" refers to the amount of deposited energy per unit mass (in Gy), "biological dose" was modelled with "Complex Lesions" (CL), clustered DNA strand breaks calculated in a previous work based on "event-by-event" track-structure simulations. The yields of complex lesions per cell and per unit dose were calculated for different radiation types and energies, and integrated into a version of FLUKA modified for this purpose, allowing us to estimate the effects of mixed fields. As an initial test simulation, the phantom was inserted into an aluminium parallelepiped and was isotropically irradiated with 500 MeV protons. Dose distributions were calculated for different values of the shielding thickness. The results were found to be organ-dependent. In most organs, with increasing shielding thickness the contribution of primary protons showed an initial flat region followed by a gradual decrease, whereas secondary particles showed an initial increase followed by a decrease at large thickness values. Secondary particles were found to provide a substantial contribution, especially to the biological dose. In particular, the decrease of their contribution occurred at larger depths than for primary protons. In addition, their contribution to biological dose was generally greater than that of primary protons.

Using graph theory to describe and model chromosome aberrations

A comprehensive description of chromosome aberrations is introduced that is suitable for all cytogenetic protocols (e.g. solid staining, banding, FISH, mFISH, SKY, bar coding) and for mathematical analyses. "Aberration multigraphs" systematically characterize and interrelate three basic aberration elements: (1) the initial configuration of chromosome breaks; (2) the exchange process, whose cycle structure helps to describe aberration complexity; and (3) the final configuration of rearranged chromosomes, which determines the observed pattern but may contain cryptic misrejoinings in addition. New aberration classification methods and a far-reaching generalization of mPAINT descriptors, applicable to any protocol, emerge. The difficult problem of trying to infer actual exchange processes from cytogenetically observed final patterns is analyzed using computer algorithms, adaptations of known theorems on cubic graphs, and some new graph-theoretical constructs. Results include the following: (1) For a painting protocol, unambiguously inferring the occurrence of a high-order cycle requires a corresponding number of different colors; (2) cycle structure can be computed by a simple trick directly from mPAINT descriptors if the initial configuration has no more than one break per homologue pair; and (3) higher-order cycles are more frequent than the obligate cycle structure specifies. Aberration multigraphs are a powerful new way to describe, classify and quantitatively analyze radiation-induced chromosome aberrations. They pinpoint (but do not eliminate) the problem that, with present cytogenetic techniques, one observed pattern corresponds to many possible initial configurations and exchange processes.

Radiation-induced chromosome aberrations: insights gained from biophysical modeling

Enzymatic misrepair of ionizing-radiation-induced DNA damage can produce large-scale rearrangements of the genome, such as translocations and dicentrics. These and other chromosome exchange aberrations can cause major phenotypic alterations, including cell death, mutation and neoplasia. Exchange formation requires that two (or more) genomic loci come together spatially. Consequently, the surprisingly rich aberration spectra uncovered by recently developed techniques, when combined with biophysically based computer modeling, help characterize large-scale chromatin architecture in the interphase nucleus. Most results are consistent with a picture whereby chromosomes are mainly confined to territories, chromatin motion is limited, and interchromosomal interactions involve mainly territory surfaces. Aberration spectra and modeling also help characterize DNA repair/misrepair mechanisms. Quantitative results for mammalian cells are best described by a breakage-and-reunion model, suggesting that the dominant recombinational mechanism during the G(0)/G(1) phase of the cell cycle is non-homologous end-joining of radiogenic DNA double strand breaks. In turn, better mechanistic and quantitative understanding of aberration formation gives new insights into health-related applications.