Global proteomic assessment of the classical protein-tyrosine phosphatome and "Redoxome".

Protein-tyrosine phosphatases (PTPs), along with protein-tyrosine kinases, play key roles in cellular signaling. All Class I PTPs contain an essential active site cysteinyl residue, which executes a nucleophilic attack on substrate phosphotyrosyl residues. The high reactivity of the catalytic cysteine also predisposes PTPs to oxidation by reactive oxygen species, such as H(2)O(2). Reversible PTP oxidation is emerging as an important cellular regulatory mechanism and might contribute to diseases such as cancer. We exploited these unique features of PTP enzymology to develop proteomic methods, broadly applicable to cell and tissue samples, that enable the comprehensive identification and quantification of expressed classical PTPs (PTPome) and the oxidized subset of the PTPome (oxPTPome). We find that mouse and human cells and tissues, including cancer cells, display distinctive PTPomes and oxPTPomes, revealing additional levels of complexity in the regulation of protein-tyrosine phosphorylation in normal and malignant cells.

The ubiquitin ligase HERC4 suppresses MafA transcriptional activity triggered by GSK3ß in myeloma by atypical K63-linked polyubiquitination.

MafA and c-Maf are close members of the Maf transcription factor family and indicators of poor prognosis of multiple myeloma (MM). Our previous study finds that the ubiquitin ligase HERC4 induces c-Maf degradation but stabilizes MafA, and the mechanism is elusive. In the present study, we find that HERC4 interacts with MafA and mediates its K63-linked polyubiquitination at K33. Moreover, HERC4 inhibits MafA phosphorylation and its transcriptional activity triggered by glycogen synthase kinase 3ß (GSK3ß). The K33R MafA variant prevents HERC4 from inhibiting MafA phosphorylation and increases MafA transcriptional activity. Further analyses reveal that MafA can also activate the STAT3 signaling, but it is suppressed by HERC4. Lastly, we demonstrate that lithium chloride, a GSK3ß inhibitor, can upregulate HERC4 and synergizes dexamethasone, a typical anti-MM drug, in inhibiting MM cell proliferation and xenograft growth in nude mice. These findings thus highlight a novel regulation of MafA oncogenic activity in MM and provide the rationale by targeting HERC4/GSK3ß/MafA for the treatment of MM.

Targeting the Otub1/c-Maf axis for the treatment of multiple myeloma.

The oncogenic transcription factor c-Maf has been proposed as an ideal therapeutic target for multiple myeloma (MM), but how to achieve it is still elusive. In the present study, we found the Otub1/c-Maf axis could be a potential target. Otub1, an OTU family deubiquitinase, was found to interact with c-Maf by mass spectrometry. Otub1 abrogates c-Maf K48-linked polyubiquitination, thus preventing its degradation and enhancing its transcriptional activity. Specifically, this deubiquitinating activity depends on its Lys71 and the N terminus but is independent of UBE2O, a known E2 of c-Maf. Otub1 promotes MM cell survival and MM tumor growth. In contrast, silence of Otub1 leads to c-Maf degradation and c-Maf-expressing MM cell apoptosis. Therefore, the Otub1/c-Maf axis could be a therapeutic target of MM. In order to explore this concept, we performed a c-Maf recognition element-driven luciferase-based screen against US Food and Drug Administration-approved drugs and natural products, from which the generic cardiac glycoside lanatoside C (LanC) is found to prevent c-Maf deubiquitination and induces its degradation by disrupting the interaction of Otub1 and c-Maf. Consequently, LanC inhibits c-Maf transcriptional activity, induces c-Maf-expressing MM cell apoptosis, and suppresses MM growth and prolongs overall survival of model mice, but without apparent toxicity. Therefore, the present study identifies Otub1 as a novel deubiquitinase of c-Maf and establishes that the Otub1/c-Maf axis is a potential therapeutic target for MM.

Oncohistone interactome profiling uncovers contrasting oncogenic mechanisms and identifies potential therapeutic targets in high grade glioma.

