Proton linear energy transfer and variable relative biological effectiveness for adolescent patients with Hodgkin lymphoma.

Objectives: Proton therapy has a theoretical dosimetric advantage due to the Bragg peak, but the linear energy transfer (LET), and therefore the relative biological effectiveness (RBE), increase at the end of range. For patients with Hodgkin lymphoma, the distal edge of beam is often located within or close to the heart, where elevated RBE would be of potential concern. The purpose of this study was to investigate the impact of RBE and the choice of beam arrangement for adolescent patients with mediastinal Hodgkin lymphoma. Methods: For three previously treated adolescent patients, proton plans with 1-3 fields were created to a prescribed dose of 19.8 Gy (RBE) in 11 fractions (Varian Eclipse v13.7), assuming an RBE of 1.1. Plans were recalculated using Monte-Carlo (Geant4 v10.3.3/Gate v8.1) to calculate dose-averaged LET. Variable RBE-weighted dose was calculated using the McNamara model, assuming an α/ß ratio of 2 Gy for organs-at-risk. Results: Although the LET decreased as the number of fields increased, the difference in RBE-weighted dose (Δdose) to organs-at-risk did not consistently decrease. Δdose values varied by patient and organ and were mostly of the order of 0-3 Gy (RBE), with a worst-case of 4.75 Gy (RBE) in near-maximum dose to the left atrium for one plan. Conclusions: RBE-weighted doses to organs-at-risk are sensitive to the choice of RBE model, which is of particular concern for the heart. Advances in knowledge: There is a need to remain cautious when evaluating proton plans for Hodgkin lymphoma, especially when near-maximum doses to organs-at-risk are considered.

Proposing a Clinical Model for RBE Based on Proton Track-End Counts.

PURPOSE: In proton therapy, the clinical application of linear energy transfer (LET) optimization remains contentious, in part because of challenges associated with the definition and calculation of LET and its exact relationship with relative biological effectiveness (RBE) because of large variation in experimental in vitro data. This has raised interest in other metrics with favorable properties for biological optimization, such as the number of proton track ends in a voxel. In this work, we propose a novel model for clinical calculations of RBE, based on proton track end counts. METHODS AND MATERIALS: We developed an effective dose concept to translate between the total proton track-end count per unit mass in a voxel and a proton RBE value. Dose, track end, and dose-averaged LET (LETd) distributions were simulated using Monte Carlo models for a series of water phantoms, in vitro radiobiological studies, and patient treatment plans. We evaluated the correlation between track ends and regions of elevated biological effectiveness in comparison to LETd-based models of RBE. RESULTS: Track ends were found to correlate with biological effects in in vitro experiments with an accuracy comparable to LETd. In patient simulations, our track end model identified the same biological hotspots as predicted by LETd-based radiobiological models of proton RBE. CONCLUSIONS: These results suggest that, for clinical optimization and evaluation, an RBE model based on proton track end counts may match LETd-based models in terms of information provided while also offering superior statistical properties.

A preclinical radiotherapy dosimetry audit using a realistic 3D printed murine phantom.

Preclinical radiation research lacks standardized dosimetry procedures that provide traceability to a primary standard. Consequently, ensuring accuracy and reproducibility between studies is challenging. Using 3D printed murine phantoms we undertook a dosimetry audit of Xstrahl Small Animal Radiation Research Platforms (SARRPs) installed at 7 UK centres. The geometrically realistic phantom accommodated alanine pellets and Gafchromic EBT3 film for simultaneous measurement of the dose delivered and the dose distribution within a 2D plane, respectively. Two irradiation scenarios were developed: (1) a 10 × 10 mm2 static field targeting the pelvis, and (2) a 5 × 5 mm2 90° arc targeting the brain. For static fields, the absolute difference between the planned dose and alanine measurement across all centres was 4.1 ± 4.3% (mean ± standard deviation), with an overall range of - 2.3 to 10.5%. For arc fields, the difference was - 1.2% ± 6.1%, with a range of - 13.1 to 7.7%. EBT3 dose measurements were greater than alanine by 2.0 ± 2.5% and 3.5 ± 6.0% (mean ± standard deviation) for the static and arc fields, respectively. 2D dose distributions showed discrepancies to the planned dose at the field edges. The audit demonstrates that further work on preclinical radiotherapy quality assurance processes is merited.

Estimating the percentage of patients who might benefit from proton beam therapy instead of X-ray radiotherapy.

