Elucidation of novel biosynthetic pathways and metabolite flux patterns by retrobiosynthetic NMR analysis

Metabolic flux analysis: a comprehensive review on sample preparation, analytical techniques, data analysis, computational modelling, and main application areas

Metabolic flux analysis (MFA) quantitatively describes cellular fluxes to understand metabolic phenotypes and functional behaviour after environmental and/or genetic perturbations. In the last decade, the application of stable isotopes became extremely important to determine and integrate in vivo measurements of metabolic reactions in systems biology. 13C-MFA is one of the most informative methods used to study central metabolism of biological systems. This review aims to outline the current experimental procedure adopted in 13C-MFA, starting from the preparation of cell cultures and labelled tracers to the quenching and extraction of metabolites and their subsequent analysis performed with very powerful software. Here, the limitations and advantages of nuclear magnetic resonance spectroscopy and mass spectrometry techniques used in carbon labelled experiments are elucidated by reviewing the most recent published papers. Furthermore, we summarise the most successful approaches used for computational modelling in flux analysis and the main application areas with a particular focus in metabolic engineering.

Total Synthesis of (-)-FD-838 and (-)-Cephalimysin A

We completed a nine-step total synthesis of (-)-FD-838 and (-)-cephalimysin A. Our synthesis features a biogenetically guided assembly of the highly oxidized spirocyclic core by Snider-type tandem epoxidations of the chiral substrate derived from an amino acid derivative. Our synthetic approach provides a general and versatile solution to access spirocyclic PKS-NRPS-based secondary fungal metabolites.

Improved riboflavin production with Ashbya gossypii from vegetable oil based on 13C metabolic network analysis with combined labeling analysis by GC/MS, LC/MS, 1D, and 2D NMR

The fungus Ashbya gossypii is an important industrial producer of riboflavin, i.e. vitamin B₂. In order to meet the constantly increasing demands for improved production processes, it appears essential to better understand the underlying metabolic pathways of the vitamin. Here, we used a highly sophisticated set-up of parallel ¹³C tracer studies with labeling analysis by GC/MS, LC/MS, 1D, and 2D NMR to resolve carbon fluxes in the overproducing strain A. gossypii B2 during growth and subsequent riboflavin production from vegetable oil as carbon source, yeast extract, and supplemented glycine. The studies provided a detailed picture of the underlying metabolism. Glycine was exclusively used as carbon-two donor of the vitamin's pyrimidine ring, which is part of its isoalloxazine ring structure, but did not contribute to the carbon-one metabolism due to the proven absence of a functional glycine cleavage system. The pools of serine and glycine were closely connected due to a highly reversible serine hydroxymethyltransferase. Transmembrane formate flux simulations revealed that the one-carbon metabolism displayed a severe bottleneck during initial riboflavin production, which was overcome in later phases of the cultivation by intrinsic formate accumulation. The transiently limiting carbon-one pool was successfully replenished by time-resolved feeding of small amounts of formate and serine, respectively. This increased the intracellular availability of glycine, serine, and formate and resulted in a final riboflavin titer increase of 45%.

Differential Off-line LC-NMR (DOLC-NMR) Metabolomics To Monitor Tyrosine-Induced Metabolome Alterations in Saccharomyces cerevisiae

A novel differential off-line LC-NMR approach (DOLC-NMR) was developed to capture and quantify nutrient-induced metabolome alterations in Saccharomyces cerevisiae. Off-line coupling of HPLC separation and ¹H NMR spectroscopy supported by automated comparative bucket analyses, followed by quantitative ¹H NMR using ERETIC 2 (electronic reference to access in vivo concentrations), has been successfully used to quantitatively record changes in the metabolome of S. cerevisiae upon intervention with the aromatic amino acid l-tyrosine. Among the 33 metabolites identified, glyceryl succinate, tyrosol acetate, tyrosol lactate, tyrosol succinate, and N-acyl-tyrosine derivatives such as N-(1-oxooctyl)-tyrosine are reported for the first time as yeast metabolites. Depending on the chain length, N-(1-oxooctyl)-, N-(1-oxodecanyl)-, N-(1-oxododecanyl)-, N-(1-oxomyristinyl)-, N-(1-oxopalmityl)-, and N-(1-oxooleoyl)-l-tyrosine imparted a kokumi taste enhancement above their recognition thresholds ranging between 145 and 1432 μmol/L (model broth). Finally, carbon module labeling (CAMOLA) and carbon bond labeling (CABOLA) experiments with ¹³C₆-glucose as the carbon source confirmed the biosynthetic pathway leading to the key metabolites; for example, the aliphatic side chain of N-(1-oxooctyl)-tyrosine could be shown to be generated via de novo fatty acid biosynthesis from four C₂-carbon modules (acetyl-CoA) originating from glucose.

Isoprenoid generating systems in plants - A handy toolbox how to assess contribution of the mevalonate and methylerythritol phosphate pathways to the biosynthetic process

Isoprenoids comprise an astonishingly diverse group of metabolites with numerous potential and actual applications in medicine, agriculture and the chemical industry. Generation of efficient platforms producing isoprenoids is a target of numerous laboratories. Such efforts are generally enhanced if the native biosynthetic routes can be identified, and if the regulatory mechanisms responsible for the biosynthesis of the compound(s) of interest can be determined. In this review a critical summary of the techniques applied to establish the contribution of the two alternative routes of isoprenoid production operating in plant cells, the mevalonate and methylerythritol pathways, with a focus on their co-operation (cross-talk) is presented. Special attention has been paid to methodological aspects of the referred studies, in order to give the reader a deeper understanding for the nuances of these powerful techniques. This review has been designed as an organized toolbox, which might offer the researchers comments useful both for project design and for interpretation of results obtained.

