# State-Dependent Impulsive Control Strategies for a Tumor-Immune Model

Authors: Kwang Su Kim; Giphil Cho; Lin-Fei Nie; Il Hyo Jung
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

Controlling the number of tumor cells leads us to expect more efficient strategies for treatment of tumor. Towards this goal, a tumor-immune model with state-dependent impulsive treatments is established. This model may give an efficient treatment schedule to control tumor’s abnormal growth. By using the Poincaré map and analogue of Poincaré criterion, some conditions for the existence and stability of a positive order-1 periodic solution of this model are obtained. Moreover, we carry out numerical simulations to illustrate the feasibility of our main results and compare fixed-time impulsive treatment effects with state-dependent impulsive treatment effects. The results of our simulations say that, in determining optimal treatment timing, the model with state-dependent impulsive control is more efficient than that with fixed-time impulsive control.


---

## Body

## 1. Introduction

Cancer is a class of diseases characterized by out of control cell growth. Abnormal cells divide without control and are able to invade other tissues. It can also spread to other parts of the body through the blood and lymph systems. Cancer is a leading cause of death worldwide. In the past decades, doctors are trying to cure cancer but it has been difficult. Among many treatments for cancer, immunotherapy is a treatment method by using immune system in a body to fight cancer cells. Immune system recognizes cancer cells and leads to destruction of cancer cells before cancer cells are big enough to see. However, since the mechanism of immune system in the body is very complex, it is difficult to find an efficient treatment schedule to eradicate some cancers. So, in order to find such schedule, a mathematical model describing dynamics of tumor-immune interactions would provide a new strategy for treatment of cancer.In the past, mathematical modeling about tumor-immune interactions has been studied by many scholars (see [1–3]). Particularly, Gałach [4] provided three kinds of different models which describe a competition between the tumor and immune cells by using a mathematical model with delay. Actually they investigated the stability of the steady state and observed a state of the “returning" tumor in the model with time delay. In several papers, there are several mathematical models that describe treatment effects: immunotherapy, chemotherapy, radiation therapy, and so on [1, 5].In real world, such treatments need not once but periodic schedules. In particular, the impulsive differential equations can describe well such periodic treatments and pulse vaccination. Disease mathematical models using time impulsive differential equations have been studied by many researchers. For example, Meng and Chen [6] suggested a delayed SIR epidemic model with impulsive effect and analyzed the dynamic behaviors of the model. For SI model, Pei et al. [7] analyzed the dynamic behaviors by an impulsive delayed SI model. In Gao et al. [8], they proposed the periodic therapy from antiretroviral drugs for HIV by using impulsive differential equations. Qiao et al. [9] suggested a Hepatitis B virus infection model with impulsive vaccination and time delay. Also they showed the existence of a periodic solution of the model and derived sufficient conditions that Hepatitis B virus will be eliminated or persistent. Huang et al. [10] suggested two novel mathematical models by using an impulsive control term for injection of insulin and showed the existence and stability of the periodic solution of the model.A noticed fact is that the treatment of different cancers has different ways. For some cancers that cannot be cured completely, a natural idea is therefore to keep the density of tumor cells at a low level to slow the progression of cancers. Based on this idea, Tang et al. [11] proposed an SIR epidemic model with state-dependent pulse control strategies and demonstrated that the combination of pulse vaccination and treatment is optimal in terms of cost under certain conditions, and the existence and stability of periodic solution with the maximum value of the infective no greater than threshold value RL are studied. Nie et al. [12] proposed an SIR epidemic model with state-dependent pulse vaccination and proved the existence and the stability of a positive periodic solution for the model. Additionally, in tumor treating, people have high death rate if they have tumor cells too much. So it is important that we control to reduce the number of tumor cells below the value of hazardous. Therefore, to describe tumor-immune systems with treatments, state-dependent impulsive differential equation is more appropriate than fixed-time dependent impulsive differential equation. However, the mathematical models for tumor-immune interactions with state-dependent impulsive treatments are rare. Additionally, state-dependent impulsive control strategies are also used in the treatment and control of viral diseases. For example, Nie et al. [13] proposed two novel virus dynamics models with state-dependent control strategies and analytically proved the existence and orbit stability of semitrivial periodic solution and positive periodic solution. Theoretical results showed that the density of infected cells or free virus can be kept within a low level over a long period of time by adjusting control strength. Tang et al. [14] constructed a mathematical model with state-dependent impulsive interventions where the comprehensive therapy involving combining surgery with immunotherapy is considered. In it, authors examined the global dynamics of the immune tumour system with state-dependent feedback control, including the existence and stability of the semitrivial order-1 periodic solution and the positive order-k periodic solution.Inspired by the above discussion, the objective of this study is to develop a tumor-immune model with state-dependent impulsive treatment: immunotherapy and chemotherapy by using state-dependent impulsive differential equations, find some conditions for the existence and stability of periodic solution, and compare between state-dependent impulsive control and fixed-time impulsive control for control of tumor.The rest of this paper is organized as follows: In Section2, we give a mathematical model for the tumor-immune system with state-dependent impulsive control strategies, which is based on the Gałach’s model [4], and describe the parameters used in this model. In Section 3, for the main results, we show the existence and stability of a positive order-1 periodic solution. In Section 4, we carry out some numerical simulations to illustrate the main results and compare a fixed-time pulse treatment with a state-dependent pulse treatment. Finally, the conclusion is given in Section 5.


