The evolution of genome-scale models of cancer metabolism

MACAW: a method for semi-automatic detection of errors in genome-scale metabolic models

Genome-scale metabolic models (GSMMs) are used to predict metabolic fluxes, with applications ranging from identifying novel drug targets to engineering microbial metabolism. Erroneous or missing reactions, scattered throughout densely interconnected networks, are a limiting factor in these applications. We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a suite of algorithms that helps to identify and visualize errors at the level of connected pathways, rather than individual reactions. We show how MACAW highlights inaccuracies of varying severity in manually curated and automatically generated GSMMs for humans, yeast, and bacteria and helps to identify systematic issues to be addressed in future model construction efforts.

Semi-Automatic Detection of Errors in Genome-Scale Metabolic Models

Background

Genome-Scale Metabolic Models (GSMMs) are used for numerous tasks requiring computational estimates of metabolic fluxes, from predicting novel drug targets to engineering microbes to produce valuable compounds. A key limiting step in most applications of GSMMs is ensuring their representation of the target organism's metabolism is complete and accurate. Identifying and visualizing errors in GSMMs is complicated by the fact that they contain thousands of densely interconnected reactions. Furthermore, many errors in GSMMs only become apparent when considering pathways of connected reactions collectively, as opposed to examining reactions individually.

Results

We present Metabolic Accuracy Check and Analysis Workflow (MACAW), a collection of algorithms for detecting errors in GSMMs. The relative frequencies of errors we detect in manually curated GSMMs appear to reflect the different approaches used to curate them. Changing the method used to automatically create a GSMM from a particular organism's genome can have a larger impact on the kinds of errors in the resulting GSMM than using the same method with a different organism's genome. Our algorithms are particularly capable of identifying errors that are only apparent at the pathway level, including loops, and nontrivial cases of dead ends.

Conclusions

MACAW is capable of identifying inaccuracies of varying severity in a wide range of GSMMs. Correcting these errors can measurably improve the predictive capacity of a GSMM. The relative prevalence of each type of error we identify in a large collection of GSMMs could help shape future efforts for further automation of error correction and GSMM creation.

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

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

Metabolic Connectome and Its Role in the Prediction, Diagnosis, and Treatment of Complex Diseases

The interconnectivity of advanced biological systems is essential for their proper functioning. In modern connectomics, biological entities such as proteins, genes, RNA, DNA, and metabolites are often represented as nodes, while the physical, biochemical, or functional interactions between them are represented as edges. Among these entities, metabolites are particularly significant as they exhibit a closer relationship to an organism's phenotype compared to genes or proteins. Moreover, the metabolome has the ability to amplify small proteomic and transcriptomic changes, even those from minor genomic changes. Metabolic networks, which consist of complex systems comprising hundreds of metabolites and their interactions, play a critical role in biological research by mediating energy conversion and chemical reactions within cells. This review provides an introduction to common metabolic network models and their construction methods. It also explores the diverse applications of metabolic networks in elucidating disease mechanisms, predicting and diagnosing diseases, and facilitating drug development. Additionally, it discusses potential future directions for research in metabolic networks. Ultimately, this review serves as a valuable reference for researchers interested in metabolic network modeling, analysis, and their applications.

Control Theory and Systems Biology: Potential Applications in Neurodegeneration and Search for Therapeutic Targets

Control theory, a well-established discipline in engineering and mathematics, has found novel applications in systems biology. This interdisciplinary approach leverages the principles of feedback control and regulation to gain insights into the complex dynamics of cellular and molecular networks underlying chronic diseases, including neurodegeneration. By modeling and analyzing these intricate systems, control theory provides a framework to understand the pathophysiology and identify potential therapeutic targets. Therefore, this review examines the most widely used control methods in conjunction with genomic-scale metabolic models in the steady state of the multi-omics type. According to our research, this approach involves integrating experimental data, mathematical modeling, and computational analyses to simulate and control complex biological systems. In this review, we find that the most significant application of this methodology is associated with cancer, leaving a lack of knowledge in neurodegenerative models. However, this methodology, mainly associated with the Minimal Dominant Set (MDS), has provided a starting point for identifying therapeutic targets for drug development and personalized treatment strategies, paving the way for more effective therapies.

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.

Cytosolic malic enzyme and glucose-6-phosphate dehydrogenase modulate redox balance in NSCLC with acquired drug resistance

Lung cancer cells often show elevated levels of reactive oxygen species (ROS) and nicotinamide adenine dinucleotide phosphate (NADPH). However, the connections between deregulated redox homeostasis in different subtypes of lung cancer and acquired drug resistance in lung cancer have not yet been fully established. Herein, we analyzed different subtypes of lung cancer data reported in the Cancer Cell Line Encyclopedia (CCLE) database, the Cancer Genome Atlas program (TCGA), and the sequencing data obtained from a gefitinib-resistant non-small-cell lung cancer (NSCLC) cell line (H1975GR). Using flux balance analysis (FBA) model integrated with multiomics data and gene expression profiles, we identified cytosolic malic enzyme 1 (ME1) and glucose-6-phosphate dehydrogenase as the major contributors to the significantly upregulated NADPH flux in NSCLC tissues as compared with normal lung tissues, and gefitinib-resistant NSCLC cell line as compared with the parental cell line. Silencing the gene expression of either of these two enzymes in two osimertinib-resistant NSCLC cell lines (H1975OR and HCC827OR) exhibited strong antiproliferative effects. Our findings not only underscored the pivotal roles of cytosolic ME1 and glucose-6-phosphate dehydrogenase in regulating redox states in NSCLC cells but also provided novel insights into their potential roles in drug-resistant NSCLC cells with disturbed redox states.

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.

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.

