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

Characterization of winged helix domain fusion endonucleases as N6-methyladenine-dependent type IV restriction systems

Winged helix (wH) domains, also termed winged helix-turn-helix (wHTH) domains, are widespread in all kingdoms of life and have diverse roles. In the context of DNA binding and DNA modification sensing, some eukaryotic wH domains are known as sensors of non-methylated CpG. In contrast, the prokaryotic wH domains in DpnI and HhiV4I act as sensors of adenine methylation in the 6mApT (N6-methyladenine, 6mA, or N6mA) context. DNA-binding modes and interactions with the probed dinucleotide are vastly different in the two cases. Here, we show that the role of the wH domain as a sensor of adenine methylation is widespread in prokaryotes. We present previously uncharacterized examples of PD-(D/E)XK-wH (FcyTI, Psp4BI), PUA-wH-HNH (HtuIII), wH-GIY-YIG (Ahi29725I, Apa233I), and PLD-wH (Aba4572I, CbaI) fusion endonucleases that sense adenine methylation in the Dam+ Gm6ATC sequence contexts. Representatives of the wH domain endonuclease fusion families with the exception of the PLD-wH family could be purified, and an in vitro preference for adenine methylation in the Dam context could be demonstrated. Like most other modification-dependent restriction endonucleases (MDREs, also called type IV restriction systems), the new fusion endonucleases except those in the PD-(D/E)XK-wH family cleave close to but outside the recognition sequence. Taken together, our data illustrate the widespread combinatorial use of prokaryotic wH domains as adenine methylation readers. Other potential 6mA sensors in modified DNA are also discussed.

Orientation-Guided Immobilization of Probe DNA on swCNT-FET for Enhancing Sensitivity of EcoRV Detection

We present a novel approach that integrates electrical measurements with molecular dynamics (MD) simulations to assess the activity of type-II restriction endonucleases, specifically EcoRV. Our approach employs a single-walled carbon nanotube field-effect transistor (swCNT-FET) functionalized with the EcoRV substrate DNA, enabling the detection of enzymatic cleavage events. Notably, we leveraged the methylene blue (MB) tag as an "orientation guide" to immobilize the EcoRV substrate DNA in a specific direction, thereby enhancing the proximity of the DNA cleavage reaction to the swCNT surface and consequently improving the sensitivity in EcoRV detection. We conducted computational modeling to compare the conformations and electrostatic potential (ESP) of MB-tagged DNA with its MB-free counterpart, providing strong support for our electrical measurements. Both conformational and ESP simulations exhibited robust agreement with our experimental data. The inhibitory efficacy of the EcoRV inhibitor aurintricarboxylic acid (ATA) was also evaluated, and the selectivity of the sensing device was examined.

The crystal structures of Sau3AI with and without bound DNA suggest a self-activation-based DNA cleavage mechanism

The type II restriction endonuclease Sau3AI cleaves the sequence 5'-GATC-3' in double-strand DNA producing two sticky ends. Sau3AI cuts both DNA strands regardless of methylation status. Here, we report the crystal structures of the active site mutant Sau3AI-E64A and the C-terminal domain Sau3AI-C with a bound GATC substrate. Interestingly, the catalytic site of the N-terminal domain (Sau3AI-N) is spatially blocked by the C-terminal domain, suggesting a potential self-inhibition of the enzyme. Interruption of Sau3AI-C binding to substrate DNA disrupts Sau3AI function, suggesting a functional linkage between the N- and C-terminal domains. We propose that Sau3AI-C behaves as an allosteric effector binding one GATC substrate, which triggers a conformational change to open the N-terminal catalytic site, resulting in the subsequent GATC recognition by Sau3AI-N and cleavage of the second GATC site. Our data indicate that Sau3AI and UbaLAI might represent a new subclass of type IIE restriction enzymes.

Mitochondrial DNA breaks activate an integrated stress response to reestablish homeostasis

Mitochondrial DNA double-strand breaks (mtDSBs) lead to the degradation of circular genomes and a reduction in copy number; yet, the cellular response in human cells remains elusive. Here, using mitochondrial-targeted restriction enzymes, we show that a subset of cells with mtDSBs exhibited defective mitochondrial protein import, reduced respiratory complexes, and loss of membrane potential. Electron microscopy confirmed the altered mitochondrial membrane and cristae ultrastructure. Intriguingly, mtDSBs triggered the integrated stress response (ISR) via the phosphorylation of eukaryotic translation initiation factor 2α (eIF2α) by DELE1 and heme-regulated eIF2α kinase (HRI). When ISR was inhibited, the cells experienced intensified mitochondrial defects and slower mtDNA recovery post-breakage. Lastly, through proteomics, we identified ATAD3A-a membrane-bound protein interacting with nucleoids-as potentially pivotal in relaying signals from impaired genomes to the inner mitochondrial membrane. In summary, our study delineates the cascade connecting damaged mitochondrial genomes to the cytoplasm and highlights the significance of the ISR in maintaining mitochondrial homeostasis amid genome instability.

Identification and characterization of a new HNH restriction endonuclease with unusual properties

Restriction-modification (R-M) systems form a large superfamily constituting bacterial innate immunity mechanism. The restriction endonucleases (REases) are very diverse in subunit structure, DNA recognition, co-factor requirement, and mechanism of action. Among the different catalytic motifs, HNH active sites containing REases are the second largest group distinguished by the presence of the ββα-metal finger fold. KpnI is the first member of the HNH-family REases whose homologs are present in many bacteria of Enterobacteriaceae having varied degrees of sequence similarity between them. Considering that the homologs with a high similarity may have retained KpnI-like properties, while those with a low similarity could be different, we have characterized a distant KpnI homolog present in a pathogenic Klebsiella pneumoniae NTUH K2044. A comparison of the properties of KpnI and KpnK revealed that despite their similarity and the HNH motif, these two enzymes have different properties viz oligomerization, cleavage pattern, metal ion requirement, recognition sequence, and sequence specificity. Unlike KpnI, KpnK is a monomer in solution, nicks double-stranded DNA, recognizes degenerate sequence, and catalyses the degradation of DNA into smaller products after the initial cleavage at preferred sites. Due to several distinctive properties, it can be classified as a variant of the Type IIS enzyme having nicking endonuclease activity. KEY POINTS: • KpnK is a distant homolog of KpnI and belongs to the ββα-metal finger superfamily. • Both KpnI and KpnK have widespread occurrence in K. pneumoniae strains. • KpnK is a Type IIS restriction endonuclease with a single-strand nicking property.

