Profile of the DNA recognition site of the archaeal homing endonuclease I-DmoI

F-CphI represents a new homing endonuclease family using the Endo VII catalytic motif

Background

There are six known families of homing endonucleases, LAGLIDADG, GIY-YIG, HNH, His-Cys box, PD-(D/E)-XK, and EDxHD, which are characterized by their conserved residues. Previously, we discovered a novel homing endonuclease F-CphI encoded by ORF177 of cyanophage S-PM2. F-CphI does not resemble any characterized homing endonucleases. Instead, the C-terminus of F-CphI aligns well with the N-terminal catalytic domain of a Holliday junction DNA resolvase, phage T4 endonuclease VII (Endo VII).

Results

A PSI-BLAST search resulted in a total of 313 Endo VII motif–containing sequences in sequenced genomes. Multiple sequence alignment showed that the catalytically important residues of T4 Endo VII were all well conserved in these proteins. Our site-directed mutagenesis studies further confirmed that the catalytically important residues of T4 Endo VII were also essential for F-CphI activity, and thus F-CphI might use a similar protein fold as Endo VII for DNA cleavage. A phylogenetic tree of the Endo VII motif–containing sequences showed that putative resolvases grouped into one clade while putative homing endonucleases and restriction endonucleases grouped into another clade.

Conclusions

Based on the unique conserved residues, we proposed that F-CphI represents a new homing endonuclease family, which was named the DHHRN family. Our phylogenetic analysis could be used to predict the functions of many previously unknown proteins.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0132-5) contains supplementary material, which is available to authorized users.

A group II intron encodes a functional LAGLIDADG homing endonuclease and self-splices under moderate temperature and ionic conditions

A group II intron encoding a protein belonging to the LAGLIDADG family of homing endonucleases was identified in the mitochondrial rns gene of the filamentous fungus Leptographium truncatum, and the catalytic activities of both the intron and its encoded protein were characterized. A model of the RNA secondary structure indicates that the intron is a member of the IIB1 subclass and the open reading frame is inserted in ribozyme domain III. In vitro assays carried out with two versions of the intron, one in which the open reading frame was removed and the other in which it was present, demonstrate that both versions of the intron readily self-splice at 37 degrees C and at a concentration of MgCl(2) as low as 6 mM. The open reading frame encodes a functional LAGLIDADG homing endonuclease that cleaves 2 (top strand) and 6 (bottom strand) nucleotides (nt) upstream of the intron insertion site, generating 4 nt 3' OH overhangs. In vitro splicing assays carried out in the absence and presence of the intron-encoded protein indicate that the protein does not enhance intron splicing, and RNA-binding assays show that the protein does not appear to bind to the intron RNA precursor transcript. These findings raise intriguing questions concerning the functional and evolutionary relationships of the two components of this unique composite element.

Homing endonucleases: from basics to therapeutic applications

Homing endonucleases (HE) are double-stranded DNAses that target large recognition sites (12-40 bp). HE-encoding sequences are usually embedded in either introns or inteins. Their recognition sites are extremely rare, with none or only a few of these sites present in a mammalian-sized genome. However, these enzymes, unlike standard restriction endonucleases, tolerate some sequence degeneracy within their recognition sequence. Several members of this enzyme family have been used as templates to engineer tools to cleave DNA sequences that differ from their original wild-type targets. These custom HEs can be used to stimulate double-strand break homologous recombination in cells, to induce the repair of defective genes with very low toxicity levels. The use of tailored HEs opens up new possibilities for gene therapy in patients with monogenic diseases that can be treated ex vivo. This review provides an overview of recent advances in this field.

Phylogenetic screening of a bacterial, metagenomic library using homing endonuclease restriction and marker insertion

Metagenomics provides access to the uncultured majority of the microbial world. The approaches employed in this field have, however, had limited success in linking functional genes to the taxonomic or phylogenetic origin of the organism they belong to. Here we present an efficient strategy to recover environmental DNA fragments that contain phylogenetic marker genes from metagenomic libraries. Our method involves the cleavage of 23S ribsosmal RNA (rRNA) genes within pooled library clones by the homing endonuclease I-CeuI followed by the insertion and selection of an antibiotic resistance cassette. This approach was applied to screen a library of 6500 fosmid clones derived from the microbial community associated with the sponge Cymbastela concentrica. Several fosmid clones were recovered after the screen and detailed phylogenetic and taxonomic assignment based on the rRNA gene showed that they belong to previously unknown organisms. In addition, compositional features of these fosmid clones were used to classify and taxonomically assign a dataset of environmental shotgun sequences. Our approach represents a valuable tool for the analysis of rapidly increasing, environmental DNA sequencing information.

