Computational Approaches on Stoichiometric and Kinetic Modeling for Efficient Strain Design

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multilocus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Self-controlled in silico gene knockdown strategies to enhance the sustainable production of heterologous terpenoid by Saccharomyces cerevisiae

Microbial bioengineering is a growing field for producing plant natural products (PNPs) in recent decades, using heterologous metabolic pathways in host cells. Once heterologous metabolic pathways have been introduced into host cells, traditional metabolic engineering techniques are employed to enhance the productivity and yield of PNP biosynthetic routes, as well as to manage competing pathways. The advent of computational biology has marked the beginning of a novel epoch in strain design through in silico methods. These methods utilize genome-scale metabolic models (GEMs) and flux optimization algorithms to facilitate rational design across the entire cellular metabolic network. However, the implementation of in silico strategies can often result in an uneven distribution of metabolic fluxes due to the rigid knocking out of endogenous genes, which can impede cell growth and ultimately impact the accumulation of target products. In this study, we creatively utilized synthetic biology to refine in silico strain design for efficient PNPs production. OptKnock simulation was performed on the GEM of Saccharomyces cerevisiae OA07, an engineered strain for oleanolic acid (OA) bioproduction that has been reported previously. The simulation predicted that the single deletion of fol1, fol2, fol3, abz1, and abz2, or a combined knockout of hfd1, ald2 and ald3 could improve its OA production. Consequently, strains EK1∼EK7 were constructed and cultivated. EK3 (OA07△fol3), EK5 (OA07△abz1), and EK6 (OA07△abz2) had significantly higher OA titers in a batch cultivation compared to the original strain OA07. However, these increases were less pronounced in the fed-batch mode, indicating that gene deletion did not support sustainable OA production. To address this, we designed a negative feedback circuit regulated by malonyl-CoA, a growth-associated intermediate whose synthesis served as a bypass to OA synthesis, at fol3, abz1, abz2, and at acetyl-CoA carboxylase-encoding gene acc1, to dynamically and autonomously regulate the expression of these genes in OA07. The constructed strains R_3A, R_5A and R_6A had significantly higher OA titers than the initial strain and the responding gene-knockout mutants in either batch or fed-batch culture modes. Among them, strain R_3A stand out with the highest OA titer reported to date. Its OA titer doubled that of the initial strain in the flask-level fed-batch cultivation, and achieved at 1.23 ± 0.04 g L-1 in 96 h in the fermenter-level fed-batch mode. This indicated that the integration of optimization algorithm and synthetic biology approaches was efficiently rational for PNP-producing strain design.

Computational analysis into the potential of azo dyes as a feedstock for actinorhodin biosynthesis in Pseudomonas putida

Fermentation-based biosynthesis in synthetic biology relies heavily on sugar-derived feedstocks, a limited and carbon-intensive commodity. Unconventional feedstocks from less-noble sources such as waste are being utilized to produce high-value chemical products. Azo dyes, a major pollutant commonly discharged by food, textile, and pharmaceutical industries, present significant health and environmental risks. We explore the potential of engineering Pseudomonas putida KT2440 to utilize azo dyes as a substrate to produce a polyketide, actinorhodin (ACT). Using the constrained minimal cut sets (cMCS) approach, we identified metabolic interventions that optimize ACT biosynthesis and compare the growth-coupling solutions attainable on an azo dye compared to glucose. Our results predicted that azo dyes could perform better as a feedstock for ACT biosynthesis than glucose as it allowed growth-coupling regimes that are unfeasible with glucose and generated an 18.28% higher maximum ACT flux. By examining the flux distributions enabled in different carbon sources, we observed that carbon fluxes from aromatic compounds like azo dyes have a unique capability to leverage gluconeogenesis to support both growth and production of secondary metabolites that produce excess NADH. Carbon sources are commonly chosen based on the host organism, availability, cost, and environmental implications. We demonstrated that careful selection of carbon sources is also crucial to ensure that the resulting flux distribution is suitable for further metabolic engineering of microbial cell factories.

