Comparison of Spatial Structures and Packaging of Phosphorybosil Pyrophosphate Synthetase 2 from Thermus thermophilus HB27 in Rhombohedral and Tetragonal Crystals

: We report the spatial structure of phosphoribosyl pyrophosphate synthetase 2 from the thermophilic bacterium Thermus thermophilus HB27 (TthPRPPS2) obtained at a 1.85 Å resolution using a diffraction set collected from rhombohedral crystals (space group R3 2 -h ), grown with lithium sulfate as a precipitant. This crystal structure was compared with the structure of TthPRPPS2, previously obtained at a 2.2 Å resolution using diffraction sets from the tetragonal crystals (space group P4 1 2 1 2 ), grown with ammonium sulfate as a precipitant. The comparison of these structures allows the study of the differences between protein molecules in both crystalline structures, as well as the packaging of enzyme molecules in crystals of both spatial groups. Our results may contribute to the research of the structural basis of catalytic activity and substrate speciﬁcity of this enzyme.


Introduction
Phosphoribosyl pyrophosphate synthetases (PRPPS) catalyze the synthesis of 5-phosphorybosyl-α-1-pyrophosphate (5PRPP), using ribose-5-phosphate (R5P) and ATP as substrates. Enzymes of the PRPPS family are widespread in prokaryotic and eukaryotic organisms [1,2]. The functioning of PRPPS accompanies the most important metabolic processes and is vitally important for a living cell. These enzymes are involved in metabolic processes, including the synthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, and pyridine coenzymes. PRPP synthetases are the objects of a number of our studies in the multienzymatic cascade conversion of D-pentoses into purine nucleotides [3,4]. One of the ways to improve the cascade and expand the possibilities of its utilization is the use of enzymes with altered substrate specificity.
As mentioned earlier, the product of the PRPPS reaction is PRPP. This molecule is mainly involved in the de novo and salvage biosynthesis of purine and pyrimidine nucleotides and acts as a regulator of both processes [5][6][7]. The significant role of PRPPS in cellular metabolism requires strict regulation of its activity. The enzyme activity is controlled by two mechanisms: feedback inhibition by excess substrate and allosteric inhibition. The molecule responsible for the latter inhibition mechanism is adenosine diphosphate (ADP), which binds with the allosteric center of the enzyme [8][9][10].
There are three classes of PRPP synthetases, which are determined by the type of their quaternary structure, the mechanism of regulatory activity, and specificity with respect to the pyrophosphate residue donor [11][12][13]. The object of research is the thermophilic enzyme PRPPS2 from Thermus thermophilus HB27 (TthPRPPS2). This enzyme belongs to the first-class PRPPS. Enzymes of this subfamily have a hexameric quaternary structure, and their activity depends on phosphate ions and is regulated by allosterically binding ADP or GDF molecules in the allosteric center [14]. Only ATP can serve as the pyrophosphate donor for this class of PRPPS.
Previously, we have obtained the structure of TthPRPPS2 at a resolution of 2.2 Å using protein crystals grown with ammonium sulfate at pH 4.9 as a precipitant [15]. The crystals belonged to the space group P4 1 2 1 2.
In this study, we used a protein solution of the same composition as when growing tetragonal crystals but with lithium sulfate as a precipitant at pH 7.5. This resulted in the growth of the rhombohedral crystals of TthPRPPS2, belonging to the R3 2 h space group. A diffraction set was collected and used to solve the spatial structure of the enzyme at a resolution of 1.85 Å. This study compares the crystal structures of TthPRPPS2 molecules obtained at a resolution of 2.2 and 1.85 Å, and the packing of enzyme molecules in tetragonal and rhombohedral crystals.
The spatial structure of tetragonal crystals is described in our previous paper [15] (PDB ID: 5T3O). The spatial structures of both modifications were solved by the molecular replacement method using the Phaser program [19] and E. coli PRPPS coordinates (PDB ID: 1DKR) as a starting model [20], for tetragonal crystals (space group P4 1 2 1 2), and TthPRPPS2 coordinates (PDB ID: 5T3O) as a starting model for rhombohedral crystals [15]. The diffraction set from rhombohedral crystals was collected at a temperature of 100 K at the BL41 station of the SPring-8 synchrotron facility (Japan) at the resolution of 1.85 Å.
The crystal structure (space group R3 2 -h) was refined at the resolution of 1.85 Å. The refined coordinates (REFMAC program [21]) of the atomic model were deposited in the PDB, PDB_ID: 7PN0. In both structures, ADP molecules and sulfate ions bounded with the enzyme are localized on the electron density map. The statistics of the X-ray diffraction set and the refinement of the structure are given in Table 1. A fragment of the electron density map is shown in Figure S1. The 5T3O and 7PN0 structures were compared by superposing them with the PyMol program [22]. The analysis of structures, intramolecular, and intermolecular contacts in the crystal structures of TthPRPPS2 was conducted using the PISA web service [23].

