Segregated flux balance analysis constrained by population structure/function data: The case of PHA production by mixed microbial cultures

In this study we developed a segregated flux balance analysis (FBA) method to calculate metabolic flux distributions of the individual populations present in a mixed microbial culture (MMC). Population specific flux data constraints were derived from the raw data typically obtained by the fluorescence in situ hybridization (FISH) and microautoradiography (MAR)‐FISH techniques. This method was applied to study the metabolic heterogeneity of a MMC that produces polyhydroxyalkanoates (PHA) from fermented sugar cane molasses. Three populations were identified by FISH, namely Paracoccus sp., Thauera sp., and Azoarcus sp. The segregated FBA method predicts a flux distribution for each of the identified populations. The method is shown to predict with high accuracy the average PHA storage flux and the respective monomeric composition for 16 independent experiments. Moreover, flux predictions by segregated FBA were slightly better than those obtained by nonsegregated FBA, and also highly concordant with metabolic flux analysis (MFA) estimated fluxes. The segregated FBA method can be of high value to assess metabolic heterogeneity in MMC systems and to derive more efficient eco‐engineering strategies. For the case of PHA‐producing MMC considered in this work, it becomes apparent that the PHA average monomeric composition might be controlled not only by the volatile fatty acids (VFA) feeding profile but also by the population composition present in the MMC. Biotechnol. Bioeng. 2013; 110: 2267–2276. © 2013 Wiley Periodicals, Inc.     

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and ¹³C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.     

Large‐scale identification of genetic design strategies using local search

In the past decade, computational methods have been shown to be well suited to unraveling the complex web of metabolic reactions in biological systems. Methods based on flux–balance analysis (FBA) and bi‐level optimization have been used to great effect in aiding metabolic engineering. These methods predict the result of genetic manipulations and allow for the best set of manipulations to be found computationally. Bi‐level FBA is, however, limited in applicability because the required computational time and resources scale poorly as the size of the metabolic system and the number of genetic manipulations increase. To overcome these limitations, we have developed Genetic Design through Local Search (GDLS), a scalable, heuristic, algorithmic method that employs an approach based on local search with multiple search paths, which results in effective, low‐complexity search of the space of genetic manipulations. Thus, GDLS is able to find genetic designs with greater in silico production of desired metabolites than can feasibly be found using a globally optimal search and performs favorably in comparison with heuristic searches based on evolutionary algorithms and simulated annealing.     

Improving the accuracy of flux balance analysis through the implementation of carbon availability constraints for intracellular reactions

Constraint-based modeling methods, such as Flux Balance Analysis (FBA), have been extensively used to decipher complex, information rich -omics datasets to elicit system-wide behavioral patterns of cellular metabolism. FBA has been successfully used to gain insight in a wide range of applications, such as range of substrate utilization, product yields and to design metabolic engineering strategies to improve bioprocess performance. A well-known challenge associated with large genome-scale metabolic networks is that they result in underdetermined problem formulations. Consequently, rather than unique solutions, FBA and related methods examine ranges of reaction flux values that are consistent with the studied physiological conditions. The wider the reported flux ranges, the higher the uncertainty in the determination of basic reaction properties, limiting interpretability of and confidence in the results. Herein, we propose a new, computationally efficient approach that refines flux range predictions by constraining reaction fluxes on the basis of the elemental balance of carbon. We compared carbon constraint FBA (ccFBA) against experimentally-measured intracellular fluxes using the latest CHO GEM (iCHO1766) and were able to substantially improve the accuracy of predicted flux values compared with FBA. ccFBA can be used as a stand-alone method but is also compatible with and complimentary to other constraint-based approaches.     

ΔFBA—Predicting metabolic flux alterations using genome-scale metabolic models and differential transcriptomic data

Genome-scale metabolic models (GEMs) provide a powerful framework for simulating the entire set of biochemical reactions in a cell using a constraint-based modeling strategy called flux balance analysis (FBA). FBA relies on an assumed metabolic objective for generating metabolic fluxes using GEMs. But, the most appropriate metabolic objective is not always obvious for a given condition and is likely context-specific, which often complicate the estimation of metabolic flux alterations between conditions. Here, we propose a new method, called ΔFBA (deltaFBA), that integrates differential gene expression data to evaluate directly metabolic flux differences between two conditions. Notably, ΔFBA does not require specifying the cellular objective. Rather, ΔFBA seeks to maximize the consistency and minimize inconsistency between the predicted flux differences and differential gene expression. We showcased the performance of ΔFBA through several case studies involving the prediction of metabolic alterations caused by genetic and environmental perturbations in Escherichia coli and caused by Type-2 diabetes in human muscle. Importantly, in comparison to existing methods, ΔFBA gives a more accurate prediction of flux differences.     

