Protonentherapie ist eine wichtige Behandlungsmodalität in der Radioonkologie. Aufgrund einer vorteilhaften Dosisverteilung im bestrahlten Volumen kann diese Bestrahlungsmethode das tumorumgebende Normalgewebe schützen. Dadurch können Nebenwirkungen in bestimmten Patientenpopulationen, zum Beispiel Kindern oder Patienten mit Gehirntumoren, verringert werden. Trotzdem können nach Protonenbestrahlung von Gehirntumorpatienten Normalgewebsschäden auftreten. Gründe dafür können der notwendige klinische Sicherheitssaum im Normalgewebe, der Einfluss der relativen biologischen Wirksamkeit RBE sowie eine erhöhte Strahlensensitivität bestimmter Gehirnregionen sein. Um diese Aspekte zu beleuchten, werden geeignete präklinische Modelle für die Normalgewebsreaktion im Gehirn nach Protonenbestrahlung benötigt. Darüber hinaus kann eine Risikostratifizierung der Patienten durch die Vorhersage von Nebenwirkungswahrscheinlichkeiten oder der Tumorantwort den Behandlungserfolg erhöhen. Auch hier können präklinische Modelle helfen, um neue prädiktive Biomarker zu finden und um die zugrunde liegenden Mechanismen strahleninduzierter Gehirnschäden besser zu verstehen. Das Ziel dieser Dissertation war die Etablierung und Charakterisierung von adäquaten präklinischen Modellen für die Untersuchung von strahleninduzierten Normalgewebsschäden im Gehirn. Diese Modelle bilden die Grundlage für zukünftige Studien zur Untersuchung von RBE Effekten, der spezifische Strahlensensitivität einzelner Gehirnregionen und neuer Biomarker. Die getesteten Modellsysteme waren in vitro Kulturen von adulten organotypischen Gehirnschnitten, Tumorschnittkultur sowie in vivo Bestrahlung von Gehirnsubvolumina, jeweils mit dem Modellorganismus Maus. Die Etablierung eines Bestrahlungssetups in der experimentellen Protonenanlage und dessen dosimetrische Charakterisierung waren von großer Bedeutung für die Durchführung der biologischen Experimente. Ein weiteres Hauptziel war die Definition klinisch relevanter Endpunkte für frühe und späte Nebenwirkungen. Die Gewebsschnitte wurden durch Messungen des Zellüberlebens und der Entzündungsreaktion, sowie mittels in situ Analyse von Zellmorphologie und DNA Schäden untersucht. Als ergänzendes Modell wurde die Tumorschnittkultur etabliert und ähnliche Endpunkte analysiert. Adulte Gehirnschnitte stellten sich als ungeeignet für präklinische Experimente in der Radioonkologie heraus. Die Messungen von Zelltod und Entzündungswerten zeigten eine starke Zellreaktion auf die Inkulturnahme, aber keine auf die Protonenbestrahlung. In der Histologie wurden gestörte Zellmorphologie, reduzierte Vitalität und eingeschränkte Reparaturfähigkeit von DNA Schäden beobachtet. Daher sollten für strahlenbiologische Experimente andere 3D Zellkulturmodelle in Betracht gezogen werden, wie zum Beispiel Organoide oder durch Tissue Engineering hergestellte Kulturen. Durch die Publikation der Daten leistet diese Dissertation einen wichtigen Beitrag zur aktuellen Forschung, da so künftig die limitierten Ressourcen, die für strahlenbiologische Experimente mit Protonen zur Verfügung stehen, auf relevantere Modelle verwendet werden können. Die Bestrahlung von Gehirnsubvolumina in Mäusen wurde mit dem Ziel etabliert, klinisch vergleichbare Felder zu erreichen. Das gewählte Zielvolumen war der rechte Hippocampus; der Protonenstrahl sollte in der Mitte des Gehirns stoppen. Im Rahmen des Projekts wurde ein Arbeitsablauf für präzise und reproduzierbare Bestrahlung entwickelt. Zur Verifizierung wurde der induzierte DNA Schaden ausgewertet und anschließend mit Monte-Carlos Dosissimulationen korreliert. Die Maushirnbestrahlung lieferte wertvolle Ergebnisse für frühe Zeitpunkte (d.h. innerhalb 24 h nach Bestrahlung). Im Verlauf des Projekts wurde ein Algorithmus erstellt, der schnell und zuverlässig die räumliche Verteilung des DNA Schadens in Relation zur Gesamtzellzahl analysiert. Diese Auswertung zeigte, wie bei der Bestrahlungsplanung vorgesehen, ein Stoppen des Protonenstrahls im Gehirn. Eine anschließende Korrelation der Schadensverteilung mit der applizierten Dosis weist nach, dass das Modell einen wichtigen Beitrag zur Untersuchung des RBE leisten kann. In einer darauf folgenden Studie wurde der Dosis-Zeitverlauf der beobachteten Strahlenreaktion des Normalgewebes genauer beleuchtet. Dafür wurden Untersuchungen des Allgemeinzustands der Versuchstiere, regelmäßige Magnetresonanztomografie (MRI) Messungen über einen Zeitraum von sechs Monaten, sowie abschließende Histologie korreliert. Die Volumenzunahme des Kontrastmittelaustritts, die den Zusammenbruch der Blut-Hirn-Schranke anzeigt, wurde konturiert; aus diesen Daten entstand ein prädiktives Dosis-Volumen Modell. Die Pilotstudie konnte eine dosisabhängige Strahlenreaktion nachweisen, die sich im Zusammenbruch der Blut-Hirn-Schranke, einer Hautreaktion mit vorrübergehender Alopezie, Gewichtsabnahme und zelluläre Veränderung äußerte. Das von den MRI Messungen abgeleitete Modell konnte zuverlässig das Eintreten der Nebenwirkungen, den Krankheitsverlauf, sowie die geschätzte Überlebensdauer der Mäuse vorhersagen. Zusätzlich konnte ein Zusammenhang zwischen den MRI Bildänderungen und den pathologischen Gewebsveränderungen beobachtet werden. Durch die außerordentlich homogene Strahlenreaktion der Tiere können aus den vorliegenden Daten künftig zuverlässig geeignete Dosen für spezifische experimentelle Endpunkte bestimmt werden. Zusammenfassend wurden in dieser Arbeit zwei präklinische Modelle für die Protonengehirnbestrahlung etabliert, nämlich organotypische Gewebsschnitte als 3D Zellkulturmodell sowie in vivo Bestrahlung von Gehirnsubvolumina in Mäusen. Während Zellkulturexperimente die Erwartungen nicht erfüllen konnten, stellen sich die Tierexperimente als hervorragendes Modell für translationale Radioonkologie heraus, welches zusätzlich für andere Strahlenqualitäten eingesetzt werden kann. Darauf basierend können aktuelle und zukünftige Studien die Ursachen von strahleninduzierten Normalgewebsschäden im Gehirn beleuchten, RBE Effekte untersuchen und neue prädiktive Biomarker erforschen.:Contents Abstract i Zusammenfassung v Publications ix List of Figures xiii List of Acronyms and Abbreviations xiv 1 Introduction 3 2 Background 5 2.1 Proton therapy for brain cancer treatment 5 2.1.1 Fundamentals of radiobiology 5 2.1.2 Proton therapy 6 2.1.3 Tumors of the central nervous system 8 2.2 Radiation effects on brain cells 8 2.2.1 Neurons and myelin 9 2.2.2 Blood-brain barrier 9 2.2.3 Astrocytes 10 2.2.4 Microglia 10 2.3 Principles of histology 11 2.3.1 Hematoxylin & eosin staining 12 2.3.2 Immunohistochemistry 13 2.3.3 Bioimage analysis 13 2.4 Techniques in medical imaging 14 2.4.1 Projectional radiography 14 2.4.2 Computed tomography 14 2.4.3 Magnetic resonance imaging 15 2.5 Preclinical models for radiation injury 17 2.5.1 Technical requirements 17 2.5.2 In vitro models 17 2.5.3 Small animal models 18 3 Applying Tissue Slice Culture in Cancer Research – Insights from Preclinical Proton Radiotherapy 19 3.1 Aim of the study 19 3.2 Conclusion 19 3.3 Author’s contribution 19 3.4 Publication 21 4 High-precision image-guided proton irradiation of mouse brain sub-volumes 41 4.1 Aim of the study 41 4.2 Conclusion 41 4.3 Author’s contribution 41 4.4 Publication 43 5 Late side effects in normal mouse brain tissue after proton irradiation 51 5.1 Aim of the study 51 5.2 Conclusion 51 5.3 Author’s contribution 52 5.4 Publication 53 6 Discussion 71 6.1 Establishment of preclinical models for radiooncology 71 6.1.1 3D cell culture 71 6.1.2 In vivo irradiation of brain subvolumes 73 6.2 Current applications of the mouse model 75 6.2.1 Ongoing data analysis 75 6.2.2 Innovating on-site imaging 76 6.2.3 RBE investigations 77 6.3 Future studies of radiation-induced brain tissue toxicities 79 Acknowledgement XV Supplementary Material XVII 1 Applying Tissue Slice Culture in Cancer Research – Insights from Preclinical Proton Radiotherapy XVII 2 High-precision image-guided proton irradiation of mouse brain sub-volumes XXVI 3 Late side effects in normal mouse brain tissue after proton irradiation XXXI / Proton therapy is an important modality in radiation oncology. Due to a favorable dose distribution in the irradiated volume, this treatment allows to spare tumor-surrounding normal tissue. Although this protection can lead to reduced side effects in certain patient populations, such as brain tumor or pediatric patients, normal tissue toxicities can occur to some extend. This could be due to clinical safety margins around the tumor that lead to dose deposition in the normal tissue. The underlying causes might also be related to relative biological effectiveness (RBE) variations or elevated radiosensitivity of certain brain regions. To address these issues, suitable preclinical models for normal brain tissue reaction after proton therapy are needed. In addition, patient stratification to predict the tumor response or the probability of side effects will contribute to increased treatment effectiveness. Preclinical models can improve the process of finding new predictive biomarkers and help to understand underlying mechanisms of radiation-induced brain injury. The aim of this thesis was to establish and characterize suitable preclinical models of brain tissue irradiation effects and set the base for future studies designed to reveal RBE effects, brain region specific radiation sensitivities, and novel biomarkers. The tested model systems were in vitro organotypic brain slice culture (OBSC) and in vivo irradiation of brain subvolumes, both on mouse brain tissue. Setup establishment at the experimental proton beam line and subsequent dosimetry built the foundation for conducting the biological experiments. Additionally, one main goal was defining clinically relevant endpoints for both short- and long-term effects. For OBSC, assays for cell death and inflammation, as well as in situ analysis of cell morphology and DNA damage induction were tested. As comparative model to OBSC, tumor slice culture was established and the results were also used for proton investigation. Adult OBSC turned out as inadequate model for preclinical experiments in radiation oncology. The assays measuring cell death and inflammation indicated a severe reaction during the first days in culture, but no response to irradiation. Histology revealed deficient cell morphology, reduced vitality and impaired DNA damage repair. In conclusion, other 3D cell culture models, such as organoids or tissue engineered constructs, should be considered for radiobiological experiments with protons. By publishing the observations, this thesis contributes to conserving the limited resources of proton radiobiology for more meaningful models. A methodology for irradiation of mouse brain subvolumes was established with a focus on creating fields comparable to clinical practice. The chosen target was the right hippocampus and the goal was to stop the proton beam in the middle of the brain. The project included a workflow for this precise irradiation in a robust and reproducible manner. Evaluation of the induced DNA damage and its correlation to Monte Carlo dose simulations were used for verification. Irradiation of mouse brain subvolumes yielded valuable results for early (i.e. within 24 h after irradiation) time points. An evaluation algorithm was designed for fast and robust analysis of spatial DNA damage distribution in relation to the total cell count. This ratio showed that the beam stopped in the brain tissue, in accordance to the treatment planning. Furthermore, the DNA damage could be reliably correlated with the dose simulation, which proves the value of the presented model for future RBE studies. In a follow-up experiment, the dose-time relationship of induced normal tissue reactions was analysed. For this, scoring of the animals' health status was combined with regular MRI measurements over the course of up to 6 months, and final histopathology. The volume increase of contrast agent leakage - representing breakdown of the blood brain barrier (BBB) - was contoured and the data was used to create a dose-volume response model. This pilot study on long-term radiation effects revealed dose-dependent normal tissue toxicities, including breakdown of the BBB, a skin reaction with temporary alopecia, weight reduction and changes on the cellular level. The model derived from MRI data reliably predicts onset of side effects, volume of brain damage as well as the expected animal survival. In addition, MRI image changes could be correlated to underlying tissue alterations by histopathology. Due to the uniform radiation response of the animals this data set enables to determine endpoint-specific dose values in future experiments. In conclusion, two preclinical models for proton brain irradiation were established, namely OBSC as 3D cell culture model and in vivo irradiation of mouse brain subvolumes. While the former could not yield the anticipated results, the latter emerged as excellent model for translational radiooncology, which can also be applied for experiments with other radiation types. Ongoing and future studies will focus on revealing the causes of normal brain tissue toxicities, studying RBE effects, and investigating new predictive biomarkers.:Contents Abstract i Zusammenfassung v Publications ix List of Figures xiii List of Acronyms and Abbreviations xiv 1 Introduction 3 2 Background 5 2.1 Proton therapy for brain cancer treatment 5 2.1.1 Fundamentals of radiobiology 5 2.1.2 Proton therapy 6 2.1.3 Tumors of the central nervous system 8 2.2 Radiation effects on brain cells 8 2.2.1 Neurons and myelin 9 2.2.2 Blood-brain barrier 9 2.2.3 Astrocytes 10 2.2.4 Microglia 10 2.3 Principles of histology 11 2.3.1 Hematoxylin & eosin staining 12 2.3.2 Immunohistochemistry 13 2.3.3 Bioimage analysis 13 2.4 Techniques in medical imaging 14 2.4.1 Projectional radiography 14 2.4.2 Computed tomography 14 2.4.3 Magnetic resonance imaging 15 2.5 Preclinical models for radiation injury 17 2.5.1 Technical requirements 17 2.5.2 In vitro models 17 2.5.3 Small animal models 18 3 Applying Tissue Slice Culture in Cancer Research – Insights from Preclinical Proton Radiotherapy 19 3.1 Aim of the study 19 3.2 Conclusion 19 3.3 Author’s contribution 19 3.4 Publication 21 4 High-precision image-guided proton irradiation of mouse brain sub-volumes 41 4.1 Aim of the study 41 4.2 Conclusion 41 4.3 Author’s contribution 41 4.4 Publication 43 5 Late side effects in normal mouse brain tissue after proton irradiation 51 5.1 Aim of the study 51 5.2 Conclusion 51 5.3 Author’s contribution 52 5.4 Publication 53 6 Discussion 71 6.1 Establishment of preclinical models for radiooncology 71 6.1.1 3D cell culture 71 6.1.2 In vivo irradiation of brain subvolumes 73 6.2 Current applications of the mouse model 75 6.2.1 Ongoing data analysis 75 6.2.2 Innovating on-site imaging 76 6.2.3 RBE investigations 77 6.3 Future studies of radiation-induced brain tissue toxicities 79 Acknowledgement XV Supplementary Material XVII 1 Applying Tissue Slice Culture in Cancer Research – Insights from Preclinical Proton Radiotherapy XVII 2 High-precision image-guided proton irradiation of mouse brain sub-volumes XXVI 3 Late side effects in normal mouse brain tissue after proton irradiation XXXI

