Abstract:Post-translational modification by SUMO (small ubiquitin-related modifier) is an important mechanism that regulates a wide variety of cellular functions, such as nuclear transport, gene expression, stress response, cell cycle control, oncogenesis, and response to virus infection.1-3 SUMO proteins are conjugated to target proteins by an enzymatic cascade involving SUMO activating enzyme E1 (Aos1/Uba2, also named SAE1/SAE2),4-6 SUMO conjugating enzyme E2 (Ubc9),7, 8 and E3 ligase. Trypanosoma brucei, the causative agent of African sleeping sickness, is a unicellular protozoan. The sequences of SUMO in Trypanosoma brucei (Tb-SUMO) are 37% and 33% identical with that of human SUMO-1 and yeast Smt3, respectively (Supporting Information Figure S1A). Here, the solution structure of truncated Tb-SUMO (residues 7–108) from T. brucei was determined by NMR spectroscopy. Its interaction surface with human Ubc9 was also determined by chemical shift perturbation using NMR spectroscopy. This is the first structure of SUMO determined in protist. The 3D structure of Tb-SUMO is highly conserved with that of SUMO-1 and Smt3. Tb-SUMO interacts with human Ubc9 through residues located on the β-sheet of Tb-SUMO, which is also similar to that of SUMO-1 and Smt3. All these results suggest the evolutionary conservation of SUMO and sumoylation in eukaryotes. The DNA encoding residues 7–108 of SUMO in T. brucei was cloned into the expression vector pET22b(+) (Novgen). The recombinant vector was transformed into E.coli BL21 (DE3) for expression. Recombinant Tb-SUMO has 110 amino acid residues, including a C-terminal His tag (LEHHHHHH). The Tb-SUMO proteins were purified with a Ni-chelating column (Qiagen). NMR samples contained 0.6 mM Tb-SUMO, 50 mM phosphate buffer (pH 6.25), 100 mM sodium chloride, 2 mM EDTA in either 90% H2O/10%D2O or 100% D2O.9 Human GST-Ubc9 fusion protein was produced in E. coli DH5α cells.10 The fusion protein was purified on a GSH-Sepharose 4B column (Pharmacia) according to the manufacturer's specification. The fusion protein was digested with 50 units/mL thrombin (Pharmacia) for 16 h at 20°C. All NMR data were collected at 298 K on a Bruker DMX500 spectrometer with self-shielded z-axis gradients. NMR data processing was carried out using NMRPipe and NMRDraw. [1H-15N] HSQC, HNCACB, CACB (CO)NH, HNHA, HCC(CO)NH, H(CC)CONH, 3D 15N-edited NOESY, and 13C-edited NOESY spectra were required for structure analysis. The data were analyzed with Sparky 3 software. All software was run on a Linux system. Structures were calculated using CYANA.11 Initially, only typical medium- and long-range NOEs defining secondary structures from 13C-edited NOESY spectrum were used together with NOEs obtained from 15N-edited NOESY for structure calculation. Dihedral angle restraints, hydrogen bonds, and all other NOEs were introduced in consecutive steps. The statistics on experimental constraints, coordinate precision, and stereochemical quality of 20 structures with lowest energy were analyzed with MOLMOL and PROCHECK. To investigate the human Ubc9 binding site on Tb-SUMO, 0.4 mM 15N-labeled Tb-SUMO was prepared. Unlabeled human Ubc9 was titrated into the sample containing Tb-SUMO. Six different molar ratios of Ubc9 and Tb-SUMO (1:6, 1:3, 1:2, 2:3, 5:6, and 1:1) were designed. The concentration of both proteins at the end of the titration was ∼0.25 mM. The final molar ratio of Ubc9 and Tb-SUMO is about 1:1. The HSQC spectra for the binding surface on Tb-SUMO were recorded at 298 K. The 3D structure was calculated using CYANA 2.1.11 About 97% chemical shifts in backbone and 95% chemical shifts in total were assigned. The best-fit superposition of the ensemble of the 20 lowest-energy calculated structures is shown in Figure 1A. The atomic coordinates have been deposited in the Protein Data Bank as entry 2K8H. The statistical parameters in Table I showed a high-quality NMR structure of Tb-SUMO. The NMR solution structure of Tb-SUMO consists of two α-helices (α1, 56–67; α2, 88–91) and one mixed β-sheet (β1, 33–39; β2, 43–49; β3, 73–76; β4, 80–82; β5, 96–102) [Fig. 1(B)]. The β-sheet consists of five β-strands composed of three antiparallel β-strands and one parallel β-strand. The structure represents the typical ubiquitin-like fold. Pairwise backbone (residues L32–G108) Cα RMSD between the average NMR structure of Tb-SUMO and the average NMR structure of human SUMO-1 is 2.5 Å. For the same residue range, the RMSD value between Tb-SUMO and Smt3 is 3.1 Å. The structural comparison indicates that structure of Tb-SUMO is highly similar with that of human SUMO-1 and yeast Smt313 (Supporting Information Figure S1B and 1C), suggesting the conservation of SUMO protein among eukaryotes. NMR structure of Tb-SUMO: A: Backbone superposition of 20 selected conformers (residue range, 31–108) from the final CYANA 2.111 calculation; B: ribbon representation of the minimized averaged structure with the secondary structure elements highlighted. This figure was produced with MOLMOL.12 Previous studies indicate that Ubc9 can interact with SUMO proteins through noncovalent binding.14-16 No Ubc9 homologue has been identified in T. brucei yet. Therefore, the human Ubc9 was chosen to investigate the interface between Tb-SUMO and Ubc9 in this study. The chemical shift perturbation was used to map the Ubc9-binding interface on Tb-SUMO. Residues 31, 38–43, 75–80, 92–102, and 104 have the most significant changes in chemical shifts on complex formation with Ubc9 (Figure S2). Similar to other SUMO, these residues are located in β strands 1, 3, 5, the loop between β1 and β2 and the region between β3 and β5 [Fig. 2(A); Supporting Information Figure S3A]. Among the surface residues affected by interaction between Ubc9 and Tb-SUMO, six negatively charged residues, D41, D78, E94, D95, I99, and D100, are highly conserved among SUMO family members [Fig. 2(B)]. These residues, together with E89, D96, D97, and E104, constitute a large negative potential surface (Supporting Information Figure S3B and S3C). The result suggests an important role for electrostatic attraction in the interaction between Tb-SUMO and human Ubc9. The binding surface on Ubc9 has been reported previously.15 These residues are mainly on the helix α1, β-strands 1–3, and loop between α1 and β1. These highly conserved and positive residues on this region constitute a large positive potential surface. This surface is apparently compatible with the negative potential surface of Tb-SUMO. Besides, there are some hydrophobic residues affected significantly in this region. The complementary electrostatic potential plus hydrophobic interaction contribute to the specific binding of Tb-SUMO to human Ubc9. Ubc9 binding region and conserved residues on the surface of Tb-SUMO: A: residues whose combined chemical shift changes were more than the mean value are colored in red; B: residues which are conserved on the Ubc9 binding surface are colored in dark pink. Another interesting observation in this study comes from the N-terminal region of Tb-SUMO. Residue T30, which is located on the N-terminal region near to the β1-sheet, exhibits an obvious change on Ubc9 binding. Residue T30 is a polar residue and might contribute to the hydrophobic interaction between SUMO and Ubc9. This is unique for Tb-SUMO since N-terminal region of SUMO from other organisms has never been found to interact with Ubc9. One obvious discrepancy between SUMO and ubiquitin lies in the N-terminal extension in SUMO proteins. Whether the extension has special function is not clear. Tb-SUMO has the longer N-terminal extension than human and yeast SUMO proteins. The unique feature of N-terminal extension of Tb-SUMO regarding its interaction with Ubc9 might reflect the divergence during its evolution in eukaryotes. Gene expression in eukaryotes is regulated at several different levels including transcriptional level and post-translational level. However, in T. brucei, the lack of transcriptional factors implies the deficiency of regulation at transcriptional level which makes the post-translational regulation such as sumoylation become extremely more important in regulating the gene expression. A large body of evidence has indicated that SUMO is essential for many important nuclear events, including DNA replication and repair, kinetochore-tubulin attachment, genome integrity maintenance.1, 13, 17 Independent studies in our laboratory indicated that Tb-SUMO plays an essential role in the mitosis of cell cycle in T. brucei (Liao et al., unpublished dada). Our study, based on the sequence alignment, the similarity of Tb-SUMO structure compared with other SUMO structure and the conservation of molecular mechanism of interaction between SUMO and E2 conjugating enzyme (Ubc9), confirms the existence of SUMO and imply a SUMO system in T. brucei. In conclusion, the first solution structure of Tb-SUMO from protist has been determined. The 3D structure of Tb-SUMO and its interaction with human Ubc9 exhibit highly similarity among all SUMO from eukaryotes. All these accumulated data confirm the conservation of SUMO and sumoylation from protist to vertebrate. The authors thank F. Delaglio and A. Bax for providing NMRPipe and NMRDraw software. They also thank Goddard and Kneller for providing SPARKY software. Additional Supporting Information may be found in the online version of this article. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Abstract:Ubiquitin-like protein modifiers (Ubls) are involved in diverse biological processes and regulate the activity and function of target proteins. When they are conjugated to target proteins, they are either as a regulatory post-translational modifier, or tag for target protein degradation in the proteosome by an enzymatic cascade involving an activating enzyme, a conjugating enzyme, and a ligase.1, 2 Many Ubls have been detected in eukaryotes.3 However, no protein modifier has been detected in prokaryotes. The evolutionary progress of protein modifier remains a puzzle. Urm1 (ubiquitin-related modifier 1, PF09138) belongs to Ubls family. It was found to form a thioester with Uba4, the E1-like enzyme in the Urm1 conjugation pathway.4 Urm1 plays a few important roles, including oxidant-stress response, budding in vegetative growth, invasive growth into agar in the haploid state, pseudohyphal growth, and cell elongation under starvation conditions in the diploid state in Saccharomyces cerevisiae.5, 6 Meanwhile, the solution structure of Urm1 from yeast has been determined. Structural comparison reveals that Urm1 has the similar structure with MoaD (molybdopterin synthase small subunit) and ThiS (involved in thiamin biosynthesis).7 In Escherichia coli, MoaD and ThiS are two sulfur carrier proteins that are activated in an ATP-dependent manner by sulfur transferase, MoaB, and ThiF.8, 9 The mechanism is similar to the activation of ubiquitin by the ubiquitin-activating enzyme E1.10, 11 Additionally, ubiquitin, MoaD, and ThiS all have conventional β-grasp fold. These accumulated evidences imply that MoaD and ThiS might be prokaryotic homologs of ubiquitin.12-14 Structural comparison of yeast Urm1 and MoaD and ThiS, combining phylogenetic analysis of the ubiquitin superfamily, suggests that Urm1 might act as a unique "molecular fossil" connecting evolutionary course between ATP-dependent protein conjunction in eukaryotes and ATP-dependent cofactor sulfuration.7 Trypanosoma brucei is the most ancient and evolutionarily divergent eukaryote with many unique biological features. It causes sleeping sickness in human and nagana in cattle in sub-Saharan Africa.15 Although it belongs to Ubls family, TbUrm1 (Urm1 in T. brucei) has only 11% sequence identity with ubiquitin from human. In this article, the solution structure of TbUrm1 is determined by NMR. It shares the typical β-fold structure of ubiquitin superfamily. Structural comparison reveals that TbUrm1 has the most similar structure with Urm1 from Saccharomyces cerevisiae and MoaD from Escherichia coli. These results further confirm an important role of Urm1 in the evolutionary history of ubiquitin superfamily. The gene encoding TbUrm1 was obtained from T. brucei DNA gene library by PCR and cloned into the NdeI/XhoI-cleaved plasmid PET22b(+) (Novagen), providing the C-terminal His-tagged (LEHHHHHH) protein. The recombinant vector was transformed into expression host Rosetta (DE3). The supernatant of lysed cells was collected and purified with a Ni-chelating column (Qiagen) as described before.16 NMR sample contained 0.6 mM Urm1, 50 mM phosphate (pH 6.5), 100 mM sodium chloride, 2 mM DTT in either 90% H2O/10% 2H2O or 100% 2H2O. All NMR data were collected at 298 k on a Bruker DMX600 spectrometer. A set of standard triple-resonance spectra was recorded for backbone and side chain assignments. NOE distance restraints were obtained from 3D 15N- and 13C-edited NOESY spectra acquired with a mixing time of 130 ms. After these experiments, the sample was lyophilized and redissolved in 99.96%2 H2O. A series of 15N-HSQC experiments was performed to monitor the disappearance of NH signals to obtain the hydrogen bond information. NMR data were processed with NMRPipe and analyzed with Sparky 3 software. The interproton restraints were classified into four categories: 1.8–3.0 Å, 1.8–4.0 Å, 1.8–5.0 Å, and 1.8–6.0 Å according to NOE intensities. Based on analysis of Cα, Cβ, C′, and Hα chemical shifts using the chemical shift index, the information on the φ and ψ backbone dihedral angles was obtained using TALOS program.17 Hydrogen bond restraints were obtained by assignment of slow-exchange amide protons located in regular Secondary structural elements (SSEs). All calculations were performed with CNS.18, 19 The CNS program20 was used to calculate 3D structure of Urm1 by distance restraints using the ARIA setup and protocols. Short-range NOEs and long-range NOEs were first used to determine SSEs of Urm1. φ and ψ backbone dihedral angles and hydrogen bond restraints were added in consecutive steps to constraint 3D structure of Urm1 better. The 20 structures with lowest energy were analyzed with MOLMOL.21 The Ramachandran plot was analyzed with PROCHECK.22 The NMR data used for structure calculations are summarized in Table I. The 20 structures with lowest energy from 200 calculated structures were selected for structural statistic calculation. The chemical shift assignments of TbUrm1 had been deposited in the Biological Magnetic Resonance Data Bank (accession number 15587), and the structure of TbUrm1 had been deposited in Protein Data Bank (PDB ID code 2K9X). Yeast and mouse Urm1 were used to search potential Urm1 gene of T. brucei in the GeneDB (http://www.genedb.org/). Only one homologue (Tb 927. 4. 1830) was identified. Sequence identity and similarity between TbUrm1 and yeast Urm1 are 30 and 51%, respectively. In addition, some features of TbUrm1 primary sequence are similar to that of Urm1 from other species, such as double Gly motif in their C-terminal. The structure of TbUrm1 has four β-stands and five α-helices. Strand 1 (residues 7–12) and strand 4 (residues 92–97) are parallel to each other, whereas strand 2 (residues 24–25) and strand 3 (residues 68–72) are anti-parallel. Five helices, α1 (residues 15–18), α2 (residues 32–34), α3 (residues 37–47), α4 (residues 52–57), and α5 (residues 77–81) are on the concave side of the curved β-sheet. The interaction of the inner face of α2 with the concave face of the β-sheet forms the hydrophobic core of TbUrm1, which is essential to maintain the compact structure (see Fig. 1). This type of hydrophobic core is very similar to that from other ubiquitin superfamily proteins, suggesting Urm1 shares the highly conserved structure in this family. NMR structure of TbUrm1. (A) Backbone superposition of 20 selected conformers with the lowest energy from the final CNS calculation. (B) Ribbon representation of the minimized averaged structure of TbUrm1 with the secondary structure elements highlighted. This figure was produced with MOLMOL. The minimum energy structure of TbUrm1 was submitted to the fold recognition programs DALI23 to search for its similar structure. The result showed that TbUrm1 has the similar structure with ubiquitin from human and yeast Urm1 (see Fig. 2), further indicating that the structure of Urm1 is highly conserved in its ubiquitin superfamily. Notably, the structural comparison revealed that yeast Urm1 protein (PDB ID code 2PKO) has the most similar structure with that of TbUrm1. The Cα RMSD between TbUrm1 and yeast Urm1 was 1.9 Å with corresponding Z-scores 13. Meanwhile, the Cα RMSD between TbUrm1 and MoaD (PDB ID code 1V8C-A), ThiS (PDB ID code 3CWI-A) were 2.3 Å, 2.5 Å with corresponding Z-scores 10.0, 4.9 (see Fig. 2). The ubiquitin superfamily can be classified into three clusters by different structural features including ubiquitin-related fold, ThiS-related fold, and MoaD-related fold. In the ubiquitin-related fold, the 310 helix and the core helix interact through their N-terminal residues, which are perpendicular to each other. In the ThiS-related fold, the short α-helix interacts with the core helix N terminus residues by its C-terminal region residues following the β-4. The MoaD-related fold has different feature. In the C-terminal region, two short helixes are packed together. TbUrm1 has two additional short α helixes like the MoaD-related fold and short α helix immediately following β4 similar to the ThiS-related fold. Therefore, the structure of TbUrm1 can be classified into MoaD-related fold or ThiS-related fold by structural comparison. Although distributional regions of similar residues have some difference, overall electrostatic surface distributions of TbUrm1, yeast Urm1, MoaD, and ThiS are similar (Supporting Information Fig. 1). TbUrm1 shows about 24 and 23% sequence identity to MoaD from Escherichia coli and ThiS from Escherichia coli, respectively. The sequence and structural similarity between TbUrm1 and two sulfur carrier proteins demonstrates that TbUrm1 retains structural features of the common ancestor of ubiquitin superfamily. (A) Sequence alignment of TbUrm1 with ubiquitin, Urm1 from yeast and MoaD by ClustalW24 and decorated using ESPript.25 (B) Ribbon structures of four protein modifiers. Left top: TbUrm1, right top: yeast Urm1 (2AX5), bottom left: MoaD from Thermus thermophilus (1V8C), bottom right: ubiquitin from human (1UBI). In conclusion, the NMR structure of Urm1 from T. brucei, a most ancient eukaryote, has been determined. The structure of TbUrm1 shares the highly conserved structure among its ubiquitin superfamily. Furthermore, structural comparison reveals that it has the most similar structure with MoaD from prokaryote, supporting the hypothesis that Urm1 is the unique "molecular fossil" establishing an evolutionary link between ATP-dependent protein conjunction in eukaryotes and ATP-dependent cofactor sulfuration in prokaryotes. This information will undoubtedly shed more light on the road exploring evolution of ubiquitin superfamily. The authors thank F. Delaglio and A. Bax for providing NMRPipe and NMRDraw, T. D. Goddard and D. Kneller for Sparky, A. T. Brünger for CNS, and R. Koradi and K. Wuthrich for MOLMOL. Additional Supporting Information may be found in the online version of this article. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Solution Structure of a Ubiquitin-Like Protein from Trypanosoma Brucei

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

Identification of enzymes involved in SUMOylation in Trypanosoma brucei

Solution Structure of the N-terminal Domain of DC-UbP/UBTD2 and Its Interaction with Ubiquitin.

Solution Structure of Tbtfiis2-1 Pwwp Domain from Trypanosoma Brucei

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

Solution structure of the ubiquitin-like domain of human DC-UbP from dendritic cells.

Solution Structure of Urm1 from Trypanosoma Brucei

Solution Structure of Q388A3 PDZ Domain from Trypanosoma Brucei

The Small Ubiquitin-Like Modifier (SUMO) is Essential in Cell Cycle Regulation in Trypanosoma Brucei

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

Solution Structure of TbTFIIS2-2 PWWP Domain from Trypanosoma Brucei and Its Binding to H4K17me3 and H3K32me3

Novel Bilobe Components in Trypanosoma brucei Identified Using Proximity-Dependent Biotinylation

Structural Analysis of Recombinant Human Ubiquitin-Conjugating Enzyme UbcH5c.

Crystal structure of the ubiquitin-like domain of human TBK1

High-resolution cryo-electron microscopy structure of the Trypanosoma brucei ribosome

Structure and Ubiquitin Interactions of the Conserved Zinc Finger Domain of Npl4

Crystal structure of TbAlba1 from Trypanosoma brucei

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

The Protein Neddylation Pathway in Trypanosoma Brucei

Elucidating the Mechanism Underlying UBA7-UBE2L6 Disulfide Complex Formation