1
|
Terekhov KM, Butakov ID, Danilov AA, Vassilevski YV. Dynamic adaptive moving mesh finite-volume method for the blood flow and coagulation modeling. Int J Numer Method Biomed Eng 2023; 39:e3731. [PMID: 38018385 DOI: 10.1002/cnm.3731] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Revised: 04/05/2023] [Accepted: 04/27/2023] [Indexed: 11/30/2023]
Abstract
In this work, we develop numerical methods for the solution of blood flow and coagulation on dynamic adaptive moving meshes. We consider the blood flow as a flow of incompressible Newtonian fluid governed by the Navier-Stokes equations. The blood coagulation is introduced through the additional Darcy term, with a permeability coefficient dependent on reactions. To this end, we introduce moving mesh collocated finite-volume methods for the Navier-Stokes equations, advection-diffusion equations, and a method for the stiff cascade of reactions. A monolithic nonlinear system is solved to advance the solution in time. The finite volume method for the Navier-Stokes equations features collocated arrangement of pressure and velocity unknowns and a coupled momentum and mass flux. The method is conservative and inf-sup stable despite the saddle point nature of the system. It is verified on a series of analytical problems and applied to the blood flow problem in the deforming domain of the right ventricle, reconstructed from a time series of computed tomography scans. At last, we demonstrate the ability to model the coagulation process in deforming microfluidic capillaries.
Collapse
Affiliation(s)
- Kirill M Terekhov
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russia
- Sirius University of Science and Technology, Sochi, Russia
| | - Ivan D Butakov
- Sirius University of Science and Technology, Sochi, Russia
| | - Alexander A Danilov
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russia
- Sirius University of Science and Technology, Sochi, Russia
- Sechenov University, Moscow, Russia
| | - Yuri V Vassilevski
- Marchuk Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow, Russia
- Sirius University of Science and Technology, Sochi, Russia
- Sechenov University, Moscow, Russia
| |
Collapse
|
2
|
Ratto N, Bouchnita A, Chelle P, Marion M, Panteleev M, Nechipurenko D, Tardy-Poncet B, Volpert V. Patient-Specific Modelling of Blood Coagulation. Bull Math Biol 2021; 83:50. [PMID: 33772645 PMCID: PMC7998098 DOI: 10.1007/s11538-021-00890-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 03/12/2021] [Indexed: 10/24/2022]
Abstract
Blood coagulation represents one of the most studied processes in biomedical modelling. However, clinical applications of this modelling remain limited because of the complexity of this process and because of large inter-patient variation of the concentrations of blood factors, kinetic constants and physiological conditions. Determination of some of these patients-specific parameters is experimentally possible, but it would be related to excessive time and material costs impossible in clinical practice. We propose in this work a methodological approach to patient-specific modelling of blood coagulation. It begins with conventional thrombin generation tests allowing the determination of parameters of a reduced kinetic model. Next, this model is used to study spatial distributions of blood factors and blood coagulation in flow, and to evaluate the results of medical treatment of blood coagulation disorders.
Collapse
Affiliation(s)
- N Ratto
- UMR 5208 CNRS, Institute Camille Jordan, Ecole Centrale de Lyon, Ecully, France
| | - A Bouchnita
- University of Texas at Austin, Austin, TX, 78712, USA
| | - P Chelle
- Center for Health Engineering, UMR 5307, Ecole Nationale Superieure des Mines de Saint-Etienne, 2023, Saint-Étienne, France.,EA3065, University Jean Monnet, 42023, Saint-Étienne, France
| | - M Marion
- UMR 5208 CNRS, Institute Camille Jordan, Ecole Centrale de Lyon, Ecully, France
| | - M Panteleev
- Faculty of Physics, Lomonosov Moscow State University, Moscow, Russia.,Center for Theoretical Problems of Physicochemical Pharmacology of the Russian Academy of Sciences, Moscow, Russia.,National Medical Research Center of Pediatric Hematology, Oncology and Immunology named after Dmitry Rogachev, Moscow, Russia
| | - D Nechipurenko
- Faculty of Physics, Lomonosov Moscow State University, Moscow, Russia.,Center for Theoretical Problems of Physicochemical Pharmacology of the Russian Academy of Sciences, Moscow, Russia.,National Medical Research Center of Pediatric Hematology, Oncology and Immunology named after Dmitry Rogachev, Moscow, Russia
| | - B Tardy-Poncet
- EA3065, University Jean Monnet, 42023, Saint-Étienne, France.,Inserm CIC1408, CHU de Saint-Etienne, 42023, Saint-Étienne, France
| | - V Volpert
- UMR 5208 CNRS, Institut Camille Jordan, University Lyon 1, 69622, Villeurbanne, France. .,INRIA Team Dracula, INRIA Lyon La Doua, 69603, Villeurbanne, France. .,Peoples' Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow, Russia, 117198.
| |
Collapse
|
3
|
Bouchnita A, Terekhov K, Nony P, Vassilevski Y, Volpert V. A mathematical model to quantify the effects of platelet count, shear rate, and injury size on the initiation of blood coagulation under venous flow conditions. PLoS One 2020; 15:e0235392. [PMID: 32726315 PMCID: PMC7390270 DOI: 10.1371/journal.pone.0235392] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 06/16/2020] [Indexed: 11/18/2022] Open
Abstract
Platelets upregulate the generation of thrombin and reinforce the fibrin clot which increases the incidence risk of venous thromboembolism (VTE). However, the role of platelets in the pathogenesis of venous cardiovascular diseases remains hard to quantify. An experimentally validated model of thrombin generation dynamics is formulated. The model predicts that a high platelet count increases the peak value of generated thrombin as well as the endogenous thrombin potential (ETP) as reported in experimental data. To investigate the effects of platelets density, shear rate, and wound size on the initiation of blood coagulation, we calibrate a previously developed model of venous thrombus formation and implement it in 3D using a novel cell-centered finite-volume solver. We conduct numerical simulations to reproduce in vitro experiments of blood coagulation in microfluidic capillaries. Then, we derive a reduced one-equation model of thrombin distribution from the previous model under simplifying hypotheses and we use it to determine the conditions of clotting initiation on the platelet count, the shear rate, and the plasma composition. The initiation of clotting also exhibits a threshold response to the size of the wounded region in good agreement with the reported experimental findings.
Collapse
Affiliation(s)
| | - Kirill Terekhov
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
| | - Patrice Nony
- Services de Pharmacologie Clinique, Hospices Civils de Lyon, Lyon, France
| | - Yuri Vassilevski
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
- Sechenov University, Moscow, Russia
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Vitaly Volpert
- Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences, Moscow, Russia
- Institut Camille Jordan, Université Lyon 1, Villeurbanne, France
- INRIA team Dracula, INRIA Lyon La Doua, Villeurbanne, France
- Peoples’ Friendship University of Russia (RUDN University), Moscow, Russia
| |
Collapse
|
4
|
Ratto N, Tokarev A, Chelle P, Tardy-Poncet B, Volpert V. Clustering of Thrombin Generation Test Data Using a Reduced Mathematical Model of Blood Coagulation. Acta Biotheor 2020; 68:21-43. [PMID: 31853681 DOI: 10.1007/s10441-019-09372-w] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Accepted: 12/04/2019] [Indexed: 11/28/2022]
Abstract
Correct interpretation of the data from integral laboratory tests, including Thrombin Generation Test (TGT), requires biochemistry-based mathematical models of blood coagulation. The purpose of this study is to describe the experimental TGT data from healthy donors and hemophilia A (HA) and B (HB) patients. We derive a simplified ODE model and apply it to analyze the TGT data from healthy donors and HA/HB patients with in vitro added tissue factor pathway inhibitor (TFPI) antibody. This model allows the characterization of hemophilia patients in the space of three most important model parameters. The proposed approach may provide a new quantitative tool for the analysis of experimental TGT. Also, it gives a reduced model of coagulation verified against clinical data to be used in future theoretical large-scale modeling of thrombosis in flow.
Collapse
Affiliation(s)
- N Ratto
- Institute Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France.
| | - A Tokarev
- Peoples' Friendship University of Russia (RUDN University), Moscow, 117198, Russian Federation
| | - P Chelle
- Center for Health Engineering, UMR 5307, Ecole Nationale Supérieure des Mines de Saint-Etienne, 42023, Saint-Etienne, France
- Unite Inserm Sainbiose U1059, Université Jean Monnet, 42023, Saint-Etienne, France
| | - B Tardy-Poncet
- Unite Inserm Sainbiose U1059, Université Jean Monnet, 42023, Saint-Etienne, France
- Inserm CIC1408, CHU de Saint-Etienne, 42023, Saint-Etienne, France
- CRC Hémophilie CHU St Etienne, 42023, Saint-Etienne, France
| | - V Volpert
- Institute Camille Jordan, UMR 5208 CNRS, University Lyon 1, 69622, Villeurbanne, France
- Institute Camille Jordan, INRIA, Université de Lyon, Université Lyon 1, 69200, Villeurbanne, France
- Peoples' Friendship University of Russia (RUDN University), Moscow, 117198, Russian Federation
| |
Collapse
|
5
|
Kushchenko YK, Belyaev AV. Effects of hydrophobicity, tethering and size on flow-induced activation of von Willebrand factor multimers. J Theor Biol 2019; 485:110050. [PMID: 31618612 DOI: 10.1016/j.jtbi.2019.110050] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2019] [Revised: 09/12/2019] [Accepted: 10/12/2019] [Indexed: 01/14/2023]
Abstract
Von Willebrand factor (VWF) is a multimeric protein of blood plasma that mediates platelet adhesion to injury under strong hemodynamic flows in arterias and microvasvulature. We present a 3D coarse-grained computer model of VWF multimers in flowing viscous fluid that explicitely grasps the dynamics, the conformational changes and the hydrodynamics-induced activation of adhesivity of these protein concatamers to GPIb receptor of blood platelets. The model is based on the fluctuating Lattice Boltzmann method for modelling the hydrodynamics in the simulation box and the Lagrangian particle dynamics coupled to the fluid by a viscous drag force. The model has been validated by the comparison with the experimental data found in literature. We studied the effect of hydrophobic interactions on the conformational dynamics of VWF multimers. The simulations suggest that the contour length is an important parameter that controls the functionality of VWF multimers in blood. We also demonstrate that tethering to the surface of a vessel wall promoted the flow-induced activation of VWF, while those multimers that remain untethered and move freely in the blood plasma require a stronger shearing to get activated.
Collapse
Affiliation(s)
- Yulia K Kushchenko
- Lomonosov Moscow State University, Faculty of Physics, Moscow 119991, Russia
| | - Aleksey V Belyaev
- Lomonosov Moscow State University, Faculty of Physics, Moscow 119991, Russia; S.M. Nikol'skii Mathematical Institute, RUDN University, Moscow 115419, Russia.
| |
Collapse
|
6
|
Chen J, Diamond SL. Reduced model to predict thrombin and fibrin during thrombosis on collagen/tissue factor under venous flow: Roles of γ'-fibrin and factor XIa. PLoS Comput Biol 2019; 15:e1007266. [PMID: 31381558 PMCID: PMC6695209 DOI: 10.1371/journal.pcbi.1007266] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 08/15/2019] [Accepted: 07/08/2019] [Indexed: 01/29/2023] Open
Abstract
During thrombosis, thrombin generates fibrin, however fibrin reversibly binds thrombin with low affinity E-domain sites (KD = 2.8 μM) and high affinity γ’-fibrin sites (KD = 0.1 μM). For blood clotting on collagen/tissue factor (1 TF-molecule/μm2) at 200 s-1 wall shear rate, high μM-levels of intraclot thrombin suggest robust prothrombin penetration into clots. Setting intraclot zymogen concentrations to plasma levels (and neglecting cofactor rate limitations) allowed the linearization of 7 Michaelis-Menton reactions between 6 species to simulate intraclot generation of: Factors FXa (via TF/VIIa or FIXa), FIXa (via TF/FVIIa or FXIa), thrombin, fibrin, and FXIa. This reduced model [7 rates, 2 KD’s, enzyme half-lives~1 min] predicted the measured clot elution rate of thrombin-antithrombin (TAT) and fragment F1.2 in the presence and absence of the fibrin inhibitor Gly-Pro-Arg-Pro. To predict intraclot fibrin reaching 30 mg/mL by 15 min, the model required fibrinogen penetration into the clot to be strongly diffusion-limited (actual rate/ideal rate = 0.05). The model required free thrombin in the clot (~100 nM) to have an elution half-life of ~2 sec, consistent with measured albumin elution, with most thrombin (>99%) being fibrin-bound. Thrombin-feedback activation of FXIa became prominent and reached 5 pM FXIa at >500 sec in the simulation, consistent with anti-FXIa experiments. In predicting intrathrombus thrombin and fibrin during 15-min microfluidic experiments, the model revealed “cascade amplification” from 30 pM levels of intrinsic tenase to 15 nM prothrombinase to 15 μM thrombin to 90 μM fibrin. Especially useful for multiscale simulation, this reduced model predicts thrombin and fibrin co-regulation during thrombosis under flow. During blood clotting events, a complex series of reaction are involved. Simulation gives insights to the concentration of different enzymes which are at too low of concentration to be detected. However, the models are often large and difficult to solve for clotting under flow conditions. With a thin film approximation, we were able to simplify clotting under flow with parameters from literature, with only 3 adjusted in order to fit the experimental data. This model gave insights into the dynamics of the species involved, and the roles of γ’-fibrin and thrombin feedback activation. This reduced model may be useful in further multiscale simulations.
Collapse
Affiliation(s)
- Jason Chen
- Department of Chemical and Biomolecular Engineering, Institute for Medicine and Engineering, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Scott L. Diamond
- Department of Chemical and Biomolecular Engineering, Institute for Medicine and Engineering, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
- * E-mail:
| |
Collapse
|
7
|
Belyaev AV, Dunster JL, Gibbins JM, Panteleev MA, Volpert V. Modeling thrombosis in silico: Frontiers, challenges, unresolved problems and milestones. Phys Life Rev 2018; 26-27:57-95. [PMID: 29550179 DOI: 10.1016/j.plrev.2018.02.005] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Revised: 02/21/2018] [Accepted: 02/24/2018] [Indexed: 12/24/2022]
Abstract
Hemostasis is a complex physiological mechanism that functions to maintain vascular integrity under any conditions. Its primary components are blood platelets and a coagulation network that interact to form the hemostatic plug, a combination of cell aggregate and gelatinous fibrin clot that stops bleeding upon vascular injury. Disorders of hemostasis result in bleeding or thrombosis, and are the major immediate cause of mortality and morbidity in the world. Regulation of hemostasis and thrombosis is immensely complex, as it depends on blood cell adhesion and mechanics, hydrodynamics and mass transport of various species, huge signal transduction networks in platelets, as well as spatiotemporal regulation of the blood coagulation network. Mathematical and computational modeling has been increasingly used to gain insight into this complexity over the last 30 years, but the limitations of the existing models remain profound. Here we review state-of-the-art-methods for computational modeling of thrombosis with the specific focus on the analysis of unresolved challenges. They include: a) fundamental issues related to physics of platelet aggregates and fibrin gels; b) computational challenges and limitations for solution of the models that combine cell adhesion, hydrodynamics and chemistry; c) biological mysteries and unknown parameters of processes; d) biophysical complexities of the spatiotemporal networks' regulation. Both relatively classical approaches and innovative computational techniques for their solution are considered; the subjects discussed with relation to thrombosis modeling include coarse-graining, continuum versus particle-based modeling, multiscale models, hybrid models, parameter estimation and others. Fundamental understanding gained from theoretical models are highlighted and a description of future prospects in the field and the nearest possible aims are given.
Collapse
|