A Predictive Model of High Shear Thrombus Growth

N-Acetyl Cysteine Prevents Arterial Thrombosis in a Dose-Dependent Manner In Vitro and in Mice

BACKGROUND

Platelet-rich thrombi occlude arteries causing fatal infarcts like heart attacks and strokes. Prevention of thrombi by current antiplatelet agents can cause major bleeding. Instead, we propose using N-acetyl cysteine (NAC) to act against the protein VWF (von Willebrand factor), and not platelets, to prevent arterial thrombi from forming.

METHODS

NAC was assessed for its ability to prevent arterial thrombosis by measuring platelet accumulation rate and occlusion time using a microfluidic model of arterial thrombosis with human blood. Acute clot formation, clot stability, and tail bleeding were measured in vivo with the murine modified Folts model. The effect of NAC in the murine model after 6 hours was also measured to determine any persistent effects of NAC after it has been cleared from the blood.

RESULTS

We demonstrate reduction of thrombi formation following treatment with NAC in vitro and in vivo. Human whole blood treated with 3 or 5 mmol/L NAC showed delayed thrombus formation 2.0× and 3.7× longer than control, respectively (P<0.001). Blood treated with 10 mmol/L NAC did not form an occlusive clot, and no macroscopic platelet aggregation was visible (P<0.001). In vivo, a 400-mg/kg dose of NAC prevented occlusive clots from forming in mice without significantly affecting tail bleeding times. A lower dose of NAC significantly reduced clot stability. Mice given multiple injections showed that NAC has a lasting and cumulative effect on clot stability, even after being cleared from the blood (P<0.001).

CONCLUSIONS

Both preclinical models demonstrate that NAC prevents thrombus formation in a dose-dependent manner without significantly affecting bleeding time. This work highlights a new pathway for preventing arterial thrombosis, different from antiplatelet agents, using an amino acid derivative as an antithrombotic therapeutic.

Uncertainty quantification of a thrombosis model considering the clotting assay PFA-100®

Mathematical models of thrombosis are currently used to study clinical scenarios of pathological thrombus formation. As these models become more complex to predict thrombus formation dynamics high computational cost must be alleviated and inherent uncertainties must be assessed. Evaluating model uncertainties allows to increase the confidence in model predictions and identify avenues of improvement for both thrombosis modeling and anti-platelet therapies. In this work, an uncertainty quantification analysis of a multi-constituent thrombosis model is performed considering a common assay for platelet function (PFA-100®). The analysis is facilitated thanks to time-evolving polynomial chaos expansions used as a parametric surrogate for the full thrombosis model considering two quantities of interest; namely, thrombus volume and occlusion percentage. The surrogate is thoroughly validated and provides a straightforward access to a global sensitivity analysis via computation of Sobol' coefficients. Six out of 15 parameters linked to thrombus consitution, vWF activity, and platelet adhesion dynamics were found to be most influential in the simulation variability considering only individual effects; while parameter interactions are highlighted when considering the total Sobol' indices. The influential parameters are related to thrombus constitution, vWF activity, and platelet to platelet adhesion dynamics. The surrogate model allowed to predict realistic PFA-100® closure times of 300,000 virtual cases that followed the trends observed in clinical data. The current methodology could be used including common anti-platelet therapies to identify scenarios that preserve the hematological balance.

Structure of shear-induced platelet aggregated clot formed in an in vitro arterial thrombosis model

The structure of occlusive arterial thrombi is described in this paper. Macroscopic thrombi are made from whole blood in a collagen-coated, large-scale stenosis with high shear flow conditions similar to an atherosclerotic artery. The millimeter-sized thrombi are harvested for histology and scanning electron microscopy. Histological images showed 3 distinctive structures of the thrombus. A) The upstream region showed string-like platelet aggregates growing out from the wall to protrude into the central lumen, while RBCs were trapped between the strings. The strings were >10x as long as they were wide and reached out to join the strings from the opposite wall. B) Near the apex, the platelet strings coalesced into a dense mass with microchannels that effectively occludes the lumen. C) In the expansion region, the thrombus ended abruptly with an annulus of free blood in the flow separation zone. Scanning electron microscopy showed dense clusters of spherical platelets upstream and downstream, with amorphous platelets in the occluded throat consistent with prior activation. The total clot is estimated to contain 1.23 billion platelets with pores on the order of 10-100 microns. The results reveal a complex structure of arterial thrombi that grow from their tips under high shear stress to bridge the 2.5 mm lumen quickly with VWF-platelet strings. The occlusion leaves many microchannels that allow some flow through the bulk of the thrombus. This architecture can create occlusion or hemostasis rapidly with minimal material, yet remain porous for potential delivery of lytic agents to the core of the thrombus.

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.

