Low Exposure to Direct Oral Anticoagulants Is Associated with Ischemic Stroke and Its Severity

Article information

J Stroke. 2022;24(1):88-97
Publication date (electronic) : 2022 January 31
doi : https://doi.org/10.5853/jos.2020.04952
aDepartment of Neurology, University of Heidelberg, Im Neuenheimer Feld, Heidelberg, Germany
bDepartment of Clinical Pharmacology and Pharmacoepidemiology, University of Heidelberg, Im Neuenheimer Feld, Heidelberg, Germany
cFaculty of Health/School of Medicine, Witten/Herdecke University, Witten, Germany
dDepartment for Child and Adolescent Psychiatry, Johannes Gutenberg-University, Mainz, Germany
Correspondence: Timolaos Rizos Department of Neurology, University Heidelberg, Im Neuenheimer Feld 400, 69120 Heidelberg, Germany Tel: +49-6221-56-7504 Fax: +49-6221-56-5654 E-mail: timolaos.rizos@med.uniheidelberg.de
*These authors contributed equally to the manuscript as first author.
Received 2020 December 13; Revised 2021 July 17; Accepted 2021 September 7.

Abstract

Background and purpose

In acute stroke patients, plasma concentrations of direct oral anticoagulants (DOAC) at hospital admission only poorly mirror DOAC exposure or the coagulation status at the time of the event. Here, we evaluated whether DOAC exposure and DOAC plasma concentration at the time of transient ischemic attacks (TIA) and ischemic strokes correlate with their likelihood of occurrence.

Methods

Prospectively, consecutive DOAC patients with acute ischemic stroke or TIA were included. Admission DOAC plasma concentrations were measured by ultraperformance liquid chromatography– tandem mass spectrometry. Individual DOAC exposure (area under the curve) and DOAC concentrations at event onset were derived from population pharmacokinetic analyses.

Results

DOAC exposure was successfully modeled in 211 patients (ischemic stroke 74.4%, TIA 25.6%). Compared to published values, 63.0% had relatively lower DOAC exposure and they more often received lower DOAC doses than recommended (odds ratio [OR], 2.125; 95% confidence interval [CI], 1.039 to 4.560; P=0.044). These patients more likely suffered ischemic strokes than TIA (OR, 2.411; 95% CI, 1.254 to 4.638; P=0.008) and their strokes were more severe (slope, 3.161; 95% CI, 0.741 to 5.58; P=0.011). Low relative DOAC concentrations at event onset were likewise associated with ischemic strokes (OR, 4.123; 95% CI, 1.834 to 9.268; P=0.001), but not to stroke severity (P=0.272). DOAC exposure had a higher explanatory value for stroke severity than concentrations at event.

Conclusions

Low DOAC exposure is strongly associated to ischemic stroke and its severity. By monitoring DOAC plasma concentrations, patients prone to ischemic stroke might be identified.

Introduction

Increasing numbers of patients exposed to direct oral anticoagulants (DOAC) are admitted to emergency departments with ischemic stroke and transient ischemic attack (TIA) [1-4]. Assessing the coagulation status and DOAC plasma concentrations at hospital admission is important in these patients because their results can influence acute treatment decisions [1-3]. Although non-specific coagulation tests are available, many centers still do not measure DOAC concentrations in these patients at all. Reportedly, specific coagulation tests are performed in only about one quarter of patients with ischemic stroke [2]. Therefore, current knowledge on DOAC plasma concentrations of patients in the acute phase of stroke is still limited, but available evidence suggests high exposure variability [2]. As an example, DOAC admission concentrations in the multicenter Registry of Acute Stroke Under New Oral Anticoagulants (RASUNOA) registry were below expected trough concentrations in 25% of patients [2]. Similar findings were observed in a single-center registry, which observed low plasma concentrations at admission in 27.7% of acute ischemic stroke patients [5] and in another registry-based study, only including rivaroxaban-treated acute ischemic stroke patients, DOAC plasma concentrations were low in 66.3% [6]. However, these reports cannot be readily compared because they applied different definitions for low DOAC plasma ranges, used different calibrator assays, and did not standardize the timing of the measurements with respect to last drug intake. Moreover, DOAC-specific coagulation tests used in emergency situations (most commonly chromogenic anti-Xa assays) can considerably differ in terms of diagnostic performance when compared to ultraperformance liquid chromatography–tandem mass spectrometry (UPLC-MS/MS), which is the gold standard for accurate DOAC quantification and which is much more accurate in the low concentration range [7].

Finally, although plasma concentrations at the time of hospital admission can guide treatment decisions, these measurements do not mirror DOAC exposure or the coagulation status at onset of an ischemic cerebrovascular event, in particular when hospital admission is considerably delayed. Knowledge of actual DOAC plasma concentrations at onset of an ischemic cerebrovascular event would enable evaluating and establishing a concentration-effect relationship and help deciding whether DOAC exposure of affected patients is within the target range (signifying non-response to treatment) or whether patients are underdosed or non-adherent.

In the present study, we addressed these questions by measuring DOAC plasma exposure at admission and estimating concentrations at the onset of ischemic symptoms.

Methods

Patients and clinical data

Consecutive patients aged >18 years, suspected to take a DOAC, and admitted to the neurological department of Heidelberg University Hospital due to symptoms of acute ischemic stroke or TIA between September 4, 2016 and June 12, 2018, were prospectively considered for inclusion into this observational study. Routine work-up for all patients encompassed a clinical examination by board-certified neurologists, brain imaging (computed tomography [CT] or magnetic resonance imaging [MRI]), and assessment of the National Institutes of Health Stroke Scale (NIHSS) for stroke severity. Occlusions of brain-supplying vessels were diagnosed by CT, MRI, or digital subtraction angiography. Occlusions of the distal internal carotid artery, the middle cerebral artery, the anterior cerebral artery, the posterior cerebral artery, or the basilar artery were defined as large vessel occlusion (LVO). Moreover, demographic variables, cardiovascular risk factors, time of symptom onset (if known), and time of hospital admission were recorded.

