Modeling and Numerical Solution of a Cancer Therapy Optimal Control Problem
- 1 Institut für Mathematik, Universität Würzburg, Würzburg, Germany
- 2 Institut für Mathematik, Universität Würzburg, Würzburg, Germany
- 3 Lehrstuhl für Wissenschaftliches Rechnen, Fakultüt für Mathematik, Physik und Informatik, Universität Bayreuth, Bayreuth, Germany
- 4 Institut für Mathematik, Universität Würzburg, Würzburg, Germany
Abstract
A mathematical optimal-control tumor therapy framework consisting of radio- and anti-angiogenesis control strategies that are included in a tumor growth model is investigated. The governing system, resulting from the combination of two well established models, represents the differential constraint of a non-smooth optimal control problem that aims at reducing the volume of the tumor while keeping the radio- and anti-angiogenesis chemical dosage to a minimum. Existence of optimal solutions is proved and necessary conditions are formulated in terms of the Pontryagin maximum principle. Based on this principle, a so-called sequential quadratic Hamiltonian (SQH) method is discussed and benchmarked with an “interior point optimizer—a mathematical programming language” (IPOPT-AMPL) algorithm. Results of numerical experiments are presented that successfully validate the SQH solution scheme. Further, it is shown how to choose the optimisation weights in order to obtain treatment functions that successfully reduce the tumor volume to zero.
- Azizi, M. and Hesameddini, E. (2017) Grünwald-Letnikov Scheme for System of Chronic Myelogenous Leukemia Fractional Differential Equations and Its Optimal Control of Drug Treatment. Journal of Mahani Mathematical Research Center, 5, 51-57.
- Moore, H. (2018) How to Mathematically Optimize Drug Regimens Using Optimal Control. Journal of Pharmacokinetics and Pharmacodynamics, 45, 127-137. https://doi.org/10.1007/s10928-018-9568-y
- Schättler, H. and Ledzewicz, U. (2015) Optimal Control for Mathematical Models of Cancer Therapies. Springer, Berlin/Heidelberg. https://doi.org/10.1007/978-1-4939-2972-6
- Chudej, K., Huebner, D. and Pesch, H.J. (2016) Numerische Lösung eines mathematischen Modells für eine optimale Krebskombinationstherapie aus Anti-Angiogenese und Strahlentherapie. In: Tagungsband ASIM 2016-23 Symposium Simulationstechnik, Dresden, Volume 52 of Argesim Report, ARGESIM Verlag, Wien, Austria, 169-176.
- Hahnfeldt, P., Panigrahy, D., Folkman, J. and Hlatky, L. (1999) Tumor Development under Angiogenic Signaling. Cancer Research, 59, 4770-4775.
- Ergun, A., Camphausen, K. and Wein, L.M. (2003) Optimal Scheduling of Radiotherapy and Angiogenic Inhibitors. Bulletin of Mathematical Biology, 65, 407-424. https://doi.org/10.1016/S0092-8240(03)00006-5
- Wein, L.M., Cohen, J.E. and Wu, J.T. (2000) Dynamic Optimization of a Linear-Quadratic Model with Incomplete Repair and Volume-Dependent Sensitivity and Repopulation. International Journal of Radiation Oncology, Biology, Physics, 47, 1073-1083. https://doi.org/10.1016/S0360-3016(00)00534-4
- Bellomo, N. and Preziosi, L. (2000) Modelling and Mathematical Problems Related to Tumor Evolution and Its Interaction with the Immune System. Mathematical and Computer Modelling, 32, 413-452. https://doi.org/10.1016/S0895-7177(00)00143-6
- Chudej, K., Wagner, L. and Pesch, H.J. (2015) Numerical Solution of an Optimal Control Problem in Cancer Treatment: Combined Radio and Anti-Angiogenic Therapy. IFAC-PapersOnLine, 48, 665-666. https://doi.org/10.1016/j.ifacol.2015.05.082
- Sontag, E.D. (2013) Mathematical Control Theory: Deterministic Finite Dimensional Systems. Volume 6. Springer Science & Business Media, Berlin/Heidelberg.
- Tröltzsch, F. (2010) Optimal Control of Partial Differential Equations. Graduate Studies in Mathematics. American Mathematical Society, Providence.
- Adams, R.A. and Fournier, J.J.F. (2003) Sobolev Spaces. Volume 140. Academic Press, Cambridge.