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.

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.

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 to4He-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-like4He-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.

Novelty Improves the Formation and Persistence of Memory in a Naturalistic School Scenario.

One of the top challenges in education and neuroscience consists in translating laboratory results into strategies to improve learning and memory in teaching environments. In that sense, during the last two decades, researchers have discovered specific temporal windows around learning, during which the intervention with some experiences induces modulatory effects on the formation and/or persistence of memory. Based on these results, the aim of the present study was to design a specific strategy to improve the memory of students in a high-school scenario, by assessing the effect of a novel situation experienced close to learning. We found that the long-term memory about a geometrical figure was more precise in the group of students that faced a novel situation 1 h before or after learning the figure than the control group of students who did not face the novelty. This enhancement was probably triggered by processes acting on memory formation mechanisms that remained evident 45 days after learning, indicating that the improvement was sustained over time. In addition, our results showed that novelty no longer improved the memory if it was experienced 4 h before or after learning. However, far beyond this window of efficacy, when it was faced around 10 h after learning, the novel experience improved the memory persistence tested 7 days later. In summary, our findings characterized different temporal windows of the effectiveness of novelty acting on memory processing, providing a simple and inexpensive strategy that could be used to improve memory formation and persistence in high-school students.

DNA fragmentation induced in human fibroblasts by 56Fe ions: experimental data and Monte Carlo simulations.

We studied the DNA fragmentation induced in human fibroblasts by iron-ion beams of two different energies: 115 MeV/nucleon and 414 MeV/nucleon. Experimental data were obtained in the fragment size range 1-5700 kbp; Monte Carlo simulations were performed with the PARTRAC code; data analysis was also performed through the Generalized Broken Stick (GBS) model. The comparison between experimental and simulated data for the number of fragments produced in two different size ranges, 1-23 kbp and 23-5700 kbp, gives a satisfactory agreement for both radiation qualities. The Monte Carlo simulations also allow the counting of fragments outside the experimental range: The number of fragments smaller than 1 kbp is large for both beams, although with a strong difference between the two cases. As a consequence, we can compute different RBEs depending on the size range considered for the fragment counting. The PARTRAC evaluation takes into account fragments of all sizes, while the evaluation from the experimental data considers only the fragments in the range of 1-5700 kbp. When the PARTRAC evaluation is restricted to this range, the agreement between experimental and computed RBE values is again good. When fragments smaller than 1 kbp are also considered, the RBE increases considerably, since gamma rays produce a small number of such fragments. The analysis performed with the GBS model proved to be quite sensitive to showing, with a phenomenological single parameter, variations in double-strand break (DSB) correlation.

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.

A New Standard DNA Damage (SDD) Data Format.

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 understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Exams at classroom have bidirectional effects on the long-term memory of an unrelated graphical task.

The influence of a given event on long-term memory formation of another one has been a relevant topic of study in the neuroscience field in recent years. Students at school learn contents which are usually tested in exam format. However, exam elevates the arousal state of the students acting as a mild stressor that could influence another memory formation ongoing process. Thus, in this study we examine in high school students the effect of exams on long-term retention of unrelated information, learned at different times before or after the exams. Our results show that exams are not innocuous and that they could improve or reduce the retention of temporarily associated content. These effects did not show gender differences. Our findings should alert teachers about the side effects of exams on the learning of other content within the same school day.

Modelling radiation-induced bystander effect and cellular communication.

In the last 10 years evidence has accumulated on the so-called radiation-induced 'non-targeted effects' and in particular on bystander effects, consisting of damage induction in non-irradiated cells most likely following the release of soluble factors by the irradiated ones. These phenomena were observed for different biological endpoints, both lethal and non-lethal for the cell. Although the underlying mechanisms are largely unknown, it is now widely recognised that two types of cellular communication (i.e. via gap junctions and/or release of molecular messengers into the extracellular environment) play a pivotal role. Furthermore, the effects can be significantly modulated by parameters such as cell type and cell-cycle stage, cell density, time after irradiation etc. Theoretical models and simulation codes can be of help to improve our knowledge of the mechanisms, as well as to investigate the possible role of these effects in determining deviations from the linear relationship between dose and risk which is generally applied in radiation protection. In this paper three models, including an approach under development at the University of Pavia, will be presented in detail. The focus will be on the various adopted assumptions, together with their implications in terms of non-targeted radiobiological damage and, more generally, low-dose radiation risk. Comparisons with experimental data will also be discussed.

