Mechanisms by which thrombolytic therapy results in nonuniform lysis and residual thrombus after reperfusion

Computational modeling of blood clot lysis considering the effect of vessel wall and pulsatile blood flow

Stroke is one of the major causes of global death, which can occur due to blockage in a blood vessel by a clot. The immediate dissolving of the clot is essential to restore the blood flow and prevent tissue necrosis. Clot dissolution can be achieved via thrombolytic therapy using plasminogen activators. In this study, a clot dissolution model is developed for a three-dimensional patient-specific carotid artery that investigates the effect of different vessel wall models on clot dissolution. The lysis pattern of the clot and hemodynamics of blood flow are evaluated using three different models of the vessel wall, namely, rigid, linear elastic, and Mooney-Rivlin hyperelastic. The effect of flow condition is considered by solving the Navier-Stokes equations for the free flow domain and the Brinkman equation for the clot domain with the same pressure and velocity fields. This will result in continuous pressure and velocity over the interfaces of the free flow and clot domains. The blood inflow is assumed to be pulsatile. In addition, the species transport driven by diffusion and convection is considered to be different in the porous medium and plasma. The obtained results show that in all models, the starting time of clot volume decrease is almost the same and the clot starts dissolving from the inner curvature of the artery. However, in the hyperelastic model, dissolving the clot takes longer compared to the other two models. By monitoring the vessel wall deformation, the exact time of vessel recanalization is determined.

Blood clot contraction: Mechanisms, pathophysiology, and disease

A State of the Art lecture titled "Blood Clot Contraction: Mechanisms, Pathophysiology, and Disease" was presented at the International Society on Thrombosis and Haemostasis (ISTH) Congress in 2022. This was a systematic description of blood clot contraction or retraction, driven by activated platelets and causing compaction of the fibrin network along with compression of the embedded erythrocytes. The consequences of clot contraction include redistribution of the fibrin-platelet meshwork toward the periphery of the clot and condensation of erythrocytes in the core, followed by their deformation from the biconcave shape into polyhedral cells (polyhedrocytes). These structural signatures of contraction have been found in ex vivo thrombi derived from various locations, which indicated that clots undergo intravital contraction within the blood vessels. In hemostatic clots, tightly packed polyhedrocytes make a nearly impermeable seal that stems bleeding and is impaired in hemorrhagic disorders. In thrombosis, contraction facilitates the local blood flow by decreasing thrombus obstructiveness, reducing permeability, and changing susceptibility to fibrinolytic enzymes. However, in (pro)thrombotic conditions, continuous background platelet activation is followed by platelet exhaustion, refractoriness, and impaired intravital clot contraction, which is associated with weaker thrombi predisposed to embolization. Therefore, assays that detect imperfect in vitro clot contraction have potential diagnostic and prognostic values for imminent or ongoing thrombosis and thrombotic embolism. Collectively, the contraction of blood clots and thrombi is an underappreciated and understudied process that has a pathogenic and clinical significance in bleeding and thrombosis of various etiologies. Finally, we have summarized relevant new data on this topic presented during the 2022 ISTH Congress.

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.

Microstructure aware modeling of biochemical transport in arterial blood clots

Flow-mediated transport of biochemical species is central to thrombotic phenomena. Comprehensive three-dimensional modeling of flow-mediated transport around realistic macroscale thrombi poses challenges owing to their arbitrary heterogeneous microstructure. Here, we develop a microstructure aware model for species transport within and around a macroscale thrombus by devising a custom preconditioned fictitious domain formulation for thrombus-hemodynamics interactions, and coupling it with a fictitious domain advection-diffusion formulation for transport. Microstructural heterogeneities are accounted through a hybrid discrete particle-continuum approach for the thrombus interior. We present systematic numerical investigations on unsteady arterial flow within and around a three-dimensional macroscale thrombus; demonstrate the formation of coherent flow structures around the thrombus which organize advective transport; illustrate the role of the permeation processes at the thrombus boundary and subsequent intra-thrombus transport; and characterize species transport from bulk flow to the thrombus boundary and vice versa.

