Probabilistic Approach to Predicting Substrate Specificity of Methyltransferases

We present a general probabilistic framework for predicting the substrate specificity of enzymes. We designed this approach to be easily applicable to different organisms and enzymes. Therefore, our predictive models do not rely on species-specific properties and use mostly sequence-derived data. Maximum Likelihood optimization is used to fine-tune model parameters and the Akaike Information Criterion is employed to overcome the issue of correlated variables. As a proof-of-principle, we apply our approach to predicting general substrate specificity of yeast methyltransferases (MTases). As input, we use several physico-chemical and biological properties of MTases: structural fold, isoelectric point, expression pattern and cellular localization. Our method accurately predicts whether a yeast MTase methylates a protein, RNA or another molecule. Among our experimentally tested predictions, 89% were confirmed, including the surprising prediction that YOR021C is the first known MTase with a SPOUT fold that methylates a substrate other than RNA (protein). Our approach not only allows for highly accurate prediction of functional specificity of MTases, but also provides insight into general rules governing MTase substrate specificity.     

Alanine-scanning mutagenesis of the predicted rRNA-binding domain of ErmC' redefines the substrate-binding site and suggests a model for protein-RNA interactions

The Erm family of adenine-N(6) methyltransferases (MTases) is responsible for the development of resistance to macrolide-lincosamide-streptogramin B antibiotics through the methylation of 23S ribosomal RNA. Hence, these proteins are important potential drug targets. Despite the availability of the NMR and crystal structures of two members of the family (ErmAM and ErmC', respectively) and extensive studies on the RNA substrate, the substrate-binding site and the amino acids involved in RNA recognition by the Erm MTases remain unknown. It has been proposed that the small C-terminal domain functions as a target-binding module, but this prediction has not been tested experimentally. We have undertaken structure-based mutational analysis of 13 charged or polar residues located on the predicted rRNA-binding surface of ErmC' with the aim to identify the area of protein-RNA interactions. The results of in vivo and in vitro analyses of mutant protein suggest that the key RNA-binding residues are located not in the small domain, but in the large catalytic domain, facing the cleft between the two domains. Based on the mutagenesis data, a preliminary three-dimensional model of ErmC' complexed with the minimal substrate was constructed. The identification of the RNA-binding site of ErmC' may be useful for structure-based design of novel drugs that do not necessarily bind to the cofactor-binding site common to many S-adenosyl-L- methionine-dependent MTases, but specifically block the substrate-binding site of MTases from the Erm family.     

Sequence analysis and structure prediction of aminoglycoside-resistance 16S rRNA:m7G methyltransferases

Methylation of G1405 within bacterial 16S ribosomal RNA results in high-level resistance to specific combinations of aminoglycoside antibiotics. Only a few closely related methyltransferases (MTases), which carry out the respective modification (here dubbed "Agr", for aminoglycoside resistance), are known. It is not clear, whether they are related to "typical" S-adenosylmethionine (AdoMet)-dependent MTases or not. Demydchuk et al., 1998 proposed that the cofactor-binding region is localized at the C-terminus of Agr MTases, which implies an interesting case of sequence permutation. Since the Agr MTases lack significant sequence similarity to other proteins, we tested that hypothesis using more sensitive sequence/structure threading approach. Structure prediction confirmed the presence of a putative AdoMet-binding site in these proteins, albeit at a distinct location, resembling that of "typical", non-permuted MTases. Additionally, a small alpha-helical domain dissimilar to other proteins in the database was identified in the N-terminal region of Agr MTases. Comparison of a three-dimensional model of the Agr family member with a recently solved structure of reovirus mRNA capping MTase suggests that the mechanism of guanine-N7 methylation in rRNA and mRNA may be different.     

Identification and preliminary characterization of an O6-methylguanine DNA repair methyltransferase in the yeast Saccharomyces cerevisiae

Saccharomyces cerevisiae contains a DNA repair methyltransferase (MTase) that repairs O6-methylguanine. Methyl groups are irreversibly transferred from O6-methylguanine in DNA to a 25-kilodalton protein in S. cerevisiae cell extracts, and methyl transfer is accompanied by the formation of S-methylcysteine. The yeast MTase is expressed at approximately 150 molecules/cell in exponentially growing yeast cultures but is not detectable in stationary phase cells. Unlike mammalian and bacterial MTases, the yeast MTase is very temperature-sensitive, having a half-life of about 4 min at 37 degrees C, which may explain why others have failed to detect it. Like other DNA repair MTases, the S. cerevisiae MTase repairs O6-methylguanine more efficiently in double-stranded DNA than in single-stranded DNA. Synthesis of the yeast DNA MTase is apparently not inducible by sublethal exposures to alkylating agent, but rather MTase activity is depleted in cells exposed to low doses of alkylating agent. Judging from its molecular weight and substrate specificity, the yeast DNA MTase is more closely related to mammalian MTases than to Escherichia coli MTases.     

Protein methylation in mitochondria

Many proteins are modified by posttranslational methylation, introduced by a number of methyltransferases (MTases). Protein methylation plays important roles in modulating protein function and thus in optimizing and regulating cellular and physiological processes. Research has mainly focused on nuclear and cytosolic protein methylation, but it has been known for many years that also mitochondrial proteins are methylated. During the last decade, significant progress has been made on identifying the MTases responsible for mitochondrial protein methylation and addressing its functional significance. In particular, several novel human MTases have been uncovered that methylate lysine, arginine, histidine, and glutamine residues in various mitochondrial substrates. Several of these substrates are key components of the bioenergetics machinery, e.g., respiratory Complex I, citrate synthase, and the ATP synthase. In the present review, we report the status of the field of mitochondrial protein methylation, with a particular emphasis on recently discovered human MTases. We also discuss evolutionary aspects and functional significance of mitochondrial protein methylation and present an outlook for this emergent research field.     

In silico identification, structure prediction and phylogenetic analysis of the 2'-O-ribose (cap 1) methyltransferase domain in the large structural protein of ssRNA negative-strand viruses

The Escherichia coli RrmJ gene product has recently been shown to be the 23S rRNA:U2552 specific 2'-O-ribose methyltransferase (MTase) (RrmJ). Its structure has been solved and refined to 1.5 A resolution, demonstrating conservation of the three-dimensional fold and key catalytic side chains with the vaccinia virus VP39 protein, which functions as an mRNA 5'm(7)G-cap-N-specific 2'-O-ribose MTase. Using the amino acid sequence of RrmJ as an initial probe in an iterative search of sequence databases, we identified a homologous domain in the sequence of the L protein of non-segmented, negative-sense, single-stranded RNA viruses. The plausibility of the prediction was confirmed by homology modeling and checking whether important residues at substrate/ligand-binding sites were conserved. The predicted structural compatibility and the conservation of the active site between the novel putative MTase domain and genuine 2'-O-ribose MTases, together with the available results of biochemical studies, strongly suggest that this domain is a 5'm(7)G-cap-N-specific 2'-O-ribose MTase (i.e. the cap 1 MTase). Evolutionary relationships between these proteins are also discussed.     

Characterization of DNA methyltransferase specificities using single-molecule, real-time DNA sequencing

DNA methylation is the most common form of DNA modification in prokaryotic and eukaryotic genomes. We have applied the method of single-molecule, real-time (SMRT®) DNA sequencing that is capable of direct detection of modified bases at single-nucleotide resolution to characterize the specificity of several bacterial DNA methyltransferases (MTases). In addition to previously described SMRT sequencing of N6-methyladenine and 5-methylcytosine, we show that N4-methylcytosine also has a specific kinetic signature and is therefore identifiable using this approach. We demonstrate for all three prokaryotic methylation types that SMRT sequencing confirms the identity and position of the methylated base in cases where the MTase specificity was previously established by other methods. We then applied the method to determine the sequence context and methylated base identity for three MTases with unknown specificities. In addition, we also find evidence of unanticipated MTase promiscuity with some enzymes apparently also modifying sequences that are related, but not identical, to the cognate site.     

Characterization of Protein Methyltransferases Rkm1, Rkm4, Efm4, Efm7, Set5 and Hmt1 Reveals Extensive Post-Translational Modification

Protein methylation is one of the major post-translational modifications (PTMs) in the cell. In Saccharomyces cerevisiae, over 20 protein methyltransferases (MTases) and their respective substrates have been identified. However, the way in which these MTases are modified and potentially subject to regulation remains poorly understood. Here, we investigated six overexpressed S. cerevisiae protein MTases (Rkm1, Rkm4, Efm4, Efm7, Set5 and Hmt1) to identify PTMs of potential functional relevance. We identified 48 PTM sites across the six MTases, including phosphorylation, acetylation and methylation. Forty-two sites are novel. We contextualized the PTM sites in structural models of the MTases and revealed that many fell in catalytic pockets or enzyme–substrate interfaces. These may regulate MTase activity. Finally, we compared PTMs on Hmt1 with those on its human homologs PRMT1, PRMT3, CARM1, PRMT6 and PRMT8. This revealed that several PTMs are conserved from yeast to human, whereas others are only found in Hmt1. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD006767.     

Structure and mechanism of the rebeccamycin sugar 4'-O-methyltransferase RebM

The 2.65-angstroms crystal structure of the rebeccamycin 4'-O-methyltransferase RebM in complex with S-adenosyl-l-homocysteine revealed RebM to adopt a typical S-adenosylmethionine-binding fold of small molecule O-methyltransferases (O-MTases) and display a weak dimerization domain unique to MTases. Using this structure as a basis, the RebM substrate binding model implicated a predominance of nonspecific hydrophobic interactions consistent with the reported ability of RebM to methylate a wide range of indolocarbazole surrogates. This model also illuminated the three putative RebM catalytic residues (His140/141 and Asp166) subsequently found to be highly conserved among sequence-related natural product O-MTases from GC-rich bacteria. Interrogation of these residues via site-directed mutagenesis in RebM demonstrated His140 and Asp166 to be most important for catalysis. This study reveals RebM to be a member of the general acid/base-dependent O-MTases and, as the first crystal structure for a sugar O-MTase, may also present a template toward the future engineering of natural product MTases for combinatorial applications.     

Structural and functional insights into the molecular mechanism of rRNA m6A methyltransferase RlmJ

RlmJ catalyzes the m⁶A2030 methylation of 23S rRNA during ribosome biogenesis in Escherichia coli. Here, we present crystal structures of RlmJ in apo form, in complex with the cofactor S-adenosyl-methionine and in complex with S-adenosyl-homocysteine plus the substrate analogue adenosine monophosphate (AMP). RlmJ displays a variant of the Rossmann-like methyltransferase (MTase) fold with an inserted helical subdomain. Binding of cofactor and substrate induces a large shift of the N-terminal motif X tail to make it cover the cofactor binding site and trigger active-site changes in motifs IV and VIII. Adenosine monophosphate binds in a partly accommodated state with the target N6 atom 7 Å away from the sulphur of AdoHcy. The active site of RlmJ with motif IV sequence ₁₆₄DPPY₁₆₇ is more similar to DNA m⁶A MTases than to RNA m⁶₂A MTases, and structural comparison suggests that RlmJ binds its substrate base similarly to DNA MTases T4Dam and M.TaqI. RlmJ methylates in vitro transcribed 23S rRNA, as well as a minimal substrate corresponding to helix 72, demonstrating independence of previous modifications and tertiary interactions in the RNA substrate. RlmJ displays specificity for adenosine, and mutagenesis experiments demonstrate the critical roles of residues Y4, H6, K18 and D164 in methyl transfer.     

Predicting protein-protein interactions from protein domains using a set cover approach

One goal of contemporary proteome research is the elucidation of cellular protein interactions. Based on currently available protein-protein interaction and domain data, we introduce a novel method, maximum specificity set cover (MSSC), for the prediction of protein-protein interactions. In our approach, we map the relationship between interactions of proteins and their corresponding domain architectures to a generalized weighted set cover problem. The application of a greedy algorithm provides sets of domain interactions which explain the presence of protein interactions to the largest degree of specificity. Utilizing domain and protein interaction data of S. cerevisiae, MSSC enables prediction of previously unknown protein interactions, links that are well supported by a high tendency of coexpression and functional homogeneity of the corresponding proteins. Focusing on concrete examples, we show that MSSC reliably predicts protein interactions in well-studied molecular systems, such as the 26S proteasome and RNA polymerase II of S. cerevisiae. We also show that the quality of the predictions is comparable to the maximum likelihood estimation while MSSC is faster. This new algorithm and all data sets used are accessible through a Web portal at http://ppi-cse.nd.edu     

Role of conserved active site residues in catalysis by phospholipase B1 from Cryptococcus neoformans

Phospholipase B1 (PLB1), secreted by the pathogenic yeast Cryptococcus neoformans, has an established role in virulence. Although the mechanism of its phospholipase B, lysophospholipase, and lysophospholipase transacylase activities is unknown, it possesses lipase, subtilisin protease aspartate, and phospholipase motifs containing putative catalytic residues S146, D392, and R108, respectively, conserved in fungal PLBs and essential for human cytosolic phospholipase A2 (cPLA2) catalysis. To determine the role of these residues in PLB1 catalysis, each was substituted with alanine, and the mutant cDNAs were expressed in Saccharomyces cerevisiae. The mutant PLB1s were deficient in all three enzymatic activities. As the active site structure of PLB1 is unknown, a homology model was developed, based on the X-ray structure of the cPLA2 catalytic domain. This shows that the two proteins share a closely related fold, with the three catalytic residues located in identical positions as part of a single active site, with S146 and D392 forming a catalytic dyad. The model suggests that PLB1 lacks the "lid" region which occludes the cPLA2 active site and provides a mechanism of interfacial activation. In silico substrate docking studies with cPLA2 reveal the binding mode of the lipid headgroup, confirming the catalytic dyad mechanism for the cleavage of the sn-2 ester bond within one of two separate binding tracts for the lipid acyl chains. Residues specific for binding arachidonic and palmitic acids, preferred substrates for cPLA2 and PLB1, respectively, are identified. These results provide an explanation for differences in substrate specificity between lipases sharing the cPLA2 catalytic domain fold and for the differential effect of inhibitors on PLB1 enzymatic activities.     

Structure, function and mechanism of exocyclic DNA methyltransferases

DNA MTases (methyltransferases) catalyse the transfer of methyl groups to DNA from AdoMet (S-adenosyl-L-methionine) producing AdoHcy (S-adenosyl-L-homocysteine) and methylated DNA. The C5 and N4 positions of cytosine and N6 position of adenine are the target sites for methylation. All three methylation patterns are found in prokaryotes, whereas cytosine at the C5 position is the only methylation reaction that is known to occur in eukaryotes. In general, MTases are two-domain proteins comprising one large and one small domain with the DNA-binding cleft located at the domain interface. The striking feature of all the structurally characterized DNA MTases is that they share a common core structure referred to as an 'AdoMet-dependent MTase fold'. DNA methylation has been reported to be essential for bacterial virulence, and it has been suggested that DNA adenine MTases (Dams) could be potential targets for both vaccines and antimicrobials. Drugs that block Dam could slow down bacterial growth and therefore drug-design initiatives could result in a whole new generation of antibiotics. The transfer of larger chemical entities in a MTase-catalysed reaction has been reported and this represents an interesting challenge for bio-organic chemists. In general, amino MTases could therefore be used as delivery systems for fluorescent or other reporter groups on to DNA. This is one of the potential applications of DNA MTases towards developing non-radioactive DNA probes and these could have interesting applications in molecular biology. Being nucleotide-sequence-specific, DNA MTases provide excellent model systems for studies on protein-DNA interactions. The focus of this review is on the chemistry, enzymology and structural aspects of exocyclic amino MTases.     

High plasticity of multispecific DNA methyltransferases in the region carrying DNA target recognizing enzyme modules.

Multispecific cytosine C5 DNA methyltransferases (MTases) methylate more than one specific DNA target. This is due to the presence of several target recognizing domains (TRDs) in these enzymes. Such TRDs form part of a variable centre in the MTase primary sequence, which separates conserved enzyme core sequences responsible for general steps in the methylation reaction. By deleting, rearranging and exchanging several TRDs of multispecific MTases, we demonstrate their modular character; they mediate target recognition independent of a particular TRD or core sequence context. We show also that multispecific MTases can accommodate inert material of non-MTase origin within their variable region without losing their activity. The remarkable plasticity with respect to the material that can be integrated into this region suggests that the enzyme core sequences preceding or following it form separable functional domains. In spite of the documented flexibility multispecific MTases could not be endowed with novel specificities by integration of putative TRDs of monospecific MTases, pointing to differences between multi- and monospecific MTases in the way their core and TRD sequences interact.     

Crystal structure of tRNA m1G9 methyltransferase Trm10: insight into the catalytic mechanism and recognition of tRNA substrate

Transfer RNA (tRNA) methylation is necessary for the proper biological function of tRNA. The N¹ methylation of guanine at Position 9 (m¹G9) of tRNA, which is widely identified in eukaryotes and archaea, was found to be catalyzed by the Trm10 family of methyltransferases (MTases). Here, we report the first crystal structures of the tRNA MTase spTrm10 from Schizosaccharomyces pombe in the presence and absence of its methyl donor product S-adenosyl-homocysteine (SAH) and its ortholog scTrm10 from Saccharomyces cerevisiae in complex with SAH. Our crystal structures indicated that the MTase domain (the catalytic domain) of the Trm10 family displays a typical SpoU-TrmD (SPOUT) fold. Furthermore, small angle X-ray scattering analysis reveals that Trm10 behaves as a monomer in solution, whereas other members of the SPOUT superfamily all function as homodimers. We also performed tRNA MTase assays and isothermal titration calorimetry experiments to investigate the catalytic mechanism of Trm10 in vitro. In combination with mutational analysis and electrophoretic mobility shift assays, our results provide insights into the substrate tRNA recognition mechanism of Trm10 family MTases.     

Identification and Characterization of a Novel Evolutionarily Conserved Lysine-specific Methyltransferase Targeting Eukaryotic Translation Elongation Factor 2 (eEF2)

The components of the cellular protein translation machinery, such as ribosomal proteins and translation factors, are subject to numerous post-translational modifications. In particular, this group of proteins is frequently methylated. However, for the majority of these methylations, the responsible methyltransferases (MTases) remain unknown. The human FAM86A (family with sequence similarity 86) protein belongs to a recently identified family of protein MTases, and we here show that FAM86A catalyzes the trimethylation of eukaryotic elongation factor 2 (eEF2) on Lys-525. Moreover, we demonstrate that the Saccharomyces cerevisiae MTase Yjr129c, which displays sequence homology to FAM86A, is a functional FAM86A orthologue, modifying the corresponding residue (Lys-509) in yeast eEF2, both in vitro and in vivo. Finally, Yjr129c-deficient yeast cells displayed phenotypes related to eEF2 function (i.e. increased frameshifting during protein translation and hypersensitivity toward the eEF2-specific drug sordarin). In summary, the present study establishes the function of the previously uncharacterized MTases FAM86A and Yjr129c, demonstrating that these enzymes introduce a functionally important lysine methylation in eEF2. Based on the previous naming of similar enzymes, we have redubbed FAM86A and Yjr129c as eEF2-KMT and Efm3, respectively.     

Sequence motifs specific for cytosine methyltransferases

Using a new alignment method, the sequences of 13 m5C methyltransferases (MTases) have been examined. Five extremely well-conserved blocks of sequence have been detected and have been used as fixed points for the alignment of the 13 sequences. Following this initial alignment, five further blocks of similarity have been identified to give a total of ten recognizable blocks of sequence homology that are all arranged in a common order. The structures of these MTases consist of a variable-length N-terminal arm followed by eight well-conserved blocks each separated by small variable-length regions. A large variable-length segment of 90 to 270 amino acids (aa) then follows. After this are two blocks, and a variable-length C-terminal segment completes the sequence. Within the final alignment, 20 aa in the protein sequences, and 86 nucleotides in the nucleotide sequences are invariant. The strongest conservation is found in proximity to a suspected functional site that contains the dipeptide proline-cysteine. Consensus patterns can be defined for the five best conserved blocks and, when used as search motifs, are able to clearly distinguish between the m5C MTases and all other identified proteins in the PIR database. This suggests they may be of use in identifying putative MTases among protein sequences of unknown function.     

Crystal structure of RlmM, the 2��O-ribose methyltransferase for C2498 of Escherichia coli 23S rRNA

RlmM (YgdE) catalyzes the S-adenosyl methionine (AdoMet)-dependent 2′O methylation of C2498 in 23S ribosomal RNA (rRNA) of Escherichia coli. Previous experiments have shown that RlmM is active on 23S rRNA from an RlmM knockout strain but not on mature 50S subunits from the same strain. Here, we demonstrate RlmM methyltransferase (MTase) activity on in vitro transcribed 23S rRNA and its domain V. We have solved crystal structures of E. coli RlmM at 1.9 Å resolution and of an RlmM–AdoMet complex at 2.6 Å resolution. RlmM consists of an N-terminal THUMP domain and a C-terminal catalytic Rossmann-like fold MTase domain in a novel arrangement. The catalytic domain of RlmM is closely related to YiiB, TlyA and fibrillarins, with the second K of the catalytic tetrad KDKE shifted by two residues at the C-terminal end of a beta strand compared with most 2′O MTases. The AdoMet-binding site is open and shallow, suggesting that RNA substrate binding may be required to form a conformation needed for catalysis. A continuous surface of conserved positive charge indicates that RlmM uses one side of the two domains and the inter-domain linker to recognize its RNA substrate.     

Cellular localization and phosphorylation of Hrb1p is independent of Sky1p

Protein phosphorylation plays a major role in regulating cellular functions. We have previously demonstrated that Sky1p, the SR protein kinase of the budding yeast Saccharomyces cerevisiae, is a regulator of polyamine transport and ion homeostasis. Since its kinase activity was demonstrated essential for fulfilling these roles, we assumed that Sky1p function via substrates phosphorylation. Using an in vitro phosphorylation assay, we have identified Hrb1p as a putative Sky1p substrate. However, phosphorylation analysis in WT and sky1Delta cells and localization studies disproved Hrb1p as a true Sky1p substrate, although a segment of the RS domain is required for determining its subcellular localization. Furthermore, we demonstrate that Hrb1p and additional putative Sky1p substrates, identified by computational approach, are not involved in mediating the spermine tolerant phenotype of sky1Delta cells.     

Is the HemK family of putative S-adenosylmethionine-dependent methyltransferases a "missing" zeta subfamily of adenine methyltransferases? A hypothesis

Previous comparative studies revealed close similarity among various groups of S-adenosyl-L-methionine (AdoMet)-dependent methyltransferases (MTases), indicating their common evolutionary origin. We present evidence for a remarkable similarity between the sequence and predicted structure of HemK (a widespread family of putative proteins encoded in genomes from bacteria to humans) and the catalytic domain of the gamma-subfamily of adenine-specific DNA MTases (N6mA MTases). We predict the structure and function of the putative catalytic domain of HemK proteins and speculate that the target-recognizing function may be conferred by the N-terminal variable region.