Probabilistic assignment of formulas to mass peaks in metabolomics experiments

<b>Motivation</b>: High-accuracy mass spectrometry is a popular technology for high-throughput measurements of cellular metabolites (metabolomics). One of the major challenges is the correct identification of the observed mass peaks, including the assignment of their empirical formula, based on the measured mass.<p></p> <b>Results</b>: We propose a novel probabilistic method for the assignment of empirical formulas to mass peaks in high-throughput metabolomics mass spectrometry measurements. The method incorporates information about possible biochemical transformations between the empirical formulas to assign higher probability to formulas that could be created from other metabolites in the sample. In a series of experiments, we show that the method performs well and provides greater insight than assignments based on mass alone. In addition, we extend the model to incorporate isotope information to achieve even more reliable formula identification.<p></p&gt

Metabolomics to unveil and understand phenotypic diversity between pathogen populations

Visceral leishmaniasis is caused by a parasite called Leishmania donovani, which every year infects about half a million people and claims several thousand lives. Existing treatments are now becoming less effective due to the emergence of drug resistance. Improving our understanding of the mechanisms used by the parasite to adapt to drugs and achieve resistance is crucial for developing future treatment strategies. Unfortunately, the biological mechanism whereby Leishmania acquires drug resistance is poorly understood. Recent years have brought new technologies with the potential to increase greatly our understanding of drug resistance mechanisms. The latest mass spectrometry techniques allow the metabolome of parasites to be studied rapidly and in great detail. We have applied this approach to determine the metabolome of drug-sensitive and drug-resistant parasites isolated from patients with leishmaniasis. The data show that there are wholesale differences between the isolates and that the membrane composition has been drastically modified in drug-resistant parasites compared with drug-sensitive parasites. Our findings demonstrate that untargeted metabolomics has great potential to identify major metabolic differences between closely related parasite strains and thus should find many applications in distinguishing parasite phenotypes of clinical relevance

MetaboSearch: Tool for Mass-Based Metabolite Identification Using Multiple Databases

Searching metabolites against databases according to their masses is often the first step in metabolite identification for a mass spectrometry-based untargeted metabolomics study. Major metabolite databases include Human Metabolome DataBase (HMDB), Madison Metabolomics Consortium Database (MMCD), Metlin, and LIPID MAPS. Since each one of these databases covers only a fraction of the metabolome, integration of the search results from these databases is expected to yield a more comprehensive coverage. However, the manual combination of multiple search results is generally difficult when identification of hundreds of metabolites is desired. We have implemented a web-based software tool that enables simultaneous mass-based search against the four major databases, and the integration of the results. In addition, more complete chemical identifier information for the metabolites is retrieved by cross-referencing multiple databases. The search results are merged based on IUPAC International Chemical Identifier (InChI) keys. Besides a simple list of m/z values, the software can accept the ion annotation information as input for enhanced metabolite identification. The performance of the software is demonstrated on mass spectrometry data acquired in both positive and negative ionization modes. Compared with search results from individual databases, MetaboSearch provides better coverage of the metabolome and more complete chemical identifier information. Availability: The software tool is available at http://omics.georgetown.edu/MetaboSearch.html

Ultra-fast searching assists in evaluating sub-ppm mass accuracy enhancement in U-HPLC/Orbitrap MS data

A strategy, detailed methodology description and software are given with which the mass accuracy of U-HPLC-Orbitrap data (resolving power 50,000 FWHM) can be enhanced by an order of magnitude to sub-ppm levels. After mass accuracy enhancement all 211 reference masses have mass errors within 0.5 ppm; only 14 of these are outside the 0.2 ppm error margin. Further demonstration of mass accuracy enhancement is shown on a pre-concentrated urine sample in which evidence for 89 (342 ions) potential hydroxylated and glucuronated DHEA-metabolites is found. Although most DHEA metabolites have low-intensity mass signals, only 11 out of 342 are outside the ±1 ppm error envelop; 272 mass signals have errors below 0.5 ppm (142 below 0.2 ppm). The methodology consists of: (a) a multiple internal lock correction (here ten masses; no identity of internal lock masses is required) to avoid suppression problems of a single internal lock mass as well as to increase lock precision, (b) a multiple external mass correction (here 211 masses) to correct for calibration errors, (c) intensity dependant mass correction, (d) file averaging. The strategy is supported by ultra-fast file searching of baseline corrected, noise-reduced metAlign output. The output and efficiency of ultra-fast searching is essential in obtaining the required information to visualize the distribution of mass errors and isotope ratio deviations as a function of mass and intensity