Seit 1998 werden am Helmholtz-Zentrum Berlin für Materialien und Energie GmbH (HZB) in Kooperation mit der Charité-Universitätsmedizin Berlin Augentumore mit Protonen bestrahlt. Dazu beschleunigt ein Zyklotron, welches hauptsächlich für die Therapie genutzt wird, die Protonen auf eine Energie von 68MeV. Neben dem Therapieplatz existiert auch ein Experimentierplatz zu Forschungszwecken. Die Strahlaufbereitung erfolgt in beiden Anwendungsgebieten primär durch passives Aufstreuen des Strahls mit einer einzelnen Streufolie. Dies führt jedoch zu erheblichen Strahlverlusten, die die erreichbare Strahlintensität bzw. Dosisleistung im Feld limitieren. Ziel dieser Arbeit ist es, einen Bestrahlungsplatz mit aktiver, dreidimensionaler Strahlführung aufzubauen, um eine möglichst effiziente Strahlstromnutzung zu ermöglichen und flexibel auf Anforderungen an das Bestrahlungsfeld eingehen zu können. Zu diesem Zweck wurde ein Modell entwickelt, welches unter Berücksichtigung verschiedener Parameter der Messgeometrie, Strahlcharakteristik und Felderzeugung Strahlprofile berechnet. Dieses Modell wurde im Vergleich mit Monte-Carlo Simulationen und Messungen getestet und verifiziert. Anschließend konnten Untersuchungen zur Positionierung der verwendeten Komponenten durchgeführt sowie deren Einflüsse auf Penumbra, Homogenität und Transmission der Strahlfelder analysiert werden. Es zeigte sich, dass die hohen klinischen Anforderungen für die Anwendung in der Augentumortherapie am HZB erfüllt werden können, sofern stärkere Magnetspulen verwendet werden. Mit Hilfe des Modells wurden für einen optimierten Aufbau Strahlprofile berechnet und charakterisiert. Es konnte so gezeigt werden, dass unter diesen optimierten Bedingungen eine Anwendbarkeit des aktiven Systems in der Augentumortherapie möglich ist. Des Weiteren wurde im Verlauf dieser Arbeit ein statischer Modulator für den Bestrahlungsplatz am HZB entwickelt, mit Hilfe eines 3D-Druckers hergestellt und erfolgreich getestet. / Since 1998, ocular tumours have been irradiated with protons at the Helmholtz- Zentrum Berlin für Materialien und Energie GmbH (HZB) in cooperation with Charité-Universitätsmedizin Berlin. For this purpose a cyclotron is operated which accelerates protons to an energy of 68MeV and is mainly used for therapy purposes. In addition to the therapy station there is also an experimental station for research purposes. In both application areas, beam delivery is primarily achieved by passive scattering using a single scattering foil. However, this leads to considerable beam losses, which in turn limit the maximum beam intensity or dose rate in the irradiation field. The aim of this work is to construct an irradiation setup with active three-dimensional beam delivery in order to enable the most efficient possible use of beam current and to be able to respond flexibly to different irradiation field requirements. Therefore, a model was developed which calculates beam profiles taking into account various parameters of the measurement geometry, beam characteristics and field generation. This model was tested and verified in comparison with Monte Carlo simulations and measurements. Subsequently investigations into the positioning of the components used could be carried out and their influences on penumbra, homogeneity and transmission of the beam fields analysed. It was found that the high clinical requirements in terms of maximum penumbra and field homogeneity for the application in eye tumour therapy at HZB can be met if stronger magnetic coils were used. Model calculations were used to calculate and characterize different beam profiles for an optimized setup. It could be shown that under these optimized conditions an applicability of the active system in eye tumour therapy is possible. Furthermore, a static modulator for the irradiation site at HZB was developed in the course of this work, manufactured using a 3D printer and successfully tested.

