Commissioning of a clinical pencil beam scanning proton therapy unit for ultra-high dose rates (FLASH).

PURPOSE: The purpose of this work was to provide a flexible platform for FLASH research with protons by adapting a former clinical pencil beam scanning gantry to irradiations with ultra-high dose rates. METHODS: PSI Gantry 1 treated patients until December 2018. We optimized the beamline parameters to transport the 250 MeV beam extracted from the PSI COMET accelerator to the treatment room, maximizing the transmission of beam intensity to the sample. We characterized a dose monitor on the gantry to ensure good control of the dose, delivered in spot-scanning mode. We characterized the beam for different dose rates and field sizes for transmission irradiations. We explored scanning possibilities in order to enable conformal irradiations or transmission irradiations of large targets (with transverse scanning). RESULTS: We achieved a transmission of 86% from the cyclotron to the treatment room. We reached a peak dose rate of 9000 Gy/s at 3 mm water equivalent depth, along the central axis of a single pencil beam. Field sizes of up to 5 × 5 mm2 were achieved for single-spot FLASH irradiations. Fast transverse scanning allowed to cover a field of 16 × 1.2 cm2 . With the use of a nozzle-mounted range shifter, we are able to span depths in water ranging from 19.6 to 37.9 cm. Various dose levels were delivered with precision within less than 1%. CONCLUSIONS: We have realized a proton FLASH irradiation setup able to investigate continuously a wide dose rate spectrum, from 1 to 9000 Gy/s in single-spot irradiation as well as in the pencil beam scanning mode. As such, we have developed a versatile test bench for FLASH research.

Towards FLASH proton therapy: the impact of treatment planning and machine characteristics on achievable dose rates.

Background: This study aimed at evaluating spatially varying instantaneous dose rates for different intensity-modulated proton therapy (IMPT) planning strategies and delivery scenarios, and comparing these with FLASH dose rates (>40 Gy/s). Material and methods: In order to quantify dose rates in three-dimensions, we proposed the 'dose-averaged dose rate' (DADR) metric, defined for each voxel as the dose-weighted mean of the instantaneous dose rates of all spots (i.e., pencil beams). This concept was applied to four head-and-neck cases, each planned with clinical (4 fields) and various spot-reduced IMPT techniques: 'standard' (4 fields), 'arc' (120 fields) and 'arc-shoot-through' (120 fields; 229 MeV only). For all plans, different delivery scenarios were simulated: constant beam intensity, variable beam intensity for a clinical Varian ProBeam system, varied per energy layer or per spot, and theoretical spot-wise variable beam intensity (i.e., no monitor/safety limitations). DADR distributions were calculated assuming 2-Gy or 6-Gy fractions. Results: Spot-reduced plans contained 17-52 times fewer spots than clinical plans, with no deterioration of plan quality. For the clinical plans, the mean DADR in normal tissue for 2-Gy fractionation was 1.7 Gy/s (median over all patients) at maximum, whereas in standard spot-reduced plans it was 0.7, 4.4, 7.1, and 12.1 Gy/s, for the constant, energy-layer-wise, spot-wise, and theoretical spot-wise delivery scenarios, respectively. Similar values were observed for arc plans. Arc-shoot-through planning resulted in DADR values of 3.0, 6.0, 14.1, and 24.4 Gy/s, for the abovementioned scenarios. Hypofractionation (3×) generally resulted in higher dose rates, up to 73.2 Gy/s for arc-shoot-through plans. The DADR was inhomogeneously distributed with highest values at beam entrance and at the Bragg peak. Conclusion: FLASH dose rates were not achieved for conventional planning and clinical spot-scanning machines. As such, increased spot-wise beam intensities, spot-reduced planning, hypofractionation and arc-shoot-through plans were required to achieve FLASH compatible dose rates.

Sparing the region of the salivary gland containing stem cells preserves saliva production after radiotherapy for head and neck cancer.

Each year, 500,000 patients are treated with radiotherapy for head and neck cancer, resulting in relatively high survival rates. However, in 40% of patients, quality of life is severely compromised because of radiation-induced impairment of salivary gland function and consequent xerostomia (dry mouth). New radiation treatment technologies enable sparing of parts of the salivary glands. We have determined the parts of the major salivary gland, the parotid gland, that need to be spared to ensure that the gland continues to produce saliva after irradiation treatment. In mice, rats, and humans, we showed that stem and progenitor cells reside in the region of the parotid gland containing the major ducts. We demonstrated in rats that inclusion of the ducts in the radiation field led to loss of regenerative capacity, resulting in long-term gland dysfunction with reduced saliva production. Then we showed in a cohort of patients with head and neck cancer that the radiation dose to the region of the salivary gland containing the stem/progenitor cells predicted the function of the salivary glands one year after radiotherapy. Finally, we showed that this region of the salivary gland could be spared during radiotherapy, thus reducing the risk of post-radiotherapy xerostomia.

