CERN-based experiments and Monte-Carlo studies on focused dose delivery with very high energy electron (VHEE) beams for radiotherapy applications.

Very High Energy Electron (VHEE) beams are a promising alternative to conventional radiotherapy due to their highly penetrating nature and their applicability as a modality for FLASH (ultra-high dose-rate) radiotherapy. The dose distributions due to VHEE need to be optimised; one option is through the use of quadrupole magnets to focus the beam, reducing the dose to healthy tissue and allowing for targeted dose delivery at conventional or FLASH dose-rates. This paper presents an in depth exploration of the focusing achievable at the current CLEAR (CERN Linear Electron Accelerator for Research) facility, for beam energies >200 MeV. A shorter, more optimal quadrupole setup was also investigated using the TOPAS code in Monte Carlo simulations, with dimensions and beam parameters more appropriate to a clinical situation. This work provides insight into how a focused VHEE radiotherapy beam delivery system might be achieved.

Characterisation of the UK high energy proton research beamline for high and ultra-high dose rate (FLASH) irradiation.

Objective. This work sets out the capabilities of the high energy proton research beamline developed in the Christie proton therapy centre for Ultra-High Dose Rate (UHDR) irradiation and FLASH experiments. It also characterises the lower limits of UHDR operation for this Pencil Beam Scanning (PBS) proton hardware.Approach. Energy dependent nozzle transmission was measured using a Faraday Cup beam collector. Spot size was measured at the reference plane using a 2D scintillation detector. Integrated depth doses (IDDs) were measured. EBT3 Gafchromic film was used to compare UHDR and conventional dose rate spots. Our beam monitor calibration methodolgy for UHDR is described. A microDiamond detector was used to determine dose rates at zref. Instantaneous depth dose rates were calculated for 70-245 MeV. PBS dose rate distributions were calculated using Folkerts and Van der Water definitions.Main results. Transmission of 7.05 ± 0.1% is achieveable corresponding to a peak instantaneous dose rate of 112.7 Gy s-1. Beam parameters are comparable in conventional and UHDR mode with a spot size ofσx= 4.6 mm,σy= 6.6 mm. Dead time in the beam monitoring electonics warrants a beam current dependent MU correction in the present configuration. Fast beam scanning of 26.4 m s-1(X) and 12.1 m s-1(Y) allows PBS dose rates of the order tens of Grays per second.Significance. UHDR delivery is possible for small field sizes and high energies enabling research into the FLASH effect with PBS protons at our facility. To our knowledge this is also the first thorough characterisation of UHDR irradiation using the hardware of this clinical accelerator at energies less than 250 MeV. The data set out in this publication can be used for designing experiments at this UK research facility and inform the possible future clinical translation of UHDR PBS proton therapy.

Focused VHEE (very high energy electron) beams and dose delivery for radiotherapy applications.

This paper presents the first demonstration of deeply penetrating dose delivery using focused very high energy electron (VHEE) beams using quadrupole magnets in Monte Carlo simulations. We show that the focal point is readily modified by linearly changing the quadrupole magnet strength only. We also present a weighted sum of focused electron beams to form a spread-out electron peak (SOEP) over a target region. This has a significantly reduced entrance dose compared to a proton-based spread-out Bragg peak (SOBP). Very high energy electron (VHEE) beams are an exciting prospect in external beam radiotherapy. VHEEs are less sensitive to inhomogeneities than proton and photon beams, have a deep dose reach and could potentially be used to deliver FLASH radiotherapy. The dose distributions of unfocused VHEE produce high entrance and exit doses compared to other radiotherapy modalities unless focusing is employed, and in this case the entrance dose is considerably improved over existing radiations. We have investigated both symmetric and asymmetric focusing as well as focusing with a range of beam energies.

Determining the parameter space for effective oxygen depletion for FLASH radiation therapy.

There has been a recent revival of interest in the FLASH effect, after experiments have shown normal tissue sparing capabilities of ultra-high-dose-rate radiation with no compromise on tumour growth restraint. A model has been developed to investigate the relative importance of a number of fundamental parameters considered to be involved in the oxygen depletion paradigm of induced radioresistance. An example eight-dimensional parameter space demonstrates the conditions under which radiation may induce sufficient depletion of oxygen for a diffusion-limited hypoxic cellular response. Initial results support experimental evidence that FLASH sparing is only achieved for dose rates on the order of tens of Gy s-1or higher, for a sufficiently high dose, and only for tissue that is slightly hypoxic at the time of radiation. We show that the FLASH effect is the result of a number of biological, radiochemical and delivery parameters. Also, the threshold dose for a FLASH effect occurring would be more prominent when the parameterisation was optimised to produce the maximum effect. The model provides a framework for further FLASH-related investigation and experimental design. An understanding of the mechanistic interactions producing an optimised FLASH effect is essential for its translation into clinical practice.