Role of DNA/chromatin organisation and scavenging capacity in USX- and proton- induced DNA damage.

DNA higher-order structures and (non-histonic) *;OH radical scavengers have well known protective effects in the induction of single- and double-strand breaks by ionising radiation. In a previous work, such protective roles have been quantified for gamma radiation (Valota et al., Int. J. Radiat. Biol. 79, 2003). As a starting base for the simulations, we used the PARTRAC Monte Carlo code, developed within a collaboration involving the University of Pavia and the GSF institute. The code can reproduce the track structure of photons, electrons, protons and heavier ions in liquid water, and it can simulate the DNA content of a human cell at different organisation levels, based on an atom-by-atom approach. In this work we extended the calculations to Ultra-Soft X rays (USX) and protons, separately analysing the effects of different radiation types on various DNA structures (i.e. linear DNA, SV40 'minichromosomes' and compact chromatin) as a function of the *OH scavenging capacity (SC). Both for USX and protons, the calculated damage yields decreased by increasing the SC for the three considered target types. Such decrease can be ascribed to the competition between the reactions *OH-DNA and *OH-scavenger, which becomes more and more likely by increasing the SC. Furthermore, linear DNA was found to be more radiosensitive than SV40 'minichromosomes', which in turn were more radiosensitive than compact chromatin, which is protected by histones. Comparisons with experimental data by Fulford et al. (Int. J. Radiat. Biol. 77, 2001) relative to USX irradiation showed very good agreement. The dependence of the modulating role played by DNA organisation and scavenging capacity on radiation quality is presented and discussed.

Gamma ray-induced bystander effect in tumour glioblastoma cells: a specific study on cell survival, cytokine release and cytokine receptors.

Recent experimental evidence has challenged the paradigm according to which radiation traversal through the nucleus of a cell is a prerequisite for producing genetic changes or biological responses. Thus, unexposed cells in the vicinity of directly irradiated cells or recipient cells of medium from irradiated cultures can also be affected. The aim of the present study was to evaluate, by means of the medium transfer technique, whether interleukin-8 and its receptor (CXCR1) may play a role in the bystander effect after gamma irradiation of T98G cells in vitro. In fact the cell specificity in inducing the bystander effect and in receiving the secreted signals that has been described suggests that not only the ability to release the cytokines but also the receptor profiles are likely to modulate the cell responses and the final outcome. The dose and time dependence of the cytokine release into the medium, quantified using an enzyme linked immunosorbent assay, showed that radiation causes alteration in the release of interleukin-8 from exposed cells in a dose-independent but time-dependent manner. The relative receptor expression was also affected in exposed and bystander cells.

Human exposure to space radiation: role of primary and secondary particles.

Human exposure to space radiation implies two kinds of risk, both stochastic and deterministic. Shielding optimisation therefore represents a crucial goal for long-term missions, especially in deep space. In this context, the use of radiation transport codes coupled with anthropomorphic phantoms allows to simulate typical radiation exposures for astronauts behind different shielding, and to calculate doses to different organs. In this work, the FLUKA Monte Carlo code and two phantoms, a mathematical model and a voxel model, were used, taking the Galactic Cosmic Rays (GCR) spectra from the model of Badhwar and O'Neill. The time integral spectral proton fluence of the August 1972 Solar Particle Event (SPE) was represented by an exponential function. For each aluminium shield thickness, besides total doses the contributions from primary and secondary particles for different organs and tissues were calculated separately. More specifically, organ-averaged absorbed doses, dose equivalents and a form of 'biological dose', defined on the basis of initial (clustered) DNA damage, were calculated. As expected, the SPE doses dramatically decreased with increasing shielding, and doses in internal organs were lower than in skin. The contribution of secondary particles to SPE doses was almost negligible; however it is of note that, at high shielding (10 g cm(-2)), most of the secondaries are neutrons. GCR organ doses remained roughly constant with increasing Al shielding. In contrast to SPE results, for the case of cosmic rays, secondary particles accounted for a significant fraction of the total dose.

