Next Article in Journal
Divergent Effect of Cigarette Smoke on Innate Immunity in Inflammatory Bowel Disease: A Nicotine-Infection Interaction
Previous Article in Journal
Production of Injectable Marine Collagen-Based Hydrogel for the Maintenance of Differentiated Chondrocytes in Tissue Engineering Applications
Previous Article in Special Issue
The Impact of Hypoxia on Neutrophil Degranulation and Consequences for the Host
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phosphoproteomic Analysis of Rat Neutrophils Shows the Effect of Intestinal Ischemia/Reperfusion and Preconditioning on Kinases and Phosphatases

1
Laboratory of Protein Chemistry and Biochemistry, Department of Cell Biology, University of Brasilia, Brasilia 70910-900, Brazil
2
Department of Biochemistry and Molecular Biology, University of Southern Denmark, DK-5230 Odense M, Denmark
3
Laboratory of Surgical Physiopathology (LIM-62), Faculty of Medicine, University of São Paulo, São Paulo 01246903, Brazil
4
Department of Biochemistry, Albert Einstein College of Medicine, Bronx, NY 10461, USA
*
Author to whom correspondence should be addressed.
Contributed equally to this work.
Present address: Department of Pathology, Harvard Medical School; Department of Lab Medicine, Children’s Hospital Boston, Boston, MA 02115, USA.
Int. J. Mol. Sci. 2020, 21(16), 5799; https://doi.org/10.3390/ijms21165799
Submission received: 9 February 2020 / Revised: 11 March 2020 / Accepted: 17 April 2020 / Published: 13 August 2020
(This article belongs to the Special Issue The Impact of Hypoxia on Neutrophil Signaling and Function)

Abstract

:
Intestinal ischemia reperfusion injury (iIRI) is a severe clinical condition presenting high morbidity and mortality worldwide. Some of the systemic consequences of IRI can be prevented by applying ischemic preconditioning (IPC), a series of short ischemia/reperfusion events preceding the major ischemia. Although neutrophils are key players in the pathophysiology of ischemic injuries, neither the dysregulation presented by these cells in iIRI nor the protective effect of iIPC have their regulation mechanisms fully understood. Protein phosphorylation, as well as the regulation of the respective phosphatases and kinases are responsible for regulating a large number of cellular functions in the inflammatory response. Moreover, in previous work we found hydrolases and transferases to be modulated in iIR and iIPC, suggesting the possible involvement of phosphatases and kinases in the process. Therefore, in the present study, we analyzed the phosphoproteome of neutrophils from rats submitted to mesenteric ischemia and reperfusion, either submitted or not to IPC, compared to quiescent controls and sham laparotomy. Proteomic analysis was performed by multi-step enrichment of phosphopeptides, isobaric labeling, and LC-MS/MS analysis. Bioinformatics was used to determine phosphosite and phosphopeptide abundance and clustering, as well as kinases and phosphatases sites and domains. We found that most of the phosphorylation-regulated proteins are involved in apoptosis and migration, and most of the regulatory kinases belong to CAMK and CMGC families. An interesting finding revealed groups of proteins that are modulated by iIR, but such modulation can be prevented by iIPC. Among the regulated proteins related to the iIPC protective effect, Vamp8 and Inpp5d/Ship are discussed as possible candidates for control of the iIR damage.

Graphical Abstract

1. Introduction

Intestinal ischemia reperfusion (iIR) injury can be caused by various clinical conditions triggering a severe clinical syndrome and even death. These conditions include acute mesenteric ischemia, incarcerated hernia, small intestine transplantation, and intestinal obstruction [1,2]. Moreover, mesenteric ischemia may also be a consequence of hemorrhagic, cardiogenic, or septic shock [3]. The overall mortality rate due to acute mesenteric ischemia is 60% to 80% with increasing incidences over time. The iIR injury not only results in local acute inflammatory response but also pulmonary injury due to the infiltrating neutrophils and systemic inflammatory changes leading to development of multiple organ failure (MOF) [4,5,6,7,8]. Various studies have reported the key role of neutrophils in the pathophysiology of iIR as a result of leukocytes and endothelial cells (ECs) interactions in the intestinal mucosa [9]. The clinical outcome from iIR conditions is related to an interplay among neutrophils and many other cell types in the inflammatory tissue microenvironment, such as the homeostatic balance between PMN and mucosal cells in HIF regulation [10], or the correlation of neutrophil-to-lymphocyte ratio with mesenteric thrombosis [11]. Depletion of neutrophils from blood before reperfusion in human small bowel decreases the effect of IR [12]. Neutrophils contribute to the tissue damage by proteolytic enzyme secretion from cytoplasmic granules and production of free radicals by respiratory burst, causing damage to the microcirculation and prolonging ischemia [13,14,15,16].
Various strategies including pharmacological approaches were tested to reduce neutrophil infiltration to the tissues that resulted in reduced ischemia and reperfusion injury [6,17]. Neutrophil blocking or depletion results in less organ failure in the pro-inflammatory phase; however, increased organ failure incidence can be seen by severe infection in the anti-inflammatory phase. Instead of shutting down an important defense mechanism, it seems more favorable to regulate the neutrophil behavior. A phenomenon called ischemic preconditioning (IPC) was reported about 30 years ago that protects the organs against IR injury [18]. The IPC is a phenomenon triggered by short episodes of ischemia and reperfusion before a long ischemia and reperfusion. The beneficial role of IPC is clearly evident from the literature, and it has been tested in different organs including intestine, brain, spinal cord, kidney, lung, liver, and human heart [19,20,21,22,23,24,25,26]. There are some studies explaining the changes in the neutrophil behavior after IPC. Ischemia reperfusion causes systemic neutrophil activation with induced CD11b expression and neutrophil-platelet complex formation. Consequently, profound endothelial dysfunction can be observed; however, the IPC attenuates both the effects [27], even though it can be considered as a first stimulation of the inflammatory response [28,29]. The protein abundance profile of the neutrophil exposed to iIRI changes drastically, presenting alterations in phagocytosis and ROS production enzymes [30]. Some of these regulated proteins have their abundance profile kept at control levels if the IPC is applied, and most of the enzymes presenting this profile are hydrolases or transferases [31]. Less frequent, but still significant abundance changes were found for proteins related to the protein biosynthesis and proteasome degradation processes in mild surgical trauma [32]. Since protein phosphatases are a subclass of hydrolases and protein kinases are a subclass of transferases, the analysis of protein phosphorylation was developed as a follow-up study from our previous work.
Protein phosphorylation is an essential and one of the most common post-translational modifications (PTM) affecting cell cycle progression, gene expression, signal transduction, and other cellular activities [33]. Protein kinases, over 500 members, phosphorylate proteins at serine (S), threonine (T), and tyrosine (Y) residues. It has been predicted that kinases encompass 1.7% of the human genome [34]. The estimated relative abundance of S/T/Y phosphorylation found in human proteome is 86.4%/11.8%/1.8%, respectively [35]. Protein phosphatases catalyze the enzymatic removal of the phosphate group from these residues (S/T/Y) and return the proteins to their non-phosphorylated state [36]. In addition to the protein function regulation due to conformational changes, phosphorylation sometimes disrupts surfaces for protein–substrate interactions and vice versa [37,38,39]. Kinases and phosphatases regulate the complex dynamics of protein phosphorylation that can be disrupted in various diseases including diabetes, cancer, neurodegenerative, and autoimmune diseases [40,41].
In that sense, the main objective of the present work was to analyze the phosphoproteome of neutrophils from rats submitted to mesenteric ischemia and reperfusion, either submitted or not to IPC, compared to quiescent controls and sham laparotomy. In such dataset, we sought for evidences of phosphorylation-based regulation of events related to the iIRI and its prevention by iIPC. Among our findings, it is interesting to mention that most of the phosphorylation-regulated proteins are involved in neutrophil apoptosis, migration, and degranulation, and most of the regulatory kinases belong to CAMK and CMGC families. Some relevant regulatory proteins discussed are PKCδ and NDR1/Stk38, and some of them, such as Inpp5d/Ship and Vamp8, present an abundance profile of phosphopeptides that are compatible with modulation by iIR and prevention of such modulation by iIPC.

2. Results

2.1. Phosphoproteome Analysis of Rat Neutrophils

A total of 40 rats were randomly allocated into the four experimental groups. These groups included the control group (Ctrl, without any surgery), sham laparotomy group (Lap, without clamping of any artery but received the same surgical procedure), intestinal ischemia and reperfusion group (iIR, had 45 min of superior mesenteric artery occlusion (SMAO) followed by 120 min of reperfusion), and intestinal ischemic preconditioning (iIPC, subjected to 10 min of SMAO followed by 10 min of reperfusion immediately before 45 min of ischemia and 120 min of reperfusion). About 7–10 mL of blood was collected from the heart after the surgical procedures, and rat neutrophils from the four different groups (Ctrl, Lap, iIR, and iIPC) were isolated and processed for protein extraction. After the tryptic digestion, the peptides were iTRAQ labeled and multiplexed, followed by the enrichment for phosphorylated peptides. The purified peptides were HILIC fractionated and analyzed by mass spectrometry (MS) (Figure 1). The MS analysis of the peptides resulted in the identification of a total of 2151 proteins (Supplementary Table S1), out of which 549 were phosphorylated proteins (Figure 2A). For protein identification, a standard criterion of two peptides per protein was used, whereas the site localization of the phosphorylation was accepted with phosphoRS score of > 95%. The Venn diagram shown in Figure 2B,C represents regulated phosphorylated proteins and peptides. Unique and overlapping phosphosites for serine, threonine, and tyrosine (S/T/Y) were observed with phosphorylation abundance order of S > T > Y. The number of phosphorylated peptides with Sph, Tph, and Yph showed a pattern distribution of 75%, 22%, and 3%, respectively (Figure 2D), and a similar pattern was followed by the phosphosites (Figure 2E). The most common phosphorylated residue was serine (Sph) in all of the phosphorylated proteins followed by threonine (Tph) and tyrosine (Yph) phosphorylation at 60%, 33%, and 7%, respectively (Figure 2F).
A cluster-based analysis of the total proteins and phosphopeptides was performed using an in-house developed R script as described in the methodology section (Supplementary Tables S2 and S3). A total of six different clusters resulted after analysis based on the abundance profiles in the Ctrl, Lap, iIR, and iIPC groups (Figure 3A). Clusters 1 and 6 show a progressive up- and down-regulation of the phosphorylation profiles, respectively. Such a profile could be related to modifications progressively induced by long traumatic surgical procedures. Cluster 2 represents down-regulation in iIPC compared to Lap and iIR, reflecting the possible effects of iIPC alone, regardless of previous stimulation. However, clusters 3 and 5 show a pronounced effect with sharp down- or up-regulation (respectively) in iIR, and an opposite regulation in iIPC, showing proteins that might be related to the protective effect of iIPC. Cluster 4 represented a sharp decline in the expression profile in Lap. In cluster 1, 378 phosphopeptides from 188 proteins were grouped; cluster 2 carried 621 phosphopeptides from 229 proteins; cluster 3 had 541 phosphopeptides from 232 proteins; 565 phosphopeptides from 224 proteins were in cluster 4; 754 phosphopeptides from 249 proteins were in cluster 5; cluster 6 had 651 phosphopeptides from 251 proteins (Figure 3B). A different number of Sph, Tph, and Yph sites were seen across the clusters. Most of the phosphorylation regulation was observed in cluster 5, followed by clusters 6, 2, and 3.

2.2. Kinases and Phosphatases Phosphorylation in Rat Neutrophils

A total of 12 protein kinases and 4 phosphatases were identified in our analysis with significantly differential phosphorylation in at least one of the residues. Most of them contain phosphorylation in catalytic domains of the proteins. Five of the kinases identified were found with significant regulation in iIR vs. Ctrl, and two kinases were significantly regulated (FDR < 0.05) in iIR vs. iIPC. Significant differences in iIR vs. Ctrl, but not to the Lap group, might indicate proteins that were affected by tissue damage, regardless of the damage intensity, or gradually from Ctrl to Lap to iIR, but only showing significance between Ctrl and iIR. While significant differences in iIR vs. iIPC might indicate proteins involved in the protective effect of iIPC. The identified kinases and phosphatases with cluster-wise abundance profile, number and position of phosphosites, and their presence in domain region are given in Table 1. Two tyrosine-protein kinases important for the neutrophil function, fgr and syk, were found with significant regulation in 7 peptides and belonged to cluster 2. The most important phosphorylation was found in the catalytic domain at tyr400. The phosphorylation was higher in iIR compared to Ctrl, suggesting a decrease in domain phosphorylation in iIPC as well. The syk protein was found with one phosphorylation regulation in iIPC compared to Ctrl in our study and not in any domain region. Another kinase, beta-adrenergic receptor kinase 1, was found with significantly higher phosphorylation at S670 in iIR compared to Ctrl. The peptide was assigned to cluster 2, and a phosphorylation site was found in the PH domain region. Another important kinase, protein kinase c delta, was found with significant regulation in phosphorylation at position S643 and S645. The peptide was assigned to cluster 6 (lower in iIR and iIPC compared to Ctrl and Lap), and phosphorylation sites were in the catalytic domain of the protein. Another kinase with regulated phosphorylation in the catalytic domain at position T452 identified was serine/threonine-protein kinase 38. The phosphorylation was up-regulated in iIR compared to iIPC.
There are some phosphatases identified in our proteomic analysis of the rat neutrophils with significant regulation in phosphorylation in the domain regions. Tyrosine-protein phosphatase non-receptor type 6 was found with 2 phosphorylation sites with up-regulation in iIR compared to iIPC and assigned to cluster 2. The up-regulation in phosphorylation was found at Ser12 and in N-SH2 domain. Another phosphatidylinositol, 3,4,5-triphosphate 5-phosphatase 1, was identified with four up-regulated phosphorylation sites in iIR compared to iIPC and assigned to cluster 3. The phosphorylation at T519 was found in a catalytic region called SH2 domain and putative catalytic site, putative active site, putative Mg-binding site, and also in putative PI/IP-binding sites of this phosphatase.
A number of other kinases and phosphatases with significant regulation in their phosphorylation were also identified, but their phosphorylation was not found to be in any of the domain regions. These kinases and phosphatases are listed in Table 1.

2.3. Motif-x Enrichment Analysis of Phosphopeptides

The Motif X algorithm was used for the analysis of the sequence motif enrichment of phosphorylation sites (Table 2). The amino acids window of 17 and Rattus norvegicus background database were selected as parameters. Proteins with highly significant phosphorylation site at serine in their motifs showed more abundance (138) compared to the highly phosphorylated threonine residue in their motifs (8). The logos observed in Motif X are shown in the Supplementary Figure S1. Some of the proteins were observed to have both serine and threonine in their motifs with high regulation (Figure 4A).
After the motif-x analysis, distribution of the enriched serine and threonine motifs were checked across the six clusters (Figure 4B). As presented, the abundance of S.S motif was somewhat homogeneous through clusters 2, 3, 5, and 6, whereas SP motif was more prominent in these clusters. Comparatively, clusters 1 and 4 have the least motif containing phosphopeptides. The phosphopeptides with threonine motifs were prominently found in cluster 2, 3, and 5, whereas no enrichment was found for the tyrosine motifs, as the tyrosine containing peptides were low in number compared to serine- and threonine-containing phosphopeptides.

2.4. Prediction of Kinases Responsible for Regulated Phosphorylation in Domain Regions

A group-based prediction system (GPS) software package [42] was used for the phosphoproteome data. A prediction analysis was performed for the in vivo site-specific kinase-substrate relation for the proteins with regulated phosphorylation in the domain regions Table 3. The distribution of all the predicted kinase families responsible for the phosphorylation event in the serine and threonine motifs was grouped in pie charts (Figure 5). It is clear from Figure 5A that the kinase families with higher percentage for phosphorylation of serine motifs are CAMK, CMGC, and AGC, whereas for threonine motifs, they belong to CMGC among other families (Figure 5B). A list of commercially available inhibitors for these predicted kinases is given in the Supplementary Table S4.

2.5. Pathways Analysis of the Proteins with Significantly Regulated Phosphorylation

To analyze the effect of the iIPC on the regulated pathways in rat neutrophil proteome, the phosphoproteome data was searched against KEGG pathways database. The most significantly regulated pathways observed included spliceosome, chemokine signaling pathway, Fc-gamma R-mediated phagocytosis, and focal adhesion. Each pathway involved a different number of the phosphorylated proteins (Table 4).

3. Discussion

Protein phosphorylation is an essential post-translation modification, and the eukaryotic cells rely extensively on it for their basic cellular processes like metabolism, motility, organelle trafficking, muscle contraction, immunity, learning, and memory [34,43]. Signal transduction depends on the phosphorylation of the hydroxyl group on the side chain of serine, threonine, and tyrosine [44,45]. Various studies have documented different ratios of S/T/Y phosphorylation in different cell lines; Macek et al. reported a ratio of 70:20:10 in B. subtilis, and in another study, the same group found 86:12:2 for S/T/Y in human cell culture [35,46]. The serine–threonine–tyrosine ratio in our study for the rat neutrophil phosphoproteome was 75:22:3 (Figure 2B), whereas the phosphoproteomics analysis ratio of 14 different rat organs and tissues by Lundby et al. was 88.1:11.4:1.5 for S/T/Y, respectively [47]. In the present study, we analyzed the phosphoproteome of the rat neutrophil. For the relative quantification, an iTRAQ-based approach was adopted. Although presenting some limitations [48], the method is very robust for the large number of biological replicates and conditions, especially when metabolic labeling is not feasible due to the use of non-dividing cells. This labelling approach, in addition to others like TMT, allows not only the reproducible enrichment for phosphorylated peptides but also for other post-translational modifications in parallel in one experiment [49,50].
Phosphorylation events are regulated by protein kinases and phosphatases. The protein kinases play a vital role in cell signaling cascades and make up to 2% of the genome [43,51,52]. It is clear now that short linear motifs around the site of phosphorylation provide primary specificity for kinase substrate recognition [53]. The cluster 4 represents the effect of laparotomy on neutrophil phosphoproteome (Figure 3A). During our analysis, the laparotomy group did not have a statistically significant effect on phosphoproteome of the kinases and phosphatases, compared to the iIR and iIPC. There are very few proteins with significantly higher phosphorylation regulation in the laparotomy group compared to the control group, like the up-regulation of 60S acidic ribosomal protein P2 (P02401) and the down-regulation of 26S proteasome non-ATPase regulatory subunit 1 (O88761), in agreement with previous findings of higher rates of mRNA processing and protein synthesis for the neutrophil after surgical trauma [32].

3.1. Phosphorylation Events Related to the Ischemia-Reperfusion Injury

A considerable number of regulated phosphosites was associated to differences between the iIR and the Lap or the Ctrl groups, suggesting possible molecular mechanisms related to the ischemia-reperfusion injury. An important tyrosine-protein kinase, Fgr, a member of Scr family, was found with seven significantly regulated phosphorylation sites, one of them (tyr400), found in the catalytic domain, was assigned to cluster 2 with up-regulation after ischemia, compared to Ctrl. Phosphorylation in the activation loop of the kinase domain at this position switches the kinase to its active conformation and regulates migratory ability and chemokine secretion of the neutrophil [54,55,56]. The up-regulation of the phosphorylation in the catalytic domain in iIR suggests Fgr activation in neutrophils regulating its migration and chemokine secretion after iIR. Interestingly, in our analysis, five regulated phosphopeptides were clustered in cluster 5 containing three motifs, SP, S.S, and S…….S (Table 1 and Table 2). This makes the Fgr kinase a strong candidate for modulation of the neutrophil function with multiple catalytic functions contributed by the distant motif sequences of the protein [57]. Proto-oncogene tyrosine-protein kinase, Src, is responsible for the specific phosphorylation of tyrosine residues in other proteins [58]. In our search analysis, it is the predicted kinase with the highest score that might be a possible candidate responsible for the Fgr phosphorylation at Y400 (Table 3). Another relevant kinase for the possible Y400 phosphorylation is tyrosine-protein kinase, Tec, which plays an important role in human neutrophil responses to monosodium urate [59]. The high score of the predicted kinase indicates a high possibility of the residue to be phosphorylated by that kinase. Although some of the predicted kinases like Tec do not have proper inhibitors, others, such as Src, can be further tested by inhibition of the phosphorylation event to study the neutrophil migration and induced damage to the host tissues or organs after ischemia and reperfusion (Supplementary Table S4).
The β-adrenergic receptor kinase-1 (Grk2/P26817) has a C-terminal PH domain that possesses binding specificity for Gβγ subunits and for the membrane phospholipid PIP2 [60,61,62]. We identified a serine residue at S670 with SP motif in the PH domain of this protein with significantly higher phosphorylation in iIR compared to Ctrl. Its higher expression in neutrophils desensitization to chemoattractants was observed in sepsis and in neutrophil treatment with cytokines plus LPS [63]. In our analysis, its regulation in phosphorylation may affect the sensation of neutrophil towards chemoattractant and hence their migration. Among the five possible different kinases responsible for this specific S670 phosphorylation in Grk2, CMGC/MAPK has the highest score and known inhibitors. It has been shown previously that S670 phosphorylation by ERK1/2 kinases impaired strongly Grk2/Gβγ interaction, inhibiting its catalytic activity and kinase translocation towards membrane receptor substrates. However, phosphorylation of the same serine residue by MAPK kinase disrupts the complex allowing Grk2/GIT1 association and favors cell migration preceded by the localized activation of Rac/PAK/MEK/ERK signaling [64]. Our finding shows that Grk2 phosphorylation might be because of MAPK or ERK1/2 that can alter Grk2 catalytic activity in neutrophils after ischemia towards membrane receptor substrates and consequently activation of various cascades via GPCR signaling.
Protein kinase C delta (PKCδ) is a novel PKC that plays an important role in neutrophil transmigration mediated by IL-1 beta and fMLP by regulating neutrophil adherence and was identified in our analysis with two significantly regulated phosphorylation sites (S643 and S 645) in the catalytic domain in iIPC versus control. [65]. S643 phosphorylation has already been shown to be required for catalytic maturation of PKCδ, whereas in neutrophils the role is not known. Enhanced EC survival was observed after decreased activity of this enzyme by hypoxia [66]. In the present work, two kinases, CAMK/RAD53 and CAMK/CAMK1, were predicted with the highest score for the regulated phosphosites S643 and S645, respectively. These multifunctional CAMKs belong to a family of serine/threonine kinases that is sensitive to intracellular Ca2+ and affects various cellular activities like gene expression, cell cycle, differentiation, and ischemic tolerance [67].
Serine/threonine-protein kinase (OSR1) and serine/threonine-protein kinase 38 (Stk38) are the next two important kinases we identified with significantly higher phosphorylation in iIR compared to iIPC. The Stk38, also known as Nuclear Dbf2-Related kinase 1 (NDR1), has a regulated phosphosite on threonine residue T452 in the catalytic domain region (Table 1). Selimoglu et al. found that osmotic pressure and oxidative stress result in NDR1 phosphorylation on T444 by MAP4k4 kinase and induce apoptosis [68,69]. The possible kinases for T452 site phosphorylation are tyrosine kinase like (TKL) and PKC, and the exact role of this phosphosite is not known yet. AGC kinase family including PKC mainly regulates different physiological processes important for metabolism, cell growth, proliferation, and survival [70]. In neutrophils, phosphorylation of NDR1 in the catalytic domain can be a good target during ischemia reperfusion that consequently might lower tissue damage, whereas in iIPC, the phosphorylation was significantly lower compared to iIR, validating its protective role.
Significantly regulated phosphopeptides were identified in some other kinases in our analysis, but the respective phosphosites were not localized in any domain region of the kinases. Serine/threonine kinase N1 (Pkn1) is one of the members of protein kinase C superfamily of serine/threonine kinases and was found with up-regulation in phosphorylation in iIR compared to Ctrl at serine residue S377 with S.S motif. The serine/threonine kinase family is one of the first identified effectors for RhoA-GTPase that exists in integral plasma membrane and cytosolic pool [71]. The spleen tyrosine kinase (Syk) was found with significant down-regulation in S.S motif region at serine S291 residue in iIR compared to Ctrl group during our analysis. It has been found that CD18 integrins downstream signaling is dependent on Syk, and loss of Syk kinase mediates integrin signaling impaired activation of leukocytes resulting in reduced host defense responses [72]. For neutrophils, Syk is an important player and required for phagocytosis, and its inhibition lowers the phagocytosis of IgG coated particles [73]. Some studies show its role in making the actin filament cup during FCγR-mediated phagocytosis [74]. The significant down-regulation of the phosphorylation on S291 site might affect actin filament formation or defense response after iIPC. Another kinase, STE20-like serine/threonine-protein kinase (Slk), found with a significantly regulated phosphosite, exacerbates apoptotic activity and may regulate survival of the cell during repair or injury [75]. There are some studies showing the apoptotic correlation of Slk expression in kidney cells subjected to in vitro ischemia reperfusion injury [76]. However, other findings show its effect on cytoskeletal reorganization as part of anti-apoptotic cascades [77,78]. Oxsr1 is an oxidative-stress responsive protein and belongs to the STE20 family. It regulates the Na+/H+ exchanger activity [79]. In our analysis, significantly up-regulated phosphorylation was found for Oxsr1 at serine residues S324 and S359 in SP motif in iIR surgical group (Table 1), suggesting its role in reducing the production of reactive oxygen species (ROS) and NADPH oxidase inhibition [80]. Pak2 was also identified with two phosphorylated peptides regulated only in iIR compared to Ctrl and Lap, and it has been shown as one of the most regulated proteins for actin cytoskeleton regulation pathway playing a key role in directional migration [31].
It has been well documented that phosphatases cause de-phosphorylation of the phosphorylated proteins. Some important protein phosphatases with significant regulation in their phosphorylation have been identified in our analysis (Table 1). Phosphatases have regulatory sequences that flank the catalytic domain and control the activity by interactions at the active sites [81,82]. Protein tyrosine phosphatase receptor type C (Ptprc/CD45) is a phosphatase identified in our analysis with one phosphosite on S1209 significantly regulated in iIR compared to Lap and iIPC groups. It is commonly known as LCA (leukocyte common antigen) and has an important role in neutrophil adhesion, chemotaxis, ROS production, phagocytosis, and bacterial killing [83]. Another phosphatase called tyrosine-protein phosphatase non-receptor type 6 (Ptpn6/Shp-1) was found with an up-regulated phosphosite on S12 in SH2_N-SH2_SHP_like domain and assigned to cluster 2. The N-terminal SH2 domain regulates Shp-1 catalytic activity negatively by binding directly to Shp-1 catalytic domain. Shp-1 also modulates apoptosis signaling cascades in neutrophils [84]. Specific surface immune inhibitory receptors containing immunoreceptor tyrosine-based inhibitory motifs (ITIM) in their cytoplasmic domain mediate inhibitory signaling [85]. Following receptor activation, the Src family kinases phosphorylate a tyrosine residue of the ITIM domain and recruit SH2 domain-containing tyrosine phosphatases (Shp-1). Their subsequent phosphorylation results in the deactivation of tyrosine kinases inhibiting the survival signaling, especially TNF-α, TRAIL, Fas ligand, G-CSF, GM-CSF, or IFN-γ in neutrophils [86]. Shp-1 activation in inhibition of survival signaling leads to longer lifespan found in activated neutrophils after iIR surgical procedures. The identified significant phosphosite regulation on S12 in N-Sh2 domain of Shp-1 can be an important possible candidate for regulating these mechanisms and needs further validation. Different kinases have been reported to bind to Shp-1 through a novel phosphorylation-independent kinase tyrosyl inhibitory motif [87]. We identified significant regulation in phosphosite S12, in the R..S motif having binding possibility for kinases from CAMK and CMGC families (Table 1, Table 2 and Table 3). The role of CAMK and CMGC kinases dependent Ptpn6 has been shown to positively regulate TCR signaling cascades in T cells; however, the role of CAMKs in neutrophils to regulate phosphatases like Ptpn6 has not yet been described and requires further investigations [88]. Phosphatidylinositol 3, 4, 5-triphosphate 5-phosphatase 1 (Inpp5d/Ship) was identified with significant up-regulation in iIR compared to iIPC surgical group. It was assigned to cluster 3 with four regulated phosphosites that were localized in catalytic inositol polyphosphate 5 phosphatase (Inpp5d) domain, SH2 domain, putative active and catalytic sites, and putative Mg and PI/IP binding sites. Ship modulates most of the phosphatidylinositol 5’-phosphatase activity and plays a key role in neutrophils for their efficient cytoskeleton polarization and directional migration. However, spontaneous neutrophil death may result from its activation and consequently down-regulating Akt [89,90]. Ship has been well known for negatively regulating most of the immune and cellular responses like phagocytosis, degranulation, proliferation, cell migration, survival, and sensitivity to chemokines [91,92]. The phosphosite regulation in iIR in cluster 3 might negatively regulate neutrophil dependent inflammatory processes like in acute lung injury resulting in spontaneous death, chemotaxis, and directional migration after iIPC [93]. Therefore, Ship activation can be an important step in protection against iIR, found associated with iIPC. The predicted Tousled-like kinase (Tlk) (Table 3) belongs to a family of serine/threonine kinases that are conserved in plants and animals and have an important role in cell cycle progression by regulating chromatin dynamics [94]. Its proper interaction and regulation for Ship is not well known yet.

3.2. Molecular Events Associated to iIPC

Importantly, clusters 3 and 5 show abundance profiles of proteins and phosphopeptides that are profoundly affected by iIR but such effect is not observed in iIPC (Figure 3A), indicating a possible correlation to the protective effect of iIPC. The unsupervised pathways analysis searches for the rat neutrophil phosphorylated proteins with significantly regulated phosphorylation in any of the two conditions from clusters 3 and 5 were performed against KEGG database (Table 4). The most significantly enriched pathway was spliceosome with four proteins, u2af2, srsf2 (SR family), sf3b1, and usp39. In humans, about 200 spliceosomal and splicing-associated proteins regulate alternative splicing (AS) altering cytokine signaling cascades and modulating protein functions [95]. Some of the proteins assigned to spliceosome with regulated phosphosites were u2af2, sf3b1, and usp39. Changes in these proteins were found related to myeloid malignancy [96] (u2af2), clinical monoclonal B-cell lymphocytosis [97] (sf3b1), and prostate cancer [98] (usp39). The srsf2 protein in stimulated T cells was found to regulate CD45 [99]. AS can regulate protein expression in a cell-specific and tissue-specific manner in response to environmental and developmental cues; however, the exact knowledge of spliceosome dependent gene expression regulation in neutrophil is incomplete.
Both the insulin and mTOR-signaling pathways play a pivotal role in a variety of neutrophil functions. Insulin is known as an anabolic hormone that inhibits the degradation of carbohydrates, lipids, and proteins and helps in increasing their synthesis and storage [100,101]. The insulin-signaling pathway maintains glucose homeostasis. It has also been reported that insulin plays a key role in dose-dependent chemotaxis in monocytes and macrophages during wound healing [102]. Another study found an indirect confirmation that priming of neutrophil respiratory burst is a result of crosstalk of signaling pathways of insulin and fMLP receptors with the participation of tyrosine phosphorylation [103]. The identified phosphorylated proteins in insulin pathway during our analysis are Cbl, Rps6, and Inpp5d. The Cbl protein, regulated in iIR compared to Ctrl, has been known for its role in apoptosis, focal adhesion, and that its inhibition protects the heart against IR injury and minimizes the ubiquitination and subsequent degradation of EGFR and monocyte death in response to injury [104,105].
Another phosphorylated protein of the pathway, Ribosomal protein S6, regulated in iIR compared to control, Lap, and iIPC, was found common in both the insulin and mTOR pathways during our analysis. The phosphorylation of RpS6 proteins controls the cell size and also regulates the glucose homeostasis [106]. It has been found that IPC results in the phosphorylation of RpS6 protein, which interacts with components of mTOR complex 2 (mTORC2), increases mTORC2 activity, and plays a key role in cardioprotection after ischemia reperfusion [107]. The mTOR-signaling pathway plays an important role in various neutrophil functions. It has been reported that the loss of mTORC2 in neutrophils caused by silencing of Rictor results in failure of neutrophil polarization towards the chemoattractants [108]. It plays a key function in the actin cytoskeleton and also controls the chemotaxis and directional movement of the neutrophils [109,110,111].
Chemokine signaling pathway was found with significant enrichment for the four phosphorylated proteins from clusters 3 and 5 with regulated phosphosites after iIR (compared to Ctrl) and iIPC (compared to iIR). One of the assigned proteins, dock2, has been reported for its role in chemokine dependent activation in F-actin polymerization present at the leading edge of the migrating cell [112]. It also affects other cellular activities of the cell, like chemotaxis, production of superoxide, and formation of extracellular traps [113]. Paxillin, another protein of the chemokine-regulated pathway, was identified in this study with regulated phosphorylation at tyrosine and has been previously documented for its role to be involved in adhesion due to a domain interacting with the C-terminus of focal adhesion kinase (FAK) [114]. Activation of neutrophils due to tyrosine phosphorylation of Fgr domain has already been discussed, causing respiratory burst and polymerization of F-actin [115]. The rolling and recruitment of leukocytes depends on E-selectin, and it has been found to be slow and was associated with rasgrp2 [116]. Our analysis confirms the key role of phosphorylation events in dock2, Fgr, paxillin, and rasgrp2 in neutrophil activation and highlights the need of a keen understanding of the involved processes. Different pro-inflammatory mediators and cytokines are released as a result of Fc receptor interaction with Ig ligands that triggers a process of phagocytosis called antibody-dependent cell-mediated cytotoxicity [117] (ADCC). Another interesting pathway that was found significantly enriched for the phosphorylated proteins from cluster 3 and 5 is Fc-γ R-mediated phagocytosis. The enriched proteins from the two clusters, 3 and 5, to the Fc-γ R-mediated phagocytosis are dock2, Inpp5d, and marcks. Dock2 has already been known for its key role in F-actin architecture and motility and polarity during chemotaxis in neutrophil [118,119]. Marcks (myristoylated alanine-rich C kinase substrate) is a PKC target and has a role in cross-linking actin filaments which is involved in macrophage phagocytosis. Marcks deficiency resulted in the prevention of such processes [120].
Focal adhesion is an important phenomenon for the immune cells, and the proteins that affect this process are also very important in neutrophils. In our pathway analysis, the phosphorylated proteins enriched for the focal adhesion include Zyx and Pxn, regulated in iIR and iIPC vs. Ctrl and Ppp1r12a regulated in iIR and iIPC vs. Ctrl and in iIR vs. iIPC. Studies have shown acute inflammatory response as a result of high protein efflux and adhesion of PMNs in postcapillary venules [121]. Intestinal ischemia and reperfusion in mouse overexpresses P-selectin and E-selectin on neutrophils and ECs, respectively. Reduced neutrophil rolling and adhesion as result of P-selectin blockage attenuates the injury [122]. Focal adhesion-like structures have been observed in the primed neutrophils of severely injured patients that regulate various neutrophil functions [123]. A large number of adhesion associated components including various proteins are linked to actin like vinculin, alpha-actinin, talin, and zyxin among others [124]. The regulation of the phosphorylated proteins that are related to adhesion and migration reflects modulation of the major events within neutrophils after iIR and iIPC.
The SNARE protein vesicle-associated membrane protein 8 (Vamp8 - Q9WUF4) was identified in cluster 3 in its phosphorylated form. Phosphorylation has already been described as a suppressor of its function in mast cells [125]. SNARE proteins play a major role on extracellular release of neutrophil granular contents [126], which is critical for innate immune responses and the outburst of inflammation but can also cause huge organ injury during chronic or acute inflammatory conditions like that triggered by iIR [127]. Vamp8, also known as endobrevin was initially characterized as involved in fusion between early and late endosomes and can interact with syntaxin-7 (also identified as phosphorylated in cluster 3), syntaxin-8, and Vtib to form an endosomal fusion complex [128]. More recently, Vamp8 was demonstrated to be involved in regulated exocytosis by different cell types, including granulocytes, and to be localized to azurophilic granules and in low-density membrane-enriched fractions [129]. Vamp8 phosphorylation allows vesicle docking but suppresses secretion by preventing fusion with the plasma membrane [125]. Interestingly, analysis of Vamp8-deficient mice revealed a major physiological role for this protein in inflammatory conditions such as pancreatitis, as the Vamp8-/- mice showed a significantly attenuated inflammatory response compared to wild-type animals [130]. It is interesting to note that in the present study, the down-regulation of phosphorylated Vamp8 by iIR was prevented by iIPC, thus phosphorylated Vamp8 may be an important target to investigate IPC protection against IR injury. Furthermore, we found phosphorylated forms of Rab7a, Rab1a, Rab27a, and Rab31 up-regulated in iIR compared to other groups (cluster 5). The interaction between Rab proteins and SNAREs plays a pivotal role in neutrophil exocytosis [131]. As iIR seems to increase phosphorylation of Rabs, one could speculate a regulatory mechanism of Vamp8 and Rabs by phosphorylation that connects them to iIR injury and iIPC protection.

4. Materials and Methods

4.1. Chemicals and Reagents

Inhibitor Mix and Quant-iT™ Protein Assay Kit were from GE Healthcare Biosciences (Pittsburgh, USA). TEAB was from Sigma Aldrich, Switzerland. Protease PhosSTOP from Roche, Germany. Iodoacetamide was from GE Healthcare (Amersham, UK). Trypsin was from Promega (Madison, WI). DTT and SDS were from GE Healthcare Bio-Sciences (Uppsala, Sweden). iTRAQ reagents from Applied Biosystems (Foster City, CA). Vivacon 30KDa filters were from Sartorius Stedim Biotech (Goettingen, Germany). TSKGel Amide-80 2 mm, 3 mm particle size was from Tosoh Bioscience (Stuttgart, Germany). Empore C8 extraction disk was purchased from 3 M Bioanalytical Technologies (St. Paul, MN). TSKGel Amide-80 2 mm, 3 mm particle size was from Tosoh Bioscience (Stuttgart, Germany). Poros Oligo R2 and R3 reversed-phase materials from PerSeptive Biosystems (Framingham, MA). ReproSil-Pur C18 AQ 3 µm material was from Dr. Maisch (Ammerbuch-Entringen, Germany). Low-binding polypropylene microtubes were from Sorenson Bioscience (Salt Lake City, UT). TFA was from Fluka (St. Louis, MO). All other solvents and chemicals used in this work were HPLC grade or higher.

4.2. Preparation of Experimental Subjects and Sample Collection

A total of 40 male Wistar rats (weighing 250–350 g) were collected from the animal house facility, Faculty of Medicine, University of São Paulo, Brazil. The rats were randomly allocated into the four groups, i.e., control group (Ctrl, without any surgical procedure), sham laparotomy group (Lap, without clamping of any artery but receiving the same surgical procedure), intestinal ischemia and reperfusion group (iIR, subjected to 120 min of reperfusion preceded by 45 min of superior mesenteric artery occlusion (SMAO)), and intestinal ischemic preconditioning group (iIPC, subjected to 10 min of SMAO followed by 10 min of reperfusion immediately prior to 45 min of ischemia and 120 min of reperfusion as in iIR). The project was approved by the ethical committee of FMUSP (Protocol No. 8186), all the surgical procedures and sample preparation were performed in the Laboratory of Surgical Physiopathology (LIM-62), Department of Surgery, FMUSP. To ensure a high sample quality and absence of pre-existing inflammatory responses, hemocytometric analysis using granulocytes count as selection criterion was performed before the surgical procedures and sample collection as previously reported [132].

4.3. Neutrophil Isolation, Protein Digestion and iTRAQ Labelling

About 10–12 mL of blood was collected from the heart and processed for the neutrophil isolation by using Ficoll density gradient method as reported [31]. After the neutrophil isolation and counting, 3 × 106 neutrophils were tip sonicated for protein extraction in a 200 μL of lysis buffer (2% SDS, 20mM TEAB, 10 mM DTT and protease inhibitor Mix). After sonication, the protein samples were heated at 80 °C for 10 min in water bath and centrifuged, and the protein concentration was determined by using Quant-iT Protein Assay Kit (Thermo Scientific, Waltham, MA, USA). The samples were combined in pools of two to get the required protein concentration, therefore, from ten rats we got five biological replicates. To remove SDS from the protein samples, Vivacon spin filters (Sartorius Stedim, CA, USA) of 30 KDa were used according to Wisniewski et al. [133]. The DTT reduced samples were alkylated with 20 mM of Iodoacetamide in the dark, in 400 μL of 1% SDC solution at room temperature. The samples were washed with 400 μL of 1% SDC solution followed by on filter trypsin digestion in 1:50 (trypsin:protein) ratio at 37 °C overnight. After digestion, the samples were acidified with TFA to a final concentration of 0.1% and centrifuged at 15,000× g for 10 min. The peptides were purified using homemade microcolumns packed with Poros Oligo R2/R3 resins (~1 cm long) in p200 tips and the concentration was determined by amino acid analysis [134,135]. After quantification, 100 μg of purified peptides was picked into new low binding tubes from all the four conditions and replicates and dried down. The peptides were reconstituted in 20 μL of iTRAQ dissolution buffer in accordance with the manufacturer protocol. Peptides from Ctrl were labeled with 114, Lap with 115, iIR with 116, and iIPC with 117. Labeling was confirmed by MALDI and samples were multiplexed in 1:1:1:1 ratio.

4.4. Enrichment of Phosphopeptides (TiO2-SIMAC-HILIC Procedure)

Enrichment of the phosphopeptides from the iTRAQ labeled samples was performed using the following steps:
1:
First TiO2 purification
2:
IMAC purification
3:
Second TiO2 purification

4.4.1. First TiO2 Purification

The enrichment of phosphorylated peptides was performed using the titanium dioxide chromatography as previously reported [136,137,138]. Briefly, a 400 μg of the labeled peptides was reconstituted in 800 μL of loading buffer (5% TFA (v/v), 1M glycolic acid, and 80% acetonitrile (v/v) (ACN)) in low binding tubes. The samples were incubated with 0.6 mg of the TiO2 beads per 100 μg of peptides with constant high agitation at room temperature for 15 min. The samples were centrifuged (6000 rpm for 30 sec), and the supernatant was incubated with half of the amount of TiO2 used in the first incubation in a new tube. The process was repeated to recover the maximum of the phosphopeptides possible. The flow-through (FT) from the second TiO2 incubation was transferred to a new low binding tube for further downstream analysis of the total proteome. The TiO2 beads from two incubations were pooled together by using 100 μL of the loading buffer. The beads were washed with 100 μL of washing buffer-1 (1% TFA, 80% ACN) and 100 μL washing buffer-2 (0.2% TFA and 10% ACN) and dried down for 10–15 min. The bound peptides from the TiO2 surface were eluted with 100 μL of the elution buffer (60 μL ammonia solution (28%) in 940 μL Milli-Q, pH 11.3) with high agitation for 15 min at room temperature. To maximize the phosphopeptides elution and recovery, the elution step was repeated for the second time. The peptides solution was passed through a C8 stage tip to remove the TiO2 beads and the samples were lyophilized completely.

4.4.2. Enzymatic Deglycosylation

As the TiO2 chromatography also enriches the sialylated glycopeptides in parallel to the phosphorylated peptides, in order to get rid of the sialylated glycopeptides, the peptides were processed for the deglycosylation. The lyophilized peptides were re-solubilized in 50 μL of 20mM of TEAB and 2 μL of 1U/μL PNGase F, and 0.5 μL of 1U Sialidase A was added to the samples. The enzymatic reaction was performed at 37 °C over night in a wet chamber [139].

4.4.3. IMAC Purification of Multi-Phosphorylated Peptides

The identification of the multi-phosphorylated peptides is usually compromised compared to the monophosphorylated peptides due to their low abundance. To increase identification of the multi-phosphorylated peptides, IMAC purification was performed for the enrichment. The peptides solution from the deglycosylation was acidified by adding 1 μL of 10% TFA and diluted with 200 μL of 50% ACN/0.1% TFA, and the pH was adjusted to 1.6–1.8. Then 200 μL of IMAC beads were picked and washed twice with 50% ACN/0.1% TFA. The washed beads were added to the peptide solution and incubated for 30 min at room temperature. The samples were centrifuged, and half of the supernatant was transferred to a new tube whereas the remaining supernatant with the IMAC beads was shifted to the 200 μL GeLoader tip (constricted at the tip to hold the beads). The remaining supernatant was passed through with help of a syringe and the IMAC beads column was packed. The beads were washed with 100 μL of 50% ACN/0.1% TFA. The beads were washed for the second time with 70 μL of 20% ACN/1% TFA and the FT was collected with the supernatant containing monophosphorylated peptides. This step was performed slowly to allow the monophosphorylated peptides to fall off the IMAC beads and still allow the multiphosphorylated peptides to be retained. The multiphosphorylated peptides were eluted with 100 μL of ammonia (as for the first TiO2 elution) directly in a p200 stage tip with Poros Oligo R3 materials and acidified with 2 μL of 10% TFA prior to the sample touching the R3 material of the column. The multiphosphorylated peptides were R3 purified and lyophilized for analysis by mass spectrometry.

4.4.4. Second TiO2 Purification

The SIMAC FT was dried completely and re-suspended in 200 μL of 70% ACN/2% TFA solution. The peptides were incubated twice with the same amount of the TiO2 beads as in the first time. The beads were washed with 100 μL washing buffer-1 (80% ACN/1% TFA) and 100 μL of washing buffer-2 (10% ACN/0.2% TFA). All the supernatants/FTs were pooled together containing the sialylated glycopeptides (de-glycosylated). The beads were dried for 10 min, and the monophosphorylated peptides were eluted with ammonia. The peptides were passed through a C8 plug to remove the TiO2 spheres. The sample was R3 purified and HILIC fractionated.

4.4.5. HILIC Fractionation of Total and Mono-Phosphorylated Peptides

The HILIC fractionation of the purified total peptides (first TiO2 FT) and mono-phosphorylated peptides (second TiO2 eluate) was performed as previously documented [140]. Briefly, the dried peptides were re-solubilized in a buffer (90% ACN/0.1% TFA) and 40 μL of the sample was loaded onto an in-house packed TSK gel Amide-80 HILIC 320 μm × 170 mm capillary HPLC column connected to an Agilent 1200 HPLC system. The gradient used for the elution of the peptides was from 90% ACN/0.1% TFA to 60% ACN/0.1% TFA for 35 min at a flow rate of 6 μL/min. The fractions were automatically collected at 1 min interval in a microwell plate after UV detection at 210 nm. The fractions were dried completely and analyzed by the mass spectrometer.

4.5. Nano-Liquid Chromatography Tandem Mass Spectrometry (nLC-MS/MS)

The multi-phosphorylated and fractionated mono-phosphorylated peptides were analyzed by a Proxeon Easy-nLC system (Thermo Fischer Scientific, Odense, Denmark), coupled with LTQ-Orbitrap Velos mass spectrometer (Thermo Fischer Scientific). The samples were loaded onto an 18 cm long home-made reversed phase analytical column (75 μm inner diameter) packed with ReproSil-Pur C18 AQ 3 μm material (Dr. Maisch, Ammerbuch Entringen, Germany). Peptidic fractions eluted from the nLC directly into the ion source of the mass spectrometer, that was operated in data-dependent acquisition mode. The first MS scan (MS1) was carried out in the m/z range of 400–1200, at a resolution of 30,000. From each MS1 scan, the seven most intense components above a signal of 20,000 were submitted to HCD fragmentation at NCE of 36% and analyzed in the Orbitrap detector at a resolution of 7500. The spectra were saved as raw files for further analyses and were deposited at theProteomeXchange [141] consortium via the Mass Spectrometry Interactive Virtual Environment (MassIVE) platform [142] under the ID PXD016182, MSV000084545.

4.6. Database Searching and Bioinformatics Analysis

The raw files resulted from the mass spectrometer were analyzed by Proteome Discoverer version 1.4.0.288 (Thermo Fischer Scientific). The .mgf files were created from the tandem MS/MS spectra, and searches were made against the UniProt rodents’ database using Mascot (v2.3.2, Matrix Science, London, UK). The parameters used for searching of the database were: precursor mass tolerance 10 ppm, fragment (MS/MS) mass tolerance 0.05 Da, two missed cleavages were allowed, trypsin was selected as digestion enzyme, variable modifications included oxidation (M), deamidation, and for phosphopeptides serine, threonine and tyrosine were added to the searches. The carbamidomethylation of cysteine residue was selected as a fixed modification. The results were filtered for 1% false discovery rate (FDR) using Percolator as validator [143]. Further stringent filters were applied for the phosphopeptides by excluding all with phosphoRS 3.0 probability lower than 95%.
For statistical analysis, an in-house script developed under statistical package R software was used. The iTRAQ intensities from all the five biological replicates were log2 transformed and median normalized. One peptide measurement was allowed by choosing mean instead of median in RRollup function of DanteR package for the multiple measurements of the same peptide [144]. Limma and Rank Products were applied, providing sufficient power to deal with low replicate numbers and additional missing values [145,146]. Both statistical tests were carried out against label 114 (Ctrl) and corrected for multiple testing. The protein and phosphopeptides with q-value below 0.05 (5% FDR) were considered as regulated.
Cluster analysis was performed, by calculating the mean over all five replicated values for each condition. A single dataset was created after merging the phosphorylated peptides and proteins data. We applied Fuzzy c-mean clustering after determining the value of the fuzzifier and obtaining the number of clusters according to Schwämmle [147]. A standard principle component analysis (PCA) analysis was performed using R to check similarity among the replicates of the same group and variability between different conditions.
ProteinCenter (Thermo Fischer Scientific, Waltham, USA) was used to interpret the results at protein level (statistical GO slim functional categorization with 5% FDR and number of domains). WebGestalt was used with default parameters for KEGG pathways analysis [148]. Using GPS software, we performed a kinase recognition site analysis on protein sequences, and Rattus norvegicus was selected as database. The kinase-substrate prediction analysis of the significantly regulated phosphopeptides was performed in iGPS 1.0 [42] against M. musculus database using high prediction confidence, and string/experiment was selected to find the interactions. MotifX algorithm was used with windows of 17 amino acids for the analysis of sequence motif enrichment of the phosphorylation sites [149]. Phosphorylation in the regulating domains of kinases and phosphatases were checked manually in the NCBI database (www.nvbi.nlm.nih.gov/cdd).

5. Conclusions

Our deep analyses of the rat neutrophil phosphoproteome revealed many phosphorylation events that occur in various domains and motifs regions of interesting kinases and phosphatases. Most of these proteins regulate apoptosis, degranulation, and migration in the neutrophils. Our experimental design validated the role of previously studied kinases and extended the knowledge of neutrophil biology with addition of new possible players identified for the first time in our analysis. Further investigation of these new candidates with their regulating pathways will open new windows to mimic the tissues injury after iIR by the activated neutrophils. Furthermore, changes in neutrophil behavior might be due to phosphorylation events taking place in PKCδ, NDR1/Stk38, or proteins that present a phosphorylation profile compatible with modulation by iIR and prevention of such modulation by iIPC, such as Inpp5d/Ship and Vamp8. Our analysis shows that most of the important predicted kinases that might regulate other kinases identified here in neutrophils from iIR and iIPC belong to the CAMK and CMGC families. The understanding of possible interaction among different kinases and the use of documented inhibitors will help to define possible therapeutic targets to decrease the neutrophil activation and tissue injury in iIR surgical procedures.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/16/5799/s1.

Author Contributions

Conceptualization, B.F., P.R. and W.F.; Data curation, M.T., S.S., V.S. and I.S.L.; Formal analysis, M.T., M.S.C., V.S., W.F. and I.S.L.; Funding acquisition, M.S.C. and W.F.; Investigation, M.T. and P.R.; Methodology, M.T., S.A. and S.S.; Project administration, M.S.C. and W.F.; Resources, M.S.C. and P.R.; Software, M.T. and S.S.; Supervision, B.F., P.R. and W.F.; Validation, M.T. and S.A.; Writing—original draft, M.T. and S.A.; Writing—review & editing, B.F., M.S.C., W.F. and I.S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The World Academy of Sciences (TWAS), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Fundação de Apoio à Pesquisa do Distrito Federal (FAP-DF), Fundação de Empreendimentos Científicos e Tecnológicos (Finatec), and Universidade de Brasília (FUB-UnB).

Conflicts of Interest

The authors declare no conflict of interest.

Study Limitations

The present study was designed as a discovery analysis and the proteomics approach is not appropriate for target validation. The proteins and pathways highlighted here should be considered as target candidates that should be further validated by specific biochemical assays (e.g., inhibition assays) and/or by absolute quantitation (e.g., SRM assays).

Abbreviations

iIRIntestinal ischemia reperfusion
iIRIIntestinal ischemia reperfusion injury
MOFMultiple organ failure
ECsEndothelial cells
IRIschemia and reperfusion
IPCIschemic preconditioning
iIPCIntestinal ischemic preconditioning
PTMPost-translational modification
SMAOSuperior mesenteric artery occlusion
FMUSPMedical Faculty of the University of São Paulo
DTTDithiotreitol
SDCSodiu deoxycholate
iTRAQIsobaric tags for relative and absolute quantitation
SIMACSequential elution from immobilized metal affinity chromatography
HILICHydrophilic interaction liquid chromatography
ACNAcetonitrile
FTFlow-through
TFATrifluoroacetic acid
nLC-MS/MSNano-liquid chromatography tandem mass spectrometry
FDRFalse discovery rate
TiSHTiO2 SIMAC HILIC fractionation
LC-MS/MSLiquid chromatography - tandem mass spectrometry
CAMKCa2+/calmodulin-dependent protein kinase
CMGCincludes cyclin-dependent kinases (CDKs), mitogen-activated protein kinases (MAP kinases), glycogen synthase kinases (GSK) and CDK-like kinases
PMNPolymorphonuclear leukocytes
HIFHypoxia-inducible factor
LTQLinear trap quadrupole
TMTTandem mass tag
LPSLipopolysaccharides
ECEnzyme Commission
SRMSelected reaction monitoring

