Energy deposition in small cylindrical targets by monoenergetic electrons

Application of a simple DNA damage model developed for electrons to proton irradiation

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keV μm−1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Impact of DNA Geometry and Scoring on Monte Carlo Track-Structure Simulations of Initial Radiation-Induced Damage

Track structure Monte Carlo simulations are a useful tool to investigate the damage induced to DNA by ionizing radiation. These simulations usually rely on simplified geometrical representations of the DNA subcomponents. DNA damage is determined by the physical and physicochemical processes occurring within these volumes. In particular, damage to the DNA backbone is generally assumed to result in strand breaks. DNA damage can be categorized as direct (ionization of an atom part of the DNA molecule) or indirect (damage from reactive chemical species following water radiolysis). We also consider quasi-direct effects, i.e., damage originated by charge transfers after ionization of the hydration shell surrounding the DNA. DNA geometries are needed to account for the damage induced by ionizing radiation, and different geometry models can be used for speed or accuracy reasons. In this work, we use the Monte Carlo track structure tool TOPAS-nBio, built on top of Geant4-DNA, for simulation at the nanometer scale to evaluate differences among three DNA geometrical models in an entire cell nucleus, including a sphere/spheroid model specifically designed for this work. In addition to strand breaks, we explicitly consider the direct, quasi-direct, and indirect damage induced to DNA base moieties. We use results from the literature to determine the best values for the relevant parameters. For example, the proportion of hydroxyl radical reactions between base moieties was 80%, and between backbone, moieties was 20%, the proportion of radical attacks leading to a strand break was 11%, and the expected ratio of base damages and strand breaks was 2.5-3. Our results show that failure to update parameters for new geometric models can lead to significant differences in predicted damage yields.

A double-strand-break model for the relative biological effectiveness of electrons based on ionization clustering

BACKGROUND

The effectiveness of ionizing radiation regarding DNA damage induction depends on its spatial energy deposition pattern. For electrons an increased effectiveness is observed at low kinetic energies due to the enhanced density of energy deposition events at electron track ends.

PURPOSE

A model is presented, which enables the calculation of the double-strand-break (DSB) yield and the relative biological effectiveness (RBE) for DSB induction of electrons.

METHODS

The model applies the mean free path between two ionizations and the assumption that two ionizations within a certain threshold distance are necessary to potentially lead to a DSB. Next to an expression for the electron RBE according to its common definition, a local RBE is determined, which describes the electrons' local effectiveness at a defined point on their track.

RESULTS

This local RBE allows a better understanding of microscopic processes resulting from radiation and can be used, for instance, to describe the mean effectiveness of the mixed electron radiation field as a function of the radial distance to the center of an ion track.

CONCLUSIONS

The presented model reflects the experimentally observed increased effectiveness of low-energetic electrons. It will be used in a future work to improve RBE predictions for ions performed with the local effect model.

Effect of spatial distribution of boron and oxygen concentration on DNA damage induced from boron neutron capture therapy using Monte Carlo simulations

PURPOSE

This paper aims to investigate how the spatial distribution of boron in cells and oxygen concentration affect the DNA damage induced by charged particles in boron neutron capture therapy (BNCT) by Monte Carlo simulations, and further to evaluate the relative biological effectiveness (RBE) of DNA double-strand breaks (DSBs) induction.

MATERIALS AND METHODS

The kinetic energy spectra of α, ⁷Li particles in BNCT arriving at the nucleus surface were obtained from GEANT4 (Geant4 10.05.p01). The DNA damage caused by BNCT was then evaluated using MCDS (MCDS 3.10A).

RESULTS

When α or ⁷Li particles were distributed in the cytomembrane or cytoplasm, the difference in DNA damage of the same types was less than 0.5%. Taking the 137Cs photons as the reference radiation, when the oxygen concentration varied from 0% to 50%, the RBE of 0.54MeV protons and recoil protons varied from 5 to 2, whereas it decreased from 10 to 3 for α or ⁷Li particles.

CONCLUSION

The RBE of DSB induction all charged particles in BNCT decreased with the increase of oxygen concentration. This work indicated that the RBE of different radiation particles of BNCT might be affected by many factors, which should be paid attention to in theoretical research or clinical application.

Simulation of Proton-Induced DNA Damage Patterns Using an Improved Clustering Algorithm

Simulations of deoxyribonucleic acid (DNA) molecular damage use the traversal algorithm that has the disadvantages of being time-consuming, slowly converging, and requiring high-performance computer clusters. This work presents an improved version of the algorithm, "density-based spatial clustering of applications with noise" (DBSCAN), using a KD-tree approach to find neighbors of each point for calculating clustered DNA damage. The resulting algorithm considers the spatial distributions for sites of energy deposition and hydroxyl radical attack, yielding the statistical probability of (single and double) DNA strand breaks. This work achieves high accuracy and high speed at calculating clustered DNA damage that has been induced by proton treatment at the molecular level while running on an i7 quad-core CPU. The simulations focus on the indirect effect generated by hydroxyl radical attack on DNA. The obtained results are consistent with those of other published experiments and simulations. Due to the array of chemical processes triggered by proton treatment, it is possible to predict the effects that different track structures of various energy protons produce on eliciting direct and indirect damage of DNA.

