Context-Specific Genome-Scale Metabolic Modelling and Its Application to the Analysis of COVID-19 Metabolic Signatures

Epigenetic control of metabolic identity across cell types

Constraint-based network modelling is a powerful tool for analysing cellular metabolism at genomic scale. Here, we conducted an integrative analysis of metabolic networks reconstructed from RNA-seq data with paired epigenomic data from the EpiATLAS resource of the International Human Epigenome Consortium (IHEC). Applying a state-of-the-art contextualisation algorithm, we reconstructed metabolic networks across 1,555 samples corresponding to 58 tissues and cell types. Analysis of these networks revealed the distribution of metabolic functionalities across human cell types and provides a compendium of human metabolic activity. This integrative approach allowed us to define, across tissues and cell types, i) reactions that fulfil the basic metabolic processes (core metabolism), and ii) cell type-specific functions (unique metabolism), that shape the metabolic identity of a cell or a tissue. Integration with EpiATLAS-derived cell type-specific gene-level chromatin states and enhancer-gene interactions identified enhancers, transcription factors, and key nodes controlling core and unique metabolism. Transport and first reactions of pathways were enriched for high expression, active chromatin state, and Polycomb-mediated repression in cell types where pathways are inactive, suggesting that key nodes are targets of repression. This integrative analysis forms the basis for identifying regulation points that control metabolic identity in human cells.

Identification of Anticancer Enzymes and Biomarkers for Hepatocellular Carcinoma through Constraint-Based Modeling

Hepatocellular carcinoma (HCC) results in the abnormal regulation of cellular metabolic pathways. Constraint-based modeling approaches can be utilized to dissect metabolic reprogramming, enabling the identification of biomarkers and anticancer targets for diagnosis and treatment. In this study, two genome-scale metabolic models (GSMMs) were reconstructed by employing RNA sequencing expression patterns of hepatocellular carcinoma (HCC) and their healthy counterparts. An anticancer target discovery (ACTD) framework was integrated with the two models to identify HCC targets for anticancer treatment. The ACTD framework encompassed four fuzzy objectives to assess both the suppression of cancer cell growth and the minimization of side effects during treatment. The composition of a nutrient may significantly affect target identification. Within the ACTD framework, ten distinct nutrient media were utilized to assess nutrient uptake for identifying potential anticancer enzymes. The findings revealed the successful identification of target enzymes within the cholesterol biosynthetic pathway using a cholesterol-free cell culture medium. Conversely, target enzymes in the cholesterol biosynthetic pathway were not identified when the nutrient uptake included a cholesterol component. Moreover, the enzymes PGS1 and CRL1 were detected in all ten nutrient media. Additionally, the ACTD framework comprises dual-group representations of target combinations, pairing a single-target enzyme with an additional nutrient uptake reaction. Additionally, the enzymes PGS1 and CRL1 were identified across the ten-nutrient media. Furthermore, the ACTD framework encompasses two-group representations of target combinations involving the pairing of a single-target enzyme with an additional nutrient uptake reaction. Computational analysis unveiled that cell viability for all dual-target combinations exceeded that of their respective single-target enzymes. Consequently, integrating a target enzyme while adjusting an additional exchange reaction could efficiently mitigate cell proliferation rates and ATP production in the treated cancer cells. Nevertheless, most dual-target combinations led to lower side effects in contrast to their single-target counterparts. Additionally, differential expression of metabolites between cancer cells and their healthy counterparts were assessed via parsimonious flux variability analysis employing the GSMMs to pinpoint potential biomarkers. The variabilities of the fluxes and metabolite flow rates in cancer and healthy cells were classified into seven categories. Accordingly, two secretions and thirteen uptakes (including eight essential amino acids and two conditionally essential amino acids) were identified as potential biomarkers. The findings of this study indicated that cancer cells exhibit a higher uptake of amino acids compared with their healthy counterparts.

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.

