Kinase inhibitor pulldown assay (KiP) for clinical proteomics.

Protein kinases are frequently dysregulated and/or mutated in cancer and represent essential targets for therapy. Accurate quantification is essential. For breast cancer treatment, the identification and quantification of the protein kinase ERBB2 is critical for therapeutic decisions. While immunohistochemistry (IHC) is the current clinical diagnostic approach, it is only semiquantitative. Mass spectrometry-based proteomics offers quantitative assays that, unlike IHC, can be used to accurately evaluate hundreds of kinases simultaneously. The enrichment of less abundant kinase targets for quantification, along with depletion of interfering proteins, improves sensitivity and thus promotes more effective downstream analyses. Multiple kinase inhibitors were therefore deployed as a capture matrix for kinase inhibitor pulldown (KiP) assays designed to profile the human protein kinome as broadly as possible. Optimized assays were initially evaluated in 16 patient derived xenograft models (PDX) where KiP identified multiple differentially expressed and biologically relevant kinases. From these analyses, an optimized single-shot parallel reaction monitoring (PRM) method was developed to improve quantitative fidelity. The PRM KiP approach was then reapplied to low quantities of proteins typical of yields from core needle biopsies of human cancers. The initial prototype targeting 100 kinases recapitulated intrinsic subtyping of PDX models obtained from comprehensive proteomic and transcriptomic profiling. Luminal and HER2 enriched OCT-frozen patient biopsies subsequently analyzed through KiP-PRM also clustered by subtype. Finally, stable isotope labeled peptide standards were developed to define a prototype clinical method. Data are available via ProteomeXchange with identifiers PXD044655 and PXD046169.

Modulation of the proteostasis network promotes tumor resistance to oncogenic KRAS inhibitors.

Despite substantial advances in targeting mutant KRAS, tumor resistance to KRAS inhibitors (KRASi) remains a major barrier to progress. Here, we report proteostasis reprogramming as a key convergence point of multiple KRASi-resistance mechanisms. Inactivation of oncogenic KRAS down-regulated both the heat shock response and the inositol-requiring enzyme 1α (IRE1α) branch of the unfolded protein response, causing severe proteostasis disturbances. However, IRE1α was selectively reactivated in an ER stress-independent manner in acquired KRASi-resistant tumors, restoring proteostasis. Oncogenic KRAS promoted IRE1α protein stability through extracellular signal-regulated kinase (ERK)-dependent phosphorylation of IRE1α, leading to IRE1α disassociation from 3-hydroxy-3-methylglutaryl reductase degradation (HRD1) E3-ligase. In KRASi-resistant tumors, both reactivated ERK and hyperactivated AKT restored IRE1α phosphorylation and stability. Suppression of IRE1α overcame resistance to KRASi. This study reveals a druggable mechanism that leads to proteostasis reprogramming and facilitates KRASi resistance.

Kinome-wide screening uncovers a role for Bromodomain Protein 3 in DNA double-stranded break repair.

Double-stranded breaks (DSBs) are toxic DNA damage and a serious threat to genomic integrity. Thus, all living organisms have evolved multiple mechanisms of DNA DSB repair, the two principal ones being classical-non homologous end joining (C-NHEJ), and homology dependent recombination (HDR). In mammals, C-NHEJ is the predominate DSB repair pathway, but how a cell chooses to repair a particular DSB by a certain pathway is still not mechanistically clear. To uncover novel regulators of DSB repair pathway choice, we performed a kinome-wide screen in a human cell line engineered to express a dominant-negative C-NHEJ factor. The intellectual basis for such a screen was our hypothesis that a C-NHEJ-crippled cell line might need to upregulate other DSB repair pathways, including HDR, in order to survive. This screen identified Bromodomain-containing Protein 3 (BRD3) as a protein whose expression was almost completely ablated specifically in a C-NHEJ-defective cell line. Subsequent experimentation demonstrated that BRD3 is a negative regulator of HDR as BRD3-null cell lines proved to be hyper-recombinogenic for gene conversion, sister chromatid exchanges and gene targeting. Mechanistically, BRD3 appears to be working at the level of Radiation Sensitive 51 (RAD51) recruitment. Overall, our results demonstrate that BRD3 is a novel regulator of human DSB repair pathway choice.

Exploiting DNA Replication Stress as a Therapeutic Strategy for Breast Cancer.

