Evaluating proteome allocation of Saccharomyces cerevisiae phenotypes with resource balance analysis

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Engineering of Saccharomyces cerevisiae for enhanced metabolic robustness and L-lactic acid production from lignocellulosic biomass

Metabolic engineering for high productivity and increased robustness is needed to enable sustainable biomanufacturing of lactic acid from lignocellulosic biomass. Lactic acid is an important commodity chemical used for instance as a monomer for production of polylactic acid, a biodegradable polymer. Here, rational and model-based optimization was used to engineer a diploid, xylose fermenting Saccharomyces cerevisiae strain to produce L-lactic acid. The metabolic flux was steered towards lactic acid through the introduction of multiple lactate dehydrogenase encoding genes while deleting ERF2, GPD1, and CYB2. A production of 93 g/L of lactic acid with a yield of 0.84 g/g was achieved using xylose as the carbon source. To increase xylose utilization and reduce acetic acid synthesis, PHO13 and ALD6 were also deleted from the strain. Finally, CDC19 encoding a pyruvate kinase was overexpressed, resulting in a yield of 0.75 g lactic acid/g sugars consumed, when the substrate used was a synthetic lignocellulosic hydrolysate medium, containing hexoses, pentoses and inhibitors such as acetate and furfural. Notably, modeling also provided leads for understanding the influence of oxygen in lactic acid production. High lactic acid production from xylose, at oxygen-limitation could be explained by a reduced flux through the oxidative phosphorylation pathway. On the contrast, higher oxygen levels were beneficial for lactic acid production with the synthetic hydrolysate medium, likely as higher ATP concentrations are needed for tolerating the inhibitors therein. The work highlights the potential of S. cerevisiae for industrial production of lactic acid from lignocellulosic biomass.

KETCHUP: Parameterizing of large-scale kinetic models using multiple datasets with different reference states

Large-scale kinetic models provide the computational means to dynamically link metabolic reaction fluxes to metabolite concentrations and enzyme levels while also conforming to substrate level regulation. However, the development of broadly applicable frameworks for efficiently and robustly parameterizing models remains a challenge. Challenges arise due to both the heterogeneity, paucity, and difficulty in obtaining flux and/or concentration data but also due to the computational difficulties of the underlying parameter identification problem. Both the computational demands for parameterization, degeneracy of obtained parameter solutions and interpretability of results has so far limited widespread adoption of large-scale kinetic models despite their potential. Herein, we introduce the Kinetic Estimation Tool Capturing Heterogeneous Datasets Using Pyomo (KETCHUP), a flexible parameter estimation tool that leverages a primal-dual interior-point algorithm to solve a nonlinear programming (NLP) problem that identifies a set of parameters capable of recapitulating the (non)steady-state fluxes and concentrations in wild-type and perturbed metabolic networks. KETCHUP is benchmarked against previously parameterized large-scale kinetic models demonstrating an at least an order of magnitude faster convergence than the tool K-FIT while at the same time attaining better data fits. This versatile toolbox accepts different kinetic descriptions, metabolic fluxes, enzyme levels and metabolite concentrations, under either steady-state or instationary conditions to enable robust kinetic model construction and parameterization. KETCHUP supports the SBML format and can be accessed at https://github.com/maranasgroup/KETCHUP.

StressME: Unified computing framework of Escherichia coli metabolism, gene expression, and stress responses

Generalist microbes have adapted to a multitude of environmental stresses through their integrated stress response system. Individual stress responses have been quantified by E. coli metabolism and expression (ME) models under thermal, oxidative and acid stress, respectively. However, the systematic quantification of cross-stress & cross-talk among these stress responses remains lacking. Here, we present StressME: the unified stress response model of E. coli combining thermal (FoldME), oxidative (OxidizeME) and acid (AcidifyME) stress responses. StressME is the most up to date ME model for E. coli and it reproduces all published single-stress ME models. Additionally, it includes refined rate constants to improve prediction accuracy for wild-type and stress-evolved strains. StressME revealed certain optimal proteome allocation strategies associated with cross-stress and cross-talk responses. These stress-optimal proteomes were shaped by trade-offs between protective vs. metabolic enzymes; cytoplasmic vs. periplasmic chaperones; and expression of stress-specific proteins. As StressME is tuned to compute metabolic and gene expression responses under mild acid, oxidative, and thermal stresses, it is useful for engineering and health applications. The modular design of our open-source package also facilitates model expansion (e.g., to new stress mechanisms) by the computational biology community.