High-resolution profiling of homing endonuclease binding and catalytic specificity using yeast surface display

Experimental analysis and manipulation of protein–DNA interactions pose unique biophysical challenges arising from the structural and chemical homogeneity of DNA polymers. We report the use of yeast surface display for analytical and selection-based applications for the interaction between a LAGLIDADG homing endonuclease and its DNA target. Quantitative flow cytometry using oligonucleotide substrates facilitated a complete profiling of specificity, both for DNA-binding and catalysis, with single base pair resolution. These analyses revealed a comprehensive segregation of binding specificity and affinity to one half of the pseudo-dimeric interaction, while the entire interface contributed specificity at the level of catalysis. A single round of targeted mutagenesis with tandem affinity and catalytic selection steps provided mechanistic insights to the origins of binding and catalytic specificity. These methods represent a dynamic new approach for interrogating specificity in protein–DNA interactions.

Generation of a nicking enzyme that stimulates site-specific gene conversion from the I-AniI LAGLIDADG homing endonuclease

Homing endonucleases stimulate gene conversion by generating site-specific DNA double-strand breaks that are repaired by homologous recombination. These enzymes are potentially valuable tools for targeted gene correction and genome engineering. We have engineered a variant of the I-AniI homing endonuclease that nicks its cognate target site. This variant contains a mutation of a basic residue essential for proton transfer and solvent activation in one active site. The cleavage mechanism, DNA-binding affinity, and substrate specificity profile of the nickase are similar to the wild-type enzyme. I-AniI nickase stimulates targeted gene correction in human cells, in cis and in trans, at approximately 1/4 the efficiency of the wild-type enzyme. The development of sequence-specific nicking enzymes like the I-AniI nickase will facilitate comparative analyses of DNA repair and mutagenesis induced by single- or double-strand breaks.

Optimization of in vivo activity of a bifunctional homing endonuclease and maturase reverses evolutionary degradation

The LAGLIDADG homing endonuclease (LHE) I-AniI has adopted an extremely efficient secondary RNA splicing activity that is beneficial to its host, balanced against inefficient DNA cleavage. A selection experiment identified point mutations in the enzyme that act synergistically to improve endonuclease activity. The amino-acid substitutions increase target affinity, alter the thermal cleavage profile and significantly increase targeted recombination in transfected cells. The RNA splicing activity is not affected by these mutations. The improvement in DNA cleavage activity is largely focused on one of the enzyme's two active sites, corresponding to a rearrangement of a lysine residue hypothesized to act as a general base. Most of the constructs isolated in the screen contain one or more mutations that revert an amino-acid identity to a residue found in one or more close homologues of I-AniI. This implies that mutations that have previously reduced the endonuclease activity of I-AniI are identified and reversed, sometimes in combination with additional ‘artificial’ mutations, to optimize its in vivo activity.

Crystal structure of I-DmoI in complex with its target DNA provides new insights into meganuclease engineering

Homing endonucleases, also known as meganucleases, are sequence-specific enzymes with large DNA recognition sites. These enzymes can be used to induce efficient homologous gene targeting in cells and plants, opening perspectives for genome engineering with applications in a wide series of fields, ranging from biotechnology to gene therapy. Here, we report the crystal structures at 2.0 and 2.1 A resolution of the I-DmoI meganuclease in complex with its substrate DNA before and after cleavage, providing snapshots of the catalytic process. Our study suggests that I-DmoI requires only 2 cations instead of 3 for DNA cleavage. The structure sheds light onto the basis of DNA binding, indicating key residues responsible for nonpalindromic target DNA recognition. In silico and in vivo analysis of the I-DmoI DNA cleavage specificity suggests that despite the relatively few protein-base contacts, I-DmoI is highly specific when compared with other meganucleases. Our data open the door toward the generation of custom endonucleases for targeted genome engineering using the monomeric I-DmoI scaffold.

