BACKGROUND
mass spectrometry  has emerged in recent years as one of the most powerful tools for protein analysis available to proteomics research.
ms-based protein identification strategies typically involve the digestion of protein samples prior to introduction into the mass spectrometer by a site-specific protease such as trypsin.
the derived peptides are subsequently ionized at entry into the mass spectrometer and measured as intact fragment  ions.
subsets of these ions can then be selected on the basis of their mass-to-charge ratio  and subject to further fragmentation, most commonly using collision induced dissociation , in a process known as tandem mass spectrometry .
under the conditions utilized in cid, peptides fragment in predictable patterns resulting in a series of signature spectra.
identification of the protein components in an analyzed sample can then be achieved by correlating the observed signature spectra of individual peptides with the predicted ms/ms spectra of the amino acid sequences derived from protein databases such as swiss-prot and trembl .
over the past few years, computer-assisted database searching using mass spectrometry data has become the standard method for high-throughput protein identification.
unsurprisingly, the performance of computer search algorithms, for example sequest , mascot  and others, has a dramatic influence on the accuracy and reliability of the protein identification process  <cit> .
in general terms, such algorithms use a built-in fragmentation model to construct theoretical fragmentation spectra for candidate peptides derived from databases, and then evaluate the match of these theoretical spectra with observed spectra from ms/ms experiments using defined scoring criteria.
the candidate peptide whose predicted fragmentation spectra best matches the experimental ms/ms spectra is selected as representing the true identity of the analyzed peptide  <cit> .
unfortunately, the performance of computer algorithms currently available is still less than ideal.
generally, these algorithms tend to only utilize the positional information  contained in ms/ms spectra, whereas fail to systematically incorporate the additional intensity information available in the same spectra.
the intensity information is usually applied in an indirect way, for example in mascot, peaks are selected based on intensity for peptide matching, and in sequest peaks for y/b ions are supposed to be higher than peaks for other ions.
previous published work indicates that some efforts have been made to try to utilize spectrum intensity information more effectively  <cit> , but they have predominantly focused on the design of better scoring methods.
furthermore, the application of these previous studies was limited by the oversimplification of the peptide fragmentation models which were applied to construct the theoretical spectra.
this is probably due to our insufficient understanding of the complex mechanisms involved in peptide fragmentation during ms/ms analysis, which makes accurate prediction of spectra intensities in ms/ms very difficult.
recently, a number of research groups have proposed novel fragmentation models in attempts to better understand the mechanisms involved in ms/ms.
for example, wysocki, et al.
proposed the "mobile proton" hypothesis in which protons added to a peptide can transfer along its backbone from the initial site of protonation and subsequently induce fragmentation  <cit> .
according to the hypothesis, peptides can be classified as "mobile" or "non-mobile" by the ratio of charge to arginine number.
they also statistically examined the effect of particular amino acid residues such as asparagine, proline and histidine on fragmentation patterns, with the aim of deducing rules for the influence of these residues on spectra intensities  <cit> .
the "mobile proton" model was later expanded by kapp, et al.
into the "relative mobile proton"  model  <cit> , in which peptides are classified as "mobile", "non-mobile" and "partial mobile" based on their charge number and basic residue number.
schutz, et al.
used a linear model based on rmp hypothesis to calculate the influence of sequence context effects on fragmentation  <cit> .
a kinetic model was developed by zhang  <cit>  to simulate the fragmentation process of a peptide undergoing low-energy cid, and further used to predict the spectra intensity patterns of given peptides.
machine-learning techniques such as bayesian decision trees have also been used to investigate peptide fragmentation behaviour  <cit>  and from this work a group of features that may influence peptide fragmentation have been proposed.
this was the first attempt to our knowledge to systematically utilize intensity information of peptide fragmentation.
however, the machine learning approach used in  <cit>  discovered only a limited number of features to have significant effect on peptide fragmentation and many of these features have already been revealed by other researchers, for example, the presence of basic residues in a peptide sequence, the charge state of the peptide and the presence of proline residue in peptide sequence, etc  <cit> .
whether many other putative determinants are of relevance, and the extent of their influence on fragmentation, is still in question.
in this work, we present a probabilistic machine-learning approach designed to analyze the intensity information contained in ms/ms data, with the aim of developing a better understanding of the rules involved in peptide fragmentation events.
a library of peptide-relevant features as listed in table  <dig> was examined and a score was assigned to each feature to represent the magnitude of its influence on peptide fragmentation.
this information was then used to develop a more sophisticated model to predict the intensity patterns of spectra generated in ms/ms with the expectation that this will improve the reliability of peptide identification.
overall, we attempted to find answers for three basic questions: what factors influence peptide fragmentation during cid?
what is the relationship between the features that influence peptide fragmentation and the resulting intensity pattern of fragmentation spectra?
and finally, how can we accurately predict the spectrum intensity pattern of a given peptide and use this information to improve peptide identification?
all features listed above are supposed to exert influences on the gas-phase fragmentation of peptides.
they are subject to further examination by the bayesian neural network model.
RESULTS
the experimental design following development of our peptide fragmentation model comprised two phases: a feature selection stage for the determination of peptide characteristics that have significant influence on fragmentation, and a model development stage that trained a bayesian neural network with features identified from the first stage.
the performance of the model was then tested by using it to predict spectra intensity patterns for given peptides and subsequently compared with experimental data.
different data and data filtering algorithms were applied during the different phases.
experiment stage 1
in this part of the study, ms/ms spectra data as described in  <cit>  were acquired from wysocki vh.
the intensity information contained within the spectra was then used to verify a library of features that are supposed to influence peptide fragmentation .
the values of relevant amino acid properties that were used for calculating these features can be found in table  <dig>  this feature set is a modified version of what was used by elias, et al.
<cit> .
we aimed to determine a group of features that genuinely influence the intensity patterns of ms/ms spectra.
for this purpose, a bayesian neural network model was developed.
the structure of the network model is illustrated in figure  <dig> and more details can be found in the method section of the paper.
the values of different peptide property used in the study are listed in the table.
values for all the features listed in table  <dig> are calculated with these property values during network training.
the values for mass, hydrophobicity, helicity and basicity are cited from  <cit>  and the values for pi are cited from
in brief the data comprised peptide ms/ms spectra from two micro organisms, shewanella oneidensis and deinococcus radiodurans.
the datasets were derived using lc/ms analysis with ion trap instrumentation .
peptide sequences were assigned to these spectra using the sequest algorithm with a minimum xcorr score of  <dig>  for peptides with molecular weight <  <dig>  and  <dig>  for all other peptides.
using the same chromatographic conditions, accurate masses of the precursor ions detected at the same retention time by ft-icr were used to confirm the assigned sequences.
finally, a total of  <dig> spectra of unique sequence and charge state  were acquired and subject to further analyses.
in our work, we wished to analyse only spectra with non-biased peptide intensities so that the genuine influences of all the features can be determined.
for this purpose, the following filtering criteria were applied to the available  <dig> spectra:
<dig>  only doubly charged peptide spectra were retained for the study.
<dig>  for a given peptide, the intensities of detected b/y ions  according to the assigned sequence, should be no less than 25% of the total intensities of all peaks within the particular spectrum.
this criterion came from our belief that a correctly identified peptide should be able to explain all peaks in the corresponding spectra reasonably well.
accordingly all spectra with this correlation lower than an arbitrary threshold were considered to be either mismatches or biased spectra due to undetected degradation/modification events.
<dig>  for a given peptide, the total intensities of the detected b/y ions should be no less than the intensity of the parent ion of the peptide.
application of this criterion was intended to ensure that all selected peptides are fully fragmented.
<dig>  finally, all candidate peptides were classified according to the "relative mobile proton"  hypothesis  <cit> .
applying the rmp model as a classification criterion enables us to analyze peptides with different relative mobility separately, and also makes it easier for the machine learning algorithm to identify correct rules involved in peptide fragmentation.
as a result, a total number of  <dig> spectra were analysed in this study, comprising  <dig> mobile peptides,  <dig> partially mobile peptides and  <dig> non-mobile peptides.
the length of these peptides ranges from  <dig> to  <dig> with a mean value of  <dig>  the data provided  <dig> input patterns  for the training of our network model.
the first stage of our experiment, a feature selection stage as described above, began with training the bayesian neural network  <dig> epochs using the features listed in table  <dig>  details of network structure and training method can be found in the method section.
the importance of each individual feature was evaluated by updating its "relevance coefficient" α as in eq.
<dig>  the results of coefficients were ranked and normalized, with their mean values defined as 'irrelevance' scores.
the greater an irrelevance score is, the less significant the influence of the corresponding feature is.
the irrelevance scores of each feature of different peptide mobility status are compared in figure  <dig>  and the values of the original scores can be found in additional file  <dig>
as shown in figure  <dig>  the features that influence the fragmentation pathway of peptides vary considerably depending on peptide mobility status.
peptides of mobile or partial-mobile status generally share similar influential feature sets, but for peptides of non-mobile status, the features that influence fragmentation appear to be completely different.
such an observation implies that peptides of mobile- and partial-mobile status do not have fundamental differences in their fragmentation mechanism, whereas non-mobile peptides appear to possess their own unique method of fragmentation.
the results in figure  <dig> indicate that the ion intensity pattern under non-mobile status depends highly on the sequence context of the fragmented peptide.
it is well known that the identities of residues at the either side of a cleavage site play a very important role in determining whether cleavage can occur at this site, and the extent of this cleavage; but their influence are especially prominent for non-mobile peptides, who's spectra are often observed to be dominated by a limited number of ions.
for mobile and partial-mobile peptides, however, fragmentation pathways appear to be determined by a mixture of factors including sequence context, position of cleavage site, mass and length of fragmented peptide, and many others.
the results show that cleavage is more likely to occur at the middle of a peptide rather than at the two ends, as mentioned before by kapp, et al.
<cit> .
we speculate that the specificity of tryptic digestion may contribute to this.
it is also conceivable that the low mass cut-off inherent in ion trap mass spectrometers play a role in this position-selective phenomenon.
it is observed that this phenomenon is less significant for non-mobile peptides, most likely because of the dominant residue-specific fragmentation pathway.
our analysis also reveals that the presence of basic residues can hinder fragmentation at peptide bonds close to them, as reported in other publications  <cit> .
the influence of individual residue will be discussed in details in the next section.
it does not appear from our results that the basic nature of specific residues can influence the fragmentation pathway directly.
although the presence of basic residues within a peptide can result in marked changes in spectra intensity patterns, the basic nature of a particular residue  appears to have little relevance to the fragmentation pattern .
however, the basic characteristic of the whole peptide  does appear to play an important role in fragmentation irrespective of peptide mobility status, and the basic characteristics of fragmented y ions  can influence peptides of mobile and partial-mobile status.
in general, the tendency of amino acid residues to contribute to the helicity nature of a peptide correlates with medium to high irrelevance scores, indicating that these characteristics do not have significant influence on peptide fragmentation, especially for non-mobile peptides .
specific hydrophobicity-related features , however, appear to be important in the fragmentation of both mobile and partial-mobile peptides, but they show little influence on peptides of non-mobile status.
to the best of our knowledge there is no published theory suggesting a mechanism to explain how peptide hydrophobicity may influence fragmentation events in ms/ms, and, we are unsure whether our results stem from a causal relationship or simply a numerical correspondence.
indeed, this may be a topic worthy of future study.
the pi values of residues show little direct influence on fragmentation of peptides .
the features "number of basic residues in the whole peptide"  and "number of basic residues in the fragmentation ion"  were unsurprisingly ranked as having little influence on mobile peptides, because in the great majority of cases doubly charged mobile peptides contain only one single basic residue, which will be located at the c-terminus given the sequence specificity of trypsin.
accordingly these features are of little relevance for mobile peptides.
in contrast, these two features do appear to be influential on peptides under other mobility status because variable numbers of basic residues are usually present in those cases.
it is also apparent that the distance from the fragmented bond to basic residues has little influence on fragmentation pathways .
this feature indeed appears to influence mobile peptides, but such an effect is more likely to be a numerical correspondence only, because the sole basic residue in a doubly charged mobile peptide is located at its c-terminus, making this feature effectively synonymous with the feature "distance from fragmented bond to peptide c-terminus"  that has been shown to influence peptide fragmentation.
finally, we find that the mass and length of fragmented ions/whole peptide can influence the overall fragmentation pattern .
comparatively, the ratio of mass/length are more influential  than absolute values of the two.
this result agrees with findings reported elsewhere  <cit> .
many studies have been conducted to find out how the presence of a particular residue influences the subsequent fragmentation pathway of a whole peptide.
a series of rules has been derived from both statistical analysis and manual interpretation of ms/ms spectra  <cit> .
in our model, every residue has  <dig> separate nodes which represent its presence on the n- or c-terminus of a peptide bond.
we are able to determine the influence of each residue by evaluating the weight values assigned to these nodes.
the results are illustrated in figure  <dig>  as can be seen, many of the defined features appear to influence fragmentation, and most of them conform to the established rules.
this correlation lends credence to the effectiveness of our approach, and supports the validity of the influence of the features as we suggest above.
when a free proton is available within a peptide , we unsurprisingly find that proline  has a significant influence on fragmentation.
as has previously been extensively documented  <cit> , proline markedly enhances cleavage at its n-terminal peptide bond while greatly inhibiting c-terminal cleavage.
conversely, aspartic acid  and glutamic acid  residues appear to inhibit cleavage at their n-termini, and similarly, asparagine  is found to have the same inhibitory effect on peptide fragmentation but to a less significant degree.
isoleucine  and valine  are found to promote c-terminal cleavage most, whereas glycine  and asparagine  residues have the greatest inhibitory effect  on cleavage at the c-terminus.
however with non-mobile peptides, for example those containing multiple arginine  residues, protons are sequestered by the basic amino acids, and as a result the peptide fragments in a totally different manner .
in this situation proline still has the greatest influence on cleavage on n- terminal cleavage, but in comparison to the situation in mobile peptides, this effect is much reduced.
aspartic acid is now the most influential residue in respect to enhanced c-terminal cleavage , although its ability to inhibit cleavage at its n-terminal peptide bond is reduced.
it is clear from the figure that the influence of aspartic acid is almost twice as much as that of proline, so even if they appear in the same peptide, the resulting spectra will be dominated by ions derived from aspartic acid-derived fragmentation.
glutamic acid  favours peptide cleavage at its c-terminus, a characteristic which probably results from the presence of a similar side chain to that of aspartic acid.
glycine-dependent inhibition of cleavage at its c-terminus is observed in all mobility status.
arginine  is observed to strongly promote cleavage at its c-terminus, and the other two basic residues lysine  and histidine  also present the same favour but in a less significant way.
the rules defined above have also been reported previously by wysocki group in work using a statistical method and the same ms/ms spectra dataset  <cit> , and by zhang, using his kinetic model  <cit> .
we also observed a number of novel peptide sequence-context effects.
firstly arginine  residues show a markedly inhibited cleavage at their n-termini in non-mobile peptides.
secondly histidine  appears to favour cleavage at its n-terminus, and such effect is observed in all mobility status.
besides these, previous studies have proposed that leucine  residues promote cleavage to their c-terminal peptide bonds, irrespective of the mobility status of the peptide  <cit> .
this effect is not apparent from our study, with the presence of leucine only having a relatively minor effect  on c-terminal cleavage
in the classical proton mobility theory peptides are classified into  <dig> distinct groups, designated as either mobile or non-mobile, according to the number of arginine residues present within the peptide.
in addition, kapp, et al.
<cit>  have since proposed another class: an intermediate or 'partial' mobility state.
we have also analyzed the fragmentation behaviour of peptides belonging to this third mobility class, and the results are indicated in figure 3-b.
we find that peptides falling into this notional group fragment according to a combination of rules predominant only to either mobile or non-mobile peptides.
effectively the mechanism of fragmentation in partially mobile peptides appears to obey a hybrid rule set.
in this rule set, the influence of residues on fragmentation at their n-terminal peptide bonds is similar to that for peptides of mobile status, in which proline has a dominant enhancing effect, and aspartic acid and glutamic acid inhibit cleavage most.
in marked contrast, the influence of amino acid residues on peptide bonds at their c-terminus more closely resembles that occurring in peptides of non-mobile status, where aspartic acid has the most profound effect.
at the same time, isoleucine and valine enhance cleavage at their c-terminal peptide bonds in partial mobility peptides, as they do in mobile status peptides.
exceptionally, lysine  is observed to enhance c-terminal cleavage in partially mobile peptides.
such an effect was not observed under any other mobility status.
glycine and proline have the most marked inhibitory effect on c-terminal cleavage as they do in both other peptide mobility groups.
it is worth noting that in the earlier work of elias, et al.
<cit>  using a bayesian decision tree method, a similar feature set was examined to analyse the influence of each component on peptide fragmentation.
in that study however, only the 'proline effect' was observed, and the influences of other residues were suggested to be insignificant.
by revealing a considerably larger set of valid fragmentation rules using a similar feature set, it appears that our machine learning model has abetter learning capacity and is capable of identifying more subtle, yet nevertheless significant differences, in the contribution of different amino acid residues to peptide fragmentation during cid.
having determined the irrelevance scores for all features examined, a new feature set can be defined containing only those found to markedly influence peptide fragmentation.
to this end, we sequentially discarded the features with highest scores as listed in figure  <dig>  and then retrained the network with the reduced feature set.
comparison of the training results for all networks is illustrated in figure  <dig>  taking non-mobile peptides as an example, the training error increases significantly when  <dig> less relevant features are removed, indicating that at most  <dig> features can be removed.
the remaining features are indicated with filled circles in figure  <dig>
experiment stage 2
with the reduced feature set derived in the first stage, a new network was trained to predict the intensity patterns of fragmentation spectra for given peptides.
benefiting from the application of the bayesian theory, the network can not only predict the absolute values of spectra intensities, but also assign variances for the predictions.
the obtained results are thus more robust against noise and system errors that unavoidably appear in the experimental ms/ms data.
details of the prediction method used can be found in the methods section.
a new ms/ms dataset was applied to evaluate the performance of the bayesian intensity model.
the dataset is a controlled dataset containing  <dig> different proteins as described in  <cit> .
the details of how the dataset is generated can be found in the original paper.
there are in total  <dig> doubly charged spectra that have been verified to be correctly identified.
applying the same filtering method as described in experimental stage  <dig>  we finally obtained  <dig> doubly charged peptides for model testing.
the theoretical spectra of these peptides were predicted by the bayesian intensity model and then compared with the experimental counterparts to evaluate the accuracy of the model.
in order to compare an experimental spectrum with its predicted counterpart, a score capable of evaluating the similarity of two spectra has to be defined.
as described in the method section, it is assumed that the log-transformed intensities of a given spectrum are gaussian distributed with mean values and variances as predicted by our model.
accordingly, experimental spectra were normalized using the following method:
<dig>  all peaks related to parent ions  are removed.
<dig>  divide each spectrum with its intensity of total ion current  and then times  <dig>
<dig>  log-transformed and then normalized to  <cit>  scale.
it is then straightforward to design the following scoring system  to measure the degree of similarity between the two spectra:
score=∑k=1npexp⁡{−22σpk2}
where np is the number of peptide bonds within peptide p, ipk-predict is the predicted mean intensity value of the peak at peptide bond k, ipk-real is the observed intensity of the peak at peptide bond k in the experimental spectrum, and σpk is the standard deviation  of the peak intensity at peptide bond k predicted by the intensity model.
the more similar the two spectra, the higher the resulting score is.
an example of spectrum predicted by the bayesian intensity model can be found in figure  <dig> for the peptide gysfvttaer.
the prediction for this peptide achieves one of the highest scores .
it can be seen that the predicted spectrum matches its experimental counterpart very well, and the small differences between the two spectra are well within the variance.
further examples using peptides of different mobility status can be found in additional file  <dig>
similarly, the prediction for peptide vlypndnffegk is illustrated in figure  <dig>  this peptide attained one of the lowest scores , indicating a probable failure of spectrum prediction.
indeed, as can be seen in figure 6-c, none of the peaks lie within the expected ranges.
it is apparent from figure 6-a that the experimental spectrum of the peptide is dominated by the y <dig> ion resulting from cleavage at the y-p bond of the peptide, the other ions of the expected y-ion series are either very low or even below the level of detection.
this pattern is characteristic of the type of spectrum that often leads to random  matches in database searching using current m/z based peptide identification algorithms.
however, as shown in figure 6-b, our model did correctly predict the general pattern of the spectrum, i.e.
that y <dig> and y <dig> are the highest two peaks and the others peaks are lower and of relatively equal height.
the experimental spectrum therefore represents a greatly exaggerated version of the predicted pattern.
the similarity scores were firstly calculated for spectra predicted by the bayesian intensity model as illustrated in figure  <dig>  the same scores were subsequently recalculated with intensity information excluded, i.e.
after assigning the same intensity value to each peak within a spectrum.
such an approach, using intensity-free spectra, is typical of most peptide identification algorithms in current use, e.g.
sequest.
in order to compare spectra with/without intensity information on an impartial basis, the similarity scores for intensity-free spectra were firstly calculated using the same variance values as used for scores with intensity information, and then recalculated with the influence of variance eliminated .
the former case was illustrated in figure 7-a and the latter one in figure 7-b.
it can be clearly seen that the scores derived using intensity information are consistently higher than those derived without.
this result indicates that our network model can accurately predict fragmentation spectra for given peptides, and the predicted spectra fit experimental spectra much better than those generated using intensity-free information.
in order to further validate our bayesian intensity model, programs for kinetic intensity model published by zhang  <cit>  were used to predict the intensity patterns of the same test dataset as mentioned above.
the differences between experimental spectra and predicted counterparts are calculated and compared with those from our bayesian intensity model.
it should be noted that zhang's kinetic model is able to predict intensity of b, y ions and their degradations.
so only y related ions were pick out for comparison.
at the same time, the variance information predicted by the bayesian model was also ignored.
as listed in table  <dig>  the bayesian model has a smaller prediction error in  <dig> spectra out of the total  <dig>  showing a slightly higher accuracy than the kinetic model.
however, the mean values and sd values of prediction errors from the two approaches are rather close.
it is reasonable to conclude that the two approaches have similar performance in predicting spectra intensity values for given peptides, and our bayesian intensity model can potentially be more informative because extra variances can be assigned to the predictions to tolerate the prediction error.
the intensity patterns of  <dig> test data are predicted with zhang's kinetic model and our bayesian intensity model.
the accuracy of predictions and mean/sd of prediction errors are listed in the table.
discussion
in this work, a novel bayesian neural network approach was applied to examine features that were thought to influence peptide fragmentation.
the benefit of this approach includes making the features numerically analysable so that large number of features regarding various characteristics of fragmented peptides can be compared directly at one time.
in the experiment stage  <dig>  a new network was trained to predict intensity patterns of fragmentation spectra for given peptides.
only a limited number of features with significant influence were applied in this stage, and the others were discarded.
it is worth noting that the discarded features are not necessarily irrelevant to peptide fragmentation.
indeed, they may still influence fragmentation pathways but in a less significant or indirect way.
however they were discarded given that results indicated that the accuracy of spectra intensity predictions was not significantly affected by the elimination of these features.
the ms/ms data used in this work is dominated by spectra from mobile and partial-mobile peptides, whereas the number of non-mobile data is relatively small.
this is mainly because these peptides were identified by sequest, who rely heavily on m/z information to make identification.
non-mobile peptides, unfortunately, usually fail to present enough amounts of ions in the spectra, and therefore fail to be identified.
it is worth noting that although many reasonable rules in fragmentation have been derived for non-mobile peptides, it is very likely that some important fragmentation rules are missing.
to the best of our knowledge, most algorithms available for peptide identification to date make their identification based on spectra m/z information only.
the importance of spectra intensity information on peptide identification have been realised by many researchers, but successful applications in published literature are still rare:  <cit>  is the only one to the best of our knowledge, and it is simply an application of zhang's kinetic model.
this phenomenon is partially because it is difficult to predict intensity patterns of fragmentation spectra accurately, and even if intensity patterns are successfully predicted, how to compare the predicted spectra with experimental ones still remains a problem.
the large volumes of noise and system errors inevitably contained in experimental spectra make it difficult to apply conventional comparison methods to evaluate the true degree of similarity between spectra.
our network model provided a good way to solve the problem by assigning variances to the predictions to obtain certain degree of tolerance to the fluctuation of spectra intensities.
in this study, we proposed a scoring method  that combined the predicted variance information to compare two spectra under total ion current normalization.
this method worked well in validating the intensity model and could be used for peptide identification.
however, there were still cases in which the predicted spectra failed to match the experimental counterpart even if the general pattern of spectra intensity was predicted correctly .
future work will involve development of more robust scoring methods to further improve the performance of our intensity model in peptide identification, and allow for the effect of post-translational modifications  on peptide fragmentation pathways.
the influence of ptms is unpredictable at present, as modified peptides may fragment in a different manner to unmodified ones, thereby making predictions using current fragmentation models less reliable.
it is important to note that the intensity pattern of ms/ms spectra for the same peptide can differ depending on the nature of the mass spectrometer platform used for analysis, e.g.
whether it is an ion trap or a q-tof system etc., and the method of ionisation employed .
in this situation, the intensity model needs to be retrained each time to adapt to the different machine types and peptide dissociation methods
CONCLUSIONS
in this work, we have shown that the intensity patterns of fragmentation spectra are informative and can be used to analyze the influence of various characteristics of fragmented peptides on their fragmentation pathway.
the features with significant influence can be used in turn to predict spectra intensities given the sequences of peptides.
it has been demonstrated that the intensity pattern of fragmentation spectra predicted by our model fits experimental data reasonably well.
it is suggested that such intensity predictions can be used with current peptide and protein identification algorithms to make them more reliable in high-throughput proteomics experiments.

 mass spectrometry has emerged in recent years as one of the most powerful tools for protein analysis available to proteomics research. ms-based protein identification strategies typically involve the digestion of protein samples prior to introduction into the mass spectrometer by a site-specific protease such as trypsin. the derived peptides are subsequently ionized at entry into the mass spectrometer and measured as intact fragment ions. subsets of these ions can then be selected on the basis of their mass-to-charge ratio and subject to further fragmentation, most commonly using collision induced dissociation, in a process known as tandem mass spectrometry. under the conditions utilized in cid, peptides fragment in predictable patterns resulting in a series of signature spectra. identification of the protein components in an analyzed sample can then be achieved by correlating the observed signature spectra of individual peptides with the predicted ms/ms spectra of the amino acid sequences derived from protein databases such as swiss-prot and trembl. over the past few years, computer-assisted database searching using mass spectrometry data has become the standard method for high-throughput protein identification. unsurprisingly, the performance of computer search algorithms, for example sequest, mascot and others, has a dramatic influence on the accuracy and reliability of the protein identification process. in general terms, such algorithms use a built-in fragmentation model to construct theoretical fragmentation spectra for candidate peptides derived from databases, and then evaluate the match of these theoretical spectra with observed spectra from ms/ms experiments using defined scoring criteria. the candidate peptide whose predicted fragmentation spectra best matches the experimental ms/ms spectra is selected as representing the true identity of the analyzed peptide. unfortunately, the performance of computer algorithms currently available is still less than ideal. generally, these algorithms tend to only utilize the positional information contained in ms/ms spectra, whereas fail to systematically incorporate the additional intensity information available in the same spectra. the intensity information is usually applied in an indirect way, for example in mascot, peaks are selected based on intensity for peptide matching, and in sequest peaks for y/b ions are supposed to be higher than peaks for other ions. previous published work indicates that some efforts have been made to try to utilize spectrum intensity information more effectively, but they have predominantly focused on the design of better scoring methods. furthermore, the application of these previous studies was limited by the oversimplification of the peptide fragmentation models which were applied to construct the theoretical spectra. this is probably due to our insufficient understanding of the complex mechanisms involved in peptide fragmentation during ms/ms analysis, which makes accurate prediction of spectra intensities in ms/ms very difficult. recently, a number of research groups have proposed novel fragmentation models in attempts to better understand the mechanisms involved in ms/ms. for example, wysocki, et al. proposed the "mobile proton" hypothesis in which protons added to a peptide can transfer along its backbone from the initial site of protonation and subsequently induce fragmentation. according to the hypothesis, peptides can be classified as "mobile" or "non-mobile" by the ratio of charge to arginine number. they also statistically examined the effect of particular amino acid residues such as asparagine, proline and histidine on fragmentation patterns, with the aim of deducing rules for the influence of these residues on spectra intensities. the "mobile proton" model was later expanded by kapp, et al. into the "relative mobile proton" model, in which peptides are classified as "mobile", "non-mobile" and "partial mobile" based on their charge number and basic residue number. schutz, et al. used a linear model based on rmp hypothesis to calculate the influence of sequence context effects on fragmentation. a kinetic model was developed by zhang to simulate the fragmentation process of a peptide undergoing low-energy cid, and further used to predict the spectra intensity patterns of given peptides. machine-learning techniques such as bayesian decision trees have also been used to investigate peptide fragmentation behaviour and from this work a group of features that may influence peptide fragmentation have been proposed. this was the first attempt to our knowledge to systematically utilize intensity information of peptide fragmentation. however, the machine learning approach used in discovered only a limited number of features to have significant effect on peptide fragmentation and many of these features have already been revealed by other researchers, for example, the presence of basic residues in a peptide sequence, the charge state of the peptide and the presence of proline residue in peptide sequence, etc. whether many other putative determinants are of relevance, and the extent of their influence on fragmentation, is still in question. in this work, we present a probabilistic machine-learning approach designed to analyze the intensity information contained in ms/ms data, with the aim of developing a better understanding of the rules involved in peptide fragmentation events. a library of peptide-relevant features as listed in table was examined and a score was assigned to each feature to represent the magnitude of its influence on peptide fragmentation. this information was then used to develop a more sophisticated model to predict the intensity patterns of spectra generated in ms/ms with the expectation that this will improve the reliability of peptide identification. overall, we attempted to find answers for three basic questions: what factors influence peptide fragmentation during cid? what is the relationship between the features that influence peptide fragmentation and the resulting intensity pattern of fragmentation spectra? and finally, how can we accurately predict the spectrum intensity pattern of a given peptide and use this information to improve peptide identification? all features listed above are supposed to exert influences on the gas-phase fragmentation of peptides. they are subject to further examination by the bayesian neural network model.  the experimental design following development of our peptide fragmentation model comprised two phases: a feature selection stage for the determination of peptide characteristics that have significant influence on fragmentation, and a model development stage that trained a bayesian neural network with features identified from the first stage. the performance of the model was then tested by using it to predict spectra intensity patterns for given peptides and subsequently compared with experimental data. different data and data filtering algorithms were applied during the different phases. experiment stage 1 in this part of the study, ms/ms spectra data as described in were acquired from wysocki vh. the intensity information contained within the spectra was then used to verify a library of features that are supposed to influence peptide fragmentation. the values of relevant amino acid properties that were used for calculating these features can be found in table this feature set is a modified version of what was used by elias, et al.. we aimed to determine a group of features that genuinely influence the intensity patterns of ms/ms spectra. for this purpose, a bayesian neural network model was developed. the structure of the network model is illustrated in figure and more details can be found in the method section of the paper. the values of different peptide property used in the study are listed in the table. values for all the features listed in table are calculated with these property values during network training. the values for mass, hydrophobicity, helicity and basicity are cited from and the values for pi are cited from in brief the data comprised peptide ms/ms spectra from two micro organisms, shewanella oneidensis and deinococcus radiodurans. the datasets were derived using lc/ms analysis with ion trap instrumentation. peptide sequences were assigned to these spectra using the sequest algorithm with a minimum xcorr score of for peptides with molecular weight < and for all other peptides. using the same chromatographic conditions, accurate masses of the precursor ions detected at the same retention time by ft-icr were used to confirm the assigned sequences. finally, a total of spectra of unique sequence and charge state were acquired and subject to further analyses. in our work, we wished to analyse only spectra with non-biased peptide intensities so that the genuine influences of all the features can be determined. for this purpose, the following filtering criteria were applied to the available spectra: only doubly charged peptide spectra were retained for the study. for a given peptide, the intensities of detected b/y ions according to the assigned sequence, should be no less than 25% of the total intensities of all peaks within the particular spectrum. this criterion came from our belief that a correctly identified peptide should be able to explain all peaks in the corresponding spectra reasonably well. accordingly all spectra with this correlation lower than an arbitrary threshold were considered to be either mismatches or biased spectra due to undetected degradation/modification events. for a given peptide, the total intensities of the detected b/y ions should be no less than the intensity of the parent ion of the peptide. application of this criterion was intended to ensure that all selected peptides are fully fragmented. finally, all candidate peptides were classified according to the "relative mobile proton" hypothesis. applying the rmp model as a classification criterion enables us to analyze peptides with different relative mobility separately, and also makes it easier for the machine learning algorithm to identify correct rules involved in peptide fragmentation. as a result, a total number of spectra were analysed in this study, comprising mobile peptides, partially mobile peptides and non-mobile peptides. the length of these peptides ranges from to with a mean value of the data provided input patterns for the training of our network model. the first stage of our experiment, a feature selection stage as described above, began with training the bayesian neural network epochs using the features listed in table details of network structure and training method can be found in the method section. the importance of each individual feature was evaluated by updating its "relevance coefficient" α as in eq. the results of coefficients were ranked and normalized, with their mean values defined as'irrelevance' scores. the greater an irrelevance score is, the less significant the influence of the corresponding feature is. the irrelevance scores of each feature of different peptide mobility status are compared in figure and the values of the original scores can be found in additional file as shown in figure the features that influence the fragmentation pathway of peptides vary considerably depending on peptide mobility status. peptides of mobile or partial-mobile status generally share similar influential feature sets, but for peptides of non-mobile status, the features that influence fragmentation appear to be completely different. such an observation implies that peptides of mobile- and partial-mobile status do not have fundamental differences in their fragmentation mechanism, whereas non-mobile peptides appear to possess their own unique method of fragmentation. the results in figure indicate that the ion intensity pattern under non-mobile status depends highly on the sequence context of the fragmented peptide. it is well known that the identities of residues at the either side of a cleavage site play a very important role in determining whether cleavage can occur at this site, and the extent of this cleavage; but their influence are especially prominent for non-mobile peptides, who's spectra are often observed to be dominated by a limited number of ions. for mobile and partial-mobile peptides, however, fragmentation pathways appear to be determined by a mixture of factors including sequence context, position of cleavage site, mass and length of fragmented peptide, and many others. the results show that cleavage is more likely to occur at the middle of a peptide rather than at the two ends, as mentioned before by kapp, et al.. we speculate that the specificity of tryptic digestion may contribute to this. it is also conceivable that the low mass cut-off inherent in ion trap mass spectrometers play a role in this position-selective phenomenon. it is observed that this phenomenon is less significant for non-mobile peptides, most likely because of the dominant residue-specific fragmentation pathway. our analysis also reveals that the presence of basic residues can hinder fragmentation at peptide bonds close to them, as reported in other publications. the influence of individual residue will be discussed in details in the next section. it does not appear from our results that the basic nature of specific residues can influence the fragmentation pathway directly. although the presence of basic residues within a peptide can result in marked changes in spectra intensity patterns, the basic nature of a particular residue appears to have little relevance to the fragmentation pattern. however, the basic characteristic of the whole peptide does appear to play an important role in fragmentation irrespective of peptide mobility status, and the basic characteristics of fragmented y ions can influence peptides of mobile and partial-mobile status. in general, the tendency of amino acid residues to contribute to the helicity nature of a peptide correlates with medium to high irrelevance scores, indicating that these characteristics do not have significant influence on peptide fragmentation, especially for non-mobile peptides. specific hydrophobicity-related features, however, appear to be important in the fragmentation of both mobile and partial-mobile peptides, but they show little influence on peptides of non-mobile status. to the best of our knowledge there is no published theory suggesting a mechanism to explain how peptide hydrophobicity may influence fragmentation events in ms/ms, and, we are unsure whether our results stem from a causal relationship or simply a numerical correspondence. indeed, this may be a topic worthy of future study. the pi values of residues show little direct influence on fragmentation of peptides. the features "number of basic residues in the whole peptide" and "number of basic residues in the fragmentation ion" were unsurprisingly ranked as having little influence on mobile peptides, because in the great majority of cases doubly charged mobile peptides contain only one single basic residue, which will be located at the c-terminus given the sequence specificity of trypsin. accordingly these features are of little relevance for mobile peptides. in contrast, these two features do appear to be influential on peptides under other mobility status because variable numbers of basic residues are usually present in those cases. it is also apparent that the distance from the fragmented bond to basic residues has little influence on fragmentation pathways. this feature indeed appears to influence mobile peptides, but such an effect is more likely to be a numerical correspondence only, because the sole basic residue in a doubly charged mobile peptide is located at its c-terminus, making this feature effectively synonymous with the feature "distance from fragmented bond to peptide c-terminus" that has been shown to influence peptide fragmentation. finally, we find that the mass and length of fragmented ions/whole peptide can influence the overall fragmentation pattern. comparatively, the ratio of mass/length are more influential than absolute values of the two. this result agrees with findings reported elsewhere. many studies have been conducted to find out how the presence of a particular residue influences the subsequent fragmentation pathway of a whole peptide. a series of rules has been derived from both statistical analysis and manual interpretation of ms/ms spectra. in our model, every residue has separate nodes which represent its presence on the n- or c-terminus of a peptide bond. we are able to determine the influence of each residue by evaluating the weight values assigned to these nodes. the results are illustrated in figure as can be seen, many of the defined features appear to influence fragmentation, and most of them conform to the established rules. this correlation lends credence to the effectiveness of our approach, and supports the validity of the influence of the features as we suggest above. when a free proton is available within a peptide, we unsurprisingly find that proline has a significant influence on fragmentation. as has previously been extensively documented, proline markedly enhances cleavage at its n-terminal peptide bond while greatly inhibiting c-terminal cleavage. conversely, aspartic acid and glutamic acid residues appear to inhibit cleavage at their n-termini, and similarly, asparagine is found to have the same inhibitory effect on peptide fragmentation but to a less significant degree. isoleucine and valine are found to promote c-terminal cleavage most, whereas glycine and asparagine residues have the greatest inhibitory effect on cleavage at the c-terminus. however with non-mobile peptides, for example those containing multiple arginine residues, protons are sequestered by the basic amino acids, and as a result the peptide fragments in a totally different manner. in this situation proline still has the greatest influence on cleavage on n- terminal cleavage, but in comparison to the situation in mobile peptides, this effect is much reduced. aspartic acid is now the most influential residue in respect to enhanced c-terminal cleavage, although its ability to inhibit cleavage at its n-terminal peptide bond is reduced. it is clear from the figure that the influence of aspartic acid is almost twice as much as that of proline, so even if they appear in the same peptide, the resulting spectra will be dominated by ions derived from aspartic acid-derived fragmentation. glutamic acid favours peptide cleavage at its c-terminus, a characteristic which probably results from the presence of a similar side chain to that of aspartic acid. glycine-dependent inhibition of cleavage at its c-terminus is observed in all mobility status. arginine is observed to strongly promote cleavage at its c-terminus, and the other two basic residues lysine and histidine also present the same favour but in a less significant way. the rules defined above have also been reported previously by wysocki group in work using a statistical method and the same ms/ms spectra dataset, and by zhang, using his kinetic model. we also observed a number of novel peptide sequence-context effects. firstly arginine residues show a markedly inhibited cleavage at their n-termini in non-mobile peptides. secondly histidine appears to favour cleavage at its n-terminus, and such effect is observed in all mobility status. besides these, previous studies have proposed that leucine residues promote cleavage to their c-terminal peptide bonds, irrespective of the mobility status of the peptide. this effect is not apparent from our study, with the presence of leucine only having a relatively minor effect on c-terminal cleavage in the classical proton mobility theory peptides are classified into distinct groups, designated as either mobile or non-mobile, according to the number of arginine residues present within the peptide. in addition, kapp, et al. have since proposed another class: an intermediate or'partial' mobility state. we have also analyzed the fragmentation behaviour of peptides belonging to this third mobility class, and the results are indicated in figure 3-b. we find that peptides falling into this notional group fragment according to a combination of rules predominant only to either mobile or non-mobile peptides. effectively the mechanism of fragmentation in partially mobile peptides appears to obey a hybrid rule set. in this rule set, the influence of residues on fragmentation at their n-terminal peptide bonds is similar to that for peptides of mobile status, in which proline has a dominant enhancing effect, and aspartic acid and glutamic acid inhibit cleavage most. in marked contrast, the influence of amino acid residues on peptide bonds at their c-terminus more closely resembles that occurring in peptides of non-mobile status, where aspartic acid has the most profound effect. at the same time, isoleucine and valine enhance cleavage at their c-terminal peptide bonds in partial mobility peptides, as they do in mobile status peptides. exceptionally, lysine is observed to enhance c-terminal cleavage in partially mobile peptides. such an effect was not observed under any other mobility status. glycine and proline have the most marked inhibitory effect on c-terminal cleavage as they do in both other peptide mobility groups. it is worth noting that in the earlier work of elias, et al. using a bayesian decision tree method, a similar feature set was examined to analyse the influence of each component on peptide fragmentation. in that study however, only the'proline effect' was observed, and the influences of other residues were suggested to be insignificant. by revealing a considerably larger set of valid fragmentation rules using a similar feature set, it appears that our machine learning model has abetter learning capacity and is capable of identifying more subtle, yet nevertheless significant differences, in the contribution of different amino acid residues to peptide fragmentation during cid. having determined the irrelevance scores for all features examined, a new feature set can be defined containing only those found to markedly influence peptide fragmentation. to this end, we sequentially discarded the features with highest scores as listed in figure and then retrained the network with the reduced feature set. comparison of the training results for all networks is illustrated in figure taking non-mobile peptides as an example, the training error increases significantly when less relevant features are removed, indicating that at most features can be removed. the remaining features are indicated with filled circles in figure experiment stage 2 with the reduced feature set derived in the first stage, a new network was trained to predict the intensity patterns of fragmentation spectra for given peptides. benefiting from the application of the bayesian theory, the network can not only predict the absolute values of spectra intensities, but also assign variances for the predictions. the obtained results are thus more robust against noise and system errors that unavoidably appear in the experimental ms/ms data. details of the prediction method used can be found in the methods section. a new ms/ms dataset was applied to evaluate the performance of the bayesian intensity model. the dataset is a controlled dataset containing different proteins as described in. the details of how the dataset is generated can be found in the original paper. there are in total doubly charged spectra that have been verified to be correctly identified. applying the same filtering method as described in experimental stage we finally obtained doubly charged peptides for model testing. the theoretical spectra of these peptides were predicted by the bayesian intensity model and then compared with the experimental counterparts to evaluate the accuracy of the model. in order to compare an experimental spectrum with its predicted counterpart, a score capable of evaluating the similarity of two spectra has to be defined. as described in the method section, it is assumed that the log-transformed intensities of a given spectrum are gaussian distributed with mean values and variances as predicted by our model. accordingly, experimental spectra were normalized using the following method: all peaks related to parent ions are removed. divide each spectrum with its intensity of total ion current and then times log-transformed and then normalized to scale. it is then straightforward to design the following scoring system to measure the degree of similarity between the two spectra: score=∑k=1npexp⁡{−22σpk2} where np is the number of peptide bonds within peptide p, ipk-predict is the predicted mean intensity value of the peak at peptide bond k, ipk-real is the observed intensity of the peak at peptide bond k in the experimental spectrum, and σpk is the standard deviation of the peak intensity at peptide bond k predicted by the intensity model. the more similar the two spectra, the higher the resulting score is. an example of spectrum predicted by the bayesian intensity model can be found in figure for the peptide gysfvttaer. the prediction for this peptide achieves one of the highest scores. it can be seen that the predicted spectrum matches its experimental counterpart very well, and the small differences between the two spectra are well within the variance. further examples using peptides of different mobility status can be found in additional file similarly, the prediction for peptide vlypndnffegk is illustrated in figure this peptide attained one of the lowest scores, indicating a probable failure of spectrum prediction. indeed, as can be seen in figure 6-c, none of the peaks lie within the expected ranges. it is apparent from figure 6-a that the experimental spectrum of the peptide is dominated by the y ion resulting from cleavage at the y-p bond of the peptide, the other ions of the expected y-ion series are either very low or even below the level of detection. this pattern is characteristic of the type of spectrum that often leads to random matches in database searching using current m/z based peptide identification algorithms. however, as shown in figure 6-b, our model did correctly predict the general pattern of the spectrum, i.e. that y and y are the highest two peaks and the others peaks are lower and of relatively equal height. the experimental spectrum therefore represents a greatly exaggerated version of the predicted pattern. the similarity scores were firstly calculated for spectra predicted by the bayesian intensity model as illustrated in figure the same scores were subsequently recalculated with intensity information excluded, i.e. after assigning the same intensity value to each peak within a spectrum. such an approach, using intensity-free spectra, is typical of most peptide identification algorithms in current use, e.g. sequest. in order to compare spectra with/without intensity information on an impartial basis, the similarity scores for intensity-free spectra were firstly calculated using the same variance values as used for scores with intensity information, and then recalculated with the influence of variance eliminated. the former case was illustrated in figure 7-a and the latter one in figure 7-b. it can be clearly seen that the scores derived using intensity information are consistently higher than those derived without. this result indicates that our network model can accurately predict fragmentation spectra for given peptides, and the predicted spectra fit experimental spectra much better than those generated using intensity-free information. in order to further validate our bayesian intensity model, programs for kinetic intensity model published by zhang were used to predict the intensity patterns of the same test dataset as mentioned above. the differences between experimental spectra and predicted counterparts are calculated and compared with those from our bayesian intensity model. it should be noted that zhang's kinetic model is able to predict intensity of b, y ions and their degradations. so only y related ions were pick out for comparison. at the same time, the variance information predicted by the bayesian model was also ignored. as listed in table the bayesian model has a smaller prediction error in spectra out of the total showing a slightly higher accuracy than the kinetic model. however, the mean values and sd values of prediction errors from the two approaches are rather close. it is reasonable to conclude that the two approaches have similar performance in predicting spectra intensity values for given peptides, and our bayesian intensity model can potentially be more informative because extra variances can be assigned to the predictions to tolerate the prediction error. the intensity patterns of test data are predicted with zhang's kinetic model and our bayesian intensity model. the accuracy of predictions and mean/sd of prediction errors are listed in the table. discussion in this work, a novel bayesian neural network approach was applied to examine features that were thought to influence peptide fragmentation. the benefit of this approach includes making the features numerically analysable so that large number of features regarding various characteristics of fragmented peptides can be compared directly at one time. in the experiment stage a new network was trained to predict intensity patterns of fragmentation spectra for given peptides. only a limited number of features with significant influence were applied in this stage, and the others were discarded. it is worth noting that the discarded features are not necessarily irrelevant to peptide fragmentation. indeed, they may still influence fragmentation pathways but in a less significant or indirect way. however they were discarded given that results indicated that the accuracy of spectra intensity predictions was not significantly affected by the elimination of these features. the ms/ms data used in this work is dominated by spectra from mobile and partial-mobile peptides, whereas the number of non-mobile data is relatively small. this is mainly because these peptides were identified by sequest, who rely heavily on m/z information to make identification. non-mobile peptides, unfortunately, usually fail to present enough amounts of ions in the spectra, and therefore fail to be identified. it is worth noting that although many reasonable rules in fragmentation have been derived for non-mobile peptides, it is very likely that some important fragmentation rules are missing. to the best of our knowledge, most algorithms available for peptide identification to date make their identification based on spectra m/z information only. the importance of spectra intensity information on peptide identification have been realised by many researchers, but successful applications in published literature are still rare: is the only one to the best of our knowledge, and it is simply an application of zhang's kinetic model. this phenomenon is partially because it is difficult to predict intensity patterns of fragmentation spectra accurately, and even if intensity patterns are successfully predicted, how to compare the predicted spectra with experimental ones still remains a problem. the large volumes of noise and system errors inevitably contained in experimental spectra make it difficult to apply conventional comparison methods to evaluate the true degree of similarity between spectra. our network model provided a good way to solve the problem by assigning variances to the predictions to obtain certain degree of tolerance to the fluctuation of spectra intensities. in this study, we proposed a scoring method that combined the predicted variance information to compare two spectra under total ion current normalization. this method worked well in validating the intensity model and could be used for peptide identification. however, there were still cases in which the predicted spectra failed to match the experimental counterpart even if the general pattern of spectra intensity was predicted correctly. future work will involve development of more robust scoring methods to further improve the performance of our intensity model in peptide identification, and allow for the effect of post-translational modifications on peptide fragmentation pathways. the influence of ptms is unpredictable at present, as modified peptides may fragment in a different manner to unmodified ones, thereby making predictions using current fragmentation models less reliable. it is important to note that the intensity pattern of ms/ms spectra for the same peptide can differ depending on the nature of the mass spectrometer platform used for analysis, e.g. whether it is an ion trap or a q-tof system etc., and the method of ionisation employed. in this situation, the intensity model needs to be retrained each time to adapt to the different machine types and peptide dissociation methods  in this work, we have shown that the intensity patterns of fragmentation spectra are informative and can be used to analyze the influence of various characteristics of fragmented peptides on their fragmentation pathway. the features with significant influence can be used in turn to predict spectra intensities given the sequences of peptides. it has been demonstrated that the intensity pattern of fragmentation spectra predicted by our model fits experimental data reasonably well. it is suggested that such intensity predictions can be used with current peptide and protein identification algorithms to make them more reliable in high-throughput proteomics experiments. 
