An analytic and systematic framework for estimating metabolic flux ratios from 13C tracer experiments

Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease

Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.

13C metabolic flux analysis: Classification and characterization from the perspective of mathematical modeling and application in physiological research of neural cell

¹³C metabolic flux analysis (¹³C-MFA) has emerged as a forceful tool for quantifying in vivo metabolic pathway activity of different biological systems. This technology plays an important role in understanding intracellular metabolism and revealing patho-physiology mechanism. Recently, it has evolved into a method family with great diversity in experiments, analytics, and mathematics. In this review, we classify and characterize the various branch of ¹³C-MFA from a unified perspective of mathematical modeling. By linking different parts in the model to each step of its workflow, the specific technologies of ¹³C-MFA are put into discussion, including the isotope labeling model (ILM), isotope pattern measuring technique, optimization algorithm and statistical method. Its application in physiological research in neural cell has also been reviewed.

Effects of Acute Subdural Hematoma-Induced Brain Injury On Energy Metabolism in Peripheral Blood Mononuclear Cells

ABSTRACT

In activated immune cells, differentiation and function are determined by cell type-specific modifications of metabolic patterns. After traumatic brain injury both immune cell activation and suppression were reported. Therefore, we sought to explore immune cell energy metabolism in a long-term, resuscitated porcine model of acute subdural hematoma (ASDH)-induced acute brain injury devoid of impaired systemic hemodynamics and oxygen transport.Before and up to 50 h after induction of ASDH, peripheral blood mononuclear cells (PBMCs) were separated by density gradient centrifugation, and cell metabolism was analyzed using high-resolution respirometry for mitochondrial respiration and electron spin resonance for reactive oxygen species production. After incubation with stable isotope-labeled 1,2-13C2-glucose or 13C5-glutamine, distinct labeling patterns of intermediates of glycolysis or tricarboxylic acid (TCA) cycle and 13CO2 production were measured by gas chromatography-mass spectroscopy. Principal component analysis was followed by a varimax rotation on the covariance across all measured variables and all measured time points.After ASDH induction, average PBMC metabolic activity remained unaffected, possibly because strict adherence to intensive care unit guidelines limited trauma to ASDH induction without any change in parameters of systemic hemodynamics, oxygen transport, and whole-body metabolism. Despite decreased glycolytic activity fueling the TCA cycle, the principal component analysis indicated a cell type-specific activation pattern with biosynthetic and proliferative characteristics.

Workflow Based on the Combination of Isotopic Tracer Experiments to Investigate Microbial Metabolism of Multiple Nutrient Sources

Studies in the field of microbiology rely on the implementation of a wide range of methodologies. In particular, the development of appropriate methods substantially contributes to providing extensive knowledge of the metabolism of microorganisms growing in chemically defined media containing unique nitrogen and carbon sources. In contrast, the management through metabolism of multiple nutrient sources, despite their broad presence in natural or industrial environments, remains virtually unexplored. This situation is mainly due to the lack of suitable methodologies, which hinders investigations. We report an experimental strategy to quantitatively and comprehensively explore how metabolism operates when a nutrient is provided as a mixture of different molecules, i.e., a complex resource. Here, we describe its application for assessing the partitioning of multiple nitrogen sources through the yeast metabolic network. The workflow combines information obtained during stable isotope tracer experiments using selected ¹³C- or ¹⁵N-labeled substrates. It first consists of parallel and reproducible fermentations in the same medium, which includes a mixture of N-containing molecules; however,a selected nitrogen source is labeled each time. A combination of analytical procedures (HPLC, GC-MS) is implemented to assess the labeling patterns of targeted compounds and to quantify the consumption and recovery of substrates in other metabolites. An integrated analysis of the complete dataset provides an overview of the fate of consumed substrates within cells. This approach requires an accurate protocol for the collection of samples-facilitated by a robot-assisted system for online monitoring of fermentations-and the achievement of numerous time-consuming analyses. Despite these constraints, it allowed understanding, for the first time, the partitioning of multiple nitrogen sources throughout the yeast metabolic network. We elucidated the redistribution of nitrogen from more abundant sources toward other N-compounds and determined the metabolic origins of volatile molecules and proteinogenic amino acids.

Estimation of flux ratios without uptake or release data: Application to serine and methionine metabolism

