Search

Summary: Saccharomyces cerevisiae adjusts its metabolism based on nutrient availability, typically transitioning from glucose fermentation to ethanol respiration as glucose becomes limiting. However, our understanding of the regulation of metabolism is largely based on population averages, whereas nutrient transitions may cause heterogeneous responses. Here we introduce iCRAFT, a method that couples the ATP Förster resonance energy transfer (FRET)-based biosensor yAT1.03 with Antimycin A to differentiate fermentative and respiratory metabolisms in individual yeast cells. Upon Antimycin A addition, respiratory cells experienced a sharp decrease of the normalized FRET ratio, while respiro-fermentative cells showed no response. Next, we tracked changes in metabolism during the diauxic shift of a glucose pre-grown culture. Following glucose exhaustion, the entire cell population experienced a progressive rise in cytosolic ATP produced via respiration, suggesting a gradual increase in respiratory capacity. Overall, iCRAFT is a robust tool to distinguish fermentation from respiration, offering a new single-cell opportunity to study yeast metabolism.

ABSTRACT Bacteria adapt to nutrient availability by regulating the synthesis of enzymes. Transcriptome- and multi-sugar growth studies suggest that Lactococcus cremoris represses genes involved in the catabolization of lower growth rate-supporting (lower quality) sugars in a hierarchical order. Furthermore, L. cremoris appears to always express genes involved in the catabolization of higher growth rate-supporting sugars (higher quality) relative to the sugar it is growing on. Here, we unraveled the sugar catabolization hierarchy by determining the sugar catabolizing capacity and the proteome of cells exponentially growing on glucose (µmax = 0.72 h−1), lactose (µmax = 0.6 h−1), galactose (µmax = 0.44 h−1), or maltose (µmax = 0.43 h−1). We found that L. cremoris can grow on 14 of the 96 sugars in a Biolog plate, with µmax ranging from 0.32 to 0.72 h−1. Proteome and catabolization rate measurements show that L. cremoris consistently prepares for the catabolization of higher-quality sugars, except trehalose. While cells were not prepared for the catabolization of most lower-quality sugars, some proteins related to fructose and lactose consumption were always present. Moreover, reducing the growth rate of glucose through salt stress had only a minor influence on the sugars that L. cremoris could catabolize. These findings demonstrate that the catabolization hierarchy is not strictly linked to absolute growth rate or sugar quality. Cells instantly catabolizing a higher-quality sugar require enhanced expression of ribosomal and nucleotide metabolism functions for growth rate maximization, whereas transitioning to lower-quality sugars requires enhanced synthesis of proteins related to arginine catabolism and mixed acid fermentation, besides sugar-specific catabolic proteins. IMPORTANCE The availability of nutrients to microorganisms varies considerably between different environments, and changes can occur rapidly. As a general rule, a fast growth rate—typically growth on glucose—is associated with the repression of other carbohydrate utilization genes, but it is not clear to what extent catabolite repression is exerted by other sugars. We investigated the hierarchy of sugar utilization after substrate transitions in Lactococcus cremoris. For this, we determined the proteome and carbohydrate utilization capacity after growth on different sugars. The results show that the preparedness of cells for the utilization of “slower” sugars is not strictly determined by the growth rate. The data point to individual proteins relevant for various sugar transitions and suggest that the evolutionary history of the organism might be responsible for deviations from a strictly growth rate-related sugar catabolization hierarchy.

Metabolically active organelles compete for cytosolic space and resources during metabolism rewiring. Here, the authors develop a computational model of yeast metabolism and resource allocation to predict condition- and compartment-specific proteome constraints that govern metabolic strategies.

