Generation and analysis of context-specific genome-scale metabolic models derived from single-cell RNA-Seq data.

Single-cell RNA sequencing combined with genome-scale metabolic models (GEMs) has the potential to unravel the differences in metabolism across both cell types and cell states but requires new computational methods. Here, we present a method for generating cell-type-specific genome-scale models from clusters of single-cell RNA-Seq profiles. Specifically, we developed a method to estimate the minimum number of cells required to pool to obtain stable models, a bootstrapping strategy for estimating statistical inference, and a faster version of the task-driven integrative network inference for tissues algorithm for generating context-specific GEMs. In addition, we evaluated the effect of different RNA-Seq normalization methods on model topology and differences in models generated from single-cell and bulk RNA-Seq data. We applied our methods on data from mouse cortex neurons and cells from the tumor microenvironment of lung cancer and in both cases found that almost every cell subtype had a unique metabolic profile. In addition, our approach was able to detect cancer-associated metabolic differences between cancer cells and healthy cells, showcasing its utility. We also contextualized models from 202 single-cell clusters across 19 human organs using data from Human Protein Atlas and made these available in the web portal Metabolic Atlas, thereby providing a valuable resource to the scientific community. With the ever-increasing availability of single-cell RNA-Seq datasets and continuously improved GEMs, their combination holds promise to become an important approach in the study of human metabolism.

Metastatic triple negative breast cancer adapts its metabolism to destination tissues while retaining key metabolic signatures.

Triple negative breast cancer (TNBC) metastases are assumed to exhibit similar functions in different organs as in the original primary tumor. However, studies of metastasis are often limited to a comparison of metastatic tumors with primary tumors of their origin, and little is known about the adaptation to the local environment of the metastatic sites. We therefore used transcriptomic data and metabolic network analyses to investigate whether metastatic tumors adapt their metabolism to the metastatic site and found that metastatic tumors adopt a metabolic signature with some similarity to primary tumors of their destinations. The extent of adaptation, however, varies across different organs, and metastatic tumors retain metabolic signatures associated with TNBC. Our findings suggest that a combination of anti-metastatic approaches and metabolic inhibitors selected specifically for different metastatic sites, rather than solely targeting TNBC primary tumors, may constitute a more effective treatment approach.

Metabolic collaboration between cells in the tumor microenvironment has a negligible effect on tumor growth.

The tumor microenvironment is composed of a complex mixture of different cell types interacting under conditions of nutrient deprivation, but the metabolism therein is not fully understood due to difficulties in measuring metabolic fluxes and exchange of metabolites between different cell types in vivo. Genome-scale metabolic modeling enables estimation of such exchange fluxes as well as an opportunity to gain insight into the metabolic behavior of individual cell types. Here, we estimated the availability of nutrients and oxygen within the tumor microenvironment using concentration measurements from blood together with a metabolite diffusion model. In addition, we developed an approach to efficiently apply enzyme usage constraints in a comprehensive metabolic model of human cells. The combined modeling reproduced severe hypoxic conditions and the Warburg effect, and we found that limitations in enzymatic capacity contribute to cancer cells' preferential use of glutamine as a substrate to the citric acid cycle. Furthermore, we investigated the common hypothesis that some stromal cells are exploited by cancer cells to produce metabolites useful for the cancer cells. We identified over 200 potential metabolites that could support collaboration between cancer cells and cancer-associated fibroblasts, but when limiting to metabolites previously identified to participate in such collaboration, no growth advantage was observed. Our work highlights the importance of enzymatic capacity limitations for cell behaviors and exemplifies the utility of enzyme-constrained models for accurate prediction of metabolism in cells and tumor microenvironments.

Ginsenoside Rg3 Reduces the Toxicity of Graphene Oxide Used for pH-Responsive Delivery of Doxorubicin to Liver and Breast Cancer Cells.

Doxorubicin (DOX) is extensively used in chemotherapy, but it has serious side effects and is inefficient against some cancers, e.g., hepatocarcinoma. To ameliorate the delivery of DOX and reduce its side effects, we designed a pH-responsive delivery system based on graphene oxide (GO) that is capable of a targeted drug release in the acidic tumor microenvironment. GO itself disrupted glutathione biosynthesis and induced reactive oxygen species (ROS) accumulation in human cells. It induced IL17-directed JAK-STAT signaling and VEGF gene expression, leading to increased cell proliferation as an unwanted effect. To counter this, GO was conjugated with the antioxidant, ginsenoside Rg3, prior to loading with DOX. The conjugation of Rg3 to GO significantly reduced the toxicity of the GO carrier by abolishing ROS production. Furthermore, treatment of cells with GO-Rg3 did not induce IL17-directed JAK-STAT signaling and VEGF gene expression-nor cell proliferation-suggesting GO-Rg3 as a promising drug carrier. The anticancer activity of GO-Rg3-DOX conjugates was investigated against Huh7 hepatocarcinoma and MDA-MB-231 breast cancer cells. GO-Rg3-DOX conjugates significantly reduced cancer cell viability, primarily via downregulation of transcription regulatory genes and upregulation of apoptosis genes. GO-Rg3 is an effective, biocompatible, and pH responsive DOX carrier with potential to improve chemotherapy-at least against liver and breast cancers.

Augmentation of oxidative stress-induced apoptosis in MCF7 cells by ascorbate-tamoxifen and/or ascorbate-juglone treatments.

Since reactive oxygen species (ROS) play diverse roles in cancer, modulating the redox status of cancerous cells seems to be a promising therapeutic approach. Oxidant-targeted therapy appears logical for intervention with the acquired adaptive response to oxidative stress in cancer. In this study, we investigated the cytotoxic effects of juglone (J) and tamoxifen (T) and also the combination of each with ascorbate (A): tamoxifen/ascorbate (TA) and/or juglone/ascorbate (JA) on MCF7 cancerous cells. The results revealed that the growth inhibitory effects of juglone and tamoxifen were each associated with enhanced levels of ROS production and lipid peroxidation. These effects were markedly intensified in tamoxifen/ascorbate and juglone/ascorbate co-treatments. On the other hand, the intracellular anti-oxidant components such as reduced glutathione (GSH), catalase, superoxide dismutase (SOD), and glutathione peroxidase significantly declined in cells subjected to combination treatments compared to that in cells exposed solely to tamoxifen, juglone, and the untreated control cells. In addition, ascorbate association induced more apoptotic and necrotic or necrotic-like cell death than cells treated with each drug alone. These results were further confirmed by comparing the Bax/Bcl2 ratio in combination-treated cells. Additionally, ascorbate was able to potentiate the cytotoxic effects of combination therapy via activation of ROS-responsive factors including Foxo family members.