Development of a flow-through system to create occluding thrombus

Occlusive thrombosis accounts for many heart attacks and strokes. These acute events are difficult to catch in patients and animal test methods may be misleading because anti-thrombotic therapeutics often do not cross-react with different species. This paper presents a new flow-through system that leads to rapid occlusive thrombosis in arterial flow conditions. Whole porcine blood is perfused through a tubular test section. The growing thrombus is visualized in real time from early platelet attachment, through accumulation, to occlusion. The progression of flow rate reduction provides a clear distinguishing parameter between thrombus formation and embolization. Thrombus growth rate is a linear function of very high shear rate beyond 40,000 s(-1). The histology of the thrombus reveals predominantly platelet accumulation and growth as a rough surface with tendrils. This flow-through system may be useful for the economic testing of new anti-thrombosis therapies.     

A computational model for false lumen thrombosis in type B aortic dissection following thoracic endovascular repair

Thoracic endovascular repair (TEVAR) has recently been established as the preferred treatment option for complicated type B dissection. This procedure involves covering the primary entry tear to stimulate aortic remodelling and promote false lumen thrombosis thereby restoring true lumen flow. However, complications associated with incomplete false lumen thrombosis, such as aortic dilatation and stent graft induced new entry tears, can arise after TEVAR. This study presents the application and validation of a recently developed mathematical model for patient-specific prediction of thrombus formation and growth under physiologically realistic flow conditions. The model predicts thrombosis through the evaluation of shear rates, fluid residence time and platelet distribution, based on convection-diffusion-reaction transport equations. The model was applied to 3 type B aortic dissection patients: two TEVAR cases showing complete and incomplete false lumen thrombosis respectively, and one medically treated dissection with no signs of thrombosis. Predicted thrombus growth over time was validated against follow-up CT scans, showing good agreement with in vivo data in all cases with a maximum difference between predicted and measured false lumen reduction below 8%. Our results demonstrate that TEVAR-induced thrombus formation in type B aortic dissection can be predicted based on patient-specific anatomy and physiologically realistic boundary conditions. Our model can be used to identify anatomical or stent graft related factors that are associated with incomplete false lumen thrombosis following TEVAR, which may help clinicians develop personalised treatment plans for dissection patients in the future.     

von Willebrand factor unfolding mediates platelet deposition in a model of high-shear thrombosis

Thrombosis under high-shear conditions is mediated by the mechanosensitive blood glycoprotein von Willebrand factor (vWF). vWF unfolds in response to strong flow gradients and facilitates rapid recruitment of platelets in flowing blood. While the thrombogenic effect of vWF is well recognized, its conformational response in complex flows has largely been omitted from numerical models of thrombosis. We recently presented a continuum model for the unfolding of vWF, where we represented vWF transport and its flow-induced conformational change using convection-diffusion-reaction equations. Here, we incorporate the vWF component into our multi-constituent model of thrombosis, where the local concentration of stretched vWF amplifies the deposition rate of free-flowing platelets and reduces the shear cleaning of deposited platelets. We validate the model using three benchmarks: in vitro model of atherothrombosis, a stagnation point flow, and the PFA-100, a clinical blood test commonly used for screening for von Willebrand disease (vWD). The simulations reproduced the key aspects of vWF-mediated thrombosis observed in these experiments, such as the thrombus location, thrombus growth dynamics, and the effect of blocking platelet-vWF interactions. The PFA-100 simulations closely matched the reported occlusion times for normal blood and several hemostatic deficiencies, namely, thrombocytopenia, vWD type 1, and vWD type 3. Overall, this multi-constituent model of thrombosis enables macro-scale 3D simulations of thrombus formation in complex geometries over a wide range of shear rates and accounts for qualitative and quantitative hemostatic deficiencies in patient blood.     

Determining possible thrombus sites in an extracorporeal device,using computational fluid dynamics-derived relative residence time

