Genome-scale metabolic modelling predicts biomarkers and therapeutic targets for neuropsychiatric disorders

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.

Transcriptome-based biomarker prediction for Parkinson's disease using genome-scale metabolic modeling

Parkinson's disease (PD) is the second most common neurodegenerative disease in the world. Identification of PD biomarkers is crucial for early diagnosis and to develop target-based therapeutic agents. Integrative analysis of genome-scale metabolic models (GEMs) and omics data provides a computational approach for the prediction of metabolite biomarkers. Here, we applied the TIMBR (Transcriptionally Inferred Metabolic Biomarker Response) algorithm and two modified versions of TIMBR to investigate potential metabolite biomarkers for PD. To this end, we mapped thirteen post-mortem PD transcriptome datasets from the substantia nigra region onto Human-GEM. We considered a metabolite as a candidate biomarker if its production was predicted to be more efficient by a TIMBR-family algorithm in control or PD case for the majority of the datasets. Different metrics based on well-known PD-related metabolite alterations, PD-associated pathways, and a list of 25 high-confidence PD metabolite biomarkers compiled from the literature were used to compare the prediction performance of the three algorithms tested. The modified algorithm with the highest prediction power based on the metrics was called TAMBOOR, TrAnscriptome-based Metabolite Biomarkers by On-Off Reactions, which was introduced for the first time in this study. TAMBOOR performed better in terms of capturing well-known pathway alterations and metabolite secretion changes in PD. Therefore, our tool has a strong potential to be used for the prediction of novel diagnostic biomarkers for human diseases.

Boltzmann's Theorem Revisited: Inaccurate Time-to-Action Clocks in Affective Disorders

Timely goal-oriented behavior is essential for survival and is shaped by experience. In this paper, a multileveled approach was employed, ranging from the polymorphic level through thermodynamic molecular, cellular, intracellular, extracellular, non-neuronal organelles and electrophysiological waves, attesting for signal variability. By adopting Boltzmann's theorem as a thermodynamic conceptualization of brain work, we found deviations from excitation-inhibition balance and wave decoupling, leading to wider signal variability in affective disorders compared to healthy individuals. Recent evidence shows that the overriding on-off design of clock genes paces the accuracy of the multilevel parallel sequencing clocks and that the accuracy of the time-to-action is more crucial for healthy behavioral reactions than their rapidity or delays. In affective disorders, the multilevel clocks run free and lack accuracy of responsivity to environmentally triggered time-to-action as the clock genes are not able to rescue mitochondria organelles from oxidative stress to produce environmentally-triggered energy that is required for the accurate time-to-action and maintenance of the thermodynamic equilibrium. This maintenance, in turn, is dependent on clock gene transcription of electron transporters, leading to higher signal variability and less signal accuracy in affective disorders. From a Boltzmannian thermodynamic and energy-production perspective, the option of reversibility to a healthier time-toaction, reducing entropy is implied. We employed logic gates to show deviations from healthy levelwise communication and the reversed conditions through compensations implying the role of nonneural cells and the extracellular matrix in return to excitation-inhibition balance and accuracy in the time-to-action signaling.

Multilayer Analysis of RNA Sequencing Data in Alzheimer's Disease to Unravel Molecular Mysteries

Alzheimer's disease (AD) is a complex disease, and numerous cellular events may be involved in etiology. RNAseq-based transcriptome data hold multilayer information content, which could be crucial in unraveling molecular mysteries of AD. It enables quantification of gene expression levels, identification of genomic variants, and elucidation of splicing anomalies such as exon skipping and intron retention. Additional integration of this information into protein-protein interaction networks and genome-scale metabolic models from the literature has potential to decipher functional modules and affected mechanisms for complex scenarios such as AD. In this chapter, we review the application areas of the multilayer content of RNAseq and associated integrative approaches available, with a special focus on AD.

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

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

A contribution of metabolic engineering to addressing medical problems: Metabolic flux analysis

Metabolic engineering has served as a systematic discipline for industrial biotechnology as it has offered systematic tools and methods for strain development and bioprocess optimization. Because these metabolic engineering tools and methods are concerned with the biological network of a cell with emphasis on metabolic network, they have also been applied to a range of medical problems where better understanding of metabolism has also been perceived to be important. Metabolic flux analysis (MFA) is a unique systematic approach initially developed in the metabolic engineering community, and has proved its usefulness and potential when addressing a range of medical problems. In this regard, this review discusses the contribution of MFA to addressing medical problems. For this, we i) provide overview of the milestones of MFA, ii) define two main branches of MFA, namely constraint-based reconstruction and analysis (COBRA) and isotope-based MFA (iMFA), and iii) present successful examples of their medical applications, including characterizing the metabolism of diseased cells and pathogens, and identifying effective drug targets. Finally, synergistic interactions between metabolic engineering and biomedical sciences are discussed with respect to MFA.

