Monte Carlo dosimetry for forthcoming clinical trials in x-ray microbeam radiation therapy

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.

Tolerance of Normal Rabbit Facial Bones and Teeth to Synchrotron X-Ray Microbeam Irradiation

Microbeam radiation therapy, an alternative radiosurgical treatment under preclinical investigation, aims to safely treat muzzle tumors in pet animals. This will require data on the largely unknown radiation toxicity of microbeam arrays for bones and teeth. To this end, the muzzle of six young adult New Zealand rabbits was irradiated by a lateral array of microplanar beamlets with peak entrance doses of 200, 330 or 500 Gy. The muzzles were examined 431 days postirradiation by computed microtomographic imaging (micro-CT) ex vivo, and extensive histopathology. The boundaries of the radiation field were identified histologically by microbeam tracks in cartilage and other tissues. There was no radionecrosis of facial bones in any rabbit. Conversely, normal incisor teeth exposed to peak entrance doses of 330 Gy or 500 Gy developed marked caries-like damage, whereas the incisors of the two rabbits exposed to 200 Gy remained unscathed. A single, unidirectional array of microbeams with a peak entrance dose ≤200 Gy (valley dose14 Gy) did not damage normal bone, teeth and soft tissues of the muzzle of normal rabbits longer than one year after irradiation. Because of that, Microbeam radiation therapy of muzzle tumors in pet animals is unlikely to cause sizeable damage to normal teeth, bone and soft tissues, if a single array as used here delivers a limited entrance dose of 200 Gy and a valley dose of ≤14 Gy.

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.

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.

Experimental optimisation of the X-ray energy in microbeam radiation therapy

Microbeam radiation therapy has demonstrated superior normal tissue sparing properties compared to broadbeam radiation fields. The ratio of the microbeam peak dose to the valley dose (PVDR), which is dependent on the X-ray energy/spectrum and geometry, should be maximised for an optimal therapeutic ratio. Simulation studies in the literature report the optimal energy for MRT based on the PVDR. However, most of these studies have considered different microbeam geometries to that at the Imaging and Medical Beamline (50 μm beam width with a spacing of 400 μm). We present the first fully experimental investigation of the energy dependence of PVDR and microbeam penumbra. Using monochromatic X-ray energies in the range 40-120 keV the PVDR was shown to increase with increasing energy up to 100 keV before plateauing. PVDRs measured for pink beams were consistently higher than those for monochromatic energies similar or equivalent to the average energy of the spectrum. The highest PVDR was found for a pink beam average energy of 124 keV. Conversely, the microbeam penumbra decreased with increasing energy before plateauing for energies above 90 keV. The effect of bone on the PVDR was investigated at energies 60, 95 and 120 keV. At depths greater than 20 mm beyond the bone/water interface there was almost no effect on the PVDR. In conclusion, the optimal energy range for MRT at IMBL is 90-120 keV, however when considering the IMBL flux at different energies, a spectrum with 95 keV weighted average energy was found to be the best compromise.

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.

Carbon and oxygen minibeam radiation therapy: An experimental dosimetric evaluation

PURPOSE

To perform dosimetric characterization of a minibeam collimator in both carbon and oxygen ion beams to guide optimal setup geometry and irradiation for future radiobiological studies.

METHODS

Carbon and oxygen minibeams were generated using a prototype tungsten multislit collimator presenting line apertures 700 μm wide, which are spaced 3500 μm centre-to-centre distance apart. Several radiation beam spots generated the desired field size of 15 × 15 mm² and production of a 50 mm long spread out Bragg peak (SOBP) centered at 80 mm depth in water. Dose evaluations were performed with two different detectors: a PTW microDiamond® single crystal diamond detector and radiochromic films (EBT3). Peak-to-valley dose ratio (PVDR) values, output factors (OF), penumbras, and full width at half maximum (FWHM) were measured.

RESULTS

Measured lateral dose profiles exhibited spatial fractionation of dose at depth in a water phantom in the expected form of peaks and valleys for both carbon and oxygen radiation fields. The diamond detector and radiochromic film provided measurements of PVDR in good agreement. PVDR values at shallow depth were about 60 and decreased to about 10 at 80 mm depth in water. OF in the center of the SOBP was about 0.4; this value is larger than the corresponding one in proton minibeam radiation therapy measured using a comparable collimator due to a reduced lateral scattering for carbon and oxygen minibeams.

CONCLUSIONS

Carbon and oxygen minibeams may be produced by a mechanical collimator. PVDR values and output factors measured in this first study of these minibeam radiation types indicate there is potential for their therapeutic use. Optimization of minibeam collimator design and the number and size of focal spots for irradiation are advocated to improve PDVR values and dose distributions for each specific applied use.

