A Monte Carlo study for the calculation of the average linear energy transfer (LET) distributions for a clinical proton beam line and a radiobiological carbon ion beam line

A linear energy transfer distributions computation method for inhomogeneous medium by using the water equivalent ratio approximation

Dose-averaged linear energy transfer (LET), LETd is widely used in proton therapy. Compared with analytical models, Monte Carlo (MC) simulations are more accurate in obtaining LETd distributions, but they are time-consuming. This study used the 3D LETd distributions of proton beam spots in water by MC simulations as a benchmark data set. Subsequently, by combining the water equivalent ratio of various human tissues, the 3D LETd distributions of clinical cases could be quickly obtained. Our method was applied to a single spot of 160 MeV proton beam in a water-bone phantom and a pelvic case. We also computed the 3D LETd distributions for multiple proton beam spots in the pelvic case and a lung case. The results of our method were compared with the results of MC simulations, demonstrating that our method can rapidly provide 3D LETd distributions of clinical cases with acceptable differences from MC simulations.

Sensitivity of a mini-TEPC to radiation quality variations in clinical proton beams

PURPOSE

This work aims at studying the sensitivity of a miniaturized Tissue-Equivalent Proportional Counter to variations of beam quality in clinical radiation fields, to further investigate its performances as radiation quality monitor.

METHODS

Measurements were taken at the CATANA facility (INFN-LNS, Catania, Italy), in a monoenergetic and an energy-modulated proton beam with the same initial energy of 62 MeV. PMMA layers were placed in front of the detector to measure at different depths along the depth-dose profile. The frequency- and dose-mean lineal energy were compared to the track- and dose-averaged LET calculated by Monte Carlo simulations. A microdosimetric evaluation of the Relative Biological Effectiveness (RBE) was performed and compared with cell survival experiments.

RESULTS

Microdosimetric distributions measured at identical depths in the two beams show spectral differences reflecting their different radiation quality. Discrepancies are most evident at depths corresponding to the Spread-Out Bragg Peak, while spectra at the entrance and in the dose fall-off regions are similar. This can be explained by the different energy components that compose the pristine and spread-out peaks at each depth. The trend of microdosimetric mean values matches that of calculated LET averages along the entire penetration depth, and the microdosimetric estimation of RBE is consistent with radiobiological data not only at 2 Gy but also at lower dose levels, such as those absorbed by healthy tissues.

CONCLUSIONS

The mini-TEPC is sensitive to differences in radiation quality resulting from different modulations of the proton beam, confirming its potential for beam quality monitoring in proton therapy.

Computational approaches in the estimation of radiobiological damage for human-malignant cells irradiated with clinical proton and carbon beams

PURPOSE

The use of Monte Carlo (MC) simulations capable of reproducing radiobiological effects of ionising radiation on human cell lines is of great importance, especially for cases involving protons and heavier ion beams. In the latter, huge uncertainties can arise mainly related to the effects of the secondary particles produced in the beam-tissue interaction. This paper reports on a detailed MC study performed using Geant4-based approach on three cancer cell lines, the HTB-177, CRL-5876 and MCF-7, that were previously irradiated with therapeutic proton and carbon ion beams.

METHODS

A Geant4-based approach used jointly with analytical calculations has been developed to provide a more realistic estimation of the radiobiological damage produced by proton and carbon beams in tissues, reproducing available data obtained from in vitro cell irradiations. The MC "Hadrontherapy" Geant4 application and the Local Effect Model: LEM I, LEM II and LEM III coupled with the different numerical approaches: RapidRusso (RR) and RapidScholz (RS) were used in the study.

RESULTS

Experimental survival curves are compared with those evaluated using the highlighted Geant4 MC-based approach via chi-square statistical analysis, for the combinations of radiobiological models and numerical approaches, as outlined above.

CONCLUSION

This study has presented a comparison of the survival data from MC simulations to experimental survival data for three cancer cell lines. An overall best level of agreement was obtained for the HTB-177 cells.

Uncertainty in the range of the protons and C-ions in particle therapy due to a hydration level of a human body model

Cancer treatment with protons and carbon ions relies on the property of the accelerated charged particles to deposit most of their energy in the vicinity of their range (around the Bragg peak). The level of hydration in a cancer patient's body may vary within hours. Some patients may be heavy to moderately dehydrated, and some may be well and even excessively hydrated. In this research, we aim to estimate the uncertainty of the protons and C-ion ranges because of the different hydration levels of the human body. For the study of the impact of body hydration level on the particle's ranges, we have designed a new phantom model - a homogeneous mixture of an Average HUuman BOdy constituting elements (AHUBO) in three states of hydration: normal (n), dehydrated (d), and excessively hydrated (e) by applying corresponding recalibration in the "atomic-stoichiometry model" due to the water sufficiency/deficiency. The purpose of the study is to estimate the shift in the ranges depending on the hydration level, possibly suggest particle beam energy adjustments to overcome the range uncertainties, to deliver the prescribed dose to the tumour while sparing the healthy tissue. Herein we present the results of the FLUKA-Flair simulations of the therapeutic range of energies of protons (50-105 MeV) and C-ions (30-210 MeV) respectively, into an AHUBO head phantom model at three levels of hydration (normal, dehydrated, and excessively hydrated). The range uncertainty was estimated via the shifts of the Bragg-peaks position for the three different hydration levels. The estimations showed that the range uncertainty (ΔR) due to body hydration for the maximum energy in the range for protons at 105 MeV is about 0.04 mm and for C-ions at 190 MeV/u is about 0.06 mm.

An iterative convex relaxation method for proton LET optimization

