Systems biology analysis of drivers underlying hallmarks of cancer cell metabolism

Direct Infusion Mass Spectrometry to Rapidly Map Metabolic Flux of Substrates Labeled with Stable Isotopes

Direct infusion-high-resolution mass spectrometry (DI-HRMS) allows for rapid profiling of complex mixtures of metabolites in blood, cerebrospinal fluid, tissue samples and cultured cells. Here, we present a DI-HRMS method suitable for the rapid determination of metabolic fluxes of isotopically labeled substrates in cultured cells and organoids. We adapted an automated annotation pipeline by selecting labeled adducts that best represent the majority of 13C and/or 15N-labeled glycolytic and tricarboxylic acid cycle intermediates as well as a number of their derivatives. Furthermore, valine, leucine and several of their degradation products were included. We show that DI-HRMS can determine anticipated and unanticipated alterations in metabolic fluxes along these pathways that result from the genetic alteration of single metabolic enzymes, including pyruvate dehydrogenase (PDHA1) and glutaminase (GLS). In addition, it can precisely pinpoint metabolic adaptations to the loss of methylmalonyl-CoA mutase in patient-derived liver organoids. Our results highlight the power of DI-HRMS in combination with stable isotopically labeled compounds as an efficient screening method for fluxomics.

Aged breast matrix bound vesicles promote breast cancer invasiveness

Aging is one of the inherent risk factors for breast cancer. Although the influence of age-related cellular alterations on breast cancer development has been extensively explored, little is known about the alterations in the aging breast tissue microenvironment, specifically the extracellular matrix (ECM). Here, for the first time in literature, we have identified tissue resident matrix bound vesicles (MBVs) within the healthy mouse breast ECM, investigated and compared their characteristics in young and aged healthy breast tissues, and studied the effects of these MBVs on normal (KTB21) and cancerous (MDA-MB-231) human mammary epithelial cells with respect to the tissue age that they are extracted from. Using vesicle labeling technology, we were able to visualize cellular uptake of the MBVs directly from the native decellularized tissue sections, showing that these MBVs have regulatory roles in the tissue microenvironment. We mimicked the ECM by embedding the MBVs in collagen gels, and showed that MBVs could be taken up by the cells. The miRNA and cytokine profiling showed that MBVs shifted towards a more tumorigenic and invasive phenotype with age, as evidenced by the more pronounced presence of cancer-associated cytokines, and higher expression levels of oncomiRs miR-10b, miR-30e, and miR-210 in MBVs isolated from aged mice. When treated with MBVs or these upregulated factors, KTB21 and MDA-MB-231 cells showed significantly higher motility and invasion compared to untreated controls. Treatment of cells with a cocktail of miRNAs (miR-10b, miR-30e, and miR-210) or with the agonist of adiponectin (AdipoRon), which both were enriched in the aged MBVs, recapitulated the effect of aged MBVs on cells. This study shows for the first time that the MBVs have a regulatory role in the tissue microenvironment and that the MBV contents change towards cancer-promoting upon aging. Studying the effects of MBVs and their cargos on cellular behavior could lead to a better understanding of the critical roles of MBVs played in breast cancer progression and metastasis.

Mitochondrial ATP generation is more proteome efficient than glycolysis

Metabolic efficiency profoundly influences organismal fitness. Nonphotosynthetic organisms, from yeast to mammals, derive usable energy primarily through glycolysis and respiration. Although respiration is more energy efficient, some cells favor glycolysis even when oxygen is available (aerobic glycolysis, Warburg effect). A leading explanation is that glycolysis is more efficient in terms of ATP production per unit mass of protein (that is, faster). Through quantitative flux analysis and proteomics, we find, however, that mitochondrial respiration is actually more proteome efficient than aerobic glycolysis. This is shown across yeast strains, T cells, cancer cells, and tissues and tumors in vivo. Instead of aerobic glycolysis being valuable for fast ATP production, it correlates with high glycolytic protein expression, which promotes hypoxic growth. Aerobic glycolytic yeasts do not excel at aerobic growth but outgrow respiratory cells during oxygen limitation. We accordingly propose that aerobic glycolysis emerges from cells maintaining a proteome conducive to both aerobic and hypoxic growth.

Reconstruction of cell-specific models capturing the influence of metabolism on DNA methylation in cancer

The imbalance of epigenetic regulatory mechanisms such as DNA methylation, which can promote aberrant gene expression profiles without affecting the DNA sequence, may cause the deregulation of signaling, regulatory, and metabolic processes, contributing to a cancerous phenotype. Since some metabolites are substrates and cofactors of epigenetic regulators, their availability can be affected by characteristic cancer cell metabolic shifts, feeding cancer onset and progression through epigenetic deregulation. Hence, there is a need to study the influence of cancer metabolic reprogramming in DNA methylation to design new effective treatments. In this study, a generic Genome-Scale Metabolic Model (GSMM) of a human cell, integrating DNA methylation or demethylation reactions, was obtained and used for the reconstruction of Genome-Scale Metabolic Models enhanced with Enzymatic Constraints using Kinetic and Omics data (GECKOs) of 31 cancer cell lines. Furthermore, cell-line-specific DNA methylation levels were included in the models, as coefficients of a DNA composition pseudo-reaction, to depict the influence of metabolism over global DNA methylation in each of the cancer cell lines. Flux simulations demonstrated the ability of these models to provide simulated fluxes of exchange reactions similar to the equivalent experimentally measured uptake/secretion rates and to make good functional predictions. In addition, simulations found metabolic pathways, reactions and enzymes directly or inversely associated with the gene promoter methylation. Two potential candidates for targeted cancer epigenetic therapy were identified.