Identifying essential genes in genome-scale metabolic models of consensus molecular subtypes of colorectal cancer

Identifying essential targets in the genome-scale metabolic networks of cancer cells is a time-consuming process. The present study proposed a fuzzy hierarchical optimization framework for identifying essential genes, metabolites and reactions. On the basis of four objectives, the present study developed a framework for identifying essential targets that lead to cancer cell death and evaluating metabolic flux perturbations in normal cells that have been caused by cancer treatment. Through fuzzy set theory, a multiobjective optimization problem was converted into a trilevel maximizing decision-making (MDM) problem. We applied nested hybrid differential evolution to solve the trilevel MDM problem to identify essential targets in genome-scale metabolic models for five consensus molecular subtypes (CMSs) of colorectal cancer. We used various media to identify essential targets for each CMS and discovered that most targets affected all five CMSs and that some genes were CMS-specific. We obtained experimental data on the lethality of cancer cell lines from the DepMap database to validate the identified essential genes. The results reveal that most of the identified essential genes were compatible with the colorectal cancer cell lines obtained from DepMap and that these genes, with the exception of EBP, LSS, and SLC7A6, could generate a high level of cell death when knocked out. The identified essential genes were mostly involved in cholesterol biosynthesis, nucleotide metabolisms, and the glycerophospholipid biosynthetic pathway. The genes involved in the cholesterol biosynthetic pathway were also revealed to be determinable, if a cholesterol uptake reaction was not induced when the cells were in the culture medium. However, the genes involved in the cholesterol biosynthetic pathway became non-essential if such a reaction was induced. Furthermore, the essential gene CRLS1 was revealed as a medium-independent target for all CMSs.

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.

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

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

Expanding biochemical knowledge and illuminating metabolic dark matter with ATLASx

