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.

Brain energy metabolism is optimized to minimize the cost of enzyme synthesis and transport.

The energy metabolism of the brain is poorly understood partly due to the complex morphology of neurons and fluctuations in ATP demand over time. To investigate this, we used metabolic models that estimate enzyme usage per pathway, enzyme utilization over time, and enzyme transportation to evaluate how these parameters and processes affect ATP costs for enzyme synthesis and transportation. Our models show that the total enzyme maintenance energy expenditure of the human body depends on how glycolysis and mitochondrial respiration are distributed both across and within cell types in the brain. We suggest that brain metabolism is optimized to minimize the ATP maintenance cost by distributing the different ATP generation pathways in an advantageous way across cell types and potentially also across synapses within the same cell. Our models support this hypothesis by predicting export of lactate from both neurons and astrocytes during peak ATP demand, reproducing results from experimental measurements reported in the literature. Furthermore, our models provide potential explanation for parts of the astrocyte-neuron lactate shuttle theory, which is recapitulated under some conditions in the brain, while contradicting other aspects of the theory. We conclude that enzyme usage per pathway, enzyme utilization over time, and enzyme transportation are important factors for defining the optimal distribution of ATP production pathways, opening a broad avenue to explore in brain metabolism.

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.

Enhanced metabolism and negative regulation of ER stress support higher erythropoietin production in HEK293 cells.

Recombinant protein production can cause severe stress on cellular metabolism, resulting in limited titer and product quality. To investigate cellular and metabolic characteristics associated with these limitations, we compare HEK293 clones producing either erythropoietin (EPO) (secretory) or GFP (non-secretory) protein at different rates. Transcriptomic and functional analyses indicate significantly higher metabolism and oxidative phosphorylation in EPO producers compared with parental and GFP cells. In addition, ribosomal genes exhibit specific expression patterns depending on the recombinant protein and the production rate. In a clone displaying a dramatically increased EPO secretion, we detect higher gene expression related to negative regulation of endoplasmic reticulum (ER) stress, including upregulation of ATF6B, which aids EPO production in a subset of clones by overexpression or small interfering RNA (siRNA) knockdown. Our results offer potential target pathways and genes for further development of the secretory power in mammalian cell factories.

Genome-scale metabolic network reconstruction of model animals as a platform for translational research.

Genome-scale metabolic models (GEMs) are used extensively for analysis of mechanisms underlying human diseases and metabolic malfunctions. However, the lack of comprehensive and high-quality GEMs for model organisms restricts translational utilization of omics data accumulating from the use of various disease models. Here we present a unified platform of GEMs that covers five major model animals, including Mouse1 (Mus musculus), Rat1 (Rattus norvegicus), Zebrafish1 (Danio rerio), Fruitfly1 (Drosophila melanogaster), and Worm1 (Caenorhabditis elegans). These GEMs represent the most comprehensive coverage of the metabolic network by considering both orthology-based pathways and species-specific reactions. All GEMs can be interactively queried via the accompanying web portal Metabolic Atlas. Specifically, through integrative analysis of Mouse1 with RNA-sequencing data from brain tissues of transgenic mice we identified a coordinated up-regulation of lysosomal GM2 ganglioside and peptide degradation pathways which appears to be a signature metabolic alteration in Alzheimer's disease (AD) mouse models with a phenotype of amyloid precursor protein overexpression. This metabolic shift was further validated with proteomics data from transgenic mice and cerebrospinal fluid samples from human patients. The elevated lysosomal enzymes thus hold potential to be used as a biomarker for early diagnosis of AD. Taken together, we foresee that this evolving open-source platform will serve as an important resource to facilitate the development of systems medicines and translational biomedical applications.

Single-cell BCR and transcriptome analysis after influenza infection reveals spatiotemporal dynamics of antigen-specific B cells.

B cell responses are critical for antiviral immunity. However, a comprehensive picture of antigen-specific B cell differentiation, clonal proliferation, and dynamics in different organs after infection is lacking. Here, by combining single-cell RNA and B cell receptor (BCR) sequencing of antigen-specific cells in lymph nodes, spleen, and lungs after influenza infection in mice, we identify several germinal center (GC) B cell subpopulations and organ-specific differences that persist over the course of the response. We discover transcriptional differences between memory cells in lungs and lymphoid organs and organ-restricted clonal expansion. Remarkably, we find significant clonal overlap between GC-derived memory and plasma cells. By combining BCR-mutational analyses with monoclonal antibody (mAb) expression and affinity measurements, we find that memory B cells are highly diverse and can be selected from both low- and high-affinity precursors. By linking antigen recognition with transcriptional programming, clonal proliferation, and differentiation, these finding provide important advances in our understanding of antiviral immunity.

