Robust treatment planning with conditional value at risk chance constraints in intensity‐modulated proton therapy

Background and purpose: Intensity‐modulated proton therapy (IMPT) is highly sensitive to range uncertainties and uncertainties caused by setup variation. The conventional inverse treatment planning of IMPT based on the planning target volume (PTV) is not often sufficient to ensure robustness of treatment plans. We applied a probabilistic framework (chance‐constrained optimization) in IMPT planning to hedge against the influence of uncertainties. Material and methods: We retrospectively selected one patient with lung cancer, one patient with head and neck (H&N) cancer, and one with prostate cancer for this analysis. Using their original images and prescriptions, we created new IMPT plans using two methods: (1) a robust chance‐constrained treatment planning method with the clinical target volume (CTV) as the target; (2) the margin‐based method with PTV as the target, which was solved by commercial software, CPLEX, using linear programming. For the first method, we reformulated the model into a tractable mixed‐integer programming problem and sped up the calculation using Benders decomposition. The dose‐volume histograms (DVHs) from the nominal and perturbed dose distributions were used to assess and compare plan quality. DVHs for all uncertain scenarios along with the nominal DVH were plotted. The width of the “bands” of DVHs was used to quantify the plan sensitivity to uncertainty. The newly developed Benders decomposition method was compared with a commercial solution to demonstrate its computational efficiency. The trade‐off between nominal plan quality and plan robustness was investigated. Results: Our chance‐constrained model outperformed the PTV method in terms of tumor coverage, tumor dose homogeneity, and plan robustness. Our model was shown to produce IMPT plans to meet the dose‐volume constraints of organs at risk (OARs) and had better sparing of OARs than the PTV method in the three clinical cases included in this study. The chance‐constrained model provided a flexible tool for users to balance between plan robustness and plan quality. In addition, our in‐house developed method was found to be much faster than the commercial solution. Conclusion: With explicit control of plan robustness, the chance‐constrained robust optimization model generated superior IMPT plans compared to the PTV‐based method.

[1]  Wei Liu,et al.  Preliminary evaluation of multifield and single-field optimization for the treatment planning of spot-scanning proton therapy of head and neck cancer. , 2013, Medical physics.

[2]  Wei Liu,et al.  Dosimetric benefits of robust treatment planning for intensity modulated proton therapy for base-of-skull cancers. , 2014, Practical radiation oncology.

[3]  Wei Chen,et al.  Including robustness in multi-criteria optimization for intensity-modulated proton therapy , 2011, Physics in medicine and biology.

[4]  E. Pedroni,et al.  Intensity modulated proton therapy: a clinical example. , 2001, Medical physics.

[5]  Anders Forsgren,et al.  Minimax optimization for handling range and setup uncertainties in proton therapy. , 2011, Medical physics.

[6]  Wei Liu,et al.  Influence of robust optimization in intensity-modulated proton therapy with different dose delivery techniques. , 2012, Medical physics.

[7]  Wei Liu,et al.  Exploratory Study of 4D versus 3D Robust Optimization in Intensity Modulated Proton Therapy for Lung Cancer. , 2016, International journal of radiation oncology, biology, physics.

[8]  Radhe Mohan,et al.  Intensity-modulated proton therapy reduces the dose to normal tissue compared with intensity-modulated radiation therapy or passive scattering proton therapy and enables individualized radical radiotherapy for extensive stage IIIB non-small-cell lung cancer: a virtual clinical study. , 2010, International journal of radiation oncology, biology, physics.

[9]  Xiaodong Zhang,et al.  Beyond Gaussians: a study of single-spot modeling for scanning proton dose calculation , 2012, Physics in medicine and biology.

[10]  Christodoulos A. Floudas Generalized Benders Decomposition , 2009, Encyclopedia of Optimization.

[11]  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.

[12]  Albin Fredriksson,et al.  A characterization of robust radiation therapy treatment planning methods-from expected value to worst case optimization. , 2012, Medical physics.

[13]  Rasmus Bokrantz,et al.  A critical evaluation of worst case optimization methods for robust intensity-modulated proton therapy planning. , 2014, Medical physics.

[14]  A. Lomax,et al.  Intensity modulated proton therapy and its sensitivity to treatment uncertainties 1: the potential effects of calculational uncertainties , 2008, Physics in medicine and biology.

[15]  Bo Zeng,et al.  Chance Constrained Mixed Integer Program: Bilinear and Linear Formulations, and Benders Decomposition , 2014, 1403.7875.

[16]  Steve B. Jiang,et al.  An experimental investigation on intra-fractional organ motion effects in lung IMRT treatments. , 2003, Physics in medicine and biology.

[17]  Radhe Mohan,et al.  Proton stereotactic body radiation therapy for clinically challenging cases of centrally and superiorly located stage I non-small-cell lung cancer. , 2011, International journal of radiation oncology, biology, physics.

[18]  A J Lomax,et al.  Intensity modulated proton therapy and its sensitivity to treatment uncertainties 2: the potential effects of inter-fraction and inter-field motions , 2008, Physics in medicine and biology.

[19]  Thomas A. Henzinger,et al.  Probabilistic programming , 2014, FOSE.

[20]  Thomas Bortfeld,et al.  Reducing the sensitivity of IMPT treatment plans to setup errors and range uncertainties via probabilistic treatment planning. , 2008, Medical physics.

[21]  Radhe Mohan,et al.  Robust optimization of intensity modulated proton therapy. , 2012, Medical physics.

[22]  Wei Liu,et al.  Impact of respiratory motion on worst-case scenario optimized intensity modulated proton therapy for lung cancers. , 2015, Practical radiation oncology.

[23]  Gino J. Lim,et al.  A two-phase method for selecting IMRT treatment beam angles: Branch-and-Prune and local neighborhood search , 2012, Eur. J. Oper. Res..

[24]  Radhe Mohan,et al.  Effectiveness of robust optimization in intensity-modulated proton therapy planning for head and neck cancers. , 2013, Medical physics.

[25]  Thomas Bortfeld,et al.  Evaluation of Dosimetric Gain and Uncertainties in Proton Therapy Delivery with Scanned Pencil Beam in Treatment of Base-of-skull and Spinal Tumors , 2010 .

[26]  Arvind Kumar,et al.  A New Linear Programming Approach to Radiation Therapy Treatment Planning Problems , 2006, Oper. Res..

[27]  Steven J Frank,et al.  PTV-based IMPT optimization incorporating planning risk volumes vs robust optimization. , 2013, Medical physics.

[28]  Wei Liu,et al.  Uncertainty incorporated beam angle optimization for IMPT treatment planning. , 2012, Medical physics.

[29]  Joe Y. Chang,et al.  4D Proton treatment planning strategy for mobile lung tumors. , 2007, International journal of radiation oncology, biology, physics.

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

[31]  Houra Mahmoudzadeh,et al.  A robust-CVaR optimization approach with application to breast cancer therapy , 2014, Eur. J. Oper. Res..

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