GPS-ARM: computational analysis of the APC/C recognition motif by predicting D-boxes and KEN-boxes

A Review of Machine Learning and Algorithmic Methods for Protein Phosphorylation Site Prediction

Post-translational modifications (PTMs) have key roles in extending the functional diversity of proteins and, as a result, regulating diverse cellular processes in prokaryotic and eukaryotic organisms. Phosphorylation modification is a vital PTM that occurs in most proteins and plays a significant role in many biological processes. Disorders in the phosphorylation process lead to multiple diseases, including neurological disorders and cancers. The purpose of this review is to organize this body of knowledge associated with phosphorylation site (p-site) prediction to facilitate future research in this field. At first, we comprehensively review all related databases and introduce all steps regarding dataset creation, data preprocessing, and method evaluation in p-site prediction. Next, we investigate p-site prediction methods, which are divided into two computational groups: algorithmic and machine learning (ML). Additionally, it is shown that there are basically two main approaches for p-site prediction by ML: conventional and end-to-end deep learning methods, both of which are given an overview. Moreover, this review introduces the most important feature extraction techniques, which have mostly been used in p-site prediction. Finally, we create three test sets from new proteins related to the released version of the database of protein post-translational modifications (dbPTM) in 2022 based on general and human species. Evaluating online p-site prediction tools on newly added proteins introduced in the dbPTM 2022 release, distinct from those in the dbPTM 2019 release, reveals their limitations. In other words, the actual performance of these online p-site prediction tools on unseen proteins is notably lower than the results reported in their respective research papers.

Counting Degrons: Lessons From Multivalent Substrates for Targeted Protein Degradation

E3s comprise a structurally diverse group of at least 800 members, most of which target multiple substrates through specific and regulated protein-protein interactions. These interactions typically rely on short linear motifs (SLiMs), called "degrons", in an intrinsically disordered region (IDR) of the substrate, with variable rules of engagement governing different E3-docking events. These rules of engagement are of importance to the field of targeted protein degradation (TPD), where substrate ubiquitination and destruction require tools to effectively harness ubiquitin ligases (E3s). Substrates are often found to contain multiple degrons, or multiple copies of a degron, contributing to the affinity and selectivity of the substrate for its E3. One important paradigm for E3-substrate docking is presented by the Anaphase-Promoting Complex/Cyclosome (APC/C), a multi-subunit E3 ligase that targets hundreds of proteins for destruction during mitotic exit. APC/C substrate targeting takes place in an ordered manner thought to depend on tightly regulated interactions of substrates, with docking sites provided by the substoichiometric APC/C substrate adaptors and coactivators, Cdc20 or Cdh1/FZR1. Both structural and functional studies of individual APC/C substrates indicate that productive ubiquitination usually requires more than one degron, and that degrons are of different types docking to distinct sites on the coactivators. However, the dynamic nature of APC/C substrate recruitment, and the influence of multiple degrons, remains poorly understood. Here we review the significance of multiple degrons in a number of E3-substrate interactions that have been studied in detail, illustrating distinct kinetic effects of multivalency and allovalency, before addressing the role of multiple degrons in APC/C substrates, key to understanding ordered substrate destruction by APC/C. Lastly, we consider how lessons learnt from these studies can be applied in the design of TPD tools.

Most non-canonical proteins uniquely populate the proteome or immunopeptidome

Combining RNA sequencing, ribosome profiling, and mass spectrometry, we elucidate the contribution of non-canonical translation to the proteome and major histocompatibility complex (MHC) class I immunopeptidome. Remarkably, of 14,498 proteins identified in three human B cell lymphomas, 2,503 are non-canonical proteins. Of these, 28% are novel isoforms and 72% are cryptic proteins encoded by ostensibly non-coding regions (60%) or frameshifted canonical genes (12%). Cryptic proteins are translated as efficiently as canonical proteins, have more predicted disordered residues and lower stability, and critically generate MHC-I peptides 5-fold more efficiently per translation event. Translating 5' "untranslated" regions hinders downstream translation of genes involved in transcription, translation, and antiviral responses. Novel protein isoforms show strong enrichment for signaling pathways deregulated in cancer. Only a small fraction of cryptic proteins detected in the proteome contribute to the MHC-I immunopeptidome, demonstrating the high preferential access of cryptic defective ribosomal products to the class I pathway.

pCysMod: Prediction of Multiple Cysteine Modifications Based on Deep Learning Framework