A guide to metabolic flux analysis in metabolic engineering: Methods,tools and applications

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, “Metabolic fluxes and metabolic engineering” (Metabolic Engineering, 1: 1–11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. ¹³C-metabolic flux analysis (¹³C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.     

A generalized computational framework to streamline thermodynamics and kinetics analysis of metabolic pathways

Metabolic engineering is a critical biotechnological approach in addressing global energy and environment challenges. Most engineering efforts, however, consist of laborious and inefficient trial-and-error of target pathways, due in part to the lack of methodologies that can comprehensively assess pathway properties in thermodynamics and kinetics. Metabolic engineering can benefit from computational tools that evaluate feasibility, expense and stability of non-natural metabolic pathways. Such tools can also help us understand natural pathways and their regulation at systems level. Here we introduce a computational toolbox, PathParser, which, for the first time, integrates multiple important functions for pathway analysis including thermodynamics analysis, kinetics-based protein cost optimization and robustness analysis. Specifically, PathParser enables optimization of the driving force of a pathway by minimizing the Gibbs free energy of least thermodynamically favorable reaction. In addition, based on reaction thermodynamics and enzyme kinetics, it can compute the minimal enzyme protein cost that supports metabolic flux, and evaluate pathway stability and flux in response to enzyme concentration perturbations. In a demo analysis of the Calvin–Benson–Bassham cycle and photorespiration pathway in the model cyanobacterium Synechocystis PCC 6803, the computation results are corroborated by experimental proteomics data as well as metabolic engineering outcomes. This toolbox may have broad application in metabolic engineering and systems biology in other microbial systems.     

Energy balance for analysis of complex metabolic networks

Predicting behavior of large-scale biochemical networks represents one of the greatest challenges of bioinformatics and computational biology. Computational tools for predicting fluxes in biochemical networks are applied in the fields of integrated and systems biology, bioinformatics, and genomics, and to aid in drug discovery and identification of potential drug targets. Approaches, such as flux balance analysis (FBA), that account for the known stoichiometry of the reaction network while avoiding implementation of detailed reaction kinetics are promising tools for the analysis of large complex networks. Here we introduce energy balance analysis (EBA)--the theory and methodology for enforcing the laws of thermodynamics in such simulations--making the results more physically realistic and revealing greater insight into the regulatory and control mechanisms operating in complex large-scale systems. We show that EBA eliminates thermodynamically infeasible results associated with FBA.     

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol⁻¹, and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.     

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures—Systems biology‐based interpretation using genome‐scale metabolic flux balance model and multivariate data analysis

Genome‐scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745–753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome‐scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163–1173, 2016     

Development and experimental verification of a genome-scale metabolic model for Corynebacterium glutamicum

Background  

In silico genome-scale metabolic models enable the analysis of the characteristics of metabolic systems of organisms. In this study, we reconstructed a genome-scale metabolic model of Corynebacterium glutamicum on the basis of genome sequence annotation and physiological data. The metabolic characteristics were analyzed using flux balance analysis (FBA), and the results of FBA were validated using data from culture experiments performed at different oxygen uptake rates.     

Error propagation in constraint-based modeling of Chinese hamster ovary cells

Chinese hamster ovary (CHO) cells are the most popular mammalian cell factories for the production of glycosylated biopharmaceuticals. To further increase titer and productivity and ensure product quality, rational system-level engineering strategies based on constraint-based metabolic modeling, such as flux balance analysis (FBA), have gained strong interest. However, the quality of FBA predictions depends on the accuracy of the experimental input data, especially on the exchange rates of extracellular metabolites. Yet, it is not standard practice to devote sufficient attention to the accurate determination of these rates. In this work, we investigated to what degree the sampling frequency during a batch culture and the measurement errors of metabolite concentrations influence the accuracy of the calculated exchange rates and further, how this error then propagates into FBA predictions of growth rates. We determined that accurate measurements of essential amino acids with low uptake rates are crucial for the accuracy of FBA predictions, followed by a sufficient number of analyzed time points. We observed that the measured difference in growth rates of two cell lines can only be reliably predicted when both high measurement accuracy and sampling frequency are ensured.     

Urinary C-peptide is not an accurate bioindicator of energy balance in humans

The apprehension of the factors that affect long term regulation of energy balance is indispensable to understand the rise in obesity prevalence as well as to delineate levers to prevent it. Accurate measurements of energy balance are however challenging during free‐living conditions. Recent studies proposed urinary C‐peptide, a metabolic byproduct of insulin synthesis, as reliable noninvasive assessment of energy balance. These studies were in fact essentially based on correlations between urinary C‐peptide and energy intake and only focused on nonhuman primates. During a bed‐rest study conducted in 16 healthy women in a controlled environment, we tested the existence of a relationship between 24 h‐urinary C‐peptide and energy balance in humans. Daily energy intake and body mass, body composition (dual‐energy X‐ray absorptiometry (DXA)) and total energy expenditure (doubly labeled water (DLW) method) was measured and energy balance was calculated as the difference between energy intake and expenditure. Urinary C‐peptide was positively correlated with bed‐rest‐induced changes in fat mass (r² = 0.285; P = 0.03) and energy balance assessed at the end of the bed‐rest (r² = 0.302; P = 0.027). However, in this tightly controlled environment, urinary C‐peptide only accounted for 30% of variations in energy balance. No relationship was noted between urinary C‐peptide and body or fat mass both at baseline and at the end of the bed‐rest. These results indicate that urinary C‐peptide cannot be used as an accurate biomarker of energy balance in the general human population in free‐living conditions.     

A Method to Constrain Genome-Scale Models with 13C Labeling Data

Current limitations in quantitatively predicting biological behavior hinder our efforts to engineer biological systems to produce biofuels and other desired chemicals. Here, we present a new method for calculating metabolic fluxes, key targets in metabolic engineering, that incorporates data from ¹³C labeling experiments and genome-scale models. The data from ¹³C labeling experiments provide strong flux constraints that eliminate the need to assume an evolutionary optimization principle such as the growth rate optimization assumption used in Flux Balance Analysis (FBA). This effective constraining is achieved by making the simple but biologically relevant assumption that flux flows from core to peripheral metabolism and does not flow back. The new method is significantly more robust than FBA with respect to errors in genome-scale model reconstruction. Furthermore, it can provide a comprehensive picture of metabolite balancing and predictions for unmeasured extracellular fluxes as constrained by ¹³C labeling data. A comparison shows that the results of this new method are similar to those found through ¹³C Metabolic Flux Analysis (¹³C MFA) for central carbon metabolism but, additionally, it provides flux estimates for peripheral metabolism. The extra validation gained by matching 48 relative labeling measurements is used to identify where and why several existing COnstraint Based Reconstruction and Analysis (COBRA) flux prediction algorithms fail. We demonstrate how to use this knowledge to refine these methods and improve their predictive capabilities. This method provides a reliable base upon which to improve the design of biological systems.     

Modeling <Emphasis Type="Italic">Lactococcus lactis</Emphasis> using a genome-scale flux model

Background  

Genome-scale flux models are useful tools to represent and analyze microbial metabolism. In this work we reconstructed the metabolic network of the lactic acid bacteria Lactococcus lactis and developed a genome-scale flux model able to simulate and analyze network capabilities and whole-cell function under aerobic and anaerobic continuous cultures. Flux balance analysis (FBA) and minimization of metabolic adjustment (MOMA) were used as modeling frameworks.     

MapMaker and PathTracer for tracking carbon in genome‐scale metabolic models

Constraint‐based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome‐scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint‐based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO₂ fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint‐based models to identify feasible, active, and high flux paths between metabolites.     

A hybrid kinetic and constraint-based model of leaf metabolism allows predictions of metabolic fluxes in different environments

While flux balance analysis (FBA) provides a framework for predicting steady-state leaf metabolic network fluxes, it does not readily capture the response to environmental variables without being coupled to other modelling formulations. To address this, we coupled an FBA model of 903 reactions of soybean (Glycine max) leaf metabolism with e-photosynthesis, a dynamic model that captures the kinetics of 126 reactions of photosynthesis and associated chloroplast carbon metabolism. Successful coupling was achieved in an iterative formulation in which fluxes from e-photosynthesis were used to constrain the FBA model and then, in turn, fluxes computed from the FBA model used to update parameters in e-photosynthesis. This process was repeated until common fluxes in the two models converged. Coupling did not hamper the ability of the kinetic module to accurately predict the carbon assimilation rate, photosystem II electron flux, and starch accumulation of field-grown soybean at two CO₂ concentrations. The coupled model also allowed accurate predictions of additional parameters such as nocturnal respiration, as well as analysis of the effect of light intensity and elevated CO₂ on leaf metabolism. Predictions included an unexpected decrease in the rate of export of sucrose from the leaf at high light, due to altered starch–sucrose partitioning, and altered daytime flux modes in the tricarboxylic acid cycle at elevated CO₂. Mitochondrial fluxes were notably different between growing and mature leaves, with greater anaplerotic, tricarboxylic acid cycle and mitochondrial ATP synthase fluxes predicted in the former, primarily to provide carbon skeletons and energy for protein synthesis.     

Plankton metabolic balance at the margins of very large lakes: temporal variability and evidence for dominance of autochthonous processes

1. Planktonic metabolic balance (PMB_m) of the surface mixed layer (SML) was measured as the ratio of areal rates of gross photosynthesis (AGP) to community respiration (AR) to test the idea that previously neglected allochthonous inputs of organic matter may support chronic excess respiration relative to photosynthesis even in very large lakes during the summer (May–October) season. Four Laurentian Great Lakes coastal sites of varying trophic status, physical structure and dissolved organic carbon (DOC) concentration were studied with oxygen light‐and dark bottle and ¹⁴C methods, with excess respiration anticipated in the higher DOC sites. 2. Planktonic metabolic balance was net autotrophic in 73% of the observations. The calculated mixing depth at which respiration would predominate over photosynthesis was greater than typically observed mixing depths, varying from 11 to 25 m in the more transparent, low DOC (<3 g m⁻³) sites to 8–15 m in the higher DOC (4–6 g m⁻³) sites. Biweekly measurements at one higher and one lower DOC site over two successive summer seasons showed that seasonal gross photosynthesis (ΣAGP) exceeded seasonal community respiration (ΣAR). Despite the location of the sites at the periphery of the lakes, where allochthonous influences should be strongest, the measurements indicated prevailing conditions of net autotrophy in the SML. 3. Individual measurements of AR from this study and the literature were correlated with AGP but season average values were more tightly correlated, suggesting a tighter coupling of metabolic rates on a larger scale and a looser coupling on a shorter scale. The observed temporal variability was variable in pattern among years, and likely to confound inferences based on limited sampling. 4. It is shown that accepted formulations for AGP and AR lead to the conclusion that PMB_m should be largely predictable from knowledge of a biological properties ratio (light‐saturated gross photosynthesis to plankton community respiration, P_max/R) and a physical properties ratio (euphotic to mixing depths, Z_eu/Z_m) and this prediction was confirmed using data from this study and from the literature. The evident success of this model points to the pre‐eminent importance of plankton biomass and physical conditions in determining metabolic balance. Variation in these fundamental factors appears capable of explaining the diversity of PMB_m reported for different Great Lakes.     

Flux balance analysis of Chlorella sp. FC2 IITG under photoautotrophic and heterotrophic growth conditions

Quantification of carbon flux distribution in the metabolic network of microalgae remains important to understand the complex interplay between energy metabolism, carbon fixation, and assimilation pathways. This is even more relevant with respect to cyclic metabolism of microalgae under light–dark cycle. In the present study, flux balance analysis (FBA) was carried out for an indigenous isolate Chlorella sp. FC2 IITG under photoautotrophic and heterotrophic growth conditions. A shift in intracellular flux distribution was predicted during transition from nutrient sufficient phase to nutrient starvation phase of growth. Further, dynamic flux analysis (dFBA) was carried out to capture light–dark metabolism over discretized pseudo steady state time intervals. Our key findings include the following: (i) unlike heterotrophic condition, oxidative pentose phosphate (PP) pathway, and Krebs cycle were relatively inactive under photoautotrophic growth; (ii) in both growth conditions, while transhydrogenation reaction was highly active, glyoxalate shunt was found to be nonoperative; (iii) flux distribution during transition period was marked with up regulation of carbon flux toward nongrowth associated (NGA) maintenance energy, oxidative phosphorylation, and photophosphorylation; (iv) redirection of carbon flux from polysaccharide and neutral lipid resulted in up regulation of Krebs cycle flux in the dark phase; (v) elevated glycolytic and acetyl-CoA flux were coupled with induction of neutral lipid during light cycle of the growth; (vi) significantly active photophosphorylation in the light phase was able to satisfy cellular energy requirement without need of oxidative PP pathway; and (vi) unlike static FBA, dFBA predicted an unaltered NGA maintenance energy of 1.5 mmol g^?1 DCW h^?1.