A Simplified Cluster Analysis of Electron Track Structure for Estimating Complex DNA Damage Yields

Complex DNA damage, defined as at least two vicinal lesions within 10-20 base pairs (bp), induced after exposure to ionizing radiation, is recognized as fatal damage to human tissue. Due to the difficulty of directly measuring the aggregation of DNA damage at the nano-meter scale, many cluster analyses of inelastic interactions based on Monte Carlo simulation for radiation track structure in liquid water have been conducted to evaluate DNA damage. Meanwhile, the experimental technique to detect complex DNA damage has evolved in recent decades, so both approaches with simulation and experiment get used for investigating complex DNA damage. During this study, we propose a simplified cluster analysis of ionization and electronic excitation events within 10 bp based on track structure for estimating complex DNA damage yields for electron and X-ray irradiations. We then compare the computational results with the experimental complex DNA damage coupled with base damage (BD) measured by enzymatic cleavage and atomic force microscopy (AFM). The computational results agree well with experimental fractions of complex damage yields, i.e., single and double strand breaks (SSBs, DSBs) and complex BD, when the yield ratio of BD/SSB is assumed to be 1.3. Considering the comparison of complex DSB yields, i.e., DSB + BD and DSB + 2BD, between simulation and experimental data, we find that the aggregation degree of the events along electron tracks reflects the complexity of induced DNA damage, showing 43.5% of DSB induced after 70 kVp X-ray irradiation can be classified as a complex form coupled with BD. The present simulation enables us to quantify the type of complex damage which cannot be measured through in vitro experiments and helps us to interpret the experimental detection efficiency for complex BD measured by AFM. This simple model for estimating complex DNA damage yields contributes to the precise understanding of the DNA damage complexity induced after X-ray and electron irradiations.

Microdosimetric calculations of the direct DNA damage induced by low energy electrons using the Geant4-DNA Monte Carlo code

To calculate the yield of direct DNA damage induced by low energy electrons using Monte Carlo generated microdosimetric spectra at the nanometer scale and examine the influence of various simulation inputs. The potential of classical microdosimetry to offer a viable and simpler alternative to more elaborate mechanistic approaches for practical applications is discussed. Track-structure simulations with the Geant4-DNA low-energy extension of the Geant4 Monte Carlo toolkit were used for calculating lineal energy spectra in spherical volumes with dimensions relevant to double-strand-break (DSB) induction. The microdosimetric spectra were then used to calculate the yield of simple and clustered DSB based on literature values of the threshold energy of DNA damage. The influence of the different implementations of the dielectric function of liquid water available in Geant4-DNA (Option 2 and Option 4 constructors), as well as the effect of particle tracking cutoff energy and target size are examined. Frequency- and dose-mean lineal energies in liquid-water spheres of 2, 2.3, 2.6, and 3.4 nm diameter, as well as, number of simple and clustered DSB/Gy/cell are presented for electrons over the 100 eV to 100 keV energy range. Results are presented for both the 'default' (Option 2) and 'Ioannina' (Option 4) physics models of Geant4-DNA applying several commonly used tracking cutoff energies (10, 20, 50, 100 eV). Overall, the choice of the physics model and target diameter has a moderate effect (up to ~10%-30%) on the DSB yield whereas the effect of the tracking cutoff energy may be significant (>100%). Importantly, the yield of both simple and clustered DSB was found to vary significantly (by a factor of 2 or more) with electron energy over the examined range. The yields of electron-induced simple and clustered DSB exhibit a strong energy dependence over the 100 eV-100 keV range with implications to radiation quality issues. It is shown that a classical microdosimetry approach for the calculation of DNA damage based on lineal energy spectra in nanometer-size targets predicts comparable results to computationally intensive mechanistic approaches which use detailed atomistic DNA geometries, thus, offering a relatively simple and robust alternative for some practical applications.

Computational approach to determine the relative biological effectiveness of fast neutrons using the Geant4-DNA toolkit and a DNA atomic model from the Protein Data Bank

This study proposes an innovative approach to estimate relative biological effectiveness (RBE) of fast neutrons using the Geant4 toolkit. The Geant4-DNA version cannot track heavy ions below 0.5 MeV/nucleon. In order to explore the impact of this issue, secondary particles are simulated instead of the primary low-energy neutrons. The Evaluated Nuclear Data File library is used to determine the cross sections for the elastic and inelastic interactions of neutrons with water and to find the contribution of each secondary particle spectrum. Two strategies are investigated in order to find the best possible approach and results. The first one takes into account only light particles, protons produced from elastic scattering, and α particles from inelastic scattering. Geantino particles are shot instead of heavy ions; hence all heavy ions are considered in the simulations, though their physical effects on DNA not. The second strategy takes into account all the heavy and light ions, although heavy ions cannot be tracked down to very low energies (E<0.5 MeV/nucleon). Our model is based on the combination of an atomic resolution DNA geometrical model and a Monte Carlo simulation toolkit for tracking particles. The atomic coordinates of the DNA double helix are extracted from the Protein Data Bank. Since secondary particle spectra are used instead of simulating the interaction of neutrons explicitly, this method reduces the computation times dramatically. Double-strand break induction is used as the end point for the estimation of the RBE of fast neutrons. ^{60}Co γ rays are used as the reference radiation quality. Both strategies succeed in reproducing the behavior of the RBE_{max} as a function of the incident neutron energy ranging from 0.1 to 14 MeV, including the position of its peak. A comparison of the behavior of the two strategies shows that for neutrons with energies less than 0.7 MeV, the effect of heavy ions would not be very significant, but above 0.7 MeV, heavy ions have an important role in neutron RBE.

