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.     

A model of restriction endonuclease MvaI in complex with DNA: a template for interpretation of experimental data and a guide for specificity engineering

R.MvaI is a Type II restriction enzyme (REase), which specifically recognizes the pentanucleotide DNA sequence 5'-CCWGG-3' (W indicates A or T). It belongs to a family of enzymes, which recognize related sequences, including 5'-CCSGG-3' (S indicates G or C) in the case of R.BcnI, or 5'-CCNGG-3' (where N indicates any nucleoside) in the case of R.ScrFI. REases from this family hydrolyze the phosphodiester bond in the DNA between the 2nd and 3rd base in both strands, thereby generating a double strand break with 5'-protruding single nucleotides. So far, no crystal structures of REases with similar cleavage patterns have been solved. Characterization of sequence-structure-function relationships in this family would facilitate understanding of evolution of sequence specificity among REases and could aid in engineering of enzymes with new specificities. However, sequences of R.MvaI or its homologs show no significant similarity to any proteins with known structures, thus precluding straightforward comparative modeling. We used a fold recognition approach to identify a remote relationship between R.MvaI and the structure of DNA repair enzyme MutH, which belongs to the PD-(D/E)XK superfamily together with many other REases. We constructed a homology model of R.MvaI and used it to predict functionally important amino acid residues and the mode of interaction with the DNA. In particular, we predict that only one active site of R.MvaI interacts with the DNA target at a time, and the cleavage of both strands (5'-CCAGG-3' and 5'-CCTGG-3') is achieved by two independent catalytic events. The model is in good agreement with the available experimental data and will serve as a template for further analyses of R.MvaI, R.BcnI, R.ScrFI and other related enzymes.     

Structural and evolutionary classification of Type II restriction enzymes based on theoretical and experimental analyses

For a very long time, Type II restriction enzymes (REases) have been a paradigm of ORFans: proteins with no detectable similarity to each other and to any other protein in the database, despite common cellular and biochemical function. Crystallographic analyses published until January 2008 provided high-resolution structures for only 28 of 1637 Type II REase sequences available in the Restriction Enzyme database (REBASE). Among these structures, all but two possess catalytic domains with the common PD-(D/E)XK nuclease fold. Two structures are unrelated to the others: R.BfiI exhibits the phospholipase D (PLD) fold, while R.PabI has a new fold termed 'half-pipe'. Thus far, bioinformatic studies supported by site-directed mutagenesis have extended the number of tentatively assigned REase folds to five (now including also GIY-YIG and HNH folds identified earlier in homing endonucleases) and provided structural predictions for dozens of REase sequences without experimentally solved structures. Here, we present a comprehensive study of all Type II REase sequences available in REBASE together with their homologs detectable in the nonredundant and environmental samples databases at the NCBI. We present the summary and critical evaluation of structural assignments and predictions reported earlier, new classification of all REase sequences into families, domain architecture analysis and new predictions of three-dimensional folds. Among 289 experimentally characterized (not putative) Type II REases, whose apparently full-length sequences are available in REBASE, we assign 199 (69%) to contain the PD-(D/E)XK domain. The HNH domain is the second most common, with 24 (8%) members. When putative REases are taken into account, the fraction of PD-(D/E)XK and HNH folds changes to 48% and 30%, respectively. Fifty-six characterized (and 521 predicted) REases remain unassigned to any of the five REase folds identified so far, and may exhibit new architectures. These enzymes are proposed as the most interesting targets for structure determination by high-resolution experimental methods. Our analysis provides the first comprehensive map of sequence-structure relationships among Type II REases and will help to focus the efforts of structural and functional genomics of this large and biotechnologically important class of enzymes.     

A homology model of restriction endonuclease SfiI in complex with DNA

Background

Restriction enzymes (REases) are commercial reagents commonly used in recombinant DNA technologies. They are attractive models for studying protein-DNA interactions and valuable targets for protein engineering. They are, however, extremely divergent: the amino acid sequence of a typical REase usually shows no detectable similarities to any other proteins, with rare exceptions of other REases that recognize identical or very similar sequences. From structural analyses and bioinformatics studies it has been learned that some REases belong to at least four unrelated and structurally distinct superfamilies of nucleases, PD-DxK, PLD, HNH, and GIY-YIG. Hence, they are extremely hard targets for structure prediction and homology-based inference of sequence-function relationships and the great majority of REases remain structurally and evolutionarily unclassified.

Results

SfiI is a REase which recognizes the interrupted palindromic sequence 5'GGCCNNNN^NGGCC3' and generates 3 nt long 3' overhangs upon cleavage. SfiI is an archetypal Type IIF enzyme, which functions as a tetramer and cleaves two copies of the recognition site in a concerted manner. Its sequence shows no similarity to other proteins and nothing is known about the localization of its active site or residues important for oligomerization. Using the threading approach for protein fold-recognition, we identified a remote relationship between SfiI and BglI, a dimeric Type IIP restriction enzyme from the PD-DxK superfamily of nucleases, which recognizes the 5'GCCNNNN^NGGC3' sequence and whose structure in complex with the substrate DNA is available. We constructed a homology model of SfiI in complex with its target sequence and used it to predict residues important for dimerization, tetramerization, DNA binding and catalysis.

Conclusions

The bioinformatics analysis suggest that SfiI, a Type IIF enzyme, is more closely related to BglI, an "orthodox" Type IIP restriction enzyme, than to any other REase, including other Type IIF REases with known structures, such as NgoMIV. NgoMIV and BglI belong to two different, very remotely related branches of the PD-DxK superfamily: the α-class (EcoRI-like), and the β-class (EcoRV-like), respectively. Thus, our analysis provides evidence that the ability to tetramerize and cut the two DNA sequences in a concerted manner was developed independently at least two times in the evolution of the PD-DxK superfamily of REases. The model of SfiI will also serve as a convenient platform for further experimental analyses.     

Restriction endonuclease BpuJI specific for the 5'-CCCGT sequence is related to the archaeal Holliday junction resolvase family

Type IIS restriction endonucleases (REases) recognize asymmetric DNA sequences and cleave both DNA strands at fixed positions downstream of the recognition site. REase BpuJI recognizes the asymmetric sequence 5′-CCCGT, however it cuts at multiple sites in the vicinity of the target sequence. We show that BpuJI is a dimer, which has two DNA binding surfaces and displays optimal catalytic activity when bound to two recognition sites. BpuJI is cleaved by chymotrypsin into an N-terminal domain (NTD), which lacks catalytic activity but binds specifically to the recognition sequence as a monomer, and a C-terminal domain (CTD), which forms a dimer with non-specific nuclease activity. Fold recognition approach reveals that the CTD of BpuJI is structurally related to archaeal Holliday junction resolvases (AHJR). We demonstrate that the isolated catalytic CTD of BpuJI possesses end-directed nuclease activity and preferentially cuts 3nt from the 3′-terminus of blunt-ended DNA. The nuclease activity of the CTD is repressed in the apo-enzyme and becomes activated upon specific DNA binding by the NTDs. This leads to a complicated pattern of specific DNA cleavage in the vicinity of the target site. Bioinformatics analysis identifies the AHJR-like domain in the putative Type III enzymes and functionally uncharacterized proteins.     

Dual role for Zn2+ in maintaining structural integrity and inducing DNA sequence specificity in a promiscuous endonuclease

We describe two uncommon roles for Zn2+ in enzyme KpnI restriction endonuclease (REase). Among all of the REases studied, KpnI REase is unique in its DNA binding and cleavage characteristics. The enzyme is a poor discriminator of DNA sequences, cleaving DNA in a promiscuous manner in the presence of Mg2+. Unlike most Type II REases, the active site of the enzyme comprises an HNH motif, which can accommodate Mg2+, Mn2+, or Ca2+. Among these metal ions, Mg2+ and Mn2+ induce promiscuous cleavage by the enzyme, whereas Ca2+-bound enzyme exhibits site-specific cleavage. Examination of the sequence of the protein revealed the presence of a zinc finger CCCH motif rarely found in proteins of prokaryotic origin. The zinc binding motif tightly coordinates zinc to provide a rigid structural framework for the enzyme needed for its function. In addition to this structural scaffold, another atom of zinc binds to the active site to induce high fidelity cleavage and suppress the Mg2+- and Mn2+-mediated promiscuous behavior of the enzyme. This is the first demonstration of distinct structural and catalytic roles for zinc in an enzyme, suggesting the distinct origin of KpnI REase.     

Type II restriction endonuclease R.KpnI is a member of the HNH nuclease superfamily

The restriction endonuclease (REase) R.KpnI is an orthodox Type IIP enzyme, which binds to DNA in the absence of metal ions and cleaves the DNA sequence 5′-GGTAC^C-3′ in the presence of Mg²⁺ as shown generating 3′ four base overhangs. Bioinformatics analysis reveals that R.KpnI contains a ββα-Me-finger fold, which is characteristic of many HNH-superfamily endonucleases, including homing endonuclease I-HmuI, structure-specific T4 endonuclease VII, colicin E9, sequence non-specific Serratia nuclease and sequence-specific homing endonuclease I-PpoI. According to our homology model of R.KpnI, D148, H149 and Q175 correspond to the critical D, H and N or H residues of the HNH nucleases. Substitutions of these three conserved residues lead to the loss of the DNA cleavage activity by R.KpnI, confirming their importance. The mutant Q175E fails to bind DNA at the standard conditions, although the DNA binding and cleavage can be rescued at pH 6.0, indicating a role for Q175 in DNA binding and cleavage. Our study provides the first experimental evidence for a Type IIP REase that does not belong to the PD…D/EXK superfamily of nucleases, instead is a member of the HNH superfamily.     

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.     

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.     

Structural integrity of the beta beta alpha-Metal finger motif is required for DNA binding and stable protein-DNA complex formation in R.KpnI

Restriction endonuclease (REase) R.KpnI from Klebsiella pneumoniae is a homodimeric enzyme, which recognizes palindromic sequence GGTAC|C and cleaves generating 4 base 3' end overhangs. R.KpnI belongs to the HNH superfamily of nucleases, which are characterized by the presence of the beta beta alpha-Me finger motif. Structurally, this motif consists of a twisted beta-hairpin followed by an alpha-helix, and serves as a scaffold for side chains of residues involved in co-ordination of a divalent metal ion that is required for catalysis. Homology modeling studies of R.KpnI suggested a crossover structure for the alpha-helix, which could possibly form dimeric interface and/or structural scaffold for the active site. We have evaluated the role of the residues present in this alpha-helix in intersubunit interactions and/or stabilization of the active site. We show here that mutations of residues in the alpha-helix lead to a loss of the enzyme activity, but not dimerization ability. Intrinsic fluorescence and circular dichroism studies revealed that the loss of function phenotype was due to the structural perturbation of the beta beta alpha-Me finger motif. The results of mutational analysis suggest that the alpha-helix of the beta beta alpha-Me finger of R.KpnI plays an important role for the stability of the protein-DNA complex.     

The structure and enzymatic properties of a novel RNase II family enzyme from Deinococcus radiodurans

Exoribonucleases are vital in nearly all aspects of RNA metabolism, including RNA maturation, end-turnover, and degradation. RNase II and RNase R are paralogous members of the RNR superfamily of nonspecific, 3'→5', processive exoribonucleases. In Escherichia coli, RNase II plays a primary role in mRNA decay and has a preference for unstructured RNA. RNase R, in contrast, is capable of digesting structured RNA and plays a role in the degradation of both mRNA and stable RNA. Deinococcus radiodurans, a radiation-resistant bacterium, contains two RNR family members. The shorter of these, DrR63, includes a sequence signature typical of RNase R, but we show here that this enzyme is an RNase II-type exonuclease and cannot degrade structured RNA. We also report the crystal structure of this protein, now termed DrII. The DrII structure reveals a truncated RNA binding region in which the N-terminal cold shock domains, typical of most RNR family nucleases, are replaced by an unusual winged helix-turn-helix domain, where the "wing" is contributed by the C-terminal S1 domain. Consistent with its truncated RNA binding region, DrII is able to remove 3' overhangs from RNA molecules closer to duplexes than do other RNase II-type enzymes. DrII also displays distinct sensitivity to pyrimidine-rich regions of single-stranded RNA and is able to process tRNA precursors with adenosine-rich 3' extensions in vitro. These data indicate that DrII is the RNase II of D. radiodurans and that its structure and catalytic properties are distinct from those of other related enzymes.     

