Substrate-binding loop interactions with pseudouridine trigger conformational changes that promote catalytic efficiency of pseudouridine kinase PUKI

Pseudouridine, one major RNA modification, is catabolized into uracil and ribose-5′-phosphate by two sequential enzymatic reactions. In the first step, pseudouridine kinase (PUKI) phosphorylates pseudouridine to pseudouridine 5′-monophosphate. High-fidelity catalysis of pseudouridine by PUKI prevents possible disturbance of in vivo pyrimidine homeostasis. However, the molecular basis of how PUKI selectively phosphorylates pseudouridine over uridine with >100-fold greater efficiency despite minor differences in their Km values has not been elucidated. To investigate this selectivity, in this study we determined the structures of PUKI from Escherichia coli strain B (EcPUKI) in various ligation states. The structure of EcPUKI was determined to be similar to PUKI from Arabidopsis thaliana, including an α/β core domain and β-stranded small domain, with dimerization occurring via the β-stranded small domain. In a binary complex, we show that Ser30 in the substrate-binding loop of the small domain mediates interactions with the hallmark N1 atom of pseudouridine nucleobase, causing conformational changes in its quaternary structure. Kinetic and fluorescence spectroscopic analyses also showed that the Ser30-mediated interaction is a prerequisite for conformational changes and subsequent catalysis by EcPUKI. Furthermore, S30A mutation or EcPUKI complexed with other nucleosides homologous to pseudouridine but lacking the pseudouridine-specific N1 atom did not induce such conformational changes, demonstrating the catalytic significance of the proposed Ser30-mediated interaction. These analyses provide structural and functional evidence for a pseudouridine-dependent conformational change of EcPUKI and its functional linkage to catalysis.

RNA modifications have been recognized in noncoding RNAs and tRNAs and are also present in various coding RNAs (1,2). They have important roles in the maintenance of stability for tRNAs (3)(4)(5) and mRNAs (5) and in the regulation of the translation efficiency of mRNAs (6,7). Their biological effects are diverse. For example, RNA modifications in humans and yeast are regulated in response to changes in the environment (8,9). Among more than 140 RNA modifications identified in all organisms, N 6 -methylated adenine and pseudouridine (Ψ), a C5-glycoside isomer of uridine (Fig. 1A), are the most prevalent (10,11). The biogenesis of these modifications and their biological effects are reviewed elsewhere (1,2,5).
The catabolic fates of N 6 -methylated adenine (12) and pseudouridine are emerging from studies of noncanonical nucleotides derived from modification-containing RNAs. In particular, pseudouridine catabolism was first suggested in pyrimidine auxotrophic Escherichia coli mutants (13). Thirtyeight years later, two enzymes in the pathway-YeiC and YeiN-were identified (14). Details regarding pseudouridine catabolism, with its physiological consequences, were gleaned from studies of Arabidopsis (15). A lack of high sequence conservation from E. coli to eukaryotes hampered the identification of genes involved in pseudouridine catabolism within Arabidopsis thaliana. The metabolic turnover of pseudouridine-containing RNAs in vacuole produces pseudouridine monophosphate, which is dephosphorylated and exported to the peroxisome via the cytosol. In the peroxisome, pseudouridine is converted into uracil and ribose-5 0 -phosphate by two consecutive reactions that correspond to the activities of YeiC and YeiN, respectively, from E. coli. Therefore, the pseudouridine catabolic pathway recycles metabolic intermediates of uracil and ribose-5 0 -phosphate for use in other metabolic pathways.
Pseudouridine kinase (PUKI), a member of the phosphofructokinase B (PfkB) family (16), catalyzes the first of the two consecutive reactions (15). PUKI phosphorylates pseudouridine to pseudouridine 5 0 -monophosphate in an ATP-dependent manner (Fig. 1A); in the second step, pseudouridine 5 0 -monophosphate is hydrolyzed into uracil and ribose-5 0 -phosphate by pseudouridine monophosphate glycosylase. Recently, we performed a structural analysis of PUKI from A. thaliana (AtPUKI) and disclosed the structural basis for its substrate specificity toward pseudouridine (17). The enzyme efficiency (i.e., k cat /K m ) was 100-fold greater for pseudouridine than for structurally homologous uridine, despite minor differences in their K m values. These observations led us to postulate that AtPUKI catalytically distinguishes pseudouridine from uridine possibly by pseudouridineinduced conformational changes.
Here, we conducted a structural analysis of PUKI from E. coli strain B (EcPUKI) in the absence or the presence of pseudouridine, uridine, or cytidine. We also performed kinetic and fluorescence spectroscopic analyses. These findings provide structural and functional evidence concerning the pseudouridine-dependent conformational change in PUKI and its functional linkage to catalysis.

Phylogenetic analysis of PUKIs in the PfkB proteins
The sequence identity of full-length EcPUKI and AtPUKI is low (21%). We therefore investigated the evolutionary relationship of PUKIs in PfkB proteins using local structural information for AtPUKI (17). Before the construction of a phylogenetic tree, AtPUKI homologs were identified using the BLAST (18) by querying the full-length AtPUKI. The top 20 proteins ranked by similarity score were subjected to homology modeling in SWISS-MODEL (19) using AtPUKI (Protein Data Bank [PDB] ID: 7C1X) as the template. This modeling allowed us to identify the ligand-binding region with two hydrophilic active site residues; those residues correspond to positions 137 and 166 in AtPUKI, which are crucial for pseudouridine affinity (17). Multiple sequence alignment (MSA) was performed using amino acid sequences containing the two active site residues (20). This demonstrated the conservation of N (Asn) and K (Lys) in EcPUKI and AtPUKI (Fig. 1B). Notably, phylogenetic analysis revealed that EcPUKI and AtPUKI, together with E. coli YeiI, whose function is unknown, are evolutionarily closer than other PfkB proteins, such as ribokinases (RBSKs) (Fig. 1C). Pairwise distances for the alignments confirmed that functionally equivalent PfkB family proteins (ortholog) exhibit high sequence identities (Fig. 1D). The active site regions of PUKI were conserved between E. coli and A. thaliana, although the full sequences were phylogenetically distant. Therefore, EcPUKI and AtPUKI are homologous based on the similarity of two active site residues and the local structure.

Structure of EcPUKI in the absence of substrate
In the crystal structure of unliganded EcPUKI at 2.15 Å resolution, the asymmetric unit of the space group P3 has two dimers ( Table 1). The dimerization of EcPUKI in solution was Figure 1. Reaction scheme and phylogenetic comparison of PUKIs. A, in Escherichia coli, PUKI catalyzes the phosphorylation of pseudouridine using ATP as a phosphate donor, yielding pseudouridine 5 0 -monophosphate and ADP. B, MSA of local amino acid sequences including two hydrophilic active site residues for the top 20 highly homologous proteins to AtPUKI. Q94AT3_ARATH and PSUK_ECOLI refer to AtPUKI and EcPUKI, respectively. The function of E. coli YeiI is uncharacterized. Other homologs are described in the Experimental procedures section. C, phylogenetic tree of homologs of MSA constructed using the neighbor-joining method. D, heat map of pairwise distances of MSA. AtPUKI, PUKI from Arabidopsis thaliana; EcPUKI, PUKI from Escherichia coli strain B; MSA, multiple sequence alignment; PUKI, pseudouridine kinase. characterized by size-exclusion chromatography ( Fig. 2A). Therefore, the dimer is a biological functional unit. From a structural perspective, the four protomers in the two dimers are essentially identical, with an RMSD of 0.19 to 0.21 Å for about 307 Cα atoms.
These structural features are comparable to the features of AtPUKI, including the substrate-binding loop. A structural homology search using program DALI (22) indicated that a monomeric EcPUKI exhibited high similarity to AtPUKI (PDB ID: 7C1Y; Z-score, 39.2; sequence identity, 21%; and RMSD of 1.76 Å for the 267 Cα atoms) (17), EcRBSK (PDB ID: 1RKA; Zscore, 39.4; sequence identity, 22%; and RMSD of 1.57 Å for the 281 Cα atoms) (23), and RBSK from A. thaliana (PDB ID: 6ILT; Z-score, 38.7; sequence identity, 21%; and RMSD of 1.80 Å for the 290 Cα atoms) (24). In a previous structural and kinetic study of AtPUKI (17), we proposed that pseudouridineinduced conformational changes are required for catalysis. Those changes in conformation likely cause two substrates, pseudouridine and ATP, to be properly oriented for a kinase reaction. However, the proposed structural features were not unraveled in a ternary complex of AtPUKI with pseudouridine and ADP; only minor changes were identified (17), compared with the large conformational changes in EcRBSK (21,23). Analogous to EcRBSK and AtPUKI, we also adopted the Cαinteratomic distance as an indicator of conformational changes. Ser200, which is located at the edge of the core domain, was selected as the indicator residue (Fig. 2D). In unliganded EcPUKI, the Cα-interatomic distance between Ser200 and Ser200* (asterisk indicates a residue or an element from an adjacent subunit) in the two subunits was approximately 66 Å.
where I h is the observed intensity and <I h > is the average intensity. c The CC 1/2 is the Pearson correlation coefficient (CC) calculated from each subset containing a random half of the measurements of unique reflection. d R work = Σ||F obs | − |F cal ||/Σ|F obs |. e R free is the same as R obs for a selected subset (5%) of the reflections that was not included in prior refinement calculations. f R free of the EcPUKI-Ψ binary complex is relatively high because some of the monomers are suffered from poor density.

