Genome-scale Metabolic Models Research Articles

GMCSpy: Efficient and accurate computation of genetic minimal cut sets in python.

The identification of minimal genetic interventions that modulate metabolic processes constitutes one of the most relevant applications of genome-scale metabolic models (GEMs). The concept of Minimal Cut Sets (MCSs) and its extension at the gene level, genetic Minimal Cut Sets (gMCSs), have attracted increasing interest in the field of Systems Biology to address this task. Different computational tools have been developed to calculate MCSs and gMCSs using both commercial and open-source software. Here, we present gMCSpy, an efficient Python package to calculate gMCSs in GEMs using both commercial and non-commercial optimization solvers. We show that gMCSpy substantially overperforms our previous computational tool GMCS, which exclusively relied on commercial software. Moreover, we compared gMCSpy with recently published competing algorithms in the literature, finding significant improvements in both accuracy and computation time. All these advances make gMCSpy an attractive tool for researchers in the field of Systems Biology for different applications in health and biotechnology. The Python package gMCSpy can be accessed at: https://github.com/PlanesLab/gMCSpy. Supplementary data are available at Bioinformatics online.

Construction of an enzyme-constrained metabolic network model for Myceliophthora thermophila using machine learning-based kcat data

BackgroundGenome-scale metabolic models (GEMs) serve as effective tools for understanding cellular phenotypes and predicting engineering targets in the development of industrial strain. Enzyme-constrained genome-scale metabolic models (ecGEMs) have emerged as a valuable advancement, providing more accurate predictions and unveiling new engineering targets compared to models lacking enzyme constraints. In 2022, a stoichiometric GEM, iDL1450, was reconstructed for the industrially significant fungus Myceliophthora thermophila. To enhance the GEM’s performance, an ecGEM was developed for M. thermophila in this study.ResultsInitially, the model iDL1450 underwent refinement and updates, resulting in a new version named iYW1475. These updates included adjustments to biomass components, correction of gene-protein-reaction (GPR) rules, and a consensus on metabolites. Subsequently, the first ecGEM for M. thermophila was constructed using machine learning-based kcat data predicted by TurNuP within the ECMpy framework. During the construction, three versions of ecGEMs were developed based on three distinct kcat collection methods, namely AutoPACMEN, DLKcat and TurNuP. After comparison, the ecGEM constructed using TurNuP-predicted kcat values performed better in several aspects and was selected as the definitive version of ecGEM for M. thermophila (ecMTM). Comparing ecMTM to iYW1475, the solution space was reduced and the growth simulation results more closely resembled realistic cellular phenotypes. Metabolic adjustment simulated by ecMTM revealed a trade-off between biomass yield and enzyme usage efficiency at varying glucose uptake rates. Notably, hierarchical utilization of five carbon sources derived from plant biomass hydrolysis was accurately captured and explained by ecMTM. Furthermore, based on enzyme cost considerations, ecMTM successfully predicted reported targets for metabolic engineering modification and introduced some new potential targets for chemicals produced in M. thermophila.ConclusionsIn this study, the incorporation of enzyme constraint to iYW1475 not only improved prediction accuracy but also broadened the model’s applicability. This research demonstrates the effectiveness of integrating of machine learning-based kcat data in the construction of ecGEMs especially in situations where there is limited measured enzyme kinetic parameters for a specific organism.

Revealing microbiota characteristics and predicting flavor-producing sub-communities in Nongxiangxing baijiu pit mud through metagenomic analysis and metabolic modeling

The microorganisms of the pit mud (PM) of Nongxiangxing baijiu (NXXB) have an important role in the synthesis of flavor substances, and they determine attributes and quality of baijiu. Herein, we utilize metagenomics and genome-scale metabolic models (GSMMs) to investigate the microbial composition, metabolic functions in PM microbiota, as well as to identify microorganisms and communities linked to flavor compounds. Metagenomic data revealed that the most prevalent assembly of bacteria and archaea was Proteiniphilum, Caproicibacterium, Petrimonas, Lactobacillus, Clostridium, Aminobacterium, Syntrophomonas, Methanobacterium, Methanoculleus, and Methanosarcina. The important enzymes ofPMwere in bothGH and GT familymetabolism. A total of 38 high-quality metagenome-assembled genomes (MAGs) were obtained, including those at the family level (n = 13), genus level (n = 17), and species level (n = 8). GSMMs of the 38 MAGs were then constructed. From the GSMMs, individual and community capabilities respectively were predicted to be able to produce 111 metabolites and 598 metabolites. Twenty-three predicted metabolites were consistent with the metabonomics detected flavors and served as targets. Twelve sub-community of were screened by cross-feeding of 38 GSMMs. Of them, Methanobacterium, Sphaerochaeta, Muricomes intestini, Methanobacteriaceae, Synergistaceae, and Caloramator were core microorganisms for targets in each sub-community. Overall, this study of metagenomic and target-community screening could help our understanding of the metabolite–microbiome association and further bioregulation of baijiu.

