Phenotype-based cell-specific metabolic modeling reveals metabolic liabilities of cancer

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.

Towards personalized genome-scale modeling of inborn errors of metabolism for systems medicine applications

Inborn errors of metabolism (IEMs) are a group of more than 1000 inherited diseases that are individually rare but have a cumulative global prevalence of 50 per 100,000 births. Recently, it has been recognized that like common diseases, patients with rare diseases can greatly vary in the manifestation and severity of symptoms. Here, we review omics-driven approaches that enable an integrated, holistic view of metabolic phenotypes in IEM patients. We focus on applications of Constraint-based Reconstruction and Analysis (COBRA), a widely used mechanistic systems biology approach, to model the effects of inherited diseases. Moreover, we review evidence that the gut microbiome is also altered in rare diseases. Finally, we outline an approach using personalized metabolic models of IEM patients for the prediction of biomarkers and tailored therapeutic or dietary interventions. Such applications could pave the way towards personalized medicine not just for common, but also for rare diseases.

Identification of mitochondrial-related signature and molecular subtype for the prognosis of osteosarcoma

Mitochondria play a vital role in osteosarcoma. Therefore, the purpose of this study was to investigate the potential role of mitochondrial-related genes (MRGs) in osteosarcoma. Based on 92 differentially expressed MRGs, osteosarcoma samples were divided into two subtypes using the nonnegative matrix factorization (NMF). Ultimately, a univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox analysis were performed to construct a prognostic risk model. The single-sample gene set enrichment analysis assessed the immune infiltration characteristics of osteosarcoma patients. Finally, we identified an osteosarcoma biomarker, malonyl-CoA decarboxylase (MLYCD), which showed downregulation. Osteosarcoma cells proliferation, migration, and invasion were effectively inhibited by the overexpression of MLYCD. Our findings will help us to further understand the molecular mechanisms of osteosarcoma and contribute to the discovery of new diagnostic biomarkers and therapeutic targets.

The Warburg Effect Explained: Integration of Enhanced Glycolysis with Heterogeneous Mitochondria to Promote Cancer Cell Proliferation

The Warburg effect is the long-standing riddle of cancer biology. How does aerobic glycolysis, inefficient in producing ATP, confer a growth advantage to cancer cells? A new evaluation of a large set of literature findings covering the Warburg effect and its yeast counterpart, the Crabtree effect, led to an innovative working hypothesis presented here. It holds that enhanced glycolysis partially inactivates oxidative phosphorylation to induce functional rewiring of a set of TCA cycle enzymes to generate new non-canonical metabolic pathways that sustain faster growth rates. The hypothesis has been structured by constructing two metabolic maps, one for cancer metabolism and the other for the yeast Crabtree effect. New lines of investigation, suggested by these maps, are discussed as instrumental in leading toward a better understanding of cancer biology in order to allow the development of more efficient metabolism-targeted anticancer drugs.

SL-scan identifies synthetic lethal interactions in cancer using metabolic networks

Exploiting synthetic lethality is a promising strategy for developing targeted cancer therapies. However, identifying clinically significant synthetic lethal (SL) interactions among a large number of gene combinations is a challenging computational task. In this study, we developed the SL-scan pipeline based on metabolic network modeling to discover SL interaction. The SL-scan pipeline identifies the association between simulated Flux Balance Analysis knockout scores and mutation data across cancer cell lines and predicts putative SL interactions. We assessed the concordance of the SL pairs predicted by SL-scan with those of obtained from analysis of the CRISPR, shRNA, and PRISM datasets. Our results demonstrate that the SL-scan pipeline outperformed existing SL prediction approaches based on metabolic networks in identifying SL pairs in various cancers. This study emphasizes the importance of integrating multiple data sources, particularly mutation data, when identifying SL pairs for targeted cancer therapies. The findings of this study may lead to the development of novel targeted cancer therapies.

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.

MCD Inhibits Lipid Deposition in Goat Intramuscular Preadipocytes

Malonyl-CoA decarboxylase (MCD) is a major regulator of fatty acid oxidation catalyzing the decarboxylation of malonyl coenzyme A (malonyl-CoA). Although its involvement in human diseases has been well studied, its role in intramuscular fat (IMF) deposition remains unknown. In this present study, 1726 bp of MCD cDNA was cloned (OM937122) from goat liver, including 5'UTR of 27 bp, 3'UTR of 199 bp, and CDS of 1500 bp, encoding 499 amino acids. In this present study, although the overexpression of MCD increased the mRNA expression of FASN and DGAT2, the expression of ATGL and ACOX1 was also activated significantly and resulted in a decrease in cellular lipid deposition in goat intramuscular preadipocytes. Meanwhile, the silencing of MCD increased the cellular lipid deposition and was accompanied by the expression activation of DGAT2 and the expression suppression of ATGL and HSL, despite the expression suppression of genes related to fatty acid synthesis, including ACC and FASN. However, the expression of DGAT1 was not affected significantly (p > 0.05) by the expression alteration of MCD in this present study. Furthermore, 2025 bp of MCD promoter was obtained and predicted to be regulated by C/EBPα, SP1, SREBP1, and PPARG. In summary, although different pathways may respond to the expression alteration of MCD, the expression of MCD was negatively correlated with the cellular lipid deposition in goat intramuscular preadipocytes. These data may be beneficial for enhancing our understanding of the regulation of IMF deposition in goats.

Lysine Malonylation and Its Links to Metabolism and Diseases

Malonylation is a recently identified post-translational modification with malonyl-coenzyme A as the donor. It conserved both in prokaryotes and eukaryotes. Recent advances in the identification and quantification of lysine malonylation by bioinformatic analysis have improved our understanding of its role in the regulation of protein activity, interaction, and localization and have elucidated its involvement in many biological processes. Malonylation has been linked to diverse physiological processes, including metabolic disorders, inflammation, and immune regulation. This review discusses malonylation in theory, describes the underlying mechanism, and summarizes the recent progress in malonylation research. The latest findings point to novel functions of malonylation and highlight the mechanisms by which malonylation regulates a variety of cellular processes. Our review also marks the association between lysine malonylation, the enzymes involved, and various diseases, and discusses promising diagnostic and therapeutic biomolecular targets for future clinical applications.

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

