Genome-Scale Metabolic Modeling of Glioblastoma Reveals Promising Targets for Drug Development

From Data to Cure: A Comprehensive Exploration of Multi-omics Data Analysis for Targeted Therapies

In the dynamic landscape of targeted therapeutics, drug discovery has pivoted towards understanding underlying disease mechanisms, placing a strong emphasis on molecular perturbations and target identification. This paradigm shift, crucial for drug discovery, is underpinned by big data, a transformative force in the current era. Omics data, characterized by its heterogeneity and enormity, has ushered biological and biomedical research into the big data domain. Acknowledging the significance of integrating diverse omics data strata, known as multi-omics studies, researchers delve into the intricate interrelationships among various omics layers. This review navigates the expansive omics landscape, showcasing tailored assays for each molecular layer through genomes to metabolomes. The sheer volume of data generated necessitates sophisticated informatics techniques, with machine-learning (ML) algorithms emerging as robust tools. These datasets not only refine disease classification but also enhance diagnostics and foster the development of targeted therapeutic strategies. Through the integration of high-throughput data, the review focuses on targeting and modeling multiple disease-regulated networks, validating interactions with multiple targets, and enhancing therapeutic potential using network pharmacology approaches. Ultimately, this exploration aims to illuminate the transformative impact of multi-omics in the big data era, shaping the future of biological research.

Targeting metabolic fluxes reverts metastatic transitions in ovarian cancer

The formation of spheroids during epithelial ovarian cancer progression is correlated with peritoneal metastasis, disease recurrence, and poor prognosis. Although metastasis has been demonstrated to be driven by metabolic changes in transformed cells, mechanistic associations between metabolism and phenotypic transitions remain ill-explored. We performed quantitative proteomics to identify protein signatures associated with three distinct phenotypic morphologies (2D monolayers and two geometrically distinct three-dimensional spheroidal states) of the high-grade serous ovarian cancer line OVCAR-3. We obtained disease-driving phenotype-specific metabolic reaction modules and elucidated gene knockout strategies to reduce metabolic alterations that could drive phenotypic transitions. Exploring the DrugBank database, we identified and evaluated drugs that could impair such transitions and, hence, cancer progression. Finally, we experimentally validated our predictions by confirming the ability of one of our predicted drugs, the neuraminidase inhibitor oseltamivir, to inhibit spheroidogenesis in three ovarian cancer cell lines without any cytotoxic effects on untransformed stromal mesothelia.

Novel sights on therapeutic, prognostic, and diagnostics aspects of non-coding RNAs in glioblastoma multiforme

Glioblastoma Multiforme (GBM) is the primary brain tumor and accounts for 200,000 deaths each year worldwide. The standard therapy includes surgical resection followed by temozolomide (TMZ)-based chemotherapy and radiotherapy. The survival period of GBM patients is only 12-15 months. Therefore, novel treatment modalities for GBM treatment are urgently needed. Mounting evidence reveals that non-coding RNAs (ncRNAs) were involved in regulating gene expression, the pathophysiology of GBM, and enhancing therapeutic outcomes. The combinatory use of ncRNAs, chemotherapeutic drugs, and tumor suppressor gene expression induction might provide an innovative, alternative therapeutic approach for managing GBM. Studies have highlighted the role of Long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) in prognosis and diagnosis. Dysregulation of ncRNAs is observed in virtually all tumor types, including GBMs. Studies have also indicated the blood-brain barrier (BBB) as a crucial factor that hinders chemotherapy. Although several nanoparticle-mediated drug deliveries were degrading effectively against GBM in vitro conditions. However, the potential to cross the BBB and optimum delivery of oligonucleotide RNA into GBM cells in the brain is currently under intense clinical trials. Despite several advances in molecular pathogenesis, GBM remains resistant to chemo and radiotherapy. Targeted therapies have less clinical benefit due to high genetic heterogeneity and activation of alternative pathways. Thus, identifying GBM-specific prognostic pathways, essential genes, and genomic aberrations provide several potential benefits as subtypes of GBM. Also, these approaches will provide insights into new strategies to overcome the heterogenous nature of GBM, which will eventually lead to successful therapeutic interventions toward precision medicine and precision oncology.

Genome-scale metabolic modeling of Aspergillus fumigatus strains reveals growth dependencies on the lung microbiome

Aspergillus fumigatus, an opportunistic human pathogen, frequently infects the lungs of people with cystic fibrosis and is one of the most common causes of infectious-disease death in immunocompromised patients. Here, we construct 252 strain-specific, genome-scale metabolic models of this important fungal pathogen to study and better understand the metabolic component of its pathogenic versatility. The models show that 23.1% of A. fumigatus metabolic reactions are not conserved across strains and are mainly associated with amino acid, nucleotide, and nitrogen metabolism. Profiles of non-conserved reactions and growth-supporting reaction fluxes are sufficient to differentiate strains, for example by environmental or clinical origin. In addition, shotgun metagenomics analysis of sputum from 40 cystic fibrosis patients (15 females, 25 males) before and after diagnosis with an A. fumigatus colonization suggests that the fungus shapes the lung microbiome towards a more beneficial fungal growth environment associated with aromatic amino acid availability and the shikimate pathway. Our findings are starting points for the development of drugs or microbiome intervention strategies targeting fungal metabolic needs for survival and colonization in the non-native environment of the human lung.

