A kinetic model of the branch-point between the methionine and threonine biosynthesis pathways in Arabidopsis thaliana

The Arabidopsis electron-transfer flavoprotein:ubiquinone oxidoreductase is required during normal seed development and germination

The importance of the alternative donation of electrons to the ubiquinol pool via the electron-transfer flavoprotein/electron-transfer flavoprotein:ubiquinone oxidoreductase (ETF/ETFQO) complex has been demonstrated. However, the functional significance of this pathway during seed development and germination remains to be elucidated. To assess the function of this pathway, we performed a detailed metabolic and transcriptomic analysis of Arabidopsis mutants to test the molecular consequences of a dysfunctional ETF/ETFQO pathway. We demonstrate that the disruption of this pathway compromises seed germination in the absence of an external carbon source and also impacts seed size and yield. Total protein and storage protein content is reduced in dry seeds, whilst sucrose levels remain invariant. Seeds of ETFQO and related mutants were also characterized by an altered fatty acid composition. During seed development, lower levels of fatty acids and proteins accumulated in the etfqo-1 mutant as well as in mutants in the alternative electron donors isovaleryl-CoA dehydrogenase (ivdh-1) and d-2-hydroxyglutarate dehydrogenase (d2hgdh1-2). Furthermore, the content of several amino acids was increased in etfqo-1 mutants during seed development, indicating that these mutants are not using such amino acids as alternative energy source for respiration. Transcriptome analysis revealed alterations in the expression levels of several genes involved in energy and hormonal metabolism. Our findings demonstrated that the alternative pathway of respiration mediated by the ETF/ETFQO complex affects seed germination and development by directly adjusting carbon storage during seed filling. These results indicate a role for the pathway in the normal plant life cycle to complement its previously defined roles in the response to abiotic stress.

The level of threonine in tobacco seeds is limited by substrate availability, while the level of methionine is limited also by the activity of cystathionine γ-synthase

Methionine and threonine are two essential amino acids whose low levels limit the nutritional quality of seeds. The current objective was to define factors that regulate and might increase their levels in seeds. Feeding experiments carried out on receptacles of developing tobacco (Nicotiana tabacum) capsules showed that 1 mM of S-methylmethionine increased the level of methionine to contents similar to 2.5 mM of homoserine, an intermediate metabolite of the aspartate family of amino acids. The latter also increased the level of threonine. Based on these findings, we generated tobacco seeds that expressed a combination of bacterial feedback-insensitive aspartate kinase (bAK), which was previously reported to have a high level of threonine/methionine, and feedback-insensitive cystathionine γ-synthase (CGS), the regulatory enzyme of the methionine biosynthesis pathway. Plants expressing this latter gene previously showed having higher levels of methionine. The results of total amino acids analysis showed that the level of threonine was highest in the bAK line, which has moderate levels of methionine and lysine, while the highest level of methionine was found in seeds expressing both heterologous genes. The results suggest that the level of threonine in tobacco seeds is limited by the substrate, while that of methionine is limited also by the activity of CGS.

Distribution of control in the sulfur assimilation in Arabidopsis thaliana depends on environmental conditions

Sulfur assimilation is central to the survival of plants and has been studied under different environmental conditions. Multiple studies have been published trying to determine rate-limiting or controlling steps in this pathway. However, the picture remains inconclusive with at least two different enzymes proposed to represent such rate-limiting steps. Here, we used computational modeling to gain an integrative understanding of the distribution of control in the sulfur assimilation pathway of Arabidopsis thaliana. For this purpose, we set up a new ordinary differential equation (ODE)-based, kinetic model of sulfur assimilation encompassing all biochemical reactions directly involved in this process. We fitted the model to published experimental data and produced a model ensemble to deal with parameter uncertainties. The ensemble was validated against additional published experimental data. We used the model ensemble to subsequently analyse the control pattern and robustly identified a set of processes that share the control in this pathway under standard conditions. Interestingly, the pattern of control is dynamic and not static, that is it changes with changing environmental conditions. Therefore, while adenosine-5'-phosphosulfate reductase (APR) and sulfite reductase (SiR) share control under standard laboratory conditions, APR takes over an even more dominant role under sulfur starvation conditions.

Using metabolome data for mathematical modeling of plant metabolic systems

Plant metabolism is characterized by a wide diversity of metabolites, with systems far more complicated than those of microorganisms. Mathematical modeling is useful for understanding dynamic behaviors of plant metabolic systems for metabolic engineering. Time-series metabolome data has great potential for estimating kinetic model parameters to construct a genome-wide metabolic network model. However, data obtained by current metabolomics techniques does not meet the requirement for constructing accurate models. In this article, we highlight novel strategies and algorithms to handle the underlying difficulties and construct dynamic in vivo models for large-scale plant metabolic systems. The coarse but efficient modeling enables the prediction of unknown mechanisms regulating plant metabolism.

Inhibition of homoserine dehydrogenase by formation of a cysteine-NAD covalent complex

Homoserine dehydrogenase (EC 1.1.1.3, HSD) is an important regulatory enzyme in the aspartate pathway, which mediates synthesis of methionine, threonine and isoleucine from aspartate. Here, HSD from the hyperthermophilic archaeon Sulfolobus tokodaii (StHSD) was found to be inhibited by cysteine, which acted as a competitive inhibitor of homoserine with a Ki of 11 μM and uncompetitive an inhibitor of NAD and NADP with Ki's of 0.55 and 1.2 mM, respectively. Initial velocity and product (NADH) inhibition analyses of homoserine oxidation indicated that StHSD first binds NAD and then homoserine through a sequentially ordered mechanism. This suggests that feedback inhibition of StHSD by cysteine occurs through the formation of an enzyme-NAD-cysteine complex. Structural analysis of StHSD complexed with cysteine and NAD revealed that cysteine situates within the homoserine binding site. The distance between the sulfur atom of cysteine and the C4 atom of the nicotinamide ring was approximately 1.9 Å, close enough to form a covalent bond. The UV absorption-difference spectrum of StHSD with and without cysteine in the presence of NAD, exhibited a peak at 325 nm, which also suggests formation of a covalent bond between cysteine and the nicotinamide ring.

