Perpetuating the homing endonuclease life cycle: identification of mutations that modulate and change I-TevI cleavage preference

Homing endonucleases are sequence-tolerant DNA endonucleases that act as mobile genetic elements. The ability of homing endonucleases to cleave substrates with multiple nucleotide substitutions suggests a high degree of adaptability in that changing or modulating cleavage preference would require relatively few amino acid substitutions. Here, using directed evolution experiments with the GIY-YIG homing endonuclease I-TevI that targets the thymidylate synthase gene of phage T4, we readily isolated variants that dramatically broadened I-TevI cleavage preference, as well as variants that fine-tuned cleavage preference. By combining substitutions, we observed an ∼10 000-fold improvement in cleavage on some substrates not cleaved by the wild-type enzyme, correlating with a decrease in readout of information content at the cleavage site. Strikingly, we were able to change the cleavage preference of I-TevI to that of the isoschizomer I-BmoI which targets a different cleavage site in the thymidylate synthase gene, recapitulating the evolution of cleavage preference in this family of homing endonucleases. Our results define a strategy to isolate GIY-YIG nuclease domains with distinct cleavage preferences, and provide insight into how homing endonucleases may escape a dead-end life cycle in a population of saturated target sites by promoting transposition to different target sites.     

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.     

Importance of a single base pair for discrimination between intron-containing and intronless alleles by endonuclease I-BmoI

Homing endonucleases initiate mobility of their host group I introns by binding to and cleaving lengthy recognition sequences that are typically centered on the intron insertion site (IS) of intronless alleles. Because the intron interrupts the endonucleases' recognition sequence, intron-containing alleles are immune to cleavage by their own endonuclease. I-TevI and I-BmoI are related GIY-YIG endonucleases that bind a homologous stretch of thymidylate synthase (TS)-encoding DNA but use different strategies to distinguish intronless from intron-containing substrates. I-TevI discriminates between substrates at the level of DNA binding, as its recognition sequence is centered on the intron IS. I-BmoI, in contrast, possesses a very asymmetric recognition sequence with respect to the intron IS, binds both intron-containing and intronless TS-encoding substrates, but efficiently cleaves only intronless substrate. Here, we show that I-BmoI is extremely tolerant of multiple substitutions around its cleavage sites and has a low specific activity. However, a single G-C base pair, at position -2 of a 39-base pair recognition sequence, is a major determinant for cleavage efficiency and distinguishes intronless from intron-containing alleles. Strikingly, this G-C base pair is universally conserved in phylogenetically diverse TS-coding sequences; this finding suggests that I-BmoI has evolved exquisite cleavage requirements to maximize the potential to spread to variant intronless alleles, while minimizing cleavage at its own intron-containing allele.     

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.     

SegH and Hef: two novel homing endonucleases whose genes replace the mobC and mobE genes in several T4-related phages

T4 contains two groups of genes with similarity to homing endonucleases, the seg-genes (similarity to endonucleases encoded by group I introns) containing GIY-YIG motifs and the mob-genes (similarity to mobile endonucleases) containing H-N-H motifs. The four seg-genes characterized to date encode homing endonucleases with cleavage sites close to their respective gene loci while none of the mob-genes have been shown to cleave DNA. Of 18 phages screened, only T4 was found to have mobC while mobE genes were found in five additional phages. Interestingly, three phages encoded a seg-like gene (hereby called segH) with a GIY-YIG motif in place of mobC. An additional phage has an unrelated gene called hef (homing endonuclease-like function) in place of the mobE gene. The gene products of both novel genes displayed homing endonuclease activity with cleavage site specificity close to their respective genes. In contrast to intron encoded homing endonucleases, both SegH and Hef can cleave their own DNA as well as DNA from phages without the genes. Both segH and mobE (and most likely hef) can home between phages in mixed infections. We discuss why it might be a selective advantage for phage freestanding homing endonucleases to cleave both HEG-containing and HEG-less genomes.     

Strand-specific contacts and divalent metal ion regulate double-strand break formation by the GIY-YIG homing endonuclease I-BmoI

GIY-YIG homing endonucleases are modular enzymes consisting of a well-defined N-terminal catalytic domain connected to a variable C-terminal DNA-binding domain. Previous studies have revealed that the role of the DNA-binding domain is to recognize and bind intronless DNA substrate, positioning the N-terminal catalytic domain such that it is poised to generate a staggered double-strand break by an unknown mechanism. Interactions of the N-terminal catalytic domain with intronless substrate are therefore a critical step in the reaction pathway but have been difficult to define. Here, we have taken advantage of the reduced activity of I-BmoI, an isoschizomer of the well-studied bacteriophage T4 homing endonuclease I-TevI, to examine double-strand break formation by I-BmoI. We present evidence demonstrating that I-BmoI generates a double-strand break by two sequential but chemically independent nicking reactions where divalent metal ion is a limiting factor in top-strand nicking. We also show by in-gel footprinting that contacts by the I-BmoI catalytic domain induce significant minor groove DNA distortions that occur independently of bottom-strand nicking. Bottom-strand contacts are critical for accurate top-strand nicking, whereas top-strand contacts have little influence on the accuracy of bottom-strand nicking. We discuss our results in the context of current models of GIY-YIG endonuclease function, with emphasis on the role of divalent metal ion and strand-specific contacts in regulating the activity of a single active site to generate a staggered double-strand break.     

