KRSA: An R package and R Shiny web application for an end-to-end upstream kinase analysis of kinome array data.

Phosphorylation by serine-threonine and tyrosine kinases is critical for determining protein function. Array-based platforms for measuring reporter peptide signal levels allow for differential phosphorylation analysis between conditions for distinct active kinases. Peptide array technologies like the PamStation12 from PamGene allow for generating high-throughput, multi-dimensional, and complex functional proteomics data. As the adoption rate of such technologies increases, there is an imperative need for software tools that streamline the process of analyzing such data. We present Kinome Random Sampling Analyzer (KRSA), an R package and R Shiny web-application for analyzing kinome array data to help users better understand the patterns of functional proteomics in complex biological systems. KRSA is an All-In-One tool that reads, formats, fits models, analyzes, and visualizes PamStation12 kinome data. While the underlying algorithm has been experimentally validated in previous publications, we demonstrate KRSA workflow on dorsolateral prefrontal cortex (DLPFC) in male (n = 3) and female (n = 3) subjects to identify differential phosphorylation signatures and upstream kinase activity. Kinase activity differences between males and females were compared to a previously published kinome dataset (11 female and 7 male subjects) which showed similar global phosphorylation signals patterns.

Phosphoproteomic Comparison of Four Eimeria tenella Life Cycle Stages.

Protein phosphorylation is an important post-translational modification (PTM) involved in diverse cellular functions. It is the most prevalent PTM in both Toxoplasma gondii and Plasmodium falciparum, but its status in Eimeria tenella has not been reported. Herein, we performed a comprehensive, quantitative phosphoproteomic profile analysis of four stages of the E. tenella life cycle: unsporulated oocysts (USO), partially sporulated (7 h) oocysts (SO7h), sporulated oocysts (SO), and sporozoites (S). A total of 15,247 phosphorylation sites on 9514 phosphopeptides corresponding to 2897 phosphoproteins were identified across the four stages. In addition, 456, 479, and 198 differentially expressed phosphoproteins (DEPPs) were identified in the comparisons SO7h vs. USO, SO vs. SO7h, and S vs. SO, respectively. Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEPPs suggested that they were involved in diverse functions. For SO7h vs. USO, DEPPs were mainly involved in cell division, actin cytoskeleton organization, positive regulation of transport, and pyruvate metabolism. For SO vs. SO7h, they were related to the peptide metabolic process, translation, and RNA transport. DEPPs in the S vs. SO comparison were associated with the tricarboxylic acid metabolic process, positive regulation of ATPase activity, and calcium ion binding. Time course sequencing data analysis (TCseq) identified six clusters with similar expression change characteristics related to carbohydrate metabolism, cytoskeleton organization, and calcium ion transport, demonstrating different regulatory profiles across the life cycle of E. tenella. The results revealed significant changes in the abundance of phosphoproteins during E. tenella development. The findings shed light on the key roles of protein phosphorylation and dephosphorylation in the E. tenella life cycle.

Respiratory Syncytial Virus Phosphoprotein Residue S156 Plays a Role in Regulating Genome Transcription and Replication.

Respiratory syncytial virus (RSV) is a single-stranded, negative-sense RNA virus in the family Pneumoviridae and genus Orthopneumovirus that can cause severe disease in infants, immunocompromised adults, and the elderly. The RSV viral RNA-dependent RNA polymerase (vRdRp) complex is composed of the phosphoprotein (P) and the large polymerase protein (L). The P protein is constitutively phosphorylated by host kinases and has 41 serine (S) and threonine (T) residues as potential phosphorylation sites. To identify important phosphorylation residues in the P protein, we systematically and individually mutated all S and T residues to alanine (A) and analyzed their effects on genome transcription and replication by using a minigenome system. We found that the mutation of eight residues resulted in minigenome activity significantly lower than that of wild-type (WT) P. We then incorporated these mutations (T210A, S203A, T151A, S156A, T160A, S23A, T188A, and T105A) into full-length genome cDNA to rescue recombinant RSV. We were able to recover four recombinant viruses (with T151A, S156A, T160A, or S23A), suggesting that RSV-P residues T210, S203, T188, and T105 are essential for viral RNA replication. Among the four recombinant viruses rescued, rRSV-T160A caused a minor growth defect relative to its parental virus while rRSV-S156A had severely restricted replication due to decreased levels of genomic RNA. During infection, P-S156A phosphorylation was decreased, and when passaged, the S156A virus acquired a known compensatory mutation in L (L795I) that enhanced both WT-P and P-S156A minigenome activity and was able to partially rescue the S156A viral growth defect. This work demonstrates that residues T210, S203, T188, and T105 are critical for RSV replication and that S156 plays a critical role in viral RNA synthesis. IMPORTANCE RSV-P is a heavily phosphorylated protein that is required for RSV replication. In this study, we identified several residues, including P-S156, as phosphorylation sites that play critical roles in efficient viral growth and genome replication. Future studies to identify the specific kinase(s) that phosphorylates these residues can lead to kinase inhibitors and antiviral drugs for this important human pathogen.