ABSTRACT The fission yeast, Schizosaccharomyces pombe, is a popular eukaryal model organism for cell division and cell cycle studies. With this extensive knowledge of its cell and molecular biology, S. pombe also holds promise for use in metabolism research and industrial applications. However, unlike the baker’s yeast, Saccharomyces cerevisiae, a major workhorse in these areas, cell physiology and metabolism of S. pombe remain less explored. One way to advance understanding of organism-specific metabolism is construction of computational models and their use for hypothesis testing. To this end, we leverage existing knowledge of S. cerevisiae to generate a manually curated high-quality reconstruction of S. pombe’s metabolic network, including a proteome-constrained version of the model. Using these models, we gain insights into the energy demands for growth, as well as ribosome kinetics in S. pombe. Furthermore, we predict proteome composition and identify growth-limiting constraints that determine optimal metabolic strategies under different glucose availability regimes and reproduce experimentally determined metabolic profiles. Notably, we find similarities in metabolic and proteome predictions of S. pombe with S. cerevisiae, which indicate that similar cellular resource constraints operate to dictate metabolic organization. With these cases, we show, on the one hand, how these models provide an efficient means to transfer metabolic knowledge from a well-studied to a lesser-studied organism, and on the other, how they can successfully be used to explore the metabolic behavior and the role of resource allocation in driving different strategies in fission yeast. IMPORTANCE Our understanding of microbial metabolism relies mostly on the knowledge we have obtained from a limited number of model organisms, and the diversity of metabolism beyond the handful of model species thus remains largely unexplored in mechanistic terms. Computational modeling of metabolic networks offers an attractive platform to bridge the knowledge gap and gain new insights into physiology of lesser-studied organisms. Here we showcase an example of successful knowledge transfer from the budding yeast Saccharomyces cerevisiae to a popular model organism in molecular and cell biology, fission yeast Schizosaccharomyces pombe, using computational models.

Abstract Cells adapt to different conditions via gene expression that tunes metabolism for maximal fitness. Constraints on cellular proteome may limit such expression strategies and introduce trade‐offs. Resource allocation under proteome constraints has explained regulatory strategies in bacteria. It is unclear, however, to what extent these constraints can predict evolutionary changes, especially for microorganisms that evolved under nutrient‐rich conditions, i.e., multiple available nitrogen sources, such as Lactococcus lactis. Here, we present a proteome‐constrained genome‐scale metabolic model of L. lactis (pcLactis) to interpret growth on multiple nutrients. Through integration of proteomics and flux data, in glucose‐limited chemostats, the model predicted glucose and arginine uptake as dominant constraints at low growth rates. Indeed, glucose and arginine catabolism were found upregulated in evolved mutants. At high growth rates, pcLactis correctly predicted the observed shutdown of arginine catabolism because limited proteome availability favored lactate for ATP production. Thus, our model‐based analysis is able to identify and explain the proteome constraints that limit growth rate in nutrient‐rich environments and thus form targets of fitness improvement.

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

ABSTRACT Yeasts constitute over 1,500 species with great potential for biotechnology. Still, the yeast Saccharomyces cerevisiae dominates industrial applications, and many alternative physiological capabilities of lesser-known yeasts are not being fully exploited. While comparative genomics receives substantial attention, little is known about yeasts’ metabolic specificity in batch cultures. Here, we propose a multiphase multiobjective dynamic genome-scale model of yeast batch cultures that describes the uptake of carbon and nitrogen sources and the production of primary and secondary metabolites. The model integrates a specific metabolic reconstruction, based on the consensus Yeast8, and a kinetic model describing the time-varying culture environment. In addition, we proposed a multiphase multiobjective flux balance analysis to compute the dynamics of intracellular fluxes. We then compared the metabolism of S. cerevisiae and Saccharomyces uvarum strains in a rich medium fermentation. The model successfully explained the experimental data and brought novel insights into how cryotolerant strains achieve redox balance. The proposed model (along with the corresponding code) provides a comprehensive picture of the main steps occurring inside the cell during batch cultures and offers a systematic approach to prospect or metabolically engineering novel yeast cell factories. IMPORTANCE Nonconventional yeast species hold the promise to provide novel metabolic routes to produce industrially relevant compounds and tolerate specific stressors, such as cold temperatures. This work validated the first multiphase multiobjective genome-scale dynamic model to describe carbon and nitrogen metabolism throughout batch fermentation. To test and illustrate its performance, we considered the comparative metabolism of three yeast strains of the Saccharomyces genus in rich medium fermentation. The study revealed that cryotolerant Saccharomyces species might use the γ-aminobutyric acid (GABA) shunt and the production of reducing equivalents as alternative routes to achieve redox balance, a novel biological insight worth being explored further. The proposed model (along with the provided code) can be applied to a wide range of batch processes started with different yeast species and media, offering a systematic and rational approach to prospect nonconventional yeast species metabolism and engineering novel cell factories.

