Performance of LGAD strip detectors for particle counting of therapeutic proton beams

Objective. The performance of silicon detectors with moderate internal gain, named low-gain avalanche diodes (LGADs), was studied to investigate their capability to discriminate and count single beam particles at high fluxes, in view of future applications for beam characterization and on-line beam monitoring in proton therapy.Approach. Dedicated LGAD detectors with an active thickness of 55μm and segmented in 2 mm2strips were characterized at two Italian proton-therapy facilities, CNAO in Pavia and the Proton Therapy Center of Trento, with proton beams provided by a synchrotron and a cyclotron, respectively. Signals from single beam particles were discriminated against a threshold and counted. The number of proton pulses for fixed energies and different particle fluxes was compared with the charge collected by a compact ionization chamber, to infer the input particle rates.Main results. The counting inefficiency due to the overlap of nearby signals was less than 1% up to particle rates in one strip of 1 MHz, corresponding to a mean fluence rate on the strip of about 5 × 107p/(cm2·s). Count-loss correction algorithms based on the logic combination of signals from two neighboring strips allow to extend the maximum counting rate by one order of magnitude. The same algorithms give additional information on the fine time structure of the beam.Significance. The direct counting of the number of beam protons with segmented silicon detectors allows to overcome some limitations of gas detectors typically employed for beam characterization and beam monitoring in particle therapy, providing faster response times, higher sensitivity, and independence of the counts from the particle energy.

Radiation measurements in the International Space Station, Columbus module, in 2020-2022 with the LIDAL detector

The Light Ion Detector for ALTEA (LIDAL) is a new instrument designed to measure flux, energy spectra and Time of Flight of ions in a space habitat. It was installed in the International Space Station (Columbus) on January 19, 2020 and it is still operating. This paper presents the results of LIDAL measurements in the first 17 months of operation (01/2020-05/2022). Particle flux, dose rate, Time of Flight and spectra are presented and studied in the three ISS orthogonal directions and in the different geomagnetic regions (high latitude, low latitude, and South Atlantic Anomaly, SAA). The results are consistent with previous measurements. Dose rates range between 1.8 nGy/s and 2.4 nGy/s, flux between 0.21 particles/(sr cm2 s) and 0.32 particles/(sr cm2 s) as measured across time and directions during the full orbit. These data offer insights concerning the radiation measurements in the ISS and demonstrate the capabilities of LIDAL as a unique tool for the measurement of space radiation in space habitats, also providing novel information relevant to assess radiation risks for astronauts.

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.

Calibration method and performance of a time-of-flight detector to measure absolute beam energy in proton therapy

BACKGROUND

The beam energy is one of the most significant parameters in particle therapy since it is directly correlated to the particles' penetration depth inside the patient. Nowadays, the range accuracy is guaranteed by offline routine quality control checks mainly performed with water phantoms, 2D detectors with PMMA wedges, or multi-layer ionization chambers. The latter feature low sensitivity, slow collection time, and response dependent on external parameters, which represent limiting factors for the quality controls of beams delivered with fast energy switching modalities, as foreseen in future treatments. In this context, a device based on solid-state detectors technology, able to perform a direct and absolute beam energy measurement, is proposed as a viable alternative for quality assurance measurements and beam commissioning, paving the way for online range monitoring and treatment verification.

PURPOSE

This work follows the proof of concept of an energy monitoring system for clinical proton beams, based on Ultra Fast Silicon Detectors (featuring tenths of ps time resolution in 50 μm active thickness, and single particle detection capability) and time-of-flight techniques. An upgrade of such a system is presented here, together with the description of a dedicated self-calibration method, proving that this second prototype is able to assess the mean particles energy of a monoenergetic beam without any constraint on the beam temporal structure, neither any a priori knowledge of the beam energy for the calibration of the system.

METHODS

A new detector geometry, consisting of sensors segmented in strips, has been designed and implemented in order to enhance the statistics of coincident protons, thus improving the accuracy of the measured time differences. The prototype was tested on the cyclotron proton beam of the Trento Protontherapy Center (TPC). In addition, a dedicated self-calibration method, exploiting the measurement of monoenergetic beams crossing the two telescope sensors for different flight distances, was introduced to remove the systematic uncertainties independently from any external reference.

RESULTS

The novel calibration strategy was applied to the experimental data collected at TPC (Trento) and CNAO (Pavia). Deviations between measured and reference beam energies in the order of a few hundreds of keV with a maximum uncertainty of 0.5 MeV were found, in compliance with the clinically required water range accuracy of 1 mm.

CONCLUSIONS

The presented version of the telescope system, minimally perturbative of the beam, relies on a few seconds of acquisition time to achieve the required clinical accuracy and therefore represents a feasible solution for beam commission, quality assurance checks, and online beam energy monitoring.

The INFN proton computed tomography system for relative stopping power measurements: calibration and verification