Neuroinflammation, Energy and Sphingolipid Metabolism Biomarkers Are Revealed by Metabolic Modeling of Autistic Brains

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders generally characterized by repetitive behaviors and difficulties in communication and social behavior. Despite its heterogeneous nature, several metabolic dysregulations are prevalent in individuals with ASD. This work aims to understand ASD brain metabolism by constructing an ASD-specific prefrontal cortex genome-scale metabolic model (GEM) using transcriptomics data to decipher novel neuroinflammatory biomarkers. The healthy and ASD-specific models are compared via uniform sampling to identify ASD-exclusive metabolic features. Noticeably, the results of our simulations and those found in the literature are comparable, supporting the accuracy of our reconstructed ASD model. We identified that several oxidative stress, mitochondrial dysfunction, and inflammatory markers are elevated in ASD. While oxidative phosphorylation fluxes were similar for healthy and ASD-specific models, and the fluxes through the pathway were nearly undisturbed, the tricarboxylic acid (TCA) fluxes indicated disruptions in the pathway. Similarly, the secretions of mitochondrial dysfunction markers such as pyruvate are found to be higher, as well as the activities of oxidative stress marker enzymes like alanine and aspartate aminotransferases (ALT and AST) and glutathione-disulfide reductase (GSR). We also detected abnormalities in the sphingolipid metabolism, which has been implicated in many inflammatory and immune processes, but its relationship with ASD has not been thoroughly explored in the existing literature. We suggest that important sphingolipid metabolites, such as sphingosine-1-phosphate (S1P), ceramide, and glucosylceramide, may be promising biomarkers for the diagnosis of ASD and provide an opportunity for the adoption of early intervention for young children.

Context-Specific Genome-Scale Metabolic Modelling and Its Application to the Analysis of COVID-19 Metabolic Signatures

Genome-scale metabolic models (GEMs) have found numerous applications in different domains, ranging from biotechnology to systems medicine. Herein, we overview the most popular algorithms for the automated reconstruction of context-specific GEMs using high-throughput experimental data. Moreover, we describe different datasets applied in the process, and protocols that can be used to further automate the model reconstruction and validation. Finally, we describe recent COVID-19 applications of context-specific GEMs, focusing on the analysis of metabolic implications, identification of biomarkers and potential drug targets.

Inferring miRNA sponge modules across major neuropsychiatric disorders

The role of non-coding RNAs in neuropsychiatric disorders (NPDs) is an emerging field of study. The long non-coding RNAs (lncRNAs) are shown to sponge the microRNAs (miRNAs) from interacting with their target mRNAs. Investigating the sponge activity of lncRNAs in NPDs will provide further insights into biological mechanisms and help identify disease biomarkers. In this study, a large-scale inference of the lncRNA-related miRNA sponge network of pan-neuropsychiatric disorders, including autism spectrum disorder (ASD), schizophrenia (SCZ), and bipolar disorder (BD), was carried out using brain transcriptomic (RNA-Seq) data. The candidate miRNA sponge modules were identified based on the co-expression pattern of non-coding RNAs, sharing of miRNA binding sites, and sensitivity canonical correlation. miRNA sponge modules are associated with chemical synaptic transmission, nervous system development, metabolism, immune system response, ribosomes, and pathways in cancer. The identified modules showed similar and distinct gene expression patterns depending on the neuropsychiatric condition. The preservation of miRNA sponge modules was shown in the independent brain and blood-transcriptomic datasets of NPDs. We also identified miRNA sponging lncRNAs that may be potential diagnostic biomarkers for NPDs. Our study provides a comprehensive resource on miRNA sponging in NPDs.

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.

Gut Microbial Dysbiosis and Cognitive Impairment in Bipolar Disorder: Current Evidence

Recent studies have reported that the gut microbiota influences mood and cognitive function through the gut-brain axis, which is involved in the pathophysiology of neurocognitive and mental disorders, including Parkinson’s disease, Alzheimer’s disease, and schizophrenia. These disorders have similar pathophysiology to that of cognitive dysfunction in bipolar disorder (BD), including neuroinflammation and dysregulation of various neurotransmitters (i.e., serotonin and dopamine). There is also emerging evidence of alterations in the gut microbial composition of patients with BD, suggesting that gut microbial dysbiosis contributes to disease progression and cognitive impairment in BD. Therefore, microbiota-centered treatment might be an effective adjuvant therapy for BD-related cognitive impairment. Given that studies focusing on connections between the gut microbiota and BD-related cognitive impairment are lagging behind those on other neurocognitive disorders, this review sought to explore the potential mechanisms of how gut microbial dysbiosis affects cognitive function in BD and identify potential microbiota-centered treatment.

