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.

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.

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.

Quantitative evaluation of the effect of changes in effective energy on the image quality in X-ray computed tomography

The X-ray effective energy differs for each computed tomography (CT) scanner even at the same tube voltage because of differences in the bow-tie filter and additional filter. Even when scanning with the same tube voltage and dose setting, these differences in effective energy result in different image noise levels. Although this qualitative change is known, the related quantitative changes have not been clarified. In this study, using two CT scanners with the same geometric specifications and detector configurations, we quantitatively assessed the reduction in image noise accompanying the increase in effective energy. We also clarified the fluctuations in CT number. For both CT scanners, the effective energy, the standard deviation (SD) of the noise image when using two water phantoms with diameters of 240 mm and 320 mm, and CT numbers of the sensitometry module were measured. Further, the dose required to obtain the same image noise level in each CT scanner was calculated. The effective energy difference was 5.5 keV to 10.7 keV, and the difference tended to be larger when the scan field of view was larger. The SD differences were 24% and 14% for the 320-mm and 240-mm phantoms, respectively. For converting to the dose required to obtain the same SD, the dose can be reduced by 42% and 24% for the 320-mm and 240-mm phantoms, respectively. The CT number difference of both CT scanners was small. Therefore, higher effective energy contributes to the reduction of image noise.

Technical advances in x-ray microbeam radiation therapy

In the last 25 years microbeam radiation therapy (MRT) has emerged as a promising alternative to conventional radiation therapy at large, third generation synchrotrons. In MRT, a multi-slit collimator modulates a kilovoltage x-ray beam on a micrometer scale, creating peak dose areas with unconventionally high doses of several hundred Grays separated by low dose valley regions, where the dose remains well below the tissue tolerance level. Pre-clinical evidence demonstrates that such beam geometries lead to substantially reduced damage to normal tissue at equal tumour control rates and hence drastically increase the therapeutic window. Although the mechanisms behind MRT are still to be elucidated, previous studies indicate that immune response, tumour microenvironment, and the microvasculature may play a crucial role. Beyond tumour therapy, MRT has also been suggested as a microsurgical tool in neurological disorders and as a primer for drug delivery. The physical properties of MRT demand innovative medical physics and engineering solutions for safe treatment delivery. This article reviews technical developments in MRT and discusses existing solutions for dosimetric validation, reliable treatment planning and safety. Instrumentation at synchrotron facilities, including beam production, collimators and patient positioning systems, is also discussed. Specific solutions reviewed in this article include: dosimetry techniques that can cope with high spatial resolution, low photon energies and extremely high dose rates of up to 15 000 Gy s-1, dose calculation algorithms-apart from pure Monte Carlo Simulations-to overcome the challenge of small voxel sizes and a wide dynamic dose-range, and the use of dose-enhancing nanoparticles to combat the limited penetrability of a kilovoltage energy spectrum. Finally, concepts for alternative compact microbeam sources are presented, such as inverse Compton scattering set-ups and carbon nanotube x-ray tubes, that may facilitate the transfer of MRT into a hospital-based clinical environment. Intensive research in recent years has resulted in practical solutions to most of the technical challenges in MRT. Treatment planning, dosimetry and patient safety systems at synchrotrons have matured to a point that first veterinary and clinical studies in MRT are within reach. Should these studies confirm the promising results of pre-clinical studies, the authors are confident that MRT will become an effective new radiotherapy option for certain patients.

Validation of a Monte Carlo simulation for Microbeam Radiation Therapy on the Imaging and Medical Beamline at the Australian Synchrotron

Microbeam Radiation Therapy (MRT) is an emerging cancer treatment modality characterised by the use of high-intensity synchrotron-generated x-rays, spatially fractionated by a multi-slit collimator (MSC), to ablate target tumours. The implementation of an accurate treatment planning system, coupled with simulation tools that allow for independent verification of calculated dose distributions are required to ensure optimal treatment outcomes via reliable dose delivery. In this article we present data from the first Geant4 Monte Carlo radiation transport model of the Imaging and Medical Beamline at the Australian Synchrotron. We have developed the model for use as an independent verification tool for experiments in one of three MRT delivery rooms and therefore compare simulation results with equivalent experimental data. The normalised x-ray spectra produced by the Geant4 model and a previously validated analytical model, SPEC, showed very good agreement using wiggler magnetic field strengths of 2 and 3 T. However, the validity of absolute photon flux at the plane of the Phase Space File (PSF) for a fixed number of simulated electrons was unable to be established. This work shows a possible limitation of the G4SynchrotronRadiation process to model synchrotron radiation when using a variable magnetic field. To account for this limitation, experimentally derived normalisation factors for each wiggler field strength determined under reference conditions were implemented. Experimentally measured broadbeam and microbeam dose distributions within a Gammex RMI457 Solid Water® phantom were compared to simulated distributions generated by the Geant4 model. Simulated and measured broadbeam dose distributions agreed within 3% for all investigated configurations and measured depths. Agreement between the simulated and measured microbeam dose distributions agreed within 5% for all investigated configurations and measured depths.

