Difference between revisions of "De Pillis chemotherapy model"

From mintOC
Jump to: navigation, search
(→‎References)
(→‎Source Code)
Line 77: Line 77:
 
==Source Code==
 
==Source Code==
  
−
[[:Category:C | C code]]  
+
* [[:Category:C | C code]] at [[De Pillis chemotherapy model (C)]]
 
<source lang="cpp">
 
<source lang="cpp">
 
/* volume of tumor in absolute cell count*/
 
/* volume of tumor in absolute cell count*/

Revision as of 18:40, 27 January 2016

De Pillis chemotherapy model
State dimension: 1
Differential states: 6
Continuous control functions: 3


The model by de Pillis combines chemotherapy with immunotherapy.

Mathematical formulation

For t \in [t_0, t_f] and D = d \frac{(x_2/x_0)^l}{s+(x_2/x_0)^l} the optimal control problem is given by

\begin{array}{llcl} \displaystyle \min_{x, u} & p_0 x_0(t_f) & + & \int_{t_0}^{t_f} p_1 x_0(t)^2 \mbox{d}t + \sum_{i=0}^{3} \int_{t_0}^{t_f} p_{i+2} u_i(t)^2 \mbox{d}t \\[1.5ex] \mbox{s.t.} & \dot{x}_0(t) & = & a x_0 (1-b x_0) -c x_1 x_0 - D x_0 - K_T (1- \mbox{e}^{- x_4}) x_0 ,\\ & \dot{x}_1(t) & = & e x_3 - f x_1 + g \frac{x_0^2}{h+x_0^2}-p x_1 x_0 - K_N (1- \mbox{e}^{- x_4}) x_1, \\ & \dot{x}_2(t) & = & -m x_2 + j \frac{D^2 x_0^2}{k+ D^2 x_0^2} x_2 - q x_1 x_2 + (r_1 x_1 + r_2 x_3) x_0 \\ & & & - v x_1 x_2^2 - K_L (1- \mbox{e}^{- x_4}) x_2 + \frac{p_I x_2 x_5}{g_I + x_5} + u_2, \\ & \dot{x}_3(t) & = & \alpha - \beta x_3 - K_C (1- \mbox{e}^{- x_4}) x_3,\\ & \dot{x}_4(t) & = & - \gamma x_4 + u_0,\\ & \dot{x}_5(t) & = & - \mu_I x_5 + u_1. \end{array}

The states x_0 to x_3 are measured in absolute cell counts, where x_0 describes the number of tumor cells, x_1 of unspecific immune cells, x_2 of tumor-specific cytotoxic T-cells (CD8^+ T) and x_3 of circulating lymphocytes. The chemotherapeutic drug concentration is given by x_4 and the immunotherapeutic by x_5 (Interleukin-2) respectively.

Parameters

This set of parameters can be found as “patient 9” in [].

\begin{array}{lll} a = 4.31 \times 10^{-1} & b = 1.02 \times 10^{-9} & c = 6.41 \times 10^{-11}\\ d = 2.34 & e = 2.08 \times 10^{-7} & f = 4.12 \times 10^{-2}\\ g = 1.25 \times 10^{-2} & h = 2.02 \times 10^{7} & j = 2.49 \times 10^{-2}\\ k = 3.66 \times 10^{7} & l = 2.09 & m = 2.04 \times 10^{-1}\\ q = 1.42 \times 10^{-6} & p = 3.42 \times 10^{-6} & s = 8.39 \times 10^{-2}\\ r_1 = 1.01 \times 10^{-7} & r_2 = 6.50 \times 10^{-11} & u = 3.00 \times 10^{-10}\\ \alpha = 7.50 \times 10^{8} & \beta = 1.20 \times 10^{-2} & \gamma = 9.00 \times 10^{-1}\\ p_I = 1.25 \times 10^{-1} & g_I = 2.00 \times 10^{7} & \mu_I = 1.00 \times 10^{1}\\ K_T = 9.00 \times 10^{-1} & K_N = K_L = K_C = 6 \times 10^{-1} \end{array}

Reference Solutions

The problem can be solved with the [multiple shooting method]. For the following solutions the control functions and states are discretized on the same grid, with 100 nodes. For the objective function parameters have been chosen from the following sets to grand the results shown in the graphics below.

Objective function 1 \begin{array}{rclr} p_0 &=& 1,&\\ p_2 &=& 1,&\\ p_i &=& 0,&otherwise.\\ \end{array}

Objective function 2 \begin{array}{rclr} p_0 &=& -1,&\\ p_2 &=& 1,&\\ p_i &=& 0,&otherwise.\\ \end{array}

In both objective function the amount of chemotherapeutic drugs is penalized. The objective function 2 describes the worst case scenario of the tumor growth at end time.

  • Reference solution plots

  • Optimal controls and states computed with a multiple shooting approach on the same discretization grid with 100 nodes and objective function 1.

  • Optimal controls and states with objective function 2.

Source Code

/* volume of tumor in absolute cell count*/
rhs[0] = - a*x[0] (1-b*x[0]) -c*x[1]*x[0] - D*x[0] - K_T*(1- exp(- x[4]))*x[0];
 
/* volume of unspecific immune cells */
rhs[1] = e*x[3] - f*x[1] + g*x[0]*x[0]/(h+x[0]*x[0])-p*x[1]*x[0] - K_N*(1- exp(- x[4]))*x[1];
 
/* volume of tumor-specific cytotoxic T-cells */
rhs[2] = -m*x[2] + j*D*D*x[0]*x[0]/(k+ D*D*x[0]*x[0])*x[2] - q*x[1]*x[2] + (r_1*x[1] + r_2*x[3])*x[0] ...
- v*x[1]* x[2]*x[2] - K_L* (1- exp(- x[4]))*x[2] + p_I*x[2]*x[5]/(g_I + x[5]) + u[2];
 
/* volume of circulating lymphocytes */
rhs[3] = alpha - beta *x[3] - K_C *(1- exp(- x[4])) *x[3];
 
/* amount of chemotherapeutic drug */
rhs[4] = u[0] - gamma*x[4];
 
/* amount of immunotherapeutic drug */
rhs[5] = u[1] - mu_I*x[5];

References

Retrieved from "https://mintoc.de/index.php?title=De_Pillis_chemotherapy_model&oldid=1514"