MAPRes: an efficient method to analyze protein sequence around post-translational modification sites

Functional switches are often regulated by dynamic protein modifications. Assessing protein functions, in vivo, and their functional switches remains still a great challenge in this age of development. An alternative methodology based on in silico procedures may facilitate assessing the multifunctionality of proteins and, in addition, allow predicting functions of those proteins that exhibit their functionality through transitory modifications. Extensive research is ongoing to predict the sequence of protein modification sites and analyze their dynamic nature. This study reports the analysis performed on phosphorylation, Phospho.ELM (version 3.0) and glycosylation, OGlycBase (version 6.0) data for mining association patterns utilizing a newly developed algorithm, MAPRes. This method, MAPRes (Mining Association Patterns among preferred amino acid residues in the vicinity of amino acids targeted for post-translational modifications), is based on mining association among significantly preferred amino acids of neighboring sequence environment and modification sites themselves. Association patterns arrived at by association pattern/rule mining were in significant conformity with the results of different approaches. However, attempts to analyze substrate sequence environment of phosphorylation sites catalyzed for Tyr kinases and the sequence data for O-GlcNAc modification were not successful, due to the limited data available. Using the MAPRes algorithm for developing an association among PTM site with its vicinal amino acids is a valid method with many potential uses: this is indeed the first method ever to apply the association pattern mining technique to protein post-translational modification data.     

Influence of the sequence environment and properties of neighboring amino acids on amino‐acetylation: Relevance for structure‐function analysis

Proteins function is regulated by co‐translational modifications and post‐translational modifications (PTMs) such as phosphorylation, glycosylation, and acetylation, which induce proteins to perform multiple tasks in a specified environment. Acetylation takes place post‐translationally on the ε‐amino group of Lys in histone proteins, allowing regulation of gene expression. Furthermore, amino group acetylation also occurs co‐translationally on Ser, Thr, Gly, Met, and Ala, possibly contributing to the stability of proteins. In this work, the influence of amino acids next to acetylated sites has been investigated by using MAPRes (Mining Association Patterns among preferred amino acid residues in the vicinity of amino acids targeted for PTMs). MAPRes was utilized to examine the sequence patterns vicinal to modified and non‐modified residues, taking into account their charge and polarity. The PTMs data were further sub‐divided according to their sub‐cellular location (nuclear, mitochondrial, and cytoplasmic), and their association patterns were mined. The association patterns mined by MAPRes for acetylated and non‐acetylated residues are consistent with the existing literature but also revealed novel patterns. These rules have been utilized to describe the acetylation and its effects on the protein structure‐function relationship. J. Cell. Biochem. 114: 874–887, 2013. © 2012 Wiley Periodicals, Inc.     

MAPRes: Mining association patterns among preferred amino acid residues in the vicinity of amino acids targeted for post-translational modifications

Post-translational modification (PTM) of a protein is an important event in regulating cellular functions. An algorithm, MAPRes, has been developed for mining associations among PTM sites and the preferred amino acids in their vicinity. The algorithm has been implemented to O-glycosylation and O-phosphorylation data (phosphorylated/glycosylated Ser/Thr/Tyr). The association patterns mined by MAPRes demonstrate significant correlations and the results are in conformity with the existing methods. These association rules/patterns will be helpful in predicting the sequences/motifs involved for specific PTMs in proteins.     

Pim kinase substrate identification and specificity

The Pim family of Ser/Thr kinases has been implicated in the process of lymphomagenesis and cell survival. Known substrates of Pim kinases are few and poorly characterized. In this study we set out to identify novel Pim-2 substrates using the Kinase Substrate Tracking and Elucidation (KESTREL) approach. Two potential substrates, eukaryotic initiation factor 4B (eIF4B) and apoptosis inhibitor 5 (API-5), were identified from rat thymus extracts. Sequence comparison of the Pim-2 kinase phosphorylation sites of eIF4B and mouse BAD, the only other known Pim-2 substrate, revealed conserved amino acids preceding the phosphorylated serine residue. Stepwise replacement of the conserved residues produced a consensus sequence for Pim kinase recognition: RXRHXS. Pim-1 and Pim-2 catalyzed the phosphorylation of this recognition sequence 20-fold more efficiently than the original (K/R-K/R-R-K/R-L-S/T-a; a = small chain amino acid) Pim-1 phosphorylation site. The identification of the novel Pim kinase consensus sequence provides a more sensitive and versatile peptide based assay for screening modulators of Pim kinase activity.     