Computational fluid dynamic (CFD) modeling of simultaneous extraction and fermentation process in a single sugar beet cossette

For simulations of flow and microbial conversion reactions, related to modeling of simultaneous extraction and fermentation process in a single sugar beet cossette a software package OpenFOAM was used. The mass transfer of the components (sucrose, glucose, fructose and ethanol) in the studied system was controlled by the convection and diffusion processes. Microbial conversion rates and yield coefficients were experimentally determined and/or estimated by mathematical simulation. Dimensions of the model sugar beet cossette (SBC) were: average length of cosettes 40.10 mm, average thickness 3.32 mm and average width 3.5 mm, and represented in the model as a square-shape cross-section mathematical simulation. Dimensions of the model sugar beet cossette (SBC) were: average length of cosettes 40.10 mm, average thickness 3.32 mm and average width 3.5 mm, and represented in the model as a square-shape cross-section used to study the mass transfer and microbial conversion rates on the scale of single sugar beet cossette in the short time scales (up to 25 s). This model can be used for simulation of extractant flow around single sugar beet cossette as well as for description of simultaneous extraction and fermentation process in the studied system.

Improving Nutritional Quality of Plant Proteins Through Genetic Engineering

Humans and animals are unable to synthesize essential amino acids such as branch chain amino acids methionine (Met), lysine (Lys) and tryptophan (Trp). Therefore, these amino acids need to be supplied through the diets. Several essential amino acids are deficient or completely lacking among crops used for human food and animal feed. For example, soybean is deficient in Met; Lys and Trp are lacking in maize. In this mini review, we will first summarize the roles of essential amino acids in animal nutrition. Next, we will address the question: “What are the amino acids deficient in various plants and their biosynthesis pathways?” And: “What approaches are being used to improve the availability of essential amino acids in plants?” The potential targets for metabolic engineering will also be discussed, including what has already been done and what remains to be tested.

Methionine and S-methylmethionine exhibit temporal and spatial accumulation patterns during the Arabidopsis life cycle

Methionine is a nutritionally essential sulfur-containing amino acid found at low levels in plant tissues. Yet, the factors that regulate its synthesis and accumulation in seeds are not fully known. Recent genetic studies demonstrate that Arabidopsis seeds are able to synthesize methionine de novo through the aspartate family pathway similarly to vegetative tissues; however, additional biochemical studies suggest that the S-methylmethionine (SMM) cycle also plays a major role in methionine synthesis in seeds. To better understand the contribution of these two pathways to methionine synthesis, we have sampled various vegetative and reproductive tissues during the Arabidopsis life cycle and determined the contents of soluble and protein-incorporated methionine, SMM, as well as the expression levels of the key genes involved in these two pathways. Our results strengthen the hypothesis that SMM that is produced in the rosette leaves from methionine contributes to methionine accumulation in seeds. However, the SMM cycle may have additional functions in plant tissues since its key genes were expressed in all of the examined tissues, although at different rates. The accumulation patterns of soluble and protein-incorporated methionine during the Arabidopsis life cycle were found to be similar to most of the other amino acids, especially to those belonging to the branched-chain and aromatic amino acids that are produced in chloroplasts together with methionine. This indicates that similar factors regulate the levels of amino acids during development.

Predictive sulfur metabolism - a field in flux

The key role of sulfur metabolites in response to biotic and abiotic stress in plants, as well as their importance in diet and health has led to a significant interest and effort in trying to understand and manipulate the production of relevant compounds. Metabolic engineering utilizes a set of theoretical tools to help rationally design modifications that enhance the production of a desired metabolite. Such approaches have proven their value in bacterial systems, however, the paucity of success stories to date in plants, suggests that challenges remain. Here, we review the most commonly used methods for understanding metabolic flux, focusing on the sulfur assimilatory pathway. We highlight known issues with both experimental and theoretical approaches, as well as presenting recent methods for integrating different modeling strategies, and progress toward an understanding of flux at the whole plant level.

Analytical kinetic modeling: a practical procedure

This chapter describes a practical procedure to dissect metabolic systems, simplify them, and use or derive enzyme rate equations in order to build a mathematical model of a metabolic system and run simulations. We first deal with a simple example, modeling a single enzyme that follows Michaelis-Menten kinetics and operates in the middle of an unbranched metabolic pathway. Next we describe the rules that can be followed to isolate sub-systems from their environment to simulate their behavior. Finally we use examples to show how to derive suitable rate equations, simpler than those needed for mechanistic studies, though adequate to describe the behavior over the physiological range of conditions.Many of the general characteristics of kinetic models will be obvious to readers familiar with the theory of metabolic control analysis (Cornish-Bowden, Fundamentals of Enzyme Kinetics, Wiley-Blackwell, Weinheim, 327-380, 2012), but here we shall not assume such knowledge, as the chapter is directed toward practical application rather than theory.

Application of metabolic flux analysis to plants

This volume compiles a series of chapters that cover the major aspects of plant metabolic flux analysis, such as but not limited to labeling of plant material, acquisition of labeling data, mathematical modeling of metabolic network at the cell, tissue, and plant level. A short revue, including methodological points and applications of flux analysis to plants, is presented in this introductory chapter.

