Integration of clinical data with a genome-scale metabolic model of the human adipocyte

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

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.

Reconstruction of cell-specific models capturing the influence of metabolism on DNA methylation in cancer

The imbalance of epigenetic regulatory mechanisms such as DNA methylation, which can promote aberrant gene expression profiles without affecting the DNA sequence, may cause the deregulation of signaling, regulatory, and metabolic processes, contributing to a cancerous phenotype. Since some metabolites are substrates and cofactors of epigenetic regulators, their availability can be affected by characteristic cancer cell metabolic shifts, feeding cancer onset and progression through epigenetic deregulation. Hence, there is a need to study the influence of cancer metabolic reprogramming in DNA methylation to design new effective treatments. In this study, a generic Genome-Scale Metabolic Model (GSMM) of a human cell, integrating DNA methylation or demethylation reactions, was obtained and used for the reconstruction of Genome-Scale Metabolic Models enhanced with Enzymatic Constraints using Kinetic and Omics data (GECKOs) of 31 cancer cell lines. Furthermore, cell-line-specific DNA methylation levels were included in the models, as coefficients of a DNA composition pseudo-reaction, to depict the influence of metabolism over global DNA methylation in each of the cancer cell lines. Flux simulations demonstrated the ability of these models to provide simulated fluxes of exchange reactions similar to the equivalent experimentally measured uptake/secretion rates and to make good functional predictions. In addition, simulations found metabolic pathways, reactions and enzymes directly or inversely associated with the gene promoter methylation. Two potential candidates for targeted cancer epigenetic therapy were identified.

Thyroid cancer cell metabolism: A glance into cell culture system-based metabolomics approaches

Thyroid cancer is the most common malignancy of the endocrine system and the seventh most prevalent cancer in women worldwide. It is a complex and diverse disease characterized by heterogeneity, underscoring the importance of understanding the underlying metabolic alterations within tumor cells. Metabolomics technologies offer a powerful toolset to explore and identify endogenous and exogenous biochemical reaction products, providing crucial insights into the intricate metabolic pathways and processes within living cells. Metabolism plays a central role in cell function, making metabolomics a valuable reflection of a cell's phenotype. In the OMICs era, metabolomics analysis of cells brings numerous advantages over existing methods, propelling cell metabolomics as an emerging field with vast potential for investigating metabolic pathways and their perturbation in pathophysiological conditions. This review article aims to look into recent developments in applying metabolomics for characterizing and interpreting the cellular metabolome in thyroid cancer cell lines, exploring their unique metabolic characteristics. Understanding the metabolic alterations in tumor cells can lead to the identification of critical nodes in the metabolic network that could be targeted for therapeutic intervention.

Integrating glycolysis, citric acid cycle, pentose phosphate pathway, and fatty acid beta-oxidation into a single computational model

The metabolic network of a living cell is highly intricate and involves complex interactions between various pathways. In this study, we propose a computational model that integrates glycolysis, the pentose phosphate pathway (PPP), the fatty acids beta-oxidation, and the tricarboxylic acid cycle (TCA cycle) using queueing theory. The model utilizes literature data on metabolite concentrations and enzyme kinetic constants to calculate the probabilities of individual reactions occurring on a microscopic scale, which can be viewed as the reaction rates on a macroscopic scale. However, it should be noted that the model has some limitations, including not accounting for all the reactions in which the metabolites are involved. Therefore, a genetic algorithm (GA) was used to estimate the impact of these external processes. Despite these limitations, our model achieved high accuracy and stability, providing real-time observation of changes in metabolite concentrations. This type of model can help in better understanding the mechanisms of biochemical reactions in cells, which can ultimately contribute to the prevention and treatment of aging, cancer, metabolic diseases, and neurodegenerative disorders.

Integrating Omics Data in Genome-Scale Metabolic Modeling: A Methodological Perspective for Precision Medicine

Recent advancements in omics technologies have generated a wealth of biological data. Integrating these data within mathematical models is essential to fully leverage their potential. Genome-scale metabolic models (GEMs) provide a robust framework for studying complex biological systems. GEMs have significantly contributed to our understanding of human metabolism, including the intrinsic relationship between the gut microbiome and the host metabolism. In this review, we highlight the contributions of GEMs and discuss the critical challenges that must be overcome to ensure their reproducibility and enhance their prediction accuracy, particularly in the context of precision medicine. We also explore the role of machine learning in addressing these challenges within GEMs. The integration of omics data with GEMs has the potential to lead to new insights, and to advance our understanding of molecular mechanisms in human health and disease.

A holistic view of muscle metabolic reprogramming through personalized metabolic modeling in newly diagnosed diabetic patients

Type 2 diabetes mellitus (T2DM) is a challenging and progressive metabolic disease caused by insulin resistance. Skeletal muscle is the major insulin-sensitive tissue that plays a pivotal role in blood sugar homeostasis. Dysfunction of muscle metabolism is implicated in the disturbance of glucose homeostasis, the development of insulin resistance, and T2DM. Understanding metabolism reprogramming in newly diagnosed patients provides opportunities for early diagnosis and treatment of T2DM as a challenging disease to manage. Here, we applied a system biology approach to investigate metabolic dysregulations associated with the early stage of T2DM. We first reconstructed a human muscle-specific metabolic model. The model was applied for personalized metabolic modeling and analyses in newly diagnosed patients. We found that several pathways and metabolites, mainly implicating in amino acids and lipids metabolisms, were dysregulated. Our results indicated the significance of perturbation of pathways implicated in building membrane and extracellular matrix (ECM). Dysfunctional metabolism in these pathways possibly interrupts the signaling process and develops insulin resistance. We also applied a machine learning method to predict potential metabolite markers of insulin resistance in skeletal muscle. 13 exchange metabolites were predicted as the potential markers. The efficiency of these markers in discriminating insulin-resistant muscle was successfully validated.

