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

The nuclease-associated short prokaryotic Argonaute system nonspecifically degrades DNA upon activation by target recognition

Prokaryotic Argonautes (pAgos) play a vital role in host defense by utilizing short nucleic acid guides to recognize and target complementary nucleic acids. Despite being the majority of pAgos, short pAgos have only recently received attention. Short pAgos are often associated with proteins containing an APAZ domain and a nuclease domain including DUF4365, SMEK, or HNH domain. In contrast to long pAgos that specifically cleave the target DNA, our study demonstrates that the short pAgo from Thermocrispum municipal, along with its associated DUF4365-APAZ protein, forms a heterodimeric complex. Upon RNA-guided target DNA recognition, this complex is activated to nonspecifically cleave DNA. Additionally, we found that the TmuRE-Ago complex shows a preference for 5'-OH guide RNA, specifically requires a uridine nucleotide at the 5' end of the guide RNA, and is sensitive to single-nucleotide mismatches between the guide RNA and target DNA. Based on its catalytic properties, our study has established a novel nucleic acid detection method and demonstrated its feasibility. This study not only expands our understanding of the defense mechanism employed by short pAgo systems but also suggests their potential applications in nucleic acid detection.

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4's. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Post-genomic era in agriculture and veterinary science: successful and proposed application of genetic targeting technologies

Gene editing tools have become an indispensable part of research into the fundamental aspects of cell biology. With a vast body of literature having been generated based on next generation sequencing technologies, keeping track of this ever-growing body of information remains challenging. This necessitates the translation of genomic data into tangible applications. In order to address this objective, the generated Next Generation Sequencing (NGS) data forms the basis for targeted genome editing strategies, employing known enzymes of various cellular machinery, in generating organisms with specifically selected phenotypes. This review focuses primarily on CRISPR/Cas9 technology in the context of its advantages over Zinc finger proteins (ZNF) and Transcription activator-like effector nucleases (TALEN) and meganucleases mutagenesis strategies, for use in agricultural and veterinary applications. This review will describe the application of CRISPR/Cas9 in creating modified organisms with custom-made properties, without the undesired non-targeted effects associated with virus vector vaccines and bioactive molecules produced in bacterial systems. Examples of the successful and unsuccessful applications of this technology to plants, animals and microorganisms are provided, as well as an in-depth look into possible future trends and applications in vaccine development, disease resistance and enhanced phenotypic traits will be discussed.

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4’s. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Genome-Wide Identification and Transcriptional Analysis of Arabidopsis DUF506 Gene Family

The Domain of unknown function 506 (DUF506) family, which belongs to the PD-(D/E)XK nuclease superfamily, has not been functionally characterized. In this study, 266 DUF506 domain-containing genes were identified from algae, mosses, and land plants showing their wide occurrence in photosynthetic organisms. Bioinformatics analysis identified 211 high-confidence DUF506 genes across 17 representative land plant species. Phylogenetic modeling classified three groups of plant DUF506 genes that suggested functional preservation among the groups based on conserved gene structure and motifs. Gene duplication and Ka/Ks evolutionary rates revealed that DUF506 genes are under purifying positive selection pressure. Subcellular protein localization analysis revealed that DUF506 proteins were present in different organelles. Transcript analyses showed that 13 of the Arabidopsis DUF506 genes are ubiquitously expressed in various tissues and respond to different abiotic stresses and ABA treatment. Protein-protein interaction network analysis using the STRING-DB, AtPIN (Arabidopsis thaliana Protein Interaction Network), and AI-1 (Arabidopsis Interactome-1) tools indicated that AtDUF506s potentially interact with iron-deficiency response proteins, salt-inducible transcription factors, or calcium sensors (calmodulins), implying that DUF506 genes have distinct biological functions including responses to environmental stimuli, nutrient-deficiencies, and participate in Ca(2+) signaling. Current results provide insightful information regarding the molecular features of the DUF506 family in plants, to support further functional characterizations.

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.

Genome editing approaches: manipulating of lovastatin and taxol synthesis of filamentous fungi by CRISPR/Cas9 system

Filamentous fungi are prolific repertoire of structurally diverse secondary metabolites of remarkable biological activities such as lovastatin and paclitaxel that have been approved by FDA as drugs for hypercholesterolemia and cancer treatment. The clusters of genes encoding lovastatin and paclitaxel are cryptic at standard laboratory cultural conditions (Kennedy et al. Science 284:1368-1372, 1999; Bergmann et al. Nature Chem Biol 3:213-217, 2007). The expression of these genes might be triggered in response to nutritional and physical conditions; nevertheless, the overall yield of these metabolites does not match the global need. Consequently, overexpression of the downstream limiting enzymes and/or blocking the competing metabolic pathways of these metabolites could be the most successful technologies to enhance their yield. This is the first review summarizing the different strategies implemented for fungal genome editing, molecular regulatory mechanisms, and prospective of clustered regulatory interspaced short palindromic repeat/Cas9 system in metabolic engineering of fungi to improve their yield of lovastatin and taxol to industrial scale. Thus, elucidating the putative metabolic pathways in fungi for overproduction of lovastatin and taxol was the ultimate objective of this review.

