Metabolic Control Analysis of glycolysis in tuber tissue of potato (Solanum tuberosum): explanation for the low control coefficient of phosphofructokinase over respiratory flux

A new selective force driving metabolic gene clustering

The evolution of operons has puzzled evolutionary biologists since their discovery, and many theories exist to explain their emergence, spreading, and evolutionary conservation. In this work, we suggest that DNA replication introduces a selective force for the clustering of functionally related genes on chromosomes, which we interpret as a preliminary and necessary step in operon formation. Our reasoning starts from the observation that DNA replication produces copy number variations of genomic regions, and we propose that such changes perturb metabolism. The formalization of this effect by exploiting concepts from metabolic control analysis suggests that the minimization of such perturbations during evolution could be achieved through the formation of gene clusters and operons. We support our theoretical derivations with simulations based on a realistic metabolic network, and we confirm that present-day genomes have a degree of compaction of functionally related genes, which is significantly correlated to the proposed perturbations introduced by replication. The formation of clusters of functionally related genes in microbial genomes has puzzled microbiologists since their first discovery. Here, we suggest that replication, and the copy number variations due to the replisome passage, might play a role in the process through a perturbation in metabolite homeostasis. We provide theoretical support to this hypothesis, and we found that both simulations and genomic analysis support our hypothesis.

IMPORTANCE

The formation of clusters of functionally related genes in microbial genomes has puzzled microbiologists since their discovery. Here, we suggest that replication, and the copy number variations due to the replisome passage, might play a role in the process through a perturbation in metabolite homeostasis. We provide theoretical support to this hypothesis, and we found that both simulations and genomic analysis support our hypothesis.

Enzyme-based kinetic modelling of ASC-GSH cycle during tomato fruit development reveals the importance of reducing power and ROS availability

The ascorbate-glutathione (ASC-GSH) cycle is at the heart of redox metabolism, linking the major redox buffers with central metabolism through the processing of reactive oxygen species (ROS) and pyridine nucleotide metabolism. Tomato fruit development is underpinned by changes in redox buffer contents and their associated enzyme capacities, but interactions between them remain unclear. Based on quantitative data obtained for the core redox metabolism, we built an enzyme-based kinetic model to calculate redox metabolite concentrations with their corresponding fluxes and control coefficients. Dynamic and associated regulations of the ASC-GSH cycle throughout the whole fruit development were analysed and pointed to a sequential metabolic control of redox fluxes by ASC synthesis, NAD(P)H and ROS availability depending on the developmental phase. Furthermore, we highlighted that monodehydroascorbate reductase and the availability of reducing power were found to be the main regulators of the redox state of ASC and GSH during fruit growth under optimal conditions. Our kinetic modelling approach indicated that tomato fruit development displayed growth phase-dependent redox metabolism linked with central metabolism via pyridine nucleotides and H2 O2 availability, while providing a new tool to the scientific community to investigate redox metabolism in fruits.

Delving deeper: Relating the behaviour of a metabolic system to the properties of its components using symbolic metabolic control analysis

High-level behaviour of metabolic systems results from the properties of, and interactions between, numerous molecular components. Reaching a complete understanding of metabolic behaviour based on the system’s components is therefore a difficult task. This problem can be tackled by constructing and subsequently analysing kinetic models of metabolic pathways since such models aim to capture all the relevant properties of the system components and their interactions. Symbolic control analysis is a framework for analysing pathway models in order to reach a mechanistic understanding of their behaviour. By providing algebraic expressions for the sensitivities of system properties, such as metabolic flux or steady-state concentrations, in terms of the properties of individual reactions it allows one to trace the high level behaviour back to these low level components. Here we apply this method to a model of pyruvate branch metabolism in Lactococcus lactis in order to explain a previously observed negative flux response towards an increase in substrate concentration. With this method we are able to show, first, that the sensitivity of flux towards changes in reaction rates (represented by flux control coefficients) is determined by the individual metabolic branches of the pathway, and second, how the sensitivities of individual reaction rates towards their substrates (represented by elasticity coefficients) contribute to this flux control. We also quantify the contributions of enzyme binding and mass-action to enzyme elasticity separately, which allows for an even finer-grained understanding of flux control. These analytical tools allow us to analyse the control properties of a metabolic model and to arrive at a mechanistic understanding of the quantitative contributions of each of the enzymes to this control. Our analysis provides an example of the descriptive power of the general principles of symbolic control analysis.

Control and regulation of pathways via negative feedback

The biochemical networks found in living organisms include a huge variety of control mechanisms at multiple levels of organization. While the mechanistic and molecular details of many of these control mechanisms are understood, their exact role in driving cellular behaviour is not. For example, yeast glycolysis has been studied for almost 80 years but it is only recently that we have come to understand the systemic role of the multitude of feedback and feed-forward controls that exist in this pathway. In this article, control theory is discussed as an approach to dissect the control logic of complex pathways. One of the key issues is distinguishing between the terms control and regulation and how these concepts are applied to regulated enzymes such as phosphofructokinase. In doing so, one of the paradoxes in metabolic regulation can be resolved where enzymes such as phosphofructokinase have little control but, nevertheless, possess significant regulatory influence.

Quantitative Multilevel Analysis of Central Metabolism in Developing Oilseeds of Oilseed Rape during in Vitro Culture

Seeds provide the basis for many food, feed, and fuel products. Continued increases in seed yield, composition, and quality require an improved understanding of how the developing seed converts carbon and nitrogen supplies into storage. Current knowledge of this process is often based on the premise that transcriptional regulation directly translates via enzyme concentration into flux. In an attempt to highlight metabolic control, we explore genotypic differences in carbon partitioning for in vitro cultured developing embryos of oilseed rape (Brassica napus). We determined biomass composition as well as 79 net fluxes, the levels of 77 metabolites, and 26 enzyme activities with specific focus on central metabolism in nine selected germplasm accessions. Overall, we observed a tradeoff between the biomass component fractions of lipid and starch. With increasing lipid content over the spectrum of genotypes, plastidic fatty acid synthesis and glycolytic flux increased concomitantly, while glycolytic intermediates decreased. The lipid/starch tradeoff was not reflected at the proteome level, pointing to the significance of (posttranslational) metabolic control. Enzyme activity/flux and metabolite/flux correlations suggest that plastidic pyruvate kinase exerts flux control and that the lipid/starch tradeoff is most likely mediated by allosteric feedback regulation of phosphofructokinase and ADP-glucose pyrophosphorylase. Quantitative data were also used to calculate in vivo mass action ratios, reaction equilibria, and metabolite turnover times. Compounds like cyclic 3',5'-AMP and sucrose-6-phosphate were identified to potentially be involved in so far unknown mechanisms of metabolic control. This study provides a rich source of quantitative data for those studying central metabolism.

