Crosstalk between AML and stromal cells triggers acetate secretion through the metabolic rewiring of stromal cells.

Acute myeloid leukaemia (AML) cells interact and modulate components of their surrounding microenvironment into their own benefit. Stromal cells have been shown to support AML survival and progression through various mechanisms. Nonetheless, whether AML cells could establish beneficial metabolic interactions with stromal cells is underexplored. By using a combination of human AML cell lines and AML patient samples together with mouse stromal cells and a MLL-AF9 mouse model, here we identify a novel metabolic crosstalk between AML and stromal cells where AML cells prompt stromal cells to secrete acetate for their own consumption to feed the tricarboxylic acid cycle (TCA) and lipid biosynthesis. By performing transcriptome analysis and tracer-based metabolic NMR analysis, we observe that stromal cells present a higher rate of glycolysis when co-cultured with AML cells. We also find that acetate in stromal cells is derived from pyruvate via chemical conversion under the influence of reactive oxygen species (ROS) following ROS transfer from AML to stromal cells via gap junctions. Overall, we present a unique metabolic communication between AML and stromal cells and propose two different molecular targets, ACSS2 and gap junctions, that could potentially be exploited for adjuvant therapy.

Quantitative Isotopomer Rates in Real-Time Metabolism of Cells Determined by NMR Methods.

Tracer-based metabolism is becoming increasingly important for studying metabolic mechanisms in cells. NMR spectroscopy offers several approaches to measure label incorporation in metabolites, including 13 C- and 1 H-detected spectra. The latter are generally more sensitive, but quantification depends on the proton-carbon 1 JCH coupling constant, which varies significantly between different metabolites. It is therefore not possible to have one experiment optimised for all metabolites, and quantification of 1 H-edited spectra such as HSQCs requires precise knowledge of coupling constants. Increasing interest in tracer-based and metabolic flux analysis requires robust analyses with reasonably small acquisition times. Herein, we compare 13 C-filtered and 13 C-edited methods for quantification and show the applicability of the methods for real-time NMR spectroscopy of cancer-cell metabolism, in which label incorporations are subject to constant flux. We find an approach using a double filter to be most suitable and sufficiently robust to reliably obtain 13 C incorporations from difference spectra. This is demonstrated for JJN3 multiple myeloma cells processing glucose over 24 h. The proposed method is equally well suited for calculating the level of label incorporation in labelled cell extracts in the context of metabolic flux analysis.

A framework for tracer-based metabolism in mammalian cells by NMR.

Metabolism changes extensively during the normal proliferation and differentiation of mammalian cells, and in cancer and inflammatory diseases. Since changes in the metabolic network reflect interactions between genetic, epigenetic and environmental changes, it is helpful to study the flow of label from isotopically labelled precursors into other metabolites rather than static metabolite levels. For this Nuclear Magnetic Resonance (NMR) spectroscopy is an attractive technique as it can quantify site-specific label incorporation. However, for applications using human cells and cell lines, the challenge is to optimize the process to maximize sensitivity and reproducibility. Here we present a new framework to analyze metabolism in mammalian cell lines and primary cells, covering the workflow from the preparation of cells to the acquisition and analysis of NMR spectra. We have applied this new approach in hematological and liver cancer cell lines and confirm the feasibility of tracer-based metabolism in primary liver cells.

PhenoMeNal: processing and analysis of metabolomics data in the cloud.

BACKGROUND: Metabolomics is the comprehensive study of a multitude of small molecules to gain insight into an organism's metabolism. The research field is dynamic and expanding with applications across biomedical, biotechnological, and many other applied biological domains. Its computationally intensive nature has driven requirements for open data formats, data repositories, and data analysis tools. However, the rapid progress has resulted in a mosaic of independent, and sometimes incompatible, analysis methods that are difficult to connect into a useful and complete data analysis solution. FINDINGS: PhenoMeNal (Phenome and Metabolome aNalysis) is an advanced and complete solution to set up Infrastructure-as-a-Service (IaaS) that brings workflow-oriented, interoperable metabolomics data analysis platforms into the cloud. PhenoMeNal seamlessly integrates a wide array of existing open-source tools that are tested and packaged as Docker containers through the project's continuous integration process and deployed based on a kubernetes orchestration framework. It also provides a number of standardized, automated, and published analysis workflows in the user interfaces Galaxy, Jupyter, Luigi, and Pachyderm. CONCLUSIONS: PhenoMeNal constitutes a keystone solution in cloud e-infrastructures available for metabolomics. PhenoMeNal is a unique and complete solution for setting up cloud e-infrastructures through easy-to-use web interfaces that can be scaled to any custom public and private cloud environment. By harmonizing and automating software installation and configuration and through ready-to-use scientific workflow user interfaces, PhenoMeNal has succeeded in providing scientists with workflow-driven, reproducible, and shareable metabolomics data analysis platforms that are interfaced through standard data formats, representative datasets, versioned, and have been tested for reproducibility and interoperability. The elastic implementation of PhenoMeNal further allows easy adaptation of the infrastructure to other application areas and 'omics research domains.