The prediction of conditions that may result in thrombus formation is a useful application of computational fluid dynamics. A number of techniques exist, based on the consideration of wall shear stress and regions of low blood flow; however, no clear guideline exists for the best practice of their use. In this paper, the sensitivity of each parameter and the specific mechanical forces are explained, before the optimal indicator of thrombosis risk is outlined. An extracorporeal access device cavity provides a suitable geometry to test the methodology. The recommended method for thrombus prediction considers areas with a calculated residence time (RT) and shear strain rate (SSR) thresholds, here set to RT>1 and SSR < 10 s^? 1. Evidence of thrombosis was found for physiological waveforms with an absence of reverse flow, which is expected to ‘wash out’ the cavity. The predicted thrombosis sites compare well with evidence collected from explanted devices.     

Platelet Transport Rates and Binding Kinetics at High Shear over a?Thrombus

Thrombus formation over a ruptured atherosclerotic plaque cap can occlude an artery with fatal consequences. We describe a computational model of platelet transport and binding to interpret rate-limiting steps seen in experimental thrombus formation over a collagen-coated stenosis. The model is used to compute shear rates in stenoses with growing boundaries. In the model, moving erythrocytes influence platelet transport based on shear-dependent enhanced diffusivity and a nonuniform platelet distribution. Adhesion is modeled as platelet-platelet binding kinetics. The results indicate that observed thrombus growth rates are limited by platelet transport to the wall for shear rates up to 6000 s⁻¹. Above 7000 s⁻¹, the thrombus growth rate is likely limited by binding kinetics (10⁻⁴ m/s). Thrombus growth computed from these rate-limiting steps match the thrombus location and occlusion times for experimental conditions if a lag time for platelet activation is included. Using fitted parameters, the model is then used to predict thrombus size and shape at a higher Reynolds number flow consistent with coronary artery disease.     

Microfluidic Thrombosis under Multiple Shear Rates and Antiplatelet Therapy Doses

The mainstay of treatment for thrombosis, the formation of occlusive platelet aggregates that often lead to heart attack and stroke, is antiplatelet therapy. Antiplatelet therapy dosing and resistance are poorly understood, leading to potential incorrect and ineffective dosing. Shear rate is also suspected to play a major role in thrombosis, but instrumentation to measure its influence has been limited by flow conditions, agonist use, and non-systematic and/or non-quantitative studies.In this work we measured occlusion times and thrombus detachment for a range of initial shear rates (500, 1500, 4000, and 10000 s⁻¹) and therapy concentrations (0–2.4 µM for eptifibatide, 0–2 mM for acetyl-salicylic acid (ASA), 3.5–40 Units/L for heparin) using a microfluidic device. We also measured complete blood counts (CBC) and platelet activity using whole blood impedance aggregometry. Effects of shear rate and dose were analyzed using general linear models, logistic regressions, and Cox proportional hazards models.Shear rates have significant effects on thrombosis/dose-response curves for all tested therapies. ASA has little effect on high shear occlusion times, even at very high doses (up to 20 times the recommended dose). Under ASA therapy, thrombi formed at high shear rates were 4 times more prone to detachment compared to those formed under control conditions. Eptifibatide reduced occlusion when controlling for shear rate and its efficacy increased with dose concentration. In contrast, the hazard of occlusion from ASA was several orders of magnitude higher than that of eptifibatide. Our results show similar dose efficacy to our low shear measurements using whole blood aggregometry. This quantitative and statistically validated study of the effects of a wide range of shear rate and antiplatelet therapy doses on occlusive thrombosis contributes to more accurate understanding of thrombosis and to models for optimizing patient treatment.     

Wall shear over high degree stenoses pertinent to atherothrombosis

Atherothrombosis can induce acute myocardial infarction and stroke by progressive stenosis of a blood vessel lumen to full occlusion. Since thrombus formation and embolization may be shear-dependent, we quantify the magnitude of shear rates in idealized severely stenotic coronary arteries (≥75% by diameter) using computational fluid dynamics to characterize the shear environment that may exist during atherothrombosis. Maximum shear rates in severe short stenoses were found to exceed 250,000 s^?1 (9500 dynes/cm²) and can reach a peak value of 425,000 s^?1 for a 98% stenosis. These high shear rates exceed typical shear used for in vitro blood flow experiments by an order of magnitude, indicating the need to examine thrombosis at very high shear rates. Pulsatility and stenosis eccentricity were found to have minor effects on the maximum wall shear rates in severe stenoses. In contrast, increases in the stenosis length reduced the maximum shear to 107,000 s^?1 (98% stenosis), while surface roughness could increase focal wall shear rates to a value reaching 610,000 s^?1 (90% stenosis). The “shear histories” of circulating platelets in these stenoses are far below reported activation thresholds. Platelets may be required to form bonds in 5 μs and resist shear forces reaching 8000 pN per platelet. Arterial thrombosis occurs in the face of pathological high shear stress, creating rapid and strong bonds without prior activation of circulating platelets.     

