Metabolic modeling of dynamic brain ¹³C NMR multiplet data: concepts and simulations with a two-compartment neuronal-glial model

Challenges of Investigating Compartmentalized Brain Energy Metabolism Using Nuclear Magnetic Resonance Spectroscopy in vivo

Brain function requires continuous energy supply. Thus, unraveling brain metabolic regulation is critical not only for our basic understanding of overall brain function, but also for the cellular basis of functional neuroimaging techniques. While it is known that brain energy metabolism is exquisitely compartmentalized between astrocytes and neurons, the metabolic and neuro-energetic basis of brain activity is far from fully understood. 1H nuclear magnetic resonance (NMR) spectroscopy has been widely used to detect variations in metabolite levels, including glutamate and GABA, while 13C NMR spectroscopy has been employed to study metabolic compartmentation and to determine metabolic rates coupled brain activity, focusing mainly on the component corresponding to excitatory glutamatergic neurotransmission. The rates of oxidative metabolism in neurons and astrocytes are both associated with the rate of the glutamate-glutamine cycle between neurons and astrocytes. However, any possible correlation between energy metabolism pathways and the inhibitory GABAergic neurotransmission rate in the living brain remains to be experimentally demonstrated. That is due to low GABA levels, and the consequent challenge of determining GABAergic rates in a non-invasive manner. This brief review surveys the state-of-the-art analyses of energy metabolism in neurons and astrocytes contributing to glutamate and GABA synthesis using 13C NMR spectroscopy in vivo, and identifies limitations that need to be overcome in future studies.

Magnetic resonance spectroscopy in the rodent brain: Experts' consensus recommendations

In vivo MRS is a non-invasive measurement technique used not only in humans, but also in animal models using high-field magnets. MRS enables the measurement of metabolite concentrations as well as metabolic rates and their modifications in healthy animals and disease models. Such data open the way to a deeper understanding of the underlying biochemistry, related disturbances and mechanisms taking place during or prior to symptoms and tissue changes. In this work, we focus on the main preclinical 1H, 31P and 13C MRS approaches to study brain metabolism in rodent models, with the aim of providing general experts' consensus recommendations (animal models, anesthesia, data acquisition protocols). An overview of the main practical differences in preclinical compared with clinical MRS studies is presented, as well as the additional biochemical information that can be obtained in animal models in terms of metabolite concentrations and metabolic flux measurements. The properties of high-field preclinical MRS and the technical limitations are also described.

How Energy Metabolism Supports Cerebral Function: Insights from 13C Magnetic Resonance Studies In vivo

Cerebral function is associated with exceptionally high metabolic activity, and requires continuous supply of oxygen and nutrients from the blood stream. Since the mid-twentieth century the idea that brain energy metabolism is coupled to neuronal activity has emerged, and a number of studies supported this hypothesis. Moreover, brain energy metabolism was demonstrated to be compartmentalized in neurons and astrocytes, and astrocytic glycolysis was proposed to serve the energetic demands of glutamatergic activity. Shedding light on the role of astrocytes in brain metabolism, the earlier picture of astrocytes being restricted to a scaffold-associated function in the brain is now out of date. With the development and optimization of non-invasive techniques, such as nuclear magnetic resonance spectroscopy (MRS), several groups have worked on assessing cerebral metabolism in vivo. In this context, ¹H MRS has allowed the measurements of energy metabolism-related compounds, whose concentrations can vary under different brain activation states. ¹H-[¹³C] MRS, i.e., indirect detection of signals from ¹³C-coupled ¹H, together with infusion of ¹³C-enriched glucose has provided insights into the coupling between neurotransmission and glucose oxidation. Although these techniques tackle the coupling between neuronal activity and metabolism, they lack chemical specificity and fail in providing information on neuronal and glial metabolic pathways underlying those processes. Currently, the improvement of detection modalities (i.e., direct detection of ¹³C isotopomers), the progress in building adequate mathematical models along with the increase in magnetic field strength now available render possible detailed compartmentalized metabolic flux characterization. In particular, direct ¹³C MRS offers more detailed dataset acquisitions and provides information on metabolic interactions between neurons and astrocytes, and their role in supporting neurotransmission. Here, we review state-of-the-art MR methods to study brain function and metabolism in vivo, and their contribution to the current understanding of how astrocytic energy metabolism supports glutamatergic activity and cerebral function. In this context, recent data suggests that astrocytic metabolism has been underestimated. Namely, the rate of oxidative metabolism in astrocytes is about half of that in neurons, and it can increase as much as the rate of neuronal metabolism in response to sensory stimulation.

Experimental strategies for in vivo13C NMR spectroscopy

In vivo carbon-13 (¹³C) MRS opens unique insights into the metabolism of intact organisms, and has led to major advancements in the understanding of cellular metabolism under normal and pathological conditions in various organs such as skeletal muscles, the heart, the liver and the brain. However, the technique comes at the expense of significant experimental difficulties. In this review we focus on the experimental aspects of non-hyperpolarized ¹³C MRS in vivo. Some of the enrichment strategies which have been proposed so far are described; the various MRS acquisition paradigms to measure ¹³C labeling are then presented. Finally, practical aspects of ¹³C spectral quantification are discussed.

Mechanism of antineoplastic activity of lonidamine

Lonidamine (LND) was initially introduced as an antispermatogenic agent. It was later found to have anticancer activity sensitizing tumors to chemo-, radio-, and photodynamic-therapy and hyperthermia. Although the mechanism of action remained unclear, LND treatment has been known to target metabolic pathways in cancer cells. It has been reported to alter the bioenergetics of tumor cells by inhibiting glycolysis and mitochondrial respiration, while indirect evidence suggested that it also inhibited l-lactic acid efflux from cells mediated by members of the proton-linked monocarboxylate transporter (MCT) family and also pyruvate uptake into the mitochondria by the mitochondrial pyruvate carrier (MPC). Recent studies have demonstrated that LND potently inhibits MPC activity in isolated rat liver mitochondria (K_i 2.5μM) and cooperatively inhibits l-lactate transport by MCT1, MCT2 and MCT4 expressed in Xenopus laevis oocytes with K_0.5 and Hill coefficient values of 36-40μM and 1.65-1.85, respectively. In rat heart mitochondria LND inhibited the MPC with similar potency and uncoupled oxidation of pyruvate was inhibited more effectively (IC₅₀~7μM) than other substrates including glutamate (IC₅₀~20μM). LND inhibits the succinate-ubiquinone reductase activity of respiratory Complex II without fully blocking succinate dehydrogenase activity. LND also induces cellular reactive oxygen species through Complex II and has been reported to promote cell death by suppression of the pentose phosphate pathway, which resulted in inhibition of NADPH and glutathione generation. We conclude that MPC inhibition is the most sensitive anti-tumour target for LND, with additional inhibitory effects on MCT-mediated l-lactic acid efflux, Complex II and glutamine/glutamate oxidation.

(13)C MRS and LC-MS Flux Analysis of Tumor Intermediary Metabolism

We present the first validated metabolic network model for analysis of flux through key pathways of tumor intermediary metabolism, including glycolysis, the oxidative and non-oxidative arms of the pentose pyrophosphate shunt, the TCA cycle as well as its anaplerotic pathways, pyruvate-malate shuttling, glutaminolysis, and fatty acid biosynthesis and oxidation. The model that is called Bonded Cumomer Analysis for application to (13)C magnetic resonance spectroscopy ((13)C MRS) data and Fragmented Cumomer Analysis for mass spectrometric data is a refined and efficient form of isotopomer analysis that can readily be expanded to incorporate glycogen, phospholipid, and other pathways thereby encompassing all the key pathways of tumor intermediary metabolism. Validation was achieved by demonstrating agreement of experimental measurements of the metabolic rates of oxygen consumption, glucose consumption, lactate production, and glutamate pool size with independent measurements of these parameters in cultured human DB-1 melanoma cells. These cumomer models have been applied to studies of DB-1 melanoma and DLCL2 human diffuse large B-cell lymphoma cells in culture and as xenografts in nude mice at 9.4 T. The latter studies demonstrate the potential translation of these methods to in situ studies of human tumor metabolism by MRS with stable (13)C isotopically labeled substrates on instruments operating at high magnetic fields (≥7 T). The melanoma studies indicate that this tumor line obtains 51% of its ATP by mitochondrial metabolism and 49% by glycolytic metabolism under both euglycemic (5 mM glucose) and hyperglycemic conditions (26 mM glucose). While a high level of glutamine uptake is detected corresponding to ~50% of TCA cycle flux under hyperglycemic conditions, and ~100% of TCA cycle flux under euglycemic conditions, glutaminolysis flux and its contributions to ATP synthesis were very small. Studies of human lymphoma cells demonstrated that inhibition of mammalian target of rapamycin (mTOR) signaling produced changes in flux through the glycolytic, pentose shunt, and TCA cycle pathways that were evident within 8 h of treatment and increased at 24 and 48 h. Lactate was demonstrated to be a suitable biomarker of mTOR inhibition that could readily be monitored by (1)H MRS and perhaps also by FDG-PET and hyperpolarized (13)C MRS methods.

Refined Analysis of Brain Energy Metabolism Using In Vivo Dynamic Enrichment of 13C Multiplets

Carbon-13 nuclear magnetic resonance spectroscopy in combination with the infusion of ¹³C-labeled precursors is a unique approach to study in vivo brain energy metabolism. Incorporating the maximum information available from in vivo localized ¹³C spectra is of importance to get broader knowledge on cerebral metabolic pathways. Metabolic rates can be quantitatively determined from the rate of ¹³C incorporation into amino acid neurotransmitters such as glutamate and glutamine using suitable mathematical models. The time course of multiplets arising from ¹³C-¹³C coupling between adjacent carbon atoms was expected to provide additional information for metabolic modeling leading to potential improvements in the estimation of metabolic parameters.

The aim of the present study was to extend two-compartment neuronal/glial modeling to include dynamics of ¹³C isotopomers available from fine structure multiplets in ¹³C spectra of glutamate and glutamine measured in vivo in rats brain at 14.1 T, termed bonded cumomer approach. Incorporating the labeling time courses of ¹³C multiplets of glutamate and glutamine resulted in elevated precision of the estimated fluxes in rat brain as well as reduced correlations between them.

Bonded Cumomer Analysis of Human Melanoma Metabolism Monitored by 13C NMR Spectroscopy of Perfused Tumor Cells

A network model for the determination of tumor metabolic fluxes from (13)C NMR kinetic isotopomer data has been developed and validated with perfused human DB-1 melanoma cells carrying the BRAF V600E mutation, which promotes oxidative metabolism. The model generated in the bonded cumomer formalism describes key pathways of tumor intermediary metabolism and yields dynamic curves for positional isotopic enrichment and spin-spin multiplets. Cells attached to microcarrier beads were perfused with 26 mm [1,6-(13)C2]glucose under normoxic conditions at 37 °C and monitored by (13)C NMR spectroscopy. Excellent agreement between model-predicted and experimentally measured values of the rates of oxygen and glucose consumption, lactate production, and glutamate pool size validated the model. ATP production by glycolytic and oxidative metabolism were compared under hyperglycemic normoxic conditions; 51% of the energy came from oxidative phosphorylation and 49% came from glycolysis. Even though the rate of glutamine uptake was ∼ 50% of the tricarboxylic acid cycle flux, the rate of ATP production from glutamine was essentially zero (no glutaminolysis). De novo fatty acid production was ∼ 6% of the tricarboxylic acid cycle flux. The oxidative pentose phosphate pathway flux was 3.6% of glycolysis, and three non-oxidative pentose phosphate pathway exchange fluxes were calculated. Mass spectrometry was then used to compare fluxes through various pathways under hyperglycemic (26 mm) and euglycemic (5 mm) conditions. Under euglycemic conditions glutamine uptake doubled, but ATP production from glutamine did not significantly change. A new parameter measuring the Warburg effect (the ratio of lactate production flux to pyruvate influx through the mitochondrial pyruvate carrier) was calculated to be 21, close to upper limit of oxidative metabolism.

Metabolic Modeling of Dynamic (13)C NMR Isotopomer Data in the Brain In Vivo: Fast Screening of Metabolic Models Using Automated Generation of Differential Equations

Most current brain metabolic models are not capable of taking into account the dynamic isotopomer information available from fine structure multiplets in (13)C spectra, due to the difficulty of implementing such models. Here we present a new approach that allows automatic implementation of multi-compartment metabolic models capable of fitting any number of (13)C isotopomer curves in the brain. The new automated approach also makes it possible to quickly modify and test new models to best describe the experimental data. We demonstrate the power of the new approach by testing the effect of adding separate pyruvate pools in astrocytes and neurons, and adding a vesicular neuronal glutamate pool. Including both changes reduced the global fit residual by half and pointed to dilution of label prior to entry into the astrocytic TCA cycle as the main source of glutamine dilution. The glutamate-glutamine cycle rate was particularly sensitive to changes in the model.

