Predicted Flux Distributions Research Articles

Genome-scale flux balance analysis reveals redox trade-offs in the metabolism of the thermoacidophile Methylacidiphilum fumariolicum under auto-, hetero-and methanotrophic conditions

Members of the genus Methylacidiphilum are thermoacidophile methanotrophs with optimal growth temperatures between 50°C and 60°C, and pH between 1.0 and 3.0. These microorganisms, as well as other extremophile bacteria, offer an attractive platform for environmental and industrial biotechnology because of their robust operating conditions and capacity to grow using low-cost substrates. In this study, we isolated Methylacidiphilum fumariolicum str. Pic from a crater lake located in the state of Chiapas, Mexico. We sequenced the genome and built a genome-scale metabolic model. The manually curated model contains 667 metabolites, 729 reactions, and 473 genes. Predicted flux distributions using flux balance analysis identified changes in redox trade-offs under methanotrophic and autotrophic conditions (H2+CO2). This was also predicted under heterotrophic conditions (acetone, isopropanol, and propane). Model validation was performed by testing the capacity of the strains to grow using four substrates: CH4, acetone, isopropanol, and LP-Gas. The results suggest that the metabolism of M. fumariolicum str. Pic is limited by the regeneration of redox equivalents such as NAD(P)H and reduced cytochromes.

Probabilistic thermodynamic analysis of metabolic networks.

MotivationRandom sampling of metabolic fluxes can provide a comprehensive description of the capabilities of a metabolic network. However, current sampling approaches do not model thermodynamics explicitly, leading to inaccurate predictions of an organism’s potential or actual metabolic operations.ResultsWe present a probabilistic framework combining thermodynamic quantities with steady-state flux constraints to analyze the properties of a metabolic network. It includes methods for probabilistic metabolic optimization and for joint sampling of thermodynamic and flux spaces. Applied to a model of Escherichia coli, we use the methods to reveal known and novel mechanisms of substrate channeling, and to accurately predict reaction directions and metabolite concentrations. Interestingly, predicted flux distributions are multimodal, leading to discrete hypotheses on E.coli’s metabolic capabilities.Availability and implementationPython and MATLAB packages available at https://gitlab.com/csb.ethz/pta.Supplementary informationSupplementary data are available at Bioinformatics online.

Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content.

Marine cyanobacteria are promising microbes to capture and convert atmospheric CO2 and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.

A genome-scale metabolic network reconstruction of extremely halophilic bacterium Salinibacter ruber.

A genome-scale metabolic network reconstruction of Salinibacter ruber DSM13855 is presented here. To our knowledge, this is the first metabolic model of an organism in the phylum Rhodothermaeota. This model, which will be called iMB631, was reconstructed based on genomic and biochemical data available on the strain Salinibacter ruber DSM13855. This network consists of 1459 reactions, 1363 metabolites and 631 genes. Model evaluation was performed based on existing biochemical data in the literature and also by performing laboratory experiments. For growth on different carbon sources, we show that iMB631 is able to correctly predict the growth in 91% of cases where growth has been observed experimentally and 83% of conditions in which S. ruber did not grow. The F-score was 93%, demonstrating a generally acceptable performance of the model. Based on the predicted flux distributions, we found that under certain autotrophic condition, a reductive tricarboxylic acid cycle (rTCA) has fluxes in all necessary reactions to support autotrophic growth. To include special metabolites of the bacterium, salinixanthin biosynthesis pathway was modeled based on the pathway proposed recently. For years, main glucose consumption pathway has been under debates in S. ruber. Using flux balance analysis, iMB631 predicts pentose phosphate pathway, rather than glycolysis, as the active glucose consumption method in the S. ruber.

Predicting Dynamic Metabolic Demands in the Photosynthetic Eukaryote Chlorella vulgaris.

Phototrophic organisms exhibit a highly dynamic proteome, adapting their biomass composition in response to diurnal light/dark cycles and nutrient availability. Here, we used experimentally determined biomass compositions over the course of growth to determine and constrain the biomass objective function (BOF) in a genome-scale metabolic model of Chlorella vulgaris UTEX 395 over time. Changes in the BOF, which encompasses all metabolites necessary to produce biomass, influence the state of the metabolic network thus directly affecting predictions. Simulations using dynamic BOFs predicted distinct proteome demands during heterotrophic or photoautotrophic growth. Model-driven analysis of extracellular nitrogen concentrations and predicted nitrogen uptake rates revealed an intracellular nitrogen pool, which contains 38% of the total nitrogen provided in the medium for photoautotrophic and 13% for heterotrophic growth. Agreement between flux and gene expression trends was determined by statistical comparison. Accordance between predicted flux trends and gene expression trends was found for 65% of multisubunit enzymes and 75% of allosteric reactions. Reactions with the highest agreement between simulations and experimental data were associated with energy metabolism, terpenoid biosynthesis, fatty acids, nucleotides, and amino acid metabolism. Furthermore, predicted flux distributions at each time point were compared with gene expression data to gain new insights into intracellular compartmentalization, specifically for transporters. A total of 103 genes related to internal transport reactions were identified and added to the updated model of C. vulgaris, iCZ946, thus increasing our knowledgebase by 10% for this model green alga.

