Modeling hybridoma cell metabolism using a generic genome-scale metabolic model of Mus musculus

The reconstructed cellular metabolic network of Mus musculus, based on annotated genomic data, pathway databases, and currently available biochemical and physiological information, is presented. Although incomplete, it represents the first attempt to collect and characterize the metabolic network of a mammalian cell on the basis of genomic data. The reaction network is generic in nature and attempts to capture the carbon, energy, and nitrogen metabolism of the cell. The metabolic reactions were compartmentalized between the cytosol and the mitochondria, including transport reactions between the compartments and the extracellular medium. The reaction list consists of 872 internal metabolites involved in a total of 1220 reactions, whereof 473 relate to known open reading frames. Initial in silico analysis of the reconstructed model is presented.     

Flux balance analysis of photoautotrophic metabolism

Photosynthesis is the principal process responsible for fixation of inorganic carbon dioxide into organic molecules with sunlight as the energy source. Potentially, many chemicals could be inexpensively produced by photosynthetic organisms. Mathematical modeling of photoautotrophic metabolism is therefore important to evaluate maximum theoretical product yields and to deeply understand the interactions between biochemical energy, carbon fixation, and assimilation pathways. Flux balance analysis based on linear programming is applied to photoautotrophic metabolism. The stoichiometric network of a model photosynthetic prokaryote, Synechocystis sp. PCC 6803, has been reconstructed from genomic data and biochemical literature and coupled with a model of the photophosphorylation processes. Flux map topologies for the hetero-, auto-, and mixotrophic modes of metabolism under conditions of optimal growth were determined and compared. The roles of important metabolic reactions such as the glyoxylate shunt and the transhydrogenase reaction were analyzed. We also theoretically evaluated the effect of gene deletions or additions on biomass yield and metabolic flux distributions.     

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.     

13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches

Conventional metabolic flux analysis uses the information gained from determination of measurable fluxes and a steady-state assumption for intracellular metabolites to calculate the metabolic fluxes in a given metabolic network. The determination of intracellular fluxes depends heavily on the correctness of the assumed stoichiometry including the presence of all reactions with a noticeable impact on the model metabolite balances. Determination of fluxes in complex metabolic networks often requires the inclusion of NADH and NADPH balances, which are subject to controversial debate. Transhydrogenation reactions that transfer reduction equivalents from NADH to NADPH or vice versa can usually not be included in the stoichiometric model, because they result in singularities in the stoichiometric matrix. However, it is the NADPH balance that, to a large extent, determines the calculated flux through the pentose phosphate pathway. Hence, wrong assumptions on the presence or activity of transhydrogenation reactions will result in wrong estimations of the intracellular flux distribution. Using 13C tracer experiments and NMR analysis, flux analysis can be performed on the basis of only well established stoichiometric equations and measurements of the labeling state of intracellular metabolites. Neither NADH/NADPH balancing nor assumptions on energy yields need to be included to determine the intracellular fluxes. Because metabolite balancing methods and the use of 13C labeling measurements are two different approaches to the determination of intracellular fluxes, both methods can be used to verify each other or to discuss the origin and significance of deviations in the results. Flux analysis based entirely on metabolite balancing and flux analysis, including labeling information, have been performed independently for a wild-type strain of Aspergillus oryzae producing alpha-amylase. Two different nitrogen sources, NH4+ and NO3-, have been used to investigate the influence of the NADPH requirements on the intracellular flux distribution. The two different approaches to the calculation of fluxes are compared and deviations in the results are discussed. Copyright 1998 John Wiley & Sons, Inc.     

A Strategy to Calculate the Patterns of Nutrient Consumption by Microorganisms Applying a Two-Level Optimisation Principle to Reconstructed Metabolic Networks

Bacterial responses to environmental changes rely on a complex network of biochemical reactions. The properties of the metabolic network determining these responses can be divided into two groups: the stoichiometric properties, given by the stoichiometry matrix, and the kinetic/thermodynamic properties, given by the rate equations of the reaction steps. The stoichiometry matrix represents the maximal metabolic capabilities of the organism, and the regulatory mechanisms based on the rate laws could be considered as being responsible for the administration of these capabilities. Post-genomic reconstruction of metabolic networks provides us with the stoichiometry matrix of particular strains of microorganisms, but the kinetic aspects of in vivo rate laws are still largely unknown. Therefore, the validity of predictions of cellular responses requiring detailed knowledge of the rate equations is difficult to assert. In this paper, we show that by applying optimisation criteria to the core stoichiometric network of the metabolism of Escherichia coli, and including information about reversibility/irreversibility only of the reaction steps, it is possible to calculate bacterial responses to growth media with different amounts of glucose and galactose. The target was the minimisation of the number of active reactions (subject to attaining a growth rate higher than a lower limit) and subsequent maximisation of the growth rate (subject to the number of active reactions being equal to the minimum previously calculated). Using this two-level target, we were able to obtain by calculation four fundamental behaviours found experimentally: inhibition of respiration at high glucose concentrations in aerobic conditions, turning on of respiration when glucose decreases, induction of galactose utilisation when the system is depleted of glucose and simultaneous use of glucose and galactose as carbon sources when both sugars are present in low concentrations. Preliminary results of the coarse pattern of sugar utilisation were also obtained with a genome-scale E. coli reconstructed network, yielding similar qualitative results.     

