Constraint-based Flux Analysis Research Articles

Genome sequencing and metabolic network reconstruction of a novel sulfur-oxidizing bacterium Acidithiobacillus Ameehan.

Sulfur-oxidizing bacteria play a crucial role in various processes, including mine bioleaching, biodesulfurization, and treatment of sulfur-containing wastewater. Nevertheless, the pathway involved in sulfur oxidation is highly intricate, making it complete comprehension a formidable and protracted undertaking. The mechanisms of sulfur oxidation within the Acidithiobacillus genus, along with the process of energy production, remain areas that necessitate further research and elucidation. In this study, a novel strain of sulfur-oxidizing bacterium, Acidithiobacillus Ameehan, was isolated. Several physiological characteristics of the strain Ameehan were verified and its complete genome sequence was presented in the study. Besides, the first genome-scale metabolic network model (AMEE_WP1377) was reconstructed for Acidithiobacillus Ameehan to gain a comprehensive understanding of the metabolic capacity of the strain.The characteristics of Acidithiobacillus Ameehan included morphological size and an optimal growth temperature range of 37-45°C, as well as an optimal growth pH range of pH 2.0-8.0. The microbe was found to be capable of growth when sulfur and K2O6S4 were supplied as the energy source and electron donor for CO2 fixation. Conversely, it could not utilize Na2S2O3, FeS2, and FeSO4·7H2O as the energy source or electron donor for CO2 fixation, nor could it grow using glucose or yeast extract as a carbon source. Genome annotation revealed that the strain Ameehan possessed a series of sulfur oxidizing genes that enabled it to oxidize elemental sulfur or various reduced inorganic sulfur compounds (RISCs). In addition, the bacterium also possessed carbon fixing genes involved in the incomplete Calvin-Benson-Bassham (CBB) cycle. However, the bacterium lacked the ability to oxidize iron and fix nitrogen. By implementing a constraint-based flux analysis to predict cellular growth in the presence of 71 carbon sources, 88.7% agreement with experimental Biolog data was observed. Five sulfur oxidation pathways were discovered through model simulations. The optimal sulfur oxidation pathway had the highest ATP production rate of 14.81 mmol/gDW/h, NADH/NADPH production rate of 5.76 mmol/gDW/h, consumed 1.575 mmol/gDW/h of CO2, and 1.5 mmol/gDW/h of sulfur. Our findings provide a comprehensive outlook on the most effective cellular metabolic pathways implicated in sulfur oxidation within Acidithiobacillus Ameehan. It suggests that the OMP (outer-membrane proteins) and SQR enzymes (sulfide: quinone oxidoreductase) have a significant impact on the energy production efficiency of sulfur oxidation, which could have potential biotechnological applications.

Genome-scale metabolic modeling and in silico analysis of opportunistic skin pathogen Cutibacterium acnes.

Cutibacterium acnes, one of the most abundant skin microbes found in the sebaceous gland, is known to contribute to the development of acne vulgaris when its strains become imbalanced. The current limitations of acne treatment using antibiotics have caused an urgent need to develop a systematic strategy for selectively targeting C. acnes, which can be achieved by characterizing their cellular behaviors under various skin environments. To this end, we developed a genome-scale metabolic model (GEM) of virulent C. acnes, iCA843, based on the genome information of a relevant strain from ribotype 5 to comprehensively understand the pathogenic traits of C. acnes in the skin environment. We validated the model qualitatively by demonstrating its accuracy prediction of propionate and acetate production patterns, which were consistent with experimental observations. Additionally, we identified unique biosynthetic pathways for short-chain fatty acids in C. acnes compared to other GEMs of acne-inducing skin pathogens. By conducting constraint-based flux analysis under endogenous carbon sources in human skin, we discovered that the Wood-Werkman cycle is highly activated under acnes-associated skin condition for the regeneration of NAD, resulting in enhanced propionate production. Finally, we proposed potential anti-C. acnes targets by using the model-guided systematic framework based on gene essentiality analysis and protein sequence similarity search with abundant skin microbiome taxa.

