Experimental and Monte Carlo characterization of radionuclidic impurities originated from proton irradiation of [18O]H2O in a modern medical cyclotron.

In this work we present a characterization of the radionuclidic impurities originated by proton irradiation of enriched water [18O]H2O in a medical cyclotron through Monte Carlo simulations and experimental measurements. A set of standard samples of enriched water loaded in the cyclotron target cell have been irradiated at 30 µA proton current for 1 h each and, after an appropriate cooling time, measured by HPGe gamma spectrometry. In this way it was possible to study the direct release of radionuclidic impurities from target components as well as the release as a function of target ageing. Previously to experimental measurements, Monte Carlo calculations with the PHITS Code have been carried out to estimate the radionuclides generated within the target components (in particular Havar® foil) with the aim to identify the nuclides expected to be found in the irradiated water due to cell-to-water transmission mechanisms. Comparison between simulations data and experimental measurements by gamma spectrometry showed that only a very small amount of the radionuclides produced in the target window are released in the enriched water through corrosion/erosion effects, while the release decreases with increasing aging of the target.

Radiation resistance of elastomeric O-rings in mixed neutron and gamma fields: Testing methodology and experimental results.

Materials and components employed in the presence of intense neutron and gamma fields are expected to absorb high dose levels that may induce deep modifications of their physical and mechanical properties, possibly causing loss of their function. A protocol for irradiating elastomeric materials in reactor mixed neutron and gamma fields and for testing the evolution of their main mechanical and physical properties with absorbed dose has been developed. Four elastomeric compounds used for vacuum O-rings, one fluoroelastomer polymer (FPM) based and three ethylene propylene diene monomer rubber (EPDM) based, presently available on the market have been selected for the test. One EPDM is rated as radiation resistant in gamma fields, while the other elastomers are general purpose products. Particular care has been devoted to dosimetry calculations, since absorbed dose in neutron fields, unlike pure gamma fields, is strongly dependent on the material composition and, in particular, on the hydrogen content. The products have been tested up to about 2 MGy absorbed dose. The FPM based elastomer, in spite of its lower dose absorption in fast neutron fields, features the largest variations of properties, with a dramatic increase in stiffness and brittleness. Out of the three EPDM based compounds, one shows large and rapid changes in the main mechanical properties, whereas the other two feature more stable behaviors. The performance of the EPDM rated as radiation resistant in pure gamma fields does not appear significantly better than that of the standard product. The predictive capability of the accelerated irradiation tests performed as well as the applicable concepts of threshold of radiation damage is discussed in view of the use of the examined products in the selective production of exotic species facility, now under construction at the Legnaro National Laboratories of the Italian Istituto Nazionale di Fisica Nucleare. It results that a careful account of dose rate effects and oxygen penetration in the material, both during test irradiations and in operating conditions, is needed to obtain reliable predictions.

Characterisation of the secondary neutron field generated by a compact PET cyclotron with MCNP6 and experimental measurements.

The production of the most common used PET radioisotope Fluorine-18 with commercial cyclotrons is obtained from the 18O(p,n)18F nuclear reaction when 18O-enriched water is bombarded with a proton beam. We present the characterization of the secondary neutron field spectra produced by this reaction in different locations around the cyclotron, through a comparison between MCNP6 Monte Carlo simulation results and experimental data obtained with Neutron Activation Analysis (NAA) of thin target foils of different materials.

The origin of neutron biological effectiveness as a function of energy.

The understanding of the impact of radiation quality in early and late responses of biological targets to ionizing radiation exposure necessarily grounds on the results of mechanistic studies starting from physical interactions. This is particularly true when, already at the physical stage, the radiation field is mixed, as it is the case for neutron exposure. Neutron Relative Biological Effectiveness (RBE) is energy dependent, maximal for energies ~1 MeV, varying significantly among different experiments. The aim of this work is to shed light on neutron biological effectiveness as a function of field characteristics, with a comprehensive modeling approach: this brings together transport calculations of neutrons through matter (with the code PHITS) and the predictive power of the biophysical track structure code PARTRAC in terms of DNA damage evaluation. Two different energy dependent neutron RBE models are proposed: the first is phenomenological and based only on the characterization of linear energy transfer on a microscopic scale; the second is purely ab-initio and based on the induction of complex DNA damage. Results for the two models are compared and found in good qualitative agreement with current standards for radiation protection factors, which are agreed upon on the basis of RBE data.

Energy dependence of the complexity of DNA damage induced by carbon ions.

