Next Article in Journal
The Importance of the Concentration of Selected Cytokines (IL-6, IL-10, IL-12, IL-15, TNF-α) and Inflammatory Markers (CRP, NLR, PLR, LMR, SII) in Predicting the Course of Rehabilitation for Patients after COVID-19 Infection
Previous Article in Journal
Resting-State EEG Alterations of Practice-Related Spectral Activity and Connectivity Patterns in Depression
Previous Article in Special Issue
Basivertebral Nerve Ablation for Treatment of Lower Back Pain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predictive Model for Opioid Use Disorder in Chronic Pain: A Development and Validation Study

1
Pharmacogenetics Unit, Clinical Pharmacology Department, Alicante Institute for Health and Biomedical Research (ISABIAL), Diagnostic Center, Gray Building, 5th Floor, Avda Pintor Baeza, 12, 03010 Alicante, Spain
2
Department of Pharmacology, Paediatrics and Organic Chemistry, Miguel Hernández University of Elche, 03550 Sant Joan de Alicante, Spain
3
Pain Unit, Dr. Balmis General University Hospital, ISABIAL, 03010 Alicante, Spain
4
Department of Psychology and Cognitive Science, University of Trento, 38060 Trento, Italy
5
Operations Research Center, University Miguel Hernández de Elche, 03202 Elche, Spain
6
Institute of Bioengineering, Miguel Hernández University, 03202 Elche, Spain
*
Author to whom correspondence should be addressed.
Biomedicines 2024, 12(9), 2056; https://doi.org/10.3390/biomedicines12092056
Submission received: 23 August 2024 / Revised: 2 September 2024 / Accepted: 4 September 2024 / Published: 10 September 2024

Abstract

:
Background/Objective: There are several questionnaires for the challenge of anticipating opioid use disorder (OUD). However, many are not specific for chronic non-cancer pain (CNCP) or have been developed in the American population, whose sociodemographic factors are very different from the Spanish population, leading to scarce translation into clinical practice. Thus, the aim of this study is to prospectively validate a predictive model for OUD in Spanish patients under long-term opioids. Methods: An innovative two-stage predictive model was developed from retrospective (n = 129) and non-overlapping prospective (n = 100) cohorts of real-world CNCP outpatients. All subjects used prescribed opioids for 6 or more months. Sociodemographic, clinical and pharmacological covariates were registered. Mu-opioid receptor 1 (OPRM1, A118G, rs1799971) and catechol-O-methyltransferase (COMT, G472A, rs4680) genetic variants plus cytochrome P450 2D6 (CYP2D6) liver enzyme phenotypes were also analyzed. The model performance and diagnostic accuracy were calculated. Results: The two-stage model comprised risk factors related to OUD (younger age, work disability and high daily opioid dose) and provided new useful information about other risk factors (low quality of life, OPRM-G allele and CYP2D6 extreme phenotypes). The validation showed a satisfactory accuracy (70% specificity and 75% sensitivity) for our predictive model with acceptable discrimination and goodness of fit. Conclusions: Our study presents the results of an innovative model for predicting OUD in our setting. After external validation, it could represent a change in the paradigm of opioid treatment, helping clinicians to better identify and manage the risks and reduce the side effects and complications.

1. Introduction

Opioid analgesics are recognized as a legitimate medical therapy for selected patients with severe chronic non-cancer pain (CNCP) that does not respond to other therapies [1]. In Spain, the use of opioid analgesics has increased by almost 54% overall from 2013 to 2019. However, opioids are associated with risks that include aberrant drug-related behaviors and opioid use disorder (OUD) in up to 20% of regular opioid prescriptions [2]. OUD is detailed in DSM-5 as an unsuccessful effort to cut down or control use, social problems and a failure to fulfill major role obligations [3]. Thus, clinical tools, greater understanding and better assessment are needed to properly screen and monitor patients under long-term opioid prescriptions.
There is not enough understanding of the inter-individual variability in analgesic administration and aberrant opioid-related behaviors [4]. The scientific literature has indicated that some genetic polymorphisms (mu-opioid receptor 1-OPRM1 (A118G, rs1799971), CYP2D6 (10 single nucleotide polymorphisms and copy number variation) and cathecol-O-methyltransferase-COMT (G472A, rs4680), among others) may as well contribute to inter-individual differences in morphine consumption [5,6].
In our Pain Unit (PU), a predictive model for OUD [7] was developed from a retrospective cohort of patients under chronic use of opioids [8,9,10,11]. This model included these actionable pharmacogenetic markers and other covariates such as age, employment status and equivalent morphine daily dose (MEDD). In this way, the aim of this study was to internally validate the model, in real-world CNCP outpatients, for future clinical translation. We also examined the characteristics of routinely ambulatory CNCP patients under long-term use of opioids.

2. Materials and Methods

2.1. Patients

An observational cross-sectional study was conducted on patients who routinely attended the PU of Dr. Balmis General University Hospital. This study was approved on 15 May 2020, by the Institutional Review Board (code PI2020-047) where procedures were carried out in accordance with the Declaration of Helsinki. Enrolment began in September 2021 and ended in July 2022. All patients included were ≥18 years old with CNCP (moderate or severe pain lasting for six or more months) under long-term opioids (≥6 months). The exclusion criteria were oncologic pain, opioid prescription <6 months and prior inclusion in the retrospective cohort (model development). Informed consent was obtained from all subjects involved in this study. This manuscript adheres to the applicable STROBE guidelines.

