Tissue-specific Metabolic Models Research Articles

New workflow predicts drug targets against SARS-CoV-2 via metabolic changes in infected cells.

COVID-19 is one of the deadliest respiratory diseases, and its emergence caught the pharmaceutical industry off guard. While vaccines have been rapidly developed, treatment options for infected people remain scarce, and COVID-19 poses a substantial global threat. This study presents a novel workflow to predict robust druggable targets against emerging RNA viruses using metabolic networks and information of the viral structure and its genome sequence. For this purpose, we implemented pymCADRE and PREDICATE to create tissue-specific metabolic models, construct viral biomass functions and predict host-based antiviral targets from more than one genome. We observed that pymCADRE reduces the computational time of flux variability analysis for internal optimizations. We applied these tools to create a new metabolic network of primary bronchial epithelial cells infected with SARS-CoV-2 and identified enzymatic reactions with inhibitory effects. The most promising reported targets were from the purine metabolism, while targeting the pyrimidine and carbohydrate metabolisms seemed to be promising approaches to enhance viral inhibition. Finally, we computationally tested the robustness of our targets in all known variants of concern, verifying our targets' inhibitory effects. Since laboratory tests are time-consuming and involve complex readouts to track processes, our workflow focuses on metabolic fluxes within infected cells and is applicable for rapid hypothesis-driven identification of potentially exploitable antivirals concerning various viruses and host cell types.

Metabolic modeling of host\u2013microbe interactions for therapeutics in colorectal cancer

The onset of colorectal cancer (CRC) is often attributed to gut bacterial dysbiosis, and thus gut microbiota are highly relevant in devising treatment strategies. Certain gut microbes, like Enterococcus spp., exhibit remarkable anti-neoplastic and probiotic properties, which can aid in silver nanoparticle (AgNPs) induced reactive oxygen species (ROS)-based CRC treatment. However, the effects of AgNPs on gut microbial metabolism have not been reported thus far. In this study, a detailed systems-level understanding of ROS metabolism in Enterococcus durans (E. durans), a representative gut microbe, was gained using constraint-based modeling, wherein, the critical association between ROS and folate metabolism was established. Experimental studies involving low AgNP concentration treatment of E. durans cultures confirmed these modeling predictions (an increased extracellular folate concentration by 52%, at the 9th h of microbial growth, was observed). Besides, the computational studies established various metabolic pathways involving amino acids, energy metabolites, nucleotides, and SCFAs as the key players in elevating folate levels on ROS exposure. The anti-cancer potential of E. durans was also studied through MTT analysis of HCT 116 cells treated with microbial culture (AgNP treated) supernatant. A decrease in cell viability by 19% implicated the role of microbial metabolites (primarily folate) in causing cell death. The genome-scale modeling approach was then extended to extensively model CRC metabolism, as well as CRC–E. durans interactions in the context of CRC treatment, using tissue-specific metabolic models of CRC and healthy colon. These findings on further validation can facilitate the development of robust and effective cancer therapy.

Can we capture an accurate view of tissue metabolism from an expression profile?

Genome-scale models of metabolism have been constructed with the aim of illuminating the molecular basis of cell phenotypes. Since some enzymes are only active in specific cell or tissue types, several algorithms use various omics data types (e.g., transcriptomics, proteomics, metabolomics, genomic variant calls) to construct cell line and tissue-specific metabolic models from genome-scale models. However, it is unclear how algorithm and parameter selection (e.g., gene expression thresholds, metabolic constraints) affects model content and predictive accuracy. To investigate this, we built hundreds of models of four different cancer cell lines using transcriptomics, metabolomics, various algorithms and different parameter sets. Model content varied substantially across different parameter sets, but all algorithms generally increased accuracy in identifying cancer cell vulnerabilities. However, the choice of algorithm had the largest impact on model accuracy. Finally, we further developed a framework to predict, de novo, the metabolic capabilities of a cell or tissue, based on its gene expression and metabolomic profiles. These insights will guide further development of tissue and cell type specific models, and enable us to predict phenotype from genotype. Support or Funding Information The Novo Nordisk Foundation through the Center for Biosustainability at the Technical University of Denmark (NNF10CC1016517 and NNF16CC0021858) and from the NIGMS (R35 GM119850). Tissue specific metabolic models can be constructed using various omics data, but careful decisions on algorithms and parameter choice are necessary to predict a tissue's metabolic functions de novo. This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Analysing Algorithms and Data Sources for the Tissue-Specific Reconstruction of Liver Healthy and Cancer Cells.

