In-Silico Trials for Treatment of Acute Ischemic Stroke

Modelling midline shift and ventricle collapse in cerebral oedema following acute ischaemic stroke

In ischaemic stroke, a large reduction in blood supply can lead to the breakdown of the blood-brain barrier and to cerebral oedema after reperfusion therapy. The resulting fluid accumulation in the brain may contribute to a significant rise in intracranial pressure (ICP) and tissue deformation. Changes in the level of ICP are essential for clinical decision-making and therapeutic strategies. However, the measurement of ICP is constrained by clinical techniques and obtaining the exact values of the ICP has proven challenging. In this study, we propose the first computational model for the simulation of cerebral oedema following acute ischaemic stroke for the investigation of ICP and midline shift (MLS) relationship. The model consists of three components for the simulation of healthy blood flow, occluded blood flow and oedema, respectively. The healthy and occluded blood flow components are utilized to obtain oedema core geometry and then imported into the oedema model for the simulation of oedema growth. The simulation results of the model are compared with clinical data from 97 traumatic brain injury patients for the validation of major model parameters. Midline shift has been widely used for the diagnosis, clinical decision-making, and prognosis of oedema patients. Therefore, we focus on quantifying the relationship between ICP and midline shift (MLS) and identify the factors that can affect the ICP-MLS relationship. Three major factors are investigated, including the brain geometry, blood-brain barrier damage severity and the types of oedema (including rare types of oedema). Meanwhile, the two major types (stress and tension/compression) of mechanical brain damage are also presented and the differences in the stress, tension, and compression between the intraparenchymal and periventricular regions are discussed. This work helps to predict ICP precisely and therefore provides improved clinical guidance for the treatment of brain oedema.

MRI-based parameter inference for cerebral perfusion modelling in health and ischaemic stroke

Cerebral perfusion modelling is a promising tool to predict the impact of acute ischaemic stroke treatments on the spatial distribution of cerebral blood flow (CBF) in the human brain. To estimate treatment efficacy based on CBF, perfusion simulations need to become suitable for group-level investigations and thus account for physiological variability between individuals. However, computational perfusion modelling to date has been restricted to a few patient-specific cases. This study set out to establish automated parameter inference for perfusion modelling based on neuroimaging data and thus enable CBF simulations of groups. Magnetic resonance imaging (MRI) data from 75 healthy senior adults were utilised. Brain geometries were computed from healthy reference subjects' T1-weighted MRI. Haemodynamic model parameters were determined from spatial CBF maps measured by arterial spin labelling (ASL) perfusion MRI. Thereafter, perfusion simulations were conducted in 75 healthy cases followed by 150 acute ischaemic stroke cases representing an occlusion and CBF cessation in the left and right middle cerebral arteries. The anatomical fitness of the brain geometries was evaluated by comparing the simulated grey (GM) and white matter (WM) volumes to measurements in healthy reference subjects. Strong positive correlations were found in both tissue types (GM: Pearson's r 0.74, P<0.001; WM: Pearson's r 0.84, P<0.001). Haemodynamic parameter tuning was verified by comparing the total volumetric blood flow rate to the brain in healthy reference subjects and simulations (Pearson's r 0.89, P<0.001). In acute ischaemic stroke cases, the simulated infarct volume using a perfusion-based estimate was 197±25 ml. Computational predictions were in agreement with anatomical and haemodynamic values from the literature concerning T1-weighted, T2-weighted, and phase-contrast MRI measurements in healthy scenarios and acute ischaemic stroke cases. The acute stroke simulations did not capture small infarcts (left tail of the distribution), which could be explained by neglected compensatory mechanisms, e.g. collaterals. The proposed parameter inference method provides a foundation for group-level CBF simulations and for in silico clinical stroke trials which could assist in medical device and drug development.

Perfusion estimation using synthetic MRI-based measurements and a porous media flow model

The measurement of perfusion and filtration of blood in biological tissue give rise to important clinical parameters used in diagnosis, follow-up, and therapy. In this paper, we address techniques for perfusion analysis using processed contrast agent concentration data from dynamic MRI acquisitions. A new methodology for analysis is evaluated and verified using synthetic data generated on a tissue geometry.

A systematic review of cardiac in-silico clinical trials

Computational models of the heart are now being used to assess the effectiveness and feasibility of interventions through in-silico clinical trials (ISCTs). As the adoption and acceptance of ISCTs increases, best practices for reporting the methodology and analysing the results will emerge. Focusing in the area of cardiology, we aim to evaluate the types of ISCTs, their analysis methods and their reporting standards. To this end, we conducted a systematic review of cardiac ISCTs over the period of 1 January 2012-1 January 2022, following the preferred reporting items for systematic reviews and meta-analysis (PRISMA). We considered cardiac ISCTs of human patient cohorts, and excluded studies of single individuals and those in which models were used to guide a procedure without comparing against a control group. We identified 36 publications that described cardiac ISCTs, with most of the studies coming from the US and the UK. In 75% of the studies, a validation step was performed, although the specific type of validation varied between the studies. ANSYS FLUENT was the most commonly used software in 19% of ISCTs. The specific software used was not reported in 14% of the studies. Unlike clinical trials, we found a lack of consistent reporting of patient demographics, with 28% of the studies not reporting them. Uncertainty quantification was limited, with sensitivity analysis performed in only 19% of the studies. In 97% of the ISCTs, no link was provided to provide easy access to the data or models used in the study. There was no consistent naming of study types with a wide range of studies that could potentially be considered ISCTs. There is a clear need for community agreement on minimal reporting standards on patient demographics, accepted standards for ISCT cohort quality control, uncertainty quantification, and increased model and data sharing.

In silico thrombectomy trials for acute ischemic stroke

BACKGROUND AND OBJECTIVE

In silico trials aim to speed up the introduction of new devices in clinical practice by testing device design and performance in different patient scenarios and improving patient stratification for optimizing clinical trials. In this paper, we demonstrate an in silico trial framework for thrombectomy treatment of acute ischemic stroke and apply this framework to compare treatment outcomes in different subpopulations and with different thrombectomy stent-retriever devices. We employ a novel surrogate thrombectomy model to evaluate the thrombectomy success in the in silico trial.

METHODS

The surrogate thrombectomy model, built using data from a fine-grained finite-element model, is a device-specific binary classifier (logistic regression), to estimate the probability of successful recanalization, the outcome of interest. We incorporate this surrogate model within our previously developed in silico trial framework and demonstrate its use with three examples of in silico clinical trials. The first trial is a validation trial for the surrogate thrombectomy model. We then present two exploratory trials: one evaluating the performance of a commercially available device based on the fibrin composition in the occluding thrombus and one comparing the performance of two commercially available stent retrievers.

RESULTS

The Validation Trial showed the surrogate thrombectomy model was able to reproduce a similar recanalization rate as the real-life MR CLEAN trial (p=0.6). Results from the first exploratory trial showed that the chance of successful thrombectomy increases with higher blood cell concentrations in the thrombi, which is in line with observations from clinical data. The second exploratory trial showed improved recanalization success with a newer stent retriever device; however, these results require further investigation as the surrogate model for the newer stent retriever device has not yet been validated.

CONCLUSIONS

In this novel study, we have shown that in silico trials have the potential to help inform medical device developers on the performance of a new device and may also be used to select populations of interest for a clinical trial. This would reduce the time and costs involved in device development and traditional clinical trials.

Computational modelling of cerebral oedema and osmotherapy following ischaemic stroke

In ischaemic stroke, a large reduction in blood supply can lead to the breakdown of the blood brain barrier and to cerebral oedema after reperfusion therapy. Cerebral oedema is marked by elevated intracranial pressure (ICP), tissue herniation and reduced cerebral perfusion pressure. In clinical settings, osmotherapy has been a common practice to decrease ICP. However, there are no guidelines on the choice of administration protocol parameters such as injection doses, infusion time and retention time. Most importantly, the effects of osmotherapy have been proven controversial since the infusion of osmotic agents can lead to a range of side effects. Here, a new Finite Element model of brain oedema and osmotherapy is thus proposed to predict treatment outcome. The model consists of three components that simulate blood perfusion, oedema, and osmotherapy, respectively. In the perfusion model (comprising arteriolar, venous, and capillary blood compartments), an anatomically accurate brain geometry is used to identify regions with a perfusion reduction and potential oedema occurrence in stroke. The oedema model is then used to predict ICP using a porous circulation model with four fluid compartments (arteriolar blood, venular blood, capillary blood, and interstitial fluid). In the osmotherapy model, the osmotic pressure is varied and the changes in ICP during different osmotherapy episodes are quantified. The simulation results of the model show excellent agreement with available clinical data and the model is employed to study osmotherapy under various parameters. Consequently, it is demonstrated how therapeutic strategies can be proposed for patients with different pathological parameters based on simulations.

Modelling collateral flow and thrombus permeability during acute ischaemic stroke

The presence of collaterals and high thrombus permeability are associated with good functional outcomes after an acute ischaemic stroke. We aim to understand the combined effect of the collaterals and thrombus permeability on cerebral blood flow during an acute ischaemic stroke. A cerebral blood flow model including the leptomeningeal collateral circulation is used to simulate cerebral blood flow during an acute ischaemic stroke. The collateral circulation is varied to capture the collateral scores: absent, poor, moderate and good. Measurements of the transit time, void fraction and thrombus length in acute ischaemic stroke patients are used to estimate thrombus permeability. Estimated thrombus permeability ranges between 10^-7 and 10^-4 mm². Measured flow rates through the thrombus are small and the effect of a permeable thrombus on brain perfusion during stroke is small compared with the effect of collaterals. Our simulations suggest that the collaterals are a dominant factor in the resulting infarct volume after a stroke.

Quantification of hypoxic regions distant from occlusions in cerebral penetrating arteriole trees

The microvasculature plays a key role in oxygen transport in the mammalian brain. Despite the close coupling between cerebral vascular geometry and local oxygen demand, recent experiments have reported that microvascular occlusions can lead to unexpected distant tissue hypoxia and infarction. To better understand the spatial correlation between the hypoxic regions and the occlusion sites, we used both in vivo experiments and in silico simulations to investigate the effects of occlusions in cerebral penetrating arteriole trees on tissue hypoxia. In a rat model of microembolisation, 25 μm microspheres were injected through the carotid artery to occlude penetrating arterioles. In representative models of human cortical columns, the penetrating arterioles were occluded by simulating the transport of microspheres of the same size and the oxygen transport was simulated using a Green’s function method. The locations of microspheres and hypoxic regions were segmented, and two novel distance analyses were implemented to study their spatial correlation. The distant hypoxic regions were found to be present in both experiments and simulations, and mainly due to the hypoperfusion in the region downstream of the occlusion site. Furthermore, a reasonable agreement for the spatial correlation between hypoxic regions and occlusion sites is shown between experiments and simulations, which indicates the good applicability of in silico models in understanding the response of cerebral blood flow and oxygen transport to microemboli.

The brain function depends on the continuous oxygen supply through the bloodstream inside the microvasculature. Occlusions in the microvascular network will disturb the oxygen delivery in the brain and result in hypoxic tissues that can lead to infarction and cognitive dysfunction. To aid in understanding the formation of hypoxic tissues caused by micro-occlusions in the penetrating arteriole trees, we use rodent experiments and simulations of human vascular networks to study the spatial correlations between the hypoxic regions and the occlusion locations. Our results suggest that hypoxic regions can form distally from the occlusion site, which agrees with the previous observations in the rat brain. These distant hypoxic regions are primarily due to the lack of blood flow in the brain tissues downstream of the occlusion. Moreover, a reasonable agreement of the spatial relationship is found between the experiments and the simulations, which indicates the applicability of in silico models to study the effects of microemboli on the brain tissue.

A Review of the Advancements in the in-vitro Modelling of Acute Ischemic Stroke and Its Treatment

In-vitro neurovascular models of large vessel occlusions (LVOs) causing acute ischemic stroke (AIS) are used extensively for pre-clinical testing of new treatment devices. They enable physicians and engineers to examine device performance and the response of the occlusion to further advance design solutions for current unmet clinical needs. These models also enable physicians to train on basic skills, to try out new devices and new procedural approaches, and for the stroke team to practice workflows together in the comfort of a controlled environment in a non-clinical setting. Removal of the occlusive clot in its entirety is the primary goal of the endovascular treatment of LVOs via mechanical thrombectomy (MT) and the medical treatment via thrombolysis. In MT, recanalization after just one pass is associated with better clinical outcomes than procedures that take multiple passes to achieve the same level of recanalization, commonly known as first pass effect (FPE). To achieve this, physicians and engineers are continually investigating new devices and treatment approaches. To distinguish between treatment devices in the pre-clinical setting, test models must also be optimized and expanded become more nuanced and to represent challenging patient cohorts that could be improved through new technology or better techniques. The aim of this paper is to provide a perspective review of the recent advancements in the in-vitro modeling of stroke and to outline how these models need to advance further in future. This review provides an overview of the various in-vitro models used for the modeling of AIS and compares the advantages and limitations of each. In-vitro models remain an extremely useful tool in the evaluation and design of treatment devices, and great strides have been made to improve replication of physiological conditions. However, further advancement is still required to represent the expanding indications for thrombectomy and thrombolysis, and the generation of new thrombectomy devices, to ensure that smaller treatment effects are captured.

Comparison of stent retriever thrombectomy using 3-dimensional patient-specific models of intracranial circulation with actual middle cerebral artery occlusion thrombectomy cases

BACKGROUND AND PURPOSE

Stent retriever (SR) thrombectomy is commonly used for the treatment of emergent large vessel occlusion (ELVO) in acute ischemic stroke. Clot imaging parameters such as clot length, diameter, distance to the internal carotid artery terminus, and vessel angle where the SR is deployed may predict the likelihood of achieving first pass effect (FPE). Most of the proposed factors that seem to affect recanalization success have been studied individually, and conflicting data derived from clinical versus in vitro studies using 3-dimensional printed models of intracranial circulation currently exist.

METHODS

Using patient-specific 3-dimensional phantoms of the cervical and intracranial circulation, we simulated middle cerebral arteries (MCA) M1 and M2 occlusions treated with SR thrombectomy using Solitaire (Medtronic) or Trevo (Styker). Our primary outcome was FPE, defined as Thrombolysis in Cerebral Infarction score of 2c-3 achieved after a single thrombectomy attempt. We also performed retrospective analysis of same clot imaging characteristics of consecutive cases of MCA occlusion and its association with FPE matching the 3-dimensional in vitro experiments. Analysis was conducted using IBM SPSS Statistics Version 25 (IBM Corp., Armonk, NY). Chi-square tests and bivariate logistic regressions were the main statistical tests used in analysis. A p-value of less than .05 was considered to indicate statistical significance. Ninety-five confidence intervals (95% CI) were generated.

RESULTS

We compared 41 thrombectomy experiments performed using patient-specific 3-dimensional in vitro models with a retrospective cohort of 41 patients treated with SR thrombectomy. We found that in the in vitro cohort, higher MCA angulation was associated with a lower likelihood of FPE (odds ratio [OR] = 0.967, 95% CI = 0.944-0.991, p = .008). Meanwhile in the in vivo cohort, higher MCA angulation was associated with a higher likelihood of FPE (OR = 1.039, 95% CI = 1.003-1.077, p = .033). Neither clot length nor location of clot (M1 vs. M2) was associated with a difference in FPE rates in either cohort.

DISCUSSION

Comparison of SR thrombectomy performed during actual MCA occlusion cases versus patient-specific 3-dimensional replicas revealed MCA angulation as an independent predictor of procedure success or failure. However, the opposite direction of effect was observed between the two studied environments, indicating potential limitations of studying SR thrombectomy using 3-dimensional models of LVO.

Thrombolysis: Observations and numerical models

This perspective paper considers thrombolysis in the context of ischemic strokes, intending to build eventually a numerical model capable of simulating the thrombolytic treatment and predicting patient outcomes. Numerical modeling is a scientific methodology based on an abstraction of a system but requires understanding their spatio-temporal interactions. However, although important, the current knowledge on thrombolysis is fragmented in contributions from which it is difficult to obtain a complete picture of the process, especially in a clinically relevant setup. This paper discusses, from a general point of view, how to develop a numerical model to describe the evolution of a patient clot under the action of a thrombolytic drug. We will present critical, yet fundamental, open questions that have emerged during this elaboration and discuss original experimental observations that challenge some of our current knowledge of thrombolysis.

Aspects of ischemic stroke biomechanics derived using ex-vivo and in-vitro methods relating to mechanical thrombectomy

Establishing the underlying biomechanics of acute ischemic stroke (AIS) and its treatment is fundamental to developing more effective clinical treatments for one of society's most impactful diseases. Recent changes in AIS management, driven by clinical evidence of improved treatments, has already led to a rapid rate of innovation, which is likely to be sustained for many years to come. These unprecedented AIS triage and treatment innovations provide a great opportunity to better understand the disease. In this article we provide a perspective on the recreation of AIS in the laboratory to inform contemporary device design and procedural techniques in mechanical thrombectomy. Presentation of these findings, which have been used to solve the applied problem of designing mechanical thrombectomy devices, is intended to help inform the development of basic biomechanics solutions for AIS.