Genome-scale reconstruction and analysis of the metabolic network in the hyperthermophilic archaeon sulfolobus solfataricus

We describe the reconstruction of a genome-scale metabolic model of the crenarchaeon Sulfolobus solfataricus, a hyperthermoacidophilic microorganism. It grows in terrestrial volcanic hot springs with growth occurring at pH 2-4 (optimum 3.5) and a temperature of 75-80°C (optimum 80°C). The genome of Sulfolobus solfataricus P2 contains 2,992,245 bp on a single circular chromosome and encodes 2,977 proteins and a number of RNAs. The network comprises 718 metabolic and 58 transport/exchange reactions and 705 unique metabolites, based on the annotated genome and available biochemical data. Using the model in conjunction with constraint-based methods, we simulated the metabolic fluxes induced by different environmental and genetic conditions. The predictions were compared to experimental measurements and phenotypes of S. solfataricus. Furthermore, the performance of the network for 35 different carbon sources known for S. solfataricus from the literature was simulated. Comparing the growth on different carbon sources revealed that glycerol is the carbon source with the highest biomass flux per imported carbon atom (75% higher than glucose). Experimental data was also used to fit the model to phenotypic observations. In addition to the commonly known heterotrophic growth of S. solfataricus, the crenarchaeon is also able to grow autotrophically using the hydroxypropionate-hydroxybutyrate cycle for bicarbonate fixation. We integrated this pathway into our model and compared bicarbonate fixation with growth on glucose as sole carbon source. Finally, we tested the robustness of the metabolism with respect to gene deletions using the method of Minimization of Metabolic Adjustment (MOMA), which predicted that 18% of all possible single gene deletions would be lethal for the organism.     

Constraint-Based Model of Shewanella oneidensis MR-1 Metabolism: A Tool for Data Analysis and Hypothesis Generation

Shewanellae are gram-negative facultatively anaerobic metal-reducing bacteria commonly found in chemically (i.e., redox) stratified environments. Occupying such niches requires the ability to rapidly acclimate to changes in electron donor/acceptor type and availability; hence, the ability to compete and thrive in such environments must ultimately be reflected in the organization and utilization of electron transfer networks, as well as central and peripheral carbon metabolism. To understand how Shewanella oneidensis MR-1 utilizes its resources, the metabolic network was reconstructed. The resulting network consists of 774 reactions, 783 genes, and 634 unique metabolites and contains biosynthesis pathways for all cell constituents. Using constraint-based modeling, we investigated aerobic growth of S. oneidensis MR-1 on numerous carbon sources. To achieve this, we (i) used experimental data to formulate a biomass equation and estimate cellular ATP requirements, (ii) developed an approach to identify cycles (such as futile cycles and circulations), (iii) classified how reaction usage affects cellular growth, (iv) predicted cellular biomass yields on different carbon sources and compared model predictions to experimental measurements, and (v) used experimental results to refine metabolic fluxes for growth on lactate. The results revealed that aerobic lactate-grown cells of S. oneidensis MR-1 used less efficient enzymes to couple electron transport to proton motive force generation, and possibly operated at least one futile cycle involving malic enzymes. Several examples are provided whereby model predictions were validated by experimental data, in particular the role of serine hydroxymethyltransferase and glycine cleavage system in the metabolism of one-carbon units, and growth on different sources of carbon and energy. This work illustrates how integration of computational and experimental efforts facilitates the understanding of microbial metabolism at a systems level.     

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.     

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

【目的】通过挖掘实验性文献,建立巨大芽胞杆菌事实型代谢网络模型,以详尽解析生理特性,优化其生理功能。【方法】从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合成途径和氨基酸代谢途径,并在一定程度上反映了营养底物与基因对巨大芽胞杆菌生长性能的影响。     

Assessment of the metabolic capabilities of Haemophilus influenzae Rd through a genome-scale pathway analysis

