Phenotypic characterization of <Emphasis Type="Italic">Corynebacterium glutamicum</Emphasis> using elementary modes towards synthesis of amino acids

Elementary flux mode (EFM) analysis is a powerful tool to represent the metabolic network structure and can be further utilized for flux analysis. The method enables characterization and quantification of feasible phenotypes in microbes. EFM analysis was employed to characterize the phenotype of Corynebacterium glutamicum to yield various amino acids. The metabolic network of C. glutamicum yielded 62 elementary modes by incorporating the accumulation of amino acids namely, lysine, alanine, valine, glutamine and glutamate. The analysis also allowed us to compute the maximum theoretical yield for the synthesis of various amino acids. These 62 elementary modes were further used to obtain optimal phenotypic space towards accumulation of biomass and lysine. The study indicated that the optimal solution space from 62 elementary modes forms a super space which incorporates various mutants including lysine producing strain of C. glutamicum. The analysis was also extended to obtain sensitivity of the network to variation in the stoichiometry of NADP in the definition of biomass.     

Thermodynamic limitations of PHB production from formate and fructose in Cupriavidus necator

The chemolithotroph Cupriavidus necator H16 is known as a natural producer of the bioplastic-polymer PHB, as well as for its metabolic versatility to utilize different substrates, including formate as the sole carbon and energy source. Depending on the entry point of the substrate, this versatility requires adjustment of the thermodynamic landscape to maintain sufficiently high driving forces for biological processes. Here we employed a model of the core metabolism of C. necator H16 to analyze the thermodynamic driving forces and PHB yields from formate for different metabolic engineering strategies. For this, we enumerated elementary flux modes (EFMs) of the network and evaluated their PHB yields as well as thermodynamics via Max-min driving force (MDF) analysis and random sampling of driving forces. A heterologous ATP:citrate lyase reaction was predicted to increase driving force for producing acetyl-CoA. A heterologous phosphoketolase reaction was predicted to increase maximal PHB yields as well as driving forces. These enzymes were then verified experimentally to enhance PHB titers between 60 and 300% in select conditions. The EFM analysis also revealed that PHB production from formate may be limited by low driving forces through citrate lyase and aconitase, as well as cofactor balancing, and identified additional reactions associated with low and high PHB yield. Proteomics analysis of the engineered strains confirmed an increased abundance of aconitase and cofactor balancing. The findings of this study aid in understanding metabolic adaptation. Furthermore, the outlined approach will be useful in designing metabolic engineering strategies in other non-model bacteria.     

Study of metabolic network of Cupriavidus necator DSM 545 growing on glycerol by applying elementary flux modes and yield space analysis

A metabolic network consisting of 48 reactions was established to describe intracellular processes during growth and poly-3-hydroxybutyrate (PHB) production for Cupriavidus necator DSM 545. Glycerol acted as the sole carbon source during exponential, steady-state cultivation conditions. Elementary flux modes were obtained by the program Metatool and analyzed by using yield space analysis. Four sets of elementary modes were obtained, depending on whether the pair NAD/NADH or FAD/FADH₂ contributes to the reaction of glycerol-3-phosphate dehydrogenase (GLY-3-P DH), and whether 6-phosphogluconate dehydrogenase (6-PG DH) is present or not. Established metabolic network and the related system of equations provide multiple solutions for the simultaneous synthesis of PHB and biomass; this number of solutions can be further increased if NAD/NADH or FAD/FADH₂ were assumed to contribute in the reaction of GLY-3-P DH. As a major outcome, it was demonstrated that experimentally determined yields for biomass and PHB with respect to glycerol fit well to the values obtained in silico when the Entner–Doudoroff pathway (ED) dominates over the glycolytic pathway; this is also the case if the Embden–Meyerhof–Parnas pathway dominates over the ED.     

Design, construction and performance of the most efficient biomass producing E. coli bacterium

