Mechanobiology of shear-induced platelet aggregation leading to occlusive arterial thrombosis: A multiscale in silico analysis

Recent advances in microfluidic technology of arterial thrombosis investigations

Microfluidic technology has emerged as a powerful tool in studying arterial thrombosis, allowing researchers to construct artificial blood vessels and replicate the hemodynamics of blood flow. This technology has led to significant advancements in understanding thrombosis and platelet adhesion and aggregation. Microfluidic models have various types and functions, and by studying the fabrication methods and working principles of microfluidic chips, applicable methods can be selected according to specific needs. The rapid development of microfluidic integrated system and modular microfluidic system makes arterial thrombosis research more diversified and automated, but its standardization still needs to be solved urgently. One key advantage of microfluidic technology is the ability to precisely control fluid flow in microchannels and to analyze platelet behavior under different shear forces and flow rates. This allows researchers to study the physiological and pathological processes of blood flow, shedding light on the underlying mechanisms of arterial thrombosis. In conclusion, microfluidic technology has revolutionized the study of arterial thrombosis by enabling the construction of artificial blood vessels and accurately reproducing hemodynamics. In the future, microfluidics will place greater emphasis on versatility and automation, holding great promise for advancing antithrombotic therapeutic and prophylactic measures.

Numerical study of ultra-large von Willebrand factor multimers in coagulopathy

An excessive von Willebrand factor (VWF) secretion, coupled with a moderate to severe deficiency of ADAMTS13 activity, serves as a linking mechanism between inflammation to thrombosis. The former facilitates platelet adhesion to the vessel wall and the latter is required to cleave VWF multimers. As a result, the ultra-large VWF (UL-VWF) multimers released by Weibel-Palade bodies remain uncleaved. In this study, using a computational model based on first principles, we quantitatively show how the uncleaved UL-VWF multimers interact with the blood cells to initiate microthrombosis. We observed that platelets first adhere to unfolded and stretched uncleaved UL-VWF multimers anchored to the microvessel wall. By the end of this initial adhesion phase, the UL-VWF multimers and platelets make a mesh-like trap in which the red blood cells increasingly accumulate to initiate a gradually growing microthrombosis. Although high-shear rate and blood flow velocity are required to activate platelets and unfold the UL-VWFs, during the initial adhesion phase, the blood velocity drastically drops after thrombosis, and as a result, the wall shear stress is elevated near UL-VWF roots, and the pressure drops up to 6 times of the healthy condition. As the time passes, these trends progressively continue until the microthrombosis fully develops and the effective size of the microthrombosis and these flow quantities remain almost constant. Our findings quantitatively demonstrate the potential role of UL-VWF in coagulopathy.

Illustrated State-of-the-Art Capsules of the ISTH 2024 Congress

Here, we present a series of illustrated capsules from the State of the Art (SOA) speakers at the 2024 International Society on Thrombosis and Haemostasis Congress in Bangkok, Thailand. This year's Congress marks the first time that the International Society on Thrombosis and Haemostasis has held its flagship scientific meeting in Southeast Asia and is the first to be organized by an international Planning Committee. The Bangkok program will feature innovative science and clinical updates from around the world, reflecting the diversity and multidisciplinary growth of our field. In these illustrated SOA capsules, you will find an exploration of novel models of thrombosis and bleeding and biomaterial discoveries that can trigger or block coagulation. Thromboinflammation is now understood to drive many disease states, and the SOA speakers cover cellular and coagulation responses to COVID-19 and other infections. The theme of crosstalk between coagulation and inflammation expands with capsules on protein S signaling, complement, and fibrinolytic inhibitors. Novel agents for hemophilia and thrombosis prevention are introduced. Challenging clinical conditions are also covered, such as inherited platelet disorders and antiphospholipid antibody syndrome. The scientific program in Bangkok will also showcase the work of clinicians and scientists from all parts of the world and chronicle real-world challenges. For example, 2 SOA capsules address the diagnosis and management of von Willebrand disease in low-income settings. Take some time to browse through these short illustrated reviews; we're sure that you'll be entertained, educated, and inspired to further explore the world of thrombosis and hemostasis.

