Optimization of beam arrangements in proton minibeam radiotherapy by cell survival simulations

Tuning spatially fractionated radiotherapy dose profiles using the moiré effect

Spatially Fractionated Radiotherapy (SFRT) has demonstrated promising potential in cancer treatment, combining the advantages of reduced post-radiation effects and enhanced local control rates. Within this paradigm, proton minibeam radiotherapy (pMBRT) was suggested as a new treatment modality, possibly producing superior normal tissue sparing to conventional proton therapy, leading to improvements in patient outcomes. However, an effective and convenient beam generation method for pMBRT, capable of implementing various optimum dose profiles, is essential for its real-world application. Our study investigates the potential of utilizing the moiré effect in a dual collimator system (DCS) to generate pMBRT dose profiles with the flexibility to modify the center-to-center distance (CTC) of the dose distribution in a technically simple way.We employ the Geant4 Monte Carlo simulations tool to demonstrate that the angle between the two collimators of a DCS can significantly impact the dose profile. Varying the DCS angle from 10∘ to 50∘ we could cover CTC ranging from 11.8 mm to 2.4 mm, respectively. Further investigations reveal the substantial influence of the multi-slit collimator's (MSC) physical parameters on the spatially fractionated dose profile, such as period (CTC), throughput, and spacing between MSCs. These findings highlight opportunities for precision dose profile adjustments tailored to specific clinical scenarios.The DCS capacity for rapid angle adjustments during the energy transition stages of a spot scanning system can facilitate dynamic alterations in the irradiation profile, enhancing dose contrast in normal tissues. Furthermore, its unique attribute of spatially fractionated doses in both lateral directions could potentially improve normal tissue sparing by minimizing irradiated volume. Beyond the realm of pMBRT, the dual MSC system exhibits remarkable versatility, showing compatibility with different types of beams (X-rays and electrons) and applicability across various SFRT modalities.Our study illuminates the dual MSC system's potential as an efficient and adaptable tool in the refinement of pMBRT techniques. By enabling meticulous control over irradiation profiles, this system may expedite advancements in clinical and experimental applications, thereby contributing to the evolution of SFRT strategies.

TVL1-IMPT: Optimization of Peak-to-Valley Dose Ratio Via Joint Total-Variation and L1 Dose Regularization for Spatially Fractionated Pencil-Beam-Scanning Proton Therapy

PURPOSE

Proton minibeam radiation therapy (pMBRT) is a novel proton modality of spatially fractionated RT. pMBRT can reduce the radiation damage to normal tissues via biological dose sparing of high peak-to-valley dose ratio (PVDR). This work will develop a new spatially fractionated IMPT treatment planning method for pMBRT that jointly optimizes the plan quality and maximizes the PVDR.

METHODS

The new optimization method simultaneously maximizes the normal-tissue PVDR and optimizes the dose distribution at tumor targets and organs at risk. The PVDR maximization is through the joint total variation (TV) and L1 regularization with respect to the normal-tissue dose. That is, the beam-eye view projects dose slices of several depths for each beam angle; the TV of dose is maximized, corresponding to the PVDR maximization; and the L1 of dose is minimized, corresponding to the minimization of the organs-at-risk dose and maximization of survival fraction (SF).

RESULTS

The new IMPT method with TV and L1 regularization was validated in comparison with the conventional IMPT method for pMBRT in several clinical cases. The results show that TVL1 provided larger PVDR and SF than the conventional IMPT method for biological sparing of normal tissues, with preserved plan quality in terms of physical dose distribution.

CONCLUSIONS

A new spatially fractionated IMPT treatment planning method was developed for pMBRT that can optimize and improve normal-tissue PVDR and SF by incorporating TV and L1 dose regularization with properly chosen regularization parameters into IMPT.

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).

Monte Carlo evaluation of high-gradient magnetically focused planar proton minibeams in a passive nozzle

