Carbon and oxygen minibeam radiation therapy: An experimental dosimetric evaluation

Dosimetric study for breathing-induced motion effects in an abdominal pancreas phantom for carbon ion mini-beam radiotherapy

BACKGROUND

Particle mini-beam therapy exhibits promise in sparing healthy tissue through spatial fractionation, particularly notable for heavy ions, further enhancing the already favorable differential biological effectiveness at both target and entrance regions. However, breathing-induced organ motion affects particle mini-beam irradiation schemes since the organ displacements exceed the mini-beam structure dimensions, decreasing the advantages of spatial fractionation.

PURPOSE

In this study, the impact of breathing-induced organ motion on the dose distribution was examined at the target and organs at risk(OARs) during carbon ion mini-beam irradiation for pancreatic cancer.

METHODS

As a first step, the carbon ion mini-beam pattern was characterized with Monte Carlo simulations. To analyze the impact of breathing-induced organ motion on the dose distribution of a virtual pancreas tumor as target and related OARs, the anthropomorphic Pancreas Phantom for Ion beam Therapy (PPIeT) was irradiated with carbon ions. A mini-beam collimator was used to deliver a spatially fractionated dose distribution. During irradiation, varying breathing motion amplitudes were induced, ranging from 5 to 15 mm. Post-irradiation, the 2D dose pattern was analyzed, focusing on the full width at half maximum (FWHM), center-to-center distance (ctc), and the peak-to-valley dose ratio (PVDR).

RESULTS

The mini-beam pattern was visible within OARs, while in the virtual pancreas tumor a more homogeneous dose distribution was achieved. Applied motion affected the mini-beam pattern within the kidney, one of the OARs, reducing the PVDR from 3.78  ± $\pm$  0.12 to 1.478  ± $\pm$  0.070 for the 15 mm motion amplitude. In the immobile OARs including the spine and the skin at the back, the PVDR did not change within 3.4% comparing reference and motion conditions.

CONCLUSIONS

This study provides an initial understanding of how breathing-induced organ motion affects spatial fractionation during carbon ion irradiation, using an anthropomorphic phantom. A decrease in the PVDR was observed in the right kidney when breathing-induced motion was applied, potentially increasing the risk of damage to OARs. Therefore, further studies are needed to explore the clinical viability of mini-beam radiotherapy with carbon ions when irradiating abdominal regions.

Investigating the biochemical response of proton minibeam radiation therapy by means of synchrotron-based infrared microspectroscopy

The biology underlying proton minibeam radiation therapy (pMBRT) is not fully understood. Here we aim to elucidate the biological effects of pMBRT using Fourier Transform Infrared Microspectroscopy (FTIRM). In vitro (CTX-TNA2 astrocytes and F98 glioma rat cell lines) and in vivo (healthy and F98-bearing Fischer rats) irradiations were conducted, with conventional proton radiotherapy and pMBRT. FTIRM measurements were performed at ALBA Synchrotron, and multivariate data analysis methods were employed to assess spectral differences between irradiation configurations and doses. For astrocytes, the spectral regions related to proteins and nucleic acids were highly affected by conventional irradiations and the high-dose regions of pMBRT, suggesting important modifications on these biomolecules. For glioma, pMBRT had a great effect on the nucleic acids and carbohydrates. In animals, conventional radiotherapy had a remarkable impact on the proteins and nucleic acids of healthy rats; analysis of tumour regions in glioma-bearing rats suggested major nucleic acid modifications due to pMBRT.

Spatially fractionated radiation therapy: a critical review on current status of clinical and preclinical studies and knowledge gaps

Spatially fractionated radiation therapy (SFRT) is a therapeutic approach with the potential to disrupt the classical paradigms of conventional radiation therapy. The high spatial dose modulation in SFRT activates distinct radiobiological mechanisms which lead to a remarkable increase in normal tissue tolerances. Several decades of clinical use and numerous preclinical experiments suggest that SFRT has the potential to increase the therapeutic index, especially in bulky and radioresistant tumors. To unleash the full potential of SFRT a deeper understanding of the underlying biology and its relationship with the complex dosimetry of SFRT is needed. This review provides a critical analysis of the field, discussing not only the main clinical and preclinical findings but also analyzing the main knowledge gaps in a holistic way.

Monte Carlo study of the free radical yields in minibeam radiation therapy

BACKGROUND

Minibeam radiation therapy (MBRT) is a novel technique which has been shown to widen the therapeutic window through significant normal tissue sparing. Despite the heterogeneous dose distributions, tumor control is still ensured. Nevertheless the exact radiobiological mechanisms responsible for MBRT efficacy are not fully understood.

PURPOSE

