Within-host Model Research Articles

Modelling and investigating memory immune responses in infectious disease. Application to influenza a virus and sars-cov-2 reinfections

Understanding effector and memory immune responses against influenza A virus (IAV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections and re-infections is extremely important, given that they are now endemic in the community. The goal of this study is to investigate the role of memory cells and antibodies in the immune responses against IAV and SARS-CoV-2 re-infections. To this end, we adapt a previously-published within-host mathematical model (Sadria & Layton, 2021) for the primary immune response against SARS-CoV-2 infections, by including two types of memory immune cells, i.e., memory CD8+ T-cells and memory B-cells, and by parametrising the new model with values specific to the two viruses. We first investigate the long-term dynamics of the model by identifying the virus-free steady states and studying the conditions that ensure the stability of these states. Then, we investigate the transient dynamics of this in-host model by simulating different viral reinfection times: 20 days, 60 days and 400 days after the first encounter with the pathogen. This allows us to highlight which memory immune components have the greatest impact on the viral elimination depending on the time of reinfection. Our results suggest that memory immune responses have a greater impact in the case of IAV infections compared to SARS-CoV-2 infections. Moreover, we observe that the immune response after a secondary infection is more efficient when the reinfection occurs at a shorter time.

Modelling the mechanisms of antibody mixtures in viral infections: the cases of sequential homologous and heterologous dengue infections.

Antibodies play an essential role in the immune response to viral infections, vaccination or antibody therapy. Nevertheless, they can be either protective or harmful during the immune response. Moreover, competition or cooperation between mixed antibodies can enhance or reduce this protective or harmful effect. Using the laws of chemical reactions, we propose a new approach to modelling the antigen-antibody complex activity. The resulting expression covers not only purely competitive or purely independent binding but also synergistic binding which, depending on the antibodies, can promote either neutralization or enhancement of viral activity. We then integrate this expression of viral activity in a within-host model and investigate the existence of steady-states and their asymptotic stability. We complete our study with numerical simulations to illustrate different scenarios: firstly, where both antibodies are neutralizing and secondly, where one antibody is neutralizing and the other enhancing. The results indicate that efficient viral neutralization is associated with purely independent antibody binding, whereas strong viral activity enhancement is expected in the case of purely competitive antibody binding. Finally, data collected during a secondary dengue infection were used to validate the model. The dataset includes sequential measurements of virus and antibody titres during viremia in patients. Data fitting shows that the two antibodies are in strong competition, as the synergistic binding is low. This contributes to the high levels of virus titres and may explain the antibody-dependent enhancement phenomenon. Besides, the mortality of infected cells is almost twice as high as that of susceptible cells, and the heterogeneity of viral kinetics in patients is associated with variability in antibody responses between individuals. Other applications of the model may be considered, such as the efficacy of vaccines and antibody-based therapies.

Modeling virus-stimulated proliferation of CD4[formula omitted] T-cell, cell-to-cell transmission and viral loss in HIV infection dynamics

Human immunodeficiency virus (HIV) can persist in infected individuals despite prolonged antiretroviral therapy and it may spread through two modes: virus-to-cell and cell-to-cell transmissions. Understanding viral infection dynamics is pivotal for elucidating HIV pathogenesis. In this study, we incorporate the loss term of virions, and both virus-to-cell and cell-to-cell infection modes into a within-host HIV model, which also takes into consideration the proliferation of healthy target cells stimulated by free viruses. By constructing suitable Lyapunov function and applying geometric methods, we establish global stability results of the infection free equilibrium and the infection persistent equilibrium, respectively. Our findings highlight the crucial role of the basic reproduction number in the threshold dynamics. Moreover, we use the loss rate of virions as the bifurcation parameter to investigate stability switches of the positive equilibrium, local Hopf bifurcation, and its global continuation. Numerical simulations validate our theoretical results, revealing rich viral dynamics including backward bifurcation, saddle–node bifurcation, and bistability phenomenon in the sense that the infection free equilibrium and a limit cycle are both locally asymptotically stable. These insights contribute to a deeper understanding of HIV dynamics and inform the development of effective therapeutic strategies.

Identifiability and Parameter Estimation of Within-Host Model of HIV with Immune Response

Identifiability and Parameter Estimation of Within-Host Model of HIV with Immune Response

Quantifying trade-offs between therapeutic efficacy and resistance dissemination for enrofloxacin dose regimens in cattle

The use of antimicrobial drugs in food-producing animals contributes to the selection pressure on pathogenic and commensal bacteria to become resistant. This study aims to evaluate the existence of trade-offs between treatment effectiveness, cost, and the dynamics of resistance in gut commensal bacteria. We developed a within-host ordinary differential equation model to track the dynamics of antimicrobial drug concentrations and bacterial populations in the site of infection (lung) and the gut. The model was parameterized to represent enrofloxacin treatment for bovine respiratory disease (BRD) caused by Pastereulla multocida in cattle. Three approved enrofloxacin dosing regimens were compared for their effects on resistance on P. multocida and commensal E. coli: 12.5 mg/kg and 7.5 mg/kg as a single dose, and 5 mg/kg as three doses. Additionally, we explored non-FDA-approved regimes. Our results indicated that both 12.5 mg/kg and 7.5 mg/kg as a single dose scenario increased the most the treatment costs and prevalence of P. multocida resistance in the lungs, while 5 mg/kg as three doses increased resistance in commensal E. coli bacteria in the gut the most out of the approved scenarios. A proposed non-FDA-approved scenario (7.5 mg/kg, two doses 24 h apart) showed low economic costs, minimal P. multocida, and moderate effects on resistant E. coli. Overall, the scenarios that decrease P. multocida, including resistant P. multocida did not coincide with those that decrease resistant E. coli the most, suggesting a trade-off between both outcomes. The sensitivity analysis suggests that bacterial populations were the most sensitive to drug conversion factors into plasma (β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\beta}$$\\end{document}), elimination of the drug from the colon (ϑ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\vartheta$$\\end{document}), fifty percent sensitive bacteria (P. multocida) killing effect (Ls50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{L}}_{\ ext{s50}}$$\\end{document}), fifty percent of bacteria (E. coli) above ECOFF killing effect (Cr50\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ ext{C}}_{\ ext{r50}}$$\\end{document}), and net drug transfer rate in the lung (γ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\gamma$$\\end{document}) parameters.