Computational Simulations of Thrombolytic Therapy in Acute Ischaemic Stroke

Ischaemic stroke can occur when an artery to the brain is blocked by a blood clot. The use of thrombolytic agents, such as tissue plasminogen activator (tPA), to dissolve the occluding clot is limited by the risk of intracerebral haemorrhage (ICH), a known side effect associated with tPA. We developed a computational thrombolysis model for a 3D patient-specific artery coupled with a compartmental model for temporal concentrations of tPA and lysis proteins during intravenous infusion of tPA, in order to evaluate the effects of tPA dose on the efficacy of thrombolytic therapy and the risk of ICH. The model was applied to a 3-mm-long fibrin clot with two different fibrin fibre radii in the middle cerebral artery (MCA) - a setting relevant to ischaemic stroke, and results for different tPA dose levels and fibrin fibre radii were compared. Our simulation results showed that clot lysis was accelerated at higher tPA doses at the expense of a substantial increase in the risk of ICH. It was also found that a fine clot with a smaller fibre radius dissolved much slowly than a coarse clot due to a slower tPA penetration into the clots.

Systems Analysis of Thrombus Formation

The systems analysis of thrombosis seeks to quantitatively predict blood function in a given vascular wall and hemodynamic context. Relevant to both venous and arterial thrombosis, a Blood Systems Biology approach should provide metrics for rate and molecular mechanisms of clot growth, thrombotic risk, pharmacological response, and utility of new therapeutic targets. As a rapidly created multicellular aggregate with a polymerized fibrin matrix, blood clots result from hundreds of unique reactions within and around platelets propagating in space and time under hemodynamic conditions. Coronary artery thrombosis is dominated by atherosclerotic plaque rupture, complex pulsatile flows through stenotic regions producing high wall shear stresses, and plaque-derived tissue factor driving thrombin production. In contrast, venous thrombosis is dominated by stasis or depressed flows, endothelial inflammation, white blood cell-derived tissue factor, and ample red blood cell incorporation. By imaging vessels, patient-specific assessment using computational fluid dynamics provides an estimate of local hemodynamics and fractional flow reserve. High-dimensional ex vivo phenotyping of platelet and coagulation can now power multiscale computer simulations at the subcellular to cellular to whole vessel scale of heart attacks or strokes. In addition, an integrated systems biology approach can rank safety and efficacy metrics of various pharmacological interventions or clinical trial designs.

Effect of encapsulation on plasminogen activator delivery to the microcirculation and its implications for bleeding

BACKGROUND AND PURPOSE

It is known that encapsulation can alter the delivery of plasminogen activators by flow to accelerate fibrinolysis while other experimental studies suggest encapsulation may reduce the risk of hemorrhage with administration of the agent. The aim of this research is to resolve the effect of encapsulation on fibrinolysis and bleeding in the microcirculation.

METHODS

An established rabbit model of fibrinolytic hemorrhage was utilized to explore the potential of encapsulation to limit bleeding. Equal dosages of free or microencapsulated streptokinase (MESK) were infused to initiate thrombolysis of small vessel clots while tracking blood loss.

RESULTS

Compared to free streptokinase, significant improvements in bleeding were observed with MESK as demonstrated by (1) delayed onset of bleeding, (2) shortened duration, and (3) reduction in the volume of lost blood, consistent with less systemic fibrinogen degradation.

CONCLUSIONS

Findings demonstrate that encapsulation of streptokinase can inhibit clot lysis in small vessels. Combined with prior work on accelerated thrombolysis, results suggest a time-based regimen for avoiding bleeding complications during thrombolytic therapy with encapsulated agent.

Towards a multi-physics modelling framework for thrombolysis under the influence of blood flow

