A community-driven global reconstruction of human metabolism

Big data in IBD: big progress for clinical practice

IBD is a complex multifactorial inflammatory disease of the gut driven by extrinsic and intrinsic factors, including host genetics, the immune system, environmental factors and the gut microbiome. Technological advancements such as next-generation sequencing, high-throughput omics data generation and molecular networks have catalysed IBD research. The advent of artificial intelligence, in particular, machine learning, and systems biology has opened the avenue for the efficient integration and interpretation of big datasets for discovering clinically translatable knowledge. In this narrative review, we discuss how big data integration and machine learning have been applied to translational IBD research. Approaches such as machine learning may enable patient stratification, prediction of disease progression and therapy responses for fine-tuning treatment options with positive impacts on cost, health and safety. We also outline the challenges and opportunities presented by machine learning and big data in clinical IBD research.

Non-coding RNAs, metabolic stress and adaptive mechanisms in cancer

Metabolic reprogramming in cancer describes the multifaceted alterations in metabolism that contribute to tumorigenesis. Major determinants of metabolic phenotypes are the changes in signalling pathways associated with oncogenic activation together with cues from the tumor microenvironment. Therein, depleted oxygen and nutrient levels elicit metabolic stress, requiring cancer cells to engage adaptive mechanisms. Non-coding RNAs (ncRNAs) act as regulatory elements within metabolic pathways and their widespread dysregulation in cancer contributes to altered metabolic phenotypes. Indeed, ncRNAs are the regulatory accomplices of many prominent effectors of metabolic reprogramming including c-MYC and HIFs that are activated by metabolic stress. By example, this review illustrates the range of ncRNAs mechanisms impacting these effectors throughout their DNA-RNA-protein lifecycle along with presenting the mechanistic roles of ncRNAs in adaptive responses to glucose, glutamine and lipid deprivation. We also discuss the facultative activation of metabolic enzymes by ncRNAs, a phenomenon which may reflect a broad but currently invisible level of metabolic regulation. Finally, the translational challenges associated with ncRNA discoveries are discussed, emphasizing the gaps in knowledge together with importance of understanding the molecular basis of ncRNA regulatory mechanisms.

Bioenergetic deficits in Huntington's disease iPSC-derived neural cells and rescue with glycolytic metabolites

Altered cellular metabolism is believed to be an important contributor to pathogenesis of the neurodegenerative disorder Huntington's disease (HD). Research has primarily focused on mitochondrial toxicity, which can cause death of the vulnerable striatal neurons, but other aspects of metabolism have also been implicated. Most previous studies have been carried out using postmortem human brain or non-human cells. Here, we studied bioenergetics in an induced pluripotent stem cell-based model of the disease. We found decreased adenosine triphosphate (ATP) levels in HD cells compared to controls across differentiation stages and protocols. Proteomics data and multiomics network analysis revealed normal or increased levels of mitochondrial messages and proteins, but lowered expression of glycolytic enzymes. Metabolic experiments showed decreased spare glycolytic capacity in HD neurons, while maximal and spare respiratory capacities driven by oxidative phosphorylation were largely unchanged. ATP levels in HD neurons could be rescued with addition of pyruvate or late glycolytic metabolites, but not earlier glycolytic metabolites, suggesting a role for glycolytic deficits as part of the metabolic disturbance in HD neurons. Pyruvate or other related metabolic supplements could have therapeutic benefit in HD.

Systems Biology Approaches to Understand the Host-Microbiome Interactions in Neurodegenerative Diseases

Neurodegenerative diseases (NDDs) comprise a broad range of progressive neurological disorders with multifactorial etiology contributing to disease pathophysiology. Evidence of the microbiome involvement in the gut-brain axis urges the interest in understanding metabolic interactions between the microbiota and host physiology in NDDs. Systems Biology offers a holistic integrative approach to study the interplay between the different biologic systems as part of a whole, and may elucidate the host–microbiome interactions in NDDs. We reviewed direct and indirect pathways through which the microbiota can modulate the bidirectional communication of the gut-brain axis, and explored the evidence of microbial dysbiosis in Alzheimer’s and Parkinson’s diseases. As the gut microbiota being strongly affected by diet, the potential approaches to targeting the human microbiota through diet for the stimulation of neuroprotective microbial-metabolites secretion were described. We explored the potential of Genome-scale metabolic models (GEMs) to infer microbe-microbe and host-microbe interactions and to identify the microbiome contribution to disease development or prevention. Finally, a systemic approach based on GEMs and ‘omics integration, that would allow the design of sustainable personalized anti-inflammatory diets in NDDs prevention, through the modulation of gut microbiota was described.

Nutrient Sensing by Histone Marks: Reading the Metabolic Histone Code Using Tracing, Omics, and Modeling

Several metabolites serve as substrates for histone modifications and communicate changes in the metabolic environment to the epigenome. Technologies such as metabolomics and proteomics have allowed us to reconstruct the interactions between metabolic pathways and histones. These technologies have shed light on how nutrient availability can have a dramatic effect on various histone modifications. This metabolism-epigenome cross talk plays a fundamental role in development, immune function, and diseases like cancer. Yet, major challenges remain in understanding the interactions between cellular metabolism and the epigenome. How the levels and fluxes of various metabolites impact epigenetic marks is still unclear. Discussed herein are recent applications and the potential of systems biology methods such as flux tracing and metabolic modeling to address these challenges and to uncover new metabolic-epigenetic interactions. These systems approaches can ultimately help elucidate how nutrients shape the epigenome of microbes and mammalian cells.

From bag-of-genes to bag-of-genomes: metabolic modelling of communities in the era of metagenome-assembled genomes

Metagenomic sequencing of complete microbial communities has greatly enhanced our understanding of the taxonomic composition of microbiotas. This has led to breakthrough developments in bioinformatic disciplines such as assembly, gene clustering, metagenomic binning of species genomes and the discovery of an incredible, so far undiscovered, taxonomic diversity. However, functional annotations and estimating metabolic processes from single species - or communities - is still challenging. Earlier approaches relied mostly on inferring the presence of key enzymes for metabolic pathways in the whole metagenome, ignoring the genomic context of such enzymes, resulting in the 'bag-of-genes' approach to estimate functional capacities of microbiotas. Here, we review recent developments in metagenomic bioinformatics, with a special focus on emerging technologies to simulate and estimate metabolic information, that can be derived from metagenomic assembled genomes. Genome-scale metabolic models can be used to model the emergent properties of microbial consortia and whole communities, and the progress in this area is reviewed. While this subfield of metagenomics is still in its infancy, it is becoming evident that there is a dire need for further bioinformatic tools to address the complex combinatorial problems in modelling the metabolism of large communities as a 'bag-of-genomes'.

Integration of transcriptomic profile of SARS-CoV-2 infected normal human bronchial epithelial cells with metabolic and protein-protein interaction networks