Motivation: Radiotherapy is an important modality in cancer treatment commonly using photon beams from compact electron linear accelerators. However, due to the inverse depth dose profile (Bragg peak) with maximum dose deposition at the end of their path, proton beams allow a dose escalation within the target volume and reduction in surrounding normal tissue. Up to 20% of all radiotherapy patients could benefit from proton therapy (PT). Conventional accelerators are utilized to obtain proton beams with therapeutic energies of 70 – 250 MeV. These beams are then transported to the patient via magnetic transferlines and a rotatable beamline, called gantry, which are large and bulky. PT requires huge capex, limiting it to only a few big centres worldwide treating much less than 1% of radiotherapy patients. The new particle acceleration by ultra-intense laser pulses occurs on micrometer scales, potentially enabling more compact PT facilities and increasing their widespread. These laser-accelerated proton (LAP) bunches have been observed recently with energies of up to 90 MeV and scaling models predict LAP with therapeutic energies with the next generation petawatt laser systems. Challenges: Intense pulses with maximum 10 Hz repetition rate, broad energy spectrum, large divergence and short duration characterize LAP beams. In contrast, conventional accelerators generate mono-energetic, narrow, quasi-continuous beams. A new multifunctional gantry is needed for LAP beams with a capture and collimation system to control initial divergence, an energy selection system (ESS) to filter variable energy widths and a large acceptance beam shaping and scanning system. An advanced magnetic technology is also required for a compact and light gantry design. Furthermore, new dose deposition models and treatment planning systems (TPS) are needed for high quality, efficient dose delivery. Materials and Methods: In conventional dose modelling, mono-energetic beams with decreasing energies are superimposed to deliver uniform spread-out Bragg peak (SOBP). The low repetition rate of LAP pulses puts a critical constraint on treatment time and it is highly inefficient to utilize conventional dose models. It is imperative to utilize unique LAP beam properties to reduce total treatment times. A new 1D Broad Energy Assorted depth dose Deposition (BEAD) model was developed. It could deliver similar SOBP by superimposing several LAP pulses with variable broad energy widths. The BEAD model sets the primary criteria for the gantry, i.e. to filter and transport pulses with up to 20 times larger energy widths than conventional beams for efficient dose delivery. Air-core pulsed magnets can reach up to 6 times higher peak magnetic fields than conventional iron-core magnets and the pulsed nature of laser-driven sources allowed their use to reduce the size and weight of the gantry. An isocentric gantry was designed with integrated laser-target assembly, beam capture and collimation, variable ESS and large acceptance achromatic beam transport. An advanced clinical gantry was designed later with a novel active beam shaping and scanning system, called ELPIS. The filtered beam outputs via the advanced gantry simulations were implemented in an advanced 3D TPS, called LAPCERR. A LAP beam gantry and TPS were brought together for the first time, and clinical feasibility was studied for the advanced gantry via tumour conformal dose calculations on real patient data. Furthermore, for realization of pulsed gantry systems, a first pulsed beamline section consisting of prototypes of a capturing solenoid and a sector magnet was designed and tested at tandem accelerator with 10MeV pulsed proton beams. A first air-core pulsed quadrupole was also designed. Results: An advanced gantry with the new ELPIS system was designed and simulated. Simulated results show that achromatic beams with actively selectable beam sizes in the range of 1 – 20 cm diameter with selectable energy widths ranging from 19 – 3% can be delivered via the advanced gantry. ELPIS can also scan these large beams to a 20 × 10 cm2 irradiation field. This gantry is about 2.5 m in height and about 3.5 m in length, which is about 4 times smaller in volume than the conventional PT gantries. The clinical feasibility study on a head and neck tumour patient shows that these filtered beams can deliver state-of-the-art 3D intensity modulated treatment plans. Experimental characterization of a prototype pulsed beamline section was performed successfully and the synchronization of proton pulse with peak magnetic field in the individual magnets was established. This showed the practical applicability and feasibility of pulsed beamlines. The newly designed pulsed quadrupole with three times higher field gradients than iron-core quadrupoles is already manufactured and will be tested in near future. Conclusion: The main hurdle towards laser-driven PT is a laser accelerator providing beams of therapeutic quality, i.e. energy, intensity, stability, reliability. Nevertheless, the presented advanced clinical gantry design presents a complete beam transport solution for future laser-driven sources and shows the prospect and limitations of a compact laser-driven PT facility. Further development in the LAP-CERR is needed as it has the potential to utilize advanced beam controls from the ELPIS system and optimize doses on the basis of advanced dose schemes, like partial volume irradiation, to bring treatment times further down. To realize the gantry concept, further research, development and testing in higher field and higher (up to 10 Hz) repetition rate pulsed magnets to cater therapeutic proton beams is crucial.

Die Protonentherapie, welche zur Behandlung von Tumoren mittlerweile weltweit eingesetzt wird, könnte von einem System zur Reichweitekontrolle des Protonenstrahls bedeutend profitieren. Einen Rückschluss auf die Position der Protonen in Echtzeit ermöglicht die prompte Gammastrahlung. Diese Sekundärstrahlung entsteht durch nukleare Wechselwirkungen der Protonen mit den Atomkernen im Gewebe. Derzeit existieren verschiedene Methoden die aus einer Messung der zeitlichen oder räumlichen Verteilung dieser hochenergetischen Photonen versuchen eine Reichweiteinformation zu gewinnen. Die Methoden zur Nutzung prompter Gammastrahlung beruhen auf den verlässlichen Voraussagen von Teilchentransportrechnungen. Im Vergleich zu Messungen zeigen diese jedoch Diskrepanzen in den Photonenproduktionsquerschnitten aufgrund mangelnder experimenteller Stützstellen. Um einen Beitrag zu der Datenlage zu leisten, soll das Emissionsspektrum prompter Gammastrahlung am Kohlenstoff bestimmt werden. Dazu werden mithilfe von Entfaltungsalgorithmen aus den aufbereiteten Messdaten und der simulierten Detektorantwort vorhandene Energielinien und deren Intensität extrahiert. Über eine Normierung auf die Anzahl einfallender Protonen erfolgt die Bestimmung von Ausbeuten prompter Gammastrahlung bei der Bestrahlung homogener Kohlenstofftargets mit Protonen.:1 Einleitung und Motivation 2 Grundlagen der Entfaltung von Gammaspektren 2.1 Inverses Problem 2.1.1 Gold-Dekonvolution 2.1.2 Spektrum-Stripping 3 Experimentelle Bestimmung des prompten Gammaspektrums 3.1 Vorbetrachtungen und Aufbau des Experimentes 3.2 Verarbeitung der Messdaten 3.3 Optimierung des Messaufbaus und Maßnahmen der Untergrundreduktion 3.3.1 Abschirmung der Detektoren und Untergrundkorrektur 3.3.2 Flugzeitdiskriminierung 4 Entfaltung der gewonnenen Gammaspektren 4.1 Bestimmung der Detektorantwort mit GEANT 4 4.1.1 Bestimmung der detektorspezifischen Energieauflösung und des Ansprechvermögens über Quellmessungen 4.2 Entfaltung 4.2.1 Ergebnisse der Gold-Dekonvolution 4.2.2 Ergebnisse des Spektrum-Stripping 4.3 Protonenfluenznormierung 5 Zusammenfassung und Ausblick Anhang A Depositionsspektren der Detektoren 2 (125°) und 3 (55°) 69 Literaturverzeichnis Abbildungsverzeichnis Tabellenverzeichnis Abkürzungsverzeichnis Danksagung Selbstständigkeitserklärung / Proton therapy, which has become a clinically established technology for cancer treatment worldwide, would greatly benefit from an appropriate range verification system. A signature, which can be used to trace the range of the primary protons in real time during the treatment, is the so-called prompt gamma radiation. It is produced in nuclear reactions of the protons with the nuclei in tissue. Currently, there are different approaches which try to decode a range information from the spatial or time distribution of these high energetic photons. Methods utilizing prompt gamma-rays rely on the prediction of particle transport simulations. In comparison to measurements these simulations show severe deviations in prompt gamma-ray yield due to insufficient availability of experimental data. To make a contribution to the data situation, the prompt gamma-ray emission spectra for carbon should be investigated. Therefore deconvolution algorithms are applied to measured spectra with knowledge of the detector response function. Unfolded energy lines and their intensities are examined. With help of proton fluence some yields for the irradiation of homogeneous graphite targets are determined.:1 Einleitung und Motivation 2 Grundlagen der Entfaltung von Gammaspektren 2.1 Inverses Problem 2.1.1 Gold-Dekonvolution 2.1.2 Spektrum-Stripping 3 Experimentelle Bestimmung des prompten Gammaspektrums 3.1 Vorbetrachtungen und Aufbau des Experimentes 3.2 Verarbeitung der Messdaten 3.3 Optimierung des Messaufbaus und Maßnahmen der Untergrundreduktion 3.3.1 Abschirmung der Detektoren und Untergrundkorrektur 3.3.2 Flugzeitdiskriminierung 4 Entfaltung der gewonnenen Gammaspektren 4.1 Bestimmung der Detektorantwort mit GEANT 4 4.1.1 Bestimmung der detektorspezifischen Energieauflösung und des Ansprechvermögens über Quellmessungen 4.2 Entfaltung 4.2.1 Ergebnisse der Gold-Dekonvolution 4.2.2 Ergebnisse des Spektrum-Stripping 4.3 Protonenfluenznormierung 5 Zusammenfassung und Ausblick Anhang A Depositionsspektren der Detektoren 2 (125°) und 3 (55°) 69 Literaturverzeichnis Abbildungsverzeichnis Tabellenverzeichnis Abkürzungsverzeichnis Danksagung Selbstständigkeitserklärung