Objective. To investigate the potential of using a single quadrupole magnet with a high magnetic field gradient to create planar minibeams suitable for clinical applications of proton minibeam radiation therapy.Approach. We performed Monte Carlo simulations involving single quadrupole Halbach cylinders in a passively scattered nozzle in clinical use for proton therapy. Pencil beams produced by the nozzle of 10-15 mm initial diameters and particle range of ∼10-20 cm in water were focused by magnets with field gradients of 225-350 T m^-1and cylinder lengths of 80-110 mm to produce very narrow elongated (planar) beamlets. The corresponding dose distributions were scored in a water phantom. Composite minibeam dose distributions composed from three beamlets were created by laterally shifting copies of the single beamlet distribution to either side of a central beamlet. Modulated beamlets (with 18-30 mm nominal central SOBP) and corresponding composite dose distributions were created in a similar manner. Collimated minibeams were also compared with beams focused using one magnet/particle range combination.Main results. The focusing magnets produced planar beamlets with minimum lateral FWHM of ∼1.1-1.6 mm. Dose distributions composed from three unmodulated beamlets showed a high degree of proximal spatial fractionation and a homogeneous target dose. Maximal peak-to-valley dose ratios (PVDR) for the unmodulated beams ranged from 32 to 324, and composite modulated beam showed maximal PVDR ranging from 32 to 102 and SOBPs with good target dose coverage.Significance.Advantages of the high-gradient magnets include the ability to focus beams with phase space parameters that reflect beams in operation today, and post-waist particle divergence allowing larger beamlet separations and thus larger PVDR. Our results suggest that high gradient quadrupole magnets could be useful to focus beams of moderate emittance in clinical proton therapy.

Magnetically focused 70 MeV proton minibeams for preclinical experiments combining a tandem accelerator and a 3 GHz linear post-accelerator

PURPOSE

Radiotherapy plays an important role for the treatment of tumor diseases in two-thirds of all cases, but it is limited by side effects in the surrounding healthy tissue. Proton minibeam radiotherapy (pMBRT) is a promising option to widen the therapeutic window for tumor control at reduced side effects. An accelerator concept based on an existing tandem Van de Graaff accelerator and a linac enables the focusing of 70 MeV protons to form minibeams with a size of only 0.1 mm for a preclinical small animal irradiation facility, while avoiding the cost of an RFQ injector.

METHODS

The tandem accelerator provides a 16 MeV proton beam with a beam brightness of B = 4 nA mm 2 mrad 2 as averaged from 5 µs long pulses with a flat top current of 17 µA at 200 Hz repetition rate. Subsequently, the protons are accelerated to 70 MeV by a 3 GHz linear post-accelerator consisting of two Side Coupled Drift Tube Linac (SCDTL) structures and four Coupled Cavity Linac (CCL) structures [design: AVO-ADAM S.A (Geneva, Switzerland)]. A 3 GHz buncher and four magnetic quadrupole lenses are placed between the tandem and the post-accelerator to maximize the transmission through the linac. A quadrupole triplet situated downstream of the linac structure focuses the protons into an area of (0.1 × 0.1) mm² . The beam dynamics of the facility is optimized using the particle optics code TRACE three-dimensional (3D). Proton transmission through the facility is elaborated using the particle tracking code TRAVEL.

RESULTS

A study about buncher amplitude and phase shift between buncher and linac is showing that 49% of all protons available from the tandem can be transported through the post-accelerator. A mean beam current up to 19 nA is expected within an area of (0.1 × 0.1) mm² at the beam focus.

CONCLUSION

An extension of existing tandem accelerators by commercially available 3 GHz structures is able to deliver a proton minibeam that serves all requirements to obtain proton minibeams to perform preclinical minibeam irradiations as it would be the case for a complete commercial 3 GHz injector-RFQ-linac combination. Due to the modularity of the linac structure, the irradiation facility can be extended to clinically relevant proton energies up to or above 200 MeV.

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.

Normal Tissue Response of Combined Temporal and Spatial Fractionation in Proton Minibeam Radiation Therapy

PURPOSE

Proton minibeam radiation therapy, a spatial fractionation concept, widens the therapeutic window. By reducing normal tissue toxicities, it allows a temporally fractionated regime with high daily doses. However, an array shift between daily fractions can affect the tissue-sparing effect by decreasing the total peak-to-valley dose ratio. Therefore, combining temporal fractions with spatial fractionation raises questions about the impact of daily applied dose modulations, reirradiation accuracies, and total dose modulations.

