## Abstract

In large-scale proteomic experiments, multiple peptide precursors are often cofragmented simultaneously in the same *mixture* tandem mass (MS/MS) spectrum. These spectra tend to elude current computational tools because of the ubiquitous assumption that each spectrum is generated from only one peptide. Therefore, tools that consider multiple peptide matches to each MS/MS spectrum can potentially improve the relatively low spectrum identification rate often observed in proteomics experiments. More importantly, data independent acquisition protocols *promoting* the cofragmentation of multiple precursors are emerging as alternative methods that can greatly improve the throughput of peptide identifications but their success also depends on the availability of algorithms to identify multiple peptides from each MS/MS spectrum. Here we address a fundamental question in the identification of mixture MS/MS spectra: determining the statistical significance of multiple peptides matched to a given MS/MS spectrum. We propose the MixGF generating function model to rigorously compute the statistical significance of peptide identifications for mixture spectra and show that this approach improves the sensitivity of current mixture spectra database search tools by a ≈30–390%. Analysis of multiple data sets with MixGF reveals that in complex biological samples the number of identified mixture spectra can be as high as 20% of all the identified spectra and the number of unique peptides identified only in mixture spectra can be up to 35.4% of those identified in single-peptide spectra.

The advancement of technology and instrumentation has made tandem mass (MS/MS)^{1} spectrometry the leading high-throughput method to analyze proteins (1, 2, 3). In typical experiments, tens of thousands to millions of MS/MS spectra are generated and enable researchers to probe various aspects of the proteome on a large scale. Part of this success hinges on the availability of computational methods that can analyze the large amount of data generated from these experiments. The classical question in computational proteomics asks: given an MS/MS spectrum, what is the *peptide* that generated the spectrum? However, it is increasingly being recognized that this assumption that each MS/MS spectrum comes from *only one* peptide is often not valid. Several recent analyses show that as many as 50% of the MS/MS spectra collected in typical proteomics experiments come from more than one peptide precursor (4, 5). The presence of multiple peptides in mixture spectra can decrease their identification rate to as low as one half of that for MS/MS spectra generated from only one peptide (6, 7, 8). In addition, there have been numerous developments in data independent acquisition (DIA) technologies where multiple peptide precursors are intentionally selected to cofragment in each MS/MS spectrum (9, 10, 11, 12, 13, 14, 15). These emerging technologies can address some of the enduring disadvantages of traditional data-dependent acquisition (DDA) methods (*e.g.* low reproducibility (16)) and potentially increase the throughput of peptide identification 5–10 fold (4, 17). However, despite the growing importance of mixture spectra in various contexts, there are still only a few computational tools that can analyze mixture spectra from more than one peptide (18, 19, 20, 21, 8, 22). Our recent analysis indicated that current database search methods for mixture spectra still have relatively low sensitivity compared with their single-peptide counterpart and the main bottleneck is their limited ability to separate true matches from false positive matches (8). Traditionally problem of peptide identification from MS/MS spectra involves two sub-problems: 1) define a Peptide-Spectrum-Match (PSM) scoring function that assigns each MS/MS spectrum to the peptide sequence that most likely generated the spectrum; and 2) given a set of top-scoring PSMs, select a subset that corresponds to statistical significance PSMs. Here we focus on the second problem, which is still an ongoing research question even for the case of single-peptide spectra (23, 24, 25, 26). Intuitively the second problem is difficult because one needs to consider spectra across the whole data set (instead of comparing different peptide candidates against one spectrum as in the first problem) and PSM scoring functions are often not well-calibrated across different spectra (*i.e.* a PSM score of 50 may be good for one spectrum but poor for a different spectrum). Ideally, a scoring function will give high scores to all true PSMs and low scores to false PSMs regardless of the peptide or spectrum being considered. However, in practice, some spectra may receive higher scores than others simply because they have more peaks or their precursor mass results in more peptide candidates being considered from the sequence database (27, 28). Therefore, a scoring function that accounts for spectrum or peptide-specific effects can make the scores more comparable and thus help assess the confidence of identifications across different spectra. The MS-GF solution to this problem is to compute the per-spectrum statistical significance of each top-scoring PSM, which can be defined as the probability that a random peptide (out of all possible peptide within parent mass tolerance) will match to the spectrum with a score at least as high as that of the top-scoring PSM. This measures how good the current best match is in relation to all possible peptides matching to the same spectrum, normalizing any spectrum effect from the scoring function. Intuitively, our proposed MixGF approach extends the MS-GF approach to now calculate the statistical significance of the top *pair* of peptides matched from the database to a given mixture spectrum M (*i.e.* the significance of the top peptide–peptide spectrum match (PPSM)). As such, MixGF determines the probability that a random pair of peptides (out of all possible peptides within parent mass tolerance) will match a given mixture spectrum with a score at least as high as that of the top-scoring PPSM.

