Analytical Biochemistry 475 (2015) 4–13
Contents lists available at ScienceDirect
Analytical Biochemistry journal homepage: www.elsevier.com/locate/yabio
Alkaline conditions in hydrophilic interaction liquid chromatography for intracellular metabolite quantification using tandem mass spectrometry Attila Teleki 1, Andrés Sánchez-Kopper 1, Ralf Takors ⇑ Institute of Biochemical Engineering, University of Stuttgart, 70569 Stuttgart, Germany
a r t i c l e
i n f o
Article history: Received 17 September 2014 Received in revised form 8 December 2014 Accepted 5 January 2015 Available online 16 January 2015 Keywords: Metabolic profiling HILIC LC–MS/MS
a b s t r a c t Modeling of metabolic networks as part of systems metabolic engineering requires reliable quantitative experimental data of intracellular concentrations. The hydrophilic interaction liquid chromatography–e lectrospray ionization–tandem mass spectrometry (HILIC–ESI–MS/MS) method was used for quantitative profiling of more than 50 hydrophilic key metabolites of cellular metabolism. Without prior derivatization, sugar phosphates, organic acids, nucleotides, and amino acids were measured under alkaline and acidic mobile phase conditions with pre-optimized multiple reaction monitoring (MRM) transitions. Irrespective of the polarity mode of the acquisition method used, alkaline conditions achieved the best quantification limits and linear dynamic ranges. Fully 90% of the analyzed metabolites presented detection limits better than 0.5 pmol (on column), and 70% presented 1.5-fold higher signal intensities under alkaline mobile phase conditions. The quality of the method was further demonstrated by absolute quantification of selected metabolites in intracellular extracts of Escherichia coli. In addition, quantification bias caused by matrix effects was investigated by comparison of calibration strategies: standard-based external calibration, isotope dilution, and standard addition with internal standards. Here, we recommend the use of alkaline mobile phase with polymer-based zwitterionic hydrophilic interaction chromatography (ZIC–pHILIC) as the most sensitive scenario for absolute quantification for a broad range of metabolites. Ó 2015 Elsevier Inc. All rights reserved.
Metabolomics addresses the identification and quantification of small molecule metabolites, in essence the reactants of biological systems, and enables the comprehensive analysis of complex biochemical networks systemic studies. Whereas the genome and proteome represent upstream biochemical events, the metabolome reveals underlying regulatory mechanisms and therefore reflects, most closely, the actual cellular physiological state [1–3]. As such, metabolome analysis can be subdivided into targeted and nontargeted approaches. Nontargeted high-throughput applications such as metabolic fingerprinting aim at the rapid and qualitative or semiquantitative analysis of whole-cell metabolic patterns, reducing the analytical effort considerably [4–6]. On the contrary, targeted applications such as quantitative metabolic profiling are focused on the nonbiased absolute quantification of a subset of metabolites of predefined metabolic pathways or classes of compounds [2,4,6]. In this context, monitoring of intracellular ⇑ Corresponding author. Fax: +49 0711 685 65164. 1
E-mail address: [email protected] (R. Takors). Contributed equally to this manuscript.
http://dx.doi.org/10.1016/j.ab.2015.01.002 0003-2697/Ó 2015 Elsevier Inc. All rights reserved.
metabolism dynamics as a result of rapid changes in the extracellular environment within stimulus–response experiments can be used to investigate in vivo enzyme kinetics offering insights into underlying regulatory mechanisms [7–10]. Considering that the enormous number of intracellular metabolites within various biological matrices and concentrations range over several orders of magnitude (pM–mM), a variety of different analytical methods have been investigated [2,6]. Used analytical methods vary widely, from enzymatic assays [11] for targeted measurements to modern and more powerful soft ionization mass spectrometry (SIMS)2 coupled to liquid chromatography (LC) [12–19], gas chromatography
2 Abbreviations used: SIMS, soft ionization mass spectrometry; LC, liquid chromatography; GC, gas chromatography; CE, capillary electrophoresis; NMR, nuclear magnetic resonance; LC–MS, liquid chromatography combined with mass spectrometry; HILIC, hydrophilic interaction liquid chromatography; LC–ESI–MS/MS, liquid chromatography–electrospray ionization–tandem mass spectrometry; MRM, multiple reaction monitoring; ZIC–pHILIC, polymer-based zwitterionic hydrophilic interaction chromatography; LB, Luria–Bertani; MDL, method detection limit; ANOVA, analysis of variance; RPLC, reversed phase liquid chromatography; WWR, wrong-way-round.
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
(GC) [20–25], or capillary electrophoresis (CE) [26–28] as well as nuclear magnetic resonance (NMR) [29,30]. For comprehensive metabolite quantification studies, a combination of multiple analytical methods is usually indispensable, resulting in data sets with varying precision and accuracy [1,2,6,31]. However, for dynamic modeling of cellular metabolism, in connection with stimulus–response studies, reliable quantifications of absolute intracellular concentrations in defined cellular states and time points are of crucial importance [7,10,25]. Compared with GC, CE, and NMR, liquid chromatography combined with mass spectrometry (LC–MS) offers inherent advantages; LC–MS does not require further derivatization like GC, is able to separate polar and nonpolar compounds unlike CE with selectivity depending on column material, and offers a higher sensibility compared with NMR. Consequently, LC–MS approaches are becoming highly accepted as the most universal analytical platform for reliable quantification of intracellular metabolites [2,3,25]. Recently, hydrophilic interaction liquid chromatography (HILIC) has become increasingly popular because it enables the separation of charged and polar analytes and exhibits excellent compatibility with MS detection [32–34]. Related studies focused mainly on silica-based stationary HILIC phases, modified by various functional groups, enabling comprehensive quantitative analysis of metabolites under acidic and neutral mobile phase conditions [12,13,15,16]. Motivated by these encouraging studies, we developed the approach further for exploiting its potential, concentrating on so far nonaddressed modes of operation. In this study, we present an analytical platform based on two liquid chromatography–electrospray ionization–tandem mass spectrometry (LC–ESI–MS/MS) methods for quantitative profiling of more than 50 key metabolites of the cellular metabolism, comprising amino acids, sugar phosphates, organic acids, nucleotides, and coenzymes. Metabolites were detected by MS/MS on a triple quadrupole instrument in the multiple reaction monitoring (MRM) mode with high selectivity based on pre-optimized MRM transitions. Chromatographic separation was performed by polymer-based zwitterionic hydrophilic interaction chromatography (ZIC–pHILIC), and the mobile phase pH effect with respect to analyte sensibility was evaluated. Special focus placed on method quantification limits, linear dynamic range, and repeatability for a representative cross section of 56 common intracellular metabolites in acidic and alkaline mobile phase conditions. Finally, the applicability of the method was demonstrated by absolute quantification of selected metabolites in intracellular extracts of Escherichia coli showing consistency with previously published endogenous steady-state concentrations. In addition, occurrence of quantification bias caused by matrix effects is presented, comparing quantitative results obtained by standards-based external calibration, isotope dilution, and standard addition as calibration strategies. Materials and methods Chemicals Metabolite standards, reagents, and U-13C-labeled lyophilized algal cells were supplied by Sigma–Aldrich (Taufkirchen, Germany). MS-grade water was purchased from VWR (Darmstadt, Germany). MS-grade acetonitrile was purchased from Carl Roth (Essen, Germany). Standard stock solutions were prepared in LC–MS water and stored at 70 °C.
5
(NH4)2HPO4, 20.3 mM (NH4)2SO4, 6.2 mM Na2SO4, 1.0 mM MgSO4, 0.1 mM CaCl2, and 0.01 mM thiamine hydrochloride. Overnight precultures (12 h, 37 °C, 120 rpm) were inoculated from a cryoculture (50% glycerol and Luria–Bertani [LB] medium) and were grown in 100-ml baffled shake flasks with 20 ml of minimal medium (12 h, 37 °C, 120 rpm). Main cultures (37 °C, 150 rpm) were inoculated in 1:100 dilution and grown in triplicates in 500-ml baffled shake flasks with 60 ml of minimal medium (10 h, 37 °C, 150 rpm). Sampling, quenching, and metabolite extraction Cells were sampled during the exponential growth phase at a biomass concentration of approximately 2 g L 1. Amounts of approximately 4 mg biomass were sampled by fast centrifugation (20,000g, 20 s) and washed with 2 ml of isotonic 0.9% (w/v) sodium chloride solution (20,000g, 20 s). Biomasses were quenched by liquid nitrogen ( 196 °C) and temporarily stored at 70 °C. Subsequently, a defined amount of water (MS grade) was added to obtain an extraction concentration of 15 g L 1. Resulting suspensions were immediately pre-incubated for 1 min at 100 °C in a water bath for enzymatic inactivation and resuspended by short-time vortexing. Subsequently, samples were incubated for 5 min at 100 °C in a water bath and afterward chilled on ice water. Metabolite extracts were separated from cell debris by centrifugation (20,000g, 10 min) and stored at 70 °C [35]. Preparation of U-13C-labeled internal standard Commercially available U-13C-labeled lyophilized algal cells (>99 atom% 13C, lot no. 487945, Sigma–Aldrich) were weighed into 2-ml reaction vessels in small amounts. Preheated water (100 °C) was added aiming for an extraction concentration of 90 g L 1. The resulting suspensions were incubated at 100 °C for 2 min in a water bath and resuspended by short-time vortexing. The procedure was repeated two more times, and the resulting samples were chilled on ice water. 13C-labeled metabolite extracts were separated from algal cells by centrifugation (20,000g, 10 min) and stored at 70 °C [25]. Optimization of chromatographic and source conditions Development of chromatography was performed on a Sequant ZIC–pHILIC column (150 2.1 mm, 5 lm) with guard column (Sequant ZIC–pHILIC, 20 2.1 mm, 5 lm). The optimization of chromatographic conditions was based on a selection of highly polar intracellular metabolites within the central metabolism (6-phophogluconate, fructose 1,6-bisphosphate, glucose 6-phosphate, phosphoenolpyruvate, malate, succinate, and citrate). Bicratic chromatographic runs were evaluated in regard to selectivity and sensitivity of applied standard mixtures. Optimization was focused on pH of mobile phases (5.0, 7.0, and 9.0), buffer concentration (2–20 mM ammonium acetate), flow rate (0.05–0.20 ml/min), column temperature (20–60 °C), and gradient slope of polar eluent B (2,5–5% B/min) within a bicratic elution mode. In addition, source parameters (ESI) were optimized with respect to maximal signal intensities, focusing on the nitrogen gas flow (7–13 L/min), nebulizer pressure (15–60 psi), and capillary voltage (±3000–4000 V). An overview of optimized chromatographic and source parameters is shown further below.
Strains and growth condition Sample processing and optimized chromatographic conditions E. coli K-12 MG1655 (purchased from DSMZ, Leipzig, Germany) was grown in minimal medium with 0.75% (w/v) D-glucose as sole carbon source. The minimal medium contained the following components: 63.1 mM NaH2PO4, 26.9 mM K2HPO4, 10.0 mM
Standards and endogenous cellular extracts were analyzed on an Agilent 1200 high-performance liquid chromatography (HPLC) system (Agilent, Waldbronn, Germany) consisting of a degasser, a
6
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
binary pump, a thermostated column compartment, and a bio-inert thermostated autosampler (Agilent 1260 Infinity series) maintained at 5 °C. Samples and standards were prepared in 60% (v/v) acetonitrile and 10 mM ammonium acetate buffer (adjusted to pH 5.6 or 9.2). An injection of 5 ll was separated on a Sequant ZIC–pHILIC column (150 2.1 mm, 5 lm) with guard column (Sequant ZIC–pHILIC, 20 2.1 mm, 5 lm) kept at 40 °C and with a flow rate of 0.2 ml/min. Mobile phases were composed of 10% (v/v) aqueous buffer solution (10 mM ammonium acetate) and 90% (v/v) acetonitrile for eluent A and 90% (v/v) aqueous buffer solution and 10% (v/v) acetonitrile for eluent B, both adjusted to pH 5.6 (with acetic acid) or pH 9.2 (with 25% [v/v] ammonium hydroxide). Gradient elution was carried out by the following program: isocratic hold 0% B for 1 min, next a linear gradient from 0% B to 75% B for 30 min, followed by a linear gradient from 75% B to 100% B for 4 min. At the end of the run, the column was washed with 100% B for 5 min and afterward equilibrated to starting conditions by a linear gradient from 100% B to 0% B for 10 min and 0% B for 15 min. Instrumentation and ESI–MS/MS Data were acquired with an Agilent 6410B Triple Quad mass spectrometer (Agilent) with ESI ion source. Negative (ESI ) and positive (ESI+) ionization modes were performed in separated runs in MRM mode depending on metabolite group. Source conditions were as follows: nitrogen gas flow rate of 10 L min 1 with a temperature of 350 °C, a nebulizer pressure of 30 psi, and a capillary voltage of 4000 V. Fragmentor voltages and collision energies were optimized in continuous direct injection mode for each metabolite precursor-to-product ion transition. Dwell time was adjusted to 100 ms for each MRM transition. System control, acquisition, and analysis of data were performed by use of commercial MassHunter B.04.00 software. Linearity range and sensibility tests Linearity ranges of both method configurations (acidic and alkaline mobile phase conditions) were evaluated with regard to the amount of respective injected analytes. Multicomponent standard solutions were prepared from properly adjusted metabolite mixtures and were assigned to three different groups depending on ionization mode and retention time. Focusing on small concentrations (65 lM) and depending on sensitivities of measured metabolite groups, a calibration range of 1 nM to 5 lM based on 12 levels was applied. Additional linearities of several selected metabolites were evaluated by use of an extended calibration range (0.5– 200 lM) based on 12 levels. Sensibility tests were implemented by measuring 5-lM standards for each metabolite under acidic and alkaline mobile phase conditions. Absolute quantification in bacterial metabolite extracts (quantification strategies) The suitability of the method was demonstrated by quantification of absolute concentrations of selected metabolites in intracellular extracts of wild-type E. coli biomasses. In addition, quantification bias caused by matrix effects was evaluated by comparison of quantitative results by external calibration, isotope dilution, and standard addition. Multicomponent standard solutions for external calibration of defined nonlabeled metabolite mixtures were obtained. Depending on previously estimated concentrations (obtained by spiking experiments), metabolite-specific calibration ranges with seven levels were adjusted. Quantifications based on isotope dilution were performed by constant addition of U-13C-labeled algal extracts as internal standard to samples and
external calibration standards. External calibration curves were prepared by linear regression of ratios of analyte peak areas and internal standard areas against the concentration levels of the respective analyte. Standard addition was performed by addition of defined amounts of respective standard directly to sample aliquots. Applied multicomponent metabolite mixtures were adjusted based on previously estimated concentrations of selected metabolites. Internal calibration curves were prepared by linear regression of nonspiked and incrementally increased analyte peak areas against concentration levels of the respective added standards. Results The objective of the current work was to develop a robust LC–MS/MS method for quantitative analysis of a broad range of analytes with high selectivity and sensitivity in complex matrix samples such as intracellular metabolite extracts. To overcome matrix effects and coelution of isomeric compounds, frequently present in intracellular matrices, a suitable chromatographic method is required. Recent developments in polymer-based HILIC stationary phases enable chromatographic separation of complex samples under alkaline conditions, improving sensitivity and chromatographic performance. Optimization of MS/MS parameters Measurements were performed on an HPLC–ESI–QQQ (triple quadrupole) instrument in MRM mode with a metal-free bio-inert HPLC system. Precursor-to-product ion transitions were optimized using continuous direct injection of pure standard solutions in positive and negative ionization modes. Suitable product ions, fragmentor voltages, and collision energies were selected in regard to maximal signal intensities. Specific MRM transitions for a representative cross section of 56 common intracellular metabolites, including amino acids, carboxylic acids, sugar phosphates, nucleotides, and coenzymes as well as their U-13C-labeled analogues for isotope dilution, are provided in the online Supplementary material. Chromatographic performance under alkaline and acidic mobile phases Chromatographic performance was evaluated on retention time reproducibility, peak width, and peak shape. To avoid potential column effects due to pH switches, two different columns from the same sorbent lot were used for acidic and alkaline conditions. Even when similar chromatographic quality was obtained for organic acids and amino acids, the peak shape was significantly improved for nucleotides under basic conditions. Only citrate and isocitrate presented better resolution in higher concentration levels under acidic conditions. A selection of extracted ion chromatograms of isobaric metabolites is shown in Fig. 1. Retention time reproducibility was fairly high: Analyzing individual series, 98% of the 5-lM standards showed absolute standard deviations less than 0.03 min under alkaline conditions, whereas 88% of the metabolites fulfill the likewise criterion under acidic conditions. The highest variability was observed for cytidine triphosphate (CTP, 0.17 min) and guanosine triphosphate (GTP, 0.12 min) under acidic conditions. The comparison of multiple measurement series confirmed this finding, with 98% of the metabolites showing absolute time shifts less than 0.04 min under alkaline conditions and 80% under acidic conditions. An overview of analyte retention times is given in Table 1. A compilation of MRM chromatograms of all 56 metabolites measured under acidic and alkaline mobile phase conditions is included in the Supplementary material.
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
7
Fig.1. Extracted ion chromatograms obtained with optimized HILIC method under alkaline (pH 9.2) and acidic (pH 5.6) mobile phase conditions focusing on isobaric compounds frequently present in intracellular matrices. Citrate/Isocitrate mixtures of 5 and 50 lM were injected at pHs 9.2 and 5.6, respectively.
Evaluation of linearity and method detection limits for different mobile phase pHs The linearity range for each metabolite was determined under acidic (pH 5.6) and alkaline (pH 9.2) mobile phase conditions. A range of metabolite concentrations was considered linear when the squared correlation coefficient (R2) was better than 0.98 for the average of three calibration curves. To increase the predictive significance of the measurements, internal standards were considered for monitoring of instrumental fluctuations but not for normalization of obtained peak areas. Extending the calibration range from 0.5 to 200 lM, the linearity of instrument response with respect to the amount of the injected analyte was revealed to be specific for each metabolite, sample context, and mobile phase condition. Fig. 2 shows two different metabolite standards with exemplary deviations from linearity of applied calibration lines. Obtained metabolic-specific linearity ranges differ significantly at lower concentrations (1 nM to 5 lM) between both conditions. The median for the lower linearity limits was approximately 30 nM for alkaline and 120 nM for acidic mobile phase conditions. Within alkaline conditions, 40% of the metabolites show a linearity range between 5 and 20 nM as lower boundary and 5 lM as upper limit. In contrast to this, only 12% of the metabolites fulfilled the same criterion under acidic conditions. An overview of metabolic-specific linearity ranges for acidic and alkaline mobile phase conditions is shown in Table 1. The method detection limit (MDL) was calculated as the amount of metabolite that results at statistically significant peak areas above background noise. MDL was based on the standard deviation of triplicates at the lower linearity boundary (lowest concentration of the linearity range) multiplied by the expansion coefficient (t = 6.96) for defining the 99% confidence level [36]. In column MDLs were 0.01 to 3.7 pmol for alkaline conditions and 0.02 to 14.7 pmol for acidic conditions. Fully 89% of the analyzed metabolites showed MDLs lower than 100 nM (0.5 pmol on column) under alkaline conditions, whereas just 46% of the metabolites fulfilled this criterion under acidic conditions (Fig. 3). The lowest detection limit was achieved for malic acid in positive ionization mode (2.23 nM, 0.011 pmol on column) using basic conditions. An overview of the metabolite-specific MDLs for acidic and alkaline mobile phase conditions is shown in Table 1. Comparing both conditions, the most significant differences were observed for the di- and triphosphate nucleotides and NADP as well as for
fumaric acid, isocitrate, aspartic acid, and glutamine, showing approximately 20 to 50 times lower detection limits under alkaline conditions. ESI–MS responsivity under alkaline and acidic mobile phases Next, the impact of pH of mobile phases on ESI–MS responsivity was studied based on comparison of peak areas obtained by injection of 5-lM metabolite standards under both conditions. Special care was given to the preparation of 10-mM ammonium acetate buffer solutions to afford comparable ion strength under acidic and alkaline conditions. Either acetic acid or ammonium hydroxide was added in equimolar amounts to adjust the buffer system pH to 5.6 or 9.2. To compare both scenarios, averaged peak areas obtained by alkaline conditions were normalized by those obtained by acidic conditions. Strikingly, 70% of the metabolites presented more than 1.5-fold higher signal intensities under alkaline mobile phase conditions (Fig. 4). Absolute quantification in bacterial metabolite extracts For challenging the applicability of the approach, absolute concentrations of selected metabolites were quantified in intracellular extracts of E. coli biomasses. Because matrix effects are suspected to be a major source for inaccuracies and unreliability in ESI– MS-based quantitative metabolomics [6,16,18], related impacts were evaluated as well. Quantification strategies were adjusted based on previously estimated concentrations of the targeted intracellular metabolites. Intracellular pool sizes were detected by standard-based external calibration, isotope dilution, and standard addition for subsequent comparison. Standard-based external calibration and isotope dilution were applied by metabolic-specific calibration ranges in multicomponent mixtures. Standard addition was performed for each metabolite at concentrations corresponding to approximately 2-, 3-, and 4-fold the estimated values. Results obtained from exemplary measurements of six intracellular metabolites with three different approaches are summarized in Fig. 5. Corresponding correlation coefficients are shown in Table 2. Good (R2 P 0.99) to excellent (R2 P 0.999) linearity coefficients were found for all three calibration strategies. Four replicates were made for sample analysis, where calculated relative standard deviations of the measured concentrations were less than 9% for each metabolite. As depicted in Fig. 5, individual sample analysis
Metabolite
Polarity mode
Linearity range
Method detection limit
Alkaline
Acid
Alkaline
Acid
(min)
(min)
4.2 606.0 138.0 8.5 59.6 513.3 454.1 41.6 181.9 1131.1 1713.6 67.5 67.5 53.2 79.3 2949.3 423.1 532.7 148.4 148.4 13.1 1495.0 144.1 1428.0 1130.4 198.5 2274.7 1525.2 160.8 1030.1 2004.1 123.4 1094.0 1711.8 375.8 3.6 18.2 22.0
0.10 0.17 0.59 0.01 0.05 2.20 0.11 0.02 0.15 3.72 – 0.77 0.77 0.18 0.17 1.80 0.27 0.43 0.44 0.44 0.04 0.16 0.07 0.20 0.30 0.25 0.35 0.51 0.25 0.28 0.26 0.02 0.20 0.36 0.30 0.21 0.06 0.05
0.02 3.03 0.69 0.04 0.30 2.57 2.27 0.21 0.91 5.66 8.57 0.34 0.34 0.27 0.40 14.75 2.12 2.66 0.74 0.74 0.07 7.47 0.72 7.14 5.65 0.99 11.37 7.63 0.80 5.15 10.02 0.62 5.47 8.56 1.88 0.02 0.09 0.11
8.02 ± 0.04 22.07 ± 0.02 20.59 ± 0.02 21.83 ± 0.01 21.42 ± 0.02 24.15 ± 0.01 24.83 ± 0.01 23.08 ± 0.01 21.15 ± 0.01 24.03 ± 0.01 – 21.11 ± 0.02 21.27 ± 0.01 21.63 ± 0.02 22.49 ± 0.02 24.99 ± 0.01 20.77 ± 0.01 23.61 ± 0.02 22.86 ± 0.02 23.11 ± 0.02 19.64 ± 0.01 23.15 ± 0.01 19.41 ± 0.01 21.31 ± 0.01 22.80 ± 0.02 22.07 ± 0.01 23.95 ± 0.02 25.45 ± 0.01 21.38 ± 0.01 23.21 ± 0.01 24.66 ± 0.01 20.54 ± 0.01 22.57 ± 0.01 24.11 ± 0.01 15.35 ± 0.01 15.86 ± 0.01 17.64 ± 0.01 17.51 ± 0.01
8.63 ± 0.07 22.57 ± 0.01 20.12 ± 0.01 22.14 ± 0.01 21.77 ± 0.01 23.97 ± 0.02 24.75 ± 0.02 22.76 ± 0.01 21.14 ± 0.01 24.53 ± 0.01 21.73 ± 0.11 – – 21.81 ± 0.01 22.72 ± 0.01 25.93 ± 0.07 20.78 ± 0.01 24.07 ± 0.02 23.32c ± 0.02 23.61c ± 0.02 20.03 ± 0.01 23.51 ± 0.01 19.33 ± 0.01 22.08 ± 0.03 23.75 ± 0.05 22.11 ± 0.01 24.91 ± 0.03 26.73c ± 0.12 21.34 ± 0.01 24.02 ± 0.01 25.90 ± 0.17 20.40 ± 0.01 23.33 ± 0.02 25.24 ± 0.01 15.65 ± 0.01 16.26 ± 0.01 18.02 ± 0.01 17.92 ± 0.01
33.3 14.0 4.9 28.5 41.0
286.0 38.3 22.6 99.2 167.5
0.17 0.07 0.02 0.14 0.20
1.43 0.19 0.11 0.50 0.84
15.76 ± 0.01 18.69 ± 0.01 21.45 ± 0.01 17.96 ± 0.01 16.91 ± 0.01
15.00 ± 0.01 19.08 ± 0.01 22.10 ± 0.01 18.37 ± 0.01 17.32 ± 0.01
0.9964
9.9
140.9
0.05
0.70
18.45 ± 0.01
18.90 ± 0.01
0.9968
9.1
66.8
0.05
0.33
14.37 ± 0.01
14.76 ± 0.01
5000
0.9964
26.4
49.7
0.13
0.25
14.54 ± 0.01
14.94 ± 0.01
5000
0.9967
43.1
97.6
0.22
0.49
17.14 ± 0.01
17.58 ± 0.02
20
5000
0.9965
7.9
35.9
0.04
0.18
12.74 ± 0.01
13.15 ± 0.01
0.9990
20
5000
0.9966
17.2
22.2
0.09
0.11
13.22 ± 0.01
13.64 ± 0.01
0.9991
100
5000
0.9964
15.0
100.6
0.08
0.50
18.35 ± 0.01
18.79 ± 0.01
R2
Min (nM)
Max (nM)
R2
(nM)
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (+) (+) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (+) (+) (+)
50 50 500 20 50 50 20 20 50 1000 – 20 20 50 50 500 200 50 100 100 10 20 50 100 100 100 100 500 100 50 200 100 50 100 100 20 10 20
5000 5000 5000 5000 5000 5000 5000 1000 5000 5000 – 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 2000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000
1.0000 0.9999 0.9999 0.9998 0.9974 0.9938 0.9953 0.9899 0.9884 0.9903 – 0.9973 0.9973 0.9954 0.9972 0.9998 0.9990 0.9904 0.9874 0.9874 0.9982 0.9901 0.9995 0.9928 0.9925 0.9990 0.9844 0.9915 0.9990 0.9905 0.9915 0.9991 0.9858 0.9968 0.9970 0.9996 0.9994 0.9990
50 500 500 20 100 1000 500 50 200 1000 2000 50 50 100 100 2000 500 500 500 500 20 500 100 1000 1000 100 1000 2000 100 1000 1000 100 500 2000 100 50 20 50
5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000
0.9998 0.9951 0.9976 0.9998 0.9916 0.9838 0.9945 0.9854 0.9953 0.9950 0.9972 0.9956 0.9956 0.9990 0.9992 0.9980 0.9985 0.9994 0.9984 0.9984 0.9963 0.9847 0.9987 0.9832 0.9836 0.9960 0.9975 0.9971 0.9974 0.9938 0.9842 0.9985 0.9840 0.9850 0.9988 0.9938 0.9920 0.9963
20.7 33.8 117.1 2.2 9.4 439.6 21.5 3.9 30.1 744.3 – 154.5 154.5 35.1 33.5 360.8 53.9 86.5 87.8 87.8 7.6 31.4 13.4 40.1 59.6 49.9 70.9 102.4 49.9 56.3 51.9 4.7 39.2 72.8 59.3 41.7 12.7 10.2
L-Ala
(+) (+) (+) (+) (+)
20 10 10 200 100
5000 5000 5000 5000 5000
0.9996 0.9991 0.9953 0.9992 0.9990
200 50 50 500 200
5000 5000 5000 5000 5000
0.9986 0.9956 0.9938 0.9963 0.9956
L-Ser
(+)
100
5000
0.9988
200
5000
L-Pro
(+)
10
5000
0.9992
50
5000
L-Val
(+)
20
5000
0.9991
100
L-Thr
(+)
50
5000
0.9993
100
L-Leu
(+)
20
5000
0.9993
L-Ile
(+)
20
5000
L-Asn
(+)
50
5000
L-Hom-Ser
O-Ac-Ser OSHS lL-Cysta Gly
Alkaline
Acid
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
Alkaline
Max (nM)
aKB Fum Suc Mal aKG Cit IsoCit cis-Aco Xu5P 6PG E4Pb R5Pa Ru5Pa F6P G6P FbP DHAP PEP 2PGa 3PGa NAD NADP AMP ADP ATP GMP GDP GTP CMP CDP CTP UMP UDP UTP cAMP AIBA GABA
Retention time
(pmol, on column)
Min (nM)
Acid
8
Table 1 Overview of obtained linearity ranges (
Contents lists available at ScienceDirect
Analytical Biochemistry journal homepage: www.elsevier.com/locate/yabio
Alkaline conditions in hydrophilic interaction liquid chromatography for intracellular metabolite quantification using tandem mass spectrometry Attila Teleki 1, Andrés Sánchez-Kopper 1, Ralf Takors ⇑ Institute of Biochemical Engineering, University of Stuttgart, 70569 Stuttgart, Germany
a r t i c l e
i n f o
Article history: Received 17 September 2014 Received in revised form 8 December 2014 Accepted 5 January 2015 Available online 16 January 2015 Keywords: Metabolic profiling HILIC LC–MS/MS
a b s t r a c t Modeling of metabolic networks as part of systems metabolic engineering requires reliable quantitative experimental data of intracellular concentrations. The hydrophilic interaction liquid chromatography–e lectrospray ionization–tandem mass spectrometry (HILIC–ESI–MS/MS) method was used for quantitative profiling of more than 50 hydrophilic key metabolites of cellular metabolism. Without prior derivatization, sugar phosphates, organic acids, nucleotides, and amino acids were measured under alkaline and acidic mobile phase conditions with pre-optimized multiple reaction monitoring (MRM) transitions. Irrespective of the polarity mode of the acquisition method used, alkaline conditions achieved the best quantification limits and linear dynamic ranges. Fully 90% of the analyzed metabolites presented detection limits better than 0.5 pmol (on column), and 70% presented 1.5-fold higher signal intensities under alkaline mobile phase conditions. The quality of the method was further demonstrated by absolute quantification of selected metabolites in intracellular extracts of Escherichia coli. In addition, quantification bias caused by matrix effects was investigated by comparison of calibration strategies: standard-based external calibration, isotope dilution, and standard addition with internal standards. Here, we recommend the use of alkaline mobile phase with polymer-based zwitterionic hydrophilic interaction chromatography (ZIC–pHILIC) as the most sensitive scenario for absolute quantification for a broad range of metabolites. Ó 2015 Elsevier Inc. All rights reserved.
Metabolomics addresses the identification and quantification of small molecule metabolites, in essence the reactants of biological systems, and enables the comprehensive analysis of complex biochemical networks systemic studies. Whereas the genome and proteome represent upstream biochemical events, the metabolome reveals underlying regulatory mechanisms and therefore reflects, most closely, the actual cellular physiological state [1–3]. As such, metabolome analysis can be subdivided into targeted and nontargeted approaches. Nontargeted high-throughput applications such as metabolic fingerprinting aim at the rapid and qualitative or semiquantitative analysis of whole-cell metabolic patterns, reducing the analytical effort considerably [4–6]. On the contrary, targeted applications such as quantitative metabolic profiling are focused on the nonbiased absolute quantification of a subset of metabolites of predefined metabolic pathways or classes of compounds [2,4,6]. In this context, monitoring of intracellular ⇑ Corresponding author. Fax: +49 0711 685 65164. 1
E-mail address: [email protected] (R. Takors). Contributed equally to this manuscript.
http://dx.doi.org/10.1016/j.ab.2015.01.002 0003-2697/Ó 2015 Elsevier Inc. All rights reserved.
metabolism dynamics as a result of rapid changes in the extracellular environment within stimulus–response experiments can be used to investigate in vivo enzyme kinetics offering insights into underlying regulatory mechanisms [7–10]. Considering that the enormous number of intracellular metabolites within various biological matrices and concentrations range over several orders of magnitude (pM–mM), a variety of different analytical methods have been investigated [2,6]. Used analytical methods vary widely, from enzymatic assays [11] for targeted measurements to modern and more powerful soft ionization mass spectrometry (SIMS)2 coupled to liquid chromatography (LC) [12–19], gas chromatography
2 Abbreviations used: SIMS, soft ionization mass spectrometry; LC, liquid chromatography; GC, gas chromatography; CE, capillary electrophoresis; NMR, nuclear magnetic resonance; LC–MS, liquid chromatography combined with mass spectrometry; HILIC, hydrophilic interaction liquid chromatography; LC–ESI–MS/MS, liquid chromatography–electrospray ionization–tandem mass spectrometry; MRM, multiple reaction monitoring; ZIC–pHILIC, polymer-based zwitterionic hydrophilic interaction chromatography; LB, Luria–Bertani; MDL, method detection limit; ANOVA, analysis of variance; RPLC, reversed phase liquid chromatography; WWR, wrong-way-round.
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
(GC) [20–25], or capillary electrophoresis (CE) [26–28] as well as nuclear magnetic resonance (NMR) [29,30]. For comprehensive metabolite quantification studies, a combination of multiple analytical methods is usually indispensable, resulting in data sets with varying precision and accuracy [1,2,6,31]. However, for dynamic modeling of cellular metabolism, in connection with stimulus–response studies, reliable quantifications of absolute intracellular concentrations in defined cellular states and time points are of crucial importance [7,10,25]. Compared with GC, CE, and NMR, liquid chromatography combined with mass spectrometry (LC–MS) offers inherent advantages; LC–MS does not require further derivatization like GC, is able to separate polar and nonpolar compounds unlike CE with selectivity depending on column material, and offers a higher sensibility compared with NMR. Consequently, LC–MS approaches are becoming highly accepted as the most universal analytical platform for reliable quantification of intracellular metabolites [2,3,25]. Recently, hydrophilic interaction liquid chromatography (HILIC) has become increasingly popular because it enables the separation of charged and polar analytes and exhibits excellent compatibility with MS detection [32–34]. Related studies focused mainly on silica-based stationary HILIC phases, modified by various functional groups, enabling comprehensive quantitative analysis of metabolites under acidic and neutral mobile phase conditions [12,13,15,16]. Motivated by these encouraging studies, we developed the approach further for exploiting its potential, concentrating on so far nonaddressed modes of operation. In this study, we present an analytical platform based on two liquid chromatography–electrospray ionization–tandem mass spectrometry (LC–ESI–MS/MS) methods for quantitative profiling of more than 50 key metabolites of the cellular metabolism, comprising amino acids, sugar phosphates, organic acids, nucleotides, and coenzymes. Metabolites were detected by MS/MS on a triple quadrupole instrument in the multiple reaction monitoring (MRM) mode with high selectivity based on pre-optimized MRM transitions. Chromatographic separation was performed by polymer-based zwitterionic hydrophilic interaction chromatography (ZIC–pHILIC), and the mobile phase pH effect with respect to analyte sensibility was evaluated. Special focus placed on method quantification limits, linear dynamic range, and repeatability for a representative cross section of 56 common intracellular metabolites in acidic and alkaline mobile phase conditions. Finally, the applicability of the method was demonstrated by absolute quantification of selected metabolites in intracellular extracts of Escherichia coli showing consistency with previously published endogenous steady-state concentrations. In addition, occurrence of quantification bias caused by matrix effects is presented, comparing quantitative results obtained by standards-based external calibration, isotope dilution, and standard addition as calibration strategies. Materials and methods Chemicals Metabolite standards, reagents, and U-13C-labeled lyophilized algal cells were supplied by Sigma–Aldrich (Taufkirchen, Germany). MS-grade water was purchased from VWR (Darmstadt, Germany). MS-grade acetonitrile was purchased from Carl Roth (Essen, Germany). Standard stock solutions were prepared in LC–MS water and stored at 70 °C.
5
(NH4)2HPO4, 20.3 mM (NH4)2SO4, 6.2 mM Na2SO4, 1.0 mM MgSO4, 0.1 mM CaCl2, and 0.01 mM thiamine hydrochloride. Overnight precultures (12 h, 37 °C, 120 rpm) were inoculated from a cryoculture (50% glycerol and Luria–Bertani [LB] medium) and were grown in 100-ml baffled shake flasks with 20 ml of minimal medium (12 h, 37 °C, 120 rpm). Main cultures (37 °C, 150 rpm) were inoculated in 1:100 dilution and grown in triplicates in 500-ml baffled shake flasks with 60 ml of minimal medium (10 h, 37 °C, 150 rpm). Sampling, quenching, and metabolite extraction Cells were sampled during the exponential growth phase at a biomass concentration of approximately 2 g L 1. Amounts of approximately 4 mg biomass were sampled by fast centrifugation (20,000g, 20 s) and washed with 2 ml of isotonic 0.9% (w/v) sodium chloride solution (20,000g, 20 s). Biomasses were quenched by liquid nitrogen ( 196 °C) and temporarily stored at 70 °C. Subsequently, a defined amount of water (MS grade) was added to obtain an extraction concentration of 15 g L 1. Resulting suspensions were immediately pre-incubated for 1 min at 100 °C in a water bath for enzymatic inactivation and resuspended by short-time vortexing. Subsequently, samples were incubated for 5 min at 100 °C in a water bath and afterward chilled on ice water. Metabolite extracts were separated from cell debris by centrifugation (20,000g, 10 min) and stored at 70 °C [35]. Preparation of U-13C-labeled internal standard Commercially available U-13C-labeled lyophilized algal cells (>99 atom% 13C, lot no. 487945, Sigma–Aldrich) were weighed into 2-ml reaction vessels in small amounts. Preheated water (100 °C) was added aiming for an extraction concentration of 90 g L 1. The resulting suspensions were incubated at 100 °C for 2 min in a water bath and resuspended by short-time vortexing. The procedure was repeated two more times, and the resulting samples were chilled on ice water. 13C-labeled metabolite extracts were separated from algal cells by centrifugation (20,000g, 10 min) and stored at 70 °C [25]. Optimization of chromatographic and source conditions Development of chromatography was performed on a Sequant ZIC–pHILIC column (150 2.1 mm, 5 lm) with guard column (Sequant ZIC–pHILIC, 20 2.1 mm, 5 lm). The optimization of chromatographic conditions was based on a selection of highly polar intracellular metabolites within the central metabolism (6-phophogluconate, fructose 1,6-bisphosphate, glucose 6-phosphate, phosphoenolpyruvate, malate, succinate, and citrate). Bicratic chromatographic runs were evaluated in regard to selectivity and sensitivity of applied standard mixtures. Optimization was focused on pH of mobile phases (5.0, 7.0, and 9.0), buffer concentration (2–20 mM ammonium acetate), flow rate (0.05–0.20 ml/min), column temperature (20–60 °C), and gradient slope of polar eluent B (2,5–5% B/min) within a bicratic elution mode. In addition, source parameters (ESI) were optimized with respect to maximal signal intensities, focusing on the nitrogen gas flow (7–13 L/min), nebulizer pressure (15–60 psi), and capillary voltage (±3000–4000 V). An overview of optimized chromatographic and source parameters is shown further below.
Strains and growth condition Sample processing and optimized chromatographic conditions E. coli K-12 MG1655 (purchased from DSMZ, Leipzig, Germany) was grown in minimal medium with 0.75% (w/v) D-glucose as sole carbon source. The minimal medium contained the following components: 63.1 mM NaH2PO4, 26.9 mM K2HPO4, 10.0 mM
Standards and endogenous cellular extracts were analyzed on an Agilent 1200 high-performance liquid chromatography (HPLC) system (Agilent, Waldbronn, Germany) consisting of a degasser, a
6
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
binary pump, a thermostated column compartment, and a bio-inert thermostated autosampler (Agilent 1260 Infinity series) maintained at 5 °C. Samples and standards were prepared in 60% (v/v) acetonitrile and 10 mM ammonium acetate buffer (adjusted to pH 5.6 or 9.2). An injection of 5 ll was separated on a Sequant ZIC–pHILIC column (150 2.1 mm, 5 lm) with guard column (Sequant ZIC–pHILIC, 20 2.1 mm, 5 lm) kept at 40 °C and with a flow rate of 0.2 ml/min. Mobile phases were composed of 10% (v/v) aqueous buffer solution (10 mM ammonium acetate) and 90% (v/v) acetonitrile for eluent A and 90% (v/v) aqueous buffer solution and 10% (v/v) acetonitrile for eluent B, both adjusted to pH 5.6 (with acetic acid) or pH 9.2 (with 25% [v/v] ammonium hydroxide). Gradient elution was carried out by the following program: isocratic hold 0% B for 1 min, next a linear gradient from 0% B to 75% B for 30 min, followed by a linear gradient from 75% B to 100% B for 4 min. At the end of the run, the column was washed with 100% B for 5 min and afterward equilibrated to starting conditions by a linear gradient from 100% B to 0% B for 10 min and 0% B for 15 min. Instrumentation and ESI–MS/MS Data were acquired with an Agilent 6410B Triple Quad mass spectrometer (Agilent) with ESI ion source. Negative (ESI ) and positive (ESI+) ionization modes were performed in separated runs in MRM mode depending on metabolite group. Source conditions were as follows: nitrogen gas flow rate of 10 L min 1 with a temperature of 350 °C, a nebulizer pressure of 30 psi, and a capillary voltage of 4000 V. Fragmentor voltages and collision energies were optimized in continuous direct injection mode for each metabolite precursor-to-product ion transition. Dwell time was adjusted to 100 ms for each MRM transition. System control, acquisition, and analysis of data were performed by use of commercial MassHunter B.04.00 software. Linearity range and sensibility tests Linearity ranges of both method configurations (acidic and alkaline mobile phase conditions) were evaluated with regard to the amount of respective injected analytes. Multicomponent standard solutions were prepared from properly adjusted metabolite mixtures and were assigned to three different groups depending on ionization mode and retention time. Focusing on small concentrations (65 lM) and depending on sensitivities of measured metabolite groups, a calibration range of 1 nM to 5 lM based on 12 levels was applied. Additional linearities of several selected metabolites were evaluated by use of an extended calibration range (0.5– 200 lM) based on 12 levels. Sensibility tests were implemented by measuring 5-lM standards for each metabolite under acidic and alkaline mobile phase conditions. Absolute quantification in bacterial metabolite extracts (quantification strategies) The suitability of the method was demonstrated by quantification of absolute concentrations of selected metabolites in intracellular extracts of wild-type E. coli biomasses. In addition, quantification bias caused by matrix effects was evaluated by comparison of quantitative results by external calibration, isotope dilution, and standard addition. Multicomponent standard solutions for external calibration of defined nonlabeled metabolite mixtures were obtained. Depending on previously estimated concentrations (obtained by spiking experiments), metabolite-specific calibration ranges with seven levels were adjusted. Quantifications based on isotope dilution were performed by constant addition of U-13C-labeled algal extracts as internal standard to samples and
external calibration standards. External calibration curves were prepared by linear regression of ratios of analyte peak areas and internal standard areas against the concentration levels of the respective analyte. Standard addition was performed by addition of defined amounts of respective standard directly to sample aliquots. Applied multicomponent metabolite mixtures were adjusted based on previously estimated concentrations of selected metabolites. Internal calibration curves were prepared by linear regression of nonspiked and incrementally increased analyte peak areas against concentration levels of the respective added standards. Results The objective of the current work was to develop a robust LC–MS/MS method for quantitative analysis of a broad range of analytes with high selectivity and sensitivity in complex matrix samples such as intracellular metabolite extracts. To overcome matrix effects and coelution of isomeric compounds, frequently present in intracellular matrices, a suitable chromatographic method is required. Recent developments in polymer-based HILIC stationary phases enable chromatographic separation of complex samples under alkaline conditions, improving sensitivity and chromatographic performance. Optimization of MS/MS parameters Measurements were performed on an HPLC–ESI–QQQ (triple quadrupole) instrument in MRM mode with a metal-free bio-inert HPLC system. Precursor-to-product ion transitions were optimized using continuous direct injection of pure standard solutions in positive and negative ionization modes. Suitable product ions, fragmentor voltages, and collision energies were selected in regard to maximal signal intensities. Specific MRM transitions for a representative cross section of 56 common intracellular metabolites, including amino acids, carboxylic acids, sugar phosphates, nucleotides, and coenzymes as well as their U-13C-labeled analogues for isotope dilution, are provided in the online Supplementary material. Chromatographic performance under alkaline and acidic mobile phases Chromatographic performance was evaluated on retention time reproducibility, peak width, and peak shape. To avoid potential column effects due to pH switches, two different columns from the same sorbent lot were used for acidic and alkaline conditions. Even when similar chromatographic quality was obtained for organic acids and amino acids, the peak shape was significantly improved for nucleotides under basic conditions. Only citrate and isocitrate presented better resolution in higher concentration levels under acidic conditions. A selection of extracted ion chromatograms of isobaric metabolites is shown in Fig. 1. Retention time reproducibility was fairly high: Analyzing individual series, 98% of the 5-lM standards showed absolute standard deviations less than 0.03 min under alkaline conditions, whereas 88% of the metabolites fulfill the likewise criterion under acidic conditions. The highest variability was observed for cytidine triphosphate (CTP, 0.17 min) and guanosine triphosphate (GTP, 0.12 min) under acidic conditions. The comparison of multiple measurement series confirmed this finding, with 98% of the metabolites showing absolute time shifts less than 0.04 min under alkaline conditions and 80% under acidic conditions. An overview of analyte retention times is given in Table 1. A compilation of MRM chromatograms of all 56 metabolites measured under acidic and alkaline mobile phase conditions is included in the Supplementary material.
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
7
Fig.1. Extracted ion chromatograms obtained with optimized HILIC method under alkaline (pH 9.2) and acidic (pH 5.6) mobile phase conditions focusing on isobaric compounds frequently present in intracellular matrices. Citrate/Isocitrate mixtures of 5 and 50 lM were injected at pHs 9.2 and 5.6, respectively.
Evaluation of linearity and method detection limits for different mobile phase pHs The linearity range for each metabolite was determined under acidic (pH 5.6) and alkaline (pH 9.2) mobile phase conditions. A range of metabolite concentrations was considered linear when the squared correlation coefficient (R2) was better than 0.98 for the average of three calibration curves. To increase the predictive significance of the measurements, internal standards were considered for monitoring of instrumental fluctuations but not for normalization of obtained peak areas. Extending the calibration range from 0.5 to 200 lM, the linearity of instrument response with respect to the amount of the injected analyte was revealed to be specific for each metabolite, sample context, and mobile phase condition. Fig. 2 shows two different metabolite standards with exemplary deviations from linearity of applied calibration lines. Obtained metabolic-specific linearity ranges differ significantly at lower concentrations (1 nM to 5 lM) between both conditions. The median for the lower linearity limits was approximately 30 nM for alkaline and 120 nM for acidic mobile phase conditions. Within alkaline conditions, 40% of the metabolites show a linearity range between 5 and 20 nM as lower boundary and 5 lM as upper limit. In contrast to this, only 12% of the metabolites fulfilled the same criterion under acidic conditions. An overview of metabolic-specific linearity ranges for acidic and alkaline mobile phase conditions is shown in Table 1. The method detection limit (MDL) was calculated as the amount of metabolite that results at statistically significant peak areas above background noise. MDL was based on the standard deviation of triplicates at the lower linearity boundary (lowest concentration of the linearity range) multiplied by the expansion coefficient (t = 6.96) for defining the 99% confidence level [36]. In column MDLs were 0.01 to 3.7 pmol for alkaline conditions and 0.02 to 14.7 pmol for acidic conditions. Fully 89% of the analyzed metabolites showed MDLs lower than 100 nM (0.5 pmol on column) under alkaline conditions, whereas just 46% of the metabolites fulfilled this criterion under acidic conditions (Fig. 3). The lowest detection limit was achieved for malic acid in positive ionization mode (2.23 nM, 0.011 pmol on column) using basic conditions. An overview of the metabolite-specific MDLs for acidic and alkaline mobile phase conditions is shown in Table 1. Comparing both conditions, the most significant differences were observed for the di- and triphosphate nucleotides and NADP as well as for
fumaric acid, isocitrate, aspartic acid, and glutamine, showing approximately 20 to 50 times lower detection limits under alkaline conditions. ESI–MS responsivity under alkaline and acidic mobile phases Next, the impact of pH of mobile phases on ESI–MS responsivity was studied based on comparison of peak areas obtained by injection of 5-lM metabolite standards under both conditions. Special care was given to the preparation of 10-mM ammonium acetate buffer solutions to afford comparable ion strength under acidic and alkaline conditions. Either acetic acid or ammonium hydroxide was added in equimolar amounts to adjust the buffer system pH to 5.6 or 9.2. To compare both scenarios, averaged peak areas obtained by alkaline conditions were normalized by those obtained by acidic conditions. Strikingly, 70% of the metabolites presented more than 1.5-fold higher signal intensities under alkaline mobile phase conditions (Fig. 4). Absolute quantification in bacterial metabolite extracts For challenging the applicability of the approach, absolute concentrations of selected metabolites were quantified in intracellular extracts of E. coli biomasses. Because matrix effects are suspected to be a major source for inaccuracies and unreliability in ESI– MS-based quantitative metabolomics [6,16,18], related impacts were evaluated as well. Quantification strategies were adjusted based on previously estimated concentrations of the targeted intracellular metabolites. Intracellular pool sizes were detected by standard-based external calibration, isotope dilution, and standard addition for subsequent comparison. Standard-based external calibration and isotope dilution were applied by metabolic-specific calibration ranges in multicomponent mixtures. Standard addition was performed for each metabolite at concentrations corresponding to approximately 2-, 3-, and 4-fold the estimated values. Results obtained from exemplary measurements of six intracellular metabolites with three different approaches are summarized in Fig. 5. Corresponding correlation coefficients are shown in Table 2. Good (R2 P 0.99) to excellent (R2 P 0.999) linearity coefficients were found for all three calibration strategies. Four replicates were made for sample analysis, where calculated relative standard deviations of the measured concentrations were less than 9% for each metabolite. As depicted in Fig. 5, individual sample analysis
Metabolite
Polarity mode
Linearity range
Method detection limit
Alkaline
Acid
Alkaline
Acid
(min)
(min)
4.2 606.0 138.0 8.5 59.6 513.3 454.1 41.6 181.9 1131.1 1713.6 67.5 67.5 53.2 79.3 2949.3 423.1 532.7 148.4 148.4 13.1 1495.0 144.1 1428.0 1130.4 198.5 2274.7 1525.2 160.8 1030.1 2004.1 123.4 1094.0 1711.8 375.8 3.6 18.2 22.0
0.10 0.17 0.59 0.01 0.05 2.20 0.11 0.02 0.15 3.72 – 0.77 0.77 0.18 0.17 1.80 0.27 0.43 0.44 0.44 0.04 0.16 0.07 0.20 0.30 0.25 0.35 0.51 0.25 0.28 0.26 0.02 0.20 0.36 0.30 0.21 0.06 0.05
0.02 3.03 0.69 0.04 0.30 2.57 2.27 0.21 0.91 5.66 8.57 0.34 0.34 0.27 0.40 14.75 2.12 2.66 0.74 0.74 0.07 7.47 0.72 7.14 5.65 0.99 11.37 7.63 0.80 5.15 10.02 0.62 5.47 8.56 1.88 0.02 0.09 0.11
8.02 ± 0.04 22.07 ± 0.02 20.59 ± 0.02 21.83 ± 0.01 21.42 ± 0.02 24.15 ± 0.01 24.83 ± 0.01 23.08 ± 0.01 21.15 ± 0.01 24.03 ± 0.01 – 21.11 ± 0.02 21.27 ± 0.01 21.63 ± 0.02 22.49 ± 0.02 24.99 ± 0.01 20.77 ± 0.01 23.61 ± 0.02 22.86 ± 0.02 23.11 ± 0.02 19.64 ± 0.01 23.15 ± 0.01 19.41 ± 0.01 21.31 ± 0.01 22.80 ± 0.02 22.07 ± 0.01 23.95 ± 0.02 25.45 ± 0.01 21.38 ± 0.01 23.21 ± 0.01 24.66 ± 0.01 20.54 ± 0.01 22.57 ± 0.01 24.11 ± 0.01 15.35 ± 0.01 15.86 ± 0.01 17.64 ± 0.01 17.51 ± 0.01
8.63 ± 0.07 22.57 ± 0.01 20.12 ± 0.01 22.14 ± 0.01 21.77 ± 0.01 23.97 ± 0.02 24.75 ± 0.02 22.76 ± 0.01 21.14 ± 0.01 24.53 ± 0.01 21.73 ± 0.11 – – 21.81 ± 0.01 22.72 ± 0.01 25.93 ± 0.07 20.78 ± 0.01 24.07 ± 0.02 23.32c ± 0.02 23.61c ± 0.02 20.03 ± 0.01 23.51 ± 0.01 19.33 ± 0.01 22.08 ± 0.03 23.75 ± 0.05 22.11 ± 0.01 24.91 ± 0.03 26.73c ± 0.12 21.34 ± 0.01 24.02 ± 0.01 25.90 ± 0.17 20.40 ± 0.01 23.33 ± 0.02 25.24 ± 0.01 15.65 ± 0.01 16.26 ± 0.01 18.02 ± 0.01 17.92 ± 0.01
33.3 14.0 4.9 28.5 41.0
286.0 38.3 22.6 99.2 167.5
0.17 0.07 0.02 0.14 0.20
1.43 0.19 0.11 0.50 0.84
15.76 ± 0.01 18.69 ± 0.01 21.45 ± 0.01 17.96 ± 0.01 16.91 ± 0.01
15.00 ± 0.01 19.08 ± 0.01 22.10 ± 0.01 18.37 ± 0.01 17.32 ± 0.01
0.9964
9.9
140.9
0.05
0.70
18.45 ± 0.01
18.90 ± 0.01
0.9968
9.1
66.8
0.05
0.33
14.37 ± 0.01
14.76 ± 0.01
5000
0.9964
26.4
49.7
0.13
0.25
14.54 ± 0.01
14.94 ± 0.01
5000
0.9967
43.1
97.6
0.22
0.49
17.14 ± 0.01
17.58 ± 0.02
20
5000
0.9965
7.9
35.9
0.04
0.18
12.74 ± 0.01
13.15 ± 0.01
0.9990
20
5000
0.9966
17.2
22.2
0.09
0.11
13.22 ± 0.01
13.64 ± 0.01
0.9991
100
5000
0.9964
15.0
100.6
0.08
0.50
18.35 ± 0.01
18.79 ± 0.01
R2
Min (nM)
Max (nM)
R2
(nM)
( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (+) (+) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (+) (+) (+)
50 50 500 20 50 50 20 20 50 1000 – 20 20 50 50 500 200 50 100 100 10 20 50 100 100 100 100 500 100 50 200 100 50 100 100 20 10 20
5000 5000 5000 5000 5000 5000 5000 1000 5000 5000 – 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 2000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000
1.0000 0.9999 0.9999 0.9998 0.9974 0.9938 0.9953 0.9899 0.9884 0.9903 – 0.9973 0.9973 0.9954 0.9972 0.9998 0.9990 0.9904 0.9874 0.9874 0.9982 0.9901 0.9995 0.9928 0.9925 0.9990 0.9844 0.9915 0.9990 0.9905 0.9915 0.9991 0.9858 0.9968 0.9970 0.9996 0.9994 0.9990
50 500 500 20 100 1000 500 50 200 1000 2000 50 50 100 100 2000 500 500 500 500 20 500 100 1000 1000 100 1000 2000 100 1000 1000 100 500 2000 100 50 20 50
5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000 5000
0.9998 0.9951 0.9976 0.9998 0.9916 0.9838 0.9945 0.9854 0.9953 0.9950 0.9972 0.9956 0.9956 0.9990 0.9992 0.9980 0.9985 0.9994 0.9984 0.9984 0.9963 0.9847 0.9987 0.9832 0.9836 0.9960 0.9975 0.9971 0.9974 0.9938 0.9842 0.9985 0.9840 0.9850 0.9988 0.9938 0.9920 0.9963
20.7 33.8 117.1 2.2 9.4 439.6 21.5 3.9 30.1 744.3 – 154.5 154.5 35.1 33.5 360.8 53.9 86.5 87.8 87.8 7.6 31.4 13.4 40.1 59.6 49.9 70.9 102.4 49.9 56.3 51.9 4.7 39.2 72.8 59.3 41.7 12.7 10.2
L-Ala
(+) (+) (+) (+) (+)
20 10 10 200 100
5000 5000 5000 5000 5000
0.9996 0.9991 0.9953 0.9992 0.9990
200 50 50 500 200
5000 5000 5000 5000 5000
0.9986 0.9956 0.9938 0.9963 0.9956
L-Ser
(+)
100
5000
0.9988
200
5000
L-Pro
(+)
10
5000
0.9992
50
5000
L-Val
(+)
20
5000
0.9991
100
L-Thr
(+)
50
5000
0.9993
100
L-Leu
(+)
20
5000
0.9993
L-Ile
(+)
20
5000
L-Asn
(+)
50
5000
L-Hom-Ser
O-Ac-Ser OSHS lL-Cysta Gly
Alkaline
Acid
Alkaline HILIC MS metabolite quantification / A. Teleki et al. / Anal. Biochem. 475 (2015) 4–13
Alkaline
Max (nM)
aKB Fum Suc Mal aKG Cit IsoCit cis-Aco Xu5P 6PG E4Pb R5Pa Ru5Pa F6P G6P FbP DHAP PEP 2PGa 3PGa NAD NADP AMP ADP ATP GMP GDP GTP CMP CDP CTP UMP UDP UTP cAMP AIBA GABA
Retention time
(pmol, on column)
Min (nM)
Acid
8
Table 1 Overview of obtained linearity ranges (
