Treatment planning aspects and Monte Carlo methods in proton therapy

Over the last years, the interest in proton radiotherapy is rapidly increasing. Protons provide superior physical properties compared with conventional radiotherapy using photons. These properties result in depth dose curves with a large dose peak at the end of the proton track and the finite proton range allows sparing the distally located healthy tissue. These properties offer an increased flexibility in proton radiotherapy, but also increase the demand in accurate dose estimations. To carry out accurate dose calculations, first an accurate and detailed characterization of the physical proton beam exiting the treatment head is necessary for both currently available delivery techniques: scattered and scanned proton beams. Since Monte Carlo (MC) methods follow the particle track simulating the interactions from first principles, this technique is perfectly suited to accurately model the treatment head. Nevertheless, careful validation of these MC models is necessary. While for the dose estimation pencil beam algorithms provide the advantage of fast computations, they are limited in accuracy. In contrast, MC dose calculation algorithms overcome these limitations and due to recent improvements in efficiency, these algorithms are expected to improve the accuracy of the calculated dose distributions and to be introduced in clinical routine in the near future.

Hypomagnesemia Induced by Long-Term Treatment with Proton-Pump Inhibitors.

In 2006, hypomagnesemia was first described as a complication of proton-pump inhibitors. To address this issue, we systematically reviewed the literature. Hypomagnesemia, mostly associated with hypocalcemic hypoparathyroidism and hypokalemia, was reported in 64 individuals on long-term proton-pump inhibitors. Hypomagnesemia recurred following replacement of one proton-pump inhibitor with another but not with a histamine type-2 receptor antagonist. The association between proton-pump inhibitors and magnesium metabolism was addressed in 14 case-control, cross-sectional studies. An association was found in 11 of them: 6 reports found that the use of proton-pump inhibitors is associated per se with a tendency towards hypomagnesemia, 2 found that this tendency is more pronounced in patients concurrently treated with diuretics, carboplatin, or cisplatin, and 2 found a relevant tendency to hypomagnesemia in patients with poor renal function. Finally, findings likely reflecting decreased intestinal magnesium uptake were observed on treatment with proton-pump inhibitors. Three studies did not disclose any relationship between magnesium metabolism and treatment with histamine type-2 receptor antagonists. In conclusion, proton-pump inhibitors may cause hypomagnesemia. In these cases, switching to a histamine type-2 receptor antagonist is advised.

Robustness of the Voluntary Breath-Hold Approach for the Treatment of Peripheral Lung Tumors Using Hypofractionated Pencil Beam Scanning Proton Therapy.

PURPOSE The safe clinical implementation of pencil beam scanning (PBS) proton therapy for lung tumors is complicated by the delivery uncertainties caused by breathing motion. The purpose of this feasibility study was to investigate whether a voluntary breath-hold technique could limit the delivery uncertainties resulting from interfractional motion. METHODS AND MATERIALS Data from 15 patients with peripheral lung tumors previously treated with stereotactic radiation therapy were included in this study. The patients had 1 computed tomographic (CT) scan in voluntary breath-hold acquired before treatment and 3 scans during the treatment course. PBS proton treatment plans with 2 fields (2F) and 3 fields (3F), respectively, were calculated based on the planning CT scan and subsequently recalculated on the 3 repeated CT scans. Recalculated plans were considered robust if the V95% (volume receiving ≥95% of the prescribed dose) of the gross target volume (GTV) was within 5% of what was expected from the planning CT data throughout the simulated treatment. RESULTS A total of 14/15 simulated treatments for both 2F and 3F met the robustness criteria. Reduced V95% was associated with baseline shifts (2F, P=.056; 3F, P=.008) and tumor size (2F, P=.025; 3F, P=.025). Smaller tumors with large baseline shifts were also at risk for reduced V95% (interaction term baseline/size: 2F, P=.005; 3F, P=.002). CONCLUSIONS The breath-hold approach is a realistic clinical option for treating lung tumors with PBS proton therapy. Potential risk factors for reduced V95% are small targets in combination with large baseline shifts. On the basis of these results, the baseline shift of the tumor should be monitored (eg, through image guided therapy), and appropriate measures should be taken accordingly. The intrafractional motion needs to be investigated to confirm that the breath-hold approach is robust.

