Impact of respiratory motion on variable relative biological effectiveness in 4D-dose distributions of proton therapy.

BACKGROUND: Organ motion during radiation therapy with scanned protons leads to deviations between the planned and the delivered physical dose. Using a constant relative biological effectiveness (RBE) of 1.1 linearly maps these deviations into RBE-weighted dose. However, a constant value cannot account for potential nonlinear variations in RBE suggested by variable RBE models. Here, we study the impact of motion on recalculations of RBE-weighted dose distributions using a phenomenological variable RBE model. MATERIAL AND METHODS: 4D-dose calculation including variable RBE was implemented in the open source treatment planning toolkit matRad. Four scenarios were compared for one field and two field proton treatments for a liver cancer patient assuming (α∕ß)x = 2 Gy and (α∕ß)x = 10 Gy: (A) the optimized static dose distribution with constant RBE, (B) a static recalculation with variable RBE, (C) a 4D-dose recalculation with constant RBE and (D) a 4D-dose recalculation with variable RBE. For (B) and (D), the variable RBE was calculated by the model proposed by McNamara. For (C), the physical dose was accumulated with direct dose mapping; for (D), dose-weighted radio-sensitivity parameters of the linear quadratic model were accumulated to model synergistic irradiation effects on RBE. RESULTS: Dose recalculation with variable RBE led to an elevated biological dose at the end of the proton field, while 4D-dose recalculation exhibited random deviations everywhere in the radiation field depending on the interplay of beam delivery and organ motion. For a single beam treatment assuming (α∕ß)x = 2 Gy, D95% was 1.98 Gy (RBE) (A), 2.15 Gy (RBE) (B), 1.81 Gy (RBE) (C) and 1.98 Gy (RBE) (D). The homogeneity index was 1.04 (A), 1.08 (B), 1.23 (C) and 1.25 (D). CONCLUSION: For the studied liver case, intrafractional motion did not reduce the modulation of the RBE-weighted dose postulated by variable RBE models for proton treatments.

Development of the open-source dose calculation and optimization toolkit matRad.

PURPOSE: We report on the development of the open-source cross-platform radiation treatment planning toolkit matRad and its comparison against validated treatment planning systems. The toolkit enables three-dimensional intensity-modulated radiation therapy treatment planning for photons, scanned protons and scanned carbon ions. METHODS: matRad is entirely written in Matlab and is freely available online. It re-implements well-established algorithms employing a modular and sequential software design to model the entire treatment planning workflow. It comprises core functionalities to import DICOM data, to calculate and optimize dose as well as a graphical user interface for visualization. matRad dose calculation algorithms (for carbon ions this also includes the computation of the relative biological effect) are compared against dose calculation results originating from clinically approved treatment planning systems. RESULTS: We observe three-dimensional γ-analysis pass rates ≥ 99.67% for all three radiation modalities utilizing a distance to agreement of 2 mm and a dose difference criterion of 2%. The computational efficiency of matRad is evaluated in a treatment planning study considering three different treatment scenarios for every radiation modality. For photons, we measure total run times of 145 s-1260 s for dose calculation and fluence optimization combined considering 4-72 beam orientations and 2608-13597 beamlets. For charged particles, we measure total run times of 63 s-993 s for dose calculation and fluence optimization combined considering 9963-45574 pencil beams. Using a CT and dose grid resolution of 0.3 cm3 requires a memory consumption of 1.59 GB-9.07 GB and 0.29 GB-17.94 GB for photons and charged particles, respectively. CONCLUSION: The dosimetric accuracy, computational performance and open-source character of matRad encourages a future application of matRad for both educational and research purposes.

Worst case optimization for interfractional motion mitigation in carbon ion therapy of pancreatic cancer.