useful rules involved in peptide fragmentation are found and subsets of features which have significant influence on fragmentation pathway of peptides are characterised. the intensity patterns of fragmentation spectra are informative and can be used to analyze the influence of various characteristics of fragmented peptides on their fragmentation pathway. a better understanding of the mechanisms involved in gas-phase fragmentation of peptides is essential for the development of more reliable algorithms for high-throughput protein identification using mass spectrometry. current methodologies depend predominantly on the use of derived m/z values of fragment ions, and, the knowledge provided by the intensity information present in ms/ms spectra has not been fully exploited. indeed spectrum intensity information is very rarely utilized in the algorithms currently in use for high-throughput protein identification. in this work, a bayesian neural network approach is employed to analyze ion intensity information present in different ms/ms spectra. the predictions include not only the mean values of spectra intensity but also the variances that can be used to tolerate noises and system biases within experimental ms/ms spectra. an intensity model is built based on the selected features and the model can make an accurate prediction of the intensity patterns for given ms/ms spectra. the features with significant influence can be used in turn to predict spectra intensities. the influence of a library of features on peptide fragmentation is examined under different proton mobility conditions.
 a better understanding of the mechanisms involved in gas-phase fragmentation of peptides is essential for the development of more reliable algorithms for high-throughput protein identification using mass spectrometry. current methodologies depend predominantly on the use of derived m/z values of fragment ions, and, the knowledge provided by the intensity information present in ms/ms spectra has not been fully exploited. indeed spectrum intensity information is very rarely utilized in the algorithms currently in use for high-throughput protein identification.  in this work, a bayesian neural network approach is employed to analyze ion intensity information present in different ms/ms spectra. the influence of a library of features on peptide fragmentation is examined under different proton mobility conditions. useful rules involved in peptide fragmentation are found and subsets of features which have significant influence on fragmentation pathway of peptides are characterised. an intensity model is built based on the selected features and the model can make an accurate prediction of the intensity patterns for given ms/ms spectra. the predictions include not only the mean values of spectra intensity but also the variances that can be used to tolerate noises and system biases within experimental ms/ms spectra.  the intensity patterns of fragmentation spectra are informative and can be used to analyze the influence of various characteristics of fragmented peptides on their fragmentation pathway. the features with significant influence can be used in turn to predict spectra intensities. such information can help develop more reliable algorithms for peptide and protein identification. 