Fortifying plants with the essential amino acids lysine and methionine to improve nutritional quality

Humans, as well as farm animals, cannot synthesize a number of essential amino acids, which are critical for their survival. Hence, these organisms must obtain these essential amino acids from their diets. Cereal and legume crops, which represent the major food and feed sources for humans and livestock worldwide, possess limiting levels of some of these essential amino acids, particularly Lys and Met. Extensive efforts were made to fortify crop plants with these essential amino acids using traditional breeding and mutagenesis. However, aside from some results obtained with maize, none of these approaches was successful. Therefore, additional efforts using genetic engineering approaches concentrated on increasing the synthesis and reducing the catabolism of these essential amino acids and also on the expression of recombinant proteins enriched in them. In the present review, we discuss the basic biological aspects associated with the synthesis and accumulation of these amino acids in plants and also describe recent developments associated with the fortification of crop plants with essential amino acids by genetic engineering approaches.

Transcripts of sulphur metabolic genes are co-ordinately regulated in developing seeds of common bean lacking phaseolin and major lectins

The lack of phaseolin and phytohaemagglutinin in common bean (dry bean, Phaseolus vulgaris) is associated with an increase in total cysteine and methionine concentrations by 70% and 10%, respectively, mainly at the expense of an abundant non-protein amino acid, S-methyl-cysteine. Transcripts were profiled between two genetically related lines differing for this trait at four stages of seed development using a high density microarray designed for common bean. Transcripts of multiple sulphur-rich proteins were elevated, several previously identified by proteomics, including legumin, basic 7S globulin, albumin-2, defensin, albumin-1, the Bowman-Birk type proteinase inhibitor, the double-headed trypsin inhibitor, and the Kunitz trypsin inhibitor. A co-ordinated regulation of transcripts coding for sulphate transporters, sulphate assimilatory enzymes, serine acetyltransferases, cystathionine β-lyase, homocysteine S-methyltransferase and methionine gamma-lyase was associated with changes in cysteine and methionine concentrations. Differential gene expression of sulphur-rich proteins preceded that of sulphur metabolic enzymes, suggesting a regulation by demand from the protein sink. Up-regulation of SERAT1;1 and -1;2 expression revealed an activation of cytosolic O-acetylserine biosynthesis. Down-regulation of SERAT2;1 suggested that cysteine and S-methyl-cysteine biosynthesis may be spatially separated in different subcellular compartments. Analysis of free amino acid profiles indicated that enhanced cysteine biosynthesis was correlated with a depletion of O-acetylserine. These results contribute to our understanding of the regulation of sulphur metabolism in developing seed in response to a change in the composition of endogenous proteins.

Characterization of chloroplastic fructose 1,6-bisphosphate aldolases as lysine-methylated proteins in plants

In pea (Pisum sativum), the protein-lysine methyltransferase (PsLSMT) catalyzes the trimethylation of Lys-14 in the large subunit (LS) of ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco), the enzyme catalyzing the CO(2) fixation step during photosynthesis. Homologs of PsLSMT, herein referred to as LSMT-like enzymes, are found in all plant genomes, but methylation of LS Rubisco is not universal in the plant kingdom, suggesting a species-specific protein substrate specificity of the methyltransferase. In this study, we report the biochemical characterization of the LSMT-like enzyme from Arabidopsis thaliana (AtLSMT-L), with a focus on its substrate specificity. We show that, in Arabidopsis, LS Rubisco is not naturally methylated and that the physiological substrates of AtLSMT-L are chloroplastic fructose 1,6-bisphosphate aldolase isoforms. These enzymes, which are involved in the assimilation of CO(2) through the Calvin cycle and in chloroplastic glycolysis, are trimethylated at a conserved lysyl residue located close to the C terminus. Both AtLSMT-L and PsLSMT are able to methylate aldolases with similar kinetic parameters and product specificity. Thus, the divergent substrate specificity of LSMT-like enzymes from pea and Arabidopsis concerns only Rubisco. AtLSMT-L is able to interact with unmethylated Rubisco, but the complex is catalytically unproductive. Trimethylation does not modify the kinetic properties and tetrameric organization of aldolases in vitro. The identification of aldolases as methyl proteins in Arabidopsis and other species like pea suggests a role of protein lysine methylation in carbon metabolism in chloroplasts.

Development of a continuous assay and steady-state characterization of Escherichia coli threonine synthase

Threonine synthase (TS) catalyzes the hydrolysis of O-phospho-L-homoserine (OPHS) to produce L-threonine (L-Thr) and inorganic phosphate. Here, we report a simplified purification protocol for the OPHS substrate and a continuous, coupled-coupled, spectrophotometric TS assay. The sequential actions of threonine deaminase (TD) and hydroxyisocaproate dehydrogenase (HO-HxoDH) convert the L-Thr product of TS to α-ketobutyrate (α-KB) and then to 2-hydroxybutyrate, respectively, and are monitored as the decrease in absorbance at 340 nm resulting from the concomitant oxidation of β-nicotinamide adenine dinucleotide (NADH) to NAD(+) by HO-HxoDH. The effect of pH on the activities of Escherichia coli TD and Lactobacillus delbrueckii HO-HxoDH was determined to establish this continuous assay as suitable for steady-state characterization and to facilitate the optimization of coupling enzyme concentrations under different assay conditions to enable studies of TS across phyla. To validate this assay, TS from E. coli was characterized. The kinetic parameters (k(cat)=4s(-1) and K(m)=0.34 mM) and the pH optimum of 8.7, determined using the continuous assay, are consistent with values reported for this enzyme based on the discontinuous malachite green assay. The k(cat)/K(m)(OPHS) versus pH profile of E. coli TS is bell-shaped, and the apparent pK(a) values for the acidic and basic limbs are 7.1 and 10.4, respectively.