A three-dimensional multiscale model for the prediction of thrombus growth under flow with single-platelet resolution

Modeling thrombus growth in pathological flows allows evaluation of risk under patient-specific pharmacological, hematological, and hemodynamical conditions. We have developed a 3D multiscale framework for the prediction of thrombus growth under flow on a spatially resolved surface presenting collagen and tissue factor (TF). The multiscale framework is composed of four coupled modules: a Neural Network (NN) that accounts for platelet signaling, a Lattice Kinetic Monte Carlo (LKMC) simulation for tracking platelet positions, a Finite Volume Method (FVM) simulator for solving convection-diffusion-reaction equations describing agonist release and transport, and a Lattice Boltzmann (LB) flow solver for computing the blood flow field over the growing thrombus. A reduced model of the coagulation cascade was embedded into the framework to account for TF-driven thrombin production. The 3D model was first tested against in vitro microfluidics experiments of whole blood perfusion with various antiplatelet agents targeting COX-1, P₂Y₁, or the IP receptor. The model was able to accurately capture the evolution and morphology of the growing thrombus. Certain problems of 2D models for thrombus growth (artifactual dendritic growth) were naturally avoided with realistic trajectories of platelets in 3D flow. The generalizability of the 3D multiscale solver enabled simulations of important clinical situations, such as cylindrical blood vessels and acute flow narrowing (stenosis). Enhanced platelet-platelet bonding at pathologically high shear rates (e.g., von Willebrand factor unfolding) was required for accurately describing thrombus growth in stenotic flows. Overall, the approach allows consideration of patient-specific platelet signaling and vascular geometry for the prediction of thrombotic episodes.     

Structural and hydrodynamic simulation of an acute stenosis-dependent thrombosis model in mice

Platelet activation under blood flow is thought to be critically dependent on the autologous secretion of soluble platelet agonists (chemical activators) such as ADP and thromboxane. However, recent evidence challenging this model suggests that platelet activation can occur independent of soluble agonist signalling, in response to the mechanical effects of micro-scale shear gradients. A key experimental tool utilized to define the effect of shear gradients on platelet aggregation is the murine intravital microscopy model of platelet thrombosis under conditions of acute controlled arteriolar stenosis. This paper presents a computational structural and hydrodynamic simulation of acute stenotic blood flow in the small bowel mesenteric vessels of mice. Using a homogeneous fluid at low Reynolds number (0.45) we investigated the relationship between the local hydrodynamic strain-rates and the severity of arteriolar stensosis. We conclude that the critical rates of blood flow acceleration and deceleration at sites of artificially induced stenosis (vessel side-wall compression or ligation) are a function of tissue elasticity. By implementing a structural simulation of arteriolar side wall compression, we present a mechanistic model that provides accurate simulations of stenosis in vivo and allows for predictions of the effects on local haemodynamics in the murine small bowel mesenteric thrombosis model.     

Modeling the effect of blood vessel bifurcation ratio on occlusive thrombus formation

Vascular geometry is a major determinant of the hemodynamics that promote or prevent unnecessary vessel occlusion from thrombus formation. Bifurcations in the vascular geometry are repeating structures that introduce flow separation between parent and daughter vessels. We modelled the blood flow and shear rate in a bifurcation during thrombus formation and show that blood vessel bifurcation ratios determine the maximum shear rate on the surface of a growing thrombus. We built an analytical model that may aid in predicting microvascular bifurcation ratios that are prone to occlusive thrombus formation. We also observed that bifurcation ratios that adhere to Murray’s law of bifurcations may be protected from occlusive thrombus formation. These results may be useful in the rational design of diagnostic microfluidic devices and microfluidic blood oxygenators.     