Inverse metabolic engineering based on elementary mode analysis was applied to maximize the biomass yield of Escherchia coli MG1655. Elementary mode analysis was previously employed to identify among 1691 possible pathways for cell growth the most efficient pathway with maximum biomass yield. The metabolic network analysis predicted that deletion of only 6 genes reduces the number of possible elementary modes to the most efficient pathway. We have constructed a strain containing these gene deletions and we evaluated its properties in batch and in chemostat growth experiments. The results show that the theoretical predictions are closely matched by the properties of the designed strain.     

基于Petri网T不变量的PHB新陈代谢建模分析(英文)

从拓扑结构的角度分析生化反应网络是生物信息学研究中的一个热点问题。通过将两种传统的途径分析方法(基元模式和极端途径)与Petri网的T不变量分析进行了比较,结果表明:它们本质上是一致的,但是采用Petri网的T不变量分析更便捷。然后,利用Petri网技术构建了PHB代谢模型。对该模型作了结构分析,将计算得到的23个T不变量进行了分组:I组表示简单的可逆反应,II组表示循环的反应,III组可用于调控ATP/ADP比率,IV组是与PHB生产直接相关的反应,可用于代谢工程以提高PHB的产率。最后讨论了Petri网的T不变量分析在这个领域中的应用。     

Kinetic studies and biochemical pathway analysis of anaerobic poly-(R)-3-hydroxybutyric acid synthesis in Escherichia coli

Poly-(R)-3-hydroxybutyric acid (PHB) was synthesized anaerobically in recombinant Escherichia coli. The host anaerobically accumulated PHB to more than 50% of its cell dry weight during cultivation in either growth or nongrowth medium. The maximum specific PHB production rate during growth-associated synthesis was approximately 2.3 +/- 0.2 mmol of PHB/g of residual cell dry weight/h. The by-product secretion profiles differed significantly between the PHB-synthesizing strain and the control strain. PHB production decreased acetate accumulation for both growth and nongrowth-associated PHB synthesis. For instance under nongrowth cultivation, the PHB-synthesizing culture produced approximately 66% less acetate on a glucose yield basis as compared to a control culture. A theoretical biochemical network model was used to provide a rational basis to interpret the experimental results like the fermentation product secretion profiles and to study E. coli network capabilities under anaerobic conditions. For example, the maximum theoretical carbon yield for anaerobic PHB synthesis in E. coli is 0.8. The presented study is expected to be generally useful for analyzing, interpreting, and engineering cellular metabolisms.     

Kinetic Studies and Biochemical Pathway Analysis of Anaerobic Poly-(R)-3-Hydroxybutyric Acid Synthesis in Escherichia coli

Poly-(R)-3-hydroxybutyric acid (PHB) was synthesized anaerobically in recombinant Escherichia coli. The host anaerobically accumulated PHB to more than 50% of its cell dry weight during cultivation in either growth or nongrowth medium. The maximum specific PHB production rate during growth-associated synthesis was approximately 2.3 ± 0.2 mmol of PHB/g of residual cell dry weight/h. The by-product secretion profiles differed significantly between the PHB-synthesizing strain and the control strain. PHB production decreased acetate accumulation for both growth and nongrowth-associated PHB synthesis. For instance under nongrowth cultivation, the PHB-synthesizing culture produced approximately 66% less acetate on a glucose yield basis as compared to a control culture. A theoretical biochemical network model was used to provide a rational basis to interpret the experimental results like the fermentation product secretion profiles and to study E. coli network capabilities under anaerobic conditions. For example, the maximum theoretical carbon yield for anaerobic PHB synthesis in E. coli is 0.8. The presented study is expected to be generally useful for analyzing, interpreting, and engineering cellular metabolisms.     

Combinatorial complexity of pathway analysis in metabolic networks