Measuring the effect of thrombosis, thrombus maturation and thrombolysis on clot mechanical properties in an in-vitro model

Changes in acute ischemic stroke thrombi structure and composition may result in significant differences in treatment responsiveness. Ischemic stroke patients are often treated with a thrombolytic agent to dissolve thrombi, however these patients may subsequently undergo mechanical thrombectomy to remove the occlusive clot. We set out to determine if rt-PA thrombolysis treatment of blood clots changes their mechanical properties, which in turn may impact mechanical thrombectomy. Using a design-of-experiment approach, ovine clot analogues were prepared with varying composition and further exposed to different levels of compaction force to simulate the effect of arterial blood pressure. Finally, clots were treated with three r-tPA doses for different durations. Clot mass and mechanical behaviour was analysed to assess changes due to (i) Platelet driven contraction (ii) Compaction force and (iii) Thrombolysis. Clots that were exposed to r-tPA for longer duration showed significant reduction in clot mass (p < 0.001). Exposure time to r-tPA (p < 0.001) was shown to be an independent predictor of lower clot stiffness. A decrease in energy dissipation ratio during mechanical compression was associated with longer exposure time in r-tPA (p = 0.001) and a higher platelet concentration ratio (p = 0.018). The dose of r-tPA was not a significant factor in reducing clot mass or changing mechanical properties of the clots. Fibrinolysis reduces clot stiffness which may explain increased distal clot migration observed in patients treated with r-tPA and should be considered as a potential clot modification factor before mechanical thrombectomy.

Quantitative 3D analysis of tissue damage in a rat model of microembolization

There is a discrepancy between successful recanalization and good clinical outcome after endovascular treatment (EVT) in acute ischemic stroke patients. During removal of a thrombus, a shower of microemboli may release and lodge to the distal circulation. The objective of this study was to determine the extent of damage on brain tissue caused by microemboli. In a rat model of microembolization, a mixture of microsphere (MS) sizes (15, 25 and 50 µm diameter) was injected via the left internal carotid artery. A 3D image of the left hemisphere was reconstructed and a point-pattern spatial analysis was applied based on G- and K-functions to unravel the spatial correlation between MS and the induced hypoxia or infarction. We show a spatial correlation between MS and hypoxia or infarction spreading up to a distance of 1000-1500 µm. These results imply that microemboli, which individually may not always be harmful, can interact and result in local areas of hypoxia or even infarction when lodged in large numbers.

A Thrombectomy Model Based on Ex Vivo Whole Human Brains

BACKGROUND AND PURPOSE

The persistent challenges in thrombectomy for large-vessel occlusion, such as suboptimal complete recanalization and first-pass effect imply an insufficient understanding of the artery-clot-device interaction. In this study, we present a thrombectomy model using fresh human brains, which can capture the artery-clot-device interaction through concurrent transmural and angiographic visualizations.

MATERIALS AND METHODS

Fresh nonfrozen whole adult human brains were collected and connected to a customized pump system tuned to deliver saline flow at a physiologic flow rate and pressure. Angiography was performed to verify the flow in the anterior-posterior and vertebrobasilar circulations and collaterals. Large-vessel occlusion was simulated by embolizing a radiopaque clot analog. Thrombectomy was tested, and the artery-clot-device interactions were recorded by transmural and angiographic videos.

RESULTS

Baseline cerebral angiography revealed excellent penetration of contrast in the anterior-posterior and vertebrobasilar circulations without notable arterial cutoffs and with robust collaterals. Small branches (<0.5 mm) and perforating arteries were consistently opacified with good patency. Three device passes were performed to achieve recanalization, with failure modes including elongation, fragmentation, and distal embolization.

CONCLUSIONS

This model enables concurrent transmural and angiographic analysis of artery-clot-device interaction in a human brain and provides critical insights into the action mechanism and failure modes of current and upcoming thrombectomy devices.

In silico trials for treatment of acute ischemic stroke: Design and implementation

