Molecular Isotopic Distribution Analysis (MIDAs) with adjustable mass accuracy.

In this paper, we present Molecular Isotopic Distribution Analysis (MIDAs), a new software tool designed to compute molecular isotopic distributions with adjustable accuracies. MIDAs offers two algorithms, one polynomial-based and one Fourier-transform-based, both of which compute molecular isotopic distributions accurately and efficiently. The polynomial-based algorithm contains few novel aspects, whereas the Fourier-transform-based algorithm consists mainly of improvements to other existing Fourier-transform-based algorithms. We have benchmarked the performance of the two algorithms implemented in MIDAs with that of eight software packages (BRAIN, Emass, Mercury, Mercury5, NeutronCluster, Qmass, JFC, IC) using a consensus set of benchmark molecules. Under the proposed evaluation criteria, MIDAs's algorithms, JFC, and Emass compute with comparable accuracy the coarse-grained (low-resolution) isotopic distributions and are more accurate than the other software packages. For fine-grained isotopic distributions, we compared IC, MIDAs's polynomial algorithm, and MIDAs's Fourier transform algorithm. Among the three, IC and MIDAs's polynomial algorithm compute isotopic distributions that better resemble their corresponding exact fine-grained (high-resolution) isotopic distributions. MIDAs can be accessed freely through a user-friendly web-interface at http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/midas/index.html.

Introduction

Most biomolecules are composed of hydrogen, carbon, nitrogen, oxygen, and sulphur. It is known that the natural isotopes of these elements occur with different probabilities [1, 2], and in some experiments the relative abundances of an element’s isotopes can be manipulated by using a technique known as stable isotopic labeling [3, 4]. The relative abundances of isotopes determine a molecule’s isotopic distribution (ID), which can be measured experimentally using a mass spectrometer. The measured ID constrains the elemental composition when compared with the in-silico computed ID and, hence, helps in identifying the underlying molecule. The realization of this goal, however, demands accurate in-silico ID prediction [5-11].

The information content in an experimentally measured ID depends on the resolution of the mass spectrometer. An ID generated by a low resolution instrument contains less information than that by an ultra-high resolution instrument [12-15]. Based on the instrument resolution, three different types of IDs are commonly mentioned in the literature: the aggregated, the fine structure, and the hyper-fine structure IDs [16]. The aggregated ID is computed by merging isotopic variants that have the same nucleon number into one aggregated isotopic variant [17, 18] whose corresponding molecular mass (MM) and occurrence probability are computed respectively from the probability-weighted sum of masses and from the sum of the probabilities of the isotopic variants merged. The fine and hyper-fine structure IDs are computed similarly to the aggregated ID, except that one merges only isotopic variants whose molecular mass differences are within some pre-specified mass accuracy.

To make practical use of experimentally measured IDs, it is imperative to have methods that can compute in-silico IDs when given molecular formulas. Rockwood et al. [19] mentioned several criteria for a sound ID-computing method (IDCM): an IDCM must accurately compute in a very short time the masses and intensities without consuming much computational resource. We propose a few additional criteria by which to assess an IDCM’s application value: to handle experimentally generated IDs from both low-resolution and high-resolution instruments, an IDCM should allow adjustable mass accuracy; given that customized isotopic labeling has become a common experimental technique for quantitative analyses, an IDCM should be able to handle customized (or user-specified) isotopic abundances (or occurrence probabilities) of all chemical elements considered; finally, an IDCM should be able to compute IDs for a wide mass range and be user-friendly. Although there are several available methods [16] that can compute an aggregated ID [17, 18, 20–23], fine structure ID [19], and hyper-fine structure ID [23-26], there are not many methods that can satisfy all the requirements mentioned above.

In this manuscript, we present MIDAs, a software tool satisfying all the requirements above. MIDAs provides users with two accurate and efficient algorithms to compute IDs: the first algorithms belongs to the class of polynomial methods [27, 28], whereas the other algorithm belongs to the class of Fourier transform methods [29, 30]. The latter consists mainly of changes made to the existing Fourier transform method [19], and the changes made are shown to improve significantly the accuracy of the computed ID. Both algorithms can compute low and high resolution IDs, referred to as the coarse-grained isotopic distribution (CGID) and the fine-grained isotopic distribution (FGID), respectively, for the remainder of this manuscript. Also both algorithms implemented in MIDAs are capable of computing CGID and FGID with adjustable mass accuracy.

To evaluate the performance of MIDAs, we have benchmarked it against eight methods: four of these methods—Mercury [19], NeutronCluster (NC) [17], Emass [21], and BRAIN [18, 31]—are the four best performing methods taking from a recent publication by Claesen et al. [18]; four other methods included are Mercury5 (a new version of Mercury2) [32], Qmass [20], Isotope Calculator (IC) [33], and a Fourier-transform-based method recently published [34], which we refer to as JFC. JFC is an improved version of Isotopica [35], which incorporates BRAIN’s generating function. The program of JFC was downloaded from http://bioinformatica.cigb.edu.cu/isotopica/centermass.html. The BRAIN code was downloaded from http://www.bioconductor.org/packages/release/bioc/html/BRAIN.html. The program IC was downloaded from http://agarlabs.com/. The rest of the programs were provided by the code authors, whom we acknowledge in the Acknowledgment section.

The performance evaluation was conducted using 25 molecules. Ten of these molecules are benchmark proteins previously used to evaluate the accuracy of computed CGIDs [17, 18]. Another 10 are hydrocarbon molecules whose CGIDs and FGIDs can be exactly computed, making them ideal for evaluating the accuracy of computed IDs. The remaining five molecules, made of a combination of sulfur, mercury, carbon and hydrogen, are used together with some of the other 20 molecules to evaluate the computational time of MIDAs’s algorithms. Results from our investigation show that MIDAs [both the polynomial-based algorithm (MIDAs) and the Fourier-transform-based algorithm (MIDAs)], Emass, and JFC compute CGIDs with equivalent accuracy and are more accurate than the other methods evaluated. When computing the FGIDs, IC and MIDAs yield FGIDs that are closest to the exact FGIDs. The results also show that MIDAs and MIDAs satisfy all aforementioned requirements to be considered a valuable tool, providing the community with two new options for computing accurate IDs.

Methods

In the subsections below we explain in detail the two algorithms implemented in MIDAs. The first subsection explains MIDAs, a polynomial-based algorithm. The second subsection describes MIDAs, a fast Fourier transform (FFT) based algorithm. Both algorithms can be used to compute CGIDs and FGIDs.

MIDAs Polynomial Multiplication Algorithm (MIDAs)

It is well known that the ID of a molecule can be obtained by expanding the corresponding product of polynomials: each expanded term corresponds to an isotopic composition of the molecule’s elements. For example, the ID of a molecule having molecular formula (MF) x y is given by expandingwhere I is an indicator variable, x and y are the isotopes of elements x and y, respectively, p(x ) and p(y ) are normalized probabilities of occurrence, and m(x ) and m(y ) are the exact atomic masses.