A comprehensive mechanistic model of adipocyte signaling with layers of confidence

Adipocyte signaling, normally and in type 2 diabetes, is far from fully understood. We have earlier developed detailed dynamic mathematical models for several well-studied, partially overlapping, signaling pathways in adipocytes. Still, these models only cover a fraction of the total cellular response. For a broader coverage of the response, large-scale phosphoproteomic data and systems level knowledge on protein interactions are key. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. We have developed a method to first establish a core model by connecting existing models of adipocyte cellular signaling for: (1) lipolysis and fatty acid release, (2) glucose uptake, and (3) the release of adiponectin. Next, we use publicly available phosphoproteome data for the insulin response in adipocytes together with prior knowledge on protein interactions, to identify phosphosites downstream of the core model. In a parallel pairwise approach with low computation time, we test whether identified phosphosites can be added to the model. We iteratively collect accepted additions into layers and continue the search for phosphosites downstream of these added layers. For the first 30 layers with the highest confidence (311 added phosphosites), the model predicts independent data well (70-90% correct), and the predictive capability gradually decreases when we add layers of decreasing confidence. In total, 57 layers (3059 phosphosites) can be added to the model with predictive ability kept. Finally, our large-scale, layered model enables dynamic simulations of systems-wide alterations in adipocytes in type 2 diabetes.

Sexual dimorphism in mitochondrial dysfunction and diabetes mellitus: evidence from a population-based cohort study

BACKGROUND

Pathophysiological mechanisms underlying sex-based differences in diabetes remain poorly understood. Mitochondrial metabolite methylmalonic acid (MMA) accumulation reflects mitochondrial dysfunction which is involved in sex-specific pathophysiological responses biologically. We aimed to investigate the sex-specific associations between mortality risk and MMA in adults with the presence or absence of type 2 diabetes.

METHODS

This cohort study included 24,164 adults (12,123 females and 12,041 males) from the NHANES study during 1999-2014. Both sexes were separately categorized as those with no diabetes, prediabetes, undiagnosed diabetes, and diagnosed diabetes. Circulating MMA level was measured at baseline by mass-spectrometric detection. Mortality status was ascertained from baseline until December 31, 2015.

RESULTS

During a median follow-up of 11.1 years, 3375 deaths were documented. Males had a particularly higher mortality than females in adults with diagnosed diabetes compared to differences in those with no diabetes, prediabetes and undiagnosed diabetes (sex differences in mortality rate per 1000 person-years across diabetic status: 0.62, 1.44, 5.78, and 9.77, p < 0.001). Notably, the sex-specific difference in associations between MMA and mortality was significant only in adults with diagnosed diabetes (p for interaction = 0.028), not in adults with no diabetes and prediabetes. Adjusted HRs (95%CIs) per doubling of MMA for all-cause mortality were 1.19 (1.04-1.37) in females with diagnosed diabetes versus 1.58 (1.36-1.86) in male counterparts. In addition, MMA levels had an insignificant or weak correlation with sex hormone profiles at baseline, regardless of diabetes status and sex.

CONCLUSIONS

Sex difference in mortality risk was especially significant in diagnosed type 2 diabetes. Increasing equivalent exposure to mitochondrial metabolite MMA was associated with a greater excess risk of future mortality in males with diabetes than in females.

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.

Systematic diet composition swap in a mouse genome-scale metabolic model reveals determinants of obesogenic diet metabolism in liver cancer

Dietary nutrient availability and gene expression, together, influence tissue metabolic activity. Here, we explore whether altering dietary nutrient composition in the context of mouse liver cancer suffices to overcome chronic gene expression changes that arise from tumorigenesis and western-style diet (WD). We construct a mouse genome-scale metabolic model and estimate metabolic fluxes in liver tumors and non-tumoral tissue after computationally varying the composition of input diet. This approach, called Systematic Diet Composition Swap (SyDiCoS), revealed that, compared to a control diet, WD increases production of glycerol and succinate irrespective of specific tissue gene expression patterns. Conversely, differences in fatty acid utilization pathways between tumor and non-tumor liver are amplified with WD by both dietary carbohydrates and lipids together. Our data suggest that combined dietary component modifications may be required to normalize the distinctive metabolic patterns that underlie selective targeting of tumor metabolism.

Combined metabolic activators improve metabolic functions in the animal models of neurodegenerative diseases

BACKGROUND

Neurodegenerative diseases (NDDs), including Alzheimer's disease (AD) and Parkinson's disease (PD), are associated with metabolic abnormalities. Integrative analysis of human clinical data and animal studies have contributed to a better understanding of the molecular and cellular pathways involved in the progression of NDDs. Previously, we have reported that the combined metabolic activators (CMA), which include the precursors of nicotinamide adenine dinucleotide and glutathione can be utilized to alleviate metabolic disorders by activating mitochondrial metabolism.