OBJECTIVES: High-energy Proton Beam Therapy (PBT) commenced in England in 2018 and NHS England commissions PBT for 1.5% of patients receiving radical radiotherapy. We sought expert opinion on the level of provision. METHODS: Invitations were sent to 41 colleagues working in PBT, most at one UK centre, to contribute by completing a spreadsheet. 39 responded: 23 (59%) completed the spreadsheet; 16 (41%) declined, arguing that clinical outcome data are lacking, but joined six additional site-specialist oncologists for two consensus meetings. The spreadsheet was pre-populated with incidence data from Cancer Research UK and radiotherapy use data from the National Cancer Registration and Analysis Service. 'Mechanisms of Benefit' of reduced growth impairment, reduced toxicity, dose escalation and reduced second cancer risk were examined. RESULTS: The most reliable figure for percentage of radical radiotherapy patients likely to benefit from PBT was that agreed by 95% of the 23 respondents at 4.3%, slightly larger than current provision. The median was 15% (range 4-92%) and consensus median 13%. The biggest estimated potential benefit was from reducing toxicity, median benefit to 15% (range 4-92%), followed by dose escalation median 3% (range 0 to 47%); consensus values were 12 and 3%. Reduced growth impairment and reduced second cancer risk were calculated to benefit 0.5% and 0.1%. CONCLUSIONS: The most secure estimate of percentage benefit was 4.3% but insufficient clinical outcome data exist for confident estimates. The study supports the NHS approach of using the evidence base and developing it through randomised trials, non-randomised studies and outcomes tracking. ADVANCES IN KNOWLEDGE: Less is known about the percentage of patients who may benefit from PBT than is generally acknowledged. Expert opinion varies widely. Insufficient clinical outcome data exist to provide robust estimates. Considerable further work is needed to address this, including international collaboration; much is already underway but will take time to provide mature data.

A Monte Carlo study of different LET definitions and calculation parameters for proton beam therapy.

The strongin vitroevidence that proton Relative Biological Effectiveness (RBE) varies with Linear Energy Transfer (LET) has led to an interest in applying LET within treatment planning. However, there is a lack of consensus on LET definition, Monte Carlo (MC) parameters or clinical methodology. This work aims to investigate how common variations of LET definition may affect potential clinical applications. MC simulations (GATE/GEANT4) were used to calculate absorbed dose and different types of LET for a simple Spread Out Bragg Peak (SOBP) and for four clinical PBT plans covering a range of tumour sites. Variations in the following LET calculation methods were considered: (i) averaging (dose-averaged LET (LETd) & track-averaged LET); (ii) scoring (LETdto water, to medium and to mass density); (iii) particle inclusion (LETdto all protons, to primary protons and to particles); (iv) MC settings (hit type and Maximum Step Size (MSS)). LET distributions were compared using: qualitative comparison, LET Volume Histograms (LVHs), single value criteria (maximum and mean values) and optimised LET-weighted dose models. Substantial differences were found between LET values in averaging, scoring and particle type. These differences depended on the methodology, but for one patient a difference of ∼100% was observed between the maximum LETdfor all particles and maximum LETdfor all protons within the brainstem in the high isodose region (4 keVµm-1and 8 keVµm-1respectively). An RBE model using LETdincluding heavier ions was found to predict substantially different LET-weighted dose compared to those using other LET definitions. In conclusion, the selection of LET definition may affect the results of clinical metrics considered in treatment planning and the results of an RBE model. The authors' advocate for the scoring of dose-averaged LET to water for primary and secondary protons using a random hit type and automated MSS.

Automated Monte-Carlo re-calculation of proton therapy plans using Geant4/Gate: implementation and comparison to plan-specific quality assurance measurements.

OBJECTIVES: Software re-calculation of proton pencil beam scanning plans provides a method of verifying treatment planning system (TPS) dose calculations prior to patient treatment. This study describes the implementation of AutoMC, a Geant4 v10.3.3/Gate v8.1 (Gate-RTion v1.0)-based Monte-Carlo (MC) system for automated plan re-calculation, and presents verification results for 153 patients (730 fields) planned within year one of the proton service at The Christie NHS Foundation Trust. METHODS: A MC beam model for a Varian ProBeam delivery system with four range-shifter options (none, 2 cm, 3 cm, 5 cm) was derived from beam commissioning data and implemented in AutoMC. MC and TPS (Varian Eclipse v13.7) calculations of 730 fields in solid-water were compared to physical plan-specific quality assurance (PSQA) measurements acquired using a PTW Octavius 1500XDR array and PTW 31021 Semiflex 3D ion chamber. RESULTS: TPS and MC showed good agreement with array measurements, evaluated using γ analyses at 3%, 3 mm with a 10% lower dose threshold:>94% of fields calculated by the TPS and >99% of fields calculated by MC had γ ≤ 1 for>95% of measurement points within the plane. TPS and MC also showed good agreement with chamber measurements of absolute dose, with systematic differences of <1.5% for all range-shifter options. CONCLUSIONS: Reliable independent verification of the TPS dose calculation is a valuable complement to physical PSQA and may facilitate reduction of the physical PSQA workload alongside a thorough delivery system quality assurance programme. ADVANCES IN KNOWLEDGE: A Gate/Geant4-based MC system is thoroughly validated against an extensive physical PSQA dataset for 730 clinical fields, showing that clinical implementation of MC for PSQA is feasible.

