Integrative Proteogenomics for Differential Expression and Splicing Variation in a DM1 Mouse Model.

Dysregulated mRNA splicing is involved in the pathogenesis of many diseases including cancer, neurodegenerative diseases, and muscular dystrophies such as myotonic dystrophy type 1 (DM1). Comprehensive assessment of dysregulated splicing on the transcriptome and proteome level has been methodologically challenging, and thus investigations have often been targeting only few genes. Here, we performed a large-scale coordinated transcriptomic and proteomic analysis to characterize a DM1 mouse model (HSALR) in comparison to wild type. Our integrative proteogenomics approach comprised gene- and splicing-level assessments for mRNAs and proteins. It recapitulated many known instances of aberrant mRNA splicing in DM1 and identified new ones. It enabled the design and targeting of splicing-specific peptides and confirmed the translation of known instances of aberrantly spliced disease-related genes (e.g., Atp2a1, Bin1, Ryr1), complemented by novel findings (Flnc and Ywhae). Comparative analysis of large-scale mRNA and protein expression data showed quantitative agreement of differentially expressed genes and splicing patterns between disease and wild type. We hence propose this work as a suitable blueprint for a robust and scalable integrative proteogenomic strategy geared toward advancing our understanding of splicing-based disorders. With such a strategy, splicing-based biomarker candidates emerge as an attractive and accessible option, as they can be efficiently asserted on the mRNA and protein level in coordinated fashion.

Deuterium induces a distinctive Escherichia coli proteome that correlates with the reduction in growth rate.

Substitution of protium (H) for deuterium (D) strongly affects biological systems. Whereas higher eukaryotes such as plants and mammals hardly survive a deuterium content of >30%, many microorganisms can grow on fully deuterated media, albeit at reduced rates. Very little is known about how the H/D replacement influences life at the systems level. Here, we used MS-based analysis to follow the adaptation of a large part of the Escherichia coli proteome from growth on a protonated full medium, over a protonated minimal medium, to a completely deuterated minimal medium. We could quantify >1800 proteins under all conditions, several 100 of which exhibited strong regulation during both adaptation processes. The adaptation to minimal medium strongly up-regulated amino acid synthesis and sugar metabolism and down-regulated translational proteins on average by 9%, concomitant with a reduction in growth rate from 1.8 to 0.67 h-1 In contrast, deuteration caused a very wide proteomic response over many cell functional categories, together with an additional down-regulation of the translational proteins by 5%. The latter coincided with a further reduction in growth rate to 0.37 h-1, revealing a clear linear correlation between growth rate and abundance of translational proteins. No significant morphological effects are observed under light and electron microscopies. Across all protein categories, about 80% of the proteins up-regulated under deuteration are enzymes with hydrogen transfer functions. Thus, the H/D kinetic isotope effect appears as the major limiting factor of cellular functions under deuteration.

Global Ion Suppression Limits the Potential of Mass Spectrometry Based Phosphoproteomics.

Mass spectrometry based proteomics has become the method of choice for pinpointing and monitoring thousands of post-translational modifications, predominately phosphorylation sites, in cellular signaling studies. Critical for achieving this analytical depth is the enrichment of phosphorylated peptides prior to liquid chromatography-mass spectrometry (MS) analysis. Despite the high prevalence of this modification, the numbers of identified phosphopeptides lag behind those achieved for unmodified peptides, and the cause for this still remains controversial. Here, we use an effective phosphatase protocol that considerably improves global ionization efficiency and, therefore, the overall sensitivity and coverage of standard phosphoproteomics studies. We demonstrate the power of our method on the model system of Salmonella-infected macrophages by extending the current quantitative picture of immune signaling pathways involved in infection. In combination with sensitive, label-free targeted MS for phosphorylation site validation, our approach is ideally suited to exploring cellular phosphorylation based signaling networks in high detail.

Evaluation and Improvement of Quantification Accuracy in Isobaric Mass Tag-Based Protein Quantification Experiments.

The multiplexing capabilities of isobaric mass tag-based protein quantification, such as Tandem Mass Tags or Isobaric Tag for Relative and Absolute Quantitation have dramatically increased the scope of mass spectrometry-based proteomics studies. Not only does the technology allow for the simultaneous quantification of multiple samples in a single MS injection, but its seamless compatibility with extensive sample prefractionation methods allows for comprehensive studies of complex proteomes. However, reporter ion-based quantification has often been criticized for limited quantification accuracy due to interference from coeluting peptides and peptide fragments. In this study, we investigate the extent of this problem and propose an effective and easy-to-implement remedy that relies on spiking a 6-protein calibration mixture to the samples. We evaluated our ratio adjustment approach using two large scale TMT 10-plex data sets derived from a human cancer and noncancer cell line as well as E. coli cells grown at two different conditions. Furthermore, we analyzed a complex 2-proteome artificial sample mixture and investigated the precision of TMT and precursor ion intensity-based label free quantification. Studying the protein set identified by both methods, we found that differentially abundant proteins were assigned dramatically higher statistical significance when quantified using TMT. Data are available via ProteomeXchange with identifier PXD003346.

