Regular articleTarget Cell Limited and Immune Control Models of HIV Infection: A Comparison
Abstract
We develop various mathematical models of the clinical latency stage of HIV-1 infection assuming that HIV-1 infection is limited either by the availability of cells that HIV can infect or by a specific anti-HIV cellular immune response. The former models we call “target-cell-limited”. Comparing the models by phase plane analysis we find that they all belong to the class of predator-prey models. In the target-cell limited models the virus is a predator feeding upon target cell prey, while in the immune-control models the virus is a prey that is controlled by an immune response predator. Because both classes of models are of predator-prey type they behave similarly in most circumstances. We find that both types of model can account for the generic picture of disease progression in which the CD4 T cell count slowly decreases and the viral load slowly increases. Additionally, we find that both types of models can adequately describe the clinically observed changes in the plasma HIV-1 RNA loads in response to retroviral therapies.
References (0)
Cited by (297)
Dynamic behaviors of a stochastic virus infection model with Beddington–DeAngelis incidence function, eclipse-stage and Ornstein–Uhlenbeck process
2024, Mathematical BiosciencesIn this paper, we present a virus infection model that incorporates eclipse-stage and Beddington–DeAngelis function, along with perturbation in infection rate using logarithmic Ornstein–Uhlenbeck process. Rigorous analysis demonstrates that the stochastic model has a unique global solution. Through construction of appropriate Lyapunov functions and a compact set, combined with the strong law of numbers and Fatou’s lemma, we obtain the existence of the stationary distribution under a critical condition, which indicates the long-term persistence of T-cells and virions. Moreover, a precise probability density function is derived around the quasi-equilibrium of the model, and spectral radius analysis is employed to identify critical condition for elimination of the virus. Finally, numerical simulations are presented to validate theoretical results, and the impact of some key parameters such as the speed of reversion, volatility intensity and mean infection rate are investigated.
What explains the poor contraction of the viral load during paediatric HIV infection?
2023, Journal of Theoretical BiologyAn acute HIV infection in young children differs markedly from that in adults: Children have higher viral loads (VL), and a poor contraction to a setpoint VL that is not much lower than the peak VL. As a result, children progress faster towards AIDS in the absence of treatment. We used a classical ordinary differential equation model for viral infection dynamics to study why children have a lower viral contraction ratio than adults. We performed parameter sweeps to identify factors explaining the observed difference between children and adults. We grouped parameters associated with the host, the infection, or the immune response. Based on paediatric data available from datasets within the EPIICAL project (https://www.epiical.org/), we refuted that viral replication rates differ between young children and adults, and therefore these cannot be responsible for the low VL contraction ratios seen in children. The major differences in lowering VL contraction ratio resulted from sweeping the parameters linked to the immune response. Thus, we postulate that an ”ineffective” (late and/or weak) immune response is the most parsimonious explanation for the higher setpoint VL in young children, and hence the reason for their fast disease progression.
Modeling the contribution of antibody attack rates to single and dual helminth infections in a natural system
2023, Mathematical BiosciencesWithin-host models of infection can provide important insights into the processes that affect parasite spread and persistence in host populations. However, modeling can be limited by the availability of empirical data, a problem commonly encountered in natural systems. Here, we used six years of immune-infection observations of two gastrointestinal helminths (Trichostrongylus retortaeformis and Graphidium strigosum) from a population of European rabbits (Oryctolagus cuniculus) to develop an age-dependent, mathematical model that explicitly included species-specific and cross-reacting antibody (IgA and IgG) responses to each helminth in hosts with single or dual infections. Different models of single infection were formally compared to test alternative mechanisms of parasite regulation. The two models that best described single infections of each helminth species were then coupled through antibody cross-immunity to examine how the presence of one species could alter the host immune response to, and the within-host dynamics of, the other species. For both single infections, model selection suggested that either IgA or IgG responses could equally explain the observed parasite intensities by host age. However, the antibody attack rate and affinity level changed between the two helminths, it was stronger against T. retortaeformis than against G. strigosum and caused contrasting age-intensity profiles. When the two helminths coinfect the same host, we found variation of the species-specific antibody response to both species together with an asymmetric cross-immune response driven by IgG. Lower attack rate and affinity of antibodies in dual than single infections contributed to the significant increase of both helminth intensities. By combining mathematical modeling with immuno-infection data, our work provides a tractable model framework for disentangling some of the complexities generated by host-parasite and parasite–parasite interactions in natural systems.
Excitable dynamics in a molecularly-explicit model of cell motility: Mixed-mode oscillations and beyond: Mixed-mode oscillations in cell motility
2023, Journal of Theoretical BiologyMesenchymal cell motility is mainly regulated by two members of the Rho-family of GTPases, called Rac and Rho. The mutual inhibition exerted by these two proteins on each other’s activation and the promotion of Rac activation by an adaptor protein called paxillin have been implicated in driving cellular polarization comprised of front (high active Rac) and back (high active Rho) during cell migration. Mathematical modeling of this regulatory network has previously shown that bistability is responsible for generating a spatiotemporal pattern underscoring cellular polarity called wave-pinning when diffusion is included. We previously developed a 6V reaction–diffusion model of this network to decipher the role of Rac, Rho and paxillin (along with other auxiliary proteins) in generating wave-pinning. In this study, we simplify this model through a series of steps into an excitable 3V ODE model comprised of one fast variable (the scaled concentration of active Rac), one slow variable (the maximum paxillin phosphorylation rate — turned into a variable) and a very slow variable (a recovery rate — also turned into a variable). We then explore, through slow–fast analysis, how excitability is manifested by showing that the model can exhibit relaxation oscillations (ROs) as well as mixed-mode oscillations (MMOs) whose underlying dynamics are consistent with a delayed Hopf bifurcation with a canard explosion. By reintroducing diffusion and the scaled concentration of inactive Rac into the model, we obtain a 4V PDE model that generates several unique spatiotemporal patterns that are relevant to cell motility. These patterns are then characterized and their impact on cell motility are explored by employing the cellular potts model (CPM). Our results reveal that wave pinning produces purely very directed motion in CPM, while MMOs allow for meandering and non-motile behaviors to occur. This highlights the role of MMOs as a potential mechanism for mesenchymal cell motility.
Modeling the effects of drugs of abuse on within-host dynamics of two HIV species
2023, Journal of Theoretical BiologyInjection drug use is one of the most significant risk factors associated with contracting human immunodeficiency virus (HIV), and drug users infected with HIV suffer from a higher viral load and rapid disease progression. While replication of HIV may result in many mutant viruses that can escape recognition of the host’s immune response, the presence of morphine (a drug of abuse) can decrease the viral mutation rate and cellular immune responses. This study develops a mathematical model to explore the effects of morphine-altered mutation and cellular immune response on the within-host dynamics of two HIV species, a wild-type and a mutant. Our model predicts that the morphine-altered mutation rate and cellular immune response allow the wild-type virus to outcompete the mutant virus, resulting in a higher set point viral load and lower CD4 count. We also compute the basic reproduction numbers and show that the dominant species is determined by morphine concentration, with the mutant dominating below and the wild-type dominating above a threshold. Furthermore, we identified three biologically relevant equilibria, infection-free, mutant-only, and coexistence, which are completely characterized by the fitness cost of mutation, mutant escape rate, and morphine concentration.
Modeling the virus-induced tumor-specific immune response with delay in tumor virotherapy
2022, Communications in Nonlinear Science and Numerical SimulationIt is urgently required to design novel cancer therapies which lead to permanent cancer eradication or cancer control. Oncolytic virotherapy is a promising cancer treatment strategy using genetically engineered viruses which can selectively infect, replicate in and kill tumor cells without harming normal cells. Due to the variable interactions between tumor cells and oncolytic viruses (OVs) as well as the accordingly immune response, the impact of tumor virotherapy on tumor control is prominent. We propose a novel mathematical modeling framework based on delay differential equation to study tumor virotherapy with mediated antitumor immunity by OVs, which incorporates complex tumor–virus–immune system interactions. We initially study the existence and local stability of equilibria, and theoretically and numerically investigate the local Hopf bifurcation from the positive equilibrium by considering the time delay as a bifurcation parameter. Further, we analyze the direction of Hopf bifurcation and the stability of bifurcating periodic solution. Main results show that the time delay can induce Hopf bifurcation and result in periodic oscillations, indicating that the delayed tumor-specific CTL response induced by OVs leads to complex dynamics and may significantly influence the development process of tumor growth. Our findings also provide insight into important aspects of the virotherapy, including the dependence of the efficacy on key factors. We find that enhancing the induction of antitumor immunity by OVs can reduce the complexity of dynamics by strengthening the stability and show considerable effects on combating tumor. Further, proper choice of OVs with relative strong ability to lyse tumor cells is beneficial to ultimately control the development process of tumor growth.
- f1
Author to whom correspondence should be addressed. E-mail: [email protected]