Nature and dose of DOAC and the time of the last DOAC intake were recorded. If the last intake was not reported or unknown (n=58/211; 27.5%), it was assumed that the DOAC was taken regularly and that the last ingestion occurred at 8:00 AM (once a day and twice a day regimens) and 8:00 PM (twice a day schedules). Physicians´ adherence to approved doses as recommended in the corresponding product characteristics (SmPCs) by the European Medicines Agency was evaluated by comparing the actual prescription scheme with the prescribing information of the respective DOAC in each patient. The degree of disability before stroke was assessed by the premorbid modified Rankin Scale (mRS). All patients with proven ischemic stroke or TIA, current DOAC prescription, and at least one immediate blood sampling at admission were considered for inclusion. Functional outcome after 3 months was assessed by the mRS during a standardized telephone interview.

Laboratory data and DOAC plasma concentrations

As part of our standard of care, immediately after admission, standardized routine laboratory testing was performed. This included a full blood count, glucose, electrolytes, urea, creatinine, prothrombin time, activated partial thromboplastin time, international normalized ratio (INR), and DOAC plasma concentrations. DOAC plasma concentrations were measured from plasma samples drawn at admission (t1) and 6 hours later (t2) by UPLC-MS/MS (lower limit of quantification [LLOQ] of 1 ng/mL), developed and validated according to U.S. Food and Drug Administration and European Medicines Agency guidelines for bioanalytical method validation [8]. DOAC exposure in the dosing interval and DOAC concentration at the onset of the ischemic event were derived from population pharmacokinetic analyses, based on established models for apixaban [9], dabigatran [10], edoxaban [11], and rivaroxaban [12]. Individual pharmacokinetic parameters were obtained as empirical Bayes estimates (conditional mode) from the respective model applied to each patient together with complete profiles (data points every 2 hours) of 50 simulated individuals with identical covariates as the respective patient. Of note, we used a simplified covariate-free adaptation of the complex edoxaban model [11] and generally set missing body weight measurements to the median value of the literature model in nine patients. For the simulated 50 virtual individuals serving as a reference, we assumed that steady-state was reached on the 5th treatment day.

Estimation of individual pharmacokinetic parameters of our patients enabled us to estimate individual concentration-time profiles (for an example, see Figure 1). Upon visual inspection of the profiles, we had to manually exclude 25 patients whose sparse data points did not yield clinically plausible profiles (Figure 2). Reasons for exclusion included (1) measurements or predicted profiles far beyond the 95% confidence range that were to expect from regular intake at steady-state (n=11); (2) extreme profiles that were not supported by any measurement within the regular dosing interval (n=15); and (3) implausible measurement pairs (e.g., second measurement with higher value than first measurement, n=6; please note: more than one reason was present in single patients). Based on the remaining 192 informative profiles, we calculated the area under the curve (AUC) to describe exposures with the respective DOAC and estimate DOAC concentration at onset (if the time of the event was known).

Figure 1.

Model-based concentration-time profile of an exemplary patient treated with apixaban (red line). Shaded areas (purple) visualize the typical range of 95% virtual patients as predicted from the underlying population pharmacokinetic model for apixaban while accounting for available covariates of the particular patient or assuming median values [9]. Red dots represent the actual measurements in the particular patient and the red line shows the a posteriori estimates of the individual concentration-time profile based on individual pharmacokinetic parameters of the particular patient. The concentration at the time of event (+black cross) was directly derived from the profile.

Figure 2.

Selection of the study population. DOAC, direct oral anticoagulant; TIA, transient ischemic attack; AUC, area under the curve; NA, non-adherence.

As a measure of internal validity, we checked the agreement between model-based extrapolated concentrations and non-parametric extrapolations in a sample of patients with more than one measurement. In order to estimate individual elimination rate constants (λ), two DOAC UPLC-MS/MS measurements had to be available that were collected ≥6 hours apart within the same dosing interval, of which the first had to be drawn not earlier than 2 hours after the last drug intake, i.e., after expected peak concentrations. The slope of the line connecting the logarithms of the two measured concentrations over time was defined as λ [λ=(log C1–log C2)/(t2–t1)], where C1 is the concentration in the first-drawn admission sample, C2 the concentration in the sample collected ≥6 hours after C1 within the same elimination phase, and t2–t1 the numeric interval between the corresponding concentrations.