A raison d'être for two distinct pathways in the early steps of plant isoprenoid biosynthesis?

When compared to other organisms, plants are atypical with respect to isoprenoid biosynthesis: they utilize two distinct and separately compartmentalized pathways to build up isoprene units. The co-existence of these pathways in the cytosol and in plastids might permit the synthesis of many vital compounds, being essential for a sessile organism. While substrate exchange across membranes has been shown for a variety of plant species, lack of complementation of strong phenotypes, resulting from inactivation of either the cytosolic pathway (growth and development defects) or the plastidial pathway (pigment bleaching), seems to be surprising at first sight. Hundreds of isoprenoids have been analyzed to determine their biosynthetic origins. It can be concluded that in angiosperms, under standard growth conditions, C₂₀-phytyl moieties, C₃₀-triterpenes and C₄₀-carotenoids are made nearly exclusively within compartmentalized pathways, while mixed origins are widespread for other types of isoprenoid-derived molecules. It seems likely that this coexistence is essential for the interaction of plants with their environment. A major purpose of this review is to summarize such observations, especially within an ecological and functional context and with some emphasis on regulation. This latter aspect still requires more work and present conclusions are preliminary, although some general features seem to exist.

Tools for the tract: understanding the functionality of the gastrointestinal tract

The human gastrointestinal tract comprises a series of complex and dynamic organs ranging from the stomach to the distal colon, which harbor immense microbial assemblages that are known to be vital for human health. Until recently, most of the details concerning our gut microbiota remained obscure. Over the past several years, however, a number of crucial technological and conceptual innovations have been introduced to shed more light on the composition and functionality of human gut microbiota. Recently developed high throughput approaches, including next-generation sequencing technologies and phylogenetic microarrays targeting ribosomal RNA gene sequences, allow for comprehensive analysis of the diversity and dynamics of the gut microbiota composition. Nevertheless, most of the microbes especially in the human large intestine still remain uncultured, and the in situ functions of distinct groups of the gut microbiota are therefore largely unknown, but pivotal to the understanding of their role in human physiology. Apart from functional and metagenomics approaches, stable isotope probing is a promising tool to link the metabolic activity and diversity of microbial communities, including yet uncultured microbes, in a complex environment. Advancements in current stable isotope probing approaches integrated with the application of high-throughput diagnostic microarray-based phylogenetic profiling and metabolic flux analysis should facilitate the understanding of human microbial ecology and will enable the development of innovative strategies to treat or prevent intestinal diseases of as yet unknown etiology.

On the origins and biosynthesis of tetrodotoxin

The potent neurotoxin tetrodotoxin (TTX) has been identified from taxonomically diverse marine organisms. TTX possesses a unique cage-like structure, however, its biosynthesis has yet to be elucidated. Biosynthetic studies in the TTX-producing newt Taricha torosa, and in bacterial genera, including Vibrio, have proven inconclusive. Indeed, very few studies have been performed that address the cellular production of TTX. Here we review the sources of TTX described to date and provide evidence for the biosynthesis of TTX by symbiotic microorganisms in higher taxa. Chemical and genetic based biosynthesis studies of TTX undertaken thus far are discussed and we outline approaches which may be useful for expanding upon the current body of knowledge. The complex biosynthesis of structurally similar toxins, that reveal clues into the biosynthetic pathway of TTX, is also presented.

Construction of an E. Coli genome-scale atom mapping model for MFA calculations

Metabolic flux analysis (MFA) has so far been restricted to lumped networks lacking many important pathways, partly due to the difficulty in automatically generating isotope mapping matrices for genome-scale metabolic networks. Here we introduce a procedure that uses a compound matching algorithm based on the graph theoretical concept of pattern recognition along with relevant reaction information to automatically generate genome-scale atom mappings which trace the path of atoms from reactants to products for every reaction. The procedure is applied to the iAF1260 metabolic reconstruction of Escherichia coli yielding the genome-scale isotope mapping model imPR90068. This model maps 90,068 non-hydrogen atoms that span all 2,077 reactions present in iAF1260 (previous largest mapping model included 238 reactions). The expanded scope of the isotope mapping model allows the complete tracking of labeled atoms through pathways such as cofactor and prosthetic group biosynthesis and histidine metabolism. An EMU representation of imPR90068 is also constructed and made available.

Biosynthesis of andrographolide in Andrographis paniculata

Andrographolide, a diterpene lactone, is isolated from Andrographis paniculata which is well known for its medicinal properties. The biosynthetic route to andrographolide was studied using [1-(13)C]acetate, [2-(13)C]acetate and [1,6-(13)C(2)]glucose. The peak enrichment of eight carbon atoms in the (13)C NMR spectra of andrographolide suggested that deoxyxylulose pathway (DXP) is the major biosynthetic pathway to this diterpene. The contribution of the mevalonic acid pathway (MVA) is indicated by the observed (13)C-labeling pattern, and because the labeling patterns indicate a simultaneous contribution of both methyl erythritol phosphate (MEP) and MVA pathways it can be deduced that cross-talk occurs between plastids and cytoplasm.

Carbon metabolism of intracellular bacterial pathogens and possible links to virulence

New technologies such as high-throughput methods and 13C-isotopologue-profiling analysis are beginning to provide us with insight into the in vivo metabolism of microorganisms, especially in the host cell compartments that are colonized by intracellular bacterial pathogens. In this Review, we discuss the recent progress made in determining the major carbon sources and metabolic pathways used by model intracellular bacterial pathogens that replicate either in the cytosol or in vacuoles of infected host cells. Furthermore, we highlight the possible links between intracellular carbon metabolism and the expression of virulence genes.

Carbon metabolism of Listeria monocytogenes growing inside macrophages

The intracellular metabolism of Listeria monocytogenes was studied by (13)C-isotopologue profiling using murine J774A.1 macrophages as host cells. Six hours after infection, bacteria were separated from the macrophages and hydrolyzed. Amino acids were converted into tert-butyl-dimethylsilyl derivatives and subjected to gas chromatography/mass spectrometry. When the macrophages were supplied with [U-(13)C(6)]glucose prior to infection, but not during infection, label was detected only in Ala, Asp and Glu of the macrophage and bacterial protein with equal isotope distribution. When [U-(13)C(6)]glucose was provided during the infection period, (13)C label was found again in Ala, Asp and Glu from host and bacterial protein, but also in Ser, Gly, Thr and Val from the bacterial fraction. Mutants of L. monocytogenes defective in the uptake and catabolism of the C(3)-metabolites, glycerol and/or dihydroxyacetone, showed reduced incorporation of [U-(13)C(6)]glucose into bacterial amino acids under the same experimental settings. The (13)C pattern suggests that (i) significant fractions (50-100%) of bacterial amino acids were provided by the host cell, (ii) a C(3)-metabolite can serve as carbon source for L. monocytogenes under intracellular conditions and (iii) bacterial biosynthesis of Asp, Thr and Glu proceeds via oxaloacetate by carboxylation of pyruvate.

Gaining insight into microbial physiology in the large intestine: a special role for stable isotopes

The importance of the human large intestine for nutrition, health, and disease, is becoming increasingly realized. There are numerous indications of a distinct role for the gut in such important issues as immune disorders and obesity-linked diseases. Research on this long-neglected organ, which is colonized by a myriad of bacteria, is a rapidly growing field that is currently providing fascinating new insights into the processes going on in the colon, and their relevance for the human host. This review aims to give an overview of studies dealing with the physiology of the intestinal microbiota as it functions within and in interaction with the host, with a special focus on approaches involving stable isotopes. We have included general aspects of gut microbial life as well as aspects specifically relating to genomic, proteomic, and metabolomic studies. A special emphasis is further laid on reviewing relevant methods and applications of stable isotope-aided metabolic flux analysis (MFA). We argue that linking MFA with the '-omics' technologies using innovative modeling approaches is the way to go to establish a truly integrative and interdisciplinary approach. Systems biology thus actualized will provide key insights into the metabolic regulations involved in microbe-host mutualism and their relevance for health and disease.

13CO2 as a universal metabolic tracer in isotopologue perturbation experiments

A tobacco plant was illuminated for 5h in an atmosphere containing (13)CO(2) and then maintained for 10 days under standard greenhouse conditions. Nicotine, glucose, and amino acids from proteins were isolated chromatographically. Isotopologue abundances of isolated metabolites were determined quantitatively by NMR spectroscopy and mass spectrometry. The observed non-stochastic isotopologue patterns indicate (i) formation of multiply labeled photosynthetic carbohydrates during the (13)CO(2) pulse phase followed by (ii) partial catabolism of the primary photosynthetic products, and (iii) recombination of the (13)C-labeled fragments with unlabeled intermediary metabolites during the chase period. The detected and simulated isotopologue profiles of glucose and amino acids reflect carbon partitioning that is dominated by the Calvin cycle and glycolysis/glucogenesis. Retrobiosynthetic analysis of the nicotine pattern is in line with its known formation from nicotinic acid and putrescine via aspartate, glyceraldehyde phosphate and alpha-ketoglutarate as basic building blocks. The study demonstrates that pulse/chase labeling with (13)CO(2) as precursor is a powerful tool for the analysis of quantitative aspects of plant metabolism in completely unperturbed whole plants.

Metabolic flux elucidation for large-scale models using 13C labeled isotopes

A key consideration in metabolic engineering is the determination of fluxes of the metabolites within the cell. This determination provides an unambiguous description of metabolism before and/or after engineering interventions. Here, we present a computational framework that combines a constraint-based modeling framework with isotopic label tracing on a large scale. When cells are fed a growth substrate with certain carbon positions labeled with (13)C, the distribution of this label in the intracellular metabolites can be calculated based on the known biochemistry of the participating pathways. Most labeling studies focus on skeletal representations of central metabolism and ignore many flux routes that could contribute to the observed isotopic labeling patterns. In contrast, our approach investigates the importance of carrying out isotopic labeling studies using a more comprehensive reaction network consisting of 350 fluxes and 184 metabolites in Escherichia coli including global metabolite balances on cofactors such as ATP, NADH, and NADPH. The proposed procedure is demonstrated on an E. coli strain engineered to produce amorphadiene, a precursor to the antimalarial drug artemisinin. The cells were grown in continuous culture on glucose containing 20% [U-(13)C]glucose; the measurements are made using GC-MS performed on 13 amino acids extracted from the cells. We identify flux distributions for which the calculated labeling patterns agree well with the measurements alluding to the accuracy of the network reconstruction. Furthermore, we explore the robustness of the flux calculations to variability in the experimental MS measurements, as well as highlight the key experimental measurements necessary for flux determination. Finally, we discuss the effect of reducing the model, as well as shed light onto the customization of the developed computational framework to other systems.