Proliferating cells rely on DNA replication to ensure accurate genome duplication. Cancer cells, including breast cancer cells, exhibit elevated replication stress (RS) due to the uncontrolled oncogenic activation, loss of key tumor suppressors, and defects in the DNA repair machinery. This intrinsic vulnerability provides a great opportunity for therapeutic exploitation. An increasing number of drug candidates targeting RS in breast cancer are demonstrating promising efficacy in preclinical and early clinical trials. However, unresolved challenges lie in balancing the toxicity of these drugs while maintaining clinical efficacy. Furthermore, biomarkers of RS are urgently required to guide patient selection. In this review, we introduce the concept of targeting RS, detail the current therapies that target RS, and highlight the integration of RS with immunotherapies for breast cancer treatment. Additionally, we discuss the potential biomarkers to optimizing the efficacy of these therapies. Together, the continuous advances in our knowledge of targeting RS would benefit more patients with breast cancer.

Alpha-single chains of collagen type VI inhibit the fibrogenic effects of triple helical collagen VI in hepatic stellate cells.

The interaction of extracellular matrix (ECM) components with hepatic stellate cells (HSCs) is thought to perpetuate fibrosis by stimulating signaling pathways that drive HSC activation, survival and proliferation. Consequently, disrupting the interaction between ECM and HSCs is considered a therapeutical avenue although respective targets and underlying mechanisms remain to be established. Here we have interrogated the interaction between type VI collagen (CVI) and HSCs based on the observation that CVI is 10-fold upregulated during fibrosis, closely associates with HSCs in vivo and promotes cell proliferation and cell survival in cancer cell lines. We exposed primary rat HSCs and a rat hepatic stellate cell line (CFSC) to soluble CVI and determined the rate of proliferation, apoptosis and fibrogenesis in the absence of any additional growth factors. We find that CVI in nanomolar concentrations prevents serum starvation-induced apoptosis. This potent anti-apoptotic effect is accompanied by induction of proliferation and acquisition of a pronounced pro-fibrogenic phenotype characterized by increased α-smooth muscle actin, TGF-ß, collagen type I and TIMP-1 expression and diminished proteolytic MMP-13 expression. The CVI-HSC interaction can be disrupted with the monomeric α2(VI) and α3(VI) chains and abrogates the activating CVI effects. Further, functional relevant α3(VI)-derived 30 amino acid peptides lead to near-complete inhibition of the CVI effect. In conclusion, CVI serves as a potent mitogen and activating factor for HSCs. The antagonistic effects of the CVI monomeric chains and peptides point to linear peptide sequences that prevent activation of CVI receptors which may allow a targeted antifibrotic therapy.

Fructose Causes Liver Damage, Polyploidy, and Dysplasia in the Setting of Short Telomeres and p53 Loss.

Studies in humans and model systems have established an important role of short telomeres in predisposing to liver fibrosis through pathways that are incompletely understood. Recent studies have shown that telomere dysfunction impairs cellular metabolism, but whether and how these metabolic alterations contribute to liver fibrosis is not well understood. Here, we investigated whether short telomeres change the hepatic response to metabolic stress induced by fructose, a sugar that is highly implicated in non-alcoholic fatty liver disease. We find that telomere shortening in telomerase knockout mice (TKO) imparts a pronounced susceptibility to fructose as reflected in the activation of p53, increased apoptosis, and senescence, despite lower hepatic fat accumulation in TKO mice compared to wild type mice with long telomeres. The decreased fat accumulation in TKO is mediated by p53 and deletion of p53 normalizes hepatic fat content but also causes polyploidy, polynuclearization, dysplasia, cell death, and liver damage. Together, these studies suggest that liver tissue with short telomers are highly susceptible to fructose and respond with p53 activation and liver damage that is further exacerbated when p53 is lost resulting in dysplastic changes.

RON signalling promotes therapeutic resistance in ESR1 mutant breast cancer.

