# Optimal Treatment Strategy for a Tumor Model under Immune Suppression

Authors: Kwang Su Kim; Giphil Cho; Il Hyo Jung
Journal: Computational and Mathematical Methods in Medicine
 (2014)

---

## Abstract

We propose a mathematical model describing tumor-immune interactions under immune suppression. These days evidences indicate that the immune suppression related to cancer contributes to its progression. The mathematical model for tumor-immune interactions would provide a new methodology for more sophisticated treatment options of cancer. To do this we have developed a system of 11 ordinary differential equations including the movement, interaction, and activation of NK cells, CD8+T-cells, CD4+T cells, regulatory T cells, and dendritic cells under the presence of tumor and cytokines and the immune interactions. In addition, we apply two control therapies, immunotherapy and chemotherapy to the model in order to control growth of tumor. Using optimal control theory and numerical simulations, we obtain appropriate treatment strategies according to the ratio of the cost for two therapies, which suggest an optimal timing of each administration for the two types of models, without and with immunosuppressive effects. These results mean that the immune suppression can have an influence on treatment strategies for cancer.


---

## Body

## 1. Introduction

Cancer is a leading cause of death worldwide. Cancer is a term used for diseases in which abnormal cells divide without control and are able to invade some tissues. The possibility that cancers can be eradicated by specific immune responses has been the impetus for the field of tumor immunology. The existence of immune surveillance, which was proposed by Macfarlane Burnet in the 1950s, has been demonstrated by the incidence of some types of tumors in immunocompromised experimental animals and humans [1].Chemotherapy and immunotherapy are typical treatment methods of cancer. Chemotherapy directly targets the transformed tumor cell. But chemotherapy has some general side effects such as hair loss, a sore mouth, vomiting, and diarrhea [2]. Immunotherapy is treatment method that uses body’s own immune system to help fighting cancer. There are such many treatment methods for cancer. It is important to know how to combine these treatment methods as well as to find treatment methods for cancer treatment. Although the immune system is very complex, we may suggest more effective treatment strategies for cancer control through mathematical models associated with the immune system. So, identifying a mathematical model of tumor-immune interactions that mediate the immune responses and immune suppression would provide a new strategy for more sophisticated treatment methods.There are amount of papers that deal with mathematical models about tumor-immune interaction. In Kuznetsov et al. [3], a mathematical model describing the cytotoxic T lymphocyte response to the growth of an immunogenic tumor was proposed. Kirschner and Panetta [4] presented a mathematical model describing the dynamics between tumor cells, immune-effector cells, and IL-2. They explained both short-tumor oscillations in tumor sizes and long-term tumor relapse. In de Pillis et al. [5], they proposed a mathematical model which was based on de Pillis and Radunskaya’s model [6] and includes tumor cells and three immune cells as well as two drug concentrations in the bloodstream and the model described bifurcation-like behavior to be reproduced under particular realistic conditions. In addition, they showed that combination therapy has more effects than only one therapy for tumor control through numerical simulations. The model in [7] extended de Pillis’ model [6] and identified appropriate values for the parameters of the new model according to recent empirical data. de Pillis’ paper [8] dealt with optimal control problem for the model of tumor-immune interactions with chemotherapy. They characterized the optimal controls related to drug therapy, including a quadratic control, a linear control, and state constraint. Also, Engelhart et al. [9] introduced four different cancer mathematical models of chemotherapy from the literature and compared with results of optimal control on their models. Here we note that the above-mentioned all models describe tumor-immune interactions without immune suppression. Although the immune system response tumor cells, the fact that each year many people die from cancer suggests that immune response to tumor cells is often ineffective. In fact, an immunosuppressive effect is known to be one of the main causes of such phenomenon so that there are several mechanisms by which tumor cells appear to evade the immune system [1, 10]. The paper of Vaage [11] gave that the immune system plays a significant role in removing the tumors of cancer and pointed out that it does not always work. Chouaib et al. [12] found out that immune suppression obstructs the immune system in removing tumors. Recently, Robertson-Tessi et al. [13] proposed a mathematical model of tumor-immune interactions. They showed only existence of optimal antigenicity maximizing the immune response.So even though considering immunosuppressive effects in the tumor model is important, mathematical modeling for tumor-immune interactions under immune suppression is still rare. Here we focus on finding optimal treatment strategies for tumor-immune models under immune suppression.For this goal, combining the ideas of [7, 13], we consider an improved tumor dynamic model including three types of immune response: NK cell as innate immune response; CD8+T cell, CD4+T cell, and IL-2 as adaptive immune response; and regulatory T cell, TGF-β, and IL-10 as a role of immune suppression. Actually we have developed a system of 11 ordinary differential equations describing the movement, interaction, and activation of tumor-immune system under immune suppression. We set up a controlled tumor immune system incorporating two control measures: immunotherapy and chemotherapy. In addition, we establish the existence of an optimal control for the control system and also drive an optimal control and the optimality system using optimal control techniques. Using the numerical simulations, we suggest the optimal treatment strategies on the base of time for the control system under immunosuppressive effects. We find out that if the cost of chemotherapy is more expensive than that of immunotherapy, then the optimal treatment strategy for the model without immunosuppressive effects needs to be taken for a longer time for chemotherapy comparing immunotherapy. On the other hand, when the model has immunosuppressive effects and the cost of chemotherapy is more expensive than that of immunotherapy, the optimal treatment strategy for the model needs to be taken for a longer time for immunotherapy comparing chemotherapy.The rest of this paper is organized as follows. In Section2, we give the mathematical model for the tumor-immune system and describe the parameters used in the equations with optimal control problem. The necessary conditions for an optimal control and the corresponding states are derived using Pontryagin’s maximum principle. In Section 3, we show that the efficacy of immunotherapy and chemotherapy using numerical simulation and the resulting optimality system and the parameter sensitivity analysis are numerically solved. Finally, discussion and conclusions are given in Section 4.


