Journal of Applied Mathematics

Volume 2013, Article ID 381286, 8 pages

http://dx.doi.org/10.1155/2013/381286

## Extinction of Disease Pathogenesis in Infected Population and Its Subsequent Recovery: A Stochastic Approach

^{1}Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata, West Bengal 700032, India^{2}Department of Mathematics, Barasat College, Kolkata 700126, India^{3}Department of Chemistry, Narula Institute of Technology, Kolkata 700109, India^{4}Agricultural and Ecological Research Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India^{5}Department of Stochastics, Institute of Mathematics, Budapest University of Technology and Economics, Budapest 1521, Hungary

Received 5 February 2013; Accepted 21 May 2013

Academic Editor: Xinyu Song

Copyright © 2013 Priti Kumar Roy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A stochastic mathematical model of host-pathogen interaction has been developed to estimate the time to extinction of infected population. It has been assumed in the model that the infected host does not grow or reproduce but can recover from pathogenic infection and move to add to the susceptible host population using various drugs or vaccination. Extinction of infected population in host-pathogen interaction depends significantly upon the total population. Here, we consider an extension of our previous work with the stochastic approach to predict the time to extinction of disease pathogenesis. The optimal control approach helped in designing an innovative, safe therapeutic regimen where the susceptible host population enhanced with simultaneous decrease in the infected population. By means of an optimal control theory paradigm, it has also been shown in our preceding research paper that the cost-effective combination of treatment may depend on the population size. In this research paper, we have studied an approximation which is derived in favor of quasi-stationary distribution along with the expected time to extinction for the model of host-pathogen interactions. The complete study has been roofed through the stochastic approach in context that disease pathogenesis is to be extinct and infected population are going to be recovered. Numerical simulation is also done to confirm the analysis.

#### 1. Introduction

The modeling of epidemic diseases is an age-old problem [1] and makes a sincere effort in understanding the development of mathematical models for epidemics from the 18th century to the present day. These models are shown to be of used in predicting and controlling the spread of infection. But in the recent context of epidemiological research, microbial pathogenesis reflects the interaction between two entities, host and pathogen, which is somewhat related to the predator-prey model [2, 3]. These classes of models are also relevant for host-parasite type of models. Host-pathogen models are mathematical prototypes pertaining to epidemiology, and are of immense importance in view of the emergence and reemergence of epidemiological diseases in the present day global scenario [4, 5].

Application of mathematical concepts and techniques to analyze host-pathogen interactions was done by various researchers [6–8]. A very recent study throws light on the antiviral drug treatment along with immune activator IL-2 and optimal control in disease pathogenesis by deterministic model formulation [9]. In this research paper, a conventional host-pathogen model has been considered including the recovery of the infected individuals to the healthy susceptible organisms. In this case, the host population is divided into two classes: susceptible , that is, healthy organisms, and infected individuals . Pathogens cause infection to host population by transforming to . Over the last several years, many researchers focused their attention on the mathematical and biological aspects of host-pathogen interactions. Beltrami and Carroll [10] as well as Venturino, worked on the role of viral disease considering a three-species model of susceptible and infected phytoplankton as well as their predator [4]. Hethcote and Driessche [11] formulated SIS epidemiological models where delay has been incorporated corresponding to the infectious period and disease-related deaths. Another pioneering work was reported by [12] where the transport of coevolved host-pathogen systems into new environment leads to the evolution of altered levels of pathogen aggressiveness, if transmission rates are different in the new environment. With the growing research in the prey-predator and other prototypical systems, it is apparent that the pathogen or viral growth through replication influences the model dynamics. This has been emphasized by Bairagi et al. in a subsequent communication [13]. In another pioneering communication in recent years, Bairagi et al. carried a comparative study of the prey-predator model with several response functions [14].

Epidemiological modeling uses both deterministic and stochastic approaches in host-pathogen interactions [15]. Both model types have their respective advantages and weaknesses. Deterministic models may be considered as an approximation of a corresponding stochastic model. An important difference between the two is that the stochastic model deals with a finite population size, while the deterministic model deals only with proportion [16]. The deterministic version of the model can be derived as an approximation of the stochastic version using the state variables which are susceptible host , infected host , and pathogen where the total population comprising the susceptible and infected host is . Recurrence can be explained by the combined influence of epidemic and demographic forces. In stochastic models, infection will eventually become extinct, and time to extinction is an important measure of the persistence of the infection [17].

