Structure and Mechanism of a Cyclic Trinucleotide-Activated Bacterial Endonuclease Mediating Bacteriophage Immunity.

Bacteria possess an array of defenses against foreign invaders, including a broadly distributed bacteriophage defense system termed CBASS (cyclic oligonucleotide-based anti-phage signaling system). In CBASS systems, a cGAS/DncV-like nucleotidyltransferase synthesizes cyclic di- or tri-nucleotide second messengers in response to infection, and these molecules activate diverse effectors to mediate bacteriophage immunity via abortive infection. Here, we show that the CBASS effector NucC is related to restriction enzymes but uniquely assembles into a homotrimer. Binding of NucC trimers to a cyclic tri-adenylate second messenger promotes assembly of a NucC homohexamer competent for non-specific double-strand DNA cleavage. In infected cells, NucC activation leads to complete destruction of the bacterial chromosome, causing cell death prior to completion of phage replication. In addition to CBASS systems, we identify NucC homologs in over 30 type III CRISPR/Cas systems, where they likely function as accessory nucleases activated by cyclic oligoadenylate second messengers synthesized by these systems' effector complexes.

Structural analysis of the BisI family of modification dependent restriction endonucleases.

The BisI family of restriction endonucleases is unique in requiring multiple methylated or hydroxymethylated cytosine residues within a short recognition sequence (GCNGC), and in cleaving directly within this sequence, rather than at a distance. Here, we report that the number of modified cytosines that are required for cleavage can be tuned by the salt concentration. We present crystal structures of two members of the BisI family, NhoI and Eco15I_Ntd (N-terminal domain of Eco15I), in the absence of DNA and in specific complexes with tetra-methylated GCNGC target DNA. The structures show that NhoI and Eco15I_Ntd sense modified cytosine bases in the context of double-stranded DNA (dsDNA) without base flipping. In the co-crystal structures of NhoI and Eco15I_Ntd with DNA, the internal methyl groups (G5mCNGC) interact with the side chains of an (H/R)(V/I/T/M) di-amino acid motif near the C-terminus of the distal enzyme subunit and arginine residue from the proximal subunit. The external methyl groups (GCNG5mC) interact with the proximal enzyme subunit, mostly through main chain contacts. Surface plasmon resonance analysis for Eco15I_Ntd shows that the internal and external methyl binding pockets contribute about equally to sensing of cytosine methyl groups.

A sensitive approach to map genome-wide 5-hydroxymethylcytosine and 5-formylcytosine at single-base resolution.

Mapping genome-wide 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC) at single-base resolution is important to understand their biological functions. We present a cost-efficient mapping method that combines 5hmC-specific restriction enzyme PvuRts1I with a 5hmC chemical labeling enrichment method. The sensitive method enables detection of low-abundance 5hmC sites, providing a more complete 5hmC landscape than available bisulfite-based methods. This method generated a genome-wide 5fC map at single-base resolution. Parallel analyses revealed that 5hmC and 5fC in non-CpG context exhibit lower abundance, more dynamically, than those in CpG context. In the genic region, distribution of 5hmCpG and 5fCpG differed from 5hmCH and 5fCH (H = A, T, C). 5hmC and 5fC were distributed distinctly at regulatory protein-DNA binding sites, depleted in permissive transcription factor binding sites, and enriched at active and poised enhancers. This sensitive bisulfite conversion-free method can be applied to biological samples with limited starting material or low-abundance cytosine modifications.

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.

DNA backbone interactions impact the sequence specificity of DNA sulfur-binding domains: revelations from structural analyses.

The sulfur atom of phosphorothioated DNA (PT-DNA) is coordinated by a surface cavity in the conserved sulfur-binding domain (SBD) of type IV restriction enzymes. However, some SBDs cannot recognize the sulfur atom in some sequence contexts. To illustrate the structural determinants for sequence specificity, we resolved the structure of SBDSpr, from endonuclease SprMcrA, in complex with DNA of GPSGCC, GPSATC and GPSAAC contexts. Structural and computational analyses explained why it binds the above PT-DNAs with an affinity in a decreasing order. The structural analysis of SBDSpr-GPSGCC and SBDSco-GPSGCC, the latter only recognizes DNA of GPSGCC, revealed that a positively charged loop above the sulfur-coordination cavity electrostatically interacts with the neighboring DNA phosphate linkage. The structural analysis indicated that the DNA-protein hydrogen bonding pattern and weak non-bonded interaction played important roles in sequence specificity of SBD protein. Exchanges of the positively-charged amino acid residues with the negatively-charged residues in the loop would enable SBDSco to extend recognization for more PT-DNA sequences, implying that type IV endonucleases can be engineered to recognize PT-DNA in novel target sequences.