Metabolically engineered soybean seed with enhanced threonine levels: biochemical characterization and seed-specific expression of lysine-insensitive variants of aspartate kinases from the enteric bacterium Xenorhabdus bovienii

Threonine (Thr) is one of a few limiting essential amino acids (EAAs) in the animal feed industry, and its level in feed rations can impact production of important meat sources, such as swine and poultry. Threonine as well as EAAs lysine (Lys) and methionine (Met) are all synthesized via the aspartate family pathway. Here, we report a successful strategy to produce high free threonine soybean seed via identification of a feedback-resistant aspartate kinase (AK) enzyme that can be over-expressed in developing soybean seed. Towards this goal, we have purified and biochemically characterized AK from the enteric bacterium Xenorhabdus bovienii (Xb). Site-directed mutagenesis of XbAK identified two key regulatory residues Glu-257 and Thr-359 involved in lysine inhibition. Three feedback-resistant alleles, XbAK_T359I, XbAK_E257K and XbAK_E257K/T359I, have been generated. This study is the first to kinetically characterize the XbAK enzyme and provide biochemical and transgenic evidence that Glu-257 near the catalytic site is a critical residue for the allosteric regulation of AK. Furthermore, seed-specific expression of the feedback-resistant XbAK_T359I or XbAK_E257K allele results in increases of free Thr levels of up to 100-fold in R(1) soybean seed when compared to wild-type. Expression of feedback-sensitive wild-type AK did not substantially impact seed Thr content. In addition to high Thr, transgenic seed also showed substantial increases in other major free amino acid (FAA) levels, resulting in an up to 3.5-fold increase in the total FAA content. The transgenic seed was normal in appearance and germinated well under greenhouse conditions.

Network analysis of enzyme activities and metabolite levels and their relationship to biomass in a large panel of Arabidopsis accessions

Natural genetic diversity provides a powerful resource to investigate how networks respond to multiple simultaneous changes. In this work, we profile maximum catalytic activities of 37 enzymes from central metabolism and generate a matrix to investigate species-wide connectivity between metabolites, enzymes, and biomass. Most enzyme activities change in a highly coordinated manner, especially those in the Calvin-Benson cycle. Metabolites show coordinated changes in defined sectors of metabolism. Little connectivity was observed between maximum enzyme activities and metabolites, even after applying multivariate analysis methods. Measurements of posttranscriptional regulation will be required to relate these two functional levels. Individual enzyme activities correlate only weakly with biomass. However, when they are used to estimate protein abundances, and the latter are summed and expressed as a fraction of total protein, a significant positive correlation to biomass is observed. The correlation is additive to that obtained between starch and biomass. Thus, biomass is predicted by two independent integrative metabolic biomarkers: preferential investment in photosynthetic machinery and optimization of carbon use.

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.

SBRML: a markup language for associating systems biology data with models

MOTIVATION

Research in systems biology is carried out through a combination of experiments and models. Several data standards have been adopted for representing models (Systems Biology Markup Language) and various types of relevant experimental data (such as FuGE and those of the Proteomics Standards Initiative). However, until now, there has been no standard way to associate a model and its entities to the corresponding datasets, or vice versa. Such a standard would provide a means to represent computational simulation results as well as to frame experimental data in the context of a particular model. Target applications include model-driven data analysis, parameter estimation, and sharing and archiving model simulations.

RESULTS

We propose the Systems Biology Results Markup Language (SBRML), an XML-based language that associates a model with several datasets. Each dataset is represented as a series of values associated with model variables, and their corresponding parameter values. SBRML provides a flexible way of indexing the results to model parameter values, which supports both spreadsheet-like data and multidimensional data cubes. We present and discuss several examples of SBRML usage in applications such as enzyme kinetics, microarray gene expression and various types of simulation results.

AVAILABILITY AND IMPLEMENTATION

The XML Schema file for SBRML is available at http://www.comp-sys-bio.org/SBRML under the Academic Free License (AFL) v3.0.

Recent progress in deciphering the biosynthesis of aspartate-derived amino acids in plants

Plants are either directly or indirectly the source of most of the essential amino acids in animal diets. Four of these essential amino acids-methionine, threonine, isoleucine, and lysine-are all produced from aspartate via a well studied biosynthesis pathway. Given the nutritional interest in essential amino acids, the aspartate-derived amino acid pathway has been the subject of extensive research. Additionally, several pathway enzymes serve as targets for economically important herbicides, and some of the downstream products are biosynthetic precursors for other essential plant metabolites such as ethylene and S-adenosylmethionine. Recent and ongoing research on the aspartate-derived family of amino acids has identified new enzyme activities, regulatory mechanisms, and in vivo metabolic functions. Together, these discoveries will open up new possibilities for plant metabolic engineering.

Understanding the regulation of aspartate metabolism using a model based on measured kinetic parameters

The aspartate-derived amino-acid pathway from plants is well suited for analysing the function of the allosteric network of interactions in branched pathways. For this purpose, a detailed kinetic model of the system in the plant model Arabidopsis was constructed on the basis of in vitro kinetic measurements. The data, assembled into a mathematical model, reproduce in vivo measurements and also provide non-intuitive predictions. A crucial result is the identification of allosteric interactions whose function is not to couple demand and supply but to maintain a high independence between fluxes in competing pathways. In addition, the model shows that enzyme isoforms are not functionally redundant, because they contribute unequally to the flux and its regulation. Another result is the identification of the threonine concentration as the most sensitive variable in the system, suggesting a regulatory role for threonine at a higher level of integration.

Computational modeling of biochemical networks using COPASI

