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.

A genome-scale metabolic network reconstruction of tomato (Solanum lycopersicum L.) and its application to photorespiratory metabolism.

Tomato (Solanum lycopersicum L.) has been studied extensively due to its high economic value in the market, and high content in health-promoting antioxidant compounds. Tomato is also considered as an excellent model organism for studying the development and metabolism of fleshy fruits. However, the growth, yield and fruit quality of tomatoes can be affected by drought stress, a common abiotic stress for tomato. To investigate the potential metabolic response of tomato plants to drought, we reconstructed iHY3410, a genome-scale metabolic model of tomato leaf, and used this metabolic network to simulate tomato leaf metabolism. The resulting model includes 3410 genes and 2143 biochemical and transport reactions distributed across five intracellular organelles including cytosol, plastid, mitochondrion, peroxisome and vacuole. The model successfully described the known metabolic behaviour of tomato leaf under heterotrophic and phototrophic conditions. The in silico investigation of the metabolic characteristics for photorespiration and other relevant metabolic processes under drought stress suggested that: (i) the flux distributions through the mevalonate (MVA) pathway under drought were distinct from that under normal conditions; and (ii) the changes in fluxes through core metabolic pathways with varying flux ratio of RubisCO carboxylase to oxygenase may contribute to the adaptive stress response of plants. In addition, we improved on previous studies of reaction essentiality analysis for leaf metabolism by including potential alternative routes for compensating reaction knockouts. Altogether, the genome-scale model provides a sound framework for investigating tomato metabolism and gives valuable insights into the functional consequences of abiotic stresses.

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.

Substrate Utilisation and Energy Metabolism in Non-Growing Campylobacter jejuni M1cam.

Campylobacter jejuni, the major cause of bacterial foodborne illness, is also a fastidious organism that requires strict growth requirements in the laboratory. Our aim was to study substrate utilisation and energy metabolism in non-growing C. jejuni to investigate the ability of these bacteria to survive so effectively in the food chain. We integrated phenotypic microarrays and genome-scale metabolic modelling (GSM) to investigate the survival of C. jejuni on 95 substrates. We further investigated the underlying metabolic re-adjustment associated with varying energy demands on each substrate. We identified amino acids, organic acids and H2, as single substrates supporting survival without growth. We identified several different mechanisms, which were used alone or in combination, for ATP production: substrate-level phosphorylation via acetate kinase, the TCA cycle, and oxidative phosphorylation via the electron transport chain that utilised alternative electron donors and acceptors. The benefit of ATP production through each of these mechanisms was associated with the cost of enzyme investment, nutrient availability and/or O2 utilisation. C. jejuni can utilise a wide range of substrates as energy sources, including organic acids commonly used for marination or preservation of ingredients, which might contribute to the success of their survival in changing environments.

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.

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.

Global responses to oxytetracycline treatment in tetracycline-resistant Escherichia coli.

We characterized the global transcriptome of Escherichia coli MG1655:: tetA grown in the presence of ½ MIC (14 mg/L) of OTC, and for comparison WT MG1655 strain grown with 1//2 MIC of OTC (0.25 mg/L OTC). 1646 genes changed expression significantly (FDR > 0.05) in the resistant strain, the majority of which (1246) were also regulated in WT strain. Genes involved in purine synthesis and ribosome structure and function were top-enriched among up-regulated genes, and anaerobic respiration, nitrate metabolism and aromatic amino acid biosynthesis genes among down-regulated genes. Blocking of the purine-synthesis- did not affect resistance phenotypes (MIC and growth rate with OTC), while blocking of protein synthesis using low concentrations of chloramphenicol or gentamicin, lowered MIC towards OTC. Metabolic-modeling, using a novel model for MG1655 and continuous weighing factor that reflected the degree of up or down regulation of genes encoding a reaction, identified 102 metabolic reactions with significant change in flux in MG1655:: tetA when grown in the presence of OTC compared to growth without OTC. These pathways could not have been predicted by simply analyzing functions of the up and down regulated genes, and thus this work has provided a novel method for identification of reactions which are essential in the adaptation to growth in the presence of antimicrobials.

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.

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.

The intertwined metabolism during symbiotic nitrogen fixation elucidated by metabolic modelling.

Genome-scale metabolic network models can be used for various analyses including the prediction of metabolic responses to changes in the environment. Legumes are well known for their rhizobial symbiosis that introduces nitrogen into the global nutrient cycle. Here, we describe a fully compartmentalised, mass and charge-balanced, genome-scale model of the clover Medicago truncatula, which has been adopted as a model organism for legumes. We employed flux balance analysis to demonstrate that the network is capable of producing biomass components in experimentally observed proportions, during day and night. By connecting the plant model to a model of its rhizobial symbiont, Sinorhizobium meliloti, we were able to investigate the effects of the symbiosis on metabolic fluxes and plant growth and could demonstrate how oxygen availability influences metabolic exchanges between plant and symbiont, thus elucidating potential benefits of inter organism amino acid cycling. We thus provide a modelling framework, in which the interlinked metabolism of plants and nodules can be studied from a theoretical perspective.

A Genome Scale Model of Geobacillus thermoglucosidasius (C56-YS93) reveals its biotechnological potential on rice straw hydrolysate.

Rice straw is a major crop residue which is burnt in many countries, creating significant air pollution. Thus, alternative routes for disposal of rice straw are needed. Biotechnological treatment of rice straw hydrolysate has potential to convert this agriculture waste into valuable biofuel(s) and platform chemicals. Geobacillus thermoglucosidasius is a thermophile with properties specially suited for use as a biocatalyst in lignocellulosic bioprocesses, such as high optimal temperature and tolerance to high levels of ethanol. However, the capabilities of G. thermoglucosidasius to utilise sugars in rice straw hydrolysate for making bioethanol and other platform chemicals have not been fully explored. In this work, we have created a genome scale metabolic model (denoted iGT736) of the organism containing 736 gene products, 1159 reactions and 1163 metabolites. The model was validated both by purely theoretical approaches and by comparing the behaviour of the model to previously published experimental results. The model was then used to determine the yields of a variety of platform chemicals from glucose and xylose - two primary sugars in rice straw hydrolysate. A comparison with results from a model of Escherichia coli shows that G. thermoglucosidasius is capable of producing a wider range of products, and that for the products also produced by E. coli, the yields are comparable. We also discuss strategies to utilise arabinose, a minor component of rice straw hydrolysate, and propose additional reactions to lead to the synthesis of xylitol, not currently produced by G. thermoglucosidasius. Our results provide additional motivation for the current exploration of the industrial potential of G. thermoglucosidasius and we make our model publicly available to aid the development of metabolic engineering strategies for this organism.

Metabolic trade-offs between biomass synthesis and photosynthate export at different light intensities in a genome-scale metabolic model of rice.

Previously we have used a genome scale model of rice metabolism to describe how metabolism reconfigures at different light intensities in an expanding leaf of rice. Although this established that the metabolism of the leaf was adequately represented, in the model, the scenario was not that of the typical function of the leaf-to provide material for the rest of the plant. Here we extend our analysis to explore the transition to a source leaf as export of photosynthate increases at the expense of making leaf biomass precursors, again as a function of light intensity. In particular we investigate whether, when the leaf is making a smaller range of compounds for export to the phloem, the same changes occur in the interactions between mitochondrial and chloroplast metabolism as seen in biomass synthesis for growth when light intensity increases. Our results show that the same changes occur qualitatively, though there are slight quantitative differences reflecting differences in the energy and redox requirements for the different metabolic outputs.

Modular decomposition of metabolic systems via null-space analysis.

We describe a method by which the reactions in a metabolic system may be grouped hierarchically into sets of modules to form a metabolic reaction tree. In contrast to previous approaches, the method described here takes into account the fact that, in a viable network, reactions must be capable of sustaining a steady-state flux. In order to achieve this decomposition we introduce a new concept--the reaction correlation coefficient, phi, and show that this is a logical extension of the concept of enzyme (or reaction) subsets. In addition to their application to modular decomposition, reaction correlation coefficients have a number of other interesting properties, including a convenient means for identifying disconnected subnetworks in a system and potential applications to metabolic engineering. The method computes reaction correlation coefficients from an orthonormal basis of the null-space of the stoichiometry matrix. We show that reaction correlation coefficients are uniquely defined, even though the basis of the null-space is not. Once a complete set of reaction correlation coefficients is calculated, a metabolic reaction tree can be determined through the application of standard programming techniques. Computation of the reaction correlation coefficients, and the subsequent construction of the metabolic reaction tree is readily achievable for genome-scale models using a commodity desk-top PC.

Elementary modes analysis of photosynthate metabolism in the chloroplast stroma.

We briefly review the metabolism of the chloroplast stroma, and describe the structural modelling technique of elementary modes analysis. The technique is applied to a model of chloroplast metabolism to investigate viable pathways in the light, in the dark, and in the dark with the addition of sedoheptulose-1,7-bisphosphatase (normally inactive in the dark). The results of the analysis show that it is possible for starch degradation to enhance photosynthetic triose phosphate export in the light, but the reactions of the Calvin cycle alone are not capable of providing a sustainable flux from starch to triose phosphate in the dark. When reactions of the oxidative pentose phosphate pathway are taken into consideration, triose phosphate export in the dark becomes possible by the operation of a cyclic pathway not previously described. The effect of introducing sedoheptulose-1,7-bisphosphatase to the system are relatively minor and, we predict, innocuous in vivo. We conclude that, in contrast with the traditional view of the Calvin cycle and oxidative pentose phosphate pathway as separate systems, they are in fact, in the context of the chloroplast, complementary and overlapping components of the same system.