Identification of major ERK-related phosphorylation sites in Gab1

Gab1 (Grb2-associated binder1) belongs to a family of multifunctional docking proteins that play a central role in the integration of receptor tyrosine kinase (RTK) signaling, i.e., mediating cellular growth response, transformation, and apoptosis. In addition to RTK-specific tyrosine phosphorylation, these docking proteins also can be phosphorylated on serine/threonine residues affecting signal transduction. Since serine and threonine phosphorylation are capable of modulating the initial signal one major task to elucidate signal transduction via Gab1 is to determine the exact localization of distinct phosphorylation sites. To address this question in this report we examined extracellular signal-regulated kinases 1/2 (ERK) specific serine/threonine phosphorylation of the entire Gab1 engaged in insulin signaling in more detail in vitro. To elucidate the ERK1/2-specific phosphorylation pattern of Gab1, we used phosphopeptide mapping by two-dimensional HPLC analysis. Subsequently, phosphorylated serine/threonine residues were identified by sequencing the separated phosphopeptides using matrix assisted laser desorption ionization mass spectrometry (MALDI-MS) and Edman degradation. Our results demonstrate that ERK1/2 phosphorylate Gab1 at six serine/threonine residues (T312, S381, S454, T476, S581, S597) in consensus motifs for MAP kinase phosphorylation. Serine residues S454, S581, S597, and threonine residue T476 represent nearly 80% of overall incorporated phosphate. These sites are located adjacent to src homology region-2 (SH2) binding motifs (YVPM-motif: Y447, Y472, Y619) specific for the phosphatidylinositol 3kinase (PI3K). The biological role of identified phosphorylation sites was proven by PI3K and Akt activity in intact cells. These data demonstrate that ERK1/2 modulate insulin action via Gab1 by targeting serine and threonine residues beside YXXM motifs. Accordingly, insulin signaling is blocked at the level of PI3K.     

Evolved to be active: sulfate ions define substrate recognition sites of CK2alpha and emphasise its exceptional role within the CMGC family of eukaryotic protein kinases

CK2alpha is the catalytic subunit of protein kinase CK2 and a member of the CMGC family of eukaryotic protein kinases like the cyclin-dependent kinases, the MAP kinases and glycogen-synthase kinase 3. We present here a 1.6 A resolution crystal structure of a fully active C-terminal deletion mutant of human CK2alpha liganded by two sulfate ions, and we compare this structure systematically with representative structures of related CMGC kinases. The two sulfate anions occupy binding pockets at the activation segment and provide the structural basis of the acidic consensus sequence S/T-D/E-X-D/E that governs substrate recognition by CK2. The anion binding sites are conserved among those CMGC kinases. In most cases they are neutralized by phosphorylation of a neighbouring threonine or tyrosine side-chain, which triggers conformational changes for regulatory purposes. CK2alpha, however, lacks both phosphorylation sites at the activation segment and structural plasticity. Here the anion binding sites are functionally changed from regulation to substrate recognition. These findings underline the exceptional role of CK2alpha as a constitutively active enzyme within a family of strictly controlled protein kinases.     

Specificity determinants of substrate recognition by the protein kinase DYRK1A

