# The Dynamical Behaviors for a Class of Immunogenic Tumor Model with Delay

Authors: Ping Bi; Zijian Liu; Mutei Damaris Muthoni; Jianhua Pang
Journal: Computational and Mathematical Methods in Medicine
 (2017)

---

## Abstract

This paper aims at studying the model proposed by Kuznetsov and Taylor in 1994. Inspired by Mayer et al., time delay is introduced in the general model. The dynamic behaviors of this model are studied, which include the existence and stability of the equilibria and Hopf bifurcation of the model with discrete delays. The properties of the bifurcated periodic solutions are studied by using the normal form on the center manifold. Numerical examples and simulations are given to illustrate the bifurcation analysis and the obtained results.


---

## Body

## 1. Introduction

For the longest time, malignant tumors have posed a threat to human lives. Effective strategies based on the immune system have been championed to complement traditional methods of cancer treatment, such as surgery, radiation, and chemotherapy. Cancer immunotherapy is the use of the immune system to treat cancer. Immunotherapy is used to provoke the immune system into attacking the tumor cells by using these cancer antigens as targets. Cell types that can be used in this way are natural killer cells, lymphokine-activated killer cells, cytotoxic T cells, and dendritic cells. There has been much interest in mathematical models describing the interaction between tumor cells and effector cells of the immune system; see [1–6]. An ideal model can provide insight into the dynamics of interactions of the immune response with the tumor and may play an important role in understanding of cancer and developing effective drug therapies. Developing ideal models for such complex processes is not easy. Simple models which display some useful immunological phenomena have been proposed and studied extensively. See Bi and Xiao [7], Galach [4], and Yafia [8–10], and the references cited therein. In 1994, Kuznetsov et al. [11] introduced a model, which describes competition between the tumor and immune cells. It also describes the response of effector cells (ECs) to the growth of tumor cells (TCs). It is assumed that the tumor cells are “immunogenic” and thus subject to immune attack by cytotoxic (CT) or natural killer (NK) cells. This model studies the infiltration of TCs by ECs and also the possibility of ECs inactivation.It is assumed that interactions between ECs and TCs in vitro can be described by the following kinetic scheme describing interactions between ECs and TCs:(1)E+T⇄nnk-1nnnnk1nnC⟶k2E∗+T⟶k3T∗+Ewhere E,T,C,T∗,E∗ are the local concentrations of ECs, TCs, EC-TC complexes, inactivated ECs, and TCs, respectively. The parameters k1,k-1,k2 and k3 are nonnegative constants, which describe the conversion rates of differential cells. Then Kuznetsov and Taylor’s model is as follows:(2)dEdt=s+FC,T-d1E-k1ET+k-1+k3C,dTdt=aT1-bTtot-k1ET+k-1+k2C,dCdt=k1ET-k-1+k2+k3C,dT∗dt=k3C-d2T∗,dE∗dt=k2C-d3E∗,where C is the normal rate of the flow of adult ECs into the tumor site, F(C,T) describes the accumulation of effector cells in the tumor cells localization due to the presence of the tumor. F(C,T)=fC/(g+T), (f,g are constants). If Ttot≈T(Ttot=T+C), it is reasonable to assume dC/dt≈0, that is, C≈kET(k=k1/(k-1+k2+k3)), then F(C,T)=F(E,T). Then we only need to analyze the dynamical behavior of the first two equations.In2003, Galach [4] suggested that the function F is in the Michaelis-Menten form F(E,T)=θET, and thus (2) can be simplified as(3)dxdt=c+θxy-m1xy-d1x,dydt=ay1-by-m2xy,where x is the local concentrations of ECs, and y is the local concentrations of TCs, m1=kk2, m2=kk3, and all the parameters are positive. The properties of the model (3) is studied in [7–9]. The dynamical behaviors and the bifurcation of the model with delay are also studied by [7, 10].In this paper, we analyze the dynamics of an immune response function with Michaelis-Menten formF(E,T)=pET/(g+T), where p and g are positive constants, that is,(4)dEdt=s+pETg+T-mET-dE,dTdt=aT1-bT-nET.In order to simplify the original model, we nondimensionalize (4) by choosing scale for E and T cell population, respectively. Let τ=nT0t, and replace τ with t; thus the model can be written as follows:(5)dxdt=σ+αxy1+y-μxy-δx,dydt=fy1-βy-xy,where y=T/g, x=E/g, τ=tng,σ=s/ng2,α=p/ng,μ=m/n,δ=d/ng,f=a/ng,β=bg and all the above parameters are positive.This paper is organized as follows: In Section2, the model without delay is considered, and the conditions for existence and stability of equilibria are given. In Section 3, the model with delay is studied. The existence of Hopf bifurcation is obtained; the direction and stability of bifurcated periodic solutions are also given with the help of center manifold and bifurcation theories. At the end of this paper, numerical results are given to illustrate the main results of this paper.