## 2. Mathematical Model of Tumor-Immune Dynamics

Gałach’s autonomous tumor-immune dynamics model is given by the following differential equations [4]: (1)
d
E
t
d
t
=
σ
+
ω
T
t
E
t
-
δ
E
t
,
d
T
t
d
t
=
α
T
t
1
-
β
T
t
-
T
t
E
t
.Gałach changed the model proposed in [3] by replacing the Michaelis-Menten form of function with a Lotka-Volterra form. Here E
(
t
), T
(
t
) represent the numbers of immune cells and tumor cells, respectively, and parameter σ
>
0 denotes the rate of immune cell proliferation, parameter ω
>
0 describes the immune response to the appearance of the tumor cells, parameter δ
>
0 is a rate of immune cell turnover, parameter α
>
0 is a growth rate of tumor cells, and parameter β
>
0 denotes the inverse of tumor cell’s carrying capacity; 1
/
β is the mean tumor cell’s carrying capacity.On the stability of an equilibrium, Gałach [4] has the following results.Lemma 1.
(i) Ifω
>
0 and α
δ
<
σ, then model (1) has a trivial equilibrium P
0
(
σ
/
δ
,
0
) which is stable.
(ii) Ifω
>
0 and α
δ
>
σ, then model (1) has an unstable equilibrium P
0 and a stable positive equilibrium P
1
(
E
1
,
T
1
), where (2)
E
1
=
α
β
δ
+
ω
-
Δ
2
α
β
ω
,
T
1
=
-
α
β
δ
-
ω
+
Δ
2
ω
,
Δ
=
α
2
β
δ
-
ω
2
+
4
α
β
σ
ω
>
0
.
(iii) Ifω
<
0, α
δ
>
σ and α
(
β
δ
-
ω
)
2
+
4
β
σ
ω
>
0, then model (1) has an unstable equilibrium P
0 and a stable positive equilibrium P
1
(
E
1
,
T
1
), where (3)
E
1
=
α
β
δ
+
ω
-
Δ
2
α
β
ω
,
T
1
=
-
α
β
δ
-
ω
+
Δ
2
ω
,
Δ
=
α
2
β
δ
-
ω
2
+
4
α
β
σ
ω
>
0
.It is well known that immunotherapy and chemotherapy are effective to control the tumor cells. Chemotherapy attacks not only tumor cells but also normal cells and immunotherapy is applied by injection of immune cells. Therefore, cancer patients need periodic treatments for a long time. For the reality, a state-dependent impulsive control strategy for tumor abnormal growth is proposed here, rather than the usual continuous control or pulse fixed-time control strategy. The treatment strategy is based on the following assumption.(
A
)
When the amount of tumor cell reaches the hazardous threshold valueH
(
>
0
) at time t
i
(
H
) at the ith time, immunotherapy and chemotherapy are taken and the amounts of immune cells and tumor cells abruptly turn to (
1
-
p
)
E
(
t
i
-
(
H
)
)
+
τ and (
(
1
-
q
)
T
(
t
i
-
(
H
)
)
), respectively, where p
,
q
∈
(
0,1
) are the chemotherapy intensity and τ
∈
(
0,1
) is the amounts of immunotherapy, respectively.Under assumption(
A
), we build the control process that is modeled by the following differential equation with state-dependent impulsive effects:(4)
d
E
t
d
t
=
σ
+
ω
T
t
E
t
-
δ
E
t
d
T
t
d
t
=
α
T
t
1
-
β
T
t
-
T
t
E
t
T
<
H
,
Δ
E
t
=
E
t
-
E
t
-
=
-
p
E
t
-
+
τ
Δ
T
t
=
T
t
-
T
t
-
=
-
q
T
t
-
T
=
H
.Remark 2.
We indicate that we do not know a priori time of pulse treatments of chemotherapy and immunotherapy since it is taken at the time when the number of tumor cellsT reaches the threshold value H. So the pulse treatments of chemotherapy and immunotherapy time obviously depend on the solution and our model adapts state-dependent pulse treatments of chemotherapy and immunotherapy.Remark 3.
The control parametersp, q, τ and hazardous threshold H rely heavily on the characteristics of various cancer types. The choices of these values are very important, which are closely geared to the growth or elimination of cancer cells.Lemma 4.
Each component of solution(
E
(
t
)
,
T
(
t
)
) of model (4) with nonnegative initial condition (
E
0
,
T
0
)
=
(
E
(
t
0
)
,
T
(
t
0
)
) is positive for every t
≥
0.Proof.
For any nonnegative initial condition(
E
0
,
T
0
), suppose that (
E
(
t
)
,
T
(
t
)
) is a solution of model (4) with the positive initial conditions (
E
0
,
T
0
); then relation of solution (
E
(
t
)
,
T
(
t
)
) with line K
:
T
=
H is one of the following cases.
(i) The Solution
(
E
(
t
)
,
T
(
t
)
)
Intersects Line
K
Finite Times. Since (
E
(
t
)
,
T
(
t
)
) intersects line K finite times and the positive equilibrium P
1
(
E
1
,
T
1
) is stable, E
(
t
)
≥
0, T
(
t
)
≥
0 for every t
≥
0.
(ii) The Solution
(
E
(
t
)
,
T
(
t
)
)
Intersects Line
K
Infinite Times. Suppose that solution (
E
(
t
)
,
T
(
t
)
) intersects with line K
:
T
=
H at times t
1
,
t
2
,
t
3
,
…
,
t
k
,
…, so that l
i
m
k
→
∞
t
k
=
∞. If the conclusion of Lemma 4 is not true, then there is a t
∗
>
t
0 such that min
⁡
{
E
(
t
∗
)
,
T
(
t
∗
)
}
<
0, and E
(
t
)
>
0, T
(
t
)
>
0 for every t
0
≤
t
<
t
∗. For this t
∗, there is a positive integer n such that t
n
-
1
≤
t
∗
≤
t
n. (i)
IfE
(
t
∗
)
<
0, T
(
t
∗
)
≥
0, then from the first and third equations of model (4), we obtain (5)
E
t
∗
≥
1
-
p
n
-
1
E
t
0
exp
⁡
∫
t
0
t
∗
ω
T
s
-
δ
d
s
≥
0
,which leads to a contradiction to E
(
t
∗
)
<
0.
(ii)
IfT
(
t
∗
)
<
0, E
(
t
∗
)
≥
0, then from the second and fourth equations of model (4), we have (6)
T
t
∗
≥
T
0
1
-
β
T
0
e
-
α
t
∗
+
β
T
0
≥
0
,which leads to a contradiction to T
(
t
∗
)
<
0.So we know that t
∗ does not exist such that min
⁡
{
E
(
t
∗
)
,
T
(
t
∗
)
}
<
0. Therefore, we can prove that E
(
t
)
≥
0
,
T
(
t
)
≥
0 for every t
≥
0.Next, for the convenience of statement in the remainder of this paper, we introduce the following definitions. LetΓ be an arbitrary set in R
2 and Y be an arbitrary point in R
2. The distance between Y and Γ is defined by d
(
Y
,
Γ
)
=
inf
Y
0
∈
Γ
⁡
Y
-
Y
0
.Definition 5 (orbital stability [15]).
TrajectoryO
+
(
X
0
,
t
0
) is said to be orbitally stable if, given ε
>
0, there exists a constant δ
=
δ
(
ε
)
>
0 such that, for any solution X
∗
(
t
) of model (4), when d
(
X
∗
(
t
0
)
,
O
+
(
X
0
,
t
0
)
)
<
δ, one has d
(
X
∗
(
t
)
,
O
+
(
X
0
,
t
0
)
)
<
ε for all t
>
t
0.Definition 6 (orbitally asymptotical stability [15]).
TrajectoryO
+
(
X
0
,
t
0
) is said to be orbitally asymptotically stable if it is orbitally stable; there exists a constant η
>
0 such that, for any other solution X
∗
(
t
) of model (4), when d
(
X
∗
(
t
0
)
,
O
+
(
X
0
,
t
0
)
)
<
η, then lim
t
→
∞
⁡
d
(
X
∗
(
t
)
,
O
+
(
X
0
,
t
0
)
)
=
0.To analyze the dynamic behaviors of model (4), we define two sections to the phase space of model (4) by (7)
Γ
q
≔
E
,
T
:
E
≥
0
,
T
=
1
-
q
H
,
Γ
H
≔
E
,
T
:
E
≥
0
,
T
=
H
.For any point a
n
(
E
n
,
H
) on section Γ
H, suppose that trajectory O
+
(
a
n
,
t
n
) starting from point a
n intersects section Γ
H infinitely many times. That is, trajectory O
+
(
a
n
,
t
n
) jumps to point a
n
+
(
(
1
-
p
)
E
n
+
τ
,
(
1
-
q
)
H
) on section Γ
q due to impulsive effects Δ
E
(
t
)
=
-
p
E
(
t
-
)
+
τ and Δ
T
(
t
)
=
-
q
H. Furthermore, trajectory O
+
(
a
n
,
t
n
) intersects section Γ
H at point a
n
+
1
(
E
n
+
1
,
H
) and then jumps to point a
n
+
1
+
(
(
1
-
p
)
E
n
+
1
+
τ
,
(
1
-
q
)
H
) on section Γ
q again (see Figure 1(b)). Repeating this process, we have two impulsive point sequences a
n
+
1
-
p
E
n
+
τ
,
1
-
q
H and a
n
E
n
,
H, where E
n
+
1 is determined by E
n, p, q, τ, H and n
=
1,2
,
…. We then define the Poincaré map of section Γ
H as follows:(8)
E
n
+
1
=
P
E
n
,
p
,
q
,
τ
,
H
.Figure 1
The difference of dynamical behaviors for models (1) and (4).
(a)
(b)Definition 7.
A trajectoryO
+
(
a
n
,
t
n
) of model (4) is said to be order-k periodic if there exists a positive integer k
≥
1 such that k is the smallest integer for E
n
+
k
=
E
n.Next, we consider the following autonomous model with impulsive effects:(9)
d
x
d
t
=
f
x
,
y
,
d
y
d
t
=
g
x
,
y
,
φ
x
,
y
≠
0
,
Δ
x
=
ξ
x
,
y
,
Δ
y
=
η
x
,
y
,
φ
x
,
y
=
0
,where f
(
x
,
y
) and g
(
x
,
y
) are continuous differential functions defined on R
2 and φ
(
x
,
y
) is a sufficiently smooth function with ∇
φ
(
x
,
y
)
≠
0. Let (
μ
(
t
)
,
ν
(
t
)
) be a positive T-periodic solution of the model (9). The following result comes from Corollary  2 of [16].Lemma 8 (analogue of Poincaré criterion).
If the Floquet multiplierμ satisfies condition μ
<
1, where (10)
μ
=
∏
j
=
1
n
κ
j
exp
⁡
∫
0
T
∂
f
μ
t
,
ν
t
∂
x
+
∂
g
μ
t
,
ν
t
∂
y
d
twith (11)
κ
j
=
∂
η
/
∂
y
∂
φ
/
∂
x
-
∂
η
/
∂
x
∂
φ
/
∂
y
+
∂
φ
/
∂
x
f
+
+
∂
ξ
/
∂
x
∂
φ
/
∂
y
-
∂
ξ
/
∂
y
∂
φ
/
∂
x
+
∂
φ
/
∂
y
g
+
∂
φ
/
∂
x
f
+
∂
φ
/
∂
y
gand f, g, ∂
ξ
/
∂
x, ∂
ξ
/
∂
y, ∂
η
/
∂
x, ∂
η
/
∂
y, ∂
φ
/
∂
x, and ∂
φ
/
∂
y have been calculated at the points (
μ
(
τ
j
)
,
ν
(
τ
j
)
), f
+
=
f
(
μ
(
τ
j
+
)
,
ν
(
τ
j
+
)
), g
+
=
g
(
μ
(
τ
j
+
)
,
ν
(
τ
j
+
)
), and τ
j
(
j
∈
N
) is the time of the jth jump, then, (
μ
(
t
)
,
ν
(
t
)
) is orbitally asymptotically stable.


