The threonine degradation pathway of the Trypanosoma brucei procyclic form: the main carbon source for lipid biosynthesis is under metabolic control

The Trypanosoma brucei procyclic form resides within the digestive tract of its insect vector, where it exploits amino acids as carbon sources. Threonine is the amino acid most rapidly consumed by this parasite, however its role is poorly understood. Here, we show that the procyclic trypanosomes grown in rich medium only use glucose and threonine for lipid biosynthesis, with threonine's contribution being ～ 2.5 times higher than that of glucose. A combination of reverse genetics and NMR analysis of excreted end‐products from threonine and glucose metabolism, shows that acetate, which feeds lipid biosynthesis, is also produced primarily from threonine. Interestingly, the first enzymatic step of the threonine degradation pathway, threonine dehydrogenase (TDH, EC 1.1.1.103), is under metabolic control and plays a key role in the rate of catabolism. Indeed, a trypanosome mutant deleted for the phosphoenolpyruvate decarboxylase gene (PEPCK, EC 4.1.1.49) shows a 1.7‐fold and twofold decrease of TDH protein level and activity, respectively, associated with a 1.8‐fold reduction in threonine‐derived acetate production. We conclude that TDH expression is under control and can be downregulated in response to metabolic perturbations, such as in the PEPCK mutant in which the glycolytic metabolic flux was redirected towards acetate production.     

Controle de la chaine de biosynthese de la threonine chezE. coli

Resume Le travail décrit dans cet article s'appuie sur la théorie du contrôle du métabolisme et a pour but la mesure des coefficients de contrôle des différentes étapes sur le flux de production de thréonine. Le coefficient de contrôle d'une étape sur un flux mesure quantitativement la réponse du flux aux variations de l'étape. Cette notion est donc particulièrement importante aussi bien dans les situations pathologiques (diminution de l'activité d'une étape) qu'en biotechnologies où au contraire les étapes sont amplifiées.La mesure des coefficients de contrôle des étapes d'une chaîne métabolique permet donc de connaître celle(s) dont l'amplification doit entraîner une augmentation concomitante du flux.Nous avons appliqué ces concepts à l'étude de la voie de biosynthèse de la thréonine à partir de l'aspartate.La voie de la biosynthèse de la thréonine à partir de l'aspartate est constituée de cinq étapes catalysées par cinq activités enzymatiques: l'aspartokinase (AK), l'aspartate semi-aldéhyde déshydrogénase (ASA-DH), l'homosérine déshydrogénase (HDH), l'homosérine kinase (HK) et la thréonine synthase (TS).La mesure du coefficient de contrôle de la première étape (AK, insensible à la rétro-inhibition par la thréonine dans la souche étudiée) a montré qu'elle était faiblement contrôlante. L'étude a révélé la présence d'une inhibition jusqu'alors inconnue de l'homosérine kinase par la lysine.Un début de modélisation de cette chaîne de biosynthèse permet d'expliquer les résultats expérimentaux.

---

This paper deals with the application of the metabolic control theory, especially the measurement of control coefficients, to the threonine pathway inE. coli. The control coefficient of a step on a metabolic flux quantitatively assesses the flux response to the step variations. This concept is particularly relevant both in pathological situations (decrease in the activity of an enzymatic step in the metabolism) and in biotechnologies, where, on the contrary steps are amplified.Measurement of the control coefficients of the steps of a metabolic network makes it possible to know those whose amplification should lead to a simultaneous increase in the fluxes.We have applied these concepts to threonine biosynthesis from aspartate inE. coli. The threonine pathway starting from aspartate involves five steps catalyzed by five enzyme activities: aspartokinase (AK), aspartate-semialdehyde-dehydrogenase (ASA-DH), homoserine dehydrogenase (HDH), homoserine kinase (HK) and hreonine synthetase activity (TS).Measurement of the control coefficient of the first step (AK, insensitive to threonine inhibition in the studied strain) has shown that it controls threonine production weakly. Our study has revealed a hitherto unknown inhibition of homoserine kinase activity by lysine.Mathematical modeling of this metabolic pathway has been undertaken to better understand our experimental results.

     

