A genome-scale metabolic model of Cupriavidus necator H16 integrated with TraDIS and transcriptomic data reveals metabolic insights for biotechnological applications.

Exploiting biological processes to recycle renewable carbon into high value platform chemicals provides a sustainable and greener alternative to current reliance on petrochemicals. In this regard Cupriavidus necator H16 represents a particularly promising microbial chassis due to its ability to grow on a wide range of low-cost feedstocks, including the waste gas carbon dioxide, whilst also naturally producing large quantities of polyhydroxybutyrate (PHB) during nutrient-limited conditions. Understanding the complex metabolic behaviour of this bacterium is a prerequisite for the design of successful engineering strategies for optimising product yields. We present a genome-scale metabolic model (GSM) of C. necator H16 (denoted iCN1361), which is directly constructed from the BioCyc database to improve the readability and reusability of the model. After the initial automated construction, we have performed extensive curation and both theoretical and experimental validation. By carrying out a genome-wide essentiality screening using a Transposon-directed Insertion site Sequencing (TraDIS) approach, we showed that the model could predict gene knockout phenotypes with a high level of accuracy. Importantly, we indicate how experimental and computational predictions can be used to improve model structure and, thus, model accuracy as well as to evaluate potential false positives identified in the experiments. Finally, by integrating transcriptomics data with iCN1361 we create a condition-specific model, which, importantly, better reflects PHB production in C. necator H16. Observed changes in the omics data and in-silico-estimated alterations in fluxes were then used to predict the regulatory control of key cellular processes. The results presented demonstrate that iCN1361 is a valuable tool for unravelling the system-level metabolic behaviour of C. necator H16 and can provide useful insights for designing metabolic engineering strategies.

Overexpression of phospholipid: diacylglycerol acyltransferase in Brassica napus results in changes in lipid metabolism and oil accumulation.

The regulation of lipid metabolism in oil seeds is still not fully understood and increasing our knowledge in this regard is of great economic, as well as intellectual, importance. Oilseed rape (Brassica napus) is a major global oil crop where increases in triacylglycerol (TAG) accumulation have been achieved by overexpression of relevant biosynthetic enzymes. In this study, we expressed Arabidopsis phospholipid: diacylglycerol acyltransferase (PDAT1), one of the two major TAG-forming plant enzymes in B. napus DH12075 to evaluate its effect on lipid metabolism in developing seeds and to estimate its flux control coefficient. Despite several-fold increase in PDAT activity, seeds of three independently generated PDAT transgenic events showed a small but consistent decrease in seed oil content and had altered fatty acid composition of phosphoglycerides and TAG, towards less unsaturation. Mass spectrometry imaging of seed sections confirmed the shift in lipid compositions and indicated that PDAT overexpression altered the distinct heterogeneous distributions of phosphatidylcholine (PC) molecular species. Similar, but less pronounced, changes in TAG molecular species distributions were observed. Our data indicate that PDAT exerts a small, negative, flux control on TAG biosynthesis and could have under-appreciated effects in fine-tuning of B. napus seed lipid composition in a tissue-specific manner. This has important implications for efforts to increase oil accumulation in similar crops.

Stoichiometric analysis of the energetics and metabolic impact of photorespiration in C3 plants.

Analysis of the impact of photorespiration on plant metabolism is usually based on manual inspection of small network diagrams. Here we create a structural metabolic model that contains the reactions that participate in photorespiration in the plastid, peroxisome, mitochondrion and cytosol, and the metabolite exchanges between them. This model was subjected to elementary flux modes analysis, a technique that enumerates all the component, minimal pathways of a network. Any feasible photorespiratory metabolism in the plant will be some combination of the elementary flux modes (EFMs) that contain the Rubisco oxygenase reaction. Amongst the EFMs we obtained was the classic photorespiratory cycle, but there were also modes that involve photorespiration coupled with mitochondrial metabolism and ATP production, the glutathione-ascorbate cycle and nitrate reduction to ammonia. The modes analysis demonstrated the underlying basis of the metabolic linkages with photorespiration that have been inferred experimentally. The set of reactions common to all the elementary modes showed good agreement with the gene products of mutants that have been reported to have a defective phenotype in photorespiratory conditions. Finally, the set of modes provided a formal demonstration that photorespiration itself does not impact on the CO2 :O2 ratio (assimilation quotient), except in those modes associated with concomitant nitrate reduction.