Metabolomic Evidence for a Field Effect in Histologically Normal and Metaplastic Tissues in Patients with Esophageal Adenocarcinoma.

Patients with Barrett's esophagus (BO) are at increased risk of developing esophageal adenocarcinoma (EAC). Most Barrett's patients, however, do not develop EAC, and there is a need for markers that can identify those most at risk. This study aimed to see if a metabolic signature associated with the development of EAC existed. For this, tissue extracts from patients with EAC, BO, and normal esophagus were analyzed using 1H nuclear magnetic resonance. Where possible, adjacent histologically normal tissues were sampled in those with EAC and BO. The study included 46 patients with EAC, 7 patients with BO, and 68 controls who underwent endoscopy for dyspeptic symptoms with normal appearances. Within the cancer cohort, 9 patients had nonneoplastic Barrett's adjacent to the cancer suitable for biopsy. It was possible to distinguish between histologically normal, BO, and EAC tissue in EAC patients [area under the receiver operator curve (AUROC) 1.00, 0.86, and 0.91] and between histologically benign BO in the presence and absence of EAC (AUROC 0.79). In both these cases, sample numbers limited the power of the models. Comparison of histologically normal tissue proximal to EAC versus that from controls (AUROC 1.00) suggests a strong field effect which may develop prior to overt EAC and hence be useful for identifying patients at high risk of developing EAC. Excellent sensitivity and specificity were found for this model to distinguish histologically normal squamous esophageal mucosa in EAC patients and healthy controls, with 8 metabolites being very significantly altered. This may have potential diagnostic value if a molecular signature can detect tissue from which neoplasms subsequently arise.

Malonate as a ROS product is associated with pyruvate carboxylase activity in acute myeloid leukaemia cells.

BACKGROUND: The role of anaplerotic nutrient entry into the Krebs cycle via pyruvate carboxylase has been the subject of increased scrutiny and in particular whether this is dysregulated in cancer. Here, we use a tracer-based NMR analysis involving high-resolution (1)H-(13)C-HSQC spectra to assess site-specific label incorporation into a range of metabolite pools, including malate, aspartate and glutamate in the acute myeloid leukaemia cell line K562. We also determine how this is affected following treatment with the redeployed drug combination of the lipid-regulating drug bezafibrate and medroxyprogesterone (BaP). RESULTS: Using the tracer-based approach, we assessed the contribution of pyruvate carboxylase (PC) vs. pyruvate dehydrogenase (PDH) activity in the derivation of Krebs cycle intermediates. Our data show that PC activity is indeed high in K562 cells. We also demonstrate a branched entry to the Krebs cycle of K562 cells with one branch running counterclockwise using PC-derived oxaloacetate and the other clockwise from the PDH activity. Finally, we show that the PC activity of K562 cells exclusively fuels the ROS-induced decarboxylation of oxaloacetate to malonate in response to BaP treatment; resulting in further Krebs cycle disruption via depletion of oxaloacetate and malonate-mediated inhibition of succinate dehydrogenase (SDH) resulting in a twofold reduction of fumarate. CONCLUSIONS: This study extends the interest in the PC activity in solid cancers to include leukaemias and further demonstrates the value of tracer-based NMR approaches in generating a more accurate picture of the flow of carbons and metabolites within the increasingly inappropriately named Krebs cycle. Moreover, our studies indicate that the PC activity in cancer cells can be exploited as an Achilles heel by using treatments, such as BaP, that elevate ROS production.