Integration of genome-scale metabolic model with biorefinery process model reveals market-competitive carbon-negative sustainable aviation fuel utilizing microbial cell mass lipids and biogenic CO2

Producing scalable, economically viable, low-carbon biofuels or biochemicals hinges on more efficient bioconversion processes. While microbial conversion can offer robust solutions, the native microbial growth process often redirects a large fraction of carbon to CO2 and cell mass. By integrating genome-scale metabolic models with techno-economic and life cycle assessment models, this study analyzes the effects of converting cell mass lipids to hydrocarbon fuels, and CO2 to methanol on the facility’s costs and life-cycle carbon footprint. Results show that upgrading microbial lipids or both microbial lipids and CO2 using renewable hydrogen produces carbon-negative bisabolene. Additionally, on-site electrolytic hydrogen production offers a supply of pure oxygen to use in place of air for bioconversion and fuel combustion in the boiler. To reach cost parity with conventional jet fuel, renewable hydrogen needs to be produced at less than $2.2 to $3.1/kg, with a bisabolene yield of 80% of the theoretical yield, along with cell mass and CO2 yields of 22 wt% and 54 wt%, respectively. The economic combination of cell mass, CO2, and bisabolene yields demonstrated in this study provides practical insights for prioritizing research, selecting suitable hosts, and determining necessary engineered production levels.

Genome-scale model of Rothia mucilaginosa predicts gene essentialities and reveals metabolic capabilities.

Cystic fibrosis (CF), an inherited genetic disorder caused by mutations in the cystic fibrosis transmembrane conductance regulator gene, results in sticky and thick mucosal fluids. This environment facilitates the colonization of various microorganisms, some of which can cause acute and chronic lung infections, while others may positively impact the disease. Rothia mucilaginosa, an oral commensal, is relatively abundant in the lungs of CF patients. Recent studies have unveiled its anti-inflammatory properties using in vitro three-dimensional lung epithelial cell cultures and in vivo mouse models relevant to chronic lung diseases. Apart from this, R. mucilaginosa has been associated with severe infections. However, its metabolic capabilities and genotype-phenotype relationships remain largely unknown. To gain insights into its cellular metabolism and genetic content, we developed the first manually curated genome-scale metabolic model, iRM23NL. Through growth kinetics and high-throughput phenotypic microarray testings, we defined its complete catabolic phenome. Subsequently, we assessed the model's effectiveness in accurately predicting growth behaviors and utilizing multiple substrates. We used constraint-based modeling techniques to formulate novel hypotheses that could expedite the development of antimicrobial strategies. More specifically, we detected putative essential genes and assessed their effect on metabolism under varying nutritional conditions. These predictions could offer novel potential antimicrobial targets without laborious large-scale screening of knockouts and mutant transposon libraries. Overall, iRM23NL demonstrates a solid capability to predict cellular phenotypes and holds immense potential as a valuable resource for accurate predictions in advancing antimicrobial therapies. Moreover, it can guide metabolic engineering to tailor R. mucilaginosa's metabolism for desired performance.IMPORTANCECystic fibrosis (CF) is a genetic disorder characterized by thick mucosal secretions, leading to chronic lung infections. Rothia mucilaginosa is a common bacterium found in various parts of the human body, acting as a normal part of the flora. In people with weakened immune systems, it can become an opportunistic pathogen, while it is prevalent and active in CF airways. Recent studies have highlighted its anti-inflammatory properties in the lower pulmonary system, indicating the intricate relationship between microbes and human health. Herein, we have developed the first manually curated metabolic model of R. mucilaginosa. Our study examined the previously unknown relationships between the bacterium's genotype and phenotype and identified essential genes that impact the metabolism under various conditions. With this, we opt for paving the way for developing new strategies in antimicrobial therapy and metabolic engineering, leading to enhanced therapeutic outcomes in cystic fibrosis and related conditions.

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.