Model-based metabolic flux analysis (MFA) using isotope-labeled substrates has provided great insight into intracellular metabolic activities across a host of organisms. One challenge with applying MFA in mammalian systems, however, is the need for absolute quantification of nutrient uptake, biomass composition, and byproduct release fluxes. Such measurements are often not feasible in complex culture systems or in vivo. One way to address this issue is to estimate flux ratios, the fractional contribution of a flux to a metabolite pool, which are independent of absolute measurements and yet informative for cellular metabolism. Prior work has focused on "local" estimation of a handful of flux ratios for specific metabolites and reactions. Here, we perform systematic, model-based estimation of all flux ratios in a metabolic network using isotope labeling data, in the absence of uptake/release data. In a series of examples, we investigate what flux ratios can be well estimated with reasonably tight confidence intervals, and contrast this with confidence intervals on normalized fluxes. We find that flux ratios can provide useful information on the metabolic state, and is complementary to normalized fluxes: for certain metabolic reactions, only flux ratios can be well estimated, while for others normalized fluxes can be obtained. Simulation studies of a large human metabolic network model suggest that estimation of flux ratios is technically feasible for complex networks, but additional studies on data from actual isotopomer labeling experiments are needed to validate these results. Finally, we experimentally study serine and methionine metabolism in cancer cells using flux ratios. We find that, in these cells, the methionine cycle is truncated with little remethylation from homocysteine, and polyamine synthesis in the absence of methionine salvage leads to loss of 5-methylthioadenosine, suggesting a new mode of overflow metabolism in cancer cells. This work highlights the potential for flux ratio analysis in the absence of absolute quantification, which we anticipate will be important for both in vitro and in vivo studies of cancer metabolism.

Recent advances in high-throughput 13C-fluxomics

The rise of high throughput (HT) strain engineering tools accompanying the area of synthetic biology is supporting the generation of a large number of microbial cell factories. A current bottleneck in process development is our limited capacity to rapidly analyze the metabolic state of the engineered strains, and in particular their intracellular fluxes. HT ¹³C-fluxomics workflows have not yet become commonplace, despite the existence of several HT tools at each of the required stages. This includes cultivation and sampling systems, analytics for isotopic analysis, and software for data processing and flux calculation. Here, we review recent advances in the field and highlight bottlenecks that must be overcome to allow the emergence of true HT ¹³C-fluxomics workflows.

SUMOFLUX: A Generalized Method for Targeted 13C Metabolic Flux Ratio Analysis

Metabolic fluxes are a cornerstone of cellular physiology that emerge from a complex interplay of enzymes, carriers, and nutrients. The experimental assessment of in vivo intracellular fluxes using stable isotopic tracers is essential if we are to understand metabolic function and regulation. Flux estimation based on 13C or 2H labeling relies on complex simulation and iterative fitting; processes that necessitate a level of expertise that ordinarily preclude the non-expert user. To overcome this, we have developed SUMOFLUX, a methodology that is broadly applicable to the targeted analysis of 13C-metabolic fluxes. By combining surrogate modeling and machine learning, we trained a predictor to specialize in estimating flux ratios from measurable 13C-data. SUMOFLUX targets specific flux features individually, which makes it fast, user-friendly, applicable to experimental design and robust in terms of experimental noise and exchange flux magnitude. Collectively, we predict that SUMOFLUX's properties realistically pave the way to high-throughput flux analyses.

Natural isotope correction of MS/MS measurements for metabolomics and (13)C fluxomics

Fluxomics and metabolomics are crucial tools for metabolic engineering and biomedical analysis to determine the in vivo cellular state. Especially, the application of (13)C isotopes allows comprehensive insights into the functional operation of cellular metabolism. Compared to single MS, tandem mass spectrometry (MS/MS) provides more detailed and accurate measurements of the metabolite enrichment patterns (tandem mass isotopomers), increasing the accuracy of metabolite concentration measurements and metabolic flux estimation. MS-type data from isotope labeling experiments is biased by naturally occurring stable isotopes (C, H, N, O, etc.). In particular, GC-MS(/MS) requires derivatization for the usually non-volatile intracellular metabolites introducing additional natural isotopes leading to measurements that do not directly represent the carbon labeling distribution. To make full use of LC- and GC-MS/MS mass isotopomer measurements, the influence of natural isotopes has to be eliminated (corrected). Our correction approach is analyzed for the two most common applications; (13)C fluxomics and isotope dilution mass spectrometry (IDMS) based metabolomics. Natural isotopes can have an impact on the calculated flux distribution which strongly depends on the substrate labeling and the actual flux distribution. Second, we show that in IDMS based metabolomics natural isotopes lead to underestimated concentrations that can and should be corrected with a nonlinear calibration. Our simulations indicate that the correction for natural abundance in isotope based fluxomics and quantitative metabolomics is essential for correct data interpretation.