Metabolic networks in motion: 13C-based flux analysis

Many properties of complex networks cannot be understood from monitoring the components—not even when comprehensively monitoring all protein or metabolite concentrations—unless such information is connected and integrated through mathematical models. The reason is that static component concentrations, albeit extremely informative, do not contain functional information per se. The functional behavior of a network emerges only through the nonlinear gene, protein, and metabolite interactions across multiple metabolic and regulatory layers. I argue here that intracellular reaction rates are the functional end points of these interactions in metabolic networks, hence are highly relevant for systems biology. Methods for experimental determination of metabolic fluxes differ fundamentally from component concentration measurements; that is, intracellular reaction rates cannot be detected directly, but must be estimated through computer model-based interpretation of stable isotope patterns in products of metabolism.

13C isotopologue perturbation studies of Listeria monocytogenes carbon metabolism and its modulation by the virulence regulator PrfA

The carbon metabolism of Listeria monocytogenes (Lm) EGD and the two isogenic mutant strains LmDeltaprfA and LmDeltaprfApPRFA* (showing no or enhanced expression, respectively, of the virulence factor PrfA) was determined by 13C isotopologue perturbation. After growth of the bacteria in a defined medium containing a mixture of [U-13C6]glucose and glucose with natural 13C abundance (1:25, wt/wt), 14 amino acids were isolated and analyzed by NMR spectroscopy. Multiply 13C-labeled isotopologues were determined quantitatively by signal deconvolution. The 13C enrichments and isotopologue patterns allowed the reconstruction of most amino acid biosynthesis pathways and illustrated that overproduced PrfA may strongly influence the synthesis of some amino acids, notably that of the branched amino acids (Val, Ile, and Leu). Retrobiosynthetic analysis of the isotopologue compositions showed that degradation of glucose occurs to a large extent via the pentose phosphate pathway and that the citrate cycle is incomplete because of the absence of 2-oxoglutarate dehydrogenase activity. The reconstructed labeling pattern of oxaloacetate indicated its formation by carboxylation of pyruvate. This metabolic reaction seems to have a strong impact on the growth requirement in defined minimal medium. Bioinformatical steady-state network analyses and flux distribution predictions confirmed the experimental data and predicted metabolite fluxes through the enzymes of the pathways under study.

Measuring multiple fluxes through plant metabolic networks

Fluxes through metabolic networks are crucial for cell function, and a knowledge of these fluxes is essential for understanding and manipulating metabolic phenotypes. Labeling provides the key to flux measurement, and in network flux analysis the measurement of multiple fluxes allows a flux map to be superimposed on the metabolic network. The principles and practice of two complementary methods, dynamic and steady-state labeling, are described, emphasizing best practice and illustrating their contribution to network flux analysis with examples taken from the plant and microbial literature. The principal analytical methods for the detection of stable isotopes are also described, as well as the procedures for obtaining flux maps from labeling data. A series of boxes summarizing the key concepts of network flux analysis is provided for convenience.

Beyond diversity: functional microbiomics of the human colon

Molecular tools have revealed wide microbial diversity in the human alimentary tract. Most intestinal microorganisms have not been cultured and the in situ functions of distinct groups of the intestinal microbiota are largely unknown but pivotal to understanding the role of these microorganisms in health and disease. Promising strategies to gain more insight into the functionality of the complex microbial communities in the human alimentary tract, including fermentation processes in the colon, are discussed. These research approaches could provide a basis for the definition of a healthy gut based on key properties of microbial functionality. This will also enable the development of direct nutritional strategies for intestinal disease prevention and health promotion.

Studies of the intermediary metabolism in cultured cells of the insect Spodoptera frugiperda using 13C- or 15N-labelled tracers

Background

Insect cells can serve as host systems for the recombinant expression of eukaryotic proteins. Using this platform, the controlled expression of ¹⁵N/¹³C labelled proteins requires the analysis of incorporation paths and rates of isotope-labelled precursors present in the medium into amino acids. For this purpose, Spodoptera frugiperda cells were grown in a complex medium containing [U-¹³C₆]glucose. In a second experiment, cultures of S. frugiperda were grown in the presence of ¹⁵N-phenylalanine.

Results

Quantitative NMR analysis showed incorporation of the proffered [U-¹³C₆]glucose into the ribose moiety of ribonucleosides (40 – 45%) and into the amino acids, alanine (41%), glutamic acid/glutamine (C-4 and C-5, 30%) and aspartate/asparagine (15%). Other amino acids and the purine ring of nucleosides were not formed from exogenous glucose in significant amounts (> 5%). Prior to the incorporation into protein the proffered ¹⁵N-phenylalanine lost about 70% of its label by transamination and the labelled compound was not converted into tyrosine to a significant extent.

Conclusion

Growth of S. frugiperda cells in the presence of [U-¹³C₆]glucose is conducive to the fractional labelling of ribonucleosides, alanine, glutamic acid/glutamine and aspartic acid/asparagine. The isotopolog compositions of the ribonucleosides and of alanine indicate considerable recycling of carbohydrate intermediates in the reductive branch of the pentose phosphate pathway. The incorporation of ¹⁵N-labelled amino acids may be hampered by loss of the ¹⁵N-label by transamination.

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.

Biosynthesis of benzofuran derivatives in root cultures of Tagetes patula via phenylalanine and 1-deoxy-D-xylulose 5-phosphate