Thrombolytic therapy is an effective means of treating thromboembolic diseases but can also give rise to life-threatening side effects. The infusion of a high drug concentration can provoke internal bleeding while an insufficient dose can lead to artery reocclusion. It is hoped that mathematical modelling of the process of clot lysis can lead to a better understanding and improvement of thrombolytic therapy. To this end, a multi-physics continuum model has been developed to simulate the dissolution of clot over time upon the addition of tissue plasminogen activator (tPA). The transport of tPA and other lytic proteins is modelled by a set of reaction–diffusion–convection equations, while blood flow is described by volume-averaged continuity and momentum equations. The clot is modelled as a fibrous porous medium with its properties being determined as a function of the fibrin fibre radius and voidage of the clot. A unique feature of the model is that it is capable of simulating the entire lytic process from the initial phase of lysis of an occlusive thrombus (diffusion-limited transport), the process of recanalization, to post-canalization thrombolysis under the influence of convective blood flow. The model has been used to examine the dissolution of a fully occluding clot in a simplified artery at different pressure drops. Our predicted lytic front velocities during the initial stage of lysis agree well with experimental and computational results reported by others. Following canalization, clot lysis patterns are strongly influenced by local flow patterns, which are symmetric at low pressure drops, but asymmetric at higher pressure drops, which give rise to larger recirculation regions and extended areas of intense drug accumulation.

An In Vitro Thrombolysis Study Using a Mixture of Fast-Acting and Slower Release Microspheres

PURPOSE

To test the hypothesis that a mixture combining fast and slower release rate microspheres can restore blood flow rapidly and prevent formation of another blockage in thrombolysis.

METHODS

We used polyethylene glycol (PEG) microspheres which provide the release of the encapsulated streptokinase (SK) on the scale of minutes, and Eudragit FS30D (Eud), a polymethacrylate polymer, for development of delayed release microspheres which were desirable to prevent a putative second thrombus. Eud microspheres were coated with chitosan (CS) to further extend half-life. Experiments included the development, characterization of Eud/SK and CS-Eud/SK microspheres, and in vitro thrombolytic studies of the mixtures of PEG/SK and Eud /SK microspheres and of PEG/SK and CS-Eud/SK microspheres.

RESULTS

CS-Eud/SK microspheres have slightly lower encapsulation efficiency, reduced activity of SK, and a much slower release of SK when compared with microspheres of Eud/SK microspheres. Counter-intuitively, slower release leads to faster thrombolysis after reocclusion as a result of greater retention of agent and the mechanism of distributed intraclot thrombolysis.

CONCLUSIONS

A mixture of PEG/SK and CS-Eud/SK microspheres could break up the blood clot rapidly while providing clot-lytic efficacy in prevention of a second blockage up to 4 h.

Evaluation of thrombolysis by using ultrasonic imaging: an in vitro study

The hematocrit of a thrombus is a key factor associated with the susceptibility to thrombolysis. Ultrasonic imaging is currently the first-line screening tool for thrombus examinations. Different hematocrits result in different acoustical structures of thrombi, which alter the behavior of ultrasonic backscattering. This study explored the relationships among thrombolytic efficiencies, hematocrits, and ultrasonic parameters (the echo intensity and backscattered statistics). Porcine thrombi with different hematocrits, ranging from 0% to 50%, were induced in vitro. An ultrasonic scanner was used to scan thrombi and acquire raw image data for B-mode (echo intensity measurements) and Nakagami imaging (backscattered statistics analysis). Experiments on thrombolysis were performed using urokinase to explore the effect of the hematocrit on thrombolytic efficiency. Results showed that the weight loss ratio of thrombi exponentially decreased as the hematocrit increased from 0% to 50%. Compared with the echo intensity obtained from the conventional B-scan, the Nakagami parameter predicts the weight loss ratio, increasing from 0.6 to 1.2 as the weight loss ratio decreased from 0.67 to 0.26. The current findings suggest that using Nakagami imaging characterizing thrombi provides information of backscattered statistics, which may be associated with the thrombolytic efficiency.