High-quality genome-scale metabolic network reconstruction of probiotic bacterium Escherichia coli Nissle 1917

BackgroundEscherichia coli Nissle 1917 (EcN) is a probiotic bacterium used to treat various gastrointestinal diseases. EcN is increasingly being used as a chassis for the engineering of advanced microbiome therapeutics. To aid in future engineering efforts, our aim was to construct an updated metabolic model of EcN with extended secondary metabolite representation.ResultsAn updated high-quality genome-scale metabolic model of EcN, iHM1533, was developed based on comparison with 55 E. coli/Shigella reference GEMs and manual curation, including expanded secondary metabolite pathways (enterobactin, salmochelins, aerobactin, yersiniabactin, and colibactin). The model was validated and improved using phenotype microarray data, resulting in an 82.3% accuracy in predicting growth phenotypes on various nutrition sources. Flux variability analysis with previously published 13C fluxomics data validated prediction of the internal central carbon fluxes. A standardised test suite called Memote assessed the quality of iHM1533 to have an overall score of 89%. The model was applied by using constraint-based flux analysis to predict targets for optimisation of secondary metabolite production. Modelling predicted design targets from across amino acid metabolism, carbon metabolism, and other subsystems that are common or unique for influencing the production of various secondary metabolites.ConclusioniHM1533 represents a well-annotated metabolic model of EcN with extended secondary metabolite representation. Phenotype characterisation and the iHM1533 model provide a better understanding of the metabolic capabilities of EcN and will help future metabolic engineering efforts.

Comprehensive analysis of a metabolic model for lipid production in Rhodosporidium toruloides

The yeast Rhodosporidium toruloides has been extensively studied for its application in biolipid production. The knowledge of its metabolism capabilities and the application of constraint-based flux analysis methodology provide useful information for process prediction and optimization. The accuracy of the resulting predictions is highly dependent on metabolic models. A metabolic reconstruction for R. toruloides metabolism has been recently published. On the basis of this model, we developed a curated version that unblocks the central nitrogen metabolism and, in addition, completes charge and mass balances in some reactions neglected in the former model. Then, a comprehensive analysis of network capability was performed with the curated model and compared with the published metabolic reconstruction. The flux distribution obtained by lipid optimization with flux balance analysis was able to replicate the internal biochemical changes that lead to lipogenesis in oleaginous microorganisms. These results motivate the development of a genome-scale model for complete elucidation of R. toruloides metabolism.

Genome-scale modeling and transcriptome analysis of Leuconostoc mesenteroides unravel the redox governed metabolic states in obligate heterofermentative lactic acid bacteria

Obligate heterofermentative lactic acid bacteria (LAB) are well-known for their beneficial health effects in humans. To delineate the incompletely characterized metabolism that currently limits their exploitation, at systems-level, we developed a genome-scale metabolic model of the representative obligate heterofermenting LAB, Leuconostoc mesenteroides (iLME620). Constraint-based flux analysis was then used to simulate several qualitative and quantitative phenotypes of L. mesenteroides, thereby evaluating the model validity. With established predictive capabilities, we subsequently employed iLME620 to elucidate unique metabolic characteristics of L. mesenteroides, such as the limited ability to utilize amino acids as energy source, and to substantiate the role of malolactic fermentation (MLF) in the reduction of pH-homeostatic burden on F0F1-ATPase. We also reported new hypothesis on the MLF mechanism that could be explained via a substrate channelling-like phenomenon mainly influenced by intracellular redox state rather than the intermediary reactions. Model simulations further revealed possible proton-symporter dependent activity of the energy efficient glucose-phosphotransferase system in obligate heterofermentative LAB. Moreover, integrated transcriptomic analysis allowed us to hypothesize transcriptional regulatory bias affecting the intracellular redox state. The insights gained here about the low ATP-yielding metabolism of L. mesenteroides, dominantly controlled by the cellular redox state, could potentially aid strain design for probiotic and cell factory applications.

Construction and simulation of the Bradyrhizobium diazoefficiens USDA110 metabolic network: a comparison between free-living and symbiotic states.