Dose-painting intensity-modulated proton therapy for intermediate- and high-risk meningioma.

BACKGROUND Newly diagnosed WHO grade II-III or any WHO grade recurrent meningioma exhibit an aggressive behavior and thus are considered as high- or intermediate risk tumors. Given the unsatisfactory rates of disease control and survival after primary or adjuvant radiation therapy, optimization of treatment strategies is needed. We investigated the potential of dose-painting intensity-modulated proton beam-therapy (IMPT) for intermediate- and high-risk meningioma. MATERIAL AND METHODS Imaging data from five patients undergoing proton beam-therapy were used. The dose-painting target was defined using [68]Ga-[1,4,7,10-tetraazacyclododecane tetraacetic acid]- d-Phe(1),Tyr(3)-octreotate ([68]Ga-DOTATATE)-positron emission tomography (PET) in target delineation. IMPT and photon intensity-modulated radiation therapy (IMRT) treatment plans were generated for each patient using an in-house developed treatment planning system (TPS) supporting spot-scanning technology and a commercial TPS, respectively. Doses of 66 Gy (2.2 Gy/fraction) and 54 Gy (1.8 Gy/fraction) were prescribed to the PET-based planning target volume (PTVPET) and the union of PET- and anatomical imaging-based PTV, respectively, in 30 fractions, using simultaneous integrated boost. RESULTS Dose coverage of the PTVsPET was equally good or slightly better in IMPT plans: dose inhomogeneity was 10 ± 3% in the IMPT plans vs. 13 ± 1% in the IMRT plans (p = 0.33). The brain Dmean and brainstem D50 were small in the IMPT plans: 26.5 ± 1.5 Gy(RBE) and 0.002 ± 0.0 Gy(RBE), respectively, vs. 29.5 ± 1.5 Gy (p = 0.001) and 7.5 ± 11.1 Gy (p = 0.02) for the IMRT plans, respectively. The doses delivered to the optic structures were also decreased with IMPT. CONCLUSIONS Dose-painting IMPT is technically feasible using currently available planning tools and resulted in dose conformity of the dose-painted target comparable to IMRT with a significant reduction of radiation dose delivered to the brain, brainstem and optic apparatus. Dose escalation with IMPT may improve tumor control and decrease radiation-induced toxicity.

Intrinsic characteristics of the proton pump in the luminal membrane of a tight urinary epithelium. The relation between transport rate and delta mu H.

A number of tight urinary epithelia, as exemplified by the turtle bladder, acidify the luminal solution by active transport of H+ across the luminal cell membrane. The rate of active H+ transport (JH) decreases as the electrochemical potential difference for H+ [delta mu H = mu H(lumen) - mu H(serosa)] across the epithelium is increased. The luminal cell membrane has a low permeability for H+ equivalents and a high electrical resistance compared with the basolateral cell membrane. Changes in JH thus reflect changes in active H+ transport across the luminal membrane. To examine the control of JH by delta mu H in the turtle bladder, transepithelial electrical potential differences (delta psi) were imposed at constant acid-base conditions or the luminal pH was varied at delta psi = 0 and constant serosal PCO2 and pH. When the luminal compartment was acidified from pH 7 to 4 or was made electrically positive, JH decreased as a linear function of delta mu H as previously described. When the luminal compartment was made alkaline from pH 7 to 9 or was made electrically negative, JH reached a maximal value, which was the same whether the delta mu H was imposed as a delta pH or a delta psi. The nonlinear JH vs. delta mu H relation does not result from changes in the number of pumps in the luminal membrane or from changes in the intracellular pH, but is a characteristic of the H+ pumps themselves. We propose a general scheme, which, because of its structural features, can account for the nonlinearity of the JH vs. delta mu H relations and, more specifically, for the kinetic equivalence of the effects of the chemical and electrical components of delta mu H. According to this model, the pump complex consists of two components: a catalytic unit at the cytoplasmic side of the luminal membrane, which mediates the ATP-driven H+ translocation, and a transmembrane channel, which mediates the transfer of H+ from the catalytic unit to the luminal solution. These two components may be linked through a buffer compartment for H+ (an antechamber).

