Modeling cancer metabolism on a genome scale

Targeting tumour metabolism in melanoma to enhance response to immune checkpoint inhibition: A balancing act

Immune checkpoint inhibition has transformed the treatment landscape of advanced melanoma and long-term survival of patients is now possible. However, at least half of the patients do not benefit sufficiently. Metabolic reprogramming is a hallmark of cancer cells and may contribute to both tumour growth and immune evasion by the tumour. Preclinical studies have indeed demonstrated that modulating tumour metabolism can reduce tumour growth while improving the functionality of immune cells. Since metabolic pathways are commonly shared between immune and tumour cells, it is essential to understand how modulating tumour metabolism in patients influences the intricate balance of pro-and anti-tumour immune effects in the tumour microenvironment. The key question is whether modulating tumour metabolism can inhibit tumour cell growth as well as facilitate an anti-tumour immune response. Here, we review current knowledge on the effect of tumour metabolism on the immune response in melanoma. We summarise metabolic pathways in melanoma and non-cancerous cells in the tumour microenvironment and discuss models and techniques available to study the metabolic-immune interaction. Finally, we discuss clinical use of these techniques to improve our understanding of how metabolic interventions can tip the balance towards a favourable, immune permissive microenvironment in melanoma patients.

Genome-scale metabolic network of human carotid plaque reveals the pivotal role of glutamine/glutamate metabolism in macrophage modulating plaque inflammation and vulnerability

BACKGROUND

Metabolism is increasingly recognized as a key regulator of the function and phenotype of the primary cellular constituents of the atherosclerotic vascular wall, including endothelial cells, smooth muscle cells, and inflammatory cells. However, a comprehensive analysis of metabolic changes associated with the transition of plaque from a stable to a hemorrhaged phenotype is lacking.

METHODS

In this study, we integrated two large mRNA expression and protein abundance datasets (BIKE, n = 126; MaasHPS, n = 43) from human atherosclerotic carotid artery plaque to reconstruct a genome-scale metabolic network (GEM). Next, the GEM findings were linked to metabolomics data from MaasHPS, providing a comprehensive overview of metabolic changes in human plaque.

RESULTS

Our study identified significant changes in lipid, cholesterol, and inositol metabolism, along with altered lysosomal lytic activity and increased inflammatory activity, in unstable plaques with intraplaque hemorrhage (IPH+) compared to non-hemorrhaged (IPH-) plaques. Moreover, topological analysis of this network model revealed that the conversion of glutamine to glutamate and their flux between the cytoplasm and mitochondria were notably compromised in hemorrhaged plaques, with a significant reduction in overall glutamate levels in IPH+ plaques. Additionally, reduced glutamate availability was associated with an increased presence of macrophages and a pro-inflammatory phenotype in IPH+ plaques, suggesting an inflammation-prone microenvironment.

CONCLUSIONS

This study is the first to establish a robust and comprehensive GEM for atherosclerotic plaque, providing a valuable resource for understanding plaque metabolism. The utility of this GEM was illustrated by its ability to reliably predict dysregulation in the cholesterol hydroxylation, inositol metabolism, and the glutamine/glutamate pathway in rupture-prone hemorrhaged plaques, a finding that may pave the way to new diagnostic or therapeutic measures.

Cell-specific genome-scale metabolic modeling of SARS-CoV-2-infected lung to identify antiviral enzymes

Computational systems biology plays a key role in the discovery of suitable antiviral targets. We designed a cell-specific, constraint-based modeling technique for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected lungs. We used the gene sequence of the alpha variant of SARS-CoV-2 to build a viral biomass reaction (VBR). We also used the mass proportion of lipids between the viral biomass and its host cell to estimate the stoichiometric coefficients of viral lipids in the reaction. We then integrated the VBR, the gene expression of the alpha variant of SARS-CoV-2, and the generic human metabolic network Recon3D to reconstruct a cell-specific genome-scale metabolic model. An antiviral target discovery (AVTD) platform was introduced using this model to identify therapeutic drug targets for combating COVID-19. The AVTD platform not only identified antiviral genes for eliminating viral replication but also predicted side effects of treatments. Our computational results revealed that knocking out dihydroorotate dehydrogenase (DHODH) might reduce the synthesis rate of cytidine-5'-triphosphate and uridine-5'-triphosphate, which terminate the viral building blocks of DNA and RNA for SARS-CoV-2 replication. Our results also indicated that DHODH is a promising antiviral target that causes minor side effects, which is consistent with the results of recent reports. Moreover, we discovered that the genes that participate in the de novo biosynthesis of glycerophospholipids and ceramides become unidentifiable if the VBR does not involve the stoichiometry of lipids.

Genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer

Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment (TME), which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the TME. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

Flux estimation analysis systematically characterizes the metabolic shifts of the central metabolism pathway in human cancer

Introduction

Glucose and glutamine are major carbon and energy sources that promote the rapid proliferation of cancer cells. Metabolic shifts observed on cell lines or mouse models may not reflect the general metabolic shifts in real human cancer tissue.

Method

In this study, we conducted a computational characterization of the flux distribution and variations of the central energy metabolism and key branches in a pan-cancer analysis, including the glycolytic pathway, production of lactate, tricarboxylic acid (TCA) cycle, nucleic acid synthesis, glutaminolysis, glutamate, glutamine, and glutathione metabolism, and amino acid synthesis, in 11 cancer subtypes and nine matched adjacent normal tissue types using TCGA transcriptomics data.

Result

Our analysis confirms the increased influx in glucose uptake and glycolysis and decreased upper part of the TCA cycle, i.e., the Warburg effect, in almost all the analyzed cancer. However, increased lactate production and the second half of the TCA cycle were only seen in certain cancer types. More interestingly, we failed to detect significantly altered glutaminolysis in cancer tissues compared to their adjacent normal tissues. A systems biology model of metabolic shifts through cancer and tissue types is further developed and analyzed. We observed that (1) normal tissues have distinct metabolic phenotypes; (2) cancer types have drastically different metabolic shifts compared to their adjacent normal controls; and (3) the different shifts in tissue-specific metabolic phenotypes result in a converged metabolic phenotype through cancer types and cancer progression.

