skip to main content

Title: 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 more » HIV.

« less
; ;
Award ID(s):
Publication Date:
Journal Name:
Mathematical Biosciences and Engineering
Page Range or eLocation-ID:
3636 to 3672
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Drugs of abuse, such as opiates, have been widely associated with enhancing HIV replication, accelerating disease progression and diminishing host-immune responses, thereby making it harder to effectively manage HIV infection. It is thus important to study the effects of drugs of abuse on HIV-infection and immune responses. Here, we develop mathematical models that incorporate the effects of morphine-altered antibody responses on HIV/SIV dynamics. Based on fitting the model to experimental data from simian immunodeficiency virus (SIV) infections in control and morphine-addicted macaques, we found that two of the most significant effects of virus specific antibodies are neutralizing viral particles and enhancing viral clearance. Using our model, we quantified how morphine alters virus-specific antibody responses, and how this alteration affects the key components of virus dynamics such as infection rate, virus clearance, viral load, CD4+T cell count, and CD4+T cell loss in SIV-infected macaques under conditioning with morphine. We found that in a subpopulation of SIV-infected morphine addicted macaques, the presence of drugs of abuse may cause significantly diminished antibody responses, resulting in more severe infection with increased SIV infectivity, a decreased viral clearance rate, increased viral load, and higher CD4+T cell loss.

  2. Abstract In 2016, more than 11 million Americans abused prescription opioids. The National Institute on Drug Abuse considers the opioid crisis a national addiction epidemic, as an increasing number of people are affected each year. Using the framework developed in mathematical modelling of infectious diseases, we create and analyse a compartmental opioid-abuse model consisting of a system of ordinary differential equations. Since $40\%$ of opioid overdoses are caused by prescription opioids, our model includes prescription compartments for the four most commonly prescribed opioids, as well as for the susceptible, addicted and recovered populations. While existing research has focused on drug abuse models in general and opioid models with one prescription compartment, no previous work has been done comparing the roles that the most commonly prescribed opioids have had on the crisis. By combining data from the Substance Abuse and Mental Health Services Administration (which tracked the proportion of people who used or misused one of the four individual opioids) with data from the Centers of Disease Control and Prevention (which counted the total number of prescriptions), we estimate prescription rates and probabilities of addiction for the four most commonly prescribed opioids. Additionally, we perform a sensitivity analysis and reallocate prescriptionsmore »to determine which opioid has the largest impact on the epidemic. Our results indicate that oxycodone prescriptions are both the most likely to lead to addiction and have the largest impact on the size of the epidemic, while hydrocodone prescriptions had the smallest impact.« less
  3. Drug-resistant HIV-1 has caused a growing concern in clinic and public health. Although combination antiretroviral therapy can contribute massively to the suppression of viral loads in patients with HIV-1, it cannot lead to viral eradication. Continuing viral replication during sub-optimal therapy (due to poor adherence or other reasons) may lead to the accumulation of drug resistance mutations, resulting in an increased risk of disease progression. Many studies also suggest that events occurring during the early stage of HIV-1 infection (i.e., the first few hours to days following HIV exposure) may determine whether the infection can be successfully established. However, the numbers of infected cells and viruses during the early stage are extremely low and stochasticity may play a critical role in dictating the fate of infection. In this paper, we use stochastic models to investigate viral infection and the emergence of drug resistance of HIV-1. The stochastic model is formulated by a continuous-time Markov chain (CTMC), which is derived based on an ordinary differential equation model proposed by Kitayimbwa et al. that includes both forward and backward mutations. An analytic estimate of the probability of the clearance of HIV infection of the CTMC model near the infection-free equilibrium is obtainedmore »by a multitype branching process approximation. The analytical predictions are validated by numerical simulations. Unlike the deterministic dynamics where the basic reproduction number $ \mathcal{R}_0 $ serves as a sharp threshold parameter (i.e., the disease dies out if $ \mathcal{R}_0 < 1 $ and persists if $ \mathcal{R}_0 > 1 $), the stochastic models indicate that there is always a positive probability for HIV infection to be eradicated in patients. In the presence of antiretroviral therapy, our results show that the chance of clearance of the infection tends to increase although drug resistance is likely to emerge.

    « less
  4. Abstract This project is funded by the US National Science Foundation (NSF) through their NSF RAPID program under the title “Modeling Corona Spread Using Big Data Analytics.” The project is a joint effort between the Department of Computer & Electrical Engineering and Computer Science at FAU and a research group from LexisNexis Risk Solutions. The novel coronavirus Covid-19 originated in China in early December 2019 and has rapidly spread to many countries around the globe, with the number of confirmed cases increasing every day. Covid-19 is officially a pandemic. It is a novel infection with serious clinical manifestations, including death, and it has reached at least 124 countries and territories. Although the ultimate course and impact of Covid-19 are uncertain, it is not merely possible but likely that the disease will produce enough severe illness to overwhelm the worldwide health care infrastructure. Emerging viral pandemics can place extraordinary and sustained demands on public health and health systems and on providers of essential community services. Modeling the Covid-19 pandemic spread is challenging. But there are data that can be used to project resource demands. Estimates of the reproductive number (R) of SARS-CoV-2 show that at the beginning of the epidemic, each infectedmore »person spreads the virus to at least two others, on average (Emanuel et al. in N Engl J Med. 2020, Livingston and Bucher in JAMA 323(14):1335, 2020). A conservatively low estimate is that 5 % of the population could become infected within 3 months. Preliminary data from China and Italy regarding the distribution of case severity and fatality vary widely (Wu and McGoogan in JAMA 323(13):1239–42, 2020). A recent large-scale analysis from China suggests that 80 % of those infected either are asymptomatic or have mild symptoms; a finding that implies that demand for advanced medical services might apply to only 20 % of the total infected. Of patients infected with Covid-19, about 15 % have severe illness and 5 % have critical illness (Emanuel et al. in N Engl J Med. 2020). Overall, mortality ranges from 0.25 % to as high as 3.0 % (Emanuel et al. in N Engl J Med. 2020, Wilson et al. in Emerg Infect Dis 26(6):1339, 2020). Case fatality rates are much higher for vulnerable populations, such as persons over the age of 80 years (> 14 %) and those with coexisting conditions (10 % for those with cardiovascular disease and 7 % for those with diabetes) (Emanuel et al. in N Engl J Med. 2020). Overall, Covid-19 is substantially deadlier than seasonal influenza, which has a mortality of roughly 0.1 %. Public health efforts depend heavily on predicting how diseases such as those caused by Covid-19 spread across the globe. During the early days of a new outbreak, when reliable data are still scarce, researchers turn to mathematical models that can predict where people who could be infected are going and how likely they are to bring the disease with them. These computational methods use known statistical equations that calculate the probability of individuals transmitting the illness. Modern computational power allows these models to quickly incorporate multiple inputs, such as a given disease’s ability to pass from person to person and the movement patterns of potentially infected people traveling by air and land. This process sometimes involves making assumptions about unknown factors, such as an individual’s exact travel pattern. By plugging in different possible versions of each input, however, researchers can update the models as new information becomes available and compare their results to observed patterns for the illness. In this paper we describe the development a model of Corona spread by using innovative big data analytics techniques and tools. We leveraged our experience from research in modeling Ebola spread (Shaw et al. Modeling Ebola Spread and Using HPCC/KEL System. In: Big Data Technologies and Applications 2016 (pp. 347-385). Springer, Cham) to successfully model Corona spread, we will obtain new results, and help in reducing the number of Corona patients. We closely collaborated with LexisNexis, which is a leading US data analytics company and a member of our NSF I/UCRC for Advanced Knowledge Enablement. The lack of a comprehensive view and informative analysis of the status of the pandemic can also cause panic and instability within society. Our work proposes the HPCC Systems Covid-19 tracker, which provides a multi-level view of the pandemic with the informative virus spreading indicators in a timely manner. The system embeds a classical epidemiological model known as SIR and spreading indicators based on causal model. The data solution of the tracker is built on top of the Big Data processing platform HPCC Systems, from ingesting and tracking of various data sources to fast delivery of the data to the public. The HPCC Systems Covid-19 tracker presents the Covid-19 data on a daily, weekly, and cumulative basis up to global-level and down to the county-level. It also provides statistical analysis for each level such as new cases per 100,000 population. The primary analysis such as Contagion Risk and Infection State is based on causal model with a seven-day sliding window. Our work has been released as a publicly available website to the world and attracted a great volume of traffic. The project is open-sourced and available on GitHub. The system was developed on the LexisNexis HPCC Systems, which is briefly described in the paper.« less
  5. Co-infections of hosts by multiple pathogen species are ubiquitous, but predicting their impact on disease remains challenging. Interactions between co-infecting pathogens within hosts can alter pathogen transmission, with the impact on transmission typically dependent on the relative arrival order of pathogens within hosts (within-host priority effects). However, it is unclear how these within-host priority effects influence multi-pathogen epidemics, particularly when the arrival order of pathogens at the host-population scale varies. Here, we combined models and experiments with zooplankton and their naturally co-occurring fungal and bacterial pathogens to examine how within-host priority effects influence multi-pathogen epidemics. Epidemiological models parametrized with within-host priority effects measured at the single-host scale predicted that advancing the start date of bacterial epidemics relative to fungal epidemics would decrease the mean bacterial prevalence in a multi-pathogen setting, while models without within-host priority effects predicted the opposite effect. We tested these predictions with experimental multi-pathogen epidemics. Empirical dynamics matched predictions from the model including within-host priority effects, providing evidence that within-host priority effects influenced epidemic dynamics. Overall, within-host priority effects may be a key element of predicting multi-pathogen epidemic dynamics in the future, particularly as shifting disease phenology alters the order of infection within hosts.