Identification of Therapeutic Targets for Medulloblastoma by Tissue-Specific Genome-Scale Metabolic Model

Medulloblastoma (MB), occurring in the cerebellum, is the most common childhood brain tumor. Because conventional methods decline life quality and endanger children with detrimental side effects, computer models are needed to imitate the characteristics of cancer cells and uncover effective therapeutic targets with minimum toxic effects on healthy cells. In this study, metabolic changes specific to MB were captured by the genome-scale metabolic brain model integrated with transcriptome data. To determine the roles of sphingolipid metabolism in proliferation and metastasis in the cancer cell, 79 reactions were incorporated into the MB model. The pathways employed by MB without a carbon source and the link between metastasis and the Warburg effect were examined in detail. To reveal therapeutic targets for MB, biomass-coupled reactions, the essential genes/gene products, and the antimetabolites, which might deplete the use of metabolites in cells by triggering competitive inhibition, were determined. As a result, interfering with the enzymes associated with fatty acid synthesis (FAs) and the mevalonate pathway in cholesterol synthesis, suppressing cardiolipin production, and tumor-supporting sphingolipid metabolites might be effective therapeutic approaches for MB. Moreover, decreasing the activity of succinate synthesis and GABA-catalyzing enzymes concurrently might be a promising strategy for metastatic MB.

Whole-cell energy modeling reveals quantitative changes of predicted energy flows in RAS mutant cancer cell lines

Cellular utilization of available energy flows to drive a multitude of forms of cellular "work" is a major biological constraint. Cells steer metabolism to address changing phenotypic states but little is known as to how bioenergetics couples to the richness of processes in a cell as a whole. Here, we outline a whole-cell energy framework that is informed by proteomic analysis and an energetics-based gene ontology. We separate analysis of metabolic supply and the capacity to generate high-energy phosphates from a representation of demand that is built on the relative abundance of ATPases and GTPases that deliver cellular work. We employed mouse embryonic fibroblast cell lines that express wild-type KRAS or oncogenic mutations and with distinct phenotypes. We observe shifts between energy-requiring processes. Calibrating against Seahorse analysis, we have created a whole-cell energy budget with apparent predictive power, for instance in relation to protein synthesis.

A Risk Signature Consisting of Eight m6A Methylation Regulators Predicts the Prognosis of Glioma

Glioma progression seriously correlates to the epigenetic context. This study aims to identify glioma subtypes by clustering analysis of patients using the multi-omics data of N6-methyladenosine (m⁶A) methylation regulators and to construct a risk signature for investigating the role of m⁶A methylation regulators in the prognosis of glioma. Multi-omics data of glioma and normal control tissues were obtained through The Cancer Genome Atlas (TCGA) database. The clustering analysis of multi-omics data of patients was conducted using the R package iClusterPlus software. The risk model was constructed by univariate and multivariate Cox analysis, and the glioma expression data and related clinical data were obtained by Chinese Glioma Genome Atlas (CGGA) datasets to verify the risk model. By analyzing the glioma data in TCGA, we found that the risk signature could be constructed according to the eight genes with m⁶A methylation modification function, including ALKBH5, HNRNPA2B1, IGF2BP2, IGF2BP3, RBM15, WTAP, YTHDF1, and YTHDF2. Meanwhile, we found that IGF2BP2 and IGF2BP3 were highly expressed in glioma subtypes with high-risk scores and closely related to the prognosis of glioma patients. m⁶A methylation regulators, especially IGF2BP2 and IGF2BP3, play important roles in the malignant progression of glioma. The risk signature constructed by eight m⁶A methylation regulators can predict the prognosis of glioma. IGF2BP2 and IGF2BP3 may be the key regulatory factors of m⁶A methylation regulators involved in the occurrence and development of glioma, and can serve as molecular markers for the prognosis of glioma.

In silico analysis of metabolic effects of bipolar disorder on prefrontal cortex identified altered GABA, glutamate-glutamine cycle, energy metabolism and amino acid synthesis pathways