Recent approaches in computational modelling for controlling pathogen threats.

In this review, we assess the status of computational modelling of pathogens. We focus on three disparate but interlinked research areas that produce models with very different spatial and temporal scope. First, we examine antimicrobial resistance (AMR). Many mechanisms of AMR are not well understood. As a result, it is hard to measure the current incidence of AMR, predict the future incidence, and design strategies to preserve existing antibiotic effectiveness. Next, we look at how to choose the finite number of bacterial strains that can be included in a vaccine. To do this, we need to understand what happens to vaccine and non-vaccine strains after vaccination programmes. Finally, we look at within-host modelling of antibody dynamics. The SARS-CoV-2 pandemic produced huge amounts of antibody data, prompting improvements in this area of modelling. We finish by discussing the challenges that persist in understanding these complex biological systems.

Modeling and analysis of a two-strain immuno-epidemiological model with reinfection

In this paper, we formulate a two-strain model with reinfection that combines immunological and epidemiological dynamics across scales, using the COVID-19 pandemic as a case study. Firstly, we conduct a qualitative analysis of both within-host and between-host models. For the within-host model, we prove the existence and stability of equilibria, and Hopf bifurcation occurs from the infection equilibrium with immune response. This implies that, under specific immune states, the virus within an infected individual may persist, and its concentration may also oscillate periodically. For the between-host model, the disease-free equilibrium always exists and is locally asymptotically stable when the epidemiological basic reproduction number ℜ0<1. In addition, the model can have boundary equilibria of strain 1 or strain 2, which are locally asymptotically stable under specific conditions. However, the co-existence equilibrium does not exist. Secondly, to explore the infection and transmission mechanisms of two strain models and obtain reliable parameter values, we utilize statistical data to fit the immuno-epidemiological model. Simultaneously, we conduct an identifiability analysis of the immuno-epidemiological model to ensure the robustness of the fitted parameters. The results demonstrate the reliable estimation of parameter ranges for structurally unidentifiable parameters with minor measurement errors using the affine invariant ensemble Markov Chain Monte Carlo algorithm (GWMCMC). Moreover, simulations illustrate that enhancing treatment of patients infected with BA.2 strains to inhibit the number of viruses released by infected cells can significantly reduce disease spread.

Validation of a Multi-Strain HIV Within-Host Model with AIDS Clinical Studies

Validation of a Multi-Strain HIV Within-Host Model with AIDS Clinical Studies

