Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily,but exhibits an unusual active site

Background  

Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds.     

Oligomeric Structure Diversity within the GIY-YIG Nuclease Family

The GIY-YIG nuclease domain has been identified in homing endonucleases, DNA repair and recombination enzymes, and restriction endonucleases. The Type II restriction enzyme Eco29kI belongs to the GIY-YIG nuclease superfamily and, like most of other family members, including the homing endonuclease I-TevI, is a monomer. It recognizes the palindromic sequence 5′-CCGC/GG-3′ (“/” marks the cleavage position) and cuts it to generate 3′-staggered ends. The Eco29kI monomer, which contains a single active site, either has to nick sequentially individual DNA strands or has to form dimers or even higher-order oligomers upon DNA binding to make a double-strand break at its target site. Here, we provide experimental evidence that Eco29kI monomers dimerize on a single cognate DNA molecule forming the catalytically active complex. The mechanism described here for Eco29kI differs from that of Cfr42I isoschisomer, which also belongs to the GIY-YIG family but is functional as a tetramer. This novel mechanism may have implications for the function of homing endonucleases and other enzymes of the GIY-YIG family.     

Amino acid residues in the GIY-YIG endonuclease II of phage T4 affecting sequence recognition and binding as well as catalysis

Phage T4 endonuclease II (EndoII), a GIY-YIG endonuclease lacking a carboxy-terminal DNA-binding domain, was subjected to site-directed mutagenesis to investigate roles of individual amino acids in substrate recognition, binding, and catalysis. The structure of EndoII was modeled on that of UvrC. We found catalytic roles for residues in the putative catalytic surface (G49, R57, E118, and N130) similar to those described for I-TevI and UvrC; in addition, these residues were found to be important for substrate recognition and binding. The conserved glycine (G49) and arginine (R57) were essential for normal sequence recognition. Our results are in agreement with a role for these residues in forming the DNA-binding surface and exposing the substrate scissile bond at the active site. The conserved asparagine (N130) and an adjacent proline (P127) likely contribute to positioning the catalytic domain correctly. Enzymes in the EndoII subfamily of GIY-YIG endonucleases share a strongly conserved middle region (MR, residues 72 to 93, likely helical and possibly substituting for heterologous helices in I-TevI and UvrC) and a less strongly conserved N-terminal region (residues 12 to 24). Most of the conserved residues in these two regions appeared to contribute to binding strength without affecting the mode of substrate binding at the catalytic surface. EndoII K76, part of a conserved NUMOD3 DNA-binding motif of homing endonucleases found to overlap the MR, affected both sequence recognition and catalysis, suggesting a more direct involvement in positioning the substrate. Our data thus suggest roles for the MR and residues conserved in GIY-YIG enzymes in recognizing and binding the substrate.     

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.     

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.     

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.     

Catalytic domain structure and hypothesis for function of GIY-YIG intron endonuclease I-TevI

I-TevI, a member of the GIY-YIG family of homing endonucleases, consists of an N-terminal catalytic domain and a C-terminal DNA-binding domain joined by a flexible linker. The GIY-YIG motif is in the N-terminal domain of I-TevI, which corresponds to a phylogenetically widespread catalytic cartridge that is often associated with mobile genetic elements. The crystal structure of the catalytic domain of I-TevI, the first of any GIY-YIG endonuclease, reveals a novel alpha/beta-fold with a central three-stranded antiparallel beta-sheet flanked by three helices. The most conserved and putative catalytic residues are located on a shallow, concave surface and include a metal coordination site. Similarities in the three-dimensional arrangement of the catalytically important residues and the cation-binding site with those of the His-Cys box endonuclease I-PpoI suggest the possibility of mechanistic relationships among these different families of homing endonucleases despite completely different folds.     

Nucleotide flipping by restriction enzymes analyzed by 2-aminopurine steady-state fluorescence

Many DNA modification and repair enzymes require access to DNA bases and therefore flip nucleotides. Restriction endonucleases (REases) hydrolyze the phosphodiester backbone within or in the vicinity of the target recognition site and do not require base extrusion for the sequence readout and catalysis. Therefore, the observation of extrahelical nucleotides in a co-crystal of REase Ecl18kI with the cognate sequence, CCNGG, was unexpected. It turned out that Ecl18kI reads directly only the CCGG sequence and skips the unspecified N nucleotides, flipping them out from the helix. Sequence and structure conservation predict nucleotide flipping also for the complexes of PspGI and EcoRII with their target DNAs (/CCWGG), but data in solution are limited and indirect. Here, we demonstrate that Ecl18kI, the C-terminal domain of EcoRII (EcoRII-C) and PspGI enhance the fluorescence of 2-aminopurines (2-AP) placed at the centers of their recognition sequences. The fluorescence increase is largest for PspGI, intermediate for EcoRII-C and smallest for Ecl18kI, probably reflecting the differences in the hydrophobicity of the binding pockets within the protein. Omitting divalent metal cations and mutation of the binding pocket tryptophan to alanine strongly increase the 2-AP signal in the Ecl18kI–DNA complex. Together, our data provide the first direct evidence that Ecl18kI, EcoRII-C and PspGI flip nucleotides in solution.     

