MacroSEQUEST: efficient candidate-centric searching and high-resolution correlation analysis for large-scale proteomics data sets

Software Options for the Analysis of MS-Proteomic Data

Mass spectrometry (MS)-based proteomics is currently the most successful approach to measure and compare peptides and proteins in a large variety of biological samples. Modern mass spectrometers, equipped with high-resolution analyzers, provide large amounts of data output. This is the case of shotgun/bottom-up proteomics, which consists in the enzymatic digestion of protein into peptides that are then measured by MS-instruments through a data dependent acquisition (DDA) mode. Dedicated bioinformatic tools and platforms have been developed to face the increasing size and complexity of raw MS data that need to be processed and interpreted for large-scale protein identification and quantification. This chapter illustrates the most popular bioinformatics solution for the analysis of shotgun MS-proteomics data. A general description will be provided on the data preprocessing options and the different search engines available, including practical suggestions on how to optimize the parameters for peptide search, based on hands-on experience.

CHD7 and Runx1 interaction provides a braking mechanism for hematopoietic differentiation

Hematopoietic stem and progenitor cell (HSPC) formation and lineage differentiation involve gene expression programs orchestrated by transcription factors and epigenetic regulators. Genetic disruption of the chromatin remodeler chromodomain-helicase-DNA-binding protein 7 (CHD7) expanded phenotypic HSPCs, erythroid, and myeloid lineages in zebrafish and mouse embryos. CHD7 acts to suppress hematopoietic differentiation. Binding motifs for RUNX and other hematopoietic transcription factors are enriched at sites occupied by CHD7, and decreased RUNX1 occupancy correlated with loss of CHD7 localization. CHD7 physically interacts with RUNX1 and suppresses RUNX1-induced expansion of HSPCs during development through modulation of RUNX1 activity. Consequently, the RUNX1:CHD7 axis provides proper timing and function of HSPCs as they emerge during hematopoietic development or mature in adults, representing a distinct and evolutionarily conserved control mechanism to ensure accurate hematopoietic lineage differentiation.

Anatomy and evolution of database search engines-a central component of mass spectrometry based proteomic workflows

Sequence database search engines are bioinformatics algorithms that identify peptides from tandem mass spectra using a reference protein sequence database. Two decades of development, notably driven by advances in mass spectrometry, have provided scientists with more than 30 published search engines, each with its own properties. In this review, we present the common paradigm behind the different implementations, and its limitations for modern mass spectrometry datasets. We also detail how the search engines attempt to alleviate these limitations, and provide an overview of the different software frameworks available to the researcher. Finally, we highlight alternative approaches for the identification of proteomic mass spectrometry datasets, either as a replacement for, or as a complement to, sequence database search engines.

Direct regulation of Chk1 protein stability by E3 ubiquitin ligase HUWE1

The HECT E3 ubiquitin ligase HUWE1 is required for a wide array of important functions in cell biology. Although HUWE1 is known to play a role in DNA damage signaling, the mechanism(s) that underlie this function remain elusive. HUWE1 regulates effectors of DNA replication and genotoxic stress tolerance. However, the loss of HUWE1 can also result in the accrual of significant endogenous DNA damage due to insufficient remediation of replication stress induced by an overabundance of key substrates. We discovered that HUWE1 depletion leads to a significant increase in levels of the single-strand break effector kinase Chk1, independent of the DNA damage response, activation of apical DNA damage repair (DDR) signaling kinases (ATM and ATR), and the tumor suppressor p53. We also identified multiple lysine residues on Chk1 that are polyubiquitinated by HUWE1 in vitro, many of which are within the kinase domain. HUWE1 knockdown also markedly prolonged the protein half-life of Chk1 in steady-state conditions and resulted in greater stabilization of Chk1 protein than depletion of Cul4A, an E3 ubiquitin ligase previously described to control Chk1 abundance. Moreover, prolonged replication stress induced by hydroxyurea or camptothecin resulted in a reduction of Chk1 protein levels, which was rescued by HUWE1 knockdown. Our study indicates that HUWE1 plays a significant role in the regulation of the DDR signaling pathway by directly modulating the abundance of Chk1 protein.

