Deisher, A.

BACKGROUND: Cardiac radioablation is emerging as a treatment modality for refractory ventricular tachycardia. This study aimed to evaluate the effects of radiation on the coronary arteries in a swine model of proton beam cardiac radioablation. METHODS: Eighteen swine underwent single-fraction 30 to 40 Gy pencil-beam scanning proton therapy targeting the left ventricle and were euthanized 12 to 40 weeks later. Treatment planning was performed without restricting the dose to the coronary arteries. The maximum point dose (D max ) to the epicardial coronary arteries was calculated. In secondary analyses, the mean (D mean ) and minimum dose received by the highest irradiated 0.01 cm 3 (D 0.01 ) of each coronary artery were also calculated. Coronary artery segments were harvested from the D max sites for histological analysis, and the D max was correlated with stenosis severity. RESULTS: Ninety-six coronary arteries were analyzed. No stenoses were observed by computed tomography imaging preirradiation. By histological analysis posteuthanasia, 25/96 (26%) coronaries sampled at their D max sites had ≥75% stenosis. The median D max was 4.7 Gy for the <75% stenosis group and 29.7 Gy for the ≥75% stenosis group ( P <0.001). The AUC-ROC for the association between D max and stenosis ≥75% was 92.2%. A D max value of 20.1 Gy best predicted stenosis ≥75%, with sensitivity 92.3% and specificity 87.1%. The AUC-ROCs for the associations of D mean and D 0.01 with stenosis ≥75% were 84.8% and 91.6%, respectively. In histopathologic analysis, intimal hyperplasia was the most common coronary artery abnormality at the D max sites, and it was present in 61.5% of all arteries and in 93.9% of arteries with D max ≥20 Gy. CONCLUSIONS: In this preclinical model of proton beam cardiac radioablation, coronary stenoses occurred in a dose-dependent manner, with D max showing the closest correlation with stenosis ≥75%. These data provide for the first time a framework for dose constraint considerations for the coronary arteries during treatment planning for cardiac radioablation and thoracic malignancy radiation.

Background: Cardiac radioablation (CRA) is emerging as a treatment modality for refractory ventricular tachycardia (VT). This study aimed to evaluate the effects of radiation on the coronary arteries in a swine model of proton beam CRA. Methods: Eighteen domestic swine underwent 30–40 Gy pencil-beam scanning proton therapy in a single-fraction targeting the left ventricle with no attempt to spare coronary arteries and were euthanized 12-40 weeks later. Contrast-enhanced cardiac computed tomography (CT) was performed at baseline for treatment planning and at study end. Using dose-volume histograms, the maximum point (D max ), mean (D mean ) doses and the minimal doses received by the highest irradiated volumes of 0.01 cc (D 0.01cc ) for the left main, left anterior descending, circumflex, and right coronary artery, as well as the large diagonal and obtuse marginal branches were calculated. The D max location for each artery was mapped from the baseline planning CT onto the follow-up CT using deformable registration. Coronary artery segments were harvested from the D max sites for histological analysis. Receiver operating characteristic (ROC) analysis was performed to identify dose thresholds for predicting histologically severe coronary stenoses (≥75%). Results: Ninety-six coronary artery segments were analyzed. No stenoses were observed by CT imaging before irradiation. By histological analysis post-euthanasia, 25/96 (26%) coronary segments sampled at the Dmax site per artery had ≥75% stenosis. The median D max was 4.7 Gy for the <75% stenosis group and 29.7 Gy for the ≥75% stenosis group (p<0.001). The AUC-ROC for the association between D max and severe stenosis was 92.2%. A D max value of 20.1 Gy in a single fraction best predicted severe stenosis, with sensitivity of 92.3% and specificity of 87.1%. The AUC-ROC for the association between severe stenosis versus D mean and D 0.01cc was 84.8% and 91.6%, respectively. Based on histopathologic analysis, intimal hyperplasia was the most common coronary artery abnormality at the D max sites. Conclusion: In this pre-clinical model of proton beam CRA, coronary stenoses occurred in a dose-dependent manner, with D max being a better dosimetric predictor of stenosis than D mean and D 0.01cc . These data provide for the first time radiation dose constraint information for the epicardial coronary arteries that can be considered in treatment planning for CRA and radiation for thoracic malignancies.