Genome-scale metabolic models (GEMs) have found numerous applications in different domains, ranging from biotechnology to systems medicine. Herein, we overview the most popular algorithms for the automated reconstruction of context-specific GEMs using high-throughput experimental data. Moreover, we describe different datasets applied in the process, and protocols that can be used to further automate the model reconstruction and validation. Finally, we describe recent COVID-19 applications of context-specific GEMs, focusing on the analysis of metabolic implications, identification of biomarkers and potential drug targets.

Cholangiocarcinoma progression depends on the uptake and metabolization of extracellular lipids

BACKGROUND AND AIMS

Cholangiocarcinoma (CCA) includes a heterogeneous group of biliary cancers with a dismal prognosis. We investigated if lipid metabolism is disrupted in CCA and its role in tumor proliferation.

APPROACH AND RESULTS

The in vitro and in vivo tumorigenic capacity of five human CCA cell lines was analyzed. Proteome, lipid content, and metabolic fluxes were evaluated in CCA cells and compared with normal human cholangiocytes (NHC). The Akt1/NOTCH1 intracellular cytoplasmic domain (Nicd1)-driven CCA mouse model was also evaluated. The proteome of CCA cells was enriched in pathways involved in lipid and lipoprotein metabolism. The EGI1 CCA cell line presented the highest tumorigenic capacity. Metabolic studies in high (EGI1) versus low (HUCCT1) proliferative CCA cells in vitro showed that both EGI1 and HUCCT1 incorporated more fatty acids (FA) than NHC, leading to increased triglyceride storage, also observed in Akt1/Nicd1-driven CCA mouse model. The highly proliferative EGI1 CCA cells showed greater uptake of very-low-density and HDLs than NHC and HUCCT1 CCA cells and increased cholesteryl ester content. The FA oxidation (FAO) and related proteome enrichment were specifically up-regulated in EGI1, and consequently, pharmacological blockade of FAO induced more pronounced inhibition of their tumorigenic capacity compared with HUCCT1. The expression of acyl-CoA dehydrogenase ACADM, the first enzyme involved in FAO, was increased in human CCA tissues and correlated with the proliferation marker PCNA.

CONCLUSIONS

Highly proliferative human CCA cells rely on lipid and lipoprotein uptake to fuel FA catabolism, suggesting that inhibition of FAO and/or lipid uptake could represent a therapeutic strategy for this CCA subclass.

Clinical stratification improves the diagnostic accuracy of small omics datasets within machine learning and genome-scale metabolic modelling methods

BACKGROUND

Recently, multi-omic machine learning architectures have been proposed for the early detection of cancer. However, for rare cancers and their associated small datasets, it is still unclear how to use the available multi-omics data to achieve a mechanistic prediction of cancer onset and progression, due to the limited data available. Hepatoblastoma is the most frequent liver cancer in infancy and childhood, and whose incidence has been lately increasing in several developed countries. Even though some studies have been conducted to understand the causes of its onset and discover potential biomarkers, the role of metabolic rewiring has not been investigated in depth so far.

METHODS

Here, we propose and implement an interpretable multi-omics pipeline that combines mechanistic knowledge from genome-scale metabolic models with machine learning algorithms, and we use it to characterise the underlying mechanisms controlling hepatoblastoma.

RESULTS AND CONCLUSIONS

While the obtained machine learning models generally present a high diagnostic classification accuracy, our results show that the type of omics combinations used as input to the machine learning models strongly affects the detection of important genes, reactions and metabolic pathways linked to hepatoblastoma. Our method also suggests that, in the context of computer-aided diagnosis of cancer, optimal diagnostic accuracy can be achieved by adopting a combination of omics that depends on the patient's clinical characteristics.

Cellular Genome-Scale Metabolic Modeling Identifies New Potential Drug Targets Against Hepatocellular Carcinoma

Genome-scale metabolic modeling (GEM) is one of the key approaches to unpack cancer metabolism and for discovery of new drug targets. In this study, we report the Transcriptional Regulated Flux Balance Analysis-CORE (TRFBA-), an algorithm for GEM using key growth-correlated reactions using hepatocellular carcinoma (HCC), an important global health burden, as a case study. We generated a HepG2 cell-specific GEM by integrating this cell line transcriptomic data with a generic human metabolic model to forecast potential drug targets for HCC. A total of 108 essential genes for growth were predicted by the TRFBA-CORE. These genes were enriched for metabolic pathways involved in cholesterol, sterol, and steroid biosynthesis. Furthermore, we silenced a predicted essential gene, 11-beta dehydrogenase hydroxysteroid type 2 (HSD11B2), in HepG2 cells resulting in a reduction in cell viability. To further identify novel potential drug targets in HCC, we examined the effect of nine drugs targeting the essential genes, and observed that most drugs inhibited the growth of HepG2 cells. Some of these drugs in this model performed better than Sorafenib, the first-line therapeutic against HCC. A HepG2 cell-specific GEM highlights sterol metabolism to be essential for cell growth. HSD11B2 downregulation results in lower cell growth. Most of the compounds, selected by drug repurposing approach, show a significant inhibitory effect on cell growth in a wide range of concentrations. These findings offer new molecular leads for drug discovery for hepatic cancer while also illustrating the importance of GEM and drug repurposing in cancer therapeutics innovation.

Review of Current Human Genome-Scale Metabolic Models for Brain Cancer and Neurodegenerative Diseases