## 3. Main Results

The first result is on the existence of a positive order-1 periodic solution for model (4).Theorem 9.
For anyq
,
p
∈
(
0,1
), model (4) admits a positive order-1 periodic solution.Proof.
Let the pointc
1
(
ε
,
(
1
-
q
)
H
)
∈
Γ
q for sufficiently small ε with ε
≤
τ. In view of the geometrical structure of the phase space of model (4), trajectory O
+
(
c
1
,
t
0
) of model (4) starting from initial point c
1 intersects section Γ
H at point d
1
(
E
1
,
H
). At point d
1, the trajectory O
+
(
c
1
,
t
0
) jumps to point c
2
(
(
1
-
p
)
E
1
+
τ
,
(
1
-
q
)
H
) on section Γ
q due to impulsive effects Δ
E
(
t
)
=
-
p
E
(
t
-
)
+
τ and Δ
T
(
t
)
=
-
q
T
(
t
-
). Furthermore, the trajectory O
+
(
c
1
,
t
0
) intersects section Γ
H at point d
2
(
E
2
,
H
).
Sinceε
≤
τ, then ε
<
(
1
-
p
)
E
1
+
τ. It follows that c
1 is on the left of c
2. We claim that d
1 is also on the left of d
2. In fact, if d
1 is on the right of d
2, then orbits c
1
d
1
~ and c
2
d
2
~ intersect at the point (
E
0
,
T
0
). This shows that there are two different solutions which start from the point (
E
0
,
T
0
). This contradicts with the uniqueness of solutions for model (4). So by (8), we have E
2
=
P
(
E
1
,
p
,
q
,
τ
,
H
) and (12)
P
E
1
,
p
,
q
,
τ
,
H
-
E
1
=
E
2
-
E
1
>
0
.
On the other hand, suppose that the curvel
1
:
α
-
α
β
T
-
E
=
0 intersects section Γ
q at point a
0
(
α
-
(
1
-
q
)
α
β
H
,
(
1
-
q
)
H
). Then trajectory O
+
(
a
0
,
t
0
) starting from the initial point a
0 intersects section Γ
H at point b
1
(
E
^
1
,
H
) and then jumps to point a
1
(
(
1
-
p
)
E
^
1
+
τ
,
(
1
-
q
)
H
) on section Γ
q and finally reaches a point b
2
(
E
^
2
,
H
) in section Γ
H again. If there is a positive constant τ
∗ such that (
1
-
p
)
E
^
1
+
τ
∗
=
α
-
(
1
-
q
)
α
β
H, then a
1 coincides with a
0 for τ
=
τ
∗
∈
(
0,1
). That is, b
1 coincides with b
2. Otherwise, a
1 is on the left of a
0 for (
1
-
p
)
E
^
1
+
τ
<
α
-
(
1
-
q
)
α
β
H and a
1 is on the right of a
0 for (
1
-
p
)
E
^
1
+
τ
>
α
-
(
1
-
q
)
α
β
H. However, from the geometrical structure of a phase space of model (4), b
2 is on the left of b
1 for any τ
∈
(
0
,
τ
∗
)
∪
(
τ
∗
,
+
∞
).
To sum up, we get, from the above discussions, that(i)
whenE
^
1
=
E
^
2, model (4) has a positive order-1 periodic solution;
(ii)
whenE
^
2
<
E
^
1, (13)
P
E
^
1
,
p
,
q
,
τ
,
H
-
E
^
1
=
E
^
2
-
E
^
1
<
0
.
By (12) and (13), it follows that the Poincaré map (8) has a fixed point. This means that model (4) has a positive order-1 periodic solution. The proof is complete.The next result is on the orbital stability of a positive order-1 periodic solution for model (4).Theorem 10.
Let(
ϕ
(
t
)
,
ψ
(
t
)
) be a positive order-1 periodic solution of model (4) with period ϖ. If(14)
μ
=
κ
exp
⁡
∫
0
ϖ
ω
ψ
t
-
δ
d
t
<
1
,where (15)
κ
=
1
-
p
1
-
q
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
ϖ
-
τ
α
-
α
β
H
-
ϕ
ϖ
,then (
ϕ
(
t
)
,
ψ
(
t
)
) is locally orbitally asymptotically stable.Proof.
Suppose that(
ϕ
(
t
)
,
ψ
(
t
)
) intersects sections Γ
p and Γ
H at points K
+
(
(
1
-
p
)
ϕ
(
ϖ
)
+
τ
,
(
1
-
p
)
H
) and E
(
ϕ
(
ϖ
)
,
H
), respectively. Comparing with model (9), we have (16)
f
E
,
T
=
σ
+
ω
T
E
-
δ
E
,
g
E
,
T
=
α
T
1
-
β
T
-
T
E
,and ξ
(
E
,
T
)
=
-
p
E
+
τ, η
(
E
,
T
)
=
-
q
T, φ
(
E
,
T
)
=
T
-
H, (
ϕ
(
ϖ
)
,
ψ
(
ϖ
)
)
=
(
ϕ
(
ϖ
)
,
H
), and (
ϕ
(
T
∗
+
)
,
ψ
(
T
∗
+
)
)
=
(
(
1
-
p
)
ϕ
(
ϖ
)
+
τ
,
(
1
-
q
)
H
). Thus(17)
∂
f
∂
E
=
ω
T
-
δ
,
∂
g
∂
T
=
α
-
2
α
β
T
-
E
,
∂
ξ
∂
E
=
-
p
,
∂
η
∂
T
=
-
q
,
∂
φ
∂
T
=
1
,
∂
ξ
∂
T
=
∂
η
∂
E
=
∂
φ
∂
E
=
0
.Furthermore, it follows from (17) that(18)
κ
=
∂
η
/
∂
T
∂
φ
/
∂
E
-
∂
η
/
∂
E
∂
φ
/
∂
T
+
∂
φ
/
∂
E
f
+
+
∂
ξ
/
∂
E
∂
φ
/
∂
T
-
∂
ξ
/
∂
T
∂
φ
/
∂
E
+
∂
φ
/
∂
T
g
+
∂
φ
/
∂
E
f
+
∂
φ
/
∂
I
g
=
1
-
p
g
+
ϕ
T
∗
+
,
ψ
T
∗
+
g
ϕ
ϖ
,
ψ
ϖ
=
1
-
p
1
-
q
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
ϖ
-
τ
α
-
α
β
H
-
ϕ
ϖ
,
μ
=
κ
exp
⁡
∫
0
T
ω
ψ
t
-
δ
+
α
-
2
α
β
ψ
t
-
ϕ
t
d
t
.
On the other hand, integrating both sides of the second equation of model (4) along the orbit c
+
c
~, we have(19)
ln
⁡
1
1
-
q
=
∫
1
-
q
H
H
d
T
T
=
∫
0
T
∗
α
-
2
α
β
T
-
E
d
t
=
∫
0
ϖ
α
-
2
α
β
ψ
t
-
ϕ
t
d
t
.From (18)-(19), we can obtain (20)
μ
=
1
-
p
1
-
q
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
ϖ
-
τ
α
-
α
β
H
-
ϕ
ϖ
1
1
-
q
exp
⁡
∫
0
ϖ
ω
ψ
t
-
δ
d
t
=
1
-
p
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
ϖ
-
τ
α
-
α
β
H
-
ϕ
ϖ
exp
⁡
∫
0
ϖ
ω
ψ
t
-
δ
d
t
.
By condition (14), we see that model (4) satisfies all conditions of Lemma 8. It then follows from Lemma 8 that the order-1 periodic solution (
ϕ
(
t
)
,
ψ
(
t
)
) of model (4) is locally orbitally asymptotically stable and has asymptotic phase property. This completes the proof.From the proof of Theorem10, integrating both sides of the first equation of model (4) along orbit c
+
c
~, it follows that (21)
exp
⁡
∫
0
ϖ
ω
ψ
t
-
δ
d
t
≤
exp
⁡
∫
0
ϖ
ω
H
-
δ
d
t
≤
exp
⁡
ω
H
-
ω
ϖ
.Therefore, we have the consequence of Theorem10 as Corollary 11 below.Corollary 11.
Let(
ϕ
(
t
)
,
ψ
(
t
)
) be a positive order-1 periodic solution of model (4) with periodic ϖ. If (22)
μ
=
1
-
p
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
ϖ
-
τ
α
-
α
β
H
-
ϕ
ϖ
exp
⁡
ω
H
-
δ
ϖ
<
1
,then (
ϕ
(
t
)
,
ψ
(
t
)
) is locally orbitally asymptotically stable.Theorem 12.
Ifτ
<
α
-
α
β
(
1
-
q
)
H, then model (4) has a positive order-1 periodic solution which is orbitally asymptotical stable.Proof.
Suppose that the curvel
1
:
α
-
α
β
T
-
E
=
0 intersects section Γ
q at point a
0
(
α
-
(
1
-
q
)
α
β
H
,
(
1
-
q
)
H
). Let α
-
(
1
-
q
)
α
β
H
=
λ
0. Then, since τ
<
λ
0, there is a positive constant λ
∗ such that (
1
-
p
)
λ
∗
+
τ
=
λ
0. For any two points a
i
(
E
i
,
H
) and a
j
(
E
j
,
H
) on section γ
H, E
i
,
E
j
∈
(
0
,
λ
∗
) and E
i
<
E
j. In view of the impulsive treatment effects, point a
i
+
(
(
1
-
p
)
E
i
+
τ
,
(
1
-
q
)
H
) is on the left of point a
j
+
(
(
1
-
p
)
E
j
+
τ
,
(
1
-
q
)
H
). Therefore, from the geometrical construction of the phase space of model (4) we have (23)
E
i
+
1
<
E
j
+
1
<
λ
∗
.Then, for any E
0
∈
(
0
,
λ
∗
), from the Poincaré map (8) of section Γ
H we have E
1
=
F
(
E
0
,
p
,
q
,
τ
,
H
), E
2
=
F
(
E
1
,
p
,
q
,
τ
,
H
), and E
n
+
1
=
F
(
E
n
,
p
,
q
,
τ
,
H
)
(
n
=
2,3
,
…
). In particular, if E
0
=
E
1, then model (4) has a positive order-1 periodic solution.
Next, we discuss the general circumstance; that is,E
0
≠
E
1
≠
E
2
≠
E
3
≠
⋯
≠
E
k
(
k
≥
2
)
.
(1)
IfE
1
<
E
0, then it follows that E
2
<
E
1, E
3
<
E
2. Repeating the above process, we have (24)
0
<
⋯
<
E
5
<
E
4
<
E
3
<
E
2
<
E
1
<
E
0
.
(2)
IfE
0
<
E
1, then it follows that E
1
<
E
2, E
2
<
E
3. Repeating the above process, we have (25)
E
0
<
E
1
<
E
2
<
E
3
<
E
4
<
E
5
<
⋯
<
λ
∗
.
Further in Case (1) it followslim
k
→
∞
⁡
E
k
=
δ
1, where 0
<
δ
1
<
E
0. Similarly, in Case (
2
)
lim
k
→
∞
⁡
S
k
=
δ
2, where S
0
<
δ
2
<
λ
∗. So model (4) has orbitally asymptotically stable positive order-1 periodic solution.Remark 13.
Obviously, the treatment strengthp of chemotherapy is closely related to the stable region of an order-1 solution. If p is close to 1, then λ
∗ is large, so by Theorem 12, the local asymptotic orbital stability region is large. Otherwise, if p is close to 0, then the local asymptotic orbital stability region is small. That is, the local asymptotic orbital stability region depends on value of p.Remark 14.
We study the existence and the stability of a positive order-1 periodic solution for the case (ii) of Lemma1. For the case (ii) of Lemma 1, model (4) has an unstable equilibrium P
0 and a stable positive equilibrium P
1, so that we can obtain similar results of case (ii) of Lemma 1 by using the same methods. Therefore, we only consider numerical simulations for case (iii) of Lemma 1 in next section.


