Integration of transcriptomic data in a genome-scale metabolic model to investigate the link between obesity and breast cancer

ICON-GEMs: integration of co-expression network in genome-scale metabolic models, shedding light through systems biology

BACKGROUND

Flux Balance Analysis (FBA) is a key metabolic modeling method used to simulate cellular metabolism under steady-state conditions. Its simplicity and versatility have led to various strategies incorporating transcriptomic and proteomic data into FBA, successfully predicting flux distribution and phenotypic results. However, despite these advances, the untapped potential lies in leveraging gene-related connections like co-expression patterns for valuable insights.

RESULTS

To fill this gap, we introduce ICON-GEMs, an innovative constraint-based model to incorporate gene co-expression network into the FBA model, facilitating more precise determination of flux distributions and functional pathways. In this study, transcriptomic data from both Escherichia coli and Saccharomyces cerevisiae were integrated into their respective genome-scale metabolic models. A comprehensive gene co-expression network was constructed as a global view of metabolic mechanism of the cell. By leveraging quadratic programming, we maximized the alignment between pairs of reaction fluxes and the correlation of their corresponding genes in the co-expression network. The outcomes notably demonstrated that ICON-GEMs outperformed existing methodologies in predictive accuracy. Flux variabilities over subsystems and functional modules also demonstrate promising results. Furthermore, a comparison involving different types of biological networks, including protein-protein interactions and random networks, reveals insights into the utilization of the co-expression network in genome-scale metabolic engineering.

CONCLUSION

ICON-GEMs introduce an innovative constrained model capable of simultaneous integration of gene co-expression networks, ready for board application across diverse transcriptomic data sets and multiple organisms. It is freely available as open-source at https://github.com/ThummaratPaklao/ICOM-GEMs.git .

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.

TumorMet: A repository of tumor metabolic networks derived from context-specific Genome-Scale Metabolic Models

Studies about the metabolic alterations during tumorigenesis have increased our knowledge of the underlying mechanisms and consequences, which are important for diagnostic and therapeutic investigations. In this scenario and in the era of systems biology, metabolic networks have become a powerful tool to unravel the complexity of the cancer metabolic machinery and the heterogeneity of this disease. Here, we present TumorMet, a repository of tumor metabolic networks extracted from context-specific Genome-Scale Metabolic Models, as a benchmark for graph machine learning algorithms and network analyses. This repository has an extended scope for use in graph classification, clustering, community detection, and graph embedding studies. Along with the data, we developed and provided Met2Graph, an R package for creating three different types of metabolic graphs, depending on the desired nodes and edges: Metabolites-, Enzymes-, and Reactions-based graphs. This package allows the easy generation of datasets for downstream analysis.

Measurement(s)	gene expression, metabolic relationships
Technology Type(s)	Genome Scale Metabolic Models; Computational network biology
Sample Characteristic - Organism	Homo sapiens

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.

All Driven by Energy Demand? Integrative Comparison of Metabolism of Enterococcus faecalis Wildtype and a Glutamine Synthase Mutant