DNA-mediated coupling of ATPase, translocase and nuclease activities of a Type ISP restriction-modification enzyme.

Enzymes involved in nucleic acid transactions often have a helicase-like ATPase coordinating and driving their functional activities, but our understanding of the mechanistic details of their coordination is limited. For example, DNA cleavage by the antiphage defense system Type ISP restriction-modification enzyme requires convergence of two such enzymes that are actively translocating on DNA powered by Superfamily 2 ATPases. The ATPase is activated when the enzyme recognizes a DNA target sequence. Here, we show that the activation is a two-stage process of partial ATPase stimulation upon recognition of the target sequence by the methyltransferase and the target recognition domains, and complete stimulation that additionally requires the DNA to interact with the ATPase domain. Mutagenesis revealed that a ß-hairpin loop and motif V of the ATPase couples DNA translocation to ATP hydrolysis. Deletion of the loop inhibited translocation, while mutation of motif V slowed the rate of translocation. Both the mutations inhibited the double-strand (ds) DNA cleavage activity of the enzyme. However, a translocating motif V mutant cleaved dsDNA on encountering a translocating wild-type enzyme. Based on these results, we conclude that the ATPase-driven translocation not only brings two nucleases spatially close to catalyze dsDNA break, but that the rate of translocation influences dsDNA cleavage.

Structure of HhaI endonuclease with cognate DNA at an atomic resolution of 1.0 Å.

HhaI, a Type II restriction endonuclease, recognizes the symmetric sequence 5'-GCG↓C-3' in duplex DNA and cleaves ('↓') to produce fragments with 2-base, 3'-overhangs. We determined the structure of HhaI in complex with cognate DNA at an ultra-high atomic resolution of 1.0 Å. Most restriction enzymes act as dimers with two catalytic sites, and cleave the two strands of duplex DNA simultaneously, in a single binding event. HhaI, in contrast, acts as a monomer with only one catalytic site, and cleaves the DNA strands sequentially, one after the other. HhaI comprises three domains, each consisting of a mixed five-stranded ß sheet with a defined function. The first domain contains the catalytic-site; the second contains residues for sequence recognition; and the third contributes to non-specific DNA binding. The active-site belongs to the 'PD-D/EXK' superfamily of nucleases and contains the motif SD-X11-EAK. The first two domains are similar in structure to two other monomeric restriction enzymes, HinP1I (G↓CGC) and MspI (C↓CGG), which produce fragments with 5'-overhangs. The third domain, present only in HhaI, shifts the positions of the recognition residues relative to the catalytic site enabling this enzyme to cleave the recognition sequence at a different position. The structure of M.HhaI, the biological methyltransferase partner of HhaI, was determined earlier. Together, these two structures represent the first natural pair of restriction-modification enzymes to be characterized in atomic detail.

The structure of the Thermococcus gammatolerans McrB N-terminal domain reveals a new mode of substrate recognition and specificity among McrB homologs.

McrBC is a two-component, modification-dependent restriction system that cleaves foreign DNA-containing methylated cytosines. Previous crystallographic studies have shown that Escherichia coli McrB uses a base-flipping mechanism to recognize these modified substrates with high affinity. The side chains stabilizing both the flipped base and the distorted duplex are poorly conserved among McrB homologs, suggesting that other mechanisms may exist for binding modified DNA. Here we present the structures of the Thermococcus gammatolerans McrB DNA-binding domain (TgΔ185) both alone and in complex with a methylated DNA substrate at 1.68 and 2.27 Å resolution, respectively. The structures reveal that TgΔ185 consists of a YT521-B homology (YTH) domain, which is commonly found in eukaryotic proteins that bind methylated RNA and is structurally unrelated to the E. coli McrB DNA-binding domain. Structural superposition and co-crystallization further show that TgΔ185 shares a conserved aromatic cage with other YTH domains, which forms the binding pocket for a flipped-out base. Mutational analysis of this aromatic cage supports its role in conferring specificity for the methylated adenines, whereas an extended basic surface present in TgΔ185 facilitates its preferential binding to duplex DNA rather than RNA. Together, these findings establish a new binding mode and specificity among McrB homologs and expand the biological roles of YTH domains.

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.

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).

Crystal structure of the EcoKMcrA N-terminal domain (NEco): recognition of modified cytosine bases without flipping.

EcoKMcrA from Escherichia coli restricts CpG methylated or hydroxymethylated DNA, and may act as a barrier against host DNA. The enzyme consists of a novel N-terminal specificity domain that we term NEco, and a C-terminal catalytic HNH domain. Here, we report that NEco and full-length EcoKMcrA specificities are consistent. NEco affinity to DNA increases more from hemi- to full-methylation than from non- to hemi-methylation, indicating cooperative binding of the methyl groups. We determined the crystal structures of NEco in complex with fully modified DNA containing three variants of the Y5mCGR EcoKMcrA target sequence: C5mCGG, T5mCGA and T5hmCGA. The structures explain the specificity for the two central base pairs and one of the flanking pairs. As predicted based on earlier biochemical experiments, NEco does not flip any DNA bases. The proximal and distal methyl groups are accommodated in separate pockets. Changes to either pocket reduce DNA binding by NEco and restriction by EcoKMcrA, confirming the relevance of the crystallographically observed binding mode in solution.

Systematic microsatellite repeat expansion cloning and validation.

Approximately 3% of the human genome is composed of short tandem repeat (STR) DNA sequence known as microsatellites, which can be found in both coding and non-coding regions. When associated with genic regions, expansion of microsatellite repeats beyond a critical threshold causes dozens of neurological repeat expansion disorders. To better understand the molecular pathology of repeat expansion disorders, precise cloning of microsatellite repeat sequence and expansion size is highly valuable. Unfortunately, cloning repeat expansions is often challenging and presents a significant bottleneck to practical investigation. Here, we describe a clear method for seamless and systematic cloning of practically any microsatellite repeat expansion. We use cloning and expansion of GGGGCC repeats, which are the leading genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), as an example. We employ a recursive directional ligation (RDL) technique to build multiple GGGGCC repeat-containing vectors. We describe methods to validate repeat expansion cloning, including diagnostic restriction digestion, PCR across the repeat, and next-generation long-read MinION nanopore sequencing. Validated cloning of microsatellite repeats beyond the critical expansion threshold can facilitate step-by-step characterization of disease mechanisms at the cellular and molecular level.

RNA-DNA hybrid (R-loop) immunoprecipitation mapping: an analytical workflow to evaluate inherent biases.

The impact of R-loops on the physiology and pathology of chromosomes has been demonstrated extensively by chromatin biology research. The progress in this field has been driven by technological advancement of R-loop mapping methods that largely relied on a single approach, DNA-RNA immunoprecipitation (DRIP). Most of the DRIP protocols use the experimental design that was developed by a few laboratories, without paying attention to the potential caveats that might affect the outcome of RNA-DNA hybrid mapping. To assess the accuracy and utility of this technology, we pursued an analytical approach to estimate inherent biases and errors in the DRIP protocol. By performing DRIP-sequencing, qPCR, and receiver operator characteristic (ROC) analysis, we tested the effect of formaldehyde fixation, cell lysis temperature, mode of genome fragmentation, and removal of free RNA on the efficacy of RNA-DNA hybrid detection and implemented workflows that were able to distinguish complex and weak DRIP signals in a noisy background with high confidence. We also show that some of the workflows perform poorly and generate random answers. Furthermore, we found that the most commonly used genome fragmentation method (restriction enzyme digestion) led to the overrepresentation of lengthy DRIP fragments over coding ORFs, and this bias was enhanced at the first exons. Biased genome sampling severely compromised mapping resolution and prevented the assignment of precise biological function to a significant fraction of R-loops. The revised workflow presented herein is established and optimized using objective ROC analyses and provides reproducible and highly specific RNA-DNA hybrid detection.

Binary control of enzymatic cleavage of DNA origami by structural antideterminants.

Controlling DNA nanostructure interaction with protein is essential in developing nanodevices with programmable function, reactivity, and stability for biological and medical applications. Here, we show that the sequence-specific action of restriction endonucleases towards sharp triangular or rectangular DNA origami exhibits a novel, binary 'on/off' behaviour, as canonical recognition sites are either essentially fully reactive, or strongly resistant to enzymatic cutting. Moreover, introduction of structural defects in the sharp triangle can activate an otherwise unreactive site, with a site-to-defect distance of ∼50 nm. We argue that site reactivity is dependent upon programmable, mechanical coupling in the two-dimensional DNA origami, with specific structural elements, including DNA nicks and branches proximal to the sites that can function as negative(anti) determinants of reactivity. Empirically modelling the constraints to DNA degrees of freedom associated with each recognition site in the sharp triangle can rationalize the pattern of suppressed reactivity towards nine restriction endonucleases, in substantial agreement with the experimental results. These results provide a basis for a predictive understanding of structure-reactivity correlates of specific DNA nanostructures, which will allow a better understanding of the behaviour of nucleic acids under nanoscale confinement, as well as in the rational design of functional nanodevices based on self-assembling nucleic acids.

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.

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.

The H-subunit of the restriction endonuclease CglI contains a prototype DEAD-Z1 helicase-like motor.

CglI is a restriction endonuclease from Corynebacterium glutamicum that forms a complex between: two R-subunits that have site specific-recognition and nuclease domains; and two H-subunits, with Superfamily 2 helicase-like DEAD domains, and uncharacterized Z1 and C-terminal domains. ATP hydrolysis by the H-subunits catalyses dsDNA translocation that is necessary for long-range movement along DNA that activates nuclease activity. Here, we provide biochemical and molecular modelling evidence that shows that Z1 has a fold distantly-related to RecA, and that the DEAD-Z1 domains together form an ATP binding interface and are the prototype of a previously undescribed monomeric helicase-like motor. The DEAD-Z1 motor has unusual Walker A and Motif VI sequences those nonetheless have their expected functions. Additionally, it contains DEAD-Z1-specific features: an H/H motif and a loop (aa 163-aa 172), that both play a role in the coupling of ATP hydrolysis to DNA cleavage. We also solved the crystal structure of the C-terminal domain which has a unique fold, and demonstrate that the Z1-C domains are the principal DNA binding interface of the H-subunit. Finally, we use small angle X-ray scattering to provide a model for how the H-subunit domains are arranged in a dimeric complex.

Engineering altered protein-DNA recognition specificity.

Protein engineering is used to generate novel protein folds and assemblages, to impart new properties and functions onto existing proteins, and to enhance our understanding of principles that govern protein structure. While such approaches can be employed to reprogram protein-protein interactions, modifying protein-DNA interactions is more difficult. This may be related to the structural features of protein-DNA interfaces, which display more charged groups, directional hydrogen bonds, ordered solvent molecules and counterions than comparable protein interfaces. Nevertheless, progress has been made in the redesign of protein-DNA specificity, much of it driven by the development of engineered enzymes for genome modification. Here, we summarize the creation of novel DNA specificities for zinc finger proteins, meganucleases, TAL effectors, recombinases and restriction endonucleases. The ease of re-engineering each system is related both to the modularity of the protein and the extent to which the proteins have evolved to be capable of readily modifying their recognition specificities in response to natural selection. The development of engineered DNA binding proteins that display an ideal combination of activity, specificity, deliverability, and outcomes is not a fully solved problem, however each of the current platforms offers unique advantages, offset by behaviors and properties requiring further study and development.

Structure and operation of the DNA-translocating type I DNA restriction enzymes.

Type I DNA restriction/modification (RM) enzymes are molecular machines found in the majority of bacterial species. Their early discovery paved the way for the development of genetic engineering. They control (restrict) the influx of foreign DNA via horizontal gene transfer into the bacterium while maintaining sequence-specific methylation (modification) of host DNA. The endonuclease reaction of these enzymes on unmethylated DNA is preceded by bidirectional translocation of thousands of base pairs of DNA toward the enzyme. We present the structures of two type I RM enzymes, EcoKI and EcoR124I, derived using electron microscopy (EM), small-angle scattering (neutron and X-ray), and detailed molecular modeling. DNA binding triggers a large contraction of the open form of the enzyme to a compact form. The path followed by DNA through the complexes is revealed by using a DNA mimic anti-restriction protein. The structures reveal an evolutionary link between type I RM enzymes and type II RM enzymes.

Iteratively improving Hi-C experiments one step at a time.

The 3D organization of eukaryotic chromosomes affects key processes such as gene expression, DNA replication, cell division, and response to DNA damage. The genome-wide chromosome conformation capture (Hi-C) approach can characterize the landscape of 3D genome organization by measuring interaction frequencies between all genomic regions. Hi-C protocol improvements and rapid advances in DNA sequencing power have made Hi-C useful to study diverse biological systems, not only to elucidate the role of 3D genome structure in proper cellular function, but also to characterize genomic rearrangements, assemble new genomes, and consider chromatin interactions as potential biomarkers for diseases. Yet, the Hi-C protocol is still complex and subject to variations at numerous steps that can affect the resulting data. Thus, there is still a need for better understanding and control of factors that contribute to Hi-C experiment success and data quality. Here, we evaluate recently proposed Hi-C protocol modifications as well as often overlooked variables in sample preparation and examine their effects on Hi-C data quality. We examine artifacts that can occur during Hi-C library preparation, including microhomology-based artificial template copying and chimera formation that can add noise to the downstream data. Exploring the mechanisms underlying Hi-C artifacts pinpoints steps that should be further optimized in the future. To improve the utility of Hi-C in characterizing the 3D genome of specialized populations of cells or small samples of primary tissue, we identify steps prone to DNA loss which should be considered to adapt Hi-C to lower cell numbers.