Decomposition of Flux Distributions into Metabolic Pathways

Genome-scale reconstructions are often used for studying relationships between fundamental components of a metabolic system. In this study, we develop a novel computational method for analyzing predicted flux distributions for metabolic reconstructions. Because chemical reactions may have multiple reactants and products, a directed hypergraph where hyperarcs may have multiple tail vertices and head vertices is a more appropriate representation of the metabolic network than a conventional network. We use this view to represent predicted flux distributions by maximum generalized flows on hypergraphs. We then demonstrate that the generalized hyperflow problem may be transformed to an equivalent network flow problem with side constraints. This transformation allows a flux to be decomposed into chains of reactions. Subsequent analysis of these chains helps to characterize active pathways in a flux distribution. Such characterizations facilitate comparisons of flux distributions for different environmental conditions. The proposed method is applied to compare predicted flux distributions for Salmonella typhimurium to study changes in metabolism that cause enhanced virulence during a space flight. The differences between flux distributions corresponding to normal and enhanced virulence states confirm previous observations concerning infection mechanisms and suggest new pathways for exploration.

A Genome-Scale Metabolic Model Accurately Predicts Fluxes in Central Carbon Metabolism under Stress Conditions

Flux is a key measure of the metabolic phenotype. Recently, complete (genome-scale) metabolic network models have been established for Arabidopsis (Arabidopsis thaliana), and flux distributions have been predicted using constraints-based modeling and optimization algorithms such as linear programming. While these models are useful for investigating possible flux states under different metabolic scenarios, it is not clear how close the predicted flux distributions are to those occurring in vivo. To address this, fluxes were predicted for heterotrophic Arabidopsis cells and compared with fluxes estimated in parallel by (13)C-metabolic flux analysis (MFA). Reactions of the central carbon metabolic network (glycolysis, the oxidative pentose phosphate pathway, and the tricarboxylic acid [TCA] cycle) were independently analyzed by the two approaches. Net fluxes in glycolysis and the TCA cycle were predicted accurately from the genome-scale model, whereas the oxidative pentose phosphate pathway was poorly predicted. MFA showed that increased temperature and hyperosmotic stress, which altered cell growth, also affected the intracellular flux distribution. Under both conditions, the genome-scale model was able to predict both the direction and magnitude of the changes in flux: namely, increased TCA cycle and decreased phosphoenolpyruvate carboxylase flux at high temperature and a general decrease in fluxes under hyperosmotic stress. MFA also revealed a 3-fold reduction in carbon-use efficiency at the higher temperature. It is concluded that constraints-based genome-scale modeling can be used to predict flux changes in central carbon metabolism under stress conditions.

Integration of Metabolic Modeling and Phenotypic Data in Evaluation and Improvement of Ethanol Production Using Respiration-Deficient Mutants of Saccharomyces cerevisiae

Flux balance analysis and phenotypic data were used to provide clues to the relationships between the activities of gene products and the phenotypes resulting from the deletion of genes involved in respiratory function in Saccharomyces cerevisiae. The effect of partial or complete respiratory deficiency on the ethanol production and growth characteristics of hap4Delta/hap4Delta, mig1Delta/mig1Delta, qdr3Delta/qdr3Delta, pdr3Delta/pdr3Delta, qcr7Delta/qcr7Delta, cyt1Delta/cyt1Delta, and rip1Delta/rip1Delta mutants grown in microaerated chemostats was investigated. The study provided additional evidence for the importance of the selection of a physiologically relevant objective function, and it may improve quantitative predictions of exchange fluxes, as well as qualitative estimations of changes in intracellular fluxes. Ethanol production was successfully predicted by flux balance analysis in the case of the qdr3Delta/qdr3Delta mutant, with maximization of ethanol production as the objective function, suggesting an additional role for Qdr3p in respiration. The absence of similar changes in estimated intracellular fluxes in the qcr7Delta/qcr7Delta mutant compared to the rip1Delta/rip1Delta and cyt1Delta/cyt1Delta mutants indicated that the effect of the deletion of this subunit of complex III was somehow compensated for. Analysis of predicted flux distributions indicated self-organization of intracellular fluxes to avoid NAD(+)/NADH imbalance in rip1Delta/rip1Delta and cyt1Delta/cyt1Delta mutants, but not the qcr7Delta/qcr7Delta mutant. The flux through the glycerol efflux channel, Fps1p, was estimated to be zero in all strains under the investigated conditions. This indicates that previous strategies for improving ethanol production, such as the overexpression of the glutamate synthase gene GLT1 in a GDH1 deletion background or deletion of the glycerol efflux channel gene FPS1 and overexpression of GLT1, are unnecessary in a respiration-deficient background.

