Identifying optimal clinical scenarios for synchrotron microbeam radiation therapy: A treatment planning study

Inhomogeneously distributed ferroptosis with a high peak-to-valley ratio may improve the antitumor immune response

Combined radiotherapy (RT) and mild hyperthermia have been used clinically for decades to increase local control. Both modalities tend to achieve a homogeneous dose distribution within treatment targets to induce immunogenic cell death. However, marked, and long-lasting abscopal effects have not usually been observed. We proposed a hypothesis to emphasize the importance of the peak-to-valley ratio of the dose distribution inside the tumor to induce immunogenic ferrroptosis in peak area while avoid nonimmunogenic ferroptosis in valley area. Although inhomogeneous distributed energy absorption has been noted in many anticancer medical fields, the idea of sedulously created dose inhomogeneity related to antitumor immunity has not been discussed. To scale up the peak-to-valley ratio, we proposed possible implications by the combination of nanoparticles (NP) with conventional RT or hyperthermia, or the use of a high modulation depth of extremely low frequency hyperthermia or high resolution spatially fractionated radiotherapy (SFRT) to enhance the antitumor immune reactions.

Accurate and Fast Deep Learning Dose Prediction for a Preclinical Microbeam Radiation Therapy Study Using Low-Statistics Monte Carlo Simulations

Microbeam radiation therapy (MRT) utilizes coplanar synchrotron radiation beamlets and is a proposed treatment approach for several tumor diagnoses that currently have poor clinical treatment outcomes, such as gliosarcomas. Monte Carlo (MC) simulations are one of the most used methods at the Imaging and Medical Beamline, Australian Synchrotron to calculate the dose in MRT preclinical studies. The steep dose gradients associated with the 50μm-wide coplanar beamlets present a significant challenge for precise MC simulation of the dose deposition of an MRT irradiation treatment field in a short time frame. The long computation times inhibit the ability to perform dose optimization in treatment planning or apply online image-adaptive radiotherapy techniques to MRT. Much research has been conducted on fast dose estimation methods for clinically available treatments. However, such methods, including GPU Monte Carlo implementations and machine learning (ML) models, are unavailable for novel and emerging cancer radiotherapy options such as MRT. In this work, the successful application of a fast and accurate ML dose prediction model for a preclinical MRT rodent study is presented for the first time. The ML model predicts the peak doses in the path of the microbeams and the valley doses between them, delivered to the tumor target in rat patients. A CT imaging dataset is used to generate digital phantoms for each patient. Augmented variations of the digital phantoms are used to simulate with Geant4 the energy depositions of an MRT beam inside the phantoms with 15% (high-noise) and 2% (low-noise) statistical uncertainty. The high-noise MC simulation data are used to train the ML model to predict the energy depositions in the digital phantoms. The low-noise MC simulations data are used to test the predictive power of the ML model. The predictions of the ML model show an agreement within 3% with low-noise MC simulations for at least 77.6% of all predicted voxels (at least 95.9% of voxels containing tumor) in the case of the valley dose prediction and for at least 93.9% of all predicted voxels (100.0% of voxels containing tumor) in the case of the peak dose prediction. The successful use of high-noise MC simulations for the training, which are much faster to produce, accelerates the production of the training data of the ML model and encourages transfer of the ML model to different treatment modalities for other future applications in novel radiation cancer therapies.

A high resolution dose calculation engine for x-ray microbeams radiation therapy

BACKGROUND

Microbeam radiation therapy (MRT) is a treatment modality based on spatial fractionation of synchrotron generated x-rays into parallel, high dose, microbeams of a few microns width. MRT is still an under-development radiosurgery technique for which, promising preclinical results on brain tumors and epilepsy encourages its clinical transfer.

PURPOSE

A safe clinical transfer of MRT needs a specific treatment planning system (TPS) that provides accurate dose calculations in human patients, taking into account the MRT beams properties (high dose gradients, spatial fractionation, polarization effects). So far, the most advanced MRT treatment planning system, based on a hybrid dose calculation algorithm, is limited to a macroscopic rendering of the dose and does not account for the complex dose distribution inherent to MRT if delivered as conformal irradiations with multiple incidences. For overcoming these limitations, a multi-scale full Monte-Carlo calculation engine called penMRT has been developed and benchmarked against two general purpose Monte Carlo codes: penmain based on PENELOPE and Gate based on Geant4.