Preclinical dosimetry: exploring the use of small animal phantoms.

Preclinical radiotherapy studies using small animals are an indispensable step in the pathway from in vitro experiments to clinical implementation. As radiotherapy techniques advance in the clinic, it is important that preclinical models evolve to keep in line with these developments. The use of orthotopic tumour sites, the development of tissue-equivalent mice phantoms and the recent introduction of image-guided small animal radiation research platforms has enabled similar precision treatments to be delivered in the laboratory.These technological developments, however, are hindered by a lack of corresponding dosimetry standards and poor reporting of methodologies. Without robust and well documented preclinical radiotherapy quality assurance processes, it is not possible to ensure the accuracy and repeatability of dose measurements between laboratories. As a consequence current RT-based preclinical models are at risk of becoming irrelevant.In this review we explore current standardization initiatives, focusing in particular on recent developments in small animal irradiation equipment, 3D printing technology to create customisable tissue-equivalent dosimetry phantoms and combining these phantoms with commonly used detectors.

A Monte Carlo study on the collimation of pencil beam scanning proton therapy beams.

PURPOSE: The lateral edge of a proton therapy beam is commonly used to achieve conformality to the treatment volume where critical structures reside close to the target. However, when treating shallow depths, the lateral edge of a pencil beam scanning (PBS) system may be broader than that of a double scattered (DS) system. Use of a range-shifter to degrade the beam and allow treatment of very shallow depths further blurs the lateral edge. The authors investigate the potential use of a collimator with a PBS system for delivery of 3D uniform dose-volumes to a water-tank phantom, identifying the key factors controlling the width of the lateral edge. METHODS: The geant4 application for tomographic emission (gate) Monte Carlo (MC) environment was used, following validation against previously published data. Key parameters for PBS beams were investigated to assess their impact on the lateral edge of both monoenergetic beams and uniform dose-volumes. These parameters included nozzle-to-surface distance (NSD), vacuum window-to-surface distance (VSD), use of a range-shifter, and spot optimization parameters. RESULTS: The lateral edge of an uncollimated PBS beam is particularly sensitive to VSD and NSD. While use of a range-shifter blurs the lateral edge, collimation allows the edge to be sharpened to between 2 and 4 mm depending on the depth of the target. Optimization of the spot weightings alone can provide a penumbral width close to that of a single spot, but also leads to poorer uniformity near the edge of the target volume. CONCLUSIONS: Collimation of PBS beams should be considered for superficial targets particularly for beams delivered through a range-shifter, since the resultant sharpening of the lateral edge will allow improved sparing of adjacent normal tissues. Further work is needed to develop collimators which are integrated into both nozzle designs and planning system optimization algorithms.

Parametrized rectal dose and associations with late toxicity in prostate cancer radiotherapy.

OBJECTIVE: We investigated possible associations between planned dose-volume parameters and rectal late toxicity in 170 patients having radical prostate cancer radiotherapy. METHODS: For each patient, the rectum was outlined from anorectal junction to sigmoid colon, and rectal dose was parametrized using dose-volume (DVH), dose-surface (DSH) and dose-line (DLH) histograms. Generation of DLHs differed from previous studies in that the rectal dose was parametrized without first unwrapping onto 2-dimensional dose-surface maps. Patient-reported outcomes were collected using a validated Later Effects in Normal Tissues Subjective, Objective, Management and Analytic questionnaire. Associations between dose and toxicity were assessed using a one-sided Mann-Whitney U test. RESULTS: Associations (p < 0.05) were found between equieffective dose (EQD23) and late toxicity as follows: overall toxicity with DVH and DSH at 13-24 Gy; proctitis with DVH and DSH at 25-36 Gy and with DVH, DSH and DLH at 61-67 Gy; bowel urgency with DVH and DSH at 10-20 Gy. None of these associations met statistical significance following the application of a Bonferroni correction. CONCLUSION: Independently confirmed associations between rectal dose and late toxicity remain elusive. Future work to increase the accuracy of the knowledge of the rectal dose, either by accounting for interfraction and intrafraction rectal motion or via stabilization of the rectum during treatment, may be necessary to allow for improved dose-toxicity comparisons. ADVANCES IN KNOWLEDGE: This study is the first to use parametrized DLHs to study associations with patient-reported toxicity for prostate radiotherapy showing that it is feasible to model rectal dose mapping in three dimensions.

Simulation of realistic linac motion improves the accuracy of a Monte Carlo based VMAT plan QA system.

PURPOSE: To investigate the use of a software-based pre-treatment QA system for VMAT, which incorporates realistic linac motion during delivery. METHODS: A beam model was produced using the GATE platform for GEANT4 Monte Carlo dose calculations. Initially validated against static measurements, the model was then integrated with a VMAT delivery emulator, which reads plan files and generates a set of dynamic delivery instructions analogous to the linac control system. Monte Carlo simulations were compared to measurements on dosimetric phantoms for prostate and head and neck VMAT plans. Comparisons were made between calculations using fixed control points, and simulations of continuous motion utilising the emulator. For routine use, the model was incorporated into an automated pre-treatment QA system. RESULTS: The model showed better agreement with measurements when incorporating linac motion: mean gamma pass (Γ<1) over 5 prostate plans was 100.0% at 3%/3mm and 97.4% at 2%/2mm when compared to measurement. For the head and neck plans, delivered to the anatomical phantom, gamma passes were 99.4% at 4%/4mm and 94.94% at 3%/3mm. For example simulations within patient CT data, gamma passes were observed which are within our centre's tolerance for pre-treatment QA. CONCLUSIONS: Through comparison to phantom measurements, it was found that the incorporation of a realistic linac motion improves the accuracy of the model compared to the simulation of fixed control points. The ability to accurately calculate dose as a second check of the planning system, and determine realistic delivery characteristics, may allow for the reduction of machine-based pre-treatment plan QA for VMAT.

An analysis of the origin of differences between measured and simulated fields produced by a 15-element ultrasound phased array.

Modeling provides an attractive approach for the design of phased array ultrasound transducers for hyperthermia. However, measurements on physical transducers reveal differences from the idealized field profiles predicted by simulation. In this paper we report a method of analyzing the origins of these differences. The measured performance of a 15-element sparse phased array is described and compared with simulated fields calculated using the point source method. It highlighted two notable differences: First, that the focal region was located closer to the surface of the physical transducer than in the simulated fields; and second, that numerous intensity maxima were present between the surface of the transducer and the focal zone in the experimental data, but not in the simulated fields. We identified six factors that could potentially affect the field but were not taken into account by the default simulations, and we performed a sensitivity analysis on these: (i) Variation in the amplitude of the output from each element, (ii) the presence of square-wave harmonics in the drive signals, (iii) nonpistonlike vibration of elements, (iv) quantization of the applied phases, (v) errors in the spatial positioning of each element; and (vi) interelement cross-coupling. Both the independent impact of each factor and the interactions between multiple factors were analyzed by using a full-factorial experimental design composed of 64 (2(6)) simulations. The results indicated that nonpistonlike motion of elements is likely to be the primary cause of differences between the measured and modelled fields. Determination of the precise vibrational modes of elements in an array is complex and would require full finite element analysis. However, the simple vibrational mode considered within the present work, corresponding to the addition of a surface Rayleigh wave originating at the element center and propagating radially, produced simulation results that were in good agreement with the measured data.

The design and characterization of an ultrasound phased array suitable for deep tissue hyperthermia.

In this paper we describe the design and evaluation of a planar phased-array ultrasound transducer suitable for producing localized hyperthermia in solid tumors deep within the body. Simulation using a customized version of Ultrasim has been used to determine the relationship between the size and position of the focus and parameters of the array. These parameters include the overall size of the array and the size, shape and distribution of the individual elements. A 15-element prototype array has been constructed using the results of the simulation. Beam profile measurements on this transducer made in an acoustic tank were compared with the beam profile predicted by simulation. The results showed good agreement in the shape of the focal region, but with the focus closer to the surface of the physical transducer when compared with the simulation and with small high-intensity areas between the surface of the transducer and the focus in the measured profile. A sensitivity analysis using a simulated factorial experiment indicated that the presence of a secondary vibrational mode within the elements of the array was the principal cause for both the shift in the position of the focus and for the unwanted maxima close to the surface of the array. The results also showed that the array was tolerant of a large variation in output intensity of the individual elements in the array in producing a focal region. Extrapolation of the results obtained in this study indicate that an array of 60 elements, based on the design described, driven by 550 V peak-to-peak pulses would be capable of producing a peak focal intensity of 50 Wcm(-2) at a depth of 60 mm in tissue, which would be appropriate for hyperthermia used as an adjunct to radiotherapy or chemotherapy.