1
|
Kim YR, Min Y, Okogun-Odompley JN. A mathematical model of COVID-19 with multiple variants of the virus under optimal control in Ghana. PLoS One 2024; 19:e0303791. [PMID: 38954691 PMCID: PMC11218976 DOI: 10.1371/journal.pone.0303791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 04/30/2024] [Indexed: 07/04/2024] Open
Abstract
In this paper, we suggest a mathematical model of COVID-19 with multiple variants of the virus under optimal control. Mathematical modeling has been used to gain deeper insights into the transmission of COVID-19, and various prevention and control strategies have been implemented to mitigate its spread. Our model is a SEIR-based model for multi-strains of COVID-19 with 7 compartments. We also consider the circulatory structure to account for the termination of immunity for COVID-19. The model is established in terms of the positivity and boundedness of the solution and the existence of equilibrium points, and the local stability of the solution. As a result of fitting data of COVID-19 in Ghana to the model, the basic reproduction number of the original virus and Delta variant was estimated to be 1.9396, and the basic reproduction number of the Omicron variant was estimated to be 3.4905, which is 1.8 times larger than that. We observe that even small differences in the incubation and recovery periods of two strains with the same initial transmission rate resulted in large differences in the number of infected individuals. In the case of COVID-19, infections caused by the Omicron variant occur 1.5 to 10 times more than those caused by the original virus. In terms of the optimal control strategy, we formulate three control strategies focusing on social distancing, vaccination, and testing-treatment. We have developed an optimal control model for the three strategies outlined above for the multi-strain model using the Pontryagin's Maximum Principle. Through numerical simulations, we analyze three optimal control strategies for each strain and also consider combinations of the two control strategies. As a result of the simulation, all control strategies are effective in reducing disease spread, in particular, vaccination strategies are more effective than the other two control strategies. In addition the combination of the two strategies also reduces the number of infected individuals by 1/10 compared to implementing one strategy, even when mild levels are implemented. Finally, we show that if the testing-treatment strategy is not properly implemented, the number of asymptomatic and unidentified infections may surge. These results could help guide the level of government intervention and prevention strategy formulation.
Collapse
Affiliation(s)
- Young Rock Kim
- Major in Mathematics Education, Graduate School of Education, Hankuk University of Foreign Studies, Seoul, Republic of Korea
| | - Youngho Min
- Major in Mathematics Education, Graduate School of Education, Hankuk University of Foreign Studies, Seoul, Republic of Korea
| | | |
Collapse
|
2
|
Ciupe SM, Dahari H, Ploss A. Mathematical Models of Early Hepatitis B Virus Dynamics in Humanized Mice. Bull Math Biol 2024; 86:53. [PMID: 38594319 PMCID: PMC11003933 DOI: 10.1007/s11538-024-01284-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 03/15/2024] [Indexed: 04/11/2024]
Abstract
Analyzing the impact of the adaptive immune response during acute hepatitis B virus (HBV) infection is essential for understanding disease progression and control. Here we developed mathematical models of HBV infection which either lack terms for adaptive immune responses, or assume adaptive immune responses in the form of cytolytic immune killing, non-cytolytic immune cure, or non-cytolytic-mediated block of viral production. We validated the model that does not include immune responses against temporal serum hepatitis B DNA (sHBV) and temporal serum hepatitis B surface-antigen (HBsAg) experimental data from mice engrafted with human hepatocytes (HEP). Moreover, we validated the immune models against sHBV and HBsAg experimental data from mice engrafted with HEP and human immune system (HEP/HIS). As expected, the model that does not include adaptive immune responses matches the observed high sHBV and HBsAg concentrations in all HEP mice. By contrast, while all immune response models predict reduction in sHBV and HBsAg concentrations in HEP/HIS mice, the Akaike Information Criterion cannot discriminate between non-cytolytic cure (resulting in a class of cells refractory to reinfection) and antiviral block functions (of up to 99 % viral production 1-3 weeks following peak viral load). We can, however, reject cytolytic killing, as it can only match the sHBV and HBsAg data when we predict unrealistic levels of hepatocyte loss.
Collapse
Affiliation(s)
- Stanca M Ciupe
- Department of Mathematics, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA.
| | - Harel Dahari
- Division of Hepatology, Department of Medicine, Loyola University, Chicago, IL, USA
| | - Alexander Ploss
- Department of Molecular Biology, Princeton University, Princeton, NJ, USA
| |
Collapse
|
3
|
Zhang S, Li F, Xu X. Dynamics and control strategy for a delayed viral infection model. JOURNAL OF BIOLOGICAL DYNAMICS 2022; 16:44-63. [PMID: 35075973 DOI: 10.1080/17513758.2022.2028024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2020] [Accepted: 01/04/2022] [Indexed: 06/14/2023]
Abstract
In this paper, we derive a delayed epidemic model to describe the characterization of cytotoxic T lymphocyte (CTL)-mediated immune response against virus infection. The stability of equilibria and the existence of Hopf bifurcation are analysed. Theoretical results reveal that if the basic reproductive number is greater than 1, the positive equilibrium may lose its stability and the bifurcated periodic solution occurs when time delay is taken as the bifurcation parameter. Furthermore, we investigate an optimal control problem according to the delayed model based on the available therapy for hepatitis B infection. With the aim of minimizing the infected cells and viral load with consideration for the treatment costs, the optimal solution is discussed analytically. For the case when periodic solution occurs, numerical simulations are performed to suggest the optimal therapeutic strategy and compare the model-predicted consequences.
Collapse
Affiliation(s)
- Suxia Zhang
- School of Science, Xi'an University of Technology, Xi'an, People's Republic of China
| | - Fei Li
- School of Science, Xi'an University of Technology, Xi'an, People's Republic of China
| | - Xiaxia Xu
- School of Science, Xi'an University of Technology, Xi'an, People's Republic of China
| |
Collapse
|
4
|
de los Reyes AA, Kim Y. Optimal regulation of tumour-associated neutrophils in cancer progression. ROYAL SOCIETY OPEN SCIENCE 2022; 9:210705. [PMID: 35127110 PMCID: PMC8808100 DOI: 10.1098/rsos.210705] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 11/19/2021] [Indexed: 06/14/2023]
Abstract
In a tumour microenvironment, tumour-associated neutrophils could display two opposing differential phenotypes: anti-tumour (N1) and pro-tumour (N2) effector cells. Converting N2 to N1 neutrophils provides innovative therapies for cancer treatment. In this study, a mathematical model for N1-N2 dynamics describing the cancer survival and immune inhibition in response to TGF-β and IFN-β is considered. The effects of exogenous intervention of TGF-β inhibitor and IFN-β are examined in order to enhance N1 recruitment to combat tumour progression. Our approach employs optimal control theory to determine drug infusion protocols that could minimize tumour volume with least administration cost possible. Four optimal control scenarios corresponding to different therapeutic strategies are explored, namely, TGF-β inhibitor control only, IFN-β control only, concomitant TGF-β inhibitor and IFN-β controls, and alternating TGF-β inhibitor and IFN-β controls. For each scheme, different initial conditions are varied to depict different pathophysiological condition of a cancer patient, leading to adaptive treatment schedule. TGF-β inhibitor and IFN-β drug dosages, total drug amount, infusion times and relative cost of drug administrations are obtained under various circumstances. The control strategies achieved could guide in designing individualized therapeutic protocols.
Collapse
Affiliation(s)
- Aurelio A. de los Reyes
- Biomedical Mathematics Group, Pioneer Research Center for Mathematical and Computational Sciences, Institute for Basic Science, Daejeon 34126, Republic of Korea
- Institute of Mathematics, University of the Philippines Diliman, Quezon City 1101, Philippines
| | - Yangjin Kim
- Department of Mathematics, Konkuk University, Seoul 05029, Republic of Korea
- Mathematical Biosciences Institute, Columbus, OH 43210, USA
| |
Collapse
|
5
|
Optimal Allocation of Vaccine and Antiviral Drugs for Influenza Containment over Delayed Multiscale Epidemic Model considering Time-Dependent Transmission Rate. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2021; 2021:4348910. [PMID: 34707682 PMCID: PMC8545521 DOI: 10.1155/2021/4348910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Revised: 07/19/2021] [Accepted: 08/26/2021] [Indexed: 11/17/2022]
Abstract
In this study, two types of epidemiological models called "within host" and "between hosts" have been studied. The within-host model represents the innate immune response, and the between-hosts model signifies the SEIR (susceptible, exposed, infected, and recovered) epidemic model. The major contribution of this paper is to break the chain of infectious disease transmission by reducing the number of susceptible and infected people via transferring them to the recovered people group with vaccination and antiviral treatment, respectively. Both transfers are considered with time delay. In the first step, optimal control theory is applied to calculate the optimal final time to control the disease within a host's body with a cost function. To this end, the vaccination that represents the effort that converts healthy cells into resistant-to-infection cells in the susceptible individual's body is used as the first control input to vaccinate the susceptible individual against the disease. Moreover, the next control input (antiviral treatment) is applied to eradicate the concentrations of the virus and convert healthy cells into resistant-to-infection cells simultaneously in the infected person's body to treat the infected individual. The calculated optimal time in the first step is considered as the delay of vaccination and antiviral treatment in the SEIR dynamic model. Using Pontryagin's maximum principle in the second step, an optimal control strategy is also applied to an SEIR mathematical model with a nonlinear transmission rate and time delay, which is computed as optimal time in the first step. Numerical results are consistent with the analytical ones and corroborate our theoretical results.
Collapse
|
6
|
A quantitative systems pharmacology model for acute viral hepatitis B. Comput Struct Biotechnol J 2021; 19:4997-5007. [PMID: 34589180 PMCID: PMC8449028 DOI: 10.1016/j.csbj.2021.08.052] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Revised: 08/31/2021] [Accepted: 08/31/2021] [Indexed: 12/25/2022] Open
Abstract
Mechanistic model characterizing acute immune response and HBV system interactions. Key role of the cellular and regulatory response triggering hepatitis B chronicity. Modelling framework to easily incorporate and explore additional biological mechanisms.
Hepatitis B liver infection is caused by hepatitis B virus (HBV) and represents a major global disease problem when it becomes chronic, as is the case for 80–90% of vertical or early life infections. However, in the vast majority (>95%) of adult exposures, the infected individuals are capable of mounting an effective immune response leading to infection resolution. A good understanding of HBV dynamics and the interaction between the virus and immune system during acute infection represents an essential step to characterize and understand the key biological processes involved in disease resolution, which may help to identify potential interventions to prevent chronic hepatitis B. In this work, a quantitative systems pharmacology model for acute hepatitis B characterizing viral dynamics and the main components of the innate, adaptive, and tolerant immune response has been successfully developed. To do so, information from multiple sources and across different organization levels has been integrated in a common mechanistic framework. The final model adequately describes the chronology and plausibility of an HBV-triggered immune response, as well as clinical data from acute patients reported in the literature. Given the holistic nature of the framework, the model can be used to illustrate the relevance of the different immune pathways and biological processes to ultimate response, observing the negligible contribution of the innate response and the key contribution of the cellular response on viral clearance. More specifically, moderate reductions of the proliferation of activated cytotoxic CD8+ lymphocytes or increased immunoregulatory effects can drive the system towards chronicity.
Collapse
Key Words
- AHB, acute hepatitis B
- ALT, alanine aminotransferase
- CHB, chronic hepatitis B
- CTL*, activated CTL
- CTL, antigen-specific cytotoxic T lymphocytes
- CTLm, memory CTL
- DC*, activated dendritic cells
- DC, dendritic cells
- HB, Hepatitis B
- HBV, hepatitis B virus, HBV DNA, circulating DNA levels of HBV
- HBsAg, hepatitis B surface antigen
- Hep, hepatocytes
- Hepatitis B
- Heptot, total hepatocytes
- IFN, interferon
- Immune system dynamics
- LN, lymph node
- LPC, long-lived plasma cells
- LV, liver
- MDSC, myeloid-derived suppressor cells
- Mechanistic modeling
- NK*, activated NK
- NK, natural killer cells
- ODE, ordinary differential equations
- PB, plasmablasts
- PC, plasma cells
- PL, plasma
- QSP, quantitative systems pharmacology
- Quantitative systems pharmacology
- SPC, short-lived plasma cells
- TRAIL, tumor necrosis factor–related apoptosis-inducing ligand
- Th0, naïve T cells
- Treg, regulatory T cells
- Viral dynamics
- anti-HBc, specific antibodies against core hepatitis B antigen
- anti-HBs, specific antibodies against surface hepatitis B antigen
- dHep, debris hepatocytes
- iHep, infected hepatocytes
- pDC, plasmacytoid DC
Collapse
|
7
|
Ciupe SM, Vaidya NK, Forde JE. Early events in hepatitis B infection: the role of inoculum dose. Proc Biol Sci 2021; 288:20202715. [PMID: 33563115 DOI: 10.1098/rspb.2020.2715] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The relationship between the inoculum dose and the ability of the pathogen to invade the host is poorly understood. Experimental studies in non-human primates infected with different inoculum doses of hepatitis B virus have shown a non-monotonic relationship between dose magnitude and infection outcome, with high and low doses leading to 100% liver infection and intermediate doses leading to less than 0.1% liver infection, corresponding to CD4 T-cell priming. Since hepatitis B clearance is CD8 T-cell mediated, the question of whether the inoculum dose influences CD8 T-cell dynamics arises. To help answer this question, we developed a mathematical model of virus-host interaction following hepatitis B virus infection. Our model explains the experimental data well, and predicts that the inoculum dose affects both the timing of the CD8 T-cell expansion and the quality of its response, especially the non-cytotoxic function. We find that a low-dose challenge leads to slow CD8 T-cell expansion, weak non-cytotoxic functions, and virus persistence; high- and medium-dose challenges lead to fast CD8 T-cell expansion, strong cytotoxic and non-cytotoxic function, and virus clearance; while a super-low-dose challenge leads to delayed CD8 T-cell expansion, strong cytotoxic and non-cytotoxic function, and virus clearance. These results are useful for designing immune cell-based interventions.
Collapse
Affiliation(s)
- Stanca M Ciupe
- Department of Mathematics, Virginia Tech, Blacksburg, 24060 VA, USA
| | - Naveen K Vaidya
- Department of Mathematics and Statistics, San Diego State University, San Diego, CA 92182, USA.,Computational Science Research Center, San Diego State University, San Diego, CA 92182, USA.,Viral Information Institute, San Diego State University, San Diego, CA 92182, USA
| | - Jonathan E Forde
- Department of Mathematics and Computer Science, Hobart and William Smith Colleges, Geneva, New York 14456, USA
| |
Collapse
|
8
|
Abstract
Hepatitis B virus infection is the cause of liver diseases such as cirrhosis and liver cancer. Understanding the host-virus mechanisms that mediate virus pathogenesis can help design better preventive measures for disease control. Mathematical models have been used alongside experimental data to provide insight into the role of immune responses during the acute and chronic hepatitis B infections as well as virus dynamics following administration of combined drug therapy. In this paper, we review several modeling studies on virus-host interactions during acute infection, the virus-host characteristics responsible for transition to chronic disease, and the efficacy and optimal control measures of drug therapy. We conclude by presenting our opinion on the future directions of the field.
Collapse
Affiliation(s)
- Stanca M Ciupe
- Department of Mathematics, Virginia Tech, Blacksburg, VA, USA
| |
Collapse
|
9
|
Ciupe SM, Heffernan JM. In-host modeling. Infect Dis Model 2017; 2:188-202. [PMID: 29928736 PMCID: PMC6001971 DOI: 10.1016/j.idm.2017.04.002] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2017] [Revised: 04/24/2017] [Accepted: 04/26/2017] [Indexed: 01/14/2023] Open
Abstract
Understanding the mechanisms governing host-pathogen kinetics is important and can guide human interventions. In-host mathematical models, together with biological data, have been used in this endeavor. In this review, we present basic models used to describe acute and chronic pathogenic infections. We highlight the power of model predictions, the role of drug therapy, and advantage of considering the dynamics of immune responses. We also present the limitations of these models due in part to the trade-off between the complexity of the model and their predictive power, and the challenges a modeler faces in determining the appropriate formulation for a given problem.
Collapse
Affiliation(s)
- Stanca M. Ciupe
- Department of Mathematics, Virginia Tech, Blacksburg, VA, USA
| | - Jane M. Heffernan
- Centre for Disease Modelling, Department of Mathematics & Statistics, York University, Toronto, ON, Canada
| |
Collapse
|