Background: Proton beam irradiation to the ventricular myocardium may enable non-invasive targeting of ventricular tachycardia. Previous work has demonstrated the feasibility of imaging ablation lesion with late gadolinium enhanced (LGE) cardiac magnetic resonance imaging (MRI). Hypothesis: We hypothesize that proton beam irradiation will lead to progressively increasing density of LGE corresponding to lesion formation during follow-up after treatment delivery. Aim: The aim of the current study was to determine the time course of LGE intensity in ventricular myocardium following proton beam therapy. Methods: Eleven swine were irradiated with 40 Gy proton beam in the anterior or lateral left ventricular (LV) followed by LGE cardiac MRI every 4 weeks for up to 24 weeks. After manual contouring of the endocardial and epicardial myocardium, a region of interest (ROI) in the area irradiated with 0-5 Gy was contoured (remote ROI). The mean and standard deviation (SD) of the signal intensity (SI) of the remote ROI were calculated. Light LGE was defined as voxels with SI > mean+2SD but < mean+3SD of the remote ROI, while voxels with SI > mean+3SD of the remote ROI were classified as having dense LGE. Voxels identified as bright due to image artifacts or noise were manually removed. The analysis was conducted in a fashion blinded to swine ID and time point after irradiation. Results: A total of 33 time points were analyzed (3.1 ± 1.4 MRI studies per swine). Figure 1A shows a representative time course of LGE in the same swine, progressing from diffuse light LGE at 8 weeks to dense LGE at 16-24 weeks. Figure 1B demonstrates the proportion of light versus dense LGE across timepoints in irradiated myocardial regions. At 8 weeks, light and dense LGE areas were evenly distributed. At 12 to 24 weeks, dense LGE increased to approximately 75%. The proportion of dense versus light LGE was significantly higher at 12, 16, and 24 weeks, as compared to 8 weeks. Conclusion: In ventricular myocardium, diffuse light LGE evident at 8 weeks progresses to denser LGE at 12-24 weeks following proton beam irradiation corresponding with lesion formation.

Background: Particle therapy is a noninvasive, catheter-free modality for cardiac ablation. We previously demonstrated the efficacy for creating ablation lesions in the porcine heart. Despite several earlier studies, the exact mechanism of early biophysical effects of proton and photon beam delivery on the myocardium remain incompletely resolved. Methods: Ten normal and 9 infarcted in situ porcine hearts received proton beam irradiation (40 Gy) delivered to the left ventricular myocardium with follow-up for 8 weeks. High-resolution electroanatomical mapping of the left ventricular was performed at baseline and follow-up. Bipolar voltage amplitude, conduction velocity, and connexin-43 were determined within the irradiated and nonirradiated areas. Results: The irradiated area in normal hearts showed a significant reduction of bipolar voltage amplitude (10.1±4.9 mV versus 5.7±3.2, P <0.0001) and conduction velocity (85±26 versus 55±13 cm/s, P =0.03) beginning at 4 weeks after irradiation. In infarcted myocardium after irradiation, bipolar voltage amplitude of the infarct scar (2.0±2.9 versus 0.8±0.7 mV, P =0.008) was significantly reduced as well as the conduction velocity in the infarcted heart (43.7±15.7 versus 26.3±11.4 cm/s, P =0.02). There were no significant changes in bipolar voltage amplitude and conduction velocity in nonirradiated myocardium. Myocytolysis, capillary hyperplasia, and dilation were seen in the irradiated myocardium 8 weeks after irradiation. Active caspase-3 and reduction of connexin-43 expression began in irradiated myocardium 1 week after irradiation and decreased over 8 weeks. Conclusions: Irradiation of the myocardium with proton beams reduce connexin-43 expression, conduction velocity, and bipolar conducted electrogram amplitude in a large porcine model. The changes in biomarkers preceded electrophysiological changes after proton beam therapy.