Generation and analysis of mesophilic variants of the thermostable archaeal I-DmoI homing endonuclease

The hyperthermophilic archaeon Desulfurococcus mobilis I-DmoI protein belongs to the family of proteins known as homing endonucleases (HEs). HEs are highly specific DNA-cleaving enzymes that recognize long stretches of DNA and are powerful tools for genome engineering. Because of its monomeric nature, I-DmoI is an ideal scaffold for generating mutant enzymes with novel DNA specificities, similarly reported for homodimeric HEs, but providing single chain endonucleases instead of dimers. However, this would require the use of a mesophilic variant cleaving its substrate at temperatures of 37 degrees C and below. We have generated mesophilic mutants of I-DmoI, using a single round of directed evolution that relies on a functional assay in yeast. The effect of mutations identified in the novel proteins has been investigated. These mutations are located distant to the DNA-binding site and cause changes in the size and polarity of buried residues, suggesting that they act by destabilizing the protein. Two of the novel proteins have been produced and analyzed in vitro. Their overall structures are similar to that of the parent protein, but they are destabilized against thermal and chemical denaturation. The temperature-dependent activity profiles for the mutants shifted toward lower temperatures with respect to the wild-type activity profile. However, the most destabilized mutant was not the most active at low temperatures, suggesting that other effects, like local structural distortions and/or changes in the protein dynamics, also influence their activity. These mesophilic I-DmoI mutants form the basis for generating new variants with tailored DNA specificities.

Coevolution of a homing endonuclease and its host target sequence

We have determined the specificity profile of the homing endonuclease I-AniI and compared it to the conservation of its host gene. Homing endonucleases are encoded within intervening sequences such as group I introns. They initiate the transfer of such elements by cleaving cognate alleles lacking the intron, leading to their transfer via homologous recombination. Each structural homing endonuclease family has arrived at an appropriate balance of specificity and fidelity that avoids toxicity while maximizing target recognition and invasiveness. I-AniI recognizes a strongly conserved target sequence in a host gene encoding apocytochrome B and has fine-tuned its specificity to correlate with wobble versus nonwobble positions across that sequence and to the amount of degeneracy inherent in individual codons. The physiological target site in the host gene is not the optimal substrate for recognition and cleavage: at least one target variant identified during a screen is bound more tightly and cleaved more rapidly. This is a result of the periodic cycle of intron homing, which at any time can present nonoptimal combinations of endonuclease specificity and insertion site sequences in a biological host.

Evolution of Pleopsidium (Lichenized Ascomycota) S943 Group I Introns and the Phylogeography of an Intron-Encoded Putative Homing Endonuclease

The sporadic distribution of nuclear group I introns among different fungal lineages can be explained by vertical inheritance of the introns followed by successive losses, or horizontal transfers from one lineage to another through intron homing or reverse splicing. Homing is mediated by an intron-encoded homing endonuclease (HE) and recent studies suggest that the introns and their associated HE gene (HEG) follow a recurrent cyclical model of invasion, degeneration, loss, and reinvasion. The purpose of this study was to compare this model to the evolution of HEGs found in the group I intron at position S943 of the nuclear ribosomal DNA of the lichen-forming fungus Pleopsidium. Forty-eight S943 introns were found in the 64 Pleopsidium samples from a worldwide screen, 22 of which contained a full-length HEG that encodes a putative 256-amino acid HE, and 2 contained HE pseudogenes. The HEGs are divided into two closely related types (as are the introns that encode them) that differ by 22.6% in their nucleotide sequences. The evolution of the Pleopsidium intron-HEG element shows strong evidence for a cyclical model of evolution. The intron was likely acquired twice in the genus and then transmitted via two or three interspecific horizontal transfers. Close geographical proximity plays an important role in intron-HEG horizontal transfer because most of these mobile elements were found in Europe. Once acquired in a lineage, the intron-HEG element was also vertically transmitted, and occasionally degenerated or was lost.

Degenerated recognition property of a mitochondrial homing enzyme in the unicellular green alga Chlamydomonas smithii

Target sequence cleavage is the essential step for intron invasion into an intronless allele. DNA cleavage at a specific site is performed by an endonuclease, termed a homing enzyme, which is encoded by an open reading frame within the intron. The recognition properties of them have only been analyzed in vitro, using purified, recombinant homing enzyme and various mutated DNA substrates, but it is unclear whether the homing enzyme behaves similarly in vivo. To answer this question, we determined the recognition properties of I-CsmI in vivo. I-CsmI is a homing enzyme encoded by the open reading frame of the alpha-group I-intron, located in the mitochondrial apocytochrome b gene of the green alga Chlamydomonas smithii. The in vivo recognition properties of it were determined as the frequency of intron invasion into a mutated target site. For this purpose, we utilized hybrid diploid cells developed by crossing alpha-intron-plus C. smithii to intron-minus C. reinhardtii containing mutated target sequences. The intron invasion frequency was much higher than the expected from the in vitro cleavage frequency of the respective mutated substrates. Even the substrates that had very little cleavage in the in vitro experiment were efficiently invaded in vivo, and were accompanied by a large degree of coconversion. Considering the ease of the homing enzyme invading into various mutated target sequences, we propose that the principle bottleneck for lateral intron transmission is not the sequence specificity of the homing enzyme, but instead is limited by the rare occurrence of inter-specific cell fusion.

From monomeric to homodimeric endonucleases and back: engineering novel specificity of LAGLIDADG enzymes

Monomeric homing endonucleases of the LAGLIDADG family recognize DNA in a bipartite manner, reflecting the underlying structural assembly of two protein domains (A and B) related by pseudo 2-fold symmetry. This architecture allows for changes in DNA specificity via the distinct combination of these half-site domains. The key to engineering such hybrid proteins lies in the LAGLIDADG two-helix bundle that forms both the domain interface and the endonuclease active site. In this study, we utilize domain A of the monomeric I-DmoI to demonstrate the feasibility of generating functional homodimeric endonucleases that recognize palindromic DNA sequences derived from the original, non-palindromic target. Wild-type I-DmoI domain A is capable of forming a homodimer (H-DmoA) that binds tightly to, but does not cleave efficiently, its anticipated DNA target. Partial restoration of DNA cleavage ability was obtained by re-engineering the LAGLIDADG dimerization interface (H-DmoC). Upon fusing two copies of H-DmoC via a short peptide linker, a novel, site-specific DNA endonuclease was created (H-DmoC2). Like I-DmoI, H-DmoC2 is thermostable and cleaves the new target DNA to generate the predicted 4 nt 3'-OH overhangs but, unlike I-DmoI, H-DmoC2 retains stringent cleavage specificity when substituting Mn2+ for Mg2+ as co-factor. This novel endonuclease allows speculation regarding specificity of monomeric LAGLIDADG proteins, while it supports the evolutionary genesis of these proteins by a gene duplication event.

Lateral gene transfer and phylogenetic assignment of environmental fosmid clones

Metagenomic data, especially sequence data from large insert clones, are most useful when reasonable inferences about phylogenetic origins of inserts can be made. Often, clones that bear phylotypic markers (usually ribosomal RNA genes) are sought, but sometimes phylogenetic assignments have been based on the preponderance of blast hits obtained with predicted protein coding sequences (CDSs). Here we use a cloning method which greatly enriches for ribosomal RNA-bearing fosmid clones to ask two questions: (i) how reliably can we judge the phylogenetic origin of a clone (that is, its RNA phylotype) from the sequences of its CDSs? and (ii) how much lateral gene transfer (LGT) do we see, as assessed by CDSs of different phylogenetic origins on the same fosmid? We sequenced 12 rRNA containing fosmid clones, obtained from libraries constructed using DNA isolated from Baltimore harbour sediments. Three of the clones are from bacterial candidate divisions for which no cultured representatives are available, and thus represent the first protein coding sequences from these major bacterial lineages. The amount of LGT was assessed by making phylogenetic trees of all the CDSs in the fosmid clones and comparing the phylogenetic position of the CDS to the rRNA phylotype. We find that the majority of CDSs in each fosmid, 57-96%, agree with their respective rRNA genes. However, we also find that a significant fraction of the CDSs in each fosmid, 7-44%, has been acquired by LGT. In several cases, we can infer co-transfer of functionally related genes, and generate hypotheses about mechanism and ecological significance of transfer.

Adaptation of intronic homing endonuclease for successful horizontal transmission

Group I introns are thought to be self-propagating mobile elements, and are distributed over a wide range of organisms through horizontal transmission. Intron invasion is initiated through cleavage of a target DNA by a homing endonuclease encoded in an open reading frame (ORF) found within the intron. The intron is likely of no benefit to the host cell and is not maintained over time, leading to the accumulation of mutations after intron invasion. Therefore, regular invasional transmission of the intron to a new species at least once before its degeneration is likely essential for its evolutionary long-term existence. In many cases, the target is in a protein-coding region which is well conserved among organisms, but contains ambiguity at the third nucleotide position of the codon. Consequently, the homing endonuclease might be adapted to overcome sequence polymorphisms at the target site. To address whether codon degeneracy affects horizontal transmission, we investigated the recognition properties of a homing enzyme, I-CsmI, that is encoded in the intronic ORF of a group I intron located in the mitochondrial COB gene of the unicellular green alga Chlamydomonas smithii. We successfully expressed and purified three types of N-terminally truncated I-CsmI polypeptides, and assayed the efficiency of cleavage for 81 substrates containing single nucleotide substitutions. We found a slight but significant tendency that I-CsmI cleaves substrates containing a silent or tolerated amino acid change more efficiently than nonsilent or nontolerated ones. The published recognition properties of I-SpomI, I-ScaI, and I-SceII were reconsidered from this point of view, and we detected proficient adaptation of I-SpomI, I-ScaI, and I-SceII for target site sequence degeneracy. Based on the results described above, we propose that intronic homing enzymes are adapted to cleave sequences that might appear at the target region in various species, however, such adaptation becomes less prominent in proportion to the time elapsed after intron invasion into a new host.

Analysis of the LAGLIDADG interface of the monomeric homing endonuclease I-DmoI

The general structural fold of the LAGLIDADG endonuclease family consists of two similar alpha/beta domains (alphabetabetaalphabetabetaalpha) that assemble either as homodimers or monomers with the domains related by pseudo-two-fold symmetry. At the center of this symmetry is the closely packed LAGLIDADG two-helix bundle that forms the main inter- or intra-molecular contact region between the domains of single- or double-motif proteins, respectively. In this work, we further examine the role of the LAGLIDADG residues involved in the helix-helix interaction. The interchangeability of the LAGLIDADG helix interaction was explored by grafting interfacial residues from the homodimeric I-CreI into the corresponding positions in the monomeric I-DmoI. The resulting LAGLIDADG exchange mutant is partially active, preferring to nick dsDNA rather than making the customary double-strand break. A series of partial revertants within the mutated LAGLIDADG region are shown to restore cleavage activity to varying degrees resulting in one I-DmoI mutant that is more active than wild-type I-DmoI. The phenotype of some of these mutants was reconciled on the basis of similarity to the GxxxG helix interaction found in transmembrane proteins. Additionally, a split variant of I-DmoI was created, demonstrating that the LAGLIDADG helices of I-DmoI are capable of forming and maintaining the protein-protein interface in trans to create an active heterodimer.

Degeneration of a homing endonuclease and its target sequence in a wild yeast strain

Mobile introns and inteins self-propagate by 'homing', a gene conversion process initiated by site-specific homing endonucleases. The VMA intein, which encodes the PI-SceI endonuclease in Saccharomyces cerevisiae, is present in several different yeast strains. Surprisingly, a wild wine yeast (DH1-1A) contains not only the intein(+) allele, but also an inteinless allele that has not undergone gene conversion. To elucidate how these two alleles co-exist, we characterized the endonuclease encoded by the DH1-1A intein(+) allele and the target site in the intein(-) allele. Sequence analysis reveals seven mutations in the 31 bp recognition sequence, none of which occurs at positions that are individually critical for activity. However, binding and cleavage of the sequence by PI-SceI is reduced 10-fold compared to the S.cerevisiae target. The PI-SceI analog encoded by the DH1-1A intein(+) allele contains 11 mutations at residues in the endonuclease and protein splicing domains. None affects protein splicing, but one, a R417Q substitution, accounts for most of the decrease in DNA cleavage and DNA binding activity of the DH1-1A protein. Loss of activity in the DH1-1A endonuclease and target site provides one explanation for co-existence of the intein(+) and intein(-) alleles.

