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.

[1]  J. Antaki,et al.  Study of blood flow in several benchmark micro-channels using a two-fluid approach. , 2015, International journal of engineering science.

[2]  M. Massoudi,et al.  On the coefficients of the interaction forces in a two-phase flow of a fluid infused with particles , 2014 .

[3]  D. Monroe,et al.  Platelets and thrombin generation. , 2002, Arteriosclerosis, thrombosis, and vascular biology.

[4]  B. Savage,et al.  Mechanisms of platelet aggregation. , 2001, Current opinion in hematology.

[5]  Hai-Chao Han,et al.  Effect of Red Blood Cells on Platelet Activation and Thrombus Formation in Tortuous Arterioles , 2013, Front. Bioeng. Biotechnol..

[6]  David Gladstone Review article , 2005, Health Care Analysis.

[7]  Donghui Zhu,et al.  Mathematical modeling of blood coagulation cascade: kinetics of intrinsic and extrinsic pathways in normal and deficient conditions , 2007, Blood coagulation & fibrinolysis : an international journal in haemostasis and thrombosis.

[8]  M. Frojmovic,et al.  Dynamic measurements of the platelet membrane glycoprotein IIb-IIIa receptor for fibrinogen by flow cytometry. II. Platelet size-dependent subpopulations. , 1991, Biophysical journal.

[9]  B. Samuelsson,et al.  Effect of Leucotrienes C4, D4, Prostacyclin and Thromboxane A2 on isolated human cerebral arteries , 2005, Acta Neurochirurgica.

[10]  Drake D Pedersen,et al.  Visualization and analysis of biomaterial-centered thrombus formation within a defined crevice under flow. , 2016, Biomaterials.

[11]  M. Simon,et al.  Activation of G12/G13 Results in Shape Change and Rho/Rho-Kinase–mediated Myosin Light Chain Phosphorylation in Mouse Platelets , 1999, The Journal of cell biology.

[12]  A. Fogelson,et al.  Computational model of whole blood exhibiting lateral platelet motion induced by red blood cells , 2010, International journal for numerical methods in biomedical engineering.

[13]  Zhiliang Xu,et al.  A multiscale model of thrombus development , 2008, Journal of The Royal Society Interface.

[14]  Nadine Aubry,et al.  A numerical study of blood flow using mixture theory. , 2014, International journal of engineering science.

[15]  S. Jackson,et al.  Signaling events underlying thrombus formation , 2003, Journal of thrombosis and haemostasis : JTH.

[16]  James F. Antaki,et al.  Computational Simulation of Platelet Deposition and Activation: I. Model Development and Properties , 1999, Annals of Biomedical Engineering.

[17]  C. O’Leary,et al.  Bleeding and Thrombosis , 2021, Oncologic Emergency Medicine.

[18]  S. Coughlin,et al.  Protease-activated receptors 1 and 4 mediate activation of human platelets by thrombin. , 1999, The Journal of clinical investigation.

[19]  G. V. R. Born,et al.  Activation time of blood platelets , 1980, The Journal of Membrane Biology.

[20]  Nader Moazami,et al.  An analysis of pump thrombus events in patients in the HeartWare ADVANCE bridge to transplant and continued access protocol trial. , 2014, The Journal of heart and lung transplantation : the official publication of the International Society for Heart Transplantation.

[21]  D. Ku,et al.  A low-volume, single pass in-vitro system of high shear thrombosis in a stenosis. , 2013, Thrombosis research.

[22]  S. Diamond,et al.  Pathological von Willebrand factor fibers resist tissue plasminogen activator and ADAMTS13 while promoting the contact pathway and shear‐induced platelet activation , 2015, Journal of thrombosis and haemostasis : JTH.

[23]  H. Engelberg,et al.  Plasma Heparin Levels in Normal Man , 1961, Circulation.

[24]  James F. Antaki,et al.  Computational Simulation of Platelet Deposition and Activation: II. Results for Poiseuille Flow over Collagen , 1999, Annals of Biomedical Engineering.

[25]  H. Weiss Platelets: Pathophysiology and antiplatelet drug therapy , 1982 .

[26]  Kumbakonam R. Rajagopal,et al.  A Model Incorporating Some of the Mechanical and Biochemical Factors Underlying Clot Formation and Dissolution in Flowing Blood , 2003 .

[27]  Scott L. Diamond,et al.  Thrombus Growth and Embolism on Tissue Factor-Bearing Collagen Surfaces Under Flow: Role of Thrombin With and Without Fibrin , 2012, Arteriosclerosis, thrombosis, and vascular biology.

[28]  James F. Antaki,et al.  Micro-Flow Visualization of Red Blood Cell-Enhanced Platelet Concentration at Sudden Expansion , 2008, Annals of Biomedical Engineering.

