1
|
Deng Q, Guo T, Qiu Z, Chen Y. A mathematical model for HIV dynamics with multiple infections: implications for immune escape. J Math Biol 2024; 89:6. [PMID: 38762831 DOI: 10.1007/s00285-024-02104-w] [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: 08/26/2023] [Revised: 12/15/2023] [Accepted: 04/25/2024] [Indexed: 05/20/2024]
Abstract
Multiple infections enable the recombination of different strains, which may contribute to viral diversity. How multiple infections affect the competition dynamics between the two types of strains, the wild and the immune escape mutant, remains poorly understood. This study develops a novel mathematical model that includes the two strains, two modes of viral infection, and multiple infections. For the representative double-infection case, the reproductive numbers are derived and global stabilities of equilibria are obtained via the Lyapunov direct method and theory of limiting systems. Numerical simulations indicate similar viral dynamics regardless of multiplicities of infections though the competition between the two strains would be the fiercest in the case of quadruple infections. Through sensitivity analysis, we evaluate the effect of parameters on the set-point viral loads in the presence and absence of multiple infections. The model with multiple infections predict that there exists a threshold for cytotoxic T lymphocytes (CTLs) to minimize the overall viral load. Weak or strong CTLs immune response can result in high overall viral load. If the strength of CTLs maintains at an intermediate level, the fitness cost of the mutant is likely to have a significant impact on the evolutionary dynamics of mutant viruses. We further investigate how multiple infections alter the viral dynamics during the combination antiretroviral therapy (cART). The results show that viral loads may be underestimated during cART if multiple-infection is not taken into account.
Collapse
Affiliation(s)
- Qi Deng
- School of Mathematics and Statistics, Nanjing University of Science and Technology, Nanjing, 210094, People's Republic of China
- Department of Mathematics, Wilfrid Laurier University, Waterloo, N2L 3C5, Canada
| | - Ting Guo
- Aliyun School of Big Data, Changzhou University, Changzhou, 213164, People's Republic of China
| | - Zhipeng Qiu
- School of Mathematics and Statistics, Nanjing University of Science and Technology, Nanjing, 210094, People's Republic of China
| | - Yuming Chen
- Department of Mathematics, Wilfrid Laurier University, Waterloo, N2L 3C5, Canada.
| |
Collapse
|
2
|
Guo T, Deng Q, Qiu Z, Rong L. HIV infection dynamics and viral rebound: Modeling results from humanized mice. J Theor Biol 2023; 567:111490. [PMID: 37054969 DOI: 10.1016/j.jtbi.2023.111490] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Revised: 03/22/2023] [Accepted: 04/05/2023] [Indexed: 04/15/2023]
Abstract
Despite years of combined antiretroviral therapy (cART), HIV persists in infected individuals. The virus also rebounds after the cessation of cART. The sources contributing to viral persistence and rebound are not fully understood. When viral rebound occurs, what affects the time to rebound and how to delay the rebound remain unclear. In this paper, we started with the data fitting of an HIV infection model to the viral load data in treated and untreated humanized myeloid-only mice (MoM) in which macrophages serve as the target of HIV infection. By fixing the parameter values for macrophages from the MoM fitting, we fit a mathematical model including the infection of two target cell populations to the viral load data from humanized bone marrow/liver/thymus (BLT) mice, in which both CD4+ T cells and macrophages are the target of HIV infection. Data fitting suggests that the viral load decay in BLT mice under treatment has three phases. The loss of infected CD4+ T cells and macrophages is a major contributor to the first two phases of viral decay, and the last phase may be due to the latent infection of CD4+ T cells. Numerical simulations using parameter estimates from the data fitting show that the pre-ART viral load and the latent reservoir size at treatment cessation can affect viral growth rate and predict the time to viral rebound. Model simulations further reveal that early and prolonged cART can delay the viral rebound after cessation of treatment, which may have implications in the search for functional control of HIV infection.
Collapse
Affiliation(s)
- Ting Guo
- Aliyun School of Big Data, Changzhou University, Changzhou, 213164, China
| | - Qi Deng
- School of Science, Nanjing University of Science and Technology, Nanjing 210094, China
| | - Zhipeng Qiu
- Center for Basic Teaching and Experiment, Nanjing University of Science and Technology, Jiangyin 214443, China
| | - Libin Rong
- Department of Mathematics, University of Florida, Gainesville, FL 32611, USA.
| |
Collapse
|
3
|
Mathematical Analysis of Epidemic Models with Treatment in Heterogeneous Networks. Bull Math Biol 2023; 85:11. [PMID: 36602636 PMCID: PMC9813916 DOI: 10.1007/s11538-022-01116-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Accepted: 12/26/2022] [Indexed: 01/06/2023]
Abstract
In this paper, we formulate two different network-based epidemic models to investigate the effect of partly effective treatment on disease dynamics. The first network model represents the individuals with heterogeneous number of contacts in a population as choosing a new partner at each moment, whereas the second one assumes the individuals have fixed or stable neighbors. The basic reproduction number [Formula: see text] is computed for each model, using the next generation matrix method. In particular, the critical treatment rate is defined for the model, above which the disease can be eliminated through the treatment. The final epidemic size relations are derived, and the solvability of these implicit equations is studied. In particular, a unique solution of the implicit equation for the final epidemic size is determined, and by rewriting the implicit equation as a suitable fixed point problem, it is proved that the iteration of the fixed point problem converges to the unique solution. Stochastic simulations and numerical simulations, including in comparison with the model outputs and the joint influence of network topology and treatment on the final epidemic size, are conducted to illustrate the theoretical results.
Collapse
|
4
|
Towards a new combination therapy with vectored immunoprophylaxis for HIV: Modeling "shock and kill" strategy. Math Biosci 2023; 355:108954. [PMID: 36525996 DOI: 10.1016/j.mbs.2022.108954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 09/23/2022] [Accepted: 12/06/2022] [Indexed: 12/15/2022]
Abstract
Latently infected cells are considered as a major barrier to curing Human Immunodeficiency Virus (HIV) infection. Reactivation of latently infected cells followed by killing the actively infected cells may be a potential strategy ("shock and kill") to purge the latent reservoir. Based on vectored immunoprophylaxis (VIP) experiment that can elicit bNAbs, in this paper a mathematical model is formulated to explore the efficacy of "shock and kill" strategy with VIP. We derive the basic reproduction number R0 of the model and show that R0 completely determines the dynamics of the model: if R0<1, the disease-free equilibrium is globally asymptotically stable; if R0>1, the system is uniformly persistent. Numerical simulations suggest that the "shock and kill" strategy with VIP can effectively control HIV infection while this strategy cannot eradicate the reservoir without VIP although it can alleviate the HIV infection. To model the administration of drugs and vaccine more realistically, pharmacokinetics and pulse vaccination are incorporated into the model of ordinary differential equations. The resultants are described by impulsive differential equations. The thresholds are obtained for the frequency and strength of the vaccination to eliminate the viruses. Furthermore, the most appropriate times are numerically investigated for starting a short-term latency-reversing agents (LRAs) treatment relative to ART considering the toxicity of LRAs. The results show that LRAs treatment at the beginning of ART might be a better option. These results have important implications for the design of HIV cure-related clinical trials.
Collapse
|
5
|
Salehipour A, Bagheri M, Sabahi M, Dolatshahi M, Boche D. Combination Therapy in Alzheimer’s Disease: Is It Time? J Alzheimers Dis 2022; 87:1433-1449. [DOI: 10.3233/jad-215680] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Alzheimer’s disease (AD) is the most common cause of dementia globally. There is increasing evidence showing AD has no single pathogenic mechanism, and thus treatment approaches focusing only on one mechanism are unlikely to be meaningfully effective. With only one potentially disease modifying treatment approved, targeting amyloid-β (Aβ), AD is underserved regarding effective drug treatments. Combining multiple drugs or designing treatments that target multiple pathways could be an effective therapeutic approach. Considering the distinction between added and combination therapies, one can conclude that most trials fall under the category of added therapies. For combination therapy to have an actual impact on the course of AD, it is likely necessary to target multiple mechanisms including but not limited to Aβ and tau pathology. Several challenges have to be addressed regarding combination therapy, including choosing the correct agents, the best time and stage of AD to intervene, designing and providing proper protocols for clinical trials. This can be achieved by a cooperation between the pharmaceutical industry, academia, private research centers, philanthropic institutions, and the regulatory bodies. Based on all the available information, the success of combination therapy to tackle complicated disorders such as cancer, and the blueprint already laid out on how to implement combination therapy and overcome its challenges, an argument can be made that the field has to move cautiously but quickly toward designing new clinical trials, further exploring the pathological mechanisms of AD, and re-examining the previous studies with combination therapies so that effective treatments for AD may be finally found.
Collapse
Affiliation(s)
- Arash Salehipour
- Neurosurgery Research Group (NRG), Student Research Committee, Hamadan University of Medical Sciences, Hamadan, Iran
- NeuroImaging Network (NIN), Universal Scientific Education and Research Network (USERN), Tehran, Iran
| | - Motahareh Bagheri
- Neurosurgery Research Group (NRG), Student Research Committee, Hamadan University of Medical Sciences, Hamadan, Iran
| | - Mohammadmahdi Sabahi
- Neurosurgery Research Group (NRG), Student Research Committee, Hamadan University of Medical Sciences, Hamadan, Iran
- NeuroImaging Network (NIN), Universal Scientific Education and Research Network (USERN), Tehran, Iran
| | - Mahsa Dolatshahi
- NeuroImaging Network (NIN), Universal Scientific Education and Research Network (USERN), Tehran, Iran
- Students’ Scientific Research Center (SSRC), Tehran University of Medical Sciences, Tehran, Iran
| | - Delphine Boche
- Clinical Neurosciences, Clinical and Experimental Sciences, Faculty of Medicine, University of Southampton, United Kingdom
| |
Collapse
|
6
|
Optimal Control of an HIV Model with Gene Therapy and Latency Reversing Agents. MATHEMATICAL AND COMPUTATIONAL APPLICATIONS 2021. [DOI: 10.3390/mca26040077] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In this paper, we study the dynamics of HIV under gene therapy and latency reversing agents. While previous works modeled either the use of gene therapy or latency reversing agents, we consider the effects of a combination treatment strategy. For constant treatment controls, we establish global stability of the disease-free equilibrium and endemic equilibrium based on the value of R0. We then consider time-dependent controls and formulate an associated optimal control problem that emphasizes reduction of the latent reservoir. Characterizations for the optimal control profiles are found using Pontryagin’s Maximum Principle. We perform numerical simulations of the optimal control model using the fourth-order Runge–Kutta forward-backward sweep method. We find that a combination treatment of gene therapy with latency reversing agents provides better remission times than gene therapy alone. We conclude with a discussion of our findings and future work.
Collapse
|
7
|
Bai N, Xu R. Mathematical analysis of an HIV model with latent reservoir, delayed CTL immune response and immune impairment. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:1689-1707. [PMID: 33757205 DOI: 10.3934/mbe.2021087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In this paper, an in-host HIV infection model with latent reservoir, delayed CTL immune response and immune impairment is investigated. By using suitable Lyapunov functions and LaSalle's invariance principle, it is shown that when time delay is equal to zero, the immunity-inactivated reproduction ratio is a threshold determining the global dynamics of the model. By means of the persistence theory for infinite dimensional systems, it is proven that if the immunity-inactivated reproduction ratio is greater than unity, the model is permanent. Choosing time delay as the bifurcation parameter and analyzing the corresponding characteristic equation of the linearized system, the existence of a Hopf bifurcation at the immunity-activated equilibrium is established. Numerical simulations are carried out to illustrate the theoretical results and reveal the effects of some key parameters on viral dynamics.
Collapse
Affiliation(s)
- Ning Bai
- Complex Systems Research Center, Shanxi University, Taiyuan 030006, China
| | - Rui Xu
- Complex Systems Research Center, Shanxi University, Taiyuan 030006, China
| |
Collapse
|
8
|
Mutua JM, Wang FB, Vaidya NK. Effects of periodic intake of drugs of abuse (morphine) on HIV dynamics: Mathematical model and analysis. Math Biosci 2020; 326:108395. [PMID: 32485213 DOI: 10.1016/j.mbs.2020.108395] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2019] [Revised: 05/27/2020] [Accepted: 05/27/2020] [Indexed: 11/27/2022]
Abstract
Drugs of abuse, such as opiates, have been widely associated with diminishing host-immune responses, including suppression of HIV-specific antibody responses. In particular, periodic intake of the drugs of abuse can result in time-varying periodic antibody level within HIV-infected individuals, consequently altering the HIV dynamics. In this study, we develop a mathematical model to analyze the effects of periodic intake of morphine, a widely used opiate. We consider two routes of morphine intake, namely, intravenous morphine (IVM) and slow-release oral morphine (SROM), and integrate several morphine pharmacodynamic parameters into HIV dynamics model. Using our non-autonomous model system we formulate the infection threshold, Ri, for global stability of infection-free equilibrium, which provides a condition for avoiding viral infection in a host. We demonstrate that the infection threshold highly depends on the morphine pharmacodynamic parameters. Such information can be useful in the design of antibody-based vaccines. In addition, we also thoroughly evaluate how alteration of the antibody level due to periodic intake of morphine can affect the viral load and the CD4 count in HIV infected drug abusers.
Collapse
Affiliation(s)
- Jones M Mutua
- Department of Computer Science, Mathematics, & Physics, Missouri Western State University, St. Joseph, MO, USA
| | - Feng-Bin Wang
- Department of Natural Science in the Center for General Education, Chang Gung University, Guishan, Taoyuan 333, Taiwan; Community Medicine Research Center, Chang Gung Memorial Hospital, Keelung Branch, Keelung 204, Taiwan
| | - Naveen K Vaidya
- Department of Mathematics and Statistics, San Diego State University, San Diego, CA, USA; Computational Science Research Center, San Diego State University, San Diego, CA, USA; Viral Information Institute, San Diego State University, San Diego, CA, USA.
| |
Collapse
|