Thrombosis and Hemodynamics: external and intrathrombus gradients

Distinct from dilute, isotropic, and homogeneous reaction systems typically used in laboratory kinetic assays, blood is concentrated, two-phase, flowing, and highly anisotropic when clotting on a surface. This review focuses on spatial gradients that are generated and can dictate thrombus structure and function. Novel experimental and computational tools have recently emerged to explore reaction-transport coupling during clotting. Multiscale simulations help bridge tissue length scales (the coronary arteries) to millimeter scales of a growing clot to the microscopic scale of single-cell signaling and adhesion. Microfluidic devices help create and control pathological velocity profiles, albeit at a low Reynolds number. Since rate processes and force loading are often coupled, this review highlights prevailing convective-diffusive transport physics that modulate cellular and molecular processes during thrombus formation.

Global Thrombosis Test: Occlusion by Coagulation or SIPA?

The global thrombosis test (GTT) is a point of care device that tests thrombotic and thrombolytic status. The device exposes whole blood flow to a combination of both high and low shear stress past and between ball bearings potentially causing thrombin and fibrin formation. The question arises as to whether thrombosis in the GTT is dominated by coagulation-triggered red clot or high shear-induced white clot. We investigated the nature of the thrombus formed in the GTT, the device efficacy, human factors use, and limitations. The GTT formed clots that were histologically fibrin-rich with trapped red blood cells. The occlusion time (OT) was more consistent with coagulation than high shear white clot and was strongly lengthened by heparin and citrate, two common anticoagulants. The clot was lysed by tissue plasminogen activator (tPA), also consistent with a fibrin-rich red clot. Changing the bead to a collagen-coated surface and eliminating the low shear zone between the beads induced a rapid OT consistent with a platelet-rich thrombus that was relatively resistant to heparin or tPA. The evidence points to the GTT as occluding primarily due to fibrin-rich red clot from coagulation rather than high shear platelet aggregation and occlusion associated with arterial thrombosis.

Blood Flow Velocimetry in a Microchannel During Coagulation Using Particle Image Velocimetry and Wavelet-Based Optical Flow Velocimetry

This article describes novel measurements of the velocity of whole blood flow in a microchannel during coagulation. The blood is imaged volumetrically using a simple optical setup involving a white light source and a microscope camera. The images are processed using particle image velocimetry (PIV) and wavelet-based optical flow velocimetry (wOFV), both of which use images of individual blood cells as flow tracers. Measurements of several clinically relevant parameters such as the clotting time, decay rate, and blockage ratio are computed. The high-resolution wOFV results yield highly detailed information regarding thrombus formation and corresponding flow evolution that is the first of its kind.

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.

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

Occlusive thrombosis in arteries causes heart attacks and strokes. The rapid growth of thrombus at elevated shear rates (~10,000 1/s) relies on shear-induced platelet aggregation (SIPA) thought to come about from the entanglement of von Willebrand factor (VWF) molecules. The mechanism for SIPA is not yet understood in terms of cell- and molecule-level dynamics in fast flowing bloodstreams. Towards this end, we develop a multiscale computational model to recreate SIPA in silico, where the suspension dynamics and interactions of individual platelets and VWF multimers are resolved directly. The platelet-VWF interaction via GP1b-A1 bonds is prescribed with intrinsic binding rates theoretically derived and informed by single-molecule measurements. The model is validated against existing microfluidic SIPA experiments, showing good agreement with the in vitro observations in terms of the morphology, traveling distance and capture time of the platelet aggregates. Particularly, the capture of aggregates can occur in a few milliseconds, comparable to the platelet transit time through pathologic arterial stenotic sections and much shorter than the time for shear-induced platelet activation. The multiscale SIPA simulator provides a cross-scale tool for exploring the biophysical mechanisms of SIPA in silico that are difficult to access with single-molecule measurements or micro-/macro-fluidic assays only.

Flow Dynamics and Mixing in Extracorporeal Support: A Study of the Return Cannula