Proton minibeam radiation therapy: Experimental dosimetry evaluation

PURPOSE

Proton minibeam radiation therapy (pMBRT) is a new radiotherapy (RT) approach that allies the inherent physical advantages of protons with the normal tissue preservation observed when irradiated with submillimetric spatially fractionated beams. This dosimetry work aims at demonstrating the feasibility of the technical implementation of pMBRT. This has been performed at the Institut Curie - Proton Therapy Center in Orsay.

METHODS

Proton minibeams (400 and 700 μm-width) were generated by means of a brass multislit collimator. Center-to-center distances between consecutive beams of 3200 and 3500 μm, respectively, were employed. The (passive scattered) beam energy was 100 MeV corresponding to a range of 7.7 cm water equivalent. Absolute dosimetry was performed with a thimble ionization chamber (IBA CC13) in a water tank. Relative dosimetry was carried out irradiating radiochromic films interspersed in a IBA RW3 slab phantom. Depth dose curves and lateral profiles at different depths were evaluated. Peak-to-valley dose ratios (PVDR), beam widths, and output factors were also assessed as a function of depth.

RESULTS

A pattern of peaks and valleys was maintained in the transverse direction with PVDR values decreasing as a function of depth until 6.7 cm. From that depth, the transverse dose profiles became homogeneous due to multiple Coulomb scattering. Peak-to-valley dose ratio values extended from 8.2 ± 0.5 at the phantom surface to 1.08 ± 0.06 at the Bragg peak. This was the first time that dosimetry in such small proton field sizes was performed. Despite the challenge, a complete set of dosimetric data needed to guide the first biological experiments was achieved.

CONCLUSIONS

pMBRT is a novel strategy in order to reduce the side effects of RT. This works provides the experimental proof of concept of this new RT method: clinical proton beams might allow depositing a (high) uniform dose in a brain tumor located in the center of the brain (7.5 cm depth, the worst scenario), while a spatial fractionation of the dose is retained in the normal tissues in the beam path, potentially leading to a gain in tissue sparing. This is the first complete experimental implementation of this promising technique. Biological experiments are needed in order to confirm the clinical potential of pMBRT.

Dosimetric evaluation of new approaches in GRID therapy using nonconventional radiation sources

PURPOSE

Spatial fractionation of the dose has proven to be a promising approach to increase the tolerance of healthy tissue, which is the main limitation of radiotherapy. A good example of that is GRID therapy, which has been successfully used in the management of large tumors with low toxicity. The aim of this work is to explore new avenues using nonconventional sources: GRID therapy by using kilovoltage (synchrotron) x-rays, the use of very high-energy electrons, and proton GRID therapy. They share in common the use of the smallest possible grid sizes in order to exploit the dose-volume effects.

METHODS

Monte Carlo simulations (penelope/peneasy and geant4/GATE codes) were used as a method to study dose distributions resulting from irradiations in different configurations of the three proposed techniques. As figure of merit, percentage (peak and valley) depth dose curves, penumbras, and central peak-to-valley dose ratios (PVDR) were evaluated. As shown in previous biological experiments, high PVDR values are requested for healthy tissue sparing. A superior tumor control may benefit from a lower PVDR.

RESULTS

High PVDR values were obtained in the healthy tissue for the three cases studied. When low energy photons are used, the treatment of deep-seated tumors can still be performed with submillimetric grid sizes. Superior PVDR values were reached with the other two approaches in the first centimeters along the beam path. The use of protons has the advantage of delivering a uniform dose distribution in the tumor, while healthy tissue benefits from the spatial fractionation of the dose. In the three evaluated techniques, there is a net reduction in penumbra with respect to radiosurgery.

CONCLUSIONS

The high PVDR values in the healthy tissue and the use of small grid sizes in the three presented approaches might constitute a promising alternative to treat tumors with such spatially fractionated radiotherapy techniques. The dosimetric results presented here support the interest of performing radiobiology experiments in order to evaluate these new avenues.

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.

Effects of high-dose microbeam irradiation on tumor microvascular function and angiogenesis