Pseudouridine-induced conformational changes
The 2.75 Å resolution crystal structure of EcPUKI complexed with pseudouridine provides structural evidence for conformational changes upon binding of pseudouridine (Fig. 3A). This binary structure, which has eight monomers in the asymmetric unit of the space group P2 1 (Table 1), exhibits two features. In each dimer, pseudouridine is present in only one subunit (gray in Fig. 3A). The other subunit shows a highly disordered structure in some regions. In particular, two pseudouridine-unbound subunits among the eight monomers had a disordered electron density. However, the Cα tracing in the other subunits was reliable. In an overall structure of the binary complex, the Cα-interatomic distance indicator was 51 Å (Fig. 3A), compared with 66 Å for unliganded EcPUKI, revealing global conformational changes upon binding of pseudouridine. These changes in quaternary structure are attributable to the β-stranded small domain, rather than the α/β core domain. The structural superposition of the corresponding subunit in unliganded EcPUKI and its binary complex resulted in an RMSD of 0.5 Å for 251 Cα atoms for the α/β core domain. Large structural differences are associated with the β-stranded small domain, with positional changes of up to approximately 9 Å in the substrate-binding loop between the corresponding Cα atoms (Fig. S2). Thus, the β-stranded small domain in a binary complex is more bent toward the α/β core domain, compared with that of unliganded EcPUKI. In dimeric EcPUKI, these dynamic features of the β-stranded small domain result in global conformational changes upon binding of pseudouridine, effectively sealing off pseudouridine for catalysis. 2). The elution peak for EcPUKI indicates a molecular weight of approximately 75 kDa, corresponding to dimeric EcPUKI. Left inset, 6.5 to 158 kDa molecular weight markers (Cytiva). Right inset, SDS-PAGE of monomeric EcPUKI. B, monomeric structure of unliganded EcPUKI is shown. An α/β core domain is shown, with the central β-sheet in yellow, flanking αhelices in gray, and the β-stranded small domain in orange. The monovalent cation-binding site is indicated by a green circle with the letter M. The secondary structural elements in Fig. S1 are labeled. C, topology of EcPUKI. The color codes and secondary structural elements are identical to those of B. D, side view of dimeric unliganded EcPUKI with a surface representation. The β-stranded small domain of each subunit mediates dimerization; a substratebinding loop (green) of the small domain transverses into an adjacent subunit (gray). The Cα-interatomic distance between Ser200 in each monomer was approximately 66 Å. EcPUKI, pseudouridine kinase from Escherichia coli strain B.
Pseudouridine in the binary complex is located in a cleft between the β-stranded small domain and the concave side of the α/β core domain (Fig. 3, B and C), with its nucleobase (hereafter referred to as uracil-Ψ) pointing toward the substrate-binding loop (green in Fig. 3B). Specifically, the ribosyl moiety forms several interactions within 3.3 Å with nearby residues: Asp16 with O2 0 and O3'; Asn14 with O2'; and Asp256, a proposed catalytic base, with O5 0 , to which the γ-phosphate group of ATP is attached. The uracil-Ψ moiety is surrounded by Tyr97, Asn112, and Met114 on one side and by the hydrophilic residues Ser30*, Asn143, and Lys170 (Fig. 3B). There are possible hydrogen bonds within 3.5 Å around the uracil-Ψ moiety: N3 and O4 of the base with Asn143 via bidentate interactions, O2 with Lys170, and N1 with Ser30*, a residue of the substrate-binding loop of the adjacent subunit.
Considering crystal packing interactions could affect the conformational state(s) of EcPUKI, we analyzed those packing interactions of unliganded EcPUKI in the P3 space group and in a binary complex with pseudouridine in the P2 1 space group (Table 1). In the PISA (Protein Interfaces, Surfaces and Assemblies) analysis (25), the interface area for crystal packing around a particular monomer was calculated. Although the two space groups had similar interface area values (2894 Å 2 for P3 and 2737 Å 2 for P2 1 ), a small number of crystal contacts in the critical regions could affect the conformational freedom of molecules. Therefore, we still cannot rule out the possible effects of crystal contacts on the changes in EcPUKI conformation. Consistent with the crystal structure of EcPUKI complexed with pseudouridine, further structural analysis, mutagenesis, and fluorescence spectroscopy of various EcPUKIs in the current study unanimously support the pseudouridine-induced conformational changes of EcPUKI.

Catalysis of EcPUKI
The binding environment of pseudouridine in EcPUKI is comparable to that of pseudouridine in AtPUKI (17). In AtPUKI, there are variations in the hydrophobic residues around the uracil-Ψ moiety, but three hydrophilic residues were identified: Asn137 AtPUKI , Lys166 AtPUKI , and Thr26 AtPUKI (superscript AtPUKI indicates residue(s) in AtPUKI) (Fig. 1B). Three residues, which differ in Thr26 AtPUKI and Ser30 EcPUKI , are hallmarks of PUKI family proteins.
EcPUKI residues interacting with a uracil-Ψ have crucial roles in enzyme activity. For a functional analysis, we validated that the monovalent cation K + , the divalent cation Mg 2+ , and a phosphate donor ATP are required for EcPUKI activity, with high specificity for pseudouridine and no requirement for inorganic phosphate as an activator (Fig. 4, A-E). Specific activity measurements of mutants demonstrated that residues around pseudouridine have various effects on enzyme activity (Figs. 3B and 5A). Among them, the mutation of three hydrophilic residues around the uracil-Ψ moiety-Ser30, Asn143, and Lys170-resulted in catalytically incompetent enzymes, with activities of 0.1 to 5.5% of the wildtype enzyme (Fig. 5A). Therefore, the possible hydrogen bonds mediated by these three residues have central roles in the kinase activity of EcPUKI. Steady-state kinetic analyses yielded supporting evidence (Fig. 5, B and C). Mutation of the uracil-Ψ-interacting residues (Y97A, N112A, and M114A) and the three hydrophilic residues (S30A, N143A, and K170A) increased the K m value by approximately 1.4-to 55.1-fold, and their k cat values are in a range of 0.087 to 109%, compared with wildtype EcPUKI. The enzyme efficiency (i.e., k cat /K m ) of the mutants was therefore 0.06 to 31% of wildtype EcPUKI. Mutant S30A showed the lowest efficiency for pseudouridine (0.06% of wildtype EcPUKI). Unlike the N143A and K170A mutants showing 55and 11-fold increased K m values, respectively, the S30A mutant had a significantly reduced k cat value (0.087% of wildtype EcPUKI), but its K m value is only a 1.4-fold increase.