References

  1. Guneli, E.; Cavdar, Z.; Islekel, H.; Sarioglu, S.; Erbayraktar, S.; Kiray, M.; Sokmen, S.; Yilmaz, O.; Gokmen, N. Erythropoietin Protects the Intestine Against Ischemia/Reperfusion Injury in Rats. Mol. Med. 2007, 13, 509–517. [Google Scholar] [CrossRef] [PubMed]
  2. Mojzis, J.; Hviscová, K.; Germanova, D.; Bukovicová, D.; Mirossay, L. Protective effect of quercetin on ischemia/reperfusion-induced gastric mucosal injury in rats. Physiol. Res. 2001, 50, 501–506. [Google Scholar] [PubMed]
  3. Trompeter, M.; Brazda, T.; Remy, C.T.; Vestring, T.; Reimer, P. Non-occlusive mesenteric ischemia: Etiology, diagnosis, and interventional therapy. Eur. Radiol. 2002, 12, 1179–1187. [Google Scholar] [CrossRef] [PubMed]
  4. Bradbury, A.W.; Brittenden, J.; McBride, K.; Ruckley, C. V Mesenteric ischaemia: A multidisciplinary approach. Br. J. Surg. 1995, 82, 1446–1459. [Google Scholar] [CrossRef]
  5. Heys, S.D.; Brittenden, J.; Crofts, T.J. Acute mesenteric ischaemia: The continuing difficulty in early diagnosis. Postgrad. Med. J. 1993, 69, 48–51. [Google Scholar] [CrossRef]
  6. Souza, D.G.; Cara, D.C.; Cassali, G.D.; Coutinho, S.F.; Silveira, M.R.; Andrade, S.P.; Poole, S.P.; Teixeira, M.M. Effects of the PAF receptor antagonist UK74505 on local and remote reperfusion injuries following ischaemia of the superior mesenteric artery in the rat. Br. J. Pharmacol. 2000, 131, 1800–1808. [Google Scholar] [CrossRef] [Green Version]
  7. Ceppa, E.P.; Fuh, K.C.; Bulkley, G.B. Mesenteric hemodynamic response to circulatory shock. Curr. Opin. Crit. Care 2003, 9, 127–132. [Google Scholar] [CrossRef] [Green Version]
  8. Fontes, W.; Sousa, M.V.; Aragão, J.B.; Morhy, L. Determination of the amino acid sequence of the plant cytolysin enterolobin. Arch. Biochem. Biophys. 1997, 347, 201–207. [Google Scholar] [CrossRef]
  9. Massberg, S.; Gonzalez, A.P.; Leiderer, R.; Menger, M.D.; Messmer, K. Messmer In vivo assessment of the influence of cold preservation time on microvascular reperfusion injury after experimental small bowel transplantation. Br. J. Surg. 1998, 85, 127–133. [Google Scholar] [CrossRef]
  10. Campbell, E.L.; Kao, D.J.; Colgan, S.P. Neutrophils and the inflammatory tissue microenvironment in the mucosa. Immunol. Rev. 2016, 273, 112–120. [Google Scholar] [CrossRef]
  11. Wang, S.; Liu, H.; Wang, Q.; Cheng, Z.; Sun, S.; Zhang, Y.; Sun, X.; Wang, Z.; Ren, L. Neutrophil-to-Lymphocyte Ratio and Platelet-to-Lymphocyte Ratio Are Effective Predictors of Prognosis in Patients with Acute Mesenteric Arterial Embolism and Thrombosis. Ann. Vasc. Surg. 2018, 49, 115–122. [Google Scholar] [CrossRef] [PubMed]
  12. Sisley, A.C.; Desai, T.; Harig, J.M.; Gewertz, B.L. Neutrophil Depletion Attenuates Human Intestinal Reperfusion Injury. J. Surg. Res. 1994, 57, 192–196. [Google Scholar] [CrossRef] [PubMed]
  13. Carden, D.L.; Smith, J.K.; Korthuis, R.J. Neutrophil-mediated microvascular dysfunction in postischemic canine skeletal muscle. Role of granulocyte adherence. Circ. Res. 1990, 66, 1436–1444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Epstein, F.H.; Weiss, S.J. Tissue Destruction by Neutrophils. N. Engl. J. Med. 1989, 320, 365–376. [Google Scholar] [CrossRef]
  15. Bagge, U.L.F.; Amundson, B.; Lauritzen, C. White blood cell deformability and plugging of skeletal muscle capillaries in hemorrhagic shock. Acta Physiol. Scand. 1980, 108, 159–163. [Google Scholar] [CrossRef]
  16. Nascimento, A.; Chapeaurouge, A.; Perales, J.; Sebben, A.; Sousa, M.V.; Fontes, W.; Castro, M.S. Purification, characterization and homology analysis of ocellatin 4, a cytolytic peptide from the skin secretion of the frog Leptodactylus ocellatus. Toxicon 2007, 50, 1095–1104. [Google Scholar] [CrossRef]
  17. Willerson, J.T. Pharmacologic Approaches to Reperfusion Injury. In Advances in Pharmacology; Academic Press: Cambridge, MA, USA, 1997; pp. 291–312. [Google Scholar]
  18. Murry, C.E.; Jennings, R.B.; Reimer, K.A. Preconditioning with ischemia: A delay of lethal cell injury in ischemic myocardium. Circulation 1986, 74, 1124–1136. [Google Scholar] [CrossRef] [Green Version]
  19. Heurteaux, C.; Lauritzen, I.; Widmann, C.; Lazdunski, M. Essential role of adenosine, adenosine A1 receptors, and ATP-sensitive K+ channels in cerebral ischemic preconditioning. Proc. Natl. Acad. Sci. USA 1995, 92, 4666–4670. [Google Scholar] [CrossRef] [Green Version]
  20. Matsuyama, K.; Chiba, Y.; Ihaya, A.; Kimura, T.; Tanigawa, N.; Muraoka, R. Effect of spinal cord preconditioning on paraplegia during cross-clamping of the thoracic aorta. Ann. Thorac. Surg. 1997, 63, 1315–1320. [Google Scholar] [CrossRef]
  21. Turman, M.A.; Bates, C.M. Susceptibility of Human Proximal Tubular Cells to Hypoxia: Effect of Hypoxic Preconditioning and Comparison to Glomerular Cells. Ren. Fail. 1997, 19, 47–60. [Google Scholar] [CrossRef] [Green Version]
  22. Hotter, G.; Closa, D.; Prados, M.; Fernández-Cruz, L.; Prats, N.; Gelpí, E.; Roselló-Catafau, J. Intestinal Preconditioning Is Mediated by a Transient Increase in Nitric Oxide. Biochem. Biophys. Res. Commun. 1996, 222, 27–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Du, Z.Y.; Hicks, M.; Winlaw, D.; Spratt, P.; Macdonald, P. Ischemic preconditioning enhances donor lung preservation in the rat. J. Heart Lung Transplant. 1996, 15, 1258–1267. [Google Scholar] [PubMed]
  24. Prieto-Moure, B.; Lloris-Carsí, J.M.; Barrios-Pitarque, C.; Toledo-Pereyra, L.-H.; Lajara-Romance, J.M.; Berda-Antolí, M.; Lloris-Cejalvo, J.M.; Cejalvo-Lapeña, D. Pharmacology of Ischemia–Reperfusion. Translational Research Considerations. J. Investig. Surg. 2016, 29, 234–249. [Google Scholar] [CrossRef] [PubMed]
  25. Jenkins, D.P.; Pugsley, W.B.; Alkhulaifi, A.M.; Kemp, M.; Hooper, J.; Yellon, D.M. Ischaemic preconditioning reduces troponin T release in patients undergoing coronary artery bypass surgery. Heart 1997, 77, 314–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Clavien, P.-A.; Selzner, M.; Rüdiger, H.A.; Graf, R.; Kadry, Z.; Rousson, V.; Jochum, W. A Prospective Randomized Study in 100 Consecutive Patients Undergoing Major Liver Resection With Versus Without Ischemic Preconditioning. Ann. Surg. 2003, 238, 843–852. [Google Scholar] [CrossRef] [Green Version]
  27. Kharbanda, R.K.; Peters, M.; Walton, B.; Kattenhorn, M.; Mullen, M.; Klein, N.; Vallance, P.; Deanfield, J.; MacAllister, R. Ischemic preconditioning prevents endothelial injury and systemic neutrophil activation during ischemia-reperfusion in humans in vivo. Circulation 2001, 103, 1624–1630. [Google Scholar] [CrossRef] [Green Version]
  28. Morris, C.F.M.; Tahir, M.; Arshid, S.; Castro, M.S.; Fontes, W. Reconciling the IPC and two-hit models: Dissecting the underlying cellular and molecular mechanisms of two seemingly opposing frameworks. J. Immunol. Res. 2015, 2015. [Google Scholar] [CrossRef] [Green Version]
  29. Morris, C.F.M.; Castro, M.S.; Fontes, W. Neutrophil proteome: Lessons from different standpoints. Protein Pept. Lett. 2008, 15. [Google Scholar] [CrossRef]
  30. Tahir, M.; Arshid, S.; Fontes, B.; Castro, M.S.; Luz, I.S.; Botelho, K.L.R.; Sidoli, S.; Schwämmle, V.; Roepstorff, P.; Fontes, W. Analysis of the Effect of Intestinal Ischemia and Reperfusion on the Rat Neutrophils Proteome. Front. Mol. Biosci. 2018, 5, 89. [Google Scholar] [CrossRef] [Green Version]
  31. Arshid, S.; Tahir, M.; Fontes, B.; Montero, E.F.S.; Castro, M.S.; Sidoli, S.; Schwämmle, V.; Roepstorff, P.; Fontes, W. Neutrophil proteomic analysis reveals the participation of antioxidant enzymes, motility and ribosomal proteins in the prevention of ischemic effects by preconditioning. J. Proteom. 2017, 151, 162–173. [Google Scholar] [CrossRef]
  32. Arshid, S.; Tahir, M.; Fontes, B.; de Souza Montero, E.F.; Castro, M.S.; Sidoli, S.; Roepstorff, P.; Fontes, W. High performance mass spectrometry based proteomics reveals enzyme and signaling pathway regulation in neutrophils during the early stage of surgical trauma. Proteom. Clin. Appl. 2017, 11, 1600001. [Google Scholar] [CrossRef] [PubMed]
  33. Cohen, P. The regulation of protein function by multisite phosphorylation–a 25 year update. Trends Biochem. Sci. 2000, 25, 596–601. [Google Scholar] [CrossRef]
  34. Manning, G.; Plowman, G.D.; Hunter, T.; Sudarsanam, S. Evolution of protein kinase signaling from yeast to man. Trends Biochem. Sci. 2002, 27, 514–520. [Google Scholar] [CrossRef]
  35. Olsen, J.V.; Blagoev, B.; Gnad, F.; Macek, B.; Kumar, C.; Mortensen, P.; Mann, M. Resource Global, In Vivo, and Site-Specific Phosphorylation Dynamics in Signaling Networks. Cell 2006, 2, 635–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhang, Z.-Y. Protein Tyrosine Phosphatases: Structure and Function, Substrate Specificity, and Inhibitor Development. Annu. Rev. Pharmacol. Toxicol. 2002, 42, 209–234. [Google Scholar] [CrossRef] [PubMed]
  37. Hurley, J.H.; Dean, A.M.; Thorsness, P.E.; Koshlandjunior, D.E.; Stroud, R.M. Regulation of Isocitrate Dehydrogenase By Phosphorylation Involves no Long-Range Conformational Change in the Free enzyme. J. Biol. Chem. 1990, 265, 3599–3602. [Google Scholar]
  38. Miller, M.L.; Jensen, L.J.; Diella, F.; Jorgensen, C.; Tinti, M.; Li, L.; Hsiung, M.; Parker, S.A.; Bordeaux, J.; Sicheritz-Ponten, T.; et al. Linear Motif Atlas for Phosphorylation-Dependent Signaling. Sci. Signal. 2008, 1, ra2. [Google Scholar] [CrossRef]
  39. Campos, V.; Melo, R.C.N.; Silva, L.P.; Aquino, E.N.; Castro, M.S.; Fontes, W. Characterization of neutrophil adhesion to different titanium surfaces. Bull. Mater. Sci. 2014, 37, 157–166. [Google Scholar] [CrossRef] [Green Version]
  40. Cohen, P. The role of protein phosphorylation in human health and disease. Eur. J. Biochem. 2001, 268, 5001–5010. [Google Scholar] [CrossRef]
  41. Libério, M.S.; Joanitti, G.A.; Fontes, W.; Castro, M.S. Anticancer peptides and proteins: A panoramic view. Protein Pept. Lett. 2013, 20, 380–391. [Google Scholar]
  42. Song, C.; Ye, M.; Liu, Z.; Cheng, H.; Jiang, X.; Han, G.; Songyang, Z.; Tan, Y.; Wang, H.; Ren, J.; et al. Systematic Analysis of Protein Phosphorylation Networks From Phosphoproteomic Data. Mol. Cell. Proteom. 2012, 11, 1070–1083. [Google Scholar] [CrossRef] [Green Version]
  43. Manning, G. The Protein Kinase Complement of the Human Genome. Science 2002, 298, 1912–1934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hunter, T. Protein kinases and phosphatases: The Yin and Yang of protein phosphorylation and signaling. Cell 1995, 80, 225–236. [Google Scholar] [CrossRef] [Green Version]
  45. Schlessinger, J. Cell Signaling by Receptor Tyrosine Kinases. Cell 2000, 103, 211–225. [Google Scholar] [CrossRef] [Green Version]
  46. Macek, B.; Mijakovic, I.; Olsen, J.V.; Gnad, F.; Kumar, C.; Jensen, P.R.; Mann, M. The Serine/Threonine/Tyrosine Phosphoproteome of the Model BacteriumBacillus subtilis. Mol. Cell. Proteom. 2007, 6, 697–707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Lundby, A.; Secher, A.; Lage, K.; Nordsborg, N.B.; Dmytriyev, A.; Lundby, C.; Olsen, J. V Quantitative maps of protein phosphorylation sites across 14 different rat organs and tissues. Nat. Commun. 2012, 3, 1–10. [Google Scholar] [CrossRef] [Green Version]
  48. Karp, N.A.; Huber, W.; Sadowski, P.G.; Charles, P.D.; Hester, S.V.; Lilley, K.S. Addressing accuracy and precision issues in iTRAQ quantitation. Mol. Cell. Proteom. 2010, 9, 1885–1897. [Google Scholar] [CrossRef] [Green Version]
  49. Glibert, P.; Meert, P.; Van Steendam, K.; Van Nieuwerburgh, F.; De Coninck, D.; Martens, L.; Dhaenens, M.; Deforce, D. Phospho-iTRAQ: Assessing isobaric labels for the large-scale study of phosphopeptide stoichiometry. J. Proteome Res. 2015, 14, 839–849. [Google Scholar] [CrossRef]
  50. Palmisano, G.; Parker, B.L.; Engholm-Keller, K.; Lendal, S.E.; Kulej, K.; Schulz, M.; Schwämmle, V.; Graham, M.E.; Saxtorph, H.; Cordwell, S.J.; et al. A novel method for the simultaneous enrichment, identification, and quantification of phosphopeptides and sialylated glycopeptides applied to a temporal profile of mouse brain development. Mol. Cell. Proteom. 2012, 11, 1191–1202. [Google Scholar] [CrossRef] [Green Version]
  51. Caenepeel, S.; Charydczak, G.; Sudarsanam, S.; Hunter, T.; Manning, G. The mouse kinome: Discovery and comparative genomics of all mouse protein kinases. Proc. Natl. Acad. Sci. USA 2004, 101, 11707–11712. [Google Scholar] [CrossRef] [Green Version]
  52. Zhu, H.; Klemic, J.F.; Chang, S.; Bertone, P.; Casamayor, A.; Klemic, K.G.; Smith, D.; Gerstein, M.; Reed, M.A.; Snyder, M. Analysis of yeast protein kinases using protein chips. Nat. Genet. 2000, 26, 283–289. [Google Scholar] [CrossRef] [PubMed]
  53. Kobe, B.; Kampmann, T.; Forwood, J.K.; Listwan, P.; Brinkworth, R.I. Substrate specificity of protein kinases and computational prediction of substrates. Biochim. Biophys. Acta Proteins Proteom. 2005, 1754, 200–209. [Google Scholar] [CrossRef] [PubMed]
  54. Yan, S.R.; Fumagalli, L.; Berton, G. Activation of SRC family kinases in human neutrophils. Evidence that p58c-frg and p53/56lyn redistributed to a Triton X-100-insoluble cytoskeletal fraction, also enriched in the caveolar protein Caveolin, display an enhanced kinase activity. FEBS Lett. 1996, 380, 198–203. [Google Scholar] [PubMed] [Green Version]
  55. Zhang, H.; Meng, F.; Chu, C.-L.; Takai, T.; Lowell, C.A. The Src Family Kinases Hck and Fgr Negatively Regulate Neutrophil and Dendritic Cell Chemokine Signaling via PIR-B. Immunity 2005, 22, 235–246. [Google Scholar] [CrossRef] [Green Version]
  56. Mazzi, P.; Caveggion, E.; Lapinet-Vera, J.A.; Lowell, C.A.; Berton, G. The Src-Family Kinases Hck and Fgr Regulate Early Lipopolysaccharide-Induced Myeloid Cell Recruitment into the Lung and Their Ability To Secrete Chemokines. J. Immunol. 2015, 195, 2383–2395. [Google Scholar] [CrossRef] [Green Version]
  57. Mariño-Ramírez, L.; Hu, J.C. Isolation and mapping of self-assembling protein domains encoded by theSaccharomyces cerevisiaegenome using λ repressor fusions. Yeast 2002, 19, 641–650. [Google Scholar] [CrossRef]
  58. Wheeler, D.L.; Iida, M.; Dunn, E.F. The Role of Src in Solid Tumors. Oncologist 2009, 14, 667–678. [Google Scholar] [CrossRef]
  59. Gagné, V.; Marois, L.; Levesque, J.-M.; Galarneau, H.; Lahoud, M.H.; Caminschi, I.; Naccache, P.H.; Tessier, P.; Fernandes, M.J.G. Modulation of monosodium urate crystal-induced responses in neutrophils by the myeloid inhibitory C-type lectin-like receptor: Potential therapeutic implications. Arthritis Res. Ther. 2013, 15, R73. [Google Scholar] [CrossRef] [Green Version]
  60. Lodowski, D.T. Keeping G Proteins at Bay: A Complex between G Protein-Coupled Receptor Kinase 2 and Gbetagamma. Science 2003, 300, 1256–1262. [Google Scholar] [CrossRef]
  61. Tesmer, J.J.G.; Tesmer, V.M.; Lodowski, D.T.; Steinhagen, H.; Huber, J. Structure of Human G Protein-Coupled Receptor Kinase 2 in Complex with the Kinase Inhibitor Balanol. J. Med. Chem. 2010, 53, 1867–1870. [Google Scholar] [CrossRef] [Green Version]
  62. Ribas, C.; Penela, P.; Murga, C.; Salcedo, A.; García-Hoz, C.; Jurado-Pueyo, M.; Aymerich, I.; Mayor, F. The G protein-coupled receptor kinase (GRK) interactome: Role of GRKs in GPCR regulation and signaling. Biochim. Biophys. Acta Biomembr. 2007, 1768, 913–922. [Google Scholar] [CrossRef] [Green Version]
  63. Arraes, S.M.A.; Freitas, M.S.; da Silva, S.V.; de Paula Neto, H.A.; Alves-Filho, J.C.; Martins, M.A.; Basile-Filho, A.; Tavares-Murta, B.M.; Barja-Fidalgo, C.; Cunha, F.Q. Impaired neutrophil chemotaxis in sepsis associates with GRK expression and inhibition of actin assembly and tyrosine phosphorylation. Blood 2006, 108, 2906–2913. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Penela, P.; Ribas, C.; Mayor, F. Mechanisms of regulation of the expression and function of G protein-coupled receptor kinases. Cell. Signal. 2003, 15, 973–981. [Google Scholar] [CrossRef]
  65. Xu, X.; Jin, T. The Novel Functions of the PLC/PKC/PKD Signaling Axis in G Protein-Coupled Receptor-Mediated Chemotaxis of Neutrophils. J. Immunol. Res. 2015, 2015, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Shizukuda, Y.; Helisch, A.; Yokota, R.; Ware, J.A. Downregulation of Protein Kinase Cδ Activity Enhances Endothelial Cell Adaptation to Hypoxia. Circulation 1999, 100, 1909–1916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Soderling, T.R. The Ca2+–calmodulin-dependent protein kinase cascade. Trends Biochem. Sci. 1999, 24, 232–236. [Google Scholar] [CrossRef]
  68. Camonis, J.H.; Hergovich, A. RalA GTPase and MAP4K4 Function through NDR1 Activation in Stress Response and Apoptotic Signaling. Cell Biol. Cell Metab. 2014, 1, 1–11. [Google Scholar] [CrossRef] [Green Version]
  69. Carvalho, L.O.; Aquino, E.N.; Neves, A.C.D.; Fontes, W. The Neutrophil Nucleus and Its Role in Neutrophilic Function. J. Cell. Biochem. 2015, 116, 1831–1836. [Google Scholar] [CrossRef] [PubMed]
  70. Hergovich, A.; Stegert, M.R.; Schmitz, D.; Hemmings, B.A. NDR kinases regulate essential cell processes from yeast to humans. Nat. Rev. Mol. Cell Biol. 2006, 7, 253–264. [Google Scholar] [CrossRef] [PubMed]
  71. Zhu, Y.; Stolz, D.B.; Guo, F.; Ross, M.A.; Watkins, S.C.; Tan, B.E.E.J.E.N.; Qi, R.Z.; Manser, E.D.; Li, Q.I.U.T.; Bay, B.H.; et al. Signaling via a novel integral plasma membrane pool of a serine/threonine protein kinase PRK1 in mammalian cells. FASEB J. 2004, 18, 1722–1724. [Google Scholar] [CrossRef] [PubMed]
  72. Van Ziffle, J.A.; Lowell, C.A. Neutrophil-specific deletion of Syk kinase results in reduced host defense to bacterial infection. Blood 2009, 114, 4871–4882. [Google Scholar] [CrossRef] [Green Version]
  73. Raeder, E.M.; Mansfield, P.J.; Hinkovska-Galcheva, V.; Shayman, J.A.; Boxer, L.A. Syk activation initiates downstream signaling events during human polymorphonuclear leukocyte phagocytosis. J. Immunol. 1999, 163, 6785–6793. [Google Scholar] [PubMed]
  74. Cox, D.; Chang, P.; Kurosaki, T.; Greenberg, S. Syk Tyrosine Kinase Is Required for Immunoreceptor Tyrosine Activation Motif-dependent Actin Assembly. J. Biol. Chem. 1996, 271, 16597–16602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Luhovy, A.Y.; Jaberi, A.; Papillon, J.; Guillemette, J.; Cybulsky, A. V Regulation of the Ste20-like Kinase, SLK. J. Biol. Chem. 2011, 287, 5446–5458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Cybulsky, A.V.; Takano, T.; Papillon, J.; Khadir, A.; Bijian, K.; Chien, C.-C.; Alpers, C.E.; Rabb, H. Renal expression and activity of the germinal center kinase SK2. Am. J. Physiol. Physiol. 2004, 286, F16–F25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Cybulsky, A.V.; Takano, T.; Guillemette, J.; Papillon, J.; Volpini, R.A.; Di Battista, J.A. The Ste20-like kinase SLK promotes p53 transactivation and apoptosis. Am. J. Physiol. Renal Physiol. 2009, 297, F971–F980. [Google Scholar] [CrossRef] [PubMed]
  78. Hao, W.; Takano, T.; Guillemette, J.; Papillon, J.; Ren, G.; Cybulsky, A. V Induction of Apoptosis by the Ste20-like Kinase SLK, a Germinal Center Kinase That Activates Apoptosis Signal-regulating Kinase and p38. J. Biol. Chem. 2005, 281, 3075–3084. [Google Scholar] [CrossRef] [Green Version]
  79. Pedersen, S.F. The Na+/H+ exchanger NHE1 in stress-induced signal transduction: Implications for cell proliferation and cell death. Pflügers Arch. Eur. J. Physiol. 2006, 452, 249–259. [Google Scholar] [CrossRef]
  80. Henderson, L.M.; Chappell, J.B.; Jones, O.T.G. Internal pH changes associated with the activity of NADPH oxidase of human neutrophils. Further evidence for the presence of an H+ conducting channel. Biochem. J. 1988, 251, 563–567. [Google Scholar] [CrossRef] [Green Version]
  81. Felberg, J.; Johnson, P. Characterization of Recombinant CD45 Cytoplasmic Domain Proteins. J. Biol. Chem. 1998, 273, 17839–17845. [Google Scholar] [CrossRef] [Green Version]
  82. Streuli, M.; Krueger, N.X.; Thai, T.; Tang, M.; Saito, H. Distinct functional roles of the two intracellular phosphatase like domains of the receptor-linked protein tyrosine phosphatases LCA and LAR. EMBO J. 1990, 9, 2399–2407. [Google Scholar] [CrossRef] [PubMed]
  83. Zhu, J.W.; Doan, K.; Park, J.; Chau, A.H.; Zhang, H.; Lowell, C.A.; Weiss, A. Receptor-like Tyrosine Phosphatases CD45 and CD148 Have Distinct Functions in Chemoattractant-Mediated Neutrophil Migration and Response to S. aureus. Immunity 2011, 35, 757–769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Pei, D.; Lorenz, U.; Klingmuller, U.; Neel, B.G.; Walsh, C.T. Intramolecular Regulation of Protein Tyrosine Phosphatase SH-PTP1: A New Function for Src Homology 2 Domains. Biochemistry 1994, 33, 15483–15493. [Google Scholar] [CrossRef] [PubMed]
  85. Ravetch, J. V Immune Inhibitory Receptors. Science 2000, 290, 84–89. [Google Scholar] [CrossRef]
  86. Yousefi, S.; Simon, H.-U. SHP-1: A regulator of neutrophil apoptosis. Semin. Immunol. 2003, 15, 195–199. [Google Scholar] [CrossRef]
  87. Abu-Dayyeh, I.; Shio, M.T.; Sato, S.; Akira, S.; Cousineau, B.; Olivier, M. Leishmania-induced IRAK-1 inactivation is mediated by SHP-1 interacting with an evolutionarily conserved KTIM motif. PLoS Negl. Trop. Dis. 2008, 2, e305. [Google Scholar] [CrossRef]
  88. Chylek, L.A.; Akimov, V.; Dengjel, J.; Rigbolt, K.T.G.; Hu, B.; Hlavacek, W.S.; Blagoev, B. Phosphorylation Site Dynamics of Early T-cell Receptor Signaling. PLoS ONE 2014, 9, e104240. [Google Scholar] [CrossRef] [Green Version]
  89. Nishio, M.; Watanabe, K.; Sasaki, J.; Taya, C.; Takasuga, S.; Iizuka, R.; Balla, T.; Yamazaki, M.; Watanabe, H.; Itoh, R.; et al. Control of cell polarity and motility by the PtdIns(3,4,5)P3 phosphatase SHIP1. Nat. Cell Biol. 2006, 9, 36–44. [Google Scholar] [CrossRef]
  90. Webb, P.R.; Wang, K.Q.; Scheel-Toellner, D.; Pongracz, J.; Salmon, M.; Lord, J.M. Regulation of neutrophil apoptosis: A role for protein kinase C and phosphatidylinositol-3-kinase. Apoptosis 2000, 5, 451–458. [Google Scholar] [CrossRef]
  91. Leung, W.-H.; Tarasenko, T.; Bolland, S. Differential roles for the inositol phosphatase SHIP in the regulation of macrophages and lymphocytes. Immunol. Res. 2008, 43, 243–251. [Google Scholar] [CrossRef] [Green Version]
  92. Mondal, S.; Subramanian, K.K.; Sakai, J.; Bajrami, B.; Luo, H.R. Phosphoinositide lipid phosphatase SHIP1 and PTEN coordinate to regulate cell migration and adhesion. Mol. Biol. Cell 2012, 23, 1219–1230. [Google Scholar] [CrossRef] [PubMed]
  93. Strassheim, D.; Kim, J.-Y.; Park, J.-S.; Mitra, S.; Abraham, E. Involvement of SHIP in TLR2-Induced Neutrophil Activation and Acute Lung Injury. J. Immunol. 2005, 174, 8064–8071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Carrera, P. Tousled-like kinase functions with the chromatin assembly pathway regulating nuclear divisions. Genes Dev. 2003, 17, 2578–2590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Sahoo, A.; Im, S.-H. Interleukin and Interleukin Receptor Diversity: Role of Alternative Splicing. Int. Rev. Immunol. 2010, 29, 77–109. [Google Scholar] [CrossRef]
  96. Graubert, T.A.; Shen, D.; Ding, L.; Okeyo-Owuor, T.; Lunn, C.L.; Shao, J.; Krysiak, K.; Harris, C.C.; Koboldt, D.C.; Larson, D.E.; et al. Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes. Nat. Genet. 2011, 44, 53–57. [Google Scholar] [CrossRef] [Green Version]
  97. Quesada, V.; Conde, L.; Villamor, N.; Ordóñez, G.R.; Jares, P.; Bassaganyas, L.; Ramsay, A.J.; Beà, S.; Pinyol, M.; Martínez-Trillos, A.; et al. Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia. Nat. Genet. 2011, 44, 47–52. [Google Scholar] [CrossRef]
  98. Wen, D.; Xu, Z.; Xia, L.; Liu, X.; Tu, Y.; Lei, H.; Wang, W.; Wang, T.; Song, L.; Ma, C.; et al. Important Role of SUMOylation of Spliceosome Factors in Prostate Cancer Cells. J. Proteome Res. 2014, 13, 3571–3582. [Google Scholar] [CrossRef]
  99. Tong, A.; Nguyen, J.; Lynch, K.W. Differential Expression of CD45 Isoforms Is Controlled by the Combined Activity of Basal and Inducible Splicing-regulatory Elements in Each of the Variable Exons. J. Biol. Chem. 2005, 280, 38297–38304. [Google Scholar] [CrossRef] [Green Version]
  100. Krogh-Madsen, R.; Plomgaard, P.; Keller, P.; Keller, C.; Pedersen, B.K. Insulin stimulates interleukin-6 and tumor necrosis factor-alpha gene expression in human subcutaneous adipose tissue. Am. J. Physiol. Endocrinol. Metab. 2004, 286, E234–E238. [Google Scholar] [CrossRef]
  101. Iida, K.T.; Shimano, H.; Kawakami, Y.; Sone, H.; Toyoshima, H.; Suzuki, S.; Asano, T.; Okuda, Y.; Yamada, N. Insulin up-regulates tumor necrosis factor-alpha production in macrophages through an extracellular-regulated kinase-dependent pathway. J. Biol. Chem. 2001, 276, 32531–32537. [Google Scholar] [CrossRef] [Green Version]
  102. Liu, Y.; Dhall, S.; Castro, A.; Chan, A.; Alamat, R.; Martins-Green, M. Insulin regulates multiple signaling pathways leading to monocyte/macrophage chemotaxis into the wound tissue. Biol. Open 2018, 7, 026187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Safronova, V.G.; Gabdoulkhakova, A.G.; Miller, A.V.; Kosarev, I.V.; Vasilenko, R.N. Variations of the effect of insulin on neutrophil respiratory burst. The role of Tyrosine kinases and phosphatases. Biochem. 2001, 66, 840–849. [Google Scholar]
  104. Rafiq, K.; Kolpakov, M.A.; Seqqat, R.; Guo, J.; Guo, X.; Qi, Z.; Yu, D.; Mohapatra, B.; Zutshi, N.; An, W.; et al. c-Cbl inhibition improves cardiac function and survival in response to myocardial ischemia. Circulation 2014, 129, 2031–2043. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Rafiq, K.; Guo, J.; Vlasenko, L.; Guo, X.; Kolpakov, M.A.; Sanjay, A.; Houser, S.R.; Sabri, A. c-Cbl ubiquitin ligase regulates focal adhesion protein turnover and myofibril degeneration induced by neutrophil protease cathepsin G. J. Biol. Chem. 2012, 287, 5327–5339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Ruvinsky, I.; Sharon, N.; Lerer, T.; Cohen, H.; Stolovich-Rain, M.; Nir, T.; Dor, Y.; Zisman, P.; Meyuhas, O. Ribosomal protein S6 phosphorylation is a determinant of cell size and glucose homeostasis. Genes Dev. 2005, 19, 2199–2211. [Google Scholar] [CrossRef] [Green Version]
  107. Yano, T.; Ferlito, M.; Aponte, A.; Kuno, A.; Miura, T.; Murphy, E.; Steenbergen, C. Pivotal role of mTORC2 and involvement of ribosomal protein S6 in cardioprotective signaling. Circ. Res. 2014, 114, 1268–1280. [Google Scholar] [CrossRef] [Green Version]
  108. Liu, L.; Das, S.; Losert, W.; Parent, C.A. MTORC2 Regulates Neutrophil Chemotaxis in a cAMP- and RhoA-Dependent Fashion. Dev. Cell 2010, 19, 845–857. [Google Scholar] [CrossRef] [Green Version]
  109. Jacinto, E.; Loewith, R.; Schmidt, A.; Lin, S.; Rüegg, M.A.; Hall, A.; Hall, M.N. Mammalian TOR complex 2 controls the actin cytoskeleton and is rapamycin insensitive. Nat. Cell Biol. 2004, 6, 1122–1128. [Google Scholar] [CrossRef]
  110. He, Y.; Li, D.; Cook, S.L.; Yoon, M.-S.; Kapoor, A.; Rao, C.V.; Kenis, P.J.A.; Chen, J.; Wang, F. Mammalian target of rapamycin and Rictor control neutrophil chemotaxis by regulating Rac/Cdc42 activity and the actin cytoskeleton. Mol. Biol. Cell 2013, 24, 3369–3380. [Google Scholar] [CrossRef] [Green Version]
  111. Bao, Y.; Ledderose, C.; Graf, A.F.; Brix, B.; Birsak, T.; Lee, A.; Zhang, J.; Junger, W.G. mTOR and differential activation of mitochondria orchestrate neutrophil chemotaxis. J. Cell Biol. 2015, 210, 1153–1164. [Google Scholar] [CrossRef]
  112. Nishikimi, A.; Fukuhara, H.; Su, W.; Hongu, T.; Takasuga, S.; Mihara, H.; Cao, Q.; Sanematsu, F.; Kanai, M.; Hasegawa, H.; et al. Sequential Regulation of DOCK2 Dynamics by Two Phospholipids During Neutrophil Chemotaxis. Science 2009, 324, 384–387. [Google Scholar] [CrossRef] [Green Version]
  113. Watanabe, M.; Terasawa, M.; Miyano, K.; Yanagihara, T.; Uruno, T.; Sanematsu, F.; Nishikimi, A.; Côté, J.-F.; Sumimoto, H.; Fukui, Y. DOCK2 and DOCK5 Act Additively in Neutrophils To Regulate Chemotaxis, Superoxide Production, and Extracellular Trap Formation. J. Immunol. 2014, 193, 5660–5667. [Google Scholar] [CrossRef] [Green Version]
  114. Tachibana, K. Direct association of pp125FAK with paxillin, the focal adhesion- targeting mechanism of pp125FAK. J. Exp. Med. 1995, 182, 1089–1099. [Google Scholar] [CrossRef]
  115. Fumagalli, L.; Zhang, H.; Baruzzi, A.; Lowell, C.A.; Berton, G. The Src Family Kinases Hck and Fgr Regulate Neutrophil Responses to N-Formyl-Methionyl-Leucyl-Phenylalanine. J. Immunol. 2007, 178, 3874–3885. [Google Scholar] [CrossRef]
  116. Stadtmann, A.; Brinkhaus, L.; Mueller, H.; Rossaint, J.; Bolomini-Vittori, M.; Bergmeier, W.; Van Aken, H.; Wagner, D.D.; Laudanna, C.; Ley, K.; et al. Rap1a activation by CalDAG-GEFI and p38 MAPK is involved in E-selectin-dependent slow leukocyte rolling. Eur. J. Immunol. 2011, 41, 2074–2085. [Google Scholar] [CrossRef] [Green Version]
  117. Sánchez-Mejorada, G.; Rosales, C. Signal transduction by immunoglobulin Fc receptors. J. Leukoc. Biol. 1998, 63, 521–533. [Google Scholar] [CrossRef]
  118. Floc’h, A.L.; Tanaka, Y.; Bantilan, N.S.; Voisinne, G.; Altan-Bonnet, G.; Fukui, Y.; Huse, M. Annular PIP3 accumulation controls actin architecture and modulates cytotoxicity at the immunological synapse. J. Exp. Med. 2013, 210, 2721–2737. [Google Scholar] [CrossRef] [Green Version]
  119. Kunisaki, Y.; Nishikimi, A.; Tanaka, Y.; Takii, R.; Noda, M.; Inayoshi, A.; Watanabe, K.; Sanematsu, F.; Sasazuki, T.; Sasaki, T.; et al. DOCK2 is a Rac activator that regulates motility and polarity during neutrophil chemotaxis. J. Cell Biol. 2006, 174, 647–652. [Google Scholar] [CrossRef]
  120. Carballo, E.; Pitterle, D.M.; Stumpo, D.J.; Sperling, R.T.; Blackshear, P.J. Phagocytic and macropinocytic activity in MARCKS-deficient macrophages and fibroblasts. Am. J. Physiol. Physiol. 1999, 277, C163–C173. [Google Scholar] [CrossRef]
  121. Granger, D.N.; Korthuis, R.J. Physiologic Mechanisms of Postischemic Tissue Injury. Annu. Rev. Physiol. 1995, 57, 311–332. [Google Scholar] [CrossRef]
  122. Riaz, A.A.; Wan, M.X.; Schaefer, T.; Schramm, R.; Ekberg, H.; Menger, M.D.; Jeppsson, B.; Thorlacius, H. Fundamental and Distinct Roles of P-Selectin and LFA-1 in Ischemia/Reperfusion-Induced Leukocyte-Endothelium Interactions in the Mouse Colon. Ann. Surg. 2002, 236, 777–784. [Google Scholar] [CrossRef]
  123. Aquino, E.N.; Neves, A.C.; Santos, K.C.; Uribe, C.E.; De Souza, P.E.; Correa, J.R.; Castro, M.S.; Fontes, W. Proteomic Analysis of Neutrophil Priming by PAF. Protein Pept. Lett. 2015, 23, 142–151. [Google Scholar] [CrossRef]
  124. Yoshigi, M.; Hoffman, L.M.; Jensen, C.C.; Yost, H.J.; Beckerle, M.C. Mechanical force mobilizes zyxin from focal adhesions to actin filaments and regulates cytoskeletal reinforcement. J. Cell Biol. 2005, 171, 209–215. [Google Scholar] [CrossRef] [Green Version]
  125. Malmersjö, S.; Palma, S.D.; Diao, J.; Lai, Y.; Pfuetzner, R.A.; Wang, A.L.; Mcmahon, M.A.; Hayer, A.; Porteus, M.; Bodenmiller, B.; et al. Phosphorylation of residues inside the SNARE complex suppresses secretory vesicle fusion. EMBO J. 2016, 35, 1810–1821. [Google Scholar] [CrossRef] [PubMed]
  126. Mollinedo, F.; Calafat, J.; Janssen, H.; Martín-Martín, B.; Canchado, J.; Nabokina, S.M.; Gajate, C. Combinatorial SNARE Complexes Modulate the Secretion of Cytoplasmic Granules in Human Neutrophils. J. Immunol. 2006, 177, 2831–2841. [Google Scholar] [CrossRef] [Green Version]
  127. Tapper, H. The secretion of preformed granules by macrophages and neutrophils. J. Leukoc. Biol. 1996, 59, 613–622. [Google Scholar] [CrossRef]
  128. Antonin, W. A SNARE complex mediating fusion of late endosomes defines conserved properties of SNARE structure and function. EMBO J. 2000, 19, 6453–6464. [Google Scholar] [CrossRef]
  129. Logan, M.R.; Lacy, P.; Odemuyiwa, S.O.; Steward, M.; Davoine, F.; Kita, H.; Moqbel, R. A critical role for vesicle-associated membrane protein-7 in exocytosis from human eosinophils and neutrophils. Allergy 2006, 61, 777–784. [Google Scholar] [CrossRef]
  130. Wang, C.C.; Ng, C.P.; Lu, L.; Atlashkin, V.; Zhang, W.; Seet, L.F.; Hong, W. A role of VAMP8/Endobrevin in regulated exocytosis of pancreatic acinar cells. Dev. Cell 2004, 7, 359–371. [Google Scholar] [CrossRef] [Green Version]
  131. Lacy, P. The role of Rho GTPases and SNAREs in mediator release from granulocytes. Pharmacol. Ther. 2005, 107, 358–376. [Google Scholar] [CrossRef]
  132. Tahir, M.; Arshid, S.; Heimbecker, A.M.C.; Castro, M.S.; de Souza Montero, E.F.; Fontes, B.; Fontes, W. Evaluation of the effects of ischemic preconditioning on the hematological parameters of rats subjected to intestinal ischemia and reperfusion. Clinics 2015, 70, 61–68. [Google Scholar] [CrossRef]
  133. Wiśniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef] [PubMed]
  134. Gobom, J.; Nordhoff, E.; Mirgorodskaya, E.; Ekman, R.; Roepstorff, P. Sample purification and preparation technique based on nano-scale reversed-phase columns for the sensitive analysis of complex peptide mixtures by matrix-assisted laser desorption/ionization mass spectrometry. J. Mass Spectrom. 1999, 34, 105–116. [Google Scholar] [CrossRef]
  135. Laursen, I.; Højrup, P.; Houen, G.; Christiansen, M. Characterisation of the 1st SSI purified MBL standard. Clin. Chim. Acta 2008, 395, 159–161. [Google Scholar] [CrossRef]
  136. Engholm-Keller, K.; Hansen, T.A.; Palmisano, G.; Larsen, M.R. Multidimensional Strategy for Sensitive Phosphoproteomics Incorporating Protein Prefractionation Combined with SIMAC, HILIC, and TiO2Chromatography Applied to Proximal EGF Signaling. J. Proteome Res. 2011, 10, 5383–5397. [Google Scholar] [CrossRef] [PubMed]
  137. Thingholm, T.E.; Jørgensen, T.J.D.; Jensen, O.N.; Larsen, M.R. Highly selective enrichment of phosphorylated peptides using titanium dioxide. Nat. Protoc. 2006, 1, 1929–1935. [Google Scholar] [CrossRef]
  138. Larsen, M.R.; Thingholm, T.E.; Jensen, O.N.; Roepstorff, P.; Jørgensen, T.J.D. Highly selective enrichment of phosphorylated peptides from peptide mixtures using titanium dioxide microcolumns. Mol. Cell. Proteom. 2005, 4, 873–886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Larsen, M.R.; Jensen, S.S.; Jakobsen, L.A.; Heegaard, N.H.H. Exploring the Sialiome Using Titanium Dioxide Chromatography and Mass Spectrometry. Mol. Cell. Proteom. 2007, 6, 1778–1787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  140. McNulty, D.E.; Annan, R.S. Hydrophilic Interaction Chromatography Reduces the Complexity of the Phosphoproteome and Improves Global Phosphopeptide Isolation and Detection. Mol. Cell. Proteom. 2008, 7, 971–980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  141. Perez-Riverol, Y.; Alpi, E.; Wang, R.; Hermjakob, H.; Vizcaino, J.A. Making proteomics data accessible and reusable: Current state of proteomics databases and repositories. Proteomics 2015, 15, 930–949. [Google Scholar] [CrossRef] [Green Version]
  142. Jarnuczak, A.F.; Vizcaino, J.A. Using the PRIDE Database and ProteomeXchange for Submitting and Accessing Public Proteomics Datasets. Curr. Protoc. Bioinform. 2017, 59, 13–31. [Google Scholar] [CrossRef] [PubMed]
  143. Spivak, M.; Weston, J.; Bottou, L.; Käll, L.; Noble, W.S. Improvements to the Percolator Algorithm for Peptide Identification from Shotgun Proteomics Data Sets. J. Proteome Res. 2009, 8, 3737–3745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  144. Taverner, T.; Karpievitch, Y.V.; Polpitiya, A.D.; Brown, J.N.; Dabney, A.R.; Anderson, G.A.; Smith, R.D. DanteR: An extensible R-based tool for quantitative analysis of -omics data. Bioinformatics 2012, 28, 2404–2406. [Google Scholar] [CrossRef] [PubMed]
  145. Breitling, R.; Armengaud, P.; Amtmann, A.; Herzyk, P. Rank products: A simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573, 83–92. [Google Scholar] [CrossRef]
  146. Schwämmle, V.; León, I.R.; Jensen, O.N. Assessment and Improvement of Statistical Tools for Comparative Proteomics Analysis of Sparse Data Sets with Few Experimental Replicates. J. Proteome Res. 2013, 12, 3874–3883. [Google Scholar] [CrossRef]
  147. Schwämmle, V.; Jensen, O.N. A simple and fast method to determine the parameters for fuzzy c–means cluster analysis. Bioinformatics 2010, 26, 2841–2848. [Google Scholar] [CrossRef]
  148. Wang, J.; Duncan, D.; Shi, Z.; Zhang, B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): Update 2013. Nucleic Acids Res. 2013, 41, W77–W83. [Google Scholar] [CrossRef] [Green Version]
  149. Chou, M.F.; Schwartz, D. Biological Sequence Motif Discovery Usingmotif-x. In Current Protocols in Bioinformatics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2011; ISBN 0471250953. [Google Scholar]
