A novel within-host model of HIV and nutrition.

In this paper we develop a four compartment within-host model of nutrition and HIV. We show that the model has two equilibria: an infection-free equilibrium and infection equilibrium. The infection free equilibrium is locally asymptotically stable when the basic reproduction number $ \mathcal{R}_0 < 1 $, and unstable when $ \mathcal{R}_0 > 1 $. The infection equilibrium is locally asymptotically stable if $ \mathcal{R}_0 > 1 $ and an additional condition holds. We show that the within-host model of HIV and nutrition is structured to reveal its parameters from the observations of viral load, CD4 cell count and total protein data. We then estimate the model parameters for these 3 data sets. We have also studied the practical identifiability of the model parameters by performing Monte Carlo simulations, and found that the rate of clearance of the virus by immunoglobulins is practically unidentifiable, and that the rest of the model parameters are only weakly identifiable given the experimental data. Furthermore, we have studied how the data frequency impacts the practical identifiability of model parameters.

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.

Optimal control of a multi-scale HIV-opioid model.

In this study, we apply optimal control theory to an immuno-epidemiological model of HIV and opioid epidemics. For the multi-scale model, we used four controls: treating the opioid use, reducing HIV risk behaviour among opioid users, entry inhibiting antiviral therapy, and antiviral therapy which blocks the viral production. Two population-level controls are combined with two within-host-level controls. We prove the existence and uniqueness of an optimal control quadruple. Comparing the two population-level controls, we find that reducing the HIV risk of opioid users has a stronger impact on the population who is both HIV-infected and opioid-dependent than treating the opioid disorder. The within-host-level antiviral treatment has an effect not only on the co-affected population but also on the HIV-only infected population. Our findings suggest that the most effective strategy for managing the HIV and opioid epidemics is combining all controls at both within-host and between-host scales.

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.

Modeling the interplay between albumin-globulin metabolism and HIV infection.

Human immunodeficiency virus (HIV) infection is a major public health concern with 1.2 million people living with HIV in the United States. The role of nutrition in general, and albumin/globulin in particular in HIV progression has long been recognized. However, no mathematical models exist to describe the interplay between HIV and albumin/globulin. In this paper, we present a family of models of HIV and the two protein components albumin and globulin. We use albumin, globulin, viral load and target cell data from simian immunodeficiency virus (SIV)-infected monkeys to perform model selection on the family of models. We discover that the simplest model accurately and uniquely describes the data. The selection of the simplest model leads to the observation that albumin and globulin do not impact the infection rate of target cells by the virus and the clearance of the infected target cells by the immune system. Moreover, the recruitment of target cells and immune cells are modeled independently of globulin in the selected model. Mathematical analysis of the selected model reveals that the model has an infection-free equilibrium and a unique infected equilibrium when the immunological reproduction number is above one. The infection-free equilibrium is locally stable when the immunological reproduction number is below one, and unstable when the immunological reproduction number is greater than one. The infection equilibrium is locally stable whenever it exists. To determine the parameters of the best fitted model we perform structural and practical identifiability analysis. The structural identifiability analysis reveals that the model is identifiable when the immune cell infection rate is fixed at a value obtained from the literature. Practical identifiability reveals that only seven of the sixteen parameters are practically identifiable with the given data. Practical identifiability of parameters performed with synthetic data sampled a lot more frequently reveals that only two parameters are practically unidentifiable. We conclude that experiments that will improve the quality of the data can help improve the parameter estimates and lead to better understanding of the interplay of HIV and albumin-globulin metabolism.

Structural identifiability analysis of epidemic models based on differential equations: a tutorial-based primer.

The successful application of epidemic models hinges on our ability to estimate model parameters from limited observations reliably. An often-overlooked step before estimating model parameters consists of ensuring that the model parameters are structurally identifiable from the observed states of the system. In this tutorial-based primer, intended for a diverse audience, including students training in dynamic systems, we review and provide detailed guidance for conducting structural identifiability analysis of differential equation epidemic models based on a differential algebra approach using differential algebra for identifiability of systems (DAISY) and Mathematica (Wolfram Research). This approach aims to uncover any existing parameter correlations that preclude their estimation from the observed variables. We demonstrate this approach through examples, including tutorial videos of compartmental epidemic models previously employed to study transmission dynamics and control. We show that the lack of structural identifiability may be remedied by incorporating additional observations from different model states, assuming that the system's initial conditions are known, using prior information to fix some parameters involved in parameter correlations, or modifying the model based on existing parameter correlations. We also underscore how the results of structural identifiability analysis can help enrich compartmental diagrams of differential-equation models by indicating the observed state variables and the results of the structural identifiability analysis.

