Article Text

Original research
Expediting workflow in the acute stroke pathway for endovascular thrombectomy in the northern Netherlands: a simulation model
  1. Willemijn J Maas1,2,
  2. Maarten M H Lahr2,
  3. Maarten Uyttenboogaart1,3,
  4. Erik Buskens2,4,
  5. Durk-Jouke van der Zee2,4
  6. On behalf of CONTRAST investigators
  1. 1Department of Neurology, University of Groningen, University Medical Centre Groningen, Groningen, The Netherlands
  2. 2Health Technology Assessment, Department of Epidemiology, University of Groningen, University Medical Centre Groningen, Groningen, The Netherlands
  3. 3Department of Radiology, Medical Imaging Centre, University of Groningen, University Medical Centre Groningen, Groningen, The Netherlands
  4. 4Department of Operations, Faculty of Economics & Business, University of Groningen, Groningen, The Netherlands
  1. Correspondence to Willemijn J Maas; w.j.maas{at}umcg.nl

Abstract

Objective The objective of this study is to identify barriers for the timely delivery of endovascular thrombectomy (EVT) and to investigate the effects of potential workflow improvements in the acute stroke pathway.

Design Hospital data prospectively collected in the MR CLEAN Registry were linked to emergency medical services data for each EVT patient and used to build two Monte Carlo simulation models. The ‘mothership (MS) model’, reflecting patients who arrived directly at the comprehensive stroke centre (CSC); and the ‘drip and ship’ (DS) model, reflecting patients who were transferred to the CSC from primary stroke centres (PSCs).

Setting Northern region of the Netherlands. One CSC provides EVT, and its catchment area includes eight PSCs.

Participants 248 patients who were treated with EVT between July 2014 and November 2017.

Outcome measures The main outcome measures were total delay from stroke onset until groin puncture, functional independence at 90 days (modified Rankin Scale 0–2) and mortality.

Results Barriers identified included fast-track emergency department routing, prealert for transfer to the CSC, reduced handover time between PSC and ambulance, direct transfer from CSC arrival to angiography suite entry, and reducing time to groin puncture. Taken together, all workflow improvements could potentially reduce the time from onset to groin puncture by 59 min for the MS model and 61 min for the DS model. These improvements could thus result in more patients—3.7% MS and 7.4% DS—regaining functional independence after 90 days, in addition to decreasing mortality by 3.0% and 5.0%, respectively.

Conclusions In our region, the proposed workflow improvements might reduce time to treatment by about 1 hour and increase the number of patients regaining functional independence by 6%. Simulation modelling is useful for assessing the potential effects of interventions aimed at reducing time from onset to EVT.

  • stroke
  • organisation of health services
  • organisational development

Data availability statement

No data are available. The data for this sub-study from the MR CLEAN Registry and the data of the EMS are not publicly available, as they allow for the identification of individual centres. The sharing of such data is in conflict with the privacy regulations in the Netherlands.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Data were collected on time delays along the acute stroke pathway for patients treated with endovascular thrombectomy (EVT), thereby allowing the identification, analysis and simulation of barriers from onset to treatment.

  • An extensive set of workflow improvements is suggested based on data analysis, expert opinion and literature.

  • A simulation model of the acute stroke pathway is developed, enabling the effective and efficient assessment of workflow improvements, relying on realistic in-silico modelling.

  • The simulation model includes only patients treated with EVT in a region with one comprehensive stroke centre, but it could be extended to all suspected patients who had a stroke, thereby allowing a more comprehensive assessment of stroke care.

Introduction

Acute ischaemic stroke places a large burden on society, and the overall incidence has increased by 78% since 1990.1 The main reperfusion treatments for acute ischaemic stroke due to large vessel occlusion are intravenous thrombolysis (IVT) and endovascular thrombectomy (EVT). The phrase ‘time is brain’ applies to both treatments. For EVT, the probability of regaining functional independence at 90 days after stroke declines by 5%–6% for each additional hour delay from onset to groin puncture (OTG).2 3

