Person:

Paganetti, Harald

Background: To precisely calculate skin dose and thus to evaluate the relationship between the skin dose and permanent alopecia for pediatric medulloblastoma patients treated with proton beams. Methods: The dosimetry and alopecia outcomes of 12 children with medulloblastoma (ages 4-15 years) comprise the study cohort. Permanent alopecia was assessed and graded after completion of the entire therapy. Skin threshold doses of permanent alopecia were calculated based on the skin dose from the craniospinal irradiation (CSI) plan using the concept of generalized equivalent uniform dose (gEUD) and accounting for chemotherapy intensity. Monte Carlo simulations were employed to accurately assess uncertainties due to beam range prediction and secondary particles. Results: Increasing the dose of the CSI field or the dose given by the boost field to the posterior fossa increased total skin dose delivered in that region. It was found that permanent alopecia could be correlated with CSI dose with a threshold of about 21 Gy (relative biological effectiveness, RBE) with high dose chemotherapy and 30 Gy (RBE) with conventional chemotherapy. Conclusions: Our results based on 12 patients provide a relationship between the skin dose and permanent alopecia for pediatric medulloblastoma patients treated with protons. The alopecia risk as assessed with gEUD could be predicted based on the treatment plan information.

Purpose

To evaluate the lifespan and risk of cancer following whole-body exposure of mice to neutrons generated by a passively scattered clinical SOBP proton beam.

Methods and Materials

Three hundred young adult female FVB/N mice, 152 test and 148 control, were entered into the experiment. Mice were placed in an annular cassette around a cylindrical phantom, which was positioned lateral to the mid SOBP of a 165 MeV, clinical proton beam. The average distance from the edge of the mid SOBP to the conscious active mice was 21.5 cm. The phantom was irradiated with once daily fractions of 25 Gy, 4 days per week, for 6 weeks. The age at death and cause of death, i.e., cancer and type vs. non-cancer causes, were assessed over the lifespan of the mice.

Results

Exposure of mice to a dose of 600 Gy of proton beam generated neutrons, reduced the median lifespan of the mice by 4.2% (Kaplan-Meier cumulative survival, P = 0.053). The relative risk of death from cancer in neutron exposed vs. control mice was 1.40 for cancer of all types (P = 0.0006) and 1.22 for solid cancers (P = 0.09). For a typical 60 Gy dose of clinical protons, the observed 22% increased risk of solid cancer would be expected to decrease by a factor of 10.

Conclusions

Exposure of mice to neutrons generated by a proton dose which exceeds a typical course of radiotherapy by a factor of 10, resulted in a statistically significant increase in the background incidence of leukemia and a marginally significant increase in solid cancer. The results indicate that the risk of out-of-field 2nd solid cancers from SOBP proton generated neutrons and typical treatment schedules, is 6 - 10 times less than is suggested by current neutron risk estimates.

Monte Carlo models of proton therapy treatment heads are being used to improve beam delivery systems and to calculate the radiation field for patient dose calculations. The achievable accuracy of the model depends on the exact knowledge of the treatment head geometry and time structure, the material characteristics, and the underlying physics. This work aimed at studying the uncertainties in treatment head simulations for passive scattering proton therapy. The sensitivities of spread-out Bragg peak (SOBP) dose distributions on material densities, mean ionization potentials, initial proton beam energy spread and spot size were investigated. An improved understanding of the nature of these parameters may help to improve agreement between calculated and measured SOBP dose distributions and to ensure that the range, modulation width, and uniformity are within clinical tolerance levels. Furthermore, we present a method to make small corrections to the uniformity of spread-out Bragg peaks by utilizing the time structure of the beam delivery. In addition, we re-commissioned the models of the two proton treatment heads located at our facility using the aforementioned correction methods presented in this paper.

The goal of this work was to develop a mathematical model to predict Kaplan–Meier survival curves for chemotherapy combined with radiation in Non-Small Cell Lung Cancer patients for use in clinical trial design. The Gompertz model was used to describe tumor growth, radiation effect was simulated by the linear-quadratic model with an α/β-ratio of 10, and chemotherapy effect was based on the log-cell kill model. To account for repopulation during treatment, we considered two independent methods: 1) kickoff-repopulation using exponential growth with a decreased volume doubling time, or 2) Gompertz-repopulation using the gradually accelerating growth rate with tumor shrinkage. The input parameters were independently estimated by fitting to the SEER database for untreated tumors, RTOG-8808 for radiation only, and RTOG-9410 for sequential chemo-radiation. Applying the model, the benefit from concurrent chemo-radiation comparing to sequential for stage III patients was predicted to be a 6.6% and 6.2% improvement in overall survival for 3 and 5-years respectively, comparing well to the 5.3% and 4.5% observed in RTOG-9410. In summary, a mathematical model was developed to model tumor growth over extended periods of time, and can be used for the optimization of combined chemo-radiation scheduling and sequencing.

Several studies have shown that pediatric patients have an increased risk of developing a secondary malignancy several decades after treatment with radiotherapy and chemotherapy. In this work, we use a biologically motivated mathematical formalism to estimate the relative risks of breast, lung and thyroid cancers in childhood cancer survivors due to concurrent therapy regimen. This model specifically includes possible organ-specific interaction between radiotherapy and chemotherapy. The model predicts relative risks for developing secondary cancers after chemotherapy in breast, lung and thyroid tissues, and compared with the epidemiological data. For a concurrent therapy protocol, our model predicted relative risks of 3.2, 9.3, 4.5 as compared to the clinical data, i.e., 1.4, 8.0, 2.3 for secondary breast, lung and thyroid cancer risks, respectively. The extracted chemotherapy mutation induction rates for breast, lung and thyroid are 10−9, 0.5 × 10−6, 0.9 × 10−7 respectively. We found that there exists no synergistic interaction between radiation and chemotherapy for neither mutation induction nor cell kill in lung tissue, but there is an interaction in cell kill for the breast and thyroid organs. These findings help understand the risks of current clinical protocols and might provide rational guidance to develop future multi-modality treatment protocols to minimize secondary cancer risks.

Abstract Purpose Our group previously introduced an in vivo proton range verification methodology in which a silicon diode array system is used to correlate the dose rate profile per range modulation wheel cycle of the detector signal to the water‐equivalent path length (WEPL) for passively scattered proton beam delivery. The implementation of this system requires a set of calibration data to establish a beam‐specific response to WEPL fit for the selected ‘scout’ beam (a 1 cm overshoot of the predicted detector depth with a dose of 4 cGy) in water‐equivalent plastic. This necessitates a separate set of measurements for every ‘scout’ beam that may be appropriate to the clinical case. The current study demonstrates the use of Monte Carlo simulations for calibration of the time‐resolved diode dosimetry technique. Methods: Measurements for three ‘scout’ beams were compared against simulated detector response with Monte Carlo methods using the Tool for Particle Simulation (TOPAS). The ‘scout’ beams were then applied in the simulation environment to simulated water‐equivalent plastic, a CT of water‐equivalent plastic, and a patient CT data set to assess uncertainty. Results: Simulated detector response in water‐equivalent plastic was validated against measurements for ‘scout’ spread out Bragg peaks of range 10 cm, 15 cm, and 21 cm (168 MeV, 177 MeV, and 210 MeV) to within 3.4 mm for all beams, and to within 1 mm in the region where the detector is expected to lie. Conclusion: Feasibility has been shown for performing the calibration of the detector response for three ‘scout’ beams through simulation for the time‐resolved diode dosimetry technique in passive scattered proton delivery.