# American Institute of Mathematical Sciences

• Previous Article
Newton's method for nonlinear stochastic wave equations driven by one-dimensional Brownian motion
• MBE Home
• This Issue
• Next Article
Local controllability and optimal control for\newline a model of combined anticancer therapy with control delays
February  2017, 14(1): 217-235. doi: 10.3934/mbe.2017014

## On drug resistance and metronomic chemotherapy: A mathematical modeling and optimal control approach

 1 Institute of Mathematics, Lodz University of Technology, 90-924 Lodz, Poland 2 Dept. of Mathematics and Statistics, Southern Illinois University, Edwardsville, Il, 62026-1653, USA 3 Dept. of Electrical and Systems Engineering, Washington University, St. Louis, Mo, 63130, USA 4 Service d'hématologie et Oncologie Pédiatrique, Centre Hospitalo-Universitaire Timone Enfants, AP-HM, Marseille, UMR S_911 CRO2 Aix Marseille Université, Marseille, France 5 Metronomics Global Health Initiative, Marseille, France 6 Childrens Cancer Institute Australia, Lowy Cancer Research Centre, University of New South Wales, Randwick, NSW, Australia

* Corresponding author

Received  December 2015 Accepted  June 30, 2016 Published  October 2016

Effects that tumor heterogeneity and drug resistance have on the structure of chemotherapy protocols are discussed from a mathematical modeling and optimal control point of view. In the case when two compartments consisting of sensitive and resistant cells are considered, optimal protocols consist of full dose chemotherapy as long as the relative proportion of sensitive cells is high. When resistant cells become more dominant, optimal controls switch to lower dose regimens defined by so-called singular controls. The role that singular controls play in the structure of optimal therapy protocols for cell populations with a large number of traits is explored in mathematical models.

Citation: Urszula Ledzewicz, Shuo Wang, Heinz Schättler, Nicolas André, Marie Amélie Heng, Eddy Pasquier. On drug resistance and metronomic chemotherapy: A mathematical modeling and optimal control approach. Mathematical Biosciences & Engineering, 2017, 14 (1) : 217-235. doi: 10.3934/mbe.2017014
##### References:

show all references

