Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

Understanding the adaptive laboratory evolution of multiple stress-resistant yeast strains by genome-scale modeling.

Stress responses triggered by external exposures in adaptive laboratory evolution studies alter the ordinary behavior of cells, and the identification of the differences between the starting and the evolved strains would provide ideal strategies to obtain the desired strains. Metabolic networks are one of the most useful tools to analyze data for this purpose. This study integrates differential expression profiles of multiple Saccharomyces cerevisiae strains that have evolved in eight different stress conditions (ethanol, caffeine, coniferyl aldehyde, iron, nickel, phenylethanol, and silver) and enzyme kinetics into a genome-scale metabolic model of yeast, following a new enhanced method. Flux balance analysis, flux variability analysis, robustness, phenotype phase plane, minimization of metabolic adjustment, survivability, sensitivity analyses, and random sampling are conducted to identify the most common and divergent points within strains. Results were examined both individually and comparatively, and the target reactions, metabolites, and enzymes were identified. Our results showed that the models reconstructed by our methodology were able to simulate experimental conditions where efficient protein allocation was the main goal for survival under stressful conditions, and most of the metabolic changes in the adaptation process mainly arose from the differences in the metabolic reactions of energy maintenance (through coenzyme-A and FAD utilization), cell division (folate requirement of DNA synthesis), and cell wall formation (through sterol and ergosterol biosynthesis).