Microbeam radiation therapy (MRT) is a form of cancer treatment in which a single large dose of radiation is spatially fractionated in-line or grid-like patterns. Preclinical studies have demonstrated that MRT is capable of eliciting high levels of tumor response while sparing normal tissue that is exposed to the same radiation field. Since a large fraction of the MRT-treated tumor is in the dose valley region that is not directly irradiated, tumor response may be driven by radiation bystander effects, which in turn elicit a microvascular response. Differential alterations in hemodynamics between the tumor and normal tissue may explain the therapeutic advantages of MRT. Direct observation of these dynamic responses presents a challenge for conventional ex vivo analysis. Furthermore, knowledge gleaned from in vitro studies of radiation bystander response has not been widely incorporated into in vivo models of tumor radiotherapy, and the biological contribution of the bystander effect within the tumor microenvironment is unknown. In this study, we employed noninvasive, serial observations of the tumor microenvironment to address the question of how tumor vasculature and HIF-1 expression are affected by microbeam radiotherapy. Tumors (approximately 4 mm in diameter) grown in a dorsal window chamber were irradiated in a single fraction using either a single, microplanar beam (300 micron wide swath) or a wide-field setup (whole-window chamber) to a total dose of 50 Gy. The tumors were optically observed daily for seven days postirradiation. Microvascular changes in the tumor and surrounding normal tissue differed greatly between the wide-field and microbeam treatments. We present evidence that these changes may be due to dissimilar spatial and temporal patterns of HIF-1 expression induced through radiation bystander effects.

Double-strand breaks on F98 glioma rat cells induced by minibeam and broad-beam synchrotron radiation therapy

PURPOSE

To assess the DNA damage induced by MBRT and BB radiations on glioma cells.

METHODS

The analysis of fluorescent intensity emitted per nucleus was plotted versus DNA content 2 and 17 h after irradiations. At around cell-doubling time (17 h) after exposures, the remaining DNA radiation damage could be correlated with cellular death.

RESULTS

A higher γH2AX IF intensity per cell could be detected 2 and 17 h after MBRT when compared with BB. 17 h after MBRT, misrepaired damaged cells remained arrested in both G1 and G2 phases.

CONCLUSIONS

A pronounced G2 phase arrest was detected at 17 h after MBRT and BB. However, only after MBRT, a dose-dependent increasing number of damaged cells appeared arrested also in the G1 phase, and a higher amount of cells more prone to undergo apoptosis were detected. The threshold dose required to enhance the effectiveness of both synchrotron radiation techniques was 12 Gy.

Proton-minibeam radiation therapy: a proof of concept

PURPOSE

This Monte Carlo simulation work aims at studying a new radiotherapy approach called proton-minibeam radiation therapy (pMBRT). The main objective of this proof of concept was the evaluation of the possible gain in tissue sparing, thanks to the spatial fractionation of the dose, which could be used to deposit higher and potentially curative doses in clinical cases where tissue tolerances are a limit for conventional methods.

METHODS

Monte Carlo simulations (GATE v.6) have been used as a method to calculate the ratio of the peak-to-valley doses (PVDR) for arrays of proton minibeams of 0.7 mm width and several center-to-center distances, at different depths in a water phantom. The beam penumbras were also evaluated as an important parameter for tissue sparing, for example, in the treatment of non-cancer diseases like epilepsy. Two proton energies were considered in this study: a clinically relevant energy (105 MeV) and a very high energy (1 GeV), to benefit from a reduced lateral scattering. For the latter case, an interlaced geometry was also evaluated.

RESULTS

Higher or similar PVDR than the ones obtained in x-rays minibeam radiation therapy were achieved in several pMBRT configurations. In addition, for the two energies studied, the beam penumbras are smaller than in the case of Gamma Knife radiosurgery.

CONCLUSIONS

The high PVDR obtained for some configurations and the small penumbras in comparison with existing radiosurgery techniques, suggest a potential gain in healthy tissue sparing in this new technique. Biological studies are warranted to assess the effects of pMBRT on both normal and tumoral tissues.

Monte Carlo-based dose calculation engine for minibeam radiation therapy

Minibeam radiation therapy (MBRT) is an innovative radiotherapy approach based on the well-established tissue sparing effect of arrays of quasi-parallel micrometre-sized beams. In order to guide the preclinical trials in progress at the European Synchrotron Radiation Facility (ESRF), a Monte Carlo-based dose calculation engine has been developed and successfully benchmarked with experimental data in anthropomorphic phantoms. Additionally, a realistic example of treatment plan is presented. Despite the micron scale of the voxels used to tally dose distributions in MBRT, the combination of several efficiency optimisation methods allowed to achieve acceptable computation times for clinical settings (approximately 2 h). The calculation engine can be easily adapted with little or no programming effort to other synchrotron sources or for dose calculations in presence of contrast agents.

Synchrotron-generated microbeam sensorimotor cortex transections induce seizure control without disruption of neurological functions

