Next Article in Journal / Special Issue
Quantitative Microdialysis: Experimental Protocol and Software for Small Molecule Protein Affinity Determination and for Exclusion of Compounds with Poor Physicochemical Properties
Previous Article in Journal
Unexpectedly High Levels of Inverted Re-Insertions Using Paired sgRNAs for Genomic Deletions
Previous Article in Special Issue
Synthesis of Quenchbodies for One-Pot Detection of Stimulant Drug Methamphetamine
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
MPs
Volume 3
Issue 3
10.3390/mps3030054
mps-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

IsoSearch: An Untargeted and Unbiased Metabolite and Lipid Isotopomer Tracing Strategy from HR-LC-MS/MS Datasets

by
He Huang
1,2,3,
Min Yuan
1,
Phillip Seitzer
4,5,
Susan Ludwigsen
4 and
John M. Asara
1,2,*
1
Division of Signal Transduction, Beth Israel Deaconess Medical Center, Boston, MA 02115, USA
2
Department of Medicine, Harvard Medical School, Boston, MA 02115, USA
3
Institute of Metabolism and Integrative Biology, Fudan University, Shanghai 200438, China
4
Proteome Software, Inc., Portland, OR 97219, USA
5
Calico Life Sciences, San Francisco, CA 94080, USA
*
Author to whom correspondence should be addressed.
Methods Protoc. 2020, 3(3), 54; https://0-doi-org.brum.beds.ac.uk/10.3390/mps3030054
Submission received: 25 June 2020 / Revised: 27 July 2020 / Accepted: 28 July 2020 / Published: 30 July 2020
(This article belongs to the Special Issue Feature Issue on the Development of Novel Cutting Edge Protocols)
Browse Figures
Versions Notes

Abstract

:
Stable isotopic tracer analysis is a technique used to determine carbon or nitrogen atom incorporation into biological systems. A number of mass spectrometry based approaches have been developed for this purpose, including high-resolution tandem mass spectrometry (HR-LC-MS/MS), selected reaction monitoring (SRM) and parallel reaction monitoring (PRM). We have developed an approach for analyzing untargeted metabolomic and lipidomic datasets using high-resolution mass spectrometry with polarity switching and implemented our approach in the open-source R script IsoSearch and in Scaffold Elements software. Using our strategy, which requires an unlabeled reference dataset and isotope labeled datasets across various biological conditions, we traced metabolic isotopomer alterations in breast cancer cells (MCF-7) treated with the metabolic drugs 2-deoxy-glucose, 6-aminonicotinamide, compound 968, and rapamycin. Metabolites and lipids were first identified by the commercial software Scaffold Elements and LipidSearch, then IsoSearch successfully profiled the 13C-isotopomers extracted metabolites and lipids from 13C-glucose labeled MCF-7 cells. The results interpreted known models, such as glycolysis and pentose phosphate pathway inhibition, but also helped to discover new metabolic/lipid flux patterns, including a reactive oxygen species (ROS) defense mechanism induced by 6AN and triglyceride accumulation in rapamycin treated cells. The results suggest the IsoSearch/Scaffold Elements platform is effective for studying metabolic tracer analysis in diseases, drug metabolism, and metabolic engineering for both polar metabolites and non-polar lipids.

1. Introduction

