key: cord-0017518-lfv3mveq
authors: Ratigan, Amanda Rondinelli; Michaud, Veronique; Turgeon, Jacques; Bikmetov, Ravil; Gaona Villarreal, Gabriela; Anderson, Heather D.; Pulver, Gerald; Pace, Wilson D.
title: Longitudinal Association of a Medication Risk Score With Mortality Among Ambulatory Patients Acquired Through Electronic Health Record Data
date: 2021-03-31
journal: J Patient Saf
DOI: 10.1097/pts.0000000000000829
sha: 3324021ea2a1cde3799870026a03dcd8afaf9d0c
doc_id: 17518
cord_uid: lfv3mveq

The use of electronic health records allows for the application of a novel medication risk score for the rapid identification of ambulatory patients at risk of adverse drug events. We sought to examine the longitudinal association of medication risk score with mortality. This retrospective cohort study included patients whose data were available through electronic health records from multiple health care organizations in the United States that provided data as part of a Patient Safety Organization. Patients were included if they had ≥1 visit and ≥1 medication in their record between January 1, 2011, to June 30, 2017. Cox proportional hazards regression was used to examine the association between continuous and categorized medication risk score with all-cause mortality. Among 427,103 patients, the median age was 50 years (interquartile range, 29–64 years); 61% were female; 50% were White, 11% were Black, and 38% were Hispanic; and 6873 had a death date recorded. Patients 30 to 49 years old had the highest hazard ratios (HRs), followed by the 50- to 64-year-olds and lastly those 65 years or older. Controlling for all covariates, 30- to 49-year-olds with a score of 20 to 30 (versus <10) had a 604% increase in the hazard of death (HR, 7.04; 95% confidence interval [CI], 3.86–12.85), 50- to 64-year-olds had a 254% increase (HR, 3.54; 95% CI, 2.71–4.63), and ≥65-year-olds had an 87% increase (HR, 1.87; 95% CI, 1.67–2.09). The medication risk score was independently associated with death, adjusting for multimorbidities and other conditions. Risk was found to vary by age group and score. Results suggest that pharmaceutical interventions among those with elevated scores could improve medication safety for patients taking multiple medications.

T he act of prescribing a medication is the most common health care intervention in the United States, with approximately 19% of Americans aged 45 to 64 years and almost 40% of adults 65 years and older taking 5 or more prescription drugs in 2015 to 2016. [1] [2] [3] [4] Moreover, prescribing, dispensing, and ingestion of a medication by a patient have long been recognized as the riskiest ambulatory activity in terms of absolute number of people harmed. 5, 6 This results both from the high number of medications prescribed 3 and the complex processes involved between clinical decision making, communication to the pharmacy, dispensing, and actual use. 1, 2 The widespread adoption of electronic health records (EHRs) has been considered the primary driver to more adequately address medication risk. [7] [8] [9] Electronic health records have improved legibility issues and perhaps dosing and diseasedrug and drug interactions [10] [11] [12] [13] ; however, medication safety, and the prevention of medication errors and subsequent adverse drug events (ADEs) continues to be a significant global public health challenge. [14] [15] [16] [17] In the United States, it is estimated that ADEs account for approximately 3.5 million physician office visits, 1 million emergency department (ED) visits, and 125,000 hospital admissions annually. 18, 19 Prevention of ADEs is a key patient safety priority of many national health care-regulating agencies, [20] [21] [22] and practical, evidence-based tools for the early identification of patients at risk of ADEs are needed to significantly improve medication safety. 23 The increase in the average number of medications used by a single individual has highlighted the need to develop new methods of addressing medication risk. 3, 24 As process errors in the order, dispensing, and use chain improve, more complex drug and disease interactions can be prioritized to reduce harm from therapeutic medication use. [25] [26] [27] Multimorbidity is associated with higher complexity of medication regimens and a greater chance that complex interactions between medications and an individual's conditions may lead to harm. 28, 29 In an attempt to manage medication safety, many larger health care systems have implemented Medication Therapy Management (MTM) programs, and the U.S. Centers for Medicare & Medicaid Services requires that MTM programs be included as part of the provision of Part D benefits. [30] [31] [32] These programs typically identify individuals based on the number of drugs they are using, particular disease states (e.g., diabetes mellitus or hypertension) and/or an event (e.g., hospitalization). [33] [34] [35] A clinical pharmacist is then needed to review the medication regimens and to talk with the patients to determine if adjustments are warranted. This process relies heavily on the pharmacists' clinical knowledge of the various conditions encountered in each patient. More advanced techniques are needed as knowledge increases regarding medication metabolism, side effects, and interactions with developed and fixed (i.e., genetic) patient characteristics. In 2017, the U.S. Centers for Medicare & Medicaid Services introduced a Part D Enhanced MTM 5-year pilot program where advanced clinical decision support systems and pharmacist-directed interventions are used. 36 Computer-based risk algorithms hold the potential to efficiently and effectively support medication safety programs. 23, 37 Tabula Rasa HealthCare has developed a proprietary advanced clinical decision support system and medication risk score (MRS), which incorporates a series of algorithms that calculate an overall score, as well as subcomponent scores on sedation burden, [38] [39] [40] anticholinergic cognitive burden, 41 a competitive CYP450 drug interaction burden, 24 drug regimen relative odds ratio for ADE using the Food and Drug Administration adverse event reporting system, 42 and a drug-induced risk of creating or worsening a prolonged QT cardiac repolarization interval. 37, 42, 43 An observational study using claims data among participants of the Program of All-Inclusive Care for the Elderly found the total MRS to be positively correlated with increased ADEs, health care utilization, hospitalization, ED visits, and hospital length of stay. 44 However, this study did not assess an independent risk of death and did not include other clinical variables that may have a strong impact on the overall score.