There are several polynomial-based methods designed to compute an ID from the MF. Methods such as the stepwise procedure and its improvement [36, 37], symbolic expansion [4], and multinomial expansion [28, 38] have been proposed to compute the expansion of the above polynomial. Although these methods have been shown to perform well for small molecules, they fail to handle large molecules, yielding inaccurate IDs, requiring a significant amount of computer memory, and taking a considerable amount of computational time [16].

Here we present MIDAs, a polynomial-based algorithm that is simple and easy to understand. Our algorithm computes the molecule’s CGID by directly performing polynomial multiplication. To simplify the explanation, define the polynomial between the brackets in Equation (1) containing the probabilities and atomic masses of an element’s isotopes as the element fundamental polynomial (EFP). Let us represent the EFP of element x by P , and also define the following recursion operation that multiplies together polynomials Q and P and assigns the resulting polynomial back to Q as Q  ← (Q  × P ). Substituting these definitions in Equation (1) with Q initialized to one giveswhere ⌊z⌋ represents the integer part of z for any positive number z. Using the recursion operation mentioned earlier, all the x-element related polynomials finally merge into Q and all the y-element related polynomials finally merge into Q as shown in algorithms 1 and 2. By first computing in Equation (2), one considerably reduces the computational time needed to obtain the polynomial expansion of an EFP. The logic in computing (or ) and not [P ] (or [P y]) is that the former requires a smaller number of arithmetic operations. This is due to two heuristic procedures of MIDAS, prune and merge, which reduce the number of retained terms in the expanded polynomial . These heuristics are similar to the hyperatom concept [37] and the superatom concept [10], except that the number of atoms ⌊N/10⌋ in a superstructure is not fixed in our case. The choice of using 10 in Equation (2) was somewhat arbitrary but seemed to generate an accurate ID for each molecule used in our investigation in a short amount of computational time. Evidently, one may use a number other than 10. Choosing a smaller number, however, means that we need a larger memory to hold Q. Choosing a larger number, on the other hand, results in a longer computation time. We find using the number 10 seems to provide a good balance between the two.

The first heuristic employed by the MIDAs algorithm prunes terms from the polynomial Q that have probability smaller than a pre-set probability value (η). The second heuristic procedure merges polynomial terms from Q that are within some user specified mass accuracy (∈) of each other into a new polynomial term. The new polynomial term is assigned a new mass () that is equal to the probability-weighted sum of the merged termswhere m and p stand for the mass and probability of the merged terms, respectively. This new term associated with is then assigned a probability equal to the sum of the probabilities of the merged terms. The pseudo-code for computing a CGID is given by algorithm 1, which is used by MIDAs.

Algorithm 1. Computes Coarse-Grained Isotopic Distribution

To compute the FGID for an MF, for every element x, MIDAs first computes the expected number of occurrences μ[x ] for each isotope x of x. MIDAs then computes σ 2[x ], the variance of the number of occurrences. As an example, for the molecular formula x y , the expectation and variance in the number of atoms for a given isotope of element x is given byand

Using the computed expectation and variance values, we denote the range as allowable for , the number of atoms of isotope x . The upper bound and the lower bound are given by

For isotope x , we choose to be the span of sum as this quantity is guaranteed to be simultaneously larger then 10σ [x ] and 10 daltons. For each element x, the and are used to construct a polynomial, , by means of the multinomial expansion formula

By summing only the contributions bounded by and , we direct the calculations to the relevant part of the ID. It has counter-part in FT based method, namely the heterodyning of in [24].

For numerical accuracy and efficiency we employ the following simple identityand ln(n !) = ∑  ln k. This representation reduces computational time of Equation (7) since by tabulation one only enumerates all the logarithmic terms in Equation (8) once. Once all have been computed, they are used together with a user specified ∈ to compute a FGID using algorithm 2.

Algorithm 2. Computes Fine-Grained Isotopic Distribution 2

MIDAs Fast Fourier Transform Algorithm (MIDAs)

The MIDAs algorithm is similar to an early FFT algorithm by Rockwood et al. [19], which was implemented in a computer program called Mercury. These two algorithms differ, however, in a few aspects. First, using the exact isotopic masses in discrete FFT (DFFT) [39, 40], Mercury produces IDs with leakages (assigning nonzero probabilities to masses where exactly zero probability is expected) and employs an apodization function to minimize leakage [41]. On the other hand, by assigning each isotope mass to a point on a fixed grid, MIDAs avoids the leakage problem. Using discrete masses to avoid leakage is not new: Rockwood and Van Orden [32] have written a computer program, whose latest version is called Mercury5, to compute IDs based on the nucleon numbers (or roughly using one dalton mass grid). The improvement we made was to allow the users to specify the mass accuracy other than 1 Da. Second, Mercury uses a fixed number of sample points with the DFFT, whereas in MIDAs the number of sample points used depends on the mass accuracy, which is a parameter adjustable by the user.

Every FFT based method relies on the convolution theorem, which states that a convolution can be performed as multiplication in the Fourier domain.

As we shall discuss in the Appendix, there are two key conditions in order for the convolution theorem to be used in the discrete case while computing IDs. The first one is that the masses of each isotope must lie on grid points. Using a mass that is not on the grid causes the “leakage" phenomenon [41]. If the masses considered all reside on grid points, the leakage problem no longer exists. The second important condition is that the mass domain must be large enough so that the “folded-back" phenomenon (which is also known as “aliasing”, “fold over”, or “wrap around” in the signal processing community) near the tail of the distribution is negligible (see Appendix).

Prior to delving into detail constructs of MIDAs, let us first describe how one may compute the theoretical molecular mass variance σ 2. Using our example molecule x y , one note that the molecular mass variance of this molecule can be rigorously written as σ 2 = Nσ 2[m ] + Mσ 2[m ], where σ 2[m ] is the molecular mass variance associated with element x. Explicitly, one may calculate σ 2[m ] as followswhere the index i runs over all isotopes of element x and p(x ) again represents the occurrence probability of isotope x .

A key constraint of DFFT based ID method is that the total number of sample points, denoted by S, must be an integral power of two [42]. For a given molecule and specified mass accuracy ∈, the total number of sample points S used in MIDAs’s DFFT is given byand σ 2 is the theoretical variance in MM due to the elements’ isotopes [32]. The quantity ⌈z⌉ represents the smallest integer that is larger than z for any positive number z. Again the quantity is chosen so that S covers on both ends more than 7.5 standard deviations from the mean molecular mass, which prevents folded-back mass regions from having significant probabilities.

In order to avoid the problem of leakage, instead of using exact masses of isotopes and then applying filtering windows, we pin all isotopic masses to grid points. For each isotope mass m(x ), we first find a corresponding grid index n(x ) by the following formula

Using this discrete approach, the probability function of the mass of element x becomeswhere n(x )∈ is the approximate expression for the exact mass m(x ), and the Kronecker delta function takes value one if its two indices coincide and zero otherwise.

