Improving the dose distributions in minibeam radiation therapy: Helium ions vs protons

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.

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.

Longitudinally Heterogeneous Tumor Dose Optimizes Proton Broadbeam, Interlaced Minibeam, and FLASH Therapy

Simple Summary

The aim of any kind of external radiation therapy is to control a tumor with the highest possible probability of the lowest possible side effects. Here, we study further opportunities of reducing the side effects of proton therapy by applying longitudinally heterogeneous dose distributions in the tumor respecting the delivery of a minimum prescribed dose. In our simulations, the longitudinally heterogeneous dose distributions show a reduced dose in the healthy tissue already in the case of proton broadbeam irradiations, but a much higher (calculated) mean cell survival in the case of proton minibeam irradiation. This demonstrates its potential to substantially reduce side effects at a simultaneously higher tumor control probability, opening new opportunities of easier application when striving for high dose-rate applications of proton beams (>~10 Gy/s), in order to additionally profit from the so-called FLASH effects.

Abstract

The prerequisite of any radiation therapy modality (X-ray, electron, proton, and heavy ion) is meant to meet at least a minimum prescribed dose at any location in the tumor for the best tumor control. In addition, there is also an upper dose limit within the tumor according to the International Commission on Radiation Units (ICRU) recommendations in order to spare healthy tissue as well as possible. However, healthy tissue may profit from the lower side effects when waving this upper dose limit and allowing a larger heterogeneous dose deposition in the tumor, but maintaining the prescribed minimum dose level, particularly in proton minibeam therapy. Methods: Three different longitudinally heterogeneous proton irradiation modes and a standard spread-out Bragg peak (SOBP) irradiation mode are simulated for their depth-dose curves under the constraint of maintaining a minimum prescribed dose anywhere in the tumor region. Symmetric dose distributions of two opposing directions are overlaid in a 25 cm-thick water phantom containing a 5 cm-thick tumor region. Interlaced planar minibeam dose distributions are compared to those of a broadbeam using the same longitudinal dose profiles. Results and Conclusion: All longitudinally heterogeneous proton irradiation modes show a dose reduction in the healthy tissue compared to the common SOBP mode in the case of broad proton beams. The proton minibeam cases show eventually a much larger mean cell survival and thus a further reduced equivalent uniform dose (EUD) in the healthy tissue than any broadbeam case. In fact, the irradiation mode using only one proton energy from each side shows better sparing capabilities in the healthy tissue than the common spread-out Bragg peak irradiation mode with the option of a better dose fall-off at the tumor edges and an easier technical realization, particularly in view of proton minibeam irradiation at ultra-high dose rates larger than ~10 Gy/s (so-called FLASH irradiation modes).

Radiation quality correction factors for improved dosimetry in preclinical minibeam radiotherapy

BACKGROUND

In reference dosimetry, radiation quality correction factors are used in order to account for changes in the detector's response among different radiation qualities, improving dosimetric accuracy.

PURPOSE

Reference dosimetry radiation quality corrections factors for the PTW microDiamond were calculated for preclinical X-ray and proton minibeams, and their impact in dosimetric accuracy was evaluated.

METHODS

A formalism for the calculation of radiation quality correction factors for absolute dosimetry in minibeam fields was developed. Following our formalism, radiation quality correction factors were calculated for the PTW microDiamond detector, using the Monte Carlo method. Models of the detector, and X-ray and proton irradiation platform, were imported into the TOPAS Monte Carlo simulation toolkit. The radiation quality correction factors were calculated in the following scenarios: (i) reference dosimetry open field to minibeam center of the central peak, (ii) different positions at the minibeam profile (along the peaks and valleys direction) to the center of the central minibeam, and (iii) some representative depth positions. In addition, the radiation quality correction factors needed for the calculation of the peak-to-valley dose ratio at different depths were calculated.

RESULTS

An important overestimation of the dose (about 10%) was found in the case of the open to minibeam field for both X-rays and proton beams, when the correction factors were used. Smaller differences were observed in the other cases.

CONCLUSIONS

The usage of the PTW microDiamond detector requires radiation quality correction factors in order to be used in minibeam reference dosimetry.

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.

Proton Minibeam Radiation Therapy and Arc Therapy: Proof of Concept of a Winning Alliance

Simple Summary

Normal tissue’s morbidity continues to limit the increase in the therapeutic index in radiation therapy. This study explores the potential advantages of combining proton arc therapy and proton minibeam radiation therapy, which have already individually shown a significant normal tissue’s sparing. This alliance aims to integrate the benefits of those techniques in a single approach.

Abstract

(1) Background: Proton Arc Therapy and Proton Minibeam Radiation Therapy are two novel therapeutic approaches with the potential to lower the normal tissue complication probability, widening the therapeutic window for radioresistant tumors. While the benefits of both modalities have been individually evaluated, their combination and its potential advantages are being assessed in this proof-of-concept study for the first time. (2) Methods: Monte Carlo simulations were employed to evaluate the dose and LET distributions in brain tumor irradiations. (3) Results: a net reduction in the dose to normal tissues (up to 90%), and the preservation of the spatial fractionation of the dose were achieved for all configurations evaluated. Additionally, Proton Minibeam Arc Therapy (pMBAT) reduces the volumes exposed to high-dose and high-LET values at expense of increased low-dose and intermediate-LET values. (4) Conclusions: pMBAT enhances the individual benefits of proton minibeams while keeping those of conventional proton arc therapy. These results might facilitate the path towards patients’ treatments since lower peak doses in normal tissues would be needed than in the case of a single array of proton minibeams.