A novel coronavirus (SARS-CoV-2, formerly known as nCoV-2019) that causes an acute respiratory disease has emerged in Wuhan, China and spread globally in early 2020. On January the 30th, the World Health Organization (WHO) declared spread of this virus as an epidemic and a public health emergency. With its highly contagious characteristic and long incubation time, confinement of SARS-CoV-2 requires drastic lock-down measures to be taken and therefore early diagnosis is crucial. We analysed transcriptome of SARS-CoV-2 infected human lung epithelial cells, compared it with mock-infected cells, used network-based reporter metabolite approach and integrated the transcriptome data with protein-protein interaction network to elucidate the early cellular response. Significantly affected metabolites have the potential to be used in diagnostics while pathways of protein clusters have the potential to be used as targets for supportive or novel therapeutic approaches. Our results are in accordance with the literature on response of IL6 family of cytokines and their importance, in addition, we find that matrix metalloproteinase 2 (MMP2) and matrix metalloproteinase 9 (MMP9) with keratan sulfate synthesis pathway may play a key role in the infection. We hypothesize that MMP9 inhibitors have potential to prevent "cytokine storm" in severely affected patients.

Metabolomics of Endurance Capacity in World Tour Professional Cyclists

The study of elite athletes provides a unique opportunity to define the upper limits of human physiology and performance. Across a variety of sports, these individuals have trained to optimize the physiological parameters of their bodies in order to compete on the world stage. To characterize endurance capacity, techniques such as heart rate monitoring, indirect calorimetry, and whole blood lactate measurement have provided insight into oxygen utilization, and substrate utilization and preference, as well as total metabolic capacity. However, while these techniques enable the measurement of individual, representative variables critical for sports performance, they lack the molecular resolution that is needed to understand which metabolic adaptations are necessary to influence these metrics. Recent advancements in mass spectrometry-based analytical approaches have enabled the measurement of hundreds to thousands of metabolites in a single analysis. Here we employed targeted and untargeted metabolomics approaches to investigate whole blood responses to exercise in elite World Tour (including Tour de France) professional cyclists before and after a graded maximal physiological test. As cyclists within this group demonstrated varying blood lactate accumulation as a function of power output, which is an indicator of performance, we compared metabolic profiles with respect to lactate production to identify adaptations associated with physiological performance. We report that numerous metabolic adaptations occur within this physically elite population (n = 21 males, 28.2 ± 4.7 years old) in association with the rate of lactate accumulation during cycling. Correlation of metabolite values with lactate accumulation has revealed metabolic adaptations that occur in conjunction with improved endurance capacity. In this population, cycling induced increases in tricarboxylic acid (TCA) cycle metabolites and Coenzyme A precursors. These responses occurred proportionally to lactate accumulation, suggesting a link between enhanced mitochondrial networks and the ability to sustain higher workloads. In association with lactate accumulation, altered levels of amino acids before and after exercise point to adaptations that confer unique substrate preference for energy production or to promote more rapid recovery. Cyclists with slower lactate accumulation also have higher levels of basal oxidative stress markers, suggesting long term physiological adaptations in these individuals that support their premier competitive status in worldwide competitions.

Analysis of human metabolism by reducing the complexity of the genome-scale models using redHUMAN

Altered metabolism is associated with many human diseases. Human genome-scale metabolic models (GEMs) were reconstructed within systems biology to study the biochemistry occurring in human cells. However, the complexity of these networks hinders a consistent and concise physiological representation. We present here redHUMAN, a workflow for reconstructing reduced models that focus on parts of the metabolism relevant to a specific physiology using the recently established methods redGEM and lumpGEM. The reductions include the thermodynamic properties of compounds and reactions guaranteeing the consistency of predictions with the bioenergetics of the cell. We introduce a method (redGEMX) to incorporate the pathways used by cells to adapt to the medium. We provide the thermodynamic curation of the human GEMs Recon2 and Recon3D and we apply the redHUMAN workflow to derive leukemia-specific reduced models. The reduced models are powerful platforms for studying metabolic differences between phenotypes, such as diseased and healthy cells.

An analysis of the mechanism underlying photocatalytic disinfection based on integrated metabolic networks and transcriptional data

Photocatalytic disinfection has long been used to combat pathogenic bacteria. However, the specific mechanism underlying photocatalytic disinfection and its corresponding targets remain unclear. In this study, an analysis of the potential mechanism underlying photocatalytic disinfection was performed based on integrated metabolic networks and transcriptional data. Two sets of RNA-seq data (wild type and a photocatalysis-resistant mutant mediated by titanium dioxide (TiO₂)) were processed to constrain the genome scale metabolic models (GSMM) of E. coli. By analyzing the metabolic network, the differential metabolic flux of every reaction was computed in constrained GSMM, and several significantly differential metabolic fluxes in reactions were extracted and analyzed. Most of these reactions were involved in the transmembrane transport of substances and occurred on the inner membrane or were an important component of the cell membrane. These results, which are consistent with the reported information, validated our analysis process. In addition, our work also identified other new and valuable metabolic pathways, such as the reaction ALCD2x, which has a great effect on the energy production process under bacterial anaerobic conditions. The DHAK reaction is also related to the metabolic process of ATP. These reactions with large differential metabolic fluxes merit further research. Additionally, to provide a strategy to address photocatalysis-resistant mutant bacteria, a metabolic compensation analysis was also performed. The metabolic compensation analysis results provided suggestions for a combined method that can effectively combat resistant bacteria. This method could also be used to explore the mechanisms of drug resistance in other microorganisms.

Extending the small-molecule similarity principle to all levels of biology with the Chemical Checker

Small molecules are usually compared by their chemical structure, but there is no unified analytic framework for representing and comparing their biological activity. We present the Chemical Checker (CC), which provides processed, harmonized and integrated bioactivity data on ~800,000 small molecules. The CC divides data into five levels of increasing complexity, from the chemical properties of compounds to their clinical outcomes. In between, it includes targets, off-targets, networks and cell-level information, such as omics data, growth inhibition and morphology. Bioactivity data are expressed in a vector format, extending the concept of chemical similarity to similarity between bioactivity signatures. We show how CC signatures can aid drug discovery tasks, including target identification and library characterization. We also demonstrate the discovery of compounds that reverse and mimic biological signatures of disease models and genetic perturbations in cases that could not be addressed using chemical information alone. Overall, the CC signatures facilitate the conversion of bioactivity data to a format that is readily amenable to machine learning methods.

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease

A multitude of factors contribute to complex diseases and can be measured with ‘omics’ methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources ‘Human metabolism’, ‘Gut microbiome’, ‘Disease’, ‘Nutrition’, and ‘ReconMaps’. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH’s unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

Enzyme capacity-based genome scale modelling of CHO cells

