The yeast metabolome addressed by electrospray ionization mass spectrometry: Initiation of a mass spectral library and its applications for metabolic footprinting by direct infusion mass spectrometry

Mass spectrometry (MS) has been a major driver for metabolomics, and gas chromatography (GC)-MS has been one of the primary techniques used for microbial metabolomics. The use of liquid chromatography (LC)-MS has however been limited, but electrospray ionization (ESI) is very well suited for ionization of microbial metabolites without any previous derivatization needed. To address the capabilities of ESI-MS in detecting the metabolome of Saccharomyces cerevisiae, the in silico metabolome of this organism was used as a template to present a theoretical metabolome. This showed that in combination with the specificity of MS up to 84% of the metabolites can be identified in a high mass accuracy ESI-spectrum. A total of 66 metabolites were systematically analyzed by positive and negative ESI-MS/MS with the aim of initiating a spectral library for ESI of microbial metabolites. This systematic analysis gave insight into the ionization and fragmentation characteristics of the different metabolites. With this insight, a small study of metabolic footprinting with ESI-MS demonstrated that biological information can be extracted from footprinting spectra. Statistical analysis of the footprinting data revealed discriminating ions, which could be assigned using the in silico metabolome. By this approach metabolic footprinting can advance from a classification method that is used to derive biological information based on guilt-by-association, to a tool for extraction of metabolic differences, which can guide new targeted biological experiments. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mass spectrometry for metabolic flux analysis

Mass spectrometry in combination with tracer experiments based on 13C substrates can serve as a powerful tool for the modeling and analysis of intracellular fluxes and the investigation of biochemical networks. The theoretical background for the application of mass spectrometry to metabolic flux analysis is discussed. Mass spectrometry methods are especially useful to determine mass distribution of metabolites. Additional information gained from fragmentation of metabolites, e.g., by electron impact ionization, allows further localization of labeling positions, up to complete resolution of isotopomer pools. To effectively handle mass distributions in simulation experiments, a matrix based general methodology is formulated. The natural isotope distribution of carbon, oxygen, hydrogen and nitrogen in the target metabolites is considered by introduction of correction matrices. It is shown by simulation results for the central carbon metabolism that neglecting natural isotope distributions causes significant errors in intracellular flux distributions. By varying relative fluxes into pentosephosphate pathway and pyruvate carboxylation reaction, marked changes in the mass distributions of metabolites result, which are illustrated for pyruvate, oxaloacetate, and alpha-ketoglutarate. In addition mass distributions of metabolites are significantly influenced over a broad range by the degree of reversibility of transaldolase and transketolase reactions in the pentosephosphate pathway. The mass distribution of metabolites is very sensitive towards intracellular flux patterns and can be measured with high accuracy by routine mass spectrometry methods. Copyright 1999 John Wiley & Sons, Inc.     

MapMaker and PathTracer for tracking carbon in genome‐scale metabolic models

Constraint‐based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome‐scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint‐based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO₂ fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint‐based models to identify feasible, active, and high flux paths between metabolites.     

Profiling of Arabidopsis secondary metabolites by capillary liquid chromatography coupled to electrospray ionization quadrupole time-of-flight mass spectrometry

Large-scale metabolic profiling is expected to develop into an integral part of functional genomics and systems biology. The metabolome of a cell or an organism is chemically highly complex. Therefore, comprehensive biochemical phenotyping requires a multitude of analytical techniques. Here, we describe a profiling approach that combines separation by capillary liquid chromatography with the high resolution, high sensitivity, and high mass accuracy of quadrupole time-of-flight mass spectrometry. About 2000 different mass signals can be detected in extracts of Arabidopsis roots and leaves. Many of these originate from Arabidopsis secondary metabolites. Detection based on retention times and exact masses is robust and reproducible. The dynamic range is sufficient for the quantification of metabolites. Assessment of the reproducibility of the analysis showed that biological variability exceeds technical variability. Tools were optimized or established for the automatic data deconvolution and data processing. Subtle differences between samples can be detected as tested with the chalcone synthase deficient tt4 mutant. The accuracy of time-of-flight mass analysis allows to calculate elemental compositions and to tentatively identify metabolites. In-source fragmentation and tandem mass spectrometry can be used to gain structural information. This approach has the potential to significantly contribute to establishing the metabolome of Arabidopsis and other model systems. The principles of separation and mass analysis of this technique, together with its sensitivity and resolving power, greatly expand the range of metabolic profiling.     

