==== Front Nat Commun Nat Commun Nature Communications 2041-1723 Nature Publishing Group UK London 20035 10.1038/s41467-020-20035-0 Article Reverse C-glycosidase reaction provides C-nucleotide building blocks of xenobiotic nucleic acids Pfeiffer Martin 12 http://orcid.org/0000-0002-5030-2643Nidetzky Bernd bernd.nidetzky@tugraz.at 12 1 grid.410413.30000 0001 2294 748XInstitute of Biotechnology and Biochemical Engineering, Graz University of Technology, NAWI Graz, Petersgasse 12, 8010 Graz, Austria 2 grid.432147.70000 0004 0591 4434Austrian Centre of Industrial Biotechnology (acib), Petersgasse 14, 8010 Graz, Austria 8 12 2020 8 12 2020 2020 11 627020 5 2020 9 11 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.C-Analogues of the canonical N-nucleosides have considerable importance in medicinal chemistry and are promising building blocks of xenobiotic nucleic acids (XNA) in synthetic biology. Although well established for synthesis of N-nucleosides, biocatalytic methods are lacking in C-nucleoside synthetic chemistry. Here, we identify pseudouridine monophosphate C-glycosidase for selective 5-β-C-glycosylation of uracil and derivatives thereof from pentose 5-phosphate (d-ribose, 2-deoxy-d-ribose, d-arabinose, d-xylose) substrates. Substrate requirements of the enzymatic reaction are consistent with a Mannich-like addition between the pyrimidine nucleobase and the iminium intermediate of enzyme (Lys166) and open-chain pentose 5-phosphate. β-Elimination of the lysine and stereoselective ring closure give the product. We demonstrate phosphorylation-glycosylation cascade reactions for efficient, one-pot synthesis of C-nucleoside phosphates (yield: 33 – 94%) from unprotected sugar and nucleobase. We show incorporation of the enzymatically synthesized C-nucleotide triphosphates into nucleic acids by RNA polymerase. Collectively, these findings implement biocatalytic methodology for C-nucleotide synthesis which can facilitate XNA engineering for synthetic biology applications. C-nucleosides are analogues of the canonical N-nucleosides and, despite their synthetic value, biocatalysis has not targeted them yet. Here, the authors report a pseudouridine monophosphate C-glycosidase enzyme for selective 5-β-C-glycosylation of uracil and its derivatives from pentose 5- phosphate substrates. Subject terms BiocatalysisGlycosidesEnzyme mechanismsHydrolasesSynthetic chemistry methodologyhttps://doi.org/10.13039/501100002428Austrian Science Fund (Fonds zur Förderung der Wissenschaftlichen Forschung)DK Molecular Enzymology W901Nidetzky Bernd issue-copyright-statement© The Author(s) 2020 ==== Body Introduction Most natural and synthetic nucleosides (e.g., uridine, U; Fig. 1a) have their sugar and nucleobase parts connected via a hemiaminal ether (N-glycosidic) linkage. Nucleoside analogs replacing the canonical N- by a rare-to-nature C-glycosidic linkage resist the common routes of biological degradation by N-glycoside hydrolases and phosphorylases1–4. Since in vivo stability is crucial for compound efficacy and applicability in therapy, synthetic C-nucleosides have considerable importance in medicinal chemistry. Nucleoside analogs generally represent promising lead structures in drug design for antiviral, antibacterial and antitumor activity5–9. Synthetic C-nucleosides, such as BCX4430 (Galidesivir) and GS-6620 (Supplementary Fig. 1), can be remarkably effective against filovirus infections and hepatitis C virus, respectively7,10–13. Originally developed for use against filoviruses, GS-5734 (Remdesivir; Supplementary Fig. 1) is now also evaluated against infections of the SARS-CoV-2 virus10,13,14. Other C-nucleosides, such as showdomycin, formycin, pseudouridimycin and minimycin (Supplementary Fig. 1), constitute an emerging class of microbial natural products that exhibit unusual structures and diverse biological activities, including antibiotic efficacy15–21.Fig. 1 Proposed reaction mechanism of YeiN and the enzymatic phosphorylation-glycosylation cascade used in this study. The N-nucleoside uridine (U) (a) compared to the C-nucleoside pseudouridine (Ψ) (b) and reverse reaction of YeiN for synthesis of Ψ 5’-phosphate (ΨMP) from d-Rib5P (1) and Ura (12) (c). In panel c, B is the base residue from the enzyme. The key catalytic steps are the following: formation of a covalent iminium ion intermediate between the enzyme (Lys166) and the open-chain d-Rib5P, assisted by Glu31 (I); Mannich-like reaction for C-C coupling (II); β-elimination of the amine of Lys166 (III); cyclization through an oxa-Michael addition-like reaction (IV). For further discussion of the enzyme mechanism, see the text under the Results and the Discussion. d Phosphorylation-glycosylation cascades for one-pot synthesis of ΨMP and Ψ 5’-triphosphate (ΨTP) from unphosphorylated d-Rib are shown. Dephosphorylation of ΨMP yields Ψ. RbsK d-Rib-5-kinase, PK pyruvate kinase, CIP calf intestine phosphatase, NDK nucleoside diphosphate kinase. PEP phosphoenolpyruvate, PYR pyruvate. C-nucleosides are furthermore important as natural or synthetic building blocks of nucleic acids, RNA in particular. The U isomer pseudouridine (Ψ; Fig. 1b) occurs naturally in all types of RNA (tRNA, rRNA, snRNA, snoRNA, scaRNA, and mRNA)22–26. The Ψ formation is catalyzed by Ψ-synthase, which post-translationally isomerizes U to Ψ 27,28. Owing to additional hydrogen bonding from the nucleobase N1 in Ψ compared to U (Fig. 1a, b), positional substitution of U with Ψ can increase RNA stability by locally enhancing the macromolecular rigidity23. Chemically synthesized Ψ and derivatives thereof are important building blocks for study of C-nucleoside biology as well as for numerous synthetic biology applications based on xenobiotic nucleic acids (XNA)29–36. For example, the uniform substitution of U by Ψ can render synthetic mRNA non-immunogenic and was shown to enhance the translation efficiency30–32. Additionally, Ψ-modified mRNA can boost the gene-cutting efficiency in CRISPR-based gene editing37,38. To advance C-nucleosides for their emerging applications in XNA engineering or in drug design, the target compounds must be made available through efficient and structurally precise synthetic methods. Chemical routes to C-nucleosides typically follow the synthetic rationale of (1) reconstructing the heteroaryl nucleobase on a C1′-functionalized sugar or (2) coupling the prepared heteroaryl with a suitable glycosyl reagent, such as a hydroxy group-protected glycal or pentonolactone5,39–43. However, despite much progress from research over decades, a purely chemical approach to defined C-nucleosides remains challenging to establish44–50. A synthesis that exploits the high stereo- and regiocontrol of enzymatic reactions for installing the C-glycosidic linkage would be highly desirable. Natural enzymes of C-glycoside biosynthesis, such as tRNA pseudouridine synthase27,28,51,52 and C-glycosynthase15,18–21,53,54 have complex substrate requirements, leading to perceived limitations on their applicability. Therefore, unlike N-nucleosides for which N-nucleoside phosphorylases present excellent tools of glycoside synthesis55–60, biocatalytic methods are lacking for C-nucleoside synthetic chemistry. Here, we identify pseudouridine monophosphate (ΨMP) C-glycosidase for C-nucleoside synthesis. The enzyme’s natural role is the recycling of Ψ via degradation of its 5′-phosphate into uracil (Ura(12)) and d-ribose 5-phosphate (d-Rib5P(1)), as shown in Supplementary Fig. 261,62. Despite early evidence suggesting applicability63–65, the reverse C-glycosidase reaction (Fig. 1c) is unexplored synthetically. Multiplicity of enzymes in C-nucleoside turnover and overlapped enzyme naming used previously complicate clear assignment of enzyme type to function. For this study, we use the biochemically and structurally well-characterized ΨMP C-glycosidase from E. coli (YeiN). We show that YeiN catalyzes the selective 5-β-C-glycosylation of Ura (12) and derivatives thereof from several pentose 5-phosphate substrates, including the 5-phosphate of d-ribose, 2-deoxy-d-ribose, 2-deoxy-2-fluoro-d-ribose, d-arabinose or d-xylose. The ΨMP and its corresponding sugar or nucleobase variants are obtained in a single-step, highly efficient transformation. Remarkably, although formally a reversion of hydrolysis, the overall enzymatic conversion benefits from the favorable thermodynamics of its central catalytic step, a carbon-carbon bond-forming (Mannich-like) addition. It thus proceeds with excellent synthetic yields (≥90%), even in water. Finally, we show incorporation of the enzymatically synthesized C-nucleotide triphosphates into RNA. Collectively, therefore, our study implements biocatalytic methodology for C-nucleotide synthesis, which can facilitate XNA engineering for synthetic biology applications and establish new routes to natural product C-glycosides. Results The substrate scope of YeiN in reverse C-glycosidase reaction Figure 2 shows the substrates used. Besides the native substrates d-Rib5P (1) and Ura (12), we tested a series of sugar phosphates (in combination with Ura, 12) and nucleobase analogs (in combination with d-Rib5P, 1) as alternative substrates (Fig. 2) of C-glycosylation by YeiN. Based on product formation detected with HPLC, 2-deoxy-d-Rib5P (d-dRib5P, 2), d-xylose 5-phosphate (d-Xyl5P, 3), d-arabinose 5-phosphate (d-Ara5P, 4) and 2-deoxy-2-fluoro-d-ribose 5-phosphate (2F-d-Rib5PP, 5) are active. l-Ribose 5-phosphate (l-Rib5P, 6), d-glucose 6-phosphate (d-Glc6P, 7), d-mannose 6-phosphate (d-Man6P, 8), d-fructose 6-phosphate (d-Fru6P, 9), d-ribulose 5-phosphate (d-Rul5P, 10) and d-erythrose 4-phosphate (d-Ery4P, 11) are inactive. We also considered the 5-phosphates of 2-amino-2-deoxy- and 3-amino-3-deoxy-d-ribose, l-arabinose, d-lyxose, and l-xylose (Supplementary Fig. 3), but enzymatic phosphorylation of these pentoses from ATP proved not possible. The nucleobase analogs used for reaction with YeiN feature substitutions at each Ura carbon (Fig. 2). 3-Methyl-Ura (13), 6-amino-Ura (14), 4-thio-Ura (15) and 2-thio-Ura (16) are active. Absence of activity with 1-methyl-Ura (17), 5-fluoro-Ura (18) and 6-trifluoro-methyl-Ura (19) can be mechanistically relevant, as discussed below. Additionally, thymine (5-methyl-Ura, 20), 6-aza-Ura (21), 6-Cl-Ura (22), cytosine (23), 2,4-thio-Ura (24), 5-ethynyl-Ura (25), maleimide (26), adenine (27), guanine (28), 4-pyrimidone (29) and pyrimidine (30) were evaluated but found to be inactive. Combinatorial assessment of all sugar 5-phosphates active with Ura (12) against the full panel of Ura analogs used with d-Rib5P (1) confirms the substrate scope of YeiN (Supplementary Table 1).Fig. 2 Substrate scope of YeiN. Substrates tested for C-glycosylation by YeiN. The compounds framed in blue are active while those framed in magenta are not. We obtained kinetic parameters for the substrates (Fig. 2 and Supplementary Table 2). Compared with d-Rib5P (1) and Ura (12), catalytic efficiencies (kcat/KM) for alternative substrates are decreased between ~10 to ~103-fold, reflecting substantial (≥5-fold) changes in both kcat (decrease) and KM (increase). Characterization of the YeiN reaction The substrate scope of YeiN is rationalized in terms of structural and mechanistic requirements of the enzymatic reaction. Molecular-docking results show that, with the terminal phosphate group accommodated in the enzyme binding pocket (Fig. 3a), a five-carbon chain is needed to precisely position the substrates carbonyl for iminium ion formation with the ε-amino group of Lys166. Docking of d-Rib5P (1) results in a potentially catalytic complex, with the ε-amino group of Lys166 positioned (3.3 Å) for nucleophilic attack at the C1 carbonyl carbon. Thr139 (2.6 Å) and Glu31 (3.2 Å) are placed suitably to stabilize negative charge developing on the C1 oxygen atom during carbinol intermediate formation (Fig. 3a, c). The carbinol intermediate may be stabilized by Thr131 and the main chain carbonyl oxygen of Gly132 (Fig. 3c). Glu31 is furthermore positioned to promote expulsion of water from the carbinol to give the iminium ion (protonated Schiff base). In terms of structural arrangement and proposed catalytic function, Thr139 and Glu31 represent a clear analogy to aldolase active sites66–72. In terms of sequence and three-dimensional structure, however, YeiN shows hardly any similarity with well-studied aldolases, such as deoxy-ribose phosphate aldolase67,71, d-fructose bisphosphate aldolase69, and 2-dehydro-3-deoxy-phosphogluconate aldolase70,72. Mechanistically, YeiN differs from these aldolases in that it requires the iminium ion as an electrophile for C–C coupling (Fig. 1c). The aldolases convert the iminium ion to an enamine (carbanion iminium by resonance) so that it can function as a nucleophile in their aldol reactions66–70,72.Fig. 3 The substrate specificity of YeiN analyzed with molecular docking. Docking poses of a d-Rib5P (1, cyan) and b Ura (12, white). The enzyme bound Mn2+ is shown as a violet sphere. Distances are shown in Å. c Suggested mechanism of formation of the key covalent iminium ion intermediate that serves as electrophile in the subsequent Mannich-like addition, as shown in Fig. 1c. Contrary to d-Rib5P (1), the docking poses for d-Ery4P (11), d-Rul5P (10), and d-Fru6P (9) have Thr139 and Glu31 not properly aligned to the carbonyl oxygen, geometrically (d-Ery4P, 11; Supplementary Fig. 4a) or in terms of interaction distance (d-Rul5P, 10; d-Fru6P, 9; Supplementary Fig. 4b, d, respectively). Additionally, the poses of d-Ery4P (11) and d-Rul5P (10) suggest loss of the stabilizing interaction with Gly132. The d-Glc6P (7) pose (Supplementary Fig. 4c) shows an enlarged distance between the ε-amino of Lys166 and the C1 carbon. Lack of activity of these sugar phosphates is thus plausibly explained. Docking of Ura (12) into the structure of the covalent enzyme-d-Rib5P adduct (iminium ion form) places the reactive sugar C1′ and nucleobase C5 at a distance (3.4 Å) and orientation (91° angle between C5-C1′ and pyrimidine ring plane; Fig. 3b) well suited for carbon-carbon bond formation. Mechanistically, the catalytic reaction on the enzyme resembles a Mannich-like addition, for it involves the addition of a nucleophile (the cyclic enamine of the nucleobase) to an iminium ion (the activated aldehyde of d-Rib5P, 1) (Figs. 1c and 3c). Except for positioning, there seems to be limited catalytic facilitation from the enzyme to the agency of the Ura (12) C5 as a nucleophile of the reaction. The crystal structure of a covalent intermediate between YeiN and ΨMP shows just a single hydrogen bond between the C2 oxygen atom of Ura and Asn289 (Supplementary Fig. 5)62. The docking pose of Ura (12) suggests additional polar contacts (Fig. 3b) that can contribute to the precise positioning of the C5 of Ura above the C1 of the covalently bound d-Rib5P (1). Further docking analysis summarized in Supplementary Fig. 6 shows that lack of reactivity with d-Rib5P (1) can reflect steric restraint on the productive binding of certain nucleobases, such as 6-trifluoromethyl-Ura (19), 6-Cl-Ura (22) and 1-methyl-Ura (17) (Supplementary Fig. 6a, b, d). Other Ura (12) analogs, including 5-fluoro-Ura (18) as well as 5-methyl-Ura (20) and 6-aza-Ura (21) (Supplementary Fig. 6e, f, c) are accommodated in the enzyme binding pocket almost identically as the parent Ura (12). Positioned properly for catalysis, these Ura (12) analogs seem to be rendered inactive intrinsically. The proposed path for the catalytic reaction of YeiN (Fig. 1c) suggests mechanistic rationale for why the particular substitutions are incompatible with the enzyme activity. The 1-methyl substituent precludes formation of the mechanistically important cyclic imine tautomer in the catalytic step of C–C coupling (Fig. 1c). Owing to electronic effect, the 6-aza-Ura (21) is an unlikely nucleophile for reaction with the iminium ion of the covalent enzyme intermediate. In the subsequent catalytic step, as shown in Fig. 1c, the Mannich adduct undergoes β-elimination of the amine of the lysine (Fig. 1c). Ura (12) analogs bearing substituents at the C5 are chemically incompetent to release the lysine analogously. In the last step of the synthetic reaction, cyclization to a β-C-pentofuranosyl ring occurs and this proceed chemically through a kind of oxa-Michael addition (Fig. 1c). The cyclization could happen on the enzyme, in which case it is likely to be stereospecific, or in solution. Anomeric stereoselectivity can be a distinguishing feature of enzymes that catalyze C-ribosylation. The C-glycosynthase AlnA produces a mixture of α- (16%) and β-glycosidic (84%) forms during C-ribosylation of the naphthoquinone prealnumycin from d-Rib5P (1), suggesting ribosyl ring closure without enzymatic assistance (Supplementary Fig. 7)54. Here, we used in situ proton NMR in 2H2O solvent to monitor ΨMP formation from, and degradation to, d-Rib5P (1) and Ura (12). The time-resolved NMR spectra for the synthesis direction (Fig. 4a–c; 14 mM Ura, 12 and 8 mM d-Rib5P, 1) show that only a single anomer (β-C-riboside) is formed. They also show that solvent deuterium is incorporated at the Ura (12) C5 at a rate (0.07 mM/min) about 3-fold lower than the rate of conversion of Ura (12) into ΨMP (0.22 mM/min). This suggests that the deuterium uptake does not happen immediately during ΨMP synthesis but occurs only later when ΨMP and Ura (12) equilibrate dynamically. Suitable controls show that deuterium uptake into Ura (12) necessitates ΨMP formation (Supplementary Fig. 8). Inactive sugar phosphates (d-Ery4P, 11; d-Rul5P, 10; d-Glc6P, 7; and d-Fru6P, 9) do not promote deuterium uptake into the C5 of Ura (12). Analyzing the hydrolysis of ΨMP (Fig. 4d–f), we confirm that the Ura (12) released from ΨMP incorporates solvent deuterium at its C5. The hydrolysis reaction approaches equilibrium after about 100 min and an equilibrium constant of 8.9 × 10−5 M can be calculated from the data (0.9 mM Ura, 12 and d-Rib5P, 1; 9.1 mM ΨMP).Fig. 4 Enzymatic synthesis and hydrolysis of ΨMP monitored by 1H-NMR. a Scheme of enzymatic synthesis of ΨMP. Solvent deuterium is incorporated into Ura (12) C5 due to enzymatic hydrolysis of ΨMP. b 1H-NMR traces of ΨMP synthesis from 14 mM Ura (12) and 8 mM d-Rib5P (1), using 3 µM YeiN. c Time course of ΨMP synthesis as determined from the NMR spectra (n = 1). d Scheme of solvent deuterium incorporation into Ura C5 due to the enzymatic hydrolysis of ΨMP e 1H-NMR traces of hydrolysis of ΨMP (10 mM) by 3 µM YeiN. f Time course of ΨMP hydrolysis as determined from the NMR spectra (n = 1). Relative amounts of β-d-Rib5P and α-d-Rib5P are displayed on the secondary axis. The NMR analysis (Fig. 4e, f) additionally shows that at all times, the d-Rib5P (1) released from ΨMP involves α-furanose (33%) and β-furanose (67%) forms at anomeric equilibrium. Plausible interpretation is that the open-chain form of d-Rib5P is released into solution where the furanose structures are formed spontaneously. However, the mutarotation of d-Rib5P (1)73 occurs at the limits of temporal resolution of the experiment, leaving the conformation of the released d-Rib5P (1) undecided. Phosphorylation-glycosylation cascade for the synthesis of C-nucleoside monophosphates We considered that limitation on the synthetic utility of the YeiN-catalyzed formation of C-nucleoside monophosphates can arise from the requirement for a pentose 5-phosphate as the substrate. The unphosphorylated pentose, however, is an attractive starting material for synthesis. We therefore designed an enzymatic cascade reaction (Fig. 1d) in which pentose phosphorylation74 and C-glycosylation of the nucleobase are telescoped into a one-pot transformation. We show efficient enzymatic synthesis of ΨMP from d-Rib and Ura (12) and demonstrate applicability of the biocatalytic system to the synthesis of all sugar and nucleobase variants of the C-glycoside that were shown in Fig. 2 to be accessible via the YeiN reaction. As shown in Supplementary Fig. 9 and summarized in Table 1, the yield of synthesis ranges between useful and excellent (33-95%), and product (10–60 mg) can be recovered generally with high purity from small-scale reactions (15 mL). Product structures are confirmed by NMR (Supplementary Fig. 10–29) and agree well with published spectra75.Table 1 Enzymatic syntheses with yields and product purities indicated. n.d. no product formation detected. n.i. not isolated. *NMR spectra not recorded. Purity based on HPLC. §Based on Ura (12) converted. #Based on pentose converted. i5-(2-deoxy-2-fluoro-β-d-ribofuranosyl 5-phosphate)uracil, ii5-(β-d-arabinofuranosyl 5-phosphate)uracil, iii5-(β-d-xylofuranosyl 5-phosphate)uracil. Dephosphorylation and phosphorylation for structural diversification of the C-nucleoside triphosphates We considered that chemical biology and medicinal chemistry applications can require the unphosphorylated C-nucleoside or the C-nucleoside triphosphate. We therefore show enzymatic dephosphorylation of ΨMP with non-specific phosphatase to obtain the C-nucleoside in excellent yield (Table 1 and Supplementary Fig. 30). Dephosphorylation can be analogously applied to the other C-glycoside monophosphates. To prepare the C-nucleoside triphosphate (e.g., ΨTP), we set up a cascade reaction in which the monophosphate is phosphorylated by nucleotide monophosphate kinase and pyruvate kinase in the presence of ATP and PEP (Fig. 1d). Interestingly, the addition of a nucleotide diphosphate kinase did not improve efficiency of the overall phosphorylation (Supplementary Fig. 31)76. We show synthesis of ΨTP, 2-deoxy-ΨTP (dΨTP) and 6-amino ΨTP in a one-pot reaction starting from unphosphorylated sugar and Ura (12) (Fig. 5 and Supplementary Fig. 32). Products were recovered in mg amounts and in suitable purity (≥85%) for characterization and further studies.Fig. 5 Enzymatic cascade synthesis of ΨTP and dΨTP. Panels a and b show ΨTP and dΨTP, respectively. Pentoses (100 mM) were phosphorylated applying RbsK (15 µM for d-Rib; 30 µM for d-dRib) and pyruvate kinase (4 µM) using 100 mM PEP and 2 mM ATP as phosphate donors (see Supplementary Fig. 33). ΨMP synthesis was initiated by addition of 15 mM Ura, 1 mM MnCl2, YeiN (1 µM for d-Rib, 3 µM for d-dRib) and 10 mM HEPES buffer, resulting in a concentration of 25 mM pentose 5-phosphate (time point zero a, b. The arrow indicates the start of the phosphorylation reaction by addition of 15 mM PEP, 15 µM CMPK, and 1 mM ATP, resulting in a volume increase of 33% (n = 1). C-Nucleoside triphosphates for enzymatic synthesis of xenobiotic RNA Using the YeiN-encoding gene as the template, we examined amplification of the corresponding RNA and DNA using ΨTP, 6-amino-ΨTP and dΨTP in place of, or together with, the native UTP and dTTP, respectively. Applying T7 RNA polymerase as shown in Fig. 6a, we succeeded in synthesizing RNA built solely from Ψ instead of U.Fig. 6 Polymerase-catalyzed amplification of RNA and DNA. Panels a and b show RNA and DNA, respectively. RNA was transcribed from YeiN-encoding template DNA (n = 1). a In vitro translation was performed using standard conditions, 30 U of T7 RNA polymerase, 700 ng template DNA, transcription buffer (40 mM Tris-HCl, pH 7.9, at 25 °C; 6 mM MgCl2, 10 mM DTT, 10 mM NaCl and 2 mM spermidine) and 2 mM of each NTP (ATP, GTP, UTP, CTP). Modified RNA was synthesized either by addition of 2 mM of ΨTP to the NTPs (NTPs +  ΨTP) or by fully replacing UTP by 2 mM ΨTP (NTPs + ΨTP − UTP). b The yein gene was amplified from 50 ng expression vector (pet15_yein) using 1.25 U of Dream-Taq DNA polymerase, 0.2 mM of each dNTP (dATP, dCTP, dGTP, dTTP) and 0.5 µM of forward and reverse primer. To incorporate dΨ into the DNA either 0.2 mM of ΨTP were additionally added to dNTP mix (dNTPs + dΨTP) or TTP was fully replaced by dΨTP (dNTPs + YTP − TTP). Further details are described in “Methods”. Densitometric analysis of gel bands suggested that RNA polymerase transcription efficiency (expressed as gel band intensity) was similar when ΨTP was used to substitute UTP. Unexpectedly, the band intensity decreased to about half when ΨTP and UTP were added at the same time. Overloading of the RNA polymerase with NTP substrate may have caused inhibition, perhaps due to limitation in the concentration of Mg2+ used. 6-Amino-ΨTP and dΨTP were not accepted by the used RNA polymerase. Applying DreamTaq-polymerase that is known for its relaxed dNTP substrate specificity, we observed that dΨTP cannot substitute for dTTP (Fig. 6b). However, when used together with dTTP, the dΨTP appeared to inhibit the amplification slightly (∼25%). Discussion We here demonstrate the reverse reaction of ΨMP C-glycosidase for efficient biocatalytic synthesis of C-nucleoside monophosphates. Besides the natural ΨMP, a panel of sugar and nucleobase analogs of the C-nucleoside monophosphate are synthesized with excellent regiocontrol and stereoselectivity. The substrate requirements of the enzymatic reaction are consistent with the mechanistic proposal (Fig. 1c) of a Mannich-like C-C addition between the iminium ion of an enzyme bound open-chain pentose 5-phosphate intermediate and the pyrimidine nucleobase. β-Elimination of the enzyme residue and stereospecific cyclization via oxa-Michael-like addition give the C-nucleoside monophosphate product. Use of the iminium ion intermediate as an electrophile of the reaction distinguishes YeiN from well-studied aldolases that form essentially the same covalent intermediate but convert it to an enamine nucleophile for aldol coupling. Chemical studies of C-glycoside synthesis have shown that the condensation product is favored thermodynamically, even in water77,78, and the synthetic use of the chemical reaction was demonstrated73–77. Equilibrium on the YeiN reaction far on the side of the C-nucleoside monophosphate, which would seem unusual for a nominal glycosidase reaction, is thus explained. Enzymatic C–C coupling reactions are synthetically important in the carbohydrate sciences 65,78–86 and the current study adds new application-relevant feature to them. The development of enzymatic phosphorylation-glycosylation cascade reactions enable expedient, one-pot biocatalytic syntheses of the C-nucleoside 5′-phosphates from unphosphorylated pentoses as the substrates (Fig. 1d). We show that the C-nucleoside 5′-phosphates are readily dephosphorylated to give the C-nucleoside and that they can also be converted into the corresponding 5′-triphosphates. The non-phosphorylated C-nucleoside and its 5′-triphosphate are promising substrates for chemical and polymerase-catalyzed XNA synthesis, respectively. Phosphorylation-glycosylation cascades extend the synthetic utility of the ΨMP C-glycosidase reaction, enabling C-nucleoside phosphates to be obtained from unprotected sugar and nucleobase. The ΨMP C-glycosidases are related evolutionary to a class of C-glycosynthases, represented by the well-characterized AlnA from Streptomyces sp. CM020, that catalyze β-C-riboside formation from d-Rib5P (1) and different naphthoquinones (Supplementary Fig. 7)53,54,87. The herein developed idea of modular phosphorylation-glycosylation cascades could be analogously applied to AlnA-type C-glycosynthases, to also harness these enzymes’ catalytic proficiency and substrate promiscuity for practical synthesis. Methods Materials ATP ( > 95% purity) was from Roth (Karlsruhe, Germany). Deuterium oxide (99.96%) was from Euriso-Top (Saint-Aubin Cedex, France). Calf intestine phosphatase (CIP) was from New England Biolabs (Frankfurt am Main, Germany). Rabbit muscle pyruvate kinase (PK), adenosine triphosphate (ATP), d-ribose (d-Rib), d-2-deoxy-ribose (d-dRib), d-xylose (d-Xyl), d-arabinose (d-Ara), d-mannose (d-Man), d-glucosamine (d-GlcN), N-acetyl-d-glucosamine (d-GlcNAc), d-galactose (d-Gal), d-galactosamine (d-GalN), d-lyxose (d-Lyx), l-arabinose (l-Ara), l-xylose (l-Xyl), d-ribose-5-phosphate (d-Rib5P(1)), d-fructose 6-phosphate d-Fru6P(9), d-mannose 6-phosphate (d-Man6P (8)), d-glucose 6-phosphate (d-Glc6-P (7)), d-erythrose 4-phosphate (d-Ery4P (11)), d-ribulose 5-phosphate (d-Rul5P (10)), uracil (Ura (12)), 6-amino-Ura (14), 4-thio-Ura (15), 2-thio-Ura (16), 6-trifluoromethyl-Ura (19), 1-methyl-uracil (17), pyrimidine (29), pyrimidone (30), maleimide (26), adenine (27), guanin (28), thymine (20), and cytosine (23) were from Sigma Aldrich (Vienna, Austria). 3-Methyl-Ura (13), 6-aza-Ura (21), 6-Cl-Ura (22), and l-ribose 5-phosphate (6) were from Carbosynth (Compton, UK). 2-Amino-2-deoxy-d-ribose (2A-d-Rib), 3-amino-3-deoxy-d-ribose (3A-d-Rib), and 2-deoxy-2-fluoro-d-ribose (2F-d-Rib)) were from Synthose (Concord, Ontario, Canada). Expression vectors (pet15b) containing the genes for ribokinase (RbsK) or β-pseudouridine 5′-phosphate glycosidase (YeiN) were obtained from Genescript (Leiden, The Netherlands). Restriction enzymes, Phusion polymerase, T7 RNA Polymerase, Dream Taq DNA polymerase, and GenJET® PCR clean-up kit were from Thermo Fisher Scientific (Altham, MA, USA). DNA sequencing was done at LGC Genomics (Berlin, Germany). Ψ-5’-phosphate glycosidase (YeiN) An N-terminally His6-tagged E. coli YeiN was used. The expression vector (pet15_yein) as obtained from Genescript was transformed into the E. coli BL21_DE3 (pLys) expression strain by electroporation and transformants were selected on 0.1 mg/mL ampicillin LB-agar plates. YeiN was expressed in 1-L baffled shaken flasks at 37 °C and 110 rpm using 250 mL Lennox-media containing 0.1 mg/mL of ampicillin. Inoculation was to an initial OD600 of 0.1. At an OD600 of 0.8, the temperature was decreased to 18 °C and gene expression was induced with 0.4 mM IPTG (isopropyl β-d-1-thiogalactopyranoside) for 20 h. Cells were harvested by centrifugation at 2000 × g at 4 °C for 30 min using a Sorvall RC-5B refrigerated superspeed centrifuge (Du Pont Instruments, Newtown, CT, USA). The supernatant was discarded, and the pellet was suspended in 50 4-(2-hydroxyethyl) piperazine-1-ethanesulfonic acid (HEPES) buffer (pH 8.0). Cells were disrupted using sonication (Fisherbrand Sonic Dismembrator, model Ultrasonic Processor FB-505; Fisher Scientific, Vienna, Austria; 5 min) on ice and the cell-free supernatant was recovered by centrifugation at 4 °C and 2000 × g for 30 min. Pre-treated cell lysate (20 mL) was loaded on 2 × 5 mL HisTrap FF column (GE Healthcare, Little Chalfont, UK), equilibrated with buffer A (50 mM HEPES, pH 8.0 containing 125 mM NaCl, 20 mM imidazole) and mounted on an ÄKTA prime plus (GE Healthcare) system. Protein was eluted using an imidazole gradient from 0% to 80% buffer B (50 mM HEPES, pH 8.0, containing 125 mM NaCl and 500 mM imidazole). The temperature was 10 °C and a flow rate of 3 mL/min was used. Fractions containing the target protein were pooled, concentrated and buffer-exchanged with Amicon Ultra-15 Centrifugal Filter Units (Millipore; Billerica, MA, USA). The final protein concentration was 30 mg/mL in 50 mM HEPES buffer containing 2 mM MgCl2 (pH 7.0). Enzyme was stored at −20 °C until further use. Protein purification was monitored by sodium dodecyl sulphate–polyacrylamide gel electrophoresis (SDS PAGE) (Supplementary Fig. 34). The YeiN activity assay was performed at 37 °C in 50 mM HEPES buffer supplemented with 2 mM MnCl2 (pH 7.0). About 1.5 µM YeiN was used in the reaction. The substrates used were d-Rib5P (1) (2 mM) and Ura (12) (1 mM). One unit of YeiN activity is the amount of enzyme producing 1 μmol of ΨMP/min under the conditions used. Enzyme preparations had typically a specific activity of ∼3 U/mg. Ribokinase (RbsK) The expression vector (pet15_RbsK) as obtained from Genescript was transformed into the E. coli BL21_DE3 (pLys) expression strain by electroporation and transformants were selected on 0.1 mg/mL ampicillin LB-agar plates. RbsK was expressed and purified as described for YeiN, except the protein was stored in 50 mM HEPES supplemented with 2 mM MgCl2 (pH 7.0). An N-terminally His6-tagged E. coli RbsK was used. The RbsK activity assay was performed at 37 °C in 50 mM HEPES buffer (pH 7.0, pH adjusted with KOH) containing 2 mM MgCl2 and 3 µM RbsK. The substrates used were ATP (10 mM) and d-Rib (10 mM). One unit of RbsK activity is the amount of enzyme producing 1 μmol of d-Rib5P (1)/min under the conditions used. The enzyme preparation had a specific activity of ∼50 U/mg. Pentose-5-phosphate phosphatase BT4131 The sugar phosphate phosphatase BT4131 from Bacteroides thetaiotaomicron VPI-5482 was used as an N-terminal fusion with the positively charged binding module Zbasic288. Expression, harvest and cell disruption was performed as described for YeiN except 50 µg/mL kanamycin was used instead of ampicillin. Cell-free extract (20 mL) was loaded on 2 × 5 mL HiTrap SP FF column (GE Healthcare), equilibrated with buffer A (50 mM MES buffer, pH 7.0, containing 100 mM NaCl) and mounted on an ÄKTA prime plus (GE Healthcare) system. Protein was eluted using an NaCl gradient from 0% to 100% buffer B (50 mM HEPES buffer supplemented with 2 M NaCl, pH 7.0) at 10 °C using a flow rate of 3 mL/min. Fractions containing the target protein were pooled, concentrated and buffer-exchanged with Amicon Ultra-15 Centrifugal Filter Units (Millipore; Billerica, MA, USA). The final protein concentration was 10 mg/mL in 50 mM HEPES containing 2 mM MgCl2 (pH 7.0). The enzyme was stored at 4 °C until further use. Protein purification was monitored using SDS PAGE. The specific activity was determined in 50 mM HEPES buffer supplemented with 5 mM MgCl2 (pH 7.0) at 37 °C using 20 mM d-Glc6P as the substrate. The reaction was started with 1.5 µM BT4131 and samples taken were quenched by heat treatment. Inorganic phosphate was measured at 850 nm88. One unit of BT4131 activity is the enzyme amount releasing 1 μmol of phosphate from d-Glc6P/min under the conditions applied. The specific activity of the phosphatase preparation was typically ∼16 U/mg. Nucleoside monophosphate kinases (NPK) The uridine monophosphate kinase (umpk) and cytidine monophosphate kinase (cmpk) genes were amplified from genomic E. coli BL21_DE3 (pLys) DNA by PCR using Phusion polymerase and the following set of primers (restriction sites underlined): UMPK_NdeI_fw (CATATGGCTACCAATGCAAAACCCGTC), UMPK_XhoI_rv (CTCGAGTTATTCCGTGATTAAAGTCCCTTCTTTTTC); CMPK_NdeI_fw (CATATGACGGCAATTGCCCCGGTTATTACC), CMPK_XhoI_rv (CTCGAGTTATGCGAGAGCCAATTTCTGGCGC). Doubly digested (NdeI and XhoI) PCR products were purified from agarose gels treated with shrimp alkaline phosphatase (NEB, Ipswich, MA, USA) and ligated (T4 ligase) into the linearized (NdeI and XhoI) pET-STREP3 expression vector. The constructed vector was transformed into electro-competent E. coli BL21_DE3 (pLys) and transformants were selected on LB-agar plates (0.05 mg/mL kanamycin). Isolated vectors were sequenced (LGC Genomics). NPK expression, cell harvesting and cell disruption was performed as described for YeiN, except that 0.05 mg/mL kanamycin instead of ampicillin was used and cells were dissolved in wash buffer (100 mM Tris-Cl pH 8.0, 150 mM NaCl, 1 mM EDTA). Cell-free extract (20 mL) was loaded (2 mL/min) on 1 × 5 mL StrepTrap HP column (GE Healthcare), equilibrated with wash buffer and mounted on an ÄKTA prime plus (GE Healthcare) system. Unspecifically bound protein was washed off with 10 column volumes of wash buffer at a flow of 3 mL/min. Protein was eluted using elution buffer (100 mM Tris-Cl pH 8.0, 150 mM NaCl, 1 mM EDTA, 2.5 mM desthiobiotin). Fractions containing the target protein were pooled, concentrated and buffer-exchanged with Amicon Ultra-15 Centrifugal Filter Units (Millipore). The final protein concentration was 15 mg/mL in 50 mM HEPES containing 2 mM MgCl2 (pH 7.0). Enzyme was stored at −20 °C until further use. Protein purification was monitored using SDS PAGE (Supplementary Fig. 34). The UMPK was only weakly expressed and most of the recombinant enzyme was found in inclusion bodies. An N-terminally strep-tagged CMPK was used. The NPK-assay was performed at 37 °C in 50 mM HEPES buffer (pH 7.0) containing 1 mM MgCl2 and 5 µM CMPK. The substrates used were 3 mM UMP and 3 mM ATP. One unit of NPK activity is the amount of enzyme producing 1 μmol of UDP/min under the conditions used. CMPK preparation typically exhibited a specific activity of 20 U/mg. Nucleoside diphosphate kinase (NDK) An N-terminally Strep-tagged variant was used. The ndk gene was amplified from genomic E. coli BL21_DE3 (pLys) DNA by PCR using Phusion polymerase and the following set of primers (overlapping 5'-overhangs in lower case letters): NDK_fwd (ATGGCTATTGAACGTACTTTTTCCATC), NDK_rv (TTAACGGGTGCGCGGGC). The PCR product was cloned into the pET-STREP3 expression vector by circular polymerase extension cloning. Briefly, using the primer pairs CEPEC_NDK_fwd (GTTCGAGAAAGGCTTAATTAACCATATGgctattgaacgtactttttccatc), CPEC-NDK_rv (GTGGTGGTGGTGGTGCTCGAGttaacgggtgcgcgggc) and CPEC-NDK-Strep_fwd (ggaaaaagtacgttcaatagcCATATGGTTAATTAAGCCTTTCTCGAAC) CPEC-NDK-Strep_rv (gcccgcgcacccgttaaCTCGAGCACCACCACCACCAC) in two independent PCR reactions, complementary overlaps are introduced into the ndk gene and into the target vector. The purified PCR products were used in an overlap-extension-PCR to prime each other and to create the final expression vector, which was transformed into E. coli XL1 Blue. Transformants were selected on LB-agar plates containing 0.1 mg/mL of ampicillin, and the isolated plasmid was sequenced. The sequenced vector was transformed into the expression strain E. coli BL21_DE3 (pLys) and transformants were selected on LB-agar plates containing 0.1 mg/mL of ampicillin. NDK expression and purification were as described for NDKs (Supplementary Fig. 34). The NDK activity assay was performed at 37 °C in 50 mM HEPES buffer (pH 7.0) containing 1 mM MgCl2 and 5 µM NDK. The substrates used were 3 mM uridine diphosphate (UDP) and 3 mM ATP. One unit of NDK activity is the amount of enzyme producing 1 μmol of UTP/min under the conditions used. The enzyme preparation had a specific activity of 21 U/mg. Determination of protein concentration Protein concentrations were measured using DS-11 Spectrophotometer (DeNovix; Wilmington, DE, USA) at 280 nm. Protein concentration was calculated using the corresponding molar extinction coefficient and the molecular weight computed by protparam (https://web.expasy.org/protparam/). RbsK (12,490 M−1cm−1, 32,290 Da), YeiN (10,095 M−1cm−1, 32,909 Da), CMPK (13,075 M−1cm−1, 25,970 Da), NDK (4470 M−1cm−1, 15,463 Da), BT4131 (14,690 M−1cm−1, 28,856 Da). Determination of nucleotide concentration Uracil content was measured photometrically using a DS-11 Spectrophotometer (DeNovix). Ura (12) (260 nm, 10 absorbance units (AU) mM−1 cm−1), 4-thio-Ura (15) (270 nm, 22 AU mM−1 cm−1), 2-thio-Ura (16) (273 nm, 12 AU mM−1 cm−1), thymidine (20) (260 nm, 8.9 AU mM−1 cm−1), guanine (28) (246 nm, 10.7 mM−1 cm−1), cytosine (23) (267 nm, 6.1 mM−1 cm−1), adenine (27) (260 nm, 13 mM−1 cm−1). High-performance liquid chromatography (HPLC) Nucleosides were quantified by reversed-phase ion-pairing HPLC. Typically, 5 µL of sample containing an overall nucleoside concentration of about 1 mM were loaded on a Kinetex C18 column (Phenomenex, Aschaffenburg, Germany; 5 µm, 100 Å, 50 × 4.6 mm). Analytes were separated in 5-min long isocratic runs using 20 mM phosphate buffer, pH 5.9, containing 40 mM tetra-n-butylammonium bromide (TBAB) and 12.5% acetonitrile. The flow rate was 2 mL/min and the temperature set to 35 °C. Eluents were detected at 260 nm. A 5-min HPLC trace is shown in Supplementary Fig. 35. Typical retention times were as follows: Ura (12) (0.3 min), ΨMP (0.4 min), ΨDP (0.9 min), ADP (1.1 min), ΨTP (2.7 min), and ATP (4 min). Thin-layer chromatography (TLC) About 1–2 µL sample were spotted onto a TLC plate (Merk, Darmstadt, Germany). The mobile phase was composed of 2-BuOH:AcOH:H2O (2:1:1). The development solution consisted of 0.5 g thymol, 95 mL EtOH and 5 mL H2SO4. Reducing sugar spots were developed by heating the plate to approximately 70 °C for 2 min using a heat-gun (Supplementary Fig. 33). Synthes is and isolation of pentose 5-phosphates Prior to synthesis, the substrate scope of RibK was evaluated using mixture of 100 mM sugar (d-Rib, d-dRib, d-Xyl, d-Ara, 2F-d-Rib, d-Glc, 2F-d-Rib, 3A-d-Rib, l-Xyl, l-Ara, d-Lyx, d-Man, d-GlcN, d-GlcNAc, d-GlcA, d-Gal, d-GalN), 100 mM PEP, 2 mM ATP and 10 mM HEPES buffer (pH 7.0) containing 1 mM MgCl2 and enzymes (60 µM RbsK, 0.4 µM PK) in a total volume of 100 µL. The reactions were incubated at 37 °C without agitation for 20 h and 5 µL sample were analyzed by TLC. Of the tested sugars, d-Rib, d-dRib, d-Xyl, d-Ara,2F-d-Rib, d-Glc were accepted. Synthesis was performed at 5 mL scale, except for 2-deoxy-2-fluoro-d-Rib 5-phosphate (2F-D-Rib5P (5) that was prepared in a 1 mL reaction. Reaction mixtures contained 100 mM pentose (d-Rib, d-dRib, d-Xyl, d-Ara, 2F-d-Rib), 100 mM PEP, 2 mM ATP in 10 mM HEPES buffer (pH 7.0) with 1 mM MgCl2 and 0.4 µM PK. For the phosphorylation of d-Rib, 15 µM RbsK, for d-dRib 30 µM RbsK, and for other pentoses 60 µM RbsK were used. The reactions were incubated at 37 °C without agitation for 20 h and 5 µL sample were analyzed by TLC. The reaction was stopped by heat treatment at 99 °C for 10 min and precipitated protein was removed by centrifugation for 10 min at 4 °C and 16,000 × g. The reactions were diluted with twice volume of water and applied (20 mL; 2 mL/min) to a FliQ FPLC column (10 mL; 104 × 11.0 mm; Generon, Maidenhead, U.K.) packed with SuperQ-650M (Tosh, Tokyo, Japan) anion exchange resin mounted onto an ÄKTA prime plus FPLC system (GE Healthcare). The column was equilibrated in 10 mM sodium acetate (NaAc), pH 4.5. Products were eluted using a linear gradient of 1 M NaAc (pH 4.5), from 0% to 50% in 120 mL with a flow of 4 mL/min. Eluting fractions were evaluated for the presence of sugar phosphates using TLC. Following the product elution, the column was regenerated with 100% of 1 M NaAc buffer for 10 min. The fractions containing pentose 5-phosphates were pooled and concentrated under reduced pressure (20 mbar, 40 °C) to a total volume of 1 mL using a rotary evaporator (Laborta 500-efficient; Heidolph Instruments, Schwabach, Germany). The concentration of pentose 5-phosphates were evaluated as follows: sugar phosphates were completely hydrolyzed by applying 10 µM of sugar phosphate phosphatase (AGP)88 for 20 h. The released inorganic phosphate was determined in a colorimetric assay based on the formation of a blue phospho-molybdate complex at 850 nm89. Screening sugar phosphate and nucleobase substrates All reactions were performed at 37 °C in 50 mM HEPES buffer (pH 7.0) containing 2 mM MnCl2 in a total volume of 100 μL. Nucleobase screening was performed by incubating 5 mM d-Rib5P (1) and 1 mM of nucleobase (Ura (12), 3-methyl-Ura (13), 6-amino-Ura (14), 2-thio-Ura (15), 4-thio-Ura (16), 1-methyl-Ura (17), 5-fluoro-Ura (18), 6-trifluoro-methyl-Ura (19), 6-aza-Ura (21), 6-Cl-Ura (22), 2,4-dithio-Ura (24), 5-ethynyl-Ura (25), maleimide (26), cytosine(23), thymine (20), adenine (27), guanine (28), 4-pyrimidone (29) and pyrimidine (30)) with 30 µM YeiN. Sugar phosphates were screened by incubating 1 mM Ura and 5 mM of sugar phosphate (d-Rib5P (1), d-dRib5P (2), d-Xyl5P (3), d-Ara5P (4), d-Fru6P (9), d-Man6P (8), d-Glc6P (7), d-Rul (10), d-Ery4P (11), 2F-d-Rib5P (5), l-Rib5P (6) with 30 µM YeiN. Samples (20 µL) were withdrawn after 0, 3, 24 h and quenched by methanol addition (100 µL). Precipitated protein was removed by centrifugation at 20,000 × g for 5 min and supernatant was analyzed by HPLC. Five millimolar of pentose phosphates (previously isolated as described above) accepted by YeiN in the initial screen (d-Rib5P (1), d-dRib5P (2), d-Xyl5P (3), d-Ara5P (4), d-2F-Rib5P(5)) were incubated with 2 mM of each nucleobase (Ura (12), 3-methyl-Ura (13), 6-amino-Ura (14), 2-thio-Ura (15), 4-thio-Ura (16), 1-methyl-Ura (17), 5-fluoro-Ura (18), 6-trifluoro-methyl-Ura (19), 6-aza-Ura (21), 6-Cl-Ura (22), 2,4-dithio-Ura (24), 5-ethynyl-Ura (25), maleimide (26), cytosine(23), thymine (20), adenine (27), guanin (28), 4-pyrimidone (29) and pyrimidine (30)) and 30 µM of YeiN for 3 h. The reactions were quenched by methanol addition (500 µL) and precipitated protein was removed by centrifugation at 20,000 × g for 5 min prior to HPLC analysis. YeiN kinetics The rate constant kcat and the Michaelis constant KM were determined by measuring the initial rate of ΨMP synthesis catalyzed by YeiN (1 µM–120 µM) in 50 mM HEPES (pH 7.0) supplemented with 2 mM MnCl2. Reactions were incubated at 37 °C on a Thermomixer comfort (Eppendorf, Hamburg, Germany) with agitation of 300 rpm. Initial rates for ΨMP synthesis with nucleobases (Ura (12), 6-amino-Ura (13), 4-thio-Ura (14), 2-thio-Ura (15), 3-methy-Ura(16)) were determined at a constant d-Rib5P (1) concentration of 5 mM and at nucleobase concentrations between 0.05 mM and 1.00 mM. Initial rates for ΨMP synthesis with pentose 5-phosphates (d-Rib5P (1), d-dRib5P (2), d-Ara5P (4), d-Xyl5P(3) were determined at a constant Ura (10) concentration of 5 mM and 1.00−0.05 mM of the corresponding pentose 5-phosphate. Reactions were started by YeiN (1.00 µM−100 µM) added in a concentration adjusted to the enzyme’s specific activity with the substrate used. The reaction mixture was incubated for 20 min and quenched by methanol addition (1:1, v:v). Precipitated protein was removed by centrifugation at 20,000 × g for 5 min and samples were analyzed by HPLC. It was ensured that conversion was below 20%. Kinetic parameters were obtained from non-linear least squares fits (SigmaPlot; Systat, Erkrath, Germany) of initial ΨMP-synthesis rates to Michaelis–Menten Eq. (1). V is the initial rate, Vmax is the maximum initial rate, KM is the Michaelis constant and [S] is the initial substrate concentration. The kcat was determined from the relationship kcat = Vmax/[E], where [E] is the molar enzyme concentration. [E] was determined from the protein concentration and enzyme molecular weight. 1 V=Vmax[S](KM+[S]) Synthesis of ΨMP derivatives Synthesis of ΨMP derivatives was performed in two steps in one pot. Reactions were started by addition of enzyme. All reactions were incubated at 37 °C without agitation. Step one: pentose phosphorylation. The reaction mixture contained 100 mM pentose (d-Rib, d-dRib, d-Xyl, d-Ara, 2F-d-Rib), 100 mM PEP, 2 mM ATP and 10 mM HEPES buffer (pH 7.0), 1 mM MgCl2 and 0.4 µM pyruvate kinase. For the phosphorylation of d-Rib 15 µM RbsK, for d-dRib 30 µM RbsK and for d-Ara and d-Xyl 60 µM RbsK was used. Reactions were performed in a total volume of 5 mL. Samples (5 µL) were taken at reaction start and then after 2 h and 20 h. They were analyzed by TLC. Step two: C-glycosylation using YeiN. After phosphorylation, the reaction mixture was diluted to a final concentration of 25 mM pentose 5-phosphate (d-Rib5P (1), d-dRib5P (2), d-Xyl5P (3), d-Ara5P(4)) using 10 mM HEPES (pH 7.0). Ten mM of nucleobase (Ura (12), 6-amino-Ura (13), 2-thio-Ura (14), 4-thio-Ura (15), 3-methyl-Ura (16)) were added to the phosphorylation reaction. ΨMP-synthesis was started with 15 µM YeiN. To compensate for the lower activity of YeiN towards different pentose 5-phosphates and nucleobases, the YeiN concentration was adjusted as follows: for the synthesis of dΨMP 30 µM YeiN was used; for synthesis of Ara-ΨMP, Xyl-ΨMP, 6-amino-ΨMP, 2-thio-ΨMP and 4-thio-ΨMP, 90 µM YeiN was used. Samples (5 µL) were taken at start and after 5, 15, 30, 60, 120, 180, 240, 300, 360 min, and 10 h for 2-thio-ΨMP, Ara-ΨMP, Xyl-ΨMP. Each sample was quenched by 100 µL methanol. Precipitated protein was removed by centrifugation for 20,000 × g for 5 min and the sample was analyzed by HPLC. After completion of the reaction or latest after 10 h of incubation, the reactions were stopped by heat inactivation for 5 min at 99 °C. Precipitated protein was removed by centrifugation for 10 min at 4 °C and 16,000 ×  g. ΨMP dephosphorylation ΨMP (5 mM) was dissolved in 15 mL 10 mM HEPES buffer, pH 8.0, supplemented with 1 mM MgCl2 and incubated with 2 U/mL of CIP at 37 °C for 600 min. Ten microliters samples were withdrawn at start, and after 30, 90, 600 min, and quenched in 100 µL methanol. Precipitated enzyme was removed by centrifugation (16,000 × g, 4 °C, 5 min) and supernatant was analyzed by HPLC. CMPK, NDK, PK, phosphorylation cascade The reaction set-up consisted of 0.5 mM ATP, 5 mM purified ΨMP, 12 µM CMPK, 1 mM MgCl2, 30 µM PK, 12 mM PEP, 50 mM HEPES, pH 7.0, in a total volume of 250 µL. Ten µL samples taken at certain times were quenched in 100 µL methanol. The precipitated enzyme was removed by centrifugation (20,000 × g, 4 °C, 5 min) and supernatant was analyzed by HPLC. To test the effect of NDK removal for ΨTP-synthesis the same reaction set up without addition of 40 µM NDK was used (Supplementary Fig. 31). One-pot multistep enzymatic synthesis of Ψ-triphosphates Pentose phosphorylation was performed as described above under “Step one: pentose phosphorylation”. For ΨMP synthesis, 15 mM nucleobase (Ura (12), 6-amino-Ura (13), 2-thio-Ura (14), 4-thio-Ura (15), 3-methyl-Ura (16)), 1 mM MnCl2, and 10 mM HEPES buffer were added to the reaction mixture to a final concentration of 25 mM pentose 5-phosphate. The reaction was started with YeiN (1 µM for the ΨMP synthesis, 3 µM for the dΨMP synthesis and 90 µM for the synthesis using Ura derivatives). ΨMP phosphorylation was started with 15 mM PEP, 15 µM CMPK, and 1 mM ATP, resulting in a volume increase of 33%. Ten microliters of samples were quenched by addition of 100 μL methanol and precipitated protein was removed by centrifugation for 20,000 × g for 5 min prior to HPLC analysis. Enzymes were heat inactivated (5 min, 99 °C) after 21 h reaction and precipitated protein was removed by centrifugation for 10 min at 4 °C and 16,000 × g. Dephosphorylation of pentose 5-phosphates BT4131 (10 µM) and 1 mM MgCl2 were added to the reaction mixtures and incubated for 60 min of at 37 °C. The reaction was stopped by heat inactivation (5 min, 99 °C) and precipitated protein was removed by centrifugation for 10 min at 4 °C and 16,000 ×  g. Anion exchange chromatography Enzyme-free reaction mixture was applied (20 mL; 2 mL/min) to a FliQ FPLC column (10 mL; 104 × 11.0 mm; Generon, Maidenhead, U.K.) packed with SuperQ-650M (Tosh, Tokyo, Japan) anion exchange resin mounted onto an ÄKTA prime plus FPLC system (GE Healthcare.). The column was equilibrated in 10 mM NaAc, pH 4.5. Products were eluted using a linear gradient of 1 M NaAc, pH 4.5, from 0% to 50% in 120 mL with a flow of 4 mL/min and monitored at 260 nm. Following product elution, buffer B was set to 100% for 10 min to elute unspecifically bound substances. The fractions containing the desired C-nucleoside were pooled and concentrated under reduced pressure (20 mbar, 40 °C) to a total volume of 10 mL using a rotary evaporator (Laborta 500-efficient; Heidolph Instruments, Schwabach, Germany). Size exclusion chromatography To remove excess NaAc, concentrated products were split in 5 × 2 mL portions and loaded onto a size exclusion column (16 × 1000 mm; XR 16/100 column; GE Healthcare) packed with Sephadex G10 (exclusion limit <700 Da) equilibrated in double distilled H2O and mounted onto an ÄKTA prime plus FPLC system (GE Healthcare). Compound elusion was performed using doubly distilled water and monitored at 260 nm. Fractions containing salt-free product were collected, pooled and concentrated under reduced pressure (20 mbar, 40 °C) to a total volume of 10 mL using a rotary evaporator (Laborta 500-efficient; Heidolph Instruments) and lyophilized. Docking Molecular docking was used to structurally evaluate the substrate scope of YeinN. AutoDock1390 as implemented in Yasara v. 18.2.7 was used for enzyme-ligand docking. The AmberFB-1591 force field and the default parameters provided by the standard docking macro were used. The structure of YeiN (pdb code: 4GIL)62 was used as a receptor in molecular-docking experiments. For dockings of sugar phosphates (d-Rib5P (1), d-Rul5P (10), d-Fru6P (9), d-Glu6P (7), d-Ery4P(11)), the ring-opened intermediate was deleted and a simulation cell of 15 Å × 10 Å × 10 Å was placed at the former location of d-Rib5P (1). Sugar phosphates were docked in the ring-open conformation. The protonation state of key catalytic residues (Glu31 protonated; Lys166 deprotonated) was set prior to the docking experiment. The protonation states of all other protein residues were set automatically for a pH of 7.00. For dockings of nucleobases (Ura (12), 6-trifluoro-methyl-Ura (19), 6-Cl-Ura (22), 6-aza-Ura (21), 1-methyl-Ura (17), 5-fluoro-Ura (18), 5-methyl-Ura(20)), the uracil moiety was deleted and the charge of lysine was set to fit the covalent iminium intermediate. The protonation states of all other protein residues were set automatically for a pH of 7.00. A simulation cell of 15 Å × 10 Å × 10 Å was placed at the former location of Ura (11). Following the docking runs, the obtained poses were optimized by energy using the standard macro included in Yasara. All ligands were generated using the Grade Web Server (http://grade.globalphasing.org/cgi-bin/grade/server.cgi)92. Docking poses were evaluated by their associated free energy and mechanistic plausibility. Synthesis of pseudouridinated RNA In vitro translation (IVT) using T7 RNA polymerase (Thermo Fisher Scientific) was performed following the manufacturer’s protocol. Shortly, linear DNA template containing the T7 promoter and the yein gene was amplified from the YeiN expression vector (peT15_yein) by PCR using Phusion polymerase and a set of primers complementary to the T7 promoter (TAATACGACTCACTATAGGG) or terminator (GCTAGTTATTGCTCAGCGG). The PCR product was purified via column purification (GenJET® PCR clean up; Thermo Fisher Scientific) and eluted in DEPC-water. IVT was performed using 30 U of T7 RNA polymerase, 600 ng template DNA, Transcription Buffer (40 mM Tris-HCl pH 7.9, at 25 °C), 6 mM MgCl2, 10 mM DTT, 10 mM NaCl and 2 mM spermidine) and 2 mM of each NTP (ATP, GTP, UTP, CTP) in a total volume of 50 µL. Synthesis of modified RNA was performed in two ways; either 2 mM of ΨTP (or: dΨTP, 6-amino-ΨTP) was additionally added to the conventional NTPs or, UTP was fully replaced by 2 mM of either ΨTP, dΨTP or 6-amino-ΨTP. The reaction was incubated at 37 °C in the thermocycler (Doppio-Dual; VWR, Radnor, PA, USA) and the resulting PCR products were evaluated by agarose gel electrophoresis. Synthesis of pseudouridinated DNA To evaluate the potential of PCR catalyzed generation of pseudouridinated DNA, the yein gene was amplified from 50 ng expression vector (pet15_yein) using 1.25 U of Dream-Taq DNA polymerase, 0.2 mM of dNTPs (dATP, dCTP, dGTP, dTTP), 0.5 µM T7_fwd, 0.5 µM T7_bwd primers and DreamTaq Buffer as recommended by the supplier (Thermo Fisher Scientific). To incorporate Ψ into the DNA either 0.2 mM of ΨTP were additionally added to dNTP mix or fully replaced dTTP by dΨTP. The temperature profile consisted of a preheating step at 95 °C for 1 min, followed by 30 reaction cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and elongation at 72 °C for 2 min. The final extension step was carried out at 72 °C for 5 min. PCRs were performed using the Doppio-Dual thermocycler (VWR). The resulting PCR products were evaluated by agarose gel electrophoresis. In situ 1H-NMR of ΨMP synthesis or hydrolysis YeiN-catalyzed synthesis of ΨMP from d-Rib5P (1) and Ura (12), as well as hydrolysis of ΨMP, were analyzed by 1H-NMR. All reagents including the enzyme were dissolved in D2O and the reaction was set up in a 5 mm high-precision NMR sample tube (Promochem, Wesel, Germany). YeiN was buffered by washing the enzyme with 10-fold excess of 10 mM potassium phosphate buffer (pD 7.5) supplemented with 0.1 mM MnCl2, using centrifugation (16,000 × g, 15 min, 4 °C) with Vivaspin 500 centrifugal concentrators (10 kDa cut-off). The synthesis reaction contained 14 mM Ura (12) and 8 mM d-Rib5P (1), while 10 mM ΨMP were used for the hydrolysis reaction. Each reaction was initiated by adding 3 µM enzyme. The reactions were incubated at 30 °C without agitation in the electromagnetic field. Measurements were performed with a Varian INOVA 500-MHz NMR spectrometer (Agilent Technologies; Santa Clara, CA, USA). The VNMRJ 2.2D software was used for all measurements. 1H-NMR spectra (499.98 MHz) were recorded on a 5 mm indirect detection PFG-probe with pre-saturation of the water signal by a shaped pulse. The following standard pre-saturation sequence was used: 2 s relaxation delay; 90° proton pulse; 2.048 s acquisition time; 8 kHz spectral width; number of points 32 k. 1H-NMR analysis of YeiN-catalyzed deuterium exchange The potential of YeiN-catalyzed deuterium wash-in into Ura (12) without transient formation of ΨMP was evaluated by 1H-NMR. All reagents including the enzyme were dissolved in D2O and prepared as described above. The 600 µL reactions contained 5 mM Ura (12) and 5 mM of sugar phosphate (d-Rib5P (1), d-Ery4P (11), d-Rul5P (10), d-Glc6P (7), d-Fru6P(9)) in 10 mM potassium phosphate buffer (pD 7.5) supplemented with 0.1 mM MnCl2. Conversions were started by addition of 3 µM enzyme and incubated at 30 °C in a Thermomixer comfort (Eppendorf, Hamburg, Germany) with agitation of 300 rpm for 3 h. The reactions were stopped by heat treatment (5 min, 99 °C) and precipitated protein was removed by centrifugation for 10 min at 4 °C and 16,000 × g. The supernatant was transferred into a 5 mm NMR sample tube and analyzed by 1H-NMR. NMR analysis of purified reaction products Lyophilized products were dissolved to 3−10 mM in 50 mM potassium phosphate buffer (pD 7.9) in a total volume of 600 μL and transferred to a 5 mm high-precision NMR sample tube. They were analyzed using a Varian INOVA 500-MHz NMR spectrometer (Agilent Technologies). The VNMRJ 2.2D software was used for measurements. 1H-NMR spectra (499.98 MHz) were measured on a 5 mm indirect detection PFG probe. Standard pre-saturation sequence was used: relaxation delay 2 s; 908 proton pulse; 2.048 s acquisition time; spectral width 8 kHz; number of points 32,000 13C NMR spectra (125.71 MHz). Standard pre-saturation sequence was used: relaxation delay 2 s; 908 proton pulse; 2.048 s acquisition time; spectral width 8 kHz; number of points 32 k. MNova 12.0.3 (Mestrelab Research S.L. Bajo, Santiago de Compostela, Spain) was used for data processing and analysis (Supplementary Figs. 10–29). Recorded spectra were in good agreement with published data75. Reporting summary Further information on research design is available in the Nature Research Reporting Summary linked to this article. Supplementary information Supplementary Information Peer Review File Reporting Summary Source data Source Data Peer review information Nature Communications thanks Pere Clapes and the other, anonymous reivewer(s) for their contribution to the peer review of this work. Peer review reports are available. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary information Supplementary information is available for this paper at 10.1038/s41467-020-20035-0. Acknowledgements We thank Alma Memic and Stefan Drießler (Institute of Biotechnology and Biochemical Engineering, Graz University of Technology) for experiments; Prof. Hansjörg Weber (Institute of Organic Chemistry, Graz University of Technology) for NMR measurements. Financial support from the Austrian Science Funds (FWF, DK Molecular Enzymology W901; to B.N. and M.P.) is gratefully acknowledged. Author contributions M.P. design of study; all experiments and data analysis; writing of paper. B.N. design of study; funding acquisition; discussion; writing of paper. Data availability All relevant data are reported in the manuscript and in the associated Supplementary Information. All data are available from the corresponding author upon reasonable request. There is no restriction on data availability. Source data are provided with this paper. Competing interests The authors declare no competing interests. ==== Refs References 1. Evans GB 8-Aza-immucillins as transition-state analogue inhibitors of purine nucleoside phosphorylase and nucleoside hydrolases J. Med. Chem. 2003 46 155 160 10.1021/jm0203332 12502369 2. Kicska GA Immucillin H, a powerful transition-state analog inhibitor of purine nucleoside phosphorylase, selectively inhibits human T lymphocytes Proc. Natl Acad. Sci. 2001 98 4593 4598 10.1073/pnas.071050798 11287638 3. Miles RW Tyler PC Furneaux RH Bagdassarian CK Schramm VL One-third-the-sites transition-state inhibitors for purine nucleoside phosphorylase Biochemistry 1998 37 8615 8621 10.1021/bi980658d 9628722 4. Berg M Evaluation of nucleoside hydrolase inhibitors for treatment of african trypanosomiasis Antimicrob. Agents Chemother. 2010 54 1900 1908 10.1128/AAC.01787-09 20194690 5. Jan Š Hocek M Pavel K C-Nucleosides: synthetic strategies and biological applications Chem. Rev. 2009 109 6729 6764 10.1021/cr9002165 19761208 6. Eyer, L., Nencka, R., de Clercq, E., Seley-Radtke, K. & Růžek, D. Nucleoside analogs as a rich source of antiviral agents active against arthropod-borne flaviviruses. Antivir. Chem. Chemother. 26, 2040206618761299 (2018). 7. Warren TK Protection against filovirus diseases by a novel broad-spectrum nucleoside analogue BCX4430 Nature 2014 508 402 405 10.1038/nature13027 24590073 8. De Clercq E C-Nucleosides to be revisited J. Med. Chem. 2016 59 2301 2311 10.1021/acs.jmedchem.5b01157 26513594 9. Rosenqvist P Characterization of C-nucleoside antimicrobials from Streptomyces albus DSM 40763: strepturidin is pseudouridimycin Sci. Rep. 2019 9 1 9 10.1038/s41598-019-45375-w 30626917 10. Li G Clercq ED Therapeutic options for the 2019 novel coronavirus (2019-nCoV) Nat. Rev. Drug Discov. 2020 19 149 150 10.1038/d41573-020-00016-0 32127666 11. Zumla A Chan JFW Azhar EI Hui DSC Yuen K-Y Coronaviruses—drug discovery and therapeutic options Nat. Rev. Drug Discov. 2016 15 327 347 10.1038/nrd.2015.37 26868298 12. Warren TK Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys Nature 2016 531 381 385 10.1038/nature17180 26934220 13. Wang M Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro Cell Res. 2020 30 269 271 10.1038/s41422-020-0282-0 32020029 14. Li Z Rapid review for the anti-coronavirus effect of remdesivir Drug Discov. Ther. 2020 14 73 76 10.5582/ddt.2020.01015 32378648 15. Kong L Divergent biosynthesis of C-nucleoside minimycin and indigoidine in bacteria iScience 2019 22 430 440 10.1016/j.isci.2019.11.037 31816530 16. Roy-Burman S Roy-Burman P Visser DW Showdomycin, a new nucleoside antibiotic Cancer Res. 1968 28 1605 1610 5691758 17. Maffioli SI Antibacterial nucleoside-analog inhibitor of bacterial RNA polymerase Cell 2017 169 1240 1248 10.1016/j.cell.2017.05.042 28622509 18. Sosio M Analysis of the pseudouridimycin biosynthetic pathway provides insights into the formation of C-nucleoside antibiotics Cell Chem. Biol. 2018 25 540 549 10.1016/j.chembiol.2018.02.008 29551347 19. Ren D Identification of the C-glycoside synthases during biosynthesis of the pyrazole-C-nucleosides formycin and pyrazofurin Angew. Chem. Int. Ed. 2019 58 16512 16516 10.1002/anie.201910356 20. Palmu K Discovery of the showdomycin gene cluster from Streptomyces showdoensis ATCC 15227 yields Insight into the biosynthetic logic of C-nucleoside antibiotics ACS Chem. Biol. 2017 12 1472 1477 10.1021/acschembio.7b00078 28418235 21. Gao S Uncovering the chemistry of C–C bond formation in C-nucleoside biosynthesis: crystal structure of a C-glycoside synthase/PRPP complex Chem. Commun. 2020 56 7617 7620 10.1039/D0CC02834G 22. Carlile TM Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells Nature 2014 515 143 146 10.1038/nature13802 25192136 23. Harcourt EM Kietrys AM Kool ET Chemical and structural effects of base modifications in messenger RNA Nature 2017 541 339 346 10.1038/nature21351 28102265 24. Charette M Gray MW Pseudouridine in RNA: what, where, how, and why IUBMB Life 2000 49 341 351 10.1080/152165400410182 10902565 25. Li X Ma S Yi C Pseudouridine: the fifth RNA nucleotide with renewed interests Curr. Opin. Chem. Biol. 2016 33 108 116 10.1016/j.cbpa.2016.06.014 27348156 26. Spenkuch F Motorin Y Helm M Pseudouridine: still mysterious, but never a fake (uridine)! RNA Biol. 2014 11 1540 1554 10.4161/15476286.2014.992278 25616362 27. Hamma T Ferré-D’Amaré AR Pseudouridine synthases Chem. Biol. 2006 13 1125 1135 10.1016/j.chembiol.2006.09.009 17113994 28. Veerareddygari GR Singh SK Mueller EG The pseudouridine synthases proceed through a glycal intermediate J. Am. Chem. Soc. 2016 138 7852 7855 10.1021/jacs.6b04491 27292228 29. Mandal PK Rossi DJ Reprogramming human fibroblasts to pluripotency using modified mRNA Nat. Protoc. 2013 8 568 582 10.1038/nprot.2013.019 23429718 30. Uchida S Kataoka K Itaka K Screening of mRNA chemical modification to maximize protein expression with Rreduced immunogenicity Pharmaceutics 2015 7 137 151 10.3390/pharmaceutics7030137 26213960 31. Karikó K Incorporation of pseudouridine into mRNA yields superior nonimmunogenic vector with increased translational capacity and biological stability Mol. Ther. 2008 16 1833 1840 10.1038/mt.2008.200 18797453 32. Karikó K Muramatsu H Keller JM Weissman D Increased erythropoiesis in mice injected with submicrogram quantities of pseudouridine-containing mRNA encoding erythropoietin Mol. Ther. 2012 20 948 953 10.1038/mt.2012.7 22334017 33. Eyler DE Pseudouridinylation of mRNA coding sequences alters translation Proc. Natl Acad. Sci. 2019 116 23068 23074 10.1073/pnas.1821754116 31672910 34. Karijolich J Yu Y-T Converting nonsense codons into sense codons by targeted pseudouridylation Nature 2011 474 395 398 10.1038/nature10165 21677757 35. Hoernes TP Nucleotide modifications within bacterial messenger RNAs regulate their translation and are able to rewire the genetic code Nucleic Acids Res. 2016 44 852 862 10.1093/nar/gkv1182 26578598 36. Fernández IS Unusual base pairing during the decoding of a stop codon by the ribosome Nature 2013 500 107 110 10.1038/nature12302 23812587 37. Li B Engineering CRISPR–Cpf1 crRNAs and mRNAs to maximize genome editing efficiency Nat. Biomed. Eng. 2017 1 1 10 10.1038/s41551-016-0001 38. Vaidyanathan S Uridine depletion and chemical modification increase Cas9 mRNA activity and reduce Iimmunogenicity without HPLC purification Mol. Ther. Nucleic Acids 2018 12 530 542 10.1016/j.omtn.2018.06.010 30195789 39. Temburnikar K Seley-Radtke KL Recent advances in synthetic approaches for medicinal chemistry of C-nucleosides Beilstein J. Org. Chem. 2018 14 772 785 10.3762/bjoc.14.65 29719574 40. Daves, G. D. & Cheng, C. C. in Progress in Medicinal Chemistry (eds. Ellis, G. P. & West, G. B.) Vol. 13, 303–349 (Elsevier, 1976). 41. Hacksell, U. & Daves, G. D. in Progress in Medicinal Chemistry (eds. Ellis, G. P. & West, G. B.) Vol. 22, 1–65 (Elsevier, 1985). 42. Meanwell M A short de novo synthesis of nucleoside analogs Science 2020 369 725 730 10.1126/science.abb3231 32764073 43. Xue F Improvement of the C-glycosylation step for the synthesis of remdesivir Org. Process Res. Dev. 2020 24 1772 1777 10.1021/acs.oprd.0c00310 44. Shapiro R Chambers R W. synthesis of pseudouridine J. Am. Chem. Soc. 1961 83 3920 3921 10.1021/ja01479a057 45. van Rijssel ER Stereoselectivity in the lewis acid mediated reduction of ketofuranoses J. Org. Chem. 2015 80 4553 4565 10.1021/acs.joc.5b00419 25826382 46. Hanessian S Machaalani R A highly stereocontrolled and efficient synthesis of α- and β-pseudouridines Tetrahedron Lett. 2003 44 8321 8323 10.1016/j.tetlet.2003.09.103 47. Yu C-P Chang H-Y Chien T-C Total synthesis of pseudouridine via Heck-type C-glycosylation N. J. Chem. 2019 43 8796 8803 10.1039/C9NJ01012B 48. Grohar PJ Chow CS A practical synthesis of the modified RNA nucleoside pseudouridine Tetrahedron Lett. 1999 40 2049 2052 10.1016/S0040-4039(99)00162-8 49. Lerch U Burdon MG Moffatt JG C-Glycosyl nucleosides. synthesis of pseudouridine and related compounds J. Org. Chem. 1971 36 1507 1513 10.1021/jo00810a016 5574816 50. Brown, D. M. & Ogden, R. C. A synthesis of pseudouridine. J. Chem. Soc. Perkin 1 723–725 (1981). 10.1039/P19810000723. 51. Hoang C Ferré-D’Amaré AR Cocrystal structure of a tRNA Ψ55 pseudouridine synthase: nucleotide flipping by an RNA-modifying enzyme Cell 2001 107 929 939 10.1016/S0092-8674(01)00618-3 11779468 52. Foster PG Huang L Santi DV Stroud RM The structural basis for tRNA recognition and pseudouridine formation by pseudouridine synthase I Nat. Struct. Biol. 2000 7 23 27 10.1038/71219 10625422 53. Blauenburg B Oja T Klika KD Metsä-Ketelä M Chemoenzymatic synthesis of novel C-ribosylated naphthoquinones ACS Chem. Biol. 2013 8 2377 2382 10.1021/cb400384c 24015959 54. Oja T Structural basis for C-ribosylation in the alnumycin A biosynthetic pathway Proc. Natl Acad. Sci. 2013 110 1291 1296 10.1073/pnas.1207407110 23297194 55. Burgess K Cook D Syntheses of nucleoside triphosphates Chem. Rev. 2000 100 2047 2060 10.1021/cr990045m 11749283 56. Prasad AK Trikha S Parmar VS Nucleoside synthesis mediated by glycosyl transferring enzymes Bioorg. Chem. 1999 27 135 154 10.1006/bioo.1998.1127 57. Hennen WJ Wong CH Design of an in vitro biocatalytic cascade for the manufacture of islatravir Science 2019 366 1255 1259 10.1126/science.aay8484 31806816 58. Kaspar F General principles for yield optimization of nucleoside phosphorylase-catalyzed transglycosylations ChemBioChem 2020 21 1428 1432 10.1002/cbic.201900740 31820837 59. Yehia H Efficient biocatalytic synthesis of dihalogenated purine nucleoside analogues applying thermodynamic calculations Molecules 2020 25 934 10.3390/molecules25040934 60. Vichier-Guerre S Ku TC Pochet S Seley-Radtke KL An expedient synthesis of flexible nucleosides through enzymatic glycosylation of proximal and distal fleximer bases ChemBioChem 2020 21 1412 1417 10.1002/cbic.201900714 31899839 61. Preumont A Snoussi K Stroobant V Collet J-F Van Schaftingen E Molecular identification of pseudouridine-metabolizing enzymes J. Biol. Chem. 2008 283 25238 25246 10.1074/jbc.M804122200 18591240 62. Huang S Mahanta N Begley TP Ealick SE Pseudouridine monophosphate glycosidase: a new glycosidase mechanism Biochemistry 2012 51 9245 9255 10.1021/bi3006829 23066817 63. Remsen JF Matsushita T Chirikjian JG Davis FF Enzymatic synthesis of deoxypseudouridylic acid and a study of certain of its properties Biochim. Biophys. Acta BBA-Nucleic Acids Protein Synth. 1972 281 481 487 10.1016/0005-2787(72)90148-7 64. Goldberg IH Rabinowitz M The enzymic synthesis of pseudouridine triphosphate Biochim. Biophys. Acta 1961 54 202 204 10.1016/0006-3002(61)90959-3 13899697 65. Goldberg IH Rabinowitz M Comparative utilization of pseudouridine triphosphate and uridine triphosphate by ribonucleic acid polymerase J. Biol. Chem. 1963 238 1793 1800 13948670 66. Frey, P. A. & Hegeman, A. D. Enzymatic Reaction Mechanisms. 25–30, 631–634 (Oxford University Press, 2007). 67. Heine A Observation of covalent intermediates in an enzyme mechanism at atomic resolution Science 2001 294 369 374 10.1126/science.1063601 11598300 68. Lehwess-Litzmann A Twisted schiff base intermediates and substrate locale revise transaldolase mechanism Nat. Chem. Biol. 2011 7 678 684 10.1038/nchembio.633 21857661 69. Lorentzen E Siebers B Hensel R Pohl E Mechanism of the schiff base forming fructose-1,6-bisphosphate aldolase: structural analysis of reaction intermediates Biochemistry 2005 44 4222 4229 10.1021/bi048192o 15766250 70. Allard J Grochulski P Sygusch J Covalent intermediate trapped in 2-keto-3-deoxy-6- phosphogluconate (KDPG) aldolase structure at 1.95-Å resolution Proc. Natl Acad. Sci. 2001 98 3679 3684 10.1073/pnas.071380898 11274385 71. Heine A Luz JG Wong C-H Wilson IA Analysis of the class I aldolase binding site architecture based on the crystal structure of 2-deoxyribose-5-phosphate aldolase at 0.99A resolution J. Mol. Biol. 2004 343 1019 1034 10.1016/j.jmb.2004.08.066 15476818 72. Fullerton SWB Mechanism of the class I KDPG aldolase Bioorg. Med. Chem. 2006 14 3002 3010 10.1016/j.bmc.2005.12.022 16403639 73. Pierce J Serianni AS Barker R Anomerization of furanose sugars and sugar phosphates J. Am. Chem. Soc. 1985 107 2448 2456 10.1021/ja00294a041 74. Chuvikovsky DV Ribokinase from E. coli : expression, purification, and substrate specificity Bioorg. Med. Chem. 2006 14 6327 6332 10.1016/j.bmc.2006.05.057 16784868 75. Chu CK El-kabbani FM Thompson BB Determination of the anomeric configuration of C-nucleosides by 1 H and 13 C NMR spectroscopy Nucleosides Nucleotides 1984 3 1 31 10.1080/07328318408079416 76. Wong CH Haynie SL Whitesides GM Preparation of a mixture of nucleoside triphosphates from yeast RNA: use in enzymic synthesis requiring nucleoside triphosphate regeneration and conversion to nucleoside diphosphate sugars J. Am. Chem. Soc. 1983 105 115 117 10.1021/ja00339a021 77. Rodrigues F Canac Y Lubineau A A convenient, one-step, synthesis of β-C-glycosidic ketones in aqueous media Chem. Commun. 2000 20 2049 2050 10.1039/b006642g 78. Bigi F Chesini L Maggi R Sartori G Montmorillonite KSF as an inorganic, water stable, and reusable catalyst for the knoevenagel synthesis of coumarin-3-carboxylic acids J. Org. Chem. 1999 64 1033 1035 10.1021/jo981794r 11674183 79. Mahrwald R The long underestimated carbonyl function of carbohydrates–an organocatalyzed shot into carbohydrate chemistry Chem. Commun. 2015 51 13868 13877 10.1039/C5CC04386G 80. Villadsen K Martos-Maldonado MC Jensen KJ Thygesen MB Chemoselective reactions for the synthesis of glycoconjugates from unprotected carbohydrates Chembiochem. Eur. J. Chem. Biol. 2017 18 574 612 10.1002/cbic.201600582 81. Frantz M-C Divergent entry to C-glycosides from unprotected sugars Org. Lett. 2019 21 2684 2687 10.1021/acs.orglett.9b00666 30916969 82. Ranoux A Lemiègre L Benoit M Guégan J-P Benvegnu T Horner–Wadsworth–Emmons reaction of unprotected sugars in water or in the absence of any solvent: one‐step access to C‐glycoside amphiphiles Eur. J. Org. Chem. 2010 2010 1314 1323 10.1002/ejoc.200901286 83. Cavezza A Synthesis of pro-xylane: A new biologically active C-glycoside in aqueous media Bioorg. Med. Chem. Lett. 2009 19 845 849 10.1016/j.bmcl.2008.12.037 19135365 84. Clapés P Fessner W-D Sprenger GA Samland AK Recent progress in stereoselective synthesis with aldolases Curr. Opin. Chem. Biol. 2010 14 154 167 10.1016/j.cbpa.2009.11.029 20071212 85. Garrabou X Asymmetric self- and cross-aldol reactions of glycolaldehyde catalyzed by d-fructose-6-phosphate Aldolase Angew. Chem. Int. Ed. 2009 48 5521 5525 10.1002/anie.200902065 86. Güclü D Minimalist protein engineering of an aldolase provokes unprecedented substrate promiscuity ACS Catal. 2016 6 1848 1852 10.1021/acscatal.5b02805 87. Thapa K Oja T Metsä-Ketelä M Molecular evolution of the bacterial pseudouridine-5’-phosphate glycosidase protein family FEBS J. 2014 281 4439 4449 10.1111/febs.12950 25081867 88. Wildberger P Pfeiffer M Brecker L Nidetzky B Diastereoselective synthesis of glycosyl phosphates by using a phosphorylase–phosphatase combination catalyst Angew. Chem. Int. Ed. 2015 54 15867 15871 10.1002/anie.201507710 89. Saheki S Takeda A Shimazu T Assay of inorganic phosphate in the mild pH range, suitable for measurement of glycogen phosphorylase activity Anal. Biochem. 1985 148 277 281 10.1016/0003-2697(85)90229-5 4061809 90. Morris GM AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility J. Comput. Chem. 2009 30 2785 2791 10.1002/jcc.21256 19399780 91. Wang L-P Building a more predictive protein force field: a systematic and reproducible route to AMBER-FB15 J. Phys. Chem. B 2017 121 4023 4039 10.1021/acs.jpcb.7b02320 28306259 92. Grade Web Server. http://grade.globalphasing.org/cgi-bin/grade/server.cgi. (2020).