BACKGROUND: Oestrogen Receptor 1 (ESR1) mutations are frequently acquired in oestrogen receptor (ER)-positive metastatic breast cancer (MBC) patients who were treated with aromatase inhibitors (AI) in the metastatic setting. Acquired ESR1 mutations are associated with poor prognosis and there is a lack of effective therapies that selectively target these cancers. METHODS: We performed a proteomic kinome analysis in ESR1 Y537S mutant cells to identify hyperactivated kinases in ESR1 mutant cells. We validated Recepteur d'Origine Nantais (RON) and PI3K hyperactivity through phospho-immunoblot analysis, organoid growth assays, and in an in vivo patient-derived xenograft (PDX) metastatic model. RESULTS: We demonstrated that RON was hyperactivated in ESR1 mutant models, and in acquired palbociclib-resistant (PalbR) models. RON and insulin-like growth factor 1 receptor (IGF-1R) interacted as shown through pharmacological and genetic inhibition and were regulated by the mutant ER as demonstrated by reduced phospho-protein expression with endocrine therapies (ET). We show that ET in combination with a RON inhibitor (RONi) decreased ex vivo organoid growth of ESR1 mutant models, and as a monotherapy in PalbR models, demonstrating its therapeutic efficacy. Significantly, ET in combination with the RONi reduced metastasis of an ESR1 Y537S mutant PDX model. CONCLUSIONS: Our results demonstrate that RON/PI3K pathway inhibition may be an effective treatment strategy in ESR1 mutant and PalbR MBC patients. Clinically our data predict that ET resistance mechanisms can also contribute to CDK4/6 inhibitor resistance.

Neurofibromin Is an Estrogen Receptor-α Transcriptional Co-repressor in Breast Cancer.

We report that neurofibromin, a tumor suppressor and Ras-GAP (GTPase-activating protein), is also an estrogen receptor-α (ER) transcriptional co-repressor through leucine/isoleucine-rich motifs that are functionally independent of GAP activity. GAP activity, in turn, does not affect ER binding. Consequently, neurofibromin depletion causes estradiol hypersensitivity and tamoxifen agonism, explaining the poor prognosis associated with neurofibromin loss in endocrine therapy-treated ER+ breast cancer. Neurofibromin-deficient ER+ breast cancer cells initially retain sensitivity to selective ER degraders (SERDs). However, Ras activation does play a role in acquired SERD resistance, which can be reversed upon MEK inhibitor addition, and SERD/MEK inhibitor combinations induce tumor regression. Thus, neurofibromin is a dual repressor for both Ras and ER signaling, and co-targeting may treat neurofibromin-deficient ER+ breast tumors.

Microscaled proteogenomic methods for precision oncology.

Cancer proteogenomics promises new insights into cancer biology and treatment efficacy by integrating genomics, transcriptomics and protein profiling including modifications by mass spectrometry (MS). A critical limitation is sample input requirements that exceed many sources of clinically important material. Here we report a proteogenomics approach for core biopsies using tissue-sparing specimen processing and microscaled proteomics. As a demonstration, we analyze core needle biopsies from ERBB2 positive breast cancers before and 48-72 h after initiating neoadjuvant trastuzumab-based chemotherapy. We show greater suppression of ERBB2 protein and both ERBB2 and mTOR target phosphosite levels in cases associated with pathological complete response, and identify potential causes of treatment resistance including the absence of ERBB2 amplification, insufficient ERBB2 activity for therapeutic sensitivity despite ERBB2 amplification, and candidate resistance mechanisms including androgen receptor signaling, mucin overexpression and an inactive immune microenvironment. The clinical utility and discovery potential of proteogenomics at biopsy-scale warrants further investigation.

DPYSL3 modulates mitosis, migration, and epithelial-to-mesenchymal transition in claudin-low breast cancer.

A Clinical Proteomic Tumor Analysis Consortium (CPTAC) proteogenomic analysis prioritized dihydropyrimidinase-like-3 (DPYSL3) as a multilevel (RNA/protein/phosphoprotein) expression outlier specific to the claudin-low (CLOW) subset of triple-negative breast cancers. A PubMed informatics tool indicated a paucity of data in the context of breast cancer, which further prioritized DPYSL3 for study. DPYSL3 knockdown in DPYSL3-positive ([Formula: see text]) CLOW cell lines demonstrated reduced proliferation, yet enhanced motility and increased expression of epithelial-to-mesenchymal transition (EMT) markers, suggesting that DPYSL3 is a multifunctional signaling modulator. Slower proliferation in DPYSL3-negative ([Formula: see text]) CLOW cells was associated with accumulation of multinucleated cells, indicating a mitotic defect that was associated with a collapse of the vimentin microfilament network and increased vimentin phosphorylation. DPYSL3 also suppressed the expression of EMT regulators SNAIL and TWIST and opposed p21 activated kinase 2 (PAK2)-dependent migration. However, these EMT regulators in turn induce DPYSL3 expression, suggesting that DPYSL3 participates in negative feedback on EMT. In conclusion, DPYSL3 expression identifies CLOW tumors that will be sensitive to approaches that promote vimentin phosphorylation during mitosis and inhibitors of PAK signaling during migration and EMT.

gpGrouper: A Peptide Grouping Algorithm for Gene-Centric Inference and Quantitation of Bottom-Up Proteomics Data.

In quantitative mass spectrometry, the method by which peptides are grouped into proteins can have dramatic effects on downstream analyses. Here we describe gpGrouper, an inference and quantitation algorithm that offers an alternative method for assignment of protein groups by gene locus and improves pseudo-absolute iBAQ quantitation by weighted distribution of shared peptide areas. We experimentally show that distributing shared peptide quantities based on unique peptide peak ratios improves quantitation accuracy compared with conventional winner-take-all scenarios. Furthermore, gpGrouper seamlessly handles two-species samples such as patient-derived xenografts (PDXs) without ignoring the host species or species-shared peptides. This is a critical capability for proper evaluation of proteomics data from PDX samples, where stromal infiltration varies across individual tumors. Finally, gpGrouper calculates peptide peak area (MS1) based expression estimates from multiplexed isobaric data, producing iBAQ results that are directly comparable across label-free, isotopic, and isobaric proteomics approaches.

Functional Annotation of ESR1 Gene Fusions in Estrogen Receptor-Positive Breast Cancer.

RNA sequencing (RNA-seq) detects estrogen receptor alpha gene (ESR1) fusion transcripts in estrogen receptor-positive (ER+) breast cancer, but their role in disease pathogenesis remains unclear. We examined multiple ESR1 fusions and found that two, both identified in advanced endocrine treatment-resistant disease, encoded stable and functional fusion proteins. In both examples, ESR1-e6>YAP1 and ESR1-e6>PCDH11X, ESR1 exons 1-6 were fused in frame to C-terminal sequences from the partner gene. Functional properties include estrogen-independent growth, constitutive expression of ER target genes, and anti-estrogen resistance. Both fusions activate a metastasis-associated transcriptional program, induce cellular motility, and promote the development of lung metastasis. ESR1-e6>YAP1- and ESR1-e6>PCDH11X-induced growth remained sensitive to a CDK4/6 inhibitor, and a patient-derived xenograft (PDX) naturally expressing the ESR1-e6>YAP1 fusion was also responsive. Transcriptionally active ESR1 fusions therefore trigger both endocrine therapy resistance and metastatic progression, explaining the association with fatal disease progression, although CDK4/6 inhibitor treatment is predicted to be effective.

Combinatorial inhibition of PTPN12-regulated receptors leads to a broadly effective therapeutic strategy in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer diagnosed in more than 200,000 women each year and is recalcitrant to targeted therapies. Although TNBCs harbor multiple hyperactive receptor tyrosine kinases (RTKs), RTK inhibitors have been largely ineffective in TNBC patients thus far. We developed a broadly effective therapeutic strategy for TNBC that is based on combined inhibition of receptors that share the negative regulator PTPN12. Previously, we and others identified the tyrosine phosphatase PTPN12 as a tumor suppressor that is frequently inactivated in TNBC. PTPN12 restrains several RTKs, suggesting that PTPN12 deficiency leads to aberrant activation of multiple RTKs and a co-dependency on these receptors. This in turn leads to the therapeutic hypothesis that PTPN12-deficient TNBCs may be responsive to combined RTK inhibition. However, the repertoire of RTKs that are restrained by PTPN12 in human cells has not been systematically explored. By methodically identifying the suite of RTK substrates (MET, PDGFRß, EGFR, and others) inhibited by PTPN12, we rationalized a combination RTK-inhibitor therapy that induced potent tumor regression across heterogeneous models of TNBC. Orthogonal approaches revealed that PTPN12 was recruited to and inhibited these receptors after ligand stimulation, thereby serving as a feedback mechanism to limit receptor signaling. Cancer-associated mutation of PTPN12 or reduced PTPN12 protein levels diminished this feedback mechanism, leading to aberrant activity of these receptors. Restoring PTPN12 protein levels restrained signaling from RTKs, including PDGFRß and MET, and impaired TNBC survival. In contrast with single agents, combined inhibitors targeting the PDGFRß and MET receptors induced the apoptosis in TNBC cells in vitro and in vivo. This therapeutic strategy resulted in tumor regressions in chemo-refractory patient-derived TNBC models. Notably, response correlated with PTPN12 deficiency, suggesting that impaired receptor feedback may establish a combined addiction to these proto-oncogenic receptors. Taken together, our data provide a rationale for combining RTK inhibitors in TNBC and other malignancies that lack receptor-activating mutations.

The Histone Variant MacroH2A1 Is a BRCA1 Ubiquitin Ligase Substrate.

The breast- and ovarian-cancer-specific tumor suppressor BRCA1 and its heterodimeric partner BARD1 contain RING domains that implicate them as E3 ubiquitin ligases. Despite extensive efforts, the bona fide substrates of BRCA1/BARD1 remain elusive. Here, we used recombinant GST fused to four UBA domains to enrich ubiquitinated proteins followed by a Lys-ε-Gly-Gly (diGly) antibody to enrich ubiquitinated tryptic peptides. This tandem affinity purification method coupled with mass spectrometry identified 101 putative BRCA1/BARD1 E3 substrates. We identified the histone variant macroH2A1 from the screen and showed that BRCA1/BARD1 ubiquitinates macroH2A1 at lysine 123 in vitro and in vivo. Primary human fibroblasts stably expressing a ubiquitination-deficient macroH2A1 mutant were defective in cellular senescence compared to their wild-type counterpart. Our study demonstrates that BRCA1/BARD1 is a macroH2A1 E3 ligase and implicates a role for macroH2A1 K123 ubiquitination in cellular senescence.

The Dual Estrogen Receptor α Inhibitory Effects of the Tissue-Selective Estrogen Complex for Endometrial and Breast Safety.

The conjugated estrogen /: bazedoxifene tissue-selective estrogen complex (TSEC) is designed to minimize the undesirable effects of estrogen in the uterus and breast tissues and to allow the beneficial effects of estrogen in other estrogen-target tissues, such as the bone and brain. However, the molecular mechanism underlying endometrial and breast safety during TSEC use is not fully understood. Estrogen receptor α (ERα)-estrogen response element (ERE)-DNA pull-down assays using HeLa nuclear extracts followed by mass spectrometry-immunoblotting analyses revealed that, upon TSEC treatment, ERα interacted with transcriptional repressors rather than coactivators. Therefore, the TSEC-mediated recruitment of transcriptional repressors suppresses ERα-mediated transcription in the breast and uterus. In addition, TSEC treatment also degraded ERα protein in uterine tissue and breast cancer cells, but not in bone cells. Interestingly, ERα-ERE-DNA pull-down assays also revealed that, upon TSEC treatment, ERα interacted with the F-box protein 45 (FBXO45) E3 ubiquitin ligase. The loss-of- and gain-of-FBXO45 function analyses indicated that FBXO45 is involved in TSEC-mediated degradation of the ERα protein in endometrial and breast cells. In preclinical studies, these synergistic effects of TSEC on ERα inhibition also suppressed the estrogen-dependent progression of endometriosis. Therefore, the endometrial and breast safety effects of TSEC are associated with synergy between the selective recruitment of transcriptional repressors to ERα and FBXO45-mediated degradation of the ERα protein.

Proteomic analysis of coregulators bound to ERα on DNA and nucleosomes reveals coregulator dynamics.

Chromatin immunoprecipitation studies have mapped protein occupancies at many genomic loci. However, a detailed picture of the complexity of coregulators (CoRs) bound to a defined enhancer along with a transcription factor is missing. To address this, we used biotin-DNA pull-down assays coupled with mass spectrometry-immunoblotting to identify at least 17 CoRs from nuclear extracts bound to 17ß-estradiol (E2)-liganded estrogen receptor-α on estrogen response elements (EREs). Unexpectedly, these complexes initially are biochemically stable and contain certain atypical corepressors. Addition of ATP dynamically converts these complexes to an "activated" state by phosphorylation events, primarily mediated by DNA-dependent protein kinase. Importantly, a "natural" ERE-containing enhancer and nucleosomal EREs recruit similar complexes. We further discovered the mechanism whereby H3K4me3 stimulates ERα-mediated transcription as compared with unmodified nucleosomes. H3K4me3 templates promote specific CoR dynamics in the presence of ATP and AcCoA, as manifested by CBP/p300 and SRC-3 dismissal and SAGA and TFIID stabilization/recruitment.

Proteome-wide profiling of activated transcription factors with a concatenated tandem array of transcription factor response elements.

Transcription factors (TFs) are families of proteins that bind to specific DNA sequences, or TF response elements (TFREs), and function as regulators of many cellular processes. Because of the low abundance of TFs, direct quantitative measurement of TFs on a proteome scale remains a challenge. In this study, we report the development of an affinity reagent that permits identification of endogenous TFs at the proteome scale. The affinity reagent is composed of a synthetic DNA containing a concatenated tandem array of the consensus TFREs (catTFRE) for the majority of TF families. By using catTFRE to enrich TFs from cells, we were able to identify as many as 400 TFs from a single cell line and a total of 878 TFs from 11 cell types, covering more than 50% of the gene products that code for the DNA-binding TFs in the genome. We further demonstrated that catTFRE pull-downs could quantitatively measure proteome-wide changes in DNA binding activity of TFs in response to exogenous stimulation by using a label-free MS-based quantification approach. Applying catTFRE on the evaluation of drug effects, we described a panoramic view of TF activations and provided candidates for the elucidation of molecular mechanisms of drug actions. We anticipate that the catTFRE affinity strategy will find widespread applications in biomedical research.

Identification and characterization of nardilysin as a novel dimethyl H3K4-binding protein involved in transcriptional regulation.

Histone methylation on lysine residues is believed to function primarily as docking sites to recruit specific proteins termed as histone code "readers" or "effectors." Each lysine residue can be mono-, di, and tri-methylated and different methylation states can have different effect on chromatin function. While an increasing number of proteins have been identified and characterized as specific effectors for methylated histones, very few of the proteins are known to recognize a particular state of methylation. In this study, we identified nardilysin (NRDc), a member of M16 family metalloendopeptidases, as a novel dimethyl-H3K4 (H3K4me2)-binding protein. Among three methylated states, NRDc binds preferentially H3K4me2 both in vitro and in vivo. Biochemical purification demonstrated that NRDc interacts with the NCoR/SMRT corepressor complex. We identified target genes repressed by NRDc through microarray. We showed that NRDc is physically associated with and recruits the NCoR complex to some of the repressed genes and this association correlates with binding of H3K4me2. Thus, our study has identified a novel H3K4me2-binding protein and revealed a role of NRDc in transcriptional regulation.

Analysis of the human endogenous coregulator complexome.

Elucidation of endogenous cellular protein-protein interactions and their networks is most desirable for biological studies. Here we report our study of endogenous human coregulator protein complex networks obtained from integrative mass spectrometry-based analysis of 3290 affinity purifications. By preserving weak protein interactions during complex isolation and utilizing high levels of reciprocity in the large dataset, we identified many unreported protein associations, such as a transcriptional network formed by ZMYND8, ZNF687, and ZNF592. Furthermore, our work revealed a tiered interplay within networks that share common proteins, providing a conceptual organization of a cellular proteome composed of minimal endogenous modules (MEMOs), complex isoforms (uniCOREs), and regulatory complex-complex interaction networks (CCIs). This resource will effectively fill a void in linking correlative genomic studies with an understanding of transcriptional regulatory protein functions within the proteome for formulation and testing of future hypotheses.

A data set of human endogenous protein ubiquitination sites.

Lysine ubiquitination is an important and versatile protein post-translational modification. Numerous cellular functions are regulated by ubiquitination, suggesting that extensive numbers of proteins, if not all, are modified with ubiquitin at certain times. However, proteome-wide profiling of ubiquitination sites in the mammalian system is technically challenging. We report the design and characterization of an engineered protein affinity reagent for the isolation of ubiquitinated proteins and the identification of ubiquitination sites with mass spectrometry. This recombinant protein consists of four tandem repeats of ubiquitin-associated domain from UBQLN1 fused to a GST tag. We used this GST-qUBA reagent to isolate polyubiquitinated proteins and identified 294 endogenous ubiquitination sites on 223 proteins from human 293T cells without proteasome inhibitors or overexpression of ubiquitin. Mitochondrial proteins constitute 14.7% of this data set, implicating ubiquitination in a wide range of mitochondrial functions.