Geometrical structures for radiation biology research as implemented in the TOPAS-nBio toolkit

Computational simulations, such as Monte Carlo track structure simulations, offer a powerful tool for quantitatively investigating radiation interactions within cells. The modelling of the spatial distribution of energy deposition events as well as diffusion of chemical free radical species, within realistic biological geometries, can help provide a comprehensive understanding of the effects of radiation on cells. Track structure simulations, however, generally require advanced computing skills to implement. The TOPAS-nBio toolkit, an extension to TOPAS (TOol for PArticle Simulation), aims to provide users with a comprehensive framework for radiobiology simulations, without the need for advanced computing skills. This includes providing users with an extensive library of advanced, realistic, biological geometries ranging from the micrometer scale (e.g. cells and organelles) down to the nanometer scale (e.g. DNA molecules and proteins). Here we present the geometries available in TOPAS-nBio.

Estimation of the radiation-induced DNA double-strand breaks number by considering cell cycle and absorbed dose per cell nucleus

DNA double-strand breaks (DSBs) are thought to be the main cause of cell death after irradiation. In this study, we estimated the probability distribution of the number of DSBs per cell nucleus by considering the DNA amount in a cell nucleus (which depends on the cell cycle) and the statistical variation in the energy imparted to the cell nucleus by X-ray irradiation. The probability estimation of DSB induction was made following these procedures: (i) making use of the Chinese Hamster Ovary (CHO)-K1 cell line as the target example, the amounts of DNA per nucleus in the logarithmic and the plateau phases of the growth curve were measured by flow cytometry with propidium iodide (PI) dyeing; (ii) the probability distribution of the DSB number per cell nucleus for each phase after irradiation with 1.0 Gy of 200 kVp X-rays was measured by means of γ-H2AX immunofluorescent staining; (iii) the distribution of the cell-specific energy deposition via secondary electrons produced by the incident X-rays was calculated by WLTrack (in-house Monte Carlo code); (iv) according to a mathematical model for estimating the DSB number per nucleus, we deduced the induction probability density of DSBs based on the measured DNA amount (depending on the cell cycle) and the calculated dose per nucleus. The model exhibited DSB induction probabilities in good agreement with the experimental results for the two phases, suggesting that the DNA amount (depending on the cell cycle) and the statistical variation in the local energy deposition are essential for estimating the DSB induction probability after X-ray exposure.

Investigation on the correlation between energy deposition and clustered DNA damage induced by low-energy electrons

This study presents the correlation between energy deposition and clustered DNA damage, based on a Monte Carlo simulation of the spectrum of direct DNA damage induced by low-energy electrons including the dissociative electron attachment. Clustered DNA damage is classified as simple and complex in terms of the combination of single-strand breaks (SSBs) or double-strand breaks (DSBs) and adjacent base damage (BD). The results show that the energy depositions associated with about 90% of total clustered DNA damage are below 150 eV. The simple clustered DNA damage, which is constituted of the combination of SSBs and adjacent BD, is dominant, accounting for 90% of all clustered DNA damage, and the spectra of the energy depositions correlating with them are similar for different primary energies. One type of simple clustered DNA damage is the combination of a SSB and 1-5 BD, which is denoted as SSB + BD. The average contribution of SSB + BD to total simple clustered DNA damage reaches up to about 84% for the considered primary energies. In all forms of SSB + BD, the SSB + BD including only one base damage is dominant (above 80%). In addition, for the considered primary energies, there is no obvious difference between the average energy depositions for a fixed complexity of SSB + BD determined by the number of base damage, but average energy depositions increase with the complexity of SSB + BD. In the complex clustered DNA damage constituted by the combination of DSBs and BD around them, a relatively simple type is a DSB combining adjacent BD, marked as DSB + BD, and it is of substantial contribution (on average up to about 82%). The spectrum of DSB + BD is given mainly by the DSB in combination with different numbers of base damage, from 1 to 5. For the considered primary energies, the DSB combined with only one base damage contributes about 83% of total DSB + BD, and the average energy deposition is about 106 eV. However, the energy deposition increases with the complexity of clustered DNA damage, and therefore, the clustered DNA damage with high complexity still needs to be considered in the study of radiation biological effects, in spite of their small contributions to all clustered DNA damage.

Measurements and simulations of microscopic damage to DNA in water by 30 keV electrons: A general approach applicable to other radiation sources and biological targets

