Clustering and Filtering Tandem Mass Spectra Acquired in Data-Independent Mode

Data-independent mass spectrometry activates all ion species isolated within a given mass-to-charge window (m/z) regardless of their abundance. This acquisition strategy overcomes the traditional data-dependent ion selection boosting data reproducibility and sensitivity. However, several tandem mass (MS/MS) spectra of the same precursor ion are acquired during chromatographic elution resulting in large data redundancy. Also, the significant number of chimeric spectra and the absence of accurate precursor ion masses hamper peptide identification. Here, we describe an algorithm to preprocess data-independent MS/MS spectra by filtering out noise peaks and clustering the spectra according to both the chromatographic elution profiles and the spectral similarity. In addition, we developed an approach to estimate the m/z value of precursor ions from clustered MS/MS spectra in order to improve database search performance. Data acquired using a small 3 m/z units precursor mass window and multiple injections to cover a m/z range of 400–1400 was processed with our algorithm. It showed an improvement in the number of both peptide and protein identifications by 8 % while reducing the number of submitted spectra by 18 % and the number of peaks by 55 %. We conclude that our clustering method is a valid approach for data analysis of these data-independent fragmentation spectra. The software including the source code is available for the scientific community.

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.     

Electron Transfer Dissociation of Milk Oligosaccharides

For structural identification of glycans, the classic collision-induced dissociation (CID) spectra are dominated by product ions that derived from glycosidic cleavages, which provide only sequence information. The peaks from cross-ring fragmentation are often absent or have very low abundances in such spectra. Electron transfer dissociation (ETD) is being applied to structural identification of carbohydrates for the first time, and results in some new and detailed information for glycan structural studies. A series of linear milk sugars was analyzed by a variety of fragmentation techniques such as MS/MS by CID and ETD, and MS(3) by sequential CID/CID, CID/ETD, and ETD/CID. In CID spectra, the detected peaks were mainly generated via glycosidic cleavages. By comparison, ETD generated various types of abundant cross-ring cleavage ions. These complementary cross-ring cleavages clarified the different linkage types and branching patterns of the representative milk sugar samples. The utilization of different MS(3) techniques made it possible to verify initial assignments and to detect the presence of multiple components in isobaric peaks. Fragment ion structures and pathways could be proposed to facilitate the interpretation of carbohydrate ETD spectra, and the main mechanisms were investigated. ETD should contribute substantially to confident structural analysis of a wide variety of oligosaccharides.     

Screening non-colored phenolics in red wines using liquid chromatography/ultraviolet and mass spectrometry/mass spectrometry libraries

Liquid chromatography/ultraviolet (LC/UV) and mass spectrometry/mass spectrometry (MS/MS) libraries containing 39 phenolic compounds were established by coupling a LC and an ion trap MS with an electrospray ionization (ESI) source, operated in negative ion mode. As a result, the deprotonated [M-H]- molecule was observed for all the analyzed compounds. Using MS/MS hydroxybenzoic acid and hydroxycinnamic acids showed a loss of CO2 and production of a [M-H-44]- fragment and as expected, the UV spectra of these two compounds were affected by their chemical structures. For flavonol and flavonol glycosides, the spectra of their glycosides and aglycones produced deprotonated [M-H]- and [A-H]- species, respectively, and their UV spectra each presented two major absorption peaks. The UV spectra and MS/MS data of flavan-3-ols and stilbenes were also investigated. Using the optimized LC/MS/MS analytical conditions, the phenolic extracts from six representative wine samples were analyzed and 31 phenolic compounds were detected, 26 of which were identified by searching the LC/UV and MS/MS libraries. Finally, the presence of phenolic compounds was confirmed in different wine samples using the LC/UV and LC/MS/MS libraries.     

SpectraBank: an open access tool for rapid microbial identification by MALDI-TOF MS fingerprinting

