Canalization effect in the coagulation cascade and the interindividual variability of oral anticoagulant response. a simulation Study

<p>Abstract</p> <p>Background</p> <p>Increasing the predictability and reducing the rate of side effects of oral anticoagulant treatment (OAT) requires further clarification of the cause of about 50% of the interindividual variability of OAT response that is currently unaccounted for. We explore numerically the hypothesis that the effect of the interindividual expression variability of coagulation proteins, which does not usually result in a variability of the coagulation times in untreated subjects, is unmasked by OAT.</p> <p>Results</p> <p>We developed a stochastic variant of the Hockin-Mann model of the tissue factor coagulation pathway, using literature data for the variability of coagulation protein levels in the blood of normal subjects. We simulated <it>in vitro </it>coagulation and estimated the Prothrombin Time and the INR across a model population. In a model of untreated subjects a "canalization effect" can be observed in that a coefficient of variation of up to 33% of each protein level results in a simulated INR of 1 with a clinically irrelevant dispersion of 0.12. When the mean and the standard deviation of vitamin-K dependent protein levels were reduced by 80%, corresponding to the usual Warfarin treatment intensity, the simulated INR was 2.98 ± 0.48, a clinically relevant dispersion, corresponding to a reduction of the canalization effect.</p> <p>Then we combined the Hockin-Mann stochastic model with our previously published model of population response to Warfarin, that takes into account the genetical and the phenotypical variability of Warfarin pharmacokinetics and pharmacodynamics. We used the combined model to evaluate the coagulation protein variability effect on the variability of the Warfarin dose required to reach an INR target of 2.5. The dose variance when removing the coagulation protein variability was 30% lower. The dose was mostly related to the pretreatment levels of factors VII, X, and the tissue factor pathway inhibitor (TFPI).</p> <p>Conclusions</p> <p>It may be worth exploring in experimental studies whether the pretreatment levels of coagulation proteins, in particular VII, X and TFPI, are predictors of the individual warfarin dose, even though, maybe due to a canalization-type effect, their effect on the INR variance in untreated subjects appears low.</p

Factor Xa Generation by Computational Modeling: An Additional Discriminator to Thrombin Generation Evaluation

Factor (f)Xa is a critical enzyme in blood coagulation that is responsible for the initiation and propagation of thrombin generation. Previously we have shown that analysis of computationally generated thrombin profiles is a tool to investigate hemostasis in various populations. In this study, we evaluate the potential of computationally derived time courses of fXa generation as another approach for investigating thrombotic risk. Utilizing the case (n = 473) and control (n = 426) population from the Leiden Thrombophilia Study and each individual's plasma protein factor composition for fII, fV, fVII, fVIII, fIX, fX, antithrombin and tissue factor pathway inhibitor, tissue factor-initiated total active fXa generation was assessed using a mathematical model. FXa generation was evaluated by the area under the curve (AUC), the maximum rate (MaxR) and level (MaxL) and the time to reach these, TMaxR and TMaxL, respectively. FXa generation was analyzed in the entire populations and in defined subgroups (by sex, age, body mass index, oral contraceptive use). The maximum rates and levels of fXa generation occur over a 10- to 12- fold range in both cases and controls. This variation is larger than that observed with thrombin (3–6 fold) in the same population. The greatest risk association was obtained using either MaxR or MaxL of fXa generation; with an ∼2.2 fold increased risk for individuals exceeding the 90th percentile. This risk was similar to that of thrombin generation(MaxR OR 2.6). Grouping defined by oral contraceptive (OC) use in the control population showed the biggest differences in fXa generation; a >60% increase in the MaxR upon OC use. FXa generation can distinguish between a subset of individuals characterized by overlapping thrombin generation profiles. Analysis of fXa generation is a phenotypic characteristic which may prove to be a more sensitive discriminator than thrombin generation among all individuals

Defining the Boundaries of Normal Thrombin Generation: Investigations into Hemostasis

In terms of its soluble precursors, the coagulation proteome varies quantitatively among apparently healthy individuals. The significance of this variability remains obscure, in part because it is the backdrop against which the hemostatic consequences of more dramatic composition differences are studied. In this study we have defined the consequences of normal range variation of components of the coagulation proteome by using a mechanism-based computational approach that translates coagulation factor concentration data into a representation of an individual's thrombin generation potential. A novel graphical method is used to integrate standard measures that characterize thrombin generation in both empirical and computational models (e.g max rate, max level, total thrombin, time to 2 nM thrombin (“clot time”)) to visualize how normal range variation in coagulation factors results in unique thrombin generation phenotypes. Unique ensembles of the 8 coagulation factors encompassing the limits of normal range variation were used as initial conditions for the computational modeling, each ensemble representing “an individual” in a theoretical healthy population. These “individuals” with unremarkable proteome composition was then compared to actual normal and “abnormal” individuals, i.e. factor ensembles measured in apparently healthy individuals, actual coagulopathic individuals or artificially constructed factor ensembles representing individuals with specific factor deficiencies. A sensitivity analysis was performed to rank either individual factors or all possible pairs of factors in terms of their contribution to the overall distribution of thrombin generation phenotypes. Key findings of these analyses include: normal range variation of coagulation factors yields thrombin generation phenotypes indistinguishable from individuals with some, but not all, coagulopathies examined; coordinate variation of certain pairs of factors within their normal ranges disproportionately results in extreme thrombin generation phenotypes, implying that measurement of a smaller set of factors may be sufficient to identify individuals with aberrant thrombin generation potential despite normal coagulation proteome composition