The determination of the microscopic dose-damage relationship for DNA in an aqueous environment is of a fundamental interest for dosimetry and applications in radiation therapy and protection. We combine geant4 particle-scattering simulations in water with calculations concerning the movement of biomolecules to obtain the energy deposit in the biologically relevant nanoscopic volume. We juxtaposition these results to the experimentally determined damage to obtain the dose-damage relationship at a molecular level. This approach is tested for an experimentally challenging system concerning the direct irradiation of plasmid DNA (pUC19) in water with electrons as primary particles. Here a microscopic target model for the plasmid DNA based on the relation of lineal energy and radiation quality is used to calculate the effective target volume. It was found that on average fewer than two ionizations within a 7.5-nm radius around the sugar-phosphate backbone are sufficient to cause a single strand break, with a corresponding median lethal energy deposit being E_{1/2}=6±4 eV. The presented method is applicable for ionizing radiation (e.g., γ rays, x rays, and electrons) and a variety of targets, such as DNA, proteins, or cells.

Track structure model of microscopic energy deposition by protons and heavy ions in segments of neuronal cell dendrites represented by cylinders or spheres

Changes to cognition, including memory, following radiation exposure are a concern for cosmic ray exposures to astronauts and in Hadron therapy with proton and heavy ion beams. The purpose of the present work is to develop computational methods to evaluate microscopic energy deposition (ED) in volumes representative of neuron cell structures, including segments of dendrites and spines, using a stochastic track structure model. A challenge for biophysical models of neuronal damage is the large sizes (> 100µm) and variability in volumes of possible dendritic segments and pre-synaptic elements (spines and filopodia). We consider cylindrical and spherical microscopic volumes of varying geometric parameters and aspect ratios from 0.5 to 5 irradiated by protons, and ³He and ¹²C particles at energies corresponding to a distance of 1cm to the Bragg peak, which represent particles of interest in Hadron therapy as well as space radiation exposure. We investigate the optimal axis length of dendritic segments to evaluate microscopic ED and hit probabilities along the dendritic branches at a given macroscopic dose. Because of large computation times to analyze ED in volumes of varying sizes, we developed an analytical method to find the mean primary dose in spheres that can guide numerical methods to find the primary dose distribution for cylinders. Considering cylindrical segments of varying aspect ratio at constant volume, we assess the chord length distribution, mean number of hits and ED profiles by primary particles and secondary electrons (δ-rays). For biophysical modeling applications, segments on dendritic branches are proposed to have equal diameters and axes lengths along the varying diameter of a dendritic branch.

Comments on potential health effects of MRI-induced DNA lesions: quality is more important to consider than quantity

Magnetic resonance imaging (MRI) is increasingly being used in cardiology to detect heart disease and guide therapy. It is mooted to be a safer alternative to imaging techniques, such as computed tomography (CT) or coronary angiographic imaging. However, there has recently been an increased interest in the potential long-term health risks of MRI, especially in the light of the controversy resulting from a small number of research studies reporting an increase in DNA damage following exposure, with calls to limit its use and avoid unnecessary examination, according to the precautionary principle. Overall the published data are somewhat limited and inconsistent; the ability of MRI to produce DNA lesions has yet to be robustly demonstrated and future experiments should be carefully designed to optimize sensitivity and benchmarked to validate and assess reproducibility. The majority of the current studies have focussed on the initial induction of DNA damage, and this has led to comparisons between the reported induction of γH2AX and implied double-strand break (DSB) yields produced following MRI with induction by imaging techniques using ionizing radiation. However, γH2AX is not only a marker of classical double-ended DSB, but also a marker of stalled replication forks and in certain circumstances stalled DNA transcription. Additionally, ionizing radiation is efficient at producing complex DNA damage, unique to ionizing radiation, with an associated reduction in repairability. Even if the fields associated with MRI are capable of producing DNA damage, the lesions produced will in general be simple, similar to those produced by endogenous processes. It is therefore inappropriate to try and infer cancer risk by simply comparing the yields of γH2AX foci or DNA lesions potentially produced by MRI to those produced by a given exposure of ionizing radiation, which will generally be more biologically effective and have a greater probability of leading to long-term health effects. As a result, it is important to concentrate on more relevant downstream end points (e.g. chromosome aberration production), along with potential mechanisms by which MRI may lead to DNA lesions. This could potentially involve a perturbation in homeostasis of oxidative stress, modifying the background rate of endogenous DNA damage induction. In summary, what the field needs at the moment is more research and less fear mongering.

Monte Carlo calculations of energy deposition distributions of electrons below 20 keV in protein

The distributions of energy depositions of electrons in semi-infinite bulk protein and the radial dose distributions of point-isotropic mono-energetic electron sources [i.e., the so-called dose point kernel (DPK)] in protein have been systematically calculated in the energy range below 20 keV, based on Monte Carlo methods. The ranges of electrons have been evaluated by extrapolating two calculated distributions, respectively, and the evaluated ranges of electrons are compared with the electron mean path length in protein which has been calculated by using electron inelastic cross sections described in this work in the continuous-slowing-down approximation. It has been found that for a given energy, the electron mean path length is smaller than the electron range evaluated from DPK, but it is large compared to the electron range obtained from the energy deposition distributions of electrons in semi-infinite bulk protein. The energy dependences of the extrapolated electron ranges based on the two investigated distributions are given, respectively, in a power-law form. In addition, the DPK in protein has also been compared with that in liquid water. An evident difference between the two DPKs is observed. The calculations presented in this work may be useful in studies of radiation effects on proteins.

Health risks of space exploration: targeted and nontargeted oxidative injury by high-charge and high-energy particles

SIGNIFICANCE

During deep space travel, astronauts are often exposed to high atomic number (Z) and high-energy (E) (high charge and high energy [HZE]) particles. On interaction with cells, these particles cause severe oxidative injury and result in unique biological responses. When cell populations are exposed to low fluences of HZE particles, a significant fraction of the cells are not traversed by a primary radiation track, and yet, oxidative stress induced in the targeted cells may spread to nearby bystander cells. The long-term effects are more complex because the oxidative effects persist in progeny of the targeted and affected bystander cells, which promote genomic instability and may increase the risk of age-related cancer and degenerative diseases.

RECENT ADVANCES

Greater understanding of the spatial and temporal features of reactive oxygen species bursts along the tracks of HZE particles, and the availability of facilities that can simulate exposure to space radiations have supported the characterization of oxidative stress from targeted and nontargeted effects.

CRITICAL ISSUES

The significance of secondary radiations generated from the interaction of the primary HZE particles with biological material and the mitigating effects of antioxidants on various cellular injuries are central to understanding nontargeted effects and alleviating tissue injury.

FUTURE DIRECTIONS

Elucidation of the mechanisms underlying the cellular responses to HZE particles, particularly under reduced gravity and situations of exposure to additional radiations, such as protons, should be useful in reducing the uncertainty associated with current models for predicting long-term health risks of space radiation. These studies are also relevant to hadron therapy of cancer.

Radioactivity and lung cancer-mathematical models of radionuclide deposition in the human lungs

The human respiratory tract is regarded as pathway for radionuclides and other hazardous airborne materials to enter the body. Radioactive particles inhaled and deposited in the lungs cause an irradiation of bronchial/alveolar tissues. At the worst, this results in a malignant cellular transformation and, as a consequence of that, the development of lung cancer. In general, naturally occurring radionuclides (e.g., (222)Rn, (40)K) are attached to so-called carrier aerosols. The aerodynamic diameters of such radioactively labeled particles generally vary between several nanometers (ultrafine particles) and few micrometers, whereby highest particle fractions adopt sizes around 100 nm. Theoretical simulations of radioactive particle deposition in the human lungs were based on a stochastic lung geometry and a particle transport/deposition model using the random-walk algorithm. Further a polydisperse carrier aerosol (diameter: 1 nm-10 µm, ρ ≈ 1 g cm(-3)) with irregularly shaped particles and the effect of breathing characteristics and certain respiratory parameters on the transport of radioactive particles to bronchial/alveolar tissues were considered. As clearly shown by the results of deposition modeling, distribution patterns of radiation doses mainly depend on the size of the carrier aerosol. Ultrafine (< 10 nm) and large (> 2 µm) aerosol particles are preferentially deposited in the extrathoracic and upper bronchial region, whereas aerosol particles with intermediate size (10 nm-2 µm) may penetrate to deeper lung regions, causing an enhanced damage of the alveolar tissue by the attached radionuclides.

Microdosimetry of low-energy electrons

PURPOSE

To investigate differences in energy depositions and microdosimetric parameters of low-energy electrons in liquid and gaseous water using Monte Carlo track structure simulations.

MATERIALS AND METHODS

KURBUC-liq (Kyushu University and Radiobiology Unit Code for liquid water) was used for simulating electron tracks in liquid water. The inelastic scattering cross sections of liquid water were obtained from the dielectric response model of Emfietzoglou et al. (Radiation Research 2005;164:202-211). Frequencies of energy deposited in nanometre-size cylindrical targets per unit absorbed dose and associated lineal energies were calculated for 100-5000 eV monoenergetic electrons and the electron spectrum of carbon K edge X-rays. The results for liquid water were compared with those for water vapour.

RESULTS

Regardless of electron energy, there is a limit how much energy electron tracks can deposit in a target. Phase effects on the frequencies of energy depositions are largely visible for the targets with diameters and heights smaller than 30 nm. For the target of 2.3 nm by 2.3 nm (similar to dimension of DNA segments), the calculated frequency- and dose-mean lineal energies for liquid water are up to 40% smaller than those for water vapour. The corresponding difference is less than 12% for the targets with diameters ≥ 30 nm.

CONCLUSIONS

Condensed-phase effects are non-negligible for microdosimetry of low-energy electrons for targets with sizes smaller than a few tens of nanometres, similar to dimensions of DNA molecular structures and nucleosomes.

Formulation of the multi-hit model with a non-Poisson distribution of hits

PURPOSE

We proposed a formulation of the multi-hit single-target model in which the Poisson distribution of hits was replaced by a combination of two distributions: one for the number of particles entering the target and one for the number of hits a particle entering the target produces. Such an approach reflects the fact that radiation damage is a result of two different random processes: particle emission by a radiation source and interaction of particles with matter inside the target.

METHODS AND MATERIALS

Poisson distribution is well justified for the first of the two processes. The second distribution depends on how a hit is defined. To test our approach, we assumed that the second distribution was also a Poisson distribution. The two distributions combined resulted in a non-Poisson distribution. We tested the proposed model by comparing it with previously reported data for DNA single- and double-strand breaks induced by protons and electrons, for survival of a range of cell lines, and variation of the initial slopes of survival curves with radiation quality for heavy-ion beams.

RESULTS

Analysis of cell survival equations for this new model showed that they had realistic properties overall, such as the initial and high-dose slopes of survival curves, the shoulder, and relative biological effectiveness (RBE) In most cases tested, a better fit of survival curves was achieved with the new model than with the linear-quadratic model. The results also suggested that the proposed approach may extend the multi-hit model beyond its traditional role in analysis of survival curves to predicting effects of radiation quality and analysis of DNA strand breaks.

CONCLUSIONS

Our model, although conceptually simple, performed well in all tests. The model was able to consistently fit data for both cell survival and DNA single- and double-strand breaks. It correctly predicted the dependence of radiation effects on parameters of radiation quality.

Ionizing radiation-induced metabolic oxidative stress and prolonged cell injury

Cellular exposure to ionizing radiation leads to oxidizing events that alter atomic structure through direct interactions of radiation with target macromolecules or via products of water radiolysis. Further, the oxidative damage may spread from the targeted to neighboring, non-targeted bystander cells through redox-modulated intercellular communication mechanisms. To cope with the induced stress and the changes in the redox environment, organisms elicit transient responses at the molecular, cellular and tissue levels to counteract toxic effects of radiation. Metabolic pathways are induced during and shortly after the exposure. Depending on radiation dose, dose-rate and quality, these protective mechanisms may or may not be sufficient to cope with the stress. When the harmful effects exceed those of homeostatic biochemical processes, induced biological changes persist and may be propagated to progeny cells. Physiological levels of reactive oxygen and nitrogen species play critical roles in many cellular functions. In irradiated cells, levels of these reactive species may be increased due to perturbations in oxidative metabolism and chronic inflammatory responses, thereby contributing to the long-term effects of exposure to ionizing radiation on genomic stability. Here, in addition to immediate biological effects of water radiolysis on DNA damage, we also discuss the role of mitochondria in the delayed outcomes of ionization radiation. Defects in mitochondrial functions lead to accelerated aging and numerous pathological conditions. Different types of radiation vary in their linear energy transfer (LET) properties, and we discuss their effects on various aspects of mitochondrial physiology. These include short and long-term in vitro and in vivo effects on mitochondrial DNA, mitochondrial protein import and metabolic and antioxidant enzymes.

Applications of the microdosimetric function implemented in the macroscopic particle transport simulation code PHITS

PURPOSE

Microdosimetric quantities such as lineal energy are generally considered to be better indices than linear energy transfer (LET) for expressing the relative biological effectiveness (RBE) of high charge and energy particles. To calculate their probability densities (PD) in macroscopic matter, it is necessary to integrate microdosimetric tools such as track-structure simulation codes with macroscopic particle transport simulation codes.

METHODS

As an integration approach, the mathematical model for calculating the PD of microdosimetric quantities developed based on track-structure simulations was incorporated into the macroscopic particle transport simulation code PHITS (Particle and Heavy Ion Transport code System). The improved PHITS enables the PD in macroscopic matter to be calculated within a reasonable computation time, while taking their stochastic nature into account.

APPLICATIONS

The microdosimetric function of PHITS was applied to biological dose estimation for charged-particle therapy and risk estimation for astronauts. The former application was performed in combination with the microdosimetric kinetic model, while the latter employed the radiation quality factor expressed as a function of lineal energy.

CONCLUSION

Owing to the unique features of the microdosimetric function, the improved PHITS has the potential to establish more sophisticated systems for radiological protection in space as well as for the treatment planning of charged-particle therapy.

Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC

This review describes the PARTRAC suite of comprehensive Monte Carlo simulation tools for calculations of track structures of a variety of ionizing radiation qualities and their biological effects. A multi-scale target model characterizes essential structures of the whole genomic DNA within human fibroblasts and lymphocytes in atomic resolution. Calculation methods and essential results are recapitulated regarding the physical, physico-chemical and chemical stage of track structure development of radiation damage induction. Recent model extension towards DNA repair processes extends the time dimension by about 12 orders of magnitude and paves the way for superior predictions of radiation risks.

The role of hydration in the distribution of free radical trapping in directly ionized DNA

The purpose of this study was to elucidate the role of hydration (Gamma) in the distribution of free radical trapping in directly ionized DNA. Solid-state films of pUC18 (2686 bp) plasmids were hydrated to Gamma in the range 2.5 < or = Gamma < or = 22.5 mol water/mol nucleotide. Free radical yields, G(Sigmafr), measured by EPR at 4 K are seen to increase from 0.28 +/- 0.01 micromol/J at Gamma = 2.5 to 0.63 +/- 0.01 micromol/J at Gamma= 22.5, respectively. Based on a semi-empirical model of the free radical trapping events that follow the initial ionizations of the DNA components, we conclude that two-thirds of the holes formed on the inner solvation shell (Gamma < 10) transfer to the sugar-phosphate backbone. Likewise, of the holes produced by direct ionization of the sugar-phosphate, about one-third are trapped by deprotonation as neutral sugar-phosphate radical species, while the remaining two-thirds are found to transfer to the bases. This analysis provides the best measure to date for the probability of hole transfer (approximately 67%) into the base stack. It can thus be predicted that the distribution of holes formed in fully hydrated DNA at 4 K will be 78% on the bases and 22% on the sugar-phosphate. Adding the radicals due to electron attachment (confined to the pyrimidine bases), the distribution of all trapped radicals will be 89% on the bases and 11% on the sugar-phosphate backbone. This prediction is supported by partitioning results obtained from the high dose-response curves fitted to the two-component model. These results not only add to our understanding of how the holes redistribute after ionization but are also central to predicting the yield and location of strand breaks in DNA exposed to the direct effects of ionizing radiation.

A contribution to the linear no-threshold discussion

The paper approaches the linear no-threshold (LNT) hypothesis, currently used as the basis for recommendations in radiological protection, from the point of view of the radiation mechanism. All considerations of the validity of the LNT hypothesis based on experiment or epidemiology are dismissed because of the impossibility of deriving statistically significant data at very low doses. Instead, the LNT hypothesis is assessed from a consideration of the mechanism of radiation action. The DNA double-strand break is proposed to be the crucial radiation-induced molecular lesion. A trace is made using a series of correlations that link the DNA double-strand break to effects at the cellular level and these cellular effects are linked to the induction of cancer. Multistep modelling of carcinogenesis is used to take the link through to a consideration of radiation risk. It is concluded that, from the point of view of radiation mechanism, at very low doses the LNT hypothesis of radiation action is valid, that is, the risk function has a positive slope from zero dose.

Cross sections for low-energy (10-50 eV) electron damage to DNA

We report direct measurements of the formation of single-, double- and multiple strand breaks in pure plasmid DNA as a function of exposure to 10-50 eV electrons. The effective cross sections to produce these different types of DNA strand breaks were determined and were found to range from approximately 10(-17) to 3 x 10(-15) cm(2). The total effective cross section and the effective range for destruction of supercoiled DNA extend from 3.4 to 4.4 x 10(-15) cm(2) and 12 to 14 nm, respectively, over the range 10-50 eV. The variation of the effective cross sections with electron energy is discussed in terms of the electron's inelastic mean free path, penetration depth, and dissociation mechanisms, including resonant electron capture; the latter is found to dominate the effective cross sections for single- and double-strand breaks at 10 eV. The most striking observations are that (1) supercoiled DNA is approximately one order of magnitude more sensitive to the formation of double-strand breaks by low-energy electrons than is relaxed circular DNA, and (2) the dependence of the effective cross sections on the incident electron energy is unrelated to the corresponding ionization cross sections. This finding suggests that the traditional notion that radiobiological damage is related to the number of ionization events would not apply at very low energies.

Model for radial dependence of frequency distributions for energy imparted in nanometer volumes from HZE particles

This paper develops a deterministic model of frequency distributions for energy imparted (total energy deposition) in small volumes similar to DNA molecules from high-energy ions of interest for space radiation protection and cancer therapy. Frequency distributions for energy imparted are useful for considering radiation quality and for modeling biological damage produced by ionizing radiation. For high-energy ions, secondary electron (delta-ray) tracks originating from a primary ion track make dominant contributions to energy deposition events in small volumes. Our method uses the distribution of electrons produced about an ion's path and incorporates results from Monte Carlo simulation of electron tracks to predict frequency distributions for ions, including their dependence on radial distance. The contribution from primary ion events is treated using an impact parameter formalism of spatially restricted linear energy transfer (LET) and energy-transfer straggling. We validate our model by comparing it directly to results from Monte Carlo simulations for proton and alpha-particle tracks. We show for the first time frequency distributions of energy imparted in DNA structures by several high-energy ions such as cosmic-ray iron ions. Our comparison with results from Monte Carlo simulations at low energies indicates the accuracy of the method.

Issues in risk assessment from solar particle events

Uncertainties in risk assessment from solar particle events (SPE) include the role of high linear energy transfer (LET) secondary ions, the assessment of dose-rate effects as they relate to acute injury, the risk of cancer mortality and the modification of health effects due to the stress of spaceflight. We discuss several issues where new knowledge is required for improving estimates of radiation risk from SPE's. Secondary particles such as neutrons and low energy and charge ions (LZE) may dominate radiation risk behind a storm shelter and their biological effects are poorly understood, especially at low dose-rate. Dose-rate modulation of radiation response is also related to genetic pre-disposition an important determinant of radiation sensitivity. Molecular pathways that control cell death and tissue response have been elucidated in recent years and should provide new understanding of dose-rate effects for risk assessment. We consider some of these factors and discuss calculations using radiation transport codes, track structure models of energy deposition, and a molecular kinetics approach to model radiation response.