Stabilization of threonine production in an Escherichia coli overproducing strain using tightly regulated T7 promoter system

The industrial production of compounds by E. coli strains is often accompanied with variability in yield levels. To investigate the mechanism of such instability the over-production of threonine was used as a model. The instability in this strain appears to be caused by a metabolic burden resulting in occurrence of low producing revertants. A successful application of the tightly regulated T7 expression system is presented as a possible solution providing a substantial stabilization of the threonine production.     

A phenomenological model to represent the kinetics of growth by <Emphasis Type="Italic">Corynebacterium glutamicum</Emphasis> for lysine production

Corynebacterium glutamicum is commonly used for lysine production. In the last decade, several metabolic engineering approaches have been successfully applied to C. glutamicum. However, only few studies have been focused on the kinetics of growth and lysine production. Here, we present a phenomenological model that captures the growth and lysine production during different phases of fermentation at various initial dextrose concentrations. The model invokes control coefficients to capture the dynamics of lysine and trehalose synthesis. The analysis indicated that maximum lysine productivity can be obtained using 72 g/L of initial dextrose concentration in the media, while growth was optimum at 27 g/L of dextrose concentration. The predictive capability was demonstrated through a two-stage fermentation strategy to enhance the productivity of lysine by 1.5 times of the maximum obtained in the batch fermentation. Two-stage fermentation indicated that the kinetic model could be further extended to predict the optimal feeding strategy for fed-batch fermentation.     

Dynamic metabolic flux analysis (DMFA): a framework for determining fluxes at metabolic non-steady state

Metabolic flux analysis (MFA) is a key tool for measuring in vivo metabolic fluxes in systems at metabolic steady state. Here, we present a new method for dynamic metabolic flux analysis (DMFA) of systems that are not at metabolic steady state. The advantages of our DMFA method are: (1) time-series of metabolite concentration data can be applied directly for estimating dynamic fluxes, making data smoothing and estimation of average extracellular rates unnecessary; (2) flux estimation is achieved without integration of ODEs, or iterations; (3) characteristic metabolic phases in the fermentation data are identified automatically by the algorithm, rather than selected manually/arbitrarily. We demonstrate the application of the new DMFA framework in three example systems. First, we evaluated the performance of DMFA in a simple three-reaction model in terms of accuracy, precision and flux observability. Next, we analyzed a commercial glucose-limited fed-batch process for 1,3-propanediol production. The DMFA method accurately captured the dynamic behavior of the fed-batch fermentation and identified characteristic metabolic phases. Lastly, we demonstrate that DMFA can be used without any assumed metabolic network model for data reconciliation and detection of gross measurement errors using carbon and electron balances as constraints.     

Enzymology in vivo using NMR and molecular genetics

Models of metabolic flux regulation are frequently based on an extrapolation of the kinetic properties of enzymes measured in vitro to the intact cell. Such an extrapolation assumes a detailed knowledge of the intracellular environment of these enzymes in terms of their free substates and effectors concentrations and possible interaction with other cellular macromolecules, which may modify their kinetic properties. These is a considerable incentive, therefore, to study the properties of enzymes directly in vivo. We have been using non-invasive NMR techniques, in conjunction with molecular genetic manipulation of enzyme levels, to study the kinetic properties of individual enzymes in vivo. We have also developed a novel strategy which has allowed us to monitor, by NMR, the ligand binding properties and mobilities of enzymes in the intact cell. This technique may also allow us to measured the diffusion coefficients of these proteins in the cell. These studies should give new insight into the properties of enzymes in vivo     

Metabolic design based on a coupled gene expression-metabolic network model of tryptophan production in Escherichia coli