Platelet Biorheology and Mechanobiology in Thrombosis and Hemostasis: Perspectives from Multiscale Computation

Thrombosis is the pathological clot formation under abnormal hemodynamic conditions, which can result in vascular obstruction, causing ischemic strokes and myocardial infarction. Thrombus growth under moderate to low shear (<1000 s-1) relies on platelet activation and coagulation. Thrombosis at elevated high shear rates (>10,000 s-1) is predominantly driven by unactivated platelet binding and aggregating mediated by von Willebrand factor (VWF), while platelet activation and coagulation are secondary in supporting and reinforcing the thrombus. Given the molecular and cellular level information it can access, multiscale computational modeling informed by biology can provide new pathophysiological mechanisms that are otherwise not accessible experimentally, holding promise for novel first-principle-based therapeutics. In this review, we summarize the key aspects of platelet biorheology and mechanobiology, focusing on the molecular and cellular scale events and how they build up to thrombosis through platelet adhesion and aggregation in the presence or absence of platelet activation. In particular, we highlight recent advancements in multiscale modeling of platelet biorheology and mechanobiology and how they can lead to the better prediction and quantification of thrombus formation, exemplifying the exciting paradigm of digital medicine.

Influence of Hematocrit Level and Integrin αIIbβIII Function on vWF-Mediated Platelet Adhesion and Shear-Induced Platelet Aggregation in a Sudden Expansion

Purpose

Shear-mediated thrombosis is a clinically relevant phenomenon that underlies excessive arterial thrombosis and device-induced thrombosis. Red blood cells are known to mechanically contribute to physiological hemostasis through margination of platelets and vWF, facilitating the unfurling of vWF multimers, and increasing the fraction of thrombus-contacting platelets. Shear also plays a role in this phenomenon, increasing both the degree of margination and the near-wall forces experienced by vWF and platelets leading to unfurling and activation. Despite this, the contribution of red blood cells in shear-induced platelet aggregation has not been fully investigated-specifically the effect of elevated hematocrit has not yet been demonstrated.

Methods

Here, a microfluidic model of a sudden expansion is presented as a platform for investigating platelet adhesion at hematocrits ranging from 0 to 60% and shear rates ranging from 1000 to 10,000 s-1. The sudden expansion geometry models nonphysiological flow separation characteristic to mechanical circulatory support devices, and the validatory framework of the FDA benchmark nozzle. PDMS microchannels were fabricated and coated with human collagen. Platelets were fluorescently tagged, and blood was reconstituted at variable hematocrit prior to perfusion experiments. Integrin function of selected blood samples was inhibited by a blocking antibody, and platelet adhesion and aggregation over the course of perfusion was monitored.

Results

Increasing shear rates at physiological and elevated hematocrit levels facilitate robust platelet adhesion and formation of large aggregates. Shear-induced platelet aggregation is demonstrated to be dependent on both αIIbβIII function and the presence of red blood cells. Inhibition of αIIbβIII results in an 86.4% reduction in overall platelet adhesion and an 85.7% reduction in thrombus size at 20-60% hematocrit. Hematocrit levels of 20% are inadequate for effective platelet margination and subsequent vWF tethering, resulting in notable decreases in platelet adhesion at 5000 and 10,000 s-1 compared to 40% and 60%. Inhibition of αIIbβIII triggered dramatic reductions in overall thrombus coverage and large aggregate formation. Stability of platelets tethered by vWF are demonstrated to be αIIbβIII-dependent, as adhesion of single platelets treated with A2A9, an anti-αIIbβIII blocking antibody, is transient and did not lead to sustained thrombus formation.

Conclusions

