Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling

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.

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

BACKGROUND

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

METHODS

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

RESULTS AND CONCLUSIONS

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

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

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

Genetically personalised organ-specific metabolic models in health and disease

Understanding how genetic variants influence disease risk and complex traits (variant-to-function) is one of the major challenges in human genetics. Here we present a model-driven framework to leverage human genome-scale metabolic networks to define how genetic variants affect biochemical reaction fluxes across major human tissues, including skeletal muscle, adipose, liver, brain and heart. As proof of concept, we build personalised organ-specific metabolic flux models for 524,615 individuals of the INTERVAL and UK Biobank cohorts and perform a fluxome-wide association study (FWAS) to identify 4312 associations between personalised flux values and the concentration of metabolites in blood. Furthermore, we apply FWAS to identify 92 metabolic fluxes associated with the risk of developing coronary artery disease, many of which are linked to processes previously described to play in role in the disease. Our work demonstrates that genetically personalised metabolic models can elucidate the downstream effects of genetic variants on biochemical reactions involved in common human diseases.

Genome-Scale Metabolic Modeling Reveals Sequential Dysregulation of Glutathione Metabolism in Livers from Patients with Alcoholic Hepatitis

Alcoholic hepatitis (AH) is the most severe form of alcoholic liver disease for which there is no efficacious treatment aiding most patients. AH manifests differently in individuals, with some patients showing debilitating symptoms more so than others. Previous studies showed significant metabolic dysregulation associated with AH. Therefore, we sought to analyze how the activity of metabolic pathways differed in the liver of patients with varying degrees of AH severity. We utilized a genome-scale metabolic modeling approach that allowed for integration of a generic human cellular metabolic model with specific RNA-seq data corresponding to healthy and multiple liver disease states to predict the metabolic fluxes within each disease state. Additionally, we performed a systems-level analysis of the transcriptomic data and predicted metabolic flux data to identify the regulatory and functional differences in liver metabolism with increasing severity of AH. Our results provide unique insights into the sequential dysregulation of the solute transport mechanisms underlying the glutathione metabolic pathway with increasing AH disease severity. We propose targeting of the solute transporters in the glutathione pathway to mimic the flux activity of the healthy liver state as a potential therapeutic intervention for AH.

Combining denoising of RNA-seq data and flux balance analysis for cluster analysis of single cells

Background

Sophisticated methods to properly pre-process and analyze the increasing collection of single-cell RNA sequencing (scRNA-seq) data are increasingly being developed. On the contrary, the best practices to integrate these data into metabolic networks, aiming at describing metabolic phenotypes within a heterogeneous cell population, have been poorly investigated. In this regard, a critical factor is the presence of false zero values in reactions essential for a fundamental metabolic function, such as biomass or energy production. Here, we investigate the role of denoising strategies in mitigating this problem.

Methods

We applied state-of-the-art denoising strategies - namely MAGIC, ENHANCE, and SAVER - on three public scRNA-seq datasets. We then associated a metabolic flux distribution with every single cell by embedding its noise-free transcriptomics profile in the constraints of the optimization of a core metabolic model. Finally, we used the obtained single-cell optimal metabolic fluxes as features for cluster analysis. We compared the results obtained with different techniques, and with or without the use of denoising. We also investigated the possibility of applying denoising directly on the Reaction Activity Scores, which are metabolic features extracted from the read counts, rather than on the read counts.

Results

We show that denoising of transcriptomics data improves the clustering of single cells. We also illustrate that denoising restores important metabolic properties, such as the correlation between cell cycle phase and biomass accumulation, and between the RAS scores of reactions belonging to the same metabolic pathway. We show that MAGIC performs better than ENHANCE and SAVER, and that, denoising applied directly on the RAS matrix could be an effective alternative in removing false zero values from essential metabolic reactions.

Conclusions

Our results indicate that including denoising as a pre-processing operation represents a milestone to integrate scRNA-seq data into Flux Balance Analysis simulations and to perform single-cell cluster analysis with a focus on metabolic phenotypes.

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

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

Genome-scale metabolic models for natural and long-term drug-induced viral control in HIV infection

Genome-scale metabolic models (GSMMs) can provide novel insights into metabolic reprogramming during disease progression and therapeutic interventions. We developed a context-specific system-level GSMM of people living with HIV (PLWH) using global RNA sequencing data from PBMCs with suppressive viremia either by natural (elite controllers, PLWHEC) or drug-induced (PLWHART) control. This GSMM was compared with HIV-negative controls (HC) to provide a comprehensive systems-level metabo-transcriptomic characterization. Transcriptomic analysis identified up-regulation of oxidative phosphorylation as a characteristic of PLWHART, differentiating them from PLWHEC with dysregulated complexes I, III, and IV. The flux balance analysis identified altered flux in several intermediates of glycolysis including pyruvate, α-ketoglutarate, and glutamate, among others, in PLWHART The in vitro pharmacological inhibition of OXPHOS complexes in a latent lymphocytic cell model (J-Lat 10.6) suggested a role for complex IV in latency reversal and immunosenescence. Furthermore, inhibition of complexes I/III/IV induced apoptosis, collectively indicating their contribution to reservoir dynamics.

