Applicability assessment of a stent-retriever thrombectomy finite-element model

A hyper-viscoelastic uniaxial characterization of collagenous embolus analogs in acute ischemic stroke

PURPOSE

Acute ischemic stroke is a leading cause of death and morbidity worldwide. Despite advances in medical technology, nearly 30% of strokes result in incomplete vessel recanalization. Recent studies have demonstrated that clot composition correlates with success rates of mechanical thrombectomy procedures. To understand clot behavior during thrombectomy, which exerts considerable strains on thrombi, in vitro studies must characterize the rate-dependent high-strain behavior of embolus analogs (EAs) with different formation conditions, which can be used to fit models of hyper-viscoelasticity.

METHODS

In this study, the effect of collagen infiltration as a carotid-induced collagen-rich thrombosis surrogate is considered as a contributor to embolus analog high-strain stiffness, when compared to 40% hematocrit EAs.

RESULTS

EA high-strain stiffnesses, characterized on a uniaxial load frame, increase by an order of magnitude for collagenous clot analogs. Chandler loop analogs show high-strain stiffnesses and clot compositions commensurate with previous reports of stroke patient clots, and collagenous clots show significant increase in stiffness when compared to stroke patient clots. Finally, hyper-viscoelastic curve fitting demonstrates the asymmetry between tension and compression. Nonlinear, rate-dependent models that consider clot-stiffening behavior match the high strain stiffness of clots fairly well. Furthermore, we demonstrate that the stability of the elastic energy needs to be considered to obtain optimal curve fits for high-strain, rate dependent data.

CONCLUSION

This study provides a framework for the development of dynamically formed EAs that mimic the mechanical and structural properties of in vivo clots and provides parameters for numerical simulation of clot behavior with hyper-viscoelastic models.

Flow reduction due to arterial catheterization during stroke treatment - A computational study using a distributed compartment model

The effectiveness of various stroke treatments depends on the anatomical variability of the cerebral vasculature, particularly the collateral blood vessel network. Collaterals at the level of the Circle of Willis and distal collaterals, such as the leptomeningeal arteries, serve as alternative avenues of flow when the primary pathway is obstructed during an ischemic stroke. Stroke treatment typically involves catheterization of the primary pathway, and the potential risk of further flow reduction to the affected brain area during this treatment has not been previously investigated. To address this clinical question, we derived the lumped parameters for catheterized blood vessels and implemented a corresponding distributed compartment (0D) model. This 0D model was validated against an experimental model and benchmark test cases solved using a 1D model. Additionally, we compared various off-center catheter trajectories modeled using a 3D solver to this 0D model. The differences between them were minimal, validating the simplifying assumption of the central catheter placement in the 0D model. The 0D model was then used to simulate blood flows in realistic cerebral arterial networks with different collateralization characteristics. Ischemic strokes were modeled by occlusion of the M1 segment of the middle cerebral artery in these networks. Catheters of different diameters were inserted up to the obstructed segment and flow alterations in the network were calculated. Results showed up to 45% maximum blood flow reduction in the affected brain region. These findings suggest that catheterization during stroke treatment may have a further detrimental effect for some patients with poor collateralization.

Energy barriers govern catheter herniation during endovascular procedures: a 2.5D vascular flow model analysis

Endovascular procedures rely on navigating guidewires, catheters and other devices through tortuous vasculature to treat disease. A critical challenge in these procedures is catheter herniation, in which the device deviates from its intended path, often irrecoverably. To elucidate the mechanics of herniation, we developed a physical flow model of the aortic arch that enables direct measurement of device curvature during experimentally simulated neuroendovascular procedures conducted from an upper arterial access. Combined with measurements of initial, unstressed device shapes and flexural rigidities, the method enables the experimental estimation of the device bending energies during these simulated procedures. Characteristic energy profiles revealed distinct stages in both herniation and successful navigation, governed by the interplay between device properties and vascular anatomy. A deterministic progression from successful navigation to herniation was identified, with catheter systems following paths determined by measurable energy barriers. Increasing guidewire stiffness or decreasing catheter stiffness reduced the energy barrier for successful navigation while increasing that for herniation. This framework enables the prediction of endovascular herniation risk and offers unique insight into improved device design and clinical decision-making.

