Inhibition of platelet-surface-bound proteins during coagulation under flow II: Antithrombin and heparin.

Blood coagulation is a self-repair process regulated by activated platelet surfaces, clotting factors, and inhibitors. Antithrombin (AT) is one such inhibitor that impedes coagulation by targeting and inactivating several key coagulation enzymes. The effect of AT is greatly enhanced in the presence of heparin, a common anticoagulant drug. When heparin binds to AT, it either bridges with the target enzyme or induces allosteric changes in AT leading to more favorable binding with the target enzyme. AT inhibition of fluid-phase enzymes caused little suppression of thrombin generation in our previous mathematical models of blood coagulation under flow. This is because in that model, flow itself was a greater inhibitor of the fluid-phase enzymes than AT. From clinical observations, it is clear that AT and heparin should have strong inhibitory effects on thrombin generation, and thus we hypothesized that AT could be inhibiting enzymes bound to activated platelet surfaces that are not subject to being washed away by flow. We extended our mathematical model to include the relevant reactions of AT inhibition at the activated platelet surfaces as well as those for unfractionated heparin and a low molecular weight heparin. Our results show that AT alone is only an effective inhibitor at low tissue factor densities, but in the presence of heparin, it can greatly alter, and in some cases shut down, thrombin generation. Additionally, we studied each target enzyme separately and found that inactivation of no single enzyme could substantially suppress thrombin generation.

Inhibition of platelet-surface-bound proteins during coagulation under flow I: TFPI.

Blood coagulation is a self-repair process regulated by activated platelet surfaces, clotting factors, and inhibitors. Tissue factor pathway inhibitor (TFPI) is one such inhibitor, well known for its inhibitory action on the active enzyme complex comprising tissue factor (TF) and activated clotting factor VII. This complex forms when TF embedded in the blood vessel wall is exposed by injury and initiates coagulation. A different role for TFPI, independent of TF:VIIa, has recently been discovered whereby TFPI binds a partially cleaved form of clotting factor V (FV-h) and impedes thrombin generation on activated platelet surfaces. We hypothesized that this TF-independent inhibitory mechanism on platelet surfaces would be a more effective platform for TFPI than the TF-dependent one. We examined the effects of this mechanism on thrombin generation by including the relevant biochemical reactions into our previously validated mathematical model. Additionally, we included the ability of TFPI to bind directly to and inhibit platelet-bound FXa. The new model was sensitive to TFPI levels and, under some conditions, TFPI could completely shut down thrombin generation. This sensitivity was due entirely to the surface-mediated inhibitory reactions. The addition of the new TFPI reactions increased the threshold level of TF needed to elicit a strong thrombin response under flow, but the concentration of thrombin achieved, if there was a response, was unchanged. Interestingly, we found that direct binding of TFPI to platelet-bound FXa had a greater anticoagulant effect than did TFPI binding to FV-h alone, but that the greatest effects occurred if both reactions were at play. The model includes activated platelets' release of FV species, and we explored the impact of varying the FV/FV-h composition of the releasate. We found that reducing the zymogen FV fraction of this pool, and thus increasing the fraction that is FV-h, led to acceleration of thrombin generation.

Computationally Driven Discovery in Coagulation.

Bleeding frequency and severity within clinical categories of hemophilia A are highly variable and the origin of this variation is unknown. Solving this mystery in coagulation requires the generation and analysis of large data sets comprised of experimental outputs or patient samples, both of which are subject to limited availability. In this review, we describe how a computationally driven approach bypasses such limitations by generating large synthetic patient data sets. These data sets were created with a mechanistic mathematical model, by varying the model inputs, clotting factor, and inhibitor concentrations, within normal physiological ranges. Specific mathematical metrics were chosen from the model output, used as a surrogate measure for bleeding severity, and statistically analyzed for further exploration and hypothesis generation. We highlight results from our recent study that employed this computationally driven approach to identify FV (factor V) as a key modifier of thrombin generation in mild to moderate hemophilia A, which was confirmed with complementary experimental assays. The mathematical model was used further to propose a potential mechanism for these observations whereby thrombin generation is rescued in FVIII-deficient plasma due to reduced substrate competition between FV and FVIII for FXa (activated factor X).