The monomeric GIY-YIG homing endonuclease I-BmoI uses a molecular anchor and a flexible tether to sequentially nick DNA

The GIY-YIG nuclease domain is found within protein scaffolds that participate in diverse cellular pathways and contains a single active site that hydrolyzes DNA by a one-metal ion mechanism. GIY-YIG homing endonucleases (GIY-HEs) are two-domain proteins with N-terminal GIY-YIG nuclease domains connected to C-terminal DNA-binding and they are thought to function as monomers. Using I-BmoI as a model GIY-HE, we test mechanisms by which the single active site is used to generate a double-strand break. We show that I-BmoI is partially disordered in the absence of substrate, and that the GIY-YIG domain alone has weak affinity for DNA. Significantly, we show that I-BmoI functions as a monomer at all steps of the reaction pathway and does not transiently dimerize or use sequential transesterification reactions to cleave substrate. Our results are consistent with the I-BmoI DNA-binding domain acting as a molecular anchor to tether the GIY-YIG domain to substrate, permitting rotation of the GIY-YIG domain to sequentially nick each DNA strand. These data highlight the mechanistic differences between monomeric GIY-HEs and dimeric or tetrameric GIY-YIG restriction enzymes, and they have implications for the use of the GIY-YIG domain in genome-editing applications.     

Divalent metal ion differentially regulates the sequential nicking reactions of the GIY-YIG homing endonuclease I-BmoI

Homing endonucleases are site-specific DNA endonucleases that function as mobile genetic elements by introducing double-strand breaks or nicks at defined locations. Of the major families of homing endonucleases, the modular GIY-YIG endonucleases are least understood in terms of mechanism. The GIY-YIG homing endonuclease I-BmoI generates a double-strand break by sequential nicking reactions during which the single active site of the GIY-YIG nuclease domain must undergo a substantial reorganization. Here, we show that divalent metal ion plays a significant role in regulating the two independent nicking reactions by I-BmoI. Rate constant determination for each nicking reaction revealed that limiting divalent metal ion has a greater impact on the second strand than the first strand nicking reaction. We also show that substrate mutations within the I-BmoI cleavage site can modulate the first strand nicking reaction over a 314-fold range. Additionally, in-gel DNA footprinting with mutant substrates and modeling of an I-BmoI-substrate complex suggest that amino acid contacts to a critical GC-2 base pair are required to induce a bottom-strand distortion that likely directs conformational changes for reaction progress. Collectively, our data implies mechanistic roles for divalent metal ion and substrate bases, suggesting that divalent metal ion facilitates the re-positioning of the GIY-YIG nuclease domain between sequential nicking reactions.     

102 Genome engineering nucleases derived from GIY-YIG homing endonucleases

Efficient targeted manipulation of complex genomes requires highly specific endonucleases to generate double-strand breaks at defined locations (Bibikova et al., 2003; Bogdanove and Voytas, 2011). The predominantly engineered nucleases, zinc-finger nucleases (ZFNs), and TAL effector nucleases (TALENs) use the catalytic domain of FokI as the nuclease portion. This domain, however, functions as a dimer to nonspecifically cleave DNA meaning that ZFNs and TALENs must be designed in head-to-head pairs to target a desired sequence. To overcome this limitation and expand the toolbox of genome editing reagents, we used the N-terminal catalytic domain and interdomain linker of the monomeric GIY-YIG homing endonuclease I-TevI to create I-TevI-zinc-fingers (Tev-ZFEs), and I-TevI-TAL effectors (Tev-TALs) (Kleinstiver et al. 2012). We also made I-TevI fusions to LAGLIDADGs homing endonucleases (I-Tev-LHEs). All the three fusions showed activity on model substrates on par with ZFNs and TALENs in yeast-based recombination assays. These proof-of-concept experiments demonstrate that the catalytic domain of GIY-YIG homing endonucleases can be targeted to relevant loci by fusing the domain to characterize DNA-binding platforms. Recent efforts have focused on improving the Tev-TAL platform by (1) understanding the spacing requirements between the nuclease cleavage site and the DNA binding site, (2) probing the DNA binding requirements of the I-TevI linker domain, and (3) demonstrating activity in mammalian systems.     