A network immuno-epidemiological model of HIV and opioid epidemics.

In this paper, we introduce a novel multi-scale network model of two epidemics: HIV infection and opioid addiction. The HIV infection dynamics is modeled on a complex network. We determine the basic reproduction number of HIV infection, $ \mathcal{R}_{v} $, and the basic reproduction number of opioid addiction, $ \mathcal{R}_{u} $. We show that the model has a unique disease-free equilibrium which is locally asymptotically stable when both $ \mathcal{R}_{u} $ and $ \mathcal{R}_{v} $ are less than one. If $ \mathcal{R}_{u} > 1 $ or $ \mathcal{R}_{v} > 1 $, then the disease-free equilibrium is unstable and there exists a unique semi-trivial equilibrium corresponding to each disease. The unique opioid only equilibrium exist when the basic reproduction number of opioid addiction is greater than one and it is locally asymptotically stable when the invasion number of HIV infection, $ \mathcal{R}^{1}_{v_i} $ is less than one. Similarly, the unique HIV only equilibrium exist when the basic reproduction number of HIV is greater than one and it is locally asymptotically stable when the invasion number of opioid addiction, $ \mathcal{R}^{2}_{u_i} $ is less than one. Existence and stability of co-existence equilibria remains an open problem. We performed numerical simulations to better understand the impact of three epidemiologically important parameters that are at the intersection of two epidemics: $ q_v $ the likelihood of an opioid user being infected with HIV, $ q_u $ the likelihood of an HIV-infected individual becoming addicted to opioids, and $ \delta $ recovery from opioid addiction. Simulations suggest that as the recovery from opioid use increases, the prevalence of co-affected individuals, those who are addicted to opioids and are infected with HIV, increase significantly. We demonstrate that the dependence of the co-affected population on $ q_u $ and $ q_v $ are not monotone.

Identifiability of parameters in mathematical models of SARS-CoV-2 infections in humans.

Determining accurate estimates for the characteristics of the severe acute respiratory syndrome coronavirus 2 in the upper and lower respiratory tracts, by fitting mathematical models to data, is made difficult by the lack of measurements early in the infection. To determine the sensitivity of the parameter estimates to the noise in the data, we developed a novel two-patch within-host mathematical model that considered the infection of both respiratory tracts and assumed that the viral load in the lower respiratory tract decays in a density dependent manner and investigated its ability to match population level data. We proposed several approaches that can improve practical identifiability of parameters, including an optimal experimental approach, and found that availability of viral data early in the infection is of essence for improving the accuracy of the estimates. Our findings can be useful for designing interventions.

Parameter identifiability and optimal control of an SARS-CoV-2 model early in the pandemic.

We fit an SARS-CoV-2 model to US data of COVID-19 cases and deaths. We conclude that the model is not structurally identifiable. We make the model identifiable by prefixing some of the parameters from external information. Practical identifiability of the model through Monte Carlo simulations reveals that two of the parameters may not be practically identifiable. With thus identified parameters, we set up an optimal control problem with social distancing and isolation as control variables. We investigate two scenarios: the controls are applied for the entire duration and the controls are applied only for the period of time. Our results show that if the controls are applied early in the epidemic, the reduction in the infected classes is at least an order of magnitude higher compared to when controls are applied with 2-week delay. Further, removing the controls before the pandemic ends leads to rebound of the infected classes.

Immuno-epidemiological co-affection model of HIV infection and opioid addiction.