Histone H3 mutations at amino acids 27 (H3K27M) and 34 (H3G34R) are recurrent drivers of pediatric-type high-grade glioma (pHGG). H3K27M mutations lead to global disruption of H3K27me3 through dominant negative PRC2 inhibition, while H3G34R mutations lead to local losses of H3K36me3 through inhibition of SETD2. However, their broader oncogenic mechanisms remain unclear. We characterized the H3.1K27M, H3.3K27M and H3.3G34R interactomes, finding that H3K27M is associated with epigenetic and transcription factor changes; in contrast H3G34R removes a break on cryptic transcription, limits DNA methyltransferase access, and alters mitochondrial metabolism. All 3 mutants had altered interactions with DNA repair proteins and H3K9 methyltransferases. H3K9me3 was reduced in H3K27M-containing nucleosomes, and cis-H3K9 methylation was required for H3K27M to exert its effect on global H3K27me3. H3K9 methyltransferase inhibition was lethal to H3.1K27M, H3.3K27M and H3.3G34R pHGG cells, underscoring the importance of H3K9 methylation for oncohistone-mutant gliomas and suggesting it as an attractive therapeutic target.

A drug discovery platform to identify compounds that inhibit EGFR triple mutants.

Receptor tyrosine kinases (RTKs) are transmembrane receptors of great clinical interest due to their role in disease. Historically, therapeutics targeting RTKs have been identified using in vitro kinase assays. Due to frequent development of drug resistance, however, there is a need to identify more diverse compounds that inhibit mutated but not wild-type RTKs. Here, we describe MaMTH-DS (mammalian membrane two-hybrid drug screening), a live-cell platform for high-throughput identification of small molecules targeting functional protein-protein interactions of RTKs. We applied MaMTH-DS to an oncogenic epidermal growth factor receptor (EGFR) mutant resistant to the latest generation of clinically approved tyrosine kinase inhibitors (TKIs). We identified four mutant-specific compounds, including two that would not have been detected by conventional in vitro kinase assays. One of these targets mutant EGFR via a new mechanism of action, distinct from classical TKI inhibition. Our results demonstrate how MaMTH-DS is a powerful complement to traditional drug screening approaches.

The deubiquitinase USP7 stabilizes Maf proteins to promote myeloma cell survival.

The Maf proteins, including c-Maf, MafA, and MafB, are critical transcription factors in myelomagenesis. Previous studies demonstrated that Maf proteins are processed by the ubiquitin-proteasome pathway, but the mechanisms remain elusive. This study applied MS to identify MafB ubiquitination-associated proteins and found that the ubiquitin-specific protease USP7 was present in the MafB interactome. Moreover, USP7 also interacted with c-Maf and MafA and blocked their polyubiquitination and degradation. Consistently, knockdown of USP7 resulted in Maf protein degradation along with increased polyubiquitination levels. The action of USP7 thus promoted Maf transcriptional activity as evidenced by luciferase assays and by the up-regulation of the expression of Maf-modulated genes. Furthermore, USP7 was up-regulated in myeloma cells, and it was negatively associated with the survival of myeloma patients. USP7 promoted myeloma cell survival, and when it was inhibited by its specific inhibitor P5091, myeloma cell lines underwent apoptosis. These results therefore demonstrated that USP7 is a deubiquitinase of Maf proteins and promotes MM cell survival in association with Maf stability. Given the significance of USP7 and Maf proteins in myeloma genesis, targeting the USP7/Maf axle is a potential strategy to the precision therapy of MM.

Refined RIP-seq protocol for epitranscriptome analysis with low input materials.

N6-Methyladenosine (m6A) accounts for approximately 0.2% to 0.6% of all adenosine in mammalian mRNA, representing the most abundant internal mRNA modifications. m6A RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq) is a powerful technique to map the m6A location transcriptome-wide. However, this method typically requires 300 µg of total RNA, which limits its application to patient tumors. In this study, we present a refined m6A MeRIP-seq protocol and analysis pipeline that can be applied to profile low-input RNA samples from patient tumors. We optimized the key parameters of m6A MeRIP-seq, including the starting amount of RNA, RNA fragmentation, antibody selection, MeRIP washing/elution conditions, methods for RNA library construction, and the bioinformatics analysis pipeline. With the optimized immunoprecipitation (IP) conditions and a postamplification rRNA depletion strategy, we were able to profile the m6A epitranscriptome using 500 ng of total RNA. We identified approximately 12,000 m6A peaks with a high signal-to-noise (S/N) ratio from 2 lung adenocarcinoma (ADC) patient tumors. Through integrative analysis of the transcriptome, m6A epitranscriptome, and proteome data in the same patient tumors, we identified dynamics at the m6A level that account for the discordance between mRNA and protein levels in these tumors. The refined m6A MeRIP-seq method is suitable for m6A epitranscriptome profiling in a limited amount of patient tumors, setting the ground for unraveling the dynamics of the m6A epitranscriptome and the underlying mechanisms in clinical settings.