Collateral flow and pulsatility during large vessel occlusions: insights from a quantitative in vitro study

Acute ischemic stroke caused by large vessel occlusions is being increasingly treated with neurovascular interventions. The hemodynamics within the collateral system of the circle of Willis (CoW) hemodynamics play a fundamental role in therapy success. However, transient in vivo data on pathological collateral flow during large vessel occlusions are not available. Moreover, there are no flow models that accurately simulate the hemodynamic conditions in the CoW during large vessel occlusions. We used a circulatory loop to generate highly reproducible cerebrovascular-like flows and pressures and used non-invasive flow visualization and high-resolution flow and pressure measurements to acquire detailed, time-dependent hemodynamics inside an anatomical phantom of the CoW. After calibrating a physiological reference case, we induced occlusions in the 1. middle cerebral artery, 2. terminal carotid artery, and 3. basilar artery; and measured the left posterior communicating artery flow. Mean arterial pressure and pulse pressure remained unchanged in the different occlusion cases compared to the physiological reference case, while total cerebral flow decreased by up to 19%. In all three occlusion cases, reversed flow was found in the left posterior communicating artery compared to the reference case with different flow magnitudes and pulsatility index values. The experimental results were compared with clinical findings, demonstrating the capability of this realistic cerebrovascular flow setup. This novel cerebrovascular flow setup opens the possibility for investigating different topics of neurovascular interventions under various clinical conditions in controlled preclinical laboratory studies.

A one-dimensional computational model for blood flow in an elastic blood vessel with a rigid catheter

Strokes are one of the leading causes of death in the United States. Stroke treatment involves removal or dissolution of the obstruction (usually a clot) in the blocked artery by catheter insertion. A computer simulation to systematically plan such patient-specific treatments needs a network of about 105 blood vessels including collaterals. The existing computational fluid dynamic (CFD) solvers are not employed for stroke treatment planning as they are incapable of providing solutions for such big arterial trees in a reasonable amount of time. This work presents a novel one-dimensional mathematical formulation for blood flow modeling in an elastic blood vessel with a centrally placed rigid catheter. The governing equations are first-order hyperbolic partial differential equations, and the hypergeometric function needs to be computed to obtain the characteristic system of these hyperbolic equations. We employed the Discontinuous Galerkin method to solve the hyperbolic system and validated the implementation by comparing it against a well-established 3D CFD solver using idealized vessels and a realistic truncated arterial network. The results showed clinically insignificant differences in steady flow cases, with overall variations between 1D and 3D models remaining below 10%. Additionally, the solver accurately captured wave reflection phenomena at domain discontinuities in unsteady cases. A primary advantage of this model over 3D solvers is its ease in obtaining a discretized geometry of complex vasculatures with multiple arterial branches. Thus, the 1D computational model offers good accuracy and applicability in simulating complex vasculatures, demonstrating promising potential for investigating patient-specific endovascular interventions in strokes.

Development of a clot-adhesive coating to improve the performance of thrombectomy devices

BACKGROUND

The first-pass complete recanalization by mechanical thrombectomy (MT) for the treatment of stroke remains limited due to the poor integration of the clot within current devices. Aspiration can help retrieval of the main clot but fails to prevent secondary embolism in the distal arterial territory. The dense meshes of extracellular DNA, recently described in stroke-related clots, might serve as an anchoring platform for MT devices. We aimed to evaluate the potential of a DNA-reacting surface to aid the retention of both the main clot and small fragments within the thrombectomy device to improve the potential of MT procedures.

METHODS