Identification of the expressed protein and the impact of change in ascorbate peroxidase activity related to endodormancy breaking in Pyrus pyrifolia

Endodormancy is an important feature of perennial deciduous fruit trees that survive in the extreme climates brought about by seasonal variation. To acquire a comprehensive knowledge of the biochemical processes occurring just before endodormancy breaking, the buds collected in the pre-breaking period (PP) phase were used as samples to identify the proteins related to the breaking of endodormancy in the Japanese pear (Pyrus pyrifolia Nakai). Using nano-ESI-LC-MS/MS analysis, 96 proteins were overlapped by analyses of three times and identified as expressed proteins at the PP stage. Among these proteins, dehydrin, several classes of heat shock proteins (HSP), auxin-binding protein, and auxin-induced protein were identified in the floral bud in the PP stage. The majority of these proteins were involved primarily in the oxidation-reduction process. We focused on catalase (CAT), peroxidase (POD), and ascorbate peroxidase (APX) as enzymes regulating the levels of hydrogen peroxide (H2O2) in the bud. From measurements taken during the deepest period (DP), PP, mid-breaking period (MP), and late-breaking period (LP) of endodormancy, CAT activity decreased gradually, while APX activity also decreased from DP to MP, but then increased rapidly during LP. Protein data for PP and the rapid increase in APX activity observed in LP provided knowledge of the biochemical processes that regulate the consecutive transition from endodormancy breaking to ecodormancy induction in the Japanese pear.

Model-assisted analysis of sugar metabolism throughout tomato fruit development reveals enzyme and carrier properties in relation to vacuole expansion

A kinetic model combining enzyme activity measurements and subcellular compartmentation was parameterized to fit the sucrose, hexose, and glucose-6-P contents of pericarp throughout tomato (Solanum lycopersicum) fruit development. The model was further validated using independent data obtained from domesticated and wild tomato species and on transgenic lines. A hierarchical clustering analysis of the calculated fluxes and enzyme capacities together revealed stage-dependent features. Cell division was characterized by a high sucrolytic activity of the vacuole, whereas sucrose cleavage during expansion was sustained by both sucrose synthase and neutral invertase, associated with minimal futile cycling. Most importantly, a tight correlation between flux rate and enzyme capacity was found for fructokinase and PPi-dependent phosphofructokinase during cell division and for sucrose synthase, UDP-glucopyrophosphorylase, and phosphoglucomutase during expansion, thus suggesting an adaptation of enzyme abundance to metabolic needs. In contrast, for most enzymes, flux rates varied irrespectively of enzyme capacities, and most enzymes functioned at <5% of their maximal catalytic capacity. One of the major findings with the model was the high accumulation of soluble sugars within the vacuole together with organic acids, thus enabling the osmotic-driven vacuole expansion that was found during cell division.

Label-free quantitative proteomics analysis of dormant terminal buds of poplar

Induction and break of bud dormancy are important features for perennial plants surviving extreme seasonal variations in climate. However, the molecular mechanism of the dormancy regulation, still remain poorly understood. To better understand the molecular basis of poplar bud dormancy, we used a label-free quantitative proteomics method based on nanoscale ultra performance liquid chromatography-ESI-MS(E) for investigation of differential protein expression during dormancy induction, dormancy, and dormancy break in apical buds of poplar (Populus simonii × P. nigra). Among these identified over 300 proteins during poplar bud dormancy, there are 74 significantly altered proteins, most of which involved in carbohydrate metabolism (22 %), redox regulation (19 %), amino acid transport and metabolism (10 %), and stress response (8 %). Thirty-one of these proteins were up-regulated, five were down-regulated during three phase, and thirty-eight were expressed specifically under different conditions. Pathway analysis suggests that there are still the presence of various physiological activities and a particular influence on photosynthesis and energy metabolism during poplar bud dormancy. Differential expression patterns were identified for key enzymes involved in major metabolic pathways such as glycolysis and the pentose phosphate pathway, thus manifesting the interplay of intricate molecular events in energy generation for new protein synthesis in the dormant buds. Furthermore, there are significant changes present in redox regulation and defense response proteins, for instance in peroxidase and ascorbate peroxidase. Overall, this study provides a better understanding of the possible regulation mechanisms during poplar bud dormancy.

Metabolic control and regulation of the tricarboxylic acid cycle in photosynthetic and heterotrophic plant tissues

The tricarboxylic acid (TCA) cycle is a crucial component of respiratory metabolism in both photosynthetic and heterotrophic plant organs. All of the major genes of the tomato TCA cycle have been cloned recently, allowing the generation of a suite of transgenic plants in which the majority of the enzymes in the pathway are progressively decreased. Investigations of these plants have provided an almost complete view of the distribution of control in this important pathway. Our studies suggest that citrate synthase, aconitase, isocitrate dehydrogenase, succinyl CoA ligase, succinate dehydrogenase, fumarase and malate dehydrogenase have control coefficients flux for respiration of -0.4, 0.964, -0.123, 0.0008, 0.289, 0.601 and 1.76, respectively; while 2-oxoglutarate dehydrogenase is estimated to have a control coefficient of 0.786 in potato tubers. These results thus indicate that the control of this pathway is distributed among malate dehydrogenase, aconitase, fumarase, succinate dehydrogenase and 2-oxoglutarate dehydrogenase. The unusual distribution of control estimated here is consistent with specific non-cyclic flux mode and cytosolic bypasses that operate in illuminated leaves. These observations are discussed in the context of known regulatory properties of the enzymes and some illustrative examples of how the pathway responds to environmental change are given.

Adaptation to HIF-1 deficiency by upregulation of the AMP/ATP ratio and phosphofructokinase activation in hepatomas

