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.

Mathematical modeling to understand the role of bivalent thrombin-fibrin binding during polymerization.

Thrombin is an enzyme produced during blood coagulation that is crucial to the formation of a stable clot. Thrombin cleaves soluble fibrinogen into fibrin, which polymerizes and forms an insoluble, stabilizing gel around the growing clot. A small fraction of circulating fibrinogen is the variant γA/γ', which has been associated with high-affinity thrombin binding and implicated as a risk factor for myocardial infarctions, deep vein thrombosis, and coronary artery disease. Thrombin is also known to be strongly sequestered by polymerized fibrin for extended periods of time in a way that is partially regulated by γA/γ'. However, the role of γA/γ'-thrombin interactions during fibrin polymerization is not fully understood. Here, we present a mathematical model of fibrin polymerization that considered the interactions between thrombin, fibrinogen, and fibrin, including those with γA/γ'. In our model, bivalent thrombin-fibrin binding greatly increased thrombin residency times and allowed for thrombin-trapping during fibrin polymerization. Results from the model showed that early in fibrin polymerization, γ' binding to thrombin served to localize the thrombin to the fibrin(ogen), which effectively enhanced the enzymatic conversion of fibrinogen to fibrin. When all the fibrin was fully generated, however, the fibrin-thrombin binding persisted but the effect of fibrin on thrombin switched quickly to serve as a sink, essentially removing all free thrombin from the system. This dual role for γ'-thrombin binding during polymerization led to a paradoxical decrease in trapped thrombin as the amount of γ' was increased. The model highlighted biochemical and biophysical roles for fibrin-thrombin interactions during polymerization and agreed well with experimental observations.

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.

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.

A Mathematical Model of Bivalent Binding Suggests Physical Trapping of Thrombin within Fibrin Fibers.

Thrombin is an enzyme that plays many important roles in the blood clotting process; it activates platelets, cleaves coagulation proteins within feedback loops, and cleaves fibrinogen into fibrin, which polymerizes into fibers to form a stabilizing gel matrix in and around growing clots. Thrombin also binds to the formed fibrin matrix, but this interaction is not well understood. Thrombin-fibrin binding is often described as two independent, single-step binding events, one high-affinity and one low-affinity. However, kinetic schemes describing these single-step binding events do not explain experimentally-observed residency times of fibrin-bound thrombin. In this work, we study a bivalent, sequential-step binding scheme as an alternative to the high-affinity event and, in addition to the low-affinity one. We developed mathematical models for the single- and sequential-step schemes consisting of reaction-diffusion equations to compare to each other and to experimental data. We then used Bayesian inference, in the form of Markov chain Monte Carlo, to learn model parameter distributions from previously published experimental data. For the model to best fit the data, we made an additional assumption that thrombin was irreversibly sequestered; we hypothesized that this could be due to thrombin becoming physically trapped within fibrin fibers as they formed. We further estimated that ∼30% of thrombin in the experiments to which we compare our model output became physically trapped. The notion of physically trapped thrombin may provide new insights into conflicting observations regarding the speed of fibrinolysis. Finally, we show that our new model can be used to further probe scenarios dealing with thrombin allostery.

Assessing the impact of product inhibition in a chromogenic assay.

Chromogenic substrates (CS) are synthetic substrates used to monitor the activity of a target enzyme. It has been reported that some CSs display competitive product inhibition with their target enzyme. Thus, in assays where enzyme activity is continuously monitored over long periods of time, the product inhibition may significantly interfere with the reactions being monitored. Despite this knowledge, it is rare for CSs to be directly incorporated into mathematical models that simulate these assays. This devalues the predictive power of the models. In this study, we examined the interactions between a single enzyme, coagulation factor Xa, and its chromogenic substrate. We developed, and experimentally validated, a mathematical model of a chromogenic assay for factor Xa that explicitly included product inhibition from the CS. We employed Bayesian inference, in the form of Markov-Chain Monte Carlo, to estimate the strength of the product inhibition and other sources of uncertainty such as pipetting error and kinetic rate constants. Our model, together with carefully calibrated biochemistry experiments, allowed for full characterization of the strength and impact of product inhibition in the assay. The effect of CS product inhibition in more complex reaction mixtures was further explored using mathematical models.

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.

Modeling the distribution of enzymes on lipid vesicles: A novel framework for surface-mediated reactions in coagulation.

Blood coagulation is a network of biochemical reactions wherein dozens of proteins act collectively to initiate a rapid clotting response. Coagulation reactions are lipid-surface dependent, and this dependence is thought to help localize coagulation to the site of injury and enhance the association between reactants. Current mathematical models of coagulation either do not consider lipid as a variable or do not agree with experiments where lipid concentrations were varied. Since there is no analytic rate law that depends on lipid, only apparent rate constants can be derived from enzyme kinetic experiments. We developed a new mathematical framework for modeling enzymes reactions in the presence of lipid vesicles. Here the concentrations are such that only a fraction of the vesicles harbor bound enzymes and the rest remain empty. We call the lipid vesicles with and without enzyme TF:VIIa+ and TF:VIIa- lipid, respectively. Since substrate binds to both TF:VIIa+ and TF:VIIa- lipid, our model shows that excess empty lipid acts as a strong sink for substrate. We used our framework to derive an analytic rate equation and performed constrained optimization to estimate a single, global set of intrinsic rates for the enzyme-substrate pair. Results agree with experiments and reveal a critical lipid concentration where the conversion rate of the substrate is maximized, a phenomenon known as the template effect. Next, we included product inhibition of the enzyme and derived the corresponding rate equations, which enables kinetic studies of more complex reactions. Our combined experimental and mathematical study provides a general framework for uncovering the mechanisms by which lipid mediated reactions impact coagulation processes.

Stoichiometry for entry and binding properties of the Env protein of R5 T cell-tropic HIV-1 and its evolutionary variant of macrophage-tropic HIV-1.

Human immunodeficiency virus type 1 typically requires a high density of CD4 for efficient entry as a mechanism to target CD4+ T cells (T-tropic), with CCR5 being used most often as the coreceptor. When target T cells are limiting, the virus can evolve to infect cells with a low density of CD4 such as macrophages (M-tropic). The entry phenotype is known to be encoded in the viral Env protein on the surface of the virus particle. Using data showing a dose response for infectivity based on CD4 surface density, we built a model consistent with T-tropic viruses requiring multiple CD4 molecules to mediate infection, whereas M-tropic viruses can infect cells using a single CD4 receptor molecule interaction. We also found that T-tropic viruses bound to the surface of cells with a low density of CD4 are released more slowly than M-tropic viruses which we modeled to be due to multiple interactions of the T-tropic virus with multiple CD4 molecules to allow the initial stable binding. Finally, we found that some M-tropic Env proteins, as the gp120 subunit, possess an enhanced affinity for CD4 compared with their T-tropic pair, indicating that the evolution of macrophage tropism can be reflected both in the closed Env trimer conformation on the virion surface and, in some cases, also in the open confirmation of gp120 Env. Collectively, these studies reveal differences in the stoichiometry of interaction of T-tropic and M-tropic viruses with CD4 and start to identify the basis of binding differences at the biochemical level. IMPORTANCE: Human immunodeficiency virus type 1 normally targets CD4+ T cells for viral replication. When T cells are limiting, the virus can evolve to infect myeloid cells. The evolutionary step involves a change from requiring a high surface density of CD4 for entry to being able to infect cells with a low density of CD4, as is found on myeloid lineage cells such as macrophage and microglia. Viruses able to infect macrophages efficiently are most often found in the CNS late in the disease course, and such viruses may contribute to neurocognitive impairment. Here, we examine the CD4 binding properties of the viral Env protein to explore these two different entry phenotypes.

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.

clotFoam: An open-source framework to simulate blood clot formation under arterial flow.