Model-assisted metabolic engineering of Escherichia coli for long chain alkane and alcohol production.

Biologically-derived hydrocarbons are considered to have great potential as next-generation biofuels owing to the similarity of their chemical properties to contemporary diesel and jet fuels. However, the low yield of these hydrocarbons in biotechnological production is a major obstacle for commercialization. Several genetic and process engineering approaches have been adopted to increase the yield of hydrocarbon, but a model driven approach has not been implemented so far. Here, we applied a constraint-based metabolic modeling approach in which a variable demand for alkane biosynthesis was imposed, and co-varying reactions were considered as potential targets for further engineering of an E. coli strain already expressing cyanobacterial enzymes towards higher chain alkane production. The reactions that co-varied with the imposed alkane production were found to be mainly associated with the pentose phosphate pathway (PPP) and the lower half of glycolysis. An optimal modeling solution was achieved by imposing increased flux through the reaction catalyzed by glucose-6-phosphate dehydrogenase (zwf) and iteratively removing 7 reactions from the network, leading to an alkane yield of 94.2% of the theoretical maximum conversion determined by in silico analysis at a given biomass rate. To validate the in silico findings, we first performed pathway optimization of the cyanobacterial enzymes in E. coli via different dosages of genes, promoting substrate channelling through protein fusion and inducing substantial equivalent protein expression, which led to a 36-fold increase in alka(e)ne production from 2.8 mg/L to 102 mg/L. Further, engineering of E. coli based on in silico findings, including biomass constraint, led to an increase in the alka(e)ne titer to 425 mg/L (major components being 249 mg/L pentadecane and 160 mg/L heptadecene), a 148.6-fold improvement over the initial strain, respectively; with a yield of 34.2% of the theoretical maximum. The impact of model-assisted engineering was also tested for the production of long chain fatty alcohol, another commercially important molecule sharing the same pathway while differing only at the terminal reaction, and a titer of 1506 mg/L was achieved with a yield of 86.4% of the theoretical maximum. Moreover, the model assisted engineered strains had produced 2.54 g/L and 12.5 g/L of long chain alkane and fatty alcohol, respectively, in the bioreactor under fed-batch cultivation condition. Our study demonstrated successful implementation of a combined in silico modeling approach along with the pathway and process optimization in achieving the highest reported titers of long chain hydrocarbons in E. coli.

Evolution of central carbon metabolism.

Organisms share a common core to their metabolic networks. But what determined this: chance, chemical necessity, or evolutionary optimization? In this issue of Molecular Cell, Noor et al. (2010) provide new evidence for selection of a network with optimal features from a broader set of possibilities.

A Diel Flux Balance Model Captures Interactions between Light and Dark Metabolism during Day-Night Cycles in C3 and Crassulacean Acid Metabolism Leaves.

Although leaves have to accommodate markedly different metabolic flux patterns in the light and the dark, models of leaf metabolism based on flux-balance analysis (FBA) have so far been confined to consideration of the network under continuous light. An FBA framework is presented that solves the two phases of the diel cycle as a single optimization problem and, thus, provides a more representative model of leaf metabolism. The requirement to support continued export of sugar and amino acids from the leaf during the night and to meet overnight cellular maintenance costs forces the model to set aside stores of both carbon and nitrogen during the day. With only minimal constraints, the model successfully captures many of the known features of C3 leaf metabolism, including the recently discovered role of citrate synthesis and accumulation in the night as a precursor for the provision of carbon skeletons for amino acid synthesis during the day. The diel FBA model can be applied to other temporal separations, such as that which occurs in Crassulacean acid metabolism (CAM) photosynthesis, allowing a system-level analysis of the energetics of CAM. The diel model predicts that there is no overall energetic advantage to CAM, despite the potential for suppression of photorespiration through CO2 concentration. Moreover, any savings in enzyme machinery costs through suppression of photorespiration are likely to be offset by the higher flux demand of the CAM cycle. It is concluded that energetic or nitrogen use considerations are unlikely to be evolutionary drivers for CAM photosynthesis.

