Abstract:Ubiquitins and ubiquitin-like modifier proteins regulate numerous cellular processes via covalent conjugation to the target substrates.1-3 To date, five ubiquitin-like pathways have been identified in the yeast Saccharomyces cerevisiae.4-6 A protein of 99 residues encoded by yeast open reading frame YIL008W was characterized as one of the ubiquitin-related modifiers (Urm1). Conjugation of Urm1 to its target proteins depends on the E1-like activating enzyme Uba4.7 A recent article on the solution structure and phylogenetic analyses of Urm1 indicated it might be a common ancestor of the entire ubiquitin superfamily.8 To date, little is known about the molecular function of urmylation pathway, but is required for normal growth, particularly at high temperature.7 Urm1 was also found to play important roles in budding and haploid invasive growth,9 posttranslational modification/proteolysis of the elongator subunit Tot1/YLR384C10 and oxidative stress response via modifying the alkyl hydroperoxide reductase Ahp1.11 Very recently, urmylation was also reported to participate in regulating the expression of genes involved in sensing and controlling amino acids levels.12 Oligomerization have been reported to be a common feature of ubiquitins and related modifiers.13 Ubiquitin forms polyubiquitin conjugation via isopeptide bonds between the active Lys of one monomer and the C-terminal Gly of another.3 SUMO-2/-3 forms oligomers in vitro via a specific consensus SUMOylation site, VKXE,14 which was further proved by the crystal structures of the human SUMO-2 protein in a possible assembly form of trimer.15 Nevertheless, the oligomerization of Urm1 has never been reported previously. Here, we reported the crystal structure of Urm1 from the yeast S. cerevisiae at 1.44 Å resolution, dimerized via hydrophobic and polar interactions, mainly contributed by the C-terminal residues. The stable dimer was also found to be the form of majority in solution. The dimerization of Urm1 might provide an autoprotection to the highly active C-terminal Gly residue before attacking its substrates. The URM1/YIL008W open reading frame was amplified by polymerase chain reaction (PCR), using S. cerevisiae S288c genomic DNA as the template, and cloned into a pET29a-derived plasmid (Novagen) without introducing any additional residues. The recombinant protein was overexpressed in E. coli Rosetta (DE3) strain at 37°C for 5 h. Pellets were collected by centrifugation at 8,000g for 10 min at 4°C and resuspended in 50 mM NaCl, 20 mM Tris-HCl, pH 7.0, followed by sonication. After centrifugation at 25,000g for 30 min at 4°C, the supernatant was pooled and applied to an ion-exchange Q-FF column (Amersham Biosciences) equilibrated with 20 mM Tris-HCl, pH 7.0. Elution of 40 mL volume was carried out by a linear gradient of NaCl from 0 to 1.0M at a flow rate of 1 mL/min. Urm1 was eluted at 250 mM NaCl and then pooled and loaded to a Superdex 200 16/60 column (Amersham Biosciences), equilibrated with 200 mM NaCl, 20 mM Tris-HCl, pH 7.0. Fractions containing Urm1 were pooled, desalted, and concentrated to 10 mg/mL in a final buffer of 50 mM NaCl, 20 mM Tris-HCl, pH 7.0, used for crystallization and analytical size exclusion chromatography. Protein concentration was determined by its absorbance at 280 nm with the theoretical molar extinction coefficient of 6,990 M−1 cm−1 (http://www.expasy.org). Urm1 was also cloned into pET15b plasmid (Novagen) with additional residues of MGSSHHHHHHSSGLVPRGH at the N-terminus. The recombinant protein Urm1-6xHis-tag was overexpressed by using Rosetta (DE3) strain of E. coli at 16°C for 20 h. Pellets were collected by centrifugation at 8000g for 10 min at 4°C and resuspended in 200 mM NaCl, 20 mM Tris-HCl, pH 7.0, followed by sonication. After centrifugation at 25,000g for 30 min 4°C, the supernatant was collected and applied to a Ni2+-NTA column (Amersham Biosciences). The target protein was eluted with 300 mM imidazole in 200 mM NaCl, 20 mM Tris-HCl, pH 7.0, and loaded to a Superdex 75 16/60 column, equilibrated with 200 mM NaCl, 20 mM Tris-HCl, pH 7.0. Fractions containing the recombinant Urm1-6xHis-tag were pooled, desalted, and concentrated to 10 mg/mL in a final buffer of 50 mM NaCl, 20 mM Tris-HCl. Protein concentration was determined by its absorbance at 280 nm with a theoretical molar extinction coefficient of 6990 M−1cm−1 (http://www.expasy.org). The crystals were grown at 289 K using the sitting-drop vapor-diffusion technique, with initial conditions formed by mixing equal volumes (1 μL) of the protein sample with mother liquor [0.2M MgCl2, 0.1M Tris, pH8.5, 30% polyethylene glycol (PEG) 4000]. Typically, crystals of about 0.30 × 0.30 × 0.05 mm3 appeared in a week. The crystal was transferred to the cryoprotectant of the reservoir solution supplemented with 25% glycerol and flash cooled with liquid nitrogen. The X-ray diffraction data were collected at 100 K in a liquid nitrogen stream using a Rigaku MM007 X-ray generator (λ = 1.5418 Å) with a MarResearch 345 image-plate detector at School of Life Sciences, University of Science and Technology of China (USTC, Hefei, China). Data were processed using the program AUTOMAR 1.2.16 The crystal structure of Urm1 was solved by the molecular replacement method with the program Phaser17 using the mean solution structure of Urm1 (PDB code 2AX5) without the flexible residues at both termini as the search model. According to the 2Fo-Fc and Fo-Fc maps, the initial model was fitted and rebuilt with the program O18 and refined with program CNS19 using reflections between 30 and 3.0 Å with 5% of the data set reserved for Rfree. When the Rfactor and Rfree were decreased to 43.8% and 44.7%, respectively, the resultant model was applied to restrained refinement with REFMAC5.18, 20 As refinement progressed, the resolution was extended to 1.44 Å and anisotropic temperature factors were applied in refinement. Finally, 248 water molecules, two magnesium ions, and one ethylene glycol (EDO) were added according to the Fo-Fc map. The side chains of eight residues were modeled in two discrete alternative conformations. The stereochemical quality of the final structure was verified using the program PROCHECK V 3.4.4.21 The data collection parameters, final refinement statistics, and the overall quality parameters of the structure were listed in Table I. All figures were prepared with the program PyMOL.22 An aliquot of 20 μL protein mixture of molecular weight (MW) maker was loaded to HPLC size exclusion column Bio-Sil SEC 400-5 (Bio-Rad), equilibrated with 50 mM NaCl, 50 mM sodium phosphate buffer, pH 6.5, followed by elution at a flow rate of 1 mL/min, and absorbance measured at 280 nm. The protein markers are bovine thyroglobulin (670 kDa), bovine γ-globulin (158 kDa), chicken ovalbumin (44 kDa), horse myoglobin (17 kDa), yeast recombinant Grx1 (13.4 kDa), and Grx2 (14.1 kDa). Using the peak elution volumes, the standard curve equation was determined as Log MW= −0.6369Ve + 10.948 with the R2 value of 0.922. Urm1 was diluted to 10 μM in the same sodium phosphate buffer and then loaded to the column. The molecular weight of Urm1 was calculated from the equation with its peak elution volume. Urm1 (11.0 kDa) and Urm1-6xHis-tag (13.1 kDa) were mixed with a molar ratio of 10:1. The mixture was added to 10 mL buffer of 200 mM NaCl, 20 mM Tris-HCl, pH 7.0 with 8M urea, followed by slow dialysis at 4°C overnight to a final buffer of 200 mM NaCl, 20 mM Tris-HCl, pH 7.0 with 1.5M urea. The sample was then collected, desalted to remove the residual urea and concentrated in a final buffer of 50 mM NaCl, 20 mM Tris-HCl, pH 7.0. Then, 50 μL protein sample was added with 20 μL Ni2+-NTA gel and incubated at room temperature for 5 min. Supernatant was removed after centrifugation at 10,000g for 1 min. The sediment was resuspended with 60 μL washing buffer of 50 mM NaCl, 20 mM Tris-HCl, pH 7.0, followed by centrifugation as above to remove supernatant containing unbound Urm1. After the gel was washed three times, Urm1-6xHis-tag was then eluted with 60 μL washing buffer containing 300 mM imidazole, followed by centrifugation. All supernatants were applied to SDS-PAGE to check the protein. The overall fold of Urm1 crystal structure resembles its solution structure (Fig. 1), with a five-strand β-sheet and four helices on the concave side of the curved β-sheet. Nevertheless, the main chain superposition between the crystal structure and the mean solution structure gave a root-mean-square deviation (r.m.s.d.) of 1.192 Å, mainly due to the region from Pro49 to Ile61 and the C terminus. The former region consists of helix α3′ flanked by a loop on each side in the solution structure [Fig. 1(B)], compared with a region of helices α2 and η2 followed by a small antiparallel β-sheet composed of β3 and β4 in the crystal structure [Fig. 1(A)]. In the solution structure strand β4′ is absent from the cartoon representation but identified by nuclear overhauser effect (NOE) connectivity,8 while in the crystal structure this β-strand was numbered as strand β6. Besides, helices η3 and η4 were identified in the crystal structure, which is corresponding to helix α4′ in the solution structure. Cartoon representation of the overall fold of Urm1 in (A) the crystal structure and (B) the mean solution structure, colored and labeled according to the secondary structures. Strand β4′ in the solution structure was identified by NOE connectivity. In Urm1 solution structure, the C terminus from Thr95 to Gly99 extends as a flexible tail protruding to the solvent, probably due to the C-terminal tag LEHHHHHH introduced during cloning.8 In contrast, the C-terminal tail twists and covers the region between strands β4 and β5 in the crystal structure [Fig. 1(A)]. Although there is only one Urm1 molecule in an asymmetric unit of the crystal structure, crystal packing produces two possible dimeric assemblies. Calculated with the program AreaIMol in CCP4,23 the contact area of two dimeric assemblies is 297 and 692 Å2, respectively. The larger one accounts for 13.7% of the total molecular surface area [Fig. 2(A)], which is assumed to be the majority of Urm1 dimer. LIGPLOT26 diagram represents that the interface is formed by residues of Glu7, Gly10, Asp13, Arg20, Ile66, Leu68, Asp71, Asp73, Glu75, and Thr93 to Gly99, and a group of water molecules [Fig. 2(B)]. Hydrophobic interactions are found contributed by Gly10, Ile66, Leu68, Asp71, Asp 73, Thr93, Leu96, His97, and Gly99. Beside Ile66, Leu68, and Gly99, other residues are also involved in polar interactions, some of which forming hydrogen bonds mediated by water molecules. These residues are highly conserved among the Urm1 homologues, and C-terminal residues are almost identical [Fig. 2(C)]. The structure of Urm1 crystal obtained from the crystallization condition of 30% PEG8000, 0.1M sodium cacodylate, pH 6.5, and 0.2M NaAc also showed the same crystal packing and conformation of the C-terminal tail (data not shown). This indicated that Urm1 is of great possibility to form a dimer at the pH range from 6.5 to 8.5, which is close to the physiological pH value. A: Dimeric Urm1 in the crystal structure. Monomers are colored in yellow and lightblue, respectively. The N terminus and C terminus are labeled, accordingly. B: LIGPLOT diagram of the interface of dimeric Urm1. Residues at the first surface are colored in blue, residues at the second surface in brown; hydrogen bonds are shown as green dashed lines, labeled in Angstrom; residues in hydrophobic contacts are presented by semicircle with radiating spokes. C: Multiple sequence alignment of Urm1 from Saccharomyces cerevisiae to its homologs from other species. The alignment was performed using MultAlin and ESPript (http://prodes.toulouse.inra.fr/multalin).24, 25 The secondary structural elements are identified from Urm1 (PDB code 2QJL) using ESPript and displayed at the top of the alignment. The α-helices, η-helices, β-sheets, and strict β-turns are denoted α, η, β, and TT, correspondingly. Residues of >70% consensus level are in blue rectangle. Completely conserved residues are indicated by white lettering on a red background. The residues involved in intermolecular interactions are indicated by asterisks at the bottom of the alignment. In an attempt to verify if Urm1 exists as a dimer in solution as well, we determined its molecular weight (MW) with 10 μM protein in the solution of 50 mM NaCl, 50 mM sodium phosphate buffer, pH 6.5. Calculated from the standard curve [Fig. 3(A)], the MW of Urm1 is 21.9 kD, which is very close to twice of the theoretical MW of the monomer (11.0 kDa). In addition, no elution peak of the monomer was found, which suggested the majority of Urm1 exists as dimer in the solution. Furthermore, renaturation between Urm1 (11.0 kDa, without any additional tag) and Urm1-6xHis-tag (13.1 kDa) was applied to check the existence of dimeric Urm1 in solution. During renaturation, heterodimeric Urm1 was speculated to be formed between Urm1 and Urm1-6xHis-tag. Urm1-6xHis-tag was then chelated with Ni2+-NTA gel and the redundant Urm1 was removed completely by centrifugation and washing. As shown in Figure 3(B), the imidazole elution sample in lane 7 contained two bands corresponding to Urm1 and Urm1-6xHis-tag, respectively, which validated the existence of dimeric Urm1 in solution. A: The molecular weight of Urm1 in the buffer of 50 mM NaCl, 50 mM Tirs, pH 7.0, determined from the standard curve of analytical size exclusion chromatography. The protein markers are bovine thyroglobulin (670 kDa), bovine γ-globulin (158 kDa), chicken ovalbumin (44 kDa), horse myoglobin (17 kDa), recombinant yeast Grx1 (13.4 kDa), and Grx2 (14.1 kDa). B: SDS-PAGE analysis of renaturation samples applied to Ni2+-NTA gel. Lane 1, renaturation sample; lane 2, marker; lane 3, unbound Urm1 to Ni2+-NTA gel; lane 4, first gel washing; lane 5, second gel washing; lane 6, third gel washing; lane 7, elution with buffer containing 300 mM imidazole. Numbers on the left indicate the positions of MW standards (in kDa). Taken together, Urm1 forms a dimer either in solution or in the crystal packing. The dimeric form will block the mobile C-terminal Gly-Gly, which could provide an autoprotection of this active site. As we know, the C-terminus of Urm1 should be recognized and activated by E1-like Uba4 before conjugation, as strongly indicated from the bacterial counterpart MoaD-MoeB complex.27 Thus, the dipeptide Gly-Gly of Urm1 should be exposed to let it accessible to its protein partners, not only for activation but also for conjugation. But what triggers the dissociation of Urm1 dimer remains an open question. We also thank Mr. Zhiqiang Zhu at USTC for advice on structure solution.

Abstract:The TM1506 gene of Thermotoga maritima encodes a protein with a molecular weight of 15.7 kDa (residues 1–139) and a calculated isoelectric point (pI) of 8.9. TM1506 was selected for structure determination to extend the structural coverage of the T. maritima proteome (coverage data: http://ffas.burnham.org/ffas-cgi/cgi/tm_cov.pl), which is currently one of the highest for a single organism.1 TM1506 is a member of an uncharacterized protein family (PF08973, DUF1893) but, in contrast to other members of this family, has an additional Lys/Arg-rich N-terminal α-helix and a basic pI. The function of TM1506 is unknown, although its genomic neighborhood indicates possible functional associations with S12 and other ribosomal proteins. Here, we report the crystal structure of TM1506, which was determined using the semiautomated, high-throughput pipeline of the Joint Center for Structural Genomics (JCSG),2 as part of the National Institute of General Medical Sciences' Protein Structure Initiative (PSI). The gene encoding TM1506 (TIGR: TM1506, Swiss-Prot: Q9X1J4) was amplified by polymerase chain reaction (PCR) from genomic DNA using PfuTurbo DNA polymerase (Stratagene) and primers corresponding to the predicted 5′ and 3′ ends. The PCR product was cloned into plasmid pMH1, which encodes an expression and purification tag (MGSDKIHHHHHH) at the amino terminus of the full-length protein. The cloning junctions were confirmed by DNA sequencing. Protein expression was performed in a selenomethionine-containing medium using the Escherichia coli methionine auxotrophic strain DL41. At the end of fermentation, lysozyme was added to the culture to a final concentration of 250 μg/mL, and the cells were harvested. After one freeze/thaw cycle, the cells were sonicated in Lysis Buffer [50 mM Tris pH 7.9, 50 mM NaCl, 10 mM imidazole, and 0.25 mM Tris(2-carboxyethyl)phosphine hydrochloride (TCEP)], and the lysate was centrifuged at 3400 × g for 60 min. The supernatant was applied to nickel-chelating resin (GE Healthcare) pre-equilibrated with Lysis Buffer, the resin was washed with Wash Buffer [50 mM potassium phosphate pH 7.8, 300 mM NaCl, 40 mM imidazole, 10% (v/v) glycerol, and 0.25 mM TCEP], and the protein was eluted with Elution Buffer [20 mM Tris pH 7.9, 300 mM imidazole, 10% (v/v) glycerol, 0.25 mM TCEP]. The eluate was buffer exchanged with Buffer Q [20 mM Tris pH 7.9, 50 mM NaCl, 5% (v/v) glycerol, 0.25 mM TCEP] and applied to a RESOURCE Q column (GE Healthcare) pre-equilibrated with the same buffer. The flow-through fraction, which contained TM1506, was buffer exchanged with Crystallization Buffer [20 mM Tris pH 7.9, 150 mM NaCl, and 0.25 mM TCEP] and concentrated for crystallization assays to 15 mg/mL by centrifugal ultrafiltration (Millipore). TM1506 was crystallized using the nanodroplet vapor diffusion method3 with standard JCSG crystallization protocols.2 The crystallization reagent that produced the crystal used for structure solution contained 40% (v/v) polyethylene glycol 600, 0.1 M imidazole pH 8.0, and 0.2 M Zn(OAc)2 (final pH 5.8). No additional cryoprotectant was added to the crystal. Initial screening for diffraction was carried out using the Stanford Automated Mounting system (SAM)4 and an X-ray microsource5 installed in a Stanford Synchrotron Radiation Laboratory beamline (SSRL, Menlo Park, CA). The crystal was indexed in hexagonal space group P6222 (Table I).6, 7 Molecular weight and oligomeric state of TM1506 were determined using a 1 cm × 30 cm Superdex 200 column (GE Healthcare) in combination with static light scattering (Wyatt Technology). The mobile phase consisted of 20 mM Tris pH 8.0 and 150 mM NaCl. Mutation of Asp56 to Ala and cloning of the D56A mutant was achieved using the polymerase incomplete primer extension (PIPE) method.8 The appropriate substitution and the absence of unwanted mutations were confirmed by sequencing the insert from both directions. Native, mutant (D56A), and selenomethionine-labeled proteins, expressed with an N-terminal tag (MGSDKI HHHHHH), were analyzed by LC-MS by coupling a reversed-phase Protein Captrap column (Michrom Bioresources) to a Q-Tof 2 mass spectrometer (Micromass) operating in V mode. The proteins were desalted on the column with 0.2% (v/v) formic acid and 4% (v/v) acetonitrile in water and then eluted with 0.2% (v/v) formic acid and 65% acetonitrile (v/v) in water. Deconvolution of the spectra using the MaxEnt1 algorithm was performed using MassLynx Software (Waters). Native and selenomethionine-labeled proteins were reduced with 10 mM dithiothreitol in 100 mM ammonium bicarbonate for 1 h at room temperature and then thioalkylated by the addition of 20 mM iodoacetamide for 30 min in the dark. The proteins were then digested overnight at 37°C by adding trypsin at a 1:40 mass ratio (trypsin/protein). Native protein was further digested for 8 h at 37°C with endoproteinase Glu-C at a mass ratio of 1:40 (protease/protein). The digested proteins were analyzed by LC-MS/MS by interfacing a reversed-phase C18-packed capillary column to a Q-Tof Ultima mass spectrometer (Micromass) operating in data-dependent MS/MS switching mode. Multi-wavelength anomalous diffraction (MAD) data were collected at the Advanced Light Source (ALS, Berkeley, CA) on beamline 8.2.1 at wavelengths corresponding to the inflection (λ1), low energy remote (λ2), and peak (λ3) of a selenium MAD experiment. The data sets were collected at 100 K using an ADSC Q210 CCD detector. The MAD data were integrated and reduced using Mosflm9 and then scaled with the program SCALA from the CCP4 suite.6 Data statistics are summarized in Table I. Phasing was performed with SHELXD10 and autoSHARP,11 and automated model building was performed with ARP/wARP12 and RESOLVE.13 Model completion and refinement were performed with Xfit14 and REFMAC5.15 Refinement statistics are summarized in Table I. Analysis of the stereochemical quality of the model was accomplished using AutoDepInputTool,16 MolProbity,17 SFcheck 4.0,6 and WHATIF 5.0.18 Protein quaternary structure analysis was performed using the PQS (Protein Quaternary Structure) server,19 the PISA (Protein Interfaces, Surfaces and Assemblies) server,20 and PITA (Protein InTerfaces and Assemblies) software.21 Figure 1(B) was adapted from an analysis using PDBsum.22 Figures 1(A) and 2(B,D) were prepared with PyMOL (DeLano Scientific) and ProtSkin (http://www.mcgnmr.ca/ProtSkin/). Sequence conservation in Figure 2(C) was calculated using rate4site.23 Atomic coordinates and experimental structure factors for TM1506 from T. maritima at 2.70 Å resolution have been deposited in the PDB and are accessible under the code 1vk9. Crystal structure of TM1506 from Thermotoga maritima. (A) Stereo ribbon diagram of the TM1506 monomer color-coded from N-terminus (blue) to C-terminus (red). Helices H1–H6 and β-strands β1–β5 are indicated. (B) Diagram showing the secondary structural elements of TM1506 superimposed on its primary sequence. The α-helices, β-strands, and β-turns are indicated. The β-sheet is indicated by a red A and the β-hairpin is depicted as a red loop. Dashed lines indicate regions that are not included in the protein structure. ADP-ribosylation reaction, model of TM1506 with conserved residues and ADP-ribose shown on its surface, and sequence alignment. (A) Proposed reaction from NAD+ to ADP-ribosylated Asp56. (B) Residues conserved in TM1506 are shown on its surface, where the intensity of red color indicates the degree of conservation. The ADP-ribose (shown as a stick model colored by atom type: carbon, green; nitrogen, blue; oxygen, red; phosphorus, orange) was manually fitted into the extra electron density in the active site; metals and waters are shown as spheres (zinc, blue; magnesium, green; water, red). (C) Multiple sequence alignment of TM1506 and its homologs. The intensity of red color indicates the degree of sequence conservation. The conservation scores (0, not conserved; 4, highly conserved) were derived from an analysis using rate4site software using all available homologs of TM1506 and not just those shown in the alignment. The positions of the most conserved residues in TM1506 are indicated above the alignment. Homologs of TM1506 used in the alignment include: ZP_01189901, conserved hypothetical protein from Halothermothrix orenii H 168; YP_001244876, hypothetical protein Tpet_1286 from Thermotoga petrophila RKU-1; ZP_01354436, putative TonB-dependent outer membrane receptor from Clostridium phytofermentans ISDg; YP_001319085, domain of unknown function DUF1893 from Alkaliphilus metalliredigenes QYMF; YP_101571, putative TonB-dependent outer membrane receptor from Bacteroides fragilis YCH46; NP_809874, hypothetical protein BT0961 from Bacteroides thetaiotaomicron VPI-5482. (D) The site of ADP-ribosylation. A plausible model of ADP-ribose that fits the density well is shown as a stick model (carbon, grey; nitrogen, blue; oxygen, red; phosphorus, orange); metals and waters are shown as spheres (zinc, blue; magnesium, green; water, red). The experimental map is shown after solvent flattening. Side-chain atoms of Asp56 and Cys113 are shown as sticks (carbon, yellow; oxygen, red). The crystal structure of TM1506 was determined to 2.70 Å resolution using the MAD method (see Fig. 1). Data collection, model, and refinement statistics are summarized in Table I. The final model includes one monomer (residues 1–136 plus the entire N-terminal expression and purification tag except for the initial Met residue, whose absence was confirmed by LC-MS), six zinc ions, an unknown ligand (UNL), and nine water molecules in the asymmetric unit. It is assumed that the Met from the tag was removed by an endogenous methionine aminopeptidase. No electron density was observed for residues 137–139 or the side chains of Glu2, Lys3, Arg7, Lys37, Lys107, Lys110, Glu126, and Glu127. The Matthews' coefficient (Vm)24 for TM1506 is 4.9 Å3/Da, and the estimated solvent content is 74.9%. The Ramachandran plot produced by MolProbity25 shows that 99.3% of the residues are in favored regions with no outliers. TM1506 is composed of five β-strands (β1–β5) and six α-helices (Lys3–Lys15, Lys37–Arg45, Lys61–Met70, Lys82–Glu90, Pro114–Leu120, and Pro125–Leu133; Fig. 1). The total β-sheet and α-helical content is 21.1% and 38.8%, respectively. The TM1506 monomer comprises a single domain and its secondary structural elements are arranged in a three-layer core (α/β/α) with a mixed β-sheet of five strands (order: 21345), where strand 1 is antiparallel to the rest. As per the SCOP classification scheme,26 TM1506 belongs to the cytidine deaminase-like fold and is classified as a member of the cytidine deaminase-like SCOP superfamily (sunid: 53927). A DALI27 structural similarity search of TM1506 found several hits to proteins with a cytidine deaminase-like fold (PDB codes 1uaq, 1g8m, 2b3j, 2b3z, 2g84, 1teo, and 1tiy) with significant Z-scores (8.5, 8.5, 7.3, 7.3, 6.9, 6.1, and 5.2, respectively), but residues important for zinc binding and catalysis in those enzymes are not conserved in TM1506. In addition, the position of the C-terminal helices (H5 and H6) in TM1506 occludes the typical dimerization interface that, in many cases, is functionally important for proteins with a cytidine deaminase-like fold.28, 29 The structurally similar (DALI Z-score = 5.8) T4 bacteriophage 2′-deoxycytidylate deaminase (PDB code 1vq2) was crystallized with an inhibitor located in a deep grove that is not present in TM1506. Comparison of TM1506 with tRNA adenosine deaminase TadA (PDB code 2b3j, crystallized with tRNA) indicates that TM1506 may be a tRNA-binding protein, consistent with the finding that several nucleotide-binding proteins are targeted for ADP-ribosylation.30 Analytical size-exclusion chromatography coupled with static light scattering suggests a monomer as the likely oligomeric state of TM1506, while crystallographic packing analysis using the PQS server predicts a homotetramer. This prediction is likely artifactual because of the presence of the expression/purification tag, which makes up a significant part of the tetramerization interface. Analysis of the TM1506 structure with the tag removed (as well as the associated zinc atoms) using the PISA server and PITA software suggests a monomer. The actual observed mass of the native protein was 17,576.6 Da, as determined by LC-MS, and corresponded to a delta mass of +540.7 Da based on the assumption that the N-terminal methionine from the tag was removed by an endogenous methionine aminopeptidase. Unmodified protein lacking this methionine has a calculated average mass of 17,035.9 Da. Analysis of selenomethionine-labeled protein revealed a mass of 17,811.6 Da, which also corresponded to a delta mass of +540.7 Da. Protein devoid of the methionine from the tag and having all five methionine residues substituted with selenomethionine has a theoretical average mass of 17,270.9 Da. LC-MS data suggest that TM1506 is covalently modified with ADP-ribose, since a protein with this modification would have an additional mass of 541 Da. The ADP-ribosylated peptides 45-RFDNLEGSLVIDK-57 and 46-FDNLEGSLVIDK-57 were identified by LC-MS/MS analysis of tryptic digests of native and selenomethionine-labeled proteins. Several tandem MS spectra revealed intense signal at m/z 136 ([adenine+H]1+), m/z 250 ([adenosine+H-H2O]1+), m/z 348 ([AMP+H]1+), and m/z 428 ([ADP+H]1+), indicative of fragmentation of the covalently bound ADP-ribose. The ADP-ribosylated peptide 51-GSLVIDK-57 was detected by analysis of the tryptic/endoproteinase Glu-C dual digest of native protein. The tandem MS spectrum suggested that Asp56 was modified, given that the C-terminal backbone fragments (y ions) y2–y5 were found ADP-ribosylated or as fragments that had undergone the neutral loss of AMP from the ADP-ribose moiety. Furthermore, as determined by LC-MS, a D56A mutant of TM1506 lacked the additional mass of 541 Da, supporting our hypothesis that Asp56 is the site of modification. An analysis of the TM1506 structure using the CASTp server31 revealed a deep cavity of 366 Å3 formed by Ser19, Leu20, Asp33, Ser34, Gly35, Leu36, Pro38, Val39, Asp56, Lys57, Met58, Gly60, Ala62, Val104, Cys113, Phe115, and Glu116. Location of electron density in this cavity corresponds to the presence of a ligand covalently bound to the conserved Asp56 (see Fig. 2). As indicated earlier, mass spectrometry revealed that the protein is ADP-ribosylated and suggested that the invariant Asp56 is the site of modification. As expected, the electron density closely approximates an ADP-ribose. However, the adenosine moiety, which is exposed to solvent, is disordered and appears to be partly displaced by a Zn2+ ion that is coordinated to Glu99. The Zn2+ ion was confirmed by anomalous difference Fourier. Because of the limited resolution, this ambiguity cannot be fully resolved. As a result, the ligand was assigned as a UNL in the coordinates deposited in the PDB. Figure 2(A) illustrates the proposed reaction32 leading from NAD+ to ADP-ribosylated Asp56 and a possible conformation of ADP-ribose (manually docked into the extra density). Figure 2(D) illustrates one model where ADP-ribose fits reasonably well into the density. Additionally, two metals are likely to be located in the pocket to provide stabilization for the phosphoryl groups. Based on anomalous difference maps, a heavy metal (possibly Zn2+) is likely located near Cys113. A second, lighter metal, such as Mg2+, is probably adjacent to Ser34. The information presented here, in combination with further biochemical and biophysical studies, should (1) yield valuable insights into the functional role of this protein, (2) help to determine whether NAD+ or its derivatives can bind in the conserved cleft, (3) reveal the function of ADP-ribosylation (ADP-ribosylation reaction has been proposed as a potential new drug target33), (4) confirm tRNA binding, and (5) elucidate if a high proportion of basic residues is important for the function of TM1506. The JCSG has developed The Open Protein Structure Annotation Network (TOPSAN), a wiki-based community project to collect, share, and distribute information about protein structures determined at PSI centers. TOPSAN offers a combination of automatically generated, as well as comprehensive, expert-curated annotations, provided by JCSG personnel and members of the research community. Additional information about TM1506 is available at https://www.topsan.org/1VK9. Portions of this research were performed at the Advanced Light Source (ALS) and Stanford Synchrotron Radiation Laboratory (SSRL). The SSRL is a national user facility operated by Stanford University on behalf of the United States Department of Energy, Office of Basic Energy Sciences. The SSRL Structural Molecular Biology Program is supported by the Department of Energy, Office of Biological and Environmental Research, and by the National Institutes of Health (National Center for Research Resources, Biomedical Technology Program, and the National Institute of General Medical Sciences).

Solution Structure of Urm1 from Trypanosoma Brucei

Solution structure of SUMO from Trypanosoma brucei and its interaction with Ubc9.

Crystal Structure Of Human Spermidine/Spermine N'-Acetyltransferase (Hssat): The First Structure Of A New Sequence Family Of Transferase Homologous Superfamily

Solution Structure of TbUfm1 from Trypanosoma Brucei and Its Binding to TbUba5

Solution Structure of Urm1 and Its Implications for the Origin of Protein Modifiers.

Crystal Structure Of The Dimeric Urm1 From The Yeast Saccharomyces Cerevisiae

The dual role of ubiquitin-like protein Urm1 as a protein modifier and sulfur carrier

Crystal structure of the copper homeostasis protein (CutCm) from Shigella flexneri at 1.7 Å resolution: The first structure of a new sequence family of TIM barrels

Solution Structure of a Ubiquitin-Like Protein from Trypanosoma Brucei

Monomeric Trna (M(7)G46) Methyltransferase From Escherichia Coli Presents A Novel Structure At The Function-Essential Insertion

1 H, 13 C and 15 N Resonance Assignment of Urm1 from Trypanosoma Brucei

Crystal Structure of the Hypothetical Protein Ttha1013 from Thermus Thermophilus Hb8

Complex Structure Of Ospi And Ubc13: The Molecular Basis Of Ubc13 Deamidation And Convergence Of Bacterial And Host E2 Recognition

The emerging roles of ubiquitin-like protein Urm1 in eukaryotes

Solution Structure of Human SUMO-3 C47S and Its Binding Surface for Ubc9,

Crystal Structures of the Apo and Gdp-Bound Forms of A Cupin-Like Protein Bbduf985 from Branchiostoma Belcheri Tsingtauense

Crystal Structure of B-Subtilis Yjcg Characterizing the Yjcg-Like Group of 2h Phosphoesterase Superfamily

Solution structure of the ubiquitin-associated domain of human BMSC-UbP and its complex with ubiquitin.

Molecular Basis of Threonine ADP-ribosylation of Ubiquitin by Bacterial ARTs.

Molecular Mechanism for the Substrate Recognition of USP7

Crystal Structure of an ADP‐ribosylated Protein with a Cytidine Deaminase‐like Fold, but Unknown Function (TM1506), from Thermotoga Maritima at 2.70 Å Resolution