The Art and Science of Building a Computational Model to Understand Hemostasis.

Computational models of various facets of hemostasis and thrombosis have increased substantially in the last decade. These models have the potential to make predictions that can uncover new mechanisms within the complex dynamics of thrombus formation. However, these predictions are only as good as the data and assumptions they are built upon, and therefore model building requires intimate coupling with experiments. The objective of this article is to guide the reader through how a computational model is built and how it can inform and be refined by experiments. This is accomplished by answering six questions facing the model builder: (1) Why make a model? (2) What kind of model should be built? (3) How is the model built? (4) Is the model a "good" model? (5) Do we believe the model? (6) Is the model useful? These questions are answered in the context of a model of thrombus formation that has been successfully applied to understanding the interplay between blood flow, platelet deposition, and coagulation and in identifying potential modifiers of thrombin generation in hemophilia A.

Clot Permeability, Agonist Transport, and Platelet Binding Kinetics in Arterial Thrombosis.

The formation of wall-adherent platelet aggregates is a critical process in arterial thrombosis. A growing aggregate experiences frictional drag forces exerted on it by fluid moving over or through the aggregate. The magnitude of these forces is strongly influenced by the permeability of the developing aggregate; the permeability depends on the aggregate's porosity. Aggregation is mediated by formation of ensembles of molecular bonds; each bond involves a plasma protein bridging the gap between specific receptors on the surfaces of two different platelets. The ability of the bonds existing at any time to sustain the drag forces on the aggregate determines whether it remains intact or sheds individual platelets or larger fragments (emboli). We investigate platelet aggregation in coronary-sized arteries using both computational simulations and in vitro experiments. The computational model tracks the formation and breaking of bonds between platelets and treats the thrombus as an evolving porous, viscoelastic material, which moves differently from the background fluid. This relative motion generates drag forces which the fluid and thrombus exert on one another. These forces are computed from a permeability-porosity relation parameterized from experimental measurements. Basing this relation on measurements from occlusive thrombi formed in our flow chamber experiments, along with other physiological parameter values, the model produced stable dense thrombi on a similar timescale to the experiments. When we parameterized the permeability-porosity relation using lower permeabilities reported by others, bond formation was insufficient to balance drag forces on an early thrombus and keep it intact. Under high shear flow, soluble agonist released by platelets was limited to the thrombus and a boundary layer downstream, thus restricting thrombus growth into the vessel lumen. Adding to the model binding and activation of unactivated platelets through von Willebrand-factor-mediated processes allowed greater growth and made agonist-induced activation more effective.

A MATHEMATICAL MODEL OF PLATELET AGGREGATION IN AN EXTRAVASCULAR INJURY UNDER FLOW.

We present the first mathematical model of flow-mediated primary hemostasis in an extravascular injury which can track the process from initial deposition to occlusion. The model consists of a system of ordinary differential equations (ODEs) that describe platelet aggregation (adhesion and cohesion), soluble-agonist-dependent platelet activation, and the flow of blood through the injury. The formation of platelet aggregates increases resistance to flow through the injury, which is modeled using the Stokes-Brinkman equations. Data from analogous experimental (microfluidic flow) and partial differential equation models informed parameter values used in the ODE model description of platelet adhesion, cohesion, and activation. This model predicts injury occlusion under a range of flow and platelet activation conditions. Simulations testing the effects of shear and activation rates resulted in delayed occlusion and aggregate heterogeneity. These results validate our hypothesis that flow-mediated dilution of activating chemical adenosine diphosphate hinders aggregate development. This novel modeling framework can be extended to include more mechanisms of platelet activation as well as the addition of the biochemical reactions of coagulation, resulting in a computationally efficient high throughput screening tool of primary and secondary hemostasis.

Elevated hematocrit enhances platelet accumulation following vascular injury.

