Optimization of Dose Fractionation for Radiotherapy of a Solid Tumor with Account of Oxygen Effect and Proliferative Heterogeneity

A spatially-distributed continuous mathematical model of solid tumor growth and treatment by fractionated radiotherapy is presented. The model explicitly accounts for three time and space-dependent factors that influence the efficiency of radiotherapy fractionation schemes—tumor cell repopulation, reoxygenation and redistribution of proliferative states. A special algorithm is developed, aimed at finding the fractionation schemes that provide increased tumor cure probability under the constraints of maximum normal tissue damage and maximum fractional dose. The optimization procedure is performed for varied radiosensitivity of tumor cells under the values of model parameters, corresponding to different degrees of tumor malignancy. The resulting optimized schemes consist of two stages. The first stages are aimed to increase the radiosensitivity of the tumor cells, remaining after their end, sparing the caused normal tissue damage. This allows to increase the doses during the second stages and thus take advantage of the obtained increased radiosensitivity. Such method leads to significant expansions in the curative ranges of the values of tumor radiosensitivity parameters. Overall, the results of this study represent the theoretical proof of concept that non-uniform radiotherapy fractionation schemes may be considerably more effective that uniform ones, due to the time and space-dependent effects.