Cannulation strategies in medical treatment such as in extracorporeal life support along with the associated cannula position, orientation and design, affects the mixing and the mechanical shear stress appearing in the flow field. This in turn influences platelet activation state and blood cell destruction. In this study, a co-flowing confined jet similar to a return cannula flow configuration found in extracorporeal membrane oxygenation was investigated experimentally. Cannula diameters, flow rate ratios between the jet and the co-flow and cannula position were studied using Particle Image Velocimetry and Planar Laser Induced Fluorescence. The jet was turbulent for all but two cases, in which a transitional regime was observed. The mixing, governed by flow entrainment, shear layer induced vortices and a backflow along the vessel wall, was found to require 9–12 cannula diameters to reach a fully homogeneous mixture. This can be compared to the 22–30 cannula diameters needed to obtain a fully developed flow. Although not significantly affecting mixing characteristics, cannula position altered the development of the flow structures, and hence the shear stress characteristics.

Lysis of arterial thrombi by perfusion of N,N'-Diacetyl-L-cystine (DiNAC)

The search persists for a safe and effective agent to lyse arterial thrombi in the event of acute heart attacks or strokes due to thrombotic occlusion. The culpable thrombi are composed either primarily of platelets and von Willebrand Factor (VWF), or polymerized fibrin, depending on the mechanism of formation. Current thrombolytics were designed to target red fibrin-rich clots, but may be not be efficacious on white VWF-platelet-rich arterial thrombi. We have developed an in vitro system to study the efficacy of known and proposed thrombolytic agents on white clots formed from whole blood in a stenosis with arterial conditions. The agents and adjuncts tested were tPA, ADAMTS-13, abciximab, N-acetyl cysteine, and N,N'-Diacetyl-L-cystine (DiNAC). Most of the agents, including tPA, had little thrombolytic effect on the white clots. In contrast, perfusion of DiNAC lysed thrombi as quickly as 1.5 min, which ranged up to 30 min at lower concentrations, and resulted in an average reduction in surface area of 71 ± 20%. The clot burden was significantly reduced compared to both tPA and a saline control (p<0.0001). We also tested the efficacy of all agents on red fibrinous clots formed in stagnant conditions. DiNAC did not lyse red clots, whereas tPA significantly lysed red clot over 48 h (p<0.01). These results lead to a novel use for DiNAC as a possible thrombolytic agent against acute arterial occlusions that could mitigate the risk of hyper-fibrinolytic bleeding.

Haemodynamic flow conditions at the initiation of high-shear platelet aggregation: a combined in vitro and cellular in silico study

The influence of the flow environment on platelet aggregation is not fully understood in high-shear thrombosis. The objective of this study is to investigate the role of a high shear rate in initial platelet aggregation. The haemodynamic conditions in a microfluidic device are studied using cell-based blood flow simulations. The results are compared with in vitro platelet aggregation experiments performed with porcine whole blood (WB) and platelet-rich-plasma (PRP). We studied whether the cell-depleted layer in combination with high shear and high platelet flux can account for the distribution of platelet aggregates. High platelet fluxes at the wall were found in silico. In WB, the platelet flux was about twice as high as in PRP. Additionally, initial platelet aggregation and occlusion were observed in vitro in the stenotic region. In PRP, the position of the occlusive thrombus was located more downstream than in WB. Furthermore, the shear rates and stresses in cell-based and continuum simulations were studied. We found that a continuum simulation is a good approximation for PRP. For WB, it cannot predict the correct values near the wall.

Establishment of an in vitro thrombogenicity test system with cyclic olefin copolymer substrate for endothelial layer formation

In vitro thrombogenicity test systems require co-cultivation of endothelial cells and platelets under blood flow-like conditions. Here, a commercially available perfusion system is explored using plasma-treated cyclic olefin copolymer (COC) as a substrate for the endothelial cell layer. COC was characterized prior to endothelialization and co-cultivation with platelets under static or flow conditions. COC exhibits a low roughness and a moderate hydrophilicity. Flow promoted endothelial cell growth and prevented platelet adherence. These findings show the suitability of COC as substrate and the importance of blood flow-like conditions for the assessment of the thrombogenic risk of drugs or cardiovascular implant materials.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1557/s43579-021-00072-6.

Biorheology of occlusive thrombi formation under high shear: in vitro growth and shrinkage