Inhibition of Mitochondrial Complex II by the Anticancer Agent Lonidamine

The antitumor agent lonidamine (LND; 1-(2,4-dichlorobenzyl)-1H-indazole-3-carboxylic acid) is known to interfere with energy-yielding processes in cancer cells. However, the effect of LND on central energy metabolism has never been fully characterized. In this study, we report that a significant amount of succinate is accumulated in LND-treated cells. LND inhibits the formation of fumarate and malate and suppresses succinate-induced respiration of isolated mitochondria. Utilizing biochemical assays, we determined that LND inhibits the succinate-ubiquinone reductase activity of respiratory complex II without fully blocking succinate dehydrogenase activity. LND also induces cellular reactive oxygen species through complex II, which reduced the viability of the DB-1 melanoma cell line. The ability of LND to promote cell death was potentiated by its suppression of the pentose phosphate pathway, which resulted in inhibition of NADPH and glutathione generation. Using stable isotope tracers in combination with isotopologue analysis, we showed that LND increased glutaminolysis but decreased reductive carboxylation of glutamine-derived α-ketoglutarate. Our findings on the previously uncharacterized effects of LND may provide potential combinational therapeutic approaches for targeting cancer metabolism.

How to measure metabolic fluxes: a taxonomic guide for (13)C fluxomics

Metabolic reaction rates (fluxes) contribute fundamentally to our understanding of metabolic phenotypes and mechanisms of cellular regulation. Stable isotope-based fluxomics integrates experimental data with biochemical networks and mathematical modeling to 'measure' the in vivo fluxes within an organism that are not directly observable. In recent years, (13)C fluxomics has evolved into a technology with great experimental, analytical, and mathematical diversity. This review aims at establishing a unified taxonomy by means of which the various fluxomics methods can be compared to each other. By linking the developed modeling approaches to recent studies, their challenges and opportunities are put into perspective. The proposed classification serves as a guide for scientific 'travelers' who are striving to resolve research questions with the currently available (13)C fluxomics toolset.

Characterization of the usage of the serine metabolic network in human cancer

The serine, glycine, one-carbon (SGOC) metabolic network is implicated in cancer pathogenesis, but its general functions are unknown. We carried out a computational reconstruction of the SGOC network and then characterized its expression across thousands of cancer tissues. Pathways including methylation and redox metabolism exhibited heterogeneous expression indicating a strong context dependency of their usage in tumors. From an analysis of coexpression, simultaneous up- or downregulation of nucleotide synthesis, NADPH, and glutathione synthesis was found to be a common occurrence in all cancers. Finally, we developed a method to trace the metabolic fate of serine using stable isotopes, high-resolution mass spectrometry, and a mathematical model. Although the expression of single genes didn't appear indicative of flux, the collective expression of several genes in a given pathway allowed for successful flux prediction. Altogether, these findings identify expansive and heterogeneous functions for the SGOC metabolic network in human cancer.

(13)C NMR spectroscopy applications to brain energy metabolism

(13)C nuclear magnetic resonance (NMR) spectroscopy is the method of choice for studying brain metabolism. Indeed, the most convincing data obtained to decipher metabolic exchanges between neurons and astrocytes have been obtained using this technique, thus illustrating its power. It may be difficult for non-specialists, however, to grasp thefull implication of data presented in articles written by spectroscopists. The aim of the review is, therefore, to provide a fundamental understanding of this topic to facilitate the non-specialists in their reading of this literature. In the first part of this review, we present the metabolic fate of (13)C-labeled substrates in the brain in a detailed way, including an overview of some general neurochemical principles. We also address and compare the various spectroscopic strategies that can be used to study brain metabolism. Then, we provide an overview of the (13)C NMR experiments performed to analyze both intracellular and intercellular metabolic fluxes. More particularly, the role of lactate as a potential energy substrate for neurons is discussed in the light of (13)C NMR data. Finally, new perspectives and applications offered by (13)C hyperpolarization are described.