The Arabidopsis Root Tip (Phospho)Proteomes at Growth-Promoting versus Growth-Repressing Conditions Reveal Novel Root Growth Regulators.

Auxin plays a dual role in growth regulation and, depending on the tissue and concentration of the hormone, it can either promote or inhibit division and expansion processes in plants. Recent studies have revealed that, beyond transcriptional reprogramming, alternative auxin-controlled mechanisms regulate root growth. Here, we explored the impact of different concentrations of the synthetic auxin NAA that establish growth-promoting and -repressing conditions on the root tip proteome and phosphoproteome, generating a unique resource. From the phosphoproteome data, we pinpointed (novel) growth regulators, such as the RALF34-THE1 module. Our results, together with previously published studies, suggest that auxin, H+-ATPases, cell wall modifications and cell wall sensing receptor-like kinases are tightly embedded in a pathway regulating cell elongation. Furthermore, our study assigned a novel role to MKK2 as a regulator of primary root growth and a (potential) regulator of auxin biosynthesis and signalling, and suggests the importance of the MKK2 Thr31 phosphorylation site for growth regulation in the Arabidopsis root tip.

HLA molecule expression on the surface of cells and microparticles in platelet concentrates.

BACKGROUND: Platelet (PLT) transfusions are an essential treatment for bleeding disorders. However, immunologic complications can occur, including alloantibody production against Class I HLA molecules. The principal source of HLA molecules in PLT concentrates (PCs) is the PLTs themselves. However, extracellular microparticles (MPs) present in PCs may express HLA molecules. STUDY DESIGN AND METHODS: We used nanoscale flow cytometry to explore the expression of HLA-A2, HLA-B7, and HLA-B57 on the surface of cells, PLT-derived MPs (PMPs), lymphocyte-derived MPs (LMPs), and monocyte-derived MPs (MMPs) present in PCs. Expression was studied during 7 days of storage. RESULTS: Platelets were not the only source of HLA molecules in PCs. HLA molecules were present on PMPs, LMPs, and MMPs. The level of HLA Class I molecule expression varied between haplotypes and MPs of different origins and during storage. CONCLUSION: Platelets or residual cells remaining after leukoreduction are not the only source of HLA Class I molecules in PCs, highlighting the contribution of MPs to alloimmunization mechanisms. These data may be relevant for the development of new transfusion guidelines.

Characterization of Phosphorylated Proteins Using Mass Spectrometry.

Phosphorylation is arguably the most important post-translational modification that occurs within proteins. Phosphorylation is used as a signal to control numerous physiological activities ranging from gene expression to metabolism. Identifying phosphorylation sites within proteins was historically a challenge as it required either radioisotope labeling or the use of phospho-specific antibodies. The advent of mass spectrometry (MS) has had a major impact on the ability to qualitatively and quantitatively characterize phosphorylated proteins. In this article, we describe MS methods for characterizing phosphorylation sites within individual proteins as well as entire proteome samples. The utility of these methods is illustrated in examples that show the information that can be gained using these MS techniques.