Figure 1. Schematic representation of the experimental workflow. Isolation of neutrophils by Ficoll density gradient from control (Ctrl), laparotomy (Lap), intestinal ischemia/reperfusion (iIR), and intestinal ischemic preconditioning (iIPC) groups followed by protein extraction and digestion. The iTRAQ labeling 114 for Ctrl, 115 for Lap, 116 for iIR, and 117 for iIPC was used. Mono- and multi-phosphopeptides enrichment was performed using TiSH (TiO2 SIMAC HILIC fractionation). LTQ Orbitrap mass spectrometry (MS) analysis followed bioinformatics tools.
Figure 1. Schematic representation of the experimental workflow. Isolation of neutrophils by Ficoll density gradient from control (Ctrl), laparotomy (Lap), intestinal ischemia/reperfusion (iIR), and intestinal ischemic preconditioning (iIPC) groups followed by protein extraction and digestion. The iTRAQ labeling 114 for Ctrl, 115 for Lap, 116 for iIR, and 117 for iIPC was used. Mono- and multi-phosphopeptides enrichment was performed using TiSH (TiO2 SIMAC HILIC fractionation). LTQ Orbitrap mass spectrometry (MS) analysis followed bioinformatics tools.
Ijms 21 05799 g001
Figure 2. Statistical overview of total and phosphoproteome. (A) Overlap between proteome and phosphoproteome. (B) Overlap between phosphorylated serine, threonine, and tyrosine containing phosphoproteome. (C) Overlap between serine, threonine, and tyrosine phosphorylated peptides. (D) Phosphorylated peptides distribution. (E) Phosphosites distribution. (F) Phosphorylated proteins distribution.
Figure 2. Statistical overview of total and phosphoproteome. (A) Overlap between proteome and phosphoproteome. (B) Overlap between phosphorylated serine, threonine, and tyrosine containing phosphoproteome. (C) Overlap between serine, threonine, and tyrosine phosphorylated peptides. (D) Phosphorylated peptides distribution. (E) Phosphosites distribution. (F) Phosphorylated proteins distribution.
Ijms 21 05799 g002
Figure 3. (A) Cluster-based relative abundance profiles of the identified proteins and phosphopeptides in Ctrl, Lap, iIR, and iIPC from rat neutrophil. (B) Phosphosites, phosphopeptides, and phosphoproteins distribution across the six clusters.
Figure 3. (A) Cluster-based relative abundance profiles of the identified proteins and phosphopeptides in Ctrl, Lap, iIR, and iIPC from rat neutrophil. (B) Phosphosites, phosphopeptides, and phosphoproteins distribution across the six clusters.
Ijms 21 05799 g003
Figure 4. Analysis of phosphorylation motifs grouped by regulatory trends. (A) Venn diagram of phosphorylated proteins with enriched serine and threonine motifs. (B) Distribution of the enriched serine and threonine motifs abundance across six clusters.
Figure 4. Analysis of phosphorylation motifs grouped by regulatory trends. (A) Venn diagram of phosphorylated proteins with enriched serine and threonine motifs. (B) Distribution of the enriched serine and threonine motifs abundance across six clusters.
Ijms 21 05799 g004
Figure 5. Distribution of the predicted kinase families responsible for the phosphorylation event in the serine (A) and threonine (B) motifs.
Figure 5. Distribution of the predicted kinase families responsible for the phosphorylation event in the serine (A) and threonine (B) motifs.
Ijms 21 05799 g005
Table 1. Kinases and Phosphatases with regulated phosphorylation.
Table 1. Kinases and Phosphatases with regulated phosphorylation.
Protein IDDescriptionGene SymbolRegulated Phosphopeptide in DomainPhospho Reg. ClusterSig. RegulationPhosphosite in ProteinEnzyme CodeDomain Containing the PhosphopeptideDomain Description
Kinases
Q6P6U0Tyrosine-protein kinase FgrFgrLIVDDEYphNPQQGTKFPIK2IR vs. Ctrl Y400EC: 2.7.10.2PTKc_Src_Fyn_likeCatalytic domain of a subset of Src kinase-like Protein Tyrosine Kinases
Pkinase_TyrProtein tyrosine kinase
Q64303Serine/threonine-protein kinase PAK 2Pak2--------EC: 2.7.11.1----
E9PTG8Serine/threonine-protein kinase 10Stk10--------EC: 2.7.11.1----
Q63531Ribosomal protein S6 kinase alpha-1Rps6ka1--------EC: 2.7.11.1----
Q9WUT3Ribosomal protein S6 kinase alpha-2Rps6ka2--------EC: 2.7.11.1----
P26817Beta-adrenergic receptor kinase 1Adrbk1NKPRSphPVVELSK2IR vs. Ctrl S670EC: 2.7.11.15PH_GRK2_subgroupG Protein-Coupled Receptor Kinase 2 subgroup pleckstrin homology (PH) domain
G-beta gamma binding site
Q63433Serine/threonine-protein kinase N1Pkn1--------EC: 2.7.11.13----
P09215Protein kinase C delta typePrkcdSPSDYSNFDPEFLNEKPQLSphFSDK6IR vs. Ctrl S643EC: 2.7.11.13STKc_nPKC_deltaCatalytic domain of the Serine/Threonine Kinase, Novel Protein Kinase C delta
P09215Protein kinase C delta typePrkcdSPSDYSNFDPEFLNEKPQLSFSphDK6IR vs. Ctrl S645EC: 2.7.11.13STKc_nPKC_deltaCatalytic domain of the Serine/Threonine Kinase, Novel Protein Kinase C delta
Q64725Tyrosine-protein kinase SYKSyk--------EC: 2.7.10.2----
O08815STE20-like serine/threonine-protein kinaseSlk--------EC: 2.7.11.1----
Q6P9R2Serine/threonine-protein kinase OSR1Oxsr1--------EC: 2.7.11.1----
Q91VJ4Serine/threonine-protein kinase 38Stk38FEGLTphAR1IR vs. IPC T452EC: 2.7.11.1STKc_NDR1Catalytic domain of the Serine/Threonine Kinase, Nuclear Dbf2-Related kinase 1
Phosphatases
P81718Tyrosine-protein phosphatase non-receptor type 6Ptpn6DLSphGPDAETLLK2IR vs. Ctrl S12EC: 3.1.3.48SH2_N-SH2_SHP_likeN-terminal Src homology 2 (N-SH2) domain found in SH2 domain Phosphatases (SHP) proteins
P97573Phosphatidylinositol 3,4,5-trisphosphate 5-phosphatase 1Inpp5dTGIANTphLGNK3IR vs. IPC T519EC: 3.1.3.86INPP5D_SHIP1-INPP5DCatalytic inositol polyphosphate 5-phosphatase (INPP5d) domain of SH2 domain
Putative catalytic site
Putative active site
Putative Mg binding site
Putative PI/IP binding site
B2GV87Receptor-type tyrosine-protein phosphatase epsilon precursorPtpre--------EC: 3.1.3.48----
P04157Receptor-type tyrosine-protein phosphatase C isoform 4 precursorPtprc--------EC: 3.1.3.48----
Table 2. Enriched motifs found in the regulated phosphoproteins.
Table 2. Enriched motifs found in the regulated phosphoproteins.
Class/IDSignificantly Regulated PhosphopeptideSignificant RegulationModified pos. in PeptideData Source 1Position in ProteinClusterMotif
Kinases
Q6P6U0EDVGLEGDFRSphQGAEERIR vs. CtrlS11RATS255
LIVDDEYphNPQQGTKFPIKIR vs. CtrlY7RATY4002
SphSSISPQPISPAFLNVGNIRIR vs. CtrlS1RATS412S.S
SphSSISPQPISphPAFLNVGNIRIR vs. CtrlS1RATS415S.S
SSphSISPQPISphPAFLNVGNIRIR vs. CtrlS2RATS425S.......S
SSSISphPQPISphPAFLNVGNIRIR vs. CtrlS5RATS455SP
SSSphISphPQPISPAFLNVGNIRIR vs. CtrlS3RATS435S.S
SSSphISPQPISphPAFLNVGNIRIR vs. CtrlS3RATS435S.S
EDVGLEGDFRSphQGAEERIR vs. LapS11RATS255
Q64303FYDSphNTVKIR vs. CtrlS4RATS1322
YLSphFTPPEKIR vs. LapS3RATS1416
E9PTG8ILRLSphTFEKIR vs. CtrlS5RATS136
LSphTFEKIR vs. CtrlS2RATS136
LSphTFEKIR vs. LapS2RATS136
LSTphFEKIR vs. IPCT3RATT141
Q63531KLPSphTTLIR vs. CtrlS4RATS7322
Q9WUT3LEPVLSphSSLAQRIR vs. CtrlS6MOUSES7163S.S
LEPVLSphSSLAQRIR vs. IPCS6MOUSES7163S.S
P26817NKPRSphPVVELSKIR vs. CtrlS5RATS6702SP
Q63433SGSphLSGRIR vs. CtrlS3RATS3772S.S
P09215SPSDYSNFDPEFLNEKPQLSFSphDKIR vs. CtrlS22RATS6456S.S
SPSDYSNFDPEFLNEKPQLSphFSDKIR vs. CtrlS20RATS6436S.S
SPSDYSNFDPEFLNEKPQLSphFSDKIR vs. LapS20RATS6436S.S
SPSDYSNFDPEFLNEKPQLSphFSDKIR vs. IPCS20RATS6436S.S
Q64725SYSphFPKPGHKIR vs. CtrlS3RATS2916S.S
O08815TKDSGSphVSLQETRIR vs. LapS6RATS7783S.S
Q6P9R2AAISQLRSphPRIR vs. IPCS8MOUSES3595SP
RVPGSphSphGRLHKIR vs. IPCS5MOUSES3243
Q91VJ4FEGLTphARIR vs. IPCT5MOUSET4521
Phosphatases
P81718DLSphGPDAETLLKIR vs. CtrlS3RATS122R..S
TphSSKHKEEVYENVHSKIR vs. CtrlT1RATT5575
DLSphGPDAETLLKIR vs. LapS3RATS122R..S
P97573DSSLGPGRGEGPPTphPPSQPPLSPKIR vs. CtrlT14RATT9632TP
GEGPPTphPPSQPPLSphPKIR vs. CtrlT6RATT9635TP
GEGPPTphPPSQPPLSphPKKIR vs. CtrlT6RATT9635TP
KEQESphPKIR vs. CtrlS5RATS10372SP
KEQESphPKIR vs. LapS5RATS10372SP
GEGPPTphPPSQPPLSPKIR vs. IPCT6RATT9632TP
TGIANTphLGNKIR vs. IPCT6RATT5193
B2GV87SPSphGPKKIR vs. CtrlS3RATS1063S.S
SPSphGPKKIR vs. LapS3RATS1063S.S
P04157ANSphQDKIEFHNEVDGAKIR vs. LapS3RATS12096
KANSphQDKIR vs. LapS4RATS12096
KANSphQDKIEFHNEVDGAKIR vs. LapS4RATS12096
ANSphQDKIEFHNEVDGAKIR vs. IPCS3RATS12096
KANSphQDKIR vs. IPCS4RATS12096
KANSphQDKIEFHNEVDGAKIR vs. IPCS4RATS12096
Table 3. Enriched motifs found in the regulated phosphoproteins.
Table 3. Enriched motifs found in the regulated phosphoproteins.
Acc. No.PhosphositePredicted KinasePhosphorylated PeptideScoreCutoffMotifCluster
Kinases
Q6P6U0Y400TK/SrcLIVDDEYphNPQQGTKFPIK24.6451.63--6
Y400TK/Tec*LIVDDEYphNPQQGTKFPIK23.3413.584--
Y400TK/JakLIVDDEYphNPQQGTKFPIK18.2428.154--
Y400TK/FAK*LIVDDEYphNPQQGTKFPIK14.8085.968--
Y400TK/DDR*LIVDDEYphNPQQGTKFPIK112.683--
Y400TK/SykLIVDDEYphNPQQGTKFPIK8.8252.436--
Y400TK/Csk*LIVDDEYphNPQQGTKFPIK6.7784.886--
Y400TK/Abl*LIVDDEYphNPQQGTKFPIK5.4224.747--
Y400TK/Met*LIVDDEYphNPQQGTKFPIK4.6072.379--
Y400TK/VEGFR*LIVDDEYphNPQQGTKFPIK4.5563.654--
Y400TK/AlkLIVDDEYphNPQQGTKFPIK3.6673.333--
Y400TK/PDGFR*LIVDDEYphNPQQGTKFPIK3.6462.352--
Y400TK*LIVDDEYphNPQQGTKFPIK2.5160.166--
P26817S670CMGC/MAPKNKPRSphPVVELSK42.73535.046SP4
S670AGC/PDK1NKPRSphPVVELSK5.2222.257SP
S670CMGC*NKPRSphPVVELSK1.3560.963SP
S670CMGC/DYRK*NKPRSphPVVELSK1.3331.276SP
S670AGC/PKCNKPRSphPVVELSK0.3840.236SP
P09215S643CAMK/RAD53*SPSDYSNFDPEFLNEKPQLSphFSDK13.0757.385S.S5
S643AGC/PKCSPSDYSNFDPEFLNEKPQLSphFSDK0.9050.236S.S
S643Atypical/PDHK*SPSDYSNFDPEFLNEKPQLSphFSDK4.4053.075S.S
S645CAMK/CAMK1*SPSDYSNFDPEFLNEKPQLSFSphDK3.2592.488S.S
S645TKL/STKR*SPSDYSNFDPEFLNEKPQLSFSphDK2.9382.562S.S
Q91VJ4T452TKL*FEGLTphAR4.6484.354--4
T452AGC/PKCFEGLTphAR0.6080.236--
Phosphatases
P81718S12CAMK/PHK*DLSphGPDAETLLK22.2699.527R..S4
S12CMGC/CK2*DLSphGPDAETLLK12.4679.894R..S
S12CMGC/CLK*DLSphGPDAETLLK5.3754.3R..S
P97573T519Other/TLK*TGIANTphLGNK6.255.775--5
T519Other/TTK*TGIANTphLGNK5.1885.009--
Predicted Kinases with * shows the absence of inhibitors in the database.
Table 4. Overrepresented KEGG pathways associated to regulated phosphoproteins.
Table 4. Overrepresented KEGG pathways associated to regulated phosphoproteins.
KEGG PathwayAccession No.Gene SymbolDescriptionNo. of Phosphopeptides
Spliceosome (C = 135; O = 4; E = 0.26; R = 15.49; rawP = 0.0001; adjP = 0.0017)
P26369U2af2U2 small nuclear ribonucleoprotein auxiliary factor (U2AF) 21
Q62093Srsf2serine/arginine-rich splicing factor 22
Q99NB9Sf3b1splicing factor 3b, subunit 12
Q3TIX9Usp39ubiquitin specific peptidase 391
Chemokine signaling pathway(C = 178; O = 4; E = 0.34; R = 11.75; rawP = 0.0004; adjP = 0.0034)
Q8C3J5Dock2dedicator of cytokinesis 24
Q66H76Pxnpaxillin2
Q6P6U0FgrGardner–Rasheed feline sarcoma viral (v-fgr) oncogene homolog6
P0C643Rasgrp2RAS guanyl releasing protein 2 (calcium and DAG-regulated)1
Fc gamma R-mediated phagocytosis(C = 91; O = 3; E = 0.17; R = 17.24; rawP = 0.0007; adjP = 0.0040)
Q8C3J5Dock2dedicator of cytokinesis 24
P97573Inpp5dinositol polyphosphate-5-phosphatase D2
P30009Marcksmyristoylated alanine rich protein kinase C substrate8
Insulin signaling pathway(C = 131; O = 3; E = 0.25; R = 11.97; rawP = 0.0021; adjP = 0.0089)
P22682CblCbl proto-oncogene, E3 ubiquitin protein ligase2
P62754Rps6ribosomal protein S61
P97573Inpp5dinositol polyphosphate-5-phosphatase D2
Focal adhesion(C = 186; O = 3; E = 0.36; R = 8.43; rawP = 0.0055; adjP = 0.0134)
Q62523Zyxzyxin1
Q10728Ppp1r12aprotein phosphatase 1, regulatory subunit 12A4
Q66H76Pxnpaxillin2
Phagosome(C = 185; O = 3; E = 0.35; R = 8.48; rawP = 0.0055; adjP = 0.0134)
O70257Stx7syntaxin 72
P35278Rab5cRAB5C, member RAS oncogene family1
Q91ZN1Coro1acoronin, actin binding protein 1A1
mTOR signaling pathway(C = 51; O = 2; E = 0.10; R = 20.50; rawP = 0.0043; adjP = 0.0134)
Q9WUT3Rps6ka2ribosomal protein S6 kinase polypeptide 21
P62754Rps6ribosomal protein S61
Endocytosis(C = 230; O = 3; E = 0.44; R = 6.82; rawP = 0.0099; adjP = 0.0153)
Q5FVC7Acap2ArfGAP with coiled coil, ankyrin repeat and PH domains 21
P22682CblCbl proto-oncogene, E3 ubiquitin protein ligase2
P35278Rab5cRAB5C, member RAS oncogene family1
B cell receptor signaling pathway(C = 75; O = 2; E = 0.14; R = 13.94; rawP = 0.0092; adjP = 0.0153)
P97573Inpp5dinositol polyphosphate-5-phosphatase D2
P81718Ptpn6protein tyrosine phosphatase, non-receptor type 63
Long-term potentiation(C = 69; O = 2; E = 0.13; R = 15.16; rawP = 0.0078; adjP = 0.0153)
Q10728Ppp1r12aprotein phosphatase 1, regulatory subunit 12A4
Q9WUT3Rps6ka2ribosomal protein S6 kinase polypeptide 21
C: the number of reference genes in the category, O: the number of genes in the gene set and also in the category, E: the expected number in the category, R: ratio of enrichment, rawP: p-value from hypergeometric test, adjP: p-value adjusted for multiple tests.