Dissecting Holistic Metabolic Acclimatization of Mucor circinelloides WJ11 Defective in Carotenoid Biosynthesis.

Mucor circinelloides WJ11 is a lipid-producing strain with industrial potential. A holistic approach using gene manipulation and bioprocessing development has improved lipid production and the strain's economic viability. However, the systematic regulation of lipid accumulation and carotenoid biosynthesis in M. circinelloides remains unknown. To dissect the metabolic mechanism underlying lipid and carotenoid biosynthesis, transcriptome analysis and reporter metabolites identification were implemented between the wild-type (WJ11) and ΔcarRP WJ11 strains of M. circinelloides. As a result, transcriptome analysis revealed 10,287 expressed genes, with 657 differentially expressed genes (DEGs) primarily involved in amino acid, carbohydrate, and energy metabolism. Integration with a genome-scale metabolic model (GSMM) identified reporter metabolites in the ΔcarRP WJ11 strain, highlighting metabolic pathways crucial for amino acid, energy, and nitrogen metabolism. Notably, the downregulation of genes associated with carotenoid biosynthesis and acetyl-CoA generation suggests a coordinated relationship between the carotenoid and fatty acid biosynthesis pathways. Despite disruptions in the carotenoid pathway, lipid production remains stagnant due to reduced acetyl-CoA availability, emphasizing the intricate metabolic interplay. These findings provide insights into the coordinated relationship between carotenoid and fatty acid biosynthesis in M. circinelloides that are valuable in applied research to design optimized strains for producing desired bioproducts through emerging technology.

Decipher syntrophies within C2-C4 organic acids-degrading anaerobic microbiomes: A multi-omic exploration

In methanogenic ecosystems, carboxylic acid degradation is a crucial process facilitated by specialized bacteria working in syntrophy with methanogens. This study explores the transcriptomic responses in 178 metagenome-assembled genomes upon changes in feeding composition, specifically lactate, butyrate, propionate, and acetate, in four lab-scale bioreactors. Methanothrix soehngenii and Methanoculleus sp. emerged as key biomethanation contributors. Through metatranscriptomics, machine learning-based functional annotation, and flux balance analysis the underlying microbial interactions and dynamics were deciphered. Syntrophomonadaceae and Mesotoga species maintained mutualistic metabolite exchanges with hydrogenotrophic methanogens, degrading primarily butyrate and lactate, respectively. Acetate was mainly consumed by Smithellaceae sp., in competition with M. soehngenii in all reactors. Furthermore, the acetoclastic archaeon exhibited a previously undocumented capability to metabolize lactate, thereby confirming the prevalence of acetotrophic pathway. Transcriptomic profiles revealed an additional layer of complexity, where propionate dismutation and beta-oxidation pathways are interconnected by the exchanges of butyrate between putative syntrophs, including Syntrophomonadaceae and Smithellaceae species. Integrating multi-omics data with genome-scale metabolic modeling enabled the accurate reconstruction of dynamics within the controlled ecosystem. This composite novel approach was applied to the AD system to unveil the intricate relationships operating within the microbiome to promote thriving. Results elucidated the still poorly explored organization of the so-called microbial dark matter.

Bioinformatic Analysis of Metabolomic Data: From Raw Spectra to Biological Insight

Metabolites are at the end of the gene–transcript–protein–metabolism cascade. As such, metabolomics is the omics approach that offers the most direct correlation with phenotype. This allows, where genomics, transcriptomics and proteomics fail to explain a trait, metabolomics to possibly provide an answer. Complex phenotypes, which are determined by the influence of multiple small-effect alleles, are an example of these situations. Consequently, the interest in metabolomics has increased exponentially in recent years. As a newer discipline, metabolomic bioinformatic analysis pipelines are not as standardized as in the other omics approaches. In this review, we synthesized the different steps that need to be carried out to obtain biological insight from annotated metabolite abundance raw data. These steps were grouped into three different modules: preprocessing, statistical analysis, and metabolic pathway enrichment. We included within each one of them the different state-of-the-art procedures and tools that can be used depending on the characteristics of the study, providing details about each method’s characteristics and the issues the reader might encounter. Finally, we introduce genome-scale metabolic modeling as a tool for obtaining pseudo-metabolomic data in situations where their acquisition is difficult, enabling the analysis of the resulting data with the modules of the described workflow.

ECMpy 2.0: A Python package for automated construction and analysis of enzyme-constrained models

Genome-scale metabolic models (GEMs) have been widely employed to predict microorganism behaviors. However, GEMs only consider stoichiometric constraints, leading to a linear increase in simulated growth and product yields as substrate uptake rates rise. This divergence from experimental measurements prompted the creation of enzyme-constrained models (ecModels) for various species, successfully enhancing chemical production. Building upon studies that allocate macromolecule resources, we developed a Python-based workflow (ECMpy) that constructs an enzyme-constrained model. This involves directly imposing an enzyme amount constraint in GEM and accounting for protein subunit composition in reactions. However, this procedure demands manual collection of enzyme kinetic parameter information and subunit composition details, making it rather user-unfriendly. In this work, we've enhanced the ECMpy toolbox to version 2.0, broadening its scope to automatically generate ecGEMs for a wider array of organisms. ECMpy 2.0 automates the retrieval of enzyme kinetic parameters and employs machine learning for predicting these parameters, which significantly enhances parameter coverage. Additionally, ECMpy 2.0 introduces common analytical and visualization features for ecModels, rendering computational results more user accessible. Furthermore, ECMpy 2.0 seamlessly integrates three published algorithms that exploit ecModels to uncover potential targets for metabolic engineering. ECMpy 2.0 is available at https://github.com/tibbdc/ECMpy or as a pip package (https://pypi.org/project/ECMpy/).

Genome-scale model development and genomic sequencing of the oleaginous clade Lipomyces.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

Logic programming-based Minimal Cut Sets reveal consortium-level therapeutic targets for chronic wound infections

Minimal Cut Sets (MCSs) identify sets of reactions which, when removed from a metabolic network, disable certain cellular functions. The traditional search for MCSs within genome-scale metabolic models (GSMMs) targets cellular growth, identifies reaction sets resulting in a lethal phenotype if disrupted, and retrieves a list of corresponding gene, mRNA, or enzyme targets. Using the dual link between MCSs and Elementary Flux Modes (EFMs), our logic programming-based tool aspefm was able to compute MCSs of any size from GSMMs in acceptable run times. The tool demonstrated better performance when computing large-sized MCSs than the mixed-integer linear programming methods. We applied the new MCSs methodology to a medically-relevant consortium model of two cross-feeding bacteria, Staphylococcus aureus and Pseudomonas aeruginosa. aspefm constraints were used to bias the computation of MCSs toward exchanged metabolites that could complement lethal phenotypes in individual species. We found that interspecies metabolite exchanges could play an essential role in rescuing single-species growth, for instance inosine could complement lethal reaction knock-outs in the purine synthesis, glycolysis, and pentose phosphate pathways of both bacteria. Finally, MCSs were used to derive a list of promising enzyme targets for consortium-level therapeutic applications that cannot be circumvented via interspecies metabolite exchange.

Review and meta-analysis of the genetic Minimal Cut Set approach for gene essentiality prediction in cancer metabolism.

Cancer metabolism is a marvellously complex topic, in part, due to the reprogramming of its pathways to self-sustain the malignant phenotype in the disease, to the detriment of its healthy counterpart. Understanding these adjustments can provide novel targeted therapies that could disrupt and impair proliferation of cancerous cells. For this very purpose, genome-scale metabolic models (GEMs) have been developed, with Human1 being the most recent reconstruction of the human metabolism. Based on GEMs, we introduced the genetic Minimal Cut Set (gMCS) approach, an uncontextualized methodology that exploits the concepts of synthetic lethality to predict metabolic vulnerabilities in cancer. gMCSs define a set of genes whose knockout would render the cell unviable by disrupting an essential metabolic task in GEMs, thus, making cellular proliferation impossible. Here, we summarize the gMCS approach and review the current state of the methodology by performing a systematic meta-analysis based on two datasets of gene essentiality in cancer. First, we assess several thresholds and distinct methodologies for discerning highly and lowly expressed genes. Then, we address the premise that gMCSs of distinct length should have the same predictive power. Finally, we question the importance of a gene partaking in multiple gMCSs and analyze the importance of all the essential metabolic tasks defined in Human1. Our meta-analysis resulted in parameter evaluation to increase the predictive power for the gMCS approach, as well as a significant reduction of computation times by only selecting the crucial gMCS lengths, proposing the pertinency of particular parameters for the peak processing of gMCS.

Investigating the Role of β-Disodium Glycerophosphate and Urea in Promoting Growth of Streptococcus thermophilus from Omics-Integrated Genome-Scale Models.

This study investigates the impact of urea and β-GP on the growth of Streptococcus thermophilus S-3, a bacterium commonly used in industrial fermentation processes. Through a series of growth experiments, transcriptome, metabolome, and omics-based analyses, the research demonstrates that both urea and β-GP can enhance the biomass of S. thermophilus, with urea showing a more significant effect. The optimal urea concentration for growth was determined to be 3 g/L in M17 medium. The study also highlights the metabolic pathways influenced by urea and β-GP, particularly the galactose metabolism pathway, which is crucial for cell growth when lactose is the substrate. The integration of omics data into the genome-scale metabolic model of S. thermophilus, iCH502, allowed for a more accurate prediction of metabolic fluxes and growth rates. The study concludes that urea can serve as a viable substitute for β-GP in the cultivation of S. thermophilus, offering potential cost and efficiency benefits in industrial fermentation processes. The findings are supported by validation experiments with 11 additional strains of S. thermophilus, which showed increased biomass in UM17 medium.

Reconstruction of the genome-scale metabolic network model of Sinorhizobium fredii CCBAU45436 for free-living and symbiotic states.

Sinorhizobium fredii CCBAU45436 is an excellent rhizobium that plays an important role in agricultural production. However, there still needs more comprehensive understanding of the metabolic system of S. fredii CCBAU45436, which hinders its application in agriculture. Therefore, based on the first-generation metabolic model iCC541 we developed a new genome-scale metabolic model iAQY970, which contains 970 genes, 1,052 reactions, 942 metabolites and is scored 89% in the MEMOTE test. Cell growth phenotype predicted by iAQY970 is 81.7% consistent with the experimental data. The results of mapping the proteome data under free-living and symbiosis conditions to the model showed that the biomass production rate in the logarithmic phase was faster than that in the stable phase, and the nitrogen fixation efficiency of rhizobia parasitized in cultivated soybean was higher than that in wild-type soybean, which was consistent with the actual situation. In the symbiotic condition, there are 184 genes that would affect growth, of which 94 are essential; In the free-living condition, there are 143 genes that influence growth, of which 78 are essential. Among them, 86 of the 94 essential genes in the symbiotic condition were consistent with the prediction of iCC541, and 44 essential genes were confirmed by literature information; meanwhile, 30 genes were identified by DEG and 33 genes were identified by Geptop. In addition, we extracted four key nitrogen fixation modules from the model and predicted that sulfite reductase (EC 1.8.7.1) and nitrogenase (EC 1.18.6.1) as the target enzymes to enhance nitrogen fixation by MOMA, which provided a potential focus for strain optimization. Through the comprehensive metabolic model, we can better understand the metabolic capabilities of S. fredii CCBAU45436 and make full use of it in the future.

Abstract 4673: Genome-scale metabolic models identify gene targets that enhance drug response of drug-resistant breast cancer cells

Abstract Adaptive drug resistance poses a significant therapeutic challenge in treatment of breast cancer patients. Metabolic reprogramming in cancer cells is one of the major mechanisms associated with drug resistance. In this study, we aim to identify gene targets that enhance the drug response of drug-resistant breast cancer cells through downregulation. We hypothesized that metabolic perturbation in drug-resistant breast cancer cells could revert their metabolic features to those of drug-sensitive cells. To predict such gene targets, we used a computational model known as a genome-scale metabolic model (GEM). A GEM simulates the entire metabolic reaction fluxes of a cell under varying genetic and/or environmental conditions. We specifically targeted two engineered MCF7 breast cancer cell lines: one resistant to doxorubicin and the other to paclitaxel. First, we generated proteome data from drug-sensitive and drug-resistant MCF7 cell lines, and used them to construct the cell-specific GEMs. Subsequently, single-gene knockout simulation was conducted using these cell-specific GEMs. The initial list of the gene knockout targets was subjected to pathway analysis involving proteome and metabolome data to finalize the gene targets. The final gene targets were expected to lead to metabolic features similar to those of drug-sensitive MCF7 upon the downregulation. In vitro experiments using inhibitors were conducted to validate these predicted gene targets. This study provides a computational approach that can also be applied to other types of drug-resistant cancer cells. Citation Format: Hae Deok Jung, JinA Lim, Han Suk Ryu, Yoosik Kim, Hyun Uk Kim. Genome-scale metabolic models identify gene targets that enhance drug response of drug-resistant breast cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4673.

Uncovering the Effect of SARS-CoV-2 on Liver Metabolism via Genome-Scale Metabolic Modeling for Reprogramming and Therapeutic Strategies.

Genome-scale metabolic models (GEMs) are promising computational tools that contribute to elucidating host-virus interactions at the system level and developing therapeutic strategies against viral infection. In this study, the effect of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) on liver metabolism was investigated using integrated GEMs of human hepatocytes and SARS-CoV-2. They were generated for uninfected and infected hepatocytes using transcriptome data. Reporter metabolite analysis resulted in significant transcriptional changes around several metabolites involved in xenobiotics, drugs, arachidonic acid, and leukotriene metabolisms due to SARS-CoV-2 infection. Flux balance analysis and minimization of metabolic adjustment approaches unraveled possible virus-induced hepatocellular reprogramming in fatty acid, glycerophospholipid, sphingolipid cholesterol, and folate metabolisms, bile acid biosynthesis, and carnitine shuttle among others. Reaction knockout analysis provided critical reactions in glycolysis, oxidative phosphorylation, purine metabolism, and reactive oxygen species detoxification subsystems. Computational analysis also showed that administration of dopamine, glucosamine, D-xylose, cysteine, and (R)-3-hydroxybutanoate contributes to alleviating viral infection. In essence, the reconstructed host-virus GEM helps us understand metabolic programming and develop therapeutic strategies to battle SARS-CoV-2.

Abstract 7036: Integrative biomarker profiling for predicting immunotherapy outcomes in gynecologic cancers

Abstract Immunotherapy with immune checkpoint inhibitors (ICI) has proven effective for various solid tumors, such as those in the lungs, urinary tract, and skin. However, its efficacy in gynecologic (GYN) cancers has been limited, with response rates of 11-17% in recurrent cases. Only a minority of patients with DNA mismatch repair deficiencies exhibit favorable responses to PD-1 inhibitors. The predictability of PD-L1 expression in GYN cancers remains inconsistent. Our research aimed to identify and prioritize biomarkers and mechanisms linked to ICI resistance in GYN oncology, targeting those who could benefit from these treatments.We gathered clinical and pathological data from 49 patients with ovarian, endometrial, and cervical cancers at various stages, all treated with ICI. We performed comprehensive analyses, including transcriptome-based immune cell abundance estimation, signaling pathway enrichment, and correlation with ICI treatment responses and survival outcomes. Our previous work highlighted TGF-β and CD47's roles in immunotherapy resistance, facilitating an immunosuppressive environment within our GYN cancer cohort. Moreover, we conducted a stoichiometry constraint-based genome-scale metabolic model and computed correlation-grouped metabolic activity scores (cgScore) using metabolic enzyme gene expression as constraints. Among 1533 metabolic reactions, pyruvate transport and folate metabolism stood out as the most upregulated in responders versus those who progressed.Upon identifying several promising biomarkers, we used logistic regression classifiers to decipher the predictive value of each. Our model included clinical-demographic (age, BMI, cancer type, stage) and molecular-genomic features (immune marker gene expression, TopTGF-β score, immune cell abundance from transcriptomic deconvolution, and the top 2 cgScores). Analysis divided into training and testing sets highlighted features such as stage, BMI, eosinophils, CD274, CD47, and the pyruvate transport reaction as significant predictors of ICI response. The model achieved an 80% accuracy rate on the test set, confirmed by permutation importance from scikit-learn. We also used a Regularized Cox model to link predictors with progression-free survival, emphasizing the roles of BMI, plasma cells, CD276, pyruvate transport reaction, and folate metabolism. In conclusion, there are immunological and metabolic markers have the potential to serve as predictive indicators for ICI treatment response and survival in GYN cancers. Understanding these predictors can profoundly influence clinical decision-making in immunotherapy, ultimately leading to enhanced patient outcomes and personalized treatment approaches. Citation Format: Gideon Idumah, Lin Li, Haider Mahdi, Ying Ni. Integrative biomarker profiling for predicting immunotherapy outcomes in gynecologic cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7036.

Design of a proteolytic module for improved metabolic modeling of Bacteroides caccae.

The gut microbiota plays a crucial role in health and is significantly modulated by human diets. In addition to Western diets which are rich in proteins, high-protein diets are used for specific populations or indications, mainly weight loss. In this study, we investigated the effect of protein supplementation on Bacteroides caccae, a Gram-negative gut symbiont. The supplementation with whey proteins led to a significant increase in growth rate, final biomass, and short-chain fatty acids production. A comprehensive genomic analysis revealed that B. caccae possesses a set of 156 proteases with putative intracellular and extracellular localization and allowed to identify amino acid transporters and metabolic pathways. We developed a fully curated genome-scale metabolic model of B. caccae that incorporated its proteolytic activity and simulated its growth and production of fermentation-related metabolites in response to the different growth media. We validated the model by comparing the predicted phenotype to experimental data. The model accurately predicted B. caccae's growth and metabolite production (R2 = 0.92 for the training set and R2 = 0.89 for the validation set). We found that accounting for both ATP consumption related to proteolysis, and whey protein accessibility is necessary for accurate predictions of metabolites production. These results provide insights into B. caccae's adaptation to a high-protein diet and its ability to utilize proteins as a source of nutrition. The proposed model provides a useful tool for understanding the feeding mechanism of B. caccae in the gut microbiome.IMPORTANCEMicrobial proteolysis is understudied despite the availability of dietary proteins for the gut microbiota. Here, the proteolytic potential of the gut symbiont Bacteroides caccae was analyzed for the first time using pan-genomics. This sketches a well-equipped bacteria for protein breakdown, capable of producing 156 different proteases with a broad spectrum of cleavage targets. This functional potential was confirmed by the enhancement of growth and metabolic activities at high protein levels. Proteolysis was included in a B. caccae metabolic model which was fitted with the experiments and validated on external data. This model pinpoints the links between protein availability and short-chain fatty acids production, and the importance for B. caccae to gain access to glutamate and asparagine to promote growth. This integrated approach can be generalized to other symbionts and upscaled to complex microbiota to get insights into the ecological impact of proteins on the gut microbiota.

Exposure to environmental toxicants is associated with gut microbiome dysbiosis, insulin resistance and obesity

Environmental toxicants (ETs) are associated with adverse health outcomes. Here we hypothesized that exposures to ETs are linked with obesity and insulin resistance partly through a dysbiotic gut microbiota and changes in the serum levels of secondary bile acids (BAs). Serum BAs, per- and polyfluoroalkyl substances (PFAS) and additional twenty-seven ETs were measured by mass spectrometry in 264 Danes (121 men and 143 women, aged 56.6 ± 7.3 years, BMI 29.7 ± 6.0 kg/m2) using a combination of targeted and suspect screening approaches. Bacterial species were identified based on whole-genome shotgun sequencing (WGS) of DNA extracted from stool samples. Personalized genome-scale metabolic models (GEMs) of gut microbial communities were developed to elucidate regulation of BA pathways. Subsequently, we compared findings from the human study with metabolic implications of exposure to perfluorooctanoic acid (PFOA) in PPARα-humanized mice. Serum levels of twelve ETs were associated with obesity and insulin resistance. High chemical exposure was associated with increased abundance of several bacterial species (spp.) of genus (Anaerotruncus, Alistipes, Bacteroides, Bifidobacterium, Clostridium, Dorea, Eubacterium, Escherichia, Prevotella, Ruminococcus, Roseburia, Subdoligranulum, and Veillonella), particularly in men. Conversely, females in the higher exposure group, showed a decrease abundance of Prevotella copri. High concentrations of ETs were correlated with increased levels of secondary BAs including lithocholic acid (LCA), and decreased levels of ursodeoxycholic acid (UDCA). In silico causal inference analyses suggested that microbiome-derived secondary BAs may act as mediators between ETs and obesity or insulin resistance. Furthermore, these findings were substantiated by the outcome of the murine exposure study. Our combined epidemiological and mechanistic studies suggest that multiple ETs may play a role in the etiology of obesity and insulin resistance. These effects may arise from disruptions in the microbial biosynthesis of secondary BAs.