METHODS

PenMRT, is based on the PENELOPE (2018) Monte Carlo (MC) code, modified to take into account the voxelized geometry of the patients (CT-scans) and offering an adaptive micrometric dose calculation grid independent to the CT size, location and orientation. The implementation of the dynamic memory allocation in penMRT, makes the simulations feasible within a huge number of dose scoring bins. The possibility of using a source replication approach to simulate arrays of microbeams, and the parallelization using OpenMPI have been added to penMRT in order to increase the calculation speed for clinical usages. This engine can be implemented in a TPS as a dose calculation core.

RESULTS

The performance tests highlight the reliability of penMRT to be used for complex irradiation conditions in MRT. The benchmarking against a standard PENELOPE code did not show any significant difference for calculations in centimetric beams, for a single microbeam and for a microbeam array. The comparisons between penMRT and Gate as an independent MC code did not show any difference in the beam paths, whereas in valley regions, relative differences between the two codes rank from 1 to 7.5% which are probably due to the differences in physics lists that are used in these two codes. The reliability of the source replication approach has also been tested and validated with an underestimation of no more than 0.6% in low dose areas.

CONCLUSIONS

Good agreements (a relative difference between 0 to 8%) were found when comparing calculated peak to valley dose ratio (PVDR) values using penMRT, for irradiations with a full microbeam array, with calculated values in the literature. The high-resolution calculated dose maps obtained with penMRT are used to extract differential and cumulative dose-volume histograms (DVHs) and analyze treatment plans with much finer metrics regarding the irradiation complexity. To our knowledge, these are the first high-resolution dose maps and associated DVHs ever obtained for cross-fired microbeams irradiation, which is bringing a significant added value to the field of treatment planning in spatially fractionated radiation therapy. This article is protected by copyright. All rights reserved.

Treatment Planning Study for Microbeam Radiotherapy Using Clinical Patient Data

Microbeam radiotherapy (MRT) is a novel, still preclinical dose delivery technique. MRT has shown reduced normal tissue effects at equal tumor control rates compared to conventional radiotherapy. Treatment planning studies are required to permit clinical application. The aim of this study was to establish a dose comparison between MRT and conventional radiotherapy and to identify suitable clinical scenarios for future applications of MRT. We simulated MRT treatment scenarios for clinical patient data using an inhouse developed planning algorithm based on a hybrid Monte Carlo dose calculation and implemented the concept of equivalent uniform dose (EUD) for MRT dose evaluation. The investigated clinical scenarios comprised fractionated radiotherapy of a glioblastoma resection cavity, a lung stereotactic body radiotherapy (SBRT), palliative bone metastasis irradiation, brain metastasis radiosurgery and hypofractionated breast cancer radiotherapy. Clinically acceptable treatment plans were achieved for most analyzed parameters. Lung SBRT seemed the most challenging treatment scenario. Major limitations comprised treatment plan optimization and dose calculation considering the tissue microstructure. This study presents an important step of the development towards clinical MRT. For clinical treatment scenarios using a sophisticated dose comparison concept based on EUD and EQD2, we demonstrated the capability of MRT to achieve clinically acceptable dose distributions.

A commercial treatment planning system with a hybrid dose calculation algorithm for synchrotron radiotherapy trials

Synchrotron Radiotherapy (SyncRT) is a preclinical radiation treatment which delivers synchrotron x-rays to cancer targets. SyncRT allows for novel treatments such as Microbeam Radiotherapy, which has been shown to have exceptional healthy tissue sparing capabilities while maintaining good tumour control. Veterinary trials in SyncRT are anticipated to take place in the near future at the Australian Synchrotron's Imaging and Medical Beamline (IMBL). However, before veterinary trials can commence, a computerised treatment planning system (TPS) is required, which can quickly and accurately calculate the synchrotron x-ray dose through patient CT images. Furthermore, SyncRT TPS's must be familiar and intuitive to radiotherapy planners in order to alleviate necessary training and reduce user error. We have paired an accurate and fast Monte Carlo (MC) based SyncRT dose calculation algorithm with Eclipse^TM, the most widely implemented commercial TPS in the clinic. Using Eclipse^TM, we have performed preliminary SyncRT trials on dog cadavers at the IMBL, and verified calculated doses against dosimetric measurement to within 5% for heterogeneous tissue-equivalent phantoms. We have also performed a validation of the TPS against a full MC simulation for constructed heterogeneous phantoms in Eclipse^TM, and showed good agreement for a range of water-like tissues to within 5%-8%. Our custom Eclipse^TM TPS for SyncRT is ready to perform live veterinary trials at the IMBL.

Investigating the Dosimetric Characteristics of Microbeam Radiation Treatment

BACKGROUND

High-radiation therapeutic gain could be achieved by the modern technique of microbeam radiation treatment (MRT). The aim of this study was to investigate the dosimetric properties of MRT.

METHODS

The EGSnrc Monte Carlo (MC) code system was used to transport photons and electrons in MRT. The mono-energetic beams (1 cm × 1 cm array) of 50, 100, and 150 keV and the spectrum photon beam (European Synchrotron Radiation Facility [ESRF]) were modeled to transport through multislit collimators with the aperture's widths of 25 and 50 μm and the center-to-center (c-t-c) distance between two adjacent microbeams (MBs) of 200 μm. The calculated phase spaces at the upper surface of water phantom (1 cm × 1 cm) were implemented in DOSXYZnrc code to calculate the percentage depth dose (PDD), the dose profile curves (in depths of 0-1, 1-2, and 3-4 cm), and the peak-to-valley dose ratios (PVDRs).

RESULTS

The PDD, dose profile curves, and PVDRs were calculated for different effective parameters. The more flatness of lateral dose profile was obtained for the ESRF spectrum MB. With constant c-t-c distance, an increase in the MB size increased the peak and valley dose; simultaneously, the PVDR was larger for the 25 μm MB (33.5) compared to 50 μm MB (21.9) beam, due to the decreased scattering photons followed to the lower overlapping of the adjacent MBs. An increase in the depth decreased the PVDRs (i.e., 54.9 in depth of 0-1 cm).

CONCLUSION

Our MC model of MRT successfully calculated the effect of dosimetric parameters including photon's energy, beam width, and depth to estimate the dose distribution.

Iodine nanoparticle radiotherapy of human breast cancer growing in the brains of athymic mice

About 30% of breast cancers metastasize to the brain; those widely disseminated are fatal typically in 3-4 months, even with the best available treatments, including surgery, drugs, and radiotherapy. To address this dire situation, we have developed iodine nanoparticles (INPs) that target brain tumors after intravenous (IV) injection. The iodine then absorbs X-rays during radiotherapy (RT), creating free radicals and local tumor damage, effectively boosting the local RT dose at the tumor. Efficacy was tested using the very aggressive human triple negative breast cancer (TNBC, MDA-MB-231 cells) growing in the brains of athymic nude mice. With a well-tolerated non-toxic IV dose of the INPs (7 g iodine/kg body weight), tumors showed a heavily iodinated rim surrounding the tumor having an average uptake of 2.9% iodine by weight, with uptake peaks at 4.5%. This is calculated to provide a dose enhancement factor of approximately 5.5 (peaks at 8.0), the highest ever reported for any radiation-enhancing agents. With RT alone (15 Gy, single dose), all animals died by 72 days; INP pretreatment resulted in longer-term remissions with 40% of mice surviving 150 days and 30% surviving > 280 days.

A Monte Carlo model of synchrotron radiotherapy shows good agreement with experimental dosimetry measurements: Data from the imaging and medical beamline at the Australian Synchrotron

Experimental measurement of Synchrotron Radiotherapy (SyncRT) doses is challenging, especially for Microbeam Radiotherapy (MRT), which is characterised by very high dynamic ranges with spatial resolutions on the micrometer scale. Monte Carlo (MC) simulation is considered a gold standard for accurate dose calculation in radiotherapy, and is therefore routinely relied upon to produce verification data. We present a MC model for Australian Synchrotron's Imaging and Medical Beamline (IMBL), which is capable of generating accurate dosimetry data to inform and/or verify SyncRT experiments. Our MC model showed excellent agreement with dosimetric measurement for Synchrotron Broadbeam Radiotherapy (SBBR). Our MC model is also the first to achieve validation for MRT, using two methods of dosimetry, to within clinical tolerances of 5% for a 20×20 mm² field size, except for surface measurements at 5 mm depth, which remained to within good agreement of 7.5%. Our experimental methodology has allowed us to control measurement uncertainties for MRT doses to within 5-6%, which has also not been previously achieved, and provides a confidence which until now has been lacking in MRT validation studies. The MC model is suitable for SyncRT dose calculation of clinically relevant field sizes at the IMBL, and can be extended to include medical beamlines at other Synchrotron facilities as well. The presented MC model will be used as a validation tool for treatment planning dose calculation algorithms, and is an important step towards veterinary SyncRT trials at the Australian Synchrotron.

Dose profiles and x-ray energy optimization for microbeam radiation therapy by high-dose, high resolution dosimetry using Sm-doped fluoroaluminate glass plates and Monte Carlo transport simulation

Microbeam radiation therapy (MRT) utilizes highly collimated synchrotron generated x-rays to create narrow planes of high dose radiation for the treatment of tumors. Individual microbeams have a typical width of 30-50 µm and are separated by a distance of 200-500 µm. The dose delivered at the center of the beam is lethal to cells in the microbeam path, on the order of hundreds of Grays (Gy). The tissue between each microbeam is spared and helps aid in the repair of adjacent damaged tissue. Radiation interactions within the peak of the microbeam, such as the photoelectric effect and incoherent (atomic Compton) scattering, cause some dose to be delivered to the valley areas adjacent to the microbeams. As the incident x-ray energy is modified, radiation interactions within a material change and affect the probability of interactions, as well as the directionality and energy of ionizing particles (electrons) that deposit energy in the valley regions surrounding the microbeam peaks. It is crucial that the valley dose between microbeams be minimal to maintain the effectiveness of MRT. Using a monochromatic x-ray source with x-ray energies ranging from 30 to 150 keV, a detailed investigation into the effect of incident x-ray energy on the dose profiles of microbeams was performed using samarium doped fluoroaluminate (FA) glass as the medium. All dosimetric measurements were carried out using a purpose-built fluorescence confocal microscope dosimetric technique that used Sm-doped FA glass plates as the irradiated medium. Dose profiles are measured over a very a wide range of x-ray energies at micrometer resolution and dose distribution in the microbeam are mapped. The measured microbeam profiles at different energies are compared with the MCNP6 radiation transport code, a general transport code which can calculate the energy deposition of electrons as they pass through a given material. The experimentally measured distributions can be used to validate the results for electron energy deposition in fluoroaluminate glass. Code validation is necessary for using transport codes in future treatment planning for MRT and other radiation therapies. It is shown that simulated and measured micro beam-profiles are in good agreement, and micrometer level changes can be observed using this high-resolution dosimetry technique. Full width at 10% of the maximum peak (FW@10%) was used to quantify the microbeam width. Experimental measurements on FA glasses and simulations on the dependence of the FW@10% at various energies are in good agreement. Simulations on energy deposited in water indicate that FW@10% reaches a local minimum around energies 140 keV. In addition, variable slit width experiments were carried out at an incident x-ray energy of 100 keV in order to determine the effect of the narrowing slit width on the delivered peak dose. The microbeam width affects the peak dose, which decreases with the width of the microbeam. Experiments suggest that a typical microbeam width for MRT is likely to be between 20-50 µm based on this work.

Animal Models in Microbeam Radiation Therapy: A Scoping Review

BACKGROUND

Microbeam Radiation Therapy (MRT) is an innovative approach in radiation oncology where a collimator subdivides the homogeneous radiation field into an array of co-planar, high-dose beams which are tens of micrometres wide and separated by a few hundred micrometres.

OBJECTIVE

This scoping review was conducted to map the available evidence and provide a comprehensive overview of the similarities, differences, and outcomes of all experiments that have employed animal models in MRT.

METHODS

We considered articles that employed animal models for the purpose of studying the effects of MRT. We searched in seven databases for published and unpublished literature. Two independent reviewers screened citations for inclusion. Data extraction was done by three reviewers.

RESULTS

After screening 5688 citations and 159 full-text papers, 95 articles were included, of which 72 were experimental articles. Here we present the animal models and pre-clinical radiation parameters employed in the existing MRT literature according to their use in cancer treatment, non-neoplastic diseases, or normal tissue studies.

CONCLUSIONS

The study of MRT is concentrated in brain-related diseases performed mostly in rat models. An appropriate comparison between MRT and conventional radiotherapy (instead of synchrotron broad beam) is needed. Recommendations are provided for future studies involving MRT.