Es wird erwartet, dass die Integration der Magnetresonanztomografie (MRT) in die Protonentherapie die Treffgenauigkeit bei der Strahlentherapie für Krebserkrankungen deutlich verbessern wird. Besonders für Tumoren in beweglichen Organen des Thorax oder des Abdomens könnte die MRT-integrierte Protonentherapie (MRiPT) eine Synchronisierung der Bestrahlung mit der Tumorposition ermöglichen, was zu einer verminderten Normalgewebsdosis und weniger Nebenwirkungen führen könnte. Bis heute ist solch eine Integration jedoch aufgrund fehlender Studien zu potenziellen gegenseitigen Störeinflüssen dieser beiden Systeme nicht vollzogen worden. Diese Arbeit widmete sich zwei solcher Störeinflüsse, und zwar der Ablenkung des Protonenstrahls im Magnetfeld des MRT- Scanners, und umgekehrt, dem Einfluss der elekromagnetischen Felder der Protonentherapieanlage und des Protonenstrahls selbst auf die MRT-Bilder. Obwohl vorangegangene Studien den derzeitigen Konsens aufgezeigt haben, dass die Trajektorie eines abgebremsten Protonenstrahls im homogenen Phantom in einem transversalen Magnetfeld vorhersagbar ist, zeigte sich im quantitativen Vergleich der publizierten Modelle, der im ersten Teil dieser Arbeit vorgestellt wurde, dass die Vorhersagen dieser Modelle nur für eine begrenzte Anzahl von Kombinationen aus Magnetfeldstärke und Protonenenergie übereinstimmen. Die Schwächen bestehender analytischer Modelle wurden deshalb analysiert und quantifiziert. Kritische Annahmen und die mangelnde Anwendbarkeit auf realistische, d.h. inhomogene Magnetfeldstärken und Patientengeometrien wurden als Hauptprobleme identifiziert. Um diese zu überwinden, wurde ein neues semianalytisches Modell namens RAMDIM entwickelt. Es wurde gezeigt, dass dieses auf realistischere Fälle anwendbar und genauer ist als existierende analytische Modelle und dabei schneller als Monte-Carlo-basierte Teilchenspursimulationen. Es wird erwartet, dass dieses Modell in der MRiPT Anwendung findet zur schnellen und genauen Ablenkungsberechnung, zur Betrahlungsplanoptimierung und bei der MRT-geführten Strahlnachführung. In einem zweiten Schritt wurde die magnetfeldinduzierte Protonenstrahlablenkung in einem gewebeähnlichen Material durch Filmdosimetrie erstmalig gemessen und mit Monte-Carlo-Simulationen verglichen. In einem transversalen Magnetfeld einer Flussdichte von 0,95 T wurde experimentell gezeigt, dass die laterale Versetzung des Bragg-Peaks für Protonenenergien zwischen 80 und 180 MeV in PMMA zwischen 1 und 10 mm liegt. Die Retraktion des Bragg-Peaks war ≤ 0,5 mm. Es wurde gezeigt, dass die gemessene Versetzung des Bragg-Peaks innerhalb von 0,8 mm mit Monte-Carlo-basierten Vorhersagen übereinstimmt. Diese Ergebnisse weisen darauf hin, dass die Protonenstrahlablenkung durch Monte-Carlo-Simulationen genau vorhersagbar ist und damit der Realisierbarkeit der MRiPT nicht im Wege steht. Im zweiten Teil dieser Arbeit wurde erstmalig ein MRT-Scanner in eine Protonenstrahlführung integriert. Hierfür wurde ein offener Niederfeld-MRT-Scanner am Ende einer statischen Forschungsstrahlführung einer Protonentherapieanlage platziert. Die durch das statische Magnetfeld des MRT-Scanners hervorgerufene Strahlablenkung wurde bei der Ausrichtung des MRT-Scanners berücksichtigt. Die sequenzabhängigen, veränderlichen Gradientenfelder hatten keinen messbaren Einfluss auf das transversale Strahlprofil hinter dem MRT-Scanner. Die Magnetfeldhomogenität des Scanners lag innerhalb der Herstellervorgaben und zeigte keinen relevanten Einfluss von Rotationen der Protonengantry im benachbarten Bestrahlungsraum. Eine magnetische Abschirmung war zum gleichzeitigen Betrieb des MRT-Scanners und der Protonentherapieanlage nicht notwendig. Dies beweist die Machbarkeit gleichzeitiger Bestrahlung und Bildgebung in einem ersten MRiPT Aufbau. Die MRT-Bildqualität des Aufbaus wurde darauffolgend anhand eines angepassten Standardprotokolls aus Spin-Echo- und Gradienten-Echo-Sequenzen quantifiziert und es wurde gezeigt, dass die Bildqualität sowohl ohne als auch mit gleichzeitiger Bestrahlung hinreichend ist. Alle bestimmten geometrischen Parameter stimmten mit den physikalischen Abmessungen des verwendeten Phantoms innerhalb eines Bildpixels überein. Wie es für Niederfeld-MRT-Scanner üblich ist, war das Signal-Rausch-Verhältnis (SNR) der MRT-Bilder gering, was im Vergleich zu den Standardkriterien zu einer geringen Bildhomogenität und zu einem hohen Geisterbildanteil im Bild führte. Außerdem wurde aufgrund von Unsicherheiten in der Hochfrequenzkalibrierung des MRT-Scanners eine starke Schwankung der vertikalen Phantomposition mit einem Interquartilabstand von bis zu 1,5 mm beobachtet. T2*-gewichtete Gradientenechosequenzen zeigten zudem aufgrund von Magnetfeldinho- mogenitäten relevante ortsabhängige Bildverzerrungen. Es wurde gezeigt, dass die meisten Bildqualitätsparameter mit und ohne gleichzeitige Betrahlung äquivalent sind. Es wurde jedoch ein signifikanter Betrahlungseinfluss in Form von einer vertikalen Bildverschiebung und einer Verminderung des SNR beobachtet, die durch eine Änderung im Magnetfeld des MRT-Scanners erklärt werden können, welche durch zu diesem Feld parallel ausgerichtete Komponenten im Fernfeld der Strahlführungsmagneten hervorgerufen wird. Während das verminderte SNR vermutlich irrelevant ist (Dif- ferenz im Median ≤ 1,5), ist die sequenzabhängige Bildverschiebung (Differenz im Median bis zu 0,7 mm) nicht immer vernachlässigbar. Diese Ergebisse zeigen, dass die MRT-Bilder durch gleichzeitige Bildgebung nicht schwerwiegend verfälscht werden, dass aber eine dedizierte Optimierung der Hochfrequenzkalibrierung und der MRT-Bildsequenzen notwendig ist. Im letzten Teil der Arbeit wurde gezeigt, dass ein stromabhängiger Einfluss des Protonenstrahls auf MRT-Bilder eines Wasserphantoms durch zwei verschiedene MRT-Sequenzen messbar gemacht und zur Reichweiteverifikation genutzt werden kann. Der Effekt war in verschiedenen Flüssigkeiten, jedoch nicht in viskosen und festen Materialen, nachweisbar und wurde auf Hitzekonvektion zurückgeführt. Es wird erwartet, dass diese Methode in der MRiPT für Konstanztests der Protonenreichweite bei der Maschinenqualitätssicherung nützlich sein wird. Zusammenfassend hat diese Arbeit die Genauigkeit der Vorhersage der Strahlablenkung quantifiziert und verbessert, sowie Potenzial und Realisierbarkeit einer gleichzeitigen MRT-Bildgebung und Protonenbestrahlung gezeigt. Die weitere Entwicklung eines ersten MRiPT-Prototyps ist demnach gerechtfertigt.:List of Figures v List of Tables vii 1 General Introduction 1 2 State of the Art: Proton Therapy and Magnetic Resonance Imaging 3 2.1 Proton Therapy 4 2.1.1 Physical Principle 4 2.1.2 Beam Delivery 7 2.1.3 Motion Management and the Role of Image Guidance 10 2.2 Magnetic Resonance Imaging 14 2.2.1 Physical Principle 14 2.2.2 Image Generation by Pulse Sequences 18 2.2.3 Image Quality 21 2.3 MR-Guided Radiotherapy 24 2.3.1 Offline MR Guidance 24 2.3.2 On-line MR Guidance 25 2.4 MR-Integrated Proton Therapy 28 2.4.1 Aims of this Thesis 32 3 Magnetic Field-Induced Beam Deflection and Bragg Peak Displacement 35 3.1 Analytical Description 36 3.1.1 Review of Analytical Models 36 3.1.2 New Model Formulation 41 3.1.3 Evaluation of Analytical and Numerical Models 44 3.1.4 Discussion 51 3.2 Monte Carlo Simulation and Experimental Verification 54 3.2.1 Verification Setup 54 3.2.2 Monte Carlo Simulation 56 3.2.3 Experimental Verification 60 3.2.4 Discussion 61 3.3 Summary 63 4 Integrated In-Beam MR System: Proof of Concept 65 4.1 Integration of a Low-Field MR Scanner and a Static Research Beamline 65 4.1.1 Proton Therapy System 66 4.1.2 MR Scanner 66 4.1.3 Potential Sources of Interference 67 4.1.4 Integration of Both Systems 68 4.2 Beam and Image Quality in the Integrated Setup 70 4.2.1 Beam Profile 70 4.2.2 MR Magnetic Field Homogeneity 72 4.2.3 MR Image Quality - Qualitative In Vivo and Ex Vivo Test 74 4.2.4 MR Image Quality - Quantitative Phantom Tests 77 4.3 Feasibility of MRI-based Range Verification 86 4.3.1 MR Sequences 86 4.3.2 Proton Beam Parameters 88 4.3.3 Target Material Dependence 91 4.3.4 Discussion 92 4.4 Summary 96 5 Discussion and Future Perspectives 99 6 Summary/Zusammenfassung 105 6.1 Summary 105 6.2 Zusammenfassung 108 Bibliography I Supplementary Information XXIX A Beam Deflection: Experimental Measurements XXIX A.1 Setup XXIX A.2 Film Handling and Evaluation XXX A.3 Uncertainty Estimation XXX B Beam Deflection: Monte Carlo Simulations XXXIII B.1 Magnetic Field Model XXXIII B.2 Uncertainty Estimation XXXIV C Integrated MRiPT Setup XXXVI C.1 Magnetic Field Map XXXVI C.2 Sequence Parameters XXXVI C.3 Image Quality Parameters XLII C.4 Range Verification Sequences XLII / The integration of magnetic resonance imaging (MRI) into proton therapy is expected to strongly increase the targeting accuracy in radiation therapy for cancerous diseases. Especially for tumours situated in mobile organs in the thorax and abdomen, MR-integrated proton therapy (MRiPT) could enable the synchronisation of irradiation to the tumour position, resulting in less dose to normal tissue and reduced side effects. However, such an integration has been hindered so far by a lack of scientific studies on the potential mutual interference between the two components. This thesis was dedicated to two of these sources of interference, namely the deflection of the proton beam by the magnetic field of the MR scanner and, vice versa, alterations of the MR image induced by the electromagnetic fields of the proton therapy facility and by the beam itself. Although previous work has indicated that there is general consensus that the trajectory of a slowing down proton beam in a homogeneous phantom inside a transverse magnetic field is predictable, a quantitative comparison of the published methods, as presented in the first part of this thesis, has shown that predictions of different models only agree for certain proton beam energies and magnetic flux densities. Therefore, shortcomings of previously published analytical methods have been analysed and quantified. The inclusion of critical assumptions and the lack of applicability to realistic, i.e. non-uniform, magnetic flux densities and patient anatomies have been identified as main problems. To overcome these deficiencies, a new semi-analytical model called RAMDIM has been developed. It was shown that this model is both applicable to more realistic setups and less assumptive than existing analytical approaches, and faster than Monte Carlo based particle tracking simulations. This model is expected to be useful in MRiPT for fast and accurate deflection estimations, treatment plan optimisation, and MR-guided beam tracking. In a second step, the magnetic field-induced proton beam deflection has been measured for the first time in a tissue-mimicking medium by film dosimetry and has been compared against Monte Carlo simulations. In a transverse magnetic field of 0.95 T, it was experimentally shown that the lateral Bragg peak displacement ranges between 1 mm and 10 mm for proton energies between 80 and 180 MeV in PMMA. Range retraction was found to be ≤ 0.5 mm. The measured Bragg peak displacement was shown to agree within 0.8 mm with Monte Carlo simulations. These results indicate that proton beam deflection in a homogeneous medium is accurately predictable for intermediate proton beam energies and magnetic flux densities by Monte Carlo simulations and therefore not impeding the feasibility of MRiPT. In the second part of this thesis, an MR scanner has been integrated into a proton beam line for the first time. For this purpose, an open low-field MR scanner has been placed at the end of a fixed horizontal proton research beam line in a proton therapy facility. The beam deflection induced by the static magnetic field of the scanner was taken into account for alignment of the beam and the FOV of the scanner. The pulse sequence-dependent dynamic gradient fields did not measurably affect the transverse beam profile behind the MR scanner. The MR magnetic field homogeneity was within the vendor’s specifications and not relevantly influenced by the rotation of the proton gantry in the neighbouring treatment room. No magnetic field compensation system was required for simultaneous operation of the MR scanner and the proton therapy system. These results proof that simultaneous irradiation and imaging is feasible in an in-beam MR setup. The MR image quality of the in-beam MR scanner was then quantified by an adapted standard protocol comprising spin and gradient echo imaging and shown to be acceptable both with and without simultaneous proton beam irradiation. All geometrical parameters agreed with the mechanical dimensions of the used phantom within one pixel width. As common for low-field MR scanners, the signal-to-noise ratio (SNR) of the MR images was low, which resulted in a low image uniformity and a high ghosting ratio in comparison to the standardised test criteria. Furthermore, a strong fluctuation of the vertical phantom position due to uncertainties in the pre-scan frequency calibration was observed, with an interquartile range of up to 1.5 mm. T2*-weighted gradient echo images showed relevant nonuniform deformations due to magnetic field inhomogeneities. Most image quality parameters were shown to be equivalent with and without simultaneous proton beam irradiation. However, a significant influence of simultaneous irradiation was observed as a shift of the vertical phantom position and a decrease in the SNR, both of which can be explained by a change in the B0 field of the MR scanner induced by components of the fringe field of the beam line magnets directed parallel to B0 . While the decrease in SNR is not expected to be relevant (median differences were within 1.5 ), the sequence-dependent phantom shift (median differences of up to 0.7 mm) can become non-negligible. These results show that the MR images are not severely distorted by simultaneous irradiation, but a dedicated optimisation of the pre-scan RF calibration and the MR sequences is required for MRiPT. Lastly, a current-dependent influence of the proton beam on the MR image was shown to be measurable in water in two different MR sequences, which allowed for range verification measurements. The effect was observed in different liquids but not in highly viscose and solid materials, and most probably induced by heat convection. This method is expected to be useful in MRiPT for consistency tests of the proton range during machine-specific quality assurance. In conclusion, this work has improved and quantified the accuracy of beam deflection predictions and shown the feasibility and potential of in-beam MR imaging, justifying further research towards a first MRiPT prototype.:List of Figures v List of Tables vii 1 General Introduction 1 2 State of the Art: Proton Therapy and Magnetic Resonance Imaging 3 2.1 Proton Therapy 4 2.1.1 Physical Principle 4 2.1.2 Beam Delivery 7 2.1.3 Motion Management and the Role of Image Guidance 10 2.2 Magnetic Resonance Imaging 14 2.2.1 Physical Principle 14 2.2.2 Image Generation by Pulse Sequences 18 2.2.3 Image Quality 21 2.3 MR-Guided Radiotherapy 24 2.3.1 Offline MR Guidance 24 2.3.2 On-line MR Guidance 25 2.4 MR-Integrated Proton Therapy 28 2.4.1 Aims of this Thesis 32 3 Magnetic Field-Induced Beam Deflection and Bragg Peak Displacement 35 3.1 Analytical Description 36 3.1.1 Review of Analytical Models 36 3.1.2 New Model Formulation 41 3.1.3 Evaluation of Analytical and Numerical Models 44 3.1.4 Discussion 51 3.2 Monte Carlo Simulation and Experimental Verification 54 3.2.1 Verification Setup 54 3.2.2 Monte Carlo Simulation 56 3.2.3 Experimental Verification 60 3.2.4 Discussion 61 3.3 Summary 63 4 Integrated In-Beam MR System: Proof of Concept 65 4.1 Integration of a Low-Field MR Scanner and a Static Research Beamline 65 4.1.1 Proton Therapy System 66 4.1.2 MR Scanner 66 4.1.3 Potential Sources of Interference 67 4.1.4 Integration of Both Systems 68 4.2 Beam and Image Quality in the Integrated Setup 70 4.2.1 Beam Profile 70 4.2.2 MR Magnetic Field Homogeneity 72 4.2.3 MR Image Quality - Qualitative In Vivo and Ex Vivo Test 74 4.2.4 MR Image Quality - Quantitative Phantom Tests 77 4.3 Feasibility of MRI-based Range Verification 86 4.3.1 MR Sequences 86 4.3.2 Proton Beam Parameters 88 4.3.3 Target Material Dependence 91 4.3.4 Discussion 92 4.4 Summary 96 5 Discussion and Future Perspectives 99 6 Summary/Zusammenfassung 105 6.1 Summary 105 6.2 Zusammenfassung 108 Bibliography I Supplementary Information XXIX A Beam Deflection: Experimental Measurements XXIX A.1 Setup XXIX A.2 Film Handling and Evaluation XXX A.3 Uncertainty Estimation XXX B Beam Deflection: Monte Carlo Simulations XXXIII B.1 Magnetic Field Model XXXIII B.2 Uncertainty Estimation XXXIV C Integrated MRiPT Setup XXXVI C.1 Magnetic Field Map XXXVI C.2 Sequence Parameters XXXVI C.3 Image Quality Parameters XLII C.4 Range Verification Sequences XLII

