Simulation of thrombosis in a stenotic microchannel: The effects of vWF-enhanced shear activation of platelets

Simulation of the effect of hemolysis on thrombosis in blood-contacting medical devices

Heart failure, broadly characterized by the gradual decline of the ability of the heart to maintain adequate blood flow throughout the body's vascular network of veins and arteries, is one of the leading causes of death worldwide. Mechanical Circulatory Support is one of the few available alternative interventions for late-stage heart failure with reduced ejection fraction. A ventricular assist device is surgically implanted and connected to the left and or right heart ventricles to provide additional bloodflow, off-loading the work required by the heart to maintain circulation. Modern mechanical circulatory support devices generate non-physiological flow conditions that can lead to the damage and rupture of blood cells (hemolysis), and the formation of blood clots (thrombosis), which pose severe health risks to the patient. It is essential to improve prediction tools for blood damage to reduce the risk of hemolysis and thrombosis. A simulation-based approach examines the interaction between hemolysis and thrombosis. Incompressible finite-volume computational fluid dynamics simulations are executed on an open-hub axial flow ventricular assist device. A continuum model of thrombosis and the intrinsic coagulation process is extended to include the effect of hemolysis. The model accounts for the effect of activation of platelets by shear stress, paracrine signaling, adhesion, and hemoglobin and ADP released during hemolysis. The effect of hemolysis with thrombosis is modelled by accounting for the hyper-adhesivity of von-Willebrand Factor on extracellular hemoglobin, and the increased rate of platelet activation induced by ADP release. Thrombosis is assessed at varying inflow rates and rotor speeds, and cases are executed where thrombosis is affected by ADP release and Hb-induced hyper-adhesivity. It is found that there is a non-negligible effect from hemolysis on thrombosis across a range of rotor speeds, and that hyperadhesivity plays a dominant role in thrombus formation in the presence of hemolysis.

A continuum model for the elongation and orientation of Von Willebrand factor with applications in arterial flow

The blood protein Von Willebrand factor (VWF) is critical in facilitating arterial thrombosis. At pathologically high shear rates, the protein unfolds and binds to the arterial wall, enabling the rapid deposition of platelets from the blood. We present a novel continuum model for VWF dynamics in flow based on a modified viscoelastic fluid model that incorporates a single constitutive relation to describe the propensity of VWF to unfold as a function of the scalar shear rate. Using experimental data of VWF unfolding in pure shear flow, we fix the parameters for VWF's unfolding propensity and the maximum VWF length, so that the protein is half unfolded at a shear rate of approximately 5000 s - 1 . We then use the theoretical model to predict VWF's behaviour in two complex flows where experimental data are challenging to obtain: pure elongational flow and stenotic arterial flow. In pure elongational flow, our model predicts that VWF is 50% unfolded at approximately 2000 s - 1 , matching the established hypothesis that VWF unfolds at lower shear rates in elongational flow than in shear flow. We demonstrate the sensitivity of this elongational flow prediction to the value of maximum VWF length used in the model, which varies significantly across experimental studies, predicting that VWF can unfold between 2000 and 3200 s - 1 depending on the selected value. Finally, we examine VWF dynamics in a range of idealised arterial stenoses, predicting the relative extension of VWF in elongational flow structures in the centre of the artery compared to high shear regions near the arterial walls.

A homogenized two-phase computational framework for meso- and macroscale blood flow simulations

BACKGROUND AND OBJECTIVE

Owing to the complexity of physics linked with blood flow and its associated phenomena, appropriate modeling of the multi-constituent rheology of blood is of primary importance. To this effect, various kinds of computational fluid dynamic models have been developed, each with merits and limitations. However, when additional physics like thrombosis and embolization is included within the framework of these models, computationally efficient scalable translation becomes very difficult. Therefore, this paper presents a homogenized two-phase blood flow framework with similar characteristics to a single fluid model but retains the flow resolution of a classical two-fluid model. The presented framework is validated against four different sets of experiments.

METHODS