Monte Carlo and analytical dose calculations for ocular proton therapy

Uveal melanoma is a rare but life-threatening form of ocular cancer. Contemporary treatment techniques include proton therapy, which enables conservation of the eye and its useful vision. Dose to the proximal structures is widely believed to play a role in treatment side effects, therefore, reliable dose estimates are required for properly evaluating the therapeutic value and complication risk of treatment plans. Unfortunately, current simplistic dose calculation algorithms can result in errors of up to 30% in the proximal region. In addition, they lack predictive methods for absolute dose per monitor unit (D/MU) values. ^ To facilitate more accurate dose predictions, a Monte Carlo model of an ocular proton nozzle was created and benchmarked against measured dose profiles to within ±3% or ±0.5 mm and D/MU values to within ±3%. The benchmarked Monte Carlo model was used to develop and validate a new broad beam dose algorithm that included the influence of edgescattered protons on the cross-field intensity profile, the effect of energy straggling in the distal portion of poly-energetic beams, and the proton fluence loss as a function of residual range. Generally, the analytical algorithm predicted relative dose distributions that were within ±3% or ±0.5 mm and absolute D/MU values that were within ±3% of Monte Carlo calculations. Slightly larger dose differences were observed at depths less than 7 mm, an effect attributed to the dose contributions of edge-scattered protons. Additional comparisons of Monte Carlo and broad beam dose predictions were made in a detailed eye model developed in this work, with generally similar findings. ^ Monte Carlo was shown to be an excellent predictor of the measured dose profiles and D/MU values and a valuable tool for developing and validating a broad beam dose algorithm for ocular proton therapy. The more detailed physics modeling by the Monte Carlo and broad beam dose algorithms represent an improvement in the accuracy of relative dose predictions over current techniques, and they provide absolute dose predictions. It is anticipated these improvements can be used to develop treatment strategies that reduce the incidence or severity of treatment complications by sparing normal tissue. ^

Proton therapy versus intensity modulated x-ray therapy in the treatment of prostate cancer: Estimating secondary cancer risks

External beam radiation therapy is used to treat nearly half of the more than 200,000 new cases of prostate cancer diagnosed in the United States each year. During a radiation therapy treatment, healthy tissues in the path of the therapeutic beam are exposed to high doses. In addition, the whole body is exposed to a low-dose bath of unwanted scatter radiation from the pelvis and leakage radiation from the treatment unit. As a result, survivors of radiation therapy for prostate cancer face an elevated risk of developing a radiogenic second cancer. Recently, proton therapy has been shown to reduce the dose delivered by the therapeutic beam to normal tissues during treatment compared to intensity modulated x-ray therapy (IMXT, the current standard of care). However, the magnitude of stray radiation doses from proton therapy, and their impact on this incidence of radiogenic second cancers, was not known. ^ The risk of a radiogenic second cancer following proton therapy for prostate cancer relative to IMXT was determined for 3 patients of large, median, and small anatomical stature. Doses delivered to healthy tissues from the therapeutic beam were obtained from treatment planning system calculations. Stray doses from IMXT were taken from the literature, while stray doses from proton therapy were simulated using a Monte Carlo model of a passive scattering treatment unit and an anthropomorphic phantom. Baseline risk models were taken from the Biological Effects of Ionizing Radiation VII report. A sensitivity analysis was conducted to characterize the uncertainty of risk calculations to uncertainties in the risk model, the relative biological effectiveness (RBE) of neutrons for carcinogenesis, and inter-patient anatomical variations. ^ The risk projections revealed that proton therapy carries a lower risk for radiogenic second cancer incidence following prostate irradiation compared to IMXT. The sensitivity analysis revealed that the results of the risk analysis depended only weakly on uncertainties in the risk model and inter-patient variations. Second cancer risks were sensitive to changes in the RBE of neutrons. However, the findings of the study were qualitatively consistent for all patient sizes and risk models considered, and for all neutron RBE values less than 100. ^