Tempest: Accelerated MS/MS Database Search Software for Heterogeneous Computing Platforms

MS/MS database search algorithms derive a set of candidate peptide sequences from in silico digest of a protein sequence database, and compute theoretical fragmentation patterns to match these candidates against observed MS/MS spectra. The original Tempest publication described these operations mapped to a CPU-GPU model, in which the CPU (central processing unit) generates peptide candidates that are asynchronously sent to a discrete GPU (graphics processing unit) to be scored against experimental spectra in parallel. The current version of Tempest expands this model, incorporating OpenCL to offer seamless parallelization across multicore CPUs, GPUs, integrated graphics chips, and general-purpose coprocessors. Three protocols describe how to configure and run a Tempest search, including discussion of how to leverage Tempest's unique feature set to produce optimal results. © 2016 by John Wiley & Sons, Inc.

Visualizing the Mechanism of Epoxide Hydrolysis by the Bacterial Virulence Enzyme Cif

The CFTR inhibitory factor (Cif) is an epoxide hydrolase (EH) virulence factor secreted by the bacterium Pseudomonas aeruginosa. Sequence alignments reveal a pattern of Cif-like substitutions that proved to be characteristic of a new subfamily of bacterial EHs. At the same time, crystallographic and mutagenetic data suggest that EH activity is required for virulence and that Cif's active site remains generally compatible with a canonical two-step EH mechanism. A hallmark of this mechanism is the formation of a covalent hydroxyalkyl-enzyme intermediate by nucleophilic attack. In several well-studied EHs, this intermediate has been captured at near stoichiometric levels, presumably reflecting rate-limiting hydrolysis. Here we show by mass spectrometry that only minimal levels of the expected intermediate can be trapped with WT Cif. In contrast, substantial amounts of intermediate are recovered from an active-site mutant (Cif-E153Q) that selectively targets the second, hydrolytic release step. Utilizing Cif-E153Q and a previously reported nucleophile mutant (Cif-D129S), we then captured Cif in the substrate-bound, hydroxyalkyl-intermediate, and product-bound states for 1,2-epoxyhexane, yielding the first crystallographic snapshots of an EH at these key stages along the reaction coordinate. Taken together, our data illuminate the proposed two-step hydrolytic mechanism of a new class of bacterial virulence factor. They also suggest that the failure of WT Cif to accumulate a covalent hydroxyalkyl-enzyme intermediate reflects an active-site chemistry in which hydrolysis is no longer the rate-limiting step, a noncanonical kinetic regime that may explain similar observations with a number of other EHs.

The SEQUEST family tree

Since its introduction in 1994, SEQUEST has gained many important new capabilities, and a host of successor algorithms have built upon its successes. This Account and Perspective maps the evolution of this important tool and charts the relationships among contributions to the SEQUEST legacy. Many of the changes represented improvements in computing speed by clusters and graphics cards. Mass spectrometry innovations in mass accuracy and activation methods led to shifts in fragment modeling and scoring strategies. These changes, as well as the movement of laboratories and lab members, have led to great diversity among the members of the SEQUEST family. Graphical Abstract ᅟ.

A deeper look into Comet--implementation and features

The Comet database search software was initially released as an open source project in late 2012. Prior to that, Comet existed as the University of Washington's academic version of the SEQUEST database search tool. Despite its availability and widespread use over the years, some details about its implementation have not been previously disseminated or are not well understood. We address a few of these details in depth and highlight new features available in the latest release. Comet is freely available for download at http://comet-ms.sourceforge.net or it can be accessed as a component of a number of larger software projects into which it has been incorporated. Graphical Abstract ᅟ.

Quantitative phosphoproteomics reveals new roles for the protein phosphatase PP6 in mitotic cells