Bradyrhizobium diazoefficiens is a rhizobium able to convert atmospheric nitrogen into ammonium by establishing mutualistic symbiosis with soybean. It has been recognized as an important parent strain for microbial agents and is widely applied in agricultural and environmental fields. In order to study the metabolic properties of symbiotic nitrogen fixation and the differences between a free-living cell and a symbiotic bacteroid, a genome-scale metabolic network of B. diazoefficiens USDA110 was constructed and analyzed. The metabolic network, iYY1101, contains 1031 reactions, 661 metabolites, and 1101 genes in total. Metabolic models reflecting free-living and symbiotic states were determined by defining the corresponding objective functions and substrate input sets, and were further constrained by high-throughput transcriptomic and proteomic data. Constraint-based flux analysis was used to compare the metabolic capacities and the effects on the metabolic targets of genes and reactions between the two physiological states. The results showed that a free-living rhizobium possesses a steady state flux distribution for sustaining a complex supply of biomass precursors while a symbiotic bacteroid maintains a relatively condensed one adapted to nitrogen-fixation. Our metabolic models may serve as a promising platform for better understanding the symbiotic nitrogen fixation of this species.

FOCuS: a metaheuristic algorithm for computing knockouts from genome-scale models for strain optimization.

Although handful tools are available for constraint-based flux analysis to generate knockout strains, most of these are either based on bilevel-MIP or its modifications. However, metaheuristic approaches that are known for their flexibility and scalability have been less studied. Moreover, in the existing tools, sectioning of search space to find optimal knocks has not been considered. Herein, a novel computational procedure, termed as FOCuS (Flower-pOllination coupled Clonal Selection algorithm), was developed to find the optimal reaction knockouts from a metabolic network to maximize the production of specific metabolites. FOCuS derives its benefits from nature-inspired flower pollination algorithm and artificial immune system-inspired clonal selection algorithm to converge to an optimal solution. To evaluate the performance of FOCuS, reported results obtained from both MIP and other metaheuristic-based tools were compared in selected case studies. The results demonstrated the robustness of FOCuS irrespective of the size of metabolic network and number of knockouts. Moreover, sectioning of search space coupled with pooling of priority reactions based on their contribution to objective function for generating smaller search space significantly reduced the computational time.

Genome-scale metabolic modeling and in silico analysis of lipid accumulating yeast Candida tropicalis for dicarboxylic acid production.

Recently, the bio-production of α,ω-dicarboxylic acids (DCAs) has gained significant attention, which potentially leads to the replacement of the conventional petroleum-based products. In this regard, the lipid accumulating yeast Candida tropicalis, has been recognized as a promising microbial host for DCA biosynthesis: it possess the unique ω-oxidation pathway where the terminal carbon of α-fatty acids is oxidized to form DCAs with varying chain lengths. However, despite such industrial importance, its cellular physiology and lipid accumulation capability remain largely uncharacterized. Thus, it is imperative to better understand the metabolic behavior of this lipogenic yeast, which could be achieved by a systems biological approach. To this end, herein, we reconstructed the genome-scale metabolic model of C. tropicalis, iCT646, accounting for 646 unique genes, 945 metabolic reactions, and 712 metabolites. Initially, the comparative network analysis of iCT646 with other yeasts revealed several distinctive metabolic reactions, mainly within the amino acid and lipid metabolism including the ω-oxidation pathway. Constraints-based flux analysis was, then, employed to predict the in silico growth rates of C. tropicalis which are highly consistent with the cellular phenotype observed in glucose and xylose minimal media chemostat cultures. Subsequently, the lipid accumulation capability of C. tropicalis was explored in comparison with Saccharomyces cerevisiae, indicating that the formation of "citrate pyruvate cycle" is essential to the lipid accumulation in oleaginous yeasts. The in silico flux analysis also highlighted the enhanced ability of pentose phosphate pathway as NADPH source rather than malic enzyme during lipogenesis. Finally, iCT646 was successfully utilized to highlight the key directions of C. tropicalis strain design for the whole cell biotransformation application to produce long-chain DCAs from alkanes. Biotechnol. Bioeng. 2016;113: 1993-2004. © 2016 Wiley Periodicals, Inc.

