Abstract
BACKGROUND: Because understanding of the inventory, connectivity and dynamics of the components characterizing the process of coagulation is relatively mature, it has become an attractive target for physiochemical modeling. Such models can potentially improve the design of therapeutics. The prothrombinase complex (composed of the protease factor (F)Xa and its cofactor FVa) plays a central role in this network as the main producer of thrombin, which catalyses both the activation of platelets and the conversion of fibrinogen to fibrin, the main substances of a clot. A key negative feedback loop that prevents clot propagation beyond the site of injury is the thrombin-dependent generation of activated protein C (APC), an enzyme that inactivates FVa, thus neutralizing the prothrombinase complex. APC inactivation of FVa is complex, involving the production of partially active intermediates and "protection" of FVa from APC by both FXa and prothrombin. An empirically validated mathematical model of this process would be useful in advancing the predictive capacity of comprehensive models of coagulation. RESULTS: A model of human APC inactivation of prothrombinase was constructed in a stepwise fashion by analyzing time courses of FVa inactivation in empirical reaction systems with increasing number of interacting components and generating corresponding model constructs of each reaction system. Reaction mechanisms, rate constants and equilibrium constants informing these model constructs were initially derived from various research groups reporting on APC inactivation of FVa in isolation, or in the presence of FXa or prothrombin. Model predictions were assessed against empirical data measuring the appearance and disappearance of multiple FVa degradation intermediates as well as prothrombinase activity changes, with plasma proteins derived from multiple preparations. Our work integrates previously published findings and through the cooperative analysis of in vitro experiments and mathematical constructs we are able to produce a final validated model that includes 24 chemical reactions and interactions with 14 unique rate constants which describe the flux in concentrations of 24 species. CONCLUSION: This study highlights the complexity of the inactivation process and provides a module of equations describing the Protein C pathway that can be integrated into existing comprehensive mathematical models describing tissue factor initiated coagulation.