Rpn (YhgA-Like) Proteins of Escherichia coli K-12 and Their Contribution to RecA-Independent Horizontal Transfer

Bacteria use a variety of DNA-mobilizing enzymes to facilitate environmental niche adaptation via horizontal gene transfer. This has led to real-world problems, like the spread of antibiotic resistance, yet many mobilization proteins remain undefined. In the study described here, we investigated the uncharacterized family of YhgA-like transposase_31 (Pfam PF04754) proteins. Our primary focus was the genetic and biochemical properties of the five Escherichia coli K-12 members of this family, which we designate RpnA to RpnE, where Rpn represents recombination-promoting nuclease. We employed a conjugal system developed by our lab that demanded RecA-independent recombination following transfer of chromosomal DNA. Overexpression of RpnA (YhgA), RpnB (YfcI), RpnC (YadD), and RpnD (YjiP) increased RecA-independent recombination, reduced cell viability, and induced the expression of reporter of DNA damage. For the exemplar of the family, RpnA, mutational changes in proposed catalytic residues reduced or abolished all three phenotypes in concert. In vitro, RpnA displayed magnesium-dependent, calcium-stimulated DNA endonuclease activity with little, if any, sequence specificity and a preference for double-strand cleavage. We propose that Rpn/YhgA-like family nucleases can participate in gene acquisition processes.IMPORTANCE Bacteria adapt to new environments by obtaining new genes from other bacteria. Here, we characterize a set of genes that can promote the acquisition process by a novel mechanism. Genome comparisons had suggested the horizontal spread of the genes for the YhgA-like family of proteins through bacteria. Although annotated as transposase_31, no member of the family has previously been characterized experimentally. We show that four Escherichia coli K-12 paralogs contribute to a novel RecA-independent recombination mechanism in vivo For RpnA, we demonstrate in vitro action as a magnesium-dependent, calcium-stimulated nonspecific DNA endonuclease. The cleavage products are capable of providing priming sites for DNA polymerase, which can enable DNA joining by primer-template switching.

Genomics Approaches For Improving Salinity Stress Tolerance in Crop Plants

Salinity is one of the major factors which reduces crop production worldwide. Plant responses to salinity are highly complex and involve a plethora of genes. Due to its multigenicity, it has been difficult to attain a complete understanding of how plants respond to salinity. Genomics has progressed tremendously over the past decade and has played a crucial role towards providing necessary knowledge for crop improvement. Through genomics, we have been able to identify and characterize the genes involved in salinity stress response, map out signaling pathways and ultimately utilize this information for improving the salinity tolerance of existing crops. The use of new tools, such as gene pyramiding, in genetic engineering and marker assisted breeding has tremendously enhanced our ability to generate stress tolerant crops. Genome editing technologies such as Zinc finger nucleases, TALENs and CRISPR/Cas9 also provide newer and faster avenues for plant biologists to generate precisely engineered crops.

Single-Cell (Meta-)Genomics of a Dimorphic Candidatus Thiomargarita nelsonii Reveals Genomic Plasticity

The genus Thiomargarita includes the world's largest bacteria. But as uncultured organisms, their physiology, metabolism, and basis for their gigantism are not well understood. Thus, a genomics approach, applied to a single Candidatus Thiomargarita nelsonii cell was employed to explore the genetic potential of one of these enigmatic giant bacteria. The Thiomargarita cell was obtained from an assemblage of budding Ca. T. nelsonii attached to a provannid gastropod shell from Hydrate Ridge, a methane seep offshore of Oregon, USA. Here we present a manually curated genome of Bud S10 resulting from a hybrid assembly of long Pacific Biosciences and short Illumina sequencing reads. With respect to inorganic carbon fixation and sulfur oxidation pathways, the Ca. T. nelsonii Hydrate Ridge Bud S10 genome was similar to marine sister taxa within the family Beggiatoaceae. However, the Bud S10 genome contains genes suggestive of the genetic potential for lithotrophic growth on arsenite and perhaps hydrogen. The genome also revealed that Bud S10 likely respires nitrate via two pathways: a complete denitrification pathway and a dissimilatory nitrate reduction to ammonia pathway. Both pathways have been predicted, but not previously fully elucidated, in the genomes of other large, vacuolated, sulfur-oxidizing bacteria. Surprisingly, the genome also had a high number of unusual features for a bacterium to include the largest number of metacaspases and introns ever reported in a bacterium. Also present, are a large number of other mobile genetic elements, such as insertion sequence (IS) transposable elements and miniature inverted-repeat transposable elements (MITEs). In some cases, mobile genetic elements disrupted key genes in metabolic pathways. For example, a MITE interrupts hupL, which encodes the large subunit of the hydrogenase in hydrogen oxidation. Moreover, we detected a group I intron in one of the most critical genes in the sulfur oxidation pathway, dsrA. The dsrA group I intron also carried a MITE sequence that, like the hupL MITE family, occurs broadly across the genome. The presence of a high degree of mobile elements in genes central to Thiomargarita's core metabolism has not been previously reported in free-living bacteria and suggests a highly mutable genome.