MALDI-TOF MS has proved to be an accurate, rapid, and cost-effective technique for microbial identification in which the spectral fingerprint of an unknown strain can be compared to a database of spectra from reference strains. Most of the existing databases are private and often costly to access, and little spectral information is shared among researchers. The objective of the present communication is to introduce the SpectraBank database (http://www.spectrabank.org), which provides open access MALDI-TOF mass spectra from a variety of microorganisms. This work aims to familiarize readers with the SpectraBank database, from the sample preparation, data collection, and data analysis to how the spectral reference data can be used for microbial species identification. The database currently includes more than 200 MALDI-TOF MS spectra from more than 70 bacterial species and links to the freely available web-based application SPECLUST (http://bioinfo.thep.lu.se/speclust.html) to allow comparisons of the obtained peak mass lists and evaluate phyloproteomic relationships. The SpectraBank database is intended to be expanded by the addition of new spectra from microbial strains, obtained in our laboratory and by other researchers.     

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.     

Semi-automated alignment and quantification of peaks using parallel factor analysis for comprehensive two-dimensional liquid chromatography–diode array detector data sets

Parallel factor analysis was used to quantify the relative concentrations of peaks within four-way comprehensive two dimensional liquid chromatography–diode array detector data sets. Since parallel factor analysis requires that the retention times of peaks between each injection are reproducible, a semi-automated alignment method was developed that utilizes the spectra of the compounds to independently align the peaks without the need for a reference injection. Peak alignment is achieved by shifting the optimized chromatographic component profiles from a three-way parallel factor analysis model applied to each injection. To ensure accurate shifting, components are matched up based on their spectral signature and the position of the peak in both chromatographic dimensions. The degree of shift, for each peak, is determined by calculating the distance between the median data point of the respective dimension (in either the second or first chromatographic dimension) and the maximum data point of the peak furthest from the median. All peaks that were matched to this peak are then aligned to this common retention data point. Target analyte recoveries for four simulated data sets were within 2% of 100% recovery in all cases. Two different experimental data sets were also evaluated. Precision of quantification of two spectrally similar and partially coeluting peaks present in urine was as good as or better than 4%. Good results were also obtained for a challenging analysis of phenytoin in waste water effluent, where the results of the semi-automated alignment method agreed with the reference LC–LC MS/MS method within the precision of the methods.     

ESIprot: a universal tool for charge state determination and molecular weight calculation of proteins from electrospray ionization mass spectrometry data

Electrospray ionization (ESI) ion trap mass spectrometers with relatively low resolution are frequently used for the analysis of natural products and peptides. Although ESI spectra of multiply charged protein molecules also can be measured on this type of devices, only average spectra are produced for the majority of naturally occurring proteins. Evaluating such ESI protein spectra would provide valuable information about the native state of investigated proteins. However, no suitable and freely available software could be found which allows the charge state determination and molecular weight calculation of single proteins from average ESI‐MS data. Therefore, an algorithm based on standard deviation optimization (scatter minimization) was implemented for the analysis of protein ESI‐MS data. The resulting software ESIprot was tested with ESI‐MS data of six intact reference proteins between 12.4 and 66.7 kDa. In all cases, the correct charge states could be determined. The obtained absolute mass errors were in a range between ?0.2 and 1.2 Da, the relative errors below 30 ppm. The possible mass accuracy allows for valid conclusions about the actual condition of proteins. Moreover, the ESIprot algorithm demonstrates an extraordinary robustness and allows spectral interpretation from as little as two peaks, given sufficient quality of the provided m/z data, without the necessity for peak intensity data. ESIprot is independent from the raw data format and the computer platform, making it a versatile tool for mass spectrometrists. The program code was released under the open‐source GPLv3 license to support future developments of mass spectrometry software. Copyright © 2010 John Wiley & Sons, Ltd.     

i-RUBY: a novel software for quantitative analysis of highly accurate shotgun-proteomics liquid chromatography/tandem mass spectrometry data obtained without stable-isotope labeling of proteins

We developed a novel software named i-RUBY (identification-Related qUantification-Based strategY algorithm for liquid chromatography/tandem mass spectrometry (LC/MS/MS) data) that enables us to perform fully automatic ion current-based spectral feature analysis of highly accurate data obtained by LC/MS/MS. At the 1st step, this software utilizes accurate peptide/protein identification information for peak detection and peak matching among measurements. Then, at the 2nd step, it picks yet unidentified peaks and matches them to the peaks identified at the 1st step by a linear interpolation algorithm. The analysis of human plasma externally spiked with a known amount of yeast alcohol dehydrogenase 1 showed a good linear relationship between the amount of protein spiked and the quantitative values obtained by i-RUBY analysis. Experiment using human plasma digests spiked with a mixture of known amounts of synthetic peptides derived from two yeast proteins, alcohol dehydrogenase 1 and glucose-6-phospate isomerase, showed the expansion by the 2nd step of i-RUBY of the lower quantification limits to 1/10 to 1/1000 of those reached only by identified peaks at the 1st step. Good correlations between the i-RUBY results and the amount of proteins were confirmed by the analysis of real samples, i.e., sera of normal subjects and cancer patients, by comparing quantitative values of acute-phase proteins obtained by i-RUBY analysis of LC/MS/MS data with those obtained by an immunological method using Bio-Plex. These results taken together show that i-RUBY is a useful tool for obtaining dependable quantitative information from highly accurate shotgun-proteomics LC/MS/MS data.     

Weighted two-band target entropy minimization for the reconstruction of pure component mass spectra: Simulation studies and the application to real systems

A method is proposed, on the basis of a recently developed algorithm--Band Target Entropy Minimization (BTEM)--to reconstruct mass spectra of pure components from mixture spectra. This method is particular useful in dealing with spectral data with discrete features (like mass spectra). Compared to the original BTEM, which has been applied to differentiable spectroscopies such as Fourier-transfer infrared spectroscopy (FTIR), ultraviolet (UV), Raman, and nuclear magnetic resonance (NMR), the latest modifications were obtained through: (1) Reformulating the objective function using the peak heights instead of their derivatives; (2) weighting the abstract vector VT to reduce the effect of noise; (3) using a two-peak targeting strategy (tBTEM) to deal with strongly overlapping peaks; and (4) using exhaustive search to locate all the component spectra. A set of 50 multi-component mass spectra was generated from ten reference experimental pure component spectra. Many of the compounds chosen have common MS fragments and therefore, many of the pure component spectra have considerable intensity in same data channels. In addition, a set of MS spectra from a real system with four components was used to examine the newly developed algorithm. Successful reconstruction of the ten component spectra of the simulated system and the four component spectra of the real system was rapidly achieved using the new tBTEM algorithm. The advantages of the new algorithm and its implication for rapid system identification of unknown mixtures are readily apparent.     

Mass spectral studies on aryl-substituted N-carbamoyl/N-thiocarbamoyl narcotine and related compounds

Positive ion mass spectral fragmentation of new N-carbamoyl/N-thiocarbamoyl derivatives of narcotine and compounds closely related to it are reported and discussed. The techniques used include electron impact (EI), fast-atom bombardment (FAB), matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) and electrospray ionization tandem mass spectrometry (ESI-MS/MS). Prominent peaks in the mass spectra of these compounds appear to involve C-C bond cleavage beta to the amine nitrogen with loss of the 4,5-dimethoxy(1H)isobenzofuranone moiety from their molecular ions, along with another prominent peak at m/z 382. No molecular ion peaks of these compounds were recorded in EI, whereas intense [M + H]+ ion peaks were observed in FAB and ESI spectra. MALDI also yielded [M + H]+ ion peaks in good agreement with FAB and ESI studies.     

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.     