Objective. This paper describes the procedure to calibrate the three-dimensional (3D) proton stopping power relative to water (SPR) maps measured by the proton computed tomography (pCT) apparatus of the Istituto Nazionale di Fisica Nucleare (INFN, Italy). Measurements performed on water phantoms are used to validate the method. The calibration allowed for achieving measurement accuracy and reproducibility to levels below 1%.Approach. The INFN pCT system is made of a silicon tracker for proton trajectory determination followed by a YAG:Ce calorimeter for energy measurement. To perform the calibration, the apparatus has been exposed to protons of energies ranging from 83 to 210 MeV. Using the tracker, a position-dependent calibration has been implemented to keep the energy response uniform across the calorimeter. Moreover, correction algorithms have been developed to reconstruct the proton energy when this is shared in more than one crystal and to consider the energy loss in the non-uniform apparatus material. To verify the calibration and its reproducibility, water phantoms have been imaged with the pCT system during two data-taking sessions.Main results. The energy resolution of the pCT calorimeter resulted to beσEE≅0.9%at 196.5 MeV. The average values of the water SPR in fiducial volumes of the control phantoms have been calculated to be 0.995±0.002. The image non-uniformities were below 1%. No appreciable variation of the SPR and uniformity values between the two data-taking sessions could be identified.Significance. This work demonstrates the accuracy and reproducibility of the calibration of the INFN pCT system at a level below 1%. Moreover, the uniformity of the energy response keeps the image artifacts at a low level even in the presence of calorimeter segmentation and tracker material non-uniformities. The implemented calibration technique allows the INFN-pCT system to face applications where the precision of the SPR 3D maps is of paramount importance.

Characterization of the INFN proton CT scanner for cross-calibration of x-ray CT

OBJECTIVE

The goal of this study was to assess the imaging performances of the pCT system developed in the framework of INFN-funded (Italian National Institute of Nuclear Physics) research projects. The spatial resolution, noise power spectrum and RSP accuracy has been investigated, as a preliminary step to implement a new cross-calibration method for x-ray CT (xCT). 
Approach: The INFN pCT apparatus, made of four planes of silicon micro-strip detectors and a YAG:Ce scintillating calorimeter, reconstructs 3D RSP maps by a filtered-back projection algorithm. The imaging performances (i.e. spatial resolution, noise power spectrum and RSP accuracy) of the pCT system were assessed on a custom-made phantom, made of plastic materials with different densities ([0.66, 2.18] g/cm3). For comparison, the same phantom was acquired with a clinical xCT system.
Main results: The spatial resolution analysis revealed the non-linearity of the imaging system, showing different imaging responses in air or water phantom background. Applying the Hann filter in the pCT reconstruction, it was possible to investigate the imaging potential of the system. Matching the spatial resolution value of the xCT (0.54 lp/mm) and acquiring both with the same dose level (11.6 mGy), the pCT appeared to be less noisy than xCT, with an RSP standard deviation of 0.0063. Concerning the RSP accuracy, the measured Mean Absolute Percentage Errors were (0.23+-0.09)% in air and (0.21+-0.07)% in water.
Significance: The obtained performances confirm that the INFN pCT system provides a very accurate RSP estimation, appearing to be a feasible clinical tool for verification and correction of xCT calibration in proton treatment planning.

Does variable RBE affect toxicity risks for mediastinal lymphoma patients? NTCP-based evaluation after proton therapy treatment

INTRODUCTION

Mediastinal lymphoma (ML) is a solid malignancy affecting young patients. Modern combined treatments allow obtaining good survival probability, together with a long life expectancy, and therefore with the need to minimize treatment-related toxicities. We quantified the expected toxicity risk for different organs and endpoints in ML patients treated with intensity-modulated proton therapy (IMPT) at our centre, accounting also for uncertainties related to variable RBE.

METHODS

Treatment plans for ten ML patients were recalculated with a TOPAS-based Monte Carlo code, thus retrieving information on LET and allowing the estimation of variable RBE. Published NTCP models were adopted to calculate the toxicity risk for hypothyroidism, heart valve defects, coronary heart disease and lung fibrosis. NTCP was calculated assuming both constant (i.e. 1.1) and variable RBE. The uncertainty associated with individual radiosensitivity was estimated by random sampling α/β values before RBE evaluation.

RESULTS

Variable RBE had a minor impact on hypothyroidism risk for 7 patients, while it led to significant increase for the remaining three (+24% risk maximum increase). Lung fibrosis was slightly affected by variable RBE, with a maximum increase of ≅ 1%. This was similar for heart valve dysfunction, with the exception of one patient showing an about 10% risk increase, which could be explained by means of large heart volume and D₁ increase.

DISCUSSION

The use of NTCP models allows for identifying those patients associated with a higher toxicity risk. For those patients, it might be worth including variable RBE in plan evaluation.

Investigation of In-Field and Out-of-Field Radiation Quality With Microdosimetry and Its Impact on Relative Biological Effectiveness in Proton Therapy

PURPOSE

Using microdosimetry, this study investigated the relative biological effectiveness (RBE) and quality factor (Q¯) variations in field and out of field as a function of radiation quality for clinical protons.

METHODS AND MATERIALS

A water phantom with a spread-out Bragg peak (SOBP) was irradiated to acquire microdosimetric spectra at several distal and lateral depths with a tissue equivalent proportional counter. The measurements were used as inputs to microdosimetric kinetic and Loncol models to determine the RBE spatial distribution and compare it with predictions from the dose-averaged linear energy transfer-based McNamara model. Q¯ values and biological and dose equivalent values were also calculated.

RESULTS

The data demonstrated that radiation quality changed more rapidly with depth than lateral distance from the SOBP. In beam, y_D ranged from approximately 4 keV/μm at the entrance to 8 keV/μm at the SOBP far end, reaching approximately 15 keV/μm at the penumbra. Out of field, the overall highest value of 23 ± 2 keV/μm was observed at the beam-edge penumbra. Radiation quality changes caused RBE deviations from the clinical value of 1.1, whose extent depends on the approach used for assessing radiation quality as well as on the radiobiological model. For RBE₁₀, microdosimetry-based models appeared to better reproduce the radiobiological data than the dose-averaged linear energy transfer model. Out of field, both the RBE and Q¯ values appeared to have limitations in describing the radiation biological effectiveness. This research also presents a first comprehensive benchmark of TOPAS code against in-field and out-of-field microdosimetric spectra of therapeutic protons.

CONCLUSIONS

Further investigation will be necessary to evaluate the quantitative effects of RBE variations on treatment planning and assess the clinical consequences in terms of both tumor control and normal-tissue toxicity. The achievement of this goal calls for accurate radiobiological data to validate the RBE models.

LIDAL, a Time-of-Flight Radiation Detector for the International Space Station: Description and Ground Calibration

LIDAL (Light Ion Detector for ALTEA, Anomalous Long-Term Effects on Astronauts) is a radiation detector designed to measure the flux, the energy spectra and, for the first time, the time-of-flight of ions in a space habitat. It features a combination of striped silicon sensors for the measurement of deposited energy (using the ALTEA device, which operated from 2006 to 2012 in the International Space Station) and fast scintillators for the time-of-flight measurement. LIDAL was tested and calibrated using the proton beam line at TIFPA (Trento Institute for Fundamental Physics Application) and the carbon beam line at CNAO (National Center for Oncology Hadron-therapy) in 2019. The performance of the time-of-flight system featured a time resolution (sigma) less than 100 ps. Here, we describe the detector and the results of these tests, providing ground calibration curves along with the methodology established for processing the detector's data. LIDAL was uploaded in the International Space Station in November 2019 and it has been operative in the Columbus module since January 2020.

Diagnosis of Langerhans cell histiocytosis on cytological examination of cerebrospinal fluid: Report of the first case

Langerhans cell histiocytosis (LCH) is a disease of unknown etiology characterized by a proliferation of histiocytic cells resembling dendritic Langerhans cells. LCH can be unifocal or multifocal, with one- or many-organ involvement. The serous fluids are rarely involved. Cytological diagnosis of LCH is possible and relies on recognition of the typical cytomorphological features and subsequent immunocytochemical confirmation. Given the possibility of multisystem involvement, after diagnosing LCH it is necessary to carry out staging exams such as a bone survey, abdominal ultrasound, complete blood count, screening for diabetes insipidus and pulmonary function tests. We present the first case of LCH where the diagnosis was reached on cytological material from the cerebrospinal fluid. To the best of our knowledge, this is the first such case reported in the international literature to date. The morphological and immunocytochemical characteristics of our case are described, and the relevant literature is reviewed.

Modelling the HPRT-gene mutation induction of particle beams: systematic in vitro data collection, analysis and microdosimetric kinetic model implementation

Objective. Since the early years, particle therapy treatments have been associated with concerns for late toxicities, especially secondary cancer risk (SCR). Nowadays, this concern is related to patients for whom long-term survival is expected (e.g. breast cancer, lymphoma, paediatrics). In the aim to contribute to this research, we present a dedicated statistical and modelling analysis aiming at improving our understanding of the RBE for mutation induction ( RBE M ˜ ) for different particle species. Approach. We built a new database based on a systematic collection of RBE data for mutation assays of the gene encoding for the purine salvage enzyme hypoxanthine-guanine phosphoribosyltransferase from literature (105 entries, distributed among 3 cell lines and 16 particle species). The data were employed to perform statistical and modelling analysis. For the latter, we adapted the microdosimetric kinetic model (MKM) to describe the mutagenesis in analogy to lethal lesion induction. Main results. Correlation analysis between RBE for survival (RBES) and RBE M ˜ reveals significant correlation between these two quantities (ρ = 0.86, p < 0.05). The correlation gets stronger when looking at subsets of data based on cell line and particle species. We also show that the MKM can be successfully employed to describe RBE M ˜ , obtaining comparably good agreement with the experimental data. Remarkably, to improve the agreement with experimental data the MKM requires, consistently in all the analysed cases, a reduced domain size for the description of mutation induction compared to that adopted for survival. Significance. We were able to show that RBES and RBE M ˜ are strongly related quantities. We also showed for the first time that the MKM could be successfully applied to the description of mutation induction, representing an endpoint different from the more traditional cell killing. In analogy to the RBES, RBE M ˜ can be implemented into treatment planning system evaluations.

Quantification of biological range uncertainties in patients treated at the Krakow proton therapy centre

Background

Variable relative biological effectiveness (vRBE) in proton therapy might significantly modify the prediction of RBE-weighted dose delivered to a patient during proton therapy. In this study we will present a method to quantify the biological range extension of the proton beam, which results from the application of vRBE approach in RBE-weighted dose calculation.

Methods and materials

The treatment plans of 95 patients (brain and skull base patients) were used for RBE-weighted dose calculation with constant and the McNamara RBE model. For this purpose the Monte Carlo tool FRED was used. The RBE-weighted dose distributions were analysed using indices from dose-volume histograms. We used the volumes receiving at least 95% of the prescribed dose (V95) to estimate the biological range extension resulting from vRBE approach.

Results

The vRBE model shows higher median value of relative deposited dose and D95 in the planning target volume by around 1% for brain patients and 4% for skull base patients. The maximum doses in organs at risk calculated with vRBE was up to 14 Gy above dose limit. The mean biological range extension was greater than 0.4 cm.

Discussion

Our method of estimation of biological range extension is insensitive for dose inhomogeneities and can be easily used for different proton plans with intensity-modulated proton therapy (IMPT) optimization. Using volumes instead of dose profiles, which is the common method, is more universal. However it was tested only for IMPT plans on fields arranged around the tumor area.

Conclusions

Adopting a vRBE model results in an increase in dose and an extension of the beam range, which is especially disadvantageous in cancers close to organs at risk. Our results support the need to re-optimization of proton treatment plans when considering vRBE.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-022-02022-5.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach

Simple Summary

Proton therapy is now an established external radiotherapy modality for cancer treatment. Clinical routine currently neglects the radiobiological impact of nuclear target fragments even if experimental evidences show a significant enhancement in cell-killing effect due to secondary particles. This paper quantifies the contribution of proton target fragments of different charge in different irradiation scenarios and compares the computationally predicted corrections to the overall biological dose with experimental data.

Abstract

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments’ contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/β values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments’ contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

Study of relationship between dose, LET and the risk of brain necrosis after proton therapy for skull base tumors

PURPOSE

We investigated the relationship between RBE-weighted dose (DRBE) calculated with constant (cRBE) and variable RBE (vRBE), dose-averaged linear energy transfer (LETd) and the risk of radiographic changes in skull base patients treated with protons.

METHODS

Clinical treatment plans of 45 patients were recalculated with Monte Carlo tool FRED. Radiographic changes (i.e. edema and/or necrosis) were identified by MRI. Dosimetric parameters for cRBE and vRBE were computed. Biological margin extension and voxel-based analysis were employed looking for association of DRBE(vRBE) and LETd with brain edema and/or necrosis.

RESULTS

When using vRBE, Dmax in the brain was above the highest dose limits for 38% of patients, while such limit was never exceeded assuming cRBE. Similar values of Dmax were observed in necrotic regions, brain and temporal lobes. Most of the brain necrosis was in proximity to the PTV. The voxel-based analysis did not show evidence of an association with high LETd values.

CONCLUSIONS

When looking at standard dosimetric parameters, the higher dose associated with vRBE seems to be responsible for an enhanced risk of radiographic changes. However, as revealed by a voxel-based analysis, the large inter-patient variability hinders the identification of a clear effect for high LETd.

Technical Note: CT calibration for proton treatment planning by cross-calibration with proton CT data

PURPOSE

This study explores the possibility of a new method for x-ray computed tomography (CT) calibration by means of cross-calibration with proton CT (pCT) data. The proposed method aims at a more accurate conversion of CT Hounsfield Units (HU) into proton stopping power ratio (SPR) relative to water to be used in proton-therapy treatment planning.

METHODS

X-ray CT scan was acquired on a synthetic anthropomorphic phantom, composed of different tissue equivalent materials (TEMs). A pCT apparatus was instead adopted to obtain a reference three-dimensional distribution of the phantom's SPR values. After rigid registration, the x-ray CT was artificially blurred to the same resolution of pCT. Then a scatter plot showing voxel-by-voxel SPR values as a function of HU was employed to link the two measurements and thus obtaining a cross-calibrated x-ray CT calibration curve. The cross-calibration was tested at treatment planning system and then compared with a conventional calibration based on exactly the same TEMs constituting the anthropomorphic phantom.

RESULTS

Cross-calibration provided an accurate SPR mapping, better than by conventional TEMs calibration. The dose distribution of single beams optimized on the reference SPR map was recomputed on cross-calibrated CT, showing, with respect to conventional calibration, minor deviation at the dose fall-off (lower than 1%).

CONCLUSIONS

