A computational model of human blood clottings: simulation, analysis, control, and validation

Complex biological systems pose many challenges to researchers, including, inter alia, choice of computational model, with its consequences for simulation and analysis; methods of manipulating the system exogenously (control); and model validation. I attempt to address these issues for human blood clotting. By treating the system as comprising interacting discrete and continuous aspects, i.e. as hybrid, the entire coagulation cascade may be simulated: Blood proteins, elemental ions, and other state elements are modeled either as real-valued concentrations or as binary variables (present/absent); interactions are rendered either as ODEs (per their chemical equations) or as discrete events. Techniques from nonlinear control theory are then used to devise drug therapies for diseased patients. Finally, the model is used to warp variations in the input parameters—rate constants and initial conditions—into an output space where pathologies and healthy clotting are cleanly separated by a semi-supervised clustering analysis. This serves to validate the model as well as to summarize efficiently the predicted clinical consequences of individual variations.

[1]  H C Hemker,et al.  The limits of simulation of the clotting system , 2006, Journal of thrombosis and haemostasis : JTH.

[2]  S. Olson,et al.  Transient kinetics of heparin-catalyzed protease inactivation by antithrombin III. Linkage of protease-inhibitor-heparin interactions in the reaction with thrombin. , 1988, The Journal of biological chemistry.

[3]  K. Mann,et al.  "Normal" thrombin generation. , 1999, Blood.

[4]  J. Francis,et al.  Mechanisms in Blood Coagulation, Fibrinolysis and the Complement System , 1992 .

[5]  Gary J. Balas,et al.  Hybrid Systems: Review and Recent Progress , 2003 .

[6]  M N Liebman,et al.  Qualitative modeling of normal blood coagulation and its pathological states using stochastic activity networks. , 1997, International journal of biological macromolecules.

[7]  H Matsuno,et al.  Hybrid Petri net representation of gene regulatory network. , 1999, Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing.

[8]  S. Sastry Nonlinear Systems: Analysis, Stability, and Control , 1999 .

[9]  K. Mann,et al.  What is all that thrombin for? , 2003, Journal of thrombosis and haemostasis : JTH.

[10]  K. Mann,et al.  A model describing the inactivation of factor Va by APC: bond cleavage, fragment dissociation, and product inhibition. , 1999, Biochemistry.

[11]  Sharene D Bungay,et al.  A mathematical model of lipid-mediated thrombin generation. , 2003, Mathematical medicine and biology : a journal of the IMA.

[12]  M A Khanin,et al.  Analysis of the activated partial thromboplastin time test using mathematical modeling. , 2001, Thrombosis research.

[13]  Jean-Luc Gouzé,et al.  Hybrid Modeling and Simulation of Genetic Regulatory Networks: A Qualitative Approach , 2003, HSCC.

[14]  S. Sastry,et al.  The Art of Hybrid Systems , 2005 .

[15]  P. Kokotovic,et al.  Nonlinear control via approximate input-output linearization: the ball and beam example , 1992 .

[16]  Stephen Neuendorffer,et al.  Modeling real-world control systems: beyond hybrid systems , 2004, Proceedings of the 2004 Winter Simulation Conference, 2004..

[17]  S. Everse,et al.  Predicting the pharmacology of thrombin inhibitors , 2003, Journal of thrombosis and haemostasis : JTH.

[18]  Mikhail A Panteleev,et al.  Tissue factor pathway inhibitor: a possible mechanism of action. , 2002, European journal of biochemistry.

[19]  Ádám M. Halász,et al.  Understanding the Bacterial Stringent Response Using Reachability Analysis of Hybrid Systems , 2004, HSCC.

[20]  Deyan Luan,et al.  Computationally Derived Points of Fragility of a Human Cascade Are Consistent with Current Therapeutic Strategies , 2007, PLoS Comput. Biol..

[21]  Pravin Varaiya,et al.  What's decidable about hybrid automata? , 1995, STOC '95.

[22]  Rainer Drath Description of Hybrid Systems by Modified Petri Nets , 2002 .

[23]  Georges Bastin,et al.  On state accessibility in reaction systems , 1993 .

[24]  R. Horowitz,et al.  Control design of an automated highway system , 2000, Proceedings of the IEEE.

[25]  W. Blair Series solution to the general linear time varying system , 1971 .

[26]  A. Schmaier,et al.  Activation of the plasma kallikrein/kinin system on endothelial cell membranes. , 1999, Immunopharmacology.

[27]  F. Rosner,et al.  Medicine in the Bible and the Talmud: Selections from Classical Jewish Sources , 1977 .

[28]  A. Fersht Structure and mechanism in protein science , 1998 .

[29]  M. Mosesson,et al.  Fibrinogen functions and fibrin assembly , 2000 .

[30]  Vijay Kumar,et al.  Hybrid control of formations of robots , 2001, Proceedings 2001 ICRA. IEEE International Conference on Robotics and Automation (Cat. No.01CH37164).

[31]  L.C.G.J.M. Habets,et al.  Control of multi-affine systems on rectangles with applications to hybrid biomolecular networks , 2002, Proceedings of the 41st IEEE Conference on Decision and Control, 2002..