Abstract Background Several genome-scale metabolic reconstruction software platforms have been developed and are being continuously updated. These tools have been widely applied to reconstruct metabolic models for hundreds of microorganisms ranging from important human pathogens to species of industrial relevance. However, these platforms, as yet, have not been systematically evaluated with respect to software quality, best potential uses and intrinsic capacity to generate high-quality, genome-scale metabolic models. It is therefore unclear for potential users which tool best fits the purpose of their research. Results In this work, we perform a systematic assessment of current genome-scale reconstruction software platforms. To meet our goal, we first define a list of features for assessing software quality related to genome-scale reconstruction. Subsequently, we use the feature list to evaluate the performance of each tool. To assess the similarity of the draft reconstructions to high-quality models, we compare each tool’s output networks with that of the high-quality, manually curated, models of Lactobacillus plantarum and Bordetella pertussis, representatives of gram-positive and gram-negative bacteria, respectively. We additionally compare draft reconstructions with a model of Pseudomonas putida to further confirm our findings. We show that none of the tools outperforms the others in all the defined features. Conclusions Model builders should carefully choose a tool (or combinations of tools) depending on the intended use of the metabolic model. They can use this benchmark study as a guide to select the best tool for their research. Finally, developers can also benefit from this evaluation by getting feedback to improve their software.

Abstract Background A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized. Results We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments. Conclusion This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Summary: The metabolic capabilities of cells determine their biotechnological potential, fitness in ecosystems, pathogenic threat levels, and function in multicellular organisms. Their comprehensive experimental characterization is generally not feasible, particularly for unculturable organisms. In principle, the full range of metabolic capabilities can be computed from an organism's annotated genome using metabolic network reconstruction. However, current computational methods cannot deal with genome-scale metabolic networks. Part of the problem is that these methods aim to enumerate all metabolic pathways, while computation of all (elementally balanced) conversions between nutrients and products would suffice. Indeed, the elementary conversion modes (ECMs, defined by Urbanczik and Wagner) capture the full metabolic capabilities of a network, but the use of ECMs has not been accessible until now. We explain and extend the theory of ECMs, implement their enumeration in ecmtool, and illustrate their applicability. This work contributes to the elucidation of the full metabolic footprint of any cell. The Bigger Picture: Understanding the metabolic capabilities of cells is of profound importance. Microbial metabolism shapes global cycles of elements and cleans polluted soils. Human and pathogen metabolism affects our health. Recent advances allow for automatic reconstruction of reaction networks for any organism, which is already used in synthetic biology, (food) microbiology, and agriculture to compute optimal yields from resources to products. However, computational tools are limited to optimal states or subnetworks, leaving many capabilities of organisms hidden. Our program, ecmtool, creates a blueprint of any organism's metabolic functionalities, drastically improving insights obtained from genome sequences. Ecmtool may become essential in exploratory research, especially for studying cells that are not culturable in laboratory conditions. Ideally, elementary conversion mode enumeration will someday be a standard step after metabolic network reconstruction, achieving the metabolic characterization of all known organisms.

Amino acids are attractive metabolites for the pharmaceutical and food industry field. On one hand, the construction of microbial cell factories for large-scale production aims to satisfy the demand for amino acids as bulk biochemical. On the other hand, amino acids enhance flavor formation in fermented foods. Concerning the latter, flavor formation in dairy products, such as cheese is associated with the presence of lactic acid bacteria (LAB). In particular, Lactococcus lactis, one of the most important LAB, is used as a starter culture in fermented foods. The proteolytic activity of some L. lactis strains results in peptides and amino acids, which are flavor compounds or flavor precursors. However, it is still a challenge to isolate bacterial cells with enhanced amino acid production and secretion activity. In this work, we developed a growth-based sensor strain to detect the essential amino acids isoleucine, leucine, valine, histidine and methionine. Amino acids are metabolites that can be secreted by some bacteria. Therefore, our biosensor allowed us to identify wild-type L. lactis strains that naturally secrete amino acids, by using co-cultures of the biosensor strain with potential amino acid producing strains. Subsequently, we used this biosensor in combination with a droplet-based screening approach, and isolated three mutated L. lactis IPLA838 strains with 5–10 fold increased amino acid-secretion compared to the wild type. Genome re-sequencing revealed mutations in genes encoding proteins that participate in peptide uptake and peptide degradation. We argue that an unbalance in the regulation of amino acid levels as a result of these gene mutations may drive the accumulation and secretion of these amino acids. This biosensing system tackles the problem of selection for overproduction of secreted molecules, which requires the coupling of the product to the producing cell in the droplets.