Elementary flux mode analysis is a promising approach for a pathway-oriented perspective of metabolic networks. However, in larger networks it is hampered by the combinatorial explosion of possible routes. In this work we give some estimations on the combinatorial complexity including theoretical upper bounds for the number of elementary flux modes in a network of a given size. In a case study, we computed the elementary modes in the central metabolism of Escherichia coli while utilizing four different substrates. Interestingly, although the number of modes occurring in this complex network can exceed half a million, it is still far below the upper bound. Hence, to a certain extent, pathway analysis of central catabolism is feasible to assess network properties such as flexibility and functionality.     

Metabolic pathway structures for recombinant protein synthesis in Escherichia coli

Escherichia coli is a valuable commercial host for the production of heterologous proteins. We used elementary mode analysis to identify all possible genetically independent pathways for the production of three specific recombinant proteins, green fluorescent protein, savinase and an artificial protein consisting of repeating units of a five-amino-acid cassette. Analysis of these pathways led to the identification of the most efficient pathways for the production of each of these proteins. The results indicate that the amino acid composition of expressed proteins has a profound effect on the number and identity of possible pathways for the production of these proteins. We show that several groups of elementary modes produce the same ratio of biomass and recombinant protein. The pattern of occurrence of these modes is dependent on the amino acid composition of the specific foreign protein produced. These pathways are formed as systemic combinations of other pathways that produce biomass or foreign protein alone after the elimination of fluxes in specific internal reversible reactions or the reversible carbon dioxide exchange reaction. Since these modes represent pathway options that enable the cell to produce biomass and protein without utilizing these reactions, removal of these reactions would constrain the cells to utilize these modes for producing biomass and foreign protein at constant ratios.     

Supply-Side Analysis of Growth of Bacillus subtilis on Glucose-Citrate Medium: Feasible Network Alternatives and Yield Optimality

Our prior work revealed that compared to the case for glucose metabolism, increased carbon yield and nil acid formation result when Bacillus subtilis grows on glucose medium containing citrate. To scrutinize further how citrate addition may alter metabolic flux regulation and the degree that the observed carbon yield corresponds to the maximal value, experimental (by least-squares analysis) and optimal (by linear programming) fluxes and yields were contrasted. Networks with differing reaction routes, directionality constraints, and transhydrogenase activities were examined. To attain an elevated carbon yield, citrate-glucose utilization need not alleviate any stoichiometric constraints that can sometimes interfere with the attainment of network objectives. Rather, the high carbon yield and nil acid formation attained may be linked to restriction of glycolytic capacity, particularly at the level of pyruvate kinase, which is consistent with a hypothesized effect of coupled metal-citrate uptake. Allowing for malic enzyme activity, hexose monophosphate pathway cycling, and transhydrogenase activity may also lead to the flux distributions underlying the high carbon yield observed. Finally, the observed carbon yield corresponded well to the maximum yield provided by all the network alternatives examined. Collectively, these results suggest that (i) the observed carbon yield is essentially equal to the maximal values associated with plausible networks and (ii), as suggested by others, nonoptimal flux regulation may contribute significantly to apparent cellular maintenance requirements.     

Metabolic flux distribution and thermodynamic analysis of green fluorescent protein production in recombinant Escherichia coli: The effect of carbon source and CO 2 partial pressure

Increasing recombinant protein production yields from bacterial cultures remains an important challenge in biotechnology. Acetate accumulation due to high dissolved carbon dioxide (pCO₂) concentrations in the medium has been identified as a factor that negatively affects such yields. Under appropriate culture conditions, acetate could be re-assimilated by bacterial cells to maintain heterologous proteins production. In this work, we developed a simplified metabolic network aiming to establish a reaction rate analysis for a recombinant Escherichia coli when producing green fluorescent protein (GFP) under controlled pCO₂ concentrations. Because E. coli is able to consume both glucose and acetate, the analysis was performed in two stages. Our results indicated that GFP synthesis is an independent process of cellular growth in some culture phases. Additionally, recombinant protein production is influenced by the available carbon source and the amount of pCO₂ in the culture medium. When growing on glucose, the increase in the pCO₂ concentration produced a down-regulation of central carbon metabolism by directing the carbon flux toward acetate accumulation; as a result, cellular growth and the overall GFP yield decreased. However, the maximum specific rate of GFP synthesis occurred with acetate as the main available carbon source, despite the low activity in the other metabolic pathways. To maintain cellular functions, including GFP synthesis, carbon flux was re-distributed toward the tricarboxylic acid cycle and the pentose phosphate pathway to produce ATP and NADH. The thermodynamic analysis allowed demonstrating the feasibility of the simplified network for describing the metabolic state of a recombinant system.     