Integrated proteomics and phosphoproteomics reveal perturbed regulative pathways in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) has a dismal prognosis largely owing to its inefficient diagnosis, rapid progress, and tenacious drug resistance. Here, we aimed to analyze the expressive patterns of proteins and phosphorylation in PDAC tissue samples and compare them to normal pancreatic tissue to investigate the underlying mechanisms and to reveal potential protein targets for diagnosis and drug development. Liquid chromatography coupled to mass spectrometry (LC-MS) based proteomics and phosphoproteomics analyses were performed using 20 pairs of patient-derived PDAC tissue and normal pancreatic tissue samples. Protein identification and quantification were conducted using MaxQuant software. Bioinformatics analysis was used to retrieve PDAC-relevant pathways and gene ontology (GO) terms. 4985 proteins and 3643 phosphoproteins were identified with high confidence; of these, 322 proteins and 235 phosphoproteins were dysregulated in PDAC. Several pathways, including several extracellular matrix-related pathways, were found to be strongly associated with PDAC. Further, the expression levels of filamin A (FLNA), integrin alpha-V (ITGAV), thymidine phosphorylase (TYMP), medium-chain specific acyl-CoA dehydrogenase, mitochondrial (ACADM), short-chain specific acyl-CoA dehydrogenase, mitochondrial (ACADS), and acetyl-CoA acetyltransferase, mitochondrial (ACAT1) were examined through western blot and immunohistochemistry analysis, and the results confirmed the validity of the proteomics analysis results. These findings provide comprehensive insight into the dysregulated regulative networks in PDAC tissue samples at the protein and phosphorylation levels, and they provide clues for further pathological studies and drug-target development.

Neuromedin U (NMU) regulates osteoblast differentiation and activity.

Osteoporosis is a disease of low bone mass that places individuals at enhanced risk for fracture, disability, and death. Osteoporosis rates are expected to rise significantly in the coming decades yet there are limited pharmacological treatment options, particularly for long-term management of this chronic condition. The drug development pipeline is relatively bereft of new strategies, causing an urgent and unmet need for developing new strategies and targets for treating osteoporosis. Here, we examine a lesser-studied bone remodeling pathway, Neuromedin U (NMU), which is expressed in the bone microenvironment along with its cognate receptors NMU receptor 1 (NMUR1) and 2 (NMUR2). We independently corroborate a prior report that global loss of NMU expression leads to high bone mass and test the hypothesis that NMU negatively regulates osteoblast differentiation. Consistent with this, in vitro studies reveal NMU represses osteoblastic differentiation of osteogenic precursors but, in contrast, promotes osteoblastic marker expression, proliferation and activity of osteoblast-like cells. Phospho-profiling arrays were used to detail differential signaling outcomes that may underlie the opposite responses of these cell types. Collectively, our findings indicate that NMU exerts cell-type-specific responses to regulate osteoblast differentiation and activity.

Enrichment of ATP Binding Proteins Unveils Proteomic Alterations in Human Macrophage Cell Death, Inflammatory Response, and Protein Synthesis after Interaction with Candida albicans.

Macrophages are involved in the primary human response to Candida albicans. After pathogen recognition, signaling pathways are activated, leading to the production of cytokines, chemokines, and antimicrobial peptides. ATP binding proteins are crucial for this regulation. Here, a quantitative proteomic and phosphoproteomic approach was carried out for the study of human macrophage ATP-binding proteins after interaction with C. albicans. From a total of 547 nonredundant quantified proteins, 137 were ATP binding proteins and 59 were detected as differentially abundant. From the differentially abundant ATP-binding proteins, 6 were kinases (MAP2K2, SYK, STK3, MAP3K2, NDKA, and SRPK1), most of them involved in signaling pathways. Furthermore, 85 phosphopeptides were quantified. Macrophage proteomic alterations including an increase of protein synthesis with a consistent decrease in proteolysis were observed. Besides, macrophages showed changes in proteins of endosomal trafficking together with mitochondrial proteins, including some involved in the response to oxidative stress. Regarding cell death mechanisms, an increase of antiapoptotic over pro-apoptotic signals is suggested. Furthermore, a high pro-inflammatory response was detected, together with no upregulation of key mi-RNAs involved in the negative feedback of this response. These findings illustrate a strategy to deepen the knowledge of the complex interactions between the host and the clinically important pathogen C. albicans.

Comparative Analysis of the Transcriptome and Proteome during Mouse Placental Development.

The condition of the placenta is a determinant of the short- and long-term health of the mother and the fetus. However, critical processes occurring in early placental development, such as trophoblast invasion and establishment of placental metabolism, remain poorly understood. To gain a better understanding of the genes involved in regulating these processes, we utilized a multiomics approach, incorporating transcriptome, proteome, and phosphoproteome data generated from mouse placental tissue collected at two critical developmental time points. We found that incorporating information from both the transcriptome and proteome identifies genes associated with time point-specific biological processes, unlike using the proteome alone. We further inferred genes upregulated on the basis of the proteome data but not the transcriptome data at each time point, leading us to identify 27 genes that we predict to have a role in trophoblast migration or placental metabolism. Finally, using the phosphoproteome data set, we discovered novel phosphosites that may play crucial roles in the regulation of placental transcription factors. By generating the largest proteome and phosphoproteome data sets in the developing placenta, and integrating transcriptome analysis, we uncovered novel aspects of placental gene regulation.