In this paper we try to describe all possible molecular states (phenotypes) for a cell that fabricates itself at a constant rate, given its enzyme kinetics and the stoichiometry of all reactions. For this, we must understand the process of cellular growth: steady-state self-fabrication requires a cell to synthesize all of its components, including metabolites, enzymes and ribosomes, in proportions that match its own composition. Simultaneously, the concentrations of these components affect the rates of metabolism and biosynthesis, and hence the growth rate. We here derive a theory that describes all phenotypes that solve this circular problem. All phenotypes can be described as a combination of minimal building blocks, which we call Elementary Growth Modes (EGMs). EGMs can be used as the theoretical basis for all models that explicitly model self-fabrication, such as the currently popular Metabolism and Expression models. We then use our theory to make concrete biological predictions. We find that natural selection for maximal growth rate drives microorganisms to states of minimal phenotypic complexity: only one EGM will be active when growth rate is maximised. The phenotype of a cell is only extended with one more EGM whenever growth becomes limited by an additional biophysical constraint, such as a limited solvent capacity of a cellular compartment. The theory presented here extends recent results on Elementary Flux Modes: the minimal building blocks of cellular growth models that lack the self-fabrication aspect. Our theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs.

Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.

Acetobacter pasteurianus 386B is a candidate functional starter culture for the cocoa bean fermentation process. To allow in silico simulations of its related metabolism in response to different environmental conditions, a genome-scale metabolic model for A. pasteurianus 386B was reconstructed. This is the first genome-scale metabolic model reconstruction for a member of the genus Acetobacter. The metabolic network reconstruction process was based on extensive genome re-annotation and comparative genomics analyses. The information content related to the functional annotation of metabolic enzymes and transporters was placed in a metabolic context by exploring and curating a Pathway/Genome Database of A. pasteurianus 386B using the Pathway Tools software. Metabolic reactions and curated gene-protein-reaction associations were bundled into a genome-scale metabolic model of A. pasteurianus 386B, named iAp386B454, containing 454 genes, 322 reactions, and 296 metabolites embedded in two cellular compartments. The reconstructed model was validated by performing growth experiments in a defined medium, which revealed that lactic acid as the sole carbon source could sustain growth of this strain. Further, the reconstruction of the A. pasteurianus 386B genome-scale metabolic model revealed knowledge gaps concerning the metabolism of this strain, especially related to the biosynthesis of its cell envelope and the presence or absence of metabolite transporters.

Microbial life usually takes place in a community where individuals interact, by competition for nutrients, cross-feeding, inhibition by end-products, but also by their spatial distribution. Lactic acid bacteria are prominent members of microbial communities responsible for food fermentations. Their niche in a community depends on their own properties as well as those of the other species. Here, we apply a computational approach, which uses only genomic and metagenomic information and functional annotation of genes, to find properties that distinguish a species from others in the community, as well as to follow individual species in a community. We analyzed isolated and sequenced strains from a kefir community, and metagenomes from wine fermentations. We demonstrate how the distinguishing properties of an organism lead to experimentally testable hypotheses concerning the niche and the interactions with other species. We observe, for example, that L. kefiranofaciens, a dominant organism in kefir, stands out among the Lactobacilli because it potentially has more amino acid auxotrophies. Using metagenomic analysis of industrial wine fermentations we investigate the role of an inoculated L. plantarum in malolactic fermentation. We observed that L. plantarum thrives better on white than on red wine fermentations and has the largest number of phosphotransferase system among the bacteria observed in the wine communities. Also, L. plantarum together with Pantoea, Erwinia, Asaia, Gluconobacter, and Komagataeibacter genera had the highest number of genes involved in biosynthesis of amino acids.

