High-throughput sequencing of EcoWI restriction fragments maps the genome-wide landscape of phosphorothioate modification at base resolution.

Phosphorothioation (PT), in which a non-bridging oxygen is replaced by a sulfur, is one of the rare modifications discovered in bacteria and archaea that occurs on the sugar-phosphate backbone as opposed to the nucleobase moiety of DNA. While PT modification is widespread in the prokaryotic kingdom, how PT modifications are distributed in the genomes and their exact roles in the cell remain to be defined. In this study, we developed a simple and convenient technique called EcoWI-seq based on a modification-dependent restriction endonuclease to identify genomic positions of PT modifications. EcoWI-seq shows similar performance than other PT modification detection techniques and additionally, is easily scalable while requiring little starting material. As a proof of principle, we applied EcoWI-seq to map the PT modifications at base resolution in the genomes of both the Salmonella enterica cerro 87 and E. coli expressing the dnd+ gene cluster. Specifically, we address whether the partial establishment of modified PT positions is a stochastic or deterministic process. EcoWI-seq reveals a systematic usage of the same subset of target sites in clones for which the PT modification has been independently established.

Mechanism of DNA cleavage by the endonuclease SauUSI: a major barrier to horizontal gene transfer and antibiotic resistance in Staphylococcus aureus.

Acquisition of foreign DNA by Staphylococcus aureus, including vancomycin resistance genes, is thwarted by the ATP-dependent endonuclease SauUSI. Deciphering the mechanism of action of SauUSI could unravel the reason how it singularly plays a major role in preventing horizontal gene transfer (HGT) in S. aureus. Here, we report a detailed biochemical and structural characterization of SauUSI, which reveals that in the presence of ATP, the enzyme can cleave DNA having a single or multiple target site/s. Remarkably, in the case of multiple target sites, the entire region of DNA flanked by two target sites is shred into smaller fragments by SauUSI. Crystal structure of SauUSI reveals a stable dimer held together by the nuclease domains, which are spatially arranged to hydrolyze the phosphodiester bonds of both strands of the duplex. Thus, the architecture of the dimeric SauUSI facilitates cleavage of either single-site or multi-site DNA. The structure also provides insights into the molecular basis of target recognition by SauUSI. We show that target recognition activates ATP hydrolysis by the helicase-like ATPase domain, which powers active directional movement (translocation) of SauUSI along the DNA. We propose that a pile-up of multiple translocating SauUSI molecules against a stationary SauUSI bound to a target site catalyzes random double-stranded breaks causing shredding of the DNA between two target sites. The extensive and irreparable damage of the foreign DNA by shredding makes SauUSI a potent barrier against HGT.

Crystal structures of the EVE-HNH endonuclease VcaM4I in the presence and absence of DNA.

Many modification-dependent restriction endonucleases (MDREs) are fusions of a PUA superfamily modification sensor domain and a nuclease catalytic domain. EVE domains belong to the PUA superfamily, and are present in MDREs in combination with HNH nuclease domains. Here, we present a biochemical characterization of the EVE-HNH endonuclease VcaM4I and crystal structures of the protein alone, with EVE domain bound to either 5mC modified dsDNA or to 5mC/5hmC containing ssDNA. The EVE domain is moderately specific for 5mC/5hmC containing DNA according to EMSA experiments. It flips the modified nucleotide, to accommodate it in a hydrophobic pocket of the enzyme, primarily formed by P24, W82 and Y130 residues. In the crystallized conformation, the EVE domain and linker helix between the two domains block DNA binding to the catalytic domain. Removal of the EVE domain and inter-domain linker, but not of the EVE domain alone converts VcaM4I into a non-specific toxic nuclease. The role of the key residues in the EVE and HNH domains of VcaM4I is confirmed by digestion and restriction assays with the enzyme variants that differ from the wild-type by changes to the base binding pocket or to the catalytic residues.

A protein architecture guided screen for modification dependent restriction endonucleases.

Modification dependent restriction endonucleases (MDREs) often have separate catalytic and modification dependent domains. We systematically looked for previously uncharacterized fusion proteins featuring a PUA or DUF3427 domain and HNH or PD-(D/E)XK catalytic domain. The enzymes were clustered by similarity of their putative modification sensing domains into several groups. The TspA15I (VcaM4I, CmeDI), ScoA3IV (MsiJI, VcaCI) and YenY4I groups, all featuring a PUA superfamily domain, preferentially cleaved DNA containing 5-methylcytosine or 5-hydroxymethylcytosine. ScoA3V, also featuring a PUA superfamily domain, but of a different clade, exhibited 6-methyladenine stimulated nicking activity. With few exceptions, ORFs for PUA-superfamily domain containing endonucleases were not close to DNA methyltransferase ORFs, strongly supporting modification dependent activity of the endonucleases. DUF3427 domain containing fusion proteins had very little or no endonuclease activity, despite the presence of a putative PD-(D/E)XK catalytic domain. However, their expression potently restricted phage T4gt in Escherichia coli cells. In contrast to the ORFs for PUA domain containing endonucleases, the ORFs for DUF3427 fusion proteins were frequently found in defense islands, often also featuring DNA methyltransferases.

Crystal structure of the modification-dependent SRA-HNH endonuclease TagI.

TagI belongs to the recently characterized SRA-HNH family of modification-dependent restriction endonucleases (REases) that also includes ScoA3IV (Sco5333) and TbiR51I (Tbis1). Here, we present a crystal structure of dimeric TagI, which exhibits a DNA binding site formed jointly by the nuclease domains, and separate binding sites for modified DNA bases in the two protomers. The nuclease domains have characteristic features of HNH/ßßα-Me REases, and catalyze nicks or double strand breaks, with preference for /RY and RYN/RY sites, respectively. The SRA domains have the canonical fold. Their pockets for the flipped bases are spacious enough to accommodate 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC), but not glucosyl-5-hydroxymethylcytosine (g5hmC). Such preference is in agreement with the biochemical determination of the TagI modification dependence and the results of phage restriction assays. The ability of TagI to digest plasmids methylated by Dcm (C5mCWGG), M.Fnu4HI (G5mCNGC) or M.HpyCH4IV (A5mCGT) suggests that the SRA domains of the enzyme are tolerant to different sequence contexts of the modified base.

Activity and structure of EcoKMcrA.

Escherichia coli McrA (EcoKMcrA) acts as a methylcytosine and hydroxymethylcytosine dependent restriction endonuclease. We present a biochemical characterization of EcoKMcrA that includes the first demonstration of its endonuclease activity, small angle X-ray scattering (SAXS) data, and a crystal structure of the enzyme in the absence of DNA. Our data indicate that EcoKMcrA dimerizes via the anticipated C-terminal HNH domains, which together form a single DNA binding site. The N-terminal domains are not homologous to SRA domains, do not interact with each other, and have separate DNA binding sites. Electrophoretic mobility shift assay (EMSA) and footprinting experiments suggest that the N-terminal domains can sense the presence and sequence context of modified cytosines. Pyrrolocytosine fluorescence data indicate no base flipping. In vitro, EcoKMcrA DNA endonuclease activity requires Mn2+ ions, is not strictly methyl dependent, and is not observed when active site variants of the enzyme are used. In cells, EcoKMcrA specifically restricts DNA that is modified in the correct sequence context. This activity is impaired by mutations of the nuclease active site, unless the enzyme is highly overexpressed.

Restriction endonuclease AgeI is a monomer which dimerizes to cleave DNA.

Although all Type II restriction endonucleases catalyze phosphodiester bond hydrolysis within or close to their DNA target sites, they form different oligomeric assemblies ranging from monomers, dimers, tetramers to higher order oligomers to generate a double strand break in DNA. Type IIP restriction endonuclease AgeI recognizes a palindromic sequence 5΄-A/CCGGT-3΄ and cuts it ('/' denotes the cleavage site) producing staggered DNA ends. Here, we present crystal structures of AgeI in apo and DNA-bound forms. The structure of AgeI is similar to the restriction enzymes that share in their target sites a conserved CCGG tetranucleotide and a cleavage pattern. Structure analysis and biochemical data indicate, that AgeI is a monomer in the apo-form both in the crystal and in solution, however, it binds and cleaves the palindromic target site as a dimer. DNA cleavage mechanism of AgeI is novel among Type IIP restriction endonucleases.

Engineering nicking enzymes that preferentially nick 5-methylcytosine-modified DNA.

N.ϕGamma is a strand-specific and site-specific DNA nicking enzyme (YCG↓GT or AC↑CGR). Here we describe the isolation of single and double mutants of N.ϕGamma with attenuated activity. The nicking domains (NDs) of E59A and 11 double mutants were fused to the 5mCG-binding domain of MBD2 and generated fusion enzymes that preferentially nick 5mCG-modified DNA. The CG dinucleotide can be modified by C5 methyltransferases (MTases) such as M.SssI, M.HhaI or M.HpaII to create composite sites AC↑YGG N(8-15) 5mCG. We also constructed a fusion enzyme 2xMBD2-ND(N.BceSVIII) targeting more frequent composite sites AS↑YS N(5-12) 5mCG in Mn2+ buffer. 5mCG-dependent nicking requires special digestion conditions in high salt (0.3 M KCl) or in Ni2+ buffer. The fusion enzyme can be used to nick and label 5mCG-modified plasmid and genomic DNAs with fluorescently labeled Cy3-dUTP and potentially be useful for diagnostic applications, DNA sequencing and optical mapping of epigenetic markers. The importance of the predicted catalytic residues D89, H90, N106 and H115 in N.ϕGamma was confirmed by mutagenesis. We found that the wild-type enzyme N.ϕGamma prefers to nick 5mCG-modified DNA in Ni2+ buffer even though the nicking activity is sub-optimal compared to the activity in Mg2+ buffer.

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](2). 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.

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(2+) ion higher than 500 µM mediates promiscuous activity, Ca(2+) 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(2+)-mediated cleavage but imparting high cleavage fidelity with Mg(2+). 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(2+)-mediated cleavage activity by altering the geometry of the Ca(2+)-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.

The MspJI family of modification-dependent restriction endonucleases for epigenetic studies.

MspJI is a novel modification-dependent restriction endonuclease that cleaves at a fixed distance away from the modification site. Here, we present the biochemical characterization of several MspJI homologs, including FspEI, LpnPI, AspBHI, RlaI, and SgrTI. All of the enzymes specifically recognize cytosine C5 modification (methylation or hydroxymethylation) in DNA and cleave at a constant distance (N(12)/N(16)) away from the modified cytosine. Each displays its own sequence context preference, favoring different nucleotides flanking the modified cytosine. By cleaving on both sides of fully modified CpG sites, they allow the extraction of 32-base long fragments around the modified sites from the genomic DNA. These enzymes provide powerful tools for direct interrogation of the epigenome. For example, we show that RlaI, an enzyme that prefers (m)CWG but not (m)CpG sites, generates digestion patterns that differ between plant and mammalian genomic DNA, highlighting the difference between their epigenomic patterns. In addition, we demonstrate that deep sequencing of the digested DNA fragments generated from these enzymes provides a feasible method to map the modified sites in the genome. Altogether, the MspJI family of enzymes represent appealing tools of choice for method development in DNA epigenetic studies.

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.

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.

A type IV modification-dependent restriction enzyme SauUSI from Staphylococcus aureus subsp. aureus USA300.

A gene encoding a putative DNA helicase from Staphylococcus aureus USA300 was cloned and expressed in Escherichia coli. The protein was purified to over 90% purity by chromatography. The purified enzyme, SauUSI, predominantly cleaves modified DNA containing 5mC and 5-hydroxymethylcytosine. Cleavage of 5mC-modified plasmids indicated that the sites S5mCNGS (S = C or G) are preferentially digested. The endonuclease activity requires the presence of adenosine triphosphate (ATP) or dATP whereas the non-hydrolyzable γ-S-ATP does not support activity. SauUSI activity was inhibited by ethylenediaminetetraacetic acid. It is most active in Mg(++) buffers. No companion methylase gene was found near the SauUSI restriction gene. The absence of a cognate methylase and cleavage of modified DNA indicate that SauUSI belongs to type IV restriction endonucleases, a group that includes EcoK McrBC and Mrr. SauUSI belongs to a family of highly similar homologs found in other sequenced S. aureus, S. epidermidis and S. carnosus genomes. More distant SauUSI orthologs can be found in over 150 sequenced bacterial/archaea genomes. Finally, we demonstrated the biological function of the type IV REase in restricting 5mC-modified plasmid DNA by transformation into clinical S. aureus strain SA564, and in restricting phage λ infection when the endonuclease is expressed in E. coli.

Asymmetric DNA recognition by the OkrAI endonuclease, an isoschizomer of BamHI.

Restriction enzymes share little or no sequence homology with the exception of isoschizomers, or enzymes that recognize and cleave the same DNA sequence. We present here the structure of a BamHI isoschizomer, OkrAI, bound to the same DNA sequence (TATGGATCCATA) as that cocrystallized with BamHI. We show that OkrAI is a more minimal version of BamHI, lacking not only the N- and C-terminal helices but also an internal 3(10) helix and containing ß-strands that are shorter than those in BamHI. Despite these structural differences, OkrAI recognizes the DNA in a remarkably similar manner to BamHI, including asymmetric contacts via C-terminal 'arms' that appear to 'compete' for the minor groove. However, the arms are shorter than in BamHI. We observe similar DNA-binding affinities between OkrAI and BamHI but OkrAI has higher star activity (at 37°C) compared to BamHI. Together, the OkrAI and BamHI structures offer a rare opportunity to compare two restriction enzymes that work on exactly the same DNA substrate.

Characterization and crystal structure of the type IIG restriction endonuclease RM.BpuSI.

A type IIG restriction endonuclease, RM.BpuSI from Bacillus pumilus, has been characterized and its X-ray crystal structure determined at 2.35Å resolution. The enzyme is comprised of an array of 5-folded domains that couple the enzyme's N-terminal endonuclease domain to its C-terminal target recognition and methylation activities. The REase domain contains a PD-x(15)-ExK motif, is closely superimposable against the FokI endonuclease domain, and coordinates a single metal ion. A helical bundle domain connects the endonuclease and methyltransferase (MTase) domains. The MTase domain is similar to the N6-adenine MTase M.TaqI, while the target recognition domain (TRD or specificity domain) resembles a truncated S subunit of Type I R-M system. A final structural domain, that may form additional DNA contacts, interrupts the TRD. DNA binding and cleavage must involve large movements of the endonuclease and TRD domains, that are probably tightly coordinated and coupled to target site methylation status.

Universal NicE-Seq: A Simple and Quick Method for Accessible Chromatin Detection in Fixed Cells.

Genome-wide accessible chromatin sequencing and identification has enabled deciphering the epigenetic information encoded in chromatin, revealing accessible promoters, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding. The starting biological materials are often fixed using formaldehyde crosslinking. Here, we describe accessible chromatin library preparation from low numbers of formaldehyde-crosslinked cells using a modified nick translation method, where a nicking enzyme nicks one strand of DNA and DNA polymerase incorporates biotin-conjugated dATP, dCTP, and methyl-dCTP. Once the DNA is labeled, it can be isolated for NGS library preparation. We termed this method as universal NicE-seq (nicking enzyme-assisted sequencing). We also demonstrate a single tube method that enables direct NGS library preparation from low cell numbers without DNA purification. Furthermore, we demonstrated universal NicE-seq on FFPE tissue section sample.

Characterization of type II and III restriction-modification systems from Bacillus cereus strains ATCC 10987 and ATCC 14579.

The genomes of two Bacillus cereus strains (ATCC 10987 and ATCC 14579) have been sequenced. Here, we report the specificities of type II/III restriction (R) and modification (M) enzymes. Found in the ATCC 10987 strain, BceSI is a restriction endonuclease (REase) with the recognition and cut site CGAAG 24-25/27-28. BceSII is an isoschizomer of AvaII (G/GWCC). BceSIII cleaves at ACGGC 12/14. The BceSIII C terminus resembles the catalytic domains of AlwI, MlyI, and Nt.BstNBI. BceSIV is composed of two subunits and cleaves on both sides of GCWGC. BceSIV activity is strongly stimulated by the addition of cofactor ATP or GTP. The large subunit (R1) of BceSIV contains conserved motifs of NTPases and DNA helicases. The R1 subunit has no endonuclease activity by itself; it strongly stimulates REase activity when in complex with the R2 subunit. BceSIV was demonstrated to hydrolyze GTP and ATP in vitro. BceSIV is similar to CglI (GCSGC), and homologs of R1 are found in 11 sequenced bacterial genomes, where they are paired with specificity subunits. In addition, homologs of the BceSIV R1-R2 fusion are found in many sequenced microbial genomes. An orphan methylase, M.BceSV, was found to modify GCNGC, GGCC, CCGG, GGNNCC, and GCGC sites. A ParB-methylase fusion protein appears to nick DNA nonspecifically. The ATCC 14579 genome encodes an active enzyme Bce14579I (GCWGC). BceSIV and Bce14579I belong to the phospholipase D (PLD) family of endonucleases that are widely distributed among Bacteria and Archaea. A survey of type II and III restriction-modification (R-M) system genes is presented from sequenced B. cereus, Bacillus anthracis, and Bacillus thuringiensis strains.

Characterization of MspNI (G/GWCC) and MspNII (R/GATCY), novel thermostable Type II restriction endonucleases from Meiothermus sp., isoschizomers of AvaII and BstYI.

MspNI and MspNII, isoschizomers of prototype Type II restriction endonucleases AvaII and BstYI, were extracted from an extreme thermophile bacterium belonging to the genus Meiothermus, isolated from the hot sulphur springs in north Himalayan region of India where temperature and pH ranged from 60 to 80°C and 7.5 to 8.5, respectively. The two enzymes were purified to homogeneity using Cibacron-Blue 3GA Agarose, Q-Sepharose and SP-Sepharose chromatography and were homodimers with subunit molecular weights of 27 and 45 kDa, respectively. Restriction mapping and run-off sequencing of MspNI and MspNII cleaved pBR322 DNA showed that they recognized and cleaved 5'-G/GWCC-3' and 5'-R/GATCY-3' sites, respectively. MspNI and MspNII worked optimally at 60 and 70°C, 6 and 5 mM MgCl(2), respectively and showed no star activity in organic solvents. Both were resistant to sequence methylation and were stable up to 25 PCR cycles.

Engineering Nt.BtsCI and Nb.BtsCI nicking enzymes and applications in generating long overhangs.

Type IIS restriction endonuclease BtsCI (GGATG 2/0) is a neoschizomer of FokI (GGATG 9/13) and cleaves closer to the recognition sequence. Although M.BtsCI shows 62% amino acid sequence identity to M.FokI, BtsCI and FokI restriction endonucleases do not share significant amino acid sequence similarity. BtsCI belongs to a group of Type IIS restriction endonucleases, BsmI, Mva1269I and BsrI, that carry two different catalytic sites in a single polypeptide. By inactivating one of the catalytic sites through mutagenesis, we have generated nicking variants of BtsCI that specifically nick the bottom-strand or the top-strand of the target site. By treating target DNA sequentially with the appropriate combinations of FokI and BtsCI nicking variants, we are able to generate long overhangs suitable for fluorescent labeling through end-filling or other techniques based on annealing of complementary DNA sequences.