Streptomyces coelicolor is a model organism for the Actinobacteria, a phylum known to produce an extensive range of different bioactive compounds that include antibiotics currently used in the clinic. Biosynthetic gene clusters discovered in genomes of other Actinobacteria can be transferred to and expressed in S. coelicolor, making it a factory for heterologous production of secondary metabolites. Genome‐scale metabolic reconstructions have successfully been used in several biotechnology applications to facilitate the over‐production of target metabolites. Here, the authors present iKS1317, the most comprehensive and accurate reconstructed genome‐scale metabolic model (GEM) for S. coelicolor. The model reconstruction is based on previous models, publicly available databases, and published literature and includes 1317 genes, 2119 reactions, and 1581 metabolites. It correctly predicts wild‐type growth in 96.5% of the evaluated growth environments and gene knockout predictions in 78.4% when comparing with observed mutant growth phenotypes, with a total accuracy of 83.3%. However, using a minimal nutrient environment for the gene knockout predictions, iKS1317 has an accuracy of 87.1% in predicting mutant growth phenotypes. Furthermore, we used iKS1317 and existing strain design algorithms to suggest robust gene‐knockout strategies to increase the production of acetyl‐CoA. Since acetyl‐CoA is the most important precursor for polyketide antibiotics, the suggested strategies may be implemented in vivo to improve the function of S. coelicolor as a heterologous expression host.

Summary Many biosynthetic gene clusters (BGCs) require heterologous expression to realize their genetic potential, including silent and metagenomic BGCs. Although the engineered Streptomyces coelicolor M1152 is a widely used host for heterologous expression of BGCs, a systemic understanding of how its genetic modifications affect the metabolism is lacking and limiting further development. We performed a comparative analysis of M1152 and its ancestor M145, connecting information from proteomics, transcriptomics, and cultivation data into a comprehensive picture of the metabolic differences between these strains. Instrumental to this comparison was the application of an improved consensus genome-scale metabolic model (GEM) of S. coelicolor . Although many metabolic patterns are retained in M1152, we find that this strain suffers from oxidative stress, possibly caused by increased oxidative metabolism. Furthermore, precursor availability is likely not limiting polyketide production, implying that other strategies could be beneficial for further development of S. coelicolor for heterologous production of novel compounds.

Genome-scale metabolic modeling is an important tool in the study of metabolism by enhancing the collation of knowledge, interpretation of data, and prediction of metabolic capabilities. A frequent assumption in the use of genome-scale models is that the in vivo organism is evolved for optimal growth, where growth is represented by flux through a biomass objective function (BOF). While the specific composition of the BOF is crucial, its formulation is often inherited from similar organisms due to the experimental challenges associated with its proper determination. A cell’s macro-molecular composition is not fixed and it responds to changes in environmental conditions. As a consequence, initiatives for the high-fidelity determination of cellular biomass composition have been launched. Thus, there is a need for a mathematical and computational framework capable of using multiple measurements of cellular biomass composition in different environments. Here, we propose two different computational approaches for directly addressing this challenge: Biomass Trade-off Weighting (BTW) and Higher-dimensional-plane InterPolation (HIP). In lieu of experimental data on biomass composition-variation in response to changing nutrient environment, we assess the properties of BTW and HIP using three hypothetical, yet biologically plausible, BOFs for the Escherichia coli genome-scale metabolic model iML1515. We find that the BTW and HIP formulations have a significant impact on model performance and phenotypes. Furthermore, the BTW method generates larger growth rates in all environments when compared to HIP. Using acetate secretion and the respiratory quotient as proxies for phenotypic changes, we find marked differences between the methods as HIP generates BOFs more similar to a reference BOF than BTW. We conclude that the presented methods constitute a conceptual step in developing genome-scale metabolic modelling approaches capable of addressing the inherent dependence of cellular biomass composition on nutrient environments.

Streptomyces coelicolor M1152 is a widely used host strain for the heterologous production of novel small molecule natural products, genetically engineered for this purpose through e.g. deletion of four of its native biosynthetic gene clusters (BGCs) for improved precursor supply. Regardless of its potential, a systems understanding of its tight regulatory network and the effects of the significant genomic changes in M1152 is missing. In this study, we compare M1152 to its ancestor M145, thereby connecting observed phenotypic differences to changes on transcription and translation. Measured protein levels are connected to predicted metabolic fluxes, facilitated by an enzyme-constrained genome-scale model (GEM), that by itself is a consensus result of a community effort. This approach connects observed differences in growth rate and glucose consumption to changes in central carbon metabolism, accompanied by differential expression of important regulons. Results suggest that precursors supply is not limiting secondary metabolism, informing that alternative strategies will be beneficial for further development of S. coelicolor for heterologous production of novel compounds.

Genome-scale metabolic modeling is an important tool in understanding metabolism, by enhancing collation of knowledge, interpretation of data, and prediction of metabolic capabilities. A central assumption in the construction and use of genome-scale models is that the in vivo organism is evolved for optimal growth, where growth is represented by flux through a biomass objective function (BOF). While the specific composition of the BOF is crucial, its formulation is often inherited from similar organisms due to the experimental challenges associated with its proper determination.However, a cell’s macro-molecular composition is not fixed and it responds to changes in environmental conditions. As a consequence, initiatives for the high-fidelity determination of cellular biomass composition have been launched. Thus, there is a need for a mathematical and computational framework capable of using multiple measurements of cellular biomass composition in different environments. Here, we propose two different computational approaches for directly addressing this challenge: Biomass Trade-off Weighting (BTW) and Higher-dimensional-plane InterPolation (HIP).In lieu of experimental data on biomass composition-variation in response to changing nutrient environment, we assess the properties of BTW and HIP using three hypothetical, yet biologically plausible, BOFs for the Escherichia coli genome-scale metabolic model iML1515. We find that the BTW and HIP formulations have a significant impact on model performance and phenotypes. Furthermore, the BTW method generates larger growth rates in all environments when compared to HIP. Using acetate secretion and the respiratory quotient as proxies for phenotypic changes, we find marked differences between the methods as HIP generates BOFs more similar to a reference BOF than BTW. We conclude that the presented methods constitute a first conceptual step in developing genome-scale metabolic modelling approaches capable of addressing the inherent dependence of cellular biomass composition on nutrient environments.Author summaryChanges in the environment promote changes in an organism’s metabolism. To achieve balanced growth states for near-optimal function, cells respond through metabolic rearrangements, which may influence the biosynthesis of metabolic precursors for building a cell’s molecular constituents. Therefore, it is necessary to take the dependence of biomass composition on environmental conditions into consideration. While measuring the biomass composition for some environments is possible, and should be done, it cannot be completed for all possible environments.In this work, we propose two main approaches, BTW and HIP, for addressing the challenge of estimating biomass composition in response to environmental changes. We evaluate the phenotypic consequences of BTW and HIP by characterizing their effect on growth, secretion potential, respiratory efficiency, and gene essentiality of a cell.Our work constitutes a first conceptual step in accounting for the influence of growth conditions on biomass composition, and in turn the biomass composition’s effect on metabolic phenotypic traits, within constraint-based modelling. As such, we believe it will improve the relevance of constraint-based methods in metabolic engineering and drug discovery, since the biosynthetic potential of microbes for generating industrially relevant products or drugs often is closely linked to their biomass composition.