Exploiting the multiplexing capabilities of tandem mass tags for high-throughput estimation of cellular protein abundances by mass spectrometry.

The generation of dynamic models of biological processes critically depends on the determination of precise cellular concentrations of biomolecules. Measurements of system-wide absolute protein levels are particularly valuable information in systems biology. Recently, mass spectrometry based proteomics approaches have been developed to estimate protein concentrations on a proteome-wide scale. However, for very complex proteomes, fractionation steps are required, increasing samples number and instrument analysis time. As a result, the number of full proteomes that can be routinely analyzed is limited. Here we combined absolute quantification strategies with the multiplexing capabilities of isobaric tandem mass tags to determine cellular protein abundances in a high throughput and proteome-wide scale even for highly complex biological systems, such as a whole human cell line. We generated two independent data sets to demonstrate the power of the approach regarding sample throughput, dynamic range, quantitative precision and accuracy as well as proteome coverage in comparison to existing mass spectrometry based strategies.

A mammalian transcription factor-specific peptide repository for targeted proteomics.

Site-specific transcription factors (TFs) play an essential role in mammalian development and function as they are vital for the majority of cellular processes. Despite their biological importance, TF proteomic data is scarce in the literature, likely due to difficulties in detecting peptides as the abundance of TFs in cells tends to be low. In recent years, significant improvements in MS-based technologies in terms of sensitivity and specificity have increased the interest in developing quantitative methodologies specifically targeting relatively lowly abundant proteins such as TFs in mammalian models. Such efforts would be greatly aided by the availability of TF peptide-specific information as such data would not only enable improvements in speed and accuracy of protein identifications, but also ameliorate cross-comparisons of quantitative proteomics data and allow for a more efficient development of targeted proteomics assays. However, to date, no comprehensive TF proteotypic peptide database has been developed. To address this evident lack of TF peptide data in public repositories, we are generating a comprehensive, experimentally derived TF proteotypic peptide spectral library dataset based on in vitro protein expression. Our library currently contains peptide information for 89 TFs and this number is set to increase in the near future. All MS data have been deposited in the ProteomeXchange with identifier PXD001212 (http://proteomecentral.proteomexchange.org/dataset/PXD001212).

Comparison of Different Sample Preparation Protocols Reveals Lysis Buffer-Specific Extraction Biases in Gram-Negative Bacteria and Human Cells.

We evaluated different in-solution and FASP-based sample preparation strategies for absolute protein quantification. Label-free quantification (LFQ) was employed to compare different sample preparation strategies in the bacterium Pseudomonas aeruginosa and human embryonic kidney cells (HEK), and organismal-specific differences in general performance and enrichment of specific protein classes were noted. The original FASP protocol globally enriched for most proteins in the bacterial sample, whereas the sodium deoxycholate in-solution strategy was more efficient with HEK cells. Although detergents were found to be highly suited for global proteome analysis, higher intensities were obtained for high-abundant nucleic acid-associated protein complexes, like the ribosome and histone proteins, using guanidine hydrochloride. Importantly, we show for the first time that the observable total proteome mass of a sample strongly depends on the sample preparation protocol, with some protocols resulting in a significant underestimation of protein mass due to incomplete protein extraction of biased protein groups. Furthermore, we demonstrate that some of the observed abundance biases can be overcome by incorporating a nuclease treatment step or, alternatively, a correction factor for complementary sample preparation approaches.

Systems-level overview of host protein phosphorylation during Shigella flexneri infection revealed by phosphoproteomics.

The enteroinvasive bacterium Shigella flexneri invades the intestinal epithelium of humans. During infection, several injected effector proteins promote bacterial internalization, and interfere with multiple host cell responses. To obtain a systems-level overview of host signaling during infection, we analyzed the global dynamics of protein phosphorylation by liquid chromatography-tandem MS and identified several hundred of proteins undergoing a phosphorylation change during the first hours of infection. Functional bioinformatic analysis revealed that they were mostly related to the cytoskeleton, transcription, signal transduction, and cell cycle. Fuzzy c-means clustering identified six temporal profiles of phosphorylation and a functional module composed of ATM-phosphorylated proteins related to genotoxic stress. Pathway enrichment analysis defined mTOR as the most overrepresented pathway. We showed that mTOR complex 1 and 2 were required for S6 kinase and AKT activation, respectively. Comparison with a published phosphoproteome of Salmonella typhimurium-infected cells revealed a large subset of coregulated phosphoproteins. Finally, we showed that S. flexneri effector OspF affected the phosphorylation of several hundred proteins, thereby demonstrating the wide-reaching impact of a single bacterial effector on the host signaling network.

Extensive mass spectrometry-based analysis of the fission yeast proteome: the Schizosaccharomyces pombe PeptideAtlas.

We report a high quality and system-wide proteome catalogue covering 71% (3,542 proteins) of the predicted genes of fission yeast, Schizosaccharomyces pombe, presenting the largest protein dataset to date for this important model organism. We obtained this high proteome and peptide (11.4 peptides/protein) coverage by a combination of extensive sample fractionation, high resolution Orbitrap mass spectrometry, and combined database searching using the iProphet software as part of the Trans-Proteomics Pipeline. All raw and processed data are made accessible in the S. pombe PeptideAtlas. The identified proteins showed no biases in functional properties and allowed global estimation of protein abundances. The high coverage of the PeptideAtlas allowed correlation with transcriptomic data in a system-wide manner indicating that post-transcriptional processes control the levels of at least half of all identified proteins. Interestingly, the correlation was not equally tight for all functional categories ranging from r(s) >0.80 for proteins involved in translation to r(s) <0.45 for signal transduction proteins. Moreover, many proteins involved in DNA damage repair could not be detected in the PeptideAtlas despite their high mRNA levels, strengthening the translation-on-demand hypothesis for members of this protein class. In summary, the extensive and publicly available S. pombe PeptideAtlas together with the generated proteotypic peptide spectral library will be a useful resource for future targeted, in-depth, and quantitative proteomic studies on this microorganism.

Evaluation of data-dependent and -independent mass spectrometric workflows for sensitive quantification of proteins and phosphorylation sites.

In recent years, directed and, particularly, targeted mass spectrometric workflows have gained momentum as alternative techniques to conventional data-dependent acquisition (DDA) LC-MS/MS approaches. By focusing on specific peptide species, these methods allow hypothesis-driven analysis of selected proteins of interest, and they have been shown to be suited to monitor low-abundance proteins within complex mixtures. Despite their growing popularity, no study has systematically evaluated these various MS strategies in terms of quantification, detection, and identification limits when they are applied to complex samples. Here, we systematically compared the performance of conventional DDA, directed, and various targeted MS approaches on two different instruments, namely, a hybrid linear ion trap--Orbitrap and a triple quadrupole instrument. We assessed the limits of identification, quantification, and detection for each method by analyzing a dilution series of 20 unmodified and 10 phosphorylated synthetic heavy-labeled reference peptides, respectively, covering 6 orders of magnitude in peptide concentration with and without a complex human cell digest background. We found that all methods performed similarly in the absence of background proteins; however, when analyzing whole-cell lysates, targeted methods were at least 5-10 times more sensitive than that of the directed or DDA method. In particular, higher stage fragmentation (MS3) of the neutral loss peak using a linear ion trap increased the dynamic quantification range of some phosphopeptides up to 100-fold. We illustrate the power of this targeted MS3 approach for phosphopeptide monitoring by successfully quantifying nine phosphorylation sites of the kinetochore and spindle assembly checkpoint component Mad1 over different cell cycle states from nonenriched pull-down samples.

Small molecule induced STING degradation facilitated by the HECT ligase HERC4.

Stimulator of interferon genes (STING) is a central component of the cytosolic nucleic acids sensing pathway and as such master regulator of the type I interferon response. Due to its critical role in physiology and its' involvement in a variety of diseases, STING has been a focus for drug discovery. Targeted protein degradation (TPD) has emerged as a promising pharmacology for targeting previously considered undruggable proteins by hijacking the cellular ubiquitin proteasome system (UPS) with small molecules. Here, we identify AK59 as a STING degrader leveraging HERC4, a HECT-domain E3 ligase. Additionally, our data reveals that AK59 is effective on the common pathological STING mutations, suggesting a potential clinical application of this mechanism. Thus, these findings introduce HERC4 to the fields of TPD and of compound-induced degradation of STING, suggesting potential therapeutic applications.

DCAF1-based PROTACs with activity against clinically validated targets overcoming intrinsic- and acquired-degrader resistance.

Targeted protein degradation (TPD) mediates protein level through small molecule induced redirection of E3 ligases to ubiquitinate neo-substrates and mark them for proteasomal degradation. TPD has recently emerged as a key modality in drug discovery. So far only a few ligases have been utilized for TPD. Interestingly, the workhorse ligase CRBN has been observed to be downregulated in settings of resistance to immunomodulatory inhibitory drugs (IMiDs). Here we show that the essential E3 ligase receptor DCAF1 can be harnessed for TPD utilizing a selective, non-covalent DCAF1 binder. We confirm that this binder can be functionalized into an efficient DCAF1-BRD9 PROTAC. Chemical and genetic rescue experiments validate specific degradation via the CRL4DCAF1 E3 ligase. Additionally, a dasatinib-based DCAF1 PROTAC successfully degrades cytosolic and membrane-bound tyrosine kinases. A potent and selective DCAF1-BTK-PROTAC (DBt-10) degrades BTK in cells with acquired resistance to CRBN-BTK-PROTACs while the DCAF1-BRD9 PROTAC (DBr-1) provides an alternative strategy to tackle intrinsic resistance to VHL-degrader, highlighting DCAF1-PROTACS as a promising strategy to overcome ligase mediated resistance in clinical settings.

Critical assessment of proteome-wide label-free absolute abundance estimation strategies.

There is a great interest in reliable ways to obtain absolute protein abundances at a proteome-wide scale. To this end, label-free LC-MS/MS quantification methods have been proposed where all identified proteins are assigned an estimated abundance. Several variants of this quantification approach have been presented, based on either the number of spectral counts per protein or MS1 peak intensities. Equipped with several datasets representing real biological environments, containing a high number of accurately quantified reference proteins, we evaluate five popular low-cost and easily implemented quantification methods (Absolute Protein Expression, Exponentially Modified Protein Abundance Index, Intensity-Based Absolute Quantification Index, Top3, and MeanInt). Our results demonstrate considerably improved abundance estimates upon implementing accurately quantified reference proteins; that is, using spiked in stable isotope labeled standard peptides or a standard protein mix, to generate a properly calibrated quantification model. We show that only the Top3 method is directly proportional to protein abundance over the full quantification range and is the preferred method in the absence of reference protein measurements. Additionally, we demonstrate that spectral count based quantification methods are associated with higher errors than MS1 peak intensity based methods. Furthermore, we investigate the impact of miscleaved, modified, and shared peptides as well as protein size and the number of employed reference proteins on quantification accuracy.

Large-scale quantitative assessment of different in-solution protein digestion protocols reveals superior cleavage efficiency of tandem Lys-C/trypsin proteolysis over trypsin digestion.

The complete and specific proteolytic cleavage of protein samples into peptides is crucial for the success of every shotgun LC-MS/MS experiment. In particular, popular peptide-based label-free and targeted mass spectrometry approaches rely on efficient generation of fully cleaved peptides to ensure accurate and sensitive protein quantification. In contrast to previous studies, we globally and quantitatively assessed the efficiency of different digestion strategies using a yeast cell lysate, label-free quantification, and statistical analysis. Digestion conditions include double tryptic, surfactant-assisted, and tandem-combinatorial Lys-C/trypsin digestion. In comparison to tryptic digests, Lys-C/trypsin digests were found most efficient to yield fully cleaved peptides while reducing the abundance of miscleaved peptides. Subsequent sequence context analysis revealed improved digestion performances of Lys-C/trypsin for miscleaved sequence stretches flanked by charged basic and particulary acidic residues. Furthermore, targeted MS analysis demonstrated a more comprehensive protein cleavage only after Lys-C/trypsin digestion, resulting in a more accurrate absolute protein quantification and extending the number of peptides suitable for SRM assay development. Therefore, we conclude that a serial Lys-C/trypsin digestion is highly attractive for most applications in quantitative MS-based proteomics building on in-solution digestion schemes.

An improved method for the construction of decoy peptide MS/MS spectra suitable for the accurate estimation of false discovery rates.

The relevance of libraries of annotated MS/MS spectra is growing with the amount of proteomic data generated in high-throughput experiments. These reference libraries provide a fast and accurate way to identify newly acquired MS/MS spectra. In the context of multiple hypotheses testing, the control of the number of false-positive identifications expected in the final result list by means of the calculation of the false discovery rate (FDR). In a classical sequence search where experimental MS/MS spectra are compared with the theoretical peptide spectra calculated from a sequence database, the FDR is estimated by searching randomized or decoy sequence databases. Despite on-going discussion on how exactly the FDR has to be calculated, this method is widely accepted in the proteomic community. Recently, similar approaches to control the FDR of spectrum library searches were discussed. We present in this paper a detailed analysis of the similarity between spectra of distinct peptides to set the basis of our own solution for decoy library creation (DeLiberator). It differs from the previously published results in some key points, mainly in implementing new methods that prevent decoy spectra from being too similar to the original library spectra while keeping important features of real MS/MS spectra. Using different proteomic data sets and library creation methods, we evaluate our approach and compare it with alternative methods.

QuickMod: A tool for open modification spectrum library searches.

MS2 library spectra are rich in reproducible information about peptide fragmentation patterns compared to theoretical spectra modeled by a sequence search tool. So far, spectrum library searches are mostly applied to detect peptides as they are present in the library. However, they also allow finding modified variants of the library peptides if the search is done with a large precursor mass window and an adapted Spectrum-Spectrum Match (SSM) scoring algorithm. We perform a thorough evaluation on the use of library spectra as opposed to theoretical peptide spectra for the identification of PTMs, analyzing spectra of a well-annotated modification-rich test data set compiled from public data repositories. These initial studies motivate the development of our modification tolerant spectrum library search tool QuickMod, designed to identify modified variants of the peptides listed in the spectrum library without any prior input from the user estimating the modifications present in the sample. We built the search algorithm of QuickMod after carefully testing different SSM similarity scores. The final spectrum scoring scheme uses a support vector machine (SVM) on a selection of scoring features to classify correct and incorrect SSM. After identification of a list of modified peptides at a given False Discovery Rate (FDR), the modifications need to be positioned on the peptide sequence. We present a rapid modification site assignment algorithm and evaluate its positioning accuracy. Finally, we demonstrate that QuickMod performs favorably in terms of speed and identification rate when compared to other software solutions for PTM analysis.

PAX8 and MECOM are interaction partners driving ovarian cancer.

The transcription factor PAX8 is critical for the development of the thyroid and urogenital system. Comprehensive genomic screens furthermore indicate an additional oncogenic role for PAX8 in renal and ovarian cancers. While a plethora of PAX8-regulated genes in different contexts have been proposed, we still lack a mechanistic understanding of how PAX8 engages molecular complexes to drive disease-relevant oncogenic transcriptional programs. Here we show that protein isoforms originating from the MECOM locus form a complex with PAX8. These include MDS1-EVI1 (also called PRDM3) for which we map its interaction with PAX8 in vitro and in vivo. We show that PAX8 binds a large number of genomic sites and forms transcriptional hubs. At a subset of these, PAX8 together with PRDM3 regulates a specific gene expression module involved in adhesion and extracellular matrix. This gene module correlates with PAX8 and MECOM expression in large scale profiling of cell lines, patient-derived xenografts (PDXs) and clinical cases and stratifies gynecological cancer cases with worse prognosis. PRDM3 is amplified in ovarian cancers and we show that the MECOM locus and PAX8 sustain in vivo tumor growth, further supporting that the identified function of the MECOM locus underlies PAX8-driven oncogenic functions in ovarian cancer.

Unrestricted identification of modified proteins using MS/MS.

Proteins undergo PTM, which modulates their structure and regulates their function. Estimates of the PTM occurrence vary but it is safe to assume that there is an important gap between what is currently known and what remains to be discovered. The highest throughput and most comprehensive efforts to catalogue protein mixtures have so far been using MS-based shotgun proteomics. The standard approach to analyse MS/MS data is to use Peptide Fragment Fingerprinting tools such as Sequest, MASCOT or Phenyx. These tools commonly identify 5-30% of the spectra in an MS/MS data set while only a limited list of predefined protein modifications can be screened. An important part of the unidentified spectra is likely to be spectra of peptides carrying modifications not considered in the search. Bioinformatics for PTM discovery is an active area of research. In this review we focus on software solutions developed for unrestricted identification of modifications in MS/MS data, here referred to as open modification search tools. We give an overview of the conceptually different algorithmic solutions to evaluate the large number of candidate peptides per spectrum when accounting for modifications of unrestricted size and demonstrate the value of results of large-scale open modification search studies. Efficient and easy-to-use tools for protein modification discovery should prove valuable in the quest for mapping the dynamics of proteomes.

A simple workflow to increase MS2 identification rate by subsequent spectral library search.

Searching a spectral library for the identification of protein MS/MS data has proven to be a fast and accurate method, while yielding a high identification rate. We investigated the potential to increase peptide discovery rate, with little increase in computational time, by constructing a workflow based on a sequence search with Phenyx followed by a library search with SpectraST. Searching a consensus library compiled from the search results of the prior Phenyx search increased the number of confidently matched spectra by up to 156%. Additionally matched spectra by SpectraST included noisy spectra, spectra representing missed cleaved peptides as well as spectra from post-translationally modified peptides.