Cell Survival Computation via the Generalized Stochastic Microdosimetric Model (GSM2); Part II: Numerical Results.

In the present paper we numerically investigate, using Monte Carlo simulation, the theoretical results predicted by the Generalized Stochastic Microdosimetric Model (GSM2), as shown in the published companion paper. Taking advantage of the particle irradiation data ensemble (PIDE) dataset, we calculated GSM2 biological parameters of human salivary gland (HSG) and V79 cell lines. Further, exploiting the TOPAS-microdosimetric extension, we simulated the microdosimetric spectra of different radiation fields of therapeutic interest generated by four different ions (protons, helium-4, carbon-12 and oxygen-16) each at three different residual ranges. We investigated the properties of the initial damage distributions as well as the cell survival curve predicted by GSM2, focusing especially on the non-Poissonian effects naturally included in the model. GSM2 successfully computed cell survival curves, accurately describing experimental behavior even under challenging LET and dose conditions.

Microdosimetric analysis of the radiation quality of two different proton beams in the distal edge of the depth-dose profile.

Proton-therapy exploits the advantageous depth-dose profile of protons to induce the highest damage to tumoral cells in the last millimetres of their range in sharp Bragg Peak. To cover the whole tumoral volume, beams of different energies are combined to create the Spread Out Bragg Peak (SOBP). In passive modulated beams, the energy spread is created with modulators in which the highest energy beam is degraded through different thicknesses of calibrated plastic materials. The highest energy is chosen depending on the deepest point that needs to be treated. This study aims to investigate differences in the radiation quality in the distal edge of SOBP beams with different initial energy and modulation techniques based on microdosimetric measurements with mini Tissue-Equivalent Proportional Counters. The beams investigated are the 62 MeV proton SOBP of the clinical facility of CATANA and the 148 MeV proton SOBP of the research beam line of the proton-therapy centre of Trento.

Interaction of therapeutic12C ions with bone-like targets: physical characterization and dosimetric effect at material interfaces.

Carbon therapy is a promising treatment option for cancer. The physical and biological properties of carbon ions can theoretically allow for the delivery of curative doses to the tumor, while simultaneously limiting risks of toxicity to adjacent healthy structures. The treatment effectiveness can be further improved by decreasing the uncertainties stemming from several sources, including the modeling of tissue heterogeneity. Current treatment plans employ density-based conversion methods to translate patient-specific anatomy into a water system, where dose distribution is calculated. This approach neglects differences in nuclear interactions stemming from the elemental composition of each tissue. In this work, we investigated the interaction of therapeutic carbon ions with bone-like materials. The study concentrated on nuclear interactions and included attenuation curves of 200 and 400 AMeV beams in different types of bones, as well as kinetic energy spectra of all charged fragments produced up to 29 degrees from the beam direction. The comparison between measurements and calculations of the treatment planning system TRiP98 indicated that bone tissue causes less fragmentation of carbon ions than water. Overall, hydrogen and helium particles were found to be the most abundant species, while heavier fragments were mostly detected within 5 degrees from the beam direction. We also investigated how the presence of a soft tissue-bone interface could affect the depth-dose profile. The results revealed a dose spike in the transition region, that extended from the entry channel to the target volume. The findings of this work indicated that the tissue-to-water conversion method based only on density considerations can result in dose inaccuracies. Tissue heterogeneity regions containing bones can potentially produce dose spikes, whose magnitude will depend on the patient anatomy. Dose uncertainties can be decreased by modeling nuclear interactions directly in bones, without applying the tissue-to-water conversion.

Generalized stochastic microdosimetric model: The main formulation.

The present work introduces a rigorous stochastic model, called the generalized stochastic microdosimetric model (GSM^{2}), to describe biological damage induced by ionizing radiation. Starting from the microdosimetric spectra of energy deposition in tissue, we derive a master equation describing the time evolution of the probability density function of lethal and potentially lethal DNA damage induced by a given radiation to a cell nucleus. The resulting probability distribution is not required to satisfy any a priori conditions. After the initial assumption of instantaneous irradiation, we generalized the master equation to consider damage induced by a continuous dose delivery. In addition, spatial features and damage movement inside the nucleus have been taken into account. In doing so, we provide a general mathematical setting to fully describe the spatiotemporal damage formation and evolution in a cell nucleus. Finally, we provide numerical solutions of the master equation exploiting Monte Carlo simulations to validate the accuracy of GSM^{2}. Development of GSM^{2} can lead to improved modeling of radiation damage to both tumor and normal tissues, and thereby impact treatment regimens for better tumor control and reduced normal tissue toxicities.

