Iterative database searches demonstrate that glycoside hydrolase families 27, 31, 36 and 66 share a common evolutionary origin with family 13

A catalytic sequence motif PDX10-30(E/D)XK is found in many restriction enzymes. On the basis of sequence similarities and mapping of the conserved residues to the crystal structure of NgoMIV we suggest that residues D160, K182, R186, R188 and E195 contribute to the catalytic/DNA binding site of the Ecl18kI restriction endonuclease. Mutational analysis confirms the functional significance of the conserved residues of Ecl18kI. Therefore, we conclude that the active site motif 159VDX21KX12E of Ecl18kI differs from the canonical PDX10-30(E/D)XK motif characteristic for most of the restriction enzymes. Moreover, we propose that two subfamilies of endonucleases Ecl18kI/PspGI/EcoRII and Cfr10I/Bse634I/NgoMIV, specific, respectively, for CCNGG/CCWGG and RCCGGY/GCCGGC sites, share conserved active site architecture and DNA binding elements.     

An Mrr-family nuclease motif in the single polypeptide restriction–modification enzyme LlaGI

Bioinformatic analysis of the putative nuclease domain of the single polypeptide restriction–modification enzyme LlaGI reveals amino acid motifs characteristic of the Escherichia coli methylated DNA-specific Mrr endonuclease. Using mutagenesis, we examined the role of the conserved residues in both DNA translocation and cleavage. Mutations in those residues predicted to play a role in DNA hydrolysis produced enzymes that could translocate on DNA but were either unable to cleave the polynucleotide track or had reduced nuclease activity. Cleavage by LlaGI is not targeted to methylated DNA, suggesting that the conserved motifs in the Mrr domain are a conventional sub-family of the PD-(D/E)XK superfamily of DNA nucleases.     

Unusual evolutionary history of the tRNA splicing endonuclease EndA: relationship to the LAGLIDADG and PD-(D/E)XK deoxyribonucleases

The tRNA splicing endoribonuclease EndA from Methanococcus jannaschii is a homotetramer formed via heterologous interaction between the two pairs of homodimers. Each monomer consists of two alpha/beta domains, the N-terminal domain (NTD) and the C-terminal domain (CTD) containing the RNase A-like active site. Comparison of the EndA coordinates with the publicly available protein structure database revealed the similarity of both domains to site-specific deoxyribonucleases: the NTD to the LAGLIDADG family and the CTD to the PD-(D/E)XK family. Superposition of the NTD on the catalytic domain of LAGLIDADG homing endonucleases allowed a suggestion to be made about which amino acid residues of the tRNA splicing nuclease might participate in formation of a presumptive cryptic deoxyribonuclease active site. On the other hand, the CTD and PD-(D/E)XK endonucleases, represented by restriction enzymes and a phage lambda exonuclease, were shown to share extensive similarities of the structural framework, to which entirely different active sites might be attached in two alternative locations. These findings suggest that EndA evolved from a fusion protein with at least two distinct endonuclease activities: the ribonuclease, which made it an essential "antitoxin" for the cells whose RNA genes were interrupted by introns, and the deoxyribonuclease, which provided the means for homing-like mobility. The residues of the noncatalytic CTDs from the positions corresponding to the catalytic side chains in PD-(D/E)XK deoxyribonucleases map to the surface at the opposite side to the tRNA binding site, for which no function has been implicated. Many restriction enzymes from the PD-(D/E)XK superfamily might have the potential to maintain an additional active or binding site at the face opposite the deoxyribonuclease active site, a property that can be utilized in protein engineering.     

Activity and Specificity of the Bacterial PD-(D/E)XK Homing Endonuclease I-Ssp6803I

The restriction endonuclease fold [a three-layer α-β sandwich containing variations of the PD-(D/E)XK nuclease motif] has been greatly diversified during evolution, facilitating its use for many biological functions. Here we characterize DNA binding and cleavage by the PD-(D/E)XK homing endonuclease I-Ssp6803I. Unlike most restriction endonucleases harboring the same core fold, the specificity profile of this enzyme extends over a long (17 bp) target site. The DNA binding and cleavage specificity profiles of this enzyme were independently determined and found to be highly correlated. However, the DNA target sequence contains several positions where binding and cleavage activities are not tightly coupled: individual DNA base-pair substitutions at those positions that significantly decrease cleavage activity have minor effects on binding affinity. These changes in the DNA target sequence appear to correspond to substitutions that uniquely increase the free energy change between the ground state and the transition state, rather than simply decreasing the overall DNA binding affinity. The specificity of the enzyme reflects constraints on its host gene and limitations imposed by the enzyme's quaternary structure and illustrate the highly diverse repertoire of DNA recognition specificities that can be adopted by the related folds surrounding the PD-(D/E)XK nuclease motif.     

Mva1269I: a monomeric type IIS restriction endonuclease from Micrococcus varians with two EcoRI- and FokI-like catalytic domains

Type II restriction endonuclease Mva1269I recognizes an asymmetric DNA sequence 5'-GAATGCN / -3'/5'-NG / CATTC-3' and cuts top and bottom DNA strands at positions, indicated by the "/" symbol. Most restriction endonucleases require dimerization to cleave both strands of DNA. We found that Mva1269I is a monomer both in solution and upon binding of cognate DNA. Protein fold-recognition analysis revealed that Mva1269I comprises two "PD-(D/E)XK" domains. The N-terminal domain is related to the 5'-GAATTC-3'-specific restriction endonuclease EcoRI, whereas the C-terminal one resembles the nonspecific nuclease domain of restriction endonuclease FokI. Inactivation of the C-terminal catalytic site transformed Mva1269I into a very active bottom strand-nicking enzyme, whereas mutants in the N-terminal domain nicked the top strand, but only at elevated enzyme concentrations. We found that the cleavage of the bottom strand is a prerequisite for the cleavage of the top strand. We suggest that Mva1269I evolved the ability to recognize and to cleave its asymmetrical target by a fusion of an EcoRI-like domain, which incises the bottom strand within the target, and a FokI-like domain that completes the cleavage within the nonspecific region outside the target sequence. Our results have implications for the molecular evolution of restriction endonucleases, as well as for perspectives of engineering new restriction and nicking enzymes with asymmetric target sites.     

Identification of a single HNH active site in type IIS restriction endonuclease Eco31I

Type IIS restriction endonuclease Eco31I is a "short-distance cutter", which cleaves DNA strands close to its recognition sequence, 5'-GGTCTC(1/5). Previously, it has been proposed that related endonucleases recognizing a common sequence core GTCTC possess two active sites for cleavage of both strands in the DNA substrate. Here, we present bioinformatic identification and experimental evidence for a single nuclease active site. We identified a short region of homology between Eco31I and HNH nucleases, constructed a three-dimensional model of the putative catalytic domain and validated our predictions by random and site-specific mutagenesis. The restriction mechanism of Eco31I is suggested by analogy to the mechanisms of phage T4 endonuclease VII and homing endonuclease I-PpoI. We propose that residues D311 and N334 coordinate the cofactor. H312 acts as a general base-activating water molecule for the nucleophilic attack. K337 together with R340 and D345 are located in close proximity to the active center and are essential for correct folding of catalytic motif, while D345 together with R264 and D273 could be directly involved in DNA binding. We also predict that the Eco31I catalytic domain contains a putative Zn-binding site, which is essential for its structural integrity. Our results suggest that the HNH-like active site is involved in the cleavage of both strands in the DNA substrate. On the other hand, analysis of site-specific mutants in the region, previously suggested to harbor the second active site, revealed its irrelevance to the nuclease activity. Thus, our data argue against the earlier prediction and indicate the presence of a single conserved active site in type IIS restriction endonucleases that recognize common sequence core GTCTC.     

Protein stability indicates divergent evolution of PD-(D/E)XK type II restriction endonucleases