##### References:
Schematic illustration of the distribution of traits around dominant steady states
Dose rates (left) for the control $u_{\tau}$ defined by equation (26) for $\tau=1$ (top row), $\tau=1.75$ (middle row) and $\tau =2$ (bottom row), and corresponding time evolutions of the concentrations $c$ (middle) and the states $S$ and $R$ (right)
Evolution of the states S and R if the therapy horizon is extended to 42 days for the control strategies uτ for τ = 1:75 (left) and τ = 2 (right)
An administration protocol which gives full dose until time $\tau$ and then switches to the dose given by the singular control over the interval $[\tau, 10]$ optimized over $\tau$. The weights in the objective are chosen equal using $\bar{\alpha} = \bar{\beta} = \frac{1}{21}$ and $\gamma = 100$. Shown are the control (top, left), corresponding evolution of the total population $T=\sum_{i=1}^{21} N_i(t)$ (top, right), evolution of the traits $N_i(t)$ for $i=1, \ldots, 21$ (middle) and a comparison of the initial density $n(0, x)\equiv \frac{200}{21}$ and the terminal density $n(10, x)$ shown as red curve (bottom)
An administration protocol which gives full dose until time $\tau$ and then switches to the dose given by the singular control over the interval $[\tau, 10]$ optimized over $\tau$. The weights in the objective are chosen equal using $\bar{\alpha} = \frac{1}{21}$, $\bar{\beta} = \frac{8}{21}$ and $\gamma = 100$. Shown are the control (top, left), corresponding evolution of the total population $T=\sum_{i=1}^{21} N_i(t)$ (top, right), evolution of the traits $N_i(t)$ for $i=1, \ldots, 21$ (middle) and a comparison of the initial density $n(0, x)\equiv \frac{200}{21}$ and the terminal density $n(10, x)$ shown as red curve (bottom)
Values for the initial data and parameters used in numerical computations
 parameters interpretation value S0 initial condition of sensitive cells 9:4051 × 109 R0 initial condition of resistant cells 0:5949 × 109 r1 growth rate of sensitive population 3:5 r2 growth rate of resistant population 1 θ1 rate at which sensitive cells become resistant 0:15 θ2 rate at which resistant cells become resensitized 0:02 ψ1 log-kill coefficient for sensitive population 5 ψ2 log-kill coefficient for resistant population 1 k clearance rate of drug [Taxol] 2:9706
 parameters interpretation value S0 initial condition of sensitive cells 9:4051 × 109 R0 initial condition of resistant cells 0:5949 × 109 r1 growth rate of sensitive population 3:5 r2 growth rate of resistant population 1 θ1 rate at which sensitive cells become resistant 0:15 θ2 rate at which resistant cells become resensitized 0:02 ψ1 log-kill coefficient for sensitive population 5 ψ2 log-kill coefficient for resistant population 1 k clearance rate of drug [Taxol] 2:9706
 [1] Urszula Ledzewicz, Heinz Schättler. Drug resistance in cancer chemotherapy as an optimal control problem. Discrete & Continuous Dynamical Systems - B, 2006, 6 (1) : 129-150. doi: 10.3934/dcdsb.2006.6.129 [2] Shuo Wang, Heinz Schättler. Optimal control of a mathematical model for cancer chemotherapy under tumor heterogeneity. Mathematical Biosciences & Engineering, 2016, 13 (6) : 1223-1240. doi: 10.3934/mbe.2016040 [3] Urszula Ledzewicz, Heinz Schättler, Shuo Wang. On the role of tumor heterogeneity for optimal cancer chemotherapy. Networks & Heterogeneous Media, 2019, 14 (1) : 131-147. doi: 10.3934/nhm.2019007 [4] 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 [5] Urszula Ledzewicz, Heinz Schättler, Mostafa Reisi Gahrooi, Siamak Mahmoudian Dehkordi. On the MTD paradigm and optimal control for multi-drug cancer chemotherapy. Mathematical Biosciences & Engineering, 2013, 10 (3) : 803-819. doi: 10.3934/mbe.2013.10.803 [6] Arturo Alvarez-Arenas, Konstantin E. Starkov, Gabriel F. Calvo, Juan Belmonte-Beitia. Ultimate dynamics and optimal control of a multi-compartment model of tumor resistance to chemotherapy. Discrete & Continuous Dynamical Systems - B, 2019, 24 (5) : 2017-2038. doi: 10.3934/dcdsb.2019082 [7] Wei Feng, Shuhua Hu, Xin Lu. Optimal controls for a 3-compartment model for cancer chemotherapy with quadratic objective. Conference Publications, 2003, 2003 (Special) : 544-553. doi: 10.3934/proc.2003.2003.544 [8] Urszula Ledzewicz, Heinz Schättler. The Influence of PK/PD on the Structure of Optimal Controls in Cancer Chemotherapy Models. Mathematical Biosciences & Engineering, 2005, 2 (3) : 561-578. doi: 10.3934/mbe.2005.2.561 [9] Urszula Ledzewicz, Behrooz Amini, Heinz Schättler. Dynamics and control of a mathematical model for metronomic chemotherapy. Mathematical Biosciences & Engineering, 2015, 12 (6) : 1257-1275. doi: 10.3934/mbe.2015.12.1257 [10] Piotr Bajger, Mariusz Bodzioch, Urszula Foryś. Singularity of controls in a simple model of acquired chemotherapy resistance. Discrete & Continuous Dynamical Systems - B, 2019, 24 (5) : 2039-2052. doi: 10.3934/dcdsb.2019083 [11] Luis A. Fernández, Cecilia Pola. Catalog of the optimal controls in cancer chemotherapy for the Gompertz model depending on PK/PD and the integral constraint. Discrete & Continuous Dynamical Systems - B, 2014, 19 (6) : 1563-1588. doi: 10.3934/dcdsb.2014.19.1563 [12] Urszula Ledzewicz, Mozhdeh Sadat Faraji Mosalman, Heinz Schättler. Optimal controls for a mathematical model of tumor-immune interactions under targeted chemotherapy with immune boost. Discrete & Continuous Dynamical Systems - B, 2013, 18 (4) : 1031-1051. doi: 10.3934/dcdsb.2013.18.1031 [13] Craig Collins, K. Renee Fister, Bethany Key, Mary Williams. Blasting neuroblastoma using optimal control of chemotherapy. Mathematical Biosciences & Engineering, 2009, 6 (3) : 451-467. doi: 10.3934/mbe.2009.6.451 [14] Cristian Tomasetti, Doron Levy. An elementary approach to modeling drug resistance in cancer. Mathematical Biosciences & Engineering, 2010, 7 (4) : 905-918. doi: 10.3934/mbe.2010.7.905 [15] Joseph Malinzi, Rachid Ouifki, Amina Eladdadi, Delfim F. M. Torres, K. A. Jane White. Enhancement of chemotherapy using oncolytic virotherapy: Mathematical and optimal control analysis. Mathematical Biosciences & Engineering, 2018, 15 (6) : 1435-1463. doi: 10.3934/mbe.2018066 [16] Urszula Ledzewicz, Heinz Schättler. Controlling a model for bone marrow dynamics in cancer chemotherapy. Mathematical Biosciences & Engineering, 2004, 1 (1) : 95-110. doi: 10.3934/mbe.2004.1.95 [17] Urszula Ledzewicz, Helmut Maurer, Heinz Schättler. Optimal and suboptimal protocols for a mathematical model for tumor anti-angiogenesis in combination with chemotherapy. Mathematical Biosciences & Engineering, 2011, 8 (2) : 307-323. doi: 10.3934/mbe.2011.8.307 [18] Urszula Ledzewicz, Omeiza Olumoye, Heinz Schättler. On optimal chemotherapy with a strongly targeted agent for a model of tumor-immune system interactions with generalized logistic growth. Mathematical Biosciences & Engineering, 2013, 10 (3) : 787-802. doi: 10.3934/mbe.2013.10.787 [19] Ami B. Shah, Katarzyna A. Rejniak, Jana L. Gevertz. Limiting the development of anti-cancer drug resistance in a spatial model of micrometastases. Mathematical Biosciences & Engineering, 2016, 13 (6) : 1185-1206. doi: 10.3934/mbe.2016038 [20] Clara Rojas, Juan Belmonte-Beitia, Víctor M. Pérez-García, Helmut Maurer. Dynamics and optimal control of chemotherapy for low grade gliomas: Insights from a mathematical model. Discrete & Continuous Dynamical Systems - B, 2016, 21 (6) : 1895-1915. doi: 10.3934/dcdsb.2016028

2018 Impact Factor: 1.313