Metabolic network analysis ofpenicillium chrysogenumusing13c-labeled glucose

Efforts to install a heterologous Wood-Ljungdahl pathway in Clostridium acetobutylicum enable the identification of the native tetrahydrofolate (THF) cycle and result in early induction of solvents

Here, we report the construction of a Clostridium acetobutylicum strain ATCC 824 (pCD07239) by heterologous expression of carbonyl branch genes (CD630_0723∼CD630_0729) from Clostridium difficile, aimed at installing a heterologous Wood-Ljungdahl pathway (WLP). As part of this effort, in order to validate the methyl branch of the WLP in the C. acetobutylicum, we performed ¹³C-tracing analysis on knockdown mutants of four genes responsible for the formation of 5-methyl-tetrahydrofolate (5-methyl-THF) from formate: CA_C3201, CA_C2310, CA_C2083, and CA_C0291. While C. acetobutylicum 824 (pCD07239) could not grow autotrophically, in heterotrophic fermentation, it began producing butanol at the early growth phase (OD₆₀₀ of 0.80; 0.162 g/L butanol). In contrast, solvent production in the parent strain did not begin until the early stationary phase (OD₆₀₀ of 7.40). This study offers valuable insights for future research on biobutanol production during the early growth phase.

1H NMR-based chemometric metabolomics characterization of soymilk fermented by Bacillus subtilis BSNK-5

Microbial fermentation can endow food with unique flavors, increase its nutritional value and enhance functional characteristics. Our previous research has shown that liquid fermentation of soymilk by Bacillus subtilis BSNK-5 imparted new functional properties of to the fermented product via production of nattokinase. In this study, in order to further investigate the changes in the flavor, nutritional quality and functional characteristics of soymilk during fermentation using proton nuclear magnetic resonance (¹H NMR) metabolomics to monitor the metabolite profile of BSNK-5-fermented soymilk. A total of 44 differential metabolites were identified between BSNK-5-fermented soymilk and uninoculated/unfermented soymilk, among which the levels of flavor-related substances (acetate, isovalerate and 2-methylbutyrate), nutrient-related substances (12 free amino acids), and functional substances (taurine, GABA and genistein) significantly increased after fermentation. These metabolites were closely associated with eight potential metabolic pathways. This work highlighted the significance of BSNK-5 strain in improving the nutritional quality and functional characteristics of fermented soymilk; however, the use of the strain also caused flavor deterioration. This study lays a theoretical foundation for the improvement and development of fermented soy products via liquid fermentation with B. subtilis.

Biosynthetic Incorporation of Site-Specific Isotopes in β-Lactam Antibiotics Enables Biophysical Studies

A biophysical understanding of the mechanistic, chemical, and physical origins underlying antibiotic action and resistance is vital to the discovery of novel therapeutics and the development of strategies to combat the growing emergence of antibiotic resistance. The site-specific introduction of stable-isotope labels into chemically complex natural products is particularly important for techniques such as NMR, IR, mass spectrometry, imaging, and kinetic isotope effects. Toward this goal, we developed a biosynthetic strategy for the site-specific incorporation of ¹³C labels into the canonical β-lactam carbonyl of penicillin G and cefotaxime, the latter via cephalosporin C. This was achieved through sulfur-replacement with 1-¹³C-l-cysteine, resulting in high isotope incorporations and milligram-scale yields. Using ¹³C NMR and isotope-edited IR difference spectroscopy, we illustrate how these molecules can be used to interrogate interactions with their protein targets, e.g., TEM-1 β-lactamase. This method provides a feasible route to isotopically labeled penicillin and cephalosporin precursors for future biophysical studies.

Flux prediction using artificial neural network (ANN) for the upper part of glycolysis

The selection of optimal enzyme concentration in multienzyme cascade reactions for the highest product yield in practice is very expensive and time-consuming process. The modelling of biological pathways is a difficult process because of the complexity of the system. The mathematical modelling of the system using an analytical approach depends on the many parameters of enzymes which rely on tedious and expensive experiments. The artificial neural network (ANN) method has been successively applied in different fields of science to perform complex functions. In this study, ANN models were trained to predict the flux for the upper part of glycolysis as inferred by NADH consumption, using four enzyme concentrations i.e., phosphoglucoisomerase, phosphofructokinase, fructose-bisphosphate-aldolase, triose-phosphate-isomerase. Out of three ANN algorithms, the neuralnet package with two activation functions, “logistic” and “tanh” were implemented. The prediction of the flux was very efficient: RMSE and R² were 0.847, 0.93 and 0.804, 0.94 respectively for logistic and tanh functions using a cross validation procedure. This study showed that a systemic approach such as ANN could be used for accurate prediction of the flux through the metabolic pathway. This could help to save a lot of time and costs, particularly from an industrial perspective. The R-code is available at: https://github.com/DSIMB/ANN-Glycolysis-Flux-Prediction.

13C-assisted metabolic flux analysis to investigate heterotrophic and mixotrophic metabolism in Cupriavidus necator H16

INTRODUCTION

Cupriavidus necator H16 is a gram-negative bacterium, capable of lithoautotrophic growth by utilizing hydrogen as an energy source and fixing carbon dioxide (CO2) through Calvin-Benson-Bassham (CBB) cycle. The potential to utilize synthesis gas (Syngas) and the prospects of rerouting carbon from polyhydroxybutyrate synthesis to value-added compounds makes C. necator an excellent chassis for industrial application.

OBJECTIVES

In the context of lack of sufficient quantitative information of the metabolic pathways and to advance in rational metabolic engineering for optimized product synthesis in C. necator H16, we carried out a metabolic flux analysis based on steady-state 13C-labelling.

METHODS

In this study, steady-state carbon labelling experiments, using either d-[1-13C]fructose or [1,2-13C]glycerol, were undertaken to investigate the carbon flux through the central carbon metabolism in C. necator H16 under heterotrophic and mixotrophic growth conditions, respectively.