Device-suitable alloy samples were coated with 15 different compounds and put in contact with extracellular DNA or with human peripheral whole blood, to compare their binding to DNA versus blood elements in vitro. Clinical-grade MT devices were coated with two selected compounds and evaluated in functional bench tests to study clot retrieval efficacy and quantify distal emboli using an M1 occlusion model.

RESULTS

Binding properties of samples coated with all compounds were increased for DNA (≈3-fold) and decreased (≈5-fold) for blood elements, as compared with the bare alloy samples in vitro. Functional testing showed that surface modification with DNA-binding compounds improved clot retrieval and significantly reduced distal emboli during experimental MT of large vessel occlusion in a three-dimensional model.

CONCLUSION

Our results suggest that clot retrieval devices coated with DNA-binding compounds can considerably improve the outcome of the MT procedures in stroke patients.

In vitro models for the experimental evaluation of mechanical thrombectomy devices in acute ischemic stroke

Mechanical thrombectomy has become an important method for the treatment of acute ischemic stroke for large vessel occlusions. The current hotspots of mechanical thrombectomy are optimizing the treatment methods, improving the recanalization rate and reducing complications. The in vitro model has become a common and convenient method for mechanical thrombectomy research. This review summarizes the in vitro model in the following aspects: the preparation of clot analogues; the experimental platform; the application of the in vitro model in the testing of thrombectomy devices; and the advantages, limitations and future trends of the in vitro experimental model. This review describes the characteristics and applications of the in vitro experimental model with the hope that the in vitro experimental model will be further improved and play a more effective role in the study of mechanical thrombectomy.

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.

Recommendations for finite element modelling of nickel-titanium stents-Verification and validation activities

The objective of this study is to present a credibility assessment of finite element modelling of self-expanding nickel-titanium (Ni-Ti) stents through verification and validation (VV) activities, as set out in the ASME VV-40 standard. As part of the study, the role of calculation verification, model input sensitivity, and model validation is examined across three different application contexts (radial compression, stent deployment in a vessel, fatigue estimation). A commercially available self-expanding Ni-Ti stent was modelled, and calculation verification activities addressed the effects of mesh density, element integration and stable time increment on different quantities of interests, for each context of use considered. Sensitivity analysis of the geometrical and material input parameters and validation of deployment configuration with in vitro comparators were investigated. Results showed similar trends for global and local outputs across the contexts of use in response to the selection of discretization parameters, although with varying sensitivities. Mesh discretisation showed substantial variability for less than 4 × 4 element density across the strut cross-section in radial compression and deployment cases, while a finer grid was deemed necessary in fatigue estimation for reliable predictions of strain/stress. Element formulation also led to substantial variation depending on the chosen integration options. Furthermore, for explicit analyses, model results were highly sensitive to the chosen target time increment (e.g., mass scaling parameters), irrespective of whether quasistatic conditions were ensured (ratios of kinetic and internal energies below 5%). The higher variability was found for fatigue life simulation, with the estimation of fatigue safety factor varying up to an order of magnitude depending on the selection of discretization parameters. Model input sensitivity analysis highlighted that the predictions of outputs such as radial force and stresses showed relatively low sensitivity to Ni-Ti material parameters, which suggests that the calibration approaches used in the literature to date appear reasonable, but a higher sensitivity to stent geometry, namely strut thickness and width, was found. In contrast, the prediction of vessel diameter following deployment was least sensitive to numerical parameters, and its validation with in vitro comparators offered a simple and accurate (error ~ 1-2%) method when predicting diameter gain, and lumen area, provided that the material of the vessel is appropriately characterized and modelled.

Microcatheter tracking in thrombectomy procedures: A finite-element simulation study

BACKGROUND AND OBJECTIVE

Mechanical thrombectomy is a minimally invasive procedure that aims at removing the occluding thrombus from the vasculature of acute ischemic stroke patients. Thrombectomy success and failure can be studied using in-silico thrombectomy models. Such models require realistic modeling steps to be effective. We here present a new approach to model microcatheter tracking during thrombectomy.

