Predicting microbial interactions with approaches based on flux balance analysis: an evaluation.

BACKGROUND: Given a genome-scale metabolic model (GEM) of a microorganism and criteria for optimization, flux balance analysis (FBA) predicts the optimal growth rate and its corresponding flux distribution for a specific medium. FBA has been extended to microbial consortia and thus can be used to predict interactions by comparing in-silico growth rates for co- and monocultures. Although FBA-based methods for microbial interaction prediction are becoming popular, a systematic evaluation of their accuracy has not yet been performed. RESULTS: Here, we evaluate the accuracy of FBA-based predictions of human and mouse gut bacterial interactions using growth data from the literature. For this, we collected 26 GEMs from the semi-curated AGORA database as well as four previously published curated GEMs. We tested the accuracy of three tools (COMETS, Microbiome Modeling Toolbox and MICOM) by comparing growth rates predicted in mono- and co-culture to growth rates extracted from the literature and also investigated the impact of different tool settings and media. We found that except for curated GEMs, predicted growth rates and their ratios (i.e. interaction strengths) do not correlate with growth rates and interaction strengths obtained from in vitro data. CONCLUSIONS: Prediction of growth rates with FBA using semi-curated GEMs is currently not sufficiently accurate to predict interaction strengths reliably.

A strategy to detect metabolic changes induced by exposure to chemicals from large sets of condition-specific metabolic models computed with enumeration techniques.

