Model-driven discovery of long-chain fatty acid metabolic reprogramming in heterogeneous prostate cancer cells

Applications of genome-scale metabolic models to the study of human diseases: A systematic review

BACKGROUND AND OBJECTIVES

Genome-scale metabolic networks (GEMs) represent a valuable modeling and computational tool in the broad field of systems biology. Their ability to integrate constraints and high-throughput biological data enables the study of intricate metabolic aspects and processes of different cell types and conditions. The past decade has witnessed an increasing number and variety of applications of GEMs for the study of human diseases, along with a huge effort aimed at the reconstruction, integration and analysis of a high number of organisms. This paper presents a systematic review of the scientific literature, to pursue several important questions about the application of constraint-based modeling in the investigation of human diseases. Hopefully, this paper will provide a useful reference for researchers interested in the application of modeling and computational tools for the investigation of metabolic-related human diseases.

METHODS

This systematic review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Elsevier Scopus®, National Library of Medicine PubMed® and Clarivate Web of Science™ databases were enquired, resulting in 566 scientific articles. After applying exclusion and eligibility criteria, a total of 169 papers were selected and individually examined.

RESULTS

The reviewed papers offer a thorough and up-to-date picture of the latest modeling and computational approaches, based on genome-scale metabolic models, that can be leveraged for the investigation of a large variety of human diseases. The numerous studies have been categorized according to the clinical research area involved in the examined disease. Furthermore, the paper discusses the most typical approaches employed to derive clinically-relevant information using the computational models.

CONCLUSIONS

The number of scientific papers, utilizing GEM-based approaches for the investigation of human diseases, suggests an increasing interest in these types of approaches; hopefully, the present review will represent a useful reference for scientists interested in applying computational modeling approaches to investigate the aetiopathology of human diseases; we also hope that this work will foster the development of novel applications and methods for the discovery of clinically-relevant insights on metabolic-related diseases.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

A Multi-Level Systems Biology Analysis of Aldrin's Metabolic Effects on Prostate Cancer Cells

Although numerous studies support a dose-effect relationship between Endocrine disruptors (EDs) and the progression and malignancy of tumors, the impact of a chronic exposure to non-lethal concentrations of EDs in cancer remains unknown. More specifically, a number of studies have reported the impact of Aldrin on a variety of cancer types, including prostate cancer. In previous studies, we demonstrated the induction of the malignant phenotype in DU145 prostate cancer (PCa) cells after a chronic exposure to Aldrin (an ED). Proteins are pivotal in the regulation and control of a variety of cellular processes. However, the mechanisms responsible for the impact of ED on PCa and the role of proteins in this process are not yet well understood. Here, two complementary computational approaches have been employed to investigate the molecular processes underlying the acquisition of malignancy in prostate cancer. First, the metabolic reprogramming associated with the chronic exposure to Aldrin in DU145 cells was studied by integrating transcriptomics and metabolomics via constraint-based metabolic modeling. Second, gene set enrichment analysis was applied to determine (i) altered regulatory pathways and (ii) the correlation between changes in the transcriptomic profile of Aldrin-exposed cells and tumor progression in various types of cancer. Experimental validation confirmed predictions revealing a disruption in metabolic and regulatory pathways. This alteration results in the modification of protein levels crucial in regulating triacylglyceride/cholesterol, linked to the malignant phenotype observed in Aldrin-exposed cells.

A Dual Role for FADD in Human Precursor T-Cell Neoplasms

A reduction in FADD levels has been reported in precursor T-cell neoplasms and other tumor types. Such reduction would impact on the ability of tumor cells to undergo apoptosis and has been associated with poor clinical outcomes. However, FADD is also known to participate in non-apoptotic functions, but these mechanisms are not well-understood. Linking FADD expression to the severity of precursor T-cell neoplasms could indicate its use as a prognostic marker and may open new avenues for targeted therapeutic strategies. Using transcriptomic and clinical data from patients with precursor T-cell neoplasms, complemented by in vitro analysis of cellular functions and by high-throughput interactomics, our results allow us to propose a dual role for FADD in precursor T-cell neoplasms, whereby resisting cell death and chemotherapy would be a canonical consequence of FADD deficiency in these tumors, whereas deregulation of the cellular metabolism would be a relevant non-canonical function in patients expressing FADD. These results reveal that evaluation of FADD expression in precursor T-cell neoplasms may aid in the understanding of the biological processes that are affected in the tumor cells. The altered biological processes can be of different natures depending on the availability of FADD influencing its ability to exert its canonical or non-canonical functions. Accordingly, specific therapeutic interventions would be needed in each case.

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.