In Silico Models of DNA Damage and Repair in Proton Treatment Planning: A Proof of Concept.

There is strong in vitro cell survival evidence that the relative biological effectiveness (RBE) of protons is variable, with dependence on factors such as linear energy transfer (LET) and dose. This is coupled with the growing in vivo evidence, from post-treatment image change analysis, of a variable RBE. Despite this, a constant RBE of 1.1 is still applied as a standard in proton therapy. However, there is a building clinical interest in incorporating a variable RBE. Recently, correlations summarising Monte Carlo-based mechanistic models of DNA damage and repair with absorbed dose and LET have been published as the Manchester mechanistic (MM) model. These correlations offer an alternative path to variable RBE compared to the more standard phenomenological models. In this proof of concept work, these correlations have been extended to acquire RBE-weighted dose distributions and calculated, along with other RBE models, on a treatment plan. The phenomenological and mechanistic models for RBE have been shown to produce comparable results with some differences in magnitude and relative distribution. The mechanistic model found a large RBE for misrepair, which phenomenological models are unable to do. The potential of the MM model to predict multiple endpoints presents a clear advantage over phenomenological models.

Mechanistic modelling supports entwined rather than exclusively competitive DNA double-strand break repair pathway.

Following radiation induced DNA damage, several repair pathways are activated to help preserve genome integrity. Double Strand Breaks (DSBs), which are highly toxic, have specified repair pathways to address them. The main repair pathways used to resolve DSBs are Non-Homologous End Joining (NHEJ) and Homologous Recombination (HR). Cell cycle phase determines the availability of HR, but the repair choice between pathways in the G2 phases where both HR and NHEJ can operate is not clearly understood. This study compares several in silico models of repair choice to experimental data published in the literature, each model representing a different possible scenario describing how repair choice takes place. Competitive only scenarios, where initial protein recruitment determines repair choice, are unable to fit the literature data. In contrast, the scenario which uses a more entwined relationship between NHEJ and HR, incorporating protein co-localisation and RNF138-dependent removal of the Ku/DNA-PK complex, is better able to predict levels of repair similar to the experimental data. Furthermore, this study concludes that co-localisation of the Mre11-Rad50-Nbs1 (MRN) complexes, with initial NHEJ proteins must be modeled to accurately depict repair choice.

Clinically relevant nanodosimetric simulation of DNA damage complexity from photons and protons.

Relative Biological Effectiveness (RBE), the ratio of doses between radiation modalities to produce the same biological endpoint, is a controversial and important topic in proton therapy. A number of phenomenological models incorporate variable RBE as a function of Linear Energy Transfer (LET), though a lack of mechanistic description limits their applicability. In this work we take a different approach, using a track structure model employing fundamental physics and chemistry to make predictions of proton and photon induced DNA damage, the first step in the mechanism of radiation-induced cell death. We apply this model to a proton therapy clinical case showing, for the first time, predictions of DNA damage on a patient treatment plan. Our model predictions are for an idealised cell and are applied to an ependymoma case, at this stage without any cell specific parameters. By comparing to similar predictions for photons, we present a voxel-wise RBE of DNA damage complexity. This RBE of damage complexity shows similar trends to the expected RBE for cell kill, implying that damage complexity is an important factor in DNA repair and therefore biological effect.

Image Guidance for Proton Therapy.

Image-guided radiotherapy has an established role in all forms of radiotherapy treatment delivery. Proton therapy seeks to deliver superior dose distributions through utilising the Bragg peak to target tumour and avoid sensitive normal tissue. The Bragg peak and sharp falloff in dose delivered by proton therapy necessitate careful treatment planning and treatment delivery. The dose distribution delivered by proton therapy is particularly sensitive to uncertainty in the prediction of proton range during treatment planning and deviations from the planned delivery during the course of the fractionated treatment. Realising the superior dose distribution of proton therapy requires increased diligence and image guidance has a key role in ensuring that treatments are planned and delivered. This article will outline the current status of image guidance for proton therapy, particularly highlighting differences with regard to high-energy X-ray therapy, and will look at a number of future improvements in image-guided proton therapy.

In Silico Non-Homologous End Joining Following Ion Induced DNA Double Strand Breaks Predicts That Repair Fidelity Depends on Break Density.

This work uses Monte Carlo simulations to investigate the dependence of residual and misrepaired double strand breaks (DSBs) at 24 hours on the initial damage pattern created during ion therapy. We present results from a nanometric DNA damage simulation coupled to a mechanistic model of Non-Homologous End Joining, capable of predicting the position, complexity, and repair of DSBs. The initial damage pattern is scored by calculating the average number of DSBs within 70 nm from every DSB. We show that this local DSB density, referred to as the cluster density, can linearly predict misrepair regardless of ion species. The models predict that the fraction of residual DSBs is constant, with 7.3% of DSBs left unrepaired following 24 hours of repair. Through simulation over a range of doses and linear energy transfer (LET) we derive simple correlations capable of predicting residual and misrepaired DSBs. These equations are applicable to ion therapy treatment planning where both dose and LET are scored. This is demonstrated by applying the correlations to an example of a clinical proton spread out Bragg peak. Here we see a considerable biological effect past the distal edge, dominated by residual DSBs.

Nanodosimetric Simulation of Direct Ion-Induced DNA Damage Using Different Chromatin Geometry Models.

Monte Carlo based simulation has proven useful in investigating the effect of proton-induced DNA damage and the processes through which this damage occurs. Clustering of ionizations within a small volume can be related to DNA damage through the principles of nanodosimetry. For simulation, it is standard to construct a small volume of water and determine spatial clusters. More recently, realistic DNA geometries have been used, tracking energy depositions within DNA backbone volumes. Traditionally a chromatin fiber is built within the simulation and identically replicated throughout a cell nucleus, representing the cell in interphase. However, the in vivo geometry of the chromatin fiber is still unknown within the literature, with many proposed models. In this work, the Geant4-DNA toolkit was used to build three chromatin models: the solenoid, zig-zag and cross-linked geometries. All fibers were built to the same chromatin density of 4.2 nucleosomes/11 nm. The fibers were then irradiated with protons (LET 5-80 keV/µm) or alpha particles (LET 63-226 keV/µm). Nanodosimetric parameters were scored for each fiber after each LET and used as a comparator among the models. Statistically significant differences were observed in the double-strand break backbone size distributions among the models, although nonsignificant differences were noted among the nanodosimetric parameters. From the data presented in this article, we conclude that selection of the solenoid, zig-zag or cross-linked chromatin model does not significantly affect the calculated nanodosimetric parameters. This allows for a simulation-based cell model to make use of any of these chromatin models for the scoring of direct ion-induced DNA damage.

Marvin: an anatomical phantom for dosimetric evaluation of complex radiotherapy of the head and neck.

We report on the design of Marvin, a Model Anatomy for Radiotherapy Verification and audit In the head and Neck and present results demonstrating its use in the development of the Elekta volumetric modulated arc therapy (VMAT) technique at the Christie, and in the audit of TomoTherapy and Varian RapidArc at other institutions. The geometry of Marvin was generated from CT datasets of eight male and female patients lying in the treatment position, with removable inhomogeneities modelling the sinuses and mandible. A modular system allows the phantom to be used with a range of detectors, with the locations of the modules being based on an analysis of a range of typical treatment plans (27 in total) which were mapped onto the phantom geometry. Results demonstrate the use of Gafchromic EBT2/EBT3 film for measurement of relative dose in a plane through the target and organs-at-risk, and the use of a small-volume ionization chamber for measurement of absolute dose in the target and spinal cord. Measurements made during the development of the head and neck VMAT protocol at the Christie quantified the improvement in plan delivery resulting from the installation of the Elekta Integrity upgrade (which permits an effectively continuously variable dose rate), with plans delivered before and after the upgrade having 88.5 ± 9.4% and 98.0 ± 2.2% respectively of points passing a gamma analysis (at 4%, 4 mm, global). Audits of TomoTherapy and Varian RapidArc neck techniques at other institutions showed a similar quality of plan delivery as for post-Integrity Elekta VMAT: film measurements for both techniques had >99% of points passing a gamma analysis at the clinical criteria of 4%, 4 mm, global, and >95% of points passing at tighter criteria of 3%, 3 mm, global; and absolute dose measurements in the PTV and spinal cord were within 1.5% and 3.5% of the planned doses respectively for both techniques. The results demonstrate that Marvin is an efficient and effective means of assessing the quality of delivery of complex radiotherapy in the head and neck, and is a useful tool to assist development and audit of these techniques.

Modelling the throughput capacity of a single-accelerator multitreatment room proton therapy centre.

OBJECTIVE: We describe a model for evaluating the throughput capacity of a single-accelerator multitreatment room proton therapy centre with the aims of (1) providing quantitative estimates of the throughput and waiting times and (2) providing insight into the sensitivity of the system to various physical parameters. METHODS: A Monte Carlo approach was used to compute various statistics about the modelled centre, including the throughput capacity, fraction times for different groups of patients and beam waiting times. A method of quantifying the saturation level is also demonstrated. RESULTS: Benchmarking against the MD Anderson Cancer Center showed good agreement between the modelled (140 ± 4 fractions per day) and reported (133 ± 35 fractions per day) throughputs. A sensitivity analysis of that system studied the impact of beam switch time, the number of treatment rooms, patient set-up times and the potential benefit of having a second accelerator. Finally, scenarios relevant to a potential UK facility were studied, finding that a centre with the same four-room, single-accelerator configuration as the MD Anderson Cancer Center but handling a more complex UK-type caseload would have a throughput reduced by approximately 19%, but still be capable of treating in excess of 100 fractions per 16-h treatment day. CONCLUSIONS: The model provides a useful tool to aid in understanding the operating dynamics of a proton therapy facility, and for investigating potential scenarios for prospective centres. ADVANCES IN KNOWLEDGE: The model helps to identify which technical specifications should be targeted for future improvements.

Radiotherapy physics research in the UK: challenges and proposed solutions.

In 2011, the Clinical and Translational Radiotherapy Research Working Group (CTRad) of the National Cancer Research Institute brought together UK radiotherapy physics leaders for a think tank meeting. Following a format that CTRad had previously and successfully used with clinical oncologists, 23 departments were asked to complete a pre-meeting evaluation of their radiotherapy physics research infrastructure and the strengths, weaknesses, opportunities and threats within their own centre. These departments were brought together with the CTRad Executive Group and research funders to discuss the current state of radiotherapy physics research, perceived barriers and possible solutions. In this Commentary, we summarise the submitted materials, presentations and discussions from the meeting and propose an action plan. It is clear that there are challenges in both funding and staffing of radiotherapy physics research. Programme and project funding streams sometimes struggle to cater for physics-led work, and increased representation on research funding bodies would be valuable. Career paths for academic radiotherapy physicists need to be examined and an academic training route identified within Modernising Scientific Careers; the introduction of formal job plans may allow greater protection of research time, and should be considered. Improved access to research facilities, including research linear accelerators, would enhance research activity and pass on developments to patients more quickly; research infrastructure could be benchmarked against centres in the UK and abroad. UK National Health Service departments wishing to undertake radiotherapy research, with its attendant added value for patients, need to develop a strategy with their partner higher education institution, and collaboration between departments may provide enhanced opportunities for funded research.

The use of a realistic VMAT delivery emulator to optimize dynamic machine parameters for improved treatment efficiency.

The delivery of volumetric modulated arc therapy (VMAT) requires the simultaneous movement of the linear accelerator gantry, multi-leaf collimators and jaws while the dose rate is varied. In this study, a VMAT delivery emulator was developed to accurately predict the characteristics of a given treatment plan, incorporating realistic parameters for gantry inertia and the variation in leaf speed with respect to gravity. The emulator was used to assess the impact of dynamic machine parameters on the delivery efficiency, using a set of prostate and head and neck VMAT plans. Initially, assuming a VMAT system with fixed dose rate bins, the allowable leaf and jaw speeds were increased and a significant improvement in treatment time and average dose rate was observed. The software was then adapted to simulate a VMAT system with continuously varying dose rate, and the increase in delivery efficiency was quantified, along with the impact of an increased leaf and jaw speed. Finally, a set of optimal dynamic machine parameters was derived assuming an idealized scenario in which the treatment is delivered in a single arc at constant maximum gantry speed.

UK guidelines for the safe delivery of intensity-modulated radiotherapy.

This paper outlines the guidelines for the development of intensity-modulated radiotherapy (IMRT) in the UK. The guidelines are designed to cover the complete implementation of IMRT, with guidelines in the following categories: commissioning, quality, clinical, audit, and training and education. These guidelines have been compiled by the Radiotherapy Development Board of the Royal College of Radiologists and will support the safe application of IMRT in the UK.

A novel imaging technique for fusion of high-quality immobilised MR images of the head and neck with CT scans for radiotherapy target delineation.

Uncertainty and inconsistency are observed in target volume delineation in the head and neck for radiotherapy treatment planning based only on CT imaging. Alternative modalities such as MRI have previously been incorporated into the delineation process to provide additional anatomical information. This work aims to improve on previous studies by combining good image quality with precise patient immobilisation in order to maintain patient position between scans. MR images were acquired using quadrature coils placed over the head and neck while the patient was immobilised in the treatment position using a five-point thermoplastic shell. The MR image and CT images were automatically fused in the Pinnacle treatment planning system using Syntegra software. Image quality, distortion and accuracy of the image registration using patient anatomy were evaluated. Image quality was found to be superior to that acquired using the body coil, while distortion was < 1.0 mm to a radius of 8.7 cm from the scan centre. Image registration accuracy was found to be 2.2 mm (+/- 0.9 mm) and < 3.0 degrees (n = 6). A novel MRI technique that combines good image quality with patient immobilization has been developed and is now in clinical use. The scan duration of approximately 15 min has been well tolerated by all patients.

The cost effectiveness of in vivo dosimetry is not proven.

In vivo dosimetry is a check of the radiotherapy treatment process. Recently, it has been argued that in vivo dosimetry at the start of every patient's radiotherapy is cost effective. We have examined and extended the cost effectiveness argument and show that the model does not deliver the stated benefits.