Impact of Tissue Factor Localization on Blood Clot Structure and Resistance under Venous Shear

The structure and growth of a blood clot depend on the localization of tissue factor (TF), which can trigger clotting during the hemostatic process or promote thrombosis when exposed to blood under pathological conditions. We sought to understand how the growth, structure, and mechanical properties of clots under flow are shaped by the simultaneously varying TF surface density and its exposure area. We used an eight-channel microfluidic device equipped with a 20- or 100-μm-long collagen surface patterned with lipidated TF of surface densities ～0.1 and ～2 molecules/μm². Human whole blood was perfused at venous shear, and clot growth was continually measured. Using our recently developed computational model of clot formation, we performed simulations to gain insights into the clot’s structure and its resistance to blood flow. An increase in TF exposure area resulted not only in accelerated bulk platelet, thrombin, and fibrin accumulation, but also in increased height of the platelet mass and increased clot resistance to flow. Moreover, increasing the TF surface density or exposure area enhanced platelet deposition by approximately twofold, and thrombin and fibrin generation by greater than threefold, thereby increasing both clot size and its viscous resistance. Finally, TF effects on blood flow occlusion were more pronounced for the longer thrombogenic surface than for the shorter one. Our results suggest that TF surface density and its exposure area can independently enhance both the clot’s occlusivity and its resistance to blood flow. These findings provide, to our knowledge, new insights into how TF affects thrombus growth in time and space under flow.     

Computational Fluid Dynamics Simulations at Micro-Scale Stenosis for Microfluidic Thrombosis Model Characterization

Platelet aggregation plays a central role in pathological thrombosis, preventing healthy physiological blood flow within the circulatory system. For decades, it was believed that platelet aggregation was primarily driven by soluble agonists such as thrombin, adenosine diphosphate and thromboxane A2. However, recent experimental findings have unveiled an intriguing but complementary biomechanical mechanism—the shear rate gradients generated from flow disturbance occurring at sites of blood vessel narrowing, otherwise known as stenosis, may rapidly trigger platelet recruitment and subsequent aggregation. In our Nature Materials 2019 paper [1], we employed microfluidic devices which incorporated micro-scale stenoses to elucidate the molecular insights underlying the prothrombotic effect of blood flow disturbance. Nevertheless, the rheological mechanisms associated with this stenotic microfluidic device are poorly characterized. To this end, we developed a computational fluid dynamics (CFD) simulation approach to systematically analyze the hemodynamic influence of bulk flow mechanics and flow medium. Grid sensitivity studies were performed to ensure accurate and reliable results. Interestingly, the peak shear rate was significantly reduced with the device thickness, suggesting that fabrication of microfluidic devices should retain thicknesses greater than 50 µm to avoid unexpected hemodynamic aberration, despite thicker devices raising the cost of materials and processing time of photolithography. Overall, as many groups in the field have designed microfluidic devices to recapitulate the effect of shear rate gradients and investigate platelet aggregation, our numerical simulation study serves as a guideline for rigorous design and fabrication of microfluidic thrombosis models.     

A volumetric model for growth of arterial walls with arbitrary geometry and loads

Stress and deformation in arterial wall tissue are factors which may influence significantly its response and evolution. In this work we develop models based on nonlinear elasticity and finite element numerical solutions for the mechanical behaviour and the remodelling of the soft tissue of arteries, including anisotropy induced by the presence of collagen fibres. Remodelling and growth in particular constitute important features in order to interpret stenosis and atherosclerosis. The main object of this work is to model accurately volumetric growth, induced by fluid shear stress in the intima and local wall stress in arteries with patient-specific geometry and loads. The model is implemented in a nonlinear finite element setting which may be applied to realistic 3D geometries obtained from in vivo measurements. The capabilities of this method are demonstrated in several examples. Firstly a stenotic process on an idealised geometry induced by a non-uniform shear stress distribution is considered. Following the growth of a right coronary artery from an in vivo reconstructed geometry is presented. Finally, experimental measurements for growth under hypertension for rat carotid arteries are modelled.     

An integrated fluid–structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid–structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection–diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young’s modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

     