RESULTS

We found that the CBB cycle is active even under heterotrophic condition, and growth is indeed mixotrophic. While Entner-Doudoroff (ED) pathway is shown to be the major route for sugar degradation, tricarboxylic acid (TCA) cycle is highly active in mixotrophic condition. Enhanced flux is observed in reductive pentose phosphate pathway (redPPP) under the mixotrophic condition to supplement the precursor requirement for CBB cycle. The flux distribution was compared to the mRNA abundance of genes encoding enzymes involved in key enzymatic reactions of the central carbon metabolism.

CONCLUSION

This study leads the way to establishing 13C-based quantitative fluxomics for rational pathway engineering in C. necator H16.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.

Computing the various pathways of penicillin synthesis and their molar yields

More than 80 years after its discovery, penicillin is still a widely used and commercially highly important antibiotic. Here, we analyse the metabolic network of penicillin synthesis in Penicillium chrysogenum based on the concept of elementary flux modes. In particular, we consider the synthesis of the invariant molecular core of the various subtypes of penicillin and the two major ways of incorporating sulfur: transsulfuration and direct sulfhydrylation. 66 elementary modes producing this invariant core are obtained. These show four different yields with respect to glucose, notably ½, 2/5, 1/3, and 2/7, with the highest yield of ½ occurring only when direct sulfhydrylation is used and α-aminoadipate is completely recycled. In the case of no recycling of this intermediate, we find the maximum yield to be 2/7. We compare these values with earlier literature values. Our analysis provides a systematic overview of the redundancy in penicillin synthesis and a detailed insight into the corresponding routes. Moreover, we derive suggestions for potential knockouts that could increase the average yield.

Metabolic flux analysis of Arthrobacter sp. CGMCC 3584 for cAMP production based on 13C tracer experiments and gas chromatography-mass spectrometry

Arthrobacter sp. CGMCC 3584 are able to produce cAMP from glucose by the purine synthesis pathway via de novo or salvage biosynthesis. In order to gain an improved understanding of its metabolism, (13)C-labeling experiment and gas chromatography-mass spectrometry (GC-MS) analysis were employed to determine the metabolic network structure and estimate the intracellular fluxes. GC-MS analysis helps to reflect the activity of the intracellular pathways and reactions. The metabolic network mainly contains glycolytic and pentose phosphate pathways, the tricarboxylic acid cycle, and the inactive glyoxylate shunt. Hypoxanthine as a precursor of cAMP and sodium fluoride as an inhibitor of glycolysis were found to increase the cAMP production, as well as the flux through the PP pathway. The effects of adding hypoxanthine and sodium fluoride are discussed based on the enzyme assays and metabolic flux analysis. In conclusion, our results provide quantitative insights into how cells manipulate the metabolic network under different culture conditions and this may be of value in metabolic regulation for desirable production.

Quantitative quenching evaluation and direct intracellular metabolite analysis in Penicillium chrysogenum

Sustained progress in metabolic engineering methodologies has stimulated new efforts toward optimizing fungal production strains such as through metabolite analysis of Penicillium chrysogenum industrial-scale processes. Accurate intracellular metabolite quantification requires sampling procedures that rapidly stop metabolism (quenching) and avoid metabolite loss via the cell membrane (leakage). When sampling protocols are validated, the quenching efficiency is generally not quantitatively assessed. For fungal metabolomics, quantitative biomass separation using centrifugation is a further challenge. In this study, P. chrysogenum intracellular metabolites were quantified directly from biomass extracts using automated sampling and fast filtration. A master/slave bioreactor concept was applied to provide industrial production conditions. Metabolic activity during sampling was monitored by 13C tracing. Enzyme activities were efficiently stopped and metabolite leakage was absent. This work provides a reliable method for P. chrysogenum metabolomics and will be an essential base for metabolic engineering of industrial processes.

Highlighting the tricarboxylic acid cycle: liquid and gas chromatography-mass spectrometry analyses of (13)C-labeled organic acids

The tricarboxylic acid (TCA) cycle is involved in the complete oxidation of organic acids to carbon dioxide in aerobic cells. It not only uses the acetyl-CoA derived from glycolysis but also uses breakdown products of proteins, fatty acids, and nucleic acids. Therefore, the TCA cycle involves numerous carbon fluxes through central metabolism to produce reductant power and transfer the generated electrons to the aerobic electron transport system where energy is formed by oxidative phosphorylation. Although the TCA cycle plays a crucial role in aerobic organisms and tissues, the lack of direct isotopic labeling information in its intermediates (organic acids) makes the quantification of its metabolic fluxes rather approximate. This is the major technical gap that this study intended to fill. In this work, we established and validated liquid and gas chromatography-mass spectrometry methods to determine (13)C labeling in organic acids involved in the TCA cycle using scheduled multiple reaction monitoring and single ion monitoring modes, respectively. Labeled samples were generated using maize embryos cultured with [(13)C]glucose or [(13)C]glutamine. Once steady-state labeling was reached, (13)C-labeled organic acids were extracted and purified. When applying our mass spectrometric methods to those extracts, mass isotopomer abundances of seven major organic acids were successfully determined.

Substrate cycles in Penicillium chrysogenum quantified by isotopic non-stationary flux analysis

Background

Penicillium chrysogenum, the main production strain for penicillin-G, has a high content of intracellular carbohydrates, especially reduced sugars such as mannitol, arabitol, erythritol, as well as trehalose and glycogen. In previous steady state ¹³C wash-in experiments a delay of labeling enrichments in glycolytic intermediates was observed, which suggests turnover of storage carbohydrates. The turnover of storage pools consumes ATP which is expected to reduce the product yield for energy demanding production pathways like penicillin-G.

Results