Consider the mass distribution of our example molecule x y . By the convolution theorem, the Fourier transform of the mass distribution, denoted by Ψ(v), can be written aswhere v takes S discrete values: 0,1,…, S − 1. The sample function Ψ(v) is heterodyned to have zero average mass by multiplying it by , where n is equal to (the molecule’s probability-weighted average grid index computed using n(x ) and p(x )) rounded to the nearest integer.

Once the function Ψ(v) has been calculated, three other operations are performed in order to generate the final FGID and CGID. The first operation performed is the inverse discrete fast Fourier transform (IDFFT), which transforms the sample function Ψ(v) to Φ(n) on the mass grid. Second, we apply a denoising procedure to remove small amplitudes due to rounding errors that occur during IDFFT. The rounding errors are expected to create small positive and negative amplitudes of equal amounts in the mass domain. MIDAs thus removes all amplitudes whose absolute magnitude are smaller than that of the most negative amplitude. As a matter of fact, to be more conservative, MIDAs uses an amplitude cutoff value that is twice the absolute value of the most negative amplitude. This means that only terms having amplitude greater than the cutoff value are reported in a computed FGID and CGID, with the amplitude values renormalized to sum to one. Figure 1 shows an example of the overlap between the positive amplitude histogram and the negative amplitude histogram. Right below the cutoff absolute amplitude, we see that the two histograms resemble each other, reflecting the fact that rounding errors have equal probability to be positive and negative. Following Rockwood and Van Orden [32], in the third step, MIDAs applies a linear transformation to rescale the masses associated with the IDs to ensure a good agreement between the theoretically calculated and the numerically computed mean molecular mass as well as standard deviation of the molecular mass. The procedure described above is summarized in the pseudo code give by algorithm 3.

Figure 1

Example of rounding errors. The curves plotted above are the histogram for the logarithm of absolute value of the positive (green solid line histogram) and negative (blue long-dashed line histogram) amplitudes obtained after applying the discrete Fourier transform to compute an isotopic distribution for molecule C2023H3208N524O619S20 using a mass accuracy of 0.01 Da (≈0.22 ppm). Absent the leakage, the negative amplitude can only come from rounding errors, among which equal amounts of small positive amplitudes and negative amplitudes are expected. The above histograms overlap for terms that have magnitude in amplitude less than 4.2e–10, displayed above by a dash black line, and at this point is the rounding error cutoff value used by MIDAs

Algorithm 3. Computes Fine-Grained and Coarse-Grained Isotopic Distribution

Results and Discussion

All methods used in our investigation were evaluated using their default parameter settings, except for a few parameter changes made to ensure that the atomic masses and abundances of elements’ isotopes were the same for all methods (see Table 1). To conduct the evaluation, we used the 22 biomolecules and three inorganic compounds, all of which are listed in Table 2. Ten of these biomolecules, (1)–(10), are benchmark proteins previously used to evaluate the computed CGID [17, 18]; 10 biomolecules, (11)–(20), are hydrocarbon molecules whose FGID can be exactly computed and were employed to evaluate the FGID; the remaining five molecules, (21)–(25), made of a combination of sulfur, mercury, carbon, and hydrogen, are used together with the some of the other 20 biomolecules to evaluate the computational time of MIDAs’s algorithms.

Table 1

Atomic Masses and Abundances used for Benchmark Test in this Paper

Isotope	Atomic mass Da	Abundance (%)
Atomic masses and naturally occurring isotopic abundances [1]
12C	12.0000000000	98.9300
13C	13.0033548378	1.0700
1H	1.0078250321	99.9885
2H	2.0141017780	0.0115
14 N	14.0030740052	99.6320
15 N	15.0001088984	0.3680
160	15.9949146	99.7570
170	16.9991312	0.0380
180	17.9991603	0.2050
32S	31.97207070	94.9300
33S	32.97145843	0.7600
34S	33.96786665	4.2900
36S	35.96708062	0.0200
196Hg	195.965833	0.0015
198Hg	197.966769	0.0997
199Hg	198.968279	0.1687
200Hg	199.968326	0.2310
201Hg	200.970302	0.1318
202Hg	201.970643	0.2986
204Hg	203.973493	0.0687
Atomic Masses and Enriched Carbon’s Isotopic Abundances
12C	12.0000000000	1.0000
13C	13.0033548378	99.0000

Table 2

Molecules for which the Isotopic Distribution was Computed by Various Methods

No. 1	Molecular formula	Lightest Mass (Da)2	Average Mass (Da)
(1)	C50H71N13O12	1045.5345145467	1046.1811074558
(2)	C254H377N65O75S6	5729.6008666397	5733.5107592120
(3)	C520H817N139O147S8	11616.8493497485	11624.4487510271
(4)	C744H1224N210O222S5	16812.9547750824	16823.3213522608
(5)	C2023H3208N524O619S20	45387.0070331016	45415.6793695079
(6)	C2934H4615N781O897S39	66389.8624747027	66432.4555603617
(7)	C5047H8014N1338O1495S48	112823.8795468070	112895.1259319964
(8)	C8574H13378N2092O2392S77	186386.7992654122	186506.0525933526
(9)	C17600H2674N4752O5486S197	398470.3669960258	398722.9724824960
(10)	C23832H37816N6528O7031S170	533403.4750914392	533735.2146493989
(11)	C5H5	65.0391251605	65.0933832534
(12)	C10H10	130.0782503209	130.1867665069
(13)	C50H50	650.3912516049	650.9338325345
(14)	C100H100	1300.7825032099	1301.8676650690
(15)	C1000H1000	13007.8250320999	13018.6766506902
(16)	C10000H10000	130078.2503209999	130186.7665069023
(17)	C20000H20000	260156.5006419999	260373.5330138047
(18)	C30000H30000	390234.7509629999	390560.2995207072
(19)	C40000H40000	520313.0012839999	520747.0660276095
(20)	C50000H50000	650391.2516049999	650933.8325345119
(21)	S20000	639441.4139999999	641321.6938997399
(22)	Hg5000	159860.3534999999	160330.4234749349
(23)	Hg1000S1000	227937.9037000000	232665.2510595869
(24)	S1000C1000H1000	44979.8957320999	45084.7613456772
(25)	Hg1000C1000H1000	208973.6580321000	213617.8430152902

1Reference number associated with a molecular formula (biomolecule or inorganic compound).

2The unified atomic mass unit dalton (Da).

Overview of Methods Benchmarked

MIDAs’s performance was evaluated against eight published methods: Mercury [19], Mercury5 [32], JFC [34], Isotope Calculator (IC) [33], Qmass [20], BRAIN [18, 31, 43], NeutronCluster (NC) [17], and Emass [21]. The first three published methods are Fourier-transform-based, IC utilizes a divide-and-recursively-combine algorithm, Qmass has its core based on FFT, BRAIN and NeutronCluster are polynomial-based, whereas Emass is based on a direct convolution approach related to the stepwise procedure and its improvement [36, 37]. BRAIN, Qmass, NC, Emass, JFC, and Mercury5 all use nucleon numbers to classify molecule’s isotopic variants, while all but the last assign to a given nucleon number the average isotopic mass of all variants of that nucleon number.