In this paper, we present a multi-scale co-affection model of HIV infection and opioid addiction. The population scale epidemiological model is linked to the within-host model which describes the HIV and opioid dynamics in a co-affected individual. CD4 cells and viral load data obtained from morphine addicted SIV-infected monkeys are used to validate the within-host model. AIDS diagnoses, HIV death and opioid mortality data are used to fit the between-host model. When the rates of viral clearance and morphine uptake are fixed, the within-host model is structurally identifiable. If in addition the morphine saturation and clearance rates are also fixed the model becomes practical identifiable. Analytical results of the multi-scale model suggest that in addition to the disease-addiction-free equilibrium, there is a unique HIV-only and opioid-only equilibrium. Each of the boundary equilibria is stable if the invasion number of the other epidemic is below one. Elasticity analysis suggests that the most sensitive number is the invasion number of opioid epidemic with respect to the parameter of enhancement of HIV infection of opioid-affected individual. We conclude that the most effective control strategy is to prevent opioid addicted individuals from getting HIV, and to treat the opioid addiction directly and independently from HIV.

Sensitivity Analysis in an Immuno-Epidemiological Vector-Host Model.

Sensitivity Analysis (SA) is a useful tool to measure the impact of changes in model parameters on the infection dynamics, particularly to quantify the expected efficacy of disease control strategies. SA has only been applied to epidemic models at the population level, ignoring the effect of within-host virus-with-immune-system interactions on the disease spread. Connecting the scales from individual to population can help inform drug and vaccine development. Thus the value of understanding the impact of immunological parameters on epidemiological quantities. Here we consider an age-since-infection structured vector-host model, in which epidemiological parameters are formulated as functions of within-host virus and antibody densities, governed by an ODE system. We then use SA for these immuno-epidemiological models to investigate the impact of immunological parameters on population-level disease dynamics such as basic reproduction number, final size of the epidemic or the infectiousness at different phases of an outbreak. As a case study, we consider Rift Valley Fever Disease utilizing parameter estimations from prior studies. SA indicates that [Formula: see text] increase in within-host pathogen growth rate can lead up to [Formula: see text] increase in [Formula: see text] up to [Formula: see text] increase in steady-state infected host abundance, and up to [Formula: see text] increase in infectiousness of hosts when the reproduction number [Formula: see text] is larger than one. These significant increases in population-scale disease quantities suggest that control strategies that reduce the within-host pathogen growth can be important in reducing disease prevalence.

Determining reliable parameter estimates for within-host and within-vector models of Zika virus.

In this paper, we introduce three within-host and one within-vector models of Zika virus. The within-host models are the target cell limited model, the target cell limited model with natural killer (NK) cells class, and a within-host-within-fetus model of a pregnant individual. The within-vector model includes the Zika virus dynamics in the midgut and salivary glands. The within-host models are not structurally identifiable with respect to data on viral load and NK cell counts. After rescaling, the scaled within-host models are locally structurally identifiable. The within-vector model is structurally identifiable with respect to viremia data in the midgut and salivary glands. Using Monte Carlo Simulations, we find that target cell limited model is practically identifiable from data on viremia; the target cell limited model with NK cell class is practically identifiable, except for the rescaled half saturation constant. The within-host-within-fetus model has all fetus-related parameters not practically identifiable without data on the fetus, as well as the rescaled half saturation constant is also not practically identifiable. The remaining parameters are practically identifiable. Finally we find that none of the parameters of the within-vector model is practically identifiable.

Effects of social-distancing on infectious disease dynamics: an evolutionary game theory and economic perspective.

We propose two models inspired by the COVID-19 pandemic: a coupled disease-human behaviour (or disease-game theoretic), and a coupled disease-human behaviour-economic model, both of which account for the impact of social-distancing on disease control and economic growth. The models exhibit rich dynamical behaviour including multistable equilibria, a backward bifurcation, and sustained bounded periodic oscillations. Analyses of the first model suggests that the disease can be eliminated if everybody practices full social-distancing, but the most likely outcome is some level of disease coupled with some level of social-distancing. The same outcome is observed with the second model when the economy is weaker than the social norms to follow health directives. However, if the economy is stronger, it can support some level of social-distancing that can lead to disease elimination.

A Network Immuno-Epidemiological HIV Model.

