Metabolite signal identification in accurate mass metabolomics data with MZedDB, an interactive m/z annotation tool utilising predicted ionisation behaviour 'rules'
- John Draper1Email author,
- David P Enot†1, 2,
- David Parker†1, 3,
- Manfred Beckmann†1,
- Stuart Snowdon1,
- Wanchang Lin1 and
- Hassan Zubair1
© Draper et al; licensee BioMed Central Ltd. 2009
Received: 20 February 2009
Accepted: 21 July 2009
Published: 21 July 2009
Metabolomics experiments using Mass Spectrometry (MS) technology measure the mass to charge ratio (m/z) and intensity of ionised molecules in crude extracts of complex biological samples to generate high dimensional metabolite 'fingerprint' or metabolite 'profile' data. High resolution MS instruments perform routinely with a mass accuracy of < 5 ppm (parts per million) thus providing potentially a direct method for signal putative annotation using databases containing metabolite mass information. Most database interfaces support only simple queries with the default assumption that molecules either gain or lose a single proton when ionised. In reality the annotation process is confounded by the fact that many ionisation products will be not only molecular isotopes but also salt/solvent adducts and neutral loss fragments of original metabolites. This report describes an annotation strategy that will allow searching based on all potential ionisation products predicted to form during electrospray ionisation (ESI).
Metabolite 'structures' harvested from publicly accessible databases were converted into a common format to generate a comprehensive archive in MZedDB. 'Rules' were derived from chemical information that allowed MZedDB to generate a list of adducts and neutral loss fragments putatively able to form for each structure and calculate, on the fly, the exact molecular weight of every potential ionisation product to provide targets for annotation searches based on accurate mass. We demonstrate that data matrices representing populations of ionisation products generated from different biological matrices contain a large proportion (sometimes > 50%) of molecular isotopes, salt adducts and neutral loss fragments. Correlation analysis of ESI-MS data features confirmed the predicted relationships of m/z signals. An integrated isotope enumerator in MZedDB allowed verification of exact isotopic pattern distributions to corroborate experimental data.
We conclude that although ultra-high accurate mass instruments provide major insight into the chemical diversity of biological extracts, the facile annotation of a large proportion of signals is not possible by simple, automated query of current databases using computed molecular formulae. Parameterising MZedDB to take into account predicted ionisation behaviour and the biological source of any sample improves greatly both the frequency and accuracy of potential annotation 'hits' in ESI-MS data.
Changes in the overall metabolite composition of living cells (metabolome) reflect a key end point in gene expression and make a major contribution to organism phenotype . Although, no single analytical platform can offer a fully comprehensive survey of the chemical diversity representing the metabolome, continuous improvements in mass spectrometry (MS) instrumentation have allowed development of relatively standardised metabolite profiling or fingerprinting procedures . A fundamental principle of mass spectrometry is the representation of metabolite features in any biological matrix by measurement of the spectrum of signals reflecting the mass to charge ratios (m/z) of their ionisation products. One advantage of MS over alternative spectroscopic methods such as Nuclear Magnetic Resonance (NMR) and Fourier Transform Infrared (FT-IR) is the opportunity to putatively annotate directly a spectral component by virtue of its atomic mass. In the context of a metabolomics experiment these 'first pass' annotations may be used to develop hypotheses relating to metabolite identity which are then tested by subsequent, more targeted, analytical chemistry methods. Traditional hyphenated MS profiling approaches provide simultaneous detection and quantification of discrete metabolite-derived peaks after chromatographic separation. In gas chromatography MS (GC-MS) a couple of hundred well resolved metabolite peaks are identified where possible by matching their positively charged ion spectrum (following fragmentation by electron impact in the gas phase) and column retention time to those of known standards . GC-MS methods are well established and robust but are limited to analysis only of volatile metabolites. Although derivatisation can increase the volatility of a wide range of metabolite classes such chemical modification further increases the complexity of any annotation process based on atomic mass.
Alternative profiling methods utilising liquid chromatography coupled to mass spectrometry (LC-MS) provide sensitive tools for the analysis of a wider range of metabolites with higher polarity, lower volatility and much larger mass range without a need for derivatisation . Resolved metabolite peaks can only be efficiently ionised when outside of the liquid phase under atmospheric pressure. Typical approaches include atmospheric pressure chemical ionisation (APCI) and electrospray ionisation (ESI). In contrast to electron impact, both methods allow 'soft' ionisation with little fragmentation in which a major product may be a pseudo-molecular ion comprising the protonated (+ ve ion data) or de-protonated (-ve ion data) parent molecule.
Peak detection and particularly spectral deconvolution in GC-MS and LC-MS are both technically challenging, time-consuming and very difficult to automate; LC-MS particularly has been hampered by poor analyte peak resolution and retention time variability which confound metabolite peak alignment [5, 6]. As a consequence, recent approaches to LC-MS profiling have concentrated on the non-targeted quantification of all detected ionisation products above a pre-set noise threshold which are compiled subsequently as several thousand un-annotated variables [7–9]. The variables in this high dimensional data consist simply of an ion mass and its retention time. Such data tables will contain much redundancy as the data acquisition rate of most instruments will be sufficient to generate at least 3–5 scans across each eluting peak. In addition, even with the extra resolving power of ultra high pressure liquid chromatography (UHPLC) a large number of peaks will still overlap or even co-elute in many regions of the chromatogram which demands the use of powerful signal alignment software for data pre-processing.
An alternative approach to LC-MS profiling is to generate a metabolite fingerprint representation without using the chromatographic dimension in which data variables are simply the detected mass values . In such approaches [2, 10–16] the sample may be dissolved in an appropriate solvent and injected directly into the ion source (Direct Injection Mass Spectrometry; DIMS) or infused as a 'plug' flow using a HPLC system without a chromatography column (Flow Injection Electrospray Ionisation Mass Spectrometry; FIE-MS). Data representations may take two forms depending mainly on the accuracy and subsequent resolving power (mass/mass accuracy) of the system mass analyser. In nominal mass fingerprinting methods the intensities of all masses from all scans are integrated within pre-defined ranges (mass bins). Linear ion traps with quadrupole (Q) detectors are extremely robust and offer rapid scanning over large m/z ranges. However, as full-scan mass resolution is often less than 1000 all signals are binned to the nearest nominal mass value from around m/z 50 up to as high as m/z 2000, depending on instrument; under these circumstances metabolites are considered to be mono-isotopic and each mass bin could obviously contain signals derived from several metabolites. Despite these limitations nominal mass fingerprinting methods have attracted considerable interest as a first pass investigation tool as data pre-processing is quick and trivial with minimal likelihood of error and with a short cycle time (typically less than 5 minutes) they are suitable for experimental designs which require a high throughput of samples . Time-of-flight (TOF) and hybrid Q-TOF mass analysers exhibit an average resolving power as high as 10,000 which, with internal calibration, may be sufficient to allow binning of signals in finer mass ranges (0.1 amu). This level of resolution however is not sufficient to distinguish between a large proportion of metabolites of different elemental composition which remain effectively isobaric at the maximum mass resolution achievable. With the advent of ultra-high accuracy mass analyzers such as the Orbitrap (resolution ~100,000) and Fourier Transform Ion Cyclotron Resonance (FT-ICR) analyzer (resolution potentially in excess of 1,000,000) this latter problem is greatly reduced [17–19]. Particularly using FT-ICR-MS instruments there is the opportunity to generate an accurate mass fingerprint (often referred to as a profile as variables are now considered effectively discrete) of an infused sample extract [20–22]. In such data representations mass 'peaks' are detected in the raw data and the centroid mass intensity calculated. The resolving power and sensitivity of these ultra-high accuracy mass analyzers is such that metabolite signals from molecules containing naturally abundant elemental isotopes (e.g. 13C, 41K, 15N, 17O,34S and 37Cl) are visible in the data [17, 18].
The past few years has witnessed a rapid increase in the number of publications describing experiments using particularly electrospray ionisation mass spectrometry (ESI-MS) profiling or fingerprinting at a range of mass accuracies. This activity has been accompanied by the recognition that the majority of m/z variables in most complex biological matrices cannot be matched unambiguously to metabolites of identical elemental composition found in publicly accessible databases [9, 19–27]. With the ultra-high mass accuracy achievable using FT-ICR-MS methodology [17, 18] this might seem surprising for metabolites with a mass < 500 Da for which the number of isobaric molecules based on calculated elemental composition should be very limited. Several factors contribute to this problem. In the simplest situation, given that many natural metabolites (particularly from plant sources) remain to be structurally identified, the relevant molecules may simply not be present in current databases. A more complex issue stems from the fact that all metabolites have a chance of containing one or more natural isotopes [19, 23, 24] of constituent atoms (e.g. 13C, 41K, 15N, 17O,34S and 37Cl); metabolite signals thus are rarely present as a single mono-isotopic peak and the spectral isotope distributions often results in measured masses of significant intensity that by chance are undistinguishable (ie not resolvable) from completely unrelated chemistry [18, 24]. A final factor relates to the diversity of potential ionisation products generated in the ion source in addition to parental pseudo-molecular ions (i.e. [M+H]1+ and [M - H]1-). Thus, in both fingerprinting and profiling approaches it is common to find adducts with either sample matrix components (e.g. in positive ion data: [M + Na]1+, [M + K]1+ and [M + 2Na - H]1+ and negative ion data: [M + Cl]1- and [M + Na - 2H]1-) or mobile phase solvents (e.g. [M + NH4]1+ in positive ion data). Additionally, although ESI is traditionally considered a 'soft' ionisation method a number of more sensitive metabolites will rearrange and fragment with neutral loss of common moieties such as water from a hydroxyl group [M - H2O + H]1+ or ammonia from an amine group [M - NH2 + H]1+ or formate from a carboxyl group [M - CO2H + H]1+. Clearly, it is possible also that a range of potential ionisation products will contain various combinations of isotopes, adducts and neutral losses. Finally, ionisation products also can carry more than one charge (e.g. [M + 2H]2+) or can, through non-covalent interactions, produce stable molecular clusters, for example dimers, in which only one component is charged (e.g. [2M + H]1+). The potential to form different ionisation products will be dependent on the structure and physical properties of individual metabolites.
It can be concluded that a key aspect of any variable annotation strategy is to assess the relationships between signals and to base putative annotations on models that best fit all of the data. The presence of isotopic atoms and either loss, or gain, of selected chemical moieties all afford specific changes in accurate mass which mean that the exact mathematical relationships between detected ions are informative. Natural isotope abundances in metabolites reflect the percentages at which isotopes of an element occur naturally in the environment. Thus analysis of isotope relative abundance patterns provide a useful approach to 'filter' unlikely annotations based on molecular formulae calculated from accurate mass measurements [19, 23]. Annotation strategies have been reported recently incorporating algorithms interrogating some, but not all, of these atomic mass mathematical relationships [19, 22, 23, 25, 26]. In most instances the adopted strategy centred on an analysis of the full spectrum of aligned peaks resolved at ultra-high mass accuracy in the biological matrices under comparison. A less computationally intensive strategy more suitable for high throughput experiments is to focus the annotation exercise on m/z signals which make significant contributions to models explaining the metabolome differences between two biologically different sample classes [10, 14–16, 22, 28]. In nominal mass metabolite fingerprint data [10, 14–16, 28] such experiments make the simple assumption that 'explanatory' metabolites will be represented by several mathematically-related ionisation products (e.g. pseudo-molecular parental ion, 13C isotope, salt adducts) which additionally exhibit tight co-variance in a simple correlation analysis; where LC-MS profiling is used ionisation products related to the same metabolite should additionally be found at the same retention time [14, 29]. Similar approaches are clearly possible using accurate mass data [7, 8, 24, 27].
It is time consuming and often computationally intensive to pre-process LC-MS profile/fingerprint data to a level where the researcher is absolutely confident of understanding the origin of any ionisation product (i.e. whether signals of interest are parental pseudo-ions, salt adducts, solvent adducts, neutral losses, clusters). This is a particular problem when high-throughput is required or for exploratory analyses of a new biological matrix. An alternative approach described here is the use of a database strategy that will allow searching based on all potential ionisation products that may be predicted to form during electrospray ionisation.
Searching for likely annotation candidates based on accurate mass information in publicly accessible databases is in itself time consuming as individual database coverage of natural chemistry varies and so a comprehensive search requires query of several. Unfortunately, with few exceptions, databases with appropriate metabolite mass information can contain much redundancy, resources for curation are often limited (consequently it is not uncommon to find mistakes relating to mass values, molecular formulae and structure) and some of the entries relate to ionic states, often from interactions with salts. To overcome some of these problems we have developed MZedDB  a database which uses an archive in a common format of all metabolite 'structures' derived from several widely used and publicly accessible databases. The value of such an approach to avoid ontology problems was recognised recently by a consortium of yeast researchers . Using a set of 'rules' derived from structural information and physical properties (such as number of H-bond acceptor/donors, number of OH/COOH/NH2 groups, number of acidic H or basic O in molecule) MZedDB will generate a list of potential adducts and neutral loss fragments that are likely to form for each structure and calculates on the fly the accurate mass of every potential ionisation product which provide targets for searches based on accurate mass. Starting with a list of m/z signals MZedDB supports a range of manual or semi-automated annotation strategies based on either m/z mass or predicted molecular formulae at a range of mass resolutions. In the present article we describe the development of MZedDB and illustrate several typical applications relating to annotation of ESI-MS data which take into account both predicted and measured ionisation behaviour.
Results and discussion
Developing MZedDB construction strategy
Comparison of data fields useful for LC-MS m/z signal annotation in a selection of online databases
Major Data Class
Databases with field
H, MT, Kp
K, P, Ch
H, P, Ch, Cy, MB, Cs
P, Ch, Cs
H, P, Ch, Cy, MB, ML, MC, MT, K, Cs, Kp
H, P, Ch, Cy, MB, ML, MC, MT, K, Cs, Kp
H, MB, Kp
K, H, P, MB, MT, Kp
ML, MB, MT
ML, MB, MT, Cs
H, P, Cs
H, P, Cy
K, H, Cy, MC
K, Cy, MC
K, Cy, MC
Development of 'rules' to link metabolite structures represented in MZedDB to potential ionisation products
SMILES, a unique ID and hyperlinks are the only information permanently stored in MZedDB for each metabolite entry accessed from external databases. SMILES can be used to generate a structural representation of the metabolite in question using the Java applet Jmol if required. Automated analysis of SMILES following a set of simple 'rules' can be used to detect the presence of specific moieties that are likely to effect potential ionisation behaviour; for example the possession of NH2 or COOH groups which could be lost as ammonia or formate respectively, or the presence, for example, of hydrogen bond acceptor or donor groups which could affect adduct formation. Additionally, the SMILES can be used to create a molecular formula from which an accurate monoisotopic accurate mass can be calculated.
Example default ionisation product mass calculation rules
Nch = 1
Nacc>0 AND Nch = 0
Nacc>0 AND Nch = 0
Nacc>0 AND Nch = 0
Nacc>0 AND Nch = 0
Nnhh>0 AND Nch = 0
Ncooh>0 AND Nch = 0
Noh>0 AND Nch = 0
Nch = -1
Ndon>0 AND Nch = 0
Ndon>1 AND Nacc>0 AND Nch = 0
Nacc>0 AND Nch = 0
Ndon>1 AND Nacc>0 AND Nch = 0
pMZedBD architecture and functionalities
Parameterisation of MZedDB searches on the basis of potential ionisation behaviour
Number of putative annotations of FT-ICR-MS signals using MZedDB
No. of potential hits at 3 ppm
Potential m/z annotations
Additional m/z annotations
[M + H]1+ m/z 159.0764:
4-methylene-L-glutamine; [M+NH4]1+ m/z 159.0764:
OR Gentianaine OR 6-oxo-
[M + Na]1+ m/z 166.0839:
[M + H]1+ m/z 172.0005:
[M + K]1+ m/z 172.0007:
Aspartic acid OR 2-hydroxy-succinamate
No annotation hits
No annotation hits
No annotation hits
[M + 2K - H]1+ m/z 268.9461:
Diketogulonic acid OR 7 other metabolites
Prevalence of potential common isotopes and adducts signals in FT-ICR-MS data derived from analysis of extracts of various biological tissues.
Signal relationship as a percentage of signals as common adducts/isotopes
Potato tuber (polar)
Potato tuber (non polar)
Ideally, MZedDB should be customised intelligently by selection of a subset of PIPs known to be abundant in the matrix under study and searches constrained to databases representing the organism in question. In circumstances where new biological matrices are investigated it may be helpful initially to use the MZedDB default settings in which three parameterisation boxes are automatically 'checked' (Adducts default selection [positive mode]; Apply adducts formation rules; Only include C, H, N, O, P, S) and a search of all database entries is activated from the drop down menu. If negative ion data is under investigation then Adducts default selection should be changed to negative mode. It may be helpful to change the default mass accuracy of 1 ppm to 5 ppm if accurate mass data from Time-of-flight (TOF) or hybrid Q-TOF instruments is used. If large numbers of annotation suggestions are returned then searches can be customised initially by restricting to entries in specific external databases by selecting one or more (using shift key) from the dropdown menu. Further reduction in numbers of annotation suggestions can be achieved by 'unchecking' the default selection boxes and then selecting a customised list of adducts from the pull down menu. Typically, in positive ion data selected adducts would include [M+H]1+, [M+Na]1+, [M+K]1+ and [M-H2O+H]+ and in negative ion data [M-H]1-, [M+Cl]1-, [M+Na-2H]1- and [M+K-2H]1- and [M-H2O-H]1-. Customisation could be extended sequentially as desired to include less common neutral losses (e.g. [M-NH2+H]1+ or [M-CO2H+H]1+), ion clusters (e.g. [2M+K]1+) or ions with multiple charges.
Guiding m/z annotation decisions by examination of m/z signal relationships
In high throughput LC-MS fingerprinting/profiling techniques utilising electrospray ionisation (ESI-MS) it is a common precaution to use analytical instruments at lower than maximal mass resolution to avoid problems with data alignment [21, 24]. As a result, m/z annotation suggestions may often include a large number of structurally diverse but effectively isobaric candidates when only parental pseudo-ions (i.e. [M + H]+ or [M - H]- m/z signals) are considered as ionisation products. Although recent papers have stressed the importance of recognising isotopes to avoid miss-annotation of m/z signals in accurate mass LC-MS data [9, 18, 21, 24–26] the present paper shows clearly that ionisation products other than parental pseudo-ions in reality can account for upwards of 40% of the signals in many accurate mass spectra. Against this finding it is unsurprising that many accurate mass m/z signals lack annotation 'hits' if structures other than parental pseudo-ions are not included in the searchable databases. Thus we conclude that annotation decisions can be usefully guided by both determining from the outset whether any signal in question is actually an adduct, isotope or neutral loss fragment or indeed any combination of all three. Takahashi et al.  report the ability to modify a search of the species-specific KNApSAcK database dependent upon the presence of sodium, potassium or ammonium adduct ions derived from the solvent used for sample preparation. Other researchers have downloaded the general metabolite database KEGG LIGAND and calculated the exact mass of 7 possible adduct ions for each entry to provide query targets in a customised database . MZedDB consolidates and expands on this strategy by harvesting SMILES and calculating the accurate mass of a much more comprehensive list of metabolites from a range of databases, as well as including hyperlinks to rapidly access detailed information from entries in the original databases. By restricting MZedDB searches to entries derived from only a selected number of the original databases m/z annotations can be refined in a species specific context. An additional advantage of MZedDB centres on the ability to extend annotation searches to include a much larger range of adducts, neutral loss fragments and clusters by implementation (using on the fly calculations) of a simple set of Potential Ionisation Product 'rules' which effectively filter searches to exclude unlikely candidates based on SMILES information and which can be customised in line with any prior knowledge concerning the matrix under analysis (e.g. salt content of tissue, extraction or HPLC solvent). Ultimately, metabolite identity assignment is crucial to derive biological meaning from metabolome modelling experiments and thus data annotation tools such as MZedDB which allow more intelligent searching need to be expanded to cover more comprehensively a wider range of natural chemistry and take into account the characteristics of new biological matrices. For example, to improve the analysis of human urine samples it is anticipated that future implementations of MZedDB will be extended to include any common biotransformations of parent metabolites (e.g. glucuronidation, sulphation, and de-methylation) in attempts to link ionisation products to original metabolites.
Generation of FT-ICR mass spectrometry data
The majority of data used to evaluate the functionality of the MZedDB annotation systems was generated on a Thermo LTQ instrument fitted with a 7-Telsa FT-ICR mass analyser. Extract preparation and basic instrument set up has been described previously  but the present study utilised an Advion Nanomate chip-based direct infusion nanospray ionisation to introduce the sample . Typical nanospray conditions comprised 200 nl/min flow rate, 0.5 psi back pressure, and +1.6 kV (positive ion data) or -1.6 kV (negative ion data) electrospray voltage, controlled by Chipsoft software (Advion). FT-ICR-MS parameters included an automatic control gain setting of 1 × 105 and a mass resolution of 100,000 and data was recorded for 6 min per replicate infusion using Xcalibur software (Thermo Scientific). Blank quality control samples comprising extraction solvents were interspersed at random into the run sequence to monitor instrument performance and detect system peaks. Prior to any statistical analysis the data was log transformed to reduce the chance of the higher intensity peaks dominating in any multivariate analyses.
An FT-ICR-MS data set representing flounder liver extract generated by the SIM-stitching method  was kindly supplied by Prof. Mark Viant and Dr Andrew Southam (School of Biosciences, University of Birmingham, UK).
Data Processing and Statistical Analysis
Data within each biological matrix class (ie leaf tissue, urine, potato tuber etc) were aligned and any peaks not represented in 70% of replicates were removed from the matrix. Two basic types of data sets were available for the project. The first type comprised a population of biological replicates representing a typical 'tissue' class; to provide some diversity the samples included fish liver (flounder), human urine, human blood plasma, plant leaves (Brachypodium distachyon) and potato tubers. The second type of samples were characterised by being part of a larger experiment that contained several biologically related sample matrices. For example a series of leaf sections harvested at different time points following infection of the grass B. distachyon with the rice blast pathogen  provided a complex sample series involving the interaction of two organisms. All sample classes could be used to generate FT-ICR-MS spectra containing thousands of resolved accurate mass signals (full matrix analysis) to provide typical annotation targets. FT-ICR-MS matrices derived from the infected Brachypodium samples additionally allowed both a detailed analysis of signal correlation behaviour and provided sample populations for further supervised multivariate data mining  in order to identify a small sub-set of 'explanatory' features (m/z) able to discriminate between the different diseased states.
All statistical tests were carried out in the R environment using the FIEMSpro metabolomics data analysis package  which is web accessible . Explanatory feature selection was performed using Random Forest [15, 16, 27]. Signal correlation analysis was carried out by the Pearson correlation method and performed on the explanatory m/z obtained after feature selection. Hierarchical cluster analysis based on the correlation coefficient was employed to identify the set of clusters, which satisfy some setting, for example, a signal correlation coefficient larger than 0.75. Determining the mathematical relationships between m/z is performed in R as described . This code searches for operator predetermined mass differences between measured accurate masses at an adjustable sensitivity (examples of mass difference searches are shown in Figure 4). In theory any mass difference can be searched for providing the operator knows the exact expected mass difference between the measured masses. This process is important to indicate possible isotope signals present in the matrix for which a prediction would not be wanted, and as an indication of the relationship between ionisation products within the matrix.
Developing MZedDB database
MZedDB is designed as a relational database which is implemented in a MySQL server running on PowerPC 1.8 ghz with 2 GB ram. The web-based interface is publicly available . All calculations were performed within the same environment using Perl and Unix shell commands to generate MZedDB tables and the ChemAxon software suite  to manipulate chemical structures. Interaction between MySQL database, user queries, results display and external applications was implemented using PHP. The Isotope Pattern Calculator (IPC v 1.3)  was employed to simulate isotopic patterns based on element atomic masses and probabilities. The integrated molecular formula generator  is a modified version of the tool available from Tobias Kind .
Overview of MZedDB table structure
MZedDB comprises a set of several tables containing relevant information from a range of web-accessible databases containing molecular information relating to metabolites (see Table 1 in Results). Due to the lack of a universal ontology for representing metabolite chemical information  and the non existence of a standalone platform where metabolite information is stored, our strategy relies on harvesting molecular information from web accessible metabolite databases. This information is then aggregated according to the actual chemical structures contained in each database, rather than using ambiguous annotation based on preferred names, synonyms or even CAS registry number. At the time of writing, only a few metabolite databases have been deposited and regularly updated to PubChem, making the use of its standard chemical information for searching PIP fairly restrictive. Finally, links between entries in different databases were not used because they may either be incomplete or somehow erroneous, irrespective of their level of curation.
Generating MZedDB dbcompound table
Chemical structures and associated pathway and reactions when available were downloaded from the corresponding websites (see Table 1 in Results section). Molecular structures, database entry ID, names and synonyms were parsed using bespoke Perl scripts. Because of the various formats used to represent molecular structure in each source (e.g. InChI http://www.inchi.info/, SMILES http://www.daylight.com/smiles and SDF http://www.mdli.com/solutions) and the inherent difficulties associated with matching chemical information from different chemoinformatics software, we adopted a common SMILES (Simplified Molecular Input Line Entry System) format generated with MolConverter (v. 5.l.4 ).
To enable comparisons across databases, we followed a multiple step strategy to clean up each entry. Ambiguities related to salts and complexes were removed by keeping the largest fragment. Structures containing less than 6 atoms, elements other than (C, H, N, O, P, S, K, Cl, Br, F) and ambiguous fragments (e.g. R) were then systematically excluded. A third stage consisted of reducing a major source of discrepancies between databases resulting from the dual representation of neutral and ionic forms of the same entity by removal (positive ion) or addition (negative ion) of a hydrogen atom. Finally, each molecular entity was aggregated according to: (1) Any exact match identification number; (2) A selected name from the list of names given by the original source; (3) Aggregated list of all synonyms; (4) Entry SMILES; (5) Canonical SMILES. Whenever available, stereochemical information was retained during the curation process, giving rise to MZedDB redundancies (e.g. glucose has several entries corresponding to the full or partial description of the configuration of its 5 asymmetric centres). The provision of Canonical SMILES (i.e. representation without specified stereochemistry) allow proposition of alternative entries with similar chemical structure.
Other MZedDB tables
Alongside the MZedDB compound table, dblinks functions to regroup the ID in the primary metabolite database for each MZedDB entry in order to generate dynamic hyperlinks to the original database. The table dbmetrule1 contains properties of each MZedDB entry that are required to perform the Potential Ionisation Product search. This table encapsulates molecular formula generated from the SMILES and the accurate monoisotopic mass computed to seven decimal places to avoid inaccuracies while rounding at a lower precision level. The remaining columns in dbmetrule consist of the overall charge and cardinalities of a predefined set of chemical groups necessary to apply adduct formation rule restrictions during Potential Ionisation Product (PIP) search. Metabolic pathway related information is contained in two tables. The table reac1 regroups the reaction description from each primary database (i.e. identifier, names, enzyme annotations, pathway identifiers and names). The table reac2 provides links for each MZedDB entry to reactions in which the compound is known to be involved.
MZedDB search capabilities
As the initial goal of MZedDB was not to provide a chemoinformatics environment for metabolomics, any metabolite-related queries were reduced to elementary searches based on chemical composition, molecular weights and synonyms. A metabolite entry reflects chemical properties relevant to m/z annotation and provides links to primary database sources for known pathway/enzymatic information to help biological interpretation. Links to entries with the same skeleton (canonical SMILES) and molecular formula are also provided. A typical output file based on a query for information on pipecolic acid is presented [see Additional file 1]. The user can interactively check its potential ionisation products and verify their exact isotopic pattern distributions to corroborate with experimental data.
Primarily, MZedDB focuses on extensive PIP search capabilities  using an extendable list of potential adducts (in both positive and negative ionisation modes) and chemistry based formation rules to avoid impossible or highly improbable PIPs entering the list of potential candidates . For example the possession of an amine group that will allow the neutral loss of ammonia or the presence of strong hydrogen bond donors or acceptors which affect interactions with salts and solvents to generate molecular adducts. As MZedDB allows either a default selection, or user-specified selection (including no rules) of adducts, neutral losses, clusters, charged states and number of charges, only a subset of 'rules' are applied 'on the fly' depending on which options are 'checked' when setting up a query. The end result is that when MZedDB is queried with one or more accurate masses only the masses of parent metabolite derivatives that satisfy the chosen set of PIP rules are available as annotation targets.
The annotation search can be parameterised by defining the level of mass resolution (in ppm), simply an m/z range or nominal mass given the truncation value for binning. Further parameterisation of MZedDB is achieved by limiting to metabolites found only in specific organisms, groups of organisms or databases. Databases restrictions are obviously straightforward. Limiting search to metabolites found in specific organisms is achieved by linkage to KEGG information. Additionally, the search can be extended to the more general generation of acceptable PIP molecular formula given the overall charge, the number of non covalently bound units and a mass precision. Potential molecular formula of the parent compound deduced from the PIP table are used to search a local installation of PubChem (to avoid heavy traffic to the server).
Availability and requirements
Project name: MZedDB
Project home page: http://maltese.dbs.aber.ac.uk:8888/hrmet/index.html
Operating system(s): Multiple platform (tested on Windows and Mac OS X).
Programming language: php, MySQL, Perl, C++
Other requirements: Java 1.3.1 or higher
Any restrictions to use by non-academics: None
We thank Kathleen Tailliart for technical support for FT-ICR-MS analysis and Robert Darby for maintenance of laboratory infrastructure. We thank Andrew Southam and Mark Viant (University of Birmingham, UK) for supplying a flounder liver FT-ICR-MS data matrix and for commenting on the utility of MZedDB for environmental data sets. We thank Augustin Scalbert and colleagues (Unité de Nutrition Humaine, INRA, Centre de Recherche de Clermont-Ferrand/Theix 63122 Saint-Genes-Champanelle, France) for commenting on the utility of MZedDB for human nutrition studies. We thank Tobias Kind (UCL Davis, USA) for discussions regarding refinement of the Molecular Formula Generator. JD and DE were supported by Aberystwyth University. MB was supported jointly by a Research Councils UK Fellowship and Aberystwyth University. DP, HZ and WL were supported by the UK Biotechnology and Biological Sciences Research Council grant BB/D006953/1. SS was supported by a Ph. D studentship from the UK Food Standards Agency.
- Fiehn O: Metabolomics – the link between genotypes and phenotypes. Plant Mol Biol 2002, 48: 155–171.View ArticlePubMedGoogle Scholar
- Dunn WB, Bailey NJC, Johnson HE: Measuring the metabolome: current analytical technologies. Analyst 2005, 130: 606–625.View ArticlePubMedGoogle Scholar
- Lisec J, Schauer N, Kopka J, Willmitzer L, Fernie AR: Gas chromatography mass spectrometry-based metabolite profiling in plants. Nat Protocols 2006, 1: 387–396.View ArticlePubMedGoogle Scholar
- De Vos RC, Moco S, Lommen A, Keurentjes JJ, Bino RJ, Hall RD: Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry. Nat Protocols 2007, 22: 778–791.View ArticleGoogle Scholar
- Smith CA, Want EJ, O'Maille G, Abagyan R, Siuzdak G: XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem 2006, 78: 779–787.View ArticlePubMedGoogle Scholar
- Lange E, Tautenhahn R, Neumann S, Gropl C: Critical assessment of alignment procedures for LC-MS proteomics and metabolomics measurements. BMC Bioinformatics 2008, 9: 375.PubMed CentralView ArticlePubMedGoogle Scholar
- Idborg H, Zamani L, Edlund PO, Schuppe-Koistinen I, Jacobsson SP: Metabolic fingerprinting of rat urine by LC/MS Part 2. Data pretreatment methods for handling of complex data. J Chromat B-Analytical Technologies in the Biomedical and Life Sciences 2005, 828: 14–20.View ArticleGoogle Scholar
- Keurentjes JJ, Fu J, de Vos CH, Lommen A, Hall RD, Bino RJ, van der Plas LH, Jansen RC, Vreugdenhil D, Kornneef M: The genetics of plant metabolism. Nat Genet 2006, 38: 842–849.View ArticlePubMedGoogle Scholar
- Sana TR, Roark JC, Li X, Waddell K, Fischer SM: Molecular formula and METLIN personal metabolite database matching applied to the identification of compounds generated by LC/TOF-MS. J Biomol Tech 2008, 19: 258–266.PubMed CentralPubMedGoogle Scholar
- Beckmann M, Parker D, Enot D, Duval E, Draper J: High throughput metabolome fingerprinting using Flow Injection Electrospray Mass Spectrometry. Nat Protocols 2008, 3: 486–504.View ArticlePubMedGoogle Scholar
- Aharoni A, De Vos CHR, Verhoeven HA, Maliepaard CA, Kruppa G, Bino R, Goodenowe DB, et al.: Nontargeted metabolome analysis by use of Fourier Transform Ion Cyclotron Mass Spectrometry. OMICS 2002, 6: 217–234.View ArticlePubMedGoogle Scholar
- Allen J, Davey HM, Broadhurst D, Heald JK, Rowland JJ, Oliver SG, Kell DB: Metabolic footprinting: a high throughput, high-information approach to cellular characterisation and functional genomics. Nat Biotech 2003, 21: 692–696.View ArticleGoogle Scholar
- Scholz M, Gatzek S, Sterling A, Fiehn O, Selbig J: Metabolite fingerprinting: detecting biological features by independent component analysis. Bioinformatics 2004, 20: 2447–2454.View ArticlePubMedGoogle Scholar
- Catchpole GS, Beckmann M, Enot DP, Mondhe M, Zywicki B, Taylor J, Hardy N, Smith A, King RD, Kell DB, Fiehn O, Draper J: Hierarchical metabolomics demonstrates substantial compositional similarity between genetically-modified and conventional potato crops. Proc Natl Acad Sci USA 2005, 102: 14458–14462.PubMed CentralView ArticlePubMedGoogle Scholar
- Enot DP, Beckmann M, Overy DP, Draper J: Predicting interpretability of metabolome models based on behavior, putative identity and biological relevance of explanatory signals. Proc Natl Acad Sci USA 2006, 103: 14718–14723.View ArticleGoogle Scholar
- Beckmann M, Enot DP, Overy DP, Draper J: Representation, comparison, and interpretation of metabolome fingerprint data for total composition analysis and quality trait investigation in potato cultivars. J Agric Food Chem 2007, 55: 3444–3451.View ArticlePubMedGoogle Scholar
- Alan G, Marshall AG, Hendrickson CL: High-Resolution Mass Spectrometers. Ann Rev Anal Chem 2008, 1: 579–599.View ArticleGoogle Scholar
- Breitling R, Pitt AR, Barrett MP: Precision mapping of the metabolome. Trends Biotech 2006, 24: 543–548.View ArticleGoogle Scholar
- Kind T, Fiehn O: Metabolomics database annotations via query of elemental compositions: Mass accuracy is insufficient even at less than 1 ppm. BMC Bioinformatics 2006, 7: 234.PubMed CentralView ArticlePubMedGoogle Scholar
- Madalinski G, Godat E, Alves S, Lesage D, Genin E, Levi P, Labarre J, Tabet JC, Ezan E, Junot C: Direct introduction of biological samples into a LTQ-Oritrap hybrid mass spectrometer as a tool for fast metabolome analysis. Anal Chem 2008, 80: 3291–3303.View ArticlePubMedGoogle Scholar
- Han J, Danell RM, Pate JR, Gumerov DR, Scarlett CO, Speir JP, Parker CE, Rusyn I, Zeisel S, Borchers C: Towards high-throughput metabolomics using ultrahigh-field Fourier transform ion cyclotron resonance mass spectrometery. Metabolomics 2008, 4: 128–140.PubMed CentralView ArticlePubMedGoogle Scholar
- Taylor NS, Weber RJ, Southam AD, Payne TG, Hrydziuszko O, Arvanitis TN, Viant MR: A new approach to toxicity testing in Daphnia magna : application of high throughput FT-ICR mass spectrometry metabolomics. Metabolomics 2008, 5: 44–58.View ArticleGoogle Scholar
- Stoll N, Schmidt E, Thurow K: Isotope pattern evaluation for the reduction of elemental compositions assigned to high-resolution mass spectral data from electrospray ionisation Fourier transform ion cyclotron resonance mass spectrometry. J Am Soc Mass Spectrom 2006, 17: 1692–1699.View ArticlePubMedGoogle Scholar
- Tautenhahn R, Bottcher C, Neumann S: Annotation of LC/ESI-MS mass signals. In Lecture Notes in Computer Science, Bioinformatics Research and Development. Volume 4414. Springer Berlin/Heidelberg; 2007:371–380.View ArticleGoogle Scholar
- Iijima Y, Nakamura Y, Ogata Y, Tanaka K, Sakurai N, Suzuki T, Suzuki H, Okazaki K, Kitayama M, Kanaya S, Aoki K, Shibata D: Metabolite annotations based on the integration of mass spectral information. Plant J 2008, 54: 949–962.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakamura Y, Kanaya S, Sakurai N, Iijima Y, Aoki K, Okazaki K, Suzuki H, Kitayama M, Shibata D: A tool for high-throughput prediction of molecular formulas and identification of isotopic peaks from large scale mass spectrometry data. Plant Biotech 2008, 25: 377–380.View ArticleGoogle Scholar
- Takahashi H, Kai K, Shinbo Y, Tanaka K, Ohta D, Oshima T, Altaf-Ul-Amin M, Kurokawa K, Ogasawara N, Kanaya S: Metabolomics approach for determining growth-specific metabolites based on fourier transform ion cyclotron resonance mass spectrometry. Anal Bioanal Chem 2008, 391: 2769–2782.PubMed CentralView ArticlePubMedGoogle Scholar
- Enot DP, Lin W, Beckmann M, Parker D, Overy DP, Draper J: Pre-processing, classification modelling and feature selection using Flow Injection Electrospray Mass Spectrometry (FIE-MS) metabolome fingerprint data. Nat Protocols 2008, 3: 446–470.View ArticlePubMedGoogle Scholar
- Gipson GT, Tatsuoka KS, Sokhansanj BA, Ball RJ, Connor SC: Assignment of MS-based metabolomic datasets via compound interaction pair mapping. Metabolomics 2008, 4: 94–103.View ArticleGoogle Scholar
- Herrgård MJ, Swainston N, Dobson P, et al.: A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. Nat Biotech 2008, 26: 1155–1160.View ArticleGoogle Scholar
- Sud M, Fahy E, Cotter D, Brown HA, Dennis EA, Glass CK, Merrill AH Jr, Murphy RC, Raetz RC, Russell DW, Subramanian S: LMSD: LIPID MAPS structure database. Nuc Acids Res 2007, 35: 527–532.View ArticleGoogle Scholar
- Moco S, Bino RJ, Vorst O, Verhoeven HA, de Groot J, van Beek TA, Vervoort J, de Vos CHS: A liquid chromatography-mass spectrometry-based metabolome database for tomato. Plant Physiol 2006, 141: 1205–1218.PubMed CentralView ArticlePubMedGoogle Scholar
- Grafahrend-Belau E, Stephan S, Koschützki D, Scholz U, Junker BH, Schreiber F: MetaCrop: a detailed database of crop plant metabolism. Nuc Acids Res 2008, 36: 954–958.View ArticleGoogle Scholar
- Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, Paley S, Rhee SY, Shearer AG, Tissier C, Walk C, Zhang P, Karp PD: The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases. Nuc Acids Res 2008, 36: 623–631.View ArticleGoogle Scholar
- Wishart DS, Knox C, Guo AC, et al.: HMDB: a knowledgebase for the human metabolome. Nucleic Acids Res 2009, 37: 603–610.View ArticleGoogle Scholar
- Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, et al.: Database resources of the National Centre for Biotechnology Informatio. Nuc Acids Res 2007, 35: 5–12.View ArticleGoogle Scholar
- Ionisation 'rules'[http://maltese.dbs.aber.ac.uk:8888/hrmet/search/disprules.php]
- Hau J: Formula calculator (elemental composition calculator) HiRes MS vs. 20050617.[http://homepage.sunrise.ch/mysunrise/joerg.hau/sci/index.html]
- Kind T, Fiehn O: Seven Golden Rules for heuristic filtering of molecular formulas obtained by accurate mass spectrometry. BMC Bioinformatics 2007, 8: 105.PubMed CentralView ArticlePubMedGoogle Scholar
- Parker D, Beckmann M, Enot D, Overy D, Caracuel-Rios Z, Gilbert M, Talbot N, Draper J: Rice blast infection of Brachypodium distachyon as a model system to study dynamic host/pathogen interactions. Nat Protocols 2008, 3: 435–445.View ArticlePubMedGoogle Scholar
- Southam AD, Payne TG, Cooper HJ, Arvanitis TN, Viant MR: Dynamic range and mass accuracy of wide-scan direct infusion nanoelectrospray fourier transform ion cyclotron resonance mass spectrometry-based metabolomics increased by the spectral stitching method. Anal Chem 2007, 79: 4595–4602.View ArticlePubMedGoogle Scholar
- Mass mathematical relationships[http://maltese.dbs.aber.ac.uk:8888/hrmet/supp/askMZedDBworkflow.r]
- Isotope calculator[http://isotopatcalc.sourceforge.net]
- Seven Golden Rules[http://fiehnlab.ucdavis.edu/projects/Seven_Golden_Rules/]
- Potential Adducts list[http://maltese.dbs.aber.ac.uk:8888/hrmet/search/addsearch0.php]
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.