Nature Biotechnology: doi: /nbt Supplementary Figure 1. The workflow of Open-pFind.

Transcription

1 Supplementary Figure 1 The workflow of Open-pFind. The MS data are first preprocessed by pparse, and then the MS/MS data are searched by the open search module. Next, the MS/MS data are re-searched by the restricted search module against a refined search space based on the learned information in the reranking step. Finally, the results obtained from both the open and restricted searches are merged, reranked again and reported.

2 Supplementary Figure 2 An overview of the Open-pFind results for the Dong-Ecoli-QE dataset. This figure is from a screenshot of the embedded software tool pbuild, which was designed to graphically represent pfind Studio results. a) Distributions of mass deviations (in ppm) of the identified PSMs in the five raw files. The number of PSMs for each raw (N) is 16,007 (10mM), 13,418 (30mM), 16,931 (60mM), 18,575 (150mM3) and 24,523 (1000mM), respectively. Box-plot elements: center line, median; box limits, first and third quartile (Q1 and Q3); whiskers, from Q1 1.5 IQR to Q3+1.5 IQR. b) Score distributions of the target (blue) and decoy (red) PSMs. c) Distribution of the enzymatic specificity of identified peptides. d) Percentages of the highly abundant modifications in all identified peptides. e) Distribution of the number of missed cleavage sites. f) Distribution of the number of peptides identified from one tandem mass spectrum, e.g., two peptides can be identified from 12.65% of the total spectra.

3 Supplementary Figure 3 Precision evaluation of search engines using the Dong-Ecoli-QE dataset. Metabolically labeled datasets are searched against the protein database, in which only the unlabeled peptides are considered, and pquant is then independently used for check the quantitation ratio of each PSM.

4 Supplementary Figure 4 Percentages of NaN-ratio PSMs identified in the Dong-Ecoli-QE dataset (I). For the open search engines, only PSMs with no or common modifications specified in the restricted search engines were considered in the analysis. a) Consistently identified PSMs from the comparison of every two search engines are considered. b) Separately identified PSMs from the comparison of every two search engines are considered. The 15 N- and unlabeled peptides are used to calculate the quantitation values in both a) and b). c) Similar to a), but the 13 C- and unlabeled peptides are used to calculate the quantitation values. d) Similar to b), but the 13 C- and unlabeled peptides are used to calculate the quantitation values.

5 Supplementary Figure 5 Percentages of NaN-ratio PSMs identified in the Xu-Yeast-QEHF dataset. For the open search engines, only PSMs with no or common modifications specified in the restricted search engines were considered in the analysis. a) Consistently identified PSMs from the comparison of every two search engines are considered. b) Separately identified PSMs from the comparison of every two search engines are considered.

6 Supplementary Figure 6 Percentages of NaN-ratio PSMs identified in the Dong-Ecoli-QE dataset (II). All PSMs (both modified and unmodified ones) were considered in the analysis. a) Consistently identified PSMs from the comparison of every two search engines are considered. b) Separately identified PSMs from the comparison of every two search engines are considered. The 15 N- and unlabeled peptides are used to calculate the quantitation values in both a) and b). c) Similar to a), but the 13 C- and unlabeled peptides are used to calculate the quantitation values. d) Similar to b), but the 13 C- and unlabeled peptides are used to calculate the quantitation values.

7 Supplementary Figure 7 Analysis of NaN-ratio proteins in the Dong-Ecoli-QE dataset. Open-pFind is compared with each of the other seven engines. a) The number of proteins that are consistently identified by two engines (&), separately identified by Open-pFind (+) or the other engine ( ), and the percentages of NaN-ratio proteins in those three sets (&, + and, respectively). The quantitation value of a protein is the NaN ratio if no supporting protein-unique peptides are assigned with normal quantitation values. b) Similar to a), but proteins supported with at least two protein-unique peptides are considered. The quantitation value of a protein is the NaN ratio if the number of the supporting protein-unique peptides with normal quantitation values is less than two. 15 N-labeled peptides were used to calculate the quantitation values in both a) and b). c) Similar to a), but 13 C-labeled peptides were used for calculating the quantitation values. d) Similar to b), but 13 C-labeled peptides were used for calculating the quantitation values.