A method for accounting for maintenance costs in flux balance analysis improves the prediction of plant cell metabolic phenotypes under stress conditions.

Flux balance models of metabolism generally utilize synthesis of biomass as the main determinant of intracellular fluxes. However, the biomass constraint alone is not sufficient to predict realistic fluxes in central heterotrophic metabolism of plant cells because of the major demand on the energy budget due to transport costs and cell maintenance. This major limitation can be addressed by incorporating transport steps into the metabolic model and by implementing a procedure that uses Pareto optimality analysis to explore the trade-off between ATP and NADPH production for maintenance. This leads to a method for predicting cell maintenance costs on the basis of the measured flux ratio between the oxidative steps of the oxidative pentose phosphate pathway and glycolysis. We show that accounting for transport and maintenance costs substantially improves the accuracy of fluxes predicted from a flux balance model of heterotrophic Arabidopsis cells in culture, irrespective of the objective function used in the analysis. Moreover, when the new method was applied to cells under control, elevated temperature and hyper-osmotic conditions, only elevated temperature led to a substantial increase in cell maintenance costs. It is concluded that the hyper-osmotic conditions tested did not impose a metabolic stress, in as much as the metabolic network is not forced to devote more resources to cell maintenance.

Identification of potential drug targets in Salmonella enterica sv. Typhimurium using metabolic modelling and experimental validation.

Salmonella enterica sv. Typhimurium is an established model organism for Gram-negative, intracellular pathogens. Owing to the rapid spread of resistance to antibiotics among this group of pathogens, new approaches to identify suitable target proteins are required. Based on the genome sequence of S. Typhimurium and associated databases, a genome-scale metabolic model was constructed. Output was based on an experimental determination of the biomass of Salmonella when growing in glucose minimal medium. Linear programming was used to simulate variations in the energy demand while growing in glucose minimal medium. By grouping reactions with similar flux responses, a subnetwork of 34 reactions responding to this variation was identified (the catabolic core). This network was used to identify sets of one and two reactions that when removed from the genome-scale model interfered with energy and biomass generation. Eleven such sets were found to be essential for the production of biomass precursors. Experimental investigation of seven of these showed that knockouts of the associated genes resulted in attenuated growth for four pairs of reactions, whilst three single reactions were shown to be essential for growth.

Responses to light intensity in a genome-scale model of rice metabolism.

We describe the construction and analysis of a genome-scale metabolic model representing a developing leaf cell of rice (Oryza sativa) primarily derived from the annotations in the RiceCyc database. We used flux balance analysis to determine that the model represents a network capable of producing biomass precursors (amino acids, nucleotides, lipid, starch, cellulose, and lignin) in experimentally reported proportions, using carbon dioxide as the sole carbon source. We then repeated the analysis over a range of photon flux values to examine responses in the solutions. The resulting flux distributions show that (1) redox shuttles between the chloroplast, cytosol, and mitochondrion may play a significant role at low light levels, (2) photorespiration can act to dissipate excess energy at high light levels, and (3) the role of mitochondrial metabolism is likely to vary considerably according to the balance between energy demand and availability. It is notable that these organelle interactions, consistent with many experimental observations, arise solely as a result of the need for mass and energy balancing without any explicit assumptions concerning kinetic or other regulatory mechanisms.

Kinetic modelling of the Zymomonas mobilis Entner-Doudoroff pathway: insights into control and functionality.