Brain disorders represent 32% of the global disease burden, with 169 million Europeans affected. Constraint-based metabolic modelling and other approaches have been applied to predict new treatments for these and other diseases. Many recent studies focused on enhancing, among others, drug predictions by generating generic metabolic models of brain cells and on the contextualisation of the genome-scale metabolic models with expression data. Experimental flux rates were primarily used to constrain or validate the model inputs. Bi-cellular models were reconstructed to study the interaction between different cell types. This review highlights the evolution of genome-scale models for neurodegenerative diseases and glioma. We discuss the advantages and drawbacks of each approach and propose improvements, such as building bi-cellular models, tailoring the biomass formulations for glioma and refinement of the cerebrospinal fluid composition.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

System-level metabolic modeling facilitates unveiling metabolic signature in exceptional longevity

Although it is well known that metabolic control plays a crucial role in regulating the health span and life span of various organisms, little is known for the systems metabolic profile of centenarians, the paradigm of human healthy aging and longevity. Meanwhile, how to well characterize the system‐level metabolic states in an organism of interest remains to be a major challenge in systems metabolism research. To address this challenge and better understand the metabolic mechanisms of healthy aging, we developed a method of genome‐wide precision metabolic modeling (GPMM) which is able to quantitatively integrate transcriptome, proteome and kinetome data in predictive modeling of metabolic networks. Benchmarking analysis showed that GPMM successfully characterized metabolic reprogramming in the NCI‐60 cancer cell lines; it dramatically improved the performance of the modeling with an R² of 0.86 between the predicted and experimental measurements over the performance of existing methods. Using this approach, we examined the metabolic networks of a Chinese centenarian cohort and identified the elevated fatty acid oxidation (FAO) as the most significant metabolic feature in these long‐lived individuals. Evidence from serum metabolomics supports this observation. Given that FAO declines with normal aging and is impaired in many age‐related diseases, our study suggests that the elevated FAO has potential to be a novel signature of healthy aging of humans.

Compartmentalization of metabolism between cell types in multicellular organisms: a computational perspective

In multicellular organisms, metabolism is compartmentalized at many levels, including tissues and organs, different cell types, and subcellular compartments. Compartmentalization creates a coordinated homeostatic system where each compartment contributes to the production of energy and biomolecules the organism needs to carrying out specific metabolic tasks. Experimentally studying metabolic compartmentalization and metabolic interactions between cells and tissues in multicellular organisms is challenging at a systems level. However, recent progress in computational modeling provides an alternative approach to this problem. Here we discuss how integrating metabolic network modeling with omics data offers an opportunity to reveal metabolic states at the level of organs, tissues and, ultimately, individual cells. We review the current status of genome-scale metabolic network models in multicellular organisms, methods to study metabolic compartmentalization in silico, and insights gained from computational analyses. We also discuss outstanding challenges and provide perspectives for the future directions of the field.

Drug Target Prediction Using Context-Specific Metabolic Models Reconstructed from rFASTCORMICS

Metabolic modeling is a powerful computational tool to analyze metabolism. It has not only been used to identify metabolic rewiring strategies in cancer but also to predict drug targets and candidate drugs for repurposing. Here, we will elaborate on the reconstruction of context-specific metabolic models of cancer using rFASTCORMICS and the subsequent prediction of drugs for repurposing using our drug prediction workflow.

Discovering Essential Multiple Gene Effects Through Large Scale Optimization: An Application to Human Cancer Metabolism

Computational modelling of metabolic processes has proven to be a useful approach to formulate our knowledge and improve our understanding of core biochemical systems that are crucial to maintaining cellular functions. Towards understanding the broader role of metabolism on cellular decision-making in health and disease conditions, it is important to integrate the study of metabolism with other core regulatory systems and omics within the cell, including gene expression patterns. After quantitatively integrating gene expression profiles with a genome-scale reconstruction of human metabolism, we propose a set of combinatorial methods to reverse engineer gene expression profiles and to find pairs and higher-order combinations of genetic modifications that simultaneously optimize multi-objective cellular goals. This enables us to suggest classes of transcriptomic profiles that are most suitable to achieve given metabolic phenotypes. We demonstrate how our techniques are able to compute beneficial, neutral or "toxic" combinations of gene expression levels. We test our methods on nine tissue-specific cancer models, comparing our outcomes with the corresponding normal cells, identifying genes as targets for potential therapies. Our methods open the way to a broad class of applications that require an understanding of the interplay among genotype, metabolism, and cellular behaviour, at scale.

Computer-Aided Design for Identifying Anticancer Targets in Genome-Scale Metabolic Models of Colon Cancer

Simple Summary

Discovery of anticancer targets with minimal side effects is a major challenge in drug discovery and development. This study developed a fuzzy optimization framework for identifying anticancer targets. The framework was applied to identify not only gene regulator targets but also metabolite- and reaction-centric targets. The computational results show that the combination of a carbon metabolism target and any one-target gene that participates in the sphingolipid, glycerophospholipid, nucleotide, cholesterol biosynthesis, or pentose phosphate pathways is more effective for treatment than one-target inhibition is, and a two-target combination of 5-FU and folate supplement can improve cell viability, reduce metabolic deviation, and reduce side effects of normal cells.

Abstract

The efficient discovery of anticancer targets with minimal side effects is a major challenge in drug discovery and development. Early prediction of side effects is key for reducing development costs, increasing drug efficacy, and increasing drug safety. This study developed a fuzzy optimization framework for Identifying AntiCancer Targets (IACT) using constraint-based models. Four objectives were established to evaluate the mortality of treated cancer cells and to minimize side effects causing toxicity-induced tumorigenesis on normal cells and smaller metabolic perturbations. Fuzzy set theory was applied to evaluate potential side effects and investigate the magnitude of metabolic deviations in perturbed cells compared with their normal counterparts. The framework was applied to identify not only gene regulator targets but also metabolite- and reaction-centric targets. A nested hybrid differential evolution algorithm with a hierarchical fitness function was applied to solve multilevel IACT problems. The results show that the combination of a carbon metabolism target and any one-target gene that participates in the sphingolipid, glycerophospholipid, nucleotide, cholesterol biosynthesis, or pentose phosphate pathways is more effective for treatment than one-target inhibition is. A clinical antimetabolite drug 5-fluorouracil (5-FU) has been used to inhibit synthesis of deoxythymidine-5′-triphosphate for treatment of colorectal cancer. The computational results reveal that a two-target combination of 5-FU and a folate supplement can improve cell viability, reduce metabolic deviation, and reduce side effects of normal cells.