This slope (l) was used to back-extrapolate concentrations (Ce) expected at the time (te) of the ischemic event [log Ce=λ((e–t1)+log C1], provided that symptom onset occurred after drug absorption (i.e., >2 hours after the last drug intake). In 19 patients all measurements were below the analytical LLOQ; these values were set to zero and the patients were considered non-adherent.

The responsible independent Ethics Committee of the Medical Faculty of Heidelberg University approved this study, written informed consent was obtained from patients or their legal representatives.

Statistical analysis

To facilitate standardized comparisons across the four DOACs, we calculated fold changes from dividing our estimated values by the mean concentration and AUC at steady-state extracted from external reference populations (for apixaban [13], dabigatran [14], edoxaban [15], and rivaroxaban [16]). Considering measurements set to zero in non-adherent patients, the Yeo-Johnson transformation yielded power-transformed values following a normal distribution around the center peak of 0.5 [17]. Thus, normalized fold changes above this normalized threshold indicated relatively higher concentrations or exposure than expected from the external reference (henceforth denoted as groups with relatively high, or vice versa low, relative DOAC AUCs or concentrations).

Standard statistical methods were applied for univariate comparisons. Multivariate logistic and linear regression was performed to explore factors explaining differences between patients with TIA and ischemic stroke, stroke severity (NIHSS), the occurrence of LVO in ischemic stroke patients, and the severity of clinical outcome (mRS) after 3 months (dichotomized at the level of 2).

To determine whether overall exposure (AUC) or the estimated concentration at the time of the event had a higher explanatory value to predict the nature of events (ischemic stroke or TIA) or stroke severity (NIHSS), we compared subpopulations for which both measures were available. Therefore, separate models were fitted including either normalized AUC or normalized concentration values as independent variables (in addition to clinical covariates for confounding adjustment). A formal model comparison of these non-nested models was based on an appropriate likelihood ratio test assessing the working hypothesis that the AUC model fits better than the respective concentration model [18].

All tests were two-sided and a P-value of ≤0.05 was considered significant. Data were analysed using the SPSS version 26.0 (IBM Co., Armonk, NY, USA), the R software/environment version 3.6.1 (R Foundation for Statistical Computing, Vienna, Austria), Monolix version 2019R2 (LIXOFT, Antony, France), and GraphPad Prism version 8.4.2 (GraphPad, La Jolla, CA, USA).

Results

Selection and characteristics of patients

The selection of the study population is illustrated in Figure 2. The median age of the 211 patients of our primary analysis set with evaluable DOAC exposure was 80 years (interquartile range [IQR], 75 to 86), 50.7% were female (Supplementary Table 1). Ischemic stroke was present in 74.4%, TIA in 25.6%. Etiology of ischemic stroke according to Trial of Org 10172 in Acute Stroke Treatment (TOAST) was cardioembolic in 79.0% (Supplementary Table 1) and in 33.6% LVO was observed. Most patients were either taking rivaroxaban (n=88; 41.7%) or apixaban (n=79; 37.4%); dabigatran was used by 24 (11.4%) and edoxaban by 20 (9.5%) (Supplementary Table 1). Expectedly, the predominant reason for oral anticoagulation was atrial fibrillation (AF; 91.0%) and the prevalence of cardiovascular risk factors was high, resulting in a median CHA2DS2-VASc Score prior to the actual event of 5 (IQR, 4 to 6).

DOAC plasma concentrations at admission and DOAC exposure

The median time between reported last drug administration and admission was 582 minutes (IQR, 282 to 891). DOAC concentrations at admission for individual DOACs are listed in Supplementary Table 1. The ratios of modeled exposure of individual patients and published average AUC values as a measure of individual DOAC exposure are presented in Figure 3 (residuals of underlying pharmacokinetic models are provided in Supplementary Figure 1). Non-adherence (UPLC-MS/MS measurements below LLOQ) was observed in 19 patients (9%).

Figure 3.

Box plots of the ratio of modeled area under the curve (AUC) and published figures corresponding to the respective doses with superimposed individual AUC values of all included patients. The broken line indicates the threshold between relatively higher and lower direct oral anticoagulant (DOAC) exposure than the published expected average. Each boxplot contains the median (horizontal line in the box), the upper quartile (75th percentile, top of box), the lower quartile (25th percentile, bottom of box). The whiskers plot the minimal and maximal DOAC exposure. Solid circles: Dot plots of categorized individual AUC results of all included patients (lower limit of quantification for ultraperformance liquid chromatography–tandem mass spectrometry 1 ng/mL). Open circles: non-adherent patients (concentration <1 ng/mL; n=19).

Lower relative DOAC exposure (AUC) than expected was observed in almost two-thirds (63.0%) of our patients. Clinical and DOAC-specific variables in patients with relatively lower and higher DOAC exposure are given in Supplementary Table 2 as are the results of univariate correlation analyses.

Multivariate logistic regression revealed a substantially higher risk for ischemic stroke when DOAC exposure was low (odds ratio [OR], 2.411; 95% confidence interval [CI], 1.254 to 4.638; P=0.008) (Table 1). Adjusted for sex, age, and creatinine, low DOAC exposure was associated with physicians’ non-adherence (underdosage) to SmPCs (OR, 2.125; 95% CI, 1.039 to 4.560; P=0.044) (Supplementary Table 3) and, in these ischemic stroke patients, low DOAC exposure was associated with greater stroke severity (NIHSS; P=0.011; multivariate linear regression) (Table 2), but not with the presence of LVO (OR, 1.360; 95% CI, 0.660 to 2.803; P=0.404) (Supplementary Table 4).

Multivariate logistic regression between ischemic stroke and TIA patients in all included patients and in all patients with modeled DOAC concentration at the event

Multivariate linear regression of stroke severity (NIHSS) in all included ischemic stroke patients and in all ischemic stroke patients with modeled DOAC concentration at the event

Estimated DOAC concentrations at the time of event

Information on the specific time of ischemic stroke or TIA was available in 179/211 patients thus allowing to extrapolate DOAC concentration at the time of the event (Figure 2). The median time between last drug intake and event was 480 minutes (IQR, 210 to 690) and neither differed between individual substances (PANOVA =0.173) nor between substances with once (rivaroxaban and edoxaban) or twice daily use (apixaban and dabigatran, P=0.221). Moreover, non-parametric comparison of the intervals between the time of the last drug intake and time of the event revealed no difference between TIA and ischemic stroke patients (PMann-Whitney-U test=0.832).

Univariate comparisons between patients with relatively lower and higher DOAC concentrations are provided in Supplementary Table 5. Multivariate logistic regression analysis revealed that patient with low relative onset concentrations were also more likely to suffer ischemic stroke than TIA (OR, 4.123; 95% CI, 1.834 to 9.268; P=0.001) (Table 1), but no association with stroke severity (NIHSS) was present (P=0.222) (Table 2). In ischemic stroke, lower DOAC concentrations at event onset were not associated with the presence of LVO (OR, 0.619; 95% CI, 0.287 to 1.335; P=0.222) (Supplementary Table 4).

Sensitivity analyses and substance-specific influences

Sensitivity analyses with alternating event definitions revealed robust results for the influence of DOAC exposure and DOAC plasma concentration at the time of the event on the risk of ischemic stroke also when data were restricted to patients with AF or patients with cardioembolic ischemic strokes (Supplementary Figure 2). No substance-specific modulation of the overall DOAC effect of exposure (AUC) or event concentration was present (Supplementary Figure 3). Hence, substance-specific influences did not appear to be present in our cohort.

Significance of DOAC exposure and DOAC concentration at event

In ischemic stroke patients, DOAC exposure (AUC) better predicted the stroke severity than estimated DOAC concentration at the event (P=0.049). However, AUC did not better explain the nature of observed ischemic events than estimated DOAC concentration (P=0.939).

Functional outcome

Overall, in 199 patients of our cohort (94.3%), follow-up information after 3 months was available. In these patients, median functional outcome (mRS) was 3 (1 to 5) (Supplementary Table 1). Of all included patients with ischemic stroke (n=157), follow-up information was available in 149 patients (94.9%). An unfavorable outcome (mRS 3–6) was observed in 92 of them (61.7%). Differences between ischemic stroke patients with favorable (mRS 0–2) and unfavorable (mRS 3–6) outcome are presented in Supplementary Table 6. In univariate analyses, neither DOAC exposure nor estimated concentration at the time of stroke were associated with the 3-month outcome (Supplementary Table 6). Multivariate logistic regression revealed that only stroke severity (OR, 1.174; 95% CI, 1.079 to 1.278; P<0.001) and the premorbid functional status (OR, 2.881; 95% CI, 1.720 to 4.828; P<0.001) were independently associated with an unfavorable outcome after 3 months (Supplementary Table 7).

Discussion

The main findings of our study are that (1) a large proportion of patients hospitalized with acute ischemic cerebrovascular events had lower than expected DOAC exposure; (2) these patients were more likely to have a stroke than a TIA and their strokes were more severe; and (3) low DOAC exposure was more likely when the prescribed dosage regimens contradicted approved standards. Finally (4) DOAC exposure had a higher explanatory value for stroke severity than single concentrations at the time of the event.

In contrast to vitamin K antagonists (VKAs), the anticoagulatory effect of all DOACs is rapid and clearly concentration-dependent [19,20], suggesting that current plasma concentrations might reflect efficacy. In fact, concentration monitoring in pivotal DOAC trials revealed that important clinical endpoints, including ischemic stroke, correlated with DOAC trough concentrations [21,22]. Although, ischemic stroke in patients taking DOACs has become an important and common scenario in clinical practice, data on actual DOAC plasma concentrations in the hyperacute phase of stroke are still limited, not least because these events typically occur outside direct medical care. In particular, it is unknown whether DOAC plasma concentrations reflect clinical endpoints as well as they mirror the inhibition of coagulation factors.

Because of the importance of supporting acute treatment decisions, studies in DOAC-pretreated patients with acute cerebrovascular events have so far used the results of (faster) non-specific or specific coagulation tests as surrogates for DOAC activity at the time of hospital admission, which occurs often hours after event onset. However, such single activity testing of coagulation parameters represents an approximate snap-shot of the actual DOAC plasma concentration only, the quality of anticoagulation over time cannot be assessed and it gives no information on DOAC concentrations at the event time. To close this gap, we meticulously collected information on dose and time of individual drug intake, modeled individual concentration-time profiles, and correlated these data with the reported occurrence of stroke symptoms. This approach is independent of time of admission and considers the specifics of individual dosing regimens and drug intake. To the best of our knowledge, this is the first study designed to evaluate the significance of DOAC plasma concentrations at the event for the occurrence of TIA, ischemic strokes, and their severity.

Irrespective of the prescribed DOAC, exposure of our patients was lower than expected and in a substantial proportion of patients this was a result of physicians´ non-adherence to guidelines (underdosing: 24.2%) or non-adherence of patients (9%). Poor guideline adherence with a trend to lower doses has been reported repeatedly and is associated with less favorable outcomes [23]. Similarly, also poor patient adherence and persistence is alarmingly frequent and associated with increased frequencies of thromboembolic events and death [24]. In accordance with the known relationship between DOAC dose and efficacy, our study revealed that patients with DOAC plasma exposure below the population mean were 2.4 times more likely to have a stroke than a TIA and their strokes were more severe. This difference was maintained when comparing the DOAC concentrations at the time of the ischemic event but the latter did not correlate with stroke severity. Because no study evaluated DOAC plasma concentrations at the time of ischemic stroke before, comparison of our results to former reports is not possible. In descriptive analyses, a single-center registry study assessing associations between specific coagulation tests and clinical and imaging characteristics reported lower functional plasma levels at hospital admission to be associated with stroke severity, higher risk of persisting neurological deficits, and cerebral infarction on brain imaging [5]. Furthermore, in a multivariate model, low results of such coagulation tests were observed to be an independent predictor of LVO [5]. However, differences between TIA and ischemic stroke patients were not analysed and time-intervals between last DOAC administration and laboratory assessments was not reported [5]. By considering DOAC dosing schemes prior to ischemic stroke, another analysis of this registry revealed that only one in three patients with ischemic stroke followed appropriately dosed DOAC regimes and underdosing was associated to greater stroke severity and worse functional outcome after 3 months [25]. However, only patients using coagulation factor Xa inhibitors were included and no DOAC plasma concentrations were measured at admission thus leaving open whether exposure was indeed low. A further retrospective evaluation of patients with recurrent ischemic stroke at a tertiary stroke center reported that missed DOAC use for ≥48 hours or DOAC dosing below current recommendations were associated with LVO [26]. However again, no information on DOAC plasma concentrations was included [26].

Our study enabled to determine whether DOAC exposure over time (AUC) or the concentration at event onset better predicted ischemic stroke or its severity and found that AUC was superior to explain stroke severity than DOAC concentration at the time of the event. Our findings therefore suggest that lower DOAC exposure (and not the concentration at stroke onset) is more relevant for treatment efficacy. These results are in accordance to other reports, suggesting that AUC values might more closely reflect actual drug action than either peak [27,28] or trough DOAC concentrations [29]; however, this appears not to be true for all DOACs [29].

When transferring the broad knowledge on treatment with VKA to DOAC patients, our observation that stroke severity was better predicted by DOAC exposure than by plasma concentration at the event is plausible. Exposure to VKA can be expressed by INR values and it is well known that INR values of ≥2.0 reduce the frequency of ischemic stroke and its severity as well as the risk of death from stroke in patients with AF [30]. Moreover, the proportion of treatment time outside the therapeutic range has been clearly associated with adverse outcomes, including ischemic stroke [31,32]. Furthermore, admission INR values are inversely correlated with infarct volume on diffusion-weighted MRI imaging [33] and the intensity of anticoagulation with VKAs in cardioembolic ischemic stroke is inversely associated with the severity of neurological deficits and clinical outcome [34]. Such observations are the foundation for the generally known and established recommendation to thoroughly monitor INR values in VKA patients to gain optimal prevention for thromboembolic events by achieving adequate drug exposure.

Taken together, our findings have important clinical implications for the future management of patients using DOACs. By now, recommendations to perform regular controls of plasma concentrations in DOACs are absent. In contrast, our data suggest that monitoring these patients by thorough control of dosing schemes and long-term drug exposure might help identifying patients with inadequately low exposure and subsequent dose adaptions could considerably improve stroke prevention. Whether this strategy will help improving outcomes will now have to be tested in a prospective clinical trial.

An obvious strength of the present study is the prospective inclusion of a large number of anticoagulated acute ischemic stroke and TIA patients, irrespective of a prespecified DOAC treatment. Moreover, external standardization of DOAC plasma concentrations provided sound knowledge on associations between DOAC plasma concentrations and ischemic cerebrovascular events. However, modeled exposure in our study is based on the existent literature and limited to the last hours before stroke. Longitudinal exposure information was not available, the numbers of patients for single DOACs are limited and the single-center approach might restrict generalizing our results. It could be speculated that assuming a regular drug use in patients in whom the last drug intake was not reported or unknown may have generated bias. However, excluding these patients yielded even higher risks for stroke occurrence due to low DOAC plasma exposure (OR, 2.99, data not shown). This indicates that the main approach of approximating drug intake times in those patients can be considered as rather conservative as main results appeared to be biased towards the null, if at all. Due to the observational design, we did not perform brain imaging in a predefined structured manner and therefore cannot adequately report stroke volumes. Moreover, we did not consider details of recanalization therapies or changes in renal function when determining the functional outcome. Longitudinal information on exposure was not available, the number of patients with individual DOACs was small, as was the number of 1–2 observations per patient available for extrapolation of individual pharmacokinetics. The lack of a control group and of patients with intracerebral hemorrhages is a further limitation.

Conclusions

Our data reveal that low DOAC exposure is strongly associated to ischemic stroke and its severity. In consequence, monitoring plasma concentrations in DOAC patients on a regular basis might identify patients prone to ischemic stroke and subsequent dose adaptions could considerably improve preventing ischemic cerebrovascular events by achieving adequate DOAC exposure. Monitoring these patients longitudinally would also enable to examine associations between the impact of duration of suboptimal anticoagulation on the occurrence of cerebrovascular events. This should now be assessed in an adequately powered prospective trial.

Supplementary materials

Supplementary materials related to this article can be found online at https://doi.org/10.5853/jos.2020.04952.

Supplementary Table 1

Demographic findings, cardiovascular risk factors, and DOAC-specific findings of all included patients (n=211)

jos-2020-04952-suppl1.pdf
Supplementary Table 2.

Differences between low and high DOAC exposure (n=211)

jos-2020-04952-suppl2.pdf
Supplementary Table 3.

Multivariate logistic regression between high and low DOAC exposure (n=211)

jos-2020-04952-suppl3.pdf
Supplementary Table 4.

Multivariate logistic regression between patients with and without large vessel occlusions of all included patients and of patients with modeled DOAC concentration at the event

jos-2020-04952-suppl4.pdf
Supplementary Table 5.

Differences between low and high DOAC concentrations at onset of the event (n=179)

jos-2020-04952-suppl5.pdf
Supplementary Table 6.

Differences between favorable and unfavorable 3 months outcome in ischemic stroke patients (n=149)

jos-2020-04952-suppl6.pdf
Supplementary Table 7.

Multivariate logistic regression between favorable and unfavorable 3 months outcome in ischemic stroke patients with follow-up information (n=149)

jos-2020-04952-suppl7.pdf
Supplementary Figure 1.

Diagnostic plots to explore deviations of predicted concentrations from observed measurements at the respective times of measurement. (A, C) Panels visualize the bivariate relationship between predicted concentrations und observed measurements with distinct pairs of values colour-coding individual direct oral anticoagulants (DOACs) or whether the indicator was derived from a pharmacokinetic profile estimated by one or two measurements per patient. Residual plots illustrate residuals over time after the last DOAC administration (B) or residuals stratified by their origin from pharmacokinetic profiles with one-point or two-point estimation (D).

jos-2020-04952-suppl8.pdf
Supplementary Figure 2.

Forest plot of effect estimates (expressed as odds ratios) of either low direct oral anticoagulant (DOAC) plasma exposure (A) or low DOAC plasma concentration at the time of the event (B) on the risk of ischemic stroke. Different outcome definitions are plotted on the discrete y-axis with the results from the main analysis, a restricted sample including only patients with atrial fibrillation (AF), and a sample including only cardioembolic ischemic strokes. Whiskers around the point estimates visualize 95% confidence intervals. AUC, area under the curve.

jos-2020-04952-suppl9.pdf
Supplementary Figure 3.

Substance-specific modulation of overall direct oral anticoagulant (DOAC) effects expressed as odds ratio for the probability of stroke occurrence (left side, A, C) or as a linear increase in the National Institutes of Health Stroke Scale (NIHSS) score for stroke severity (right side, B, D). The top row indicates substance-specific effects around the DOAC class effect (dashed vertical line) in terms of exposure (area under the curve [AUC]), while the bottom row indicated substance-specific effects in term of concentration at the time of the event. Whiskers visualize 95% confidence intervals of the respective estimates.

jos-2020-04952-suppl10.pdf

Notes

Disclosure

Timolaos Rizos received consulting honoraria, speakers’ honoraria and travel support from Bristol-Myers Squibb/Pfizer, Boehringer Ingelheim, Bayer HealthCare, Portola Pharmaceuticals and Daiichi Sankyo, outside of the present work; Jan Purrucker received consulting honoraria and travel support from Abbott, Akcea, Boehringer Ingelheim, Daiichi Sankyo and Pfizer, outside of the present work; David Czock received honoraria and travel support from Biogen GmbH, Chiesi GmbH and Daiichi Sankyo, outside of the present work; Peter A. Ringleb and Walter E. Haefeli received honoraria from Bayer, Boehringer Ingelheim, Daiichi Sankyo and Pfizer, outside of the present work.

Andreas D. Meid, Chris Dumschat, Andrea Huppertz, Kathrin I. Foerster, Jürgen Burhenne, and Ekkehart Jenetzky report no disclosures.

Acknowledgements

The work of Andreas D. Meid was supported by a grant from the Medical Faculty of Heidelberg University (Physician Scientist program).

References

1. Purrucker JC, Steiner T. Management of acute stroke in patients on oral anticoagulants. Curr Opin Neurol 2017;30:1–7.
2. Purrucker JC, Haas K, Rizos T, Khan S, Poli S, Kraft P, et al. Coagulation testing in acute ischemic stroke patients taking non-vitamin k antagonist oral anticoagulants. Stroke 2017;48:152–158.
3. Seiffge DJ, Traenka C, Polymeris A, Hert L, Fisch U, Peters N, et al. Feasibility of rapid measurement of rivaroxaban plasma levels in patients with acute stroke. J Thromb Thrombolysis 2017;43:112–116.
4. Volbers B, Köhrmann M, Kallmünzer B, Kurka N, Breuer L, Ringwald J, et al. Dabigatran plasma levels in acute cerebrovascular events. J Stroke Cerebrovasc Dis 2016;25:877–882.
5. Macha K, Marsch A, Siedler G, Breuer L, Strasser EF, Engelhorn T, et al. Cerebral ischemia in patients on direct oral anticoagulants. Stroke 2019;50:873–879.
6. Seiffge DJ, Kägi G, Michel P, Fischer U, Béjot Y, Wegener S, et al. Rivaroxaban plasma levels in acute ischemic stroke and intracerebral hemorrhage. Ann Neurol 2018;83:451–459.
7. Gosselin RC, Adcock DM, Bates SM, Douxfils J, Favaloro EJ, Gouin-Thibault I, et al. International Council for Standardization in Haematology (ICSH) recommendations for laboratory measurement of direct oral anticoagulants. Thromb Haemost 2018;118:437–450.
8. Foerster KI, Huppertz A, Müller OJ, Rizos T, Tilemann L, Haefeli WE, et al. Simultaneous quantification of direct oral anticoagulants currently used in anticoagulation therapy. J Pharm Biomed Anal 2018;148:238–244.
9. Byon W, Sweeney K, Frost C, Boyd RA. Population pharmacokinetics, pharmacodynamics, and exploratory exposure-response analyses of apixaban in subjects treated for venous thromboembolism. CPT Pharmacometrics Syst Pharmacol 2017;6:340–349.
10. Liesenfeld KH, Lehr T, Dansirikul C, Reilly PA, Connolly SJ, Ezekowitz MD, et al. Population pharmacokinetic analysis of the oral thrombin inhibitor dabigatran etexilate in patients with non-valvular atrial fibrillation from the RE-LY trial. J Thromb Haemost 2011;9:2168–2175.
11. Krekels EH, Niebecker R, Karlsson MO, Miller R, Shimizu T, Karlsson KE, et al. Population pharmacokinetics of edoxaban in patients with non-valvular atrial fibrillation in the ENGAGE AF-TIMI 48 study, a phase III clinical trial. Clin Pharmacokinet 2016;55:1079–1090.
12. Girgis IG, Patel MR, Peters GR, Moore KT, Mahaffey KW, Nessel CC, et al. Population pharmacokinetics and pharmacodynamics of rivaroxaban in patients with non-valvular atrial fibrillation: results from ROCKET AF. J Clin Pharmacol 2014;54:917–927.
13. Cirincione B, Kowalski K, Nielsen J, Roy A, Thanneer N, Byon W, et al. Population pharmacokinetics of apixaban in subjects with nonvalvular atrial fibrillation. CPT Pharmacometrics Syst Pharmacol 2018;7:728–738.
14. Stangier J. Clinical pharmacokinetics and pharmacodynamics of the oral direct thrombin inhibitor dabigatran etexilate. Clin Pharmacokinet 2008;47:285–295.
15. Yin OQ, Tetsuya K, Miller R. Edoxaban population pharmacokinetics and exposure-response analysis in patients with non-valvular atrial fibrillation. Eur J Clin Pharmacol 2014;70:1339–1351.
16. Assessment report Xarelto VHF European Medicines Agency; https://www.ema.europa.eu/en/documents/variation-report/xarelto-h-c-944-ii-0012-epar-assessment-report-variation_en.pdf. 2012. Accessed November 24, 2021.
17. Yeo IK, Johnson RA. A new family of power transformations to improve normality or symmetry. Biometrika 2000;87:954–959.
18. Vuong QH. Likelihood ratio tests for model selection and non-nested hypotheses. Econometrica 1989;57:307–333.
19. Wang X, Tirucherai G, Marbury TC, Wang J, Chang M, Zhang D, et al. Pharmacokinetics, pharmacodynamics, and safety of apixaban in subjects with end-stage renal disease on hemodialysis. J Clin Pharmacol 2016;56:628–636.
20. Kubitza D, Becka M, Wensing G, Voith B, Zuehlsdorf M. Safety, pharmacodynamics, and pharmacokinetics of BAY 59-7939–an oral, direct Factor Xa inhibitor: after multiple dosing in healthy male subjects. Eur J Clin Pharmacol 2005;61:873–880.
21. Reilly PA, Lehr T, Haertter S, Connolly SJ, Yusuf S, Eikelboom JW, et al. The effect of dabigatran plasma concentrations and patient characteristics on the frequency of ischemic stroke and major bleeding in atrial fibrillation patients: the RE-LY Trial (Randomized Evaluation of Long-Term Anticoagulation Therapy). J Am Coll Cardiol 2014;63:321–328.
22. Ruff CT, Giugliano RP, Braunwald E, Morrow DA, Murphy SA, Kuder JF, et al. Association between edoxaban dose, concentration, anti-factor Xa activity, and outcomes: an analysis of data from the randomised, double-blind ENGAGE AF-TIMI 48 trial. Lancet 2015;385:2288–2295.
23. Steinberg BA, Shrader P, Pieper K, Thomas L, Allen LA, Ansell J, et al. Frequency and outcomes of reduced dose non-vitamin K antagonist anticoagulants: results from ORBIT-AF II (The Outcomes Registry for Better Informed Treatment of Atrial Fibrillation II). J Am Heart Assoc 2018;7:e007633.
24. Borne RT, O‘Donnell C, Turakhia MP, Varosy PD, Jackevicius CA, Marzec LN, et al. Adherence and outcomes to direct oral anticoagulants among patients with atrial fibrillation: findings from the veterans health administration. BMC Cardiovasc Disord 2017;17:236.
25. Stoll S, Macha K, Marsch A, Gerner ST, Siedler G, Fröhlich K, et al. Ischemic stroke and dose adjustment of oral Factor Xa inhibitors in patients with atrial fibrillation. J Neurol 2020;267:2007–2012.
26. Woo HG, Chung I, Gwak DS, Kim BK, Kim BJ, Bae HJ, et al. Recurrent ischemic stroke in atrial fibrillation with non-vitamin K antagonist oral anticoagulation. J Clin Neurosci 2019;64:127–133.
27. Herink MC, Zhuo YF, Williams CD, DeLoughery TG. Clinical management of pharmacokinetic drug interactions with direct oral anticoagulants (DOACs). Drugs 2019;79:1625–1634.
28. Leil TA, Feng Y, Zhang L, Paccaly A, Mohan P, Pfister M. Quantification of apixaban’s therapeutic utility in prevention of venous thromboembolism: selection of phase III trial dose. Clin Pharmacol Ther 2010;88:375–382.
29. Sakaguchi T, Osanai H, Murase Y, Ishii H, Nakashima Y, Asano H, et al. Monitoring of anti-Xa activity and factors related to bleeding events: a study in Japanese patients with nonvalvular atrial fibrillation receiving rivaroxaban. J Cardiol 2017;70:244–249.
30. Hylek EM, Go AS, Chang Y, Jensvold NG, Henault LE, Selby JV, et al. Effect of intensity of oral anticoagulation on stroke severity and mortality in atrial fibrillation. N Engl J Med 2003;349:1019–1026.
31. Gallagher AM, Setakis E, Plumb JM, Clemens A, van Staa TP. Risks of stroke and mortality associated with suboptimal anticoagulation in atrial fibrillation patients. Thromb Haemost 2011;106:968–977.
32. Connolly SJ, Pogue J, Eikelboom J, Flaker G, Commerford P, Franzosi MG, et al. Benefit of oral anticoagulant over antiplatelet therapy in atrial fibrillation depends on the quality of international normalized ratio control achieved by centers and countries as measured by time in therapeutic range. Circulation 2008;118:2029–2037.
33. Ay H, Arsava EM, Gungor L, Greer D, Singhal AB, Furie KL, et al. Admission international normalized ratio and acute infarct volume in ischemic stroke. Ann Neurol 2008;64:499–506.
34. Nakamura A, Ago T, Kamouchi M, Hata J, Matsuo R, Kuroda J, et al. Intensity of anticoagulation and clinical outcomes in acute cardioembolic stroke: the Fukuoka Stroke Registry. Stroke 2013;44:3239–3242.

Article information Continued

Figure 1.

Model-based concentration-time profile of an exemplary patient treated with apixaban (red line). Shaded areas (purple) visualize the typical range of 95% virtual patients as predicted from the underlying population pharmacokinetic model for apixaban while accounting for available covariates of the particular patient or assuming median values [9]. Red dots represent the actual measurements in the particular patient and the red line shows the a posteriori estimates of the individual concentration-time profile based on individual pharmacokinetic parameters of the particular patient. The concentration at the time of event (+black cross) was directly derived from the profile.

Figure 2.

Selection of the study population. DOAC, direct oral anticoagulant; TIA, transient ischemic attack; AUC, area under the curve; NA, non-adherence.

Figure 3.

Box plots of the ratio of modeled area under the curve (AUC) and published figures corresponding to the respective doses with superimposed individual AUC values of all included patients. The broken line indicates the threshold between relatively higher and lower direct oral anticoagulant (DOAC) exposure than the published expected average. Each boxplot contains the median (horizontal line in the box), the upper quartile (75th percentile, top of box), the lower quartile (25th percentile, bottom of box). The whiskers plot the minimal and maximal DOAC exposure. Solid circles: Dot plots of categorized individual AUC results of all included patients (lower limit of quantification for ultraperformance liquid chromatography–tandem mass spectrometry 1 ng/mL). Open circles: non-adherent patients (concentration <1 ng/mL; n=19).

Table 1.

Multivariate logistic regression between ischemic stroke and TIA patients in all included patients and in all patients with modeled DOAC concentration at the event

Variable Ischemic stroke vs. TIA
All included patients (n=211)
Patients with extrapolated DOAC concentration at the event (n=179)
OR 95% CI P OR 95% CI P
Low DOAC exposure* 2.411 1.254–4.638 0.008 - - -
Low DOAC concentration at event* - - - 4.123 1.834–9.268 0.001
Female sex 0.610 0.294–1.270 0.187 0.528 0.244–1.145 0.106
Age 0.961 0.917–1.008 0.100 0.983 0.935–1.034 0.506
Hypertension 1.132 0.365–3.516 0.830 0.722 0.204–2.555 0.614
Diabetes mellitus 0.487 0.235–1.008 0.053 0.340 0.149–0.774 0.010
Hypercholesterolemia 0.657 0.320–1.347 0.251 0.726 0.324–1.624 0.435
Previous stroke/TIA 1.066 0.531–2.142 0.857 0.956 0.440–2.078 0.956
Congestive heart failure 1.428 0.576–3.540 0.442 1.436 0.534–3.857 0.473
Vascular disease 0.940 0.445–1.989 0.872 1.047 0.460–2.385 0.912
Atrial fibrillation 0.880 0.254–3.049 0.076 0.907 0.229–3.600 0.890

TIA, transient ischemic attack; DOAC, direct oral anticoagulant; OR, odds ratio; CI, confidence interval.

*

Derived from area under the curve (AUC) ratios normalized to reference populations;

No published data available for comparison.

Table 2.

Multivariate linear regression of stroke severity (NIHSS) in all included ischemic stroke patients and in all ischemic stroke patients with modeled DOAC concentration at the event

Variable NIHSS
All included ischemic stroke patients (n=157)
Ischemic stroke patients with extrapolated DOAC concentration at the event (n=131)
Slope 95% CI P Slope 95% CI P
Low DOAC exposure* 3.161 0.741 to 5.581 0.011 - - -
Low DOAC concentration at event* - - - 1.570 –1.320 to 4.334 0.293
Female sex –2.708 –5.125 to –0.291 0.028 –2.559 –5.398 to 0.281 0.077
Age 0.076 –0.053 to 0.205 0.246 0.099 –0.060 to 0.257 0.220
Hypertension 0.103 –3.717 to 3.924 0.957 0.379 –3.875 to 4.633 0.860
Diabetes mellitus 0.380 –2.312 to 3.072 0.781 0.154 –3.214 to 3.522 0.928
Hypercholesterolemia –1.629 –4.039 to 0.780 0.183 –1.454 –4.318 to 1.409 0.317
Previous stroke/TIA 1.737 –0.769 to 4.243 0.173 1.661 –1.322 to 4.644 0.272
Congestive heart failure 1.001 –1.924 to 3.927 0.500 1.119 –2.383 to 4.622 0.528
Vascular disease 1.212 –1.292 to 3.717 0.340 1.354 –1.648 to 4.356 0.374
Large vessel occlusion 7.773 5.449 to 10.098 <0.001 7.597 4.796 to 10.397 <0.001
Atrial fibrillation 2.285 –1.687 to 6.257 0.257 1.936 –2.977 to 6.850 0.437

NIHSS, National Institutes of Health Stroke Scale; DOAC, direct oral anticoagulant; CI, confidence interval; TIA, transient ischemic attack.

*

Derived from ratios normalized to reference populations;

No published data available for comparison.