Lactic acid bacteria (LAB) play a significant role in biotechnology, e.g., food industry and also in human health. Many LAB genera have developed a multidrug resistance in the past few years, causing a serious problem in controlling hospital germs worldwide. Enterococcus faecalis accounts for a large part of the human infections caused by LABs. Therefore, studying its adaptive metabolism under various environmental conditions is particularly important to promote the development of new therapeutic approaches. In this study, we investigated the effect of glutamine auxotrophy (ΔglnA mutant) on metabolic and proteomic adaptations of E. faecalis in response to a changing pH in its environment. Changing pH values are part of the organism's natural environment in the human body and play a role in the food industry. We compared the results with those of the wildtype. Using a genome-scale metabolic model constrained by metabolic and proteomic data, our integrative method allows us to understand the bigger picture of the adaptation strategies of this bacterium. The study showed that energy demand is the decisive factor in adapting to a new environmental pH. The energy demand of the mutant was higher at all conditions. It has been reported that ΔglnA mutants of bacteria are energetically less effective. With the aid of our data and model we are able to explain this phenomenon as a consequence of a failure to regulate glutamine uptake and the costs for the import of glutamine and the export of ammonium. Methodologically, it became apparent that taking into account the nonspecificity of amino acid transporters is important for reproducing metabolic changes with genome-scale models because it affects energy balance. IMPORTANCE The integration of new pH-dependent experimental data on metabolic uptake and release fluxes, as well as of proteome data with a genome-scale computational model of a glutamine synthetase mutant of E. faecalis is used and compared with those of the wildtype to understand why glutamine auxotrophy results in a less efficient metabolism and how-in comparison with the wildtype-the glutamine synthetase knockout impacts metabolic adjustments during acidification or simply exposure to lower pH. We show that forced glutamine auxotrophy causes more energy demand and that this is likely due to a disregulated glutamine uptake. Proteome changes during acidification observed for the mutant resemble those of the wildtype with the exception of glycolysis-related genes, as the mutant is already energetically stressed at a higher pH and the respective proteome changes were in effect.

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.

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

Endothelial cell plasticity at the single-cell level

The vascular endothelium is characterized by a remarkable level of plasticity, which is the driving force not only of physiological repair/remodeling of adult tissues but also of pathological angiogenesis. The resulting heterogeneity of endothelial cells (ECs) makes targeting the endothelium challenging, no less because many EC phenotypes are yet to be identified and functionally inventorized. Efforts to map the vasculature at the single-cell level have been instrumental to capture the diversity of EC types and states at a remarkable depth in both normal and pathological states. Here, we discuss new EC subtypes and functions emerging from recent single-cell studies in health and disease. Interestingly, such studies revealed distinct metabolic gene signatures in different EC phenotypes, which deserve further consideration for therapy. We highlight how this metabolic targeting strategy could potentially be used to promote (for tissue repair) or block (in tumor) angiogenesis in a tissue or even vascular bed-specific manner.

Mapping the Metabolic Networks of Tumor Cells and Cancer-Associated Fibroblasts

Metabolism is considered to be the core of all cellular activity. Thus, extensive studies of metabolic processes are ongoing in various fields of biology, including cancer research. Cancer cells are known to adapt their metabolism to sustain high proliferation rates and survive in unfavorable environments with low oxygen and nutrient concentrations. Hence, targeting cancer cell metabolism is a promising therapeutic strategy in cancer research. However, cancers consist not only of genetically altered tumor cells but are interwoven with endothelial cells, immune cells and fibroblasts, which together with the extracellular matrix (ECM) constitute the tumor microenvironment (TME). Cancer-associated fibroblasts (CAFs), which are linked to poor prognosis in different cancer types, are one important component of the TME. CAFs play a significant role in reprogramming the metabolic landscape of tumor cells, but how, and in what manner, this interaction takes place remains rather unclear. This review aims to highlight the metabolic landscape of tumor cells and CAFs, including their recently identified subtypes, in different tumor types. In addition, we discuss various in vitro and in vivo metabolic techniques as well as different in silico computational tools that can be used to identify and characterize CAF-tumor cell interactions. Finally, we provide our view on how mapping the complex metabolic networks of stromal-tumor metabolism will help in finding novel metabolic targets for cancer treatment.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Leveraging genome-scale metabolic models for human health applications

Genome-scale metabolic modeling is a scalable and extensible computational method for analyzing and predicting biological function. With the ongoing improvements in computational methods and experimental capabilities, genome-scale metabolic models (GEMs) are demonstrating utility in addressing human health applications. The initial areas of highest impact are likely to be health applications where disease states involve metabolic changes. In this review, we focus on recent application of GEMs to studying cancer and the human microbiome by describing the enabling methodologies and outcomes of these studies. We conclude with proposing some areas of research that are likely to arise as a result of recent methodological advances.

Computational approaches to understanding nutrient metabolism and metabolic disorders

Computational methods are becoming more and more essential to elucidate biological systems. Many different approaches exist with pros and cons. This paper reviews the most useful technologies focusing on nutrient metabolism and metabolic disorders. Space limitation prevents from exploring the examples in details, but pointers to the relevant papers are reported.