FLUKA simulation of target fragmentation in proton therapy.

In proton therapy, secondary fragments are created in nuclear interactions of the beam with the target nuclei. The secondary fragments have low kinetic energies and high atomic numbers as compared to primary protons. Fragments have a high LET and deposit all their energy close to the generation point. For their characteristics, secondary fragments can alter the dose distribution and lead to an increase of RBE for the same delivered physical dose. Moreover, the radiobiological impact of target fragmentation is significant mostly in the region before the Bragg peak, where generally healthy tissues are present, and immediately after Bragg peak. Considering the high biological impact of those particles, especially in the case of healthy tissues or organs at risk, the inclusion of target fragmentation processes in the dose calculation of a treatment planning system can be relevant to improve the treatment accuracy and for this reason it is one of the major tasks of the MoVe IT project. In this study, Monte Carlo simulations were employed to fully characterize the mixed radiation field generated by target fragmentation in proton therapy. The dose averaged LET has been evaluated in case of a Spread Out Bragg Peak (SOBP). Starting from LET distribution, RBE has been evaluated with two different phenomenological models. In order to characterize the mixed radiation field, the production cross section has been evaluated by means of the FLUKA code. The future development of present work is to generate a MC database of fragments fluence to be included in TPS.

Microdosimetric measurements as a tool to assess potential in-field and out-of-field toxicity regions in proton therapy.

Relative biological effectiveness (RBE) variations are thought to be one of the primary causes of unexpected normal-tissue toxicities during tumor treatments with charged particles. Unlike carbon therapy, where treatment planning is optimized on the basis of the RBE-weighted dose, a constant RBE value of 1.1 is currently used in proton therapy. Assuming a uniform value can lead to under- or over-dosage, not just to the tumor but also to surrounding normal tissue. RBE changes have been linked with dose/fraction, the biological endpoint and beam properties. Understanding radiation quality and the associated RBE can improve the prediction of normal-tissue toxicities. In this study, we exploited microdosimetry for characterizing radiation quality in proton therapy in-field, and off-beam at 20 (beam edge), 50 (close out-of-field) and 100 (far out-of-field) mm from the beam center. We measured the lineal energy y spectra in a water phantom irradiated with 152 MeV protons, from which beam quality as well as the physical dose could be obtained. Taking advantage of the linear quadratic model and a modified version of the microdosimetric kinetic model, the microdosimetric data were combined with radiobiological parameters (α and ß) of human salivary gland tumor cells for assessing cell survival RBE and RBE-weighted dose. The results indicate that if a dose of 60 Gy is delivered to the peak, the beam edge receives up to 6 Gy while the close and far out-of-field regions receive doses on the order of 10-3 Gy and 10-4 Gy, respectively. The RBE estimate in-beam shows large variations, ranging from 1.0 ± 0.2 at the entrance channel to 2.51 ± 0.15 at the tail. The beam edge follows a similar trend but the RBE calculated at the Bragg peak depth is 2.27 ± 0.17, i.e. twice the RBE in-beam (1.05 ± 0.15). Out-of-field, the estimated RBE is always significantly higher than 1.1 and increases with increasing lateral distance, reaching the overall highest value of 3.4 ± 0.3 at a depth of 206 mm and a lateral distance of 10 mm. The combination of RBE and dose into the biological dose points to the beam edge and the end-of-range in-beam as the areas with the highest risk of potential toxicities.

Transversal dose profile reconstruction for clinical proton beams: A detectors inter-comparison.

PURPOSE: The main purpose of this work is the inter-comparison between different devices devoted to the transversal dose profile recostruction for daily QA tests in proton therapy. METHODS: The results obtained with the EBT3 radiochromic films, used as a reference, and other common quality control devices, have been compared with those obtained with a beam profiling system developed at the "Laboratori Nazionali del Sud" of Italian Institute for Nuclear Physics (INFN-LNS, Catania, Italy). It consists of a plastic scintillator screen (thickness 1 mm), mounted perpendicularly to the beam axis and coupled with a highly sensitive CCD detector in a light-tight box. RESULTS AND CONCLUSION: The tests, carried out both at the INFN-LNS and Trento Proton Therapy Center facilities, show, in general, a good agreement between the different detectors. The beam profiling system, in particular, appears to be a promising quality control device for 2-D relative dosimetry, because of its linear response in a dose rate range useful for proton therapy treatments, its high spatial resolution and its short acquisition and processing time.

Systematic quantification of nanoscopic dose enhancement of gold nanoparticles in ion beams.