Functional role of Ser30 in the substrate-binding loop
The substrate-binding loop containing Ser30 has dynamic structural features that enable the pseudouridine-induced Activities were compared with wildtype EcPUKI (48.2 μmol min −1 mg −1 protein set to 100%; Fig. 5A). Consistent with PUKI of E. coli K-12 (14), EcPUKI showed high specificity for pseudouridine, with negligible activities for other potential substrates, including similar pyrimidine nucleosides (uridine, cytidine, 5 0methyluridine, thymidine, and 2 0 -deoxycytidine). EcPUKI, pseudouridine kinase from Escherichia coli strain B; n.d., not detected. Figure 5. Mutational and kinetic analyses of EcPUKIs. A, specific activity of EcPUKIs. Red, wildtype; blue, pseudouridine-binding pocket mutants; and yellow, mutants for possible catalytic residues. The specific activities in parentheses were compared with wildtype EcPUKI (48.2 μmol min −1 mg −1 protein set to 100%). Each measurement was conducted in triplicate; error bars indicate standard deviations. Activity assays were performed as described (17). ADP production was monitored for 30 s under substrate saturation conditions (4 mM ATP and 1.25 mM pseudouridine) for wildtype EcPUKI. B, steady-state kinetic analyses of EcPUKI with pseudouridine and ATP as substrates. Error bars are standard deviations (n = 3). Kinetic analyses involved 4 mM ATP or 1.25 mM pseudouridine with 100 nM EcPUKI. C-E, steady-state kinetic parameters of the mutants. Kinetic assays of phosphate acceptors (pseudouridine, uridine, and cytidine) were conducted in the presence of 4 mM ATP. The kinetic assay for ATP in Figure 5C was performed in the presence of 1.25 mM pseudouridine. Concentrations of EcPUKI are in a range of 100 nM-24 μM. EcPUKI, pseudouridine kinase from Escherichia coli strain B; n.d., not detected.
conformational changes (Figs. 3A and S2). Considering the significant reduction of the k cat , but not the K m , value of the S30A mutant (Fig. 5C), we propose that Ser30-mediated interaction between a substrate-binding loop and uracil-Ψ constitutes a prerequisite for conformational changes and catalysis by EcPUKI. Therefore, the effects of the S30A mutation likely result from an inability to undergo the conformational changes required for activity.
We performed a structural analysis of EcPUKI to validate the proposed function of Ser30. Indeed, the 1.90 Å resolution crystal structure of S30A complexed with pseudouridine did not show global conformational changes ( Fig. 3D and Table 1; distance indicator of 65 Å, similar to unliganded EcPUKI). Pseudouridine is present in both subunits, unlike the wildtype enzyme and its binding environment in the S30A mutant is reminiscent of that of wildtype EcPUKI (Fig. 3, E and F). However, because of the Ser30 to Ala30 mutation, there is no interaction within 7 Å between Ala30 and the N1 of uracil-Ψ, restricting movement of the β-stranded small domain and its adoption of an orientation identical to unliganded EcPUKI. Next, we determined the structure of wildtype EcPUKI complexed with uridine or cytidine ( Fig. S3 and Table 1) at 2.15 and 2.20 Å resolution, respectively, which do not interact with Ser30 because of the absence of the hallmark uracil-Ψ N1 atom in the nucleobases (Fig. 5, D and E). Those ligands, which are structurally homologous to pseudouridine, bound within a pseudouridine-binding cleft in an almost identical manner to wildtype EcPUKI complexed with pseudouridine. However, they did not induce conformational changes (Fig. S3). These structural features are correlated with the kinetic parameters of wildtype EcPUKI and S30A mutant for uridine and cytidine (Fig. 5, D and E), with a 1500-to 7000-fold reduction in k cat and a 4-to 14-fold increase in K m , compared with wildtype EcPUKI for pseudouridine. Therefore, the conformational changes in EcPUKI induced by the binding of pseudouridine to a substrate-binding loop via Ser30 are essential for activity.