Despite the theoretical attractiveness of computing statistical significance, it is generally prohibitive for any database search methods to score all possible peptides against a spectrum. Therefore, earlier works in this direction focus on approximating this probability by assuming the score distribution of all PSMs follows certain analytical form such as the normal, Poisson or hypergeometric distributions (29, 30, 31). In practice, because score distributions are highly data-dependent and spectrum-specific, these model assumptions do not always hold. Other approaches tried to learn the score distribution empirically from the data (29, 27). However, one is most interested in the region of the score distribution where only a small fraction of false positives are allowed (typically at 1% FDR). This usually corresponds to the extreme tail of the distribution where *p* values are on the order of 10^{−9} or lower and thus there is typically lack of sufficient data points to accurately model the tail of the score distribution (32). More recently, Kim *et al.* (24) and Alves *et al.* (33), in parallel, proposed a generating function approach to compute the exact score distribution of random peptide matches for any spectra without explicitly matching all peptides to a spectrum. Because it is an exact computation, no assumption is made about the form of score distribution and the tail of the distribution can be computed very accurately. As a result, this approach substantially improved the ability to separate true matches from false positive ones and lead to a significant increase in sensitivity of peptide identification over state-of-the-art database search tools in single-peptide spectra (24).

For mixture spectra, it is expected that the scores for the top-scoring match will be even less comparable across different spectra because now more than one peptide and different numbers of peptides can be matched to each spectrum at the same time. We extend the generating function approach (24) to rigorously compute the statistical significance of multiple-Peptide-Spectrum Matches (mPSMs) and demonstrate its utility toward addressing the peptide identification problem in mixture spectra. In particular, we show how to extend the generating approach for mixture from two peptides. We focus on this relatively simple case of mixture spectra because it accounts for a large fraction of mixture spectra presented in traditional DDA workflows (5). This allows us to test and develop algorithmic concepts using readily-available DDA data because data with more complex mixture spectra such as those from DIA workflows (11) is still not widely available in public repositories.

## MATERIALS AND METHODS

##### Spectral Probability for a Mixture Spectrum

For single-peptide spectra, the statistical significance of a particular peptide *P* matched to a spectrum *S* with score *T* is determined by the probability that a random peptide *R* (out of all possible peptides) when matched to S has a score greater or equal to *T*: Pr(*Score*(*R*, *S*)≥*T*) where *Score*(*R*, *S*) is a scoring function for a peptide-spectrum-match. From here on, we will refer to this as the *Single-peptide probability* in order to distinguish it from the other definitions introduced below. Analogously, to compute the statistical significance of a particular peptide pair (*P*, *Q*) matched to a mixture spectrum (*M*) with a score of *T*, we are interested in two statistical questions: 1) *Joint probability* ≡ Pr(*Score*(*R*1,*R*2,*M*) ≥*T*): the probability that a random peptide pair (*R*_{1}, *R*_{2}) (out of all possible peptide pairs) when matched to *M* yields a score greater or equal to *T* and 2) *Conditional probability* ≡ Pr(*Score*(*R*1,*R*2,*M*) ≥*T* | *R*1 = *P*): given a peptide *P*, the probability that a random peptide *R*_{2} (out of all possible peptides) together with *P* when matched to *M* yields a score greater or equal to *T*. Intuitively a peptide–peptide spectrum match (PPSM) can fall into three categories: (1) *Correct-match*: both peptides are correct matches; (2) *Half-correct match*: one peptide is correct and the other peptide is an incorrect match; and 3) *Incorrect-match*: both peptides are incorrect matches. We are interested in separating the correct matches from incorrect and half-correct matches. The definitions above address this question in two steps. The joint probability assesses the chance that two random peptides have the same or higher score than a given match. When this probability is very low, this means that at least one peptide is a statistically significant match to the spectrum (*i.e.* it is a correct or half-correct match). Once we assume that at least one peptide is a true match, the conditional probability assesses whether the second peptide is also a statistically significant match (*i.e.* correct matches). In summary, one is looking for PPSMs with both low joint probability and conditional probability.