Structural characterization of MepB from Staphylococcus aureus reveals homology to endonucleases

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

Homing endonucleases from mobile group I introns: discovery to genome engineering

Homing endonucleases are highly specific DNA cleaving enzymes that are encoded within genomes of all forms of microbial life including phage and eukaryotic organelles. These proteins drive the mobility and persistence of their own reading frames. The genes that encode homing endonucleases are often embedded within self-splicing elements such as group I introns, group II introns and inteins. This combination of molecular functions is mutually advantageous: the endonuclease activity allows surrounding introns and inteins to act as invasive DNA elements, while the splicing activity allows the endonuclease gene to invade a coding sequence without disrupting its product. Crystallographic analyses of representatives from all known homing endonuclease families have illustrated both their mechanisms of action and their evolutionary relationships to a wide range of host proteins. Several homing endonucleases have been completely redesigned and used for a variety of genome engineering applications. Recent efforts to augment homing endonucleases with auxiliary DNA recognition elements and/or nucleic acid processing factors has further accelerated their use for applications that demand exceptionally high specificity and activity.

Homing endonucleases: from genetic anomalies to programmable genomic clippers

Homing endonucleases are strong drivers of genetic exchange and horizontal transfer of both their own genes and their local genetic environment. The mechanisms that govern the function and evolution of these genetic oddities have been well documented over the past few decades at the genetic, biochemical, and structural levels. This wealth of information has led to the manipulation and reprogramming of the endonucleases and to their exploitation in genome editing for use as therapeutic agents, for insect vector control and in agriculture. In this chapter we summarize the molecular properties of homing endonucleases and discuss their strengths and weaknesses in genome editing as compared to other site-specific nucleases such as zinc finger endonucleases, TALEN, and CRISPR-derived endonucleases.

Group I introns in Staphylococcus bacteriophages

The abundance of group I introns, intragenic RNA sequences capable of self-splicing, in Gram-positive bacteriophage genomes, is illustrated by various new group I introns recently described in Staphylococcus phage genomes. These introns were found to interrupt DNA metabolism genes as well as late genes. These group I introns often code for homing endonucleases, which promote lateral transfer of group I introns, thereby enabling spread through a population. Homing endonucleases encoded by group I introns in Staphylococcus phage genomes were predicted to belong to the GIY-YIG, LAGLIDADG, HNH or EDxHD family of endonucleases. The group I intron distribution in Staphylococcus phage genomes exemplifies the homology between these introns as well as the encoded endonucleases. Despite several suggested functions, the role of group I introns in bacteriophages remains unclear or might be nonexistent. However, transcriptome analysis might provide additional information to elucidate the possible purpose of group I introns in phage genomes.

The mitochondrial LSU rRNA group II intron of Ustilago maydis encodes an active homing endonuclease likely involved in intron mobility

Background

The a2 mating type locus gene lga2 is critical for uniparental mitochondrial DNA inheritance during sexual development of Ustilago maydis. Specifically, the absence of lga2 results in biparental inheritance, along with efficient transfer of intronic regions in the large subunit rRNA gene between parental molecules. However, the underlying role of the predicted LAGLIDADG homing endonuclease gene I-UmaI located within the group II intron LRII1 has remained unresolved.

Methodology/Principal Findings

We have investigated the enzymatic activity of I-UmaI in vitro based on expression of a tagged full-length and a naturally occurring mutant derivative, which harbors only the N-terminal LAGLIDADG domain. This confirmed Mg²⁺-dependent endonuclease activity and cleavage at the LRII1 insertion site to generate four base pair extensions with 3′ overhangs. Specifically, I-UmaI recognizes an asymmetric DNA sequence with a minimum length of 14 base pairs (5′-GACGGGAAGACCCT-3′) and tolerates subtle base pair substitutions within the homing site. Enzymatic analysis of the mutant variant indicated a correlation between the activity in vitro and intron homing. Bioinformatic analyses revealed that putatively functional or former functional I-UmaI homologs are confined to a few members within the Ustilaginales and Agaricales, including the phylogenetically distant species Lentinula edodes, and are linked to group II introns inserted into homologous positions in the LSU rDNA.