Endothelial cells metabolically regulate breast cancer invasion toward a microvessel

Breast cancer metastasis is initiated by invasion of tumor cells into the collagen type I-rich stroma to reach adjacent blood vessels. Prior work has identified that metabolic plasticity is a key requirement of tumor cell invasion into collagen. However, it remains largely unclear how blood vessels affect this relationship. Here, we developed a microfluidic platform to analyze how tumor cells invade collagen in the presence and absence of a microvascular channel. We demonstrate that endothelial cells secrete pro-migratory factors that direct tumor cell invasion toward the microvessel. Analysis of tumor cell metabolism using metabolic imaging, metabolomics, and computational flux balance analysis revealed that these changes are accompanied by increased rates of glycolysis and oxygen consumption caused by broad alterations of glucose metabolism. Indeed, restricting glucose availability decreased endothelial cell-induced tumor cell invasion. Our results suggest that endothelial cells promote tumor invasion into the stroma due, in part, to reprogramming tumor cell metabolism.

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.

Systems Biology and Cytokines Potential Role in Lung Cancer Immunotherapy Targeting Autophagic Axis

Lung cancer accounts for the highest number of deaths among men and women worldwide. Although extensive therapies, either alone or in conjunction with some specific drugs, continue to be the principal regimen for evolving lung cancer, significant improvements are still needed to understand the inherent biology behind progressive inflammation and its detection. Unfortunately, despite every advancement in its treatment, lung cancer patients display different growth mechanisms and continue to die at significant rates. Autophagy, which is a physiological defense mechanism, serves to meet the energy demands of nutrient-deprived cancer cells and sustain the tumor cells under stressed conditions. In contrast, autophagy is believed to play a dual role during different stages of tumorigenesis. During early stages, it acts as a tumor suppressor, degrading oncogenic proteins; however, during later stages, autophagy supports tumor cell survival by minimizing stress in the tumor microenvironment. The pivotal role of the IL6-IL17-IL23 signaling axis has been observed to trigger autophagic events in lung cancer patients. Since the obvious roles of autophagy are a result of different immune signaling cascades, systems biology can be an effective tool to understand these interconnections and enhance cancer treatment and immunotherapy. In this review, we focus on how systems biology can be exploited to target autophagic processes that resolve inflammatory responses and contribute to better treatment in carcinogenesis.

TFAP2A activates HMGA1 to promote glycolysis and lung adenocarcinoma progression

BACKGROUND

Lung cancer is the most common cancer in the world. High Mobility Group AT-Hook 1 (HMGA1) is found to be associated with the glycolytic pathway in a variety of cancers, and abnormal glycolysis function is one of the important characteristics of cancer cells. Therefore, this paper discusses the effect of HMGA1 on glycolysis of lung adenocarcinoma (LUAD) cells METHODS: The mRNA expression data were downloaded from TCGA-LUAD database. Groups were set according to the median expression of HMGA1, followed by GSEA enrichment analysis. The upstream transcriptional regulators of HMGA1 were predicted by bioinformatics. The correlation between HMGA1 and Transcription Factor AP-2 Alpha (TFAP2A) and their expression in LUAD tissues were analyzed as well. mRNA expression levels of HMGA1 and TFAP2A were detected by qRT-PCR. The binding of HMGA1 and TFAP2A was demonstrated by ChIP and dual luciferase reporter assays. Cell function experiments were utilized to assay proliferation, apoptosis, glycolysis ability of LUAD cells, and glycolysis-related protein expression in each treatment group.

RESULTS

HMGA1 was highly expressed in LUAD patients' tissues and enriched in the glycolytic pathway. Additionally, silencing HMGA1 markedly hampered cell proliferation and glycolysis, and promoted cell apoptosis. The upstream transcriptional regulator TFAP2A was predicted to be highly expressed in LUAD. ChIP and dual luciferase reporter assays confirmed the targeted relationship between HMGA1 and TFAP2A. Cell rescue assay confirmed that TFAP2A promoted glycolysis and LUAD progression by activating HMGA1.

CONCLUSION

TFAP2A promotes glycolysis, proliferation and hampers apoptosis of LUAD cells by stimulating HMGA1. Hence, TFAP2A/HMGA1 may be a feasible therapeutic target for LUAD.

AVAILABILITY OF DATA AND MATERIALS

All the data within this manuscript could be gotten from corresponding author at reasonable request.

Uses of Multi-Objective Flux Analysis for Optimization of Microbial Production of Secondary Metabolites

Secondary metabolites are not essential for the growth of microorganisms, but they play a critical role in how microbes interact with their surroundings. In addition to this important ecological role, secondary metabolites also have a variety of agricultural, medicinal, and industrial uses, and thus the examination of secondary metabolism of plants and microbes is a growing scientific field. While the chemical production of certain secondary metabolites is possible, industrial-scale microbial production is a green and economically attractive alternative. This is even more true, given the advances in bioengineering that allow us to alter the workings of microbes in order to increase their production of compounds of interest. This type of engineering requires detailed knowledge of the "chassis" organism's metabolism. Since the resources and the catalytic capacity of enzymes in microbes is finite, it is important to examine the tradeoffs between various bioprocesses in an engineered system and alter its working in a manner that minimally perturbs the robustness of the system while allowing for the maximum production of a product of interest. The in silico multi-objective analysis of metabolism using genome-scale models is an ideal method for such examinations.

