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

The FLUKA Code: An Accurate Simulation Tool for Particle Therapy

Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both ⁴He and ¹²C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth–dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features.

The relative contributions of DNA strand breaks, base damage and clustered lesions to the loss of DNA functionality induced by ionizing radiation

The majority of studies on lethal radiobiological damage have focused on double-strand breaks (DSBs), a type of clustered DNA damage and the evaluation of their toxicity, while other types of clustered DNA damage have received much less attention. The main purpose of this study is to evaluate the contribution of different lesions induced by ionizing radiation to the loss of plasmid DNA functionality. We employed a simple model system comprising E. coli transformed with an irradiated plasmid [pGEM-3Zf (-)] to determine the effect of DSBs and other lesions including base damage and clustered lesions on the functionality ("viability") of the plasmid. The yields of γ-radiation-induced single-strand breaks (SSBs) and DSBs were measured by gel electrophoresis. We found that the transformation efficiency decreases with radiation dose, but this decrease cannot be explained by the formation of DSBs. For example, at doses of 500 and 700 Gy, the relative transformation efficiency falls from 100% to 53% and 26%, respectively, while only 5.7% and 9.1% of the plasmids contain a DSB. In addition, it is also unlikely that randomly distributed base lesions could explain the loss of functionality of the plasmid, since cells can repair them efficiently. However, clustered lesions other than DSBs, which are difficult to repair and result in the loss of information on both DNA strands, have the potential to induce the loss of plasmid functionality. We therefore measured the yields of γ-radiation-induced base lesions and cluster damage, which are respectively converted into SSBs and DSBs by the base excision repair enzymes endonuclease III (Nth) and formamidopyrimidine-DNA glycosylase (Fpg). Our data demonstrate that the yield of cluster damage (i.e., lesions that yield DSBs following digestion) is 31 times higher than that of frank DSBs. This finding suggests that frank DSBs make a relatively minor contribution to the loss of DNA functionality induced by ionizing radiation, while other toxic lesions formed at a much higher frequencies than DSBs must be responsible for the loss of plasmid functionality. These lesions may be clustered lesions/locally multiply damaged sites (LMDS), including base damage, SSBs and/or intrastrand and interstrand crosslinks, leading to the loss of vital information in the DNA. Using a mathematical model, we estimate that at least three toxic lesions are required for the inactivation of plasmid functionality, in part because even these complex lesions can be repaired.

Effects of shielding on the induction of 53BP1 foci and micronuclei after Fe ion exposures

High atomic number and high-energy (HZE) particles in deep space are of low abundance but substantially contribute to the biological effects of space radiation. Shielding is so far the most effective way to partially protect astronauts from these highly penetrating particles. However, simulated calculations and measurements have predicted that secondary particles resulting from the shielding of cosmic rays produce a significant fraction of the total dose and dose equivalent. In this study, we investigated the biological effects of secondary radiation with two cell types, and with cells exposed in different phases of the cell cycle, by comparing the biological effects of a 200 MeV/u iron beam with a shielded beam in which the energy of the iron ion beam was decreased from 500 MeV/u to 200 MeV/u with PMMA, polyethylene (PE), or aluminum. We found that beam shielding resulted in increased induction of 53BP1 foci and micronuclei in a cell-type-dependent manner compared with the unshielded 200 MeV/u Fe ion beam. These findings provide experimental proof that the biological effects of secondary particles resulting from the interaction between HZE particles and shielding materials should be considered in shielding design.

Effect of SPE-like proton or photon radiation on the kinetics of mouse peripheral blood cells and radiation biological effectiveness determinations

Exploration missions outside low-Earth orbit are being planned; therefore, it is critical to understand the risk astronauts would be exposed to in the space environment, especially during extravehicular activities (EVAs). Reductions in white blood cell (WBC) numbers can occur as a result of exposure to solar particle event (SPE) radiation. The aim of the present study was to determine the duration of the effects on blood cell numbers from exposure to a single whole-body dose of SPE-like proton radiation or photon radiation as well as to determine the radiation biological effectiveness (RBE) values at those times when radiation exposure causes blood cell numbers to experience the most critical effects when using mice as a model. Our results indicate that both types of radiation cause significant reductions in the numbers of all blood cell types at different times post-irradiation. The RBE values were not significantly different from 1.0. These results indicate that the risk estimations for astronauts from exposure of mice to SPE-like proton radiation are comparable to those previously made for doses of standard reference radiations, suggesting that countermeasures should be developed for the decreases in blood cell counts observed following the exposure of mice to SPE radiation.

The effect of anatomical modeling on space radiation dose estimates: a comparison of doses for NASA phantoms and the 5th, 50th, and 95th percentile male and female astronauts

The National Aeronautics and Space Administration (NASA) performs organ dosimetry and risk assessment for astronauts using model-normalized measurements of the radiation fields encountered in space. To determine the radiation fields in an organ or tissue of interest, particle transport calculations are performed using self-shielding distributions generated with the computer program CAMERA to represent the human body. CAMERA mathematically traces linear rays (or path lengths) through the computerized anatomical man (CAM) phantom, a computational stylized model developed in the early 1970s with organ and body profiles modeled using solid shapes and scaled to represent the body morphometry of the 1950 50th percentile (PCTL) Air Force male. With the increasing use of voxel phantoms in medical and health physics, a conversion from a mathematical-based to a voxel-based ray-tracing algorithm is warranted. In this study, the voxel-based ray tracer (VoBRaT) is introduced to ray trace voxel phantoms using a modified version of the algorithm first proposed by Siddon (1985 Med. Phys. 12 252-5). After validation, VoBRAT is used to evaluate variations in body self-shielding distributions for NASA phantoms and six University of Florida (UF) hybrid phantoms, scaled to represent the 5th, 50th, and 95th PCTL male and female astronaut body morphometries, which have changed considerably since the inception of CAM. These body self-shielding distributions are used to generate organ dose equivalents and effective doses for five commonly evaluated space radiation environments. It is found that dosimetric differences among the phantoms are greatest for soft radiation spectra and light vehicular shielding.

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.

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.