Whole breast proton irradiation for maximal reduction of heart dose in breast cancer patients.

PURPOSE: In left-sided breast cancer radiotherapy, tangential intensity modulated radiotherapy combined with breath-hold enables a dose reduction to the heart and left anterior descending (LAD) coronary artery. Aim of this study was to investigate the added value of intensity modulated proton therapy (IMPT) with regard to decreasing the radiation dose to these structures. METHODS: In this comparative planning study, four treatment plans were generated in 20 patients: an IMPT plan and a tangential IMRT plan, both with breath-hold and free-breathing. At least 97 % of the target volume had to be covered by at least 95 % of the prescribed dose in all cases. Specifically with respect to the heart, the LAD, and the target volumes, we analyzed the maximum doses, the mean doses, and the volumes receiving 5-30 Gy. RESULTS: As compared to IMRT, IMPT resulted in significant dose reductions to the heart and LAD-region even without breath-hold. In the majority of the IMPT cases, a reduction to almost zero to the heart and LAD-region was obtained. IMPT treatment plans yielded the lowest dose to the lungs. CONCLUSIONS: With IMPT the dose to the heart and LAD-region could be significantly decreased compared to tangential IMRT with breath-hold. The clinical relevance should be assessed individually based on the baseline risk of cardiac complications in combination with the dose to organs at risk. However, as IMPT for breast cancer is currently not widely available, IMPT should be reserved for patients remaining at high risk for major coronary events.

Emerging technologies in proton therapy.

An increasing number of proton therapy facilities are being planned and built at hospital based centers. Most facilities are employing traditional dose delivery methods. A second generation of dose application techniques, based on pencil beam scanning, is slowly being introduced into the commercially available proton therapy systems. New developments in accelerator physics are needed to accommodate and fully exploit these new techniques. At the same time new developments such as the development of small cyclotrons, Dielectric Wall Accelerator (DWA) and laser driven systems, aim for smaller, single room treatment units. In general the benefits of proton therapy could be exploited optimally when achieving a higher level in accuracy, beam energy, beam intensity, safety and system reliability. In this review an overview of the current developments will be given followed by a discussion of upcoming new technologies and needs, like increase of energy, on-line MRI and proton beam splitting for independent uses of treatment rooms.

Volume-dependent expression of in-field and out-of-field effects in the proton-irradiated rat lung.

PURPOSE: To investigate whether occurrence of early radiation effects in lung tissue depends on local dose only. METHODS AND MATERIALS: Twenty-five percent, 50%, 66%, 88%, or 100% of the rat lung was irradiated using single fractions of 150-MeV protons. For all volumes, in-field and out-of-field dose-response curves were obtained 8 weeks after irradiation. The pathohistology of parenchymal inflammation, infiltrates, fibrosis, and vascular damage and the relative expression of proinflammatory cytokines interleukin (IL)-1α, transforming growth factor-ß, IL-6, and tumor necrosis factor-α were assessed. RESULTS: For all histologic endpoints, irradiated dose- and volume-dependent in-field and out-of-field effects were observed, albeit with different dynamics. Of note, the out-of-field effects for vascular damage were very similar to the in-field effects. Interestingly, only IL-6 showed a clear dose-dependent increase in expression both in-field and out-of-field, whereas the expression levels of IL-1α, transforming growth factor-ß, and tumor necrosis factor-α were either very low or without a clear dose-volume relation. As such, none of the radiation effects studied depended only on local dose to the tissue. CONCLUSION: The effects of radiation to lung tissue do not only depend on local dose to that tissue. Especially at high-volume irradiation, lung damage seems to present globally rather than locally. The accuracy of predictive modeling may be improved by including nonlocal effects.

Quantifying local radiation-induced lung damage from computed tomography.

