Development of a computational model for acute ischemic stroke recanalization through cyclic aspiration

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.

Systems-level computational modeling in ischemic stroke: from cells to patients

Ischemic stroke, a significant threat to human life and health, refers to a class of conditions where brain tissue damage is induced following decreased cerebral blood flow. The incidence of ischemic stroke has been steadily increasing globally, and its disease mechanisms are highly complex and involve a multitude of biological mechanisms at various scales from genes all the way to the human body system that can affect the stroke onset, progression, treatment, and prognosis. To complement conventional experimental research methods, computational systems biology modeling can integrate and describe the pathogenic mechanisms of ischemic stroke across multiple biological scales and help identify emergent modulatory principles that drive disease progression and recovery. In addition, by running virtual experiments and trials in computers, these models can efficiently predict and evaluate outcomes of different treatment methods and thereby assist clinical decision-making. In this review, we summarize the current research and application of systems-level computational modeling in the field of ischemic stroke from the multiscale mechanism-based, physics-based and omics-based perspectives and discuss how modeling-driven research frameworks can deliver insights for future stroke research and drug development.

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.

Current Approaches and Methods to Understand Acute Ischemic Stroke Treatment Using Aspiration Thrombectomy

Acute ischemic stroke occurs when a blood clot occludes a cerebral artery. Mechanical interventions, primarily stent retrievers and aspiration thrombectomy, are used currently for removing the occluding clot and restoring blood flow. Aspiration involves using a long catheter to traverse the cerebral vasculature to reach the blood clot, followed by application of suction through the catheter bore. Aspiration is also used in conjunction with other techniques such as stent retrievers and balloon guide catheters. Despite the wide use of aspiration, our physical understanding of the process and the causes of the failure of aspiration to retrieve cerebral clots in certain scenarios is not well understood. Experimental and computational studies can help develop the capability to provide deeper insights into the procedure and enable development of new devices and more effective treatment methods. We recapitulate the aspiration-based thrombectomy techniques in clinical practice and provide a perspective of existing engineering methods for aspiration. We articulate the current knowledge gap in the understanding of aspiration and highlight possible directions for future engineering studies to bridge this gap, help clinical translation of engineering studies, and develop new patient-specific stroke therapy.

Decoding thrombosis through code: a review of computational models

From the molecular level up to a blood vessel, thrombosis and hemostasis involves many interconnected biochemical and biophysical processes over a wide range of length and time scales. Computational modeling has gained eminence in offering insights into these processes beyond what can be obtained from in vitro or in vivo experiments, or clinical measurements. The multiscale and multiphysics nature of thrombosis has inspired a wide range of modeling approaches that aim to address how a thrombus forms and dismantles. Here, we review recent advances in computational modeling with a focus on platelet-based thrombosis. We attempt to summarize the diverse range of modeling efforts straddling the wide-spectrum of physical phenomena, length scales, and time scales; highlighting key advancements and insights from existing studies. Potential information gleaned from models is discussed, ranging from identification of thrombus-prone regions in patient-specific vasculature to modeling thrombus deformation and embolization in response to fluid forces. Furthermore, we highlight several limitations of current models, future directions in the field, and opportunities for clinical translation, to illustrate the state-of-the-art. There are a plethora of opportunity areas for which models can be expanded, ranging from topics of thromboinflammation to platelet production and clearance. Through successes demonstrated in existing studies described here, as well as continued advancements in computational methodologies and computer processing speeds and memory, in silico investigations in thrombosis are poised to bring about significant knowledge growth in the years to come.

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.

Computational analysis of effects of clot length on Acute ischemic stroke recanalization under different cyclic aspiration loading conditions

Acute ischemic stroke, the second leading cause of death worldwide, results from occlusion of a cerebral artery by a blood clot. Application of cyclic aspiration using an aspiration catheter is a current therapy for the removal of lodged clots. In this study, we perform finite element simulations to analyze deformation of long clots, having length to radius ratio of 2-10, which corresponds to clot-length of 2.85-14.25 mm, under peak-to-peak cyclic aspiration pressures of 10-50 mmHg, and frequencies of 0.5, 1, and 2 Hz. Our computational system comprises of a nonlinear viscoelastic solid clot, a hyperelastic artery, and a nonlinear viscoelastic cohesive zone, the latter modeling the clot-artery interface. We observe that clots having length-to-radius ratio approximately greater than two separate from the inner arterial surface somewhere between the axial and distal ends, irrespective of the cyclic aspiration loading conditions. The stress distribution within the clot shows large tensile stresses in the clot interior, indicating the possibility of simultaneous fragmentation of the clot. Thus, this study shows us the various failure mechanisms simultaneously present in the clot during cyclic aspiration. Similarly, the stress distribution within the artery implies a possibility of endothelial damage to the arterial wall near the end where the aspiration pressure is applied. This framework provides a foundation for further investigation to clot fracture and adhesion characterization.

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.

Preclinical modeling of mechanical thrombectomy

Mechanical thrombectomy to treat large vessel occlusions (LVO) causing a stroke is one of the most effective treatments in medicine, with a number needed to treat to improve clinical outcomes as low as 2.6. As the name implies, it is a mechanical solution to a blocked artery and modeling these mechanics preclinically for device design, regulatory clearance and high-fidelity physician training made clinical applications possible. In vitro simulation of LVO is extensively used to characterize device performance in representative vascular anatomies with physiologically accurate hemodynamics. Embolus analogues, validated against clots extracted from patients, provide a realistic simulated use experience. In vitro experimentation produces quantitative results such as particle analysis of distal emboli generated during the procedure, as well as pressure and flow throughout the experiment. Animal modeling, used mostly for regulatory review, allows estimation of device safety. Other than one recent development, nearly all animal modeling does not incorporate the desired target organ, the brain, but rather is performed in the extracranial circulation. Computational modeling of the procedure remains at the earliest stages but represents an enormous opportunity to rapidly characterize and iterate new thrombectomy concepts as well as optimize procedure workflow. No preclinical model is a perfect surrogate; however, models available can answer important questions during device development and have to date been successful in delivering efficacious and safe devices producing excellent clinical outcomes. This review reflects on the developments of preclinical modeling of mechanical thrombectomy with particular focus on clinical translation, as well as articulate existing gaps requiring additional research.

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.

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.

Modeling acute ischemic stroke recanalization through cyclic aspiration

We model the deformation of a thromboembolus lodged in a cerebral artery under the application of aspiration pressure as it would be provided by an aspiration catheter during a mechanical thrombectomy procedure. The system considered consists of (i) a clot modeled as a viscoelastic solid; (ii) an artery modeled as a hyperelastic solid; and (iii) a viscoelastic cohesive interface between the clot and the artery. For the chosen system and geometry, we show that the application of aspiration pressure results in the impingement of the thrombus against the inner arterial wall near the aspiration location. Conditions leading to interfacial failure are nucleated at the distal end of the clot and, depending on the details of the loading conditions, propagate toward the proximal end. The results provide useful information in identifying the circumstances that play a decisive role for clot removal by aspiration alone.

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.

Hydrodynamics in Acute Ischemic Stroke Catheters Under Static and Cyclic Aspiration Conditions

PURPOSE

Previous studies have suggested improved recanalization efficiency by using cyclic aspiration as opposed to static aspiration for mechanical thrombectomy in the treatment of acute ischemic stroke (AIS). However, there has not been an investigation into which parameters optimize this theoretical approach. Therefore, this study was designed to investigate the application of static and cyclic pressure in AIS aspiration catheters.

METHODS

An experimental flow system was designed to apply and measure both static and cyclic aspiration pressures in four commercial aspiration catheters. A medical aspirator and solenoid-valve were used to generate a range of applied pump pressures and cyclic valve frequencies. The measured mean and oscillatory pressure amplitudes were analyzed and compared across catheters and operating conditions.

RESULTS

The greatest mean pressures and forces were achieved under contact-aspiration conditions using the largest diameter catheter, while the largest oscillating pressure (ΔP) and force (ΔF) values were generated in the least compliant catheter, determined to be proportional to both its inner radius and wall thickness. Furthermore, lower frequency cyclic aspiration was able to generate the largest values of ΔP and ΔF.

CONCLUSION

For future investigations of cyclic aspiration AIS therapy, the best options are to use large inner diameter and low compliance catheters operating at lower cyclic frequencies to generate the greatest mean and oscillating forces on lodged thromboemboli.