Characterizing cancer metabolism from bulk and single-cell RNA-seq data using METAFlux

Cells often alter metabolic strategies under nutrient-deprived conditions to support their survival and growth. Characterizing metabolic reprogramming in the tumor microenvironment (TME) is of emerging importance in cancer research and patient care. However, recent technologies only measure a subset of metabolites and cannot provide in situ measurements. Computational methods such as flux balance analysis (FBA) have been developed to estimate metabolic flux from bulk RNA-seq data and can potentially be extended to single-cell RNA-seq (scRNA-seq) data. However, it is unclear how reliable current methods are, particularly in TME characterization. Here, we present a computational framework METAFlux (METAbolic Flux balance analysis) to infer metabolic fluxes from bulk or single-cell transcriptomic data. Large-scale experiments using cell-lines, the cancer genome atlas (TCGA), and scRNA-seq data obtained from diverse cancer and immunotherapeutic contexts, including CAR-NK cell therapy, have validated METAFlux's capability to characterize metabolic heterogeneity and metabolic interaction amongst cell types.

Convergent Approaches to Delineate the Metabolic Regulation of Tumor Invasion by Hyaluronic Acid Biosynthesis

Metastasis is the leading cause of breast cancer-related deaths and is often driven by invasion and cancer-stem like cells (CSCs). Both the CSC phenotype and invasion are associated with increased hyaluronic acid (HA) production. How these independent observations are connected, and which role metabolism plays in this process, remains unclear due to the lack of convergent approaches integrating engineered model systems, computational tools, and cancer biology. Using microfluidic invasion models, metabolomics, computational flux balance analysis, and bioinformatic analysis of patient data, the functional links between the stem-like, invasive, and metabolic phenotype of breast cancer cells as a function of HA biosynthesis are investigated. These results suggest that CSCs are more invasive than non-CSCs and that broad metabolic changes caused by overproduction of HA play a role in this process. Accordingly, overexpression of hyaluronic acid synthases (HAS) 2 or 3 induces a metabolic phenotype that promotes cancer cell stemness and invasion in vitro and upregulates a transcriptomic signature predictive of increased invasion and worse patient survival. This study suggests that HA overproduction leads to metabolic adaptations to satisfy the energy demands for 3D invasion of breast CSCs highlighting the importance of engineered model systems and multidisciplinary approaches in cancer research.

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.

The functions and molecular mechanisms of Tribbles homolog 3 (TRIB3) implicated in the pathophysiology of cancer

Currently, cancer ranks as the second leading cause of death worldwide, and at the same time, the burden of cancer continues to increase. The underlying molecular pathways involved in the initiation and development of cancer are the subject of considerable research worldwide. Further understanding of these pathways may lead to new cancer treatments. Growing data suggest that Tribble's homolog 3 (TRIB3) is essential in oncogenesis in many types of cancer. The mammalian tribbles family's proteins regulate various cellular and physiological functions, such as the cell cycle, stress response, signal transduction, propagation, development, differentiation, immunity, inflammatory processes, and metabolism. To exert their activities, Tribbles proteins must alter key signaling pathways, including the mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3 kinase (PI3K)/AKT pathways. Recent evidence supports that TRIB3 dysregulation has been linked to various diseases, including tumor development and chemoresistance. It has been speculated that TRIB3 may either promote or inhibit the onset and development of cancer. However, it is still unclear how TRIB3 performs this dual function in cancer. In this review, we present and discuss the most recent data on the role of TRIB3 in cancer pathophysiology and chemoresistance. Furthermore, we describe in detail the molecular mechanism TRIB3 regulates in cancer.

Guidelines for extracting biologically relevant context-specific metabolic models using gene expression data

Genome-scale metabolic models comprehensively describe an organism's metabolism and can be tailored using omics data to model condition-specific physiology. The quality of context-specific models is impacted by (i) choice of algorithm and parameters and (ii) alternate context-specific models that equally explain the -omics data. Here we quantify the influence of alternate optima on microbial and mammalian model extraction using GIMME, iMAT, MBA, and mCADRE. We find that metabolic tasks defining an organism's phenotype must be explicitly and quantitatively protected. The scope of alternate models is strongly influenced by algorithm choice and the topological properties of the parent genome-scale model with fatty acid metabolism and intracellular metabolite transport contributing much to alternate solutions in all models. mCADRE extracted the most reproducible context-specific models and models generated using MBA had the most alternate solutions. There were fewer qualitatively different solutions generated by GIMME in E. coli, but these increased substantially in the mammalian models. Screening ensembles using a receiver operating characteristic plot identified the best-performing models. A comprehensive evaluation of models extracted using combinations of extraction methods and expression thresholds revealed that GIMME generated the best-performing models in E. coli, whereas mCADRE is better suited for complex mammalian models. These findings suggest guidelines for benchmarking -omics integration algorithms and motivate the development of a systematic workflow to enumerate alternate models and extract biologically relevant context-specific models.

Saturation of the mitochondrial NADH shuttles drives aerobic glycolysis in proliferating cells