Technical note: Providing proton fields down to the few-MeV level at clinical pencil beam scanning facilities for radiobiological experiments

PURPOSE

The adequate performance of radiobiological experiments using clinical proton beams typically requires substantial preparations to provide the appropriate setup for specific experiments. Providing radiobiologically interesting low-energy protons is a particular challenge, due to various physical effects that become more pronounced with larger absorber thickness and smaller proton energy. This work demonstrates the generation of decelerated low-energy protons from a clinical proton beam.

METHODS

Monte Carlo simulations of proton energy spectra were performed for energy absorbers with varying thicknesses to reduce the energy of the clinical proton beam down to the few-MeV level corresponding to μ m-ranges. In this way, a setup with an optimum thickness of the absorber with a maximum efficiency of the proton fluence for the provisioning of low-energy protons is supposed to be found. For the specific applications of 2.5-3.3 MeV protons and α -particle range equivalent protons, the relative depth dose was measured and simulated together with the dose-averaged linear energy transfer (LETd) distribution.

RESULTS

The resulting energy spectra from Monte Carlo simulations indicate an optimal absorber thickness for providing low-energy protons with maximum efficiency of proton fluence at an user-requested energy range for experiments. For instance, providing energies lower than 5 MeV, an energy spectrum with a relative total efficiency of 38.6 % to the initial spectrum was obtained with the optimal setup. The measurements of the depth dose, compared to the Monte Carlo simulations, showed that the dosimetry of low-energy protons works and protons with high LETd down to the range of α -particles can be produced.

CONCLUSIONS

This work provides a method for generating all clinically and radiobiologically relevant energies - especially down to the few-MeV level - at one clinical facility with pencil beam scanning. Thereby, it enables radiobiological experiments under environmentally uniform conditions.

Optimizing proton minibeam radiotherapy by interlacing and heterogeneous tumor dose on the basis of calculated clonogenic cell survival

Proton minibeam radiotherapy (pMBRT) is a spatial fractionation method using sub-millimeter beams at center-to-center (ctc) distances of a few millimeters to widen the therapeutic index by reduction of side effects in normal tissues. Interlaced minibeams from two opposing or four orthogonal directions are calculated to minimize side effects. In particular, heterogeneous dose distributions applied to the tumor are investigated to evaluate optimized sparing capabilities of normal tissues at the close tumor surrounding. A 5 cm thick tumor is considered at 10 cm depth within a 25 cm thick water phantom. Pencil and planar minibeams are interlaced from two (opposing) directions as well as planar beams from four directions. An initial beam size of σ0 = 0.2 mm (standard deviation) is assumed in all cases. Tissue sparing potential is evaluated by calculating mean clonogenic cell survival using a linear-quadratic model on the calculated dose distributions. Interlacing proton minibeams for homogeneous irradiation of the tumor has only minor benefits for the mean clonogenic cell survival compared to unidirectional minibeam irradiation modes. Enhanced mean cell survival, however, is obtained when a heterogeneous dose distribution within the tumor is permitted. The benefits hold true even for an elevated mean tumor dose, which is necessary to avoid cold spots within the tumor in concerns of a prescribed dose. The heterogeneous irradiation of the tumor allows for larger ctc distances. Thus, a high mean cell survival of up to 47% is maintained even close to the tumor edges for single fraction doses in the tumor of at least 10 Gy. Similar benefits would result for heavy ion minibeams with the advantage of smaller minibeams in deep tissue potentially offering even increased tissue sparing. The enhanced mean clonogenic cell survival through large ctc distances for interlaced pMBRT with heterogeneous tumor dose distribution results in optimum tissue sparing potential. The calculations show the largest enhancement of the mean cell survival in normal tissue for high-dose fractions. Thus, hypo-fractionation or even single dose fractions become possible for tumor irradiation. A widened therapeutic index at big cost reductions is offered by interlaced proton or heavy ion minibeam therapy.

First proton minibeam radiation therapy treatment plan evaluation

Proton minibeam radiation therapy (pMBRT) is a novel dose delivery method based on spatial dose fractionation. pMBRT has been shown to be promising in terms of reduced side effects and superior tumour control in high-grade glioma-bearing rats compared to standard irradiation. These findings, together with the recent optimized implementation of pMBRT in a clinical pencil beam scanning system, have triggered reflection on the possible application to patient treatments. In this context, the present study was designed to conduct a first theoretical investigation of the clinical potential of this technique. For this purpose, a dedicated dose engine was developed and used to evaluate two clinically relevant patient treatment plans (high-grade glioma and meningioma). Treatment plans were compared with standard proton therapy plans assessed by means of a commercial treatment planning system (ECLIPSE-Varian Medical systems) and Monte Carlo simulations. A multislit brass collimator consisting of 0.4 mm wide slits separated by a centre-to-centre distance of 4 or 6 mm was placed between the nozzle and the patient to shape the planar minibeams. For each plan, spread-out Bragg peaks and homogeneous dose distributions (±7% dose variations) can be obtained in target volumes. The Peak-to-Valley Dose Ratios (PVDR) were evaluated between 9.2 and 12.8 at a depth of 20 mm for meningioma and glioma, respectively. Dose volume histograms (DVHs) for target volumes and organs at risk were quantitatively compared, resulting in a slightly better target homogeneity with standard PT than with pMBRT plans, but similar DVHs for deep-seated organs-at-risk and lower average dose for shallow organs. The proposed delivery method evaluated in this work opens the way to an effective treatment for radioresistant tumours and will support the design of future clinical research.