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

Evolution of Restriction–Modification Systems Consisting of One Restriction Endonuclease and Two DNA Methyltransferases

Some restriction-modification systems contain two DNA methyltransferases. In the present work, we have classified such systems according to the families of catalytic domains present in the restriction endonucleases and both DNA methyltransferases. Evolution of the restriction-modification systems containing an endonuclease with a NOV_C family domain and two DNA methyltransferases, both with DNA_methylase family domains, was investigated in detail. Phylogenetic tree of DNA methyltransferases from the systems of this class consists of two clades of the same size. Two DNA methyltransferases of each restriction-modification system of this class belong to the different clades. This indicates independent evolution of the two methyltransferases. We detected multiple cross-species horizontal transfers of the systems as a whole, as well as the cases of gene transfer between the systems.

Toward Reducing Phylostratigraphic Errors and Biases

Phylostratigraphy is a method for estimating gene age, usually applied to large numbers of genes in order to detect nonrandom age-distributions of gene properties that could shed light on mechanisms of gene origination and evolution. However, phylostratigraphy underestimates gene age with a nonnegligible probability. The underestimation is severer for genes with certain properties, creating spurious age distributions of these properties and those correlated with these properties. Here we explore three strategies to reduce phylostratigraphic error/bias. First, we test several alternative homology detection methods (PSIBLAST, HMMER, PHMMER, OMA, and GLAM2Scan) in phylostratigraphy, but fail to find any that noticeably outperforms the commonly used BLASTP. Second, using machine learning, we look for predictors of error-prone genes to exclude from phylostratigraphy, but cannot identify reliable predictors. Finally, we remove from phylostratigraphic analysis genes exhibiting errors in simulation, which by definition minimizes error/bias if the simulation is sufficiently realistic. Using this last approach, we show that some previously reported phylostratigraphic trends (e.g., younger proteins tend to evolve more rapidly and be shorter) disappear or even reverse, reconfirming the necessity of controlling phylostratigraphic error/bias. Taken together, our analyses demonstrate that phylostratigraphic errors/biases are refractory to several potential solutions but can be controlled at least partially by the exclusion of error-prone genes identified via realistic simulations. These results are expected to stimulate the judicious use of error-aware phylostratigraphy and reevaluation of previous phylostratigraphic findings.

Natural zinc ribbon HNH endonucleases and engineered zinc finger nicking endonuclease

Many bacteriophage and prophage genomes encode an HNH endonuclease (HNHE) next to their cohesive end site and terminase genes. The HNH catalytic domain contains the conserved catalytic residues His-Asn-His and a zinc-binding site [CxxC]₂. An additional zinc ribbon (ZR) domain with one to two zinc-binding sites ([CxxxxC], [CxxxxH], [CxxxC], [HxxxH], [CxxC] or [CxxH]) is frequently found at the N-terminus or C-terminus of the HNHE or a ZR domain protein (ZRP) located adjacent to the HNHE. We expressed and purified 10 such HNHEs and characterized their cleavage sites. These HNHEs are site-specific and strand-specific nicking endonucleases (NEase or nickase) with 3- to 7-bp specificities. A minimal HNH nicking domain of 76 amino acid residues was identified from Bacillus phage γ HNHE and subsequently fused to a zinc finger protein to generate a chimeric NEase with a new specificity (12–13 bp). The identification of a large pool of previously unknown natural NEases and engineered NEases provides more ‘tools’ for DNA manipulation and molecular diagnostics. The small modular HNH nicking domain can be used to generate rare NEases applicable to targeted genome editing. In addition, the engineered ZF nickase is useful for evaluation of off-target sites in vitro before performing cell-based gene modification.

Chemically-induced affinity star restriction specificity: a novel TspGWI/sinefungin endonuclease with theoretical 3-bp cleavage frequency

The type IIS/IIC restriction endonuclease TspGWI recognizes the sequence 5'-ACGGA-3', cleaving DNA 11/9 nucleotides downstream. Here we show that sinefungin, a cofactor analog of S-adenosyl methionine, induces a unique type of relaxation in DNA recognition specificity. In the presence of sinefungin, TspGWI recognizes and cleaves at least 12 degenerate variants of the original recognition sequence that vary by single base pair changes from the original 5-bp restriction site with only a single degeneracy per variant appearing to be allowed. In addition, sinefungin was found to have a stimulatory effect on cleavage at these nondegenerate TspGWI recognition sites, irrespective of their number or the DNA topology. Interestingly, no fixed "core" could be identified among the new recognition sequences. Theoretically, TspGWI cleaves DNA every 1024 bp, while sinefungin-induced activity cleaves every 78.8 bp, corresponding to a putative 3-bp long recognition site. Thus, the combination of sinefungin and TspGWI represents a novel frequent cutter, next only to CviJI/CviJI*, that should prove useful in DNA cloning methodologies.