The two-phase model of blood presented here is based on the classical diffusion-flux framework. Diffusion flux models are known to be less computationally expensive than two-fluid multiphase models since the numerical implementation resembles single-phase flow models. Diffusion flux models typically use empirical slip velocity correlations to resolve the motion between phases. However, such correlations do not exist for blood. Therefore, a modified slip velocity equation is proposed, derived rigorously from the two-fluid governing equations. An additional drag law for red blood cells (RBCs) as a function of volume fraction is evaluated using a previously published cell-resolved solver. A new hematocrit-dependent expression for lift force on RBCs is proposed. The final governing equations are discretized and solved using the open-source software OpenFOAM.

RESULTS

The framework is validated against four sets of experiments: (i) flow through a rectangular microchannel to validate RBC velocity profiles against experimental measurements and compare computed hematocrit distributions against previously reported simulation results (ii) flow through a sudden expansion microchannel for comparing experimentally obtained contours of hematocrit distributions and normalized cell-free region length obtained at different flowrates and inlet hematocrits, (iii) flow through two hyperbolic channels to evaluate model predictions of cell-free layer thickness, and (iv) flow through a microchannel that mimics crevices of a left ventricular assist device to predict hematocrit distributions observed experimentally. The simulation results exhibit good agreement with the results of all four experiments.

CONCLUSION

The computational framework presented in this paper has the advantage of resolving the multiscale physics of blood flow while still leveraging numerical techniques used for solving single-phase flows. Therefore, it becomes an excellent candidate for addressing more complicated problems related to blood flow, such as modeling mechanical entrapment of RBCs within blood clots, predicting thrombus composition, and visualizing clot embolization.

Shear Histories Alter Local Shear Effects on Thrombus Nucleation and Growth

Our goal was to determine the impact of physiological and pathological shear histories on platelet nucleation and thrombus growth at various local shear rates. We designed and characterized a microfluidic device capable of subjecting platelets to shear histories reaching as high as 6700 s- 1 in a single passage. Time-lapse videos of platelets and thrombi are captured using fluorescence microscopy. Thrombi are tracked, and the degree of thrombosis is evaluated through surface coverage, platelet nucleation maps, and ensemble-averaged aggregate areas and intensities. Surface coverage rates were the lowest when platelets deposited at high shear rates following a pathological shear history and were highest at low shear rates following a pathological shear history. Early aggregate area growth rates were significantly larger for thrombi developing at high shear following physiological shear history than at high shear following a pathological shear history. Aggregate vertical growth was restricted when depositing at low shear following a pathological shear history. In contrast, thrombi grew faster vertically following physiological shear histories. These results show that physiological shear histories pose thrombotic risks via volumetric growth, and pathological shear histories drastically promote nucleation. These findings may inform region-based geometries for biomedical devices and refine thrombosis simulations.

Multi-Method Investigation of Blood Damage Induced By Blood Pumps in Different Clinical Support Modes

To investigate the effects of blood pumps operated in different modes on nonphysiologic flow patterns, cell and protein function, and the risk of bleeding, thrombosis, and hemolysis, an extracorporeal blood pump (CentriMag) was operated in three clinical modalities including heart failure (HF), venous-venous (V-V) extracorporeal membrane oxygenation (ECMO), and venous-arterial (V-A) ECMO. Computational fluid dynamics (CFD) methods and coupled hemolysis models as well as recently developed bleeding and thrombosis models associated with changes in platelet and von Willebrand factor (vWF) function were used to predict hydraulic performance and hemocompatibility. The V-A ECMO mode had the highest flow losses and shear stress levels, the V-V ECMO mode was intermediate, and the HF mode was the lowest. Different nonphysiologic flow patterns altered cell/protein morphology and function. The V-A ECMO mode resulted in the highest levels of platelet activation, receptor shedding, vWF unfolding, and high molecular weight multimers vWF (HMWM-vWF) degradation, leading to the lowest platelet adhesion and the highest vWF binding capacity, intermediate in the V-V ECMO mode, and opposite in the HF mode. The V-A ECMO mode resulted in the highest risk of bleeding, thrombosis, and hemolysis, with the V-V ECMO mode intermediate and the HF mode lowest. These findings are supported by published experimental or clinical statistics. Further studies found that secondary blood flow passages resulted in the highest risk of blood damage. Nonphysiologic blood flow patterns were strongly associated with cell and protein function changing, blood damage, and complications.

