Quantitative Site-specific Phosphorylation Dynamics of Human Protein Kinases during Mitotic Progression

Reversible protein phosphorylation is a key regulatory mechanism of mitotic progression. Importantly, protein kinases themselves are also regulated by phosphorylation-dephosphorylation processes; hence, phosphorylation dynamics of kinases hold a wealth of information about phosphorylation networks. Here, we investigated the site-specific phosphorylation dynamics of human kinases during mitosis using synchronization of HeLa suspension cells, kinase enrichment, and high resolution mass spectrometry. In biological triplicate analyses, we identified 206 protein kinases and more than 900 protein kinase phosphorylation sites, including 61 phosphorylation sites on activation segments, and quantified their relative abundances across three specific mitotic stages. Around 25% of the kinase phosphorylation site ratios were found to be changed by at least 50% during mitotic progression. Further network analysis of jointly regulated kinase groups suggested that Cyclin-dependent kinase- and mitogen-activated kinase-centered interaction networks are coordinately down- and up-regulated in late mitosis, respectively. Importantly, our data cover most of the already known mitotic kinases and, moreover, identify attractive candidates for future studies of phosphorylation-based mitotic signaling. Thus, the results of this study provide a valuable resource for cell biologists and provide insight into the system properties of the mitotic phosphokinome.Reversible phosphorylation is a ubiquitous posttranslational protein modification that is involved in the regulation of almost all biological processes (1–3). In human, 518 protein kinases have been identified in the genome that phosphorylate the majority of cellular proteins and increase the diversity of the proteome by severalfold (4). Addition of a phosphate group to a protein can alter its structural, catalytic, and functional properties; hence, kinases require tight regulation to avoid unspecific phosphorylation, which can be deleterious to cells (5–7). As a result, cells use a variety of mechanisms to ensure proper regulation of kinase activities (8). Importantly, most kinases are also in turn regulated through autophosphorylation and phosphorylation by other kinases, thus generating complex phosphorylation networks. In particular, phosphorylation on activation segments is a common mechanism to modulate kinase activities (9–11), but additional phosphorylation sites are also frequently required for fine tuning of kinase localizations and functions (12). Some kinases contain phosphopeptide binding domains that recognize prephosphorylated sites on other kinases, resulting in processive phosphorylation and/or targeting of kinases to distinct cellular locations (13–16). Because such priming phosphorylation events depend on the activities of the priming kinases, these motifs act as conditional docking sites and restrict the interaction with docking kinases to a particular point in time and physiological state. In addition, phosphorylation sites may act through combinatorial mechanisms or through cross-talk with other posttranslational modifications (PTMs)¹ (17, 18), thus further increasing the complexity of kinase regulatory networks.Regulation of kinases is of particular interest in mitosis as most of the mitotic events are regulated by reversible protein phosphorylation (19). During mitosis, error-free segregation of sister chromatids into the two daughter cells is essential to ensure genomic stability. Physically, this process is carried out by the mitotic spindle, a highly dynamic microtubule-based structure. After entry into mitosis, the major microtubule-organizing centers in animal cells, the centrosomes, start to increase microtubule nucleation and move to opposite poles of the cell. Throughout prometaphase, microtubules emanating from centrosomes are captured by kinetochores, protein complexes assembled on centromeric chromosomal DNA. This eventually leads to the alignment of all chromosomes in a metaphase plate. Because proper bipolar attachment of chromosomes to spindle microtubules is essential for the correct segregation of chromosomes, this critical step is monitored by a signaling pathway known as the spindle assembly checkpoint (SAC) (20). This checkpoint is silenced only after all chromosomes have attached to the spindle in a bioriented fashion, resulting in the synchronous segregation of sister chromatids during anaphase. Simultaneously, a so-called central spindle is formed between the separating chromatids, and the formation of a contractile ring initiates cytokinesis. Finally, in telophase, the chromosomes decondense and reassemble into nuclei, whereas remnants of the central spindle form the midbody, marking the site of abscission. Cyclin-dependent kinase 1 (Cdk1), an evolutionarily conserved master mitotic kinase, is activated prior to mitosis and initiates most of the mitotic events. Cdk1 works in close association with other essential mitotic kinases such as Plk1, Aurora A, and Aurora B for the regulation of mitotic progression (19, 21–24). Plk1 and Aurora kinases dynamically localize to different subcellular locations to perform multiple functions during mitosis and are phosphorylated at several conserved sites. Although little is known about the precise roles of these phosphorylation sites, emerging data indicate that they are involved in regulating localization-specific functions (25, 26). Furthermore, the kinases Bub1, BubR1, and TTK (Mps1) and kinases of the Nek family play important roles in maintaining the fidelity and robustness of mitosis (19). Recently, a genome-wide RNA-mediated interference screen identified M phase phenotypes for many kinases that have not previously been implicated in cell cycle functions, indicating that additional kinases have important mitotic functions (27).Although protein phosphorylation plays a pivotal role in the regulation of cellular networks, many phosphorylation events remain undiscovered mainly because of technical limitations (28). The advent of mass spectrometry-based proteomics along with developments in phosphopeptide enrichment methods has enabled large scale global phosphoproteomics studies (29, 30). However, the number of phosphorylation sites identified on kinases is limited compared with other proteins because of their frequently low expression levels. To overcome this problem, small inhibitor-based kinase enrichment strategies were developed, resulting in the identification of more than 200 kinases from HeLa cell lysates (31, 32). This method was also used recently to compare the phosphokinomes during S phase and M phase of the cell cycle, resulting in the identification of several hundreds of M phase-specific kinase phosphorylation sites (31). In the present study, we address the dynamics of the phosphokinome during mitotic progression using large scale cell synchronization at three distinct mitotic stages, small inhibitor-based kinase enrichment, and stable isotope labeling by amino acids in cell culture (SILAC)-based quantitative mass spectrometry. Thus, we determined the mitotic phosphorylation dynamics of more than 900 kinase phosphorylation sites and identified distinctly regulated kinase interaction networks. Our results provide a valuable resource for the dynamics of the kinome during mitotic progression and give insight into the system properties of kinase interaction networks.     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Large-scale Proteomics Analysis of the Human Kinome