BACKGROUND

HIF-1 deficiency has marked effects on tumour glycolysis and growth. We therefore investigated the consequences of HIF-1 deficiency in mice, using the well established Hepa-1 wild-type (WT) and HIF-1β-deficient (c4) model. These mechanisms could be clinically relevant, since HIF-1 is now a therapeutic target.

METHODS

Hepa-1 WT and c4 tumours grown in vivo were analysed by 18FDG-PET and 19FDG Magnetic Resonance Spectroscopy for glucose uptake; by HPLC for adenine nucleotides; by immunohistochemistry for GLUTs; by immunoblotting and by DIGE followed by tandem mass spectrometry for protein expression; and by classical enzymatic methods for enzyme activity.

RESULTS

HIF-1β deficient Hepa-1 c4 tumours grew significantly more slowly than WT tumours, and (as expected) showed significantly lower expression of many glycolytic enzymes. However, HIF-1β deficiency caused no significant change in the rate of glucose uptake in c4 tumours compared to WT when assessed in vivo by measuring fluoro-deoxyglucose (FDG) uptake. Immunohistochemistry demonstrated less GLUT-1 in c4 tumours, whereas GLUT-2 (liver type) was similar to WT. Factors that might upregulate glucose uptake independently of HIF-1 (phospho-Akt, c-Myc) were shown to have either lower or similar expression in c4 compared to WT tumours. However the AMP/ATP ratio was 4.5 fold higher (p < 0.01) in c4 tumours, and phosphofructokinase-1 (PFK-1) activity, measured at prevailing cellular ATP and AMP concentrations, was up to two-fold higher in homogenates of the deficient c4 cells and tumours compared to WT (p < 0.001), suggesting that allosteric PFK activation could explain their normal level of glycolysis. Phospho AMP-Kinase was also higher in the c4 tumours.

CONCLUSIONS

Despite their defective HIF-1 and consequent down-regulation of glycolytic enzyme expression, Hepa-1 c4 tumours maintain glucose uptake and glycolysis because the resulting low [ATP] high [AMP] allosterically activate PFK-1. This mechanism of resistance would keep glycolysis functioning and also result in activation of AMP-Kinase and growth inhibition; it may have major implications for the therapeutic activity of HIF inhibitors in vivo. Interestingly, this control mechanism does not involve transcriptional control or proteomics, but rather the classical activation and inhibition mechanisms of glycolytic enzymes.

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.

Switch between life history strategies due to changes in glycolytic enzyme gene dosage in Saccharomyces cerevisiae

Adaptation is the process whereby a population or species becomes better fitted to its habitat through modifications of various life history traits which can be positively or negatively correlated. The molecular factors underlying these covariations remain to be elucidated. Using Saccharomyces cerevisiae as a model system, we have investigated the effects on life history traits of varying the dosage of genes involved in the transformation of resources into energy. Changing gene dosage for each of three glycolytic enzyme genes (hexokinase 2, phosphoglucose isomerase, and fructose-1,6-bisphosphate aldolase) resulted in variation in enzyme activities, glucose consumption rate, and life history traits (growth rate, carrying capacity, and cell size). However, the range of effects depended on which enzyme was expressed differently. Most interestingly, these changes revealed a genetic trade-off between carrying capacity and cell size, supporting the discovery of two extreme life history strategies already described in yeast populations: the "ants," which have lower glycolytic gene dosage, take up glucose slowly, and have a small cell size but reach a high carrying capacity, and the "grasshoppers," which have higher glycolytic gene dosage, consume glucose more rapidly, and allocate it to a larger cell size but reach a lower carrying capacity. These results demonstrate antagonist pleiotropy for glycolytic genes and show that altered dosage of a single gene drives a switch between two life history strategies in yeast.

Combining genomics, metabolome analysis, and biochemical modelling to understand metabolic networks

Now that complete genome sequences are available for a variety of organisms, the elucidation of gene functions involved in metabolism necessarily includes a better understanding of cellular responses upon mutations on all levels of gene products, mRNA, proteins, and metabolites. Such progress is essential since the observable properties of organisms – the phenotypes – are produced by the genotype in juxtaposition with the environment. Whereas much has been done to make mRNA and protein profiling possible, considerably less effort has been put into profiling the end products of gene expression, metabolites. To date, analytical approaches have been aimed primarily at the accurate quantification of a number of pre-defined target metabolites, or at producing fingerprints of metabolic changes without individually determining metabolite identities. Neither of these approaches allows the formation of an in-depth understanding of the biochemical behaviour within metabolic networks. Yet, by carefully choosing protocols for sample preparation and analytical techniques, a number of chemically different classes of compounds can be quantified simultaneously to enable such understanding. In this review, the terms describing various metabolite-oriented approaches are given, and the differences among these approaches are outlined. Metabolite target analysis, metabolite profiling, metabolomics, and metabolic fingerprinting are considered. For each approach, a number of examples are given, and potential applications are discussed.

A kinetic model describes metabolic response to perturbations and distribution of flux control in the benzenoid network of Petunia hybrida

In recent years there has been much interest in the genetic enhancement of plant metabolism; however, attempts at genetic modification are often unsuccessful due to an incomplete understanding of network dynamics and their regulatory properties. Kinetic modeling of plant metabolic networks can provide predictive information on network control and response to genetic perturbations, which allow estimation of flux at any concentration of intermediate or enzyme in the system. In this research, a kinetic model of the benzenoid network was developed to simulate whole network responses to different concentrations of supplied phenylalanine (Phe) in petunia flowers and capture flux redistributions caused by genetic manipulations. Kinetic parameters were obtained by network decomposition and non-linear least squares optimization of data from petunia flowers supplied with either 75 or 150 mm(2)H(5)-Phe. A single set of kinetic parameters simultaneously accommodated labeling and pool size data obtained for all endogenous and emitted volatiles at the two concentrations of supplied (2)H(5)-Phe. The generated kinetic model was validated using flowers from transgenic petunia plants in which benzyl CoA:benzyl alcohol/phenylethanol benzoyltransferase (BPBT) was down-regulated via RNAi. The determined in vivo kinetic parameters were used for metabolic control analysis, in which flux control coefficients were calculated for fluxes around the key branch point at Phe and revealed that phenylacetaldehyde synthase activity is the primary controlling factor for the phenylacetaldehyde branch of the benzenoid network. In contrast, control of flux through the beta-oxidative and non-beta-oxidative pathways is highly distributed.

Decreased expression of cytosolic pyruvate kinase in potato tubers leads to a decline in pyruvate resulting in an in vivo repression of the alternative oxidase

The aim of this work was to investigate the effect of decreased cytosolic pyruvate kinase (PKc) on potato (Solanum tuberosum) tuber metabolism. Transgenic potato plants with strongly reduced levels of PKc were generated by RNA interference gene silencing under the control of a tuber-specific promoter. Metabolite profiling showed that decreased PKc activity led to a decrease in the levels of pyruvate and some other organic acids involved in the tricarboxylic acid cycle. Flux analysis showed that this was accompanied by changes in carbon partitioning, with carbon flux being diverted from glycolysis toward starch synthesis. However, this metabolic shift was relatively small and hence did not result in enhanced starch levels in the tubers. Although total respiration rates and the ATP to ADP ratio were largely unchanged, transgenic tubers showed a strong decrease in the levels of alternative oxidase (AOX) protein and a corresponding decrease in the capacity of the alternative pathway of respiration. External feeding of pyruvate to tuber tissue or isolated mitochondria resulted in activation of the AOX pathway, both in the wild type and the PKc transgenic lines, providing direct evidence for the regulation of AOX by changes in pyruvate levels. Overall, these results provide evidence for a crucial role of PKc in the regulation of pyruvate levels as well as the level of the AOX in heterotrophic plant tissue, and furthermore reveal that these parameters are interlinked in vivo.

Metabolic control analysis: a tool for designing strategies to manipulate metabolic pathways

The traditional experimental approaches used for changing the flux or the concentration of a particular metabolite of a metabolic pathway have been mostly based on the inhibition or over-expression of the presumed rate-limiting step. However, the attempts to manipulate a metabolic pathway by following such approach have proved to be unsuccessful. Metabolic Control Analysis (MCA) establishes how to determine, quantitatively, the degree of control that a given enzyme exerts on flux and on the concentration of metabolites, thus substituting the intuitive, qualitative concept of rate limiting step. Moreover, MCA helps to understand (i) the underlying mechanisms by which a given enzyme exerts high or low control and (ii) why the control of the pathway is shared by several pathway enzymes and transporters. By applying MCA it is possible to identify the steps that should be modified to achieve a successful alteration of flux or metabolite concentration in pathways of biotechnological (e.g., large scale metabolite production) or clinical relevance (e.g., drug therapy). The different MCA experimental approaches developed for the determination of the flux-control distribution in several pathways are described. Full understanding of the pathway properties when is working under a variety of conditions can help to attain a successful manipulation of flux and metabolite concentration.

SBML-SAT: a systems biology markup language (SBML) based sensitivity analysis tool

Background

It has long been recognized that sensitivity analysis plays a key role in modeling and analyzing cellular and biochemical processes. Systems biology markup language (SBML) has become a well-known platform for coding and sharing mathematical models of such processes. However, current SBML compatible software tools are limited in their ability to perform global sensitivity analyses of these models.

Results

This work introduces a freely downloadable, software package, SBML-SAT, which implements algorithms for simulation, steady state analysis, robustness analysis and local and global sensitivity analysis for SBML models. This software tool extends current capabilities through its execution of global sensitivity analyses using multi-parametric sensitivity analysis, partial rank correlation coefficient, SOBOL's method, and weighted average of local sensitivity analyses in addition to its ability to handle systems with discontinuous events and intuitive graphical user interface.

Conclusion

SBML-SAT provides the community of systems biologists a new tool for the analysis of their SBML models of biochemical and cellular processes.

Decreased expression of plastidial adenylate kinase in potato tubers results in an enhanced rate of respiration and a stimulation of starch synthesis that is attributable to post-translational redox-activation of ADP-glucose pyrophosphorylase

Adenine nucleotides are of general importance for many aspects of cell function, but their role in the regulation of biosynthetic processes is still unclear. It was previously reported that decreased expression of plastidial adenylate kinase, catalysing the interconversion of ATP and AMP to ADP, leads to increased adenylate pools and starch content in transgenic potato tubers. However, the underlying mechanisms were not elucidated. Here, it is shown that decreased expression of plastidial adenylate kinase in growing tubers leads to increased rates of respiratory oxygen consumption and increased carbon fluxes into starch. Increased rates of starch synthesis were accompanied by post-translational redox-activation of ADP-glucose pyrophosphorylase (AGPase), catalysing the key regulatory step of starch synthesis in the plastid, while there were no substantial changes in metabolic intermediates or sugar levels. A similar increase in post-translational redox-activation of AGPase was found after supplying adenine to wild-type potato tuber discs to increase adenine nucleotide levels. Results provide first evidence for a link between redox-activation of AGPase and adenine nucleotide levels in plants.

Kinetic modeling can describe in vivo glycolysis in Entamoeba histolytica

Glycolysis in the human parasite Entamoeba histolytica is characterized by the absence of cooperative modulation and the prevalence of pyrophosphate-dependent (over ATP-dependent) enzymes. To determine the flux-control distribution of glycolysis and understand its underlying control mechanisms, a kinetic model of the pathway was constructed by using the software gepasi. The model was based on the kinetic parameters determined in the purified recombinant enzymes, and the enzyme activities, and steady-state fluxes and metabolite concentrations determined in amoebal trophozoites. The model predicted, with a high degree of accuracy, the flux and metabolite concentrations found in trophozoites, but only when the pyrophosphate concentration was held constant; at variable pyrophosphate, the model was not able to completely account for the ATP production/consumption balance, indicating the importance of the pyrophosphate homeostasis for amoebal glycolysis. Control analysis by the model revealed that hexokinase exerted the highest flux control (73%), as a result of its low cellular activity and strong AMP inhibition. 3-Phosphoglycerate mutase also exhibited significant flux control (65%) whereas the other pathway enzymes showed little or no control. The control of the ATP concentration was also mainly exerted by ATP consuming processes and 3-phosphoglycerate mutase and hexokinase (in the producing block). The model also indicated that, in order to diminish the amoebal glycolytic flux by 50%, it was required to decrease hexokinase or 3-phosphoglycerate mutase by 24% and 55%, respectively, or by 18% for both enzymes. By contrast, to attain the same reduction in flux by inhibiting the pyrophosphate-dependent enzymes pyrophosphate-phosphofructokinase and pyruvate phosphate dikinase, they should be decreased > 70%. On the basis of metabolic control analysis, steps whose inhibition would have stronger negative effects on the energy metabolism of this parasite were identified, thus becoming alternative targets for drug design.

Metabolic flux analysis in plants using dynamic labeling technique: application to tryptophan biosynthesis in cultured rice cells

The concept and methodology of using dynamic labeling for the MFA of plant metabolic pathways are described, based on a case study to develop a method for the MFA of the tryptophan biosynthetic pathway in cultured rice cells. Dynamic labeling traces the change in the labeling level of a metabolite in a metabolic pathway after the application of a stable isotope-labeled compound. In this study, [1-(13)C] l-serine was fed as a labeling precursor and the labeling level of Trp was determined by using the LC-MS/MS. The value of metabolic flux is determined by fitting a model describing the labeling dynamics of the pathway to the observed labeling data. The biosynthetic flux of Trp in rice suspension cultured cell was determined to be 6.0+/-1.1 nmol (gFWh)(-1). It is also demonstrated that an approximately sixfold increase in the biosynthetic flux of Trp in transgenic rice cells expressing the feedback-insensitive version of anthranilate synthase alpha-subunit gene (OASA1D) resulted in a 45-fold increase in the level of Trp. In this article, the basic workflow for the experiment is introduced and the details of the actual experimental procedures are explained. Future perspectives are also discussed by referring recent advances in the dynamic labeling approach.

Experimental and mathematical approaches to modeling plant metabolic networks

To support their sessile and autotrophic lifestyle higher plants have evolved elaborate networks of metabolic pathways. Dynamic changes in these metabolic networks are among the developmental forces underlying the functional differentiation of organs, tissues and specialized cell types. They are also important in the various interactions of a plant with its environment. Further complexity is added by the extensive compartmentation of the various interconnected metabolic pathways in plants. Thus, although being used widely for assessing the control of metabolic flux in microbes, mathematical modeling approaches that require steady-state approximations are of limited utility for understanding complex plant metabolic networks. However, considerable progress has been made when manageable metabolic subsystems were studied. In this article, we will explain in general terms and using simple examples the concepts underlying stoichiometric modeling (metabolic flux analysis and metabolic pathway analysis) and kinetic approaches to modeling (including metabolic control analysis as a special case). Selected studies demonstrating the prospects of these approaches, or combinations of them, for understanding the control of flux through particular plant pathways are discussed. We argue that iterative cycles of (dry) mathematical modeling and (wet) laboratory testing will become increasingly important for simulating the distribution of flux in plant metabolic networks and deriving rational experimental designs for metabolic engineering efforts.

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.

A bioinformatician's view of the metabolome

The study of a collection of metabolites as a whole (metabolome), as opposed to isolated small molecules, is a fast-growing field promising to take us one step further towards understanding cell biology, and relating the genetic capabilities of an organism to its observed phenotype. The new sciences of metabolomics and metabonomics can exploit a variety of existing experimental and computational methods, but they also require new technology that can deal with both the amount and the diversity of the data relating to the rich world of metabolites. More specifically, the collaboration between bioinformaticians and chemoinformaticians promises to advance our view of cognate molecules, by shedding light on their atomic structure and properties. Modelling of the interactions of metabolites with other entities in the cell, and eventually complete modelling of reaction pathways will be essential for analysis of the experimental data, and prediction of an organism's response to environmental challenges.

Determining and understanding the control of glycolysis in fast-growth tumor cells. Flux control by an over-expressed but strongly product-inhibited hexokinase

Control analysis of the glycolytic flux was carried out in two fast-growth tumor cell types of human and rodent origin (HeLa and AS-30D, respectively). Determination of the maximal velocity (V(max)) of the 10 glycolytic enzymes from hexokinase to lactate dehydrogenase revealed that hexokinase (153-306 times) and phosphofructokinase-1 (PFK-1) (22-56 times) had higher over-expression in rat AS-30D hepatoma cells than in normal freshly isolated rat hepatocytes. Moreover, the steady-state concentrations of the glycolytic metabolites, particularly those of the products of hexokinase and PFK-1, were increased compared with hepatocytes. In HeLa cells, V(max) values and metabolite concentrations for the 10 glycolytic enzyme were also significantly increased, but to a much lesser extent (6-9 times for both hexokinase and PFK-1). Elasticity-based analysis of the glycolytic flux in AS-30D cells showed that the block of enzymes producing Fru(1,6)P2 (i.e. glucose transporter, hexokinase, hexosephosphate isomerase, PFK-1, and the Glc6P branches) exerted most of the flux control (70-75%), whereas the consuming block (from aldolase to lactate dehydrogenase) exhibited the remaining control. The Glc6P-producing block (glucose transporter and hexokinase) also showed high flux control (70%), which indicated low flux control by PFK-1. Kinetic analysis of PFK-1 showed low sensitivity towards its allosteric inhibitors citrate and ATP, at physiological concentrations of the activator Fru(2,6)P2. On the other hand, hexokinase activity was strongly inhibited by high, but physiological, concentrations of Glc6P. Therefore, the enhanced glycolytic flux in fast-growth tumor cells was still controlled by an over-produced, but Glc6P-inhibited hexokinase.

Integrative analysis of metabolic networks: from peaks to flux models?

Recent developments in genomic and post-genomic technologies have led to the amassment of data describing genome sequences, transcript, protein and metabolite abundances, protein modifications, and protein-protein and protein-DNA interactions. Such technologies have vastly expanded the inventory of detectable molecular species and can be used to describe their interdependence, but they have yet to fulfill their promise in enhancing our knowledge of how flux through metabolic pathways is regulated. A convergence of traditional reductionistic and novel holistic experimental approaches could aid in elucidating flux control.

Cytosolic pyruvate kinase: subunit composition, activity, and amount in developing castor and soybean seeds, and biochemical characterization of the purified castor seed enzyme