Genome-Scale Metabolic Models (GSMMs), mathematical representations of the cell metabolism in different organisms including humans, are resourceful tools to simulate metabolic phenotypes and understand associated diseases, such as obesity, diabetes and cancer. In the last years, different algorithms have been developed to generate tissue-specific metabolic models that simulate different phenotypes for distinct cell types. Hepatocyte cells are one of the main sites of metabolic conversions, mainly due to their diverse physiological functions. Most of the liver's tissue is formed by hepatocytes, being one of the largest and most important organs regarding its biological functions. Hepatocellular carcinoma is, also, one of the most important human cancers with high mortality rates. In this study, we will analyze four different algorithms (MBA, mCADRE, tINIT and FASTCORE) for tissue-specific model reconstruction, based on a template model and two types of data sources: transcriptomics and proteomics. These methods will be applied to the reconstruction of metabolic models for hepatocyte cells and HepG2 cancer cell line. The models will be analyzed and compared under different perspectives, emphasizing their functional analysis considering a set of metabolic liver tasks. The results show that there is no "ideal" algorithm. However, with the current analysis, we were able to retrieve knowledge about the metabolism of the liver.

A Systematic Evaluation of Methods for Tailoring Genome-Scale Metabolic Models

Genome-scale models of metabolism can illuminatethe molecular basis of cell phenotypes. Since some enzymes are only active in specific cell types, several algorithms use omics data to construct cell-line- and tissue-specific metabolic models from genome-scale models. However, these methods are often not rigorously benchmarked, and it is unclear how algorithm and parameter selection (e.g., gene expression thresholds, metabolic constraints) affects model content and predictive accuracy. To investigate this, we built hundreds of models of four different cancer cell lines using six algorithms, four gene expressionthresholds, and three sets of metabolic constraints.Model content varied substantially across different parameter sets, but the algorithms generally increased accuracy in gene essentiality predictions. However, model extraction method choice had the largest impact on model accuracy. We further highlight howassumptions during model development influence model prediction accuracy. These insights will guide further development of context-specific models, thus more accurately resolving genotype-phenotype relationships.

Systematic investigation of metabolic reprogramming in different cancers based on tissue-specific metabolic models.

Cancer cells have different metabolism in contrast to normal cells. The advancement in omics measurement technology enables the genome-wide characterization of altered cellular processes in cancers, but the metabolic flux landscape of cancer is still far from understood. In this study, we compared the well-reconstructed tissue-specific models of five cancers, including breast, liver, lung, renal, and urothelial cancer, and their corresponding normal cells. There are similar patterns in majority of significantly regulated pathways and enriched pathways in correlated reaction sets. But the differences among cancers are also explicit. The renal cancer demonstrates more dramatic difference with other cancer models, including the smallest number of reactions, flux distribution patterns, and specifically correlated pathways. We also validated the predicted essential genes and revealed the Warburg effect by in silico simulation in renal cancer, which are consistent with the measurements for renal cancer. In conclusion, the tissue-specific metabolic model is more suitable to investigate the cancer metabolism. The similarity and heterogenicity of metabolic reprogramming in different cancers are crucial for understanding the aberrant mechanisms of cancer proliferation, which is fundamental for identifying drug targets and biomarkers.

Multi-tissue computational modeling analyzes pathophysiology of type 2 diabetes in MKR mice.