AbstractBackground and purposeProton therapy is a key modality used in the treatment of ocular melanoma. Traditionally ocular sites are treated using a dedicated eyeline with a passively scattered proton beam and a brass aperture. This work aims to design and characterize a beam‐collimating aperture to treat ocular targets with a gantry‐based spot scanning proton beam.MethodsA plastic aperture system that slides into the gantry nozzle of a spot scanning proton beam was designed and constructed. It consists of an intermediate scraper layer to attenuate stray protons and a 3D‐printed patient‐specific aperture positioned 5.7 cm from the surface of the eye. The aperture system was modeled in TOPAS and Monte Carlo simulations were validated with film measurements. Two different spot configurations were investigated for treatment planning and characterized based on lateral penumbra, central axis (CAX) dose and relative efficiency. Alignment and leakage were investigated through experimental film measurements. Range was verified using a multi‐layer ionization chamber. Reference dose measurements were made with a PinPoint 3D ion chamber. Neutron dose was evaluated through Monte Carlo simulations.ResultsAperture alignment with radiation isocenter was determined to be within 0.31 mm at a gantry angle of 0°. A single‐spot configuration with a 10 mm diameter aperture yielded film‐measured lateral penumbras of 1 mm to 1.25 mm, depending on depth in the spread‐out Bragg peak. TOPAS simulations found that a single spot configuration results in a flat dose distribution for a 10 mm diameter aperture and provides a CAX dose of less than 106% for apertures less than 14 mm in diameter. For larger targets, adding four corner spots to fill in the dose distribution is beneficial. Trade‐offs between lateral penumbra, CAX dose and relative efficiency were characterized for different spot configurations and can be used for future clinical decision‐making. The aperture was experimentally determined to not affect proton beam range, and no concerning leakage radiation or neutron dose was identified. Reference dose measurements with a PinPoint ion chamber were within 2.1% of Monte Carlo calculated doses.ConclusionThe aperture system developed in this work provides a method of treating ocular sites on a gantry‐based spot scanning proton system. Additional work to develop compatible gaze tracking and gating infrastructure is ongoing.

Purpose: This study evaluates beam angles used to generate highly individualized proton therapy treatment plans for patients eligible for carbon ion radiotherapy (CIRT). Methods and Materials: We retrospectively evaluated patients treated with pencil beam scanning intensity modulated proton therapy from 2015 to 2020 who had indications for CIRT. Patients were treated with a 190° rotating gantry with a robotic patient positioning system. Treatment plans were individualized to provide maximal prescription dose delivery to the tumor target volume while sparing organs at risk. The utilized beam angles were grouped, and anatomic sites with at least 10 different beam angles were sorted into histograms. Results: A total of 467 patients with 484 plans and 1196 unique beam angles were evaluated and characterized by anatomic treatment site and the number of beam angles utilized. The most common beam angles used were 0° and 180°. A wide range of beam angles were used in treating almost all anatomic sites. Only esophageal cancers had a predominantly unimodal grouping of beam angles. Pancreas cancers showed a modest grouping of beam angles. Conclusions: The wide distribution of beam angles used to treat CIRT-eligible patients suggests that a rotating gantry is optimal to provide highly individualized beam arrangements.

AbstractWith complex head and neck diseases that closely approximate crucial organs or have advanced spread and cannot be surgically resected, intensity modulated proton beam radiotherapy is favored due to its superior dose targeting through magnetic manipulation of the proton energy. A radiation mask and an oral positioning device immobilize craniofacial, cervical, and oral structures for accurate and reliable radiation delivery. Widely available prefabricated thermoplastic oral positioning devices are manufactured in standardized forms and materials that unpredictably affect the proton beams’ pathways and range. This technique article showcases a workflow that combines analog and digital dental techniques to produce a customized 3D‐printed oral positioning device within 2 appointments.

Background: Proton beam therapy offers radiophysical properties that are appealing for noninvasive arrhythmia elimination. This study was conducted to use scanned proton beams for ablation of cardiac tissue, investigate electrophysiological outcomes, and characterize the process of lesion formation in a porcine model using particle therapy. Methods: Twenty-five animals received scanned proton beam irradiation. ECG-gated computed tomography scans were acquired at end-expiration breath hold. Structures (atrioventricular junction or left ventricular myocardium) and organs at risk were contoured. Doses of 30, 40, and 55 Gy were delivered during expiration to the atrioventricular junction (n=5) and left ventricular myocardium (n=20) of intact animals. Results: In this study, procedural success was tracked by pacemaker interrogation in the atrioventricular junction group, time-course magnetic resonance imaging in the left ventricular group, and correlation of lesion outcomes displayed in gross and microscopic pathology. Protein extraction (active caspase-3) was performed to investigate tissue apoptosis. Doses of 40 and 55 Gy caused slowing and interruption of cardiac impulse propagation at the atrioventricular junction. In 40 left ventricular irradiated targets, all lesions were identified on magnetic resonance after 12 weeks, being consistent with outcomes from gross pathology. In the majority of cases, lesion size plateaued between 12 and 16 weeks. Active caspase-3 was seen in lesions 12 and 16 weeks after irradiation but not after 20 weeks. Conclusions: Scanned proton beams can be used as a tool for catheter-free ablation, and time-course of tissue apoptosis was consistent with lesion maturation.

