Mathematical Techniques for Understanding Platelet Regulation and the Development of New Pharmacological Approaches

Dissecting thrombus-directed chemotaxis and random movement in neutrophil near-thrombus motion in flow chambers

BACKGROUND

Thromboinflammation is caused by mutual activation of platelets and neutrophils. The site of thromboinflammation is determined by chemoattracting agents release by endothelium, immune cells, and platelets. Impaired neutrophil chemotaxis contributes to the pathogenesis of Shwachman-Diamond syndrome (SDS). In this hereditary disorder, neutrophils are known to have aberrant chemoattractant-induced F-actin properties. Here, we aim to determine whether neutrophil chemotaxis could be analyzed using our previously developed ex vivo assay of the neutrophils crawling among the growing thrombi.

METHODS

Adult and pediatric healthy donors, alongside with pediatric patients with SDS, were recruited for the study. Thrombus formation and granulocyte movement in hirudinated whole blood were visualized by fluorescent microscopy in fibrillar collagen-coated parallel-plate flow chambers. Alternatively, fibrinogen, fibronectin, vWF, or single tumor cells immobilized on coverslips were used. A computational model of chemokine distribution in flow chamber with a virtual neutrophil moving in it was used to analyze the observed data.

RESULTS

The movement of healthy donor neutrophils predominantly occurred in the direction and vicinity of thrombi grown on collagen or around tumor cells. For SDS patients or on coatings other than collagen, the movement was characterized by randomness and significantly reduced velocities. Increase in wall shear rates to 300-500 1/s led to an increase in the proportion of rolling neutrophils. A stochastic algorithm simulating leucocyte chemotaxis movement in the calculated chemoattractant field could reproduce the experimental trajectories of moving neutrophils for 72% of cells.

CONCLUSIONS

In samples from healthy donors, but not SDS patients, neutrophils move in the direction of large, chemoattractant-releasing platelet thrombi growing on collagen.

Chemical Adjustment of Fibrinolysis

Fibrinolysis is the process of the fibrin-platelet clot dissolution initiated after bleeding has been stopped. It is regulated by a cascade of proteolytic enzymes with plasmin at its core. In pathological cases, the balance of normal clot formation and dissolution is replaced by a too rapid lysis, leading to bleeding, or an insufficient one, leading to an increased thrombotic risk. The only approved therapy for emergency thrombus lysis in ischemic stroke is recombinant tissue plasminogen activator, though streptokinase or urokinase-type plasminogen activators could be used for other conditions. Low molecular weight compounds are of great interest for long-term correction of fibrinolysis dysfunctions. Their areas of application might go beyond the hematology field because the regulation of fibrinolysis could be important in many conditions, such as fibrosis. They enhance or weaken fibrinolysis without significant effects on other components of hemostasis. Here we will describe and discuss the main classes of these substances and their mechanisms of action. We will also explore avenues of research for the development of new drugs, with a focus on the use of computational models in this field.

In Silico Hemostasis Modeling and Prediction

Computational physiology, i.e., reproduction of physiological (and, by extension, pathophysiological) processes in silico, could be considered one of the major goals in computational biology. One might use computers to simulate molecular interactions, enzyme kinetics, gene expression, or whole networks of biochemical reactions, but it is (patho)physiological meaning that is usually the meaningful goal of the research even when a single enzyme is its subject. Although exponential rise in the use of computational and mathematical models in the field of hemostasis and thrombosis began in the 1980s (first for blood coagulation, then for platelet adhesion, and finally for platelet signal transduction), the majority of their successful applications are still focused on simulating the elements of the hemostatic system rather than the total (patho)physiological response in situ. Here we discuss the state of the art, the state of the progress toward the efficient "virtual thrombus formation," and what one can already get from the existing models.

Control of Platelet CLEC-2-Mediated Activation by Receptor Clustering and Tyrosine Kinase Signaling

Platelets are blood cells responsible for vascular integrity preservation. The activation of platelet receptor C-type lectin-like receptor II-type (CLEC-2) could partially mediate the latter function. Although this receptor is considered to be of importance for hemostasis, the rate-limiting steps of CLEC-2-induced platelet activation are not clear. Here, we aimed to investigate CLEC-2-induced platelet signal transduction using computational modeling in combination with experimental approaches. We developed a stochastic multicompartmental computational model of CLEC-2 signaling. The model described platelet activation beginning with CLEC-2 receptor clustering, followed by Syk and Src family kinase phosphorylation, determined by the cluster size. Active Syk mediated linker adaptor for T cell protein phosphorylation and membrane signalosome formation, which resulted in the activation of Bruton's tyrosine kinase, phospholipase and phosphoinositide-3-kinase, calcium, and phosphoinositide signaling. The model parameters were assessed from published experimental data. Flow cytometry, total internal reflection fluorescence and confocal microscopy, and western blotting quantification of the protein phosphorylation were used for the assessment of the experimental dynamics of CLEC-2-induced platelet activation. Analysis of the model revealed that the CLEC-2 receptor clustering leading to the membrane-based signalosome formation is a critical element required for the accurate description of the experimental data. Both receptor clustering and signalosome formation are among the rate-limiting steps of CLEC-2-mediated platelet activation. In agreement with these predictions, the CLEC-2-induced platelet activation, but not activation mediated by G-protein-coupled receptors, was strongly dependent on temperature conditions and cholesterol depletion. Besides, the model predicted that CLEC-2-induced platelet activation results in cytosolic calcium spiking, which was confirmed by single-platelet total internal reflection fluorescence microscopy imaging. Our results suggest a refined picture of the platelet signal transduction network associated with CLEC-2. We show that tyrosine kinase activation is not the only rate-limiting step in CLEC-2-induced activation of platelets. Translocation of receptor-agonist complexes to the signaling region and linker adaptor for T cell signalosome formation in this region are limiting CLEC-2-induced activation as well.

Interspecies differences in protein expression do not impact the spatiotemporal regulation of glycoprotein VI mediated activation

BACKGROUND

Accurate protein quantification is a vital prerequisite for generating meaningful predictions when using systems biology approaches, a method that is increasingly being used to unravel the complexities of subcellular interactions and as part of the drug discovery process. Quantitative proteomics, flow cytometry, and western blotting have been extensively used to define human platelet protein copy numbers, yet for mouse platelets, a model widely used for platelet research, evidence is largely limited to a single proteomic dataset in which the total amount of proteins was generally comparatively higher than those found in human platelets.

OBJECTIVES

To investigate the functional implications of discrepancies between levels of mouse and human proteins in the glycoprotein VI (GPVI) signalling pathway using a systems pharmacology model of GPVI.

METHODS

The protein copy number of mouse platelet receptors was determined using flow cytometry. The Virtual Platelet, a mathematical model of GPVI signalling, was used to determine the consequences of protein copy number differences observed between human and mouse platelets.

RESULTS AND CONCLUSION

Despite the small size of mouse platelets compared to human platelets they possessed a greater density of surface receptors alongside a higher concentration of intracellular signalling proteins. Surprisingly the predicted temporal profile of Syk activity was similar in both species with predictions supported experimentally. Super resolution microscopy demonstrates that the spatial distribution of Syk is similar between species, suggesting that the spatial distribution of receptors and signalling molecules in activated platelets, rather than their copy number, is important for signalling pathway regulation.