PeakSelect: preprocessing tandem mass spectra for better peptide identification

We present a new preprocessing method, PeakSelect, to improve the accuracy and efficiency of Tandem Mass-Spec peptide (protein) identification. The fundamental difference between noise and fragment ions in spectra is that ions have isotopes but noise does not. We propose a new and important concept of an Isotope Pattern Vector (IPV) which characterizes the isotope cluster of fragment ions. Then the noise and real peaks can be distinguished by the quantitative IPV values. PeakSelect first uses a new method of the Gaussian Mixture Model and Expectation-Maximization (EM) algorithm to find the base intensity level (baseline) in a spectrum. Then PeakSelect selects features based on the IPV and baseline, and constructs a decision tree to automatically classify the peaks into different categories such as noise, single ion peaks, and overlapping peaks. Experiments show that PeakSelect can help to reduce the Mascot searching time and increase the reliability of peptide identifications. In particular, PeakSelect performs well on complex spectra with a large number of peaks from large peptides, and supports more sequence identification than other well-known systems.     

Competing fragmentation processes in tandem mass spectra of heparin-like glycosaminoglycans

Heparin-like glycosaminoglycans (HLGAGs) are highly sulfated, linear carbohydrates attached to proteoglycan core proteins and expressed on cell surfaces and in basement membranes. These carbohydrates bind several families of growth factors and growth factor receptors and act as coreceptors for these molecules. Tandem mass spectrometry has the potential to increase our understanding of the biological significance of HLGAG expression by providing a facile means for sequencing these molecules without the need for time-consuming total purification. The challenge for tandem mass spectrometric analysis of HLGAGs is to produce abundant ions derived via glycosidic bond cleavages while minimizing the abundances of ions produced from elimination of the fragile sulfate groups. This work describes the competing fragmentation pathways that result from dissociation of high negative charge state ions generated from HLGAGs. Glycosidic bond cleavage ion formation competes with losses of equivalents of H2SO4, resulting in complex ion patterns. For the most highly sulfated structure examined, an octasulfated tetramer, an unusual loss of charge from the precursor ion was observed, accompanied by low abundance ions originating from subsequent backbone cleavages. These results demonstrate that fragmentation processes competing with glycosidic bond cleavages are more favored for highly sulfated HLGAG ions. In conclusion, reduction of charge-charge repulsions, such as is achieved by pairing the HLGAG ions with metal cations, is necessary in order to minimize the abundances of ions derived via fragmentation processes that compete with glycosidic bond cleavages.     

SILVER helps assign peptides to tandem mass spectra using intensity-based scoring

Tandem mass spectrometry is commonly used to identify peptides (and thereby proteins) that are present in complex mixtures. Peptide identification from tandem mass spectra is partially automated, but still requires human curation to resolve "borderline" peptide-spectrum matches (PSMs). SILVER is web-based software that assists manual curation of tandem mass spectra, using a recently developed intensity-based machine-learning approach to scoring PSMs, Elias et al. In this method, a large training set of peptide, fragment, and peak-intensity properties for both matched and mismatched PSMs was used to develop a score measuring consistency between each predicted fragment ion of a candidate peptide and its corresponding observed spectral peak intensity. The SILVER interface provides a visual representation of match quality between each candidate fragment ion and the observed spectrum, thereby expediting manual curation of tandem mass spectra. SILVER is available online at http://llama.med.harvard.edu/Software.html.     

Fast-atom bombardment-induced condensation of glycerol with ammonium surfactants II: Time dependence of mass spectra and tandem mass spectra

Studies of fast-atom bombardment (FAB)-induced condensation between trimethyltetradecylammonium cations and glycerol have been extended to consider spectral time dependence. To enhance reproducibility of time dependence, a modified FAB target was used. FAB mass spectrometry of deuterium-labeled surfactants and FAB/collision-induced dissociation (CID) of nonlabeled material demonstrate that products of condensation at the surfactant “head group” predominate early in the analysis, while tail adducts become prominent later. This time dependence correlates with the expected surface activity of the products. It is incompatible with gas-phase reaction, but consistent with reaction in the condensed phase. Subtle variations in the surface activities of various condensation products (derived from changes in the number of hydroxyls from the reactive glycerol radical or in the position of attack along the surfactant chain) are reflected in the time dependence of FAB and CID spectra. CID spectra of deuterium-labeled cations provide evidence for intramolecular hydrogen transfer from the surfactant tail to the head within a surfactant radical. This transfer shows no significant kinetic isotope effect.     

