Jacob L.
Gorenflos López
ab,
Peter
Schmieder
a,
Kristin
Kemnitz-Hassanin
a,
Hatice Ceyda
Asikoglu
ac,
Arif
Celik
ab,
Christian E.
Stieger
ab,
Dorothea
Fiedler
ab,
Stephan
Hinderlich
*c and
Christian P. R.
Hackenberger
*ab
aLeibniz-Institut für Molekulare Pharmakologie, Robert-Roessle-Strasse 10, 13125 Berlin, Germany. E-mail: hackenbe@fmp-berlin.de
bHumboldt Universität zu Berlin, Department Chemie, Brook-Taylor-Strasse 2, 12489, Germany
cBerliner Hochschule für Technik, Department Life Sciences & Technology, Seestrasse 64, 13347 Berlin, Germany
First published on 21st February 2023
Sialic acids are part of the outermost component of the glycocalyx of all vertebrates; as such, they are fundamental markers in physiological and pathological processes. In this study, we introduce a real-time assay to monitor individual enzymatic steps of sialic acid biosynthesis, either with recombinant enzymes, in particular using UDP-N-acetylglucosamine 2-epimerase (GNE) or N-acetylmannosamine kinase (MNK), or in cytosolic rat liver extract. Using state-of-the-art NMR techniques, we are able to follow the characteristic signal of the N-acetyl methyl group, which displays different chemical shifts for the biosynthesis intermediates UDP-N-acetylglucosamine, N-acetylmannosamine (and its 6-phosphate) and N-acetylneuraminic acid (and its 9-phosphate). Pseudo 2- and 3-D NMR demonstrated that in rat liver cytosolic extract, the phosphorylation reaction of MNK is exclusive for N-acetylmannosamine generated by GNE. Thus, we speculate that phosphorylation of this sugar from other sources (e.g. external application to cells) or N-acetylmannosamine derivatives often applied in metabolic glycoengineering is not conducted by MNK but by a yet unknown sugar kinase. Competition experiments with the most prevalent neutral carbohydrates demonstrated that of these, only N-acetylglucosamine slowed N-acetylmannosamine phosphorylation kinetics, suggesting an N-acetylglucosamine-preferring kinase as the acting enzyme.
In vertebrates, de novo sialic acid biosynthesis proceeds in four consecutive steps in the cytosol. First, uridine diphosphate N-acetylglucosamine (UDP-GlcNAc), produced from fructose 6-phosphate by the hexosamine biosynthetic pathway,3 is epimerised to produce N-acetylmannosamine (ManNAc). Subsequently, using ATP, it is phosphorylated to ManNAc-6-phosphate (ManNAc-6P). Both steps are catalysed by the bifunctional UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE/MNK).4 Then, sialic acid synthase (SAS) catalyses the condensation reaction of phosphoenolpyruvic acid (PEP) and ManNAc-6P to Neu5Ac 9-phosphate (Neu5Ac-9P). Neu5Ac-9P is dephosphorylated by sialic acid phosphatase (SAP), the only non-essential enzyme of this pathway.5 Neu5Ac is activated to CMP-Neu5Ac by CMP-Neu5Ac synthetase (CMP-SA) in the nucleus. Afterwards, it is transported by the specific nucleotide sugar transporter solute carrier 35A1 (SLC35A1) into the Golgi,6 where it serves as a substrate for sialyltransferases (Scheme 1).
Metabolic glycan engineering (MGE) is a methodology used to introduce specific chemical entities into cells via glycobiological pathways.7 Because of the terminal position of sialic acids and the promiscuity of the involved enzymes, many MGE approaches use ManNAc derivatives, which are directly metabolised to sialic acids.7,8 The most widely used of these is N-azidoacetylmannosamine (ManNAz), usually in its per-acetylated form to ensure cell permeability, which is metabolised to N-azidoacetylneuraminic acid (Neu5Az) and can be used as an easily accessible bioorthogonal handle targeting the azide on cell surface glycoconjugates9,10 and soluble glycoproteins.11,12 The first metabolic step of this and related N-acyl-modified ManNAc derivatives, after their acetyl esters have been cleaved off by unspecific esterases, is the phosphorylation to ManNAc-6P derivatives. First, it was speculated that MNK catalyses this step. Still, low conversion efficiencies of the recombinant kinase13 and studies with GNE/MNK knock-down cells14 indicated that another kinase might be responsible for this enzymatic step. Furthermore, experiments with external ManNAc (we differentiate between biosynthetic ManNAc produced in cells and external ManNAc, which is added artificially) rescues the cell surface sialic acid content of lymphoma BJA-B cells with a GNE/MNK knock-out.15
Analysis of sialic acid biosynthesis metabolites is a challenging approach. Recombinant and purified enzymes often allow colourimetric assays or optical tests for the characterisation of single steps of the pathway, including recording of enzyme kinetics.16,17 Analysis of the metabolites in a cellular context was initially performed by radioactive compounds, which were separated by paper or thin-layer chromatography.18 Improved methods allowed high-performance liquid chromatography (HPLC) assays for labelled Neu5Ac and its derivatives19 or HPLC-mass spectrometry (MS), thereby often employing 13C or 1H labelled compounds.20,21 However, all of these methods are endpoint assays, which so far did not allow real-time monitoring of one metabolite or even several compounds in parallel. Therefore, our aim in the current paper was to develop a real-time assay to monitor key steps in sialic acid biosynthesis in a cell-free cytosolic system.
Real-time NMR (RT-NMR) monitoring of metabolic pathways in physiological environments has received considerable attention in recent years. Due to the non-invasive, non-disruptive and quantitative nature of RT-NMR, with NMR signal intensities reflecting the absolute concentrations of components,22–26 we considered RT-NMR to be the ideal method to follow different biosynthetic intermediates of the sialic acid biosynthesis pathway, such as UDP-GlcNAc, ManNAc, ManNAc-6P and Neu5Ac, even in the complex environment of a tissue homogenate. In particular, we envisioned using this approach to obtain more information on the metabolic pathway flux and on the phosphorylation of biosynthetic vs. external ManNAc, due to its importance for MGE.
The α to β ratios of ManNAc and ManNAc-6P were determined by integrating the peaks to roughly be 1.2:1 (ESI† Fig. 3). A combination of HMBC and HMQC revealed, in tune with the literature, the downfield peak of ManNAc to be the β-anomer (ESI† Fig. 4).29 The methyl groups of GlcNAc and GalNAc only present one peak for both anomers. Surprisingly, both anomers of ManNAc-6P show the same chemical shift difference of 0.03 ppm for the N-Ac methyl in 1H-NMR compared to ManNAc. This was not the case for the GlcNAc and GlcNAc-6P pair (Fig. 1B). We speculated that the axial orientation of the N-Ac methyl in ManNAc and ManNAc-6P positions the methyl group to interact with the anomeric hydroxyl, thus causing peak separation between α- and β-anomer. In the case of the chemical shift between the phosphorylated and the unphosphorylated ManNAc, we considered that the phosphate interacts with the axial amide either sterically, through a salt bridge or a hydrogen bond (ESI† Fig. 5).
In parallel to measuring the chemical shifts of the methyl group of N-acetyl hexosamines, 1H-NMR also enabled the monitoring of many other proton-containing molecules, such as ATP, ADP, AMP, by the aromatic protons of the nucleotides, and PEP and pyruvate by their vinyl and methyl protons, respectively. With the help of standards, characteristic peaks for these molecules could be attributed to their origin, which provided further highly valuable information on coenzymes and substrates from our assay (ESI† Fig. 6).
At 10 μg mL−1 (215 nM) GNE converted 400 μM UDP-GlcNAc in about 5 h to ManNAc (half-life 1.1 h, ESI† Table 1). 1H-NMR spectra were recorded in hourly intervals. The stacking of these spectra enabled the quantification of the epimerization of UDP-GlcNAc to ManNAc. The UDP-GlcNAc signal was normalised to the first measurement. The ManNAc signals were normalised to the averaged signal of the plateau after 7 h (Fig. 2A). This and all following reactions were analysed using a one-phase decay model (ESI† Fig. 10).
Fig. 2 Topographic view (oblique mode) of the pseudo-2D NMR; the monitored reactions and the graphs resulting from the integration of the raw data are displayed here in triplicates. (A) Real-time 1H-NMR of the epimerisation of UDP-GlcNAc (400 μM) to ManNAc by recombinant GNE (10 μg mL−1) – nonlinear fit data in ESI† Table 1. (B) Real-time 1H-NMR of the phosphorylation of ManNAc to ManNAc-6p by recombinant MNK (1 μg mL−1) with ATP (5 mM). (C) Real-time 1H-NMR of the conversion of UDP-GlcNAc (400 μM) to ManNAc-6p by recombinant GNE/MNK (10 μg mL−1) with ATP (5 mM) – nonlinear fit data in ESI† Table 2. |
Recombinant MNK proved to be significantly more active than recombinant GNE, 1 μg mL−1 (29 nM) MNK was determined to be a more suitable concentration for the 1 h measuring intervals (ESI† Fig. 11). Real-time 1H-NMR spectroscopy enabled monitoring of the phosphorylation of ManNAc to ManNAc-6P by MNK, but due to poor peak separation the ManNAc and ManNAc-6P peaks could not be integrated separately, which made separate quantification of educt and product unreliable (Fig. 2B).
To demonstrate that we could monitor epimerisation and phosphorylation a cocktail of both enzymes at 10 μg mL−1 of each, GNE (215 nM) and MNK (290 nM), was employed to convert UDP-GlcNAc to ManNAc-6P (Fig. 2C). To mimic the production of Neu5Ac by the sialic acid synthase, commercially available sialic acid aldolase from E. coli (1 mg mL−1; 29 μM) and pyruvate (200 mM) were used. ManNAc was generated in situ using GNE and UDP-GlcNAc. Real-time 1H-NMR spectroscopy proved to be efficient at monitoring the conversion reactions (ESI† Fig. 12A). Neu5Ac production was corroborated by the characteristic signal of the axial ring protons at C3. Part of the ring proton signal overlapped with unidentified compounds (ESI† Fig. 12B).
Fig. 3 (A) Topographic view (oblique mode) of the pseudo 2D NMR spectrum of the phosphorylation of ManNAc (400 μM) to ManNAc-6p by recombinant MNK (1 μg mL−1) with ATP (5 mM). The individual plane (contour mode) demonstrates shows the lack of peak separation of ManNAc and ManNAc-6P in 1H-NMR. (B) Schematic representation of SOFAST-HMQC coupling of ManNAc/-6P-2-13C and the resulting 2D spectra. These 2D spectra can be stacked into a cube, introducing two pseudo planes, each containing separate peaks; (C) pseudo planes of β-ManNAc/-6P described in (B). The ManNAc and ManNAc-6P peaks can be integrated separately. (D) 31P-HMBC spectrum of ManNAc-6P. (E) 31P-HMBC spectra can be stacked into cubes, allowing the introduction of pseudo planes displayed here. (F) Integration of the pseudo planes of the SOFAST-HMQC and 31P-HMBC – non-linear fit data in ESI† Table 4. |
The phosphorylation of ManNAc-2-13C by MNK was monitored using SOFAST-HMQC,32 which was run using non uniform sampling (NUS)33 measuring only 25% of all data points and reducing the spectral width in the 13C-dimension to 0.6 ppm around the ManNAc/-6P signals, which differed on the two utilized spectrometers (22.0 ± 0.3 ppm). The resulting spectra were stacked in a cube, and the peak separation was achieved through the introduction of pseudo-planes (Fig. 3B), which enabled the separation of the β-ManNAc and the β-ManNAc-6P peaks and the integration (Fig. 3C). The same was possible with the α-ManNAc signal (ESI† Fig. 18). All following experiments were analysed using the β-ManNAc and the β-ManNAc-6P peaks because the α signals in part overlapped with other N-acetyl methyl signals, like that of GlcNAc (Fig. 1A/B). To further validate this method, a 31P-HMBC approach was designed, which was also facilitated by NUS of 25% of the data points. ManNAc-6P has a characteristic 31P-HMBC signal (Fig. 3D) that can be stacked like the SOFAST-HMQC spectra. Again time-resolved analysis became possible by introducing a pseudo plane and the integration (Fig. 3E). Multiple measurements were run in sequence. First, the 1H-NMR was run, then the SOFAST-HMQC and last, the 31P-HMBC. All three experiments in a row had been optimised to take less than 5 min, which included moving the probe in the autosampler. While the 1H-NMR provided an initial overview of the data (Fig. 3A), which can also be used to monitor parallel reactions of ATP and PEP, the SOFAST-HMQC offered a powerful method to monitor ManNAc phosphorylation by recombinant enzymes. For recombinant MNK, the 31P-HMBC method could validate the rate constants from the SOFAST-HMQC experiments (Fig. 3F). 31P-HMBC therefore proved to be a robust method to monitor the phosphorylation of monosaccharides by recombinant enzymes.
First, we performed a GNE assay with RLCE. In the absence of ATP and PEP, UDP-GlcNAc was converted in roughly equal parts to ManNAc and GlcNAc. This is likely due to the activity of GlcNAc 2-epimerase.35,36 However, we could not rule out that UDP-GlcNAc is converted to GlcNAc by a non-physiological activity of GNE/MNK in the absence of ATP, or an independent UDP-GlcNAc hydrolysing enzyme. Because of overlapping signals with GlcNAc, we assumed an α to β ratio of 1.2:1 for ManNAc based on our initial experiments (ESI† Table 5). 27% of the initial UDP-GlcNAc signal was not regained by the sum of ManNAc and GlcNAc signals (Fig. 4A). It is likely, that both amino sugars were further converted by catabolic pathways, including release of the N-acetyl group. This cleaved group is not detected by our NMR assay, leading to loss of signal.
Fig. 4 Topographic view (oblique mode) of the pseudo-2D NMR; the monitored reactions and the graphs resulting from the integration of the raw data are displayed here in triplicates. (A) Real-time 1H-NMR of the conversion of UDP-GlcNAc (400 μM) in RLCE (1.6 mg mL−1) without further cofactors. The percentages are given relative to the initial UDP-GlcNAc signal – nonlinear fit data in ESI† Table 5. (B) Real-time 1H-NMR of the phosphorylation of ManNAc-2-13C (400 μM) in RLCE (1.6 mg mL−1) – nonlinear fit data in ESI† Table 6. (C) Real-time 1H-NMR of the conversion of UDP-GlcNAc (400 μM) to ManNAc-6P with ATP (5 mM) and PEP (25 mM) in RLCE (1.6 mg mL−1) – nonlinear fit data in ESI† Table 7. |
External ManNAc (400 μM) was phosphorylated to ManNAc-6P in RLCE (1.6 mg mL−1). The half-life of the reaction was 0.84 h (ESI† Table 6). After 3 h of reaction, ManNAc was completely converted to ManNAc-6P (Fig. 4B). Although significant amounts of PEP were present, ManNAc-6P was not converted to Neu5Ac. It is notable that under the same conditions the half-life of the phosphorylation of external ManNAc to ManNAc-6P (0.84 h) was more than an order of magnitude higher than the conversion of UDP-GlcNAc to ManNAc-6P 18.83 h (ESI† Table 7), indicating the activity of another kinase than MNK phosphorylating ManNAc.
In the presence of ATP and PEP, RLCE converts UDP-GlcNAc to ManNAc-6P (Fig. 4C). Though, the half-life of UDP-GlcNAc increases by an order of magnitude from 1.66 h (ESI† Table 5) to 18.83 h (ESI† Table 7). Under these conditions, the production of ManNAc was faster than the production of ManNAc-6P (ESI† Fig. 22).
Commercially available ManNAc-6P was converted with PEP to Neu5Ac/-9P in RLCE (6 mg mL−1). Due to the unavailability of Neu5Ac-9P standards, we could not differentiate whether the detected signal stemmed from the phosphorylated or the dephosphorylated Neu5Ac. Because of overlapping signals with Neu5Ac, we again could not measure the α-anomer signals of ManNAc-6P. Therefore, we calculated Neu5Ac concentrations assuming an α to β ratio of 1.2:1 (ESI† Fig. 23 and Table 8). Of all the recorded reactions of the sialic acid biosynthesis pathway in RLCE, the synthase reaction was the slowest – RLCE concentration had to be increased almost 4-fold to obtain a system in which significant concentrations of substrate were converted to product. Under these conditions, the half-life of ManNAc-6P was 6.46 h (ESI† Table 8). GNE/MNK and SAS are expressed in comparable levels in RLCE (ESI† Fig. 20). Therefore, specific activity of SAS must be significantly lower than the ones of GNE/MNK. Former studies revealed that in RLCE both specific activities of GNE/MNK4 are more than one order of magnitude higher than the specific activity of SAS.37 In conclusion, SAS could be suggested as the bottleneck for sialic acid biosynthesis in RLCE. This is in agreement with the essential role of SAS in sialic acid biosynthesis,5 which is also true for glycoengineering approaches to establish the pathway in sialic acid-deficient cells from insects38 or plants.39
Taken together, we could demonstrate that our RT-NMR approach using 1H-NMR and SOFAST-HMQC experiments could be employed to measure substrate conversions of GNE/MNK and SAS in RLCE in the complex environment of RLCE. Furthermore, the flexibility of the assay to measure unpredicted side products (e.g. GlcNAc) was demonstrated.
Fig. 5 (A) Simultaneous conversion of UDP-GlcNAc (400 μM) and ManNAc (400 μM) in RLCE (1.6 mg mL−1). The decreasing UDP-GlcNAc signals were monitored and the conversion of ManNAc to ManNAc-6P. The respective systems with either only UDP-GlcNAc or only ManNAc were used as control experiments, to reference the conversion velocities – nonlinear fit data in and ESI† Table 10. This experiment was conducted with a different RLCE than all other experiments. (B) Highly abundant monosaccharides ordered by structural similarity to ManNAc. (C) Competition experiments of the conversion of ManNAc (200 μM) to ManNAc-6P with other carbohydrates. The influence of GlcNAc was assessed at increasing ratios (1:2 to 1:10). Other relevant carbohydrates were only used at the ratio of 1:10. The decreasing ManNAc signal and the increasing ManNAc-6P signals were modelled with a one-phase decay model (nonlinear fit data in ESI† Table 11) and the resulting half-life times were plotted. The significance is given relative to the experiment without competing carbohydrate; data are represented as means ± SEM (n = 45; 15 time points per experiment conducted in triplicates) analysed by multiple t-test (P-value: 0.1234 (ns), 0.0332 (*), 0.0021 (**), 0.0002 (***) and <0.0001). |
We assumed that the responsible kinase for external ManNAc phosphorylation uses highly abundant and structurally related monosaccharides as substrates. These encompass primarily the N-acetyl hexosamines GlcNAc and GalNAc, and the hexoses mannose (Man), glucose (Glc) and galactose (Gal) (Fig. 5B). Only GlcNAc had a significant effect on the half-life of ManNAc-6P production from ManNAc, which increased from 1.20 h to nearly 3.76 h with a 10-fold excess of GlcNAc over ManNAc (Fig. 5C and ESI† Table 12) by use of a decreased RCLE concentration to 1 mg mL−1. This makes a kinase with a primary activity for GlcNAc a prime candidate to phosphorylate external ManNAc, instead of a kinase with a primary activity for ManNAc phosphorylation. Previously, significant MNK activity was already observed for partially purified N-acetylglucosamine kinase (NAGK) from rat liver by Allen and Walker.40 This observation was confirmed by homogenous purification of rat liver NAGK; MNK activity of NAGK could clearly be distinguished from the MNK activity of GNE/MNK by chromatographic separation.41 Our data were further underlined by human cell lines with low or completely absent GNE/MNK expression, where MNK activity could almost completely assigned to NAGK.15 Treatment of GNE myopathy, a hereditary inclusion body myopathy caused by mutations in GNE/MNK,42 with ManNAc, is currently in a phase 2 clinical study. Understanding the nature of the kinase responsible for the first step of metabolisation will enhance the knowledge with regard to failure and success of this trials.43,44
Based on the plate-reader experiments, two concentrations were tested for the GlcNAc phosphorylation by NAGK in NMR experiments. Whereas 10 μg mL−1 NAGK did not allow quantification of enzyme activity because the conversion was finished before the first measurement, 0.1 μg mL−1 NAGK showed complete phosphorylation of GlcNAc after 6 h (Fig. 6A). The phosphorylation of ManNAc (400 μM) by 10 μg mL−1 MNK with ATP (5 mM) was finished after 48 min. At 1 μg mL−1, complete transformation was reached at 16 h (Fig. 6B). These data confirm the activities of the two recombinant kinases observed in the plate-reader assay.
Fig. 6 (A and B) Real-time 31P-HMBC-NMR of the phosphorylation (ATP 5 mM) of GlcNAc and ManNAc (400 μM) by NAGK and MNK at different concentrations - nonlinear fit data in ESI† Table 13 and 14. (C and D) Real-time 31P-HMBC-NMR of the phosphorylation of different carbohydrates (400 μM) by NAGK and MNK with ATP (5 mM) - nonlinear fit data in ESI† Tables 15 and 16. |
GlcNAc-6P standard revealed a similar signal pattern in a 31P-HMBC as ManNAc-6P (ESI† Fig. 26). From this, we extrapolated that all hexoses phosphorylated at the hydroxy groups at carbon 6 show similar patterns in 31P-HMBC NMR spectra. As such, we used the 31P-HMBC real-time NMR assay to do studies over 18 h on the activity of NAGK and MNK with different substrates. At 10 μg mL−1 NAGK phosphorylates ManNAc, but also Glc and Man with similar rates within the detection limits of this assay. These data are in agreement with the identification of NAGK as “extrahepatic glucokinase”,46 suggesting that NAGK is an enzyme with high substrate promiscuity, at least in in vitro assays, with high Km values for ManNAc (Km = 0.95 mM) and Glc (Km = 600 mM).47 On the other hand, NAGK reveals no activity with Gal and GalNAc (Fig. 6C), revealing a physiologically relevant selection of substrates. MNK does not appear to be promiscuous over longer incubation periods (Fig. 6D), underlining its specificity for ManNAc as a substrate.
In the context of MGE the question is still up to debate, which kinase phosphorylates the ManNAc derivatives. Earlier studies had demonstrated the ability of GNE/MNK knock down cells to still metabolize ManNAz,13 the most commonly used ManNAc derivative for MGE.7 Previous studies had further shown that recombinant rat GNE/MNK could phosphorylate certain ManNAc derivatives with low efficiency in dependence of the length of their N-acyl side chains.13 We therefore investigated both kinases for their ability to phosphorylate ManNAz. Surprisingly, MNK showed no activity for this artificial sugar, whereas NAGK was able to phosphorylate ManNAz with the same activity as shown for the other alternative substrates (Fig. 6C and D), again pointing towards NAGK as a responsible kinase in MGE of N-acyl-mannosamine derivatives.
GNE myopathy is presumably caused by a lack of sialic acid production due to a malfunctioning GNE/MNK. As such, our study supports the clinical application of external ManNAc, because external ManNAc phosphorylation was shown functioning independently of GNE/MNK. Still, it also highlights a problem with this treatment: the deregulation of sialic acid biosynthesis. The utilization of external ManNAc circumvents the master regulator of sialic acid biosynthesis: GNE/MNK.48 Without the feedback inhibition of the end-product CMP-Neu5Ac sialic acid production is deregulated and could cause symptoms akin to Sialuria, a hereditary developmental disorder caused by mutations in the CMP-Neu5Ac binding pocket of GNE/MNK.49,50 Furthermore, we think that this deregulation of sialic acid biosynthesis is the reason for why MGE with ManNAc derivatives has been employed so successfully. A recent study showed that azide modifications of carbohydrates were hampering metabolisation rates.51,52 Nonetheless, in a HEK 293T cell model, 65% of the sialome after 24 h of treatment with Ac4ManNAz carried an azide.53 We speculate that the kinase responsible for the phosphorylation of external ManNAc, also phosphorylates ManNAc derivatives.
In summary, with this study we highlight the power of a real-time NMR assay, which, in contrast to other endpoint methods, can differentiate between forming monosaccharide species in situ. As such, we envision this method to be a singularly useful approach to study the conversion velocities of biologically highly relevant carbohydrates in complex environments to evaluate the performance of potential inhibitors and elucidate metabolic pathways.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2sc06986e |
This journal is © The Royal Society of Chemistry 2023 |