Conclusions/Significance

The present data provide strong evidence that intron homing efficiently operates under conditions of biparental inheritance in U. maydis. Conversely, uniparental inheritance may be critical to restrict the transmission of mobile introns. Bioinformatic analyses suggest that I-UmaI-associated introns have been acquired independently in distant taxa and are more widespread than anticipated from available genomic data.

I-PfoP3I: a novel nicking HNH homing endonuclease encoded in the group I intron of the DNA polymerase gene in Phormidium foveolarum phage Pf-WMP3

Homing endonucleases encoded in a group I self-splicing intron in a protein-coding gene in cyanophage genomes have not been reported, apart from some free-standing homing edonucleases. In this study, a nicking DNA endonuclease, I-PfoP3I, encoded in a group IA2 intron in the DNA polymerase gene of a T7-like cyanophage Pf-WMP3, which infects the freshwater cyanobacterium Phormidium foveolarum is described. The Pf-WMP3 intron splices efficiently in vivo and self-splices in vitro simultaneously during transcription. I-PfoP3I belongs to the HNH family with an unconventional C-terminal HNH motif. I-PfoP3I nicks the intron-minus Pf-WMP3 DNA polymerase gene more efficiently than the Pf-WMP4 DNA polymerase gene that lacks any intervening sequence in vitro, indicating the variable capacity of I-PfoP3I. I-PfoP3I cleaves 4 nt upstream of the intron insertion site on the coding strand of EXON 1 on both intron-minus Pf-WMP3 and Pf-WMP4 DNA polymerase genes. Using an in vitro cleavage assay and scanning deletion mutants of the intronless target site, the minimal recognition site was determined to be a 14 bp region downstream of the cut site. I-PfoP3I requires Mg²⁺, Ca²⁺ or Mn²⁺ for nicking activity. Phylogenetic analysis suggests that the intron and homing endonuclease gene elements might be inserted in Pf-WMP3 genome individually after differentiation from Pf-WMP4. To our knowledge, this is the first report of the presence of a group I self-splicing intron encoding a functional homing endonuclease in a protein-coding gene in a cyanophage genome.

Homing endonucleases: DNA scissors on a mission

Buried within the genomes of many microorganisms are genetic elements that encode rare-cutting homing endonucleases that assist in the mobility of the elements that encode them, such as the self-splicing group I and II introns and in some cases inteins. There are several different families of homing endonucleases and their ability to initiate and target specific sequences for lateral transfers makes them attractive reagents for gene targeting. Homing endonucleases have been applied in promoting DNA modification or genome editing such as gene repair or "gene knockouts". This review examines the categories of homing endonucleases that have been described so far and their possible applications to biotechnology. Strategies to engineer homing endonucleases to alter target site specificities will also be addressed. Alternatives to homing endonucleases such as zinc finger nucleases, transcription activator-like effector nucleases, triplex forming oligonucleotide nucleases, and targetrons are also briefly discussed.

Gene flow and biological conflict systems in the origin and evolution of eukaryotes

The endosymbiotic origin of eukaryotes brought together two disparate genomes in the cell. Additionally, eukaryotic natural history has included other endosymbiotic events, phagotrophic consumption of organisms, and intimate interactions with viruses and endoparasites. These phenomena facilitated large-scale lateral gene transfer and biological conflicts. We synthesize information from nearly two decades of genomics to illustrate how the interplay between lateral gene transfer and biological conflicts has impacted the emergence of new adaptations in eukaryotes. Using apicomplexans as example, we illustrate how lateral transfer from animals has contributed to unique parasite-host interfaces comprised of adhesion- and O-linked glycosylation-related domains. Adaptations, emerging due to intense selection for diversity in the molecular participants in organismal and genomic conflicts, being dispersed by lateral transfer, were subsequently exapted for eukaryote-specific innovations. We illustrate this using examples relating to eukaryotic chromatin, RNAi and RNA-processing systems, signaling pathways, apoptosis and immunity. We highlight the major contributions from catalytic domains of bacterial toxin systems to the origin of signaling enzymes (e.g., ADP-ribosylation and small molecule messenger synthesis), mutagenic enzymes for immune receptor diversification and RNA-processing. Similarly, we discuss contributions of bacterial antibiotic/siderophore synthesis systems and intra-genomic and intra-cellular selfish elements (e.g., restriction-modification, mobile elements and lysogenic phages) in the emergence of chromatin remodeling/modifying enzymes and RNA-based regulation. We develop the concept that biological conflict systems served as evolutionary “nurseries” for innovations in the protein world, which were delivered to eukaryotes via lateral gene flow to spur key evolutionary innovations all the way from nucleogenesis to lineage-specific adaptations.