Phosphoenolpyruvate carboxylase dentified as a key enzyme in erythrocytic Plasmodium falciparum carbon metabolism

Phospoenolpyruvate carboxylase (PEPC) is absent from humans but encoded in thePlasmodium falciparum genome, suggesting that PEPC has a parasite-specific function. To investigate its importance in P. falciparum, we generated a pepc null mutant (D10Δpepc), which was only achievable when malate, a reduction product of oxaloacetate, was added to the growth medium. D10Δpepc had a severe growth defect in vitro, which was partially reversed by addition of malate or fumarate, suggesting that pepc may be essential in vivo. Targeted metabolomics using 13C-U-D-glucose and 13C-bicarbonate showed that the conversion of glycolytically-derived PEP into malate, fumarate, aspartate and citrate was abolished in D10Δpepc and that pentose phosphate pathway metabolites and glycerol 3-phosphate were present at increased levels. In contrast, metabolism of the carbon skeleton of 13C,15N-U-glutamine was similar in both parasite lines, although the flux was lower in D10Δpepc; it also confirmed the operation of a complete forward TCA cycle in the wild type parasite. Overall, these data confirm the CO2 fixing activity of PEPC and suggest that it provides metabolites essential for TCA cycle anaplerosis and the maintenance of cytosolic and mitochondrial redox balance. Moreover, these findings imply that PEPC may be an exploitable target for future drug discovery

Benznidazole biotransformation and multiple targets in <i>Trypanosoma</i> cruzi revealed by metabolomics

&lt;b&gt;Background&lt;/b&gt;&lt;p&gt;&lt;/p&gt; The first line treatment for Chagas disease, a neglected tropical disease caused by the protozoan parasite Trypanosoma cruzi, involves administration of benznidazole (Bzn). Bzn is a 2-nitroimidazole pro-drug which requires nitroreduction to become active, although its mode of action is not fully understood. In the present work we used a non-targeted MS-based metabolomics approach to study the metabolic response of T. cruzi to Bzn.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Methodology/Principal findings&lt;/b&gt;&lt;p&gt;&lt;/p&gt; Parasites treated with Bzn were minimally altered compared to untreated trypanosomes, although the redox active thiols trypanothione, homotrypanothione and cysteine were significantly diminished in abundance post-treatment. In addition, multiple Bzn-derived metabolites were detected after treatment. These metabolites included reduction products, fragments and covalent adducts of reduced Bzn linked to each of the major low molecular weight thiols: trypanothione, glutathione, γ-glutamylcysteine, glutathionylspermidine, cysteine and ovothiol A. Bzn products known to be generated in vitro by the unusual trypanosomal nitroreductase, TcNTRI, were found within the parasites, but low molecular weight adducts of glyoxal, a proposed toxic end-product of NTRI Bzn metabolism, were not detected.&lt;p&gt;&lt;/p&gt; &lt;b&gt;Conclusions/significance&lt;/b&gt;&lt;p&gt;&lt;/p&gt; Our data is indicative of a major role of the thiol binding capacity of Bzn reduction products in the mechanism of Bzn toxicity against T. cruzi

Assessment of Metabolome Annotation Quality: A Method for Evaluating the False Discovery Rate of Elemental Composition Searches

