A multimorbidity model for estimating health outcomes from the syndemic of injection drug use and associated infections in the United States

Background: Fatal drug overdoses and serious injection-related infections are rising in the US. Multiple concurrent infections in people who inject drugs (PWID) exacerbate poor health outcomes, but little is known about how the synergy among infections compounds clinical outcomes and costs. Injection drug use (IDU) converges multiple epidemics into a syndemic in the US, including opioid use and HIV. Estimated rates of new injection-related infections in the US are limited due to widely varying estimates of the number of PWID in the US, and in the absence of clinical trials and nationally representative longitudinal observational studies of PWID, simulation models provide important insights to policymakers for informed decisions.

Methods: We developed and validated a MultimorbiditY model to Reduce Infections Associated with Drug use (MYRIAD). This microsimulation model of drug use and associated infections (HIV, hepatitis C virus [HCV], and severe bacterial infections) uses inputs derived from published data to estimate national level trends in the US. We used Latin hypercube sampling to calibrate model output against published data from 2015 to 2019 for fatal opioid overdose rates. We internally validated the model for HIV and HCV incidence and bacterial infection hospitalization rates among PWID. We identified best fitting parameter sets that met pre-established goodness-of-fit targets using the Pearson's chi-square test. We externally validated the model by comparing model output to published fatal opioid overdose rates from 2020.

Results: Out of 100 sample parameter sets for opioid use, the model produced 3 sets with well-fitting results to key calibration targets for fatal opioid overdose rates with Pearson's chi-square test ranging from 1.56E-5 to 2.65E-5, and 2 sets that met validation targets. The model produced well-fitting results within validation targets for HIV and HCV incidence and serious bacterial infection hospitalization rates. From 2015 to 2019, the model estimated 120,000 injection-related overdose deaths, 17,000 new HIV infections, and 144,000 new HCV infections among PWID.

Conclusions: This multimorbidity microsimulation model, populated with data from national surveillance data and published literature, accurately replicated fatal opioid overdose, incidence of HIV and HCV, and serious bacterial infections hospitalization rates. The MYRIAD model of IDU could be an important tool to assess clinical and economic outcomes related to IDU behavior and infections with serious morbidity and mortality for PWID.