In vivo quantification of neuro-glial metabolism and glial glutamate concentration using 1H-[13C] MRS at 14.1T

Astrocytes have recently become a major center of interest in neurochemistry with the discoveries on their major role in brain energy metabolism. An interesting way to probe this glial contribution is given by in vivo (13) C NMR spectroscopy coupled with the infusion labeled glial-specific substrate, such as acetate. In this study, we infused alpha-chloralose anesthetized rats with [2-(13) C]acetate and followed the dynamics of the fractional enrichment (FE) in the positions C4 and C3 of glutamate and glutamine with high sensitivity, using (1) H-[(13) C] magnetic resonance spectroscopy (MRS) at 14.1T. Applying a two-compartment mathematical model to the measured time courses yielded a glial tricarboxylic acid (TCA) cycle rate (Vg ) of 0.27 ± 0.02 μmol/g/min and a glutamatergic neurotransmission rate (VNT ) of 0.15 ± 0.01 μmol/g/min. Glial oxidative ATP metabolism thus accounts for 38% of total oxidative metabolism measured by NMR. Pyruvate carboxylase (VPC ) was 0.09 ± 0.01 μmol/g/min, corresponding to 37% of the glial glutamine synthesis rate. The glial and neuronal transmitochondrial fluxes (Vx (g) and Vx (n) ) were of the same order of magnitude as the respective TCA cycle fluxes. In addition, we estimated a glial glutamate pool size of 0.6 ± 0.1 μmol/g. The effect of spectral data quality on the fluxes estimates was analyzed by Monte Carlo simulations. In this (13) C-acetate labeling study, we propose a refined two-compartment analysis of brain energy metabolism based on (13) C turnover curves of acetate, glutamate and glutamine measured with state of the art in vivo dynamic MRS at high magnetic field in rats, enabling a deeper understanding of the specific role of glial cells in brain oxidative metabolism. In addition, the robustness of the metabolic fluxes determination relative to MRS data quality was carefully studied.

Computational approaches for understanding energy metabolism

There has been a surge of interest in understanding the regulation of metabolic networks involved in disease in recent years. Quantitative models are increasingly being used to interrogate the metabolic pathways that are contained within this complex disease biology. At the core of this effort is the mathematical modeling of central carbon metabolism involving glycolysis and the citric acid cycle (referred to as energy metabolism). Here, we discuss several approaches used to quantitatively model metabolic pathways relating to energy metabolism and discuss their formalisms, successes, and limitations.

Metabolic Flux and Compartmentation Analysis in the Brain In vivo

Through significant developments and progresses in the last two decades, in vivo localized nuclear magnetic resonance spectroscopy (MRS) became a method of choice to probe brain metabolic pathways in a non-invasive way. Beside the measurement of the total concentration of more than 20 metabolites, (1)H MRS can be used to quantify the dynamics of substrate transport across the blood-brain barrier by varying the plasma substrate level. On the other hand, (13)C MRS with the infusion of (13)C-enriched substrates enables the characterization of brain oxidative metabolism and neurotransmission by incorporation of (13)C in the different carbon positions of amino acid neurotransmitters. The quantitative determination of the biochemical reactions involved in these processes requires the use of appropriate metabolic models, whose level of details is strongly related to the amount of data accessible with in vivo MRS. In the present work, we present the different steps involved in the elaboration of a mathematical model of a given brain metabolic process and its application to the experimental data in order to extract quantitative brain metabolic rates. We review the recent advances in the localized measurement of brain glucose transport and compartmentalized brain energy metabolism, and how these reveal mechanistic details on glial support to glutamatergic and GABAergic neurons.

Skeletal muscle mitochondrial uncoupling in a murine cancer cachexia model