Statistical evaluation of electrospray tandem mass spectra for optimized peptide fragmentation

Liquid chromatography-tandem mass spectrometry (LC-MS/MS) has become the method of choice for the analysis of complex peptide mixtures. It combines the separation power of nanoflow LC with highly specific sequence analysis, allowing automated peptide sequencing with high resolution and throughput. For peptide fragmentation, the current experimental setup uses predefined parameters based on the mass-to-charge ratio of the individual precursor. Suitable parameters are typically established by empirical evaluation of fragment spectra of individual peptides used as standards. As a result, nonoptimal fragment spectra are obtained if peptides show fragmentation behavior different from these standards, which often result in the loss of sequence-specific fragment ion information. Here we describe a statistical approach for the systematic evaluation of the quality of individual peptide fragment spectra based on the calculation of their arithmetic mean and standard deviation. The method utilizes the dependence of these parameters on the difference in electric potential across the collision cell to determine the value that results in maximum information content. We show that the method is applicable to fragment spectra generated from a variety of multiply-charged tryptic peptides, over a wide concentration range, and on different types of mass analyzers. We also show how this novel approach can be used to define optimized collision energy settings over a wide mass-to-charge range.     

Automatic recalibration and processing of tandem mass spectra using formula annotation

High accuracy, high resolution tandem mass spectrometry (MS/MS) is becoming more common in analytical applications, yet databases of these spectra remain limited. Databases require good quality spectra with sufficient compound information, but processing, calibration, noise reduction and retrieval of compound information are time‐consuming tasks that prevent many contributions. We present a comprehensive workflow for the automatic processing of MS/MS using formula annotation for recalibration and cleanup to generate high quality spectra of standard compounds for upload to MassBank ( www.massbank.jp ). Compound information is retrieved via Internet services. Reference standards of 70 pesticides were measured at various collision energies on an LTQ‐Orbitrap XL to develop and evaluate the workflow. A total of 944 resulting spectra are now available on MassBank. Evidence of nitrogen adduct formation during MS/MS fragmentation processes was found, highlighting the benefits high accuracy MS/MS offers for spectral interpretation. A database of recalibrated, cleaned‐up spectra resulted in the most correct spectra ranked in first place, regardless of whether the search spectra were recalibrated or not, whereas the average rank of the correct molecular formula was improved from 2.55 (uncalibrated) to 1.53 when using recalibrated MS/MS data. The workflow is available as an R package RMassBank capable of generating MassBank records from raw MS and MS/MS data and can be adjusted to process data acquired with different settings and instruments. This workflow is a vital step towards addressing the need for more high quality, high accuracy MS/MS spectra in spectral databases and provides important information for spectral interpretation. Copyright © 2012 John Wiley & Sons, Ltd.     

Generation of substructure identification rules using feature-combinations from tandem mass spectra

Software to interpret tandem mass spectra, entitled Method for Analyzing Patterns in Spectra (MAPS), has been developed to provide substructure information for an automated compound identification system, This software consists of several program modules which manipulate databases of tandem mass spectra and substructure information, generate substructure identification rules, and apply these rules to the tandem mass spectra of unknown compounds to identify components of their structure. The MAPS rule generation program has been modified to generate rules based on specific combinations of spectral features that occur concertedly. False positives are drastically reduced by searching for “feature-combinations” that have 100% uniqueness with respect to a reference database of compounds. Recall is increased by the determination of multiple feature-combinations indicative of the presence of a given substructure. Strategies were developed in the algorithm for the discovery of feature-combinations that avoid the computation “explosion” that occurs when working with a large number of spectral features. The rules developed have the form: “IF feature-eombination a (FC a) or FC b,..., or FC x, THEN substructure SSn is present.”     

Improved peptide sequencing using isotope information inherent in tandem mass spectra

We demonstrate here the use of natural isotopic 'labels' in peptides to aid in the identification of peptides with a de novo algorithm. Using data from ion trap tandem mass spectrometric (MS/MS) analysis of 102 tryptic peptides, we have analyzed multiple series of peaks within LCQ MS/MS spectra that 'spell' peptide sequences. Isotopic peaks from naturally abundant isotopes are particularly prominent even after peak centroiding on y- and b-series ions and lead to increased confidence in the identification of the precursor peptides. Sequence analysis of the MS/MS data is accomplished by finding sequences and subsequences in a hierarchical manner within the spectra.     