Prediction of LC-MS/MS Properties of Peptides from Sequence by Deep Learning.

Deep learning models for prediction of three key LC-MS/MS properties from peptide sequences were developed. The LC-MS/MS properties or behaviors are indexed retention times (iRT), MS1 or survey scan charge state distributions, and sequence ion intensities of HCD spectra. A common core deep supervised learning architecture, bidirectional long-short term memory (LSTM) recurrent neural networks was used to construct the three prediction models. Two featurization schemes were proposed and demonstrated to allow for efficient encoding of modifications. The iRT and charge state distribution models were trained with on order of 105 data points each. An HCD sequence ion prediction model was trained with 2 × 106 experimental spectra. The iRT prediction model and HCD sequence ion prediction model provide improved accuracies over the start-of-the-art models available in literature. The MS1 charge state distribution prediction model offers excellent performance. The prediction models can be used to enhance peptide identification and quantification in data-dependent acquisition and data-independent acquisition (DIA) experiments as well as to assist MRM (multiple reaction monitoring) and PRM (parallel reaction monitoring) experiment design.

Repeat-Preserving Decoy Database for False Discovery Rate Estimation in Peptide Identification.

The sequence database searching method is widely used in proteomics for peptide identification. To control the false discovery rate (FDR) of the searching results, the target-decoy method generates and searches a decoy database together with the target database. A known problem is that the target protein sequence database may contain numerous repeated peptides. The structures of these repeats are not preserved by most existing decoy generation algorithms. Previous studies suggest that such discrepancy between the target and decoy databases may lead to an inaccurate FDR estimation. Based on the de Bruijn graph model, we propose a new repeat-preserving algorithm to generate decoy databases. We prove that this algorithm preserves the structures of the repeats in the target database to a great extent. The de Bruijn method has been compared with a few other commonly used methods and demonstrated superior FDR estimation accuracy and an improved number of peptide identification.

Data Dependent-Independent Acquisition (DDIA) Proteomics.

Data dependent acquisition (DDA) and data independent acquisition (DIA) are traditionally separate experimental paradigms in bottom-up proteomics. In this work, we developed a strategy combining the two experimental methods into a single LC-MS/MS run. We call the novel strategy data dependent-independent acquisition proteomics, or DDIA for short. Peptides identified from DDA scans by a conventional and robust DDA identification workflow provide useful information for interrogation of DIA scans. Deep learning based LC-MS/MS property prediction tools, developed previously, can be used repeatedly to produce spectral libraries facilitating DIA scan extraction. A complete DDIA data processing pipeline, including the modules for iRT vs RT calibration curve generation, DIA extraction classifier training, and false discovery rate control, has been developed. Compared to another spectral library-free method, DIA-Umpire, the DDIA method produced a similar number of peptide identifications, but nearly twice as many protein group identifications. The primary advantage of the DDIA method is that it requires minimal information for processing its data.

Evosep One Enables Robust Deep Proteome Coverage Using Tandem Mass Tags while Significantly Reducing Instrument Time.