METHODS AND MATERIALS

Healthy mouse ear pinnae were irradiated with 4 daily fractions of 30 Gy mean dose, applying proton pencil minibeams (pMB) of Gaussian σ = 222 μm in 3 different schemes: a 16 pMB array with a center-to-center distance of 1.8 mm irradiated the same position in all sessions (FS1) or was shifted by 0.9 mm to never hit the previously irradiated tissue in each session (FS2), or a 64 pMB array with a center-to-center distance of 0.9 mm irradiated the same position in all sessions (FS3), resulting in the same total dose distribution as FS2. Reirradiation positioning and its accuracy were obtained from image guidance using the unique vessel structure of ears. Acute toxicities (swelling, erythema, and desquamation) were evaluated for 153 days after the first fraction. Late toxicities (fibrous tissue, inflammation) were analyzed on day 153.

RESULTS

Reirradiation of highly dose-modulated arrays at a positioning accuracy of 110 ± 52 μm induced the least severe acute and late toxicities. A shift of the same array in FS2 led to significantly inducted acute toxicities, a higher otitis score, and a slight increase in fibrous tissue. FS3 led to the strongest increase in acute and late toxicities.

CONCLUSIONS

The highest normal-tissue sparing is achieved after accurate reirradiation of a highly dose modulated pMB array, although high positioning accuracies are challenging in a clinical environment. Nevertheless, the same integral dose applied in highly dose-modulated fractions is superior to low daily dose-modulated fractions.

Proton pencil minibeam irradiation of an in-vivo mouse ear model spares healthy tissue dependent on beam size

Proton radiotherapy using minibeams of sub-millimeter dimensions reduces side effects in comparison to conventional proton therapy due to spatial fractionation. Since the proton minibeams widen with depth, the homogeneous irradiation of a tumor can be ensured by adjusting the beam distances to tumor size and depth to maintain tumor control as in conventional proton therapy. The inherent advantages of protons in comparison to photons like a limited range that prevents a dosage of distal tissues are maintained by proton minibeams and can even be exploited for interlacing from different beam directions. A first animal study was conducted to systematically investigate and quantify the tissue-sparing effects of proton pencil minibeams as a function of beam size and dose distributions, using beam widths between σ = 95, 199, 306, 411, 561 and 883 μm (standard deviation) at a defined center-to-center beam distance (ctc) of 1.8 mm. The average dose of 60 Gy was distributed in 4x4 minibeams using 20 MeV protons (LET ~ 2.7 keV/μm). The induced radiation toxicities were measured by visible skin reactions and ear swelling for 90 days after irradiation. The largest applied beam size to ctc ratio (σ/ctc = 0.49) is similar to a homogeneous irradiation and leads to a significant 3-fold ear thickness increase compared to the control group. Erythema and desquamation was also increased significantly 3–4 weeks after irradiation. With decreasing beam sizes and thus decreasing σ/ctc, the maximum skin reactions are strongly reduced until no ear swelling or other visible skin reactions should occur for σ/ctc < 0.032 (extrapolated from data). These results demonstrate that proton pencil minibeam radiotherapy has better tissue-sparing for smaller σ/ctc, corresponding to larger peak-to-valley dose ratios PVDR, with the best effect for σ/ctc < 0.032. However, even quite large σ/ctc (e.g. σ/ctc = 0.23 or 0.31, i.e. PVDR = 10 or 2.7) show less acute side effects than a homogeneous dose distribution. This suggests that proton minibeam therapy spares healthy tissue not only in the skin but even for dose distributions appearing in deeper layers close to the tumor enhancing its benefits for clinical proton therapy.

Beam size limit for pencil minibeam radiotherapy determined from side effects in an in-vivo mouse ear model