Multi-omics personalized network analyses highlight progressive disruption of central metabolism associated with COVID-19 severity

The clinical outcome and disease severity in coronavirus disease 2019 (COVID-19) are heterogeneous, and the progression or fatality of the disease cannot be explained by a single factor like age or comorbidities. In this study, we used system-wide network-based system biology analysis using whole blood RNA sequencing, immunophenotyping by flow cytometry, plasma metabolomics, and single-cell-type metabolomics of monocytes to identify the potential determinants of COVID-19 severity at personalized and group levels. Digital cell quantification and immunophenotyping of the mononuclear phagocytes indicated a substantial role in coordinating the immune cells that mediate COVID-19 severity. Stratum-specific and personalized genome-scale metabolic modeling indicated monocarboxylate transporter family genes (e.g., SLC16A6), nucleoside transporter genes (e.g., SLC29A1), and metabolites such as α-ketoglutarate, succinate, malate, and butyrate could play a crucial role in COVID-19 severity. Metabolic perturbations targeting the central metabolic pathway (TCA cycle) can be an alternate treatment strategy in severe COVID-19.

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.

ELIMINATOR: essentiality analysis using multisystem networks and integer programming

A gene is considered as essential when it is indispensable for cells to grow and replicate in a certain environment. However, gene essentiality is not a structural property but rather a contextual one, which depends on the specific biological conditions affecting the cell. This circumstantial essentiality of genes is what brings the attention of scientist since we can identify genes essential for cancer cells but not essential for healthy cells. This same contextuality makes their identification extremely challenging. Huge experimental efforts such as Project Achilles where the essentiality of thousands of genes is measured together with a plethora of molecular data (transcriptomics, copy number, mutations, etc.) in over one thousand cell lines can shed light on the causality behind the essentiality of a gene in a given environment. Here, we present an in-silico method for the identification of patient-specific essential genes using constraint-based modelling (CBM). Our method expands the ideas behind traditional CBM to accommodate multisystem networks. In essence, it first calculates the minimum number of lowly expressed genes required to be activated by the cell to sustain life as defined by a set of requirements; and second, it performs an exhaustive in-silico gene knockout to find those that lead to the need of activating additional lowly expressed genes. We validated the proposed methodology using a set of 452 cancer cell lines derived from the Cancer Cell Line Encyclopedia where an exhaustive experimental large-scale gene knockout study using CRISPR (Achilles Project) evaluates the impact of each removal. We also show that the integration of different essentiality predictions per gene, what we called Essentiality Congruity Score, reduces the number of false positives. Finally, we explored our method in a breast cancer patient dataset, and our results showed high concordance with previous publications. These findings suggest that identifying genes whose activity is fundamental to sustain cellular life in a patient-specific manner is feasible using in-silico methods. The patient-level gene essentiality predictions can pave the way for precision medicine by identifying potential drug targets whose deletion can induce death in tumour cells.

Genome-scale metabolic network models: from first-generation to next-generation

Over the last two decades, thousands of genome-scale metabolic network models (GSMMs) have been constructed. These GSMMs have been widely applied in various fields, ranging from network interaction analysis, to cell phenotype prediction. However, due to the lack of constraints, the prediction accuracy of first-generation GSMMs was limited. To overcome these limitations, the next-generation GSMMs were developed by integrating omics data, adding constrain condition, integrating different biological models, and constructing whole-cell models. Here, we review recent advances of GSMMs from the first generation to the next generation. Then, we discuss the major application of GSMMs in industrial biotechnology, such as predicting phenotypes and guiding metabolic engineering. In addition, human health applications, including understanding biological mechanisms, discovering biomarkers and drug targets, are also summarized. Finally, we address the challenges and propose new trend of GSMMs. KEY POINTS: •This mini-review updates the literature on almost all published GSMMs since 1999. •Detailed insights into the development of the first- and next-generation GSMMs. •The application of GSMMs is summarized, and the prospects of integrating machine learning are emphasized.

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.

Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into such models was first enabled by the GECKO toolbox, allowing the study of phenotypes constrained by protein limitations. Here, we upgrade the toolbox in order to enhance models with enzyme and proteomics constraints for any organism with a compatible GEM reconstruction. With this, enzyme-constrained models for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus are generated to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions reveal that upregulation and high saturation of enzymes in amino acid metabolism are common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO is expanded with an automated framework for continuous and version-controlled update of enzyme-constrained GEMs, also producing such models for Escherichia coli and Homo sapiens. In this work, we facilitate the utilization of enzyme-constrained GEMs in basic science, metabolic engineering and synthetic biology purposes.