Chinese hamster ovary (CHO) cells are most prevalently used for producing recombinant therapeutics in biomanufacturing. Recently, more rational and systems approaches have been increasingly exploited to identify key metabolic bottlenecks and engineering targets for cell line engineering and process development based on the CHO genome-scale metabolic model which mechanistically characterizes cell culture behaviours. However, it is still challenging to quantify plausible intracellular fluxes and discern metabolic pathway usages considering various clonal traits and bioprocessing conditions. Thus, we newly incorporated enzyme kinetic information into the updated CHO genome-scale model (iCHO2291) and added enzyme capacity constraints within the flux balance analysis framework (ecFBA) to significantly reduce the flux variability in biologically meaningful manner, as such improving the accuracy of intracellular flux prediction. Interestingly, ecFBA could capture the overflow metabolism under the glucose excess condition where the usage of oxidative phosphorylation is limited by the enzyme capacity. In addition, its applicability was successfully demonstrated via a case study where the clone- and media-specific lactate metabolism was deciphered, suggesting that the lactate-pyruvate cycling could be beneficial for CHO cells to efficiently utilize the mitochondrial redox capacity. In summary, iCHO2296 with ecFBA can be used to confidently elucidate cell cultures and effectively identify key engineering targets, thus guiding bioprocess optimization and cell engineering efforts as a part of digital twin model for advanced biomanufacturing in future.

Computational models applied to metabolomics data hints at the relevance of glutamine metabolism in breast cancer

Background

Metabolomics has a great potential in the development of new biomarkers in cancer and it has experiment recent technical advances.

Methods

In this study, metabolomics and gene expression data from 67 localized (stage I to IIIB) breast cancer tumor samples were analyzed, using (1) probabilistic graphical models to define associations using quantitative data without other a priori information; and (2) Flux Balance Analysis and flux activities to characterize differences in metabolic pathways.

Results

On the one hand, both analyses highlighted the importance of glutamine in breast cancer. Moreover, cell experiments showed that treating breast cancer cells with drugs targeting glutamine metabolism significantly affects cell viability. On the other hand, these computational methods suggested some hypotheses and have demonstrated their utility in the analysis of metabolomics data and in associating metabolomics with patient’s clinical outcome.

Conclusions

Computational analyses applied to metabolomics data suggested that glutamine metabolism is a relevant process in breast cancer. Cell experiments confirmed this hypothesis. In addition, these computational analyses allow associating metabolomics data with patient prognosis.

Transcriptome analysis reveals gender-specific differences in overall metabolic response of male and female patients in lung adenocarcinoma

Background

Evidence from multiple studies suggests metabolic abnormalities play an important role in lung cancer. Lung adenocarcinoma (LUAD) is the most common subtype of lung cancer. The present study aimed to explore differences in the global metabolic response between male and female patients in LUAD and to identify the metabolic genes associated with lung cancer susceptibility.

Methods

Transcriptome and clinical LUAD data were acquired from The Cancer Genome Atlas (TCGA) database. Information on metabolic genes and metabolic subsystems were collected from the Recon3D human metabolic model. Two validation datasets (GSE68465 and GSE72094) were downloaded from the Gene Expression Omnibus (GEO) database. Differential expression analysis, gene set enrichment analysis and protein-protein interaction networks were used to identified key metabolic pathways and genes. Functional experiments were used to verify the effects of genes on proliferation, migration, and invasion in lung cancer cells in vitro.

Results

Samples of tumors and adjacent non-tumor tissue from both male and female patients exhibited distinct global patterns of gene expression. In addition, we found large differences in methionine and cysteine metabolism, pyruvate metabolism, cholesterol metabolism, nicotinamide adenine dinucleotide (NAD) metabolism, and nuclear transport between male and female LUAD patients. We identified 34 metabolic genes associated with lung cancer susceptibility in males and 15 in females. Most of the metabolic cancer-susceptibility genes had high prediction accuracy for lung cancer (AUC > 0.9). Furthermore, both bioinformatics analysis and experimental results showed that TAOK2 was down-regulated and ASAH1 was up-regulated in male tumor tissue and female tumor tissue in LUAD. Functional experiments showed that inhibiting ASAH1 suppressed the proliferation, migration, and invasion of lung cancer cells.

Conclusions

Metabolic cancer-susceptibility genes may be used alone or in combination as diagnostic markers for LUAD. Further studies are required to elucidate the functions of these genes in LUAD.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

A census of pathway maps in cancer systems biology

A key goal of cancer systems biology is to use big data to elucidate the molecular networks by which cancer develops. However, to date there has been no systematic evaluation of how far these efforts have progressed. In this Analysis, we survey six major systems biology approaches for mapping and modelling cancer pathways with attention to how well their resulting network maps cover and enhance current knowledge. Our sample of 2,070 systems biology maps captures all literature-curated cancer pathways with significant enrichment, although the strong tendency is for these maps to recover isolated mechanisms rather than entire integrated processes. Systems biology maps also identify previously underappreciated functions, such as a potential role for human papillomavirus-induced chromosomal alterations in ovarian tumorigenesis, and they add new genes to known cancer pathways, such as those related to metabolism, Hippo signalling and immunity. Notably, we find that many cancer networks have been provided only in journal figures and not for programmatic access, underscoring the need to deposit network maps in community databases to ensure they can be readily accessed. Finally, few of these findings have yet been clinically translated, leaving ample opportunity for future translational studies. Periodic surveys of cancer pathway maps, such as the one reported here, are critical to assess progress in the field and identify underserved areas of methodology and cancer biology.

Current Status and Future Prospects of Genome-Scale Metabolic Modeling to Optimize the Use of Mesenchymal Stem Cells in Regenerative Medicine

Mesenchymal stem cells are a promising source for externally grown tissue replacements and patient-specific immunomodulatory treatments. This promise has not yet been fulfilled in part due to production scaling issues and the need to maintain the correct phenotype after re-implantation. One aspect of extracorporeal growth that may be manipulated to optimize cell growth and differentiation is metabolism. The metabolism of MSCs changes during and in response to differentiation and immunomodulatory changes. MSC metabolism may be linked to functional differences but how this occurs and influences MSC function remains unclear. Understanding how MSC metabolism relates to cell function is however important as metabolite availability and environmental circumstances in the body may affect the success of implantation. Genome-scale constraint based metabolic modeling can be used as a tool to fill gaps in knowledge of MSC metabolism, acting as a framework to integrate and understand various data types (e.g., genomic, transcriptomic and metabolomic). These approaches have long been used to optimize the growth and productivity of bacterial production systems and are being increasingly used to provide insights into human health research. Production of tissue for implantation using MSCs requires both optimized production of cell mass and the understanding of the patient and phenotype specific metabolic situation. This review considers the current knowledge of MSC metabolism and how it may be optimized along with the current and future uses of genome scale constraint based metabolic modeling to further this aim.

Fibroblast-specific genome-scale modelling predicts an imbalance in amino acid metabolism in Refsum disease