The annotated full DNA sequence is becoming available for a growing number of organisms. This information along with additional biochemical and strain-specific data can be used to define metabolic genotypes and reconstruct cellular metabolic networks. The first free-living organism for which the entire genomic sequence was established was Haemophilus influenzae. Its metabolic network is reconstructed herein and contains 461 reactions operating on 367 intracellular and 84 extracellular metabolites. With the metabolic reaction network established, it becomes necessary to determine its underlying pathway structure as defined by the set of extreme pathways. The H. influenzae metabolic network was subdivided into six subsystems and the extreme pathways determined for each subsystem based on stoichiometric, thermodynamic, and systems-specific constraints. Positive linear combinations of these pathways can be taken to determine the extreme pathways for the complete system. Since these pathways span the capabilities of the full system, they could be used to address a number of important physiological questions. First, they were used to reconcile and curate the sequence annotation by identifying reactions whose function was not supported in any of the extreme pathways. Second, they were used to predict gene products that should be co-regulated and perhaps co-expressed. Third, they were used to determine the composition of the minimal substrate requirements needed to support the production of 51 required metabolic products such as amino acids, nucleotides, phospholipids, etc. Fourth, sets of critical gene deletions from core metabolism were determined in the presence of the minimal substrate conditions and in more complete conditions reflecting the environmental niche of H. influenzae in the human host. In the former case, 11 genes were determined to be critical while six remained critical under the latter conditions. This study represents an important milestone in theoretical biology, namely the establishment of the first extreme pathway structure of a whole genome.     

Stoichiometric model of Escherichia coli metabolism: incorporation of growth-rate dependent biomass composition and mechanistic energy requirements

A stoichiometric model of metabolism was developed to describe the balance of metabolic reactions during steady-state growth of Escherichia coli on glucose (or metabolic intermediates) and mineral salts. The model incorporates 153 reversible and 147 irreversible reactions and 289 metabolites from several metabolic data bases for the biosynthesis of the macromolecular precursors, coenzymes, and prosthetic groups necessary for synthesis of all cellular macromolecules. Correlations describing how the cellular composition changes with growth rate were developed from experimental data and were used to calculate the drain of precursors to macromolecules, coenzymes, and prosthetic groups from the metabolic network for the synthesis of those macromolecules at a specific growth rate. Energy requirements for macromolecular polymerization and proofreading, transport of metabolites, and maintenance of transmembrane gradients were included in the model rather than a lumped maintenance energy term. The underdetermined set of equations was solved using the Simplex algorithm, employing realistic objective functions and constraints; the drain of precursors, coenzymes, and prosthetic groups and the energy requirements for the synthesis of macromolecules served as the primary set of constraints. The model accurately predicted experimentally determined metabolic fluxes for aerobic growth on acetate or acetate plus glucose. In addition, the model predicted the genetic and metabolic regulation that must occur for growth under different conditions, such as the opening of the glyoxylate shunt during growth on acetate and the branching of the tricarboxylic acid cycle under anaerobic growth. Sensitivity analyses were performed to determine the flexibility of pathways and the effects of different rates and growth conditions on the distribution of fluxes. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 56: 398-421, 1997.     

Metabolic network reconstruction and flux variability analysis of storage synthesis in developing oilseed rape (Brassica napus L.) embryos

Computational simulation of large‐scale biochemical networks can be used to analyze and predict the metabolic behavior of an organism, such as a developing seed. Based on the biochemical literature, pathways databases and decision rules defining reaction directionality we reconstructed bna572, a stoichiometric metabolic network model representing Brassica napus seed storage metabolism. In the highly compartmentalized network about 25% of the 572 reactions are transport reactions interconnecting nine subcellular compartments and the environment. According to known physiological capabilities of developing B. napus embryos, four nutritional conditions were defined to simulate heterotrophy or photoheterotrophy, each in combination with the availability of inorganic nitrogen (ammonia, nitrate) or amino acids as nitrogen sources. Based on mathematical linear optimization the optimal solution space was comprehensively explored by flux variability analysis, thereby identifying for each reaction the range of flux values allowable under optimality. The range and variability of flux values was then categorized into flux variability types. Across the four nutritional conditions, approximately 13% of the reactions have variable flux values and 10–11% are substitutable (can be inactive), both indicating metabolic redundancy given, for example, by isoenzymes, subcellular compartmentalization or the presence of alternative pathways. About one‐third of the reactions are never used and are associated with pathways that are suboptimal for storage synthesis. Fifty‐seven reactions change flux variability type among the different nutritional conditions, indicating their function in metabolic adjustments. This predictive modeling framework allows analysis and quantitative exploration of storage metabolism of a developing B. napus oilseed.     

Metabolic networks: enzyme function and metabolite structure