METHODS

We first analysed the brain transcriptomics data from AD patients and controls using a brain-specific genome-scale metabolic model (GEM). Then, we investigated the effect of CMA administration in animal models of AD and PD. We evaluated pathological and immunohistochemical findings of brain and liver tissues. Moreover, PD rats were tested for locomotor activity and apomorphine-induced rotation.

FINDINGS

Analysis of transcriptomics data with GEM revealed that mitochondrial dysfunction is involved in the underlying molecular pathways of AD. In animal models of AD and PD, we showed significant damage in the high-fat diet groups' brain and liver tissues compared to the chow diet. The histological analyses revealed that hyperemia, degeneration and necrosis in neurons were improved by CMA administration in both AD and PD animal models. These findings were supported by immunohistochemical evidence of decreased immunoreactivity in neurons. In parallel to the improvement in the brain, we also observed dramatic metabolic improvement in the liver tissue. CMA administration also showed a beneficial effect on behavioural functions in PD rats.

INTERPRETATION

Overall, we showed that CMA administration significantly improved behavioural scores in parallel with the neurohistological outcomes in the AD and PD animal models and is a promising treatment for improving the metabolic parameters and brain functions in NDDs.

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.

Quantitative modeling of human liver reveals dysregulation of glycosphingolipid pathways in nonalcoholic fatty liver disease

Nonalcoholic fatty liver disease (NAFLD) is an increasingly prevalent disease that is associated with multiple metabolic disturbances, yet the metabolic pathways underlying its progression are poorly understood. Here, we studied metabolic pathways of the human liver across the full histological spectrum of NAFLD. We analyzed whole liver tissue transcriptomics and serum metabolomics data obtained from a large, prospectively enrolled cohort of 206 histologically characterized patients derived from the European NAFLD Registry and developed genome-scale metabolic models (GEMs) of human hepatocytes at different stages of NAFLD. We identified several metabolic signatures in the liver and blood of these patients, specifically highlighting the alteration of vitamins (A, E) and glycosphingolipids, and their link with complex glycosaminoglycans in advanced fibrosis. Furthermore, we derived GEMs and identified metabolic signatures of three common NAFLD-associated gene variants (PNPLA3, TM6SF2, and HSD17B13). The study demonstrates dysregulated liver metabolic pathways which may contribute to the progression of NAFLD.

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.

Robust identification of temporal biomarkers in longitudinal omics studies

MOTIVATION

Longitudinal studies increasingly collect rich 'omics' data sampled frequently over time and across large cohorts to capture dynamic health fluctuations and disease transitions. However, the generation of longitudinal omics data has preceded the development of analysis tools that can efficiently extract insights from such data. In particular, there is a need for statistical frameworks that can identify not only which omics features are differentially regulated between groups but also over what time intervals. Additionally, longitudinal omics data may have inconsistencies, including non-uniform sampling intervals, missing data points, subject dropout and differing numbers of samples per subject.

RESULTS

In this work, we developed OmicsLonDA, a statistical method that provides robust identification of time intervals of temporal omics biomarkers. OmicsLonDA is based on a semi-parametric approach, in which we use smoothing splines to model longitudinal data and infer significant time intervals of omics features based on an empirical distribution constructed through a permutation procedure. We benchmarked OmicsLonDA on five simulated datasets with diverse temporal patterns, and the method showed specificity greater than 0.99 and sensitivity greater than 0.87. Applying OmicsLonDA to the iPOP cohort revealed temporal patterns of genes, proteins, metabolites and microbes that are differentially regulated in male versus female subjects following a respiratory infection. In addition, we applied OmicsLonDA to a longitudinal multi-omics dataset of pregnant women with and without preeclampsia, and OmicsLonDA identified potential lipid markers that are temporally significantly different between the two groups.

AVAILABILITY AND IMPLEMENTATION