Protein phosphorylation is an important regulatory mechanism controlling mitotic progression. Protein phosphatase 6 (PP6) is an essential enzyme with conserved roles in chromosome segregation and spindle assembly from yeast to humans. We applied a baculovirus-mediated gene silencing approach to deplete HeLa cells of the catalytic subunit of PP6 (PP6c) and analyzed changes in the phosphoproteome and proteome in mitotic cells by quantitative mass spectrometry-based proteomics. We identified 408 phosphopeptides on 272 proteins that increased and 298 phosphopeptides on 220 proteins that decreased in phosphorylation upon PP6c depletion in mitotic cells. Motif analysis of the phosphorylated sites combined with bioinformatics pathway analysis revealed previously unknown PP6c-dependent regulatory pathways. Biochemical assays demonstrated that PP6c opposed casein kinase 2-dependent phosphorylation of the condensin I subunit NCAP-G, and cellular analysis showed that depletion of PP6c resulted in defects in chromosome condensation and segregation in anaphase, consistent with dysregulation of condensin I function in the absence of PP6 activity.

LRRK2 Facilitates tau Phosphorylation through Strong Interaction with tau and cdk5

Leucine-rich repeat kinase 2 (LRRK2) and tau have been identified as risk factors of Parkinson's disease (PD). As LRRK2 is a kinase and tau is hyperphosphorylated in some LRRK2 mutation carriers of PD patients, the obvious hypothesis is that tau could be a substrate of LRRK2. Previous reports that LRRK2 phosphorylates free tau or tubulin-associated tau provide direct support for this proposition. By comparing LRRK2 with cdk5, we show that wild-type LRRK2 and the G2019S mutant phosphorylate free recombinant full-length tau protein with specific activity 480- and 250-fold lower than cdk5, respectively. More strikingly tau binds to wt LRRK2 or the G2019S mutant 140- or 200-fold more strongly than cdk5. The extremely low activity of LRRK2 but strong binding affinity with tau suggests that LRRK2 may facilitate tau phosphorylation as a scaffold protein rather than as a major tau kinase. This hypothesis is further supported by the observation that (i) cdk5 or tau coimmunoprecipitates with endogenous LRRK2 in SH-SY5Y cells, in mouse brain tissue, and in human PBMCs; (ii) knocking down endogenous LRRK2 by its siRNA in SH-SY5Y cells reduces tau phosphorylation at Ser396 and Ser404; (iii) inhibiting LRRK2 kinase activity by its inhibitors has no effect on tau phosphorylation at these two sites; and (iv) overexpressing wt LRRK2, the G2019S mutant, or the D1994A kinase-dead mutant in SH-SY5Y cells has no effect on tau phosphorylation. Our results suggest that LRRK2 facilitates tau phosphorylation indirectly by recruiting tau or cdk5 rather than by directly phosphorylating tau.

Quantitative phosphoproteomics reveals pathways for coordination of cell growth and division by the conserved fission yeast kinase pom1

Complex phosphorylation-dependent signaling networks underlie the coordination of cellular growth and division. In the fission yeast Schizosaccharomyces pombe, the Dual specificity tyrosine-(Y)-phosphorylation regulated kinase (DYRK) family protein kinase Pom1 regulates cell cycle progression through the mitotic inducer Cdr2 and controls cell polarity through unknown targets. Here, we sought to determine the phosphorylation targets of Pom1 kinase activity by SILAC-based phosphoproteomics. We defined a set of high-confidence Pom1 targets that were enriched for cytoskeletal and cell growth functions. Cdr2 was the only cell cycle target of Pom1 kinase activity that we identified in cells. Mutation of Pom1-dependent phosphorylation sites in the C terminus of Cdr2 inhibited mitotic entry but did not impair Cdr2 localization. In addition, we found that Pom1 phosphorylated multiple substrates that function in polarized cell growth, including Tea4, Mod5, Pal1, the Rho GAP Rga7, and the Arf GEF Syt22. Purified Pom1 phosphorylated these cell polarity targets in vitro, confirming that they are direct substrates of Pom1 kinase activity and likely contribute to regulation of polarized growth by Pom1. Our study demonstrates that Pom1 acts in a linear pathway to control cell cycle progression while regulating a complex network of cell growth targets.

