Spatiotemporal optimisation of prostate intensity modulated proton therapy (IMPT) treatments

Objective. In intensity modulated particle therapy (IMPT), the adoption of spatially and temporally heterogeneous dose distributions allows to decouple the fractionation scheme from the patient anatomy, so that an hypofractionated schedule can be selectively created inside the tumour, while simultaneously exploiting the fractionation effect in the healthy tissues. In this paper, the authors show the reproducibility of the method on a set of prostate patients, quantifying the dependencies of the achievable benefit with respect to conventional and hypofractionated schemes and the sensitivity of the method to setup errors and range uncertainty. Approach. On a cohort of 9 patients, non-uniform IMPT plans were optimised and compared to conventional and hypofractionated schedules. For each patient, the comparison of the three strategies has been based on the output of the cost function used to optimise the treatments. The analysis has been repeated considering different α/β ratios for the tumour, namely 1.5, 3 and 4.5 Gy. For a single patient, setup errors and beam range uncertainty have been analysed: the plans, for each optimisation strategy, have been iteratively forward planned 500 times with randomly varying the patient position in each fraction, and 200 times for systematic range shift. Main results. An average 10% benefit has been shown for the lowest α/β ratio considered for the tumour, where the non-uniform schedule generally converges to hypofractionation; the benefit decreases to 5%–7% for higher α/β ratios, for which the non-uniform schedule always showed better outcomes with respect to the other fractionation schedules. An increased sensitivity to uncertainty, especially for setup errors, has been shown, which can be associated to the spatial non-uniformity of the dose distributions peculiar of the spatiotemporal plans. Significance. This work represents the first investigation of spatiotemporal fractionation for prostate cancer and the beginning of further investigations before clinical implementation can be considered.

TOPAS-nBio validation for simulating water radiolysis and DNA damage under low-LET irradiation

The chemical stage of the Monte Carlo track-structure simulation code Geant4-DNA has been revised and validated. The root-mean-square (RMS) empirical parameter that dictates the displacement of water molecules after an ionization and excitation event in Geant4-DNA has been shortened to better fit experimental data. The pre-defined dissociation channels and branching ratios were not modified, but the reaction rate coefficients for simulating the chemical stage of water radiolysis were updated. The evaluation of Geant4-DNA was accomplished with TOPAS-nBio. For that, we compared predicted time-dependentGvalues in pure liquid water for^·OH, e^-_aq, and H₂with published experimental data. For H₂O₂and H^·, simulation of added scavengers at different concentrations resulted in better agreement with measurements. In addition, DNA geometry information was integrated with chemistry simulation in TOPAS-nBio to realize reactions between radiolytic chemical species and DNA. This was used in the estimation of the yield of single-strand breaks (SSB) induced by¹³⁷Csγ-ray radiolysis of supercoiled pUC18 plasmids dissolved in aerated solutions containing DMSO. The efficiency of SSB induction by reaction between radiolytic species and DNA used in the simulation was chosen to provide the best agreement with published measurements. An RMS displacement of 1.24 nm provided agreement with measured data within experimental uncertainties for time-dependentGvalues and under the presence of scavengers. SSB efficiencies of 24% and 0.5% for^·OH and H^·, respectively, led to an overall agreement of TOPAS-nBio results within experimental uncertainties. The efficiencies obtained agreed with values obtained with published non-homogeneous kinetic model and step-by-step Monte Carlo simulations but disagreed by 12% with published direct measurements. Improvement of the spatial resolution of the DNA damage model might mitigate such disagreement. In conclusion, with these improvements, Geant4-DNA/TOPAS-nBio provides a fast, accurate, and user-friendly tool for simulating DNA damage under low linear energy transfer irradiation.

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

A New Standard DNA Damage (SDD) Data Format

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Monte Carlo simulation of chemistry following radiolysis with TOPAS-nBio