Successful and timely EVT largely depends on the regional organisation of acute stroke care delivery. Delays that can occur during prehospital and intrahospital processes, as well as along each step in the acute stroke pathway, have the potential to worsen patient outcomes or even rule out the possibility of acute treatment. Pathway elements that have been identified as having the potential to cause treatment delays include prehospital stroke management, in-hospital patient transfer, anaesthetic management, teamwork and inter-hospital patient transfer.4

Most studies of interventions aimed at improving workflow processes have focused on specific interventions, examining bits and pieces of the acute stroke pathway separately. The joint analysis of several improvements might lead to the identification of actual improvements. Simulation modelling has been suggested as a means of supporting such comprehensive analyses, and it has been performed within the context of IVT based on a variety of organisational models.5 6

The objectives of this study are (1) to assess delays in the workflow of acute stroke care, based on patient-level data; and (2) to estimate the impact of reducing delays throughout the process, from work-up to EVT treatment, based on simulation modelling.

Methods

Setting

This study is based on prospective data collected in the MR CLEAN Registry7 from patients treated with EVT in one comprehensive stroke centre (CSC), which provides EVT for eligible patients in the northern part of the Netherlands (1.7 million inhabitants). Its catchment area includes eight primary stroke centres (PSCs), spaced at distances of 6–84 km, as shown in online supplemental figure S1.

Participants and data collection

Between July 2014 and November 2017, 285 patients were included. According to the emergency medical services (EMS) protocol,8 patients suspected of acute stroke were routed to the nearest IVT-capable hospital. The patients were either sent directly to a CSC (mothership (MS) model) or first presented at a PSC and subsequently transferred to the CSC for EVT (drip and ship (DS) model). In the eastern part of the province of Groningen, patients were routed directly to the CSC, reflecting a centralised organisational model.9

Patient data on clinical characteristics, diagnostic processes, time delays and ambulance routing patterns were used as input for simulation modelling. In-hospital time delays included onset or time last seen well, CT, IVT initiation, CT angiogram (CTA), arrival at the angiography suite and the time of groin puncture. In-hospital (PSC or CSC) patients were routed through the emergency department (ED) according to three routes: (1) CT to IVT to CTA; (2) CT to CTA to IVT and (3) CT to CTA (patients ineligible for IVT). Following secondary transfer, DS patients arriving at the CSC could undergo additional diagnostics (eg, CT and/or CTA).

Prehospital data from three EMS organisations were collected retrospectively and linked to the MR CLEAN Registry data for each patient. Time-delay items collected included 911 notification, EMS arrival at the stroke-onset location, departure to hospital and arrival at hospital. Additional data collected for DS patients included the timestamps for EMS transfer notification, arrival at PSC, departure to CSC and arrival at CSC.

Patients were excluded from analyses in case of a prior modified Rankin Scale (mRS) >2 and when OTG exceeded 390 min, as EVT based on perfusion CT beyond 6 hours was not indicated at that time. Missing values were excluded from analyses.

Patient and public involvement

No patients involved.

Simulation

Separate Monte Carlo simulation models were developed for the MS and DS organisation models.10 Prior to model building, conceptual modelling was performed in order to abstract real-world acute stroke pathways, as shown in figure 1. Conceptual models were validated using expert opinion (MU), combined with literature observations and input from stroke experts participating in the national collaboration for new treatments of acute stroke (CONTRAST) consortium.11

Figure 1

Conceptual models of the acute stroke pathway: ‘mothership’ and ‘drip and ship’. CT, computed tomography; CTA, CT angiography; EMS, emergency medical services; EVT, endovascular thrombectomy; IVT, intravenous thrombolysis; POC, point of care.

Both simulation models were developed using Plant Simulation.12 Distributions for the individual time-delay variables were based on patient data and obtained using ExpertFit.13 Details are presented as online supplemental tables 1 and 2.

Modelling scenarios

We identified barriers along the acute stroke pathway by analysing patient data, relevant literature and expert opinion (MU). These barriers were used to create hypothetical scenarios, which we tested ‘in silico’ using the simulation model developed for this purpose.

Outcome measures

Outcome measures include OTG, likelihood of functional independence (mRS 0–2) and mortality (mRS 6) at 90 days.

Analysis

The simulation models were validated numerically by comparing mean, median, SD, minimum and maximum time values of real-world patient data and observations to model data and outputs.

Within the simulation model, ordinal logistic regression was used to estimate the likelihood of each of the seven scales belonging to the mRS score, ranging from 0 (no symptoms) to 6 (death). Known prognostic variables were OTG (continuous), age (continuous), National Institutes of Health Stroke Scale score (continuous) and CTA collateral grading score in four categories (absence of collaterals, less than 50% filling of the occluded area, more than 50% filling but less than 100% filling of the occluded area and 100% filling of the occluded area). The likelihood of functional independence (mRS 0–2) was calculated from the formulas obtained by ordinal logistic regression, using IBM SPSS Statistics V.23 software. Details are presented as online supplemental material 1.

For each scenario, we calculated the clinical benefits in terms of reduction in OTG and the likelihood of regaining functional independence and reducing mortality. Significance testing was inappropriate, as the goal was to assess the potential gain expected based on 100 000 hypothetical patients, rather than to test a hypothesis as in an actual experiment.

Results

In all, 248 patients met the inclusion criteria. Of these patients, 27 were excluded because of a prestroke mRS>2, and/or an unknown OTG of >390 min (12 patients). Patient characteristics, diagnostics and median time delays for each model are presented in table 1. For MS patients (n=83), the median (IQR) OTG was 205 (160–260) min; 51.8% regained functional independence after 90 days and mortality was 26.5%. For DS patients (n=165), the respective figures were 230 (198–275) min, 52.1% and 22.4%. To obtain the likelihood formulas for each of the seven mRSs, data from 80 MS patients and 154 DS patients were used. Despite faster OTG, the MS patients had a lower likelihood of functional independence and a higher likelihood of mortality after 90 days compared with DS patients.

Table 1

Characteristics, diagnostics and time delays of the MS and DS models

Identified delays

We identified multiple opportunities for improving workflow for both the DS and MS models.

DS model, PSC workflow

The door-in-door-out (DIDO) time was used to estimate the entire PSC workflow, defined as time from PSC arrival until departure to the CSC. The DIDO time of patients routed through the ED according to route 2 (CT to CTA to IVT) was less than that of patients routed according to route 1 (CT to IVT to CTA), with a mean (SD) of 82 (25) min versus 100 (37) min, respectively.

We also assessed the handover time from PSC to ambulance for transfer to the CSC. The lowest median (IQR) handover time in one of the PSCs was 11 (8–14) min, as compared with an overall median time of 14 (10–16) min.

DS model, CSC workflow

If no additional diagnostics are required, DS patients arriving at the CSC should be transferred directly to the angiography suite.14 The observed median (IQR) transfer time from CSC arrival to angiography suite was 26 (16–38) min, and from angiography suite arrival to groin puncture 30 (24–35) min.

MS model, CSC workflow

We assessed the time from CSC presentation to arrival at the angiography suite for each route through the ED. Patients who were routed according to route 2 (CT to CTA to IVT) had shorter delays compared with those who were routed according to route 1 (CT to IVT to CTA); with a mean (SD) of 103 (46) min compared with 113 (42) min, respectively. The observed median (IQR) time from the last examination at the ED to angiography suite arrival was 58 (44–82) min, and between angiography suite arrival and groin puncture 28 (25–35) min.

Modelling scenarios

The following scenarios were defined, based on the barriers identified for the DS model (online supplemental table S3): routing all patients without contraindication for IVT through the ED according to route 2 (CT to CTA to IVT) (scenario 1a); EMS prealert is used, thus reducing the ambulance response time to 0 min (scenario 1b); reducing the handover time from PSC to ambulance to 11 min (scenario 1c); and combining all three experiments (scenario 1d).

The following scenarios were considered for the CSC optimised workflow improvements (DS model): direct transfer from CSC arrival to the angiography suite (maximum of 5 min, scenario 2a); reducing the time from angiography suite arrival to groin puncture to 10 min, based on expert opinion, analysis of the MR CLEAN Registry dataset for all hospitals in the Netherlands, and a previously published study15 (scenario 2b); and combining the two experiments (scenario 2c). In addition, the PSC and CSC workflow improvements were combined into one experiment (scenario 3).