The balance between comprehensively analyzing the proteome and using valuable mass spectrometry time is a genuine challenge in the field of proteomics. Multidimensional fractionation strategies have significantly increased proteome coverage, but often at the cost of increased mass analysis time, despite advances in mass spectrometer acquisition rates. Recently, the Evosep One liquid chromatography system was shown to analyze peptide samples in a high-throughput manner without sacrificing in-depth proteomics coverage. We demonstrate the incorporation of Evosep One technology into our multiplexing workflow for analysis of tandem mass tag (TMT)-labeled nonsmall cell lung carcinoma (NSCLC) patient-derived xenografts (PDXs). By the use of a 30 samples per day Evosep workflow, >12 000 proteins were identified in 48 h of mass spectrometry time, which is comparable to the number of proteins identified by our conventional concatenated EASY-nLC workflow in 60 h. Shorter Evosep gradient lengths reduced the number of protein identifications by 10% while decreasing the mass analysis time by 50%. This Evosep workflow will enable quantitative analysis of multiplexed samples in less time without conceding depth of proteome coverage.

The ubiquitin ligase HERC4 mediates c-Maf ubiquitination and delays the growth of multiple myeloma xenografts in nude mice.

The transcription factor c-Maf is extensively involved in the pathophysiology of multiple myeloma (MM), a fatal malignancy of plasma cells. In the present study, affinity chromatography and mass spectrometry were used to identify c-Maf ubiquitination-associated proteins, from which the E3 ligase HERC4 was found to interact with c-Maf and catalyzed its polyubiquitination and subsequent proteasome-mediated degradation. HERC4 mediated polyubiquitination at K85 and K297 in c-Maf, and this polyubiquitination could be prevented by the isopeptidase USP5. Further analysis on the NCI-60 cell line collection revealed that RPMI 8226, a MM-derived cell line, expressed the lowest level of HERC4. Primary bone marrow analysis revealed HERC4 expression was high in normal bone marrow, but was steadily decreased during myelomagenesis. These findings suggested HERC4 played an important role in MM progression. Moreover, ectopic HERC4 expression decreased MM proliferation in vitro, and delayed xenograft tumor growth in vivo. Therefore, modulation of c-Maf ubiquitination by targeting HERC4 may represent a new therapeutic modality for MM.

Protein-phosphotyrosine proteome profiling by superbinder-SH2 domain affinity purification mass spectrometry, sSH2-AP-MS.

Recently, "superbinder" SH2 domain variants with three amino acid substitutions (sSH2) were reported to have 100-fold or greater affinity for protein-phosphotyrosine (pY) than natural SH2 domains. Here we report a protocol in which His-tagged Src sSH2 efficiently captures pY-peptides from protease-digested HeLa cell total protein extracts. Affinity purification of pY-peptides by this method shows little bias for pY-proximal amino acid sequences, comparable to that achieved by using antibodies to pY, but with equal or higher yield. Superbinder-SH2 affinity purification mass spectrometry (sSH2-AP-MS) therefore provides an efficient and economical approach for unbiased pY-directed phospho-proteome profiling without the use of antibodies.

Integrated analysis of proteome, phosphotyrosine-proteome, tyrosine-kinome, and tyrosine-phosphatome in acute myeloid leukemia.

Reversible protein-tyrosine phosphorylation is catalyzed by the antagonistic actions of protein-tyrosine kinases (PTKs) and phosphatases (PTPs), and represents a major form of cell regulation. Acute myeloid leukemia (AML) is an aggressive hematological malignancy that results from the acquisition of multiple genetic alterations, which in some instances are associated with deregulated protein-phosphotyrosine (pY) mediated signaling networks. However, although individual PTKs and PTPs have been linked to AML and other malignancies, analysis of protein-pY networks as a function of activated PTKs and PTPs has not been done. In this study, MS was used to characterize AML proteomes, and phospho-proteome-subsets including pY proteins, PTKs, and PTPs. AML proteomes resolved into two groups related to high or low degrees of maturation according to French-American-British classification, and reflecting differential expression of cell surface antigens. AML pY proteomes reflect canonical, spatially organized signaling networks, unrelated to maturation, with heterogeneous expression of activated receptor and nonreceptor PTKs. We present the first integrated analysis of the pY-proteome, activated PTKs, and PTPs. Every PTP and most PTKs have both positive and negative associations with the pY-proteome. pY proteins resolve into groups with shared PTK and PTP correlations. These findings highlight the importance of pY turnover and the PTP phosphatome in shaping the pY-proteome in AML.

Molecular heterogeneity of non-small cell lung carcinoma patient-derived xenografts closely reflect their primary tumors.