Antibodies against Brassica napus cytosolic pyruvate kinase (PKc) (EC 2.7.1.40) were employed to examine PKc subunit composition and developmental profiles in castor and soybean seeds. A 56-kDa immunoreactive polypeptide was uniformly detected on immunoblots of clarified extracts from developing castor endosperm or soybean embryos. Maximal PKc activities occurred early in castor oil seed (COS) and soybean development (7.1 and 5.5 (micromol of pyruvate produced/min) g(-1) FW, respectively) and were up to 25-fold greater than those of fully mature seeds. Time-course studies revealed a close correlation between extractable PKc activity and the relative amount of the immunoreactive 56-kDa PKc polypeptide. PKc from developing COS was purified 1,874-fold to homogeneity and a final specific activity of 73.1 (micromol of pyruvate produced/min) mg(-1) protein. Gel filtration and SDS-PAGE indicated that this PKc exists as a 230-kDa homotetramer composed of 56-kDa subunits. The mass fingerprint of tryptic peptides of the 56-kDa COS PKc subunit best matched three putative PK(c)s from Arabidopsis thaliana. The purified enzyme was relatively heat-stable and displayed a broad pH optimum of 6.4. However, more efficient substrate utilization (in terms of Vmax /Km for phosphoenolpyruvate or ADP) was observed at pH 7.4. Glutamate was the most effective inhibitor, whereas aspartate functioned as an activator by partially relieving glutamate inhibition. Together with our previous studies, the results: (1) allow a model to be formulated regarding the coordinate allosteric control of PKc and phosphoenolpyruvate carboxylase by aspartate and glutamate in developing COS, and (2) provide further biochemical evidence that castor plant PKc exists as tissue-specific isozymes that exhibit substantial differences in their respective physical and regulatory properties.

Tuber-specific cytosolic expression of a bacterial phosphoglucomutase in potato (Solanum tuberosum L.) dramatically alters carbon partitioning

Constitutive antisense inhibition of the cytosolic isoform of phosphoglucomutase in the potato (Solanum tuberosum L.) results in restriction of photosynthesis, growth inhibition and modified tuber morphology, and a severe restriction of tuber starch synthesis. Here we describe the consequences of the tuber-specific expression of an Escherichia coli phosphoglucomutase in the cytosol. Analysis of [14C]glucose metabolism by tuber discs isolated from wild type and transformants revealed that the rates of sucrose and starch synthesis were unaltered but that the rate of glycolysis was depressed in the transgenics. The transformant tubers also contained dramatically reduced amino acid content and significantly higher levels of ADP, but were characterized by elevated levels of Krebs cycle intermediates and an unaltered rate of respiration. In addition to these metabolic consequences of the overexpression of the E. coli enzyme, we observed morphological changes in tubers, with the transformants having a smaller number of larger tubers which exhibited delayed rates of sprouting with respect to the wild type. These results are discussed with respect to current models of the regulation of central plant metabolism and tuber dormancy.

Purification and characterization of pyrophosphate- and ATP-dependent phosphofructokinases from banana fruit

Pyrophosphate-dependent phosphofructokinase (PFP; EC 2.7.1.90) and two isoforms of ATP-dependent phosphofructokinase (PFK I and PFK II; EC 2.7.1.11) from ripened banana ( Musa cavendishii L. cv. Cavendish) fruits were resolved via hydrophobic interaction fast protein liquid chromatography (FPLC), and further purified using anion-exchange and gel filtration FPLC. PFP was purified 1,158-fold to a final specific activity of 13.9 micromol fructose 1,6-bisphosphate produced (mg protein)(-1) x min(-1). Gel filtration FPLC and immunoblot analyses indicated that this PFP exists as a 490-kDa heterooctomer composed of equal amounts of 66- (alpha) and 60-kDa (beta) subunits. PFP displayed hyperbolic saturation kinetics for fructose 6-phosphate (Fru 6-P), PPi, fructose 1,6-bisphosphate, and Pi ( K(m) values = 32, 9.7, 25, and 410 microM, respectively) in the presence of saturating (5 microM) fructose 2,6-bisphosphate, which elicited a 24-fold enhancement of glycolytic PFP activity ( K(a)=8 nM). PFK I and PFK II were each purified about 350-fold to final specific activities of 5.5-6.0 micromol fructose 1,6-bisphosphate produced (mg protein)(-1) x min(-1). Analytical gel filtration yielded respective native molecular masses of 210 and 160 kDa for PFK I and PFK II. Several properties of PFK I and PFK II were consistent with their respective designation as plastid and cytosolic PFK isozymes. PFK I and PFK II exhibited: (i) pH optima of 8.0 and 7.3, respectively; (ii) hyperbolic saturation kinetics for ATP ( K(m)=34 and 21 microM, respectively); and (iii) sigmoidal saturation kinetics for Fru 6-P ( S0.5=540 and 90 microM, respectively). Allosteric effects of phospho enolpyruvate (PEP) and Pi on the activities of PFP, PFK I, and PFK II were characterized. Increasing concentrations of PEP or Pi progressively disrupted fructose 2,6-bisphosphate binding by PFP. PEP potently inhibited PFK I and to a lesser extent PFK II ( I50=2.3 and 900 microM, respectively), while Pi activated PFK I by reducing its sensitivity to PEP inhibition. Our results are consistent with: (i) the respiratory climacteric being regulated by fine (allosteric) control of pre-existing enzymes; and (ii) primary and secondary glycolytic flux control being exerted at the levels of PEP and Fru 6-P metabolism, respectively.

Metabolomics--the link between genotypes and phenotypes

Metabolites are the end products of cellular regulatory processes, and their levels can be regarded as the ultimate response of biological systems to genetic or environmental changes. In parallel to the terms 'transcriptome' and proteome', the set of metabolites synthesized by a biological system constitute its 'metabolome'. Yet, unlike other functional genomics approaches, the unbiased simultaneous identification and quantification of plant metabolomes has been largely neglected. Until recently, most analyses were restricted to profiling selected classes of compounds, or to fingerprinting metabolic changes without sufficient analytical resolution to determine metabolite levels and identities individually. As a prerequisite for metabolomic analysis, careful consideration of the methods employed for tissue extraction, sample preparation, data acquisition, and data mining must be taken. In this review, the differences among metabolite target analysis, metabolite profiling, and metabolic fingerprinting are clarified, and terms are defined. Current approaches are examined, and potential applications are summarized with a special emphasis on data mining and mathematical modelling of metabolism.