METHODS

For 3 patient-specific vessel geometries, we performed finite-element simulations of the microcatheter tracking (1) following the vessel centerline (centerline method) and (2) as a one-step insertion simulation, where the microcatheter tip was advanced along the vessel centerline while its body was free to interact with the vessel wall (tip-dragging method). Qualitative validation of the two tracking methods was performed with the patient's digital subtraction angiography (DSA) images. In addition, we compared simulated thrombectomy outcomes (successful vs unsuccessful thrombus retrieval) and maximum principal stresses on the thrombus between the centerline and tip-dragging method.

RESULTS

Qualitative comparison with the DSA images showed that the tip-dragging method more realistically resembles the patient-specific microcatheter-tracking scenario, where the microcatheter approaches the vessel walls. Although the simulated thrombectomy outcomes were similar in terms of thrombus retrieval, the thrombus stress fields (and the associated fragmentation of the thrombus) were strongly different between the two methods, with local differences in the maximum principal stress curves up to 84%.

CONCLUSIONS

Microcatheter positioning with respect to the vessel affects the stress fields of the thrombus during retrieval, and therefore, may influence thrombus fragmentation and retrieval in-silico thrombectomy.

Investigation of Cerebral Hemodynamics During Endovascular Aspiration: Development of an Experimental and Numerical Setup

PURPOSE

Acute ischemic stroke is a life-threatening emergency caused by an occlusion of a cerebral artery through a blood clot. Aspiration thrombectomy is an endovascular therapy for the removal of vessel occlusions. However, open questions regarding the hemodynamics during the intervention remain, motivating investigations of blood flow within cerebral arteries. In this study, we present a combined experimental and numerical approach to analyze hemodynamics during endovascular aspiration.

METHODS

We have developed an in vitro setup for investigations of hemodynamic changes during endovascular aspiration within a compliant model of patient-specific cerebral arteries. Pressures, flows, and locally resolved velocities were obtained. In addition, we established a computational fluid dynamics (CFD) model and compared the simulations during physiological conditions and in two aspiration scenarios with different occlusions.

RESULTS

Flow redistribution within cerebral arteries after ischemic stroke is strongly dependent on the severity of the occlusion and on the volume flow extracted by endovascular aspiration. Numerical simulations exhibit an excellent correlation of R = 0.92 for flow rates and a good correlation of R = 0.73 for pressures. Further on, the local velocity field inside the basilar artery had a good agreement between CFD model and particle image velocimetry (PIV) data.

CONCLUSION

The presented setup allows for in vitro investigations of artery occlusions and endovascular aspiration techniques on arbitrary patient-specific cerebrovascular anatomies. The in silico model provides consistent predictions of flows and pressures in several aspiration scenarios.

A low dimensional surrogate model for a fast estimation of strain in the thrombus during a thrombectomy procedure

BACKGROUND

Intra-arterial thrombectomy is the main treatment for acute ischemic stroke due to large vessel occlusions and can consist in mechanically removing the thrombus with a stent-retriever. A cause of failure of the procedure is the fragmentation of the thrombus and formation of micro-emboli, difficult to remove. This work proposes a methodology for the creation of a low-dimensional surrogate model of the mechanical thrombectomy procedure, trained on realizations from high-fidelity simulations, able to estimate the evolution of the maximum first principal strain in the thrombus.

METHOD

A parametric finite-element model was created, composed of a tapered vessel, a thrombus, a stent-retriever and a catheter. A design of experiments was conducted to sample 100 combinations of the model parameters and the corresponding thrombectomy simulations were run and post-processed to extract the maximum first principal strain in the thrombus during the procedure. Then, a surrogate model was built with a combination of principal component analysis and Kriging.

RESULTS

The surrogate model was chosen after a sensitivity analysis on the number of principal components and was tested with 10 additional cases. The model provided predictions of the strain curves with correlation above 0.9 and a maximum error of 28%, with an error below 20% in 60% of the test cases.