Metabolic flux analysis (MFA) using stable isotope tracers is a technique used to investigate the intracellular metabolic rate of cells and organisms [1,2]. Using nuclear magnetic resonance (NMR) or mass spectrometry (MS), steady-state metabolite flux in biological systems has been widely studied in the past few decades [3,4,5]. Although unlabeled steady-state metabolomics profiling can provide information on metabolite level alterations, the directionality associated with metabolic fluxes is not easily obtained [6]. Stable isotope labeled nutrients, such as 13C-glucose and 15N-glutamine, are frequently used in targeted mass spectrometry (GC-MS or LC-MS) based approaches to analyze metabolite isotopomer incorporations [1,6,7]. In particular, selected reaction monitoring (SRM) is often used to detect the heavy isotope labeled metabolites identified by their unique precursor ion and fragment ion transitions based on the number of labeled carbons. The SRM based targeted flux analysis built in our lab has been used to profile central carbon metabolism in various diseases. In references [8,9,10,11,12,13] However, like most existing MFA approaches, targeted fluxomics are limited to selected metabolite reactions in defined metabolic pathways (such as central carbon metabolism, amino acid metabolism, nucleotides, etc.); thus, a global metabolic fluxomics method is needed to discover new metabolic reactions from alternative pathways including minor metabolic contributions from reaction intermediates. Many laboratories have developed various software tools to interpret the MFA results generated by NMR [14], GC-MS [15,16,17], or LC-MS, such as SUMOFLUX [18], 13C-FLUX [14,19], X13CMS [20], and Omix [21]. SUMOFLUX is a MATLAB based program that uses surrogate modeling and machine learning techniques to analyze targeted 13C-metabolic flux models. 13C-FLUX is a toolbox that takes GC-MS data or converted NMR data for targeted metabolic flux analysis [15]. The 13C-FLUX results can also be visualized by the Omix software. Although both SUMOFLUX and 13C-FLUX can process targeted metabolic flux results, they are less powerful for managing untargeted fluxomics results and discovering new metabolic flux mechanisms. X13CMS aims to analyze untargeted metabolic isotopic data by requiring a paired isotope labeled experiment with an unlabeled sample, and differentiates the isotopomers based on the mass differences. Using X13CMS, several new carbon source fluxes were discovered, such as itaconate flux in the polarized macrophages [22], acetoacetate shuttles between the hepatocyte and macrophage [23], and mammalian mitochondrial lactate flux [24]. However, it is not always possible to obtain the paired isotopic experiment samples required by X13CMS due to acquisition conditions. Moreover, the dependence on XCMS outputs makes X13CMS unlikely to incorporate results from other software such as Proteome Software’s Scaffold Elements [25] and Thermo’s LipidSearch.
Here, we introduce an untargeted metabolic isotopomer tracing strategy which uses high resolution LC-MS/MS in polarity switching mode to trace the stable isotopic tracer’s (i.e., 13C[6]-glucose) fate and profile the results using an in-house R program named IsoSearch. As the core portion of the untargeted isotopomer tracing workflow (Figure 1), IsoSearch can take untargeted metabolomics or lipidomics results generated by high-resolution MS and deconvolute the heavy isotopic metabolite or lipid spectra, and eventually used for flux analysis. Identities of metabolites and lipids are obtained via data-dependent acquisition (DDA) and Scaffold Elements software (Proteome Software) matching against the (National Institute of Standards and Technology (NIST) 2017 and Human Metabolome Database (HMDB) spectral libraries, or by using LipidSearch [26] from the unlabeled dataset. Two sample sets (unlabeled and 13C or 15N labeled) are acquired by the untargeted metabolomic or lipidomic platform with a positive/negative polarity switch high-resolution MS. The isotopomers of the metabolites and lipids are interpreted from 13C isotopic labeled samples using the IsoSearch program. Using IsoSearch, every identified metabolite and lipid molecule can produce an isotopomer pattern since IsoSearch pools unlabeled samples and compares each [13] C isotopic labeled sample to the unlabeled pool. One of the core features of the IsoSearch algorithm is to screen out potential false positives using a score that incorporates both m/z and retention time (RT) simultaneously. We also have introduced this method as a feature in the commercial Scaffold Elements 2.1 software (http://www.proteomesoftware.com/products/elements/) which provides a user-friendly interface. As an example of the technology, we show MCF-7 human breast cancer cells after treatment with known inhibitors of central carbon metabolism and compare how universally labeled 13C[6]-glucose carbon atoms are tracked throughout metabolic processes. Our specific focus was concentrated on aerobic and anaerobic glycolysis, pentose phosphate pathway (PPP), tricarboxylic acid cycle (TCA), and glucose associated phosphatidylinositol 3-kinase/Akt kinase/mammalian target of rapamycin (PI3K/Akt/mTOR) signaling pathways to examine the well-known glucose involved central carbon metabolism, but also to reveal new metabolic fluxes in alternative and potentially unknown pathways.

2. Results

2.1. Development of IsoSearch Strategy-Identification and Peak Picking

Two sets of RAW untargeted high-resolution data were acquired, unlabeled (12C) data from control samples and 13C labeled datasets from various experimental cell conditions. The 12C-unlabeled control files were processed by Scaffold Elements 2.1 software. Using a 0.7 ID score and 0.5 MS2 score cut-off against NIST 2017 database (652,475 MS/MS spectra), more than 470 metabolites were selected. The identified metabolites were then exported as a pooled reference peak list with all related information including the m/z and retention time (RT). The 12C-unlabeled lipid. RAW files were processed similarly to obtain the pooled reference peak list with lipid m/z and RT information. However, instead of using Scaffold Elements, lipids were identified using Thermo LipidSearch 4.1.3 software. Lipids graded A-C and m-score ≥5.0 were accepted. MS2 data is necessary only for the identification process. Further isotope tracer analysis is performed on the MS1 level.
The 13C-labeled RAW files from both metabolite and lipid datasets were converted to mzXML files and MS1 peak lists were created using the R package enviPick(v1.0) (https://rdrr.io/cran/enviPick/). The parameters of enviPick() functions to generate the peak lists are listed in Table 1 and enviPick parameters used for our high resolution data from the QEaxctive Plus/HF Orbitrap was as follows: enviPickwrap (filepath.mzXML, MSlevel = c(1), dmzgap = 15, dmzdens = 4, ppm = TRUE, drtgap = 300, drtsmall = 20, drtdens = 250, drtfill = 10, drttotal = 200, minpeak = 3, recurs = 3, weight = 3, SB = 2, SN = 3, minint = 1E4, maxint = 1E7, ended = 2, ion_mode = positive/negative, progbar = FALSE).

2.2. IsoSearch Analysis

The 13C-labeled peak lists were then compared against the 12C-unlabeled pooled reference peak list using an in-house R program “IsoSearch”. It is important to note that only well annotated identified lipids or metabolites from the unlabeled reference samples can be used for further isotope tracer studies. The IsoSearch program performs the untargeted metabolomic/lipidomic isotope tracing analysis in two steps using multiple functions (Table 1). In the first step, IsoSearch reconstructs the 12C-unlabeled reference peak list by adding isotopic mass difference (13C/z = 1.0033548378/z) to each metabolite/lipid feature’s m/z value. Then, the 13C-labeled samples are searched against the reconstructed 12C-unlabeled reference peak list to sort out the metabolite isotopomers that have the closest m/z and RT. The IsoSearch program introduced a scoring algorithm to adjust the metabolite’s isotopomer accuracy and screen out the false positive identifications using the equation:
S c o r e = 1 − ( Δ   m / z × β ) 2 +   ( Δ   R T ) 2
where in the Δ m/z and Δ RT are the m/z and RT differences between isotopic labeled feature list vs. reference list, respectively. β (s ∙ C/g) is a unit-fixed parameter, which can justify the m/z and RT weight coefficient and normalize unit differences. In this article, we set β value to 0.1 for both the metabolite and lipid isotopomer searches. For high resolution and high mass accuracy data with Orbitrap technology, the Δ m/z is~1.4 ppm and Δ RT is~0.14 sec. For polar metabolites, the maximum number of possible isotopomers is 10, while the number of non-polar lipid isotopomers maxes out at 40 because of the relative number of carbon atoms in each molecule type.
The results generated by the IsoSearch:Flux_result function were summarized as a spreadsheet (.csv) whose attributes were listed in the Table 2.
The above rules and parameters for IsoSearch have been adopted and incorporated into Scaffold Elements metabolomics software v2.1 distributed by Proteome Software, Inc. (Portland, OR, USA), (http://www.proteomesoftware.com/products/elements/).

2.3. LC-MS/MS Data Preprocessing

The SRM data generated from targeted metabolomics were integrated to peak area values using MultiQuant 3.0 software (AB/SCIEX) and text output. The targeted data was acquired only to test the competency of the untargeted metabolomics IsoSearch workflow. More experimental details, including the full isotopomer datasets from both lipids and metabolites, are included in the Supplementary Materials Section.
Since we have a well established method for targeted metabolic flux analysis in our laboratory [8], we developed a methodology for unbiased and untargeted metabolic isotopomer tracing analysis suitable for any biological source whereby 13C or 15N atoms can be incorporated into the organism’s metabolism. Our goal was to utilize high resolution untargeted metabolomics and lipidomics as a way to trace the isotopomers and to discover more metabolite features. We designed the technique and software parameters using MCF-7 breast cancer cells. In order to test the power of the flux methodology for untargeted metabolomics and lipidomics data, we validated known metabolites in central carbon metabolism against targeted SRM technology [8].

3. Discussion

3.1. Principle of IsoSearch Based Metabolic/Lipid Isotopomer Tracing

Two major components of untargeted 13C-isotopomer tracing should be considered in the workflow; (1) discovering unbiased metabolites for flux analysis; and (2) minimizing false positive isotopomer IDs. The IsoSearch software implements untargeted heavy isotopic pattern interpretation with low false positive rate. To demonstrate the capabilities of this untargeted isotopomer-tracing platform, we used four drugs that inhibit metabolic or kinase signaling pathways in the MCF-7 breast cancer cell line, which was cultured in 13C[6]-glucose media, and LC-MS/MS profiled the 13C isotopomers at three time points after drug treatments (2 h, 16 h, and 24 h). Metabolites and lipids were extracted and run via both untargeted metabolomic and lipidomics with positive/negative polarity switching. The analysis of all metabolites and lipids as well as their associated 13C isotopomers was unbiased. The metabolic/lipid heavy isotopomers are the features that elute at the same time as the unlabeled metabolites with m/z values that increase by M + 1, M + 2, M + 3, etc. IsoSearch selects all the 13C isotopomers that have the identical retention time (RT) and theoretical m/z, and also filters off the unrelated features (i.e., the false positives). Due to the high degree of carbon incorporation of lipid isotopomers, falsely identified features occur more frequently in the lipid isotopomer profiling. It should be noted that IsoSearch is capable of accurately picking out the correct lipid isotopomers and screening out all of the false positives.
Another mission of untargeted 13C-fluxomics using IsoSearch technology is to identify novel metabolic and lipid features and discover 13C incorporation of metabolites and lipids detected in both positive and negative modes. The metabolites or lipids in 12C unlabeled samples were identified using Scaffold Elements (NIST 2017 MS/MS spectral library) or LipidSearch (internal lipid library), respectively. The identified metabolites or lipids from the unlabeled samples were then pooled as an identification list. The IsoSearch algorithm uses the pooled identification list as the reference for identifying and quantifying 13C isotopomers, and 13C-labeled samples are compared to the pooled metabolite/lipid reference list with a tight m/z and RT match. The identification of metabolite/lipid 13C isotopomer distribution is dependent upon the metabolite/lipid candidates that are identified via DDA from the reference samples; thus, the more molecules that are discovered in the reference samples, the more metabolites and their associated isotopomer changes can be revealed in the 13C-labeled samples by the untargeted approach. Since multiple LC-MS/MS runs typically result in additional small molecule identifications and increased metabolite coverage, we performed multiple unlabeled runs over various treatments to expand our identification reference pool for subsequent 13C isotopic pattern discovering [27,28]. For untargeted metabolic fluxomics, Scaffold Elements identified a total of 340 unique metabolites, and IsoSearch found that 241 of those metabolites were paired with at least one 13C isotopomer for a total of 969 associated 13C isotopomers (Figure S1A). In addition, LipidSearch identified 844 total unique lipid ions and 485 of them contained a total of 2369 13C isotopomers (Figure S1B).

3.2. Untargeted Isotopic Patterns Are Consistent with the Targeted Isotopomer Detection

We set out to validate the results of IsoSearch software by comparing the isotopomer results of several overlapping metabolites with a well-established targeted SRM 13C isotopic analysis technology [8,9,10]. As an example, 12C-unlabeled and 13C-labeled uridine triphosphate (UTP) was detected using both the targeted and untargeted methods (Figure 2A). The 13C[6]-glucose labeled cells can flux 13C atoms through UTP, which accumulates in the (M+5) isotopomer peak in UTP, likely in the pentose sugar [29]. The high-resolution MS spectra detected in the untargeted metabolomic analysis interpreted by IsoSearch also indicate the same isotopic pattern as the targeted method. Several other metabolites including alanine and lactate are also detected the same isotopic patterns in both untargeted and targeted approaches (Figure S2), thus, further validating the IsoSearch algorithm. In addition, the same metabolite (S-adenosylmethionine) isotopomer profile can be obtained with both IsoSearch and Scaffold Elements 2.1 (Figure S3), which contains a 13C metabolic flux algorithm built from the parameters used by IsoSearch.
In addition to IsoSearch’s interpretation of the metabolite isotopomer patterns, intact lipid isotopomers can also be profiled. Most phospholipids contain a large number of carbon atoms (>30) which makes isotopomer interpretation difficult since lipids of similar size of the same class are more likely to have a similar retention time. Like metabolites, lipids are first identified using the M+0 (12C only) ion from a MS/MS based identification software (Thermo LipidSearch). Since lipids are larger in size than metabolites with many more carbon atoms, peak overlap can often occur. The phosphatidylcholine lipid, (PC(16:0/16:0)) (m/z 778.56) was detected in both 12C-unlabeled and 13C-labeled samples (Figure 2B), and co-eluted with a heavier lipid at m/z 804.58 [PC(16:0/18:1) + (HCOO-)]. Although the raw spectra indicate that the (M+30) to (M+40) peaks of PC(16:0/16:0) are overlapping with the (M+0) to (M+10) peaks of the heavier lipid (Figure 2B, left panel), using the score screening algorithm d, IsoSearch can differentiate the overlapping peaks from PC(16:0/16:0) isotopomers (Figure 2B, right panel).

3.3. Untargeted Metabolic Time-Course Analysis Reveals Drug Stimulated 13C[6]-Glucose Flux Alteration

Cancer signaling pathways are frequently associated with metabolic regulations [30,31,32]. In order to demonstrate the IsoSearch strategy and study metabolic regulation in cancer cells, the estrogen receptor (ER) positive MCF-7 breast cancer cell line, was cultured in either 12 C regular glucose DMEM or with the addition of universally 13C labeled (U-13C6) glucose DMEM to investigate the glucose metabolic flux. Four metabolism or kinase inhibitors, 2-deoxy-glucose (2DG), 6-aminonicotinamide (6AN), compound 968 (C968) and rapamycin (Rapa), were added to the cell culture media to alter the flux of glucose in MCF-7 for different purposes.
To investigate the influence of metabolic inhibitors in MCF-7 cells, 2DG was used to inhibit the production of glucose-6-phosphate and the downstream glycolysis pathway [33] (Figure S4A). 6AN was used to inhibit the pentose phosphate pathway (PPP) by blocking 6-phosphogluconate dehydrogenase [34] (Figure S4B). C968 is a glutaminase inhibitor used to block glutamine to glutamate production in the TCA cycle though it does not directly affect glucose metabolism [35,36]. Finally, Rapa was used as an mTOR inhibitor to test its effect on glucose metabolism [37]. The polar metabolite isotopomer profiles of MCF-7 breast cancer cells were analyzed by both targeted and untargeted metabolomics (Figures S5–S10). The targeted data was used as a control for some key central carbon metabolites while the untargeted data was used for discovery. As expected, 2DG induced glucose phosphorylation inhibition resulted the reduction of glycolytic flux compared to the vehicle (DMSO) and other inhibitors (Figures S5–S7). The 13C isotopomer incorporations of the glycolysis intermediates (fructose-6-phosphate, fructose-1,6-bisphosphate, glyceraldehyde-3-phosphate, 3-phosphoglycerate, pyruvate, etc.) were effectively reduced at the 2-h time point (Figure S5) and remained low over 24 h (Figures S6 and S7). The 6AN treatment behaved as expected by inhibiting 6-phosphogluconate dehydrogenase and the accumulation of 13C isotopomers of 6-phosphogluconolactone and 6-phosphogluconate over 24 h [34] (Figures S5–S7). This inhibition resulted in a decrease in 13C isotopomer incorporation into ribose-5-phosphate and PPP metabolism. The targeted isotopomer tracing results also showed that C968 had little effect on glucose metabolism, whereby most of the glycolytic intermediates displayed high levels of [13] C incorporation as compared to the other three metabolic inhibitors. The inhibition of mTOR kinase by Rapa decreased the 13C incorporation into phosphoribosyl pyrophosphate (PRPP) over 24 h (Figures S5–S7). This depletion of the PRPP 13C isotopomers resulted in reduced nucleic acids biosynthesis flux at 2 h and 16 h, which indicated mTOR inhibition can decrease pyrimidine and purine biosynthesis [9,38].
Simultaneously, we collected untargeted metabolomic results for the time-course analysis at 2-h, 16-h, and 24-h drug treatments (Figures S8–S10). The metabolites and their associated 13C isotopomers revealed the same flux alterations as the targeted methods, such as the glycolysis inhibition by 2DG and PPP inhibition via 6-phosphogluconolactonate accumulation by 6AN. Moreover, untargeted metabolic 13C isotope tracer analysis using IsoSearch can be useful to discover novel or unexpected metabolic and lipidomic pathways which are not profiled by targeted methods. In general, the untargeted metabolomics results in MCF-7 breast cancer cells revealed that 13C isotopomers of some phospholipids, amino acids, 1-carbon metabolites, and signaling intermediates were highly dysregulated by the metabolic inhibitors (Figure 3A). The exceedingly increased metabolites implied metabolic isotopomers were altered by drug treatments in MCF-7 cells (Figure 3B–D). For instance, the dramatically elevated 13C-labeled modified amino acids (13C5-3-hydroxy-3-methylglutarate, 13C5-N-acetyl-alanine and 13C5-ornithine) in C968 treated cells indicated that 13C atoms fluxed to the amino acid via urea cycle [39] and amino acid catabolism [40,41] due to glutaminase inhibition; 3-hydroxy-3-methylglutarate (HMG) and N-acetyl-alanine are two amino acids–valine and alanine derivatives which can be catabolized via the urea cycle and increase ornithine 13C isotopomers was observed after C968 treatment; the hepatic glutaminase dysfunction can provide the urea cycle with ammonia, as previously described [42,43]. The highly accumulated fully 13C labeled 6-phosphogluconate by 6AN treatment was due to the inhibition of 6-phosphogluconate dehydrogenase [34] since 6-phosphogluconate links the glycolysis and PPP, and the 6AN inhibits the 13C atom flow from glucose-6-phosphate to pentose phosphate were observed in myelogenous leukemia and pheochromocytoma [44,45].

3.4. Drugs Induce Oxidative Responses in MCF-7 Breast Cancer Cells

The levels of an oxidative stress biomarker, L-cysteine-glutathione disulfide, were different in 2DG and 6AN treated MCF-7 cells (Figure 3A), which could be attributed to metabolic pathways, such as aerobic glycolysis, electron transport chain (ETC), and redox reactions [46]. 2DG inhibited glycolysis so that 13C labeled glucose was blocked from flowing to downstream intermediates (Figure 4A,B). It has been reported that cancer cells are more likely to catabolize glucose through the fermentation pathway and produce large quantities of lactate [47,48]. However, due to 2DG inhibition, 13C[6]-glucose (Figure 4E) was transported into the MCF-7 cells, but did not flux the 13C to lactate (Figure 4F) or TCA cycle intermediates (Figure 4I,J). Nevertheless, the untargeted fluxomics results indicated that 6AN increased the levels of 13C3-lactate (Figure 4F), 13C5-NAD + and 13C10-NAD + at the three time points (Figure 4L). Once 6AN blocked the 13C flux to PPP, flux to the other glucose catabolism pathways were upregulated, and 13C6-glucose was converted to fully labeled lactate with the help of NAD + in MCF-7 cells (Figure 4C,D). The 6AN also increased the flux to fumarate in the TCA cycle over 24 h (Figure 4I) while altering succinate flux differently which may indicate a succinate dehydrogenase (SDH) involved electron transport chain (ETC) disturbance (Figure 4J). SDH catalyzes the conversion of succinate to fumarate and is an important respiratory enzyme involved in redox reactions in the mitochondria [49,50]. The 13C isotopic distribution of succinate by 6AN treatment indicated that succinate was highly labeled at 2 h with predominantly four 13C atoms, but the 13C atoms flowed to fumarate after 16 h and 24 h. In metabolic homeostasis, the energy intake and biomass exchange in cells are maintained at a balanced rate. However, it has been reported that an altered TCA cycle network can induce reactive oxidative stress [51]. Using a DCFDA/H2DCFDA fluorescent assay (see Supplemental Materials section), ROS production in MCF-7 cells was detected at 2-h, 16-h, and 24-h time intervals after drug treatments (Figure S14). At 2 h, the ROS levels in the 6AN treated cells were dramatically increased, while they slowly dropped over the next 22 h (Figure S14A). This indicated that 6AN can stimulate ROS generation in MCF-7 cells over a short time period, but the cells later produce antioxidant chemicals to counteract intracellular ROS. The results indicated that 6AN induced a flux imbalance from succinate to fumarate that can result in the perturbation of ETC and accumulation of ROS, which was also observed in ischemic reperfusion injury in the mouse by invoking oxidative stress [52]. One of the oxidative stress biomarkers [53,54,55], the 13C isotopomers of L-cysteine-glutathione disulfide were increased with 6AN treatment (Figure 4G). In order to defend against intracellular ROS generation, MCF-7 cells utilized reductive species to eliminate oxidative stress and maintain the redox homeostasis. Two antioxidant compounds, carnitine and carnosine, were detected with a dramatically higher 13C isotopomer incorporation with 6AN treatment (Figure 4K,H). The 13C atoms from 13C[6]-glucose flowed to carnosine and carnitine at later time points (16 h and 24 h) which implied that oxidative stress defense was stimulated by the ROS.

3.5. Lipid Isotopomer Regulation by Rapamycin Treatment

Using the 13C glucose isotopic tracer labeling, untargeted time-course isotopomer analysis can also reveal dynamic changes in lipid metabolism. Tracing lipid isotopomers is one of the IsoSearch applications in addition to metabolic isotopomer profiling. The untargeted lipid LC-MS/MS platform is nearly identical to the polar metabolite platform except that C18 reversed phase chromatography is used rather than hydrophilic interaction chromatography (HILIC). Lipids were acquired with their associated isotopomers and IsoSearch was used to profile the lipid fluxes over 2, 16, and 24-h time points (Figures S11–S13). It takes a relatively long time until 13C flows to lipids from 13C[6]-glucose; therefore, we selected 24-h treatment at to investigate how the 2DG and Rapa affect lipid flux. Previous work in our lab has shown that triglycerides (TG) accumulate to high levels with overnight Rapa treatment [26]. Various lipid classes including phospholipids, acylglycerols, sphingolipids, and ceramides and triglycerides were identified using LipidSearch identification software and their 13C isotopomer incorporations were profiled by IsoSearch. The common lipid class isotopomers were altered after 24-h treatment with 2DG or Rapa (Figure 5A). The IsoSearch results illustrated that most of the 13C carbon atoms that flowed into lipids were reduced by 2DG treatment via glycolysis inhibition [56,57] (Figure 5A,B, and Figures S8 and S9). Most phospholipid classes such as phosphatidylcholine (PC) were slightly up-regulated by Rapa treatment through 16 h and then down-regulated after 24-h while the flux of lysophosphatidylcholine (LPC) and phosphatidylglycerol (PG) were down-regulated and triglyceride (TG) was up-regulated at 24 h of treatment [26,58]. Three lipid classes (PC, LPC, and TG) were selected to represent the different regulation with Rapa. Lipids with different fatty acid chain lengths can behave differently even within the same lipid class. Using IsoSearch, the isotopomer ratio of three lipid classes (PC, LPC, and TG) were profiled based on different fatty acid chain lengths (Figure 5B). Lipids containing the same fatty acid chain lengths were selected and their associated isotopomers were summarized separately. The ratio of 13C isotopomer intensity (sum of 13C isotopomer intensities/sum of 13C isotopomer intensities and 12C isotopomer intensity) was used to calculate the 13C isotopic incorporation value. The reduced PC flux is likely a combination of two effects; (1) the altered short fatty acid chains such as palmitate (16:0) and stearate (18:0) and (2) down-regulated isotopomer incorporation of longer fatty acid chains including arachidate (C20:0), lignocerate (C24:0), montanate (C28:0), etc.).
Using IsoSearch, more than thirty 13C carbon atom isotopomers were detected. The isotopic distributions of PC(16:0/16:1), LPC(16:0), and TG(16:0/18:1/18:1) were different in vehicle, 2DG or Rapa treated cells (Figure 5C–E), but all three lipids had a relatively high (M+3) (three 13C atoms) isotopic peak. The (M+3) 13C isotope peak is derived from the lipid head group from glycerol [59,60]. In the 13C[6]-glucose treated experiment, the 13C3-glycerol is the product of glycerol-3-phosphate, which is generated via glycolysis. 2DG blocked glycolysis and also impeded the biosynthesis of many lipids in MCF-7 cells leading to lower levels of the (M+ 3) isotopomer in PC(16:0/16:1, LPC(16:0), and TG(16:0/18:1/18:1). Rapamycin treated cells altered (PC(16:0/16:1), LPC(16:0) and TG(16:0/18:1/18:1)) differently; the 13C isotope distribution pattern of PC(16:0/16:1) was similar when the isotopomers were less than (M+16), the LPC(16:0) isotopomer levels were decreased by Rapa treatment while TG(16:0/18:1/18:1) increased. The higher isotopomers (>ten 13C atoms) decrease with an unchanged (M+3) head group peak in PC(16:0/16:1) implying that Rapa affects the phosphatidylcholine salvage pathway where the fatty acid chains are transferred from diacylglycerol(16:0/16:1) to CDP-choline [61]. Similar to the 2DG treated cells, the isotopomers including the (M+3) head group in LPC(16:0) was decreased with Rapa treatment. LPC is a partial hydrolysis product of phosphatidylcholines, which are catalyzed by phospholipase A2. Considering the unchanged level of (M+3) head group in PC(16:0/16:1), the decrease of (M+3) in LPC(16:0) in Rapa treated cells was the effect of phosphatidylcholine hydrolysis [62,63]. Once the PC hydrolysis reaction was altered by Rapa, the LPC lipids with various fatty acid chain lengths were downregulated in MCF-7 cells. Rapa treated cells showed high intensity of TG(16:0/18:1/18:1) isotopomers, which indicated the TG storage was suppressed by mTOR expression in MCF-7 breast cancer cells. The elevation of TG lipids has also been found in Rapa treated hypertriglyceridemia patients as well as the tsc2 mutated mouse embryonic fibroblast (MEF) cells, which indicated that mTOR suppressed TG accumulation, can also occur in some other cells and patients [26,64,65] We also attempted to incorporate triple labeled stable isotope labeling by amino acids in cell culture (SILAC) proteomics to correlate any identifiable metabolic or lipidomic enzymes to our isotope tracing results (see Supplemental Materials Section, Figures S15 and S16) for multi-omics analyses.

4. Materials and Methods

4.1. Chemicals and Reagents

Methyl tertiary-butyl ether (MTBE) used for lipid and metabolite extraction was purchased from Sigma-Aldrich. LC-MS grade water, LC-MS grade acetonitrile (ACN), and HPLC grade isopropanol (IPA) were purchased from Fisher Chemical. HPLC grade methanol (MeOH) was purchased from Pharmco-Aaper. Mobile phase buffer LC-MS grade formic acid and formic acetate were purchased from Sigma-Aldrich. Fully carbon labeled 13C6 glucose for isotopic labeling was purchased from Cambridge Isotope Laboratories. High glucose DMEM and glucose-free Dulbecco’s modified eagle medium (DMEM) were purchased from Gibco. For the drug treatments, 6-aminonicotinamide, 2-deoxyglucose, compound 968 and rapamycin were purchased from Sigma-Aldrich. A Cadenza CD-C18 column (150 × 2 mm) was purchased from Imtakt, and an XBridge Amide HILIC column (3.5 µm, 4.6 × 100 mm) was purchased from Waters.

4.2. Cell Culture

MCF-7 cells were seeded on 10-cm2 cell culture plates and grown in high glucose DMEM media at 37 °C overnight and reached ~30% confluency. The media was replaced by 13C6 glucose DMEM (13C6 labeled) which prepared by adding 4.5 g/L 13C6-glucose into the glucose-free DMEM or regular 4.5 g/L high glucose DMEM (12C-unlabeled) next day, and four drugs including 2DG, 6AN, Rapa, C968, and DMSO (vehicle) were simultaneously added to the cell culture plates, respectively. The cells were harvested at 2-h, 16-h, or 24-h after the drug treatments.

4.3. Lipid and Metabolite Co-Extraction

The harvested cells were extracted using the methanol/MTBE based extraction method for both metabolites and lipids according to previously published studies [26,66]. The upper phase was collected to obtain the non-polar lipids, and the middle aqueous phase was collected to obtain the polar metabolites; both phases were dried in a SpeedVac at room temperature separately. Polar metabolites and non-polar lipids were stored at −80 °C until mass spectrometry analysis [26,66]. Metabolites were reconstituted in 20 µL of water and lipids were reconstituted in 35 µL of MeOH/IPA (v:v 50:50) just prior to LC-MS/MS analyses.

4.4. LC-MS/MS Based Metabolomic and Lipidomic Analyses

The reconstituted polar metabolite samples split into two equal parts for either targeted or untargeted metabolomics analysis. For targeted metabolomics, 5 µL sample was injected into the 5500 or 6500 QTRAP hybrid triple quadrupole MS (AB/SCIEX) using amide HILIC chromatography and acquired via SRM with a positive/negative ion-switching targeted analysis with 300 unlabeled targets [6,8]. For untargeted metabolomics and lipidomics, the same volume of sample was injected into a high resolution QExactive HF Orbitrap MS (Thermo Fisher Scientific) with positive/negative switching. The source voltages were 4.25 kV in both modes, the source heater temperature was 300 °C and the capillary temperature was set to 250 °C. The resolution was set to 70,000 in MS1 mode and 35,000 in MS2 mode. The top eight ions in each mode were selected for MS2 via data-dependent acquisition (DDA). In MS2 mode, an automatic gain control (AGC) target of 1e6 and a max IT of 85 msec was used. In MS1 mode, the AGC was 5e5 and IT was set to 65 msec. The untargeted metabolomics method adopted the same HPLC conditions as previously published targeted metabolomics, using amide HILIC XBridge (Waters) column with a gradient of 85% B to 2% B over 15 min. A buffer (pH = 9.0: 95% (vol/vol) water, 5% (vol/vol) acetonitrile, 20 mM ammonium hydroxide, 20 mM ammonium acetate); B buffer (100% acetonitrile) via DDA [66]. The 5 µL of the reconstituted lipid samples were injected onto the QExactive Plus LC-MS/MS system using polarity switching / DDA mode with C18 reversed phase chromatography [26,66] as previously published with a gradient of 32% B to 97% B over 25 min. A buffer (pH = 3.5: 39.9% (v/v) water, 60% (v/v) acetonitrile, 10 mM ammonium formate, 0.1% formic acid); B buffer (89.9% (v/v) isopropanol, 10.0% (v/v) acetonitrile, 10 mM ammonium formate, 0.1% formic acid). The untargeted lipidomics MS acquisition parameters were the same as above for untargeted metabolomics. We found that amide HILIC was the most robust and reproducible chromatographic media for polar small molecule metabolites, while intact non-polar lipids were best separated and detected using C18 reversed phase chromatography columns [6,8,26].

5. Conclusions

Isotopic tracer flux analyses typically report the detection of molecules present in well-known metabolic pathways. In this report, we have described a strategy to interpret metabolite and lipid isotopomers of polar metabolites and non-polar lipids identified from untargeted high resolution LC-MS/MS analyses, and subsequently re-identified and quantified in 13C isotopically labeled samples. The IsoSearch, strategy profiles the distribution of each identified molecule’s isotope profile, and providing information for subsequent polar metabolic flux modeling. We also applied our technique to a set of non-polar lipids to reveal lipid class regulation, fatty acid chain composition changes, and individual lipid isotopic distribution alterations. This R script derived isotope labeling strategy is also released as a user-friendly component of the commercial software Scaffold Elements versions 2.0 and later (Proteome Software) and requires high resolution and high mass accuracy MS and MS/MS data. The ability to perform isotope tracer analyses using polarity switching has the added benefits of conserving precious sample while also improving sample throughput and a longer chromatographic run can compensate for any potential data loss due to MS acquisition cycling. In general, IsoSearch provides metabolite/lipid isotopomer tracing capabilities to unbiased untargeted metabolomics and lipidomics to facilitate the large-scale global exploration of isotopic fluxomics from any biological source, whereby metabolites and lipids can be extracted and profiled by LC-MS/MS. Prior to isotope tracer analysis, our LC-MS/MS polarity switching platform can identify and annotate ~300–450 metabolites and ~ 800–1500 lipids for many biological sample types. The ability to trace isotopomers from any identified molecule is a unique feature of IsoSearch (and Scaffold Elements). In addition, only molecules that have been confidently identified and annotated based on MS/MS data can be further used for IsoSearch isotopomer profile tracing.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2409-9279/3/3/54/s1: Figures S1–S17, ROS detection assay and SILAC methodology and multi-omics results. Metabolite and Lipid results files have been published on the Zenodo data repository, 10.5281/zenodo.3908617.

Author Contributions

H.H., J.M.A., and S.L. conceptualized the methodology. H.H. and P.S. designed and implemented the algorithm and software. H.H., M.Y., and J.M.A. performed mass spectrometry analyses. H.H. and J.M.A. analyzed data and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

We acknowledge the National Institutes of Health (NIH) grants 5P01CA120964 (J.M.A.), 5P30CA006516 (J.M.A.). We also acknowledge support from the BIDMC Capital Equipment fund (J.M.A.).

Acknowledgments

We thank Gerburg Wulf at BIDMC/HMS in Boston for providing MCF-7 breast cancer cells.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sauer, U. Metabolic networks in motion: 13C-based flux analysis. Mol. Syst. Biol. 2006, 2, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Buescher, J.M.; Antoniewicz, M.R.; Boros, L.G.; Burgess, S.C.; Brunengraber, H.; Clish, C.B.; DeBerardinis, R.J.; Feron, O.; Frezza, C.; Ghesquiere, B.; et al. A roadmap for interpreting (13)C metabolite labeling patterns from cells. Curr. Opin. Biotechnol. 2015, 34, 189–201. [Google Scholar] [CrossRef] [PubMed]
  3. Schwender, J.; Goffman, F.; Ohlrogge, J.B.; Shachar-Hill, Y. Rubisco without the calvin cycle improves the carbon efficiency of developing green seeds. Nature 2004, 432, 779–782. [Google Scholar] [CrossRef] [PubMed]
  4. Sriram, G.; Fulton, D.B.; Iyer, V.V.; Peterson, J.M.; Zhou, R.; Westgate, M.E.; Spalding, M.H.; Shanks, J.V. Quantification of compartmented metabolic fluxes in developing soybean embryos by employing biosynthetically directed fractional (13)C labeling, two-dimensional [(13)C, (1)H] nuclear magnetic resonance, and comprehensive isotopomer balancing. Plant. Physiol. 2004, 136, 3043–3057. [Google Scholar] [CrossRef] [Green Version]
  5. Hellerstein, M.K.; Murphy, E. Stable isotope-mass spectrometric measurements of molecular fluxes in vivo: Emerging applications in drug development. Curr. Opin. Mol. Ther. 2004, 6, 249–264. [Google Scholar]
  6. Yuan, M.; Breitkopf, S.B.; Yang, X.; Asara, J.M. A positive/negative ion-switching, targeted mass spectrometry-based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nat. Protoc. 2012, 7, 872–881. [Google Scholar] [CrossRef] [Green Version]
  7. Reisz, J.A.; D’Alessandro, A. Measurement of metabolic fluxes using stable isotope tracers in whole animals and human patients. Curr. Opin. Clin. Nutr. Metab. Care 2017, 20, 366–374. [Google Scholar] [CrossRef]
  8. Yuan, M.; Kremer, D.M.; Huang, H.; Breitkopf, S.B.; Ben-Sahra, I.; Manning, B.D.; Lyssiotis, C.A.; Asara, J.M. Ex vivo and in vivo stable isotope labelling of central carbon metabolism and related pathways with analysis by LC-MS/MS. Nat. Protoc. 2019, 14, 313–330. [Google Scholar] [CrossRef]
  9. Ben-Sahra, I.; Howell, J.J.; Asara, J.M.; Manning, B.D. Stimulation of de novo pyrimidine synthesis by growth signaling through MTOR and S6K1. Science 2013, 339, 1323–1328. [Google Scholar] [CrossRef] [Green Version]
  10. Son, J.; Lyssiotis, C.A.; Ying, H.; Wang, X.; Hua, S.; Ligorio, M.; Perera, R.M.; Ferrone, C.R.; Mullarky, E.; Shyh-Chang, N.; et al. Glutamine supports pancreatic cancer growth through a KRAS-regulated metabolic pathway. Nature 2013, 496, 101–105. [Google Scholar] [CrossRef]
  11. Ying, H.; Kimmelman, A.C.; Lyssiotis, C.A.; Hua, S.; Chu, G.C.; Fletcher-Sananikone, E.; Locasale, J.W.; Son, J.; Zhang, H.; Coloff, J.L.; et al. Oncogenic kras maintains pancreatic tumors through regulation of anabolic glucose metabolism. Cell 2012, 149, 656–670. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Sousa, C.M.; Biancur, D.E.; Wang, X.; Halbrook, C.J.; Sherman, M.H.; Zhang, L.; Kremer, D.; Hwang, R.F.; Witkiewicz, A.K.; Ying, H.; et al. Pancreatic stellate cells support tumour metabolism through autophagic alanine secretion. Nature 2016, 536, 479–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Nicolay, B.N.; Gameiro, P.A.; Tschop, K.; Korenjak, M.; Heilmann, A.M.; Asara, J.M.; Stephanopoulos, G.; Iliopoulos, O.; Dyson, N.J. Loss of RBF1 changes glutamine catabolism. Genes Dev. 2013, 27, 182–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wiechert, W. 13C metabolic flux analysis. Metab. Eng. 2001, 3, 195–206. [Google Scholar] [CrossRef]
  15. Zamboni, N.; Fendt, S.-M.; Ruhl, M.; Sauer, U. (13)C-based metabolic flux analysis. Nat. Protoc. 2009, 4, 878–892. [Google Scholar] [CrossRef]
  16. Nanchen, A.; Fuhrer, T.; Sauer, U. Determination of metabolic flux ratios from 13C-experiments and gas chromatography-mass spectrometry data: Protocol and principles. Methods Mol. Biol. 2007, 358, 177–197. [Google Scholar]
  17. Crown, S.B.; Long, C.P.; Antoniewicz, M.R. Optimal tracers for parallel labeling experiments and (13)C metabolic flux analysis: A new precision and synergy scoring system. Metab. Eng. 2016, 38, 10–18. [Google Scholar] [CrossRef] [Green Version]
  18. Kogadeeva, M.; Zamboni, N. SUMOFLUX: A generalized method for targeted 13C metabolic flux ratio analysis. PLoS Comput. Biol. 2016, 12, e1005109. [Google Scholar] [CrossRef] [Green Version]
  19. Weitzel, M.; Nöh, K.; Dalman, T.; Niedenführ, S.; Stute, B.; Wiechert, W. 13CFLUX2—High-performance software suite for (13)C-metabolic flux analysis. Bioinformatics 2013, 29, 143–145. [Google Scholar] [CrossRef]
  20. Huang, X.; Chen, Y.-J.; Cho, K.; Nikolskiy, I.; Crawford, P.A.; Patti, G.J. X13CMS: Global tracking of isotopic labels in untargeted metabolomics. Anal. Chem. 2014, 86, 1632–1639. [Google Scholar] [CrossRef]
  21. Noh, K.; Droste, P.; Wiechert, W. Visual workflows for 13C-metabolic flux analysis. Bioinformatics 2015, 31, 346–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Puchalska, P.; Huang, X.; Martin, S.E.; Han, X.; Patti, G.J.; Crawford, P.A. Isotope tracing untargeted metabolomics reveals macrophage polarization-state-specific metabolic coordination across intracellular compartments. iScience 2018, 9, 298–313. [Google Scholar] [CrossRef] [Green Version]
  23. Puchalska, P.; Martin, S.E.; Huang, X.; Lengfeld, J.E.; Daniel, B.; Graham, M.J.; Han, X.; Nagy, L.; Patti, G.J.; Crawford, P.A. Hepatocyte-macrophage acetoacetate shuttle protects against tissue fibrosis. Cell Metab. 2019, 29, 383–398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chen, Y.-J.; Mahieu, N.G.; Huang, X.; Singh, M.; Crawford, P.A.; Johnson, S.L.; Gross, R.W.; Schaefer, J.; Patti, G.J. Lactate metabolism is associated with mammalian mitochondria. Nat. Chem. Biol. 2016, 12, 937. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Seitzer, P.M.; Searle, B.C. Incorporating in-source fragment information improves metabolite identification accuracy in untargeted LC–MS data sets. J. Proteome Res. 2019, 18, 791–796. [Google Scholar] [CrossRef]
  26. Breitkopf, S.B.; Ricoult, S.J.H.; Yuan, M.; Xu, Y.; Peake, D.A.; Manning, B.D.; Asara, J.M. A relative quantitative positive/negative ion switching method for untargeted lipidomics via high resolution LC-MS/MS from any biological source. Metabolomics 2017, 13, 30. [Google Scholar] [CrossRef] [Green Version]
  27. Gibbons, H.; Brennan, L. Metabolomics as a tool in the identification of dietary biomarkers. Proc. Nutr. Soc. 2017, 76, 42–53. [Google Scholar] [CrossRef]
  28. Schrimpe-Rutledge, A.C.; Codreanu, S.G.; Sherrod, S.D.; McLean, J.A. Untargeted metabolomics strategies-challenges and emerging directions. J. Am. Soc. Mass Spectrom. 2016, 27, 1897–1905. [Google Scholar] [CrossRef] [Green Version]
  29. Nam, I.; Lee, J.K.; Nam, H.G.; Zare, R.N. Abiotic production of sugar phosphates and uridine ribonucleoside in aqueous microdroplets. Proc. Natl. Acad. Sci. USA 2017, 114, 12396–12400. [Google Scholar] [CrossRef] [Green Version]
  30. Kaushik, A.K.; DeBerardinis, R.J. Applications of metabolomics to study cancer metabolism. Biochim. Biophys. Acta Rev. Cancer 2018, 1870, 2–14. [Google Scholar] [CrossRef]
  31. Kang, Y.P.; Ward, N.P.; DeNicola, G.M. Recent advances in cancer metabolism: A technological perspective. Exp. Mol. Med. 2018, 50, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Sousa, C.M.; Kimmelman, A.C. The complex landscape of pancreatic cancer metabolism. Carcinogenesis 2014, 35, 1441–1450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wick, A.N.; Drury, D.R.; Nakada, H.I.; Wolfe, J.B. Localization of the primary metabolic block produced by 2-deoxyglucose. J. Biol. Chem. 1957, 224, 963–969. [Google Scholar] [PubMed]
  34. Salimi, A.; Paeezi, M.; Yousefsani, B.S.; Shadnia, S.; Hassanian-Moghaddam, H.; Zamani, N.; Pourahmad, J. Inhibition of glucose-6-phosphate dehydrogenase protects hepatocytes from aluminum phosphide-induced toxicity. Pestic. Biochem. Physiol. 2017, 143, 141–146. [Google Scholar] [CrossRef] [PubMed]
  35. Yuan, L.; Sheng, X.; Clark, L.H.; Zhang, L.; Guo, H.; Jones, H.M.; Willson, A.K.; Gehrig, P.A.; Zhou, C.; Bae-Jump, V.L. Glutaminase inhibitor compound 968 inhibits cell proliferation and sensitizes paclitaxel in ovarian cancer. Am. J. Transl. Res. 2016, 8, 4265–4277. [Google Scholar]
  36. Wang, D.; Meng, G.; Zheng, M.; Zhang, Y.; Chen, A.; Wu, J.; Wei, J. The glutaminase-1 inhibitor 968 enhances dihydroartemisinin-mediated antitumor efficacy in hepatocellular carcinoma cells. PLoS ONE 2016, 11, e0166423. [Google Scholar] [CrossRef]
  37. Caron, A.; Richard, D.; Laplante, M. The roles of MTOR complexes in lipid metabolism. Annu. Rev. Nutr. 2015, 35, 321–348. [Google Scholar] [CrossRef]
  38. Ben-Sahra, I.; Hoxhaj, G.; Ricoult, S.J.H.; Asara, J.M.; Manning, B.D. MTORC1 induces purine synthesis through control of the mitochondrial tetrahydrofolate cycle. Science 2016, 351, 728–733. [Google Scholar] [CrossRef] [Green Version]
  39. Wright, P.; Felskie, A.; Anderson, P. Induction of ornithine-urea cycle enzymes and nitrogen metabolism and excretion in rainbow trout (oncorhynchus mykiss) during early life stages. J. Exp. Biol. 1995, 198, 127–135. [Google Scholar]
  40. Camardella, L.; Damonte, G.; Carratore, V.; Benatti, U.; Tonetti, M.; Moneti, G. Glucose 6-phosphate dehydrogenase from human erythrocytes: Identification of N-acetyl-alanine at the N-terminus of the mature protein. Biochem. Biophys. Res. Commun. 1995, 207, 331–338. [Google Scholar] [CrossRef]
  41. Dena, R.; Fabbro, M.; Rigoni, F. Formation and utilization of 3-hydroxy-3-methylglutarate in liver mitochondria of starved and streptozotocin-diabetic rats. Biochem. J. 1978, 172, 371–375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Moorman, A.F.; de Boer, P.A.; Watford, M.; Dingemanse, M.A.; Lamers, W.H. Hepatic glutaminase MRNA is confined to part of the urea cycle domain in the adult rodent liver lobule. FEBS Lett. 1994, 356, 76–80. [Google Scholar] [CrossRef] [Green Version]
  43. Häussinger, D.; Gerok, W. Regulation of hepatic glutamate metabolism. Eur. J. Biochem. 1984, 143, 491–497. [Google Scholar] [PubMed]
  44. Chen, Y.; Xu, Q.; Ji, D.; Wei, Y.; Chen, H.; Li, T.; Wan, B.; Yuan, L.; Huang, R.; Chen, G. Inhibition of pentose phosphate pathway suppresses acute myelogenous leukemia. Tumour Biol. 2016, 37, 6027–6034. [Google Scholar] [CrossRef] [PubMed]
  45. Davis, L.H.; Kauffman, F.C. Metabolism via the pentose phosphate pathway in rat pheochromocytoma PC12 cells: Effects of nerve growth factor and 6-aminonicotinamide. Neurochem. Res. 1987, 12, 521–527. [Google Scholar] [CrossRef]
  46. Panieri, E.; Santoro, M.M. ROS homeostasis and metabolism: A dangerous liason in cancer cells. Cell Death Dis. 2016, 7, e2253. [Google Scholar] [CrossRef]
  47. Liberti, M.V.; Locasale, J.W. The warburg effect: How does it benefit cancer cells? Trends Biochem. Sci. 2016, 41, 211–218. [Google Scholar] [CrossRef] [Green Version]
  48. Burns, J.S.; Manda, G. Metabolic pathways of the warburg effect in health and disease: Perspectives of choice, chain or chance. Int. J. Mol. Sci. 2017, 18, 2755. [Google Scholar] [CrossRef] [Green Version]
  49. Mills, E.L.; Kelly, B.; Logan, A.; Costa, A.S.H.; Varma, M.; Bryant, C.E.; Tourlomousis, P.; Dabritz, J.H.M.; Gottlieb, E.; Latorre, I.; et al. Succinate dehydrogenase supports metabolic repurposing of mitochondria to drive inflammatory macrophages. Cell 2016, 167, 457–470. [Google Scholar] [CrossRef] [Green Version]
  50. Bezawork-Geleta, A.; Rohlena, J.; Dong, L.; Pacak, K.; Neuzil, J. Mitochondrial complex II: At the crossroads. Trends Biochem. Sci. 2017, 42, 312–325. [Google Scholar] [CrossRef]
  51. Mailloux, R.J.; Beriault, R.; Lemire, J.; Singh, R.; Chenier, D.R.; Hamel, R.D.; Appanna, V.D. The tricarboxylic acid cycle, an ancient metabolic network with a novel twist. PLoS ONE 2007, 2, e690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Chouchani, E.T.; Pell, V.R.; Gaude, E.; Aksentijevic, D.; Sundier, S.Y.; Robb, E.L.; Logan, A.; Nadtochiy, S.M.; Ord, E.N.J.; Smith, A.C.; et al. Ischaemic accumulation of succinate controls reperfusion injury through mitochondrial ROS. Nature 2014, 515, 431–435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Heman-Ackah, S.E.; Juhn, S.K.; Huang, T.C.; Wiedmann, T.S. A combination antioxidant therapy prevents age-related hearing loss in C57BL/6 mice. Otolaryngol. Head. Neck Surg. 2010, 143, 429–434. [Google Scholar] [CrossRef] [PubMed]
  54. Gao, X.-H.; Qanungo, S.; Pai, H.V.; Starke, D.W.; Steller, K.M.; Fujioka, H.; Lesnefsky, E.J.; Kerner, J.; Rosca, M.G.; Hoppel, C.L.; et al. Aging-dependent changes in rat heart mitochondrial glutaredoxins—Implications for redox regulation. Redox Biol. 2013, 1, 586–598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Berkeley, L.I.; Cohen, J.F.; Crankshaw, D.L.; Shirota, F.N.; Nagasawa, H.T. Hepatoprotection by L-cysteine-glutathione mixed disulfide, a sulfhydryl-modified prodrug of glutathione. J. Biochem. Mol. Toxicol. 2003, 17, 95–97. [Google Scholar] [CrossRef]
  56. Hong, S.Y.; Yu, F.-X.; Luo, Y.; Hagen, T. Oncogenic activation of the PI3K/Akt pathway promotes cellular glucose uptake by downregulating the expression of thioredoxin-interacting protein. Cell. Signal. 2016, 28, 377–383. [Google Scholar] [CrossRef]
  57. Cao, B.-Y.; Li, R.; Tian, H.-H.; Ma, Y.-J.; Hu, X.-G.; Jia, N.; Wang, Y.-Y. PI3K-GLUT4 signal pathway associated with effects of EX-B3 electroacupuncture on hyperglycemia and insulin resistance of T2DM rats. Evid. Based Complement. Alternat. Med. 2016, 2016, 7914387. [Google Scholar] [CrossRef] [Green Version]
  58. Ricoult, S.J.H.; Yecies, J.L.; Ben-Sahra, I.; Manning, B.D. Oncogenic PI3K and K-Ras stimulate de novo lipid synthesis through MTORC1 and SREBP. Oncogene 2016, 35, 1250–1260. [Google Scholar] [CrossRef] [Green Version]
  59. Wang, X.; Xiong, X.; Sa, N.; Roje, S.; Chen, S. Metabolic engineering of enhanced Glycerol-3-phosphate synthesis to increase lipid production in synechocystis sp. PCC 6803. Appl. Microbiol. Biotechnol. 2016, 100, 6091–6101. [Google Scholar] [CrossRef]
  60. Xue, L.-L.; Chen, H.-H.; Jiang, J.-G. Implications of glycerol metabolism for lipid production. Prog. Lipid Res. 2017, 68, 12–25. [Google Scholar] [CrossRef]
  61. Moessinger, C.; Klizaite, K.; Steinhagen, A.; Philippou-Massier, J.; Shevchenko, A.; Hoch, M.; Ejsing, C.S.; Thiele, C. Two different pathways of phosphatidylcholine synthesis, the kennedy pathway and the lands cycle, differentially regulate cellular triacylglycerol storage. BMC Cell Biol. 2014, 15, 43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Lim, C.W.; Kim, B.H.; Kim, I.-H.; Lee, M.-W. Modeling and optimization of phospholipase A(1)-catalyzed hydrolysis of phosphatidylcholine using response surface methodology for lysophosphatidylcholine production. Biotechnol. Prog. 2015, 31, 35–41. [Google Scholar] [CrossRef] [PubMed]
  63. Tsao, F.H.; Zachman, R.D. Phosphatidylcholine-lysophosphatidylcholine cycle pathway enzymes in rabbit lung. II. marked differences in the effect of gestational age on activity compared to the CDP-choline pathway. Pediatr. Res. 1977, 11, 858–861. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kasiske, B.L.; de Mattos, A.; Flechner, S.M.; Gallon, L.; Meier-Kriesche, H.-U.; Weir, M.R.; Wilkinson, A. Mammalian target of rapamycin inhibitor dyslipidemia in kidney transplant recipients. Am. J. Transplant. 2008, 8, 1384–1392. [Google Scholar] [CrossRef]
  65. Morrisett, J.D.; Abdel-Fattah, G.; Hoogeveen, R.; Mitchell, E.; Ballantyne, C.M.; Pownall, H.J.; Opekun, A.R.; Jaffe, J.S.; Oppermann, S.; Kahan, B.D. Effects of sirolimus on plasma lipids, lipoprotein levels, and fatty acid metabolism in renal transplant patients. J. Lipid Res. 2002, 43, 1170–1180. [Google Scholar] [CrossRef] [Green Version]
  66. Breitkopf, S.B.; Taveira, M.D.O.; Yuan, M.; Wulf, G.M.; Asara, J.M. Serial-omics of P53-/-, Brca1-/- mouse breast tumor and normal mammary gland. Sci. Rep. 2017, 7, 14503. [Google Scholar] [CrossRef] [Green Version]
Mps 03 00054 g001 550
Figure 1. The workflow of untargeted fluxomic analysis. Two sets of cells are cultured in 12C-glucose or 13C[6]-glucose containing media for isotope labeling. Both non-polar lipids and polar metabolites were extracted using MTBE/methanol, and samples were acquired by high resolution LC-MS/MS via data-dependent acquisition (DDA) with polarity switching. The metabolites and lipids of unlabeled samples were then identified using commercial software and serve as references for labeled samples. The isotopically labeled samples were matched against the references using the in-house R script IsoSearch or the commercial software Scaffold Elements.
Figure 1. The workflow of untargeted fluxomic analysis. Two sets of cells are cultured in 12C-glucose or 13C[6]-glucose containing media for isotope labeling. Both non-polar lipids and polar metabolites were extracted using MTBE/methanol, and samples were acquired by high resolution LC-MS/MS via data-dependent acquisition (DDA) with polarity switching. The metabolites and lipids of unlabeled samples were then identified using commercial software and serve as references for labeled samples. The isotopically labeled samples were matched against the references using the in-house R script IsoSearch or the commercial software Scaffold Elements.
Mps 03 00054 g001
Mps 03 00054 g002 550
Figure 2. Validation of untargeted high-resolution fluxomics. (A) Uridine triphosphate (UTP) and its associated isotopomers were detected using both targeted and untargeted metabolomics. The abundant M+5 peak of UTP indicates the ribose ring comes from the 13C[6]-glucose. Both the targeted fluxomics via selected reaction monitoring (SRM) (left) and untargeted high-resolution isotopomer tracing (IsoSearch) (right) show a similar isotopomer pattern after 13C[6]-glucose tracing for UTP. (B) The phospholipid PC(16:0/16:0) and its associated isotopomers were profiled using untargeted lipid fluxomics. The overlap of lipid isotopomer peaks with high carbon numbers and other lipid peaks are distinguished using the IsoSearch strategy.
Figure 2. Validation of untargeted high-resolution fluxomics. (A) Uridine triphosphate (UTP) and its associated isotopomers were detected using both targeted and untargeted metabolomics. The abundant M+5 peak of UTP indicates the ribose ring comes from the 13C[6]-glucose. Both the targeted fluxomics via selected reaction monitoring (SRM) (left) and untargeted high-resolution isotopomer tracing (IsoSearch) (right) show a similar isotopomer pattern after 13C[6]-glucose tracing for UTP. (B) The phospholipid PC(16:0/16:0) and its associated isotopomers were profiled using untargeted lipid fluxomics. The overlap of lipid isotopomer peaks with high carbon numbers and other lipid peaks are distinguished using the IsoSearch strategy.
Mps 03 00054 g002
Mps 03 00054 g003 550
Figure 3. An overview of untargeted fluxomics analysis from drug treated MCF-7 breast cancer cells reveals the changes in metabolites and their associated isotopomers. (A) A heat map of untargeted fluxomics results displays various metabolites and their isotopomer alterations after drug treatments with 2DG, Rapa, 6AN and DMSO vehicle control. The MA scatterplot of (B) the top 10 increased (red) and decreased (blue) metabolites induced by 2DG at 24 h; and (C) the top 10 increased (red) and decreased (blue) metabolites induced by 6AN at 24 h; and (D) the top 10 increased (red) and decreased (blue) metabolites induced by rapamycin at 24 h.
Figure 3. An overview of untargeted fluxomics analysis from drug treated MCF-7 breast cancer cells reveals the changes in metabolites and their associated isotopomers. (A) A heat map of untargeted fluxomics results displays various metabolites and their isotopomer alterations after drug treatments with 2DG, Rapa, 6AN and DMSO vehicle control. The MA scatterplot of (B) the top 10 increased (red) and decreased (blue) metabolites induced by 2DG at 24 h; and (C) the top 10 increased (red) and decreased (blue) metabolites induced by 6AN at 24 h; and (D) the top 10 increased (red) and decreased (blue) metabolites induced by rapamycin at 24 h.
Mps 03 00054 g003
Mps 03 00054 g004 550
Figure 4. The 2DG and 6AN affect metabolic fluxes of 13C-glucose in MCF-7 cells differently. The bar graphs indicate the ratio of 13C incorporation for the isotopomers of each metabolite via IsoSearch. (A) 2DG inhibits the formation of glucose 6-phosphate and 13C flux to glycolysis. The circles represent the carbon atoms of each polar metabolite, and 13C-labeled carbons are filled in red. (B) 2DG inhibits 13C[6]-glucose metabolic flux and its downstream pathways. The arrow thickness represents the relative flux ratio through the metabolic pathways. (C) The 6AN inhibits 6-phosphogluconate dehydrogenase and 13C flux to PPP. (D) The 6AN drives the 13C[6]-glucose metabolic flux to aerobic glycolysis and oxidative stress. (E) Glucose, (F) Lactate, (G) Cysteine-glutathione disulfide, (H) carnosine, (I) fumarate, (J) succinate, (K) carnitine, (L) NAD+ isotopomer incorporations are altered by 2DG or 6AN at 2 h, 16 h, and 24 h.
Figure 4. The 2DG and 6AN affect metabolic fluxes of 13C-glucose in MCF-7 cells differently. The bar graphs indicate the ratio of 13C incorporation for the isotopomers of each metabolite via IsoSearch. (A) 2DG inhibits the formation of glucose 6-phosphate and 13C flux to glycolysis. The circles represent the carbon atoms of each polar metabolite, and 13C-labeled carbons are filled in red. (B) 2DG inhibits 13C[6]-glucose metabolic flux and its downstream pathways. The arrow thickness represents the relative flux ratio through the metabolic pathways. (C) The 6AN inhibits 6-phosphogluconate dehydrogenase and 13C flux to PPP. (D) The 6AN drives the 13C[6]-glucose metabolic flux to aerobic glycolysis and oxidative stress. (E) Glucose, (F) Lactate, (G) Cysteine-glutathione disulfide, (H) carnosine, (I) fumarate, (J) succinate, (K) carnitine, (L) NAD+ isotopomer incorporations are altered by 2DG or 6AN at 2 h, 16 h, and 24 h.
Mps 03 00054 g004
Mps 03 00054 g005 550
Figure 5. Lipid 13C isotopic flux changes are profiled using untargeted LC-MS/MS with polarity switching. (A) The 2DG and Rapa dysregulate various lipid classes compared to the control (DMSO). (B) The 2DG and Rapa alter the phosphatidylcholine (PC), lysophosphatidylcholine (LPC), and triglyceride (TG) lipid classes with various fatty acid chains differently. (C) The 2DG decreases all of the PC(16:0/16:1) isotopomer levels, but Rapa only decreases PC(16:0/16:1) flux through higher 13C isotopes (>M+11). (D) Both 2DG and Rapa down-regulated the levels of LPC(16:0) 13C isotopomers. (E) TG(16:0/18:1/18:1) isotopomer levels decreased with 2DG treatment, but increased with Rapa.
Figure 5. Lipid 13C isotopic flux changes are profiled using untargeted LC-MS/MS with polarity switching. (A) The 2DG and Rapa dysregulate various lipid classes compared to the control (DMSO). (B) The 2DG and Rapa alter the phosphatidylcholine (PC), lysophosphatidylcholine (LPC), and triglyceride (TG) lipid classes with various fatty acid chains differently. (C) The 2DG decreases all of the PC(16:0/16:1) isotopomer levels, but Rapa only decreases PC(16:0/16:1) flux through higher 13C isotopes (>M+11). (D) Both 2DG and Rapa down-regulated the levels of LPC(16:0) 13C isotopomers. (E) TG(16:0/18:1/18:1) isotopomer levels decreased with 2DG treatment, but increased with Rapa.
Mps 03 00054 g005
Table
Table 1. List of IsoSearch() functions.
Table 1. List of IsoSearch() functions.
FunctionParametersDescription
Flux_result()input_negative, input_positive, referInput, score = 0.61 (default)main function to generate the final result of the flux analysis
ref_13C_neg()refInput, refOutputcreate the reference list using negative mode unlabeled samples (12C)
ref_13C_pos()refInput, refOutputcreate the reference list using positive mode unlabeled samples (12C)
msMatch()mzFile, refilesearching process
resultWrap_neg()inputFile_neg, refInputwrapping the flux 13C-labeling results of negative ion mode samples
resultWrap_pos()inputFile_pos, refInputwrapping the flux 13C-labeling results of positive ion mode samples
Sgrade()inputResultscoring function using Equation (1) to calculate the filtering scores
refLibPos_ls()refInputpreserved function for lipidSearch
refLibNeg_ls()refInputpreserved function for lipidSearch
refLibPos_ele()refInputpreserved function for Elements
refLibNeg_ele()refInputpreserved function for Elements
IsoSearch program built-in functions. IsoSearch() uses the 13C-labeling peak lists and 12C-unlabeled reference list for flux search processing.
Table
Table 2. Description of IsoSearch() output table (.csv).
Table 2. Description of IsoSearch() output table (.csv).
Column HeadingDescription
mz1mass to charge ratio of the searched file
rt1Retention time value of the data file
Intensityfeature peak intensity
mz2m/z value of the 13C-labeled experimental file
rt2retention time value of the data file
Metabolite/Lipidname of the feature
fattyAcid, lipidClass, lipidFormfatty acid chain, lipid class and chemical formula of the lipid (lipid fluxomics only)
Accessionaccession number assigned by the database
Theoretical_mzisotopomer theoretical m/z value
Adductadduct ion (+H, −H, etc.)
Charge (z)ion charge mode
Annotationisotopomer notation (M+1, M+2, M+3, etc.)
Scorescore used for feature screening
Δmz_ppmthe difference between experimental and theoretical m/z in ppm
GradesQuality associated with the score where A is best
The fluxomics results generated by IsoSearch() are presented in a spreadsheet with headings.

Share and Cite

MDPI and ACS Style

Huang, H.; Yuan, M.; Seitzer, P.; Ludwigsen, S.; Asara, J.M. IsoSearch: An Untargeted and Unbiased Metabolite and Lipid Isotopomer Tracing Strategy from HR-LC-MS/MS Datasets. Methods Protoc. 2020, 3, 54. https://0-doi-org.brum.beds.ac.uk/10.3390/mps3030054

AMA Style

Huang H, Yuan M, Seitzer P, Ludwigsen S, Asara JM. IsoSearch: An Untargeted and Unbiased Metabolite and Lipid Isotopomer Tracing Strategy from HR-LC-MS/MS Datasets. Methods and Protocols. 2020; 3(3):54. https://0-doi-org.brum.beds.ac.uk/10.3390/mps3030054

Chicago/Turabian Style

Huang, He, Min Yuan, Phillip Seitzer, Susan Ludwigsen, and John M. Asara. 2020. "IsoSearch: An Untargeted and Unbiased Metabolite and Lipid Isotopomer Tracing Strategy from HR-LC-MS/MS Datasets" Methods and Protocols 3, no. 3: 54. https://0-doi-org.brum.beds.ac.uk/10.3390/mps3030054

Article Metrics

Back to TopTop