Machine learning-based investigation of the cancer protein secretory pathway.

Deregulation of the protein secretory pathway (PSP) is linked to many hallmarks of cancer, such as promoting tissue invasion and modulating cell-cell signaling. The collection of secreted proteins processed by the PSP, known as the secretome, is often studied due to its potential as a reservoir of tumor biomarkers. However, there has been less focus on the protein components of the secretory machinery itself. We therefore investigated the expression changes in secretory pathway components across many different cancer types. Specifically, we implemented a dual approach involving differential expression analysis and machine learning to identify PSP genes whose expression was associated with key tumor characteristics: mutation of p53, cancer status, and tumor stage. Eight different machine learning algorithms were included in the analysis to enable comparison between methods and to focus on signals that were robust to algorithm type. The machine learning approach was validated by identifying PSP genes known to be regulated by p53, and even outperformed the differential expression analysis approach. Among the different analysis methods and cancer types, the kinesin family members KIF20A and KIF23 were consistently among the top genes associated with malignant transformation or tumor stage. However, unlike most cancer types which exhibited elevated KIF20A expression that remained relatively constant across tumor stages, renal carcinomas displayed a more gradual increase that continued with increasing disease severity. Collectively, our study demonstrates the complementary nature of a combined differential expression and machine learning approach for analyzing gene expression data, and highlights key PSP components relevant to features of tumor pathophysiology that may constitute potential therapeutic targets.

Sources of variation in cell-type RNA-Seq profiles.

Cell-type specific gene expression profiles are needed for many computational methods operating on bulk RNA-Seq samples, such as deconvolution of cell-type fractions and digital cytometry. However, the gene expression profile of a cell type can vary substantially due to both technical factors and biological differences in cell state and surroundings, reducing the efficacy of such methods. Here, we investigated which factors contribute most to this variation. We evaluated different normalization methods, quantified the variance explained by different factors, evaluated the effect on deconvolution of cell type fractions, and examined the differences between UMI-based single-cell RNA-Seq and bulk RNA-Seq. We investigated a collection of publicly available bulk and single-cell RNA-Seq datasets containing B and T cells, and found that the technical variation across laboratories is substantial, even for genes specifically selected for deconvolution, and this variation has a confounding effect on deconvolution. Tissue of origin is also a substantial factor, highlighting the challenge of using cell type profiles derived from blood with mixtures from other tissues. We also show that much of the differences between UMI-based single-cell and bulk RNA-Seq methods can be explained by the number of read duplicates per mRNA molecule in the single-cell sample. Our work shows the importance of either matching or correcting for technical factors when creating cell-type specific gene expression profiles that are to be used together with bulk samples.

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.

UPLC-ESI-MRM/MS for Absolute Quantification and MS/MS Structural Elucidation of Six Specialized Pyranonaphthoquinone Metabolites From Ventilago harmandiana.

Pyranonaphthoquinones (PNQs) are important structural scaffolds found in numerous natural products. Research interest in these specialized metabolites lies in their natural occurrence and therapeutic activities. Nonetheless, research progress has thus far been hindered by the lack of analytical standards and analytical methods for both qualitative and quantitative analysis. We report here that various parts of Ventilago harmandiana are rich sources of PNQs. We developed an ultraperformance liquid chromatography-electrospray ionization multiple reaction monitoring/mass spectrometry method to quantitatively determine six PNQs from leaves, root, bark, wood, and heartwood. The addition of standards in combination with a stable isotope of salicylic acid-D6 was used to overcome the matrix effect with average recovery of 82% ± 1% (n = 15). The highest concentration of the total PNQs was found in the root (11,902 µg/g dry weight), whereas the lowest concentration was found in the leaves (28 µg/g dry weight). Except for the root, PNQ-332 was found to be the major compound in all parts of V. harmandiana, accounting for ∼48% of the total PNQs quantified in this study. However, PNQ-318A was the most abundant PNQ in the root sample, accounting for 27% of the total PNQs. Finally, we provide novel MS/MS spectra of the PNQs at different collision induction energies: 10, 20, and 40 eV (POS and NEG). For structural elucidation purposes, we propose complete MS/MS fragmentation pathways of PNQs using MS/MS spectra at collision energies of 20 and 40 eV. The MS/MS spectra along with our discussion on structural elucidation of these PNQs should be very useful to the natural products community to further exploring PNQs in V. harmandiana and various other sources.