Red blood cells (RBCs) demonstrate procoagulant properties in vitro, and elevated hematocrit is associated with reduced bleeding and increased thrombosis risk in humans. These observations suggest RBCs contribute to thrombus formation. However, effects of RBCs on thrombosis are difficult to assess because humans and mice with elevated hematocrit typically have coexisting pathologies. Using an experimental model of elevated hematocrit in healthy mice, we measured effects of hematocrit in 2 in vivo clot formation models. We also assessed thrombin generation, platelet-thrombus interactions, and platelet accumulation in thrombi ex vivo, in vitro, and in silico. Compared with controls, mice with elevated hematocrit (RBCHIGH) formed thrombi at a faster rate and had a shortened vessel occlusion time. Thrombi in control and RBCHIGH mice did not differ in size or fibrin content, and there was no difference in levels of circulating thrombin-antithrombin complexes. In vitro, increasing the hematocrit increased thrombin generation in the absence of platelets; however, this effect was reduced in the presence of platelets. In silico, direct numerical simulations of whole blood predicted elevated hematocrit increases the frequency and duration of interactions between platelets and a thrombus. When human whole blood was perfused over collagen at arterial shear rates, elevating the hematocrit increased the rate of platelet deposition and thrombus growth. These data suggest RBCs promote arterial thrombosis by enhancing platelet accumulation at the site of vessel injury. Maintaining a normal hematocrit may reduce arterial thrombosis risk in humans.

A physics-based model for maintenance of the pH gradient in the gastric mucus layer.

It is generally accepted that the gastric mucus layer provides a protective barrier between the lumen and the mucosa, shielding the mucosa from acid and digestive enzymes and preventing autodigestion of the stomach epithelium. However, the precise mechanisms that contribute to this protective function are still up for debate. In particular, it is not clear what physical processes are responsible for transporting hydrogen protons, secreted within the gastric pits, across the mucus layer to the lumen without acidifying the environment adjacent to the epithelium. One hypothesis is that hydrogen may be bound to the mucin polymers themselves as they are convected away from the mucosal surface and eventually degraded in the stomach lumen. It is also not clear what mechanisms prevent hydrogen from diffusing back toward the mucosal surface, thereby lowering the local pH. In this work we investigate a physics-based model of ion transport within the mucosal layer based on a Nernst-Planck-like equation. Analysis of this model shows that the mechanism of transporting protons bound to the mucus gel is capable of reproducing the trans-mucus pH gradients reported in the literature. Furthermore, when coupled with ion exchange at the epithelial surface, our analysis shows that bicarbonate secretion alone is capable of neutralizing the epithelial pH, even in the face of enormous diffusive gradients of hydrogen. Maintenance of the pH gradient is found to be robust to a wide array of perturbations in both physiological and phenomenological model parameters, suggesting a robust physiological control mechanism.NEW & NOTEWORTHY This work combines modeling techniques based on physical principles, as well as novel numerical simulations to test the plausibility of one hypothesized mechanism for proton transport across the gastric mucus layer. Results show that this mechanism is able to maintain the extreme pH gradient seen in in vivo experiments and suggests a highly robust regulation mechanism to maintain this gradient in the face of dynamic lumen composition.

Synergy Between Tissue Factor and Exogenous Factor XIa in Initiating Coagulation.

OBJECTIVE: Recent evidence suggests involvement of coagulation factor XIa (FXIa) in thrombotic event development. This study was conducted to explore possible synergies between tissue factor (TF) and exogenous FXIa (E-FXIa) in thrombin generation. APPROACH AND RESULTS: In thrombin generation assays, for increasing concentrations of E-FXIa with low, but not with high TF concentrations, peak thrombin significantly increased whereas lag time and time to peak significantly decreased. Similar dependencies of lag times and rates of thrombin generation were found in mathematical model simulations. In both in vitro and in silico experiments that included E-FXIa, thrombin bursts were seen for TF levels much lower than those required without E-FXIa. For in silico thrombin bursts initiated by the synergistic action of TF and E-FXIa, the mechanisms leading to the burst differed substantially from those for bursts initiated by high TF alone. For the synergistic case, sustained activation of platelet-bound FIX by E-FXIa, along with the feedback-enhanced activation of platelet-bound FVIIIa and FXa, was needed to elicit a thrombin burst. Furthermore, the initiation of thrombin bursts by high TF levels relied on different platelet FIX/FIXa binding sites than those involved in bursts initiated by low TF levels with E-FXIa. CONCLUSIONS: Low concentrations of TF and exogenous FXIa, each too low to elicit a burst in thrombin production alone, act synergistically when in combination to cause substantial thrombin production. The observation about FIX/FIXa binding sites may have therapeutic implications.

A Mathematical Model of Venous Thrombosis Initiation.

We present a mathematical model for the initiation of venous thrombosis (VT) due to slow flow and the consequent activation of the endothelial cells (ECs) lining the vein, in the absence of overt mechanical disruption of the EC layer. It includes all reactions of the tissue factor (TF) pathway of coagulation through fibrin formation, incorporates the accumulation of blood cells on activated ECs, accounts for the flow-mediated delivery and removal of coagulation proteins and blood cells from the locus of the reactions, and accounts for the activity of major inhibitors including heparan-sulfate-accelerated antithrombin and activated protein C. The model reveals that the occurrence of robust thrombin generation (a thrombin burst) depends in a threshold manner on the density of TF on the activated ECs and on the concentration of thrombomodulin and the degree of heparan-sulfate accelerated antithrombin activity on those cells. Small changes in any of these in appropriate narrow ranges switches the response between "no burst" and "burst." The model predicts synergies among the inhibitors, both in terms of each inhibitor's multiple targets, and in terms of interactions between the different inhibitors. The model strongly suggests that the rate and extent of accumulation of activated monocytes, platelets, and MPs that can support the coagulation reactions has a powerful influence on whether a thrombin burst occurs and the thrombin response when it does. The slow rate of accumulation of cells supporting coagulation is one reason that the progress of VT is so much slower than that of arterial thrombosis initiated by subendothelial exposure.

A computational investigation of occlusive arterial thrombosis.

The generation of occlusive thrombi in stenotic arteries involves the rapid deposition of millions of circulating platelets under high shear flow. The process is mediated by the formation of molecular bonds of several distinct types between platelets; the bonds capture the moving platelets and stabilize the growing thrombi under flow. We investigated the mechanisms behind occlusive thrombosis in arteries with a two-phase continuum model. The model explicitly tracks the formation and rupture of the two types of interplatelet bonds, the rates of which are coupled with the local flow conditions. The motion of platelets in the thrombi results from competition between the viscoelastic forces generated by the interplatelet bonds and the fluid drag. Our simulation results indicate that stable occlusive thrombi form only under specific combinations for the ranges of model parameters such as rates of bond formation and rupture, platelet activation time, and number of bonds required for platelet attachment.

Platelet motion near a vessel wall or thrombus surface in two-dimensional whole blood simulations.

Computational simulations using a two-dimensional lattice-Boltzmann immersed boundary method were conducted to investigate the motion of platelets near a vessel wall and close to an intravascular thrombus. Physiological volume fractions of deformable red blood cells and rigid platelet-size elliptic particles were studied under arteriolar flow conditions. Tumbling of platelets in the red-blood-cell depleted zone near the vessel walls was strongly influenced by nearby red blood cells. The thickness of the red-blood-cell depleted zone was greatly reduced near a thrombus, and platelets in this zone were pushed close to the surface of the thrombus to distances that would facilitate their cohesion to it. The distance, nature, and duration of close platelet-thrombus encounters were influenced by the porosity of the thrombus. The strong influence on platelet-thrombus encounters of red-blood-cell motion and thrombus porosity must be taken into account to understand the dynamics of platelet attachment to a growing thrombus.

The influence of hindered transport on the development of platelet thrombi under flow.

Vascular injury triggers two intertwined processes, platelet deposition and coagulation, and can lead to the formation of an intravascular clot (thrombus) that may grow to occlude the vessel. Formation of the thrombus involves complex biochemical, biophysical, and biomechanical interactions that are also dynamic and spatially-distributed, and occur on multiple spatial and temporal scales. We previously developed a spatial-temporal mathematical model of these interactions and looked at the interplay between physical factors (flow, transport to the clot, platelet distribution within the blood) and biochemical ones in determining the growth of the clot. Here, we extend this model to include reduction of the advection and diffusion of the coagulation proteins in regions of the clot with high platelet number density. The effect of this reduction, in conjunction with limitations on fluid and platelet transport through dense regions of the clot can be profound. We found that hindered transport leads to the formation of smaller and denser clots compared to the case with no protein hindrance. The limitation on protein transport confines the important activating complexes to small regions in the interior of the thrombus and greatly reduces the supply of substrates to these complexes. Ultimately, this decreases the rate and amount of thrombin production and leads to greatly slowed growth and smaller thrombus size. Our results suggest a possible physical mechanism for limiting thrombus growth.

A study of different modeling choices for simulating platelets within the immersed boundary method.

The Immersed Boundary (IB) method is a widely-used numerical methodology for the simulation of fluid-structure interaction problems. The IB method utilizes an Eulerian discretization for the fluid equations of motion while maintaining a Lagrangian representation of structural objects. Operators are defined for transmitting information (forces and velocities) between these two representations. Most IB simulations represent their structures with piecewise linear approximations and utilize Hookean spring models to approximate structural forces. Our specific motivation is the modeling of platelets in hemodynamic flows. In this paper, we study two alternative representations - radial basis functions (RBFs) and Fourier-based (trigonometric polynomials and spherical harmonics) representations - for the modeling of platelets in two and three dimensions within the IB framework, and compare our results with the traditional piecewise linear approximation methodology. For different representative shapes, we examine the geometric modeling errors (position and normal vectors), force computation errors, and computational cost and provide an engineering trade-off strategy for when and why one might select to employ these different representations.

Towards understanding the effect of fibrinogen interactions on fibrin gel structure.

Fibrin gelation involves the enzymatic conversion of the plasma protein fibrinogen to fibrin monomers which then polymerize to form the gel that is a major structural component of a blood clot. Because fibrinogen provides the material from which fibrin is made, it is generally regarded as promoting the gelation process. However, fibrinogen can bind to a site on a fibrin oligomer, preventing another fibrin oligomer from binding there, thus slowing the polymerization process. "Soluble fibrin oligomers," which are mixtures of fibrin and fibrinogen, are found in the blood plasma and serve as biomarkers for various clotting disorders, so understanding the interplay between fibrin and fibrinogen during fibrin polymerization may have medical importance. We present a kinetic gelation model of fibrin polymerization which accounts for the dual and antagonistic roles of fibrinogen. It builds on our earlier model of fibrin polymerization that proposed a novel mechanism for branch formation, which is a necessary component of gelation. This previous model captured salient experimental observations regarding the determinants of the structure of the gel, but did not include fibrinogen binding. Here, we add to that model reactions between fibrinogen and fibrin, so oligomers are now mixtures of fibrin and fibrinogen, and characterizing their dynamics leads to equations of substantially greater complexity than previously. Using a moment generating function approach, we derive a closed system of moment equations and we track their dynamics until the finite time blow-up of specific second moments indicates that a gel has formed. In simulations begun with an initial mixture of fibrin and fibrinogen monomers, a sufficiently high relative concentration of fibrinogen prevents gelation; the critical concentration increases with the branch formation rate. In simulations begun with only fibrinogen monomers that are converted to fibrin at a specified rate, the rates of conversion, fibrinogen binding to oligomers, and branch formation together determine whether a gel forms, how long it takes to form, and the structural properties of the gel that results.

A model of fluid-structure and biochemical interactions for applications to subclinical leaflet thrombosis.

Subclinical leaflet thrombosis (SLT) is a potentially serious complication of aortic valve replacement with a bioprosthetic valve in which blood clots form on the replacement valve. SLT is associated with increased risk of transient ischemic attacks and strokes and can progress to clinical leaflet thrombosis. SLT following aortic valve replacement also may be related to subsequent structural valve deterioration, which can impair the durability of the valve replacement. Because of the difficulty in clinical imaging of SLT, models are needed to determine the mechanisms of SLT and could eventually predict which patients will develop SLT. To this end, we develop methods to simulate leaflet thrombosis that combine fluid-structure interaction and a simplified thrombosis model that allows for deposition along the moving leaflets. Additionally, this model can be adapted to model deposition or absorption along other moving boundaries. We present convergence results and quantify the model's ability to realize changes in valve opening and pressures. These new approaches are an important advancement in our tools for modeling thrombosis because they incorporate both adhesion to the surface of the moving leaflets and feedback to the fluid-structure interaction.