Conformational changes of EcPUKI in solution
We further investigated substrate-induced conformational changes of EcPUKI in solution using intrinsic fluorescence measurements. Using an excitation wavelength of 289 nm, the emission spectra of EcPUKI in the presence of various ligands showed that intrinsic fluorescence at 343 nm was enhanced by up to 35% by an authentic phosphate acceptor pseudouridine; it was not enhanced by uridine, cytidine, or the authentic phosphate donor ATP (Fig. 6A). Considering that pseudouridine is not fluorescent under the experimental conditions (Fig. 6B), the enhanced fluorescence emission is related to changes in EcPUKI associated with the binding of pseudouridine.
Subsequent experiments indicated that the enhanced fluorescence emission was caused by pseudouridine-induced conformational changes. In the structural transition between unliganded EcPUKI and a binary complex with pseudouridine, the side chain of Trp169 between the small and core domains undergoes large positional displacements, with notable changes in its binding environments (Fig. S4, A-C). The side chain of Trp169, with disordered density, is in a fully exposed orientation in the unliganded state, although it is buried in the binary complex by conformational changes. Considering that the fluorescence of tryptophan is affected by its binding environment, we evaluated whether Trp169 could serve as a fluorescent probe for conformational changes in EcPUKI. Using a W169A mutant, pseudouridine-dependent fluorescence at 343 nm was enhanced only by 13%, compared with the unliganded W169A mutant (Fig. 6C). The W169A mutant is catalytically competent; its specific activity and kinetic parameters for pseudouridine are comparable with wildtype EcPUKI (Figs. 5, A and C, S4D). Therefore, Trp169 is a major contributor to the enhanced pseudouridine-dependent fluorescence that is caused by pseudouridine-dependent conformational changes. Pseudouridine did not affect the fluorescence of the S30A mutant (Fig. 6D). Consistent with the structural analysis, these findings suggest that the S30A mutant does not undergo conformational changes in solution. Therefore, the pseudouridine-induced conformational changes are required for catalysis by EcPUKI in solution.

Ligand-induced conformational changes and catalysis
The EcPUKI-dependent kinase activity requires ATP as a phosphate donor and pseudouridine as a phosphate acceptor, with cofactors of monovalent and divalent cations. A proposed hydrogen bond between the side-chain hydroxyl group of Ser30 and the N1 of uracil-Ψ in pseudouridine triggers conformational changes via the β-stranded small domain.
In a PUKI-dependent kinase reaction, a catalytic base (e.g., Asp256 EcPUKI ) deprotonates O5 0 of the ribosyl moiety in pseudouridine (Figs. 3B and S5A). The resulting deprotonated O5 0 of the ribosyl moiety performs a nucleophilic attack on the γ-phosphate group of ATP, thereby producing pseudouridine 5 0 -monophosphate and ADP. Liberation of the terminal γphosphate of ATP requires a catalytic acid. Structural and functional analyses of EcPUKI (this study), AtPUKI (17), and RBSK (21,24,26,27) indicated that the aspartate residue in the immediate vicinity of the O5 0 of the ribosyl moiety constituted the catalytic base (Fig. S5, B-D). However, the catalytic acid is diverse: Lys291 in RBSK from A. thaliana (Fig. S5B) (24) and possibly a water molecule in the Mg 2+ -coordinating shell near ATP in AtPUKI (Fig. S5C) (17). Our activity and structural analyses identified Asp256 as the catalytic base (Figs. 3B, 5, A and C). Consistent with a catalytic role for Asp256, the D256A mutant showed a 3280-fold decreased k cat and 1.4-fold increased K m for ATP (Fig. 5C). Moreover, Lys185 in EcPUKI could be located in the vicinity of the γ-phosphate of ATP (Fig. S5E) and be crucial for activity (Fig. 5A). However, a role for Lys185 as a catalytic acid is controversial; the K185A variant showed a 54-fold increased K m and 100-fold decreased k cat for ATP (Fig. 5C).
We were unable to evaluate the effects of conformational changes on catalysis because a ternary complex of EcPUKI with ATP and pseudouridine is unavailable. However, fluorescence measurements suggested that conformational changes require pseudouridine; ATP was insufficient (Fig. 6A). A structural comparison indicated that the conformational changes do not affect the structural integrity of the ATP-binding site in EcPUKI (Fig. S6). The function of pseudouridine-induced conformational changes is unclear; it could include orienting the O5 0 of the ribosyl moiety of pseudouridine and the γ-phosphate group of ATP in a position suitable for nucleophilic attack, enhancing affinity for ATP, or providing the catalytic acid. Our structural and kinetic analyses suggest a functional linkage between pseudouridine-dependent conformational changes and catalysis.

Conclusion
We provide structural evidence for pseudouridine-dependent conformational changes of dimeric EcPUKI. Kinetic and fluorescence analyses further indicated that those conformational changes are required for the activity of EcPUKI. Ser30 in the substrate-binding loop of EcPUKI forms a hydrogen bond with a pseudouridine nucleobase, triggering pseudouridinedependent conformational changes. Therefore, a serine (or threonine) residue in the substrate-binding loop is conserved among PUKI family proteins, as are two key asparagine and lysine residues, which mediate hydrogen bonding with a pseudouridine nucleobase. The presence of the active site residues Asn143 and Lys170 for maintaining the substrate affinity and Ser30 for inducing conformational changes for catalysis provides mechanistic insight into the high fidelity of PUKI for pseudouridine, without disturbing pyrimidine homeostasis.