A Systematic Investigation of the Malignant Functions and Diagnostic Potential of the Cancer Secretome.

The collection of proteins secreted from a cell-the secretome-is of particular interest in cancer pathophysiology due to its diagnostic potential and role in tumorigenesis. However, cancer secretome studies are often limited to one tissue or cancer type or focus on biomarker prediction without exploring the associated functions. We therefore conducted a pan-cancer analysis of secretome gene expression changes to identify candidate diagnostic biomarkers and to investigate the underlying biological function of these changes. Using transcriptomic data spanning 32 cancer types and 30 healthy tissues, we quantified the relative diagnostic potential of secretome proteins for each cancer. Furthermore, we offer a potential mechanism by which cancer cells relieve secretory pathway stress by decreasing the expression of tissue-specific genes, thereby facilitating the secretion of proteins promoting invasion and proliferation. These results provide a more systematic understanding of the cancer secretome, facilitating its use in diagnostics and its targeting for therapeutic development.

Targeting CDK2 overcomes melanoma resistance against BRAF and Hsp90 inhibitors.

Novel therapies are undergoing clinical trials, for example, the Hsp90 inhibitor, XL888, in combination with BRAF inhibitors for the treatment of therapy-resistant melanomas. Unfortunately, our data show that this combination elicits a heterogeneous response in a panel of melanoma cell lines including PDX-derived models. We sought to understand the mechanisms underlying the differential responses and suggest a patient stratification strategy. Thermal proteome profiling (TPP) identified the protein targets of XL888 in a pair of sensitive and unresponsive cell lines. Unbiased proteomics and phosphoproteomics analyses identified CDK2 as a driver of resistance to both BRAF and Hsp90 inhibitors and its expression is regulated by the transcription factor MITF upon XL888 treatment. The CDK2 inhibitor, dinaciclib, attenuated resistance to both classes of inhibitors and combinations thereof. Notably, we found that MITF expression correlates with CDK2 upregulation in patients; thus, dinaciclib would warrant consideration for treatment of patients unresponsive to BRAF-MEK and/or Hsp90 inhibitors and/or harboring MITF amplification/overexpression.

An integrated network analysis reveals that nitric oxide reductase prevents metabolic cycling of nitric oxide by Pseudomonas aeruginosa.

Nitric oxide (NO) is a chemical weapon within the arsenal of immune cells, but is also generated endogenously by different bacteria. Pseudomonas aeruginosa are pathogens that contain an NO-generating nitrite (NO2-) reductase (NirS), and NO has been shown to influence their virulence. Interestingly, P. aeruginosa also contain NO dioxygenase (Fhp) and nitrate (NO3-) reductases, which together with NirS provide the potential for NO to be metabolically cycled (NO→NO3-→NO2-→NO). Deeper understanding of NO metabolism in P. aeruginosa will increase knowledge of its pathogenesis, and computational models have proven to be useful tools for the quantitative dissection of NO biochemical networks. Here we developed such a model for P. aeruginosa and confirmed its predictive accuracy with measurements of NO, O2, NO2-, and NO3- in mutant cultures devoid of Fhp or NorCB (NO reductase) activity. Using the model, we assessed whether NO was metabolically cycled in aerobic P. aeruginosa cultures. Calculated fluxes indicated a bottleneck at NO3-, which was relieved upon O2 depletion. As cell growth depleted dissolved O2 levels, NO3- was converted to NO2- at near-stoichiometric levels, whereas NO2- consumption did not coincide with NO or NO3- accumulation. Assimilatory NO2- reductase (NirBD) or NorCB activity could have prevented NO cycling, and experiments with ΔnirB, ΔnirS, and ΔnorC showed that NorCB was responsible for loss of flux from the cycle. Collectively, this work provides a computational tool to analyze NO metabolism in P. aeruginosa, and establishes that P. aeruginosa use NorCB to prevent metabolic cycling of NO.

Integrative analysis of human omics data using biomolecular networks.

High-throughput '-omics' technologies have given rise to an increasing abundance of genome-scale data detailing human biology at the molecular level. Although these datasets have already made substantial contributions to a more comprehensive understanding of human physiology and diseases, their interpretation becomes increasingly cryptic and nontrivial as they continue to expand in size and complexity. Systems biology networks offer a scaffold upon which omics data can be integrated, facilitating the extraction of new and physiologically relevant information from the data. Two of the most prevalent networks that have been used for such integrative analyses of omics data are genome-scale metabolic models (GEMs) and protein-protein interaction (PPI) networks, both of which have demonstrated success among many different omics and sample types. This integrative approach seeks to unite 'top-down' omics data with 'bottom-up' biological networks in a synergistic fashion that draws on the strengths of both strategies. As the volume and resolution of high-throughput omics data continue to grow, integrative network-based analyses are expected to play an increasingly important role in their interpretation.

