Global, site-resolved analysis of ubiquitylation occupancy and turnover rate reveals systems properties.

Ubiquitylation regulates most proteins and biological processes in a eukaryotic cell. However, the site-specific occupancy (stoichiometry) and turnover rate of ubiquitylation have not been quantified. Here we present an integrated picture of the global ubiquitylation site occupancy and half-life. Ubiquitylation site occupancy spans over four orders of magnitude, but the median ubiquitylation site occupancy is three orders of magnitude lower than that of phosphorylation. The occupancy, turnover rate, and regulation of sites by proteasome inhibitors are strongly interrelated, and these attributes distinguish sites involved in proteasomal degradation and cellular signaling. Sites in structured protein regions exhibit longer half-lives and stronger upregulation by proteasome inhibitors than sites in unstructured regions. Importantly, we discovered a surveillance mechanism that rapidly and site-indiscriminately deubiquitylates all ubiquitin-specific E1 and E2 enzymes, protecting them against accumulation of bystander ubiquitylation. The work provides a systems-scale, quantitative view of ubiquitylation properties and reveals general principles of ubiquitylation-dependent governance.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy.

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Time-Resolved Analysis Reveals Rapid Dynamics and Broad Scope of the CBP/p300 Acetylome.

The acetyltransferases CBP and p300 are multifunctional transcriptional co-activators. Here, we combined quantitative proteomics with CBP/p300-specific catalytic inhibitors, bromodomain inhibitor, and gene knockout to reveal a comprehensive map of regulated acetylation sites and their dynamic turnover rates. CBP/p300 acetylates thousands of sites, including signature histone sites and a multitude of sites on signaling effectors and enhancer-associated transcriptional regulators. Time-resolved acetylome analyses identified a subset of CBP/p300-regulated sites with very rapid (<30 min) acetylation turnover, revealing a dynamic balance between acetylation and deacetylation. Quantification of acetylation, mRNA, and protein abundance after CBP/p300 inhibition reveals a kinetically competent network of gene expression that strictly depends on CBP/p300-catalyzed rapid acetylation. Collectively, our in-depth acetylome analyses reveal systems attributes of CBP/p300 targets, and the resource dataset provides a framework for investigating CBP/p300 functions and for understanding the impact of small-molecule inhibitors targeting its catalytic and bromodomain activities.

Proteome dynamics at broken replication forks reveal a distinct ATM-directed repair response suppressing DNA double-strand break ubiquitination.

Cells have evolved an elaborate DNA repair network to ensure complete and accurate DNA replication. Defects in these repair machineries can fuel genome instability and drive carcinogenesis while creating vulnerabilities that may be exploited in therapy. Here, we use nascent chromatin capture (NCC) proteomics to characterize the repair of replication-associated DNA double-strand breaks (DSBs) triggered by topoisomerase 1 (TOP1) inhibitors. We reveal profound changes in the fork proteome, including the chromatin environment and nuclear membrane interactions, and identify three classes of repair factors according to their enrichment at broken and/or stalled forks. ATM inhibition dramatically rewired the broken fork proteome, revealing that ataxia telangiectasia mutated (ATM) signalling stimulates DNA end resection, recruits PLK1, and concomitantly suppresses the canonical DSB ubiquitination response by preventing accumulation of RNF168 and BRCA1-A. This work and collection of replication fork proteomes provide a new framework to understand how cells orchestrate homologous recombination repair of replication-associated DSBs.

Enhancers are activated by p300/CBP activity-dependent PIC assembly, RNAPII recruitment, and pause release.

The metazoan-specific acetyltransferase p300/CBP is involved in activating signal-induced, enhancer-mediated transcription of cell-type-specific genes. However, the global kinetics and mechanisms of p300/CBP activity-dependent transcription activation remain poorly understood. We performed genome-wide, time-resolved analyses to show that enhancers and super-enhancers are dynamically activated through p300/CBP-catalyzed acetylation, deactivated by the opposing deacetylase activity, and kinetic acetylation directly contributes to maintaining cell identity at very rapid (minutes) timescales. The acetyltransferase activity is dispensable for the recruitment of p300/CBP and transcription factors but essential for promoting the recruitment of TFIID and RNAPII at virtually all enhancers and enhancer-regulated genes. This identifies pre-initiation complex assembly as a dynamically controlled step in the transcription cycle and reveals p300/CBP-catalyzed acetylation as the signal that specifically promotes transcription initiation at enhancer-regulated genes. We propose that p300/CBP activity uses a "recruit-and-release" mechanism to simultaneously promote RNAPII recruitment and pause release and thereby enables kinetic activation of enhancer-mediated transcription.