Expanding Metabolic Networks: Scopes of Compounds, Robustness, and Evolution

A new method for the mathematical analysis of large metabolic networks is presented. Based on the fact that the occurrence of a metabolic reaction generally requires the existence of other reactions providing its substrates, series of metabolic networks are constructed. In each step of the corresponding expansion process those reactions are incorporated whose substrates are made available by the networks of the previous generations. The method is applied to the set of all metabolic reactions included in the KEGG database. Starting with one or more seed compounds, the expansion results in a final network whose compounds define the scope of the seed. Scopes of all metabolic compounds are calculated and it is shown that large parts of cellular metabolism can be considered as the combined scope of simple building blocks. Analyses of various expansion processes reveal crucial metabolites whose incorporation allows for the increase in network complexity. Among these metabolites are common cofactors such as NAD⁺, ATP, and coenzyme A. We demonstrate that the outcome of network expansion is in general very robust against elimination of single or few reactions. There exist, however, crucial reactions whose elimination results in a dramatic reduction of scope sizes. It is hypothesized that the expansion process displays characteristics of the evolution of metabolism such as the temporal order of the emergence of metabolic pathways. [Reviewing Editor : Dr. David Pollock]     

Use of reconstituted metabolic networks to assist in metabolomic data visualization and mining

Metabolomics experiments seldom achieve their aim of comprehensively covering the entire metabolome. However, important information can be gleaned even from sparse datasets, which can be facilitated by placing the results within the context of known metabolic networks. Here we present a method that allows the automatic assignment of identified metabolites to positions within known metabolic networks, and, furthermore, allows automated extraction of sub-networks of biological significance. This latter feature is possible by use of a gap-filling algorithm. The utility of the algorithm in reconstructing and mining of metabolomics data is shown on two independent datasets generated with LC–MS LTQ-Orbitrap mass spectrometry. Biologically relevant metabolic sub-networks were extracted from both datasets. Moreover, a number of metabolites, whose presence eluded automatic selection within mass spectra, could be identified retrospectively by virtue of their inferred presence through gap filling.     

Stability of Metabolic Correlations under Changing Environmental Conditions in Escherichia coli – A Systems Approach

Background

Biological systems adapt to changing environments by reorganizing their cellular and physiological program with metabolites representing one important response level. Different stresses lead to both conserved and specific responses on the metabolite level which should be reflected in the underlying metabolic network.

Methodology/Principal Findings

Starting from experimental data obtained by a GC-MS based high-throughput metabolic profiling technology we here develop an approach that: (1) extracts network representations from metabolic condition-dependent data by using pairwise correlations, (2) determines the sets of stable and condition-dependent correlations based on a combination of statistical significance and homogeneity tests, and (3) can identify metabolites related to the stress response, which goes beyond simple observations about the changes of metabolic concentrations. The approach was tested with Escherichia coli as a model organism observed under four different environmental stress conditions (cold stress, heat stress, oxidative stress, lactose diauxie) and control unperturbed conditions. By constructing the stable network component, which displays a scale free topology and small-world characteristics, we demonstrated that: (1) metabolite hubs in this reconstructed correlation networks are significantly enriched for those contained in biochemical networks such as EcoCyc, (2) particular components of the stable network are enriched for functionally related biochemical pathways, and (3) independently of the response scale, based on their importance in the reorganization of the correlation network a set of metabolites can be identified which represent hypothetical candidates for adjusting to a stress-specific response.

Conclusions/Significance

Network-based tools allowed the identification of stress-dependent and general metabolic correlation networks. This correlation-network-based approach does not rely on major changes in concentration to identify metabolites important for stress adaptation, but rather on the changes in network properties with respect to metabolites. This should represent a useful complementary technique in addition to more classical approaches.     

细胞代谢复杂网络研究进展

复杂网络理论为细胞代谢网络研究提供了新的工具,基于复杂网络理论的细胞代谢网络研究可称细胞代谢复杂网络研究．先简要介绍了细胞代谢复杂网络的研究背景;随后详细总结和论述了细胞代谢复杂网络在建模、分析和控制三个方面的研究现状;再进一步指出了细胞代谢复杂网络在建模、分析和控制这三个方面研究中所存在的一些问题．为细胞代谢复杂网络领域的研究指出了一些有意义的方向,具有一定的参考价值。     