The scenarios for the MS model were as follows: routing all patients without contraindication for IVT through the ED according to route 2 (CT to CTA to IVT; scenario 4a); reducing time from last examination at the ED to angiography suite arrival to a maximum of 30 min (scenario 4b); and reducing the time from angiography suite arrival to groin puncture to a maximum of 10 min (scenario 4c). Scenarios 4a and 4b are based on expert opinion, analysis of the MR CLEAN Registry dataset on all hospitals in the Netherlands, and a previously published paper.2 In scenario 4d, all experiments were combined.

Simulation results

DS workflow

Implementing all workflow improvements in a PSC (scenario 1d) would imply an absolute increase of 2.2% in the number of patients regaining functional independence after 90 days, a mortality reduction of 1.5%, and a reduction in OTG of 18 min (table 2). Realising workflow improvements within the CSC (scenario 2 c) would reduce OTG by 43 min, increase the proportion of patients reaching functional independence at 90 days by 5.3% and reduce mortality by 3.6%. Combining all workflow improvements in both PSC and CSC (scenario 3) would reduce OTG by 61 min, increase the proportion of patients reaching functional independence by 7.4% and decrease mortality by 5.0%.

Table 2

Simulation results

MS Workflow

Implementing all workflow improvements (scenario 4d) would reduce OTG by 59 min, increase the number of patients regaining functional independence at 90 days by 3.7% and decrease mortality by 3.0%.

The shifts in likelihood for each mRS score when all workflow improvements are executed in the DS and MS models are displayed in figure 2.

Figure 2

Shifts in likelihood for each mRS score when all workflow improvements are executed in the DS and MS models. DS, ‘drip and ship’ model; mRS, modified Rankin Scale; MS, ‘mothership’ model.

Discussion

The results of this study demonstrate that simulation modelling can be used to identify barriers for timely EVT and to assess the impact of workflow improvements in regional acute stroke care systems. Workflow improvements (eg, ED routing of CT to CTA to IVT, prealerting the ambulance, reducing handover time between PSC and EMS, and reducing CSC workflow from hospital arrival to groin puncture) could possibly reduce the time to EVT by approximately 1 hour. For DS patients, we estimate that the suggested workflow improvements could reduce OTG by 61 min, ultimately decreasing mortality by 5.0% and increasing the number of patients regaining functional independence at 90 days by 7.4%. The implementation of all hypothetical PSC workflow improvements for DS patients could make it possible to achieve the DIDO target time value of 75 min.2 15 For MS patients, the proposed interventions could reduce OTG by 59 min, decrease mortality by 3.0% and increase the number of patients regaining functional independence at 90 days by 3.7%.

For the aforementioned improvements, we specifically considered the acute stroke pathway of our region and the potential improvements that we systematically implemented ‘in silico’. Analysis of the MR CLEAN Registry for all hospitals in the Netherlands nevertheless revealed that some hospitals have already attained the level of our proposed improvements, while others have not. This suggests that the implementation of the proposed improvements could result in even greater benefits and that the selection of policies and improvements will depend on regional set-up and characteristics of existing acute stroke care systems.

The findings for the DS model indicate slightly greater improvement than has been reported in previous studies, while those for the MS model indicate slightly less improvement, with the number of patients regaining functional independence increasing by between 5% and 6% for each hour reduction in OTG.2 3 Possible explanations for the difference between our region and other regions might have to do with the fact that data in other studies were collected shortly after the introduction of EVT was newly introduced, as well as with region-specific differences (eg, hospital infrastructure). Furthermore, the use of ordinal logistic regression revealed greater fluctuations in estimating the likelihood of mRS in the DS model, as compared with the MS model. Possible explanations include the fact that a separate ordinal logistic regression was performed for each model, the small sample size (ie, n=154 for the DS model and n=80 for the MS model), and the fact that previous studies have not analysed data in separate routing groups (ie, the DS model vs the MS model).2 3 Another striking result was the higher probability of death and poor functional outcome for MS patients, despite a decrease in OTG. One possible explanation could be that patients with highly complex comorbidity and ischaemic stroke were more likely to be transferred directly to the CSC instead of to a PSC.