A pipeline for the reconstruction and evaluation of context-specific human metabolic models at a large-scale

Constraint-based (CB) metabolic models provide a mathematical framework and scaffold for in silico cell metabolism analysis and manipulation. In the past decade, significant efforts have been done to model human metabolism, enabled by the increased availability of multi-omics datasets and curated genome-scale reconstructions, as well as the development of several algorithms for context-specific model (CSM) reconstruction. Although CSM reconstruction has revealed insights on the deregulated metabolism of several pathologies, the process of reconstructing representative models of human tissues still lacks benchmarks and appropriate integrated software frameworks, since many tools required for this process are still disperse across various software platforms, some of which are proprietary. In this work, we address this challenge by assembling a scalable CSM reconstruction pipeline capable of integrating transcriptomics data in CB models. We combined omics preprocessing methods inspired by previous efforts with in-house implementations of existing CSM algorithms and new model refinement and validation routines, all implemented in the Troppo Python-based open-source framework. The pipeline was validated with multi-omics datasets from the Cancer Cell Line Encyclopedia (CCLE), also including reference fluxomics measurements for the MCF7 cell line. We reconstructed over 6000 models based on the Human-GEM template model for 733 cell lines featured in the CCLE, using MCF7 models as reference to find the best parameter combinations. These reference models outperform earlier studies using the same template by comparing gene essentiality and fluxomics experiments. We also analysed the heterogeneity of breast cancer cell lines, identifying key changes in metabolism related to cancer aggressiveness. Despite the many challenges in CB modelling, we demonstrate using our pipeline that combining transcriptomics data in metabolic models can be used to investigate key metabolic shifts. Significant limitations were found on these models ability for reliable quantitative flux prediction, thus motivating further work in genome-wide phenotype prediction.

Genome-scale models of human metabolism are promising tools capable of contextualising large omics datasets within a framework that enables analysis and manipulation of metabolic phenotypes. Despite various successes in applying these methods to provide mechanistic hypotheses for deregulated metabolism in disease, there is no standardized workflow to extract these models using existing methods and the tools required to do so are mostly implemented using proprietary software. We have assembled a generic pipeline to extract and validate context-specific metabolic models using multi-omics datasets and implemented it using the troppo framework. We first validate our pipeline using MCF7 cell line models and assess their ability to predict lethal gene knockouts as well as flux activity using multi-omics data. We also demonstrate how this approach can be generalized for large-scale transcriptomics datasets and used to generate insights on the metabolic heterogeneity of cancer and relevant features for other data mining approaches. The pipeline is available as part of an open-source framework that is generic for a variety of applications.

Exploring the metabolic landscape of pancreatic ductal adenocarcinoma cells using genome-scale metabolic modeling

Pancreatic ductal adenocarcinoma (PDAC) is a major research focus because of its poor therapy response and dismal prognosis. PDAC cells adapt their metabolism to the surrounding environment, often relying on diverse nutrient sources. Because traditional experimental techniques appear exhaustive to find a viable therapeutic strategy, a highly curated and omics-informed PDAC genome-scale metabolic model was reconstructed using patient-specific transcriptomics data. From the model-predictions, several new metabolic functions were explored as potential therapeutic targets in addition to the known metabolic hallmarks of PDAC. Significant downregulation in the peroxisomal beta oxidation pathway, flux modulation in the carnitine shuttle system, and upregulation in the reactive oxygen species detoxification pathway reactions were observed. These unique metabolic traits of PDAC were correlated with potential drug combinations targeting genes with poor prognosis in PDAC. Overall, this study provides a better understanding of the metabolic vulnerabilities in PDAC and will lead to novel effective therapeutic strategies.

Integration of omics data to generate and analyse COVID-19 specific genome-scale metabolic models

COVID-19 presents a complex disease that needs to be addressed using systems medicine approaches that include genome-scale metabolic models (GEMs). Previous studies have used a single model extraction method (MEM) and/or a single transcriptomic dataset to reconstruct context-specific models, which proved to be insufficient for the broader biological contexts. We have applied four MEMs in combination with five COVID-19 datasets. Models produced by GIMME were separated by infection, while tINIT preserved the biological variability in the data and enabled the best prediction of the enrichment of metabolic subsystems. Vitamin D3 metabolism was predicted to be down-regulated in one dataset by GIMME, and in all by tINIT. Models generated by tINIT and GIMME predicted downregulation of retinol metabolism in different datasets, while downregulated cholesterol metabolism was predicted only by tINIT-generated models. Predictions are in line with the observations in COVID-19 patients. Our data indicated that GIMME and tINIT models provided the most biologically relevant results and should have a larger emphasis in further analyses. Particularly tINIT models identified the metabolic pathways that are a part of the host response and are potential antiviral targets. The code and the results of the analyses are available to download from https://github.com/CompBioLj/COVID_GEMs_and_MEMs.

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

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