Metabolism is one of the most complex cellular processes. Connections between biochemical reactions via substrate and product metabolites create complex metabolic networks that may be analyzed using network theory, stoichiometric analysis, and information on protein structure/function and metabolite properties. These frameworks take into consideration different aspects of enzyme chemistry, enzyme structure and metabolite structure, and demonstrate the impact of metabolic biochemistry on the systemic properties of metabolism. The integration of these approaches and the systematic classification of enzyme function and the chemical structure of metabolites will enhance our understanding of metabolism, and could improve our ability to predict enzyme function and novel metabolic pathways.     

Dynamic flux cartography of hairy roots primary metabolism

A dynamic model for plant cell and hairy root primary metabolism is presented. The model includes nutrient uptake (Pi, sugars, nitrogen sources), the glycolysis and pentose phosphate pathways, the TCA cycle, amino acid biosynthesis, respiratory chain, biosynthesis of cell building blocks (structural hexoses, organic acids, lipids, and organic phosphated molecules). The energy shuttles (ATP, ADP) and cofactors (NAD/H, NADP/H) are also included. The model describes the kinetics of 44 biochemical reactions (fluxes) of the primary metabolism of plant cells and includes 41 biochemical species (metabolites, nutrients, biomass components). Multiple Michaelis-Menten type kinetics are used to describe biochemical reaction rates. Known regulatory phenomena on metabolic pathways are included using sigmoid switch functions. A visualization framework showing fluxes and metabolite concentrations over time is presented. The visualization of fluxes and metabolites is used to analyze simulation results from Catharanthus roseus hairy root 50 d batch cultures. The visualization of the metabolic system allows analyzing split ratios between pathways and flux time-variations. For carbon metabolism, the cells were observed to have relatively high and stable fluxes for the central carbon metabolism and low and variable fluxes for anabolic pathways. For phosphate metabolism, a very high free intracellular Pi turnover rate was observed with higher flux variations than for the carbon metabolism. Nitrogen metabolism also exhibited large flux variations. The potential uses of the model are also discussed.     

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.     

Candidate metabolic network states in human mitochondria. Impact of diabetes, ischemia, and diet

The human mitochondrial metabolic network was recently reconstructed based on proteomic and biochemical data. Linear programming and uniform random sampling were applied herein to identify candidate steady states of the metabolic network that were consistent with the imposed physico-chemical constraints and available experimental data. The activity of the mitochondrion was studied under four metabolic conditions: normal physiologic, diabetic, ischemic, and dietetic. Pairwise correlations between steady-state reaction fluxes were calculated in each condition to evaluate the dependence among the reactions in the network. Applying constraints on exchange fluxes resulted in predictions for intracellular fluxes that agreed with experimental data. Analyses of the steady-state flux distributions showed that the experimentally observed reduced activity of pyruvate dehydrogenase in vivo could be a result of stoichiometric constraints and therefore would not necessarily require enzymatic inhibition. The observed changes in the energy metabolism of the mitochondrion under diabetic conditions were used to evaluate the impact of previously suggested treatments. The results showed that neither normalized glucose uptake nor decreased ketone body uptake have a positive effect on the mitochondrial energy metabolism or network flexibility. Taken together, this study showed that sampling of the steady-state flux space is a powerful method to investigate network properties under different conditions and provides a basis for in silico evaluations of effects of potential disease treatments.     

METANNOGEN: compiling features of biochemical reactions needed for the reconstruction of metabolic networks

Background  

One central goal of computational systems biology is the mathematical modelling of complex metabolic reaction networks. The first and most time-consuming step in the development of such models consists in the stoichiometric reconstruction of the network, i. e. compilation of all metabolites, reactions and transport processes relevant to the considered network and their assignment to the various cellular compartments. Therefore an information system is required to collect and manage data from different databases and scientific literature in order to generate a metabolic network of biochemical reactions that can be subjected to further computational analyses.     

Metabolic flux analysis ofBeijerinckia indica for PS-7 production

In order to investigate central metabolic changes inBeijerinckia indica, cells were grown on different carbon sources and intracellular flux distributions were studied under varying concentrations of nitrogen. Metabolic fluxes were estimated by combining material balances with extracellular substrate uptake rate, biomass formation rate, and exopolysaccharide (EPS) accumulation rate. Thirty-one metabolic reactions and 30 intracellular metabolites were considered for the flux analysis. The results revealed that most of the carbon source was directed into the Entner-Doudoroff pathway, followed by the recycling of triose-3-phosphate back to Hexose-6-phosphate. The pentose phosphate pathway was operated at a minimal level to supply the precursors for biomass formation. The different metabolic behaviors under varying nitrogen concentrations were observed with flux analysis.