##### Scoring Function for Mixture Spectrum

Before describing the computation of the different probability measures, first we review the basics of our scoring function for a peptide-spectrum-match. We represent a tandem mass (MS/MS) spectrum with parent mass N as a real-valued vector: *V* = ν* _{1}* … ν

*with*

_{N}*N*elements, where ν

*is the sum of intensity of all the peaks with mass between*

_{i}*i*− 0.5 and

*i*+ 0.5 and parent mass is defined as the sum of the masses of all amino acids in the peptide that generated the spectrum. A prefix residue mass (PRM) spectrum is a transformation of an MS/MS spectrum into a scored version

*S*=

*s*

_{1}…

*s*using a probabilistic model as described before (34). In brief, at every mass position

_{N}*i*of the PRM spectrum is a score

*s*that represents the log-likelihood that the peptide from which the spectrum was generated contains a prefix mass

_{i}*i*(35). Given a peptide

*P*, its prefix masses are defined by the amino acid masses for each peptide prefix. For a peptide

*P*of length

*n*with prefix masses

*p*

_{1}…

*p*, we define its parent mass as

_{n}*p*and the score of matching peptide

_{n}*P*to a spectrum is the sum of all the scores at its theoretical prefix masses in the PRM spectrum:

Note that the probabilistic model used to generate the PRM spectrum depends on the precursor charge state of the MS/MS spectrum (34), thus when matching *P* to *S*, the precursor charge state for *S* is determined such that the parent mass of *P* is equal to that of *S* within the specified mass error tolerance.

We define a mixture spectrum as a spectrum from two different peptides. When interpreting an MS/MS spectrum as a mixture spectrum *M*, we construct two PRM spectra, *M ^{H}* and

*M*, each generated using the corresponding scoring models for high and low-abundance peptides present in a mixture spectrum. As shown in MixDB (8), different scoring models are needed for high and low-abundance peptides because they exhibit substantially different fragmentation statistics in mixture spectra. Intuitively, this is because the low-abundance peptides will generate less intense peaks in the mixture spectrum and, in general, it also has less number of detectable peaks above noise level. For example, the median peak intensity rank (ranked by decreasing peak intensity) for a y-ion from high-abundance peptides is 19, whereas the median peak rank for a y-ion from low-abundance peptide drops to 35. Without loss of generality, when matching a mixture spectrum (

^{L}*M*) against a pair of peptides (

*P*,

*Q*) we assume that the first peptide (

*P*) is the high-abundance peptide. Thus, the score of a pair of peptides (

*P*,

*Q*) against a mixture spectrum

*M*will be the sum of scoring

*P*with

*M*and scoring

^{H}*Q*with

*M*:

^{L}To avoid double counting, when a prefix mass of *P* is the same as a prefix mass of *Q*, only the bin with the higher score is considered and the other peptide gets a score of zero for that particular mass position:

##### Computing Spectral Probabilities

In order to compute the probabilities mentioned above we need to know the score distribution for all possible peptides and peptide pairs. The original MS-GF (24) approach uses dynamic programming to efficiently compute the single-peptide probability without explicitly considering the scores for all peptides. Here we extend this generating function approach to compute the distributions for the joint and conditional probabilities. Let *J _{M}* be a three-dimensional dynamic programming matrix where each element J

_{M}(

*p, q, T*) stands for the joint probability that a pair of peptides

*P*,

*Q*with parent mass

*p*and

*q*match to

*M*with score higher than or equal to

*T*. This means

*P*matches to

*M*up to the

^{H}*p*-th bin and

*Q*matches to

*M*up to

^{L}*q*-th bin. The following recurrence can then be used to compute the joint probability:

In the equation above *a*, *a*_{1,} and *a*_{2} denote amino acids; *mass*(*a*) denotes the mass of an amino acid; *prob*(*a*) denotes the probability that a particular amino acid occurs in a peptide and recall that *M ^{H}* and

*M*are the PRM spectra defined in the previous section. When considering all possible peptide sequences this probability is uniform and has a value of for each of the 20 standard amino acids. To better reflect the amino acid composition observed in real protein sequences we can also define this probability by computing the frequency of each amino acid in the protein sequence database against which the spectra are searched. To start the computation of the recurrence, we initialize

^{L}*J*(0, 0, 0) = 1 and

_{m}*J*(

_{M}*p,q,s*) for all entries where