Approximately half of all cancer patients present with cachexia, a condition in which disease-associated metabolic changes lead to a severe loss of skeletal muscle mass. Working toward an integrated and mechanistic view of cancer cachexia, we investigated the hypothesis that cancer promotes mitochondrial uncoupling in skeletal muscle. We subjected mice to in vivo phosphorous-31 nuclear magnetic resonance (³¹P NMR) spectroscopy and subjected murine skeletal muscle samples to gas chromatography/mass spectrometry (GC/MS). The mice used in both experiments were Lewis lung carcinoma models of cancer cachexia. A novel ‘fragmented mass isotopomer’ approach was used in our dynamic analysis of ¹³C mass isotopomer data. Our ³¹P NMR and GC/MS results indicated that the adenosine triphosphate (ATP) synthesis rate and tricarboxylic acid (TCA) cycle flux were reduced by 49% and 22%, respectively, in the cancer-bearing mice (p<0.008; t-test vs. controls). The ratio of ATP synthesis rate to the TCA cycle flux (an index of mitochondrial coupling) was reduced by 32% in the cancer-bearing mice (p=0.036; t-test vs. controls). Genomic analysis revealed aberrant expression levels for key regulatory genes and transmission electron microscopy (TEM) revealed ultrastructural abnormalities in the muscle fiber, consistent with the presence of abnormal, giant mitochondria. Taken together, these data suggest that mitochondrial uncoupling occurs in cancer cachexia and thus point to the mitochondria as a potential pharmaceutical target for the treatment of cachexia. These findings may prove relevant to elucidating the mechanisms underlying skeletal muscle wasting observed in other chronic diseases, as well as in aging.

Mitochondria-targeted antioxidant promotes recovery of skeletal muscle mitochondrial function after burn trauma assessed by in vivo 31P nuclear magnetic resonance and electron paramagnetic resonance spectroscopy

Burn injury causes a major systemic catabolic response that is associated with mitochondrial dysfunction in skeletal muscle. We investigated the effects of the mitochondria-targeted peptide antioxidant Szeto-Schiller 31 (SS-31) on skeletal muscle in a mouse burn model using in vivo phosphorus-31 nuclear magnetic resonance ((31)P NMR) spectroscopy to noninvasively measure high-energy phosphate levels; mitochondrial aconitase activity measurements that directly correlate with TCA cycle flux, as measured by gas chromatography mass spectrometry (GC-MS); and electron paramagnetic resonance (EPR) to assess oxidative stress. At 6 h postburn, the oxidative ATP synthesis rate was increased 5-fold in burned mice given a single dose of SS-31 relative to untreated burned mice (P=0.002). Furthermore, SS-31 administration in burned animals decreased mitochondrial aconitase activity back to control levels. EPR revealed a recovery in redox status of the SS-31-treated burn group compared to the untreated burn group (P<0.05). Our multidisciplinary convergent results suggest that SS-31 promotes recovery of mitochondrial function after burn injury by increasing ATP synthesis rate, improving mitochondrial redox status, and restoring mitochondrial coupling. These findings suggest use of noninvasive in vivo NMR and complementary EPR offers an approach to monitor the effectiveness of mitochondrial protective agents in alleviating burn injury symptoms.

Modeling the glutamate-glutamine neurotransmitter cycle

Glutamate is the principal excitatory neurotransmitter in brain. Although it is rapidly synthesized from glucose in neural tissues the biochemical processes for replenishing the neurotransmitter glutamate after glutamate release involve the glutamate–glutamine cycle. Numerous in vivo¹³C magnetic resonance spectroscopy (MRS) experiments since 1994 by different laboratories have consistently concluded: (1) the glutamate–glutamine cycle is a major metabolic pathway with a flux rate substantially greater than those suggested by early studies of cell cultures and brain slices; (2) the glutamate–glutamine cycle is coupled to a large portion of the total energy demand of brain function. The dual roles of glutamate as the principal neurotransmitter in the CNS and as a key metabolite linking carbon and nitrogen metabolism make it possible to probe glutamate neurotransmitter cycling using MRS by measuring the labeling kinetics of glutamate and glutamine. At the same time, comparing to non-amino acid neurotransmitters, the added complexity makes it more challenging to quantitatively separate neurotransmission events from metabolism. Over the past few years our understanding of the neuronal-astroglial two-compartment metabolic model of the glutamate–glutamine cycle has been greatly advanced. In particular, the importance of isotopic dilution of glutamine in determining the glutamate–glutamine cycling rate using [1−¹³C] or [1,6-¹³C₂] glucose has been demonstrated and reproduced by different laboratories. In this article, recent developments in the two-compartment modeling of the glutamate–glutamine cycle are reviewed. In particular, the effects of isotopic dilution of glutamine on various labeling strategies for determining the glutamate–glutamine cycling rate are analyzed. Experimental strategies for measuring the glutamate–glutamine cycling flux that are insensitive to isotopic dilution of glutamine are also suggested.