IC is suitable for computing FGIDs, not CGIDs. Qmass, BRAIN, NeutronCluster, and Emass are suitable for computing CGIDs, not FGIDs. The remaining three Fourier-transform-based methods are also suitable for computing CGIDs, although Mercury is the only one that has FGID computing capacity. To benchmark the FGIDs computed by MIDAs against those of Mercury, however, would require post-processing of Mercury data files such as removing noise from leakage and rounding errors, as well as compiling output from different specified molecular masses. All of these steps may be done differently and make the benchmark test less meaningful. For these reasons, we only evaluated MIDAs’s FGIDs against that of IC, not that of Mercury.

Benchmarking of Computed CGIDs

Following previous publications [18, 19, 24], the accuracy of a method is gauged by how accurately it yields ID mean, ID standard deviation, lightest and heaviest molecular masses, while computing a CGID. In our evaluation, the lightest mass and heaviest molecular mass are defined as a molecule’s molecular mass computed using the masses of the lightest and heaviest isotopes, respectively.

Lightest masses comparisons for biomolecules, numbered (1)–(10) in Table 2, with elements having their naturally occurring isotopic abundances taken from Table 1, are displayed in Table 3. Unexpectedly, the lightest masses for the first six molecules, reported by Mercury, Mercury5, and Qmass, are even lighter than their exact lightest masses, which should be the lightest masses possible in these six IDs. This observation was also described by [18]. For Mercury5 (and Mercury), this is caused by the rounding errors (and the leakage) when applying the DFFT. (In principle, both methods can avoid this problem by not reporting any terms in the computed ID that have masses lighter than the exact lightest mass.) For Qmass, this seems to arise from computing ID terms that are outside of the allowed mass range imposed by the biomolecule’s MF. This is because in the Qmass output file the reported masses lighter than the exact lightest mass are associated with elemental compositions that differ from the biomolecule’s MF used in the evaluation.

Table 3

Coarse - Grained Isotopic Distribution Results using Naturally Occurring Isotopes

Difference in lightest mass
No.1	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	0	-3.4e - 05	-2.6e - 10	2.2e - 13	7.0	7.3	0	12.2	0
(2)	0	-1.7e - 03	-1.3e - 09	0	12.0	12.1	0	15.9	0
(3)	0	-2.6e - 03	-2.8e - 09	0	8.0	8.4	0	18.2	-1.0e - 10
(4)	0	-2.1e - 03	-4.2e - 09	0	22.0	21.6	-360	39.1	8.0e - 10
(5)	7.2e - 12	-7.4e - 03	1.0e - 08	1.4e - 11	2.9	3.3	0	0.045	-1.6e - 01
(6)	-1.4e - 11	-5.0	-1.6e - 08	0	22.0	21.3	0	65.2	-4.6
(7)	1.5e - 11	-19.1	-2.7e - 08	-8.0	-8.0	-7.2	0	-69.7	-18.1
(8)	0	-49.1	-4.1e - 08	-31.1	-55.2	-55.3	0	-90.7	-48.8
(9)	-5.8e - 11	-147.4	-1.1e - 07	-114.3	-124.3	-124.7	0	-188.3	-118.9
(10)	0	-210.6	-1.2e - 07	-172.5	-203.7	-203.7	0	-355.6	-146.4
Difference in heaviest isotopic mass
No.	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	7.4e + 01	1.4e + 02	1.5e + 02	1.4e + 02	1.5e + 02	1.5e + 02	1.5e + 02	1.4e + 02	1.4e + 02
(2)	7.2e + 02	8.4e + 02	8.7e + 02	8.4e + 02	8.5e + 02	8.5e + 02	8.6e + 02	8.4e + 02	8.4e + 02
(3)	1.6e + 03	1.8e + 03	1.8e + 03	1.7e + 02	1.8e + 03	1.8e + 03	1.8e + 03	1.8e + 03	1.8e + 03
(4)	2.5e + 03	2.6e + 03	2.6e + 03	2.6e + 03	2.6e + 03	2.6e + 03	2.3e + 03	2.6e + 03	2.6e + 03
(5)	6.8e + 03	7.0e + 03	7.0e + 03	7.0e + 03	7.0e + 03	7.0e + 04	7.0e + 03	6.9e + 03	7.0e + 03
(6)	9.9e + 03	1.0e + 04	1.0e + 04	1.0e + 04	1.0e + 04	1.0e + 04	1.0e + 04	1.0e + 04	1.0e + 04
(7)	1.7e + 04	1.7e + 04	1.7e + 04	1.7e + 04	1.7e + 04	1.7e + 04	1.7e + 04	1.8e + 04	1.7e + 04
(8)	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04	2.9e + 04
(9)	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04	6.0e + 04
(10)	8.2e + 04	8.3e + 04	8.3e + 04	8.3e + 04	8.3e + 04	8.3e + 04	8.3e + 04	8.2e + 04	8.2e + 04

1Reference number associated with a molecular formula (biomolecule).

The software NC reports correct lightest masses for nine out of the 10 molecules. For biomolecule number four, NC reports a mass that is 360 Da heavier. This same result has also been observed independently by others [17, 18].

For MIDAs, BRAIN, and Emass, the differences between exact and computed lightest masses, for small and medium size biomolecules [numbered (1)–(6)], are smaller than 1.0e–08 Da. As for JFC and MIDAs, although they do not perform as well as the polynomial-based methods above, they are not inferior to other Fourier-transform-based methods such as Mercury and Mercury5. When the biomolecules become heavier [say molecules numbered (7)–(10)], the chance of experimentally observing the exact lightest masses rapidly decreases, and the computed difference between exact and computed lightest masses becomes less important.

The evaluation of getting the correct heaviest mass is not as important under natural conditions. This is because heavy isotopes typically carry very low natural occurrence probabilities so that it is impossible to observe the exact heaviest isotopic variant of the molecule. Of course, when artificial isotopic abundances are enforced, obtaining the correct heaviest masses can become important, while obtaining the correct lightest masses can become unimportant. Since the current evaluation is using the natural isotopic abundances, we do not expect any method to provide correct heaviest masses. Indeed, because most methods are computing terms of an ID that are concentrated around a molecule’s average molecular mass, which is closer to the exact lightest mass under natural isotopic abundances, the mass range used for computing IDs usually will not include the heaviest masses. For biomolecules numbered (1)–(10), the differences between the exact heaviest masses and the heaviest masses computed by all methods considered are all of the same order of magnitude.