SegG endonuclease promotes marker exclusion and mediates co-conversion from a distant cleavage site

Bacteriophages T2 and T4 are closely related T-even phages. However, T4 genetic markers predominate in the progeny of mixed infections, a phenomenon termed marker exclusion. One region previously mapped where the frequency of T2 markers in the progeny is extremely low is located around gene 32. Here, we describe SegG, a GIY-YIG family endonuclease adjacent to gene 32 of phage T4 that is absent from phage T2. In co-infections with T2 and T4, cleavage in T2 gene 32 by T4-encoded SegG initiates a gene conversion event that results in replacement of T2 gene 32 markers with the corresponding T4 sequence. Interestingly, segG inheritance is limited, apparently because of the physical separation of its cleavage and insertion sites, which are 332 base-pairs apart. This contrasts with efficient inheritance of the phage T4 td group I intron and its endonuclease, I-TevI, for which the distance separating the I-TevI cleavage site and td insertion site is 23 base-pairs. Furthermore, we show that co-conversion tracts generated by repair of SegG and I-TevI double-strand breaks contribute to the localized exclusion of T2 markers. Our results demonstrate that the endonuclease activities of SegG and I-TevI promote the spread of these two endonucleases to progeny phage, consistent with their role as selfish genetic elements, and also provide a mechanism by which the genetic contribution of T2 markers to progeny phage is reduced.     

The nicking homing endonuclease I-BasI is encoded by a group I intron in the DNA polymerase gene of the Bacillus thuringiensis phage Bastille

Here we describe the discovery of a group I intron in the DNA polymerase gene of Bacillus thuringiensis phage Bastille. Although the intron insertion site is identical to that of the Bacillus subtilis phages SPO1 and SP82 introns, the Bastille intron differs from them substantially in primary and secondary structure. Like the SPO1 and SP82 introns, the Bastille intron encodes a nicking DNA endonuclease of the H-N-H family, I-BasI, with a cleavage site identical to that of the SPO1-encoded enzyme I-HmuI. Unlike I-HmuI, which nicks both intron-minus and intron-plus DNA, I-BasI cleaves only intron-minus alleles, which is a characteristic of typical homing endonucleases. Interestingly, the C-terminal portions of these H-N-H phage endonucleases contain a conserved sequence motif, the intron-encoded endonuclease repeat motif (IENR1) that also has been found in endonucleases of the GIY-YIG family, and which likely comprises a small DNA-binding module with a globular ββααβ fold, suggestive of module shuffling between different homing endonuclease families.     

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.     

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.     

Structural insights into the first incision reaction during nucleotide excision repair

Nucleotide excision repair is a highly conserved DNA repair mechanism present in all kingdoms of life. The incision reaction is a critical step for damage removal and is accomplished by the UvrC protein in eubacteria. No structural information is so far available for the 3' incision reaction. Here we report the crystal structure of the N-terminal catalytic domain of UvrC at 1.5 A resolution, which catalyzes the 3' incision reaction and shares homology with the catalytic domain of the GIY-YIG family of intron-encoded homing endonucleases. The structure reveals a patch of highly conserved residues surrounding a catalytic magnesium-water cluster, suggesting that the metal binding site is an essential feature of UvrC and all GIY-YIG endonuclease domains. Structural and biochemical data strongly suggest that the N-terminal endonuclease domain of UvrC utilizes a novel one-metal mechanism to cleave the phosphodiester bond.     

Role of the interdomain linker in distance determination for remote cleavage by homing endonuclease I-TevI

I-TevI is a modular intron-encoded endonuclease, consisting of an N-terminal catalytic domain and a C-terminal DNA-binding domain, joined by a 75 amino acid linker. This linker can be divided into three regions, starting at the N terminus: the deletion-intolerant (DI) region; the deletion-tolerant (DT) region; and a zinc finger, which acts as a distance determinant for cleavage. To further explore linker function, we generated deletion and substitution mutants that were tested for their preference to cleave at a particular distance or at the correct sequence. Our results demonstrate that the I-TevI linker is multi-functional, a property that sets it apart from junction sequences in most other proteins. First, the linker DI region has a role in I-TevI cleavage activity. Second, the DT linker region participates in distance determination, as evident from DT mutants that display a phenotype similar to that of the zinc-finger mutants in their selection of a cleavage site. Finally, NMR analysis of a freestanding 56 residue linker segment showed an unstructured stretch corresponding to the DI region and a portion of the DT region, followed by a β-strand corresponding to the remainder of the DT region and containing a key distance-determining arginine, R129. Mutation of this arginine to alanine abolished distance determination and disrupted the β-strand, indicating that the structure of the DT linker region has a role in cleavage at a fixed distance.     

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.     

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.