AbstractPurposeThe aim of this work is to describe the clinical implementation of respiratory‐gated spot‐scanning proton therapy (SSPT) for the treatment of thoracic and abdominal moving targets. The experience of our institution is summarized, from initial acceptance and commissioning tests to the development of standard clinical operating procedures for simulation, motion assessment, motion mitigation, treatment planning, and gated SSPT treatment delivery.Materials and methodsA custom respiratory gating interface incorporating the Real‐Time Position Management System (RPM, Varian Medical Systems, Inc., Palo Alto, CA, USA) was developed in‐house for our synchrotron‐based delivery system. To assess gating performance, a motion phantom and radiochromic films were used to compare gated vs nongated delivery. Site‐specific treatment planning protocols and conservative motion cutoffs were developed, allowing for free‐breathing (FB), breath‐holding (BH), or phase‐gating (Ph‐G). Room usage efficiency of BH and Ph‐G treatments was retrospectively evaluated using beam delivery data retrieved from our record and verify system and DICOM files from patient‐specific quality assurance (QA) procedures.ResultsMore than 70 patients were treated using active motion management between the launch of our motion mitigation program in October 2015 and the end date of data collection of this study in January 2018. During acceptance procedures, we found that overall system latency is clinically‐suitable for Ph‐G. Regarding room usage efficiency, the average number of energy layers delivered per minute was <10 for Ph‐G, 10‐15 for BH and ≥15 for FB, making Ph‐G the slowest treatment modality. When comparing to continuous delivery measured during pretreatment QA procedures, the median values of BH treatment time were extended from 6.6 to 9.3 min (+48%). Ph‐G treatments were extended from 7.3 to 13.0 min (+82%).ConclusionsActive motion management has been crucial to the overall success of our SSPT program. Nevertheless, our conservative approach has come with an efficiency cost that is more noticeable in Ph‐G treatments and should be considered in decision‐making.

Purpose Photon counting detectors (PCDs) are being introduced in advanced x‐ray computed tomography (CT) scanners. From a single PCD‐CT acquisition, multiple images can be reconstructed, each based on only a part of the original x‐ray spectrum. In this study, we investigated whether PCD‐CT can be used to estimate stopping power ratios (SPRs) for proton therapy treatment planning, both by comparing to other SPR methods proposed for single energy CT (SECT) and dual energy CT (DECT) as well as to experimental measurements. Methods A previously developed DECT‐based SPR estimation method was adapted to PCD‐CT data, by adjusting the estimation equations to allow for more energy spectra. The method was calibrated directly on noisy data to increase the robustness toward image noise. The new PCD SPR estimation method was tested in theoretical calculations as well as in an experimental setup, using both four and two energy bin PCD‐CT images, and through comparison to two other SPR methods proposed for SECT and DECT. These two methods were also evaluated on PCD‐CT images, full spectrum (one‐bin) or two‐bin images, respectively. In a theoretical framework, we evaluated the effect of patient‐specific tissue variations (density and elemental composition) and image noise on the SPR accuracy; the latter effect was assessed by applying three different noise levels (low, medium, and high noise). SPR estimates derived using real PCD‐CT images were compared to experimentally measured SPRs in nine organic tissue samples, including fat, muscle, and bone tissues. Results For the theoretical calculations, the root‐mean‐square error (RMSE) of the SPR estimation was 0.1% for the new PCD method using both two and four energy bins, compared to 0.2% and 0.7% for the DECT‐ and SECT‐based method, respectively. The PCD method was found to be very robust toward CT image noise, with a RMSE of 2.7% when high noise was added to the CT numbers. Introducing tissue variations, the RMSE only increased to 0.5%; even when adding high image noise to the changed tissues, the RMSE stayed within 3.1%. In the experimental measurements, the RMSE over the nine tissue samples was 0.8% when using two energy bins, and 1.0% for the four‐bin images. Conclusions In all tested cases, the new PCD method produced similar or better results than the SECT‐ and DECT‐based methods, showing an overall improvement of the SPR accuracy. This study thus demonstrated that PCD‐CT scans will be a qualified candidate for SPR estimations.