Quantitative comparison of late effects following photon versus proton external-beam radiation therapies: Toward an evidence-based approach to selecting a treatment modality

Radiation therapy has been used as an effective treatment for malignancies in pediatric patients. However, in many cases, the side effects of radiation diminish these patients’ quality of life. In order to develop strategies to minimize radiogenic complications, one must first quantitatively estimate pediatric patients’ relative risk for radiogenic late effects, which has not become feasible till recently because of the calculational complexity. The goals of this work were to calculate the dose delivered to tissues and organs in pediatric patients during contemporary photon and proton radiotherapies; to estimate the corresponding risk of radiogenic second cancer and cardiac toxicity based on the calculated doses and on dose-risk models from the literature; to test for the statistical significance of the difference between predicted risks after photon versus proton radiotherapies; and to provide a prototype of an evidence-based approach to selecting treatment modalities for pediatric patients, taking second cancer and cardiac toxicity into account. The results showed that proton therapy confers a lower predicted risk of radiogenic second cancer, and lower risks of radiogenic cardiac toxicities, compared to photon therapy. An uncertainty analysis revealed that the qualitative findings of this study are insensitive to changes in a wide variety of host and treatment related factors.

Verification of the clinical implementation of the respiratory gated beam delivery technique with synchrotron-based proton irradiation

The clinical advantage for protons over conventional high-energy x-rays stems from their unique depth-dose distribution, which delivers essentially no dose beyond the end of range. In order to achieve it, accurate localization of the tumor volume relative to the proton beam is necessary. For cases where the tumor moves with respiration, the resultant dose distribution is sensitive to such motion. One way to reduce uncertainty caused by respiratory motion is to use gated beam delivery. The main goal of this dissertation is to evaluate the respiratory gating technique in both passive scattering and scanning delivery mode. Our hypothesis for the study was that optimization of the parameters of synchrotron operation and respiratory gating can lead to greater efficiency and accuracy of respiratory gating for all modes of synchrotron-based proton treatment delivery. The hypothesis is tested in two specific aims. The specific aim #1 is to assess the efficiency of respiratory-gated proton beam delivery and optimize the synchrotron operations for the gated proton therapy. A simulation study was performed and introduced an efficient synchrotron operation pattern, called variable Tcyc. In addition, the simulation study estimated the efficiency in the respiratory gated scanning beam delivery mode as well. The specific aim #2 is to assess the accuracy of beam delivery in respiratory-gated proton therapy. The simulation study was extended to the passive scattering mode to estimate the quality of pulsed beam delivery to the residual motion for several synchrotron operation patterns with the gating technique. The results showed that variable Tcyc operation can offer good reproducible beam delivery to the residual motion at a certain phase of the motion. For respiratory gated scanning beam delivery, the impact of motion on the dose distributions by scanned beams was investigated by measurement. The results showed the threshold for motion for a variety of scan patterns and the proper number of paintings for normal and respiratory gated beam deliveries. The results of specific aims 1 and 2 provided supporting data for implementation of the respiratory gating beam delivery technique into both passive and scanning modes and the validation of the hypothesis.

Evaluation of Polymer Gel Dosimeters for Measurements of Dose and LET in Proton Beams

This project assessed the effectiveness of polymer gel dosimeters as tools for measuring the dose deposited by and LET of a proton beam. A total of three BANG® dosimeter formulations were evaluated: BANG®-3-Pro-2 BANGkits™ for dose measurement and two BANG®-3 variants, the LET-Baseline and LET-Meter dosimeters, for LET measurement. All dosimeters were read out using an OCT scanner. The basic characteristics of the BANGkits™ were assessed in a series of photon and electron irradiations. The dose-response relationship was found to be sigmoidal with a threshold for response of approximately 15 cGy. The active region of the dosimeter, the volume in which dosimeter response is not inhibited by oxygen, was found to make up roughly one fourth of the total dosimeter volume. Delivering a dose across multiple fractions was found to yield a greater response than delivering the same dose in a single irradiation. The dosimeter was found to accurately measure a dose distribution produced by overlapping photon fields, yielding gamma pass rates of 95.4% and 93.1% from two planar gamma analyses. Proton irradiations were performed for measurements of proton dose and LET. Initial irradiations performed through the side of a dosimeter led to OCT artifacts. Gamma pass rates of 85.7% and 89.9% were observed in two planar gamma analyses. In irradiations performed through the base of a dosimeter, gel response was found to increase with height in the dosimeter, even in areas of constant dose. After a correction was applied, gamma pass rates of 94.6% and 99.3% were observed in two planar gamma analyses. Absolute dose measurements were substantially higher (33%-100%) than the delivered doses for proton irradiations. Issues encountered while calibrating the LET-Meter gel restricted analysis of the LET measurement data to the SOBP of a proton beam. LET-Meter overresponse was found to increase linearly with track-average LET across the LET range that could be investigated (1.5 keV/micron – 3.5 keV/micron).