We provide an open-source R package (https://bioconductor.org/packages/OmicsLonDA), to enable widespread use.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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

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

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

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

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

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

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

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

Innovations in Genomics and Big Data Analytics for Personalized Medicine and Health Care: A Review

Big data in health care is a fast-growing field and a new paradigm that is transforming case-based studies to large-scale, data-driven research. As big data is dependent on the advancement of new data standards, technology, and relevant research, the future development of big data applications holds foreseeable promise in the modern day health care revolution. Enormously large, rapidly growing collections of biomedical omics-data (genomics, proteomics, transcriptomics, metabolomics, glycomics, etc.) and clinical data create major challenges and opportunities for their analysis and interpretation and open new computational gateways to address these issues. The design of new robust algorithms that are most suitable to properly analyze this big data by taking into account individual variability in genes has enabled the creation of precision (personalized) medicine. We reviewed and highlighted the significance of big data analytics for personalized medicine and health care by focusing mostly on machine learning perspectives on personalized medicine, genomic data models with respect to personalized medicine, the application of data mining algorithms for personalized medicine as well as the challenges we are facing right now in big data analytics.

Reconstruction of a generic genome-scale metabolic network for chicken: Investigating network connectivity and finding potential biomarkers

Chicken is the first sequenced avian that has a crucial role in human life for its meat and egg production. Because of various metabolic disorders, study the metabolism of chicken cell is important. Herein, the first genome-scale metabolic model of a chicken cell named iES1300, consists of 2427 reactions, 2569 metabolites, and 1300 genes, was reconstructed manually based on KEGG, BiGG, CHEBI, UNIPROT, REACTOME, and MetaNetX databases. Interactions of metabolic genes for growth were examined for E. coli, S. cerevisiae, human, and chicken metabolic models. The results indicated robustness to genetic manipulation for iES1300 similar to the results for human. iES1300 was integrated with transcriptomics data using algorithms and Principal Component Analysis was applied to compare context-specific models of the normal, tumor, lean and fat cell lines. It was found that the normal model has notable metabolic flexibility in the utilization of various metabolic pathways, especially in metabolic pathways of the carbohydrate metabolism, compared to the others. It was also concluded that the fat and tumor models have similar growth metabolisms and the lean chicken model has a more active lipid and carbohydrate metabolism.

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

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

Systems-based approaches to study immunometabolism

Technical advances at the interface of biology and computation, such as single-cell RNA-sequencing (scRNA-seq), reveal new layers of complexity in cellular systems. An emerging area of investigation using the systems biology approach is the study of the metabolism of immune cells. The diverse spectra of immune cell phenotypes, sparsity of immune cell numbers in vivo, limitations in the number of metabolites identified, dynamic nature of cellular metabolism and metabolic fluxes, tissue specificity, and high dependence on the local milieu make investigations in immunometabolism challenging, especially at the single-cell level. In this review, we define the systemic nature of immunometabolism, summarize cell- and system-based approaches, and introduce mathematical modeling approaches for systems interrogation of metabolic changes in immune cells. We close the review by discussing the applications and shortcomings of metabolic modeling techniques. With systems-oriented studies of metabolism expected to become a mainstay of immunological research, an understanding of current approaches toward systems immunometabolism will help investigators make the best use of current resources and push the boundaries of the discipline.

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

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

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.

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.

Comparison of metabolic states using genome-scale metabolic models

Genome-scale metabolic models (GEMs) are comprehensive knowledge bases of cellular metabolism and serve as mathematical tools for studying biological phenotypes and metabolic states or conditions in various organisms and cell types. Given the sheer size and complexity of human metabolism, selecting parameters for existing analysis methods such as metabolic objective functions and model constraints is not straightforward in human GEMs. In particular, comparing several conditions in large GEMs to identify condition- or disease-specific metabolic features is challenging. In this study, we showcase a scalable, model-driven approach for an in-depth investigation and comparison of metabolic states in large GEMs which enables identifying the underlying functional differences. Using a combination of flux space sampling and network analysis, our approach enables extraction and visualisation of metabolically distinct network modules. Importantly, it does not rely on known or assumed objective functions. We apply this novel approach to extract the biochemical differences in adipocytes arising due to unlimited vs blocked uptake of branched-chain amino acids (BCAAs, considered as biomarkers in obesity) using a human adipocyte GEM (iAdipocytes1809). The biological significance of our approach is corroborated by literature reports confirming our identified metabolic processes (TCA cycle and Fatty acid metabolism) to be functionally related to BCAA metabolism. Additionally, our analysis predicts a specific altered uptake and secretion profile indicating a compensation for the unavailability of BCAAs. Taken together, our approach facilitates determining functional differences between any metabolic conditions of interest by offering a versatile platform for analysing and comparing flux spaces of large metabolic networks.

Cellular metabolism is a highly complex and interconnected system. As many lifestyle diseases in humans have a strong metabolic component, it is important to understand metabolic differences between healthy and diseased states. In systems biology, metabolic behaviours are investigated using genome-scale metabolic models. In addition to the sheer size and complexity of the genome-scale metabolic models of human systems, using existing analysis methods is challenging and the parameter selection is not straightforward. Therefore, novel methodological frameworks are necessary for analysing metabolic conditions despite the challenges posed by human models. Particularly, an ongoing challenge has been that of comparing several phenotypes for identifying condition- or disease-specific metabolic signatures. We address this significant challenge by developing a scalable and model-driven approach, ComMet (Comparison of Metabolic states). ComMet enables an in-depth investigation and comparison of metabolic phenotypes in large models while also identifying the underlying functional differences. Novel hypotheses can be generated using ComMet for not only understanding known metabolic phenotypes better but also for guiding the design of new experiments to validate the processes predicted by ComMet.

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.

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

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

Argininosuccinate lyase is a metabolic vulnerability in breast development and cancer

Epithelial-to-mesenchymal transition (EMT) is fundamental to both normal tissue development and cancer progression. We hypothesized that EMT plasticity defines a range of metabolic phenotypes and that individual breast epithelial metabolic phenotypes are likely to fall within this phenotypic landscape. To determine EMT metabolic phenotypes, the metabolism of EMT was described within genome-scale metabolic models (GSMMs) using either transcriptomic or proteomic data from the breast epithelial EMT cell culture model D492. The ability of the different data types to describe breast epithelial metabolism was assessed using constraint-based modeling which was subsequently verified using ¹³C isotope tracer analysis. The application of proteomic data to GSMMs provided relatively higher accuracy in flux predictions compared to the transcriptomic data. Furthermore, the proteomic GSMMs predicted altered cholesterol metabolism and increased dependency on argininosuccinate lyase (ASL) following EMT which were confirmed in vitro using drug assays and siRNA knockdown experiments. The successful verification of the proteomic GSMMs afforded iBreast2886, a breast GSMM that encompasses the metabolic plasticity of EMT as defined by the D492 EMT cell culture model. Analysis of breast tumor proteomic data using iBreast2886 identified vulnerabilities within arginine metabolism that allowed prognostic discrimination of breast cancer patients on a subtype-specific level. Taken together, we demonstrate that the metabolic reconstruction iBreast2886 formalizes the metabolism of breast epithelial cell development and can be utilized as a tool for the functional interpretation of high throughput clinical data.

Advances in Genome-Scale Metabolic Modeling toward Microbial Community Analysis of the Human Microbiome

A genome-scale metabolic model (GEM) represents metabolic pathways of an organism in a mathematical form and can be built using biochemistry and genome annotation data. GEMs are invaluable for understanding organisms since they analyze the metabolic capabilities and behaviors quantitatively and can predict phenotypes. The development of high-throughput data collection techniques led to an immense increase in omics data such as metagenomics, which expand our knowledge on the human microbiome, but this also created a need for systematic analysis of these data. In recent years, GEMs have also been reconstructed for microbial species, including human gut microbiota, and methods for the analysis of microbial communities have been developed to examine the interaction between the organisms or the host. The purpose of this review is to provide a comprehensive guide for the applications of GEMs in microbial community analysis. Starting with GEM repositories, automatic GEM reconstruction tools, and quality control of models, this review will give insights into microbe-microbe and microbe-host interaction predictions and optimization of microbial community models. Recent studies that utilize microbial GEMs and personalized models to infer the influence of microbiota on human diseases such as inflammatory bowel diseases (IBD) or Parkinson's disease are exemplified. Being powerful system biology tools for both species-level and community-level analysis of microbes, GEMs integrated with omics data and machine learning techniques will be indispensable for studying the microbiome and their effects on human physiology as well as for deciphering the mechanisms behind human diseases.

Genome Scale Modeling to Study the Metabolic Competition between Cells in the Tumor Microenvironment

The tumor's physiology emerges from the dynamic interplay of numerous cell types, such as cancer cells, immune cells and stromal cells, within the tumor microenvironment. Immune and cancer cells compete for nutrients within the tumor microenvironment, leading to a metabolic battle between these cell populations. Tumor cells can reprogram their metabolism to meet the high demand of building blocks and ATP for proliferation, and to gain an advantage over the action of immune cells. The study of the metabolic reprogramming mechanisms underlying cancer requires the quantification of metabolic fluxes which can be estimated at the genome-scale with constraint-based or kinetic modeling. Constraint-based models use a set of linear constraints to simulate steady-state metabolic fluxes, whereas kinetic models can simulate both the transient behavior and steady-state values of cellular fluxes and concentrations. The integration of cell- or tissue-specific data enables the construction of context-specific models that reflect cell-type- or tissue-specific metabolic properties. While the available modeling frameworks enable limited modeling of the metabolic crosstalk between tumor and immune cells in the tumor stroma, future developments will likely involve new hybrid kinetic/stoichiometric formulations.

A metabolic modeling approach reveals promising therapeutic targets and antiviral drugs to combat COVID-19

In this study we have developed a method based on Flux Balance Analysis to identify human metabolic enzymes which can be targeted for therapeutic intervention against COVID-19. A literature search was carried out in order to identify suitable inhibitors of these enzymes, which were confirmed by docking calculations. In total, 10 targets and 12 bioactive molecules have been predicted. Among the most promising molecules we identified Triacsin C, which inhibits ACSL3, and which has been shown to be very effective against different viruses, including positive-sense single-stranded RNA viruses. Similarly, we also identified the drug Celgosivir, which has been successfully tested in cells infected with different types of viruses such as Dengue, Zika, Hepatitis C and Influenza. Finally, other drugs targeting enzymes of lipid metabolism, carbohydrate metabolism or protein palmitoylation (such as Propylthiouracil, 2-Bromopalmitate, Lipofermata, Tunicamycin, Benzyl Isothiocyanate, Tipifarnib and Lonafarnib) are also proposed.

Interpreting the lipidome: bioinformatic approaches to embrace the complexity

BACKGROUND

Improvements in mass spectrometry (MS) technologies coupled with bioinformatics developments have allowed considerable advancement in the measurement and interpretation of lipidomics data in recent years. Since research areas employing lipidomics are rapidly increasing, there is a great need for bioinformatic tools that capture and utilize the complexity of the data. Currently, the diversity and complexity within the lipidome is often concealed by summing over or averaging individual lipids up to (sub)class-based descriptors, losing valuable information about biological function and interactions with other distinct lipids molecules, proteins and/or metabolites.

AIM OF REVIEW

To address this gap in knowledge, novel bioinformatics methods are needed to improve identification, quantification, integration and interpretation of lipidomics data. The purpose of this mini-review is to summarize exemplary methods to explore the complexity of the lipidome.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Here we describe six approaches that capture three core focus areas for lipidomics: (1) lipidome annotation including a resolvable database identifier, (2) interpretation via pathway- and enrichment-based methods, and (3) understanding complex interactions to emphasize specific steps in the analytical process and highlight challenges in analyses associated with the complexity of lipidome data.

Lipid and glucose metabolism in white adipocytes: pathways, dysfunction and therapeutics

In mammals, the white adipocyte is a cell type that is specialized for storage of energy (in the form of triacylglycerols) and for energy mobilization (as fatty acids). White adipocyte metabolism confers an essential role to adipose tissue in whole-body homeostasis. Dysfunction in white adipocyte metabolism is a cardinal event in the development of insulin resistance and associated disorders. This Review focuses on our current understanding of lipid and glucose metabolic pathways in the white adipocyte. We survey recent advances in humans on the importance of adipocyte hypertrophy and on the in vivo turnover of adipocytes and stored lipids. At the molecular level, the identification of novel regulators and of the interplay between metabolic pathways explains the fine-tuning between the anabolic and catabolic fates of fatty acids and glucose in different physiological states. We also examine the metabolic alterations involved in the genesis of obesity-associated metabolic disorders, lipodystrophic states, cancers and cancer-associated cachexia. New challenges include defining the heterogeneity of white adipocytes in different anatomical locations throughout the lifespan and investigating the importance of rhythmic processes. Targeting white fat metabolism offers opportunities for improved patient stratification and a wide, yet unexploited, range of therapeutic opportunities.

Molecular pathways behind acquired obesity: Adipose tissue and skeletal muscle multiomics in monozygotic twin pairs discordant for BMI

Tissue-specific mechanisms prompting obesity-related development complications in humans remain unclear. We apply multiomics analyses of subcutaneous adipose tissue and skeletal muscle to examine the effects of acquired obesity among 49 BMI-discordant monozygotic twin pairs. Overall, adipose tissue appears to be more affected by excess body weight than skeletal muscle. In heavier co-twins, we observe a transcriptional pattern of downregulated mitochondrial pathways in both tissues and upregulated inflammatory pathways in adipose tissue. In adipose tissue, heavier co-twins exhibit lower creatine levels; in skeletal muscle, glycolysis- and redox stress-related protein and metabolite levels remain higher. Furthermore, metabolomics analyses in both tissues reveal that several proinflammatory lipids are higher and six of the same lipid derivatives are lower in acquired obesity. Finally, in adipose tissue, but not in skeletal muscle, mitochondrial downregulation and upregulated inflammation are associated with a fatty liver, insulin resistance, and dyslipidemia, suggesting that adipose tissue dominates in acquired obesity.

Multiomics analyses of adipose tissue and skeletal muscle in BMI-discordant twins

Excess body weight downregulates mitochondrial pathways in both tissues

Excess body weight upregulates proinflammatory pathways in both tissues

Adipose tissue alterations are associated with metabolic health in acquired obesity

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.

A dynamic multi-tissue model to study human metabolism

Metabolic modeling enables the study of human metabolism in healthy and in diseased conditions, e.g., the prediction of new drug targets and biomarkers for metabolic diseases. To accurately describe blood and urine metabolite dynamics, the integration of multiple metabolically active tissues is necessary. We developed a dynamic multi-tissue model, which recapitulates key properties of human metabolism at the molecular and physiological level based on the integration of transcriptomics data. It enables the simulation of the dynamics of intra-cellular and extra-cellular metabolites at the genome scale. The predictive capacity of the model is shown through the accurate simulation of different healthy conditions (i.e., during fasting, while consuming meals or during exercise), and the prediction of biomarkers for a set of Inborn Errors of Metabolism with a precision of 83%. This novel approach is useful to prioritize new biomarkers for many metabolic diseases, as well as for the integration of various types of personal omics data, towards the personalized analysis of blood and urine metabolites.

Preclinical techniques to investigate exercise training in vascular pathophysiology

Atherosclerosis is a dynamic process starting with endothelial dysfunction and inflammation and eventually leading to life-threatening arterial plaques. Exercise generally improves endothelial function in a dose-dependent manner by altering hemodynamics, specifically by increased arterial pressure, pulsatility, and shear stress. However, athletes who regularly participate in high-intensity training can develop arterial plaques, suggesting alternative mechanisms through which excessive exercise promotes vascular disease. Understanding the mechanisms that drive atherosclerosis in sedentary versus exercise states may lead to novel rehabilitative methods aimed at improving exercise compliance and physical activity. Preclinical tools, including in vitro cell assays, in vivo animal models, and in silico computational methods, broaden our capabilities to study the mechanisms through which exercise impacts atherogenesis, from molecular maladaptation to vascular remodeling. Here, we describe how preclinical research tools have and can be used to study exercise effects on atherosclerosis. We then propose how advanced bioengineering techniques can be used to address gaps in our current understanding of vascular pathophysiology, including integrating in vitro, in vivo, and in silico studies across multiple tissue systems and size scales. Improving our understanding of the antiatherogenic exercise effects will enable engaging, targeted, and individualized exercise recommendations to promote cardiovascular health rather than treating cardiovascular disease that results from a sedentary lifestyle.

Systems biology approaches to study lipidomes in health and disease

Lipids have many important biological roles, such as energy storage sources, structural components of plasma membranes and as intermediates in metabolic and signaling pathways. Lipid metabolism is under tight homeostatic control, exhibiting spatial and dynamic complexity at multiple levels. Consequently, lipid-related disturbances play important roles in the pathogenesis of most of the common diseases. Lipidomics, defined as the study of lipidomes in biological systems, has emerged as a rapidly-growing field. Due to the chemical and functional diversity of lipids, the application of a systems biology approach is essential if one is to address lipid functionality at different physiological levels. In parallel with analytical advances to measure lipids in biological matrices, the field of computational lipidomics has been rapidly advancing, enabling modeling of lipidomes in their pathway, spatial and dynamic contexts. This review focuses on recent progress in systems biology approaches to study lipids in health and disease, with specific emphasis on methodological advances and biomedical applications.

ReCodLiver0.9: Overcoming Challenges in Genome-Scale Metabolic Reconstruction of a Non-model Species

The availability of genome sequences, annotations, and knowledge of the biochemistry underlying metabolic transformations has led to the generation of metabolic network reconstructions for a wide range of organisms in bacteria, archaea, and eukaryotes. When modeled using mathematical representations, a reconstruction can simulate underlying genotype-phenotype relationships. Accordingly, genome-scale metabolic models (GEMs) can be used to predict the response of organisms to genetic and environmental variations. A bottom-up reconstruction procedure typically starts by generating a draft model from existing annotation data on a target organism. For model species, this part of the process can be straightforward, due to the abundant organism-specific biochemical data. However, the process becomes complicated for non-model less-annotated species. In this paper, we present a draft liver reconstruction, ReCodLiver0.9, of Atlantic cod (Gadus morhua), a non-model teleost fish, as a practicable guide for cases with comparably few resources. Although the reconstruction is considered a draft version, we show that it already has utility in elucidating metabolic response mechanisms to environmental toxicants by mapping gene expression data of exposure experiments to the resulting model.

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

Distinguishing neuropsychiatric disorders is challenging due to the overlap in symptoms and genetic risk factors. People suffering from these disorders face personal and professional challenges. Understanding the dysregulation of brain metabolism under disease condition can aid in effective diagnosis and in developing treatment strategies based on the metabolism. In this study, we reconstructed the metabolic network of three major neuropsychiatric disorders, schizophrenia (SCZ), bipolar disorder (BD) and major depressive disorder (MDD) using transcriptomic data and constrained based modelling approach. We integrated brain transcriptomic data from six independent studies with a recent comprehensive genome-scale metabolic model Recon3D. The analysis of the reconstructed network revealed the flux-level alterations in the peroxisome-mitochondria-golgi axis in neuropsychiatric disorders. We also extracted reporter metabolites and pathways that distinguish these three neuropsychiatric disorders. We found differences with respect to fatty acid oxidation, aromatic and branched chain amino acid metabolism, bile acid synthesis, glycosaminoglycans synthesis and modifications, and phospholipid metabolism. Further, we predicted network perturbations that transform the disease metabolic state to a healthy metabolic state for each disorder. These analyses provide local and global views of the metabolic changes in SCZ, BD and MDD, which may have clinical implications.

Applications of Genome-Wide Screening and Systems Biology Approaches in Drug Repositioning

Modern drug discovery through de novo drug discovery entails high financial costs, low success rates, and lengthy trial periods. Drug repositioning presents a suitable approach for overcoming these issues by re-evaluating biological targets and modes of action of approved drugs. Coupling high-throughput technologies with genome-wide essentiality screens, network analysis, genome-scale metabolic modeling, and machine learning techniques enables the proposal of new drug-target signatures and uncovers unanticipated modes of action for available drugs. Here, we discuss the current issues associated with drug repositioning in light of curated high-throughput multi-omic databases, genome-wide screening technologies, and their application in systems biology/medicine approaches.

Deep learning meets metabolomics: a methodological perspective

Deep learning (DL), an emerging area of investigation in the fields of machine learning and artificial intelligence, has markedly advanced over the past years. DL techniques are being applied to assist medical professionals and researchers in improving clinical diagnosis, disease prediction and drug discovery. It is expected that DL will help to provide actionable knowledge from a variety of 'big data', including metabolomics data. In this review, we discuss the applicability of DL to metabolomics, while presenting and discussing several examples from recent research. We emphasize the use of DL in tackling bottlenecks in metabolomics data acquisition, processing, metabolite identification, as well as in metabolic phenotyping and biomarker discovery. Finally, we discuss how DL is used in genome-scale metabolic modelling and in interpretation of metabolomics data. The DL-based approaches discussed here may assist computational biologists with the integration, prediction and drawing of statistical inference about biological outcomes, based on metabolomics data.

Elucidating potential molecular signatures through host-microbe interactions for reactive arthritis and inflammatory bowel disease using combinatorial approach

Reactive Arthritis (ReA), a rare seronegative inflammatory arthritis, lacks exquisite classification under rheumatic autoimmunity. ReA is solely established using differential clinical diagnosis of the patient cohorts, where pathogenic triggers linked to enteric and urogenital microorganisms e.g. Salmonella, Shigella, Yersinia, Campylobacter, Chlamydia have been reported. Inflammatory Bowel Disease (IBD), an idiopathic enteric disorder co-evolved and attuned to present gut microbiome dysbiosis, can be correlated to the genesis of enteropathic arthropathies like ReA. Gut microbes symbolically modulate immune system homeostasis and are elementary for varied disease patterns in autoimmune disorders. The gut-microbiota axis structured on the core host-microbe interactions execute an imperative role in discerning the etiopathogenesis of ReA and IBD. This study predicts the molecular signatures for ReA with co-evolved IBD through the enveloped host-microbe interactions and microbe-microbe 'interspecies communication', using synonymous gene expression data for selective microbes. We have utilized a combinatorial approach that have concomitant in-silico work-pipeline and experimental validation to corroborate the findings. In-silico analysis involving text mining, metabolic network reconstruction, simulation, filtering, host-microbe interaction, docking and molecular mimicry studies results in robust drug target/s and biomarker/s for co-evolved IBD and ReA. Cross validation of the target/s or biomarker/s was done by targeted gene expression analysis following a non-probabilistic convenience sampling. Studies were performed to substantiate the host-microbe disease network consisting of protein-marker-symptom/disease-pathway-drug associations resulting in possible identification of vital drug targets, biomarkers, pathways and inhibitors for IBD and ReA.Our study identified Na(+)/H(+) anti-porter (NHAA) and Kynureninase (KYNU) to be robust early and essential host-microbe interacting targets for IBD co-evolved ReA. Other vital host-microbe interacting genes, proteins, pathways and drugs include Adenosine Deaminase (ADA), Superoxide Dismutase 2 (SOD2), Catalase (CAT), Angiotensin I Converting Enzyme (ACE), carbon metabolism (folate biosynthesis) and methotrexate. These can serve as potential prognostic/theranostic biomarkers and signatures that can be extrapolated to stratify ReA and related autoimmunity patient cohorts for further pilot studies.

Systems Biology: A Multi-Omics Integration Approach to Metabolism and the Microbiome

The complex and dynamic nature of human physiology, as exemplified by metabolism, has often been overlooked due to the lack of quantitative and systems approaches. Recently, systems biology approaches have pushed the boundaries of our current understanding of complex biochemical, physiological, and environmental interactions, enabling proactive medicine in the near future. From this perspective, we review how state-of-the-art computational modelling of human metabolism, i.e., genome-scale metabolic modelling, could be used to identify the metabolic footprints of diseases, to guide the design of personalized treatments, and to estimate the microbiome contributions to host metabolism. These state-of-the-art models can serve as a scaffold for integrating multi-omics data, thereby enabling the identification of signatures of dysregulated metabolism by systems approaches. For example, increased plasma mannose levels due to decreased uptake in the liver have been identified as a potential biomarker of early insulin resistance by multi-omics approaches. In addition, we also review the emerging axis of human physiology and the human microbiome, discussing its contribution to host metabolism and quantitative approaches to study its variations in individuals.

Unraveling the molecular heterogeneity in type 2 diabetes: a potential subtype discovery followed by metabolic modeling

Background

Type 2 diabetes mellitus (T2DM) is a complex multifactorial disease with a high prevalence worldwide. Insulin resistance and impaired insulin secretion are the two major abnormalities in the pathogenesis of T2DM. Skeletal muscle is responsible for over 75% of the glucose uptake and plays a critical role in T2DM. Here, we sought to provide a better understanding of the abnormalities in this tissue.

Methods

The muscle gene expression patterns were explored in healthy and newly diagnosed T2DM individuals using supervised and unsupervised classification approaches. Moreover, the potential of subtyping T2DM patients was evaluated based on the gene expression patterns.

Results

A machine-learning technique was applied to identify a set of genes whose expression patterns could discriminate diabetic subjects from healthy ones. A gene set comprising of 26 genes was found that was able to distinguish healthy from diabetic individuals with 94% accuracy. In addition, three distinct clusters of diabetic patients with different dysregulated genes and metabolic pathways were identified.

Conclusions

This study indicates that T2DM is triggered by different cellular/molecular mechanisms, and it can be categorized into different subtypes. Subtyping of T2DM patients in combination with their real clinical profiles will provide a better understanding of the abnormalities in each group and more effective therapeutic approaches in the future.

Systems Biology Will Direct Vascular-Targeted Therapy for Obesity

Healthy adipose tissue expansion and metabolism during weight gain require coordinated angiogenesis and lymphangiogenesis. These vascular growth processes rely on the vascular endothelial growth factor (VEGF) family of ligands and receptors (VEGFRs). Several studies have shown that controlling vascular growth by regulating VEGF:VEGFR signaling can be beneficial for treating obesity; however, dysregulated angiogenesis and lymphangiogenesis are associated with several chronic tissue inflammation symptoms, including hypoxia, immune cell accumulation, and fibrosis, leading to obesity-related metabolic disorders. An ideal obesity treatment should minimize adipose tissue expansion and the advent of adverse metabolic consequences, which could be achieved by normalizing VEGF:VEGFR signaling. Toward this goal, a systematic investigation of the interdependency of vascular and metabolic systems in obesity and tools to predict personalized treatment ranges are necessary to improve patient outcomes through vascular-targeted therapies. Systems biology can identify the critical VEGF:VEGFR signaling mechanisms that can be targeted to regress adipose tissue expansion and can predict the metabolic consequences of different vascular-targeted approaches. Establishing a predictive, biologically faithful platform requires appropriate computational models and quantitative tissue-specific data. Here, we discuss the involvement of VEGF:VEGFR signaling in angiogenesis, lymphangiogenesis, adipogenesis, and macrophage specification – key mechanisms that regulate adipose tissue expansion and metabolism. We then provide useful computational approaches for simulating these mechanisms, and detail quantitative techniques for acquiring tissue-specific parameters. Systems biology, through computational models and quantitative data, will enable an accurate representation of obese adipose tissue that can be used to direct the development of vascular-targeted therapies for obesity and associated metabolic disorders.