Characterisation of the ester-substituted products of the reaction of p-t-butyl calix[4]arene and ethyl bromoacetate using LC-UV-MS and LC-DAD

A series of derivatisation reactions between p-t-butyl calix[4]arene and ethyl bromoacetate were carried out in order to prepare 1,3 diester substituted calix[4]arene. Mass spectral data, obtained from direct injection of samples, indicated that the reactions were rich in the desired product. Since the ultra violet (UV) spectra of the desired product and possible impurities are very similar, liquid chromatography (LC) chromatographic data seemed to corroborate these results. However, when on-line LC-UV-MS was carried out and each LC peak subjected to MS analysis as it eluted, a very different picture emerged. It was found that many of these reactions actually contained high levels of the monoester product which, having less affinity for sodium in the MS, is therefore seriously underestimated in any direct injection assay. LC-diode array detection (DAD) methods were also used to help successfully identify and characterise the compounds being formed in these complex reactions. The overall results obtained in this paper allowed the optimal reaction conditions to be determined for this reaction. LC-MS analysis of the chromatographic peaks also identified the presence of two isomers of the diester substituted calix[4]arene (1,3 and 1,2 diesters). The combination of LC and UV/MS detection is required for accurate analysis of the products of such reactions.     

Prediction of artifact peak intensity in linked scans for dissociations occurring in the first field-free region of sector mass spectrometers

Linked scans are commonly used on double-focusing mass spectrometers to obtain tandem mass spectrometry (MS/MS) spectra. The appearance of artifact peaks in linked scan MS/MS spectra from dissociations occurring in the first field-free region are a result of poor parent ion resolution, and they often can complicate the interpretation of the MS/MS spectra. The kinetic energy release associated with dissociation of ions of similar m/z to the “selected” parent ion is the main factor in determining the intensity of artifact peaks. A means of predicting the intensities of these artifact peaks in product ion and constant neutral loss scans is presented here. The method requires straightforward calculations based on Lacey-Macdonaldion intensity diagrams. The exact calculations require knowledge of the kinetic energy release of a particular dissociation, the kinetic energy spread of the main beam, and the parent ion and product ion mass-to-charge ratios. Adequate predictions, however, can be made by assuming a general kinetic energy release for any given reaction and a typical instrument energy resolution. Theoretical predictions are in good agreement with experimental data obtained from the product ion scans of unlabeled and isotopically labeled tirilazad and unlabeled and labeled leucine enkephalin methyl ester. There is also excellent agreement between experiment and theory in the constant neutral loss scans of rubidium bromide clusters.     

Artifacts in Fourier transform mass spectrometry

Recent work on a new, higher sensitivity preamplifier design for Fourier transform ion cyclotron resonance mass spectrometry (FTICRMS) revealed a number of artifact peaks (spectral features) which do not contain useful chemical information. In order to determine the cause of these artifacts and eliminate them, these severely distorted spectra were compared with similarly distorted signal models. The source of several common signal processing artifacts was thereby determined and correlated to radio‐frequency interference (RFI) noise and saturation of the amplifier and/or the digitizer. Under such conditions, the fast Fourier transform (FFT) generates spectral artifact peaks corresponding to harmonics and mixing frequencies of the real signal peaks and RFI frequencies. While this study was done using FTICRMS data, it is important to stress that these artifacts are inherent to the digitization and FFT process and thus are relevant to any FT‐based MS instrument, including the orbitrap and FT ion trap. Copyright © 2009 John Wiley & Sons, Ltd.     

Discrimination of Penicillium isolates by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry fingerprinting

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) was used to generate highly reproducible mass spectral 'fingerprints' for twelve Penicillium species. Prior to MALDI-TOF MS analysis, eight replicate cultures of each Penicillium species were subjected to three one-minute bead-beating cycles in an acetonitrile/trifluoroacetic acid solvent. The mass spectra contained abundant peaks in the range of m/z 5000-20 000, and allowed unambiguous discrimination between species. In addition, a biomarker common to all Penicillium mass spectra was observed at m/z 13 900. Discriminant analysis using the MALDI-TOF MS data yielded classification error rates of 0% (i.e. 100% correct identification), indicating that MALDI-TOF MS data may be a useful diagnostic tool for the objective identification of Penicillium species of environmental and clinical importance.     

Comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry peak sorting algorithm

We report a novel peak sorting method for the two-dimensional gas chromatography/time-of-flight mass spectrometry (GC x GC/TOF-MS) system. The objective of peak sorting is to recognize peaks from the same metabolite occurring in different samples from thousands of peaks detected in the analytical procedure. The developed algorithm is based on the fact that the chromatographic peaks for a given analyte have similar retention times in all of the chromatograms. Raw instrument data are first processed by ChromaTOF (Leco) software to provide the peak tables. Our algorithm achieves peak sorting by utilizing the first- and second-dimension retention times in the peak tables and the mass spectra generated during the process of electron impact ionization. The algorithm searches the peak tables for the peaks generated by the same type of metabolite using several search criteria. Our software also includes options to eliminate non-target peaks from the sorting results, e.g., peaks of contaminants. The developed software package has been tested using a mixture of standard metabolites and another mixture of standard metabolites spiked into human serum. Manual validation demonstrates high accuracy of peak sorting with this algorithm.     