Type II restriction endonucleases recognize 4-8 base-pair-long DNA sequences and catalyze their cleavage with remarkable specificity. Crystal structures of the PD-(DE)XK superfamily revealed a common alpha/beta core motif and similar active site. In contrast, these enzymes show little sequence similarity and use different strategies to interact with their substrate DNA. The intriguing question is whether this enzyme family could have evolved from a common origin. In our present work, protein structure stability elements were analyzed and compared in three parts of PD-(DE)XK type II restriction endonucleases: (1) core motif, (2) active-site residues, and (3) residues playing role in DNA recognition. High correlation was found between the active-site residues and those stabilization factors that contribute to preventing structural decay. DNA recognition sites were also observed to participate in stabilization centers. It indicates that recognition motifs and active sites in PD-(DE)XK type II restriction endonucleases should have been evolutionary more conserved than other parts of the structure. Based on this observation it is proposed that PD-(DE)XK type II restriction endonucleases have developed from a common ancestor with divergent evolution.     

I-Ssp6803I: the first homing endonuclease from the PD-(D/E)XK superfamily exhibits an unusual mode of DNA recognition

MOTIVATION: Restriction endonucleases (REases) and homing endonucleases (HEases) are biotechnologically important enzymes. Nearly all structurally characterized REases belong to the PD-(D/E)XK superfamily of nucleases, while most HEases belong to an unrelated LAGLIDADG superfamily. These two protein folds are typically associated with very different modes of protein-DNA recognition, consistent with the different mechanisms of action required to achieve high specificity. REases recognize short DNA sequences using multiple contacts per base pair, while HEases recognize very long sites using a few contacts per base pair, thereby allowing for partial degeneracy of the target sequence. Thus far, neither REases with the LAGLIDADG fold, nor HEases with the PD-(D/E)XK fold, have been found. RESULTS: Using protein fold recognition, we have identified the first member of the PD-(D/E)XK superfamily among homing endonucleases, a cyanobacterial enzyme I-Ssp6803I. We present a model of the I-Ssp6803I-DNA complex based on the structure of Type II restriction endonuclease R.BglI and predict the active site and residues involved in specific DNA sequence recognition by I-Ssp6803I. Our finding reveals a new unexpected evolutionary link between HEases and REases and suggests how PD-(D/E)XK nucleases may develop a 'HEase-like' way of interacting with the extended DNA sequence. This in turn may be exploited to study the evolution of DNA sequence specificity and to engineer nucleases with new substrate specificities.     

Tetrameric restriction enzymes: expansion to the GIY-YIG nuclease family

The GIY-YIG nuclease domain was originally identified in homing endonucleases and enzymes involved in DNA repair and recombination. Many of the GIY-YIG family enzymes are functional as monomers. We show here that the Cfr42I restriction endonuclease which belongs to the GIY-YIG family and recognizes the symmetric sequence 5′-CCGC/GG-3′ (‘/’ indicates the cleavage site) is a tetramer in solution. Moreover, biochemical and kinetic studies provided here demonstrate that the Cfr42I tetramer is catalytically active only upon simultaneous binding of two copies of its recognition sequence. In that respect Cfr42I resembles the homotetrameric Type IIF restriction enzymes that belong to the distinct PD-(E/D)XK nuclease superfamily. Unlike the PD-(E/D)XK enzymes, the GIY-YIG nuclease Cfr42I accommodates an extremely wide selection of metal-ion cofactors, including Mg²⁺, Mn²⁺, Co²⁺, Zn²⁺, Ni²⁺, Cu²⁺ and Ca²⁺. To our knowledge, Cfr42I is the first tetrameric GIY-YIG family enzyme. Similar structural arrangement and phenotypes displayed by restriction enzymes of the PD-(E/D)XK and GIY-YIG nuclease families point to the functional significance of tetramerization.     

A mutation in the Mod subunit of EcoP15I restriction enzyme converts the DNA methyltransferase to a site-specific endonuclease