Simulation of water radiolysis and the subsequent chemistry provides important information on the effect of ionizing radiation on biological material. The Geant4 Monte Carlo toolkit has added chemical processes via the Geant4-DNA project. The TOPAS tool simplifies the modeling of complex radiotherapy applications with Geant4 without requiring advanced computational skills, extending the pool of users. Thus, a new extension to TOPAS, TOPAS-nBio, is under development to facilitate the configuration of track-structure simulations as well as water radiolysis simulations with Geant4-DNA for radiobiological studies. In this work, radiolysis simulations were implemented in TOPAS-nBio. Users may now easily add chemical species and their reactions, and set parameters including branching ratios, dissociation schemes, diffusion coefficients, and reaction rates. In addition, parameters for the chemical stage were re-evaluated and updated from those used by default in Geant4-DNA to improve the accuracy of chemical yields. Simulation results of time-dependent and LET-dependent primary yields Gx (chemical species per 100 eV deposited) produced at neutral pH and 25 °C by short track-segments of charged particles were compared to published measurements. The LET range was 0.05-230 keV µm-1. The calculated Gx values for electrons satisfied the material balance equation within 0.3%, similar for protons albeit with long calculation time. A smaller geometry was used to speed up proton and alpha simulations, with an acceptable difference in the balance equation of 1.3%. Available experimental data of time-dependent G-values for [Formula: see text] agreed with simulated results within 7%  ±  8% over the entire time range; for [Formula: see text] over the full time range within 3%  ±  4%; for H2O2 from 49%  ±  7% at earliest stages and 3%  ±  12% at saturation. For the LET-dependent Gx, the mean ratios to the experimental data were 1.11  ±  0.98, 1.21  ±  1.11, 1.05  ±  0.52, 1.23  ±  0.59 and 1.49  ±  0.63 (1 standard deviation) for [Formula: see text], [Formula: see text], H2, H2O2 and [Formula: see text], respectively. In conclusion, radiolysis and subsequent chemistry with Geant4-DNA has been successfully incorporated in TOPAS-nBio. Results are in reasonable agreement with published measured and simulated data.

Comparing stochastic proton interactions simulated using TOPAS-nBio to experimental data from fluorescent nuclear track detectors

Whilst Monte Carlo (MC) simulations of proton energy deposition have been well-validated at the macroscopic level, their microscopic validation remains lacking. Equally, no gold-standard yet exists for experimental metrology of individual proton tracks. In this work we compare the distributions of stochastic proton interactions simulated using the TOPAS-nBio MC platform against confocal microscope data for Al2O3:C,Mg fluorescent nuclear track detectors (FNTDs). We irradiated [Formula: see text] mm3 FNTD chips inside a water phantom, positioned at seven positions along a pristine proton Bragg peak with a range in water of 12 cm. MC simulations were implemented in two stages: (1) using TOPAS to model the beam properties within a water phantom and (2) using TOPAS-nBio with Geant4-DNA physics to score particle interactions through a water surrogate of Al2O3:C,Mg. The measured median track integrated brightness (IB) was observed to be strongly correlated to both (i) voxelized track-averaged linear energy transfer (LET) and (ii) frequency mean microdosimetric lineal energy, [Formula: see text], both simulated in pure water. Histograms of FNTD track IB were compared against TOPAS-nBio histograms of the number of terminal electrons per proton, scored in water with mass-density scaled to mimic Al2O3:C,Mg. Trends between exposure depths observed in TOPAS-nBio simulations were experimentally replicated in the study of FNTD track IB. Our results represent an important first step towards the experimental validation of MC simulations on the sub-cellular scale and suggest that FNTDs can enable experimental study of the microdosimetric properties of individual proton tracks.

Dose enhancement effects to the nucleus and mitochondria from gold nanoparticles in the cytosol

Gold nanoparticles (GNPs) have shown potential as dose enhancers for radiation therapy. Since damage to the genome affects the viability of a cell, it is generally assumed that GNPs have to localise within the cell nucleus. In practice, however, GNPs tend to localise in the cytoplasm yet still appear to have a dose enhancing effect on the cell. Whether this effect can be attributed to stress-induced biological mechanisms or to physical damage to extra-nuclear cellular targets is still unclear. There is however growing evidence to suggest that the cellular response to radiation can also be influenced by indirect processes induced when the nucleus is not directly targeted by radiation. The mitochondrion in particular may be an effective extra-nuclear radiation target given its many important functional roles in the cell. To more accurately predict the physical effect of radiation within different cell organelles, we measured the full chemical composition of a whole human lymphocytic JURKAT cell as well as two separate organelles; the cell nucleus and the mitochondrion. The experimental measurements found that all three biological materials had similar ionisation energies  ∼70 eV, substantially lower than that of liquid water  ∼78 eV. Monte Carlo simulations for 10-50 keV incident photons showed higher energy deposition and ionisation numbers in the cell and organelle materials compared to liquid water. Adding a 1% mass fraction of gold to each material increased the energy deposition by a factor of  ∼1.8 when averaged over all incident photon energies. Simulations of a realistic compartmentalised cell show that the presence of gold in the cytosol increases the energy deposition in the mitochondrial volume more than within the nuclear volume. We find this is due to sub-micron delocalisation of energy by photoelectrons, making the mitochondria a potentially viable indirect radiation target for GNPs that localise to the cytosol.