Hydrocortisone made in yeast: metabolic engineering turns a unicellular microorganism into a drug-synthesizing factory

Inspired by the successful work of converting Saccharomyces cerevisiae into an microorganism capable of synthesizing hydrocortisone, a 27-carbon molecule, from ethanol, a 2-carbon molecule, this review provides an overview of the potential of yeast as a recombinant organism in the 21st century. Yeast has been used by man for more than 6,000 years, and is still paving the way to new discoveries. It was the first eukaryotic organism to be sequenced, in 1996, and the first to produce hydrocortisone in 2003. In addition, extensive genome-wide analyses have been performed with yeast. In this review, we discuss the pros and cons of using yeast to produce small therapeutic molecules. It is obvious that S. cerevisiae has a cutting edge advantage of being a well-known organism and time will tell if yeast "biohydrocortisone" is a unique example or the beginning of a long list of yeast bioproducts. Other organisms, such as plants and bacteria, are competing with yeast. Bacteria produce a wealth of marketed molecules and plants are capable of producing extremely complex molecules with an unbeatable yield. However, S. cerevisiae offers a unique mix of the simplicity of a recombinant organism combined with the complexity of a eukaryote.     

ACoM: A classification method for elementary flux modes based on motif finding

Elementary flux mode analysis is a powerful tool for the theoretical study of metabolic networks. However, when the networks are complex, the determination of elementary flux modes leads to combinatorial explosion of their number which prevents from drawing simple conclusions from their analysis. To deal with this problem we have developed a method based on the Agglomeration of Common Motifs (ACoM) for classifying elementary flux modes. We applied this algorithm to describe the decomposition into elementary flux modes of the central carbon metabolism in Bacillus subtilis and of the yeast mitochondrial energy metabolism. ACoM helps to give biological meaning to the different elementary flux modes and to the relatedness between reactions. ACoM, which can be viewed as a bi-clustering method, can be of general use for sets of vectors with values 0, +1 or −1.     

Rat ATP citrate-lyase. Molecular cloning and sequence analysis of a full-length cDNA and mRNA abundance as a function of diet, organ, and age

ATP citrate-lyase is the primary enzyme responsible for the synthesis of cytosolic acetyl-CoA in many tissues. We have isolated a full-length cDNA copy of 4.3 kilobase pairs encoding the ATP-citrate lyase mRNA by screening rat liver cDNA library using oligonucleotide probes designed from peptide sequences obtained from the purified rat enzyme. Expression of this cDNA in bacteria, followed by immunoblotting with antibody directed against the ATP citrate-lyase, further demonstrated the identity of this clone. Nucleic acid sequence data indicate that the cDNA contains the complete coding region for the enzyme, which is 1100 amino acids in length with a calculated molecular weight of 121,293. RNA blot analysis indicated an mRNA species of about 4.3 kilobase pairs in livers of chow-fed rats. Rats maintained on low fat, high carbohydrate diets exhibited a striking increase (50-fold) in the level of liver ATP citrate-lyase mRNA as compared with the control animals maintained on a normal diet. The tissue distribution of this mRNA in chow-fed animals revealed a relatively high abundance of the message in liver and adrenal, moderate levels were found in lung, brain, and large intestine with only trace amounts of the message in small intestine, stomach, testis, spleen, pancreas, kidney, and heart. During rat development, the ATP citrate-lyase mRNA was relatively high in the liver at parturition, followed by a reduction in its level during suckling. Higher amounts of the mRNA were detected again in adult animals. The isolation and characterization of the mRNA for ATP citrate-lyase will allow further studies on the reaction mechanism and metabolic regulation of this key enzyme in lipogenesis and cholesterogenesis.     

