Charged particle tracks in water of different phases. Monte Carlo simulation of electron tracks

Review of the Geant4-DNA Simulation Toolkit for Radiobiological Applications at the Cellular and DNA Level

Simple Summary

A brief description of the methodologies to simulate ionizing radiation transport in biologically relevant matter is presented. Emphasis is given to the physical, chemical, and biological models of Geant4-DNA that enable mechanistic radiobiological modeling at the cellular and DNA level, important to improve the efficacy of existing and novel radiotherapeutic modalities for the treatment of cancer.

Abstract

The Geant4-DNA low energy extension of the Geant4 Monte Carlo (MC) toolkit is a continuously evolving MC simulation code permitting mechanistic studies of cellular radiobiological effects. Geant4-DNA considers the physical, chemical, and biological stages of the action of ionizing radiation (in the form of x- and γ-ray photons, electrons and β^±-rays, hadrons, α-particles, and a set of heavier ions) in living cells towards a variety of applications ranging from predicting radiotherapy outcomes to radiation protection both on earth and in space. In this work, we provide a brief, yet concise, overview of the progress that has been achieved so far concerning the different physical, physicochemical, chemical, and biological models implemented into Geant4-DNA, highlighting the latest developments. Specifically, the “dnadamage1” and “molecularDNA” applications which enable, for the first time within an open-source platform, quantitative predictions of early DNA damage in terms of single-strand-breaks (SSBs), double-strand-breaks (DSBs), and more complex clustered lesions for different DNA structures ranging from the nucleotide level to the entire genome. These developments are critically presented and discussed along with key benchmarking results. The Geant4-DNA toolkit, through its different set of models and functionalities, offers unique capabilities for elucidating the problem of radiation quality or the relative biological effectiveness (RBE) of different ionizing radiations which underlines nearly the whole spectrum of radiotherapeutic modalities, from external high-energy hadron beams to internal low-energy gamma and beta emitters that are used in brachytherapy sources and radiopharmaceuticals, respectively.

A review of simulation codes and approaches for radiation chemistry

Historically, the field of radiation chemistry began shortly after the discovery of radioactivity, and its development has been closely related to discoveries in other related fields such as radiation and nuclear physics. Radiolysis of water and radiation chemistry have been very important in elucidating how radiation affects living matter and how it induces DNA damage. Nowadays, we recognize the importance of chemistry to understanding the effects of radiation on cells; however, it took several decades to obtain this insight, and much is still unknown. The radiolysis of water and aqueous solutions have been the subject of much experimental and theoretical research for many decades. One important concept closely related to radiation chemistry is radiation track structure. Track structure results from early physical and physicochemical events that lead to a highly non-homogenous distribution of radiolytic species. Because ionizing radiation creates unstable species that are distributed non-homogenously, the use of conventional reaction kinetics methods does not describe this chemistry well. In recent years, several methods have been developed for simulating radiation chemistry. In this review, we give a brief history of the field and the development of the simulation codes. We review the current methods used to simulate radiolysis of water and radiation chemistry, and we describe several radiation chemistry codes and their applications.

Lethal events in V79 cells irradiated by low-energy protons and correlations with distribution patterns of energy deposition, radical concentration and DNA damage

Published survival data of V79 cells irradiated by 0.5-5.0 MeV (7-40 keV/microm) protons have been analyzed with a detailed radiobiological model to estimate the per-track yields of lethal lesions. Their correlations with distribution patterns of deposited energy, radical concentrations and with the yields of specific classes of DNA damage have been studied. The observed correlations indicate a potential interpretation of DNA damage lethal for the cell and the initial physical and chemical processes leading to such damage.

Ionisation cluster-size formation by electrons: from macroscopic to nanometric target sizes

An indispensable prerequisite for a deeper understanding of specified physical, chemical or biological changes initiated in matter when being exposed to ionising radiation is a detailed knowledge of particle track structure. Here, the structure of electron tracks is of particular interest since electrons are set in motion in large numbers as secondary particles during the slow down of any kind of ionising radiation in matter. From the point of view of radiation induced early damage to genes and cells, which starts with the early damage to segments of the DNA molecule, the most effective secondary electrons are those at energies of a few hundred eV since the yield of double-strand breaks induced by such electrons in the DNA shows a maximum. This can be explained by the fact that in water cylinders, 2 nm in diameter and height (as a substitute to small segments of the DNA), the probability of the electron-induced formation of ionisation cluster sizes greater than or equal to two is highest also at initial electron energies of a few hundred eV. In view of this promising feature of ionisation cluster-size distributions formed by low-energy electrons in nanometric targets of liquid water for explaining particular radio-biological endpoints, it is the aim of the present work to investigate the properties of cluster-size formation by electrons as a function of target size. Here, main emphasis is laid on the behaviour of cluster-size distributions if the target size is reduced from macroscopic to nanometric volumes.

The comparison of calculated and experimental microdosimetric distributions for carbon ions

The aim of this work is to present microdosimetric characteristics of 400 MeV amu(-1) and 500 MeV amu(-1) carbon ions obtained by theoretical calculations and to analyse them with respect to experimental data obtained by tissue-equivalent proportional counter in a scope of project ICCHIBAN and by etched track detector CR 39 Page irradiated by LHE nuclotron at JINR, Dubna, Russia. Track structures provided by Monte Carlo code TRIOL are used as an input data for calculations of energy distributions. The calculations of frequency f(y) and dose d(y) distributions are performed using own developed programs.