Efficient Modeling of MS/MS Data for Metabolic Flux Analysis

Metabolic flux analysis (MFA) is a widely used method for quantifying intracellular metabolic fluxes. It works by feeding cells with isotopic labeled nutrients, measuring metabolite isotopic labeling, and computationally interpreting the measured labeling data to estimate flux. Tandem mass-spectrometry (MS/MS) has been shown to be useful for MFA, providing positional isotopic labeling data. Specifically, MS/MS enables the measurement of a metabolite tandem mass-isotopomer distribution, representing the abundance in which certain parent and product fragments of a metabolite have different number of labeled atoms. However, a major limitation in using MFA with MS/MS data is the lack of a computationally efficient method for simulating such isotopic labeling data. Here, we describe the tandemer approach for efficiently computing metabolite tandem mass-isotopomer distributions in a metabolic network, given an estimation of metabolic fluxes. This approach can be used by MFA to find optimal metabolic fluxes, whose induced metabolite labeling patterns match tandem mass-isotopomer distributions measured by MS/MS. The tandemer approach is applied to simulate MS/MS data in a small-scale metabolic network model of mammalian methionine metabolism and in a large-scale metabolic network model of E. coli. It is shown to significantly improve the running time by between two to three orders of magnitude compared to the state-of-the-art, cumomers approach. We expect the tandemer approach to promote broader usage of MS/MS technology in metabolic flux analysis. Implementation is freely available at www.cs.technion.ac.il/~tomersh/methods.html

Use of (13)C stable isotope labelling for pathway and metabolic flux analysis in Leishmania parasites

This protocol describes the combined use of metabolite profiling and stable isotope labelling to define pathways of central carbon metabolism in the protozoa parasite, Leishmania mexicana. Parasite stages are cultivated in standard or completely defined media and then rapidly transferred to chemically equivalent media containing a single (13)C-labelled nutrient. The incorporation of label can be followed over time or after establishment of isotopic equilibrium by harvesting parasites with rapid metabolic quenching. (13)C enrichment of multiple intracellular polar and apolar (lipidic) metabolites can be quantified using gas chromatography-mass spectrometry (GC-MS), while the uptake and secretion of (13)C-labelled metabolites can be measured by (13)C-NMR. Analysis of the mass isotopomer distribution of key metabolites provides information on pathway structure, while analysis of labelling kinetics can be used to infer metabolic fluxes. This protocol is exemplified using L. mexicana labelled with (13)C-U-glucose. The method can be used to measure perturbations in parasite metabolism induced by drug inhibition or genetic manipulation of enzyme levels and is broadly applicable to any cultured parasite stages.

Comparative genome-scale reconstruction of gapless metabolic networks for present and ancestral species

We introduce a novel computational approach, CoReCo, for comparative metabolic reconstruction and provide genome-scale metabolic network models for 49 important fungal species. Leveraging on the exponential growth in sequenced genome availability, our method reconstructs genome-scale gapless metabolic networks simultaneously for a large number of species by integrating sequence data in a probabilistic framework. High reconstruction accuracy is demonstrated by comparisons to the well-curated Saccharomyces cerevisiae consensus model and large-scale knock-out experiments. Our comparative approach is particularly useful in scenarios where the quality of available sequence data is lacking, and when reconstructing evolutionary distant species. Moreover, the reconstructed networks are fully carbon mapped, allowing their use in 13C flux analysis. We demonstrate the functionality and usability of the reconstructed fungal models with computational steady-state biomass production experiment, as these fungi include some of the most important production organisms in industrial biotechnology. In contrast to many existing reconstruction techniques, only minimal manual effort is required before the reconstructed models are usable in flux balance experiments. CoReCo is available at http://esaskar.github.io/CoReCo/.