Computational analysis of blood clot dissolution using a vibrating catheter tip

We developed a novel concept of endovascular thrombolysis that employs a vibrating electroactive polymer actuator. In order to predict the efficacy of thrombolysis using the developed vibrating actuator, enzyme (plasminogen activator) perfusion into a clot was analyzed by solving flow fields and species transport equations considering the fluid structure interaction. In vitro thrombolysis experiments were also performed. Computational results showed that plasminogen activator perfusion into a clot was enhanced by actuator vibration at frequencies of 1 and 5 Hz. Plasminogen activator perfusion was affected by the actuator oscillation frequencies and amplitudes that were determined by electromechanical characteristics of a polymer actuator. Computed plasminogen activator perfused volumes were compared with experimentally measured dissolved clot volumes. The computed plasminogen activator perfusion volumes with threshold concentrations of 16% of the initial plasminogen activator concentration agreed well with the in vitro experimental data. This study showed the effectiveness of actuator oscillation on thrombolysis and the validity of the computational plasminogen activator perfusion model for predicting thrombolysis in complex flow fields induced by an oscillating actuator.

Accelerating thrombolysis with chitosan-coated plasminogen activators encapsulated in poly-(lactide-co-glycolide) (PLGA) nanoparticles

Accelerating thrombolysis using plasminogen activators (PAs) encapsulated liposomes or PEG microparticles by pressure-driven permeation have been demonstrated in vitro and in vivo in animal models. However, designing and delivering PA-encapsulated nanoparticles (NPs) to enhance thrombolysis by applying electrostatic forces or ligand-receptor interactions between the NPs and blood clots has not been proposed. Therefore, without a pressure-driving factor, tissue-plasminogen activator (t-PA) encapsulated in PLGA NPs with chitosan (CS) and CS-GRGD coating and their thrombolysis capabilities in a blood clot-occluded tube model were evaluated by determining clot lysis times and the masses of the digested clots. The characteristics and release profiles of t-PA-encapsulated PLGA, PLGA/CS and PLGA/CS-GRGD NPs are determined by FT-IR, a laser particle/zeta potential analyzer and HPLC. Additionally, the permeation capacities of the NPs into flat blood clots were examined. For example, the mean particle sizes and encapsulation efficacies of t-PA for the NPs are in the ranges 260-320 nm and 65.5-70.5%, respectively. The results reveal that the NPs for the shortest clot lysis time and the highest weight percentages of digested clot are PLGA/CS (20.7 +/- 0.7 min) and PLGA/CS-GRGD (25.7 +/- 1.3 wt%), respectively. Compared with t-PA solution, the NPs can significantly shorten clot lysis times in the following order: PLGA/CS NPs (38.8 +/- 1.5%) > PLGA/CS-GRGD NPs (16.3 +/- 1.0%) > PLGA NPs (7.7 +/- 1.2%). Compared with t-PA solution, the NPs significantly increase the weight of digested clots in the order, PLGA/CS-GRGD (40.9 +/- 1.5%) > PLGA/CS (27.8 +/- 1.2%) > PLGA (8.6 +/- 0.6%). The highest release rate of t-PA in the fast release phase and the highest permeability into intra-clots of PLGA/CS and PLGA/CS-GRGD NPs, respectively, correspond to the shortest clot lysis time and the largest increase in weight of the digested clots among the NP system. In conclusion, the NPs designed based on new concepts significantly accelerate thrombolysis in vitro in this model, and may be useful in clinical study.

Modelling the effect of laminar axially directed blood flow on the dissolution of non-occlusive blood clots