Argininosuccinate lyase is a metabolic vulnerability in breast development and cancer

Epithelial-to-mesenchymal transition (EMT) is fundamental to both normal tissue development and cancer progression. We hypothesized that EMT plasticity defines a range of metabolic phenotypes and that individual breast epithelial metabolic phenotypes are likely to fall within this phenotypic landscape. To determine EMT metabolic phenotypes, the metabolism of EMT was described within genome-scale metabolic models (GSMMs) using either transcriptomic or proteomic data from the breast epithelial EMT cell culture model D492. The ability of the different data types to describe breast epithelial metabolism was assessed using constraint-based modeling which was subsequently verified using ¹³C isotope tracer analysis. The application of proteomic data to GSMMs provided relatively higher accuracy in flux predictions compared to the transcriptomic data. Furthermore, the proteomic GSMMs predicted altered cholesterol metabolism and increased dependency on argininosuccinate lyase (ASL) following EMT which were confirmed in vitro using drug assays and siRNA knockdown experiments. The successful verification of the proteomic GSMMs afforded iBreast2886, a breast GSMM that encompasses the metabolic plasticity of EMT as defined by the D492 EMT cell culture model. Analysis of breast tumor proteomic data using iBreast2886 identified vulnerabilities within arginine metabolism that allowed prognostic discrimination of breast cancer patients on a subtype-specific level. Taken together, we demonstrate that the metabolic reconstruction iBreast2886 formalizes the metabolism of breast epithelial cell development and can be utilized as a tool for the functional interpretation of high throughput clinical data.

Genome-Scale Metabolic Modeling of the Human Microbiome in the Era of Personalized Medicine

The human microbiome plays an important role in human health and disease. Meta-omics analyses provide indispensable data for linking changes in microbiome composition and function to disease etiology. Yet, the lack of a mechanistic understanding of, e.g., microbiome-metabolome links hampers the translation of these findings into effective, novel therapeutics. Here, we propose metabolic modeling of microbial communities through constraint-based reconstruction and analysis (COBRA) as a complementary approach to meta-omics analyses. First, we highlight the importance of microbial metabolism in cardiometabolic diseases, inflammatory bowel disease, colorectal cancer, Alzheimer disease, and Parkinson disease. Next, we demonstrate that microbial community modeling can stratify patients and controls, mechanistically link microbes with fecal metabolites altered in disease, and identify host pathways affected by the microbiome. Finally, we outline our vision for COBRA modeling combined with meta-omics analyses and multivariate statistical analyses to inform and guide clinical trials, yield testable hypotheses, and ultimately propose novel dietary and therapeutic interventions. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

A metabolic modeling approach reveals promising therapeutic targets and antiviral drugs to combat COVID-19

In this study we have developed a method based on Flux Balance Analysis to identify human metabolic enzymes which can be targeted for therapeutic intervention against COVID-19. A literature search was carried out in order to identify suitable inhibitors of these enzymes, which were confirmed by docking calculations. In total, 10 targets and 12 bioactive molecules have been predicted. Among the most promising molecules we identified Triacsin C, which inhibits ACSL3, and which has been shown to be very effective against different viruses, including positive-sense single-stranded RNA viruses. Similarly, we also identified the drug Celgosivir, which has been successfully tested in cells infected with different types of viruses such as Dengue, Zika, Hepatitis C and Influenza. Finally, other drugs targeting enzymes of lipid metabolism, carbohydrate metabolism or protein palmitoylation (such as Propylthiouracil, 2-Bromopalmitate, Lipofermata, Tunicamycin, Benzyl Isothiocyanate, Tipifarnib and Lonafarnib) are also proposed.

Integration of machine learning and genome-scale metabolic modeling identifies multi-omics biomarkers for radiation resistance

Resistance to ionizing radiation, a first-line therapy for many cancers, is a major clinical challenge. Personalized prediction of tumor radiosensitivity is not currently implemented clinically due to insufficient accuracy of existing machine learning classifiers. Despite the acknowledged role of tumor metabolism in radiation response, metabolomics data is rarely collected in large multi-omics initiatives such as The Cancer Genome Atlas (TCGA) and consequently omitted from algorithm development. In this study, we circumvent the paucity of personalized metabolomics information by characterizing 915 TCGA patient tumors with genome-scale metabolic Flux Balance Analysis models generated from transcriptomic and genomic datasets. Metabolic biomarkers differentiating radiation-sensitive and -resistant tumors are predicted and experimentally validated, enabling integration of metabolic features with other multi-omics datasets into ensemble-based machine learning classifiers for radiation response. These multi-omics classifiers show improved classification accuracy, identify clinical patient subgroups, and demonstrate the utility of personalized blood-based metabolic biomarkers for radiation sensitivity. The integration of machine learning with genome-scale metabolic modeling represents a significant methodological advancement for identifying prognostic metabolite biomarkers and predicting radiosensitivity for individual patients.

PIWIL1 governs the crosstalk of cancer cell metabolism and immunosuppressive microenvironment in hepatocellular carcinoma