Zymomonas mobilis, an ethanol-producing bacterium, possesses the Entner-Doudoroff (E-D) pathway, pyruvate decarboxylase and two alcohol dehydrogenase isoenzymes for the fermentative production of ethanol and carbon dioxide from glucose. Using available kinetic parameters, we have developed a kinetic model that incorporates the enzymic reactions of the E-D pathway, both alcohol dehydrogenases, transport reactions and reactions related to ATP metabolism. After optimizing the reaction parameters within likely physiological limits, the resulting kinetic model was capable of simulating glycolysis in vivo and in cell-free extracts with good agreement with the fluxes and steady-state intermediate concentrations reported in previous experimental studies. In addition, the model is shown to be consistent with experimental results for the coupled response of ATP concentration and glycolytic flux to ATPase inhibition. Metabolic control analysis of the model revealed that the majority of flux control resides not inside, but outside the E-D pathway itself, predominantly in ATP consumption, demonstrating why past attempts to increase the glycolytic flux through overexpression of glycolytic enzymes have been unsuccessful. Co-response analysis indicates how homeostasis of ATP concentrations starts to deteriorate markedly at the highest glycolytic rates. This kinetic model has potential for application in Z. mobilis metabolic engineering and, since there are currently no E-D pathway models available in public databases, it can serve as a basis for the development of models for other micro-organisms possessing this type of glycolytic pathway.

Metabolic Control Analysis of triacylglycerol accumulation in oilseed rape.

The increasing global demand for vegetable oils will only be met if there are significant improvements in the productivity of the major oil crops, such as oilseed rape. Metabolic engineering offers the prospect of further gains in yield beyond that already achieved by breeding and selection but requires guidance as to the changes that need to be made. Metabolic Control Analysis, through measurement and estimation of flux control coefficients, can indicate which enzymes have the most influence on a desired flux. Some experiments have previously reported flux control coefficients for oil accumulation in the seeds of oilseed rape, and others have measured control coefficient distributions for multi-enzyme segments of oil synthesis in seed embryo metabolism measured in vitro. In addition, other reported manipulations of oil accumulation contain results that are exploited further here to calculate previously unknown flux control coefficients. These results are then assembled within a framework that allows an integrated interpretation of the controls on oil accumulation from the assimilation of CO2 to deposition of oil in the seed. The analysis shows that the control is distributed to an extent that the gains from amplifying any single target are necessarily limited, but there are candidates for joint amplification that are likely to act synergistically to produce much more significant gains.

A genome-scale metabolic model accurately predicts fluxes in central carbon metabolism under stress conditions.

Flux is a key measure of the metabolic phenotype. Recently, complete (genome-scale) metabolic network models have been established for Arabidopsis (Arabidopsis thaliana), and flux distributions have been predicted using constraints-based modeling and optimization algorithms such as linear programming. While these models are useful for investigating possible flux states under different metabolic scenarios, it is not clear how close the predicted flux distributions are to those occurring in vivo. To address this, fluxes were predicted for heterotrophic Arabidopsis cells and compared with fluxes estimated in parallel by (13)C-metabolic flux analysis (MFA). Reactions of the central carbon metabolic network (glycolysis, the oxidative pentose phosphate pathway, and the tricarboxylic acid [TCA] cycle) were independently analyzed by the two approaches. Net fluxes in glycolysis and the TCA cycle were predicted accurately from the genome-scale model, whereas the oxidative pentose phosphate pathway was poorly predicted. MFA showed that increased temperature and hyperosmotic stress, which altered cell growth, also affected the intracellular flux distribution. Under both conditions, the genome-scale model was able to predict both the direction and magnitude of the changes in flux: namely, increased TCA cycle and decreased phosphoenolpyruvate carboxylase flux at high temperature and a general decrease in fluxes under hyperosmotic stress. MFA also revealed a 3-fold reduction in carbon-use efficiency at the higher temperature. It is concluded that constraints-based genome-scale modeling can be used to predict flux changes in central carbon metabolism under stress conditions.

Can sugars be produced from fatty acids? A test case for pathway analysis tools.