Derzeit gibt es eine intensive Debatte über die Notwendigkeit, Variationen der relativen biologischen Wirksamkeit (RBE) in der Protonentherapie zu berücksichtigen. Hier wurde die Variabilität der RBE für späte strahleninduzierte Hirnverletzungen (RIBI) untersucht, die nach einer Protonentherapie von Gliom-Hirntumorpatienten diagnostiziert wurden. Eine Gliomkohorte wurde definiert und auf späte RIBI untersucht, die in der Patienten-Nachsorge beobachtet wurden. Die RIBI Läsionen wurden mit linearen Energietransfer- und variablen RBE-Verteilungen korreliert, die unter Verwendung eines etablierten Monte-Carlo Frameworks berechnet wurden. Ein klinisches RBE Modell wurde erstellt und zur Evaluierung neuartiger Behandlungsstrategien angewendet. Die periventrikuläre Region wurde als anatomischer Risikofaktor für RIBI identifiziert und als neues Risikoorgan in die Protonenbehandlungsplanung an der Universitäts Protonen Therapie Dresden aufgenommen. Die Arbeit demonstriert die klinische Relevanz und präsentiert erste translatorische Schritte hin zur biologisch optimierten Protonentherapie unter Verwendung einer variablen RBE. / Currently, there is an intense debate on the need to consider variations in relative biological effectiveness (RBE) in proton therapy. Here, the variability of RBE was investigated for late radiation-induced brain injuries (RIBI) observed after proton therapy in glioma brain tumor patients. A glioma cohort was defined and evaluated for late RIBI observed in patient follow-up. Injury lesions were correlated to linear energy transfer and variable RBE distributions simulated using an established Monte-Carlo framework. A clinical RBE model was established and applied in treatment response assessment of novel treatment strategies. The periventricular region was identified as an anatomical risk factor for RIBI and incorporated as a new organ at risk in proton treatment planning at University Proton Therapy Dresden. The thesis demonstrates the clinical relevance and presents first translational steps towards biologically optimized proton therapy using a variable RBE.