Share and Cite

MDPI and ACS Style

Tahir, M.; Arshid, S.; Fontes, B.; S. Castro, M.; Sidoli, S.; Schwämmle, V.; Luz, I.S.; Roepstorff, P.; Fontes, W. Phosphoproteomic Analysis of Rat Neutrophils Shows the Effect of Intestinal Ischemia/Reperfusion and Preconditioning on Kinases and Phosphatases. Int. J. Mol. Sci. 2020, 21, 5799. https://doi.org/10.3390/ijms21165799

AMA Style

Tahir M, Arshid S, Fontes B, S. Castro M, Sidoli S, Schwämmle V, Luz IS, Roepstorff P, Fontes W. Phosphoproteomic Analysis of Rat Neutrophils Shows the Effect of Intestinal Ischemia/Reperfusion and Preconditioning on Kinases and Phosphatases. International Journal of Molecular Sciences. 2020; 21(16):5799. https://doi.org/10.3390/ijms21165799

Chicago/Turabian Style

Tahir, Muhammad, Samina Arshid, Belchor Fontes, Mariana S. Castro, Simone Sidoli, Veit Schwämmle, Isabelle S. Luz, Peter Roepstorff, and Wagner Fontes. 2020. "Phosphoproteomic Analysis of Rat Neutrophils Shows the Effect of Intestinal Ischemia/Reperfusion and Preconditioning on Kinases and Phosphatases" International Journal of Molecular Sciences 21, no. 16: 5799. https://doi.org/10.3390/ijms21165799

APA Style

Tahir, M., Arshid, S., Fontes, B., S. Castro, M., Sidoli, S., Schwämmle, V., Luz, I. S., Roepstorff, P., & Fontes, W. (2020). Phosphoproteomic Analysis of Rat Neutrophils Shows the Effect of Intestinal Ischemia/Reperfusion and Preconditioning on Kinases and Phosphatases. International Journal of Molecular Sciences, 21(16), 5799. https://doi.org/10.3390/ijms21165799

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop