Abstract

This paper is devoted to develop a nonlocal and time-delayed reaction-diffusion model for HIV infection within host cell-to-cell viral transmissions. In a bounded spatial domain, we study threshold dynamics in terms of basic reproduction number for the heterogeneous model. Our results show that if , the infection-free steady state is globally attractive, implying infection becomes extinct, while if , virus will persist in the host environment.

1. Introduction and Model Derivation

Infectious diseases have been threatening human health; infectious diseases such as smallpox, cholera, and acquired immune deficiency syndrome (AIDS) have brought great disaster to the national economy of a country and people’s livelihood. In order to control the spread of infectious diseases, researchers have proposed a great deal of mathematical models to study the evolution behavior of infectious diseases from different views [17]. Mathematical models have been confirmed to be an effective and valuable approach to understand the dynamical behavior of infectious diseases [811]. Many diseases are caused by a virus, for example, human immunodeficiency virus (HIV) can destroy the immune cells, reduce human immunity, and ultimately lead to AIDS. A lot of mathematical models have been built to explain the virus infection process from different views [1215]. In recent years, the virus infection dynamical models incorporating spatial dispersion [16, 17] have received widely attentions. For example, under spatial domain equipped with suitable boundary conditions, K. Wang and W. Wang [18] established the existence of travelling wave solutions via the geometric singular perturbation method and the spatial domain is assumed to be one dimensional, that is, . And then, Brauner et al. [19] extended the works [18] to a two-dimensional square domain with periodic boundary conditions by introducing a new parameter which is the largest eigenvalue of some Sturm-Liouville problem and taking the heterogeneous reproductive ratio into account. Recently, Wang et al. [20] noticed that a realistic domain should be bounded but is typically not a square, under no-flux boundary condition (homogenous Neumann boundary condition) in a general bounded domain with smooth boundary , and threshold dynamics of the model were investigated by appealing to the theory of uniform persistence and the comparison theorem. Some further developments have been performed on the virus infection dynamical models with diffusion term (see, for example, Hattaf and Yousfi [21], Wang and Ma [2224], Lai and Zou [25], Wang and Xu [26], and Wang et al. [27]).

It is widely known that virus-to-cell transmission and direct cell-to-cell transmission are two predominant infection modes in within-host environment (see, for example, [28, 29]). The virus-to-cell transmission has iterative process: the binding of virus to a receptor on the surface of CD4 T cells, the fusion of virus with host cells and then the release of genetic materials, the transcription of genetic materials in an infected cell, the assembling of virus inside the infected cells, and budding through the membrane of the infected cells [30]. The cell-to-cell transmission process is due to that a large number of viruses can be transferred from infected cells to uninfected cells through the formation of virally induced structures termed virological synapses [31]. In order to examine the effects of both diffusion and spatial heterogeneity, Wang et al. [32] proposed the following virus infection dynamical model incorporating cell-to-cell transmission under the homogeneous Neumann boundary condition:

In model (1), is the production rate of newly produced uninfected cells. , , and represent the death rates of uninfected cells, infected cells, and free viruses, respectively. denotes the transmission coefficient for virus-to-cell infection. represents the transmission coefficient for cell-to-cell infection. denotes the production rate of virus due to the lysis of infected cells. is the diffusion coefficient and is Laplacian. For model (1), the virus dynamics are fully determined by an important parameter which is called basic reproduction number , in the following sense: if , the infection-free steady state is globally asymptotically stable, while if , the model is uniformly persistent and the infection steady state is globally asymptotically stable.

The mobility of cells in the latent period will result in a delay term with spatial averaging on the spatial habitat [33]. To formulate this process with the latency properly, we introduce an infection age variable . The infected cells can be divided into two epidemiological categories: latently infected cells () and actively infected cells (). Let represent the density of cells with infection age at time and habitat . From the biological perspective, we assume that all populations remain confined to the region for all time and subject to no-flux boundary condition for : where is the outward normal to the smooth boundary . We adopt the standard model on describing age-structured population with spatial diffusion [34]. Then

In (3), is the gradient of with respect to the spatial variable and represents the divergence of . represents the diffusion rate of cells at age and habitat . is the natural death rate which is independent of the infection age . Let be the average incubation period. Then

In order to make the model mathematically tractable yet without losing the main characteristics, we make some assumptions for function as follows:

