# American Institute of Mathematical Sciences

• Previous Article
Sampled–data model predictive control: Adaptive time–mesh refinement algorithms and guarantees of stability
• DCDS-B Home
• This Issue
• Next Article
A sufficient optimality condition for delayed state-linear optimal control problems
May  2019, 24(5): 2315-2334. doi: 10.3934/dcdsb.2019097

## Optimal control for a mathematical model for anti-angiogenic treatment with Michaelis-Menten pharmacodynamics

 1 Institute of Mathematics, Lodz University of Technology, 90-924 Lodz, Poland 2 Dept. of Mathematics and Statistics, Southern Illinois University Edwardsville, Edwardsville, Il, 62026-1653, USA 3 Dept. of Electrical and Systems Engineering, Washington University, St. Louis, Mo, 63130, USA

* Corresponding author: Urszula Ledzewicz

Received  January 2018 Revised  January 2019 Published  March 2019

We analyze the structure of optimal protocols for a mathematical model of tumor anti-angiogenic treatment. The control represents the concentration of the agent and we consider the problem to administer an a priori given total amount of agents in order to achieve a minimum tumor volume/maximum tumor reduction. In earlier work, this problem was studied with a log-kill type pharmacodynamic model for drug effects which does not account for saturation of the drug concentration. Here we study the effect of incorporating a Michaelis-Menten (MM) or $E_{\max}$-type pharmacodynamic model, the most commonly used model in the field of pharmacometrics. We compare the formulations of both problems and the resulting solutions. The reformulated problem with $E_{\max}$ pharmacodynamics is no longer linear in the control. This results in qualitative changes in the structure of optimal controls which, in line with an interpretation as concentrations, now are continuous while discontinuities exist if the log-kill model is used which is more in line with an interpretation of the control as dose rates. In spite of these qualitative differences, similarities in the structures of solutions can be observed. Both aspects are discussed theoretically and illustrated numerically.

Citation: Maciej Leszczyński, Urszula Ledzewicz, Heinz Schättler. Optimal control for a mathematical model for anti-angiogenic treatment with Michaelis-Menten pharmacodynamics. Discrete & Continuous Dynamical Systems - B, 2019, 24 (5) : 2315-2334. doi: 10.3934/dcdsb.2019097
##### References:

show all references