DYRK1A is a dual-specificity protein kinase that is thought to be involved in brain development. We identified a single phosphorylated amino acid residue in the DYRK substrate histone H3 (threonine 45) by mass spectrometry, phosphoamino acid analysis, and protein sequencing. Exchange of threonine 45 for alanine abolished phosphorylation of histone H3 by DYRK1A and by the related kinases DYRK1B, DYRK2, and DYRK3 but not by CLK3. In order to define the consensus sequence for the substrate specificity of DYRK1A, a library of 300 peptides was designed in variation of the H3 phosphorylation site. Evaluation of the phosphate incorporation into these peptides identified DYRK1A as a proline-directed kinase with a phosphorylation consensus sequence (RPX(S/T)P) similar to that of ERK2 (PX(S/T)P). A peptide designed after the optimal substrate sequence (DYRKtide) was efficiently phosphorylated by DYRK1A (K(m) = 35 microM) but not by ERK2. Both ERK2 and DYRK1A phosphorylated myelin basic protein, whereas only ERK2, but not DYRK1A, phosphorylated the mitogen-activated protein kinase substrate ELK-1. This marked difference in substrate specificity between DYRK1A and ERK2 can be explained by the requirement for an arginine at the P -3 site of DYRK substrates and its presumed interaction with aspartate 247 conserved in all DYRKs.     

Positive and negative regulation of Raf kinase activity and function by phosphorylation.

Activating and inhibitory phosphorylation mechanisms play an essential role in regulating Raf kinase activity. Here we demonstrate that phosphorylation of C-Raf in the kinase activation loop (residues T491 and S494) is necessary, but not sufficient, for activation. C-Raf has additional activating phosphorylation sites at S338 and Y341. Mutating all four of these residues to acidic residues, S338D/Y341D/T491E/S494D (DDED), in C-Raf results in constitutive activity. However, acidic residue substitutions at the corresponding activation loop sites in B-Raf are sufficient to confer constitutive activity. B-Raf and C-Raf also utilize similar inhibitory phosphorylation mechanisms to regulate kinase activity. B-Raf has multiple inhibitory phosphorylation sites necessary for full kinase inhibition where C-Raf requires only one. We examined the functional significance of these inhibitory and activating phosphorylations in Caenorhabditis elegans lin-45 Raf. Eliminating the inhibitory phosphorylation or mimicking activating phosphorylation sites is sufficient to confer constitutive activity upon lin-45 Raf and induce multi-vulva phenotypes in C.elegans. Our results demonstrate that different members of the Raf family kinases have both common and distinct phosphorylation mechanisms to regulate kinase activity and biological function.     

Bioinformatic prediction and confirmation of beta-adducin as a novel substrate of glycogen synthase kinase 3

It is important to identify the true substrates of protein kinases because this illuminates the primary function of any kinase. Here, we used bioinformatics and biochemical validation to identify novel brain substrates of the Ser/Thr kinase glycogen synthase kinase 3 (GSK3). Briefly, sequence databases were searched for proteins containing a conserved GSK3 phosphorylation consensus sequence ((S/T)PXX(S/T)P or (S/T)PXXX(S/T)P), as well as other criteria of interest (e.g. brain proteins). Importantly, candidates were highlighted if they had previously been reported to be phosphorylated at these sites by large-scale phosphoproteomic studies. These criteria identified the brain-enriched cytoskeleton-associated protein β-adducin as a likely substrate of GSK3. To confirm this experimentally, it was cloned and subjected to a combination of cell culture and in vitro kinase assays that demonstrated direct phosphorylation by GSK3 in vitro and in cells. Phosphosites were mapped to three separate regions near the C terminus and confirmed using phosphospecific antibodies. Prior priming phosphorylation by Cdk5 enhanced phosphorylation by GSK3. Expression of wild type, but not non-phosphorylatable (GSK3 insensitive), β-adducin increased axon and dendrite elongation in primary cortical neurons. Therefore, phosphorylation of β-adducin by GSK3 promotes efficient neurite outgrowth in neurons.     

Multisite phosphorylation of a synthetic peptide derived from the carboxyl terminus of the ribosomal protein S6