Proliferating cells exhibit a metabolic phenotype known as "aerobic glycolysis," which is characterized by an elevated rate of glucose fermentation to lactate irrespective of oxygen availability. Although several theories have been proposed, a rationalization for why proliferating cells seemingly waste glucose carbon by excreting it as lactate remains elusive. Using the NCI-60 cell lines, we determined that lactate excretion is strongly correlated with the activity of mitochondrial NADH shuttles, but not proliferation. Quantifying the fluxes of the malate-aspartate shuttle (MAS), the glycerol 3-phosphate shuttle (G3PS), and lactate dehydrogenase under various conditions demonstrated that proliferating cells primarily transform glucose to lactate when glycolysis outpaces the mitochondrial NADH shuttles. Increasing mitochondrial NADH shuttle fluxes decreased glucose fermentation but did not reduce the proliferation rate. Our results reveal that glucose fermentation, a hallmark of cancer, is a secondary consequence of MAS and G3PS saturation rather than a unique metabolic driver of cellular proliferation.

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.

Integrating a dynamic central metabolism model of cancer cells with a hybrid 3D multiscale model for vascular hepatocellular carcinoma growth

We develop here a novel modelling approach with the aim of closing the conceptual gap between tumour-level metabolic processes and the metabolic processes occurring in individual cancer cells. In particular, the metabolism in hepatocellular carcinoma derived cell lines (HEPG2 cells) has been well characterized but implementations of multiscale models integrating this known metabolism have not been previously reported. We therefore extend a previously published multiscale model of vascular tumour growth, and integrate it with an experimentally verified network of central metabolism in HEPG2 cells. This resultant combined model links spatially heterogeneous vascular tumour growth with known metabolic networks within tumour cells and accounts for blood flow, angiogenesis, vascular remodelling and nutrient/growth factor transport within a growing tumour, as well as the movement of, and interactions between normal and cancer cells. Model simulations report for the first time, predictions of spatially resolved time courses of core metabolites in HEPG2 cells. These simulations can be performed at a sufficient scale to incorporate clinically relevant features of different tumour systems using reasonable computational resources. Our results predict larger than expected temporal and spatial heterogeneity in the intracellular concentrations of glucose, oxygen, lactate pyruvate, f16bp and Acetyl-CoA. The integrated multiscale model developed here provides an ideal quantitative framework in which to study the relationship between dosage, timing, and scheduling of anti-neoplastic agents and the physiological effects of tumour metabolism at the cellular level. Such models, therefore, have the potential to inform treatment decisions when drug response is dependent on the metabolic state of individual cancer cells.

Metabolic modeling-based drug repurposing in Glioblastoma

The manifestation of intra- and inter-tumor heterogeneity hinders the development of ubiquitous cancer treatments, thus requiring a tailored therapy for each cancer type. Specifically, the reprogramming of cellular metabolism has been identified as a source of potential drug targets. Drug discovery is a long and resource-demanding process aiming at identifying and testing compounds early in the drug development pipeline. While drug repurposing efforts (i.e., inspecting readily available approved drugs) can be supported by a mechanistic rationale, strategies to further reduce and prioritize the list of potential candidates are still needed to facilitate feasible studies. Although a variety of ‘omics’ data are widely gathered, a standard integration method with modeling approaches is lacking. For instance, flux balance analysis is a metabolic modeling technique that mainly relies on the stoichiometry of the metabolic network. However, exploring the network’s topology typically neglects biologically relevant information. Here we introduce Transcriptomics-Informed Stoichiometric Modelling And Network analysis (TISMAN) in a recombinant innovation manner, allowing identification and validation of genes as targets for drug repurposing using glioblastoma as an exemplar.

Cancer cells depend on environmental lipids for proliferation when electron acceptors are limited

Production of oxidized biomass, which requires regeneration of the cofactor NAD⁺, can be a proliferation bottleneck that is influenced by environmental conditions. However, a comprehensive quantitative understanding of metabolic processes that may be affected by NAD⁺ deficiency is currently missing. Here, we show that de novo lipid biosynthesis can impose a substantial NAD⁺ consumption cost in proliferating cancer cells. When electron acceptors are limited, environmental lipids become crucial for proliferation because NAD⁺ is required to generate precursors for fatty acid biosynthesis. We find that both oxidative and even net reductive pathways for lipogenic citrate synthesis are gated by reactions that depend on NAD⁺ availability. We also show that access to acetate can relieve lipid auxotrophy by bypassing the NAD⁺ consuming reactions. Gene expression analysis demonstrates that lipid biosynthesis strongly anti-correlates with expression of hypoxia markers across tumor types. Overall, our results define a requirement for oxidative metabolism to support biosynthetic reactions and provide a mechanistic explanation for cancer cell dependence on lipid uptake in electron acceptor-limited conditions, such as hypoxia.

Cancer Recurrence and Omics: Metabolic Signatures of Cancer Dormancy Revealed by Transcriptome Mapping of Genome-Scale Networks

A major problem in medicine and oncology is cancer recurrence through the activation of dormant cancer cells. A system scale examination of metabolic dysregulations associated with the cancer dormancy offers promise for the discovery of novel molecular targets for cancer precision medicine, and importantly, for the prevention of cancer recurrence. In this study, we mapped the total mRNA sequencing-based transcriptomic data from dormant cancer cell lines and nondormant cancer controls onto a human genome-scale metabolic network by using a graph-based approach, and two mass balance-based approaches with one based on reaction activity/inactivity and the other one on flux changes. The gene expression datasets were accessed from Gene Expression Omnibus (GSE83142 and GSE114012). This analysis included two diverse cancer types, a liquid and a solid cancer, namely, acute lymphoblastic leukemia and colorectal cancer. For the dormant cancer state, we observed changes in major adenosine triphosphate-producing pathways, including the citric acid cycle, oxidative phosphorylation, and glycolysis/gluconeogenesis, indicating a reprogramming in the metabolism of dormant cells away from Warburg-based energy metabolism. All three computational approaches unanimously predicted that folate metabolism, pyruvate metabolism, and glutamate metabolism, as well as valine/leucine/isoleucine metabolism are likely dysregulated in cancer dormancy. These findings provide new insights and molecular pathway targets on cancer dormancy, comprehensively catalog dormancy-associated metabolic pathways, and inform future research aimed at prevention of cancer recurrence in particular.