Es wird erwartet, dass die Integration der Magnetresonanztomografie (MRT) in die Protonentherapie die Treffgenauigkeit bei der Strahlentherapie für Krebserkrankungen deutlich verbessern wird. Besonders für Tumoren in beweglichen Organen des Thorax oder des Abdomens könnte die MRT-integrierte Protonentherapie (MRiPT) eine Synchronisierung der Bestrahlung mit der Tumorposition ermöglichen, was zu einer verminderten Normalgewebsdosis und weniger Nebenwirkungen führen könnte. Bis heute ist solch eine Integration jedoch aufgrund fehlender Studien zu potenziellen gegenseitigen Störeinflüssen dieser beiden Systeme nicht vollzogen worden. Diese Arbeit widmete sich zwei solcher Störeinflüsse, und zwar der Ablenkung des Protonenstrahls im Magnetfeld des MRT- Scanners, und umgekehrt, dem Einfluss der elekromagnetischen Felder der Protonentherapieanlage und des Protonenstrahls selbst auf die MRT-Bilder. Obwohl vorangegangene Studien den derzeitigen Konsens aufgezeigt haben, dass die Trajektorie eines abgebremsten Protonenstrahls im homogenen Phantom in einem transversalen Magnetfeld vorhersagbar ist, zeigte sich im quantitativen Vergleich der publizierten Modelle, der im ersten Teil dieser Arbeit vorgestellt wurde, dass die Vorhersagen dieser Modelle nur für eine begrenzte Anzahl von Kombinationen aus Magnetfeldstärke und Protonenenergie übereinstimmen. Die Schwächen bestehender analytischer Modelle wurden deshalb analysiert und quantifiziert. Kritische Annahmen und die mangelnde Anwendbarkeit auf realistische, d.h. inhomogene Magnetfeldstärken und Patientengeometrien wurden als Hauptprobleme identifiziert. Um diese zu überwinden, wurde ein neues semianalytisches Modell namens RAMDIM entwickelt. Es wurde gezeigt, dass dieses auf realistischere Fälle anwendbar und genauer ist als existierende analytische Modelle und dabei schneller als Monte-Carlo-basierte Teilchenspursimulationen. Es wird erwartet, dass dieses Modell in der MRiPT Anwendung findet zur schnellen und genauen Ablenkungsberechnung, zur Betrahlungsplanoptimierung und bei der MRT-geführten Strahlnachführung. In einem zweiten Schritt wurde die magnetfeldinduzierte Protonenstrahlablenkung in einem gewebeähnlichen Material durch Filmdosimetrie erstmalig gemessen und mit Monte-Carlo-Simulationen verglichen. In einem transversalen Magnetfeld einer Flussdichte von 0,95 T wurde experimentell gezeigt, dass die laterale Versetzung des Bragg-Peaks für Protonenenergien zwischen 80 und 180 MeV in PMMA zwischen 1 und 10 mm liegt. Die Retraktion des Bragg-Peaks war ≤ 0,5 mm. Es wurde gezeigt, dass die gemessene Versetzung des Bragg-Peaks innerhalb von 0,8 mm mit Monte-Carlo-basierten Vorhersagen übereinstimmt. Diese Ergebnisse weisen darauf hin, dass die Protonenstrahlablenkung durch Monte-Carlo-Simulationen genau vorhersagbar ist und damit der Realisierbarkeit der MRiPT nicht im Wege steht. Im zweiten Teil dieser Arbeit wurde erstmalig ein MRT-Scanner in eine Protonenstrahlführung integriert. Hierfür wurde ein offener Niederfeld-MRT-Scanner am Ende einer statischen Forschungsstrahlführung einer Protonentherapieanlage platziert. Die durch das statische Magnetfeld des MRT-Scanners hervorgerufene Strahlablenkung wurde bei der Ausrichtung des MRT-Scanners berücksichtigt. Die sequenzabhängigen, veränderlichen Gradientenfelder hatten keinen messbaren Einfluss auf das transversale Strahlprofil hinter dem MRT-Scanner. Die Magnetfeldhomogenität des Scanners lag innerhalb der Herstellervorgaben und zeigte keinen relevanten Einfluss von Rotationen der Protonengantry im benachbarten Bestrahlungsraum. Eine magnetische Abschirmung war zum gleichzeitigen Betrieb des MRT-Scanners und der Protonentherapieanlage nicht notwendig. Dies beweist die Machbarkeit gleichzeitiger Bestrahlung und Bildgebung in einem ersten MRiPT Aufbau. Die MRT-Bildqualität des Aufbaus wurde darauffolgend anhand eines angepassten Standardprotokolls aus Spin-Echo- und Gradienten-Echo-Sequenzen quantifiziert und es wurde gezeigt, dass die Bildqualität sowohl ohne als auch mit gleichzeitiger Bestrahlung hinreichend ist. Alle bestimmten geometrischen Parameter stimmten mit den physikalischen Abmessungen des verwendeten Phantoms innerhalb eines Bildpixels überein. Wie es für Niederfeld-MRT-Scanner üblich ist, war das Signal-Rausch-Verhältnis (SNR) der MRT-Bilder gering, was im Vergleich zu den Standardkriterien zu einer geringen Bildhomogenität und zu einem hohen Geisterbildanteil im Bild führte. Außerdem wurde aufgrund von Unsicherheiten in der Hochfrequenzkalibrierung des MRT-Scanners eine starke Schwankung der vertikalen Phantomposition mit einem Interquartilabstand von bis zu 1,5 mm beobachtet. T2*-gewichtete Gradientenechosequenzen zeigten zudem aufgrund von Magnetfeldinho- mogenitäten relevante ortsabhängige Bildverzerrungen. Es wurde gezeigt, dass die meisten Bildqualitätsparameter mit und ohne gleichzeitige Betrahlung äquivalent sind. Es wurde jedoch ein signifikanter Betrahlungseinfluss in Form von einer vertikalen Bildverschiebung und einer Verminderung des SNR beobachtet, die durch eine Änderung im Magnetfeld des MRT-Scanners erklärt werden können, welche durch zu diesem Feld parallel ausgerichtete Komponenten im Fernfeld der Strahlführungsmagneten hervorgerufen wird. Während das verminderte SNR vermutlich irrelevant ist (Dif- ferenz im Median ≤ 1,5), ist die sequenzabhängige Bildverschiebung (Differenz im Median bis zu 0,7 mm) nicht immer vernachlässigbar. Diese Ergebisse zeigen, dass die MRT-Bilder durch gleichzeitige Bildgebung nicht schwerwiegend verfälscht werden, dass aber eine dedizierte Optimierung der Hochfrequenzkalibrierung und der MRT-Bildsequenzen notwendig ist. Im letzten Teil der Arbeit wurde gezeigt, dass ein stromabhängiger Einfluss des Protonenstrahls auf MRT-Bilder eines Wasserphantoms durch zwei verschiedene MRT-Sequenzen messbar gemacht und zur Reichweiteverifikation genutzt werden kann. Der Effekt war in verschiedenen Flüssigkeiten, jedoch nicht in viskosen und festen Materialen, nachweisbar und wurde auf Hitzekonvektion zurückgeführt. Es wird erwartet, dass diese Methode in der MRiPT für Konstanztests der Protonenreichweite bei der Maschinenqualitätssicherung nützlich sein wird. Zusammenfassend hat diese Arbeit die Genauigkeit der Vorhersage der Strahlablenkung quantifiziert und verbessert, sowie Potenzial und Realisierbarkeit einer gleichzeitigen MRT-Bildgebung und Protonenbestrahlung gezeigt. Die weitere Entwicklung eines ersten MRiPT-Prototyps ist demnach gerechtfertigt.:List of Figures v List of Tables vii 1 General Introduction 1 2 State of the Art: Proton Therapy and Magnetic Resonance Imaging 3 2.1 Proton Therapy 4 2.1.1 Physical Principle 4 2.1.2 Beam Delivery 7 2.1.3 Motion Management and the Role of Image Guidance 10 2.2 Magnetic Resonance Imaging 14 2.2.1 Physical Principle 14 2.2.2 Image Generation by Pulse Sequences 18 2.2.3 Image Quality 21 2.3 MR-Guided Radiotherapy 24 2.3.1 Offline MR Guidance 24 2.3.2 On-line MR Guidance 25 2.4 MR-Integrated Proton Therapy 28 2.4.1 Aims of this Thesis 32 3 Magnetic Field-Induced Beam Deflection and Bragg Peak Displacement 35 3.1 Analytical Description 36 3.1.1 Review of Analytical Models 36 3.1.2 New Model Formulation 41 3.1.3 Evaluation of Analytical and Numerical Models 44 3.1.4 Discussion 51 3.2 Monte Carlo Simulation and Experimental Verification 54 3.2.1 Verification Setup 54 3.2.2 Monte Carlo Simulation 56 3.2.3 Experimental Verification 60 3.2.4 Discussion 61 3.3 Summary 63 4 Integrated In-Beam MR System: Proof of Concept 65 4.1 Integration of a Low-Field MR Scanner and a Static Research Beamline 65 4.1.1 Proton Therapy System 66 4.1.2 MR Scanner 66 4.1.3 Potential Sources of Interference 67 4.1.4 Integration of Both Systems 68 4.2 Beam and Image Quality in the Integrated Setup 70 4.2.1 Beam Profile 70 4.2.2 MR Magnetic Field Homogeneity 72 4.2.3 MR Image Quality - Qualitative In Vivo and Ex Vivo Test 74 4.2.4 MR Image Quality - Quantitative Phantom Tests 77 4.3 Feasibility of MRI-based Range Verification 86 4.3.1 MR Sequences 86 4.3.2 Proton Beam Parameters 88 4.3.3 Target Material Dependence 91 4.3.4 Discussion 92 4.4 Summary 96 5 Discussion and Future Perspectives 99 6 Summary/Zusammenfassung 105 6.1 Summary 105 6.2 Zusammenfassung 108 Bibliography I Supplementary Information XXIX A Beam Deflection: Experimental Measurements XXIX A.1 Setup XXIX A.2 Film Handling and Evaluation XXX A.3 Uncertainty Estimation XXX B Beam Deflection: Monte Carlo Simulations XXXIII B.1 Magnetic Field Model XXXIII B.2 Uncertainty Estimation XXXIV C Integrated MRiPT Setup XXXVI C.1 Magnetic Field Map XXXVI C.2 Sequence Parameters XXXVI C.3 Image Quality Parameters XLII C.4 Range Verification Sequences XLII / The integration of magnetic resonance imaging (MRI) into proton therapy is expected to strongly increase the targeting accuracy in radiation therapy for cancerous diseases. Especially for tumours situated in mobile organs in the thorax and abdomen, MR-integrated proton therapy (MRiPT) could enable the synchronisation of irradiation to the tumour position, resulting in less dose to normal tissue and reduced side effects. However, such an integration has been hindered so far by a lack of scientific studies on the potential mutual interference between the two components. This thesis was dedicated to two of these sources of interference, namely the deflection of the proton beam by the magnetic field of the MR scanner and, vice versa, alterations of the MR image induced by the electromagnetic fields of the proton therapy facility and by the beam itself. Although previous work has indicated that there is general consensus that the trajectory of a slowing down proton beam in a homogeneous phantom inside a transverse magnetic field is predictable, a quantitative comparison of the published methods, as presented in the first part of this thesis, has shown that predictions of different models only agree for certain proton beam energies and magnetic flux densities. Therefore, shortcomings of previously published analytical methods have been analysed and quantified. The inclusion of critical assumptions and the lack of applicability to realistic, i.e. non-uniform, magnetic flux densities and patient anatomies have been identified as main problems. To overcome these deficiencies, a new semi-analytical model called RAMDIM has been developed. It was shown that this model is both applicable to more realistic setups and less assumptive than existing analytical approaches, and faster than Monte Carlo based particle tracking simulations. This model is expected to be useful in MRiPT for fast and accurate deflection estimations, treatment plan optimisation, and MR-guided beam tracking. In a second step, the magnetic field-induced proton beam deflection has been measured for the first time in a tissue-mimicking medium by film dosimetry and has been compared against Monte Carlo simulations. In a transverse magnetic field of 0.95 T, it was experimentally shown that the lateral Bragg peak displacement ranges between 1 mm and 10 mm for proton energies between 80 and 180 MeV in PMMA. Range retraction was found to be ≤ 0.5 mm. The measured Bragg peak displacement was shown to agree within 0.8 mm with Monte Carlo simulations. These results indicate that proton beam deflection in a homogeneous medium is accurately predictable for intermediate proton beam energies and magnetic flux densities by Monte Carlo simulations and therefore not impeding the feasibility of MRiPT. In the second part of this thesis, an MR scanner has been integrated into a proton beam line for the first time. For this purpose, an open low-field MR scanner has been placed at the end of a fixed horizontal proton research beam line in a proton therapy facility. The beam deflection induced by the static magnetic field of the scanner was taken into account for alignment of the beam and the FOV of the scanner. The pulse sequence-dependent dynamic gradient fields did not measurably affect the transverse beam profile behind the MR scanner. The MR magnetic field homogeneity was within the vendor’s specifications and not relevantly influenced by the rotation of the proton gantry in the neighbouring treatment room. No magnetic field compensation system was required for simultaneous operation of the MR scanner and the proton therapy system. These results proof that simultaneous irradiation and imaging is feasible in an in-beam MR setup. The MR image quality of the in-beam MR scanner was then quantified by an adapted standard protocol comprising spin and gradient echo imaging and shown to be acceptable both with and without simultaneous proton beam irradiation. All geometrical parameters agreed with the mechanical dimensions of the used phantom within one pixel width. As common for low-field MR scanners, the signal-to-noise ratio (SNR) of the MR images was low, which resulted in a low image uniformity and a high ghosting ratio in comparison to the standardised test criteria. Furthermore, a strong fluctuation of the vertical phantom position due to uncertainties in the pre-scan frequency calibration was observed, with an interquartile range of up to 1.5 mm. T2*-weighted gradient echo images showed relevant nonuniform deformations due to magnetic field inhomogeneities. Most image quality parameters were shown to be equivalent with and without simultaneous proton beam irradiation. However, a significant influence of simultaneous irradiation was observed as a shift of the vertical phantom position and a decrease in the SNR, both of which can be explained by a change in the B0 field of the MR scanner induced by components of the fringe field of the beam line magnets directed parallel to B0 . While the decrease in SNR is not expected to be relevant (median differences were within 1.5 ), the sequence-dependent phantom shift (median differences of up to 0.7 mm) can become non-negligible. These results show that the MR images are not severely distorted by simultaneous irradiation, but a dedicated optimisation of the pre-scan RF calibration and the MR sequences is required for MRiPT. Lastly, a current-dependent influence of the proton beam on the MR image was shown to be measurable in water in two different MR sequences, which allowed for range verification measurements. The effect was observed in different liquids but not in highly viscose and solid materials, and most probably induced by heat convection. This method is expected to be useful in MRiPT for consistency tests of the proton range during machine-specific quality assurance. In conclusion, this work has improved and quantified the accuracy of beam deflection predictions and shown the feasibility and potential of in-beam MR imaging, justifying further research towards a first MRiPT prototype.:List of Figures v List of Tables vii 1 General Introduction 1 2 State of the Art: Proton Therapy and Magnetic Resonance Imaging 3 2.1 Proton Therapy 4 2.1.1 Physical Principle 4 2.1.2 Beam Delivery 7 2.1.3 Motion Management and the Role of Image Guidance 10 2.2 Magnetic Resonance Imaging 14 2.2.1 Physical Principle 14 2.2.2 Image Generation by Pulse Sequences 18 2.2.3 Image Quality 21 2.3 MR-Guided Radiotherapy 24 2.3.1 Offline MR Guidance 24 2.3.2 On-line MR Guidance 25 2.4 MR-Integrated Proton Therapy 28 2.4.1 Aims of this Thesis 32 3 Magnetic Field-Induced Beam Deflection and Bragg Peak Displacement 35 3.1 Analytical Description 36 3.1.1 Review of Analytical Models 36 3.1.2 New Model Formulation 41 3.1.3 Evaluation of Analytical and Numerical Models 44 3.1.4 Discussion 51 3.2 Monte Carlo Simulation and Experimental Verification 54 3.2.1 Verification Setup 54 3.2.2 Monte Carlo Simulation 56 3.2.3 Experimental Verification 60 3.2.4 Discussion 61 3.3 Summary 63 4 Integrated In-Beam MR System: Proof of Concept 65 4.1 Integration of a Low-Field MR Scanner and a Static Research Beamline 65 4.1.1 Proton Therapy System 66 4.1.2 MR Scanner 66 4.1.3 Potential Sources of Interference 67 4.1.4 Integration of Both Systems 68 4.2 Beam and Image Quality in the Integrated Setup 70 4.2.1 Beam Profile 70 4.2.2 MR Magnetic Field Homogeneity 72 4.2.3 MR Image Quality - Qualitative In Vivo and Ex Vivo Test 74 4.2.4 MR Image Quality - Quantitative Phantom Tests 77 4.3 Feasibility of MRI-based Range Verification 86 4.3.1 MR Sequences 86 4.3.2 Proton Beam Parameters 88 4.3.3 Target Material Dependence 91 4.3.4 Discussion 92 4.4 Summary 96 5 Discussion and Future Perspectives 99 6 Summary/Zusammenfassung 105 6.1 Summary 105 6.2 Zusammenfassung 108 Bibliography I Supplementary Information XXIX A Beam Deflection: Experimental Measurements XXIX A.1 Setup XXIX A.2 Film Handling and Evaluation XXX A.3 Uncertainty Estimation XXX B Beam Deflection: Monte Carlo Simulations XXXIII B.1 Magnetic Field Model XXXIII B.2 Uncertainty Estimation XXXIV C Integrated MRiPT Setup XXXVI C.1 Magnetic Field Map XXXVI C.2 Sequence Parameters XXXVI C.3 Image Quality Parameters XLII C.4 Range Verification Sequences XLII