Objective:A constant relative biological effectiveness of 1.1 in current clinical practice of proton radiotherapy (RT) is a crude approximation and may severely underestimate the biological dose from proton RT to normal tissues, especially near the treatment target at the end of Bragg peaks that exhibits high linear energy transfer (LET). LET optimization can account for biological effectiveness of protons during treatment planning, for minimizing biological proton dose and hot spots to normal tissues. However, the LET optimization is usually nonlinear and nonconvex to solve, for which this work will develop an effective optimization method based on iterative convex relaxation (ICR).Approach: In contrast to the generic nonlinear optimization method, such as Quasi-Newton (QN) method, that does not account for specific characteristics of LET optimization, ICR is tailored to LET modeling and optimization in order to effectively and efficiently solve the LET problem. Specifically, nonlinear dose-averaged LET term is iteratively linearized and becomes convex during ICR, while nonconvex dose-volume constraint and minimum-monitor-unit constraint are also handled by ICR, so that the solution for LET optimization is obtained by solving a sequence of convex and linearized convex subproblems. Since the high LET mostly occurs near the target, a 1 cm normal-tissue expansion of clinical target volume (CTV) (excluding CTV), i.e. CTV1cm, is defined to as an auxiliary structure during treatment planning, where LET is minimized.Main results: ICR was validated in comparison with QN for abdomen, lung, and head-and-neck cases. ICR was effective for LET optimization, as ICR substantially reduced the LET and biological dose in CTV1cm the ring, with preserved dose conformality to CTV. Compared to QN, ICR had smaller LET, physical and biological dose in CTV1cm, and higher conformity index values; ICR was also computationally more efficient, which was about 3 times faster than QN.Significance: A LET-specific optimization method based on ICR has been developed for solving proton LET optimization, which has been shown to be more computationally efficient than generic nonlinear optimizer via QN, with better plan quality in terms of LET, biological and physical dose conformality.

Dosimetric comparison of robust angles in carbon-ion radiation therapy for prostate cancer