Elucidating the metabolic characteristics of pancreatic β-cells from patients with type 2 diabetes (T2D) using a genome-scale metabolic modeling

Diabetes is a global health problem caused primarily by the inability of pancreatic β-cells to secrete adequate insulin. Despite extensive research, the identity of factors contributing to the dysregulated metabolism-secretion coupling in the β-cells remains elusive. The present study attempts to capture some of these factors responsible for the impaired β-cell metabolism-secretion coupling that contributes to diabetes pathogenesis. The metabolic-flux profiles of pancreatic β-cells were predicted using genome-scale metabolic modeling for ten diabetic patients and ten control subjects. Analysis of these flux states shows reduction in the mitochondrial fatty acid oxidation and mitochondrial oxidative phosphorylation pathways, that leads to decreased insulin secretion in diabetes. We also observed elevated reactive oxygen species (ROS) generation through peroxisomal fatty acid β-oxidation. In addition, cellular antioxidant defense systems were found to be attenuated in diabetes. Our analysis also uncovered the possible changes in the plasma metabolites in diabetes due to the β-cells failure. These efforts subsequently led to the identification of seven metabolites associated with cardiovascular disease (CVD) pathogenesis, thus establishing its link as a secondary complication of diabetes.

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.

Genome-scale metabolic modelling of SARS-CoV-2 in cancer cells reveals an increased shift to glycolytic energy production

Cancer is considered a high‐risk condition for severe illness resulting from COVID‐19. The interaction between severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) and human metabolism is key to elucidating the risk posed by COVID‐19 for cancer patients and identifying effective treatments, yet it is largely uncharacterised on a mechanistic level. We present a genome‐scale map of short‐term metabolic alterations triggered by SARS‐CoV‐2 infection of cancer cells. Through transcriptomic‐ and proteomic‐informed genome‐scale metabolic modelling, we characterise the role of RNA and fatty acid biosynthesis in conjunction with a rewiring in energy production pathways and enhanced cytokine secretion. These findings link together complementary aspects of viral invasion of cancer cells, while providing mechanistic insights that can inform the development of treatment strategies.

Transcriptomics Analysis Reveals Shared Pathways in Peripheral Blood Mononuclear Cells and Brain Tissues of Patients With Schizophrenia

Background: Schizophrenia is a serious mental disorder with complicated biological mechanisms. Few studies explore the transcriptional features that are shared in brain tissue and peripheral blood. In the present study, we aimed to explore the biological pathways with similar expression patterns in both peripheral blood mononuclear cells (PBMCs) and brain tissues. Methods: The present study used transcriptomics technology to detect mRNA expression of PBMCs of 10 drug-naïve patients with schizophrenia and 20 healthy controls. Transcriptome data sets of brain tissue of patients with schizophrenia downloaded from public databases were also analyzed in our study. The biological pathways with similar expression patterns in the PBMCs and brain tissues were uncovered by differential expression analysis, weighted gene co-expression network analysis (WGCNA), and pathway analysis. Finally, the expression levels of differential expressed genes (DEGs) were validated by real-time fluorescence quantitative polymerase chain reaction (qPCR) in another 12 drug-naïve patients with schizophrenia and 12 healthy controls. Results: We identified 542 DEGs, 51 DEGs, 732 DEGs, and 104 DEGs in PBMCs, dorsolateral prefrontal cortex, anterior cingulate gyrus, and nucleus accumbent, respectively. Five DEG clusters were recognized as having similar gene expression patterns in PBMCs and brain tissues by WGCNA. The pathway analysis illustrates that these DEG clusters are mainly enriched in several biological pathways that are related to phospholipid metabolism, ribosome signal transduction, and mitochondrial oxidative phosphorylation. The differential significance of PLAAT3, PLAAT4, PLD2, RPS29, RPL30, COX7C, COX7A2, NDUFAF2, and ATP5ME were confirmed by qPCR. Conclusions: This study finds that the pathways associated with phospholipid metabolism, ribosome signal transduction, and energy metabolism have similar expression patterns in PBMCs and brain tissues of patients with schizophrenia. Our results supply a novel insight for revealing the pathogenesis of schizophrenia and might offer a new approach to explore potential biological markers of peripheral blood in schizophrenia.