Growth rate is a near-universal selective pressure across microbial species. High growth rates require hundreds of metabolic enzymes, each with different nonlinear kinetics, to be precisely tuned within the bounds set by physicochemical constraints. Yet, the metabolic behaviour of many species is characterized by simple relations between growth rate, enzyme expression levels and metabolic rates. We asked if this simplicity could be the outcome of optimisation by evolution. Indeed, when the growth rate is maximized-in a static environment under mass-conservation and enzyme expression constraints-we prove mathematically that the resulting optimal metabolic flux distribution is described by a limited number of subnetworks, known as Elementary Flux Modes (EFMs). We show that, because EFMs are the minimal subnetworks leading to growth, a small active number automatically leads to the simple relations that are measured. We find that the maximal number of flux-carrying EFMs is determined only by the number of imposed constraints on enzyme expression, not by the size, kinetics or topology of the network. This minimal-EFM extremum principle is illustrated in a graphical framework, which explains qualitative changes in microbial behaviours, such as overflow metabolism and co-consumption, and provides a method for identification of the enzyme expression constraints that limit growth under the prevalent conditions. The extremum principle applies to all microorganisms that are selected for maximal growth rates under protein concentration constraints, for example the solvent capacities of cytosol, membrane or periplasmic space.

Genome-scale stoichiometric modeling methods, in particular Flux Balance Analysis (FBA) and variations thereof, are widely used to investigate cell metabolism and to optimize biotechnological processes. Given (1) a metabolic network, which can be reconstructed from an organism’s genome sequence, and (2) constraints on reaction rates, which may be based on measured nutrient uptake rates, FBA predicts which reactions maximize an objective flux, usually the production of cell components. Although FBA solutions may accurately predict the metabolic behavior of a cell, the actual flux predictions are often hard to interpret. This is especially the case for conditions with many constraints, such as for organisms growing in rich nutrient environments: it remains unclear why a certain solution was optimal. Here, we rationalize FBA solutions by explaining for which properties the optimal combination of metabolic strategies is selected. We provide a graphical formalism in which the selection of solutions can be visualized; we illustrate how this perspective provides a glimpse of the logic that underlies genome-scale modeling by applying our formalism to models of various sizes.

One of the marvels of biology is the phenotypic plasticity of microorganisms. It allows them to maintain high growth rates across conditions. Studies suggest that cells can express metabolic enzymes at tuned concentrations through adjustment of gene expression. The associated transcription factors are often regulated by intracellular metabolites. Here we study metabolite-mediated regulation of metabolic-gene expression that maximises metabolic fluxes across conditions. We developed an adaptive control theory, qORAC (for 'Specific Flux (q) Optimization by Robust Adaptive Control'), and illustrate it with several examples of metabolic pathways. The key feature of the theory is that it does not require knowledge of the regulatory network, only of the metabolic part. We derive that maximal metabolic flux can be maintained in the face of varying N environmental parameters only if the number of transcription-factor binding metabolites is at least equal to N. The controlling circuits appear to require simple biochemical kinetics. We conclude that microorganisms likely can achieve maximal rates in metabolic pathways, in the face of environmental changes.

Since the 1970s, galactose metabolism in Lactococcus lactis has been in debate. Different studies led to diverse outcomes making it difficult to conclude whether galactose uptake was PEP- or ATP- dependent and decide what the exact connection was between galactose and lactose uptake and metabolism. It was shown that some Lactococcus strains possess two galactose-specific systems – a permease and a PTS, even if they lack the lactose utilization plasmid, proving that a lactose-independent PTSGal exists. However, the PTSGal transporter was never identified. Here, with the help of transcriptome analyses and genetic knock-out mutants, we reveal the identities of two low-affinity galactose PTSs. A novel plant-niche-related PTS component Llmg_0963 forming a hybrid transporter Llmg_0963PtcBA and a glucose/mannose-specific PTS are shown to be involved in galactose transport in L. lactis MG1363.