Development and validation of a mathematical model for evaluating shear-induced damage of von Willebrand factor

PURPOSE

To develop a mathematical model for predicting shear-induced von Willebrand factor (vWF) function modification which can be used to guide ventricular assist devices (VADs) design, and evaluate the damage of high molecular weight multimers (HMWM)-vWF in VAD patients for reducing clinical complications.

METHODS

Mathematical models were constructed based on three morphological variations (globular vWF, unfolded vWF and degraded vWF) of vWF under shear stress conditions, in which parameters were obtained from previous studies or fitted by experimental data. Different clinical support modes (pediatric vs. adult mode), different VAD operating states (pulsation vs. constant mode) and different clinical VADs (HeartMate II, HeartWare and CentriMag) were utilized to analyze shear-induced damage of HMWM-vWF based on our vWF model. The accuracy and feasibility of the models were evaluated using various experimental and clinical cases, and the biomechanical mechanisms of HMWM-vWF degradation induced by VADs were further explained.

RESULTS

The mathematical model developed in this study predicted VAD-induced HMWM-vWF degradation with high accuracy (correlation with experimental data r2 > 0.99). The numerical results showed that VAD in the pediatric mode resulted in more HMWM-vWF degradation per unit time and per unit flow rate than in the adult mode. However, the total degradation of HMWM-vWF is less in the pediatric mode than in the adult mode because the pediatric mode has fewer times of blood circulation than the adult mode in the same amount of time. The ratio of HMWM-vWF degradation was lower in the pulsation mode than in the constant mode. This is due to the increased flushing of VADs in the pulsation mode, which avoids prolonged stagnation of blood in high shear regions. This study also found that the design feature, rotor size and volume of the VADs, and the superimposed regions of high shear stress and long residence time inside VADs affect the degradation of HMWM-vWF. The axial flow VADs (HeartMate II) showed higher degradation of HMWM-vWF compared to centrifugal VADs (HeartWare and CentriMag). Compared to fully magnetically suspended VADs (CentriMag), hydrodynamic suspended VADs (HeartWare) produced extremely high degradation of HWMW-vWF in its narrow hydrodynamic clearance. Finally, the study used a mathematical model of HMWM-vWF degradation to interpret the clinical statistics from a biomechanical perspective and found that minimizing the rotating speed of VADs within reasonable limits helps to reduce HWMW-vWF degradation. All predicted conclusions are supported by the experimental and clinical data.

CONCLUSION

This study provides a validated mathematical model to assess the shear-induced degradation of HMWM-vWF, which can help to evaluate the damage of HMWM-vWF in patients implanted with VADs for reducing clinical complications, and to guide the optimization of VADs for improving hemocompatibility.

Bio-inspired microfluidics: A review

Biomicrofluidics, a subdomain of microfluidics, has been inspired by several ideas from nature. However, while the basic inspiration for the same may be drawn from the living world, the translation of all relevant essential functionalities to an artificially engineered framework does not remain trivial. Here, we review the recent progress in bio-inspired microfluidic systems via harnessing the integration of experimental and simulation tools delving into the interface of engineering and biology. Development of "on-chip" technologies as well as their multifarious applications is subsequently discussed, accompanying the relevant advancements in materials and fabrication technology. Pointers toward new directions in research, including an amalgamated fusion of data-driven modeling (such as artificial intelligence and machine learning) and physics-based paradigm, to come up with a human physiological replica on a synthetic bio-chip with due accounting of personalized features, are suggested. These are likely to facilitate physiologically replicating disease modeling on an artificially engineered biochip as well as advance drug development and screening in an expedited route with the minimization of animal and human trials.

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.

Study on the antithrombotic effect and physiological mechanism of okanin

