Abstract

Opioid use disorder (OUD) has become a serious leading health issue in the USA leading to addiction, disability, or death by overdose. Research has shown that OUD can lead to a chronic lifelong disorder with greater risk for relapse and accidental overdose deaths. While the prescription opioid epidemic is a relatively new phenomenon, illicit opioid use via heroin has been around for decades. Recently, additional illicit opioids such as fentanyl have become increasingly available and problematic. We propose a mathematical model that focuses on illicit OUD and includes a class for recovered users but allows for individuals to either remain in or relapse back to the illicit OUD class. Therefore, in our model, individuals may cycle in and out of three different classes: illicit OUD, treatment, and recovered. We additionally include a treatment function with saturation, as it has been shown there is limited accessibility to specialty treatment facilities. We used 2002–2019 SAMHSA and CDC data for the US population, scaled to a medium-sized city, to obtain parameter estimates for the specific case of heroin. We found that the overdose death rate has been increasing linearly since around 2011, likely due to the increased presence of fentanyl in the heroin supply. Extrapolation of this overdose death rate, together with the obtained parameter estimates, predict that by 2038 no endemic equilibrium will exist and the only stable equilibrium will correspond to the absence of heroin use disorder in the population. There is a range of parameter values that will give rise to a backward bifurcation above a critical saturation of treatment availability. We show this for a range of overdose death rate values, thus illustrating the critical role played by the availability of specialty treatment facilities. Sensitivity analysis consistently shows the significant role of people entering treatment on their own accord, which suggests the importance of removing two of the most prevalent SAMHSA-determined reasons that individuals do not enter treatment: financial constraints and the stigma of seeking treatment for heroin use disorder.

Highlights

  • A national crisis has emerged regarding opioid use disorder (OUD) (Vivolo-Kantor et al 2018)

  • We propose an “illicit opioid use disorder” (IOUD) model to describe the role that black market opioids such as heroin, fentanyl, and other synthetic opioids play in the current opioid epidemic

  • Our paper presents a deterministic model for the dynamics of an illicit opioid use disorder (IOUD) model

Read more

Summary

Introduction

A national crisis has emerged regarding opioid use disorder (OUD) (Vivolo-Kantor et al 2018). Prescription opioid overdose, abuse, and dependence accounts for a total cost of 78.5 billion dollars a year reported by the Centers for Disease Control and Prevention (CDC) These costs are the result of elevated health care, drug abuse treatment, criminal justice, and loss of productivity expenditures (National Institutes of Health 2019; Florence et al 2016). Opioid analgesics have been overprescribed and given their effect on the brain, this has resulted in an increased risk of OUD This has influenced an increase of heroin use where multiple users (4 out of 5 reported) have switched over from opioid pain reliever prescriptions because of lower cost and accessibility issues (Kolodny et al 2015; Volkow 2018; Schuckit 2016; Connery 2015). A tool that can be used for understanding the complex issues surrounding OUD and illicit OUD, its treatment options, and methods for decreasing relapse, is a mathemat-

48 Page 4 of 49
Model Formulation and Basic Properties
48 Page 6 of 49
Model Equations
48 Page 8 of 49
Data Explanation and Parameter Estimation
48 Page 10 of 49
48 Page 12 of 49
48 Page 14 of 49
Steady-State Analysis
Calculating the Basic Reproductive Number R0
48 Page 16 of 49
Endemic Equilibria
48 Page 18 of 49
48 Page 20 of 49
Sensitivity Analysis
Discussion of the PRCC Values
48 Page 24 of 49
48 Page 30 of 49
Conclusion
48 Page 32 of 49
48 Page 34 of 49
No Backward Bifurcation for = ̨2 = 0
Findings
Additional PRCC Plots and Discussion
48 Page 48 of 49
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call