Swapping the domains of exoribonucleases RNase II and RNase R: conferring upon RNase II the ability to degrade ds RNA

RNase II and RNase R are the two E. coli exoribonucleases that belong to the RNase II super family of enzymes. They degrade RNA hydrolytically in the 3' to 5' direction in a processive and sequence independent manner. However, while RNase R is capable of degrading structured RNAs, the RNase II activity is impaired by dsRNAs. The final end-product of these two enzymes is also different, being 4 nt for RNase II and 2 nt for RNase R. RNase II and RNase R share structural properties, including 60% of amino acid sequence similarity and have a similar modular domain organization: two N-terminal cold shock domains (CSD1 and CSD2), one central RNB catalytic domain, and one C-terminal S1 domain. We have constructed hybrid proteins by swapping the domains between RNase II and RNase R to determine which are the responsible for the differences observed between RNase R and RNase II. The results obtained show that the S1 and RNB domains from RNase R in an RNase II context allow the degradation of double-stranded substrates and the appearance of the 2 nt long end-product. Moreover, the degradation of structured RNAs becomes tail-independent when the RNB domain from RNase R is no longer associated with the RNA binding domains (CSD and S1) of the genuine protein. Finally, we show that the RNase R C-terminal Lysine-rich region is involved in the degradation of double-stranded substrates in an RNase II context, probably by unwinding the substrate before it enters into the catalytic cavity.     

Elimination of inter-domain interactions increases the cleavage fidelity of the restriction endonuclease DraIII

DraIII is a type IIP restriction endonucleases (REases) that recognizes and creates a double strand break within the gapped palindromic sequence CAC↑NNN↓GTG of double-stranded DNA (↑ indicates nicking on the bottom strand; ↓ indicates nicking on the top strand). However, wild type DraIII shows significant star activity. In this study, it was found that the prominent star site is CAT↑GTT↓GTG, consisting of a star 5′ half (CAT) and a canonical 3′ half (GTG). DraIII nicks the 3′ canonical half site at a faster rate than the 5′ star half site, in contrast to the similar rate with the canonical full site. The crystal structure of the DraIII protein was solved. It indicated, as supported by mutagenesis, that DraIII possesses a ββα-metal HNH active site. The structure revealed extensive intra-molecular interactions between the N-terminal domain and the C-terminal domain containing the HNH active site. Disruptions of these interactions through sitedirected mutagenesis drastically increased cleavage fidelity. The understanding of fidelity mechanisms will enable generation of high fidelity REases.     

Hin4II, a new prototype restriction endonuclease from Haemophilus influenzae RFL4: discovery, cloning and expression in Escherichia coli

The genes encoding restriction-modification system of unknown specificity Hin4II from Haemophilus influenzae RFL4 were cloned in Escherichia coli and sequenced. The Hin4II system comprises three tandemly arranged genes coding for m6A DNA methyltransferase, m5C DNA methyltransferase and restriction endonuclease, respectively. Restriction endonuclease was expressed in E. coli and purified to apparent homogeneity. The DNA recognition sequence and cleavage positions were determined. R.Hin4II recognizes the novel non-palindromic sequence 5'-CCTTC-3' and cleaves the DNA 6 and 5 nt downstream in the top and bottom strand, respectively. The new prototype restriction endonuclease Hin4II was classified as a potential candidate of HNH nuclease family after comparison against SMART database. An amino acid sequence motif 297H-X14-N-X8-H of Hin4II was proposed as forming a putative catalytic center.     

Identification of GATC- and CCGG-recognizing Type II REases and their putative specificity-determining positions using Scan2S--a novel motif scan algorithm with optional secondary structure constraints