Mixed glucose and lactate uptake by Corynebacterium glutamicum through metabolic engineering

The Corynebacterium glutamicum ATCC 13032 lysC(fbr) strain was engineered to grow fast on racemic mixtures of lactate and to secrete lysine during growth on lactate as well as on mixtures of lactate and glucose. The wild-type C. glutamicum only grows well on L-lactate. Overexpression of D-lactate dehydrogenase (dld) achieved by exchanging the native promoter of the dld gene for the stronger promoter of the sod gene encoding superoxide dismutase in C. glutamicum resulted in a duplication of biomass yield and faster growth without any secretion of lysine. Elementary mode analysis was applied to identify potential targets for lysine production from lactate as well as from mixtures of lactate and glucose. Two targets for overexpression were pyruvate carboxylase and malic enzyme. The overexpression of these genes using again the sod promoter resulted in growth-associated production of lysine with lactate as sole carbon source with a carbon yield of 9% and a yield of 15% during growth on a lactate-glucose mixture. Both substrates were taken up simultaneously with a slight preference for lactate. As surmised from the elementary mode analysis, deletion of glucose-6-phosphate isomerase resulted in a decreased production of lysine on the mixed substrate. Elementary mode analysis together with suitable objective functions has been found a very useful tool guiding the design of strains producing lysine on mixed substrates.     

The effect of dissolved oxygen on PHB accumulation in activated sludge cultures

Nitrogen removal from wastewater is often limited by the availability of reducing power to perform denitrification, especially when treating wastewaters with a low carbon:nitrogen ratio. In the increasingly popular sequencing batch reactor (SBR), bacteria have the opportunity to preserve reducing power from incoming chemical oxygen demand (COD) as poly-beta-hydroxybutyrate (PHB). The current study uses laboratory experiments and mathematical modeling in an attempt to generate a better understanding of the effect of oxygen on microbial conversion of COD into PHB. Results from a laboratory SBR with acetate as the organic carbon source showed that the aerobic acetate uptake process was oxygen-dependent, producing higher uptake rates at higher dissolved oxygen (DO) supply rates. However, at the lower DO supply rates (k(L)a 6 to 16 h(-1), 0 mg L(-1) DO), a higher proportion of the substrate was preserved as PHB than at higher DO supply rates (k(L)a 30, 51 h(-1), DO >0.9 mg L(-1)). Up to 77% of the reducing equivalents available from acetate were converted to PHB under oxygen limitation (Y(PHB/Ac) 0.68 Cmol/Cmol), as opposed to only 54% under oxygen-excess conditions (Y(PHB/Ac) 0.48 Cmol/Cmol), where a higher fraction of acetate was used for biomass growth. It was calculated that, by oxygen management during the feast phase, the amount of PHB preserved (1.4 Cmmol L(-1) PHB) accounted for an additional denitrification potential of up to 18 mg L(-1) nitrate-nitrogen. The trends of the effect of oxygen (and hence ATP availability) on PHB accumulation could be reproduced by the simulation model, which was based on biochemical stoichiometry and maximum rates obtained from experiments. Simulated data showed that, at low DO concentrations, the limited availability of adenosine triphosphate (ATP) prevented significant biomass growth and most ATP was used for acetate transport into the cell. In contrast, high DO supply rates provided surplus ATP and hence higher growth rates, resulting in decreased PHB yields. The results suggest that oxygen management is crucial to conserving reducing power during the feast phase of SBR operation, as excessive aeration rates decrease the PHB yield and allow higher biomass growth.     

Streptomyces lividans potassium channel contains poly-(R)-3-hydroxybutyrate and inorganic polyphosphate