8 Supplementary Figure 8 Performance evaluation of the four open search engines using the entrapment strategy for the four published datasets. The datasets were chosen as shown in Supplementary Table 1. Each database search process used the target database (human or mouse) appended by the entrapment database (Arabidopsis thaliana). All databases were downloaded from UniProt ( and only the reviewed proteins were used. The other database search parameters were the same as those shown in Supplementary Table 3. In each subfigure, the x-axis denotes the number of total PSMs reported by the open search engines, and the y-axis denotes the fraction of the entrapment PSMs in the total PSMs, e.g., for the Mann-Human-Velos dataset, Open-pFind reported 23,876 PSMs at 1% FDR at the peptide level, in which 0.25% of the spectra were matched with peptides from the entrapment database.

9 Supplementary Figure 9 The distribution of the results reported by Open-pFind that were not included in the traditional search space. For example, among all distinct peptides identified by Open-pFind for Mann-Human-Velos dataset, the semi-tryptic/non-specific peptides, peptides with unexpected modifications including mutations, peptides identified with additional precursors exported by pparse, accounted for 4.3%, 17.1% and 27.7%, respectively, and a total of 42.9% peptides fall outside of the restricted search space.

10 Supplementary Figure 10 The search spaces for different search modes. The four modes includes the restricted search with no modifications (No mod), with common modifications (Common, with carbamidomethylation of C, oxidation of M, Gln pyro-glu at N-termini of peptides and acetylation at N-termini of proteins), and with phosphorylation (Phospho, with common modifications and phosphorylation of S, T, and Y), and the open search (Unimod, with at most one modification in Unimod), together with three different types of digestion. For each mode, 1,000 experiments were performed, each with 1,000 randomly selected proteins from UniProt database, which were digested in silico into peptides. Then, the distribution of the number of peptides from each experiment is shown in one box-plot (N = 1,000). Box-plot elements: center line, median; box limits, first and third quartile (Q1 and Q3); whiskers, from Q1 1.5 IQR to Q3+1.5 IQR.

11 Supplementary Figure 11 Relationship between the protein database reduction and the running time of Open-pFind. #PR denotes the number of proteins in the reduced protein database after open search, i.e., the first step of Open-pFind. #PO denotes the number of proteins in the original protein database specified before database searches. The red curve denotes the ratio of the running time of Open-pFind to the running time of pfind in the analyses of the six datasets. Each bar denotes the mean and each error bar (blue whisker) denotes the standard deviation (SD) measured among different datasets (N = 5, 22, 3, 24, 4, 24 for the number of raw files in the six datasets shown from left to right, respectively). Each grey dot denotes a data point in the corresponding dataset. For example, in the analysis of the Xu-Yeast-QEHF dataset, ~97% of the proteins in the original database were removed after the open search step on average; as a result, the speed of Open-pFind was approximately twice that of pfind.

12 Supplementary Figure 12 The comparison between the picked TDA approach and the common two-peptide TDA approach for protein inference with the Kim data. The number of proteins reported by the picked strategy was 16,133, which was larger than the number that was reported using the twopeptide rule (14,064), but seven additional olfactory receptors were detected with the picked strategy.

13 Supplementary Figure 13 An example of a mixed spectrum with three co-eluting peptides. The spectrum is entitled as Ecoli-1to1to1-un-C13-N15-60mM dta, in which three co-eluting peptides were identified. a) MS1 information for the three peptides in the isolation window around this tandem mass spectrum. b) The match between the spectrum and the three resulting peptides, taken from a screenshot of pbuild. Three isotopic clusters are clearly distinguished from each other, while the observed fragmentation sites were near-complete for all three peptides.

14 Supplementary Figure 14 Similarity between the actual and predicted fragment ions for each of the three co-eluting peptides identified from one spectrum. The spectrum is entitled as Ecoli-1to1to1-un-C13-N15-60mM dta. The predicted spectra are from pdeep algorithm. The similarity was measured with the Pearson's correlation coefficients (PCC): a) TLAEGQNVEFEIQDGQK, b) SEYLGDPDFVK and c) VSLAADPVEEIK. The number of singly charged b and y ions (N), for which the correlations were determined, is 32, 20 and 22 for a) c), respectively.