This study highlights driving factors in vWF-mediated platelet adhesion that are relevant to clinical suppression of shear-induced thrombosis and in vitro assays of platelet adhesion. Primarily, increasing hematocrit promotes platelet margination, permitting shear-induced platelet aggregation through αIIbβIII-mediated adhesion at supraphysiological shear rates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12195-024-00796-0.

Evaluating the impact of transient shear stress on platelet activation, adhesion, and aggregation with microfluidic chip technique

BACKGROUND

When nonphysiological stenosis occurs, the transient high shear stress formed in vessels increases the risk of thrombosis and is a potential factor for cardiovascular diseases. But the platelet adhesion and aggregation behavior at nonphysiological post-stenosis and its affecting factors are not fully understood yet.

METHODS

In this experiment, platelet aggregation on collagen and fibrinogen at different shear stresses and different hematocrits were observed by microfluidic technology. Platelet activation (P-selectin, glycoprotein IIb/IIIa) and monocyte-platelet aggregate (MPA) levels under different shear stresses were analyzed by flow cytometry.

RESULTS

On fibrinogen, platelets aggregate more at higher shear stress conditions. While on collagen, it becomes more difficult for platelets to form stable aggregation at higher shear stress conditions. If platelets adhere initially at low shear stress, stable platelet aggregation can be formed at subsequent high shear stress. Moreover, when the shear stress increases, platelet activity markers (P-selectin, glycoprotein IIb/IIIa and MPAs) increase significantly. Hematocrit affects the degree of platelet aggregation, and the influence of hematocrit is obvious at high shear stress.

CONCLUSION

Transient high shear stress (46 ms) can effectively activate platelets. Platelet aggregation behavior was different for coated fibrinogen and collagen protein. Stable platelet adhesion at post-stenosis is more dependent on fibrinogen and platelet aggregation is stable on both fibrinogen and collagen. Hematocrit can significantly affect the formation of platelet aggregation.

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.

Coarse-grained simulations of von Willebrand factor adsorption to collagen with consequent platelet recruitment

A multimeric glycoprotein of blood plasma-Von Willebrand factor (VWF)-mediates platelet adhesion to the fibrillar collagen of the subendothelial matrix if the blood vessel walls are damaged. The adsorption of VWF to collagen is thus essential for the initial stages of platelet hemostasis and thrombosis, as it plays a role of a molecular bridge between the injury and platelet adhesion receptors. Biomechanical complexity and sensitivity to the hydrodynamics are inherent in this system, therefore, modern computational methods supplement experimental studies of biophysical and molecular mechanisms that underlie platelet adhesion and aggregation in the blood flow. In the present paper, we propose a simulation framework for the VWF-mediated platelet adhesion to a plane wall with immobilized binding sites for VWF under the action of shear flow. VWF multimers and platelets are represented in the model by particles connected by elastic bonds and immersed in a viscous continuum fluid. This work complements the scientific field by taking into account the shape of a flattened platelet, but keeping a compromise between the detail of the description and the computational complexity of the model.

Biomechanical activation of blood platelets via adhesion to von Willebrand factor studied with mesoscopic simulations

Platelet adhesion and activation are essential initial processes of arterial and microvascular hemostasis, where high hydrodynamic forces from the bloodflow impede coagulation. The process relies on von Willebrand factor (VWF)-a linear multimeric protein of blood plasma plays a pivotal role in mechanochemical regulation of shear-induced platelet aggregation (SIPA). Adhesive interactions between VWF and glycoprotein receptors GPIb are crucial for platelet recruitment under high shear stress in fluid. Recent advances in experimental studies revealed that mechanical tension on the extracellular part of GPIb may trigger a cascade of biochemical reactions in platelets leading to activation of integrins [Formula: see text] (also known as GPIIb/IIIa) and strengthening of the adhesion. The present paper is aimed at investigation of this process by three-dimensional computer simulations of platelet adhesion to surface-grafted VWF multimers in pressure-driven flow of platelet-rich plasma. The simulations demonstrate that GPIb-mediated mechanotransduction is a feasible way of platelet activation and stabilization of platelet aggregates under high shear stress. Quantitative understanding of mechanochemical processes involved in SIPA would potentially promote the discovery of new anti-platelet medication and the development of multiscale numerical models of platelet thrombosis and hemostasis.

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.