High-Z material nanoparticles are being studied as localized dose enhancers in radiotherapeutic applications. Here, the nano-scale physical dose enhancement of proton, carbon and oxygen ion beam radiation by gold nanoparticles was studied by means of Monte Carlo track structure simulation with the TRAX code. We present 2D distributions and radial profiles of the additional dose and the dose enhancement factor for two geometries which consider an isolated and a water-embedded nanoparticle, respectively. Different nanoparticle sizes (radius of 1.2-22 nm) were found to yield qualitatively different absolute and relative dose enhancement distributions and different maximum dose enhancement factors (up to 20). Whereas the smallest nanoparticles produced the highest local dose enhancement factor close to the metal, larger ones led to lower, more diffuse dose enhancement factors that contributed more at larger distances. Differential absorption effects inside the metal were found to be responsible for those characteristics. For the energy range 15-204 MeVu-1, also a mild trend with ion E/A, regardless of the ion species, was found for embedded nanoparticles. In analogy to the width of the ion track itself, slower ions increased the enhancement at the nanoparticle surface. In contrast, no dependence on linear energy transfer was encountered. For slower ions (3-10 MeVu-1), the enhancement effect began to break down over all distances. Finally, the significance of any indirect physical effect was excluded, giving important hints especially in view of the low probabilities (at realistic concentrations and fluences) of direct ion-NP-hits. The very localized nature of the physical dose enhancement found suggests a strong action upon targets closeby, but no relevant effect at cellular distances. When pondering different possible damage enhancement mechanisms of gold nanoparticles in the context of published in vitro and in vivo experimental results, biological pathways are likely to play the key role.

Backscattered electron emission after proton impact on gold nanoparticles with and without polymer shell coating.

This work aims at measuring experimentally proton induced secondary electron energy spectra after interaction with gold nano particles (GNPs) and polymer-coated GNPs. Backscattered electron energy spectra were collected over a 0 to 1000 eV energy range using a retarding field analyzer (RFA). This paper presents the spectra obtained for proton beam energies of 0.5 and 2 MeV and diameter 2.5 and 3.8 nm GNPs. The spectra were also measured for 3.8 nm GNPs after 5 and 10 MeV proton irradiations. GNPs were deposited on a 100 nm carbon film. Each experimental spectrum was compared with dedicated simulations based on existing numerical models used in the TRAX and Geant4 Monte Carlo codes. For 100 nm carbon target, good agreement between experimental, TRAX and Geant4 simulation results can be observed. For 3.8 nm GNPs, the TRAX simulations reproduce with good agreement the electron energy spectra produced after 0.5, 2, 5 and 10 MeV proton irradiations, while Geant4 spectra display a lower secondary electron yield at low energy (<600 eV) for all the studied energies. This underestimation can mostly be explained by the 790 eV threshold applied in the condensed history model used by Geant4 which impacts the secondary electron energy distribution. Results obtained for carbon and gold targets highlight the impact of the secondary electron production threshold for proton ionization process considered in condensed history models. The experimental results demonstrate that the single interaction approach used in TRAX is adapted to reproduce secondary electron emission from GNPs. On the other hand, the standard electron generation threshold implement in G4BetheBlochModel and G4BraggModel condensed-history models used in Geant4 is not adapted to reproduce low energy electron emission in gold targets. Finally, the results highlight that the GNP coating leads to a decrease of the electron yield and mostly affects low energy electrons (<500 eV) emitted from GNPs.

Experimental measurements validate the use of the binary encounter approximation model to accurately compute proton induced dose and radiolysis enhancement from gold nanoparticles.

In protontherapy, it has been suggested that nanoparticles of high-Z material like gold (GNP) could be used as radiosensitizers. The origin of this enhancement phenomenon for proton radiation is not yet well understood and additional mechanistic insights are required. Previous works have highlighted the good capabilities of TRAX to reproduce secondary electron emission from gold material. Therefore, TRAX cross sections obtained with the binary encounter approximation (BEA) model for proton ionization were implemented within Geant4 for gold material. Based on the TRAX cross sections, improved Geant4 simulations have been developed to investigate the energy deposition and radical species production around a spherical gold nanoparticle (5 and 10 nm in diameter) placed in a water volume during proton irradiation. Simulations were performed for incident 2 MeV proton. The dose enhancement factor and the radiolysis enhancement factor were quantified. Results obtained with the BEA model were compared with results obtained with condensed-history models. Experimental irradiation of 200 nm gold films were performed to validate the secondary electron emission reproduction capabilities of physical models used in Monte Carlo (MC) simulations. TRAX simulations reproduced the experimental backscattered electron energy spectrum from gold film with better agreement than Geant4. Results on gold film obtained with the BEA model enabled to estimate the electron emission from GNPs. Results obtained in our study tend to support that the use of the BEA discrete model leads to a significant increase of the dose in the near vicinity of GNPs (<20 nm), while condensed history models used in Geant4 seem to overestimate the dose and the number of chemical species for increasing distances from the GNP. Based on discrete BEA model results, no enhancement effect due to secondary electron emitted from the GNP is expected if the GNP is not in close proximity to key cellular functional elements (DNA, mitochondria…).