Root cultures of Tagetes patula L. cv. Carmen were grown with a mixture of unlabeled glucose and [U-(13)C(6)]glucose or [1-(13)C(1)]glucose as carbon source. Isoeuparin and (-)-4-hydroxytremetone were isolated by solvent extraction of the cultured tissue, purified by chromatography and analysed by (1)H and (13)C NMR spectroscopy. Amino acids obtained by hydrolysis of protein from the same experiments were used for the reconstruction of the labelling patterns in central metabolic intermediates. These labelling patterns were used for the prediction of isotopolog compositions in the benzofuranone derivatives via different hypothetical pathways. Comparison with the experimentally observed isotopolog distributions showed that the benzenoid ring and the acetoxy group are exclusively or predominantly (>98%) derived from phenylalanine and not from acetyl-CoA via a polyketide-type biosynthesis. The isopropylidene side chain and two carbon atoms of the furan and dihydrofuran moiety, respectively, originate from an isoprenoid building block obtained exclusively or predominantly (>98%) via the deoxyxylulose phosphate pathway. The exomethylene atom of the isopropylidene side chain is biosynthetically equivalent to the (Z)-methyl group of dimethylallyl diphosphate. The data indicate that isoeuparin and (-)-4-hydroxytremetone are assembled from 4-hydroxyacetophenone and dimethylallyl diphosphate via prenyl-substituted 4-hydroxyacetophenone and dihydrobenzofurans as intermediates.

Revealing metabolic phenotypes in plants: inputs from NMR analysis

Assessing the performance of the plant metabolic network, with its varied biosynthetic capacity and its characteristic subcellular compartmentation, remains a considerable challenge. The complexity of the network is such that it is not yet possible to build large-scale predictive models of the fluxes it supports, whether on the basis of genomic and gene expression analysis or on the basis of more traditional measurements of metabolites and their interconversions. This limits the agronomic and biotechnological exploitation of plant metabolism, and it undermines the important objective of establishing a rational metabolic engineering strategy. Metabolic analysis is central to removing this obstacle and currently there is particular interest in harnessing high-throughput and/or large-scale analyses to the task of defining metabolic phenotypes. Nuclear magnetic resonance (NMR) spectroscopy contributes to this objective by providing a versatile suite of analytical techniques for the detection of metabolites and the fluxes between them. The principles that underpin the analysis of plant metabolism by NMR are described, including a discussion of the measurement options for the detection of metabolites in vivo and in vitro, and a description of the stable isotope labelling experiments that provide the basis for metabolic flux analysis. Despite a relatively low sensitivity, NMR is suitable for high-throughput system-wide analyses of the metabolome, providing methods for both metabolite fingerprinting and metabolite profiling, and in these areas NMR can contribute to the definition of plant metabolic phenotypes that are based on metabolic composition. NMR can also be used to investigate the operation of plant metabolic networks. Labelling experiments provide information on the operation of specific pathways within the network, and the quantitative analysis of steady-state labelling experiments leads to the definition of large-scale flux maps for heterotrophic carbon metabolism. These maps define multiple unidirectional fluxes between branch-points in the metabolic network, highlighting the existence of substrate cycles and discriminating in favourable cases between fluxes in the cytosol and plastid. Flux maps can be used to define a functionally relevant metabolic phenotype and the extensive application of such maps in microbial systems suggests that they could have important applications in characterising the genotypes produced by plant genetic engineering.

Quantification of compartmented metabolic fluxes in developing soybean embryos by employing biosynthetically directed fractional (13)C labeling, two-dimensional [(13)C, (1)H] nuclear magnetic resonance, and comprehensive isotopomer balancing

Metabolic flux quantification in plants is instrumental in the detailed understanding of metabolism but is difficult to perform on a systemic level. Toward this aim, we report the development and application of a computer-aided metabolic flux analysis tool that enables the concurrent evaluation of fluxes in several primary metabolic pathways. Labeling experiments were performed by feeding a mixture of U-(13)C Suc, naturally abundant Suc, and Gln to developing soybean (Glycine max) embryos. Two-dimensional [(13)C, (1)H] NMR spectra of seed storage protein and starch hydrolysates were acquired and yielded a labeling data set consisting of 155 (13)C isotopomer abundances. We developed a computer program to automatically calculate fluxes from this data. This program accepts a user-defined metabolic network model and incorporates recent mathematical advances toward accurate and efficient flux evaluation. Fluxes were calculated and statistical analysis was performed to obtain sds. A high flux was found through the oxidative pentose phosphate pathway (19.99 +/- 4.39 micromol d(-1) cotyledon(-1), or 104.2 carbon mol +/- 23.0 carbon mol per 100 carbon mol of Suc uptake). Separate transketolase and transaldolase fluxes could be distinguished in the plastid and the cytosol, and those in the plastid were found to be at least 6-fold higher. The backflux from triose to hexose phosphate was also found to be substantial in the plastid (21.72 +/- 5.00 micromol d(-1) cotyledon(-1), or 113.2 carbon mol +/-26.0 carbon mol per 100 carbon mol of Suc uptake). Forward and backward directions of anaplerotic fluxes could be distinguished. The glyoxylate shunt flux was found to be negligible. Such a generic flux analysis tool can serve as a quantitative tool for metabolic studies and phenotype comparisons and can be extended to other plant systems.

Elucidation of the biosynthetic pathway of the allelochemical sorgoleone using retrobiosynthetic NMR analysis

NMR analyses of the labeling pattern obtained using various 13C-labeled precursors indicated that both the lipid tail and the quinone head of sorgoleone, the main allelopathic component of the oily root exudate of Sorghum bicolor, were derived from acetate units, but that the two moieties were synthesized in different subcellular compartments. The 16:3 fatty acid precursor of the tail is synthesized by the combined action of fatty-acid synthase and desaturases most likely in the plastids. It is then exported out of the plastids and converted to 5-pentadecatriene resorcinol by a polyketide synthase. This resorcinol intermediate was identified in root hair extracts. The lipid resorcinol intermediate is then methylated by a S-adenosylmethionine-dependent O-methyltransferase and subsequently dihydroxylated by a P450 monooxygenase to yield the reduced form of sorgoleone.

Phenylalanine-independent biosynthesis of 1,3,5,8-tetrahydroxyxanthone. A retrobiosynthetic NMR study with root cultures of Swertia chirata

Root cultures of Swertia chirata (Gentianaceae) were grown with supplements of [1-13C]glucose, [U-13C6]glucose or [carboxy-13C]shikimic acid. 1,3,5,8-Tetrahydroxyxanthone was isolated and analysed by quantitative NMR analysis. The observed isotopomer distribution shows that 1,3,5,8-tetrahydroxyxanthone is biosynthesized via a polyketide-type pathway. The starter unit, 3-hydroxybenzoyl-CoA, is obtained from an early shikimate pathway intermediate. Phenylalanine, cinnamic acid and benzoic acid were ruled out as intermediates.

Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS

We describe here a novel methodology for rapid diagnosis of metabolic changes, which is based on probabilistic equations that relate GC-MS-derived mass distributions in proteinogenic amino acids to in vivo enzyme activities. This metabolic flux ratio analysis by GC-MS provides a comprehensive perspective on central metabolism by quantifying 14 ratios of fluxes through converging pathways and reactions from [1-13C] and [U-13C]glucose experiments. Reliability and accuracy of this method were experimentally verified by successfully capturing expected flux responses of Escherichia coli to environmental modifications and seven knockout mutations in all major pathways of central metabolism. Furthermore, several mutants exhibited additional, unexpected flux responses that provide new insights into the behavior of the metabolic network in its entirety. Most prominently, the low in vivo activity of the Entner-Doudoroff pathway in wild-type E. coli increased up to a contribution of 30% to glucose catabolism in mutants of glycolysis and TCA cycle. Moreover, glucose 6-phosphate dehydrogenase mutants catabolized glucose not exclusively via glycolysis, suggesting a yet unidentified bypass of this reaction. Although strongly affected by environmental conditions, a stable balance between anaplerotic and TCA cycle flux was maintained by all mutants in the upper part of metabolism. Overall, our results provide quantitative insight into flux changes that bring about the resilience of metabolic networks to disruption.

Biosynthesis of pteridines. Reaction mechanism of GTP cyclohydrolase I

GTP cyclohydrolase I catalyses the hydrolytic release of formate from GTP followed by cyclization to dihydroneopterin triphosphate. The enzymes from bacteria and animals are homodecamers containing one zinc ion per subunit. Replacement of Cys110, Cys181, His112 or His113 of the enzyme from Escherichia coli by serine affords catalytically inactive mutant proteins with reduced capacity to bind zinc. These mutant proteins are unable to convert GTP or the committed reaction intermediate, 2-amino-5-formylamino-6-(beta-ribosylamino)-4(3H)-pyrimidinone 5'-triphosphate, to dihydroneopterin triphosphate. The crystal structures of GTP complexes of the His113Ser, His112Ser and Cys181Ser mutant proteins determined at resolutions of 2.5A, 2.8A and 3.2A, respectively, revealed the conformation of substrate GTP in the active site cavity. The carboxylic group of the highly conserved residue Glu152 anchors the substrate GTP, by hydrogen bonding to N-3 and to the position 2 amino group. Several basic amino acid residues interact with the triphosphate moiety of the substrate. The structure of the His112Ser mutant in complex with an undefined mixture of nucleotides determined at a resolution of 2.1A afforded additional details of the peptide folding. Comparison between the wild-type and mutant enzyme structures indicates that the catalytically active zinc ion is directly coordinated to Cys110, Cys181 and His113. Moreover, the zinc ion is complexed to a water molecule, which is in close hydrogen bond contact to His112. In close analogy to zinc proteases, the zinc-coordinated water molecule is suggested to attack C-8 of the substrate affording a zinc-bound 8R hydrate of GTP. Opening of the hydrated imidazole ring affords a formamide derivative, which remains coordinated to zinc. The subsequent hydrolysis of the formamide motif has an absolute requirement for zinc ion catalysis. The hydrolysis of the formamide bond shows close mechanistic similarity with peptide hydrolysis by zinc proteases.

Dissimilation of [(13)C]methanol by continuous cultures of Bacillus methanolicus MGA3 at 50 degrees C studied by (13)C NMR and isotope-ratio mass spectrometry