DNA DSB induced in human cells by charged particles and gamma rays: experimental results and theoretical approaches.

PURPOSE: To quantify the role played by radiation track structure and background fragments in modulating DNA fragmentation in human cells exposed to gamma-rays and light ions. MATERIALS AND METHODS: Human fibroblasts were exposed in vitro to different doses (in the range from 40 - 200 Gy) of (60)Co gamma-rays and 0.84 MeV protons (Linear Energy Transfer, LET, in tissue 28.5 keV/microm). The resulting DNA fragments were scored under two electrophoretic conditions, in order to optimize separation in the size ranges 0.023 - 1.0 Mbp and 1.0 - 5.7 Mbp. In parallel, DNA fragmentation was simulated both with a phenomenological approach based on the "generalized broken-stick" model, and with a mechanistic approach based on the PARTRAC (acronym of PARticle TRACk) Monte Carlo code (1.32 MeV photons were used for the simulation of (60)Co gamma-rays). RESULTS: For both gamma-rays and protons, the experimental dose response in the range 0.023 - 5.7 Mbp could be approximated as a straight line, the slope of which provided a yield of (5.3 +/- 0.4) x 10(-9) Gy(-1) bp(-1) for gamma-rays and (7.1 +/- 0.6) x 10(-9) Gy(-1) bp(-1) for protons, leading to a Relative Biological Effectiveness (RBE) of 1.3 +/- 0.2. From both theoretical analyses it appeared that, while gamma-ray data were consistent with double-strand breaks (DSB) random induction, protons at low doses showed significant deviation from randomness, implying enhanced production of small fragments in the low molecular weight part of the experimental range. The theoretical analysis of fragment production was then extended to ranges where data were not available, i.e. to fragments larger than 5.7 Mbp and smaller than 23 kbp. The main outcome was that small fragments (<23 kbp) are produced almost exclusively via non-random processes, since their number is considerably higher than that produced by a random insertion of DSB. Furthermore, for protons the number of these small fragments is a significant fraction (about 20%) of the total number of fragments; these fragments remain undetected in these experiments. Calculations for 3.3 MeV alpha particle irradiation (for which no experimental data were available) were performed to further investigate the role of fragments smaller than 23 kbp; in this case, besides the non-random character of their production, their number resulted to be at least as much as half of the total number of fragments. CONCLUSION: Comparison between experimental data and two different theoretical approaches provided further support to the hypothesis of an important role of track structure in modulating DNA damage. According to the theoretical approaches, non-randomness of fragment production was found for proton irradiation for the smaller fragments in the experimental size range and, in a significantly larger extent, for fragments of size less than 23 kbp, both for protons and alpha particles.

The FLUKA code: new developments and application to 1 GeV/n iron beams.

The modeling of ion transport and interactions in matter is a subject of growing interest, driven by the continuous increase of possible application fields. These include hadron therapy, dosimetry, and space missions, but there are also several issues involving fundamental research, accelerator physics, and cosmic ray physics, where a reliable description of heavy ion induced cascades is important. In the present work, the capabilities of the FLUKA code for ion beams will be briefly recalled and some recent developments presented. Applications of the code to the simulation of therapeutic carbon, nitrogen and oxygen ion beams, and of iron beams, which are of direct interest for space mission related experiments, will be also presented together with interesting consideration relative to the evaluation of dosimetric quantities. Both applications involve ion beams in the AGeV range.

The application of FLUKA to dosimetry and radiation therapy.

The FLUKA Monte Carlo code has been evolving over the last several decades and is now widely used for radiation shielding calculations. In order to facilitate the use of FLUKA in dosimetry and therapy applications, supporting software has been developed to allow the direct conversion of the output files from standard CT-scans directly into a voxel geometry for transport within FLUKA. Since the CT-scan information essentially contains only the electron density information over the scanned volume, one needs the specific compositions for each voxel individually. We present here the results of a simple algorithm to assign tissues in the human body to one of four categories: soft-tissue, hard-bone, trabecular-bone and porous-lung. In addition, we explore the problem of the pathlength distributions in porous media such as trabecular bone. A mechanism will be implemented within FLUKA to allow for variable multipal fixed density materials to accommodate the pathlength distributions discovered.