This paper specifies a stochastic approach of our earlier work [18] designed for some widespread infection in closed population. It is a well-known fact that beyond a threshold value, that is, basic reproduction ratio, which is determined by the parameters of the model, the deterministic model predicts that the proportion of infected individuals will approach a positive endemic level as time approaches infinity. However, the stochastic model predicts that the infection will become extinct. Generally, stochastic version of the model is compared with prey-predator interaction or more precisely the host-pathogen interaction, and thus the concept of the quasi-stationary state is an important aspect of this model [19]. It should be mentioned here that the stochastic logistic model can be interpreted as SIS model and that SIS model is used for infection that gives no immunity. In this case, it is assumed that an individual who recovers from an infection remains susceptible to further infection. The time to extinction from quasi-stationary distribution is a measure of the persistence of the infection. The stationary distribution is to be an additional importance since the expected time to extinction from quasi-stationary can be expressed in terms of this distribution. A major goal of the analysis of the stochastic model is therefore to derive an approximation of the quasi-stationary distribution. This derivation is based on a diffusion approximation of the stochastic discrete state model.

The model is analyzed in two different avenues, analytical and numerical. In our research paper, we have found an approximation for the marginal distribution of the infected population in quasi-stationary condition and its time to extinction. The time to extinction of the infected population depends on the total population, and the study has been carried out through stochastic approach so that infected population is to be extinct and increase the susceptible population. However, being an extension of our previous work [18], where deterministic model was used and optimal control theory applied to predict the decrease of infected population along with cost-effective combination of treatment, our present work emphasizes on the time to extinction of the infected host population or the disease by stochastic model formulation and its subsequent analysis and evaluation.

#### 2. The Deterministic Model

We consider the three components of the basic three-dimensional deterministic host pathogen model [18] consisting of a host population, whose concentration is denoted by ( = number of host per designated area) and a pathogen population inflicting infection in the host population whose concentration is denoted by ( = number of pathogens per designated area). The following differential equations are formed at initial dynamics on the change of host pathogen interaction with time :

Here in the presence of pathogenic infection, the host population is divided into two disjoint classes, susceptible host and infected host . In the ideal case of no pathogen, the growth of susceptible host population follows the logistic law [19] implying that this growth is entirely controlled by an intrinsic birth rate constant with a carrying capacity . is the force of infection through contact with pathogens, and pathogens maximally infect susceptible hosts per day. This infection rate is half maximal at susceptible host population density of host. is the intensity of infection by infected host, and is the death rate constant of . Rate of cell lysis (replication of pathogens) is , and the natural death rate of pathogens is denoted as . We assume that the infected hosts do not grow or reproduce, but they can recover from pathogenic infection and move to the susceptible host population. Such recovery would stem out from immunization or vaccination. We consider a recovery rate of infected host to be denoted by .

In spirit of Bonhoefer et al., we employed a simplified system with two components, the susceptible and the infected host. It is assumed that at equilibrium point , so can be eliminated by putting . With this choice of model, (1) reduces to

where and .

An alternative deterministic formulation of the reduced model (2), one assumes birth and death rate functions and , respectively, of

where and . Here, and are intrinsic rate, and are crowding coefficient [20], and thus intrinsic growth rate and carrying capacity are defined as follows: and [20],

#### 3. The Stochastic Model Formulation

There are two state variables, namely, the number of susceptible hosts and the number of infected hosts at time . They jointly take values in the state space . The joint distribution of and at time is denoted by

We use this notation even when and/or are negative, with the convention that is equal to zero. The model is based on the following four basic events, that is, birth of a susceptible host, death of a susceptible host, infection of an uninfected host, and death or recovery of an infected host. The transition rates of the model are shown in Table 1.

#### 4. Description of the Transition States

The total number of population is increased by unity, if there is a birth of susceptible host for a small time interval . But to make the population to be the same, we should assume that there must be a natural death of susceptible host . These phenomena are captured through the first two rows of the transition matrix. On the other hand, if there is an infection in the susceptible host it can be balanced by an increase of an infected host, . The susceptible class is infected either by direct infection with the pathogen or by replication of viral generated within the infected population. This occurs one at a time, and so the increases of the infected class are reflected by the rise of unity in the transition state. If there is an infected host, the natural death should be reflected through a natural birth of susceptible host. At the end, naturally, the recovery of infected host must be reinstated to the susceptible host.