SIK2 orchestrates actin-dependent host response upon Salmonella infection.

Salmonella is an intracellular pathogen of a substantial global health concern. In order to identify key players involved in Salmonella infection, we performed a global host phosphoproteome analysis subsequent to bacterial infection. Thereby, we identified the kinase SIK2 as a central component of the host defense machinery upon Salmonella infection. SIK2 depletion favors the escape of bacteria from the Salmonella-containing vacuole (SCV) and impairs Xenophagy, resulting in a hyperproliferative phenotype. Mechanistically, SIK2 associates with actin filaments under basal conditions; however, during bacterial infection, SIK2 is recruited to the SCV together with the elements of the actin polymerization machinery (Arp2/3 complex and Formins). Notably, SIK2 depletion results in a severe pathological cellular actin nucleation and polymerization defect upon Salmonella infection. We propose that SIK2 controls the formation of a protective SCV actin shield shortly after invasion and orchestrates the actin cytoskeleton architecture in its entirety to control an acute Salmonella infection after bacterial invasion.

Automating UbiFast for High-throughput and Multiplexed Ubiquitin Enrichment.

Robust methods for deep-scale enrichment and site-specific identification of ubiquitylation sites are necessary for characterizing the myriad roles of protein ubiquitylation. To this end we previously developed UbiFast, a sensitive method for highly multiplexed ubiquitylation profiling where K-ϵ-GG peptides are enriched with anti-K-ε-GG antibody and labeled on-antibody with isobaric labeling reagents for sample multiplexing. Here, we present robotic automation of the UbiFast method using a magnetic bead-conjugated K-ε-GG antibody (mK-ε-GG) and a magnetic particle processor. We report the identification of ∼20,000 ubiquitylation sites from a TMT10-plex with 500 µg input per sample processed in ∼2 h. Automation of the UbiFast method greatly increased the number of identified and quantified ubiquitylation sites, improved reproducibility, and significantly reduced processing time. The automated method also significantly reduced variability across process replicates compared with the manual method. The workflow enables processing of up to 96 samples in a single day making it suitable to study ubiquitylation in large sample sets. Here we demonstrate the applicability of the method to profile small amounts of tissue using breast cancer patient-derived xenograft (PDX) tissue samples.

A highly multiplexed quantitative phosphosite assay for biology and preclinical studies.

Reliable methods to quantify dynamic signaling changes across diverse pathways are needed to better understand the effects of disease and drug treatment in cells and tissues but are presently lacking. Here, we present SigPath, a targeted mass spectrometry (MS) assay that measures 284 phosphosites in 200 phosphoproteins of biological interest. SigPath probes a broad swath of signaling biology with high throughput and quantitative precision. We applied the assay to investigate changes in phospho-signaling in drug-treated cancer cell lines, breast cancer preclinical models, and human medulloblastoma tumors. In addition to validating previous findings, SigPath detected and quantified a large number of differentially regulated phosphosites newly associated with disease models and human tumors at baseline or with drug perturbation. Our results highlight the potential of SigPath to monitor phosphoproteomic signaling events and to nominate mechanistic hypotheses regarding oncogenesis, response, and resistance to therapy.

TMT Labeling for the Masses: A Robust and Cost-efficient, In-solution Labeling Approach.