The role of DNA cluster damage and chromosome aberrations in radiation-induced cell killing: a theoretical approach.

The role played by DNA cluster damage and chromosome aberrations in radiation-induced cell killing was investigated, assuming that certain chromosome aberrations (dicentrics, rings and large deletions, or 'lethal aberrations') lead to clonogenic inactivation and that chromosome aberrations are due to micrometre-scale rejoining of chromosome fragments derived from DNA cluster lesions (CLs). The CL yield and the threshold distance governing fragment rejoining were left as model parameters. The model, implemented as a Monte Carlo code called BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations), provided simulated survival curves that were compared with survival data on AG1522 and V79 cells exposed to different radiation types, including heavy ions. The agreement between simulation outcomes and experimental data suggests that lethal aberrations are likely to play an important role in cell killing not only for AG1522 cells exposed to X rays, as already reported by others, but also for other radiation types and other cells. Furthermore, the results are consistent with the hypothesis that the critical DNA lesions leading to cell death and chromosome aberrations are double-strand break clusters (possibly involving the ∼1000-10 000 bp scale) and that the effects of such clusters are modulated by micrometre-scale proximity effects during DNA damage processing.

Evaluation of the dose enhancement of combined ¹°B + ¹57Gd neutron capture therapy (NCT).

An innovative molecule, GdBLDL, for boron neutron capture therapy (BNCT) has been developed and its effectiveness as a BNCT carrier is currently under evaluation using in vivo experiments on small animal tumour models. The molecule contains both (10)B (the most commonly used NCT agent) and (157)Gd nuclei. (157)Gd is the second most studied element to perform NCT, mainly thanks to its high cross section for the capture of low-energy neutrons. The main drawback of (157)Gd neutron capture reaction is the very short range and low-energy secondary charged particles (Auger electrons), which requires (157)Gd to be very close to the cellular DNA to have an appreciable biological effect. Treatment doses were calculated by Monte Carlo simulations to ensure the optimised tumour irradiation and the sparing of the healthy organs of the irradiated animals. The enhancement of the absorbed dose due to the simultaneous presence of (10)B and (157)Gd in the experimental set-up was calculated and the advantage introduced by the presence of (157)Gd was discussed.

Comparative study of the radiobiological effects induced on adherent vs suspended cells by BNCT, neutrons and gamma rays treatments.

The present work is part of a preclinical in vitro study to assess the efficacy of BNCT applied to liver or lung coloncarcinoma metastases and to limb osteosarcoma. Adherent growing cell lines can be irradiated as adherent to the culture flasks or as cell suspensions, differences in radio-sensitivity of the two modalities of radiation exposure have been investigated. Dose related cell survival and cell cycle perturbation results evidenced that the radiosensitivity of adherent cells is higher than that of the suspended ones.

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.

Contractile response of peritubular myoid cells to prostaglandin F2alpha.

Prostaglandin (PG) F2alpha, a well known agonist of smooth muscle, is produced in the male gonad. We have investigated whether PG F2alpha stimulates seminiferous tubule contractility through direct action on peritubular myoid cells. Myoid cells from prepubertal rats were highly purified through Percoll density gradient and cultured in vitro. Stimulation with PG F2alpha was observed to induce: (i) rapid and dose-dependent production of inositol phosphates; (ii) mobilization of Ca2+ from intracellular stores and (iii) cell contraction. Moreover, at a concentration of 10 microM the agonist was found to induce immediate contractile response of peritubular tissue in freshly explanted tubular fragments from both young and adult rats; the explants were examined in whole-mount preparations and the peritubular myoid cell layer was identified by selective staining for alkaline phosphatase activity. Our observations demonstrate that myoid cells are a direct target for PG F2alpha and suggest a role of the eicosanoid in the intragonadal control of seminiferous tubule contractility.