Advances in next-generation sequencing technologies are revolutionizing molecular biology. Sequencing-enabled cost-effective characterization of microbial genomes is a particularly exciting development in metabolic engineering. There, considerable effort has been put to reconstructing genome-scale metabolic networks that describe the collection of hundreds to thousands of biochemical reactions available for a microbial cell. These network models are instrumental in understanding microbial metabolism and guiding metabolic engineering efforts to improve biochemical yields. We have developed a novel computational method, CoReCo, which bridges the growing gap between the availability of sequenced genomes and respective reconstructed metabolic networks. The method reconstructs genome-scale metabolic networks simultaneously for related microbial species. It utilizes the available sequencing data from these species to correct for incomplete and missing data. We used the method to reconstruct metabolic networks for a set of 49 fungal species providing the method protein sequence data and a phylogenetic tree describing the evolutionary relationships between the species. We demonstrate the applicability of the method by comparing a metabolic reconstruction of Saccharomyces cerevisiae to the manually curated, high-quality consensus network. We also provide an easy-to-use implementation of the method, usable both in single computer and distributed computing environments.

Labelling analysis for ¹³C MFA using NMR spectroscopy

NMR spectroscopy is an efficient method for analyzing (13)C labelling of cellular metabolites. The strength of it is especially the ability to provide direct quantitative positional information on the (13)C labelling status of carbon atoms in metabolites. NMR spectroscopic methods allow also for detection of contiguously (13)C-labelled fragments in the carbon backbones of the metabolites. Furthermore, the recent developments of NMR spectroscopy hardware have substantially improved the sensitivity of the methods. In this chapter we describe a method for analyzing the (13)C labelling of the biomass amino acids for metabolic flux analysis, sample preparation for NMR spectroscopy, acquiring and processing the NMR spectra, and extracting the (13)C labelling information from the NMR data. Different NMR methods are applied depending on the (13)C labelling strategy chosen. These strategies include uniform (13)C labelling, positional (13)C labelling, or a combination of both. Not only the preparation of sample for analysis of (13)C labelling in proteinogenic amino acids in biomass is described, but also the necessary modifications to the method when analysis of (13)C labelling in free metabolic intermediates is of interest. Finally the strategies for using the different NMR-detected (13)C labelling data in (13)C MFA are discussed.

Software LS-MIDA for efficient mass isotopomer distribution analysis in metabolic modelling

BACKGROUND

The knowledge of metabolic pathways and fluxes is important to understand the adaptation of organisms to their biotic and abiotic environment. The specific distribution of stable isotope labelled precursors into metabolic products can be taken as fingerprints of the metabolic events and dynamics through the metabolic networks. An open-source software is required that easily and rapidly calculates from mass spectra of labelled metabolites, derivatives and their fragments global isotope excess and isotopomer distribution.

RESULTS

The open-source software "Least Square Mass Isotopomer Analyzer" (LS-MIDA) is presented that processes experimental mass spectrometry (MS) data on the basis of metabolite information such as the number of atoms in the compound, mass to charge ratio (m/e or m/z) values of the compounds and fragments under study, and the experimental relative MS intensities reflecting the enrichments of isotopomers in 13C- or 15 N-labelled compounds, in comparison to the natural abundances in the unlabelled molecules. The software uses Brauman's least square method of linear regression. As a result, global isotope enrichments of the metabolite or fragment under study and the molar abundances of each isotopomer are obtained and displayed.

CONCLUSIONS

The new software provides an open-source platform that easily and rapidly converts experimental MS patterns of labelled metabolites into isotopomer enrichments that are the basis for subsequent observation-driven analysis of pathways and fluxes, as well as for model-driven metabolic flux calculations.

Low mass MS/MS fragments of protonated amino acids used for distinction of their 13C-isotopomers in metabolic studies

Glu, Gln, Pro, and Ala are the main amino acids involved in ammonia detoxification in mosquitoes. In order to develop a tandem mass spectrometry method (MS(2)) to monitor each carbon of the above isotopically-labeled (13)C-amino acids for metabolic studies, the compositions and origins of atoms in fragments of the protonated amino acid should be first elucidated. Thus, various electrospray (ESI)-based MS(2) tools were employed to study the fragmentation of these unlabeled and isotopically-labeled amino acids and better understand their dissociation pathways. A broad range of fragments, including previously-undescribed low m/z fragments was revealed. The formulae of the fragments (from m/z 130 down to m/z 27) were confirmed by their accurate masses. The structures and conformations of the larger fragments of Glu were also explored by ion mobility mass spectrometry (IM-MS) and gas-phase hydrogen/deuterium exchange (HDX) experiments. It was found that some low m/z fragments (m/z 27-30) are common to Glu, Gln, Pro, and Ala. The origins of carbons in these small fragments are discussed and additional collision induced dissociation (CID) MS(2) fragmentation pathways are proposed for them. It was also found that small fragments (≤m/z 84) of protonated, methylated Glu, and methylated Gln are the same as those of the underivatized Glu and Gln. Taken together, the new approach of utilizing low m/z fragments can be applied to distinguish, identify, and quantify (13)C-amino acids labeled at various positions, either in the backbone or side chain.