## 2. Mathematical Modeling and Optimal Control Problem for a Tumor-Immune System

If a tumor cell formed in our body, firstly by innate immune response, NK cells kill the tumor cells. By the adaptive immune response, after dendritic cells knew about tumor cell’s information, dendritic cells convey information about tumor to other T cells. CD4+T cells and CD8+T cells activated each other by cytokine IL-2. By influence of immune suppression, regulatory T cells and tumor cell limit functions of other T cells by cytokine IL-10 and TGF-β. Next, CD8+T cells directly kill tumor cells by adaptive immune response [1] (see Figure 1).Figure 1
Tumor and immune response to chemotherapy and immunotherapy.In this section, we describe the above-mentioned three important immune responses and they are represented in Figure1. Our baseline model [13] describes adaptive immune response and influence of immune suppression, but they do not include effects of the innate immune response. Although de Pillis’ model [7] includes the adaptive immune response and the innate immune response, they did not consider the immune suppression to their model. Motivated by these two papers, we make mathematical model describing tumor-immune interactions which include three immune responses, chemotherapy, and immunotherapy by using mathematical model.This model is based on Robertson-Tessi’s model [13]. Their model used 12 biological variables: nine cell types and three cytokines. All three types of T cells are specific to the tumor antigen and each of the three T-cell populations is broken into three subpopulations. The memory T cell (M
E, M
H, and M
R) is a pool of T cell precursors, which are activated by dendritic cells. There is a short-lived activation phase (A
E, A
H, and A
R) where T cell proliferation occurs. The fully functional T cell phase (E, H, R) consists of the cells that kill the tumor cells, produce cytokines, and suppress the immune response. The dendritic cell population has two subpopulations. The unlicensed state, U
(
t
), is a mature dendritic cell that collects antigen from the tumor and then interacts with helper cells. After licensing, the dendritic cell, D
(
t
), is free to interact with all T cells, causing them to activate. For update to our baseline model, firstly, we assume logistic model for growth of tumor. Since experimental data of tumor growth for the kinds of tumor is different, there are many growth curves to describe tumor growth. So we use more a simple logistic growth curve than the tumor growth curve in Robertson-Tessi’s model [13]. Also, by substitute T for T
* in [13], we exclude effects of spatial access to the tumor cells. In Robertson-Tessi’s model [13], all three types of T cells are activated in the same basic way. The memory T cells enter a brief activation phase where proliferation is rapid, and then they become fully functional T cells. Thus, we may omit states of A
E, A
H, and A
R because the stages of A
E, A
H, and A
R only activate stages of E, H, and R, respectively.Secondly, we add-
c
3
E
T
T for describing CD8+T cell death by exhaustion of tumor-killing resources in the third equation of the system (1) [7]. Thirdly, we include the state of NK cells and describe innate immune response [7]. Finally, for describing death of cells due to chemotherapy, we include state of chemotherapy concentration and chemotherapy-induced cell death term [5]. The completed model is described by the following system of differential equations:
(1)
d
T
d
t
=
a
T
(
1
-
b
T
)
-
c
1
T
N
-
d
T
E
T
E
T
+
e
T
1
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
-
k
T
(
1
-
e
-
M
)
T
,
d
N
d
t
=
b
1
-
d
N
N
-
c
2
T
N
+
p
N
N
I
2
q
N
+
I
2
-
k
N
(
1
-
e
-
M
)
N
,
d
E
T
d
t
=
α
1
I
2
D
L
m
E
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
-
c
3
E
T
T
-
k
E
T
(
1
-
e
-
M
)
E
T
+
v
E
+
ω
1
u
E
(
t
)
,
d
D
U
d
t
=
p
T
(
1
+
(
I
10
/
i
2
)
)
(
1
+
(
G
T
/
g
2
)
)
-
γ
1
D
U
1
+
(
D
U
/
m
H
)
-
d
D
U
D
U
-
k
D
U
(
1
-
e
-
M
)
D
U
,
d
D
L
d
t
=
γ
1
D
U
1
+
(
D
U
/
u
1
)
-
d
D
L
D
L
,
d
H
T
d
t
=
α
2
I
2
(
D
L
+
D
U
)
m
H
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
-
γ
2
H
T
S
s
3
+
S
-
d
H
T
H
T
-
k
H
T
(
1
-
e
-
M
)
H
T
,
d
G
T
d
t
=
γ
2
H
T
S
s
3
+
S
+
α
3
I
2
D
L
m
G
(
i
1
+
I
2
)
(
d
1
+
D
L
)
-
d
G
T
G
T
-
k
G
T
(
1
-
e
-
M
)
G
T
,
d
S
d
t
=
p
1
G
T
+
p
2
T
-
d
s
S
,
d
I
2
d
t
=
α
4
H
T
(
1
+
(
I
10
/
i
3
)
)
(
1
+
(
S
/
s
4
)
)
-
d
I
2
I
2
,
d
I
10
d
t
=
p
3
G
T
+
p
4
T
-
d
I
10
I
10
,
d
M
d
t
=
-
d
M
M
+
v
M
+
ω
2
u
M
(
t
)
.HereT
(
t
) represents number of tumor cells. And immune system is partitioned into the innate immune response (N
(
t
) NK cell numbers), adaptive immune response (E
T
(
t
) (CD8+T cell numbers), D
U
(
t
) (unlicensed dendritic cell numbers), D
L (licensed dendritic cell numbers), H
T
(
t
) (CD4+T cell numbers), I
2
(
t
) (IL-2)), and immune suppression (G
T
(
t
) (regulatory T cell numbers), S
(
t
) (TGF-β), I
10
(
t
) (IL-10)).When we inject chemotherapy, tumor cells are killed but other cells are damaged also. Among the many methods of immunotherapy, we use method injecting CD8+T cells directly. M
(
t
) is chemotherapy drug concentration in the bloodstream. The parameters v
E
(
t
) and v
M
(
t
) represent amounts of chemotherapy and immunotherapy, respectively [7].For the immunochemotherapy, we follow amount of dosage suggested by the manufacturers of the drug Adria. We use the upper end of the dosing range to arrive atv
M
=
2.3869mg/L per every 21 days for chemotherapy and the amount of injected CD8+T cells number is v
L
=
1.77
×
1
0
10. Since the units of parameters are different in two models [7, 13], we have converted ng/mL into IU/L by using the specific activity of IL-10 (3
×
1
0
6 IU/mg) and TGF-β (2
×
1
0
7 IU/mg). i
1
=
0.3 ng/mL is taken from [13]. We arrive at our value for i
1 by employing the specific activity of IL-2 ((
1.8
×
1
0
7
)
/
1.1 IU/mg) to convert molar concentration to IU/L [7]. So we have
(2)
0.3
ng
1
mL
×
1
mg
1
0
6
ng
×
1000
mL
1
L
×
1.8
×
1
0
7
IU
1.1
mg
=
4909
IU/L
.
We summarized the model’s term descriptions and value of parameters in Table 1.Table 1
Equation terms and parameter descriptions.
State
Term
Parameter
Description
T
(
t
)
a
T
(
1
-
b
T
)
a
=
0.431 (Day
-
1), b
=
1.02
×
1
0
-
9 [5]
Tumor growth term
c
1
T
N
c
1
=
3.177
×
1
0
-
13 (Cells
-
1
Day
-
1) [7]
NK cells-induced tumor death
d
T
E
E
T
+
e
T
d
=
0.9 (Day
-
1) e
=
1.2 [13]
CD8+T cells-induced tumor death
1
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
g
1
=
1.2 [13]
Suppression of CD8+T cell activity  due to TGF-β and regulatory T cell
k
T
(
1
-
e
-
M
)
T
k
T
=
0.9 (Day
-
1) [7]
Chemotherapy-induced tumor death
N
(
t
)
b
1
b
1
=
3121875 (Day
-
1) [7]
Production of NK cell
d
N
N
d
N
=
0.0125 (Day
-
1) [7]
Turnover of NK cell
c
2
T
N
c
2= 3.177
×
1
0
-
13 (Cells
-
1
Day
-
1) [7]
NK death by exhaustion of tumor killing resources
p
N
N
I
2
q
N
+
I
2
p
N
=
0.0668 (Day
-
1), q
N
=
250360 (IU/L) [7]
Stimulatory effect of IL-2 on NK cell
k
N
(
1
-
e
-
M
)
N
k
N
=
0.6 (Day
-
1) [5]
Chemotherapy-induced NK cell death
E
T
(
t
)
α
1
I
2
D
L
m
E
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
α
1
=
16 (Day
-
1), s
2
=
580000 (IU/L) i
1
=
4909 (IU/L), d
1
=
579579, m
E
=
526800 [13]
Proliferation of CD8+T cells
c
3
E
T
T
c
3
=
3.42
×
1
0
-
10 (Cells
-
1
Day
-
1) [7]
CD8+T cell death by exhaustion of tumor killing resources
k
E
T
(
1
-
e
-
M
)
E
T
k
E
T
=
0.6 (Day
-
1) [5]
Chemotherapy-induced tumor death
D
U
(
t
)
p
T
(
1
+
(
I
10
/
i
2
)
)
(
1
+
(
G
T
/
g
2
)
)
p
=
0.1 (Day
-
1), i
2
=
1200 (IU/L) g
2
=
2
×
1
0
7 (cell) [13]
Proliferation of mature unlicensed dendritic cell
γ
1
D
U
1
+
(
D
U
/
m
H
)
γ
1
=
0.5 (Day
-
1), m
H
=
1053600 [13]
Licensing of dendritic cell upon encounter with CD4+T cell
d
D
U
D
U
d
D
U
=
0.14 (Day
-
1) [13]
Turnover of CD8+T cell
k
D
U
(
1
-
e
-
M
)
D
U
k
D
U
=
0.05 (Day
-
1) [5]
Chemotherapy-induced mature unlicensed dendritic cell death
D
L
(
t
)
d
D
L
D
L
d
D
L
=
0.5 (Day
-
1) [13]
Turnover of CD8+T cell
H
T
(
t
)
α
2
I
2
(
D
L
+
D
U
)
m
H
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
α
2
=
1.9 (Day
-
1) [13]
Proliferation of CD4+T cell
γ
2
H
T
S
s
3
+
S
γ
2
=
0.022 (Day
-
1), s
3
=
34000 (IU/L) [13]
Converting of CD4+T cell to regulatory T cell by TGF-β
d
H
T
H
T
d
H
T
=
0.1 (Day
-
1) [13]
Turnover of CD4+T cell
k
H
T
(
1
-
e
-
M
)
H
T
k
H
T
=
0.6 [5]
Chemotherapy-induced CD4+T cell death
G
T
(
t
)
α
3
I
2
D
L
m
G
(
i
1
+
I
2
)
(
d
1
+
D
L
)
α
3
=
3.6 (Day
-
1), m
G
=
175900 [13]
Proliferation of regulatory T cell
d
G
T
G
T
d
G
T
=
0.1 (Day
-
1) [13]
Turnover of regulatory T cell
k
G
T
(
1
-
e
-
M
)
G
T
k
G
T
=
0.6 [5]
Chemotherapy-induced regulatory T cell death
S
(
t
)
p
1
G
T
p
1
=
3.6
×
1
0
-
4 (IU/LCells
-
1
Day
-
1) [13]
Production of TGF-β by regulatory T cell
p
2
T
p
2
=
2.2
×
1
0
-
3 (IU/LCells
-
1
Day
-
1) [13]
Production of TGF-β by tumor cell
d
S
S
d
S
=
14.3 (Day
-
1) [13]
Turnover of TGF-β
I
2
(
t
)
α
4
H
T
(
1
+
(
I
10
/
i
3
)
)
(
1
+
(
S
/
s
4
)
)
α
4
=
0.278 (IU/LCells
-
1
Day
-
1), i
3
=
2250 (IU/L), s
4
=18000 (IU/L) [13]
Production of IL-2
d
I
2
I
2
d
I
2
=
12.5 (Day
-
1) [13]
Turnover of IL-2
I
10
(
t
)
p
3
G
T
p
3
=
4.2
×
1
0
-
5 (IU/LCells
-
1
Day
-
1) [13]
Production of IL-10 by regulatory T cell
p
4
T
p
4
=
3.9
×
1
0
-
7 (IU/LCells
-
1
Day
-
1) [13]
Production of IL-10 by tumor cell
d
I
10
I
10
d
I
10
=
20 (Day
-
1) [13]
Turnover of IL-10
M
(
t
)
d
M
M
d
M
=
0.9 (Day
-
1) [5]
Excretion and elimination of medicine toxicityIn this work, we introduce two control functionsu
E
(
t
) and u
M
(
t
) representing amounts of immunotherapy and chemotherapy, respectively. The parameters ω
1 and ω
2 are weight factors. We assume that number of injected CD+8 T cells for immunotherapy is 175000000 during 50 days. So we put the values of ω
1 is 3500000. And we assume that amounts of chemotherapy are 25 during 50 days. So we put the values of ω
2 to be 0.5.Next we consider a cost functional as follows:(3)
F
(
u
E
,
u
M
)
=
∫
0
t
f
‍
A
T
(
t
)
+
B
u
E
(
t
)
2
+
C
u
M
(
t
)
2
d
t
.
This functional includes the number of tumor cells, amounts of chemotherapy, and amounts of immunotherapy. In words, we are minimizing the number of tumor cells, amounts of chemotherapy, and amounts of immunotherapy. In the objective cost functional, the quantities A, B, and C represent the weight constants of tumor cell numbers, for immunotherapy and chemotherapy, respectively. The costs associated with immunotherapy and chemotherapy are described in the terms B
u
E
(
t
) and C
u
M
(
t
), respectively. Our goal is to minimize the cost functional (3), which is called the optimal control problem. That is, the optimal control problem is to seek optimal control functions (
u
E
*
(
t
)
,
u
M
*
(
t
)
) such that
(4)
F
(
u
E
*
,
u
M
*
)
=
min
⁡
{
F
(
u
E
,
u
M
)
,
(
u
E
,
u
M
)
∈
U
}
is subject to the system (1) and appropriate initial conditions are given at t
=
0, where the control set is defined as
(5)
U
=
{
u
=
(
u
E
,
u
M
)
∣
u
i
(
t
)
is
Lebesgue
measurable
,
m
m
0
≤
u
i
(
t
)
≤
1
,
t
∈
[
0
,
T
]
for
i
=
E
,
M
}
.Pontryagin’s maximum principle is used to solve this optimal control problem and the derivation of the necessary conditions. First we prove the existence of an optimal control for problem (4) and then derive the optimality system.Theorem 1.
Given that the cost functionalF
(
u
E
,
u
M
)
=
∫
0
t
f
‍
A
T
(
t
)
+
B
u
E
(
t
)
2
+
C
u
M
(
t
)
2
d
t and the control set U given by (5) is measurable, there exists on optimal control u
*
=
(
u
E
*
,
u
M
*
) such that F
(
u
E
*
,
u
M
*
)
=
min
⁡
{
F
(
u
E
,
u
M
)
,
(
u
E
,
u
M
)
∈
U
}.Proof.
To prove the existence of an optimal control, we use the result in [14]. Note that the control and the state variable are nonnegative values. In this minimizing problem, the necessary convexity of the objective functional in u
E, u
M are satisfied. The set of all the control variables (
u
E
,
u
M
)
∈
U is also convex and closed by definition. The optimal system is bounded which determines the compactness needed for the existence of the optimal control. In addition, the integrand in functional (3), A
T
(
t
)
+
B
u
E
(
t
)
2
+
C
u
M
(
t
)
2, is convex on the control set U. Also we can easily see that there exist a constant δ
>
1 and numbers ϕ
1, ϕ
2 such that
(6)
F
(
u
E
,
u
M
)
≥
ϕ
1
(
u
E
2
+
u
M
2
)
δ
/
2
-
ϕ
2
,
because the state variables are bounded, which completes the existence of an optimal control.In order to find an optimal solution of the system, first we should find the Lagrangian and Hamiltonian for the optimal control problem (4). The Lagrangian of the control problem is given by
(7)
L
=
A
T
(
t
)
+
B
u
E
(
t
)
2
+
C
u
M
(
t
)
2
.We seek for the minimal value of the Lagrangian. To do this, we define the Hamiltonian functionH for the system, where λ
i, i
=
1
,
…
,
11, are the adjoint variables:
(8)
H
=
A
T
(
t
)
+
B
u
E
(
t
)
2
+
C
u
M
(
t
)
2
+
λ
1
[
d
T
E
T
E
T
+
e
T
1
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
a
T
(
1
-
b
T
)
-
c
1
T
N
m
m
m
m
c
-
d
T
E
T
E
T
+
e
T
1
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
m
m
m
m
c
-
k
T
(
1
-
e
-
M
)
T
d
T
E
T
E
T
+
e
T
1
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
]
+
λ
2
[
b
1
-
d
N
N
-
c
2
T
N
+
p
N
N
I
2
q
N
+
I
2
-
k
N
(
1
-
e
-
M
)
N
]
+
λ
3
[
α
1
I
2
D
L
m
E
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
-
c
3
E
T
T
m
m
m
m
c
-
k
E
T
(
1
-
e
-
M
)
E
T
+
v
E
+
ω
1
u
E
(
t
)
α
1
I
2
D
L
m
E
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
]
+
λ
4
[
p
T
(
1
+
(
I
10
/
i
2
)
)
(
1
+
(
G
T
/
g
2
)
)
-
γ
1
D
U
1
+
(
D
U
/
u
1
)
m
m
m
m
c
-
d
D
U
D
U
-
k
D
U
(
1
-
e
-
M
)
D
U
p
T
(
1
+
(
I
10
/
i
2
)
)
(
1
+
(
G
T
/
g
2
)
)
]
+
λ
5
[
γ
1
D
U
1
+
(
D
U
/
u
1
)
-
d
D
L
D
L
]
+
λ
6
[
α
2
I
2
(
D
L
+
D
U
)
m
H
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
-
γ
2
H
T
S
s
3
+
S
m
m
m
m
c
-
d
H
T
H
T
-
k
H
T
(
1
-
e
-
M
)
H
T
α
2
I
2
(
D
L
+
D
U
)
m
H
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
]
+
λ
7
[
γ
2
H
T
S
s
3
+
S
+
α
3
I
2
D
L
m
G
(
i
1
+
I
2
)
(
d
1
+
D
L
)
m
m
m
m
c
-
d
G
T
G
T
-
k
G
T
(
1
-
e
-
M
)
G
T
γ
2
H
T
S
s
3
+
S
+
α
3
I
2
D
L
m
G
(
i
1
+
I
2
)
(
d
1
+
D
L
)
]
+
λ
8
[
p
1
G
T
+
p
2
T
-
d
s
S
]
+
λ
9
[
α
4
H
T
(
1
+
(
I
10
/
i
3
)
)
(
1
+
(
S
/
s
4
)
)
-
d
I
2
I
2
]
+
λ
10
[
p
3
G
T
+
p
4
T
-
d
I
10
I
10
]
+
λ
11
[
-
d
M
M
+
v
M
+
ω
2
u
M
(
t
)
]
.In order to derive the necessary conditions, we use Pontryagain’s maximum principle [15] as follows.If(
x
,
u
) is an optimal solution of an optimal control problem, then there exists a nontrivial vector function λ
=
(
λ
1
,
λ
2
,
…
,
λ
n
) satisfying the following inequalities:
(9)
d
x
d
t
=
∂
H
(
t
,
x
,
u
,
λ
)
∂
λ
,
0
=
∂
H
(
t
,
x
,
u
,
λ
)
∂
u
,
d
λ
d
t
=
-
∂
H
(
t
,
x
,
u
,
λ
)
∂
x
.We now drive the necessary conditions that optimal control functions and corresponding states must satisfy. In the following theorem, we present the adjoint system and control characterization.Theorem 2.
Given an optimal controlu
*
=
(
u
E
*
,
u
M
*
) and a solution y
*
=
(
T
*
,
N
*
,
E
T
*
,
D
U
*
,
D
L
*
,
H
T
*,G
T
*
,
S
*
,
I
2
*
,
I
10
*
,
M
*
) of the corresponding state system (1), there exists adjoint variables λ
i, i
=
1
,
…
,
11, satisfying
(10)
λ
1
′
(
t
)
=
-
∂
H
∂
T
=
-
A
-
λ
1
[
d
E
T
2
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
(
E
T
+
e
T
)
2
a
-
2
a
b
T
-
c
1
N
m
m
m
m
c
-
d
E
T
2
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
(
E
T
+
e
T
)
2
m
m
m
m
c
-
k
T
(
1
-
e
-
M
)
d
E
T
2
(
1
+
g
1
(
G
T
/
E
T
)
)
(
1
+
(
S
/
s
1
)
)
(
E
T
+
e
T
)
2
]
-
λ
2
[
-
c
2
N
]
-
λ
3
[
-
c
3
E
T
]
-
λ
4
[
p
(
1
+
(
I
10
/
i
2
)
)
(
1
+
(
G
T
/
g
2
)
)
]
-
λ
8
[
p
2
]
-
λ
10
[
p
4
]
,
λ
2
′
(
t
)
=
-
∂
H
∂
N
=
-
λ
1
[
-
c
1
T
]
-
λ
2
[
-
d
N
-
c
2
T
+
p
N
I
2
q
N
+
I
2
-
k
N
(
1
-
e
-
M
)
]
,
λ
3
′
(
t
)
=
-
∂
H
∂
E
T
=
-
λ
1
[
(
g
1
E
T
2
G
+
e
E
T
2
T
+
2
e
g
1
T
G
T
E
T
(
E
T
+
e
T
)
2
(
E
T
+
g
1
G
T
)
2
)
-
d
T
1
+
(
S
/
s
1
)
m
m
m
m
×
(
g
1
E
T
2
G
+
e
E
T
2
T
+
2
e
g
1
T
G
T
E
T
(
E
T
+
e
T
)
2
(
E
T
+
g
1
G
T
)
2
)
]
-
λ
3
[
-
c
3
T
-
k
E
T
(
1
-
e
-
M
)
]
,
λ
4
′
(
t
)
=
-
∂
H
∂
D
U
=
-
λ
4
[
-
γ
1
(
1
+
(
D
U
/
u
1
)
)
2
-
d
D
U
-
k
D
U
(
1
-
e
-
M
)
]
-
λ
5
[
γ
1
(
1
+
(
D
U
/
u
1
)
)
2
]
-
λ
6
[
α
2
I
2
m
H
d
1
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
2
]
,
λ
5
′
(
t
)
=
-
∂
H
∂
D
L
=
-
λ
3
[
α
1
I
m
E
d
1
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
2
]
-
λ
5
[
-
d
D
L
]
-
λ
6
[
α
2
I
m
H
d
1
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
2
]
-
λ
7
[
α
3
I
m
G
d
1
(
1
+
(
S
/
s
2
)
)
(
i
1
+
I
2
)
(
d
1
+
D
L
)
2
]
,
λ
6
′
(
t
)
=
-
∂
H
∂
H
T
=
-
λ
6
[
-
γ
2
S
S
3
+
S
-
d
H
T
-
k
H
T
(
1
-
e
-
M
)
]
-
λ
7
[
γ
2
S
S
3
+
S
]
-
λ
9
[
α
4
(
1
+
(
I
10
/
i
3
)
)
(
1
+
(
S
/
s
4
)
)
]
,
λ
7
′
(
t
)
=
-
∂
H
∂
G
T
=
-
λ
1
[
d
T
E
T
2
g
1
(
E
T
+
e
T
)
(
E
T
+
g
1
G
T
)
2
(
1
+
(
S
/
s
1
)
)
]
-
λ
4
[
-
p
T
g
2
(
1
+
(
I
10
/
i
2
)
)
(
g
2
+
G
T
)
2
]
-
λ
7
[
-
d
G
T
-
k
G
T
(
1
-
e
-
M
)
]
-
λ
8
[
p
1
]
-
λ
10
[
p
3
]
,
λ
8
′
(
t
)
=
-
∂
H
∂
S
=
-
λ
1
[
d
T
E
T
2
s
1
(
E
T
+
e
T
)
(
E
T
+
g
1
G
T
)
(
s
1
+
S
)
2
]
-
λ
3
[
α
1
I
2
D
L
m
E
s
2
(
i
1
+
I
2
)
(
d
1
+
D
L
)
(
s
2
+
S
)
2
]
-
λ
6
[
-
α
2
I
2
(
D
L
+
D
U
)
m
H
s
2
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
(
s
2
+
S
)
2
m
m
m
m
c
-
γ
2
H
s
3
(
s
3
+
S
)
2
α
2
I
2
(
D
L
+
D
U
)
m
H
s
2
(
i
1
+
I
2
)
(
d
1
+
D
L
+
D
U
)
(
s
2
+
S
)
2
]
-
λ
7
[
γ
2
H
T
s
3
(
s
3
+
S
)
2
]
-
λ
8
[
-
d
S
]
-
λ
9
[
-
α
4
H
T
s
4
(
1
+
(
I
10
/
i
3
)
)
(
s
4
+
S
)
2
]
,
λ
9
′
(
t
)
=
-
∂
H
∂
I
2
=
-
λ
2
[
p
N
N
q
N
(
q
N
+
I
2
)
2
]
-
λ
3
[
α
1
D
L
m
E
i
1
(
1
+
(
S
/
s
2
)
)
(
d
1
+
D
L
)
(
i
1
+
I
2
)
2
]
-
λ
6
[
α
2
(
D
L
+
D
U
)
m
H
i
1
(
1
+
(
S
/
s
2
)
)
(
d
1
+
D
L
+
D
U
)
(
i
1
+
I
2
)
2
]
-
λ
7
[
α
3
D
L
m
G
i
1
(
d
1
+
D
L
)
(
i
1
+
I
2
)
2
]
-
λ
9
[
-
d
I
2
]
,
λ
10
′
(
t
)
=
-
∂
H
∂
I
10
=
-
λ
4
[
-
p
T
i
2
(
1
+
(
G
T
/
g
2
)
)
(
i
2
+
I
10
)
2
]
-
λ
9
[
-
α
4
H
T
i
3
(
1
+
(
S
/
s
4
)
)
(
i
3
+
I
10
)
2
]
-
λ
10
[
-
d
I
10
]
,
λ
11
′
(
t
)
=
-
∂
H
∂
M
=
-
λ
1
[
-
k
T
e
-
M
T
]
-
λ
2
[
-
k
N
e
-
M
N
]
-
λ
3
[
-
k
E
T
e
-
M
E
T
]
-
λ
4
[
-
k
D
U
e
-
M
D
U
]
-
λ
6
[
-
k
H
T
e
-
M
H
T
]
-
λ
7
[
-
k
G
T
e
-
M
G
T
]
with transversality conditions
(11)
λ
i
(
t
e
n
d
)
=
0
,
i
=
1,2
,
…
,
11
.
Furthermore, the control functions u
E
*, u
M
* are given by
(12)
u
E
*
=
min
⁡
{
1
,
max
⁡
{
0
,
R
E
}
}
w
h
e
r
e
R
E
=
-
λ
3
ω
1
2
B
,
u
M
*
=
min
⁡
{
1
,
max
⁡
{
0
,
R
M
}
}
w
h
e
r
e
R
M
=
-
λ
11
ω
2
2
C
.Proof.
To determine the adjoint equations and the transversality conditions we use the Hamiltonian (8). The adjoint system results from Pontryagin’s maximum principle [15] are as follows:
(13)
λ
1
′
(
t
)
=
-
∂
H
∂
T
,
λ
2
′
(
t
)
=
-
∂
H
∂
N
,
…
,
λ
11
′
(
t
)
=
-
∂
H
∂
M
with λ
i
(
t
f
)
=
0.
To get the characterization of the optimal control given by (12), solving the equations,
(14)
∂
H
∂
u
E
=
0
,
∂
H
∂
u
M
=
0
on the interior of the control set and using the property of the control space U, we can derive the desired characterization (12).
Here we call formulas (12) for u
* the characterization of the optimal control. The optimal control and the state are found by solving the optimality system, which consists of the system (1), the adjoint system (10), initial conditions at t
=
0, boundary conditions (11), and the characterization of the optimal controls (12). To solve the optimality system, we use the initial and transversality conditions together with the characterization of the optimal control (
u
E
*
,
u
M
*
) given by (12). In addition, the second derivative of the Lagrangian with respect to u
E, u
M, respectively, are positive, which shows that the optimal problem is minimum at controls u
E
*, u
M
*. By substituting the values of u
E
*, u
M
* in the control system (1), we get the following system:
(15)
d
T
*
d
t
=
a
T
*
(
1
-
b
T
*
)
-
c
1
T
*
N
*
-
d
T
*
E
T
*
E
T
*
+
e
T
*
1
(
1
+
g
1
(
G
T
*
/
E
T
*
)
)
(
1
+
(
S
*
/
s
1
)
)
-
k
T
(
1
-
e
-
M
*
)
T
*
,
d
N
*
d
t
=
b
1
-
d
N
N
*
-
c
2
T
*
N
*
+
p
N
N
*
I
2
*
q
N
+
I
2
*
-
k
N
(
1
-
e
-
M
*
)
N
*
,
d
E
T
*
d
t
=
α
1
I
2
*
D
L
*
m
E
(
1
+
(
S
*
/
s
2
)
)
(
i
1
+
I
2
*
)
(
d
1
+
D
L
*
)
-
c
3
E
T
*
T
*
-
k
E
T
(
1
-
e
-
M
*
)
E
T
*
+
min
⁡
{
1
,
max
⁡
{
0
,
R
E
}
}
,
d
D
U
*
d
t
=
p
T
*
(
1
+
(
I
10
*
/
i
2
)
)
(
1
+
(
G
T
*
/
g
2
)
)
-
γ
1
D
U
*
1
+
(
D
U
*
/
u
1
)
-
d
D
U
*
D
U
*
-
k
D
U
*
(
1
-
e
-
M
*
)
D
U
*
,
d
D
L
*
d
t
=
γ
1
D
U
*
1
+
(
D
U
*
/
u
1
)
-
d
D
L
D
L
*
,
d
H
T
*
d
t
=
α
2
I
2
*
(
D
L
*
+
D
U
*
)
m
H
(
1
+
(
S
*
/
s
2
)
)
(
i
1
+
I
2
*
)
(
d
1
+
D
L
*
+
D
U
*
)
-
γ
2
H
T
*
S
*
s
3
+
S
*
-
d
H
T
H
T
*
-
k
H
T
(
1
-
e
-
M
*
)
H
T
*
,
d
G
T
*
d
t
=
γ
2
H
T
*
S
*
s
3
+
S
*
+
α
3
I
2
*
D
L
*
m
G
(
i
1
+
I
2
*
)
(
d
1
+
D
L
*
)
-
d
G
T
G
T
*
-
k
G
T
(
1
-
e
-
M
*
)
G
T
*
,
d
S
*
d
t
=
p
1
G
T
*
+
p
2
T
*
-
d
s
S
*
,
d
I
2
*
d
t
=
α
4
H
T
*
(
1
+
(
I
10
*
/
i
3
)
)
(
1
+
(
S
*
/
s
4
)
)
-
d
I
2
I
2
*
,
d
I
10
*
d
t
=
p
3
G
T
*
+
p
4
T
*
-
d
I
10
*
I
10
*
,
d
M
*
d
t
=
-
d
M
M
*
+
min
⁡
{
1
,
max
⁡
{
0
,
R
M
}
}
with H
* at (
t
,
T
*
,
N
*
,
E
T
*
,
D
U
*
,
D
L
*
,
H
T
*
,
G
T
*
,
S
*
,
I
2
*
,
I
10
*
,
M
*
);
(16)
H
*
=
A
T
*
+
B
min
⁡
{
1
,
max
⁡
{
0
,
R
E
}
}
*
2
+
C
min
⁡
{
1
,
max
⁡
{
0
,
R
M
}
}
*
2
+
λ
1
d
T
*
d
t
+
λ
2
d
N
*
d
t
+
λ
3
d
E
T
*
d
t
+
λ
4
d
D
U
*
d
t
+
λ
5
d
D
L
*
d
t
+
λ
6
d
H
T
*
d
t
+
λ
7
d
G
T
*
d
t
+
λ
8
d
S
*
d
t
+
λ
9
d
I
2
*
d
t
+
λ
10
d
I
10
*
d
t
+
λ
11
d
M
*
d
t
.To find out the optimal control and state, we will numerically solve the above system (15) and (16).


