Glycolysis is a defining metabolic marker in multiple biological systems. Monitoring glycolysis by measuring the extracellular flux of H+ is common, but requires correction to be quantitative and unambiguous. Here, we demonstrate how to gather and correct extracellular flux data to distinguish between respiratory and glycolytic sources of extracellular acidification.
Extracellular measurement of oxygen consumption and acid production is a simple and powerful way to monitor rates of respiration and glycolysis1. Both mitochondrial (respiration) and non-mitochondrial (other redox) reactions consume oxygen, but these reactions can be easily distinguished by chemical inhibition of mitochondrial respiration. However, while mitochondrial oxygen consumption is an unambiguous and direct measurement of respiration rate2, the same is not true for extracellular acid production and its relationship to glycolytic rate 3-6. Extracellular acid produced by cells is derived from both lactate, produced by anaerobic glycolysis, and CO2, produced in the citric acid cycle during respiration. For glycolysis, the conversion of glucose to lactate– + H+ and the export of products into the assay medium is the source of glycolytic acidification. For respiration, the export of CO2, hydration to H2CO3 and dissociation to HCO3– + H+ is the source of respiratory acidification. The proportions of glycolytic and respiratory acidification depend on the experimental conditions, including cell type and substrate(s) provided, and can range from nearly 100% glycolytic acidification to nearly 100% respiratory acidification 6. Here, we demonstrate the data collection and calculation methods needed to determine respiratory and glycolytic contributions to total extracellular acidification by whole cells in culture using C2C12 myoblast cells as a model.
The overall goal of this method is to accurately measure the glycolytic rate of cells using extracellular flux analysis. Quantitative measurement of glycolytic rate using extracellular acidification is the desired endpoint of many experiments. However, the total rate of extracellular acidification is the sum of two components: respiratory acidification, in the form of CO2 (which hydrates to H2CO3 then dissociates to HCO3– + H+), and glycolytic acidification, in the form of lactate– + H+.
The contributions of CO2 to total extracellular acidification have until recently been considered negligible in the measurement platform used here, the XF24 analyzer 7. However, it is clear in multiple other systems that CO2 can be a major contributor to extracellular acidification4-5. Multiple papers acknowledge this contribution, but do not attempt direct quantitation of CO2-derived acid 3,8,9. We recently demonstrated quantitatively that CO2 production is a significant source of extracellular acidification in this system 6. Moreover, though there are multiple metabolic pathways that generate CO2 from glucose catabolism, those carried out by matrix dehydrogenases in the citric acid cycle are the overwhelming contributors and all other sources generate amounts of CO2 that are within experimental error 6.
Without correcting for CO2 production, extracellular acidification is therefore an ambiguous indicator of glycolytic rate and cannot be used quantitatively. Our previous publication highlights several instances where respiratory CO2 comprises the bulk of the total acidification signal, even in cells generally believed to primarily use glycolysis6. Additionally, the respiratory CO2 contribution to total acidification varies widely during the course of common metabolic profiling experiments, demonstrating that correct comparison of the glycolytic rate during different parts of an experiment requires correction for CO2.
To measure the glycolytic rate of cells using the rate of extracellular acidification, it is necessary to convert pH changes to changes in total H+ generated, and to subtract the extracellular acidification caused by CO2 released during operation of the citric acid cycle. Here, we describe a straightforward method for measuring extracellular proton production rate (from extracellular changes in pH and the calibrated buffering power of the assay medium) and CO2 production (from extracellular changes in O2 concentration), and demonstrate how to calculate glycolytic rate using these measurements.
This method strengthens the utility of extracellular acidification measurement by using it to properly calculate glycolytic rate as defined by lactate production. Without correction for respiratory CO2 (or direct measurement of lactate), it is impossible to determine if and to what extent the total acidification rate reflects glycolytic rate, confounding the interpretation of experiments that use total extracellular acidification as a direct measurement of lactate production.
CALCULATIONS
CO2 and lactate are, within experimental error, the only two contributors to extracellular acid production, based on experiments with myoblast cells6. Therefore, the rate of total extracellular acidification (PPR, proton production rate) can be defined as:
PPRtot = PPRresp + PPRglyc Equation 1
where tot = total; resp = respiratory; glyc = glycolytic. Glycolytic PPR is thus:
PPRglyc = PPRtot – PPRresp Equation 2
Here,
PPRtot = ECARtot/BP Equation 3
where ECAR = extracellular acidification rate (mpH/min), and BP = buffering power (mpH/pmol H+ in 7 µl), while
PPRresp = (10pH-pK1/(1+10pH-pK1))(max H+/O2)(OCRtot – OCRrot/myx) Equation 4
where K1 = combined equilibrium constant of CO2 hydration and dissociation to HCO3– + H+; max H+/O2 = the CO2-derived acidification for a particular metabolic transformation such as complete oxidation of glucose6; OCR = oxygen consumption rate (pmol O2/min), and OCRrot/myx = non-mitochondrial OCR.
Equation 4 isolates mitochondrial OCR by subtracting any non-mitochondrial OCR (defined as OCR that is resistant to the mitochondrial respiratory poisons rotenone and myxothiazol) and accounts for the maximum H+ generated per O2 consumed for each substrate (max H+/O2) (see 6), as well as the proportion of CO2 giving rise to H+ at the experimental temperature and pH (10pH-pK1/(1+10pH-pK1). For full oxidation of glucose, mitochondrial Oxygen Consumption Rate (OCR) is exactly equal to the rate of CO2 production. In the confined assay volume of extracellular flux measurement, CO2 produced by respiration remains trapped in the assay medium. Most of the trapped CO2 is hydrated to H2CO3, which then dissociates to HCO3– + H+. A small fraction remains dissolved but not hydrated, and another small fraction is hydrated but not dissociated, as dictated thermodynamically by the combined equilibrium constant of CO2 hydration and dissociation to HCO3– + H+ at experimental temperature (37 °C) and pH (~7.4).
Thus, the complete equation for calculating PPRg by subtracting PPRresp from PPRtot is:
PPRglyc = ECARtot/BP – (10pH-pK1/(1+10pH-pK1))(max H+/O2)(OCRtot – OCRrot/myx) Equation 5
In this way, rates of respiration and glycolysis, as well as their associated ATP production rates, can be quantitatively determined from straightforward measurements (oxygen consumption, extracellular acidification, buffering capacity) and import or calculation of other required values (H+/O2, P/O, and the equilibrium constant K1) 6. The experiment described here expands on standard techniques for using the Extracellular Flux Analyzer such as Seahorse XF24 10,11; for other extracellular flux measurement formats (e.g., XFe96, or XFp), all volumes below should be scaled appropriately.
The buffering power of the assay medium can be measured by construction of a standard curve either directly in the extracellular flux platform or separately using a calibrated pH probe. Here, three options for measuring buffering by the extracellular flux assay medium are given, including using all injection ports of the extracellular flux analyzer with cell-free sample wells, or using only the last injection port in cell-containing wells (section 1) or by using an external pH measurement (section 2). See the attached spreadsheet for the full calculations of example data.
To measure buffering power using the pH-detecting capability of the extracellular flux instrument, it is safest to use cell-free wells to minimize signal variation. However, within the error, no statistical difference exists between cell-free and cell-containing wells when performing this measurement (data not shown). NOTE: The variation described in step 1.7 carries the advantage of accounting for any potential changes to buffering conferred by added compounds or by the presence of cells, with the disadvantage of noisier signal. However, as stated above, no significant differences were found in the calculated buffering power between the cell-free design shown in Table 1 and the post-experiment design in Table 2 under the experimental conditions described here.
Additionally, over small ΔpH ranges (<0.4 units; experimentally best restricted to 0.2 units), the linear slope obtained by plotting Δ mpH/pmol H+ adequately approximates the logarithmic relationship between ΔpH and [H+]. The slope of this standard curve therefore represents the buffering power of the assay medium under test in pH/nmol H+ in 7 µl, or mpH/pmol H+ in 7 µl. We recommend increasing medium buffering power or decreasing cell density for samples that exceed a 0.2 pH unit change during the measurement time. The measurement time may also be decreased, but this may shorten the steady state acidification rate and introduce error into the rate calculation.
1. Measuring Buffering Power in an Extracellular Flux Instrument: Two Variations
NOTE: the calculations and methods described here were developed using an Extracellular Flux Analyzer. For other instruments, the measurement volume must be scaled appropriately.
Table 1. Consecutive HCl injections into an extracellular flux assay well.
Table 2. Single HCl injection into an extracellular flux assay well.
2. Measuring Buffering Power Using an External pH Meter
NOTE: To measure the buffering power of a medium using an external pH probe, calibrate the probe at 37 °C and maintain this temperature for all reagents during the experiment.
Table 3. Measuring buffering power using a pH meter. Data represent a typical experiment with six 20 µl additions of 0.1 M HCl.
Figure 1. Determining buffering power. HCl standard curve measured as in Table 1, Table 2 or (here) as in Table 3. The slope of the linear curve fit gives the buffering power (pH/nmol H+ in 7 µl). Each point represents mean ± SEM of n = 9 technical replicates.
Table 4. Buffering power and buffering capacity of selected media.
3. Performing an Extracellular Flux Assay Using C2C12 Myoblast Cells
NOTE: In step 3.4.3, there were no observed differences in CO2-derived acid production dependent on the presence of carbonic anhydrase in C2C12 culture, suggesting that its presence is not required for full conversion of CO2 to HCO3– + H+. However, empirically testing this in different experimental systems is recommended before omitting carbonic anhydrase.
NOTE: For each segment of the experiment, determine the mix, wait, and measurement times desired, as well as the number of cycles per segment.
NOTE: The data in Table 5 were collected over two assay cycles of 2 min mix, 1 min wait, and 5 min measure for each segment, with three assay cycles occurring after the Port D addition of different amounts of HCl (for calibration of buffering power as in Table 2).
Table 5. Extracellular flux assay configuration.
4. Measuring End-point Lactate Concentration
NOTE: To validate the indirect assay described here in some different system, end point lactate concentration at the end of an extracellular flux experiment can be determined directly in a conventional 96-well plate by measuring the initial velocity (over 2 min) of reduction of NAD+ → NADH catalyzed by lactate dehydrogenase, described in detail in our prior publication6. For the data presented in Representative Results, the end point lactate concentration in glucose-containing assay wells was ~40 μM.
5. Measuring Protein Content
Figure 2 shows the raw data for a typical experiment. The last 10 measurement points from the point-to-point recording of both OCR and pH (shaded vertical bars) were used for the calculations. Initial concerns that the average value (middle point measurement) of each assay cycle would not provide sufficient resolution of rate for an accurate calculation, particularly as there appeared to be a slight lag between port addition and steady state acidification rate, were not borne out, as this does not appear to contribute significantly to calculation error (not shown). Alternatively, if the correct buffering capacity is entered during experimental setup, PPR can be read directly from the instrument data collection readout by displaying the PPR output in the instrument software or in the PC-compatible format available as one of the data output settings.
Figure 2. Representative extracellular flux traces of O2 and H+. OCR and pH traces for the example experiment in Table 5, containing 10 mM glucose at assay start. Port D had different HCl concentrations for calibration of buffering power (not shown in these averaged traces). Data from previous publication6. Each point represents mean ± SEM of n = 8 biological replicates. Please click here to view a larger version of this figure.
Data analysis of representative results
Using the spreadsheet shown in Table 6 and provided as an attachment, data values from individual wells may be entered in the columns shown with yellow headers. All six columns to the right are calculated from these entries. The example in Table 6 shows the calculations of PPRresp and PPRglyc using ECAR and OCR data from individual wells for the native conditions with or without added glucose, prior to Port A addition of oligomycin. Technical replicates on each biological preparation are normally averaged to give single values of the outputs in the last four columns, then data from different biological preparations are averaged with appropriate propagation of error statistics in BP and these four values.
Table 6. Calculation of respiratory and glycolytic acidification. Columns headed in yellow indicate values to be entered from calculation (e.g., BP, max H+/O2), or from data collection (e.g., ECARtot, OCR). Please click here to view a larger version of this table. | Please click here to download this table as an Excel spreadsheet.
Contributions of glycolysis and respiration to PPR after correction
Figure 3 shows the graphical output of data calculated as in Table 6 for native rates of glycolytic and respiratory acidification, rates following oligomycin addition (Port A), and rates following FCCP addition (Port B). These data clearly demonstrate how the proportions of respiratory and glycolytic acidification change with choice of substrate (glucose vs. control (ctl) with none added) and with mitochondrial status (native function vs. pharmacologically altered function).
Figure 3. Proton production rate (PPR) from glycolytic and respiratory sources. PPR from respiration (open portions) and glycolysis (filled portions) of C2C12 cells calculated using Equation 5 with added glucose (three left bars) or without added glucose (three right bars). Data from6. All data represent mean ± SEM of n = 8 biological replicates.
Extracellular acidification is an easily measured indication of cellular metabolic rate. To properly determine the rate of cellular glycolysis (as defined by lactate production) it is critical to know the buffering power of the assay medium, and to convert the extracellular flux measurements of oxygen consumption and acidification to proton production rates. By performing this calculation, the acidification resulting from CO2 released in the citric acid cycle can be subtracted, leaving the acidification that results from lactate production.
The multiple different ways given here to measure buffering power for this correction carry different advantages and disadvantages. External measurement using a pH probe is highly accurate and reproducible, but may not reflect small differences in pH detection introduced by the fluorophores contained within the assay plate, the addition of compounds during the assay, or the presence of the cells themselves. The in-plate pH measurements address these issues, but also introduce varying degrees of experimental noise.
The CO2 correction to ECAR allows for the first time the unambiguous and quantitative calculation of glycolytic rate, and reveals variation in respiratory and glycolytic contribution to total ECAR during the course of an experiment. Using Equation 5 and the measurements of OCR, ECAR, and buffering power, glycolytic rate can be calculated using the simple spreadsheet provided (Table 6). This rate can be verified by post-hoc lactate measurement if desired 6. In cells where the pentose phosphate pathway is highly active, the use of pathway inhibitors such as 6-aminonicotinamide may be useful to isolate glycolytic rate. Calculation of the contributions of both CO2– and lactate-derived H+ from the total measured Extra Cellular Acidification Rate and Oxygen Consumption Rate is an invaluable tool for using extracellular flux data to make powerful and quantitative statements about metabolic activity.
Using the procedures described here, including various modifications for measuring buffering power, and optimizing the extracellular flux experiment for the cells under investigation and data desired, the rate of glycolysis in intact cells can be quantified under a wide range of experimental conditions. This method is limited to cells that can grow in adherent culture on (or cells or organelles that can be adhered to) a polystyrene surface. It is most reliable when cultured cells are homogenous and confluent, though useful data may still be obtained over a range of these conditions. The calculations require some knowledge of the metabolism of the cells, as max H+/O2 ranges from 0.65 to 1.0 for full oxidation of different substrates and more for partial oxidation 6, however, if the cells are known to oxidize glucose, a value of 1.0 can be assumed.
Though relevant to all metabolic characterization, this method may be particularly helpful when used in systems in which a shift between respiratory and glycolytic metabolism to maintain cellular ATP supply is a critical phenotype, including the characterization of stem cells and tumor-derived cancer cells. Understanding metabolic control alterations in these and other contexts will allow a greater degree of sophistication and accuracy in the experimental design and analysis of these cell types.
The authors have nothing to disclose.
We thank David A. Ferrick and David G. Nicholls for contributing to project conception and presentation, Renata L.S. Goncalves and Akos A. Gerencser for data not shown here and for helpful discussions, Barbara Liepe for XF24 consumables, and Andy Neilson for input in developing Eq. (5).
Pherastar FS | BMG | n/a | microplate reader |
Seahorse XF-24 | Seahorse Bioscience | n/a | extracellular flux instrument |
Seahorse XF assay plate | Seahorse Bioscience | V7-PS | consumable |
XF Calibrant | Seahorse Bioscience | 100840-000 | solution |
HCl standard | Sigma | 38280 | chemical |
oligomycin | Sigma | O4876 | chemical |
FCCP | Sigma | C2920 | chemical |
Rotenone | Sigma | R8875 | chemical |
Myxothiazol | Sigma | T5580 | chemical |
DMEM | Corning | 10-013-CV | medium component |
FBS | Corning | 35-010-CV | medium component |
penicillin/streptomycin | Corning | 30-002-CI | medium component |
carbonic anhydrase | Sigma | C2624 | chemical |
96-well assay plate | Corning | CLS3991 | consumable |
NAD+ | Sigma | N7004 | chemical |
LDH | Sigma | L1254 | chemical |