Dynamic metabolic models in context: biomass backtracking

redGEM: Systematic reduction and analysis of genome-scale metabolic reconstructions for development of consistent core metabolic models

Genome-scale metabolic reconstructions have proven to be valuable resources in enhancing our understanding of metabolic networks as they encapsulate all known metabolic capabilities of the organisms from genes to proteins to their functions. However the complexity of these large metabolic networks often hinders their utility in various practical applications. Although reduced models are commonly used for modeling and in integrating experimental data, they are often inconsistent across different studies and laboratories due to different criteria and detail, which can compromise transferability of the findings and also integration of experimental data from different groups. In this study, we have developed a systematic semi-automatic approach to reduce genome-scale models into core models in a consistent and logical manner focusing on the central metabolism or subsystems of interest. The method minimizes the loss of information using an approach that combines graph-based search and optimization methods. The resulting core models are shown to be able to capture key properties of the genome-scale models and preserve consistency in terms of biomass and by-product yields, flux and concentration variability and gene essentiality. The development of these “consistently-reduced” models will help to clarify and facilitate integration of different experimental data to draw new understanding that can be directly extendable to genome-scale models.

Reduced models are used commonly to understand the metabolism of organisms and to integrate experimental data for many different studies such as physiology, fluxomics and metabolomics. Without consistent or clear criteria on how these reduced models are actually developed, it is difficult to ensure that they reflect the detailed knowledge that is kept in genome scale metabolic network models (GEMs). The redGEM algorithm presented here allows us to systematically develop consistently reduced metabolic models from their genome-scale counterparts. We applied redGEM for the construction of a core model for E. coli central carbon metabolism. We constructed the core model irJO1366 based on the latest genome-scale E. coli metabolic reconstruction (iJO1366). irJO1366 contains the central carbon pathways and other immediate pathways that must be connected to them for consistency with the iJO1366. irJO1366 can be used to understand metabolism of the organism and also to provide guidance for metabolic engineering purposes. The algorithm is also designed to be modular so that heterologous reactions or pathways can be appended to the core model akin to a “plug-and-play”, synthetic biology approach. The algorithm is applicable to any compartmentalized or non-compartmentalized GEM.

The MYpop toolbox: Putting yeast stress responses in cellular context on single cell and population scales

Systems biology holds the promise to integrate multiple sources of information in order to build ever more complete models of cellular function. To do this, the field must overcome two significant challenges. First, the current strategy to model average cells must be replaced with population based models accounting for cell-to-cell variability. Second, models must be integrated with each other and with basic cellular function. This requires a core model of cellular physiology as well as a multiscale simulation platform to support large-scale simulation of culture or tissues from single cells. Here, we present such a simulation platform with a core model of yeast physiology as scaffold to integrate and simulate SBML models. The software automates this integration helping users simulate their model of choice in context of the cell division cycle. We benchmark model merging, simulation and analysis by integrating a minimal model of osmotic stress into the core model and analyzing it. We characterize the effect of single cell differences on the dynamics of osmoadaptation, estimating when normal cell growth is resumed and obtaining an explanation for experimentally observed glycerol dynamics based on population dynamics. Hence, the platform can be used to reconcile single cell and population level data.