Availability of lung cancer models that closely mimic human tumors remains a significant gap in cancer research, as tumor cell lines and mouse models may not recapitulate the spectrum of lung cancer heterogeneity seen in patients. We aimed to establish a patient-derived tumor xenograft (PDX) resource from surgically resected non-small cell lung cancer (NSCLC). Fresh tumor tissue from surgical resection was implanted and grown in the subcutaneous pocket of non-obese severe combined immune deficient (NOD SCID) gamma mice. Subsequent passages were in NOD SCID mice. A subset of matched patient and PDX tumors and non-neoplastic lung tissues were profiled by whole exome sequencing, single nucleotide polymorphism (SNP) and methylation arrays, and phosphotyrosine (pY)-proteome by mass spectrometry. The data were compared to published NSCLC datasets of NSCLC primary and cell lines. 127 stable PDXs were established from 441 lung carcinomas representing all major histological subtypes: 52 adenocarcinomas, 62 squamous cell carcinomas, one adeno-squamous carcinoma, five sarcomatoid carcinomas, five large cell neuroendocrine carcinomas, and two small cell lung cancers. Somatic mutations, gene copy number and expression profiles, and pY-proteome landscape of 36 PDXs showed greater similarity with patient tumors than with established cell lines. Novel somatic mutations on cancer associated genes were identified but only in PDXs, likely due to selective clonal growth in the PDXs that allows detection of these low allelic frequency mutations. The results provide the strongest evidence yet that PDXs established from lung cancers closely mimic the characteristics of patient primary tumors.

Tyrosine phosphorylation of the Lyn Src homology 2 (SH2) domain modulates its binding affinity and specificity.

Src homology 2 (SH2) domains are modular protein structures that bind phosphotyrosine (pY)-containing polypeptides and regulate cellular functions through protein-protein interactions. Proteomics analysis showed that the SH2 domains of Src family kinases are themselves tyrosine phosphorylated in blood system cancers, including acute myeloid leukemia, chronic lymphocytic leukemia, and multiple myeloma. Using the Src family kinase Lyn SH2 domain as a model, we found that phosphorylation at the conserved SH2 domain residue Y(194) impacts the affinity and specificity of SH2 domain binding to pY-containing peptides and proteins. Analysis of the Lyn SH2 domain crystal structure supports a model wherein phosphorylation of Y(194) on the EF loop modulates the binding pocket that engages amino acid side chains at the pY+2/+3 position. These data indicate another level of regulation wherein SH2-mediated protein-protein interactions are modulated by SH2 kinases and phosphatases.

Comprehensive proteome analysis of fresh frozen and optimal cutting temperature (OCT) embedded primary non-small cell lung carcinoma by LC-MS/MS.

Clinical tissue samples provide valuable information for understanding human diseases. One major type of clinical tissue sample that is amenable to various kinds of analysis is fresh frozen and optimal cutting temperature (OCT)-embedded primary patient tissue. Recent advances in mass spectrometry (MS) technologies have been widely applied to study human proteomes by using clinical specimens. However, polymeric compounds such as OCT can interfere with MS analyses. Here we present methods that enable the preparation and analysis of fresh frozen and OCT embedded primary tissue samples by LC-MS/MS. A scraping method was first introduced to reduce the heterogeneity of OCT-embedded non-small cell lung carcinoma tumor sections. OCT compound was reproducibly removed by a series of washing steps involving ethanol and water prior to trypsin digestion. In data-dependent acquisition mode, optimized dynamic exclusion duration settings were established to maximize peptide identifications. These sample preparation conditions and MS parameter settings should be utilized or carefully adjusted in order to achieve optimal comprehensive proteome characterization starting from fresh frozen and OCT embedded clinical tissue specimens.

Proteomic analysis of the epidermal growth factor receptor (EGFR) interactome and post-translational modifications associated with receptor endocytosis in response to EGF and stress.