Bipolar disorder (BP) is a lifelong psychiatric condition, which often disrupts the daily life of the patients. It is characterized by unstable and periodic mood changes, which cause patients to display unusual shifts in mood, energy, activity levels, concentration and the ability to carry out day-to-day tasks. BP is a major psychiatric condition, and it is still undertreated. The causes and neural mechanisms of bipolar disorder are unclear, and diagnosis is still mostly based on psychiatric examination, furthermore the unstable character of the disorder makes diagnosis challenging. Identification of the molecular mechanisms underlying the disease may improve the diagnosis and treatment rates. Single nucleotide polymorphisms (SNP) and transcriptome profiles of patients were studied along with signalling pathways that are thought to be associated with bipolar disorder. Here, we present a computational approach that uses publicly available transcriptome data from bipolar disorder patients and healthy controls. Along with statistical analyses, data are integrated with a genome-scale metabolic model and protein-protein interaction network. Healthy individuals and bipolar disorder patients are compared based on their metabolic profiles. We hypothesize that energy metabolism alterations in bipolar disorder relate to perturbations in amino-acid metabolism and neuron-astrocyte exchange reactions. Due to changes in amino acid metabolism, neurotransmitters and their secretion from neurons and metabolic exchange pathways between neurons and astrocytes such as the glutamine-glutamate cycle are also altered. Changes in negatively charged (-1) KIV and KMV molecules are also observed, and it indicates that charge balance in the brain is highly altered in bipolar disorder. Due to this fact, we also hypothesize that positively charged lithium ions may stabilize the disturbed charge balance in neurons in addition to its effects on neurotransmission. To the best of our knowledge, our approach is unique as it is the first study using genome-scale metabolic models in neuropsychiatric research.

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.

Metabolic modeling-based drug repurposing in Glioblastoma

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

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.

Reconstruction of tissue-specific genome-scale metabolic models for human cancer stem cells

Cancer Stem Cells (CSCs) contribute to cancer aggressiveness, metastasis, chemo/radio-therapy resistance, and tumor recurrence. Recent studies emphasized the importance of metabolic reprogramming of CSCs for the maintenance and progression of the cancer phenotype through both the fulfillment of the energetic requirements and the supply of substrates fundamental for fast-cell growth, as well as through metabolite-induced epigenetic regulation. Therefore, it is of paramount importance to develop therapeutic strategies tailored to target the metabolism of CSCs. In this work, we built computational Genome-Scale Metabolic Models (GSMMs) for CSCs of different tissues. Flux simulations were then used to predict metabolic phenotypes, identify potential therapeutic targets, and spot already-known Transcription Factors (TFs), miRNAs and antimetabolites that could be used as part of drug repurposing strategies against cancer. Results were in accordance with experimental evidence, provided insights of new metabolic mechanisms for already known agents, and allowed for the identification of potential new targets and compounds that could be interesting for further in vitro and in vivo validation.

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

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.

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

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

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

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

Applications of Multi-omics Approaches for Exploring the Molecular Mechanism of Ovarian Carcinogenesis

Ovarian cancer ranks as the fifth most common cause of cancer-related death in females. The molecular mechanisms of ovarian carcinogenesis need to be explored in order to identify effective clinical therapies for ovarian cancer. Recently, multi-omics approaches have been applied to determine the mechanisms of ovarian oncogenesis at genomics (DNA), transcriptomics (RNA), proteomics (proteins), and metabolomics (metabolites) levels. Multi-omics approaches can identify some diagnostic and prognostic biomarkers and therapeutic targets for ovarian cancer, and these molecular signatures are beneficial for clarifying the development and progression of ovarian cancer. Moreover, the discovery of molecular signatures and targeted therapy strategies could noticeably improve the prognosis of ovarian cancer patients.

A Fast and Efficient Approach to Obtaining High-Purity Glioma Stem Cell Culture

Glioma is the most common and malignant primary brain tumor. Patients with malignant glioma usually have a poor prognosis due to drug resistance and disease relapse. Cancer stem cells contribute to glioma initiation, progression, resistance, and relapse. Hence, quick identification and efficient understanding of glioma stem cells (GSCs) are of profound importance for therapeutic strategies and outcomes. Ideally, therapeutic approaches will only kill cancer stem cells without harming normal neural stem cells (NSCs) that can inhibit GSCs and are often beneficial. It is key to identify the differences between cancer stem cells and normal NSCs. However, reports detailing an efficient and uniform protocol are scarce, as are comparisons between normal neural and cancer stem cells. Here, we compared different protocols and developed a fast and efficient approach to obtaining high-purity glioma stem cell by tracking observation and optimizing culture conditions. We examined the proliferative and differentiative properties confirming the identities of the GSCs with relevant markers such as Ki67, SRY-box containing gene 2, an intermediate filament protein member nestin, glial fibrillary acidic protein, and s100 calcium-binding protein (s100-beta). Finally, we identified distinct expression differences between GSCs and normal NSCs including cyclin-dependent kinase 4 and tumor protein p53. This study comprehensively describes the features of GSCs, their properties, and regulatory genes with expression differences between them and normal stem cells. Effective approaches to quickly obtaining high-quality GSCs from patients should have the potential to not only help understand the diseases and the resistances but also enable target drug screening and personalized medicine for brain tumor treatment.