PURPOSE: Optimal implementation of new radiotherapy techniques requires accurate predictive models for normal tissue complications. Since clinically used dose distributions are nonuniform, local tissue damage needs to be measured and related to local tissue dose. In lung, radiation-induced damage results in density changes that have been measured by computed tomography (CT) imaging noninvasively, but not yet on a localized scale. Therefore, the aim of the present study was to develop a method for quantification of local radiation-induced lung tissue damage using CT. METHODS AND MATERIALS: CT images of the thorax were made 8 and 26 weeks after irradiation of 100%, 75%, 50%, and 25% lung volume of rats. Local lung tissue structure (S(L)) was quantified from local mean and local standard deviation of the CT density in Hounsfield units in 1-mm(3) subvolumes. The relation of changes in S(L) (DeltaS(L)) to histologic changes and breathing rate was investigated. Feasibility for clinical application was tested by applying the method to CT images of a patient with non-small-cell lung carcinoma and investigating the local dose-effect relationship of DeltaS(L). RESULTS: In rats, a clear dose-response relationship of DeltaS(L) was observed at different time points after radiation. Furthermore, DeltaS(L) correlated strongly to histologic endpoints (infiltrates and inflammatory cells) and breathing rate. In the patient, progressive local dose-dependent increases in DeltaS(L) were observed. CONCLUSION: We developed a method to quantify local radiation-induced tissue damage in the lung using CT. This method can be used in the development of more accurate predictive models for normal tissue complications.

Bath and shower effects in the rat parotid gland explain increased relative risk of parotid gland dysfunction after intensity-modulated radiotherapy.

PURPOSE: To assess in a rat model whether adding a subtolerance dose in a region adjacent to a high-dose irradiated subvolume of the parotid gland influences its response (bath-and-shower effect). METHODS AND MATERIALS: Irradiation of the whole, cranial 50%, and/or the caudal 50% of the parotid glands of Wistar rats was performed using 150-MeV protons. To determine suitable (i.e., subtolerance) dose levels for a bath-dose, both whole parotid glands were irradiated with 5 to 25 Gy. Subsequently groups of Wistar rats received 30 Gy to the caudal 50% (shower) and 0 to 10 Gy to the cranial 50% (bath) of both parotid glands. Stimulated saliva flow rate (function) was measured before and up to 240 days after irradiation. RESULTS: Irradiation of both glands up to a dose of 10 Gy did not result in late loss of function and is thus regarded subtolerance. Addition of a dose bath of 1 to 10 Gy to a high-dose in the caudal 50% of the glands resulted in enhanced function loss. CONCLUSION: Similar to the spinal cord, the parotid gland demonstrates a bath and shower effect, which may explain the less-than-expected sparing of function after IMRT.

The impact of heart irradiation on dose-volume effects in the rat lung.

PURPOSE: To test the hypothesis that heart irradiation increases the risk of a symptomatic radiation-induced loss of lung function (SRILF) and that this can be well-described as a modulation of the functional reserve of the lung. METHODS AND MATERIALS: Rats were irradiated with 150-MeV protons. Dose-response curves were obtained for a significant increase in breathing frequency after irradiation of 100%, 75%, 50%, or 25% of the total lung volume, either including or excluding the heart from the irradiation field. A significant increase in the mean respiratory rate after 6-12 weeks compared with 0-4 weeks was defined as SRILF, based on biweekly measurements of the respiratory rate. The critical volume (CV) model was used to describe the risk of SRILF. Fits were done using a maximum likelihood method. Consistency between model and data was tested using a previously developed goodness-of-fit test. RESULTS: The CV model could be fitted consistently to the data for lung irradiation only. However, this fitted model failed to predict the data that also included heart irradiation. Even refitting the model to all data resulted in a significant difference between model and data. These results imply that, although the CV model describes the risk of SRILF when the heart is spared, the model needs to be modified to account for the impact of dose to the heart on the risk of SRILF. Finally, a modified CV model is described that is consistent to all data. CONCLUSIONS: The detrimental effect of dose to the heart on the incidence of SRILF can be described by a dose dependent decrease in functional reserve of the lung.

Influence of adjacent low-dose fields on tolerance to high doses of protons in rat cervical spinal cord.

PURPOSE: The dose-response relationship for a relatively short length (4 mm) of rat spinal cord has been shown to be significantly modified by adjacent low-dose fields. In an additional series of experiments, we have now established the dose-volume dependence of this effect. METHODS AND MATERIALS: Wistar rats were irradiated on the cervical spinal cord with single doses of unmodulated protons (150 MeV) to obtain sharp lateral penumbras, by use of the shoot-through technique, which employs the plateau of the depth-dose profile rather than the Bragg peak. Three types of inhomogeneous dose distributions were administered: Twenty millimeters of cervical spinal cord were irradiated with variable subthreshold (= bath) doses (4 and 18 Gy). At the center of the 20-mm segment, a short segment of 2 mm or 8 mm (= shower) was irradiated with variable single doses. These inhomogeneous dose distributions are referred to as symmetrical bath-and-shower experiments. An asymmetrical dose distribution was arranged by irradiation of 12 mm (= bath) of spinal cord with a dose of 4 Gy. The caudal 2 mm (= shower) of the 12-mm bath was additionally irradiated with variable single doses. This arrangement of inhomogeneous dose distribution is referred to as asymmetrical bath-and-shower experiment. The endpoint for estimation of the dose-response relationships was paralysis of the fore limbs or hind limbs and confirmation by histology. RESULTS: The 2-mm bath-and-shower experiments with a 4-Gy bath dose showed a large shift of the dose-response curves compared with the 2-mm single field, which give lower ED50 values of 61.2 Gy and 68.6 Gy for the symmetrical and asymmetrical arrangement, respectively, compared with an ED50 of 87.8 Gy after irradiation of a 2-mm field only. If the bath dose is increased to 18 Gy, the ED50 value is decreased further to 30.9 Gy. For an 8-mm field, addition of a 4-Gy bath dose did not modify the ED50 obtained for an 8-mm field only (23.2 and 23.1 Gy). CONCLUSIONS: The spinal cord tolerance of relatively small volumes (shower) is strongly affected by low-dose irradiation (= bath) of adjacent tissue. The results of all bath-and-shower experiments show the effect of a low bath dose to be highest for a field of 2 mm, less for 4 mm, and absent for 8 mm. Adding a 4-Gy bath to only 1 side of a 2-mm field still showed a large effect. Because glial progenitor cells are known to migrate over at least 2 to 3 mm, this observation indicates that interference with stem cell migration is not the most likely mechanism of a bath effect.

Radiation damage to the heart enhances early radiation-induced lung function loss.

In many thoracic cancers, the radiation dose that can safely be delivered to the target volume is limited by the tolerance dose of the surrounding lung tissue. It has been hypothesized that irradiation of the heart may be an additional risk factor for the development of early radiation-induced lung morbidity. In the current study, the dependence of lung tolerance dose on heart irradiation is determined. Fifty percent of the rat lungs were irradiated either including or excluding the heart. Proton beams were used to allow very accurate and conformal dose delivery. Lung function toxicity was scored using a breathing rate assay. We confirmed that the tolerance dose for early lung function damage depends not only on the lung region that is irradiated but also that concomitant irradiation of the heart severely reduces the tolerance of the lung. This study for the first time shows that the response of an organ to irradiation does not necessarily depend on the dose distribution in that organ alone.

Data on dose-volume effects in the rat spinal cord do not support existing NTCP models.

PURPOSE: To evaluate several existing dose-volume effect models for their ability to describe the occurrence of white matter necrosis in rat spinal cord after irradiation with small proton beams. METHODS AND MATERIALS: A large number of dose-volume effect models has been fitted to data on the occurrence of white matter necrosis after irradiation with small proton beams. The fitting was done with the maximum likelihood method. For each model, the goodness of fit was calculated. An empirical tolerance dose-volume (eTDV) model was designed to describe data obtained after uniform irradiation. RESULTS: The eTDV model, the critical element model, and critical volume model with inclusion of the repair-by-migration principle described by Shirato, were able to describe the data obtained after irradiation with uniform dose distributions of varying sizes. However, none of the models under investigation was able to describe all the data. Extension of the developed empirical model with a repair mechanism with a limited range resulted in a good description of the tolerance doses. CONCLUSIONS: In the rat spinal cord, a nonlocal repair mechanism, acting from nonirradiated to irradiated tissue, plays an important role in the (prevention of the) occurrence of white matter necrosis after irradiation. Models that take into account this effect need to be developed.

Regional differences in radiosensitivity across the rat cervical spinal cord.

PURPOSE: To study regional differences in radiosensitivity within the rat cervical spinal cord. METHODS AND MATERIALS: Three types of inhomogeneous dose distributions were applied to compare the radiosensitivity of the lateral and central parts of the rat cervical spinal cord. The left lateral half of the spinal cord was irradiated with two grazing proton beams, each with a different penumbra (20-80% isodoses): lateral wide (penumbra = 1.1 mm) and lateral tight (penumbra = 0.8 mm). In the third experiment, the midline of the cord was irradiated with a narrow proton beam with a penumbra of 0.8 mm. The irradiated spinal cord length (C1-T2) was 20 mm in all experiments. The animals were irradiated with variable single doses of unmodulated protons (150 MeV) with the shoot-through method, whereby the plateau of the depth-dose profile is used rather than the Bragg peak. The endpoint for estimating isoeffective dose (ED(50)) values was paralysis of fore and/or hind limbs within 210 days after irradiation. Histology of the spinal cords was performed to assess the radiation-induced tissue damage. RESULTS: High-precision proton irradiation of the lateral or the central part of the spinal cord resulted in a shift of dose-response curves to higher dose values compared with the homogeneously irradiated cervical cord to the same 20-mm length. The ED(50) values were 28.9 Gy and 33.4 Gy for the lateral wide and lateral tight irradiations, respectively, and as high as 71.9 Gy for the central beam experiment, compared with 20.4 Gy for the homogeneously irradiated 20-mm length of cervical cord. Histologic analysis of the spinal cords showed that the paralysis was due to white matter necrosis. The radiosensitivity was inhomogeneously distributed across the spinal cord, with a much more radioresistant central white matter (ED(50) = 71.9 Gy) compared with lateral white matter (ED(50) values = 28.9 Gy and 33.4 Gy). The gray matter did not show any noticeable lesions, such as necrosis or hemorrhage, up to 80 Gy. All lesions induced were restricted to white matter structures. CONCLUSIONS: The observed large regional differences in radiosensitivity within the rat cervical spinal cord indicate that the lateral white matter is more radiosensitive than the central part of the white matter. The gray matter is highly resistant to radiation: no lesions observable by light microscopy were induced, even after a single dose as high as 80 Gy.