Metabolism is generally required for cellular maintenance and for the generation of offspring under conditions that support growth. The rates, yields (efficiencies), adaptation time and robustness of metabolism are therefore key determinants of cellular fitness. For biotechnological applications and our understanding of the evolution of metabolism, it is necessary to figure out how the functional system properties of metabolism can be optimized, via adjustments of the kinetics and expression of enzymes, and by rewiring metabolism. The trade-offs that can occur during such optimizations then indicate fundamental limits to evolutionary innovations and bioengineering. In this paper, we review several theoretical and experimental findings about mechanisms for metabolic optimization.

In the model eukaryote Saccharomyces cerevisiae, it has long been known that a functional trehalose pathway is indispensable for transitions to high glucose conditions. Upon addition of glucose, cells with a defect in trehalose 6-phosphate synthase (Tps1), the first committed step in the trehalose pathway, display what we have termed an imbalanced glycolytic state; in this state the flux through the upper part of glycolysis outpaces that through the lower part of glycolysis. As a consequence, the intermediate fructose 1,6-bisphosphate (FBP) accumulates at low concentrations of ATP and inorganic phosphate (Pi). Despite significant research efforts, a satisfactory understanding of the regulatory role that trehalose metabolism plays during such transitions has remained infamously unresolved. In a recent study, we demonstrate that the startup of glycolysis exhibits two dynamic fates: a proper, functional, steady state or the imbalanced state described above. Both states are stable, attracting states, and the probability distribution of initial states determines the fate of a yeast cell exposed to glucose. Trehalose metabolism steers the dynamics of glycolysis towards the proper functional state through its ATP hydrolysis activity; a mechanism that ensures that the demand and supply of ATP is balanced with Pi availability under dynamic conditions. [van Heerden et al. Science (2014), DOI: 10.1126/science.1245114.]

An important challenge in microbial ecology is to infer metabolic-exchange fluxes between growing microbial species from community-level data, concerning species abundances and metabolite concentrations. Here we apply a model-based approach to integrate such experimental data and thereby infer metabolic-exchange fluxes. We designed a synthetic anaerobic co-culture of Clostridium acetobutylicum and Wolinella succinogenes that interact via interspecies hydrogen transfer and applied different environmental conditions for which we expected the metabolic-exchange rates to change. We used stoichiometric models of the metabolism of the two microorganisms that represents our current physiological understanding and found that this understanding - the model - is sufficient to infer the identity and magnitude of the metabolic-exchange fluxes and it suggested unexpected interactions. Where the model could not fit all experimental data, it indicates specific requirement for further physiological studies. We show that the nitrogen source influences the rate of interspecies hydrogen transfer in the co-culture. Additionally, the model can predict the intracellular fluxes and optimal metabolic exchange rates, which can point to engineering strategies. This study therefore offers a realistic illustration of the strengths and weaknesses of model-based integration of heterogenous data that makes inference of metabolic-exchange fluxes possible from community-level experimental data.

High-throughput data generation and genome-scale stoichiometric models have greatly facilitated the comprehensive study of metabolic networks. The computation of all feasible metabolic routes with these models, given stoichiometric, thermodynamic, and steady-state constraints, provides important insights into the metabolic capacities of a cell. How the feasible metabolic routes emerge from the interplay between flux constraints, optimality objectives, and the entire metabolic network of a cell is, however, only partially understood. We show how optimal metabolic routes, resulting from flux balance analysis computations, arise out of elementary flux modes, constraints, and optimization objectives. We illustrate our findings with a genome-scale stoichiometric model of Escherichia coli metabolism. In the case of one flux constraint, all feasible optimal flux routes can be derived from elementary flux modes alone. We found up to 120 million of such optimal elementary flux modes. We introduce a new computational method to compute the corner points of the optimal solution space fast and efficiently. Optimal flux routes no longer depend exclusively on elementary flux modes when we impose additional constraints; new optimal metabolic routes arise out of combinations of elementary flux modes. The solution space of feasible metabolic routes shrinks enormously when additional objectives---e.g. those related to pathway expression costs or pathway length---are introduced. In many cases, only a single metabolic route remains that is both feasible and optimal. This paper contributes to reaching a complete topological understanding of the metabolic capacity of organisms in terms of metabolic flux routes, one that is most natural to biochemists and biotechnologists studying and engineering metabolism.