Refsum disease (RD) is an inborn error of metabolism that is characterised by a defect in peroxisomal α‐oxidation of the branched‐chain fatty acid phytanic acid. The disorder presents with late‐onset progressive retinitis pigmentosa and polyneuropathy and can be diagnosed biochemically by elevated levels of phytanate in plasma and tissues of patients. To date, no cure exists for RD, but phytanate levels in patients can be reduced by plasmapheresis and a strict diet. In this study, we reconstructed a fibroblast‐specific genome‐scale model based on the recently published, FAD‐curated model, based on Recon3D reconstruction. We used transcriptomics (available via GEO database with identifier GSE138379), metabolomics and proteomics (available via ProteomeXchange with identifier PXD015518) data, which we obtained from healthy controls and RD patient fibroblasts incubated with phytol, a precursor of phytanic acid. Our model correctly represents the metabolism of phytanate and displays fibroblast‐specific metabolic functions. Using this model, we investigated the metabolic phenotype of RD at the genome scale, and we studied the effect of phytanate on cell metabolism. We identified 53 metabolites that were predicted to discriminate between healthy and RD patients, several of which with a link to amino acid metabolism. Ultimately, these insights in metabolic changes may provide leads for pathophysiology and therapy.

Databases

Transcriptomics data are available via GEO database with identifier GSE138379, and proteomics data are available via ProteomeXchange with identifier PXD015518.

Improving lipid mapping in Genome Scale Metabolic Networks using ontologies

INTRODUCTION

To interpret metabolomic and lipidomic profiles, it is necessary to identify the metabolic reactions that connect the measured molecules. This can be achieved by putting them in the context of genome-scale metabolic network reconstructions. However, mapping experimentally measured molecules onto metabolic networks is challenging due to differences in identifiers and level of annotation between data and metabolic networks, especially for lipids.

OBJECTIVES

To help linking lipids from lipidomics datasets with lipids in metabolic networks, we developed a new matching method based on the ChEBI ontology. The implementation is freely available as a python library and in MetExplore webserver.

METHODS

Our matching method is more flexible than an exact identifier-based correspondence since it allows establishing a link between molecules even if a different level of precision is provided in the dataset and in the metabolic network. For instance, it can associate a generic class of lipids present in the network with the molecular species detailed in the lipidomics dataset. This mapping is based on the computation of a distance between molecules in ChEBI ontology.

RESULTS

We applied our method to a chemical library (968 lipids) and an experimental dataset (32 modulated lipids) and showed that using ontology-based mapping improves and facilitates the link with genome scale metabolic networks. Beyond network mapping, the results provide ways for improvements in terms of network curation and lipidomics data annotation.

CONCLUSION

This new method being generic, it can be applied to any metabolomics data and therefore improve our comprehension of metabolic modulations.

An atlas of human metabolism