Phosphodiesterase beta is the master regulator of cAMP signalling during malaria parasite invasion.

Cyclic nucleotide signalling is a major regulator of malaria parasite differentiation. Phosphodiesterase (PDE) enzymes are known to control cyclic GMP (cGMP) levels in the parasite, but the mechanisms by which cyclic AMP (cAMP) is regulated remain enigmatic. Here, we demonstrate that Plasmodium falciparum phosphodiesterase ß (PDEß) hydrolyses both cAMP and cGMP and is essential for blood stage viability. Conditional gene disruption causes a profound reduction in invasion of erythrocytes and rapid death of those merozoites that invade. We show that this dual phenotype results from elevated cAMP levels and hyperactivation of the cAMP-dependent protein kinase (PKA). Phosphoproteomic analysis of PDEß-null parasites reveals a >2-fold increase in phosphorylation at over 200 phosphosites, more than half of which conform to a PKA substrate consensus sequence. We conclude that PDEß plays a critical role in governing correct temporal activation of PKA required for erythrocyte invasion, whilst suppressing untimely PKA activation during early intra-erythrocytic development.

Fission Yeast NDR/LATS Kinase Orb6 Regulates Exocytosis via Phosphorylation of the Exocyst Complex.

NDR/LATS kinases regulate multiple aspects of cell polarity and morphogenesis from yeast to mammals. Fission yeast NDR/LATS kinase Orb6 has been proposed to control cell polarity by regulating the Cdc42 guanine nucleotide exchange factor Gef1. Here, we show that Orb6 regulates polarity largely independently of Gef1 and that Orb6 positively regulates exocytosis. Through Orb6 inhibition in vivo and quantitative global phosphoproteomics, we identify Orb6 targets, including proteins involved in membrane trafficking. We confirm Sec3 and Sec5, conserved components of the exocyst complex, as substrates of Orb6 both in vivo and in vitro, and we show that Orb6 kinase activity is important for exocyst localization to cell tips and for exocyst activity during septum dissolution after cytokinesis. We further find that Orb6 phosphorylation of Sec3 contributes to exocyst function in concert with exocyst protein Exo70. We propose that Orb6 contributes to polarized growth by regulating membrane trafficking at multiple levels.

Highly Efficient Phosphoproteome Capture and Analysis from Urinary Extracellular Vesicles.

Analysis of protein phosphorylation in extracellular vesicles (EVs) offers an unprecedented potential for understanding cancer signaling and early stage disease diagnosis. However, prior to the phosphoproteome analysis step, the isolation of EVs from biofluids remains a challenging issue to overcome due to the low yield and impurity from current isolation methods. Here, we carry out an extensive assessment of several EV isolation methods including a novel rapid isolation method EVTRAP for highly efficient capture of extracellular vesicles from human urine sample. We demonstrate that over 95% recovery yield can be consistently achieved by EVTRAP, a significant improvement over current standard techniques. We then applied EVTRAP to identify over 16 000 unique peptides representing 2000 unique EV proteins from 200 µL urine sample, including all known EV markers with substantially increased recovery levels over ultracentrifugation. Most importantly, close to 2000 unique phosphopeptides were identified from more than 860 unique phosphoproteins using 10 mL of urine. The data demonstrated that EVTRAP is a highly effective and potentially widely implementable clinical isolation method for analysis of EV protein phosphorylation.

Exploring intrinsically disordered proteins in Chlamydomonas reinhardtii.