Metabolic models predict fotemustine and the combination of eflornithine/rifamycin and adapalene/cannabidiol for the treatment of gliomas

Gliomas are the most common type of malignant brain tumors, with glioblastoma multiforme (GBM) having a median survival of 15 months due to drug resistance and relapse. The treatment of gliomas relies on surgery, radiotherapy and chemotherapy. Only 12 anti-brain tumor chemotherapies (AntiBCs), mostly alkylating agents, have been approved so far. Glioma subtype-specific metabolic models were reconstructed to simulate metabolite exchanges, in silico knockouts and the prediction of drug and drug combinations for all three subtypes. The simulations were confronted with literature, high-throughput screenings (HTSs), xenograft and clinical trial data to validate the workflow and further prioritize the drug candidates. The three subtype models accurately displayed different degrees of dependencies toward glutamine and glutamate. Furthermore, 33 single drugs, mainly antimetabolites and TXNRD1-inhibitors, as well as 17 drug combinations were predicted as potential candidates for gliomas. Half of these drug candidates have been previously tested in HTSs. Half of the tested drug candidates reduce proliferation in cell lines and two-thirds in xenografts. Most combinations were predicted to be efficient for all three glioma types. However, eflornithine/rifamycin and cannabidiol/adapalene were predicted specifically for GBM and low-grade glioma, respectively. Most drug candidates had comparable efficiency in preclinical tests, cerebrospinal fluid bioavailability and mode-of-action to AntiBCs. However, fotemustine and valganciclovir alone and eflornithine and celecoxib in combination with AntiBCs improved the survival compared to AntiBCs in two-arms, phase I/II and higher glioma clinical trials. Our work highlights the potential of metabolic modeling in advancing glioma drug discovery, which accurately predicted metabolic vulnerabilities, repurposable drugs and combinations for the glioma subtypes.

State-of the-Art Constraint-Based Modeling of Microbial Metabolism: From Basics to Context-Specific Models with a Focus on Methanotrophs

Methanotrophy is the ability of an organism to capture and utilize the greenhouse gas, methane, as a source of energy-rich carbon. Over the years, significant progress has been made in understanding of mechanisms for methane utilization, mostly in bacterial systems, including the key metabolic pathways, regulation and the impact of various factors (iron, copper, calcium, lanthanum, and tungsten) on cell growth and methane bioconversion. The implementation of -omics approaches provided vast amount of heterogeneous data that require the adaptation or development of computational tools for a system-wide interrogative analysis of methanotrophy. The genome-scale mathematical modeling of its metabolism has been envisioned as one of the most productive strategies for the integration of muti-scale data to better understand methane metabolism and enable its biotechnological implementation. Herein, we provide an overview of various computational strategies implemented for methanotrophic systems. We highlight functional capabilities as well as limitations of the most popular web resources for the reconstruction, modification and optimization of the genome-scale metabolic models for methane-utilizing bacteria.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

Integrative computational modeling to unravel novel potential biomarkers in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a major health problem around the world. The management of this disease is complicated by the lack of noninvasive diagnostic tools and the few treatment options available. Better clinical outcomes can be achieved if HCC is detected early, but unfortunately, clinical signs appear when the disease is in its late stages. We aim to identify novel genes that can be targeted for the diagnosis and therapy of HCC. We performed a meta-analysis of transcriptomics data to identify differentially expressed genes and applied network analysis to identify hub genes. Fatty acid metabolism, complement and coagulation cascade, chemical carcinogenesis and retinol metabolism were identified as key pathways in HCC. Furthermore, we integrated transcriptomics data into a reference human genome-scale metabolic model to identify key reactions and subsystems relevant in HCC. We conclude that fatty acid activation, purine metabolism, vitamin D, and E metabolism are key processes in the development of HCC and therefore need to be further explored for the development of new therapies. We provide the first evidence that GABRP, HBG1 and DAK (TKFC) genes are important in HCC in humans and warrant further studies.