Evaluation of automatically generated substructure identification rules from tandem mass spectra

Substructure identification rules for phenothiazine and barbiturate substructures were generated by using a new version of the Method for Analyzing Patterns in Spectra (MAPS) software. This software uses tandem mass spectra and known substructure content of reference compounds to provide “feature-combination“ rules. A feature-combination is a series of tandem mass spectral features which are completely unique to compounds containing a specified substructure. The current reference databases contain over 11,000 daughter spectra of 100 compounds acquired at two different collision gas pressures (i.e., single- and multiple-collision conditions). The results of rule evaluation procedures are presented and include a comparison of the spectral features developed in rule generation to those identified in documented fragmentation pathways of the indicated substructure. Two potential sources of error due to spectral feature and substructure “cross-correlation“ were identified. If errors occur, they can be detected by calculating cross-correlation coefficients and edited from the rules. A beneficial cross-correlation involving feature-combinations was also discovered. The rules obtained by using single- and multiple-collision data were further evaluated by applying them to tandem mass spectra of 20 test compounds (compounds not in the reference database). The results of these evaluations give a good indication of the utility of the rules for use in an automated structure elucidation system for tandem mass spectrometry data.     

RScore: a peptide randomicity score for evaluating tandem mass spectra

RScore, a new criterion of randomicity for evaluating tandem mass (MS/MS) spectra, is described. RScore is defined as the relative quality in cross-correlation and matched intensity percentage of a potentially positive peptide to those of other possible candidates for the same spectrum. By utilizing RScore combined with less stringent SEQUEST score filters, the number of true positive peptides can be increased and the number of false positives in datasets from a known protein mixture can be reduced compared with current SEQUEST parameters used alone. This algorithm is simple and adds little overheads to SEQUEST computation.     

Automated structural assignment of derivatized complex N-linked oligosaccharides from tandem mass spectra

Glycoprotein function is controlled by several biological factors, one of them being the structure of carbohydrate chains (glycans) attached to specific amino acids of the protein backbone. Changes in glycan structures have been shown to modify the secondary and tertiary conformation of glycoproteins, thus their function. Powerful analytical tools are available for the characterization of sugar structures, and recently mass spectrometry (MS) has been increasingly useful for this purpose. Manual interpretation of tandem mass spectrum is possible but tedious. Automated interpretation should speed the analysis and enhance the results obtained. A new computer program for automated interpretation of tandem MS spectra of complex N-linked glycans oligosaccharides from mammals will be described. N-Linked oligosaccharides standards were derivatized with 1-phenyl-3-methyl-5-pyrazolone (PMP) and analyzed by matrix-assisted laser desorption/ionization (MALDI)-tandem MS. Simulated tandem mass spectra of other common glycans were also generated to test the algorithm. The MALDI-MS/MS spectra featured resolved isotopic distributions for the [M + H](+) and fragment ions of oligosaccharides. These isotopic distributions complicated the automated analysis of the spectra and were removed to leave only monoisotopic peaks. An algorithm was written for this purpose, yielding simplified tandem mass spectra. Another algorithm is then used to determine the structure of the oligosaccharide. A score is then given to each structure, depending on agreement with experimental results. The program successfully assigned the true structure in 24 out of the 28 cases (86%) and the true structure was among the three top scoring structures in all cases.     

Effect of electrospray ionization conditions on low-energy tandem mass spectra of peptides

The effect of electrospray ionization (ESI) conditions on low-energy tandem mass spectra of peptides in the relative molecular mass range 400–1200 was examined. For singly charged peptide ions the source skimmer potential (which determines the degree of acceleration of the ions through the intermediate pressure region in the source) can strongly influence the extent of fragmentation observed in tandem mass spectra, especially at low collision energies. For each peptide there is an optimum skimmer potential which represents a balance between generating ions with sufficient internal energy for subsequent tandem mass spectrometric experiments and inducing the onset of other processes such as source fragmentation. The fragmentation which can be achieved in tandem mass spectra with high skimmer potentials differs from ESI source fragmentation for the same peptides. We have found that fragmentation in ESI mass spectra depends both on skimmer potential and on solvent pH, presumably because the latter determines the proportion of doubly charged species generated from a given peptide. Low-energy tandem mass spectra of peptides following ESI are equally as sensitive to peptide structure and the type of adduct studied (e.g. [M + H]⁺ vs. [M + NH₄]⁺) as tandem mass spectra obtained following older ionization methods such as fast atom bombardment.     

