A two-dimensional multiscale model is introduced for studying formation of a thrombus (clot) in a blood vessel. It involves components for modelling viscous, incompressible blood plasma; non-activated and activated platelets; blood cells; activating chemicals; fibrinogen; and vessel walls and their interactions. The macroscale dynamics of the blood flow is described by the continuum Navier–Stokes equations. The microscale interactions between the activated platelets, the platelets and fibrinogen and the platelets and vessel wall are described through an extended stochastic discrete cellular Potts model. The model is tested for robustness with respect to fluctuations of basic parameters. Simulation results demonstrate the development of an inhomogeneous internal structure of the thrombus, which is confirmed by the preliminary experimental data. We also make predictions about different stages in thrombus development, which can be tested experimentally and suggest specific experiments. Lastly, we demonstrate that the dependence of the thrombus size on the blood flow rate in simulations is close to the one observed experimentally.

References

Alber M, Jiang Y& Kiskowski MLattice gas cellular automaton model for rippling and aggregation in myxobacteria. Physica D. 190, 2004a343–358.doi:10.1016/j.physd.2003.11.012. . Crossref, ISI, Google Scholar
Alber M, Kiskowski M& Jiang YTwo-stage aggregate formation via streams in myxobacteria. Phys. Rev. Lett. 93, 2004b068 301doi:10.1103/PhysRevLett.93.068102. . Crossref, ISI, Google Scholar
Alber, M., Kiskowski, M. Jiang, Y. & Newman, S. 2004c Biological lattice gas models. In Dynamics and bifurcation of patterns in dissipative systems, vol. 12 (eds G. Dangelmayr & I. Opera). World Scientific Series on Nonlinear Science, pp. 274–291. Singapore: World Scientific. Google Scholar
Alber M, Chen N, Glimm T& Lushnikov P.M. 2006Multiscale dynamics of biological cells with chemotactic interactions: from a discrete stochastic model to a continuous description. Phys. Rev. E. 73, 051 901doi:10.1103/PhysRevE.73.051901. . Crossref, ISI, Google Scholar
Alber M, Chen N, Glimm T& Lushnikov P.MTwo-dimensional multiscale model of cell motion in a chemotactic field. Single-cell-based models in biology and medicine, Anderson A.R.A, Chaplain M.A.J& Rejniak K.A. 2007app. 53–72. Eds. Basel, Switzerland:Birkhaauser. Crossref, Google Scholar
Alber M, Chen N, Newman S& Lushnikov P.MContinuous macroscopic limit of a discrete stochastic model for interaction of living cells. Phys. Rev. Lett. 99, 2007b168 102doi:10.1103/PhysRevLett.99.168102. . Crossref, ISI, Google Scholar
Anand M, Rajagopal K& Rajagopal K.R. 2005A model for the formation and lysis of blood clots. Pathophysiol. Haemost. Thromb. 34, 109–120.doi:10.1159/000089931. . Crossref, PubMed, Google Scholar
Badimon L, Badimon J.J, Galvez A& Chesebro J.H. 1986Influence of arterial damage and wall shear rate on platelet deposition. Ex vivo study in a swine model. Arteriosclerosis. 6, 312. Crossref, PubMed, Google Scholar
Baldwin S.A& Basmadjian D. 1994A mathematical model of thrombin production in blood coagulation, part I: the sparsely covered membrane case. Ann. Biomed. Eng. 22, 357–370.doi:10.1007/BF02368242. . Crossref, PubMed, ISI, Google Scholar
Barynin I.A, Starkov I.A& Khanin M.A. 1999Mathematical models in hemostasis physiology. Izv. Akad. Nauk. Ser. Biol. 1, 59–66. Google Scholar
Basmadjian D. 1990The effect of flow and mass transport in thrombogenesis. Ann. Biomed. Eng. 18, 685–709.doi:10.1007/BF02368455. . Crossref, PubMed, ISI, Google Scholar
Begent N& Born G.V. 1970Growth rate in vivo of platelet thrombi, produced by iontophoresis of ADP, as a function of mean blood flow velocity. Nature. 227, 926–930.doi:10.1038/227926a0. . Crossref, PubMed, ISI, Google Scholar
Bell, J. B., Colella, P. & Howell, H. 1991 An efficient second-order projection method for viscous incompressible flow. In Proc. 10th AIAA Computational Fluid Dynamics Conference, Honolulu, Hawaii, 24–27 June 1991. Google Scholar
Beltrami E& Jesty J. 1995Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation. Proc. Natl Acad. Sci. USA. 92, 8744–8748.doi:10.1073/pnas.92.19.8744. . Crossref, PubMed, ISI, Google Scholar
Beltrami E& Jesty J. 2001The role of membrane patch size and flow in regulating a proteolytic feedback threshold on a membrane: possible application in blood coagulation. Math. Biosci. 172, 1–13.doi:10.1016/S0025-5564(01)00064-5. . Crossref, PubMed, ISI, Google Scholar
Benoit M, Gabriel D, Gerisch G& Gaub H.E. 2000Discrete interactions in cell adhesion measured by single-molecule force spectroscopy. Nat. Cell Biol. 6, 313–317.doi:10.1038/35014000. . Crossref, ISI, Google Scholar
Casal A, Sumen C, Reddy T.E, Alber M.S& Lee P.P. 2005Agent-based modeling of the context dependency in t cell recognition. J. Theor. Biol. 236, 376–391.doi:10.1016/j.jtbi.2005.03.019. . Crossref, PubMed, ISI, Google Scholar
Chaturvedi R, et al.2005On multiscale approaches to three-dimensional modelling of morphogenesis. J. R. Soc. Interface. 2, 237–253.doi:10.1098/rsif.2005.0033. . Link, ISI, Google Scholar
Chen N, Glazier J.A, Izaguirre J& Alber M.S. 2007A parallel implementation of the cellular potts model for simulation of cell-based morphogenesis. Comp. Phys. Commun. 176, 670–681.doi:10.1016/j.cpc.2007.03.007. . Crossref, PubMed, ISI, Google Scholar
Chien S, Usami S, Taylor H.M, Lundberg J.L& Gregersen M.I. 1966Effect of hematocrit and plasma proteins on human blood rheology at low shear rates. J. Appl. Physiol. 21, 81–87. Crossref, PubMed, ISI, Google Scholar
Chien S, Sung K.L.P, Skalak R, Usami S& Tozeren A.L. 1978Theoretical and experimental studies on viscoelastic properties of erythrocyte membrane. Biophys. J. 24, 463–487. Crossref, PubMed, ISI, Google Scholar
Cickovski T.M, Huang C, Chaturvedi R, Glimm T, Hentschel H.G.E, Alber M.S, Glazier J.A, Newman S.A& Izaguirre J.A. 2005A framework for three-dimensional simulation of morphogenesis. IEEE/ACM Trans. Comput. Biol. Bioinform. 2, 273–288.doi:10.1109/TCBB.2005.46. . Crossref, PubMed, ISI, Google Scholar
Cickovski T.M, et al.2007From genes to organisms via the cell: a problem solving environment for multicellular development. Comput. Sci. Eng. 9, 50–60.doi:10.1109/MCSE.2007.74. . Crossref, PubMed, ISI, Google Scholar
Dan D, Mueller C, Chen K& Glazier J.A. 2005Solving the advection–diffusion equations in biological contexts using the cellular Potts model. Phys. Rev. E. 72, 41 909doi:10.1103/PhysRevE.72.041909. . Crossref, ISI, Google Scholar
David N.K. 1997Blood flow in arteries. Annu. Rev. Fluid Mech. 29, 399–434.doi:10.1146/annurev.fluid.29.1.399. . Crossref, ISI, Google Scholar
de Simone G, Devereux R.B, Chinali M, Best L.G, Lee E.T& Welty T.K. 2005Strong heart study. Hypertension. 45, 625–630.doi:10.1161/01.HYP.0000157526.07977.ec. . Crossref, PubMed, ISI, Google Scholar
Drasdo D& Hoehme S. 2005A single-cell based model to tumor growth in-vitro: monolayers and spheroids. Phys. Biol. 2, 133–147.doi:10.1088/1478-3975/2/3/001. . Crossref, PubMed, ISI, Google Scholar
Drasdo D, Kree R& McCaskill J. 1995Monte–Carlo approach to tissue–cell populations. Phys. Rev. E. 52, 6635doi:10.1103/PhysRevE.52.6635. . Crossref, ISI, Google Scholar
Fogelson A.L. 1992Continuum models of platelet aggregation: formulation and mechanical properties. SIAM J. Appl. Math. 52, 1089–1110.doi:10.1137/0152064. . Crossref, ISI, Google Scholar
Fogelson A.L. 1993Continuum models of platelet aggregation: mechanical properties chemically-induced phase transitions. Contemp. Math. 141, 279–294. Crossref, Google Scholar
Fogelson A.L& Guy R.D. 2004Platelet–wall interactions in continuum models of platelet thrombosis: formulation and numerical solution. Math. Med. Biol. 21, 293–334.doi:10.1093/imammb/21.4.293. . Crossref, PubMed, ISI, Google Scholar
Glazier J.A& Graner F.M.C. 1993Simulation of the differential adhesion driven rearrangement of biological cells. Phys. Rev. E. 47, 2128–2154.doi:10.1103/PhysRevE.47.2128. . Crossref, PubMed, ISI, Google Scholar
Glazier J.A, Balter A& Poplawski JMagnetization to morphogenesis: a brief history of the Glazier–Graner–Hogeweg model. Single-cell-based models in biology and medicine, Anderson A.R.A, Chaplain M.A.J& Rejniak K.A. 2007pp. 79–108. Eds. Basel, Switzerland:Birkhaauser. Crossref, Google Scholar
Graner F.M.C& Glazier J.A. 1992Simulation of biological cell sorting using a two-dimensional extended potts model. Phys. Rev. Lett. 69, 2013–2016.doi:10.1103/PhysRevLett.69.2013. . Crossref, PubMed, ISI, Google Scholar
Hirt C& Nichols B. 1981Volume of fluid (VOF) method for the dynamics of free boundaries. J. Comp. Phys. 39, 201–225.doi:10.1016/0021-9991(81)90145-5. . Crossref, ISI, Google Scholar
Jesty J, Beltrami E& Willems G. 1993Mathematical analysis of a proteolytic positive-feedback loop: dependence of lag time and enzyme yields on the initial conditions and kinetic parameters. Biochemistry. 32, 6266–6274.doi:10.1021/bi00075a021. . Crossref, PubMed, ISI, Google Scholar
Jiang Y, Pjesivac-Grbovic J, Cantrell C& Freyer J.P. 2005A multiscale model for avascular tumor growth. Biophys. J. 89, 3884–3894.doi:10.1529/biophysj.105.060640. . Crossref, PubMed, ISI, Google Scholar
Jones K.C& Mann K.G. 1994A model for the tissue factor pathway to thrombin. ii. A mathematical simulation. J. Biol. Chem. 269, 23 367–23 373. Crossref, ISI, Google Scholar
Khanin M.A& Semenov V.V. 1989A mathematical model of the kinetics of blood coagulation. J. Theor. Biol. 136, 127–134.doi:10.1016/S0022-5193(89)80220-6. . Crossref, PubMed, ISI, Google Scholar
Kiskowski M.A, Jiang Y& Alber M.S. 2004Role of streams in myxobacteria aggregate formation. Phys. Biol. 1, 173–183.doi:10.1088/1478-3967/1/3/005. . Crossref, PubMed, ISI, Google Scholar
Knewitz M.A& Mombach J.C.M. 2006Computer simulation of the influence of cellular adhesion on the morphology of the interface between tissues of proliferating and quiescent cells. Comput. Biol. Med. 36, 59–69.doi:10.1016/j.compbiomed.2004.08.002. . Crossref, PubMed, ISI, Google Scholar
Kothe D.B, Rider W.J, Mosso S.J, Brock J.S& Hochstein J.I. 1996Volume tracking of interfaces having surface tension in two and three dimensions. AIAA96–0859. Google Scholar
Kuharsky A.L& Fogelson A.L. 2001Surface-mediated control of blood coagulation: the role of binding site densities and platelet deposition. Biophys. J. 80, 1050–1074. Crossref, PubMed, ISI, Google Scholar
Lai, M. 1994 An approximate projection method for reacting flow in the zero mach number limits. PhD thesis, University of California, Berkeley. Google Scholar
Lee D& Tsuei Y. 1993A hybrid adaptive gridding procedure for rescirculating fluid flow problems. J. Comput. Phys. 108, 122–141.doi:10.1006/jcph.1993.1168. . Crossref, ISI, Google Scholar
Lobanov A.I& Starozhilova T.K. 2005The effect of convective flows on blood coagulation processes. Pathophysiol. Haemost. Thromb. 34, 121–134.doi:10.1159/000089932. . Crossref, PubMed, Google Scholar
Merks R.M.H, Brodsky S.V, Goligorksy M.S, Newman S.A& Glazier J.A. 2006Cell elongation is key to in silico replication of in vitro vasculogenesis and subsequent remodeling. Dev. Biol. 289, 44–54.doi:10.1016/j.ydbio.2005.10.003. . Crossref, PubMed, ISI, Google Scholar
Newman T.J. 2005Modeling multi-cellular systems using sub-cellular elements. Math. Biosci. Eng. 2, 611–622. Crossref, ISI, Google Scholar
Newman M& Barkema GMonte Carlo methods in statistical physics. 1999Oxford, UK:Oxford University Press. Google Scholar
Newman T.J& Grima R. 2004Many-body theory of chemotactic cell–cell interactions. Phys. Rev. E. 70, 051916doi:10.1103/PhysRevE.70.051916. . Crossref, ISI, Google Scholar
Newman S, Christley S, Glimm T, Hentschel H, Kazmierczak B, Zhang Y, Zhu J& Alber M. 2008Multiscale models for vertebrate limb development. Curr. Top. Dev. Biol. 81, 311–340.doi:10.1016/S0070-2153(07)81011-8. . Crossref, PubMed, ISI, Google Scholar
Noh, W. & Woodward, P. 1967 SLIC (simple line interface calculation). In Proc. Fifth Int. Conf. on Numerical Methods for Fluid Dynamics (eds A. I. van de Vooren & P. J. Zandbergen). Springer Lecture Notes in Physics, no. 59, pp. 330–340. New York, NY: Springer. Google Scholar
Ovanesov M.V, Krasotkina J.V, Ul'yanova L.I, Abushinova K.V, Plyushch O.P, Domogatskii S.P, Vorob'ev A.I& Ataullakhanov F.I. 2002Hemophilia a and b are associated with abnormal spatial dynamics of clot growth. Biochim. Biophys. Acta. 1572, 45–57. Crossref, PubMed, ISI, Google Scholar
Palsson E& Othmer H. 2000A model for individual and collective cell movement in Dictyostelium discoideum. Proc. Natl Acad. Sci. USA. 97, 10 448–10 453.doi:10.1073/pnas.97.19.10448. . Crossref, ISI, Google Scholar
Pivkin I.V, Richardson P.D& Karniadakis G. 2006Blood flow velocity effects and role of activation delay time on growth and form of platelet thrombi. Proc. Natl Acad. Sci. USA. 103, 17 164–17 169.doi:10.1073/pnas.0608546103. . Crossref, ISI, Google Scholar
Sorensen E.N, Burgreen G.W, Wagner W.R& Antaki J.F. 1999Computational simulation of platelet deposition and activation: I. Model development and properties. Ann. Biomed. Eng. 27, 436–448.doi:10.1114/1.200. . Crossref, PubMed, ISI, Google Scholar
Sozinova O, Jiang Y, Kaiser D& Alber M. 2005A three-dimensional model of myxobacterial aggregation by contact-mediated interactions. Proc. Natl Acad. Sci. USA. 102, 11 308–11 312.doi:10.1073/pnas.0504259102. . Crossref, ISI, Google Scholar
Sozinova O, Jiang Y, Kaiser D& Alber M. 2006A three-dimensional model of myxobacterial fruiting-body formation. Proc. Natl Acad. Sci. USA. 103, 17 255–17 259.doi:10.1073/pnas.0605555103. . Crossref, ISI, Google Scholar
Zarnitsina V.I, Pokhilko A.V& Ataullakhanov F.IA mathematical model for the spatio-temporal dynamics of intrinsic pathway of blood coagulation. i. The model description. Thromb. Res. 84, 1996a225–236.doi:10.1016/S0049-3848(96)00182-X. . Crossref, PubMed, ISI, Google Scholar
Zarnitsina V.I, Pokhilko A.V& Ataullakhanov F.IA mathematical model for the spatio-temporal dynamics of intrinsic pathway of blood coagulation. ii. Results. Thromb. Res. 84, 1996b333–344.doi:10.1016/S0049-3848(96)00197-1. . Crossref, PubMed, ISI, Google Scholar

This Issue

Journal of The Royal Society Interface cover image

06 July 2008

Volume 5Issue 24

Article Information

DOI:https://doi.org/10.1098/rsif.2007.1202
PubMed:17925274
Published by:Royal Society
Online ISSN:1742-5662

History:

Manuscript received15/08/2007
Manuscript accepted18/09/2007
Published online09/10/2007
Published in print06/07/2008

License:

Citations and impact

Keywords