#### 5. Formulation of Kolmogorov’s Forward Equation

We are supposing that in a time interval of infinitesimally little length , the probability of precisely one birth (or one death) is birth rate (or death rate) × + intuitively more than one occasion (birth and/or death) in . We as well believe the prospect , where . The Kolmogorov’s forward equations for the representation can be written as

where

Note that all proceedings consisting of more than one birth or more than one death are incorporated in the expression

#### 6. The Quasi-Stationary Distribution

We obtain primarily a deferential equation that will be used afterward in this subsection and in the discussion of time to extinction in the subsequent subsection. Place in Kolmogorov’s Forward equations (8), and put in the explicit expressions for the evolution rates given above. In addition, commence

to indicate the subsidiary allocation of the number of infected individuals at time . By summing the forward equations intended for over all -values, we acquire

After that, we develop an organization of equations for the quasi-stationary distribution. The state possibility conditioned on not being engrossed is signified and specified by

We initiate to stand for the marginal distribution for the number of infected individuals at time , habituated on not having attained any condition in the engrossing position. By differentiating the expression for in (11) and relating (10), we get hold of

By pertaining the forward equations for in (8), we achieve the following scheme of differential equations for the conditional state probabilities :

The quasi-stationary distribution is the stationary solution of this system of equations.

#### 7. The Distribution of the Time to Extinction

Two initial distributions are principally fascinating. One is the quasi-stationary distribution, and another communicates to one infected individual. The distribution of the time to extinction can be unwavering if we can resolve the Kolmogorov forward equations (8) for . The reason is that the happening that surpasses is equal to the event that the number of infected individuals at time is positive. By allowing for the balancing events we attain

Therefore, the cumulative distribution function of the time to extinction at time equals the marginal possibility in which the number of infected individuals at time equals 0. The distribution of the time to extinction is particularly straightforward when the opening distribution is equal to the quasi-stationary distribution. Let us indicate the time to extinction from quasi-stationarity by explicitly. We demonstrate that has an exponential distribution and that its predictable value is equal to

To gain this consequence, we put in (12). Thus, we are guided to the initial value problems

with solutions

By adding these expressions of over all , we obtain

The differential equation for in (10) can now be answered since the right-hand side of this equation is known from above. Remembering that we have the initial value , we acquire . This institutes the claim that has an exponential distribution with expected value agreed by (15).

#### 8. Diffusion Approximation and the Approximation of Quasi-Stationary Distribution

In this section, we have derived the diffusion approximation for the process formulated in Section 2. In order to approximate the quasi-stationary distribution, we consider the two-dimensional process which is represented by the set of differential equation (4)

The main result is the quasi-stationary distribution is approximated by a bivariate normal distribution, if is sufficiently large.

The critical point of the rescaled deterministic model that corresponds to a pathogenic infection is denoted by , where

The changes in the scaled state variables and during the time interval from to are denoted by and using Table 2: and . From the hypotheses of the original process, we determine the mean and the covariance of the vector with components and . We begin with the mean

The Jacobian matrix of the vector with respect to is denoted by

where , and .

Next, we determine the covariance matrix of the vector of changes in the state variables during the time interval

The matrix is approximated by evaluating it at the critical point corresponding to the pathogenic infection level

where and .

For large , the process is approximated by a stable bivariate Ornstein-Uhlenbeck process, with local drift matrix and local covariance matrix . Then, the stationary distribution of the Ornstein-Uhlenbeck process is bivariate normal with mean and variance through the relationship

where the superscript is used to denote the transpose. After solving the above equation, we get

Note that the parameter should also satisfy , . Thus, the diffusion approximation led to the conclusion that the marginal distribution of the infected population size in quasi-stationarity is approximately . To achieve consistency with the fact , the approximating normal distribution is modified by truncation at . Hence, we have the following approximation

where and are, respectively, the standard normal c.d.f. and the standard normal p.d.f.

##### 8.1. The Expected Time to Extinction

We find the expected time to extinction () from quasi-stationary distribution. The expected time to extinction is given by