A closer inspection of the amino acid sequence of EcoP15I DNA methyltransferase revealed a region of similarity to the PDXn(D/E)XK catalytic site of type II restriction endonucleases, except for methionine in EcoP15I DNA methyltransferase instead of proline. Substitution of methionine at position 357 by proline converts EcoP15I DNA methyltransferase to a site-specific endonuclease. EcoP15I-M357P DNA methyltransferase specifically binds to the recognition sequence 5'-CAGCAG-3' and cleaves DNA asymmetrically EcoP151-M357P.DNA methyltransferase specifically binds to the recognition sequence 5'-CAGCAG-3' and cleaves DNA asymmetrically, 5'-CAGCAG(N)(10)-3', as indicated by the arrows, in presence of magnesium ions.     

Identification of new homologs of PD-(D/E)XK nucleases by support vector machines trained on data derived from profile-profile alignments

PD-(D/E)XK nucleases, initially represented by only Type II restriction enzymes, now comprise a large and extremely diverse superfamily of proteins. They participate in many different nucleic acids transactions including DNA degradation, recombination, repair and RNA processing. Different PD-(D/E)XK families, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. This makes the identification of new superfamily members using standard homology search techniques challenging. To tackle this problem, we developed a method for the detection of PD-(D/E)XK families based on the binary classification of profile-profile alignments using support vector machines (SVMs). Using a number of both superfamily-specific and general features, SVMs were trained to identify true positive alignments of PD-(D/E)XK representatives. With this method we identified several PFAM families of uncharacterized proteins as putative new members of the PD-(D/E)XK superfamily. In addition, we assigned several unclassified restriction enzymes to the PD-(D/E)XK type. Results show that the new method is able to make confident assignments even for alignments that have statistically insignificant scores. We also implemented the method as a freely accessible web server at http://www.ibt.lt/bioinformatics/software/pdexk/.     

Mutability of an HNH nuclease imidazole general base and exchange of a deprotonation mechanism

Several unique protein folds that catalyze the hydrolysis of phosphodiester bonds have arisen independently in nature, including the PD(D/E)XK superfamily (typified by type II restriction endonucleases and many recombination and repair enzymes) and the HNH superfamily (found in an equally wide array of enzymes, including bacterial colicins and homing endonucleases). Whereas the identity and position of catalytic residues within the PD(D/E)XK superfamily are highly variable, the active sites of HNH nucleases are much more strongly conserved. In this study, the ability of an HNH nuclease to tolerate a mutation of its most conserved catalytic residue (its histidine general base), and the mechanism of the most active enzyme variant, were characterized. Conversion of this residue into several altered chemistries, glutamine, lysine, or glutamate, resulted in measurable activity. The histidine to glutamine mutant displays the highest residual activity and a pH profile similar to that of the wild-type enzyme. This activity is dependent on the presence of a neighboring imidazole ring, which has taken over as a less efficient general base for the reaction. This result implies that mutational pathways to alternative HNH-derived catalytic sites do exist but are not as extensively or successfully diverged or reoptimized in nature as variants of the PD(D/E)XK nuclease superfamily. This is possibly due to multiple steric constraints placed on the compact HNH motif, which is simultaneously involved in protein folding, DNA binding, and catalysis, as well as the use of a planar, aromatic imidazole group as a general base.     

BspRI restriction endonuclease: cloning, expression in Escherichia coli and sequential cleavage mechanism

