We analyze the effect of tumor repopulation on optimal dose delivery in radiation therapy. We are primarily motivated by accelerated tumor repopulation towards the end of radiation treatment, which is believed to play a role in treatment failure for some tumor sites. A dynamic programming framework is developed to determine an optimal fractionation scheme based on a model of cell kill due to radiation and tumor growth in between treatment days. We find that faster tumor growth suggests shorter overall treatment duration. In addition, the presence of accelerated repopulation suggests larger dose fractions later in the treatment to compensate for the increased tumor proliferation. We prove that the optimal dose fractions are increasing over time. Numerical simulations indicate potential for improvement in treatment effectiveness.
Perhaps the closest related work is @cite_36 , which considers both faster tumor proliferation and re-oxygenation during the course of treatment. While a dose intensification strategy is also suggested in @cite_36 , the primary rationale for increasing dose fractions is different: it is concluded that due to the increase in tumor sensitivity from re-oxygenation, larger fraction sizes are more effective at the end of treatment. Our work, on the other hand, suggests dose intensification (i.e., larger doses over time) as a direct consequence of a model of accelerated tumor repopulation during the course of treatment.
@cite_36
2136543444
Purpose: The linear-quadratic model typically assumes that tumor sensitivity and repopulation are constant over the time course of radiotherapy. However, evidence suggests that the growth fraction increases and the cell-loss factor decreases as the tumor shrinks. We investigate whether this evolution in tumor geometry, as well as the irregular time intervals between fractions in conventional hyperfractionation schemes, can be exploited by fractionation schedules that employ time-varying fraction sizes. Methods: We construct a mathematical model of a spherical tumor with a hypoxic core and a viable rim, which is most appropriate for a prevascular tumor, and is only a caricature of a vascularized tumor. This model is embedded into the traditional linear-quadratic model by assuming instantaneous reoxygenation. Dynamic programming is used to numerically compute the fractionation regimen that maximizes the tumor-control probability (TCP) subject to constraints on the biologically effective dose of the early and late tissues. Results: In several numerical examples that employ five or 10 fractions per week on a 1-cm or 5-cm diameter tumor, optimally varying the fraction sizes increases the TCP significantly. The optimal regimen incorporates large Friday (afternoon, if 10 fractions per week) fractions that are escalated throughout the course of treatment, and larger afternoon fractions than morning fractions. Conclusion: Numerical results suggest that a significant increase in tumor cure can be achieved by allowing the fraction sizes to vary throughout the course of treatment. Several strategies deserve further investigation: using larger fractions before overnight and weekend breaks, and escalating the dose (particularly on Friday afternoons) throughout the course of treatment. © 2000 Elsevier Science Inc. Dynamic optimization, Linear‐ quadratic, Reoxygenation, Repair, Repopulation.