The presented data demonstrated that, by means of reference pCT data, a heterogeneous phantom can be used for CT calibration, paving the way to the use of biological samples, with their accurate description of patients' tissues. This overcomes the limitations of conventional CT calibration requiring homogenous samples, only available by synthetic TEMs, which fail in accurately mimicking the properties of biological tissues. Once a heterogeneous biological sample is provided with its corresponding reference SPR maps, a cross-calibration procedure could be adopted by other PT centers, even when not equipped with a pCT system.

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.

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.

Proton or photon radiosurgery for cardiac ablation of ventricular tachycardia? Breath and ECG gated robust optimization

PURPOSE

Ventricular tachycardia (VT) is a life-threatening heart disorder. The aim of this preliminary study is to assess the feasibility of stereotactic body radiation therapy (SBRT) photon and proton therapy (PT) plans for the treatment of VT, adopting robust optimization technique for both irradiation techniques.

METHODS

ECG gated CT images (in breath hold) were acquired for one patient. Conventional planning target volume (PTV) and robust optimized plans (25GyE in single fraction) were simulated for both photon (IMRT, 5 and 9 beams) and proton (SFO, 2 beams) plans. Robust optimized plans were obtained both for protons and photons considering in the optimization setup errors (5 mm in the three orthogonal directions), range (±3.5%) and the clinical target volume (CTV) motion due to heartbeat and breath-hold variability.

RESULTS

The photon robust optimization method, compared to PTV-based plans, showed a reduction in the average dose to the heart by about 25%; robust optimization allowed also reducing the mean dose to the left lung from 3.4. to 2.8 Gy for 9-beams configuration and from 4.1 to 2.9 Gy for 5-beams configuration. Robust optimization with protons, allowed further reducing the OAR doses: average dose to the heart and to the left lung decreased from 7.3 Gy to 5.2 GyE and from 2.9 Gy to 2.2 GyE, respectively.

CONCLUSIONS

Our study demonstrates the importance of the optimization technique adopted in the treatment planning system for VT treatment. It has been shown that robust optimization can significantly reduce the dose to healthy cardiac tissues and that PT further increases this gain.

Evaluation of proton beam radiation-induced skin injury in a murine model using a clinical SOBP

The purpose of this paper is to characterize the skin deterministic damage due to the effect of proton beam irradiation in mice occurred during a long-term observational experiment. This study was initially defined to evaluate the insurgence of myelopathy irradiating spinal cords with the distal part of a Spread-out Bragg peak (SOBP). To the best of our knowledge, no study has been conducted highlighting high grades of skin injury at the dose used in this paper. Nevertheless these effects occurred. In this regard, the experimental evidence of significant insurgence of skin injury induced by protons using a SOBP configuration will be shown. Skin damages were classified into six scores (from 0 to 5) according to the severity of the injuries and correlated to ED50 (i.e. the radiation dose at which 50% of animals show a specific score) at 40 days post-irradiation (d.p.i.). The effects of radiation on the overall animal wellbeing have been also monitored and the severity of radiation-induced skin injuries was observed and quantified up to 40 d.p.i.

Clinical implementation in proton therapy of multi-field optimization by a hybrid method combining conventional PTV with robust optimization

To implement a robust multi-field optimization (MFO) technique compatible with the application of a Monte Carlo (MC) algorithm and to evaluate its robustness. Nine patients (three brain, five head-and-neck, one spine) underwent proton treatment generated by a novel robust MFO technique. A hybrid (hMFO) approach was implemented, planning dose coverage on isotropic PTV compensating for setup errors, whereas range calibration uncertainties are incorporated into PTV robust optimization process. hMFO was compared with single-field optimization (SFO) and full robust multi-field optimization (fMFO), both on the nominal plan and the worst-case scenarios assessed by robustness analysis. The SFO and the fMFO plans were normalized to hMFO on CTV to obtain iso-D95 coverage, and then the organs at risk (OARs) doses were compared. On the same OARs, in the normalized nominal plans the potential impact of variable relative biological effectiveness (RBE) was investigated. hMFO reduces the number of scenarios computed for robust optimization (from twenty-one in fMFO to three), making it practicable with the application of a MC algorithm. After normalizing on D95 CTV coverage, nominal hMFO plans were superior compared to SFO in terms of OARs sparing (p   <  0.01), without significant differences compared to fMFO. The improvement in OAR sparing with hMFO with respect to SFO was preserved in worst-case scenarios (p   <  0.01), confirming that hMFO is as robust as SFO to physical uncertainties, with no significant differences when compared to the worst case scenarios obtained by fMFO. The dose increase on OARs due to variable RBE was comparable to the increase due to physical uncertainties (i.e. 4-5 Gy(RBE)), but without significant differences between these techniques. hMFO allows improving plan quality with respect to SFO, with no significant differences with fMFO and without affecting robustness to setup, range and RBE uncertainties, making clinically feasible the application of MC-based robust optimization.

Accurate proton treatment planning for pencil beam crossing titanium fixation implants

PURPOSE

To present a planning strategy for proton pencil-beam scanning when titanium implants need to be crossed by the beam.

METHODS