Model-driven characterization of functional diversity of Pseudomonas aeruginosa clinical isolates with broadly representative phenotypes

Pseudomonas aeruginosa is a leading cause of infections in immunocompromised individuals and in healthcare settings. This study aims to understand the relationships between phenotypic diversity and the functional metabolic landscape of P. aeruginosa clinical isolates. To better understand the metabolic repertoire of P. aeruginosa in infection, we deeply profiled a representative set from a library of 971 clinical P. aeruginosa isolates with corresponding patient metadata and bacterial phenotypes. The genotypic clustering based on whole-genome sequencing of the isolates, multi-locus sequence types, and the phenotypic clustering generated from a multi-parametric analysis were compared to each other to assess the genotype-phenotype correlation. Genome-scale metabolic network reconstructions were developed for each isolate through amendments to an existing PA14 network reconstruction. These network reconstructions show diverse metabolic functionalities and enhance the collective P. aeruginosa pangenome metabolic repertoire. Characterizing this rich set of clinical P. aeruginosa isolates allows for a deeper understanding of the genotypic and metabolic diversity of the pathogen in a clinical setting and lays a foundation for further investigation of the metabolic landscape of this pathogen and host-associated metabolic differences during infection.

Insights on the Advancements of In Silico Metabolic Studies of Succinic Acid Producing Microorganisms: A Review with Emphasis on Actinobacillus succinogenes

Succinic acid (SA) is one of the top candidate value-added chemicals that can be produced from biomass via microbial fermentation. A considerable number of cell factories have been proposed in the past two decades as native as well as non-native SA producers. Actinobacillus succinogenes is among the best and earliest known natural SA producers. However, its industrial application has not yet been realized due to various underlying challenges. Previous studies revealed that the optimization of environmental conditions alone could not entirely resolve these critical problems. On the other hand, microbial in silico metabolic modeling approaches have lately been the center of attention and have been applied for the efficient production of valuable commodities including SA. Then again, literature survey results indicated the absence of up-to-date reviews assessing this issue, specifically concerning SA production. Hence, this review was designed to discuss accomplishments and future perspectives of in silico studies on the metabolic capabilities of SA producers. Herein, research progress on SA and A. succinogenes, pathways involved in SA production, metabolic models of SA-producing microorganisms, and status, limitations and prospects on in silico studies of A. succinogenes were elaborated. All in all, this review is believed to provide insights to understand the current scenario and to develop efficient mathematical models for designing robust SA-producing microbial strains.

Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.

Investigation of microbial community interactions between Lake Washington methanotrophs using -------genome-scale metabolic modeling

Background

The role of methane in global warming has become paramount to the environment and the human society, especially in the past few decades. Methane cycling microbial communities play an important role in the global methane cycle, which is why the characterization of these communities is critical to understand and manipulate their behavior. Methanotrophs are a major player in these communities and are able to oxidize methane as their primary carbon source.

Results

Lake Washington is a freshwater lake characterized by a methane-oxygen countergradient that contains a methane cycling microbial community. Methanotrophs are a major part of this community involved in assimilating methane from lake water. Two significant methanotrophic species in this community are Methylobacter and Methylomonas. In this work, these methanotrophs are computationally studied via developing highly curated genome-scale metabolic models. Each model was then integrated to form a community model with a multi-level optimization framework. The competitive and mutualistic metabolic interactions among Methylobacter and Methylomonas were also characterized. The community model was next tested under carbon, oxygen, and nitrogen limited conditions in addition to a nutrient-rich condition to observe the systematic shifts in the internal metabolic pathways and extracellular metabolite exchanges. Each condition showed variations in the methane oxidation pathway, pyruvate metabolism, and the TCA cycle as well as the excretion of formaldehyde and carbon di-oxide in the community. Finally, the community model was simulated under fixed ratios of these two members to reflect the opposing behavior in the two-member synthetic community and in sediment-incubated communities. The community simulations predicted a noticeable switch in intracellular carbon metabolism and formaldehyde transfer between community members in sediment-incubated vs. synthetic condition.

Conclusion

In this work, we attempted to predict the response of a simplified methane cycling microbial community from Lake Washington to varying environments and also provide an insight into the difference of dynamics in sediment-incubated microcosm community and synthetic co-cultures. Overall, this study lays the ground for in silico systems-level studies of freshwater lake ecosystems, which can drive future efforts of understanding, engineering, and modifying these communities for dealing with global warming issues.

Metabolic Modeling Elucidates the Transactions in the Rumen Microbiome and the Shifts Upon Virome Interactions

The complex microbial ecosystem within the bovine rumen plays a crucial role in host nutrition, health, and environmental impact. However, little is known about the interactions between the functional entities within the system, which dictates the community structure and functional dynamics and host physiology. With the advancements in high-throughput sequencing and mathematical modeling, in silico genome-scale metabolic analysis promises to expand our understanding of the metabolic interplay in the community. In an attempt to understand the interactions between microbial species and the phages inside rumen, a genome-scale metabolic modeling approach was utilized by using key members in the rumen microbiome (a bacteroidete, a firmicute, and an archaeon) and the viral phages associated with them. Individual microbial host models were integrated into a community model using multi-level mathematical frameworks. An elaborate and heuristics-based computational procedure was employed to predict previously unknown interactions involving the transfer of fatty acids, vitamins, coenzymes, amino acids, and sugars among the community members. While some of these interactions could be inferred by the available multi-omic datasets, our proposed method provides a systemic understanding of why the interactions occur and how these affect the dynamics in a complex microbial ecosystem. To elucidate the functional role of the virome on the microbiome, local alignment search was used to identify the metabolic functions of the viruses associated with the hosts. The incorporation of these functions demonstrated the role of viral auxiliary metabolic genes in relaxing the metabolic bottlenecks in the microbial hosts and complementing the inter-species interactions. Finally, a comparative statistical analysis of different biologically significant community fitness criteria identified the variation in flux space and robustness of metabolic capacities of the community members. Our elucidation of metabolite exchange among the three members of the rumen microbiome shows how their genomic differences and interactions with the viral strains shape up a highly sophisticated metabolic interplay and explains how such interactions across kingdoms can cause metabolic and compositional shifts in the community and affect the health, nutrition, and pathophysiology of the ruminant animal.

Modeling the Interplay between Photosynthesis, CO2 Fixation, and the Quinone Pool in a Purple Non-Sulfur Bacterium

Rhodopseudomonas palustris CGA009 is a purple non-sulfur bacterium that can fix carbon dioxide (CO₂) and nitrogen or break down organic compounds for its carbon and nitrogen requirements. Light, inorganic, and organic compounds can all be used for its source of energy. Excess electrons produced during its metabolic processes can be exploited to produce hydrogen gas or biodegradable polyesters. A genome-scale metabolic model of the bacterium was reconstructed to study the interactions between photosynthesis, CO₂ fixation, and the redox state of the quinone pool. A comparison of model-predicted flux values with available Metabolic Flux Analysis (MFA) fluxes yielded predicted errors of 5–19% across four different growth substrates. The model predicted the presence of an unidentified sink responsible for the oxidation of excess quinols generated by the TCA cycle. Furthermore, light-dependent energy production was found to be highly dependent on the quinol oxidation rate. Finally, the extent of CO₂ fixation was predicted to be dependent on the amount of ATP generated through the electron transport chain, with excess ATP going toward the energy-demanding Calvin-Benson-Bassham (CBB) pathway. Based on this analysis, it is hypothesized that the quinone redox state acts as a feed-forward controller of the CBB pathway, signaling the amount of ATP available.