The GGCC-specific restriction endonuclease BspRI is one of the few Type IIP restriction endonucleases, which were suggested to be a monomer. Amino acid sequence information obtained by Edman sequencing and mass spectrometry analysis was used to clone the gene encoding BspRI. The bspRIR gene is located adjacently to the gene of the cognate modification methyltransferase and encodes a 304 aa protein. Expression of the bspRIR gene in Escherichia coli was dependent on the replacement of the native TTG initiation codon with an ATG codon, explaining previous failures in cloning the gene using functional selection. A plasmid containing a single BspRI recognition site was used to analyze kinetically nicking and second-strand cleavage under steady-state conditions. Cleavage of the supercoiled plasmid went through a relaxed intermediate indicating sequential hydrolysis of the two strands. Results of the kinetic analysis of the first- and second-strand cleavage are consistent with cutting the double-stranded substrate site in two independent binding events. A database search identified eight putative restriction-modification systems in which the predicted endonucleases as well as the methyltransferases share high sequence similarity with the corresponding protein of the BspRI system. BspRI and the related putative restriction endonucleases belong to the PD-(D/E)XK nuclease superfamily.     

The crystal structure of the rare-cutting restriction enzyme SdaI reveals unexpected domain architecture

Rare-cutting restriction enzymes are important tools in genome analysis. We report here the crystal structure of SdaI restriction endonuclease, which is specific for the 8 bp sequence CCTGCA/GG ("/" designates the cleavage site). Unlike orthodox Type IIP enzymes, which are single domain proteins, the SdaI monomer is composed of two structural domains. The N domain contains a classical winged helix-turn-helix (wHTH) DNA binding motif, while the C domain shows a typical restriction endonuclease fold. The active site of SdaI is located within the C domain and represents a variant of the canonical PD-(D/E)XK motif. SdaI determinants of sequence specificity are clustered on the recognition helix of the wHTH motif at the N domain. The modular architecture of SdaI, wherein one domain mediates DNA binding while the other domain is predicted to catalyze hydrolysis, distinguishes SdaI from previously characterized restriction enzymes interacting with symmetric recognition sequences.     

Type II restriction endonuclease R.Eco29kI is a member of the GIY-YIG nuclease superfamily

Background

The majority of experimentally determined crystal structures of Type II restriction endonucleases (REases) exhibit a common PD-(D/E)XK fold. Crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI), and bioinformatics analyses supported by mutagenesis suggested that some REases belong to the HNH fold. Our previous bioinformatic analysis suggested that REase R.Eco29kI shares sequence similarities with one more unrelated nuclease superfamily, GIY-YIG, however so far no experimental data were available to support this prediction. The determination of a crystal structure of the GIY-YIG domain of homing endonuclease I-TevI provided a template for modeling of R.Eco29kI and prompted us to validate the model experimentally.

Results

Using protein fold-recognition methods we generated a new alignment between R.Eco29kI and I-TevI, which suggested a reassignment of one of the putative catalytic residues. A theoretical model of R.Eco29kI was constructed to illustrate its predicted three-dimensional fold and organization of the active site, comprising amino acid residues Y49, Y76, R104, H108, E142, and N154. A series of mutants was constructed to generate amino acid substitutions of selected residues (Y49A, R104A, H108F, E142A and N154L) and the mutant proteins were examined for their ability to bind the DNA containing the Eco29kI site 5'-CCGCGG-3' and to catalyze the cleavage reaction. Experimental data reveal that residues Y49, R104, E142, H108, and N154 are important for the nuclease activity of R.Eco29kI, while H108 and N154 are also important for specific DNA binding by this enzyme.

Conclusion

Substitutions of residues Y49, R104, H108, E142 and N154 predicted by the model to be a part of the active site lead to mutant proteins with strong defects in the REase activity. These results are in very good agreement with the structural model presented in this work and with our prediction that R.Eco29kI belongs to the GIY-YIG superfamily of nucleases. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD-(D/E)XK or HNH superfamilies of nucleases, and is instead a member of the unrelated GIY-YIG superfamily.     

The restriction fold turns to the dark side: a bacterial homing endonuclease with a PD-(D/E)-XK motif