Model-driven analysis of experimentally determined growth phenotypes for 465 yeast gene deletion mutants under 16 different conditions

An iterative approach that integrates high-throughput measurements of yeast deletion mutants and flux balance model predictions improves understanding of both experimental and computational results.

Including metabolite concentrations into flux balance analysis: thermodynamic realizability as a constraint on flux distributions in metabolic networks.

BackgroundIn recent years, constrained optimization – usually referred to as flux balance analysis (FBA) – has become a widely applied method for the computation of stationary fluxes in large-scale metabolic networks. The striking advantage of FBA as compared to kinetic modeling is that it basically requires only knowledge of the stoichiometry of the network. On the other hand, results of FBA are to a large degree hypothetical because the method relies on plausible but hardly provable optimality principles that are thought to govern metabolic flux distributions.ResultsTo augment the reliability of FBA-based flux calculations we propose an additional side constraint which assures thermodynamic realizability, i.e. that the flux directions are consistent with the corresponding changes of Gibb's free energies. The latter depend on metabolite levels for which plausible ranges can be inferred from experimental data. Computationally, our method results in the solution of a mixed integer linear optimization problem with quadratic scoring function. An optimal flux distribution together with a metabolite profile is determined which assures thermodynamic realizability with minimal deviations of metabolite levels from their expected values. We applied our novel approach to two exemplary metabolic networks of different complexity, the metabolic core network of erythrocytes (30 reactions) and the metabolic network iJR904 of Escherichia coli (931 reactions). Our calculations show that increasing network complexity entails increasing sensitivity of predicted flux distributions to variations of standard Gibb's free energy changes and metabolite concentration ranges. We demonstrate the usefulness of our method for assessing critical concentrations of external metabolites preventing attainment of a metabolic steady state.ConclusionOur method incorporates the thermodynamic link between flux directions and metabolite concentrations into a practical computational algorithm. The weakness of conventional FBA to rely on intuitive assumptions about the reversibility of biochemical reactions is overcome. This enables the computation of reliable flux distributions even under extreme conditions of the network (e.g. enzyme inhibition, depletion of substrates or accumulation of end products) where metabolite concentrations may be drastically altered.

The Determination of Angular Diameters of Stars

Ideally the determination of the angular diameter of a star would include the measurement of the distribution of intensity across the stellar disc. However, direct methods of measuring angular diameters have so far lacked adequate ‘signal to noise’ ratio to measure the intensity distribution and it has been the custom, in the first instance, to express the measured angular diameter in terms of the angular diameter of the equivalent uniform disc (θUD). Subsequent use of the angular diameter involves the assumption of a limb-darkening law and the application of an appropriate correction to θUD to find the ‘true’ angular diameter (θLD) of the star (e.g. Hanbury Brown et al., 1967). In this article we will discuss the determination of θUD for single stars and we will not refer further to the more difficult problems of determining intensity distributions involving limb-darkening and rotational effects and of measuring the angular parameters of binary systems.By itself the angular diameter of a star has no intrinsic value but when it is combined with other observational data it enables basic physical properties of the star to be determined. It is then possible to make a direct comparison of the observed properties of the star with the predictions of theoretical models of stellar atmospheres and interiors. For example, the combination of an angular diameter with the absolute monochromatic flux received from the star (ƒν), corrected for interstellar extinction, yields the absolute emergent flux at the stellar surface (). If the spectral energy distribution for the star is known it can be calibrated absolutely by and hence the effective temperature (Te) of the star can be found (this is equivalent to knowing the bolometric correction for the star and using it with the angular diameter to find Te). In addition to leading to the determination of Te, the absolute surface flux distribution may be compared directly with the predicted flux distributions for theoretical model stellar atmospheres (e.g. Davis and Webb, 1970). For O and early B. type stars a large fraction of the emergent flux is in the far ultra-violet and the effective temperatures cannot be determined from the, at present, incomplete empirical flux curves. In these cases it is possible to obtain an estimate of the effective temperatures by using the values of to calibrate a grid of model atmospheres which have Te as a parameter. In this way, by measuring the angular diameters of stars of different spectral types, it is possible to establish an effective temperature scale.