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.

Assessing Metabolic Flux in Plants with Radiorespirometry.

Carbohydrates are the dominant respiratory substrate in many plant cells. However, the route of carbohydrate oxidation varies depending on the relative cellular demands for energy, reductant, and precursors for biosynthesis. During these processes individual substrate carbon atoms are differentially released as carbon dioxide by specific reactions in the network, and this can be measured by monitoring the release of 14CO2 from a range of positionally labeled forms of [14C]glucose. Although the relative amounts of carbon dioxide produced from different carbon positions do not allow precise determination of fluxes, they are indicative of the route of carbohydrate utilization. Such information can be used to determine whether a comprehensive metabolic flux analysis is merited, and also to facilitate independent verification of flux maps generated by other techniques. This chapter describes an approach to determine and interpret the pattern of oxidation of carbohydrates by monitoring 14CO2 release during metabolism of exogenously supplied [1-14C]-, [2-14C]-, [3,4-14C]-, and [6-14C]glucose. The method is exemplified by studies on Arabidopsis cell suspension cultures, but the protocol can be easily adapted for the investigation of other plant materials.

Lipogenesis and Redox Balance in Nitrogen-Fixing Pea Bacteroids.

UNLABELLED: Within legume root nodules, rhizobia differentiate into bacteroids that oxidize host-derived dicarboxylic acids, which is assumed to occur via the tricarboxylic acid (TCA) cycle to generate NAD(P)H for reduction of N2 Metabolic flux analysis of laboratory-grown Rhizobium leguminosarum showed that the flux from [(13)C]succinate was consistent with respiration of an obligate aerobe growing on a TCA cycle intermediate as the sole carbon source. However, the instability of fragile pea bacteroids prevented their steady-state labeling under N2-fixing conditions. Therefore, comparative metabolomic profiling was used to compare free-living R. leguminosarum with pea bacteroids. While the TCA cycle was shown to be essential for maximal rates of N2 fixation, levels of pyruvate (5.5-fold reduced), acetyl coenzyme A (acetyl-CoA; 50-fold reduced), free coenzyme A (33-fold reduced), and citrate (4.5-fold reduced) were much lower in bacteroids. Instead of completely oxidizing acetyl-CoA, pea bacteroids channel it into both lipid and the lipid-like polymer poly-ß-hydroxybutyrate (PHB), the latter via a type III PHB synthase that is active only in bacteroids. Lipogenesis may be a fundamental requirement of the redox poise of electron donation to N2 in all legume nodules. Direct reduction by NAD(P)H of the likely electron donors for nitrogenase, such as ferredoxin, is inconsistent with their redox potentials. Instead, bacteroids must balance the production of NAD(P)H from oxidation of acetyl-CoA in the TCA cycle with its storage in PHB and lipids. IMPORTANCE: Biological nitrogen fixation by symbiotic bacteria (rhizobia) in legume root nodules is an energy-expensive process. Within legume root nodules, rhizobia differentiate into bacteroids that oxidize host-derived dicarboxylic acids, which is assumed to occur via the TCA cycle to generate NAD(P)H for reduction of N2 However, direct reduction of the likely electron donors for nitrogenase, such as ferredoxin, is inconsistent with their redox potentials. Instead, bacteroids must balance oxidation of plant-derived dicarboxylates in the TCA cycle with lipid synthesis. Pea bacteroids channel acetyl-CoA into both lipid and the lipid-like polymer poly-ß-hydroxybutyrate, the latter via a type II PHB synthase. Lipogenesis is likely to be a fundamental requirement of the redox poise of electron donation to N2 in all legume nodules.

The metabolic flux phenotype of heterotrophic Arabidopsis cells reveals a flexible balance between the cytosolic and plastidic contributions to carbohydrate oxidation in response to phosphate limitation.

Understanding the mechanisms that allow plants to respond to variable and reduced availability of inorganic phosphate is of increasing agricultural importance because of the continuing depletion of the rock phosphate reserves that are used to combat inadequate phosphate levels in the soil. Changes in gene expression, protein levels, enzyme activities and metabolite levels all point to a reconfiguration of the central metabolic network in response to reduced availability of inorganic phosphate, but the metabolic significance of these changes can only be assessed in terms of the fluxes supported by the network. Steady-state metabolic flux analysis was used to define the metabolic phenotype of a heterotrophic Arabidopsis thaliana cell culture grown on a Murashige and Skoog medium containing 0, 1.25 or 5 mm inorganic phosphate. Fluxes through the central metabolic network were deduced from the redistribution of (13) C into metabolic intermediates and end products when cells were labelled with [1-(13) C], [2-(13) C], or [(13) C6 ]glucose, in combination with (14) C measurements of the rates of biomass accumulation. Analysis of the flux maps showed that reduced levels of phosphate in the growth medium stimulated flux through phosphoenolpyruvate carboxylase and malic enzyme, altered the balance between cytosolic and plastidic carbohydrate oxidation in favour of the plastid, and increased cell maintenance costs. We argue that plant cells respond to phosphate deprivation by reconfiguring the flux distribution through the pathways of carbohydrate oxidation to take advantage of better phosphate homeostasis in the plastid.

Metabolic flux phenotype of tobacco hairy roots engineered for increased geraniol production.

The goal of this study was to characterise the metabolic flux phenotype of transgenic tobacco (Nicotiana tabacum) hairy roots engineered for increased biosynthesis of geraniol, an intermediate of the terpenoid indole alkaloid pathway. Steady state, stable isotope labelling was used to determine flux maps of central carbon metabolism for transgenic lines over-expressing (i) plastid-targeted geraniol synthase (pGES) from Valeriana officinalis, and (ii) pGES in combination with plastid-targeted geranyl pyrophosphate synthase from Arabidopsis thaliana (pGES+pGPPS), as well as for wild type and control-vector-transformed roots. Fluxes were constrained by the redistribution of label from [1-¹³C]-, [2-¹³C]- or [¹³C6]glucose into amino acids, sugars and organic acids at isotopic steady state, and by biomass output fluxes determined from the fractionation of [U-¹4C]glucose into insoluble polymers. No significant differences in growth and biomass composition were observed between the lines. The pGES line accumulated significant amounts of geraniol/geraniol glycosides (151±24 ng/mg dry weight) and the de novo synthesis of geraniol in pGES was confirmed by ¹³C labelling analysis. The pGES+pGPPS also accumulated geraniol and geraniol glycosides, but to lower levels than the pGES line. Although there was a distinct impact of the transgenes at the level of geraniol synthesis, other network fluxes were unaffected, reflecting the capacity of central metabolism to meet the relatively modest demand for increased precursors in the transgenic lines. It is concluded that re-engineering of the terpenoid indole alkaloid pathway will only require simultaneous manipulation of the steps producing the pathway precursors that originate in central metabolism in tissues engineered to produce at least an order of magnitude more geraniol than has been achieved so far.

Quantification of ¹³C enrichments and isotopomer abundances for metabolic flux analysis using 1D NMR spectroscopy.

The analysis of stable isotope incorporation following feeding of (13)C-labeled precursors to plant tissues provides the constraints necessary for metabolic flux analysis. This protocol describes the use of one-dimensional (1)H and (13)C nuclear magnetic resonance spectroscopy for the quantification of (13)C enrichments and isotopomer abundances in mixtures of metabolites or hydrolyzed biomass components.

Optimization of steady-state ¹³C-labeling experiments for metabolic flux analysis.

While steady-state (13)C metabolic flux analysis is a powerful method for deducing multiple fluxes in the central metabolic network of heterotrophic and mixotrophic plant tissues, it is also time-consuming and technically challenging. Key steps in the design and interpretation of steady-state (13)C labeling experiments are illustrated with a generic protocol based on applications to plant cell suspension cultures.

The metabolic flux phenotype of heterotrophic Arabidopsis cells reveals a complex response to changes in nitrogen supply.

The extent to which individual plants utilise nitrate and ammonium, the two principal nitrogen sources in the rhizosphere, is variable and many species require a balance between the two forms for optimal growth. The effects of nitrate and ammonium on gene expression, enzyme activity and metabolite composition have been documented extensively with the aim of understanding the way in which plant cells respond to the different forms of nitrogen, but ultimately the impact of these changes on the organisation and operation of the central metabolic network can only be addressed by analysing the fluxes supported by the network. Accordingly steady-state metabolic flux analysis was used to define the metabolic phenotype of a heterotrophic Arabidopsis thaliana cell culture grown in Murashige and Skoog and ammonium-free media, treatments that influenced growth and biomass composition. Fluxes through the central metabolic network were deduced from the redistribution of label into metabolic intermediates and end products observed when cells were labelled with [1-(13) C]-, [2-(13) C]- or [(13) C6 ]glucose, in tandem with (14) C-measurements of the net accumulation of biomass. Analysis of the flux maps showed that: (i) flux through the oxidative pentose phosphate pathway varied independently of the reductant demand for biosynthesis, (ii) non-plastidic processes made a significant and variable contribution to the provision of reducing power for the plastid, and (iii) the inclusion of ammonium in the growth medium increased cell maintenance costs, in agreement with the futile cycling model of ammonium toxicity. These conclusions highlight the complexity of the metabolic response to a change in nitrogen nutrition.

Subcellular flux analysis of central metabolism in a heterotrophic Arabidopsis cell suspension using steady-state stable isotope labeling.

The presence of cytosolic and plastidic pathways of carbohydrate oxidation is a characteristic feature of plant cell metabolism. Ideally, steady-state metabolic flux analysis, an emerging tool for creating flux maps of heterotrophic plant metabolism, would capture this feature of the metabolic phenotype, but the extent to which this can be achieved is uncertain. To address this question, fluxes through the pathways of central metabolism in a heterotrophic Arabidopsis (Arabidopsis thaliana) cell suspension culture were deduced from the redistribution of label in steady-state (13)C-labeling experiments using [1-(13)C]-, [2-(13)C]-, and [U-(13)C(6)]glucose. Focusing on the pentose phosphate pathway (PPP), multiple data sets were fitted simultaneously to models in which the subcellular compartmentation of the PPP was altered. The observed redistribution of the label could be explained by any one of three models of the subcellular compartmentation of the oxidative PPP, but other biochemical evidence favored the model in which the oxidative steps of the PPP were duplicated in the cytosol and plastids, with flux through these reactions occurring largely in the cytosol. The analysis emphasizes the inherent difficulty of analyzing the PPP without predefining the extent of its compartmentation and the importance of obtaining high-quality data that report directly on specific subcellular processes. The Arabidopsis flux map also shows that the potential ATP yield of respiration in heterotrophic plant cells can greatly exceed the direct metabolic requirements for biosynthesis, highlighting the need for caution when predicting flux through metabolic networks using assumptions based on the energetics of resource utilization.

Metabolic network fluxes in heterotrophic Arabidopsis cells: stability of the flux distribution under different oxygenation conditions.

Steady-state labeling experiments with [1-(13)C]Glc were used to measure multiple metabolic fluxes through the pathways of central metabolism in a heterotrophic cell suspension culture of Arabidopsis (Arabidopsis thaliana). The protocol was based on in silico modeling to establish the optimal labeled precursor, validation of the isotopic and metabolic steady state, extensive nuclear magnetic resonance analysis of the redistribution of label into soluble metabolites, starch, and protein, and a comprehensive set of biomass measurements. Following a simple modification of the cell culture procedure, cells were grown at two oxygen concentrations, and flux maps of central metabolism were constructed on the basis of replicated experiments and rigorous statistical analysis. Increased growth rate at the higher O(2) concentration was associated with an increase in fluxes throughout the network, and this was achieved without any significant change in relative fluxes despite differences in the metabolite profile of organic acids, amino acids, and carbohydrates. The balance between biosynthesis and respiration within the tricarboxylic acid cycle was unchanged, with 38% +/- 5% of carbon entering used for biosynthesis under standard O(2) conditions and 33% +/- 2% under elevated O(2). These results add to the emerging picture of the stability of the central metabolic network and its capacity to respond to physiological perturbations with the minimum of rearrangement. The lack of correlation between the change in metabolite profile, which implied significant disruption of the metabolic network following the alteration in the oxygen supply, and the unchanging flux distribution highlights a potential difficulty in the interpretation of metabolomic data.