Добавил:
kiopkiopkiop18@yandex.ru t.me/Prokururor I Вовсе не секретарь, но почту проверяю Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Ординатура / Хирургия / Библиотека им академика М.И. Перельмана / Книга_3720_Библиотеки_им_академика_М_И_Перельмана

.pdf
Скачиваний:
0
Добавлен:
31.08.2026
Размер:
30 Мб
Скачать
Mathematical Models of Thrombus Formation and Fibrinolysis Chapter | 5 77
https://t.me/med1917
The model was developed to further investigate spatial effects on the threshold behavior of thrombin, to explore how high near-wall platelet concentrations and shear rate affect thrombus growth, and to understand how solute transport within a developing clot affects its growth in space and time. The model tracks concentrations of coagulation proteins, inhibitors, and platelets throughout a two-dimensional domain that mimics a portion of an in vitro microuidic channel or blood vessel with an exposed patch of TF. Fluid motion is governed by the incompressible NaviereStokes equation with an additional porous media term to add frictional resistance to the uid due to the presence of bound platelets; the permeability of the bound platelets varies in space and time as the clotting process progresses. PDEs are used to simulate the transport of platelets and proteins by advection and diffusion through the uid and through the porous, heterogeneous clot. As in the KF model, platelets that adhere to the TF patch physically inhibit the activity there, but as an extension of the KF model, there is explicit treatment of the release of ADP, which can activate platelets. ADP that is released into the uid by activated platelets over a short time period after activation can activate more platelets or get carried away by ow. A near-wall excess of platelet s is also included to mimic the effects of platelet margination [96e98].
Model results showed that thrombin production demonstrates a threshold behavior with respect to TF. As with the ODE models described earlier, generation of a thrombin burst is dependent on there being enough IXa and VIIIa produced by TF:VIIa before the subendothelium is covered with platelets. The spatial distributions reveal something new about the threshold dependence. The thrombin concentrations plateau because there is little to no spatial overlap between tenase complexes and available binding sites to bring X to the tenase to become activated; Xa production is limited and so prothombinase levels, and therefore thrombin levels, plateau. The model showed that shear rate, thrombus permeability, and a near-wall excess of platelets each could signicantly affect the thrombus structure and heterogeneity, as well as the growth rate. This mathematical model has a physical analog [99,100] and correctly predicted experimental observations of platelet-dependent thrombin and brin generation in individuals with mild hemophilia A, including their response to treatment with rVIIa [100]; see the top of Fig. 5.8.
Leiderman and Fogelson extended this model to include a more complete hindered motion as a result of the porous thrombus [101]. The previous model hindered transport for platelet species in and around the forming thrombus because platelets were modeled as a continuum and the hindrance, in part, accounted for their physical size. In the extend ed version, advection and diffusion of solutes were hindered within a clot and resulted in slower growing, denser clots. The limited solute transport resulted in restriction of tenase and prothrombinase to areas deep within the interior of the clot and in restriction of their corresponding substrates from outside the clot from meeting them. This inside-out, outside-in hindrance slowed the production of thrombin and thus activation of platelets, and hence, clot growth slowed. This notion that clots have a dense core and a less dense periphery is in accordance with observations described in animal models [75,83], and clot proles from this model are in qualitative agreement with platelet proles extracted from ow assay data [80]; see
Fig. 5.9. The model is, as of this writing, being updated to simulate extravascular clot formation in a microuidic device;
the mathematical model and microuidic device are being developed simultaneously. This development is in the ea rly stages, but the uid dynamics simulated in a complex geometry show excellent agreement with the microuidic model (see the bottom of Fig. 5.8), which has demonstrated a hemostatic prototype sensitive to both coagulation and platelet function
[102]. We are developing these models in tandem to study coagulopathies and platelet dysfunction that result in excessive
blood loss; we refer the reader to other studies using microuidic devices to study bleeding disorders for more information
[103,104].
Govindarajan et al. extended the rst LF model to include brin formation by a simple assumption that brin forms at a rate directly proportional to the thrombin concentration [143]. In addition, they modied platelet diffusivity and allowed for shear-rate-dependent adhesion. Fibrin was assumed to contribute to the mechanical properties of the clot, and the resistance imparted by brin was found to be about 30 times higher than that imparted by platelets. Fibrin formation did not extend past the platelets in the clot, similar to other experimental studies by Colace et al. [80] and computational studies by Tosenberger et al. [105,106], but not in support of the hypothesis of a brin cap [107,108].
The models of clot formation under ow by Tosenberger and colleagues are based on the dissipative particle dynamics (DPD) method. With this method, the platelets and uid are modeled as mesoscale particles subjected to repulsive, dissipative, and random forces. Results of these types of models are sensitive to the choice of parameters describing such forces. The authors allowed platelets to aggregate prior to activation. Fibrin formed primarily in the core of the clots, aiding in rigidity and stability [105,106,109]. Karniadakis and colleagues have helped develop the eld of DPD in the context of platelet margination and aggregation in complex geometries [57,110,111].
Other studies have focused on intraclot transport in terms of the coreand shellstructure as proposed in the experimental study of Stalker and colleagues [83]. The role of a conning brin shell on clot growth was investigated by Kim and colleagues [107]. They demonstrated that the brin shell, prescribed with their own experimenta lly measured permeability, allowed thrombin to leak out of the clot and be washed away, but also limited exposure of platelets on the
78 Cardiovascular Thrombus
https://t.me/med1917
FIGURE 5.8 Mathematical modeling of bypass therapy in hemophilia and hemostasis. A patient with severe factor VIII deciency that had developed
inhibitors was treated with 90 mg/kg recombinant factor VIIa (rFVIIa). Top left: Overlay of platelet (blue) and brin (green) accumulation pre- and posttreatment on a collageneTF surface following perfusion for 5 min at a wall shear rate of 150/s. Top center: Fibrin (Fbn) deposition pre- and posttreatment in comparison with a healthy control in a microuidic assay (MFA). Top right: Predicted average thrombin concentration in the growing thrombus as a function of rFVIIa and FVIII levels using the spatial-temporal (ST) LeidermaneFogelson model [7]. (A and B) Passive tracer particles moving through the computational analog of a microuidic model of hemostasis shown in (C) using bright-eld imaging [102].
FIGURE 5.9 Limiting clot growth with hindered solute transport. Simulations of clot growth in the model of platelet deposition and blood coagulation
without (A and C) and with (B and D) hindered solute transport [7,45]. (A and B) The spatial distribution of prothrombin bound to the platelet-surface complex prothrombinase. These distributions show the regions where thrombin is produced; with hindered solute transport there is only a small regionof overlap between the substrate and its enzyme, and thus thrombin production and platelet activation slow the growth the clot. (C and D) The corresponding spatial distribution of bound platelets.
Mathematical Models of Thrombus Formation and Fibrinolysis Chapter | 5 79
https://t.me/med1917
clot surface, ultimately limiting clot growth. Using a clot geometry reconstructed from their own experiments, Voronov et al. measured thrombus permeability via simulation and reported a signicant difference between core and shell [112].In both studies, the shell and core structure was prescribed rather than an emergent property of thrombus growth. Stalker and colleagues looked further into the implications of a coreeshell structure and combined experiments and simulations in a series of three papers [78,79,113]. In these studies, they found that initial platelet accumulation has narrow gaps that produce an altered microenvironment that affects agonist distribution and subsequent thrombus growth. Through computational studies, they found that hindered solute diffusion becomes the dominant mode of transport within a thrombus, but that in the core, compared with the shell, molecules are completely restricted. Further, they demonstrated that intraclot transport regulates thrombin activity and thus platelet activation in vivo.
MODELS OF FIBRIN POLYMERIZATION
Only a few models of brin production, polymerization, and deposition have been developed as of this writing. Naski and Shafer developed a kinetic ODE model of brin production in the presence of AT [144]. Using their own previously determined kinetic parameters for release of brinopeptides [114,115], they proposed a kinetic scheme and corresponding system of ODEs to track the concentration of brinogen, brin, thrombin, and AT. Their model quantitatively matches the experimentally measured release of brinopeptides A and B and provides an estimate for the thrombin concentration necessary to produce a clot of known composition. This model was included in the ODE model of Chatterjee et al.,
and Shafer model monitors concentrations but does not track the polymerization/aggregation process. Weisel and Nagaswami formulated an ODE model of polymerization kinetics including the cleavage of brinopeptides, formation of protobrils, and lateral aggregation into bers [116]. Under certain length constraints of protobrils prior to aggregation, the model results are in good agreement with lag periods observed in turbidit y proles. By varying kinetic rates, the model can capture a wide variety of brin polymerization dynamics.
Fibrin gel structure is highly correlated with the thrombin concentration that induced its formation; as thrombin concentration increases, bers tend to be thinner with a higher degree of branching and form a gel that is less susceptible to brinolysis [8]. To better understand the branching process, Fogelson and Keener generalized a kinetic gelation model of linear aggregation from Ziff and Snell [117], to include clusters and branches [118]. Their proposed mechanism of brin branching led to different structures that depended on the supply rate of monomers, in qualitative accord with the experimental observations.
Since brin gels form under ow in vivo, the deposition of brin under ow must also be considered. The structure and density of a brin gel that forms under ow are very different from those of one that forms under static conditions; brin gels that form under ow are dense and have much lower water content than static gels [119]. Further, the rate of brin deposition is highly dependent on shear rate, as shown in in vitro ow assays [120]. Using a numerical simulation in conjunction with microuidic assays, Neeves et al. found that low shear rates (w10/s) allow for brin deposition to occur, but the deposition is limited to a thin boundary layer near the wall [121]. Guy and Fogelson used a kinetic gelation model to investigate brin gel formation under the in combines the kinetic gelation model with a one-dimensional (1D) PDE to allow for transport of prothrombin, thrombin, brinogen, and brin. Thrombin produced at the boundary (wall) converts brinogen into brin, and the brin polymerizes and feeds back onto the ow with a porous mediaelike term. At low shear rates, the gel growth is limited by the availability of thrombin. At high shear rates, brin is removed by ow before gelation can occur, similar to results described earlier. Interestingly, the transition between these two regimes is dictated by the permeability of the gel.
uence of a shear ow [122]. To predict the height of brin gels, their model
MODELS OF FIBRINOLYSIS
Understanding the interactions between all proteins involved in the brinolytic cascade is critical for effective thrombolysis, the clinical initiation of brinolysis. Current thrombolytic stroke treatment involves an infusion of recombinant tPA administered within 4.5 h of stroke onset. However, bleeding complications often occur posttreatment
[123], so safer, more effective thrombolytic therapies are desirable. Mathematical models of brinolysis, developed since
the late 1980s, are being used to understand observed phenomena like the accumulation of enzymes at the lysis front and the effects of blood ow, diffusion, and clot structure on lysis, and to propose changes to thrombolytic strategy.
Many models of brinolysis are 1D reactioneadvectionediffusion equations. Zidansek and Blinc studied the spatial distribution of plasmin within the clot and the progression of the bloodeclot boundary [124]. uPA is introduced into the homogeneous brin clot at a uniform ow velocity, and no distinction is made between brin-bound and plasma-phase
80 Cardiovascular Thrombus
https://t.me/med1917
proteins. Model results conrm experimental results: clots degrade about 2 orders of magnitude faster when exposed to pressure-induced uPA permeation compared with diffusion alone. A more detai led model by Diamond and Anand accounts for uPA and both brin-bound and uid-phase tPA, plasminogen, and plasmin [125]. The clot is a uniform porous gel composed of brin bers whose diameters shrink uniformly as degradation proceeds. As ber diameter shrinks, the local specic permeability changes, which in turn affects the overall permeation velocity across the clot. Model results suggest that thrombolytic therapy is effective over realistic time scales only if tPA is pressure driven; diffusion alone will not dissolve the clot in a reasonable time. Modications to the Diamond and Anand model include brinolytic inhibitors PAI-1 and a
-AP and solubilization of enzymes back into solution [126], adsorption and reaction of brinolytic enzymes from the
2
systemic circulation (to model intravenous thrombolytic therapy) [127], and the effect of outer convection on tPA-mediated lysis of brin clots under various ow conditions [128]. These models reasonably predict experimental lysis front ve­locities and show a high accumulation of tPA and plasmin at the lysis front [126], accurately predict measured lysis front positions and systemic tPA concentrations and show that the rate of lysis increases with increasing permeation [127], display an initial delay in tPA-mediated lysis followed by the formation of a lysis front that moves at a constant rate, and predict that outer convection modulates lysis at shear stresses below 10 dyn/cm 18 dyn/cm
2
[128].
2
, but has no effect at shear stresses above
Some models aim to study coagulation and brinolysis simultaneously. Anand and colleagues built a 1D reactioneadvectionediffusion model that includes the extrinsi c coagulation pathway and brinolysis in a viscoelastic uid
[129]. The brinolysis portion of the model is fairly simple: tPA is considered to be in steady state while plasminogen is
depleted (contrary to current knowledge); plasmin is created by tPA and is inhibited by a
-AP; brin (created through
2
the interaction of thrombin and brinogen) is degraded through an interaction with plasmin. A more biochemically comprehensive 1D reactionediffusion model, which includes both the extrinsic and the intrinsic coagulation pathway, treats brinolysis in a more realistic manner [130]. Although tPA is still assumed to be in steady state, the MichaeliseMenten activation of plasminogen to plasmin by tPA is an improvement from the previous model. Using the Anand et al. models as a starting point, Sequeira and colleagues developed a 3D model of coagulation and brinolysis in owing blood [131]. The model is simplied to neglect blood ow, however, and the stability of the continuum of equilibria is determined. It is hard to compare these types of models with experimental data, since the models that account for coagulation and brinolysis simultaneously are oversimplied.
A clot is not a homogeneous concentration of brin, as many early models assume, but rather a tangled mesh of brin
bers. Fibrin concentration is high in regions of the clot containing ber and nonexistent in regions of the clot betweenbers. To account for this brin heterogeneity, Bannish and colleagues developed a 1D reactionediffusion model in which
the brin concentration varies in space [132]. The model includes the major biochemical kinetics considered in earlier models (tPA converting plasminogen to plasmin, plasmin degrading brin, tPA and plasminogen binding to, and unbinding from, brin) and also considers plasmin-mediated exposure of lysine binding sites. While the model predicts patterns and rates of lysis similar to earlier models and experiments, it is unable to describe the observ ed differences in lysis of ne and coarse clots. Bannish and colleagues conclude that higher d imensional models are necessary to understand how clot structure affect s the patterns and rates of lysis.
To better study the geometry of clot dissolution, higher-dimensional models were created. Anand and colleagues modied the Diamond and Anand model [125] to account for pressure-driven permeation of plasmin (rather than tPA) into a2Dfibrin clot [133]. Results show dissolution fingersin the clot, consistent with experimental observations. A 2D random walk model by Zidansek and colleagues predicts similar lysing patterns [134]. This model uses the same kinetic equations for brinolytic enzymes (uPA, plasminogen, plasmin, a
-AP) as the Zidansek and Blinc model [124], but the
2
reactions are characterized phenomenologically by a single reaction time representing the time lag between arrival of uPA and clot degradation. For sufciently large values of this reaction time, nger-like patterns of clot dissolution are obtained. While Zidansek et al. modeled how a channel forms in the clot, Pleydell and colleagues created a 2D model of brinolysis of a recanalized blood vessel in low-velocity ow to study what happens to the remaining clot after the channel forms
[135]. The clot is assumed to be an annulus and a steady-state model tracks free plasminogen, tPA, plasmin, a
-AP,
2
brinogen, and brin degradation products, as well as bound tPA, plasminogen, and plasmin. Model results show that the rate of clot degradation will change as the clot dissolves because the kinetic reactions depend on ow parameters, which change as the channel widens. Sersa and colleagues extended the ideas of Pleydell et al. to study thrombolysis of nonocclusive clots in high-velocity axially directed blood ow that is either strictly laminar [136] or allowed to be turbulent
[137]. As in Pleydell et al., there is a channel through the clot, and lysis is modeled as the radial expansion of the channel.
Also, biochemical reactions are not explicitly modeled, but rather summarized phenomenologically in two parameters representing time lags in the brinolytic system. By comparing model results with experimental data, the authors show that fast ow increases the clot dissolution rate due to the work of mechanical forces of streaming blood, not just due to better
Mathematical Models of Thrombus Formation and Fibrinolysis Chapter | 5 81
https://t.me/med1917
permeation of brinolytic enzymes into the clot. As previous models treated clots homogeneously, Bajd and Sersa developed a 3D model of clot dissolution that accounts for complex microscopic structure by dening the clot to be a collection of randomly distributed spherical blood cells connected by elastic brin bonds [138]. These bonds can be corroded via biochemical degradation or eroded via mechanical forces of streaming blood. As in Sersa et al. [136], the authors nd that clot dissolution is enhanced in high-velocity ow, and also that the size of the removed clot fragments depends on ow velocity. Piebalgs and Xu developed a more complete 2D model that describes thrombol ysis of the initial occlusive clot through the recanalized mural clot [139]. The model is a higher-dimensional extension/modication of the Diamond and Anand and Wootton et al. models [125,128], which explicitly includes equations for basic brinolytic reactions and includes both inner and outer convection. Results of this model show lysis patterns similar to those in Anand et al. and Zidansek et al. [133,134], and suggest that both the lysis pattern and the lysis time depend on the pressure drop across the clot.
The most biochemically and structurally detailed model as of this writing is the 3D stochastic multiscale model from Bannish and colleagues [140,141]. A microscale model of a ber cross section is used to model degradation of a single ber. Binding sites are distributed throughout the ber cross section in accordance with known protobril spacing. tPA can convert plasminogen to plasmin, which can degrade binding sites or expose new sites. Each protobrilcross section has six total binding sites (the exposed site plus ve initial ly cryptic sites that can be exposed by plasmin), representing the six chains of two overlapping brin monomers. Plasmin can crawl to neighboring binding sites, and the ber is considered degraded when 2/3 of the binding sites within the cross section have been degraded by plasmin. In this way, single-ber degradation is accurately modeled as transverse cutting.In the macroscale model, a brin clot is modeled as a 3D square lattice, with each lattice edge representing a brin ber. The diameter of the bers and the pore size between them are adjustable. A bolus of tPA molecules is introduced into a brin-free region abutting the clot, and the model tracks each molecule as it binds to and unbinds from brin and diffuses through the clot. When a tPA molecule binds, a tPA leaving timeand ber degradation time are calculated from distributions obtained from the microscale model. With this stochastic multiscale model, it is possible to study the lysis of clots of different structure. In fact, Bannish et al. were able to reconcile the conict in the literature about whether coarse clots degrade faster or slower than ne clots [140,141]. They show that the number of tPA molecules in the system determines what type of clot degrades faster (Fig. 5.10). Coarse clots have an advantage when there are few er tPA molecules, because coarse clots have fewer bers. Even though those bers take longer to degrade, degradation can start on all bers at the clot front. When there is a lot of tPA present, ne clots degrade faster because there is enough tPA to start degradation of all bers at the clot front, and the individual bers composing the ne clot degrade more quickly than the bers in the coarse clot.
FIGURE 5.10 Lysis front velocity, in mm/min, of a ne clot (red triangles) and a coarse clot (black circles) as a function of the ratio of the number of
tissue-type plasminogen activator (tPA) molecules to the surface area of clot exposed to the brin-free region. (A) Mathematical model simulations were run with 11 different tPA-to-surface-area ratios varying from 8 to 1600 molecules/mm (B) Laboratory experiments were run with ve different tPA-to-surface-area ratios varying from 600 to 12,000 molecules/mm mean of three independent trials. Clots were formed from the pooled plasma of six donors. Reproduced from Bannish BE, Chernysh IN, Keener JP,
Fogelson AL, Weisel JW. Molecular and physical mechanisms of brinolysis and thrombolysis from mathematical modeling and experiments. Sci Rep 2017;7:6914.
2
. Each symbol is the mean of 10 independent simulations.
2
. Each symbol indicates the
82 Cardiovascular Thrombus
https://t.me/med1917
While the Bannish et al. model accounts for detailed biochemistry and brin structure, it does not include blood ow or the mechanical properties of brin bers. Models that include ow and/or mechanics do not treat the kinetics and structure in a realistic way (for example, Diamond and Anand [125] and Piebalgs and Xu [139] treat degradation as the uniform shrinking of a ber diameter, while Zidansek et al. [134] and Sersa et al. [136,137] treat the biochemical reactions phenomenologically). There is a need for a comprehensive, mechanistic model that includes detailed information about the biochemistry of lysis and the structure of brin bers and brin clots under physiological blood ow. This type of model would be an invaluable tool for innovating thrombolytic therapy.
SUMMARY
Mathematical models of blood clotting have made important contributions to the eld of hemostasis by revealing new mechanisms that regulate this complex, nonlinear process. We chose to highlight examples of such models where the real power of mathemat ical models has been exploited. Because mathematical simulations give the concentrations of every possible protein and enzyme complex at any point in space and time, one can and should use this information to systematically and seamlessly probe the system to uncover the underlying mechanisms. Too often, models are used to qualitatively reproduce experimental results without offering any mechanistic explanations of the results. We believe this is where the eld of mathematical modeling could be signicantly improved. In addition, descriptions of mathematical models can be highly technical and accessible only to the expert. It is the job of the mathematicians to analyze model results and put them into an appropriate scientic context. Better yet, mathematicians should work across disciplines to acquire other expert input on their modeling results. Multidisciplinary teams that include mathematicians, engineers, physical scientists, biologists, and clinicians are likely to generate the most successful models that will yield outputs to guide the development of therapeutic interventions and cater these interventions to individual patients.
ACKNOWLEDGMENTS
This work was supported in part by the National Institutes of Health (R01HL120728).
REFERENCES
[1] Hoffman M, Monroe DM. A cell-based model of hemostasis. Thromb Haemost 2001;85:958e65. [2] Hoffman M, Monroe DM. Coagulation 2006: a modern view of hemostasis. Hematol/Oncol Clin N 2007;21(1):1e11. [3] Monroe D, Hoffman M. What does it take to make the perfect clot? Arterioscler Thromb Vasc Biol 2006;26(1):41. [4] Monkovic DD, Tracy PB. Functional characterization of human platelet-released factor V and its activation by factor Xa and thrombin. J Biol
Chem 1990;265:17132e40. [5] Fogelson AL, Kuharsky AL. Membrane binding-site density can modulate activation thresholds in enzyme systems. J Theor Biol 1998;193:1e18. [6] Kuharsky AL, Fogelson AL. Surface-mediated control of blood coagulation: the role of binding site densities and platelet deposition. Biophys J
2001;80:1050e74. [7] Leiderman K, Fogelson AL. Grow with the ow: a spatial-temporal model of platelet deposition and blood coagulation under ow. Math Med Biol
2011;28:47e84. [8] Weisel JW. Fibrinogen and brin. In: Advances in protein chemistry [internet]. Fibrous proteins: coiled-coils, collagen and elastomers, vol. 70.
Academic Press; 2005. p. 247e99. Available from: http://www.sciencedirect.com/science/article/pii/S0065323305700085. [9] Wolberg AS. Thrombin generation and brin clot structure. Blood Rev May 2007;21(3):131e42.
[10] Blinc A, Magdic J, Fric J, Musevic I. Atomic force microscopy of brin networks and plasma clots during brinolysis. Fibrinolysis Proteolysis
2000;14(5):288e99.
[11] Veklich Y, Francis CW, White J, Weisel JW. Structural studies of brinolysis by electron microscopy. Blood 1998;92(12):4721e9. [12] Sakharov DV, Nagelkerke JF, Rijken DC. Rearrangements of the brin network and spatial distribution of brinolytic components during plasma
clot lysis. J Biol Chem 1996;271(4):2133e8.
[13] Carr ME, Alving BM. Effect of brin structure on plasmin-mediated dissolution of plasma clots. Blood Coagul Fibrinolysis 1995;6:567e73. [14] Collet JP, Soria J, Mirshahi M, Hirsch M, Dagonnet F, Caen J, et al. Dusart syndrome: a new concept of the relationship between brin clot
architecture and brin clot degradability: hypobrinolysis related to an abnormal clot structure. Blood 1993;82(8):2462e9.
[15] Kolev K, Tenekedjiev K, Komorowicz E, Machovich R. Functional evalutaion of the structural features of proteases and their substrate in brin
surface degradation. J Biol Chem 1997;272(21):13666e75.
[16] Wu JH, Siddiqui K, Diamond SL. Transport phenomena and clot dissolving therapy: an experimental investigation of diffusion-controlled and
permeation-enhanced brinolysis. Thromb Haemost 1994;72(1):105e12.
[17] Nesheim ME, Tracy RP, Mann KG. Clotspeed, a mathematical simulation of the functional properties of prothrombinase. J Biol Chem
1984;259:1447e53.
Mathematical Models of Thrombus Formation and Fibrinolysis Chapter | 5 83
https://t.me/med1917
[18] Hockin M, Jones K, Everse S, Mann K. A model for the stoichiometric regulation of blood coagulation. J Biol Chem May 2002;277(21):18322. [19] Jones KC, Mann KG. A model for the tissue factor pathway to thrombin. II. A mathematical simulation. J Biol Chem 1994;269:23367e73. [20] Butenas S, vant Veer C, Mann KG. Normalthrombin generation. Blood October 1999;94(7):2169e78. [21] Bungay SD, Gentry PA, Gentry RD. A mathematical model of lipid-mediated thrombin generation. Math Med Biol March 2003;20(1):105e29. [22] Chatterjee MS, Denney WS, Jing H, Diamond SL. Systems biology of coagulation initiation: kinetics of thrombin generation in resting and
activated human blood. PLoS Comput Biol 2010:6.
[23] Danforth CM, Orfeo T, Everse SJ, Mann KG, Brummel-Ziedins KE. Dening the boundaries of normal thrombin generation: investigations into
hemostasis. PLoS One February 2012;7(2):e30385. [24] Brummel-Ziedins K. Models for thrombin generation and risk of disease. J Thromb Haemost June 2013;11:212e23. [25] Mitrophanov AY, Rosendaal FR, Reifman J. Computational analysis of the effects of reduced temperature on thrombin generation: the
contributions of hypothermia to coagulopathy. Anesth Analg September 2013;117(3):565e74. [26] Mitrophanov AY, Rosendaal FR, Reifman J. Therapeutic correction of thrombin generation in dilution-induced coagulopathy: computational
analysis based on a data set of healthy subjects. J Trauma Acute Care Surg August 2012;73(2):S95. [27] Mitrophanov AY, Rosendaal FR, Reifman J. Computational analysis of intersubject variability and thrombin generation in dilutional coagulopathy.
Transfusion November 1, 2012;52(11):2475e86. [28] Wajima T, Isbister GK, Duffull SB. A comprehensive model for the humoral coagulation network in humans. Clin Pharmacol Ther September 1,
2009;86(3):290e8. [29] Nayak S, Lee D, Patel-Hett S, Pittman D, Martin S, Heatherington A, et al. Using a systems pharmacology model of the blood coagulation network
to predict the effects of various therapies on biomarkers. CPT Pharmacomet Syst Pharmacol July 1, 2015;4(7):396e405. [30] Zhou X, Huntjens D, Gilissen R. A systems pharmacology model for predicting effects of factor Xa inhibitors in healthy subjects: assessment of
pharmacokinetics and binding kinetics. CPT Pharmacomet Syst Pharmacol November 1, 2015;4(11):650e9. [31] Peng H, Sweeny A. Development of physiologically-based mathematical models for hemostatic resuscitation in trauma. 2016. [32] Burghaus R, Coboeken K, Gaub T, Kuepfer L, Sensse A, Siegmund H-U, et al. Evaluation of the efcacy and safety of rivaroxaban using a
computer model for blood coagulation. PLoS One April 22, 2011;6(4):e17626. [33] Willmann S, Becker C, Burghaus R, Coboeken K, Edginton A, Lippert J, et al. Development of a paediatric population-based model of the
pharmacokinetics of rivaroxaban. Clin Pharmacokinet January 1, 2014;53(1):89e102. [34] Burghaus R, Coboeken K, Gaub T, Niederalt C, Sensse A, Siegmund H-U, et al. Computational investigation of potential dosing schedules for
a switch of medication from warfarin to rivaroxabandan oral, direct Factor Xa inhibitor. Front Physiol November 7, 2014:5. Available from:
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4224077/.
[35] Diamond SL. Systems biology of coagulation. J Thromb Haemost June 1, 2013;11:224e32. [36] Brummel-Ziedins KE, Everse SJ, Mann KG, Orfeo T. Modeling thrombin generation: plasma composition based approach. J Thromb
Thrombolysis November 2013;37(1):32e44. [37] Brummel-Ziedins KE. Developing individualized coagulation proling of disease risk: thrombin generation dynamic models of the pro and
anticoagulant balance. Thromb Res May 2014;133(Suppl. 1):S9e11. [38] Panteleev MA, Sveshnikova AN, Belyaev AV, Nechipurenko DY, Gudich I, Obydenny SI, et al. Systems biology and systems pharmacology of
thrombosis. Math Model Nat Phenom January 2014;9(6):4e16. [39] Shibeko AM, Panteleev MA. Untangling the complexity of blood coagulation network: use of computational modelling in pharmacology and
diagnostics. Brief Bioinform May 1, 2016;17(3):429e39. [40] Hathcock J, Nemerson Y. Platelet deposition inhibits tissue factor activity: in vitro clots are impermeable to factor Xa. Blood 2004. [41] Okorie U, Denney WS, Chatterjee MS, Neeves KB, Diamond SL. Determination of surface tissue factor thresholds that trigger coagulation at
venous and arterial shear rates: amplication of 100 fM circulating tissue factor by ow. Blood 2008;111:3507e13. [42] Fogelson AL, Tania N. Coagulation under ow: the inuence of ow-mediated transport on the initiation and inhibition of coagulation.
Pathophysiol Haemost Thromb 2005;34(2e [43] Baugh RJ, Broze Jr GJ, Krishnaswamy S. Regulation of extrinsic pathway factor Xa formation by tissue factor pathway inhibitor. J Biol Chem
February 1998;273(8):4378e86. [44] Panteleev MA, Zarnitsina VI, Ataullakhanov FI. Tissue factor pathway inhibitor. Eur J Biochem April 1, 2002;269(8):2016e31. [45] Fogelson AL, Hussain YH, Leiderman K. Blood clot formation under ow: the importance of factor XI depends strongly on platelet count. Biophys
J January 2012;102(1):10e8. [46] Leiderman K, Chang WC, Ovanesov M, Fogelson AL. Synergy between tissue factor and exogenous factor XIa in initiating coagulation.
Arterioscler Thromb Vasc Biol January 1, 2016. Atvbaha.116.308186. [47] Roemisch JR, Kaar W, Zoechling A, Kannicht C, Putz M, Kohla G, et al. Identication of activated FXI as the major biochemical root cause in
IVIG batches associated with thromboembolic events. Analytical and experimental approaches resulting in corrective and preventive measures
implemented into the Octagam\textsuperscript [48] Ovanesov M, Shibeko AM, Woodle SA, Anderson CM, Hogwood J, Barson H, et al. Association of factor XIa with intravenous immune globulin
products implicated in thrombotic adverse events: biochemical root cause investigation. J Thromb Haemost 2011;9:272. [49] Wolberg AS, Kon RH, Monroe DM, Hoffman M. Coagulation factor XI is a contaminant in intravenous immunoglobulin preparations. Amer J
Hematol 2000;65(1):30e4.
3):91e108.
Ò
manufacturing process. 2011.
84 Cardiovascular Thrombus
https://t.me/med1917
[50] Funk MB, Gross N, Gross S, Hunfeld A, Lohmann A, Guenay S, et al. Thromboembolic events associated with immunoglobulin treatment. Vox
Sang 2013;105(1):54e64.
[51] Etscheid M, Breitner-Ruddock S, Gross S, Hunfeld A, Seitz R, Dodt J. Identication of kallikrein and FXIa as impurities in therapeutic
immunoglobulins: implications for the safety and control of intravenous blood products. Vox Sang 2012;102(1):40e6.
[52] Menis M, Sridhar G, Selvam N, Ovanesov MV, Divan HA, Liang Y, Scott D, et al. Hyperimmune globulins and same-day thrombotic adverse
events as recorded in a large healthcare database during 2008e2011. Am J Hematol 2013;88(12):1035e40. [53] Wang W, King MR. Multiscale modeling of platelet adhesion and thrombus growth. Ann Biomed Eng November 1, 2012;40(11):2345e54. [54] Fogelson AL, Neeves KB. Fluid mechanics of blood clot formation. Annu Rev Fluid Mech January 2015;47(1):377e403. [55] Hosseinzadegan H, Tafti DK. Modeling thrombus formation and growth. Biotechnol Bioeng October 1, 2017;114(10):2154e72. [56] Flamm MH, Diamond SL. Multiscale systems biology and physics of thrombosis under ow. Ann Biomed Eng November 1,
2012;40(11):2355e64. [57] Fedosov DA, Dao M, Karniadakis GE, Suresh S. Computational biorheology of human blood ow in Health and disease. Ann Biomed Eng
February 1, 2014;42(2):368e87. [58] Brass LF, Diamond SL. Transport physics and biorheology in the setting of hemostasis and thrombosis. J Thromb Haemost May 1,
2016;14(5):906e17. [59] Diamond SL. Systems analysis of thrombus formation. Circ Res April 29, 2016;118(9):1348e62. [60] Brass LF, Diamond SL, Stalker TJ. Platelets and hemostasis: a new perspective on an old subject. Blood Adv November 29, 2016;1(1):5e9. [61] Guria K, Guria GT. Spatial aspects of blood coagulation: two decades of research on the self-sustained traveling wave of thrombin. Thromb Res
March 2015;135(3):423e33. [62] Ovanesov MV, Ananyeva NM, Panteleev MA, Ataullakhanov FI, Saenko EL. Initiation and propagation of coagulation from tissue factor-bearing
cell monolayers to plasma: initiator cells do not regulate spatial growth rate. J Thromb Haemost February 1, 2005;3(2):321e31. [63] Ataullakhanov FI, Krasotkina YV, Sarbash VI, Volkova RI, Sinauridse EI, Kondratovich AY. Spatio-temporal dynamics of blood coagulation and
pattern formation: an experimental study. Int J Bifurc Chaos September 1, 2002;12(9):1969e83. [64] Zarnitsina VI, Pokhilko AV, Ataullakhanov FI. A mathematical model for the spatio-temporal dynamics of intrinsic pathway of blood coagulation.
II. Results. Thromb Res December 1, 1996;84(5):333e44. [65] Zarnitsina VI, Pokhilko AV, Ataullakhanov FI. A mathematical model for the spatio-temporal dynamics of intrinsic pathway of blood coagulation.
I. the model description. Thromb Res November 15, 1996;84(4):225e36. [66] Ataullakhanov FI, Zarnitsina VI, Pokhilko AV, Lobanov AI, Morozova OL. Spatio-temporal dynamics of blood coagulation and pattern formation:
a theoretical approach. Int J Bifurc Chaos September 1, 2002;12(9):1985e2002. [67] Panteleev MA, Ovanesov MV, Kireev DA, Shibeko AM, Sinauridze EI, Ananyeva NM, et al. Spatial propagation and localization of blood
coagulation are regulated by intrinsic and protein C pathways, respectively. Biophys J March 1, 2006;90(5):1489e500. [68] Intricate regimes of propagation of an excitation and self-organization in the blood clotting model. Phys-Uspekhi 2007;50(1):79. [69] Dashkevich NM, Ovanesov MV, Balandina AN, Karamzin SS, Shestakov PI, Soshitova NP, et al. Thrombin activity propagates in space during
blood coagulation as an excitation wave. Biophys J November 21, 2012;103(10):2233e40. [70] Beltrami E, Jesty J. The role of membrane patch size and ow in regulating a proteolytic feedback threshold on a membrane: possible application in
blood coagulation. Math Biosci July 2001;172(1):1e13. [71] Kastrup CJ, Runyon MK, Shen F, Ismagilov RF. Modular chemical mechanism predicts spatiotemporal dynamics of initiation in the complex
network of hemostasis. Proc Natl Acad Sci October 24, 2006;103(43):15747e52. [72] Kastrup CJ, Shen F, Runyon MK, Ismagilov RF. Characterization of the threshold response of initiation of blood clotting to stimulus patch size.
Biophys J October 15, 2007;93(8):2969e77. [73] Runyon MK, Kastrup CJ, Johnson-Kerner BL, Van Ha TG, Ismagilov RF. Effects of shear rate on propagation of blood clotting determined using
microuidics and numerical simulations. J Am Chem Soc March 1, 2008;130(11):3458e [74] Shen F, Kastrup CJ, Liu Y, Ismagilov RF. Threshold response of initiation of blood coagulation by tissue factor in patterned microuidic capillaries
is controlled by shear rate. Arterioscler Thromb Vasc Biol November 1, 2008;28(11):2035e41. [75] Brass LF, Wannemacher KM, Ma P, Stalker TJ. Regulating thrombus growth and stability to achieve an optimal response to injury. J Thromb
Haemost July 1, 2011;9:66e75. [76] Furie B, Furie BC. Thrombus formation in vivo. J Clin Invest December 1, 2005;115(12):3355e62. [77] Falati S, Gross P, Merrill-Skoloff G, Furie BC, Furie B. Real-time\emphin vivo imaging of platelets, tissue factor and brin during arterial
thrombus formation in a mouse. Nat Med 2002;8:1175e80. [78] Stalker TJ, Welsh JD, Tomaiuolo M, Wu J, Colace TV, Diamond SL, et al. A systems approach to hemostasis: 3. Thrombus consolidation regulates
intrathrombus solute transport and local thrombin activity. Blood June 2014. [79] Welsh JD, Stalker TJ, Voronov R, Muthard RW, Tomaiuolo M, Diamond SL, et al. A systems approach to hemostasis: 1. The interdependence of
thrombus architecture and agonist movements in the gaps between platelets. Blood June 2014. [80] Colace TV, Muthard RW, Diamond SL. Thrombus growth and embolism on tissue factor-bearing collagen surfaces under ow: role of thrombin
with and without brin. Arterioscler Thromb Vasc Biol June 1, 2012;32(6):1466e76. [81] Neeves KB, Illing DAR, Diamond SL. Thrombin ux and wall shear rate regulate brin ber deposition state during polymerization under ow.
Biophys J April 2010;98(7):1344e52.
64.
Mathematical Models of Thrombus Formation and Fibrinolysis Chapter | 5 85
https://t.me/med1917
[82] Wufsus AR, Macera NE, Neeves KB. The hydraulic permeability of blood clots as a function of brin and platelet density. Biophys J April
2013;104(8):1812e23.
[83] Stalker TJ, Traxler EA, Wu J, Wannemacher KM, Cermignano SL, Voronov R, et al. Hierarchical organization in the hemostatic response and its
relationship to the platelet-signaling network. Blood March 2013;121(10):1875e85.
[84] Welsh JD, Colace TV, Muthard RW, Stalker TJ, Brass LF, Diamond SL. Platelet-targeting sensor reveals thrombin gradients within blood clots
forming in microuidic assays and in mouse. J Thromb Haemost November 1, 2012;10(11):2344e53.
[85] Flamm MH, Colace TV, Chatterjee MS, Jing H, Zhou S, Jaeger D, et al. Multiscale prediction of patient-specic platelet function under ow.
Blood July 2012;120(1):190e8.
[86] Chatterjee MS, Purvis JE, Brass LF, Diamond SL. Pairwise agonist scanning predicts cellular signaling responses to combinatorial stimuli. Nat
Biotechnol June 2010;28(7):727e32.
[87] Flamm MH, Diamond SL, Sinno T. Lattice kinetic Monte Carlo simulations of convective-diffusive systems. J Chem Phys March
2009;130(9):94904.
[88] Flamm MH, Sinno T, Diamond SL. Simulation of aggregating particles in complex ows by the lattice kinetic Monte Carlo method. J Chem Phys
2011;134(3):34905. [89] Xu Z, Chen N, Kamocka MM, Rosen ED, Alber M. A multiscale model of thrombus development. J R Soc Interface 2008;5(24):705e22. [90] Sorensen EN, Burgreen GW, Wagner WR, Antaki JF. Computational simulation of platelet deposition and activation: I. Model development and
properties. Ann Biomed Eng 1999;27(4):436e48. [91] Graner F, Glazier JA. Simulation of biological cell sorting using a two-dimensional extended Potts model. Phys Rev Lett September 28,
1992;69(13):2013e6. [92] Glazier JA, Graner F. Simulation of the differential adhesion driven rearrangement of biological cells. Phys Rev E March 1, 1993;47(3):2128e54. [93] Xu Z, Lioi J, Alber M, Mu J, Liu X, Chen DZ, et al. Combined experimental and simulation study of blood clot formation. In: Science and
technology for humanity (TIC-STH), 2009 IEEE Toronto international conference. IEEE; 2009. p. 357e62. [94] Xu Z, Chen N, Shadden SC, Marsden JE, Kamocka MM, Rosen ED, et al. Study of blood ow impact on growth of thrombi using a multiscale
model. Soft Matter 2009;5(4):769e79. [95] Xu Z, Lioi J, Mu J, Kamocka MM, Liu X, Chen DZ, et al. A multiscale model of venous thrombus formation with surface-mediated control of
blood coagulation cascade. Biophys J 2010;98(9):1723e32. [96] Eckstein EC, Belgacem F. Model of platelet transport in owing blood with drift and diffusion terms. Biophy J 1991;60:53e69. [97] Eckstein EC, Tilles AW, Millero III FJ. Conditions for the occurrence of large near-wall excesses of small particles during blood ow. Microvasc
Res 1988;36:31e9. [98] Yeh C, Calvez AC, Eckstein EC. An estimated shape function for drift in a platelet-transport model. Biophys J 1994;67(3):1252e9. [99] Neeves KB, Maloney SF, Fong KP, Schmaier AA, Kahn ML, Brass LF, et al. Microuidic focal thrombosis model for measuring murine platelet
deposition and stability: PAR4 signaling enhances shear-resistance of platelet aggregates. J Thromb Haemost December 2008;6(12):2193e201.
[100] Onasoga-Jarvis AA, Leiderman K, Fogelson AL, Wang M, Manco-Johnson MJ, Di Paola JA, et al. The effect of factor VIII deciencies and
replacement and bypass therapies on thrombus formation under venous ow conditions in microuidic and computational models. PLoS One
November 2013;8(11):e78732.
[101] Leiderman K, Fogelson AL. The inuence of hindered transport on the development of platelet thrombi under ow. Bull Math Biol
2013;75(8):1255e83.
[102] Schoeman RM, Rana K, Danes N, Lehmann M, Paola JAD, Fogelson AL, et al. A micro
and coagulation. Cell Mol Bioeng February 1, 2017;10(1):3e15.
[103] Tasci TO, Disharoon D, Schoeman RM, Rana K, Herson PS, Marr DWM, et al. Enhanced brinolysis with magnetically powered colloidal
microwheels. Small September 1, 2017;13(36) [n/aen/a].
[104] Schoeman RM, Lehmann M, Neeves KB. Flow chamber and microuidic approaches for measuring thrombus formation in genetic bleeding
disorders. Platelets July 4, 2017;28(5):463e71.
[105] Tosenberger A, Bessonov N, Volpert V. Inuence of brinogen deciency on clot formation in ow by hybrid model. Math Model Nat Phenom
2015;10(1):36e47.
[106] Tosenberger A, Ataullakhanov F, Bessonov N, Panteleev M, Tokarev A, Volpert V. Modelling of plateletebrin clot formation in ow with a
DPDePDE method. J Math Biol February 1, 2016;72(3):649e81.
[107] Kim OV, Xu Z, Rosen ED, Alber MS. Fibrin networks regulate protein transport during thrombus development. PLoS Comput Biol June
2013;9(6):e1003095.
[108] Xu Z, Christley S, Lioi J, Kim O, Harvey C, Sun W, et al. Multiscale model of brin accumulation on the blood clot surface and platelet dynamics.
In: Asthagiri AR, Arkin AP, editors. Methods in cell biology [internet]. Computational methods in cell biology, vol. 110. Academic Press; 2012.
p. 367e88. Available from: http://www.sciencedirect.com/science/article/pii/B978012388403900014X.
[109] Tosenberger A, Ataullakhanov F, Bessonov N, Panteleev M, Tokarev A, Volpert V. Modelling of thrombus growth in ow with a DPD-PDE
method. J Theor Biol November 21, 2013;337(Suppl. C):30e41.
[110] Pivkin IV, Richardson PD, Karniadakis GE. Effect of red blood cells on platelet aggregation. IEEE Eng Med Biol Mag March 2009;28(2):32e7. [111] Yazdani A, Em Karniadakis G. Sub-cellular modeling of platelet transport in blood ow through microchannels with constriction. Soft Matter
2016;12(19):4339e51.
uidic model of hemostasis sensitive to platelet function
86 Cardiovascular Thrombus
https://t.me/med1917
[112] Voronov RS, Stalker TJ, Brass LF, Diamond SL. Simulation of intrathrombus uid and solute transport using in vivo clot structures with single
platelet resolution. Ann Biomed Eng February 2013;41(6):1297e307.
[113] Tomaiuolo M, Stalker TJ, Welsh JD, Diamond SL, Sinno T, Brass LF. A systems approach to hemostasis: 2. Computational analysis of molecular
transport in the thrombus microenvironment. Blood June 2014.
[114] Higgins DL, Lewis SD, Shafer JA. Steady state kinetic parameters for the thrombin-catalyzed conversion of human brinogen to brin. J Biol
Chem August 10, 1983;258(15):9276e82.
[115] Naski MC, Shafer JA. Alpha-thrombin-catalyzed hydrolysis of brin I. Alternative binding modes and the accessibility of the active site in brin
I-bound alpha-thrombin. J Biol Chem January 25, 1990;265(3):1401e7.
[116] Weisel JW, Nagaswami C. Computer modeling of brin polymerization kinetics correlated with electron microscope and turbidity observations:
clot structure and assembly are kinetically controlled. Biophys J 1992;63:111e28. [117] Ziff RM, Stell G. Kinetics of polymer gelation. J Chem Phys 1980;73(7):3492e9. [118] Fogelson AL, Keener JP. Toward an understanding of brin branching structure. Phys Rev E May 2010;81(5):51922. [119] Onasoga-Jarvis AA, Puls TJ, OBrien SK, Kuang L, Liang HJ, Neeves KB. Thrombin generation and brin formation under ow on biomimetic
tissue factor-rich surfaces. J Thromb Haemost March 2014;12(3):373e82. [120] Weiss HJ, Turitto VT, Baumgartner HR. Role of shear rate and platelets in promoting brin formation on rabbit subendothelium. Studies utilizing
patients with quantitative and qualitative platelet defects. J Clin Invest October 1986;78(4):1072e82. [121] Neeves KB, Illing DAR, Diamond SL. Thrombin ux and wall shear rate regulate brin ber deposition state during polymerization under ow.
Biophys J April 7, 2010;98(7):1344e52. [122] Guy RD, Fogelson AL, Keener JP. Modeling brin gel formation in a shear ow. Math Med Biol 2007;24:111e30. [123] NINDS t-PA Stroke Study. Intracerebral hemorrhage after intravenous t-PA therapy for ischemic stroke. Stroke 1997;28(11):2109e18. [124] Zidansek A, Blinc A. The inuence of transport parameters and enzyme kinetics of the brinolytic system on thrombosis: mathematical modelling
of two idealised cases. Thromb Haemost 1991;65(5):553e9. [125] Diamond SL, Anand S. Inner clot diffusion and permeation during brinolysis. Biophys J 1993;65:2622e43. [126] Anand S, Wu J, Diamond SL. Enzyme-mediated proteolysis of brous biopolymers: dissolution front movement in brin or collagen under
conditions of diffusive or convective transport. Biotechnol Bioeng 1995;48:89e107. [127] Anand S, Diamond SL. Computer simulation of systemic circulation and clot lysis dynamics during thrombolytic therapy that account for inner clot
transport and reaction. Circulation 1996;94(4):763e [128] Wootton DM, Popel AS, Alevriadou BR. An experimental and theoretical study on the dissolution of mural brin clots by tissue-type plasminogen
activator. Biotechnol Bioeng 2002;77(4):405e19. [129] Anand M, Rajagopal K, Rajagopal KR. A model for the formation and lysis of blood clots. Pathophysiol Haemost Thromb 2005;34:109e20. [130] Anand M, Rajagopal K, Rajagopal KR. A model for the formation, growth and lysis of clots in quiescent plasma. A comparison between the effects
of antithrombin III deciency and protein C deciency. J Theor Biol 2008;253:725e38. [131] Sequeira A, Santos RF, Bodnár T. Blood coagulation dynamics: mathematical modeling and stability results. Math Biosci Eng 2011;8(2):425e43. [132] Bannish BE, Keener JP, Woodbury M, Weisel JW, Fogelson AL. Modelling brinolysis: 1D continuum models. Math Med Biol
2014;31(1):45e64. [133] Anand S, Kudallur V, Pitman EB, Diamond SL. Mechanisms by which thrombolytic therapy results in nonuniform lysis and residual thrombus
after reperfusion. Ann Biomed Eng 1997;25:964e74. [134] Zidansek A, Blinc A, Lahajnar G, Keber D, Blinc R. Finger-like lysing patterns of blood clots. Biophys J 1995;69:803e9. [135] Pleydell CP, David T, Smye SW, Berridge DC. A mathematical model of post-canalization thrombolysis. Phys Med Biol 2002;47:209e24. [136] Sersa I, Vidmar J, Grobelnik B, Mikac U, Tratar G, Blinc A. Modelling the effect of laminar axially directed blood ow on the dissolution of
non-occlusive blood clots. Phys Med Biol 2007;52:2969e85. [137] Sersa I, Tratar G, Mikac U, Blinc A. A mathematical model for the dissolution of non-occlusive blood clots in fast tangential blood ow. Biorheology
2007;44:1e16. [138] Bajd F, Sersa I. Mathematical modeling of blood clot fragmentation during ow-mediated thrombolysis. Biophys J 2013;104(5):1181e90. [139] Piebalgs A, Xu XY. Towards a multi-physics modelling framework for thrombolysis under the inuence of blood ow. J R Soc Interface 2015:12. [140] Bannish BE, Chernysh IN, Keener JP, Fogelson AL, Weisel JW. Molecular and physical mechanisms of brinolysis and thrombolysis from
mathematical modeling and experiments. Sci Rep 2017;7:6914. [141] Bannish BE, Keener JP, Fogelson AL. Modelling brinolysis: a 3D stochastic multiscale model. Math Med Biol 2014;31(1):17e44. [142] Jordan SW, Chaikof EL. Simulated surface-induced thrombin generation in a ow eld. Biophysical journal 2011;101.2:276e86. [143] Govindarajan Vijay, et al. Computational study of thrombus formation and clotting factor effects under venous ow conditions. Biophysical journal
2016;110.8:1869e85. [144] Naski MC, Shafer JA. A kinetic model for the alpha-thrombin-catalyzed conversion of plasma levels of brinogen to brin in the presence of
antithrombin III. Journal of Biological Chemistry 1991;266.20:13003e10.
74.