Abstract

Two virus infection models with antibody immune response and chronically infected cells are proposed and analyzed. Bilinear incidence rate is considered in the first model, while the incidence rate is given by a saturated functional response in the second one. One main feature of these models is that it includes both short-lived infected cells and chronically infected cells. The chronically infected cells produce much smaller amounts of virus than the short-lived infected cells and die at a much slower rate. Our mathematical analysis establishes that the global dynamics of the two models are determined by two threshold parameters and . By constructing Lyapunov functions and using LaSalle's invariance principle, we have established the global asymptotic stability of all steady states of the models. We have proven that, the uninfected steady state is globally asymptotically stable (GAS) if , the infected steady state without antibody immune response exists and it is GAS if , and the infected steady state with antibody immune response exists and it is GAS if . We check our theorems with numerical simulation in the end.

1. Introduction

In recent years, many mathematical models have been proposed to study the dynamics of viral infections such as the human immunodeficiency virus (HIV), the hepatitis C virus (HCV), and the hepatitis B virus (HBV) (see, e.g., [1ā€“17]). Such virus infection models can be very useful in the control of epidemic diseases and provide insights into the dynamics of viral load in vivo. Therefore, mathematical analysis of the virus infection models can play a significant role in the development of a better understanding of diseases and various drug therapy strategies. Most of the mathematical models of viral infection presented in the literature did not differentiate between the short-lived infected cells and chronically infected cells. The chronically infected cells produce much smaller amounts of virus than the short-lived infected cells and die at a much slower rate [18]. The virus dynamics model with chronically infected cells and under the effect of antiviral drug therapy was introduced in [18] as where , , , and are the concentration of the uninfected cells, short-lived infected cells, chronically infected cells, and free virus particles, respectively. The constant is the rate at which new uninfected cells are generated and is the natural death rate constant of uninfected cells. is the infection rate constant. The fractions and with are the probabilities that, upon infection, an uninfected cell will become either short-lived infected or chronically infected. and are the death rate constants of the short-lived infected cells and chronically infected cells, respectively. and are the average number of virions produced in the lifetime of the short-lived infected and chronically infected cells, respectively. The chronically infected cells produce much smaller amounts of virus than the short-lived infected cells and die at a much slower rate (i.e., and ). The free viruses are cleared with rate constant . The drug efficacy is denoted by and .

It is observed that the basic and global properties of model (1) are not studied in the literature. Moreover, model (1) did not take into consideration the immune response. During viral infections, the host immune system reacts with antigen-specific immune response. The immune system is described as having two ā€œarmsā€: the cellular arm, which depends on T cells to mediate attacks on virally infected or cancerous cells, and the humoral arm, which depends on B cells. The B cell is a type of blood cell which belongs to a group of white blood cells (WBCs) called lymphocytes. WBCs protect the body from infection. The main job of B cells is to fight infection. B cells get activated when an infection occurs and they produce molecules called antibodies that attach to the surface of the infectious agent. These antibodies either kill the infection causing organism or make it prone to attack by other WBCs. They play a major role in the immune system, which guards the body against infection. Virus infection models with antibody immune response have been analyzed by many researchers (see [19ā€“28]). However, in all of these works, the chronically infected cells have been neglected.

In this paper, we propose two virus infection models with antibody immune response and chronically infected cells. In the first model, bilinear incidence rate which is based on the law of mass-action is considered. The second model generalizes the first one where the incidence rate is given by a saturation functional response. The global stability of all equilibria of the models is established using the method of Lyapunov function. We prove that the global dynamics of the models are determined by two threshold parameters and . If , then the infection-free equilibrium is globally asymptotically stable (GAS), if , then the infected equilibrium without antibody immune response exists and it is GAS, and if then the infected equilibrium with antibody immune response exists and it is GAS.

2. Model with Bilinear Incidence Rate

In this section we propose a viral dynamics model with antibody immune response, taking into consideration the chronically infected cells. Based on the mass-action principle, we assume that the incidence rate of infection is bilinear; that is, the infection rate per virus and per uninfected cell is constant: where is the concentration of antibody immune cells. The viruses are attacked by the antibodies with rate . The antibody immune cells are proliferated at rate and die at rate . All the other variables and parameters of the model have the same meanings as given in (1).

2.1. Positive Invariance

We note that model (2)ā€“(6) are biologically acceptable in the sense that no population goes negative. It is straightforward to check the positive invariance of the nonnegative orthant by model (2)ā€“(6) (see, e.g., [6]). In the following, we show the boundedness of the solution of model (2)ā€“(6).

Proposition 1. There exist positive numbers , , such that the compact set is positively invariant.

Proof. To show the boundedness of the solutions we let ; then where . Hence , if where . Since , , and , then , , if . Let ; then where . Hence , if , where . Since and then and if , where .

2.2. Equilibria

System (2)ā€“(6) always admits an infection-free equilibrium , where . In addition to , the system can have an infected equilibrium without antibody immune response and an infected equilibrium with antibody immune response where We discuss the local stability of the infection-free equilibrium . At the infection-free equilibrium , the system has the Jacobian matrix given by The characteristic equation of the Jacobian matrix evaluated at is where We observe that (12) has two negative eigenvalues and . By the Routh-Hurwitz criterion, the remaining three eigenvalues of (12) have negative real parts if , , and . We have and if , then and Now we define the basic reproduction number for system (2)ā€“(6) as It follows that the equilibria and can be written as We note that , , , and are positive when and that when . Now we define another threshold parameter as Clearly .

From (16) we have the following statements:(i)if , then there exists only positive equilibrium ;(ii)if , then there exist two positive equilibria and ;(iii)if , then there exist three positive equilibria , , and .

2.3. Global Stability Analysis

In this section, we study the global stability of all the equilibria of system (2)ā€“(6) employing the method of Lyapunov function.

Theorem 2. For system (2)ā€“(6), if , then is GAS.

Proof. Define a Lyapunov function as follows: where ,ā€‰ā€‰, are positive constants to be determined below. Calculating the derivative of along the solutions of the system (2)ā€“(6) and applying , we obtain Let ,ā€‰ā€‰, be chosen such as The solution of (20) is given by The values of ,ā€‰ā€‰, given by (21) will be used throughout the paper. Then If then for all . Thus the solutions of system (2)ā€“(6) limit to , the largest invariant subset of . Clearly, it follows from (22) that if and only if , , and . Noting that is invariant, for each element of we have and , and then . From (5) we derive that Since , , then . Hence if and only if , , ,ā€‰ā€‰, and . It follows from LaSalle's invariance principle that the infection-free equilibrium is GAS when .

Theorem 3. For system (2)ā€“(6), if , then is GAS.

Proof. Define the following Lyapunov function: The time derivative of along the trajectories of (2)ā€“(6) is given by Applying we get Using the following equilibrium conditions for , then we have and We have that if , then . Since the arithmetical mean is greater than or equal to the geometrical mean, then if then for all . It can be seen that if and only if , , , , and . LaSalle's invariance principle implies global stability of .

Theorem 4. For system (2)ā€“(6), if , then is GAS.

Proof. We consider a Lyapunov function Further, function along the trajectories of system (2)ā€“(6) satisfies Using the following equilibrium conditions for , we get Thus, if , then and . Since the arithmetical mean is greater than or equal to the geometrical mean, then . It can be seen that if and only if , , , and . From (5), if , then and , so and hence is equal to zero at . So, the global stability of the equilibrium follows from LaSalle's invariance principle.

3. Model with Saturation Incidence Rate

In model (2)ā€“(6), the infection process is characterized by bilinear incidence rate . However, there are a number of reasons why this bilinear incidence can be insufficient to describe infection process in detail (see, e.g., [29ā€“31]). For example, a less than linear response in could occur when the concentration of viruses becomes higher, where the infectious fraction is high so that exposure is very likely [29]. Experiments reported in [32] strongly suggested that the infection rate of microparasitic infections is an increasing function of the parasite dose and is usually sigmoidal in shape (see, e.g., [33]). In [33], to place the model on more sound biological grounds, Regoes et al. replaced the mass-action infection rate with a dose-dependent infection rates. In this section, the incidence rate is given by a saturation functional response: where is a constant, which represents the saturation infection rate constant.

All the variables and parameters have the same meanings as given in model (2)ā€“(6).

3.1. Equilibria

Similar to the previous section, we can define two threshold parameters and for system (33)ā€“(37) as Clearly . It is clear that system (33)ā€“(37) has an infection-free equilibrium , where . In addition to , the system can have an infected equilibrium without antibody immune response , where and infected equilibrium with antibody immune response , where It is clear from (39) and (40) that(i)if , then there exists only positive equilibrium ;(ii)if , then there exist two positive equilibria and ;(iii)if , then there exist three positive equilibria , , and .

3.2. Global Stability Analysis

In this section, we study the global stability of all the equilibria of system (33)ā€“(37) employing the method of Lyapunov function and LaSalle's invariance principle.

Theorem 5. For system (33)ā€“(37), if , then is GAS.

Proof. Define a Lyapunov function as follows: Calculating the derivative of along the solutions of system (33)ā€“(37) and applying , we obtain Similar to the proof of Theorem 2, one can easily show that is GAS when .

Theorem 6. For system (33)ā€“(37), if , then is GAS.

Proof. Construct a Lyapunov function as follows: The derivative of along the trajectories of system (33)ā€“(37) is given by Applying we get Using the following equilibrium conditions for , we get We have that if , then where equality occurs at . LaSalle's invariance principle implies global stability of .

Theorem 7. For system (33)ā€“(37), if , then is GAS.

Proof. We consider a Lyapunov function as follows: Further, function along the trajectories of system (33)ā€“(37) satisfies Using the following equilibrium conditions for , we get Similar to the proof of Theorem 4, one can show that is GAS.

4. Numerical Simulations