Using a continuous culture of Bacillus methanolicus MGA3 limited by 100 mM methanol in the feed and growing at a dilution rate D=0.25 h(-1), transients in dissolved methanol were studied to determine the effects of methanol toxicity and the pathway of methanol dissimilation to CO(2). Steady-state cultures were disturbed by pulses of methanol resulting in a rapid change in concentration of 6.4-12.8 mM. B. methanolicus MGA3 responded to a sudden increase in available methanol by a transient decline in the biomass concentration in the reactor. In most cases the culture returned to steady state between 4 and 12 h after pulse addition. However, at a methanol pulse of 12.8 mM, complete biomass washout occurred and the culture did not return to steady state. Integrating the response curves of the dry biomass concentration over a 12 h time period showed that a methanol pulse can cause an average transient decline in the biomass yield of up to 22%. (13)C NMR experiments using labelled methanol indicated that the transient partial or complete biomass washout was probably caused by toxic accumulation of formaldehyde in the culture. These experiments also showed accumulation of formate, indicating that B. methanolicus possesses formaldehyde dehydrogenase and formate dehydrogenase activity resulting in a methanol dissimilation pathway via formate to CO(2). Studies using isotope-ratio mass spectrometry provided further evidence of a methanol dissimilation pathway via formate. B. methanolicus MGA3, growing continuously under methanol limitation, consumed added formate at a rate of approximately 0.85 mmol l(-1) h(-1). Furthermore, significant accumulation of (13)CO(2) in the reactor exhaust gas was measured in response to a pulse addition of [(13)C]formic acid to the bioreactor. This indicates that B. methanolicus dissimilates methanol carbon to CO(2) in order to detoxify formaldehyde by both a linear pathway to formate and a cyclic mechanism as part of the RuMP pathway.

Metabolic engineering for L-lysine production by Corynebacterium glutamicum

Corynebacterium glutamicum has been used since several decades for the large-scale production of amino acids, esp. L-glutamate and L-lysine. After initial successes of random mutagenesis and screening approaches, further strain improvements now require a much more rational design, i.e. metabolic engineering. Not only recombinant DNA technology but also mathematical modelling of metabolism as well as metabolic flux analysis represent important metabolic engineering tools. This review covers as state-of-the-art examples of these techniques the genetic engineering of the L-lysine biosynthetic pathway resulting in a vectorless strain with significantly increased dihydrodipicolinate synthase activity, and the detailed metabolic flux analysis by 13C isotopomer labelling strategies of the anaplerotic enzyme activities in C. glutamicum resulting in the identification of gluconeogenic phosphoenolpyruvate carboxykinase as a limiting enzyme.

Autotrophic CO(2) fixation by Chloroflexus aurantiacus: study of glyoxylate formation and assimilation via the 3-hydroxypropionate cycle

In the facultative autotrophic organism Chloroflexus aurantiacus, a phototrophic green nonsulfur bacterium, the Calvin cycle does not appear to be operative in autotrophic carbon assimilation. An alternative cyclic pathway, the 3-hydroxypropionate cycle, has been proposed. In this pathway, acetyl coenzyme A (acetyl-CoA) is assumed to be converted to malate, and two CO(2) molecules are thereby fixed. Malyl-CoA is supposed to be cleaved to acetyl-CoA, the starting molecule, and glyoxylate, the carbon fixation product. Malyl-CoA cleavage is shown here to be catalyzed by malyl-CoA lyase; this enzyme activity is induced severalfold in autotrophically grown cells. Malate is converted to malyl-CoA via an inducible CoA transferase with succinyl-CoA as a CoA donor. Some enzyme activities involved in the conversion of malonyl-CoA via 3-hydroxypropionate to propionyl-CoA are also induced under autotrophic growth conditions. So far, no clue as to the first step in glyoxylate assimilation has been obtained. One possibility for the assimilation of glyoxylate involves the conversion of glyoxylate to glycine and the subsequent assimilation of glycine. However, such a pathway does not occur, as shown by labeling of whole cells with [1,2-(13)C(2)]glycine. Glycine carbon was incorporated only into glycine, serine, and compounds that contained C(1) units derived therefrom and not into other cell compounds.

Tracer studies with 13C-labeled carbohydrates in cultured plant cells. Retrobiosynthetic analysis of chelidonic acid biosynthesis

The biosynthesis of chelidonic acid was studied in cell suspension cultures of Leucojum aestivum. Cell cultures were supplied with [U-13C]glucose, [l-13C]glucose or [U-13Cs]ribose/ribulose in standard medium containing unlabeled glucose. 13C labeling patterns of amino acids obtained by hydrolysis of biomass were determined by NMR spectroscopy and compared to the labeling pattern of chelidonic acid. The data document the incorporation of a contiguous 4-carbon fragment derived from the pentose phosphate pool into chelidonic acid. This suggests a biosynthetic pathway involving the condensation of phosphoenolpyruvate with a pentose phosphate followed by dehydration, dehydrogenation, ring closure and decarboxylation conducive to the loss of C-5 of the pentose precursor.

Retrobiosynthetic nuclear magnetic resonance analysis of amino acid biosynthesis and intermediary metabolism. Metabolic flux in developing maize kernels

Information on metabolic networks could provide the basis for the design of targets for metabolic engineering. To study metabolic flux in cereals, developing maize (Zea mays) kernels were grown in sterile culture on medium containing [U-(13)C(6)]glucose or [1,2-(13)C(2)]acetate. After growth, amino acids, lipids, and sitosterol were isolated from kernels as well as from the cobs, and their (13)C isotopomer compositions were determined by quantitative nuclear magnetic resonance spectroscopy. The highly specific labeling patterns were used to analyze the metabolic pathways leading to amino acids and the triterpene on a quantitative basis. The data show that serine is generated from phosphoglycerate, as well as from glycine. Lysine is formed entirely via the diaminopimelate pathway and sitosterol is synthesized entirely via the mevalonate route. The labeling data of amino acids and sitosterol were used to reconstruct the labeling patterns of key metabolic intermediates (e.g. acetyl-coenzyme A, pyruvate, phosphoenolpyruvate, erythrose 4-phosphate, and Rib 5-phosphate) that revealed quantitative information about carbon flux in the intermediary metabolism of developing maize kernels. Exogenous acetate served as an efficient precursor of sitosterol, as well as of amino acids of the aspartate and glutamate family; in comparison, metabolites formed in the plastidic compartments showed low acetate incorporation.

