Challenges in Monte Carlo track structure modelling

Absolute cross-sections for DNA strand breaks and crosslinks induced by low energy electrons

Absolute cross sections (CSs) for the interaction of low energy electrons with condensed macromolecules are essential parameters to accurately model ionizing radiation induced reactions. To determine CSs for various conformational DNA damage induced by 2-20 eV electrons, we investigated the influence of the attenuation length (AL) and penetration factor (f) using a mathematical model. Solid films of supercoiled plasmid DNA with thicknesses of 10, 15 and 20 nm were irradiated with 4.6, 5.6, 9.6 and 14.6 eV electrons. DNA conformational changes were quantified by gel electrophoresis, and the respective yields were extrapolated from exposure-response curves. The absolute CS, AL and f values were generated by applying the model developed by Rezaee et al. The values of AL were found to lie between 11 and 16 nm with the maximum at 14.6 eV. The absolute CSs for the loss of the supercoiled (LS) configuration and production of crosslinks (CL), single strand breaks (SSB) and double strand breaks (DSB) induced by 4.6, 5.6, 9.6 and 14.6 eV electrons are obtained. The CSs for SSB are smaller, but similar to those for LS, indicating that SSB are the main conformational damage. The CSs for DSB and CL are about one order of magnitude smaller than those of LS and SSB. The value of f is found to be independent of electron energy, which allows extending the absolute CSs for these types of damage within the range 2-20 eV, from previous measurements of effective CSs. When comparison is possible, the absolute CSs are found to be in good agreement with those obtained from previous similar studies with double-stranded DNA. The high values of the absolute CSs of 4.6 and 9.6 eV provide quantitative evidence for the high efficiency of low energy electrons to induce DNA damage via the formation of transient anions.

Monte Carlo Electron Track Structure Calculations in Liquid Water Using a New Model Dielectric Response Function

Monte Carlo track structure codes provide valuable information for understanding radiation effects down to the DNA level, where experimental measurements are most difficult or unavailable. It is well recognized that the performance of such codes, especially at low energies and/or subcellular level, critically depends on the reliability of the interaction cross sections that are used as input in the simulation. For biological media such as liquid water, one of the most challenging issues is the role of condensed-phase effects. For inelastic scattering, such effects can be conveniently accounted for through the complex dielectric response function of the media. However, for this function to be useful it must fulfill some important sum rules and have a simple analytic form for arbitrary energy- and momentum-transfer. The Emfietzoglou-Cucinotta-Nikjoo (ECN) model offers a practical, self-consistent and fully analytic parameterization of the dielectric function of liquid water based on the best available experimental data. An important feature of the ECN model is that it includes, in a phenomenological manner, exchange and correlation effects among the screening electrons, thus, going beyond the random-phase approximation implicit in earlier models. In this work, inelastic cross sections beyond the plane wave Born approximation are calculated for low-energy electrons (10 eV-10 keV) based on the ECN model, and used for Monte Carlo track structure simulations of physical quantities relevant to the microdosimetry of low-energy electrons in liquid water. Important new developments in the physics of inelastic scattering are discussed and their effect on electron track structure is investigated by a comparison against simulations (under otherwise identical conditions) using the Born approximation and a simpler form of the dielectric function based on the Oak Ridge National Laboratory model. The results reveal that both the dielectric function and the corrections to the Born approximation may have a sizeable effect on track structure calculations at the nanometer scale (DNA level), where the details of inelastic scattering and the role of low-energy electrons are most critical.

Proton and light ion RBE for the induction of direct DNA double strand breaks

PURPOSE

To present and characterize a Monte Carlo (MC) tool for the simulation of the relative biological effectiveness for the induction of direct DNA double strand breaks (RBEDSB (direct)) for protons and light ions.

METHODS

The MC tool uses a pregenerated event-by-event tracks library of protons and light ions that are overlaid on a cell nucleus model. The cell nucleus model is a cylindrical arrangement of nucleosome structures consisting of 198 DNA base pairs. An algorithm relying on k-dimensional trees and cylindrical symmetries is used to search coincidences of energy deposition sites with volumes corresponding to the sugar-phosphate backbone of the DNA molecule. Strand breaks (SBs) are scored when energy higher than a threshold is reached in these volumes. Based on the number of affected strands, they are categorized into either single strand break (SSB) or double strand break (DSB) lesions. The number of SBs composing each lesion (i.e., its size) is also recorded. RBEDSB (direct) is obtained by taking the ratio of DSB yields of a given radiation field to a (60)Co field. The MC tool was used to obtain SSB yields, DSB yields, and RBEDSB (direct) as a function of linear energy transfer (LET) for protons ((1)H(+)), (4)He(2+), (7)Li(3+), and (12)C(6+) ions.