Computational modeling and simulation of biochemical networks is at the core of systems biology and this includes many types of analyses that can aid understanding of how these systems work. COPASI is a generic software package for modeling and simulation of biochemical networks which provides many of these analyses in convenient ways that do not require the user to program or to have deep knowledge of the numerical algorithms. Here we provide a description of how these modeling techniques can be applied to biochemical models using COPASI. The focus is both on practical aspects of software usage as well as on the utility of these analyses in aiding biological understanding. Practical examples are described for steady-state and time-course simulations, stoichiometric analyses, parameter scanning, sensitivity analysis (including metabolic control analysis), global optimization, parameter estimation, and stochastic simulation. The examples used are all published models that are available in the BioModels database in SBML format.

Aspartate-Derived Amino Acid Biosynthesis in Arabidopsis thaliana

The aspartate-derived amino acid pathway in plants leads to the biosynthesis of lysine, methionine, threonine, and isoleucine. These four amino acids are essential in the diets of humans and other animals, but are present in growth-limiting quantities in some of the world's major food crops. Genetic and biochemical approaches have been used for the functional analysis of almost all Arabidopsis thaliana enzymes involved in aspartate-derived amino acid biosynthesis. The branch-point enzymes aspartate kinase, dihydrodipicolinate synthase, homoserine dehydrogenase, cystathionine gamma synthase, threonine synthase, and threonine deaminase contain well-studied sites for allosteric regulation by pathway products and other plant metabolites. In contrast, relatively little is known about the transcriptional regulation of amino acid biosynthesis and the mechanisms that are used to balance aspartate-derived amino acid biosynthesis with other plant metabolic needs. The aspartate-derived amino acid pathway provides excellent examples of basic research conducted with A. thaliana that has been used to improve the nutritional quality of crop plants, in particular to increase the accumulation of lysine in maize and methionine in potatoes.

Overexpression of mutated forms of aspartate kinase and cystathionine gamma-synthase in tobacco leaves resulted in the high accumulation of methionine and threonine

Methionine and threonine are two essential amino acids, the levels of which limit the nutritional quality of plants. Both amino acids diverge from the same branch of the aspartate family biosynthesis pathway; therefore, their biosynthesis pathways compete for the same carbon/amino substrate. To further elucidate the regulation of methionine biosynthesis and seek ways of increasing the levels of these two amino acids, we crossed transgenic tobacco plants overexpressing the bacterial feedback-insensitive aspartate kinase (bAK), containing a significantly higher threonine level, with plants overexpressing Arabidopsis cystathionine gamma-synthase (AtCGS), the first unique enzyme of methionine biosynthesis. Plants co-expressing bAK and the full-length AtCGS (F-AtCGS) have significantly higher methionine and threonine levels compared with the levels found in wild-type plants, but the methionine level does not increase beyond that found in plants expressing F-AtCGS alone. This finding can be explained through the feedback inhibition regulation mediated by the methionine metabolite on the transcript level of AtCGS. To test this assumption, plants expressing bAK were crossed with plants expressing two mutated forms of AtCGS in which the domains responsible for the feedback regulation have been deleted. Indeed, significantly higher methionine contents and its metabolites levels accumulated in the newly produced plants, and the levels of threonine were also significantly higher than in the wild-type plants. The transcript level of the two mutated forms of AtCGS significantly increased when there was a high content of threonine in the plants, suggesting that threonine modulates, probably indirectly, the transcript level of AtCGS.

Amino acid biosynthesis: new architectures in allosteric enzymes

This review focuses on the allosteric controls in the Aspartate-derived and the branched-chain amino acid biosynthetic pathways examined both from kinetic and structural points of view. The objective is to show the differences that exist among the plant and microbial worlds concerning the allosteric regulation of these pathways and to unveil the structural bases of this diversity. Indeed, crystallographic structures of enzymes from these pathways have been determined in bacteria, fungi and plants, providing a wonderful opportunity to obtain insight into the acquisition and modulation of allosteric controls in the course of evolution. This will be examined using two enzymes, threonine synthase and the ACT domain containing enzyme aspartate kinase. In a last part, as many enzymes in these pathways display regulatory domains containing the conserved ACT module, the organization of ACT domains in this kind of allosteric enzymes will be reviewed, providing explanations for the variety of allosteric effectors and type of controls observed.

Regulation of one-carbon metabolism in Arabidopsis: the N-terminal regulatory domain of cystathionine gamma-synthase is cleaved in response to folate starvation

In all organisms, control of folate homeostasis is of vital importance to sustain the demand for one-carbon (C1) units that are essential in major metabolic pathways. In this study we induced folate deficiency in Arabidopsis (Arabidopsis thaliana) cells by using two antifolate inhibitors. This treatment triggered a rapid and important decrease in the pool of folates with significant modification in the distribution of C1-substituted folate coenzymes, suggesting an adaptive response to favor a preferential shuttling of the flux of C1 units to the synthesis of nucleotides over the synthesis of methionine (Met). Metabolic profiling of folate-deficient cells indicated important perturbation of the activated methyl cycle because of the impairment of Met synthases that are deprived of their substrate 5-methyl-tetrahydrofolate. Intriguingly, S-adenosyl-Met and Met pools declined during the initial period of folate starvation but were further restored to typical levels. Reestablishment of Met and S-adenosyl-Met homeostasis was concomitant with a previously unknown posttranslational modification that consists in the removal of 92 amino acids at the N terminus of cystathionine gamma-synthase (CGS), the first specific enzyme for Met synthesis. Rescue experiments and analysis of different stresses indicated that CGS processing is specifically associated with perturbation of the folates pool. Also, CGS processing involves chloroplastic serine-type proteases that are expressed in various plant species subjected to folate starvation. We suggest that a metabolic effector, to date unidentified, can modulate CGS activity in vivo through an interaction with the N-terminal domain of the enzyme and that removal of this domain can suppress this regulation.

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 simulation of an in vitro multi-enzyme system