Biosynthesis of the bicycloalternarenes, mixed terpenoids of Alternaria alternata

By incorporation of [2-13C]-mevalonate, [1-13C]-acetate and [1-13C]-glucose we could reveal that the phytopathogenic fungus Alternaria alternata biosynthesized the mixed terpenoids bicycloalternarenes via the classic mevalonate pathway. The polyketid pathway does not participate in the biosynthesis of bicycloalternarenes, because there is no incorporation of [13C]-acetate into the C-ring of these compounds. The labelling pattern in this nonterpenoid part of bicycloalternarenes after feeding with [1-13C]-glucose and [U-13C6]-glucose, respectively, allows the assumption that metabolites of the shikimate pathway are involved.

Pathway analysis and metabolic engineering in Corynebacterium glutamicum

The gram-positive bacterium Corynebacterium glutamicum is used for the industrial production of amino acids, e.g. of L-glutamate and L-lysine. During the last 15 years, genetic engineering and amplification of genes have become fascinating methods for studying metabolic pathways in greater detail and for the construction of strains with the desired genotypes. In order to obtain a better understanding of the central metabolism and to quantify the in vivo fluxes in C. glutamicum, the [13C]-labelling technique was combined with metabolite balancing to achieve a unifying comprehensive pathway analysis. These methods can determine the flux distribution at the branch point between glycolysis and the pentose phosphate pathway. The in vivo fluxes in the oxidative part of the pentose phosphate pathway calculated on the basis of intracellular metabolite concentrations and the kinetic constants of the purified glucose-6-phosphate and 6-phosphogluconate dehydrogenases determined in vitro were in full accordance with the fluxes measured by the [13C]-labelling technique. These data indicate that the oxidative pentose phosphate pathway in C. glutamicum is mainly regulated by the ratio of NADPH/NADP concentrations and the specific activity of glucose-6-phosphate dehydrogenase. The carbon flux via the oxidative pentose phosphate pathway correlated with the NADPH demand for L-lysine synthesis. Although it has generally been accepted that phosphoenolpyruvate carboxylase fulfills a main anaplerotic function in C. glutamicum, we recently detected that a biotin-dependent pyruvate carboxylase exists as a further anaplerotic enzyme in this bacterium. In addition to the activities of these two carboxylases three enzymes catalysing the decarboxylation of the C4 metabolites oxaloacetate or malate are also present in this bacterium. The individual flux rates at this complex anaplerotic node were investigated by using [13C]-labelled substrates. The results indicate that both carboxylation and decarboxylation occur simultaneously in C. glutamicum so that a high cyclic flux of oxaloacetate via phosphoenolpyruvate to pyruvate was found. Furthermore, we detected that in C. glutamicum two biosynthetic pathways exist for the synthesis of DL-diaminopimelate and L-lysine. As shown by NMR spectroscopy the relative use of both pathways in vivo is dependent on the ammonium concentration in the culture medium. Mutants defective in one pathway are still able to synthesise enough L-lysine for growth, but the L-lysine yields with overproducers were reduced. The luxury of having these two pathways gives C. glutamicum an increased flexibility in response to changing environmental conditions and is also related to the essential need for DL-diaminopimelate as a building block for the synthesis of the murein sacculus.

Auxin biosynthesis in maize kernels

Auxin biosynthesis was analyzed in a maize (Zea mays) kernel culture system in which the seeds develop under physiological conditions similar to the in vivo situation. This system was modified for precursor feeding experiments. Tryptophan (Trp) is efficiently incorporated into indole-3-acetic acid (IAA) with retention of the 3, 3' bond. Conversion of Trp to IAA is not competed by indole. Labeling with the general precursors [U-(13)C(6)]glucose and [1, 2-(13)C(2)]acetate followed by retrobiosynthetic analysis strongly suggest that Trp-dependent IAA synthesis is the predominant route for auxin biosynthesis in the maize kernel. The synthesis of IAA from indole glycerol phosphate and IAA formation via condensation of indole with an acetyl-coenzyme A or phosphoenolpyruvate derived metabolite can be excluded.

Biosynthesis of a neo-epi-verrucosane diterpene in the liverwort Fossombronia alaskana. A retrobiosynthetic NMR study

The biosynthesis of the diterpene 8alpha-acetoxy-13alpha-hydroxy-5-oxo-13-epi- neoverrucosane in the arctic liverwort Fossombronia alaskana was studied by incorporation experiments using [1-(13)C]- and [U-(13)C(6)]glucose as precursors. The (13)C-labeling patterns of acetyl-CoA, pyruvate, and phosphoenolpyruvate in intermediary metabolism were reconstructed from the (13)C NMR data of biosynthetic amino acids (leucine, alanine, phenylalanine) and were used to predict hypothetical labeling patterns for isopentenyl pyrophosphate formed via the mevalonate pathway and the deoxyxylulose pathway. The labeling patterns observed for the neoverrucosane diterpene were consistent with the intermediate formation of geranyllinaloyl pyrophosphate assembled from dimethylallyl pyrophosphate and three molecules of isopentenyl pyrophosphate generated predominantly or entirely via 1-deoxyxylulose 5-phosphate. The experimental data can be integrated into a detailed biosynthetic scheme involving a 1,5-hydride shift. The postulated involvement of the 1,5-hydride shift was confirmed by an incorporation experiment with [6,6-(2)H(2)]glucose.