1
|
Das S, Srivastava PK, Biswas P. Exploring Hopf-bifurcations and endemic bubbles in a tuberculosis model with behavioral changes and treatment saturation. CHAOS (WOODBURY, N.Y.) 2024; 34:013126. [PMID: 38252782 DOI: 10.1063/5.0179351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2023] [Accepted: 12/22/2023] [Indexed: 01/24/2024]
Abstract
To manage risks and minimize the transmission of contagious diseases, individuals may reduce their contact with each other and take other precautions as much as possible in their daily lives and workplaces. As a result, the transmission of the infection reduces due to the behavioral changes. These behavioral changes are incorporated into models by introducing saturation in disease incidence. In this article, we propose and analyze a tuberculosis model that incorporates saturated exogenous reinfection and treatment. The stability analysis of the model's steady states is rigorously examined. We observe that the disease-free equilibrium point and the endemic equilibrium point (EEP) are globally asymptotically stable if the basic reproduction number (R0) is less than 1 and greater than 1, respectively, only when exogenous reinfection is not present (p=0) and when treatment is available for all (ω=0). However, even when R0 is less than 1, tuberculosis may persist at a specific level in the presence of exogenous reinfection and treatment saturation, leading to a backward bifurcation in the system. The existence and direction of Hopf-bifurcations are also discussed. Furthermore, we numerically validate our analytical results using different parameter sets. In the numerical examples, we study Hopf-bifurcations for parameters such as β, p, α, and ω. In one example, we observe that increasing β leads to the loss of stability of the unique EEP through a forward Hopf-bifurcation. If β is further increased, the unique EEP restores its stability, and the bifurcation diagram exhibits an interesting structure known as an endemic bubble. The existence of an endemic bubble for the saturation constant ω is also observed.
Collapse
Affiliation(s)
- Saduri Das
- National Institute of Technology Silchar, Silchar 788010, Assam, India
| | | | - Pankaj Biswas
- National Institute of Technology Silchar, Silchar 788010, Assam, India
| |
Collapse
|
2
|
Saha P, Biswas SK, Biswas MHA, Ghosh U. An SEQAIHR model to study COVID-19 transmission and optimal control strategies in Hong Kong, 2022. NONLINEAR DYNAMICS 2023; 111:6873-6893. [PMID: 36644569 PMCID: PMC9825089 DOI: 10.1007/s11071-022-08181-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Accepted: 12/09/2022] [Indexed: 06/17/2023]
Abstract
During the COVID-19 pandemic, one of the major concerns was a medical emergency in human society. Therefore it was necessary to control or restrict the disease spreading among populations in any fruitful way at that time. To frame out a proper policy for controlling COVID-19 spreading with limited medical facilities, here we propose an SEQAIHR model having saturated treatment. We check biological feasibility of model solutions and compute the basic reproduction number ( R 0 ). Moreover, the model exhibits transcritical, backward bifurcation and forward bifurcation with hysteresis with respect to different parameters under some restrictions. Further to validate the model, we fit it with real COVID-19 infected data of Hong Kong from 19th December, 2021 to 3rd April, 2022 and estimate model parameters. Applying sensitivity analysis, we find out the most sensitive parameters that have an effect on R 0 . We estimate R 0 using actual initial growth data of COVID-19 and calculate effective reproduction number for same period. Finally, an optimal control problem has been proposed considering effective vaccination and saturated treatment for hospitalized class to decrease density of the infected class and to minimize implemented cost.
Collapse
Affiliation(s)
- Pritam Saha
- Department of Applied Mathematics, University of Calcutta, Kolkata, 700009 India
| | | | | | - Uttam Ghosh
- Department of Applied Mathematics, University of Calcutta, Kolkata, 700009 India
| |
Collapse
|
3
|
Srivastava A, Sonu, Srivastava PK. Nonlinear dynamics of a SIRI model incorporating the impact of information and saturated treatment with optimal control. EUROPEAN PHYSICAL JOURNAL PLUS 2022; 137:1028. [PMID: 36106085 PMCID: PMC9462650 DOI: 10.1140/epjp/s13360-022-03201-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 08/12/2022] [Indexed: 06/15/2023]
Abstract
In this article, we propose and analyze an infectious disease model with reinfection and investigate disease dynamics by incorporating saturated treatment and information effect. In the model, we consider the case where an individual's immunity deteriorates and they become infected again after recovering. According to our findings, multiple steady states and backward bifurcation may occur as a result of treatment saturation. Further, if treatment is available for all, the disease will be eradicated providedR 0 < 1 ; however, because limited medical resources caused saturation in treatment, the disease may persist even ifR 0 < 1 . The global stability of the unique endemic steady state is established using a geometric approach. We also establish certain conditions on the transmission rate for the occurrence of periodic oscillations in the model system. Among nonlinear dynamics, we show supercritical Hopf bifurcation, bi-stability, backward Hopf bifurcation, and double Hopf bifurcation. To illustrate and validate our theoretical results, we present numerical examples. We found that when disease information coverage is high, infection cases fall considerably, and the disease persists when the reinfection rate is high. We then extend our model by incorporating two time-dependent controls, namely inhibitory interventions and treatment. Using Pontryagin's maximum principle, we prove the existence of optimal control paths and find the optimal pair of controls. According to our numerical simulations, the second control is less effective than the first. Furthermore, while implementing a single intervention at a time may be effective, combining both interventions is most effective in reducing disease burden and cost.
Collapse
Affiliation(s)
- Akriti Srivastava
- Department of Mathematics, Indian Institute of Technology Patna, Patna, 801103 India
| | - Sonu
- Department of Mathematics, Indian Institute of Technology Patna, Patna, 801103 India
| | | |
Collapse
|
4
|
Rwezaura H, Diagne ML, Omame A, de Espindola AL, Tchuenche JM. Mathematical modeling and optimal control of SARS-CoV-2 and tuberculosis co-infection: a case study of Indonesia. MODELING EARTH SYSTEMS AND ENVIRONMENT 2022; 8:5493-5520. [PMID: 35814616 PMCID: PMC9251044 DOI: 10.1007/s40808-022-01430-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/26/2022] [Accepted: 05/16/2022] [Indexed: 01/08/2023]
Abstract
A new mathematical model incorporating epidemiological features of the co-dynamics of tuberculosis (TB) and SARS-CoV-2 is analyzed. Local asymptotic stability of the disease-free and endemic equilibria are shown for the sub-models when the respective reproduction numbers are below unity. Bifurcation analysis is carried out for the TB only sub-model, where it was shown that the sub-model undergoes forward bifurcation. The model is fitted to the cumulative confirmed daily SARS-CoV-2 cases for Indonesia from February 11, 2021 to August 26, 2021. The fitting was carried out using the fmincon optimization toolbox in MATLAB. Relevant parameters in the model are estimated from the fitting. The necessary conditions for the existence of optimal control and the optimality system for the co-infection model is established through the application of Pontryagin’s Principle. Different control strategies: face-mask usage and SARS-CoV-2 vaccination, TB prevention as well as treatment controls for both diseases are considered. Simulations results show that: (1) the strategy against incident SARS-CoV-2 infection averts about 27,878,840 new TB cases; (2) also, TB prevention and treatment controls could avert 5,397,795 new SARS-CoV-2 cases. (3) In addition, either SARS-CoV-2 or TB only control strategy greatly mitigates a significant number of new co-infection cases.
Collapse
|
5
|
BUONOMO BRUNO, DELLA MARCA ROSSELLA, SHARBAYTA SILESHISINTAYEHU. A BEHAVIORAL CHANGE MODEL TO ASSESS VACCINATION-INDUCED RELAXATION OF SOCIAL DISTANCING DURING AN EPIDEMIC. J BIOL SYST 2022. [DOI: 10.1142/s0218339022500085] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The success of mass vaccination campaigns may be jeopardized by human risky behaviors. For example, high level of vaccination coverage may induce early relaxation of social distancing. In this paper, we focus on the mutual influence between the decline in prevalence, due to the rise in the overall immunization coverage, and the consequent decrease in the compliance to social distancing measures. We consider an epidemic model where both the vaccination rate and the disease transmission rate are influenced by human behavior, which in turn depends on the current and past information about the spread of the disease. We highlight the impact of the information-related parameters on the transient and asymptotic behavior of the system that is on the early stage of the epidemic and its final outcome. Among the main results, we evidence that sustained oscillations may be triggered by the behavioral memory in the prevalence-dependent vaccination rate. However, the relaxation of social distancing may induce a switch from a cyclic regime to damped oscillations.
Collapse
Affiliation(s)
- BRUNO BUONOMO
- Department of Mathematics and Applications, University of Naples Federico II, via Cintia, I-80126 Naples, Italy
| | - ROSSELLA DELLA MARCA
- Mathematics Area, SISSA – International School for Advanced Studies, via Bonomea 265, I-34136 Trieste, Italy
| | | |
Collapse
|
6
|
Das T, Srivastava PK, Kumar A. Nonlinear dynamical behavior of an SEIR mathematical model: Effect of information and saturated treatment. CHAOS (WOODBURY, N.Y.) 2021; 31:043104. [PMID: 34251223 DOI: 10.1063/5.0039048] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2020] [Accepted: 03/15/2021] [Indexed: 06/13/2023]
Abstract
When a disease spreads in a population, individuals tend to change their behavior due to the presence of information about disease prevalence. Therefore, the infection rate is affected and incidence term in the model should be appropriately modified. In addition, a limitation of medical resources has its impact on the dynamics of the disease. In this work, we propose and analyze an Susceptible-Exposed-Infected-Recovered (SEIR) model, which accounts for the information-induced non-monotonic incidence function and saturated treatment function. The model analysis is carried out, and it is found that when R0 is below one, the disease may or may not die out due to the saturated treatment (i.e., a backward bifurcation may exist and cause multi-stability). Further, we note that in this case, disease eradication is possible if medical resources are available for all. When R0 exceeds one, there is a possibility of the existence of multiple endemic equilibria. These multiple equilibria give rise to rich and complex dynamics by showing various bifurcations and oscillations (via Hopf bifurcation). A global asymptotic stability of a unique endemic equilibrium (when it exists) is established under certain conditions. An impact of information is shown and also a sensitivity analysis of model parameters is performed. Various cases are considered numerically to provide the insight of model behavior mathematically and epidemiologically. We found that the model shows hysteresis. Our study underlines that a limitation of medical resources may cause bi(multi)-stability in the model system. Also, information plays a significant role and gives rise to a rich and complex dynamical behavior of the model.
Collapse
Affiliation(s)
- Tanuja Das
- Indian Institute of Technology Patna, Patna 801103, India
| | | | - Anuj Kumar
- Thapar Institute of Engineering and Technology, Patiala 147004, India
| |
Collapse
|
7
|
Rosales GS. Mathematical and Computational Modeling of Bacterial Infection. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11606-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
|
8
|
Ghosh I, Nadim SS, Chattopadhyay J. Zoonotic MERS-CoV transmission: modeling, backward bifurcation and optimal control analysis. NONLINEAR DYNAMICS 2021; 103:2973-2992. [PMID: 33584009 PMCID: PMC7868678 DOI: 10.1007/s11071-021-06266-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 01/27/2021] [Indexed: 05/08/2023]
Abstract
Middle East Respiratory Syndrome Coronavirus (MERS-CoV) can cause mild to severe acute respiratory illness with a high mortality rate. As of January 2020, more than 2500 cases of MERS-CoV resulting in around 860 deaths were reported globally. In the absence of neither effective treatment nor a ready-to-use vaccine, control measures can be derived from mathematical models of disease epidemiology. In this manuscript, we propose and analyze a compartmental model of zoonotic MERS-CoV transmission with two co-circulating strains. The human population is considered with eight compartments while the zoonotic camel population consist of two compartments. The expression of basic reproduction numbers are obtained for both single strain and two strain version of the proposed model. We show that the disease-free equilibrium of the system with single stain is globally asymptotically stable under some parametric conditions. We also demonstrate that both models undergo backward bifurcation phenomenon, which in turn indicates that only keeping R 0 below unity may not ensure eradication. To the best of the authors knowledge, backward bifurcation was not shown in a MERS-CoV transmission model previously. Further, we perform normalized sensitivity analysis of important model parameters with respect to basic reproduction number of the proposed model. Furthermore, we perform optimal control analysis on different combination interventions with four components namely preventive measures such as use of masks, isolation of strain-1 infected people, strain-2 infected people and infected camels. Optimal control analysis suggests that combination of preventive measures and isolation of infected camels will eventually eradicate the disease from the community.
Collapse
Affiliation(s)
- Indrajit Ghosh
- Department of Computational and Data Sciences, Indian Institute of Science, Bengalore, Karnataka 560012 India
| | - Sk Shahid Nadim
- Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, 700 108 India
| | - Joydev Chattopadhyay
- Agricultural and Ecological Research Unit, Indian Statistical Institute, Kolkata, 700 108 India
| |
Collapse
|
9
|
A T, Aggarwal R, Raj YA. A fractional order HIV-TB co-infection model in the presence of exogenous reinfection and recurrent TB. NONLINEAR DYNAMICS 2021; 104:4701-4725. [PMID: 34075277 PMCID: PMC8159726 DOI: 10.1007/s11071-021-06518-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Accepted: 05/04/2021] [Indexed: 05/13/2023]
Abstract
In this article, a novel fractional order model has been introduced in Caputo sense for HIV-TB co-infection in the presence of exogenous reinfection and recurrent TB along with the treatment for both HIV and TB. The main aim of considering the fractional order model is to incorporate the memory effect of both diseases. We have analyzed both sub-models separately with fractional order. The basic reproduction number, which measures the contagiousness of the disease, is determined. The HIV sub-model is shown to have a locally asymptotically stable disease-free equilibrium point when the corresponding reproduction number, R H , is less than unity, whereas, for R H > 1 , the endemic equilibrium point comes into existence. For the TB sub-model, the disease-free equilibrium point has been proved to be locally asymptotically stable for R T < 1 . The existence of TB endemic equilibrium points in the presence of reinfection and recurrent TB for R T < 1 justifies the existence of backward bifurcation under certain restrictions on the parameters. Further, we numerically simulate the fractional order model to verify the analytical results and highlight the role of fractional order in co-infection modeling. The fractional order derivative is shown to have a crucial role in determining the transmission dynamics of HIV-TB co-infection. It is concluded that the memory effect plays a significant role in reducing the infection prevalence of HIV-TB co-infection. An increment in the number of recovered individuals can also be observed when the memory effect is taken into consideration by introducing fractional order model.
Collapse
Affiliation(s)
- Tanvi A
- Department of Mathematics, Deshbandhu College, University of Delhi, New Delhi, 110019 India
| | - Rajiv Aggarwal
- Department of Mathematics, Deshbandhu College, University of Delhi, New Delhi, 110019 India
| | - Yashi A. Raj
- Department of Mathematics, Shaheed Rajguru College of Applied Sciences for Women, University of Delhi, New Delhi, 110096 India
| |
Collapse
|
10
|
Jafari M, Kheiri H, Jabbari A. Backward bifurcation in a fractional-order and two-patch model of tuberculosis epidemic with incomplete treatment. INT J BIOMATH 2020. [DOI: 10.1142/s1793524521500078] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this paper, we consider a tuberculosis model with incomplete treatment and extend the model to a Caputo fractional-order and two-patch version with exogenous re-infection among the treated individuals, in which only susceptible individuals can travel freely between the patches. The model has multiple equilibria. We determine conditions that lead to the appearance of a backward bifurcation. The results show that the TB model can have exogenous reinfection among the treated individuals and, at the same time, does not exhibit backward bifurcation. Also, conditions that lead to the global asymptotic stability of the disease-free equilibrium are obtained. In case without reinfection, the model has four equilibria. In this case, the global asymptotic stability of the equilibria is established using the Lyapunov function theory together with the LaSalle invariance principle for fractional differential equations (FDEs). Numerical simulations confirm the validity of the theoretical results.
Collapse
Affiliation(s)
- Mohsen Jafari
- Faculty of Mathematical Sciences, University of Tabriz, Tabriz, Iran
| | - Hossein Kheiri
- Faculty of Mathematical Sciences, University of Tabriz, Tabriz, Iran
| | - Azizeh Jabbari
- Marand Faculty of Engineering, University of Tabriz, Tabriz, Iran
| |
Collapse
|
11
|
Wangari IM, Trauer J, Stone L. Modelling heterogeneity in host susceptibility to tuberculosis and its effect on public health interventions. PLoS One 2018; 13:e0206603. [PMID: 30427891 PMCID: PMC6235601 DOI: 10.1371/journal.pone.0206603] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2018] [Accepted: 10/16/2018] [Indexed: 11/25/2022] Open
Abstract
A tuberculosis (TB) model that accounts for heterogeneity in host susceptibility to tuberculosis is proposed, with the aim of investigating the implications this may have for the effectiveness of public health interventions. The model examines the possibility that recovered individuals treated from active TB and individuals treated with preventive therapy acquire different levels of immunity. This contrasts with recent studies that assume the two cohorts acquire the same level of immunity, and therefore both groups are reinfected at the same rate. The analysis presented here examines the impact of this assumption when designing intervention strategies. Comparison of reinfection rates between cohorts treated with preventive therapy and recovered individuals who were previously treated for active TB provides important epidemiological insights. It is found that the reinfection rate of the cohort treated with preventive therapy is the one that plays the key role in qualitative changes in TB dynamics. By contrast, the reinfection rate of recovered individuals (previously treated from active TB) plays a minor role. Moreover, the study shows that preventive treatment of individuals during early latency is always beneficial regardless of the level of susceptibility to reinfection. Further, if patients have greater immunity following treatment for late latent infection, then treatment is again beneficial. However, if susceptibility increases following treatment for late latent infection, the effect of treatment depends on the epidemiological setting. That is: (i) in (very) low burden settings, the effect on reactivation predominates and the burden declines with treatment; (ii) in moderate to high burden settings the effect of reinfection predominates and burden increases with treatment. The effect is most dominant between the two reinfection thresholds, RT2 and RT1, respectively associated with individuals being treated with preventive therapy and individuals with untreated late latent TB infection.
Collapse
Affiliation(s)
- Isaac Mwangi Wangari
- Mathematical Sciences, School of Science, Royal Melbourne Institute of Technology, Melbourne, Victoria, Australia
- * E-mail:
| | - James Trauer
- School of Public Health and Preventive Medicine, Monash University, Melbourne 3004, Australia
| | - Lewi Stone
- Mathematical Sciences, School of Science, Royal Melbourne Institute of Technology, Melbourne, Victoria, Australia
- Biomathematics Unit, Department of Zoology, Faculty of Life Sciences, Tel Aviv University, Israel
| |
Collapse
|