We addressed three issues: the implementation of a CT calibration curve to assign to titanium the correct stopping power; the effect of artefacts on CT images and their reduction by a dedicated algorithm; the differences in dose computation depending on the dose engine, pencil-beam vs Monte-Carlo algorithms. We performed measurement tests on a simple cylinder phantom and on a real implant. These phantoms were irradiated with three geometries (single spots, uniform mono-energetic layer and uniform box), measuring the exit dose either by radio-chromic film or multi-layer ionization chamber. The procedure was then applied on two patients treated for chordoma.

RESULTS

We had to set in the calibration curve a mass density equal to 4.37 g/cm³ to saturated Hounsfield Units, in order to have the correct stopping power assigned to titanium in TPS. CT artefact reduction algorithm allowed a better reconstruction of the shape and size of the implant. Monte-Carlo resulted accurate in computing the dose distribution whereas the pencil-beam algorithm failed due to sharp density interfaces between titanium and the surrounding material. Finally, the treatment plans obtained on two patients showed the impact of the dose engine algorithm, with 10-20% differences between pencil-beam and Monte-Carlo in small regions distally to the titanium screws.

CONCLUSION

The described combination of CT calibration, artefacts reduction and Monte-Carlo computation provides a reliable methodology to compute dose in patients with titanium implants.

Modelling the risk of radiation induced alopecia in brain tumor patients treated with scanned proton beams

PURPOSE

To develop normal tissue complication probability (NTCP) models for radiation-induced alopecia (RIA) in brain tumor patients treated with proton therapy (PT).

METHODS AND MATERIALS

We analyzed 116 brain tumor adult patients undergoing scanning beam PT (median dose 54 GyRBE; range 36-72) for CTCAE v.4 grade 2 (G2) acute (≤90 days), late (>90 days) and permanent (>12 months) RIA. The relative dose-surface histogram (DSH) of the scalp was extracted and used for Lyman-Kutcher-Burman (LKB) modelling. Moreover, DSH metrics (Sx: the surface receiving ≥ X Gy, D_2%: near maximum dose, D_mean: mean dose) and non-dosimetric variables were included in a multivariable logistic regression NTCP model. Model performances were evaluated by the cross-validated area under the receiver operator curve (ROC-AUC).

RESULTS

Acute, late and permanent G2-RIA was observed in 52%, 35% and 19% of the patients, respectively. The LKB models showed a weak dose-surface effect (0.09 ≤ n ≤ 0.19) with relative steepness 0.29 ≤ m ≤ 0.56, and increasing tolerance dose values when moving from acute and late (22 and 24 GyRBE) to permanent RIA (44 GyRBE). Multivariable modelling selected S_21Gy for acute and S_25Gy, for late G2-RIA as the most predictive DSH factors. Younger age was selected as risk factor for acute G2-RIA while surgery as risk factor for late G2-RIA. D_2% was the only variable selected for permanent G2-RIA. Both LKB and logistic models exhibited high predictive performances (ROC-AUCs range 0.86-0.90).

CONCLUSION

We derived NTCP models to predict G2-RIA after PT, providing a comprehensive modelling framework for acute, late and permanent occurrences that, once externally validated, could be exploited for individualized scalp sparing treatment planning strategies in brain tumor patients.

Potential skin morbidity reduction with intensity-modulated proton therapy for breast cancer with nodal involvement

Background: Different modern radiation therapy treatment solutions for breast cancer (BC) and regional nodal irradiation (RNI) have been proposed. In this study, we evaluate the potential reduction in radiation-induced skin morbidity obtained by intensity modulated proton therapy (IMPT) compared with intensity modulated photon therapy (IMXT) for left-side BC and RNI. Material and Methods: Using CT scans from 10 left-side BC patients, treatment plans were generated using IMXT and IMPT techniques. A dose of 50 Gy (or Gy [RBE] for IMPT) was prescribed to the target volume (involved breast, the internal mammary, supraclavicular, and infraclavicular nodes). Two single filed optimization IMPT (IMPT₁ and IMPT₂) plans were calculated without and with skin optimization. For each technique, skin dose-metrics were extracted and normal tissue complication probability (NTCP) models from the literature were employed to estimate the risk of radiation-induced skin morbidity. NTCPs for relevant organs-at-risk (OARs) were also considered for reference. The non-parametric Anova (Friedman matched-pairs signed-rank test) was used for comparative analyses. Results: IMPT improved target coverage and dose homogeneity even if the skin was included into optimization strategy (HI_IMPT2 = 0.11 vs. HI_IMXT = 0.22 and CI_IMPT2 = 0.96 vs. CI_IMXT = 0.82, p < .05). A significant relative skin risk reduction (RR = NTCP_IMPT/NTCP_IMXT) was obtained with IMPT₂ including the skin in the optimization with a RR reduction ranging from 0.3 to 0.9 depending on the analyzed skin toxicity endpoint/model. Both IMPT plans attained significant OARs dose sparing compared with IMXT. As expected, the heart and lung doses were significantly reduced using IMPT. Accordingly, IMPT always provided lower NTCP values. Conclusions: IMPT guarantees optimal target coverage, OARs sparing, and simultaneously minimizes the risk of skin morbidity. The applied model-based approach supports the potential clinical relevance of IMPT for left-side BC and RNI and might be relevant for the setup of cost-effectiveness evaluation strategies based on NTCP predictions, as well as for establishing patient selection criteria.

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.