SPECHT - single-stage phosphopeptide enrichment and stable-isotope chemical tagging: quantitative phosphoproteomics of insulin action in muscle

The study of cellular signaling remains a significant challenge for translational and clinical research. In particular, robust and accurate methods for quantitative phosphoproteomics in tissues and tumors represent significant hurdles for such efforts. In the present work, we design, implement and validate a method for single-stage phosphopeptide enrichment and stable isotope chemical tagging, or SPECHT, that enables the use of iTRAQ, TMT and/or reductive dimethyl-labeling strategies to be applied to phosphoproteomics experiments performed on primary tissue. We develop and validate our approach using reductive dimethyl-labeling and HeLa cells in culture, and find these results indistinguishable from data generated from more traditional SILAC-labeled HeLa cells mixed at the cell level. We apply the SPECHT approach to the quantitative analysis of insulin signaling in a murine myotube cell line and muscle tissue, identify known as well as new phosphorylation events, and validate these phosphorylation sites using phospho-specific antibodies. Taken together, our work validates chemical tagging post-single-stage phosphoenrichment as a general strategy for studying cellular signaling in primary tissues.

BIOLOGICAL SIGNIFICANCE

Through the use of a quantitatively reproducible, proteome-wide phosphopeptide enrichment strategy, we demonstrated the feasibility of post-phosphopeptide purification chemical labeling and tagging as an enabling approach for quantitative phosphoproteomics of primary tissues. Using reductive dimethyl labeling as a generalized chemical tagging strategy, we compared the performance of post-phosphopeptide purification chemical tagging to the well established community standard, SILAC, in insulin-stimulated tissue culture cells. We then extended our method to the analysis of low-dose insulin signaling in murine muscle tissue, and report on the analytical and biological significance of our results.

SILAC surrogates: rescue of quantitative information for orphan analytes in spike-in SILAC experiments

Super-stable isotope labeling by amino acids in cell culture (Super-SILAC) enables the sensitive and accurate analysis of complex biological tissue and tumor samples by comparison of light peptides observed in biological samples to heavy peptides from SILAC cell culture spike-ins. However, despite the use of multiple cell lines for Super-SILAC spike-in standards, the full protein and peptide profiles of biological samples are not completely represented in these internal standards, leading to orphan analytes for which sample to standard ratios cannot be calculated. This problem is exacerbated in some biological systems, such as muscle tissue, which lack adequate cell culture lines to reflect their complex and idiosyncratic protein profiles, resulting in up to 40% of peptide analytes without heavy cognates. Furthermore, these unquantified orphan analytes may be among the most biologically interesting and significant species, since their presence is not common to cell lines cultured in vitro. Here, we report on the development of a surrogate analysis strategy to interpolate quantitative relationships between peptide species, observed across multiple biological samples, which lack representation within the spike-in standards. The precision and accuracy of this method was assessed by replicate experiments in which surrogate-derived ratios from defined mixtures of spike-in SILAC standard and tissue lysate were compared against traditional SILAC ratios for species where both light and heavy peptide cognates were observed. We demonstrate the robustness of our SILAC surrogates strategy across a variety of murine tissues, including liver, spleen, brain, and muscle. Our approach increases the quantitative coverage and precision within a biological sample by rescuing previously intractable peptide species and applying additional evidence to improve the precision of existing quantifications.

Quantitative phosphoproteomic profiling of human non-small cell lung cancer tumors