MOTIVATION: In recent years, several methods have been proposed for determining metabolic pathways in an automated way based on network topology. The aim of this work is to analyse these methods by tackling a concrete example relevant in biochemistry. It concerns the question whether even-chain fatty acids, being the most important constituents of lipids, can be converted into sugars at steady state. It was proved five decades ago that this conversion using the Krebs cycle is impossible unless the enzymes of the glyoxylate shunt (or alternative bypasses) are present in the system. Using this example, we can compare the various methods in pathway analysis. RESULTS: Elementary modes analysis (EMA) of a set of enzymes corresponding to the Krebs cycle, glycolysis and gluconeogenesis supports the scientific evidence showing that there is no pathway capable of converting acetyl-CoA to glucose at steady state. This conversion is possible after the addition of isocitrate lyase and malate synthase (forming the glyoxylate shunt) to the system. Dealing with the same example, we compare EMA with two tools based on graph theory available online, PathFinding and Pathway Hunter Tool. These automated network generating tools do not succeed in predicting the conversions known from experiment. They sometimes generate unbalanced paths and reveal problems identifying side metabolites that are not responsible for the carbon net flux. This shows that, for metabolic pathway analysis, it is important to consider the topology (including bimolecular reactions) and stoichiometry of metabolic systems, as is done in EMA.

Building and analysing genome-scale metabolic models.

Reconstructing a model of the metabolic network of an organism from its annotated genome sequence would seem, at first sight, to be one of the most straightforward tasks in functional genomics, even if the various data sources required were never designed with this application in mind. The number of genome-scale metabolic models is, however, lagging far behind the number of sequenced genomes and is likely to continue to do so unless the model-building process can be accelerated. Two aspects that could usefully be improved are the ability to find the sources of error in a nascent model rapidly, and the generation of tenable hypotheses concerning solutions that would improve a model. We will illustrate these issues with approaches we have developed in the course of building metabolic models of Streptococcus agalactiae and Arabidopsis thaliana.

A genome-scale metabolic model of Arabidopsis and some of its properties.

We describe the construction and analysis of a genome-scale metabolic model of Arabidopsis (Arabidopsis thaliana) primarily derived from the annotations in the Aracyc database. We used techniques based on linear programming to demonstrate the following: (1) that the model is capable of producing biomass components (amino acids, nucleotides, lipid, starch, and cellulose) in the proportions observed experimentally in a heterotrophic suspension culture; (2) that approximately only 15% of the available reactions are needed for this purpose and that the size of this network is comparable to estimates of minimal network size for other organisms; (3) that reactions may be grouped according to the changes in flux resulting from a hypothetical stimulus (in this case demand for ATP) and that this allows the identification of potential metabolic modules; and (4) that total ATP demand for growth and maintenance can be inferred and that this is consistent with previous estimates in prokaryotes and yeast.

Detection of stoichiometric inconsistencies in biomolecular models.

MOTIVATION: Metabolic modelling provides a mathematically rigorous basis for system-level analysis of biochemical networks. However, the growing sizes of metabolic models can lead to serious problems in their construction and validation. In this work, we describe a relatively poorly investigated type of modelling error, called stoichiometric inconsistencies. These errors are caused by incorrect definitions of reaction stoichiometries and result in conflicts between two fundamental physical constraints to be satisfied by any valid metabolic model: positivity of molecular masses of all metabolites and mass conservation in all interconversions. RESULTS: We introduce formal definitions of stoichiometric inconsistencies, inconsistent net stoichiometries, elementary leakage modes and other important fundamental properties of incorrectly defined biomolecular networks. Algorithms are described for the verification of stoichiometric consistency of a model, detection of unconserved metabolites and inconsistent minimal net stoichiometries. The usefulness of these algorithms for effective resolving of inconsistencies and for detection of input errors is demonstrated on a published genome-scale metabolic model of Saccharomyces cerevisiae and one of Streptococcus agalactiae constructed using the KEGG database. AVAILABILITY: http://mudshark.brookes.ac.uk/index.php/Albert_Gevorgyan.