2.2. Procedure and Variables

All subjects enrolled were attended by the research staff for data and saliva sample collection. The variables were collected, in both groups (retrospective and prospective), at the time of the enrolment. They were collected through patient self-report, validated scales and questionnaires (described below) and completed through Electronic Health Records (EHRs). After data collection, patients followed their routine clinical visit. Here, a medical doctor (anaesthesiologist or clinical pharmacologist expert on pain) assessed the OUD diagnosis based on DSM-5 criteria [12].
Sex (female/male), age, employment status (yes/no: active, retired, with work disability-permanent or temporary, unemployed or homemaker) were registered. The cut-off points for monthly incomes were established according to the Spanish minimum interprofessional wage (EUR 1000) and the minimum vital income (EUR 500) to facilitate the translation to other countries. Thus, data were categorized into low incomes—less than EUR 500, middle incomes—between EUR 500–1000 or upper incomes—more than EUR 1000.
The presence/absence of current and/or previous substance use disorders (SUDs) (except opioid use) related to tobacco, alcohol or other illicit drugs were collected through the review of medical diagnoses, narratives or any visit to the Addictive Behavior Unit.
Pain intensity and relief and quality of life were measured using the Visual Analogue Scale (VAS) [13] included in the Global Pain State questionnaire [14], already validated. The VAS for each indicator consists of a 100 mm horizontal line ranging from 0 (lowest) to 100 mm (highest). Specifically, quality of life was measured with the validated EuroQol-5D-3L scale (registration number: 48802, available at https://euroqol.org/, accessed on 30 May 2021), which also includes a Health Utility Score [15] (five dimensions: mobility, self-care, usual activities, pain/discomfort and anxiety/depression; scores from 0 for death to 1 for perfect health). In addition, the health resources use (any recent emergency department (ED) visit, hospitalization or drug changes due to pain and other causes) was registered.
All and only prescribed drug use was registered and contrasted with EHRs, which allows for reviewing drug prescriptions. Non-opioid analgesics (i.e., paracetamol and metamizole), non-steroidal anti-inflammatory drugs (NSAIDs), weak (i.e., tramadol and codeine) and strong opioids use (i.e., fentanyl, oxycodone, tapentadol, buprenorphine, morphine, hydromorphone and methadone) and immediate release opioids were registered. Buprenorphine, in our clinical setting, is also used for CNCP patients due to its versatility in administration and its manageable side effects. In different opioid combinations, oral MEDD was estimated using available references [16]. The prescription of antidepressants (i.e., amitriptyline, duloxetine and escitalopram), benzodiazepines and neuromodulators (i.e., pregabalin, lacosamide, gabapentin) were also collected.
What is more, patients’ reports of adverse events (AEs) were collected through a list with the most frequent adverse drug reactions (ADRs, selected according to opioid Summary of Product Characteristics frequency as “very common” and “common”) [17] and a blank space to collect any other adverse event presented. In addition, patients were asked about any depression or anxiety symptoms. They were also grouped by systems according to the Medical Dictionary for Regulatory Activities Terminology—MedDRA (available at https://www.meddra.org, accessed on 14 June 2022) [18,19].

2.3. Pharmacogenetic Analysis

Approximately 2 mL of saliva was collected in tubes containing 6 mL of PBS. Once the saliva sample was taken, it was stored at −80 °C until its processing. Genomic DNA was extracted using E.N.Z.A. forensic DNA kit (Omega Bio-Tek Inc., Norcross, GA, USA) following the manufacturer’s instructions. The following gene variants were genotyped: OPRM1 (rs1799971, A118G), COMT (rs4680, G472A) and CYP2D6*2, *3, *4, *5, *6, *10, *17, *29, *35, *41, xN using the real-time PCR rotor gene Q system (Qiagen, Hilden, Germany), through the use of specific TaqMan MGB® probes (Applied Biosystems, Waltham, MA, USA). Amplification parameters were as follows: pre-PCR section 10 min at 95 °C, 40 cycles for 15 seconds of denaturation at 92 °C and 1-min final extension at 60 °C. The specifications regarding conversion from genotype to phenotype are in Appendix A.

2.4. Statistical Methods

It was expected to have 20 cases and 100 controls from the PU for 11 months based on the inclusion rate of cases in the previous study [8]. We anticipated withdrawals, incomplete data or losses to teleassistance instead of PU visits, shrinking the targeted enrolment to 100 patients.
We compared all the variables between the retrospective cohort (Sample 1) and the new cohort (prospective cohort, Sample 2) using χ2 or Fisher’s exact test for categorical variables and t-test or U Mann–Whitney test for continuous variables, depending upon their distribution. Here, data distribution was analyzed with the Kolmogorov–Smirnov test using the Lilliefors correction method. Quantitative parametric data (pain intensity, relief and quality of life) are presented as mean and standard deviation (SD), whilst the median and interquartile range (IQR) were used for not parametric data (age, quality of life—Health Utility Score, MEDD, AEs). Categorical data (sex, employment status, income, prior SUD, health resource use, drug use, AEs and genotypes) are expressed as percentages (%). Gene frequencies were compared using the chi-squared χ2 goodness-of-fit test. For the OPRM1 genotype, the G-carriers (subjects with AG or GG genotypes) were grouped due to the limited G allelic frequency.
We applied the logistic regression model previously developed from Sample 1 (Table S1) to calculate the risk of each individual as follows: eζ/(1 + eζ), where the linear predictor ζ = b0 + b1x1 + b2x2 + … + bpxp contains five independent risk factors. In other words, ζ = 1.633 − 0.072 age + 2.012 work disability + 0.006 MEDD − 1.424 OPRM1 genotype (AG/GG) + 0.075 PM CYP2D6 phenotype + 3.172 UM CYP2D6 phenotype.
The methodology proposed for the fitting and internal validation of the model was as follows: First, the total sample was randomly divided into two parts: 80% (for model fitting) and 20% (for validation). From here, a two-stage model was proposed for the model adjustment due to the differences observed for (1) classifying patients in Sample 1 or 2 and (2) predicting OUD risk in each Sample (1 and 2). Here, variables were selected for the model on the basis of the investigators’ consensus on relevant measurable variables, the results of previous studies [8,9] and the univariate analysis (p < 0.05). Three logistic regression models were constructed based on the standards for the model-building process [20]. The model selection followed two criteria: (1) small Akaike information criterion—AIC and (2) significance of the variables. Non-significant variables were considered if their coefficients were interpretable. On the other hand, the validation was performed on the observations not used (20%) in fitting the model. This percentage was considered in order to more accurately estimate the clinical usefulness (sensitivity and specificity). Calibration (Hosmer–Lemeshow goodness-of-fit statistic) and discrimination (C-statistic, area under the receiver operating curve) were measured to assess the model performance. All statistical analyses were carried out using R (Version 3.2.0; the GNU project, Cambridge, MA, USA) and GraphPad Prism (Version 5.0, Dotmatics, Boston, MA, USA).

3. Results

3.1. Participants and Variables

Sample 1 (n = 129) corresponds to the retrospective cohort of patients used originally for the model development (7). Sample 2 (n = 100) corresponds to the new prospective cohort of patients recruited, following the detailed inclusion criteria. Three subjects were missing due to insufficient samples for the pharmacogenetic analysis (two for full analysis and one for CYP2D6 analysis) (Figure 1). All patients enrolled attended our PU for regular CNCP management due to lumbalgia (67%, mostly for disc disease pain from spinal canal stenosis with or without a radicular or myofascial pain component), knee pain (gonalgia) and other musculoskeletal pain (i.e., arthralgia and cervical joint dysfunctions).
Middle-aged (63–65 years old), predominantly female (67–70%), retired (50–40%) with 18–25% of a previous SUD (mostly tobacco 71–96%) were the main characteristics of the population included. Half (53–67%) presented middle incomes; however, a higher prevalence of upper incomes (42% vs. 13%, p = 0.04) and lower tobacco use (71% vs. 96%, p = 0.03) was evidenced in Sample 2, as seen in Table 1.
Patients of the prospective cohort (Sample 2) suffered higher pain intensity (70 (26) vs. 61 (28) mm, p = 0.02) with the near double use of non-opioid analgesics (63% vs. 34%, p < 0.001) and tramadol (45% vs. 22%, p < 0.001). However, less MEDD (median (IQR), 60 (33–108) vs. 80 (40–160) mg/day, p < 0.01) and immediate-release opioids prescription (10% vs. 24%, p < 0.01) but higher benzodiazepines use (54% vs. 35%, p < 0.01) was observed in front of the retrospective cohort (Sample 1). A different oxycodone, tapentadol and buprenorphine use was also observed between Samples. In addition, the number of AEs reported (2 (1–3) vs. 6 (3–8), p < 0.001) were significantly lower in Sample 2 (Table S3).

3.2. Model Performance

Due to the differences observed between Sample 1 and 2, the logistic regression model previously developed from Sample 1 did not present good sensitivity (fewer false negatives, 14%) in the new cohort of patients (Sample 2). In this way, a two-stage model was proposed for the adjustment of the developed model.
A total of nineteen variables were selected according to the established criteria (see in statistical analysis) and entered in the logistic regression models as candidate predictors: age, employment status (active, work disability and unemployed), prior SUD, pain intensity, quality of life, tramadol use, MEDD, strong opioids use, fentanyl use, benzodiazepines use, ED visits, vomiting, sleep disturbance, psychiatric AEs, OPRM1 genotype (AA, AG/GG), COMT genotype (GG, GA and AA) and CYP2D6 phenotypes (PM, EM and UM) (Table S2)) [9,21,22,23,24,25].
Firstly, a logistic regression model was developed to classify patients in Sample 1 or 2 (Table 2). This model included nine independent factors: ζ = 0.242 − 1.950 active − 1.740 work disability − 3.976 unemployed + 0.004 MEDD + 3.493 strong opioid use − 2.626 benzodiazepines use − 1.496 ED visits + 2.289 psychiatric AEs − 2.159 COMT genotype (GA) − 0.901 COMT genotype (GG). Here, the cut-off point (c = 0.57) presented the optimal values for specificity (0.86) and sensitivity (0.85).
Secondly, two logistic regression models were developed to estimate OUD risk in each Sample. The newly developed model for Sample 1 (Table 3) included four independent factors: ζ = −0.622 − 0.057 age + 2.859 work disability + 0.006 MEDD + 1.191 PM CYP2D6 phenotype + 3.299 UM CYP2D6 phenotype. Here, the optimal values of specificity (0.82) and sensitivity (0.94) were obtained with a cut-off point of 0.18. The C-statistic indicated a satisfactory model discrimination (0.86). The model’s ability to accurately predict the likelihood of developing OUD was measured with the test Hosmer–Lemeshow (p = 0.26), which indicated a limited model fit.
The predictive model for Sample 2 (Table 4) included three independent factors: ζ = −1.713 − 0.032 quality of life + 0.006 MEDD + 1.017 OPRM1 genotype (AG/GG) with an optimal cut-off point of 0.19 for satisfactory specificity (0.78) and sensitivity (0.73). The C-statistic indicated a satisfactory model discrimination (0.86). The Hosmer–Lemeshow (p = 0.36) showed a limited model fit.

3.3. Model Validation

The model developed for classifying patients in Sample 1 or 2 had a satisfactory sensitivity and specificity (0.78 and 0.68, respectively) with a cut-off point of 0.35. On the other hand, the models for predicting OUD risk presented adequate sensitivities (0.75 both Samples) and specificities (0.81 and 0.57, Sample 1 and Sample 2) for a cut-off point of 0.08 and 0.10, respectively. Thus, the two-stage model presented on average 70% specificity and 75% sensitivity.

4. Discussion

We have developed and internally validated a predictive model as a screening tool that can classify CNCP patients at risk for OUD when they are under long-term opioids. This tool consisted of a two-stage model that comprised well-documented risk factors related to OUD (younger age, work disability and high MEDD) and provided more useful information about other less-explored risk factors (low quality of life, OPRM-G allele and CYP2D6 extreme phenotypes).
One of the novel features of this study is the inclusion of genetic variables in the predictive model, which is lacking in many prior tools, including the Opioid Risk Tool. In this era of precision medicine and artificial intelligence, healthcare could benefit from such studies that utilize genetic predictors to stratify patients into risk categories for OUD [26]. What is more, while previous studies have developed predictive models with databases from specific populations (e.g., a veterans’ health administration database) [21,27], our model used regular CNCP patients from ambulatory PU data.
Clinical guidelines have been established and recommend CYP2D6 genotype testing prior to prescription of tramadol or codeine, as it has been associated with failure of pain treatment in PM (limited conversion to active metabolites) and a higher risk of AEs in UM [22]. In the case of OPRM1, numerous studies have associated the mutant variant (118G) with OUD risk, as they have observed a lower receptor expression in the membrane [23]. In this way, pharmacogenetic testing would allow healthcare providers to individualize prevention strategies [24]
It is well-recognized that MEDD is a major determinant for developing an OUD [25]. Experts have agreed that lower doses of opioids could reduce the risk of opioid use disorder and overdose [3]. Here, daily doses close to or greater than 100 mg/day are at higher risk than dosages < 50 mg/day. Yang S. Liu et al. [28] developed and validated an OUD risk predictive model in 316,039 patients from a national healthcare database, where MEDD was one of the ten top-ranked predictors. Nevertheless, they lacked a clear indicator for treated OUD, and the model interpretability needed further validation.
Data reveal that the opioid crisis disproportionately impacts some specific populations, such as people with low incomes, with past or current substance abuse and untreated psychiatric disorders [1]. Our data reveal that younger ages with a more vulnerable work status can condition OUD, together with a low quality of life. The latter has been previously linked with SUDs [29]. Thus, recognizing and measuring it in clinical practice should also improve the outcomes in patients with OUD.
These results need to be interpreted with caution due to their limitations. Firstly, the relatively poor incidence of OUD in our setting, being the result of clinical practice, could have prevented us from detecting the causality and other potential risk factors. Here, it is relevant to mention that patients came from a PU, where the prevalence of OUD is not even described. However, based on the sensitivity and specificity results observed, we could conclude that the predictive capacity of the model is adequate. In this study, we proposed a decision tree, which can be very useful in subpopulations with differences.
Among other limitations, it is not clear how strong the influence of genetic factors is compared to other factors in predicting OUD. Moreover, the allele frequency of CYP2D6 extreme phenotypes could have limited the detection of OUD risk with more accuracy, although the results showed an adequate sample size. In addition, we should have considered the different metabolism ways related to each opioid drug used and coadjuvants with potential inhibitory or enhancing effects, as they can have different addiction potentials and a different affinity to opioid receptors and therefore different behavioral effects. Prior SUDs were registered from the EHRs, which is limited by the missing information reported by clinicians. However, it has been reported that the prevalence of this variable is lower in prescribed opioid users [30]. Additionally, we only focused on musculoskeletal pain, avoiding other types of pain. In addition, the family history of substance use disorder together with the personality traits could have provided more comprehensive information regarding the risk profile, improving the model accuracy.
Finally, the three genes studied have biologically plausible mechanisms for affecting opioid response (CYP2D6—drug metabolism, OPRM1—drug target, COMT—pain perception). However, additional drug–gene associations have preliminary evidence (ABCB1-morphine, CYP3A4-fentanyl, CYP2B6-methadone) [31] and should be explored in the next study. Here, it is important to highlight that pharmacogenomic screening is not part of routine clinical care yet, but the increasing ease of the availability of genetic data due to the improvement and reduction in the cost of the techniques available for its analysis is progressively helping its implementation [32]. In addition, on 23 June 2023, the National Health Services proposed the inclusion of pharmacogenomic testing in the genetic services. The proposal included a minimum of 12 genes to cover a total of 65 drugs—23 for primary use and 42 for secondary use. In this way, with a blood or saliva sample, we can obtain information for more than 70% of the most common drugs, obtaining a high level of cost-effectiveness. The rest of the variables (clinical and pharmacological) included in the model are easily accessible and captured in clinical care. Thus, it must also be taken into account that there are other factors that influence the variability of the response to drugs, such as some environmental determinants. It is important to standardize the pharmacogenetic report so that it is incorporated together with the rest of the variables. In this way, the implementation of pharmacogenetics must involve multidisciplinary collaboration between health professionals, which allows the assessment of the response to the treatment of each patient taking into account their context, and other clinical variables.
The ultimate purpose of this model is to help identify patients who are prone to develop OUD and who therefore should be closely monitored or may benefit from preventive interventions. This model could only be applied in our center, where it has been developed, and for a certain time. As a future research line for external validation, it is needed to carry out a cross-validation procedure to obtain prediction intervals. At this point, it is crucial to study a generalization of the model in other patient populations and clinical settings, minimizing model overfitting.

5. Conclusions

We have developed an innovative predictive model based on real-world data from the Pain Unit routine clinical care in Spain. This model could help to focus on patients requiring monitoring or preventive interventions, optimizing medical resources and improving patients’ quality of life. In this way, external validation is crucial to ensure the clinical usefulness of the model in diverse patient populations and clinical settings.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines12092056/s1, Supplementary Data: CYP2D6-Conversion; Table S1: Original model developed from Sample 1; Table S2: Justification for the inclusion of specific predictors in the model; Table S3: Safety characteristics of Samples 1 and 2.

Author Contributions

Conceptualization, A.M.P., T.Z. and J.M.; methodology, D.M., M.E. and C.M.; software, D.M. and M.E.; validation, D.M. and M.E.; formal analysis, M.E. and J.M.; investigation, M.E., L.A., C.M. and A.P.; resources, D.M. and A.M.P.; data curation, M.E., J.M. and L.A.; writing—original draft preparation, A.M.P., J.M., M.E. and T.Z.; writing—review and editing, A.M.P., J.M., D.M., A.P., C.M., L.A., M.E. and T.Z.; visualization, A.M.P., M.E. and J.M.; supervision, A.M.P. and D.M.; project administration, A.M.P. and T.Z.; funding acquisition, M.E. and A.M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by three grants (Alicante Institute for Health and Biomedical Research, code UGP 21-116; Navarro-Tripodi Foundation, 2021 Research Grant; and National Drugs Plan, code 2023I031) and one research award from the Spanish Clinical Pharmacology Society, 2022 Research Award.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of DR BALMIS ALICANTE GENERAL HOSPITAL (PI1020-047, 15 May 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Regarding the CYP2D6 genotype, genetic analysis was based on the usual PCR methods following the instructions of Consortium of the Pharmacogenetics and Pharmacogenomics Ibero-American network 21. These amplifications were carried out in a Mastercycler 384 (Eppendorf, Hamburg, Germany).
The combination of CYP2D6 alleles is used to determine a patient’s diplotype. Each allele is assigned an activity score (AS) ranging from 0 to 1 (e.g., 0 for no function, 0.25 or 0.5 for decreased function and 1 for normal function) 22. Alleles *3, *4, *5, *6 (AS = 0) have null function, while *10 (AS = 0.25) and *17, *29, *41 (AS = 0.5) a reduced function. Alleles *1, *2, *35 (AS = 1) and duplications *1xN, *2xN, *35xN (AS = 2) have normal and greater function, respectively. If an allele contains multiple copies of a functional gene, the value is multiplied by the number of copies present. Thus, the CYP2D6 activity score is the sum of the values assigned to each allele. The CYP2D6 activity score can be translated into a standardized phenotype classification: null function (poor metabolizer, PM), normal function (extensive metabolizer, EM) and increased function (ultra-rapid metabolizers, UM) [22].

References

  1. Webster, L.R. Risk Factors for Opioid-Use Disorder and Overdose. Anesth. Analg. 2017, 125, 1741–1748. [Google Scholar] [CrossRef] [PubMed]
  2. Salsitz, E.A. Chronic Pain, Chronic Opioid Addiction: A Complex Nexus. J. Med. Toxicol. 2016, 12, 54–57. [Google Scholar] [CrossRef] [PubMed]
  3. Frieden, T.R.; Jaffe, H.W.; Cono, J.; Richards, C.L.; Iademarco, M.F. CDC Guideline for Prescribing Opioids for Chronic Pain—United States, 2016, MMWR. U.S. Department of Health and Human Services. 2016. Available online: https://www.cdc.gov/mmwr/volumes/65/rr/rr6501e1.htm (accessed on 25 August 2022).
  4. Solhaug, V.; Molden, E. Individual variability in clinical effect and tolerability of opioid analgesics—Importance of drug interactions and pharmacogenetics. Scand. J. Pain 2017, 17, 193–200. [Google Scholar] [CrossRef]
  5. Mura, E.; Govoni, S.; Racchi, M.; Carossa, V.; Ranzani, G.N.; Allegri, M.; Schaik, V.; Govoni, S.; Racchi, M.; Carossa, V.; et al. Consequences of the 118A>G polymorphism in the OPRM1 gene: Translation from bench to bedside? J. Pain Res. 2013, 6, 331–353. [Google Scholar] [CrossRef]
  6. Coller, J.K.; Beardsley, J.; Bignold, J.; Li, Y.; Merg, F.; Sullivan, T.; Cox, T.C.; Somogyi, A.A. Lack of association between the A118G polymorphism of the mu opioid receptor gene (OPRM1) and opioid dependence: A meta-analysis. Pharmacogenom. Pers. Med. 2009, 2, 9–19. [Google Scholar] [CrossRef]
  7. Escorial, M.; Muriel, J.; Agulló, L.; Zandonai, T.; Margarit, C.; Morales, D.; Peiró, A.M. Clinical prediction of opioid use disorder in chronic pain patients: A cohort-retrospective study with a pharmacogenetic approach. Minerva Anestesiol. 2024, 90, 386–396. [Google Scholar] [CrossRef] [PubMed]
  8. Muriel, J.; Margarit, C.; Planelles, B.; Serralta, M.J.; Puga, C.; Inda, M.; Cutillas, E.; Morales, D.; Horga, J.F.; Peiró, A.M. OPRM1 influence on and effectiveness of an individualized treatment plan for prescription opioid use disorder patients. Ann. N. Y. Acad. Sci. 2018, 1425, 82–93. [Google Scholar] [CrossRef]
  9. Muriel, J.; Margarit, C.; Barrachina, J.; Ballester, P.; Flor, A.; Morales, D.; Horga, J.F.; Fernández, E.; Peiró, A.M. Pharmacogenetics and prediction of adverse events in prescription opioid use disorder patients. Basic Clin. Pharmacol. Toxicol. 2019, 124, 439–448. [Google Scholar] [CrossRef]
  10. Planelles, B.; Margarit, C.; Ajo, R.; Sastre, Y.; Muriel, J.; Inda, M.; Esteban, M.D.; Peiró, A.M. Health benefits of an adverse events reporting system for chronic pain patients using long-term opioids. Acta Anaesthesiol. Scand. 2019, 63, 248–258. [Google Scholar] [CrossRef]
  11. Planelles, B.; Margarit, C.; Inda, M.-D.; Ballester, P.; Muriel, J.; Barrachina, J.; Ajo, R.; Esteban, M.-D.; Peiró, A.M. Gender based differences, pharmacogenetics and adverse events in chronic pain management. Pharmacogenom. J. 2020, 20, 320–328. [Google Scholar] [CrossRef]
  12. Boscarino, J.; Hoffman, S.N.; Han, J.J. Opioid-use disorder among patients on long-term opioid therapy: Impact of final DSM-5 diagnostic criteria on prevalence and correlates. Subst. Abus. Rehabil. 2015, 6, 83–91. [Google Scholar] [CrossRef] [PubMed]
  13. McCormack, H.M.; Horne, D.J.; Sheather, S. Clinical applications of visual analogue scales: A critical review. Psychol. Med. 1988, 18, 1007–1019. [Google Scholar] [CrossRef] [PubMed]
  14. Barrachina, J.; Muriel, J.; Margarit, C.; Planelles, B.; Ballester, P.; Richart-Martínez, M.; Cutillas, E.; Zandonai, T.; Morales, D.; Peiró, A.M. Global Pain State Questionnaire: Reliability, Validity, and Gender Gap. Arch. Intern. Med. Res. 2021, 04, 91–113. [Google Scholar] [CrossRef]
  15. Herdman, M.; Badia, X.; Berra, S. El EuroQol-5D: Una alternativa sencilla para la medición de la calidad de vida relacionada con la salud en atención primaria. Aten. Prim. 2001, 28, 425–429. [Google Scholar] [CrossRef] [PubMed]
  16. Pergolizzi, J.; Böger, R.H.; Budd, K.; Dahan, A.; Erdine, S.; Hans, G.; Kress, H.; Langford, R.; Likar, R.; Raffa, R.B.; et al. Opioids and the Management of Chronic Severe Pain in the Elderly: Consensus Statement of an International Expert Panel with Focus on the Six Clinically Most Often Used World Health Organization step III Opioids (Buprenorphine, Fentanyl, Hydromorphone, Methadone, Morphine, Oxycodone). Pain Pract. 2008, 8, 287–313. [Google Scholar] [CrossRef]
  17. Moore, R.A.; McQuay, H.J. Prevalence of opioid adverse events in chronic non-malignant pain: Systematic review of randomised trials of oral opioids. Arthritis Res. Ther. 2005, 7, R1046–R1051. [Google Scholar] [CrossRef]
  18. MedDRA®, Introductory Guide MedDRA Version 23.1. 2020. Available online: https://admin.meddra.org/sites/default/files/guidance/file/intguide_%2023_1_English.pdf (accessed on 3 September 2024).
  19. Wisher, D. Martindale: The Complete Drug Reference. 37th ed. J. Med. Libr. Assoc. 2012, 100, 75–76. [Google Scholar] [CrossRef]
  20. Steyerberg, E.W.; Vergouwe, Y. Towards better clinical prediction models: Seven steps for development and an ABCD for validation. Eur. Heart J. 2014, 35, 1925–1931. [Google Scholar] [CrossRef]
  21. Oliva, E.M.; Bowe, T.; Tavakoli, S.; Martins, S.; Lewis, E.T.; Paik, M.; Wiechers, I.; Henderson, P.; Harvey, M.; Avoundjian, T.; et al. Development and applications of the Veterans Health Administration’s Stratification Tool for Opioid Risk Mitigation (STORM) to improve opioid safety and prevent overdose and suicide. Psychol. Serv. 2017, 14, 34–49. [Google Scholar] [CrossRef]
  22. Crews, K.R.; Monte, A.A.; Huddart, R.; Caudle, K.E.; Kharasch, E.D.; Gaedigk, A.; Dunnenberger, H.M.; Leeder, J.S.; Callaghan, J.T.; Samer, C.F.; et al. Clinical Pharmacogenetics Implementation Consortium Guideline for CYP2D6, OPRM1, and COMT Genotypes and Select Opioid Therapy. Clin. Pharmacol. Ther. 2021, 110, 888–896. [Google Scholar] [CrossRef]
  23. Burns, J.A.; Kroll, D.S.; Feldman, D.E.; Liu, C.K.; Manza, P.; Wiers, C.E.; Volkow, N.D.; Wang, G.-J. Molecular Imaging of Opioid and Dopamine Systems: Insights into the Pharmacogenetics of Opioid Use Disorders. Front. Psychiatry 2019, 10, 626. [Google Scholar] [CrossRef]
  24. Ettienne, E.B.; Chapman, E.; Maneno, M.; Ofoegbu, A.; Wilson, B.; Settles-Reaves, B.; Clarke, M.; Dunston, G.; Rosenblatt, K. Pharmacogenomics-guided policy in opioid use disorder (OUD) management: An ethnically-diverse case-based approach. Addict. Behav. Rep. 2017, 6, 8–14. [Google Scholar] [CrossRef]
  25. Thielke, S.M.M.; Turner, J.A.; Shortreed, S.M.; Saunders, K.J.; LeResche, L.S.; Campbell, C.I.; Weisner, C.C.D.; Von Korff, M. Do Patient-perceived Pros and Cons of Opioids Predict Sustained Higher-Dose Use? Clin. J. Pain 2014, 30, 93–101. [Google Scholar] [CrossRef]
  26. Johnson, K.B.; Wei, W.; Weeraratne, D.; Frisse, M.E.; Misulis, K.; Rhee, K.; Zhao, J.; Snowdon, J.L. Precision Medicine, AI, and the Future of Personalized Health Care. Clin. Transl. Sci. 2021, 14, 86–93. [Google Scholar] [CrossRef]
  27. Lo-Ciganic, W.-H.; Huang, J.L.; Zhang, H.H.; Weiss, J.C.; Wu, Y.; Kwoh, C.K.; Donohue, J.M.; Cochran, G.; Gordon, A.J.; Malone, D.C.; et al. Evaluation of Machine-Learning Algorithms for Predicting Opioid Overdose Risk Among Medicare Beneficiaries with Opioid Prescriptions. JAMA Netw. Open 2019, 2, e190968. [Google Scholar] [CrossRef]
  28. Liu, Y.S.; Kiyang, L.; Hayward, J.; Zhang, Y.; Metes, D.; Wang, M.; Svenson, L.W.; Talarico, F.; Chue, P.; Li, X.-M.; et al. Individualized Prospective Prediction of Opioid Use Disorder. Can. J. Psychiatry 2023, 68, 54–63. [Google Scholar] [CrossRef]
  29. Armoon, B.; Fleury, M.-J.; Bayat, A.-H.; Bayani, A.; Mohammadi, R.; Griffiths, M.D. Quality of life and its correlated factors among patients with substance use disorders: A systematic review and meta-analysis. Arch. Public Health 2022, 80, 179. [Google Scholar] [CrossRef]
  30. Sanger, N.; Bhatt, M.; Shams, I.; Shahid, H.; Luo, C.; Tam, S.L.; Samaan, M.C.; de Souza, R.; Thabane, L.; Samaan, Z. Association Between Socio-Demographic and Health Functioning Variables Among Patients with Opioid Use Disorder Introduced by Prescription: A Prospective Cohort Study. Pain Physician 2018, 21, E623–E632. [Google Scholar]
  31. Smith, D.M.; Figg, W.D. Evidence Regarding Pharmacogenetics in Pain Management and Cancer. Oncologist 2023, 28, 189–192. [Google Scholar] [CrossRef]
  32. Morris, S.A.; Alsaidi, A.T.; Verbyla, A.; Cruz, A.; Macfarlane, C.; Bauer, J.; Patel, J.N. Cost Effectiveness of Pharmacogenetic Testing for Drugs with Clinical Pharmacogenetics Implementation Consortium (CPIC) Guidelines: A Systematic Review. Clin. Pharmacol. Ther. 2022, 112, 1318–1328. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the patients included in a real-world Pain Unit setting.
Figure 1. Flow chart of the patients included in a real-world Pain Unit setting.
Biomedicines 12 02056 g001
Table 1. Sociodemographic, clinical and pharmacological characteristics of Samples 1 and 2.
Table 1. Sociodemographic, clinical and pharmacological characteristics of Samples 1 and 2.
Retrospective
Sample 1
(n = 129)
Prospective
Sample 2
(n = 100)
Sex (% female)6770
Age (years old) (median (IQR))63 (52–72)65 (52–73)
Employment status (%)
Active1513
Retired5040
Work disability2126
Unemployed76
Homemaker715
Previous SUD (%)1825
Tobacco96 *71
Alcohol 425
Illicit substances04
Incomes (%)
Less than EUR 500205
Between EUR 500 to 10006753
More than EUR 10001342 *
Clinical outcomes (mean (SD))
Pain intensity (VAS, mm) 61 (28)70 (26) *
Pain relief (VAS, mm) 38 (31)41 (31)
Quality of life (VAS, mm)46 (24)46 (28)
Health Utility (0–1 score) (median (IQR))0.514 (0.113–0.732)0.252 (0.051–0.648)
Health resource use (%)
Emergency room visits3042
Hospitalizations1425
Medication changes5051
Drug prescription (%)
Non-opioid analgesics3463 *
NSAIDs1722
Tramadol2245 *
MEDD (mg/day) (median (IQR))80 (40–160) *60 (33–108)
Oxycodone67 *14
Fentanyl1524
Tapentadol1137 *
Buprenorphine322 *
Morphine33
Hydromorphone10
Immediate release opioids24 *10
Neuromodulators5260
Antidepressants5046
Benzodiazepines3554 *
NSAIDs: non-steroidal anti-inflammatory drugs; MEDD: morphine equivalent daily dose, VAS: visual analog scale, SUD: substance use disorder. * p-value < 0.05 comparing Sample 1 vs. Sample 2.
Table 2. The logistic regression model chosen to classify patients in Sample 1 or 2.
Table 2. The logistic regression model chosen to classify patients in Sample 1 or 2.
β-Coefficients95% CIStd. Errorz-ValuePr (>|z|) a
Intercept0.242−1.53 to 1.950.8730.2770.78
Active−1.950−3.65 to −0.350.831−2.3470.02
Work disability−1.740−3.05 to −0.540.633−2.7500.006
Unemployed−3.976−5.71 to −2.510.809−4.914<0.001
MEDD0.004−0.00 to 0.010.0031.5150.13
Strong opioids3.4932.05 to 5.240.8034.349<0.001
Benzodiazepines−2.626−3.94 to −1.500.617−4.254<0.001
ED visits−1.496−2.64 to −0.450.552−2.7070.007
Psychiatric AEs2.2891.21 to 3.520.5833.929<0.001
COMTGA−2.159−3.60 to −0.890.686−3.1470.002
GG−0.901−2.56 to 0.700.824−1.0940.27
AEs, adverse events, ED: emergency department, MEDD: morphine equivalent daily dose, a p-value associated with the z-value.
Table 3. Independent opioid use disorder risk predictors selected in Sample 1.
Table 3. Independent opioid use disorder risk predictors selected in Sample 1.
β-Coefficients95% CIStd. Errorz-ValuePr (>|z|) a
Intercept−0.622−4.77 to 2.991.933−0.3220.75
Age−0.057−0.12 to 0.000.032−1.7980.07
Work disability2.8601.33 to 4.780.8483.373<0.001
MEDD0.0060.00 to 0.010.0032.4440.02
CYP2D6 PM1.191−2.19 to 3.901.4420.8260.41
UM3.2990.82 to 5.971.2552.6280.009
MEDD: morphine equivalent daily dose, PM: poor metabolizer; UM: ultra-rapid metabolizer. a p-value associated with the z-value.
Table 4. Independent opioid use disorder risk predictors selected in Sample 2.
Table 4. Independent opioid use disorder risk predictors selected in Sample 2.
β-Coefficients95% CIStd. Errorz-ValuePr (>|z|) a
Intercept−1.713−3.35 to −0.310.759−2.2560.02
Quality of life−0.032−0.06 to −0.010.014−2.3020.02
MEDD0.005−0.00 to 0.010.0041.3940.16
OPRM1 (AG/GG)1.017−0.36 to 2.560.7271.4000.16
MEDD: morphine equivalent daily dose, a p-value associated with the z-value.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Escorial, M.; Muriel, J.; Margarit, C.; Agulló, L.; Zandonai, T.; Panadero, A.; Morales, D.; Peiró, A.M. Predictive Model for Opioid Use Disorder in Chronic Pain: A Development and Validation Study. Biomedicines 2024, 12, 2056. https://doi.org/10.3390/biomedicines12092056

AMA Style

Escorial M, Muriel J, Margarit C, Agulló L, Zandonai T, Panadero A, Morales D, Peiró AM. Predictive Model for Opioid Use Disorder in Chronic Pain: A Development and Validation Study. Biomedicines. 2024; 12(9):2056. https://doi.org/10.3390/biomedicines12092056

Chicago/Turabian Style

Escorial, Mónica, Javier Muriel, César Margarit, Laura Agulló, Thomas Zandonai, Ana Panadero, Domingo Morales, and Ana M. Peiró. 2024. "Predictive Model for Opioid Use Disorder in Chronic Pain: A Development and Validation Study" Biomedicines 12, no. 9: 2056. https://doi.org/10.3390/biomedicines12092056

APA Style

Escorial, M., Muriel, J., Margarit, C., Agulló, L., Zandonai, T., Panadero, A., Morales, D., & Peiró, A. M. (2024). Predictive Model for Opioid Use Disorder in Chronic Pain: A Development and Validation Study. Biomedicines, 12(9), 2056. https://doi.org/10.3390/biomedicines12092056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop