OpenFLUX: efficient modelling software for 13C-based metabolic flux analysis

Global 13C tracing and metabolic flux analysis of intact human liver tissue ex vivo

Liver metabolism is central to human physiology and influences the pathogenesis of common metabolic diseases. Yet, our understanding of human liver metabolism remains incomplete, with much of current knowledge based on animal or cell culture models that do not fully recapitulate human physiology. Here, we perform in-depth measurement of metabolism in intact human liver tissue ex vivo using global 13C tracing, non-targeted mass spectrometry and model-based metabolic flux analysis. Isotope tracing allowed qualitative assessment of a wide range of metabolic pathways within a single experiment, confirming well-known features of liver metabolism but also revealing unexpected metabolic activities such as de novo creatine synthesis and branched-chain amino acid transamination, where human liver appears to differ from rodent models. Glucose production ex vivo correlated with donor plasma glucose, suggesting that cultured liver tissue retains individual metabolic phenotypes, and could be suppressed by postprandial levels of nutrients and insulin, and also by pharmacological inhibition of glycogen utilization. Isotope tracing ex vivo allows measuring human liver metabolism with great depth and resolution in an experimentally tractable system.

Precise metabolic modeling in post-omics era: accomplishments and perspectives

Microbes have been extensively utilized for their sustainable and scalable properties in synthesizing desired bio-products. However, insufficient knowledge about intracellular metabolism has impeded further microbial applications. The genome-scale metabolic models (GEMs) play a pivotal role in facilitating a global understanding of cellular metabolic mechanisms. These models enable rational modification by exploring metabolic pathways and predicting potential targets in microorganisms, enabling precise cell regulation without experimental costs. Nonetheless, simplified GEM only considers genome information and network stoichiometry while neglecting other important bio-information, such as enzyme functions, thermodynamic properties, and kinetic parameters. Consequently, uncertainties persist particularly when predicting microbial behaviors in complex and fluctuant systems. The advent of the omics era with its massive quantification of genes, proteins, and metabolites under various conditions has led to the flourishing of multi-constrained models and updated algorithms with improved predicting power and broadened dimension. Meanwhile, machine learning (ML) has demonstrated exceptional analytical and predictive capacities when applied to training sets of biological big data. Incorporating the discriminant strength of ML with GEM facilitates mechanistic modeling efficiency and improves predictive accuracy. This paper provides an overview of research innovations in the GEM, including multi-constrained modeling, analytical approaches, and the latest applications of ML, which may contribute comprehensive knowledge toward genetic refinement, strain development, and yield enhancement for a broad range of biomolecules.

INCAWrapper: a Python wrapper for INCA for seamless data import, -export, and -processing

Motivation

INCA is a powerful tool for metabolic flux analysis, however, import and export of data and results can be tedious and limit the use of INCA in automated workflows.

Results

The INCAWrapper enables the use of INCA purely through Python, which allows the use of INCA in common data science workflows.

Availability and implementation

The INCAWrapper is implemented in Python and can be found at https://github.com/biosustain/incawrapper. It is freely available under an MIT License. To run INCA, the user needs their own MATLAB and INCA licenses. INCA is freely available for noncommercial use at mfa.vueinnovations.com.

Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference

Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.

Parameter Identification in Metabolic Reaction Networks by Means of Multiple Steady-State Measurements

In this work, we investigate some theoretical aspects related to the estimation approach proposed by Liebermeister and Klipp, 2006, in which general rate laws, derived from standardized enzymatic mechanisms, are exploited to kinetically describe the fluxes of a metabolic reaction network, and multiple metabolic steady-state measurements are exploited to estimate the unknown kinetic parameters. Further mathematical details are deeply investigated, and necessary conditions on the amount of information required to solve the identification problem are given. Moreover, theoretical results for the parameter identifiability are provided, and symmetrical and modular properties of the proposed approach are highlighted when the global identification problem is decoupled into smaller and simpler identification problems related to the single reactions of the network. Among the advantages of the proposed innovative approach are (i) non-restrictive conditions to guarantee the solvability of the parameter estimation problem, (ii) the unburden of the usual computational complexity for such identification problems, and (iii) the ease of obtaining the required number of measurements, which are actually steady-state data, experimentally easier to obtain with respect to the time-dependent ones. A simple example concludes the paper, highlighting the mentioned advantages of the method and the implementation of the related theoretical result.

H+-Translocating Membrane-Bound Pyrophosphatase from Rhodospirillum rubrum Fuels Escherichia coli Cells via an Alternative Pathway for Energy Generation

Inorganic pyrophosphatases (PPases) catalyze an essential reaction, namely, the hydrolysis of PP_i, which is formed in large quantities as a side product of numerous cellular reactions. In the majority of living species, PP_i hydrolysis is carried out by soluble cytoplasmic PPase (S-PPases) with the released energy dissipated in the form of heat. In Rhodospirillum rubrum, part of this energy can be conserved by proton-pumping pyrophosphatase (H⁺-PPase^Rru) in the form of a proton electrochemical gradient for further ATP synthesis. Here, the codon-harmonized gene hppa^Rru encoding H⁺-PPase^Rru was expressed in the Escherichia coli chromosome. We demonstrate, for the first time, that H⁺-PPase^Rru complements the essential native S-PPase in E. coli cells. ¹³C-MFA confirmed that replacing native PPase to H⁺-PPase^Rru leads to the re-distribution of carbon fluxes; a statistically significant 36% decrease in tricarboxylic acid (TCA) cycle fluxes was found compared with wild-type E. coli MG1655. Such a flux re-distribution can indicate the presence of an additional method for energy generation (e.g., ATP), which can be useful for the microbiological production of a number of compounds, the biosynthesis of which requires the consumption of ATP.