De novo individualized disease modules reveal the synthetic penetrance of genes and inform personalized treatment regimens

Current understandings of individual disease etiology and therapeutics are limited despite great need. To fill the gap, we propose a novel computational pipeline that collects potent disease gene cooperative pathways to envision individualized disease etiology and therapies. Our algorithm constructs individualized disease modules de novo, which enables us to elucidate the importance of mutated genes in specific patients and to understand the synthetic penetrance of these genes across patients. We reveal that importance of the notorious cancer drivers TP53 and PIK3CA fluctuate widely across breast cancers and peak in tumors with distinct numbers of mutations and that rarely mutated genes such as XPO1 and PLEKHA1 have high disease module importance in specific individuals. Furthermore, individualized module disruption enables us to devise customized singular and combinatorial target therapies that were highly varied across patients, showing the need for precision therapeutics pipelines. As the first analysis of de novo individualized disease modules, we illustrate the power of individualized disease modules for precision medicine by providing deep novel insights on the activity of diseased genes in individuals.

Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics

Colorectal cancer (CRC) is a major cause of morbidity and mortality in the United States. Tumor-stromal metabolic crosstalk in the tumor microenvironment promotes CRC development and progression, but exactly how stromal cells, in particular cancer-associated fibroblasts (CAFs), affect the metabolism of tumor cells remains unknown. Here we take a data-driven approach to investigate the metabolic interactions between CRC cells and CAFs, integrating constraint-based modeling and metabolomic profiling. Using metabolomics data, we perform unsteady-state parsimonious flux balance analysis to infer flux distributions for central carbon metabolism in CRC cells treated with or without CAF-conditioned media. We find that CAFs reprogram CRC metabolism through stimulation of glycolysis, the oxidative arm of the pentose phosphate pathway (PPP), and glutaminolysis, as well as inhibition of the tricarboxylic acid cycle. To identify potential therapeutic targets, we simulate enzyme knockouts and find that CAF-treated CRC cells are especially sensitive to inhibitions of hexokinase and glucose-6-phosphate, the rate limiting steps of glycolysis and oxidative PPP. Our work gives mechanistic insights into the metabolic interactions between CRC cells and CAFs and provides a framework for testing hypotheses towards CRC-targeted therapies.

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.

Importance of the biomass formulation for cancer metabolic modeling and drug prediction

Genome-scale metabolic reconstructions include all known biochemical reactions occurring in a cell. A typical application is the prediction of potential drug targets for cancer treatment. The precision of these predictions relies on the definition of the objective function. Generally, the biomass reaction is used to illustrate the growth capacity of a cancer cell. Today, seven human biomass reactions can be identified in published metabolic models. The impact of these differences on the metabolic model predictions has not been explored in detail. We explored this impact on cancer metabolic model predictions and showed that the metabolite composition and the associated coefficients had a large impact on the growth rate prediction accuracy, whereas gene essentiality predictions were mainly affected by the metabolite composition. Our results demonstrate the importance of defining a consensus biomass reaction compatible with most human models, which would contribute to ensuring the reproducibility and consistency of the results.

•

The definition of the biomass reaction is of utmost importance for model predictions

•

Growth rate predictions are affected by metabolite composition and their coefficients

•

Gene essentiality predictions are mainly affected by the metabolite composition

•

Need to find a standard biomass reaction for reproducibility and consistency purposes

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

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.

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.

A metabolic CRISPR-Cas9 screen in Chinese hamster ovary cells identifies glutamine-sensitive genes

Media and feed optimization have fueled many-fold improvements in mammalian biopharmaceutical production, but genome editing offers an emerging avenue for further enhancing cell metabolism and bioproduction. However, the complexity of metabolism, involving thousands of genes, makes it unclear which engineering strategies will result in desired traits. Here we present a comprehensive pooled CRISPR screen for CHO cell metabolism, including ~16,000 gRNAs against ~2500 metabolic enzymes and regulators. Using this screen, we identified a glutamine response network in CHO cells. Glutamine is particularly important since it is often over-fed to drive increased TCA cycle flux, but toxic ammonia may accumulate. With the screen we found one orphan glutamine-responsive gene with no clear connection to our network. Knockout of this novel and poorly characterized lipase, Abhd11, substantially increased growth in glutamine-free media by altering the regulation of the TCA cycle. Thus, the screen provides an invaluable targeted platform to comprehensively study genes involved in any metabolic trait, and elucidate novel regulators of metabolism.

Oncoimmunology Meets Organs-on-Chip

Oncoimmunology represents a biomedical research discipline coined to study the roles of immune system in cancer progression with the aim of discovering novel strategies to arm it against the malignancy. Infiltration of immune cells within the tumor microenvironment is an early event that results in the establishment of a dynamic cross-talk. Here, immune cells sense antigenic cues to mount a specific anti-tumor response while cancer cells emanate inhibitory signals to dampen it. Animals models have led to giant steps in this research context, and several tools to investigate the effect of immune infiltration in the tumor microenvironment are currently available. However, the use of animals represents a challenge due to ethical issues and long duration of experiments. Organs-on-chip are innovative tools not only to study how cells derived from different organs interact with each other, but also to investigate on the crosstalk between immune cells and different types of cancer cells. In this review, we describe the state-of-the-art of microfluidics and the impact of OOC in the field of oncoimmunology underlining the importance of this system in the advancements on the complexity of tumor microenvironment.