The objectives of this study were to investigate the antithrombotic effect and physiological mechanism of okanin, a flavonoid monomer in Coreopsis tinctoria Nutt. The antithrombotic effects of okanin were determined by the anticoagulant activity test in vitro and in vivo, the venous thrombosis and arterial thrombosis test in rats. To study the antithrombotic physiological mechanisms of okanin, UV spectrophotometer and enzyme-linked immunosorbent assay (ELISA) were used to determine the effects of three concentrations of okanin on the contents of 6-keto-prostaglandin F1α (6-Keto-PGF_1α), thromboxane B₂ (TXB₂), endothelin-1 (ET-1), antithrombin III (AT-Ⅲ), protein C (PC) and von willebrand factor (vWF) in the plasma of rats with arterial thrombosis; ELISA was used to detect the effects of okanin on the contents of plasminogen (PLG), tissue plasminogen activator (t-PA) and type-1 plasminogen activator inhibitor (PAI-1) in the plasma of mice and Chinese white rabbits. The results showed that okanin could prolong the coagulation time in vitro and in vivo of animals (P < 0.01 in the high dose group) and the activated partial thromboplastin time (APTT), prothrombin time (PT), and thrombin time (TT) of human venous blood (ATPP of medium dose group P < 0.01; PT, TT P < 0.05. P < 0.01 in the high dose group); inhibit the maximum platelet aggregation rate of rabbits (P < 0.05 in the low dose group; P < 0.01 in the medium and high dose groups), decrease the dry and wet weight of venous thrombosis and the wet weight of common carotid artery thrombosis in rats (low dose group, P < 0.05; medium and high dose groups, P < 0.01); increase the levels of 6-Keto-PGF_1α, AT-Ⅲ, PLG and t-PA in animal plasma; decrease the levels of TXB₂, ET-1, vWF and PAI-₁ in animal plasma. It is concluded that okanin can significantly inhibit thrombosis, and its physiological mechanisms were related to affecting the activation of related coagulation factors in endogenous and exogenous coagulation pathways, affecting the physiological characteristics of platelets, repairing damaged vascular endothelial cells and enhancing the activity of the fibrinolytic system.

Models of Shear-Induced Platelet Activation and Numerical Implementation With Computational Fluid Dynamics Approaches

Shear-induced platelet activation is one of the critical outcomes when blood is exposed to elevated shear stress. Excessively activated platelets in the circulation can lead to thrombus formation and platelet consumption, resulting in serious adverse events such as thromboembolism and bleeding. While experimental observations reveal that it is related to the shear stress level and exposure time, the underlying mechanism of shear-induced platelet activation is not fully understood. Various models have been proposed to relate shear stress levels to platelet activation, yet most are modified from the empirically calibrated power-law model. Newly developed multiscale platelet models are tested as a promising approach to capture a single platelet's dynamic shape during activation, but it would be computationally expensive to employ it for a large-scale analysis. This paper summarizes the current numerical models used to study the shear-induced platelet activation and their computational applications in the risk assessment of a particular flow pattern and clot formation prediction.

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.

The excessive formation of blood clots under flow within the circulatory system (thrombosis) is known to initiate heart attacks and strokes. Therefore, obtaining insights into the formation and progression of these clots will be useful in evaluating pharmacological options. To this end, we have developed a 3D computational model that tracks the growth of a blood clot under flow from initial platelet deposition to full vessel occlusion in the presence of soluble platelet agonists. We first validated the model against experimental predictions of blood clots formed in vitro. Due to the construction of the model in 3D, we were able to carry out simulations of clot formation under important clinical situations, namely cylindrical blood vessels and acute flow narrowings (stenoses). To our knowledge, our model is the first of its kind that can account for patient-specific platelet phenotypes to perform robust 3D simulations of thrombus growth in geometries of clinical relevance.

Mathematical and Computational Modeling of Device-Induced Thrombosis