Tracer-Based Metabolic NMR-Based Flux Analysis in a Leukaemia Cell Line.

High levels of reactive oxygen species (ROS) have a profound impact on acute myeloid leukaemia cells and can be used to specifically target these cells with novel therapies. We have previously shown how the combination of two redeployed drugs, the contraceptive steroid medroxyprogesterone and the lipid-regulating drug bezafibrate exert anti-leukaemic effects by producing ROS. Here we report a 13C-tracer-based NMR metabolic study to understand how these drugs work in K562 leukaemia cells. Our study shows that [1,2-13C]glucose is incorporated into ribose sugars, indicating activity in oxidative and non-oxidative pentose phosphate pathways alongside lactate production. There is little label incorporation into the tricarboxylic acid cycle from glucose, but much greater incorporation arises from the use of [3-13C]glutamine. The combined medroxyprogesterone and bezafibrate treatment decreases label incorporation from both glucose and glutamine into α-ketoglutarate and increased that for succinate, which is consistent with ROS-mediated conversion of α-ketoglutarate to succinate. Most interestingly, this combined treatment drastically reduced the production of several pyrimidine synthesis intermediates.

Tracer-Based Metabolic NMR-Based Flux Analysis in a Leukaemia Cell Line.

High levels of reactive oxygen species (ROS) have a profound impact on acute myeloid leukaemia cells and can be used to specifically target these cells with novel therapies. We have previously shown how the combination of two redeployed drugs, the contraceptive steroid medroxyprogesterone and the lipid-regulating drug bezafibrate exert anti-leukaemic effects by producing ROS. Here we report a 13 C-tracer-based NMR metabolic study to understand how these drugs work in K562 leukaemia cells. Our study shows that [1,2-13 C]glucose is incorporated into ribose sugars, indicating activity in oxidative and non-oxidative pentose phosphate pathways alongside lactate production. There is little label incorporation into the tricarboxylic acid cycle from glucose, but much greater incorporation arises from the use of [3-13 C]glutamine. The combined medroxyprogesterone and bezafibrate treatment decreases label incorporation from both glucose and glutamine into α-ketoglutarate and increased that for succinate, which is consistent with ROS-mediated conversion of α-ketoglutarate to succinate. Most interestingly, this combined treatment drastically reduced the production of several pyrimidine synthesis intermediates.

Structural tightening and interdomain communication in the catalytic cycle of phosphoglycerate kinase.

Changes in amide-NH chemical shift and hydrogen exchange rates as phosphoglycerate kinase progresses through its catalytic cycle have been measured to assess whether they correlate with changes in hydrogen bonding within the protein. Four representative states were compared: the free enzyme, a product complex containing 3-phosphoglyceric acid (3PG), a substrate complex containing ADP and a transition-state analogue (TSA) complex containing a 3PG-AlF(4)(-)-ADP moiety. There are an overall increases in amide protection from hydrogen exchange when the protein binds the substrate and product ligands and an additional increase when the TSA complex is formed. This is consistent with stabilisation of the protein structure by ligand binding. However, there is no correlation between the chemical shift changes and the protection factor changes, indicating that the protection factor changes are not associated with an overall shortening of hydrogen bonds in the protected ground state, but rather can be ascribed to the properties of the high-energy, exchange-competent state. Therefore, an overall structural tightening mechanism is not supported by the data. Instead, we observed that some cooperativity is exhibited in the N-domain, such that within this domain the changes induced upon forming the TSA complex are an intensification of those induced by binding 3PG. Furthermore, chemical shift changes induced by 3PG binding extend through the interdomain region to the C-domain beta-sheet, highlighting a network of hydrogen bonds between the domains that suggests interdomain communication. Interdomain communication is also indicated by amide protection in one domain being significantly altered by binding of substrate to the other, even where no associated change in the structure of the substrate-free domain is indicated by chemical shifts. Hence, the communication between domains is also manifested in the accessibility of higher-energy, exchange-competent states. Overall, the data that are consistent with structural tightening relate to defined regions and are close to the 3PG binding site and in the hinge regions of 3-phosphoglycerate kinase.

Complete assignment of hyaluronan oligosaccharides up to hexasaccharides.