Rapid evolution of the DNA-binding site in LAGLIDADG homing endonucleases

Sequence analysis of chloroplast and mitochondrial large subunit rRNA genes from over 75 green algae disclosed 28 new group I intron-encoded proteins carrying a single LAGLIDADG motif. These putative homing endonucleases form four subfamilies of homologous enzymes, with the members of each subfamily being encoded by introns sharing the same insertion site. We showed that four divergent endonucleases from the I-CreI subfamily cleave the same DNA substrates. Mapping of the 66 amino acids that are conserved among the members of this subfamily on the 3-dimensional structure of I-CreI bound to its recognition sequence revealed that these residues participate in protein folding, homodimerization, DNA recognition and catalysis. Surprisingly, only seven of the 21 I-CreI amino acids interacting with DNA are conserved, suggesting that I-CreI and its homologs use different subsets of residues to recognize the same DNA sequence. Our sequence comparison of all 45 single-LAGLIDADG proteins identified so far suggests that these proteins share related structures and that there is a weak pressure in each subfamily to maintain identical protein-DNA contacts. The high sequence variability we observed in the DNA-binding site of homologous LAGLIDADG endonucleases provides insight into how these proteins evolve new DNA specificity.

Functional alpha-fragment of beta-galactosidase can be expressed from the mobile group I intron PpLSU3 embedded in yeast pre-ribosomal RNA derived from the chromosomal rDNA locus

PpLSU3, a mobile group I intron found in the ribo-somal RNA genes of Physarum polycephalum, encodes the I-PpoI homing endonuclease. This enzyme represents one of the rare cases in nature where a protein is expressed from an RNA polymerase I transcript. Our previous results showed that the full length intron, but not a further processed species, is the messenger for I-PpoI, implying a role of the untranslated region (UTR) in gene expression. To study the function of the 3'-UTR in expression of the endonuclease and in splicing of the intron, we replaced the I-PpoI gene in PpLSU3 with the gene for the alpha-fragment of Escherichia coli beta-galactosidase, and then integrated this chimeric intron into all the chromosomal rDNA repeats of yeast. The resulting cells synthesized functional alpha-fragment, as evidenced by a complementation assay analogous to that used in E.coli. The beta-galactosidase activity thus provides an unusual and potentially valuable readout for Pol I transcription from chromosomal rDNA. This is the first example in which a eucaryotic homing endonuclease gene has been successfully replaced by a heterologous gene. Using deletion mutagenesis and a novel randomization approach with the alpha-fragment as a reporter, we found that a small segment of the 3'-UTR dramatically influences both splicing and protein expression.

Photocross-linking of the homing endonuclease PI-SceI to its recognition sequence

PI-SceI is an intein-encoded protein that belongs to the LAGLIDADG family of homing endonucleases. According to the crystal structure and mutational studies, this endonuclease consists of two domains, one responsible for protein splicing, the other for DNA cleavage, and both presumably for DNA binding. To define the DNA binding site of PI-SceI, photocross-linking was used to identify amino acid residues in contact with DNA. Sixty-three double-stranded oligodeoxynucleotides comprising the minimal recognition sequence and containing single 5-iodopyrimidine substitutions in almost all positions of the recognition sequence were synthesized and irradiated in the presence of PI-SceI with a helium/cadmium laser (325 nm). The best cross-linking yield (approximately 30%) was obtained with an oligodeoxynucleotide with a 5-iododeoxyuridine at position +9 in the bottom strand. The subsequent analysis showed that cross-linking had occurred with amino acid His-333, 6 amino acids after the second LAGLIDADG motif. With the H333A variant of PI-SceI or in the presence of excess unmodified oligodeoxynucleotide, no cross-linking was observed, indicating the specificity of the cross-linking reaction. Chemical modification of His residues in PI-SceI by diethylpyrocarbonate leads to a substantial reduction in the binding and cleavage activity of PI-SceI. This inactivation can be suppressed by substrate binding. This result further supports the finding that at least one His residue is in close contact to the DNA. Based on these and published results, conclusions are drawn regarding the DNA binding site of PI-SceI.

Crystal structure of the thermostable archaeal intron-encoded endonuclease I-DmoI