Can sugars be produced from fatty acids? A test case for pathway analysis tools.

MOTIVATION: In recent years, several methods have been proposed for determining metabolic pathways in an automated way based on network topology. The aim of this work is to analyse these methods by tackling a concrete example relevant in biochemistry. It concerns the question whether even-chain fatty acids, being the most important constituents of lipids, can be converted into sugars at steady state. It was proved five decades ago that this conversion using the Krebs cycle is impossible unless the enzymes of the glyoxylate shunt (or alternative bypasses) are present in the system. Using this example, we can compare the various methods in pathway analysis. RESULTS: Elementary modes analysis (EMA) of a set of enzymes corresponding to the Krebs cycle, glycolysis and gluconeogenesis supports the scientific evidence showing that there is no pathway capable of converting acetyl-CoA to glucose at steady state. This conversion is possible after the addition of isocitrate lyase and malate synthase (forming the glyoxylate shunt) to the system. Dealing with the same example, we compare EMA with two tools based on graph theory available online, PathFinding and Pathway Hunter Tool. These automated network generating tools do not succeed in predicting the conversions known from experiment. They sometimes generate unbalanced paths and reveal problems identifying side metabolites that are not responsible for the carbon net flux. This shows that, for metabolic pathway analysis, it is important to consider the topology (including bimolecular reactions) and stoichiometry of metabolic systems, as is done in EMA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Contribution of NADH increases to ethanol's inhibition of retinol oxidation by human ADH isoforms.

BACKGROUND: A decrease in retinoic acid levels due to alcohol consumption has been proposed as a contributor to such conditions as fetal alcohol spectrum diseases and ethanol-induced cancers. One molecular mechanism, competitive inhibition by ethanol of the catalytic activity of human alcohol dehydrogenase (EC 1.1.1.1) (ADH) on all-trans-retinol oxidation has been shown for the ADH7 isoform. Ethanol metabolism also causes an increase in the free reduced nicotinamide adenine dinucleotide (NADH) in cells, which might reasonably be expected to decrease the retinol oxidation rate by product inhibition of ADH isoforms. METHODS: To understand the relative importance of these two mechanisms by which ethanol decreases the retinol oxidation in vivo we need to assess them quantitatively. We have built a model system of 4 reactions: (1) ADH oxidation of ethanol and NAD(+), (2) ADH oxidation of retinol and NAD(+), (3) oxidation of ethanol by a generalized Ethanol(oxidase) that uses NAD(+), (4) NADH(oxidase) which carries out NADH turnover. RESULTS: Using the metabolic modeling package ScrumPy, we have shown that the ethanol-induced increase in NADH contributes from 0% to 90% of the inhibition by ethanol, depending on (ethanol) and ADH isoform. Furthermore, while the majority of flux control of retinaldehyde production is exerted by ADH, Ethanol(oxidase) and the NADH(oxidase) contribute as well. CONCLUSIONS: Our results show that the ethanol-induced increase in NADH makes a contribution of comparable importance to the ethanol competitive inhibition throughout the range of conditions likely to occur in vivo, and must be considered in the assessment of the in vivo mechanism of ethanol interference with fetal development and other diseases.

Metabolic pathways in the post-genome era.

Metabolic pathways are a central paradigm in biology. Historically, they have been defined on the basis of their step-by-step discovery. However, the genome-scale metabolic networks now being reconstructed from annotation of genome sequences demand new network-based definitions of pathways to facilitate analysis of their capabilities and functions, such as metabolic versatility and robustness, and optimal growth rates. This demand has led to the development of a new mathematically based analysis of complex, metabolic networks that enumerates all their unique pathways that take into account all requirements for cofactors and byproducts. Applications include the design of engineered biological systems, the generation of testable hypotheses regarding network structure and function, and the elucidation of properties that can not be described by simple descriptions of individual components (such as product yield, network robustness, correlated reactions and predictions of minimal media). Recently, these properties have also been studied in genome-scale networks. Thus, network-based pathways are emerging as an important paradigm for analysis of biological systems.