Mathematical modeling of plant metabolic pathways

The understanding of the control of metabolic flux in plants requires integrated mathematical formulations of gene and protein expression, enzyme kinetics, and developmental biology. Plants have a large number of metabolically active compartments, and non-steady-state conditions are frequently encountered. Consequently steady-state metabolic flux balance and isotopic flux balance modeling approaches have limited utility in probing plant metabolic systems. Transient isotopic flux analysis and kinetic modeling are powerful proven techniques for the quantification of metabolic fluxes in compartmentalized, dynamic metabolic systems. These tools are now widely used to address metabolic flux responses to environmental and genetic perturbations in plant metabolism. Continued developments in isotopic and kinetic modeling, quantifying metabolite exchange between compartments, and transcriptional and posttranscriptional regulatory mechanisms governing enzyme level and activity will enable simulation of large sections of plant metabolism under non-steady-state conditions. Metabolic control analysis will continue to make substantial contributions to the understanding of quantitative distribution of control of flux. From the synergy between mathematical models and experiments, creative methods for controlling the distribution of flux by genetic or environmental means will be discovered and rationally implemented.

Metabolomics--the link between genotypes and phenotypes

Metabolites are the end products of cellular regulatory processes, and their levels can be regarded as the ultimate response of biological systems to genetic or environmental changes. In parallel to the terms 'transcriptome' and proteome', the set of metabolites synthesized by a biological system constitute its 'metabolome'. Yet, unlike other functional genomics approaches, the unbiased simultaneous identification and quantification of plant metabolomes has been largely neglected. Until recently, most analyses were restricted to profiling selected classes of compounds, or to fingerprinting metabolic changes without sufficient analytical resolution to determine metabolite levels and identities individually. As a prerequisite for metabolomic analysis, careful consideration of the methods employed for tissue extraction, sample preparation, data acquisition, and data mining must be taken. In this review, the differences among metabolite target analysis, metabolite profiling, and metabolic fingerprinting are clarified, and terms are defined. Current approaches are examined, and potential applications are summarized with a special emphasis on data mining and mathematical modelling of metabolism.

Control of the threonine-synthesis pathway in Escherichia coli: a theoretical and experimental approach

A computer simulation of the threonine-synthesis pathway in Escherichia coli Tir-8 has been developed based on our previous measurements of the kinetics of the pathway enzymes under near-physiological conditions. The model successfully simulates the main features of the time courses of threonine synthesis previously observed in a cell-free extract without alteration of the experimentally determined parameters, although improved quantitative fits can be obtained with small parameter adjustments. At the concentrations of enzymes, precursors and products present in cells, the model predicts a threonine-synthesis flux close to that required to support cell growth. Furthermore, the first two enzymes operate close to equilibrium, providing an example of a near-equilibrium feedback-inhibited enzyme. The predicted flux control coefficients of the pathway enzymes under physiological conditions show that the control of flux is shared between the first three enzymes: aspartate kinase, aspartate semialdehyde dehydrogenase and homoserine dehydrogenase, with no single activity dominating the control. The response of the model to the external metabolites shows that the sharing of control between the three enzymes holds across a wide range of conditions, but that the pathway flux is sensitive to the aspartate concentration. When the model was embedded in a larger model to simulate the variable demands for threonine at different growth rates, it showed the accumulation of free threonine that is typical of the Tir-8 strain at low growth rates. At low growth rates, the control of threonine flux remains largely with the pathway enzymes. As an example of the predictive power of the model, we studied the consequences of over-expressing different enzymes in the pathway.

Effects of overexpression of the liver subunit of 6-phosphofructo-1-kinase on the metabolism of a cultured mammalian cell line

Overexpression of the liver subunit of 6-phosphofructo-1-kinase in Chinese hamster ovary K1 cells was shown to increase the steady-state level of the enzyme's product, fructose 1, 6-bisphosphate, and to produce a small but significant decrease in the concentration of fructose 2,6-bisphosphate, which is an allosteric activator of the enzyme. However, overexpression of the enzyme had no effect on glycolytic flux under a variety of different substrate conditions. This latter observation is consistent with similar studies in fungi and in potato tubers which indicate that 6-phosphofructo-1-kinase has very little control over flux in glycolysis.

The role of the regulatory protein of glucokinase in the glucose sensory mechanism of the hepatocyte

Glucokinase has a very high flux control coefficient (greater than unity) on glycogen synthesis from glucose in hepatocytes (Agius et al., J. Biol. Chem. 271, 30479-30486, 1996). Hepatic glucokinase is inhibited by a 68-kDa glucokinase regulatory protein (GKRP) that is expressed in molar excess. To establish the relative control exerted by glucokinase and GKRP, we applied metabolic control analysis to determine the flux control coefficient of GKRP on glucose metabolism in hepatocytes. Adenovirus-mediated overexpression of GKRP (by up to 2-fold above endogenous levels) increased glucokinase binding and inhibited glucose phosphorylation, glycolysis, and glycogen synthesis over a wide range of concentrations of glucose and sorbitol. It decreased the affinity of glucokinase translocation for glucose and increased the control coefficient of glucokinase on glycogen synthesis. GKRP had a negative control coefficient of glycogen synthesis that is slightly greater than unity (-1.2) and a control coefficient on glycolysis of -0.5. The control coefficient of GKRP on glycogen synthesis decreased with increasing glucokinase overexpression (4-fold) at elevated glucose concentration (35 mM), which favors dissociation of glucokinase from GKRP, but not at 7.5 mM glucose. Under the latter conditions, glucokinase and GKRP have large and inverse control coefficients on glycogen synthesis, suggesting that a large component of the positive control coefficient of glucokinase is counterbalanced by the negative coefficient of GKRP. It is concluded that glucokinase and GKRP exert reciprocal control; therefore, mutations in GKRP affecting the expression or function of the protein may impact the phenotype even in the heterozygote state, similar to glucokinase mutations in maturity onset diabetes of the young type 2. Our results show that the mechanism comprising glucokinase and GKRP confers a markedly extended responsiveness and sensitivity to changes in glucose concentration on the hepatocyte.

