# American Institute of Mathematical Sciences

## Monotonic and nonmonotonic immune responses in viral infection systems

 1 School of Mathematics and Statistics, Bioinformatics Center of Henan University, Kaifeng 475001, Henan, China 2 Department of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, Jiangsu, China 3 Department of Mathematics, Wilfrid Laurier University, Waterloo, Ontario, N2L 3C5, Canada

* Corresponding author: Shaoli Wang

Received  September 2019 Revised  November 2020 Published  January 2021

Fund Project: This work is supported by Science and Foundation of Technology Department of Henan Province (No.192102310089), Foundation of Henan Educational Committee (No.19A110009), Natural Science Foundation of Henan (No. 202300410045) and Grant of Bioinformatics Center of Henan University (No. 2019YLXKJC02)

In this paper, we study two-dimensional, three-dimensional monotonic and nonmonotonic immune responses in viral infection systems. Our results show that the viral infection systems with monotonic immune response has no bistability appear. However, the systems with nonmonotonic immune response has bistability appear under some conditions. For immune intensity, we got two important thresholds, post-treatment control threshold and elite control threshold. When immune intensity is less than post-treatment control threshold, the virus will be rebound. The virus will be under control when immune intensity is larger than elite control threshold. While between the two thresholds is a bistable interval. When immune intensity is in the bistable interval, the system can have bistability appear. Select the rate of immune cells stimulated by the viruses as a bifurcation parameter for nonmonotonic immune responses, we prove that the system exhibits saddle-node bifurcation and transcritical bifurcation.

Citation: Shaoli Wang, Huixia Li, Fei Xu. Monotonic and nonmonotonic immune responses in viral infection systems. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021035
##### References:

show all references

##### References:
Bifurcation diagram of system (3). The solid line represents the stable equilibrium of infected CD4+ T cells and the dashed line represents the unstable equilibrium of infected CD4+ T cells. The post-treatment control threshold is $c_{2} = 0.2500$, the elite control threshold is $c^{**}_1\approx0.6505$ and the bistable interval is $(0.2500, 0.6505).$ Here, $c = 0.37\; \; \mbox{day}^{-1}$ and the values of other parameters are listed in (4)
Time histories and trajectories of system (3) with different initial conditions. The system has a stable equilibria $E^{(2)}_{1}$. Here, $c = 0.2\; \; \mbox{day}^{-1}$ is less than the post-treatment control threshold $P_I$ and other parameter values are listed in (4)
Time histories and trajectories of system (3) with different initial conditions. Here, system (3) has two different stable equilibria $E^{(2)}_{1}$ and $E_{-}^{2*}$ with $c = 0.37\; \; \mbox{day}^{-1}$. Other parameter values are listed in (4)
Time histories and trajectories of system (3) with different initial conditions. System (3) only has the positive equilibrium $E_{-}^{2*}$, which is stable with $c = 0.65\; \; \mbox{day}^{-1}$. Other parameter values are listed in (4)
Bifurcation diagram of system (6). The solid line is the stable equilibrium and the dashed line denotes the unstable equilibrium. The post-treatment control threshold is $c_{2} = 2.5000$, the elite control threshold is $c^{**}_2\approx3.5278$ and the bistable interval is $(2.5000, 3.5278).$ Here, $c = 3\; \; \mbox{day}^{-1}$ and other parameter values are listed in (7)
Time histories and phase portraits of system (6). System (6) has two different stable equilibria $E^{(4)}_{1}$ and $E_{*}^{4-}$. Here, $c = 3\; \; \mbox{day}^{-1}$ and other parameter values are listed in (7). We choose different initial values
Phase portraits of system (6). (A) Choosing $c = 2\; \; \mbox{day}^{-1}$, which is less than the post-treatment control threshold $c_{2} = 2.5000$, system (6) only has a stable equilibrium $E^{(4)}_{1}$; (B) Choosing $c = 4\; \; \mbox{day}^{-1}$, which is greater than the elite control threshold $c^{**}_2\approx3.5278$, system (6) only has the stable equilibria $E_{*}^{4-}$. Other parameter values are listed in (7)
The stabilities of the equilibria and the behaviors of system (3) in the case $1<\mathcal {R}^{(2)}_{0}<\mathcal {R}^{(1)}_{c}$
 $E^{(2)}_{0}$ $E^{(2)}_{1}$ $E_{*}^{2-}$ $E_{*}^{2+}$ System (3) $R^{(2)}_{0}<1$ GAS — — — Converges to $E^{(2)}_{0}$ $1 $ E^{(2)}_{0}  E^{(2)}_{1}  E_{*}^{2-}  E_{*}^{2+} $System (3)$ R^{(2)}_{0}<1 $GAS — — — Converges to$ E^{(2)}_{0}  1