[32]  Jacqueline Signorini,et al.  Object-Oriented Specification of Complex Bio-computing Processes: A Case Study of a Network of Proteolytic Enzymes , 2004, BioADIT.

[33]  Falko Bause,et al.  Stochastic Petri Nets: An Introduction to the Theory , 2012, PERV.

[34]  Martin Buss,et al.  Hybrid Control of Multi-fingered Dextrous Robotic Hands , 2002 .

[35]  Hans Gregersen,et al.  Dynamic model of the role of platelets in the blood coagulation system. , 2002, Medical engineering & physics.

[36]  S. Shankar Sastry,et al.  Modeling Subtilin Production in Bacillus subtilis Using Stochastic Hybrid Systems , 2004, HSCC.

[37]  Mikhail A Panteleev,et al.  Spatial propagation and localization of blood coagulation are regulated by intrinsic and protein C pathways, respectively. , 2006, Biophysical journal.

[38]  Catalina M. Lladó,et al.  PIPE v 2 . 5 : a Petri Net Tool for Performance Modeling , 2007 .

[39]  Atsushi Doi,et al.  Biopathways representation and simulation on hybrid functional Petri net , 2003, Silico Biol..

[40]  Natasha A. Neogi,et al.  Dynamic Partitioning of Large Discrete Event Biological Systems for Hybrid Simulation and Analysis , 2004, HSCC.

[41]  A. Mohan,et al.  Modeling the growth and dissolution of clots in flowing blood , 2006 .

[42]  M. Shuman,et al.  Cleavage of blood coagulation factor XIII and fibrinogen by thrombin during in vitro clotting. , 1985, The Journal of clinical investigation.

[43]  Ian A. Hiskens,et al.  Trajectory Sensitivity Analysis of Hybrid Systems , 2000 .

[44]  Ming Chen,et al.  Quantitative Petri net model of gene regulated metabolic networks in the cell , 2003, Silico Biol..

[45]  F Keller,et al.  The quick machine--a mathematical model for the extrinsic activation of coagulation. , 1994, Haemostasis.

[46]  Zia Shariat-Madar,et al.  Assembly and activation of the plasma kallikrein/kinin system: a new interpretation. , 2002, International immunopharmacology.

[47]  Atsushi Doi,et al.  Genomic Object Net: Object Oriented Representation of Biological Systems , 2000 .

[48]  Claire J. Tomlin,et al.  Lateral Inhibition through Delta-Notch Signaling: A Piecewise Affine Hybrid Model , 2001, HSCC.

[49]  John Lygeros,et al.  Lecture Notes on Hybrid Systems , 2004 .

[50]  M A Khanin,et al.  Mathematical model for the blood coagulation prothrombin time test. , 1998, Thrombosis research.

[51]  Kenneth G. Mann,et al.  PlA2 Polymorphism of &bgr;3 Integrins Is Associated With Enhanced Thrombin Generation and Impaired Antithrombotic Action of Aspirin at the Site of Microvascular Injury , 2001, Circulation.

[52]  Dale E. Seborg,et al.  Feedback linearizing control , 1997 .

[53]  Alexandre M. Bayen,et al.  Hybrid verification of an interface for an automatic landing , 2002, Proceedings of the 41st IEEE Conference on Decision and Control, 2002..

[54]  O. Stursberg,et al.  Continuous-discrete interactions in chemical processing plants , 2000, Proceedings of the IEEE.

[55]  D C Teller,et al.  Crosslinking kinetics of the human transglutaminase, factor XIII[A2], acting on fibrin gels and gamma-chain peptides. , 1997, Biochemistry.

[56]  X. H. Xu,et al.  The kinetic model and simulation of blood coagulation--the kinetic influence of activated protein C. , 2004, Medical engineering & physics.

[57]  S. Shankar Sastry,et al.  Switching through singularities , 1998 .

[58]  Olaf Czogalla,et al.  Modelling and Simulation of Controlled Road Traffic , 2002 .

[59]  P. Grant,et al.  The molecular physiology and pathology of fibrin structure/function. , 2005, Blood reviews.

[60]  S. S. Ravi,et al.  Clustering with Constraints: Feasibility Issues and the k-Means Algorithm , 2005, SDM.

[61]  V. Yee,et al.  Blood coagulation factor XIII: structure and function. , 1999, Thrombosis research.

[62]  Michael J. Kurtz,et al.  Input-output linearizing control of constrained nonlinear processes , 1997 .

[63]  S. Béguin,et al.  The Calibrated Automated Thrombogram (CAT): a universal routine test for hyper- and hypocoagulability , 2002, Pathophysiology of Haemostasis and Thrombosis.

[64]  Leo Breiman,et al.  Classification and Regression Trees , 1984 .

[65]  Alexandre M. Bayen,et al.  Validating a Hamilton-Jacobi Approximation to Hybrid System Reachable Sets , 2001, HSCC.

[66]  R J Leipold,et al.  Mathematical Model of Serine Protease Inhibition in the Tissue Factor Pathway to Thrombin (*) , 1995, The Journal of Biological Chemistry.

[67]  K. Mann,et al.  Thrombin functions during tissue factor-induced blood coagulation. , 2002, Blood.