The presumably high potential of a holistic design approach for complex biochemical reaction networks is exemplified here for the network of tryptophan biosynthesis from glucose, a system whose components have been investigated thoroughly before. A dynamic model that combines the behavior of the trp operon gene expression with the metabolic network of central carbon metabolism and tryptophan biosynthesis is investigated. This model is analyzed in terms of metabolic fluxes, metabolic control, and nonlinear optimization. We compare two models for a wild-type strain and another model for a tryptophan producer. An integrated optimization of the whole network leads to a significant increase in tryptophan production rate for all systems under study. This enhancement is well above the increase that can be achieved by an optimization of subsystems. A constant ratio of control coefficients on tryptophan synthesis rate has been identified for the models regarding or disregarding trp operon expression. Although we found some examples where flux control coefficients even contradict the trends of enzyme activity changes in an optimized profile, flux control can be used as an indication for enzymes that have to be taken into account in optimization.     

Phosphorylation of the cell division protein GpsB regulates PrkC kinase activity through a negative feedback loop in Bacillus subtilis

Although many membrane Ser/Thr‐kinases with PASTA motifs have been shown to control bacterial cell division and morphogenesis, inactivation of the Ser/Thr‐kinase PrkC does not impact Bacillus subtilis cell division. In this study, we show that PrkC localizes at the division septum. In addition, three proteins involved in cell division/elongation, GpsB, DivIVA and EzrA are required for stimulating PrkC activity in vivo. We show that GpsB interacts with the catalytic subunit of PrkC that, in turn, phosphorylates GpsB. These observations are not made with DivIVA and EzrA. Consistent with the phosphorylated residue previously detected for GpsB in a high‐throughput phosphoproteomic analysis of B. subtilis, we show that threonine 75 is the single PrkC‐mediated phosphorylation site in GpsB. Importantly, the substitution of this threonine by a phospho‐mimetic residue induces a loss of PrkC kinase activity in vivo and a reduced growth under high salt conditions as observed for gpsB and prkC null mutants. Conversely, substitution of threonine 75 by a phospho‐ablative residue does not induce such growth and PrkC kinase activity defects. Altogether, these data show that proteins of the divisome control PrkC activity and thereby phosphorylation of PrkC substrates through a negative feedback loop in B. subtilis.     

Dual autogenous regulatory role of threonine deaminase in Escherichia coli K-12.

Summary We describe the regulatory properties of two strains carrying either the ilvA624 or the ilvA625 mutations, located in the structural gene for threonine deaminase. Crude extracts of both these strains possess a threonine deaminase activity migrating on polyacrylamide gels, differently from the wild type enzyme. Growth studies demonstrate that these mutations do not cause a limitation of isoleucine biosynthesis, suggesting normal catalytic activity of deaminase.A regulatory consequence of the ilvA624 allele is a derepression of the isoleucine-valine biosynthetic enzymes, which is recessive to an ilvA ⁺ allele. The ilvA625 mutation causes a derepression which is dominant in an ilvA625/ilvA ⁺ diploid. We interpret these data assuming that threonine deaminase, previously shown to be an autogenous regulator of the ilv genes, lacks a repressor function in the ilvA624 mutant, while in the ilvA625 mutant it is a better activator than wild type threonine deaminase.The data are discussed in terms of a model requiring that threonine deaminase, or a precursor of it, is in equilibrium between two forms, one being an activator of gene expression and the other being a repressor.     

Metabolic control at the cytosol–mitochondria interface in different growth phases of CHO cells

Metabolism at the cytosol–mitochondria interface and its regulation is of major importance particularly for efficient production of biopharmaceuticals in Chinese hamster ovary (CHO) cells but also in many diseases. We used a novel systems-oriented approach combining dynamic metabolic flux analysis and determination of compartmental enzyme activities to obtain systems level information with functional, spatial and temporal resolution. Integrating these multiple levels of information, we were able to investigate the interaction of glycolysis and TCA cycle and its metabolic control. We characterized metabolic phases in CHO batch cultivation and assessed metabolic efficiency extending the concept of metabolic ratios. Comparing in situ enzyme activities including their compartmental localization with in vivo metabolic fluxes, we were able to identify limiting steps in glycolysis and TCA cycle. Our data point to a significant contribution of substrate channeling to glycolytic regulation. We show how glycolytic channeling heavily affects the availability of pyruvate for the mitochondria. Finally, we show that the activities of transaminases and anaplerotic enzymes are tailored to permit a balanced supply of pyruvate and oxaloacetate to the TCA cycle in the respective metabolic states. We demonstrate that knowledge about metabolic control can be gained by correlating in vivo metabolic flux dynamics with time and space resolved in situ enzyme activities.     

Unravelling the thioesterases responsible for propionate formation in engineered Pseudomonas putida KT2440

Pseudomonas putida KT2440 is becoming a new robust metabolic chassis for biotechnological applications, due to its metabolic versatility, low nutritional requirements and biosafety status. We have previously engineered P. putida KT2440 to be an efficient propionate producer from L-threonine, although the internal enzymes converting propionyl-CoA to propionate are not clear. In this study, we thoroughly investigated 13 genes annotated as potential thioesterases in the KT2440 mutant. One thioesterase encoded by locus tag PP_4975 was verified to be the major contributor to propionate production in vivo. Deletion of PP_4975 significantly decreased propionate production, whereas the performance was fully restored by gene complement. Compared with thioesterase HiYciA from Haemophilus influenza, thioesterase PP_4975 showed a faster substrate conversion rate in vitro. Thus, this study expands our knowledge on acyl-CoA thioesterases in P. putida KT2440 and may also reveal a new target for further engineering the strain to improve propionate production performance.     

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 ²H₅‐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 ²H₅‐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 β‐oxidative and non‐β‐oxidative pathways is highly distributed.     

Dependence of control coefficient distribution on the boundaries of a metabolic system: a generalized analysis of the effects of additional input and output reactions to a linear pathway

Both experimental and theoretical studies of metabolism are likely to relate to a segment that has been isolated for analytical purposes. In practice, it will be embedded in the whole of cellular metabolism. Thus, it is necessary to consider how conclusions about the control of an isolated pathway may be modified in this wider context where the input and output metabolites are considered as variables of cellular metabolism. Here, we analyse the effect of expanding a linear metabolic pathway by adding an extra input or an extra output. In particular, we analyse the effect of the elasticities of the extra steps on control coefficients. We derive matrix algebraic relationships for obtaining flux and concentration control coefficients from expressions depending on these extra elasticities and on parameters (elasticities and control coefficients) of the original pathway. These equations can be shown in certain cases to be generalized versions of earlier rescaling relationships and to be related to top-down analysis, but also apply where the new variable metabolite of the expanded pathway is an effector of more than one step of the original pathway. We use our relationships to analyse the dependence or independence of control coefficients upon these extra elasticities for the published analyses of the pathway of mammalian serine biosynthesis (Fell & Snell, 1988) and Escherischia coli threonine biosynthesis (Chassagnole et al., 2001). The same analysis can be applied to determine whether the transport reactions of substrates and products of a pathway in and out of a cell need to be included in estimations of the control coefficients of the enzymes.     

An In Vitro Procedure for Estimation of Lead Relative Bioavailability: With Validation

The relative bioaccessibility leaching procedure (RBALP) is a simple, reproducible, and rapid in vitro procedure for estimating the in vivo (juvenile swine) relative bioavailability (RBA) of lead in solid media. Control of pH, temperature, and agitation are the most critical parameters of this in vitro procedure. The performance of the method was evaluated by triplicate analyses of each of 19 different test substances by the author and three independent laboratories, and comparison of the results to relative bioavailability (RBA) values measured in vivo. The results indicate that the RBALP measurements are strongly correlated with the in vivo RBA values (r² = 0.924, p < 0.0001), with an average absolute error of 10% and an average predictive error of 20%. Comparison of results within and between laboratories indicates that the procedure is highly reproducible, with inter-and intra-laboratory coefficients of variation of 4% and 6%, respectively, and within-sample precision of approximately 7%. Based on the results reported here, the RBALP can be effective in providing reliable estimates of lead RBA as predicted by the immature swine model. The method may also be valuable in evaluating site-specific differences in bioaccessibility, assessing remedial technologies intended to reduce lead RBA, providing a screening mechanism for futurein vivo studies, and providing insight into the chemical and physical factors that control lead bioavailability.     

An antiphage Escherichia coli mutant for higher production of L-threonine obtained by atmospheric and room temperature plasma mutagenesis

Phage infection is common during the production of L-threonine by E. coli, and low L-threonine production and glucose conversion percentage are bottlenecks for the efficient commercial production of L-threonine. In this study, 20 antiphage mutants producing high concentration of L-threonine were obtained by atmospheric and room temperature plasma (ARTP) mutagenesis, and an antiphage E. coli variant was characterized that exhibited the highest production of L-threonine Escherichia coli ([E. coli] TRFC-AP). The elimination of fhuA expression in E. coli TRFC-AP was responsible for phage resistance. The biomass and cell growth of E. coli TRFC-AP showed no significant differences from those of the parent strain (E. coli TRFC), and the production of L-threonine (159.3 g L⁻¹) and glucose conversion percentage (51.4%) were increased by 10.9% and 9.1%, respectively, compared with those of E. coli TRFC. During threonine production (culture time of 20 h), E. coli TRFC-AP exhibited higher activities of key enzymes for glucose utilization (hexokinase, glucose phosphate dehydrogenase, phosphofructokinase, phosphoenolpyruvate carboxylase, and PYK) and threonine synthesis (glutamate synthase, aspartokinase, homoserine dehydrogenase, homoserine kinase and threonine synthase) compared to those of E. coli TRFC. The analysis of metabolic flux distribution indicated that the flux of threonine with E. coli TRFC-AP reached 69.8%, an increase of 16.0% compared with that of E. coli TRFC. Overall, higher L-threonine production and glucose conversion percentage were obtained with E. coli TRFC-AP due to increased activities of key enzymes and improved carbon flux for threonine synthesis.     

Model‐assisted identification of metabolic engineering strategies for Jatropha curcas lipid pathways

Efficient approaches to increase plant lipid production are necessary to meet current industrial demands for this important resource. While Jatropha curcas cell culture can be used for in vitro lipid production, scaling up the system for industrial applications requires an understanding of how growth conditions affect lipid metabolism and yield. Here we present a bottom‐up metabolic reconstruction of J. curcas supported with labeling experiments and biomass characterization under three growth conditions. We show that the metabolic model can accurately predict growth and distribution of fluxes in cell cultures and use these findings to pinpoint energy expenditures that affect lipid biosynthesis and metabolism. In addition, by using constraint‐based modeling approaches we identify network reactions whose joint manipulation optimizes lipid production. The proposed model and computational analyses provide a stepping stone for future rational optimization of other agronomically relevant traits in J. curcas.     

Harnessing biodiesel-producing microbes: from genetic engineering of lipase to metabolic engineering of fatty acid biosynthetic pathway

Microbial production routes, notably whole-cell lipase-mediated biotransformation and fatty-acids-derived biosynthesis, offer new opportunities for synthesizing biodiesel. They compare favorably to immobilized lipase and chemically catalyzed processes. Genetically modified whole-cell lipase-mediated in vitro route, together with in vivo and ex vivo microbial biosynthesis routes, constitutes emerging and rapidly developing research areas for effective production of biodiesel. This review presents recent advances in customizing microorganisms for producing biodiesel, via genetic engineering of lipases and metabolic engineering (including system regulation) of fatty-acids-derived pathways. Microbial hosts used include Escherichia coli, Saccharomyces cerevisiae, Pichia pastoris and Aspergillus oryzae. These microbial cells can be genetically modified to produce lipases under different forms: intracellularly expressed, secreted or surface-displayed. They can be metabolically redesigned and systematically regulated to obtain balanced biodiesel-producing cells, as highlighted in this study. Such genetically or metabolically modified microbial cells can support not only in vitro biotransformation of various common oil feedstocks to biodiesel, but also de novo biosynthesis of biodiesel from glucose, glycerol or even cellulosic biomass. We believe that the genetically tractable oleaginous yeast Yarrowia lipolytica could be developed to an effective biodiesel-producing microbial cell factory. For this purpose, we propose several engineered pathways, based on lipase and wax ester synthase, in this promising oleaginous host.