The current observational study was conducted as part of the safety work of the DARTNet Institute subcomponent Patient Safety Organization (PSO). 45 We hypothesized that the MRS was an independent predictor of death among a cohort of ambulatory patients obtained through EHR data. In comparison to claims data, EHRs include clinical variables (e.g., estimated glomerular filtration rate [eGFR], blood pressure [BP] values, and laboratory measured electrolyte levels), allowing EHR-based algorithms to have the potential for greater dissemination and more rapid response to medication changes, which may improve health outcomes and reduce health care costs.

This retrospective cohort study uses EHR data from multiple U.S. health care organizations that provide data in partnership with DARTNet's PSO. This study includes community-living, ambulatory patients 5 years and older who had data present between January 1, 2011, and June 30, 2017, and who had at least one medication present (i.e., prescribed or added during medication reconciliation) within the study period.

The dependent variable of interest is death from any cause as determined by the presence of a date of death within the patient's EHR record. The main independent variable of interest is MRS, which ranges from 0 to 53, with a higher score indicating a higher risk of ADEs due to a particular medication regime, including anticholinergic burden, sedating effects, risk of QT-interval prolongation, and the competitive inhibition on certain cytochrome P450 isoforms. 42, 46 Medication risk score was also categorized for analysis (<10, 11-14, 15-19, 20-30, >30) . Details regarding the calculation of the MRS and its components have been published elsewhere. 37, 42, 46 The calculation of the MRS considers all prescribed medications and patients' characteristics during a given period of time. Data were extracted from EHR records and compiled into eras (i.e., time frames) to account for the dynamic nature of the MRS. An era is defined by current medication exposures (including strength); laboratory results categorized as low, normal, or high (calcium (<8.5, 8.5-10.5, and >10.5 mg/dL), potassium (<3.5, 3.5-5, and >5 mg/dL), magnesium (<1.5, 1.5-2.5, and >2.5 mg/dL); and selected conditions (i.e., end-stage renal disease, atrial fibrillation, end-stage liver disease, atherosclerotic cardiovascular disease, heart failure, cardiomyopathy, prolonged QT interval, sick sinus syndrome). A change in any of these factors closes the current era and commences a new era, as long as there continues to be at least one medication exposure.

Each patient-era includes covariates measured during that era including the following: sociodemographic factors (i.e., age, sex, and race/ethnicity as identified in the EHR, categorized as White, non-Hispanic or unknown ethnicity; Black, non-Hispanic or unknown ethnicity; Hispanic/Latino; other/mixed non-Hispanic; and unknown/missing both race and ethnicity), health-related factors (body mass index [BMI], BP, kidney function based on eGFR, chronic kidney disease, and liver disease), and total days in the era and total medications listed in the EHR during that particular era. Age was divided into quartiles (ages <30, 30-49, 50-64, and ≥65 years); systolic and diastolic BP values were limited to clinically plausible ranges of 50 to 250 and 0 to 200 mm Hg, respectively, and were categorized as normal (systolic BP <120 mm Hg or diastolic BP <80 mm Hg) or elevated (systolic BP ≥140 mm Hg or diastolic BP ≥90 mm Hg); BMI values were limited to 12 to 70 kg/m 2 and were categorized for adults into underweight (<18.5 kg/m 2 ), normal (18.5-24 kg/m 2 ), overweight (25-30 kg/m 2 ), or obese (>30 kg/m 2 ); eGFR (in mL min −1 1.73 m −2 ) was categorized into stages of chronic kidney disease (stages 1-2, ≥60; stage 3a, 45-59; stage 3b, 30-44; stages 4-5, <30); and liver disease was categorized as no liver disease, chronic liver disease, or end-stage liver disease based on International Classification of Diseases, Tenth Revision diagnoses. All outlier values were set to missing for analyses. Finally, the Charlson-Deyo Score, calculated from the sum of scores for each comorbid condition listed in the Charlson Comorbidity Score Mapping Table, was used to account for multimorbidity. 47, 48 Statistical Analysis Descriptive statistics were calculated using mean and SD, or median and interquartile range (IQR) for continuous measures; and frequency and percentage for categorical measures. Differences by mortality status were examined using χ 2 tests for categorical measures and t tests for continuous measures.

Cox proportional hazards regression modeling with time-varying covariates was used to examine the unadjusted and multivariable association of MRS with all-cause mortality and to account for within-person changes over time and eras that may occur among all variables. Proportional hazards assumptions were verified by examining Schoenfeld residuals and were determined to have been met. 49 The primary outcome was number of days from the start date of an era to the date of death if it occurred during a particular era. Patients were censored at the end of each era if they did not have a date of death ending that era. Covariates were adjusted for using a stepwise approach to determine their effect on the association of MRS with death, where variables were considered to be a confounder of this association if they changed the point estimate of the main association >10% when included into the multivariable model. Hazard ratios (HRs) with 95% confidence intervals (CIs) are reported to show strength and direction of these associations. Data were analyzed using SAS 9.4 (Cary, North Carolina), and a level of P < 0.05 was used to determine statistical significance.

A total of 427,103 patients with at least one visit to a qualifying medical practice and at least one medication record during the study period were eligible for MRS calculation. A total of 6873 patients (1.6%) died by the end of the study period. The median age at baseline was 50 years (IQR, 29-64 years), and 25% were younger than 30 years (with 14.7% being 5-18 years old), 24% were 30 to 49 years old, 26% were 50 to 64 years old, and 24.5% were 65 years and older (Table 1) . More than half were female (61%); almost 50% were White, 11% were Black, and almost 38% were Hispanic. Over the entire study period, there were a total of 2,491,399 eras present for analysis, with a median number of eras per patient of 12 (IQR, 5-25; range, 1-236) and a median number of days per era of 62 (IQR, 14-189; range, 1-2372).

Compared with patients who did not have a date of death present within their EHR, those who did were significantly more likely to be older at baseline (74 versus 49 years; P < 0.01), male (52.3% versus 38.7%; P < 0.01), and White (57.2% versus 49.4%; P < 0.01) (Table 1) . Over the entire study period, they had a significantly higher median MRS (9 versus 5; P < 0.01), higher median number of medications (7 versus 4; P < 0.01), lower average BMI (28.7 versus 29.4 kg/m 2 ; P < 0.01), lower average eGFR (58.2 versus 78.1 mL min −1 1.73 m −2 ; P < 0.01), and higher average Charlson-Deyo score (4 versus 2; P < 0.01), and were more likely to have an elevated BP (26.8% versus 20.3%; P < 0.01). Furthermore, those who died had a significantly greater median number of total eras (22 versus 12; P < 0.01), with lower median days per era (36 versus 63; P < 0.01) and significantly greater average number of days in the study period (1654 versus 1530; P < 0.01). Table 2 displays the results of age-stratified unadjusted Cox proportional hazards regression models examining the association 

In this large cohort of ambulatory patients in the United States, the MRS was demonstrated to be independently associated with death after adjustments for common multimorbidities and other conditions that increase the risk of death. Furthermore, it was found that risk varies by age group and by MRS category where the highest risk was seen among patients in the 30-to 49-yearold age group who had an MRS of ≥20. Our results suggest that pharmacist-directed interventions targeted to those who have an elevated MRS could improve medication safety for ambulatory patients taking multiple medications.

The MRS encompasses a series of algorithms that consider the pharmacological characteristics of each drug's active ingredients and estimates whether the combination of drugs exhibits properties that, when taken together, put a patient at risk of ADEs. Hence, in patients with multiple chronic diseases, drugs used for each of these diseases may seem appropriate, yet their combination may predispose them to ADEs. Advanced clinical decision support systems could inform clinicians on the appropriateness of a drug regimen, indicate what elements could be improved, and offer alternatives. 46 Results from this study expand the evidence that the MRS is correlated with increased rates of ADEs and health care utilization in an elderly population 44 by demonstrating that MRS is also significantly associated with death among a general population using EHR data. Although the use of EHR data includes clinical variables, it may also overestimate the actual medication use at the individual patient level because it includes medications prescribed rather than medications filled 50 ; however, the number of medications was not found to have an effect on the association between MRS and death in this population.

Men and patients categorized as White race were at higher risk of death than women and other races/ethnicities. The reason for this is not intuitively obvious given the range of variables controlled for in this analysis. Men overall have higher cardiovascular death rates, of which arrhythmias are a contributing factor, 51 which could partially explain the higher male versus female risk; however, lower minority rates, where males also have high cardiovascular death rates, 51 are not explained by this mechanism. Furthermore, the risk varied by age group where the highest hazards were seen among those aged 30 to 49 years and the lowest among those 65 years or older, indicating that medication risk is playing a larger part in all-cause mortality among younger age groups than older age groups. Nonetheless, even the ≥65-year older age range had significantly elevated HRs. Thus, at the population level, potential medication-related deaths may be higher among older patients than among 30-to 49-year-olds given the larger number of people with high MRSs. The low number of patients 30 years and younger with high MRSs and the decreased number of deaths in this age group prohibited the estimation of HRs. Larger data sets of younger patients will be required to understand the correlation of high MRS with death in this age group.

Risk of death increased with higher MRS. Confidence intervals widened considerably at the highest risk scores, as the number of patients and therefore the number of deaths were small in these cohorts. Larger data sets with even greater multimorbidity and thus complex medication regimens would be needed to narrow the CIs in these cohorts. Nonetheless, the lower bounds of the risk of death do not approach one for any age group that had sufficient data for analysis.

Despite several strengths of this study, some limitations should be considered when interpreting these results. First, EHR data were used, which relies on data hand-entered by clinicians and documented through laboratory interfaces and medication reconciliation processes. It is unknown if the patients were actually filling and consuming the medications they were prescribed, and it is clear from other research that EHR data do not provide accurate medication use data as opposed to claims data. 50 Furthermore, EHR data provides information on clinician intent (prescribing a medication) and not patient actions (such as filling a medication). Nonetheless, our results suggest that EHR data seem adequate to predict poorer outcomes based on MRS. Electronic health record data are readily available at the clinical organization level, which may allow for more rapid recognition and intervention based on high-risk medication regimens.

It would be expected that the combination of EHR and claims data would further enhance the predictive ability of the MRS, but the finding that EHR data alone provides a reasonable estimate of possible risk allows for easier use of the risk scores and facilitates future research. Death was determined using the presence of a date of death within the EHR relying on the provider's office to enter this date, which could have led to missing death data. It was not possible to reach out to state death records to enhance death data, as many organizations strip out identifying information before submitting data to the PSO. Likewise, as far as the authors are aware, none of the organizations in the DARTNet PSO routinely collect data from their state death registries. Nonetheless, there is no reason to suspect that deaths would be differentially recorded in relationship to an MRS that was not known to the sites and was not calculated until after the deaths were recorded.

This study provides further evidence that the studied MRS is correlated with poor health outcomes across a wide range of ages. Further research is warranted to demonstrate that lowering risk scores in individuals, through changes in timing of medication use or changes in prescribed medications, results in reductions in ADEs and death. If this is demonstrated, then the application of combined, computer-generated risk scores to identify and adjust medication regimens could significantly improve health care outcomes among ambulatory patients taking medications.

Eliminating medication overload: a national action plan

Medication overload: America's other drug problem

Solving Medication Overload: American's Other Drug Problem

Advances in patient safety relationship between patient harm and reported medical errors in primary care: a report from the ASIPS Collaborative

Advances in patient safety mixed methods analysis of medical error event reports: a report from the ASIPS Collaborative

Medication errors: prevention using information technology systems

Professional, structural and organisational interventions in primary care for reducing medication errors

Committee on Identifying and Preventing Medication Errors. Preventing Medication Errors: Quality Chasm Series

Effects of computer-based clinical decision support systems on physician performance and patient outcomes: a systematic review

Effects of computer-based clinical decision support systems on clinician performance and patient outcome. A critical appraisal of research

A computer alert system to prevent injury from adverse drug events: development and evaluation in a community teaching hospital

Improving safety with information technology

Use of Medicare administrative claims to identify a population at high risk for adverse drug events and hospital use for quality improvement

Incidence and preventability of adverse drug events among older persons in the ambulatory setting

National surveillance of emergency department visits for outpatient adverse drug events

Preventable adverse drug reactions: a focus on drug interactions

National Action Plan for Adverse Drug Event Prevention

Adverse drug events in the outpatient setting: an 11-year national analysis

The Joint Commission. Ambulatory Health Care: 2020 National Patient Safety Goals

Agency for Healthcare Research and Quality. Patient Safety Network. Updates in the Role of Health IT in Patient Safety

Clinical decision support systems: great promises for better management of patients' drug therapy

Prevalence and risk of potential cytochrome P450-mediated drug-drug interactions in older hospitalized patients with polypharmacy

The evidence for the effectiveness of safety alerts in electronic patient medication record systems at the point of pharmacy order entry: a systematic review

Dispensing error rate in a highly automated mail-service pharmacy practice

A long-term follow-up evaluation of electronic health record prescribing safety

The challenge of managing drug interactions in elderly people

The effects of polypharmacy in older adults

Medicare Part D Medication Therapy Management (MTM) Programs. Available at: https:// www.cms.gov/Medicare/Prescription-Drug-Coverage/ PrescriptionDrugCovContra/Downloads/CY2019-MTM-Fact-Sheet.pdf. Published

Medication therapy management interventions in outpatient settings: a systematic review and meta-analysis

Integrating comprehensive medication management to optimize patient outcomes

Post-discharge medication reconciliation: reduction in readmissions in a geriatric primary care clinic

Pharmacist's role in hypertension management: a review of key randomized controlled trials

Improving the adherence of type 2 diabetes mellitus patients with pharmacy care: a systematic review of randomized controlled trials

Part D Enhanced Medication Therapy Management Model

Treatment Methods Having Reduced Drug-Related Toxicity and Methods of Identifying the Likelihood of Patient Harm from Prescribed Medications

A comparison of four methods to quantify the cumulative effect of taking multiple drugs with sedative properties

A model to classify the sedative load of drugs

Sedative load among community-dwelling people aged 75 years and older: a population-based study

Impact of anticholinergics on the aging brain: a review and practical application

Application of a novel medication-related risk stratification strategy to a self-funded employer population

Risk assessment of drug-induced long QT syndrome for some COVID-19 repurposed drugs

Association of a novel medication risk score with adverse drug events and other pertinent outcomes among participants of the programs of all-inclusive care for the elderly

DARTNet Component Patient Safety Organization (DCPSO)

Population-based medication risk stratification and personalized medication risk score U.S. Patent WO2019089725

A new method of classifying prognostic comorbidity in longitudinal studies: development and validation

Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases

Partial residuals for the proportional hazards regression model

Standardizing terminology and definitions of medication adherence and persistence in research employing electronic databases

Heterogeneous trends in burden of heart disease mortality by subtypes in the United States, 1999-2018: observational analysis of vital statistics

We gratefully thank Dr Elizabeth Juarez-Colunga for her contribution of statistical knowledge and interpretation of these results.