An in silico trial simulates a disease and its corresponding therapies on a cohort of virtual patients to support the development and evaluation of medical devices, drugs, and treatment. In silico trials have the potential to refine, reduce cost, and partially replace current in vivo studies, namely clinical trials and animal testing. We present the design and implementation of an in silico trial for treatment of acute ischemic stroke. We propose an event-based modelling approach for the simulation of a disease and injury, where changes to the state of the system (the events) are assumed to be instantaneous. Using this approach we are able to combine a diverse set of models, spanning multiple time scales, to model acute ischemic stroke, treatment, and resulting brain tissue injury. The in silico trial is designed to be modular to aid development and reproducibility. It provides a comprehensive framework for application to any potential in silico trial. A statistical population model is used to generate cohorts of virtual patients. Patient functional outcomes are also predicted with a statistical model, using treatment and injury results and the patient's clinical parameters. We demonstrate the functionality of the event-based modelling approach and trial framework by running proof of concept in silico trials. The proof of concept trials simulate the same cohort of patients twice: once with successful treatment (successful recanalisation) and once with unsuccessful treatment (unsuccessful treatment). Ways to overcome some of the challenges and difficulties in setting up such an in silico trial are discussed, such as validation and computational limitations.

Preclinical testing platforms for mechanical thrombectomy in stroke: a review on phantoms, in-vivo animal, and cadaveric models

Preclinical testing platforms have been instrumental in the research and development of thrombectomy devices. However, there is no single model which fully captures the complexity of cerebrovascular anatomy, physiology, and the dynamic artery-clot-device interaction. This article provides a critical review of phantoms, in-vivo animal, and human cadaveric models used for thrombectomy testing and provides insights into the strengths and limitations of each platform. Articles published in the past 10 years that reported thrombectomy testing platforms were identified. Characteristics of each test platform, such as intracranial anatomy, artery tortuosity, vessel friction, flow conditions, device-vessel interaction, and visualization, were captured and benchmarked against human cerebral vessels involved in large-vessel occlusion stroke. Thrombectomy phantoms have been constructed from silicone, direct 3D-printed polymers, and glass. These phantoms represent oversimplified patient-specific cerebrovascular geometry but enable adequate visualization of devices and clots under appropriate flow conditions. They do not realistically mimic the artery-clot interaction. For the animal models, arteries from swine, canines, and rabbits have been reported. These models can reasonably replicate the artery-clot-device interaction and have the unique value of evaluating the safety of thrombectomy devices. However, the vasculature geometries are substantially less complex and flow conditions are different from human cerebral arteries. Cadaveric models are the most accurate vascular representations but with limited access and challenges in reproducibility of testing conditions. Multiple test platforms should be likely used for comprehensive evaluation of thrombectomy devices. Interpretation of the testing results should take into consideration platform-specific limitations.

Investigating the two regimes of fibrin clot lysis: an experimental and computational approach

It has been observed in vitro that complete clot lysis is generally preceded by a slow phase of lysis during which the degradation seems to be inefficient. However, this slow regime was merely noticed, but not yet quantitatively discussed. In our experiments, we observed that the lysis ubiquitously occurred in two distinct regimes, a slow and a fast lysis regime. We quantified extensively the duration of these regimes for a wide spectrum of experimental conditions and found that on average, the slow regime lasts longer than the fast one, meaning that during most of the process, the lysis is ineffective. We proposed a computational model in which the properties of the binding of the proteins change during the lysis: first, the biochemical reactions take place at the surface of the fibrin fibers, then in the bulk, resulting in the observed fast lysis regime. This simple hypothesis appeared to be sufficient to reproduce with a great accuracy the lysis profiles obtained experimentally.

Impact of the Internal Carotid Artery Morphology on in silico Stent-Retriever Thrombectomy Outcome