Isobaric stable isotope labeling using, for example, tandem mass tags (TMTs) is increasingly being applied for large-scale proteomic studies. Experiments focusing on proteoform analysis in drug time course or perturbation studies or in large patient cohorts greatly benefit from the reproducible quantification of single peptides across samples. However, such studies often require labeling of hundreds of micrograms of peptides such that the cost for labeling reagents represents a major contribution to the overall cost of an experiment. Here, we describe and evaluate a robust and cost-effective protocol for TMT labeling that reduces the quantity of required labeling reagent by a factor of eight and achieves complete labeling. Under- and overlabeling of peptides derived from complex digests of tissues and cell lines were systematically evaluated using peptide quantities of between 12.5 and 800 µg and TMT-to-peptide ratios (wt/wt) ranging from 8:1 to 1:2 at different TMT and peptide concentrations. When reaction volumes were reduced to maintain TMT and peptide concentrations of at least 10 mm and 2 g/l, respectively, TMT-to-peptide ratios as low as 1:1 (wt/wt) resulted in labeling efficiencies of > 99% and excellent intra- and interlaboratory reproducibility. The utility of the optimized protocol was further demonstrated in a deep-scale proteome and phosphoproteome analysis of patient-derived xenograft tumor tissue benchmarked against the labeling procedure recommended by the TMT vendor. Finally, we discuss the impact of labeling reaction parameters for N-hydroxysuccinimide ester-based chemistry and provide guidance on adopting efficient labeling protocols for different peptide quantities.

SPATA2 links CYLD to the TNF-α receptor signaling complex and modulates the receptor signaling outcomes.

TNF-α is a key regulator of innate immune and proinflammatory responses. However, the composition of the TNF-α receptor-associated signaling complexes (TNF-RSC) and the architecture of the downstream signaling networks are incompletely understood. We employed quantitative mass spectrometry to demonstrate that TNF-α stimulation induces widespread protein phosphorylation and that the scope of phosphorylation expands in a temporal manner. TNF-α stimulation also induces rapid ubiquitylation of components of the TNF-RSC Temporal analysis of the TNF-RSC composition identified SPATA2 as a novel component of the TNF-RSC The predicted PUB domain in the N-terminus of SPATA2 interacts with the USP domain of CYLD, whereas the C-terminus of SPATA2 interacts with HOIP SPATA2 is required for recruitment of CYLD to the TNF-RSC Downregulation of SPATA2 augments transcriptional activation of NF-κB and inhibits TNF-α-induced necroptosis, pointing to an important function of SPATA2 in modulating the outcomes of TNF-α signaling. Taken together, our study draws a detailed map of TNF-α signaling, identifies SPATA2 as a novel component of TNF-α signaling, and provides a rich resource for further functional investigations.

Evaluation of Advanced Precursor Determination for Tandem Mass Tag (TMT)-Based Quantitative Proteomics across Instrument Platforms.

Tandem mass tag (TMT)-based quantitation is a strong modality for quantitative proteomics, as samples can be multiplexed, creating large-scale data sets with high precision and minimal missing values. However, coisolation/cofragmentation of near isobaric, coeluting precursor peptide analytes has been well-documented to show ratio compression, compromising the accuracy of peptide/protein quantitation. Advanced peak determination (APD) is a new peak-picking algorithm that shows improved identification of peak detection in survey scans (MS1) to increase the number of precursors selected for unimolecular dissociation (MS2). To increase the number of these "features" selected for MS2 APD purposefully selects multiple peptide precursors of very similar m/ z that often derive from different proteins-a major source of ratio compression in TMT quantification. Here, we evaluate the effects of various data acquisition parameters combined with APD on ratio compression. We find that data acquisition with APD enabled results in more coisolated precursors, more mixed spectra, and in turn, fewer peptide spectral matches, especially at standard on-column loads. We conclude that APD should not be utilized for isobaric tagging, MS2-based experiments.

Accurate Quantification of Site-specific Acetylation Stoichiometry Reveals the Impact of Sirtuin Deacetylase CobB on the E. coli Acetylome.