Quantitative insight into the metabolism of isoprene-producing Synechocystis sp. PCC 6803 using steady state 13C-MFA

Cyanobacteria are photosynthetic bacteria, widely studied for the conversion of atmospheric carbon dioxide to useful platform chemicals. Isoprene is one such industrially important chemical, primarily used for production of synthetic rubber and biofuels. Synechocystis sp. PCC 6803, a genetically amenable cyanobacterium, produces isoprene on heterologous expression of isoprene synthase gene, albeit in very low quantities. Rationalized metabolic engineering to re-route the carbon flux for enhanced isoprene production requires in-dept knowledge of the metabolic flux distribution in the cell. Hence, in the present study, we undertook steady state ¹³C-metabolic flux analysis of glucose-tolerant wild-type (GTN) and isoprene-producing recombinant (ISP) Synechocystis sp. to understand and compare the carbon flux distribution in the two strains. The R-values for amino acids, flux analysis predictions and gene expression profiles emphasized predominance of Calvin cycle and glycogen metabolism in GTN. Alternatively, flux analysis predicted higher activity of the anaplerotic pathway through phosphoenolpyruvate carboxylase and malic enzyme in ISP. The striking difference in the Calvin cycle, glycogen metabolism and anaplerotic pathway activity in GTN and ISP suggested a possible role of energy molecules (ATP and NADPH) in regulating the carbon flux distribution in GTN and ISP. This claim was further supported by the transcript level of selected genes of the electron transport chain. This study provides the first quantitative insight into the carbon flux distribution of isoprene-producing cyanobacterium, information critical for developing Synechocystis sp. as a single cell factory for isoprenoid/terpenoid production.

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.

CeCaFLUX: the first web server for standardized and visual instationary 13C metabolic flux analysis

SUMMARY

The number of instationary 13C-metabolic flux (INST-MFA) studies grows every year, making it more important than ever to ensure the clarity, standardization and reproducibility of each study. We proposed CeCaFLUX, the first user-friendly web server that derives metabolic flux distribution from instationary 13C-labeled data. Flux optimization and statistical analysis are achieved through an evolutionary optimization in a parallel manner. It can visualize the flux optimizing process in real time and the ultimate flux outcome. It will also function as a database to enhance the consistency and to facilitate sharing of flux studies.

AVAILABILITY AND IMPLEMENTATION

CeCaFLUX is freely available at https://www.cecaflux.net, the source code can be downloaded at https://github.com/zhzhd82/CeCaFLUX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genotypic and phenotypic diversity among Komagataella species reveals a hidden pathway for xylose utilization

BACKGROUND

The yeast genus Komagataella currently consists of seven methylotrophic species isolated from tree environments. Well-characterized strains of K. phaffii and K. pastoris are important hosts for biotechnological applications, but the potential of other species from the genus remains largely unexplored. In this study, we characterized 25 natural isolates from all seven described Komagataella species to identify interesting traits and provide a comprehensive overview of the genotypic and phenotypic diversity available within this genus.

RESULTS

Growth tests on different carbon sources and in the presence of stressors at two different temperatures allowed us to identify strains with differences in tolerance to high pH, high temperature, and growth on xylose. As Komagataella species are generally not considered xylose-utilizing yeasts, xylose assimilation was characterized in detail. Growth assays, enzyme activity measurements and 13C labeling confirmed the ability of K. phaffii to utilize D-xylose via the oxidoreductase pathway. In addition, we performed long-read whole-genome sequencing to generate genome assemblies of all Komagataella species type strains and additional K. phaffii and K. pastoris isolates for comparative analysis. All sequenced genomes have a similar size and share 83-99% average sequence identity. Genome structure analysis showed that K. pastoris and K. ulmi share the same rearrangements in difference to K. phaffii, while the genome structure of K. kurtzmanii is similar to K. phaffii. The genomes of the other, more distant species showed a larger number of structural differences. Moreover, we used the newly assembled genomes to identify putative orthologs of important xylose-related genes in the different Komagataella species.

CONCLUSIONS

By characterizing the phenotypes of 25 natural Komagataella isolates, we could identify strains with improved growth on different relevant carbon sources and stress conditions. Our data on the phenotypic and genotypic diversity will provide the basis for the use of so-far neglected Komagataella strains with interesting characteristics and the elucidation of the genetic determinants of improved growth and stress tolerance for targeted strain improvement.

GC/MS-based 13C metabolic flux analysis resolves the parallel and cyclic photomixotrophic metabolism of Synechocystis sp. PCC 6803 and selected deletion mutants including the Entner-Doudoroff and phosphoketolase pathways

BACKGROUND

Cyanobacteria receive huge interest as green catalysts. While exploiting energy from sunlight, they co-utilize sugar and CO₂. This photomixotrophic mode enables fast growth and high cell densities, opening perspectives for sustainable biomanufacturing. The model cyanobacterium Synechocystis sp. PCC 6803 possesses a complex architecture of glycolytic routes for glucose breakdown that are intertwined with the CO₂-fixing Calvin-Benson-Bassham (CBB) cycle. To date, the contribution of these pathways to photomixotrophic metabolism has remained unclear.

RESULTS

Here, we developed a comprehensive approach for ¹³C metabolic flux analysis of Synechocystis sp. PCC 6803 during steady state photomixotrophic growth. Under these conditions, the Entner-Doudoroff (ED) and phosphoketolase (PK) pathways were found inactive but the microbe used the phosphoglucoisomerase (PGI) (63.1%) and the oxidative pentose phosphate pathway (OPP) shunts (9.3%) to fuel the CBB cycle. Mutants that lacked the ED pathway, the PK pathway, or phosphofructokinases were not affected in growth under metabolic steady-state. An ED pathway-deficient mutant (Δeda) exhibited an enhanced CBB cycle flux and increased glycogen formation, while the OPP shunt was almost inactive (1.3%). Under fluctuating light, ∆eda showed a growth defect, different to wild type and the other deletion strains.

CONCLUSIONS

The developed approach, based on parallel ¹³C tracer studies with GC-MS analysis of amino acids, sugars, and sugar derivatives, optionally adding NMR data from amino acids, is valuable to study fluxes in photomixotrophic microbes to detail. In photomixotrophic cells, PGI and OPP form glycolytic shunts that merge at switch points and result in synergistic fueling of the CBB cycle for maximized CO₂ fixation. However, redirected fluxes in an ED shunt-deficient mutant and the impossibility to delete this shunt in a GAPDH2 knockout mutant, indicate that either minor fluxes (below the resolution limit of ¹³C flux analysis) might exist that could provide catalytic amounts of regulatory intermediates or alternatively, that EDA possesses additional so far unknown functions. These ideas require further experiments.

Stable Isotopes for Tracing Cardiac Metabolism in Diseases

Although metabolic remodeling during cardiovascular diseases has been well-recognized for decades, the recent development of analytical platforms and mathematical tools has driven the emergence of assessing cardiac metabolism using tracers. Metabolism is a critical component of cellular functions and adaptation to stress. The pathogenesis of cardiovascular disease involves metabolic adaptation to maintain cardiac contractile function even in advanced disease stages. Stable-isotope tracer measurements are a powerful tool for measuring flux distributions at the whole organism level and assessing metabolic changes at a systems level in vivo. The goal of this review is to summarize techniques and concepts for in vivo or ex vivo stable isotope labeling in cardiovascular research, to highlight mathematical concepts and their limitations, to describe analytical methods at the tissue and single-cell level, and to discuss opportunities to leverage metabolic models to address important mechanistic questions relevant to all patients with cardiovascular disease.

Recent advances in metabolic engineering-integration of in silico design and experimental analysis of metabolic pathways

Microorganisms are widely used to produce valuable compounds. Because thousands of metabolic reactions occur simultaneously and many metabolic reactions are related to target production and cell growth, the development of a rational design method for metabolic pathway modification to optimize target production is needed. In this paper, recent advances in metabolic engineering are reviewed, specifically considering computational pathway modification design and experimental evaluation of metabolic fluxes by ¹³C-metabolic flux analysis. Computational tools for seeking effective gene deletion targets and recruiting heterologous genes are described in flux balance analysis approaches. A kinetic model and adaptive laboratory evolution were applied to identify and eliminate the rate-limiting step in metabolic pathways. Data science-based approaches for process monitoring and control are described to maximize the performance of engineered cells in bioreactors.

The Metabolic Flux Probe (MFP)-Secreted Protein as a Non-Disruptive Information Carrier for 13C-Based Metabolic Flux Analysis

Novel cultivation technologies demand the adaptation of existing analytical concepts. Metabolic flux analysis (MFA) requires stable-isotope labeling of biomass-bound protein as the primary information source. Obtaining the required protein in cultivation set-ups where biomass is inaccessible due to low cell densities and cell immobilization is difficult to date. We developed a non-disruptive analytical concept for ¹³C-based metabolic flux analysis based on secreted protein as an information carrier for isotope mapping in the protein-bound amino acids. This “metabolic flux probe” (MFP) concept was investigated in different cultivation set-ups with a recombinant, protein-secreting yeast strain. The obtained results grant insight into intracellular protein turnover dynamics. Experiments under metabolic but isotopically nonstationary conditions in continuous glucose-limited chemostats at high dilution rates demonstrated faster incorporation of isotope information from labeled glucose into the recombinant reporter protein than in biomass-bound protein. Our results suggest that the reporter protein was polymerized from intracellular amino acid pools with higher turnover rates than biomass-bound protein. The latter aspect might be vital for ¹³C-flux analyses under isotopically nonstationary conditions for analyzing fast metabolic dynamics.

13C-metabolic flux analysis in glycerol-assimilating strains of Saccharomyces cerevisiae

Glycerol is an attractive raw material for the production of useful chemicals using microbial cells. We previously identified metabolic engineering targets for the improvement of glycerol assimilation ability in Saccharomyces cerevisiae based on adaptive laboratory evolution (ALE) and transcriptome analysis of the evolved cells. We also successfully improved glycerol assimilation ability by the disruption of the RIM15 gene encoding a Greatwall protein kinase together with overexpression of the STL1 gene encoding the glycerol/H⁺ symporter. To understand glycerol assimilation metabolism in the evolved glycerol-assimilating strains and STL1-overexpressing RIM15 disruptant, we performed metabolic flux analysis using ¹³C-labeled glycerol. Significant differences in metabolic flux distributions between the strains obtained from the culture after 35 and 85 generations in ALE were not found, indicating that metabolic flux changes might occur in the early phase of ALE (i.e., before 35 generations at least). Similarly, metabolic flux distribution was not significantly changed by RIM15 gene disruption. However, fluxes for the lower part of glycolysis and the TCA cycle were larger and, as a result, flux for the pentose phosphate pathway was smaller in the STL1-overexpressing RIM15 disruptant than in the strain obtained from the culture after 85 generations in ALE. It could be effective to increase flux for the pentose phosphate pathway to improve the glycerol assimilation ability in S. cerevisiae.

Investigating metabolism by mass spectrometry: From steady state to dynamic view

Metabolism is the set of life-sustaining reactions in organisms. These biochemical reactions are organized in metabolic pathways, in which one metabolite is converted through a series of steps catalyzed by enzymes in another chemical compound. Metabolic reactions are categorized as catabolic, the breaking down of metabolites to produce energy, and/or anabolic, the synthesis of compounds that consume energy. The balance between catabolism of the preferential fuel substrate and anabolism defines the overall metabolism of a cell or tissue. Metabolomics is a powerful tool to gain new insights contributing to the identification of complex molecular mechanisms in the field of biomedical research, both basic and translational. The enormous potential of this kind of analyses consists of two key aspects: (i) the possibility of performing so-called targeted and untargeted experiments through which it is feasible to verify or formulate a hypothesis, respectively, and (ii) the opportunity to run either steady-state analyses to have snapshots of the metabolome at a given time under different experimental conditions or dynamic analyses through the use of labeled tracers. In this review, we will highlight the most important practical (e.g., different sample extraction approaches) and conceptual steps to consider for metabolomic analysis, describing also the main application contexts in which it is used. In addition, we will provide some insights into the most innovative approaches and progress in the field of data analysis and processing, highlighting how this part is essential for the proper extrapolation and interpretation of data.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

Metabolic Flux Analysis-Linking Isotope Labeling and Metabolic Fluxes

Metabolic flux analysis (MFA) is an increasingly important tool to study metabolism quantitatively. Unlike the concentrations of metabolites, the fluxes, which are the rates at which intracellular metabolites interconvert, are not directly measurable. MFA uses stable isotope labeled tracers to reveal information related to the fluxes. The conceptual idea of MFA is that in tracer experiments the isotope labeling patterns of intracellular metabolites are determined by the fluxes, therefore by measuring the labeling patterns we can infer the fluxes in the network. In this review, we will discuss the basic concept of MFA using a simplified upper glycolysis network as an example. We will show how the fluxes are reflected in the isotope labeling patterns. The central idea we wish to deliver is that under metabolic and isotopic steady-state the labeling pattern of a metabolite is the flux-weighted average of the substrates’ labeling patterns. As a result, MFA can tell the relative contributions of converging metabolic pathways only when these pathways make substrates in different labeling patterns for the shared product. This is the fundamental principle guiding the design of isotope labeling experiment for MFA including tracer selection. In addition, we will also discuss the basic biochemical assumptions of MFA, and we will show the flux-solving procedure and result evaluation. Finally, we will highlight the link between isotopically stationary and nonstationary flux analysis.

Enzymes as Parts in Need of Replacement - and How to Extend Their Working Life

Enzymes catalyze reactions in vivo at different rates and each enzyme molecule has a lifetime limit before it is degraded and replaced to enable catalysis to continue. Considering these rates together as a unitless ratio of catalytic cycles until replacement (CCR) provides a new quantitative tool to assess the replacement schedule of and energy investment into enzymes as they relate to function. Here, we outline the challenges of determining CCRs and new approaches to overcome them and then assess the CCRs of selected enzymes in bacteria and plants to reveal a range of seven orders of magnitude for this ratio. Modifying CCRs in plants holds promise to lower cellular costs, to tailor enzymes for particular environments, and to breed enzyme improvements for crop productivity.

Parallel isotope differential modeling for instationary 13C fluxomics at the genome scale

BACKGROUND

A precise map of the metabolic fluxome, the closest surrogate to the physiological phenotype, is becoming progressively more important in the metabolic engineering of photosynthetic organisms for biofuel and biomass production. For photosynthetic organisms, the state-of-the-art method for this purpose is instationary 13C fluxomics, which has arisen as a sibling of transcriptomics or proteomics. Instationary 13C data processing requires solving high-dimensional nonlinear differential equations and leads to large computational and time costs when its scope is expanded to a genome-scale metabolic network.

RESULT

Here, we present a parallelized method to model instationary 13C labeling data. The elementary metabolite unit (EMU) framework is reorganized to allow treating individual mass isotopomers and breaking up of their networks into strongly connected components (SCCs). A variable domain parallel algorithm is introduced to process ordinary differential equations in a parallel way. 15-fold acceleration is achieved for constant-step-size modeling and ~ fivefold acceleration for adaptive-step-size modeling.

CONCLUSION

This algorithm is universally applicable to isotope granules such as EMUs and cumomers and can substantially accelerate instationary 13C fluxomics modeling. It thus has great potential to be widely adopted in any instationary 13C fluxomics modeling.

Dynamic 13C Flux Analysis Captures the Reorganization of Adipocyte Glucose Metabolism in Response to Insulin

Cellular metabolism is dynamic, but quantifying non-steady metabolic fluxes by stable isotope tracers presents unique computational challenges. Here, we developed an efficient 13C-tracer dynamic metabolic flux analysis (13C-DMFA) framework for modeling central carbon fluxes that vary over time. We used B-splines to generalize the flux parameterization system and to improve the stability of the optimization algorithm. As proof of concept, we investigated how 3T3-L1 cultured adipocytes acutely metabolize glucose in response to insulin. Insulin rapidly stimulates glucose uptake, but intracellular pathways responded with differing speeds and magnitudes. Fluxes in lower glycolysis increased faster than those in upper glycolysis. Glycolysis fluxes rose disproportionally larger and faster than the tricarboxylic acid cycle, with lactate a primary glucose end product. The uncovered array of flux dynamics suggests that glucose catabolism is additionally regulated beyond uptake to help shunt glucose into appropriate pathways. This work demonstrates the value of using dynamic intracellular fluxes to understand metabolic function and pathway regulation.

Metabolic Rearrangements Causing Elevated Proline and Polyhydroxybutyrate Accumulation During the Osmotic Adaptation Response of Bacillus megaterium

For many years now, Bacillus megaterium serves as a microbial workhorse for the high-level production of recombinant proteins in the g/L-scale. However, efficient and stable production processes require the knowledge of the molecular adaptation strategies of the host organism to establish optimal environmental conditions. Here, we interrogated the osmotic stress response of B. megaterium using transcriptome, proteome, metabolome, and fluxome analyses. An initial transient adaptation consisted of potassium import and glutamate counterion synthesis. The massive synthesis of the compatible solute proline constituted the second longterm adaptation process. Several stress response enzymes involved in iron scavenging and reactive oxygen species (ROS) fighting proteins showed higher levels under prolonged osmotic stress induced by 1.8 M NaCl. At the same time, the downregulation of the expression of genes of the upper part of glycolysis resulted in the activation of the pentose phosphate pathway (PPP), generating an oversupply of NADPH. The increased production of lactate accompanied by the reduction of acetate secretion partially compensate for the unbalanced (NADH/NAD⁺) ratio. Besides, the tricarboxylic acid cycle (TCA) mainly supplies the produced NADH, as indicated by the higher mRNA and protein levels of involved enzymes, and further confirmed by ¹³C flux analyses. As a consequence of the metabolic flux toward acetyl-CoA and the generation of an excess of NADPH, B. megaterium redirected the produced acetyl-CoA toward the polyhydroxybutyrate (PHB) biosynthetic pathway accumulating around 30% of the cell dry weight (CDW) as PHB. This direct relation between osmotic stress and intracellular PHB content has been evidenced for the first time, thus opening new avenues for synthesizing this valuable biopolymer using varying salt concentrations under non-limiting nutrient conditions.

Impact framework: A python package for writing data analysis workflows to interpret microbial physiology

Microorganisms can be genetically engineered to solve a range of challenges in diverse including health, environmental protection and sustainability. The natural complexity of biological systems makes this an iterative cycle, perturbing metabolism and making stepwise progress toward a desired phenotype through four major stages: design, build, test, and data interpretation. This cycle has been accelerated by advances in molecular biology (e.g. robust DNA synthesis and assembly techniques), liquid handling automation and scale-down characterization platforms, generating large heterogeneous data sets. Here, we present an extensible Python package for scientists and engineers working with large biological data sets to interpret, model, and visualize data: the IMPACT (Integrated Microbial Physiology: Analysis, Characterization and Translation) framework. Impact aims to ease the development of Python-based data analysis workflows for a range of stakeholders in the bioengineering process, offering open-source tools for data analysis, physiology characterization and translation to visualization. Using this framework, biologists and engineers can opt for reproducible and extensible programmatic data analysis workflows, mediating a bottleneck limiting the throughput of microbial engineering. The Impact framework is available at https://github.com/lmse/impact.

Principal metabolic flux mode analysis

Motivation

In the analysis of metabolism, two distinct and complementary approaches are frequently used: Principal component analysis (PCA) and stoichiometric flux analysis. PCA is able to capture the main modes of variability in a set of experiments and does not make many prior assumptions about the data, but does not inherently take into account the flux mode structure of metabolism. Stoichiometric flux analysis methods, such as Flux Balance Analysis (FBA) and Elementary Mode Analysis, on the other hand, are able to capture the metabolic flux modes, however, they are primarily designed for the analysis of single samples at a time, and not best suited for exploratory analysis on a large sets of samples.

Results

We propose a new methodology for the analysis of metabolism, called Principal Metabolic Flux Mode Analysis (PMFA), which marries the PCA and stoichiometric flux analysis approaches in an elegant regularized optimization framework. In short, the method incorporates a variance maximization objective form PCA coupled with a stoichiometric regularizer, which penalizes projections that are far from any flux modes of the network. For interpretability, we also introduce a sparse variant of PMFA that favours flux modes that contain a small number of reactions. Our experiments demonstrate the versatility and capabilities of our methodology. The proposed method can be applied to genome-scale metabolic network in efficient way as PMFA does not enumerate elementary modes. In addition, the method is more robust on out-of-steady steady-state experimental data than competing flux mode analysis approaches.

Availability and implementation

Matlab software for PMFA and SPMFA and dataset used for experiments are available in https://github.com/aalto-ics-kepaco/PMFA.

Supplementary information

Supplementary data are available at Bioinformatics online.

High extracellular lactate causes reductive carboxylation in breast tissue cell lines grown under normoxic conditions

In cancer tumors, lactate accumulation was initially attributed to high glucose consumption associated with the Warburg Effect. Now it is evident that lactate can also serve as an energy source in cancer cell metabolism. Additionally, lactate has been shown to promote metastasis, generate gene expression patterns in cancer cells consistent with "cancer stem cell" phenotypes, and result in treatment resistant tumors. Therefore, the goal of this work was to quantify the impact of lactate on metabolism in three breast cell lines (one normal and two breast cancer cell lines-MCF 10A, MCF7, and MDA-MB-231), in order to better understand the role lactate may have in different disease cell types. Parallel labeling metabolic flux analysis (13C-MFA) was used to quantify the intracellular fluxes under normal and high extracellular lactate culture conditions. Additionally, high extracellular lactate cultures were labelled in parallel with [U-13C] lactate, which provided qualitative information regarding the lactate uptake and metabolism. The 13C-MFA model, which incorporated the measured extracellular fluxes and the parallel labeling mass isotopomer distributions (MIDs) for five glycolysis, four tricarboxylic acid cycle (TCA), and three intracellular amino acid metabolites, predicted lower glycolysis fluxes in the high lactate cultures. All three cell lines experienced reductive carboxylation of glutamine to citrate in the TCA cycle as a result of high extracellular lactate. Reductive carboxylation previously has been observed under hypoxia and other mitochondrial stresses, whereas these cultures were grown aerobically. In addition, this is the first study to investigate the intracellular metabolic responses of different stages of breast cancer progression to high lactate exposure. These results provide insight into the role lactate accumulation has on metabolic reaction distributions in the different disease cell types while the cells are still proliferating in lactate concentrations that do not significantly decrease exponential growth rates.

Snail-Overexpression Induces Epithelial-mesenchymal Transition and Metabolic Reprogramming in Human Pancreatic Ductal Adenocarcinoma and Non-tumorigenic Ductal Cells

The zinc finger transcription factor Snail is a known effector of epithelial-to-mesenchymal transition (EMT), a process that underlies the enhanced invasiveness and chemoresistance of common to cancerous cells. Induction of Snail-driven EMT has also been shown to drive a range of pro-survival metabolic adaptations in different cancers. In the present study, we sought to determine the specific role that Snail has in driving EMT and adaptive metabolic programming in pancreatic ductal adenocarcinoma (PDAC) by overexpressing Snail in a PDAC cell line, Panc1, and in immortalized, non-tumorigenic human pancreatic ductal epithelial (HPDE) cells. Snail overexpression was able to induce EMT in both pancreatic cell lines through suppression of epithelial markers and upregulation of mesenchymal markers alongside changes in cell morphology and enhanced migratory capacity. Snail-overexpressed pancreatic cells additionally displayed increased glucose uptake and lactate production with concomitant reduction in oxidative metabolism measurements. Snail overexpression reduced maximal respiration in both Panc1 and HPDE cells, with further reductions seen in ATP production, spare respiratory capacity and non-mitochondrial respiration in Snail overexpressing Panc1 cells. Accordingly, lower expression of mitochondrial electron transport chain proteins was observed with Snail overexpression, particularly within Panc1 cells. Modelling of ¹³C metabolite flux within both cell lines revealed decreased carbon flux from glucose in the TCA cycle in snai1-overexpressing Panc1 cells only. This work further highlights the role that Snail plays in EMT and demonstrates its specific effects on metabolic reprogramming of glucose metabolism in PDAC.

The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) is the method of choice when a detailed inference of intracellular metabolic fluxes in living organisms under metabolic quasi-steady state conditions is desired. Being continuously developed since two decades, the technology made major contributions to the quantitative characterization of organisms in all fields of biotechnology and health-related research. ¹³C MFA, however, stands out from other "-omics sciences," in that it requires not only experimental-analytical data, but also mathematical models and a computational toolset to infer the quantities of interest, i.e., the metabolic fluxes. At present, these models cannot be conveniently exchanged between different labs. Here, we present the implementation-independent model description language FluxML for specifying ¹³C MFA models. The core of FluxML captures the metabolic reaction network together with atom mappings, constraints on the model parameters, and the wealth of data configurations. In particular, we describe the governing design processes that shaped the FluxML language. We demonstrate the utility of FluxML to represent many contemporary experimental-analytical requirements in the field of ¹³C MFA. The major aim of FluxML is to offer a sound, open, and future-proof language to unambiguously express and conserve all the necessary information for model re-use, exchange, and comparison. Along with FluxML, several powerful computational tools are supplied for easy handling, but also to maintain a maximum of flexibility. Altogether, the FluxML collection is an "all-around carefree package" for ¹³C MFA modelers. We believe that FluxML improves scientific productivity as well as transparency and therewith contributes to the efficiency and reproducibility of computational modeling efforts in the field of ¹³C MFA.