13C-based metabolic flux analysis: fundamentals and practice

Isotope-based metabolic flux analysis is one of the emerging technologies applied to system level metabolic phenotype characterization in metabolic engineering. Among the developed approaches, (13)C-based metabolic flux analysis has been established as a standard tool and has been widely applied to quantitative pathway characterization of diverse biological systems. To implement (13)C-based metabolic flux analysis in practice, comprehending the underlying mathematical and computational modeling fundamentals is of importance along with carefully conducted experiments and analytical measurements. Such knowledge is also crucial when designing (13)C-labeling experiments and properly acquiring key data sets essential for in vivo flux analysis implementation. In this regard, the modeling fundamentals of (13)C-labeling systems and analytical data processing are the main topics we will deal with in this chapter. Along with this, the relevant numerical optimization techniques are addressed to help implementation of the entire computational procedures aiming at (13)C-based metabolic flux analysis in vivo.

Flux-p: automating metabolic flux analysis

Quantitative knowledge of intracellular fluxes in metabolic networks is invaluable for inferring metabolic system behavior and the design principles of biological systems. However, intracellular reaction rates can not often be calculated directly but have to be estimated; for instance, via ¹³C-based metabolic flux analysis, a model-based interpretation of stable carbon isotope patterns in intermediates of metabolism. Existing software such as FiatFlux, OpenFLUX or 13CFLUX supports experts in this complex analysis, but requires several steps that have to be carried out manually, hence restricting the use of this software for data interpretation to a rather small number of experiments. In this paper, we present Flux-P as an approach to automate and standardize ¹³C-based metabolic flux analysis, using the Bio-jETI workflow framework. Exemplarily based on the FiatFlux software, it demonstrates how services can be created that carry out the different analysis steps autonomously and how these can subsequently be assembled into software workflows that perform automated, high-throughput intracellular flux analysis of high quality and reproducibility. Besides significant acceleration and standardization of the data analysis, the agile workflow-based realization supports flexible changes of the analysis workflows on the user level, making it easy to perform custom analyses.

Metabolic modelling in the development of cell factories by synthetic biology

Cell factories are commonly microbial organisms utilized for bioconversion of renewable resources to bulk or high value chemicals. Introduction of novel production pathways in chassis strains is the core of the development of cell factories by synthetic biology. Synthetic biology aims to create novel biological functions and systems not found in nature by combining biology with engineering. The workflow of the development of novel cell factories with synthetic biology is ideally linear which will be attainable with the quantitative engineering approach, high-quality predictive models, and libraries of well-characterized parts. Different types of metabolic models, mathematical representations of metabolism and its components, enzymes and metabolites, are useful in particular phases of the synthetic biology workflow. In this minireview, the role of metabolic modelling in synthetic biology will be discussed with a review of current status of compatible methods and models for the in silico design and quantitative evaluation of a cell factory.

Integrating quantitative proteomics and metabolomics with a genome-scale metabolic network model

Motivation: The availability of modern sequencing techniques has led to a rapid increase in the amount of reconstructed metabolic networks. Using these models as a platform for the analysis of high throughput transcriptomic, proteomic and metabolomic data can provide valuable insight into conditional changes in the metabolic activity of an organism. While transcriptomics and proteomics provide important insights into the hierarchical regulation of metabolic flux, metabolomics shed light on the actual enzyme activity through metabolic regulation and mass action effects. Here we introduce a new method, termed integrative omics-metabolic analysis (IOMA) that quantitatively integrates proteomic and metabolomic data with genome-scale metabolic models, to more accurately predict metabolic flux distributions. The method is formulated as a quadratic programming (QP) problem that seeks a steady-state flux distribution in which flux through reactions with measured proteomic and metabolomic data, is as consistent as possible with kinetically derived flux estimations.