Diverse metabolic response of cancer cells treated with a 213Bi-anti-EGFR-immunoconjugate

Evaluation of treatment response is among the major challenges in modern oncology. We herein used a monoclonal antibody targeting the EGF receptor (EGFR) labelled with the alpha emitter ²¹³Bi (²¹³Bi-anti-EGFR-MAb). EJ28Luc (bladder) and LN18 (glioma) cancer cells, both overexpressing EGFR, were incubated for 3 h with the radioimmunoconjugate. To assess the responses in the core carbon metabolism upon this treatment, these cancer cell lines were subsequently cultivated for 18 h in the presence of [U-¹³C₆]glucose. ¹³C-enrichment and isotopologue profiles of key amino acids were monitored by gas chromatography–mass spectrometry (GC/MS), in order to monitor the impacts of the radionuclide-treatment upon glucose metabolism. In comparison to untreated controls, treatment of EJ28Luc cells with ²¹³Bi-anti-EGFR-MAb resulted in a significantly decreased incorporation of ¹³C from [U-¹³C₆]glucose into alanine, aspartate, glutamate, glycine, proline and serine. In sharp contrast, the same amino acids did not display less ¹³C-enrichments during treatment of the LN18 cells. The data indicate early treatment response of the bladder cancer cells, but not of the glioma cells though cell lines were killed following ²¹³Bi-anti-EGFR-MAb treatment. The pilot study shows that the ¹³C-labelling approach is a valid tool to assess the responsiveness of cancer cells upon radionuclide-treatment in considerable metabolic detail.

Integrated genetic and metabolic landscapes predict vulnerabilities of temozolomide resistant glioblastoma cells

Metabolic reprogramming and its molecular underpinnings are critical to unravel the duality of cancer cell function and chemo-resistance. Here, we use a constraints-based integrated approach to delineate the interplay between metabolism and epigenetics, hardwired in the genome, to shape temozolomide (TMZ) resistance. Differential metabolism was identified in response to TMZ at varying concentrations in both the resistant neurospheroidal (NSP) and the susceptible (U87MG) glioblastoma cell-lines. The genetic basis of this metabolic adaptation was characterized by whole exome sequencing that identified mutations in signaling pathway regulators of growth and energy metabolism. Remarkably, our integrated approach identified rewiring in glycolysis, TCA cycle, malate aspartate shunt, and oxidative phosphorylation pathways. The differential killing of TMZ resistant NSP by Rotenone at low concentrations with an IC50 value of 5 nM, three orders of magnitude lower than for U87MG that exhibited an IC50 value of 1.8 mM was thus identified using our integrated systems-based approach.

TRIB3 Promotes Lung Adenocarcinoma Progression via an Enhanced Warburg Effect

Background

The pseudokinase Tribbles 3 (TRIB3) is involved in many cellular processes and various cancers. In recent years, the importance of metabolic transformation in the maintenance of malignant tumors has become increasingly prominent. Abnormal metabolism of cancer cells is considered a hallmark of cancer. However, the exact role and molecular mechanism of TRIB3 in lung adenocarcinoma (LUAD) cell reprogramming is largely unknown.

Methods

The oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) of cells were examined with a Seahorse XF Extracellular Flux Analyzer. In vitro and in vivo RT-qPCR, Western blotting, and functional assays were performed to explore the functional roles of TRIB3 in LUAD.

Results

In the present study, we demonstrated that TRIB3 is remarkably upregulated in LUAD cell lines as well as tissues. TRIB3 knockdown significantly inhibited LUAD cell growth and suppressed LUAD cell invasion, while TRIB3 overexpression conferred the opposite effects. Moreover, silencing TRIB3 suppressed the tumorigenesis and metastatic ability of LUAD cells. Mechanistically, we demonstrated that silencing TRIB3 significantly impaired aerobic glycolysis ability in LUAD cells. Furthermore, our data indicated that TRIB3 knockdown decreased hypoxia-inducible factor (HIF)1α levels and targeted the glycolytic genes regulated by HIF1α.

Conclusion

Together, our findings revealed a previously unappreciated function of TRIB3 in cancer cell metabolism and tumor progression, illustrating that TRIB3 could be considered a valuable therapeutic target for LUAD patients.

PluriMetNet: A dynamic electronic model decrypting the metabolic variations in human embryonic stem cells (hESCs) at fluctuating oxygen concentrations