EMUlator: An Elementary Metabolite Unit (EMU) Based Isotope Simulator Enabled by Adjacency Matrix

Stable isotope based metabolic flux analysis is currently the unique methodology that allows the experimental study of the integrated responses of metabolic networks. This method primarily relies on isotope labeling and modeling, which could be a challenge in both experimental and computational biology. In particular, the algorithm implementation for isotope simulation is a critical step, limiting extensive usage of this powerful approach. Here, we introduce EMUlator a Python-based isotope simulator which is developed on Elementary Metabolite Unit (EMU) algorithm, an efficient and powerful algorithm for isotope modeling. We propose a novel adjacency matrix method to implement EMU modeling and exemplify it stepwise. This method is intuitively straightforward and can be conveniently mastered for various customized purposes. We apply this arithmetic pipeline to understand the phosphoketolase flux in the metabolic network of an industrial microbe Clostridium acetobutylicum. The resulting design enables a high-throughput and non-invasive approach for estimating phosphoketolase flux in vivo. Our computational insights allow the systematic design and prediction of isotope-based metabolic models and yield a comprehensive understanding of their limitations and potentials.

Elevated temperatures do not trigger a conserved metabolic network response among thermotolerant yeasts

BACKGROUND

Thermotolerance is a highly desirable trait of microbial cell factories and has been the focus of extensive research. Yeast usually tolerate only a narrow temperature range and just two species, Kluyveromyces marxianus and Ogataea polymorpha have been described to grow at reasonable rates above 40 °C. However, the complex mechanisms of thermotolerance in yeast impede its full comprehension and the rare physiological data at elevated temperatures has so far not been matched with corresponding metabolic analyses.

RESULTS

To elaborate on the metabolic network response to increased fermentation temperatures of up to 49 °C, comprehensive physiological datasets of several Kluyveromyces and Ogataea strains were generated and used for 13C-metabolic flux analyses. While the maximum growth temperature was very similar in all investigated strains, the metabolic network response to elevated temperatures was not conserved among the different species. In fact, metabolic flux distributions were remarkably irresponsive to increasing temperatures in O. polymorpha, while the K. marxianus strains exhibited extensive flux rerouting at elevated temperatures.

CONCLUSIONS

While a clear mechanism of thermotolerance is not deducible from the fluxome level alone, the generated data can be valued as a knowledge repository for using temperature to modulate the metabolic activity towards engineering goals.

Mapping Salmonella typhimurium pathways using 13C metabolic flux analysis

In the last years, Salmonella has been extensively studied not only due to its importance as a pathogen, but also as a host to produce pharmaceutical compounds. However, the full exploitation of Salmonella as a platform for bioproduct delivery has been hampered by the lack of information about its metabolism. Genome-scale metabolic models can be valuable tools to delineate metabolic engineering strategies as long as they closely represent the actual metabolism of the target organism. In the present study, a ¹³C-MFA approach was applied to map the fluxes at the central carbon pathways of S. typhimurium LT2 growing at glucose-limited chemostat cultures. The experiments were carried out in a 2L bioreactor, using defined medium enriched with 20% ¹³C-labeled glucose. Metabolic flux distributions in central carbon pathways of S. typhimurium LT2 were estimated using OpenFLUX2 based on the labeling pattern of biomass protein hydrolysates together with biomass composition. The results suggested that pentose phosphate is used to catabolize glucose, with minor fluxes through glycolysis. In silico simulations, using Optflux and pFBA as simulation method, allowed to study the performance of the genome-scale metabolic model. In general, the accuracy of in silico simulations was improved by the superimposition of estimated intracellular fluxes to the existing genome-scale metabolic model, showing a better fitting to the experimental extracellular fluxes, whereas the intracellular fluxes of pentose phosphate and anaplerotic reactions were poorly described.

Tandem Mass Spectrometry for 13C Metabolic Flux Analysis: Methods and Algorithms Based on EMU Framework

In the past two decades, ¹³C metabolic flux analysis (¹³C-MFA) has matured into a powerful and widely used scientific tool in metabolic engineering and systems biology. Traditionally, metabolic fluxes have been determined from measurements of isotopic labeling by means of mass spectrometry (MS) or nuclear magnetic resonance (NMR). In recent years, tandem MS has emerged as a new analytical technique that can provide additional information for high-resolution quantification of metabolic fluxes in complex biological systems. In this paper, we present recent advances in methods and algorithms for incorporating tandem MS measurements into existing ¹³C-MFA approaches that are based on the elementary metabolite units (EMU) framework. Specifically, efficient EMU-based algorithms are presented for simulating tandem MS data, tracing isotopic labeling in biochemical network models and for correcting tandem MS data for natural isotope abundances.

13C-Fingerprinting and Metabolic Flux Analysis of Bacterial Metabolisms