CBFA: phenotype prediction integrating metabolic models with constraints derived from experimental data.

BackgroundFlux analysis methods lie at the core of Metabolic Engineering (ME), providing methods for phenotype simulation that allow the determination of flux distributions under different conditions. Although many constraint-based modeling software tools have been developed and published, none provides a free user-friendly application that makes available the full portfolio of flux analysis methods.ResultsThis work presents Constraint-based Flux Analysis (CBFA), an open-source software application for flux analysis in metabolic models that implements several methods for phenotype prediction, allowing users to define constraints associated with measured fluxes and/or flux ratios, together with environmental conditions (e.g. media) and reaction/gene knockouts. CBFA identifies the set of applicable methods based on the constraints defined from user inputs, encompassing algebraic and constraint-based simulation methods. The integration of CBFA within the OptFlux framework for ME enables the utilization of different model formats and standards and the integration with complementary methods for phenotype simulation and visualization of results.ConclusionsA general-purpose and flexible application is proposed that is independent of the origin of the constraints defined for a given simulation. The aim is to provide a simple to use software tool focused on the application of several flux prediction methods.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-014-0123-1) contains supplementary material, which is available to authorized users.

Abstract 3372: Constraints-based metabolic flux analysis approach links tumor stage to metabolic adaptations and survival in cancer cells

Abstract The Warburg effect has been observed in many cancers and their high glycolytic capacity has signified their dependence on glucose. More recently, glutamine has emerged not only as an important nutrient for many cancers, but also necessary for their elevated energetic requirements. Due to these high energetic demands, certain cancer cells become addicted to glutamine to maintain viability. We postulate that distinct metabolic reconfigurations of certain cancers define their dependence or independence on glutamine for survival while maintaining their proliferative propensity and redox status. An intricate picture of the metabolic profiles is to be drawn from estimating intracellular fluxes, by combining stable isotope tracer measurements and experimental metabolomics data from different cancer cell lines, which have been observed to be glutamine dependent and independent. To this extent, we describe an approach that utilizes a redox-balanced model incorporating the electron transport chain and comprehensive amino-acid metabolic reactions to elucidate the importance of oxidative phosphorylation, often overlooked in classical approaches. Diving deeper into the foray of metabolic reprogramming, we perform in silico experiments using a constraint-based multi-objective modeling approach. This methodology elucidates the switching of metabolic pathways in glutamine-dependent and -independent cancers under nutrient-available and nutrient-deprived conditions. Our approach assumes that cancer cells operate at optimal levels, maintaining multiple objectives under certain environmental conditions. Constraining the proliferative phenotype from a maximal to minimal levels of the cells under different nutrient conditions emulates the observed behavior of glutamine-dependent cells under deprivation conditions and contrasts their metabolic reprogramming against that of glutamine-independent cells. We corroborated our simulations with experimentally derived metabolic fluxes and found that glutamine anabolism over catabolism dictates adaptations and survival in invasive cancers. Our results will lead to identification of potential targets for inducing nutrient-sensitivity and enhance current therapeutic approaches. Citation Format: Abhinav Achreja, Lifeng Yang, Hongyun Zhao, Juan Marini, Deepak Nagrath. Constraints-based metabolic flux analysis approach links tumor stage to metabolic adaptations and survival in cancer cells. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3372. doi:10.1158/1538-7445.AM2014-3372

Genome-scale metabolic network reconstruction and in silico flux analysis of the thermophilic bacterium Thermus thermophilus HB27.

BackgroundThermus thermophilus, an extremely thermophilic bacterium, has been widely recognized as a model organism for studying how microbes can survive and adapt under high temperature environment. However, the thermotolerant mechanisms and cellular metabolism still remains mostly unravelled. Thus, it is highly required to consider systems biological approaches where T. thermophilus metabolic network model can be employed together with high throughput experimental data for elucidating its physiological characteristics under such harsh conditions.ResultsWe reconstructed a genome-scale metabolic model of T. thermophilus, iTT548, the first ever large-scale network of a thermophilic bacterium, accounting for 548 unique genes, 796 reactions and 635 unique metabolites. Our initial comparative analysis of the model with Escherichia coli has revealed several distinctive metabolic reactions, mainly in amino acid metabolism and carotenoid biosynthesis, producing relevant compounds to retain the cellular membrane for withstanding high temperature. Constraints-based flux analysis was, then, applied to simulate the metabolic state in glucose minimal and amino acid rich media. Remarkably, resulting growth predictions were highly consistent with the experimental observations. The subsequent comparative flux analysis under different environmental conditions highlighted that the cells consumed branched chain amino acids preferably and utilized them directly in the relevant anabolic pathways for the fatty acid synthesis. Finally, gene essentiality study was also conducted via single gene deletion analysis, to identify the conditional essential genes in glucose minimal and complex media.ConclusionsThe reconstructed genome-scale metabolic model elucidates the phenotypes of T. thermophilus, thus allowing us to gain valuable insights into its cellular metabolism through in silico simulations. The information obtained from such analysis would not only shed light on the understanding of physiology of thermophiles but also helps us to devise metabolic engineering strategies to develop T. thermophilus as a thermostable microbial cell factory.

COFACTOR MODIFICATION ANALYSIS: A COMPUTATIONAL FRAMEWORK TO IDENTIFY COFACTOR SPECIFICITY ENGINEERING TARGETS FOR STRAIN IMPROVEMENT

Cofactors, such as NAD(H) and NADP(H), play important roles in energy transfer within the cells by providing the necessary redox carriers for a myriad of metabolic reactions, both anabolic and catabolic. Thus, it is crucial to establish the overall cellular redox balance for achieving the desired cellular physiology. Of several methods to manipulate the intracellular cofactor regeneration rates, altering the cofactor specificity of a particular enzyme is a promising one. However, the identification of relevant enzyme targets for such cofactor specificity engineering (CSE) is often very difficult and labor intensive. Therefore, it is necessary to develop more systematic approaches to find the cofactor engineering targets for strain improvement. Presented herein is a novel mathematical framework, cofactor modification analysis (CMA), developed based on the well-established constraints-based flux analysis, for the systematic identification of suitable CSE targets while exploring the global metabolic effects. The CMA algorithm was applied to E. coli using its genome-scale metabolic model, iJO1366, thereby identifying the growth-coupled cofactor engineering targets for overproducing four of its native products: acetate, formate, ethanol, and lactate, and three non-native products: 1-butanol, 1,4-butanediol, and 1,3-propanediol. Notably, among several target candidates for cofactor engineering, glyceraldehyde-3-phosphate dehydrogenase (GAPD) is the most promising enzyme; its cofactor modification enhanced both the desired product and biomass yields significantly. Finally, given the identified target, we further discussed potential mutational strategies for modifying cofactor specificity of GAPD in E. coli as suggested by in silico protein docking experiments.

FastPros: screening of reaction knockout strategies for metabolic engineering

Motivation: Although constraint-based flux analysis of knockout strains has facilitated the production of desirable metabolites in microbes, current screening methods have placed a limitation on the number knockouts that can be simultaneously analyzed.Results: Here, we propose a novel screening method named FastPros. In this method, the potential of a given reaction knockout for production of a specific metabolite is evaluated by shadow pricing of the constraint in the flux balance analysis, which generates a screening score to obtain candidate knockout sets. To evaluate the performance of FastPros, we screened knockout sets to produce each metabolite in the entire Escherichia coli metabolic network. We found that 75% of these metabolites could be produced under biomass maximization conditions by adding up to 25 reaction knockouts. Furthermore, we demonstrated that using FastPros in tandem with another screening method, OptKnock, could further improve target metabolite productivity.Availability and implementation: Source code is freely available at http://www-shimizu.ist.osaka-u.ac.jp/shimizu_lab/FastPros/, implemented in MATLAB and COBRA toolbox.Contact: chikara.furusawa@riken.jp or shimizu@ist.osaka-u.ac.jpSupplementary information: Supplementary data are available at Bioinformatics online.