Die Protonentherapie zeichnet sich durch eine konformale und fokussierte Tumorbestrahlung aus, die es ermöglicht, gesundes Gewebe besser zu schonen als bei der konventionellen Strahlentherapie. Dieses Potential wird jedoch durch Unsicherheiten bei der Vorhersage der Protonenreichweite im Gewebe oder durch anatomische Veränderungen über den Verlauf der Therapie eingeschränkt. In der vorliegenden Arbeit wurde daher der klinische Nutzen eines Reichweiteverifikationssystems auf Grundlage von Prompt-Gamma-Imaging (PGI) zur Behandlungskontrolle untersucht. Dafür wurden Messungen mit einem PGI-System während Prostata- und Kopf-Hals-Tumor-Bestrahlungen durchgeführt und retrospektiv ausgewertet. Einerseits konnte dabei mittels PGI die Genauigkeit verschiedener Methoden zur Reichweitevorhersage überprüft werden. Es zeigte sich, dass die 2019 klinisch eingeführte Methode zur Reichweitevorhersage (DirectSPR) nicht von der mit PGI gemessenen Protonenreichweite in Prostata-Tumor-Bestrahlungen abweicht, wodurch die Reduktion der auf DirectSPR basierenden Reichweiteunsicherheiten unabhängig bestätigt werden konnte. Andererseits konnte die Detektionsfähigkeit von PGI bei der Erkennung relevanter und nicht relevanter anatomischer Veränderungen in applizierten Bestrahlungsfeldern nachgewiesen werden. Insbesondere wurde für die feldweise Klassifizierung der Prostata-Bestrahlungen eine Sensitivität und Spezifität von 74% bzw. 79% festgestellt. Damit konnte in dieser Dissertation erstmals systematisch das klinische Anwendungspotential eines Systems zur PGI-Reichweiteverifikation gezeigt werden. Als zusätzliche Untersuchung wurde in einer Kollaboration mit dem Massachusetts General Hospital zum ersten Mal ein Vergleich zwischen zwei verschiedenen, auf prompter Gammastrahlung basierenden Systemen zur Reichweiteverifikation durchgeführt. Dazu wurde ein standardisiertes Studienprotokoll etabliert, welches die Vergleichbarkeit und die klinische Implementierung von Reichweiteverifikationssystemen generell unterstützen könnte.:1 Einleitung 2 Strahlentherapie mit Protonen 2.1 Physikalische Grundlagen der Protonentherapie 2.2 Behandlungsablauf in der Protonentherapie 2.2.1 Bildgebung zur Therapieplanung 2.2.2 Bestrahlungsplanung 2.2.3 Strahlapplikation 2.3 Genauigkeit in der Protonentherapie 2.3.1 Ursachen für Behandlungs- und Reichweiteunsicherheiten 2.3.2 Aktueller Stand der Behandlungs- und Reichweiteverifikation 3 Methodik der Reichweiteverifikation mittels Prompt-Gamma-Bildgebung (PGI) 3.1 Funktionsprinzip der PGI-Schlitzkamera 3.2 Datenaufnahme und -verarbeitung 3.2.1 Detektoraufbau und Signalaufnahme 3.2.2 PGI-Simulation und Bestimmung der Reichweiteabweichung 3.3 Charakterisierung des PGI-Prototyps 3.3.1 Kalibrierung des Systems 3.3.2 Positionierungspräzision 3.4 Überblick zur PRIMA-Studie 3.5 Experimentelle Studien zur PGI-Simulationsgenauigkeit 3.5.1 Abhängigkeit vom PGI-Sichtfeld und der Protonenenergie 3.5.2 Validierung der erweiterten Simulationssoftware 3.5.3 Abhängigkeit von der Tumorentität 3.5.4 Schlussfolgerungen 4 Validierung der CT-basierten Reichweitevorhersage mittels PGI 4.1 Konzept der Validierung 4.2 Gesamtabschätzung der Validierungsunsicherheit 4.3 Ergebnisse der Validierung 4.4 Diskussion 5 Detektionsfähigkeit anatomischer Veränderungen mittels PGI 5.1 Prinzipieller Aufbau der Studie 5.2 Grundwahrheit auf Basis von CT- und Dosisinformationen 5.2.1 Manuelle Klassifizierung 5.2.2 Klassifizierung auf Grundlage von integrierten Tiefendosisprofilen 5.2.3 Ergebnis der Etablierung einer CT-basierten Grundwahrheit 5.3 Etablierung einer Klassifikation auf Basis von PGI-Daten 5.3.1 Verarbeitung der PGI-Daten mittels Cluster-Algorithmus 5.3.2 Definition von spot- oder clusterbasierten Klassifikationsmodellen 5.4 Ergebnisse der PGI-Detektionsfähigkeit 5.4.1 Auswertung für Patienten mit Prostata-Tumor 5.4.2 Auswertung für Patienten mit Tumoren im Kopf-Hals-Bereich 5.5 Diskussion 6 Genauigkeit zweier Reichweiteverifikationsmethoden – bizentrischer Vergleich 6.1 Material und Methoden 6.1.1 Bildgebung 6.1.2 Bestrahlungsplanung 6.1.3 Durchführung und Auswertung 6.2 Ergebnisse 6.3 Diskussion 7 Zusammenfassung 8 Summary / Proton therapy is a conformal and focused irradiation of the tumor, which allows for a better sparing of healthy tissue than with conventional radiotherapy. However, this potential is limited by uncertainties from the proton range prediction in the patient or anatomical changes over the course of the treatment. Therefore, in this work, the clinical benefit of a range verification system based on the prompt-gamma-imaging (PGI) method for treatment verification was investigated. For this purpose, measurements were carried out with a PGI system during prostate and head and neck cancer irradiations and evaluated retrospectively. On the one hand, PGI was used to review the accuracy of several range prediction methods. The results showed that a specific method for range prediction (DirectSPR), which was clinically introduced in 2019, does not deviate from the PGI-measured proton range in prostate cancer irradiations. This means that the reduction of the range uncertainties with DirectSPR could be independently confirmed. On the other hand, the detection capability of PGI in identifying relevant and non-relevant anatomical changes in delivered treatment fields was demonstrated. In particular, for the fieldwise classification of prostate irradiations a sensitivity and specificity of 74% and 79% was determined, respectively. Thus, the clinical potential of a PGI range verification system was for the first time systematically demonstrated in this thesis. Furthermore, in a collaboration with the Massachusetts General Hospital a first-time comparison of two different range verification systems based on prompt gamma radiation was conducted. Therefore, a standardized study protocol was established, which could generally foster the comparability and clinical implementation of range verification systems.:1 Einleitung 2 Strahlentherapie mit Protonen 2.1 Physikalische Grundlagen der Protonentherapie 2.2 Behandlungsablauf in der Protonentherapie 2.2.1 Bildgebung zur Therapieplanung 2.2.2 Bestrahlungsplanung 2.2.3 Strahlapplikation 2.3 Genauigkeit in der Protonentherapie 2.3.1 Ursachen für Behandlungs- und Reichweiteunsicherheiten 2.3.2 Aktueller Stand der Behandlungs- und Reichweiteverifikation 3 Methodik der Reichweiteverifikation mittels Prompt-Gamma-Bildgebung (PGI) 3.1 Funktionsprinzip der PGI-Schlitzkamera 3.2 Datenaufnahme und -verarbeitung 3.2.1 Detektoraufbau und Signalaufnahme 3.2.2 PGI-Simulation und Bestimmung der Reichweiteabweichung 3.3 Charakterisierung des PGI-Prototyps 3.3.1 Kalibrierung des Systems 3.3.2 Positionierungspräzision 3.4 Überblick zur PRIMA-Studie 3.5 Experimentelle Studien zur PGI-Simulationsgenauigkeit 3.5.1 Abhängigkeit vom PGI-Sichtfeld und der Protonenenergie 3.5.2 Validierung der erweiterten Simulationssoftware 3.5.3 Abhängigkeit von der Tumorentität 3.5.4 Schlussfolgerungen 4 Validierung der CT-basierten Reichweitevorhersage mittels PGI 4.1 Konzept der Validierung 4.2 Gesamtabschätzung der Validierungsunsicherheit 4.3 Ergebnisse der Validierung 4.4 Diskussion 5 Detektionsfähigkeit anatomischer Veränderungen mittels PGI 5.1 Prinzipieller Aufbau der Studie 5.2 Grundwahrheit auf Basis von CT- und Dosisinformationen 5.2.1 Manuelle Klassifizierung 5.2.2 Klassifizierung auf Grundlage von integrierten Tiefendosisprofilen 5.2.3 Ergebnis der Etablierung einer CT-basierten Grundwahrheit 5.3 Etablierung einer Klassifikation auf Basis von PGI-Daten 5.3.1 Verarbeitung der PGI-Daten mittels Cluster-Algorithmus 5.3.2 Definition von spot- oder clusterbasierten Klassifikationsmodellen 5.4 Ergebnisse der PGI-Detektionsfähigkeit 5.4.1 Auswertung für Patienten mit Prostata-Tumor 5.4.2 Auswertung für Patienten mit Tumoren im Kopf-Hals-Bereich 5.5 Diskussion 6 Genauigkeit zweier Reichweiteverifikationsmethoden – bizentrischer Vergleich 6.1 Material und Methoden 6.1.1 Bildgebung 6.1.2 Bestrahlungsplanung 6.1.3 Durchführung und Auswertung 6.2 Ergebnisse 6.3 Diskussion 7 Zusammenfassung 8 Summary