To assess the complexity of DNA damage induced by carbon ions as a function of their energy and LET, 2-Gy irradiations by 100 keV u(-1)-400 MeV u(-1) carbon ions were investigated using the PARTRAC code. The total number of fragments and the yield of fragments of <30 bp were calculated. The authors found a particularly important contribution of DNA fragmentation in the range of <1 kbp for specific energies of <6 MeV u(-1). They also considered the effect of different specific energies with the same LET, i.e. before and after the Bragg peak. As a first step towards a full characterisation of secondary particle production from carbon ions interacting with tissue, a comparison between DNA-damage induction by primary carbon ions and alpha particles resulting from carbon break-up is presented, for specific energies of >1 MeV u(-1).

Neutron flux characterisation of the Pavia TRIGA Mark II research reactor for radiobiological and microdosimetric applications.

Nowadays the Pavia TRIGA reactor is available for national and international collaboration in various research fields. The TRIGA Mark II nuclear research reactor of the Pavia University offers different in- and out-core neutron irradiation channels, each characterised by different neutron spectra. In the last two years a campaign of measurements and simulations has been performed in order to guarantee a better characterisation of these different fluxes and to meet the demands of irradiations that require precise information on these spectra in particular for radiobiological and microdosimetric studies. Experimental data on neutron fluxes have been collected analysing and measuring the gamma activity induced in thin target foils of different materials irradiated in different TRIGA experimental channels. The data on the induced gamma activities have been processed with the SAND II deconvolution code and finally compared with the spectra obtained with Monte Carlo simulations. The comparison between simulated and measured spectra showed a good agreement allowing a more precise characterisation of the neutron spectra and a validation of the adopted method.

Reaction mechanism interplay in determining the biological effectiveness of neutrons as a function of energy.

Neutron relative biological effectiveness (RBE) is found to be energy dependent, being maximal for energies ∼1 MeV. This is reflected in the choice of radiation weighting factors wR for radiation protection purposes. In order to trace back the physical origin of this behaviour, a detailed study of energy deposition processes with their full dependences is necessary. In this work, the Monte Carlo transport code PHITS was used to characterise main secondary products responsible for energy deposition in a 'human-sized' soft tissue spherical phantom, irradiated by monoenergetic neutrons with energies around the maximal RBE/wR. Thereafter, results on the microdosimetric characterisation of secondary protons were used as an input to track structure calculations performed with PARTRAC, thus evaluating the corresponding DNA damage induction. Within the proposed simplified approach, evidence is suggested for a relevant role of secondary protons in inducing the maximal biological effectiveness for 1 MeV neutrons.

Modeling dose deposition and DNA damage due to low-energy ß(-) emitters.

One of the main issues of low-energy internal emitters concerns the very short ranges of the beta particles, versus the dimensions of the biological targets. Depending on the chemical form, the radionuclide may be more concentrated either in the cytoplasm or in the nucleus of the target cell. Consequently, since in most cases conventional dosimetry neglects this issue it may overestimate or underestimate the dose to the nucleus and hence the biological effects. To assess the magnitude of these deviations and to provide a realistic evaluation of the localized energy deposition by low-energy internal emitters, the biophysical track-structure code PARTRAC was used to calculate nuclear doses, DNA damage yields and fragmentation patterns for different localizations of radionuclides in human interphase fibroblasts. The nuclides considered in the simulations were tritium and nickel-63, which emit electrons with average energies of 5.7 (range in water of 0.42 µm) and 17 keV (range of 5 µm), respectively, covering both very short and medium ranges of beta-decay products. The simulation results showed that the largest deviations from the conventional dosimetry occur for inhomogeneously distributed short-range emitters. For uniformly distributed radionuclides selectively in the cytoplasm but excluded from the cell nucleus, the dose in the nucleus is 15% of the average dose in the cell in the case of tritium but 64% for nickel-63. Also, the numbers of double-strand breaks (DSBs) and the distributions of DNA fragments depend on subcellular localization of the radionuclides. In the low- and medium-dose regions investigated here, DSB numbers are proportional to the nuclear dose, with about 50 DSB/Gy for both studied nuclides. In addition, DSB numbers on specific chromosomes depend on the radionuclide localization in the cell as well, with chromosomes located more peripherally in the cell nucleus being more damaged by short-ranged emitters in cytoplasm compared with chromosomes located more centrally. These results illustrate the potential for over- or underestimating the risk associated with low-energy emitters, particularly for tritium intake, when their distribution at subcellular levels is not appropriately considered.

Integration of Monte Carlo simulations with PFGE experimental data yields constant RBE of 2.3 for DNA double-strand break induction by nitrogen ions between 125 and 225 keV/µm LET.