Results: IOMA is shown to successfully predict the metabolic state of human erythrocytes (compared to kinetic model simulations), showing a significant advantage over the commonly used methods flux balance analysis and minimization of metabolic adjustment. Thereafter, IOMA is shown to correctly predict metabolic fluxes in Escherichia coli under different gene knockouts for which both metabolomic and proteomic data is available, achieving higher prediction accuracy over the extant methods. Considering the lack of high-throughput flux measurements, while high-throughput metabolomic and proteomic data are becoming readily available, we expect IOMA to significantly contribute to future research of cellular metabolism.

Contacts:kerenyiz@post.tau.ac.il; tomersh@cs.technion.ac.il

13C metabolic flux analysis in complex systems

Experimental determination of in vivo metabolic rates by methods of (13)C metabolic flux analysis is a pivotal approach to unravel structure and regulation of metabolic networks, in particular with microorganisms grown in minimal media. However, the study of real-life and eukaryotic systems calls for the quantification of fluxes also in cellular compartments, rich media, cell-wide metabolic networks, dynamic systems or single cells. These scenarios drastically increase the complexity of the task, which is only partly dealt by existing approaches that rely on rigorous simulations of label propagation through metabolic networks and require multiple labeling experiments or a priori information on pathway inactivity to simplify the problem. Albeit qualitative and largely driven by human interpretation, statistical analysis of measured (13)C-patterns remains the exclusive alternative to comprehensively handle such complex systems. In the future, this practice will be complemented by novel modeling frameworks to assay particular fluxes within a network by stable isotopic tracer for targeted validation of well-defined hypotheses.

Systems biology of recombinant protein production in Bacillus megaterium

Over the last two decades the Gram-positive bacterium Bacillus megaterium was systematically developed to a useful alternative protein production host. Multiple vector systems for high yield intra- and extracellular protein production were constructed. Strong inducible promoters were combined with DNA sequences for optimised ribosome binding sites, various leader peptides for protein export and N- as well as C-terminal affinity tags for affinity chromatographic purification of the desired protein. High cell density cultivation and recombinant protein production were successfully tested. For further system biology based control and optimisation of the production process the genomes of two B. megaterium strains were completely elucidated, DNA arrays designed, proteome, fluxome and metabolome analyses performed and all data integrated using the bioinformatics platform MEGABAC. Now, solid theoretical and experimental bases for primary modeling attempts of the production process are available.

Inferring branching pathways in genome-scale metabolic networks

Background

A central problem in computational metabolic modelling is how to find biochemically plausible pathways between metabolites in a metabolic network. Two general, complementary frameworks have been utilized to find metabolic pathways: constraint-based modelling and graph-theoretical path finding approaches. In constraint-based modelling, one aims to find pathways where metabolites are balanced in a pseudo steady-state. Constraint-based methods, such as elementary flux mode analysis, have typically a high computational cost stemming from a large number of steady-state pathways in a typical metabolic network. On the other hand, graph-theoretical approaches avoid the computational complexity of constraint-based methods by solving a simpler problem of finding shortest paths. However, while scaling well with network size, graph-theoretic methods generally tend to return more false positive pathways than constraint-based methods.

Results

In this paper, we introduce a computational method, ReTrace, for finding biochemically relevant, branching metabolic pathways in an atom-level representation of metabolic networks. The method finds compact pathways which transfer a high fraction of atoms from source to target metabolites by considering combinations of linear shortest paths. In contrast to current steady-state pathway analysis methods, our method scales up well and is able to operate on genome-scale models. Further, we show that the pathways produced are biochemically meaningful by an example involving the biosynthesis of inosine 5'-monophosphate (IMP). In particular, the method is able to avoid typical problems associated with graph-theoretic approaches such as the need to define side metabolites or pathways not carrying any net carbon flux appearing in results. Finally, we discuss an application involving reconstruction of amino acid pathways of a recently sequenced organism demonstrating how measurement data can be easily incorporated into ReTrace analysis. ReTrace is licensed under GPL and is freely available for academic use at http://www.cs.helsinki.fi/group/sysfys/software/retrace/.

Conclusion

ReTrace is a useful method in metabolic path finding tasks, combining some of the best aspects in constraint-based and graph-theoretic methods. It finds use in a multitude of tasks ranging from metabolic engineering to metabolic reconstruction of recently sequenced organisms.