A novel algorithm for tracing the interaction of a track with molecular targets--use of Delaunay triangulation

This paper describes an efficient method for tracing interactions between a radiation track and molecular targets in a cell or nucleus. The method is efficient because it ensures that a minimal subset of interactions in the track needs to be tested for each score. It is most useful for high-energy particles since the number of interactions and range of the tracks increase with particle energy. In this method, the set of interactions, i.e. the track, is considered as a collection of points in three-dimensional space. This set, together with the eight vertices of a bounding box, is discretised as a collection of tetrahedra, each of which satisfies the Delaunay criterion. Because of the geometric properties of these tetrahedra, only those points which are connected to vertices of tetrahedra through which the target passes need be tested for scoring. An efficient algorithm is used to follow the track through the cell or the nucleus. The main benefit of the method is that it eliminates the need to test those interactions which are geometrically distant from the target. In particular, in cases where there are no interactions in the target, very few interactions need to be checked.

Computational modelling of low-energy electron-induced DNA damage by early physical and chemical events

Modelling and calculations are presented as a first step towards mechanistic interpretation and prediction of radiation effects based on the spectrum of initial DNA damage produced by low energy electrons (100 eV-4.5 keV) that can be compared with experimental information. Relative yields of single and clustered strand breaks are presented in terms of complexity and source of damage, either by direct energy deposition or by reaction of OH radicals, and dependence on the activation probability of OH radicals and the amount of energy required to give a single strand break (ssb). Data show that the majority of interactions in DNA do not lead to damage in the form of strand breaks and when they do occur, they are most frequently simple ssb. However, for double-strand breaks (dsb), a high proportion (approximately 30%) are of more complex forms, even without considering additional complexity from base damage. The greater contribution is from direct interactions in the DNA but reactions of OH radicals add substantially to this, both in terms of the total number of breaks and in increasing the complexity within a cluster. It has been shown that the lengths of damaged segments of DNA from individual electron tracks tend to be short, indicating that consequent deletion length (simply by loss of a fragment between nearby dsb) would be short, very seldom exceeding a few tens of base pairs.

Modelling of radiation-induced DNA damage: the early physical and chemical event

A Monte Carlo track structure calculation of single- and double-strand breaks induced by direct energy deposition in DNA and by interacting diffusible .OH radicals with DNA has been made for low energy electrons. The .OH radicals generated within 4 nm of linear segments of DNA were diffused in order to mimic the mean diffusion distance in the cellular environment. The reactions of the radical species .OH, .H and e-aq were included in this study. The calculated values for the yield of single- and double-strand breaks have been compared with experimentally determined values from the literature. The calculations indicate, too, that the majority of dsb have additional associated damage, constituting clustered lesions of varying complexity.

DNA structure, hydration and dynamics

Although the double helical model of DNA structure is now 40 years old, there is still considerable effort being made to elucidate the range of conformations that can be adopted by this flexible molecule. We review the current state of our knowledge of DNA structure which is available from both experimental and computational approaches.

Monte Carlo track structure studies of energy deposition and calculation of initial DSB and RBE

Estimation of exposure due to environmental and other sources of radiations of high-LET and low-LET is of interest in radiobiology and radiation protection for risk assessment. To account for the differences in effectiveness of different types of radiations various parameters have been used. However, the relative inadequacy of the commonly used parameters, including dose, fluence, linear energy transfer, lineal energy, specific energy and quality factor, has been made manifest by the biological importance of the microscopic track structure and primary modes of interaction. Monte Carlo track structure simulations have been used to calculate the frequency of energy deposition by radiations of high- and low-LET in target sizes similar to DNA and higher order genomic structure. Tracks of monoenergetic heavy ions and electrons were constructed by following the molecular interaction-by-interaction histories of the particles down to 10 eV. Subsequently, geometrical models of these assumed biological targets were randomly exposed to the radiation tracks and the frequency of energy depositions obtained were normalized to unit dose in unit density liquid water (l0(3) kg m-3). From these data and a more sophisticated model of the DNA, absolute yields of both single- and double-strand breaks expressed in number of breaks per dalton per Gray were obtained and compared with the measured yields. The relative biological effectiveness (RBE) for energy depositions in cylindrical targets has been calculated using 100 keV electrons as the reference radiation assuming the electron track-ends contribution is similar to that in 250 kV X-ray or Co60 gamma-ray irradiations.

Target model of nucleosome particle for track structure calculations and DNA damage modelling

A non-homogeneous target model of nucleosome particle for track structure calculations and biophysical modelling of radiation-induced DNA damage was developed together with a sampling method for non-homogeneous targets. The model allows one to distinguish energy deposited in the histone core, primary hydration shell and nucleosome DNA. Illustrative results for low energy electrons are presented and compared with results obtained with simpler models of a nucleosome particle.