Glioblastomas are the most aggressive primary brain tumor. Despite treatment with surgery, radiation and chemotherapy, these tumors remain uncurable and few significant increases in survival have been observed over the last half-century. We recently employed a combined theoretical and experimental approach to predict the effectiveness of radiation administration schedules, identifying two schedules that led to superior survival in a mouse model of the disease (Leder et al., Cell 156 (3):603-616, 2014). Here we extended this approach to consider fractionated schedules to best minimize toxicity arising in early-and late-responding tissues. To this end, we decomposed the problem into two separate solvable optimization tasks: (i) optimization of the amount of radiation per dose, and (ii) optimization of the amount of time that 123 H. Badri et al. passes between radiation doses. To ensure clinical applicability, we then considered the impact of clinical operating hours by incorporating time constraints consistent with operational schedules of the radiology clinic. We found that there was no significant loss incurred by restricting dosage to an 8:00 a.m. to 5:00 p.m. window. Our flexible approach is also applicable to other tumor types treated with radiotherapy.Keywords Brain tumors · Radiotherapy · Nonlinear programming · Linear-quadratic model Mathematics Subject Classification 90C11 · 90C26 · 90C30 · 90C90 · 65K05

We consider the effects of parameter uncertainty on the optimal radiation schedule in the context of the linear-quadratic model. Our interest arises from the observation that if inter-patient variability in normal and tumor tissue radiosensitivity or sparing factor of the organs-at-risk (OAR) are not accounted for during radiation scheduling, the performance of the therapy may be strongly degraded or the OAR may receive a substantially larger dose than the allowable threshold. This paper proposes a stochastic radiation scheduling concept to incorporate inter-patient variability into the scheduling optimization problem. Our method is based on a probabilistic approach, where the model parameters are given by a set of random variables. Our probabilistic formulation ensures that our constraints are satisfied with a given probability, and that our objective function achieves a desired level with a stated probability. We used a variable transformation to reduce the resulting optimization problem to two dimensions. We showed that the optimal solution lies on the boundary of the feasible region and we implemented a branch and bound algorithm to find the global optimal solution. We demonstrated how the configuration of optimal schedules in the presence of uncertainty compares to optimal schedules in the absence of uncertainty (conventional schedule). We observed that in order to protect against the possibility of the model parameters falling into a region where the conventional schedule is no longer feasible, it is required to avoid extremal solutions, i.e. a single large dose or very large total dose delivered over a long period. Finally, we performed numerical experiments in the setting of head and neck tumors including several normal tissues to reveal the effect of parameter uncertainty on optimal schedules and to evaluate the sensitivity of the solutions to the choice of key model parameters.

Metastasis is the process by which cells from a primary tumor disperse and form new tumors at distant anatomical locations. The treatment and prevention of metastatic cancer remains an extremely challenging problem. This work introduces a novel biologically motivated objective function to the radiation optimization community that takes into account metastatic risk instead of the status of the primary tumor. In this work, we consider the problem of developing fractionated irradiation schedules that minimize production of metastatic cancer cells while keeping normal tissue damage below an acceptable level. A dynamic programming framework is utilized to determine the optimal fractionation scheme. We evaluated our approach on a breast cancer case using the heart and the lung as organs-at-risk (OAR). For small tumor α/β values, hypo-fractionated schedules were optimal, which is consistent with standard models. However, for relatively larger α/β values, we found the type of schedule depended on various parameters such as the time when metastatic risk was evaluated, the α/β values of the OARs, and the normal tissue sparing factors. Interestingly, in contrast to standard models, hypo-fractionated and semihypo-fractionated schedules (large initial doses with doses tapering off with time) were suggested even with large tumor α/β values. Numerical results indicate potential for significant reduction in metastatic risk.

No abstract

In this work we review past articles that have mathematically studied cancer heterogeneity and the impact of this heterogeneity on the structure of optimal therapy. We look at past works on modeling how heterogeneous tumors respond to radiotherapy, and take a particularly close look at how the optimal radiotherapy schedule is modified by the presence of heterogeneity. In addition, we review past works on the study of optimal chemotherapy when dealing with heterogeneous tumors.Reviewers: This article was reviewed by Thomas McDonald, David Axelrod, and Leonid Hanin.