(13)C-based metabolic flux analysis

Stable isotope, and in particular (13)C-based flux analysis, is the exclusive approach to experimentally quantify the integrated responses of metabolic networks. Here we describe a protocol that is based on growing microbes on (13)C-labeled glucose and subsequent gas chromatography mass spectrometric detection of (13)C-patterns in protein-bound amino acids. Relying on publicly available software packages, we then describe two complementary mathematical approaches to estimate either local ratios of converging fluxes or absolute fluxes through different pathways. As amino acids in cell protein are abundant and stable, this protocol requires a minimum of equipment and analytical expertise. Most other flux methods are variants of the principles presented here. A true alternative is the analytically more demanding dynamic flux analysis that relies on (13)C-pattern in free intracellular metabolites. The presented protocols take 5-10 d, have been used extensively in the past decade and are exemplified here for the central metabolism of Escherichia coli.

OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis

Background

The quantitative analysis of metabolic fluxes, i.e., in vivo activities of intracellular enzymes and pathways, provides key information on biological systems in systems biology and metabolic engineering. It is based on a comprehensive approach combining (i) tracer cultivation on ¹³C substrates, (ii) ¹³C labelling analysis by mass spectrometry and (iii) mathematical modelling for experimental design, data processing, flux calculation and statistics. Whereas the cultivation and the analytical part is fairly advanced, a lack of appropriate modelling software solutions for all modelling aspects in flux studies is limiting the application of metabolic flux analysis.

Results

We have developed OpenFLUX as a user friendly, yet flexible software application for small and large scale ¹³C metabolic flux analysis. The application is based on the new Elementary Metabolite Unit (EMU) framework, significantly enhancing computation speed for flux calculation. From simple notation of metabolic reaction networks defined in a spreadsheet, the OpenFLUX parser automatically generates MATLAB-readable metabolite and isotopomer balances, thus strongly facilitating model creation. The model can be used to perform experimental design, parameter estimation and sensitivity analysis either using the built-in gradient-based search or Monte Carlo algorithms or in user-defined algorithms. Exemplified for a microbial flux study with 71 reactions, 8 free flux parameters and mass isotopomer distribution of 10 metabolites, OpenFLUX allowed to automatically compile the EMU-based model from an Excel file containing metabolic reactions and carbon transfer mechanisms, showing it's user-friendliness. It reliably reproduced the published data and optimum flux distributions for the network under study were found quickly (<20 sec).

Conclusion

We have developed a fast, accurate application to perform steady-state ¹³C metabolic flux analysis. OpenFLUX will strongly facilitate and enhance the design, calculation and interpretation of metabolic flux studies. By providing the software open source, we hope it will evolve with the rapidly growing field of fluxomics.

FiD: a software for ab initio structural identification of product ions from tandem mass spectrometric data

We present FiD (Fragment iDentificator), a software tool for the structural identification of product ions produced with tandem mass spectrometric measurement of low molecular weight organic compounds. Tandem mass spectrometry (MS/MS) has proven to be an indispensable tool in modern, cell-wide metabolomics and fluxomics studies. In such studies, the structural information of the MS(n) product ions is usually needed in the downstream analysis of the measurement data. The manual identification of the structures of MS(n) product ions is, however, a nontrivial task requiring expertise, and calls for computer assistance. Commercial software tools, such as Mass Frontier and ACD/MS Fragmenter, rely on fragmentation rule databases for the identification of MS(n) product ions. FiD, on the other hand, conducts a combinatorial search over all possible fragmentation paths and outputs a ranked list of alternative structures. This gives the user an advantage in situations where the MS/MS data of compounds with less well-known fragmentation mechanisms are processed. FiD software implements two fragmentation models, the single-step model that ignores intermediate fragmentation states and the multi-step model, which allows for complex fragmentation pathways. The software works for MS/MS data produced both in positive- and negative-ion modes. The software has an easy-to-use graphical interface with built-in visualization capabilities for structures of product ions and fragmentation pathways. In our experiments involving amino acids and sugar-phosphates, often found, e.g., in the central carbon metabolism of yeasts, FiD software correctly predicted the structures of product ions on average in 85% of the cases. The FiD software is free for academic use and is available for download from www.cs.helsinki.fi/group/sysfys/software/fragid.