Lysine acetylation is a protein posttranslational modification (PTM) that occurs on thousands of lysine residues in diverse organisms from bacteria to humans. Accurate measurement of acetylation stoichiometry on a proteome-wide scale remains challenging. Most methods employ a comparison of chemically acetylated peptides to native acetylated peptides, however, the potentially large differences in abundance between these peptides presents a challenge for accurate quantification. Stable isotope labeling by amino acids in cell culture (SILAC)-based mass spectrometry (MS) is one of the most widely used quantitative proteomic methods. Here we show that serial dilution of SILAC-labeled peptides (SD-SILAC) can be used to identify accurately quantified peptides and to estimate the quantification error rate. We applied SD-SILAC to determine absolute acetylation stoichiometry in exponentially-growing and stationary-phase wild-type and Sirtuin deacetylase CobB-deficient cells. To further analyze CobB-regulated sites under conditions of globally increased or decreased acetylation, we measured stoichiometry in phophotransacetylase (ptaΔ) and acetate kinase (ackAΔ) mutant strains in the presence and absence of the Sirtuin inhibitor nicotinamide. We measured acetylation stoichiometry at 3,669 unique sites and found that the vast majority of acetylation occurred at a low stoichiometry. Manipulations that cause increased nonenzymatic acetylation by acetyl-phosphate (AcP), such as stationary-phase arrest and deletion of ackA, resulted in globally increased acetylation stoichiometry. Comparison to relative quantification under the same conditions validated our stoichiometry estimates at hundreds of sites, demonstrating the accuracy of our method. Similar to Sirtuin deacetylase 3 (SIRT3) in mitochondria, CobB suppressed acetylation to lower than median stoichiometry in WT, ptaΔ, and ackAΔ cells. Together, our results provide a detailed view of acetylation stoichiometry in E. coli and suggest an evolutionarily conserved function of Sirtuin deacetylases in suppressing low stoichiometry acetylation.

Systems-wide analysis of BCR signalosomes and downstream phosphorylation and ubiquitylation.

B-cell receptor (BCR) signaling is essential for the development and function of B cells; however, the spectrum of proteins involved in BCR signaling is not fully known. Here we used quantitative mass spectrometry-based proteomics to monitor the dynamics of BCR signaling complexes (signalosomes) and to investigate the dynamics of downstream phosphorylation and ubiquitylation signaling. We identify most of the previously known components of BCR signaling, as well as many proteins that have not yet been implicated in this system. BCR activation leads to rapid tyrosine phosphorylation and ubiquitylation of the receptor-proximal signaling components, many of which are co-regulated by both the modifications. We illustrate the power of multilayered proteomic analyses for discovering novel BCR signaling components by demonstrating that BCR-induced phosphorylation of RAB7A at S72 prevents its association with effector proteins and with endo-lysosomal compartments. In addition, we show that BCL10 is modified by LUBAC-mediated linear ubiquitylation, and demonstrate an important function of LUBAC in BCR-induced NF-κB signaling. Our results offer a global and integrated view of BCR signaling, and the provided datasets can serve as a valuable resource for further understanding BCR signaling networks.

SUMOylation of the ING1b tumor suppressor regulates gene transcription.

The INhibitor of Growth (ING) proteins are encoded as multiple isoforms in five ING genes (ING1 -5) and act as type II tumor suppressors. They are growth inhibitory when overexpressed and are frequently mislocalized or downregulated in several forms of cancer. ING1 and ING2 are stoichiometric members of histone deacetylase complexes, whereas ING3-5 are stoichiometric components of different histone acetyltransferase complexes. The INGs target these complexes to histone marks, thus acting as epigenetic regulators. ING proteins affect angiogenesis, apoptosis, DNA repair, metastasis and senescence, but how the proteins themselves are regulated is not yet clear. Here, we find a small ubiquitin-like modification (SUMOylation) of the ING1b protein and identify lysine 193 (K193) as the preferred ING1b SUMO acceptor site. We also show that PIAS4 is the E3 SUMO ligase responsible for ING1b SUMOylation on K193. Sequence alignment reveals that the SUMO consensus site on ING1b contains a phosphorylation-dependent SUMOylation motif (PDSM) and our data indicate that the SUMOylation on K193 is enhanced by the S199D phosphomimic mutant. Using an ING1b protein mutated at the major SUMOylation site (ING1b E195A), we further demonstrate that ING1b SUMOylation regulates the binding of ING1b to the ISG15 and DGCR8 promoters, consequently regulating ISG15 and DGCR8 transcription. These results suggest a role for ING1b SUMOylation in the regulation of gene transcription.