Optimal Control for Cancer Chemotherapy ODE Models: Potential of Optimal Schedules and Choice of Objective Function

In this article, four different mathematical models of chemotherapy from the literature are investigated with respect to optimal control of drug treatment schedules. The various models are based on two different sets of ordinary differential equations and contain either chemotherapy, immunotherapy, anti-angiogenic therapy or combinations of these. Optimal control problem formulations based on these models are proposed, discussed and compared. For different parameter sets, scenarios, and objective functions o ptimal control problems are solved numerically with Bock’s direct multiple shooting method. In particular, we show that an optimally controlled therapy can be the reason for the difference between a growing and a totally vanishing tumor in comparison to standard treatment schemes and untreated or wrongly treated tumors. Furthermore, we compare different objective functions. Eventually, we show that there is a high potential for optimization of chemotherapy schedules, although the currently available models are not yet appropriate for transferring the optimal therapies into medical practice due to patient- , cancer-, and therapy-specific components.

[1]  H. Bock,et al.  A Multiple Shooting Algorithm for Direct Solution of Optimal Control Problems , 1984 .

[2]  A. Diefenbach,et al.  Rae1 and H60 ligands of the NKG2D receptor stimulate tumour immunity , 2001, Nature.

[3]  Lisette G de Pillis,et al.  A validated mathematical model of cell-mediated immune response to tumor growth. , 2007, Cancer research.

[4]  Lorenz T. Biegler,et al.  On the implementation of an interior-point filter line-search algorithm for large-scale nonlinear programming , 2006, Math. Program..

[5]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[6]  Christof Büskens,et al.  Sensitivity Analysis and Real-Time Control of Parametric Optimal Control Problems Using Nonlinear Programming Methods , 2001 .

[7]  Anil V. Rao,et al.  Practical Methods for Optimal Control Using Nonlinear Programming , 1987 .

[8]  A. Radunskaya,et al.  Mixed Immunotherapy and Chemotherapy of Tumors: Modeling, Applications and Biological Interpretations , 2022 .

[9]  L. Wein,et al.  Optimal scheduling of radiotherapy and angiogenic inhibitors , 2003, Bulletin of mathematical biology.

[10]  H. Maurer,et al.  SQP-methods for solving optimal control problems with control and state constraints: adjoint variables, sensitivity analysis and real-time control , 2000 .

[11]  Jean Clairambault,et al.  Circadian timing in cancer treatments. , 2010, Annual review of pharmacology and toxicology.

[12]  A. Perelson,et al.  Nonlinear dynamics of immunogenic tumors: parameter estimation and global bifurcation analysis. , 1994, Bulletin of mathematical biology.

[13]  Irene Bauer Numerische Verfahren zur Lösung von Anfangswertaufgaben und zur Generierung von ersten und zweiten Ableitungen mit Anwendungen bei Optimierungsaufgaben in Chemie und Verfahrenstechnik , 1999 .

[14]  Ami Radunskaya,et al.  A mathematical tumor model with immune resistance and drug therapy: an optimal control approach , 2001 .

[15]  J. Folkman,et al.  Anti‐Angiogenesis: New Concept for Therapy of Solid Tumors , 1972, Annals of surgery.

[16]  L. G. de Pillis,et al.  OPTIMAL CONTROL OF MIXED IMMUNOTHERAPY AND CHEMOTHERAPY OF TUMORS , 2008 .

[17]  P. Hahnfeldt,et al.  Tumor development under angiogenic signaling: a dynamical theory of tumor growth, treatment response, and postvascular dormancy. , 1999, Cancer research.

[18]  H. Maurer,et al.  Optimization methods for the verification of second order sufficient conditions for bang–bang controls , 2005 .

[19]  H. Schättler,et al.  On optimal delivery of combination therapy for tumors. , 2009, Mathematical biosciences.

[20]  Urszula Ledzewicz,et al.  Bang-bang and singular controls in a mathematical model for combined anti-angiogenic and chemotherapy treatments , 2009, Proceedings of the 48h IEEE Conference on Decision and Control (CDC) held jointly with 2009 28th Chinese Control Conference.

[21]  Jean Clairambault,et al.  Modelling Physiological and Pharmacological Control on Cell Proliferation to Optimise Cancer Treatments , 2009 .

[22]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[23]  Alberto Gandolfi,et al.  Tumour eradication by antiangiogenic therapy: analysis and extensions of the model by Hahnfeldt et al. (1999). , 2004, Mathematical biosciences.

[24]  Johannes P. Schlöder,et al.  An efficient multiple shooting based reduced SQP strategy for large-scale dynamic process optimization. Part 1: theoretical aspects , 2003, Comput. Chem. Eng..

[25]  D. Kirschner,et al.  Modeling immunotherapy of the tumor – immune interaction , 1998, Journal of mathematical biology.

[26]  M. J. D. Powell,et al.  A fast algorithm for nonlinearly constrained optimization calculations , 1978 .

[27]  M. Raffeld,et al.  Cancer Regression and Autoimmunity in Patients After Clonal Repopulation with Antitumor Lymphocytes , 2002, Science.

[28]  Ami Radunskaya,et al.  The dynamics of an optimally controlled tumor model: A case study , 2003 .

[29]  S Chareyron,et al.  Mixed immunotherapy and chemotherapy of tumors: feedback design and model updating schemes. , 2009, Journal of theoretical biology.