Polymorphic toxin systems: Comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

BACKGROUND

Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflicts. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis.

RESULTS

Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses.

CONCLUSIONS

Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.

Sequence, structure and functional diversity of PD-(D/E)XK phosphodiesterase superfamily

Proteins belonging to PD-(D/E)XK phosphodiesterases constitute a functionally diverse superfamily with representatives involved in replication, restriction, DNA repair and tRNA-intron splicing. Their malfunction in humans triggers severe diseases, such as Fanconi anemia and Xeroderma pigmentosum. To date there have been several attempts to identify and classify new PD-(D/E)KK phosphodiesterases using remote homology detection methods. Such efforts are complicated, because the superfamily exhibits extreme sequence and structural divergence. Using advanced homology detection methods supported with superfamily-wide domain architecture and horizontal gene transfer analyses, we provide a comprehensive reclassification of proteins containing a PD-(D/E)XK domain. The PD-(D/E)XK phosphodiesterases span over 21,900 proteins, which can be classified into 121 groups of various families. Eleven of them, including DUF4420, DUF3883, DUF4263, COG5482, COG1395, Tsp45I, HaeII, Eco47II, ScaI, HpaII and Replic_Relax, are newly assigned to the PD-(D/E)XK superfamily. Some groups of PD-(D/E)XK proteins are present in all domains of life, whereas others occur within small numbers of organisms. We observed multiple horizontal gene transfers even between human pathogenic bacteria or from Prokaryota to Eukaryota. Uncommon domain arrangements greatly elaborate the PD-(D/E)XK world. These include domain architectures suggesting regulatory roles in Eukaryotes, like stress sensing and cell-cycle regulation. Our results may inspire further experimental studies aimed at identification of exact biological functions, specific substrates and molecular mechanisms of reactions performed by these highly diverse proteins.

Structural, functional and evolutionary relationships between homing endonucleases and proteins from their host organisms

Homing endonucleases (HEs) are highly specific DNA-cleaving enzymes that are encoded by invasive DNA elements (usually mobile introns or inteins) within the genomes of phage, bacteria, archea, protista and eukaryotic organelles. Six unique structural HE families, that collectively span four distinct nuclease catalytic motifs, have been characterized to date. Members of each family display structural homology and functional relationships to a wide variety of proteins from various organisms. The biological functions of those proteins are highly disparate and include non-specific DNA-degradation enzymes, restriction endonucleases, DNA-repair enzymes, resolvases, intron splicing factors and transcription factors. These relationships suggest that modern day HEs share common ancestors with proteins involved in genome fidelity, maintenance and gene expression. This review summarizes the results of structural studies of HEs and corresponding proteins from host organisms that have illustrated the manner in which these factors are related.

Towards artificial metallonucleases for gene therapy: recent advances and new perspectives

The process of DNA targeting or repair of mutated genes within the cell, induced by specifically positioned double-strand cleavage of DNA near the mutated sequence, can be applied for gene therapy of monogenic diseases. For this purpose, highly specific artificial metallonucleases are developed. They are expected to be important future tools of modern genetics. The present state of art and strategies of research are summarized, including protein engineering and artificial ‘chemical’ nucleases. From the results, we learn about the basic role of the metal ions and the various ligands, and about the DNA binding and cleavage mechanism. The results collected provide useful guidance for engineering highly controlled enzymes for use in gene therapy.

Activity, specificity and structure of I-Bth0305I: a representative of a new homing endonuclease family

Novel family of putative homing endonuclease genes was recently discovered during analyses of metagenomic and genomic sequence data. One such protein is encoded within a group I intron that resides in the recA gene of the Bacillus thuringiensis 0305ϕ8–36 bacteriophage. Named I-Bth0305I, the endonuclease cleaves a DNA target in the uninterrupted recA gene at a position immediately adjacent to the intron insertion site. The enzyme displays a multidomain, homodimeric architecture and footprints a DNA region of ∼60 bp. Its highest specificity corresponds to a 14-bp pseudopalindromic sequence that is directly centered across the DNA cleavage site. Unlike many homing endonucleases, the specificity profile of the enzyme is evenly distributed across much of its target site, such that few single base pair substitutions cause a significant decrease in cleavage activity. A crystal structure of its C-terminal domain confirms a nuclease fold that is homologous to very short patch repair (Vsr) endonucleases. The domain architecture and DNA recognition profile displayed by I-Bth0305I, which is the prototype of a homing lineage that we term the ‘EDxHD’ family, are distinct from previously characterized homing endonucleases.

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.