Altered energy metabolism of cancer cells shapes the immune cell response in the tumor microenvironment that facilitates tumor progression. Herein, we reported the novel of tumor cell-expressed Piwi Like RNA-Mediated Gene Silencing 1 (PIWIL1) in mediating the crosstalk of fatty acid metabolism and immune response of human hepatocellular carcinoma (HCC). PIWIL1 expression in HCC was increased compared to normal hepatic tissues and was positively correlated with the proliferation rate of HCC cell lines. PIWIL1 overexpression accelerated in vitro proliferation and in vivo growth of HCC tumors, while PIWIL1 knockdown showed opposite effects. PIWIL1 increased oxygen consumption and energy production via fatty acid metabolism without altering aerobic glycolysis. Inhibition of fatty acid metabolism abolished PIWIL1-induced HCC proliferation and growth. RNA-seq analysis revealed that immune system regulation might be involved, which was echoed by the experimental observation that PIWIL1-overexpressing HCC cells attracted myeloid-derived suppressor cells (MDSCs) into the tumor microenvironment. MDSCs depletion reduced the proliferation and growth of PIWIL1-overexpressing HCC tumors. Complement C3, whose secretion was induced by PIWIL1 in HCC cells, mediates the interaction of HCC cells with MDSCs by activated p38 MAPK signaling in MDSCs, which in turn initiated expression of immunosuppressive cytokine IL10. Neutralizing IL10 secretion reduced the immunosuppressive activity of MDSCs in the microenvironment of PIWIL1-overexpressing HCC. Taken together, our study unraveled the critical role of PIWIL1 in initiating the interaction of cancer cell metabolism and immune cell response in HCC. Tumor cells-expressed PIWIL1 may be a potential target for the development of novel HCC treatment.

Towards the routine use of in silico screenings for drug discovery using metabolic modelling

Currently, the development of new effective drugs for cancer therapy is not only hindered by development costs, drug efficacy, and drug safety but also by the rapid occurrence of drug resistance in cancer. Hence, new tools are needed to study the underlying mechanisms in cancer. Here, we discuss the current use of metabolic modelling approaches to identify cancer-specific metabolism and find possible new drug targets and drugs for repurposing. Furthermore, we list valuable resources that are needed for the reconstruction of cancer-specific models by integrating various available datasets with genome-scale metabolic reconstructions using model-building algorithms. We also discuss how new drug targets can be determined by using gene essentiality analysis, an in silico method to predict essential genes in a given condition such as cancer and how synthetic lethality studies could greatly benefit cancer patients by suggesting drug combinations with reduced side effects.

The level of oncogenic Ras determines the malignant transformation of Lkb1 mutant tissue in vivo

The genetic and metabolic heterogeneity of RAS-driven cancers has confounded therapeutic strategies in the clinic. To address this, rapid and genetically tractable animal models are needed that recapitulate the heterogeneity of RAS-driven cancers in vivo. Here, we generate a Drosophila melanogaster model of Ras/Lkb1 mutant carcinoma. We show that low-level expression of oncogenic Ras (RasLow) promotes the survival of Lkb1 mutant tissue, but results in autonomous cell cycle arrest and non-autonomous overgrowth of wild-type tissue. In contrast, high-level expression of oncogenic Ras (RasHigh) transforms Lkb1 mutant tissue resulting in lethal malignant tumors. Using simultaneous multiview light-sheet microcopy, we have characterized invasion phenotypes of Ras/Lkb1 tumors in living larvae. Our molecular analysis reveals sustained activation of the AMPK pathway in malignant Ras/Lkb1 tumors, and demonstrate the genetic and pharmacologic dependence of these tumors on CaMK-activated Ampk. We further show that LKB1 mutant human lung adenocarcinoma patients with high levels of oncogenic KRAS exhibit worse overall survival and increased AMPK activation. Our results suggest that high levels of oncogenic KRAS is a driving event in the malignant transformation of LKB1 mutant tissue, and uncovers a vulnerability that may be used to target this aggressive genetic subset of RAS-driven tumors.

Serine Biosynthesis Is a Metabolic Vulnerability in IDH2-Driven Breast Cancer Progression

Cancer-specific metabolic phenotypes and their vulnerabilities represent a viable area of cancer research. In this study, we explored the association of breast cancer subtypes with different metabolic phenotypes and identified isocitrate dehydrogenase 2 (IDH2) as a key player in triple-negative breast cancer (TNBC) and HER2. Functional assays combined with mass spectrometry-based analyses revealed the oncogenic role of IDH2 in cell proliferation, anchorage-independent growth, glycolysis, mitochondrial respiration, and antioxidant defense. Genome-scale metabolic modeling identified phosphoglycerate dehydrogenase (PHGDH) and phosphoserine aminotransferase (PSAT1) as the synthetic dosage lethal (SDL) partners of IDH2. In agreement, CRISPR-Cas9 knockout of PHGDH and PSAT1 showed the essentiality of serine biosynthesis proteins in IDH2-high cells. The clinical significance of the SDL interaction was supported by patients with IDH2-high/PHGDH-low tumors, who exhibited longer survival than patients with IDH2-high/PHGDH-high tumors. Furthermore, PHGDH inhibitors were effective in treating IDH2-high cells in vitro and in vivo. Altogether, our study creates a new link between two known cancer regulators and emphasizes PHGDH as a promising target for TNBC with IDH2 overexpression. SIGNIFICANCE: These findings highlight the metabolic dependence of IDH2 on the serine biosynthesis pathway, adding an important layer to the connection between TCA cycle and glycolysis, which can be translated into novel targeted therapies.

Opportunities at the Interface of Network Science and Metabolic Modeling

Metabolism plays a central role in cell physiology because it provides the molecular machinery for growth. At the genome-scale, metabolism is made up of thousands of reactions interacting with one another. Untangling this complexity is key to understand how cells respond to genetic, environmental, or therapeutic perturbations. Here we discuss the roles of two complementary strategies for the analysis of genome-scale metabolic models: Flux Balance Analysis (FBA) and network science. While FBA estimates metabolic flux on the basis of an optimization principle, network approaches reveal emergent properties of the global metabolic connectivity. We highlight how the integration of both approaches promises to deliver insights on the structure and function of metabolic systems with wide-ranging implications in discovery science, precision medicine and industrial biotechnology.

Exploring gene knockout strategies to identify potential drug targets using genome-scale metabolic models