*p*or

*q*is smaller than the smallest mass of an amino acid or

*s*is less than zero.

The computation of the conditional probability is similar to that of single-peptide probability, except that it is conditioned on the first peptide being accepted as a match. Specifically, for a peptide pair (*P*, *Q*) matched to a spectrum *M* with score *T*, we define that peptide *P* and *Q* contribute *T _{p}* and

*T*to the total score, respectively. Assuming that peptide

_{Q}*P*was matched to

*M*, we define a two-dimensional dynamic programming matrix

*C*where each element

_{M}*C*(

_{M}*q*,

*T*|

*P*) represents the conditional probability that a peptide with parent mass

*q*together with

*P*match

*M*with a score greater than or equal to

*T*. To compute this probability, we first modify

*M*by setting all the bins corresponding to a prefix mass of

^{L}*P*to zero if

*M*has a higher score at the same location. Then Conditional probability can be computed using the following recursion:

^{H}We initialize the recurrence with the base case: *C _{M}*(0,

*T*|

_{P}*P*). The base case starts at score

*T*rather than zero because the first peptide

_{P}*P*already contributes

*T*to the total score.

_{P}We note that even though the joint probability assesses whether at least one peptide is a significant match to the spectrum, it does not determine which peptide is the significant match in the case when only one peptide is a significant match. More importantly, when calculating the conditional probability one assumes that the first peptide is a true match but it is unclear which peptide is the first peptide from the joint probability assessment. In order to resolve this ambiguity, for a candidate peptide pair (*P*, *Q*) matched to a spectrum *M*, we compute their respective single-peptide probabilities and the peptide with lower (*i.e.* statistically more significant) single-peptide probability is designated as the first peptide. The dynamic programming method described above assumes that peptide fragment ions have integer masses. However, this is not appropriate for data sets with high mass accuracy in the MS/MS spectra. The details of how to extend this method for high mass accuracy data are described in the Supplementary Material. The current implementation of mixgf considers the set of all unmodified peptides or peptide pairs when computing the conditional and joint probability, however as shown in unpublished work MSGF+ (36) it is possible to extend this approach to take variable modifications into consideration.

##### Approximating Joint Probability

The dynamic programming approach described above enables the computation of the Joint probability without explicitly computing the scores for all peptide pairs. However, the computational complexity still scales exponentially with the number of peptides that possibly generated the observed spectrum (*e.g.* quadratic for two peptides), making it difficult to generalize to cases with more than two peptides. Thus, it is desirable to find a way to efficiently approximate this probability. To derive this approximation we borrow an intuition from the definition of conditional probability where the joint probability of two random events (R_{1,} R_{2}), is equal to the probability of one event times the conditional probability of the second event given the first event:

Analogously we can decompose the joint probability question into two simpler questions: (1) what is the probability P_{r}(*Score*(*R*_{1},*M*) ≥*T _{P}*) of finding a random peptide

*R*

_{1}that matches to

*M*with a score equal or better than T

_{p}=

*Score*(

*P*,

*M*)? and (2) once we find a first peptide

*P*, what is the probability Pr(

*Score*(

*R*

_{1},

*R*

_{2},

*M*) ≥

*T*|

*R*

_{1}=

*P*) of finding a random peptide

*R*

_{1}that together with

*P*scores equal or higher than

*T*when matched to

*M*? Note that the first question is just the single-peptide probability and the second question is the conditional probability. Therefore, we can define the following approximation:

From here on, we refer to this approximation as the *Product probability*. This formulation is not exactly equivalent to the definition of joint probability because it fixes *R*_{1} = *P* in the conditional probability term (where *P* is the first peptide in the PPSM) and thus does not explicitly consider the dependences between all possible *pairs* of peptides that can be matched to the mixture spectrum. However, both single-peptide probability and conditional probability can be computed efficiently in linear time and we show in the next section that this approximation is sufficiently accurate for our main use of the joint probability – to separate correct from incorrect matches to mixture spectra.

##### Classification of Matches

Because a typical proteomics data set contains both single-peptide and mixture spectra we consider three possible outcomes when searching a given query spectrum *M*: (1) *No-match: M* does not match any peptide in the database; (2) *Single-peptide match: M* matches one peptide in the database; and (3) *Mixture match: M* matches a pair of peptides in the database. Every query spectrum is initially assumed to be a putative mixture spectrum and is assigned to its top-scoring PPSM. Then a two-step procedure is used to separate true mixture matches from false mixture matches. At the first stage, all PPSMs with joint probability less than a threshold are accepted. Then PPSMs with conditional probability less than a second threshold are accepted as Mixture-matches. The probability thresholds are determined in a way such that it enforces a selected false discovery rate (FDR, see next section). Next, all the remaining spectra that do not pass either probability threshold are reconsidered as single-peptide spectra. Each PPSM is converted into a PSM by considering the first peptide as the match to the spectrum. Single-peptide probabilities are computed for all PSMs and a probability threshold is determined to enforce a selected FDR for single-peptide spectra. A graphical illustration of this classification procedure is provided in the Fig. 1.

##### Estimation of False Discovery Rates

In the classification steps of PPSMs, the probability thresholds are determined to enforce a certain FDR. For the joint probability we are interested in the FDR that an incorrect mixture match is accepted either as half-correct or correct match:

For the conditional probability we only want to accept correct matches so we are interested in the FDR where either half-correct or incorrect matches are accepted as correct mixture matches:

Each of the above FDRs is estimated by extending the Target-Decoy Approach (TDA) for single-peptide spectra (37). However, the assumptions used in TDA first need to be generalized to the case of PPSMs and their validity also needs to be tested (the detailed derivation of the TDA approach for PPSMs is given in the Supplementary Material). In brief, for a set of PPSMs if we define *TT* to be the number of PPSMs where both peptide matches are from the target database; *TD* or *DT* to be the number of cases where one peptide is from the target and the other peptide is from the decoy database and *DD* to be the cases where both peptides are from the decoy database, the two FDRs mentioned above can be computed using the following formulae:
and

Finally, for the single-peptide probability the FDR is estimated using the standard TDA approach:
where *T* is the number of PSMs from the target database and *D* is the number PSMs from the decoy database.

In summary two types of matches can be returned by MixGF: Single-match and Mixture-match. *FDR _{Single}* enforces the FDR for Single-match, whereas

*FDR*and

_{Joint}*FDR*enforce the FDR for Mixture-match. All three

_{Conditional}*FDR*operates on PSMs level (a Mixture-match is essentially treated as two PSMs, see Supplementary Material). Therefore to enforce a global FDR of 1% for all matches returned by MixGF, all three FDR thresholds were set to 1%.

##### Data sets and Data Processing