Homing endonucleases: from microbial genetic invaders to reagents for targeted DNA modification

Homing endonucleases are microbial DNA-cleaving enzymes that mobilize their own reading frames by generating double strand breaks at specific genomic invasion sites. These proteins display an economy of size, and yet recognize long DNA sequences (typically 20 to 30 base pairs). They exhibit a wide range of fidelity at individual nucleotide positions in a manner that is strongly influenced by host constraints on the coding sequence of the targeted gene. The activity of these proteins leads to site-specific recombination events that can result in the insertion, deletion, mutation, or correction of DNA sequences. Over the past fifteen years, the crystal structures of representatives from several homing endonuclease families have been solved, and methods have been described to create variants of these enzymes that cleave novel DNA targets. Engineered homing endonucleases proteins are now being used to generate targeted genomic modifications for a variety of biotech and medical applications.

Meganucleases and other tools for targeted genome engineering: perspectives and challenges for gene therapy

The importance of safer approaches for gene therapy has been underscored by a series of severe adverse events (SAEs) observed in patients involved in clinical trials for Severe Combined Immune Deficiency Disease (SCID) and Chromic Granulomatous Disease (CGD). While a new generation of viral vectors is in the process of replacing the classical gamma-retrovirus-based approach, a number of strategies have emerged based on non-viral vectorization and/or targeted insertion aimed at achieving safer gene transfer. Currently, these methods display lower efficacies than viral transduction although many of them can yield more than 1% of engineered cells in vitro. Nuclease-based approaches, wherein an endonuclease is used to trigger site-specific genome editing, can significantly increase the percentage of targeted cells. These methods therefore provide a real alternative to classical gene transfer as well as gene editing. However, the first endonuclease to be in clinic today is not used for gene transfer, but to inactivate a gene (CCR5) required for HIV infection. Here, we review these alternative approaches, with a special emphasis on meganucleases, a family of naturally occurring rare-cutting endonucleases, and speculate on their current and future potential.

The I-CreI meganuclease and its engineered derivatives: applications from cell modification to gene therapy

Meganucleases (MNs) are highly specific enzymes that can induce homologous recombination in different types of cells, including mammalian cells. Consequently, these enzymes are used as scaffolds for the development of custom gene-targeting tools for gene therapy or cell-line development. Over the past 15 years, the high resolution X-ray structures of several MNs from the LAGLIDADG family have improved our understanding of their protein-DNA interaction and mechanism of DNA cleavage. By developing and utilizing high-throughput screening methods to test a large number of variant-target combinations, we have been able to re-engineer scores of I-CreI derivatives into custom enzymes that target a specific DNA sequence of interest. Such customized MNs, along with wild-type ones, have allowed for exploring a large range of biotechnological applications, including protein-expression cell-line development, genetically modified plants and animals and therapeutic applications such as targeted gene therapy as well as a novel class of antivirals.

Natural and engineered nicking endonucleases--from cleavage mechanism to engineering of strand-specificity

Restriction endonucleases (REases) are highly specific DNA scissors that have facilitated the development of modern molecular biology. Intensive studies of double strand (ds) cleavage activity of Type IIP REases, which recognize 4–8 bp palindromic sequences, have revealed a variety of mechanisms of molecular recognition and catalysis. Less well-studied are REases which cleave only one of the strands of dsDNA, creating a nick instead of a ds break. Naturally occurring nicking endonucleases (NEases) range from frequent cutters such as Nt.CviPII (^CCD; ^ denotes the cleavage site) to rare-cutting homing endonucleases (HEases) such as I-HmuI. In addition to these bona fida NEases, individual subunits of some heterodimeric Type IIS REases have recently been shown to be natural NEases. The discovery and characterization of more REases that recognize asymmetric sequences, particularly Types IIS and IIA REases, has revealed recognition and cleavage mechanisms drastically different from the canonical Type IIP mechanisms, and has allowed researchers to engineer highly strand-specific NEases. Monomeric LAGLIDADG HEases use two separate catalytic sites for cleavage. Exploitation of this characteristic has also resulted in useful nicking HEases. This review aims at providing an overview of the cleavage mechanisms of Types IIS and IIA REases and LAGLIDADG HEases, the engineering of their nicking variants, and the applications of NEases and nicking HEases.

The crystal structure of D212 from sulfolobus spindle-shaped virus ragged hills reveals a new member of the PD-(D/E)XK nuclease superfamily