Blood clot formation under flow: the importance of factor XI depends strongly on platelet count.

A previously validated mathematical model of intravascular platelet deposition and tissue factor (TF)-initiated coagulation under flow is extended and used to assess the influence on thrombin production of the activation of factor XI (fXI) by thrombin and of the activation of factor IX (fIX) by fXIa. It is found that the importance of the thrombin-fXIa-fIXa feedback loop to robust thrombin production depends on the concentration of platelets in the blood near the injury. At a near-wall platelet concentration of ~250,000/µL, typical in vessels in which the shear rate is <200 s(-1), thrombin activation of fXI makes a significant difference only at low densities of exposed TF. If the near-wall platelet concentration is significantly higher than this, either because of a higher systemic platelet count or because of the redistribution of platelets toward the vessel walls at high shear rates, then thrombin activation of fXI makes a major difference even for relatively high densities of exposed TF. The model predicts that the effect of a severe fXI deficiency depends on the platelet count, and that fXI becomes more important at high platelet counts.

A Hybrid Semi-Lagrangian Cut Cell Method for Advection-Diffusion Problems with Robin Boundary Conditions in Moving Domains.

We present a new discretization approach to advection-diffusion problems with Robin boundary conditions on complex, time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. [8] to discretize the Laplace operator and Robin boundary condition. To overcome the small cell problem, we use a splitting scheme along with a semi-Lagrangian method to treat advection. We demonstrate second order accuracy in the L 1, L 2, and L ∞ norms for both analytic test problems and numerical convergence studies. We also demonstrate the ability of the scheme to convert one chemical species to another across a moving boundary.

DEVELOPMENT OF FIBRIN BRANCH STRUCTURE BEFORE AND AFTER GELATION.

In [Fogelson and Keener, Phys. Rev. E, 81 (2010), 051922], we introduced a kinetic model of fibrin polymerization during blood clotting that captured salient experimental observations about how the gel branching structure depends on the conditions under which the polymerization occurs. Our analysis there used a moment-based approach that is valid only before the finite time blow-up that indicates formation of a gel. Here, we extend our analyses of the model to include both pre-gel and post-gel dynamics using the PDE-based framework we introduced in [Fogelson and Keener, SIAM J. Appl. Math., 75 (2015), pp. 1346-1368]. We also extend the model to include spatial heterogeneity and spatial transport processes. Studies of the behavior of the model reveal different spatial-temporal dynamics as the time scales of the key processes of branch formation, monomer introduction, and diffusion are varied.

An Efficient High-Order Meshless Method for Advection-Diffusion Equations on Time-Varying Irregular Domains.

We present a high-order radial basis function finite difference (RBF-FD) framework for the solution of advection-diffusion equations on time-varying domains. Our framework is based on a generalization of the recently developed Overlapped RBF-FD method that utilizes a novel automatic procedure for computing RBF-FD weights on stencils in variable-sized regions around stencil centers. This procedure eliminates the overlap parameter δ, thereby enabling tuning-free assembly of RBF-FD differentiation matrices on moving domains. In addition, our framework utilizes a simple and efficient procedure for updating differentiation matrices on moving domains tiled by node sets of time-varying cardinality. Finally, advection-diffusion in time-varying domains is handled through a combination of rapid node set modification, a new high-order semi-Lagrangian method that utilizes the new tuning-free overlapped RBF-FD method, and a high-order time-integration method. The resulting framework has no tuning parameters and has O(N logN) time complexity. We demonstrate high-orders of convergence for advection-diffusion equations on time-varying 2D and 3D domains for both small and large Peclet numbers. We also present timings that verify our complexity estimates. Finally, we utilize our method to solve a coupled 3D problem motivated by models of platelet aggregation and coagulation, once again demonstrating high-order convergence rates on a moving domain.