Thiol groups on cysteines can undergo multiple post-translational modifications (PTMs), acting as a molecular switch to maintain redox homeostasis and regulating a series of cell signaling transductions. Identification of sophistical protein cysteine modifications is crucial for dissecting its underlying regulatory mechanism. Instead of a time-consuming and labor-intensive experimental method, various computational methods have attracted intense research interest due to their convenience and low cost. Here, we developed the first comprehensive deep learning based tool pCysMod for multiple protein cysteine modification prediction, including S-nitrosylation, S-palmitoylation, S-sulfenylation, S-sulfhydration, and S-sulfinylation. Experimentally verified cysteine sites curated from literature and sites collected by other databases and predicting tools were integrated as benchmark dataset. Several protein sequence features were extracted and united into a deep learning model, and the hyperparameters were optimized by particle swarm optimization algorithms. Cross-validations indicated our model showed excellent robustness and outperformed existing tools, which was able to achieve an average AUC of 0.793, 0.807, 0.796, 0.793, and 0.876 for S-nitrosylation, S-palmitoylation, S-sulfenylation, S-sulfhydration, and S-sulfinylation, demonstrating pCysMod was stable and suitable for protein cysteine modification prediction. Besides, we constructed a comprehensive protein cysteine modification prediction web server based on this model to benefit the researches finding the potential modification sites of their interested proteins, which could be accessed at http://pcysmod.omicsbio.info. This work will undoubtedly greatly promote the study of protein cysteine modification and contribute to clarifying the biological regulation mechanisms of cysteine modification within and among the cells.

MusiteDeep: a deep-learning based webserver for protein post-translational modification site prediction and visualization

MusiteDeep is an online resource providing a deep-learning framework for protein post-translational modification (PTM) site prediction and visualization. The predictor only uses protein sequences as input and no complex features are needed, which results in a real-time prediction for a large number of proteins. It takes less than three minutes to predict for 1000 sequences per PTM type. The output is presented at the amino acid level for the user-selected PTM types. The framework has been benchmarked and has demonstrated competitive performance in PTM site predictions by other researchers. In this webserver, we updated the previous framework by utilizing more advanced ensemble techniques, and providing prediction and visualization for multiple PTMs simultaneously for users to analyze potential PTM cross-talks directly. Besides prediction, users can interactively review the predicted PTM sites in the context of known PTM annotations and protein 3D structures through homology-based search. In addition, the server maintains a local database providing pre-processed PTM annotations from Uniport/Swiss-Prot for users to download. This database will be updated every three months. The MusiteDeep server is available at https://www.musite.net. The stand-alone tools for locally using MusiteDeep are available at https://github.com/duolinwang/MusiteDeep_web.

Kaposi's Sarcoma-Associated Herpesvirus Lytic Replication Is Independent of Anaphase-Promoting Complex Activity

The anaphase-promoting complex, or cyclosome (APC/C), is a large E3 ubiquitin ligase composed of 14 subunits. The activity of APC/C oscillates during the cell cycle to ensure a timely transition through each phase by promoting the degradation of important cell cycle regulators. Of the human herpesviruses, cytomegalovirus (HCMV) and Epstein-Barr virus (EBV) both impair the activity of APC/C during their lytic replication cycle through virus-encoded protein kinases. Here, we addressed whether the oncogenic Kaposi's sarcoma-associated herpesvirus (KSHV) deregulates the activity of APC/C during the lytic replication cycle. To this end, we used the well-characterized iSLK.219 cell model of KSHV infection and established a new infection model of primary lymphatic endothelial cells (LECs) infected with a lytically replicating KSHV BAC16 mutant. In contrast to those of EBV and HCMV, the KSHV lytic cycle occurs while the APC/C is active. Moreover, interfering with the activity of APC/C did not lead to major changes in the production of infectious virus. We further investigated whether rereplication stress induced by the unscheduled activation of the APC/C-CDH1 complex affects the number and integrity of KSHV viral episomes. Deep sequencing of the viral episomes and host chromosomes in iSLK.219 cells revealed that, while distinct regions in the cellular chromosomes were severely affected by rereplication stress, the integrity of the viral episomes remained unaltered.IMPORTANCE DNA viruses have evolved complex strategies to gain control over the cell cycle. Several of them target APC/C, a key cellular machinery that controls the timely progression of the cell cycle, by either blocking or enhancing its activity. Here, we investigated the activity of APC/C during the lytic replication cycle of KSHV and found that, in contrast to that of KSHV's close relatives EBV and HCMV, KSHV lytic replication occurs while the APC/C is active. Perturbing APC/C activity by depleting a core protein or the adaptor proteins of the catalytic domain, and hence interfering with normal cell-cycle progression, did not affect virus replication. This suggests that KSHV has evolved to replicate independently of the activity of APC/C and in various cell cycle conditions.

GPS 5.0: An Update on the Prediction of Kinase-specific Phosphorylation Sites in Proteins

In eukaryotes, protein phosphorylation is specifically catalyzed by numerous protein kinases (PKs), faithfully orchestrates various biological processes, and reversibly determines cellular dynamics and plasticity. Here we report an updated algorithm of Group-based Prediction System (GPS) 5.0 to improve the performance for predicting kinase-specific phosphorylation sites (p-sites). Two novel methods, position weight determination (PWD) and scoring matrix optimization (SMO), were developed. Compared with other existing tools, GPS 5.0 exhibits a highly competitive accuracy. Besides serine/threonine or tyrosine kinases, GPS 5.0 also supports the prediction of dual-specificity kinase-specific p-sites. In the classical module of GPS 5.0, 617 individual predictors were constructed for predicting p-sites of 479 human PKs. To extend the application of GPS 5.0, a species-specific module was implemented to predict kinase-specific p-sites for 44,795 PKs in 161 eukaryotes. The online service and local packages of GPS 5.0 are freely available for academic research at http://gps.biocuckoo.cn.