Effect of metal cationization on the tandem mass spectra of glycosyl dithioacetals

The effect of metal cationization on the tandem mass spectra of glycosyl dithioacetals of glucose, mannose, galactose, rhamnose, arabinose and xylose was studied by electrospray ionization mass spectrometry under ammonium and metal (Li, Na, Ag and Cu) ion cationization conditions. The ammonium-cationized glycosyl dithioacetals fragment by loss of ammonia followed by either two molecules of EtSH or one molecule of EtSH and one molecule of H2O. Lithium cationization leads to additional eliminations such as EtSEt and EtSSEt and C-C cleavages. Elimination of EtSH is not observed under sodium cationization. Silver cationization, on the other hand, leads to additional fragmentations involving the elimination of silver as AgOH and AgSEt. Copper cationization results in adducts where copper has undergone a change of oxidation state from II to I. Li+, Ag+ and Cu+ cationization seem to favour cyclization resulting in elimination of EtSH. However, the mechanisms seem to be differently affected by different metal ions. Li+ and Ag+ cationization appear to be non-specific and favour cyclization involving C2-, C4- and C5-hydroxyl hydrogens, whereas Cu+ cationization seems to favour cyclization involving C4-hydroxyl hydrogen.     

Comparison of negative ion electrospray mass spectra measured by seven tandem mass analyzers towards library formation

A library of negative ion electrospray ionization mass spectra and tandem mass spectra (MS/MS) of sulfonated dyes has been developed for fast identification purposes. The uniform protocol has been elaborated and applied to the measurements of more than 50 anionic dyes. Three collision energies are selected in our protocol which ensures that at least one of them provides a suitable ratio of product ions to the precursor ion. The robustness is investigated with altered values of tuning parameters (e.g. the pressure of the nebulizing gas, the temperature and the flow rate of drying gas, and the mobile phase composition). The results of the inter-laboratory comparison of product ion mass spectra recorded on seven different tandem mass spectrometers (three ion traps, two triple quadrupoles and two hybrid quadrupole time of flight instruments) are presented for four representative anionic dyes--azo dye Acid Red 118, anthraquinone dye Acid Violet 43, triphenylmethane dye Acid Blue 1 and Al(III) metal-complex azo dye. The fragmentation patterns are almost identical for all tandem mass analyzers, only the ratios of product ions differ somewhat which confirms the possibility of spectra transfer among different mass analyzers with the goal of library formation.     

Pattern-based algorithm for peptide sequencing from tandem high energy collision-induced dissociation mass spectra

A new strategy is reported for extracting complete and partial sequence information from collision-induced dissociation (CID) spectra of peptides, CID spectra are obtained from high energy CID of peptide molecular ions on a four-sector tandem mass spectrometer with an electro-optically coupled microchannel array detector, A peak detection routine reduces the spectrum to a list of peak masses and peak heights, which is then used for sequencing, The sequencing algorithm was designed to use spectral data to generate sequence fits directly rather than to use data to test the fit of series of sequence guesses. The peptide sequencing algorithm uses a pattern based on the polymeric nature of peptides to classify spectral peaks into sets that are related in a sequence-independent manner, It then establishes sequence relationships among these sets, Peak detection from raw data takes 10–20 s, with sequence generation requiring an additional 10–60 s on a Sun 3/60 workstation, The program is written in the C language to run on a Unix platform. The principal advantages of our method are in the speed of analysis and the potential for identifying modified or rare amino acids. The algorithm was designed to permit real-time sequencing but awaits hardware modifications to allow real-time access to CID spectra.     

Effect of tautomerization on the fast-atom bombardment tandem mass spectra of azo dyes

The effect of tautomerization on the collision-induced dissociation of negative ions produced from sulfonated azo dyes by fast-atom bombardment (FAB) was studied by analyzing the product ion spectra of several related compounds. The mechanisms by which azo dyes fragment were found to depend on the formation and stability of tautomers. The extent of tautomerization was affected by the number and location of hydroxy groups on the dye, as well as by the FAB matrix. Ions that retained a sodium were often inhibited from forming tautomers and gave different product ions. Substitution of deuterium for hydrogen on the hydroxy groups aided in the identification of ions having more than one possible structure and provided verification of proposed mechanisms. Mechanisms involving ions retaining a sodium were verified by substituting potassium for the sodium.