The objective of this study is to compare the plan robustness at various beam angles. Hence, the influence of the beam angles on robustness and linear energy transfer (LET) was evaluated in gantry-based carbon-ion radiation therapy (CIRT) for prostate cancer. 10 patients with prostate cancer were considered, and a total dose of 51.6 Gy (Relative biological effectiveness (RBE) was prescribed for the target volume in 12 fractions. Five beam field plans comprising two opposed fields with different angle pairs were characterized. Further, dose parameters were extracted, and the RBE-weighted dose and LET values for all angle pairs were compared. All plans considering the setup uncertainty satisfied the dose regimen. When a parallel beam pair was used for perturbed scenarios to take into account set-up uncertainty in the anterior direction, the LET clinical treatment volume (CTV) D _95% standard deviation was 1.5 times higher, and the standard deviation of RBE-weighted CTV D _95% was 7.9 times higher compared to an oblique pair. The oblique beam fields were superior in terms of dose sparing for the rectum compared to the dose distribution using two conventional lateral opposed fields for prostate cancer.

A phenomenological model of proton FLASH oxygen depletion effects depending on tissue vasculature and oxygen supply

Introduction

Radiation-induced oxygen depletion in tissue is assumed as a contributor to the FLASH sparing effects. In this study, we simulated the heterogeneous oxygen depletion in the tissue surrounding the vessels and calculated the proton FLASH effective-dose-modifying factor (FEDMF), which could be used for biology-based treatment planning.

Methods

The dose and dose-weighted linear energy transfer (LET) of a small animal proton irradiator was simulated with Monte Carlo simulation. We deployed a parabolic partial differential equation to account for the generalized radiation oxygen depletion, tissue oxygen diffusion, and metabolic processes to investigate oxygen distribution in 1D, 2D, and 3D solution space. Dose and dose rates, particle LET, vasculature spacing, and blood oxygen supplies were considered. Using a similar framework for the hypoxic reduction factor (HRF) developed previously, the FEDMF was derived as the ratio of the cumulative normoxic-equivalent dose (CNED) between CONV and UHDR deliveries.

Results

Dynamic equilibrium between oxygen diffusion and tissue metabolism can result in tissue hypoxia. The hypoxic region displayed enhanced radio-resistance and resulted in lower CNED under UHDR deliveries. In 1D solution, comparing 15 Gy proton dose delivered at CONV 0.5 and UHDR 125 Gy/s, 61.5% of the tissue exhibited ≥20% FEDMF at 175 μm vasculature spacing and 18.9 μM boundary condition. This percentage reduced to 34.5% and 0% for 8 and 2 Gy deliveries, respectively. Similar trends were observed in the 3D solution space. The FLASH versus CONV differential effect remained at larger vasculature spacings. A higher FLASH dose rate showed an increased region with ≥20% FEDMF. A higher LET near the proton Bragg peak region did not appear to alter the FLASH effect.

Conclusion

We developed 1D, 2D, and 3D oxygen depletion simulation process to obtain the dynamic HRF and derive the proton FEDMF related to the dose delivery parameters and the local tissue vasculature information. The phenomenological model can be used to simulate or predict FLASH effects based on tissue vasculature and oxygen concentration data obtained from other experiments.

4He dose- and track-averaged linear energy transfer: Monte Carlo algorithms and experimental verification

Objective. In the present hadrontherapy scenario, there is a growing interest in exploring the capabilities of different ion species other than protons and carbons. The possibility of using different ions paves the way for new radiotherapy approaches, such as the multi-ions treatment, where radiation could vary according to target volume, shape, depth and histologic characteristics of the tumor. For these reasons, in this paper, the study and understanding of biological-relevant quantities was extended for the case of 4He ion. Approach. Geant4 Monte Carlo based algorithms for dose- and track-averaged LET (Linear Energy Transfer) calculations, were validated for 4He ions and for the case of a mixed field characterised by the presence of secondary ions from both target and projectile fragmentation. The simulated dose and track averaged LETs were compared with the corresponding dose and frequency mean values of the lineal energy, y D ¯ and y ¯ F , derived from experimental microdosimetric spectra. Two microdosimetric experimental campaigns were carried out at the Italian eye proton therapy facility of the Laboratori Nazionali del Sud of Istituto Nazionale di Fisica Nucleare (INFN-LNS, Catania, I) using two different microdosimeters: the MicroPlus probe and the nano-TEPC (Tissue Equivalent Proportional Counter). Main results. A good agreement of L ¯ d Total and L ¯ t Total with y ¯ D and y ¯ T experimentally measured with both microdosimetric detectors MicroPlus and nano-TEPC in two configurations: full energy and modulated 4He ion beam, was found. Significance. The results of this study certify the use of a very effective tool for the precise calculation of LET, given by a Monte Carlo approach which has the advantage of allowing detailed simulation and tracking of nuclear interactions, even in complex clinical scenarios.

A systematic study of the contribution of counting statistics to the final lineal energy uncertainty in microdosimetry

Objectives. Microdosimetry is proving to be a reliable and powerful tool to be applied in different fields such as radiobiology, radiation protection and hadron therapy. However, accepted standard protocols and codes of practice are still missing. With this regard, a systematic and methodical uncertainty analysis is fundamental to build an accredited uncertainty budget of practical use. This work studied the contribution of counting statistics (i.e. number of events collected) to the final frequency-mean and dose-mean lineal energy uncertainties, aiming at providing guidelines for good experimental and simulation practice. The practical limitation of current technologies and the non-negligible probability of nuclear reactions require careful considerations and nonlinear approaches. Approach. Microdosimetric data were obtained by means of the particle tracking Monte Carlo code Geant4. The uncertainty analysis was carried out relying on a Monte Carlo based numerical analysis, as suggested by the BIPM's ‘Guide to the expression of uncertainty in measurement’. Final uncertainties were systematically investigated for proton, helium and carbon ions at an increasing number of detected events, for a range of different clinical-relevant beam energies. Main results. Rare events generated by nuclear interactions in the detector sensitive volume were found to massively degrade microdosimetric uncertainties unless a very high statistics is collected. The study showed an increasing impact of such events for increasing beam energy and lighter ions. For instance, in the entrance region of a 250 MeV proton beam, about 5 ∗ 107 events need to be collected to obtain a dose-mean lineal energy uncertainty below 10%. Significance. The results of this study help define the necessary conditions to achieve appropriate statistics in computational microdosimetry, pointing out the importance of properly taking into account nuclear interaction events. Their impact on microdosimetric quantities and on their uncertainty is significant and cannot be overlooked, particularly when characterising clinical beams and radiobiological response. This work prepared the ground for deeper investigations involving dedicated experiments and for the development of a method to properly evaluate the counting statistics uncertainty contribution in the uncertainty budget, whose accuracy is fundamental for the clinical transition of microdosimetry.

Towards harmonizing clinical linear energy transfer (LET) reporting in proton radiotherapy: a European multi-centric study

BACKGROUND

Clinical data suggest that the relative biological effectiveness (RBE) in proton therapy (PT) varies with linear energy transfer (LET). However, LET calculations are neither standardized nor available in clinical routine. Here, the status of LET calculations among European PT institutions and their comparability are assessed.

MATERIALS AND METHODS

Eight European PT institutions used suitable treatment planning systems with their center-specific beam model to create treatment plans in a water phantom covering different field arrangements and fulfilling commonly agreed dose objectives. They employed their locally established LET simulation environments and procedures to determine the corresponding LET distributions. Dose distributions D1.1 and DRBE assuming constant and variable RBE, respectively, and LET were compared among the institutions. Inter-center variability was assessed based on dose- and LET-volume-histogram parameters.

RESULTS

Treatment plans from six institutions fulfilled all clinical goals and were eligible for common analysis. D1.1 distributions in the target volume were comparable among PT institutions. However, corresponding LET values varied substantially between institutions for all field arrangements, primarily due to differences in LET averaging technique and considered secondary particle spectra. Consequently, DRBE using non-harmonized LET calculations increased inter-center dose variations substantially compared to D1.1 and significantly in mean dose to the target volume of perpendicular and opposing field arrangements (p < 0.05). Harmonizing LET reporting (dose-averaging, all protons, LET to water or to unit density tissue) reduced the inter-center variability in LET to the order of 10-15% within and outside the target volume for all beam arrangements. Consequentially, inter-institutional variability in DRBE decreased to that observed for D1.1.

CONCLUSION

Harmonizing the reported LET among PT centers is feasible and allows for consistent multi-centric analysis and reporting of tumor control and toxicity in view of a variable RBE. It may serve as basis for harmonized variable RBE dose prescription in PT.

DNA double-strand breaks in cancer cells as a function of proton linear energy transfer and its variation in time

PURPOSE

The complex relationship between linear energy transfer (LET) and cellular response to radiation is not yet fully elucidated. To better characterize DNA damage after irradiations with therapeutic protons, we monitored formation and disappearance of DNA double-strand breaks (DNA DSB) as a function of LET and time. Comparisons with conventional γ-rays and high LET carbon ions were also performed.

MATERIALS AND METHODS

In the present work, we performed immunofluorescence-based assay to determine the amount of DNA DSB induced by different LET values along the 62 MeV therapeutic proton Spread out Bragg peak (SOBP) in three cancer cell lines, i.e. HTB140 melanoma, MCF-7 breast adenocarcinoma and HTB177 non-small lung cancer cells. Time dependence of foci formation was followed as well. To determine irradiation positions, corresponding to the desired LET values, numerical simulations were carried out using Geant4 toolkit. We compared γ-H2AX foci persistence after irradiations with protons to that of γ-rays and carbon ions.

RESULTS

With the rise of LET values along the therapeutic proton SOBP, the increase of γ-H2AX foci number is detected in the three cell lines up to the distal end of the SOBP, while there is a decrease on its distal fall-off part. With the prolonged incubation time, the number of foci gradually drops tending to attain the residual level. For the maximum number of DNA DSB, irradiation with protons attain higher level than that of γ-rays. Carbon ions produce more DNA DSB than protons but not substantially. The number of residual foci produced by γ-rays is significantly lower than that of protons and particularly carbon ions. Carbon ions do not produce considerably higher number of foci than protons, as it could be expected due to their physical properties.

CONCLUSIONS

In situ visualization of γ-H2AX foci reveal creation of more lesions in the three cell lines by clinically relevant proton SOBP than γ-rays. The lack of significant differences in the number of γ-H2AX foci between the proton and carbon ion-irradiated samples suggests an increased complexity of DNA lesions and slower repair kinetics after carbon ions compared to protons. For all three irradiation types, there is no major difference between the three cell lines shortly after irradiations, while later on, the formation of residual foci starts to express the inherent nature of tested cells, therefore increasing discrepancy between them.

Measurement of therapeutic12C beam in a water phantom using CR-39

The motivation for this study was to explore a new method to test the particle spatial distribution for a therapeutic carbon beam. CR-39 plastic nuclear track detectors were irradiated to a 276.5 MeV u^-1mono-energy carbon beam at the heavy ion facility in the Shanghai Proton and Heavy Ion Center. The spatial distribution of the primary carbon beam and secondary fragments in a water phantom were systematically analyzed both in the transverse direction (perpendicular to the projection direction of the primary beam) and at different depths in the longitudinal direction (along the projection direction of the primary beam) with measured tracks on the CR-39 detectors. Meanwhile, the theoretically spatial distribution and linear energy transfer (LET) spectra of the primary beam and secondary fragments were calculated using the Monte Carlo (MC) toolkit Geant4. The results showed that the CR-39 detectors are capable of providing high lateral resolution of carbon ion at different depths. In the range of the primary carbon beam, the beam width simulated with MC is in good agreement with that of experimental measurement. The track size registered in the CR-39 has a good correlation with the particle LET. These findings indicate that the CR-39 can be used for measuring both the particle flux and its spatial distribution of carbon ions.

Evaluations of a flat-panel based compact daily quality assurance device for proton pencil beam scanning (PBS) system

PURPOSE

To evaluate the flat-panel detector quenching effect and clinical usability of a flat-panel based compact QA device for PBS daily constancy measurements.

MATERIALS & METHOD

The QA device, named Sphinx Compact, is composed of a 20x20 cm² flat-panel imager mounted on a portable frame with removable plastic modules for constancy checks of proton energy (100 MeV, 150 MeV, 200 MeV), Spread-Out-Bragg-Peak (SOBP) profile, and machine output. The potential quenching effect of the flat-panel detector was evaluated. Daily PBS QA tests of X-ray/proton isocenter coincidence, the constancy of proton spot position and sigma as well as the energy of pristine proton beam, and the flatness of SOBP proton beam through the 'transformed' profile were performed and analyzed. Furthermore, the sensitivity of detecting energy changes of pristine proton beam was also evaluated.

RESULTS

The quenching effect was observed at depths near the pristine peak regions. The flat-panel measured range of the distal 80% is within 0.9 mm to the defined ranges of the delivered proton beams. X-ray/proton isocenter coincidence tests demonstrated maximum mismatch of 0.3 mm between the two isocenters. The device can detect 0.1 mm change of spot position and 0.1 MeV energy changes of pristine proton beams. The measured transformed SOBP beam profile through the wedge module rendered as flat.

CONCLUSIONS

Even though the flat-panel detector exhibited quenching effect at the Bragg peak region, the proton range can still be accurately measured. The device can fulfill the requirements of the daily QA tests recommended by the AAPM TG224 Report.

Standardizing Monte Carlo simulation parameters for a reproducible dose-averaged linear energy transfer

OBJECTIVE

Dose-averaged linear energy transfer (LETD) is one of the factors which determines relative biological effectiveness (RBE) for treatment planning in proton therapy. It is usually determined from Monte Carlo (MC) simulation. However, no standard simulation protocols were established for sampling of LETD. Simulation parameters like maximum step length and range cut will affect secondary electrons production and have an impact on the accuracy of dose distribution and LETD. We aim to show how different combinations of step length and range cut in GEANT4 will affect the result in sampling of LETD using different MC scoring methods.

METHODS

In this work, different step length and range cut value in a clinically relevant voxel geometry were used for comparison. Different LETD scoring methods were established and the concept of covariance between energy deposition per step and step length is used to explain the differences between them.

RESULTS

We recommend a maximum step length of 0.05 mm and a range cut of 0.01 mm in MC simulation as this yields the most consistent LETD value across different scoring methods. Different LETD scoring methods are also compared and variation up to 200% can be observed at the plateau of 80 MeV proton beam. Scoring Method one has one of the lowest percentage differences compared across all simulation parameters.

CONCLUSION

We have determined a set of maximum step length and range cut parameters to be used for LETD scoring in a 1 mm voxelized geometry. LETD scoring method should also be clearly defined and standardized to facilitate cross-institutional studies.

ADVANCES IN KNOWLEDGE

Establishing a standard simulation protocol for sampling LETD would reduce the discrepancy when comparing data across different centres, and this can improve the calculation for RBE.

On the concepts of dose-mean lineal energy, unrestricted and restricted dose-averaged LET in proton therapy

To calculate 3D distributions of microdosimetric-based restricted dose-averaged LET (LETd) and dose-mean lineal energy ([Formula: see text]) in order to explore their similarities and differences between each other and with the traditional unrestricted LETd. Additionally, a new expression for optimum restricted LETd calculation is derived, allowing for disregarding straggling-associated functions in the classical microdosimetric theory. Restricted LETd and [Formula: see text] for polyenergetic beams can be obtained by integrating previously developed energy-dependent microdosimetric functions over the energetic spectrum of these beams. This calculation is extended to the entire calculation volume using an algorithm to determine spectral fluence. Equivalently, unrestricted LETd can be obtained integrating the stopping power curve on the spectrum. A new expression to calculate restricted LETd is also derived. Results for traditional and new formulas are compared for a clinical 100 MeV proton beam. Distributions of unrestricted LETd, restricted LETd and [Formula: see text] are analyzed for a prostate case, for microscopic spherical sites of 1 µm and 10 µm in diameter. Traditional and new expressions for restricted LETd remarkably agree, being the mean differences 0.05  ±  0.04 keV µm^-1 for the 1 µm site and 0.05  ±  0.02 keV µm^-1 for the 10 µm site. In the prostate case, the ratio between the maximum and the central value for central axis (CAX) profiles is around 2 for all the quantities, being the highest for restricted LETd for 1 µm (2.17) and the lowest for [Formula: see text] for 1 µm (1.78). Unrestricted LETd, restricted LETd and [Formula: see text] can be analytically computed and compared for clinical plans. Two important consequences of the calculation of [Formula: see text] are: (1) its distribution can be verified by directly measuring it in clinical beams; and (2), optimization of proton treatments based on these quantities is enabled as well as future developments of RBE models based on them.

Optimization of treatment planning for hypoxic tumours and re-modulation of radiation intensity in heavy-ion radiotherapy

Aim

The purpose of this study is to optimize treatment planning in carbon ion radiotherapy, taking into account the effect of tumour hypoxia.

Background

In conventional hadron therapy, the goal is to create a homogenous dose in the tumour area and, thus, achieve a uniform cell survival level. Since the induction of a specific damage to cells is directly influenced by the level of hypoxia in the tissue, the varying oxygen pressure in the different regions of hypoxic tumours would disrupt the uniformity of the cell survival level.

Materials and methods

Using the Geant4 Monte Carlo Code, the physical dose profile and dose-averaged linear energy transfer were calculated in the tumour. Then, the oxygen enhancement ratio in different areas of the tumour were compared with different pressures.

Results

Modulations of radiation intensities as well as energies of ion beams were calculated, both considering and disregarding the effect of hypoxia, and the required dose profiles were compared with each other. Cell survival levels were also compared between the two methods. An equation was obtained for re-modulating the beams in the presence of hypoxia, and radiation weighting factors were extracted for the beam intensities.

Conclusion

The results show that taking the effect of hypoxia into account would cause the reduction of average doses delivered to the tumour tissues up to 1.54 times. In this regard, the required dose is reduced by 1.63 times in the healthy tissues before the tumour. This will result in an effective protection of healthy tissues around the tumour.

Microdosimetry at the CATANA 62 MeV proton beam with a sealed miniaturized TEPC

A new mini-TEPC with cylindrical sensitive volume of 0.9 mm in diameter and height, and with external diameter of 2.7 mm, has been developed to work without gas flow. With such a mini counter we have measured the physical quality of the 62 MeV therapeutic proton beam of CATANA (Catania, Italy). Measurements were performed at six precise positions along the Spread-Out Bragg Peak (SOBP): 1.4, 19.4, 24.6, 29.0, 29.7 and 30.8 mm, corresponding to positions of clinical relevance (entrance, proximal, central, and distal-edge of the SOBP) or of high lineal energy transfer (LET) increment (distal-dose drop off). Without refilling the microdosimeter with new gas, the measurements were repeated at the same positions 4 months later, in order to study the stability of the response in sealed-mode operation. From the microdosimetric spectra the frequency-mean lineal energy y-_F and the dose-mean lineal energy y-_D were derived and compared with average LET values calculated by means of Geant4 simulations. The comparison points out, in particular, a good agreement between microdosimetric y-_D and the total dose-average LET¯_d, which is the average LET of the mixed radiation field, including the contribution by nuclear reactions.

MICRODOSIMETRY OF CLINICAL ION BEAMS: CONVERTING SPECTRA FROM DIAMOND SLAB TO WATER OF DIFFERENT SHAPES

Microdosemeters are frequently used today to specify the radiation quality in the framework of ion-beam therapy. The heterogeneity of the detector shapes and the materials limits the possibility of comparing directly spectra and mean lineal energies. A method was recently studied to convert the spectra obtained with unidirectional ion beams in slab detectors to those obtained with detectors of different in shape and material. The method is based on the observation that the lineal-energy spectra of slab detector, in a restricted energy interval, approximate the Linear Energy Transfer distributions at corresponding material and particle type and energies. In this study, the experimental spectra collected with a slab diamond detector are converted to the spectra that would be obtained using water detectors of spherical and cylindrical shapes.

Biological outcomes of γ-radiation induced DNA damages in breast and lung cancer cells pretreated with free radical scavengers

PURPOSE

Investigation of effects on DNA of γ-irradiated human cancer cells pretreated with free radical scavengers is aimed to create reference data which would enable assessment of the relative efficiency of high linear energy transfer (LET) radiations used in hadron therapy, i.e. protons and carbon ions.

MATERIALS AND METHODS

MCF-7 breast and HTB177 lung cancer cells are irradiated with γ-rays. To minimize indirect effects of irradiation, dimethyl sulfoxide (DMSO) or glycerol are applied as free radical scavengers. Biological response to irradiation is evaluated through clonogenic cell survival, immunocytochemical and cell cycle analysis, as well as expression of proteins involved in DNA damage response.

RESULTS

Examined cell lines reveal similar level of radioresistance. Application of scavengers leads to the rise of cell survival and decreases the number of DNA double strand breaks in irradiated cells. Differences in cell cycle and protein expression between the two cell lines are probably caused by different DNA damage repair mechanisms that are activated.

CONCLUSION

The obtained results show that DMSO and glycerol have good scavenging capacity, and may be used to minimize DNA damage induced by free radicals. Therefore, they will be used as the reference for comparison with high LET irradiations, as well as good experimental data suitable for validation of numerical simulations.

Computational models and tools

In this chapter, we describe two different methods, analytical (pencil beam) algorithms and Monte Carlo simulations, used to obtain the intended dose distributions in patients and evaluate their strengths and shortcomings. We discuss the difference between the prescribed physical dose and the biologically effective dose, the relative biological effectiveness (RBE) between ions and photons and the dependence of RBE on the linear energy transfer (LET). Lastly, we show how LET- or RBE-based optimization can be used to improve treatment plans and explore how the availability of multimodality ion beam facilities can be used to design a tumor-specific optimal treatment.

Gene expression profiling of breast cancer cell lines treated with proton and electron radiations

OBJECTIVE

Technological advances in radiation therapy are evolving with the use of hadrons, such as protons, indicated for tumors where conventional radiotherapy does not give significant advantages or for tumors located in sensitive regions, which need the maximum of dose-saving of the surrounding healthy tissues. The genomic response to conventional and non-conventional linear energy transfer exposure is a poor investigated topic and became an issue of radiobiological interest. The aim of this work was to analyze and compare molecular responses in term of gene expression profiles, induced by electron and proton irradiation in breast cancer cell lines.

METHODS

We studied the gene expression profiling differences by cDNA microarray activated in response to electron and proton irradiation with different linear energy transfer values, among three breast cell lines (the tumorigenic MCF7 and MDA-MB-231 and the non-tumorigenic MCF10A), exposed to the same sublethal dose of 9 Gy.

RESULTS

Gene expression profiling pathway analyses showed the activation of different signaling and molecular networks in a cell line and radiation type-dependent manner. MCF10A and MDA-MB-231 cell lines were found to induce factors and pathways involved in the immunological process control.

CONCLUSION

Here, we describe in a detailed way the gene expression profiling and pathways activated after electron and proton irradiation in breast cancer cells. Summarizing, although specific pathways are activated in a radiation type-dependent manner, each cell line activates overall similar molecular networks in response to both these two types of ionizing radiation. Advances in knowledge: In the era of personalized medicine and breast cancer target-directed intervention, we trust that this study could drive radiation therapy towards personalized treatments, evaluating possible combined treatments, based on the molecular characterization.

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.

Radiosensitization Effect of Talazoparib, a Parp Inhibitor, on Glioblastoma Stem Cells Exposed to Low and High Linear Energy Transfer Radiation

Despite continuous improvements in treatment of glioblastoma, tumor recurrence and therapy resistance still occur in a high proportion of patients. One underlying reason for this radioresistance might be the presence of glioblastoma cancer stem cells (GSCs), which feature high DNA repair capability. PARP protein plays an important cellular role by detecting the presence of damaged DNA and then activating signaling pathways that promote appropriate cellular responses. Thus, PARP inhibitors (PARPi) have recently emerged as potential radiosensitizing agents. In this study, we investigated the preclinical efficacy of talazoparib, a new PARPi, in association with low and high linear energy transfer (LET) irradiation in two GSC cell lines. Reduction of GSC fraction, impact on cell proliferation, and cell cycle arrest were evaluated for each condition. All combinations were compared with a reference schedule: photonic irradiation combined with temozolomide. The use of PARPi combined with photon beam and even more carbon beam irradiation drastically reduced the GSC frequency of GBM cell lines in vitro. Furthermore, talazoparib combined with irradiation induced a marked and prolonged G2/M block, and decreased proliferation. These results show that talazoparib is a new candidate that effects radiosensitization in radioresistant GSCs, and its combination with high LET irradiation, is promising.

Correction of stopping power and LET quenching for radiophotoluminescent glass dosimetry in a therapeutic proton beam

To measure the absorbed dose to water D _w in proton beams using a radiophotoluminescent glass dosimeter (RGD), a method with the correction for the change of the mass stopping power ratio (SPR) and the linear energy transfer (LET) dependence of radiophotoluminescent efficiency [Formula: see text] is proposed. The calibration coefficient in terms of D _w for RGDs (GD-302M, Asahi Techno Glass) was obtained using a ⁶⁰Co γ-ray. The SPR of water to the RGD was calculated by Monte Carlo simulation, and [Formula: see text] was investigated experimentally using a 70 MeV proton beam. For clinical usage, the residual range R _res was used as a quality index to determine the correction factor for the beam quality [Formula: see text] and the LET quenching effect of the RGD [Formula: see text]. The proposed method was evaluated by measuring D _w at different depths in a 200 MeV proton beam. For both non-modulated and modulated proton beams, [Formula: see text] decreases rapidly where R _res is less than 4 cm. The difference in [Formula: see text] between a non-modulated and a modulated proton beam is less than 0.5% for the R _res range from 0 cm to 22 cm. [Formula: see text] decreases rapidly at a LET range from 1 to 2 keV µm^-1. In the evaluation experiments, D _w using RGDs, [Formula: see text] showed good agreement with that obtained using an ionization chamber and the relative difference was within 3% where R _res was larger than 1 cm. The uncertainty budget for [Formula: see text] in a proton beam was estimated to investigate the potential of RGD postal dosimetry in proton therapy. These results demonstrate the feasibility of RGD dosimetry in a therapeutic proton beam and the general versatility of the proposed method. In conclusion, the proposed methodology for RGDs in proton dosimetry is applicable where R _res  >  1 cm and the RGD is feasible as a postal audit dosimeter for proton therapy.

Robust intensity-modulated proton therapy to reduce high linear energy transfer in organs at risk

PURPOSE

We propose a robust treatment planning model that simultaneously considers proton range and patient setup uncertainties and reduces high linear energy transfer (LET) exposure in organs at risk (OARs) to minimize the relative biological effectiveness (RBE) dose in OARs for intensity-modulated proton therapy (IMPT). Our method could potentially reduce the unwanted damage to OARs.

METHODS

We retrospectively generated plans for 10 patients including two prostate, four head and neck, and four lung cancer patients. The "worst-case robust optimization" model was applied. One additional term as a "biological surrogate (BS)" of OARs due to the high LET-related biological effects was added in the objective function. The biological surrogate was defined as the sum of the physical dose and extra biological effects caused by the dose-averaged LET. We generated nine uncertainty scenarios that considered proton range and patient setup uncertainty. Corresponding to each uncertainty scenario, LET was obtained by a fast LET calculation method developed in-house and based on Monte Carlo simulations. In each optimization iteration, the model used the worst-case BS among all scenarios and then penalized overly high BS to organs. The model was solved by an efficient algorithm (limited-memory Broyden-Fletcher-Goldfarb-Shanno) in a parallel computing environment. Our new model was benchmarked with the conventional robust planning model without considering BS. Dose-volume histograms (DVHs) of the dose assuming a fixed RBE of 1.1 and BS for tumor and organs under nominal and uncertainty scenarios were compared to assess the plan quality between the two methods.

RESULTS

For the 10 cases, our model outperformed the conventional robust model in avoidance of high LET in OARs. At the same time, our method could achieve dose distributions and plan robustness of tumors assuming a fixed RBE of 1.1 almost the same as those of the conventional robust model.

CONCLUSIONS

Explicitly considering LET in IMPT robust treatment planning can reduce the high LET to OARs and minimize the possible toxicity of high RBE dose to OARs without sacrificing plan quality. We believe this will allow one to design and deliver safer proton therapy.

Clinical and Research Activities at the CATANA Facility of INFN-LNS: From the Conventional Hadrontherapy to the Laser-Driven Approach

The CATANA proton therapy center was the first Italian clinical facility making use of energetic (62 MeV) proton beams for the radioactive treatment of solid tumors. Since the date of the first patient treatment in 2002, 294 patients have been successful treated whose majority was affected by choroidal and iris melanomas. In this paper, we report on the current clinical and physical status of the CATANA facility describing the last dosimetric studies and reporting on the last patient follow-up results. The last part of the paper is dedicated to the description of the INFN-LNS ongoing activities on the realization of a beamline for the transport of laser-accelerated ion beams for future applications. The ELIMED (ELI-Beamlines MEDical and multidisciplinary applications) project is introduced and the main scientific aspects will be described.

Response of synthetic diamond detectors in proton, carbon, and oxygen ion beams

PURPOSE

In this work, the LET-dependence of the response of synthetic diamond detectors is investigated in different particle beams.

METHOD

Measurements were performed in three nonmodulated particle beams (proton, carbon, and oxygen). The response of five synthetic diamond detectors was compared to the response of a Markus or an Advanced Markus ionization chamber. The synthetic diamond detectors were used with their axis parallel to the beam axis and without any bias voltage. A high bias voltage was applied to the ionization chambers, to minimize ion recombination, for which no correction is applied (+300 V and +400 V were applied to the Markus and Advanced Markus ionization chambers respectively).

RESULTS

The ratio between the normalized response of the synthetic diamond detectors and the normalized response of the ionization chamber shows an under-response of the synthetic diamond detectors in carbon and oxygen ion beams. No under-response of the synthetic diamond detectors is observed in protons. For each beam, combining results obtained for the five synthetic diamond detectors and considering the uncertainties, a linear fit of the ratio between the normalized response of the synthetic diamond detectors and the normalized response of the ionization chamber is determined. The response of the synthetic diamond detectors can be described as a function of LET as (-6.22E-4 ± 3.17E-3) • LET + (0.99 ± 0.01) in proton beam, (-2.51E-4 ± 1.18E-4) • LET + (1.01 ± 0.01) in carbon ion beam and (-2.77E-4 ± 0.56E-4) • LET + (1.03 ± 0.01) in oxygen ion beam. Combining results obtained in carbon and oxygen ion beams, a LET dependence of about 0.026% (±0.013%) per keV/μm is estimated.

CONCLUSIONS

Due to the high LET value, a LET dependence of the response of the synthetic diamond detector was observed in the case of carbon and oxygen beams. The effect was found to be negligible in proton beams, due to the low LET value. The under-response of the synthetic diamond detector may result from the recombination of electron/hole in the thin synthetic diamond layer, due to the high LET-values. More investigations are required to confirm this assumption.

Linear energy transfer distributions in the brainstem depending on tumour location in intensity-modulated proton therapy of paediatric cancer

BACKGROUND

For tumours near organs at risk, there is concern about unintended increase in biological dose from elevated linear energy transfer (LET) at the distal end of treatment fields. The objective of this study was therefore to investigate how different paediatric posterior fossa tumour locations impact LET and biological dose to the brainstem during intensity-modulated proton therapy (IMPT).

MATERIAL AND METHODS

Multiple IMPT plans were generated for four different simulated tumour locations relative to the brainstem for a five-year-old male patient. A prescribed dose of 59.4 Gy(RBE) was applied to the planning target volumes (PTVs). Plans with two lateral and one posterior non-coplanar fields were created, along with plans with modified field arrangements. The dose-averaged LET (LET_d) and the physical dose × RBE_LET (D × RBE_LET), where RBE_LET=1+c × LET_d, were calculated using the FLUKA Monte Carlo code. A scaling parameter c was applied to make the RBE_LET represent variations in the biological effect due to LET.

RESULTS

High LET_d values surrounded parts of the PTV and encompassed portions of the brainstem. Mean LET_d values in the brainstem were 3.2-6.6 keV/μm. The highest absolute brainstem LET_d values were seen with the tumour located most distant from the brainstem, whereas lower and more homogeneous LET_d values were seen when the tumour invaded the brainstem. In contrast, the highest mean D × RBE_LET values were found in the latter case (54.0 Gy(RBE)), while the case with largest distance between tumour and brainstem had a mean D × RBE_LET of 1.8 Gy(RBE).

CONCLUSIONS

Using IMPT to treat posterior fossa tumours may result in high LET_d values within the brainstem, particularly if the tumour volume is separated from the brainstem. However, the D × RBE_LET was greater for tumours that approached or invaded the brainstem. Changing field angles showed a reduction of LET_d and D × RBE_LET in the brainstem.

Comparison of human lung cancer cell radiosensitivity after irradiations with therapeutic protons and carbon ions

The aim of this study was to investigate effects of irradiations with the therapeutic proton and carbon ion beams in two non-small cell lung cancers, CRL5876 adenocarcinoma and HTB177 large cell lung carcinoma. The DNA damage response dynamics, cell cycle regulation, and cell death pathway activation were followed. Viability of both cell lines was lower after carbon ions compared to the therapeutic proton irradiations. HTB177 cells showed higher recovery than CRL5876 cells seven days following the treatments, but the survival rates of both cell lines were lower after exposure to carbon ions with respect to therapeutic protons. When analyzing cell cycle distribution of both CRL5876 and HTB177 cells, it was noticed that therapeutic protons predominantly induced G1 arrest, while the cells after carbon ions were arrested in G2/M phase. The results illustrated that differences in the levels of phosphorylated H2AX, a double-strand break marker, exist after therapeutic proton and carbon ion irradiations. We also observed dose- and time-dependent increase in the p53 and p21 levels after applied irradiations. Carbon ions caused larger increase in the quantity of p53 and p21 compared to therapeutic protons. These results suggested that various repair mechanisms were induced in the treated cells. Considering the fact that we have not observed any distinct change in the Bax/Bcl-2 ratio following irradiations, it seemed that different types of cell death were involved in the response to the two types of irradiations that were applied.

A new approach to integrate GPU-based Monte Carlo simulation into inverse treatment plan optimization for proton therapy

Monte Carlo (MC)-based spot dose calculation is highly desired for inverse treatment planning in proton therapy because of its accuracy. Recent studies on biological optimization have also indicated the use of MC methods to compute relevant quantities of interest, e.g. linear energy transfer. Although GPU-based MC engines have been developed to address inverse optimization problems, their efficiency still needs to be improved. Also, the use of a large number of GPUs in MC calculation is not favorable for clinical applications. The previously proposed adaptive particle sampling (APS) method can improve the efficiency of MC-based inverse optimization by using the computationally expensive MC simulation more effectively. This method is more efficient than the conventional approach that performs spot dose calculation and optimization in two sequential steps. In this paper, we propose a computational library to perform MC-based spot dose calculation on GPU with the APS scheme. The implemented APS method performs a non-uniform sampling of the particles from pencil beam spots during the optimization process, favoring those from the high intensity spots. The library also conducts two computationally intensive matrix-vector operations frequently used when solving an optimization problem. This library design allows a streamlined integration of the MC-based spot dose calculation into an existing proton therapy inverse planning process. We tested the developed library in a typical inverse optimization system with four patient cases. The library achieved the targeted functions by supporting inverse planning in various proton therapy schemes, e.g. single field uniform dose, 3D intensity modulated proton therapy, and distal edge tracking. The efficiency was 41.6  ±  15.3% higher than the use of a GPU-based MC package in a conventional calculation scheme. The total computation time ranged between 2 and 50 min on a single GPU card depending on the problem size.

Nuclear physics in particle therapy: a review

Charged particle therapy has been largely driven and influenced by nuclear physics. The increase in energy deposition density along the ion path in the body allows reducing the dose to normal tissues during radiotherapy compared to photons. Clinical results of particle therapy support the physical rationale for this treatment, but the method remains controversial because of the high cost and of the lack of comparative clinical trials proving the benefit compared to x-rays. Research in applied nuclear physics, including nuclear interactions, dosimetry, image guidance, range verification, novel accelerators and beam delivery technologies, can significantly improve the clinical outcome in particle therapy. Measurements of fragmentation cross-sections, including those for the production of positron-emitting fragments, and attenuation curves are needed for tuning Monte Carlo codes, whose use in clinical environments is rapidly increasing thanks to fast calculation methods. Existing cross sections and codes are indeed not very accurate in the energy and target regions of interest for particle therapy. These measurements are especially urgent for new ions to be used in therapy, such as helium. Furthermore, nuclear physics hardware developments are frequently finding applications in ion therapy due to similar requirements concerning sensors and real-time data processing. In this review we will briefly describe the physics bases, and concentrate on the open issues.

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the geant4 Monte Carlo code

PURPOSE

The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the geant 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from geant 4 simulations to correlate with biological effectiveness of therapeutic protons.

METHODS

The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LETt and dose-averaged LET, LETd) using geant 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LETt and LETd of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose.