DEVELOPMENT OF A BEAM-SPECIFIC PLANNING TARGET VOLUME AND A ROBUST PLAN ANALYSIS TOOL FOR PROTON THERAPY

Proton therapy is growing increasingly popular due to its superior dose characteristics compared to conventional photon therapy. Protons travel a finite range in the patient body and stop, thereby delivering no dose beyond their range. However, because the range of a proton beam is heavily dependent on the tissue density along its beam path, uncertainties in patient setup position and inherent range calculation can degrade thedose distribution significantly. Despite these challenges that are unique to proton therapy, current management of the uncertainties during treatment planning of proton therapy has been similar to that of conventional photon therapy. The goal of this dissertation research was to develop a treatment planning method and a planevaluation method that address proton-specific issues regarding setup and range uncertainties. Treatment plan designing method adapted to proton therapy: Currently, for proton therapy using a scanning beam delivery system, setup uncertainties are largely accounted for by geometrically expanding a clinical target volume (CTV) to a planning target volume (PTV). However, a PTV alone cannot adequately account for range uncertainties coupled to misaligned patient anatomy in the beam path since it does not account for the change in tissue density. In order to remedy this problem, we proposed a beam-specific PTV (bsPTV) that accounts for the change in tissue density along the beam path due to the uncertainties. Our proposed method was successfully implemented, and its superiority over the conventional PTV was shown through a controlled experiment.. Furthermore, we have shown that the bsPTV concept can be incorporated into beam angle optimization for better target coverage and normal tissue sparing for a selected lung cancer patient. Treatment plan evaluation method adapted to proton therapy: The dose-volume histogram of the clinical target volume (CTV) or any other volumes of interest at the time of planning does not represent the most probable dosimetric outcome of a given plan as it does not include the uncertainties mentioned earlier. Currently, the PTV is used as a surrogate of the CTV’s worst case scenario for target dose estimation. However, because proton dose distributions are subject to change under these uncertainties, the validity of the PTV analysis method is questionable. In order to remedy this problem, we proposed the use of statistical parameters to quantify uncertainties on both the dose-volume histogram and dose distribution directly. The robust plan analysis tool was successfully implemented to compute both the expectation value and its standard deviation of dosimetric parameters of a treatment plan under the uncertainties. For 15 lung cancer patients, the proposed method was used to quantify the dosimetric difference between the nominal situation and its expected value under the uncertainties.

Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy

DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY by James Leroy Neihart, B.S. APPROVED: ______________________________David Followill, Ph.D. ______________________________Peter Balter, Ph.D. ______________________________Narayan Sahoo, Ph.D. ______________________________Kenneth Hess, Ph.D. ______________________________Paige Summers, M.S. APPROVED: ____________________________ Dean, The University of Texas Graduate School of Biomedical Sciences at Houston DEVELOPMENT AND IMPLEMENTATION OF A DYNAMIC HETEROGENEOUS PROTON EQUIVALENT ANTHROPOMORPHIC THORAX PHANTOM FOR THE ASSESSMENT OF SCANNED PROTON BEAM THERAPY A THESIS Presented to the Faculty of The University of Texas Health Science Center at Houston andThe University of TexasMD Anderson Cancer CenterGraduate School of Biomedical Sciences in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE by James Leroy Neihart, B.S. Houston, Texas Date of Graduation August, 2013 Acknowledgments I would like to acknowledge my advisory committee members, chair David Followill, Ph.D., Peter Balter, Ph.D, Narayan Sahoo, Ph.D., Kenneth Hess, Ph.D., Paige Summers M.S. and, for their time and effort contributed to this project. I would additionally like to thank the faculty and staff at the PTC-H and the RPC who assisted in many aspects of this project. Falk Pӧnisch, Ph.D. for his breath hold proton therapy treatment expertise, Matt Palmer and Jaques Bluett for proton dosimetry assistance, Matt Kerr for verification plan assistance, Carrie Amador, Nadia Hernandez, Trang Nguyen, Andrea Molineu, Lynda McDonald for TLD and film dosimetry assistance. Finally, I would like to thank my wife and family for their support and encouragement during my research and studies. Development and implementation of a dynamic heterogeneous proton equivalent anthropomorphic thorax phantom for the assessment of scanned proton beam therapy By: James Leroy Neihart, B.S. Chair of Advisory Committee: David Followill, Ph.D Proton therapy has been gaining ground recently in radiation oncology. To date, the most successful utilization of proton therapy is in head and neck cases as well as prostate cases. These tumor locations do not suffer from the resulting difficulties of treatment delivery as a result of respiratory motion. Lung tumors require either breath hold or motion tracking, neither of which have been assessed with an end-to-end phantom for proton treatments. Currently, the RPC does not have a dynamic thoracic phantom for proton therapy procedure assessment. Additionally, such a phantom could be an excellent means of assessing quality assurance of the procedures of proton therapy centers wishing to participate in clinical trials. An eventual goal of this phantom is to have a means of evaluating and auditing institutions for the ability to start clinical trials utilizing proton therapy procedures for lung cancers. Therefore, the hypothesis of this study is that a dynamic anthropomorphic thoracic phantom can be created to evaluate end-to-end proton therapy treatment procedures for lung cancer to assure agreement between the measured and calculated dose within 5% / 5 mm with a reproducibility of 2%. Multiple materials were assessed for thoracic heterogeneity equivalency. The phantom was designed from the materials found to be in greatest agreement. The phantom was treated in an end-to-end treatment four times, which included simulation, treatment planning and treatment delivery. Each treatment plan was delivered three times to assess reproducibility. The dose measured within the phantom was compared to that of the treatment plan. The hypothesis was fully supported for three of the treatment plans, but failed the reproducibility requirement for the most aggressive treatment plan.

ATPase activity, proton gradients and proton secretion inhibition during experiments with Sepia officinalis and Sepioteuthis lessoniana

The constraints of an active life in a pelagic habitat led to numerous convergent morphological and physiological adaptations that enable cephalopod molluscs and teleost fishes to compete for similar resources. Here, we show for the first time that such convergent developments are also found in the ontogenetic progression of ion regulatory tissues; as in teleost fish, epidermal ionocytes scattered on skin and yolk sac of cephalopod embryos appear to be responsible for ionic and acid-base regulation before gill epithelia become functional. Ion and acid-base regulation is crucial in cephalopod embryos, as they are surrounded by a hypercapnic egg fluid with a Pco2 between 0.2 and 0.4 kPa. Epidermal ionocytes were characterized via immunohistochemistry, in situ hybridization, and vital dye-staining techniques. We found one group of cells that is recognized by concavalin A and MitoTracker, which also expresses Na+/H+ exchangers (NHE3) and Na+-K+-ATPase. Similar to findings obtained in teleosts, these NHE3-rich cells take up sodium in exchange for protons, illustrating the energetic superiority of NHE-based proton excretion in marine systems. In vivo electrophysiological techniques demonstrated that acid equivalents are secreted by the yolk and skin integument. Intriguingly, epidermal ionocytes of cephalopod embryos are ciliated as demonstrated by scanning electron microscopy, suggesting a dual function of epithelial cells in water convection and ion regulation. These findings add significant knowledge to our mechanistic understanding of hypercapnia tolerance in marine organisms, as it demonstrates that marine taxa, which were identified as powerful acid-base regulators during hypercapnic challenges, already exhibit strong acid-base regulatory abilities during embryogenesis.