We see that the expected time to extinction is a function of population size , which is a function when the population size is increasing (Figure 2).

#### 9. Model Modification under Immune Host

In model (4), we consider the recovery class to be joined with the susceptible host, and it may increase the population size of uninfected host. It is to be noted that sometimes for a specific disease, the recovery from infected host may not come back to the susceptible class. For such a case, after the recovery from infected host, this specific recovery class becomes immunized and can be classified into a new category. Under this assumption the modified model can be written from model (4) by just ignoring the term , in the growth equation of susceptible host. The algebraic manipulations for obtaining the quasi-stationary distribution and time to extinction for revised model are pretty similar with that of model (4).

#### 10. Numerical Illustration

In the present study, a stochastic approach has been adopted to eradicate the infected population from a system and add on to the susceptible population. This is achieved by immunization or vaccination, and infected population recovers to susceptible population using the recovery rate . The stochastic model in this case predicts the extinction of infection, and the time to extinction is the important measure of the persistence of infection.

Figure 1 represents the normalized simulated marginal distribution profile of the infected host in quasi-stationary state for , , and . It is observed that when the population size is , the quasi-stationary distribution is truncated and skewed, while when is increased to and subsequently to , the distribution displays higher kurtosis and enhanced symmetry. The population distribution of the infected class thus becomes much narrower for higher total strength of the susceptible and infected individuals. From Figure 2, we get the expected time to extinction of the infected host from quasi-stationary state when the total population varies from to . In this case, it is observed that the expected time to extinction increases with the increase in the total population . When is small, the number of infected population is also expected to be small, and at that time, the natural immune system can decelerate the rapid growth of the infected cells and add it to the susceptible class. This results in the lesser time to extinction at the initial stages of . But when the population size is high, the number of infected class is also expected to rise, and the immune system cannot control the rapid growth of the infected host population, and this leads to the high extinction time as observed. Thus, it is evident from the numerical analysis that at the onset of the time to extinction, the infected individuals get transformed and move to the susceptible population.

#### 11. Discussion and Conclusion

In this paper, we have presented a basic mathematical model of host-pathogen interaction using the stochastic approach based on the concept of quasi-stationary and the diffusion approximation result. The time to extinction has been predicted as a function of the total population size . The numerical simulation reveals that as the total population increases, the quasi-stationary distribution inclines to reduced skewness and narrower distribution of the infected class under consideration. The time to extinction of the infected class and its transformation to the susceptible class also vary with the total population size and are found to exhibit a gradual rise with increasing value of the total population. Since the deterministic version is an approximation of the stochastic model, the estimation of the time to extinction is unlikely to be feasible with the deterministic version. So, we are inclined to conclude that the deterministic version of the model gives an approximation to the stochastic model.

The diffusion approximation led to the conclusion that the marginal distribution of the infected host population size in quasi-stationarity is approximately . To achieve consistency with the fact that , the approximating normal distribution is modified by truncating at . We conclude from the above result that as the total population size of the system becomes large , the marginal distribution of is approximately , and the marginal distribution of is approximately . Further, the covariance is approximated by . The numerical simulation confirms the analysis which is relevant in Figures 1 and 2. Further, the numerical simulation reveals that as the total population increases, the quasi-stationary distribution changes from positively skewed to symmetrical nature along with increases in the time to extinction.

The entire study has been carried out in a different angle to focus on the time to extinction of disease pathogenesis, and this reflects an advancement of our earlier reported work based on the decreasing infected population by deterministic modeling. Thus, we can conclude that stochastic modeling provides a more accurate prediction in finding out the expected time to extinction of infected population and hence disease pathogenesis.

#### References