## 4. Numerical Results

In this work, we investigate a tumor-immune model with state-dependent impulsive immunotherapy and chemotherapy. By the Poincaré map, analogue of Poincaré criterion, and qualitative analysis method, we give the criteria for the existence and stability of the positive periodic solution for this model. In order to illustrate the implications of our results and the feasibility of the state-dependent pulse feedback control strategy, we consider a tumor-immune model with state-dependent impulsive immunotherapy and chemotherapy. For simulation, we use parameter values in Table1.Table 1
Equations and parameter descriptions.
State
Term
Parameter
Description
E
(
t
)
σ
σ
=
0.1181 [17]
Production of effect cell
ω
T
E
ω
=
0.04 [18]
Stimulated accumulation of effect cell by tumor cell
δ
E
δ
=
1 (assumed)
Turnover of effect cell
T
(
t
)
α
T
(
1
-
β
T
)
α
=
1.636, β
=
0.002 [18]
Tumor growth term
T
E
Effect cell induced tumor deathFirstly, it is easy to obtain that the model (4) without impulsive effects has a stable focus (
0.8333,3.1667
). The plots in Figure 1(a) show the tumor-immune dynamical behaviors in patients without therapies, which implies that cancer is deadly. Now, we apply chemotherapy and immunotherapy to control tumor cells. That is, we choose control parameters p
=
0.1, q
=
0.2, τ
=
0.22, and H
=
0.7; the plots in Figure 1(b) show that model (4) has a positive order-1 periodic solution. This indicates that tumor cells might be controlled below the threshold value by state-dependent impulsive control strategies.Nextly, we choose control parametersp
=
0.1, q
=
0.4, τ
=
0.3, and H
=
0.7 and other model parameters are fixed in Table 1. From Theorem 9, model (4) admits a positive order-1 periodic solution (
ϕ
(
t
)
,
ψ
(
t
)
) starting from point (
1.6767,0.42
)
∈
Γ
p and intersection section Γ
H at point (
1.5296,0.7
)
=
(
ϕ
(
T
)
,
ψ
(
T
)
). This is shown in Figure 2(a). Further, this can be easily calculated by Corollary 11: (26)
μ
=
1
-
p
α
-
α
β
1
-
q
H
-
1
-
p
ϕ
T
-
τ
α
-
α
β
H
-
ϕ
T
exp
⁡
ω
H
T
≈
0.4753
<
1
.This shows that the positive order-1 periodic solution is orbitally asymptotically stable. This is shown in Figure 2(b). The plots in Figures 3(a) and 3(b) show that the existence and stability of an order-1 periodic solution of model (4) for case (iii) in Lemma 1.Figure 2
The existence and stability of order-1 periodic solution for model (4) with p
=
0.1, q
=
0.4, and τ
=
0.2: (a) the existence; (b) the orbital asymptotical stability.
(a)
(b)Figure 3
The existence and stability of positive order-1 periodic solution of model (4) for case (iii) in Lemma 1.
(a)
(b)Thirdly, we investigate how the strength of chemotherapy and immunotherapy affects the control of tumor cell numbers and existence and stability of the periodic solutions. The valuep means the side effects of chemotherapy to immune cells. Now, we change only the values of p such that 0.06,0.08,0.1 but the other parameter values are fixed. In Figure 3(a), our simulation results show that if the value of p increases, then periodic T of positive order-1 periodic solution for model (4) decreases. In particular, if p
=
0.06, then we can see that the tumor cell numbers fall off more than the other cases. Secondly, let τ denote amounts of injected immune cells for immunotherapy. In Figure 3(b), we can see that the periodic T of an order-1 periodic solution decreases with decreasing value of τ. The fourth injection that the value of τ is 0.235 and the fifth injection that the value of τ is 0.195 happen simultaneously. Thirdly, the value q implies the strength of chemotherapy to tumor cells. Then there are some differences of simulation results between the value q and the values p and τ. Numerical simulation shows that the period T of positive order-1 periodic solution for model (4) increases with the increase of the strength of chemotherapy to tumor cell q (see Figure 3(c)). And if the value of p increases, tumor cell’s lower bounded line is lowered in Figure 4(a). These show that the recurrence time for tumor cells can be prolonged by increasing strength of chemotherapy to tumor cells q, amounts of injected immune cells τ, and decreasing side effects of chemotherapy to immune cells p by through simulation results (see Figures 4(b) and 4(c)).Figure 4
The effect of control parametersp, q, and τ for the dynamical behaviors of model (4): (a) p
=
0.1, 0.08, and 0.06, respectively; (b) τ
=
0.235, 0.215, and 0.195, respectively; (c) q
=
0.2, 0.3, and 0.4, respectively.
(a)
(b)
(c)Finally, we investigate the effect of the model with two types of controls: fixed-time dependent impulsive control and state-dependent impulsive control. We first consider two treatment (or injection) period scenarios (the periodsT
=
1.665,2.96) for this model with fixed-time dependent impulsive control and p
=
0.1, q
=
0.2, τ
=
0.22, E
(
0
)
=
1.58, and T
(
0
)
=
0.56 as other parameters are fixed in Table 1. When period T
=
1.665, the tumor cells are controlled in lower than the hazardous threshold value H
=
0.7 which is shown in Figure 5 by blue line. However, if period T
=
2.96, the tumor control fails. That is, the number of tumor cells is higher than the hazardous threshold value H
=
0.7 which is shown in Figure 5 by red line. On the other hand, model (4) with state-dependent impulsive control with threshold value H
=
0.7 has approximately T
=
3.0 as treatment period and we can control the number of tumor cells in the prescribed fixed boundary which is shown in Figure 5 by black line. As we know, the proper treatment schedule is very important to control tumor cells because of the costs for treatment.Figure 5
The comparison of the fixed-time impulsive effect and the state impulsive effect (p
=
0.1, q
=
0.2, τ
=
0.22, E
(
0
)
=
1.58, and T
(
0
)
=
0.56).


## 5. Conclusion

In this paper, a mathematical model of tumor and immune cells with state-dependent impulsive chemotherapy and immunotherapy was considered. This is the first attempt to investigate a tumor-immune model with state-dependent impulsive treatments. The criterion for the existence and stability of the positive periodic solution of tumor-immune model, based on Poincaré map and analogue of Poincaré criterion, was presented. We proved existence of positive order-1 periodic solutions and obtained condition of orbitally asymptotical stability of positive order-1 periodic solution for our model. From the results, we were able to see that local asymptotic orbital stability region is large as treatment strengthp of chemotherapy is increasing. Also, increasing strength of chemotherapy and immunotherapy and decreasing side effects of chemotherapy lengthened tumor recurrence time. The simulation results indicate that in determining optimal treatment timing the model with state-dependent impulsive control is more efficient than fixed-time impulsive control.



---

*Source: 2979414-2016-09-20.xml*