From fluxes and isotope labeling patterns towards in silico cells

Stereoselective accumulation and biotransformation of chiral fungicide epoxiconazole and oxidative stress, detoxification, and endogenous metabolic disturbance in earthworm (Eisenia foetida)

>80 % of applied pesticides in agriculture will enter the soil and be exposed to soil animals. Little is known about the stereoselective metabolic effects of epoxiconazole (EPO) on soil animals. In this study, EPO-mediated stereoselective enrichment, biotransformation, oxidative stress, detoxification, and global metabolic profiles in earthworms were investigated by exposure to EPO and its enantiomers at 1 mg/kg and 10 mg/kg doses. Preferential enrichment of (-)-EPO was observed, and the five transformation products (TPs) exhibited the chemically specific stereoselective accumulation with inconsistent configurations. Biochemical markers related to reactive oxygen species (ROS) and detoxification (·OH^- content, SOD, CAT, GST, and CYP450 enzymes) showed a significant stereoselective activation overall at the low-level exposure (p-value <0.05). Based on untargeted metabolomic analysis, the steroid biosynthesis and ROS-related biotransformation, glutathione metabolism, TCA cycle, amino acid metabolism, purine and pyrimidine metabolism of earthworms were significantly interfered with by EPO and its enantiomer exposure. More pronounced stereoselectivity was observed at the level of the global metabolic profile, while comparable levels of metabolic perturbations were identified at the individual metabolite level. This study provides novel insights into the stereoselective effects of the chiral fungicide EPO, and valuable evidence for soil environmental risk assessments.

Fluxomics - New Metabolomics Approaches to Monitor Metabolic Pathways

Fluxomics is an innovative -omics research field that measures the rates of all intracellular fluxes in the central metabolism of biological systems. Fluxomics gathers data from multiple different -omics fields, portraying the whole picture of molecular interactions. Recently, fluxomics has become one of the most relevant approaches to investigate metabolic phenotypes. Metabolic flux using 13C-labeled molecules is increasingly used to monitor metabolic pathways, to probe the corresponding gene-RNA and protein-metabolite interaction networks in actual time. Thus, fluxomics reveals the functioning of multi-molecular metabolic pathways and is increasingly applied in biotechnology and pharmacology. Here, we describe the main fluxomics approaches and experimental platforms. Moreover, we summarize recent fluxomic results in different biological systems.

Dynamic metabolic response of Aspergillus niger to glucose perturbation: evidence of regulatory mechanism for reduced glucoamylase production

Environmental gradient is an important common issue during scale-up process for protein production. To address the dynamic regulatory mechanism of Aspergillus niger being exposed to inhomogeneous glucose concentrations, glucose perturbation were experimented on the steady state of A. niger chemostat culture, and dynamic profiles of the intracellular metabolites in central carbon metabolism were tracked in a time scale of seconds. The upper glycolysis and pentose phosphate pathway showed sharp variations after glucose perturbation, while the lower glycolysis, TCA cycle and amino acid pools represented a moderate and prolonged response due to the allosteric regulation of enzymes and buffering function of metabolites with large pool sizes. Improved glucose-6-phosphate enhanced the metabolic flux to PP pathway remarkably, which provided not only more redox cofactors (NADPH) for protein synthesis but also more precursors (phosphoribosyl pyrophosphate and ribose-5-phosphate) for cell growth. Moreover, reduction of the total adenine nucleotides and major precursor amino acids indicated the upregulated RNA synthesis was required to produce stress proteins, and partially explained the drop of glucoamylase production when A. niger experienced a fluctuated glucose concentration environment. These findings would be valuable for improving bioreactor operation, design, and scale-up from engineering or genetic aspects.

Metabolic adjustment upon repetitive substrate perturbations using dynamic 13C-tracing in yeast

BACKGROUND

Natural and industrial environments are dynamic with respect to substrate availability and other conditions like temperature and pH. Especially, metabolism is strongly affected by changes in the extracellular space. Here we study the dynamic flux of central carbon metabolism and storage carbohydrate metabolism under dynamic feast/famine conditions in Saccharomyces cerevisiae.

RESULTS