Given the extensive and routine use of cardiovascular devices, a major limiting factor to their success is the thrombotic rate that occurs. This both poses direct risk to the patient and requires counterbalancing with anticoagulation and other treatment strategies, contributing additional risks. Developing a better understanding of the mechanisms of device-induced thrombosis to aid in device design and medical management of patients is critical to advance the ubiquitous use and durability. Thus, mathematical and computational modelling of device-induced thrombosis has received significant attention recently, but challenges remain. Additional areas that need to be explored include microscopic/macroscopic approaches, reconciling physical and numerical timescales, immune/inflammatory responses, experimental validation, and incorporating pathologies and blood conditions. Addressing these areas will provide engineers and clinicians the tools to provide safe and effective cardiovascular devices.

An integrated fluid-chemical model toward modeling the thrombus formation in an idealized model of aortic dissection

Type B aortic dissection is a major aortic catastrophe that can be acutely complicated by rapid expansion, rupture, and malperfusion syndromes. The separation of the intima from aortic walls will form a second blood-filled lumen defined as "false lumen (FL)", where the thrombus is more likely to form due to the local stasis hemodynamic conditions. Complete thrombosis of FL is associated with a beneficial outcome while patency and partial thrombosis will lead to later complications. However, the thrombosis mechanism is still unclear and little is known about the impact of chemical species transported by blood flow on this process. The proteins involved in the coagulation cascade (CC) may play an important role in the process of thrombosis, especially in the activation and stabilization of platelets. Based on this hypothesis, a reduced-order fluid-chemical model was established to simulate CC in an aortic dissection phantom with two tears. A high level of fibrin is continuously observed at the top of the FL and some time-varying areas between two tears, indicating a high likelihood of thrombus formation there. This finding is consistent with the clinical observation. The time evolution of coagulation factors is greatly affected by local hemodynamics, especially in the high disturbance zone where the evolution has characteristics of periodic changes consistent with the flow field. The ability of the proposed model to reproduce the CC response provides a potential application to integrate with a model that can simulate platelet activities, forming a biochemical-based model which would help unveil the mechanisms of thrombosis in FL and the clinical decision of appropriate treatment.

Modelling Particle Agglomeration on through Elastic Valves under Flow

This work proposes a model of particle agglomeration in elastic valves replicating the geometry and the fluid dynamics of a venous valve. The fluid dynamics is simulated with Smooth Particle Hydrodynamics, the elastic leaflets of the valve with the Lattice Spring Model, while agglomeration is modelled with a 4-2 Lennard-Jones potential. All the models are combined together within a single Discrete Multiphysics framework. The results show that particle agglomeration occurs near the leaflets, supporting the hypothesis, proposed in previous experimental work, that clot formation in deep venous thrombosis is driven by the fluid dynamics in the valve.

Influence of shear rate and surface chemistry on thrombus formation in micro-crevice

Thromboembolic complications remain a central issue in management of patients on mechanical circulatory support. Despite the best practices employed in design and manufacturing of modern ventricular assist devices, complexity and modular nature of these systems often introduces internal steps and crevices in the flow path which can serve as nidus for thrombus formation. Thrombotic potential is influenced by multiple factors including the characteristics of the flow and surface chemistry of the biomaterial. This study explored these elements in the setting of blood flow over a micro-crevice using a multi-constituent numerical model of thrombosis. The simulations reproduced the platelet deposition patterns observed experimentally and elucidated the role of flow, shear rate, and surface chemistry in shaping the deposition. The results offer insights for design and operation of blood-contacting devices.

Platelet activation via dynamic conformational changes of von Willebrand factor under shear

Shear-induced conformational changes of von Willebrand factor (VWF) play an important role in platelet activation. A novel approach describing VWF unfolding on the platelet surface under dynamic shear stress is proposed. Cumulative effect of dynamic shear on platelet activation via conformational changes of VWF is analysed. The critical condition of shear-induced platelet activation is formulated. The explicit expression for the threshold value of cumulative shear stress as a function of VWF multimer size is derived. The results open novel prospects for pharmacological regulation of shear-induced platelet activation through control of VWF multimers size distribution.