Untargeted metabolomics reveals alterations in the metabolic reprogramming of prostate cancer cells by double-stranded DNA-modified gold nanoparticles

Metabolic reprogramming plays an important role in the development of prostate cancer (PCa). However, there are few reports on the effects of nanomaterials as vectors on cancer metabolic reprogramming. Herein, a type of nanoparticle with good biocompatibility was synthesized by modifying the double-stranded of DNA containing a sulfhydryl group on the surface of gold nanoparticles (AuNPs-dsDNA) through salt-aging conjugation methods. The resultant AuNPs-dsDNA complexes possessed low toxicity to PC3 and DU145 cells in vitro. There was also no obvious hepatorenal toxicity after intravenous injection of AuNPs-dsDNA complexes in vivo, which indicated that these nanoparticles had good biological compatibilities. We investigated their biological functions using prostate cancer cells. Seahorse assay showed that AuNPs-dsDNA complexes could increase glycolysis and glycolysis capacity both in PC3 and DU145 cells. We further detected the expression of glycolysis-related genes by qPCR assay, and found that PKM2, PDHA, and LDHA were significantly upregulated. Furthermore, untargeted metabolomics revealed that PC (18:2(9Z,12Z)/18:2(9Z,12Z)) and PC (18:0/18:2 (9Z,12Z)) levels were decreased and inosinic acid level was increased in PC3 cells. Whereas (3S,6E,10E)-1,6,10,14-Phytatetraen-3-ol, Plasmenyl-PE 36:5 and Cer (d18:2/18:2) were decreased, PE 21:3 and 1-pyrrolidinecarboxaldehyde were increased in DU145 cells after co-culturing with AuNPs-dsDNA. In summary, we found that AuNPs and AuNPs-dsDNA complexes possibly regulate the metabolic reprogramming of cancer cells mainly through the lipid metabolic pathways, which could compensate for the previously mentioned phenomenon of enhanced glycolysis and glycolysis capacity. This will provide an important theoretical basis for our future research on the characteristic targeted design of nanomaterials for cancer metabolism.

[18F]FMCH PET/CT biomarkers and similarity analysis to refine the definition of oligometastatic prostate cancer

Background

The role of image-derived biomarkers in recurrent oligometastatic Prostate Cancer (PCa) is unexplored. This paper aimed to evaluate [¹⁸F]FMCH PET/CT radiomic analysis in patients with recurrent PCa after primary radical therapy. Specifically, we tested intra-patient lesions similarity in oligometastatic and plurimetastatic PCa, comparing the two most used definitions of oligometastatic disease.

Methods

PCa patients eligible for [¹⁸F]FMCH PET/CT presenting biochemical failure after first-line curative treatments were invited to participate in this prospective observational trial. PET/CT images of 92 patients were visually and quantitatively analyzed. Each patient was classified as oligometastatic or plurimetastatic according to the total number of detected lesions (up to 3 and up to 5 or > 3 and > 5, respectively). Univariate and intra-patient lesions' similarity analysis were performed.

Results

[¹⁸F]FMCH PET/CT identified 370 lesions, anatomically classified as regional lymph nodes and distant metastases. Thirty-eight and 54 patients were designed oligometastatic and plurimetastatic, respectively, using a 3-lesion threshold. The number of oligometastic scaled up to 60 patients (thus 32 plurimetastatic patients) with a 5-lesion threshold. Similarity analysis showed high lesions' heterogeneity. Grouping patients according to the number of metastases, patients with oligometastatic PCa defined with a 5-lesion threshold presented lesions heterogeneity comparable to plurimetastic patients. Lesions within patients having a limited tumor burden as defined by three lesions were characterized by less heterogeneity.

Conclusions

We found a comparable heterogeneity between patients with up to five lesions and plurimetastic patients, while patients with up to three lesions were less heterogeneous than plurimetastatic patients, featuring different cells phenotypes in the two groups. Our results supported the use of a 3-lesion threshold to define oligometastatic PCa.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13550-021-00858-8.

Agent Clustering Strategy Based on Metabolic Flux Distribution and Transcriptome Expression for Novel Drug Development

The network module-based method has been used for drug repositioning. The traditional drug repositioning method only uses the gene characteristics of the drug but ignores the drug-triggered metabolic changes. The metabolic network systematically characterizes the connection between genes, proteins, and metabolic reactions. The differential metabolic flux distribution, as drug metabolism characteristics, was employed to cluster the agents with similar MoAs (mechanism of action). In this study, agents with the same pharmacology were clustered into one group, and a total of 1309 agents from the CMap database were clustered into 98 groups based on differential metabolic flux distribution. Transcription factor (TF) enrichment analysis revealed the agents in the same group (such as group 7 and group 26) were confirmed to have similar MoAs. Through this agent clustering strategy, the candidate drugs which can inhibit (Japanese encephalitis virus) JEV infection were identified. This study provides new insights into drug repositioning and their MoAs.

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.

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.

Cysteine and Folate Metabolism Are Targetable Vulnerabilities of Metastatic Colorectal Cancer

Simple Summary

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

Abstract

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

Revisiting prostate cancer metabolism: From metabolites to disease and therapy

Prostate cancer (PCa), one of the most commonly diagnosed cancers worldwide, still presents important unmet clinical needs concerning treatment. In the last years, the metabolic reprogramming and the specificities of tumor cells emerged as an exciting field for cancer therapy. The unique features of PCa cells metabolism, and the activation of specific metabolic pathways, propelled the use of metabolic inhibitors for treatment. The present work revises the knowledge of PCa metabolism and the metabolic alterations that underlie the development and progression of the disease. A focus is given to the role of bioenergetic sources, namely, glucose, lipids, and glutamine sustaining PCa cell survival and growth. Moreover, it is described as the action of oncogenes/tumor suppressors and sex steroid hormones in the metabolic reprogramming of PCa. Finally, the status of PCa treatment based on the inhibition of metabolic pathways is presented. Globally, this review updates the landscape of PCa metabolism, highlighting the critical metabolic alterations that could have a clinical and therapeutic interest.

Low Entropy Sub-Networks Prevent the Integration of Metabolomic and Transcriptomic Data

The constantly and rapidly increasing amount of the biological data gained from many different high-throughput experiments opens up new possibilities for data- and model-driven inference. Yet, alongside, emerges a problem of risks related to data integration techniques. The latter are not so widely taken account of. Especially, the approaches based on the flux balance analysis (FBA) are sensitive to the structure of a metabolic network for which the low-entropy clusters can prevent the inference from the activity of the metabolic reactions. In the following article, we set forth problems that may arise during the integration of metabolomic data with gene expression datasets. We analyze common pitfalls, provide their possible solutions, and exemplify them by a case study of the renal cell carcinoma (RCC). Using the proposed approach we provide a metabolic description of the known morphological RCC subtypes and suggest a possible existence of the poor-prognosis cluster of patients, which are commonly characterized by the low activity of the drug transporting enzymes crucial in the chemotherapy. This discovery suits and extends the already known poor-prognosis characteristics of RCC. Finally, the goal of this work is also to point out the problem that arises from the integration of high-throughput data with the inherently nonuniform, manually curated low-throughput data. In such cases, the over-represented information may potentially overshadow the non-trivial discoveries.

Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.

Stoichiometric gene-to-reaction associations enhance model-driven analysis performance: Metabolic response to chronic exposure to Aldrin in prostate cancer

BACKGROUND

Genome-scale metabolic models (GSMM) integrating transcriptomics have been widely used to study cancer metabolism. This integration is achieved through logical rules that describe the association between genes, proteins, and reactions (GPRs). However, current gene-to-reaction formulation lacks the stoichiometry describing the transcript copies necessary to generate an active catalytic unit, which limits our understanding of how genes modulate metabolism. The present work introduces a new state-of-the-art GPR formulation that considers the stoichiometry of the transcripts (S-GPR). As case of concept, this novel gene-to-reaction formulation was applied to investigate the metabolic effects of the chronic exposure to Aldrin, an endocrine disruptor, on DU145 prostate cancer cells. To this aim we integrated the transcriptomic data from Aldrin-exposed and non-exposed DU145 cells through S-GPR or GPR into a human GSMM by applying different constraint-based-methods.

RESULTS

Our study revealed a significant improvement of metabolite consumption/production predictions when S-GPRs are implemented. Furthermore, our computational analysis unveiled important alterations in carnitine shuttle and prostaglandine biosynthesis in Aldrin-exposed DU145 cells that is supported by bibliographic evidences of enhanced malignant phenotype.

CONCLUSIONS

The method developed in this work enables a more accurate integration of gene expression data into model-driven methods. Thus, the presented approach is conceptually new and paves the way for more in-depth studies of aberrant cancer metabolism and other diseases with strong metabolic component with important environmental and clinical implications.

Current status and applications of genome-scale metabolic models