Non-small cell lung cancer (NSCLC) is the leading cause of cancer-related deaths worldwide. Within the molecular scope of NCSLC, a complex landscape of dysregulated cellular signaling has emerged, defined largely by mutations in select mediators of signal transduction, including the epidermal growth factor receptor (EGFR) and anaplastic lymphoma (ALK) kinases. Consequently, these mutant kinases become constitutively activated and targets for chemotherapeutic intervention. Encouragingly, small molecule inhibitors of these pathways have shown promise in clinical trials or are approved for clinical use. However, many protein kinases are dysregulated in NSCLC without genetic mutations. To quantify differences in tumor cell signaling that are transparent to genomic methods, we established a super-SILAC internal standard derived from NSCLC cell lines grown in vitro and labeled with heavy lysine and arginine, and deployed them in a phosphoproteomic workflow. We identified 9019 and 8753 phosphorylation sites in two separate tumors. Relative quantification of phosphopeptide abundance between tumor samples allowed for the determination of specific hubs and pathways differing between each tumor. Sites downstream of Ras showed decreased inhibitory phosphorylation (Raf/Mek) and increased activating phosphorylation (Erk1/2) in one tumor versus another. In this way, we were able to quantitatively access oncogenic kinase signaling in primary human tumors.

BIOLOGICAL SIGNIFICANCE

Through the use of quantitative proteomics, we demonstrated the feasibility and coverage that large scale mass spectrometry can leverage for understanding kinase networks in cancer. By incorporating Super-SILAC based quantitation into a typical pathology workflow, we were able to access and compare tumors from multiple patients in this analysis with high accuracy and dynamic range. We analyzed tumors from patients diagnosed with non-small cell lung cancer and were able to detect comprehensive phosphorylation networks relaying through known hubs of oncogenesis in lung cancer. We hereby show that it is possible to track changes to phosphorylation networks across multiple tumors, opening up the possibility that drug susceptibility and patient-specific stratification can be implemented downstream of classical pathology.

Thioredoxin reductase is inhibited by the carbamoylating activity of the anticancer sulfonylhydrazine drug laromustine

The thioredoxin system facilitates proliferative processes in cells and is upregulated in many cancers. The activities of both thioredoxin (Trx) and its reductase (TrxR) are mediated by oxidation/reduction reactions among cysteine residues. A common target in preclinical anticancer research, TrxR is reported here to be significantly inhibited by the anticancer agent laromustine. This agent, which has been in clinical trials for acute myelogenous leukemia and glioblastoma multiforme, is understood to be cytotoxic principally via interstrand DNA crosslinking that originates from a 2-chloroethylating species generated upon activation in situ. The spontaneous decomposition of laromustine also yields methyl isocyanate, which readily carbamoylates thiols and primary amines. Purified rat liver TrxR was inhibited by laromustine with a clinically relevant IC(50) value of 4.65 μM. A derivative of laromustine that lacks carbamoylating activity did not appreciably inhibit TrxR while another derivative, lacking only the 2-chloroethylating activity, retained its inhibitory potency. Furthermore, in assays measuring TrxR activity in murine cell lysates, a similar pattern of inhibition among these compounds was observed. These data contrast with previous studies demonstrating that glutathione reductase, another enzyme that relies on cysteine-mediated redox chemistry, was not inhibited by methylcarbamoylating agents when measured in cell lysates. Mass spectrometry of laromustine-treated enzyme revealed significant carbamoylation of TrxR, albeit not on known catalytically active residues. However, there was no evidence of 2-chloroethylation anywhere on the protein. The inhibition of TrxR is likely to contribute to the cytotoxic, anticancer mechanism of action for laromustine.

Tempest: GPU-CPU computing for high-throughput database spectral matching

Modern mass spectrometers are now capable of producing hundreds of thousands of tandem (MS/MS) spectra per experiment, making the translation of these fragmentation spectra into peptide matches a common bottleneck in proteomics research. When coupled with experimental designs that enrich for post-translational modifications such as phosphorylation and/or include isotopically labeled amino acids for quantification, additional burdens are placed on this computational infrastructure by shotgun sequencing. To address this issue, we have developed a new database searching program that utilizes the massively parallel compute capabilities of a graphical processing unit (GPU) to produce peptide spectral matches in a very high throughput fashion. Our program, named Tempest, combines efficient database digestion and MS/MS spectral indexing on a CPU with fast similarity scoring on a GPU. In our implementation, the entire similarity score, including the generation of full theoretical peptide candidate fragmentation spectra and its comparison to experimental spectra, is conducted on the GPU. Although Tempest uses the classical SEQUEST XCorr score as a primary metric for evaluating similarity for spectra collected at unit resolution, we have developed a new "Accelerated Score" for MS/MS spectra collected at high resolution that is based on a computationally inexpensive dot product but exhibits scoring accuracy similar to that of the classical XCorr. In our experience, Tempest provides compute-cluster level performance in an affordable desktop computer.

