基元模式分析在生物网络和途径分析中的应用

基元模式分析是应用最广泛的代谢途径分析方法。基元模式分析的研究对象从代谢网络发展到信号传导网络;研究尺度从细胞到生物反应器,甚至生态系统;数学描述从稳态分解到动态解析;研究领域从微生物代谢到人类疾病。以下综述了基元模式分析的算法和软件开发现状,以及其在代谢途径与鲁棒性、代谢通量分解、稳态代谢通量分析、动态模型与生物过程模拟、网络结构与调控、菌株设计和信号传导网络等方面的应用。开发新的算法解决组合爆炸问题,探索基元模式与代谢调控的关系以及提高菌株设计算法效率是今后基元模式的重要发展方向。     

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems-oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta-ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations.     

代谢网络定量分析研究进展

综述了代谢工程中代谢控制分析、代谢通量分析、生化系统理论、途径分析、控制论模型等定量分析方法的基本理论,以实例说明了这些方法的应用,并对代谢分析方法的发展进行了展望。     

Application of dynamic metabolic flux analysis for process modeling: Robust flux estimation with regularization,confidence bounds,and selection of elementary modes

In macroscopic dynamic models of fermentation processes, elementary modes (EM) derived from metabolic networks are often used to describe the reaction stoichiometry in a simplified manner and to build predictive models by parameterizing kinetic rate equations for the EM. In this procedure, the selection of a set of EM is a key step which is followed by an estimation of their reaction rates and of the associated confidence bounds. In this paper, we present a method for the computation of reaction rates of cellular reactions and EM as well as an algorithm for the selection of EM for process modeling. The method is based on the dynamic metabolic flux analysis (DMFA) proposed by Leighty and Antoniewicz (2011, Metab Eng, 13(6), 745–755) with additional constraints, regularization and analysis of uncertainty. Instead of using estimated uptake or secretion rates, concentration measurements are used directly to avoid an amplification of measurement errors by numerical differentiation. It is shown that the regularized DMFA for EM method is significantly more robust against measurement noise than methods using estimated rates. The confidence intervals for the estimated reaction rates are obtained by bootstrapping. For the selection of a set of EM for a given st oichiometric model, the DMFA for EM method is combined with a multiobjective genetic algorithm. The method is applied to real data from a CHO fed-batch process. From measurements of six fed-batch experiments, 10 EM were identified as the smallest subset of EM based upon which the data can be described sufficiently accurately by a dynamic model. The estimated EM reaction rates and their confidence intervals at different process conditions provide useful information for the kinetic modeling and subsequent process optimization.     

Estimation of metabolic fluxes, expression levels and metabolite dynamics of a secondary metabolic pathway in potato using label pulse-feeding experiments combined with kinetic network modelling and simulation

In this paper we present a method that allows dynamic flux analysis without a priori kinetic knowledge. This method was developed and validated using the pulse-feeding experimental data obtained in our previous study (Matsuda et al., 2005), in which incorporation of exogenously applied l-phenylalanine-d(5) into seven phenylpropanoid metabolites in potato tubers was determined. After identification of the topology of the metabolic network of these biosynthetic pathways, the system was described by dynamic mass balances in combination with power-law kinetics. After the first simulations, some reactions were removed from the network because they were not contributing significantly to network behaviour. As a next step, the exponents of the power-law kinetics were identified and then kept at fixed values during further analysis. The model was tested for statistical reliability using Monte Carlo simulations. Most fluxes could be identified with high accuracy. The two test cases, control and after elicitation, were clearly distinguished, and with elicitation fluxes to N-p-coumaroyloctopamine (pCO) and N-p-coumaroyltyramine (pCT) increased significantly, whereas those for chlorogenic acid (CGA) and p-coumaroylshikimate decreased significantly. According to the model, increases in the first two fluxes were caused by induction/derepression mechanisms. The decreases in the latter two fluxes were caused by decreased concentrations of their substrates, which in turn were caused by increased activity of the pCO- and pCT-producing enzymes. Flux-control analysis showed that, in most cases, flux control was changed after application of elicitor. Thus the results revealed potential targets for improving actions against tissue wounding and pathogen attack.     

A functional genomics approach using metabolomics and in silico pathway analysis

In the field of functional genomics increasing effort is being undertaken to analyze the function of orphan genes using metabolome data. Improved analytical equipment allows screening simultaneously for a high number of metabolites. Such metabolite profiles are analyzed using multivariate data analysis techniques and changes in the genotype will in many cases lead to different metabolite profiles. Here, a theoretical framework that may be applied to identify the function of orphan genes is presented. The approach is based on a combination of metabolome analysis combined with in silico pathway analysis. Pathway analysis may be carried out using convex analysis and a change in the active pathway structure of deletion mutants expressed in a different metabolite profile may disclose the function or the functional class of an orphan gene. The concept is illustrated using a simplified model for growth of Saccharomyces cerevisiae.     

代谢流量比率分析及其在代谢工程中的应用