The metabolic flux reacts fast and sensitive to cyclic perturbations in substrate availability. Compared to well-documented stimulus-response experiments using substrate pulses, different metabolic responses are observed. Especially, cells experiencing cyclic perturbations do not show a drop in ATP with the addition of glucose, but an immediate increase in energy charge. Although a high glycolytic flux of up to 5.4 mmol g DW-1  h-1 is observed, no overflow metabolites are detected. From famine to feast the glucose uptake rate increased from 170 to 4788 μmol g DW-1  h-1 in 24 s. Intracellularly, even more drastic changes were observed. Especially, the T6P synthesis rate increased more than 100-fold upon glucose addition. This response indicates that the storage metabolism is very sensitive to changes in glycolytic flux and counterbalances these rapid changes by diverting flux into large pools to prevent substrate accelerated death and potentially refill the central metabolism when substrates become scarce. Using 13C-tracer we found a dilution in the labeling of extracellular glucose, G6P, T6P and other metabolites, indicating an influx of unlabeled carbon. It is shown that glycogen and trehalose degradation via different routes could explain these observations. Based on the 13C labeling in average 15% of the carbon inflow is recycled via trehalose and glycogen. This average fraction is comparable to the steady-state turnover, but changes significantly during the cycle, indicating the relevance for dynamic regulation of the metabolic flux.

CONCLUSIONS

Comparable to electric energy grids, metabolism seems to use storage units to buffer peaks and keep reserves to maintain a robust function. During the applied fast feast/famine conditions about 15% of the metabolized carbon were recycled in storage metabolism. Additionally, the resources were distributed different to steady-state conditions. Most remarkably is a fivefold increased flux towards PPP that generated a reversed flux of transaldolase and the F6P-producing transketolase reactions. Combined with slight changes in the biomass composition, the yield decrease of 5% can be explained.

Identifying model error in metabolic flux analysis - a generalized least squares approach

BACKGROUND

The estimation of intracellular flux through traditional metabolic flux analysis (MFA) using an overdetermined system of equations is a well established practice in metabolic engineering. Despite the continued evolution of the methodology since its introduction, there has been little focus on validation and identification of poor model fit outside of identifying "gross measurement error". The growing complexity of metabolic models, which are increasingly generated from genome-level data, has necessitated robust validation that can directly assess model fit.

RESULTS

In this work, MFA calculation is framed as a generalized least squares (GLS) problem, highlighting the applicability of the common t-test for model validation. To differentiate between measurement and model error, we simulate ideal flux profiles directly from the model, perturb them with estimated measurement error, and compare their validation to real data. Application of this strategy to an established Chinese Hamster Ovary (CHO) cell model shows how fluxes validated by traditional means may be largely non-significant due to a lack of model fit. With further simulation, we explore how t-test significance relates to calculation error and show that fluxes found to be non-significant have 2-4 fold larger error (if measurement uncertainty is in the 5-10 % range).

CONCLUSIONS

The proposed validation method goes beyond traditional detection of "gross measurement error" to identify lack of fit between model and data. Although the focus of this work is on t-test validation and traditional MFA, the presented framework is readily applicable to other regression analysis methods and MFA formulations.

Achieving Metabolic Flux Analysis for S. cerevisiae at a Genome-Scale: Challenges, Requirements, and Considerations

Recent advances in 13C-Metabolic flux analysis (13C-MFA) have increased its capability to accurately resolve fluxes using a genome-scale model with narrow confidence intervals without pre-judging the activity or inactivity of alternate metabolic pathways. However, the necessary precautions, computational challenges, and minimum data requirements for successful analysis remain poorly established. This review aims to establish the necessary guidelines for performing 13C-MFA at the genome-scale for a compartmentalized eukaryotic system such as yeast in terms of model and data requirements, while addressing key issues such as statistical analysis and network complexity. We describe the various approaches used to simplify the genome-scale model in the absence of sufficient experimental flux measurements, the availability and generation of reaction atom mapping information, and the experimental flux and metabolite labeling distribution measurements to ensure statistical validity of the obtained flux distribution. Organism-specific challenges such as the impact of compartmentalization of metabolism, variability of biomass composition, and the cell-cycle dependence of metabolism are discussed. Identification of errors arising from incorrect gene annotation and suggested alternate routes using MFA are also highlighted.

OpenFLUX2: (13)C-MFA modeling software package adjusted for the comprehensive analysis of single and parallel labeling experiments

BACKGROUND

Steady-state (13)C-based metabolic flux analysis ((13)C-MFA) is the most powerful method available for the quantification of intracellular fluxes. These analyses include concertedly linked experimental and computational stages: (i) assuming the metabolic model and optimizing the experimental design; (ii) feeding the investigated organism using a chosen (13)C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the (13)C-patterns of intracellular metabolites; and (iv) computing flux parameters that minimize the differences between observed and simulated measurements, followed by evaluating flux statistics. In its early stages, (13)C-MFA was performed on the basis of data obtained in a single labeling experiment (SLE) followed by exploiting the developed high-performance computational software. Recently, the advantages of parallel labeling experiments (PLEs), where several LEs are conducted under the conditions differing only by the tracer(s) choice, were demonstrated, particularly with regard to improving flux precision due to the synergy of complementary information. The availability of an open-source software adjusted for PLE-based (13)C-MFA is an important factor for PLE implementation.

RESULTS

The open-source software OpenFLUX, initially developed for the analysis of SLEs, was extended for the computation of PLE data. Using the OpenFLUX2, in silico simulation confirmed that flux precision is improved when (13)C-MFA is implemented by fitting PLE data to the common model compared with SLE-based analysis. Efficient flux resolution could be achieved in the PLE-mediated analysis when the choice of tracer was based on an experimental design computed to minimize the flux variances from different parts of the metabolic network. The analysis provided by OpenFLUX2 mainly includes (i) the optimization of the experimental design, (ii) the computation of the flux parameters from LEs data, (iii) goodness-of-fit testing of the model's adequacy, (iv) drawing conclusions concerning the identifiability of fluxes and construction of a contribution matrix reflecting the relative contribution of the measurement variances to the flux variances, and (v) precise determination of flux confidence intervals using a fine-tunable and convergence-controlled Monte Carlo-based method.

CONCLUSIONS

The developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state (13)C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.

Isotope labeling pattern study of central carbon metabolites using GC/MS

Determination of fluxes by (13)C tracer experiments depends on monitoring the (13)C labeling pattern of metabolites during isotope experiments. In metabolome-based (13)C metabolic flux analysis, liquid chromatography combined with mass spectrometry or tandem mass spectrometry (LC/MS or LC/MS/MS, respectively) has been mainly used as an analytical platform for isotope pattern studies of central carbon metabolites. However, gas chromatography with mass spectrometry (GC/MS) has several advantages over LC/MS, such as high sensitivity, low cost, ease of operation, and availability of mass spectra databases for comparison. In this study, analysis of isotope pattern for central carbon metabolites using GC/MS was demonstrated. First, a proper set of mass ions for central carbon metabolites was selected based on carbon backbone information and structural isomers of mass fragment ions. A total of 34 mass fragment ions was selected and used for the quantification of 25 central carbon metabolites. Then, to quantify isotope fractions, a natural mass isotopomer library for selected mass fragment ions was constructed and subtracted from isotopomer mass spectra data. The results revealed a surprisingly high abundance of partially labeled (13)C intermediates, such as 56.4% of fructose 6-phosphate and 47.6% of dihydroxyacetone phosphate at isotopic steady state, which were generated in the pentose phosphate pathway. Finally, dynamic changes of isotope fragments of central metabolites were monitored with a U-(13)C glucose stimulus response experiment in Kluyveromyces marxianus. With a comprehensive study of isotope patterns of central carbon metabolites using GC/MS, 25 central carbon metabolites and their isotopic fractions were successfully quantified. Dynamic and precise acquisition of isotope pattern can then be used in combination with proper kinetic models to calculate metabolic fluxes.

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.

Exploring the metabolic state of microorganisms using metabolomics

Microorganisms depend on their ability to modulate their metabolic composition according to specific circumstances, such as different phases of the growth cycle and circadian rhythms, fluctuations in environmental conditions, as well as experimental perturbations. A thorough understanding of these metabolic adaptations requires the ability to comprehensively identify and quantify the metabolome of bacterial cells in different states. In this review, we present an overview of the diverse metabolomics approaches recently adopted to explore the metabolism of a wide variety of microorganisms. Focusing on a selection of illustrative case studies, we assess the different experimental designs used and explore the major achievements and remaining challenges in the field. We conclude by discussing the important complementary information provided by computational methods such as genome-scale metabolic modeling, which enable an integrated analysis of metabolic state changes in the context of overall cellular physiology.