Plk1 regulates the kinesin-13 protein Kif2b to promote faithful chromosome segregation

Solid tumors are frequently aneuploid, and many display high rates of ongoing chromosome missegregation in a phenomenon called chromosomal instability (CIN). The most common cause of CIN is the persistence of aberrant kinetochore-microtubule (k-MT) attachments, which manifest as lagging chromosomes in anaphase. k-MT attachment errors form during prometaphase due to stochastic interactions between kinetochores and microtubules. The kinesin-13 protein Kif2b promotes the correction of k-MT attachment errors in prometaphase, but the mechanism restricting this activity to prometaphase remains unknown. Using mass spectrometry, we identified multiple phosphorylation sites on Kif2b, some of which are acutely sensitive to inhibition of Polo-like kinase 1 (Plk1). We show that Plk1 directly phosphorylates Kif2b at threonine 125 (T125) and serine 204 (S204), and that these two sites differentially regulate Kif2b function. Phosphorylation of S204 is required for the kinetochore localization and activity of Kif2b in prometaphase, and phosphorylation of T125 is required for Kif2b activity in the correction of k-MT attachment errors. These data demonstrate that Plk1 regulates both the localization and activity of Kif2b during mitosis to promote the correction of k-MT attachment errors to ensure mitotic fidelity.

Translational research in infectious disease: current paradigms and challenges ahead

In recent years, the biomedical community has witnessed a rapid scientific and technologic evolution after the development and refinement of high-throughput methodologies. Concurrently and consequentially, the scientific perspective has changed from the reductionist approach of meticulously analyzing the fine details of a single component of biology to the "holistic" approach of broadmindedly examining the globally interacting elements of biological systems. The emergence of this new way of thinking has brought about a scientific revolution in which genomics, proteomics, metabolomics, and other "omics" have become the predominant tools by which large amounts of data are amassed, analyzed, and applied to complex questions of biology that were previously unsolvable. This enormous transformation of basic science research and the ensuing plethora of promising data, especially in the realm of human health and disease, have unfortunately not been followed by a parallel increase in the clinical application of this information. On the contrary, the number of new potential drugs in development has been decreasing steadily, suggesting the existence of roadblocks that prevent the translation of promising research into medically relevant therapeutic or diagnostic application. In this article, we will review, in a noninclusive fashion, several recent scientific advancements in the field of translational research, with a specific focus on how they relate to infectious disease. We will also present a current picture of the limitations and challenges that exist for translational research, as well as ways that have been proposed by the National Institutes of Health to improve the state of this field.

Rapid determination of multiple linear kinase substrate motifs by mass spectrometry

Kinase-substrate recognition depends on the chemical properties of the phosphorylatable residue as well as the surrounding linear sequence motif. Detailed knowledge of these characteristics increases the confidence of linking identified phosphorylation sites to kinases, predicting phosphorylation sites, and designing optimal peptide substrates. Here, we present a mass spectrometry-based approach for determining linear kinase substrate motifs by elaborating the positional and chemical preference of the kinase for a phosphorylatable residue using libraries of naturally-occurring peptides that are amenable to peptide identification by commonly used proteomics platforms. We applied this approach to a structurally and functionally diverse set of purified kinases, which recapitulated their previously described substrate motifs and discovered additional ones, including preferences of certain kinases for phosphorylatable residues adjacent to peptide termini. Furthermore, we identify specific and distinguishable motif elements for the four members of the polo-like kinase (Plk) family and verify members of these motif elements for Plk1 in vivo.