I-DmoI is a 22 kDa endonuclease encoded by an intron in the 23 S rRNA gene of the hyperthermophilic archaeon Desulfurococcus mobilis. The structure of I-DmoI has been determined to 2.2 A resolution using multi-wavelength anomalous diffraction techniques. I-DmoI, a protein of the LAGLIDADG motif family, represents the first structure of a freestanding endonuclease with two LAGLIDADG motifs, and the first of a thermostable homing endonuclease. I-DmoI consists of two similar alpha/beta domains (alphabetabetaalphabetabetaalpha) related by pseudo 2-fold symmetry. The LAGLIDADG motifs are located at the carboxy-terminal end of the first alpha-helix of each domain. These helices form a two-helix bundle at the interface between the domains and are perpendicular to a saddle-shaped DNA binding surface, formed by two four-stranded antiparallel beta-sheets. Despite substantially different sequences, the overall fold of I-DmoI is similar to that of two other LAGLIDADG proteins for which the structures are known, I-CreI and the endonuclease domain of PI-SceI. The three structures differ most in the loops connecting the beta-strands, relating to the respective DNA target site sizes and geometries. In addition, the absence of conserved residues surrounding the active site, other than those within the LAGLIDADG motif, is of mechanistic importance. Finally, the carboxy-terminal domain of I-DmoI is smaller and has a more irregular fold than the amino-terminal domain, which is more similar to I-CreI, a symmetric homodimeric endonuclease. This is reversed compared to PI-SceI, where the amino-terminal domain is more similar to carboxy-terminal domain of I-DmoI and to I-CreI, with interesting evolutionary implications.

I-PpoI and I-CreI homing site sequence degeneracy determined by random mutagenesis and sequential in vitro enrichment

Plasmid libraries containing partially randomized cleavage sites for the eukaryotic homing endonucleases I-PpoI and I-CreI were constructed, and sites that could be cleaved by I-PpoI or I-CreI were selectively recovered by successive cycles of cleavage and gel separation followed by religation and growth in Escherichia coli. Twenty-one different I-PpoI-sensitive homing sites, including the native homing site, were isolated. These sites were identical at four nucleotide positions within the 15 bp homing site, had a restricted pattern of base substitutions at the remaining 11 positions and displayed a preference for purines flanking the top strand of the homing site sequence. Twenty-one different I-CreI-sensitive homing sites, including the native site, were isolated. Ten nucleotide positions were identical in homing site variants that were I-CreI-sensitive and required the addition of SDS for efficient cleavage product release. Four of these ten positions were identical in homing sites that did not require SDS for product release. There was a preference for pyrimidines flanking the top strand of the homing site sequence. Three of the 24 I-CreI homing site nucleotide positions apparently lacked informational content, i. e. were permissive of cleavage when occupied by any nucleotide. These results suggest that I-PpoI and I-CreI make a large number of DNA-protein contacts across their homing site sequences, and that different subsets of these contacts may be sufficient to maintain a high degree of sequence-specific homing site recognition and cleavage. The sequential enrichment protocol we used should be useful for defining the sequence degeneracy and informational content of other homing endonuclease target sites.

Genetic analysis of the Chlamydomonas reinhardtii I-CreI mobile intron homing system in Escherichia coli

We have developed and used a genetic selection system in Escherichia coli to study functional requirements for homing site recognition and cleavage by a representative eukaryotic mobile intron endonuclease. The homing endonuclease, I-CreI, was originally isolated from the chloroplast of the unicellular green alga Chlamydomonas reinhardtii. I-CreI homing site mutants contained base pair substitutions or single base deletions that altered the rate of homing site cleavage and/or product release. I-CreI endonuclease mutants fell into six phenotypic classes that differed in in vivo activity, toxicity or genetic dominance. Inactivating mutations clustered in the N-terminal 60% of the I-CreI amino acid sequence, and two frameshift mutations were isolated that resulted in premature translation termination though retained partial activity. These mutations indicate that the N-terminal two-thirds of the I-CreI endonuclease is sufficient for homing site recognition and cleavage. Substitution mutations altered in four potential active site residues were examined: D20N, Q47H or R70A substitutions inactivated endonuclease activity, whereas S22A did not. The genetic approach we have taken complements phylogenetic and structural studies of mobile intron endonucleases and has provided new information on the mechanistic basis of I-CreI homing site recognition and cleavage.