Reactive oxygen species (ROS) resulting from water radiolysis were investigated given their implications not only on targeted DNA damage, but also for their role in the immune response and non-targeted cell signalling effects: two potential drivers of MBRT efficacy.

METHODS

Monte Carlo simulations were performed using TOPAS-nBio to carry out the irradiation of a water phantom with beams of protons (pMBRT), photons (xMBRT), 4 He ions (HeMBRT), and 12 C ions (CMBRT). Primary yields at the end of the chemical stage were calculated in spheres of 20 μm diameter, located in the peaks and valleys at various depths up to the Bragg peak. The chemical stage was limited to 1 ns to approximate biological scavenging, and the yield of · OH, H2 O2 , ande aq - ${\rm e}^{-}_{\rm aq}$ was recorded.

RESULTS

Beyond 10 mm, there were no substantial differences in the primary yields between peaks and valleys of the pMBRT and HeMBRT modalities. For xMBRT, there was a lower primary yield of the radical species · OH ande aq - ${\rm e}^{-}_{\rm aq}$ at all depths in the valleys compared to the peaks, and a higher primary yield of H2 O2 . Compared to the peaks, the valleys of the CMBRT modality were subject to a higher · OH ande aq - ${\rm e}^{-}_{\rm aq}$ yield, and lower H2 O2 yield. This difference between peaks and valleys became more severe in depth. Near the Bragg peak, the increase in the primary yield of the valleys over the peaks was 6% and 4% for · OH ande aq - ${\rm e}^{-}_{\rm aq}$ respectively, while there was a decrease in the yield of H2 O2 by 16%. Given the similar ROS primary yields in the peaks and valleys of pMBRT and HeMBRT, the level of indirect DNA damage is expected to be directly proportional to the peak to valley dose ratio (PVDR). The difference in the primary yields implicates a lower level of indirect DNA damage in the valleys compared to the peaks than what would be suggested by the PVDR for xMBRT, and a higher level for CMBRT.

CONCLUSIONS

These results highlight the notion that depending on the particle chosen, one can expect different levels of ROS in the peaks and valley that goes beyond what would be expected by the macroscopic PVDR. The combination of MBRT with heavier ions is shown to be particularly interesting as the primary yield in the valleys progressively diverges from the level observed in the peaks as the LET increases. While differences in the reported · OH yields of this work implicated the indirect DNA damage, H2 O2 yields particularly implicate non-targeted cell signalling effects, and therefore this work provides a point of reference for future simulations in which the distribution of this species at more biologically relevant timescales could be investigated.

Radiobiological and Treatment-Related Aspects of Spatially Fractionated Radiotherapy

The continuously evolving field of radiotherapy aims to devise and implement techniques that allow for greater tumour control and better sparing of critical organs. Investigations into the complexity of tumour radiobiology confirmed the high heterogeneity of tumours as being responsible for the often poor treatment outcome. Hypoxic subvolumes, a subpopulation of cancer stem cells, as well as the inherent or acquired radioresistance define tumour aggressiveness and metastatic potential, which remain a therapeutic challenge. Non-conventional irradiation techniques, such as spatially fractionated radiotherapy, have been developed to tackle some of these challenges and to offer a high therapeutic index when treating radioresistant tumours. The goal of this article was to highlight the current knowledge on the molecular and radiobiological mechanisms behind spatially fractionated radiotherapy and to present the up-to-date preclinical and clinical evidence towards the therapeutic potential of this technique involving both photon and proton beams.

Divide and conquer: spatially fractionated radiation therapy

Spatially fractionated radiation therapy (SFRT) challenges some of the classical dogmas in conventional radiotherapy. The highly modulated spatial dose distributions in SFRT have been shown to lead, both in early clinical trials and in small animal experiments, to a significant increase in normal tissue dose tolerances. Tumour control effectiveness is maintained or even enhanced in some configurations as compared with conventional radiotherapy. SFRT seems to activate distinct radiobiological mechanisms, which have been postulated to involve bystander effects, microvascular alterations and/or immunomodulation. Currently, it is unclear which is the dosimetric parameter which correlates the most with both tumour control and normal tissue sparing in SFRT. Additional biological experiments aiming at parametrizing the relationship between the irradiation parameters (beam width, spacing, peak-to-valley dose ratio, peak and valley doses) and the radiobiology are needed. A sound knowledge of the interrelation between the physical parameters in SFRT and the biological response would expand its clinical use, with a higher level of homogenisation in the realisation of clinical trials. This manuscript reviews the state of the art of this promising therapeutic modality, the current radiobiological knowledge and elaborates on future perspectives.

Heavy Ion Minibeam Therapy: Side Effects in Normal Brain