Proton therapy is a relatively new treatment modality for cancer, having recently been incorporated into hospitals in the last two decades. Although proton therapy has much higher start up and treatment costs than traditional methods of radiotherapy, it continues to expand in use today. One reason for this is that proton therapy has the advantage of a more precise localization of dose compared to traditional radiotherapy. Other proposed advantages of proton therapy in the treatment of cancer may lead to a faster expanse in its use if proven to be more effective than traditional radiotherapy. Therefore, much research must be done to investigate the possible negative and positive effects of using proton therapy as a treatment modality. In proton therapy, protons do account for the vast majority of dose. However, when protons travel through matter, secondary particles are created by the interactions of protons and matter en route to and within the patient. It is believed that secondary dose can lead to secondary cancer, especially in pediatric cases. Therefore, the focus of this work is determining both primary and secondary dose. In order to develop relevant simulations, the specifications of the treatment room and beam were based off of real-world facilities as closely as possible. Using available data from proton accelerators and clinical facilities, an accurate proton therapy nozzle was designed. Dose calculations were performed by MCNPX using a simple water phantom, and then beam characteristics were investigated to ensure the accuracy of the model. After validation of the beam nozzle, primary and secondary dose values were tabulated and discussed. By demonstrating the method of these calculations, the purpose of this work is to serve as a guide into the relatively recent field of Monte Carlo methods in proton therapy.

Die Bestrahlungstherapie von Aderhautmelanomen mit der Protonentherapie stellt ein sehr präzises Behandlungsverfahren dar, welches hohe Ansprüche an die zugrundliegende Bildgebung stellt. Die Magnetresonanztomographie (MRT) und Computertomographie (CT) scheinen prinzipiell in der Lage, mit geeigneten Oberflächenspulen bzw. angepaßter Untersuchungstechnik diese Ansprüche zu erfüllen. Ziel der Arbeit ist die Bereitstellung und Überprüfung eines bildgebenden Systems, das einerseits die diagnostischen und differentialdiagnostischen Aufgaben für die Differenzierung und das Staging von Aderhautmelanomen erfüllt und andererseits als Grundlage der Bestrahlungsplanung dienen kann. Erstrebenswert ist dabei die Erstellung von Schnittbildern bestehend aus isotropen Voxeln, d.h. Voxeln mit identischer Kantenlänge, mit der Magnetresonanztomographie (MRT) und der Computertomographie (CT). Zur Differenzierung der unterschiedlichen Tumorentitäten wurden verschiedene Untersuchungen mit zum Teil großen Patientenkollektiven durchgeführt und histologisch und klinisch korreliert bzw. verifiziert. Dabei konnte das MR-morphologische Erscheinungsbild der relevanten Erkrankungen studiert werden. Die Genauigkeit der MR-gestützten Tumorvolumetrie wurde in vitro im Tierversuch und in vivo am Patientenauge überprüft. Zur Überprüfung der Genauigkeit der MR-Bildgebung mit isotropen Voxeln wurden unterschiedliche Spulen getestet und verschiedene Phantomuntersuchungen durchgeführt. Entsprechend wurde auch die CT evaluiert. Mit der digitalen Bildfusion wurden die MRT- und die CT-Daten in einen gemeinsamen Datensatz überführt, der sowohl zur Bestimmung der Genauigkeit beider Verfashren als auch als Basis zur Bestrahlungsplanung verwendet werden konnte. Mit einer optimierten Untersuchungstechnik konnten Schnittbilder auf der Basis von CT und MRT zur Verfügung gestellt werden, die eine präzise Planung der Protonenbestrahlung beim Aderhautmelanom ermöglichen. / The radiation therapy of uveal melanomas using proton therapy is a very precise therapy procedure, that makes high demands on the underlying imaging modalities. In general magnetic resonance imaging (MRI) and computed tomography (CT) have the potential to fulfill these criteria using suited surface coils and optimized imaging techniques, respectively. Purpose of the investigations was the preparation and evaluation of an imaging system, that allows diagnosis and differetial diagnosis of uveal melanomas including staging on the one hand. On the other hand it should be the base for high precision radiation therapy planning. The preparation of images with isotropic voxels, i.e. the edges of the voxels have identical size, is desirable with both magnetic resonance imaging (MRI) and computed tomography (CT). Several examinations on large patient groups were performed for the differentiation of different tumor entities and compared with clinical and histopathological outcome. The precision of the MR-based tumor volumetry was evaluated with animal studies in vitro and patient eyes in vivo. The MR-imaging with isotropic voxels was evaluated using several surface coils and phantom examinations. Comparable examinations were performed regarding CT. Digital image fusion of CT and MRI data sets was used to evaluate the precision of both modalities and for planning of radiation therapy. Using optimized examination techniques both CT and especially MRI are suitable for a precise proton therapy planning of uveal melanoma.