RESULTS

For protons, the SSB yields decreased and the DSB yields increased with LET. At ≈24.5 keV μm(-1), protons generated 15% more DSBs than (12)C(6+) ions. The RBEDSB (direct) varied between 1.24 and 1.77 for proton fields between 8.5 and 30.2 keV μm(-1), and it was higher for iso-LET ions with lowest atomic number. The SSB and DSB lesion sizes showed significant differences for all radiation fields. Generally, the yields of SSB lesions of sizes ≥2 and the yields of DSB lesions of sizes ≥3 increased with LET and increased for iso-LET ions of lower atomic number. On the other hand, the ratios of SSB to DSB lesions of sizes 2-4 did not show variability with LET nor projectile atomic number, suggesting that these metrics are independent of the radiation quality. Finally, a variance of up to 8% in the DSB yields was observed as a function of the particle incidence angle on the cell nucleus. This simulation effect is due to the preferential alignment of ion tracks with the DNA nucleosomes at specific angles.

CONCLUSIONS

The MC tool can predict SSB and DSB yields for light ions of various LET and estimate RBEDSB (direct). In addition, it can calculate the frequencies of different DNA lesion sizes, which is of interest in the context of biologically relevant absolute dosimetry of particle beams.

Nanoscale radiation transport and clinical beam modeling for gold nanoparticle dose enhanced radiotherapy (GNPT) using X-rays

We review radiation transport and clinical beam modelling for gold nanoparticle dose-enhanced radiotherapy using X-rays. We focus on the nanoscale radiation transport and its relation to macroscopic dosimetry for monoenergetic and clinical beams. Among other aspects, we discuss Monte Carlo and deterministic methods and their applications to predicting dose enhancement using various metrics.

Lauriston S. Taylor Lecture on radiation protection and measurements: what makes particle radiation so effective?

The scientific basis for the physical and biological effectiveness of particle radiations has emerged from many decades of meticulous basic research. A diverse array of biologically relevant consequences at the molecular, cellular, tissue, and organism level have been reported, but what are the key processes and mechanisms that make particle radiation so effective, and what competing processes define dose dependences? Recent studies have shown that individual genotypes control radiation-regulated genes and pathways in response to radiations of varying ionization density. The fact that densely ionizing radiations can affect different gene families than sparsely ionizing radiations, and that the effects are dose- and time-dependent, has opened up new areas of future research. The complex microenvironment of the stroma and the significant contributions of the immune response have added to our understanding of tissue-specific differences across the linear energy transfer (LET) spectrum. The importance of targeted versus nontargeted effects remains a thorny but elusive and important contributor to chronic low dose radiation effects of variable LET that still needs further research. The induction of cancer is also LET-dependent, suggesting different mechanisms of action across the gradient of ionization density. The focus of this 35th Lauriston S. Taylor Lecture is to chronicle the step-by-step acquisition of experimental clues that have refined our understanding of what makes particle radiation so effective, with emphasis on the example of radiation effects on the crystalline lens of the human eye.

Absolute cross section for low-energy-electron damage to condensed macromolecules: a case study of DNA

Cross sections (CSs) for the interaction of low-energy electrons (LEE) with condensed macromolecules are essential parameters for accurate modeling of radiation-induced molecular decomposition and chemical synthesis. Electron irradiation of dry nanometer-scale macromolecular solid films has often been employed to measure CSs and other quantitative parameters for LEE interactions. Since such films have thicknesses comparable with electron thermalization distances, energy deposition varies throughout the film. Moreover, charge accumulation occurring inside the films shields a proportion of the macromolecules from electron irradiation. Such effects complicate the quantitative comparison of the CSs obtained in films of different thicknesses and limit the applicability of such measurements. Here, we develop a simple mathematical model, termed the molecular survival model, that employs a CS for a particular damage process together with an attenuation length related to the total CS, to investigate how a measured CS might be expected to vary with experimental conditions. As a case study, we measure the absolute CS for the formation of DNA strand breaks (SBs) by electron irradiation at 10 and 100 eV of lyophilized plasmid DNA films with thicknesses between 10 and 30 nm. The measurements are shown to depend strongly on the thickness and charging condition of the nanometer-scale films. Such behaviors are in accord with the model and support its validity. Via this analysis, the CS obtained for SB damage is nearly independent of film thickness and charging effects. In principle, this model can be adapted to provide absolute CSs for electron-induced damage or reactions occurring in other molecular solids across a wider range of experimental conditions.