SIPA in 10 milliseconds: VWF tentacles agglomerate and capture platelets under high shear

Agglomeration and capture of agglomerates after travelling a lag distance of >100 µm creates SIPA as fast as 10 milliseconds.

Phase diagrams of SIPA controlled by VWF length and concentration provide mechanistic insights for various thrombotic and hemostatic events.

Shear-induced platelet aggregation (SIPA) occurs under elevated shear rates (10 000 s⁻¹) found in stenotic coronary and carotid arteries. The pathologically high shear environment can lead to occlusive thrombosis by SIPA from the interaction of nonactivated platelets and von Willebrand factor (VWF) via glycoprotein Ib–A1 binding. This process under high shear rates is difficult to visualize experimentally with concurrent molecular- and cellular-resolutions. To understand this fast bonding, we employ a validated multiscale in silico model incorporating measured molecular kinetics and a thrombosis-on-a-chip device to delineate the flow-mediated biophysics of VWF and platelets assembly into mural microthrombi. We show that SIPA begins with VWF elongation, followed by agglomeration of platelets in the flow by soluble VWF entanglement before mural capture of the agglomerate by immobilized VWF. The entire SIPA process occurs on the order of 10 milliseconds with the agglomerate traveling a lag distance of a few hundred microns before capture, matching in vitro results. Increasing soluble VWF concentration by ∼20 times in silico leads to a ∼2 to 3 times increase in SIPA rates, matching the increase in occlusion rates found in vitro. The morphology of mural aggregates is primarily controlled by VWF molecular weight (length), where normal-length VWF leads to cluster or elongated aggregates and ultra-long VWF leads to loose aggregates seen by others’ experiments. Finally, we present phase diagrams of SIPA, which provides biomechanistic rationales for a variety of thrombotic and hemostatic events in terms of platelet agglomeration and capture.

Computational modeling of biomechanics and biorheology of heated red blood cells

Because of their compromised deformability, heat denatured erythrocytes have been used as labeled probes to visualize spleen tissue or to assess the ability of the spleen to retain stiff red blood cells (RBCs) for over three decades, e.g., see Looareesuwan et al. N. Engl. J. Med. (1987). Despite their good accessibility, it is still an open question how heated RBCs compare to certain diseased RBCs in terms of their biomechanical and biorheological responses, which may undermine their effective usage and even lead to misleading experimental observations. To help answering this question, we perform a systematic computational study of the hemorheological properties of heated RBCs with several physiologically relevant static and hemodynamic settings, including optical-tweezers test, relaxation of prestretched RBCs, RBC traversal through a capillary-like channel and a spleen-like slit, and a viscometric rheology test. We show that our in silico RBC models agree well with existing experiments. Moreover, under static tests, heated RBCs exhibit deformability deterioration comparable to certain disease-impaired RBCs such as those in malaria. For RBC traversal under confinement (through microchannel or slit), heated RBCs show prolonged transit time or retention depending on the level of confinement and heating procedure, suggesting that carefully heat-treated RBCs may be useful for studying splenic- or vaso-occlusion in vascular pathologies. For the rheology test, we expand the existing bulk viscosity data of heated RBCs to a wider range of shear rates (1-1000 s^-1) to represent most pathophysiological conditions in macro- or microcirculation. Although heated RBC suspension shows elevated viscosity comparable to certain diseased RBC suspensions under relatively high shear rates (100-1000 s^-1), they underestimate the elevated viscosity (e.g., in sickle cell anemia) at low shear rates (<10 s^-1). Our work provides mechanistic rationale for selective usage of heated RBC as a potentially useful model for studying the abnormal traversal dynamics and hemorheology in certain blood disorders.