¹³C-assisted metabolism analysis provides rigorous calculations of the intracellular reaction rates (i.e., fluxes) within the central metabolism of microbial hosts. This mapping of the intracellular network within microbes has proven to be essential for understanding the cell physiology. The approach is also a key to identifying central metabolic nodes, probing the rigidity of a metabolic network, revealing cofactor balances, and delineating hidden pathways. Here we present the methodology of using stable isotopic carbon substrates for both qualitative (¹³C-fingerprinting of functional pathways) and quantitative (Metabolic Flux Analysis) metabolism studies on bacterial species. In this methodology, we include step-by-step instructions to use the open source WUflux software for the steady-state flux calculations based on labeling information of amino acids or free metabolites.

A single Gal4-like transcription factor activates the Crabtree effect in Komagataella phaffii

The Crabtree phenotype defines whether a yeast can perform simultaneous respiration and fermentation under aerobic conditions at high growth rates. It provides Crabtree positive yeasts an evolutionary advantage of consuming glucose faster and producing ethanol to outcompete other microorganisms in sugar rich environments. While a number of genetic events are associated with the emergence of the Crabtree effect, its evolution remains unresolved. Here we show that overexpression of a single Gal4-like transcription factor is sufficient to convert Crabtree-negative Komagataella phaffii (Pichia pastoris) into a Crabtree positive yeast. Upregulation of the glycolytic genes and a significant increase in glucose uptake rate due to the overexpression of the Gal4-like transcription factor leads to an overflow metabolism, triggering both short-term and long-term Crabtree phenotypes. This indicates that a single genetic perturbation leading to overexpression of one gene may have been sufficient as the first molecular event towards respiro-fermentative metabolism in the course of yeast evolution.

Aerobic ethanol production, a phenomenon referred as Crabtree effect, allows yeast to outcompete other microorganisms in sugar rich environments. Here, the authors show that overexpression of a Gal4-like transcription factor can transform Komagataella phaffii from Crabtree effect negative to positive.

Analytical Considerations of Stable Isotope Labelling in Lipidomics

Over the last two decades, lipids have come to be understood as far more than merely components of cellular membranes and forms of energy storage, and are now also being implicated to play important roles in a variety of diseases, with lipid biomarker research one of the most widespread applications of lipidomic techniques both in research and in clinical settings. Stable isotope labelling has become a staple technique in the analysis of small molecule metabolism and dynamics, as it is the only experimental setup by which biosynthesis, remodelling and degradation of biomolecules can be directly measured. Using state-of-the-art analytical technologies such as chromatography-coupled high resolution tandem mass spectrometry, the stable isotope label can be precisely localized and quantified within the biomolecules. The application of stable isotope labelling to lipidomics is however complicated by the diversity of lipids and the complexity of the necessary data analysis. This article discusses key experimental aspects of stable isotope labelling in the field of mass spectrometry-based lipidomics, summarizes current applications and provides an outlook on future developments and potential.

13C Metabolic Flux Analysis of Enhanced Lipid Accumulation Modulated by Ethanolamine in Crypthecodinium cohnii

The heterotrophic microalga Crypthecodinium cohnii has attracted considerable attention due to its capability of accumulating lipids with a high fraction of docosahexaenoic acid (DHA). In our previous study, ethanolamine (ETA) was identified as an effective chemical modulator for lipid accumulation in C. cohnii. In this study, to gain a better understanding of the lipid metabolism and mechanism for the positive effects of modulator ETA, metabolic flux analysis was performed using ¹³C-labeled glucose with and without 1 mM ETA modulator. The analysis of flux distribution showed that with the addition of ETA, flux in glycolysis pathway and citrate pyruvate cycle was strengthened while flux in pentose phosphate pathway was decreased. In addition, flux in TCA cycle was slightly decreased compared with the control without ETA. The enzyme activity of malic enzyme (ME) was significantly increased, suggesting that NADP⁺-dependent ME might be the major source of NADPH for lipid accumulation. The flux information obtained by this study could be valuable for the further efforts in improving lipid accumulation and DHA production in C. cohnii.

FluxPyt: a Python-based free and open-source software for 13C-metabolic flux analyses

¹³C-Metabolic flux analysis (MFA) is a powerful approach to estimate intracellular reaction rates which could be used in strain analysis and design. Processing and analysis of labeling data for calculation of fluxes and associated statistics is an essential part of MFA. However, various software currently available for data analysis employ proprietary platforms and thus limit accessibility. We developed FluxPyt, a Python-based truly open-source software package for conducting stationary ¹³C-MFA data analysis. The software is based on the efficient elementary metabolite unit framework. The standard deviations in the calculated fluxes are estimated using the Monte-Carlo analysis. FluxPyt also automatically creates flux maps based on a template for visualization of the MFA results. The flux distributions calculated by FluxPyt for two separate models: a small tricarboxylic acid cycle model and a larger Corynebacterium glutamicum model, were found to be in good agreement with those calculated by a previously published software. FluxPyt was tested in Microsoft™ Windows 7 and 10, as well as in Linux Mint 18.2. The availability of a free and open ¹³C-MFA software that works in various operating systems will enable more researchers to perform ¹³C-MFA and to further modify and develop the package.