Advanced stoichiometric analysis of metabolic networks of mammalian systems

Metabolic engineering tools have been widely applied to living organisms to gain a comprehensive understanding about cellular networks and to improve cellular properties. Metabolic flux analysis (MFA), flux balance analysis (FBA), and metabolic pathway analysis (MPA) are among the most popular tools in stoichiometric network analysis. Although application of these tools into well-known microbial systems is extensive in the literature, various barriers prevent them from being utilized in mammalian cells. Limited experimental data, complex regulatory mechanisms, and the requirement of more complex nutrient media are some major obstacles in mammalian cell systems. However, mammalian cells have been used to produce therapeutic proteins, to characterize disease states or related abnormal metabolic conditions, and to analyze the toxicological effects of some medicinally important drugs. Therefore, there is a growing need for extending metabolic engineering principles to mammalian cells in order to understand their underlying metabolic functions. In this review article, advanced metabolic engineering tools developed for stoichiometric analysis including MFA, FBA, and MPA are described. Applications of these tools in mammalian cells are discussed in detail, and the challenges and opportunities are highlighted.

Stable isotope-resolved metabolomics and applications for drug development

Advances in analytical methodologies, principally nuclear magnetic resonance spectroscopy (NMR) and mass spectrometry (MS), during the last decade have made large-scale analysis of the human metabolome a reality. This is leading to the reawakening of the importance of metabolism in human diseases, particularly cancer. The metabolome is the functional readout of the genome, functional genome, and proteome; it is also an integral partner in molecular regulations for homeostasis. The interrogation of the metabolome, or metabolomics, is now being applied to numerous diseases, largely by metabolite profiling for biomarker discovery, but also in pharmacology and therapeutics. Recent advances in stable isotope tracer-based metabolomic approaches enable unambiguous tracking of individual atoms through compartmentalized metabolic networks directly in human subjects, which promises to decipher the complexity of the human metabolome at an unprecedented pace. This knowledge will revolutionize our understanding of complex human diseases, clinical diagnostics, as well as individualized therapeutics and drug response. In this review, we focus on the use of stable isotope tracers with metabolomics technologies for understanding metabolic network dynamics in both model systems and in clinical applications. Atom-resolved isotope tracing via the two major analytical platforms, NMR and MS, has the power to determine novel metabolic reprogramming in diseases, discover new drug targets, and facilitates ADME studies. We also illustrate new metabolic tracer-based imaging technologies, which enable direct visualization of metabolic processes in vivo. We further outline current practices and future requirements for biochemoinformatics development, which is an integral part of translating stable isotope-resolved metabolomics into clinical reality.

Metabolomics methods for the synthetic biology of secondary metabolism

Many microbial secondary metabolites are of high biotechnological value for medicine, agriculture, and the food industry. Bacterial genome mining has revealed numerous novel secondary metabolite biosynthetic gene clusters, which encode the potential to synthesize a large diversity of compounds that have never been observed before. The stimulation or "awakening" of this cryptic microbial secondary metabolism has naturally attracted the attention of synthetic microbiologists, who exploit recent advances in DNA sequencing and synthesis to achieve unprecedented control over metabolic pathways. One of the indispensable tools in the synthetic biology toolbox is metabolomics, the global quantification of small biomolecules. This review illustrates the pivotal role of metabolomics for the synthetic microbiology of secondary metabolism, including its crucial role in novel compound discovery in microbes, the examination of side products of engineered metabolic pathways, as well as the identification of major bottlenecks for the overproduction of compounds of interest, especially in combination with metabolic modeling. We conclude by highlighting remaining challenges and recent technological advances that will drive metabolomics towards fulfilling its potential as a cornerstone technology of synthetic microbiology.

Collisional fragmentation of central carbon metabolites in LC-MS/MS increases precision of ¹³C metabolic flux analysis