Point of care whole blood microfluidics for detecting and managing thrombotic and bleeding risks

Point-of-care diagnostics of platelet and coagulation function present demanding challenges. Current clinical diagnostics often use centrifuged plasmas or platelets and frozen plasma standards, recombinant protein standards, or even venoms. Almost all commercialized tests of blood do not recreate the in vivo hemodynamics where platelets accumulate to high densities and thrombin is generated from a procoagulant surface. Despite numerous drugs that target platelets, insufficient coagulation, or excess coagulation, POC blood testing is essentially limited to viscoelastic methods that provide a clotting time, clot strength, and clot lysis, while used mostly in trauma centers with specialized capabilities. Microfluidics now allows small volumes of whole blood (<1 mL) to be tested under venous or arterial shear rates with multi-color readouts to follow platelet function, thrombin generation, fibrin production, and clot stability. Injection molded chips containing pre-patterned fibrillar collagen and lipidated tissue factor can be stored dry for 6 months at 4C, thus allowing rapid blood testing on single-use disposable chips. Using only a small imaging microscope and micropump, these microfluidic devices can detect platelet inhibitors, direct oral anticoagulants (DOACs) and their reversal agents. POC microfluidics are ideal for neonatal surgical applications that involve small blood samples, rapid DOAC testing in stroke or bleeding or emergency surgery situations with patients presenting high risk cofactors for either bleeding or thrombosis.

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.

Computational studies aiming to model thrombosis often rely on an arterial wall injury. Collagen and other extracellular matrix components are exposed to the bloodstream, which facilitates platelet adhesion to the wall and subsequent clot formation. However, these models are not adequate to explain thrombosis in other settings where even in the absence of a focal lesion, clots may still form under certain flow conditions. Coronary artery aneurysm thrombosis following KD is an example of the need to understand the mechanisms of thrombus initiation in the absence of an extrinsic factor. This study provides a new framework to investigate thrombus initiation in KD from a patient-specific perspective, which integrates fluid mechanics and biochemistry and which could help quantify the pro-coagulant environment induced by the aneurysm and become a predictive tool. The work presented here has broad relevance to other clinical situations where flow stagnation and transport are driving factors in thrombus formation.

A Continuum Model for the Unfolding of von Willebrand Factor

von Willebrand Factor is a mechano-sensitive protein circulating in blood that mediates platelet adhesion to subendothelial collagen and platelet aggregation at high shear rates. Its hemostatic function and thrombogenic effect, as well as susceptibility to enzymatic cleavage, are regulated by a conformational change from a collapsed globular state to a stretched state. Therefore, it is essential to account for the conformation of the vWF multimers when modeling vWF-mediated thrombosis or vWF degradation. We introduce a continuum model of vWF unfolding that is developed within the framework of our multi-constituent model of platelet-mediated thrombosis. The model considers two interconvertible vWF species corresponding to the collapsed and stretched conformational states. vWF unfolding takes place via two regimes: tumbling in simple shear and strong unfolding in flows with dominant extensional component. These two regimes were demonstrated in a Couette flow between parallel plates and an extensional flow in a cross-slot geometry. The vWF unfolding model was then verified in several microfluidic systems designed for inducing high-shear vWF-mediated thrombosis and screening for von Willebrand Disease. The model predicted high concentration of stretched vWF in key regions where occlusive thrombosis was observed experimentally. Strong unfolding caused by the extensional flow was limited to the center axis or middle plane of the channels, whereas vWF unfolding near the channel walls relied upon the shear tumbling mechanism. The continuum model of vWF unfolding presented in this work can be employed in numerical simulations of vWF-mediated thrombosis or vWF degradation in complex geometries. However, extending the model to 3-D arbitrary flows and turbulent flows will pose considerable challenges.