Modeling the cancer transcriptome: Recent initiatives such as The Cancer Genome Atlas have mapped the genome-wide effect of individual genes on tumor growth. By unraveling genomic alterations in tumors, molecular subtypes of cancers have been identified, which is improving patient diagnostics and treatment. Uhlen et al. developed a computer-based modeling approach to examine different cancer types in nearly 8000 patients. They provide an open-access resource for exploring how the expression of specific genes influences patient survival in 17 different types of cancer. More than 900,000 patient survival profiles are available, including for tumors of colon, prostate, lung, and breast origin. This interactive data set can also be used to generate personalized patient models to predict how metabolic changes can influence tumor growth. Science , this issue p. eaan2507

Significance: Hepatocellular carcinoma (HCC) is a heterogeneous and deadly form of liver cancer. Here, we stratified and characterized HCC tumors by applying graph and control theory concepts to the topology of genome-scale metabolic networks. We identified three HCC subtypes with distinct differences in metabolic and signaling pathways and clinical survival and validated our results by performing additional experiments. We further identified HCC subtype-specific genes pivotal in controlling the entire metabolism and discovered genes that can be targeted for development of efficient treatment strategies for specific HCC subtypes. Our systems-level analyses provided a systematic way for characterization of HCC subtypes and identification of drug targets for effective treatment of HCC patients.

Viruses hijack host metabolic pathways for their replicative advantage. In this study, using patient-derived multiomics data and in vitro infection assays, we aimed to understand the role of key metabolic pathways that can regulate severe acute respiratory syndrome coronavirus-2 reproduction and their association with disease severity. We used multiomics platforms (targeted and untargeted proteomics and untargeted metabolomics) on patient samples and cell-line models along with immune phenotyping of metabolite transporters in patient blood cells to understand viral-induced metabolic modulations. We also modulated key metabolic pathways that were identified using multiomics data to regulate the viral reproduction in vitro. Coronavirus disease 2019 disease severity was characterized by increased plasma glucose and mannose levels. Immune phenotyping identified altered expression patterns of carbohydrate transporter, glucose transporter 1, in CD8⁺ T cells, intermediate and nonclassical monocytes, and amino acid transporter, xCT, in classical, intermediate, and nonclassical monocytes. In in vitro lung epithelial cell (Calu-3) infection model, we found that glycolysis and glutaminolysis are essential for virus replication, and blocking these metabolic pathways caused significant reduction in virus production. Taken together, we therefore hypothesized that severe acute respiratory syndrome coronavirus-2 utilizes and rewires pathways governing central carbon metabolism leading to the efflux of toxic metabolites and associated with disease severity. Thus, the host metabolic perturbation could be an attractive strategy to limit the viral replication and disease severity.

Genome-scale metabolic models (GEMs) offer insights into cancer metabolism and have been used to identify potential biomarkers and drug targets. Drug repositioning is a time- and cost-effective method of drug discovery that can be applied together with GEMs for effective cancer treatment. In this study, we reconstruct a prostate cancer (PRAD)-specific GEM for exploring prostate cancer metabolism and also repurposing new therapeutic agents that can be used in development of effective cancer treatment. We integrate global gene expression profiling of cell lines with >1000 different drugs through the use of prostate cancer GEM and predict possible drug-gene interactions. We identify the key reactions with altered fluxes based on the gene expression changes and predict the potential drug effect in prostate cancer treatment. We find that sulfamethoxypyridazine, azlocillin, hydroflumethiazide, and ifenprodil can be repurposed for the treatment of prostate cancer based on an in silico cell viability assay. Finally, we validate the effect of ifenprodil using an in vitro cell assay and show its inhibitory effect on a prostate cancer cell line. Our approach demonstate how GEMs can be used to predict therapeutic agents for cancer treatment based on drug repositioning. Besides, it paved a way and shed a light on the applicability of computational models to real-world biomedical or pharmaceutical problems.

Dysbiosis in the gut microbiome composition may be promoted by therapeutic drugs such as metformin, the world’s most prescribed antidiabetic drug. Under metformin treatment, disturbances of the intestinal microbes lead to increased abundance of Escherichia spp., Akkermansia muciniphila, Subdoligranulum variabile and decreased abundance of Intestinibacter bartlettii. This alteration may potentially lead to adverse effects on the host metabolism, with the depletion of butyrate producer genus. However, an increased production of butyrate and propionate was verified in metformin-treated Type 2 diabetes (T2D) patients. The mechanisms underlying these nutritional alterations and their relation with gut microbiota dysbiosis remain unclear. Here, we used Genome-scale Metabolic Models of the representative gut bacteria Escherichia spp., I. bartlettii, A. muciniphila, and S. variabile to elucidate their bacterial metabolism and its effect on intestinal nutrient pool, including macronutrients (e.g., amino acids and short chain fatty acids), minerals and chemical elements (e.g., iron and oxygen). We applied flux balance analysis (FBA) coupled with synthetic lethality analysis interactions to identify combinations of reactions and extracellular nutrients whose absence prevents growth. Our analyses suggest that Escherichia sp. is the bacteria least vulnerable to nutrient availability. We have also examined bacterial contribution to extracellular nutrients including short chain fatty acids, amino acids, and gasses. For instance, Escherichia sp. and S. variabile may contribute to the production of important short chain fatty acids (e.g., acetate and butyrate, respectively) involved in the host physiology under aerobic and anaerobic conditions. We have also identified pathway susceptibility to nutrient availability and reaction changes among the four bacteria using both FBA and flux variability analysis. For instance, lipopolysaccharide synthesis, nucleotide sugar metabolism, and amino acid metabolism are pathways susceptible to changes in Escherichia sp. and A. muciniphila. Our observations highlight important commensal and competing behavior, and their association with cellular metabolism for prevalent gut microbes. The results of our analysis have potential important implications for development of new therapeutic approaches in T2D patients through the development of prebiotics, probiotics, or postbiotics.

Redox metabolism is often considered a potential target for cancer treatment, but a systematic examination of redox responses in hepatocellular carcinoma (HCC) is missing. Here, we employed systems biology and biological network analyses to reveal key roles of genes associated with redox metabolism in HCC by integrating multi-omics data. We found that several redox genes, including 25 novel potential prognostic genes, are significantly co-expressed with liver-specific genes and genes associated with immunity and inflammation. Based on an integrative analysis, we found that HCC tumors display antagonistic behaviors in redox responses. The two HCC groups are associated with altered fatty acid, amino acid, drug and hormone metabolism, differentiation, proliferation, and NADPH-independent vs -dependent antioxidant defenses. Redox behavior varies with known tumor subtypes and progression, affecting patient survival. These antagonistic responses are also displayed at the protein and metabolite level and were validated in several independent cohorts. We finally showed the differential redox behavior using mice transcriptomics in HCC and noncancerous tissues and associated with hypoxic features of the two redox gene groups. Our integrative approaches highlighted mechanistic differences among tumors and allowed the identification of a survival signature and several potential therapeutic targets for the treatment of HCC.