Discussion

This study strongly suggests the possibility of having a unified framework for studies of cancer-inducing stressors, adaptive metabolic reprogramming, and cancerous behaviors.

Crosstalk between Metabolite Production and Signaling Activity in Breast Cancer

The reprogramming of metabolism is a recognized cancer hallmark. It is well known that different signaling pathways regulate and orchestrate this reprogramming that contributes to cancer initiation and development. However, recent evidence is accumulating, suggesting that several metabolites could play a relevant role in regulating signaling pathways. To assess the potential role of metabolites in the regulation of signaling pathways, both metabolic and signaling pathway activities of Breast invasive Carcinoma (BRCA) have been modeled using mechanistic models. Gaussian Processes, powerful machine learning methods, were used in combination with SHapley Additive exPlanations (SHAP), a recent methodology that conveys causality, to obtain potential causal relationships between the production of metabolites and the regulation of signaling pathways. A total of 317 metabolites were found to have a strong impact on signaling circuits. The results presented here point to the existence of a complex crosstalk between signaling and metabolic pathways more complex than previously was thought.

Evaluation of Lung Cancer Patient Response to First-Line Chemotherapy by Integration of Tumor Core Biopsy Metabolomics with Multiscale Modeling

The standard of care for intermediate (Stage II) and advanced (Stages III and IV) non-small cell lung cancer (NSCLC) involves chemotherapy with taxane/platinum derivatives, with or without radiation. Ideally, patients would be screened a priori to allow non-responders to be initially treated with second-line therapies. This evaluation is non-trivial, however, since tumors behave as complex multiscale systems. To address this need, this study employs a multiscale modeling approach to evaluate first-line chemotherapy response of individual patient tumors based on metabolomic analysis of tumor core biopsies obtained during routine clinical evaluation. Model parameters were calculated for a patient cohort as a function of these metabolomic profiles, previously obtained from high-resolution 2DLC-MS/MS analysis. Evaluation metrics were defined to classify patients as Disease-Control (DC) [encompassing complete-response (CR), partial-response (PR), and stable-disease (SD)] and Progressive-Disease (PD) following first-line chemotherapy. Response was simulated for each patient and compared to actual response. The results show that patient classifications were significantly separated from each other, and also when grouped as DC vs. PD and as CR/PR vs. SD/PD, by fraction of initial tumor radius metric at 6 days post simulated bolus drug injection. This study shows that patient first-line chemotherapy response can in principle be evaluated from multiscale modeling integrated with tumor tissue metabolomic data, offering a first step towards individualized lung cancer treatment prognosis.

Mathematical Modeling of Eicosanoid Metabolism in Macrophage Cells: Cybernetic Framework Combined with Novel Information-Theoretic Approaches

Cellular response to inflammatory stimuli leads to the production of eicosanoids—prostanoids (PRs) and leukotrienes (LTs)—and signaling molecules—cytokines and chemokines—by macrophages. Quantitative modeling of the inflammatory response is challenging owing to a lack of knowledge of the complex regulatory processes involved. Cybernetic models address these challenges by utilizing a well-defined cybernetic goal and optimizing a coarse-grained model toward this goal. We developed a cybernetic model to study arachidonic acid (AA) metabolism, which included two branches, PRs and LTs. We utilized a priori biological knowledge to define the branch-specific cybernetic goals for PR and LT branches as the maximization of TNFα and CCL2, respectively. We estimated the model parameters by fitting data from three experimental conditions. With these parameters, we were able to capture a novel fourth independent experimental condition as part of the model validation. The cybernetic model enhanced our understanding of enzyme dynamics by predicting their profiles. The success of the model implies that the cell regulates the synthesis and activity of the associated enzymes, through cybernetic control variables, to accomplish the chosen biological goal. The results indicated that the dominant metabolites are PGD2 (a PR) and LTB4 (an LT), aligning with their corresponding known prominent biological roles during inflammation. Using heuristic arguments, we also infer that eicosanoid overproduction can lead to increased secretion of cytokines/chemokines. This novel model integrates mechanistic knowledge, known biological understanding of signaling pathways, and data-driven methods to study the dynamics of eicosanoid metabolism.

Ensemble-based genome-scale modeling predicts metabolic differences between macrophage subtypes in colorectal cancer

1Colorectal cancer (CRC) shows high incidence and mortality, partly due to the tumor microenvironment, which is viewed as an active promoter of disease progression. Macrophages are among the most abundant cells in the tumor microenvironment. These immune cells are generally categorized as M1, with inflammatory and anti-cancer properties, or M2, which promote tumor proliferation and survival. Although the M1/M2 subclassification scheme is strongly influenced by metabolism, the metabolic divergence between the subtypes remains poorly understood. Therefore, we generated a suite of computational models that characterize the M1- and M2-specific metabolic states. Our models show key differences between the M1 and M2 metabolic networks and capabilities. We leverage the models to identify metabolic perturbations that cause the metabolic state of M2 macrophages to more closely resemble M1 cells. Overall, this work increases understanding of macrophage metabolism in CRC and elucidates strategies to promote the metabolic state of anti-tumor macrophages.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Systematic diet composition swap in a mouse genome-scale metabolic model reveals determinants of obesogenic diet metabolism in liver cancer

Dietary nutrient availability and gene expression, together, influence tissue metabolic activity. Here, we explore whether altering dietary nutrient composition in the context of mouse liver cancer suffices to overcome chronic gene expression changes that arise from tumorigenesis and western-style diet (WD). We construct a mouse genome-scale metabolic model and estimate metabolic fluxes in liver tumors and non-tumoral tissue after computationally varying the composition of input diet. This approach, called Systematic Diet Composition Swap (SyDiCoS), revealed that, compared to a control diet, WD increases production of glycerol and succinate irrespective of specific tissue gene expression patterns. Conversely, differences in fatty acid utilization pathways between tumor and non-tumor liver are amplified with WD by both dietary carbohydrates and lipids together. Our data suggest that combined dietary component modifications may be required to normalize the distinctive metabolic patterns that underlie selective targeting of tumor metabolism.

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.

TumorMet: A repository of tumor metabolic networks derived from context-specific Genome-Scale Metabolic Models

Studies about the metabolic alterations during tumorigenesis have increased our knowledge of the underlying mechanisms and consequences, which are important for diagnostic and therapeutic investigations. In this scenario and in the era of systems biology, metabolic networks have become a powerful tool to unravel the complexity of the cancer metabolic machinery and the heterogeneity of this disease. Here, we present TumorMet, a repository of tumor metabolic networks extracted from context-specific Genome-Scale Metabolic Models, as a benchmark for graph machine learning algorithms and network analyses. This repository has an extended scope for use in graph classification, clustering, community detection, and graph embedding studies. Along with the data, we developed and provided Met2Graph, an R package for creating three different types of metabolic graphs, depending on the desired nodes and edges: Metabolites-, Enzymes-, and Reactions-based graphs. This package allows the easy generation of datasets for downstream analysis.

Measurement(s)	gene expression, metabolic relationships
Technology Type(s)	Genome Scale Metabolic Models; Computational network biology
Sample Characteristic - Organism	Homo sapiens

Pharmacometabolomics-Based Translational Biomarkers: How to Navigate the Data Ocean

Metabolome is the end point of the genome-environment interplay, and enables an important holistic overview of individual adaptability and host responses to environmental, ecological, as well as endogenous changes such as disease. Pharmacometabolomics is the application of metabolome knowledge to decipher the mechanisms of interindividual and intraindividual variations in drug efficacy and safety. Pharmacometabolomics also contributes to prediction of drug treatment outcomes on the basis of baseline (predose) and postdose metabotypes through mathematical modeling. Thus, pharmacometabolomics is a strong asset for a diverse community of stakeholders interested in theory and practice of evidence-based and precision/personalized medicine: academic researchers, public health scholars, health professionals, pharmaceutical, diagnostics, and biotechnology industries, among others. In this expert review, we discuss pharmacometabolomics in four contexts: (1) an interdisciplinary omics tool and field to map the mechanisms and scale of interindividual variability in drug effects, (2) discovery and development of translational biomarkers, (3) advance digital biomarkers, and (4) empower drug repurposing, a field that is increasingly proving useful in the current era of Covid-19. As the applications of pharmacometabolomics are growing rapidly in the current postgenome era, next-generation proteomics and metabolomics follow the example of next-generation sequencing analyses. Pharmacometabolomics can also empower data reliability and reproducibility through multiomics integration strategies, which use each data layer to correct, connect with, and inform each other. Finally, we underscore here that contextual data remain crucial for precision medicine and drug development that stand the test of time and clinical relevance.

Fuzzy optimization for identifying anti-cancer targets with few side effects in constraint-based models of head and neck cancer

Computer-aided methods can be used to screen potential candidate targets and to reduce the time and cost of drug development. In most of these methods, synthetic lethality is used as a therapeutic criterion to identify drug targets. However, these methods do not consider the side effects during the identification stage. This study developed a fuzzy multi-objective optimization for identifying anti-cancer targets that not only evaluated cancer cell mortality, but also minimized side effects due to treatment. We identified potential anti-cancer enzymes and antimetabolites for the treatment of head and neck cancer (HNC). The identified one- and two-target enzymes were primarily involved in six major pathways, namely, purine and pyrimidine metabolism and the pentose phosphate pathway. Most of the identified targets can be regulated by approved drugs; thus, these drugs are potential candidates for drug repurposing as a treatment for HNC. Furthermore, we identified antimetabolites involved in pathways similar to those identified using a gene-centric approach. Moreover, HMGCR knockdown could not block the growth of HNC cells. However, the two-target combinations of (UMPS, HMGCR) and (CAD, HMGCR) could achieve cell mortality and improve metabolic deviation grades over 22% without reducing the cell viability grade.

Genome-Scale Metabolic Model Analysis of Metabolic Differences between Lauren Diffuse and Intestinal Subtypes in Gastric Cancer

Gastric cancer (GC) is one of the most lethal cancers worldwide; it has a high mortality rate, particularly in East Asia. Recently, genetic events (e.g., mutations and copy number alterations) and molecular signaling associated with histologically different GC subtypes (diffuse and intestinal) have been elucidated. However, metabolic differences among the histological GC subtypes have not been studied systematically. In this study, we utilized transcriptome-based genome-scale metabolic models (GEMs) to identify differential metabolic pathways between Lauren diffuse and intestinal subtypes. We found that diverse metabolic pathways, including cholesterol homeostasis, xenobiotic metabolism, fatty acid metabolism, the MTORC1 pathway, and glycolysis, were dysregulated between the diffuse and intestinal subtypes. Our study provides an overview of the metabolic differences between the two subtypes, possibly leading to an understanding of metabolism in GC heterogeneity.

Lipid exposure activates gene expression changes associated with estrogen receptor negative breast cancer

Improved understanding of local breast biology that favors the development of estrogen receptor negative (ER-) breast cancer (BC) would foster better prevention strategies. We have previously shown that overexpression of specific lipid metabolism genes is associated with the development of ER- BC. We now report results of exposure of MCF-10A and MCF-12A cells, and mammary organoids to representative medium- and long-chain polyunsaturated fatty acids. This exposure caused a dynamic and profound change in gene expression, accompanied by changes in chromatin packing density, chromatin accessibility, and histone posttranslational modifications (PTMs). We identified 38 metabolic reactions that showed significantly increased activity, including reactions related to one-carbon metabolism. Among these reactions are those that produce S-adenosyl-L-methionine for histone PTMs. Utilizing both an in-vitro model and samples from women at high risk for ER- BC, we show that lipid exposure engenders gene expression, signaling pathway activation, and histone marks associated with the development of ER- BC.

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.

Modeling of Tumor Growth with Input from Patient-Specific Metabolomic Data

Advances in omic technologies have provided insight into cancer progression and treatment response. However, the nonlinear characteristics of cancer growth present a challenge to bridge from the molecular- to the tissue-scale, as tumor behavior cannot be encapsulated by the sum of the individual molecular details gleaned experimentally. Mathematical modeling and computational simulation have been traditionally employed to facilitate analysis of nonlinear systems. In this study, for the first time tumor metabolomic data are linked via mathematical modeling to the tumor tissue-scale behavior, showing the capability to mechanistically simulate cancer progression personalized to omic information obtainable from patient tumor core biopsy analysis. Generally, a higher degree of metabolic dysregulation has been correlated with more aggressive tumor behavior. Accordingly, key parameters influenced by metabolomic data in this model include tumor proliferation, vascularization, aggressiveness, lactic acid production, monocyte infiltration and macrophage polarization, and drug effect. The model enables evaluating interactions of interest between these parameters which drive tumor growth based on the metabolomic data. The results show that the model can group patients consistently with the clinically observed outcomes of response/non-response to chemotherapy. This modeling approach provides a first step towards evaluation of tumor growth based on tumor-specific metabolomic data.

A mechanistic modeling framework reveals the key principles underlying tumor metabolism

While aerobic glycolysis, or the Warburg effect, has for a long time been considered a hallmark of tumor metabolism, recent studies have revealed a far more complex picture. Tumor cells exhibit widespread metabolic heterogeneity, not only in their presentation of the Warburg effect but also in the nutrients and the metabolic pathways they are dependent on. Moreover, tumor cells can switch between different metabolic phenotypes in response to environmental cues and therapeutic interventions. A framework to analyze the observed metabolic heterogeneity and plasticity is, however, lacking. Using a mechanistic model that includes the key metabolic pathways active in tumor cells, we show that the inhibition of phosphofructokinase by excess ATP in the cytoplasm can drive a preference for aerobic glycolysis in fast-proliferating tumor cells. The differing rates of ATP utilization by tumor cells can therefore drive heterogeneity with respect to the presentation of the Warburg effect. Building upon this idea, we couple the metabolic phenotype of tumor cells to their migratory phenotype, and show that our model predictions are in agreement with previous experiments. Next, we report that the reliance of proliferating cells on different anaplerotic pathways depends on the relative availability of glucose and glutamine, and can further drive metabolic heterogeneity. Finally, using treatment of melanoma cells with a BRAF inhibitor as an example, we show that our model can be used to predict the metabolic and gene expression changes in cancer cells in response to drug treatment. By making predictions that are far more generalizable and interpretable as compared to previous tumor metabolism modeling approaches, our framework identifies key principles that govern tumor cell metabolism, and the reported heterogeneity and plasticity. These principles could be key to targeting the metabolic vulnerabilities of cancer.

Machine learning-guided evaluation of extraction and simulation methods for cancer patient-specific metabolic models

Genome-scale metabolic model (GEM) has been established as an important tool to study cellular metabolism at a systems level by predicting intracellular fluxes. With the advent of generic human GEMs, they have been increasingly applied to a range of diseases, often for the objective of predicting effective metabolic drug targets. Cancer is a representative disease where the use of GEMs has proved to be effective, partly due to the massive availability of patient-specific RNA-seq data. When using a human GEM, so-called context-specific GEM needs to be developed first by using cell-specific RNA-seq data. Biological validity of a context-specific GEM highly depends on both model extraction method (MEM) and model simulation method (MSM). However, while MEMs have been thoroughly examined, MSMs have not been systematically examined, especially, when studying cancer metabolism. In this study, the effects of pairwise combinations of three MEMs and five MSMs were evaluated by examining biological features of the resulting cancer patient-specific GEMs. For this, a total of 1,562 patient-specific GEMs were reconstructed, and subjected to machine learning-guided and biological evaluations to draw robust conclusions. Noteworthy observations were made from the evaluation, including the high performance of two MEMs, namely rank-based ‘task-driven Integrative Network Inference for Tissues’ (tINIT) or ‘Gene Inactivity Moderated by Metabolism and Expression’ (GIMME), paired with least absolute deviation (LAD) as a MSM, and relatively poorer performance of flux balance analysis (FBA) and parsimonious FBA (pFBA). Insights from this study can be considered as a reference when studying cancer metabolism using patient-specific GEMs.

Development of computational models using omics data for the identification of effective cancer metabolic biomarkers

Identification of novel biomarkers has been an active area of study for the effective diagnosis, prognosis and treatment of cancers. Among various types of cancer biomarkers, metabolic biomarkers, including enzymes, metabolites and metabolic genes, deserve attention as they can serve as a reliable source for diagnosis, prognosis and treatment of cancers. In particular, efforts to identify novel biomarkers have been greatly facilitated by a rapid increase in the volume of multiple omics data generated for a range of cancer cells. These omics data in turn serve as ingredients for developing computational models that can help derive deeper insights into the biology of cancer cells, and identify metabolic biomarkers. In this review, we provide an overview of omics data generated for cancer cells, and discuss recent studies on computational models that were developed using omics data in order to identify effective cancer metabolic biomarkers.

DCcov: Repositioning of drugs and drug combinations for SARS-CoV-2 infected lung through constraint-based modeling

The 2019 coronavirus disease (COVID-19) became a worldwide pandemic with currently no approved effective antiviral drug. Flux balance analysis (FBA) is an efficient method to analyze metabolic networks. Here, FBA was applied on human lung cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to reposition metabolic drugs and drug combinations against the virus replication within the host tissue. Making use of expression datasets of infected lung tissue, genome-scale COVID-19-specific metabolic models were reconstructed. Then, host-specific essential genes and gene pairs were determined through in silico knockouts that permit reducing the viral biomass production without affecting the host biomass. Key pathways that are associated with COVID-19 severity in lung tissue are related to oxidative stress, ferroptosis, and pyrimidine metabolism. By in silico screening of Food and Drug Administration (FDA)-approved drugs on the putative disease-specific essential genes and gene pairs, 85 drugs and 52 drug combinations were predicted as promising candidates for COVID-19 (https://github.com/sysbiolux/DCcov).

Comparison of metabolic states using genome-scale metabolic models

Genome-scale metabolic models (GEMs) are comprehensive knowledge bases of cellular metabolism and serve as mathematical tools for studying biological phenotypes and metabolic states or conditions in various organisms and cell types. Given the sheer size and complexity of human metabolism, selecting parameters for existing analysis methods such as metabolic objective functions and model constraints is not straightforward in human GEMs. In particular, comparing several conditions in large GEMs to identify condition- or disease-specific metabolic features is challenging. In this study, we showcase a scalable, model-driven approach for an in-depth investigation and comparison of metabolic states in large GEMs which enables identifying the underlying functional differences. Using a combination of flux space sampling and network analysis, our approach enables extraction and visualisation of metabolically distinct network modules. Importantly, it does not rely on known or assumed objective functions. We apply this novel approach to extract the biochemical differences in adipocytes arising due to unlimited vs blocked uptake of branched-chain amino acids (BCAAs, considered as biomarkers in obesity) using a human adipocyte GEM (iAdipocytes1809). The biological significance of our approach is corroborated by literature reports confirming our identified metabolic processes (TCA cycle and Fatty acid metabolism) to be functionally related to BCAA metabolism. Additionally, our analysis predicts a specific altered uptake and secretion profile indicating a compensation for the unavailability of BCAAs. Taken together, our approach facilitates determining functional differences between any metabolic conditions of interest by offering a versatile platform for analysing and comparing flux spaces of large metabolic networks.

Cellular metabolism is a highly complex and interconnected system. As many lifestyle diseases in humans have a strong metabolic component, it is important to understand metabolic differences between healthy and diseased states. In systems biology, metabolic behaviours are investigated using genome-scale metabolic models. In addition to the sheer size and complexity of the genome-scale metabolic models of human systems, using existing analysis methods is challenging and the parameter selection is not straightforward. Therefore, novel methodological frameworks are necessary for analysing metabolic conditions despite the challenges posed by human models. Particularly, an ongoing challenge has been that of comparing several phenotypes for identifying condition- or disease-specific metabolic signatures. We address this significant challenge by developing a scalable and model-driven approach, ComMet (Comparison of Metabolic states). ComMet enables an in-depth investigation and comparison of metabolic phenotypes in large models while also identifying the underlying functional differences. Novel hypotheses can be generated using ComMet for not only understanding known metabolic phenotypes better but also for guiding the design of new experiments to validate the processes predicted by ComMet.

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.

Genome Scale Modeling to Study the Metabolic Competition between Cells in the Tumor Microenvironment

The tumor's physiology emerges from the dynamic interplay of numerous cell types, such as cancer cells, immune cells and stromal cells, within the tumor microenvironment. Immune and cancer cells compete for nutrients within the tumor microenvironment, leading to a metabolic battle between these cell populations. Tumor cells can reprogram their metabolism to meet the high demand of building blocks and ATP for proliferation, and to gain an advantage over the action of immune cells. The study of the metabolic reprogramming mechanisms underlying cancer requires the quantification of metabolic fluxes which can be estimated at the genome-scale with constraint-based or kinetic modeling. Constraint-based models use a set of linear constraints to simulate steady-state metabolic fluxes, whereas kinetic models can simulate both the transient behavior and steady-state values of cellular fluxes and concentrations. The integration of cell- or tissue-specific data enables the construction of context-specific models that reflect cell-type- or tissue-specific metabolic properties. While the available modeling frameworks enable limited modeling of the metabolic crosstalk between tumor and immune cells in the tumor stroma, future developments will likely involve new hybrid kinetic/stoichiometric formulations.

Thermodynamic Genome-Scale Metabolic Modeling of Metallodrug Resistance in Colorectal Cancer

Simple Summary

Cancer, but also its treatment, can lead to a reprogramming of cellular metabolism. These changes are observable in metabolite abundances, which can be unbiasedly measured via mass spectrometry metabolomics. However, even when the metabolome changes strongly, a (mechanistic) interpretation is difficult as metabolite levels do not necessarily directly correspond to pathway activities. Here we measure the changes of the cellular metabolome in colorectal cancer cell lines sensitive and resistant to the ruthenium-based drug BOLD-100/KP1339 and the platinum-based drug oxaliplatin. We map these changes onto a cancer-specific genome-scale metabolic model, which allows us not only to compute intracellular flux distributions, but also to disentangle drug-specific effects from growth differences from differences in metabolic adaptations due to resistance. Specifically, we find that resistance to BOLD-100/KP1339 induces more extensive reprogramming than oxaliplatin, especially with respect to fatty acid and amino acid metabolism.

Abstract

Background: Mass spectrometry-based metabolomics approaches provide an immense opportunity to enhance our understanding of the mechanisms that underpin the cellular reprogramming of cancers. Accurate comparative metabolic profiling of heterogeneous conditions, however, is still a challenge. Methods: Measuring both intracellular and extracellular metabolite concentrations, we constrain four instances of a thermodynamic genome-scale metabolic model of the HCT116 colorectal carcinoma cell line to compare the metabolic flux profiles of cells that are either sensitive or resistant to ruthenium- or platinum-based treatments with BOLD-100/KP1339 and oxaliplatin, respectively. Results: Normalizing according to growth rate and normalizing resistant cells according to their respective sensitive controls, we are able to dissect metabolic responses specific to the drug and to the resistance states. We find the normalization steps to be crucial in the interpretation of the metabolomics data and show that the metabolic reprogramming in resistant cells is limited to a select number of pathways. Conclusions: Here, we elucidate the key importance of normalization steps in the interpretation of metabolomics data, allowing us to uncover drug-specific metabolic reprogramming during acquired metal-drug resistance.

Stoichiometric Modeling of Artificial String Chemistries Reveals Constraints on Metabolic Network Structure

Uncovering the general principles that govern the structure of metabolic networks is key to understanding the emergence and evolution of living systems. Artificial chemistries can help illuminate this problem by enabling the exploration of chemical reaction universes that are constrained by general mathematical rules. Here, we focus on artificial chemistries in which strings of characters represent simplified molecules, and string concatenation and splitting represent possible chemical reactions. We developed a novel Python package, ARtificial CHemistry NEtwork Toolbox (ARCHNET), to study string chemistries using tools from the field of stoichiometric constraint-based modeling. In addition to exploring the topological characteristics of different string chemistry networks, we developed a network-pruning algorithm that can generate minimal metabolic networks capable of producing a specified set of biomass precursors from a given assortment of environmental nutrients. We found that the composition of these minimal metabolic networks was influenced more strongly by the metabolites in the biomass reaction than the identities of the environmental nutrients. This finding has important implications for the reconstruction of organismal metabolic networks and could help us better understand the rise and evolution of biochemical organization. More generally, our work provides a bridge between artificial chemistries and stoichiometric modeling, which can help address a broad range of open questions, from the spontaneous emergence of an organized metabolism to the structure of microbial communities.

Genome-scale metabolic network reconstruction of model animals as a platform for translational research

Genome-scale metabolic models (GEMs) are used extensively for analysis of mechanisms underlying human diseases and metabolic malfunctions. However, the lack of comprehensive and high-quality GEMs for model organisms restricts translational utilization of omics data accumulating from the use of various disease models. Here we present a unified platform of GEMs that covers five major model animals, including Mouse1 (Mus musculus), Rat1 (Rattus norvegicus), Zebrafish1 (Danio rerio), Fruitfly1 (Drosophila melanogaster), and Worm1 (Caenorhabditis elegans). These GEMs represent the most comprehensive coverage of the metabolic network by considering both orthology-based pathways and species-specific reactions. All GEMs can be interactively queried via the accompanying web portal Metabolic Atlas. Specifically, through integrative analysis of Mouse1 with RNA-sequencing data from brain tissues of transgenic mice we identified a coordinated up-regulation of lysosomal GM2 ganglioside and peptide degradation pathways which appears to be a signature metabolic alteration in Alzheimer's disease (AD) mouse models with a phenotype of amyloid precursor protein overexpression. This metabolic shift was further validated with proteomics data from transgenic mice and cerebrospinal fluid samples from human patients. The elevated lysosomal enzymes thus hold potential to be used as a biomarker for early diagnosis of AD. Taken together, we foresee that this evolving open-source platform will serve as an important resource to facilitate the development of systems medicines and translational biomedical applications.

Optimization of a modeling platform to predict oncogenes from genome-scale metabolic networks of non-small-cell lung cancers

Cancer cell dysregulations result in the abnormal regulation of cellular metabolic pathways. By simulating this metabolic reprogramming using constraint-based modeling approaches, oncogenes can be predicted, and this knowledge can be used in prognosis and treatment. We introduced a trilevel optimization problem describing metabolic reprogramming for inferring oncogenes. First, this study used RNA-Seq expression data of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) samples and their healthy counterparts to reconstruct tissue-specific genome-scale metabolic models and subsequently build the flux distribution pattern that provided a measure for the oncogene inference optimization problem for determining tumorigenesis. The platform detected 45 genes for LUAD and 84 genes for LUSC that lead to tumorigenesis. A high level of differentially expressed genes was not an essential factor for determining tumorigenesis. The platform indicated that pyruvate kinase (PKM), a well-known oncogene with a low level of differential gene expression in LUAD and LUSC, had the highest fitness among the predicted oncogenes based on computation. By contrast, pyruvate kinase L/R (PKLR), an isozyme of PKM, had a high level of differential gene expression in both cancers. Phosphatidylserine synthase 1 (PTDSS1), an oncogene in LUAD, was inferred to have a low level of differential gene expression, and overexpression could significantly reduce survival probability. According to the factor analysis, PTDSS1 characteristics were close to those of the template, but they were unobvious in LUSC. Angiotensin-converting enzyme 2 (ACE2) has recently garnered widespread interest as the SARS-CoV-2 virus receptor. Moreover, we determined that ACE2 is an oncogene of LUSC but not of LUAD. The platform developed in this study can identify oncogenes with low levels of differential expression and be used to identify potential therapeutic targets for cancer treatment.