Statistical modeling and analysis of the LAGLIDADG family of site-specific endonucleases and identification of an intein that encodes a site-specific endonuclease of the HNH family

The LAGLIDADG and HNH families of site-specific DNA endonucleases encoded by viruses, bacteriophages as well as archaeal, eucaryotic nuclear and organellar genomes are characterized by the sequence motifs 'LAGLIDADG' and 'HNH', respectively. These endonucleases have been shown to occur in different environments: LAGLIDADG endonucleases are found in inteins, archaeal and group I introns and as free standing open reading frames (ORFs); HNH endonucleases occur in group I and group II introns and as ORFs. Here, statistical models (hidden Markov models, HMMs) that encompass both the conserved motifs and more variable regions of these families have been created and employed to characterize known and potential new family members. A number of new, putative LAGLIDADG and HNH endonucleases have been identified including an intein-encoded HNH sequence. Analysis of an HMM-generated multiple alignment of 130 LAGLIDADG family members and the three-dimensional structure of the I- Cre I endonuclease has enabled definition of the core elements of the repeated domain (approximately 90 residues) that is present in this family of proteins. A conserved negatively charged residue is proposed to be involved in catalysis. Phylogenetic analysis of the two families indicates a lack of exchange of endonucleases between different mobile elements (environments) and between hosts from different phylogenetic kingdoms. However, there does appear to have been considerable exchange of endonuclease domains amongst elements of the same type. Such events are suggested to be important for the formation of elements of new specficity.

Purification, biochemical characterization and protein-DNA interactions of the I-CreI endonuclease produced in Escherichia coli

I- CreI is a member of the LAGLI-DADG family of homing nucleases; however, unlike most members of this family it contains only a single copy of this signature motif. I- CreI was over-expressed in Escherichia coli, and a simple purification protocol developed that gave reasonably pure protein in high yield. Size-exclusion chromatography and chemical cross-linking indicated that the protein is a dimer in solution. DNA cleavage by I- CreI was absolutely dependent on Mg2+(or Mn2+), and was inhibited by monovalent cations. I- CreI displayed a surprisingly high temperature optimum (>50 degrees C), with full activity occurring even at 70 degrees C. Interestingly, SDS was needed for efficient release of the cleavage products from the protein, indicating formation of very stable DNA-protein complexes. In contrast to these robust characteristics, purified I- CreI was unstable; however, it could be stabilized by the addition of either target or non-target DNA. Mobility shift assays revealed that I- CreI binds to DNA in the absence of Mg2+. Hydroxyl radical footprinting showed that I- CreI strongly protected the backbone of a continuous stretch of at least 12 nt on each strand that were shifted, relative to each other, by 2 bp in the 3'direction. Methylation protection and interference analyses were also performed, and together with the hydroxyl radical footprinting, indicate that I- CreI binds in both the major and minor grooves of its target DNA.

Mapping metal ions at the catalytic centres of two intron-encoded endonucleases

Divalent metal ions play a crucial role in forming the catalytic centres of DNA endonucleases. Substitution of Mg2+ ions by Fe2+ ions in two archaeal intron-encoded homing endonucleases, I-DmoI and I-PorI, yielded functional enzymes and enabled the generation of reactive hydroxyl radicals within the metal ion binding sites. Specific hydroxyl radical-induced cleavage was observed within, and immediately after, two conserved LAGLIDADG motifs in both proteins and at sites at, and near, the scissile phosphates of the corresponding DNA substrates. Titration of Fe2+-containing protein-DNA complexes with Ca2+ ions, which are unable to support endonucleolytic activity, was performed to distinguish between the individual metal ions in the complex. Mutations of single amino acids in this region impaired catalytic activity and caused the preferential loss of a subset of hydroxyl radical cleavages in both the protein and the DNA substrate, suggesting an active role in metal ion coordination for these amino acids. The data indicate that the endonucleases cleave their DNA substrates as monomeric enzymes, and contain a minimum of four divalent metal ions located at or near the catalytic centres of each endonuclease. The metal ions involved in cleaving the coding and the non-coding strand are positioned immediately after the N- and C-terminally located LAGLIDADG motifs, respectively. The dual protein/nucleic acid footprinting approach described here is generally applicable to other protein-nucleic acid complexes when the natural metal ion can be replaced by Fe2+.