Skip to main content

Calculation of partial isotope incorporation into peptides measured by mass spectrometry



Stable isotope probing (SIP) technique was developed to link function, structure and activity of microbial cultures metabolizing carbon and nitrogen containing substrates to synthesize their biomass. Currently, available methods are restricted solely to the estimation of fully saturated heavy stable isotope incorporation and convenient methods with sufficient accuracy are still missing. However in order to track carbon fluxes in microbial communities new methods are required that allow the calculation of partial incorporation into biomolecules.


In this study, we use the characteristics of the so-called 'half decimal place rule' (HDPR) in order to accurately calculate the partial13C incorporation in peptides from enzymatic digested proteins. Due to the clade-crossing universality of proteins within bacteria, any available high-resolution mass spectrometry generated dataset consisting of tryptically-digested peptides can be used as reference.

We used a freely available peptide mass dataset from Mycobacterium tuberculosis consisting of 315,579 entries. From this the error of estimated versus known heavy stable isotope incorporation from an increasing number of randomly drawn peptide sub-samples (100 times each; no repetition) was calculated. To acquire an estimated incorporation error of less than 5 atom %, about 100 peptide masses were needed. Finally, for testing the general applicability of our method, peptide masses of tryptically digested proteins from Pseudomonas putida ML2 grown on labeled substrate of various known concentrations were used and13C isotopic incorporation was successfully predicted. An easy-to-use script [1] was further developed to guide users through the calculation procedure for their own data series.


Our method is valuable for estimating13C incorporation into peptides/proteins accurately and with high sensitivity. Generally, our method holds promise for wider applications in qualitative and especially quantitative proteomics.


An important aspect of microbial ecology is to link specific microorganisms to microbially-driven processes in the natural environment [2]. Indeed, these questions can be addressed with more accuracy by using isotopically labeled substrates to follow metabolic carbon flux within microbial cultures. Subsequently, the incorporation of the label into the biomass may be followed by analysis of either fatty acids [3], deoxyribonucleic acid (DNA) [4] or ribonucleic acid (RNA) [5]. Unfortunately, those techniques are limited in their ability to resolve low levels of isotopic labeling and their absolute incorporation. Essentially, the information about low level labeling as well as the assignment of distinct incorporation levels to various species is needed in order to understand the metabolic interdependencies within complex microbial communities.

To overcome the drawbacks of other calculation methods, we recently developed a new protein-based stable isotope probing (Protein-SIP) technique [6]. A brief overview of the workflow is displayed in Figure 1. In short, the incorporation of stable e.g.13C and/or15N isotopes from a substrate was used to pinpoint the metabolically active species within a consortium. Different species incorporate different degrees of the stable isotopes into their biomass according to their physiological ability and their activity (Figure 1A), which can be detected on peptide/protein level [7]. After cell harvesting and protein extraction, samples were analyzed by various proteomic and mass-spectrometry techniques. Due to the isotopic abundance of13C and/or15N, several isotopic envelopes can be detected by high-resolution mass spectrometry (MS) [8, 9]. As expected, the isotopologues shifted to a higher mass range due to the incorporation of heavy labeled carbon into the proteins, as shown in the MS-spectrum in figure 1B. The12C-containing peptides are shown at the left hand side of each spectrum. A higher level of incorporation, as shown in the right spectra of part B, indicates a faster growth rate and/or an increased substrate usage of the carbon source of species. In order to assign these activities to different species, two steps have to be applied (i) the light peptides have to be identified by Tandem-MS (mass spectroscopy) and further assigned to the bacterial taxon and (ii) the heavy stable isotope incorporation has to be determined.

Figure 1
figure 1

Schematic overview about the workflow of analysis. (A) Assimilation of heavy stable isotopes into the biomass of various species depends on the turnover and the interaction activity of the species. The incorporation of stable13C isotopes from a substrate can be used to pinpoint the metabolically active species within a consortium. The different incorporations are indicated by various amounts of label (gray-color scale). (B) After cell harvesting and protein extraction, samples were tryptically digested and analyzed by MS. The isotopologues shifted to a higher mass range due to their incorporation of heavy labeled carbon into the proteins. A higher level of incorporation indicates a faster growth rate and/or a primary role in the degradation of the labeled substrate within the food web. (C) The incorporation of heavy stable isotopes into peptides/proteins can be estimated using the HDPR. (D) Peptides can be used to obtain phylogenetic information (in case of unique peptides), for structural analysis, and for physiological information about the actual state of the microbial cells. (E) Based on this information, C-fluxes and food web structures can be elucidated and may further help to reconstruct the interaction of microbial communities.