The homing endonuclease I-Ssp6803I causes the insertion of a group I intron into a bacterial tRNA gene-the only example of an invasive mobile intron within a bacterial genome. Using a computational fold prediction, mutagenic screen and crystal structure determination, we demonstrate that this protein is a tetrameric PD-(D/E)-XK endonuclease - a fold normally used to protect a bacterial genome from invading DNA through the action of restriction endonucleases. I-Ssp6803I uses its tetrameric assembly to promote recognition of a single long target site, whereas restriction endonuclease tetramers facilitate cooperative binding and cleavage of two short sites. The limited use of the PD-(D/E)-XK nucleases by mobile introns stands in contrast to their frequent use of LAGLIDADG and HNH endonucleases - which in turn, are rarely incorporated into restriction/modification systems.     

Mutational analysis of active‐site residues in the Mycobacterium leprae RecA intein,a LAGLIDADG homing endonuclease: Asp122 and Asp193 are crucial to the double‐stranded DNA cleavage activity whereas Asp218 is not

Mycobacterium leprae recA harbors an in‐frame insertion sequence that encodes an intein homing endonuclease (PI‐MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their active center. A common feature of LAGLIDADG‐type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI‐MleI is distinctive from other members of the family of LAGLIDADG‐type HEases for its modular structure with functionally separable domains for DNA‐binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI‐MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active‐site residues essential for DNA target site recognition and double‐stranded DNA cleavage, we performed site‐directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild‐type PI‐MleI and its variants disclosed that the two amino acid residues, Asp¹²² (in Block C) and Asp¹⁹³ (in functional Block E), are crucial to the double‐stranded DNA endonuclease activity, whereas Asp²¹⁸ (in pseudo‐Block E) is not. However, despite the reduced catalytic activity, the PI‐MleI variants, like the wild‐type PI‐MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA‐binding affinities, but abolished the double‐stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double‐stranded DNA cleavage activity, compared with the wild‐type PI‐MleI. These results provide compelling evidence that Asp¹²² and Asp¹⁹³ in DOD motif I and II, respectively, are bona fide active‐site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.     

Structural characterization of MepB from Staphylococcus aureus reveals homology to endonucleases

The MepRAB operon in Staphylococcus aureus has been identified to play a role in drug resistance. Although the functions of MepA and MepR are known, little information is available on the function of MepB. Here we report the X‐ray structure of MepB to 2.1 Å revealing its structural similarity to the PD‐(D/E)XK family of endonucleases. We further show that MepB binds DNA and RNA, with a higher affinity towards RNA and single stranded DNA than towards double stranded DNA. Notably, the PD‐(D/E)XK catalytic active site residues are not conserved in MepB. MepB's association with a drug resistance operon suggests that it plays a role in responding to antimicrobials. This role is likely carried out through MepB's interactions with nucleic acids.     

Structural analysis of the heterodimeric type IIS restriction endonuclease R.BspD6I acting as a complex between a monomeric site-specific nickase and a catalytic subunit

The heterodimeric restriction endonuclease R.BspD6I from Bacillus species D6 recognizes a pseudosymmetric sequence and cuts both DNA strands outside the recognition sequence. The large subunit, Nt.BspD6I, acts as a type IIS site-specific monomeric nicking endonuclease. The isolated small subunit, ss.BspD6I, does not bind DNA and is not catalytically active. We solved the crystal structures of Nt.BspD6I and ss.BspD6I at high resolution. Nt.BspD6I consists of three domains, two of which exhibit structural similarity to the recognition and cleavage domains of FokI. ss.BspD6I has a fold similar to that of the cleavage domain of Nt.BspD6I, each containing a PD-(D/E)XK motif and a histidine as an additional putative catalytic residue. In contrast to the DNA-bound FokI structure, in which the cleavage domain is rotated away from the DNA, the crystal structure of Nt.BspD6I shows the recognition and cleavage domains in favorable orientations for interactions with DNA. Docking models of complexes of Nt.BspD6I and R.BspD6I with cognate DNA were constructed on the basis of structural similarity to individual domains of FokI, R.BpuJI and HindIII. A three-helix bundle forming an interdomain linker in Nt.BspD6I acts as a rigid spacer adjusting the orientations of the spatially separated domains to match the distance between the recognition and cleavage sites accurately.