The purpose of this work was to investigate whether minibeam therapy with heavy ions might offer improvements of the therapeutic ratio for the treatment of human brain cancers. To assess neurotoxicity, we irradiated normal juvenile rats using 120 MeV lithium-7 ions at an absorbed integral dose of 20 Gy. Beams were configured either as a solid parallel circular beam or as an array of planar parallel minibeams having 300-micron width and 1-mm center-to-center spacing within a circular array. We followed animals for 6 months after treatment and utilized behavioral testing and immunohistochemical studies to investigate the resulting cognitive impairment and chronic pathologic changes. We found both solid-beam therapy and minibeam therapy to result in cognitive impairment compared with sham controls, with no apparent reduction in neurotoxicity using heavy ion minibeams instead of solid beams under the conditions of this study.

Biological and Mechanical Synergies to Deal With Proton Therapy Pitfalls: Minibeams, FLASH, Arcs, and Gantryless Rooms

Proton therapy has advantages and pitfalls comparing with photon therapy in radiation therapy. Among the limitations of protons in clinical practice we can selectively mention: uncertainties in range, lateral penumbra, deposition of higher LET outside the target, entrance dose, dose in the beam path, dose constraints in critical organs close to the target volume, organ movements and cost. In this review, we combine proposals under study to mitigate those pitfalls by using individually or in combination: (a) biological approaches of beam management in time (very high dose rate “FLASH” irradiations in the order of 100 Gy/s) and (b) modulation in space (a combination of mini-beams of millimetric extent), together with mechanical approaches such as (c) rotational techniques (optimized in partial arcs) and, in an effort to reduce cost, (d) gantry-less delivery systems. In some cases, these proposals are synergic (e.g., FLASH and minibeams), in others they are hardly compatible (mini-beam and rotation). Fixed lines have been used in pioneer centers, or for specific indications (ophthalmic, radiosurgery,…), they logically evolved to isocentric gantries. The present proposals to produce fixed lines are somewhat controversial. Rotational techniques, minibeams and FLASH in proton therapy are making their way, with an increasing degree of complexity in these three approaches, but with a high interest in the basic science and clinical communities. All of them must be proven in clinical applications.

Conventional dose rate spatially-fractionated radiation therapy (SFRT) treatment response and its association with dosimetric parameters-A preclinical study in a Fischer 344 rat model

Purpose

To identify key dosimetric parameters that have close associations with tumor treatment response and body weight change in SFRT treatments with a large range of spatial-fractionation scale at dose rates of several Gy/min.

Methods

Six study arms using uniform tumor radiation, half-tumor radiation, 2mm beam array radiation, 0.3mm minibeam radiation, and an untreated arm were used. All treatments were delivered on a 320kV x-ray irradiator. Forty-two female Fischer 344 rats with fibrosarcoma tumor allografts were used. Dosimetric parameters studied are peak dose and width, valley dose and width, peak-to-valley-dose-ratio (PVDR), volumetric average dose, percentage volume directly irradiated, and tumor- and normal-tissue EUD. Animal survival, tumor volume change, and body weight change (indicative of treatment toxicity) are tested for association with the dosimetric parameters using linear regression and Cox Proportional Hazards models.

Results

The dosimetric parameters most closely associated with tumor response are tumor EUD (R² = 0.7923, F-stat = 15.26*; z-test = -4.07***), valley (minimum) dose (R² = 0.7636, F-stat = 12.92*; z-test = -4.338***), and percentage tumor directly irradiated (R² = 0.7153, F-stat = 10.05*; z-test = -3.837***) per the linear regression and Cox Proportional Hazards models, respectively. Tumor response is linearly proportional to valley (minimum) doses and tumor EUD. Average dose (R² = 0.2745, F-stat = 1.514 (no sig.); z-test = -2.811**) and peak dose (R² = 0.04472, F-stat = 0.6874 (not sig.); z-test = -0.786 (not sig.)) show the weakest associations to tumor response. Only the uniform radiation arm did not gain body weight post-radiation, indicative of treatment toxicity; however, body weight change in general shows weak association with all dosimetric parameters except for valley (minimum) dose (R² = 0.3814, F-stat = 13.56**), valley width (R² = 0.2853, F-stat = 8.783**), and peak width (R² = 0.2759, F-stat = 8.382**).

Conclusions

For a single-fraction SFRT at conventional dose rates, valley, not peak, dose is closely associated with tumor treatment response and thus should be used for treatment prescription. Tumor EUD, valley (minimum) dose, and percentage tumor directly irradiated are the top three dosimetric parameters that exhibited close associations with tumor response.

Spatial fractionation of the dose in heavy ions therapy: An optimization study

PURPOSE