INTRODUCTION: The efficacy of radiation therapy treatments for pancreatic cancer is compromised by abdominal motion which limits the spatial accuracy for dose delivery - especially for particles. In this work we investigate the potential of worst case optimization for interfractional offline motion mitigation in carbon ion treatments of pancreatic cancer. METHODS: We implement a worst case optimization algorithm that explicitly models the relative biological effectiveness of carbon ions during inverse planning. We perform a comparative treatment planning study for seven pancreatic cancer patients. Treatment plans that have been generated using worst case optimization are compared against (1) conventional intensity-modulated carbon ion therapy, (2) single field uniform dose carbon ion therapy, and (3) an ideal yet impractical scenario relying on daily re-planning. The dosimetric quality and robustness of the resulting treatment plans is evaluated using reconstructions of the daily delivered dose distributions on fractional control CTs. RESULTS: Idealized daily re-planning consistently gives the best dosimetric results with regard to both target coverage and organ at risk sparing. The absolute reduction of D 95 within the gross tumor volume during fractional dose reconstruction is most pronounced for conventional intensity-modulated carbon ion therapy. Single field uniform dose optimization exhibits no substantial reduction for six of seven patients and values for D 95 for worst case optimization fall in between. The treated volume (D>95 % prescription dose) outside of the gross tumor volume is reduced by a factor of two by worst case optimization compared to conventional optimization and single field uniform dose optimization. Single field uniform dose optimization comes at an increased radiation exposure of normal tissues, e.g. ≈2 Gy (RBE) in the mean dose in the kidneys compared to conventional and worst case optimization and ≈4 Gy (RBE) in D 1 in the spinal cord compared to worst case optimization. CONCLUSION: Interfractional motion substantially deteriorates dose distributions for carbon ion treatments of pancreatic cancer patients. Single field uniform dose optimization mitigates the negative influence of motion on target coverage at an increased radiation exposure of normal tissue. Worst case optimization enables an exploration of the trade-off between robust target coverage and organ at risk sparing during inverse treatment planning beyond margin concepts.

Feasibility of dose painting using volumetric modulated arc optimization and delivery.

PURPOSE: Dose painting strategies are limited by optimization algorithms in treatment planning systems and physical constraints of the beam delivery. We investigate dose conformity using the RapidArc optimizer and beam delivery technique. Furthermore, robustness of the plans with respect to positioning uncertainties are evaluated. METHODS: A head & neck cancer patient underwent a [(61)Cu]Cu-ATSM PET/CT-scan. PET-SUVs were converted to prescribed dose with a base dose of 60 Gy, and target mean dose 90 Gy. The voxel-based prescription was converted into 3, 5, 7, 9, and 11 discrete prescription levels. Optimization was performed in Eclipse, varying the following parameters: MLC leaf width (5 mm and 2.5 mm), number of arcs (1 and 2) and collimator rotation (0, 15, 30 and 45 degrees). Dose conformity was evaluated using quality volume histograms (QVHs), and relative volumes receiving within ±5% of prescribed dose (Q(0.95-1.05)). Deliverability was tested using a Delta4(®) phantom. Robustness was tested by shifting the isocenter 1 mm and 2 mm in all directions, and recalculating the dose. RESULTS: Good conformity was obtained using MLC leaf width 2.5 mm, two arcs, and collimators 45/315 degrees, with Q(0.95-1.05)=92.8%, 91.6%, 89.7% and 84.6%. Using only one arc or increasing the MLC leaf width had a small deteriorating effect of 2-5%. Small changes in collimator angle gave small changes, but large changes in collimator angle gave a larger decrease in plan conformity; for angles of 15 and 0 degrees (two arcs, 2.5 mm leaf width), Q(0.95-1.05) decreased by up to 15%. Consistency between planned and delivered dose was good, with ∼90% of gamma values <1. For 1 mm shift, Q(0.95-1.05) was decreased by 5-15%, while for 2 mm shift, Q(0.95-1.05) was decreased to 55-60%. CONCLUSIONS: Results demonstrate feasibility of planning of prescription doses with multiple levels for dose painting using RapidArc, and plans were deliverable. Robustness to positional error was low.

Comparison of arc-modulated cone beam therapy and helical tomotherapy for three different types of cancer.

PURPOSE: Arc-modulated cone beam therapy (AMCBT) is a fast treatment technique deliverable in a single rotation with a conventional C-arm shaped linac. In this planning study, the authors assess the dosimetric properties of single-arc therapy in comparison to helical tomotherapy for three different tumor types. METHODS: Treatment plans for three patients with prostate carcinoma, three patients with anal cancer, and three patients with head and neck cancer were optimized for helical tomotherapy and AMCBT. The dosimetric comparison of the two techniques is based on physical quantities derived from dose-volume histograms. RESULTS: For prostate cancer, the quality of dose distributions calculated for AMCBT was of equal quality as that generated for tomotherapy with the additional benefits of a faster delivery and a lower integral dose. For highly complex geometries, the plan quality achievable with helical tomotherapy could not be achieved with arc-modulated cone beam therapy. CONCLUSIONS: Rotation therapy with a conventional linac in a single arc is capable to deliver a high and homogeneous dose to the target and spare organs at risk. Advantages of this technique are a fast treatment time and a lower integral dose in comparison to helical tomotherapy. For highly complex cases, e.g., with several target regions, the dose shaping capabilities of AMCBT are inferior to those of tomotherapy. However, treatment plans for AMCBT were also clinically acceptable.

Development of an optimization concept for arc-modulated cone beam therapy.

In this paper, we propose an optimization concept for a rotation therapy technique which is referred to as arc-modulated cone beam therapy (AMCBT). The aim is a reduction of the treatment time while achieving a treatment plan quality equal to or better than that of IMRT. Therefore, the complete dose is delivered in one single gantry rotation and the beam is modulated by a multileaf collimator. The degrees of freedom are the field shapes and weights for a predefined number of beam directions. In the new optimization loop, the beam weights are determined by a gradient algorithm and the field shapes by a tabu search algorithm. We present treatment plans for AMCBT for two clinical cases. In comparison to step-and-shoot IMRT treatment plans, it was possible by AMCBT to achieve dose distributions with a better dose conformity to the target and a lower mean dose for the most relevant organ at risk. Furthermore, the number of applied monitor units was reduced for AMCBT in comparison to IMRT treatment plans.

Trehalose 6-phosphate regulates starch synthesis via posttranslational redox activation of ADP-glucose pyrophosphorylase.

Trehalose is the most widespread disaccharide in nature, occurring in bacteria, fungi, insects, and plants. Its precursor, trehalose 6-phosphate (T6P), is also indispensable for the regulation of sugar utilization and growth, but the sites of action are largely unresolved. Here we use genetic and biochemical approaches to investigate whether T6P acts to regulate starch synthesis in plastids of higher plants. Feeding of trehalose to Arabidopsis leaves led to stimulation of starch synthesis within 30 min, accompanied by activation of ADP-glucose pyrophosphorylase (AGPase) via posttranslational redox modification. The response resembled sucrose but not glucose feeding and depended on the expression of SNF1-related kinase. We also analyzed transgenic Arabidopsis plants with T6P levels increased by expression of T6P synthase or decreased by expression of T6P phosphatase (TPP) in the cytosol. Compared with wild type, leaves of T6P synthase-expressing plants had increased redox activation of AGPase and increased starch, whereas TPP-expressing plants showed the opposite. Moreover, TPP expression prevented the increase in AGPase activation in response to sucrose or trehalose feeding. Incubation of intact isolated chloroplasts with 100 muM T6P significantly and specifically increased reductive activation of AGPase within 15 min. Results provide evidence that T6P is synthesized in the cytosol and acts on plastidial metabolism by promoting thioredoxin-mediated redox transfer to AGPase in response to cytosolic sugar levels, thereby allowing starch synthesis to be regulated independently of light. The discovery informs about the evolution of plant metabolism and how chloroplasts of prokaryotic origin use an intermediate of the ancient trehalose pathway to report the metabolic status of the cytosol.