Structural studies have made significant contributions to our understanding of Sulfolobus spindle-shaped viruses (Fuselloviridae), an important model system for archaeal viruses. Continuing these efforts, we report the structure of D212 from Sulfolobus spindle-shaped virus Ragged Hills. The overall fold and conservation of active site residues place D212 in the PD-(D/E)XK nuclease superfamily. The greatest structural similarity is found to the archaeal Holliday junction cleavage enzymes, strongly suggesting a role in DNA replication, repair, or recombination. Other roles associated with nuclease activity are also considered.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

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

A unified genetic, computational and experimental framework identifies functionally relevant residues of the homing endonuclease I-BmoI

Insight into protein structure and function is best obtained through a synthesis of experimental, structural and bioinformatic data. Here, we outline a framework that we call MUSE (mutual information, unigenic evolution and structure-guided elucidation), which facilitated the identification of previously unknown residues that are relevant for function of the GIY-YIG homing endonuclease I-BmoI. Our approach synthesizes three types of data: mutual information analyses that identify co-evolving residues within the GIY-YIG catalytic domain; a unigenic evolution strategy that identifies hyper- and hypo-mutable residues of I-BmoI; and interpretation of the unigenic and co-evolution data using a homology model. In particular, we identify novel positions within the GIY-YIG domain as functionally important. Proof-of-principle experiments implicate the non-conserved I71 as functionally relevant, with an I71N mutant accumulating a nicked cleavage intermediate. Moreover, many additional positions within the catalytic, linker and C-terminal domains of I-BmoI were implicated as important for function. Our results represent a platform on which to pursue future studies of I-BmoI and other GIY-YIG-containing proteins, and demonstrate that MUSE can successfully identify novel functionally critical residues that would be ignored in a traditional structure-function analysis within an extensively studied small domain of approximately 90 amino acids.

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.

Characterization of Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease, reveals a unique mode of DNA binding, helical distortion, and cleavage compared with a canonical LAGLIDADG homing endonuclease

Mycobacterium leprae, which has undergone reductive evolution leaving behind a minimal set of essential genes, has retained intervening sequences in four of its genes implicating a vital role for them in the survival of the leprosy bacillus. A single in-frame intervening sequence has been found embedded within its recA gene. Comparison of the M. leprae recA intervening sequence with the known intervening sequences indicated that it has the consensus amino acid sequence necessary for being a LAGLIDADG-type homing endonuclease. In light of massive gene decay and function loss in the leprosy bacillus, we sought to investigate whether its recA intervening sequence encodes a catalytically active homing endonuclease. Here we show that the purified M. leprae RecA intein (PI-MleI) binds to cognate DNA and displays endonuclease activity in the presence of alternative divalent cations, Mg2+ or Mn2+. A combination of approaches, including four complementary footprinting assays such as DNase I, copper-phenanthroline, methylation protection, and KMnO4, enhancement of 2-aminopurine fluorescence, and mapping of the cleavage site revealed that PI-MleI binds to cognate DNA flanking its insertion site, induces helical distortion at the cleavage site, and generates two staggered double strand breaks. Taken together, these results implicate that PI-MleI possesses a modular structure with separate domains for DNA target recognition and cleavage, each with distinct sequence preferences. From a biological standpoint, it is tempting to speculate that our findings have implications for understanding the evolution of the LAGLIDADG family of homing endonucleases.

A likely pathway for formation of mobile group I introns

Mobile group I introns are RNA splicing elements that have been invaded by endonuclease genes. These endonucleases facilitate intron mobility by a unidirectional, duplicative gene-conversion process known as homing [1]. Survival of the invading endonuclease depends upon its ability to promote intron mobility. Therefore, the endonuclease must either quickly change its cleavage specificity to match the site of intron insertion, or it must already be preadapted to cleave this sequence. Here we show that the group I intron in the DNA polymerase gene of T7-like bacteriophage PhiI is mobile, dependent upon its intronic HNH homing endonuclease gene, I-TslI. We also show that gene 5.3 of phage T3, located adjacent to its intronless DNA polymerase gene, is a homologous homing endonuclease gene whose protein product initiates efficient spread of gene 5.3 into empty sites in related phages. Both of these endonucleases cleave intronless DNA polymerase genes at identical positions. This shared feature between an intronic and free-standing endonuclease is unprecedented. Based on this evidence, we propose that introns and their homing endonucleases evolve separately to target the same highly conserved sequences, uniting afterwards to create a composite mobile element.

Selfish DNA: homing endonucleases find a home

Self-splicing group I introns come in two flavours - those with a homing endonuclease to promote mobility of the intron, and those without an endonuclease. How homing endonucleases and self-splicing introns associate to form a composite selfish genetic element is a question of long-standing interest. Recent work has revealed that a shared characteristic of both introns and endonucleases, the targeting of conserved sequences, may provide the impetus for the evolution of composite mobile genetic elements.

Fractured genes: a novel genomic arrangement involving new split inteins and a new homing endonuclease family

Inteins are genetic elements, inserted in-frame into protein-coding genes, whose products catalyze their removal from the protein precursor via a protein-splicing reaction. Intein domains can be split into two fragments and still ligate their flanks by a trans-protein-splicing reaction. A bioinformatic analysis of environmental metagenomic data revealed 26 different loci with a novel genomic arrangement. In each locus, a conserved enzyme coding region is broken in two by a split intein, with a free-standing endonuclease gene inserted in between. Eight types of DNA synthesis and repair enzymes have this ‘fractured’ organization. The new types of naturally split-inteins were analyzed in comparison to known split-inteins. Some loci include apparent gene control elements brought in with the endonuclease gene. A newly predicted homing endonuclease family, related to very-short patch repair (Vsr) endonucleases, was found in half of the loci. These putative homing endonucleases also appear in group-I introns, and as stand-alone inserts in the absence of surrounding intervening sequences. The new fractured genes organization appears to be present mainly in phage, shows how endonucleases can integrate into inteins, and may represent a missing link in the evolution of gene breaking in general, and in the creation of split-inteins in particular.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.

Activity and specificity of the bacterial PD-(D/E)XK homing endonuclease I-Ssp6803I

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

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background

Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing.

Results

Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts.

Conclusion

The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

Structures of the rare-cutting restriction endonuclease NotI reveal a unique metal binding fold involved in DNA binding

The structure of the rare-cutting restriction endonuclease NotI, which recognizes the 8 bp target 5'-GCGGCCGC-3', has been solved with and without bound DNA. Because of its specificity (recognizing a site that occurs once per 65 kb), NotI is used to generate large genomic fragments and to map DNA methylation status. NotI contains a unique metal binding fold, found in a variety of putative endonucleases, occupied by an iron atom coordinated within a tetrahedral Cys4 motif. This domain positions nearby protein elements for DNA recognition, and serves a structural role. While recognition of the central six base pairs of the target is accomplished via a saturated hydrogen bond network typical of restriction enzymes, the most peripheral base pairs are engaged in a single direct contact in the major groove, reflecting reduced pressure to recognize those positions. NotI may represent an evolutionary intermediate between mobile endonucleases (which recognize longer target sites) and canonical restriction endonucleases.

Unconventional GIY-YIG homing endonuclease encoded in group I introns in closely related strains of the Bacillus cereus group

Several group I introns have been previously found in strains of the Bacillus cereus group at three different insertion sites in the nrdE gene of the essential nrdIEF operon coding for ribonucleotide reductase. Here, we identify an uncharacterized group IA intron in the nrdF gene in 12 strains of the B. cereus group and show that the pre-mRNA is efficiently spliced. The Bacillus thuringiensis ssp. pakistani nrdF intron encodes a homing endonuclease, denoted I-BthII, with an unconventional GIY-(X)₈-YIG motif that cleaves an intronless nrdF gene 7 nt upstream of the intron insertion site, producing 2-nt 3′ extensions. We also found four additional occurrences of two of the previously reported group I introns in the nrdE gene of 25 sequenced B. thuringiensis and one B. cereus strains, and one non-annotated group I intron at a fourth nrdE insertion site in the B. thuringiensis ssp. Al Hakam sequenced genome. Two strains contain introns in both the nrdE and the nrdF genes. Phylogenetic studies of the nrdIEF operon from 39 strains of the B. cereus group suggest several events of horizontal gene transfer for two of the introns found in this operon.

Thermodynamics of DNA target site recognition by homing endonucleases

The thermodynamic profiles of target site recognition have been surveyed for homing endonucleases from various structural families. Similar to DNA-binding proteins that recognize shorter target sites, homing endonucleases display a narrow range of binding free energies and affinities, mediated by structural interactions that balance the magnitude of enthalpic and entropic forces. While the balance of ΔH and TΔS are not strongly correlated with the overall extent of DNA bending, unfavorable ΔH_binding is associated with unstacking of individual base steps in the target site. The effects of deleterious basepair substitutions in the optimal target sites of two LAGLIDADG homing endonucleases, and the subsequent effect of redesigning one of those endonucleases to accommodate that DNA sequence change, were also measured. The substitution of base-specific hydrogen bonds in a wild-type endonuclease/DNA complex with hydrophobic van der Waals contacts in a redesigned complex reduced the ability to discriminate between sites, due to nonspecific ΔS_binding.