Whole‐cell Escherichia coli lactate biosensor for monitoring mammalian cell cultures during biopharmaceutical production

Many high‐value added recombinant proteins, such as therapeutic glycoproteins, are produced using mammalian cell cultures. In order to optimize the productivity of these cultures it is important to monitor cellular metabolism, for example the utilization of nutrients and the accumulation of metabolic waste products. One metabolic waste product of interest is lactic acid (lactate), overaccumulation of which can decrease cellular growth and protein production. Current methods for the detection of lactate are limited in terms of cost, sensitivity, and robustness. Therefore, we developed a whole‐cell Escherichia coli lactate biosensor based on the lldPRD operon and successfully used it to monitor lactate concentration in mammalian cell cultures. Using real samples and analytical validation we demonstrate that our biosensor can be used for absolute quantification of metabolites in complex samples with high accuracy, sensitivity, and robustness. Importantly, our whole‐cell biosensor was able to detect lactate at concentrations more than two orders of magnitude lower than the industry standard method, making it useful for monitoring lactate concentrations in early phase culture. Given the importance of lactate in a variety of both industrial and clinical contexts we anticipate that our whole‐cell biosensor can be used to address a range of interesting biological questions. It also serves as a blueprint for how to capitalize on the wealth of genetic operons for metabolite sensing available in nature for the development of other whole‐cell biosensors. Biotechnol. Bioeng. 2017;114: 1290–1300. © 2017 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.     

13C标记代谢流分析中天然稳定性同位素修正矩阵构建方法及应用

稳定性同位素¹³C标记实验是分析细胞代谢流的一种重要手段,主要通过质谱检测胞内代谢物中¹³C标记的同位素分布,并作为胞内代谢流计算时的约束条件,进而通过代谢流分析算法得到相应代谢网络中的通量分布。然而在自然界中,并非只有C元素存在天然稳定性同位素¹³C,其他元素如O元素也有其天然稳定性同位素¹⁷O、¹⁸O等,这使得质谱方法所测得的同位素分布中会夹杂除¹³C标记之外的其他元素的同位素信息,特别是分子中含有较多其他元素的分子,这将导致很大的实验误差,因此需要在进行代谢流计算前进行质谱数据的矫正。本研究提出了一种基于Python语言的天然同位素修正矩阵的构建方法,用于修正同位素分布测量值中由于天然同位素分布引起的测定误差。文中提出的基本修正矩阵幂方法用于构建各元素修正矩阵,结构简单、易于编码实现,可直接应用于¹³C代谢流分析软件数据前处理。将该修正方法应用于¹³C标记的黑曲霉（Aspergillus niger）胞内代谢流分析,结果表明本研究提出的方法准确有效,为准确获取微生物胞内代谢流分析提供了可靠的数据修正方法。     

Precision mapping of the metabolome

The global study of the structure and dynamics of metabolic networks has been hindered by a lack of techniques that identify metabolites and their biochemical relationship in complex mixtures. The recent application of Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) to metabolomic analysis suggests a way to tackle the problem. A lower-cost alternative to high-field FTICR-MS, the Orbitrap mass analyzer, promises accelerated activity in this area. Here, we show how the ultra-high mass accuracy and resolution provided by this new generation of mass spectrometers can help to identify metabolites and connect them into metabolic networks. Data from perturbation studies and isotope-tracking experiments can complement this information to create metabolic maps de novo and chart unexplored areas of metabolism.     

Systems analysis of metabolism in the pathogenic trypanosomatid Leishmania major

Systems analyses have facilitated the characterization of metabolic networks of several organisms. We have reconstructed the metabolic network of Leishmania major, a poorly characterized organism that causes cutaneous leishmaniasis in mammalian hosts. This network reconstruction accounts for 560 genes, 1112 reactions, 1101 metabolites and 8 unique subcellular localizations. Using a systems‐based approach, we hypothesized a comprehensive set of lethal single and double gene deletions, some of which were validated using published data with approximately 70% accuracy. Additionally, we generated hypothetical annotations to dozens of previously uncharacterized genes in the L. major genome and proposed a minimal medium for growth. We further demonstrated the utility of a network reconstruction with two proof‐of‐concept examples that yielded insight into robustness of the network in the presence of enzymatic inhibitors and delineation of promastigote/amastigote stage‐specific metabolism. This reconstruction and the associated network analyses of L. major is the first of its kind for a protozoan. It can serve as a tool for clarifying discrepancies between data sources, generating hypotheses that can be experimentally validated and identifying ideal therapeutic targets.     