Genome-scale metabolic models (GEMs) computationally describe gene-protein-reaction associations for entire metabolic genes in an organism, and can be simulated to predict metabolic fluxes for various systems-level metabolic studies. Since the first GEM for Haemophilus influenzae was reported in 1999, advances have been made to develop and simulate GEMs for an increasing number of organisms across bacteria, archaea, and eukarya. Here, we review current reconstructed GEMs and discuss their applications, including strain development for chemicals and materials production, drug targeting in pathogens, prediction of enzyme functions, pan-reactome analysis, modeling interactions among multiple cells or organisms, and understanding human diseases.

Discovery of therapeutic agents for prostate cancer using genome-scale metabolic modeling and drug repositioning

BACKGROUND

Genome-scale metabolic models (GEMs) offer insights into cancer metabolism and have been used to identify potential biomarkers and drug targets. Drug repositioning is a time- and cost-effective method of drug discovery that can be applied together with GEMs for effective cancer treatment.

METHODS

In this study, we reconstruct a prostate cancer (PRAD)-specific GEM for exploring prostate cancer metabolism and also repurposing new therapeutic agents that can be used in development of effective cancer treatment. We integrate global gene expression profiling of cell lines with >1000 different drugs through the use of prostate cancer GEM and predict possible drug-gene interactions.

FINDINGS

We identify the key reactions with altered fluxes based on the gene expression changes and predict the potential drug effect in prostate cancer treatment. We find that sulfamethoxypyridazine, azlocillin, hydroflumethiazide, and ifenprodil can be repurposed for the treatment of prostate cancer based on an in silico cell viability assay. Finally, we validate the effect of ifenprodil using an in vitro cell assay and show its inhibitory effect on a prostate cancer cell line.

INTERPRETATION

Our approach demonstate how GEMs can be used to predict therapeutic agents for cancer treatment based on drug repositioning. Besides, it paved a way and shed a light on the applicability of computational models to real-world biomedical or pharmaceutical problems.

Drug Repositioning for Effective Prostate Cancer Treatment

Drug repositioning has gained attention from both academia and pharmaceutical companies as an auxiliary process to conventional drug discovery. Chemotherapeutic agents have notorious adverse effects that drastically reduce the life quality of cancer patients so drug repositioning is a promising strategy to identify non-cancer drugs which have anti-cancer activity as well as tolerable adverse effects for human health. There are various strategies for discovery and validation of repurposed drugs. In this review, 25 repurposed drug candidates are presented as result of different strategies, 15 of which are already under clinical investigation for treatment of prostate cancer (PCa). To date, zoledronic acid is the only repurposed, clinically used, and approved non-cancer drug for PCa. Anti-cancer activities of existing drugs presented in this review cover diverse and also known mechanisms such as inhibition of mTOR and VEGFR2 signaling, inhibition of PI3K/Akt signaling, COX and selective COX-2 inhibition, NF-κB inhibition, Wnt/β-Catenin pathway inhibition, DNMT1 inhibition, and GSK-3β inhibition. In addition to monotherapy option, combination therapy with current anti-cancer drugs may also increase drug efficacy and reduce adverse effects. Thus, drug repositioning may become a key approach for drug discovery in terms of time- and cost-efficiency comparing to conventional drug discovery and development process.

Development Of A Three-Gene Prognostic Signature For Hepatitis B Virus Associated Hepatocellular Carcinoma Based On Integrated Transcriptomic Analysis

Integration of public genome-wide gene expression data together with Cox regression analysis is a powerful weapon to identify new prognostic gene signatures for cancer diagnosis and prognosis. Hepatitis B virus (HBV) is a major cause of hepatocellular carcinoma (HCC), however, it remains largely unknown about the specific gene prognostic signature of HBV-associated HCC. Using Robust Rank Aggreg (RRA) method to integrate seven whole genome expression datasets, we identified 82 up-regulated genes and 577 down-regulated genes in HBV-associated HCC patients. Combination of several enrichment analysis, univariate and multivariate Cox proportional hazards regression analysis, we revealed that a three-gene (SPP2, CDC37L1, and ECHDC2) prognostic signature could act as an independent prognostic indicator for HBV-associated HCC in both the discovery cohort and the internal testing cohort. Gene set enrichment analysis showed that the high-risk group with lower expression levels of the three genes was enriched in bladder cancer and cell cycle pathway, whereas the low-risk group with higher expression levels of the three genes was enriched in drug metabolism-cytochrome P450, PPAR signaling pathway, fatty acid and histidine metabolisms. This indicates that patients of HBV-associated HCC with higher expression of these three genes may preserve relatively good hepatic cellular metabolism and function, which may also protect HCC patients from persistent drug toxicity in response to various medication. Our findings suggest a three-gene prognostic model that serves as a specific prognostic signature for HBV-associated HCC.