Abstract
Background: Illicit opioid overdose continues to rise in North America and is a leading cause of death. Mathematical modeling is a valuable tool to investigate the epidemiology of this public health issue, as it can characterize key features of population outcomes and quantify the broader effect of structural and interventional changes on overdose mortality.
The aim of this study is to quantify and predict the impact of key harm reduction strategies at differing levels of scale-up on fatal and nonfatal overdose among a population of people who inject drugs (PWID) in Toronto.
Methods: An individual-based model of PWID was built featuring demographic and behavioural variation among members. Key individual attributes known to scale the risk of fatal and nonfatal overdose were identified and incorporated into a dynamic modeling framework, wherein every member of the simulated population encompasses a set of distinct characteristics that govern demographics, intervention usage, and overdose incidence. The model was parametrized to fatal and nonfatal overdose events reported in Toronto in 2019. The interventions considered were opioid agonist therapy (OAT), supervised consumption sites (SCS), take-home naloxone (THN), drug-checking, and reducing fentanyl in the drug supply. Harm reduction scenarios were explored relative to a baseline model to examine the impact of each intervention being scaled from 0% use to 100% use on overdose events.
Results: Model simulations resulted in 3690.6 nonfatal and 295.4 fatal overdoses, coinciding with 2019 data from Toronto. From this baseline, at full scale-up, 290 deaths were averted by THN, 248 from eliminating fentanyl from the drug supply, 124 from SCS use, 173 from OAT, and 100 by drug-checking services. Drug-checking and reducing fentanyl in the drug supply were the only harm reduction strategies that reduced the number of nonfatal overdoses.
Conclusions: Within a multi-faceted harm reduction approach, scaling up take-home naloxone, and reducing fentanyl in the drug supply led to the largest reduction in overdose fatality among PWID in Toronto. Detailed model simulation studies provide an additional tool to assess and inform public health policy on harm reduction.