CONCLUSIONS

The surrogate model provides nearly instantaneous estimates and constitutes a valuable tool for evaluating the risk of thrombus rupture during pre-operative planning for the treatment of acute ischemic stroke.

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.

Combined stent-retriever and aspiration intra-arterial thrombectomy performance for fragmentable blood clots: A proof-of-concept computational study

Mechanical thrombectomy (MT) treatment of acute ischemic stroke (AIS) patients typically involves use of stent retrievers or aspiration catheters alone or in combination. For in silico trials of AIS patients, it is crucial to incorporate the possibility of thrombus fragmentation during the intervention. This study focuses on two aspects of the thrombectomy simulation: i) Thrombus fragmentation on the basis of a failure model calibrated with experimental tests on clot analogs; ii) the combined stent-retriever and aspiration catheter MT procedure is modeled by adding both the proximal balloon guide catheter and the distal access catheter. The adopted failure criterion is based on maximum principal stress threshold value. If elements of the thrombus exceed this criterion during the retrieval simulation, then they are deleted from the calculation. Comparison with in-vitro tests indicates that the simulation correctly reproduces the procedures predicting thrombus fragmentation in the case of red blood cells rich thrombi, whereas non-fragmentation is predicted for fibrin-rich thrombi. Modeling of balloon guide catheter prevents clot fragments' embolization to further distal territories during MT procedure.

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.

Self-expandable stent for thrombus removal modeling: Solid or beam finite elements?

BACKGROUND

The performance of self-expandable stents is being increasingly studied by means of finite-element analysis. As for peripheral stents, transcatheter valves and stent-grafts, there are numerous computational studies for setting up a proper model, this information is missing for stent-retrievers used in the procedure of thrombus removal in cerebral arteries. It is well known that the selection of the appropriate finite-element dimensions (topology) and formulations (typology) is a fundamental step to set up accurate and reliable computational simulations. In this context, a thorough verification analysis is here proposed, aimed at investigating how the different element typologies and topologies - available to model a stent-retriever - affect simulation results.

METHOD

Hexahedral and beam element formulations were analyzed first individually by virtually replicating a crimping test on the device, and then by replicating the thrombectomy procedure aiming at removing a thrombus from a cerebral vessel. In particular, three discretization refinements for each element type and different element formulations including both full and reduced integration were investigated and compared in terms of the resultant radial force of the stent and the stress field generated in the thrombus.

RESULTS

The sensitivity analysis on the element formulation performed with the crimping simulations allowed the identification of the optimal setting for each element family. Both setting lead to similar results in terms of stent performance in the virtual thrombectomy and should be used in future studies simulating the mechanical thrombectomy with stent-retrievers.

CONCLUSIONS

The carried out virtual thrombectomy procedures confirmed that the beam element formulation results were sufficiently accurate to model the radial force and the performance of the stent-retriever during the procedure. For different self-expandable stents, hexahedral formulation could be essential in stress analysis.

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.

Validation and Verification of High-Fidelity Simulations of Thoracic Stent-Graft Implantation

Thoracic Endovascular Aortic Repair (TEVAR) is the preferred treatment option for thoracic aortic pathologies and consists of inserting a self-expandable stent-graft into the pathological region to restore the lumen. Computational models play a significant role in procedural planning and must be reliable. For this reason, in this work, high-fidelity Finite Element (FE) simulations are developed to model thoracic stent-grafts. Experimental crimp/release tests are performed to calibrate stent-grafts material parameters. Stent pre-stress is included in the stent-graft model. A new methodology for replicating device insertion and deployment with explicit FE simulations is proposed. To validate this simulation, the stent-graft is experimentally released into a 3D rigid aortic phantom with physiological anatomy and inspected in a computed tomography (CT) scan at different time points during deployment with an ad-hoc set-up. A verification analysis of the adopted modeling features compared to the literature is performed. With the proposed methodology the error with respect to the CT is on average 0.92 ± 0.64%, while it is higher when literature models are adopted (on average 4.77 ± 1.83%). The presented FE tool is versatile and customizable for different commercial devices and applicable to patient-specific analyses.

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.