Ultra high dose rate Synchrotron Microbeam Radiation Therapy. Preclinical evidence in view of a clinical transfer

This paper reviews the current state of the art of an emerging form of radiosurgery dedicated to brain tumour treatment and which operates at very high dose rate (kGy·s^-1). Microbeam Radiation Therapy uses synchrotron-generated X-rays which triggered normal tissue sparing partially mediated by FLASH effect.

Hybrid dose calculation: a dose calculation algorithm for microbeam radiation therapy

Microbeam radiation therapy (MRT) is still a preclinical approach in radiation oncology that uses planar micrometre wide beamlets with extremely high peak doses, separated by a few hundred micrometre wide low dose regions. Abundant preclinical evidence demonstrates that MRT spares normal tissue more effectively than conventional radiation therapy, at equivalent tumour control. In order to launch first clinical trials, accurate and efficient dose calculation methods are an inevitable prerequisite. In this work a hybrid dose calculation approach is presented that is based on a combination of Monte Carlo and kernel based dose calculation. In various examples the performance of the algorithm is compared to purely Monte Carlo and purely kernel based dose calculations. The accuracy of the developed algorithm is comparable to conventional pure Monte Carlo calculations. In particular for inhomogeneous materials the hybrid dose calculation algorithm out-performs purely convolution based dose calculation approaches. It is demonstrated that the hybrid algorithm can efficiently calculate even complicated pencil beam and cross firing beam geometries. The required calculation times are substantially lower than for pure Monte Carlo calculations.

A point kernel algorithm for microbeam radiation therapy

Microbeam radiation therapy (MRT) is a treatment approach in radiation therapy where the treatment field is spatially fractionated into arrays of a few tens of micrometre wide planar beams of unusually high peak doses separated by low dose regions of several hundred micrometre width. In preclinical studies, this treatment approach has proven to spare normal tissue more effectively than conventional radiation therapy, while being equally efficient in tumour control. So far dose calculations in MRT, a prerequisite for future clinical applications are based on Monte Carlo simulations. However, they are computationally expensive, since scoring volumes have to be small. In this article a kernel based dose calculation algorithm is presented that splits the calculation into photon and electron mediated energy transport, and performs the calculation of peak and valley doses in typical MRT treatment fields within a few minutes. Kernels are analytically calculated depending on the energy spectrum and material composition. In various homogeneous materials peak, valley doses and microbeam profiles are calculated and compared to Monte Carlo simulations. For a microbeam exposure of an anthropomorphic head phantom calculated dose values are compared to measurements and Monte Carlo calculations. Except for regions close to material interfaces calculated peak dose values match Monte Carlo results within 4% and valley dose values within 8% deviation. No significant differences are observed between profiles calculated by the kernel algorithm and Monte Carlo simulations. Measurements in the head phantom agree within 4% in the peak and within 10% in the valley region. The presented algorithm is attached to the treatment planning platform VIRTUOS. It was and is used for dose calculations in preclinical and pet-clinical trials at the biomedical beamline ID17 of the European synchrotron radiation facility in Grenoble, France.

Microbeam radiation therapy - grid therapy and beyond: a clinical perspective

Microbeam irradiation is spatially fractionated radiation on a micrometer scale. Microbeam irradiation with therapeutic intent has become known as microbeam radiation therapy (MRT). The basic concept of MRT was developed in the 1980s, but it has not yet been tested in any human clinical trial, even though there is now a large number of animal studies demonstrating its marked therapeutic potential with an exceptional normal tissue sparing effect. Furthermore, MRT is conceptually similar to macroscopic grid based radiation therapy which has been used in clinical practice for decades. In this review, the potential clinical applications of MRT are analysed for both malignant and non-malignant diseases.

Preclinical radiotherapy at the Australian Synchrotron's Imaging and Medical Beamline: instrumentation, dosimetry and a small-animal feasibility study

