M.A. Bernal

THE Geant4-DNA project

The Geant4-DNA project proposes to develop an open-source simulation software based and fully included in the general-purpose Geant4 Monte-Carlo simulation toolkit. The main objective of this software is to simulate biological damages induced by ionizing radiations at the cellular and sub-cellular scale. This project was originally initiated by the European Space Agency for the prediction of the deleterious effects of radiations that may affect astronauts during future long duration space exploration missions. In this paper, the Geant4-DNA collaboration presents an overview of the whole on-going project, including its most recent developments that are available in the Geant4 toolkit since December 2009 (release 9.3), as well as an illustration example simulating the direct irradiation of a biological chromatin fiber. Expected extensions involving several research domains, such as particle physics, chemistry and cellular and molecular biology, within a fully interdisciplinary activity of the Geant4 collaboration are also discussed.

Track structure modeling in liquid water: A review of the Geant4-DNA very low energy extension of the Geant4 Monte Carlo simulation toolkit

Understanding the fundamental mechanisms involved in the induction of biological damage by ionizing radiation remains a major challenge of todays radiobiology research. The Monte Carlo simulation of physical, physicochemical and chemical processes involved may provide a powerful tool for the simulation of early damage induction. The Geant4-DNA extension of the general purpose Monte Carlo Geant4 simulation toolkit aims to provide the scientific community with an open source access platform for the mechanistic simulation of such early damage. This paper presents the most recent review of the Geant4-DNA extension, as available to Geant4 users since June 2015 (release 10.2 Beta). In particular, the review includes the description of new physical models for the description of electron elastic and inelastic interactions in liquid water, as well as new examples dedicated to the simulation of physicochemical and chemical stages of water radiolysis. Several implementations of geometrical models of biological targets are presented as well, and the list of Geant4-DNA examples is described.

Monte Carlo simulation of energy-deposit clustering for ions of the same LET in liquid water.

This work presents a Monte Carlo study of energy depositions due to protons, alpha particles and carbon ions of the same linear-energy-transfer (LET) in liquid water. The corresponding track structures were generated using the Geant4-DNA toolkit, and the energy deposition spatial distributions were analyzed using an adapted version of the DBSCAN clustering algorithm. Combining the Geant4 simulations and the clustering algorithm it was possible to compare the quality of the different radiation types. The ratios of clustered and single energy depositions are shown versus particle LET and frequency-mean lineal energies. The estimated effect of these types of radiation on biological tissues is then discussed by comparing the results obtained for different particles with the same LET.

Carbon ion fragmentation effects on the nanometric level behind the Bragg peak depth

In this study, fragmentation yields of carbon therapy beams are estimated using the Geant4 simulation toolkit version 9.5. Simulations are carried out in a step-by-step mode using the Geant4-DNA processes for each of the major contributing fragments. The energy of the initial beam is taken 400 MeV amu(-1) as this is the highest energy, which is used for medical accelerators and this would show the integral role of secondary contributions in radiotherapy irradiations. The obtained results showed that 64% of the global dose deposition is initiated by carbon ions, while up to 36% is initiated by the produced fragments including all their isotopes. The energy deposition clustering yields of each of the simulated fragments are then estimated using the DBSCAN clustering algorithm and they are compared to the yields of the incident primary beam.

Estimation of the RBE of mammography-quality beams using a combination of a Monte Carlo code with a B-DNA geometrical model

The PENELOPE code is used to determine direct strand break yields corresponding to photons from a (60)Co source and 28 and 30 kV x-ray beams impacting on a B-DNA geometrical model, which accounts for five organizational levels of the human genetic material. Direct single, double and total strand break probabilities are determined in a liquid water homogeneous medium with 1.06 g cm(-3) density. The spectra produced by the x-ray beams at various depths in the phantom have been used to study the dependence of the damage yield on the depth. The relative biological effectiveness (RBE) is also estimated using the (60)Co radiation qualities as the reference. According to this work, the damage probabilities and thus the RBE are, within the uncertainties, similar for both x-ray energies and are independent of the depth into the phantom. Furthermore, the total strand break yield is invariant with respect to the energy of the incident photons. The RBE for low-energy x-ray beams determined here (1.3 ± 0.1) is lower than that reported by Kellerer, taking into account that he used a 200 kV radiation as the reference quality. However, our RBE values are consistent with those determined by Kühne et al (2005 Radiat. Res. 164 669-76), which used the same biological endpoint and reference quality as our study. Also, our RBE values are similar to those determined by Verhaegen and Reniers (2004 Radiat. Res. 162 592-9).

The Influence of DNA Configuration on the Direct Strand Break Yield

Purpose. To study the influence of DNA configuration on the direct damage yield. No indirect effect has been accounted for. Methods. The GEANT4-DNA code was used to simulate the interactions of protons and alpha particles with geometrical models of the A-, B-, and Z-DNA configurations. The direct total, single, and double strand break yields and site-hit probabilities were determined. Certain features of the energy deposition process were also studied. Results. A slight increase of the site-hit probability as a function of the incident particle linear energy transfer was found for each DNA configuration. Each DNA form presents a well-defined site-hit probability, independently of the particle linear energy transfer. Approximately 70% of the inelastic collisions and ~60% of the absorbed dose are due to secondary electrons. These fractions are slightly higher for protons than for alpha particles at the same incident energy. Conclusions. The total direct strand break yield for a given DNA form depends weakly on DNA conformation topology. This yield is practically determined by the target volume of the DNA configuration. However, the double strand break yield increases with the packing ratio of the DNA double helix; thus, it depends on the DNA conformation.