细胞内代谢反应流量在系统理解细胞代谢特性和指导代谢工程改造等方面都起着重要的作用。由于代谢流量难以直接测量得到,在很多情况下通过跟踪稳定同位素在代谢网络中的转移并进行相应的模型计算能有效地定量代谢流量。代谢流量比率分析法能够高度体现系统的生物化学真实性、辨别细胞代谢网络的拓扑结构,并且能够相对简单快速地定量反应速率等,因此受到代谢工程研究者越来越多的重视。以下着重介绍并讨论了利用代谢物同位体分布信息分析关键代谢节点合成途径的流量比率、基于流量比率的代谢流量解析、以及应用于代谢工程等的相关原理、实验测量、数据分析、使用条件等,以期充分发挥代谢流量比率分析法的优势,并将其拓展推广至更多细胞体系的代谢特性阐明和代谢工程改造中去。     

Metabolic pathway analysis of a recombinant yeast for rational strain development

Elementary mode analysis has been used to study a metabolic pathway model of a recombinant Saccharomyces cerevisiae system that was genetically engineered to produce the bacterial storage compound poly-beta-hydroxybutyrate (PHB). The model includes biochemical reactions from the intermediary metabolism and takes into account cellular compartmentalization as well as the reversibility/irreversibility of the reactions. The reaction network connects the production and/or consumption of eight external metabolites including glucose, acetate, glycerol, ethanol, PHB, CO(2), succinate, and adenosine triphosphate (ATP). Elementary mode analysis of the wild-type S. cerevisiae system reveals 241 unique reaction combinations that balance the eight external metabolites. When the recombinant PHB pathway is included, and when the reaction model is altered to simulate the experimental conditions when PHB accumulates, the analysis reveals 20 unique elementary modes. Of these 20 modes, 7 produce PHB with the optimal mode having a theoretical PHB carbon yield of 0.67. Elementary mode analysis was also used to analyze the possible effects of biochemical network modifications and altered culturing conditions. When the natively absent ATP citrate-lyase activity is added to the recombinant reaction network, the number of unique modes increases from 20 to 496, with 314 of these modes producing PHB. With this topological modification, the maximum theoretical PHB carbon yield increases from 0.67 to 0.83. Adding a transhydrogenase reaction to the model also improves the theoretical conversion of substrate into PHB. The recombinant system with the transhydrogenase reaction but without the ATP citrate-lyase reaction has an increase in PHB carbon yield from 0.67 to 0.71. When the model includes both the ATP citrate-lyase reaction and the transhydrogenase reaction, the maximum theoretical carbon yield increases to 0.84. The reaction model was also used to explore the possibility of producing PHB under anaerobic conditions. In the absence of oxygen, the recombinant reaction network possesses two elementary modes capable of producing PHB. Interestingly, both modes also produce ethanol. Elementary mode analysis provides a means of deconstructing complex metabolic networks into their basic functional units. This information can be used for analyzing existing pathways and for the rational design of further modifications that could improve the system's conversion of substrate into product.     

基于图论的代谢网络中流通代谢物处理新方法

使用图论方法对基因组尺度代谢网络进行途径搜索是目前途径设计中最常用的方法之一.然而,由于流通代谢物(如H+、H2O、CO2和ATP等)的影响,图论搜索得到的途径在生物学上经常是不可行的.虽然前人提出了一些方法对流通代谢物进行处理,但均存在一些问题,目前还没有标准的处理方法.文中基于流通代谢物在反应中作为转移磷酸等功能基...     

Optimization-based metabolic control analysis

In this work, a novel optimization-based metabolic control analysis (OMCA) method is introduced for reducing data requirement for metabolic control analysis (MCA). It is postulated that using the optimal control approach, the fluxes in a metabolic network are correlated to metabolite concentrations and enzyme activities as a state-feedback control system that is optimal with respect to a homeostasis objective. It is then shown that the optimal feedback gains are directly related to the elasticity coefficients (ECs) of MCA. This approach requires determination of the relative "importance" of metabolites and fluxes for the system, which is possible with significantly reduced experimental data, as compared with typical MCA requirements. The OMCA approach is applied to a top-down control model of glycolysis in hepatocytes. It is statistically demonstrated that the OMCA model is capable of predicting the ECs observed experimentally with few exceptions. Further, an OMCA-based model reconciliation study shows that the modification of four assumed stoichiometric coefficients in the model can explain most of the discrepancies, with the exception of elasticities with respect to the NADH/NAD ratio.     

Inferring meaningful pathways in weighted metabolic networks