Displayed in the upper (lower) half of Table 4 are the relative differences of computed CGIDs derived molecular mass averages (standard deviations) to their theoretical values. Molecules numbered (1)–(10) in Table 2 are used with elements assuming isotopic abundances shown in Table 1. In terms of the average masses, MIDAs, JFC, and Emass have comparable errors and have slightly smaller errors than the other methods. In terms of mass standard deviations, MIDAs, JFC, and Emass have slightly smaller errors than the other methods. In principle, the accuracy of BRAIN might be improved by increasing the number of aggregated isotopic variants computed for each computed CGID. However, to accomplish this would require changing its default option. As mentioned earlier, to keep the benchmarking test simple, we only use the default option for each method considered. From Table 4, one can also infer that Qmass yields small errors for small and medium size molecules, but the error increases as the molecular mass increases.

Table 4

Coarse - Grained Isotopic Distribution Results using Naturally Occurring Isotopes

Difference in average mass
No.1	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	2.3e - 13	6.8e - 13	4.4e - 03	-4.5e - 13	9.1e - 05	-2.9e - 05	6.51e - 05	-4.6e - 13	2.3e - 12
(2)	1.8e - 12	3.5e - 11	3.2e - 01	-1.8e - 12	8.1e - 04	7.6e - 05	3.7e - 03	-7.3e - 12	3.6e - 12
(3)	-5.4e - 12	-5.4e - 12	8.2e - 02	-3.6e - 12	2.2e - 04	-5.1e - 05	5.9e - 03	5.4e - 12	0
(4)	-7.3e - 12	2.0e - 10	4.6e - 02	0	2.9e - 03	-7.3e - 04	-360	5.8e - 11	7.3e - 12
(5)	4.3e - 11	2.6e - 10	1.4e - 04	7.3e - 12	-3.1e - 03	1.8e - 04	3.7e - 03	-7.3e - 12	-2.9e - 11
(6)	0	1.3e - 10	1.7e - 06	-5.8e - 11	-4.1e - 03	3.1e - 03	-8.5e - 04	4.2e - 09	1.2e - 10
(7)	4.3e - 11	-2.4e - 09	-2.7e - 08	-1.4e - 11	-4.1e - 03	2.1e - 03	-3.9e - 03	-1.3e - 10	5.8e - 11
(8)	-2.9e - 11	1.6e - 09	-4.1e - 08	0	-5.1e - 03	-7.9e - 03	-1.0e - 02	-7.5e - 01	2.6e - 10
(9)	-3.5e - 10	7.2e - 09	-1.1e - 07	-3.5e - 10	1.7e - 02	7.7e - 03	-4.0e - 02	-9.7e - 03	7.6e - 10
(10)	-1.2e - 10	7.6e - 09	-1.2e - 07	-1.2e - 10	-1.1e - 01	3.3e - 02	-3.8e - 02	-4.4e + 02	-4.6e - 10
Difference in standard deviation
No.1	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	1.1e - 06	-4.2e - 10	1.2e - 02	1.1e - 06	1.6e - 06	-1.2e - 04	-3.6e - 04	1.1e - 06	1.1e - 06
(2)	6.5e - 06	-5.0e - 09	3.6e - 01	6.5e - 06	1.2e - 06	-1.8e - 04	-1.2e - 03	6.5e - 06	6.5e - 06
(3)	8.0e - 06	1.5e - 08	1.2e - 01	8.0e - 06	9.8e - 05	-3.3e - 05	2.2e - 03	7.9e - 06	8.0e - 06
(4)	7.2e - 06	2.5e - 08	7.3e - 02	7.2e - 06	-3.0e - 07	-4.6e - 04	-4.5e - 02	7.0e - 06	7.1e - 06
(5)	1.3e - 05	1.8e - 07	3.7e - 04	1.3e - 05	9.7e - 06	-2.7e - 04	-1.8e - 03	1.3e - 05	1.3e - 05
(6)	1.8e - 05	-1.9e - 07	2.3e - 05	1.8e - 05	-3.9e - 07	-9.0e - 04	-8.4e - 03	-2.7e - 06	1.7e - 05
(7)	2.0e - 05	-8.0e - 07	2.1e - 05	2.0e - 05	-2.7e - 07	-7.1e - 04	-7.5e - 03	2.2e - 05	2.0e - 05
(8)	2.5e - 05	2.1e - 06	2.5e - 05	2.5e - 05	4.4e - 06	-5.4e - 04	-8.7e - 03	-4.8e + 00	2.6e - 05
(9)	4.2e - 05	-7.8e - 06	4.1e - 05	4.5e - 05	-5.9e - 07	-1.5e - 03	-5.2e - 03	-9.9e - 02	3.9e - 05
(10)	4.8e - 04	-1.0e - 05	5.0e - 05	5.4e - 05	-1.2e - 05	-1.3e - 03	9.6e - 03	-1.4e + 02	3.8e - 05

1Reference number associated with a molecular formula (biomolecule).

We have also considered the possibility of deviations from the natural frequencies of occurrence of an element’s isotopes. Such customized modifications can be accomplished experimentally by a technique known as isotopic labeling [3], which is frequently employed in quantitative proteomics [44]. To mimic such a situation, we have computed CGIDs for various molecules assuming different carbon isotopic abundances: 99% 13C and 1% 12C as listed in Table 1. We then derive from the computed CGIDs the average molecular masses and standard deviations, and compare them to the corresponding theoretical values that can be analytically calculated.

The results of using the above mentioned customized carbon abundances are shown in Table 5. The differences displayed in Table 5 show that MIDAs, JFC, and Emass have the smallest errors. However, in terms of ID’s standard deviations, Mercury and Mercury5 yield comparable errors to MIDAs, JFC, and Emass. The results for Qmass are similar to the ones obtained in Table 4: in terms of average masses and standard deviations, it yields small errors for small to medium sizes biomolecules. Table 5 also shows that the current versions of BRAIN and NC are not able to compute IDs using the modified isotopic abundances for carbon. However, the developers of NC have mentioned how NC could be modified to handle stable isotope enrichment by partition of the elements of enriched isotopes away from the equatransneutronic isotopes groups [17]. This option is not currently available in NC. Also, the proposed solution reduces NC back to a polynomial method algorithm, which, if not efficiently implemented, can significantly influence the overall computation time. In BRAIN’s case, there are no reasonable IDs reported and it is difficult to speculate what might have happened.

Table 5

Coarse - Grained Isotopic Distribution Evaluation using Abundances for Carbon’s Isotopes of 99% 13C and 1% 12C