Starved Escherichia coli preserve reducing power under nitric oxide stress.

Nitric oxide (NO) detoxification enzymes, such as NO dioxygenase (NOD) and NO reductase (NOR), are important to the virulence of numerous bacteria. Pathogens use these defense systems to ward off immune-generated NO, and they do so in environments that contain additional stressors, such as reactive oxygen species, nutrient deprivation, and acid stress. NOD and NOR both use reducing equivalents to metabolically deactivate NO, which suggests that nutrient deprivation could negatively impact their functionality. To explore the relationship between NO detoxification and nutrient deprivation, we examined the ability of Escherichia coli to detoxify NO under different levels of carbon source availability in aerobic cultures. We observed failure of NO detoxification under both carbon source limitation and starvation, and those failures could have arisen from inabilities to synthesize Hmp (NOD of E. coli) and/or supply it with sufficient NADH (preferred electron donor). We found that when limited quantities of carbon source were provided, NO detoxification failed due to insufficient NADH, whereas starvation prevented Hmp synthesis, which enabled cells to maintain their NADH levels. This maintenance of NADH levels under starvation was confirmed to be dependent on the absence of Hmp. Intriguingly, these data show that under NO stress, carbon-starved E. coli are better positioned with regard to reducing power to cope with other stresses than cells that had consumed an exhaustible amount of carbon.

Discovery and dissection of metabolic oscillations in the microaerobic nitric oxide response network of Escherichia coli.

The virulence of many pathogens depends upon their ability to cope with immune-generated nitric oxide (NO·). In Escherichia coli, the major NO· detoxification systems are Hmp, an NO· dioxygenase (NOD), and NorV, an NO· reductase (NOR). It is well established that Hmp is the dominant system under aerobic conditions, whereas NorV dominates anaerobic conditions; however, the quantitative contributions of these systems under the physiologically relevant microaerobic regime remain ill defined. Here, we investigated NO· detoxification in environments ranging from 0 to 50 µM O2, and discovered a regime in which E. coli NO· defenses were severely compromised, as well as conditions that exhibited oscillations in the concentration of NO·. Using an integrated computational and experimental approach, E. coli NO· detoxification was found to be extremely impaired at low O2 due to a combination of its inhibitory effects on NorV, Hmp, and translational activities, whereas oscillations were found to result from a kinetic competition for O2 between Hmp and respiratory cytochromes. Because at least 777 different bacterial species contain the genetic requirements of this stress response oscillator, we hypothesize that such oscillatory behavior could be a widespread phenomenon. In support of this hypothesis,Pseudomonas aeruginosa, whose respiratory and NO· response networks differ considerably from those of E. coli, was found to exhibit analogous oscillations in low O2 environments. This work provides insight into how bacterial NO· defenses function under the low O2 conditions that are likely to be encountered within host environments.

Construction and Experimental Validation of a Quantitative Kinetic Model of Nitric Oxide Stress in Enterohemorrhagic Escherichia coli O157:H7.

Enterohemorrhagic Escherichia coli (EHEC) are responsible for large outbreaks of hemorrhagic colitis, which can progress to life-threatening hemolytic uremic syndrome (HUS) due to the release of Shiga-like toxins (Stx). The presence of a functional nitric oxide (NO·) reductase (NorV), which protects EHEC from NO· produced by immune cells, was previously found to correlate with high HUS incidence, and it was shown that NorV activity enabled prolonged EHEC survival and increased Stx production within macrophages. To enable quantitative study of EHEC NO· defenses and facilitate the development of NO·-potentiating therapeutics, we translated an existing kinetic model of the E. coli K-12 NO· response to an EHEC O157:H7 strain. To do this, we trained uncertain model parameters on measurements of [NO·] and [O2] in EHEC cultures, assessed parametric and prediction uncertainty with the use of a Markov chain Monte Carlo approach, and confirmed the predictive accuracy of the model with experimental data from genetic mutants lacking NorV or Hmp (NO· dioxygenase). Collectively, these results establish a methodology for the translation of quantitative models of NO· stress in model organisms to pathogenic sub-species, which is a critical step toward the application of these models for the study of infectious disease.