Unexpected changes of rat cervical spinal cord tolerance caused by inhomogeneous dose distributions.

PURPOSE: The effects of dose distribution on dose-effect relationships have been evaluated and, from this, iso-effective doses (ED(50)) established. METHODS AND MATERIALS: Wistar rats were irradiated on the cervical spinal cord with single doses of unmodulated protons (150 MeV) to obtain sharp lateral penumbras, using the shoot-through technique, which employs the plateau of the depth-dose profile rather than the Bragg peak. Two types of inhomogeneous dose distributions have been administered: (1) 2 4-mm fields with 8- or 12-mm spacing between the center of the fields (referred to as split-field) were irradiated with variable single doses and (2) cervical spinal cord was irradiated with various combinations of relatively low doses to a large volume (20 mm) combined with high doses to a small volume (4 mm) (referred to as bath and shower). The endpoint for estimating the dose-response relationships was paralysis of the fore or hind limbs. RESULTS: The split-field experiments (2 x 4 mm) showed a shift in the dose-response curves, giving significant higher ED(50) values of 45.4 Gy and 41.6 Gy for 8- and 12-mm spacing, respectively, compared with the ED(50) of 24.9 Gy for the single 8 mm (same total tissue volume irradiated). These values were closer to the ED(50) for a single 4-mm field of 53.7 Gy. The bath and shower experiments showed a large decrease of the ED(50) values from 15-22 Gy when compared with the 4-mm single field, even with a bath dose as low as 4 Gy. There were no histologic changes found in the low dose bath regions of the spinal cord at postmortem. CONCLUSIONS: Not only the integral irradiated volume is a determining factor for the ED(50) of rat cervical spinal cord, but also the shape of the dose distribution is of great importance. The high ED(50) values of a small region or shower (4 mm) decreases significantly when the adjacent tissue is irradiated with a subthreshold dose (bath), even as low as 4 Gy. The significant shift to lower ED(50) values for induction of paralysis of the limbs by adding a low-dose bath was not accompanied by changes in histologic lesions. These observations may have implications for the interpretation of complex treatment plans and normal tissue complication probability in intensity-modulated radiotherapy.

Dose-volume effects in the rat cervical spinal cord after proton irradiation.

PURPOSE: To estimate dose-volume effects in the rat cervical spinal cord with protons. METHODS AND MATERIALS: Wistar rats were irradiated on the cervical spinal cord with a single fraction of unmodulated protons (150-190 MeV) using the shoot through method, which employs the plateau of the depth-dose profile rather than the Bragg peak. Four different lengths of the spinal cord (2, 4, 8, and 20 mm) were irradiated with variable doses. The endpoint for estimating dose-volume effects was paralysis of fore or hind limbs. RESULTS: The results obtained with a high-precision proton beam showed a marginal increase of ED50 when decreasing the irradiated cord length from 20 mm (ED50 = 20.4 Gy) to 8 mm (ED50 = 24.9 Gy), but a steep increase in ED50 when further decreasing the length to 4 mm (ED50 = 53.7 Gy) and 2 mm (ED50 = 87.8 Gy). These results generally confirm data obtained previously in a limited series with 4-6-MV photons, and for the first time it was possible to construct complete dose-response curves down to lengths of 2 mm. At higher ED50 values and shorter lengths irradiated, the latent period to paralysis decreased from 125 to 60 days. CONCLUSIONS: Irradiation of variable lengths of rat cervical spinal cord with protons showed steeply increasing ED50 values for lengths of less than 8 mm. These results suggest the presence of a critical migration distance of 2-3 mm for cells involved in regeneration processes.