Realistic computer modelling of stent retriever thrombectomy: a hybrid finite-element analysis-smoothed particle hydrodynamics model

Stent retriever thrombectomy is a pre-eminent treatment modality for large vessel ischaemic stroke. Simulation of thrombectomy could help understand stent and clot mechanics in failed cases and provide a digital testbed for the development of new, safer devices. Here, we present a novel, in silico thrombectomy method using a hybrid finite-element analysis (FEA) and smoothed particle hydrodynamics (SPH). Inspired by its biological structure and components, the blood clot was modelled with the hybrid FEA-SPH method. The Solitaire self-expanding stent was parametrically reconstructed from micro-CT imaging and was modelled as three-dimensional finite beam elements. Our simulation encompassed all steps of mechanical thrombectomy, including stent packaging, delivery and self-expansion into the clot, and clot extraction. To test the feasibility of our method, we simulated clot extraction in simple straight vessels. This was compared against in vitro thrombectomies using the same stent, vessel geometry, and clot size and composition. Comparisons with benchtop tests indicated that our model was able to accurately simulate clot deflection and penetration of stent wires into the clot, the relative movement of the clot and stent during extraction, and clot fragmentation/embolus formation. In this study, we demonstrated that coupling FEA and SPH techniques could realistically model stent retriever thrombectomy.

A new constitutive model for permanent deformation of blood clots with application to simulation of aspiration thrombectomy

As a first line option in the treatment of acute ischemic stroke (AIS), direct aspiration is a fast and effective technique with promising outcomes. In silico models are widely used for design and preclinical assessment of new developed devices and therapeutic methods. Accurate modelling of the mechanical behaviour of blood clot is a key factor in the design and simulation of aspiration devices. In this study we develop a new constitutive model which incorporates the unrecoverable plastic deformation of clots. The model is developed based on the deformation-induced microstructural changes in fibrin network, including the formation and dissociation of the cross-links between fibrin fibres. The model is calibrated using previously reported experimentally measured permanent clot deformation following uniaxial stretching. The calibrated plasticity model is then used to simulate aspiration thrombectomy. Results reveal that inclusion of permanent plastic deformation results in ∼ 15 % increase in clot aspiration length at an applied aspiration pressure of 100 mmHg. The constitutive law developed in this study provides a basis for improved design and evaluation of novel aspiration catheters leading to increased first-pass revascularization rate.

Modelling Combined Intravenous Thrombolysis and Mechanical Thrombectomy in Acute Ischaemic Stroke: Understanding the Relationship between Stent Retriever Configuration and Clot Lysis Mechanisms

Background: Combined intravenous thrombolysis and mechanical thrombectomy (IVT-MT) is a common treatment in acute ischaemic stroke, however the interaction between IVT and MT from a physiological standpoint is poorly understood. In this pilot study, we conduct numerical simulations of combined IVT-MT with various idealised stent retriever configurations to evaluate performance in terms of complete recanalisation times and lysis patterns. Methods: A 3D patient-specific geometry of a terminal internal carotid artery with anterior and middle cerebral arteries is reconstructed, and a thrombus is artificially implanted in the MCA branch. Various idealised stent retriever configurations are implemented by varying stent diameter and stent placement, and a configuration without a stent retriever provides a baseline for comparison. A previously validated multi-level model of thrombolysis is used, which incorporates blood flow, drug transport, and fibrinolytic reactions within a fibrin thrombus. Results: Fastest total recanalisation was achieved in the thrombus without a stent retriever, with lysis times increasing with stent retriever diameter. Two mechanisms of clot lysis were established: axial and radial permeation. Axial permeation from the clot front was the primary mechanism of lysis in all configurations, as it facilitated increased protein binding with fibrin fibres. Introducing a stent retriever channel allowed for radial permeation, which occurred at the fluid-thrombus interface, although lysis was much slower in the radial direction because of weaker secondary velocities. Conclusions: Numerical models can be used to better understand the complex physiological relationship between IVT and MT. Two different mechanisms of lysis were established, providing a basis towards improving the efficacy of combined treatments.