Parameters often are tuned with metabolite concentration time series data to build a dynamic model of metabolism. However, such tuning may reduce the extrapolation ability (generalization capability) of the model. In this study, we determined detailed kinetic parameters of three purified Escherichia coli glycolytic enzymes using the initial velocity method for individual enzymes; i.e., the parameters were determined independently from metabolite concentration time series data. The metabolite concentration time series calculated by the model using the parameters matched the experimental data obtained in an actual multi-enzyme system consisting of the three purified E. coli glycolytic enzymes. Thus, the results indicate that kinetic parameters can be determined without using an undesirable tuning process.

Modeling of enzyme production kinetics

Models of single cells, cell populations, and cultures can be most useful in organizing information in a comprehensive system description, as well as in optimizing and controlling actual production operations. Models discussed in this article are of various degrees of biological structure and mathematical complexity. The models are developed based on the biomass formation, substrate consumption, and product formation. the potentials asn the limitations of all the models have been reported. The parameter estimation by different methods has been discussed in this communication. These parameters will be helpful to explore the areas where future-modeling studies may be especially valuable.

Allosteric monofunctional aspartate kinases fromArabidopsis

Plant monofunctional aspartate kinase is unique among all aspartate kinases, showing synergistic inhibition by lysine and S-adenosyl-l-methionine (SAM). The Arabidopsis genome contains three genes for monofunctional aspartate kinases. We show that aspartate kinase 2 and aspartate kinase 3 are inhibited only by lysine, and that aspartate kinase 1 is inhibited in a synergistic manner by lysine and SAM. In the absence of SAM, aspartate kinase 1 displayed low apparent affinity for lysine compared to aspartate kinase 2 and aspartate kinase 3. In the presence of SAM, the apparent affinity of aspartate kinase 1 for lysine increased considerably, with K(0.5) values for lysine inhibition similar to those of aspartate kinase 2 and aspartate kinase 3. For all three enzymes, the inhibition resulted from an increase in the apparent K(m) values for the substrates ATP and aspartate. The mechanism of aspartate kinase 1 synergistic inhibition was characterized. Inhibition by lysine alone was fast, whereas synergistic inhibition by lysine plus SAM was very slow. SAM by itself had no effect on the enzyme activity, in accordance with equilibrium binding analyses indicating that SAM binding to aspartate kinase 1 requires prior binding of lysine. The three-dimensional structure of the aspartate kinase 1-Lys-SAM complex has been solved [Mas-Droux C, Curien G, Robert-Genthon M, Laurencin M, Ferrer JL & Dumas R (2006) Plant Cell18, 1681-1692]. Taken together, the data suggest that, upon binding to the inactive aspartate kinase 1-Lys complex, SAM promotes a slow conformational transition leading to formation of a stable aspartate kinase 1-Lys-SAM complex. The increase in aspartate kinase 1 apparent affinity for lysine in the presence of SAM thus results from the displacement of the unfavorable equilibrium between aspartate kinase 1 and aspartate kinase 1-Lys towards the inactive form.

Methionine catabolism in Arabidopsis cells is initiated by a gamma-cleavage process and leads to S-methylcysteine and isoleucine syntheses

Despite recent progress in elucidating the regulation of methionine (Met) synthesis, little is known about the catabolism of this amino acid in plants. In this article, we present several lines of evidence indicating that the cleavage of Met catalyzed by Met gamma-lyase is the first step in this process. First, we cloned an Arabidopsis cDNA coding a functional Met gamma-lyase (AtMGL), a cytosolic enzyme catalyzing the conversion of Met into methanethiol, alpha-ketobutyrate, and ammonia. AtMGL is present in all of the Arabidopsis organs and tissues analyzed, except in quiescent dry mature seeds, thus suggesting that AtMGL is involved in the regulation of Met homeostasis in various situations. Also, we demonstrated that the expression of AtMGL was induced in Arabidopsis cells in response to high Met levels, probably to bypass the elevated Km of the enzyme for Met. Second, [13C]-NMR profiling of Arabidopsis cells fed with [13C]Met allowed us to identify labeled S-adenosylmethionine, S-methylmethionine, S-methylcysteine (SMC), and isoleucine (Ile). The unexpected production of SMC and Ile was directly associated to the function of Met gamma-lyase. Indeed, we showed that part of the methanethiol produced during Met cleavage could react with an activated form of serine to produce SMC. The second product of Met cleavage, alpha-ketobutyrate, entered the pathway of Ile synthesis in plastids. Together, these data indicate that Met catabolism in Arabidopsis cells is initiated by a gamma-cleavage process and can result in the formation of the essential amino acid Ile and a potential storage form for sulfide or methyl groups, SMC.

Simplified modelling of metabolic pathways for flux prediction and optimization: lessons from an in vitro reconstruction of the upper part of glycolysis

Explicit modelling of metabolic networks relies on well-known mathematical tools and specialized computer programs. However, identifying and estimating the values of the very numerous enzyme parameters inherent to the models remain a tedious and difficult task, and the rate equations of the reactions are usually not known in sufficient detail. A way to circumvent this problem is to use 'non-mechanistic' models, which may account for the behaviour of the systems with a limited number of parameters. Working on the first part of glycolysis reconstituted in vitro, we showed how to derive, from titration experiments, values of effective enzyme activity parameters that do not include explicitly any of the classical kinetic constants. With a maximum of only two parameters per enzyme, this approach produced very good estimates for the flux values, and enabled us to determine the optimization conditions of the system, i.e. to calculate the set of enzyme concentrations that maximizes the flux. This fast and easy method should be valuable in the context of integrative biology or for metabolic engineering, where the challenge is to deal with the dramatic increase in the number of parameters when the systems become complex.

The aspartic acid metabolic pathway, an exciting and essential pathway in plants

Aspartate is the common precursor of the essential amino acids lysine, threonine, methionine and isoleucine in higher plants. In addition, aspartate may also be converted to asparagine, in a potentially competing reaction. The latest information on the properties of the enzymes involved in the pathways and the genes that encode them is described. An understanding of the overall regulatory control of the flux through the pathways is undisputedly of great interest, since the nutritive value of all cereal and legume crops is reduced due to low concentrations of at least one of the aspartate-derived amino acids. We have reviewed the recent literature and discussed in this paper possible methods by which the concentrations of the limiting amino acids may be increased in the seeds.

Allosteric threonine synthase. Reorganization of the pyridoxal phosphate site upon asymmetric activation through S-adenosylmethionine binding to a novel site

Threonine synthase (TS) is a fold-type II pyridoxal phosphate (PLP)-dependent enzyme that catalyzes the ultimate step of threonine synthesis in plants and microorganisms. Unlike the enzyme from microorganisms, plant TS is activated by S-adenosylmethionine (AdoMet). The mechanism of activation has remained unknown up to now. We report here the crystallographic structures of Arabidopsis thaliana TS in complex with PLP (aTS) and with PLP and AdoMet (aTS-AdoMet), which show with atomic detail how AdoMet activates TS. The aTS structure reveals a PLP orientation never previously observed for a type II PLP-dependent enzyme and explains the low activity of plant TS in the absence of its allosteric activator. The aTS-AdoMet structure shows that activation of the enzyme upon AdoMet binding triggers a large reorganization of active site loops in one monomer of the structural dimer and allows the displacement of PLP to its active conformation. Comparison with other TS structures shows that activation of the second monomer may be triggered by substrate binding. This structure also discloses a novel fold for two AdoMet binding sites located at the dimer interface, each site containing two AdoMet effectors bound in tandem. Moreover, aTS-AdoMet is the first structure of an enzyme that uses AdoMet as an allosteric effector.

Towards building the silicon cell: a modular approach

Systems Biology aims to understand quantitatively how properties of biological systems can be understood as functions of the characteristics of, and interactions between their macromolecular components. Whereas, traditional biochemistry focused on isolation and characterization of cellular components, the challenge for Systems Biology lies in integration of this knowledge and the knowledge about molecular interactions. Computer models play an important role in this integration. We here discuss an approach with which we aim to link kinetic models on small parts of metabolism together, so as to form detailed kinetic models of larger chunks of metabolism, and ultimately of the entire living cell. Specifically, we will discuss techniques that can be used to model a sub-network in isolation of a larger network of which it is a part, while still maintaining the dynamics of the larger complete network. We will start by outlining the JWS online system, the silicon cell project, and the type of models we propose. JWS online is a model repository, which can be used for the storage, simulation and analysis of kinetic models. We advocate to integrate a top-down approach, where measurements on the complete system are used to derive fluxes in a detailed structural model, with a bottom-up approach, consisting of the integration of molecular mechanism-based detailed kinetic models into the structural model.

Identification of six novel allosteric effectors of Arabidopsis thaliana aspartate kinase-homoserine dehydrogenase isoforms. Physiological context sets the specificity

The Arabidopsis genome contains two genes predicted to code for bifunctional aspartate kinase-homoserine dehydrogenase enzymes (isoforms I and II). These two activities catalyze the first and the third steps toward the synthesis of the essential amino acids threonine, isoleucine, and methionine. We first characterized the kinetic and regulatory properties of the recombinant enzymes, showing that they mainly differ with respect to the inhibition of the homoserine dehydrogenase activity by threonine. A systematic search for other allosteric effectors allowed us to identify an additional inhibitor (leucine) and 5 activators (alanine, cysteine, isoleucine, serine, and valine) equally efficient on aspartate kinase I activity (4-fold activation). The six effectors of aspartate kinase I were all activators of aspartate kinase II activity (13-fold activation) and displayed a similar specificity for the enzyme. No synergy between different effectors could be observed. The activation, which resulted from a decrease in the Km values for the substrates, was detected using low substrates concentrations. Amino acid quantification revealed that alanine and threonine were much more abundant than the other effectors in Arabidopsis leaf chloroplasts. In vitro kinetics in the presence of physiological concentrations of the seven allosteric effectors confirmed that aspartate kinase I and II activities were highly sensitive to changes in alanine and threonine concentrations. Thus, physiological context rather than enzyme structure sets the specificity of the allosteric control. Stimulation by alanine may play the role of a feed forward activation of the aspartate-derived amino acid pathway in plant.

Are high-lysine cereal crops still a challenge?

The essential amino acids lysine and threonine are synthesized in higher plants via a pathway starting with aspartate that also leads to the formation of methionine and isoleucine. Lysine is one of most limiting amino acids in plants consumed by humans and livestock. Recent genetic, molecular, and biochemical evidence suggests that lysine synthesis and catabolism are regulated by complex mechanisms. Early kinetic studies utilizing mutants and transgenic plants that over-accumulate lysine have indicated that the major step for the regulation of lysine biosynthesis is at the enzyme dihydrodipicolinate synthase. Despite this tight regulation, recent strong evidence indicates that lysine catabolism is also subject to control, particularly in cereal seeds. The challenge of producing crops with a high-lysine concentration in the seeds appeared to be in sight a few years ago. However, apart from the quality protein maize lines currently commercially available, the release of high-lysine crops has not yet occurred. We are left with the question, is the production of high-lysine crops still a challenge?

The expression level of threonine synthase and cystathionine-γ-synthase is influenced by the level of both threonine and methionine in Arabidopsis plants

The biosynthesis pathways of the essential amino acids methionine and threonine diverge from O-phosphohomoserine, an intermediate metabolite in the aspartate family of amino acids. Thus, the enzymes cystathionine-gamma-synthase (CGS) in the methionine pathway and threonine synthase (TS), the last enzyme in the threonine pathway, compete for this common substrate. To study this branching point, we overexpressed TS in sense and antisense orientation in Arabidopsis plants with the aim to study its effect on the level of threonine but more importantly on the methionine content. Positive correlation was found not only between TS expression level and threonine content, but also between TS/threonine and CGS expression level. Plants expressing the sense orientation of TS showed a higher level of threonine, increased expression level of CGS, and a significantly higher level of S-methylmethionine, the transport form of methionine. By contrast, plants expressing the antisense form of TS showed lower levels of threonine and of CGS expression level. In these antisense plants, the methionine level increased up to 47-fold compared to wild-type plants. To study further the effect of threonine on CGS expression level, wild-type plants were irrigated with threonine and control plants were irrigated with methionine or water. While threonine increased the expression level of CGS but reduced that of TS, methionine reduced the expression level of CGS but increased that of TS. This data demonstrate that both methionine and threonine affect the two enzymes at the branching point, thus controlling not only their own level, but also the level of each other. This mechanism probably aids in keeping the levels of these two essential amino acids sufficiently high to support plant growth.

The Silicon Cell initiative: working towards a detailed kinetic description at the cellular level

The Silicon Cell initiative aims to understand cellular systems on the basis of the characteristics of their components. As a tool to achieve this, detailed kinetic models at the network reaction level are being constructed. Such detailed kinetic models are extremely useful for medical and biotechnological applications and form strong tools for fundamental studies. Several recently constructed detailed kinetic models on metabolism (glycolysis), signal transduction (EGF receptor), and the eukaryotic cell cycle (Saccharomyces cerevisiae) have been used to exemplify the Silicon Cell project. These models are stored and made accessible via the JWS Online Cellular Systems Modeling project, a web-based repository of kinetic models. Using a web-browser the models can be interrogated via a user-friendly graphical interface. The goal of the two projects is to combine models on parts of cellular systems and ultimately to construct detailed kinetic models at the cellular level.

Methionine and threonine synthesis are limited by homoserine availability and not the activity of homoserine kinase in Arabidopsis thaliana

Homoserine kinase (HSK) produces O-phospho-l-homoserine (HserP) used by cystathionine gamma-synthase (CGS) for Met synthesis and threonine synthase (TS) for Thr synthesis. The effects of overexpressing Arabidopsis thaliana HSK, CGS, and Escherichia coli TS (eTS), each controlled by the 35S promoter, were compared. The results indicate that in Arabidopsis Hser supply is the major factor limiting the synthesis of HserP, Met and Thr. HSK is not limiting and CGS or TS control the partitioning of HserP. HSK overexpression had no effect on the level of soluble HserP, Met or Thr, however, when treated with Hser these plants produced far more HserP than wild type. Met and Thr also accumulated markedly after Hser treatment but the increase was similar in HSK overexpressing and wild-type plants. CGS overexpression was previously shown to increase Met content, but had no effect on Thr. After Hser treatment Met accumulation increased in CGS-overexpressing plants compared with wild type, whereas HserP declined and Thr was unaffected. Arabidopsis responded differentially to eTS expression depending on the level of the enzyme. At the highest eTS level the Thr content was not increased, but the phenotype was negatively affected and the T1 plants died before reproducing. Comparatively low eTS did not affect phenotype or Thr/Met level, however after Hser treatment HserP and Met accumulation were reduced compared with wild type and Thr was increased slightly. At intermediate eTS activity seedling growth was retarded unless Met was supplied and CGS expression was induced, indicating that eTS limited HserP availability for Met synthesis.

Methionine metabolism in plants: chloroplasts are autonomous for de novo methionine synthesis and can import S-adenosylmethionine from the cytosol

The subcellular distribution of Met and S-adenosylmethionine (AdoMet) metabolism in plant cells discloses a complex partition between the cytosol and the organelles. In the present work we show that Arabidopsis contains three functional isoforms of vitamin B(12)-independent methionine synthase (MS), the enzyme that catalyzes the methylation of homocysteine to Met with 5-methyltetrahydrofolate as methyl group donor. One MS isoform is present in chloroplasts and is most likely required to methylate homocysteine that is synthesized de novo in this compartment. Thus, chloroplasts are autonomous and are the unique site for de novo Met synthesis in plant cells. The additional MS isoforms are present in the cytosol and are most probably involved in the regeneration of Met from homocysteine produced in the course of the activated methyl cycle. Although Met synthesis can occur in chloroplasts, there is no evidence that AdoMet is synthesized anywhere but the cytosol. In accordance with this proposal, we show that AdoMet is transported into chloroplasts by a carrier-mediated facilitated diffusion process. This carrier is able to catalyze the uniport uptake of AdoMet into chloroplasts as well as the exchange between cytosolic AdoMet and chloroplastic AdoMet or S-adenosylhomocysteine. The obvious function for the carrier is to sustain methylation reactions and other AdoMet-dependent functions in chloroplasts and probably to remove S-adenosylhomocysteine generated in the stroma by methyltransferase activities. Therefore, the chloroplastic AdoMet carrier serves as a link between cytosolic and chloroplastic one-carbon metabolism.