Synchrotron-generated X-ray microplanar beams (microbeams) are characterized by the ability to deliver extremely high doses of radiation to spatially restricted volumes of tissue. Minimal dose spreading outside the beam path provides an exceptional degree of protection from radio-induced damage to the neurons and glia adjacent to the microscopic slices of tissue irradiated. The preservation of cortical architecture following high-dose microbeam irradiation and the ability to induce non-invasively the equivalent of a surgical cut over the cortex is of great interest for the development of novel experimental models in neurobiology and new treatment avenues for a variety of brain disorders. Microbeams (size 100 µm/600 µm, center-to-center distance of 400 µm/1200 µm, peak entrance doses of 360-240 Gy/150-100 Gy) delivered to the sensorimotor cortex of six 2-month-old naïve rats generated histologically evident cortical transections, without modifying motor behavior and weight gain up to 7 months. Microbeam transections of the sensorimotor cortex dramatically reduced convulsive seizure duration in a further group of 12 rats receiving local infusion of kainic acid. No subsequent neurological deficit was associated with the treatment. These data provide a novel tool to study the functions of the cortex and pave the way for the development of new therapeutic strategies for epilepsy and other neurological diseases.

Evaluation of dose-volume metrics for microbeam radiation therapy dose distributions in head phantoms of various sizes using Monte Carlo simulations

This work evaluates four dose-volume metrics applied to microbeam radiation therapy (MRT) using simulated dosimetric data as input. We seek to improve upon the most frequently used MRT metric, the peak-to-valley dose ratio (PVDR), by analyzing MRT dose distributions from a more volumetric perspective. Monte Carlo simulations were used to calculate dose distributions in three cubic head phantoms: a 2 cm mouse head, an 8 cm cat head and a 16 cm dog head. The dose distribution was calculated for a 4 × 4 mm² microbeam array in each phantom, as well as a 16 × 16 mm² array in the 8 cm cat head, and a 32 × 32 mm² array in the 16 cm dog head. Microbeam widths of 25, 50 and 75 µm and center-to-center spacings of 100, 200 and 400 µm were considered. The metrics calculated for each simulation were the conventional PVDR, the peak-to-mean valley dose ratio (PMVDR), the mean dose and the percentage volume below a threshold dose. The PVDR ranged between 3 and 230 for the 2 cm mouse phantom, and between 2 and 186 for the 16 cm dog phantom depending on geometry. The corresponding ranges for the PMVDR were much smaller, being 2-49 (mouse) and 2-46 (dog), and showed a slightly weaker dependence on phantom size and array size. The ratio of the PMVDR to the PVDR varied from 0.21 to 0.79 for the different collimation configurations, indicating a difference between the geometric dependence on outcome that would be predicted by these two metrics. For unidirectional irradiation, the mean lesion dose was 102%, 79% and 42% of the mean skin dose for the 2 cm mouse, 8 cm cat and 16 cm dog head phantoms, respectively. However, the mean lesion dose recovered to 83% of the mean skin dose in the 16 cm dog phantom in intersecting cross-firing regions. The percentage volume below a 10% dose threshold was highly dependent on geometry, with ranges for the different collimation configurations of 2-87% and 33-96% for the 2 cm mouse and 16 cm dog heads, respectively. The results of this study illustrate that different dose-volume metrics exhibit different functional dependences on MRT geometry parameters, and suggest that reliance on the PVDR as a predictor of therapeutic outcome may be insufficient.

Evaluating the peak-to-valley dose ratio of synchrotron microbeams using PRESAGE fluorescence

Synchrotron-generated microbeam radiotherapy holds great promise for future treatment, but the high dose gradients present conventional dosimetry with a challenge. Measuring the important peak-to-valley dose ratio (PVDR) of a microbeam-collimated synchrotron source requires both a dosimeter and an analysis method capable of exceptional spatial resolution. The PVDR is of great interest since it is the limiting factor for potential application of the microbeam radiation therapy technique clinically for its tissue-sparing properties (i.e. the valley dose should be below the tolerance of normal tissue). In this work a new method of measuring the dose response of PRESAGE dosimeters is introduced using the fluorescence from a 638 nm laser on a confocal laser-scanning microscope. This fluorescent microscopy method produces dosimetry data at a pixel size as low as 78 nm, giving a much better spatial resolution than optical computed tomography, which is normally used for scanning PRESAGE dosimeters. Using this technique the PVDR of the BL28B2 microbeam at the SPring-8 synchrotron in Japan is estimated to be approximately 52:1 at a depth of 2.5 mm. The PVDR was also estimated with EBT2 GAFchromic films as 30.5:1 at the surface in order to compare the PRESAGE fluorescent results with a more established dosimetry system. This estimation is in good agreement with previously measured ratios using other dosimeters and Monte Carlo simulations. This means that it is possible to use PRESAGE dosimeters with confocal microscopy for the determination of PVDR.