The alliance of charged particle therapy and the spatial fractionation of the dose, as in minibeam or Grid therapy, is an innovative strategy to improve the therapeutic index in the treatment of radioresistant tumors. The aim of this work was to assess the optimum irradiation configuration in heavy ion spatially fractionated radiotherapy (SFRT) in terms of ion species, beam width, center-to-center distances, and linear energy transfer (LET), information that could be used to guide the design of the future biological experiments. The nuclear fragmentation leading to peak and valley regions composed of different secondary particles, creates the need for a more complete dosimetric description that the classical one in SFRT.

METHODS

Monte Carlo simulations (GATE 6.2) were performed to evaluate the dose distributions for different ions, beam widths, and spacings. We have also assessed the 3D-maps of dose-averaged LET and proposed a new parameter, the peak-to-valley-LET ratio, to offer a more thorough physical evaluation of the technique.

RESULTS

Our results show that beam widths larger than 400 μm are needed in order to keep a ratio between the dose in the entrance and the dose in the target of the same order as in conventional irradiations. A large ctc distance (3500 μm) would favor tissue sparing since it provides higher PVDR, it leads to a reduced contribution of the heavier nuclear fragments and a LET value in the valleys a factor 2 lower than the LET in the ctc leading to homogeneous distributions in the target.

CONCLUSIONS

Heavy ions MBRT provide advantageous dose distributions. Thanks to the reduced lateral scattering, the use of submillimetric beams still allows to keep a ratio between the dose in the entrance and the dose in the target of the same order as in conventional irradiations. Large ctc distances (3500 μm) should be preferred since they lead to valley doses composed of lighter nuclear fragments resulting in a much reduced dose-averaged LET values in normal tissue, favoring its preservation. Among the different ions species evaluated, Ne stands out as the one leading to the best balance between high PVDR and PVLR in normal tissues and high LET values (close to 100 keV/μm) and a favorable oxygen enhancement ratio in the target region.

Quantitative evaluation of potential irradiation geometries for carbon-ion beam grid therapy

PURPOSE

Radiotherapy using grids containing cm-wide beam elements has been carried out sporadically for more than a century. During the past two decades, preclinical research on radiotherapy with grids containing small beam elements, 25 μm-0.7 mm wide, has been performed. Grid therapy with larger beam elements is technically easier to implement, but the normal tissue tolerance to the treatment is decreasing. In this work, a new approach in grid therapy, based on irradiations with grids containing narrow carbon-ion beam elements was evaluated dosimetrically. The aim formulated for the suggested treatment was to obtain a uniform target dose combined with well-defined grids in the irradiated normal tissue. The gain, obtained by crossfiring the carbon-ion beam grids over a simulated target volume, was quantitatively evaluated.

METHODS

The dose distributions produced by narrow rectangular carbon-ion beams in a water phantom were simulated with the PHITS Monte Carlo code. The beam-element height was set to 2.0 cm in the simulations, while the widths varied from 0.5 to 10.0 mm. A spread-out Bragg peak (SOBP) was then created for each beam element in the grid, to cover the target volume with dose in the depth direction. The dose distributions produced by the beam-grid irradiations were thereafter constructed by adding the dose profiles simulated for single beam elements. The variation of the valley-to-peak dose ratio (VPDR) with depth in water was thereafter evaluated. The separation of the beam elements inside the grids were determined for different irradiation geometries with a selection criterion.

RESULTS

The simulated carbon-ion beams remained narrow down to the depths of the Bragg peaks. With the formulated selection criterion, a beam-element separation which was close to the beam-element width was found optimal for grids containing 3.0-mm-wide beam elements, while a separation which was considerably larger than the beam-element width was found advantageous for grids containing 0.5-mm-wide beam elements. With the single-grid irradiation setup, the VPDRs were close to 1.0 already at a distance of several cm from the target. The valley doses given to the normal tissue at 0.5 cm distance from the target volume could be limited to less than 10% of the mean target dose if a crossfiring setup with four interlaced grids was used.

CONCLUSIONS

The dose distributions produced by grids containing 0.5- and 3.0-mm wide beam elements had characteristics which could be useful for grid therapy. Grids containing mm-wide carbon-ion beam elements could be advantageous due to the technical ease with which these beams can be produced and delivered, despite the reduced threshold doses observed for early and late responding normal tissue for beams of millimeter width, compared to submillimetric beams. The treatment simulations showed that nearly homogeneous dose distributions could be created inside the target volumes, combined with low valley doses in the normal tissue located close to the target volume, if the carbon-ion beam grids were crossfired in an interlaced manner with optimally selected beam-element separations. The formulated selection criterion was found useful for the quantitative evaluation of the dose distributions produced by the different irradiation setups.