Configuration of the catalytic GIY-YIG domain of intron endonuclease I-TevI: coincidence of computational and molecular findings

I-TevI is a member of the GIY-YIG family of homing endonucleases. It is folded into two structural and functional domains, an N-terminal catalytic domain and a C-terminal DNA-binding domain, separated by a flexible linker. In this study we have used genetic analyses, computational sequence analysis andNMR spectroscopy to define the configuration of theN-terminal domain and its relationship to the flexible linker. The catalytic domain is an alpha/beta structure contained within the first 92 amino acids of the 245-amino acid protein followed by an unstructured linker. Remarkably, this structured domain corresponds precisely to the GIY-YIG module defined by sequence comparisons of 57 proteins including more than 30 newly reported members of the family. Although much of the unstructured linker is not essential for activity, residues 93-116 are required, raising the possibility that this region may adopt an alternate conformation upon DNA binding. Two invariant residues of the GIY-YIG module, Arg27 and Glu75, located in alpha-helices, have properties of catalytic residues. Furthermore, the GIY-YIG sequence elements for which the module is named form part of a three-stranded antiparallel beta-sheet that is important for I-TevI structure and function.     

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.     

Three-dimensional modeling of the I-TevI homing endonuclease catalytic domain, a GIY-YIG superfamily member, using NMR restraints and Monte Carlo dynamics.

Using a recent version of the SICHO algorithm for in silico protein folding, we made a blind prediction of the tertiary structure of the N-terminal, independently folded, catalytic domain (CD) of the I-TevI homing endonuclease, a representative of the GIY-YIG superfamily of homing endonucleases. The secondary structure of the I-TevI CD has been determined using NMR spectroscopy, but computational sequence analysis failed to detect any protein of known tertiary structure related to the GIY-YIG nucleases (Kowalski et al., Nucleic Acids Res., 1999, 27, 2115-2125). To provide further insight into the structure-function relationships of all GIY-YIG superfamily members, including the recently described subfamily of type II restriction enzymes (Bujnicki et al., Trends Biochem. Sci., 2000, 26, 9-11), we incorporated the experimentally determined and predicted secondary and tertiary restraints in a reduced (side chain only) protein model, which was minimized by Monte Carlo dynamics and simulated annealing. The subsequently elaborated full atomic model of the I-TevI CD allows the available experimental data to be put into a structural context and suggests that the GIY-YIG domain may dimerize in order to bring together the conserved residues of the active site.     

Theoretical model of restriction endonuclease HpaI in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis

Type II restriction enzymes are commercially important deoxyribonucleases and very attractive targets for protein engineering of new specificities. At the same time they are a very challenging test bed for protein structure prediction methods. Typically, enzymes that recognize different sequences show little or no amino acid sequence similarity to each other and to other proteins. Based on crystallographic analyses that revealed the same PD-(D/E)XK fold for more than a dozen case studies, they were nevertheless considered to be related until the combination of bioinformatics and mutational analyses has demonstrated that some of these proteins belong to other, unrelated folds PLD, HNH, and GIY-YIG. As a part of a large-scale project aiming at identification of a three-dimensional fold for all type II REases with known sequences (currently approximately 1000 proteins), we carried out preliminary structure prediction and selected candidates for experimental validation. Here, we present the analysis of HpaI REase, an ORFan with no detectable homologs, for which we detected a structural template by protein fold recognition, constructed a model using the FRankenstein monster approach and identified a number of residues important for the DNA binding and catalysis. These predictions were confirmed by site-directed mutagenesis and in vitro analysis of the mutant proteins. The experimentally validated model of HpaI will serve as a low-resolution structural platform for evolutionary considerations in the subgroup of blunt-cutting REases with different specificities. The research protocol developed in the course of this work represents a streamlined version of the previously used techniques and can be used in a high-throughput fashion to build and validate models for other enzymes, especially ORFans that exhibit no sequence similarity to any other protein in the database.     

Construction of an overproducing strain,purification, and biochemical characterization of the 6His-Eco29kI restriction endonuclease

We constructed a strain of Escherichia coli overproducing 6His-tagged Eco29kI by placing the coding sequence under control of a strong bacteriophage T5 promoter. The yield of 6His-Eco29kI restriction endonuclease expression could be increased to about 20% of the total cellular protein, but inclusion bodies formed consisting of insoluble 6His-Eco29kI protein. We developed a fast and effective protocol for purification of the homogeneous enzyme from both soluble and insoluble fractions and established their identity by catalytic activity assay. The isolated enzymes were tested for recognition specificity and optimal reaction conditions as a function of NaCl and KCl concentrations, temperature, and pH compared with the native Eco29kI restriction endonuclease. The 6His-tagged enzyme retained the specificity of the native protein but had an altered optimum of its catalytic reaction.     

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.     

Coincidence of cleavage sites of intron endonuclease I-TevI and critical sequences of the host thymidylate synthase gene

To maximize spread of their host intron or intein, many homing endonucleases recognize nucleotides that code for important and conserved amino acid residues of the target gene. Here, we examine the cleavage requirements for I-TevI, which binds a stretch of thymidylate synthase (TS) DNA that codes for functionally critical residues in the TS active site. Using an in vitro selection scheme, we identified two base-pairs in the I-TevI cleavage site region as important for cleavage efficiency. These were confirmed by comparison of I-TevI cleavage efficiencies on mutant and on wild-type substrates. We also showed that nicking of the bottom strand by I-TevI is not affected by mutation of residues surrounding the bottom-strand cleavage site, unlike other homing endonucleases. One of these two base-pairs is universally conserved in all TS sequences, and is identical with a previously identified cleavage determinant of I-BmoI, a related GIY-YIG endonuclease that binds a homologous stretch of TS-encoding DNA. The other base-pair is conserved only in a subset of TS genes that includes the I-TevI, but not the I-BmoI, target sequence. Both the I-TevI and I-BmoI cleavage site requirements correspond to functionally critical residues involved in an extensive hydrogen bond network within the TS active site. Remarkably, these cleavage requirements correlate with TS phylogeny in bacteria, suggesting that each endonuclease has individually adapted to efficiently cleave distinct TS substrates.     

Effects of Specific Amino Acid Substitutions on Activities of Dinitrogenase Reductase-Activating Glycohydrolase from Rhodospirillum rubrum

Site-directed mutagenesis of the draG gene was used to generate altered forms of dinitrogenase reductase-activating glycohydrolase (DRAG) with D123A, H142L, H158N, D243G, and E279R substitutions. The amino acid residues H142 and E279 are not required either for the coordination to the metal center or for catalysis since the variants H142L and E279R retained both catalytic and electron paramagnetic resonance spectral properties similar to those of the wild-type enzyme. Since DRAG-H158N and DRAG-D243G variants lost their ability to bind Mn(II) and to catalyze the hydrolysis of the substrate, H158 and D243 residues could be involved in the coordination of the binuclear Mn(II) center in DRAG.     

Related bifunctional restriction endonuclease-methyltransferase triplets: TspDTI, Tth111II/TthHB27I and TsoI with distinct specificities

Background

We previously defined a family of restriction endonucleases (REases) from Thermus sp., which share common biochemical and biophysical features, such as the fusion of both the nuclease and methyltransferase (MTase) activities in a single polypeptide, cleavage at a distance from the recognition site, large molecular size, modulation of activity by S-adenosylmethionine (SAM), and incomplete cleavage of the substrate DNA. Members include related thermophilic REases with five distinct specificities: TspGWI, TaqII, Tth111II/TthHB27I, TspDTI and TsoI.

Results

TspDTI, TsoI and isoschizomers Tth111II/TthHB27I recognize different, but related sequences: 5'-ATGAA-3', 5'-TARCCA-3' and 5'-CAARCA-3' respectively. Their amino acid sequences are similar, which is unusual among REases of different specificity. To gain insight into this group of REases, TspDTI, the prototype member of the Thermus sp. enzyme family, was cloned and characterized using a recently developed method for partially cleaving REases.

Conclusions

TspDTI, TsoI and isoschizomers Tth111II/TthHB27I are closely related bifunctional enzymes. They comprise a tandem arrangement of Type I-like domains, like other Type IIC enzymes (those with a fusion of a REase and MTase domains), e.g. TspGWI, TaqII and MmeI, but their sequences are only remotely similar to these previously characterized enzymes. The characterization of TspDTI, a prototype member of this group, extends our understanding of sequence-function relationships among multifunctional restriction-modification enzymes.     

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.     

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.