Blood clotting involves the coupled processes of platelet aggregation and coagulation. Simulating clotting under flow in complex geometries is challenging due to multiple temporal and spatial scales and high computational cost. clotFoam is an open-source software developed in OpenFOAM that employs a continuum model of platelet advection, diffusion, and aggregation in a dynamic fluid environment and a simplified coagulation model with proteins that advect, diffuse, and react within the fluid and with wall-bound species through reactive boundary conditions. Our framework provides the foundation on which one can build more complex models and perform reliable simulations in almost any computational domain.

clotFoam: An Open-Source Framework to Simulate Blood Clot Formation Under Arterial Flow.

Blood clotting involves the coupled processes of platelet aggregation and coagulation. Simulating clotting under flow in complex geometries is challenging due to multiple temporal and spatial scales and high computational cost. clotFoam is an open-source software developed in OpenFOAM that employs a continuum model of platelet advection, diffusion, and aggregation in a dynamic fluid environment and a simplified coagulation model with proteins that advect, diffuse, and react within the fluid and with wall-bound species through reactive boundary conditions. Our framework provides the foundation on which one can build more complex models and perform reliable simulations in almost any computational domain.

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.

Procoagulant activity in hemostasis and thrombosis: Virchow's triad revisited.

Virchow's triad is traditionally invoked to explain pathophysiologic mechanisms leading to thrombosis, alleging concerted roles for abnormalities in blood composition, vessel wall components, and blood flow in the development of arterial and venous thrombosis. Given the tissue-specific bleeding observed in hemophilia patients, it may be instructive to consider the principles of Virchow's triad when investigating mechanisms operant in hemostatic disorders as well. Blood composition (the function of circulating blood cells and plasma proteins) is the most well studied component of the triad. For example, increased levels of plasma procoagulant proteins such as prothrombin and fibrinogen are established risk factors for thrombosis, whereas deficiencies in plasma factors VIII and IX result in bleeding (hemophilia A and B, respectively). Vessel wall (cellular) components contribute adhesion molecules that recruit circulating leukocytes and platelets to sites of vascular damage, tissue factor, which provides a procoagulant signal of vascular breach, and a surface upon which coagulation complexes are assembled. Blood flow is often characterized by 2 key variables: shear rate and shear stress. Shear rate affects several aspects of coagulation, including transport rates of platelets and plasma proteins to and from the injury site, platelet activation, and the kinetics of fibrin monomer formation and polymerization. Shear stress modulates adhesion rates of platelets and expression of adhesion molecules and procoagulant activity on endothelial cells lining the blood vessels. That no one abnormality in any component of Virchow's triad fully predicts coagulopathy a priori suggests coagulopathies are complex, multifactorial, and interactive. In this review, we focus on contributions of blood composition, vascular cells, and blood flow to hemostasis and thrombosis, and suggest that cross-talk among the 3 components of Virchow's triad is necessary for hemostasis and determines propensity for thrombosis or bleeding. Investigative models that permit interplay among these components are necessary to understand the operant pathophysiology, and effectively treat and prevent thrombotic and bleeding disorders.

Illustrated State-of-the-Art Capsules of the ISTH 2022 Congress.

The ISTH London 2022 Congress is the first held (mostly) face-to-face again since the COVID-19 pandemic took the world by surprise in 2020. For 2 years we met virtually, but this year's in-person format will allow the ever-so-important and quintessential creativity and networking to flow again. What a pleasure and joy to be able to see everyone! Importantly, all conference proceedings are also streamed (and available recorded) online for those unable to travel on this occasion. This ensures no one misses out. The 2022 scientific program highlights new developments in hemophilia and its treatment, acquired and other inherited bleeding disorders, thromboinflammation, platelets and coagulation, clot structure and composition, fibrinolysis, vascular biology, venous thromboembolism, women's health, arterial thrombosis, pediatrics, COVID-related thrombosis, vaccine-induced thrombocytopenia with thrombosis, and omics and diagnostics. These areas are elegantly reviewed in this Illustrated Review article. The Illustrated Review is a highlight of the ISTH Congress. The format lends itself very well to explaining the science, and the collection of beautiful graphical summaries of recent developments in the field are stunning and self-explanatory. This clever and effective way to communicate research is revolutionary and different from traditional formats. We hope you enjoy this article and will be inspired by its content to generate new research ideas.

Minimal model of the hydrodynamical coupling of flagella on a spherical body with application to Volvox.

Flagella are hairlike appendages attached to microorganisms that allow the organisms to traverse their fluid environment. The algae Volvox are spherical swimmers with thousands of individual flagella on their surface, and their coordination is not fully understood. In this work, a previously developed minimal model of flagella synchronization is extended to the outer surface of a sphere submerged in a fluid. Each beating flagellum tip is modeled as a small sphere, elastically bound to a circular orbit just above the spherical surface and a regularized image system for Stokes flow outside of a sphere is used to enforce the no-slip condition. Biologically relevant distributions of rotors results in a rapidly developing and robust symplectic metachronal wave traveling from the anterior to the posterior of the spherical Volvox body.

A mathematical model of coagulation under flow identifies factor V as a modifier of thrombin generation in hemophilia A.

BACKGROUND: The variability in bleeding patterns among individuals with hemophilia A, who have similar factor VIII (FVIII) levels, is significant and the origins are unknown. OBJECTIVE: To use a previously validated mathematical model of flow-mediated coagulation as a screening tool to identify parameters that are most likely to enhance thrombin generation in the context of FVIII deficiency. METHODS: We performed a global sensitivity analysis (GSA) on our mathematical model to identify potential modifiers of thrombin generation. Candidates from the GSA were confirmed by calibrated automated thrombography (CAT) and flow assays on collagen-tissue factor (TF) surfaces at a shear rate of 100 per second. RESULTS: Simulations identified low-normal factor V (FV) (50%) as the strongest modifier, with additional thrombin enhancement when combined with high-normal prothrombin (150%). Low-normal FV levels or partial FV inhibition (60% activity) augmented thrombin generation in FVIII-inhibited or FVIII-deficient plasma in CAT. Partial FV inhibition (60%) boosted fibrin deposition in flow assays performed with whole blood from individuals with mild and moderate FVIII deficiencies. These effects were amplified by high-normal prothrombin levels in both experimental models. CONCLUSIONS: These results show that low-normal FV levels can enhance thrombin generation in hemophilia A. Further explorations with the mathematical model suggest a potential mechanism: lowering FV reduces competition between FV and FVIII for factor Xa (FXa) on activated platelet surfaces (APS), which enhances FVIII activation and rescues thrombin generation in FVIII-deficient blood.

A density-dependent FEM-FCT algorithm with application to modeling platelet aggregation.

Upon injury to a blood vessel, flowing platelets will aggregate at the injury site, forming a plug to prevent blood loss. Through a complex system of biochemical reactions, a stabilizing mesh forms around the platelet aggregate forming a blood clot that further seals the injury. Computational models of clot formation have been developed to a study intravascular thrombosis, where a vessel injury does not cause blood leakage outside the blood vessel but blocks blood flow. To model scenarios in which blood leaks from a main vessel out into the extravascular space, new computational tools need to be developed to handle the complex geometries that represent the injury. We have previously modeled intravascular clot formation under flow using a continuum approach wherein the transport of platelet densities into some spatial location is limited by the platelet fraction that already reside in that location, i.e., the densities satisfy a maximum packing constraint through the use of a hindered transport coefficient. To extend this notion to extravascular injury geometries, we have modified a finite element method flux-corrected transport (FEM-FCT) scheme by prelimiting antidiffusive nodal fluxes. We show that our modified scheme, under a variety of test problems, including mesh refinement, structured vs unstructured meshes, and for a range of reaction rates, produces numerical results that satisfy a maximum platelet-density packing constraint.