Human/SARS-CoV-2 genome-scale metabolic modeling to discover potential antiviral targets for COVID-19

BACKGROUND

Coronavirus disease 2019 (COVID-19) has caused a substantial increase in mortality and economic and social disruption. The absence of US Food and Drug Administration-approved drugs for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) highlights the need for new therapeutic drugs to combat COVID-19.

METHODS

The present study proposed a fuzzy hierarchical optimization framework for identifying potential antiviral targets for COVID-19. The objectives in the decision-making problem were not only to evaluate the elimination of the virus growth, but also to minimize side effects causing treatment. The identified candidate targets could promote processes of drug discovery and development.

SIGNIFICANT FINDINGS

Our gene-centric method revealed that dihydroorotate dehydrogenase (DHODH) inhibition could reduce viral biomass growth and metabolic deviation by 99.4% and 65.6%, respectively, and increase cell viability by 70.4%. We also identified two-target combinations that could completely block viral biomass growth and more effectively prevent metabolic deviation. We also discovered that the inhibition of two antiviral metabolites, cytidine triphosphate (CTP) and uridine-5'-triphosphate (UTP), exhibits effects similar to those of molnupiravir, which is undergoing phase III clinical trials. Our predictions also indicate that CTP and UTP inhibition blocks viral RNA replication through a similar mechanism to that of molnupiravir.

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

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

A mechanistic modeling framework reveals the key principles underlying tumor metabolism

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

Computational Models for Clinical Applications in Personalized Medicine—Guidelines and Recommendations for Data Integration and Model Validation

The future development of personalized medicine depends on a vast exchange of data from different sources, as well as harmonized integrative analysis of large-scale clinical health and sample data. Computational-modelling approaches play a key role in the analysis of the underlying molecular processes and pathways that characterize human biology, but they also lead to a more profound understanding of the mechanisms and factors that drive diseases; hence, they allow personalized treatment strategies that are guided by central clinical questions. However, despite the growing popularity of computational-modelling approaches in different stakeholder communities, there are still many hurdles to overcome for their clinical routine implementation in the future. Especially the integration of heterogeneous data from multiple sources and types are challenging tasks that require clear guidelines that also have to comply with high ethical and legal standards. Here, we discuss the most relevant computational models for personalized medicine in detail that can be considered as best-practice guidelines for application in clinical care. We define specific challenges and provide applicable guidelines and recommendations for study design, data acquisition, and operation as well as for model validation and clinical translation and other research areas.

Metabolic modeling of host-microbe interactions for therapeutics in colorectal cancer

The onset of colorectal cancer (CRC) is often attributed to gut bacterial dysbiosis, and thus gut microbiota are highly relevant in devising treatment strategies. Certain gut microbes, like Enterococcus spp., exhibit remarkable anti-neoplastic and probiotic properties, which can aid in silver nanoparticle (AgNPs) induced reactive oxygen species (ROS)-based CRC treatment. However, the effects of AgNPs on gut microbial metabolism have not been reported thus far. In this study, a detailed systems-level understanding of ROS metabolism in Enterococcus durans (E. durans), a representative gut microbe, was gained using constraint-based modeling, wherein, the critical association between ROS and folate metabolism was established. Experimental studies involving low AgNP concentration treatment of E. durans cultures confirmed these modeling predictions (an increased extracellular folate concentration by 52%, at the 9^th h of microbial growth, was observed). Besides, the computational studies established various metabolic pathways involving amino acids, energy metabolites, nucleotides, and SCFAs as the key players in elevating folate levels on ROS exposure. The anti-cancer potential of E. durans was also studied through MTT analysis of HCT 116 cells treated with microbial culture (AgNP treated) supernatant. A decrease in cell viability by 19% implicated the role of microbial metabolites (primarily folate) in causing cell death. The genome-scale modeling approach was then extended to extensively model CRC metabolism, as well as CRC-E. durans interactions in the context of CRC treatment, using tissue-specific metabolic models of CRC and healthy colon. These findings on further validation can facilitate the development of robust and effective cancer therapy.

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

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

Host-mycobiome metabolic interactions in health and disease

Fungal communities (mycobiome) have an important role in sustaining the resilience of complex microbial communities and maintenance of homeostasis. The mycobiome remains relatively unexplored compared to the bacteriome despite increasing evidence highlighting their contribution to host-microbiome interactions in health and disease. Despite being a small proportion of the total species, fungi constitute a large proportion of the biomass within the human microbiome and thus serve as a potential target for metabolic reprogramming in pathogenesis and disease mechanism. Metabolites produced by fungi shape host niches, induce immune tolerance and changes in their levels prelude changes associated with metabolic diseases and cancer. Given the complexity of microbial interactions, studying the metabolic interplay of the mycobiome with both host and microbiome is a demanding but crucial task. However, genome-scale modelling and synthetic biology can provide an integrative platform that allows elucidation of the multifaceted interactions between mycobiome, microbiome and host. The inferences gained from understanding mycobiome interplay with other organisms can delineate the key role of the mycobiome in pathophysiology and reveal its role in human disease.

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

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

Systems Biology Approaches to Decipher the Underlying Molecular Mechanisms of Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is one of the most malignant central nervous system tumors, showing a poor prognosis and low survival rate. Therefore, deciphering the underlying molecular mechanisms involved in the progression of the GBM and identifying the key driver genes responsible for the disease progression is crucial for discovering potential diagnostic markers and therapeutic targets. In this context, access to various biological data, development of new methodologies, and generation of biological networks for the integration of multi-omics data are necessary for gaining insights into the appearance and progression of GBM. Systems biology approaches have become indispensable in analyzing heterogeneous high-throughput omics data, extracting essential information, and generating new hypotheses from biomedical data. This review provides current knowledge regarding GBM and discusses the multi-omics data and recent systems analysis in GBM to identify key biological functions and genes. This knowledge can be used to develop efficient diagnostic and treatment strategies and can also be used to achieve personalized medicine for GBM.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Personalized Immunotherapy in Colorectal Cancers: Where Do We Stand?

Colorectal cancer (CRC) is the second leading cause of cancer death in the world. Immunotherapy using monoclonal antibodies, immune-checkpoint inhibitors, adoptive cell therapy, and cancer vaccines has raised great hopes for treating poor prognosis metastatic CRCs that are resistant to the conventional therapies. However, high inter-tumor and intra-tumor heterogeneity hinder the success of immunotherapy in CRC. Patients with a similar tumor phenotype respond differently to the same immunotherapy regimen. Mutation-based classification, molecular subtyping, and immunoscoring of CRCs facilitated the multi-aspect grouping of CRC patients and improved immunotherapy. Personalized immunotherapy using tumor-specific neoantigens provides the opportunity to consider each patient as an independent group deserving of individualized immunotherapy. In the recent decade, the development of sequencing and multi-omics techniques has helped us classify patients more precisely. The expansion of such advanced techniques along with the neoantigen-based immunotherapy could herald a new era in treating heterogeneous tumors such as CRC. In this review article, we provided the latest findings in immunotherapy of CRC. We elaborated on the heterogeneity of CRC patients as a bottleneck of CRC immunotherapy and reviewed the latest advances in personalized immunotherapy to overcome CRC heterogeneity.

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

Combined Metabolic Activators Decrease Liver Steatosis by Activating Mitochondrial Metabolism in Hamsters Fed with a High-Fat Diet

Although the prevalence of non-alcoholic fatty liver disease (NAFLD) continues to increase, there is no effective treatment approved for this condition. We previously showed, in high-fat diet (HFD)-fed mice, that the supplementation of combined metabolic activators (CMA), including nicotinamide riboside (NAD+ precursor) and the potent glutathione precursors serine and N-acetyl-l-cysteine (NAC), significantly decreased fatty liver by promoting fat oxidation in mitochondria. Afterwards, in a one-day proof-of-concept human supplementation study, we observed that this CMA, including also L-carnitine tartrate (LCT), resulted in increased fatty acid oxidation and de novo glutathione synthesis. However, the underlying molecular mechanisms associated with supplementation of CMA have not been fully elucidated. Here, we demonstrated in hamsters that the chronic supplementation of this CMA (changing serine for betaine) at two doses significantly decreased hepatic steatosis. We further generated liver transcriptomics data and integrated these data using a liver-specific genome-scale metabolic model of liver tissue. We systemically determined the molecular changes after the supplementation of CMA and found that it activates mitochondria in the liver tissue by modulating global lipid, amino acid, antioxidant and folate metabolism. Our findings provide extra evidence about the beneficial effects of a treatment based on this CMA against NAFLD.

Systems Analysis Reveals Ageing-Related Perturbations in Retinoids and Sex Hormones in Alzheimer's and Parkinson's Diseases