A systematic strategy using a reconstructed genome-scale metabolic network for pathogen Streptococcuspneumoniae D39 to find novel potential drug targets

Streptococcus pneumoniae is a Gram-positive bacterium that is one of the major causes of various infections such as pneumonia, meningitis, otitis media and endocarditis. Since antibiotic resistance of S. pneumoniae is pointed out as a challenge in the treatment of these infections, more studies are required to focus on disease prevention. In this research, a first manually curated genome-scale metabolic network of the pathogen S. pneumoniae D39 was reconstructed based on its genome annotation data, and biochemical knowledge from literature and databases. The model was validated by amino acid auxotrophies, gene essentiality analysis, and different carbohydrate sources. Then, a two-stage strategy was developed to find target genes for growth reduction of the pathogen and their importance in the various infection sites. In the first stage, growth-associated genes were identified by integration of transcriptomic data with the model and in the second stage, the importance of each gene in the metabolism for growth was evaluated using principal component analysis. The reports presented in the literature confirm the effect of some found genes on the growth of S. pneumoniae.

Clustering analysis of tumor metabolic networks

BACKGROUND

Biological networks are representative of the diverse molecular interactions that occur within cells. Some of the commonly studied biological networks are modeled through protein-protein interactions, gene regulatory, and metabolic pathways. Among these, metabolic networks are probably the most studied, as they directly influence all physiological processes. Exploration of biochemical pathways using multigraph representation is important in understanding complex regulatory mechanisms. Feature extraction and clustering of these networks enable grouping of samples obtained from different biological specimens. Clustering techniques separate networks depending on their mutual similarity.

RESULTS

We present a clustering analysis on tissue-specific metabolic networks for single samples from three primary tumor sites: breast, lung, and kidney cancer. The metabolic networks were obtained by integrating genome scale metabolic models with gene expression data. We performed network simplification to reduce the computational time needed for the computation of network distances. We empirically proved that networks clustering can characterize groups of patients in multiple conditions.

CONCLUSIONS

We provide a computational methodology to explore and characterize the metabolic landscape of tumors, thus providing a general methodology to integrate analytic metabolic models with gene expression data. This method represents a first attempt in clustering large scale metabolic networks. Moreover, this approach gives the possibility to get valuable information on what are the effects of different conditions on the overall metabolism.

Bioengineering tools for the production of pharmaceuticals: current perspective and future outlook

The bioengineering tools have significant advantages through less time-consuming and utilized as a promising stage for the production of pharmaceutical bioproducts under the single platform. This review highlighted the advantages and current improvement in the plant, animal and microbial bioengineering tools and outlines feasible approaches by biological and process's bioengineering levels for advancing the economic feasibility of pharmaceutical's production. The critical analysis results revealed that system biology and synthetic biology along with advanced bioengineering tools like transcriptome, proteome, metabolome and nano bioengineering tools have shown a promising impact on the development of pharmaceutical's bioproducts. Tools to overcome and resolve the accompanying encounters of pharmaceutical's production that include nano bioengineering tools are also discussed. As a summary and prospect, it also gives new insight into the challenges and possible breakthrough of the development of pharmaceutical's bioproducts through bioengineering tools.

The 2017 Network Tools and Applications in Biology (NETTAB) workshop: aims, topics and outcomes

The 17th International NETTAB workshop was held in Palermo, Italy, on October 16-18, 2017. The special topic for the meeting was "Methods, tools and platforms for Personalised Medicine in the Big Data Era", but the traditional topics of the meeting series were also included in the event. About 40 scientific contributions were presented, including four keynote lectures, five guest lectures, and many oral communications and posters. Also, three tutorials were organised before and after the workshop. Full papers from some of the best works presented in Palermo were submitted for this Supplement of BMC Bioinformatics. Here, we provide an overview of meeting aims and scope. We also shortly introduce selected papers that have been accepted for publication in this Supplement, for a complete presentation of the outcomes of the meeting.