## 3. Numerical Simulations

In this section, we give the numerical results for the effects of chemotherapy and immunotherapy, and optimal control strategy on the tumor-immune model. In our simulations, we consider two initial tumor sizes but keep all other initial state values:(17)
N
(
0
)
=
2.5
×
1
0
8
,
E
T
(
0
)
=
5.268
×
1
0
5
,
D
U
(
0
)
=
4.725
×
1
0
7
,
D
L
(
0
)
=
10
,
H
T
(
0
)
=
1.0536
×
1
0
6
,
G
T
(
0
)
=
1.795
×
1
0
5
,
I
2
(
0
)
=
1173
,
S
(
0
)
=
0
,
I
10
(
0
)
=
0
,
M
(
0
)
=
0
.
Initial conditions of NK cells, CD8+T cells, and IL-2 are from [7]. We note that initial conditions of CD4+T cells and regulatory T cells are derived by initial condition of CD8+T cells. Since the typical ratio of CD8+T cells to CD4+T cells to regulatory T cells is approximately 3 : 6 : 1 [1], we have
(18)
H
T
(
0
)
=
5.268
×
1
0
5
×
2
=
1.0536
×
1
0
6
,
G
T
(
0
)
=
5.268
×
1
0
5
×
1
3
=
1.795
×
1
0
5
.We assume initial conditions ofD
L, S, I
2, and I
10. Our normal dendritic cell counts are in agreement with those in the product insert of Miltenyi Biotec’s blood dendritic cell enumeration kit, which tabulated a total mean dendritic cell count in normal volunteers of 2.8
×
1
0
7 cells/L [16]. In [17], they report that melanoma patients have more circulating dendritic cell per milliliter of blood compared with normal controls. Their mean dendritic cell count for normal volunteers is 64 cells/mL of whole blood, whereas their melanoma patients had mean dendritic cell counts of 108 cells/mL for stage 4 diseases, respectively. So we calculated initial conditions of dendritic cells as follows:
(19)
D
U
(
0
)
=
2.8
×
1
0
7
cells
/
L
×
108
cells
/
mL
64
cells
/
mL
=
4.725
×
1
0
7
cells
/
L
.
In Figure 2(a), with no therapy, the immune system is not able to destroy the tumor cells with initial tumor size of T
(
0
)
=
1
0
8 cells as well as T
(
0
)
=
1
0
7 cells. Figure 2(b) with immunotherapy and Figure 2(c) with chemotherapy show the results of the system (1). In both Figures 2(b) and 2(c), tumor cells are only decreased by the immune system when the initial tumor cell number T
(
0
)
=
1
0
7.(a) No treatment, (b) immunotherapy, (c) chemotherapy, (d) immunotherapy, and chemotherapy.
(a)
(b)
(c)
(d)Finally, Figure2(d) displays the effects of combined therapy on initial tumor sizes T
(
0
)
=
1
0
7 cells and T
(
0
)
=
1
0
8 cells. In this case, the tumor is rapidly destroyed in two cases.For the case of exception immunosuppressive effects, we assume that initial values of regulatory T cell, IL-10, and TGF-β and values of parameters α
3, γ
2, p
1, p
2, p
3, and p
4 are all zeroes. We simulate the optimal controlled model in different scenarios. Firstly, we divide into two cases of the initial value of tumor cells, 107 in Figure 3 and 108 in Figure 4. Secondly, Figures 3 and 4 divided into two cases: the model with or without immunosuppressive effects. The optimality system is solved by using the Runge-kutta fourth-order scheme. The optimal strategy is obtained by solving the state system, the adjoint system, and the transversality conditions. We use Forward-Backward method [18–21] to solve the optimal system.Optimal controls when initial value of tumor cells is1
0
7, (a) A
=
1, B
=
1, C
=
1, (b) A
=
1, B
=
1000, C
=
10, and (c) A
=
1, B
=
10, C
=
1000.
(a)
(b)
(c)Optimal controls when initial value of tumor cells is1
0
8, (a) A
=
1, B
=
1, C
=
1, (b) A
=
1, B
=
1000, C
=
10, and (c) A
=
1, B
=
10, C
=
1000.
(a)
(b)
(c)In our numerical simulation, first we start to solve the state system (1) using the Runge-kutta fourth-order forward in time with a guess for the controls over the simulation time. Then, using the current iteration of the state equations in the system (1), the adjoint equations in the system (10) are solved by a backward method with the transversality conditions (11). We update the controls by using a convex combination of the controls in the previous iterations if the values of unknowns at the previous iteration are very close to the ones at the present iteration. In Figures 3 and 4, (a) represents scenarios for the case when cost for two therapies is the same ratio. (b) represents scenarios for the case when the cost of immunotherapy is more expensive than the cost of chemotherapy. (c) represents scenarios for the case when the cost of chemotherapy is more expensive than the cost of immunotherapy. According to these scenarios, we certify that when the model has immunosuppressive effects, the optimal treatment strategy may change based on treatment cost.In Figures3(c) and 4(c), our simulation results show that if the cost of chemotherapy is more expensive than that of immunotherapy, then the optimal treatment strategy for the model without immunosuppressive effects needs to be taken for a longer time for chemotherapy comparing immunotherapy. On the other hand, when the model has immunosuppressive effects and the cost of chemotherapy is more expensive than that of immunotherapy, the optimal treatment strategy for the model needs to be taken for a longer time for immunotherapy comparing chemotherapy. In other words, the optimal treatment strategy may be changed by immunosuppressive effects.In order to find the parameter factors that exert a strong impact on model outcome, we use the numerical parameter sensitivity analysis. For the sensitivity analysis, one parameter value in the model is increased and decreased by 20 percent and the other parameter values are fixed. After 5 days, we plot tumor sizes depending on the model parameters. In Figure5, the solid red line and blue line represent change rates of tumor cell numbers when a special parameter value was decreased 20 percent and when parameter value was increased 20 percent, respectively. From this we can check that during 5 days, tumor size is highly sensitive to parameters a (tumor growth), d (strength of immune system), and α
1 (rate of IL-2 and DC induced CD8+T cell activation) in order of list.Figure 5
Numerical parameter sensitivity.