We now use simple numerical simulations to illustrate our theoretical results for the two models. In both models we will fix the following data: ā€‰mmāˆ’3ā€‰dayāˆ’1, dayāˆ’1, ā€‰mm3ā€‰dayāˆ’1,ā€‰ā€‰ dayāˆ’1,ā€‰ā€‰,ā€‰ā€‰ day, dayāˆ’1, , , ā€‰mm3ā€‰dayāˆ’1, and ā€‰dayāˆ’1. The other parameters will be chosen below. All computations were carried out by MATLAB.

4.1. Model with Bilinear Incidence Rate

In this section, we perform simulation results for model (2)ā€“(6) to check our theoretical results given in Theorems 2ā€“4. We have the following cases.(i). We choose and ā€‰mm3ā€‰dayāˆ’1. Using these data we compute and . Figures 1, 2, 3, 4, and 5 show that the numerical results are consistent with Theorem 2. We can see that, the concentration of uninfected cells is increased and converges to its normal value ā€‰mmāˆ’3, while the concentrations of short-lived infected cells, chronically infected cells, free viruses, and antibody immune cells are decaying and tend to zero.(ii). We take and ā€‰mm3ā€‰dayāˆ’1. In this case, and . Figures 1ā€“5 show that the numerical results are consistent with Theorem 3. We can see that the trajectory of the system will tend to the infected equilibrium without antibody immune response . In this case, the infection becomes chronic but with no persistent antibody immune response.(iii). We choose and ā€‰mm3ā€‰dayāˆ’1. Then we compute and . From Figures 1ā€“5 we can see that our simulation results are consistent with the theoretical results of Theorem 4. We observe that the trajectory of the system will tend to the infected equilibrium with antibody immune response . In this case, the infection becomes chronic but with persistent antibody immune response.

We note that the values of the parameters , , and have no impact on the value of , since is independent of those parameters. This fact seems to suggest that antibodies do not play a role in eliminating the viruses. From the definition of , we can see that can be increased by increasing or decreasing .

Figures 1 and 4 show that the presence of antibody immune response (i.e., ) reduces the concentration of free viruses and increases the concentration of uninfected cells. This can be seen by comparing the virus and uninfected cell components in the equilibria and under the condition . For model (2)ā€“(6), simple calculation shows that It follows that if , then . From (2) and at any equilibrium point we have Clearly, is a decreasing function of . This yields that if , then and .

4.2. Model with Saturation Functional Response

In this section, we perform simulation results to check Theorems 5ā€“7. The parameter is chosen as ā€‰mm3. We have the following cases.(i). We take and ā€‰mm3ā€‰dayāˆ’1. Using these data, we compute and . The simulation results of this case are shown in Figures 6, 7, 8, 9, and 10. We can see that the numerical results are consistent with Theorem 5. It is observed that the viruses will be cleared and the uninfected cells will return to their normal value.(ii). To satisfy this condition, we take and ā€‰mm3ā€‰dayāˆ’1. This will give and . Figures 6ā€“10 show that the numerical results are consistent with Theorem 6. We see that the infected equilibrium is GAS, and the infection becomes chronic but with no persistent antibody immune response.(iii). This condition is satisfied by choosing and ā€‰mm3ā€‰dayāˆ’1. This yields and . Figures 6ā€“10 demonstrate the global stability of . Then, the infection becomes chronic but with persistent antibody immune response.

From the definition of the parameter , we can see that the value of the saturation infection rate constant has no impact on the value of . This means that saturation does not play a role in eliminating the virus. From the definition of , we can see that can be increased by increasing or decreasing and .

Figures 6 and 9 show that if the antibody immune response reduces the concentration of free viruses and increases the concentration of uninfected cells. For model (33)ā€“(37), simple calculation shows that As a result, if , then . From (33) and at any equilibrium point we have Then, is a decreasing function of . It follows that if then and .

5. Conclusions

In this paper, we have proposed two virus infection models with antibody immune response taking into account the chronically infected cells. In the first model we have assumed that the incidence rate of infection is bilinear while in the second model the incidence rate is given by saturation functional response. We have shown that the dynamics of the models are fully determined by two threshold parameters and . The parameter determines whether a chronic infection can be established while determines whether a persistent antibody response can be established. By constructing Lyapunov function and using LaSalle's invariance principle, we have investigated the global stability of all equilibria of the two models. We have proven that if then the infection-free equilibrium is GAS, and the viruses are cleared. If , then the infected equilibrium without antibody immune response exists and it is GAS, and the infection becomes chronic but with no persistent antibody immune response. If , then the infected equilibrium with antibody immune response exists and it is GAS, and the infection is chronic with persistent antibody immune response. Numerical simulations have been performed for the two models. Our simulation results confirm the analytic results given in Theorems 2ā€“7.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. 130-078-D1434. The authors, therefore, acknowledge with thanks DSR technical and financial support. The authors are also grateful to Professor Malay Banerjee and to the anonymous reviewers for constructive suggestions and valuable comments, which improve the quality of the paper.