key: cord-0005299-2558niqn authors: Nie, Lin-Fei; Teng, Zhi-Dong; Jung, Il Hyo title: Complex dynamic behavior in a viral model with state feedback control strategies date: 2014-04-02 journal: Nonlinear Dyn DOI: 10.1007/s11071-014-1372-7 sha: 34261678facd3ba33b2c7c7017d3aa29941de291 doc_id: 5299 cord_uid: 2558niqn With the consideration of mechanism of prevention and control for the spread of viral diseases, in this paper, we propose two novel virus dynamics models where state feedback control strategies are introduced. The first model incorporates the density of infected cells (or free virus) as control threshold value; we analytically show the existence and orbit stability of positive periodic solution. Theoretical results imply that the density of infected cells (or free virus) can be controlled within an adequate level. The other model determines the control strategies by monitoring the density of uninfected cells when it reaches a risk threshold value. We analytically prove the existence and orbit stability of semi-trivial periodic solution, which show that the viral disease dies out. Numerical simulations are carried out to illustrate the main results. Viruses are the most abundant type of biological entities and are found in almost every ecosystem on Earth, which can infect all types of life forms, from animals and plants to bacteria and archaea. Many serious diseases, such as ebola, AIDS, avian influenza, and SARS, are caused by viruses. The control and hence eradication of diseases is one of the major concerns in the study of viral epidemiology. Since viral reproduction always involves host cells and uses the cellular machinery for the synthesis of their genome and other components, there have been two methodologies used to mimic these processes by mathematical models. Perhaps the earliest and simplest classical virus dynamics model was developed by Kermack et al. in [1] and Anderson et al. in [2, 3] as follows where x denotes the density of uninfected cells and y denotes the density of infected cells, λ is the rate of production of uninfected cells, δ is their per capita death rate, β is the rate of infection of uninfected cells, and μ is the rate of disappearance of infected cells. Model (1) is very important in viral epidemiology, which has been studied by many authors. A more detailed description of model (1) and its dynamic behaviors may be found in [1] [2] [3] and the references therein. Not long after that many similar models which include the dynamics of free virus or immune response have been used to describe the short-term dynamics of virus load during drug treatment and have helped to estimate virus turnover rates in vivo. We refer to some of them in [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] and the references therein. However, Bartholdy et al. [14] and Wodarz et al. [15] found that the turnover of free virus is much faster than that of infected cells, which allowed them to make a quasisteady state assumption-that is, the amount of free virus is simply proportional to the number of infected cells. Therefore, the density of infected cells y in model (1) also can be considered as a measure of the free virus load. In recent years, many researchers studied the evolution of virus dynamics models from different perspectives. These include the existence of the threshold value, which is the index for the persistence and extinction of a viral disease; the local or global stability of virus-free equilibrium and the virus equilibrium; the existence and stability of periodic solutions, to name just a few (see, e.g., [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] ). Particularly, Nakata [26] studied the global dynamics of a viral infection model with a latent period and nonlinear function which denotes the incidence rate of the virus infection in vivo. Wang and Zhao [27] considered the investigation of the effects of periodic drug treatment on a standard within-host virus model. Ball et al. [28] examined the dynamics of strain replacement in a simple model that includes a convex trade-off between rapid virus reproduction and long-term host cell survival. Lang and Li [29] investigated the consequences of a more general CTL response and show that a sigmoidal response function gives rise to complex behaviors previously unobserved. It is well known that exploring an effective and easily implementable control measure to keep the level of spread of diseases is significant both theoretically and practically. For example, of many strategies, the efficient ways for elimination or control of malignant neoplasm are still chemotherapy, radiotherapy, and immunotherapy in clinical treatment in recent years. The clinical data show that the course of these treatments for some viral diseases is relatively short but the density of various cells in vivo changed radically, which in nature is submitted to short temporary effects that are negligible compared to the process duration. These short-time perturbations are often assumed to be in the form of impulses in the modeling process. Very recently, many works have been focused on the analysis of mathematical models described by ordinary differential equations with control effects [30] [31] [32] [33] . Particularly, the state-dependent impulsive control strategies are applied widely to the prevention of spread of infectious disease due to its economic, high efficiency, and feasibility nature. This idea can be found in many other areas like agricultural production and fishery industry; where the control measures, such as catching, poisoning, releasing the natural enemy, and harvesting, are taken only when the number of population reaches an economic threshold value. Following this idea, the dynamic behaviors of population, epidemic and language models with state-dependent impulsive effects are considered, and the existence and stability of positive periodic solution by the Poincaré map, properties of the Lambert W function, and analog of Poincaré criterion are obtained [34] [35] [36] [37] [38] . Motivated by the above works, we propose, in this paper, two novel virus dynamics models with state feedback control strategies. The main purpose is to investigate the state feedback control strategies which govern whether the viral disease dies out or not, and further to examine how the state feedback control strategies affect the prevention and control of viral disease. This paper is organized as follows. In the next section, we introduce two novel virus dynamics models where the state feedback control strategies are considered, some basic definitions, and an important lemma. The density of infected cells (or free virus) as control threshold value, and some sufficient conditions are presented in Sect. 3 for the existence and orbital stability of positive periodic solutions. In Sect. 4, the existence and stability of semi-trivial periodic solution are obtained, where the density of uninfected cells is the control threshold value. The numerical simulations are carried out in Sect. 5 for illustration. Some concluding remarks are presented in Sect. 6. The basic reproductive rate for model (1) is R 0 = λβ μδ , which is defined as the average number of secondary infected cells generated by a single infected cell placed in an uninfected cell (or free virus). It is clear that (i) if R 0 < 1, then model (1) admits only a global asymptotical stable virus-free equilibrium ( λ δ , 0); (ii) if R 0 > 1, then model (1) admits an unstable trivial equilibrium ( λ δ , 0) and a unique globally asymptotically stable positive equilibrium (x * , y * ), where Generally, the objective of treatment and control for viral disease is to reduce the density of free virus or infected cells in vivo. For all this, we propose, firstly, a novel virus dynamics model with state feedback control strategies, while regarding the density of free virus or infected cells as control threshold value. The control model is governed by the following state feedback impulsive differential equation: The meaning of model (3) is as follows: when the density of infected cells (or free virus) y reaches the critical threshold value Y c at time t i (Y c ) at the i-th time, control strategies (such as chemotherapy, radiotherapy, immunotherapy, etc.) are taken and the densities of uninfected cells x and infected cells (or free virus) y immediately become (1 − p)x(t i (Y c )) + τ and (1−q)y(t i (Y c )), respectively, where p, q ∈ (0, 1), and τ > 0. Remark 1 It is clear that a priori time of control strategies depends on the density of infected cells (or free virus), which makes our control strategy a "state feedback control." On the other hand, from the dependencies between the densities of uninfected cells and infected cells (or free virus) in model (1) , it is obvious that if the density of uninfected cells is high then infected cells (or free virus) are growing fast. Therefore, regarding the density of uninfected cells as control threshold value, we propose another novel model simulating immune boosting and viral suppressing as state feedback control strategies, and investigate its dynamical behaviors. We suppose that, therefore, when the density of uninfected cells x reaches the risk threshold value X r at time t i (X r ) at the i-th time, the control strategies are taken and the densities of uninfected cells x and infected cells (or free virus) y turn very suddenly to a great degree to (1− p) x(t i (X r ))+τ and (1−q)y(t i (X r )), respectively, where p, q ∈ (0, 1), τ > 0, and (1 − p)X r + τ < X r . Under the assumption aforementioned, we come to a control model that is governed by the following ordinary differential equation with state feedback control strategies: We assume, throughout this paper, that R 0 = λβ μδ > 1. By the biological background, we only consider models (3) and (4) in the region R 2 + = {(x, y) : x > 0, y > 0}, where the biology makes sense. Obviously, R 2 + is divided into four domains with vertical isocline dx dt = 0 and horizontal isocline dy dt = 0, followed by The global existence and uniqueness of solution for models (3) and (4) are guaranteed by the smoothness of the right-hand sides of models (3) and (4) . For more details, we refer to [39] . (3) and (4) with the initial value in the interior of R 2 + at time t = t 0 ≥ 0 are positive. The proof of Lemma 1 is obvious, hence we omit it here. Let S ⊂ R 2 + be a nonempty set and P 0 ∈ R 2 + be a point. The distance between P 0 and S is defined by ρ(P 0 , S) = inf P∈S |P − P 0 |. Let z(t) = (x(t), y(t)) be a solution of model (3) starting from initial point z 0 ∈ R 2 + at t = t 0 . We define the positive orbit as follows: Definition 2 (Orbitally asymptotical stability [40] ) A trajectory O + (z 0 , t 0 ) is said to be orbitally asymptotically stable if it is orbitally stable, and there exists a constant η > 0 such that for any other solution z Next, for model (3), let again. Therefore, we can define a Poincaré map on section Γ 2 as follows: Furthermore, we can define Similarly, we can define two sections to model (4) by and a Poincaré map on section Γ 4 Next, we consider the following autonomous system with pulse effects: where f and g are continuous differentiable functions defined on R 2 and ϕ is a sufficiently smooth function with gradϕ = 0. Let (μ(t), ν(t)) be a positive T -periodic solution of system (8) . The following result comes from Corollary 2 of Theorem 1 of [41] . is the time of the j-th jump, then (μ, ν) is orbitally asymptotically stable. Since virus equilibrium (x * , y * ) is globally asymptotically stable for model (1) , any solution of model (3) without state feedback control strategy will eventually tend to (x * , y * ). If Y c ≤ y * , we easily see that trajectory of model (3) with initial value (x 0 , y 0 ) ∈ Γ 1 will intersect section Γ 2 infinitely many times, and when Y c > y * , then trajectory of model (3) starting from point (x 0 , y 0 ) ∈ Γ 1 may intersect section Γ 2 finitely many times. Therefore, in this section, we give some sufficient conditions for the existence and stability of positive periodic solutions in two cases of Y c ≤ y * and Y c > y * , respectively. The following result is on the existence of positive order-1 periodic solution. Theorem 1 For any p, q ∈ (0, 1), and τ > 0, model (3) admits a positive order-1 periodic solution. (3) will enter into domain III or in succession enter into domain IV and finally intersect Γ 2 at point . This shows that there are two different solutions which start from this point. It contradicts with the uniqueness of solution for model (3) . Therefore, by Poincaré map (6), it follows that On the other hand, for the intersection Namely,x 2 >x 1 for this case. Thus, from the above discussion, we get that (i) ifx 1 =x 2 , then model (3) has a positive order-1 periodic solution; (ii) ifx 2 >x 1 , then By (10) and (11), it follows that Poincaré map (6) has a fixed point. This amounts to saying that model (3) has a positive order-1 periodic solution. This completes the proof. Let (φ(t), ψ(t)) be a positive order-1 periodic solution of model (3) with period T . On the orbital stability of solution (φ(t), ψ(t)) of model (3), we have the following Theorem 2. and Furthermore, it follows from (13), (14) , and (9) that and On the other hand, integrating both sides of the second equation of model (3) along the orbit C + C gives From (15)- (17), we obtain By condition (12), we know that model (3) satisfies all conditions of Lemma 2. It then follows from Lemma 2 that the order-1 periodic solution (φ(t), ψ(t)) of model (3) is orbitally asymptotically stable. This completes the proof. Remark 2 Generally, the condition (12) of Theorem 2 is not easy to test since the expression of periodic solution (φ(t), ψ(t)) is unknown. We note that, however, it is also weaker since the exponent term of condition (12) is less than 1. The following corollary is a direct consequence of Theorem 2. then (φ(t), ψ(t)) is orbitally asymptotically stable. 3.2 Case Y c > y * As already mentioned above, this case Y c > y * is complex since the virus equilibrium (x * , y * ) of model (1) is globally asymptotically stable. Additionally, since control measures are aimed at reducing the density of infected cells (or free virus) y, we suppose that (1 − q)Y c < y * . We now have the following result on the existence and orbital stability of positive orderk periodic solution of model (3) with Y c > y * . Theorem 3 For any p, q ∈ (0, 1), and τ > 0, one of the following statements is valid. has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable. has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable. x * ) will tend to virus equilibrium (x * , y * ) and not intersect with section Γ 2 . Moreover, since (1− p) λ δ +τ 2). On the relation of x 0 , x 1 , and x 2 , there are the following two different cases. (C1) x 0 < x 1 In this case, we can get from (18) that x 1 > x 2 . This results in the relation of x 0 , x 1 , and x 2 to be one of the following two cases. In this case, x 3 > x 1 > x 2 by (18). Repeating the above process, we have (ii) x 0 < x 2 < x 1 Similar to (i), we have In this case, it follows from (18) that x 1 < x 2 . This results in the relation of x 0 , x 1 , and x 2 to be one of two cases. In this case, x 2 > x 1 > x 3 by (18) . Repeating the above process, we have (iv) x 1 < x 2 < x 0 Similar to (iii), we have Further, in (i) of case (C1), we have lim k→∞ x 2k = θ * and lim k→∞ x 2k+1 = θ * , where x * < θ * < θ * < d λ . Hence θ * = F 1 (θ * ) and θ * = F 1 (θ * ). So model (3) has an orbitally asymptotically stable positive order-2 periodic solution. Similarly, in (ii) of case (C1) and (iv) of case (C2), model (3) has an orbitally asymptotically stable positive order-1 periodic solution. In (iii) of case (C2), model (3) has an orbitally asymptotically stable positive order-2 periodic solution. This proves (a). Next, we prove conclusion (b). If for any will intersect with section Γ 2 infinitely many times due to state feedback control strategies and the fact (1 − p) λ δ + τ < x * . Similar to conclusion (a), we can also obtain that model (3) has a positive order-1 or order-2 periodic solution, which is orbitally asymptotically stable. Conclusion (b) thus follows. Finally, we turn to (c). If for any x ∈ (0, x * ], the x * ] will tend to virus equilibrium (x * , y * ) and not intersect with section Γ 2 . Furthermore, any other trajectory intersects section Γ 2 at most finitely many times, and then tends to virus equilibrium (x * , y * ) due to (1 − p) λ δ + τ < x * . In this case, model (3) has no positive order-k (k ≥ 1) periodic solution and virus equilibrium (x * , y * ) is globally asymptotically stable. This is (c). This completes the proof. Remark 3 From conclusion (a) of Theorem 3, we note that (1− p) λ δ 0, and X r ≤ x * , the semi-trivial periodic solution (21) is orbitally asymptotically stable. Proof Suppose that ( x(t), y(t)) is a solution of of small-amplitude perturbation of periodic solution (φ(t), ψ(t)) with initial value ( x(0), y(0)) = ((1 − p)X r + τ, y 0 ), which first intersects section Γ 4 at point (X r , y 1 ) and then jumps to point ((1 − p)X r + τ, y + 1 ). Further, solution ( x(t), y(t)) intersects section Γ 4 at point (X r , y 2 ) again. Repeating the above process, we have two point sequences {((1 − p)X r + τ, y + n )} and {(X r , y n )}, where y + n = (1 − q) y n . Further, by Theorem 4, it is clear that lim t→∞ y(t) = 0. This shows that semi-trivial periodic solution (21) is orbitally asymptotically stable. This completes the proof. Case of x * < X r < λ δ From the properties (5) of the phase space of model (4), we know that there a point E 0 ((1 − p)X r + τ, γ ) such that trajectory O + (E 0 , t 0 ) is tangent to the section Γ 4 at the intersection E 1 (X r , λ−δ X r β X r ) of section Γ and isocline dx dt = 0, where γ = γ (X r , p, τ ). starting from point E(X r , y) (where y ≤ λ−δ X r β X r ) will intersect with Γ 4 infinitely many times due to control strategies x = −pX r + τ and y(t) = −qy(t). For any two points E i (X r , y i ) and E j (X r , y j ) on section Γ 4 , where y i , y j ∈ (0, λ−δ X r β X r ), and y i < y j , in view of control strategies, Similar to the discussion of Theorem 3, we give the following result on the existence and orbital stability of periodic solution. Theorem 6 For any p, q ∈ (0, 1), τ > 0, and X r > x * , if then model (4) admits an order-1 periodic solution which is orbitally asymptotically stable. Remark 5 We want to make clear that the order-1 periodic solution in Theorem 6 may be a positive periodic solution, or it may also be a semi-trivial periodic solution by condition (22) . To illustrate the theoretical results and the feasibility of state control strategies, we perform numerical simulations for different control parameters p, q, τ , Y c , and X r using the software MATLAB. We fixed parameters as those in papers [4, 8, 11] with λ = 50, δ = 0.1, β = 0.002, and μ = 0.5 in model (1) . It is easy to calculate that model (1) has a unique globally asymptotically stable virus equilibrium (x * , y * ) = (250, 50), which is illustrated here by the blue line in Fig. 1a . Firstly, we choose the control parameters to be p = 0.3, q = 0.5, τ = 50, and Y c = 45 < y * , respectively. Using Theorem 1, we know that model (3) has a positive order-1 periodic solution (φ(t), ψ(t)) which is shown by red line in Fig. 1a . At the same time, Fig. 1a also shows that the periodic solution starting from the initial point (φ(0), ψ(0)) = (251.8965, 22.5000). Further, by the condition of Corollary 1, it can be easily shown that Therefore, it follows that (φ(t), ψ(t)) is orbitally asymptotically stable by Corollary 1, which is shown in Fig. 1b . Further, we fixed the values of q, τ , and Y c as mentioned before and choose p to be 0.2, 0.4, 0.6, and 0.8, respectively. Numerical simulations show that the density of y can be controlled within a certain limit and the period T of order-1 periodic solution for model (3) increases with the increase of immune strength p, which is shown in Fig. 2a . The similar numerical results are presented in Fig. 2b . Secondly, let p = 0.3, q = 0.5, τ = 50, and Y c = 55 > y * ; it is easy to see that model (3) has a positive locally orbitally asymptotically stable order-1 periodic solution by parts (a) or (b) of Theorem 3, which is shown in Fig. 3 . However, we only change model parameter Y c = 56, and fixed x(0) = 250 and y(0) ∈ [5, 35] , it can be found that the attraction domains of positive order-1 periodic solution and virus equilibrium appear alternately from Fig. 4 . It is different from the case of Y c = 55. Therefore, the value of control threshold value Y c plays an important role in the treatment of various diseases. Further, we fixed Y c = 55, q = 0.5, and τ = 50 and change the control strength p to be 0. 35 Fig. 5 show that trajectories intersect with section Γ 2 finitely many times and then tend to virus equilibrium (x * , y * ) when p is small, but trajectories intersect with section Γ 2 infinitely many times and model has a orbitally asymptotically stable periodic solution with increase of parameter p. The plots in Fig. 5 show also that control strength p has only a small effect on the periodic of treatment cycle. Similar results can be obtained if we choose q as a control parameter. However, if we fixed Y c = 55, p = 0.3, q = 0.5, and change the control strength τ to be 120, 80, 40, and 0, respectively. As shown in Fig. 6 , the periodic of treatment cycle decrease as control strength τ increases. This implies that when we reduce the concentration Thirdly, we discuss the existence and stability of semi-trivial periodic solution of model (4) . Here, we choose p = 0.3, q = 0.5, τ = 5, and X r = 250 = x * . By Theorem 5, it follows that model (4) has a orbitally asymptotically stable semi-trivial periodic solution, which are shown in Fig. 7 . Theoretical results and numerical simulations show clearly that infected cells (or free virus) y dies out if we control the density of uninfected cells x within a certain range. It also pro- Fig. 6 The effect of control strength τ for the dynamical behaviors of model (3) with Y c = 55, p = 0.3, and q = 0.5 vides theoretical basis for finding a new measure to prevent the spread of virus disease. In fact, this is also consistent with the results in [42] , where Ross proposed a mathematical model to study the spread between human beings and mosquitoes for malaria in earlier 1911. A concept of threshold density is introduced and it is concluded that "· · · in order to counteract malaria anywhere we need not banish Anopheles there entirelywe need only to reduce their numbers below a certain figure. " Finally, fixing p = 0.5, τ = 10, and X r = 350 > x * , model (4) has an orbitally asymptotical stability of positive order-1 periodic solution when q = 0.4385 as shown in Fig. 8 . At the same time, (4) with p = 0.5, q = 0.4385, τ = 10, and X r = 350 > x * many numerical simulations show that the attraction domain of this periodic solution is very small; that is, the stability of this periodic solution is very sensitive to the initial value of solution. The plots in Fig. 9 show that the solution tends to virus equilibrium for q = 0.43, tends to positive order-1 periodic solution for q = 0.4385, and tends to the semi-trivial periodic solution for q = 0.44, respectively. It implies that control strength q is sensitive for the persistence and extinction of infected cells (or free virus) y for X r > x * . This is the relevant contents of Theorem 6 and Remark 5. q=0.43 q=0.4385 q=0.44 Fig. 9 The effect of control strength q for the dynamical behaviors of model (4) with p = 0.5, τ = 10, and X r = 350 > x * The dynamic complexity of two novel virus epidemic models with state feedback control strategies are analyzed systemically in this paper. Firstly, choosing the density of infected cells or free virus as a control threshold value, some sufficient conditions on the existence and orbitally asymptotical stability of positive order-1 or order-2 periodic solution of model (3) are presented by the Poincaré map, the analog of Poincaré criterion and qualitative analysis method. This amounts to that we can control the density of infected cells or free virus at a low level over a long period of time by adjusting control strength. It is concluded that the state feedback control strategies are feasible and effective. Further, choosing the density of uninfected cells x as control threshold value, we obtain that the existence and orbitally asymptotical stability of periodic solution. In particular, the existence and orbitally asymptotical stability of semi-trivial periodic are obtained. This implies that infected cells or free virus y dies out if we control the density of cells x within a certain range. It provides theoretical basis for finding a new measure to prevent and control the spread of viral disease. Contributions to the mathematical theory of epidemics Modern parasitology Population biology of infectious dieases: I Human immunodeficiency virus drug therapy and virus load Virus dynamics and drug therapy Viral dynamics in hepatitis B virus infection Population dynamics of immune responses to persistent viruses A viral infection model with periodic immune response and nonlinear CTL response Virus dynamics: a global analysis HIV-1 infection kinetics in tissue cultures Complex dynamic behavior in a viral model with delayed immune response HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection Persistent virus infection despite chronic cytotoxic Tlymphocyte activation in Gamma interferon-deficient mice infected with lymphocytic chroriomeningitis virus The importance of lytic and nonlytic immune responses in viral infections Effects of replicative fitness on competing HIV strains A model for coupling within-host and betweenhost dynamics in an infectious disease The importance of immune responses in a model of hepatitis B virus Dynamical behaviour of an HIV/AIDS epidemic model Hopf bifurcation analysis in a delayed oncolytic virus dynamics with continuous control A delayed HIV-1 infection model with Beddington-DeAngelis functional response Optimal control of drug therapy: Melding pharmacokinetics with viral dynamics Global stability of a delayed HIV infection model with nonlinear incidence rate A viral infection model with a nonlinear infection rate. Bound. Value Probl Stability and Hopf bifurcation for a virus infection model with delayed humoral immunity response Global dynamics of a viral infection model with a latent periodic and Beddington-DeAngelis response. Nonlinear Anal A within-host virus model with periodic multidrug therapy Modeling withinhost evolution of HIV: mutation, competition and strain replacement Stable and transient periodic oscillations in a mathematical model for CTL response to HTLV-I infection Biological control of a predator-prey system through provision of a super predator A mathematical study of a prey-predator model in relevance to pest control Dynamics in two nonsmooth predator-prey models with threshold harvesting Stability of discrete-time Markovian jump delay systems with delayed impulses and partly unknown transition transition probabilities A state dependent pulse control strategy for a SIRS epidemic system Dynamic analysis of an SIR epidemic model with state dependent pulse vaccination Dynamic analysis of a two-Language competitive model with control strategies The dynamics of a prey-predator model with impulsive state feedback control. Discrete Contin State-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences Theory of Impulsive Differential Equations Dynamics and Bifurcations Orbital stability of periodic solutions of autonomous systems with impulse effect The Prevention of Malaria The authors would like to thank antonymous referees for their constructive suggestions and comments that substantially improved the original manuscript. This research has been partially supported by the National Natural Sci-