Reconstruction of the central carbon metabolism of Aspergillus niger.

The topology of central carbon metabolism of Aspergillus niger was identified and the metabolic network reconstructed, by integrating genomic, biochemical and physiological information available for this microorganism and other related fungi. The reconstructed network may serve as a valuable database for annotation of genes identified in future genome sequencing projects on aspergilli. Based on the metabolic reconstruction, a stoichiometric model was set up that includes 284 metabolites and 335 reactions, of which 268 represent biochemical conversions and 67 represent transport processes between the different intracellular compartments and between the cell and the extracellular medium. The stoichiometry of the metabolic reactions was used in combination with biosynthetic requirements for growth and pseudo-steady state mass balances over intracellular metabolites for the quantification of metabolic fluxes using metabolite balancing. This framework was employed to perform an in silico characterisation of the phenotypic behaviour of A. niger grown on different carbon sources. The effects on growth of single reaction deletions were assessed and essential biochemical reactions were identified for different carbon sources. Furthermore, application of the stoichiometric model for assessing the metabolic capabilities of A. niger to produce metabolites was evaluated by using succinate production as a case study.     

Robustness analysis of the Escherichia coli metabolic network

Genomic, biochemical, and strain-specific data can be assembled to define an in silico representation of the metabolic network for a select group of single cellular organisms. Flux-balance analysis and phenotypic phase planes derived therefrom have been developed and applied to analyze the metabolic capabilities and characteristics of Escherichia coli K-12. These analyses have shown the existence of seven essential reactions in the central metabolic pathways (glycolysis, pentose phosphate pathway, tricarboxylic acid cycle) for the growth in glucose minimal media. The corresponding seven gene products can be grouped into three categories: (1) pentose phosphate pathway genes, (2) three-carbon glycolytic genes, and (3) tricarboxylic acid cycle genes. Here we develop a procedure that calculates the sensitivity of optimal cellular growth to altered flux levels of these essential gene products. The results indicate that the E. coli metabolic network is robust with respect to the flux levels of these enzymes. The metabolic flux in the transketolase and the tricarboxylic acid cycle reactions can be reduced to 15% and 19%, respectively, of the optimal value without significantly influencing the optimal growth flux. The metabolic network also exhibited robustness with respect to the ribose-5-phosphate isomerase, and the ribose-5-phosephate isomerase flux was reduced to 28% of the optimal value without significantly effecting the optimal growth flux. The metabolic network exhibited limited robustness to the three-carbon glycolytic fluxes both increased and decreased. The development presented another dimension to the use of FBA to study the capabilities of metabolic networks.     

基于约束建模法的结核菌H37Rv代谢网络分析

基于结核分枝杆菌国际标准强毒株H37Rv菌株的基因组尺度代谢网络模型iNJ661进行分析,以寻找代谢网络中培养基的关键成分和必要基因.该研究在Matlab平台上利用COBRA工具箱,采用基于约束的建模方法进行动态生长模拟、解空间抽样在酶活性水平上的具体化和基因删除模拟实验.结果发现培养基成分中铵盐、三价铁盐、磷酸盐、硫酸盐、甘油等可影响H37Rv的生长;培养基中去除磷酸盐后十种酶均在不同程度上受到抑制,其中丙糖磷酸异构酶、3-磷酸甘油醛脱氢酶、磷酸甘油酸变位酶、烯醇酶受限明显.通过基因删除得出188个必要基因以及非必要基因中的16个致死基因对.基于约束建模分析可初步了解结核杆菌H37Rv菌株代谢网络的性质,可为后续相关研究提供参考和借鉴.     

Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis

A genome-scale metabolic network reconstruction for Clostridium acetobutylicum (ATCC 824) was carried out using a new semi-automated reverse engineering algorithm. The network consists of 422 intracellular metabolites involved in 552 reactions and includes 80 membrane transport reactions. The metabolic network illustrates the reliance of clostridia on the urea cycle, intracellular L-glutamate solute pools, and the acetylornithine transaminase for amino acid biosynthesis from the 2-oxoglutarate precursor. The semi-automated reverse engineering algorithm identified discrepancies in reaction network databases that are major obstacles for fully automated network-building algorithms. The proposed semi-automated approach allowed for the conservation of unique clostridial metabolic pathways, such as an incomplete TCA cycle. A thermodynamic analysis was used to determine the physiological conditions under which proposed pathways (e.g., reverse partial TCA cycle and reverse arginine biosynthesis pathway) are feasible. The reconstructed metabolic network was used to create a genome-scale model that correctly characterized the butyrate kinase knock-out and the asolventogenic M5 pSOL1 megaplasmid degenerate strains. Systematic gene knock-out simulations were performed to identify a set of genes encoding clostridial enzymes essential for growth in silico.     