In this study, a ¹³C labeling wash-in experiment of 1 hour was performed to systematically quantify the intracellular flux distribution including eight substrate cycles. The experiments were performed using a mixed carbon source of 85% Cmol_Glc/Cmol_Glc+EtOH labeled glucose (mixture of 90% [1-¹³C₁] and 10% [U-¹³C₆]) and 15% ethanol [U-¹³C₂]. It was found, that (1) also several extracellular pools are enriched with ¹³C labeling rapidly (trehalose, mannitol, and others), (2) the intra- to extracellular metabolite concentration ratios were comparable for a large set of metabolites while for some carbohydrates (mannitol, trehalose, and glucose) the measured ratios were much higher.

Conclusions

The fast enrichment of several extracellular carbohydrates and a concentration ratio higher than the ratio expected from cell lysis (2%) indicate active (e.g. ATP consuming) transport cycles over the cellular membrane. The flux estimation indicates, that substrate cycles account for about 52% of the gap in the ATP balance based on metabolic flux analysis.

In silico aided metabolic engineering of Streptomyces roseosporus for daptomycin yield improvement

In silico metabolic network models are valuable tools for strain improvement with desired properties. In this work, based on the comparisons of each pathway flux under two different objective functions for the reconstructed metabolic network of Streptomyces roseosporus, three potential targets of zwf2 (code for glucose-6-phosphate hydrogenase), dptI (code for α-ketoglutarate methyltransferase), and dptJ (code for tryptophan oxygenase) were identified and selected for the genetic modifications. Overexpression of zwf2, dptI, and dptJ genes increased the daptomycin concentration up to 473.2, 452.5, and 489.1 mg/L, respectively. Furthermore, co-overexpression of three genes in series resulted in a 34.4% higher daptomycin concentration compared with the parental strain, which ascribed to the synergistic effect of the enzymes responsible for daptomycin biosynthesis. Finally, the engineered strain enhanced the yield of daptomycin up to 581.5 mg/L in the fed-batch culture, which was approximately 43.2% higher than that of the parental strain. These results demonstrated that the metabolic network based on in silico prediction would be accurate, reasonable, and practical for target gene identification and strain improvement.

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

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

Metabolic fluxes in the central carbon metabolism of Dinoroseobacter shibae and Phaeobacter gallaeciensis, two members of the marine Roseobacter clade

BACKGROUND

In the present work the central carbon metabolism of Dinoroseobacter shibae and Phaeobacter gallaeciensis was studied at the level of metabolic fluxes. These two strains belong to the marine Roseobacter clade, a dominant bacterial group in various marine habitats, and represent surface-associated, biofilm-forming growth (P. gallaeciensis) and symbiotic growth with eukaryotic algae (D. shibae). Based on information from recently sequenced genomes, a rich repertoire of pathways has been identified in the carbon core metabolism of these organisms, but little is known about the actual contribution of the various reactions in vivo.

RESULTS

Using 13C labelling techniques in specifically designed experiments, it could be shown that glucose-grown cells of D. shibae catabolise the carbon source exclusively via the Entner-Doudoroff pathway, whereas alternative routes of glycolysis and the pentose phosphate pathway are obviously utilised for anabolic purposes only. Enzyme assays confirmed this flux pattern and link the lack of glycolytic flux to the absence of phosphofructokinase activity. The previously suggested formation of phosphoenolpyruvate from pyruvate during mixotrophic CO2 assimilation was found to be inactive under the conditions studied. Moreover, it could be shown that pyruvate carboxylase is involved in CO2 assimilation and that the cyclic respiratory mode of the TCA cycle is utilised. Interestingly, the use of intracellular pathways was highly similar for P. gallaeciensis.

CONCLUSION

The present study reveals the first insight into pathway utilisation within the Roseobacter group. Fluxes through major intracellular pathways of the central carbon metabolism, which are closely linked to the various important traits found for the Roseobacter clade, could be determined. The close similarity of fluxes between the two physiologically rather different species might provide the first indication of more general key properties among members of the Roseobacter clade which may explain their enormous success in the marine realm.

Metabolic flux analysis in plants: coping with complexity

Theory and experience in metabolic engineering both show that metabolism operates at the network level. In plants, this complexity is compounded by a high degree of compartmentation and the synthesis of a very wide array of secondary metabolic products. A further challenge to understanding and predicting plant metabolic function is posed by our ignorance about the structure of metabolic networks even in well-studied systems. Metabolic flux analysis (MFA) provides tools to measure and model the functioning of metabolism, and is making significant contributions to coping with their complexity. This review gives an overview of different MFA approaches, the measurements required to implement them and the information they yield. The application of MFA methods to plant systems is then illustrated by several examples from the recent literature. Next, the challenges that plant metabolism poses for MFA are discussed together with ways that these can be addressed. Lastly, new developments in MFA are described that can be expected to improve the range and reliability of plant MFA in the coming years.

Studies of the production of fungal polyketides in Aspergillus nidulans by using systems biology tools

Many filamentous fungi produce polyketide molecules with great significance as human pharmaceuticals; these molecules include the cholesterol-lowering compound lovastatin, which was originally isolated from Aspergillus terreus. The chemical diversity and potential uses of these compounds are virtually unlimited, and it is thus of great interest to develop a well-described microbial production platform for polyketides. Using genetic engineering tools available for the model organism Aspergillus nidulans, we constructed two recombinant strains, one expressing the Penicillium griseofulvum 6-methylsalicylic acid (6-MSA) synthase gene and one expressing the 6-MSA synthase gene and overexpressing the native xylulose-5-phosphate phosphoketolase gene (xpkA) for increasing the pool of polyketide precursor levels. The physiology of the recombinant strains and that of a reference wild-type strain were characterized on glucose, xylose, glycerol, and ethanol media in controlled bioreactors. Glucose was found to be the preferred carbon source for 6-MSA production, and 6-MSA concentrations up to 455 mg/liter were obtained for the recombinant strain harboring the 6-MSA gene. Our findings indicate that overexpression of xpkA does not directly improve 6-MSA production on glucose, but it is possible, if the metabolic flux through the lower part of glycolysis is reduced, to obtain quite high yields for conversion of sugar to 6-MSA. Systems biology tools were employed for in-depth analysis of the metabolic processes. Transcriptome analysis of 6-MSA-producing strains grown on glucose and xylose in the presence and absence of xpkA overexpression, combined with flux and physiology data, enabled us to propose an xpkA-msaS interaction model describing the competition between biomass formation and 6-MSA production for the available acetyl coenzyme A.