Axially directed blood plasma flow can significantly accelerate thrombolysis of non-occlusive blood clots. Viscous forces caused by shearing of blood play an essential role in this process, in addition to biochemical fibrinolytic reactions. An analytical mathematical model based on the hypothesis that clot dissolution dynamics is proportional to the power of the flowing blood plasma dissipated along the clot is presented. The model assumes cylindrical non-occlusive blood clots with the flow channel in the centre, in which the flow is assumed to be laminar and flow rate constant at all times during dissolution. Effects of sudden constriction on the flow and its impact on the dissolution rate are also considered. The model was verified experimentally by dynamic magnetic resonance (MR) microscopy of artificial blood clots dissolving in an in vitro circulation system, containing plasma with a magnetic resonance imaging contrast agent and recombinant tissue-type plasminogen activator (rt-PA). Sequences of dynamically acquired 3D low resolution MR images of entire clots and 2D high resolution MR images of clots in the axial cross-section were used to evaluate the dissolution model by fitting it to the experimental data. The experimental data fitted well to the model and confirmed our hypothesis.

Numerical analysis of forced injection of enzyme during thrombolysis

Numerical analysis was performed on the enzyme transport and the flow fields in order to predict the effectiveness of forced injection in thrombolytic therapy. The species and momentum transport equations were numerically solved for the case of uniform perfusion of enzyme into the fibrin clot, and the validity of our methods were verified. In order to predict the lysis efficiency of continuous and forced intermittent injections, enzyme perfusion and clot lysis were simulated for the different injection velocities and frequencies. Intermittent injection showed faster clot lysis compared to continuous perfusion, and lysis efficiency was increased as the injection velocity and period increased.

Blood Clot Dissolution Dynamics Simulation during Thrombolytic Therapy

Nonocclusive blood clots only partially fill blood vessels and together with the adjacent vessel wall form a channel through which blood flows at usually much higher velocities than in normal vessels. Our aim was to find a theoretical explanation for the experimentally observed fact that fast flowing blood through the channel has a large effect on the increase of the clot dissolution rate compared to the dissolution rate in the absence of flow. Blood flow through the channel increases transport of dissolution agents to the clot and also exerts large forces to the surface of the clot along the channel. Proposed is a model for clot dissolution which assumes that the clot dissolution rate is proportional to the forces of flowing blood to the surface of the clot multiplied by the average blood velocity. The model has been verified by fitting to experimental magnetic resonance imaging data obtained by dynamical magnetic resonance microscopy of clots dissolved by recombinant tissue plasminogen activator in an artificial blood flow system.

Distributed intraclot thrombolysis: mechanism of accelerated thrombolysis with encapsulated plasminogen activators

BACKGROUND

The delivery of encapsulated plasminogen activators has demonstrated enhanced thrombolysis in vivo in several models. The mechanism of such improvement has not previously been established.

OBJECTIVES

We explored in vitro the mechanism by which microencapsulation of streptokinase in polymeric microparticles accelerates clot digestion and reduces reperfusion times by as much as an order of magnitude in vivo.

METHODS

The efficacy of microencapsulated streptokinase (MESK) was directly compared with identical dosages of unencapsulated streptokinase (FREE SK) at three initial pressure drops using clots formed of plasma or whole blood in 0.2-cm inner diameter glass capillary tubes.

RESULTS

MESK demonstrated accelerated flow restoration compared with FREE SK for each condition in plasma (23.8 +/- 4.5% faster) and whole blood clots (17.2 +/- 9.2% faster). Images collected by light microscopy show sites of thrombolysis internal to the clot only with MESK while the spatial distribution of fluorescently labeled streptokinase by confocal microscopy confirms greater penetration of the encapsulated agent compared with unencapsulated streptokinase. Digestion thus proceeds in three dimensions rather than restricted to a two-dimensional lysis front.

CONCLUSIONS

The improved clot penetration with MESK establishes enhanced transport with encapsulation and the concept of distributed intraclot thrombolysis as a basis for the accelerated dissolution observed with encapsulated plasminogen activators in vivo.

A mathematical model of post-canalization thrombolysis