[29]  P. Richardson Effect of Blood Flow Velocity on Growth Rate of Platelet Thrombi , 1973, Nature.

[30]  M. Acker,et al.  Pump thrombosis in the Thoratec HeartMate II device: An update analysis of the INTERMACS Registry. , 2015, The Journal of heart and lung transplantation : the official publication of the International Society for Heart Transplantation.

[31]  Kenneth A. Solen,et al.  Computational Model of Device-Induced Thrombosis and Thromboembolism , 2005, Annals of Biomedical Engineering.

[32]  S. Roberts,et al.  Chemistry, biochemistry, and pharmacological activity of prostanoids , 1979 .

[33]  James B. Young,et al.  What is the truth behind pump thrombosis in the HeartMate II device? A National Heart, Lung, and Blood Institute perspective based on data from the Interagency Registry for Mechanically Assisted Circulatory Support. , 2015, The Journal of heart and lung transplantation : the official publication of the International Society for Heart Transplantation.

[34]  Scott L Diamond,et al.  Platelet dynamics in three-dimensional simulation of whole blood. , 2014, Biophysical journal.

[35]  S. Diamond,et al.  A systems approach to hemostasis: 3. Thrombus consolidation regulates intrathrombus solute transport and local thrombin activity. , 2014, Blood.

[36]  Aaron L. Fogelson,et al.  Continuum models of platelet aggregation: formulation and mechanical properties , 1992 .

[37]  Manash S. Chatterjee,et al.  A molecular signaling model of platelet phosphoinositide and calcium regulation during homeostasis and P2Y1 activation. , 2008, Blood.

[38]  Zhiliang Xu,et al.  A multiscale model of venous thrombus formation with surface-mediated control of blood coagulation cascade. , 2010, Biophysical journal.

[39]  K. Zlobina,et al.  Platelet activation risk index as a prognostic thrombosis indicator , 2016, Scientific Reports.

[40]  A. Fogelson,et al.  Surface-mediated control of blood coagulation: the role of binding site densities and platelet deposition. , 2001, Biophysical journal.

[41]  Arie Blitz,et al.  Pump thrombosis-A riddle wrapped in a mystery inside an enigma. , 2014, Annals of cardiothoracic surgery.

[42]  N. Abraham,et al.  Additive effect of erythropoietin and heme on murine hematopoietic recovery after azidothymidine treatment. , 1993, Blood.

[43]  G. Born,et al.  Growth Rate in vivo of Platelet Thrombi, produced by Iontophoresis of ADP, as a Function of Mean Blood Flow Velocity , 1970, Nature.

[44]  M. Frojmovic,et al.  Dynamics of platelet glycoprotein IIb-IIIa receptor expression and fibrinogen binding. I. Quantal activation of platelet subpopulations varies with adenosine diphosphate concentration. , 1994, Biophysical journal.

[45]  J. Hirsh,et al.  Medical device‐induced thrombosis: what causes it and how can we prevent it? , 2015, Journal of thrombosis and haemostasis : JTH.

[46]  Mehrdad Massoudi,et al.  Flow of a fluid—solid mixture between flat plates , 1991 .

[47]  Scott L Diamond,et al.  A membrane-based microfluidic device for controlling the flux of platelet agonists into flowing blood. , 2008, Lab on a chip.

[48]  S. Offermanns The Role of Heterotrimeric G Proteins in Platelet Activation , 2000, Biological chemistry.

[49]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[50]  V. Fuster,et al.  Synergistic action of severe wall injury and shear forces on thrombus formation in arterial stenosis: definition of a thrombotic shear rate threshold. , 1994, Journal of the American College of Cardiology.

[51]  R. D'Agostino,et al.  Platelet Glycoprotein IIIa PlA Polymorphism, Fibrinogen, and Platelet Aggregability: The Framingham Heart Study , 2001, Circulation.

[52]  D. Ku,et al.  Fluid mechanics of vascular systems, diseases, and thrombosis. , 1999, Annual review of biomedical engineering.

[53]  G. Remuzzi,et al.  Reduced platelet thromboxane formation in uremia. Evidence for a functional cyclooxygenase defect. , 1983, The Journal of clinical investigation.

[54]  Scott L Diamond,et al.  Multiscale prediction of patient-specific platelet function under flow. , 2012, Blood.

[55]  M. Jamiolkowski Characterization of real time platelet deposition onto opaque surfaces under clinically-relevant flow conditions , 2015 .

[56]  B. Hodson,et al.  The effect of passage in vitro and in vivo on the properties of murine fibrosarcomas. II. Sensitivity to cell-mediated cytotoxicity in vitro. , 1985, British Journal of Cancer.