## 2. Existence and Stability of the Equilibria

It is easy to see that (5) has a tumor-free equilibrium P0=(σ/δ,0). In order to find the positive equilibria P1=(x∗,y∗), we need to solve the following equations:(6)dxdt=σ+αxy1+y-μxy-δx=0,dydt=fy1-βy-xy=0.Equation (6) is reduced to the following cubic polynomial:(7)a1y3+a2y2+a3y+a4=0,where a1=fβμ, a2=f(-αβ+βδ-μ+βμ), a3=fα-fδ+fβδ-fμ+σ, a4=-fδ+σ. Obviously, system (5) in R2+ is as “well behaved” just as in biology. We have the following lemma by qualitative analysis.Lemma 1.
Ifμ>α, then the solutions of (5) are invariable in D: (8)D=x,y∣0≤x≤2σδ,0≤y≤2β.Proof.
Ifx=0 and 0<y<2/β, then (9)dxdtx=0=σ>0,dydtx=0=fy1-βy>0for0<y<1β,<0for1β<y<2β.Similarly, if y=0 and 0<x<2σ/δ, one has (10)dxdty=0=σ-δx>0,for0<x<σδ,<0,forσδ<x<2σδ,dydty=0=0.Also, if y=2/β and 0<x<2σ/δ, we have (11)dydty=2/β,0<x<2σ/δ=yf1-βy-xy=2/β,0<x<2σ/δ<0,dxdty=2/β,0<x<2σ/δ=σ+x2αβ-β+22μ+δβββ+2y=2/β,0<x<2σ/δ<0,and, then, there exists a positive number x∗, which is given by (12)x∗=σββ+2β+22μ+δβ-2αβ.Such that dx/dt<0 when x>x∗. Similarly, dx/dt>0 when x<x∗. On the other hand, one has (13)0<σββ+2β+22μ+δβ-2αβ<2σδwith the help of μ>α. That is 0<x∗<2σ/δ.
Whenx=2σ/δ and 0<y<2/β, we can prove the result in a similar way. This finishes the proof.It is easy to see that (5) has one tumor-free equilibrium P0=(σ/δ,0). With the distribution of the eigenvalues, we can easily obtain the following results.Theorem 2.
Iff>σ/δ, then P0 is unstable. If f=σ/δ, P0 is stable. If f<σ/δ, P0 is asymptotically stable.Theorem 3.
If>fδ, δβ>μ>α, then P0 is globally asymptotically stable.Proof.
When>δf, βδ>μ>α, then it is easy to see that (14)a1>0,a4=-fδ+σ>0,a3=fα+σ-fδ+fβδ-μ>0,a2=fβμ-α+βδ-μ>0.This implies that P0 is the only positive equilibrium.
Notingσ>df, we know that P0 is locally stable with the help of Theorem 2. In the remaining part of the proof, we only need to prove the global stability. By Lemma 1, we know that the following domain D is invariable: (15)D=x,y∣0≤x≤2σδ,0≤y≤2β.Then it is easy to prove that dx/dt<0 as x>2σ/δ. In addition, when y0>2/β, (16)dydt=fy1-βy-x=-fy1-βy+xy<0.This shows that the vector fields are moving towards D as t increases.
Let Dulac functionB=1/y. Then (17)∂BP∂x+∂BQ∂y≤α-μ-fβ<0.Hence there are no closed trajectory surrounds P0 in field D. That is, the result follows.Using the original parameters, we can give the results as follows.Theorem 4.
Ifs>da/n, db>m>p/g, system (5) has only one critical equilibrium P0=(σ/δ,0). Furthermore, P0 is globally asymptotically stable.Remark 5.
Theorem4 is an instructive results to kill the tumor cells. The tumor cells will be killed out by the immune cells sooner or later under the above conditions; then we only need to take the necessary measures to control the parameters to satisfy the inequity in Theorem 4.In the following, we will study the existence and stability of the tumor-present equilibriumP1(x∗,y∗). Similar to the proof of Lemma  2.4 in [12], we can easily obtain the following results.Theorem 6.
For the number of positive equilibria, we can get the following results:(1)
If>δ>μ/β>μ, σ<fδ, and Δ<0, then (5) has three distinct positive roots.(2)
If one of the following conditions is satisfied, then (5) has two distinct positive roots:
(a)
α>δ>μ/β>μ, σ>fδ, and Δ<0.(b)
α>δ>μ/β>μ, σ>fδ and Δ=0.(3)
Assume one of the following conditions is satisfied, then (5) has one positive root:
(a)
α>δ>μ/β>μ, σ<fδ, and Δ<0.(b)
α>δ>μ/β>μ, σ<fδ, and A=B=0.(c)
δ>α, β>1, α+βδ<μ, σ<δf, and Δ>0.(d)
α>δ>μ/β>μ, σ<fδ, and Δ=0.The corresponding Jacobian matrix atP1(x∗,y∗) is(18)J=-σx∗αx∗1+y∗2-μx∗-y∗-fβy∗.Thus, we can give the following results.Theorem 7.
Ifσβ>μf, then P1 is stable. And no Hopf bifurcation appears around the equilibrium P1.Proof.
It is easy to seeTr(J(P1))=-(σ/x∗+fβy∗)<0. Noting σβ>μf, it is easy to see that (19)det⁡JP1≥y∗x∗1+y∗2σβf-μ1+y∗2+α>0.Thus the results are proved.Theorem 8.
Let(x∗,y∗) be the coordinate of the positive equilibrium P1. Then the following results hold: (1)
Iffσβ≥μx∗2, then P1 is stable.(2)
Iffσβ<μx∗2, then P1 is stable for y∗>αx∗2/(μx∗2-fσβ)-1 and unstable for 0<y∗<αx∗2/(μx∗2-fσβ)-1.Proof.
It is easy to see(20)det⁡JP1=y∗x∗1+y∗2ay∗2+2ay∗+a+α, where a=1/x∗2(σfβ-μx∗2).
Let(21)hy∗=ay∗2+2ay∗+a+α.Then the sign of the det⁡(J(P1)) is the same as that of h(y∗). Let Δ=4α(μ-σfβ/x∗2). If fσβ=μx∗2 for any y∗>0, then a=Δ=0. Hence, (22)hy∗=ay∗2+2ay∗+a+α=α>0.Then the first result can be obtained easily.
Iffσβ>μx∗2, that is, a>0, Δ<0, then h(y∗)>0 for any y∗>0. That is to say, P1 is stable as fσβ>μx∗2. On the other hand, one knows a<0 as fσβ<μx∗2. Noting α>μ, then (23)a+α=1x∗2σfβ-μx∗2+αx∗2>0.Then h(y∗)=0 has roots -αx∗2/(μx∗2-fσβ)-1<0 and αx∗2/(μx∗2-fσβ)-1>0. From y∗>0, we have h(y∗)<0 as y∗>αx∗2/(μx∗2-fσβ)-1 and h(y∗)>0 as 0<y∗<αx∗2/(μx∗2-fσβ)-1. That is to say, the second result holds.In the following, we give some simulation results of the above results. We consider the system (5) and the parameters suggested by V. A. Kuznetsov et al. Then system (5) becomes(24)dxdt=0.1181+1.131xtyt1+yt-0.00311xtyt-0.3743xt,dydt=1.636yt1-0.002yt-xtyt.Obviously, system (24) has three positive equilibria: E1(1.6348,0.366269), E2(0.6538,300.188), and E3(0.1906,441.757). By simple computation, it is easy to know that the eigenvalues of E1 are -0.0367189±0.599721i, and E2 has eigenvalues -1.45811 and 0.295245, and E3 has eigenvalues -1.68973 and -0.375447. These results are represented in Figures 1 and 2.Figure 1
(a) Stable equilibriumE1. (b) Saddle-node equilibrium E2.
(a)
(b)Figure 2
(a) Stable focus equilibriumE3. (b) Saddle node (E2) and stable focus (E3).
(a)
(b)


## 3. Dynamical Behaviors of the Model with Delay

In this section, we shall consider model (5) with delay(25)dxdt=σ+αxy1+y-μxt-τyt-τ-δx,dydt=fy1-βy-xy.The existence of the equilibria is the same as those of (5). The dynamical behaviors of the trivial equilibria and semitrivial equilibria are not difficult; we only study the positive equilibrium P1=(x∗,y∗) here. The characteristic equation of linearized system for (25) at P1 takes the following form:(26)λ2+Bλ+C+Dλ+Ee-λτ=0,where (27)B=fβy∗-αy∗1+y∗+δ=fβy∗+fβy∗2+δ+y∗δ-α1+y∗,C=fβy∗2δ-α+fβy∗2δ-α+fβδ+αx∗1+y∗2,D=μy∗,E=μy∗2fβ-μx∗y∗=μfy∗2βy∗-1. Similar to the proof in [7], then the following results can be obtained easily.Theorem 9.
(1) If(28)B+D>0,C+E>0,D2-B2+2C<0,C2-E2>0orD2-B2+2C2<4C2-E2, then all roots of (26) have negative real parts for all τ≥0.
(2) IfC2-E2<0 or D2-B2+2C>0 and (D2-B2+2C)2=4(C2-E2), then (26) has a pair of purely imaginary roots ±iω+ at τ=τj+.
(3) IfD2-B2+2C>0,C2-E2>0 and D2-B2+2C2>4(C2-E2), then (26) has a pair of purely imaginary roots ±iω+(±iω-, resp.) at τ=τj+(τ=τj-,resp.), where (29)ω2±=12D2-B2+2C±D2-B2+2C2-4C2-E2,τj±=1ω±2jπ+arccosω2-CE-DBω2E2+D2ω2,Dω2-CD+BE≥0,j=1,2,…,1ω±2i+1π-arccosω2-CE-DBω2E2+D2ω2,Dω2-CD+BE<0,i=1,2,….Theorem 10.
(1) If the conditions of (1) in Theorem9 hold, then P1 is asymptotically stable for any τ>0.
(2) IfC2<E2 or the conditions of (3) in Theorem 9 hold, then E∗ undergoes Hopf bifurcation as τ=τ0.
(3) IfD2-B2+2C>0, (D2-B2+2C)2=4(C2-E2) and B2+2DτE+τ2>D2+E2, then E∗ undergoes a Hopf bifurcation as τ=τ0.Proof.
From the analysis of the above, we only need to compute(30)SigndReλdτλ=iω=Sign±D2-B2+2C2-4C2-E2.
In third case,dReλ/dττ=τ0=0, then we need the sign of the second derivative of Reλ(τ) of (26) at the point where λ is passing through λ0.
We can easily obtain(31)SignReλ′′λ=iω=ω2E2+D2ω22E2+D2ω22ω2-2C-2DτE-τ2-2E2A-A2A3with (32)A=-Bω2-BC-DE+τωE2+D2ω2.
Sincesgn⁡Redλ/dτ-1|λ=iω=0, sgn⁡Imdλ/dτ-1|λ=iω≠0. Thus it is easy to see that E2+D2ω2(2ω2-2C-2DτE-τ2)-2E2A-A2<0. This finishes the proof.In the above, we obtained the conditions under which a family of periodic solutions bifurcated from the positive equilibrium atτ=τ0. In the following, we derive the explicit formulae for determining the properties of the Hopf bifurcated solution by using the normal form and the center manifold theory. Throughout this section, we always assume that the system (25) undergoes Hopf bifurcations at the positive equilibrium P1(x∗,y∗) for τ=τ0 and ±iω are the corresponding pure imaginary roots.Lett=τt′, x(t)=x(τt′). Then (25) becomes(33)dxdt=τσ+αxy1+y-μxt-1yt-1-δx,dydt=τfy1-βy-xy.Set τ=τ0+p, p∈R. Then p=0 is a Hopf bifurcation value for (33). Equation (5) can be written as(34)x′t=Lpxt+Fp,xt,where x(t)=(x1(t), x2(t))T∈R2, and Lp:C→R, F:R×C→R are given, respectively, by(35)Lpϕ=A1ϕ10ϕ20+A2ϕ1-1ϕ2-1,Fp,ϕ=τ0+p-αx∗1+y∗3ϕ220+α21+y∗2ϕ10ϕ20-μϕ1-1ϕ2-1-fβϕ220-ϕ10ϕ20,with ϕ(0)=(ϕ1(0),ϕ2(0))T∈C and (36)A1=τ0+pαy∗1+y∗-δαx∗1+y∗2-y∗-fβy∗,A2=τ0+p-μy∗-μx∗00.By Riesz representation theorem, there exists a functionη(θ,p) of bounded variation for θ∈[-1,0] such that(37)Lpθ=∫-10dηθ,pϕθ.For ϕ∈C1([-1,0],R3), define(38)Apϕ=dϕθdθ,θ∈-1,0,∫-10dηs,pϕs,θ=0,Rpϕ=0,θ∈-1,0,Fp,ϕ,θ=0.Then system (34) is equivalent to(39)xt′=Apxt+Rpxt,where xt(θ)=x(t+θ) for θ∈[-1,0]. For Ψ∈C1([0,1],(R3)∗), define(40)A∗ψs=dψsds,s∈0,1,∫-10dηTt,0ψ-t,s=0,and a bilinear inner product(41)ψs,ϕθ=ψ¯0ϕ0-∫-10∫ξ=0θψ¯ξ-θdηθϕξdξ,where η(θ)=η(θ,0). Then A(0) and A∗ are joint operators.Suppose thatq(θ)=(1,ρ)Teiθω0τ0 is the eigenvector of A(θ) corresponding to iω0τ0, with ρ=-y∗/(fβy∗+iω0τ0). Then q∗(s)=D(1,ρ∗)eisω0τ0 is the eigenvector of A∗ with ρ∗=(1/y∗)[αy∗/(1+y∗)-δ-μy∗eiω0τ0-iω0τ0], D=1/(1+ρ¯ρ∗+τ-μy∗-μx∗ρ¯e-iω0τ0) such that 〈q∗(s),q(θ)〉=1, 〈q∗(s),q¯(θ)〉=0.In the following, we use the ideas in Adam and Bellomo [1] to compute the coordinates describing center manifold C0 at p=0. Define(42)zt=q∗,xt,Wt,θ=xtθ-2Reztqθ.On the center manifold C0, we have(43)Wt,θ=WZt,z¯t,θ=W20θz22+W11θzz¯+W02θz¯22+⋯,where z and z¯ are local coordinates for the center manifold C0 in the direction of q∗ and q¯∗. Note that W is real if xt is real. For xt∈C0 of (39), we have(44)zt=iω0τ0z+q¯∗0F0,Wz,z¯+2Rezqθ=iω0τ0z+gz,z¯,where(45)gz,z¯=q¯∗0F0Zt,z¯t=g20θz22+g11θzz¯+g02θz¯22+g21θz2z¯2+⋯.From (43) and (44), we have (46)xt=x1tθ,x2tθ=Wt,θ+zqθ+zqθ¯and q(θ)=(1,ρ)Teiθω0τ0. Comparing the coefficients of the above equality with (45), we obtain (47)g20=2τD¯αx∗1+y∗3ρ2+α21+y∗2ρ-μρe-2iω0τ0-fβρ2-ρ,g11=2τD¯αx∗1+y∗3ρρ¯+α21+y∗22Reρ-μ2Reρ-fβρρ¯-2Reρ,g02=2τD¯αx∗1+y∗3ρ¯2+α21+y∗2ρ¯-μρ¯e2iω0τ0-fβρ¯2-ρ¯+αx∗1+y∗3W1120ρ+z2z¯α21+y∗2W2010ρ¯+W2020+W1110ρ+W1120+αx∗1+y∗3W2020ρ¯,g21=2τD¯αx∗W2020ρ¯+W1120ρ1+y∗3+αW2010ρ¯+W2020+W1110ρ+W112021+y∗2+μW201-1ρe-2iω0τ0+W202-1e-2iω0τ0+W111-1ρe-2iω0τ0+W112-1e-2iω0τ0-fβW2020ρ¯+W1120ρ-W2010ρ¯+W2020+W1110ρ+W1120. In order to determine g21, we also need to compute W20(θ) and W11(θ) as follows:(48)W˙=x˙t-z˙q-zq¯˙=AW-2Req∗¯0F0qθ,θ∈-1,0,AW-2Req∗¯0F0qθ+F0,θ=0,(49)=AW+Hz,z¯,θ,where(50)Hz,z¯,θ=H20θz22+H11θzz¯+H02θz¯22+⋯.On the center manifoldC0 near the origin,(51)W˙=Wzz˙+Wz¯z¯˙.Comparing the coefficients of z, we obtain(52)A-2iωτW20θ=-H20θ,AW11=-H11θ,H20θ=-g20qθ-g¯02θ¯,H11θ=-g11qθ-g¯11q¯θ.From (52) and the definition of A, one has(53)W˙20θ=2iω0τ0W20θ+g20q0eiω0τ0θ+g¯02q¯0e-iω0τ0θ.Solving the above ODE, then(54)W20θ=-ig20ω0τ0q0eiω0τ0θ-ig¯023ω0τ0q¯0e-iω0τ0θ+E1e2iω0τ0θ,where E1=(E1(1),E1(2))∈R2 is a constant vector. Similarly, we obtain(55)W˙11=ig11ω0τ0q0eiω0τ0θ-ig¯11ω0τ0q¯0e-iω0τ0θ+E2,where E2=(E2(1),E2(2))∈R2 is a constant vector.From the definition ofA and the above equations, we can obtain (56)E1=2Aαx∗1+y∗3ρ2+α21+y∗2ρ-μρe-2iω0τ00-fβρ¯2-ρ¯, with (57)A=2iω0+αy∗1+y∗-d-μy∗e2iω0τ0αx∗1+y∗2-μxe2iω0τ0-y∗2iω0-fβy∗-1. In the same way, it is easy to get (58)E2=2αy∗1+y∗-d-μy∗αx∗1+y∗2-μx-y∗2iω0-fβy∗-1αx∗1+y∗3ρρ¯+α21+y∗22Reρ-μ2Reρfβρ¯2-ρ¯. In order to determine E1 and E2, we can calculate W20(θ) and W11(θ) from (54) and (55). Furthermore, we can also determine g21. Thus we can be able to compute the following values:(59)c10=i2ω0τ0g20g11-2g112-g0223+g212,μ2=-Rec10Reλ′τ0,β2=2Rec10,T2=-Imc10+μ2Imλ′τ0ω0τ0.Theorem 11.
Equation (25) undergoes Hopf bifurcation as τ=τ0. If μ2<0 (μ2>0), then the bifurcated periodic solution is supercritical (subcritical); if β<0 (β<0), then the bifurcating solution is stable (unstable). The period of the bifurcated periodic solution is T2.In the following, we present some numerical results of system (25) at different values of τ. We choose the parameters the same as before; the system (25) becomes(60)dxdt=0.1181+1.131xy1+y-0.00311xt-τyt-τ-0.3743x,dydt=1.636y1-0.002y-xy.Obviously, system (60) has three positive equilibria: E1(1.6348,0.366269), E2(0.6538,300.188), and E3(0.1906,441.757). Equilibrium E1 is stable for all τ, while equilibrium E2 is unstable saddle node for all τ. When τ=0, the positive equilibrium E3 is a stable focus, as τ is growing to τ0=0.707732, λ′(τ0)=0.604929-0.176726i. When τ passes the critical value τ0, E3 loses its stability and a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcate from E1. This is clearly demonstrated by Figure 3.Figure 3
(a) Stable equilibriumE1 as τ=1. (b) Bifurcated solution around E3 as τ=0.7.
(a)
(b)


## 4. Conclusions

In 1994, Kuznetsov et al. [11] introduced a model, which describes competition between the tumor and immune cells. This mode describes the response of effector cells to the growth of tumor cells. Later, the model is simplified to a two-dimensional differential equation with Lotka-Volterra form immune response function F(E,T)=θET; the rich dynamical behaviors are studied such as stability, qualitative behaviors, and bifurcation behaviors.In this paper, we studied the model with an immune response function with Michaelis-Menten formF(E,T)=pET/(g+T). In this case, the properties of system (4) are more complicated. There are so many tedious computations even for the existence of the equilibria. The dynamical behaviors of the model are more rich.We have studied the nonlinear dynamics of a two-dimensional general differential system. We first provided linear analysis of the system at the possible equilibria, namely, the semitrivial and positive equilibria, and discussed the existence of Hopf bifurcation at the positive equilibrium. Then we consider the system with delay; we investigated the Hopf bifurcation of the system. The existence and stability of periodic solutions were given. Numerical simulations were presented to illustrate the theoretical analysis and results. Cancer immunosurveillance functions are taken as an important defense to cancer; this is the elimination process. In fact, the existence and stability of the semitrivial equilibrium correspond to the elimination process. Our results on the existence and stability of the Hopf bifurcated periodic solutions describe the equilibrium process. If a stable periodic orbit exists, then the tumor and the immune system can coexist for a long time, although the cancer cannot be eliminated. Furthermore, the parameters are important in controlling the development and progression of the tumor, which is decided by the conditions of the existence of the bifurcations.The existence of oscillatory solutions in the tummor and immune system interaction models demonstrates that the phenomenon has been observed in some related models [7, 10, 13–15]. The initial values and delay are also important in the oscillatory coexistence of the tumor cells and the effector cells. Numerical simulations indicate this information well.



---

*Source: 1642976-2017-10-25.xml*