Systems analysis unfolds the relationship between the phosphoketolase pathway and growth in Aspergillus nidulans

BACKGROUND

Aspergillus nidulans is an important model organism for studies on fundamental eukaryotic cell biology and on industrial processes due to its close relation to A. niger and A. oryzae. Here we identified the gene coding for a novel metabolic pathway in A. nidulans, namely the phosphoketolase pathway, and investigated the role of an increased phosphoketolase activity.

METHODOLOGY/PRINCIPAL FINDINGS

Over-expression of the phosphoketolase gene (phk) improved the specific growth rate on xylose, glycerol and ethanol. Transcriptome analysis showed that a total of 1,222 genes were significantly affected by over-expression of the phk, while more than half of the affected genes were carbon source specific. During growth on glucose medium, the transcriptome analysis showed that the response to phk over-expression is targeted to neutralize the effect of the over-expression by regulating the acetate metabolism and initiate a growth dampening response.

CONCLUSIONS/SIGNIFICANCE

Metabolic flux analysis using (13)C-labelled glucose, showed that over-expression of phosphoketolase added flexibility to the central metabolism. Our findings further suggests that A. nidulans is not optimized for growth on xylose, glycerol or ethanol as the sole carbon sources.

Overexpression of a novel endogenous NADH kinase in Aspergillus nidulans enhances growth