Plant food delphinidin-3-glucoside significantly inhibits platelet activation and thrombosis: novel protective roles against cardiovascular diseases

Delphinidin-3-glucoside (Dp-3-g) is one of the predominant bioactive compounds of anthocyanins in many plant foods. Although several anthocyanin compounds have been reported to be protective against cardiovascular diseases (CVDs), the direct effect of anthocyanins on platelets, the key players in atherothrombosis, has not been studied. The roles of Dp-3-g in platelet function are completely unknown. The present study investigated the effects of Dp-3-g on platelet activation and several thrombosis models in vitro and in vivo. We found that Dp-3-g significantly inhibited human and murine platelet aggregation in both platelet-rich plasma and purified platelets. It also markedly reduced thrombus growth in human and murine blood in perfusion chambers at both low and high shear rates. Using intravital microscopy, we observed that Dp-3-g decreased platelet deposition, destabilized thrombi, and prolonged the time required for vessel occlusion. Dp-3-g also significantly inhibited thrombus growth in a carotid artery thrombosis model. To elucidate the mechanisms, we examined platelet activation markers via flow cytometry and found that Dp-3-g significantly inhibited the expression of P-selectin, CD63, CD40L, which reflect platelet α- and δ-granule release, and cytosol protein secretion, respectively. We further demonstrated that Dp-3-g downregulated the expression of active integrin αIIbβ3 on platelets, and attenuated fibrinogen binding to platelets following agonist treatment, without interfering with the direct interaction between fibrinogen and integrin αIIbβ3. We found that Dp-3-g reduced phosphorylation of adenosine monophosphate-activated protein kinase, which may contribute to the observed inhibitory effects on platelet activation. Thus, Dp-3-g significantly inhibits platelet activation and attenuates thrombus growth at both arterial and venous shear stresses, which likely contributes to its protective roles against thrombosis and CVDs.     

Thrombotic risk stratification using computational modeling in patients with coronary artery aneurysms following Kawasaki disease

Kawasaki disease (KD) is the leading cause of acquired heart disease in children and can result in life-threatening coronary artery aneurysms in up to 25 % of patients. These aneurysms put patients at risk of thrombus formation, myocardial infarction, and sudden death. Clinicians must therefore decide which patients should be treated with anticoagulant medication, and/or surgical or percutaneous intervention. Current recommendations regarding initiation of anticoagulant therapy are based on anatomy alone with historical data suggesting that patients with aneurysms \(\ge \) 8 mm are at greatest risk of thrombosis. Given the multitude of variables that influence thrombus formation, we postulated that hemodynamic data derived from patient-specific simulations would more accurately predict risk of thrombosis than maximum diameter alone. Patient-specific blood flow simulations were performed on five KD patients with aneurysms and one KD patient with normal coronary arteries. Key hemodynamic and geometric parameters, including wall shear stress, particle residence time, and shape indices, were extracted from the models and simulations and compared with clinical outcomes. Preliminary fluid structure interaction simulations with radial expansion were performed, revealing modest differences in wall shear stress compared to the rigid wall case. Simulations provide compelling evidence that hemodynamic parameters may be a more accurate predictor of thrombotic risk than aneurysm diameter alone and motivate the need for follow-up studies with a larger cohort. These results suggest that a clinical index incorporating hemodynamic information be used in the future to select patients for anticoagulant therapy.     

A Novel Canine Model of Acute Vertebral Artery Occlusion

Background

The extended time window and theoretic reduction in hemorrhage make mechanical strategies an attractive approach for the treatment of patients with ischemic stroke. However, a limited availability of suitable animal models of cerebrovascular thrombosis has hampered the study of novel endovascular interventions. The aim of the present study was to develop a new technique for site-specific placement of a thrombus in a canine model that would allow for the evaluation of mechanical thrombectomy and clot retrieval methods and the visualization of thrombus dislocation or fragmentation during angiographic manipulation.

Methods

Angiography and embolization with a preformed thrombus were performed in 12 canines. Under fluoroscopic guidance, an embolism protection device (EPD) was anchored to the middle segment of the left vertebral artery (VA) via the left femoral arterial sheath. A preformed radiopaque clot was injected through the guide catheter into the left VA, via the contralateral femoral artery, proximal to the EPD. After 15 min of occlusion, the EPD was removed and persistent occlusion of the VA was documented angiographically.

