Key Points: All blood coagulation factors predominantly bind to a small “cap”-like region on procoagulant-activated platelets. Their concentration in this small region promotes acceleration of the membrane-dependent reactions of coagulation.

Blood platelets need to undergo activation to carry out their function of stopping bleeding. Different activation degrees lead to a stepped hierarchy of responses: ability to aggregate, granule release, and, in a fraction of platelets, phosphatidylserine (PS) exposure. This suggests the existence of decision-making mechanisms in the platelet intracellular signaling network. To identify and investigate them, we developed a computational model of PAR1-stimulated platelet signal transduction that included a minimal set of major players in the calcium signaling network. The model comprised three intracellular compartments: cytosol, dense tubular system (DTS) and mitochondria and extracellular space. Computer simulations showed that the stable resting state of platelets is maintained via a balance between calcium pumps and leaks through the DTS and plasma membranes. Stimulation of PAR1 induced oscillations in the cytosolic calcium concentrations, in good agreement with experimental observations. Further increase in the agonist level activated the mitochondrial uniporter leading to calcium uptake by mitochondria, which caused the collapse of mitochondrial membrane potential in a fraction of platelets leading to the PS exposure. The formation of this subpopulation was shown to be a stochastic process determined by the small number of activated PAR1 receptors and by heterogeneity in the number of ion pumps. These results demonstrate how a gradual increase of the activation degree can be converted into a stepped response hierarchy ultimately leading to formation of two distinct subpopulations from an initially homogeneous population.

Essentials Roles of the two thrombin receptors in platelet signaling are poorly understood. Computational systems biology modeling was used together with continuous flow cytometry. Dual‐receptor system has wide‐range sensitivity to thrombin and optimal response dynamics. Procoagulant platelet formation is determined by donor‐specific activities of the two receptors. Summary Background Activation of human platelets with thrombin proceeds via two protease‐activated receptors (PARs), PAR1 and PAR4, that have identical main intracellular signaling responses. Although there is evidence that they have different cleavage/inactivation kinetics (and some secondary variations in signaling), the reason for such redundancy is not clear. Methods We developed a multicompartmental stochastic computational systems biology model of dual‐receptor thrombin signaling in platelets to gain insight into the mechanisms and roles of PAR1 and PAR4 functioning. Experiments employing continuous flow cytometry of washed human platelets were used to validate the model and test its predictions. Activity of PAR receptors in donors was evaluated by mRNA measurement and by polymorphism sequencing. Results Although PAR1 activation produced rapid and short‐lived response, signaling via PAR4 developed slowly and propagated in time. Response of the dual‐receptor system was both rapid and prolonged in time. Inclusion of PAR1/PAR4 heterodimer formation promoted PAR4 signaling in the medium range of thrombin concentration (about 10 nm), with little contribution at high and low thrombin. Different dynamics and dose‐dependence of procoagulant platelet formation in healthy donors was associated with individual variations in PAR1 and PAR4 activities and particularly by the Ala120Thr polymorphism in the F2RL3 gene encoding PAR4. Conclusions The dual‐receptor combination is critical to produce a response combining three critical advantages: sensitivity to thrombin concentration, rapid onset and steady propagation; specific features of the protease‐activated receptors do not allow combination of all three in a single receptor.

Although reversible platelet aggregation observed in response to ADP stimulation in the presence of calcium is a well-known phenomenon, its mechanisms are not entirely clear. To study them, we developed a simple kinetic mass-action-law-based mathematical model to use it in combination with experiments. Light transmission platelet aggregometry (LTA) induced by ADP was performed for platelet-rich plasma or washed platelets using both conventional light transmission and aggregate size monitoring method based on optical density fluctuations. Parameter values of the model were determined by means of parameter estimation techniques implemented in COPASI software. The mathematical model was able to describe reversible platelet aggregation LTA curves without assuming changes in platelet aggregation parameters over time, but with the assumption that platelet can enter the aggregate only once. In the model, the mean size of platelet aggregates correlated with the solution transparency. This corresponded with flow cytometry analysis and with optical density fluctuations data on aggregate size. The predicted values of model parameters correlated with ADP concentration used in experiments. These data suggest that, at the start of the aggregation, when platelet integrins switch “on”, large unstable platelet aggregates are rapidly formed, which leads to an increase in light transmission. However, upon fragmentation of these aggregates, the probability of the post-aggregate platelets’ attachment to each other decreases preventing new aggregation and resulting in the reversible aggregation phenomenon.