Aberrant expression, activation, and stabilization of epidermal growth factor receptor (EGFR) are causally associated with several human cancers. Post-translational modifications and protein-protein interactions directly modulate the signaling and trafficking of the EGFR. Activated EGFR is internalized by endocytosis and then either recycled back to the cell surface or degraded in the lysosome. EGFR internalization and recycling also occur in response to stresses that activate p38 MAP kinase. Mass spectrometry was applied to comprehensively analyze the phosphorylation, ubiquitination, and protein-protein interactions of wild type and endocytosis-defective EGFR variants before and after internalization in response to EGF ligand and stress. Prior to internalization, EGF-stimulated EGFR accumulated ubiquitin at 7 K residues and phosphorylation at 7 Y sites and at S(1104). Following internalization, these modifications diminished and there was an accumulation of S/T phosphorylations. EGFR internalization and many but not all of the EGF-induced S/T phosphorylations were also stimulated by anisomycin-induced cell stress, which was not associated with receptor ubiquitination or elevated Y phosphorylation. EGFR protein interactions were dramatically modulated by ligand, internalization, and stress. In response to EGF, different E3 ubiquitin ligases became maximally associated with EGFR before (CBL, HUWE1, and UBR4) or after (ITCH) internalization, whereas CBLB was distinctively most highly EGFR associated following anisomycin treatment. Adaptin subunits of AP-1 and AP-2 clathrin adaptor complexes also became EGFR associated in response to EGF and anisomycin stress. Mutations preventing EGFR phosphorylation at Y(998) or in the S(1039) region abolished or greatly reduced EGFR interactions with AP-2 and AP-1, and impaired receptor trafficking. These results provide new insight into spatial, temporal, and mechanistic aspects of EGFR regulation.

Comprehensive identification of phosphorylation sites on the Numb endocytic adaptor protein.

Numb is an adaptor protein that functions in the endocytosis and intracellular trafficking of membrane receptors and adhesion molecules. Previous studies have indicated that Numb localization and function are regulated through phosphorylation by atypical protein kinase C at several key sites. Here, using LC-MS/MS, we report the identification of 25 serine/threonine Numb phosphorylation sites, and a single tyrosine phosphorylation site. Amino acid sequences flanking several of the sites identified conform to consensus motifs for cyclin-dependent kinase 5 (CDK5). In vitro kinase assays and immunoblotting confirmed that CDK5 phosphorylates Numb. LC-MS/MS analysis identified specific CDK5-directed phosphorylation of Numb at position S288 and at two additional regions. Therefore, Numb is likely to exist in multiple phospho-isoforms, and may be subject to phosphorylation-mediated regulation downstream of CDK5. These findings provide a basis for further investigations into the complex role of Numb phosphorylation in regulating its subcellular localization, protein interactions, and function. All MS data have been deposited in the ProteomeXchange with identifier PXD000997 (http://proteomecentral.proteomexchange.org/dataset/PXD000997).

Differential regulation of FGFR3 by PTPN1 and PTPN2.

Aberrant expression and activation of FGFR3 is associated with disease states including bone dysplasia and malignancies of bladder, cervix, and bone marrow. MS analysis of protein-phosphotyrosine in multiple myeloma cells revealed a prevalent phosphorylated motif, D/EYYR/K, derived from the kinase domain activation loops of tyrosine kinases including FGFR3 corresponding to a recognition sequence of protein-tyrosine phosphatase PTPN1. Knockdown of PTPN1 or the related enzyme PTPN2 by RNAi resulted in ligand-independent activation of FGFR3. Modulation of FGFR3 activation loop phosphorylation by both PTPN1 and PTPN2 was a function of receptor trafficking and phosphotyrosine phosphatase (PTP) compartmentalization. The FGFR3 activation loop motif DYYKK(650) is altered to DYYKE(650) in the oncogenic variant FGFR3(K650E) , and consequently it is constitutively fully activated and unaffected by activation loop phosphorylation. FGFR3(K650E) was nevertheless remarkably sensitive to negative regulation by PTPN1 and PTPN2. This suggests that in addition to modulating FGFR3 phosphorylation, PTPN1 and PTPN2 constrain the kinase domain by fostering an inactive-state. Loss of this constraint in response to ligand or impaired PTPN1/N2 may initiate FGFR3 activation. These results suggest a model wherein PTP expression levels may define conditions that select for ectopic FGFR3 expression and activation during tumorigenesis.