Metabolic “dark matter” describes currently unknown metabolic processes, which form a blind spot in our general understanding of metabolism and slow down the development of biosynthetic cell factories and naturally derived pharmaceuticals. Mapping the dark matter of metabolism remains an open challenge that can be addressed globally and systematically by existing computational solutions. In this work, we use 489 generalized enzymatic reaction rules to map both known and unknown metabolic processes around a biochemical database of 1.5 million biological compounds. We predict over 5 million reactions and integrate nearly 2 million naturally and synthetically-derived compounds into the global network of biochemical knowledge, named ATLASx. ATLASx is available to researchers as a powerful online platform that supports the prediction and analysis of biochemical pathways and evaluates the biochemical vicinity of molecule classes (https://lcsb-databases.epfl.ch/Atlas2).

“Mapping the dark matter of metabolism remains an open challenge that can be addressed globally and systematically by existing computational solutions. Here the authors present ATLASx, a repository of known and predicted enzymatic reaction, connecting millions of compounds to help synthetic biologists and metabolic engineers to design and explore metabolic pathways.”

ΔFBA-Predicting metabolic flux alterations using genome-scale metabolic models and differential transcriptomic data

Genome-scale metabolic models (GEMs) provide a powerful framework for simulating the entire set of biochemical reactions in a cell using a constraint-based modeling strategy called flux balance analysis (FBA). FBA relies on an assumed metabolic objective for generating metabolic fluxes using GEMs. But, the most appropriate metabolic objective is not always obvious for a given condition and is likely context-specific, which often complicate the estimation of metabolic flux alterations between conditions. Here, we propose a new method, called ΔFBA (deltaFBA), that integrates differential gene expression data to evaluate directly metabolic flux differences between two conditions. Notably, ΔFBA does not require specifying the cellular objective. Rather, ΔFBA seeks to maximize the consistency and minimize inconsistency between the predicted flux differences and differential gene expression. We showcased the performance of ΔFBA through several case studies involving the prediction of metabolic alterations caused by genetic and environmental perturbations in Escherichia coli and caused by Type-2 diabetes in human muscle. Importantly, in comparison to existing methods, ΔFBA gives a more accurate prediction of flux differences.

Metabolic alterations are often used as hallmarks of observable phenotypes. In this regard, reconstructed genome-scale metabolic models (GEMs) provide a rich and computable representation of the entire set of biochemical reactions in a cell. However, the performance of analytical tools for predicting metabolic reaction rates or fluxes using GEMs is sensitive to the assumed metabolic objective that is often unknown and likely context-specific. Here, we propose a novel method called ΔFBA that combines differential gene expression data and GEMs to evaluate differences in the metabolic fluxes between two conditions (perturbation vs. control) without the need for specifying a metabolic objective. In our demonstration, ΔFBA outperformed other existing methods in predicting metabolic flux alterations.

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 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.

Metabolic control of cancer progression as novel targets for therapy

Metabolism is an important part of tumorigenesis as well as progression. The various cancer metabolism pathways, such as glucose metabolism and glutamine metabolism, directly regulate the development and progression of cancer. The pathways by which the cancer cells rewire their metabolism according to their needs, surrounding environment and host tissue conditions are an important area of study. The regulation of these metabolic pathways is determined by various oncogenes, tumor suppressor genes, as well as various constituent cells of the tumor microenvironment. Expanded studies on metabolism will help identify efficient biomarkers for diagnosis and strategies for therapeutic interventions and countering ways by which cancers may acquire resistance to therapy.

Exploring the Metabolic Heterogeneity of Cancers: A Benchmark Study of Context-Specific Models

Metabolic heterogeneity is a hallmark of cancer and can distinguish a normal phenotype from a cancer phenotype. In the systems biology domain, context-specific models facilitate extracting physiologically relevant information from high-quality data. Here, to utilize the heterogeneity of metabolic patterns to discover biomarkers of all cancers, we benchmarked thousands of context-specific models using well-established algorithms for the integration of omics data into the generic human metabolic model Recon3D. By analyzing the active reactions capable of carrying flux and their magnitude through flux balance analysis, we proved that the metabolic pattern of each cancer is unique and could act as a cancer metabolic fingerprint. Subsequently, we searched for proper feature selection methods to cluster the flux states characterizing each cancer. We employed PCA-based dimensionality reduction and a random forest learning algorithm to reveal reactions containing the most relevant information in order to effectively identify the most influential fluxes. Conclusively, we discovered different pathways that are probably the main sources for metabolic heterogeneity in cancers. We designed the GEMbench website to interactively present the data, methods, and analysis results.

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.

On the Use of Topological Features of Metabolic Networks for the Classification of Cancer Samples

BACKGROUND

The increasing availability of omics data collected from patients affected by severe pathologies, such as cancer, is fostering the development of data science methods for their analysis.

INTRODUCTION

The combination of data integration and machine learning approaches can provide new powerful instruments to tackle the complexity of cancer development and deliver effective diagnostic and prognostic strategies.

METHODS

We explore the possibility of exploiting the topological properties of sample-specific metabolic networks as features in a supervised classification task. Such networks are obtained by projecting transcriptomic data from RNA-seq experiments on genome-wide metabolic models to define weighted networks modeling the overall metabolic activity of a given sample.

RESULTS

We show the classification results on a labeled breast cancer dataset from the TCGA database, including 210 samples (cancer vs. normal). In particular, we investigate how the performance is affected by a threshold-based pruning of the networks by comparing Artificial Neural Networks, Support Vector Machines and Random Forests. Interestingly, the best classification performance is achieved within a small threshold range for all methods, suggesting that it might represent an effective choice to recover useful information while filtering out noise from data. Overall, the best accuracy is achieved with SVMs, which exhibit performances similar to those obtained when gene expression profiles are used as features.

CONCLUSION

These findings demonstrate that the topological properties of sample-specific metabolic networks are effective in classifying cancer and normal samples, suggesting that useful information can be extracted from a relatively limited number of features.

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

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

Mouse Model for Efficient Simultaneous Targeting of Glycolysis, Glutaminolysis, and De Novo Synthesis of Fatty Acids in Colon Cancer

Colon cancer is a highly anabolic entity with upregulation of glycolysis, glutaminolysis, and de novo synthesis of fatty acids, which also induces a hypercatabolic state in the patient. The blockade of either cancer anabolism or host catabolism has been previously proven to be a successful anticancer experimental treatment. However, it is still unclear whether the simultaneous blockade of both metabolic counterparts can limit malignant survival and the energetic consequences of such an approach. In this chapter, by using the CT26.WT murine colon adenocarcinoma cell line as a model of study, we provide a method to simultaneously perform a pharmacological blockade of tumor anabolism and host catabolism, as a feasible therapeutic approach to treat cancer, and to limit its energetic supply.

An Integrative Gene Expression and Mathematical Flux Balance Analysis Identifies Targetable Redox Vulnerabilities in Melanoma Cells

Melanomas harboring BRAF mutations can be treated with BRAF inhibitors (BRAFi), but responses are varied and tumor recurrence is inevitable. Here we used an integrative approach of experimentation and mathematical flux balance analyses in BRAF-mutated melanoma cells to discover that elevated antioxidant capacity is linked to BRAFi sensitivity in melanoma cells. High levels of antioxidant metabolites in cells with reduced BRAFi sensitivity confirmed this conclusion. By extending our analyses to other melanoma subtypes in The Cancer Genome Atlas, we predict that elevated redox capacity is a general feature of melanomas, not previously observed. We propose that redox vulnerabilities could be exploited for therapeutic benefits and identify unsuspected combination targets to enhance the effects of BRAFi in any melanoma, regardless of mutational status. SIGNIFICANCE: An integrative bioinformatics, flux balance analysis, and experimental approach identify targetable redox vulnerabilities and show the potential for modulation of cancer antioxidant defense to augment the benefits of existing therapies in melanoma.

The combination of orlistat, lonidamine and 6-diazo-5-oxo-L-norleucine induces a quiescent energetic phenotype and limits substrate flexibility in colon cancer cells

Cancer upregulates glycolysis, glutaminolysis and lipogenesis, and induces a catabolic state in patients. The concurrent inhibition of both tumor anabolism and host catabolism, and the energetic consequences of such an approach, have not previously been fully investigated. In the present study, CT26.WT murine colon cancer cells were treated with the combination of anti-anabolic drugs orlistat, lonidamine and 6-diazo-5-oxo-L-norleucine (DON; OLD scheme), which are inhibitors of the de novo synthesis of fatty acids, glycolysis and glutaminolysis, respectively. In addition, the effects of OLD scheme sumplemented with the combination of anti-catabolic compounds, namely growth hormone, insulin and indomethacin (GII scheme), were also evaluated. The effects of the compounds used in combination on CT26.WT cell viability, clonogenicity and energetic metabolism were assessed in vitro. The results demonstrated that the anti-anabolic approach reduced cell viability, clonogenicity and cell cycle progression, and increased apoptosis. These effects were associated with decreased oxidative phosphorylation, glycolysis and fuel flexibility. Furthermore, the anti-catabolic scheme, alone or supplemented with anti-anabolic compounds, did not favor tumor growth. These findings indicated that the simultaneous pharmacological inhibition of tumor anabolism and host catabolism exhibits antitumor effects that should be further evaluated.

StanDep: Capturing transcriptomic variability improves context-specific metabolic models

Diverse algorithms can integrate transcriptomics with genome-scale metabolic models (GEMs) to build context-specific metabolic models. These algorithms require identification of a list of high confidence (core) reactions from transcriptomics, but parameters related to identification of core reactions, such as thresholding of expression profiles, can significantly change model content. Importantly, current thresholding approaches are burdened with setting singular arbitrary thresholds for all genes; thus, resulting in removal of enzymes needed in small amounts and even many housekeeping genes. Here, we describe StanDep, a novel heuristic method for using transcriptomics to identify core reactions prior to building context-specific metabolic models. StanDep clusters gene expression data based on their expression pattern across different contexts and determines thresholds for each cluster using data-dependent statistics, specifically standard deviation and mean. To demonstrate the use of StanDep, we built hundreds of models for the NCI-60 cancer cell lines. These models successfully increased the inclusion of housekeeping reactions, which are often lost in models built using standard thresholding approaches. Further, StanDep also provided a transcriptomic explanation for inclusion of lowly expressed reactions that were otherwise only supported by model extraction methods. Our study also provides novel insights into how cells may deal with context-specific and ubiquitous functions. StanDep, as a MATLAB toolbox, is available at https://github.com/LewisLabUCSD/StanDep.

Elucidating the Reprograming of Colorectal Cancer Metabolism Using Genome-Scale Metabolic Modeling

Colorectal cancer is the third most incidental cancer worldwide, and the response rate of current treatment for colorectal cancer is very low. Genome-scale metabolic models (GEMs) are systems biology platforms, and they had been used to assist researchers in understanding the metabolic alterations in different types of cancer. Here, we reconstructed a generic colorectal cancer GEM by merging 374 personalized GEMs from the Human Pathology Atlas and used it as a platform for systematic investigation of the difference between tumor and normal samples. The reconstructed model revealed the metabolic reprogramming in glutathione as well as the arginine and proline metabolism in response to tumor occurrence. In addition, six genes including ODC1, SMS, SRM, RRM2, SMOX, and SAT1 associated with arginine and proline metabolism were found to be key players in this metabolic alteration. We also investigated these genes in independent colorectal cancer patients and cell lines and found that many of these genes showed elevated level in colorectal cancer and exhibited adverse effect in patients. Therefore, these genes could be promising therapeutic targets for treatment of a specific colon cancer patient group.

Genome-scale metabolic model of the rat liver predicts effects of diet restriction

Mapping network analysis in cells and tissues can provide insights into metabolic adaptations to changes in external environment, pathological conditions, and nutrient deprivation. Here, we reconstructed a genome-scale metabolic network of the rat liver that will allow for exploration of systems-level physiology. The resulting in silico model (iRatLiver) contains 1,882 reactions, 1,448 metabolites, and 994 metabolic genes. We then used this model to characterize the response of the liver’s energy metabolism to a controlled perturbation in diet. Transcriptomics data were collected from the livers of Sprague Dawley rats at 4 or 14 days of being subjected to 15%, 30%, or 60% diet restriction. These data were integrated with the iRatLiver model to generate condition-specific metabolic models, allowing us to explore network differences under each condition. We observed different pathway usage between early and late time points. Network analysis identified several highly connected “hub” genes (Pklr, Hadha, Tkt, Pgm1, Tpi1, and Eno3) that showed differing trends between early and late time points. Taken together, our results suggest that the liver’s response varied with short- and long-term diet restriction. More broadly, we anticipate that the iRatLiver model can be exploited further to study metabolic changes in the liver under other conditions such as drug treatment, infection, and disease.

Machine and deep learning meet genome-scale metabolic modeling

Omic data analysis is steadily growing as a driver of basic and applied molecular biology research. Core to the interpretation of complex and heterogeneous biological phenotypes are computational approaches in the fields of statistics and machine learning. In parallel, constraint-based metabolic modeling has established itself as the main tool to investigate large-scale relationships between genotype, phenotype, and environment. The development and application of these methodological frameworks have occurred independently for the most part, whereas the potential of their integration for biological, biomedical, and biotechnological research is less known. Here, we describe how machine learning and constraint-based modeling can be combined, reviewing recent works at the intersection of both domains and discussing the mathematical and practical aspects involved. We overlap systematic classifications from both frameworks, making them accessible to nonexperts. Finally, we delineate potential future scenarios, propose new joint theoretical frameworks, and suggest concrete points of investigation for this joint subfield. A multiview approach merging experimental and knowledge-driven omic data through machine learning methods can incorporate key mechanistic information in an otherwise biologically-agnostic learning process.

Assessing key decisions for transcriptomic data integration in biochemical networks

To gain insights into complex biological processes, genome-scale data (e.g., RNA-Seq) are often overlaid on biochemical networks. However, many networks do not have a one-to-one relationship between genes and network edges, due to the existence of isozymes and protein complexes. Therefore, decisions must be made on how to overlay data onto networks. For example, for metabolic networks, these decisions include (1) how to integrate gene expression levels using gene-protein-reaction rules, (2) the approach used for selection of thresholds on expression data to consider the associated gene as "active", and (3) the order in which these steps are imposed. However, the influence of these decisions has not been systematically tested. We compared 20 decision combinations using a transcriptomic dataset across 32 tissues and showed that definition of which reaction may be considered as active (i.e., reactions of the genome-scale metabolic network with a non-zero expression level after overlaying the data) is mainly influenced by thresholding approach used. To determine the most appropriate decisions, we evaluated how these decisions impact the acquisition of tissue-specific active reaction lists that recapitulate organ-system tissue groups. These results will provide guidelines to improve data analyses with biochemical networks and facilitate the construction of context-specific metabolic models.

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.

Identifying and targeting cancer-specific metabolism with network-based drug target prediction

Background

Metabolic rewiring allows cancer cells to sustain high proliferation rates. Thus, targeting only the cancer-specific cellular metabolism will safeguard healthy tissues.

Methods

We developed the very efficient FASTCORMICS RNA-seq workflow (rFASTCORMICS) to build 10,005 high-resolution metabolic models from the TCGA dataset to capture metabolic rewiring strategies in cancer cells. Colorectal cancer (CRC) was used as a test case for a repurposing workflow based on rFASTCORMICS.

Findings

Alternative pathways that are not required for proliferation or survival tend to be shut down and, therefore, tumours display cancer-specific essential genes that are significantly enriched for known drug targets. We identified naftifine, ketoconazole, and mimosine as new potential CRC drugs, which were experimentally validated.

Interpretation

The here presented rFASTCORMICS workflow successfully reconstructs a metabolic model based on RNA-seq data and successfully predicted drug targets and drugs not yet indicted for colorectal cancer.

Fund

This study was supported by the University of Luxembourg (IRP grant scheme; R-AGR-0755-12), the Luxembourg National Research Fund (FNR PRIDE PRIDE15/10675146/CANBIO), the Fondation Cancer (Luxembourg), the European Union‘s Horizon2020 research and innovation programme under the Marie Sklodowska- Curie grant agreement No 642295 (MEL-PLEX), and the German Federal Ministry of Education and Research (BMBF) within the project MelanomSensitivity (BMBF/BM/7643621).

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.

Potential premalignant status of gastric portion excluded after Roux en-Y gastric bypass in obese women: A pilot study

We evaluated whether the excluded stomach (ES) after Roux-en-Y gastric bypass (RYGB) can represent a premalignant environment. Twenty obese women were prospectively submitted to double-balloon enteroscopy (DBE) with gastric juice and biopsy collection, before and 3 months after RYGB. We then evaluated morphological and molecular changes by combining endoscopic and histopathological analyses with an integrated untargeted metabolomics and transcriptomics multiplatform. Preoperatively, 16 women already presented with gastric histopathological alterations and an increased pH (≥4.0). These gastric abnormalities worsened after RYGB. A 90-fold increase in the concentration of bile acids was found in ES fluid, which also contained other metabolites commonly found in the intestinal environment, urine, and faeces. In addition, 135 genes were differentially expressed in ES tissue. Combined analysis of metabolic and gene expression data suggested that RYGB promoted activation of biological processes involved in local inflammation, bacteria overgrowth, and cell proliferation sustained by genes involved in carcinogenesis. Accumulated fluid in the ES appears to behave as a potential premalignant environment due to worsening inflammation and changing gene expression patterns that are favorable to the development of cancer. Considering that ES may remain for the rest of the patient’s life, long-term ES monitoring is therefore recommended for patients undergoing RYGB.

Metabolite systems profiling identifies exploitable weaknesses in retinoblastoma

Retinoblastoma (RB) is a childhood eye cancer. Currently, chemotherapy, local therapy, and enucleation are the main ways in which these tumors are managed. The present work is the first study that uses constraint-based reconstruction and analysis approaches to identify and explain RB-specific survival strategies, which are RB tumor specific. Importantly, our model-specific secretion profile is also found in RB1-depleted human retinal cells in vitro and suggests that novel biomarkers involved in lipid metabolism may be important. Finally, RB-specific synthetic lethals have been predicted as lipid and nucleoside transport proteins that can aid in novel drug target development.

Metabolic Flexibility as an Adaptation to Energy Resources and Requirements in Health and Disease

The ability to efficiently adapt metabolism by substrate sensing, trafficking, storage, and utilization, dependent on availability and requirement, is known as metabolic flexibility. In this review, we discuss the breadth and depth of metabolic flexibility and its impact on health and disease. Metabolic flexibility is essential to maintain energy homeostasis in times of either caloric excess or caloric restriction, and in times of either low or high energy demand, such as during exercise. The liver, adipose tissue, and muscle govern systemic metabolic flexibility and manage nutrient sensing, uptake, transport, storage, and expenditure by communication via endocrine cues. At a molecular level, metabolic flexibility relies on the configuration of metabolic pathways, which are regulated by key metabolic enzymes and transcription factors, many of which interact closely with the mitochondria. Disrupted metabolic flexibility, or metabolic inflexibility, however, is associated with many pathological conditions including metabolic syndrome, type 2 diabetes mellitus, and cancer. Multiple factors such as dietary composition and feeding frequency, exercise training, and use of pharmacological compounds, influence metabolic flexibility and will be discussed here. Last, we outline important advances in metabolic flexibility research and discuss medical horizons and translational aspects.

GSEA of mouse and human mitochondriomes reveals fatty acid oxidation in astrocytes

The prevalent view in neuroenergetics is that glucose is the main brain fuel, with neurons being mostly oxidative and astrocytes glycolytic. Evidence supporting that astrocyte mitochondria are functional has been overlooked. Here we sought to determine what is unique about astrocyte mitochondria by performing unbiased statistical comparisons of the mitochondriome in astrocytes and neurons. Using MitoCarta, a compendium of mitochondrial proteins, together with transcriptomes of mouse neurons and astrocytes, we generated cell-specific databases of nuclear genes encoding for mitochondrion proteins, ranked according to relative expression. Standard and in-house Gene Set Enrichment Analyses (GSEA) of five mouse transcriptomes revealed that genes encoding for enzymes involved in fatty acid oxidation (FAO) and amino acid catabolism are consistently more expressed in astrocytes than in neurons. FAO and oxidative-metabolism-related genes are also up-regulated in human cortical astrocytes versus the whole cortex, and in adult astrocytes versus fetal astrocytes. We thus present the first evidence of FAO in human astrocytes. Further, as shown in vitro, FAO coexists with glycolysis in astrocytes and is inhibited by glutamate. Altogether, these analyses provide arguments against the glucose-centered view of energy metabolism in astrocytes and reveal mitochondria as specialized organelles in these cells.

Crosstalk and the Dynamical Modularity of Feed-Forward Loops in Transcriptional Regulatory Networks

Network motifs, such as the feed-forward loop (FFL), introduce a range of complex behaviors to transcriptional regulatory networks, yet such properties are typically determined from their isolated study. We characterize the effects of crosstalk on FFL dynamics by modeling the cross regulation between two different FFLs and evaluate the extent to which these patterns occur in vivo. Analytical modeling suggests that crosstalk should overwhelmingly affect individual protein-expression dynamics. Counter to this expectation we find that entire FFLs are more likely than expected to resist the effects of crosstalk (≈20% for one crosstalk interaction) and remain dynamically modular. The likelihood that cross-linked FFLs are dynamically correlated increases monotonically with additional crosstalk, but is independent of the specific regulation type or connectivity of the interactions. Just one additional regulatory interaction is sufficient to drive the FFL dynamics to a statistically different state. Despite the potential for modularity between sparsely connected network motifs, Escherichia coli (E. coli) appears to favor crosstalk wherein at least one of the cross-linked FFLs remains modular. A gene ontology analysis reveals that stress response processes are significantly overrepresented in the cross-linked motifs found within E. coli. Although the daunting complexity of biological networks affects the dynamical properties of individual network motifs, some resist and remain modular, seemingly insulated from extrinsic perturbations—an intriguing possibility for nature to consistently and reliably provide certain network functionalities wherever the need arise.

Unravelling the complexity of signalling networks in cancer: A review of the increasing role for computational modelling

Cancer induction is a highly complex process involving hundreds of different inducers but whose eventual outcome is the same. Clearly, it is essential to understand how signalling pathways and networks generated by these inducers interact to regulate cell behaviour and create the cancer phenotype. While enormous strides have been made in identifying key networking profiles, the amount of data generated far exceeds our ability to understand how it all "fits together". The number of potential interactions is astronomically large and requires novel approaches and extreme computation methods to dissect them out. However, such methodologies have high intrinsic mathematical and conceptual content which is difficult to follow. This review explains how computation modelling is progressively finding solutions and also revealing unexpected and unpredictable nano-scale molecular behaviours extremely relevant to how signalling and networking are coherently integrated. It is divided into linked sections illustrated by numerous figures from the literature describing different approaches and offering visual portrayals of networking and major conceptual advances in the field. First, the problem of signalling complexity and data collection is illustrated for only a small selection of known oncogenes. Next, new concepts from biophysics, molecular behaviours, kinetics, organisation at the nano level and predictive models are presented. These areas include: visual representations of networking, Energy Landscapes and energy transfer/dissemination (entropy); diffusion, percolation; molecular crowding; protein allostery; quinary structure and fractal distributions; energy management, metabolism and re-examination of the Warburg effect. The importance of unravelling complex network interactions is then illustrated for some widely-used drugs in cancer therapy whose interactions are very extensive. Finally, use of computational modelling to develop micro- and nano- functional models ("bottom-up" research) is highlighted. The review concludes that computational modelling is an essential part of cancer research and is vital to understanding network formation and molecular behaviours that are associated with it. Its role is increasingly essential because it is unravelling the huge complexity of cancer induction otherwise unattainable by any other approach.

"Gestaltomics": Systems Biology Schemes for the Study of Neuropsychiatric Diseases

The integration of different sources of biological information about what defines a behavioral phenotype is difficult to unify in an entity that reflects the arithmetic sum of its individual parts. In this sense, the challenge of Systems Biology for understanding the “psychiatric phenotype” is to provide an improved vision of the shape of the phenotype as it is visualized by “Gestalt” psychology, whose fundamental axiom is that the observed phenotype (behavior or mental disorder) will be the result of the integrative composition of every part. Therefore, we propose the term “Gestaltomics” as a term from Systems Biology to integrate data coming from different sources of information (such as the genome, transcriptome, proteome, epigenome, metabolome, phenome, and microbiome). In addition to this biological complexity, the mind is integrated through multiple brain functions that receive and process complex information through channels and perception networks (i.e., sight, ear, smell, memory, and attention) that in turn are programmed by genes and influenced by environmental processes (epigenetic). Today, the approach of medical research in human diseases is to isolate one disease for study; however, the presence of an additional disease (co-morbidity) or more than one disease (multimorbidity) adds complexity to the study of these conditions. This review will present the challenge of integrating psychiatric disorders at different levels of information (Gestaltomics). The implications of increasing the level of complexity, for example, studying the co-morbidity with another disease such as cancer, will also be discussed.

A Systematic Evaluation of Methods for Tailoring Genome-Scale Metabolic Models

Genome-scale models of metabolism can illuminate the molecular basis of cell phenotypes. Since some enzymes are only active in specific cell types, several algorithms use omics data to construct cell-line- and tissue-specific metabolic models from genome-scale models. However, these methods are often not rigorously benchmarked, and it is unclear how algorithm and parameter selection (e.g., gene expression thresholds, metabolic constraints) affects model content and predictive accuracy. To investigate this, we built hundreds of models of four different cancer cell lines using six algorithms, four gene expression thresholds, and three sets of metabolic constraints. Model content varied substantially across different parameter sets, but the algorithms generally increased accuracy in gene essentiality predictions. However, model extraction method choice had the largest impact on model accuracy. We further highlight how assumptions during model development influence model prediction accuracy. These insights will guide further development of context-specific models, thus more accurately resolving genotype-phenotype relationships.

Systems biology analysis of drivers underlying hallmarks of cancer cell metabolism

Malignant transformation is often accompanied by significant metabolic changes. To identify drivers underlying these changes, we calculated metabolic flux states for the NCI60 cell line collection and correlated the variance between metabolic states of these lines with their other properties. The analysis revealed a remarkably consistent structure underlying high flux metabolism. The three primary uptake pathways, glucose, glutamine and serine, are each characterized by three features: (1) metabolite uptake sufficient for the stoichiometric requirement to sustain observed growth, (2) overflow metabolism, which scales with excess nutrient uptake over the basal growth requirement, and (3) redox production, which also scales with nutrient uptake but greatly exceeds the requirement for growth. We discovered that resistance to chemotherapeutic drugs in these lines broadly correlates with the amount of glucose uptake. These results support an interpretation of the Warburg effect and glutamine addiction as features of a growth state that provides resistance to metabolic stress through excess redox and energy production. Furthermore, overflow metabolism observed may indicate that mitochondrial catabolic capacity is a key constraint setting an upper limit on the rate of cofactor production possible. These results provide a greater context within which the metabolic alterations in cancer can be understood.

The Emerging Facets of Non-Cancerous Warburg Effect

The Warburg effect (WE), or aerobic glycolysis, is commonly recognized as a hallmark of cancer and has been extensively studied for potential anti-cancer therapeutics development. Beyond cancer, the WE plays an important role in many other cell types involved in immunity, angiogenesis, pluripotency, and infection by pathogens (e.g., malaria). Here, we review the WE in non-cancerous context as a "hallmark of rapid proliferation." We observe that the WE operates in rapidly dividing cells in normal and pathological states that are triggered by internal and external cues. Aerobic glycolysis is also the preferred metabolic program in the cases when robust transient responses are needed. We aim to draw attention to the potential of computational modeling approaches in systematic characterization of common metabolic features beyond the WE across physiological and pathological conditions. Identification of metabolic commonalities across various diseases may lead to successful repurposing of drugs and biomarkers.

Tumor Metabolism, the Ketogenic Diet and β-Hydroxybutyrate: Novel Approaches to Adjuvant Brain Tumor Therapy

Malignant brain tumors are devastating despite aggressive treatments such as surgical resection, chemotherapy and radiation therapy. The average life expectancy of patients with newly diagnosed glioblastoma is approximately ~18 months. It is clear that increased survival of brain tumor patients requires the design of new therapeutic modalities, especially those that enhance currently available treatments and/or limit tumor growth. One novel therapeutic arena is the metabolic dysregulation that results in an increased need for glucose in tumor cells. This phenomenon suggests that a reduction in tumor growth could be achieved by decreasing glucose availability, which can be accomplished through pharmacological means or through the use of a high-fat, low-carbohydrate ketogenic diet (KD). The KD, as the name implies, also provides increased blood ketones to support the energy needs of normal tissues. Preclinical work from a number of laboratories has shown that the KD does indeed reduce tumor growth in vivo. In addition, the KD has been shown to reduce angiogenesis, inflammation, peri-tumoral edema, migration and invasion. Furthermore, this diet can enhance the activity of radiation and chemotherapy in a mouse model of glioma, thus increasing survival. Additional studies in vitro have indicated that increasing ketones such as β-hydroxybutyrate (βHB) in the absence of glucose reduction can also inhibit cell growth and potentiate the effects of chemotherapy and radiation. Thus, while we are only beginning to understand the pluripotent mechanisms through which the KD affects tumor growth and response to conventional therapies, the emerging data provide strong support for the use of a KD in the treatment of malignant gliomas. This has led to a limited number of clinical trials investigating the use of a KD in patients with primary and recurrent glioma.

Genome scale metabolic modeling of cancer

Cancer cells reprogram metabolism to support rapid proliferation and survival. Energy metabolism is particularly important for growth and genes encoding enzymes involved in energy metabolism are frequently altered in cancer cells. A genome scale metabolic model (GEM) is a mathematical formalization of metabolism which allows simulation and hypotheses testing of metabolic strategies. It has successfully been applied to many microorganisms and is now used to study cancer metabolism. Generic models of human metabolism have been reconstructed based on the existence of metabolic genes in the human genome. Cancer specific models of metabolism have also been generated by reducing the number of reactions in the generic model based on high throughput expression data, e.g. transcriptomics and proteomics. Targets for drugs and bio markers for diagnostics have been identified using these models. They have also been used as scaffolds for analysis of high throughput data to allow mechanistic interpretation of changes in expression. Finally, GEMs allow quantitative flux predictions using flux balance analysis (FBA). Here we critically review the requirements for successful FBA simulations of cancer cells and discuss the symmetry between the methods used for modeling of microbial and cancer metabolism. GEMs have great potential for translational research on cancer and will therefore become of increasing importance in the future.

Modulating carbohydrate-protein interactions through glycoengineering of monoclonal antibodies to impact cancer physiology

Diverse glycans on proteins impact cell and organism physiology, along with drug activity. Since many protein-based biotherapeutics are glycosylated and these glycans have biological activity, there is a desire to engineer glycosylation for recombinant protein-based biotherapeutics. Engineered glycosylation can impact the recombinant protein efficacy and also influence many cell pathways by first changing glycan-protein interactions and consequently modulating disease physiologies. However, its complexity is enormous. Recent advances in glycoengineering now make it easier to modulate protein-glycan interactions. Here, we discuss how engineered glycans contribute to therapeutic monoclonal antibodies (mAbs) in the treatment of cancers, how these glycoengineered therapeutic mAbs affect the transformed phenotypes and downstream cell pathways. Furthermore, we suggest how systems biology can help in the next generation mAb glycoengineering process by aiding in data analysis and guiding engineering efforts to tailor mAb glycan and ultimately drug efficacy, safety and affordability.

Understanding the physiology of the ageing individual: computational modelling of changes in metabolism and endurance

Ageing and lifespan are strongly affected by metabolism. The maximal possible uptake of oxygen is not only a good predictor of performance in endurance sports, but also of life expectancy. Figuratively speaking, healthy ageing is a competitive sport. Although the root cause of ageing is damage to macromolecules, it is the balance with repair processes that is decisive. Reduced or intermittent nutrition, hormones and intracellular signalling pathways that regulate metabolism have strong effects on ageing. Homeostatic regulatory processes tend to keep the environment of the cells within relatively narrow bounds. On the other hand, the body is constantly adapting to physical activity and food consumption. Spontaneous fluctuations in heart rate and other processes indicate youth and health. A (homeo)dynamic aspect of homeostasis deteriorates with age. We are now in a position to develop computational models of human metabolism and the dynamics of heart rhythm and oxygen transport that will advance our understanding of ageing. Computational modelling of the connections between dietary restriction, metabolism and protein turnover may increase insight into homeostasis of the proteins in our body. In this way, the computational reconstruction of human physiological processes, the Physiome, can help prevent frailty and age-related disease.