Prospects for Biotechnology in Oil Palm (Elaeis guineensis) and Coconut (Cocos nucifera) Improvement

Abstract

Research into human metabolism is expanding rapidly due to the emergence of metabolism as a key factor in common diseases. Mathematical modeling of human cellular metabolism has traditionally been performed via kinetic approaches whose applicability for large-scale systems is limited by lack of kinetic constants data. An alternative computational approach bypassing this hurdle called constraint-based modeling (CBM) serves to analyze the function of large-scale metabolic networks by solely relying on simple physical-chemical constraints. However, while extensive research has been performed on constraint-based modeling of microbial metabolism, large-scale modeling of human metabolism is still in its infancy. Utilizing constraint-based modeling to model human cellular metabolism is significantly more complicated than modeling microbial metabolism as in multi-cellular organisms the metabolic behavior varies across cell-types and tissues. It is further complicated due to lack of data on cell type- and tissue-specific metabolite uptake from the surrounding microenvironments and tissue-specific metabolic objective functions. To overcome these problems, several studies suggested CBM methods that integrate metabolic networks with gene expression data that is easily measurable under various conditions. This specific objective functions are expected to improve the prediction accuracy of the presented methods. Such objective functions may be derived based on computational learning that would give optimal correspondence between predicted and measured metabolic phenotypes (Burgard, 2003).

The CBM methods presented here open the way for future computational investigations of metabolic disorders given the relevant expression data. A first attempt to visualize and interpret changes in gene expression data measured following gastric bypass surgery via a genome-scale metabolic network was done by Duarte et al (Duarte, 2007). Another potential application would be the prediction of diagnostic biomarkers for metabolic diseases that could be identified via biofluid metabolomics (Kell, 2007). Towards this goal, we have recently developed a CBM method for predicting metabolic biomarkers for in-born errors of metabolism by searching for changes in metabolite uptake and secretion rate due to genetic alterations (Shlomi, 2009). Incorporating cell type- and tissue-specific gene expression data within this framework can potentially improve the identification of diagnostic biomarkers. Overall, the methods presented here lay the foundation for studying normal and abnormal human cellular metabolism in tissue-specific manner based on commonly measured gene expression data.     

Determination of metabolic fluxes in a non-steady-state system

Estimation of fluxes through metabolic networks from redistribution patterns of (13)C has become a well developed technique in recent years. However, the approach is currently limited to systems at metabolic steady-state; dynamic changes in metabolic fluxes cannot be assessed. This is a major impediment to understanding the behaviour of metabolic networks, because steady-state is not always experimentally achievable and a great deal of information about the control hierarchy of the network can be derived from the analysis of flux dynamics. To address this issue, we have developed a method for estimating non-steady-state fluxes based on the mass-balance of mass isotopomers. This approach allows multiple mass-balance equations to be written for the change in labelling of a given metabolite pool and thereby permits over-determination of fluxes. We demonstrate how linear regression methods can be used to estimate non-steady-state fluxes from these mass balance equations. The approach can be used to calculate fluxes from both mass isotopomer and positional isotopomer labelling information and thus has general applicability to data generated from common spectrometry- or NMR-based analytical platforms. The approach is applied to a GC-MS time-series dataset of (13)C-labelling of metabolites in a heterotrophic Arabidopsis cell suspension culture. Threonine biosynthesis is used to demonstrate that non-steady-state fluxes can be successfully estimated from such data while organic acid metabolism is used to highlight some common issues that can complicate flux estimation. These include multiple pools of the same metabolite that label at different rates and carbon skeleton rearrangements.     

Structural analysis of metabolic networks based on flux centrality