Members of the human protein kinase superfamily are the major regulatory enzymes involved in the activity control of eukaryotic signal transduction pathways. As protein kinases reside at the nodes of phosphorylation-based signal transmission, comprehensive analysis of their cellular expression and site-specific phosphorylation can provide important insights into the architecture and functionality of signaling networks. However, in global proteome studies, low cellular abundance of protein kinases often results in rather minor peptide species that are occluded by a vast excess of peptides from other cellular proteins. These analytical limitations create a rationale for kinome-wide enrichment of protein kinases prior to mass spectrometry analysis. Here, we employed stable isotope labeling by amino acids in cell culture (SILAC) to compare the binding characteristics of three kinase-selective affinity resins by quantitative mass spectrometry. The evaluated pre-fractionation tools possessed pyrido[2,3-d]pyrimidine-based kinase inhibitors as immobilized capture ligands and retained considerable subsets of the human kinome. Based on these results, an affinity resin displaying the broadly selective kinase ligand VI16832 was employed to quantify the relative expression of more than 170 protein kinases across three different, SILAC-encoded cancer cell lines. These experiments demonstrated the feasibility of comparative kinome profiling in a compact experimental format. Interestingly, we found high levels of cytoplasmic and low levels of receptor tyrosine kinases in MV4–11 leukemia cells compared with the adherent cancer lines HCT116 and MDA-MB-435S. The VI16832 resin was further exploited to pre-fractionate kinases for targeted phosphoproteomics analysis, which revealed about 1200 distinct phosphorylation sites on more than 200 protein kinases. This hitherto largest survey of site-specific phosphorylation across the kinome significantly expands the basis for functional follow-up studies on protein kinase regulation. In conclusion, the straightforward experimental procedures described here enable different implementations of kinase-selective proteomics with considerable potential for future signal transduction and kinase drug target analysis.Reversible protein phosphorylation represents the most common type of post-translational modification (PTM)¹ in eukaryotic organisms. A plethora of studies on a large variety of proteins have established that site-specific phosphorylation events fulfill key functions in the activity control of signaling cascades and networks (1). Cellular protein phosphorylation is controlled by more than 500 members of the protein kinase superfamily, which comprises one of the largest enzyme families encoded by the human genome (2). Protein kinases represent the key elements in phosphorylation-based signal transmission. Aberrant protein kinase expression and/or activity, often because of gene amplification or mutational changes, is involved in pathological processes leading to malignant transformation and tumor development (3). Therefore, protein kinases have emerged as a major class of drug targets for therapeutic intervention (4–6). Given the diversity of molecular mechanisms related to de-regulated kinase function in human cancers, proteomic approaches could significantly enhance our understanding of disease-relevant kinase function and also help to optimize and adjust therapeutic strategies. In addition to assessing protein expression, the analysis of site-specific phosphorylations on protein kinases is of particular relevance, as these PTMs can be indicative of their cellular catalytic activities (7, 8). Protein kinases can not only modulate each other''s functions and activities through site-specific phosphorylation events, but often also undergo site-specific autophosphorylation once they get activated (9). Thus, the comprehensive assessment of kinase-derived phosphopeptides can provide important insights into the regulation of these key players in phosphorylation-controlled signaling.Regulatory enzymes such as protein kinases are often expressed at low cellular levels. This can impede their detection by LC-MS in highly complex peptide mixtures derived from total cell or tissue extracts. These analytical challenges are further aggravated in phosphoproteomic experiments due to the fact that many phosphopeptide species result from sub-stoichiometric phosphorylation events (10). Consequently, phosphopeptide isolation methods have proven to be essential. Among others, techniques such as immobilized metal affinity chromatography or enrichment by means of titanium dioxide (TiO₂)-coated beads have found widespread use in MS-based phosphoproteomics (11–13). In addition, to reduce initial sample complexity, either protein fractionation by gel electrophoresis or peptide separation by strong cation exchange chromatography is typically included in contemporary phosphoproteomics workflows (14–16). These separation techniques in combination with LC-MS on state-of-the-art mass spectrometers enabled the identification of thousands of phosphorylation sites from total cellular extracts (15, 17, 18). Despite these impressive advances, such large-scale efforts require considerable instrument time, and the current methodology is still not comprehensive across the full dynamic range of the entire phosphoproteome. This creates the rationale for sub-proteome analyses to achieve high coverage and analytical sensitivity, which is particularly relevant for members of the protein kinase enzyme family.To date, the only pre-fractionation techniques permitting the enrichment of more than a few protein kinases are affinity capture methods relying on immobilized and kinase-selective small molecule inhibitors (19–21). We and others have demonstrated that combinations of such kinase inhibitor resins efficiently pre-fractionate kinases for subsequent phosphorylation analysis (7, 22, 23). Ideally, capture molecules for kinase proteomics have two properties. First, they should exhibit high non-selectivity within the kinase superfamily. Second, they should efficiently discriminate between protein kinases and other classes of cellular proteins under the biochemical conditions of the pre-fractionation procedure.In our efforts to characterize affinity reagents fulfilling these criteria, we quantitatively compared a selection of immobilized pyrido[2,3-d]pyrimidine-based inhibitors with respect to their proteome-wide kinase binding properties. Based on this assessment, an affinity matrix displaying the small molecule VI16832 was used as an enrichment tool for the comparative expression analysis of protein kinases in different cancer cell lines. The highly efficient VI16832 affinity resin further enabled a large-scale phosphoproteomics survey resulting in the identification and confident assignment of about 1200 phosphorylation sites on more than 200 distinct protein kinases.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]     