Flux‐coupled genes and their use in metabolic flux analysis

As large volumes of omics data have become available, systems biology is playing increasingly important roles in elucidating new biological phenomena, especially through genome-scale metabolic network modeling and simulation. Much effort has been exerted on integrating omics data with metabolic flux simulation, but further development is necessary for more accurate flux estimation. To move one step forward, we adopted the concept of flux-coupled genes (FCGs), which show that their expression transition patterns upon perturbations are correlated with their corresponding flux values, as additional constraints in metabolic flux analysis. It was found that gnd, pfkB, rpe, sdhB, sdhD, sucA, and zwf genes, mostly associated with pentose phosphate pathway and TCA cycle, were the most consistent FCGs in Escherichia coli based on its transcriptome and (13) C-flux data obtained from the chemostat cultivation at five different dilution rates. Consequently, constraints-based flux analyses with FCGs as additional constraints were conducted for the seven single-gene knockout mutants, compared with those obtained without using FCGs. This strategy of constraining the metabolic flux analysis with FCGs is expected to be useful due to the relative ease in obtaining transcriptional information in the functional genomics era.

Genome-scale reconstruction and in silico analysis of Klebsiella oxytoca for 2,3-butanediol production

BackgroundKlebsiella oxytoca, a Gram-negative, rod-shaped, and facultative anaerobic bacterium, is one of the most promising 2,3-butanediol (2,3-BD) producers. In order to improve the metabolic performance of K. oxytoca as an efficient biofactory, it is necessary to assess its metabolic characteristics with a system-wide scope, and to optimize the metabolic pathways at a systems level. Provision of the complete genome sequence of K. oxytoca enabled the construction of genome-scale metabolic model of K. oxytoca and its in silico analyses.ResultsThe genome-scale metabolic model of K. oxytoca was constructed using the annotated genome with biochemical and physiological information. The stoichiometric model, KoxGSC1457, is composed of 1,457 reactions and 1,099 metabolites. The model was further refined by applying biomass composition equations and comparing in silico results with experimental data based on constraints-based flux analyses. Then, the model was applied to in silico analyses to understand the properties of K. oxytoca and also to improve its capabilities for 2,3-BD production according to genetic and environmental perturbations. Firstly, in silico analysis, which tested the effect of augmenting the metabolic flux pool of 2,3-BD precursors, elucidated that increasing the pyruvate pool is primarily important for 2,3-BD synthesis. Secondly, we performed in silico single gene knockout simulation for 2,3-BD overproduction, and investigated the changes of the in silico flux solution space of a ldhA gene knockout mutant in comparison with that of the wild-type strain. Finally, the KoxGSC1457 model was used to optimize the oxygen levels during fermentation for 2,3-BD production.ConclusionsThe genome-scale metabolic model, KoxGSC1457, constructed in this study successfully investigated metabolic characteristics of K. oxytoca at systems level. The KoxGSC1457 model could be employed as an useful tool to analyze its metabolic capabilities, to predict its physiological responses according to environmental and genetic perturbations, and to design metabolic engineering strategies to improve its metabolic performance.

Flux variability scanning based on enforced objective flux for identifying gene amplification targets