- J. Gani, “Modelling epidemic diseases,”
*Australian and New Zealand Journal of Statistics*, vol. 52, no. 3, pp. 321–329, 2010. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - H. I. Freedman, “A model of predator-prey dynamics as modified by the action of a parasite,”
*Mathematical Biosciences*, vol. 99, no. 2, pp. 143–155, 1990. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. Chattopadhyay and N. Bairagi, “Pelicans at risk in Salton sea—an eco-epidemiological model,”
*Ecological Modelling*, vol. 136, no. 2-3, pp. 103–112, 2001. View at Publisher · View at Google Scholar · View at Scopus - E. Venturino, “Epidemics in predator-prey models: disease in the predators,”
*IMA Journal of Mathematics Applied in Medicine and Biology*, vol. 19, no. 3, pp. 185–205, 2002. View at Google Scholar · View at Scopus - D. Mukherjee, “Persistence in a prey-predator system with disease in the prey,”
*Journal of Biological Systems*, vol. 11, no. 1, pp. 101–112, 2003. View at Publisher · View at Google Scholar · View at Scopus - J. Chattopadhyay, A. Mukhopadhyay, and P. K. Roy, “Effect of viral infection on the generalized Gause model of predator-prey system,”
*Journal of Biological Systems*, vol. 11, no. 1, pp. 19–26, 2003. View at Publisher · View at Google Scholar · View at Scopus - I. Nåsell, “Stochastic models of some endemic infections,”
*Mathematical Biosciences*, vol. 179, no. 1, pp. 1–19, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - P. K. Roy and B. Chattopadhyay, “Host pathogen interactions with recovery rate: a mathematical study,” in
*Proceedings of the World Congress on Engineering (WCE '10)*, Lecture Notes in Engineering and Computer Science, pp. 521–526, London, UK, June-July 2010. View at Scopus - A. N. Chatterjee and P. K. Roy, “Anti-viral drug treatment along with immune activator IL-2: a control-based mathematical approach for HIV infection,”
*International Journal of Control*, vol. 85, no. 2, pp. 220–237, 2012. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - E. Beltrami and T. O. Carroll, “Modelling the role of viral disease in recurrent phytoplankton blooms,”
*Journal of Mathematical Biology*, vol. 32, no. 8, pp. 857–863, 1994. View at Publisher · View at Google Scholar - H. W. Hethcote and P. van den Driessche, “Two SIS epidemiologic models with delays,”
*Journal of Mathematical Biology*, vol. 40, no. 1, pp. 3–26, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - R. A. Ennos, “The introduction of lodgepole pine as a major forest crop in Sweden: implications for host-pathogen evolution,”
*Forest Ecology and Management*, vol. 141, no. 1-2, pp. 85–96, 2001. View at Publisher · View at Google Scholar · View at Scopus - N. Bairagi, P. K. Roy, R. R. Sarkar, and J. Chattopadhyay, “Virus replication factor may be a controlling agent for obtaining disease-free system in a multi-species eco-epidemiological system,”
*Journal of Biological Systems*, vol. 13, no. 3, pp. 245–259, 2005. View at Publisher · View at Google Scholar · View at Scopus - N. Bairagi, P. K. Roy, and J. Chattopadhyay, “Role of infection on the stability of a predator-prey system with several response functions—a comparative study,”
*Journal of Theoretical Biology*, vol. 248, no. 1, pp. 10–25, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - T. Britton and P. Neal, “The time to extinction for a stochastic SIS-household-epidemic model,”
*Journal of Mathematical Biology*, vol. 61, no. 6, pp. 763–779, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - I. Nåsell, “A new look at the critical community size for childhood infections,”
*Theoretical Population Biology*, vol. 67, no. 3, pp. 203–216, 2005. View at Publisher · View at Google Scholar · View at Scopus - I. Nåsell, “On the time to extinction in recurrent epidemics,”
*Journal of the Royal Statistical Society B*, vol. 61, part 2, pp. 309–330, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - P. K. Roy and J. Mondal, “Host pathogen interactions: insight of delay response recovery and optimal control in disease pathogenesis,”
*Engineering Letters*, vol. 18, no. 4, 2010. View at Google Scholar · View at Scopus - M. Carletti, “On the stability properties of a stochastic model for phage-bacteria interaction in open marine environment,”
*Mathematical Biosciences*, vol. 175, no. 2, pp. 117–131, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J. H. Matis, T. R. Kiffe, E. Renshaw, and J. Hassan, “A simple saddlepoint approximation for the equilibrium distribution of the stochastic logistic model of population growth,”
*Ecological Modelling*, vol. 161, no. 3, pp. 239–248, 2003. View at Google Scholar · View at Scopus - J. H. Matis, T. R. Kiffe, and P. R. Parthasarathy, “On the cumulants of population size for the stochastic power law logistic model,”
*Theoretical Population Biology*, vol. 53, no. 1, pp. 16–29, 1998. View at Publisher · View at Google Scholar · View at Scopus