Characterizing the Potency and Impact of Carbon Ion Therapy in a Primary Mouse Model of Soft Tissue Sarcoma

Carbon ion therapy (CIT) offers several potential advantages for treating cancers compared with X-ray and proton radiotherapy, including increased biological efficacy and more conformal dosimetry. However, CIT potency has not been characterized in primary tumor animal models. Here, we calculate the relative biological effectiveness (RBE) of carbon ions compared with X-rays in an autochthonous mouse model of soft tissue sarcoma. We used Cre/loxP technology to generate primary sarcomas in KrasLSL-G12D/+; p53fl/fl mice. Primary tumors were irradiated with a single fraction of carbon ions (10 Gy), X-rays (20 Gy, 25 Gy, or 30 Gy), or observed as controls. The RBE was calculated by determining the dose of X-rays that resulted in similar time to posttreatment tumor volume quintupling and exponential growth rate as 10 Gy carbon ions. The median tumor volume quintupling time and exponential growth rate of sarcomas treated with 10 Gy carbon ions and 30 Gy X-rays were similar: 27.3 and 28.1 days and 0.060 and 0.059 mm3/day, respectively. Tumors treated with lower doses of X-rays had faster regrowth. Thus, the RBE of carbon ions in this primary tumor model is 3. When isoeffective treatments of carbon ions and X-rays were compared, we observed significant differences in tumor growth kinetics, proliferative indices, and immune infiltrates. We found that carbon ions were three times as potent as X-rays in this aggressive tumor model and identified unanticipated differences in radiation response that may have clinical implications. Mol Cancer Ther; 17(4); 858-68. ©2018 AACR.

Model-based approach for quantitative estimates of skin, heart, and lung toxicity risk for left-side photon and proton irradiation after breast-conserving surgery

BACKGROUND

Proton beam therapy represents a promising modality for left-side breast cancer (BC) treatment, but concerns have been raised about skin toxicity and poor cosmesis. The aim of this study is to apply skin normal tissue complication probability (NTCP) model for intensity modulated proton therapy (IMPT) optimization in left-side BC.

MATERIAL AND METHODS

Ten left-side BC patients undergoing photon irradiation after breast-conserving surgery were randomly selected from our clinical database. Intensity modulated photon (IMRT) and IMPT plans were calculated with iso-tumor-coverage criteria and according to RTOG 1005 guidelines. Proton plans were computed with and without skin optimization. Published NTCP models were employed to estimate the risk of different toxicity endpoints for skin, lung, heart and its substructures.

RESULTS

Acute skin NTCP evaluation suggests a lower toxicity level with IMPT compared to IMRT when the skin is included in proton optimization strategy (0.1% versus 1.7%, p < 0.001). Dosimetric results show that, with the same level of tumor coverage, IMPT attains significant heart and lung dose sparing compared with IMRT. By NTCP model-based analysis, an overall reduction in the cardiopulmonary toxicity risk prediction can be observed for all IMPT compared to IMRT plans: the relative risk reduction from protons varies between 0.1 and 0.7 depending on the considered toxicity endpoint.

CONCLUSIONS

Our analysis suggests that IMPT might be safely applied without increasing the risk of severe acute radiation induced skin toxicity. The quantitative risk estimates also support the potential clinical benefits of IMPT for left-side BC irradiation due to lower risk of cardiac and pulmonary morbidity. The applied approach might be relevant on the long term for the setup of cost-effectiveness evaluation strategies based on NTCP predictions.

Hibernation for space travel: Impact on radioprotection

Hibernation is a state of reduced metabolic activity used by some animals to survive in harsh environmental conditions. The idea of exploiting hibernation for space exploration has been proposed many years ago, but in recent years it is becoming more realistic, thanks to the introduction of specific methods to induce hibernation-like conditions (synthetic torpor) in non-hibernating animals. In addition to the expected advantages in long-term exploratory-class missions in terms of resource consumptions, aging, and psychology, hibernation may provide protection from cosmic radiation damage to the crew. Data from over half century ago in animal models suggest indeed that radiation effects are reduced during hibernation. We will review the mechanisms of increased radioprotection in hibernation, and discuss possible impact on human space exploration.

New Ions for Therapy

Purpose

Charged particle therapy (CPT) is currently based on the use of protons or carbon ions for the treatment of deep-seated and/or radioresistant tumors, which are known to return poor prognosis in photon treatments. A renovated interest has recently been observed in the possibility of extending the spectrum of ions used in CPT. The potential and limitations of different particle species will be discussed in this work, with special regard to ¹H, ⁴He, ¹²C, and ¹⁶O, that is, those presently available in the most advanced particle therapy clinical centers.

