Optimal radiotherapy dose schedules under parametric uncertainty

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.

[1]  H. Withers,et al.  Alpha/beta value and the importance of size of dose per fraction for late complications in the supraglottic larynx. , 1986, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[2]  Fatemeh Saberian,et al.  A two-variable linear program solves the standard linear-quadratic formulation of the fractionation problem in cancer radiotherapy , 2015, Oper. Res. Lett..

[3]  Albert van der Kogel,et al.  Basic Clinical Radiobiology Fourth Edition , 2009 .

[4]  Thomas Ottmann,et al.  Algorithms for Reporting and Counting Geometric Intersections , 1979, IEEE Transactions on Computers.

[5]  A. Brahme,et al.  Optimal radiation beam profiles considering the stochastic process of patient positioning in fractionated radiation therapy , 1995 .

[6]  J. M. Taylor,et al.  Late normal tissue sequelae from radiation therapy for carcinoma of the tonsil: patterns of fractionation study of radiobiology. , 1995, International journal of radiation oncology, biology, physics.

[7]  J. Unkelbach,et al.  Inclusion of organ movements in IMRT treatment planning via inverse planning based on probability distributions. , 2004, Physics in medicine and biology.

[8]  F. Saberian,et al.  Optimal fractionation in radiotherapy with multiple normal tissues. , 2016, Mathematical medicine and biology : a journal of the IMA.

[9]  Randall K Ten Haken,et al.  Physical models and simpler dosimetric descriptors of radiation late toxicity. , 2007, Seminars in radiation oncology.

[10]  T. Bortfeld,et al.  The dependence of optimal fractionation schemes on the spatial dose distribution. , 2013, Physics in medicine and biology.

[11]  R. Onimaru,et al.  A mathematical study to select fractionation regimen based on physical dose distribution and the linear-quadratic model. , 2012, International journal of radiation oncology, biology, physics.

[12]  L. Haydu,et al.  Inter- and Intra-Patient Heterogeneity of Response and Progression to Targeted Therapy in Metastatic Melanoma , 2014, PloS one.

[13]  D A Pierre,et al.  Optimization Theory with Applications , 1986 .

[14]  J. Buatti,et al.  Calculation of cranial nerve complication probability for acoustic neuroma radiosurgery. , 2000, International journal of radiation oncology, biology, physics.

[15]  H D Thames,et al.  Repair capacity and kinetics of human skin during fractionated radiotherapy: erythema, desquamation, and telangiectasia after 3 and 5 year's follow-up. , 1989, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[16]  F. Collins,et al.  The path to personalized medicine. , 2010, The New England journal of medicine.

[17]  A. Ruifrok,et al.  Fractionation sensitivity of the rat cervical spinal cord during radiation retreatment. , 1992, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[18]  C. Belka,et al.  Dying cell clearance and its impact on the outcome of tumor radiotherapy , 2012, Front. Oncol..

[19]  Steven J. M. Jones,et al.  Comprehensive molecular portraits of human breast tumors , 2012, Nature.

[20]  M. Goitein,et al.  Tolerance of normal tissue to therapeutic irradiation. , 1991, International journal of radiation oncology, biology, physics.

[21]  Fatemeh Saberian,et al.  Corrigendum to "A two-variable linear program solves the standard linear-quadratic formulation of the fractionation problem in cancer radiotherapy" [Oper. Res. Lett. 43(3) (2015) 254-258] , 2015, Oper. Res. Lett..

[22]  S. Tucker,et al.  Isoeffect models and fractionated radiation therapy. , 1987, International journal of radiation oncology, biology, physics.

[23]  Timothy C Y Chan,et al.  Accounting for range uncertainties in the optimization of intensity modulated proton therapy , 2007, Physics in medicine and biology.

[24]  U Oelfke,et al.  Worst case optimization: a method to account for uncertainties in the optimization of intensity modulated proton therapy , 2008, Physics in medicine and biology.

[25]  R. Muschel,et al.  Effects of ionizing radiation on cell cycle progression , 1995, Radiation and environmental biophysics.

[26]  F. Fazio,et al.  Significant improvement in normal tissue sparing and target coverage for head and neck cancer by means of helical tomotherapy. , 2006, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[27]  C G Orton,et al.  A unified approach to dose-effect relationships in radiotherapy. I: Modified TDF and linear quadratic equations. , 1988, International journal of radiation oncology, biology, physics.

[28]  A. Cmelak,et al.  Improved survival of patients with human papillomavirus-positive head and neck squamous cell carcinoma in a prospective clinical trial. , 2008, Journal of the National Cancer Institute.

[29]  C C Ling,et al.  Using computerized video time lapse for quantifying cell death of X-irradiated rat embryo cells transfected with c-myc or c-Ha-ras. , 1999, Cancer research.

[30]  Xavier Geets,et al.  Adaptive radiotherapy of head and neck cancer. , 2010, Seminars in radiation oncology.

[31]  E. Holland,et al.  Optimization of radiation dosing schedules for proneural glioblastoma , 2015, Journal of Mathematical Biology.

[32]  J C Stroom,et al.  Inclusion of geometrical uncertainties in radiotherapy treatment planning by means of coverage probability. , 1999, International journal of radiation oncology, biology, physics.

[33]  K. Leder,et al.  Robust and probabilistic optimization of dose schedules in radiotherapy , 2015, 1503.00399.

[34]  J. Fowler,et al.  Sensitivity of human squamous cell carcinoma of the larynx to fractionated radiotherapy. , 1993, The British journal of radiology.

[35]  M. Mackenzie,et al.  Sparing the parotid glands and surgically transferred submandibular gland with helical tomotherapy in post-operative radiation of head and neck cancer: a planning study. , 2007, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[36]  D. Haussler,et al.  The Somatic Genomic Landscape of Glioblastoma , 2013, Cell.

[37]  A. Brahme,et al.  An adaptive control algorithm for optimization of intensity modulated radiotherapy considering uncertainties in beam profiles, patient set-up and internal organ motion. , 1998, Physics in medicine and biology.

[38]  S. Henderson,et al.  Robust optimization for intensity modulated radiation therapy treatment planning under uncertainty , 2005, Physics in medicine and biology.

[39]  John N Tsitsiklis,et al.  A robust approach to IMRT optimization , 2006, Physics in medicine and biology.

[40]  J. M. Taylor,et al.  Dose fractionation and regeneration in radiotherapy for cancer of the oral cavity and oropharynx. Part 2. Normal tissue responses: acute and late effects. , 1990, International journal of radiation oncology, biology, physics.