In this paper we formulate a multi-scale nested immuno-epidemiological model of HIV on complex networks. The system is described by ordinary differential equations coupled with a partial differential equation. First, we prove the existence and uniqueness of the immunological model and then establish the well-posedness of the multi-scale model. We derive an explicit expression of the basic reproduction number [Formula: see text] of the immuno-epidemiological model. The system has a disease-free equilibrium and an endemic equilibrium. The disease-free equilibrium is globally stable when [Formula: see text] and unstable when [Formula: see text]. Numerical simulations suggest that [Formula: see text] increases as the number of nodes in the network increases. Further, we find that for a scale-free network the number of infected individuals at equilibrium is a hump-like function of the within-host reproduction number; however, the dependence becomes monotone if the network has predominantly low connectivity nodes or high connectivity nodes.

Comparison of Drug Inhibitory Effects ([Formula: see text]) in Monolayer and Spheroid Cultures.

Traditionally, the monolayer (two-dimensional) cell cultures are used for initial evaluation of the effectiveness of anticancer drugs. In particular, these experiments provide the [Formula: see text] curves that determine drug concentration that can inhibit growth of a tumor colony by half when compared to the cells grown with no exposure to the drug. Low [Formula: see text] value means that the drug is effective at low concentrations, and thus will show lower systemic toxicity when administered to the patient. However, in these experiments cells are grown in a monolayer, all well exposed to the drug, while in vivo tumors expand as three-dimensional multicellular masses, where inner cells have a limited access to the drug. Therefore, we performed computational studies to compare the [Formula: see text] curves for cells grown as a two-dimensional monolayer and a cross section through a three-dimensional spheroid. Our results identified conditions (drug diffusivity, drug action mechanisms and cell proliferation capabilities) under which these [Formula: see text] curves differ significantly. This will help experimentalists to better determine drug dosage for future in vivo experiments and clinical trials.

Efficacy of control measures in the control of Ebola, Liberia 2014-2015.

The largest outbreak of Ebola to date is the 2014 West Africa Ebola outbreak, with more than 10,000 cases and over 4000 deaths reported in Liberia alone. To control the spread of the outbreak, multiple interventions were implemented: identification and isolation of cases, contact tracing, quarantining of suspected contacts, proper personal protection, safely conducted burials, improved education, social awareness and individual protective measures. Devising rigorous methodologies for the evaluation of the effectiveness of the control measures implemented to stop an outbreak is of paramount importance. In this paper, we evaluate the effectiveness of the 2014 Ebola outbreak interventions. We rely on model selection to determine the best model that explains the 2014 Ebola outbreak data in Liberia which is the simplest model with a social distancing term. We couple structural and practical identifiability analysis with the computation of confidence intervals to pinpoint the uncertainty in the parameter estimations. Finally, we evaluate the efficacy of control measures using the Ebola model with social distancing. Among all the control measures, we find that social distancing had the most impact on the control of the 2014 Ebola epidemic in Libreria followed by isolation and quarantining.

Structural and Practical Identifiability Analysis of Zika Epidemiological Models.

The Zika virus (ZIKV) epidemic has caused an ongoing threat to global health security and spurred new investigations of the virus. Use of epidemiological models for arbovirus diseases can be a powerful tool to assist in prevention and control of the emerging disease. In this article, we introduce six models of ZIKV, beginning with a general vector-borne model and gradually including different transmission routes of ZIKV. These epidemiological models use various combinations of disease transmission (vector and direct) and infectious classes (asymptomatic and pregnant), with addition to loss of immunity being included. The disease-induced death rate is omitted from the models. We test the structural and practical identifiability of the models to find whether unknown model parameters can uniquely be determined. The models were fit to obtain time-series data of cumulative incidences and pregnant infections from the Florida Department of Health Daily Zika Update Reports. The average relative estimation errors (AREs) were computed from the Monte Carlo simulations to further analyze the identifiability of the models. We show that direct transmission rates are not practically identifiable; however, fixed recovery rates improve identifiability overall. We found ARE is low for each model (only slightly higher for those that account for a pregnant class) and help to confirm a reproduction number greater than one at the start of the Florida epidemic. Basic reproduction number, [Formula: see text], is an epidemiologically important threshold value which gives the number of secondary cases generated by one infected individual in a totally susceptible population in duration of infectiousness. Elasticity of the reproduction numbers suggests that the mosquito-to-human ratio, mosquito life span and biting rate have the greatest potential for reducing the reproduction number of Zika, and therefore, corresponding control measures need to be focused on.