[57]  D. Ku,et al.  Platelet transport rates and binding kinetics at high shear over a thrombus. , 2013, Biophysical journal.

[58]  Kai Sun,et al.  Inhibitory effect of caffeic acid on ADP-induced thrombus formation and platelet activation involves mitogen-activated protein kinases , 2015, Scientific Reports.

[59]  H. Weiss,et al.  Red blood cells: their dual role in thrombus formation. , 1980, Science.

[60]  K R Rajagopal,et al.  A model for the formation, growth, and lysis of clots in quiescent plasma. A comparison between the effects of antithrombin III deficiency and protein C deficiency. , 2008, Journal of theoretical biology.

[61]  B. Furie,et al.  Mechanisms of thrombus formation. , 2008, The New England journal of medicine.

[62]  J. Weisel,et al.  Circulating Microparticles Alter Formation, Structure, and Properties of Fibrin Clots , 2015, Scientific Reports.

[63]  J.,et al.  Whitaker Lecture : Biorheology in Thrombosis Research , 2022 .

[64]  William R Wagner,et al.  Real time visualization and characterization of platelet deposition under flow onto clinically relevant opaque surfaces. , 2015, Journal of biomedical materials research. Part A.

[65]  L. Corash,et al.  Regulation of platelet heterogeneity: effects of thrombocytopenia on platelet volume and density. , 1990, Experimental hematology.

[66]  Edward L. Snyder,et al.  Platelets: Pathophysiology and Antiplatelet Drug Therapy , 1983, The Yale Journal of Biology and Medicine.

[67]  Aaron L Fogelson,et al.  Platelet-wall interactions in continuum models of platelet thrombosis: formulation and numerical solution. , 2004, Mathematical medicine and biology : a journal of the IMA.

[68]  S. Diamond,et al.  Thrombin flux and wall shear rate regulate fibrin fiber deposition state during polymerization under flow. , 2010, Biophysical journal.

[69]  Aaron L Fogelson,et al.  Blood clot formation under flow: the importance of factor XI depends strongly on platelet count. , 2012, Biophysical journal.

[70]  K. Rajagopal,et al.  A Model for the Formation and Lysis of Blood Clots , 2006, Pathophysiology of Haemostasis and Thrombosis.

[71]  S. Diamond,et al.  Simulation of Intrathrombus Fluid and Solute Transport Using In Vivo Clot Structures with Single Platelet Resolution , 2013, Annals of Biomedical Engineering.

[72]  S. Watson,et al.  A diacylglycerol kinase inhibitor, R59022, potentiates secretion by and aggregation of thrombin-stimulated human platelets. , 1987, The Biochemical journal.

[73]  T. Mizobe,et al.  In Vivo and In Vitro Studies of the Inhibitory Effect of Propofol on Human Platelet Aggregation , 1998, Anesthesiology.

[74]  Mikhail A. Panteleev,et al.  Mathematical Modeling and Computer Simulation in Blood Coagulation , 2006, Pathophysiology of Haemostasis and Thrombosis.

[75]  W. van Oeveren,et al.  Stirred, shaken, or stagnant: What goes on at the blood-biomaterial interface. , 2017, Blood reviews.

[76]  S. Jackson,et al.  Chapter 20 – Platelet Thrombus Formation in Flowing Blood , 2013 .

[77]  David N. Ku,et al.  A Predictive Model of High Shear Thrombus Growth , 2016, Annals of Biomedical Engineering.

[78]  V. Fuster,et al.  Thrombosis: Studies under Flow Conditions , 1987, Annals of the New York Academy of Sciences.

[79]  Ming Liang,et al.  Pentamethylquercetin (PMQ) reduces thrombus formation by inhibiting platelet function , 2015, Scientific Reports.

[80]  Steven Deutsch,et al.  Development of a computational model for macroscopic predictions of device-induced thrombosis , 2016, Biomechanics and modeling in mechanobiology.

[81]  J. D. Hellums,et al.  1993 Whitaker lecture: Biorheology in thrombosis research , 1994, Annals of Biomedical Engineering.

[82]  J. Maclouf,et al.  Antiphospholipid Antibodies Enhance Thrombin-Induced Platelet Activation and Thromboxane Formation , 1993, Thrombosis and Haemostasis.

[83]  Beverly H. Koller,et al.  Decreased platelet aggregation, increased bleeding time and resistance to thromboembolism in P2Y1-deficient mice , 1999, Nature Medicine.

[84]  Amy Henderson Squilacote The Paraview Guide , 2008 .

[85]  Aaron L Fogelson,et al.  Platelet motion near a vessel wall or thrombus surface in two-dimensional whole blood simulations. , 2013, Biophysical journal.