Regulation in metabolic systems under homeostatic flux control

The general properties of metabolic systems under homeostatic flux control are analyzed. It is shown that the main characteristic point for an enzyme in such a system is a sharp transition from limitation outside the system to limitation by some enzyme inside the system. A method for the quantitative treatment of the experimental dependence of metabolic flux on enzyme content is presented. The conception of "nonlimiting," "near-limiting," and "limiting" enzymes is developed for these systems. It is pointed out that reactions close to a thermodynamic equilibrium under normal conditions can considerably limit the homeostatic fluxes. The rules for regulation of fluxes in such systems are illustrated by the data obtained for transgenic plants with reduced activities of some Calvin-cycle enzymes and further examples. A comparison is made between the developed quantitative description of metabolic fluxes under homeostatic flux control and the methods of metabolic control analysis.

Modelling photosynthesis and its control

The dynamic and steady-state behaviour of a computer simulation of the Calvin cycle reactions of the chloroplast, including starch synthesis and degradation, and triose phosphate export have been investigated. A major difference compared with previous models is that none of the reversible reactions are assumed to be at equilibrium. The model can exhibit alternate steady states of low or high carbon assimilation flux, with hysteresis in the transitions between the steady states induced by environmental factors such as phosphate and light intensity. The enzymes which have the greatest influence on the flux have been investigated by calculation of their flux control coefficients. Different patterns of control are exhibited over the assimilation flux, the flux to starch and the flux to cytosolic triose phosphate. The assimilation flux is mostly sensitive to sedoheptulose bisphosphatase and Rubisco, with the exact distribution depending on their relative activities. Other enzymes, particularly the triose phosphate translocator, become more influential when other fluxes are considered. These results are shown to be broadly consistent with observations on transgenic plants.

Top-down control analysis of ATP turnover, glycolysis and oxidative phosphorylation in rat hepatocytes

Control analysis was used to analyse the internal control of rat hepatocyte metabolism. The reactions of the cell were grouped into nine metabolic blocks linked by five key intermediates. The blocks were glycogen breakdown, glucose release, glycolysis, lactate production, NADH oxidation, pyruvate oxidation, mitochondrial proton leak, mitochondrial phosphorylation and ATP consumption. The linking intermediates were intracellular glucose-6-phosphate, pyruvate and ATP levels, cytoplasmic NADH/NAD ratio and mitochondrial membrane potential. The steady-state fluxes through the blocks and the levels of the intermediates were measured in the absence and presence of specific effectors of hepatocyte metabolism. Application of the multiple modulation approach gave the kinetic responses of each block to each intermediate (the elasticities). These were then used to calculate all of the control coefficients, which describe the degree of control each block had over the level of each intermediate, and over the rate of each process. Within this full description of control, many different interactions could be identified. One key finding was that the processes that consumed ATP had only 35% of the control over the rate of ATP consumption. Instead, the reactions that produced ATP exerted the most control over ATP consumption rate; particularly important were mitochondrial phosphorylation (30% of control) and glycolysis (19%). The rate of glycolysis was positively controlled by the glycolytic enzymes themselves (66% of control) and by ATP consumption (47%). Mitochondrial production of ATP, including oxidative, proton leak and phosphorylation processes, had negative control over glycolysis (-26%; the Pasteur effect). In contrast, glycolysis had little control over the rate of ATP production by the mitochondria (-10%; the Crabtree effect). Control over the flux through the mitochondrial phosphorylation block was shared between pyruvate oxidation (23%), ATP consumption (28%) and the mitochondrial phosphorylation block itself (64%).

Combined expression of glucokinase and invertase in potato tubers leads to a dramatic reduction in starch accumulation and a stimulation of glycolysis

The original aim of this work was to increase starch accumulation in potato tubers by enhancing their capacity to metabolise sucrose.We previously reported that specific expression of a yeast invertase in the cytosol of tubers led to a 95% reduction in sucrose content, but that this was accompanied by a larger accumulation of glucose and a reduction in starch. In the present paper we introduced a bacterial glucokinase from Zymomonas mobilis into an invertase-expressing transgenic line, with the intention of bringing the glucose into metabolism. Transgenic lines were obtained with up to threefold more glucokinase activity than in the parent invertase line and which did not accumulate glucose. Unexpectedly, there was a further dramatic reduction in starch content, down to 35% of wild-type levels. Biochemical analysis of growing tuber tissue revealed large increases in the metabolic intermediates of glycolysis, organic acids and amino acids,two- to threefold increases in the maximum catalytic activities of key enzymes in the respiratory pathways, and three- to fivefold increases in carbon dioxide production.These changes occur in the lines expressing invertase,and are accentuated following introduction of the second transgene, glucokinase. We conclude that the expression of invertase in potato tubers leads to an increased flux through the glycolytic pathway at the expense of starch synthesis and that heterologous overexpression of glucokinase enhances this change in partitioning.

Finite change analysis of glycolytic intermediates in tuber tissue of lines of transgenic potato (Solanum tuberosum) overexpressing phosphofructokinase

Genetically engineered organisms overexpressing phosphofructokinase (PFK), a supposed 'regulatory' step of glycolysis, often show little or no measurable change in glycolytic or respiratory flux, although the concentrations of glycolytic intermediates may change. We have used the finite change theory of Metabolic Control Analysis (MCA) to analyse the concentrations of glycolytic metabolites in aged disks of tuber tissue from four lines of transgenic potatoes expressing different amounts of PFK that, under aerobic conditions, showed statistically indistinguishable rates of respiration. The constancy of the metabolites' concentration deviation indices for different increases in PFK expression indicated that the metabolite changes from a graded series, excluding the possibility of anomalous behaviour that might be observed in a single transgenic line. Consequently we were able to use the finite change method to validate the results of an MCA model of tuber glycolysis [Thomas, Mooney, Burrell and Fell (1997) Biochem. J. 322, 119-127]. Furthermore the metabolite changes with PFK activity are evidence that near-equilibrium steps do not transmit increased substrate concentrations down the pathway without attenuation. Our results support the view that flux increase by activation of a single enzyme early in the pathway will, contrary to expectations, be of limited effectiveness in achieving flux increases.