Crystal structure of tRNA N2,N2-guanosine dimethyltransferase Trm1 from Pyrococcus horikoshii

Trm1 catalyzes a two-step reaction, leading to mono- and dimethylation of guanosine at position 26 in most eukaryotic and archaeal tRNAs. We report the crystal structures of Trm1 from Pyrococcus horikoshii liganded with S-adenosyl-l-methionine or S-adenosyl-l-homocysteine. The protein comprises N-terminal and C-terminal domains with class I methyltransferase and novel folds, respectively. The methyl moiety of S-adenosyl-l-methionine points toward the invariant Phe27 and Phe140 within a narrow pocket, where the target G26 might flip in. Mutagenesis of Phe27 or Phe140 to alanine abolished the enzyme activity, indicating their role in methylating G26. Structural analyses revealed that the movements of Phe140 and the loop preceding Phe27 may be involved in dissociation of the monomethylated tRNA•Trm1 complex prior to the second methylation. Moreover, the catalytic residues Asp138, Pro139, and Phe140 are in a different motif from that in DNA 6-methyladenosine methyltransferases, suggesting a different methyl transfer mechanism in the Trm1 family.     

Substrate tRNA recognition mechanism of a multisite-specific tRNA methyltransferase, Aquifex aeolicus Trm1, based on the X-ray crystal structure

Archaeal and eukaryotic tRNA (N(2),N(2)-guanine)-dimethyltransferase (Trm1) produces N(2),N(2)-dimethylguanine at position 26 in tRNA. In contrast, Trm1 from Aquifex aeolicus, a hyper-thermophilic eubacterium, modifies G27 as well as G26. Here, a gel mobility shift assay revealed that the T-arm in tRNA is the binding site of A. aeolicus Trm1. To address the multisite specificity, we performed an x-ray crystal structure study. The overall structure of A. aeolicus Trm1 is similar to that of archaeal Trm1, although there is a zinc-cysteine cluster in the C-terminal domain of A. aeolicus Trm1. The N-terminal domain is a typical catalytic domain of S-adenosyl-l-methionine-dependent methyltransferases. On the basis of the crystal structure and amino acid sequence alignment, we prepared 30 mutant Trm1 proteins. These mutant proteins clarified residues important for S-adenosyl-l-methionine binding and enabled us to propose a hypothetical reaction mechanism. Furthermore, the tRNA-binding site was also elucidated by methyl transfer assay and gel mobility shift assay. The electrostatic potential surface models of A. aeolicus and archaeal Trm1 proteins demonstrated that the distribution of positive charges differs between the two proteins. We constructed a tRNA-docking model, in which the T-arm structure was placed onto the large area of positive charge, which is the expected tRNA-binding site, of A. aeolicus Trm1. In this model, the target G26 base can be placed near the catalytic pocket; however, the nucleotide at position 27 gains closer access to the pocket. Thus, this docking model introduces a rational explanation of the multisite specificity of A. aeolicus Trm1.     

Crystal structure of archaeal tRNA(m(1)G37)methyltransferase aTrm5

Methylation of the N1 atom of guanosine at position 37 in tRNA, the position 3'-adjacent to the anticodon, generates the modified nucleoside m(1)G37. In archaea and eukaryotes, m(1)G37 synthesis is catalyzed by tRNA(m(1)G37)methyltransferase (archaeal or eukaryotic Trm5, a/eTrm5). Here we report the crystal structure of archaeal Trm5 (aTrm5) from Methanocaldococcus jannaschii (formerly known as Methanococcus jannaschii) in complex with the methyl donor analogue at 2.2 A resolution. The crystal structure revealed that the entire protein is composed of three structural domains, D1, D2, and D3. In the a/eTrm5 primary structures, D2 and D3 are highly conserved, while D1 is not conserved. The D3 structure is the Rossmann fold, which is the hallmark of the canonical class-I methyltransferases. The a/eTrm5-defining domain, D2, exhibits structural similarity to some class-I methyltransferases. In contrast, a DALI search with the D1 structure yielded no structural homologues. In the crystal structure, D3 contacts both D1 and D2. The residues involved in the D1:D3 interactions are not conserved, while those participating in the D2:D3 interactions are well conserved. D1 and D2 do not contact each other, and the linker between them is disordered. aTrm5 fragments corresponding to the D1 and D2-D3 regions were prepared in a soluble form. The NMR analysis of the D1 fragment revealed that D1 is well folded by itself, and it did not interact with either the D2-D3 fragment or the tRNA. The NMR analysis of the D2-D3 fragment revealed that it is well folded, independently of D1, and that it interacts with tRNA. Furthermore, the D2-D3 fragment was as active as the full-length enzyme for tRNA methylation. The positive charges on the surface of D2-D3 may be involved in tRNA binding. Therefore, these findings suggest that the interaction between D1 and D3 is not persistent, and that the D2-D3 region plays the major role in tRNA methylation.     

Crystal complexes of a predicted S-adenosylmethionine-dependent methyltransferase reveal a typical AdoMet binding domain and a substrate recognition domain

S-adenosyl-L-methionine-dependent methyltransferases (MTs) are abundant, and highly conserved across phylogeny. These enzymes use the cofactor AdoMet to methylate a wide variety of molecular targets, thereby modulating important cellular and metabolic activities. Thermotoga maritima protein 0872 (TM0872) belongs to a large sequence family of predicted MTs, ranging phylogenetically from relatively simple bacteria to humans. The genes for many of the bacterial homologs are located within operons involved in cell wall synthesis and cell division. Despite preliminary biochemical studies in E. coli and B. subtilis, the substrate specificity of this group of more than 150 proteins is unknown. As part of the Midwest Center for Structural Genomics initiative (www.mcsg.anl.gov), we have determined the structure of TM0872 in complexes with AdoMet and with S-adenosyl-L-homocysteine (AdoHcy). As predicted, TM0872 has a typical MT domain, and binds endogenous AdoMet, or co-crystallized AdoHcy, in a manner consistent with other known MT structures. In addition, TM0872 has a second domain that is novel among MTs in both its location in the sequence and its structure. The second domain likely acts in substrate recognition and binding, and there is a potential substrate-binding cleft spanning the two domains. This long and narrow cleft is lined with positively charged residues which are located opposite the S(+)-CH(3) bond, suggesting that a negatively charged molecule might be targeted for catalysis. However, AdoMet and AdoHcy are both buried, and access to the methyl group would presumably require structural rearrangement. These TM0872 crystal structures offer the first structural glimpses at this phylogenetically conserved sequence family.     

Structure of RsrI methyltransferase, a member of the N6-adenine beta class of DNA methyltransferases

DNA methylation is important in cellular, developmental and disease processes, as well as in bacterial restriction-modification systems. Methylation of DNA at the amino groups of cytosine and adenine is a common mode of protection against restriction endonucleases afforded by the bacterial methyltransferases. The first structure of an N:6-adenine methyltransferase belonging to the beta class of bacterial methyltransferases is described here. The structure of M. RSR:I from Rhodobacter sphaeroides, which methylates the second adenine of the GAATTC sequence, was determined to 1.75 A resolution using X-ray crystallography. Like other methyltransferases, the enzyme contains the methylase fold and has well-defined substrate binding pockets. The catalytic core most closely resembles the PVU:II methyltransferase, a cytosine amino methyltransferase of the same beta group. The larger nucleotide binding pocket observed in M. RSR:I is expected because it methylates adenine. However, the most striking difference between the RSR:I methyltransferase and the other bacterial enzymes is the structure of the putative DNA target recognition domain, which is formed in part by two helices on an extended arm of the protein on the face of the enzyme opposite the active site. This observation suggests that a dramatic conformational change or oligomerization may take place during DNA binding and methylation.     

The DNMT1 target recognition domain resides in the N terminus

DNA-cytosine-5-methyltransferase 1 (DNMT1) is the enzyme believed to be responsible for maintaining the epigenetic information encoded by DNA methylation patterns. The target recognition domain of DNMT1, the domain responsible for recognizing hemimethylated CGs, is unknown. However, based on homology with bacterial cytosine DNA methyltransferases it has been postulated that the entire catalytic domain, including the target recognition domain, is localized to 500 amino acids at the C terminus of the protein. The N-terminal domain has been postulated to have a regulatory role, and it has been suggested that the mammalian DNMT1 is a fusion of a prokaryotic methyltransferase and a mammalian DNA-binding protein. Using a combination of in vitro translation of different DNMT1 deletion mutant peptides and a solid-state hemimethylated substrate, we show that the target recognition domain of DNMT1 resides in the N terminus (amino acids 122-417) in proximity to the proliferating cell nuclear antigen binding site. Hemimethylated CGs were not recognized specifically by the postulated catalytic domain. We have previously shown that the hemimethylated substrates utilized here act as DNMT1 antagonists and inhibit DNA replication. Our results now indicate that the DNMT1-PCNA interaction can be disrupted by substrate binding to the DNMT1 N terminus. These results point toward new directions in our understanding of the structure-function of DNMT1.     

Isolation and identification by sequence homology of a putative cytosine methyltransferase from Arabidopsis thaliana.

A plant cytosine methyltransferase cDNA was isolated using degenerate oligonucleotides, based on homology between prokaryote and mouse methyltransferases, and PCR to amplify a short fragment of a methyltransferase gene. A fragment of the predicted size was amplified from genomic DNA from Arabidopsis thaliana. Overlapping cDNA clones, some with homology to the PCR amplified fragment, were identified and sequenced. The assembled nucleic acid sequence is 4720 bp and encodes a protein of 1534 amino acids which has significant homology to prokaryote and mammalian cytosine methyltransferases. Like mammalian methylases, this enzyme has a C terminal methyltransferase domain linked to a second larger domain. The Arabidopsis methylase has eight of the ten conserved sequence motifs found in prokaryote cytosine-5 methyltransferases and shows 50% homology to the murine enzyme in the methyltransferase domain. The amino terminal domain is only 24% homologous to the murine enzyme and lacks the zinc binding region that has been found in methyltransferases from both mouse and man. In contrast to mouse where a single methyltransferase gene has been identified, a small multigene family with homology to the region amplified in PCR has been identified in Arabidopsis thaliana.     

The M.AluI DNA-(cytosine C5)-methyltransferase has an unusually large, partially dispensable, variable region.

The DNA methyltransferase of the AluI restriction-modification system, from Arthrobacter luteus, converts cytosine to 5-methylcytosine in the sequence AGCT. The gene for this methyltransferase, aluIM, was cloned into Escherichia coli and sequenced. A 525-codon open reading frame was found, consistent with deletion evidence, and the deduced amino acid sequence revealed all ten conserved regions common to 5-methylcytosine methyltransferases. The aluIM sequence predicts a protein of M(r) 59.0k, in agreement with the observed M(r), making M.AluI the largest known methyltransferase from a type II restriction-modification system. M.AluI also contains the largest known variable region of any monospecific DNA methyltransferase, larger than that of most multispecific methyltransferases. In other DNA methyltransferases the variable region has been implicated as the sequence-specific target recognition domain. An in-frame deletion that removes a third of this putative target-recognition region leaves the Alu I methyltransferase still fully active.     

Crystal structure of hypothetical protein TTHB192 from Thermus thermophilus HB8 reveals a new protein family with an RNA recognition motif-like domain

We have determined the crystal structure of hypothetical protein TTHB192 from Thermus thermophilus HB8 at 1.9 A resolution. This protein is a member of the Escherichia coli ygcH sequence family, which contains approximately 15 sequence homologs of bacterial origin. These homologs have a high isoelectric point. The crystal structure reveals that TTHB192 consists of two independently folded domains, and that each domain exhibits a ferredoxin-like fold with a four-stranded antiparallel beta-sheet packed on one side by alpha-helices. These two tandem domains face each other to generate a beta-sheet platform. TTHB192 displays overall structural similarity to Sex-lethal protein and poly(A)-binding protein fragments. These proteins have RNA binding activity which is supported by a beta-sheet platform formed by two tandem repeats of an RNA recognition motif domain with signature sequence motifs on the beta-sheet surface. Although TTHB192 does not have the same signature sequence motif as the RNA recognition motif domain, the presence of an evolutionarily conserved basic patch on the beta-sheet platform could be functionally relevant for nucleic acid-binding. This report shows that TTHB192 and its sequence homologs adopt an RNA recognition motif-like domain and provides the first testable functional hypothesis for this protein family.     

N2-methylation of guanosine at position 10 in tRNA is catalyzed by a THUMP domain-containing, S-adenosylmethionine-dependent methyltransferase, conserved in Archaea and Eukaryota

In sequenced genomes, genes belonging to the cluster of orthologous group COG1041 are exclusively, and almost ubiquitously, found in Eukaryota and Archaea but never in Bacteria. The corresponding gene products exhibit a characteristic Rossmann fold, S-adenosylmethionine-dependent methyltransferase domain in the C terminus and a predicted RNA-binding THUMP (thiouridine synthases, RNA methyltransferases, and pseudouridine synthases) domain in the N terminus. Recombinant PAB1283 protein from the archaeon Pyrococcus abyssi GE5, a member of COG1041, was purified and shown to behave as a monomeric 39-kDa entity. This protein (EC 2.1.1.32), now renamed (Pab)Trm-G10, which is extremely thermostable, forms a 1:1 complex with tRNA and catalyzes the adenosylmethionine-dependent methylation of the exocyclic amino group (N(2)) of guanosine located at position 10. Depending on the experimental conditions used, as well as the tRNA substrate tested, the enzymatic reaction leads to the formation of either N(2)-monomethyl (m(2)G) or N(2)-dimethylguanosine (m(2)(2)G). Interestingly, (Pab)Trm-G10 exhibits different domain organization and different catalytic site architecture from another, earlier characterized, tRNA-dimethyltransferase from Pyrococcus furiosus ((Pfu)Trm-G26, also known as (Pfu)Trm1, a member of COG1867) that catalyzes an identical two-step dimethylation of guanosine but at position 26 in tRNAs and is also conserved among all sequenced Eukaryota and Archaea. The co-occurrence of these two guanosine dimethyltransferases in both Archaea and Eukaryota but not in Bacteria is a hallmark of distinct tRNAs maturation strategies between these domains of life.     

Mining folded proteomes in the era of accurate structure prediction

Protein structure fundamentally underpins the function and processes of numerous biological systems. Fold recognition algorithms offer a sensitive and robust tool to detect structural, and thereby functional, similarities between distantly related homologs. In the era of accurate structure prediction owing to advances in machine learning techniques and a wealth of experimentally determined structures, previously curated sequence databases have become a rich source of biological information. Here, we use bioinformatic fold recognition algorithms to scan the entire AlphaFold structure database to identify novel protein family members, infer function and group predicted protein structures. As an example of the utility of this approach, we identify novel, previously unknown members of various pore-forming protein families, including MACPFs, GSDMs and aerolysin-like proteins.     

Exploring ORFan Domains in Giant Viruses: Structure of Mimivirus Sulfhydryl Oxidase R596

The mimivirus genome contains many genes that lack homologs in the sequence database and are thus known as ORFans. In addition, mimivirus genes that encode proteins belonging to known fold families are in some cases fused to domain-sized segments that cannot be classified. One such ORFan region is present in the mimivirus enzyme R596, a member of the Erv family of sulfhydryl oxidases. We determined the structure of a variant of full-length R596 and observed that the carboxy-terminal region of R596 assumes a folded, compact domain, demonstrating that these ORFan segments can be stable structural units. Moreover, the R596 ORFan domain fold is novel, hinting at the potential wealth of protein structural innovation yet to be discovered in large double-stranded DNA viruses. In the context of the R596 dimer, the ORFan domain contributes to formation of a broad cleft enriched with exposed aromatic groups and basic side chains, which may function in binding target proteins or localization of the enzyme within the virus factory or virions. Finally, we find evidence for an intermolecular dithiol/disulfide relay within the mimivirus R596 dimer, the first such extended, intersubunit redox-active site identified in a viral sulfhydryl oxidase.     

The amino acid sequence of the eukaryotic DNA [N6-adenine]methyltransferase, M.CviBIII, has regions of similarity with the prokaryotic isoschizomer M.TaqI and other DNA [N6-adenine] methyltransferases

The sequences of the genes coding for M.CviBIII (from virus NC-1A which infects a eukaryotic alga) [Narva et al., Nucleic Acids Res. 15 (1987) 9807-9823] and M.TaqI (from the bacterium Thermus aquaticus) [Slatko et al., Nucleic Acids Res. 15 (1987) 9781-9796] have been determined recently. Both enzymes methylate adenine in the sequence TCGA. We have compared the predicted amino acid sequences of these two methyltransferases (MTases), with each other and with ten other N6 A-MTases and find regions of similarity. M.CviBIII and M.TaqI were most closely related followed by M.PaeR7, whose recognition sequence (CTCGAG) contains the M.TaqI/M.CviBIII recognition sequence TCGA, and M.PstI, whose recognition sequence is CTGCAG. All of the N6-MTases contain the sequence Asp/Asn-Pro-Pro-Tyr (B-P-P-Y) referred to by Hattman et al. [J. Bacteriol. 164 (1985) 932-937] as region IV. The predicted secondary structure of this region forms a finger-like structure ('beta finger') containing a beta-pleated sheet (...XXXB), two beta-turns (P-P) followed by another beta-pleated sheet [Y/FXXX...].     

The 2'-O-methyltransferase responsible for modification of yeast tRNA at position 4

The methylation of the ribose 2'-OH of RNA occurs widely in nature and in all stable RNAs and occurs at five positions in yeast tRNA. 2'-O-methylation of tRNA at position 4 is interesting because it occurs in the acceptor stem (which is normally undermodified), it is the only 2'-O-methylation that occurs in the middle of a duplex region in tRNA, the modification is conserved in eukaryotes, and the features of the tRNA necessary for substrate recognition are poorly defined. We show here that Saccharomyces cerevisiae ORF YOL125w (TRM13) is necessary and sufficient for 2'-O-methylation at position 4 of yeast tRNA. Biochemical analysis of the S. cerevisiae proteome shows that Trm13 copurifies with 2'-O-methylation activity, using tRNAGlyGCC as a substrate, and extracts made from a trm13-Delta strain have undetectable levels of this activity. Trm13 is necessary for activity in vivo because tRNAs isolated from a trm13-Delta strain lack the corresponding 2'-O-methylated residue for each of the three known tRNAs with this modification. Trm13 is sufficient for 2'-O-methylation at position 4 in vitro since yeast Trm13 protein purified after expression in Escherichia coli has the same activity as that produced in yeast. Trm13 protein binds substrates tRNAHis and tRNAGlyGCC with KD values of 85+/-8 and 100+/-14 nM, respectively, and has a KM for tRNAHis of 10 nM, but binds nonsubstrate tRNAs very poorly (KD>1 microM). Trm13 is conserved in eukaryotes, but there is no sequence similarity between Trm13 and other known methyltransferases.     