The Tlo proteins are stoichiometric components of Candida albicans mediator anchored via the Med3 subunit

The amplification of the TLO (for telomere-associated) genes in Candida albicans, compared to its less pathogenic, close relative Candida dubliniensis, suggests a role in virulence. Little, however, is known about the function of the Tlo proteins. We have purified the Mediator coactivator complex from C. albicans (caMediator) and found that Tlo proteins are a stoichiometric component of caMediator. Many members of the Tlo family are expressed, and each is a unique member of caMediator. Protein expression analysis of individual Tlo proteins, as well as the purification of tagged Tlo proteins, demonstrate that there is a large free population of Tlo proteins in addition to the Mediator-associated population. Coexpression and copurification of Tloα12 and caMed3 in Escherichia coli established a direct physical interaction between the two proteins. We have also made a C. albicans med3Δ/Δ strain and purified an intact Mediator from this strain. The analysis of the composition of the med3Δ Mediator shows that it lacks a Tlo subunit. Regarding Mediator function, the med3Δ/Δ strain serves as a substitute for the difficult-to-make tloΔ/Δ C. albicans strain. A potential role of the TLO and MED3 genes in virulence is supported by the inability of the med3Δ/Δ strain to form normal germ tubes. This study of caMediator structure provides initial clues to the mechanism of action of the Tlo genes and a platform for further mechanistic studies of caMediator's involvement in gene regulatory patterns that underlie pathogenesis.

Search and decoy: the automatic identification of mass spectra

In recent years, the generation and interpretation of MS/MS spectra for the identification of peptides and proteins has matured to a frequently used automatic workflow in Proteomics. Several software solutions for the automated analysis of MS/MS spectra allow for high-throughput/high-performance analyses of complex samples. Related to MS/MS searches, target-decoy approaches have gained more and more popularity: in a "decoy" part of the search database nonexistent sequences mimic real sequences (the "target" sequences). With their help, the number of falsely identified peptides/proteins can be estimated after a search and the resulting protein list can be cut at a specified false discovery rate (FDR). This is an essential prerequisite for all quantitative approaches, as they rely on correct identifications. Especially the label-free approach "spectral counting"-gaining more and more popularity due to low costs and simplicity-depends directly on the correctness of peptide-spectrum matches (PSMs). This work's aim is to describe five popular search engines-especially their general properties regarding protein identification, but also their quantification abilities, if those go beyond spectral counting. By doing so, Proteomics researchers are enabled to compare their features and to choose an appropriate solution for their specific question. Furthermore, the search engines are applied to a spectrum data set generated from a complex sample with a Thermo LTQ Velos OrbiTrap (Thermo Fisher Scientific, Waltham, MA, USA). The results of the search engines are compared, e.g., regarding time requirements, peptides and proteins found, and the search engines' behavior using the decoy approach.

Increasing phosphoproteomic coverage through sequential digestion by complementary proteases

Protein phosphorylation is a reversible post-translational modification known to regulate protein function, subcellular localization, complex formation, and protein degradation. Detailed phosphoproteomic information is critical to kinomic studies of signal transduction and for elucidation of cancer biomarkers, such as in non-small-cell lung adenocarcinoma, where phosphorylation is commonly dysregulated. However, the collection and analysis of phosphorylation data remains a difficult problem. The low concentrations of phosphopeptides in complex biological mixtures as well as challenges inherent in their chemical nature have limited phosphoproteomic characterization and some phosphorylation sites are inaccessible by traditional workflows. We developed a sequential digestion method using complementary proteases, Glu-C and trypsin, to increase phosphoproteomic coverage and supplement traditional approaches. The sequential digestion method is more productive than workflows utilizing only Glu-C and we evaluated the orthogonality of the sequential digestion method relative to replicate trypsin-based analyses. Finally, we demonstrate the ability of the sequential digestion method to access new regions of the phosphoproteome by comparison to existing public phosphoproteomic databases. Our approach increases coverage of the human lung cancer phosphoproteome by accessing both new phosphoproteins and novel phosphorylation site information.