BACKGROUND: The growing abundance of in vitro omics data, coupled with the necessity to reduce animal testing in the safety assessment of chemical compounds and even eliminate it in the evaluation of cosmetics, highlights the need for adequate computational methodologies. Data from omics technologies allow the exploration of a wide range of biological processes, therefore providing a better understanding of mechanisms of action (MoA) related to chemical exposure in biological systems. However, the analysis of these large datasets remains difficult due to the complexity of modulations spanning multiple biological processes. RESULTS: To address this, we propose a strategy to reduce information overload by computing, based on transcriptomics data, a comprehensive metabolic sub-network reflecting the metabolic impact of a chemical. The proposed strategy integrates transcriptomic data to a genome scale metabolic network through enumeration of condition-specific metabolic models hence translating transcriptomics data into reaction activity probabilities. Based on these results, a graph algorithm is applied to retrieve user readable sub-networks reflecting the possible metabolic MoA (mMoA) of chemicals. This strategy has been implemented as a three-step workflow. The first step consists in building cell condition-specific models reflecting the metabolic impact of each exposure condition while taking into account the diversity of possible optimal solutions with a partial enumeration algorithm. In a second step, we address the challenge of analyzing thousands of enumerated condition-specific networks by computing differentially activated reactions (DARs) between the two sets of enumerated possible condition-specific models. Finally, in the third step, DARs are grouped into clusters of functionally interconnected metabolic reactions, representing possible mMoA, using the distance-based clustering and subnetwork extraction method. The first part of the workflow was exemplified on eight molecules selected for their known human hepatotoxic outcomes associated with specific MoAs well described in the literature and for which we retrieved primary human hepatocytes transcriptomic data in Open TG-GATEs. Then, we further applied this strategy to more precisely model and visualize associated mMoA for two of these eight molecules (amiodarone and valproic acid). The approach proved to go beyond gene-based analysis by identifying mMoA when few genes are significantly differentially expressed (2 differentially expressed genes (DEGs) for amiodarone), bringing additional information from the network topology, or when very large number of genes were differentially expressed (5709 DEGs for valproic acid). In both cases, the results of our strategy well fitted evidence from the literature regarding known MoA. Beyond these confirmations, the workflow highlighted potential other unexplored mMoA. CONCLUSION: The proposed strategy allows toxicology experts to decipher which part of cellular metabolism is expected to be affected by the exposition to a given chemical. The approach originality resides in the combination of different metabolic modelling approaches (constraint based and graph modelling). The application to two model molecules shows the strong potential of the approach for interpretation and visual mining of complex omics in vitro data. The presented strategy is freely available as a python module ( https://pypi.org/project/manamodeller/ ) and jupyter notebooks ( https://github.com/LouisonF/MANA ).

Bdellovibrio's prey-independent lifestyle is fueled by amino acids as a carbon source.

Identifying the nutritional requirements and growth conditions of microorganisms is crucial for determining their applicability in industry and understanding their role in clinical ecology. Predatory bacteria such as Bdellovibrio bacteriovorus have emerged as promising tools for combating infections by human bacterial pathogens due to their natural killing features. Bdellovibrio's lifecycle occurs inside prey cells, using the cytoplasm as a source of nutrients and energy. However, this lifecycle supposes a challenge when determining the specific uptake of metabolites from the prey to complete the growth inside cells, a process that has not been completely elucidated. Here, following a model-based approach, we illuminate the ability of B. bacteriovorus to replicate DNA, increase biomass, and generate adenosine triphosphate (ATP) in an amino acid-based rich media in the absence of prey, keeping intact its predatory capacity. In this culture, we determined the main carbon sources used and their preference, being glutamate, serine, aspartate, isoleucine, and threonine. This study offers new insights into the role of predatory bacteria in natural environments and establishes the basis for developing new Bdellovibrio applications using appropriate metabolic and physiological methodologies. KEY POINTS: • Amino acids support axenic lifestyle of Bdellovibrio bacteriovorus. • B. bacteriovorus preserves its predatory ability when growing in the absence of prey.

Evaluating the feasibility of medium-chain oleochemical synthesis using microbial chain elongation.

Chain elongating bacteria are a unique guild of strictly anaerobic bacteria that have garnered interest for sustainable chemical manufacturing from carbon-rich wet and gaseous waste streams. They produce C6-C8 medium-chain fatty acids, which are valuable platform chemicals that can be used directly, or derivatized to service a wide range of chemical industries. However, the application of chain elongating bacteria for synthesizing products beyond C6-C8 medium-chain fatty acids has not been evaluated. In this study, we assess the feasibility of expanding the product spectrum of chain elongating bacteria to C9-C12 fatty acids, along with the synthesis of C6 fatty alcohols, dicarboxylic acids, diols, and methyl ketones. We propose several metabolic engineering strategies to accomplish these conversions in chain elongating bacteria and utilize constraint-based metabolic modelling to predict pathway stoichiometries, assess thermodynamic feasibility, and estimate ATP and product yields. We also evaluate how producing alternative products impacts the growth rate of chain elongating bacteria via resource allocation modelling, revealing a trade-off between product chain length and class versus cell growth rate. Together, these results highlight the potential for using chain elongating bacteria as a platform for diverse oleochemical biomanufacturing and offer a starting point for guiding future metabolic engineering efforts aimed at expanding their product range. ONE-SENTENCE SUMMARY: In this work, the authors use constraint-based metabolic modelling and enzyme cost minimization to assess the feasibility of using metabolic engineering to expand the product spectrum of anaerobic chain elongating bacteria.

Genome-scale metabolic models reveal determinants of phenotypic differences in non-Saccharomyces yeasts.

BACKGROUND: Use of alternative non-Saccharomyces yeasts in wine and beer brewing has gained more attention the recent years. This is both due to the desire to obtain a wider variety of flavours in the product and to reduce the final alcohol content. Given the metabolic differences between the yeast species, we wanted to account for some of the differences by using in silico models. RESULTS: We created and studied genome-scale metabolic models of five different non-Saccharomyces species using an automated processes. These were: Metschnikowia pulcherrima, Lachancea thermotolerans, Hanseniaspora osmophila, Torulaspora delbrueckii and Kluyveromyces lactis. Using the models, we predicted that M. pulcherrima, when compared to the other species, conducts more respiration and thus produces less fermentation products, a finding which agrees with experimental data. Complex I of the electron transport chain was to be present in M. pulcherrima, but absent in the others. The predicted importance of Complex I was diminished when we incorporated constraints on the amount of enzymatic protein, as this shifts the metabolism towards fermentation. CONCLUSIONS: Our results suggest that Complex I in the electron transport chain is a key differentiator between Metschnikowia pulcherrima and the other yeasts considered. Yet, more annotations and experimental data have the potential to improve model quality in order to increase fidelity and confidence in these results. Further experiments should be conducted to confirm the in vivo effect of Complex I in M. pulcherrima and its respiratory metabolism.

Quantifying beta cell function in the preclinical stages of type 1 diabetes.

Clinically symptomatic type 1 diabetes (stage 3 type 1 diabetes) is preceded by a pre-symptomatic phase, characterised by progressive loss of functional beta cell mass after the onset of islet autoimmunity, with (stage 2) or without (stage 1) measurable changes in glucose profile during an OGTT. Identifying metabolic tests that can longitudinally track changes in beta cell function is of pivotal importance to track disease progression and measure the effect of disease-modifying interventions. In this review we describe the metabolic changes that occur in the early pre-symptomatic stages of type 1 diabetes with respect to both insulin secretion and insulin sensitivity, as well as the measurable outcomes that can be derived from the available tests. We also discuss the use of metabolic modelling to identify insulin secretion and sensitivity, and the measurable changes during dynamic tests such as the OGTT. Finally, we review the role of risk indices and minimally invasive measures such as those derived from the use of continuous glucose monitoring.

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

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

Deep learning meets metabolomics: a methodological perspective.

Deep learning (DL), an emerging area of investigation in the fields of machine learning and artificial intelligence, has markedly advanced over the past years. DL techniques are being applied to assist medical professionals and researchers in improving clinical diagnosis, disease prediction and drug discovery. It is expected that DL will help to provide actionable knowledge from a variety of 'big data', including metabolomics data. In this review, we discuss the applicability of DL to metabolomics, while presenting and discussing several examples from recent research. We emphasize the use of DL in tackling bottlenecks in metabolomics data acquisition, processing, metabolite identification, as well as in metabolic phenotyping and biomarker discovery. Finally, we discuss how DL is used in genome-scale metabolic modelling and in interpretation of metabolomics data. The DL-based approaches discussed here may assist computational biologists with the integration, prediction and drawing of statistical inference about biological outcomes, based on metabolomics data.

Optogenetic control of Cdc48 for dynamic metabolic engineering in yeast.

Dynamic metabolic engineering is a strategy to switch key metabolic pathways in microbial cell factories from biomass generation to accumulation of target products. Here, we demonstrate that optogenetic intervention in the cell cycle of budding yeast can be used to increase production of valuable chemicals, such as the terpenoid ß-carotene or the nucleoside analog cordycepin. We achieved optogenetic cell-cycle arrest in the G2/M phase by controlling activity of the ubiquitin-proteasome system hub Cdc48. To analyze the metabolic capacities in the cell cycle arrested yeast strain, we studied their proteomes by timsTOF mass spectrometry. This revealed widespread, but highly distinct abundance changes of metabolic key enzymes. Integration of the proteomics data in protein-constrained metabolic models demonstrated modulation of fluxes directly associated with terpenoid production as well as metabolic subsystems involved in protein biosynthesis, cell wall synthesis, and cofactor biosynthesis. These results demonstrate that optogenetically triggered cell cycle intervention is an option to increase the yields of compounds synthesized in a cellular factory by reallocation of metabolic resources.

Multi-dimensional experimental and computational exploration of metabolism pinpoints complex probiotic interactions.

Multi-strain probiotics are widely regarded as effective products for improving gut microbiota stability and host health, providing advantages over single-strain probiotics. However, in general, it is unclear to what extent different strains would cooperate or compete for resources, and how the establishment of a common biofilm microenvironment could influence their interactions. In this work, we develop an integrative experimental and computational approach to comprehensively assess the metabolic functionality and interactions of probiotics across growth conditions. Our approach combines co-culture assays with genome-scale modelling of metabolism and multivariate data analysis, thus exploiting complementary data- and knowledge-driven systems biology techniques. To show the advantages of the proposed approach, we apply it to the study of the interactions between two widely used probiotic strains of Lactobacillus reuteri and Saccharomyces boulardii, characterising their production potential for compounds that can be beneficial to human health. Our results show that these strains can establish a mixed cooperative-antagonistic interaction best explained by competition for shared resources, with an increased individual exchange but an often decreased net production of amino acids and short-chain fatty acids. Overall, our work provides a strategy that can be used to explore microbial metabolic fingerprints of biotechnological interest, capable of capturing multifaceted equilibria even in simple microbial consortia.

Structural variability, implementational irregularities in mathematical modelling of polyhydroxyalkanoates (PHAs) production-A state-of-the-art review.

The rate and extent of microbial polyhydroxyalkanoates (PHAs) production rely on the availability of substrates, growth of microbial biomass, and intracellular accumulation of polymer under nitrogen-limited conditions. The dynamics of PHAs production captured through various structured or unstructured models can be extended to design an optimal feeding strategy for process intensification. Large variability in process assumptions, choices of kinetics, and model complexity is expected depending on substrate(s), microbial metabolism, and discretization of the process under consideration. This communication attempts to review the estimation of stoichiometric yield coefficients, metabolic modelling, and choices of unstructured kinetics in microbial PHA production. Implementational irregularities in parameter estimation and quality check in modelling exercises have also been reviewed. It is observed that the scope of the majority of the "modelling" studies is confined to the estimation of stoichiometric parameters with limited utility. In dynamic models, microbial growth is often described using either Monod or logistic variants, while PHAs production adopts a Luedeking-Piret expression with or without substrate inhibition. Though model selection, regression with experimental data, parameter estimation, and model validation are integral parts of the exercise, very few provide sufficient coverage on all those aspects. Application of the model to control or optimize the bioprocess has rarely been attempted.

EFMlrs: a Python package for elementary flux mode enumeration via lexicographic reverse search.

BACKGROUND: Elementary flux mode (EFM) analysis is a well-established, yet computationally challenging approach to characterize metabolic networks. Standard algorithms require huge amounts of memory and lack scalability which limits their application to single servers and consequently limits a comprehensive analysis to medium-scale networks. Recently, Avis et al. developed mplrs-a parallel version of the lexicographic reverse search (lrs) algorithm, which, in principle, enables an EFM analysis on high-performance computing environments (Avis and Jordan. mplrs: a scalable parallel vertex/facet enumeration code. arXiv:1511.06487 , 2017). Here we test its applicability for EFM enumeration. RESULTS: We developed EFMlrs, a Python package that gives users access to the enumeration capabilities of mplrs. EFMlrs uses COBRApy to process metabolic models from sbml files, performs loss-free compressions of the stoichiometric matrix, and generates suitable inputs for mplrs as well as efmtool, providing support not only for our proposed new method for EFM enumeration but also for already established tools. By leveraging COBRApy, EFMlrs also allows the application of additional reaction boundaries and seamlessly integrates into existing workflows. CONCLUSION: We show that due to mplrs's properties, the algorithm is perfectly suited for high-performance computing (HPC) and thus offers new possibilities for the unbiased analysis of substantially larger metabolic models via EFM analyses. EFMlrs is an open-source program that comes together with a designated workflow and can be easily installed via pip.

Understanding and mathematical modelling of cellular resource allocation in microorganisms: a comparative synthesis.

BACKGROUND: The rising consensus that the cell can dynamically allocate its resources provides an interesting angle for discovering the governing principles of cell growth and metabolism. Extensive efforts have been made in the past decade to elucidate the relationship between resource allocation and phenotypic patterns of microorganisms. Despite these exciting developments, there is still a lack of explicit comparison between potentially competing propositions and a lack of synthesis of inter-related proposals and findings. RESULTS: In this work, we have reviewed resource allocation-derived principles, hypotheses and mathematical models to recapitulate important achievements in this area. In particular, the emergence of resource allocation phenomena is deciphered by the putative tug of war between the cellular objectives, demands and the supply capability. Competing hypotheses for explaining the most-studied phenomenon arising from resource allocation, i.e. the overflow metabolism, have been re-examined towards uncovering the potential physiological root cause. The possible link between proteome fractions and the partition of the ribosomal machinery has been analysed through mathematical derivations. Finally, open questions are highlighted and an outlook on the practical applications is provided. It is the authors' intention that this review contributes to a clearer understanding of the role of resource allocation in resolving bacterial growth strategies, one of the central questions in microbiology. CONCLUSIONS: We have shown the importance of resource allocation in understanding various aspects of cellular systems. Several important questions such as the physiological root cause of overflow metabolism and the correct interpretation of 'protein costs' are shown to remain open. As the understanding of the mechanisms and utility of resource application in cellular systems further develops, we anticipate that mathematical modelling tools incorporating resource allocation will facilitate the circuit-host design in synthetic biology.

A quantitative metabolic analysis reveals Acetobacterium woodii as a flexible and robust host for formate-based bioproduction.

Cheap and renewable feedstocks such as the one-carbon substrate formate are emerging for sustainable production in a growing chemical industry. We investigated the acetogen Acetobacterium woodii as a potential host for bioproduction from formate alone and together with autotrophic and heterotrophic co-substrates by quantitatively analyzing physiology, transcriptome, and proteome in chemostat cultivations in combination with computational analyses. Continuous cultivations with a specific growth rate of 0.05 h-1 on formate showed high specific substrate uptake rates (47 mmol g-1 h-1). Co-utilization of formate with H2, CO, CO2 or fructose was achieved without catabolite repression and with acetate as the sole metabolic product. A transcriptomic comparison of all growth conditions revealed a distinct adaptation of A. woodii to growth on formate as 570 genes were changed in their transcript level. Transcriptome and proteome showed higher expression of the Wood-Ljungdahl pathway during growth on formate and gaseous substrates, underlining its function during utilization of one-carbon substrates. Flux balance analysis showed varying flux levels for the WLP (0.7-16.4 mmol g-1 h-1) and major differences in redox and energy metabolism. Growth on formate, H2/CO2, and formate + H2/CO2 resulted in low energy availability (0.20-0.22 ATP/acetate) which was increased during co-utilization with CO or fructose (0.31 ATP/acetate for formate + H2/CO/CO2, 0.75 ATP/acetate for formate + fructose). Unitrophic and mixotrophic conversion of all substrates was further characterized by high energetic efficiencies. In silico analysis of bioproduction of ethanol and lactate from formate and autotrophic and heterotrophic co-substrates showed promising energetic efficiencies (70-92%). Collectively, our findings reveal A. woodii as a promising host for flexible and simultaneous bioconversion of multiple substrates, underline the potential of substrate co-utilization to improve the energy availability of acetogens and encourage metabolic engineering of acetogenic bacteria for the efficient synthesis of bulk chemicals and fuels from sustainable one carbon substrates.

Diversity begets diversity during community assembly until ecological limits impose a diversity ceiling.

Microbial communities are hugely diverse, but we do not yet understand how species invasions and extinctions drive and limit their diversity. On the one hand, the ecological limits hypothesis posits that diversity is primarily limited by environmental resources. On the other hand, the diversity-begets-diversity hypothesis posits that such limits can be easily lifted when new ecological niches are created by biotic interactions. To find out which hypothesis better explains the assembly of microbial communities, we used metabolic modelling. We represent each microbial species by a metabolic network that harbours thousands of biochemical reactions. Together, these reactions determine which carbon and energy sources a species can use, and which metabolic by-products-potential nutrients for other species-it can excrete in a given environment. We assemble communities by modelling thousands of species invasions in a chemostat-like environment. We find that early during the assembly process, diversity begets diversity. By-product excretion transforms a simple environment into one that can sustain dozens of species. During later assembly stages, the creation of new niches slows down, existing niches become filled, successful invasions become rare, and species diversity plateaus. Thus, ecological limitations dominate the late assembly process. We conclude that each hypothesis captures a different stage of the assembly process. Species interactions can raise a community's diversity ceiling dramatically, but only within limits imposed by the environment.

Metabolic flux from the chloroplast provides signals controlling photosynthetic acclimation to cold in Arabidopsis thaliana.

Photosynthesis is especially sensitive to environmental conditions, and the composition of the photosynthetic apparatus can be modulated in response to environmental change, a process termed photosynthetic acclimation. Previously, we identified a role for a cytosolic fumarase, FUM2 in acclimation to low temperature in Arabidopsis thaliana. Mutant lines lacking FUM2 were unable to acclimate their photosynthetic apparatus to cold. Here, using gas exchange measurements and metabolite assays of acclimating and non-acclimating plants, we show that acclimation to low temperature results in a change in the distribution of photosynthetically fixed carbon to different storage pools during the day. Proteomic analysis of wild-type Col-0 Arabidopsis and of a fum2 mutant, which was unable to acclimate to cold, indicates that extensive changes occurring in response to cold are affected in the mutant. Metabolic and proteomic data were used to parameterize metabolic models. Using an approach called flux sampling, we show how the relative export of triose phosphate and 3-phosphoglycerate provides a signal of the chloroplast redox state that could underlie photosynthetic acclimation to cold.

Metabolic and molecular modelling of zebrafish gut biome to unravel antimicrobial peptides through metagenomics.

Recently efforts have been taken for unravelling mysteries between host-microbe interactions in gut microbiome studies of model organisms through metagenomics. Co-existence and the co-evolution of the microorganisms is the significant cause of the growing antimicrobial menace. There needs a novel approach to develop potential antimicrobials with capabilities to act directly on the resistant microbes with reduced side effects. One such is to tap them from the natural resources, preferably the gut of the most closely related animal model. In this study, we employed metagenomics approaches to identify the large taxonomic genomes of the zebra fish gut. About 256 antimicrobial peptides were identified using gene ontology predictions from Macrel and Pubseed servers. Upon the property predictions, the top 10 antimicrobial peptides were screened based on their action against many resistant bacterial species, including Klebsiella pneumoniae, Pseudomonas aeruginosa, Staphylococcus aureus, E. coli, and Bacillus cereus. Metabolic modelling and flux balance analysis (FBA) were computed to conclude the antibiotic such as tetracycline, cephalosporins, puromycin, neomycin biosynthesis pathways were adopted by the microbiome as protection strategies. Molecular modelling strategies, including molecular docking and dynamics, were performed to estimate the antimicrobial peptides' binding against the target-putative nucleic acid binding lipoprotein and confirm stable binding. One specific antimicrobial peptide with the sequence "MPPYLHEIQPHTASNCQTELVIKL" showed promising results with 53% hydrophobic residues and a net charge +2.5, significant for the development of antimicrobial peptides. The said peptide also showed promising interactions with the target protein and expressed stable binding with docking energy of -429.34 kcal/mol and the average root mean square deviation of 1 A0. The study is a novel approach focusing on tapping out potential antimicrobial peptides to be developed against most resistant bacterial species.

Recombinant protein production in Pichia pastoris: from transcriptionally redesigned strains to bioprocess optimization and metabolic modelling.

Pichia pastoris is one of the most widely used host for the production of recombinant proteins. Expression systems that rely mostly on promoters from genes encoding alcohol oxidase 1 or glyceraldehyde-3-phosphate dehydrogenase have been developed together with related bioreactor operation strategies based on carbon sources such as methanol, glycerol, or glucose. Although, these processes are relatively efficient and easy to use, there have been notable improvements over the last twenty years to better control gene expression from these promoters and their engineered variants. Methanol-free and more efficient protein production platforms have been developed by engineering promoters and transcription factors. The production window of P. pastoris has been also extended by using alternative feedstocks including ethanol, lactic acid, mannitol, sorbitol, sucrose, xylose, gluconate, formate or rhamnose. Herein, the specific aspects that are emerging as key parameters for recombinant protein synthesis are discussed. For this purpose, a holistic approach has been considered to scrutinize protein production processes from strain design to bioprocess optimization, particularly focusing on promoter engineering, transcriptional circuitry redesign. This review also considers the optimization of bioprocess based on alternative carbon sources and derived co-feeding strategies. Optimization strategies for recombinant protein synthesis through metabolic modelling are also discussed.

Dynamic co-culture metabolic models reveal the fermentation dynamics, metabolic capacities and interplays of cheese starter cultures.

In this study, we have investigated the cheese starter culture as a microbial community through a question: can the metabolic behaviour of a co-culture be explained by the characterized individual organism that constituted the co-culture? To address this question, the dairy-origin lactic acid bacteria Lactococcus lactis subsp. cremoris, Lactococcus lactis subsp. lactis, Streptococcus thermophilus and Leuconostoc mesenteroides, commonly used in cheese starter cultures, were grown in pure and four different co-cultures. We used a dynamic metabolic modelling approach based on the integration of the genome-scale metabolic networks of the involved organisms to simulate the co-cultures. The strain-specific kinetic parameters of dynamic models were estimated using the pure culture experiments and they were subsequently applied to co-culture models. Biomass, carbon source, lactic acid and most of the amino acid concentration profiles simulated by the co-culture models fit closely to the experimental results and the co-culture models explained the mechanisms behind the dynamic microbial abundance. We then applied the co-culture models to estimate further information on the co-cultures that could not be obtained by the experimental method used. This includes estimation of the profile of various metabolites in the co-culture medium such as flavour compounds produced and the individual organism level metabolic exchange flux profiles, which revealed the potential metabolic interactions between organisms in the co-cultures.

Multi-tissue to whole plant metabolic modelling.

Genome-scale metabolic models have been successfully applied to study the metabolism of multiple plant species in the past decade. While most existing genome-scale modelling studies have focussed on studying the metabolic behaviour of individual plant metabolic systems, there is an increasing focus on combining models of multiple tissues or organs to produce multi-tissue models that allow the investigation of metabolic interactions between tissues and organs. Multi-tissue metabolic models were constructed for multiple plants including Arabidopsis, barley, soybean and Setaria. These models were applied to study various aspects of plant physiology including the division of labour between organs, source and sink tissue relationship, growth of different tissues and organs and charge and proton balancing. In this review, we outline the process of constructing multi-tissue genome-scale metabolic models, discuss the strengths and challenges in using multi-tissue models, review the current status of plant multi-tissue and whole plant metabolic models and explore the approaches for integrating genome-scale metabolic models into multi-scale plant models.