One method of calculation is to compare12C and13C isotopologue distributions in one MS-spectrum and estimate13C incorporation by calculating the mass difference between the monoisotopic12C mass peak and the highest mass peak of the13C isotopologue. This procedure depends on a proper pre-fractionation, since the corresponding mass peaks can be rather difficult to identify. Alternative calculation methods are necessary to find the matching pairs at low incorporation levels [10].

Therefore, an improved method based on the comparison of theoretical and experimentally determined isotopic envelopes of peptides with known sequences was developed [11, 12]. However, both methods require a priori peptide identification to accurate determine the dynamic labeling of heavy stable isotopes. In order to reliable quantify partial13C incorporation into peptides, we developed a method based on the "half decimal place rule" (HDPR) (Figure 1C). Information about the dynamic incorporation levels in proteins from different microbial species can be used to elucidate the structure and function of the microbial community (Figure 1D). With this information, the carbon fluxes throughout the community may be followed and enables to deduce species interactions and activities (Figure 1E) [13, 14].

In this study, we describe a new and easy-to-apply algorithm for the determination of absolute heavy stable isotope incorporation (13C) into peptides/proteins by taking advantage of high-resolution MS data [15, 16], and the characteristic pattern arising from the decimal residuals of incorporated heavy stable carbon isotopes [17]. Annotated sequence information for measured peptides is no longer necessary for the calculation of13C incorporation that is an important consideration when analyzing uncultivable or incompletely sequenced bacteria. Moreover, we provide an easy-to-apply script set written in the statistical programming language R [18], enabling scientific researchers to calculate13C incorporation into peptides/proteins.

Theoretical background