Kill painting of hypoxic tumors with multiple ion beams.

We report on a novel method for simultaneous biological optimization of treatment plans for hypoxic tumors using multiple ion species. Our previously introduced kill painting approach, where the overall cell killing is optimized on biologically heterogeneous targets, was expanded with the capability of handling different ion beams simultaneously. The current version (MIBO) of the research treatment planning system TRiP98 has now been augmented to handle 3D (voxel-by-voxel) target oxygenation data. We present a case of idealized geometries where this method can identify optimal combinations leading to an improved peak-to-entrance effective dose ratio. This is achieved by the redistribution of particle fluences, when the heavier ions are preferentially forwarded to hypoxic target areas, while the lighter ions deliver the remaining dose to its normoxic regions. Finally, we present an in silico skull base chordoma patient case study with a combination of 4He and 16O beams, demonstrating specific indications for its potential clinical application. In this particular case, the mean dose, received by the brainstem, was reduced by 3%-5% and by 10%-12% as compared to the pure 4He and 16O plans, respectively. The new method allows a full biological optimization of different ion beams, exploiting the capabilities of actively scanned ion beams of modern particle therapy centers. The possible experimental verification of the present approach at ion beam facilities disposing of fast ion switch is presented and discussed.

STUDY FOR A PASSIVE SCATTERING LINE DEDICATED TO RADIOBIOLOGY EXPERIMENTS AT THE TRENTO PROTON THERAPY CENTER.

The recent worldwide spread of Proton Therapy centers paves the way to new opportunities for basic and applied research related to the use of accelerated proton beams. Clinical centers make use of proton beam energies up to about 230 MeV. This represents an interesting energy range for a large spectrum of applications, including detector testing, radiation shielding and space research. Additionally, radiobiology research might benefit for a larger availability of proton beams, especially in those centers where a room dedicated to research activities also exists. Here, we describe the initial activities for the setup of a radiobiology irradiation facility at the Trento Proton Therapy Center. Data referring to the characterization of the beam in air are essential to that purpose and will be presented. A basic setup for large field irradiation will be also proposed, which is needed for the majority of in vitro and in vivo radiobiology experiments.

Oxygen beams for therapy: advanced biological treatment planning and experimental verification.

Nowadays there is a rising interest towards exploiting new therapeutical beams beyond carbon ions and protons. In particular, [Formula: see text]O ions are being widely discussed due to their increased LET distribution. In this contribution, we report on the first experimental verification of biologically optimized treatment plans, accounting for different biological effects, generated with the TRiP98 planning system with [Formula: see text]O beams, performed at HIT and GSI. This implies the measurements of 3D profiles of absorbed dose as well as several biological measurements. The latter includes the measurements of relative biological effectiveness along the range of linear energy transfer values from ≈20 up to ≈750 keV µ [Formula: see text], oxygen enhancement ratio values and the verification of the kill-painting approach, to overcome hypoxia, with a phantom imitating an unevenly oxygenated target. With the present implementation, our treatment planning system is able to perform a comparative analysis of different ions, according to any given condition of the target. For the particular cases of low target oxygenation, [Formula: see text]O ions demonstrate a higher peak-to-entrance dose ratio for the same cell killing in the target region compared to [Formula: see text]C ions. Based on this phenomenon, we performed a short computational analysis to reveal the potential range of treatment plans, where [Formula: see text]O can benefit over lighter modalities. It emerges that for more hypoxic target regions (partial oxygen pressure of ≈0.15% or lower) and relatively low doses (≈4 Gy or lower) the choice of [Formula: see text]O over [Formula: see text]C or [Formula: see text]He may be justified.

Fragmentation of 120 and 200 MeV u-14He ions in water and PMMA targets.