Experimental determination of fluxes by (13)C-tracers relies on detection of (13)C-patterns in metabolites or by-products. In the field of (13)C metabolic flux analysis, the most recent developments point toward recording labeling patterns by liquid chromatography (LC)-mass spectrometry (MS)/MS directly in intermediates in central carbon metabolism (CCM) to increase temporal resolution. Surprisingly, the flux studies published so far with LC-MS measurements were based on intact metabolic intermediates-thus neglected the potential benefits of using positional information to improve flux estimation. For the first time, we exploit collisional fragmentation to obtain more fine-grained (13)C-data on intermediates of CCM and investigate their impact in (13)C metabolic flux analysis. For the case study of Bacillus subtilis grown in mineral medium with (13)C-labeled glucose, we compare the flux estimates obtained by iterative isotopologue balancing of (13)C-data obtained either by LC-MS/MS for solely intact intermediates or LC-MS/MS for intact and fragmented intermediates of CCM. We show that with LC-MS/MS data, fragment information leads to more precise estimates of fluxes in pentose phosphate pathway, glycolysis, and to the tricarboxylic acid cycle. Additionally, we present an efficient analytical strategy to rapidly acquire large sets of (13)C-patterns by tandem MS, and an in-depth analysis of the collisional fragmentation of primary intermediates. In the future, this catalogue will enable comprehensive in silico calculability analyses to identify the most sensitive measurements and direct experimental design.

Bridging the gap between fluxomics and industrial biotechnology

Metabolic flux analysis is a vital tool used to determine the ultimate output of cellular metabolism and thus detect biotechnologically relevant bottlenecks in productivity. ¹³C-based metabolic flux analysis (¹³C-MFA) and flux balance analysis (FBA) have many potential applications in biotechnology. However, noteworthy hurdles in fluxomics study are still present. First, several technical difficulties in both ¹³C-MFA and FBA severely limit the scope of fluxomics findings and the applicability of obtained metabolic information. Second, the complexity of metabolic regulation poses a great challenge for precise prediction and analysis of metabolic networks, as there are gaps between fluxomics results and other omics studies. Third, despite identified metabolic bottlenecks or sources of host stress from product synthesis, it remains difficult to overcome inherent metabolic robustness or to efficiently import and express nonnative pathways. Fourth, product yields often decrease as the number of enzymatic steps increases. Such decrease in yield may not be caused by rate-limiting enzymes, but rather is accumulated through each enzymatic reaction. Fifth, a high-throughput fluxomics tool hasnot been developed for characterizing nonmodel microorganisms and maximizing their application in industrial biotechnology. Refining fluxomics tools and understanding these obstacles will improve our ability to engineer highlyefficient metabolic pathways in microbial hosts.

Quantitative analysis of flux regulation through hierarchical regulation analysis

Regulation analysis is a methodology that quantifies to what extent a change in the flux through a metabolic pathway is regulated by either gene expression or metabolism. Two extensions to regulation analysis were developed over the past years: (i) the regulation of V(max) can be dissected into the various levels of the gene-expression cascade, such as transcription, translation, protein degradation, etc. and (ii) a time-dependent version allows following flux regulation when cells adapt to changes in their environment. The methodology of the original form of regulation analysis as well as of the two extensions will be described in detail. In addition, we will show what is needed to apply regulation analysis in practice. Studies in which the different versions of regulation analysis were applied revealed that flux regulation was distributed over various processes and depended on time, enzyme, and condition of interest. In the case of the regulation of glycolysis in baker's yeast, it appeared, however, that cells that remain under respirofermentative conditions during a physiological challenge tend to invoke more gene-expression regulation, while a shift between respirofermentative and respiratory conditions invokes an important contribution of metabolic regulation. The complexity of the regulation observed in these studies raises the question what is the advantage of this highly distributed and condition-dependent flux regulation.

13C metabolic flux analysis in complex systems

Experimental determination of in vivo metabolic rates by methods of (13)C metabolic flux analysis is a pivotal approach to unravel structure and regulation of metabolic networks, in particular with microorganisms grown in minimal media. However, the study of real-life and eukaryotic systems calls for the quantification of fluxes also in cellular compartments, rich media, cell-wide metabolic networks, dynamic systems or single cells. These scenarios drastically increase the complexity of the task, which is only partly dealt by existing approaches that rely on rigorous simulations of label propagation through metabolic networks and require multiple labeling experiments or a priori information on pathway inactivity to simplify the problem. Albeit qualitative and largely driven by human interpretation, statistical analysis of measured (13)C-patterns remains the exclusive alternative to comprehensively handle such complex systems. In the future, this practice will be complemented by novel modeling frameworks to assay particular fluxes within a network by stable isotopic tracer for targeted validation of well-defined hypotheses.