Prediction of a novel RNA 2'-O-ribose methyltransferase subfamily encoded by the Escherichia coli YgdE open reading frame and its orthologs

The amino acid sequence of the RNA 2'-O-ribose methyltranserase RrmJ was used as a probe for detecting putative homologs through iterative searches of genomic databases. We found a previously unannotated YgdE open reading frame (ORF) in the genome sequences of Escherichia coli and other gamma-Proteobacteria, which shares key features with RrmJ, despite the mutual sequence similarity of these proteins is relatively low. The predicted structural compatibility and the conservation of all functionally important residues between RrmJ and YgdE strongly suggests that the newly identified methyltranserase also modifies 2'-OH groups of ribose. The N-terminal region of YgdE, which has no counterpart in RrmJ, is predicted to form an independent domain, possibly involved in target recognition.     

Structure of pvu II DNA-(cytosine N4) methyltransferase, an example of domain permutation and protein fold assignment.

We have determined the structure of Pvu II methyltransferase (M. Pvu II) complexed with S -adenosyl-L-methionine (AdoMet) by multiwavelength anomalous diffraction, using a crystal of the selenomethionine-substituted protein. M. Pvu II catalyzes transfer of the methyl group from AdoMet to the exocyclic amino (N4) nitrogen of the central cytosine in its recognition sequence 5'-CAGCTG-3'. The protein is dominated by an open alpha/beta-sheet structure with a prominent V-shaped cleft: AdoMet and catalytic amino acids are located at the bottom of this cleft. The size and the basic nature of the cleft are consistent with duplex DNA binding. The target (methylatable) cytosine, if flipped out of the double helical DNA as seen for DNA methyltransferases that generate 5-methylcytosine, would fit into the concave active site next to the AdoMet. This M. Pvu IIalpha/beta-sheet structure is very similar to those of M. Hha I (a cytosine C5 methyltransferase) and M. Taq I (an adenine N6 methyltransferase), consistent with a model predicting that DNA methyltransferases share a common structural fold while having the major functional regions permuted into three distinct linear orders. The main feature of the common fold is a seven-stranded beta-sheet (6 7 5 4 1 2 3) formed by five parallel beta-strands and an antiparallel beta-hairpin. The beta-sheet is flanked by six parallel alpha-helices, three on each side. The AdoMet binding site is located at the C-terminal ends of strands beta1 and beta2 and the active site is at the C-terminal ends of strands beta4 and beta5 and the N-terminal end of strand beta7. The AdoMet-protein interactions are almost identical among M. Pvu II, M. Hha I and M. Taq I, as well as in an RNA methyltransferase and at least one small molecule methyltransferase. The structural similarity among the active sites of M. Pvu II, M. Taq I and M. Hha I reveals that catalytic amino acids essential for cytosine N4 and adenine N6 methylation coincide spatially with those for cytosine C5 methylation, suggesting a mechanism for amino methylation.     

Characterization of the Arabidopsis thaliana Arath;CDC25 dual-specificity tyrosine phosphatase

CDC25 enzymes are dual-specificity phosphatases involved in the regulation of the cell cycle. No CDC25 enzymes have been described in higher plant organisms. We report here the characterization of an Arabidopsis thaliana CDC25 enzyme, constituted by a sole catalytic domain and devoid of the N-terminal regulatory region found in the human CDC25. We describe the recombinant expression in Escherichia coli of the Arath;CDC25 and its purification for activity assay and structure determination by NMR. The recombinant enzyme has a tyrosine phosphatase activity towards an artificial substrate, a NMR characterization equally concludes to its correct folding. The secondary structure of the protein was predicted on the basis of the assigned chemical shift of (1)H, (15)N, and (13)C backbone atoms of the protein. The presence of a metal ion in the C-terminus of this new protein points to a zinc finger, and sequence homology indicates that this new structural element might be conserved in related plant homologs.