Stem cells are an excellent resource in translational medicine however much is known only in terms of transcriptional and epigenetic regulation of human embryonic stem cells (hESCs). Metabolic regulation of hESCs is still unexplored in many ways, particularly the role of energy metabolism, which is intrinsic to the maintenance of cell viability, however, is very little explored in the past years. Also, there exists no hESC specific core metabolic model of pluripotency as per our knowledge. Through our work, we establish such a metabolic model of hESC using combinatorial in-silico approach of genome scale model reduction and literature curation. Further, through perturbations taking oxygen as a parameter we propose that under lower levels of oxygen concentration there is a significant dynamic change in the energy metabolism of the hESC. We further investigated energy subsystem pathways and their respective reactions in order to locate the direction of energy production along with the dynamic of nutrient metabolites like glucose and glutamine. The output shows a steep increment/decrement at a certain oxygen range. These sharp increments/decrements under hypoxic conditions are termed here as a critical range for hESC metabolic pathway. The data also resonates with the previous experimental studies on hESC energy metabolism confirming the robustness of our model. The model helps to extract range for different pathways in the energy subsystem, making us a little closer in understanding the metabolism of hESC. We also demonstrated the possible range of pathway changes in hESC's energy metabolism that can serve as the crucial preliminary data for further prospective studies. The model also offers a promise in the prediction of the flux behaviour of various metabolites in hESC.Communicated by Ramaswamy H. Sarma.

On the Origin of ATP Synthesis in Cancer

ATP is required for mammalian cells to remain viable and to perform genetically programmed functions. Maintenance of the ΔG′_ATP hydrolysis of −56 kJ/mole is the endpoint of both genetic and metabolic processes required for life. Various anomalies in mitochondrial structure and function prevent maximal ATP synthesis through OxPhos in cancer cells. Little ATP synthesis would occur through glycolysis in cancer cells that express the dimeric form of pyruvate kinase M2. Mitochondrial substrate level phosphorylation (mSLP) in the glutamine-driven glutaminolysis pathway, substantiated by the succinate-CoA ligase reaction in the TCA cycle, can partially compensate for reduced ATP synthesis through both OxPhos and glycolysis. A protracted insufficiency of OxPhos coupled with elevated glycolysis and an auxiliary, fully operational mSLP, would cause a cell to enter its default state of unbridled proliferation with consequent dedifferentiation and apoptotic resistance, i.e., cancer. The simultaneous restriction of glucose and glutamine offers a therapeutic strategy for managing cancer.

Analysis of human metabolism by reducing the complexity of the genome-scale models using redHUMAN

Altered metabolism is associated with many human diseases. Human genome-scale metabolic models (GEMs) were reconstructed within systems biology to study the biochemistry occurring in human cells. However, the complexity of these networks hinders a consistent and concise physiological representation. We present here redHUMAN, a workflow for reconstructing reduced models that focus on parts of the metabolism relevant to a specific physiology using the recently established methods redGEM and lumpGEM. The reductions include the thermodynamic properties of compounds and reactions guaranteeing the consistency of predictions with the bioenergetics of the cell. We introduce a method (redGEMX) to incorporate the pathways used by cells to adapt to the medium. We provide the thermodynamic curation of the human GEMs Recon2 and Recon3D and we apply the redHUMAN workflow to derive leukemia-specific reduced models. The reduced models are powerful platforms for studying metabolic differences between phenotypes, such as diseased and healthy cells.

Quantitative analysis of amino acid metabolism in liver cancer links glutamate excretion to nucleotide synthesis

We used a combination of experimental measurements and computer simulations to understand how liver cancer cells rewire their metabolism to grow faster. We observed that glutamate is excreted by the cells, and our simulations suggest that this occurs because glutamate is formed in excess in the cytoplasm, when cells rapidly synthesize nucleotides, which are required for growth. Meanwhile, glutamate that is formed in the mitochondria is, on the other hand, not excreted. Treating glutamate as two distinct pools, a cytosolic and a mitochondrial, is useful to better understand why many cancer cells rapidly consume glutamine, the precursor of glutamate. The results point toward potential drug targets that could be used to reduce growth of liver cancer cells.

Many cancer cells consume glutamine at high rates; counterintuitively, they simultaneously excrete glutamate, the first intermediate in glutamine metabolism. Glutamine consumption has been linked to replenishment of tricarboxylic acid cycle (TCA) intermediates and synthesis of adenosine triphosphate (ATP), but the reason for glutamate excretion is unclear. Here, we dynamically profile the uptake and excretion fluxes of a liver cancer cell line (HepG2) and use genome-scale metabolic modeling for in-depth analysis. We find that up to 30% of the glutamine is metabolized in the cytosol, primarily for nucleotide synthesis, producing cytosolic glutamate. We hypothesize that excreting glutamate helps the cell to increase the nucleotide synthesis rate to sustain growth. Indeed, we show experimentally that partial inhibition of glutamate excretion reduces cell growth. Our integrative approach thus links glutamine addiction to glutamate excretion in cancer and points toward potential drug targets.

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.

An atlas of human metabolism