The number of small radiation-induced DNA fragments can be heavily underestimated when determined from measurements of DNA mass fractions by gel electrophoresis, leading to a consequent underestimation of the initial DNA damage induction. In this study we reanalyzed the experimental results for DNA fragmentation and DNA double-strand break (DSB) yields in human fibroblasts irradiated with γ rays and nitrogen ion beams with linear energy transfer (LET) equal to 80, 125, 175 and 225 keV/µm, originally measured by Höglund et al. (Radiat Res 155, 818-825, 2001 and Int J Radiat Biol 76, 539-547, 2000). In that study the authors converted the measured distributions of fragment masses into DNA fragment distributions using mid-range values of the measured fragment length intervals, in particular they assumed fragments with lengths in the interval of 0-48 kbp had the mid-range value of 24 kbp. However, our recent detailed simulations with the Monte Carlo code PARTRAC, while reasonably in agreement with the mass distributions, indicate significantly increased yields of very short fragments by high-LET radiation, so that the actual average fragment lengths, in the interval 0-48 kbp, 2.4 kbp for 225 keV/µm nitrogen ions were much shorter than the assumed mid-range value of 24 kbp. When the measured distributions of fragment masses are converted into fragment distributions using the average fragment lengths calculated by PARTRAC, significantly higher yields of DSB related to short fragments were obtained and resulted in a constant relative biological effectiveness (RBE) for DSB induction yield of 2.3 for nitrogen ions at 125-225 keV/µm LET. The previously reported downward trend of the RBE values over this LET range for DSB induction appears to be an artifact of an inadequate average fragment length in the smallest interval.

Monte Carlo evaluation of DNA fragmentation spectra induced by different radiation qualities.

The PARTRAC code has been developed constantly in the last several years. It is a Monte Carlo code based on an event-by-event description of the interactions taking place between the ionising radiation and liquid water, and in the present version simulates the transport of photons, electrons, protons, helium and heavier ions. This is combined with an atom-by-atom representation of the biological target, i.e. the DNA target model of a diploid human fibroblast in its interphase (genome of 6 Gigabase pairs). DNA damage is produced by the events of energy depositions, either directly, if they occur in the volume occupied by the sugar-phosphate backbone, or indirectly, if this volume is reached by radiation-induced radicals. This requires the determination of the probabilities of occurrence of DNA damage. Experimental data are essential for this determination. However, after the adjustment of the relevant parameters through the comparison of the simulation data with the DNA fragmentation induced by photon irradiation, the code has been used without further parameter adjustments, and the comparison with the fragmentation induced by charged particle beams has validated the code. In this paper, the results obtained for the DNA fragmentation induced by gamma rays and by charged particle beams of various LET are shown, with a particular attention to the production of very small fragments that are not detected in experiments.

Radiation-induced perturbation of cell-to-cell signalling and communication.

The investigation of the bystander phenomena (i.e. the induction of damage in cells not directly traversed by radiation) is strictly related to the study of the mechanisms of intercellular communication and of the perturbative effects of radiation. A new possible way to try to solve the bystander puzzle is through a 'systems radiation biology' approach with the total integration of experimental and theoretical activities. In particular, this contribution will focus on: (1) 'ad hoc' experiments designed to quantify key parameters involved in intercellular signalling (focusing, as a pilot study, on release, decay and internalization of interleukine-6 molecules, their modulation by radiation, and possible differences between in vivo/in vitro behaviour); (2) the implementation and the development of two different modelling approaches: a stochastic model (based on a Monte Carlo code) that takes account of the local mechanisms of release and internalization of signalling molecules (e.g. cytokines) and an analytical model where signal molecules are treated as a population and their temporal behaviour is described by differential equations. This approach provided instruments to investigate the complex phenomena of signal transmission and the role of cell communication to guarantee (maintain) the robustness of the in vitro experimental systems against the effects of perturbations.

A Monte Carlo study of the radiation quality dependence of DNA fragmentation spectra.

We simulated the irradiation of human fibroblasts with gamma rays, protons and helium, carbon and iron ions at a fixed dose of 5 Gy. The simulations were performed with the biophysical Monte Carlo code PARTRAC. From the output of the code, containing in particular the genomic positions of the radiation-induced DNA double-strand breaks (DSBs), we obtained the DNA fragmentation spectra. Very small fragments, in particular those related to "complex lesions" (few tens of base pairs), are probably very important for the late cellular consequences, but their detection is not possible with the common experimental techniques. We paid special attention to the differences among the various ions in the production of these very small fragments; in particular, we compared the fragmentation spectra for ions of the same specific energy and for ions of the same LET (linear energy transfer). As found previously for iron ions, we found that the RBE (relative biological effectiveness) for DSB production was considerably higher than 1 for all high-LET radiations considered. This is at variance with the results obtainable from experimental data, and it is due to the ability to count the contribution of small fragments. It should be noted that for a given LET this RBE decreases with increasing ion charge, due mainly to the increasing mean energy of secondary electrons. A precise quantification of the DNA initial damage can be of great importance for both radiation protection, particularly in open-space long-term manned missions, and hadrontherapy.

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.

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.