Neurodegenerative diseases, including Alzheimer's (AD) and Parkinson's diseases (PD), are complex heterogeneous diseases with highly variable patient responses to treatment. Due to the growing evidence for ageing-related clinical and pathological commonalities between AD and PD, these diseases have recently been studied in tandem. In this study, we analysed transcriptomic data from AD and PD patients, and stratified these patients into three subclasses with distinct gene expression and metabolic profiles. Through integrating transcriptomic data with a genome-scale metabolic model and validating our findings by network exploration and co-analysis using a zebrafish ageing model, we identified retinoids as a key ageing-related feature in all subclasses of AD and PD. We also demonstrated that the dysregulation of androgen metabolism by three different independent mechanisms is a source of heterogeneity in AD and PD. Taken together, our work highlights the need for stratification of AD/PD patients and development of personalised and precision medicine approaches based on the detailed characterisation of these subclasses.

Path to improving the life cycle and quality of genome-scale models of metabolism

Genome-scale models of metabolism (GEMs) are key computational tools for the systems-level study of metabolic networks. Here, we describe the "GEM life cycle," which we subdivide into four stages: inception, maturation, specialization, and amalgamation. We show how different types of GEM reconstruction workflows fit in each stage and proceed to highlight two fundamental bottlenecks for GEM quality improvement: GEM maturation and content removal. We identify common characteristics contributing to increasing quality of maturing GEMs drawing from past independent GEM maturation efforts. We then shed some much-needed light on the latent and unrecognized but pervasive issue of content removal, demonstrating the substantial effects of model pruning on its solution space. Finally, we propose a novel framework for content removal and associated confidence-level assignment which will help guide future GEM development efforts, reduce duplication of effort across groups, potentially aid automated reconstruction platforms, and boost the reproducibility of model development.

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.

Maximization of non-idle enzymes improves the coverage of the estimated maximal in vivo enzyme catalytic rates in Escherichia coli

MOTIVATION

Constraint-based modeling approaches allow the estimation of maximal in vivo enzyme catalytic rates that can serve as proxies for enzyme turnover numbers. Yet, genome-scale flux profiling remains a challenge in deploying these approaches to catalogue proxies for enzyme catalytic rates across organisms.

RESULTS

Here we formulate a constraint-based approach, termed NIDLE-flux, to estimate fluxes at a genome-scale level by using the principle of efficient usage of expressed enzymes. Using proteomics data from Escherichia coli, we show that the fluxes estimated by NIDLE-flux and the existing approaches are in excellent qualitative agreement (Pearson correlation > 0.9). We also find that the maximal in vivo catalytic rates estimated by NIDLE-flux exhibits a Pearson correlation of 0.74 with in vitro enzyme turnover numbers. However, NIDLE-flux results in a 1.4-fold increase in the size of the estimated maximal in vivo catalytic rates in comparison to the contenders. Integration of the maximum in vivo catalytic rates with publically available proteomics and metabolomics data provide a better match to fluxes estimated by NIDLE-flux. Therefore, NIDLE-flux facilitates more effective usage of proteomics data to estimate proxies for kcatomes.

AVAILABILITY

https://github.com/Rudan-X/NIDLE-flux-code.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multi-Tissue Transcriptomic-Informed In Silico Investigation of Drugs for the Treatment of Dengue Fever Disease

Transcriptomics, proteomics and pathogen-host interactomics data are being explored for the in silico–informed selection of drugs, prior to their functional evaluation. The effectiveness of this kind of strategy has been put to the test in the current COVID-19 pandemic, and it has been paying off, leading to a few drugs being rapidly repurposed as treatment against SARS-CoV-2 infection. Several neglected tropical diseases, for which treatment remains unavailable, would benefit from informed in silico investigations of drugs, as performed in this work for Dengue fever disease. We analyzed transcriptomic data in the key tissues of liver, spleen and blood profiles and verified that despite transcriptomic differences due to tissue specialization, the common mechanisms of action, “Adrenergic receptor antagonist”, “ATPase inhibitor”, “NF-kB pathway inhibitor” and “Serotonin receptor antagonist”, were identified as druggable (e.g., oxprenolol, digoxin, auranofin and palonosetron, respectively) to oppose the effects of severe Dengue infection in these tissues. These are good candidates for future functional evaluation and clinical trials.

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.

Model-based assessment of mammalian cell metabolic functionalities using omics data

Omics experiments are ubiquitous in biological studies, leading to a deluge of data. However, it is still challenging to connect changes in these data to changes in cell functions because of complex interdependencies between genes, proteins, and metabolites. Here, we present a framework allowing researchers to infer how metabolic functions change on the basis of omics data. To enable this, we curated and standardized lists of metabolic tasks that mammalian cells can accomplish. Genome-scale metabolic networks were used to define gene sets associated with each metabolic task. We further developed a framework to overlay omics data on these sets and predict pathway usage for each metabolic task. We demonstrated how this approach can be used to quantify metabolic functions of diverse biological samples from the single cell to whole tissues and organs by using multiple transcriptomic datasets. To facilitate its adoption, we integrated the approach into GenePattern (www.genepattern.org-CellFie).