Characterization of metabolism in the Fe(III)-reducing organism Geobacter sulfurreducens by constraint-based modeling

Geobacter sulfurreducens is a well-studied representative of the Geobacteraceae, which play a critical role in organic matter oxidation coupled to Fe(III) reduction, bioremediation of groundwater contaminated with organics or metals, and electricity production from waste organic matter. In order to investigate G. sulfurreducens central metabolism and electron transport, a metabolic model which integrated genome-based predictions with available genetic and physiological data was developed via the constraint-based modeling approach. Evaluation of the rates of proton production and consumption in the extracellular and cytoplasmic compartments revealed that energy conservation with extracellular electron acceptors, such as Fe(III), was limited relative to that associated with intracellular acceptors. This limitation was attributed to lack of cytoplasmic proton consumption during reduction of extracellular electron acceptors. Model-based analysis of the metabolic cost of producing an extracellular electron shuttle to promote electron transfer to insoluble Fe(III) oxides demonstrated why Geobacter species, which do not produce shuttles, have an energetic advantage over shuttle-producing Fe(III) reducers in subsurface environments. In silico analysis also revealed that the metabolic network of G. sulfurreducens could synthesize amino acids more efficiently than that of Escherichia coli due to the presence of a pyruvate-ferredoxin oxidoreductase, which catalyzes synthesis of pyruvate from acetate and carbon dioxide in a single step. In silico phenotypic analysis of deletion mutants demonstrated the capability of the model to explore the flexibility of G. sulfurreducens central metabolism and correctly predict mutant phenotypes. These results demonstrate that iterative modeling coupled with experimentation can accelerate the understanding of the physiology of poorly studied but environmentally relevant organisms and may help optimize their practical applications.     

Development of chemically defined medium for Mannheimia succiniciproducens based on its genome sequence

This study presents a novel methodology for the development of a chemically defined medium (CDM) using genome-scale metabolic network and flux balance analysis. The genome-based in silico analysis identified two amino acids and four vitamins as non-substitutable essential compounds to be supplemented to a minimal medium for the sustainable growth of Mannheimia succiniciproducens, while no substitutable essential compounds were identified. The in silico predictions were verified by cultivating the cells on a CDM containing the six non-substitutable essential compounds, and it was further demonstrated by observing no cell growth on the CDM lacking any one of the non-substitutable essentials. An optimal CDM for the enhancement of cell growth and succinic acid production, as a target product, was formulated with a single-addition technique. The fermentation on the optimal CDM increased the succinic acid productivity by 36%, the final succinic acid concentration by 17%, and the succinic acid yield on glucose by 15% compared to the cultivation using a complex medium. The optimal CDM also lowered the sum of the amounts of by-products (acetic, formic, and lactic acids) by 30%. The strategy reported in this paper should be generally applicable to the development of CDMs for other organisms, whose genome sequences are available.     

Genome-scale reconstruction of metabolic network in Bacillus subtilis based on high-throughput phenotyping and gene essentiality data

In this report, a genome-scale reconstruction of Bacillus subtilis metabolism and its iterative development based on the combination of genomic, biochemical, and physiological information and high-throughput phenotyping experiments is presented. The initial reconstruction was converted into an in silico model and expanded in a four-step iterative fashion. First, network gap analysis was used to identify 48 missing reactions that are needed for growth but were not found in the genome annotation. Second, the computed growth rates under aerobic conditions were compared with high-throughput phenotypic screen data, and the initial in silico model could predict the outcomes qualitatively in 140 of 271 cases considered. Detailed analysis of the incorrect predictions resulted in the addition of 75 reactions to the initial reconstruction, and 200 of 271 cases were correctly computed. Third, in silico computations of the growth phenotypes of knock-out strains were found to be consistent with experimental observations in 720 of 766 cases evaluated. Fourth, the integrated analysis of the large-scale substrate utilization and gene essentiality data with the genome-scale metabolic model revealed the requirement of 80 specific enzymes (transport, 53; intracellular reactions, 27) that were not in the genome annotation. Subsequent sequence analysis resulted in the identification of genes that could be putatively assigned to 13 intracellular enzymes. The final reconstruction accounted for 844 open reading frames and consisted of 1020 metabolic reactions and 988 metabolites. Hence, the in silico model can be used to obtain experimentally verifiable hypothesis on the metabolic functions of various genes.     

Probing the performance limits of the Escherichia coli metabolic network subject to gene additions or deletions.