Restriction endonucleases (REases) are DNA-cleaving enzymes that have become indispensable tools in molecular biology. Type II REases are highly divergent in sequence despite their common structural core, function and, in some cases, common specificities towards DNA sequences. This makes it difficult to identify and classify them functionally based on sequence, and has hampered the efforts of specificity-engineering. Here, we define novel REase sequence motifs, which extend beyond the PD-(D/E)XK hallmark, and incorporate secondary structure information. The automated search using these motifs is carried out with a newly developed fast regular expression matching algorithm that accommodates long patterns with optional secondary structure constraints. Using this new tool, named Scan2S, motifs derived from REases with specificity towards GATC- and CGGG-containing DNA sequences successfully identify REases of the same specificity. Notably, some of these sequences are not identified by standard sequence detection tools. The new motifs highlight potential specificity-determining positions that do not fully overlap for the GATC- and the CCGG-recognizing REases and are candidates for specificity re-engineering.     

A cDNA cloned from Physarum polycephalum encodes new type of family 3 beta-glucosidase that is a fusion protein containing a calx-beta motif

The microplasmodia of Physarum polycephalum express three types of β-glucosidases: secretory enzyme, a soluble cytoplasmic enzyme and a membrane-bound enzyme. We are interested in the physiological role of three enzymes. We report the sequence of cDNA for membrane β-glucosidase 1, which consists of 3825 nucleotides that includes an open reading frame encoding 1248 amino acids. The molecular weight of membrane β-glucosidase 1 was calculated to be 131,843 based on the predicted amino acid composition. Glycosyl hydrolase family 3 N-terminal and C-terminal domains were found within the N-terminal half of the membrane β-glucosidase 1 sequence and were highly homologous with the primary structures of fungal β-glucosidases. Notably, the C-terminal half of membrane β-glucosidase 1 contains two calx-β motifs, which are known to be Ca²⁺ binding domains in the Drosophila Na⁺/Ca²⁺ exchanger; an RGD sequence, which is known to be a cell attachment sequence; and a transmembrane region. In this way, Physarum membrane β-glucosidase 1 differs from all previously identified family 3 β-glucosidases. In addition to cDNA for membrane β-glucosidase 1, two other distinctly different mRNAs were also isolated. Two sequences were largely identical to cDNA for membrane β-glucosidase 1, but included a long insert sequence having a stop codon, leading to truncation of their products, which could account for other β-glucosidase forms occurred in Physarum poycephalum.

Thus, the membrane β-glucosidase is a new type family 3 enzyme fused with the Calx-β domain. We propose that Calx-β domain may modulate the β-glucosidase activity in response to changes in the Ca²⁺ concentration.     

Inference of relationships in the 'twilight zone' of homology using a combination of bioinformatics and site-directed mutagenesis: a case study of restriction endonucleases Bsp6I and PvuII

Thus far, identification of functionally important residues in Type II restriction endonucleases (REases) has been difficult using conventional methods. Even though known REase structures share a fold and marginally recognizable active site, the overall sequence similarities are statistically insignificant, unless compared among proteins that recognize identical or very similar sequences. Bsp6I is a Type II REase, which recognizes the palindromic DNA sequence 5′GCNGC and cleaves between the cytosine and the unspecified nucleotide in both strands, generating a double-strand break with 5′-protruding single nucleotides. There are no solved structures of REases that recognize similar DNA targets or generate cleavage products with similar characteristics. In straightforward comparisons, the Bsp6I sequence shows no significant similarity to REases with known structures. However, using a fold-recognition approach, we have identified a remote relationship between Bsp6I and the structure of PvuII. Starting from the sequence–structure alignment between Bsp6I and PvuII, we constructed a homology model of Bsp6I and used it to predict functionally significant regions in Bsp6I. The homology model was supported by site-directed mutagenesis of residues predicted to be important for dimerization, DNA binding and catalysis. Completing the picture of sequence–structure–function relationships in protein superfamilies becomes an essential task in the age of structural genomics and our study may serve as a paradigm for future analyses of superfamilies comprising strongly diverged members with little or no sequence similarity.