Checkpoint Regulation of Nuclear Tos4 Defines S Phase Arrest in Fission Yeast

From yeast to humans, the cell cycle is tightly controlled by regulatory networks that regulate cell proliferation and can be monitored by dynamic visual markers in living cells. We have observed S phase progression by monitoring nuclear accumulation of the FHA-containing DNA binding protein Tos4, which is expressed in the G1/S phase transition. We use Tos4 localization to distinguish three classes of DNA replication mutants: those that arrest with an apparent 1C DNA content and accumulate Tos4 at the restrictive temperature; those that arrest with an apparent 2C DNA content, that do not accumulate Tos4; and those that proceed into mitosis despite a 1C DNA content, again without Tos4 accumulation. Our data indicate that Tos4 localization in these conditions is responsive to checkpoint kinases, with activation of the Cds1 checkpoint kinase promoting Tos4 retention in the nucleus, and activation of the Chk1 damage checkpoint promoting its turnover. Tos4 localization therefore allows us to monitor checkpoint-dependent activation that responds to replication failure in early vs. late S phase.

Precise Prediction of Calpain Cleavage Sites and Their Aberrance Caused by Mutations in Cancer

As a widespread post-translational modification of proteins, calpain-mediated cleavage regulates a broad range of cellular processes, including proliferation, differentiation, cytoskeletal reorganization, and apoptosis. The identification of proteins that undergo calpain cleavage in a site-specific manner is the necessary foundation for understanding the exact molecular mechanisms and regulatory roles of calpain-mediated cleavage. In contrast with time-consuming and labor-intensive experimental methods, computational approaches for detecting calpain cleavage sites have attracted wide attention due to their efficiency and convenience. In this study, we established a novel computational tool named DeepCalpain (http://deepcalpain.cancerbio.info/) for predicting the potential calpain cleavage sites by adopting deep neural network and the particle swarm optimization algorithm. Through critical evaluation and comparison, DeepCalpain exhibited superior performance against other existing tools. Meanwhile, we found that protein interactions could enrich the calpain-substrate regulatory relationship. Since calpain-mediated cleavage was critical for cancer development and progression, we comprehensively analyzed the calpain cleavage associated mutations across 11 cancers with the help of DeepCalpain, which demonstrated that the calpain-mediated cleavage events were affected by mutations and heavily implicated in the regulation of cancer cells. These prediction and analysis results might provide helpful information to reveal the regulatory mechanism of calpain cleavage in biological pathways and different cancer types, which might open new avenues for the diagnosis and treatment of cancers.

SWITCH 1/DYAD is a WINGS APART-LIKE antagonist that maintains sister chromatid cohesion in meiosis

Mitosis and meiosis both rely on cohesin, which embraces the sister chromatids and plays a crucial role for the faithful distribution of chromosomes to daughter cells. Prior to the cleavage by Separase at anaphase onset, cohesin is largely removed from chromosomes by the non-proteolytic action of WINGS APART-LIKE (WAPL), a mechanism referred to as the prophase pathway. To prevent the premature loss of sister chromatid cohesion, WAPL is inhibited in early mitosis by Sororin. However, Sororin homologs have only been found to function as WAPL inhibitors during mitosis in vertebrates and Drosophila. Here we show that SWITCH 1/DYAD defines a WAPL antagonist that acts in meiosis of Arabidopsis. Crucially, SWI1 becomes dispensable for sister chromatid cohesion in the absence of WAPL. Despite the lack of any sequence similarities, we found that SWI1 is regulated and functions in a similar manner as Sororin hence likely representing a case of convergent molecular evolution across the eukaryotic kingdom.

Differential Expression of Genes for Ubiquitin Ligases in Medulloblastoma Subtypes

Using publically available datasets on gene expression in medulloblastoma (MB) subtypes, we selected genes for ubiquitin ligases and identified statistically those that best predicted each of the four major MB subgroups as separate disease entities. We identify a gene coding for an ubiquitin ligase, ZNRF3, whose overexpression alone can predict the WNT subgroup for 100% in the Pfister dataset. For the SHH subgroup, we identify a gene for a regulatory subunit of the protein phosphatase 2A (PP2A), PPP2R2C, as the major predictor among the E3 ligases genes. The ubiquitin and ubiquitin-like conjugation database (UUCD) lists PPP2R2C as coding for a Cullin Ring ubiquitin ligase adaptor. For group 3 MBs, the best ubiquitin ligase predictor was PPP2R2B, a gene which codes for another regulatory subunit of the PP2A holoenzyme. For group 4, the best E3 gene predictors were MID2, ZBTB18, and PPP2R2A, which codes for a third PP2A regulatory subunit. Heatmap analysis of the E3 gene data shows that expression of ten genes for ubiquitin ligases can be used to classify MBs into the four major consensus subgroups. This was illustrated by analysis of gene expression of ubiquitin ligases of the Pfister dataset and confirmed in the dataset of Cavalli. We conclude that genes for ubiquitin ligases can be used as genetic markers for MB subtypes and that the proteins coded for by these genes should be investigated as subtype specific therapeutic targets for MB.

APC/CFzr regulates cardiac and myoblast cell numbers and plays a crucial role during myoblast fusion

Somatic muscles are formed by the iterative fusion of myoblasts into muscle fibres. This process is driven by the recurrent recruitment of proteins to the cell membrane to induce F-actin nucleation at the fusion site. Although various proteins involved in myoblast fusion have been identified, knowledge about their sub-cellular regulation is rather elusive. We identified the anaphase-promoting complex (APC/C) adaptor Fizzy related (Fzr) as an essential regulator of heart and muscle development. We show that APC/CFzr regulates the fusion of myoblasts as well as mitotic exit of pericardial cells, cardioblasts and myoblasts. Surprisingly, over-proliferation is not causative for the observed fusion defects. Instead, fzr mutants exhibit smaller F-actin foci at the fusion site, and display reduced membrane breakdown between adjacent myoblasts. We show that lack of APC/CFzr causes the accumulation and mislocalisation of Rols and Duf, two proteins involved in the fusion process. Duf seems to serve as direct substrate of the APC/CFzr, and its destruction depends on the presence of distinct degron sequences. These novel findings indicate that protein destruction and turnover constitute major events during myoblast fusion.

Proteomics of phosphorylation and protein dynamics during fertilization and meiotic exit in the Xenopus egg

Fertilization releases the meiotic arrest and initiates the events that prepare the egg for the ensuing developmental program. Protein degradation and phosphorylation are known to regulate protein activity during this process. However, the full extent of protein loss and phosphoregulation is still unknown. We examined absolute protein and phosphosite dynamics of the fertilization response by mass spectrometry-based proteomics in electroactivated eggs. To do this, we developed an approach for calculating the stoichiometry of phosphosites from multiplexed proteomics that is compatible with dynamic, stable, and multisite phosphorylation. Overall, the data suggest that degradation is limited to a few low-abundance proteins. However, this degradation promotes extensive dephosphorylation that occurs over a wide range of abundances during meiotic exit. We also show that eggs release a large amount of protein into the medium just after fertilization, most likely related to the blocks to polyspermy. Concomitantly, there is a substantial increase in phosphorylation likely tied to calcium-activated kinases. We identify putative degradation targets and components of the slow block to polyspermy. The analytical approaches demonstrated here are broadly applicable to studies of dynamic biological systems.

The preprophase band of microtubules controls the robustness of division orientation in plants

Controlling cell division plane orientation is essential for morphogenesis in multicellular organisms. In plant cells, the future cortical division plane is marked before mitotic entry by the preprophase band (PPB). Here, we characterized an Arabidopsis trm (TON1 Recruiting Motif) mutant that impairs PPB formation but does not affect interphase microtubules. Unexpectedly, PPB disruption neither abolished the capacity of root cells to define a cortical division zone nor induced aberrant cell division patterns but rather caused a loss of precision in cell division orientation. Our results advocate for a reassessment of PPB function and division plane determination in plants and show that a main output of this microtubule array is to limit spindle rotations in order to increase the robustness of cell division.

Determinants and Regulation of Protein Turnover in Yeast

Protein turnover maintains the recycling needs of the proteome, and its malfunction has been linked to aging and age-related diseases. However, not all proteins turnover equally, and the factors that contribute to accelerate or slow down turnover are mostly unknown. We measured turnover rates for 3,160 proteins in exponentially growing yeast and analyzed their dependence on physical, functional, and genetic properties. We found that functional characteristics, including protein localization, complex membership, and connectivity, have greater effect on turnover than sequence elements. We also found that protein turnover and mRNA turnover are correlated. Analysis under nutrient perturbation and osmotic stress revealed that protein turnover highly depends on cellular state and is faster when proteins are being actively used. Finally, stress-induced changes in protein and transcript abundance correlated with changes in protein turnover. This study provides a resource of protein turnover rates and principles to understand the recycling needs of the proteome under basal conditions and perturbation.

Association of Uba6-Specific-E2 (USE1) With Lung Tumorigenesis

Background

The UBA6-specific E2 conjugating enzyme 1 (USE1) ubiquitin enzyme cascade is a poorly characterized arm of the ubiquitin-proteasome system. We investigated whether the UBA6-USE1 enzyme cascade plays a role in lung cancer tumorigenesis.

Methods

USE1 expression was assessed in tumor-normal paired samples from 106 lung cancer patients by immunoblot. USE1 was stably overexpressed and knocked down in lung cancer cell lines to evaluate cell proliferation, colony formation, and invasion. Xenograft models were used to determine the effects of USE1 on tumor growth (n = 7). Proteomics analysis was used to identify proteins interacting with USE1. The USE1 gene was sequenced in lung cancer patients, and missense mutations of USE1 were generated to evaluate its function. All statistical tests were two-sided.

Results

USE1 proteins were frequently overexpressed in lung cancer patients (92.5%) Stable overexpression of USE1 increased cell proliferation ( P = .002), migration ( P < .001), and invasion ( P < .001), whereas knockdown of USE1 reduced cell proliferation ( P < .001), migration ( P = .003), and invasion in lung cancer cells and xenograft models ( P < .001). USE1 was found to have a conserved D-box domain, and the level of the protein was regulated by the anaphase-promoting complex. Several missense mutations in USE1 identified in patients prolong the stability of the protein.

Conclusions

USE1 proteins are frequently overexpressed in lung cancer, and missense mutations in USE1 prolong the half-life of the protein, promoting tumor formation. Our findings reveal novel roles for USE1 in lung cancer and the possible use of USE1 as a novel biomarker and therapeutic target for lung cancer treatment.

New Functions of APC/C Ubiquitin Ligase in the Nervous System and Its Role in Alzheimer's Disease

The E3 ubiquitin ligase Anaphase Promoting Complex/Cyclosome (APC/C) regulates important processes in cells, such as the cell cycle, by targeting a set of substrates for degradation. In the last decade, APC/C has been related to several major functions in the nervous system, including axon guidance, synaptic plasticity, neurogenesis, and neuronal survival. Interestingly, some of the identified APC/C substrates have been related to neurodegenerative diseases. There is an accumulation of some degradation targets of APC/C in Alzheimer’s disease (AD) brains, which suggests a dysregulation of the protein complex in the disorder. Moreover, recently evidence has been provided for an inactivation of APC/C in AD. It has been shown that oligomers of the AD-related peptide, Aβ, induce degradation of the APC/C activator subunit cdh1, in vitro in neurons in culture and in vivo in the mouse hippocampus. Furthermore, in the AD mouse model APP/PS1, lower cdh1 levels were observed in pyramidal neurons in CA1 when compared to age-matched wildtype mice. In this review, we provide a complete list of APC/C substrates that are involved in the nervous system and we discuss their functions. We also summarize recent studies that show neurobiological effects in cdh1 knockout mouse models. Finally, we discuss the role of APC/C in the pathophysiology of AD.

Different modes of APC/C activation control growth and neuron-glia interaction in the developing Drosophila eye

The development of the nervous system requires tight control of cell division, fate specification and migration. The anaphase promoting complex/cyclosome (APC/C) is an E3 ubiquitin ligase that affects different steps of cell cycle progression, as well as having postmitotic functions in nervous system development. It can therefore link different developmental stages in one tissue. The two adaptor proteins Fizzy/Cdc20 and Fizzy-Related/Cdh1 confer APC/C substrate specificity. Here we show that two distinct modes of APC/C function act during Drosophila eye development. Fizzy/Cdc20 controls the early growth of the eye disc anlage and the concomitant entry of glial cells onto the disc. In contrast, fzr/cdh1 acts during neuronal patterning and photoreceptor axon growth, and subsequently affects neuron-glia interaction. To further address the postmitotic role of Fzr/Cdh1 in controlling neuron-glia interaction, we identified a series of novel APC/C candidate substrates. Four of our candidate genes are required for fzr/cdh1 dependent neuron-glia interaction, including the dynein light chain Dlc90F. Taken together, our data show how different modes of APC/C activation can couple early growth and neuron-glia interaction during eye disc development.

GPS-PAIL: prediction of lysine acetyltransferase-specific modification sites from protein sequences

Protein acetylation catalyzed by specific histone acetyltransferases (HATs) is an essential post-translational modification (PTM) and involved in the regulation a broad spectrum of biological processes in eukaryotes. Although several ten thousands of acetylation sites have been experimentally identified, the upstream HATs for most of the sites are unclear. Thus, the identification of HAT-specific acetylation sites is fundamental for understanding the regulatory mechanisms of protein acetylation. In this work, we first collected 702 known HAT-specific acetylation sites of 205 proteins from the literature and public data resources, and a motif-based analysis demonstrated that different types of HATs exhibit similar but considerably distinct sequence preferences for substrate recognition. Using 544 human HAT-specific sites for training, we constructed a highly useful tool of GPS-PAIL for the prediction of HAT-specific sites for up to seven HATs, including CREBBP, EP300, HAT1, KAT2A, KAT2B, KAT5 and KAT8. The prediction accuracy of GPS-PAIL was critically evaluated, with a satisfying performance. Using GPS-PAIL, we also performed a large-scale prediction of potential HATs for known acetylation sites identified from high-throughput experiments in nine eukaryotes. Both online service and local packages were implemented, and GPS-PAIL is freely available at: http://pail.biocuckoo.org.

MHC class I-associated peptides derive from selective regions of the human genome

MHC class I-associated peptides (MAPs) define the immune self for CD8+ T lymphocytes and are key targets of cancer immunosurveillance. Here, the goals of our work were to determine whether the entire set of protein-coding genes could generate MAPs and whether specific features influence the ability of discrete genes to generate MAPs. Using proteogenomics, we have identified 25,270 MAPs isolated from the B lymphocytes of 18 individuals who collectively expressed 27 high-frequency HLA-A,B allotypes. The entire MAP repertoire presented by these 27 allotypes covered only 10% of the exomic sequences expressed in B lymphocytes. Indeed, 41% of expressed protein-coding genes generated no MAPs, while 59% of genes generated up to 64 MAPs, often derived from adjacent regions and presented by different allotypes. We next identified several features of transcripts and proteins associated with efficient MAP production. From these data, we built a logistic regression model that predicts with good accuracy whether a gene generates MAPs. Our results show preferential selection of MAPs from a limited repertoire of proteins with distinctive features. The notion that the MHC class I immunopeptidome presents only a small fraction of the protein-coding genome for monitoring by the immune system has profound implications in autoimmunity and cancer immunology.

RNA helicase Spn-E is required to maintain Aub and AGO3 protein levels for piRNA silencing in the germline of Drosophila

Germline-specific RNA helicase Spindle-E (Spn-E) is known to be essential for piRNA silencing in Drosophila that takes place mainly in the perinuclear nuage granules. Loss-of-function spn-E mutations lead to tandem Stellate genes derepression in the testes and retrotransposon mobilization in the ovaries. However, Spn-E functions in the piRNA pathway are still obscure. Analysis of total library of short RNAs from the testes of spn-E heterozygous flies revealed the presence of abundant piRNA ping-pong pairs originating from Su(Ste) transcripts. The abundance of these ping-pong pairs were sharply reduced in the library from the testes of spn-E mutants. Thus we found that ping-pong mechanism contributed to Su(Ste) piRNA generation in the testes. The lack of Spn-E caused a significant drop of protein levels of key ping-pong participants, Aubergine (Aub) and AGO3 proteins of PIWI subfamily, in the germline of both males and females, but did not disrupt of their assembly in nuage granules. We found that observed decline of the protein expression was not caused by suppression of aub and ago3 transcription as well as total transcription, indicating possible contribution of Spn-E to post-transcriptional regulation.

A genome-scale CRISPR-Cas9 screening method for protein stability reveals novel regulators of Cdc25A

The regulation of stability is particularly crucial for unstable proteins in cells. However, a convenient and unbiased method of identifying regulators of protein stability remains to be developed. Recently, a genome-scale CRISPR-Cas9 library has been established as a genetic tool to mediate loss-of-function screening. Here, we developed a protein stability regulators screening assay (Pro-SRSA) by combining the whole-genome CRISPR-Cas9 library with a dual-fluorescence-based protein stability reporter and high-throughput sequencing to screen for regulators of protein stability. Using Cdc25A as an example, Cul4B-DDB1^DCAF8 was identified as a new E3 ligase for Cdc25A. Moreover, the acetylation of Cdc25A at lysine 150, which was acetylated by p300/CBP and deacetylated by HDAC3, prevented the ubiquitin-mediated degradation of Cdc25A by the proteasome. This is the first study to report that acetylation, as a novel posttranslational modification, modulates Cdc25A stability, and we suggest that this unbiased CRISPR-Cas9 screening method at the genome scale may be widely used to globally identify regulators of protein stability.

Degradation of Ndd1 by APC/C(Cdh1) generates a feed forward loop that times mitotic protein accumulation

Ndd1 activates the Mcm1-Fkh2 transcription factor to transcribe mitotic regulators. The anaphase-promoting complex/cyclosome activated by Cdh1 (APC/C(Cdh1)) mediates the degradation of proteins throughout G1. Here we show that the APC/C(Cdh1) ubiquitinates Ndd1 and mediates its degradation, and that APC/C(Cdh1) activity suppresses accumulation of Ndd1 targets. We confirm putative Ndd1 targets and identify novel ones, many of them APC/C(Cdh1) substrates. The APC/C(Cdh1) thus regulates these proteins in a dual manner—both pretranscriptionally and post-translationally, forming a multi-layered feedforward loop (FFL). We predict by mathematical modelling and verify experimentally that this FFL introduces a lag between APC/C(Cdh1) inactivation at the end of G1 and accumulation of genes transcribed by Ndd1 in G2. This regulation generates two classes of APC/C(Cdh1) substrates, early ones that accumulate in S and late ones that accumulate in G2. Our results show how the dual state APC/C(Cdh1) activity is converted into multiple outputs by interactions between its substrates.

Specificity of the anaphase-promoting complex: a single-molecule study

Biological processes require specific enzymatic reactions, paradoxically involving short recognition sequences. As an example, cell-cycle timing depends on a sequence of ubiquitylation events mediated by the anaphase-promoting complex (APC) based on short redundant motifs. To understand the origin of specificity, we designed single-molecule fluorescence assays that capture transient ubiquitylation reactions. We find that the APC-mediated ubiquitylation involves a highly processive initial reaction on the substrate, followed by multiple encounters and reactions at a slower rate. The initial ubiquitylation greatly enhances the substrate's binding affinity in subsequent reactions, by both increasing the on-rate and decreasing the off-rate. We postulate that these cycles of positive feedback enable high specificity for substrates with short recognition motifs in a complex cellular environment.

Outlier detection at the transcriptome-proteome interface

BACKGROUND

In high-throughput experimental biology, it is widely acknowledged that while expression levels measured at the levels of transcriptome and the corresponding proteome do not, in general, correlate well, messenger RNA levels are used as convenient proxies for protein levels. Our interest is in developing data-driven computational models that can bridge the gap between these two levels of measurement at which different mechanisms of regulation may act on different molecular species causing any observed lack of correlations. To this end, we build data-driven predictors of protein levels using mRNA levels and known proxies of translation efficiencies as covariates. Previous work showed that in such a setting, outliers with respect to the model are reliable candidates for post-translational regulation.

RESULTS

Here, we present and compare two novel formulations of deriving a protein concentration predictor from which outliers may be extracted in a systematic manner. The first approach, outlier rejecting regression, allows explicit specification of a certain fraction of the data as outliers. In a regression setting, this is a non-convex optimization problem which we solve by deriving a difference of convex functions algorithm (DCA). With post-translationally regulated proteins, one expects their concentrations to be affected primarily by disruption of protein stability. Our second algorithm exploits this observation by minimizing an asymmetric loss using quantile regression and extracts outlier proteins whose measured concentrations are lower than what a genome-wide regression would predict. We validate the two approaches on a dataset of yeast transcriptome and proteome. Functional annotation check on detected outliers demonstrate that the methods are able to identify post-translationally regulated genes with high statistical confidence.

MiCroKiTS 4.0: a database of midbody, centrosome, kinetochore, telomere and spindle

We reported an updated database of MiCroKiTS 4.0 (http://microkit.biocuckoo.org) for proteins temporally and spatially localized in distinct subcellular positions including midbody, centrosome, kinetochore, telomere and mitotic spindle during cell division/mitosis. The database was updated from our previously developed database of MiCroKit 3.0, which contained 1489 proteins mostly forming super-complexes at midbody, centrosome and kinetochore from seven eukaryotes. Since the telomere and spindle apparatus are critical for cell division, the proteins localized at the two positions were also integrated. From the scientific literature, we curated 1872 experimentally identified proteins which at least locate in one of the five positions from eight species. Then the ortholog detection was performed to identify potential MiCroKiTS proteins from 144 eukaryotic organisms, which contains 66, 45 and 33 species of animals, fungi and plants, respectively. In total, 87 983 unique proteins with corresponding localization information were integrated into the database. The primary references of experimentally identified localizations were provided and the fluorescence microscope figures for the localizations of human proteins were shown. The orthologous relations between predicted and experimental localizations were also present. Taken together, we anticipate the database can serve as a useful resource for further analyzing the molecular mechanisms during cell division.

APC/C(Cdh1) controls CtIP stability during the cell cycle and in response to DNA damage

Human cells have evolved elaborate mechanisms for responding to DNA damage to maintain genome stability and prevent carcinogenesis. For instance, the cell cycle can be arrested at different stages to allow time for DNA repair. The APC/C(C) (dh1) ubiquitin ligase mainly regulates mitotic exit but is also implicated in the DNA damage-induced G2 arrest. However, it is currently unknown whether APC/C(C) (dh1) also contributes to DNA repair. Here, we show that Cdh1 depletion causes increased levels of genomic instability and enhanced sensitivity to DNA-damaging agents. Using an integrated proteomics and bioinformatics approach, we identify CtIP, a DNA-end resection factor, as a novel APC/C(C) (dh1) target. CtIP interacts with Cdh1 through a conserved KEN box, mutation of which impedes ubiquitylation and downregulation of CtIP both during G1 and after DNA damage in G2. Finally, we find that abrogating the CtIP-Cdh1 interaction results in delayed CtIP clearance from DNA damage foci, increased DNA-end resection, and reduced homologous recombination efficiency. Combined, our results highlight the impact of APC/C(C) (dh1) on the maintenance of genome integrity and show that this is, at least partially, achieved by controlling CtIP stability in a cell cycle- and DNA damage-dependent manner.

Clock gene expression in different synovial cells of patients with rheumatoid arthritis and osteoarthritis

Patients with rheumatoid arthritis (RA) show modulated circadian rhythms of inflammatory cytokines and cortisol, which may be associated with a modified expression of clock genes. The expression of major clock genes was previously studied in synovial tissues and fibroblasts of patients with RA and osteoarthritis (OA). We therefore especially aimed to examine the localization of clock genes at the cellular level in synovial tissue. Furthermore we were interested in studying the expression of the D site of albumin promoter (albumin D-box) binding protein (DBP) at the immunohistochemical level in human samples. Methods used include the in situ expression of the clock genes Brain and muscle aryl hydrocarbon receptor nuclear translocator-like 1 (Bmal 1), Circadian Locomotor Output Cycles Kaput (Clock), Period 1 and 2 (Per 1 and Per 2), and DBP was examined by immunohistochemistry in synovial tissues of patients with RA or OA. Additionally, expression profiles of different clock genes were determined over 24h by real time PCR in synovial fibroblasts (SFs) after a 2h serum shock or TNF-α. Results show that all clock genes investigated were found to be expressed both in RA and OA synovial tissues. Double staining against cell specific markers revealed that clock proteins were especially seen in macrophages, SFs and B-lymphocytes. Cell counting showed that clock proteins were found in approximately 5-20% of cells. Additionally, preliminary cell culture experiments showed that TNF-α treatment resulted in differential 24h expression profiles between RA and OA samples and also compared to the results obtained from the serum shock experiments. From our study we conclude that the major clock genes, including DBP, are expressed in samples from RA and OA patients, especially in macrophages and synovial fibroblasts, but also in B-lymphocytes. Preliminary experiments suggest that TNF-α seems to be able to modify clock gene expression in synovial fibroblasts.

Co-regulation proteomics reveals substrates and mechanisms of APC/C-dependent degradation

Using multiplexed quantitative proteomics, we analyzed cell cycle-dependent changes of the human proteome. We identified >4,400 proteins, each with a six-point abundance profile across the cell cycle. Hypothesizing that proteins with similar abundance profiles are co-regulated, we clustered the proteins with abundance profiles most similar to known Anaphase-Promoting Complex/Cyclosome (APC/C) substrates to identify additional putative APC/C substrates. This protein profile similarity screening (PPSS) analysis resulted in a shortlist enriched in kinases and kinesins. Biochemical studies on the kinesins confirmed KIFC1, KIF18A, KIF2C, and KIF4A as APC/C substrates. Furthermore, we showed that the APC/C(CDH1)-dependent degradation of KIFC1 regulates the bipolar spindle formation and proper cell division. A targeted quantitative proteomics experiment showed that KIFC1 degradation is modulated by a stabilizing CDK1-dependent phosphorylation site within the degradation motif of KIFC1. The regulation of KIFC1 (de-)phosphorylation and degradation provides insights into the fidelity and proper ordering of substrate degradation by the APC/C during mitosis.

The Aurora B kinase in Trypanosoma brucei undergoes post-translational modifications and is targeted to various subcellular locations through binding to TbCPC1

The chromosomal passenger complex (CPC) in animals, consisting of Aurora B kinase and three evolutionarily conserved proteins, plays crucial roles in mitosis and cytokinesis. However, Trypanosoma brucei expresses an unusual CPC consisting of an Aurora-like kinase, TbAUK1, and two kinetoplastid-specific proteins, TbCPC1 and TbCPC2. Despite their essential functions, little is known about the regulation of TbAUK1 and the roles of TbCPC1 and TbCPC2. Here, we investigate the effect of post-translational modification on the activity and spatiotemporal control of TbAUK1, and demonstrate that phosphorylation of two conserved threonine residues in the activation loop of the kinase domain contributes to TbAUK1 activation and function. TbAUK1 is SUMOylated in vivo, and mutation of the SUMO-conjugation site compromises TbAUK1 function. Degradation of TbAUK1 requires two destruction boxes and is mediated by the anaphase-promoting complex/cyclosome (APC/C), whereas degradation of TbCPC1 and TbCPC2 is not dependent on the predicted destruction boxes and is APC/C-independent. Moreover, we determine the domains in CPC subunits that mediate the pairwise interactions, and show that disruption of the interaction impairs the localization of TbAUK1 and TbCPC2 but not TbCPC1. Our results demonstrate the requirement of post-translational modifications for TbAUK1 function and a crucial role of TbCPC1 in mediating TbAUK1 localization.

Bridging the gap between transcriptome and proteome measurements identifies post-translationally regulated genes

MOTIVATION

Despite much dynamical cellular behaviour being achieved by accurate regulation of protein concentrations, messenger RNA abundances, measured by microarray technology, and more recently by deep sequencing techniques, are widely used as proxies for protein measurements. Although for some species and under some conditions, there is good correlation between transcriptome and proteome level measurements, such correlation is by no means universal due to post-transcriptional and post-translational regulation, both of which are highly prevalent in cells. Here, we seek to develop a data-driven machine learning approach to bridging the gap between these two levels of high-throughput omic measurements on Saccharomyces cerevisiae and deploy the model in a novel way to uncover mRNA-protein pairs that are candidates for post-translational regulation.

RESULTS

The application of feature selection by sparsity inducing regression (l₁ norm regularization) leads to a stable set of features: i.e. mRNA, ribosomal occupancy, ribosome density, tRNA adaptation index and codon bias while achieving a feature reduction from 37 to 5. A linear predictor used with these features is capable of predicting protein concentrations fairly accurately (R² = 0.86). Proteins whose concentration cannot be predicted accurately, taken as outliers with respect to the predictor, are shown to have annotation evidence of post-translational modification, significantly more than random subsets of similar size P < 0.02. In a data mining sense, this work also shows a wider point that outliers with respect to a learning method can carry meaningful information about a problem domain.

Global subcellular characterization of protein degradation using quantitative proteomics

Protein degradation provides an important regulatory mechanism used to control cell cycle progression and many other cellular pathways. To comprehensively analyze the spatial control of protein degradation in U2OS osteosarcoma cells, we have combined drug treatment and SILAC-based quantitative mass spectrometry with subcellular and protein fractionation. The resulting data set analyzed more than 74,000 peptides, corresponding to ∼5000 proteins, from nuclear, cytosolic, membrane, and cytoskeletal compartments. These data identified rapidly degraded proteasome targets, such as PRR11 and highlighted a feedback mechanism resulting in translation inhibition, induced by blocking the proteasome. We show this is mediated by activation of the unfolded protein response. We observed compartment-specific differences in protein degradation, including proteins that would not have been characterized as rapidly degraded through analysis of whole cell lysates. Bioinformatic analysis of the entire data set is presented in the Encyclopedia of Proteome Dynamics, a web-based resource, with proteins annotated for stability and subcellular distribution.