Results

Angiography performed during the observation period confirmed the persistence of VA occlusion in each case, and displacement of the radiopaque clots did not occur during the 3-hour observation period. The technique allowed selective embolization of targeted vessels without thrombus fragmentation.

Conclusion

This study demonstrates, for the first time, a canine model of post-circulation embolism induced by autologous blood clot placement. This model can be rapidly formed and easily operated, and the site of thrombosis can be readily controlled.     

Developing pulsatile flow in a deployed coronary stent

A major consequence of stent implantation is restenosis that occurs due to neointimal formation. This patho-physiologic process of tissue growth may not be completely eliminated. Recent evidence suggests that there are several factors such as geometry and size of vessel, and stent design that alter hemodynamic parameters, including local wall shear stress distributions, all of which influence the restenosis process. The present three-dimensional analysis of developing pulsatile flow in a deployed coronary stent quantifies hemodynamic parameters and illustrates the changes in local wall shear stress distributions and their impact on restenosis. The present model evaluates the effect of entrance flow, where the stent is placed at the entrance region of a branched coronary artery. Stent geometry showed a complex three-dimensional variation of wall shear stress distributions within the stented region. Higher order of magnitude of wall shear stress of 530 dyn/cm2 is observed on the surface of cross-link intersections at the entrance of the stent. A low positive wall shear stress of 10 dyn/cm2 and a negative wall shear stress of -10 dyn/cm2 are seen at the immediate upstream and downstream regions of strut intersections, respectively. Modified oscillatory shear index is calculated which showed persistent recirculation at the downstream region of each strut intersection. The portions of the vessel where there is low and negative wall shear stress may represent locations of thrombus formation and platelet accumulation. The present results indicate that the immediate downstream regions of strut intersections are areas highly susceptible to restenosis, whereas a high shear stress at the strut intersection may cause platelet activation and free emboli formation.     

A New Near-wall Residence Time Model Applied to Three Arterio-venous Graft End-to-side Anastomoses

As a representative example of an end-to-side anastomosis, the distal portion of an arterio-venous hemodialysis graft with three graft-end designs and two outflow conditions has been considered. Arterio-venous graft failure for hemodialysis patients can be frequently related to stenotic developments, especially in the toe region. With critical blood particles such as LDL, platelets, and monocytes, being considered as a first-order approximation of; blood elements pathlines, near-wall residence times and wall stagnation points have been simulated to identify, in conjunction with the wall shear stress gradient, susceptible sites of stenotic developments. Assuming transient laminar three-dimensional non-Newtonian fluid flow, the transport equations were solved with validated finite-volume based CFD software and user-Fortran programs for the computation of fluid element pathlines and hemodynamic parameters. We conclude from this study that for the application of looped arterio-venous grafts, the graft-end geometry and the presence of retrograde outflow significantly influence the magnitude and location of the hemodynamic wall parameters considered. The new concept of near-wall blood “particle” behaviour, in conjunction with a wall shear stress based hemodynamic parameter, can be compared to similar histopathologic studies of intimal hyperplasia formation to understand better the role of hemodynamic factors in the biogenesis of the disease.     

Threshold of Microvascular Occlusion: Injury Size Defines the Thrombosis Scenario

Damage to the blood vessel triggers formation of a hemostatic plug, which is meant to prevent bleeding, yet the same phenomenon may result in a total blockade of a blood vessel by a thrombus, causing severe medical conditions. Here, we show that the physical interplay between platelet adhesion and hemodynamics in a microchannel manifests in a critical threshold behavior of a growing thrombus. Depending on the size of injury, two distinct dynamic pathways of thrombosis were found: the formation of a nonocclusive plug, if injury length does not exceed the critical value, and the total occlusion of the vessel by the thrombus otherwise. We develop a mathematical model that demonstrates that switching between these regimes occurs as a result of a saddle-node bifurcation. Our study reveals the mechanism of self-regulation of thrombosis in blood microvessels and explains experimentally observed distinctions between thrombi of different physical etiology. This also can be useful for the design of platelet-aggregation-inspired engineering solutions.