The performance of MixGF was first evaluated on a set of simulated mixture spectra (21). In brief, mixture spectra were created by linearly combining two single-peptide spectra with predefined mixture coefficients α – a parameter that reflects the relative abundance of the two peptides in the mixture spectrum. In addition, MixGF was tested on three data sets (38, 39, 40) representing typical experimental setups in proteomics studies. In brief, the *Yeast data set* (38) is from a tryptic digest of *Saccharomyces cerevisiae* that was analyzed on an LTQ Orbitrap XL mass spectrometer (Thermo Fisher Scientific) and MS/MS spectra were acquired using a data-dependent scanning mode in which each full MS scan (*m*/*z* 300–2000) was acquired on the Orbitrap at resolution 60,000, followed by eight MS/MS scans collected on the LTQ. Two Human data sets were also analyzed. The *Human-L* data set (39) (which stands for human data set with low mass accuracy MS/MS) is from a tryptic digest of HEK293 cell lysate that was fractioned using Strong Cation Exchange (SCX) and each fraction was analyzed by an LTQ Orbitrap XL ETD mass spectrometer (Thermo Fisher Scientific) in data dependent mode. MS full scans were acquired from *m*/*z* 350–1500 with a resolution of 60,000. The two most intense ions were fragmented in the linear ion trap using CID and ETD. The *Human-H* data set (40) (human data set with high mass accuracy) is from a tryptic digest of Human LOVO cell lysates. The sample was analyzed using an LTQ-Orbitrap XL mass spectrometer operated in data-dependent mode with MS full scans acquired at resolution 30,000 and MS/MS (CID and ETD) for the three most abundant peptides recorded at resolution 7500. For the Human data sets, only CID spectra were considered in the present study. All database searches were performed with precursor m/z tolerance of 3.0 Th to allow for the possible identification of co-eluting peptides in mixture spectra (the instrument does not record any parent mass for the least abundant peptide). Only tryptic peptides were considered and only precursor charge of two and three are considered; no variable modifications and only carbamidomethylation on cystein (C+57.021) was considered as a fixed modification. The MS/MS fragment mass tolerance was set to 0.5 Da for the low mass accuracy data and 0.05 Da for the high mass accuracy human data. For ProbIDtree, we separated all the spectra into two sets depending on whether ProbIDtree identified only one or multiple peptides match to the spectrum, then an FDR was enforced using the standard TDA method (37) for each subset. As shown before (8) the rationale for separate FDR determination is that an FDR calculation combining both mixture and single-peptide spectra leads to underestimation of FDR for mixture spectra. M-SPLIT searches were done with 3.0Da parent mass tolerance against the yeast (ver. 05/04/2009) and human (ver. 01/14/2010) spectral libraries downloaded from National Institute of Standards and Technology (NIST) (41). The yeast spectral library contain 86,861 unique peptide ions (*i.e.* same peptides with different charge states are counted as different ions), whereas the human spectral library contains 343,301 unique peptide ions. The protein sequence databases used were the SGD yeast protein database (*ver.5/8/2009*) and the Human protein database (downloaded from NCBI refseq, *ver.10/29/2010*).

## RESULTS

##### Separating True and False Mixture Spectrum Matches

As described above, our main goal of computing the statistical significance of PPSMs is to separate correct mixture matches from half-correct and incorrect matches. To test MixGF's ability in these tasks, we built a set of simulated mixture spectra by linearly combining pairs of single-peptide spectra same as before (21). Because we know *a priori* the peptides that generated each simulated mixture spectrum, we can extract the top-scoring correct, half-correct, and incorrect matches returned by MixDB and compute their joint and conditional probabilities. As shown in Fig. 2A, joint probability performs very well when separating correct matches from incorrect matches but there is considerable overlap between the joint probability of correct-matches and that of half-correct matches (see Fig. 2*B*). Further investigation of cases in the overlap region shows that for correct-matches usually both peptides contribute moderate scores to the final combined score but for the half-correct matches the correct peptide often contributes a very high score and thus even when paired with an incorrect match, the resulting combined high score still yields a low joint probability. Intuitively in order to separate half-correct matches from correct matches we need to look for cases that have high combined score as well as both peptides contributing significantly to the total score. The concept of conditional probability defined above aims to address exactly this question—is the score of the peptide pair (*P*, *Q*) significantly higher than that of the single peptide *P*? As illustrated in Fig. 2*C*, conditional probability is indeed better at separating correct matches from half-correct matches. Therefore, a two-step procedure is used to separate correct matches from false matches: at the first stage of MixGF, joint probability is used to filter out incorrect matches and then conditional probability is used to filter out half-correct matches.

##### Approximating Joint by the Product of Conditional Probability

To test whether the approximation of joint probability is accurate, The joint and product probability were computed for each spectrum in the simulated mixture spectra data set. As shown in Fig. 3, in most cases the joint probability is accurately approximated by product probability as most data points clustered tightly along the main diagonal. For correct matches, the product probability is sometimes lower than the true joint probability. This can be attributed to the fact that the approximation does not explicitly consider all pairs of peptides – because *P* is fixed, there are less opportunities for false positive matches to achieve high scores and thus the resulting spectral probability can be smaller in such cases. However, the range of probabilities where this underestimation occurs is well below the range where incorrect matches tend to occur. Therefore, for the purpose of separating correct matches from incorrect matches, using the approximation is nearly equivalent to computing the exact joint probability. As shown in Fig. 3*B*, correct matches and incorrect matches remain very well-separated using either the product or joint probability. In addition, the product probability can be computed much more efficiently than the joint probability. In practice the average run time for joint probability is 205.7 s per PPSM, whereas the average time it takes to compute the product probability is only 0.12 s on a Windows 7 machine with an Intel Xeon(R) E5430 CPU, resulting in a ≈1700 times speedup. This makes MixGF computation time in similar scale as searching for the top-scoring PPSMs by MixDB which takes 0.49 s per spectrum on the Yeast data set and 1.08 s on the Human data sets with a 3.0Th precursor mass tolerance.

##### Joint and Product Probability Improve the Detection of Mixture Spectra

For mixture spectra, we expect joint probability to perform better at separating correct matches from incorrect matches by explicitly considering two peptides. Intuitively we expect that single-peptide probabilities for correct peptide matches to mixture spectra to be higher (*i.e.* worse) than those for correct matches to single-peptide spectra. This is because the presence of a second peptide in mixture spectra will allow more peptides to match to the spectrum with high score. However, for false matches the single-peptide probability distribution remains comparable for both single-peptide and mixture spectra because they are random matches in either case. Therefore, the distribution of single-peptide probabilities between correct and incorrect matches should be less well-separated for mixture spectra than for single-peptide spectra. To show this we used the simulated mixture spectra where the first peptide is mixed with a second peptide at 100%, 50%, 30% of the first peptide's total intensity and then computed the single-peptide probability, joint probability and product probability for the correct matches as well as the top-scoring incorrect matches. The performance of each probability function in separating correct from incorrect matches is shown in Table I. As expected, when the second peptide is at relatively low abundance (*i.e.* 30%), the FDR-controlled performance of single-peptide probability is nearly identical to that of joint probability because the mixture spectra are more similar to single-peptide spectra than at higher second-peptide abundances. However, as we increase the relative abundance of the second peptide, joint probability performs considerably better at separating correct matches from incorrect matches. Thus, we expect that as mixture spectra with more peptides become more common in experiments, joint probability and its product probability approximation will substantially improve our ability to identify mixture spectra.

##### Identification of Mixture Spectra in Complex Biological Samples

To illustrate MixGF's ability to identify mixture spectra in realistic scenarios, we evaluated its performance on one yeast (38) and two human data sets (39, 40) that represent typical proteomics analysis of whole cell lysates. We compared the performance of MixGF with two current database search methods for identification of mixture spectra: MixDB and ProbIDtree. As shown in Table II, MixGF is able to outperform MixDB and ProbIDtree by identifying 30–76% and 160–390% more mixture spectra, respectively. To estimate their sensitivity, the various database search methods were also benchmarked against M-SPLIT, a spectral library search method able to identify mixture spectra. In general, spectral library search methods have shown to be more sensitive than database search methods (42, 21) because of their smaller search space and the known peptide fragmentation contained in the spectral library (43). Thus, the number of mixture spectra identified by M-SPLIT can be used as an estimate of the upper bound of mixture spectra that can be identified in each sample by database search methods. As shown in Table II, whereas MixDB and ProbIDtree only identify an average of 50 and 21% of the mixture spectra identified by M-SPLIT (respectively), MixGF is closing the gap between database search and spectral-library search methods by identifying up to 84% of the mixture spectra identified by library search (in the Human-H data set). Comparing search results across different data sets reflects their varying sample complexity. In the yeast data set, the number of mixture spectra identified by MixGF is only 4.71% of the single-peptide spectra. However, in the human data set the number of identified mixture spectra increased to 17.8–23.7% of the identified single-peptide spectra. This shows that as the complexity of the sample increases, more peptides with similar precursors will co-elute at same retention time and thus mixture spectra indeed constitute a significant fraction of all MS/MS spectra in the data. Such observations have been reported in studies showing that multiple precursors often coincide within the same MS/MS precursor isolation window for as many as 50% of all MS/MS spectra (5, 4). But while it remains unclear what fraction of these cofragmented peptides is identifiable, here we show that a large number of co-eluting peptides can be identified. It is worth noting that because mixture spectra contain two peptides per spectrum, they potentially contain more information in each spectrum. For example in the yeast data set even though the number of mixture spectra identified by MixGF is only 4.71% of the single-peptide spectra, the number of unique peptides identified in mixture spectra is 22.1% of the total number of peptides identified in single-peptide spectra. For the human data sets the number of unique peptides identified in mixture spectra is a strikingly 50.3–65.24% of the number of peptides identified in single-peptide spectra. Of course, while many of the peptides identified in mixture spectra were also identified in single-peptide spectra, the combined result still leads to a gain of 23.4–34.5% more unique peptides identified than single-peptide spectra. The comparison between the results in Human-L and Human-H data set also reveals that the high mass accuracy available in the Human-H data enables more mixture spectra to be identified. The total numbers of spectra identified in both human data sets are similar but the fraction of all identified spectra that are mixture spectra increased by 33.2% in the Human-H data set. Analyzing the data with a regular search engine, MSGFDB (39), shows that it was able to identify a similar number of spectra as MixGF, indicating that a large fraction of mixture spectra were also identified as single-peptide matches by MSGFDB. However, MixGF's ability to identify two peptides per mixture spectrum leads to an improvement of 13.3–15.4% in the total number of unique peptide identifications compared with a conventional analysis tool.

## DISCUSSION

It is increasingly being recognized that MS/MS spectra from more than one peptide are common and important in proteomics data (4, 5) but processing them will require the development of accurate computational tools to identify multiple peptides in each MS/MS spectrum. Two fundamental questions need to be addressed in the pursuit of this goal: (1) to separate correct multiple-peptide-spectrum matches (mPSMs) from false positive matches, and (2) to estimate the false positive or false discovery rate (FDR) in a set of mPSMs. Here we addressed these questions by computing the statistical significance of mPSMs for the special case of two peptides per spectrum. Given an MS/MS spectrum, a database search tool can always return a top-scoring peptide or multiple peptides matched to the query spectrum and by random chance, it is always possible for false peptide matches to also obtain high scores. This is especially true in the case of mixture spectra because the explosion in the search space dramatically increases the occurrence of high-scoring false matches. Thus, it is crucial to be able to compute the statistical significance of mPSMs. Here we show that for the two-peptide case, it is possible to rigorously compute the statistical significance using the generating function approach and show that the joint and conditional probabilities can be used to separate true PPSMs from false positive matches. We further show that the computationally-expensive joint probability can be accurately approximated in our range of interest using a product of probabilities that can be computed efficiently in linear time, resulting in a ≈1700 times speedup and thus potentially allowing MixGF to scale to an arbitrary number of peptides. In order to estimate the false discovery rate (FDR) for mixture spectra, we extended the traditional target-decoy approach (TDA) to perform the database search using a concatenated target-decoy sequence database allowing for the estimation half-correct matches where one peptide in the PPSM is correct and one peptide is incorrect. This is important because these matches constitute a large fraction of false positive matches in mixture spectra and, as shown in Fig. 2, these are more difficult to separate from correct matches than cases when both peptides are incorrect.

Benchmarking MixGF performance on three data sets showed that MixGF identified 30.1–390% more mixture spectra than MixDB and ProbIDtree, respectively. In addition, we found that the number of unique peptides identified in mixture spectra can be up to 65% of the number of peptides identified in single-peptide spectra with up to 34.5% of all unique peptides being identified only in mixture spectra. Though, many of the peptides in mixture spectra were also identified by MSGFDB as single-peptide matches, making MixGF overall gain in peptide identification to be 13.3–15.4% compared with MSGFDB. These observations highlight the importance of relaxing the *one-peptide-one-spectrum* assumption when designing the next generation of computational tools for identifying MS/MS spectra as instruments advances and higher complexity samples are being analyzed in less time in proteomic experiments. This also illustrates the potential of emerging data acquisition protocols (9, 10, 11, 12, 13, 14, 15, 44) where multiple peptides are intentionally cofragmented in each MS/MS spectrum. Finally, even though the proposed MixGF approach focuses on mixture spectra from two peptides, the same approach should be extensible to more than two peptides. Focusing on the case of pairs allowed for a comprehensive assessment of the various performance aspects of the proposed approach and the analysis of several data sets showed that properly interpreting mixture spectra from two peptides already greatly improves the number of spectra and peptides identified in current experimental setups. Therefore solving the problem for the case of two peptides provides a foundation toward addressing the more general scenario of mixture spectra from any number of peptides.

## Acknowledgments

We thank the National Institute of Science and Technology (NIST), ProteomeCommons, and Vanderbilt University for the public availability of the mass spectrometry data and Dr. Fedor Kryuchkov and Frank Kjeldsen for the availability of the Human-H data.

## Footnotes

Author contributions: J.W., P.E.B., and N.B. designed research; J.W. performed research; J.W. and N.B. analyzed data; J.W., P.E.B., and N.B. wrote the paper.

↵* This work was supported in part by the National Institutes of Health grant GM078596 (JW and PEB) and 8 P41 GM103485-05 (JW and NB) from the National Center for Research Resources.

↵

^{}This article contains supplemental Fig. S1 and Table S1.↵

^{1}The abbreviations used are:- MS/MS
- tandem mass spectrometry
- M-SPLIT
- mixture-spectrum partitioning using libraries of identified tandem mass spectra
- PSM
- peptide spectrum match
- PPSM
- peptide/peptide spectrum match
- FDR
- false discovery rate
- PRM
- prefix residue mass
- mPSM
- multipeptide spectrum match
- DDA
- data dependent acquisition
- DIA
- data independent acquisition
- TDA
- target decoy approach.

- Received January 9, 2014.
- Revision received September 5, 2014.

- © 2014 by The American Society for Biochemistry and Molecular Biology, Inc.