The nature of collision-induced dissociation processes of doubly protonated peptides: comparative study for the future use of matrix-assisted laser desorption/ionization on a hybrid quadrupole time-of-flight mass spectrometer in proteomics.

Comparative MS/MS studies of singly and doubly charged electrospray ionization (ESI) and matrix-assisted laser desorption/ionization (MALDI) precursor peptide ions are described. The spectra from these experiments have been evaluated with particular emphasis on the data quality for subsequent data processing and protein/amino acid sequence identification. It is shown that, once peptide ions are formed by ESI or MALDI, their charge state, as well as the collision energy, is the main parameter determining the quality of collision-induced dissociation (CID) MS/MS fragmentation spectra of a given peptide. CID-MS/MS spectra of singly charged peptides obtained on a hybrid quadrupole orthogonal time-of-flight mass spectrometer resemble very closely spectra obtained by matrix-assisted laser desorption/ionization post-source decay time-of-flight mass spectrometry (MALDI-PSD-TOFMS). On the other hand, comparison of CID-MS/MS spectra of either singly or doubly charged ion species shows no dependence on whether ions have been formed by ESI or MALDI. This observation confirms that, at the time of precursor ion selection, further mass analysis is effectively decoupled from the desorption/ionization event. Since MALDI ions are predominantly formed as singly charged species and ESI ions as doubly charged, the associated difference in the spectral quality of MS/MS spectra as described here imposes direct consequences on data processing, database searching using ion fragmentation data, and de novo sequencing when ionization techniques are changed.     

Towards a full reference library of MS(n) spectra. Testing of a library containing 3126 MS2 spectra of 1743 compounds

A library consisting of 3766 MS(n) spectra of 1743 compounds, including 3126 MS2 spectra acquired mainly using ion trap (IT) and triple-quadrupole (QqQ) instruments, was composed of numerous collections/sources. Ionization techniques were mainly electrospray ionization and also atmospheric pressure chemical ionization and chemical ionization. The library was tested for the performance in identification of unknowns, and in this context this work is believed to be the largest of all known tests of product-ion mass spectral libraries. The MS2 spectra of the same compounds from different collections were in turn divided into spectra of 'unknown' and reference compounds. For each particular compound, library searches were performed resulting in selection by taking into account the best matches for each spectral collection/source. Within each collection/source, replicate MS2 spectra differed in the collision energy used. Overall, there were up to 950 search results giving the best match factors and their ranks in corresponding hit lists. In general, the correct answers were obtained as the 1st rank in up to 60% of the search results when retrieved with (on average) 2.2 'unknown' and 6.2 reference replicates per compound. With two or more replicates of both 'unknown' and reference spectra (the average numbers of replicates were 4.0 and 7.8, respectively), the fraction of correct answers in the 1st rank increased to 77%. This value is close to the performance of established electron ionization mass spectra libraries (up to 79%) found by other workers. The hypothesis that MS2 spectra better match reference spectra acquired using the same type of tandem mass spectrometer (IT or QqQ) was neither strongly proved nor rejected here. The present work shows that MS2 spectral libraries containing sufficiently numerous different entries for each compound are sufficiently efficient for identification of unknowns and suitable for use with different tandem mass spectrometers.     

Analysis of the stochastic variation in LTQ single scan mass spectra

A better understanding of the scan-to-scan signal intensity variation can lead to more sophisticated algorithms for database searching and de novo peptide sequencing using single scan mass spectra. In this study, we systematically studied the variation in relative intensity of m/z values in the single scan product ion mass spectra (MS2) derived from five representative precursor ions (MS1) collected using an LTQ linear ion trap under constant flow direct infusion conditions with peptide concentrations held constant. We applied a matching algorithm based on a pair hidden Markov model to align the peaks from each scan belonging to the same m/z value prior to assessing the signal intensity variation. The most significant single contributor to scan-to-scan signal intensity variation for high abundance ions was centroider error. Our study also showed that the variation in signal intensity is higher than what would be expected if the ion statistics derived from the dual geometry electron multiplier detector followed a Poisson distribution.