Genome-scale integration of transcriptome and metabolome unveils squalene synthase and dihydrofolate reductase as targets against AML cells resistant to chemotherapy

The development of resistance to chemotherapeutic agents, such as Doxorubicin (DOX) and cytarabine (AraC), is one of the greatest challenges to the successful treatment of Acute Myeloid Leukemia (AML). Such acquisition is often underlined by a metabolic reprogramming that can provide a therapeutic opportunity, as it can lead to the emergence of vulnerabilities and dependencies to be exploited as targets against the resistant cells. In this regard, genome-scale metabolic models (GSMMs) have emerged as powerful tools to integrate multiple layers of data to build cancer-specific models and identify putative metabolic vulnerabilities. Here, we use genome-scale metabolic modelling to reconstruct a GSMM of the THP1 AML cell line and two derivative cell lines, one with acquired resistance to AraC and the second with acquired resistance to DOX. We also explore how, adding to the transcriptomic layer, the metabolomic layer enhances the selectivity of the resulting condition specific reconstructions. The resulting models enabled us to identify and experimentally validate that drug-resistant THP1 cells are sensitive to the FDA-approved antifolate methotrexate. Moreover, we discovered and validated that the resistant cell lines could be selectively targeted by inhibiting squalene synthase, providing a new and promising strategy to directly inhibit cholesterol synthesis in AML drug resistant cells.

Circular RNAs Regulate Glucose Metabolism in Cancer Cells

Circular RNAs (circRNAs) were originally thought to result from RNA splicing errors. However, it has been shown that circRNAs can regulate cancer onset and progression in various ways. They can regulate cancer cell proliferation, differentiation, invasion, and metastasis. Moreover, they modulate glucose metabolism in cancer cells through different mechanisms such as directly regulating glycolytic enzymes and glucose transporter (GLUT) or indirectly regulating signal transduction pathways. In this review, we elucidate on the role of circRNAs in regulating glucose metabolism in cancer cells, which partly explains the pathogenesis of malignant tumors, and provides new therapeutic targets or new diagnostic and prognostic markers for human cancers.

Stratification of patients with clear cell renal cell carcinoma to facilitate drug repositioning

Clear cell renal cell carcinoma (ccRCC) is the most common histological type of kidney cancer and has high heterogeneity. Stratification of ccRCC is important since distinct subtypes differ in prognosis and treatment. Here, we applied a systems biology approach to stratify ccRCC into three molecular subtypes with different mRNA expression patterns and prognosis of patients. Further, we developed a set of biomarkers that could robustly classify the patients into each of the three subtypes and predict the prognosis of patients. Then, we reconstructed subtype-specific metabolic models and performed essential gene analysis to identify the potential drug targets. We identified four drug targets, including SOAT1, CRLS1, and ACACB, essential in all the three subtypes and GPD2, exclusively essential to subtype 1. Finally, we repositioned mitotane, an FDA-approved SOAT1 inhibitor, to treat ccRCC and showed that it decreased tumor cell viability and inhibited tumor cell growth based on in vitro experiments.

•

Three consistent molecular ccRCC subtypes were found to guide patients' prognoses

•

REOs-based biomarker was developed to robustly classify patients at individual level

•

SOAT1 is identified as a common drug target for all ccRCC subtypes

•

Mitotane was repositioned treatment of ccRCC via inhibiting SOAT1

Guided extraction of genome-scale metabolic models for the integration and analysis of omics data

Omics data can be integrated into a reference model using various model extraction methods (MEMs) to yield context-specific genome-scale metabolic models (GEMs). How to chose the appropriate MEM, thresholding rule and threshold remains a challenge. We integrated mouse transcriptomic data from a Cyp51 knockout mice diet experiment (GSE58271) using five MEMs (GIMME, iMAT, FASTCORE, INIT an tINIT) in a combination with a recently published mouse GEM iMM1865. Except for INIT and tINIT, the size of extracted models varied with the MEM used (t-test: p-value < 0.001). The Jaccard index of iMAT models ranged from 0.27 to 1.0. Out of the three factors under study in the experiment (diet, gender and genotype), gender explained most of the variability ( > 90%) in PC1 for FASTCORE. In iMAT, each of the three factors explained less than 40% of the variability within PC1, PC2 and PC3. Among all the MEMs, FASTCORE captured the most of the true variability in the data by clustering samples by gender. Our results show that for the efficient use of MEMs in the context of omics data integration and analysis, one should apply various MEMs, thresholding rules, and thresholding values to select the MEM and its configuration that best captures the true variability in the data. This selection can be guided by the methodology as proposed and used in this paper. Moreover, we describe certain approaches that can be used to analyse the results obtained with the selected MEM and to put these results in a biological context.