An optimization-based procedure for studying the response of metabolic networks after gene knockouts or additions is introduced and applied to a linear flux balance analysis (FBA) Escherichia coli model. Both the gene addition problem of optimally selecting which foreign genes to recombine into E. coli, as well as the gene deletion problem of removing a given number of existing ones, are formulated as mixed-integer optimization problems using binary 0-1 variables. The developed modeling and optimization framework is tested by investigating the effect of gene deletions on biomass production and addressing the maximum theoretical production of the 20 amino acids for aerobic growth on glucose and acetate substrates. In the gene deletion study, the smallest gene set necessary to achieve maximum biomass production in E. coli is determined for aerobic growth on glucose. The subsequent gene knockout analysis indicates that biomass production decreases monotonically, rendering the metabolic network incapable of growth after only 18 gene deletions. In the gene addition study, the E. coli flux balance model is augmented with 3,400 non-E. coli reactions from the KEGG database to form a multispecies model. This model is referred to as the Universal model. This study reveals that the maximum theoretical production of six amino acids could be improved by the addition of only one or two genes to the native amino acid production pathway of E. coli, even though the model could choose from 3,400 foreign reaction candidates. Specifically, manipulation of the arginine production pathway showed the most promise with 8.75% and 9.05% predicted increases with the addition of genes for growth on glucose and acetate, respectively. The mechanism of all suggested enhancements is either by: 1) improving the energy efficiency and/or 2) increasing the carbon conversion efficiency of the production route.     

Systemic metabolic reactions are obtained by singular value decomposition of genome-scale stoichiometric matrices

Genome-scale metabolic networks can be reconstructed. The systemic biochemical properties of these networks can now be studied. Here, genome-scale reconstructed metabolic networks were analysed using singular value decomposition (SVD). All the individual biochemical conversions contained in a reconstructed metabolic network are described by a stoichiometric matrix (S). SVD of S led to the definition of the underlying modes that characterize the overall biochemical conversions that take place in a network and rank-ordered their importance. The modes were shown to correspond to systemic biochemical reactions and they could be used to identify the groups and clusters of individual biochemical reactions that drive them. Comparative analysis of the Escherichia coli, Haemophilus influenzae, and Helicobacter pylori genome-scale metabolic networks showed that the four dominant modes in all three networks correspond to: (1) the conversion of ATP to ADP, (2) redox metabolism of NADP, (3) proton-motive force, and (4) inorganic phosphate metabolism. The sets of individual metabolic reactions deriving these systemic conversions, however, differed among the three organisms. Thus, we can now define systemic metabolic reactions, or eigen-reactions, for the study of systems biology of metabolism and have a basis for comparing the overall properties of genome-specific metabolic networks.     

Reconstruction of the Saccharopolyspora erythraea genome-scale model and its use for enhancing erythromycin production

Genome-scale metabolic reconstructions are routinely used for the analysis and design of metabolic engineering strategies for production of primary metabolites. The use of such reconstructions for metabolic engineering of antibiotic production is not common due to the lack of simple design algorithms in the absence of a cellular growth objective function. Here, we present the metabolic network reconstruction for the erythromycin producer Saccharopolyspora erythraea NRRL23338. The model was manually curated for primary and secondary metabolism pathways and consists of 1,482 reactions (2,075 genes) and 1,646 metabolites. As part of the model validation, we explored the potential benefits of supplying amino acids and identified five amino acids “compatible” with erythromycin production, whereby if glucose is supplemented with this amino acid on a carbon mole basis, the in silico model predicts that high erythromycin yield is possible without lowering biomass yield. Increased erythromycin titre was confirmed for four of the five amino acids, namely valine, isoleucine, threonine and proline. In bioreactor experiments, supplementation with 2.5?% carbon mole of valine increased the growth rate by 20?% and simultaneously the erythromycin yield on biomass by 50?%. The model presented here can be used as a framework for the future integration of high-throughput biological data sets in S. erythraea and ultimately to realise strain designs capable of increasing erythromycin production closer to the theoretical yield.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.     

Characterization of Metabolism in the Fe(III)-Reducing Organism Geobacter sulfurreducens by Constraint-Based Modeling