Breast cancer dormancy: need for clinically relevant models to address current gaps in knowledge

Breast cancer is the most commonly diagnosed cancer in the USA. Although advances in treatment over the past several decades have significantly improved the outlook for this disease, most women who are diagnosed with estrogen receptor positive disease remain at risk of metastatic relapse for the remainder of their life. The cellular source of late relapse in these patients is thought to be disseminated tumor cells that reactivate after a long period of dormancy. The biology of these dormant cells and their natural history over a patient's lifetime is largely unclear. We posit that research on tumor dormancy has been significantly limited by the lack of clinically relevant models. This review will discuss existing dormancy models, gaps in biological understanding, and propose criteria for future models to enhance their clinical relevance.

Systematic evaluation of parameters for genome-scale metabolic models of cultured mammalian cells

Genome-scale metabolic models describe cellular metabolism with mechanistic detail. Given their high complexity, such models need to be parameterized correctly to yield accurate predictions and avoid overfitting. Effective parameterization has been well-studied for microbial models, but it remains unclear for higher eukaryotes, including mammalian cells. To address this, we enumerated model parameters that describe key features of cultured mammalian cells - including cellular composition, bioprocess performance metrics, mammalian-specific pathways, and biological assumptions behind model formulation approaches. We tested these parameters by building thousands of metabolic models and evaluating their ability to predict the growth rates of a panel of phenotypically diverse Chinese Hamster Ovary cell clones. We found the following considerations to be most critical for accurate parameterization: (1) cells limit metabolic activity to maintain homeostasis, (2) cell morphology and viability change dynamically during a growth curve, and (3) cellular biomass has a particular macromolecular composition. Depending on parameterization, models predicted different metabolic phenotypes, including contrasting mechanisms of nutrient utilization and energy generation, leading to varying accuracies of growth rate predictions. Notably, accurate parameter values broadly agreed with experimental measurements. These insights will guide future investigations of mammalian metabolism.