Difference in average mass
No.1	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	-6.8e - 13	3.9e - 12	-1.7e + 01	2.3e - 13	5.0e - 05	-4.0e - 05	3.3e - 02	-2.3e - 13	1.6e - 12
(2)	0	4.5e - 11	NR2	3.6e - 12	3.3e - 04	7.3e - 05	1.7e - 01	-1.8e - 12	4.5e - 12
(3)	1.8e - 12	2.2e - 10	NR	-1.8e - 12	-6.3e - 04	-3.1e - 04	3.1e - 01	-6.4e - 11	-2.4e - 11
(4)	-7.3e - 12	4.4e - 11	NR	0	3.9e - 03	-2.5e - 04	NR	-1.1e - 11	-1.1e - 11
(5)	-2.9e - 11	-5.8e - 11	NR	7.3e - 12	-5.2e - 03	6.4e - 04	1.8e + 03	-4.1e - 07	-7.3e - 12
(6)	5.8e - 11	-1.0e - 10	NR	2.9e - 11	-5.4e - 03	2.1e - 03	2.7e + 03	-2.9e - 11	2.9e - 11
(7)	1.4e - 11	6.8e - 10	NR	-7.3e - 11	-8.7e - 04	6.6e - 04	4.8e + 03	-4.9e - 07	1.4e - 10
(8)	3.2e - 11	-3.5e - 10	NR	8.7e - 11	-8.1e - 03	6.4e - 03	8.4e + 03	-1.1e + 02	1.2e - 10
(9)	0	4.2e - 09	NR	0	-3.7e - 03	8.7e - 03	1.7e + 04	-1.5e + 02	-4.1e - 10
(10)	2.3e - 10	7.7e - 09	NR	8.15e - 10	-1.2e - 01	5.4e - 03	2.4e + 04	-5.1e + 02	-1.2e - 10
Difference in standard deviation
No.1	MIDAsa	MIDAsb	BRAIN	Emass	Mercury	Mercury5	NC	Qmass	JFC
(1)	7.9e - 07	-3.3e - 10	-1.4e + 01	7.9e - 07	-2.6e - 08	-1.2e - 04	-4.9e + 00	7.4e - 07	7.9e - 07
(2)	5.9e - 06	5.9e - 09	NR2	5.9e - 06	2.2e - 07	-1.9e - 04	-1.1e + 01	5.9e - 06	5.9e - 06
(3)	7.6e - 06	1.4e - 08	NR	7.6e - 06	8.2e - 06	-1.3e - 04	-1.6e + 01	7.6e - 06	7.6e - 06
(4)	7.1e - 06	1.4e - 08	NR	7.1e - 06	-1.7e - 07	-4.7e - 04	NR	7.1e - 06	7.1e - 06
(5)	1.3e - 05	-2.0e - 07	NR	1.3e - 05	3.6e - 07	-2.8e - 04	5.2e + 00	1.2e - 05	1.3e - 05
(6)	1.7e - 05	-5.7e - 07	NR	1.7e - 05	-1.2e - 07	-9.2e - 04	6.6e + 00	1.8e - 05	1.8e - 05
(7)	2.1e - 05	-1.1e - 06	NR	2.1e - 05	-6.9e - 09	-7.2e - 04	8.6e + 00	1.9e - 05	2.0e - 05
(8)	2.2e - 05	6.3e - 07	NR	2.4e - 05	-1.7e - 06	-5.5e - 04	1.1e + 01	-1.1e + 02	2.5e - 05
(9)	4.0e - 05	1.0e - 06	NR	4.4e - 05	-3.5e - 06	-1.5e - 03	1.6e + 01	-2.0e + 02	5.1e - 05
(10)	3.3e - 05	1.0e - 05	NR	3.5e - 05	-1.5e - 05	-1.4e - 03	1.9e + 01	-2.3e - 04	6.8e - 05

1Reference number associated with a molecular formula (biomolecule).

2Nothing reasonable reported (NR).

Assessing Fidelity of Computed CGIDs and FGIDs

To evaluate the fidelity of CGIDs and FGIDs reported, we used 10 hydrocarbon molecules [numbered (11)–(20) in Table 2] because the “exact” CGIDs and FGIDs can be calculated for these molecules. Exact CGID is defined as follows. First, one merges isotopic variants that have the same nucleon number into one aggregated isotopic variant, whose corresponding molecular mass (MM) and occurrence probability are computed respectively from the probability-weighted sum of masses and from the sum of the probabilities of the isotopic variants merged. However, only aggregated isotopic variants having probability greater than 5e–12 were retained for accuracy evaluation. The exact FGIDs were obtained/defined similarly to the exact CGIDs, except that one merges only isotopic variants whose molecular mass differences are within some pre-specified mass accuracy, here set to 0.01 Da. The probability cutoff of 5e–12, for typical sample loads, probably already surpasses the detection capability of current mass spectrometer. Furthermore, it is also a small enough cutoff that ignoring terms below the cutoff has negligible effect in the ID profile.

Four quantities were then utilized to evaluate the fidelity of computed IDs. The first quantity is the difference in the numbers of terms (Δτ) kept by a computed ID and by its corresponding exact ID, be it the exact CGID or the exact FGID. The second quantity was the difference in the probability sums (Δχ), one from the computed ID and the other from the exact FGID (or the exact CGID). The third quantity was the root-mean-square differences of masses (σ ) between computed and the exact CGIDs (or the exact FGIDs).

In the equation above, m represents a computed mass term while each m represents a mass term in the exact FGID (or the exact CGID). N is the number of terms retained in the computed ID. That is, for every mass term in a computed ID, the closest mass term within the exact FGID (or the exact CGID) is found and their difference square is summed. The average of such sum of squares constitutes σ 2. The fourth quantity computed was the weighted correlation (ρ) between computed and exact IDs. The weighted correlation (ρ) is defined as follows. Let p(m ) and p(m ) be the terms of a computed ID and the corresponding exact FGID (or exact CGID), respectively. We first introduce the weight (w ) between a computed ID term (index i) and exact FGID (or exact CGID) term (index j) aswhere, in the above equation, min |m  − m |, is the minimum mass difference between a term (m ) from the computed ID and terms (m ) from the exact ID. The computed weights (w ) are then normalized by the normalization factor, W = ∑  w , by summing over all i terms from the computed ID that are close to a common term j in the exact FGID (or the exact CGID). The weighted correlation using the above definitions is given by

For CGIDs, ∈ is set to one Da, while for FGIDs, ∈ is set to 0.01 Da.

Using molecules numbered (11)–(20) in Table 2, we document the analysis results of the four quantities mentioned above in Table 6 (for CGIDs) and Table 7 (for FGIDs). For CGIDs, we include in Table 6 only four methods that largely satisfy the criteria for being a sound IDCM of application value. For FGIDs, only IC, MIDAs and MIDAs are included in Table 7 since they are the only methods that can do FGID-computing reasonably fast and without additional post processing.

Table 6

Coarse-Grained Isotopic Distribution (CGID) Fidelity Assessment Results τ is the Number of Terms in the Exact CGID Having Probability Greater than 5e - 12. Δτ is the Difference Between τ and the Number of Terms of a Computed CGID. Δχ is the Difference Between the Sum of Probability Terms from the Exact CGID and the Sum of Probability terms from the Computed CGID; σ is the Root-Mean-Square Differences of Masses Between Exact and Computed CGID, see Equation (11); U is the Number of Terms from the Computed CGID that are not with ± 2∈ (∈ = 1 Da) from any Terms in the Exact CGID; E is the Number of Terms in the Exact CGID that Have at Least One Corresponding Term in Computed CGID that are with ± 2∈; ρ is the Weighted Correlation Between Computed and Exact CGID