The glycosaminoglycan hyaluronan is involved in a diverse range of physiological and diseases processes and comprises repeated disaccharide units of N-acetyl-d-glucosamine (GlcNAc) and d-glucuronic acid (GlcA). A molecular description of the solution conformation of HA is required to account for this biology, which is best attained using nuclear magnetic resonance (NMR). NMR studies of the polymer, however, are frustrated by resonance overlap arising from the highly degenerate structure. In contrast, end-effects in oligosaccharides can produce some chemical shift dispersion, giving the possibility that their conformational properties can be measured and extrapolated to models of the polymer. We report the complete resolution and assignment of (1)H, (13)C and (15)N nuclei in hyaluronan oligosaccharides with seven different naturally occurring terminal rings. At 900MHz, all (1)H nuclei in the hexasaccharide GlcA-beta-(1-->3)-GlcNAc-beta-(1-->4)-GlcA-beta-(1-->3)-GlcNAc-beta-(1-->4)-GlcA-beta-(1-->3)-GlcNAc-OH were uniquely resolved and the two central rings were found to be a good model for the polymer environment. These assignments now allow resolved, unambiguous structural restraints to be acquired on this oligosaccharide and extrapolated to models for the solution conformation of the polymer.

NMR spectra of oligosaccharides at ultra-high field (900 MHz) have better resolution than expected due to favourable molecular tumbling.

Nuclear magnetic resonance (NMR) remains the most promising technique for acquiring atomic-resolution information in complex carbohydrates. Significant obstacles to the acquisition of such data are the poor chemical-shift dispersion and artifacts resultant from their degenerate chemical structures. The recent development of ultra-high-field NMR (at 900 MHz and beyond) gives new potential to overcome these problems, as we demonstrate on a hexasaccharide of the highly repetitive glycosaminoglycan hyaluronan. At 900 MHz, the expected increase in spectral dispersion due to higher resonance frequencies and reduction in strong coupling-associated distortions are observed. In addition, the fortuitous molecular tumbling rate of oligosaccharides results in longer T2-values that further significantly enhances resolution, an effect not available to proteins. Combined, the resolution enhancement can be as much as twofold relative to 600 MHz, allowing all 1H-resonances in the hexasaccharide to be unambiguously assigned using standard natural-abundance experiments. The use of ultra-high-field spectrometers is clearly advantageous and promises a new and exciting era in carbohydrate structural biology.

The denatured state under native conditions: a non-native-like collapsed state of N-PGK.

The guanidinium-denatured state of the N-domain of phosphoglycerate kinase (PGK) has been characterized using solution NMR. Rather than behaving as a homogenous ensemble of random coils, chemical shift changes for the majority of backbone amide resonances indicate that the denatured ensemble undergoes two definable equilibrium transitions upon titration with guanidinium, in addition to the major refolding event. (13)C and (15)N chemical shift changes indicate that both intermediary states have distinct helical character. At denaturant concentrations immediately above the mid-point of unfolding, size-exclusion chromatography shows N-PGK to have a compact, denatured form, suggesting that it forms a helical molten globule. Within this globule, the helices extend into some regions that become beta strands in the native state. This predisposition of the denatured state to extensive non-native-like conformation, illustrates that, rather than directing folding, conformational pre-organization in the denatured state can compete with the normal folding direction. The corresponding reduction in control of the direction of folding as proteins become larger, could thus constitute a restriction on the size of protein domains.

Effects of domain dissection on the folding and stability of the 43 kDa protein PGK probed by NMR.

The characterization of early folding intermediates is key to understanding the protein folding process. Previous studies of the N-domain of phosphoglycerate kinase (PGK) from Bacillus stearothermophilus combined equilibrium amide exchange data with a kinetic model derived from stopped-flow kinetics. Together, these implied the rapid formation of an intermediate with extensive native-like hydrogen bonding. However, there was an absence of protection in the region proximal to the C-domain in the intact protein. We now report data for the intact PGK molecule, which at 394 residues constitutes a major extension to the protein size for which such data can be acquired. The methods utilised to achieve the backbone assignment are described in detail, including a semi-automated protocol based on a simulated annealing Monte Carlo technique. A substantial increase in the stability of the contact region is observed, allowing protection to be inferred on both faces of the beta-sheet in the intermediate. Thus, the entire N-domain acts concertedly in the formation of the kinetic refolding intermediate rather than there existing a distinct local folding nucleus.