Current State of "Omics" Biomarkers in Pancreatic Cancer

Pancreatic cancer is one of the most fatal malignancies and the seventh leading cause of cancer-related deaths related to late diagnosis, poor survival rates, and high incidence of metastasis. Unfortunately, pancreatic cancer is predicted to become the third leading cause of cancer deaths in the future. Therefore, diagnosis at the early stages of pancreatic cancer for initial diagnosis or postoperative recurrence is a great challenge, as well as predicting prognosis precisely in the context of biomarker discovery. From the personalized medicine perspective, the lack of molecular biomarkers for patient selection confines tailored therapy options, including selecting drugs and their doses or even diet. Currently, there is no standardized pancreatic cancer screening strategy using molecular biomarkers, but CA19-9 is the most well known marker for the detection of pancreatic cancer. In contrast, recent innovations in high-throughput techniques have enabled the discovery of specific biomarkers of cancers using genomics, transcriptomics, proteomics, metabolomics, glycomics, and metagenomics. Panels combining CA19-9 with other novel biomarkers from different "omics" levels might represent an ideal strategy for the early detection of pancreatic cancer. The systems biology approach may shed a light on biomarker identification of pancreatic cancer by integrating multi-omics approaches. In this review, we provide background information on the current state of pancreatic cancer biomarkers from multi-omics stages. Furthermore, we conclude this review on how multi-omics data may reveal new biomarkers to be used for personalized medicine in the future.

Cysteine and Folate Metabolism Are Targetable Vulnerabilities of Metastatic Colorectal Cancer

Simple Summary

In this work, we studied the metabolic reprogramming of same-patient-derived cell lines with increasing metastatic potential to develop new therapeutic approaches against metastatic colorectal cancer. Using a novel systems biology approach to integrate multiple layers of omics data, we predicted and validated that cystine uptake and folate metabolism, two key pathways related to redox metabolism, are potential targets against metastatic colorectal cancer. Our findings indicate that metastatic cell lines are selectively dependent on redox homeostasis, paving the way for new targeted therapies.

Abstract

With most cancer-related deaths resulting from metastasis, the development of new therapeutic approaches against metastatic colorectal cancer (mCRC) is essential to increasing patient survival. The metabolic adaptations that support mCRC remain undefined and their elucidation is crucial to identify potential therapeutic targets. Here, we employed a strategy for the rational identification of targetable metabolic vulnerabilities. This strategy involved first a thorough metabolic characterisation of same-patient-derived cell lines from primary colon adenocarcinoma (SW480), its lymph node metastasis (SW620) and a liver metastatic derivative (SW620-LiM2), and second, using a novel multi-omics integration workflow, identification of metabolic vulnerabilities specific to the metastatic cell lines. We discovered that the metastatic cell lines are selectively vulnerable to the inhibition of cystine import and folate metabolism, two key pathways in redox homeostasis. Specifically, we identified the system xCT and MTHFD1 genes as potential therapeutic targets, both individually and combined, for combating mCRC.

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.

Identifying Personalized Metabolic Signatures in Breast Cancer

Cancer cells are adept at reprogramming energy metabolism, and the precise manifestation of this metabolic reprogramming exhibits heterogeneity across individuals (and from cell to cell). In this study, we analyzed the metabolic differences between interpersonal heterogeneous cancer phenotypes. We used divergence analysis on gene expression data of 1156 breast normal and tumor samples from The Cancer Genome Atlas (TCGA) and integrated this information with a genome-scale reconstruction of human metabolism to generate personalized, context-specific metabolic networks. Using this approach, we classified the samples into four distinct groups based on their metabolic profiles. Enrichment analysis of the subsystems indicated that amino acid metabolism, fatty acid oxidation, citric acid cycle, androgen and estrogen metabolism, and reactive oxygen species (ROS) detoxification distinguished these four groups. Additionally, we developed a workflow to identify potential drugs that can selectively target genes associated with the reactions of interest. MG-132 (a proteasome inhibitor) and OSU-03012 (a celecoxib derivative) were the top-ranking drugs identified from our analysis and known to have anti-tumor activity. Our approach has the potential to provide mechanistic insights into cancer-specific metabolic dependencies, ultimately enabling the identification of potential drug targets for each patient independently, contributing to a rational personalized medicine approach.

Metabolic radiogenomics in lung cancer: associations between FDG PET image features and oncogenic signaling pathway alterations

This study investigated the associations between image features extracted from tumor ¹⁸F-fluorodeoxyglucose (FDG) uptake and genetic alterations in patients with lung cancer. A total of 137 patients (age, 62.7 ± 10.2 years) who underwent FDG positron emission tomography/computed tomography (PET/CT) and targeted deep sequencing analysis for a tumor lesion, comprising 61 adenocarcinoma (ADC), 31 squamous cell carcinoma (SQCC), and 45 small cell lung cancer (SCLC) patients, were enrolled in this study. From the tumor lesions, 86 image features were extracted, and 381 genes were assessed. PET features were associated with genetic mutations: 41 genes with 24 features in ADC; 35 genes with 22 features in SQCC; and 43 genes with 25 features in SCLC (FDR < 0.05). Clusters based on PET features showed an association with alterations in oncogenic signaling pathways: Cell cycle and WNT signaling pathways in ADC (p = 0.023, p = 0.035, respectively); Cell cycle, p53, and WNT in SQCC (p = 0.045, 0.009, and 0.029, respectively); and TGFβ in SCLC (p = 0.030). In addition, SUV_peak and SUV_max were associated with a mutation of the TGFβ signaling pathway in ADC (FDR = 0.001, < 0.001). In this study, PET image features had significant associations with alterations in genes and oncogenic signaling pathways in patients with lung cancer.

