Metabolic mechanisms of a drug revealed by distortion-free <sup>13</sup>C tracer analysis
Metabolomic isotopic tracing can provide flux information useful for understanding drug mechanisms. For that, NMR has the unique advantage of giving positional isotope enrichment information, but the current 13 C 1D NMR approach suffers from low sensitivity and high overlaps. We developed a new 2D heteronuclear NMR experiment incorporating J-scaling and distortion-free elements that allows for quantitative analysis of multiplets with high sensitivity and resolution. When applied to an old chemotherapeutic drug, the approach provided a quantitative estimation of TCA-cycle turns, confirming the conventional mechanism of its mitochondrial metabolic enhancement. Additionally, the approach identified a new mechanism of the higher contribution of the pentose phosphate pathway to serine synthesis in the cytosolic compartment, possibly explaining the broad pharmacological activities of the drug. Our approach may prove beneficial in helping to find new usages or metabolic mechanisms of other drugs.
Introduction<!>Design of a new J-scaled distortion-free 2D heteronuclear NMR experiment<!>Quantitative and distortion improvement evaluation of new Jscaled HSQC<!>Estimation of TCA-cycle ux and effects of an old chemotherapeutic drug<!>Identication of the new mechanism of DCA for serine biosynthesis<!>Conclusions<!>Conflicts of interest
<p>For the last two decades, the development of metabolomics has brought about a new paradigm wherein abnormal metabolism is considered a fundamental factor underlying various diseases including cancer as well as immune and neurodegenerative diseases. [1][2][3] Whereas metabolomics measures concentrations of metabolites, it does not always provide the ux, which is the actual indicator of the efficiency of metabolic pathways. Flux information can be obtained by isotopomer analysis using stable isotopic tracers, mostly 13 C, with liquid chromatography-mass spectrometry (LC-MS) 4 or nuclear magnetic resonance (NMR) spectroscopy. 5 Related to the latter, magnetic resonance spectroscopy 6 and hyperpolarized magnetic resonance imaging 7 are being applied for ux studies. Compared with LC-MS, the merits of NMR include provision of positional information of the 13 C atom from a tracer, reproducibility, and structural information of the metabolites. However, it suffers lower sensitivity, which is even more problematic for 13 C-NMR that has been the major approach in NMR isotopomeric analysis. 5,8,9 In addition, the 1D 13 C-NMR oen has insufficient resolution for the complex mixtures used in metabolomic settings. Related hyperpolarization techniques could overcome the sensitivity issue for short-term ux analysis, but long-term and detailed isotopomeric analysis remains challenging. Therefore, an approach that can yield high-sensitivity, high-resolution positional isotope enrichment information has been highly desired.</p><p>A simple solution for the 13 C sensitivity issue would be 1 Hdetected heteronuclear 2D NMR. However, there are two major hurdles for routine use of this approach. One is the poor resolution of the indirect domain due to the nature of the detection and limited experimental time. This problem has been partly addressed using non-uniform sampling (NUS), 10,11 but only for qualitative estimation of metabolism. In addition, NUS does not give as linear a relationship as straightforward Fourier transform, thus incurring difficulties in quantitative analysis. Another problem is the phase distortion within 13 C multiplets in routine 2D heteronuclear NMR experiments with 13 C tracers; 12 this prevents quantitative analysis of individual multiplets required for detailed metabolic pathway analysis. These factors have made the 2D approach rarely usable in most quantitative isotopic tracer experiments, and 13 C-NMR has been employed instead, despite the disadvantages.</p><p>If properly applied, isotopomer analysis can provide valuable insights into mechanistic targets of drugs. Actually, most drugs have multiple targets, either known or unknown, as termed by "polypharmacology", 13 and nding new targets of existing drugs can lead to wider uses of those drugs. As drugs modulating metabolism have higher chances of multiple targets due to the intricately connected metabolic pathways, 14,15 an efficient means of investigating metabolic drug mechanisms will certainly benet drug usage and development.</p><!><p>We designed a new NMR experiment that can provide quantitative metabolic assessment while addressing the problems of 2D NMR stated above. First, the low-resolution problem was addressed with J-scaling, 16 which gives wider separation within multiplets (higher resolution) (Fig. 1a). As the J-scaling is used with Fourier transform, it retains the linear relationship, obviating the non-quantitative NUS approach for resolution enhancement. The second issue of phase distortion occurs due to the homonuclear J CC evolution during the long time block (4.4 ms) for t 1 -noise reduction and frequency discrimination (t 1 ) (Fig. 1b and S1 †). We, therefore, moved the gradient and long carbon pulses out of the period and used a pair of short adiabatic inversion pulses of only 300 ms duration (total 600 ms), thereby practically eliminating the deleterious J CC evolution (Fig. 1c). Optimally incorporating these elements into a 2D pulse sequence gave a new sensitive and distortion-free NMR experiment (Fig. S2 †). Analytical simulation of this new experiment provided essentially distortion-free and well-resolved multiplets for the usual J CC values (30-60 Hz) (Fig. S3 †).</p><!><p>We rst tested the experiment on a frequently used 13 C metabolic tracer, U-13 C 3 -lactate. A conventional 2D heteronuclear single quantum coherence (HSQC) spectrum gave a fused doublet with distorted peak shapes (Fig. 2a). In comparison, our distortion-free approach gave proportionally separated doublets according to the scaling factor, thus demonstrating the improvement in the effective resolution. Accordingly, the vertical 1D slice of the original HSQC exhibited a large negative dip in the center, whereas those from our approach exhibited a clean baseline with regular peak shapes (Fig. 2b). Therefore, it was clear that quantitative assessment of individual peaklets is possible only with those obtained with the new approach. For the optimal J-scaling factor, on one hand, a smaller value may result in insufficient separation of multiplets; on the other hand, a larger scaling factor may cause peaklet-overlaps with those from other peaks or unnecessary splitting by long-range couplings. In addition, the nal resolution of the J-scaling approach is also affected by indirect time domain points, spectral width, and the magnetic eld strength. For usual metabolomic isotopomer analysis in our NMR setting, the scaling factor of 6 resolved most of the important multiplets without pronounced overlaps or long-range splitting.</p><p>Next, actual quantitation of isotopomers was performed. For this, a 1 : 1 mixture of isotopomers of alpha-ketoglutarate (aKG) (U-13 C 5 -aKG and 1,2,3,4-13 C 4 -aKG), an important intermediate in the TCA-cycle, was used. Each isotopomer should contribute to the nal spectrum as a quartet and a doublet of C4, respectively, giving a total of 6 peaklets (Fig. 3a). The C4 doublet to quartet ratio (Q/D ratio) is an important metabolic parameter that has been conventionally obtained by 13 C NMR isotopomer analysis. 5 Therefore, we tried it using our 2D experiment with an order-of-magnitude-higher sensitivity and also compared it with a NUS approach that gives high resolution. 10 The original HSQC could not resolve the individual peaklets, hampering any  further analysis (Fig. 3b). The NUS approach, although it gave the expected 6 peaklets, exhibited negative peaks from phase distortion. In contrast, our new approach gave readily identiable in-phase 2D peaks with the expected values of peak separation. The 1D slices of these spectra along the carbon dimension gave a clearer idea of these differences (Fig. 3c). As a quantitative estimation of the Q/D ratio, the NUS gave 1 : 1.417, whereas our approach gave 1 : 1.045, much closer to the theoretical ratio of 1 : 1 (Table S1 †). Furthermore, the intensities of the individual Q peaklets in the NUS approach severely deviated from the uniformity with the ratio of 1.15 : 1 : 2.223 : 1.379, whereas that of our approach were much more homogenous (1.023 : 1 : 0.887 : 1.080). These results show that the distortion-free approach can yield much better quantitative estimation of isotopomer distribution. In addition, as the method provide well-resolved multiplets, it could be used during spiking procedures for easier metabolite identication with multiple reference peaklets.</p><!><p>Then, we applied the approach to the mechanistic study of a drug. To that end, we treated dichloroacetate (DCA) to a leukemia cell line (L1210) and analyzed the cell extract directly without any cleanup steps. DCA is an investigational drug that has been under clinical trials in several cancers. It also has activities for seemingly unrelated diseases, such as pulmonary hypertension and amyotrophic lateral sclerosis. 17 Its mechanism in cancer has been largely attributed to its inhibition of pyruvate dehydrogenase kinase (PDK), which effectively increases the ux of pyruvate into the mitochondrial TCA-cycle. However, its activities for other diverse conditions suggest that additional mechanisms may also exist.</p><p>First, we investigated whether DCA actually enhances TCAcycle ux, using the splitting patterns of glutamate C4 derived from U-13 C 6 -glucose. The high sensitivity of our approach enabled us to obtain high-quality spectra from the complex cell extracts. The spectra exhibited a mixture of a strong doublet (from 13 C-4,5; J ¼ 52 Hz) and a quartet (from 13 C-3,4,5; J ¼ 52, 34 Hz) of glutamate generated from single and multiple turns of the TCA-cycle, respectively (Fig. 4a and b). With conventional HSQC, contrastingly, it was not possible to observe the splitting. There are calculable relationships among the doublet/quartet ratio of the C4, the TCA-cycle turns, and the 13 C enrichment of acetyl-CoA that can provide quantitative estimation on average TCA-cycle turns (Fig. 4c). Based on the isotopologue abundance of acetyl-CoA (55.3% for control, 42.7% for DCA group (Fig. S4 †)) and the above relationship, the effects of DCA on the average TCA-cycle turns could be obtained from the doublet/quartet ratio of the C4. This analysis showed that the TCA-cycle turns increased by approximately two-fold in leukemia cells treated with DCA (Fig. 4d). Therefore, our results support the conventional mechanism of DCA, that is, increase in TCA-cycle efficiency in more quantitative terms. Methodology-wise, an optimization for the shorter incubation time should help more accurate estimation of the TCA-cycle turns in the le region of the curves where D/Q ratio is the most sensitive to the TCA-cycle turns.</p><!><p>The increase in average TCA ux, though consistent with the conventional anticancer mechanism of DCA, may still not account for all of the broad-reaching pharmacological effects of DCA. Therefore, another action mechanism of DCA was investigated, by analyzing the splitting patterns of other metabolites. A surprising labeling pattern was observed for an amino acid serine with a mixture of strong singlet and weak doublet (J ¼ 38.4 Hz) for its C3 carbon (Fig. 5a). The doublet can be easily explained by the conventional serine biosynthesis branching from glycolysis, as is clearly seen for the C3 of lactate, which is also synthesized through glycolysis (Fig. 5b and S5 †). By contrast, the much stronger singlet on serine C3, absent in the lactate spectrum, suggests that a large portion of the neighboring 13 C-C2 was replaced by unlabeled 12 C-carbon. Based on known serine biosynthesis, the only plausible route for this carbon reshuffling is the pentose phosphate pathway (PPP) (Fig. S6 †). The singlet-to-doublet ratio, therefore, should give the contribution of the PPP pathway relative to glycolysis in serine biosynthesis. The DCA treatment actually increased the singlet-to-doublet ratio by about 40% (Fig. 5c), indicating that DCA increases the relative contribution of PPP over glycolysis in serine synthesis. This may constitute a new action mechanism of DCA, as serine synthesis modulation recently has been linked to several diseases such as cancer, immune and neurological diseases, as noted above. [18][19][20] Our results also show that cytosol, where PPP occurs, may be a new subcellular target compartment of DCA, distinct from the mitochondrial compartment for the conventional PDK inhibition mechanism. This new target site and mechanism may help to nd new indications for DCA.</p><!><p>In conclusion, we developed a high-sensitivity 2D NMR experiment that allows for quantitative analysis of multiplets for metabolic pathway analysis. Application of this to DCA not only conrmed its conventional mitochondrial TCA-enhancing mechanism but also suggested a new mechanism of higher PPP contribution to serine synthesis in cytosol. The new approach should also facilitate identication of new metabolic mechanisms of other drugs.</p><!><p>There are no conicts to declare.</p>