Integrating both sides of (3) from to and from to , it then follows that and

For biological reasons, we assume that (see, for example, [35]), since the recruitment of newly infected cells is divided into two parts: the contact of uninfected cells and virus and the contact of uninfected cells and infected cells. Thus, adopting Beddington-DeAngelis functional response which was firstly proposed by [36, 37] leads to the following condition: where and are the measure of virus interference during infection and and determine how fast the infection rate approaches its saturation value (see, for example, [36, 37]).

In the following, we determine by the method of characteristics. Let with . For , it then follows that

Regarding as a parameter and solving the above equation, we obtain that where is the fundamental solution associated with the partial differential operator , which satisfies , . , and , (see, [38], Chapter 1).

Since , , we get the formula for for :

It then follows that and

With all these assumptions, the HIV infection dynamics modeling cell-to-cell with nonlocal infections can be described by the following nonlocal and time-delayed reaction-diffusion equations:

Since is decoupled from other equations for model (14), we consider the following dynamical model: for , . We consider a closed environment in the sense that the fluxes for each of these subpopulations are aero, and hence we propose no-flux condition on the boundary: and initial conditions

We now briefly outline the plan of this paper. In Section 2, we establish basic reproduction number . In Section 3, we study threshold dynamics in terms of basic reproduction number for the heterogeneous model. At last, we give discussions and biological implications in Section 4.

2. Basic Reproduction Number

Let be the Banach space with the supremum norm by , , where denotes the Euclidean norm in and denotes the transpose of the vector, and . We easily get that is a closed cone of .

Define with the norm and , then is an ordered Banach space. For any given continuous function for , define by , . Let and . Suppose that , , and are the strongly continuous semigroups associated with , , and subject to homogeneous Neumann boundary conditions, respectively, that is, for any , , and , . and are the Green functions associated with , , and , respectively, subject to homogenous Neumann boundary conditions. According to Section 7.1 and Corollary 7.2.3 in [39], it follows that is strongly positive and compact for each . We can also obtain that , , is a strongly continuous semigroup. Let be the generator of , respectively . Then, we get that is a semigroup generated by the operator defined on .

In view of Lemma 2.2 in [35], model (15) exists a unique infection-free steady state Linearizing model (15) at the infection-free steady state , we get the following linearized system: satisfying the following boundary conditions

We consider the following model:

Substituting and into and of model (21) which is called the infectious compartments results in the following eigenvalue problem:

By the similar argument to Theorem 7.6.1 in [39], it then follows that eigenvalue problem (22) has a principal eigenvalue with a positive eigenfunction.

We further consider the nonlocal eigenvalue problem associated with

From Theorem 2.2 in [40], we get some useful information eigenvalue problem (23).

Lemma 2.1. The eigenvalue problem (23) has a principal eigenvalue with a strictly positive eigenfunction, and for any , has the same sign as .

We adopt the same ideas as in [4144] to define basic reproduction number of model (15). For this purpose, we define the positive linear operator as follows: where

In order to compute basic reproduction number , we assume that both infected cells and virus are near the infection-free steady state and introduce initial infected cells and virus at time , where the distribution of initial infected cells and virus is described by . From model (15) and as time evolves, those distributions can reach at time . Thus, the total distribution of new infected cells is

Similarly, the total distribution of new virus can be described as

It then follows that where is the next infection operator (see, for example, [4144]). It maps the initial distribution of infected cells and virus to the total distribution of new infected cells and virus produced during the infection period.

By [4144], we define the spectral radius of as the basic reproduction number of model (15), namely,

In view of the general results developed in [42], we have the following observation.

Lemma 2.2. has the same sign as .

3. Threshold Dynamics in a Bounded Spatial Domain

In this section, we study threshold dynamics in terms of for model (15), assuming a bounded spatial domain . We assume that all the parameters , , , , , , , and are spatially dependent. Assume these functions are positive, continuous, and bounded in the domain . Let

For any and , we define by

Obviously, is locally Lipschitz. Then, model (15) can be rewritten as the following abstract functional differential equation: where .

Lemma 3.1. For any given , there exists a unique mild solution of model (15) defined on its maximal interval of existence with , where . Further, , , and are a classical solution of model (15), .

Proof 1. For any and for sufficiently small , we obtain that Thus, . This implies that for all From Corollary 4 in [45] (see, also [46], Corollary 8.1.3), we obtain the conclusion stated in Lemma 3.1. The proof is completed.

We are now in the position to study the well-posedness of model (15) in the sense of the following theorem.

Theorem 3.1. For any , model (15) has a unique solution on with . Moreover, the solution semiflow has a compact global attractor for .

Proof 2. It is easy to see that model (15) defines a semiflow by For any fixed , from the first equation of model (15), we obtain that From the comparison principle, we easily get that there exists such that for . From the second equation of model (15), we have that There exists such that for . From the third equation of model (15), we get There also exists such that for .
Consequently, the existence of solutions of model (15) claimed in Lemma 3.1 is indeed global (i.e., ). The solution semiflow is point dissipative. Moreover, by Theorem 2.2.6 in [46], we have that is compact for any . Thus, from Theorem 3.4.8 in [47], we know that has a compact global attractor in for . We complete the proof.

The following results will play an important role in establishing uniform persistence of model (15).

Lemma 3.2. Let be the solution of model (15) with , then we have the following:
(i) If there exists some such that and , then for all , .
(ii) It holds that for any , , and uniformly for .

Proof 3. From the second and third equations of model (15), it follows that with the following boundary conditions If and for some , it follows from the strong maximum principle (see, for example, [48] p. 172, Theorem 4) and the Hopf boundary lemma (see, for example, [48] p. 170, Theorem 3) that , , and for , , that is, the conclusion holds.
From the first equation of model (15), we get Let be the solution of It then follows that for all and . It follows from the standard parabolic comparison principle that uniformly for . The proof is completed.

The following theorem indicates that is a threshold quantity for virus extinction or persistence.

Theorem 3.2. Suppose is the solution of model (15) with . Then, the following statements hold: (i)If , then the infection-free steady state of model (15) is globally attractive.(ii)If , then model (15) has at least one coexistence steady state and there exists such that any nonnegative solution with and , we getuniformly for all .

Proof 4. From Lemma 2.2, we easily obtain that when . Since there exists sufficiently small such that . For fixed , there exists such that for all , . Therefore, for all , from the second and third equations of model (15), it follows that According to Lemma 2.1, we get , and there exists a positive eigenfunction corresponding to . It then obtains the following linear system: with the following boundary conditions admitting a solution . Then, for any given , there exists such that By the comparison principle, it follows that Thus, uniformly for . Hence, by [35], it then follows that uniformly for which ends the proof of Part (i).

For , we employ the persistence theory in [49] to study uniform persistence of model (15). Define

Obviously, we have that

In view of Lemma 3.2, it then follows that is positively invariant for the solution semiflow . Define

Let be the omega limit set of the orbit of through , and set .

Claim 1. .
For any given , it follows that for all . Since , we easily have that for each either or For the case where , , combined with the third equation of model (15), it can be obtained that uniformly for . Further, from the first equation of model (15) and the theory of asymptotically autonomous semiflows [50], it then follows that uniformly for . In the case where for some , by parabolic maximum principle, we have that for all and . Therefore, we get for all . From the third equation of model (15), we obtain that for all , which is a contradiction. Hence, we get . We complete the proof of Claim 1.

Claim 2. Similar to the proof of Lemma 3.3 in [17], we easily get that is a uniform weak repeller in the sense that We define a continuous function by We easily verify that . By Lemma 3.3, has the property that if either and or , then for all . Consequently, we have that is a generalized distance function for the semiflow (see [49]). From the above claims, it follows that any forward orbit of in converges to which is isolated in , and , where is the stable set of (see, for example, [49]). Further, we find that no subset of forms a cycle in from to . By Theorem 3 in [49], we prove that there exists such that for any In view of Theorem 1.3.6 in [51], we obtain that model (15) has at least one coexistence steady state if . The proof is completed.

4. Discussions and Biological Implications

This paper is devoted to develop a nonlocal and time-delayed reaction-diffusion model for HIV infection within host cell-to-cell viral transmissions. For model (15), the virus dynamics are fully determined by an important parameter which is called basic reproduction number , in the following sense: if , the infection-free steady state is globally attractive, implying infection becomes extinct; while if , the virus will persist in the host environment.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by Shandong Provincial Natural Science Foundation (no. ZR2015AQ001), the National Natural Science Foundation of China (no. 11371230), the SDUST Research Fund (2014TDJH102), Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources by Shandong Province, and SDUST.