Computational models using metabolic reconstructions for in silico simulation of metabolic disorders such as type 2 diabetes mellitus (T2DM) can provide a better understanding of disease pathophysiology and avoid high experimentation costs. There is a limited amount of computational work, using metabolic reconstructions, performed in this field for the better understanding of T2DM. In this study, a new algorithm for generating tissue-specific metabolic models is presented, along with the resulting multi-confidence level (MCL) multi-tissue model. The effect of T2DM on liver, muscle, and fat in MKR mice was first studied by microarray analysis and subsequently the changes in gene expression of frank T2DM MKR mice versus healthy mice were applied to the multi-tissue model to test the effect. Using the first multi-tissue genome-scale model of all metabolic pathways in T2DM, we found out that branched-chain amino acids' degradation and fatty acids oxidation pathway is downregulated in T2DM MKR mice. Microarray data showed low expression of genes in MKR mice versus healthy mice in the degradation of branched-chain amino acids and fatty-acid oxidation pathways. In addition, the flux balance analysis using the MCL multi-tissue model showed that the degradation pathways of branched-chain amino acid and fatty acid oxidation were significantly downregulated in MKR mice versus healthy mice. Validation of the model was performed using data derived from the literature regarding T2DM. Microarray data was used in conjunction with the model to predict fluxes of various other metabolic pathways in the T2DM mouse model and alterations in a number of pathways were detected. The Type 2 Diabetes MCL multi-tissue model may explain the high level of branched-chain amino acids and free fatty acids in plasma of Type 2 Diabetic subjects from a metabolic fluxes perspective.

Seed architecture shapes embryo metabolism in oilseed rape.

Constrained to develop within the seed, the plant embryo must adapt its shape and size to fit the space available. Here, we demonstrate how this adjustment shapes metabolism of photosynthetic embryo. Noninvasive NMR-based imaging of the developing oilseed rape (Brassica napus) seed illustrates that, following embryo bending, gradients in lipid concentration became established. These were correlated with the local photosynthetic electron transport rate and the accumulation of storage products. Experimentally induced changes in embryo morphology and/or light supply altered these gradients and were accompanied by alterations in both proteome and metabolome. Tissue-specific metabolic models predicted that the outer cotyledon and hypocotyl/radicle generate the bulk of plastidic reductant/ATP via photosynthesis, while the inner cotyledon, being enclosed by the outer cotyledon, is forced to grow essentially heterotrophically. Under field-relevant high-light conditions, major contribution of the ribulose-1,5-bisphosphate carboxylase/oxygenase-bypass to seed storage metabolism is predicted for the outer cotyledon and the hypocotyl/radicle only. Differences between in vitro- versus in planta-grown embryos suggest that metabolic heterogeneity of embryo is not observable by in vitro approaches. We conclude that in vivo metabolic fluxes are locally regulated and connected to seed architecture, driving the embryo toward an efficient use of available light and space.

Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE

BackgroundHuman tissues perform diverse metabolic functions. Mapping out these tissue-specific functions in genome-scale models will advance our understanding of the metabolic basis of various physiological and pathological processes. The global knowledgebase of metabolic functions categorized for the human genome (Human Recon 1) coupled with abundant high-throughput data now makes possible the reconstruction of tissue-specific metabolic models. However, the number of available tissue-specific models remains incomplete compared with the large diversity of human tissues.ResultsWe developed a method called metabolic Context-specificity Assessed by Deterministic Reaction Evaluation (mCADRE). mCADRE is able to infer a tissue-specific network based on gene expression data and metabolic network topology, along with evaluation of functional capabilities during model building. mCADRE produces models with similar or better functionality and achieves dramatic computational speed up over existing methods. Using our method, we reconstructed draft genome-scale metabolic models for 126 human tissue and cell types. Among these, there are models for 26 tumor tissues along with their normal counterparts, and 30 different brain tissues. We performed pathway-level analyses of this large collection of tissue-specific models and identified the eicosanoid metabolic pathway, especially reactions catalyzing the production of leukotrienes from arachidnoic acid, as potential drug targets that selectively affect tumor tissues.ConclusionsThis large collection of 126 genome-scale draft metabolic models provides a useful resource for studying the metabolic basis for a variety of human diseases across many tissues. The functionality of the resulting models and the fast computational speed of the mCADRE algorithm make it a useful tool to build and update tissue-specific metabolic models.