Functional characterization of a hypothetical protein (TTHA1873) from Thermus thermophilus

Thermus thermophilus is an extremely thermophilic organism that thrives at a temperature of 65°C. T. thermophilus genome has ~2218 genes, out of which 66% (1482 genes) have been annotated, and the remaining 34% (736 genes) are assigned as hypothetical proteins. In this work, biochemical and biophysical experiments were performed to characterize the hypothetical protein TTHA1873 from T. thermophilus. The hypothetical protein TTHA1873 acts as a nuclease, which indiscreetly cuts methylated and non-methylated DNA in divalent metal ions and relaxes the plasmid DNA in the presence of ATP. The chelation of metal ions with EDTA inhibits its activity. These results suggest that the hypothetical protein TTHA1873 would be a CRISPR-associated protein with non-specific DNase activity and ATP-dependent DNA-relaxing activity.

Structures, activity and mechanism of the Type IIS restriction endonuclease PaqCI

Type IIS restriction endonucleases contain separate DNA recognition and catalytic domains and cleave their substrates at well-defined distances outside their target sequences. They are employed in biotechnology for a variety of purposes, including the creation of gene-targeting zinc finger and TAL effector nucleases and DNA synthesis applications such as Golden Gate assembly. The most thoroughly studied Type IIS enzyme, FokI, has been shown to require multimerization and engagement with multiple DNA targets for optimal cleavage activity; however, details of how it or similar enzymes forms a DNA-bound reaction complex have not been described at atomic resolution. Here we describe biochemical analyses of DNA cleavage by the Type IIS PaqCI restriction endonuclease and a series of molecular structures in the presence and absence of multiple bound DNA targets. The enzyme displays a similar tetrameric organization of target recognition domains in the absence or presence of bound substrate, with a significant repositioning of endonuclease domains in a trapped DNA-bound complex that is poised to deliver the first of a series of double-strand breaks. PaqCI and FokI share similar structural mechanisms of DNA cleavage, but considerable differences in their domain organization and quaternary architecture, facilitating comparisons between distinct Type IIS enzymes.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

The Pet127 protein is a mitochondrial 5'-to-3' exoribonuclease from the PD-(D/E)XK superfamily involved in RNA maturation and intron degradation in yeasts

Pet127 is a mitochondrial protein found in multiple eukaryotic lineages, but absent from several taxa, including plants and animals. Distant homology suggests that it belongs to the divergent PD-(D/E)XK superfamily which includes various nucleases and related proteins. Earlier yeast genetics experiments suggest that it plays a nonessential role in RNA degradation and 5' end processing. Our phylogenetic analysis suggests that it is a primordial eukaryotic invention that was retained in diverse groups, and independently lost several times in the evolution of other organisms. We demonstrate for the first time that the fungal Pet127 protein in vitro is a processive 5'-to-3' exoribonuclease capable of digesting various substrates in a sequence nonspecific manner. Mutations in conserved residues essential in the PD-(D/E)XK superfamily active site abolish the activity of Pet127. Deletion of the PET127 gene in the pathogenic yeast Candida albicans results in a moderate increase in the steady-state levels of several transcripts and in accumulation of unspliced precursors and intronic sequences of three introns. Mutations in the active site residues result in a phenotype identical to that of the deletant, confirming that the exoribonuclease activity is related to the physiological role of the Pet127 protein. Pet127 activity is, however, not essential for maintaining the mitochondrial respiratory activity in C. albicans.

Beyond Restriction Modification: Epigenomic Roles of DNA Methylation in Prokaryotes

The amount of bacterial and archaeal genome sequence and methylome data has greatly increased over the last decade, enabling new insights into the functional roles of DNA methylation in these organisms. Methyltransferases (MTases), the enzymes responsible for DNA methylation, are exchanged between prokaryotes through horizontal gene transfer and can function either as part of restriction-modification systems or in apparent isolation as single (orphan) genes. The patterns of DNA methylation they confer on the host chromosome can have significant effects on gene expression, DNA replication, and other cellular processes. Some processes require very stable patterns of methylation, resulting in conservation of persistent MTases in a particular lineage. Other processes require patterns that are more dynamic yet more predictable than what is afforded by horizontal gene transfer and gene loss, resulting in phase-variable or recombination-driven MTase alleles. In this review, we discuss what is currently known about the functions of DNA methylation in prokaryotes in light of these evolutionary patterns. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Restriction endonuclease T.Smu451I with new cleavage specificity-neoschizomer of T.AsuI

A specific type II restriction endonuclease T.Smu451I has been purified to electrophoretic homogeneity from the frozen cells of soil bacterium Sphingobacterium multivorum 451 (formerly Flavobacterium multivorum 451), using ultrasonic grinding, nucleic acid removal by streptomycin sulfate, protein precipitation by ammonium sulfate and phosphocellulose P-11, DEAE-Cellulose DE-52, Hepharin-Sepharose CL-6B chromatography, and elucidated several characteristics of T.Smu451I. The molecular weight of the enzyme determined by gel filtration and SDS-polyacrylamide gel electrophoresis was calculated to be 45,000 ± 2000 D (dimer) and 23,000 ± 1000 D (monomer), respectively. The isoelectric point (pI) of T.Smu451I is 5.4. T.Smu451I recognizes pentanucleotide palindromic sequences 5'-GGNC↓C-3' and cleaves between C and C in position shown by arrow to produce 3'-cohesive terminus of trinucleotide. Therefore, T.Smu451I is a neoschizomer of T.AsuI.

