Search

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Over the past decade, much of the development of computational models of device-related thrombosis has focused on platelet activity. While those models have been successful in predicting thrombus formation in medical devices operating at high shear rates (> 5000 s-1), they cannot be directly applied to low-shear devices, such as blood oxygenators and catheters, where emerging information suggest that fibrin formation is the predominant mechanism of clotting and platelet activity plays a secondary role. In the current work, we augment an existing platelet-based model of thrombosis with a partial model of the coagulation cascade that includes contact activation of factor XII and fibrin production. To calibrate the model, we simulate a backward-facing-step flow channel that has been extensively characterized in-vitro. Next, we perform blood perfusion experiments through a microfluidic chamber mimicking a hollow fiber membrane oxygenator and validate the model against these observations. The simulation results closely match the time evolution of the thrombus height and length in the backward-facing-step experiment. Application of the model to the microfluidic hollow fiber bundle chamber capture both gross features such as the increasing clotting trend towards the outlet of the chamber, as well as finer local features such as the structure of fibrin around individual hollow fibers. Our results are in line with recent findings that suggest fibrin production, through contact activation of factor XII, drives the thrombus formation in medical devices operating at low shear rates with large surface area to volume ratios., Competing Interests: The authors have declared that no competing interests exist.