Metabolic reactions are fundamental to living organisms, and a large number of reactions simultaneously occur at a given time in living cells transforming diverse metabolites into each other. There has been an ongoing debate on how to classify metabolites with respect to their importance for metabolic performance, usually based on the analysis of topological properties of genome scale metabolic networks. However, none of these studies have accounted quantitatively for flux in metabolic networks, thus lacking an important component of a cell’s biochemistry.We therefore analyzed a genome scale metabolic network of Escherichia coli by comparing growth under 19 different growth conditions, using flux balance analysis and weighted network centrality investigation. With this novel concept of flux centrality we generated metabolite rankings for each particular growth condition. In contrast to the results of conventional analysis of genome scale metabolic networks, different metabolites were top-ranking dependent on the growth condition. At the same time, several metabolites were consistently among the high ranking ones. Those are associated with pathways that have been described by biochemists as the most central part of metabolism, such as glycolysis, tricarboxylic acid cycle and pentose phosphate pathway. The values for the average path length of the analyzed metabolite networks were between 10.5 and 12.6, supporting recent findings that the metabolic network of E. coli is not a small-world network.     

Robust detection and verification of linear relationships to generate metabolic networks using estimates of technical errors

Background  

The size and magnitude of the metabolome, the ratio between individual metabolites and the response of metabolic networks is controlled by multiple cellular factors. A tight control over metabolite ratios will be reflected by a linear relationship of pairs of metabolite due to the flexibility of metabolic pathways. Hence, unbiased detection and validation of linear metabolic variance can be interpreted in terms of biological control. For robust analyses, criteria for rejecting or accepting linearities need to be developed despite technical measurement errors. The entirety of all pair wise linear metabolic relationships then yields insights into the network of cellular regulation.     

Predicting functional associations from metabolism using bi-partite network algorithms

Background  

Metabolic reconstructions contain detailed information about metabolic enzymes and their reactants and products. These networks can be used to infer functional associations between metabolic enzymes. Many methods are based on the number of metabolites shared by two enzymes, or the shortest path between two enzymes. Metabolite sharing can miss associations between non-consecutive enzymes in a serial pathway, and shortest-path algorithms are sensitive to high-degree metabolites such as water and ATP that create connections between enzymes with little functional similarity.     

Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric,thermodynamic and physiological constraints

Mathematical modeling is an essential tool for the comprehensive understanding of cell metabolism and its interactions with the environmental and process conditions. Recent developments in the construction and analysis of stoichiometric models made it possible to define limits on steady-state metabolic behavior using flux balance analysis. However, detailed information on enzyme kinetics and enzyme regulation is needed to formulate kinetic models that can accurately capture the dynamic metabolic responses. The use of mechanistic enzyme kinetics is a difficult task due to uncertainty in the kinetic properties of enzymes. Therefore, the majority of recent works considered only mass action kinetics for reactions in metabolic networks. Herein, we applied the optimization and risk analysis of complex living entities (ORACLE) framework and constructed a large-scale mechanistic kinetic model of optimally grown Escherichia coli. We investigated the complex interplay between stoichiometry, thermodynamics, and kinetics in determining the flexibility and capabilities of metabolism. Our results indicate that enzyme saturation is a necessary consideration in modeling metabolic networks and it extends the feasible ranges of metabolic fluxes and metabolite concentrations. Our results further suggest that enzymes in metabolic networks have evolved to function at different saturation states to ensure greater flexibility and robustness of cellular metabolism.     

From annotated genomes to metabolic flux models and kinetic parameter fitting

Significant advances in system-level modeling of cellular behavior can be achieved based on constraints derived from genomic information and on optimality hypotheses. For steady-state models of metabolic networks, mass conservation and reaction stoichiometry impose linear constraints on metabolic fluxes. Different objectives, such as maximization of growth rate or minimization of flux distance from a reference state, can be tested in different organisms and conditions. In particular, we have suggested that the metabolic properties of mutant bacterial strains are best described by an algorithm that performs a minimization of metabolic adjustment (MOMA) upon gene deletion. The increasing availability of many annotated genomes paves the way for a systematic application of these flux balance methods to a large variety of organisms. However, such a high throughput goal crucially depends on our capacity to build metabolic flux models in a fully automated fashion. Here we describe a pipeline for generating models from annotated genomes and discuss the current obstacles to full automation. In addition, we propose a framework for the integration of flux modeling results and high throughput proteomic data, which can potentially help in the inference of whole-cell kinetic parameters.