The synthetic peptide AKRRRLSSLRASTSKSESSQK (S6-21) which corresponds to the carboxyl-terminal 21 amino acids of human ribosomal protein S6 was synthesized and tested as a substrate for S6/H4 kinase purified from human placenta. The specific activity of the enzyme with the synthetic peptide and 40 S ribosomes was 45 and 23 nmol/min/mg, respectively. The S6/H4 kinase activity with S6-21 was greater than the enzyme activity with any other substrate tested, including histones, protamine, and casein and several other synthetic peptides. The phosphorylation of the peptide was not inhibited by inhibitors of several other proteins kinases. S6/H4 kinase catalyzed the phosphorylation of three major sites in the synthetic peptide and the 40 S ribosomes. A fourth site in S6-21 was phosphorylated more slowly. The principal phosphorylation sites were serines in the acidic carboxyl-terminal domain of the peptide. A serine (Ser-7 or -8) in the amino-terminal domain was phosphorylated at approximately 25% the rate of the carboxyl-terminal domain serines. The data suggest that multiple S6 kinases may be required to phosphorylate S6 at all five sites which are modified in vivo.     

Docking sites on substrate proteins direct extracellular signal-regulated kinase to phosphorylate specific residues

Mitogen-activated protein (MAP) kinases such as extracellular signal-regulated kinase (ERK) are important signaling proteins that phosphorylate (S/T)P sites in many different protein substrates. ERK binding to substrate proteins is mediated by docking sites including the FXFP motif and the D-domain. We characterized the sequence of amino acids that can constitute the FXFP motif using peptide and protein substrates. Substitutions of the phenylalanines at positions 1 and 3 had significant effects, indicating that these phenylalanines provide substantial binding affinity, whereas substitutions of the residues at positions 2 and 4 had less effect. The FXFP and D-domain docking sites were analyzed in a variety of positions and arrangements in the proteins ELK-1 and KSR-1. Our results indicate that the FXFP and D-domain docking sites form a flexible, modular system that has two functions. First, the affinity of a substrate for ERK can be regulated by the number, type, position, and arrangement of docking sites. Second, in substrates with multiple potential phosphorylation sites, docking sites can direct phosphorylation of specific (S/T)P residues. In particular, the FQFP motif of ELK-1 is necessary and sufficient to direct phosphorylation of serine 383, whereas the D-domain directs phosphorylation of other (S/T)P sites in ELK-1.     

Detection of a variety of Ser/Thr protein kinases using a synthetic peptide with multiple phosphorylation sites

A novel peptide with multiple phosphorylation sites, which we designated as multide, was developed to detect a wide variety of protein kinases in crude cell extracts. Multide, KKRKSSLRRWSPLTPRQMSFDC, has been designed to contain consensus sequences for various Ser/Thr protein kinases including cAMP-dependent protein kinase, protein kinase C, MAP kinases, and Ca(2+)/calmodulin-dependent protein kinases in a single peptide. In-gel protein kinase assay using multide was found to be very useful for analyzing the activities of protein kinases that are altered in response to various extracellular stimuli. The substrate specificities of the protein kinases thus detected were further determined by using five multide analogs with different phosphorylation sites.     

Regulation of VASP by phosphorylation

Phosphorylations control all aspects of vasodilator-stimulated phospho-protein (VASP) function. Mapped phosphorylation sites include Y39, S157, S239, T278, and S322, and multiple kinases have been shown to mediate their phosphorylation. Recently, Protein Kinase D1 (PKD1) as a direct kinase for S157 and S322 joined this group. While S157 phosphorylation generally seems to serve as a signal for membrane localization, phosphorylations at S322 or at S239 and T278 have opposite effects on F-actin accumulation. In migrating cells, S322 phosphorylation increases filopodia numbers and length, while S239/T278 phosphorylations decrease these and also disrupt formation of focal adhesions. Therefore, the kinases mediating these phosphorylations can be seen as switches needed to facilitate cell motility.     

Construction of human activity‐based phosphorylation networks