Anticoagulants and the Propagation Phase of Thrombin Generation

The view that clot time-based assays do not provide a sufficient assessment of an individual's hemostatic competence, especially in the context of anticoagulant therapy, has provoked a search for new metrics, with significant focus directed at techniques that define the propagation phase of thrombin generation. Here we use our deterministic mathematical model of tissue-factor initiated thrombin generation in combination with reconstructions using purified protein components to characterize how the interplay between anticoagulant mechanisms and variable composition of the coagulation proteome result in differential regulation of the propagation phase of thrombin generation. Thrombin parameters were extracted from computationally derived thrombin generation profiles generated using coagulation proteome factor data from warfarin-treated individuals (N = 54) and matching groups of control individuals (N = 37). A computational clot time prolongation value (cINR) was devised that correlated with their actual International Normalized Ratio (INR) values, with differences between individual INR and cINR values shown to derive from the insensitivity of the INR to tissue factor pathway inhibitor (TFPI). The analysis suggests that normal range variation in TFPI levels could be an important contributor to the failure of the INR to adequately reflect the anticoagulated state in some individuals. Warfarin-induced changes in thrombin propagation phase parameters were then compared to those induced by unfractionated heparin, fondaparinux, rivaroxaban, and a reversible thrombin inhibitor. Anticoagulants were assessed at concentrations yielding equivalent cINR values, with each anticoagulant evaluated using 32 unique coagulation proteome compositions. The analyses showed that no anticoagulant recapitulated all features of warfarin propagation phase dynamics; differences in propagation phase effects suggest that anticoagulants that selectively target fXa or thrombin may provoke fewer bleeding episodes. More generally, the study shows that computational modeling of the response of core elements of the coagulation proteome to a physiologically relevant tissue factor stimulus may improve the monitoring of a broad range of anticoagulants

Regulation of early steps of GPVI signal transduction by phosphatases: a systems biology approach

We present a data-driven mathematical model of a key initiating step in platelet activation, a central process in the prevention of bleeding following Injury. In vascular disease, this process is activated inappropriately and causes thrombosis, heart attacks and stroke. The collagen receptor GPVI is the primary trigger for platelet activation at sites of injury. Understanding the complex molecular mechanisms initiated by this receptor is important for development of more effective antithrombotic medicines. In this work we developed a series of nonlinear ordinary differential equation models that are direct representations of biological hypotheses surrounding the initial steps in GPVI-stimulated signal transduction. At each stage model simulations were compared to our own quantitative, high-temporal experimental data that guides further experimental design, data collection and model refinement. Much is known about the linear forward reactions within platelet signalling pathways but knowledge of the roles of putative reverse reactions are poorly understood. An initial model, that includes a simple constitutively active phosphatase, was unable to explain experimental data. Model revisions, incorporating a complex pathway of interactions (and specifically the phosphatase TULA-2), provided a good description of the experimental data both based on observations of phosphorylation in samples from one donor and in those of a wider population. Our model was used to investigate the levels of proteins involved in regulating the pathway and the effect of low GPVI levels that have been associated with disease. Results indicate a clear separation in healthy and GPVI deficient states in respect of the signalling cascade dynamics associated with Syk tyrosine phosphorylation and activation. Our approach reveals the central importance of this negative feedback pathway that results in the temporal regulation of a specific class of protein tyrosine phosphatases in controlling the rate, and therefore extent, of GPVI-stimulated platelet activation

Mathematical modelling indicates that lower activity of the haemostatic system in neonates is primarily due to lower prothrombin concentration

Haemostasis is governed by a highly complex system of interacting proteins. Due to the central role of thrombin, thrombin generation and specifically the thrombin generation curve (TGC), is commonly used as an indicator of haemostatic activity. Functional characteristics of the haemostatic system in neonates and children are significantly different compared with adults; at the same time plasma levels of haemostatic proteins vary considerably with age. However, relating one to the other has been difficult, both due to significant inter-individual differences for individuals of similar age and the complexity of the biochemical reactions underlying haemostasis. Mathematical modelling has been very successful at representing the biochemistry of blood clotting. In this study we address the challenge of large inter-individual variability by parameterising the Hockin-Mann model with data from individual patients, across different age groups from neonates to adults. Calculating TGCs for each patient of a specific age group provides us with insight into the variability of haemostatic activity across that age group. From our model we observe that two commonly used metrics for haemostatic activity are significantly lower in neonates than in older patients. Because both metrics are strongly determined by prothrombin and prothrombin levels are considerably lower in neonates we conclude that decreased haemostatic activity in neonates is due to lower prothrombin availability