Recently, the use of 4He particles in cancer radiotherapy has been reconsidered as they potentially represent a good compromise between protons and 12C ions. The first step to achieve this goal is the development of a dedicated treatment planning system, for which basic physics information such as the characterization of the beam lateral scattering and fragmentation cross sections are required. In the present work, the attenuation of 4He primary particles and the build-up of secondary charged fragments at various depths in water and polymethyl methacrylate were investigated experimentally for 120 and 200 MeV u-1 beams delivered by the synchrotron at the Heidelberg Ion-Beam Therapy Center, Heidelberg. Species and isotope identification was accomplished combining energy loss and time-of-flight measurements. Differential yields and energy spectra of all fragments types were recorded between 0° and 20° with respect to the primary beam direction.

Simulations of dose enhancement for heavy atom nanoparticles irradiated by protons.

A possible dose enhancement effect by proton or electron irradiation in the vicinity of nanoparticles consisting of different high Z atomic materials has been investigated using the track structure Monte Carlo code TRAX. In the simulations, Fe, Ag, Gd, Pt and Au nanoparticles (r = 22 and 2 nm) were irradiated with monoenergetic proton beams at energies of therapeutic interest (2, 80 and 300 MeV) and 44 keV electrons. Due to the large number of electrons in atoms with high atomic numbers, many electrons can be released in Auger cascades in addition to the primary ionization process. The potential additional nanoscopic radial dose contributions in the presence of metallic nanoparticles are assessed by comparison with liquid water and water simulated with the same density as the metallic materials. We find a noticeable impact of Auger electrons emitted from the nanoparticles. Special focus has been given to the assessment of complete sets of low-energy electron cross sections for the nanoparticle materials.

Advancing the modeling in particle therapy: from track structure to treatment planning.

We present a series of implementations on Monte Carlo track structure level which might have an impact on treatment planning for particle therapy. We evaluated the effect of multiple ion scattering and radical diffusion on the nanoscopic radial dose. Our cross section database for electron interactions was extended to be able to predict the sensitizing effect of gold nanoparticles in particle therapy. We also implemented LiF as a possible target for efficiency calculations of thermoluminescent detectors (TLDs).

Including oxygen enhancement ratio in ion beam treatment planning: model implementation and experimental verification.

We present a method for adapting a biologically optimized treatment planning for particle beams to a spatially inhomogeneous tumor sensitivity due to hypoxia, and detected e.g., by PET functional imaging. The TRiP98 code, established treatment planning system for particles, has been extended for including explicitly the oxygen enhancement ratio (OER) in the biological effect calculation, providing the first set up of a dedicated ion beam treatment planning approach directed to hypoxic tumors, TRiP-OER, here reported together with experimental tests. A simple semi-empirical model for calculating the OER as a function of oxygen concentration and dose averaged linear energy transfer, generating input tables for the program is introduced. The code is then extended in order to import such tables coming from the present or alternative models, accordingly and to perform forward and inverse planning, i.e., predicting the survival response of differently oxygenated areas as well as optimizing the required dose for restoring a uniform survival effect in the whole irradiated target. The multiple field optimization results show how the program selects the best beam components for treating the hypoxic regions. The calculations performed for different ions, provide indications for the possible clinical advantages of a multi-ion treatment. Finally the predictivity of the code is tested through dedicated cell culture experiments on extended targets irradiation using specially designed hypoxic chambers, providing a qualitative agreement, despite some limits in full survival calculations arising from the RBE assessment. The comparison of the predictions resulting by using different model tables are also reported.

Ionic reactions in He nanodroplets: the [LiHHe]+ complex and its possible energy pathways into products from ab initio calculations.

Ab initio calculations at the multiconfiguration self-consistent field level followed by a multireference configuration interaction were carried out along the two possible collinear approaches of the [LiHHe]+ system, while a three-dimensional calculation of the structures of that complex with LiH+ kept at its equilibrium geometry was also completed at the same level of accuracy. The interaction forces of the lowest two electronic states indicate possible reactive behavior, with the first excited potential-energy surface clearly showing a barrierless path to HeH+ product formation. The details of the reactive pathways and their possible bearing on reaction processes, which could occur at the low temperature of a He nanodroplet holding LiH+ as an impurity, are analyzed and discussed.

Rotational quenching in ionic systems at ultracold temperatures.

The behavior of the rotational quenching for a molecular ion in collision with closed-shell neutral gases is investigated. We confirm that Wigner's threshold law for inelastic scattering holds in the presence of a long-range interaction due to polarization forces decreasing as the inverse fourth power of the distance but find that, because of the contributions of the higher angular momenta, its range of applicability is markedly reduced when compared to the scattering by neutral species. The calculations of the quenching cross sections make evident the special features of ionic systems at ultralow collision energies and yield rate coefficients of the order of 10(-9) x cm(3) x s(-1), much larger than those found for the quenching of neutral molecules.