The content of intrinsically disordered protein (IDP) is related to organism complexity, evolution, and regulation. In the Plantae, despite their high complexity, experimental investigation of IDP content is lacking. We identified by mass spectrometry 682 heat-resistant proteins from the green alga, Chlamydomonas reinhardtii. Using a phosphoproteome database, we found that 331 of these proteins are targets of phosphorylation. We analyzed the flexibility propensity of the heat-resistant proteins and their specific features as well as those of predicted IDPs from the same organism. Their mean percentage of disorder was about 20%. Most of the IDPs (~70%) were addressed to other compartments than mitochondrion and chloroplast. Their amino acid composition was biased compared to other classic IDPs. Their molecular functions were diverse; the predominant ones were nucleic acid binding and unfolded protein binding and the less abundant one was catalytic activity. The most represented proteins were ribosomal proteins, proteins associated to flagella, chaperones and histones. We also found CP12, the only experimental IDP from C. reinhardtii that is referenced in disordered protein database. This is the first experimental investigation of IDPs in C. reinhardtii that also combines in silico analysis.

Comprehensive characterization of the adult ND4 Swiss Webster mouse retina: Using discovery-based mass spectrometry to decipher the total proteome and phosphoproteome.

Purpose: Diverse groups of proteins play integral roles in both the physiology and pathophysiology of the retina. However, thorough proteomic analyses of retinas of experimental species are currently unavailable. The purpose of the present paper is providing the field with a comprehensive proteomic characterization of the retina of a commonly used laboratory mouse using a discovery-based mass spectrometry (MS) approach. Methods: Retinas from eight male and eight female 30-week-old outbred ND4 Swiss Webster mice were harvested and immediately processed for MS analysis on a Thermo Fisher (TF) Fusion Orbitrap MS. The retinal proteome and phosphoproteome were identified and subsequently analyzed using Proteome Discoverer 2.2 and Panther-GeneGo. SEQUEST-HT scoring was used for analysis, and the reference protein FASTA database was from Mus musculus. Specifically, three technical repeats were performed for each biological sample. For characterization, only high-scoring peptides were considered, with a false discovery rate (FDR) of <1%. Downstream bioinformatic analysis used Ingenuity Pathway Analysis (IPA; Qiagen). Results: Using Proteome Discoverer 2.2, 4,767 different proteins were identified and segregated into 26 major protein classes, 9 functional molecular classes, and 12 categories of biological processes. The five largest protein classes included the following: nucleic acid binding (17%), hydrolases (13%), enzyme modulators (10%), transferases (9%), and oxidoreductases (6%). "Binding" and "catalytic" proteins contributed to 81% of the molecular function class at 37% and 42%, respectively. "Cellular processing" and "metabolic processes" contributed the most to biologic activity, at 31% and 26%, respectively. Phosphopeptide enrichment yielded the identification of 610 additional unique proteins that were not originally identified. The two datasets combined produced an adult mouse retinal proteome consisting of 5,377 unique proteins. Overall, 41% of the retinal proteome was phosphorylated. The overwhelming diversity of retinal protein functionality was reflected through further analyses revealing 2,086 unique pathway hits across 241 different pathways (TF). A core analysis summary report was performed in IPA (Qiagen) to analyze the top signaling networks, protein-protein interaction (PPI) enrichments, and canonical pathways. Conclusions: Using this high-throughput technique, we have further deciphered and updated the diverse proteome of the mouse retina, including the phosphoproteome, thereby providing the most comprehensive proteomic profile for this tissue known to date. These findings, and the bioinformatic analyses we also provided, establish a platform for future studies, facilitating the elucidation of the relevance of these proteins to the molecular and cellular pathologies that underlie retinal function and disease.

Trans-biome diversity in Australian grass-specialist lizards (Diplodactylidae: Strophurus).

Comparisons of biodiversity patterns within lineages that occur across major climate gradients and biomes, can provide insights into the relative roles that lineage history, landscape and climatic variation, and environmental change have played in shaping regional biotas. In Australia, while there has been extensive research into the origins and patterns of diversity in the Australian Arid Zone (AAZ), how diversity is distributed across this biome and the Australian Monsoonal Tropics (AMT) to the north, has been less studied. We compared the timing and patterns of diversification across this broad aridity gradient in a clade of lizards (Strophurus: phasmid geckos) that only occur in association with a unique Australian radiation of sclerophyllous grasses (Triodia: spinifex). Our results indicate that overall genetic diversity is much higher, older and more finely geographically structured within the AMT, including distantly related clades endemic to the sandstone escarpments of the Kimberley and Arnhem Plateau. Niche modelling analyses also suggest that the distribution of taxa in the AMT is more strongly correlated with variation in topographic relief than in the AAZ. The two broad patterns that we recovered - (i) lineage endemism increases as latitude decreases, and (ii) endemism is tightly correlated to rocky regions - parallel and corroborate other recent studies of habitat generalists and specialised saxicoline lineages occurring across these same regions. Early Miocene diversification estimates also suggest that, soon after Triodia grasses colonised Australia and began to diversify in the Miocene, phasmid geckos with Gondwanan ancestry shifted into these grasses, and have subsequently remained closely associated with this unique vegetation type.