Simulation of stent retriever thrombectomy in acute ischemic stroke by finite element analysis

Stent retriever thrombectomy has become an effective method for treating acute ischemic stroke. Successful recanalization cannot be achieved for all patients. The outcome of mechanical thrombectomy may be associated with some mechanical factors, requiring efficient tools that are able to assess the interaction between stent retrievers and clots. In this study, a simulation by finite element analysis was developed to evaluate the outcome of stent retriever thrombectomy for clots of different sizes and frictional properties. The deployment and retrieval of the stent retriever were conducted in a middle cerebral artery model. The recanalization, deformation of the clot and stent retriever and the stress induced in the clot were assessed. The results showed that higher friction could lead to failed recanalization and increased stress in the clot. The simulation method can be used to characterize the mechanical behaviour of stent retrievers and clots, offering a potential tool for the optimization of device design and the selection of surgical strategies.

Finite element analysis of blood clots based on the nonlinear visco-hyperelastic model

Mechanical thrombectomy has become the standard treatment for patients with an acute ischemic stroke. In this approach, to remove blood clots, mechanical force is applied using thrombectomy devices, in which the interaction between the clot and the device could significantly affect the clot retrieval performance. It is expected that the finite element method (FEM) could visualize the mechanical interaction by the visualization of the stress transmission from the device to the clot. This research was aimed at verifying the constitutive theory by implementing FEM based on the visco-hyperelastic theory, using a three-dimensional clot model. We used the visco-hyperelastic FEM to reproduce the mechanical behavior of blood clots, as observed in experiments. This study is focused on the mechanical responses of clots under tensile loading and unloading because in mechanical thrombectomy, elongation is assumed to occur locally on the clots during the retrieval process. Several types of cylindrical clots were created by changing the fibrinogen dose. Tensile testing revealed that the stiffness (E_0.45-value) of clots with fibrinogen could be more than three times higher than that of clots without fibrinogen. It was also found that the stiffness was not proportional to the fibrinogen dose. By fitting to the theoretical curve, it was revealed that the Mooney-Rivlin model could reproduce the hyperelastic characteristics of clots well. From the stress-relaxation data, the three-chain Maxwell model could accurately fit the experimental viscoelastic data. FEM, taking the theoretical models into account, was then carried out, and the results matched well with the experimental visco-hyperelastic characteristics of clots under tensile load, reproducing the mechanical hysteresis during unloading, the stress dependence on the strain rate, and the time-dependent stress decrease in the stress-relaxation test.

In vitro and in silico modeling of endovascular stroke treatments for acute ischemic stroke

Acute ischemic stroke occurs when a thrombus obstructs a cerebral artery, leading to sub-optimal blood perfusion to brain tissue. A recently developed, preventive treatment is the endovascular stroke treatment (EVT), which is a minimally invasive procedure, involving the use of stent-retrievers and/or aspiration catheters. Despite its increasing use, many critical factors of EVT are not well understood. In this respect, in vitro, and in silico studies have the great potential to help us deepen our understanding of the procedure, perform further device and procedural optimization, and help in clinical training. This review paper provides an overview of the previous in vitro and in silico evaluations of EVT treatments, with a special emphasis on the four main aspects of the adopted experimental and numerical set-ups: vessel, thrombus, device, and procedural settings.

A clinical perspective on endovascular stroke treatment biomechanics

Acute ischemic stroke (AIS) is caused by blockage of an arterial blood vessel in the brain by a thrombus, which interrupts oxygen supply to the brain parenchyma. The goal of endovascular stroke treatment (mechanical thrombectomy) is to restore blood flow as quickly and completely as possible. There are numerous factors that influence endovascular treatment success. They can be broadly grouped into a) factors related to blood vessels, b) factors related to the thrombus, c) factors related to endovascular treatment technique and tools and d) operator-related factors. While blood vessel and tgthro thrombus-related factors are mostly non-modifiable in the acute setting, operator and technique-related factors can be modified, and extensive research is currently being done to investigate the complex interplay of all these variables, and to optimize the modifiable factors to the maximum possible extent. In this review, we will describe these factors and how they interact with each other in detail, and outline some of their practical implications. We will conclude with a short summary and outlook on future directions for optimizing endovascular treatment success.

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.

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.

Blood clot fracture properties are dependent on red blood cell and fibrin content

Thrombus fragmentation during endovascular stroke treatment, such as mechanical thrombectomy, leads to downstream emboli, resulting in poor clinical outcomes. Clinical studies suggest that fragmentation risk is dependent on clot composition. This current study presents the first experimental characterization of the composition-dependent fracture properties of blood clots, in addition to the development of a predictive model for blood clot fragmentation. A bespoke experimental test-rig and compact tension specimen fabrication has been developed to measure fracture toughness of thrombus material. Fracture tests are performed on three physiologically relevant clot compositions: a high-fibrin clot made from a 5% haematocrit (H) blood mixture, a medium-fibrin clot made form a 20% H blood mixture, a low-fibrin clot made from a 40% H blood mixture. Fracture toughness is observed to significantly increase with increasing fibrin content, i.e. red blood cell-rich clots are more prone to tear during loading compared to the fibrin-rich clots. Results also reveal that the mechanical behaviour of clot analogues is significantly different in compression and tension. Finite element cohesive zone modelling of clot fracture experiments show that fibrin fibres become highly aligned in the direction perpendicular to crack propagation, providing a significant toughening mechanism. The results presented in this study provide the first characterization of the composition-dependent fracture behaviour of blood clots and are of key importance for development of next-generation thrombectomy devices and clinical strategies. STATEMENT OF SIGNIFICANCE: This study provides a characterisation of the composition-dependent fracture toughness of blood clots. This entails the development of novel experimental techniques for fabrication and testing of blood clot compact tension fracture specimens. The study also develops cohesive zone models of fracture initiation and propagation in blood clots. Results reveal that the fracture resistance of fibrin-rich clots is significantly higher than red blood cell rich clots. Simulations also reveal that stretching and realignment of the fibrin network should be included in blood clot material models in order to accurately replicate compression-tension asymmetry and fibrin enhanced fracture toughness. The results of this study have potentially important clinical implications in terms of clot fracture risk and secondary embolization during mechanical thrombectomy procedures.

A new compressible hyperelastic model for the multi-axial deformation of blood clot occlusions in vessels

Mechanical thrombectomy can be significantly affected by the mechanical properties of the occluding thrombus. In this study, we provide the first characterisation of the volumetric behaviour of blood clots. We propose a new hyperelastic model for the volumetric and isochoric deformation of clot. We demonstrate that the proposed model provides significant improvements over established models in terms of accurate prediction of nonlinear stress-strain and volumetric behaviours of clots with low and high red blood cell compositions. We perform a rigorous investigation of the factors that govern clot occlusion of a tapered vessel. The motivation for such an analysis is twofold: (i) the role of clot composition on the in vivo occlusion location is an open clinical question that has significant implications for thrombectomy procedures; (ii) in vitro measurement of occlusion location in an engineered tapered tube can be used as a quick and simple methodology to assess the mechanical properties/compositions of clots. Simulations demonstrate that both isochoric and volumetric behaviours of clots are key determinants of clot lodgement location, in addition to clot-vessel friction. The proposed formulation is shown to provide accurate predictions of in vitro measurement of clot occlusion location in a silicone tapered vessel, in addition to accurately predicting the deformed shape of the clot.