The aim of this work is to propose a methodology for identifying relationships between morphological features of the cerebral vasculature and the outcome of in silico simulations of thrombectomy, the mechanical treatment for acute ischemic stroke. Fourteen patient-specific cerebral vasculature segmentations were collected and used for geometric characterization of the intracranial arteries mostly affected by large vessel occlusions, i.e., internal carotid artery (ICA), middle cerebral artery (MCA) and anterior cerebral artery (ACA). First, a set of global parameters was created, including the geometrical information commonly provided in the clinical context, namely the total length, the average diameter and the tortuosity (length over head-tail distance) of the intracranial ICA. Then, a more exhaustive geometrical analysis was performed to collect a set of local parameters. A total of 27 parameters was measured from each patient-specific vascular configuration. Fourteen virtual thrombectomy simulations were performed with a blood clot with the same length and composition placed in the middle of the MCA. The model of TREVO ProVue stent-retriever was used for all the simulations. Results from simulations produced five unsuccessful outcomes, i.e., the clot was not removed from the vessels. The geometric parameters of the successful and unsuccessful simulations were compared to find relations between the vascular geometry and the outcome. None of the global parameters alone or combined proved able to discriminate between positive and negative outcome, while a combination of local parameters allowed to correctly identify the successful from the unsuccessful simulations. Although these results are limited by the number of patients considered, this study indicates a promising methodology to relate patient-specific geometry to virtual thrombectomy outcome, which might eventually guide decision making in the treatment of acute ischemic stroke.

Applicability analysis to evaluate credibility of an in silico thrombectomy procedure

Intra-arterial thrombectomy is a minimally invasive procedure in which an obstructing thrombus (clot) is removed using a minimally-invasive device: a stent-retriever. The stent-retriever is first deployed, and then the thrombus is removed during stent-retriever retraction. This procedure can be simulated using a detailed computational model. However, to be useful for an in silico trial in a clinical setting, model credibility should be demonstrated. The aim of this work is to apply a credibility process for the validation phases to the thrombectomy procedure in order to deem it credible for use in an in silico trial. Validation evidence is proposed for the identified context of use and then used to build credibility to the numerical model. Applicability of the proposed model is justified and assessed using a rigorous step-by-step method based on the ASME V&V40 protocol.

The first virtual patient-specific thrombectomy procedure

Treatment of acute ischemic stroke has been recently improved with the introduction of endovascular mechanical thrombectomy, a minimally invasive procedure able to remove a clot using aspiration devices and/or stent-retrievers. Despite the promising and encouraging results, improvements to the procedure and to the stent design are the focus of the recent efforts. Computational studies can pave the road to these improvements, providing their ability to describe and accurately reproduce a real procedure. A patient with ischemic stroke due to intracranial large vessel occlusion was selected and after the creation of the cerebral vasculature from computed tomography images and a histologic analysis to determine the clot composition, the entire thrombectomy procedure was virtually replicated. As in the real situation, the computational replica showed that two attempts were necessary to remove the clot, as a result of the position of the stent retriever with respect to the clot. Furthermore, the results indicated that clot fragmentation did not occur as the deformations were mainly in a compressive state without the possibility for clot cracks to propagate. The accurate representation of the procedure can be used as an important step for operative optimization planning and future improvements of stent designs.

On the Sensitivity Analysis of Porous Finite Element Models for Cerebral Perfusion Estimation

Computational physiological models are promising tools to enhance the design of clinical trials and to assist in decision making. Organ-scale haemodynamic models are gaining popularity to evaluate perfusion in a virtual environment both in healthy and diseased patients. Recently, the principles of verification, validation, and uncertainty quantification of such physiological models have been laid down to ensure safe applications of engineering software in the medical device industry. The present study sets out to establish guidelines for the usage of a three-dimensional steady state porous cerebral perfusion model of the human brain following principles detailed in the verification and validation (V&V 40) standard of the American Society of Mechanical Engineers. The model relies on the finite element method and has been developed specifically to estimate how brain perfusion is altered in ischaemic stroke patients before, during, and after treatments. Simulations are compared with exact analytical solutions and a thorough sensitivity analysis is presented covering every numerical and physiological model parameter. The results suggest that such porous models can approximate blood pressure and perfusion distributions reliably even on a coarse grid with first order elements. On the other hand, higher order elements are essential to mitigate errors in volumetric blood flow rate estimation through cortical surface regions. Matching the volumetric flow rate corresponding to major cerebral arteries is identified as a validation milestone. It is found that inlet velocity boundary conditions are hard to obtain and that constant pressure inlet boundary conditions are feasible alternatives. A one-dimensional model is presented which can serve as a computationally inexpensive replacement of the three-dimensional brain model to ease parameter optimisation, sensitivity analyses and uncertainty quantification. The findings of the present study can be generalised to organ-scale porous perfusion models. The results increase the applicability of computational tools regarding treatment development for stroke and other cerebrovascular conditions.

Modelling the impact of clot fragmentation on the microcirculation after thrombectomy

Many ischaemic stroke patients who have a mechanical removal of their clot (thrombectomy) do not get reperfusion of tissue despite the thrombus being removed. One hypothesis for this 'no-reperfusion' phenomenon is micro-emboli fragmenting off the large clot during thrombectomy and occluding smaller blood vessels downstream of the clot location. This is impossible to observe in-vivo and so we here develop an in-silico model based on in-vitro experiments to model the effect of micro-emboli on brain tissue. Through in-vitro experiments we obtain, under a variety of clot consistencies and thrombectomy techniques, micro-emboli distributions post-thrombectomy. Blood flow through the microcirculation is modelled for statistically accurate voxels of brain microvasculature including penetrating arterioles and capillary beds. A novel micro-emboli algorithm, informed by the experimental data, is used to simulate the impact of micro-emboli successively entering the penetrating arterioles and the capillary bed. Scaled-up blood flow parameters-permeability and coupling coefficients-are calculated under various conditions. We find that capillary beds are more susceptible to occlusions than the penetrating arterioles with a 4x greater drop in permeability per volume of vessel occluded. Individual microvascular geometries determine robustness to micro-emboli. Hard clot fragmentation leads to larger micro-emboli and larger drops in blood flow for a given number of micro-emboli. Thrombectomy technique has a large impact on clot fragmentation and hence occlusions in the microvasculature. As such, in-silico modelling of mechanical thrombectomy predicts that clot specific factors, interventional technique, and microvascular geometry strongly influence reperfusion of the brain. Micro-emboli are likely contributory to the phenomenon of no-reperfusion following successful removal of a major clot.

Coupling one-dimensional arterial blood flow to three-dimensional tissue perfusion models for in silico trials of acute ischaemic stroke

An acute ischaemic stroke is due to the sudden blockage of an intracranial blood vessel by an embolized thrombus. In the context of setting up in silico trials for the treatment of acute ischaemic stroke, the effect of a stroke on perfusion and metabolism of brain tissue should be modelled to predict final infarcted brain tissue. This requires coupling of blood flow and tissue perfusion models. A one-dimensional intracranial blood flow model and a method to couple this to a brain tissue perfusion model for patient-specific simulations is presented. Image-based patient-specific data on the anatomy of the circle of Willis are combined with literature data and models for vessel anatomy not visible in the images, to create an extended model for each patient from the larger vessels down to the pial surface. The coupling between arterial blood flow and tissue perfusion occurs at the pial surface through the estimation of perfusion territories. The coupling method is able to accurately estimate perfusion territories. Finally, we argue that blood flow can be approximated as steady-state flow at the interface between arterial blood flow and tissue perfusion to reduce the cost of organ-scale simulations.

A porous circulation model of the human brain for in silico clinical trials in ischaemic stroke

The advancement of ischaemic stroke treatment relies on resource-intensive experiments and clinical trials. In order to improve ischaemic stroke treatments, such as thrombolysis and thrombectomy, we target the development of computational tools for in silico trials which can partially replace these animal and human experiments with fast simulations. This study proposes a model that will serve as part of a predictive unit within an in silico clinical trial estimating patient outcome as a function of treatment. In particular, the present work aims at the development and evaluation of an organ-scale microcirculation model of the human brain for perfusion prediction. The model relies on a three-compartment porous continuum approach. Firstly, a fast and robust method is established to compute the anisotropic permeability tensors representing arterioles and venules. Secondly, vessel encoded arterial spin labelling magnetic resonance imaging and clustering are employed to create an anatomically accurate mapping between the microcirculation and large arteries by identifying superficial perfusion territories. Thirdly, the parameter space of the problem is reduced by analysing the governing equations and experimental data. Fourthly, a parameter optimization is conducted. Finally, simulations are performed with the tuned model to obtain perfusion maps corresponding to an open and an occluded (ischaemic stroke) scenario. The perfusion map in the occluded vessel scenario shows promising qualitative agreement with computed tomography images of a patient with ischaemic stroke caused by large vessel occlusion. The results highlight that in the case of vessel occlusion (i) identifying perfusion territories is essential to capture the location and extent of underperfused regions and (ii) anisotropic permeability tensors are required to give quantitatively realistic estimation of perfusion change. In the future, the model will be thoroughly validated against experiments.