Structural and practical identifiability analysis of outbreak models.

Estimating the reproduction number of an emerging infectious disease from an epidemiological data is becoming more essential in evaluating the current status of an outbreak. However, these studies are lacking the fundamental prerequisite in parameter estimation problem, namely the structural identifiability of the epidemic model, which determines the possibility of uniquely determining the model parameters from the epidemic data. In this paper, we perform both structural and practical identifiability analysis to classical epidemic models such as SIR (Susceptible-Infected-Recovered), SEIR (Susceptible-Exposed-Infected-Recovered) and an epidemic model with the treatment class (SITR). We performed structural identifiability analysis on these epidemic models using a differential algebra approach to investigate the well-posedness of the parameter estimation problem. Parameters of these models are estimated from different data types, namely prevalence, cumulative incidences and treated individuals. Furthermore, we carried out practical identifiability analysis on these models using Monte Carlo simulations and Fisher's Information Matrix. Our study shows that the SIR model is both structurally and practically identifiable from the prevalence data. It is also structurally identifiable to cumulative incidence observations, but due to high correlations of the parameters, it is practically unidentifiable from the cumulative incidence data. Furthermore, we found that none of these simple epidemic models are practically identifiable from the cumulative incidence data which is the standard type of epidemiological data provided by CDC or WHO. Our analysis with simple SIR model suggest that the health agencies, if possible, should report prevalence rather than incidence data.

Backward bifurcation and oscillations in a nested immuno-eco-epidemiological model.

This paper introduces a novel partial differential equation immuno-eco-epidemiological model of competition in which one species is affected by a disease while another can compete with it directly and by lowering the first species' immune response to the infection, a mode of competition termed stress-induced competition. When the disease is chronic, and the within-host dynamics are rapid, we reduce the partial differential equation model (PDE) to a three-dimensional ordinary differential equation (ODE) model. The ODE model exhibits backward bifurcation and sustained oscillations caused by the stress-induced competition. Furthermore, the ODE model, although not a special case of the PDE model, is useful for detecting backward bifurcation and oscillations in the PDE model. Backward bifurcation related to stress-induced competition allows the second species to persist for values of its invasion number below one. Furthermore, stress-induced competition leads to destabilization of the coexistence equilibrium and sustained oscillations in the PDE model. We suggest that complex systems such as this one may be studied by appropriately designed simple ODE models.

Efficacies of prevention and control measures applied during an outbreak in Southwest Madrid, Spain.

Leishmaniasis is a vector-borne disease of worldwide distribution, currently present in 98 countries. Since late 2010, an unusual increase of human visceral and cutaneous leishmaniasis cases has been observed in the south-western Madrid region, totaling more than 600 cases until 2015. Some hosts, such as human, domestic dog and cat, rabbit (Oryctolagus cuniculus), and hare (Lepus granatensis), were found infected by the parasite of this disease in the area. Hares were described as the most important reservoir due to their higher prevalence, capacity to infect the vector, and presence of the same strains as in humans. Various measures were adopted to prevent and control the disease, and since 2013 there was a slight decline in the human sickness. We used a mathematical model to evaluate the efficacy of each measure in reducing the number of infected hosts. We identified in the present model that culling both hares and rabbits, without immediate reposition of the animals, was the best measure adopted, decreasing the proportion of all infected hosts. Particularly, culling hares was more efficacious than culling rabbits to reduce the proportion of infected individuals of all hosts. Likewise, lowering vector contact with hares highly influenced the reduction of the proportion of infected hosts. The reduction of the vector density per host in the park decreased the leishmaniasis incidence of hosts in the park and the urban areas. On the other hand, the reduction of the vector density per host of the urban area (humans, dogs and cats) decreased only their affected population, albeit at a higher proportion. The use of insecticide-impregnated collar and vaccination in dogs affected only the infected dogs' population. The parameters related to the vector contact with dog, cat or human do not present a high impact on the other hosts infected by Leishmania. In conclusion, the efficacy of each control strategy was determined, in order to direct future actions in this and in other similar outbreaks. The present mathematical model was able to reproduce the leishmaniasis dynamics in the Madrid outbreak, providing theoretical support based on successful experiences, such as the reduction of human cases in Southwest Madrid, Spain.