Geobacter sulfurreducens is a well-studied representative of the Geobacteraceae, which play a critical role in organic matter oxidation coupled to Fe(III) reduction, bioremediation of groundwater contaminated with organics or metals, and electricity production from waste organic matter. In order to investigate G. sulfurreducens central metabolism and electron transport, a metabolic model which integrated genome-based predictions with available genetic and physiological data was developed via the constraint-based modeling approach. Evaluation of the rates of proton production and consumption in the extracellular and cytoplasmic compartments revealed that energy conservation with extracellular electron acceptors, such as Fe(III), was limited relative to that associated with intracellular acceptors. This limitation was attributed to lack of cytoplasmic proton consumption during reduction of extracellular electron acceptors. Model-based analysis of the metabolic cost of producing an extracellular electron shuttle to promote electron transfer to insoluble Fe(III) oxides demonstrated why Geobacter species, which do not produce shuttles, have an energetic advantage over shuttle-producing Fe(III) reducers in subsurface environments. In silico analysis also revealed that the metabolic network of G. sulfurreducens could synthesize amino acids more efficiently than that of Escherichia coli due to the presence of a pyruvate-ferredoxin oxidoreductase, which catalyzes synthesis of pyruvate from acetate and carbon dioxide in a single step. In silico phenotypic analysis of deletion mutants demonstrated the capability of the model to explore the flexibility of G. sulfurreducens central metabolism and correctly predict mutant phenotypes. These results demonstrate that iterative modeling coupled with experimentation can accelerate the understanding of the physiology of poorly studied but environmentally relevant organisms and may help optimize their practical applications.     

Comparative metabolic capabilities for Micrococcus luteus NCTC 2665, the "Fleming" strain, and actinobacteria

Putative gene predictions of the Gram positive actinobacteria Micrococcus luteus (NCTC 2665, "Fleming strain") was used to construct a genome scale reconstruction of the metabolic network for this organism. The metabolic network comprises 586 reactions and 551 metabolites, and accounts for 21% of the genes in the genome. The reconstruction was based on the annotated genome and available biochemical information. M. luteus has one of the smallest genomes of actinobacteria with a circular chromosome of 2,501,097 base pairs and a GC content of 73%. The metabolic pathways required for biomass production in silico were determined based on earlier models of actinobacteria. The in silico network is used for metabolic comparison of M. luteus with other actinomycetes, and hence provides useful information for possible future biotechnological exploitation of this organism, e.g., for production of biofuels.     

An in silico compartmentalized metabolic model of Brassica napus enables the systemic study of regulatory aspects of plant central metabolism

Biochemical network reconstructions represent valuable tools for the computational metabolic modeling of organisms that present a great biotechnological interest. An in silico multi-compartmental model of the central metabolism of the plant Brassica napus (Rapeseed) was constructed, aiming to investigate the metabolic properties of the Brassicaceae family. This family comprises many plants with major importance for the energy and nutrition sector, including the model plant Arabidopsis thaliana. The model utilized as objective function to be subsequently optimized, the biomass production of rapeseed developing embryos, which are characterized by a very high, oil content, up to 60% of biomass weight. In order to study global network properties of seed metabolism, various methods were employed, like Flux Balance Analysis, Principal Component Analysis of the flux space and reaction deletion studies, which simulate the effect of gene knock-out experiments. The model successfully simulated seed growth during the stage of oil accumulation and provided insight, regarding certain aspects of network plasticity, with the emphasis given in lipid biosynthesis regulation.     

基于文献挖掘的巨大芽胞杆菌代谢网络模型的构建与分析

【目的】通过挖掘实验性文献,建立巨大芽胞杆菌事实型代谢网络模型,以详尽解析生理特性,优化其生理功能。【方法】从PubMed、Derwent Innovations Index、中国知网等公共文献(专利)数据库中获取与巨大芽胞杆菌(Bacillus megaterium)相关的实验性文献建立本地文献数据库。采用文献挖掘工具获取功能基因、酶、代谢物和生化反应等信息,以其为基础构建代谢网络粗模型,进一步借助KEGG等数据库修正以及Matlab程序的模拟得到精细模型(系统生物学标记语言的形式)。【结果】最终的精细模型共有292个生化反应、378个代谢物、220个酶和217个基因。以1.62 mmol/g cell/h的葡萄糖底物吸收速率为限制性条件,模拟的菌体比生长速率为0.089 h-1,略低于实验值0.11 h-1。此外,嘧啶代谢途径的单基因敲除模拟结果表明,准确率为90%。【结论】该代谢网络模型涵盖了中心代谢途径、维生素B12合成途径和氨基酸代谢途径,并在一定程度上反映了营养底物与基因对巨大芽胞杆菌生长性能的影响。     

Reconstruction and In Silico Analysis of Metabolic Network for an Oleaginous Yeast,Yarrowia lipolytica