No.1	τ	Δτ	Δχ	σ m	U	E	ρ	Method
(11)	6	0	-5.6e - 16	5.1e - 13	0	6	1.0	MIDAsa
		0	-7.4e - 15	2.2e - 04	0	6	1.0	MIDAsb
		0	-4.7e - 04	1.2e - 14	0	6	0.99999988	Emass
		0	-4.4e - 16	2.1e - 05	0	6	1.0	JFC
(12)	7	0	-8.9e - 16	7.7e - 12	0	7	1.0	MIDAsa
		0	7.8e - 16	7.5e - 05	0	7	1.0	MIDAsb
		0	-8.3e - 04	2.6e - 14	0	7	0.99999963	Emass
		0	-5.6e - 16	1.3e - 05	0	7	1.0	JFC
(13)	12	0	-5.0e - 15	1.5e - 11	0	12	1.0	MIDAsa
		0	-2.2e - 14	3.3e - 05	0	12	1.0	MIDAsb
		0	-1.6e - 03	7.3e - 14	0	12	0.99999885	Emass
		0	-5.0e - 15	8.3e - 04	0	12	1.0	JFC
(14)	15	0	-1.0e - 14	1.2e - 11	0	15	1.0	MIDAsa
		0	1.5e - 14	2.3e - 05	0	15	1.0	MIDAsb
		0	-1.3e - 03	1.9e - 13	0	15	0.99999957	Emass
		0	-1.2e - 14	3.6e - 03	0	15	1.0	JFC
(15)	40	0	-9.8e - 14	1.0e - 11	0	40	1.0	MIDAsa
		0	1.3e - 12	6.3e - 06	0	40	1.0	MIDAsb
		0	-5.6e - 04	3.2e - 12	0	40	0.99999996	Emass
		0	-9.7e - 14	2.5e - 03	0	40	1.0	JFC
(16)	139	0	-9.6e - 13	4.4e - 10	0	139	1.0	MIDAsa
		0	3.8e - 12	6.3e - 06	0	139	1.0	MIDAsb
		0	-1.9e - 04	5.2e - 11	0	139	1.0	Emass
		0	-6.6e - 13	4.1e - 02	0	139	1.0	JFC
(17)	195	0	-2.0e - 12	5.5e - 10	0	195	1.0	MIDAsa
		0	1.3e - 12	6.3e - 06	0	195	1.0	MIDAsb
		0	-1.3e - 04	1.3e - 10	0	195	1.0	Emass
		0	-1.9e - 12	6.1e - 02	0	195	1.0	JFC
(18)	238	0	-3.0e - 12	9.4e - 10	0	238	1.0	MIDAsa
		0	2.5e - 11	6.4e - 06	0	238	1.0	MIDAsb
		0	-1.1e - 04	2.3e - 10	0	238	1.0	Emass
		0	-2.6e - 12	6.0e - 02	0	238	1.0	JFC
(19)	274	0	-4.1e - 12	1.2e - 09	0	274	1.0	MIDAsa
		1	2.5e - 11	6.1e - 02	0	274	1.0	MIDAsb
		0	-9.5e - 05	3.0e - 10	0	274	1.0	Emass
		0	-5.4e - 12	5.9e - 02	0	274	1.0	JFC
(20)	306	0	-4.8e - 12	1.6e - 09	0	306	1.0	MIDAsa
		0	2.6e - 11	6.8e - 06	0	306	1.0	MIDAsb
		0	-8.5e - 05	4.2e - 10	0	306	1.0	Emass
		0	-4.6e - 12	5.6e - 02	0	306	1.0	JFC

1Reference number associated with a molecular formula (biomolecule).

Table 7

Fine - Grained Isotopic Distribution (FGID) Fidelity Assessment Results τ is the number of terms in the exact FGID having probability greater than 5e - 12; Δτ is the difference between τ and the number of terms of a computed FGID; Δχ is the difference between the sum of probability terms from the exact FGID and the sum of probability terms from the computed FGID; σ is the root-mean-square differences of masses between exact and computed FGID, see Equation (11); U is the number of terms from the computed FGID that are not with ±2∈ (∈ = 0.01 Da) from any terms in the exact FGID; E is the number of terms in the exact FGID that have at least one corresponding term in computed FGID that are with ±2∈; ρ is the weighted correlation between computed and exact FGID

No.1	(ppm)2	τ	Δτ	Δχ	σ m	U	E	ρ	Method
(11)	307.25	6	0	-5.6e - 16	1.1e - 09	0	6	1.0	MIDAsa
			0	8.9e - 14	7.3e - 05	0	6	1.0	MIDAsb
			0	1.5e - 08	1.8e - 09	0	6	1.0	IC
(12)	153.63	7	0	-2.7e - 15	4.2e - 10	0	7	1.0	MIDAsa
			0	2.3e - 12	1.6e - 05	0	7	1.0	MIDAsb
			0	3.7e - 07	2.3e - 09	0	7	1.0	IC
(13)	30.72	13	1	-5.7e - 13	3.1e - 03	0	13	0.99999591	MIDAsa
			0	3.0e - 13	4.4e - 03	0	12	0.99999592	MIDAsb
			1	-2.9e - 07	3.1e - 03	0	13	0.99999591	IC
(14)	15.36	16	3	2.7e - 12	5.3e - 03	0	15	0.99937104	MIDAsa
			3	3.6e - 12	7.1e - 03	0	15	0.99937227	MIDAsb
			4	-3.5e - 07	5.1e - 03	0	16	0.99937103	IC
(15)	1.53	65	-6	-4.6e - 11	1.7e - 03	0	58	0.99927870	MIDAsa
			3	1.2e - 11	4.0e - 03	0	64	0.98806083	MIDAsb
			5	2.4e - 05	3.6e - 03	0	65	0.98803755	IC
(16)	0.15	291	-5	2.6e - 11	5.2e - 03	0	257	0.99999001	MIDAsa
			53	6.9e - 11	5.6e - 03	1	282	0.99958599	MIDAsb
			82	1.4e - 07	5.2e - 03	0	280	0.99998237	IC
(17)	0.077	500	-18	1.8e - 10	4.0e - 03	0	453	0.99785805	MIDAsa
			126	1.3e - 10	7.3e - 03	13	488	0.95950051	MIDAsb
			124	-1.6e - 08	4.5e - 03	0	496	0.99951891	IC
(18)	0.051	715	-16	-5.4e - 10	4.5e - 03	0	636	0.99466182	MIDAsa
			242	1.5e - 10	7.0e - 03	10	681	0.71069880	MIDAsb
			19	-5.0e - 08	4.3e - 03	0	690	0.99735244	IC