Differences in virus and immune dynamics for SARS-CoV-2 Delta and Omicron infections by age and vaccination histories

Vaccination against COVID-19 was integral to controlling the pandemic that persisted with the continuous emergence of SARS-CoV-2 variants. Using a mathematical model describing SARS-CoV-2 within-host infection dynamics, we estimate differences in virus and immunity due to factors of infecting variant, age, and vaccination history (vaccination brand, number of doses and time since vaccination). We fit our model in a Bayesian framework to upper respiratory tract viral load measurements obtained from cases of Delta and Omicron infections in Singapore, of whom the majority only had one nasopharyngeal swab measurement. With this dataset, we are able to recreate similar trends in URT virus dynamics observed in past within-host modelling studies fitted to longitudinal patient data.We found that Omicron had higher R0,within values than Delta, indicating greater initial cell-to-cell spread of infection within the host. Moreover, heterogeneities in infection dynamics across patient subgroups could be recreated by fitting immunity-related parameters as vaccination history-specific, with or without age modification. Our model results are consistent with the notion of immunosenescence in SARS-CoV-2 infection in elderly individuals, and the issue of waning immunity with increased time since last vaccination. Lastly, vaccination was not found to subdue virus dynamics in Omicron infections as well as it had for Delta infections.This study provides insight into the influence of vaccine-elicited immunity on SARS-CoV-2 within-host dynamics, and the interplay between age and vaccination history. Furthermore, it demonstrates the need to disentangle host factors and changes in pathogen to discern factors influencing virus dynamics. Finally, this work demonstrates a way forward in the study of within-host virus dynamics, by use of viral load datasets including a large number of patients without repeated measurements.

Digest: Resource limitation as a mechanism for constraining the evolution of virulence in malaria parasites.

The virulence of parasites is expected to reflect an evolutionary tradeoff between increasing proliferation rates that enhance transmission and host mortality which curtails transmission. However, host resource availability may also limit parasites' proliferation rate. To understand the role of resource limitation as a driver of virulence evolution, Pak et al. (2024) use a within-host model of red blood cell (RBC) invasion by Plasmodium chabaudi. They find that within-host resource consumption limits the evolution of the parasite's proliferation rate, as the depletion of RBCs during infection results in intermediate optimal virulence. These results suggest that resource limitation, rather than host mortality, may drive the evolution of virulence.

Identifiability investigation of within-host models of acute virus infection.

Uncertainty in parameter estimates from fitting within-host models to empirical data limits the model's ability to uncover mechanisms of infection, disease progression, and to guide pharmaceutical interventions. Understanding the effect of model structure and data availability on model predictions is important for informing model development and experimental design. To address sources of uncertainty in parameter estimation, we use four mathematical models of influenza A infection with increased degrees of biological realism. We test the ability of each model to reveal its parameters in the presence of unlimited data by performing structural identifiability analyses. We then refine the results by predicting practical identifiability of parameters under daily influenza A virus titers alone or together with daily adaptive immune cell data. Using these approaches, we present insight into the sources of uncertainty in parameter estimation and provide guidelines for the types of model assumptions, optimal experimental design, and biological information needed for improved predictions.

Mathematical analysis of a within-host dengue virus dynamics model with adaptive immunity using Caputo fractional-order derivatives

AbstractDengue fever poses a significant global health threat, with over 50 million annual infections spanning more than 100 countries. Given the absence of a specific treatment, medical intervention primarily targets symptom alleviation. The present study utilizes a Caputo-type fractional-order derivative operator to investigate and analyze the dynamics of dengue virus spread within a host with adaptive immune responses. The developed model describes and analyzes the dynamics of immune cells, free dengue particles, infected monocytes, and susceptible monocytes in the presence of cytotoxic T-Lymphocytes. A range of analytical methods is employed to probe the fractional-order within-host model. The application of the generalized mean value theorem aids in investigating the model’s solutions, employing positivity and boundedness theory. Furthermore, the Banach fixed-point approach is utilized to establish the existence and uniqueness of solutions. Employing the normalized forward sensitivity approach, the fractional-order system’s response to various model parameters is scrutinized. The study reveals that the dynamics of the viral model are significantly influenced by the transmission rate and parameters representing adaptive immune responses. Numerical simulations underscore the critical role of transmission rates and adaptive immune responses in the model. Additionally, the study examines the impact of memory on the density of susceptible monocytes, infected monocytes, free dengue particles, and immune cells to optimize immune responses. Through simulations, the study illustrates the influence of memory on immune dynamics.

Proliferation in malaria parasites: How resource limitation can prevent evolution of greater virulence.

For parasites, robust proliferation within hosts is crucial for establishing the infection and creating opportunities for onward transmission. While faster proliferation enhances transmission rates, it is often assumed to curtail transmission duration by killing the host (virulence), a trade-off constraining parasite evolution. Yet in many diseases, including malaria, the preponderance of infections with mild or absent symptoms suggests that host mortality is not a sufficient constraint, raising the question of what restrains evolution toward faster proliferation. In malaria infections, the maximum rate of proliferation is determined by the burst size, the number of daughter parasites produced per infected red blood cell. Larger burst sizes should expand the pool of infected red blood cells that can be used to produce the specialized transmission forms needed to infect mosquitoes. We use a within-host model parameterized for rodent malaria parasites (Plasmodium chabaudi) to project the transmission consequences of burst size, focusing on initial acute infection where resource limitation and risk of host mortality are greatest. We find that resource limitation restricts evolution toward higher burst sizes below the level predicted by host mortality alone. Our results suggest resource limitation could represent a more general constraint than virulence-transmission trade-offs, preventing evolution towards faster proliferation.

Modeling and analysis of the effect of optimal virus control on the spread of HFMD

A within-host and between-host hand, foot and mouth disease (HFMD) mathematical model is established and the affect of optimal control in its within-host part on HFMD transmission is studied. Through define two basic reproduction numbers, by using the fast-slow system analysis method of time scale, the global stabilities of the between-host (slow) system and within-host (fast) system are researched, respectively. An optimal control problem with drug-treatment control on coupled within-host and between-host HFMD model is formulated and analysed theoretically. Finally, the purposed optimal control measures are applied to the actual HFMD epidemic analysis in Zhejiang Province, China from April 1, 2021 to June 30, 2021. The numerical results show that the drug control strategies can reduce the virus load per capita and can effectively prevent large-scale outbreaks of HFMD.

The role of long-lived plasma cells in viral clearance

ABSTRACT The adaptive immune system has two types of plasma cells (PC), long-lived plasma cells (LLPC) and short-lived plasma cells (SLPC), that differ in their lifespan. In this paper, we propose that LLPC is crucial to the clearance of viral particles in addition to reducing the viral basic reproduction number in secondary infections. We use a sequence of within-host mathematical models to show that, CD8 T cells, SLPC and memory B cells cannot achieve full viral clearance, and the viral load will reach a low positive equilibrium level because of a continuous replenishment of target cells. However, the presence of LLPC is crucial for viral clearance.

Analysis of a Delayed Multiscale AIDS/HIV-1 Model Coupling Between-Host and Within-Host Dynamics

Taking into account the effects of the immune response and delay, and complexity on HIV-1 transmission, a multiscale AIDS/HIV-1 model is formulated in this paper. The multiscale model is described by a within-host fast time model with intracellular delay and immune delay, and a between-host slow time model with latency delay. The dynamics of the fast time model is analyzed, and includes the stability of equilibria and properties of Hopf bifurcation. Further, for the coupled slow time model without an immune response, the basic reproduction number R0h is defined, which determines whether the model may have zero, one, or two positive equilibria under different conditions. This implies that the slow time model demonstrates more complex dynamic behaviors, including saddle-node bifurcation, backward bifurcation, and Hopf bifurcation. For the other case, that is, the coupled slow time model with an immune response, the threshold dynamics, based on the basic reproduction number R˜0h, is rigorously investigated. More specifically, if R˜0h&lt;1, the disease-free equilibrium is globally asymptotically stable; if R˜0h&gt;1, the model exhibits a unique endemic equilibrium that is globally asymptotically stable. With regard to the coupled slow time model with an immune response and stable periodic solution, the basic reproduction number R0 is derived, which serves as a threshold value determining whether the disease will die out or lead to periodic oscillations in its prevalence. The research results suggest that the disease is more easily controlled when hosts have an extensive immune response and the time required for new immune particles to emerge in response to antigenic stimulation is within a certain range. Finally, numerical simulations are presented to validate the main results and provide some recommendations for controlling the spread of HIV-1.

Within-host bayesian joint modeling of longitudinal and time-to-event data of Leishmania infection.

The host immune system plays a significant role in managing and clearing pathogen material during an infection, but this complex process presents numerous challenges from a modeling perspective. There are many mathematical and statistical models for these kinds of processes that take into account a wide range of events that happen within the host. In this work, we present a Bayesian joint model of longitudinal and time-to-event data of Leishmania infection that considers the interplay between key drivers of the disease process: pathogen load, antibody level, and disease. The longitudinal model also considers approximate inflammatory and regulatory immune factors. In addition to measuring antibody levels produced by the immune system, we adapt data from CD4+ and CD8+ T cell proliferation, and expression of interleukin 10, interferon-gamma, and programmed cell death 1 as inflammatory or regulatory factors mediating the disease process. The model is developed using data collected from a cohort of dogs naturally exposed to Leishmania infantum. The cohort was chosen to start with healthy infected animals, and this is the majority of the data. The model also characterizes the relationship features of the longitudinal outcomes and time-to-death due to progressive Leishmania infection. In addition to describing the mechanisms causing disease progression and impacting the risk of death, we also present the model's ability to predict individual trajectories of Canine Leishmaniosis (CanL) progression. The within-host model structure we present here provides a way forward to address vital research questions regarding the understanding of the progression of complex chronic diseases such as Visceral Leishmaniasis, a parasitic disease causing significant morbidity worldwide.

Within-host models unravelling the dynamics of dengue reinfections

Caused by four serotypes, dengue fever is a major public health concern worldwide. Current modeling efforts have mostly focused on primary and heterologous secondary infections, assuming that lifelong immunity prevents reinfections by the same serotype. However, recent findings challenge this assumption, prompting a reevaluation of dengue immunity dynamics.In this study, we develop a within-host modeling framework to explore different scenarios of dengue infections. Unlike previous studies, we go beyond a deterministic framework, considering individual immunological variability. Both deterministic and stochastic models are calibrated using empirical data on viral load and antibody (IgM and IgG) concentrations for all dengue serotypes, incorporating confidence intervals derived from stochastic realizations.With good agreement between the mean of the stochastic realizations and the mean field solution for each model, our approach not only successfully captures primary and heterologous secondary infection dynamics facilitated by antibody-dependent enhancement (ADE) but also provides, for the first time, insights into homotypic reinfection dynamics. Our study discusses the relevance of homotypic reinfections in dengue transmission at the population level, highlighting potential implications for disease prevention and control strategies.

The effect of model structure and data availability on Usutu virus dynamics at three biological scales.

Understanding the epidemiology of emerging pathogens, such as Usutu virus (USUV) infections, requires systems investigation at each scale involved in the host-virus transmission cycle, from individual bird infections, to bird-to-vector transmissions, and to USUV incidence in bird and vector populations. For new pathogens field data are sparse, and predictions can be aided by the use of laboratory-type inoculation and transmission experiments combined with dynamical mathematical modelling. In this study, we investigated the dynamics of two strains of USUV by constructing mathematical models for the within-host scale, bird-to-vector transmission scale and vector-borne epidemiological scale. We used individual within-host infectious virus data and per cent mosquito infection data to predict USUV incidence in birds and mosquitoes. We addressed the dependence of predictions on model structure, data uncertainty and experimental design. We found that uncertainty in predictions at one scale change predicted results at another scale. We proposed in silico experiments that showed that sampling every 12 hours ensures practical identifiability of the within-host scale model. At the same time, we showed that practical identifiability of the transmission scale functions can only be improved under unrealistically high sampling regimes. Instead, we proposed optimal experimental designs and suggested the types of experiments that can ensure identifiability at the transmission scale and, hence, induce robustness in predictions at the epidemiological scale.

Quantifying the relationship between within-host dynamics and transmission for viral diseases of livestock.

Understanding the population dynamics of an infectious disease requires linking within-host dynamics and between-host transmission in a quantitative manner, but this is seldom done in practice. Here a simple phenomenological model for viral dynamics within a host is linked to between-host transmission by assuming that the probability of transmission is related to log viral titre. Data from transmission experiments for two viral diseases of livestock, foot-and-mouth disease virus in cattle and swine influenza virus in pigs, are used to parametrize the model and, importantly, test the underlying assumptions. The model allows the relationship between within-host parameters and transmission to be determined explicitly through their influence on the reproduction number and generation time. Furthermore, these critical within-host parameters (time and level of peak titre, viral growth and clearance rates) can be computed from more complex within-host models, raising the possibility of assessing the impact of within-host processes on between-host transmission in a more detailed quantitative manner.