Results and Discussion
The biologically active form of phosphoribosyl pyrophosphate synthetase 2 from Thermus thermophilus is homohexamer. The hexamer PRPRS molecule consists of subunits connected by a third-order axis and a second-order axis perpendicular to it and has a point symmetry group of 32.
Rhombohedral and tetragonal crystals of Tth2PRPPS HB27 differ in the number of subunits in the independent part of the unit cell. In the asymmetric unit, tetragonal crystals contain three subunits of a hexameric molecule, and rhombohedral crystals contain two subunits. In tetragonal crystals, the hexamer molecule is formed as a result of the action of the second-order crystallographic axis (the symmetry operator y, x, −z converts a trimer into a hexamer). In a rhombohedral crystal structure, the third-order symmetry axis of the hexameric molecule coincides with the crystallographic one. A hexamer molecule is formed by applying the symmetry operator −y−1, x−y, z to two subunits contained in an independent part of the unit cell. In both structures, the point group of pseudosymmetry of the molecule is 32.
Each subunit of the hexameric enzyme molecule consists of two domains: N terminal and C terminal ( Figure 1).
Both domains have a polypeptide chain folding characteristic of PRPPS class 1. The polypeptide chain in the N-terminal and C-terminal domains forms a five-stranded βstrand surrounded by four α-helices. The N-terminal domain additionally contains a 3 10 -helix and an α-helix at the C-terminal region. The central core in both domains is covered by a double antiparallel β-layer. These β-layers (turns β3-β4 of the N-terminal domain and β9-β10 of the C-terminal domain) are called flag regions [15]. The domains are bounded by the linker (amino acid residues 147-149). The subunits in the hexameric molecule are stacked similar to a propeller with N-terminal domains at the center of the molecule and C-terminal domains on the surface of the molecule. The active center of the subunit is located in the cavity between the N-terminal and the C-terminal domains. In the N-terminal domain, it is surrounded by the amino acid residues of the flexible loop (Phe92-Ser108) and a fragment of the polypeptide chain of the neighboring subunit (Val30-Ile44). In the C-terminal site, the active center is surrounded by a pyrophosphate-binding site (Asp170-Gly175) and so-called riboso-5-phosphate (R5P) loop (220-228). The flag region of the C-terminal domain containing residues Arg196 and Lys194, which are important for catalytic activity, called catalytic β-hairpin (residues 189-208, β9-β10), closes the entrance to the active site during the transition of the enzyme to the closed conformation [15]. In each subunit of the hexameric enzyme molecule, along with the active center, there is an allosteric center where an ADP molecule can bind [16]. The amino acid residues binding the ADP molecule belong to two adjacent trimer subunits and one of the subunits of the second trimer of the hexameric molecule. For example, tetragonal crystals with three independent subunits A, B, C, forming a trimer, the SO 4 ion located in the allosteric center of subunit B is coordinated by the amino acid residues of this subunit B, subunit A and subunit C of the second trimer. Both domains have a polypeptide chain folding characteristic of PRPPS class 1. The polypeptide chain in the N-terminal and C-terminal domains forms a five-stranded βstrand surrounded by four α-helices. The N-terminal domain additionally contains a 310helix and an α-helix at the C-terminal region. The central core in both domains is covered by a double antiparallel β-layer. These β-layers (turns β3-β4 of the N-terminal domain and β9-β10 of the C-terminal domain) are called flag regions [15]. The domains are bounded by the linker (amino acid residues 147-149). The subunits in the hexameric molecule are stacked similar to a propeller with N-terminal domains at the center of the molecule and C-terminal domains on the surface of the molecule. The active center of the subunit is located in the cavity between the N-terminal and the C-terminal domains. In the N-terminal domain, it is surrounded by the amino acid residues of the flexible loop (Phe92-Ser108) and a fragment of the polypeptide chain of the neighboring subunit (Val30-Ile44). In the C-terminal site, the active center is surrounded by a pyrophosphatebinding site (Asp170-Gly175) and so-called riboso-5-phosphate (R5P) loop (220-228). The flag region of the C-terminal domain containing residues Arg196 and Lys194, which are important for catalytic activity, called catalytic β-hairpin (residues 189-208, β9-β10), closes the entrance to the active site during the transition of the enzyme to the closed conformation [15]. In each subunit of the hexameric enzyme molecule, along with the active center, there is an allosteric center where an ADP molecule can bind [16]. The amino acid residues binding the ADP molecule belong to two adjacent trimer subunits and one of the subunits of the second trimer of the hexameric molecule. For example, tetragonal crystals with three independent subunits A, B, C, forming a trimer, the SO4 ion located in the allosteric center of subunit B is coordinated by the amino acid residues of this subunit B, subunit A and subunit C of the second trimer.
In both crystal polymorphs (tetrahedral and rhombohedral), an ADP molecule and a sulfate ion are localized in the active center of each enzyme subunit, and sulfate ions are located in the allosteric centers. An additional SO4 ion is localized in rhombohedral crys- In both crystal polymorphs (tetrahedral and rhombohedral), an ADP molecule and a sulfate ion are localized in the active center of each enzyme subunit, and sulfate ions are located in the allosteric centers. An additional SO 4 ion is localized in rhombohedral crystals, coordinated by amino acid residues Thr129 (loop 128-132), Ser147, Ala148 (irregular region 146-148), and Arg177 (α-helix 174-186). For the rest of the Cα atoms of the subunits, noticeable shifts are observed in the sequences 162-216, 231-244, and in the catalytic hairpin (residues 186-208). In the independent subunits of rhombohedral crystals, the positions of the catalytic hairpins differ and, as a result, the positions of the ligands bounded in these regions: sulfate ions and ADP molecules bounded in the active center. (Figure 2a,b). The positions of sulfate ions, localized in the active centers of independent subunits, differ by 2.6 Å.
Comparison of the position of the catalytic hairpins showed that in one of the independent subunits (subunit B), the conformation of the active center is more closed (Figure 2a). Thus, the two independent subunits of the molecule have a different conformation of the active center.
Earlier [15], when comparing the conformations of the subunits of the protein molecule in the tetragonal crystal modification, it was found that the catalytic hairpin bounding the active center region has different conformations in three independent TthPRPPS2 subunits, while the closed conformation of the active site was found in one of the subunits. Comparison of conformations of the catalytic loop (residues 188-202) in three independent subunits of tetragonal crystals [15] suggested that different conformations of the catalytic When subunits with a closed conformation of the active center from different crystal modifications (subunits B with closed conformation in rhombohedral crystals and subunit B with closed conformation in tetragonal crystals) were superimposed over 255 Cα atoms, the average rmsd was 0.369 Å (Figure 3). ent subunits of tetragonal crystals [15] suggested that different conformations of the catalytic loop are realized during the transition between the closed and open conformations of the active center.
When subunits with a closed conformation of the active center from different crystal modifications (subunits B with closed conformation in rhombohedral crystals and subunit B with closed conformation in tetragonal crystals) were superimposed over 255 Сα atoms, the average rmsd was 0.369 Å (Figure 3). The greatest differences were found in the C-terminal domain, where the positions of loops 160-165, 185-189, 208-215, 239-243 differ. All these sites are located on the surface of the molecule. There are also differences in the conformation of some residues of catalytic hairpins (β9-β10). Amino acid residues, which are part of the active center of subunits and involved in the coordination of the sulfate ion and the ADP molecule, have the same conformations in the compared subunits of different crystal structures. ADP molecules localized in the active centers of the compared subunits differ in the conformation of phosphate groups (Figure 3).
Intermolecular polar contacts in crystals of two modifications were studied using the PISA program ( Table 2). The surface area of enzyme molecules accessible to the solvent is almost the same in both crystalline modifications and amounts to 61040 Å 2 for tetragonal The greatest differences were found in the C-terminal domain, where the positions of loops 160-165, 185-189, 208-215, 239-243 differ. All these sites are located on the surface of the molecule. There are also differences in the conformation of some residues of catalytic hairpins (β9-β10). Amino acid residues, which are part of the active center of subunits and involved in the coordination of the sulfate ion and the ADP molecule, have the same conformations in the compared subunits of different crystal structures. ADP molecules localized in the active centers of the compared subunits differ in the conformation of phosphate groups (Figure 3).
Intermolecular polar contacts in crystals of two modifications were studied using the PISA program ( Table 2). The surface area of enzyme molecules accessible to the solvent is almost the same in both crystalline modifications and amounts to 61040 Å 2 for tetragonal crystals and 59280 Å 2 for rhombohedral crystals. However, the number of molecules in the immediate environment and the regions of the polypeptide chain involved in intermolecular interactions are different for both modifications.
As discussed earlier, a trimer is located in the independent part of the tetragonal crystal cell, and a biologically active hexamer molecule is formed under the action of a secondorder crystallographic axis (symmetry operator y, x, −z). In the tetragonal modification, intermolecular contacts involve molecules bounded by a twofold screw axis, which is parallel to the a axis at the 3/8 c level (symmetry operator: x−1/2, −y−1/2, −z−1/4) ( Table 2). Intermolecular contacts with four neighboring molecules form four subunits of a hexamer molecule: independent subunits B, C and subunits B*, C* of the second trimer of the molecule, obtained as a result of the action of the symmetry operator y, x, −z. Hexamer molecule and contacting subunits of other molecules in the P4 1 2 1 2 crystal modification are shown in Figure 4.
The total area of contact of one molecule with four neighboring molecules is 3936 Å 2 . In tetragonal crystals, polar intermolecular interactions involve amino acid residues belonging to the regular regions of the structure: helices α4, α5, η4, and strand β14 ( Table 2).
As mentioned, in rhombohedral crystals, a dimer is located in the independent part of the cell, and a disk-shaped hexamer molecule is formed under the action of a third-order crystallographic axis. In each trimer, the conformation of its subunits is the same, but the conformations of the subunits in different trimers of molecules are different. All six subunits of the hexamer are involved in intermolecular contacts in rhombohedral crystals. Intermolecular contacts are formed between trimers that have the same conformation of subunits in contacting molecules. The total number of neighboring molecules in contact with the molecule is six. Molecular packing in crystals of the rhombohedral modification is shown in Figure 5. Amino acid residues of helices α4, α8 participate in polar contacts between subunits B and B* of neighboring molecules in rhombohedral crystals connected by a twofold axis (symmetry operator y, x, −z) ( Table 2). Between subunit A of the molecule (x, y, z) and subunit A* of the molecule associated with the original molecule by the symmetry operator (−x−2/3, −x + y−1/3, −z−1/3), twice as many polar contacts exist (Table 2). These contacts are formed by amino acid residues of the same helices and by the residues of helix α1, strand β14, loop 302-306 ( Table 2). The total surface area of a molecule in the rhombohedral crystals involved in contacts is 4215.3 Å 2 , which exceeds the contact area in a tetragonal structure (3936 Å2). The number of amino acid residues involved in the 8 of 11 contacts of each molecule is 153. The contacts mainly involve the regions of the molecule with an ordered structure ( Table 2).
As discussed earlier, a trimer is located in the independent part of the tetragonal crystal cell, and a biologically active hexamer molecule is formed under the action of a second-order crystallographic axis (symmetry operator y, x, −z). In the tetragonal modification, intermolecular contacts involve molecules bounded by a twofold screw axis, which is parallel to the a axis at the 3/8 c level (symmetry operator: x−1/2, −y−1/2, −z−1/4) ( Table  2). Intermolecular contacts with four neighboring molecules form four subunits of a hexamer molecule: independent subunits В, С and subunits В*, С* of the second trimer of the molecule, obtained as a result of the action of the symmetry operator y, x, −z. Hexamer molecule and contacting subunits of other molecules in the P41212 crystal modification are shown in Figure 4. The total area of contact of one molecule with four neighboring molecules is 3936 Å 2 . In tetragonal crystals, polar intermolecular interactions involve amino acid residues belonging to the regular regions of the structure: helices α4, α5, η4, and strand β14 ( Table 2).
As mentioned, in rhombohedral crystals, a dimer is located in the independent part of the cell, and a disk-shaped hexamer molecule is formed under the action of a thirdorder crystallographic axis. In each trimer, the conformation of its subunits is the same, but the conformations of the subunits in different trimers of molecules are different. All six subunits of the hexamer are involved in intermolecular contacts in rhombohedral crystals. Intermolecular contacts are formed between trimers that have the same conformation of subunits in contacting molecules. The total number of neighboring molecules in contact with the molecule is six. Molecular packing in crystals of the rhombohedral modification is shown in Figure 5. Amino acid residues of helices α4, α8 participate in polar contacts between subunits B and B* of neighboring molecules in rhombohedral crystals connected by a twofold axis (symmetry operator y, x, −z) ( Table 2). Between subunit A of the molecule (x, y, z) and subunit A* of the molecule associated with the original molecule by the symmetry operator (−x−2/3, −x + y−1/3, −z−1/3), twice as many polar contacts exist ( Table 2). These contacts are formed by amino acid residues of the same helices and by the residues of helix α1, strand β14, loop 302-306 ( Table 2). The total surface area of a molecule in the rhombohedral crystals involved in contacts is 4215.3 Å 2 , which exceeds the contact area in a tetragonal structure (3936 Å2). The number of amino acid residues involved in the contacts of  (Figures 6a and 7a). In rhombohedral crystals, such channels are absent, and enzyme molecules in the crystal structure form a layer perpendicular to the c axis. Matthews's coefficients [24] for crystal structures belonging to space group R32 h and P41212 are 2.73 and 3.52 Å3/Da, respectively, which indicates a higher solvent content in the tetragonal structure (65.03%), as compared to the rhombohedral one (55.00%).
Comparison of the spatial structures and conformations of phosphoribosyl pyrophosphate synthetase molecules from the thermophilic bacterium Thermus thermophilus HB27, obtained at a resolution of 1.85 Å (rhombohedral crystals of the enzyme, space group R32-h) and at a resolution of 2.2 Å (tetragonal crystals, space group Р41212), revealed a number of differences. In the rhombohedral crystal modification, the conformation of independent subunits in the region of the active center differs in location and degree of ordering of the catalytic loop, which regulates the degree of openness of the active center. One of the subunits has a closed conformation of the active center, similar to the conformation of one of the subunits in tetragonal crystals.
Analysis of molecular packing in crystal lattices of both space groups and the nature of intermolecular contacts showed that the total surface area of a molecule involved in contacts is larger in rhombohedral crystals. All six subunits of a biological molecule participate in intermolecular interactions in rhombohedral crystals; in tetragonal crystals, intermolecular contacts with neighboring molecules form four subunits of a hexameric molecule. The total number of amino acid residues involved in the intermolecular interactions of each molecule with its neighbors is 153, in rhombohedral crystals, and 134 in tetragonal.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Figure S1: The fragment of the electron density map. The 2Fo-Fc is contoured at 2 r.m.s.d.   Matthews's coefficients [24] for crystal structures belonging to space group R3 2 h and P4 1 2 1 2 are 2.73 and 3.52 Å3/Da, respectively, which indicates a higher solvent content in the tetragonal structure (65.03%), as compared to the rhombohedral one (55.00%).
Comparison of the spatial structures and conformations of phosphoribosyl pyrophosphate synthetase molecules from the thermophilic bacterium Thermus thermophilus HB27, obtained at a resolution of 1.85 Å (rhombohedral crystals of the enzyme, space group R3 2 -h) and at a resolution of 2.2 Å (tetragonal crystals, space group P4 1 2 1 2), revealed a number of differences. In the rhombohedral crystal modification, the conformation of independent subunits in the region of the active center differs in location and degree of ordering of the catalytic loop, which regulates the degree of openness of the active center. One of the subunits has a closed conformation of the active center, similar to the conformation of one of the subunits in tetragonal crystals.
Analysis of molecular packing in crystal lattices of both space groups and the nature of intermolecular contacts showed that the total surface area of a molecule involved in contacts is larger in rhombohedral crystals. All six subunits of a biological molecule participate in intermolecular interactions in rhombohedral crystals; in tetragonal crystals, intermolecular contacts with neighboring molecules form four subunits of a hexameric molecule. The total number of amino acid residues involved in the intermolecular interactions of each molecule with its neighbors is 153, in rhombohedral crystals, and 134 in tetragonal.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/cryst11091128/s1, Figure S1: The fragment of the electron density map. The 2Fo-Fc is contoured at 2 r.m.s.d. Funding: This research was funded by the Russian Science Foundation (project No. 21-13-00429) and the Ministry of Science and Higher Education within the State assignment FSRC "Crystallography and Photonics" RAS in part of collection of X-ray diffraction datasets for financial support of this work. Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.