BACKGROUND: In metabolomics researches using mass spectrometry (MS), systematic searching of high-resolution mass data against compound databases is often the first step of metabolite annotation to determine elemental compositions possessing similar theoretical mass numbers. However, incorrect hits derived from errors in mass analyses will be included in the results of elemental composition searches. To assess the quality of peak annotation information, a novel methodology for false discovery rates (FDR) evaluation is presented in this study. Based on the FDR analyses, several aspects of an elemental composition search, including setting a threshold, estimating FDR, and the types of elemental composition databases most reliable for searching are discussed. METHODOLOGY/PRINCIPAL FINDINGS: The FDR can be determined from one measured value (i.e., the hit rate for search queries) and four parameters determined by Monte Carlo simulation. The results indicate that relatively high FDR values (30-50%) were obtained when searching time-of-flight (TOF)/MS data using the KNApSAcK and KEGG databases. In addition, searches against large all-in-one databases (e.g., PubChem) always produced unacceptable results (FDR >70%). The estimated FDRs suggest that the quality of search results can be improved not only by performing more accurate mass analysis but also by modifying the properties of the compound database. A theoretical analysis indicates that FDR could be improved by using compound database with smaller but higher completeness entries. CONCLUSIONS/SIGNIFICANCE: High accuracy mass analysis, such as Fourier transform (FT)-MS, is needed for reliable annotation (FDR <10%). In addition, a small, customized compound database is preferable for high-quality annotation of metabolome data

In the Law We Trust; Some Thoughts on the ‘Legislative Gap’ in Legal Studies

The Legitimacy and Effectiveness of Law & Governance in a World of Multilevel Jurisdiction

Using a Non-Image-Based Medium-Throughput Assay for Screening Compounds Targeting N-myristoylation in Intracellular Leishmania Amastigotes

We have refined a medium-throughput assay to screen hit compounds for activity against N-myristoylation in intracellular amastigotes of Leishmania donovani. Using clinically-relevant stages of wild type parasites and an Alamar blue-based detection method, parasite survival following drug treatment of infected macrophages is monitored after macrophage lysis and transformation of freed amastigotes into replicative extracellular promastigotes. The latter transformation step is essential to amplify the signal for determination of parasite burden, a factor dependent on equivalent proliferation rate between samples. Validation of the assay has been achieved using the anti-leishmanial gold standard drugs, amphotericin B and miltefosine, with EC50 values correlating well with published values. This assay has been used, in parallel with enzyme activity data and direct assay on isolated extracellular amastigotes, to test lead-like and hit-like inhibitors of Leishmania Nmyristoyl transferase (NMT). These were derived both from validated in vivo inhibitors of Trypanosoma brucei NMT and a recent high-throughput screen against L. donovani NMT. Despite being a potent inhibitor of L. donovani NMT, the activity of the lead T. brucei NMT inhibitor (DDD85646) against L. donovani amastigotes is relatively poor. Encouragingly, analogues of DDD85646 show improved translation of enzyme to cellular activity. In testing the high-throughput L. donovani hits, we observed macrophage cytotoxicity with compounds from two of the four NMT-selective series identified, while all four series displayed low enzyme to cellular translation, also seen here with the T. brucei NMT inhibitors. Improvements in potency and physicochemical properties will be required to deliver attractive lead-like Leishmania NMT inhibitors

LC–MS-based absolute metabolite quantification:Application to metabolic flux measurement in trypanosomes

Human African trypanosomiasis is a neglected tropical disease caused by the protozoan parasite, Trypanosoma brucei. In the mammalian bloodstream, the trypanosome’s metabolism differs significantly from that of its host. For example, the parasite relies exclusively on glycolysis for energy source. Recently, computational and mathematical models of trypanosome metabolism have been generated to assist in understanding the parasite metabolism with the aim of facilitating drug development. Optimisation of these models requires quantitative information, including metabolite concentrations and/or metabolic fluxes that have been hitherto unavailable on a large scale. Here, we have implemented an LC–MS-based method that allows large scale quantification of metabolite levels by using U-13C-labelled E. coli extracts as internal standards. Known amounts of labelled E. coli extract were added into the parasite samples, as well as calibration standards, and used to obtain calibration curves enabling us to convert intensities into concentrations. This method allowed us to reliably quantify the changes of 43 intracellular metabolites and 32 extracellular metabolites in the medium over time. Based on the absolute quantification, we were able to compute consumption and production fluxes. These quantitative data can now be used to optimise computational models of parasite metabolism