##### References:
$E_{\max}$ pharmacodynamic model.
Geometric shape of the sets ${\mathcal L}_0$ and ${\mathcal L}_m$ and definition of the associated regions $R_i$ and $S_i$ for ${i}{\rm{ = 1, 2, 3, 4}}{\rm{.}}$
The region $S_2$ (inside the loop ${\mathcal L}_m$ and above the diagonal) is not positively invariant for $u = u_{\max}$.
Graph of the value $V(\varepsilon)$ for the one-parameter family ${\mathcal F}_{\varepsilon}$ near its minimum value.
Optimal control (left) and corresponding trajectory (right) for the initial condition $(p_0, q_0) = (4600, 5800)$ and $A = 12$. The graphs are shown with $UC_{50} = 50$, $U_{\max} = 1.5UC_{50} = 75$ and $A = 12UC_{50} = 600$. The optimal controlled trajectory crosses the diagonal along a full dose segment and the optimal concatenation sequence is full dose $\rightarrow$ $\tilde{u}$ $\rightarrow$ no dose. The initial full dose segment is shown in black, the intermediate interior control segment in blue and the final no dose segment in black; junctions are marked with an asterisk. The loop ${\mathcal L}_m$ is shown as the dashed red curve.
Optimal control (left) and corresponding trajectory (right) for the initial condition $(p_0, q_0) = (12000, 15000)$ and $A = 2$. The graphs are shown with $UC_{50} = 50$, $U_{\max} = 1.5UC_{50} = 75$ and $A = 2UC_{50} = 100$. The optimal controlled trajectory crosses the diagonal using intermediate doses along the control $\tilde{u}$ and the optimal concatenation sequence is of the shortened form $\tilde{u}$ $\rightarrow$ no dose. The optimal controlled trajectory is shown in red and as a comparison the corresponding full dose trajectory is shown in blue. Here there is a significant difference in the final tumor volumes.
Parameter values used in numerical computations.
 variable/coefficients interpretation numerical value dimension $p$ tumor volume $mm^{3}$ $q$ carrying capacity of the vasculature $mm^{3}$ $\xi$ tumor growth parameter $0.192$ per day $\mu$ natural loss of endothelial support $0.02$ per day $b$ stimulation parameter, 'birth' $5.85$ per day $d$ inhibition parameter, 'death' $0.00873$ $mm^{-2}$ per day $G$ anti-angiogenic killing coefficient at limiting concentration $10$ per day $UC_{50}$ concentration with 50% effect $1$ $u_{\max}$ normalized maximum dose rate $1.5$ $A$ normalized total dose $2$; $12$
 variable/coefficients interpretation numerical value dimension $p$ tumor volume $mm^{3}$ $q$ carrying capacity of the vasculature $mm^{3}$ $\xi$ tumor growth parameter $0.192$ per day $\mu$ natural loss of endothelial support $0.02$ per day $b$ stimulation parameter, 'birth' $5.85$ per day $d$ inhibition parameter, 'death' $0.00873$ $mm^{-2}$ per day $G$ anti-angiogenic killing coefficient at limiting concentration $10$ per day $UC_{50}$ concentration with 50% effect $1$ $u_{\max}$ normalized maximum dose rate $1.5$ $A$ normalized total dose $2$; $12$
 [1] Urszula Ledzewicz, Helen Moore. Optimal control applied to a generalized Michaelis-Menten model of CML therapy. Discrete & Continuous Dynamical Systems - B, 2018, 23 (1) : 331-346. doi: 10.3934/dcdsb.2018022 [2] Nasser Sweilam, Fathalla Rihan, Seham AL-Mekhlafi. A fractional-order delay differential model with optimal control for cancer treatment based on synergy between anti-angiogenic and immune cell therapies. Discrete & Continuous Dynamical Systems - S, 2020, 13 (9) : 2403-2424. doi: 10.3934/dcdss.2020120 [3] Filippo Cacace, Valerio Cusimano, Alfredo Germani, Pasquale Palumbo, Federico Papa. Closed-loop control of tumor growth by means of anti-angiogenic administration. Mathematical Biosciences & Engineering, 2018, 15 (4) : 827-839. doi: 10.3934/mbe.2018037 [4] Karl Peter Hadeler. Michaelis-Menten kinetics, the operator-repressor system, and least squares approaches. Mathematical Biosciences & Engineering, 2013, 10 (5&6) : 1541-1560. doi: 10.3934/mbe.2013.10.1541 [5] Shuo Wang, Heinz Schättler. Optimal control for cancer chemotherapy under tumor heterogeneity with Michealis-Menten pharmacodynamics. Discrete & Continuous Dynamical Systems - B, 2019, 24 (5) : 2383-2405. doi: 10.3934/dcdsb.2019100 [6] Manuel Delgado, Cristian Morales-Rodrigo, Antonio Suárez. Anti-angiogenic therapy based on the binding to receptors. Discrete & Continuous Dynamical Systems - A, 2012, 32 (11) : 3871-3894. doi: 10.3934/dcds.2012.32.3871 [7] Jagadeesh R. Sonnad, Chetan T. Goudar. Solution of the Michaelis-Menten equation using the decomposition method. Mathematical Biosciences & Engineering, 2009, 6 (1) : 173-188. doi: 10.3934/mbe.2009.6.173 [8] Ming Liu, Dongpo Hu, Fanwei Meng. Stability and bifurcation analysis in a delay-induced predator-prey model with Michaelis-Menten type predator harvesting. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020259 [9] Ellina Grigorieva, Evgenii Khailov, Andrei Korobeinikov. An optimal control problem in HIV treatment. Conference Publications, 2013, 2013 (special) : 311-322. doi: 10.3934/proc.2013.2013.311 [10] Maria do Rosário de Pinho, Helmut Maurer, Hasnaa Zidani. Optimal control of normalized SIMR models with vaccination and treatment. Discrete & Continuous Dynamical Systems - B, 2018, 23 (1) : 79-99. doi: 10.3934/dcdsb.2018006 [11] Joaquim P. Mateus, Paulo Rebelo, Silvério Rosa, César M. Silva, Delfim F. M. Torres. Optimal control of non-autonomous SEIRS models with vaccination and treatment. Discrete & Continuous Dynamical Systems - S, 2018, 11 (6) : 1179-1199. doi: 10.3934/dcdss.2018067 [12] Cristiana J. Silva, Delfim F. M. Torres. A TB-HIV/AIDS coinfection model and optimal control treatment. Discrete & Continuous Dynamical Systems - A, 2015, 35 (9) : 4639-4663. doi: 10.3934/dcds.2015.35.4639 [13] Urszula Ledzewicz, Mohammad Naghnaeian, Heinz Schättler. Dynamics of tumor-immune interaction under treatment as an optimal control problem. Conference Publications, 2011, 2011 (Special) : 971-980. doi: 10.3934/proc.2011.2011.971 [14] Yali Yang, Sanyi Tang, Xiaohong Ren, Huiwen Zhao, Chenping Guo. Global stability and optimal control for a tuberculosis model with vaccination and treatment. Discrete & Continuous Dynamical Systems - B, 2016, 21 (3) : 1009-1022. doi: 10.3934/dcdsb.2016.21.1009 [15] Cristiana J. Silva, Delfim F. M. Torres. Optimal control strategies for tuberculosis treatment: A case study in Angola. Numerical Algebra, Control & Optimization, 2012, 2 (3) : 601-617. doi: 10.3934/naco.2012.2.601 [16] Sanjukta Hota, Folashade Agusto, Hem Raj Joshi, Suzanne Lenhart. Optimal control and stability analysis of an epidemic model with education campaign and treatment. Conference Publications, 2015, 2015 (special) : 621-634. doi: 10.3934/proc.2015.0621 [17] Djamila Moulay, M. A. Aziz-Alaoui, Hee-Dae Kwon. Optimal control of chikungunya disease: Larvae reduction, treatment and prevention. Mathematical Biosciences & Engineering, 2012, 9 (2) : 369-392. doi: 10.3934/mbe.2012.9.369 [18] Holly Gaff, Elsa Schaefer. Optimal control applied to vaccination and treatment strategies for various epidemiological models. Mathematical Biosciences & Engineering, 2009, 6 (3) : 469-492. doi: 10.3934/mbe.2009.6.469 [19] Kbenesh Blayneh, Yanzhao Cao, Hee-Dae Kwon. Optimal control of vector-borne diseases: Treatment and prevention. Discrete & Continuous Dynamical Systems - B, 2009, 11 (3) : 587-611. doi: 10.3934/dcdsb.2009.11.587 [20] M. Delgado-Téllez, Alberto Ibort. On the geometry and topology of singular optimal control problems and their solutions. Conference Publications, 2003, 2003 (Special) : 223-233. doi: 10.3934/proc.2003.2003.223

2019 Impact Factor: 1.27