Pyrvinium Pamoate Induces Death of Triple-Negative Breast Cancer Stem-Like Cells and Reduces Metastases through Effects on Lipid Anabolism

Cancer stem-like cells (CSC) induce aggressive tumor phenotypes such as metastasis formation, which is associated with poor prognosis in triple-negative breast cancer (TNBC). Repurposing of FDA-approved drugs that can eradicate the CSC subcompartment in primary tumors may prevent metastatic disease, thus representing an effective strategy to improve the prognosis of TNBC. Here, we investigated spheroid-forming cells in a metastatic TNBC model. This strategy enabled us to specifically study a population of long-lived tumor cells enriched in CSCs, which show stem-like characteristics and induce metastases. To repurpose FDA-approved drugs potentially toxic for CSCs, we focused on pyrvinium pamoate (PP), an anthelmintic drug with documented anticancer activity in preclinical models. PP induced cytotoxic effects in CSCs and prevented metastasis formation. Mechanistically, the cell killing effects of PP were a result of inhibition of lipid anabolism and, more specifically, the impairment of anabolic flux from glucose to cholesterol and fatty acids. CSCs were strongly dependent upon activation of lipid biosynthetic pathways; activation of these pathways exhibited an unfavorable prognostic value in a cohort of breast cancer patients, where it predicted high probability of metastatic dissemination and tumor relapse. Overall, this work describes a new approach to target aggressive CSCs that may substantially improve clinical outcomes for patients with TNBC, who currently lack effective targeted therapeutic options. SIGNIFICANCE: These findings provide preclinical evidence that a drug repurposing approach to prevent metastatic disease in TNBC exploits lipid anabolism as a metabolic vulnerability against CSCs in primary tumors.

Nutrient Sensing by Histone Marks: Reading the Metabolic Histone Code Using Tracing, Omics, and Modeling

Several metabolites serve as substrates for histone modifications and communicate changes in the metabolic environment to the epigenome. Technologies such as metabolomics and proteomics have allowed us to reconstruct the interactions between metabolic pathways and histones. These technologies have shed light on how nutrient availability can have a dramatic effect on various histone modifications. This metabolism-epigenome cross talk plays a fundamental role in development, immune function, and diseases like cancer. Yet, major challenges remain in understanding the interactions between cellular metabolism and the epigenome. How the levels and fluxes of various metabolites impact epigenetic marks is still unclear. Discussed herein are recent applications and the potential of systems biology methods such as flux tracing and metabolic modeling to address these challenges and to uncover new metabolic-epigenetic interactions. These systems approaches can ultimately help elucidate how nutrients shape the epigenome of microbes and mammalian cells.

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.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

Computation-Driven Analysis of Model Polyextremo-tolerant Fungus Exophiala dermatitidis: Defensive Pigment Metabolic Costs and Human Applications

The polyextremotolerant black yeast Exophiala dermatitidis is a tractable model system for investigation of adaptations that support growth under extreme conditions. Foremost among these adaptations are melanogenesis and carotenogenesis. A particularly important question is their metabolic production cost. However, investigation of this issue has been hindered by a relatively poor systems-level understanding of E. dermatitidis metabolism. To address this challenge, a genome-scale model (iEde2091) was developed. Using iEde2091, carotenoids were found to be more expensive to produce than melanins. Given their overlapping protective functions, this suggests that carotenoids have an underexplored yet important role in photo-protection. Furthermore, multiple defensive pigments with overlapping functions might allow E. dermatitidis to minimize cost. Because iEde2091 revealed that E. dermatitidis synthesizes the same melanins as humans and the active sites of the key tyrosinase enzyme are highly conserved this model may enable a broader understanding of melanin production across kingdoms.

Uncovering the Underlying Mechanisms of Cancer Metabolism through the Landscapes and Probability Flux Quantifications

Cancer metabolism is critical for understanding the mechanism of tumorigenesis, yet the understanding is still challenging. We studied gene-metabolism regulatory interactions and quantified the global driving forces for cancer-metabolism dynamics as the underlying landscape and probability flux. We uncovered four steady-state attractors: a normal state attractor, a cancer OXPHOS state attractor, a cancer glycolysis state attractor, and an intermediate cancer state attractor. We identified the key regulatory interactions through global sensitivity analysis based on the landscape topography. Different landscape topographies of glycolysis switch between normal cells and cancer cells were identified. We uncovered that the normal state to cancer state transformation is associated with the peaks of the probability flux and the thermodynamic dissipation, giving dynamical and thermodynamic origin of cancer formation. We found that cancer metabolism oscillations consume more energy to support cancer malignancy. This study provides a quantitative understanding of cancer metabolism and suggests a metabolic therapeutic strategy.

MaREA4Galaxy: Metabolic reaction enrichment analysis and visualization of RNA-seq data within Galaxy

We present MaREA4Galaxy, a user-friendly tool that allows a user to characterize and to graphically compare groups of samples with different transcriptional regulation of metabolism, as estimated from cross-sectional RNA-seq data. The tool is available as plug-in for the widely-used Galaxy platform for comparative genomics and bioinformatics analyses. MaREA4Galaxy combines three modules. The Expression2RAS module, which, for each reaction of a specified set, computes a Reaction Activity Score (RAS) as a function of the expression level of genes encoding for the associated enzyme. The MaREA (Metabolic Reaction Enrichment Analysis) module that allows to highlight significant differences in reaction activities between specified groups of samples. The Clustering module which employs the RAS computed before as a metric for unsupervised clustering of samples into distinct metabolic subgroups; the Clustering tool provides different clustering techniques and implements standard methods to evaluate the goodness of the results.