With the emergence of energy scarcity, the use of renewable energy sources such as biodiesel is becoming increasingly necessary. Recently, many researchers have focused their minds on Yarrowia lipolytica, a model oleaginous yeast, which can be employed to accumulate large amounts of lipids that could be further converted to biodiesel. In order to understand the metabolic characteristics of Y. lipolytica at a systems level and to examine the potential for enhanced lipid production, a genome-scale compartmentalized metabolic network was reconstructed based on a combination of genome annotation and the detailed biochemical knowledge from multiple databases such as KEGG, ENZYME and BIGG. The information about protein and reaction associations of all the organisms in KEGG and Expasy-ENZYME database was arranged into an EXCEL file that can then be regarded as a new useful database to generate other reconstructions. The generated model iYL619_PCP accounts for 619 genes, 843 metabolites and 1,142 reactions including 236 transport reactions, 125 exchange reactions and 13 spontaneous reactions. The in silico model successfully predicted the minimal media and the growing abilities on different substrates. With flux balance analysis, single gene knockouts were also simulated to predict the essential genes and partially essential genes. In addition, flux variability analysis was applied to design new mutant strains that will redirect fluxes through the network and may enhance the production of lipid. This genome-scale metabolic model of Y. lipolytica can facilitate system-level metabolic analysis as well as strain development for improving the production of biodiesels and other valuable products by Y. lipolytica and other closely related oleaginous yeasts.     

Determining important regulatory relations of amino acids from dynamic network analysis of plasma amino acids

The changes in the concentrations of plasma amino acids do not always follow the flow-based metabolic pathway network. We have previously shown that there is a control-based network structure among plasma amino acids besides the metabolic pathway map. Based on this network structure, in this study, we performed dynamic analysis using time-course data of the plasma samples of rats fed single essential amino acid deficient diet. Using S-system model (conceptual mathematical model represented by power-law formalism), we inferred the dynamic network structure which reproduces the actual time-courses within the error allowance of 13.17%. By performing sensitivity analysis, three of the most dominant relations in this network were selected; the control paths from leucine to valine, from methionine to threonine, and from leucine to isoleucine. This result is in good agreement with the biological knowledge regarding branched-chain amino acids, and suggests the biological importance of the effect from methionine to threonine.     

基于约束的建模方法在代谢网络中的应用研究

代谢网络在各种细胞功能和生命过程中发挥着至关重要的作用。随着细胞网络重建工程的迅速发展,可用的基因组水平代谢网络越来越多,因而计算方法在这些网络的结构功能分析中越来越重要。基于约束的建模方法不像图论方法那样仅考虑代谢模型的纯拓扑结构,也不像各种动力学建模方法那样需求详尽的热力学参数,因而极具优势。采用基于约束的建模方法对一个含619个基因,655个代谢物和743个代谢反应的金黄色葡萄球菌(Staphylococcusaureus)代谢网络进行了分析,主要研究了该模型的网络结构特征,以及其最优生长率、动态生长情况和基因删除学习等。本研究提供了一个对金黄色葡萄球菌代谢网络进行约束建模分析的初步框架。     

Functional characterization of UvrD helicases from Haemophilus influenzae and Helicobacter pylori

Haemophilus influenzae and Helicobacter pylori are major bacterial pathogens that face high levels of genotoxic stress within their host. UvrD, a ubiquitous bacterial helicase that plays important roles in multiple DNA metabolic pathways, is essential for genome stability and might, therefore, be crucial in bacterial physiology and pathogenesis. In this study, the functional characterization of UvrD helicase from Haemophilus influenzae and Helicobacter pylori is reported. UvrD from Haemophilus influenzae (HiUvrD) and Helicobacter pylori (HpUvrD) exhibit strong single-stranded DNA-specific ATPase and 3'-5' helicase activities. Mutation of highly conserved arginine (R288) in HiUvrD and glutamate (E206) in HpUvrD abrogated their activities. Both the proteins were able to bind and unwind a variety of DNA structures including duplexes with strand discontinuities and branches, three- and four-way junctions that underpin their role in DNA replication, repair and recombination. HiUvrD required a minimum of 12 nucleotides, whereas HpUvrD preferred 20 or more nucleotides of 3'-single-stranded DNA tail for efficient unwinding of duplex DNA. Interestingly, HpUvrD was able to hydrolyze and utilize GTP for its helicase activity although not as effectively as ATP, which has not been reported to date for UvrD characterized from other organisms. HiUvrD and HpUvrD were found to exist predominantly as monomers in solution together with multimeric forms. Noticeably, deletion of distal C-terminal 48 amino acid residues disrupted the oligomerization of HiUvrD, whereas deletion of 63 amino acids from C-terminus of HpUvrD had no effect on its oligomerization. This study presents the characteristic features and comparative analysis of Haemophilus influenzae and Helicobacter pylori UvrD, and constitutes the basis for understanding the role of UvrD in the biology and virulence of these pathogens.