S-adenosyl homocysteine and DNA ends stimulate promiscuous nuclease activities in the Type III restriction endonuclease EcoPI

In the absence of the methyl donor S-adenosyl methionine and under certain permissive reaction conditions, EcoPI shows non-specific endonuclease activity. We show here that the cofactor analogue S-adenosyl homocysteine promotes this promiscuous DNA cleavage. Additionally, an extensive exonuclease-like processing of the DNA is also observed that can even result in digestion of non-specific DNA in trans. We suggest a model for how DNA communication events initiating from non-specific sites, and in particular free DNA ends, could produce the observed cleavage patterns.

Domain organization and functional analysis of type IIS restriction endonuclease Eco31I

Type IIS restriction endonuclease Eco31I harbors a single HNH active site and cleaves both DNA strands close to its recognition sequence, 5'-GGTCTC(1/5). A two-domain organization of Eco31I was determined by limited proteolysis. Analysis of proteolytic fragments revealed that the N-terminal domain of Eco31I is responsible for the specific DNA binding, while the C-terminal domain contains the HNH nuclease-like active site. Gel-shift and gel-filtration experiments revealed that a monomer of the N-terminal domain of Eco31I is able to bind a single copy of cognate DNA. However, in contrast to other studied type IIS enzymes, the isolated catalytic domain of Eco31I was inactive. Steady-state and transient kinetic analysis of Eco31I reactions was inconsistent with dimerization of Eco31I on DNA. Thus, we propose that Eco31I interacts with individual copies of its recognition sequence in its monomeric form and presumably remains a monomer as it cleaves both strands of double-stranded DNA. The domain organization and reaction mechanism established for Eco31I should be common for a group of evolutionary related type IIS restriction endonucleases Alw26I, BsaI, BsmAI, BsmBI and Esp3I that recognize DNA sequences bearing the common pentanucleotide 5'-GTCTC.

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

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

The recognition domain of the BpuJI restriction endonuclease in complex with cognate DNA at 1.3-A resolution

Type IIS restriction endonucleases recognize asymmetric DNA sequences and cleave both DNA strands at fixed positions downstream of the recognition site. The restriction endonuclease BpuJI recognizes the asymmetric sequence 5'-CCCGT; however, it cuts at multiple sites in the vicinity of the target sequence. BpuJI consists of two physically separate domains, with catalytic and dimerization functions in the C-terminal domain and DNA recognition functions in the N-terminal domain. Here we report the crystal structure of the BpuJI recognition domain bound to cognate DNA at 1.3-A resolution. This region folds into two winged-helix subdomains, D1 and D2, interspaced by the DL subdomain. The D1 and D2 subdomains of BpuJI share structural similarity with the similar subdomains of the FokI DNA-binding domain; however, their orientations in protein-DNA complexes are different. Recognition of the 5'-CCCGT target sequence is achieved by BpuJI through the major groove contacts of amino acid residues located on both the helix-turn-helix motifs and the N-terminal arm. The role of these interactions in DNA recognition is also corroborated by mutational analysis.

Tetrameric restriction enzymes: expansion to the GIY-YIG nuclease family

The GIY-YIG nuclease domain was originally identified in homing endonucleases and enzymes involved in DNA repair and recombination. Many of the GIY-YIG family enzymes are functional as monomers. We show here that the Cfr42I restriction endonuclease which belongs to the GIY-YIG family and recognizes the symmetric sequence 5′-CCGC/GG-3′ (‘/’ indicates the cleavage site) is a tetramer in solution. Moreover, biochemical and kinetic studies provided here demonstrate that the Cfr42I tetramer is catalytically active only upon simultaneous binding of two copies of its recognition sequence. In that respect Cfr42I resembles the homotetrameric Type IIF restriction enzymes that belong to the distinct PD-(E/D)XK nuclease superfamily. Unlike the PD-(E/D)XK enzymes, the GIY-YIG nuclease Cfr42I accommodates an extremely wide selection of metal-ion cofactors, including Mg²⁺, Mn²⁺, Co²⁺, Zn²⁺, Ni²⁺, Cu²⁺ and Ca²⁺. To our knowledge, Cfr42I is the first tetrameric GIY-YIG family enzyme. Similar structural arrangement and phenotypes displayed by restriction enzymes of the PD-(E/D)XK and GIY-YIG nuclease families point to the functional significance of tetramerization.