The complete genome sequence of the filamentous fungi Aspergillus nidulans has paved the way for fundamental research on this industrially important species. To the best of our knowledge, this is the first time a gene encoding for ATP-dependent NADH kinase (ATP:NADH 2'-phosphotransferase, EC 2.7.1.86) has been identified. The enzyme has a predicted molecular weight of 49 kDa. We characterised the role of this NADH kinase by genomic integration of the putative gene AN8837.2 under a strong constitutive promoter. The physiological effects of overexpressed NADH kinase in combination with different aeration rates were studied in well-controlled glucose batch fermentations. Metabolite profiling and metabolic network analysis with [1-(13)C] glucose were used for characterisation of the strains, and the results demonstrated that NADH kinase activity has paramount influence on growth physiology. Biomass yield on glucose and the maximum specific growth rate increased from 0.47 g/g and 0.22 h(-1) (wild type) to 0.54 g/g and 0.26 h(-1) (NADH kinase overexpressed), respectively. The results suggest that overexpression of NADH kinase improves the growth efficiency of the cell by increasing the access to NADPH. Our findings indicate that A. nidulans is not optimised for growth in nutrient-rich conditions typically found in laboratory and industrial fermentors. This conclusion may impact the design of new strains capable of generating reducing power in the form of NADPH, which is crucial for efficient production of many industrially important metabolites and enzymes.

Elucidation of an alternate isoleucine biosynthesis pathway in Geobacter sulfurreducens

The central metabolic model for Geobacter sulfurreducens included a single pathway for the biosynthesis of isoleucine that was analogous to that of Escherichia coli, in which the isoleucine precursor 2-oxobutanoate is generated from threonine. 13C labeling studies performed in G. sulfurreducens indicated that this pathway accounted for a minor fraction of isoleucine biosynthesis and that the majority of isoleucine was instead derived from acetyl-coenzyme A and pyruvate, possibly via the citramalate pathway. Genes encoding citramalate synthase (GSU1798), which catalyzes the first dedicated step in the citramalate pathway, and threonine ammonia-lyase (GSU0486), which catalyzes the conversion of threonine to 2-oxobutanoate, were identified and knocked out. Mutants lacking both of these enzymes were auxotrophs for isoleucine, whereas single mutants were capable of growth in the absence of isoleucine. Biochemical characterization of the single mutants revealed deficiencies in citramalate synthase and threonine ammonia-lyase activity. Thus, in G. sulfurreducens, 2-oxobutanoate can be synthesized either from citramalate or threonine, with the former being the main pathway for isoleucine biosynthesis. The citramalate synthase of G. sulfurreducens constitutes the first characterized member of a phylogenetically distinct clade of citramalate synthases, which contains representatives from a wide variety of microorganisms.

Production of L-carnitine by secondary metabolism of bacteria

The increasing commercial demand for L-carnitine has led to a multiplication of efforts to improve its production with bacteria. The use of different cell environments, such as growing, resting, permeabilized, dried, osmotically stressed, freely suspended and immobilized cells, to maintain enzymes sufficiently active for L-carnitine production is discussed in the text. The different cell states of enterobacteria, such as Escherichia coli and Proteus sp., which can be used to produce L-carnitine from crotonobetaine or D-carnitine as substrate, are analyzed. Moreover, the combined application of both bioprocess and metabolic engineering has allowed a deeper understanding of the main factors controlling the production process, such as energy depletion and the alteration of the acetyl-CoA/CoA ratio which are coupled to the end of the biotransformation. Furthermore, the profiles of key central metabolic activities such as the TCA cycle, the glyoxylate shunt and the acetate metabolism are seen to be closely interrelated and affect the biotransformation efficiency. Although genetically modified strains have been obtained, new strain improvement strategies are still needed, especially in Escherichia coli as a model organism for molecular biology studies. This review aims to summarize and update the state of the art in L-carnitine production using E. coli and Proteus sp, emphasizing the importance of proper reactor design and operation strategies, together with metabolic engineering aspects and the need for feed-back between wet and in silico work to optimize this biotransformation.

Science – or not? The status and dynamics of biotechnology

This review traces the emergence of biotechnology as a new scientific discipline since the 1980s, when it became a major economic force. Significant changes in theoretical perception, research strategies, aims, and experimental methods, mainly in genetic engineering techniques, occurred during this period. The article is based on an analysis of its scientific status over four decades: the 60s and 70s when work in the field proceeded in different disciplines with a low level of coherence and little integration, then a significant change during the 80s and 90s when common approaches and the merging of molecular biology and biochemical engineering created a new discipline. The analysis covers scientific highlights and outstanding technical progress, presenting two studies undertaken by scientific and governmental agencies in Germany and the USA, as well as results of interviews and a questionnaire dealing with the scientific status of biotechnology. Answers to the questionnaire were obtained from internationally known scientists and from young scientists with biotechnology degrees. The results collected trace the transition of biotechnology from heterogeneous specialties and approaches towards a scientific discipline of its own. A hypothesis is put forward suggesting a new common paradigm allowing for a coherent perception the of phenomena observed.

A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.

Engineering primary metabolic pathways of industrial micro-organisms

Metabolic engineering is a powerful tool for the optimisation and the introduction of new cellular processes. This is mostly done by genetic engineering. Since the introduction of this multidisciplinary approach, the success stories keep accumulating. The primary metabolism of industrial micro-organisms has been studied for long time and most biochemical pathways and reaction networks have been elucidated. This large pool of biochemical information, together with data from proteomics, metabolomics and genomics underpins the strategies for design of experiments and choice of targets for manipulation by metabolic engineers. These targets are often located in the primary metabolic pathways, such as glycolysis, pentose phosphate pathway, the TCA cycle and amino acid biosynthesis and mostly at major branch points within these pathways. This paper describes approaches taken for metabolic engineering of these pathways in bacteria, yeast and filamentous fungi.

Improving metabolic flux estimation via evolutionary optimization for convex solution space

MOTIVATION

Flux estimation by using (13) C-labeling pattern information of metabolites is currently the only method that can give accurate, detailed quantification of all intracellular fluxes in the central metabolism of a microorganism. In essence, it corresponds to a constrained optimization problem which minimizes a weighted distance between measured and simulated results. Characteristics, such as existence of multiple local minima, non-linear and non-differentiable make this problem a special difficulty.

RESULTS

In the present work, we propose an evolutionary-based global optimization algorithm taking advantage of the convex feature of the problem's solution space. Based on the characteristics of convex spaces, specialized initial population and evolutionary operators are designed to solve (13)C-based metabolic flux estimation problem robustly and efficiently. The algorithm was applied to estimate the central metabolic fluxes in Escherichia coli and compared with conventional optimization technique. Experimental results illustrated that our algorithm is capable of achieving fast convergence to good near-optima and maintaining the robust nature of evolutionary algorithms at the same time.

AVAILABILITY

Available from the authors upon request.

Determining Actinobacillus succinogenes metabolic pathways and fluxes by NMR and GC-MS analyses of 13C-labeled metabolic product isotopomers

Actinobacillus succinogenes is a promising candidate for industrial succinate production. However, in addition to producing succinate, it also produces formate and acetate. To understand carbon flux distribution to succinate and alternative products we fed A. succinogenes [1-(13)C]glucose and analyzed the resulting isotopomers of excreted organic acids, proteinaceous amino acids, and glycogen monomers by gas chromatography-mass spectrometry and nuclear magnetic resonance spectroscopy. The isotopomer data, together with the glucose consumption and product formation rates and the A. succinogenes biomass composition, were supplied to a metabolic flux model. Oxidative pentose phosphate pathway flux supplied, at most, 20% of the estimated NADPH requirement for cell growth. The model indicated that NADPH was instead produced primarily by the conversion of NADH to NADPH by transhydrogenase and/or by NADP-dependent malic enzyme. Transhydrogenase activity was detected in A. succinogenes cell extracts, as were formate and pyruvate dehydrogenases, which the model suggested were contributing to NADH production. Malic enzyme activity was also detected in cell extracts, consistent with the flux analysis results. Labeling patterns in amino acids and organic acids showed that oxaloacetate and malate were being decarboxylated to pyruvate. These are the first in vivo experiments to show that the partitioning of flux between succinate and alternative fermentation products can occur at multiple nodes in A. succinogenes. The implications for designing effective metabolic engineering strategies to increase A. succinogenes succinate production are discussed.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry for metabolic flux analyses using isotope-labeled ethanol

We describe a novel method for the determination of the concentration and labeling degree of ethanol originating from 1-13C-labeling experiments. This method is suitable for high-throughput metabolic flux analysis because of the possible parallel sample preparation and fast final analysis using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS). In a closed vial containing culture supernatant, ethanol is enzymatically oxidized to acetaldehyde. The acetaldehyde formed evaporates and is readily trapped in a second enclosed but open vial containing acidified 2,4-dinitrophenylhydrazine (DNPH). The 2,4-acetaldehyde dinitrophenylhydrazone (Ac-DNPH) that is formed is insoluble under these conditions. This leads to a constant conversion rate of the acetaldehyde produced from ethanol after 14 h minimum incubation time. MALDI-TOFMS was used to quantify the formed Ac-DNPH with [13C2]-ethanol as internal standard. The relative signal intensities of the unlabeled ethanol derivative as well as of [1-13C]-ethanol were linearly related to the ethanol concentration within a range of 1 to 50 mM with a limit of detection of 0.6 mM, a range which is sufficient for flux analysis in microtiter plate fermentation experiments. The method allows the estimation of the [1-13C]-ethanol originating from 1-13C-labeling experiments of Saccharomyces cerevisiae strains. In experiments where the expected flux range was exceeded, unlabeled ethanol was determined with a linear range from 30 to 500 mM. Ethanol quantification using this method was compared with enzymatic analysis and exhibited differences of less than 3.3% on average. Comparison of flux partitioning ratios between glycolysis and the pentose-phosphate pathway (PPP) based on MALDI-TOFMS and gas chromatography (GC)/MS methods showed good agreement, with differences for ethanol and alanine labeling of only 4.3%.