The stabilities of the equilibria and the behaviors of system (3) in the case $\mathcal {R}^{(2)}_{0}>\mathcal {R}^{(1)}_{c}$
 $E^{(2)}_{0}$ $E^{(2)}_{1}$ $E_{*}^{2-}$ $E_{*}^{2+}$ System (3) $R^{(2)}_{0}<1$ GAS — — — Converges to $E^{(2)}_{0}$ $R^{(2)}_{0}>1$, $0R^{(1)}_{c}>1,$ $c_{2}R^{(1)}_{c}>1,$ $c>c^{**}_1$ US US LAS — Converges to $E_{*}^{2-}$
 $E^{(2)}_{0}$ $E^{(2)}_{1}$ $E_{*}^{2-}$ $E_{*}^{2+}$ System (3) $R^{(2)}_{0}<1$ GAS — — — Converges to $E^{(2)}_{0}$ $R^{(2)}_{0}>1$, $0R^{(1)}_{c}>1,$ $c_{2}R^{(1)}_{c}>1,$ $c>c^{**}_1$ US US LAS — Converges to $E_{*}^{2-}$
The stabilities of the equilibria and the behaviors of system (6) in the case $1<\mathcal {R}^{(4)}_{0}<\mathcal {R}^{(2)}_{c}$
 $E^{(4)}_{0}$ $E^{(4)}_{1}$ $E_{*}^{4-}$ $E_{*}^{4+}$ System (6) $R^{(4)}_{0}<1$ GAS — — — Converges to $E^{(4)}_{0}$ $1 $ E^{(4)}_{0}  E^{(4)}_{1}  E_{*}^{4-}  E_{*}^{4+} $System (6)$ R^{(4)}_{0}<1 $GAS — — — Converges to$ E^{(4)}_{0}  1
The stabilities of the equilibria and the behaviors of system (6) in the case $\mathcal {R}^{(4)}_{0}>\mathcal {R}^{(2)}_{c}$
 $E^{(4)}_{0}$ $E^{(4)}_{1}$ $E_{*}^{4-}$ $E_{*}^{4+}$ System (6) $R^{(4)}_{0}<1$ GAS — — — Converges to $E^{(4)}_{0}$ $R^{(4)}_{0}>1$, $0R^{(2)}_{c}>1,$ $c_{2}R^{(2)}_{c}>1,$ $c>c^{**}_2$ US US GAS — Converges to $E_{*}^{4-}$
 $E^{(4)}_{0}$ $E^{(4)}_{1}$ $E_{*}^{4-}$ $E_{*}^{4+}$ System (6) $R^{(4)}_{0}<1$ GAS — — — Converges to $E^{(4)}_{0}$ $R^{(4)}_{0}>1$, $0R^{(2)}_{c}>1,$ $c_{2}R^{(2)}_{c}>1,$ $c>c^{**}_2$ US US GAS — Converges to $E_{*}^{4-}$
 [1] Rui Hu, Yuan Yuan. Stability, bifurcation analysis in a neural network model with delay and diffusion. Conference Publications, 2009, 2009 (Special) : 367-376. doi: 10.3934/proc.2009.2009.367 [2] Yves Dumont, Frederic Chiroleu. Vector control for the Chikungunya disease. Mathematical Biosciences & Engineering, 2010, 7 (2) : 313-345. doi: 10.3934/mbe.2010.7.313 [3] J. Frédéric Bonnans, Justina Gianatti, Francisco J. Silva. On the convergence of the Sakawa-Shindo algorithm in stochastic control. Mathematical Control & Related Fields, 2016, 6 (3) : 391-406. doi: 10.3934/mcrf.2016008 [4] Diana Keller. Optimal control of a linear stochastic Schrödinger equation. Conference Publications, 2013, 2013 (special) : 437-446. doi: 10.3934/proc.2013.2013.437 [5] Alberto Bressan, Ke Han, Franco Rampazzo. On the control of non holonomic systems by active constraints. Discrete & Continuous Dynamical Systems - A, 2013, 33 (8) : 3329-3353. doi: 10.3934/dcds.2013.33.3329 [6] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [7] Guirong Jiang, Qishao Lu. The dynamics of a Prey-Predator model with impulsive state feedback control. Discrete & Continuous Dynamical Systems - B, 2006, 6 (6) : 1301-1320. doi: 10.3934/dcdsb.2006.6.1301 [8] A. K. Misra, Anupama Sharma, Jia Li. A mathematical model for control of vector borne diseases through media campaigns. Discrete & Continuous Dynamical Systems - B, 2013, 18 (7) : 1909-1927. doi: 10.3934/dcdsb.2013.18.1909 [9] Luke Finlay, Vladimir Gaitsgory, Ivan Lebedev. Linear programming solutions of periodic optimization problems: approximation of the optimal control. Journal of Industrial & Management Optimization, 2007, 3 (2) : 399-413. doi: 10.3934/jimo.2007.3.399 [10] Shanjian Tang, Fu Zhang. Path-dependent optimal stochastic control and viscosity solution of associated Bellman equations. Discrete & Continuous Dynamical Systems - A, 2015, 35 (11) : 5521-5553. doi: 10.3934/dcds.2015.35.5521 [11] Peter Benner, Jens Saak, M. Monir Uddin. Balancing based model reduction for structured index-2 unstable descriptor systems with application to flow control. Numerical Algebra, Control & Optimization, 2016, 6 (1) : 1-20. doi: 10.3934/naco.2016.6.1

2019 Impact Factor: 1.27