RESULTS

Simulation results showed that the step limit effect was small for LETt but significant for LETd. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in geant 4 can result in incorrect LETd calculation results in the dose plateau region for small step limits. The erroneous LETd results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in geant 4. The incorrect LETd values lead to substantial differences in the calculated RBE.

CONCLUSIONS

When the geant 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LETt in the dose plateau region and LETd around the Bragg peak. For a large step limit, i.e., 500 μm, LETd is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LETd and LETt becomes positive.

Under-response of a PTW-60019 microDiamond detector in the Bragg peak of a 62 MeV/n carbon ion beam

To investigate the linear energy transfer (LET) dependence of the response of a PTW-60019 Freiburg microDiamond detector, its response was compared to the response of a plane-parallel Markus chamber in a 62 MeV/n mono-energetic carbon ion beam. Results obtained with two different experimental setups are in agreement. As recommended by IAEA TRS-398, the response of the Markus chamber was corrected for temperature, pressure, polarity effects and ion recombination. No correction was applied to the response of the microDiamond detector. The ratio of the response of the Markus chamber to the response of the microDiamond is close to unity in the plateau region. In the Bragg peak region, a significant increase of the ratio is observed, which increases to 1.2 in the distal edge region. Results indicate a correlation between the under-response of the microDiamond detector and high LET values. The combined relative standard uncertainty of the results is estimated to be 2.38% in the plateau region and 12% in the distal edge region. These values are dominated by the uncertainty of alignment in the non-uniform beam and the uncertainty of range determination.