Role of Partitioning-defective 1/Microtubule Affinity-regulating Kinases in the Morphogenetic Activity of Helicobacter pylori CagA

Helicobacter pylori CagA plays a key role in gastric carcinogenesis. Upon delivery into gastric epithelial cells, CagA binds and deregulates SHP-2 phosphatase, a bona fide oncoprotein, thereby causing sustained ERK activation and impaired focal adhesions. CagA also binds and inhibits PAR1b/MARK2, one of the four members of the PAR1 family of kinases, to elicit epithelial polarity defect. In nonpolarized gastric epithelial cells, CagA induces the hummingbird phenotype, an extremely elongated cell shape characterized by a rear retraction defect. This morphological change is dependent on CagA-deregulated SHP-2 and is thus thought to reflect the oncogenic potential of CagA. In this study, we investigated the role of the PAR1 family of kinases in the hummingbird phenotype. We found that CagA binds not only PAR1b but also other PAR1 isoforms, with order of strength as follows: PAR1b > PAR1d ≥ PAR1a > PAR1c. Binding of CagA with PAR1 isoforms inhibits the kinase activity. This abolishes the ability of PAR1 to destabilize microtubules and thereby promotes disassembly of focal adhesions, which contributes to the hummingbird phenotype. Consistently, PAR1 knockdown potentiates induction of the hummingbird phenotype by CagA. The morphogenetic activity of CagA was also found to be augmented through inhibition of non-muscle myosin II. Because myosin II is functionally associated with PAR1, perturbation of PAR1-regulated myosin II by CagA may underlie the defect of rear retraction in the hummingbird phenotype. Our findings reveal that CagA systemically inhibits PAR1 family kinases and indicate that malfunctioning of microtubules and myosin II by CagA-mediated PAR1 inhibition cooperates with deregulated SHP-2 in the morphogenetic activity of CagA.Infection with Helicobacter pylori strains bearing cagA (cytotoxin-associated gene A)-positive strains is the strongest risk factor for the development of gastric carcinoma, the second leading cause of cancer-related death worldwide (1–3). The cagA gene is located within a 40-kb DNA fragment, termed the cag pathogenicity island, which is specifically present in the genome of cagA-positive H. pylori strains (4–6). In addition to cagA, there are ∼30 genes in the cag pathogenicity island, many of which encode a bacterial type IV secretion system that delivers the cagA-encoded CagA protein into gastric epithelial cells (7–10). Upon delivery into gastric epithelial cells, CagA is localized to the plasma membrane, where it undergoes tyrosine phosphorylation at the C-terminal Glu-Pro-Ile-Tyr-Ala motifs by Src family kinases or c-Abl kinase (11–14). The C-terminal Glu-Pro-Ile-Tyr-Ala-containing region of CagA is noted for the structural diversity among distinct H. pylori isolates. Oncogenic potential of CagA has recently been confirmed by a study showing that systemic expression of CagA in mice induces gastrointestinal and hematological malignancies (15).When expressed in gastric epithelial cells, CagA induces morphological transformation termed the hummingbird phenotype, which is characterized by the development of one or two long and thin protrusions resembling the beak of the hummingbird. It has been thought that the hummingbird phenotype is related to the oncogenic action of CagA (7, 16–19). Pathophysiological relevance for the hummingbird phenotype in gastric carcinogenesis has recently been provided by the observation that infection with H. pylori carrying CagA with greater ability to induce the hummingbird phenotype is more closely associated with gastric carcinoma (20–23). Elevated motility of hummingbird cells (cells showing the hummingbird phenotype) may also contribute to invasion and metastasis of gastric carcinoma.In host cells, CagA interacts with the SHP-2 phosphatase, C-terminal Src kinase, and Crk adaptor in a tyrosine phosphorylation-dependent manner (16, 24, 25) and also associates with Grb2 adaptor and c-Met in a phosphorylation-independent manner (26, 27). Among these CagA targets, much attention has been focused on SHP-2 because the phosphatase has been recognized as a bona fide oncoprotein, gain-of-function mutations of which are found in various human malignancies (17, 18, 28). Stable interaction of CagA with SHP-2 requires CagA dimerization, which is mediated by a 16-amino acid CagA-multimerization (CM)² sequence present in the C-terminal region of CagA (29). Upon complex formation, CagA aberrantly activates SHP-2 and thereby elicits sustained ERK MAP kinase activation that promotes mitogenesis (30). Also, CagA-activated SHP-2 dephosphorylates and inhibits focal adhesion kinase (FAK), causing impaired focal adhesions. It has been shown previously that both aberrant ERK activation and FAK inhibition by CagA-deregulated SHP-2 are involved in induction of the hummingbird phenotype (31).Partitioning-defective 1 (PAR1)/microtubule affinity-regulating kinase (MARK) is an evolutionally conserved serine/threonine kinase originally isolated in C. elegans (32–34). Mammalian cells possess four structurally related PAR1 isoforms, PAR1a/MARK3, PAR1b/MARK2, PAR1c/MARK1, and PAR1d/MARK4 (35–37). Among these, PAR1a, PAR1b, and PAR1c are expressed in a variety of cells, whereas PAR1d is predominantly expressed in neural cells (35, 37). These PAR1 isoforms phosphorylate microtubule-associated proteins (MAPs) and thereby destabilize microtubules (35, 38), allowing asymmetric distribution of molecules that are involved in the establishment and maintenance of cell polarity.In polarized epithelial cells, CagA disrupts the tight junctions and causes loss of apical-basolateral polarity (39, 40). This CagA activity involves the interaction of CagA with PAR1b/MARK2 (19, 41). CagA directly binds to the kinase domain of PAR1b in a tyrosine phosphorylation-independent manner and inhibits the kinase activity. Notably, CagA binds to PAR1b via the CM sequence (19). Because PAR1b is present as a dimer in cells (42), CagA may passively homodimerize upon complex formation with the PAR1 dimer via the CM sequence, and this PAR1-directed CagA dimer would form a stable complex with SHP-2 through its two SH2 domains.Because of the critical role of CagA in gastric carcinogenesis (7, 16–19), it is important to elucidate the molecular basis underlying the morphogenetic activity of CagA. In this study, we investigated the role of PAR1 isoforms in induction of the hummingbird phenotype by CagA, and we obtained evidence that CagA-mediated inhibition of PAR1 kinases contributes to the development of the morphological change by perturbing microtubules and non-muscle myosin II.