Transcriptional reprogramming of cellular metabolism is a hallmark of cancer. However, systematic approaches to study the role of transcriptional regulators (TRs) in mediating cancer metabolic rewiring are missing. Here, we chart a genome-scale map of TR-metabolite associations in human cells using a combined computational-experimental framework for large-scale metabolic profiling of adherent cell lines. By integrating intracellular metabolic profiles of 54 cancer cell lines with transcriptomic and proteomic data, we unraveled a large space of associations between TRs and metabolic pathways. We found a global regulatory signature coordinating glucose- and one-carbon metabolism, suggesting that regulation of carbon metabolism in cancer may be more diverse and flexible than previously appreciated. Here, we demonstrate how this TR-metabolite map can serve as a resource to predict TRs potentially responsible for metabolic transformation in patient-derived tumor samples, opening new opportunities in understanding disease etiology, selecting therapeutic treatments and in designing modulators of cancer-related TRs.

The production of recombinant proteins is frequently enhanced at the levels of transcription, codon usage, protein folding and secretion. Overproduction of heterologous proteins, however, also directly affects the primary metabolism of the producing cells. By incorporation of the production of a heterologous protein into a genome scale metabolic model of the yeast Pichia pastoris, the effects of overproduction were simulated and gene targets for deletion or overexpression for enhanced productivity were predicted. Overexpression targets were localized in the pentose phosphate pathway and the TCA cycle, while knockout targets were found in several branch points of glycolysis. Five out of 9 tested targets led to an enhanced production of cytosolic human superoxide dismutase (hSOD). Expression of bacterial β-glucuronidase could be enhanced as well by most of the same genetic modifications. Beneficial mutations were mainly related to reduction of the NADP/H pool and the deletion of fermentative pathways. Overexpression of the hSOD gene itself had a strong impact on intracellular fluxes, most of which changed in the same direction as predicted by the model. In vivo fluxes changed in the same direction as predicted to improve hSOD production. Genome scale metabolic modeling is shown to predict overexpression and deletion mutants which enhance recombinant protein production with high accuracy.

No abstract

Metabolic profiling of cell line collections has become an invaluable tool to study disease etiology, drug modes of action and to select personalized treatments. However, large-scale in vitro dynamic metabolic profiling is limited by time-consuming sampling and complex measurement procedures. By adapting a mass spectrometry-based metabolomics workflow for high-throughput profiling of diverse adherent mammalian cells, we establish a framework for the rapid measurement and analysis of drug-induced dynamic changes in intracellular metabolites. This methodology is scalable to large compound libraries and is here applied to study the mechanism underlying the toxic effect of dichloroacetate in ovarian cancer cell lines. System-level analysis of the metabolic responses revealed a key and unexpected role of CoA biosynthesis in dichloroacetate toxicity and the more general importance of CoA homeostasis across diverse human cell lines. The herein-proposed strategy for high-content drug metabolic profiling is complementary to other molecular profiling techniques, opening new scientific and drug-discovery opportunities.

Metabolic flux analysis is based on the measurement of isotopologue ratios. In this work, a new GC-MS-based method was introduced enabling accurate determination of isotopologue distributions of sugar phosphates in cell extracts. A GC-TOFMS procedure was developed involving a two-step online derivatization (ethoximation followed by trimethylsilylation) offering high mass resolution, high mass accuracy and the potential of retrospective data analysis typical for TOFMS. The information loss due to fragmentation intrinsic for isotopologue analysis by electron ionization could be overcome by chemical ionization with methane. A thorough optimization regarding pressure of the reaction gas, emission current, electron energy and temperature of the ion source was carried out. For a substantial panel of sugar phosphates both of the glycolysis and the pentose phosphate pathway, sensitive determination of the protonated intact molecular ions together with low abundance fragment ions was successfully achieved. The developed method was evaluated for analysis of Pichia pastoris cell extracts. The measured isotopologue ratios were in the range of 55:1-2:1. The comparison of the experimental isotopologue fractions with the theoretical fractions was excellent, revealing a maximum bias of 4.6% and an average bias of 1.4%.

The sulfur metabolic pathway is involved in basic modes of cellular metabolism, including methylation, cell division, respiratory oscillations and stress responses. Hence, the implicated high reactivity of the sulfur pathway intermediates entails challenges for their quantitative analysis. In particular the unwanted oxidation of the thiol group-containing metabolites glutathione, cysteine, homocysteine, γ-glutamyl cysteine and cysteinyl glycine must be prevented in order to obtain accurate snapshots of this important part of cellular metabolism. Suitable analytical methodologies are therefore needed to support studies of drug metabolism and metabolic engineering. In this work, a novel sample preparation strategy targeting thiolic metabolites was established by implementing thiol group protection with N-ethyl maleimide using a cold methanol metabolite extraction procedure. It was shown that N-ethyl maleimide derivatization is compatible with typical metabolite extraction procedures and also allowed for the stabilization of the instable thiolic metabolites in a fully (13)C-labeled yeast cell extract. The stable isotope labeled metabolite analogs could be used for internal standardization to achieve metabolite quantification with high precision. Furthermore, a dedicated hydrophilic interaction liquid chromatography tandem mass spectrometry method for the separation of sulfur metabolic pathway intermediates using a sub-2 μm particle size stationary phase was developed. Coupled with tandem mass spectrometry, the presented methodology proved to be robust, and sensitive (absolute detection limits in the low femtomole range), and allowed for the quantification of cysteine, cysteinyl glycine, cystathionine, cystine, glutamic acid, glutamyl cysteine, reduced glutathione, glutathione disulfide, homocysteine, methionine, S-adenosyl homocysteine and serine in a human ovarian carcinoma cell model.

A novel on-line combination of reversed phase and porous graphitized carbon liquid chromatography increases the versatility in non-targeted metabolomics.

Molecular profiling of small-molecules offers invaluable insights into the function of compounds and allows for hypothesis generation about small molecule direct targets and secondary effects. However, current profiling methods are either limited in the number of measurable parameters or throughput. Here, we developed a multiplexed, unbiased framework that, by linking genetic to drug-induced changes in nearly a thousand metabolites, allows for high-throughput functional annotation of compound libraries in Escherichia coli. First, we generated a reference map of metabolic changes from (CRISPR) interference with 352 genes in all major essential biological processes. Next, based on the comparison of genetic with 1342 drug-induced metabolic changes we made de novo predictions of compound functionality and revealed antibacterials with unconventional Modes of Action. We show that our framework, combining dynamic gene silencing with metabolomics, can be adapted as a general strategy for comprehensive high-throughput analysis of compound functionality, from bacteria to human cell lines.