During the initial phase of lysis of an occlusive thrombus using lytic agents such as tissue plasminogen activator, blood flow through the centre of the clot is established (the process of recanalization). Following canalization, the clot remains on the vessel wall and further lysis is required. This paper develops a multi-species mathematical model to describe the bulk chemical reactions in the bloodstream and the convective and diffusive transport of chemical species to and from the clot surface in conditions following canalization. For the steady state case, the model indicates that the process of clot lysis following initial recanalization is dominated by surface chemical reactions and the bulk reactions play little role in the lytic process. Lytic rate is dependent on the clot geometry and flow conditions. The rate of clot dissolution is greatest at the upstream end of the clot and decreases steadily downstream due to lytic agent being removed from the flowing blood as it binds to the clot surface. This model may be further developed and used to simulate and compare different lytic regimes.

Engineering design of optimal strategies for blood clot dissolution

Blood clots form under hemodynamic conditions and can obstruct flow during angina, acute myocardial infarction, stroke, deep vein thrombosis, pulmonary embolism, peripheral thrombosis, or dialysis access graft thrombosis. Therapies to remove these clots through enzymatic and/or mechanical approaches require consideration of the biochemistry and structure of blood clots in conjunction with local transport phenomena. Because blood clots are porous objects exposed to local hemodynamic forces, pressure-driven interstitial permeation often controls drug penetration and the overall lysis rate of an occlusive thrombus. Reaction engineering and transport phenomena provide a framework to relate dosage of a given agent to potential outcomes. The design and testing of thrombolytic agents and the design of therapies must account for (a) the binding, catalytic, and systemic clearance properties of the therapeutic enzyme; (b) the dose and delivery regimen; (c) the biochemical and structural aspects of the thrombotic occlusion; (d) the prevailing hemodynamics and anatomical location of the thrombus; and (e) therapeutic constraints and risks of side effects. These principles also impact the design and analysis of local delivery devices.

Differential effects of c7E3 Fab on thrombus formation and rt-PA-Mediated thrombolysis under flow conditions

Although the Fab fragment of the mouse-human chimeric anti-alphaIIbbeta3 (GP IIb/IIIa) monoclonal antibody (MoAb) c7E3 facilitates recombinant tissue-type plasminogen activator (rt-PA)-mediated thrombolysis, it is not clear whether this is due to inhibition of new clot formation and/or a direct effect on the lysis rate. We employed an in vitro flow (re)circulation model to investigate how c7E3 Fab affected (a) platelet adhesion to clotted fibrin substrates under laminar flow at wall shear rates of 100 or 500 s(-1) and (b) rt-PA-induced lysis of preformed mural platelet-fibrin substrates at 500 s(-1). c7E3 Fab dose-dependently (0.5-5 microg/ml) inhibited platelet adhesion from flowing whole blood onto fibrin substrates ( approximately 14 microm thick) at each wall shear rate. When at 5 min after the onset of flow, c7E3 Fab (0.1-10 microg/ml) and rt-PA (1 microg/ml) were coinjected in flowing blood, it was found that modest fibrinolysis caused major platelet release from fibrin substrates and there was no difference in the lysis rate in the presence of rt-PA + c7E3 Fab compared to rt-PA alone. Platelet pretreatment with c7E3 Fab (10 microg/ml) had no effect on the lysis rate of thin ( approximately 40 microm), and slightly delayed the lysis rate of thick (< 250 microm), platelet-fibrin substrates containing evenly dispersed platelets (10(9)/ml). When the platelets within thick platelet-fibrin substrates were organized in platelet-rich regions ("residual thrombi"), these substrates followed a nonuniform lysis pattern, where fibrin between the thrombi lysed first and the residual thrombi lysed at a slower rate. Platelet pretreatment with c7E3 Fab (10 microg/ml) abolished the formation of the lytic-resistant residual thrombi and the associated platelet-protected fibrin zones. Hence, treatment with c7E3 Fab has no direct effect on the rate of rt-PA-mediated lysis, but is expected to block platelet-fibrin interactions that lead to clot retraction, thus maintaining a fibrin architecture that is more susceptible to lysis.