Research on new cancer drugs is performed either through gene knockout studies or phenotypic screening of drugs in cancer cell-lines. Both of these approaches are costly and time-consuming. Computational framework, e.g., genome-scale metabolic models (GSMMs), could be a good alternative to find potential drug targets. The present study aims to investigate the applicability of gene knockout strategies to be used as the finding of drug targets using GSMMs. We performed single-gene knockout studies on existing GSMMs of the NCI-60 cell-lines obtained from 9 tissue types. The metabolic genes responsible for the growth of cancerous cells were identified and then ranked based on their cellular growth reduction. The possible growth reduction mechanisms, which matches with the gene knockout results, were described. Gene ranking was used to identify potential drug targets, which reduce the growth rate of cancer cells but not of the normal cells. The gene ranking results were also compared with existing shRNA screening data. The rank-correlation results for most of the cell-lines were not satisfactory for a single-gene knockout, but it played a significant role in deciding the activity of drug against cell proliferation, whereas multiple gene knockout analysis gave better correlation results. We validated our theoretical results experimentally and showed that the drugs mitotane and myxothiazol can inhibit the growth of at least four cell-lines of NCI-60 database.

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease

A multitude of factors contribute to complex diseases and can be measured with ‘omics’ methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources ‘Human metabolism’, ‘Gut microbiome’, ‘Disease’, ‘Nutrition’, and ‘ReconMaps’. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH’s unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

Common biochemical properties of metabolic genes recurrently dysregulated in tumors

BACKGROUND

Tumor initiation and progression are associated with numerous metabolic alterations. However, the biochemical drivers and constraints that contribute to metabolic gene dysregulation are unclear.

METHODS

Here, we present MetOncoFit, a computational model that integrates 142 metabolic features that can impact tumor fitness, including enzyme catalytic activity, pathway association, network topology, and reaction flux. MetOncoFit uses genome-scale metabolic modeling and machine-learning to quantify the relative importance of various metabolic features in predicting cancer metabolic gene expression, copy number variation, and survival data.

RESULTS

Using MetOncoFit, we performed a meta-analysis of 9 cancer types and over 4500 samples from TCGA, Prognoscan, and COSMIC tumor databases. MetOncoFit accurately predicted enzyme differential expression and its impact on patient survival using the 142 attributes of metabolic enzymes. Our analysis revealed that enzymes with high catalytic activity were frequently upregulated in many tumors and associated with poor survival. Topological analysis also identified specific metabolites that were hot spots of dysregulation.

CONCLUSIONS

MetOncoFit integrates a broad range of datasets to understand how biochemical and topological features influence metabolic gene dysregulation across various cancer types. MetOncoFit was able to achieve significantly higher accuracy in predicting differential expression, copy number variation, and patient survival than traditional modeling approaches. Overall, MetOncoFit illuminates how enzyme activity and metabolic network architecture influences tumorigenesis.

Network-based metabolic characterization of renal cell carcinoma

An emerging hallmark of cancer is metabolic reprogramming, which presents opportunities for cancer diagnosis and treatment based on metabolism. We performed a comprehensive metabolic network analysis of major renal cell carcinoma (RCC) subtypes including clear cell, papillary and chromophobe by integrating transcriptomic data with the human genome-scale metabolic model to understand the coordination of metabolic pathways in cancer cells. We identified metabolic alterations of each subtype with respect to tumor-adjacent normal samples and compared them to understand the differences between subtypes. We found that genes of amino acid metabolism and redox homeostasis are significantly altered in RCC subtypes. Chromophobe showed metabolic divergence compared to other subtypes with upregulation of genes involved in glutamine anaplerosis and aspartate biosynthesis. A difference in transcriptional regulation involving HIF1A is observed between subtypes. We identified E2F1 and FOXM1 as other major transcriptional activators of metabolic genes in RCC. Further, the co-expression pattern of metabolic genes in each patient showed the variations in metabolism within RCC subtypes. We also found that co-expression modules of each subtype have tumor stage-specific behavior, which may have clinical implications.

Inferring Metabolic Flux from Time-Course Metabolomics

The metabolic activity of a mammalian cell changes dynamically over time and is tied to the changing metabolic demands of cellular processes such as cell differentiation and proliferation. While experimental tools like time-course metabolomics and flux tracing can measure the dynamics of a few pathways, they are unable to infer fluxes at the whole network level. To address this limitation, we have developed the Dynamic Flux Activity (DFA) algorithm, a genome-scale modeling approach that uses time-course metabolomics to predict dynamic flux rewiring during transitions between metabolic states. This chapter provides a protocol for applying DFA to characterize the dynamic metabolic activity of various cancer cell lines.

Genome-Scale Metabolic Modeling with Protein Expressions of Normal and Cancerous Colorectal Tissues for Oncogene Inference

Although cancer has historically been regarded as a cell proliferation disorder, it has recently been considered a metabolic disease. The first discovery of metabolic alterations in cancer cells refers to Otto Warburg’s observations. Cancer metabolism results in alterations in metabolic fluxes that are evident in cancer cells compared with most normal tissue cells. This study applied protein expressions of normal and cancer cells to reconstruct two tissue-specific genome-scale metabolic models. Both models were employed in a tri-level optimization framework to infer oncogenes. Moreover, this study also introduced enzyme pseudo-coding numbers in the gene association expression to avoid performing posterior decision-making that is necessary for the reaction-based method. Colorectal cancer (CRC) was the topic of this case study, and 20 top-ranked oncogenes were determined. Notably, these dysregulated genes were involved in various metabolic subsystems and compartments. We found that the average similarity ratio for each dysregulation is higher than 98%, and the extent of similarity for flux changes is higher than 93%. On the basis of surveys of PubMed and GeneCards, these oncogenes were also investigated in various carcinomas and diseases. Most dysregulated genes connect to catalase that acts as a hub and connects protein signaling pathways, such as those involving TP53, mTOR, AKT1, MAPK1, EGFR, MYC, CDK8, and RAS family.

Identification of Cancer-associated metabolic vulnerabilities by modeling multi-objective optimality in metabolism

BACKGROUND