Performance of a new atomistic geometrical model of the B-DNA configuration for DNA-radiation interaction simulations

We have recently developed an atomistic model of the B-DNA configuration, up to the 30-nm chromatin fiber. This model is intended to be used in Monte Carlo simulations of the DNA-radiation interaction, specifically in conjunction with the Geant4-DNA extension of the Geant4 Monte Carlo toolkit. In this work, 11449 parallel chromatin fibers have been arranged within a cube mimicking a cell nucleus containing about 6.5×109 base pairs. Each atom in the model is represented by a sphere with the corresponding van der Waals radius. Direct single, double and total DNA strand break yields due to the impact of protons and alpha particles with LET ranging from 4.57 to 207.1 keV/μm have been determined. Also, the corresponding site-hit probabilities have been calculated.

A Feasibility Study of Fricke Dosimetry as an Absorbed Dose to Water Standard for 192Ir HDR Sources

High dose rate brachytherapy (HDR) using 192Ir sources is well accepted as an important treatment option and thus requires an accurate dosimetry standard. However, a dosimetry standard for the direct measurement of the absolute dose to water for this particular source type is currently not available. An improved standard for the absorbed dose to water based on Fricke dosimetry of HDR 192Ir brachytherapy sources is presented in this study. The main goal of this paper is to demonstrate the potential usefulness of the Fricke dosimetry technique for the standardization of the quantity absorbed dose to water for 192Ir sources. A molded, double-walled, spherical vessel for water containing the Fricke solution was constructed based on the Fricke system. The authors measured the absorbed dose to water and compared it with the doses calculated using the AAPM TG-43 report. The overall combined uncertainty associated with the measurements using Fricke dosimetry was 1.4% for k = 1, which is better than the uncertainties reported in previous studies. These results are promising; hence, the use of Fricke dosimetry to measure the absorbed dose to water as a standard for HDR 192Ir may be possible in the future.

Evaluation of the mean energy deposit during the impact of charged particles on liquid water

The DNA strand break yield due to the impact of ionizing particles on living beings is closely related to the number of inelastic events per unit absorbed dose produced by these particles. The higher this number, the higher the probability of causing DNA strand breaks per unit absorbed dose. In a previous work, it was found that the total number of events produced by primary particles and the secondary electrons is almost independent of the type and energy of the incident particle (or LET). This finding could be supported by a quasi-constant mean energy deposit by inelastic event (ε). In this work, ε was defined and determined for electrons and the non-negative charge states of hydrogen (H⁰,⁺) and helium (He⁰,⁺,²⁺) species impacting on liquid water. Ionization, excitation and charge transfer (up to two-electron transfers) processes have been included in present calculations. We found that, for liquid water, ε is within 13.7 ± 4.1 eV, 14.2 ± 1.7 eV and 13.8 ± 1.4 eV for electrons, hydrogen and helium species, respectively, with impact energies changing over three orders of magnitude. Unlike the mean excitation energy, the mean energy deposit per inelastic event depends not only on the target molecule but also on the projectile features. However, this dependence is relatively weak. This fact supports the quasi-independent number of inelastic events per unit absorbed dose found previously when charged particles impact on matter.

Dosimetric evaluation of radionuclides for VCAM-1-targeted radionuclide therapy of early brain metastases.

Brain metastases develop frequently in patients with breast cancer, and present a pressing therapeutic challenge. Expression of vascular cell adhesion molecule 1 (VCAM-1) is upregulated on brain endothelial cells during the early stages of metastasis and provides a target for the detection and treatment of early brain metastases. The aim of this study was to use a model of early brain metastasis to evaluate the efficacy of α-emitting radionuclides, 149Tb, 211At, 212Pb, 213Bi and 225Ac; β-emitting radionuclides, 90Y, 161Tb and 177Lu; and Auger electron (AE)-emitters 67Ga, 89Zr, 111In and 124I, for targeted radionuclide therapy (TRT). METHODS: Histologic sections and two photon microscopy of mouse brain parenchyma were used to inform a cylindrical vessel geometry using the Geant4 general purpose Monte Carlo (MC) toolkit with the Geant4-DNA low energy physics models. Energy deposition was evaluated as a radial function and the resulting phase spaces were superimposed on a DNA model to estimate double-strand break (DSB) yields for representative β- and α-emitters, 177Lu and 212Pb. Relative biological effectiveness (RBE) values were determined by only evaluating DNA damage due to physical interactions. RESULTS: 177Lu produced 2.69 ± 0.08 DSB per GbpGy, without significant variation from the lumen of the vessel to a radius of 100 µm. The DSB yield of 212Pb included two local maxima produced by the 6.1 MeV and 8.8 MeV α-emissions from decay products, 212Bi and 212Po, with yields of 7.64 ± 0.12 and 9.15 ± 0.24 per GbpGy, respectively. Given its higher DSB yield 212Pb may be more effective for short range targeting of early micrometastatic lesions than 177Lu. CONCLUSION: MC simulation of a model of early brain metastases provides invaluable insight into the potential efficacy of α-, β- and AE-emitting radionuclides for TRT. 212Pb, which has the attributes of a theranostic radionuclide since it can be used for SPECT imaging, showed a favorable dose profile and RBE.