Therapeutic applications of synchrotron X-rays such as microbeam (MRT) and minibeam (MBRT) radiation therapy promise significant advantages over conventional clinical techniques for some diseases if successfully transferred to clinical practice. Preclinical studies show clear evidence that a number of normal tissues in animal models display a tolerance to much higher doses from MRT compared with conventional radiotherapy. However, a wide spread in the parameters studied makes it difficult to make any conclusions about the associated tumour control or normal tissue complication probabilities. To facilitate more systematic and reproducible preclinical synchrotron radiotherapy studies, a dedicated preclinical station including small-animal irradiation stage was designed and installed at the Imaging and Medical Beamline (IMBL) at the Australian Synchrotron. The stage was characterized in terms of the accuracy and reliability of the vertical scanning speed, as this is the key variable in dose delivery. The measured speed was found to be within 1% of the nominal speed for the range of speeds measured by an interferometer. Furthermore, dose measurements confirm the expected relationship between speed and dose and show that the measured dose is independent of the scan direction. Important dosimetric parameters such as peak dose, valley dose, the collimator output factor and peak-to-valley dose ratio are presented for 5 mm × 5 mm, 10 mm × 10 mm and 20 mm × 20 mm field sizes. Finally, a feasibility study on three glioma-bearing rats was performed. MRT and MBRT doses were prescribed to achieve an average dose of 65 Gy in the target, and magnetic resonance imaging follow-up was performed at various time points after irradiation to follow the tumour volume. Although it is impossible to draw conclusions on the different treatments with such a small number of animals, the feasibility of end-to-end preclinical synchrotron radiotherapy studies using the IMBL preclinical stage is demonstrated.

A kernel-based dose calculation algorithm for kV photon beams with explicit handling of energy and material dependencies

Objective:

Mimicking state-of-the-art patient radiotherapy with high-precision irradiators for small animals is expected to advance the understanding of dose–effect relationships and radiobiology in general. We work on the implementation of intensity-modulated radiotherapy-like irradiation schemes for small animals. As a first step, we present a fast analytical dose calculation algorithm for keV photon beams.

Methods:

We follow a superposition–convolution approach adapted to kV X-rays, based on previous work for microbeam therapy. We assume local energy deposition at the photon interaction point due to the short electron ranges in tissue. This allows us to separate the dose calculation into locally absorbed primary dose and the scatter contribution, calculated in a point kernel approach. We validate our dose model against Geant4 Monte Carlo (MC) simulations and compare the results to Muriplan (XStrahl Ltd, Camberley, UK).

Results:

For field sizes of (1 mm)² to (1 cm)² in water, the depth dose curves show a mean disagreement of 1.7% to MC simulations, with the largest deviations in the entrance region (4%) and at large depths (5% at 7 cm). Larger discrepancies are observed at water-to-bone boundaries, in bone and at the beam edges in slab phantoms and a mouse brain. Calculation times are in the order of 5 s for a single beam.

Conclusion:

The algorithm shows good agreement with MC simulations in an initial validation. It has the potential to become an alternative to full MC dose calculation.

Advances in knowledge:

The presented algorithm demonstrates the potential of kernel-based dose calculation for kV photon beams. It will be valuable in intensity-modulated radiotherapy and inverse treatment planning for high precision small-animal radiotherapy.

Medical physics aspects of the synchrotron radiation therapies: Microbeam radiation therapy (MRT) and synchrotron stereotactic radiotherapy (SSRT)

Stereotactic Synchrotron Radiotherapy (SSRT) and Microbeam Radiation Therapy (MRT) are both novel approaches to treat brain tumor and potentially other tumors using synchrotron radiation. Although the techniques differ by their principles, SSRT and MRT share certain common aspects with the possibility of combining their advantages in the future. For MRT, the technique uses highly collimated, quasi-parallel arrays of X-ray microbeams between 50 and 600 keV. Important features of highly brilliant Synchrotron sources are a very small beam divergence and an extremely high dose rate. The minimal beam divergence allows the insertion of so called Multi Slit Collimators (MSC) to produce spatially fractionated beams of typically ∼25-75 micron-wide microplanar beams separated by wider (100-400 microns center-to-center(ctc)) spaces with a very sharp penumbra. Peak entrance doses of several hundreds of Gy are extremely well tolerated by normal tissues and at the same time provide a higher therapeutic index for various tumor models in rodents. The hypothesis of a selective radio-vulnerability of the tumor vasculature versus normal blood vessels by MRT was recently more solidified. SSRT (Synchrotron Stereotactic Radiotherapy) is based on a local drug uptake of high-Z elements in tumors followed by stereotactic irradiation with 80 keV photons to enhance the dose deposition only within the tumor. With SSRT already in its clinical trial stage at the ESRF, most medical physics problems are already solved and the implemented solutions are briefly described, while the medical physics aspects in MRT will be discussed in more detail in this paper.

Energy spectra considerations for synchrotron radiotherapy trials on the ID17 bio-medical beamline at the European Synchrotron Radiation Facility

The aim of this study was to validate the kilovoltage X-ray energy spectrum on the ID17 beamline at the European Synchrotron Radiation Facility (ESRF). The purpose of such validation was to provide an accurate energy spectrum as the input to a computerized treatment planning system, which will be used in synchrotron microbeam radiotherapy trials at the ESRF. Calculated and measured energy spectra on ID17 have been reported previously but recent additions and safety modifications to the beamline for veterinary trials warranted a fresh investigation. The authors used an established methodology to compare X-ray attenuation measurements in copper sheets (referred to as half value layer measurements in the radiotherapy field) with the predictions of a theoretical model. A cylindrical ionization chamber in air was used to record the relative attenuation of the X-ray beam intensity by increasing thicknesses of high-purity copper sheets. The authors measured the half value layers in copper for two beamline configurations, which corresponded to differing spectral conditions. The authors obtained good agreement between the measured and predicted half value layers for the two beamline configurations. The measured first half value layer was 1.754 ± 0.035 mm Cu and 1.962 ± 0.039 mm Cu for the two spectral conditions, compared with theoretical predictions of 1.763 ± 0.039 mm Cu and 1.984 ± 0.044 mm Cu, respectively. The calculated mean energies for the two conditions were 105 keV and 110 keV and there was not a substantial difference in the calculated percentage depth dose curves in water between the different spectral conditions. The authors observed a difference between their calculated energy spectra and the spectra previously reported by other authors, particularly at energies greater than 100 keV. The validation of the beam spectrum by the copper half value layer measurements means the authors can provide an accurate spectrum as an input to a treatment planning system for the forthcoming veterinary trials of microbeam radiotherapy to spontaneous tumours in cats and dogs.

Modeling a superficial radiotherapy X-ray source for relative dose calculations

The purpose of this study was to empirically characterize and validate a kilovoltage (kV) X‐ray beam source model of a superficial X‐ray unit for relative dose calculations in water and assess the accuracy of the British Journal of Radiology Supplement 25 (BJR 25) percentage depth dose (PDD) data. We measured central axis PDDs and dose profiles using an Xstrahl 150 X‐ray system. We also compared the measured and calculated PDDs to those in the BJR 25. The Xstrahl source was modeled as an effective point source with varying spatial fluence and spectra. In‐air ionization chamber measurements were made along the x‐ and y‐axes of the X‐ray beam to derive the spatial fluence and half‐value layer (HVL) measurements were made to derive the spatially varying spectra. This beam characterization and resulting source model was used as input for our in‐house dose calculation software (kVDoseCalc) to compute radiation dose at points of interest (POIs). The PDDs and dose profiles were measured using 2, 5, and 15 cm cone sizes at 80, 120, 140, and 150 kVp energies in a scanning water phantom using IBA Farmer‐type ionization chambers of volumes 0.65 and 0.13 cc, respectively. The percent difference in the computed PDDs compared with our measurements range from −4.8% to 4.8%, with an overall mean percent difference and standard deviation of 1.5% and 0.7%, respectively. The percent difference between our PDD measurements and those from BJR 25 range from −14.0% to 15.7%, with an overall mean percent difference and standard deviation of 4.9% and 2.1%, respectively — showing that the measurements are in much better agreement with kVDoseCalc than BJR 25. The range in percent difference between kVDoseCalc and measurement for profiles was −5.9% to 5.9%, with an overall mean percent difference and standard deviation of 1.4% and 1.4%, respectively. The results demonstrate that our empirically based X‐ray source modeling approach for superficial X‐ray therapy can be used to accurately compute relative dose in a homogeneous water‐equivalent medium. They also show limitations in the accuracy of the BJR 25 PDD data.

PACS number: 87.55.kh

Use of synchrotron medical microbeam irradiation to investigate radiation-induced bystander and abscopal effects in vivo

The question of whether bystander and abscopal effects are the same is unclear. Our experimental system enables us to address this question by allowing irradiated organisms to partner with unexposed individuals. Organs from both animals and appropriate sham and scatter dose controls are tested for expression of several endpoints such as calcium flux, role of 5HT, reporter assay cell death and proteomic profile. The results show that membrane related functions of calcium and 5HT are critical for true bystander effect expression. Our original inter-animal experiments used fish species whole body irradiated with low doses of X-rays, which prevented us from addressing the abscopal effect question. Data which are much more relevant in radiotherapy are now available for rats which received high dose local irradiation to the implanted right brain glioma. The data were generated using quasi-parallel microbeams at the biomedical beamline at the European Synchrotron Radiation Facility in Grenoble France. This means we can directly compare abscopal and "true" bystander effects in a rodent tumour model. Analysis of right brain hemisphere, left brain and urinary bladder in the directly irradiated animals and their unirradiated partners strongly suggests that bystander effects (in partner animals) are not the same as abscopal effects (in the irradiated animal). Furthermore, the presence of a tumour in the right brain alters the magnitude of both abscopal and bystander effects in the tissues from the directly irradiated animal and in the unirradiated partners which did not contain tumours, meaning the type of signal was different.