Crystal structure and DNA cleavage mechanism of the restriction DNA glycosylase R.CcoLI from Campylobacter coli

While most restriction enzymes catalyze the hydrolysis of phosphodiester bonds at specific nucleotide sequences in DNA, restriction enzymes of the HALFPIPE superfamily cleave N-glycosidic bonds, similar to DNA glycosylases. Apurinic/apyrimidinic (AP) sites generated by HALFPIPE superfamily proteins are cleaved by their inherent AP lyase activities, other AP endonuclease activities or heat-promoted β-elimination. Although the HALFPIPE superfamily protein R.PabI, obtained from a hyperthermophilic archaea, Pyrococcus abyssi, shows weak AP lyase activity, HALFPIPE superfamily proteins in mesophiles, such as R.CcoLI from Campylobacter coli and R. HpyAXII from Helicobacter pylori, show significant AP lyase activities. To identify the structural basis for the AP lyase activity of R.CcoLI, we determined the structure of R.CcoLI by X-ray crystallography. The structure of R.CcoLI, obtained at 2.35-Å resolution, shows that a conserved lysine residue (Lys71), which is stabilized by a characteristic β-sheet structure of R.CcoLI, protrudes into the active site. The results of mutational assays indicate that Lys71 is important for the AP lyase activity of R.CcoLI. Our results help to elucidate the mechanism by which HALFPIPE superfamily proteins from mesophiles efficiently introduce double-strand breaks to specific sites on double-stranded DNA.

Restriction endonucleases that cleave RNA/DNA heteroduplexes bind dsDNA in A-like conformation

Restriction endonucleases naturally target DNA duplexes. Systematic screening has identified a small minority of these enzymes that can also cleave RNA/DNA heteroduplexes and that may therefore be useful as tools for RNA biochemistry. We have chosen AvaII (G↓GWCC, where W stands for A or T) as a representative of this group of restriction endonucleases for detailed characterization. Here, we report crystal structures of AvaII alone, in specific complex with partially cleaved dsDNA, and in scanning complex with an RNA/DNA hybrid. The specific complex reveals a novel form of semi-specific dsDNA readout by a hexa-coordinated metal cation, most likely Ca²⁺ or Mg²⁺. Substitutions of residues anchoring this non-catalytic metal ion severely impair DNA binding and cleavage. The dsDNA in the AvaII complex is in the A-like form. This creates space for 2′-OH groups to be accommodated without intra-nucleic acid steric conflicts. PD-(D/E)XK restriction endonucleases of known structure that bind their dsDNA targets in the A-like form cluster into structurally similar groups. Most such enzymes, including some not previously studied in this respect, cleave RNA/DNA heteroduplexes. We conclude that A-form dsDNA binding is a good predictor for RNA/DNA cleavage activity.

The key role of E418 carboxyl group in the formation of Nt.BspD6I nickase active site: Structural and functional properties of Nt.BspD6I E418A mutant

The mutated nickase Nt.BspD6I E418A has been obtained by site-directed mutagenesis. The purified protein has been crystallized, and its spatial structure has been determined at 2.45 Å resolution. An analysis of the crystal structures of the wild-type and mutated nickase have shown that the elimination of a carboxyl group due to the E418A mutation initiates marked conformational changes in both the N-terminal recognition domain and the C-terminal catalytic domain of nickase and insignificantly affects its linker domain. This is supported by changes in the functional properties of mutated nickase: an increase in the oligomerization capacity in the presence of a substrate, a reduction in the capacity to bind a substrate, and complete loss of catalytic activity.

Distortion of double-stranded DNA structure by the binding of the restriction DNA glycosylase R.PabI

R.PabI is a restriction DNA glycosylase that recognizes the sequence 5′-GTAC-3′ and hydrolyses the N-glycosidic bond of adenine in the recognition sequence. R.PabI drastically bends and unwinds the recognition sequence of double-stranded DNA (dsDNA) and flips the adenine and guanine bases in the recognition sequence into the catalytic and recognition sites on the protein surface. In this study, we determined the crystal structure of the R.PabI-dsDNA complex in which the dsDNA is drastically bent by the binding of R.PabI but the base pairs are not unwound. This structure is predicted to be important for the indirect readout of the recognition sequence by R.PabI. In the complex structure, wedge loops of the R.PabI dimer are inserted into the minor groove of dsDNA to stabilize the deformed dsDNA structure. A base stacking is distorted between the two wedge-inserted regions. R.PabI is predicted to utilize the distorted base stacking for the detection of the recognition sequence.

Functional screening of a soil metagenome for DNA endonucleases by acquired resistance to bacteriophage infection

Endonucleases play a crucial role as reagents in laboratory research and diagnostics. Here, metagenomics was used to functionally screen a fosmid library for endonucleases. A fosmid library was constructed using metagenomic DNA isolated from soil sampled from the unique environment of the Kogelberg Nature Reserve in the Western Cape of South Africa. The principle of acquired immunity against phage infection was used to develop a plate-based screening technique for the isolation of restriction endonucleases from the library. Using next-generation sequencing and bioinformatics tools, sequence data were generated and analysed, revealing 113 novel open reading frames (ORFs) encoding putative endonuclease genes and ORFs of unknown identity and function. One endonuclease designated Endo52 was selected from the putative endonuclease ORFs and was recombinantly produced in Escherichia coli Rosetta™ (DE3) pLysS. Endo52 was purified by immobilised metal affinity chromatography and yielded 0.437 g per litre of cultivation volume. Its enzyme activity was monitored by cleaving lambda DNA and pUC19 plasmid as substrates, and it demonstrated non-specific endonuclease activity. In addition to endonuclease-like genes, the screen identified several unknown genes. These could present new phage resistance mechanisms and are an opportunity for future investigations.