Materials and Methods

Literature information has been screened, as well as additional analysis has been performed, aimed at the comparison of basic physical and biological properties of several ions. The research treatment planning system TRiP98 is also employed to compare the dose distribution resulting from exposure to the different ions in different configurations, including the irradiation of hypoxic targets.

Results

Particles of intermediate charge, such as helium and lithium, offer an increased physical selectivity compared with protons, while having reduced biological effectiveness compared with carbon. The latter aspect translates into a less sensitive biological optimization of CPT treatments, though still more effective than protons in killing cancer cells. At the same time, in view of their increased linear energy transfer, heavier ions, like oxygen, are considered attractive, especially for the treatment of hypoxic tumors. While the higher biological dose released in the entrance dose represents in general a drawback for ions heavier than carbon, for oxygen beam this effect may be balanced by the lower dose increase requested to overcome hypoxia.

Conclusions

The possibility of delivering radiation quality-optimized CPT treatments seems to be the new challenge in heavy ion therapy. The potential and limitations of different particle species, according to different sensitivity and morphological scenarios, makes combined treatments of different ions an intriguing option. New ions could open new scenarios in cancer therapy, but would represent as well an opportunity for the treatment of specific non-cancer disease such as atrial fibrillation.

Application of the local effect model to predict DNA double-strand break rejoining after photon and high-LET irradiation

In the recent version of the local effect model (LEM), the biological effects of ionising radiation can be well described trough the consideration of DNA double-strand breaks (DSB) clustering at the micrometre scale. Assuming a giant-loop organisation for the chromatin higher-order structure, two classes of DSB are defined, namely isolated (iDSB) and clustered DSB (cDSB), according to whether exactly one or more than one DSB are induced in a loop, respectively. Here, a DSB kinetic rejoining model based on the LEM is applied to the description of two specific aspects of DSB rejoining, namely the dose dependence of the rejoining capacity after photon radiation and the residual damage observed at late times after ion irradiation. Based on the hypothesis that iDSB and cDSB can be associated to the fast and slow components of rejoining, the model is able to reproduce the experimental data, therefore supporting the relevance of micrometre scale clustering of damage for photon radiation as well as for high-LET radiation.

Modeling Combined Chemotherapy and Particle Therapy for Locally Advanced Pancreatic Cancer

Pancreatic ductal adenocarcinoma is the only cancer for which deaths are predicted to increase in 2014 and beyond. Combined radiochemotherapy protocols using gemcitabine and hypofractionated X-rays are ongoing in several clinical trials. Recent results indicate that charged particle therapy substantially increases local control of resectable and unresectable pancreas cancer, as predicted from previous radiobiology studies considering the high tumor hypoxia. Combination with chemotherapy improves the overall survival (OS). We compared published data on X-ray and charged particle clinical results with or without adjuvant chemotherapy calculating the biological effective dose. We show that chemoradiotherapy with protons or carbon ions results in 1 year OS significantly higher than those obtained with other treatment schedules. Further hypofractionation using charged particles may result in improved local control and survival. A comparative clinical trial using the standard X-ray scheme vs. the best current standard with carbon ions is crucial and may open new opportunities for this deadly disease.

Induction and Processing of the Radiation-Induced Gamma-H2AX Signal and Its Link to the Underlying Pattern of DSB: A Combined Experimental and Modelling Study

We present here an analysis of DSB induction and processing after irradiation with X-rays in an extended dose range based on the use of the γH2AX assay. The study was performed by quantitative flow cytometry measurements, since the use of foci counting would result in reasonable accuracy only in a limited dose range of a few Gy. The experimental data are complemented by a theoretical analysis based on the GLOBLE model. In fact, original aim of the study was to test GLOBLE predictions against new experimental data, in order to contribute to the validation of the model. Specifically, the γH2AX signal kinetics has been investigated up to 24 h after exposure to increasing photon doses between 2 and 500 Gy. The prolonged persistence of the signal at high doses strongly suggests dose dependence in DSB processing after low LET irradiation. Importantly, in the framework of our modelling analysis, this is related to a gradually increased fraction of DSB clustering at the micrometre scale. The parallel study of γH2AX dose response curves shows the onset of a pronounced saturation in two cell lines at a dose of about 20 Gy. This dose is much lower than expected according to model predictions based on the values usually adopted for the DSB induction yield (≈ 30 DSB/Gy) and for the γH2AX foci extension of approximately 2 Mbp around the DSB. We show and discuss how theoretical predictions and experimental findings can be in principle reconciled by combining an increased DSB induction yield with the assumption of a larger genomic extension for the single phosphorylated regions. As an alternative approach, we also considered in our model the possibility of a 3D spreading-mechanism of the H2AX phosphorylation around the induced DSB, and applied it to the analysis of both the aspects considered. Our results are found to be supportive for the basic assumptions on which GLOBLE is built. Apart from giving new insights into the H2AX phosphorylation process, experiments performed at high doses are of relevance in the context of radiation therapy, where hypo-fractionated schemes become increasingly popular.