BackgroundIn order to reduce time and efforts to develop microbial strains with better capability of producing desired bioproducts, genome-scale metabolic simulations have proven useful in identifying gene knockout and amplification targets. Constraints-based flux analysis has successfully been employed for such simulation, but is limited in its ability to properly describe the complex nature of biological systems. Gene knockout simulations are relatively straightforward to implement, simply by constraining the flux values of the target reaction to zero, but the identification of reliable gene amplification targets is rather difficult. Here, we report a new algorithm which incorporates physiological data into a model to improve the model’s prediction capabilities and to capitalize on the relationships between genes and metabolic fluxes.ResultsWe developed an algorithm, flux variability scanning based on enforced objective flux (FVSEOF) with grouping reaction (GR) constraints, in an effort to identify gene amplification targets by considering reactions that co-carry flux values based on physiological omics data via “GR constraints”. This method scans changes in the variabilities of metabolic fluxes in response to an artificially enforced objective flux of product formation. The gene amplification targets predicted using this method were validated by comparing the predicted effects with the previous experimental results obtained for the production of shikimic acid and putrescine in Escherichia coli. Moreover, new gene amplification targets for further enhancing putrescine production were validated through experiments involving the overexpression of each identified targeted gene under condition-controlled batch cultivation.ConclusionsFVSEOF with GR constraints allows identification of gene amplification targets for metabolic engineering of microbial strains in order to enhance the production of desired bioproducts. The algorithm was validated through the experiments on the enhanced production of putrescine in E. coli, in addition to the comparison with the previously reported experimental data. The FVSEOF strategy with GR constraints will be generally useful for developing industrially important microbial strains having enhanced capabilities of producing chemicals of interest.

Genome-scale reconstruction and in silico analysis of the Ralstonia eutropha H16 for polyhydroxyalkanoate synthesis, lithoautotrophic growth, and 2-methyl citric acid production

BackgroundRalstonia eutropha H16, found in both soil and water, is a Gram-negative lithoautotrophic bacterium that can utillize CO2 and H2 as its sources of carbon and energy in the absence of organic substrates. R. eutropha H16 can reach high cell densities either under lithoautotrophic or heterotrophic conditions, which makes it suitable for a number of biotechnological applications. It is the best known and most promising producer of polyhydroxyalkanoates (PHAs) from various carbon substrates and is an environmentally important bacterium that can degrade aromatic compounds. In order to make R. eutropha H16 a more efficient and robust biofactory, system-wide metabolic engineering to improve its metabolic performance is essential. Thus, it is necessary to analyze its metabolic characteristics systematically and optimize the entire metabolic network at systems level.ResultsWe present the lithoautotrophic genome-scale metabolic model of R. eutropha H16 based on the annotated genome with biochemical and physiological information. The stoichiometic model, RehMBEL1391, is composed of 1391 reactions including 229 transport reactions and 1171 metabolites. Constraints-based flux analyses were performed to refine and validate the genome-scale metabolic model under environmental and genetic perturbations. First, the lithoautotrophic growth characteristics of R. eutropha H16 were investigated under varying feeding ratios of gas mixture. Second, the genome-scale metabolic model was used to design the strategies for the production of poly[R-(-)-3hydroxybutyrate] (PHB) under different pH values and carbon/nitrogen source uptake ratios. It was also used to analyze the metabolic characteristics of R. eutropha when the phosphofructokinase gene was expressed. Finally, in silico gene knockout simulations were performed to identify targets for metabolic engineering essential for the production of 2-methylcitric acid in R. eutropha H16.ConclusionThe genome-scale metabolic model, RehMBEL1391, successfully represented metabolic characteristics of R. eutropha H16 at systems level. The reconstructed genome-scale metabolic model can be employed as an useful tool for understanding its metabolic capabilities, predicting its physiological consequences in response to various environmental and genetic changes, and developing strategies for systems metabolic engineering to improve its metabolic performance.

Framework for network modularization and Bayesian network analysis to investigate the perturbed metabolic network

BackgroundGenome-scale metabolic network models have contributed to elucidating biological phenomena, and predicting gene targets to engineer for biotechnological applications. With their increasing importance, their precise network characterization has also been crucial for better understanding of the cellular physiology.ResultsWe herein introduce a framework for network modularization and Bayesian network analysis (FMB) to investigate organism’s metabolism under perturbation. FMB reveals direction of influences among metabolic modules, in which reactions with similar or positively correlated flux variation patterns are clustered, in response to specific perturbation using metabolic flux data. With metabolic flux data calculated by constraints-based flux analysis under both control and perturbation conditions, FMB, in essence, reveals the effects of specific perturbations on the biological system through network modularization and Bayesian network analysis at metabolic modular level. As a demonstration, this framework was applied to the genetically perturbed Escherichia coli metabolism, which is a lpdA gene knockout mutant, using its genome-scale metabolic network model.ConclusionsAfter all, it provides alternative scenarios of metabolic flux distributions in response to the perturbation, which are complementary to the data obtained from conventionally available genome-wide high-throughput techniques or metabolic flux analysis.

Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement

BackgroundPichia pastoris has been recognized as an effective host for recombinant protein production. A number of studies have been reported for improving this expression system. However, its physiology and cellular metabolism still remained largely uncharacterized. Thus, it is highly desirable to establish a systems biotechnological framework, in which a comprehensive in silico model of P. pastoris can be employed together with high throughput experimental data analysis, for better understanding of the methylotrophic yeast's metabolism.ResultsA fully compartmentalized metabolic model of P. pastoris (iPP668), composed of 1,361 reactions and 1,177 metabolites, was reconstructed based on its genome annotation and biochemical information. The constraints-based flux analysis was then used to predict achievable growth rate which is consistent with the cellular phenotype of P. pastoris observed during chemostat experiments. Subsequent in silico analysis further explored the effect of various carbon sources on cell growth, revealing sorbitol as a promising candidate for culturing recombinant P. pastoris strains producing heterologous proteins. Interestingly, methanol consumption yields a high regeneration rate of reducing equivalents which is substantial for the synthesis of valuable pharmaceutical precursors. Hence, as a case study, we examined the applicability of P. pastoris system to whole-cell biotransformation and also identified relevant metabolic engineering targets that have been experimentally verified.ConclusionThe genome-scale metabolic model characterizes the cellular physiology of P. pastoris, thus allowing us to gain valuable insights into the metabolism of methylotrophic yeast and devise possible strategies for strain improvement through in silico simulations. This computational approach, combined with synthetic biology techniques, potentially forms a basis for rational analysis and design of P. pastoris metabolic network to enhance humanized glycoprotein production.

The genome-scale metabolic network analysis of Zymomonas mobilis ZM4 explains physiological features and suggests ethanol and succinic acid production strategies

BackgroundZymomonas mobilis ZM4 is a Gram-negative bacterium that can efficiently produce ethanol from various carbon substrates, including glucose, fructose, and sucrose, via the Entner-Doudoroff pathway. However, systems metabolic engineering is required to further enhance its metabolic performance for industrial application. As an important step towards this goal, the genome-scale metabolic model of Z. mobilis is required to systematically analyze in silico the metabolic characteristics of this bacterium under a wide range of genotypic and environmental conditions.ResultsThe genome-scale metabolic model of Z. mobilis ZM4, ZmoMBEL601, was reconstructed based on its annotated genes, literature, physiological and biochemical databases. The metabolic model comprises 579 metabolites and 601 metabolic reactions (571 biochemical conversion and 30 transport reactions), built upon extensive search of existing knowledge. Physiological features of Z. mobilis were then examined using constraints-based flux analysis in detail as follows. First, the physiological changes of Z. mobilis as it shifts from anaerobic to aerobic environments (i.e. aerobic shift) were investigated. Then the intensities of flux-sum, which is the cluster of either all ingoing or outgoing fluxes through a metabolite, and the maximum in silico yields of ethanol for Z. mobilis and Escherichia coli were compared and analyzed. Furthermore, the substrate utilization range of Z. mobilis was expanded to include pentose sugar metabolism by introducing metabolic pathways to allow Z. mobilis to utilize pentose sugars. Finally, double gene knock-out simulations were performed to design a strategy for efficiently producing succinic acid as another example of application of the genome-scale metabolic model of Z. mobilis.ConclusionThe genome-scale metabolic model reconstructed in this study was able to successfully represent the metabolic characteristics of Z. mobilis under various conditions as validated by experiments and literature information. This reconstructed metabolic model will allow better understanding of Z. mobilis metabolism and consequently designing metabolic engineering strategies for various biotechnological applications.