Cancer cells undergo global reprogramming of cellular metabolism to satisfy demands of energy and biomass during proliferation and metastasis. Computational modeling of genome-scale metabolic models is an effective approach for designing new therapeutics targeting dysregulated cancer metabolism by identifying metabolic enzymes crucial for satisfying metabolic goals of cancer cells, but nearly all previous studies neglect the existence of metabolic demands other than biomass synthesis and trade-offs between these contradicting metabolic demands. It is thus necessary to develop computational models covering multiple metabolic objectives to study cancer metabolism and identify novel metabolic targets.

METHODS

We developed a multi-objective optimization model for cancer cell metabolism at genome-scale and an integrated, data-driven workflow for analyzing the Pareto optimality of this model in achieving multiple metabolic goals and identifying metabolic enzymes crucial for maintaining cancer-associated metabolic phenotypes. Using this workflow, we constructed cell line-specific models for a panel of cancer cell lines and identified lists of metabolic targets promoting or suppressing cancer cell proliferation or the Warburg Effect. The targets were then validated using knockdown and over-expression experiments in cultured cancer cell lines.

RESULTS

We found that the multi-objective optimization model correctly predicted phenotypes including cell growth rates, essentiality of metabolic genes and cell line specific sensitivities to metabolic perturbations. To our surprise, metabolic enzymes promoting proliferation substantially overlapped with those suppressing the Warburg Effect, suggesting that simply targeting the overlapping enzymes may lead to complicated outcomes. We also identified lists of metabolic enzymes important for maintaining rapid proliferation or high Warburg Effect while having little effect on the other. The importance of these enzymes in cancer metabolism predicted by the model was validated by their association with cancer patient survival and knockdown and overexpression experiments in a variety of cancer cell lines.

CONCLUSIONS

These results confirm this multi-objective optimization model as a novel and effective approach for studying trade-off between metabolic demands of cancer cells and identifying cancer-associated metabolic vulnerabilities, and suggest novel metabolic targets for cancer treatment.

Pan-cancer analysis of the metabolic reaction network

Metabolic reprogramming is considered a hallmark of malignant transformation. However, it is not clear whether the network of metabolic reactions expressed by cancers of different origin differ from each other or from normal human tissues. In this study, we reconstructed functional and connected genome-scale metabolic models for 917 primary tumor samples across 13 types based on the probability of expression for 3765 reference metabolic genes in the sample. This network-centric approach revealed that tumor metabolic networks are largely similar in terms of accounted reactions, despite diversity in the expression of the associated genes. On average, each network contained 4721 reactions, of which 74% were core reactions (present in >95% of all models). Whilst 99.3% of the core reactions were classified as housekeeping also in normal tissues, we identified reactions catalyzed by ARG2, RHAG, SLC6 and SLC16 family gene members, and PTGS1 and PTGS2 as core exclusively in cancer. These findings were subsequently replicated in an independent validation set of 3388 genome-scale metabolic models. The remaining 26% of the reactions were contextual reactions. Their inclusion was dependent in one case (GLS2) on the absence of TP53 mutations and in 94.6% of cases on differences in cancer types. This dependency largely resembled differences in expression patterns in the corresponding normal tissues, with some exceptions like the presence of the NANP-encoded reaction in tumors not from the female reproductive system or of the SLC5A9-encoded reaction in kidney-pancreatic-colorectal tumors. In conclusion, tumors expressed a metabolic network virtually overlapping the matched normal tissues, raising the possibility that metabolic reprogramming simply reflects cancer cell plasticity to adapt to varying conditions thanks to redundancy and complexity of the underlying metabolic networks. At the same time, the here uncovered exceptions represent a resource to identify selective liabilities of tumor metabolism.

MOOMIN - Mathematical explOration of 'Omics data on a MetabolIc Network

MOTIVATION

Analysis of differential expression of genes is often performed to understand how the metabolic activity of an organism is impacted by a perturbation. However, because the system of metabolic regulation is complex and all changes are not directly reflected in the expression levels, interpreting these data can be difficult.

RESULTS

In this work, we present a new algorithm and computational tool that uses a genome-scale metabolic reconstruction to infer metabolic changes from differential expression data. Using the framework of constraint-based analysis, our method produces a qualitative hypothesis of a change in metabolic activity. In other words, each reaction of the network is inferred to have increased, decreased, or remained unchanged in flux. In contrast to similar previous approaches, our method does not require a biological objective function and does not assign on/off activity states to genes. An implementation is provided and it is available online. We apply the method to three published datasets to show that it successfully accomplishes its two main goals: confirming or rejecting metabolic changes suggested by differentially expressed genes based on how well they fit in as parts of a coordinated metabolic change, as well as inferring changes in reactions whose genes did not undergo differential expression.

AVAILABILITY AND IMPLEMENTATION

github.com/htpusa/moomin.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TCSBN: a database of tissue and cancer specific biological networks

Biological networks provide new opportunities for understanding the cellular biology in both health and disease states. We generated tissue specific integrated networks (INs) for liver, muscle and adipose tissues by integrating metabolic, regulatory and protein-protein interaction networks. We also generated human co-expression networks (CNs) for 46 normal tissues and 17 cancers to explore the functional relationships between genes as well as their relationships with biological functions, and investigate the overlap between functional and physical interactions provided by CNs and INs, respectively. These networks can be employed in the analysis of omics data, provide detailed insight into disease mechanisms by identifying the key biological components and eventually can be used in the development of efficient treatment strategies. Moreover, comparative analysis of the networks may allow for the identification of tissue-specific targets that can be used in the development of drugs with the minimum toxic effect to other human tissues. These context-specific INs and CNs are presented in an interactive website http://inetmodels.com without any limitation.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

Metabolic reprogramming dynamics in tumor spheroids: Insights from a multicellular, multiscale model