Evolutionary entanglement of mobile genetic elements and host defence systems: guns for hire

All cellular life forms are afflicted by diverse genetic parasites, including viruses and other types of mobile genetic elements (MGEs), and have evolved multiple, diverse defence systems that protect them from MGE assault via different mechanisms. Here, we provide our perspectives on how recent evidence points to tight evolutionary connections between MGEs and defence systems that reach far beyond the proverbial arms race. Defence systems incur a fitness cost for the hosts; therefore, at least in prokaryotes, horizontal mobility of defence systems, mediated primarily by MGEs, is essential for their persistence. Moreover, defence systems themselves possess certain features of selfish elements. Common components of MGEs, such as site-specific nucleases, are 'guns for hire' that can also function as parts of defence mechanisms and are often shuttled between MGEs and defence systems. Thus, evolutionary and molecular factors converge to mould the multifaceted, inextricable connection between MGEs and anti-MGE defence systems.

Role of the R349 Gene and Its Repeats in the MIMIVIRE Defense System

MIMIVIRE is a defense system described in lineage A Mimivirus (Mimiviridae family) that mediates resistance against Zamilon virophage. It is composed of putative helicase and nuclease associated with a gene of unknown function called R349, which contains four 15 bp repeats homologous to the virophage sequence. In a previous study, the silencing of such genes restored virophage susceptibility. Moreover, the CRISPR Cas-4 like activity of the nuclease has recently been characterized. In this study, a recently isolated Mimivirus of lineage A with R349 gene lacking 3 of 4 repeats was demonstrated to be susceptible to Zamilon. To reinforce the importance of the R349 gene in the MIMIVIRE system, we developed and presented, for the first time to our knowledge, a protocol for Mimivirus genomic editing. By knocking out R349 gene in a Mimivirus lineage A, we observed the replication of Zamilon, indicating that this gene is critical in the resistance against this specific group of virophages.

Structure, subunit organization and behavior of the asymmetric Type IIT restriction endonuclease BbvCI

BbvCI, a Type IIT restriction endonuclease, recognizes and cleaves the seven base pair sequence 5'-CCTCAGC-3', generating 3-base, 5'-overhangs. BbvCI is composed of two protein subunits, each containing one catalytic site. Either site can be inactivated by mutation resulting in enzyme variants that nick DNA in a strand-specific manner. Here we demonstrate that the holoenzyme is labile, with the R1 subunit dissociating at low pH. Crystallization of the R2 subunit under such conditions revealed an elongated dimer with the two catalytic sites located on opposite sides. Subsequent crystallization at physiological pH revealed a tetramer comprising two copies of each subunit, with a pair of deep clefts each containing two catalytic sites appropriately positioned and oriented for DNA cleavage. This domain organization was further validated with single-chain protein constructs in which the two enzyme subunits were tethered via peptide linkers of variable length. We were unable to crystallize a DNA-bound complex; however, structural similarity to previously crystallized restriction endonucleases facilitated creation of an energy-minimized model bound to DNA, and identification of candidate residues responsible for target recognition. Mutation of residues predicted to recognize the central C:G base pair resulted in an altered enzyme that recognizes and cleaves CCTNAGC (N = any base).

Systematic classification of the His-Me finger superfamily

The His-Me finger endonucleases, also known as HNH or ββα-metal endonucleases, form a large and diverse protein superfamily. The His-Me finger domain can be found in proteins that play an essential role in cells, including genome maintenance, intron homing, host defense and target offense. Its overall structural compactness and non-specificity make it a perfectly-tailored pathogenic module that participates on both sides of inter- and intra-organismal competition. An extremely low sequence similarity across the superfamily makes it difficult to identify and classify new His-Me fingers. Using state-of-the-art distant homology detection methods, we provide an updated and systematic classification of His-Me finger proteins. In this work, we identified over 100 000 proteins and clustered them into 38 groups, of which three groups are new and cannot be found in any existing public domain database of protein families. Based on an analysis of sequences, structures, domain architectures, and genomic contexts, we provide a careful functional annotation of the poorly characterized members of this superfamily. Our results may inspire further experimental investigations that should address the predicted activity and clarify the potential substrates, to provide more detailed insights into the fundamental biological roles of these proteins.

Evolutionary Genomics of Defense Systems in Archaea and Bacteria

Evolution of bacteria and archaea involves an incessant arms race against an enormous diversity of genetic parasites. Accordingly, a substantial fraction of the genes in most bacteria and archaea are dedicated to antiparasite defense. The functions of these defense systems follow several distinct strategies, including innate immunity; adaptive immunity; and dormancy induction, or programmed cell death. Recent comparative genomic studies taking advantage of the expanding database of microbial genomes and metagenomes, combined with direct experiments, resulted in the discovery of several previously unknown defense systems, including innate immunity centered on Argonaute proteins, bacteriophage exclusion, and new types of CRISPR-Cas systems of adaptive immunity. Some general principles of function and evolution of defense systems are starting to crystallize, in particular, extensive gain and loss of defense genes during the evolution of prokaryotes; formation of genomic defense islands; evolutionary connections between mobile genetic elements and defense, whereby genes of mobile elements are repeatedly recruited for defense functions; the partially selfish and addictive behavior of the defense systems; and coupling between immunity and dormancy induction/programmed cell death.

Alga-PrAS (Algal Protein Annotation Suite): A Database of Comprehensive Annotation in Algal Proteomes

Algae are smaller organisms than land plants and offer clear advantages in research over terrestrial species in terms of rapid production, short generation time and varied commercial applications. Thus, studies investigating the practical development of effective algal production are important and will improve our understanding of both aquatic and terrestrial plants. In this study we estimated multiple physicochemical and secondary structural properties of protein sequences, the predicted presence of post-translational modification (PTM) sites, and subcellular localization using a total of 510,123 protein sequences from the proteomes of 31 algal and three plant species. Algal species were broadly selected from green and red algae, glaucophytes, oomycetes, diatoms and other microalgal groups. The results were deposited in the Algal Protein Annotation Suite database (Alga-PrAS; http://alga-pras.riken.jp/), which can be freely accessed online.

Structure of the EndoMS-DNA Complex as Mismatch Restriction Endonuclease

Archaeal NucS nuclease was thought to degrade the single-stranded region of branched DNA, which contains flapped and splayed DNA. However, recent findings indicated that EndoMS, the orthologous enzyme of NucS, specifically cleaves double-stranded DNA (dsDNA) containing mismatched bases. In this study, we determined the structure of the EndoMS-DNA complex. The complex structure of the EndoMS dimer with dsDNA unexpectedly revealed that the mismatched bases were flipped out into binding sites, and the overall architecture most resembled that of restriction enzymes. The structure of the apo form was similar to the reported structure of Pyrococcus abyssi NucS, indicating that movement of the C-terminal domain from the resting state was required for activity. In addition, a model of the EndoMS-PCNA-DNA complex was preliminarily verified with electron microscopy. The structures strongly support the idea that EndoMS acts in a mismatch repair pathway.

Tetrameric structure of the restriction DNA glycosylase R.PabI in complex with nonspecific double-stranded DNA

R.PabI is a type II restriction enzyme that recognizes the 5′-GTAC-3′ sequence and belongs to the HALFPIPE superfamily. Although most restriction enzymes cleave phosphodiester bonds at specific sites by hydrolysis, R.PabI flips the guanine and adenine bases of the recognition sequence out of the DNA helix and hydrolyzes the N-glycosidic bond of the flipped adenine in a similar manner to DNA glycosylases. In this study, we determined the structure of R.PabI in complex with double-stranded DNA without the R.PabI recognition sequence by X-ray crystallography. The 1.9 Å resolution structure of the complex showed that R.PabI forms a tetrameric structure to sandwich the double-stranded DNA and the tetrameric structure is stabilized by four salt bridges. DNA binding and DNA glycosylase assays of the R.PabI mutants showed that the residues that form the salt bridges (R70 and D71) are essential for R.PabI to find the recognition sequence from the sea of nonspecific sequences. R.PabI is predicted to utilize the tetrameric structure to bind nonspecific double-stranded DNA weakly and slide along it to find the recognition sequence.

Intrinsic protein disorder could be overlooked in cocrystallization conditions: An SRCD case study

X-ray diffractometry dominates protein studies, as it can provide 3D structures of these diverse macromolecules or their molecular complexes with interacting partners: substrates, inhibitors, and/or cofactors. Here, we show that under cocrystallization conditions the results could reflect induced protein folds instead of the (partially) disordered original structures. The analysis of synchrotron radiation circular dichroism spectra revealed that the Im7 immunity protein stabilizes the native-like solution structure of unfolded NColE7 nuclease mutants via complex formation. This is consistent with the fact that among the several available crystal structures with its inhibitor or substrate, all NColE7 structures are virtually the same. Our results draw attention to the possible structural consequence of protein modifications, which is often hidden by compensational effects of intermolecular interactions. The growing evidence on the importance of protein intrinsic disorder thus, demands more extensive complementary experiments in solution phase with the unligated form of the protein of interest.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

A sequence-specific DNA glycosylase mediates restriction-modification in Pyrococcus abyssi

Restriction-modification systems consist of genes that encode a restriction enzyme and a cognate methyltransferase. Thus far, it was believed that restriction enzymes are sequence-specific endonucleases that introduce double-strand breaks at specific sites by catalysing the cleavages of phosphodiester bonds. Here we report that based on the crystal structure and enzymatic activity, one of the restriction enzymes, R.PabI, is not an endonuclease but a sequence-specific adenine DNA glycosylase. The structure of the R.PabI-DNA complex shows that R.PabI unwinds DNA at a 5'-GTAC-3' site and flips the guanine and adenine bases out of the DNA helix to recognize the sequence. R.PabI catalyses the hydrolysis of the N-glycosidic bond between the adenine base and the sugar in the DNA and produces two opposing apurinic/apyrimidinic (AP) sites. The opposing AP sites are cleaved by heat-promoted β elimination and/or by endogenous AP endonucleases of host cells to introduce a double-strand break.

The conserved molecular machinery in DNA mismatch repair enzyme structures

The machinery of DNA mismatch repair enzymes is highly conserved in evolution. The process is initiated by recognition of a DNA mismatch, and validated by ATP and the presence of a processivity clamp or a methylation mark. Several events in MMR promote conformational changes that lead to progression of the repair process. Here we discuss functional conformational changes in the MMR proteins and we compare the enzymes to paralogs in other systems.

Transmission of the PabI family of restriction DNA glycosylase genes: mobility and long-term inheritance

Background

R.PabI is an exceptional restriction enzyme that functions as a DNA glycosylase. The enzyme excises an unmethylated base from its recognition sequence to generate apurinic/apyrimidinic (AP) sites, and also displays AP lyase activity, cleaving the DNA backbone at the AP site to generate the 3’-phospho alpha, beta-unsaturated aldehyde end in addition to the 5’-phosphate end. The resulting ends are difficult to religate with DNA ligase. The enzyme was originally isolated in Pyrococcus, a hyperthermophilic archaeon, and additional homologs subsequently identified in the epsilon class of the Gram-negative bacterial phylum Proteobacteria, such as Helicobacter pylori.

Results

Systematic analysis of R.PabI homologs and their neighboring genes in sequenced genomes revealed co-occurrence of R.PabI with M.PabI homolog methyltransferase genes. R.PabI and M.PabI homolog genes are occasionally found at corresponding (orthologous) loci in different species, such as Helicobacter pylori, Helicobacter acinonychis and Helicobacter cetorum, indicating long-term maintenance of the gene pair. One R.PabI and M.PabI homolog gene pair is observed immediately after the GMP synthase gene in both Campylobacter and Helicobacter, representing orthologs beyond genera. The mobility of the PabI family of restriction-modification (RM) system between genomes is evident upon comparison of genomes of sibling strains/species. Analysis of R.PabI and M.PabI homologs in H. pylori revealed an insertion of integrative and conjugative elements (ICE), and replacement with a gene of unknown function that may specify a membrane-associated toxin (hrgC). In view of the similarity of HrgC with toxins in type I toxin-antitoxin systems, we addressed the biological significance of this substitution. Our data indicate that replacement with hrgC occurred in the common ancestor of hspAmerind and hspEAsia. Subsequently, H. pylori with and without hrgC were intermixed at this locus, leading to complex distribution of hrgC in East Asia and the Americas. In Malaysia, hrgC was horizontally transferred from hspEAsia to hpAsia2 strains.

Conclusions

The PabI family of RM system behaves as a mobile, selfish genetic element, similar to the other families of Type II RM systems. Our analysis additionally revealed some cases of long-term inheritance. The distribution of the hrgC gene replacing the PabI family in the subpopulations of H. pylori, hspAmerind, hspEAsia and hpAsia2, corresponds to the two human migration events, one from East Asia to Americas and the other from China to Malaysia.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2021-3) contains supplementary material, which is available to authorized users.

Investigation into the annotation of protocol sequencing steps in the sequence read archive

Background

The workflow for the production of high-throughput sequencing data from nucleic acid samples is complex. There are a series of protocol steps to be followed in the preparation of samples for next-generation sequencing. The quantification of bias in a number of protocol steps, namely DNA fractionation, blunting, phosphorylation, adapter ligation and library enrichment, remains to be determined.

Results

We examined the experimental metadata of the public repository Sequence Read Archive (SRA) in order to ascertain the level of annotation of important sequencing steps in submissions to the database. Using SQL relational database queries (using the SRAdb SQLite database generated by the Bioconductor consortium) to search for keywords commonly occurring in key preparatory protocol steps partitioned over studies, we found that 7.10%, 5.84% and 7.57% of all records (fragmentation, ligation and enrichment, respectively), had at least one keyword corresponding to one of the three protocol steps. Only 4.06% of all records, partitioned over studies, had keywords for all three steps in the protocol (5.58% of all SRA records).

Conclusions

The current level of annotation in the SRA inhibits systematic studies of bias due to these protocol steps. Downstream from this, meta-analyses and comparative studies based on these data will have a source of bias that cannot be quantified at present.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-015-0064-7) contains supplementary material, which is available to authorized users.

The SfaNI restriction-modification system from Enterococcus faecalis NEB215 is located on a putative mobile genetic element

A type IIS restriction-modification (R-M) system SfaNI from Enterococcus faecalis NEB215 has been characterized. The sfaNIM gene was cloned by the methylase selection method. Methyltransferase SfaNI, a protein of 695 amino acids, consists of two domains responsible for different DNA-strand recognition and modification, and a putative DNA-binding HTH domain located in the N-terminal part of the protein. The sfaNIR gene, located adjacent to the gene of the cognate modification methyltransferases, encodes a protein of 648 amino acids. The enzyme has been purified to apparent homogeneity and its biochemical characteristics have been described. The R-M system SfaNI is flanked by a transposase gene at its 5(') end, and a cassette chromosome recombinase (ccr) gene complex, encoding serine recombinases CcrA and CcrB, at the 3(') end. Both proteins are specifically involved in genome rearrangement and are widely distributed among staphylococcal species. These results suggested that the R-M system SfaNI is present on the putative mobile element.

Modification-dependent restriction endonuclease, MspJI, flips 5-methylcytosine out of the DNA helix

MspJI belongs to a family of restriction enzymes that cleave DNA containing 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC). MspJI is specific for the sequence 5(h)mC-N-N-G or A and cleaves with some variability 9/13 nucleotides downstream. Earlier, we reported the crystal structure of MspJI without DNA and proposed how it might recognize this sequence and catalyze cleavage. Here we report its co-crystal structure with a 27-base pair oligonucleotide containing 5mC. This structure confirms that MspJI acts as a homotetramer and that the modified cytosine is flipped from the DNA helix into an SRA-like-binding pocket. We expected the structure to reveal two DNA molecules bound specifically to the tetramer and engaged with the enzyme's two DNA-cleavage sites. A coincidence of crystal packing precluded this organization, however. We found that each DNA molecule interacted with two adjacent tetramers, binding one specifically and the other non-specifically. The latter interaction, which prevented cleavage-site engagement, also involved base flipping and might represent the sequence-interrogation phase that precedes specific recognition. MspJI is unusual in that DNA molecules are recognized and cleaved by different subunits. Such interchange of function might explain how other complex multimeric restriction enzymes act.

Type II restriction endonucleases--a historical perspective and more

This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss 'Type II' REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.

Design of a colicin E7 based chimeric zinc-finger nuclease

Colicin E7 is a natural bacterial toxin. Its nuclease domain (NColE7) enters the target cell and kills it by digesting the nucleic acids. The HNH-motif as the catalytic centre of NColE7 at the C-terminus requires the positively charged N-terminal loop for the nuclease activity—offering opportunities for allosteric control in a NColE7-based artificial nuclease. Accordingly, four novel zinc finger nucleases were designed by computational methods exploiting the special structural features of NColE7. The constructed models were subjected to MD simulations. The comparison of structural stability and functional aspects showed that these models may function as safely controlled artificial nucleases. This study was complemented by random mutagenesis experiments identifying potentially important residues for NColE7 function outside the catalytic region.

Fine tuning of the catalytic activity of colicin E7 nuclease domain by systematic N-terminal mutations

The nuclease domain of colicin E7 (NColE7) promotes the nonspecific cleavage of nucleic acids at its C-terminal HNH motif. Interestingly, the deletion of four N-terminal residues (446-449 NColE7 = KRNK) resulted in complete loss of the enzyme activity. R447A mutation was reported to decrease the nuclease activity, but a detailed analysis of the role of the highly positive and flexible N-terminus is still missing. Here, we present the study of four mutants, with a decreased activity in the following order: NColE7  >> KGNK > KGNG ∼ GGNK > GGNG. At the same time, the folding, the metal-ion, and the DNA-binding affinity were unaffected by the mutations as revealed by linear and circular dichroism spectroscopy, isothermal calorimetric titrations, and gel mobility shift experiments. Semiempirical quantum chemical calculations and molecular dynamics simulations revealed that K446, K449, and/or the N-terminal amino group are able to approach the active centre in the absence of the other positively charged residues. The results suggested a complex role of the N-terminus in the catalytic process that could be exploited in the design of a controlled nuclease.

Elimination of inter-domain interactions increases the cleavage fidelity of the restriction endonuclease DraIII

DraIII is a type IIP restriction endonucleases (REases) that recognizes and creates a double strand break within the gapped palindromic sequence CAC↑NNN↓GTG of double-stranded DNA (↑ indicates nicking on the bottom strand; ↓ indicates nicking on the top strand). However, wild type DraIII shows significant star activity. In this study, it was found that the prominent star site is CAT↑GTT↓GTG, consisting of a star 5′ half (CAT) and a canonical 3′ half (GTG). DraIII nicks the 3′ canonical half site at a faster rate than the 5′ star half site, in contrast to the similar rate with the canonical full site. The crystal structure of the DraIII protein was solved. It indicated, as supported by mutagenesis, that DraIII possesses a ββα-metal HNH active site. The structure revealed extensive intra-molecular interactions between the N-terminal domain and the C-terminal domain containing the HNH active site. Disruptions of these interactions through site-directed mutagenesis drastically increased cleavage fidelity. The understanding of fidelity mechanisms will enable generation of high fidelity REases.

Structure and mutagenesis of the DNA modification-dependent restriction endonuclease AspBHI

The modification-dependent restriction endonuclease AspBHI recognizes 5-methylcytosine (5mC) in the double-strand DNA sequence context of (C/T)(C/G)(5mC)N(C/G) (N = any nucleotide) and cleaves the two strands a fixed distance (N₁₂/N₁₆) 3′ to the modified cytosine. We determined the crystal structure of the homo-tetrameric AspBHI. Each subunit of the protein comprises two domains: an N-terminal DNA-recognition domain and a C-terminal DNA cleavage domain. The N-terminal domain is structurally similar to the eukaryotic SET and RING-associated (SRA) domain, which is known to bind to a hemi-methylated CpG dinucleotide. The C-terminal domain is structurally similar to classic Type II restriction enzymes and contains the endonuclease catalytic-site motif of DX₂₀EAK. To understand how specific amino acids affect AspBHI recognition preference, we generated a homology model of the AspBHI-DNA complex, and probed the importance of individual amino acids by mutagenesis. Ser41 and Arg42 are predicted to be located in the DNA minor groove 5′ to the modified cytosine. Substitution of Ser41 with alanine (S41A) and cysteine (S41C) resulted in mutants with altered cleavage activity. All 19 Arg42 variants resulted in loss of endonuclease activity.

Type I restriction enzymes and their relatives

Type I restriction enzymes (REases) are large pentameric proteins with separate restriction (R), methylation (M) and DNA sequence-recognition (S) subunits. They were the first REases to be discovered and purified, but unlike the enormously useful Type II REases, they have yet to find a place in the enzymatic toolbox of molecular biologists. Type I enzymes have been difficult to characterize, but this is changing as genome analysis reveals their genes, and methylome analysis reveals their recognition sequences. Several Type I REases have been studied in detail and what has been learned about them invites greater attention. In this article, we discuss aspects of the biochemistry, biology and regulation of Type I REases, and of the mechanisms that bacteriophages and plasmids have evolved to evade them. Type I REases have a remarkable ability to change sequence specificity by domain shuffling and rearrangements. We summarize the classic experiments and observations that led to this discovery, and we discuss how this ability depends on the modular organizations of the enzymes and of their S subunits. Finally, we describe examples of Type II restriction-modification systems that have features in common with Type I enzymes, with emphasis on the varied Type IIG enzymes.

Highlights of the DNA cutters: a short history of the restriction enzymes

In the early 1950's, 'host-controlled variation in bacterial viruses' was reported as a non-hereditary phenomenon: one cycle of viral growth on certain bacterial hosts affected the ability of progeny virus to grow on other hosts by either restricting or enlarging their host range. Unlike mutation, this change was reversible, and one cycle of growth in the previous host returned the virus to its original form. These simple observations heralded the discovery of the endonuclease and methyltransferase activities of what are now termed Type I, II, III and IV DNA restriction-modification systems. The Type II restriction enzymes (e.g. EcoRI) gave rise to recombinant DNA technology that has transformed molecular biology and medicine. This review traces the discovery of restriction enzymes and their continuing impact on molecular biology and medicine.

Increasing cleavage specificity and activity of restriction endonuclease KpnI

Restriction enzyme KpnI is a HNH superfamily endonuclease requiring divalent metal ions for DNA cleavage but not for binding. The active site of KpnI can accommodate metal ions of different atomic radii for DNA cleavage. Although Mg²⁺ ion higher than 500 μM mediates promiscuous activity, Ca²⁺ suppresses the promiscuity and induces high cleavage fidelity. Here, we report that a conservative mutation of the metal-coordinating residue D148 to Glu results in the elimination of the Ca²⁺-mediated cleavage but imparting high cleavage fidelity with Mg²⁺. High cleavage fidelity of the mutant D148E is achieved through better discrimination of the target site at the binding and cleavage steps. Biochemical experiments and molecular dynamics simulations suggest that the mutation inhibits Ca²⁺-mediated cleavage activity by altering the geometry of the Ca²⁺-bound HNH active site. Although the D148E mutant reduces the specific activity of the enzyme, we identified a suppressor mutation that increases the turnover rate to restore the specific activity of the high fidelity mutant to the wild-type level. Our results show that active site plasticity in coordinating different metal ions is related to KpnI promiscuous activity, and tinkering the metal ion coordination is a plausible way to reduce promiscuous activity of metalloenzymes.

To be or not to be: regulation of restriction-modification systems and other toxin-antitoxin systems

One of the simplest classes of genes involved in programmed death is that containing the toxin–antitoxin (TA) systems of prokaryotes. These systems are composed of an intracellular toxin and an antitoxin that neutralizes its effect. These systems, now classified into five types, were initially discovered because some of them allow the stable maintenance of mobile genetic elements in a microbial population through postsegregational killing or the death of cells that have lost these systems. Here, we demonstrate parallels between some TA systems and restriction–modification systems (RM systems). RM systems are composed of a restriction enzyme (toxin) and a modification enzyme (antitoxin) and limit the genetic flux between lineages with different epigenetic identities, as defined by sequence-specific DNA methylation. The similarities between these systems include their postsegregational killing and their effects on global gene expression. Both require the finely regulated expression of a toxin and antitoxin. The antitoxin (modification enzyme) or linked protein may act as a transcriptional regulator. A regulatory antisense RNA recently identified in an RM system can be compared with those RNAs in TA systems. This review is intended to generalize the concept of TA systems in studies of stress responses, programmed death, genetic conflict and epigenetics.

DNA-binding domains of plant-specific transcription factors: structure, function, and evolution

The families of the plant-specific transcription factors (TFs) are defined by their characteristic DNA-binding domains (DBDs), such as AP2/ERF, B3, NAC, SBP, and WRKY. Recently, three-dimensional structures of the DBDs, including those in complexes with DNA, were determined by NMR spectroscopy and X-ray crystallography. In this review we summarize the functional and evolutionary implications arising from structure analyses. The unexpected structural similarity between B3 and the noncatalytic DBD of the restriction endonuclease EcoRII allowed us to build structural models of the B3/DNA complex. Most of the DBDs of plant-specific TFs are likely to have originated from endonucleases associated with transposable elements. After the DBDs have been established in unicellular eukaryotes, they experienced extensive plant-specific expansion, by acquiring new functions.

Massively parallel characterization of restriction endonucleases

Restriction endonucleases are highly specific in recognizing the particular DNA sequence they act on. However, their activity is affected by sequence context, enzyme concentration and buffer composition. Changes in these factors may lead to either ineffective cleavage at the cognate restriction site or relaxed specificity allowing cleavage of degenerate ‘star’ sites. Additionally, uncharacterized restriction endonucleases and engineered variants present novel activities. Traditionally, restriction endonuclease activity is assayed on simple substrates such as plasmids and synthesized oligonucleotides. We present and use high-throughput Illumina sequencing-based strategies to assay the sequence specificity and flanking sequence preference of restriction endonucleases. The techniques use fragmented DNA from sequenced genomes to quantify restriction endonuclease cleavage on a complex genomic DNA substrate in a single reaction. By mapping millions of restriction site–flanking reads back to the Escherichia coli and Drosophila melanogaster genomes we were able to quantitatively characterize the cognate and star site activity of EcoRI and MfeI and demonstrate genome-wide decreases in star activity with engineered high-fidelity variants EcoRI-HF and MfeI-HF, as well as quantify the influence on MfeI cleavage conferred by flanking nucleotides. The methods presented are readily applicable to all type II restriction endonucleases that cleave both strands of double-stranded DNA.

Comparative genomics of defense systems in archaea and bacteria

Our knowledge of prokaryotic defense systems has vastly expanded as the result of comparative genomic analysis, followed by experimental validation. This expansion is both quantitative, including the discovery of diverse new examples of known types of defense systems, such as restriction-modification or toxin-antitoxin systems, and qualitative, including the discovery of fundamentally new defense mechanisms, such as the CRISPR-Cas immunity system. Large-scale statistical analysis reveals that the distribution of different defense systems in bacterial and archaeal taxa is non-uniform, with four groups of organisms distinguishable with respect to the overall abundance and the balance between specific types of defense systems. The genes encoding defense system components in bacterial and archaea typically cluster in defense islands. In addition to genes encoding known defense systems, these islands contain numerous uncharacterized genes, which are candidates for new types of defense systems. The tight association of the genes encoding immunity systems and dormancy- or cell death-inducing defense systems in prokaryotic genomes suggests that these two major types of defense are functionally coupled, providing for effective protection at the population level.

Identification of a novel human mitochondrial endo-/exonuclease Ddk1/c20orf72 necessary for maintenance of proper 7S DNA levels

Although the human mitochondrial genome has been investigated for several decades, the proteins responsible for its replication and expression, especially nucleolytic enzymes, are poorly described. Here, we characterized a novel putative PD-(D/E)XK nuclease encoded by the human C20orf72 gene named Ddk1 for its predicted catalytic residues. We show that Ddk1 is a mitochondrially localized metal-dependent DNase lacking detectable ribonuclease activity. Ddk1 degrades DNA mainly in a 3'-5' direction with a strong preference for single-stranded DNA. Interestingly, Ddk1 requires free ends for its activity and does not degrade circular substrates. In addition, when a chimeric RNA-DNA substrate is provided, Ddk1 can slide over the RNA fragment and digest DNA endonucleolytically. Although the levels of the mitochondrial DNA are unchanged on RNAi-mediated depletion of Ddk1, the mitochondrial single-stranded DNA molecule (7S DNA) accumulates. On the other hand, overexperssion of Ddk1 decreases the levels of 7S DNA, suggesting an important role of the protein in 7S DNA regulation. We propose a structural model of Ddk1 and discuss its similarity to other PD-(D/E)XK superfamily members.