The Streptomyces lividans KcsA potassium channel, a homotetramer of 17.6 kDa subunits, was found to contain two nonproteinaceous polymers, namely, poly-(R)-3-hydroxybutyrate (PHB) and inorganic polyphosphate (polyP). PHB and polyP are ubiquitous cellular constituents with a demonstrated capacity for cation selection and transport. PHB was detected in both tetramer and monomer species of KcsA by reaction to anti-PHB IgG on Western blots, and estimated as 28 monomer units of PHB per KcsA tetramer by a chemical assay in which PHB is converted to its unique degradation product, crotonic acid. PolyP was detected in KcsA tetramers, but not in monomers, by metachromatic reaction to o-toluidine blue stain on SDS-PAGE gels. A band of free polyP was also visible, suggesting that polyP is released when tetramers dissociate. The exopolyphosphatase of Saccharomyces cerevisiae degraded the free polyP, but tetramer-associated polyP was not affected, indicating it was inaccessible to the enzyme. PolyP in KcsA was estimated as 15 monomer units per tetramer by an enzymatic assay in which polyphosphate kinase is used to transfer phosphates from polyP to [(14)C]ADP, yielding [(14)C]ATP. The experimentally determined isoelectric point of KcsA tetramer was 6.5-7.5, substantially more acidic than the theoretical pI of 10.3, and consistent with the inclusion of a polyanion. The results suggest that PHB is covalently bound to KcsA subunits while polyP is held within tetramers by ionic forces. It is posited that KcsA protein creates an environment in which PHB/polyP is selective for K(+). The basic amino acids attenuate the negative charge density of polyP, thereby transforming the cation binding preference from multivalent to monovalent, and discrimination between K(+) and Na(+) is accomplished by adjusting the ligand geometry in cation binding cavities formed by PHB and polyP.     

Analysis of NADPH supply during xylitol production by engineered Escherichia coli

Escherichia coli strain PC09 (DeltaxylB, cAMP-independent CRP (crp*) mutant) expressing an NADPH-dependent xylose reductase from Candida boidinii (CbXR) was previously reported to produce xylitol from xylose while metabolizing glucose [Cirino et al. (2006) Biotechnol Bioeng 95(6): 1167-1176]. This study aims to understand the role of NADPH supply in xylitol yield and the contribution of key central carbon metabolism enzymes toward xylitol production. Studies in which the expression of CbXR or a xylose transporter was increased suggest that enzyme activity and xylose transport are not limiting xylitol production in PC09. A constraints-based stoichiometric metabolic network model was used to understand the roles of central carbon metabolism reactions and xylose transport energetics on the theoretical maximum molar xylitol yield (xylitol produced per glucose consumed), and xylitol yields (Y(RPG)) were measured from resting cell biotransformations with various PC09 derivative strains. For the case of xylose-proton symport, omitting the Zwf (glucose-6-phosphate dehydrogenase) or PntAB (membrane-bound transhydrogenase) reactions or TCA cycle activity from the model reduces the theoretical maximum yield from 9.2 to 8.8, 3.6, and 8.0 mol xylitol (mol glucose)(-1), respectively. Experimentally, deleting pgi (encoding phosphoglucose isomerase) from strain PC09 improves the yield from 3.4 to 4.0 mol xylitol (mol glucose)(-1), while deleting either or both E. coli transhydrogenases (sthA and pntA) has no significant effect on the measured yield. Deleting either zwf or sucC (TCA cycle) significantly reduces the yield from 3.4 to 2.0 and 2.3 mol xylitol (mol glucose)(-1), respectively. Expression of a xylose reductase with relaxed cofactor specificity increases the yield to 4.0. The large discrepancy between theoretical maximum and experimentally determined yield values suggests that biocatalysis is compromised by pathways competing for reducing equivalents and dissipating energy. The metabolic role of transhydrogenases during E. coli biocatalysis has remained largely unspecified. Our results demonstrate the importance of direct NADPH supply by NADP+-utilizing enzymes in central metabolism for driving heterologous NADPH-dependent reactions, and suggest that the pool of reduced cofactors available for biotransformation is not readily interchangeable via transhydrogenase.