Genome-scale metabolic models (GEMs) are valuable tools to study metabolism and provide a scaffold for the integrative analysis of omics data. Researchers have developed increasingly comprehensive human GEMs, but the disconnect among different model sources and versions impedes further progress. We therefore integrated and extensively curated the most recent human metabolic models to construct a consensus GEM, Human1. We demonstrated the versatility of Human1 through the generation and analysis of cell- and tissue-specific models using transcriptomic, proteomic, and kinetic data. We also present an accompanying web portal, Metabolic Atlas (https://www.metabolicatlas.org/), which facilitates further exploration and visualization of Human1 content. Human1 was created using a version-controlled, open-source model development framework to enable community-driven curation and refinement. This framework allows Human1 to be an evolving shared resource for future studies of human health and disease.

Muscle metabolomics analysis reveals potential biomarkers of exercise‑dependent improvement of the diaphragm function in chronic obstructive pulmonary disease

Decreased diaphragm function is a crucial factor leading to reduced ventilatory efficiency and worsening of quality of life in chronic obstructive pulmonary disease (COPD). Exercise training has been demonstrated to effectively improve the function of the diaphragm. However, the mechanism of this process has not been identified. The emergence of metabolomics has allowed the exploration of new ideas. The present study aimed to analyze the potential biomarkers of exercise-dependent enhancement of diaphragm function in COPD using metabolomics. Sprague Dawley rats were divided into three groups: COPD + exercise group (CEG); COPD model group (CMG); and control group (CG). The first two groups were exposed to cigarette smoke for 16 weeks to establish a COPD model. Then, the rats in the CEG underwent aerobic exercise training for 9 weeks. Following confirmation that exercise effectively improved the diaphragm function, a gas chromatography tandem time-of-flight mass spectrometry analysis system was used to detect the differential metabolites and associated pathways in the diaphragm muscles of the different groups. Following exercise intervention, the pulmonary function and diaphragm contractility of the CEG rats were significantly improved compared with those of the CMG rats. A total of 36 different metabolites were identified in the comparison between the CMG and the CG. Pathway enrichment analysis indicated that these different metabolites were involved in 17 pathways. A total of 29 different metabolites were identified in the comparison between the CMG and the CEG, which are involved in 14 pathways. Candidate biomarkers were selected, and the pathways analysis of these metabolites demonstrated that 2 types of metabolic pathways, the nicotinic acid and nicotinamide metabolism and arginine and proline metabolism pathways, were associated with exercise-induced pulmonary rehabilitation.

Oncogene inference optimization using constraint-based modelling incorporated with protein expression in normal and tumour tissues

Cancer cells are known to exhibit unusual metabolic activity, and yet few metabolic cancer driver genes are known. Genetic alterations and epigenetic modifications of cancer cells result in the abnormal regulation of cellular metabolic pathways that are different when compared with normal cells. Such a metabolic reprogramming can be simulated using constraint-based modelling approaches towards predicting oncogenes. We introduced the tri-level optimization problem to use the metabolic reprogramming towards inferring oncogenes. The algorithm incorporated Recon 2.2 network with the Human Protein Atlas to reconstruct genome-scale metabolic network models of the tissue-specific cells at normal and cancer states, respectively. Such reconstructed models were applied to build the templates of the metabolic reprogramming between normal and cancer cell metabolism. The inference optimization problem was formulated to use the templates as a measure towards predicting oncogenes. The nested hybrid differential evolution algorithm was applied to solve the problem to overcome solving difficulty for transferring the inner optimization problem into the single one. Head and neck squamous cells were applied as a case study to evaluate the algorithm. We detected 13 of the top-ranked one-hit dysregulations and 17 of the top-ranked two-hit oncogenes with high similarity ratios to the templates. According to the literature survey, most inferred oncogenes are consistent with the observation in various tissues. Furthermore, the inferred oncogenes were highly connected with the TP53/AKT/IGF/MTOR signalling pathway through PTEN, which is one of the most frequently detected tumour suppressor genes in human cancer.

Spatial modeling of prostate cancer metabolic gene expression reveals extensive heterogeneity and selective vulnerabilities

Spatial heterogeneity is a fundamental feature of the tumor microenvironment (TME), and tackling spatial heterogeneity in neoplastic metabolic aberrations is critical for tumor treatment. Genome-scale metabolic network models have been used successfully to simulate cancer metabolic networks. However, most models use bulk gene expression data of entire tumor biopsies, ignoring spatial heterogeneity in the TME. To account for spatial heterogeneity, we performed spatially-resolved metabolic network modeling of the prostate cancer microenvironment. We discovered novel malignant-cell-specific metabolic vulnerabilities targetable by small molecule compounds. We predicted that inhibiting the fatty acid desaturase SCD1 may selectively kill cancer cells based on our discovery of spatial separation of fatty acid synthesis and desaturation. We also uncovered higher prostaglandin metabolic gene expression in the tumor, relative to the surrounding tissue. Therefore, we predicted that inhibiting the prostaglandin transporter SLCO2A1 may selectively kill cancer cells. Importantly, SCD1 and SLCO2A1 have been previously shown to be potently and selectively inhibited by compounds such as CAY10566 and suramin, respectively. We also uncovered cancer-selective metabolic liabilities in central carbon, amino acid, and lipid metabolism. Our novel cancer-specific predictions provide new opportunities to develop selective drug targets for prostate cancer and other cancers where spatial transcriptomics datasets are available.

An adjustable algal chloroplast plug-and-play model for genome-scale metabolic models

The chloroplast is a central part of plant cells, as this is the organelle where the photosynthesis, fixation of inorganic carbon, and other key functions related to fatty acid synthesis and amino acid synthesis occur. Since this organelle should be an integral part of any genome-scale metabolic model for a microalgae or a higher plant, it is of great interest to generate a detailed and standardized chloroplast model. Additionally, we see the need for a novel type of sub-model template, or organelle model, which could be incorporated into a larger, less specific genome-scale metabolic model, while allowing for minor differences between chloroplast-containing organisms. The result of this work is the very first standardized chloroplast model, iGR774, consisting of 788 reactions, 764 metabolites, and 774 genes. The model is currently able to run in three different modes, mimicking the chloroplast metabolism of three photosynthetic microalgae-Nannochloropsis gaditana, Chlamydomonas reinhardtii and Phaeodactylum tricornutum. In addition to developing the chloroplast metabolic network reconstruction, we have developed multiple software tools for working with this novel type of sub-model in the COBRA Toolbox for MATLAB, including tools for connecting the chloroplast model to a genome-scale metabolic reconstruction in need of a chloroplast, for switching the model between running in different organism modes, and for expanding it by introducing more reactions either related to one of the current organisms included in the model, or to a new organism.

FastMM: an efficient toolbox for personalized constraint-based metabolic modeling

Background

Constraint-based metabolic modeling has been applied to understand metabolism related disease mechanisms, to predict potential new drug targets and anti-metabolites, and to identify biomarkers of complex diseases. Although the state-of-art modeling toolbox, COBRA 3.0, is powerful, it requires substantial computing time conducting flux balance analysis, knockout analysis, and Markov Chain Monte Carlo (MCMC) sampling, which may limit its application in large scale genome-wide analysis.

Results

Here, we rewrote the underlying code of COBRA 3.0 using C/C++, and developed a toolbox, termed FastMM, to effectively conduct constraint-based metabolic modeling. The results showed that FastMM is 2~400 times faster than COBRA 3.0 in performing flux balance analysis and knockout analysis and returns consistent outputs. When applied to MCMC sampling, FastMM is 8 times faster than COBRA 3.0. FastMM is also faster than some efficient metabolic modeling applications, such as Cobrapy and Fast-SL. In addition, we developed a Matlab/Octave interface for fast metabolic modeling. This interface was fully compatible with COBRA 3.0, enabling users to easily perform complex applications for metabolic modeling. For example, users who do not have deep constraint-based metabolic model knowledge can just type one command in Matlab/Octave to perform personalized metabolic modeling. Users can also use the advance and multiple threading parameters for complex metabolic modeling. Thus, we provided an efficient and user-friendly solution to perform large scale genome-wide metabolic modeling. For example, FastMM can be applied to the modeling of individual cancer metabolic profiles of hundreds to thousands of samples in the Cancer Genome Atlas (TCGA).

Conclusion

FastMM is an efficient and user-friendly toolbox for large-scale personalized constraint-based metabolic modeling. It can serve as a complementary and invaluable improvement to the existing functionalities in COBRA 3.0. FastMM is under GPL license and can be freely available at GitHub site: https://github.com/GonghuaLi/FastMM.

Exploring Integrative Analysis Using the BioMedical Evidence Graph

PURPOSE

The analysis of cancer biology data involves extremely heterogeneous data sets, including information from RNA sequencing, genome-wide copy number, DNA methylation data reporting on epigenetic regulation, somatic mutations from whole-exome or whole-genome analyses, pathology estimates from imaging sections or subtyping, drug response or other treatment outcomes, and various other clinical and phenotypic measurements. Bringing these different resources into a common framework, with a data model that allows for complex relationships as well as dense vectors of features, will unlock integrated data set analysis.

METHODS

We introduce the BioMedical Evidence Graph (BMEG), a graph database and query engine for discovery and analysis of cancer biology. The BMEG is unique from other biologic data graphs in that sample-level molecular and clinical information is connected to reference knowledge bases. It combines gene expression and mutation data with drug-response experiments, pathway information databases, and literature-derived associations.

RESULTS

The construction of the BMEG has resulted in a graph containing > 41 million vertices and 57 million edges. The BMEG system provides a graph query-based application programming interface to enable analysis, with client code available for Python, Javascript, and R, and a server online at bmeg.io. Using this system, we have demonstrated several forms of cross-data set analysis to show the utility of the system.

CONCLUSION

The BMEG is an evolving resource dedicated to enabling integrative analysis. We have demonstrated queries on the system that illustrate mutation significance analysis, drug-response machine learning, patient-level knowledge-base queries, and pathway level analysis. We have compared the resulting graph to other available integrated graph systems and demonstrated the former is unique in the scale of the graph and the type of data it makes available.

Computational Modeling of the Human Microbiome

The impact of microorganisms on human health has long been acknowledged and studied, but recent advances in research methodologies have enabled a new systems-level perspective on the collections of microorganisms associated with humans, the human microbiome. Large-scale collaborative efforts such as the NIH Human Microbiome Project have sought to kick-start research on the human microbiome by providing foundational information on microbial composition based upon specific sites across the human body. Here, we focus on the four main anatomical sites of the human microbiome: gut, oral, skin, and vaginal, and provide information on site-specific background, experimental data, and computational modeling. Each of the site-specific microbiomes has unique organisms and phenomena associated with them; there are also high-level commonalities. By providing an overview of different human microbiome sites, we hope to provide a perspective where detailed, site-specific research is needed to understand causal phenomena that impact human health, but there is equally a need for more generalized methodology improvements that would benefit all human microbiome research.

Multiple Pathways Involved in Palmitic Acid-Induced Toxicity: A System Biology Approach

Inflammation is a complex biological response to injuries, metabolic disorders or infections. In the brain, astrocytes play an important role in the inflammatory processes during neurodegenerative diseases. Recent studies have shown that the increase of free saturated fatty acids such as palmitic acid produces a metabolic inflammatory response in astrocytes generally associated with damaging mechanisms such as oxidative stress, endoplasmic reticulum stress, and autophagic defects. In this aspect, the synthetic neurosteroid tibolone has shown to exert protective functions against inflammation in neuronal experimental models without the tumorigenic effects exerted by sexual hormones such as estradiol and progesterone. However, there is little information regarding the specific mechanisms of tibolone in astrocytes during inflammatory insults. In the present study, we performed a genome-scale metabolic reconstruction of astrocytes that was used to study astrocytic response during an inflammatory insult by palmitate through Flux Balance Analysis methods and data mining. In this aspect, we assessed the metabolic fluxes of human astrocytes under three different scenarios: healthy (normal conditions), induced inflammation by palmitate, and tibolone treatment under palmitate inflammation. Our results suggest that tibolone reduces the L-glutamate-mediated neurotoxicity in astrocytes through the modulation of several metabolic pathways involved in glutamate uptake. We also identified a set of reactions associated with the protective effects of tibolone, including the upregulation of taurine metabolism, gluconeogenesis, cPPAR and the modulation of calcium signaling pathways. In conclusion, the different scenarios studied in our model allowed us to identify several metabolic fluxes perturbed under an inflammatory response and the protective mechanisms exerted by tibolone.

Activated Drp1-mediated mitochondrial ROS influence the gut microbiome and intestinal barrier after hemorrhagic shock

A role of the mitochondrial dynamin-related protein (Drp1) on gut microbiome composition and intestinal barrier function after hemorrhagic shock has not been identified previously and thus addressed in this study. Here, we used a combination of 16S rRNA gene sequencing and mass spectrometry-based metabolomics profiling in WT and Drp1 KO mouse models to examine the functional impact of activated Drp1 on the gut microbiome as well as mitochondrial metabolic regulation after hemorrhagic shock. Our data showed that changes in mitochondrial Drp1 activity participated in the regulation of intestinal barrier function after hemorrhagic shock. Activated Drp1 significantly perturbed gut microbiome composition in the Bacteroidetes phylum. The abundance of short-chain fatty acid (SCFA) producing microbes, such as Bacteroides, Butyricimonas and Odoribacter, was markedly decreased in mice after shock, and was inversely correlated with both the distribution of the tight junction protein ZO1 and intestinal permeability. Together, these data suggest that Drp1 activation perturbs the gut microbiome community and SCFA production in a ROS-specific manner and thereby substantially disturbs tight junctions and intestinal barrier function after hemorrhagic shock. Our findings provide novel insights for targeting Drp1-mediated mitochondrial function as well as the microbiome in the treatment of intestinal barrier dysfunction after shock.

Network-Based Metabolism-Centered Screening of Potential Drug Targets in Klebsiella pneumoniae at Genome Scale

Klebsiella pneumoniae is an opportunistic bacterial pathogen leading to life-threatening nosocomial infections. Emergence of highly resistant strains poses a major challenge in the management of the infections by healthcare-associated K. pneumoniae isolates. Thus, despite intensive efforts, the current treatment strategies remain insufficient to eradicate such infections. Failure of the conventional infection-prevention and treatment efforts explicitly indicates the requirement of new therapeutic approaches. This prompted us to systematically analyze the K. pneumoniae metabolism to investigate drug targets. Genome-scale metabolic networks (GMNs) facilitating the systematic analysis of the metabolism are promising platforms. Thus, we used a GMN of K. pneumoniae MGH 78578 to determine putative targets through gene- and metabolite-centric approaches. To develop more realistic infection models, we performed the bacterial growth simulations within different host-mimicking media, using an improved biomass formation reaction. We selected more suitable targets based on several property-based prioritization procedures. KdsA was identified as the high-ranked putative target satisfying most of the target prioritization criteria specified under the gene-centric approach. Through a structure-based virtual screening protocol, we identified potential KdsA inhibitors. In addition, the metabolite-centric approach extended the drug target list based on synthetic lethality. This revealed the importance of combined metabolic analyses for a better understanding of the metabolism. To our knowledge, this is the first comprehensive effort on the investigation of the K. pneumoniae metabolism for drug target prediction through the constraint-based analysis of its GMN in conjunction with several bioinformatic approaches. This study can guide the researchers for the future drug designs by providing initial findings regarding crucial components of the Klebsiella metabolism.

A scalable method for parameter-free simulation and validation of mechanistic cellular signal transduction network models

The metabolic modelling community has established the gold standard for bottom-up systems biology with reconstruction, validation and simulation of mechanistic genome-scale models. Similar methods have not been established for signal transduction networks, where the representation of complexes and internal states leads to scalability issues in both model formulation and execution. While rule- and agent-based methods allow efficient model definition and execution, respectively, model parametrisation introduces an additional layer of uncertainty due to the sparsity of reliably measured parameters. Here, we present a scalable method for parameter-free simulation of mechanistic signal transduction networks. It is based on rxncon and uses a bipartite Boolean logic with separate update rules for reactions and states. Using two generic update rules, we enable translation of any rxncon model into a unique Boolean model, which can be used for network validation and simulation-allowing the prediction of system-level function directly from molecular mechanistic data. Through scalable model definition and simulation, and the independence of quantitative parameters, it opens up for simulation and validation of mechanistic genome-scale models of signal transduction networks.

Pathway Commons 2019 Update: integration, analysis and exploration of pathway data

Pathway Commons (https://www.pathwaycommons.org) is an integrated resource of publicly available information about biological pathways including biochemical reactions, assembly of biomolecular complexes, transport and catalysis events and physical interactions involving proteins, DNA, RNA, and small molecules (e.g. metabolites and drug compounds). Data is collected from multiple providers in standard formats, including the Biological Pathway Exchange (BioPAX) language and the Proteomics Standards Initiative Molecular Interactions format, and then integrated. Pathway Commons provides biologists with (i) tools to search this comprehensive resource, (ii) a download site offering integrated bulk sets of pathway data (e.g. tables of interactions and gene sets), (iii) reusable software libraries for working with pathway information in several programming languages (Java, R, Python and Javascript) and (iv) a web service for programmatically querying the entire dataset. Visualization of pathways is supported using the Systems Biological Graphical Notation (SBGN). Pathway Commons currently contains data from 22 databases with 4794 detailed human biochemical processes (i.e. pathways) and ∼2.3 million interactions. To enhance the usability of this large resource for end-users, we develop and maintain interactive web applications and training materials that enable pathway exploration and advanced analysis.

Comparative Metabolic Network Flux Analysis to Identify Differences in Cellular Metabolism

Metabolic network flux analysis uses genome-scale metabolic reconstructions to integrate transcriptomics, proteomics, and/or metabolomics data to allow for comprehensive interpretation of genotype to metabolic phenotype relationships. The compilation of many Constraint-based model analysis methods into one MATLAB package, the COBRAtoolbox, has opened the possibility of using these methods to the many biologists with some knowledge of the commonly used statistical program, MATLAB. Here we outline the steps required to take a published genome-scale metabolic reconstruction and interrogate its consistency and biological feasibility. Subsequently, we demonstrate how mRNA expression data and metabolomics data, relating to one or more cell types or biological contexts, can be applied to constrain and generate metabolic models descriptive of metabolic flux phenotypes. Finally, we describe the comparison of the resulting models and model outputs with the aim of identifying metabolic biomarkers and changes in cellular metabolism.

Pathway Analysis for Targeted and Untargeted Metabolomics

Recent advances in analytical techniques, particularly LC-MS, generate increasingly large and complex metabolomics datasets. Pathway analysis tools help place the experimental observations into relevant biological or disease context. This chapter provides an overview of the general concepts and common tools for pathway analysis, including Mummichog for untargeted metabolomics. Examples of pathway mapping, MetScape, and Mummichog are explained. This serves as both a practical tutorial and a timely survey of pathway analysis for label-free metabolomics data.

Metabolic Network Reconstructions to Predict Drug Targets and Off-Target Effects

The drug development pipeline has stalled because of the difficulty in identifying new drug targets while minimizing off-target effects. Computational methods, such as the use of metabolic network reconstructions, may provide a cost-effective platform to test new hypotheses for drug targets and prevent off-target effects. Here, we summarize available methods to identify drug targets and off-target effects using either reaction-centric, gene-centric, or metabolite-centric approaches with genome-scale metabolic network reconstructions.

Compartment and hub definitions tune metabolic networks for metabolomic interpretations

BACKGROUND

Metabolic networks represent all chemical reactions that occur between molecular metabolites in an organism's cells. They offer biological context in which to integrate, analyze, and interpret omic measurements, but their large scale and extensive connectivity present unique challenges. While it is practical to simplify these networks by placing constraints on compartments and hubs, it is unclear how these simplifications alter the structure of metabolic networks and the interpretation of metabolomic experiments.

RESULTS

We curated and adapted the latest systemic model of human metabolism and developed customizable tools to define metabolic networks with and without compartmentalization in subcellular organelles and with or without inclusion of prolific metabolite hubs. Compartmentalization made networks larger, less dense, and more modular, whereas hubs made networks larger, more dense, and less modular. When present, these hubs also dominated shortest paths in the network, yet their exclusion exposed the subtler prominence of other metabolites that are typically more relevant to metabolomic experiments. We applied the non-compartmental network without metabolite hubs in a retrospective, exploratory analysis of metabolomic measurements from 5 studies on human tissues. Network clusters identified individual reactions that might experience differential regulation between experimental conditions, several of which were not apparent in the original publications.

CONCLUSIONS

Exclusion of specific metabolite hubs exposes modularity in both compartmental and non-compartmental metabolic networks, improving detection of relevant clusters in omic measurements. Better computational detection of metabolic network clusters in large data sets has potential to identify differential regulation of individual genes, transcripts, and proteins.

Proteomic and Transcriptomic Changes in Hibernating Grizzly Bears Reveal Metabolic and Signaling Pathways that Protect against Muscle Atrophy

Muscle atrophy is a physiological response to disuse and malnutrition, but hibernating bears are largely resistant to this phenomenon. Unlike other mammals, they efficiently reabsorb amino acids from urine, periodically activate muscle contraction, and their adipocytes differentially responds to insulin. The contribution of myocytes to the reduced atrophy remains largely unknown. Here we show how metabolism and atrophy signaling are regulated in skeletal muscle of hibernating grizzly bear. Metabolic modeling of proteomic changes suggests an autonomous increase of non-essential amino acids (NEAA) in muscle and treatment of differentiated myoblasts with NEAA is sufficient to induce hypertrophy. Our comparison of gene expression in hibernation versus muscle atrophy identified several genes differentially regulated during hibernation, including Pdk4 and Serpinf1. Their trophic effects extend to myoblasts from non-hibernating species (including C. elegans), as documented by a knockdown approach. Together, these changes reflect evolutionary favored adaptations that, once translated to the clinics, could help improve atrophy treatment.

Genome-Scale Metabolic Modeling with Protein Expressions of Normal and Cancerous Colorectal Tissues for Oncogene Inference

Although cancer has historically been regarded as a cell proliferation disorder, it has recently been considered a metabolic disease. The first discovery of metabolic alterations in cancer cells refers to Otto Warburg’s observations. Cancer metabolism results in alterations in metabolic fluxes that are evident in cancer cells compared with most normal tissue cells. This study applied protein expressions of normal and cancer cells to reconstruct two tissue-specific genome-scale metabolic models. Both models were employed in a tri-level optimization framework to infer oncogenes. Moreover, this study also introduced enzyme pseudo-coding numbers in the gene association expression to avoid performing posterior decision-making that is necessary for the reaction-based method. Colorectal cancer (CRC) was the topic of this case study, and 20 top-ranked oncogenes were determined. Notably, these dysregulated genes were involved in various metabolic subsystems and compartments. We found that the average similarity ratio for each dysregulation is higher than 98%, and the extent of similarity for flux changes is higher than 93%. On the basis of surveys of PubMed and GeneCards, these oncogenes were also investigated in various carcinomas and diseases. Most dysregulated genes connect to catalase that acts as a hub and connects protein signaling pathways, such as those involving TP53, mTOR, AKT1, MAPK1, EGFR, MYC, CDK8, and RAS family.

Bioenergetics and translational metabolism: implications for genetics, physiology and precision medicine

It is now becoming clear that human metabolism is extremely plastic and varies substantially between healthy individuals. Understanding the biochemistry that underlies this physiology will enable personalized clinical interventions related to metabolism. Mitochondrial quality control and the detailed mechanisms of mitochondrial energy generation are central to understanding susceptibility to pathologies associated with aging including cancer, cardiac and neurodegenerative diseases. A precision medicine approach is also needed to evaluate the impact of exercise or caloric restriction on health. In this review, we discuss how technical advances in assessing mitochondrial genetics, cellular bioenergetics and metabolomics offer new insights into developing metabolism-based clinical tests and metabolotherapies. We discuss informatics approaches, which can define the bioenergetic-metabolite interactome and how this can help define healthy energetics. We propose that a personalized medicine approach that integrates metabolism and bioenergetics with physiologic parameters is central for understanding the pathophysiology of diseases with a metabolic etiology. New approaches that measure energetics and metabolomics from cells isolated from human blood or tissues can be of diagnostic and prognostic value to precision medicine. This is particularly significant with the development of new metabolotherapies, such as mitochondrial transplantation, which could help treat complex metabolic diseases.

n-Butylamine for Improving the Efficiency of Untargeted Mass Spectrometry Analysis of Plasma Metabolite Composition

A comparative study of the impact of n-butylamine and traditionally used additives (ammonium hydroxide and formic acid) on the efficiency of the electrospray ionization (ESI) process for the enhancement of metabolite coverage was performed by direct injection mass spectrometry (MS) analysis in negative mode. Evaluation of obtained MS data showed that n-butylamine is one of the most effective additives for the analysis of metabolite composition in ESI in negative ion mode (ESI(-)) The limitations of the use of n-butylamine and other alkylamines in the analysis of metabolic composition and a decontamination procedure that can reduce MS device contamination after their application are discussed. The proposed procedure allows the performance of high-sensitivity analysis of low-molecular-weight compounds on the same MS device in both polarities.

Network Modeling Approaches and Applications to Unravelling Non-Alcoholic Fatty Liver Disease

Non-alcoholic fatty liver disease (NAFLD) is a progressive condition of the liver encompassing a range of pathologies including steatosis, non-alcoholic steatohepatitis (NASH), cirrhosis, and hepatocellular carcinoma. Research into this disease is imperative due to its rapid growth in prevalence, economic burden, and current lack of FDA approved therapies. NAFLD involves a highly complex etiology that calls for multi-tissue multi-omics network approaches to uncover the pathogenic genes and processes, diagnostic biomarkers, and potential therapeutic strategies. In this review, we first present a basic overview of disease pathogenesis, risk factors, and remaining knowledge gaps, followed by discussions of the need and concepts of multi-tissue multi-omics approaches, various network methodologies and application examples in NAFLD research. We highlight the findings that have been uncovered thus far including novel biomarkers, genes, and biological pathways involved in different stages of NAFLD, molecular connections between NAFLD and its comorbidities, mechanisms underpinning sex differences, and druggable targets. Lastly, we outline the future directions of implementing network approaches to further improve our understanding of NAFLD in order to guide diagnosis and therapeutics.

MetaboRank: network-based recommendation system to interpret and enrich metabolomics results

Motivation

Metabolomics has shown great potential to improve the understanding of complex diseases, potentially leading to therapeutic target identification. However, no single analytical method allows monitoring all metabolites in a sample, resulting in incomplete metabolic fingerprints. This incompleteness constitutes a stumbling block to interpretation, raising the need for methods that can enrich those fingerprints. We propose MetaboRank, a new solution inspired by social network recommendation systems for the identification of metabolites potentially related to a metabolic fingerprint.

Results

MetaboRank method had been used to enrich metabolomics data obtained on cerebrospinal fluid samples from patients suffering from hepatic encephalopathy (HE). MetaboRank successfully recommended metabolites not present in the original fingerprint. The quality of recommendations was evaluated by using literature automatic search, in order to check that recommended metabolites could be related to the disease. Complementary mass spectrometry experiments and raw data analysis were performed to confirm these suggestions. In particular, MetaboRank recommended the overlooked α-ketoglutaramate as a metabolite which should be added to the metabolic fingerprint of HE, thus suggesting that metabolic fingerprints enhancement can provide new insight on complex diseases.

Availability and implementation

Method is implemented in the MetExplore server and is available at www.metexplore.fr. A tutorial is available at https://metexplore.toulouse.inra.fr/com/tutorials/MetaboRank/2017-MetaboRank.pdf.

Supplementary information

Supplementary data are available at Bioinformatics online.

Comprehensive Flux Modeling of Chlamydia trachomatis Proteome and qRT-PCR Data Indicate Biphasic Metabolic Differences Between Elementary Bodies and Reticulate Bodies During Infection

Metabolic adaptation to the host cell is important for obligate intracellular pathogens such as Chlamydia trachomatis (Ct). Here we infer the flux differences for Ct from proteome and qRT-PCR data by comprehensive pathway modeling. We compare the comparatively inert infectious elementary body (EB) and the active replicative reticulate body (RB) systematically using a genome-scale metabolic model with 321 metabolites and 277 reactions. This did yield 84 extreme pathways based on a published proteomics dataset at three different time points of infection. Validation of predictions was done by quantitative RT-PCR of enzyme mRNA expression at three time points. Ct’s major active pathways are glycolysis, gluconeogenesis, glycerol-phospholipid (GPL) biosynthesis (support from host acetyl-CoA) and pentose phosphate pathway (PPP), while its incomplete TCA and fatty acid biosynthesis are less active. The modeled metabolic pathways are much more active in RB than in EB. Our in silico model suggests that EB and RB utilize folate to generate NAD(P)H using independent pathways. The only low metabolic flux inferred for EB involves mainly carbohydrate metabolism. RB utilizes energy -rich compounds to generate ATP in nucleic acid metabolism. Validation data for the modeling include proteomics experiments (model basis) as well as qRT-PCR confirmation of selected metabolic enzyme mRNA expression differences. The metabolic modeling is made fully available here. Its detailed insights and models on Ct metabolic adaptations during infection are a useful modeling basis for future studies.

Identification of Cancer-associated metabolic vulnerabilities by modeling multi-objective optimality in metabolism

BACKGROUND

Cancer cells undergo global reprogramming of cellular metabolism to satisfy demands of energy and biomass during proliferation and metastasis. Computational modeling of genome-scale metabolic models is an effective approach for designing new therapeutics targeting dysregulated cancer metabolism by identifying metabolic enzymes crucial for satisfying metabolic goals of cancer cells, but nearly all previous studies neglect the existence of metabolic demands other than biomass synthesis and trade-offs between these contradicting metabolic demands. It is thus necessary to develop computational models covering multiple metabolic objectives to study cancer metabolism and identify novel metabolic targets.

METHODS

We developed a multi-objective optimization model for cancer cell metabolism at genome-scale and an integrated, data-driven workflow for analyzing the Pareto optimality of this model in achieving multiple metabolic goals and identifying metabolic enzymes crucial for maintaining cancer-associated metabolic phenotypes. Using this workflow, we constructed cell line-specific models for a panel of cancer cell lines and identified lists of metabolic targets promoting or suppressing cancer cell proliferation or the Warburg Effect. The targets were then validated using knockdown and over-expression experiments in cultured cancer cell lines.

RESULTS

We found that the multi-objective optimization model correctly predicted phenotypes including cell growth rates, essentiality of metabolic genes and cell line specific sensitivities to metabolic perturbations. To our surprise, metabolic enzymes promoting proliferation substantially overlapped with those suppressing the Warburg Effect, suggesting that simply targeting the overlapping enzymes may lead to complicated outcomes. We also identified lists of metabolic enzymes important for maintaining rapid proliferation or high Warburg Effect while having little effect on the other. The importance of these enzymes in cancer metabolism predicted by the model was validated by their association with cancer patient survival and knockdown and overexpression experiments in a variety of cancer cell lines.

CONCLUSIONS

These results confirm this multi-objective optimization model as a novel and effective approach for studying trade-off between metabolic demands of cancer cells and identifying cancer-associated metabolic vulnerabilities, and suggest novel metabolic targets for cancer treatment.