In the past, mass spectrometers suffered from their low resolution and accuracy, restricting the possibility to predict chemical element contents based on the mass signals. Higher accuracy allows to create exact fitting curves describing the behavior of a mass signal as a function of both the mass and the e.g. sulfur content [19]. With up-coming of high resolution MS devices, a linear relationship between tryptic peptide masses and the decimal residuals (=digits behind the mass' decimal point) has been firstly observed by Mann [20] (for additional references see [19, 21]). Closer investigation of this mass mapping phenomenon by Schmidt et al. [17], led to the definition of the linear relationship of tryptic peptides and the corresponding decimal places, the so-called "half decimal place rule" (HDRP).

This rule declares that the decimal place of a tryptic peptide is near the half of the first digit for tryptic peptides in the range of 500-1,000 Dalton (Da), near the half of the first two digits for the range of 1,000-1,999 Da and again near the half of the first digit for masses from 2,000-3,000 Da. The rule was found to be helpful for e.g. the detection of non-peptide contaminants in mass spectrometric measurements and the recalibration of peptide masses [17]. Since only the digits behind the decimal point are necessary as information in order to identify some artifacts of tryptic peptides.

During the metabolic incorporation of13C amino acid precursors into proteins,13C atoms are gradually replaced12C atoms. This elemental substitution increases each peptide mass by exactly 1.003355 Da (= the difference between12C and13C atomic mass) per substituted carbon atom. In terms of decimal places for peptides, each incorporated heavy carbon atom result in a mass shift of exactly 0.003355 Da. Applying the HDRP to partly or fully labeled peptides, the linear relationship between peptide masses and their decimal places will result in an increase of the corresponding slope. The steepness of the slope increases in proportion with the amount of incorporated heavy stable isotopes. This effect can be used for the exact estimation of the amount of incorporation for any given set of measured peptide masses. The detection of these small mass shifts, however, requires highly accurate measurements. This is now possible with modern mass spectrometers such as a fourier transformation ion cyclotron resonance (FT-ICR), linear ion trap with an Orbitrap (LTQ-Orbitrap) or quadruple time-of-flight (Q-TOF) instrument.


Calculation of isotopic composition of peptides

Processing of reference dataset

Since proteogenic amino acids are universal within all bacterial species, any high-resolution dataset can be used for referencing the theoretical decimal place calculation (as reference data and for calibration). An existing tryptic peptide dataset of Mycobacterium tuberculosis H37Rv was used [22] that originally contained ~4,000 encoding genes of which 3,924 were identified as proteins [23]. The complete FASTA protein dataset is freely available from the Sanger Institute [24]. After a tryptic in-silico digest by MS-Digest [25], peptides between the mass range m/z 300-6,000 were considered for further analysis, resulting in 315,579 peptide fragment sequences. Further, we restricted the dataset to include only those tryptic peptides with a ChemScore ≥ 10 (ChemScore = sum of total free protein binding energy (for further details see [2628]) and without missing cleavage sites or modifications. After this screening, the dataset was reduced to 90,637 remaining peptide sequences containing lengths between 2 and 40 amino acids. The monoisotopic mass of each sequence was calculated with given atomic masses of12C = 12.000000 Da,14N = 14.003074 Da,16O = 15.994915 Da,1H = 1.007825 Da,32S = 31.972071 Da to obtain the masses of 'light' (0 atom %13C incorporation) amino acids. Subsequently, by counting the numbers of carbons for each sequence and replacing the12C mass for the heavier13C (13.003355 Da), a dataset for the theoretical 100 atom %13C incorporation was obtained. The advantage of calculating the weight by counting C atom numbers is that the same procedure can be used to easily calculate a wide range of different13C incorporation levels.

Data classification and modification

Direct plotting of peptide masses against their decimal residuals would result in regular, diagonal linear shaped patterns (Figure 2). With increasing peptide masses, also their decimal residuals increase until they approach values close to 1. With further increasing, the decimal residuals start with values of 0 again. In order to obtain a continuous plot, the points of the second band had to be added on top of the first, and those of the third band on top of the second etc. by subsequently adding 0, 1, 2, etc Dalton to the residuals of each band.

Figure 2
figure 2

A scatter plot of masses from 90,637 peptides ( m/z ) and decimal residuals (digits behind the decimal point) of Mycobacterium tuberculosis within a mass range of m/z 0-5,000. (A) Unlabeled sequences (0 atomic %13C incorporation), (B) fully labeled (100 atomic %13C incorporation). Lines depict the overlap of scatter groups.

However, due to slight deviations from linearity caused by the amino acid residues, an overlap of the bands occurs (indicated by vertical lines in Figure 2A), disabling a direct distinct and straightforward separation. Moreover, increasing substitution of peptides with heavy stable13C isotopes does not only result in an alteration of the slope becoming steeper, but also produces a continuous and gradual shift of the band along the m/z axis. Therefore, we developed a semi-automatic and more flexible approach. Each point cloud was separated by a standard classification method using k-means clustering. Prior to clustering, we transformed the original mass values in relation to their decimal residual values by the following equation,


where PM Trans = transposed peptide masses, PM = peptide masses, and DR = decimal residuals producing a plot as given in Figure 3. The value of 1,800 within the formula (1) was iteratively estimated. Gradually increase of this value made the bands steeper until reaching a maximal vertical position at a value of about 1,800. A value greater than 1,800 made the bands tilt towards the opposite direction. The following classification was conducted using the Hartigan and Wong algorithm for k-means clustering [29] with three (0; 2,000; 4,000 Da) and four (0; 1,600; 3,200; 4,800 Da) pre-set clustering centers for the 'light' (no13C incorporation) and 'heavy' (complete13C substitution) dataset. For the clustering, pre-set center mass values do not have to be overly accurate since the precise numbers are automatically determined during the clustering procedure from actual group means [30]. Finally, the original peptide masses were taken and ranked according to cluster affiliation of the corresponding transformed mass values. In a next step, 0 Da was added to the peptide mass values belonging to the first cluster, 1 Da was added to the values of the second cluster and so on. Plotting these new mass values versus their decimal residuals resulted in the final straight plot (Figure 4).

Figure 3
figure 3

A scatter plot of temporally transposed masses from 90,637 peptides ( m/z ) and decimal residuals of a Mycobacterium tuberculosis dataset (details see text)). Grey dots indicate cluster affiliation after classification by k-means clustering for (A) Unlabeled sequences (0 atomic %13C incorporation), (B) fully labeled (100 atomic %13C incorporation).

Figure 4
figure 4

A scatter plot of 90,637 peptide masses ( m/z ) and their decimal residuals of a Mycobacterium tuberculosis dataset after classification and data modification for unlabeled (0 atom %13C incorporation; light grey) and fully labeled (100 atom %13C incorporation; dark grey) peptide masses resulted. The two lines represent the slopes for the calculation of13C incorporation. The values of the two scatters are for the unlabeled b0% = 5.1357e-4 and b100% = 6.3347e-4 for fully labeled dataset. Numbers in each segment represent Daltons added to decimal residuals to achieve straight plot from figure 3 (for details see text).

Data plotting and linear slope estimation

Following this principle for peptide masses relating to 0 atom % and 100 atom %13C incorporation, two uniformly increasing point clouds with dissimilar slopes were produced (Figure 4). For both datasets, a fit on a linear model given by D R = b*P M with b equal to the slope value and an intercept of zero was carried out. The slope value b was estimated by adjusting the slope and minimizing the error sum of squares (SSE) of the decimal residual deviations:


with slope b calculated as:


with DR = decimal residual, PM = peptide mass, and n = number of peptide masses [31]. Slopes for 0 atom % and 100 atom %13C incorporation were estimated with an accuracy of 6 decimal digits and served as references for partial13C incorporation calculation.

All methods described above are summarized as easy-to-follow R-scripts which can be run under the statistical platform R [18]. 'Script 1' enables calculation of own reference dataset for 0% atom and 100% atom13C from any other protein mass dataset while reference slope calculation of can be done with help of 'Script 2'. Both scripts can be downloaded under [1].

Calculation of13C incorporation from practical mass spectrometric measurements

In order to estimate13C incorporation from experimentally derived datasets, processing occurred as described above. However, k-means clustering with very small datasets sometimes failed, because the algorithm required that at least one data point in each cluster to be present. In order to circumvent this difficulty, a set of dummy data containing one data point with values at the exact centers were added. Prior to13C incorporation calculation, those dummy points are removed again and, thus, have no influence on the following calculation.

In a second step, a linear fitting was conducted. While the standard linear fit by minimizing SSE is sufficient for large datasets, it is not feasible for smaller datasets like many user data consisting of relatively few peptide masses. Moreover, these data are often irregularly distributed and can contain outliers with extreme values characterized by strong deviations from linearity. In such cases, standard linear fitting is neither appropriate nor robust. Moreover, commonly applied SSE-based methods rely on the assumption that data residuals are normally distributed that, within user data is rarely found in practice. With datasets containing high variability and/or outliers, classical linear approximation methods often perform poorly [32, 33]. Therefore, we applied the more accurate and robust linear fitting algorithm (rlm) provided by the MASS package [30] in R that are better suited for small and highly variable datasets.

The rlm fitting approximation is using an iterative process examining which points cause strong shifts in the slope, and inversely re-weighting the data points during the least squares fitting process (IWLS-method [32, 34]). Using likelihood estimators, statistically robust methods minimize small deviations from model assumptions [30, 34].

Using reference slope values estimated (Figure 4) the relative13C incorporation of practical datasets is finally calculated by the relative fraction of the practical dataset


with13CUser % = relative isotope incorporation values from user data, buser = slope of regression estimated for practical data, b12C and b13C =12C and13C regression slope value for 0 atom %13C and 100 atom %13C incorporation, respectively.

We developed an easy-to-use R script for the estimation of13C atom % incorporation for user datasets which can be downloaded under [1] as 'script 3'.

User data have to be provided as an ASCII .txt file. The file should have the values of the peptide masses containing one peptide mass (z = 1 after deconvolution) per line with no header and/or row names. Measurements can be provided with comma or point decimal markers.

During the run of the script, a window opens and the user has to indicate where the ASCII.txt file is located. The file should have the values of the peptide masses containing one peptide mass (z = 1 after deconvolution) per line with no header and/or row names. After selecting the file, the calculation starts immediately and results in a scatter-plot with lines displaying the 0 atom % and 100 atom % incorporation reference lines, the user's peptide masses, the fitting line of the data with the standard error as quality proxy, and the estimated incorporation including the standard error estimation (Figure 5). In a loop process, the estimation of multiple measurement campaigns can be carried out.

Figure 5
figure 5

Typical plot resulting from user experimental measurements (dark circles), reference slope lines for unlabeled (0 atom %13C incorporation; dark grey line) and fully labeled (100 atom %13C incorporation; light grey line) and the slope calculated from user data measurements (thick black line). The13C incorporation for user data is directly calculated and displayed in the plot. This example consists of 150 practical peptides from Pseudomonas putida grown on fully labeled benzene.

Accuracy of the method

In order to determine the minimal number of needed peptide masses with an estimated incorporation accuracy of less than 5%, we additionally calculated an in silico peptide dataset with 50 atom % isotope incorporation of M. tuberculosis. We sub-sampled then the two already existing 0 and 100 atom, and the newly created 50%13C datasets by randomly drawing groups of 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 500, and 1,000 theoretical peptides across all peptide lengths with no repetition for 100 replications. For each within-group sub-sample, we calculated the median and the standard error of each subgroup. A resulting box-and-whisker plot for each drawing shows the accuracy of the incorporation as a function of number of peptide masses (Figure 6).

Figure 6
figure 6

Box-and-whisker plot for approximate accuracy of13C incorporation for groups of 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 500, and 1000 peptides from 100 random draws (no repetition) for Mycobacterium tuberculosis peptide dataset of 0, 50 and 100 atom %13C incorporation. Boxes represent lower and upper quartile, whiskers the upper and lower (95%) confidence interval and thick lines in boxes the median within the sub-sample groups. The two thin black lines depict 5% incorporation limits and the thick black line depicts the number of samples necessary to cross the 5% incorporation limit.

General applicability of the method

The general applicability of our method was tested using a dataset from Pseudomonas putida ML2 growing on fully labeled13C-benzene (0.6 mM) until the stationary growth phase. The detailed experimental conditions are explained by Jehmlich et al., [6]. After protein extraction, peptides were analyzed by nano-LC-LTQ Orbitrap-MS. From these dataset, information form 150 tryptically-digested peptide masses were taken and the heavy stable isotopes incorporation was estimated.

Results and Discussion

Data processing

For the calibration of the unlabeled and completely labeled incorporation, a dataset of M. tuberculosis containing atomic masses of > 90,000 peptide fragments were used. Plotting unlabeled peptide masses directly against their decimal residuals. As a result, up to m/z values of ~2,000 a first point cloud increased linearly with corresponding decimal residuals approaching values close to 1 (Figure 2A). When the decimal residuals reaching 1, the next residual digits start to fill up and thus start at 0 again. Therefore, increases across values of 1 leads to the development of the next point cloud with values close to 0. For m/z ranges between 300 and 5,000, new point clouds begin. Especially, for fully labeled peptide masses, these breaks occur at slightly lower values of ~1,800 Da, ~3,200 Da, and ~4,800 Da, since overall slopes of scatter plots here become slightly steeper (Figure 2B). However, variations in the molecular composition of each theoretical peptide mass impeded a simple separation into groups in order to create a continuous linear plot. Therefore, transient transformation using formula (1) was required. During the transformation process, peptide masses with small residual values (i.e. those close to 0) are comparatively altered very little, while those with big values (close to 1) experience strong shifts towards the ordinate and thus result in the expected increase of slope steepness for each group (Figures 3A and 3B). Subsequent to transformation with formula (1), allows straightforward separation of the point clouds for both the unlabeled and the fully labeled dataset, and the clustering algorithm was able to separate all data points into distinct classes. Misgrouping occasionally occurred with higher peptide masses, where points become more dispersed within the group. However, misgrouping occurred only in less than 0.1% of all data points.

After separation into groups, and the addition of 0, 1, 2 and 3 Da according to the group affiliation of each data point, two separated linear point clusters with distinctive slopes were obtained (Figure 4). Two characteristics can be observed: (i) the two point clouds display heteroscedastic characteristics. The initially compact plot becomes more dispersed at higher peptide masses, with higher inclination from the ideal linearity. However, peptide masses have increased variance; (ii) since the used database also contains short-chained amino acid sequences, the point clouds also become less dense with increasing peptide masses.

Linear slopes and estimation of relative isotope incorporation

The accurate slope value for the unlabeled dataset was estimated as b0% = 5.1357e-4 and for the fully labeled dataset b100% = 6.3347e-4. These two values were used in further calculations as our standards for the 0 atomic % and 100 atomic % incorporation.

By applying formula (4), we were able to calcuate incomplete heavy stable isotope incorporation into peptides of user data. However, standard linear fitting depending on minimizing SSE turned out to be strongly influenced by outliers. Strong deviation of the expected slope caused by outliers occurred most frequently in scenarios with inhomogeneous distribution of the measurements along the m/z range. In the presence of outliers, especially with high values, traditional methods are inefficient and biased because the least squares predictions are often heavily dragged towards these outliers. While every added point in the traditional SSE method has a direct influence on the slope, SSE has a breakdown point value of 0 and signifying that every added point has the same influences on the slope value. Our recommended robust linear model algorithm has a breakdown point value of 0.5, meaning that at least 50% of the data need to be altered to cause the slope estimation to change. Thus, when peptide masses with high mass strongly deviated from the linear reference slope, the final slope completely affected by these values. Using the robust linear fitting (rlm) algorithm, the slope for the calculation can only altered if more than 50% of all used data points varied. Therefore, the rlm algorithm is more robust and stable if outliers are used in the dataset.

Accuracy of method

The accuracy of incorporation is strongly depends on the considered number of peptide masses. As expected, the accuracy was very poor with small datasets (when using only 10-30 peptide masses). The prediction accuracy was only < 25% at known 0, 50, and 100 atom %13C incorporation (Figure 6). However, after including more than 50 peptides, the prediction asymptotically approached less than 10% precision. After incorporation of about 100 peptide masses, the accuracy became better than 5%. However, the inclusion of more peptide masses (> 200) did not improve the accuracy further.

General applicability of the method

The method can be applied to a broad range of bacterial taxa. Data from other bacteria than the M. tuberculosis reference strain proved also to be very successful. In order to validate the method, we used peptides of the bacterium P. putida ML2 grown on fully labeled substrate in a batch culture. Obtained mass spectra were analyzed by a LTQ-Orbitrap and the incorporation of13C into the peptides was calculated applying our method. After data points fitting a graph with a slope of 6.2910e-4 was generated that corresponds to 97.3% ± 0.3%13C incorporation. This value of13C incorporation is quite comparable with other MS-based methods where the labeling efficiency was estimated as 98.6% ± 0.2% for13C labeled peptides [35]. However, in comparison of our algorithm with other MS-based methods, generally for the latter a much higher effort for estimating incorporation values is needed.


We demonstrated that our method can be used to estimate the13C incorporation into peptides/proteins accurately and sensitively. The method requires about 100 peptide masses in order to achieve an accuracy of less than 5%.

Generally, our method holds promise for wider applications in qualitative and quantitative proteomics. Comprehensive proteomic analyses of mixed communities would enable microbial ecologist to link community composition, physiology, function, ecology, interaction, and evolutionary processes [36]. Metaproteome analyses were performed by either gel-based [37] or by LC-MS based approaches [38] and would go one-step further with the combination of metabolic labeling and proteome analysis (Protein-SIP) by obtaining the metabolic activity of various species in microbial communities.


  1. UFZ Leipzig. []

  2. Neufeld JD, Wagner M, Murrell JC: Who eats what where and when? Isotope-labelling experiments are coming of age. Isme J. 2007, 1 (2): 103-110. 10.1038/ismej.2007.30.

    Article  CAS  PubMed  Google Scholar 

  3. Boschker HTS, Nold SC, Wellsbury P, Bos D, de Graaf W, Pel R, Parkes RJ, Cappenberg TE: Direct linking of microbial populations to specific biogeochemical processes by C-13-labelling of biomarkers. Nature. 1998, 392 (6678): 801-805. 10.1038/33900.

    Article  CAS  Google Scholar 

  4. Radajewski S, McDonald IR, Murrell JC: Stable-isotope probing of nucleic acids: a window to the function of uncultured microorganisms. Curr Opin Biotechnol. 2003, 14 (3): 296-302. 10.1016/S0958-1669(03)00064-8.

    Article  CAS  PubMed  Google Scholar 

  5. Manefield M, Whiteley AS, Griffiths RI, Bailey MJ: RNA stable isotope probing a novel means of linking microbial community function to phylogeny. Appl Environ Microbiol. 2002, 68 (11): 5367-5373. 10.1128/AEM.68.11.5367-5373.2002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Jehmlich N, Schmidt F, Hartwich M, von Bergen M, Richnow HH, Vogt C: Incorporation of carbon and nitrogen atoms into proteins measured by protein-based stable isotope probing (Protein-SIP). Rapid Commun Mass Spectrom. 2008, 22 (18): 2889-2897. 10.1002/rcm.3684.

    Article  CAS  PubMed  Google Scholar 

  7. Jehmlich N, Schmidt F, Taubert M, Seifert J, von Bergen M, Richnow HH, Vogt C: Comparison of methods for simultaneous identification of bacterial species and determination of metabolic activity by protein-based stable isotope probing (Protein-SIP) experiments. Rapid Commun Mass Spectrom. 2009, 23 (12): 1871-1878. 10.1002/rcm.4084.

    Article  CAS  PubMed  Google Scholar 

  8. Cargile BJ, Bundy JL, Grunden AM, Stephenson JL: Synthesis/degradation ratio mass spectrometry for measuring relative dynamic protein turnover. Analytical Chemistry. 2004, 76 (1): 86-97. 10.1021/ac034841a.

    Article  CAS  PubMed  Google Scholar 

  9. Papageorgopoulos C, Caldwell K, Shackleton C, Schweingrubber H, Hellerstein MK: Measuring protein synthesis by mass isotopomer distribution analysis (MIDA). Analytical Biochemistry. 1999, 267 (1): 1-16. 10.1006/abio.1998.2958.

    Article  CAS  PubMed  Google Scholar 

  10. McIlwain S, Page D, Huttlin EL, Sussman MR: Using dynamic programming to create isotopic distribution maps from mass spectra. Bioinformatics. 2007, 23 (13): i328-336. 10.1093/bioinformatics/btm198.

    Article  CAS  PubMed  Google Scholar 

  11. Choudhary K, Spicer VL, Donald LJ, Duckworth HW, Ens W, Loewen PC, Standing KG: Method for estimating the isotopic distributions of metabolically labeled proteins by MALDI-TOFMS: Application to NMR samples. Analytical Chemistry. 2006, 78 (15): 5419-5423. 10.1021/ac060507d.

    Article  CAS  PubMed  Google Scholar 

  12. Snijders APL, de Koning B, Wright PC: Perturbation and interpretation of nitrogen isotope distribution patterns in proteomics. Journal of Proteome Research. 2005, 4 (6): 2185-2191. 10.1021/pr050260l.

    Article  CAS  PubMed  Google Scholar 

  13. Jacob U, Brey T, Fetzer I, Kaehler S, Mintenbeck K, Dunton K, Beyer K, Struck U, Pakhomov EA, Arntz WE: Towards the trophic structure of the Bouvet Island marine ecosystem. Polar Biology. 2006, 29 (2): 106-113. 10.1007/s00300-005-0071-8.

    Article  Google Scholar 

  14. Jehmlich N, Schmidt F, von Bergen M, Richnow HH, Vogt C: Protein-based stable isotope probing (Protein-SIP) reveals active species within anoxic mixed cultures. Isme J. 2008, 2 (11): 1122-1133. 10.1038/ismej.2008.64.

    Article  CAS  PubMed  Google Scholar 

  15. Brun V, Dupuis A, Adrait A, Marcellin M, Thomas D, Court M, Vandenesch F, Garin J: Isotope-labeled protein standards: toward absolute quantitative proteomics. Mol Cell Proteomics. 2007, 6 (12): 2139-2149. 10.1074/mcp.M700163-MCP200.

    Article  CAS  PubMed  Google Scholar 

  16. Domon B, Aebersold R: Mass spectrometry and protein analysis. Science. 2006, 312 (5771): 212-217. 10.1126/science.1124619.

    Article  CAS  PubMed  Google Scholar 

  17. Schmidt F, Schmid M, Jungblut PR, Mattow J, Facius A, Pleissner KP: Iterative data analysis is the key for exhaustive analysis of peptide mass fingerprints from proteins separated by two-dimensional electrophoresis. J Am Soc Mass Spectrom. 2003, 14 (9): 943-956. 10.1016/S1044-0305(03)00345-3.

    Article  CAS  PubMed  Google Scholar 

  18. R project. []

  19. Gay S, Binz PA, Hochstrasser DF, Appel RD: Modeling peptide mass fingerprinting data using the atomic composition of peptides. Electrophoresis. 1999, 20 (18): 3527-3534. 10.1002/(SICI)1522-2683(19991201)20:18<3527::AID-ELPS3527>3.0.CO;2-9.

    Article  CAS  PubMed  Google Scholar 

  20. Mann M: 43rd ASMS Conference on Mass Spectrometry and Applied Topics: 1995; Atlanta GA. 1995: 639-640.

  21. Karty JA, Ireland MM, Brun YV, Reilly JP: Artifacts and unassigned masses encountered in peptide mass mapping. J Chromatogr B Analyt Technol Biomed Life Sci. 2002, 782 (1-2): 363-383. 10.1016/S1570-0232(02)00550-0.

    Article  CAS  PubMed  Google Scholar 

  22. Mattow J, Schmidt F, Hohenwarter W, Siejak F, Schaible UE, Kaufmann SH: Protein identification and tracking in two-dimensional electrophoretic gels by minimal protein identifiers. Proteomics. 2004, 4 (10): 2927-2941. 10.1002/pmic.200400908.

    Article  CAS  PubMed  Google Scholar 

  23. Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, Gordon SV, Eiglmeier K, Gas S, Barry CE: Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998, 393 (6685): 537-544. 10.1038/31159. 3

    Article  Google Scholar 

  24. Sanger Institute. []

  25. ProteinProspector - MS Digest. []

  26. Baxter CA, Murray CW, Clark DE, Westhead DR, Eldridge MD: Flexible docking using Tabu search and an empirical estimate of binding affinity. Proteins. 1998, 33 (3): 367-382. 10.1002/(SICI)1097-0134(19981115)33:3<367::AID-PROT6>3.0.CO;2-W.

    Article  CAS  PubMed  Google Scholar 

  27. Eldridge MD, Murray CW, Auton TR, Paolini GV, Mee RP: Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes. J Comput Aided Mol Des. 1997, 11 (5): 425-445. 10.1023/A:1007996124545.

    Article  CAS  PubMed  Google Scholar 

  28. Parker KC: Scoring methods in MALDI peptide mass fingerprinting: ChemScore and the ChemApplex program. J Am Soc Mass Spectrom. 2002, 13 (1): 22-39. 10.1016/S1044-0305(01)00320-8.

    Article  CAS  PubMed  Google Scholar 

  29. Hartigan JA, Wong MA: A k-means clustering algorithm. Applied Statistics. 1979, 28 (1): 100-108. 10.2307/2346830.

    Article  Google Scholar 

  30. Venables W, Ripley B, (eds.): Modern Applied Statistics with S. 2002, New York: Springer

  31. Crawley MJ, (ed.): Statistical Computing - An Introduction to Data Analysis Using S-Plus. 2002, Wiley-VCH

  32. Hampel FR, Ronchetti EM, Rousseeuw PJ, Stahel WA, (eds.): Robust Statistics - The Approach Based on Influence Functions. 2005, Wiley-VCH

  33. Rousseeuw PJ, Leroy AM, (eds.): Robust Regression and Outlier Detection. 2003, Wiley-VCH

  34. Huber PJ, (eds): Robust Statistics. 1981, New York: Wiley-VCH

  35. Snijders APL, de Vos MGJ, Wright PC: Novel approach for peptide quantitation and sequencing based on N-15 and C-13 metabolic labeling. Journal of Proteome Research. 2005, 4 (2): 578-585. 10.1021/pr0497733.

    Article  CAS  PubMed  Google Scholar 

  36. Wilmes P, Bond PL: Microbial community proteomics: elucidating the catalysts and metabolic mechanisms that drive the Earth's biogeochemical cycles. Curr Opin Microbiol. 2009, 12 (3): 310-317. 10.1016/j.mib.2009.03.004.

    Article  CAS  PubMed  Google Scholar 

  37. Benndorf D, Balcke GU, Harms H, von Bergen M: Functional metaproteome analysis of protein extracts from contaminated soil and groundwater. Isme Journal. 2007, 1 (3): 224-234. 10.1038/ismej.2007.39.

    Article  CAS  PubMed  Google Scholar 

  38. Ram RJ, VerBerkmoes NC, Thelen MP, Tyson GW, Baker BJ, Blake RC, Shah M, Hettich RL, Banfield JF: Community proteomics of a natural microbial biofilm. Science. 2005, 308 (5730): 1915-1920. 10.1126/science. 1109070.

    Article  CAS  PubMed  Google Scholar 

Download references


We would like to thank Christine Schumann and Michaela Öhler for excellent technical assistance. We gratefully want to thank Brandon E. Morris for kindly revising the English of this manuscript. Part of the work was funded by a grant within the priority programme SPP1319 provided by the German Research Foundation (DFG).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ingo Fetzer.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

IF participated in the conceptual design of study, conducted the computational implementation of the algorithm, conducted all calculations and statistical analyses. He also implemented the method into 'R' and wrote the scripts for user application of the method. Together with NJ, he drafted the manuscript. NJ contributed to the acquisition and interpretation of protein data and isotopic measurements used for manuscript and protein data analysis. CV conducted the necessary experimental work in the lab, carried out the lab data acquisition and final interpretation of the experimental data. He also critically revised the manuscript and helped to improve it substantially. H-HR contributed in the isotopic interpretation of the data, participated in the critical revision of the manuscript and helped in the final improvement of the version for publishing. JS performed protein data analysis and did the protein data acquisition and interpretation. HH revised the MS critically and helped to restructure the MS to its final publishable version. MvB was involved in the design of the study, mainly contributed to the introductory section of the manuscript, revised the manuscript critically and participated in the final improvement of this version. FS came up with the original idea for the algorithm, participated in the design of the study and was responsible for the coordination of the study. Additionally, he substantially contributed to the writing and improvement of the manuscript.

All authors read and approved the revised final mansucript

Ingo Fetzer, Nico Jehmlich contributed equally to this work.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fetzer, I., Jehmlich, N., Vogt, C. et al. Calculation of partial isotope incorporation into peptides measured by mass spectrometry. BMC Res Notes 3, 178 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: