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

Experimental flux measurements on a network scale

Metabolic flux is a fundamental property of living organisms. In recent years, methods for measuring metabolic flux in plants on a network scale have evolved further. One major challenge in studying flux in plants is the complexity of the plant's metabolism. In particular, in the presence of parallel pathways in multiple cellular compartments, the core of plant central metabolism constitutes a complex network. Hence, a common problem with the reliability of the contemporary results of (13)C-Metabolic Flux Analysis in plants is the substantial reduction in complexity that must be included in the simulated networks; this omission partly is due to limitations in computational simulations. Here, I discuss recent emerging strategies that will better address these shortcomings.

Systems approaches to modelling pathways and networks

It has become commonly accepted that systems approaches to biology are of outstanding importance to gain understanding from the vast amount of data which is presently being generated by advancing high-throughput technologies. The diversity of methods to model pathways and networks has significantly expanded over the past two decades. Modern and traditional approaches are equally important and recent activities aim at integrating the advantages of both. While traditional methods, based on differential equations, are useful to study the dynamics of small systems, modern constraint-based models can be applied to genome-scale systems, but are not able to capture dynamic features. Integrating different approaches is important to develop consistent theoretical descriptions encompassing various scales of biological information. The rapid progress of the field of theoretical systems biology, however, demonstrates how our fundamental theoretical understanding of biology is gaining momentum. The scientific community has apparently accepted the challenge to truly understand the principles of life.

Flux-balance modeling of plant metabolism

Flux-balance modeling of plant metabolic networks provides an important complement to (13)C-based metabolic flux analysis. Flux-balance modeling is a constraints-based approach in which steady-state fluxes in a metabolic network are predicted by using optimization algorithms within an experimentally bounded solution space. In the last 2 years several flux-balance models of plant metabolism have been published including genome-scale models of Arabidopsis metabolism. In this review we consider what has been learnt from these models. In addition, we consider the limitations of flux-balance modeling and identify the main challenges to generating improved and more detailed models of plant metabolism at tissue- and cell-specific scales. Finally we discuss the types of question that flux-balance modeling is well suited to address and its potential role in metabolic engineering and crop improvement.

Computational analysis of storage synthesis in developing Brassica napus L. (oilseed rape) embryos: flux variability analysis in relation to ¹³C metabolic flux analysis

Plant oils are an important renewable resource, and seed oil content is a key agronomical trait that is in part controlled by the metabolic processes within developing seeds. A large-scale model of cellular metabolism in developing embryos of Brassica napus (bna572) was used to predict biomass formation and to analyze metabolic steady states by flux variability analysis under different physiological conditions. Predicted flux patterns are highly correlated with results from prior ¹³C metabolic flux analysis of B. napus developing embryos. Minor differences from the experimental results arose because bna572 always selected only one sugar and one nitrogen source from the available alternatives, and failed to predict the use of the oxidative pentose phosphate pathway. Flux variability, indicative of alternative optimal solutions, revealed alternative pathways that can provide pyruvate and NADPH to plastidic fatty acid synthesis. The nutritional values of different medium substrates were compared based on the overall carbon conversion efficiency (CCE) for the biosynthesis of biomass. Although bna572 has a functional nitrogen assimilation pathway via glutamate synthase, the simulations predict an unexpected role of glycine decarboxylase operating in the direction of NH₄⁺ assimilation. Analysis of the light-dependent improvement of carbon economy predicted two metabolic phases. At very low light levels small reductions in CO₂ efflux can be attributed to enzymes of the tricarboxylic acid cycle (oxoglutarate dehydrogenase, isocitrate dehydrogenase) and glycine decarboxylase. At higher light levels relevant to the ¹³C flux studies, ribulose-1,5-bisphosphate carboxylase activity is predicted to account fully for the light-dependent changes in carbon balance.

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.

Comparison between elementary flux modes analysis and 13C-metabolic fluxes measured in bacterial and plant cells

BACKGROUND

(13)C metabolic flux analysis is one of the pertinent ways to compare two or more physiological states. From a more theoretical standpoint, the structural properties of metabolic networks can be analysed to explore feasible metabolic behaviours and to define the boundaries of steady state flux distributions. Elementary flux mode analysis is one of the most efficient methods for performing this analysis. In this context, recent approaches have tended to compare experimental flux measurements with topological network analysis.

RESULTS

Metabolic networks describing the main pathways of central carbon metabolism were set up for a bacteria species (Corynebacterium glutamicum) and a plant species (Brassica napus) for which experimental flux maps were available. The structural properties of each network were then studied using the concept of elementary flux modes. To do this, coefficients of flux efficiency were calculated for each reaction within the networks by using selected sets of elementary flux modes. Then the relative differences - reflecting the change of substrate i.e. a sugar source for C. glutamicum and a nitrogen source for B. napus - of both flux efficiency and flux measured experimentally were compared. For both organisms, there is a clear relationship between these parameters, thus indicating that the network structure described by the elementary flux modes had captured a significant part of the metabolic activity in both biological systems. In B. napus, the extension of the elementary flux mode analysis to an enlarged metabolic network still resulted in a clear relationship between the change in the coefficients and that of the measured fluxes. Nevertheless, the limitations of the method to fit some particular fluxes are discussed.

CONCLUSION

This consistency between EFM analysis and experimental flux measurements, validated on two metabolic systems allows us to conclude that elementary flux mode analysis could be a useful tool to complement (13)C metabolic flux analysis, by allowing the prediction of changes in internal fluxes before carbon labelling experiments.

The remarkable diversity of plant PEPC (phosphoenolpyruvate carboxylase): recent insights into the physiological functions and post-translational controls of non-photosynthetic PEPCs

PEPC [PEP (phosphoenolpyruvate) carboxylase] is a tightly controlled enzyme located at the core of plant C-metabolism that catalyses the irreversible β-carboxylation of PEP to form oxaloacetate and Pi. The critical role of PEPC in assimilating atmospheric CO2 during C4 and Crassulacean acid metabolism photosynthesis has been studied extensively. PEPC also fulfils a broad spectrum of non-photosynthetic functions, particularly the anaplerotic replenishment of tricarboxylic acid cycle intermediates consumed during biosynthesis and nitrogen assimilation. An impressive array of strategies has evolved to co-ordinate in vivo PEPC activity with cellular demands for C4–C6 carboxylic acids. To achieve its diverse roles and complex regulation, PEPC belongs to a small multigene family encoding several closely related PTPCs (plant-type PEPCs), along with a distantly related BTPC (bacterial-type PEPC). PTPC genes encode ~110-kDa polypeptides containing conserved serine-phosphorylation and lysine-mono-ubiquitination sites, and typically exist as homotetrameric Class-1 PEPCs. In contrast, BTPC genes encode larger ~117-kDa polypeptides owing to a unique intrinsically disordered domain that mediates BTPC's tight interaction with co-expressed PTPC subunits. This association results in the formation of unusual ~900-kDa Class-2 PEPC hetero-octameric complexes that are desensitized to allosteric effectors. BTPC is a catalytic and regulatory subunit of Class-2 PEPC that is subject to multi-site regulatory phosphorylation in vivo. The interaction between divergent PEPC polypeptides within Class-2 PEPCs adds another layer of complexity to the evolution, physiological functions and metabolic control of this essential CO2-fixing plant enzyme. The present review summarizes exciting developments concerning the functions, post-translational controls and subcellular location of plant PTPC and BTPC isoenzymes.

Redefining plant systems biology: from cell to ecosystem

Molecular biologists typically restrict systems biology to cellular levels. By contrast, ecologists define biological systems as communities of interacting individuals at different trophic levels that process energy, nutrient and information flows. Modern plant breeding needs to increase agricultural productivity while decreasing the ecological footprint. This requires a holistic systems biology approach that couples different aggregation levels while considering the variables that affect these biological systems from cell to community. The challenge is to generate accurate experimental data that can be used together with modelling concepts and techniques that allow experimentally verifying in silico predictions. The coupling of aggregation levels in plant sciences, termed Integral Quantification of Biological Organization (IQ(BiO)), might enhance our abilities to generate new desired plant phenotypes.

Metabolomic approaches toward understanding nitrogen metabolism in plants

Plants can assimilate inorganic nitrogen (N) sources to organic N such as amino acids. N is the most important of the mineral nutrients required by plants and its metabolism is tightly coordinated with carbon (C) metabolism in the fundamental processes that permit plant growth. Increased understanding of N regulation may provide important insights for plant growth and improvement of quality of crops and vegetables because N as well as C metabolism are fundamental components of plant life. Metabolomics is a global biochemical approach useful to study N metabolism because metabolites not only reflect the ultimate phenotypes (traits), but can mediate transcript levels as well as protein levels directly and/or indirectly under different N conditions. This review outlines analytical and bioinformatic techniques particularly used to perform metabolomics for studying N metabolism in higher plants. Examples are used to illustrate the application of metabolomic techniques to the model plants Arabidopsis and rice, as well as other crops and vegetables.