Mathematical modeling provides the predictive ability to understand the metabolic reprogramming and complex pathways that mediate cancer cells’ proliferation. We present a mathematical model using a multiscale, multicellular approach to simulate avascular tumor growth, applied to pancreatic cancer. The model spans three distinct spatial and temporal scales. At the extracellular level, reaction diffusion equations describe nutrient concentrations over a span of seconds. At the cellular level, a lattice-based energy driven stochastic approach describes cellular phenomena including adhesion, proliferation, viability and cell state transitions, occurring on the timescale of hours. At the sub-cellular level, we incorporate a detailed kinetic model of intracellular metabolite dynamics on the timescale of minutes, which enables the cells to uptake and excrete metabolites and use the metabolites to generate energy and building blocks for cell growth. This is a particularly novel aspect of the model. Certain defined criteria for the concentrations of intracellular metabolites lead to cancer cell growth, proliferation or death. Overall, we model the evolution of the tumor in both time and space. Starting with a cluster of tumor cells, the model produces an avascular tumor that quantitatively and qualitatively mimics experimental measurements of multicellular tumor spheroids. Through our model simulations, we can investigate the response of individual intracellular species under a metabolic perturbation and investigate how that response contributes to the response of the tumor as a whole. The predicted response of intracellular metabolites under various targeted strategies are difficult to resolve with experimental techniques. Thus, the model can give novel predictions as to the response of the tumor as a whole, identifies potential therapies to impede tumor growth, and predicts the effects of those therapeutic strategies. In particular, the model provides quantitative insight into the dynamic reprogramming of tumor cells at the intracellular level in response to specific metabolic perturbations. Overall, the model is a useful framework to study targeted metabolic strategies for inhibiting tumor growth.

Cancer cells expertly alter their metabolism in order to sustain growth, a hallmark of cancer. Quantitative details about this metabolic reprogramming are difficult to obtain without the use of predictive mathematical models. Here, we present a robust computational model of avascular tumor growth. The novel aspect of this work lies in the incorporation of a detailed model of the dynamics of metabolism within each individual cell, which directly influence growth of the multicellular tumor as a whole. We apply the model to simulate how the tumor grows in space and time and to predict how the tumor responds to targeted inhibition of specific intracellular metabolic reactions. Our results show, first-hand, the dynamic metabolic reprogramming that occurs in cancer cells. Specifically, the model provides insight into how the cells alter their metabolism to compensate for the loss of a nutrient by exploiting alternative pathways for continued tumor growth. Our work provides a quantitative tool for identifying the impact of cellular and sub-cellular features on the evolution of a tumor. This framework is useful for developing potential cancer therapies, complementing experimental studies.

Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.

A Prob(e)able Route to Lysine Acylation

Non-enzymatic modification of proteins by acyl-CoA species involved in intermediary metabolism is a possible explanation for widespread protein acylation. In this issue, Kulkarni et al. (2017) develop a set of chemoproteomic probes to interrogate the role of malonyl-CoA in mediating protein malonylation and find malonylation influences glycolysis in cancer cells.

A benchmark-driven approach to reconstruct metabolic networks for studying cancer metabolism

Genome-scale metabolic modeling has emerged as a promising way to study the metabolic alterations underlying cancer by identifying novel drug targets and biomarkers. To date, several computational methods have been developed to integrate high-throughput data with existing human metabolic reconstructions to generate context-specific cancer metabolic models. Despite a number of studies focusing on benchmarking the context-specific algorithms, no quantitative assessment has been made to compare the predictive performance of these methods. Here, we integrated various and different datasets used in previous works to design a quantitative platform to examine functional and consistency performance of several existing genome-scale cancer modeling approaches. Next, we used the results obtained here to develop a method for the reconstruction of context-specific metabolic models. We then compared the predictive power and consistency of networks generated by our method to other computational approaches investigated here. Our results showed a satisfactory performance of the developed method in most of the benchmarks. This benchmarking platform is of particular use in algorithm selection and assessing the performance of newly developed algorithms. More importantly, it can serve as guidelines for designing and developing new methods focusing on weaknesses and strengths of existing algorithms.

Current Status and Future Prospects of Clinically Exploiting Cancer-specific Metabolism-Why Is Tumor Metabolism Not More Extensively Translated into Clinical Targets and Biomarkers?

Tumor cells exhibit a specialized metabolism supporting their superior ability for rapid proliferation, migration, and apoptotic evasion. It is reasonable to assume that the specific metabolic needs of the tumor cells can offer an array of therapeutic windows as pharmacological disturbance may derail the biochemical mechanisms necessary for maintaining the tumor characteristics, while being less important for normally proliferating cells. In addition, the specialized metabolism may leave a unique metabolic signature which could be used clinically for diagnostic or prognostic purposes. Quantitative global metabolic profiling (metabolomics) has evolved over the last two decades. However, despite the technology’s present ability to measure 1000s of endogenous metabolites in various clinical or biological specimens, there are essentially no examples of metabolomics investigations being translated into actual utility in the cancer clinic. This review investigates the current efforts of using metabolomics as a tool for translation of tumor metabolism into the clinic and further seeks to outline paths for increasing the momentum of using tumor metabolism as a biomarker and drug target opportunity.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

Exploring candidate biomarkers for lung and prostate cancers using gene expression and flux variability analysis

Genome-scale metabolic models have provided valuable resources for exploring changes in metabolism under normal and cancer conditions. However, metabolism itself is strongly linked to gene expression, so integration of gene expression data into metabolic models might improve the detection of genes involved in the control of tumor progression. Herein, we considered gene expression data as extra constraints to enhance the predictive powers of metabolic models. We reconstructed genome-scale metabolic models for lung and prostate, under normal and cancer conditions to detect the major genes associated with critical subsystems during tumor development. Furthermore, we utilized gene expression data in combination with an information theory-based approach to reconstruct co-expression networks of the human lung and prostate in both cohorts. Our results revealed 19 genes as candidate biomarkers for lung and prostate cancer cells. This study also revealed that the development of a complementary approach (integration of gene expression and metabolic profiles) could lead to proposing novel biomarkers and suggesting renovated cancer treatment strategies which have not been possible to detect using either of the methods alone.