Microdosimetry distributions for 40-200 MeV protons

Theoretical calculations have been performed to obtain microdosimetrical characteristics for protons in energy range from 40 to 200 MeV. This energy range is a representative of proton energies in tissue during radiation therapy and it also represents a large portion of the proton fluency in the South Atlantic Anomaly. Distributions of deposited energy calculated using Monte Carlo track structure code TRIOL and own-made programs were compared with experimental data obtained using spherical tissue-equivalent proportional counter. A good agreement between calculated and experimentally obtained microdosimetry spectra has been found.

Significance of 8-oxoG in the spectrum of DNA damages caused by ionising radiation of different quality

The yields and composition of DNA damages caused by ionising radiation depends on radiation quality. With increasing light energy transfer (LET), the proportion of isolated DNA damages with respect to cluster damaged sites decreases. Non-double strand break complex damages are induced by gamma radiation in mammalian cells at least four times more frequently that prompt DSB. The most important product of oxidative damage to DNA bases is 8-oxo-7,8-dihydroguanine (8-oxoG). The modelling of DNA damage induced by ionising radiation of different qualities was performed to assess frequencies and composition of complex damages containing 8-oxoG. The occurrence of clusters containing 8-oxoG increases from 6 to 11% for LET in the range 0.4-160 keV microm(-1). Distributions of single strand break (SSB) on opposite DNA strand around induced 8-oxoG have similar shape for different ionising radiations, but differ in their occurrence in the whole spectrum of DNA damages. The most probable configuration is a strand break localised at position +/-3 bases from 8-oxoG.

The effect of model approximations on single-collision distributions of low-energy electrons in liquid water

The development of cross sections for the inelastic interaction of low-energy electrons with condensed tissue-like media is best accomplished within the framework of the dielectric theory. In this work we investigate the degree to which various model approximations, used in the above methodology, influence electron single-collision distributions. These distributions are of major importance to Monte Carlo track structure codes, namely, the energy-loss spectrum, the inelastic inverse mean free path, and the ionization efficiency. In particular, we make quantitative assessment of the influence of (1) the optical data set, (2) the dispersion algorithm, and (3) the perturbation and exchange Born corrections. It is shown that, although the shape and position of the energy-loss spectrum remains almost fixed, its peak height may vary by up to a factor of 1.5. Discrepancies in the calculated inelastic inverse mean free path are largely within 20-30% above 100 eV; they increase drastically, though, at lower energies. Exchange and perturbation Born corrections increase gradually below 1 keV leading to a approximately 30 to 40% reduction of the inverse mean free path at 100 eV. The perturbation effect contributes more than the exchange effect to this reduction. Similar to the dispersion situation, the effect of Born corrections at lower energies is also unclear since the models examined disagree strongly below 100 eV. In comparison, the vapor data are higher than the liquid calculations by 20 to 50% as the energy decreases from 1 to 0.1 keV, respectively. The excitation contribution is the main cause of this difference, since the ionization efficiency in the liquid levels off at approximately 90%, whereas the plateau value for the vapor is approximately 70%. It is concluded that electron inelastic distributions for liquid water, although in some respects distinctively different from the vapor phase, have associated uncertainties that are comparable in magnitude to the phase differences. The situation below 100 eV is uncertain.

Monte Carlo simulation of the energy loss of low-energy electrons in liquid water

A Monte Carlo code that performs detailed (i.e. event-by-event) simulation of the transport and energy loss of low-energy electrons (approximately 50-10 000 eV) in water in the liquid phase is presented. The inelastic model for energy loss is based on a semi-empirical dielectric-response function for the valence-shells of the liquid whereas an exchange corrected semi-classical formula was used for K-shell ionization. Following a methodology widely used for the vapour phase, we succeeded in parametrizing the dielectric cross-sections of the liquid in accordance with the Bethe asymptote, thus providing a unified approach for both phases of water and greatly facilitating the computations. Born-corrections at lower energies have been implemented in terms of a second-order perturbation term with a simple Coulomb-field correction and the use of a Mott-type exchange modification. Angular deflections were determined by empirical schemes established from vapour data. Electron tracks generated by the code were used to calculate energy- and interaction-point-kernel distributions at low electron energies in liquid water. The effect of various model assumptions (e.g., dispersion, Born-corrections, phase) on both the single-collision and slowing-down distributions is examined.

La radiolyse de l'eau et des solutions aqueuses : historique et actualité

Experiments showing that water is decomposed by the action of high-energy radiations date back to the first days of the discovery of radioactivity, a century ago. On the occasion of this anniversary, we have attempted to give a comprehensive account of the radiation chemistry of water and its solutions since its origin, with special emphasis on the various physical and chemical stages that led to the present state of this science. To this aim, we describe the effect of different intervening factors on the molecular and radical yields, including dissolved solute concentration, pH, radiation intensity (or dose rate), type and energy of the radiation, presence of oxygen, temperature, phase, and pressure. We also discuss briefly the chemical behavior of the free radicals produced in radiolyzed aqueous solutions. A good, albeit incomplete, description of the phenomena is obtained that leads to various perspectives concerning, on the one hand, the development of this science and, on the other hand, its potential for applications.Key words : radical chemistry, dilution curve, water, hydrated electron, hydroxyl and superoxide radicals, free radicals, radiolysis, chain reactions, molecular and radical yields, cell survival, linear energy transfer.