Side effects caused by radiation are a limiting factor to the amount of dose that can be applied to a tumor volume. A novel method to reduce side effects in radiotherapy is the use of spatial fractionation, in which a pattern of sub-millimeter beams (minibeams) is applied to spare healthy tissue. In order to determine the skin reactions in dependence of single beam sizes, which are relevant for spatially fractionated radiotherapy approaches, single pencil beams of submillimeter to 6 millimeter size were applied in BALB/c mice ears at a Small Animal Radiation Research Platform (SARRP) with a plateau dose of 60 Gy. Radiation toxicities in the ears were observed for 25 days after irradiation. Severe radiation responses were found for beams ≥ 3 mm diameter. The larger the beam diameter the stronger the observed reactions. No ear swelling and barely reddening or desquamation were found for the smallest beam sizes (0.5 and 1 mm). The findings were confirmed by histological sections. Submillimeter beams are preferred in minibeam therapy to obtain optimized tissue sparing. The gradual increase of radiation toxicity with beam size shows that also larger beams are capable of healthy tissue sparing in spatial fractionation.

A Current Review of Spatial Fractionation: Back to the Future?

Spatially fractionated radiation therapy represents a significant departure from canonical thinking in radiation oncology despite having origins in the early 1900s. The original and most common implementation of spatially fractionated radiation therapy uses commercially available blocks or multileaf collimators to deliver a nonconfluent, sieve-like pattern of radiation to the target volume in a nonuniform dose distribution. Dosimetrically, this is parameterized by the ratio of the valley dose in cold spots to the peak dose in hot spots, or the valley-to-peak dose ratio. The radiobiologic mechanisms are postulated to involve radiation-induced bystander effects, microvascular alterations, and/or immunomodulation. Current indications include bulky or locally advanced disease that would not be amenable to conventional radiation or that has proved refractory to chemoradiation. Early-phase clinical trials have shown remarkable success, with some response rates >90% and minimal toxicity. This has promoted technological developments in 3-dimensional formats (LATTICE), micron-size beams (microbeam), and proton arrays. Nevertheless, more clinical and biological data are needed to specify ideal dosimetry parameters and to formulate robust clinical indications and guidelines for optimal standardized care.

Spatially fractionated proton minibeams

Extraordinary normal tissue response to highly spatially fractionated X-ray beams has been explored for over 25 years. More recently, alternative radiation sources have been developed and utilized with the aim to evoke comparable effects. These include protons, which lend themselves well for this endeavour due to their physical depth dose characteristics as well as corresponding variable biological effectiveness. This paper addresses the motivation for using protons to generate spatially fractionated beams and reviews the technological implementations and experimental results to date. This includes simulation and feasibility studies, collimation and beam characteristics, dosimetry and biological considerations as well as the results of in vivo and in vitro studies. Experimental results are emerging indicating an extraordinary normal tissue sparing effect analogous to what has been observed for synchrotron generated X-ray microbeams. The potential for translational research and feasibility of spatially modulated proton beams in clinical settings is discussed.

Challenges and Contradictions of Metal Nano-Particle Applications for Radio-Sensitivity Enhancement in Cancer Therapy

From the very beginnings of radiotherapy, a crucial question persists with how to target the radiation effectiveness into the tumor while preserving surrounding tissues as undamaged as possible. One promising approach is to selectively pre-sensitize tumor cells by metallic nanoparticles. However, though the "physics" behind nanoparticle-mediated radio-interaction has been well elaborated, practical applications in medicine remain challenging and often disappointing because of limited knowledge on biological mechanisms leading to cell damage enhancement and eventually cell death. In the present study, we analyzed the influence of different nanoparticle materials (platinum (Pt), and gold (Au)), cancer cell types (HeLa, U87, and SKBr3), and doses (up to 4 Gy) of low-Linear Energy Transfer (LET) ionizing radiation (γ- and X-rays) on the extent, complexity and reparability of radiation-induced γH2AX + 53BP1 foci, the markers of double stand breaks (DSBs). Firstly, we sensitively compared the focus presence in nuclei during a long period of time post-irradiation (24 h) in spatially (three-dimensionally, 3D) fixed cells incubated and non-incubated with Pt nanoparticles by means of high-resolution immunofluorescence confocal microscopy. The data were compared with our preliminary results obtained for Au nanoparticles and recently published results for gadolinium (Gd) nanoparticles of approximately the same size (2⁻3 nm). Next, we introduced a novel super-resolution approach-single molecule localization microscopy (SMLM)-to study the internal structure of the repair foci. In these experiments, 10 nm Au nanoparticles were used that could be also visualized by SMLM. Altogether, the data show that different nanoparticles may or may not enhance radiation damage to DNA, so multi-parameter effects have to be considered to better interpret the radiosensitization. Based on these findings, we discussed on conclusions and contradictions related to the effectiveness and presumptive mechanisms of the cell radiosensitization by nanoparticles. We also demonstrate that SMLM offers new perspectives to study internal structures of repair foci with the goal to better evaluate potential differences in DNA damage patterns.