Isotopomer analysis of myocardial substrate metabolism: a systems biology approach

The increasing accessibility of mass isotopomer data via GC-MS and NMR technology has necessitated the use of a systematic and reliable method to take advantage of such data for flux analysis. Here we applied a nonlinear, optimization-based method to study substrate metabolism in cardiomyocytes using (13)C data from perfused mouse hearts. The myocardial metabolic network used in this study accounts for 257 reactions and 240 metabolites, which are further compartmentalized into extracellular space, cytosol, and mitochondrial matrix. Analysis of the perfused mouse heart showed that the steady-state ATP production rate was 16.6 +/- 2.3 micromol/min . gww, with 30% of the ATP coming from glycolysis. Of the four substrates available in the perfusate (glucose, pyruvate, lactate, and oleate), exogenous glucose forms the majority of cytosolic pyruvate. Pyruvate decaboxylation is significantly higher than carboxylation, suggesting that anaplerosis is low in the perfused heart. Exchange fluxes were predicted to be high for reversible enzymes in the citric acid cycle (CAC), but low in the glycolytic pathway. Pseudoketogenesis amounted to approximately 50% of the net ketone body uptake. Sensitivity analysis showed that the estimated flux distributions were relatively insensitive to experimental errors. The application of isotopomer data drastically improved the estimation of reaction fluxes compared to results computed with respect to reaction stoichiometry alone. Further study of 12 commonly used (13)C glucose mixtures showed that the mixtures of 20% [U-(13)C(6)] glucose, 80% [3 (13)C] glucose and 20% [U-(13)C(6)] glucose, 80% [4 (13)C] were best for resolving fluxes in the current network.

13C-labeled gluconate tracing as a direct and accurate method for determining the pentose phosphate pathway split ratio in Penicillium chrysogenum

In this study we developed a new method for accurately determining the pentose phosphate pathway (PPP) split ratio, an important metabolic parameter in the primary metabolism of a cell. This method is based on simultaneous feeding of unlabeled glucose and trace amounts of [U-13C]gluconate, followed by measurement of the mass isotopomers of the intracellular metabolites surrounding the 6-phosphogluconate node. The gluconate tracer method was used with a penicillin G-producing chemostat culture of the filamentous fungus Penicillium chrysogenum. For comparison, a 13C-labeling-based metabolic flux analysis (MFA) was performed for glycolysis and the PPP of P. chrysogenum. For the first time mass isotopomer measurements of 13C-labeled primary metabolites are reported for P. chrysogenum and used for a 13C-based MFA. Estimation of the PPP split ratio of P. chrysogenum at a growth rate of 0.02 h(-1) yielded comparable values for the gluconate tracer method and the 13C-based MFA method, 51.8% and 51.1%, respectively. A sensitivity analysis of the estimated PPP split ratios showed that the 95% confidence interval was almost threefold smaller for the gluconate tracer method than for the 13C-based MFA method (40.0 to 63.5% and 46.0 to 56.5%, respectively). From these results we concluded that the gluconate tracer method permits accurate determination of the PPP split ratio but provides no information about the remaining cellular metabolism, while the 13C-based MFA method permits estimation of multiple fluxes but provides a less accurate estimate of the PPP split ratio.

Engineering of the redox imbalance of Fusarium oxysporum enables anaerobic growth on xylose

Dissimilatory nitrate reduction metabolism, of the natural xylose-fermenting fungus Fusarium oxysporum, was used as a strategy to achieve anaerobic growth and ethanol production from xylose. Beneficial alterations of the redox fluxes and thereby of the xylose metabolism were obtained by taking advantage of the regeneration of the cofactor NAD(+) during the denitrification process. In batch cultivations, nitrate sustained growth under anaerobic conditions (1.21 g L(-1) biomass) and simultaneously a maximum yield of 0.55 moles of ethanol per mole of xylose was achieved, whereas substitution of nitrate with ammonium limited the growth significantly (0.15 g L(-1) biomass). Using nitrate, the maximum acetate yield was 0.21 moles per mole of xylose and no xylitol excretion was observed. Furthermore, the network structure in the central carbon metabolism of F. oxysporum was characterized in steady state. F. oxysporum grew anaerobically on [1-(13)C] labelled glucose and unlabelled xylose in chemostat cultivation with nitrate as nitrogen source. The use of labelled substrate allowed the precise determination of the glucose and xylose contribution to the carbon fluxes in the central metabolism of this poorly described microorganism. It was demonstrated that dissimilatory nitrate reduction allows F. oxysporum to exhibit typical respiratory metabolic behaviour with a highly active TCA cycle and a large demand for NADPH.

MET-IDEA: Data Extraction Tool for Mass Spectrometry-Based Metabolomics

A current and significant limitation to metabolomics is the large-scale, high-throughput conversion of raw chromatographically coupled mass spectrometry datasets into organized data matrices necessary for further statistical processing and data visualization. This article describes a new data extraction tool, MET-IDEA (Metabolomics Ion-based Data Extraction Algorithm) which surmounts this void. MET-IDEA is compatible with a diversity of chromatographically coupled mass spectrometry systems, generates an output similar to traditional quantification methods, utilizes the sensitivity and selectivity associated with selected ion quantification, and greatly reduces the time and effort necessary to obtain large-scale organized datasets by several orders of magnitude. The functionality of MET-IDEA is illustrated using metabolomics data obtained for elicited cell culture exudates from the model legume, Medicago truncatula. The results indicate that MET-IDEA is capable of rapidly extracting semiquantitative data from raw data files, which allows for more rapid biological insight. MET-IDEA is freely available to academic users upon request.

Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements

Metabolic fluxes, estimated from stable isotope studies, provide a key to quantifying physiology in fields ranging from metabolic engineering to the analysis of human metabolic diseases. A serious drawback of the flux estimation method in current use is that it does not produce confidence limits for the estimated fluxes. Without this information it is difficult to interpret flux results and expand the physiological significance of flux studies. To address this shortcoming we derived analytical expressions of flux sensitivities with respect to isotope measurements and measurement errors. These tools allow the determination of local statistical properties of fluxes and relative importance of measurements. Furthermore, we developed an efficient algorithm to determine accurate flux confidence intervals and demonstrated that confidence intervals obtained with this method closely approximate true flux uncertainty. In contrast, confidence intervals approximated from local estimates of standard deviations are inappropriate due to inherent system nonlinearities. We applied these methods to analyze the statistical significance and confidence of estimated gluconeogenesis fluxes from human studies with [U-13C]glucose as tracer and found true limits for flux estimation in specific human isotopic protocols.

CreA influences the metabolic fluxes of Aspergillus nidulans during growth on glucose and xylose

The physiological phenotype of Aspergillus nidulans was investigated for different genetic and environmental conditions of glucose repression through the quantification of in vivo fluxes in the central carbon metabolism using (13)C-metabolic-flux analysis. The particular focus was the role of the carbon repressor CreA, which is the major regulatory protein mediating carbon repression in many fungal species, in the primary metabolism of A. nidulans. Batch cultivations were performed with a reference strain and a deletion mutant strain (creADelta4) using [1-(13)C]glucose as carbon source. The mutant strain was also grown on a mixture of [1-(13)C]glucose and unlabelled xylose. Fractional enrichment data were measured by gas chromatography-mass spectrometry. A model describing the central metabolism of A. nidulans was used in combination with fractional enrichment data, and measurements of extracellular rates and biomass composition for the estimation of the in vivo metabolic fluxes. The creA-mutant strain showed a lower maximum specific growth rate than the reference strain when grown on glucose (0.11 and 0.25 h(-1), respectively), whereas the specific growth rate of the mutant strain grown on the glucose/xylose mixture was identical to that on glucose (0.11 h(-1)). Different patterns and increased levels of extracellular polyols were observed both upon deletion of the creA gene and upon addition of xylose to the growth medium of the mutant strain. Concerning metabolic fluxes, the major change observed in the flux distribution of A. nidulans upon deletion of the creA gene was a 20 % decrease in the flux through the oxidative part of the pentose-phosphate pathway. Addition of xylose to the growth medium of the mutant resulted in an increase of about 40 % in the activity of the oxidative part of the pentose-phosphate pathway, as well as decreases in the fluxes through the Embden-Meyerhof-Parnas pathway and the tricarboxylic acid cycle (in the range of 20-30 %). The derepression of key pathways leads to alterations in the demands for cofactors, thereby imposing changes in the central metabolism due to the coupling of the many different reactions via the redox and energy metabolism of the cells.

Revisiting the 13C-label distribution of the non-oxidative branch of the pentose phosphate pathway based upon kinetic and genetic evidence

The currently applied reaction structure in stoichiometric flux balance models for the nonoxidative branch of the pentose phosphate pathway is not in accordance with the established ping-pong kinetic mechanism of the enzymes transketolase (EC 2.2.1.1) and transaldolase (EC 2.2.1.2). Based upon the ping-pong mechanism, the traditional reactions of the nonoxidative branch of the pentose phosphate pathway are replaced by metabolite specific, reversible, glycolaldehyde moiety (C(2)) and dihydroxyacetone moiety (C(3)) fragments producing and consuming half-reactions. It is shown that a stoichiometric model based upon these half-reactions is fundamentally different from the currently applied stoichiometric models with respect to the number of independent C(2) and C(3) fragment pools in the pentose phosphate pathway and can lead to different label distributions for (13)C-tracer experiments. To investigate the actual impact of the new reaction structure on the estimated flux patterns within a cell, mass isotopomer measurements from a previously published (13)C-based metabolic flux analysis of Saccharomyces cerevisiae were used. Different flux patterns were found. From a genetic point of view, it is well known that several micro-organisms, including Escherichia coli and S. cerevisiae, contain multiple genes encoding isoenzymes of transketolase and transaldolase. However, the extent to which these gene products are also actively expressed remains unknown. It is shown that the newly proposed stoichiometric model allows study of the effect of isoenzymes on the (13)C-label distribution in the nonoxidative branch of the pentose phosphate pathway by extending the half-reaction based stoichiometric model with two distinct transketolase enzymes instead of one. Results show that the inclusion of isoenzymes affects the ensuing flux estimates.

Comparative metabolic network analysis of two xylose fermenting recombinant Saccharomyces cerevisiae strains

The recombinant xylose fermenting strain Saccharomyces cerevisiae TMB3001 can grow on xylose, but the xylose utilisation rate is low. One important reason for the inefficient fermentation of xylose to ethanol is believed to be the imbalance of redox co-factors. In the present study, a metabolic flux model was constructed for two recombinant S. cerevisiae strains: TMB3001 and CPB.CR4 which in addition to xylose metabolism have a modulated redox metabolism, i.e. ammonia assimilation was shifted from being NADPH to NADH dependent by deletion of gdh1 and over-expression of GDH2. The intracellular fluxes were estimated for both strains in anaerobic continuous cultivations when the growth limiting feed consisted of glucose (2.5 g L-1) and xylose (13 g L-1). The metabolic network analysis with 13C labelled glucose showed that there was a shift in the specific xylose reductase activity towards use of NADH as co-factor rather than NADPH. This shift is beneficial for solving the redox imbalance and it can therefore partly explain the 25% increase in the ethanol yield observed for CPB.CR4. Furthermore, the analysis indicated that the glyoxylate cycle was activated in CPB.CR4.

High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints

13C-constrained flux balancing analysis based on gas chromatography-mass spectrometry data is presented here as a simple and robust method for the estimation of intracellular carbon fluxes. In this approach, the underdetermined system of metabolite balances deduced from stoichiometric relations and measured extracellular rates is complemented with 13C constraints from metabolic flux ratio analysis. Fluxes in central carbon metabolism of exponentially growing Escherichia coli were estimated by 13C-constrained flux balancing from three different 13C-labeled glucose experiments. The best resolution of the network was achieved using 13C constraints derived from [U-13C]glucose and [1-13C]glucose experiments. The corresponding flux estimate was in excellent agreement with a solution that was independently obtained with a comprehensive isotopomer model. This new methodology was also demonstrated to faithfully capture the intracellular flux distribution in E. coli shake flasks and 1-ml deep-well microtiter plates. Due to its simplicity, speed, and robustness, 13C-constrained metabolic flux balancing is promising for routine and high-throughput analysis on a miniaturized scale.

Metabolic network analysis of Streptomyces tenebrarius, a Streptomyces species with an active entner-doudoroff pathway

Streptomyces tenebrarius is an industrially important microorganism, producing an antibiotic complex that mainly consists of the aminoglycosides apramycin, tobramycin carbamate, and kanamycin B carbamate. When S. tenebrarius is used for industrial tobramycin production, kanamycin B carbamate is an unwanted by-product. The two compounds differ only by one hydroxyl group, which is present in kanamycin carbamate but is reduced during biosynthesis of tobramycin. (13)C metabolic flux analysis was used for elucidating connections between the primary carbon metabolism and the composition of the antibiotic complex. Metabolic flux maps were constructed for the cells grown on minimal medium with glucose or with a glucose-glycerol mixture as the carbon source. The addition of glycerol, which is more reduced than glucose, led to a three-times-greater reduction of the kanamycin portion of the antibiotic complex. The labeling indicated an active Entner-Doudoroff (ED) pathway, which was previously considered to be nonfunctional in Streptomyces. The activity of the pentose phosphate (PP) pathway was low (10 to 20% of the glucose uptake rate). The fluxes through Embden-Meyerhof-Parnas (EMP) and ED pathways were almost evenly distributed during the exponential growth on glucose. During the transition from growth phase to production phase, a metabolic shift was observed, characterized by a decreased flux through the ED pathway and increased fluxes through the EMP and PP pathways. Higher specific NADH and NADPH production rates were calculated in the cultivation on glucose-glycerol, which was associated with a lower percentage of nonreduced antibiotic kanamycin B carbamate.

Glucose metabolism in the antibiotic producing actinomycete Nonomuraea sp. ATCC 39727

The actinomycete Nonomuraea sp. ATCC 39727, producer of the glycopeptide A40926 that is used as precursor for the novel antibiotic dalbavancin, has an unusual carbon metabolism. Glucose is primarily metabolized via the Entner-Doudoroff (ED) pathway, although the energetically more favorable Embden-Meyerhof-Parnas (EMP) pathway is present in this organism. Moreover, Nonomuraea utilizes a PPi-dependent phosphofructokinase, an enzyme that has been connected with anaerobic metabolism in eukaryotes and higher plants, but recently has been recognized in several actinomycetes. In order to study its primary carbon metabolism in further detail, Nonomuraea was cultivated with [1-13C] glucose as the only carbon source and the 13C-labeling patterns of proteinogenic amino acids were determined by GC-MS analysis. Through this method, the fluxes in the central carbon metabolism during balanced growth were estimated. Moreover, a shift in the label incorporation pattern was observed in connection with phosphate limitation and increased antibiotic productivity in Nonomuraea. The shift indicated an increased flux through the EMP pathway at the expense of the flux through the ED pathway, a suggestion that was supported by alterations in intracellular metabolite levels during phosphate limitation. In contrast, expression levels of genes encoding enzymes in the ED and EMP pathways were not affected by phosphate limitation.

Metabolic network analysis on Phaffia rhodozyma yeast using 13C-labeled glucose and gas chromatography-mass spectrometry

Carotenoid production by microorganisms, as opposed to chemical synthesis, could fulfill an ever-increasing demand for 'all natural' products. The yeast Phaffia rhodozyma has received considerable attention because it produces the red pigment astaxanthin, commonly used as an animal feed supplement. In order to have a better understanding of its metabolism, labeling experiments with [1-(13)C]glucose were conducted with the wildtype strain (CBS5905T) and a hyper-producing carotenoid strain (J4-3) in order to determine their metabolic network structure and estimate intracellular fluxes. Amino acid labeling patterns, as determined by GC-MS, were in accordance with a metabolic network consisting of the Embden-Meyerhof-Parnas pathway, the pentose phosphate pathway, and the TCA cycle. Glucose was mainly consumed along the pentose phosphate pathway ( approximately 65% for wildtype strain), which reflected high NADPH requirements for lipid biosynthesis. Although common to other oleaginous yeast, there was no, or very little, malic enzyme activity for carbon-limited growth. In addition, there was no evidence of phosphoketolase activity. The central carbon metabolism of the mutant strain was similar to that of the wildtype strain, though the relative pentose phosphate flux was lower and the TCA cycle flux in accordance with the biomass yield being lower.

Manipulation of malic enzyme in Saccharomyces cerevisiae for increasing NADPH production capacity aerobically in different cellular compartments

The yeast Saccharomyces cerevisiae is an attractive cell factory, but in many cases there are constraints related with balancing the formation and consumption of redox cofactors. In this work, we studied the effect of having an additional source of NADPH in the cell. In order to do this, two strains were engineered by overexpression of malic enzyme. In one of them, malic enzyme was overexpressed as its wild-type mitochondrial form, and in the other strain a short form lacking the mitochondrial targeting sequence was overexpressed. The recombinant strains were analyzed in aerobic batch and continuous cultivations, and the basic growth characteristics were generally not affected to a great extent, even though pleiotropic effects of the manipulations could be seen by the altered in vitro activities of selected enzymes of the central metabolism. Moreover, the decreased pentose-phosphate pathway flux and the ratios of redox cofactors showed that a net transhydrogenase effect was obtained, which can be directed to the cytosol or the mitochondria. This may find application in redirecting fluxes for improving specific biotechnological applications.