Genome-scale metabolic models (GEMs) are valuable tools to study metabolism and provide a scaffold for the integrative analysis of omics data. Researchers have developed increasingly comprehensive human GEMs, but the disconnect among different model sources and versions impedes further progress. We therefore integrated and extensively curated the most recent human metabolic models to construct a consensus GEM, Human1. We demonstrated the versatility of Human1 through the generation and analysis of cell- and tissue-specific models using transcriptomic, proteomic, and kinetic data. We also present an accompanying web portal, Metabolic Atlas (https://www.metabolicatlas.org/), which facilitates further exploration and visualization of Human1 content. Human1 was created using a version-controlled, open-source model development framework to enable community-driven curation and refinement. This framework allows Human1 to be an evolving shared resource for future studies of human health and disease.

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.

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.

Probabilistic controllability approach to metabolic fluxes in normal and cancer tissues

Recent research has shown that many types of cancers take control of specific metabolic processes. We compiled metabolic networks corresponding to four healthy and cancer tissues, and analysed the healthy–cancer transition from the metabolic flux change perspective. We used a Probabilistic Minimum Dominating Set (PMDS) model, which identifies a minimum set of nodes that act as driver nodes and control the entire network. The combination of control theory with flux correlation analysis shows that flux correlations substantially increase in cancer states of breast, kidney and urothelial tissues, but not in lung. No change in the network topology between healthy and cancer networks was observed, but PMDS analysis shows that cancer states require fewer controllers than their corresponding healthy states. These results indicate that cancer metabolism is characterised by more streamlined flux distributions, which may be focused towards a reduced set of objectives and controlled by fewer regulatory elements.

Metabolic rewiring is a feature of many cancers. Here, the authors combine control theory and flux correlation analysis to study the transition of healthy metabolic networks to cancer states, and find that cancer metabolism is characterized by more streamlined flux distributions.

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.

UHRF1 promotes aerobic glycolysis and proliferation via suppression of SIRT4 in pancreatic cancer

UHRF1 (ubiquitin like with plant homeodomain and ring finger domains 1) is an epigenetic modifier that is overexpressed in some cancers, including pancreatic cancer, and mediates silencing of tumor suppressor genes. However, the role of UHRF1 in regulating pancreatic cancer metabolism and metastasis is not clear. In the present study, we demonstrated that silencing UHRF1 significantly inhibited aerobic glycolysis in pancreatic cancer cells. Furthermore, we demonstrated that UHRF1 knockdown decreased hypoxia inducible factor (HIF)1α levels and HIF1α targeted glycolytic genes. The Cancer Genome Atlas dataset analysis supported this observation. The Sirtuin (SIRT) family members regulate aerobic glycolysis in many cancers. We analyzed the correlation between UHRF1 and SIRT3-5 expression and found a significant negative correlation between UHRF1 and SIRT4. Further transcriptional and functional analysis demonstrates that SIRT4 is a downstream target of UHRF1 and negatively regulated aerobic glycolysis, cell proliferation and tumor growth. Our study identified a novel UHRF1/SIRT4 axis in regulation of pancreatic cancer cell proliferation, metabolism, and metastasis.

Hypo-osmotic-like stress underlies general cellular defects of aneuploidy

Aneuploidy, which refers to unbalanced chromosome numbers, represents a class of genetic variation that is associated with cancer, birth defects and eukaryotic micro-organisms1-4. Whereas it is known that each aneuploid chromosome stoichiometry can give rise to a distinct pattern of gene expression and phenotypic profile4,5, it remains a fundamental question as to whether there are common cellular defects that are associated with aneuploidy. Here we show the existence in budding yeast of a common aneuploidy gene-expression signature that is suggestive of hypo-osmotic stress, using a strategy that enables the observation of common transcriptome changes of aneuploidy by averaging out karyotype-specific dosage effects in aneuploid yeast-cell populations with random and diverse chromosome stoichiometry. Consistently, aneuploid yeast exhibited increased plasma-membrane stress that led to impaired endocytosis, and this defect was also observed in aneuploid human cells. Thermodynamic modelling showed that hypo-osmotic-like stress is a general outcome of the proteome imbalance that is caused by aneuploidy, and also predicted a relationship between ploidy and cell size that was observed in yeast and aneuploid cancer cells. A genome-wide screen uncovered a general dependency of aneuploid cells on a pathway of ubiquitin-mediated endocytic recycling of nutrient transporters. Loss of this pathway, coupled with the endocytic defect inherent to aneuploidy, leads to a marked alteration of intracellular nutrient homeostasis.

Systems biology in inflammatory bowel diseases: on the way to precision medicine

Inflammatory bowel diseases (IBD) are chronic and recurrent inflammatory disorders of the gastrointestinal tract. The elucidation of their etiopathology requires complex and multiple approaches. Systems biology has come to fulfill this need in approaching the pathogenetic mechanisms of IBD and its etiopathology, in a comprehensive way, by combining data from different scientific sources. In combination with bioinformatics and network medicine, it uses principles from computer science, mathematics, physics, chemistry, biology, medicine and computational tools to achieve its purposes. Systems biology utilizes scientific sources that provide data from omics studies (e.g., genomics, transcriptomics, etc.) and clinical observations, whose combined analysis leads to network formation and ultimately to a more integrative image of disease etiopathogenesis. In this review, we analyze the current literature on the methods and the tools utilized by systems biology in order to cover an innovative and exciting field: IBD-omics.

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.

Increasing consensus of context-specific metabolic models by integrating data-inferred cell functions

Genome-scale metabolic models provide a valuable context for analyzing data from diverse high-throughput experimental techniques. Models can quantify the activities of diverse pathways and cellular functions. Since some metabolic reactions are only catalyzed in specific environments, several algorithms exist that build context-specific models. However, these methods make differing assumptions that influence the content and associated predictive capacity of resulting models, such that model content varies more due to methods used than cell types. Here we overcome this problem with a novel framework for inferring the metabolic functions of a cell before model construction. For this, we curated a list of metabolic tasks and developed a framework to infer the activity of these functionalities from transcriptomic data. We protected the data-inferred tasks during the implementation of diverse context-specific model extraction algorithms for 44 cancer cell lines. We show that the protection of data-inferred metabolic tasks decreases the variability of models across extraction methods. Furthermore, resulting models better capture the actual biological variability across cell lines. This study highlights the potential of using biological knowledge, inferred from omics data, to obtain a better consensus between existing extraction algorithms. It further provides guidelines for the development of the next-generation of data contextualization methods.

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.