Implementation of planar proton minibeam radiation therapy using a pencil beam scanning system: A proof of concept study

PURPOSE

Proton minibeam radiation therapy (pMBRT) is an innovative approach that combines the advantages of minibeam radiation therapy with the more precise ballistics of protons to further reduce the side effects of radiation. One of the main challenges of this approach is the generation of very narrow proton pencil beams with an adequate dose-rate to treat patients within a reasonable treatment time (several minutes) in existing clinical facilities. The aim of this study was to demonstrate the feasibility of implementing pMBRT by combining the pencil beam scanning (PBS) technique with the use of multislit collimators. This proof of concept study of pMBRT with a clinical system is intended to guide upcoming biological experiments.

METHODS

Monte Carlo simulations (TOPAS v3.1.p2) were used to design a suitable multislit collimator to implement planar pMBRT for conventional pencil beam scanning settings. Dose distributions (depth-dose curves, lateral profiles, Peak-to-Valley Dose Ratio (PVDR) and dose-rates) for different proton beam energies were assessed by means of Monte Carlo simulations and experimental measurements in a water tank using commercial ionization chambers and a new p-type silicon diode, the IBA RAZOR. An analytical intensity-modulated dose calculation algorithm designed to optimize the weight of individual Bragg peaks composing the field was also developed and validated.

RESULTS

Proton minibeams were then obtained using a brass multislit collimator with five slits measuring 2 cm × 400 μm in width with a center-to-center distance of 4 mm. The measured and calculated dose distributions (depth-dose curves and lateral profiles) showed a good agreement. Spread-out Bragg peaks (SOBP) and homogeneous dose distributions around the target were obtained by means of intensity modulation of Bragg peaks, while maintaining spatial fractionation at shallow depths. Mean dose-rates of 0.12 and 0.09 Gy/s were obtained for one iso-energy layer and a SOBP conditions in the presence of multislit collimator.

CONCLUSIONS

This study demonstrates the feasibility of implementing pMBRT on a PBS system. It also confirms the reliability of RAZOR detector for pMBRT dosimetry. This newly developed experimental methodology will support the design of future preclinical research with pMBRT.

Luminescence imaging of water during uniform-field irradiation by spot scanning proton beams

Luminescence was found during pencil-beam proton irradiation to water phantom and range could be estimated from the luminescence images. However, it is not yet clear whether the luminescence imaging is applied to the uniform fields made of spot-scanning proton-beam irradiations. For this purpose, imaging was conducted for the uniform fields having spread out Bragg peak (SOBP) made by spot scanning proton beams. We designed six types of the uniform fields with different ranges, SOBP widths and irradiation fields. One of the designed fields was irradiated to water phantom and a cooled charge coupled device camera was used to measure the luminescence image during irradiations. We estimated the ranges, field widths, and luminescence intensities from the luminescence images and compared those with the dose distribution calculated by a treatment planning system. For all types of uniform fields, we could obtain clear images of the luminescence showing the SOBPs. The ranges and field widths evaluated from the luminescence were consistent with those of the dose distribution calculated by a treatment planning system within the differences of  -4 mm and  -11 mm, respectively. Luminescence intensities were almost proportional to the SOBP widths perpendicular to the beam direction. The luminescence imaging could be applied to uniform fields made of spot scanning proton beam irradiations. Ranges and widths of the uniform fields with SOBP could be estimated from the images. The luminescence imaging is promising for the range and field width estimations in proton therapy.