Stability and bifurcations in an epidemic model with varying immunity period

An epidemic model with distributed time delay is derived to describe the dynamics of infectious diseases with varying immunity. It is shown that solutions are always positive, and the model has at most two steady states: disease-free and endemic. It is proved that the disease-free equilibrium is locally and globally asymptotically stable. When an endemic equilibrium exists, it is possible to analytically prove its local and global stability using Lyapunov functionals. Bifurcation analysis is performed using DDE-BIFTOOL and traceDDE to investigate different dynamical regimes in the model using numerical continuation for different values of system parameters and different integral kernels.Comment: 16 pages, 5 figure

Time-delayed model of autoimmune dynamics

Among various environmental factors associated with triggering or exacerbating autoimmune response, an important role is played by infections. A breakdown of immune tolerance as a byproduct of immune response against these infections is one of the major causes of autoimmune disease. In this paper we analyse the dynamics of immune response with particular emphasis on the role of time delays characterising the infection and the immune response, as well as on interactions between different types of T cells and cytokines that mediate their behaviour. Stability analysis of the model provides insights into how different model parameters affect the dynamics. Numerical stability analysis and simulations are performed to identify basins of attraction of different dynamical states, and to illustrate the behaviour of the model in different regime

Viral Dynamics of Delayed CTL-inclusive HIV-1 Infection Model With Both Virus-to-cell and Cell-to-cell Transmissions

We consider a mathematical model that describes a viral infection of HIV-1 with both virus-tocell and cell-to-cell transmission, CTL response immune and four distributed delays, describing intracellular delays and immune response delay. One of the main features of the model is that it includes a constant production rate of CTLs export from thymus, and an immune response delay. We derive the basic reproduction number and show that if the basic reproduction number is less than one, then the infection free equilibrium is globally asymptotically stable; whereas, if the basic reproduction number is greater than one, then there exist a chronic infection equilibrium, which is globally asymptotically stable in absence of immune response delay. Furthermore, for the special case with only immune response delay, we determine some conditions for stability switches of the chronic infection equilibrium. Numerical simulations indicate that the intracellular delays and immune response delay can stabilize and/or destabilize the chronic infection equilibrium

Study of Virus Dynamics by Mathematical Models

This thesis studies virus dynamics within host by mathematical models, and topics discussed include viral release strategies, viral spreading mechanism, and interaction of virus with the immune system. Firstly, we propose a delay differential equation model with distributed delay to investigate the evolutionary competition between budding and lytic viral release strategies. We find that when antibody is not established, the dynamics of competition depends on the respective basic reproduction numbers of the two viruses. If the basic reproductive ratio of budding virus is greater than that of lytic virus and one, budding virus can survive. When antibody is established for both strains but the neutralization capacities are the same for both strains, consequence of the competition also depends only on the basic reproduction numbers of the budding and lytic viruses. Using two concrete forms of the viral production functions, we are also able to conclude that budding virus will outcompete if the rates of viral production, death rates of infected cells and neutralizing capacities of the antibodies are the same for budding and lytic viruses. In this case, budding strategy would have evolutionary advantage. However, if the antibody neutralization capacity for the budding virus is larger than that for the lytic virus, lytic virus can outcompete provided that its reproductive ratio is very high. An explicit threshold is derived. Secondly, we consider model containing two modes for viral infection and spread, one is the diffusion-limited free virus transmission and the other is the direct cell-to-cell transfer of viral particles. By incorporating infection age, a rigorous analysis of the model shows that the model demonstrates a global threshold dynamics, fully described by the basic reproduction number, which is identified explicitly. The formula for the basic reproduction number of our model reveals the effects of various model parameters including the transmission rates of the two modes, and the impact of the infection age. We show that basic reproduction number is underestimated in the existing models that only consider the cell-free virus transmission, or the cell-to-cell infection, ignoring the other. Assuming logistic growth for target cells, we find that if the basic reproduction number is greater than one, the infection can persist and Hopf bifurcation can occur from the positive equilibrium within certain parameter ranges. Thirdly, the repulsion effect of superinfecting virion by infected cells is studied by a reaction diffusion equation model for virus infection dynamics. In this model, the diffusion of virus depends not only on its concentration gradient but also on the concentration of infected cells. The basic reproduction number, linear stability of steady states, spreading speed, and existence of traveling wave solutions for the model are discussed. It is shown that viruses spread more rapidly with the repulsion effect of infected cells on superinfecting virions, than with random diffusion only. For our model, the spreading speed of free virus is not consistent with the minimal traveling wave speed. With our general model, numerical computations of the spreading speed shows that the repulsion of superinfecting vision promotes the spread of virus, which confirms, not only qualitatively but also quantitatively, some recent experimental results. Finally, the effect of chemotactic movement of CD8+ cytotoxic T lymphocytes (CTLs) on HIV-1 infection dynamics is studied by a reaction diffusion model with chemotaxis. Choosing a typical chemosensitive function, we find that chemoattractive movement of CTLs due to HIV infection does not change stability of the positive steady state of the model. However, chemorepulsion movement of CTLs destabilizes the positive steady state as the strength of the chemotactic sensitivity increases. In this case, Turing instability occurs, which can be Hopf bifurcation or steady state bifurcation, and spatial heterogeneous patterns may form

Dynamical Models of Biology and Medicine

Mathematical and computational modeling approaches in biological and medical research are experiencing rapid growth globally. This Special Issue Book intends to scratch the surface of this exciting phenomenon. The subject areas covered involve general mathematical methods and their applications in biology and medicine, with an emphasis on work related to mathematical and computational modeling of the complex dynamics observed in biological and medical research. Fourteen rigorously reviewed papers were included in this Special Issue. These papers cover several timely topics relating to classical population biology, fundamental biology, and modern medicine. While the authors of these papers dealt with very different modeling questions, they were all motivated by specific applications in biology and medicine and employed innovative mathematical and computational methods to study the complex dynamics of their models. We hope that these papers detail case studies that will inspire many additional mathematical modeling efforts in biology and medicin