An approach is presented for computing meaningful pathways in the network of small molecule metabolism comprising the chemical reactions characterized in all organisms. The metabolic network is described as a weighted graph in which all the compounds are included, but each compound is assigned a weight equal to the number of reactions in which it participates. Path finding is performed in this graph by searching for one or more paths with lowest weight. Performance is evaluated systematically by computing paths between the first and last reactions in annotated metabolic pathways, and comparing the intermediate reactions in the computed pathways to those in the annotated ones. For the sake of comparison, paths are computed also in the un-weighted raw (all compounds and reactions) and filtered (highly connected pool metabolites removed) metabolic graphs, respectively. The correspondence between the computed and annotated pathways is very poor (<30%) in the raw graph; increasing to approximately 65% in the filtered graph; reaching approximately 85% in the weighted graph. Considering the best-matching path among the five lightest paths increases the correspondence to 92%, on average. We then show that the average distance between pairs of metabolites is significantly larger in the weighted graph than in the raw unfiltered graph, suggesting that the small-world properties previously reported for metabolic networks probably result from irrelevant shortcuts through pool metabolites. In addition, we provide evidence that the length of the shortest path in the weighted graph represents a valid measure of the "metabolic distance" between enzymes. We suggest that the success of our simplistic approach is rooted in the high degree of specificity of the reactions in metabolic pathways, presumably reflecting thermodynamic constraints operating in these pathways. We expect our approach to find useful applications in inferring metabolic pathways in newly sequenced genomes.     

Modes and cuts in metabolic networks: complexity and algorithms

Constraint-based approaches recently brought new insight into our understanding of metabolism. By making very simple assumptions such as that the system is at steady-state and some reactions are irreversible, and without requiring kinetic parameters, general properties of the system can be derived. A central concept in this methodology is the notion of an elementary mode (EM for short) which represents a minimal functional subsystem. The computation of EMs still forms a limiting step in metabolic studies and several algorithms have been proposed to address this problem leading to increasingly faster methods. However, although a theoretical upper bound on the number of elementary modes that a network may possess has been established, surprisingly, the complexity of this problem has never been systematically studied. In this paper, we give a systematic overview of the complexity of optimisation problems related to modes. We first establish results regarding network consistency. Most consistency problems are easy, i.e., they can be solved in polynomial time. We then establish the complexity of finding and counting elementary modes. We show in particular that finding one elementary mode is easy but that this task becomes hard when a specific EM (i.e. an EM containing some specified reactions) is sought. We then show that counting the number of elementary modes is musical sharpP-complete. We emphasize that the easy problems can be solved using currently existing software packages. We then analyse the complexity of a closely related task which is the computation of so-called minimum reaction cut sets and we show that this problem is hard. We then present two positive results which both allow to avoid computing EMs as a prior to the computation of reaction cuts. The first one is a polynomial approximation algorithm for finding a minimum reaction cut set. The second one is a test for verifying whether a set of reactions constitutes a reaction cut; this test can be readily included in existing algorithms to improve their performance. Finally, we discuss the complexity of other cut-related problems.     

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome‐wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single‐nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype‐to‐phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme‐coding genes. To this end, we also constructed a genome‐scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth‐determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.     

Metabolomics assisted metabolic network modeling and network wide analysis of metabolites in microbiology

Metabolomics is the science of qualitatively and quantitatively analyzing low molecular weight metabolites occur in a given biological system. It provides valuable information to elucidate the functional roles and relations of different metabolites in a metabolic pathway. In recent years, a large amount of research on microbial metabolomics has been conducted. It has become a useful tool for achieving highly efficient synthesis of target metabolites. At the same time, many studies have been conducted over the years in order to integrate metabolomics data into metabolic network modeling, which has yielded many exciting results. Additionally, metabolomics also shows great advantages in analyzing the relationship of metabolites network wide. Integrating metabolomics data into metabolic network construction and applying it in network wide analysis of cell metabolism would further improve our ability to control cellular metabolism and optimize the design of cell factories for the overproduction of valuable biochemicals. This review will examine recent progress in the application of metabolomics approaches in metabolic network modeling and network wide analysis of microbial cell metabolism.     

A review on metabolic pathway analysis with emphasis on isotope labeling approach

The recent progress on metabolic systems engineering was reviewed based on our recent research results in terms of (1) metabolic signal flow diagram approach, (2) metabolic flux analysis (MFA) in particular with intracellular isotopomer distribution using NMR and/or GC-MS, (3) synthesis and optimization of metabolic flux distribution (MFD), (4) modification of MFD by gene manipulation and by controlling culture environment, (5) metabolic control analysis (MCA), (6) design of metabolic regulation structure, and (7) identification of unknown pathways with isotope tracing by NMR. The main characteristics of metabolic engineering is to treat metabolism as a network or entirety instead of individual reactions. The applications were made for poly-3-hydroxybutyrate (PHB) production usingRalstonia eutropha and recombinantEscherichia coli, lactate production by recombinantSaccharomyces cerevisiae, pyruvate production by vitamin auxotrophic yeastToluropsis glabrata, lysine production usingCorynebacterium glutamicum, and energetic analysis of photosynthesic microorganisms such as Cyanobateria. The characteristics of each approach were reviewed with their applications. The approach based on isotope labeling experiments gives reliable and quantitative results for metabolic flux analysis. It should be recognized that the next stage should be toward the investigation of metabolic flux analysis with gene and protein expressions to uncover the metabolic regulation in relation to genetic modification and/or the change in the culture condition.