Experimental procedures
Phylogenetic analysis BLASTP 2.9.0+ was used to identify homologs of AtPUKI in the UniProtKB/Swiss-Prot Database (https://www.ebi.ac.uk/ Tools/sss/ncbiblast/) (18,28). To identify regions with two active site residues in the homolog proteins, homology modeling was performed using SWISS-MODEL (https://swissmodel.expasy.org/) (19). MSA was conducted for the amino acid sequences of the regions using the R package "msa" (20). Next, the squared root of the pairwise distance was calculated by the method "identity" of the R package "seqinr" (29). Finally, clustering for phylogenetic tree construction was performed by neighbor-joining method implemented in the R package "ape" (30).

Cloning and purification of EcPUKI
The PUKI gene from E. coli strain B (National Center for Biotechnology Information reference sequence: WP_ 001208132.1) was amplified by polymerase chain reaction. The Pseudouridine-induced conformational changes resulting amplification product was cloned into a modified pET28b vector (Merck) that contained a tobacco etch virus protease cleavage site. The vector harboring the gene for EcPUKI and an N-terminal His-tag was transformed into E. coli BL21 (DE3) cells. The cells were grown at 37 C in Luria-Bertani medium until an absorbance of 0.5 at 600 nm. Recombinant EcPUKI was then induced with 0.5 mM IPTG for 16 h at 20 C. The E. coli cells were harvested, suspended, and sonicated in buffer A (50 mM Tris-HCl at pH 8.0, 300 mM KCl, and 5% [w/v] glycerol). Recombinant EcPUKI with an N-terminal His-tag was purified by immobilized metal affinity chromatography using a HisTrap HP column (GE Healthcare) equilibrated with buffer A and subsequently eluted with buffer A plus 0.5 M imidazole. The N-terminal His-tag was removed by treatment with tobacco etch virus protease overnight at 4 C. EcPUKI was further subjected to size-exclusion chromatography on a Superdex 200 column (GE Healthcare) equilibrated with buffer A.
The mutants for structural and functional analyses were constructed using mutagenic primers (Table S1). Vector construction and protein purification were performed as described previously. CD measurements of several enzymes, including wildtype EcPUKI and the S30A, N143A, W169A, and K170A mutants, indicate that mutations do not cause global structural perturbations (Fig. S7). To form a binary complex with pseudouridine, EcPUKI was purified with buffer B (50 mM potassium phosphate at pH 7.4, 300 mM KCl, and 5% [w/v] glycerol).

Crystallization and structure determination
Tag-free EcPUKI was concentrated to 15 mg/ml and then subjected to crystallization using the sitting drop vapordiffusion method at 22 C. Initial crystal screening was performed using the commercially available kits. Crystals of unliganded EcPUKI in buffer A were produced in 0.2 M magnesium formate, 18 to 22% PEG 3350, and 5% glycerol. The crystallization of a possible ternary complex was carried out, with ADP and pseudouridine, uridine, or cytidine as a phosphate acceptor. Specifically, EcPUKI in buffer B was reacted at 4 C with 2 mM ADP and 2 mM ligand of interest for 1 h and then subjected to crystallization. However, these crystallization conditions resulted in a structure with a phosphate acceptor but not ADP. Therefore, we considered these conditions to be sufficient for the binary complex. The crystallization of wildtype EcPUKI with uridine or cytidine and the S30A mutant with pseudouridine was successful under conditions identical to the conditions for the unliganded EcPUKI. However, crystals of wildtype EcPUKI with pseudouridine were produced in 1.0 M lithium chloride, 0.1 M MES (pH 6.0), and 18% PEG 6000. Prior to data collection, crystals for the binary complex with various ligand(s) were soaked in crystallization solution plus 4 mM ADP and 4 mM phosphate acceptor. Even with these soaking experiments, the enzyme formed a binary complex only with a phosphate acceptor. Cryoprotection was achieved using 20% glycerol plus crystallization solution.
Diffraction data collection was conducted at 100 K with a 0.5 oscillation angle on beamlines 7A and 11C at the Pohang Accelerator Laboratory. The collected data were indexed, integrated, and scaled using HKL2000 program (31). The highresolution cutoff was based on a CC 1/2 statistical value of approximately 0.4 (32,33). All crystals belonged to the space group P3, with four monomers in the asymmetric unit, except for the EcPUKI-pseudouridine complex, which contained eight monomers in an asymmetric unit of the space group P2 1 .
The structure of unliganded wildtype EcPUKI was solved by molecular replacement using Phaser in PHENIX program (34), with the structure of a putative sugar kinase from Clostridium perfringens (PDB ID: 3KZH; sequence identity, 34%) as a search model. The initial model was manually rebuilt and refined using COOT (35) and PHENIX, respectively. Initial refinement by PHENIX yielded a model with a high R work and an R free value of >30%. At this stage, only residues or regions of the model with a reasonable electron density were included into the working structure. Subsequent cycles of PHENIX with anisotropic B-factor refinement were performed with incorporation of more residues. Cycles of refinement and extensive model rebuilding improved the structure by reducing the R work and R free values. At that stage, the refined model of unliganded EcPUKI was used as a search model for the molecular replacement of binary complexes with uridine, cytidine, or pseudouridine. In the middle of refinement, uridine, cytidine, and pseudouridine were present in the structure based on the Fo-Fc map of the active site. The data collection and refinement statistics are provided in Table 1. In the EcPUKI-pseudouridine complex, two of eight monomers had highly disordered features, possibly leading to a high R free value (Table 1). MolProbity (36) in PHENIX program was used for structure validation.

Activity assay
N-terminal His-tagged enzymes were purified in buffer A by immobilized metal affinity chromatography as described previously and subjected to desalting on a HiPrep 26/10 (GE Healthcare) with buffer A.
Two different assays, a direct assay (37) and an enzymecoupled assay (14,24,38), were employed to measure EcPUKI activity, as reported previously for the assays of AtPUKI (17). In the direct assay, the reaction mixture contained 40 mM Tris-HCl (pH 7.5), 20 mM divalent cation, 50 mM monovalent cation, 0.003% phenol red, 4 mM ATP, and 95 nM wildtype EcPUKI. After 2 min of incubation at 25 C, the enzyme reaction was triggered by adding 1.25 mM pseudouridine, and absorbance at 430 nm was monitored for 30 s. The monovalent cations tested were LiCl, NaCl, KCl, RbCl, and CsCl, and the divalent cations were MgCl 2 , MnCl 2 , and CaCl 2 . The phosphate donor candidates were ATP, GTP, CTP, and UTP at 1 mM. All assays were conducted in triplicate.
In the enzyme-coupled assay, the reaction mixture consisted of 40 mM Tris-HCl (pH 7.5), 20 mM MgCl 2 , 50 mM KCl, 0.2 mM NADH, 2 mM phosphoenolpyruvate, the indicated concentrations of EcPUKI, 18 to 28 units of lactate dehydrogenase and 12 to 20 units of pyruvate kinase per milliliter, and pseudouridine or ATP as the substrates. After 2 min of incubation at 25 C, the enzyme reaction was initiated by adding the remaining substrate, and the decrease in absorbance at 340 nm was monitored. The initial velocity was determined between 15 and 45 s and expressed as the corresponding ADP concentration change per minute, using an extinction coefficient of 6220 M −1 cm −1 at 340 nm for NADH. Sigmaplot (Systat Software) was used to calculate the K m and V max .

Fluorescence measurement
Fluorescence spectra were collected at room temperature from 1.0 mg/ml N-terminal His-tagged EcPUKI and mutants using the QuantaMaster 400 (Horiba-PTI). The emission spectra of EcPUKIs were obtained at an excitation wavelength of 289 nm. For fluorescence measurement in the presence of ligand(s), EcPUKI was incubated at 4 C for 4 h with 312.5 μM nucleoside (pseudouridine, uridine, or cytidine) or 750 μM ATP.

Data availability
The atomic coordinates and structural factors have been deposited in the PDB (http://www.rcsb.org) under ID code 7VTD, 7VTE, 7VTF, 7VTG, and 7VVA.