Occlusive thrombi formed under high flow shear rates develop very rapidly in arteries and may lead to myocardial infarction or stroke. Rapid platelet accumulation (RPA) and occlusion of platelet-rich thrombi and clot shrinkage have been studied after flow arrest. However, the influence of margination and shear rate on occlusive clot formation is not fully understood yet. In this study, the influence of flow on the growth and shrinkage of a clot is investigated. Whole blood (WB) and platelet-rich plasma (PRP) were perfused at high shear rates (> 3,000 s⁻¹) through two microfluidic systems with a stenotic section under constant pressure. The stenotic section of the two devices are different in stenotic length (1,000 vs 150 μm) and contraction angle of the stenosis (15° vs 80°). In all experiments, the flow chamber occluded in the stenotic section. Besides a significantly increased lag time and decreased RPA rate for PRP compared to WB (p < 0.01), the device with a shorter stenotic section and steeper contraction angle showed a shear-dependent occlusion and lag time for both PRP and WB. This shear-dependent behavior of the platelet aggregate formation might be caused by the stenotic geometry.

Platelet α-granules are required for occlusive high-shear-rate thrombosis

von Willebrand factor (VWF) is essential for the induction of arterial thrombosis. In this study, we investigated the critical role of platelet VWF in occlusive thrombosis formation at high shear in mice that do not express platelet VWF (Nbeal2-/-). Using in silico modeling, in vitro high-shear microfluidics, and an in vivo Folts model of arterial thrombosis we reproduced the platelet dynamics that occur under pathological flow in a stenosed vessel. Computational fluid dynamics (CFDs) simulated local hemodynamics in a stenosis based on arterial geometries. The model predicted shear rates, time course of platelet adhesion, and time to occlusion. These predictions were validated in vitro and in vivo. Occlusive thrombosis developed in wild-type control mice that had normal levels of plasma VWF and platelet VWF in vitro and in vivo. Occlusive thrombosis did not form in the Nbeal2-/- mice that had normal plasma VWF and an absence of platelet VWF. Occlusive thrombosis was corrected in Nbeal2-/- microfluidic assays by the addition of exogenous normal platelets with VWF. Combining model and experimental data, we demonstrated the necessary requirement of platelet VWF in α-granules in forming an occlusive thrombus under high shear. These results could inspire new pharmacological targets specific to pathological conditions and prevent arterial thrombosis.

Occlusive thrombosis in arteries

Thrombus formation in major arteries is life threatening. In this review article, we discuss how an arterial thrombus can form under pathologically high shear stresses, with bonding rates estimated to be the fastest Kon values in biochemistry. During occlusive thrombosis in arteries, the growth rate of the thrombus explodes to capture a billion platelets in about 10 min. Close to 100% of all platelets passing the thrombus are captured by long von Willebrand factor (vWF) strands that quickly form tethered nets. The nets grow in patches where shear stress is high, and the local concentration of vWF is elevated due to α-granule release by previously captured platelets. This rapidly formed thrombus has few red blood cells and so has a white appearance and is much stronger and more porous than clots formed through coagulation. Understanding and modeling the biophysics of this event can predict totally new approaches to prevent and treat heart attacks and strokes.

Identifying the start of a platelet aggregate by the shear rate and the cell-depleted layer

Computer simulations were performed to study the transport of red blood cells and platelets in high shear flows, mimicking earlier published in vitro experiments in microfluidic devices with high affinity for platelet aggregate formation. The goal is to understand and predict where thrombus formation starts. Additionally, the need of cell-based modelling in these microfluidic devices is demonstrated by comparing our results with macroscopic models, wherein blood is modelled as a continuous fluid. Hemocell, a cell-based blood flow simulation framework is used to investigate the transport physics in the microfluidic devices. The simulations show an enlarged cell-depleted layer at the site where a platelet aggregate forms in the experiments. In this enlarged cell-depleted layer, the probability to find a platelet is higher than in the rest of the microfluidic device. In addition, the shear rates are sufficiently high to allow for the von Willebrand factor to elongate in this region. We hypothesize that the enlarged cell-depleted layer combined with a sufficiently large platelet flux and sufficiently high shear rates result in an haemodynamic environment that is a preferred location for initial platelet aggregation.

Shear-induced platelet aggregation: 3D-grayscale microfluidics for repeatable and localized occlusive thrombosis

Atherothrombosis leads to complications of myocardial infarction and stroke as a result of shear-induced platelet aggregation (SIPA). Clinicians and researchers may benefit from diagnostic and benchtop microfluidic assays that assess the thrombotic activity of an individual. Currently, there are several different proposed point-of-care diagnostics and microfluidic thrombosis assays with different design parameters and end points. The microfluidic geometry, surface coatings, and anticoagulation may strongly influence the precision of these assays. Variability in selected end points also persists, leading to ambiguous results. This study aims to assess the effects of three physiologically relevant extrinsic design factors on the variability of a single end point to provide a quantified rationale for design parameter and end-point standardization. Using a design of experiments approach, we show that the methods of channel fabrication and collagen surface coating significantly impact the variability of occlusion time from porcine whole blood, while anticoagulant selection between heparin and citrate did not significantly impact the variability. No factor was determined to significantly impact the mean occlusion time within the assay. Occlusive thrombus was found to consistently form in the first third (333 μm) of the high shear zone and not in the shear gradient regions. The selection of these factors in the design of point-of-care diagnostics and experimental SIPA assays may lead to increased precision and specificity in high shear thrombosis studies.

The Margination of Particles in Areas of Constricted Blood Flow

Stroke is a leading cause of death globally and is caused by stenoses, abnormal narrowings of blood vessels. Recently, there has been an increased interest in shear-activated particle clusters for the treatment of stenosis, but there is a lack of literature investigating the impact of different stenosis geometries on particle margination. Margination refers to the movement of particles toward the blood vessel wall and is desirable for drug delivery. The current study investigated ten different geometries and their effects on margination. Microfluidic devices with a constricted area were fabricated to mimic a stenosed blood vessel with different extent of occlusion, constricted length, and eccentricity (gradualness of the constriction and expansion). Spherical fluorescent particles with a diameter of 2.11 μm were suspended in blood and tracked as they moved into, through, and out of the constricted area. A margination parameter, M, was used to quantify margination based on the particle distribution after velocity normalization. Experimental results suggested that a constriction leads to an enhanced margination, whereas an expansion is responsible for a decrease in margination. Further, margination was found to increase with increasing percent occlusion and constriction length, likely a result of higher shear rate and longer residence time, respectively. Margination decreases as the stenosis geometry becomes more gradual (eccentricity increases) with the exception of a sudden constriction/expansion geometry. The findings demonstrate the importance of geometric effects on margination and call for detailed numerical modeling and geometric characterization of the stenosed areas to fully understand the underlying physics.

Prediction of Thrombus Growth: Effect of Stenosis and Reynolds Number

Shear stresses play a major role in platelet-substrate interactions and thrombus formation and growth in blood flow, where under both pathological and physiological conditions platelet adhesion and accumulation occur. In this study, a shear-dependent continuum model for platelet activation, adhesion and aggregation is presented. The model was first verified under three different shear conditions and at two heparin levels. Three-dimensional simulations were then carried out to evaluate the performance of the model for severely damaged (stripped) aortas with mild and severe stenosis degrees in laminar flow regime. For these cases, linear shear-dependent functions were developed for platelet-surface and platelet-platelet adhesion rates. It was confirmed that the platelet adhesion rate is not only a function of Reynolds number (or wall shear rate) but also the stenosis severity of the vessel. General correlations for adhesion rates of platelets as functions of stenosis and Reynolds number were obtained based on these cases. Finally using the new platelet adhesion rates, the model was applied to different experimental systems and shown to agree well with measured platelet deposition.

Multi-Constituent Simulation of Thrombus Deposition

In this paper, we present a spatio-temporal mathematical model for simulating the formation and growth of a thrombus. Blood is treated as a multi-constituent mixture comprised of a linear fluid phase and a thrombus (solid) phase. The transport and reactions of 10 chemical and biological species are incorporated using a system of coupled convection-reaction-diffusion (CRD) equations to represent three processes in thrombus formation: initiation, propagation and stabilization. Computational fluid dynamic (CFD) simulations using the libraries of OpenFOAM were performed for two illustrative benchmark problems: in vivo thrombus growth in an injured blood vessel and in vitro thrombus deposition in micro-channels (1.5 mm × 1.6 mm × 0.1 mm) with small crevices (125 μm × 75 μm and 125 μm × 137 μm). For both problems, the simulated thrombus deposition agreed very well with experimental observations, both spatially and temporally. Based on the success with these two benchmark problems, which have very different flow conditions and biological environments, we believe that the current model will provide useful insight into the genesis of thrombosis in blood-wetted devices, and provide a tool for the design of less thrombogenic devices.

A General Shear-Dependent Model for Thrombus Formation

Modeling the transport, activation, and adhesion of platelets is crucial in predicting thrombus formation and growth following a thrombotic event in normal or pathological conditions. We propose a shear-dependent platelet adhesive model based on the Morse potential that is calibrated by existing invivo and invitro experimental data and can be used over a wide range of flow shear rates (100<γ˙<28,000s-1). We introduce an Eulerian-Lagrangian model where hemodynamics is solved on a fixed Eulerian grid, while platelets are tracked using a Lagrangian framework. A force coupling method is introduced for bidirectional coupling of platelet motion with blood flow. Further, we couple the calibrated platelet aggregation model with a tissue-factor/contact pathway coagulation cascade, representing the relevant biology of thrombin generation and the subsequent fibrin deposition. The range of shear rates covered by the proposed model encompass venous and arterial thrombosis, ranging from low-shear-rate conditions in abdominal aortic aneurysms and thoracic aortic dissections to thrombosis in stenotic arteries following plaque rupture, where local shear rates are extremely high.

Hemostasis (thrombus formation) is the normal physiological response that prevents significant blood loss after vascular injury. The resulting clots can form under different flow conditions in the veins as well as the arteries. The excessive and undesirable formation of clots (i.e., thrombosis) in our circulatory system may lead to significant morbidity and mortality. Some of these pathologies are deep vein thrombosis and pulmonary embolism and atherothrombosis (thrombosis triggered by plaque rupture) in coronary arteries, to name a few. The process of clot formation and growth at a site on a blood vessel wall involves a number of simultaneous processes including: multiple chemical reactions in the coagulation cascade, species transport and platelet adhesion all of which are strongly influenced by the hydrodynamic forces. Numerical models for blood clotting normally focus on one of the processes under a specific flow condition. Here, we propose a general numerical model that encompass a wide range of hemodynamic conditions in the veins and arteries, with individual platelets and their adhesive dynamics included explicitly in the models. Further, we include the biochemistry of coagulation cascade, which is essential to modeling thrombus formation, and couple that to our platelet aggregation model. The simulation results—tested against three different experiments—demonstrate that the proposed model is effective in capturing the invivo and invitro experimental observations.

Choice of a hemodynamic model for occlusive thrombosis in arteries

Intravascular thrombosis can lead to heart attacks and strokes that together are the leading causes of death in the US (Kochanek, K.D., Murphy, S.L., Xu, J.Q., 2014). The ability to identify the offending biofluid mechanical conditions and predict the timescale of thrombotic occlusion in vessels and devices may improve patient outcomes. A computational model was developed to describe the growth of thrombus based on the local hemodynamic shear rate. The model predicts thrombus deposition based on initial geometric and fluid mechanical conditions, which are updated throughout the simulation to reflect the changing lumen dimensions. Thrombus growth and occlusion from whole blood was measured in in vitro experiments using stenotic glass capillary tubes, a PDMS microfluidic channel, and a PTFE stenotic aorto-iliac graft. Comparison of the predicted occlusion times to experimental results shows excellent agreement. The results indicate that local shear rate plays a critical role in acute thrombosis, and that hemodynamic characterization may have clinical utility.

High fidelity computational simulation of thrombus formation in Thoratec HeartMate II continuous flow ventricular assist device

Continuous flow ventricular assist devices (cfVADs) provide a life-saving therapy for severe heart failure. However, in recent years, the incidence of device-related thrombosis (resulting in stroke, device-exchange surgery or premature death) has been increasing dramatically, which has alarmed both the medical community and the FDA. The objective of this study was to gain improved understanding of the initiation and progression of thrombosis in one of the most commonly used cfVADs, the Thoratec HeartMate II. A computational fluid dynamics simulation (CFD) was performed using our recently updated mathematical model of thrombosis. The patterns of deposition predicted by simulation agreed well with clinical observations. Furthermore, thrombus accumulation was found to increase with decreased flow rate, and can be completely suppressed by the application of anticoagulants and/or improvement of surface chemistry. To our knowledge, this is the first simulation to explicitly model the processes of platelet deposition and thrombus growth in a continuous flow blood pump and thereby replicate patterns of deposition observed clinically. The use of this simulation tool over a range of hemodynamic, hematological, and anticoagulation conditions could assist physicians to personalize clinical management to mitigate the risk of thrombosis. It may also contribute to the design of future VADs that are less thrombogenic.