The landscape of human phosphorylation networks has not been systematically explored, representing vast, unchartered territories within cellular signaling networks. Although a large number of in vivo phosphorylated residues have been identified by mass spectrometry (MS)‐based approaches, assigning the upstream kinases to these residues requires biochemical analysis of kinase‐substrate relationships (KSRs). Here, we developed a new strategy, called CEASAR, based on functional protein microarrays and bioinformatics to experimentally identify substrates for 289 unique kinases, resulting in 3656 high‐quality KSRs. We then generated consensus phosphorylation motifs for each of the kinases and integrated this information, along with information about in vivo phosphorylation sites determined by MS, to construct a high‐resolution map of phosphorylation networks that connects 230 kinases to 2591 in vivo phosphorylation sites in 652 substrates. The value of this data set is demonstrated through the discovery of a new role for PKA downstream of Btk (Bruton's tyrosine kinase) during B‐cell receptor signaling. Overall, these studies provide global insights into kinase‐mediated signaling pathways and promise to advance our understanding of cellular signaling processes in humans.     

Substrate specificities of tyrosine-specific protein kinases toward cytoskeletal proteins in vitro

We have previously reported that fodrin (beta subunit), tubulin (alpha subunit) and microtubule-associated proteins (MAPs; MAP2 and tau) are good substrates for the purified insulin receptor kinase (Kadowaki, T., Nishida, E., Kasuga, M., Akiyama, T., Takaku, F., Ishikawa, M., Sakai, H., Kathuria, S., and Fujita-Yamaguchi, Y. (1985) Biochem. Biophys. Res. Commun. 127, 493-500 and Kadowaki, T., Fujita-Yamaguchi, Y., Nishida, E., Takaku, F., Akiyama, T., Kathuria, S., Akanuma, Y., and Kasuga, M. (1985) J. Biol. Chem. 260, 4016-4020). In this study, to investigate the substrate specificities of tyrosine kinases, we have examined the actions of the purified epidermal growth factor (EGF) receptor kinase and Rous sarcoma virus src kinase on purified microfilament- and microtubule-related proteins. Among microfilament-related proteins examined, the purified EGF receptor kinase phosphorylated the beta subunit, but not the alpha subunit, of fodrin on tyrosine residues with a Km below the micromolar range. The fodrin phosphorylation by the EGF receptor kinase was markedly inhibited by F-actin. In contrast, the purified src kinase preferentially phosphorylated the alpha subunit of fodrin on tyrosine residues. Fodrin phosphorylation by the src kinase was not inhibited by F-actin. Among microtubule proteins examined, MAP2 was the best substrate for the EGF receptor kinase. By contrast, src kinase favored phosphorylation of tubulin as compared to MAP2. The peptide mapping of MAP2 phosphorylated by the EGF receptor kinase and by the insulin receptor kinase produced very similar patterns of phosphopeptides, while that of MAP2 phosphorylated by the src kinase gave a distinctly different pattern. When the phosphorylation of the tubulin subunits was examined, the EGF receptor kinase preferred beta subunit to alpha subunit, but the src kinase phosphorylated both alpha and beta subunits to a similar extent. These results, together with our previous results, indicate that the substrate specificities of the EGF receptor kinase and the insulin receptor kinase are very similar, but not identical, while that of the src kinase is distinctly different from that of these growth factor receptor kinases.     

Comparison of Peptide Array Substrate Phosphorylation of c-Raf and Mitogen Activated Protein Kinase Kinase Kinase 8

Kinases are pivotal regulators of cellular physiology. The human genome contains more than 500 putative kinases, which exert their action via the phosphorylation of specific substrates. The determinants of this specificity are still only partly understood and as a consequence it is difficult to predict kinase substrate preferences from the primary structure, hampering the understanding of kinase function in physiology and prompting the development of technologies that allow easy assessment of kinase substrate consensus sequences. Hence, we decided to explore the usefulness of phosphorylation of peptide arrays comprising of 1176 different peptide substrates with recombinant kinases for determining kinase substrate preferences, based on the contribution of individual amino acids to total array phosphorylation. Employing this technology, we were able to determine the consensus peptide sequences for substrates of both c-Raf and Mitogen Activated Protein Kinase Kinase Kinase 8, two highly homologous kinases with distinct signalling roles in cellular physiology. The results show that although consensus sequences for these two kinases identified through our analysis share important chemical similarities, there is still some sequence specificity that could explain the different biological action of the two enzymes. Thus peptide arrays are a useful instrument for deducing substrate consensus sequences and highly homologous kinases can differ in their requirement for phosphorylation events.     

Phosphorylation of human microtubule-associated protein tau by protein kinases of the AGC subfamily

Intraneuronal inclusions made of hyperphosphorylated microtubule-associated protein tau are a defining neuropathological characteristic of Alzheimer's disease, and of several other neurodegenerative disorders. Many phosphorylation sites in tau are S/TP sites that flank the microtubule-binding repeats. Others are KXGS motifs in the repeats. One site upstream of the repeats lies in a consensus sequence for AGC kinases. This site (S214) is believed to play an important role in the events leading from normal, soluble to filamentous, insoluble tau. Here, we show that all AGC kinases tested phosphorylated S214. RSK1 and p70 S6 kinase also phosphorylated the neighbouring T212, a TP site that conforms weakly to the AGC kinase consensus sequence. MSK1 phosphorylated S214, as well as S262, a KXGS site in the first repeat, and S305 in the second repeat.     

Regulatory roles of conserved phosphorylation sites in the activation T-loop of the MAP kinase ERK1

The catalytic domains of most eukaryotic protein kinases are highly conserved in their primary structures. Their phosphorylation within the well-known activation T-loop, a variable region between protein kinase catalytic subdomains VII and VIII, is a common mechanism for stimulation of their phosphotransferase activities. Extracellular signal–regulated kinase 1 (ERK1), a member of the extensively studied mitogen-activated protein kinase (MAPK) family, serves as a paradigm for regulation of protein kinases in signaling modules. In addition to the well-documented T202 and Y204 stimulatory phosphorylation sites in the activation T-loop of ERK1 and its closest relative, ERK2, three additional flanking phosphosites have been confirmed (T198, T207, and Y210 from ERK1) by high-throughput mass spectrometry. In vitro kinase assays revealed the functional importance of T207 and Y210, but not T198, in negatively regulating ERK1 catalytic activity. The Y210 site could be important for proper conformational arrangement of the active site, and a Y210F mutant could not be recognized by MEK1 for phosphorylation of T202 and Y204 in vitro. Autophosphorylation of T207 reduces the catalytic activity and stability of activated ERK1. We propose that after the activation of ERK1 by MEK1, subsequent slower phosphorylation of the flanking sites results in inhibition of the kinase. Because the T207 and Y210 phosphosites of ERK1 are highly conserved within the eukaryotic protein kinase family, hyperphosphorylation within the kinase activation T-loop may serve as a general mechanism for protein kinase down-regulation after initial activation by their upstream kinases.     

Identification of a CK2 phosphorylation site in mdm2.

Mdm2 is a cellular oncoprotein the most obvious function of which is the down-regulation of the growth suppressor protein p53. It represents a highly phosphorylated protein but only little is yet known about the sites phosphorylated in vivo, the kinases that are responsible for the phosphorylation or the functional relevance of the phosphorylation status. Recently, we have shown that mdm2 is a good substrate for protein kinase CK2 at least in vitro. Computer analysis of the primary amino acid sequence of mdm2 revealed 19 putative CK2 phosphorylation sites. By using deletion mutants of mdm2 and a peptide library we identified the serine residue at position 269 which lies within a canonical CK2 consensus sequence (EGQELSDEDDE) as the most important CK2 phosphorylation site. Moreover, by using the mdm2 S269A mutant for in vitro phosphorylation assays this site was shown to be phosphorylated by CK2. Binding studies revealed that phosphorylation of mdm2 at S269 does not have any influence on the binding of p53 to mdm2.