1
|
Arora G, Bairagi N, Chatterjee S. A mathematical model to study low-dose metronomic scheduling for chemotherapy. Math Biosci 2024; 372:109186. [PMID: 38580078 DOI: 10.1016/j.mbs.2024.109186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 02/21/2024] [Accepted: 03/27/2024] [Indexed: 04/07/2024]
Abstract
Metronomic chemotherapy refers to the frequent administration of chemotherapeutic agents at a lower dose and presents an attractive alternative to conventional chemotherapy with encouraging response rates. However, the schedule of the therapy, including the dosage of the drug, is usually based on empiricism. The confounding effects of tumor-endothelial-immune interactions during metronomic administration of drugs have not yet been explored in detail, resulting in an incomplete assessment of drug dose and frequency evaluations. The present study aimed to gain a mechanistic understanding of different actions of metronomic chemotherapy using a mathematical model. We have established an analytical condition for determining the dosage and frequency of the drug depending on its clearance rate for complete tumor elimination. The model also brings forward the immune-mediated clearance of the tumor during the metronomic administration of the chemotherapeutic agent. The results from the global sensitivity analysis showed an increase in the sensitivity of drug and immune-mediated killing factors toward the tumor population during metronomic scheduling. Our results emphasize metronomic scheduling over the maximum tolerated dose (MTD) and define a model-based approach for approximating the optimal schedule of drug administration to eliminate tumors while minimizing harm to the immune cells and the patient's body.
Collapse
Affiliation(s)
- Garhima Arora
- Complex Analysis Group, Translational Health Science and Technology Institute, NCR Biotech Science Cluster, Faridabad, 121001, India
| | - Nandadulal Bairagi
- Department of Mathematics, Centre for Mathematical Biology and Ecology, Jadavpur University, Kolkata, 700032, India
| | - Samrat Chatterjee
- Complex Analysis Group, Translational Health Science and Technology Institute, NCR Biotech Science Cluster, Faridabad, 121001, India.
| |
Collapse
|
2
|
Majumder A, Sardar S, Bairagi N. The effect of noise in an HIV infection model with cytotoxic T-lymphocyte impairment. CHAOS (WOODBURY, N.Y.) 2022; 32:113131. [PMID: 36456349 DOI: 10.1063/5.0105770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Accepted: 10/06/2022] [Indexed: 06/17/2023]
Abstract
The human immunodeficiency virus (HIV) interacts with the immune cells within the human body, where the environment is uncertain and noisy. Stochastic models can successfully encapsulate the effect of such a noisy environment compared to their deterministic counterparts. The human immune system is complex but well-coordinated with various immune cells like C D 4T cells, dendritic cells, and cytotoxic T-lymphocyte (CTL) cells, among many others. The CTL can kill the antigenic cells after its recognition. However, the efficacy of CTL in removing the infected C D 4T cells is progressively compromised in HIV-infected individuals. This paper considers a noise-induced HIV-immune cell interaction model with immune impairment. A multiplicative white noise is introduced in the infection rate parameter to represent the fluctuations around the average value of the rate parameter as a causative effect of the noise. We analyzed the deterministic and stochastic models and prescribed sufficient conditions for infection eradication and persistence. It is determined under what parametric restrictions the asymptotic solutions of the noise-induced system will be a limiting case of the deterministic solutions. Simulation results revealed that the solutions of the deterministic system either converge to a CTL-dominated interior equilibrium or a CTL-free immunodeficient equilibrium, depending on the initial values of the system. Stochastic analysis divulged that higher noise might be helpful in the infection removal process. The extinction time of infected C D 4T cells for some fixed immune impairment gradually decreases with increasing noise intensity and follows the power law.
Collapse
Affiliation(s)
- Abhijit Majumder
- Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata 700032, India
| | - Shibani Sardar
- Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata 700032, India
| | - Nandadulal Bairagi
- Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata 700032, India
| |
Collapse
|
3
|
Backward bifurcation in within-host HIV models. Math Biosci 2021; 335:108569. [PMID: 33636199 DOI: 10.1016/j.mbs.2021.108569] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 02/16/2021] [Accepted: 02/16/2021] [Indexed: 11/21/2022]
Abstract
The activation and proliferation of naive CD4 T cells produce helper T cells, and increase the susceptible population in the presence of HIV. This may cause backward bifurcation. To verify this, we construct a simple within-host HIV model that includes the key variables, namely healthy naive CD4 T cells, helper T cells, infected CD4 T cells and virus. When the viral basic reproduction number R0 is less than unity, we show theoretically and numerically that bistability for RC<R0<1 can be caused by a backward bifurcation due to a new susceptible population produced by activation of healthy naive CD4 T cells that become helper T cells. An extended model including the CTL dynamics may also show this backward bifurcation. In the case that the homeostatic source of healthy naive CD4 T cells is large, RC is approximately the threshold for HIV to persist independent of initial conditions. The backward bifurcation may still occur even when we consider latent infections of naive CD4 T cells. Thus to control the spread of within-host HIV, it may be necessary for treatment to reduce the reproduction number below RC.
Collapse
|
4
|
Abstract
The interplay between immune response and HIV is intensely studied via mathematical modeling, with significant insights but few direct answers. In this short review, we highlight advances and knowledge gaps across different aspects of immunity. In particular, we identify the innate immune response and its role in priming the adaptive response as ripe for modeling. The latter have been the focus of most modeling studies, but we also synthesize key outstanding questions regarding effector mechanisms of cellular immunity and development of broadly neutralizing antibodies. Thus far, most modeling studies aimed to infer general immune mechanisms; we foresee that significant progress will be made next by detailed quantitative fitting of models to data, and prediction of immune responses.
Collapse
Affiliation(s)
- Jessica M Conway
- Department of Mathematics and Center for Infectious Disease Dynamics, Pennsylvania State University, University Park PA 16802, USA
| | - Ruy M Ribeiro
- Laboratorio de Biomatematica, Faculdade de Medicina da Universidade de Lisboa, Portugal and Theoretical Biology and Biophysics, Los Alamos National Laboratory, Los Alamos, NM 87545, USA
| |
Collapse
|
5
|
Abstract
Some HIV-infected patients (the so-called post-treatment controllers) can control the virus after cessation of antiretroviral therapy. A small fraction of patients can even naturally maintain undetectable viral load without therapy (they are called elite controllers). The immune response may play an important role in viral control in these patients. In this paper, we analyze a within-host model including immune response to study the virus dynamics in HIV-infected patients. We derived two threshold values for the immune cell proliferation parameter. Below the lower immune proliferation rate, the model has a stable immune-free steady state, which predicts that patients have a high viral load. Above the higher immune proliferation rate, the model has a stable low infected steady state, which indicates that patients are under elite control. Between the two immune thresholds, the model exhibits the dynamic behavior of bistability, which suggests that patients either undergo viral rebound after treatment termination or achieve the post-treatment control. These results may explain the different virus dynamics in HIV-infected patients.
Collapse
Affiliation(s)
- SHAOLI WANG
- School of Mathematics and Statistics, Henan University, Kaifeng, Henan 475001, P. R. China
| | - FEI XU
- Department of Mathematics, Wilfrid Laurier University, Waterloo, Ontario N2L 3C5, Canada
| | - LIBIN RONG
- Department of Mathematics, University of Florida, Gainesville FL 32611, USA
| |
Collapse
|
6
|
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: 3.9] [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
|
7
|
Sun GQ, Jusup M, Jin Z, Wang Y, Wang Z. Pattern transitions in spatial epidemics: Mechanisms and emergent properties. Phys Life Rev 2016; 19:43-73. [PMID: 27567502 PMCID: PMC7105263 DOI: 10.1016/j.plrev.2016.08.002] [Citation(s) in RCA: 81] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2016] [Accepted: 08/04/2016] [Indexed: 12/19/2022]
Abstract
Infectious diseases are a threat to human health and a hindrance to societal development. Consequently, the spread of diseases in both time and space has been widely studied, revealing the different types of spatial patterns. Transitions between patterns are an emergent property in spatial epidemics that can serve as a potential trend indicator of disease spread. Despite the usefulness of such an indicator, attempts to systematize the topic of pattern transitions have been few and far between. We present a mini-review on pattern transitions in spatial epidemics, describing the types of transitions and their underlying mechanisms. We show that pattern transitions relate to the complexity of spatial epidemics by, for example, being accompanied with phenomena such as coherence resonance and cyclic evolution. The results presented herein provide valuable insights into disease prevention and control, and may even be applicable outside epidemiology, including other branches of medical science, ecology, quantitative finance, and elsewhere.
Collapse
Affiliation(s)
- Gui-Quan Sun
- Complex Systems Research Center, Shanxi University, Taiyuan, Shanxi 030006, PR China; School of Mathematical Sciences, Fudan University, Shanghai 200433, PR China.
| | - Marko Jusup
- Department of Vector Ecology and Environment, Nagasaki University Institute of Tropical Medicine (NEKKEN), Nagasaki 852-8523, Japan; Center of Mathematics for Social Creativity, Hokkaido University, Sapporo 060-0812, Japan.
| | - Zhen Jin
- Complex Systems Research Center, Shanxi University, Taiyuan, Shanxi 030006, PR China.
| | - Yi Wang
- Department of Mathematics, Southeast University, Nanjing 210096, PR China; Department of Mathematics and Statistics, University of Victoria, Victoria BC V8W 3R4, Canada
| | - Zhen Wang
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Fukuoka, 816-8580, Japan.
| |
Collapse
|
8
|
Bielejec F, Baele G, Rodrigo AG, Suchard MA, Lemey P. Identifying predictors of time-inhomogeneous viral evolutionary processes. Virus Evol 2016; 2:vew023. [PMID: 27774306 PMCID: PMC5072463 DOI: 10.1093/ve/vew023] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Various factors determine the rate at which mutations are generated and fixed in viral genomes. Viral evolutionary rates may vary over the course of a single persistent infection and can reflect changes in replication rates and selective dynamics. Dedicated statistical inference approaches are required to understand how the complex interplay of these processes shapes the genetic diversity and divergence in viral populations. Although evolutionary models accommodating a high degree of complexity can now be formalized, adequately informing these models by potentially sparse data, and assessing the association of the resulting estimates with external predictors, remains a major challenge. In this article, we present a novel Bayesian evolutionary inference method, which integrates multiple potential predictors and tests their association with variation in the absolute rates of synonymous and non-synonymous substitutions along the evolutionary history. We consider clinical and virological measures as predictors, but also changes in population size trajectories that are simultaneously inferred using coalescent modelling. We demonstrate the potential of our method in an application to within-host HIV-1 sequence data sampled throughout the infection of multiple patients. While analyses of individual patient populations lack statistical power, we detect significant evidence for an abrupt drop in non-synonymous rates in late stage infection and a more gradual increase in synonymous rates over the course of infection in a joint analysis across all patients. The former is predicted by the immune relaxation hypothesis while the latter may be in line with increasing replicative fitness during the asymptomatic stage.
Collapse
Affiliation(s)
- Filip Bielejec
- Department of Microbiology and Immunology, Rega Institute, KU Leuven, Leuven, Belgium
| | - Guy Baele
- Department of Microbiology and Immunology, Rega Institute, KU Leuven, Leuven, Belgium
| | - Allen G Rodrigo
- Research School of Biology, Australian National University, Canberra, ACT, Australia
| | - Marc A Suchard
- Department of Biomathematics and Human Genetics, David Geffen School of Medicine at UCLA, University of California, Los Angeles, CA 90095, USA; Department of Biostatistics, UCLA Fielding School of Public Health, University of California, Los Angeles, CA 90095, USA
| | - Philippe Lemey
- Department of Microbiology and Immunology, Rega Institute, KU Leuven, Leuven, Belgium
| |
Collapse
|
9
|
Abstract
Antiretroviral therapy (ART) for HIV is not a cure. However, recent studies suggest that ART, initiated early during primary infection, may induce post-treatment control (PTC) of HIV infection with HIV RNA maintained at <50 copies per mL. We investigate the hypothesis that ART initiated early during primary infection permits PTC by limiting the size of the latent reservoir, which, if small enough at treatment termination, may allow the adaptive immune response to prevent viral rebound (VR) and control infection. We use a mathematical model of within host HIV dynamics to capture interactions among target cells, productively infected cells, latently infected cells, virus, and cytotoxic T lymphocytes (CTLs). Analysis of our model reveals a range in CTL response strengths where a patient may show either VR or PTC, depending on the size of the latent reservoir at treatment termination. Below this range, patients will always rebound, whereas above this range, patients are predicted to behave like elite controllers. Using data on latent reservoir sizes in patients treated during primary infection, we also predict population-level VR times for noncontrollers consistent with observations.
Collapse
|
10
|
Fan R, Dong Y, Huang G, Takeuchi Y. Apoptosis in virus infection dynamics models. JOURNAL OF BIOLOGICAL DYNAMICS 2014; 8:20-41. [PMID: 24963975 PMCID: PMC4220821 DOI: 10.1080/17513758.2014.895433] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/29/2013] [Accepted: 02/04/2014] [Indexed: 06/03/2023]
Abstract
In this paper, on the basis of the simplified two-dimensional virus infection dynamics model, we propose two extended models that aim at incorporating the influence of activation-induced apoptosis which directly affects the population of uninfected cells. The theoretical analysis shows that increasing apoptosis plays a positive role in control of virus infection. However, after being included the third population of cytotoxic T lymphocytes immune response in HIV-infected patients, it shows that depending on intensity of the apoptosis of healthy cells, the apoptosis can either promote or comfort the long-term evolution of HIV infection. Further, the discrete-time delay of apoptosis is incorporated into the pervious model. Stability switching occurs as the time delay in apoptosis increases. Numerical simulations are performed to illustrate the theoretical results and display the different impacts of a delay in apoptosis.
Collapse
Affiliation(s)
- Ruili Fan
- School of Mathematics and Physics, China University of Geosciences, Wuhan430074, China
| | - Yueping Dong
- Graduate School of Science and Technology, Shizuoka University, Hamamatsu432-8561, Japan
| | - Gang Huang
- School of Mathematics and Physics, China University of Geosciences, Wuhan430074, China
| | - Yasuhiro Takeuchi
- Department of Physics and Mathematics, Aoyama Gakuin University, Sagamihara252-5258, Japan
| |
Collapse
|
11
|
Huang G, Takeuchi Y, Korobeinikov A. HIV evolution and progression of the infection to AIDS. J Theor Biol 2012; 307:149-59. [PMID: 22634206 DOI: 10.1016/j.jtbi.2012.05.013] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2011] [Revised: 05/11/2012] [Accepted: 05/14/2012] [Indexed: 12/12/2022]
Abstract
In this paper, we propose and discuss a possible mechanism, which, via continuous mutations and evolution, eventually enables HIV to break from immune control. In order to investigate this mechanism, we employ a simple mathematical model, which describes the relationship between evolving HIV and the specific CTL response and explicitly takes into consideration the role of CD4(+)T cells (helper T cells) in the activation of the CTL response. Based on the assumption that HIV evolves towards higher replication rates, we quantitatively analyze the dynamical properties of this model. The model exhibits the existence of two thresholds, defined as the immune activation threshold and the immunodeficiency threshold, which are critical for the activation and persistence of the specific cell-mediated immune response: the specific CTL response can be established and is able to effectively control an infection when the virus replication rate is between these two thresholds. If the replication rate is below the immune activation threshold, then the specific immune response cannot be reliably established due to the shortage of antigen-presenting cells. Besides, the specific immune response cannot be established when the virus replication rate is above the immunodeficiency threshold due to low levels of CD4(+)T cells. The latter case implies the collapse of the immune system and beginning of AIDS. The interval between these two thresholds roughly corresponds to the asymptomatic stage of HIV infection. The model shows that the duration of the asymptomatic stage and progression of the disease are very sensitive to variations in the model parameters. In particularly, the rate of production of the naive lymphocytes appears to be crucial.
Collapse
Affiliation(s)
- Gang Huang
- School of Mathematics and Physics, China University of Geosciences, Wuhan 430074, PR China
| | | | | |
Collapse
|