(19)	0.038	881	57	7.6e - 11	4.8e - 03	0	824	0.95007936	MIDAsa
			437	1.9e - 10	7.8e - 03	33	866	0.59270671	MIDAsb
			-26	-1.7e - 06	5.9e - 03	0	713	0.97935224	IC
(20)	0.031	1143	93	-1.4e - 10	5.7e - 03	0	1011	0.85638390	MIDAsa
			498	2.2e - 10	8.4e - 03	47	1092	0.63960000	MIDAsb
			-173	-1.7e - 05	4.6e - 03	0	838	0.98564325	IC

1Reference number associated with a molecular formula (biomolecule).

2Using a mass accuracy of 2∈ the equivalent resolution in parts per million.

For fidelity assessment of CGIDs, all four methods shown in Table 6 yield small Δτ and ρ values close to one. In terms of σ and Δχ, more differences are revealed. Emass always yields small σ , reflecting good fidelity in terms of mass locations, but seems to give a larger |Δχ|, reflecting less accuracy in amplitudes. JFC and MIDAs seem to yield less precise mass locations, evidenced by a larger σ , but seem to provide more accurate amplitudes, evidenced by a smaller |Δχ|. MIDAs yields both accurate mass locations and accurate amplitudes.

The values of Δχ and σ in Table 7 indicate that IC, MIDAs, and MIDAs report FGID terms with similar mass accuracy and with probability sums that are close to the expected value. For small to medium molecules, numbered (11)–(15), IC, MIDAs, and MIDAs have equivalently accurate results. For molecules numbered (16)–(20), IC and MIDAs have comparable performances, both slightly better than MIDAs. The values for Δτ indicates that MIDAs reports many more terms than expected in its computed FGID. Not expecting any leakage, MIDAs gains these extra terms mainly due to rounding errors associated with the DFFT numerical procedure.

The difference observed in Δτ for MIDAs is caused by the pruning and merging procedures employed by the algorithm. All the FGID terms computed by IC and MIDAs are within 2∈ from the exact FGID terms, which is shown by the number of unexplained term (U) being zero in Table 7. It is also true that most of the terms computed from MIDAs are within 2∈ from the exact FGID terms with the exception of molecules (17)–(20) where the number U ranges from 1 to 47. The computed weighted correlation also shows that for heavier molecules, (18)–(20), both IC and MIDAs produce FGIDs that are more similar to the exact FGIDs than MIDAs.

What causes MIDAs to perform worse here might be related to the fact that pinning the elemental masses to grid points may introduce appreciable mass errors while computing IDs for larger molecules. In the worst case scenario, the mass error introduced is apparently proportional to the number of atoms contained in the molecule. Even though MIDAs employs a mass rescaling [32] to bring the computed average masses and standard deviations close to their theoretical values, the linear mass rescaling is not sufficient to guarantee the full profile resemblance between the computed ID and the exact ID. The non-negligible discrepancy (indicated by the weighted correlation ρ not very close to one) between the computed FGID and the exact FGID for molecules (18)–(20) is reflecting this problem.

MIDAs Web Interface

MIDAs web interface http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/midas/index.html is user-friendly, but at the same time offers considerable flexibility. For example, in terms of the input molecule, the user may type in the box an elemental composition, a molecular formula, a peptide, or even a protein sequence. The program recognizes the input molecule in all formats above and extracts the corresponding elemental compositions for computing CGIDs and FGIDs. The isotopic abundances and elements’ masses can also be customized within the web interface. The user simply clicks on the “change” button to edit the abundance table of all elements. Other fields that can be easily customized and specified by the user are the charge of the input molecule and the cutoff probability. MIDAs displays both CGID and FGID together using user-specified accuracies, one for each. The “algorithm” drop down box allows the user to select either the FFT or the polynomial algorithms. The output, including the lightest mass, theoretical average mass, theoretical mass standard deviation, computed average mass, computed mass standard deviation, FGID peak list, and CGID peak list can be exported to a flat file by clicking on the “download output” button on the result page. There is also a contextual help for every functional button.

Conclusion and Outlook

The two algorithms introduced here, MIDAs and MIDAs, for the 25 molecules tested, seem to be able to compute IDs quickly and accurately. Between the two, MIDAs seems slightly more accurate. For CGIDs MIDAs appears to be faster (see Table 8), whereas for FGIDs they are of comparable speed. Both algorithms benchmark well with existing methods and stand out because of their ability to compute CGIDs and FGIDs using a user-specified accuracy. These two algorithms were also shown to accurately compute IDs for molecules labeled with stable isotopes, which was not the case for some of the methods evaluated. In summary, in terms of CGIDs derived average masses, MIDAs, MIDAs, JFC, and Emass yield smaller errors than other methods. In terms of CGIDs derived standard deviation, our investigation shows that MIDAs, MIDAs, JFC, and Emass yield smaller errors than other methods. When computing the FGID, MIDAs computes a FGID that better resembles the exact FGID than MIDAs using our evaluation gauges. Both algorithms described here were coded using the C++ programming language in a computer program called MIDAs that is available for download at http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/downloads/MIDAs.html. To make these algorithms widely accessible, we have made them available through a user-friendly web-interface at http://www.ncbi.nlm.nih.gov/CBBresearch/Yu/midas/index.html.

Table 8

Computation Time in Seconds (s) and Number of Terms Reported with MIDAs’s Computed Coarse-Grained (CG) and Fine-Grained (FG) Isotopic Distributions (ID) Using 1.0 Da and 0.01 Da Mass Accuracy, Respectively

MIDAsa
No.1	Number of terms CGID	CGID time(s)	Number of terms FGID	FGID time(s)
(1)	85	0.006	42	0.001
(2)	154	0.01	151	0.002
(3)	186	0.02	246	0.001
(4)	203	0.02	301	0.001
(5)	290	0.04	809	0.006
(6)	341	0.05	1269	0.01
(7)	423	0.1	1945	0.05
(8)	540	0.13	3145	0.06
(9)	820	0.14	6579	0.3
(10)	956	0.2	7850	0.4
(21)	3022	0.35	13834	0.8
(22)	5908	1.0	74994	1.0
(23)	2706	0.23	28508	2.1
(24)	617	0.05	2805	0.01
(25)	2623	0.21	18261	0.5
MIDAsb
No.	Number of terms CGID	CGID time(s)	Number of terms FGID	FGID time(s)
(1)	15	0.0006	29	0.025
(2)	29	0.001	114	0.041
(3)	38	0.001	193	0.043
(4)	42	0.001	241	0.08
(5)	78	0.002	740	0.08
(6)	95	0.002	1166	0.14
(7)	123	0.004	1784	0.14
(8)	157	0.004	2953	0.14
(9)	230	0.004	6527	0.3
(10)	257	0.01	7818	0.3
(21)	794	0.005	12405	0.4
(22)	1500	0.01	74994	1.0
(23)	706	0.01	23367	0.7
(24)	188	0.003	2209	0.2
(25)	698	0.01	16384	0.8

1Reference number associated with a molecular formula (biomolecule or inorganic compound).