The results of our study can be generalised in part to other regions. Suggested improvements for the acute stroke pathway may be related to a generic conceptual model of care delivery that is consistent with many existing regional pathways and that faces similar challenges. While the impact of these improvements within specific regions will differ, they can jointly create a relevant starting point for optimising stroke systems. The most important benefit of the proposed simulation modelling study is that it allows the testing of potential improvements and the estimation of their impact for specific regions. As suggested by guidelines, and taking regional and patient characteristics into account,16 simulation modelling may be particularly useful for re-populating the generic model (ie, using conceptual models and patient data from other regions). In addition, simulation modelling might be an attractive option in terms of efficiency, as it starts with hypothetical improvements without immediately requiring investments and costs associated with hardware and organisation. Although it cannot completely replace RCTs, simulation modelling can be useful as a precursor to clinical studies, as a tool for organisational learning, and as a design approach (eg, for acute stroke care).17 18

Limitations

Our study is subject to several limitations. The simulation model includes only the acute stroke pathway for patients with large vessel occlusion. Ideally, a simulation model should take all suspected patients who had a stroke into account, thereby allowing a more comprehensive assessment of stroke care.

In addition, as a consequence of identifying the optimal ED routing for timely EVT, additional delays for administering IVT were not taken into account. For patients with large vessel occlusion, rapid IVT administration is associated with less disability at 90 days.19 Furthermore, many questions remain unanswered with regard to the most beneficial treatment for these occlusion patients: faster IVT and fast EVT; faster EVT with increased delay for IVT; or direct EVT without IVT. Direct EVT is currently being studied in the MR CLEAN NO-IV (ISRCTN80619088)20 and the SWIFT DIRECT (NCT03192332)21 trials. The recently published DIRECT-MT study reports that direct EVT was non-inferior compared with IVT and EVT.22 Until this question is answered, it will be necessary to balance the relative benefits of both treatments.

Conclusions

Simulation is useful in assessing the potential effects of reducing region-specific delays from OTG. In our region, potential workflow improvements could reduce the time to treatment by 1 hour, thereby increasing the number of patients regaining functional independence after 90 days by 8% (DS model) and 4% (MS model), in addition to decreasing mortality by 5% (DS model) and 3% (MS model).

Data availability statement

No data are available. The data for this sub-study from the MR CLEAN Registry and the data of the EMS are not publicly available, as they allow for the identification of individual centres. The sharing of such data is in conflict with the privacy regulations in the Netherlands.

Ethics statements

Patient consent for publication

Ethics approval

The MR CLEAN Registry data collection has been approved for the Netherlands by the central medical ethics committee and research board (MEC-2014-235). The need for individual patient consent was waived. A Data Transfer Agreement was drafted and implemented for purposes of linking hospital patient data to the corresponding EMS data.

Acknowledgments

We acknowledge the support of the Cardiovascular Research Initiative, part of the Dutch Heart Foundation (CVON2015-01: CONTRAST), the Brain Foundation Netherlands (HA2015.01.06), Health~Holland, Top Sector Life Sciences & Health (LSHM17016), Medtronic and Cerenovus. We also acknowledge the UMCG Emergency Medical Services, Kijlstra Emergency Medical Services and Emergency Medical Services Groningen.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors All authors designed the study. WJM, MMHL and MU gathered data. WJM and D-JvdZ analysed the data and made the simulation models. WJM wrote the draft of the manuscript and MMHL, MU, EB and D-JvdZ revised the manuscript for important intellectual content. WJM is the guarantor of this research and accepts full responsibility for the finished work and the conduct of the study, had access to the data and controlled the decision to publish.

  • Funding The CONTRAST consortium is supported by Netherlands Cardiovascular Research Initiative, an initiative of the Dutch Heart Foundation (CVON2015-01: CONTRAST) and by the Brain Foundation Netherlands. It is powered by Health~Holland, Top Sector Life Sciences, and it receives unrestricted funding from Medtronic and Cerenovus. Additional funding for this collaborative project is provided by the Netherlands Ministry of Economic Affairs through a PPP Allowance made available by the Top Sector Life Sciences & Health to stimulate public-private partnerships.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.