## 4. Discussion and Conclusions

We constructed a mathematical model describing tumor-immune interactions under immune suppression. From this model, we suggested a treatment protocol for each of chemotherapy and immunotherapy. We investigated the dynamics which effectiveness and efficiency of two therapies by changing of initial tumor cell numbers. We use optimal control techniques and numerical simulation to find a combined therapy strategy for treatment of the tumor. By using Pontryagin’s maximum principle, we derived the necessary conditions of optimality for the control system.To analyse the parameter sensitivity, we plotted the percentage change in tumor size from day zero to day five as a result of changing each of the model parameters by 20% in both directions. From such sensitivity analysis, we found out some special parameters that have a strong influence on tumor growth. Our optimal control experiments demonstrated how chemotherapy and immunotherapy might be combined for more effective treatment. We showed that combined therapy is more effective than each therapy; that is, the number of tumor cells decreases in special parameter sets. In addition, we found out that if the cost of chemotherapy is more expensive than that of immunotherapy, then the optimal treatment strategy for the model without immunosuppressive effects needs to be taken for a longer time for chemotherapy comparing immunotherapy. On the other hand, when the model has immunosuppressive effects and the cost of chemotherapy is more expensive than that of immunotherapy, the optimal treatment strategy for the model needs to be taken for a longer time for immunotherapy comparing chemotherapy.Even though the methodology of this paper is standard, we provided a key process to develop the optimal control problem related to the cancer model. Based on the parameter sensitive analysis, we formulated an optimal control problem related to the tumor-immune interaction under immune suppression, which is a fresh idea in the optimal control problems. Moreover, this paper gives theoretical and experimental results in the sense of mathematical analysis but if we would have any field data for some cancer treatment and patients, then our results will be applied to the cancer model as well as some other disease models very well.



---

*Source: 206287-2014-07-23.xml*