iPhos-PseEvo: Identifying Human Phosphorylated Proteins by Incorporating Evolutionary Information into General PseAAC via Grey System Theory.

Protein phosphorylation plays a critical role in human body by altering the structural conformation of a protein, causing it to become activated/deactivated, or functional modification. Given an uncharacterized protein sequence, can we predict whether it may be phosphorylated or may not? This is no doubt a very meaningful problem for both basic research and drug development. Unfortunately, to our best knowledge, so far no high throughput bioinformatics tool whatsoever has been developed to address such a very basic but important problem due to its extremely complexity and lacking sufficient training data. Here we proposed a predictor called iPhos-PseEvo by (1) incorporating the protein sequence evolutionary information into the general pseudo amino acid composition (PseAAC) via the grey system theory, (2) balancing out the skewed training datasets by the asymmetric bootstrap approach, and (3) constructing an ensemble predictor by fusing an array of individual random forest classifiers thru a voting system. Rigorous jackknife tests have indicated that very promising success rates have been achieved by iPhos-PseEvo even for such a difficult problem. A user-friendly web-server for iPhos-PseEvo has been established at http://www.jci-bioinfo.cn/iPhos-PseEvo, by which users can easily obtain their desired results without the need to go through the complicated mathematical equations involved. It has not escaped our notice that the formulation and approach presented here can be used to analyze many other problems in protein science as well.

Differential Phosphoproteomic Analysis of Recombinant Chinese Hamster Ovary Cells Following Temperature Shift.

Phosphorylation is one of the most important post-translational modifications, playing a crucial role in regulating many cellular processes, including transcription, cytoskeletal rearrangement, cell proliferation, differentiation, apoptosis, and signal transduction. However, to date, little work has been carried out on the phosphoproteome in CHO cells. In this study we have carried out a large scale differential phosphoproteomic analysis of recombinant CHO cells following a reduction of culture temperature (temperature shift). The reduction of culture temperature during the exponential phase of growth is commonly employed by the biopharmaceutical industry to increase product yield; however, the molecular mechanisms of temperature shift in CHO cells remain poorly understood. We have identified 700 differentially expressed phosphopeptides using quantitative label-free LC-MS/MS phosphoproteomic analysis in conjunction with IMAC and TiO2 phosphopeptide enrichment strategies, following a reduction in temperature from 37 to 31 °C. Functional assessment of the phosphoproteomic data using gene ontology analysis showed a significant enrichment of biological processes related to growth (e.g., cell cycle, cell division), ribosomal biogenesis, and cytoskeleton organization, and molecular functions related to RNA binding, transcription factor activity, and protein serine/threonine kinase activity. Differential phosphorylation of two proteins, ATF2 and NDRG1, was confirmed by Western blotting. This data suggests the importance of including the post-translational layer of regulation, such as phosphorylation, in CHO "omics" studies. This study also has the potential to identify phosphoprotein targets that could be modified using cell line engineering approaches to improve the efficiency of recombinant protein production.

Glucose-regulated and drug-perturbed phosphoproteome reveals molecular mechanisms controlling insulin secretion.

Insulin-secreting beta cells play an essential role in maintaining physiological blood glucose levels, and their dysfunction leads to the development of diabetes. To elucidate the signalling events regulating insulin secretion, we applied a recently developed phosphoproteomics workflow. We quantified the time-resolved phosphoproteome of murine pancreatic cells following their exposure to glucose and in combination with small molecule compounds that promote insulin secretion. The quantitative phosphoproteome of 30,000 sites clustered into three main groups in concordance with the modulation of the three key kinases: PKA, PKC and CK2A. A high-resolution time course revealed key novel regulatory sites, revealing the importance of methyltransferase DNMT3A phosphorylation in the glucose response. Remarkably a significant proportion of these novel regulatory sites is significantly downregulated in diabetic islets. Control of insulin secretion is embedded in an unexpectedly broad and complex range of cellular functions, which are perturbed by drugs in multiple ways.