Simultaneous Integration of Gene Expression and Nutrient Availability for Studying the Metabolism of Hepatocellular Carcinoma Cell Lines

How cancer cells utilize nutrients to support their growth and proliferation in complex nutritional systems is still an open question. However, it is certainly determined by both genetics and an environmental-specific context. The interactions between them lead to profound metabolic specialization, such as consuming glucose and glutamine and producing lactate at prodigious rates. To investigate whether and how glucose and glutamine availability impact metabolic specialization, we integrated computational modeling on the genome-scale metabolic reconstruction with an experimental study on cell lines. We used the most comprehensive human metabolic network model to date, Recon3D, to build cell line-specific models. RNA-Seq data was used to specify the activity of genes in each cell line and the uptake rates were quantitatively constrained according to nutrient availability. To integrated both constraints we applied a novel method, named Gene Expression and Nutrients Simultaneous Integration (GENSI), that translates the relative importance of gene expression and nutrient availability data into the metabolic fluxes based on an observed experimental feature(s). We applied GENSI to study hepatocellular carcinoma addiction to glucose/glutamine. We were able to identify that proliferation, and lactate production is associated with the presence of glucose but does not necessarily increase with its concentration when the latter exceeds the physiological concentration. There was no such association with glutamine. We show that the integration of gene expression and nutrient availability data into genome-wide models improves the prediction of metabolic phenotypes.

Building patient-specific models for receptor tyrosine kinase signaling networks

Cancer progresses due to changes in the dynamic interactions of multidimensional factors associated with gene mutations. Cancer research has actively adopted computational methods, including data-driven and mathematical model-driven approaches, to identify causative factors and regulatory rules that can explain the complexity and diversity of cancers. A data-driven, statistics-based approach revealed correlations between gene alterations and clinical outcomes in many types of cancers. A model-driven mathematical approach has elucidated the dynamic features of cancer networks and identified the mechanisms of drug efficacy and resistance. More recently, machine learning methods have emerged that can be used for mining omics data and classifying patient. However, as the strengths and weaknesses of each method becoming apparent, new analytical tools are emerging to combine and improve the methodologies and maximize their predictive power for classifying cancer subtypes and prognosis. Here, we introduce recent advances in cancer systems biology aimed at personalized medicine, with focus on the receptor tyrosine kinase signaling network.

Metabolomics, metabolic flux analysis and cancer pharmacology

Metabolic reprogramming is a hallmark of cancer and increasing evidence suggests that reprogrammed cell metabolism supports tumor initiation, progression, metastasis and drug resistance. Understanding metabolic dysregulation may provide therapeutic targets and facilitate drug research and development for cancer therapy. Metabolomics enables the high-throughput characterization of a large scale of small molecule metabolites in cells, tissues and biofluids, while metabolic flux analysis (MFA) tracks dynamic metabolic activities using stable isotope tracer methods. Recent advances in metabolomics and MFA technologies make them powerful tools for metabolic profiling and characterizing metabolic activities in health and disease, especially in cancer research. In this review, we introduce recent advances in metabolomics and MFA analytical technologies, and provide the first comprehensive summary of the most commonly used isotope tracing methods. In addition, we highlight how metabolomics and MFA are applied in cancer pharmacology studies particularly for discovering targetable metabolic vulnerabilities, understanding the mechanisms of drug action and drug resistance, exploring potential strategies with dietary intervention, identifying cancer biomarkers, as well as enabling precision treatment with pharmacometabolomics.

Identifying functional metabolic shifts in heart failure with the integration of omics data and a heart-specific, genome-scale model

In diseased states, the heart can shift to use different carbon substrates, measured through changes in uptake of metabolites by imaging methods or blood metabolomics. However, it is not known whether these measured changes are a result of transcriptional changes or external factors. Here, we explore transcriptional changes in late-stage heart failure using publicly available data integrated with a model of heart metabolism. First, we present a heart-specific genome-scale metabolic network reconstruction (GENRE), iCardio. Next, we demonstrate the utility of iCardio in interpreting heart failure gene expression data by identifying tasks inferred from differential expression (TIDEs), which represent metabolic functions associated with changes in gene expression. We identify decreased gene expression for nitric oxide (NO) and N-acetylneuraminic acid (Neu5Ac) synthesis as common metabolic markers of heart failure. The methods presented here for constructing a tissue-specific model and identifying TIDEs can be extended to multiple tissues and diseases of interest.

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

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