The PD-(D/E)XK superfamily revisited: identification of new members among proteins involved in DNA metabolism and functional predictions for domains of (hitherto) unknown function

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.

High-resolution temporal profiling of E. coli transcriptional response

Understanding how cells dynamically adapt to their environment is a primary focus of biology research. Temporal information about cellular behavior is often limited by both small numbers of data time-points and the methods used to analyze this data. Here, we apply unsupervised machine learning to a data set containing the activity of 1805 native promoters in E. coli measured every 10 minutes in a high-throughput microfluidic device via fluorescence time-lapse microscopy. Specifically, this data set reveals E. coli transcriptome dynamics when exposed to different heavy metal ions. We use a bioinformatics pipeline based on Independent Component Analysis (ICA) to generate insights and hypotheses from this data. We discovered three primary, time-dependent stages of promoter activation to heavy metal stress (fast, intermediate, and steady). Furthermore, we uncovered a global strategy E. coli uses to reallocate resources from stress-related promoters to growth-related promoters following exposure to heavy metal stress.

Molecular Biology of Cytoplasmic Incompatibility Caused by Wolbachia Endosymbionts

Among endosymbiotic bacteria living within eukaryotic cells, Wolbachia is exceptionally widespread, particularly in arthropods. Inherited through the female germline, it has evolved ways to increase the fraction of bacterially infected offspring by inducing parthenogenesis, feminization, male killing, or, most commonly, cytoplasmic incompatibility (CI). In CI, Wolbachia infection of males causes embryonic lethality unless they mate with similarly infected females, creating a relative reproductive advantage for infected females. A set of related Wolbachia bicistronic operons encodes the CI-inducing factors. The downstream gene encodes a deubiquitylase or nuclease and is responsible for CI induction by males, while the upstream product when expressed in females binds its sperm-introduced cognate partner and rescues viability. Both toxin-antidote and host-modification mechanisms have been proposed to explain CI. Interestingly, male killing by either Spiroplasma or Wolbachia endosymbionts involves deubiquitylases as well. Interference with the host ubiquitin system may therefore be a common theme among endosymbiont-mediated reproductive alterations.

The novel anti-phage system Shield co-opts an RmuC domain to mediate phage defense across Pseudomonas species

Competitive bacteria-bacteriophage interactions have resulted in the evolution of a plethora of bacterial defense systems preventing phage propagation. In recent years, computational and bioinformatic approaches have underpinned the discovery of numerous novel bacterial defense systems. Anti-phage systems are frequently encoded together in genomic loci termed defense islands. Here we report the identification and characterisation of a novel anti-phage system, that we have termed Shield, which forms part of the Pseudomonas defensive arsenal. The Shield system comprises the core component ShdA, a membrane-bound protein harboring an RmuC domain. Heterologous production of ShdA alone is sufficient to mediate bacterial immunity against several phages. We demonstrate that Shield and ShdA confer population-level immunity and that they can also decrease transformation efficiency. We further show that ShdA homologues can degrade DNA in vitro and, when expressed in a heterologous host, can alter the organisation of the host chromosomal DNA. Use of comparative genomic approaches identified how Shield can be divided into four subtypes, three of which contain additional components that in some cases can negatively affect the activity of ShdA and/or provide additional lines of phage defense. Collectively, our results identify a new player within the Pseudomonas bacterial immunity arsenal that displays a novel mechanism of protection, and reveals a role for RmuC domains in phage defense.

A broadly distributed predicted helicase/nuclease confers phage resistance via abortive infection

There is strong selection for the evolution of systems that protect bacterial populations from viral attack. We report a single phage defense protein, Hna, that provides protection against diverse phages in Sinorhizobium meliloti, a nitrogen-fixing alpha-proteobacterium. Homologs of Hna are distributed widely across bacterial lineages, and a homologous protein from Escherichia coli also confers phage defense. Hna contains superfamily II helicase motifs at its N terminus and a nuclease motif at its C terminus, with mutagenesis of these motifs inactivating viral defense. Hna variably impacts phage DNA replication but consistently triggers an abortive infection response in which infected cells carrying the system die but do not release phage progeny. A similar host cell response is triggered in cells containing Hna upon expression of a phage-encoded single-stranded DNA binding protein (SSB), independent of phage infection. Thus, we conclude that Hna limits phage spread by initiating abortive infection in response to a phage protein.

Comparative Genomic Analysis and Physiological Properties of Limosilactobacillus fermentum SMFM2017-NK2 with Ability to Inflammatory Bowel Disease

The objective of this study was to evaluate the anti-inflammatory effect of Latilactobacillus sakei SMFM2017-NK1 (LS1), L. sakei SMFM2017-NK3 (LS2), and Limosilactobacillus fermentum SMFM2017-NK2 (LF) on colitis using an animal model. DSS (dextran sulfate sodium salt) was orally injected into C57BL/6N mice to induce inflammation in the colon for seven days. Colitis mice models were treated with LS1, LS2, and LF, respectively, and Lacticaseibacillus rhamnosus GG (LGG) was used as a positive control. During oral administration of lactic acid bacteria, the weights of the mice were measured, and the disease activity index (DAI) score was determined by judging the degree of diarrhea and bloody stool. When comparing the differences between the minimum weight after DSS administration and the maximum weight after lactic acid bacteriaadministration were compared, the LF-treated group showed the highest weight gain at 8.91%. The DAI scores of the LF, LS2, and LGG groups were lower than that of the control group. After sacrifice, mRNA expression levels for proinflammatory cytokines (TNF-α, IL-1β, IL-6, and IFN-γ) and mediators (iNOS and COX-2) in the colon were measured. LF was selected as a superior strain for anti-inflammation in the colon. It was further analyzed to determine its biochemical characteristics, cytotoxicity, and thermal stability. Catalase and oxidase activities for LF were negative. In cytotoxicity and heat stability tests, the LF group had higher cell viability than the LGG group. The genome of LF was obtained, and 5682 CDS, 114 tRNA, 2 RNA, and 5 repeat regions were predicted. Especially, LF could be distinguished from the other three L. fermentum strains based on taxonomic profiling, specific orthologous genes of the strain, and genomic variants. The results of this study suggest that L. fermentum SMFM2017-NK2 is a novel strain with an anti-inflammatory effect on colitis.

Beyond a platform protein for the degradosome assembly: The Apoptosis-Inducing Factor as an efficient nuclease involved in chromatinolysis

The Apoptosis-Inducing Factor (AIF) is a moonlighting flavoenzyme involved in the assembly of mitochondrial respiratory complexes in healthy cells, but also able to trigger DNA cleavage and parthanatos. Upon apoptotic-stimuli, AIF redistributes from the mitochondria to the nucleus, where upon association with other proteins such as endonuclease CypA and histone H2AX, it is proposed to organize a DNA-degradosome complex. In this work, we provide evidence for the molecular assembly of this complex as well as for the cooperative effects among its protein components to degrade genomic DNA into large fragments. We have also uncovered that AIF has nuclease activity that is stimulated in the presence of either Mg²⁺ or Ca²⁺. Such activity allows AIF by itself and in cooperation with CypA to efficiently degrade genomic DNA. Finally, we have identified TopIB and DEK motifs in AIF as responsible for its nuclease activity. These new findings point, for the first time, to AIF as a nuclease able to digest nuclear dsDNA in dying cells, improving our understanding of its role in promoting apoptosis and opening paths for the development of new therapeutic strategies.

Genomic Insights into the Radiation-Resistant Capability of Sphingomonas qomolangmaensis S5-59T and Sphingomonas glaciei S8-45T, Two Novel Bacteria from the North Slope of Mount Everest

Mount Everest provides natural advantages to finding radiation-resistant extremophiles that are functionally mechanistic and possess commercial significance. (1) Background: Two bacterial strains, designated S5-59T and S8-45T, were isolated from moraine samples collected from the north slope of Mount Everest at altitudes of 5700m and 5100m above sea level. (2) Methods: The present study investigated the polyphasic features and genomic characteristics of S5-59^T and S8-45^T. (3) Results: The major fatty acids and the predominant respiratory menaquinone of S5-59^T and S8-45^T were summed as feature 3 (comprising C16:1 ω6c and/or C16:1 ω7c) and ubiquinone-10 (Q-10). Phylogenetic analyses based on 16S rRNA sequences and average nucleotide identity values among these two strains and their reference type strains were below the species demarcation thresholds of 98.65% and 95%. Strains S5-59^T and S8-45^T harbored great radiation resistance. The genomic analyses showed that DNA damage repair genes, such as mutL, mutS, radA, radC, recF, recN, etc., were present in the S5-59^T and S8-45^T strains. Additionally, strain S5-59^T possessed more genes related to DNA protection proteins. The pan-genome analysis and horizontal gene transfers revealed that strains of Sphingomonas had a consistently homologous genetic evolutionary radiation resistance. Moreover, enzymatic antioxidative proteins also served critical roles in converting ROS into harmless molecules that resulted in resistance to radiation. Further, pigments and carotenoids such as zeaxanthin and alkylresorcinols of the non-enzymatic antioxidative system were also predicted to protect them from radiation. (4) Conclusions: Type strains S5-59^T (=JCM 35564T =GDMCC 1.3193T) and S8-45^T (=JCM 34749T =GDMCC 1.2715T) represent two novel species of the genus Sphingomonas with the proposed name Sphingomonas qomolangmaensis sp. nov. and Sphingomonas glaciei sp. nov. The type strains, S5-59^T and S8-45^T, were assessed in a deeply genomic study of their radiation-resistant mechanisms and this thus resulted in a further understanding of their greater potential application for the development of anti-radiation protective drugs.

Muramidase, nuclease, or hypothetical protein genes intervene between paired genes encoding DNA packaging terminase and portal proteins in Wolbachia phages and prophages

Genomes of the obligate intracellular alpha proteobacterium Wolbachia pipientis often encode prophage-like regions, and in a few cases, purified particles have been recovered. Because the structure of a conserved WO phage genome has been difficult to establish, we examined paired terminase and portal genes in Wolbachia phages and prophages, relative to those encoded by the gene transfer agent RcGTA from the free-living alpha proteobacterium Rhodobacter capsulatus. Terminase and portal proteins from Wolbachia have higher similarity to orthologs encoded by RcGTA than to orthologs encoded by bacteriophage lambda. In lambdoid phages, these proteins play key roles in assembly of mature phage particles, while in less well-studied gene transfer agents, terminase and portal proteins package random fragments of bacterial DNA, which could confound elucidation of WO phage genomes. In WO phages and prophages, terminase genes followed by a short gpW gene may be separated from the downstream portal gene by open-reading frames encoding a GH_25 hydrolase/muramidase, a PD-(D/E)XK nuclease, a hypothetical protein and/or a RelE/ParE toxin-antitoxin module. These aspects of gene organization, coupled with evidence for a low, non-inducible yield of WO phages, and the small size of WO phage particles described in the literature raise the possibility that Wolbachia prophage regions participate in processes that extend beyond conventional bacteriophage lysogeny and lytic replication. These intervening genes, and their possible relation to functions associated with GTAs, may contribute to variability among WO phage genomes recovered from physical particles and impact the ability of WO phages to act as transducing agents.

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.

LC_Glucose-Inhibited Division Protein Is Required for Motility, Biofilm Formation, and Stress Response in Lysobacter capsici X2-3

Glucose-inhibited division protein (GidA) plays a critical role in the growth, stress response, and virulence of bacteria. However, how gidA may affect plant growth-promoting bacteria (PGPB) is still not clear. Our study aimed to describe the regulatory function of the gidA gene in Lysobacter capsici, which produces a variety of lytic enzymes and novel antibiotics. Here, we generated an LC_GidA mutant, MT16, and an LC_GidA complemented strain, Com-16, by plasmid integration. The deletion of LC_GidA resulted in an attenuation of the bacterial growth rate, motility, and biofilm formation of L. capsici. Root colonization assays demonstrated that the LC_GidA mutant showed reduced colonization of wheat roots. In addition, disruption of LC_GidA showed a clear diminution of survival in the presence of high temperature, high salt, and different pH conditions. The downregulated expression of genes related to DNA replication, cell division, motility, and biofilm formation was further validated by real-time quantitative PCR (RT–qPCR). Together, understanding the regulatory function of GidA is helpful for improving the biocontrol of crop diseases and has strong potential for biological applications.

Crystal structure and initial characterization of a novel archaeal-like Holliday junction-resolving enzyme from Thermus thermophilus phage Tth15-6

This study describes the production, characterization and structure determination of a novel Holliday junction-resolving enzyme. The enzyme, termed Hjc_15-6, is encoded in the genome of phage Tth15-6, which infects Thermus thermophilus. Hjc_15-6 was heterologously produced in Escherichia coli and high yields of soluble and biologically active recombinant enzyme were obtained in both complex and defined media. Amino-acid sequence and structure comparison suggested that the enzyme belongs to a group of enzymes classified as archaeal Holliday junction-resolving enzymes, which are typically divalent metal ion-binding dimers that are able to cleave X-shaped dsDNA–Holliday junctions (Hjs). The crystal structure of Hjc_15-6 was determined to 2.5 Å resolution using the selenomethionine single-wavelength anomalous dispersion method. To our knowledge, this is the first crystal structure of an Hj-resolving enzyme originating from a bacteriophage that can be classified as an archaeal type of Hj-resolving enzyme. As such, it represents a new fold for Hj-resolving enzymes from phages. Characterization of the structure of Hjc_15-6 suggests that it may form a dimer, or even a homodimer of dimers, and activity studies show endonuclease activity towards Hjs. Furthermore, based on sequence analysis it is proposed that Hjc_15-6 has a three-part catalytic motif corresponding to E–SD–EVK, and this motif may be common among other Hj-resolving enzymes originating from thermophilic bacteriophages.

The key players of parthanatos: opportunities for targeting multiple levels in the therapy of parthanatos-based pathogenesis

Parthanatos is a form of regulated cell death involved in the pathogenesis of many diseases, particularly neurodegenerative disorders, such as Parkinson's disease, Alzheimer's disease, Huntington's disease, and amyotrophic lateral sclerosis. Parthanatos is a multistep cell death pathway cascade that involves poly (ADP-ribose) polymerase 1 (PARP-1) overactivation, PAR accumulation, PAR binding to apoptosis-inducing factor (AIF), AIF release from the mitochondria, nuclear translocation of the AIF/macrophage migration inhibitory factor (MIF) complex, and MIF-mediated large-scale DNA fragmentation. All the key players in the parthanatos pathway are pleiotropic proteins with diverse functions. An in-depth understanding of the structure-based activity of the key factors, and the biochemical mechanisms of parthanatos, is crucial for the development of drugs and therapeutic strategies. In this review, we delve into the key players of the parthanatos pathway and reveal the multiple levels of therapeutic opportunities for treating parthanatos-based pathogenesis.

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

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

Improving mobilization of foreign DNA into Zymomonas mobilis ZM4 by removal of multiple restriction systems

Zymomonas mobilis has emerged as a promising candidate for production of high-value bioproducts from plant biomass. However, a major limitation in equipping Z. mobilis with novel pathways to achieve this goal is restriction of heterologous DNA. Here, we characterized the contribution of several defense systems of Z. mobilis strain ZM4 to impeding heterologous gene transfer from an Escherichia coli donor. Bioinformatic analysis revealed that Z. mobilis ZM4 encodes a previously described mrr-like type IV restriction modification (RM) system, a type I-F CRISPR system, a chromosomal type I RM system (hsdMS_c), and a previously uncharacterized type I RM system, located on an endogenous plasmid (hsdRMS_p). The DNA recognition motif of HsdRMS_p was identified by comparing the methylated DNA sequence pattern of mutants lacking one or both of the hsdMS_c and hsdRMS_p systems to that of the parent strain. The conjugation efficiency of synthetic plasmids containing single or combinations of the HsdMS_c and HsdRMS_p recognition sites indicated that both systems are active and decrease uptake of foreign DNA. In contrast, deletions of mrr and cas3 led to no detectable improvement in conjugation efficiency for the exogenous DNA tested. Thus, the suite of markerless restriction-negative strains that we constructed and the knowledge of this new restriction system and its DNA recognition motif provide the necessary platform to flexibly engineer the next generation of Z. mobilis strains for synthesis of valuable products.

IMPORTANCEZymomonas mobilis is equipped with a number of traits that make it a desirable platform organism for metabolic engineering to produce valuable bioproducts. Engineering strains equipped with synthetic pathways for biosynthesis of new molecules requires integration of foreign genes. In this study, we developed an all-purpose strain, devoid of known host restriction systems and free of any antibiotic resistance markers, which dramatically improves the uptake efficiency of heterologous DNA into Z. mobilis ZM4. We also confirmed the role of a previously known restriction system as well as identifying a previously unknown type I RM system on an endogenous plasmid. Elimination of the barriers to DNA uptake as shown here will allow facile genetic engineering of Z. mobilis.

Comparative Genomics Reveals Factors Associated with Phenotypic Expression of Wolbachia

Wolbachia is a widespread, vertically transmitted bacterial endosymbiont known for manipulating arthropod reproduction. Its most common form of reproductive manipulation is cytoplasmic incompatibility (CI), observed when a modification in the male sperm leads to embryonic lethality unless a compatible rescue factor is present in the female egg. CI attracts scientific attention due to its implications for host speciation and in the use of Wolbachia for controlling vector-borne diseases. However, our understanding of CI is complicated by the complexity of the phenotype, whose expression depends on both symbiont and host factors. In the present study, we perform a comparative analysis of nine complete Wolbachia genomes with known CI properties in the same genetic host background, Drosophila simulans STC. We describe genetic differences between closely related strains and uncover evidence that phages and other mobile elements contribute to the rapid evolution of both genomes and phenotypes of Wolbachia. Additionally, we identify both known and novel genes associated with the modification and rescue functions of CI. We combine our observations with published phenotypic information and discuss how variability in cif genes, novel CI-associated genes, and Wolbachia titer might contribute to poorly understood aspects of CI such as strength and bidirectional incompatibility. We speculate that high titer CI strains could be better at invading new hosts already infected with a CI Wolbachia, due to a higher rescue potential, and suggest that titer might thus be a relevant parameter to consider for future strategies using CI Wolbachia in biological control.

The Card1 nuclease provides defence during type III CRISPR immunity

In the type III CRISPR-Cas immune response of prokaryotes, infection triggers the production of cyclic oligoadenylates that bind and activate proteins that contain a CARF domain^1,2. Many type III loci are associated with proteins in which the CRISPR-associated Rossman fold (CARF) domain is fused to a restriction  endonuclease-like domain^3,4. However, with the exception of the well-characterized Csm6 and Csx1 ribonucleases^5,6, whether and how these inducible effectors provide defence is not known. Here we investigated a type III CRISPR accessory protein, which we name cyclic-oligoadenylate-activated single-stranded ribonuclease and single-stranded deoxyribonuclease 1 (Card1). Card1 forms a symmetrical dimer that has a large central cavity between its CRISPR-associated Rossmann fold and restriction endonuclease domains that binds cyclic tetra-adenylate. The binding of ligand results in a conformational change comprising the rotation of individual monomers relative to each other to form a more compact dimeric scaffold, in which a manganese cation coordinates the catalytic residues and activates the cleavage of single-stranded-but not double-stranded-nucleic acids (both DNA and RNA). In vivo, activation of Card1 induces dormancy of the infected hosts to provide immunity against phage infection and plasmids. Our results highlight the diversity of strategies used in CRISPR systems to provide immunity.

Active site residue identity regulates cleavage preference of LAGLIDADG homing endonucleases

LAGLIDADG homing endonucleases (meganucleases) are site-specific mobile endonucleases that can be adapted for genome-editing applications. However, one problem when reprogramming meganucleases on non-native substrates is indirect readout of DNA shape and flexibility at the central 4 bases where cleavage occurs. To understand how the meganuclease active site regulates DNA cleavage, we used functional selections and deep sequencing to profile the fitness landscape of 1600 I-LtrI and I-OnuI active site variants individually challenged with 67 substrates with central 4 base substitutions. The wild-type active site was not optimal for cleavage on many substrates, including the native I-LtrI and I-OnuI targets. Novel combinations of active site residues not observed in known meganucleases supported activity on substrates poorly cleaved by the wild-type enzymes. Strikingly, combinations of E or D substitutions in the two metal-binding residues greatly influenced cleavage activity, and E184D variants had a broadened cleavage profile. Analyses of I-LtrI E184D and the wild-type proteins co-crystallized with the non-cognate AACC central 4 sequence revealed structural differences that correlated with kinetic constants for cleavage of individual DNA strands. Optimizing meganuclease active sites to enhance cleavage of non-native central 4 target sites is a straightforward addition to engineering workflows that will expand genome-editing applications.

A universal fluorescence-based toolkit for real-time quantification of DNA and RNA nuclease activity

DNA and RNA nucleases play a critical role in a growing number of cellular processes ranging from DNA repair to immune surveillance. Nevertheless, many nucleases have unknown or poorly characterized activities. Elucidating nuclease substrate specificities and co-factors can support a more definitive understanding of cellular mechanisms in physiology and disease. Using fluorescence-based methods, we present a quick, safe, cost-effective, and real-time versatile nuclease assay, which uniquely studies nuclease enzyme kinetics. In conjunction with a substrate library we can now analyse nuclease catalytic rates, directionality, and substrate preferences. The assay is sensitive enough to detect kinetics of repair enzymes when confronted with DNA mismatches or DNA methylation sites. We have also extended our analysis to study the kinetics of human single-strand DNA nuclease TREX2, DNA polymerases, RNA, and RNA:DNA nucleases. These nucleases are involved in DNA repair, immune regulation, and have been associated with various diseases, including cancer and immune disorders.

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

Evolutionary Genetics of Cytoplasmic Incompatibility Genes cifA and cifB in Prophage WO of Wolbachia

The bacterial endosymbiont Wolbachia manipulates arthropod reproduction to facilitate its maternal spread through host populations. The most common manipulation is cytoplasmic incompatibility (CI): Wolbachia-infected males produce modified sperm that cause embryonic mortality, unless rescued by embryos harboring the same Wolbachia. The genes underlying CI, cifA and cifB, were recently identified in the eukaryotic association module of Wolbachia’s prophage WO. Here, we use transcriptomic and genomic approaches to address three important evolutionary facets of the cif genes. First, we assess whether or not cifA and cifB comprise a classic toxin–antitoxin operon in wMel and show that the two genes exhibit striking, transcriptional differences across host development. They can produce a bicistronic message despite a predicted hairpin termination element in their intergenic region. Second, cifA and cifB strongly coevolve across the diversity of phage WO. Third, we provide new domain and functional predictions across homologs within Wolbachia, and show that amino acid sequences vary substantially across the genus. Finally, we investigate conservation of cifA and cifB and find frequent degradation and loss of the genes in strains that no longer induce CI. Taken together, we demonstrate that cifA and cifB exhibit complex transcriptional regulation in wMel, provide functional annotations that broaden the potential mechanisms of CI induction, and report recurrent erosion of cifA and cifB in non-CI strains, thus expanding our understanding of the most widespread form of reproductive parasitism.

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

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

A nuclease that mediates cell death induced by DNA damage and poly(ADP-ribose) polymerase-1

Inhibition or genetic deletion of poly(ADP-ribose) (PAR) polymerase-1 (PARP-1) is protective against toxic insults in many organ systems. The molecular mechanisms underlying PARP-1-dependent cell death involve release of mitochondrial apoptosis-inducing factor (AIF) and its translocation to the nucleus, which results in chromatinolysis. We identified macrophage migration inhibitory factor (MIF) as a PARP-1-dependent AIF-associated nuclease (PAAN). AIF was required for recruitment of MIF to the nucleus, where MIF cleaves genomic DNA into large fragments. Depletion of MIF, disruption of the AIF-MIF interaction, or mutation of glutamic acid at position 22 in the catalytic nuclease domain blocked MIF nuclease activity and inhibited chromatinolysis, cell death induced by glutamate excitotoxicity, and focal stroke. Inhibition of MIF's nuclease activity is a potential therapeutic target for diseases caused by excessive PARP-1 activation.

Functional Diversity of Cytotoxic tRNase/Immunity Protein Complexes from Burkholderia pseudomallei

Contact-dependent growth inhibition (CDI) is a widespread mechanism of inter-bacterial competition. CDI(+) bacteria deploy large CdiA effector proteins, which carry variable C-terminal toxin domains (CdiA-CT). CDI(+) cells also produce CdiI immunity proteins that specifically neutralize cognate CdiA-CT toxins to prevent auto-inhibition. Here, we present the crystal structure of the CdiA-CT/CdiI(E479) toxin/immunity protein complex from Burkholderia pseudomallei isolate E479. The CdiA-CT(E479) tRNase domain contains a core α/β-fold that is characteristic of PD(D/E)XK superfamily nucleases. Unexpectedly, the closest structural homolog of CdiA-CT(E479) is another CDI toxin domain from B. pseudomallei 1026b. Although unrelated in sequence, the two B. pseudomallei nuclease domains share similar folds and active-site architectures. By contrast, the CdiI(E479) and CdiI(1026b) immunity proteins share no significant sequence or structural homology. CdiA-CT(E479) and CdiA-CT(1026b) are both tRNases; however, each nuclease cleaves tRNA at a distinct position. We used a molecular docking approach to model each toxin bound to tRNA substrate. The resulting models fit into electron density envelopes generated by small-angle x-ray scattering analysis of catalytically inactive toxin domains bound stably to tRNA. CdiA-CT(E479) is the third CDI toxin found to have structural homology to the PD(D/E)XK superfamily. We propose that CDI systems exploit the inherent sequence variability and active-site plasticity of PD(D/E)XK nucleases to generate toxin diversity. These findings raise the possibility that many other uncharacterized CDI toxins may belong to the PD(D/E)XK superfamily.

Sin Nombre hantavirus nucleocapsid protein exhibits a metal-dependent DNA-specific endonucleolytic activity

We demonstrate that the nucleocapsid protein of Sin Nombre hantavirus (SNV-N) has a DNA-specific endonuclease activity. Upon incubation of SNV-N with DNA in the presence of magnesium or manganese, we observed DNA digestion in sequence-unspecific manner. In contrast, RNA was not affected under the same conditions. Moreover, pre-treatment of SNV-N with RNase before DNA cleavage increased the endonucleolytic activity. Structure-based protein fold prediction using known structures from the PDB database revealed that Asp residues in positions 88 and 103 of SNV-N show sequence similarity with the active site of the restriction endonuclease HindIII. Crystal structure of HindIII predicts that residues Asp93 and Asp108 are essential for coordination of the metal ions required for HindIII DNA cleavage. Therefore, we hypothesized that homologous residues in SNV-N, Asp88 and Asp103, may have a similar function. Replacing Asp88 and Asp103 by alanine led to an SNV-N protein almost completely abrogated for endonuclease activity.

Endonuclease domain of non-LTR retrotransposons: loss-of-function mutants and modeling of the R2Bm endonuclease

Non-LTR retrotransposons are an important class of mobile elements that insert into host DNA by target-primed reverse transcription (TPRT). Non-LTR retrotransposons must bind to their mRNA, recognize and cleave their target DNA, and perform TPRT at the site of DNA cleavage. As DNA binding and cleavage are such central parts of the integration reaction, a better understanding of the endonuclease encoded by non-LTR retrotransposons is needed. This paper explores the R2 endonuclease domain from Bombyx mori using in vitro studies and in silico modeling. Mutations in conserved sequences located across the putative PD-(D/E)XK endonuclease domain reduced DNA cleavage, DNA binding and TPRT. A mutation at the beginning of the first α-helix of the modeled endonuclease obliterated DNA cleavage and greatly reduced DNA binding. It also reduced TPRT when tested on pre-cleaved DNA substrates. The catalytic K was located to a non-canonical position within the second α-helix. A mutation located after the fourth β-strand reduced DNA binding and cleavage. The motifs that showed impaired activity form an extensive basic region. The R2 biochemical and structural data are compared and contrasted with that of two other well characterized PD-(D/E)XK endonucleases, restriction endonucleases and archaeal Holliday junction resolvases.

Novel recA-Independent Horizontal Gene Transfer in Escherichia coli K-12

In bacteria, mechanisms that incorporate DNA into a genome without strand-transfer proteins such as RecA play a major role in generating novelty by horizontal gene transfer. We describe a new illegitimate recombination event in Escherichia coli K-12: RecA-independent homologous replacements, with very large (megabase-length) donor patches replacing recipient DNA. A previously uncharacterized gene (yjiP) increases the frequency of RecA-independent replacement recombination. To show this, we used conjugal DNA transfer, combining a classical conjugation donor, HfrH, with modern genome engineering methods and whole genome sequencing analysis to enable interrogation of genetic dependence of integration mechanisms and characterization of recombination products. As in classical experiments, genomic DNA transfer begins at a unique position in the donor, entering the recipient via conjugation; antibiotic resistance markers are then used to select recombinant progeny. Different configurations of this system were used to compare known mechanisms for stable DNA incorporation, including homologous recombination, F'-plasmid formation, and genome duplication. A genome island of interest known as the immigration control region was specifically replaced in a minority of recombinants, at a frequency of 3 X 10(-12) CFU/recipient per hour.

Genome-wide DNA methylation analysis of Haloferax volcanii H26 and identification of DNA methyltransferase related PD-(D/E)XK nuclease family protein HVO_A0006

Restriction-modification (RM) systems have evolved to protect the cell from invading DNAs and are composed of two enzymes: a DNA methyltransferase and a restriction endonuclease. Although RM systems are present in both archaeal and bacterial genomes, DNA methylation in archaea has not been well defined. In order to characterize the function of RM systems in archaeal species, we have made use of the model haloarchaeon Haloferax volcanii. A genomic DNA methylation analysis of H. volcanii strain H26 was performed using PacBio single molecule real-time (SMRT) sequencing. This analysis was also performed on a strain of H. volcanii in which an annotated DNA methyltransferase gene HVO_A0006 was deleted from the genome. Sequence analysis of H26 revealed two motifs which are modified in the genome: C(m4)TAG and GCA(m6)BN6VTGC. Analysis of the ΔHVO_A0006 strain indicated that it exhibited reduced adenine methylation compared to the parental strain and altered the detected adenine motif. However, protein domain architecture analysis and amino acid alignments revealed that HVO_A0006 is homologous only to the N-terminal endonuclease region of Type IIG RM proteins and contains a PD-(D/E)XK nuclease motif, suggesting that HVO_A0006 is a PD-(D/E)XK nuclease family protein. Further bioinformatic analysis of the HVO_A0006 gene demonstrated that the gene is rare among the Halobacteria. It is surrounded by two transposition genes suggesting that HVO_A0006 is a fragment of a Type IIG RM gene, which has likely been acquired through gene transfer, and affects restriction-modification activity by interacting with another RM system component(s). Here, we present the first genome-wide characterization of DNA methylation in an archaeal species and examine the function of a DNA methyltransferase related gene HVO_A0006.

Crystal structure of a Fanconi anemia-associated nuclease homolog bound to 5' flap DNA: basis of interstrand cross-link repair by FAN1

Fanconi anemia (FA) is an autosomal recessive genetic disorder caused by defects in FA genes responsible for processing DNA interstrand cross-links (ICLs). FA-associated nuclease (FAN1) is recruited to lesions by a monoubiquitinated FANCI–FANCD2 (ID) complex and participates in ICL repair. Here, Gwon et al. determined the crystal structure of Pseudomonas aeruginosa FAN1 (PaFAN1) lacking the UBZ (ubiquitin-binding zinc) domain in complex with 5′ flap DNA. The PaFAN1 structure provides insights into how FAN1 integrates with the FA complex to participate in ICL repair.

Fanconi anemia (FA) is an autosomal recessive genetic disorder caused by defects in any of 15 FA genes responsible for processing DNA interstrand cross-links (ICLs). The ultimate outcome of the FA pathway is resolution of cross-links, which requires structure-selective nucleases. FA-associated nuclease 1 (FAN1) is believed to be recruited to lesions by a monoubiquitinated FANCI–FANCD2 (ID) complex and participates in ICL repair. Here, we determined the crystal structure of Pseudomonas aeruginosa FAN1 (PaFAN1) lacking the UBZ (ubiquitin-binding zinc) domain in complex with 5′ flap DNA. All four domains of the right-hand-shaped PaFAN1 are involved in DNA recognition, with each domain playing a specific role in bending DNA at the nick. The six-helix bundle that binds the junction connects to the catalytic viral replication and repair (VRR) nuclease (VRR nuc) domain, enabling FAN1 to incise the scissile phosphate a few bases distant from the junction. The six-helix bundle also inhibits the cleavage of intact Holliday junctions. PaFAN1 shares several conserved features with other flap structure-selective nucleases despite structural differences. A clamping motion of the domains around the wedge helix, which acts as a pivot, facilitates nucleolytic cleavage. The PaFAN1 structure provides insights into how archaeal Holliday junction resolvases evolved to incise 5′ flap substrates and how FAN1 integrates with the FA complex to participate in ICL repair.

Genomic diversification in strains of Rickettsia felis Isolated from different arthropods

Rickettsia felis (Alphaproteobacteria: Rickettsiales) is the causative agent of an emerging flea-borne rickettsiosis with worldwide occurrence. Originally described from the cat flea, Ctenocephalides felis, recent reports have identified R. felis from other flea species, as well as other insects and ticks. This diverse host range for R. felis may indicate an underlying genetic variability associated with host-specific strains. Accordingly, to determine a potential genetic basis for host specialization, we sequenced the genome of R. felis str. LSU-Lb, which is an obligate mutualist of the parthenogenic booklouse Liposcelis bostrychophila (Insecta: Psocoptera). We also sequenced the genome of R. felis str. LSU, the second genome sequence for cat flea-associated strains (cf. R. felis str. URRWXCal2), which are presumably facultative parasites of fleas. Phylogenomics analysis revealed R. felis str. LSU-Lb diverged from the flea-associated strains. Unexpectedly, R. felis str. LSU was found to be divergent from R. felis str. URRWXCal2, despite sharing similar hosts. Although all three R. felis genomes contain the pRF plasmid, R. felis str. LSU-Lb carries an additional unique plasmid, pLbaR (plasmid of L. bostrychophila associated Rickettsia), nearly half of which encodes a unique 23-gene integrative conjugative element. Remarkably, pLbaR also encodes a repeats-in-toxin-like type I secretion system and associated toxin, heretofore unknown from other Rickettsiales genomes, which likely originated from lateral gene transfer with another obligate intracellular parasite of arthropods, Cardinium (Bacteroidetes). Collectively, our study reveals unexpected genomic diversity across three R. felis strains and identifies several diversifying factors that differentiate facultative parasites of fleas from obligate mutualists of booklice.

Distinct facilitated diffusion mechanisms by E. coli Type II restriction endonucleases

The passive search by proteins for particular DNA sequences involving nonspecific DNA is essential for gene regulation, DNA repair, phage defense, and diverse epigenetic processes. Distinct mechanisms contribute to these searches, and it remains unresolved as to which mechanism or blend of mechanisms best suits a particular protein and, more importantly, its biological role. To address this, we compare the translocation properties of two well-studied bacterial restriction endonucleases (ENases), EcoRI and EcoRV. These dimeric, magnesium-dependent enzymes hydrolyze related sites (EcoRI ENase, 5'-GAATTC-3'; EcoRV ENase, 5'-GATATC-3'), leaving overhangs and blunt DNA segments, respectively. Here, we demonstrate that the extensive sliding by EcoRI ENase, involving sliding up to ∼600 bp prior to dissociating from the DNA, contrasts with a larger reliance on hopping mechanism(s) by EcoRV ENase. The mechanism displayed by EcoRI ENase results in a highly thorough search of DNA, whereas the EcoRV ENase mechanism results in an extended, yet less rigorous, interrogation of DNA sequence space. We describe how these mechanistic distinctions are complemented by other aspects of these endonucleases, such as the 10-fold higher in vivo concentrations of EcoRI ENase compared to that of EcoRV ENase. Further, we hypothesize that the highly diverse enzyme arsenal that bacteria employ against foreign DNA involves seemingly similar enzymes that rely on distinct but complementary search mechanisms. Our comparative approach reveals how different proteins utilize distinct site-locating strategies.

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.

Structural characterization of MepB from Staphylococcus aureus reveals homology to endonucleases

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

Cofactor analogue-induced chemical reactivation of endonuclease activity in a DNA cleavage/methylation deficient TspGWI N₄₇₃A variant in the NPPY motif

We reported previously that TspGWI, a prototype enzyme of a new Thermus sp. family of restriction endonucleases-methyltransferases (REases-MTases), undergoes the novel phenomenon of sinefungin (SIN)-caused specificity transition. Here we investigated mutant TspGWI N₄₇₃A, containing a single amino acid (aa) substitution in the NPPY motif of the MTase. Even though the aa substitution is located within the MTase polypeptide segment, DNA cleavage and modification are almost completely abolished, indicating that the REase and MTase are intertwined. Remarkably, the TspGWI N₄₇₃A REase functionality can be completely reconstituted by the addition of SIN. We hypothesize that SIN binds specifically to the enzyme and restores the DNA cleavage-competent protein tertiary structure. This indicates the significant role of allosteric effectors in DNA cleavage in Thermus sp. enzymes. This is the first case of REase mutation suppression by an S-adenosylmethionine (SAM) cofactor analogue. Moreover, the TspGWI N₄₇₃A clone strongly affects E. coli division control, acting as a ‘selfish gene’. The mutant lacks the competing MTase activity and therefore might be useful for applications in DNA manipulation. Here we present a case study of a novel strategy for REase activity/specificity alteration by a single aa substitution, based on the bioinformatic analysis of active motif locations, combining (a) aa sequence engineering (b) the alteration of protein enzymatic properties, and (c) the use of cofactor–analogue cleavage reconstitution and stimulation.

Structure and cleavage activity of the tetrameric MspJI DNA modification-dependent restriction endonuclease

The MspJI modification-dependent restriction endonuclease recognizes 5-methylcytosine or 5-hydroxymethylcytosine in the context of CNN(G/A) and cleaves both strands at fixed distances (N₁₂/N₁₆) away from the modified cytosine at the 3′-side. We determined the crystal structure of MspJI of Mycobacterium sp. JLS at 2.05-Å resolution. Each protein monomer harbors two domains: an N-terminal DNA-binding domain and a C-terminal endonuclease. The N-terminal domain is structurally similar to that of the eukaryotic SET and RING-associated domain, which is known to bind to a hemi-methylated CpG dinucleotide. Four protein monomers are found in the crystallographic asymmetric unit. Analytical gel-filtration and ultracentrifugation measurements confirm that the protein exists as a tetramer in solution. Two monomers form a back-to-back dimer mediated by their C-terminal endonuclease domains. Two back-to-back dimers interact to generate a tetramer with two double-stranded DNA cleavage modules. Each cleavage module contains two active sites facing each other, enabling double-strand DNA cuts. Biochemical, mutagenesis and structural characterization suggest three different monomers of the tetramer may be involved respectively in binding the modified cytosine, making the first proximal N₁₂ cleavage in the same strand and then the second distal N₁₆ cleavage in the opposite strand. Both cleavage events require binding of at least a second recognition site either in cis or in trans.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

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

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

Related bifunctional restriction endonuclease-methyltransferase triplets: TspDTI, Tth111II/TthHB27I and TsoI with distinct specificities

BACKGROUND

We previously defined a family of restriction endonucleases (REases) from Thermus sp., which share common biochemical and biophysical features, such as the fusion of both the nuclease and methyltransferase (MTase) activities in a single polypeptide, cleavage at a distance from the recognition site, large molecular size, modulation of activity by S-adenosylmethionine (SAM), and incomplete cleavage of the substrate DNA. Members include related thermophilic REases with five distinct specificities: TspGWI, TaqII, Tth111II/TthHB27I, TspDTI and TsoI.

RESULTS

TspDTI, TsoI and isoschizomers Tth111II/TthHB27I recognize different, but related sequences: 5'-ATGAA-3', 5'-TARCCA-3' and 5'-CAARCA-3' respectively. Their amino acid sequences are similar, which is unusual among REases of different specificity. To gain insight into this group of REases, TspDTI, the prototype member of the Thermus sp. enzyme family, was cloned and characterized using a recently developed method for partially cleaving REases.

CONCLUSIONS

TspDTI, TsoI and isoschizomers Tth111II/TthHB27I are closely related bifunctional enzymes. They comprise a tandem arrangement of Type I-like domains, like other Type IIC enzymes (those with a fusion of a REase and MTase domains), e.g. TspGWI, TaqII and MmeI, but their sequences are only remotely similar to these previously characterized enzymes. The characterization of TspDTI, a prototype member of this group, extends our understanding of sequence-function relationships among multifunctional restriction-modification enzymes.

Restriction endonucleases: natural and directed evolution

Type II restriction endonucleases (REs) are highly sequence-specific compared with other classes of nucleases. PD-(D/E)XK nucleases, initially represented by only type II REs, now comprise a large and extremely diverse superfamily of proteins and, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. Sequence similarity can only be observed in methylases and few isoschizomers. As a consequence, REs are classified according to combinations of functional properties rather than on the basis of genetic relatedness. New alignment matrices and classification systems based on structural core connectivity and cleavage mechanisms have been developed to characterize new REs and related proteins. REs recognizing more than 300 distinct specificities have been identified in RE database (REBASE: http://rebase.neb.com/cgi-bin/statlist ) but still the need for newer specificities is increasing due to the advancement in molecular biology and applications. The enzymes have undergone constant evolution through structural changes in protein scaffolds which include random mutations, homologous recombinations, insertions, and deletions of coding DNA sequences but rational mutagenesis or directed evolution delivers protein variants with new functions in accordance with defined biochemical or environmental pressures. Redesigning through random mutation, addition or deletion of amino acids, methylation-based selection, synthetic molecules, combining recognition and cleavage domains from different enzymes, or combination with domains of additional functions change the cleavage specificity or substrate preference and stability. There is a growing number of patents awarded for the creation of engineered REs with new and enhanced properties.

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

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

Viral RNase Involvement in Strategies of Infection

The overwhelming majority of RNase activity is engaged in catabolic processes. Viruses have no metabolism of their own, but rely completely on host cellular energy and substrate provision to support the biochemical processes necessary for virus replication. It is therefore obvious that RNA hydrolysis does not represent an obligate step in the viral life cycle that would have to be governed by viral proteins. Accordingly, RNases are found only rarely in the viral proteomes and serve special functions. In this chapter, several virus-specific RNases will be described and their role in the viral life cycle discussed. The text will concentrate on RNases of members of the nidoviruses, herpesviruses, pestiviruses, and several viruses with segmented negative-strand RNA genome including influenza virus. These enzymes are involved in specific steps of viral gene expression, viral genome replication, shutoff of host cellular gene expression, and interference with the host’s immune response to virus infection.

KIAA1018/FAN1 nuclease protects cells against genomic instability induced by interstrand cross-linking agents

Fanconi anemia (FA) is a rare genetic disease characterized by congenital defects, bone marrow failure, chromosomal instability, and cancer susceptibility. One hallmark of cells from FA patients is hypersensitivity to interstrand cross-linking agents, such as the chemotherapeutics cisplatin and mitomycin C (MMC). We have recently characterized a FANCD2/FANCI-associated nuclease, KIAA1018/FAN1, the depletion of which sensitizes human cells to these agents. However, as the down-regulation of FAN1 in human cells was mediated by siRNA and thus only transient, we were unable to study the long-term effects of FAN1 loss on chromosomal stability. We now describe the generation of chicken DT40 B cells, in which the FAN1 locus was disrupted by gene targeting. FAN1-null cells are highly sensitive to cisplatin and MMC, but not to ionizing or UV radiation, methyl methanesulfonate, or camptothecin. The cells do not display elevated sister chromatid exchange frequencies, either sporadic or MMC-induced. Interestingly, MMC treatment causes chromosomal instability that is quantitatively, but not qualitatively, comparable to that seen in FA cells. This finding, coupled with evidence showing that DT40 cells deficient in both FAN1 and FANCC, or FAN1 and FANCJ, exhibited increased sensitivity to cisplatin compared with cells lacking only FAN1, suggests that, despite its association with FANCD2/FANCI, FAN1 in DT40 cells participates in the processing of damage induced by interstrand cross-linking-generating agents also independently of the classical FA pathway.

Identification of new homologs of PD-(D/E)XK nucleases by support vector machines trained on data derived from profile-profile alignments

PD-(D/E)XK nucleases, initially represented by only Type II restriction enzymes, now comprise a large and extremely diverse superfamily of proteins. They participate in many different nucleic acids transactions including DNA degradation, recombination, repair and RNA processing. Different PD-(D/E)XK families, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. This makes the identification of new superfamily members using standard homology search techniques challenging. To tackle this problem, we developed a method for the detection of PD-(D/E)XK families based on the binary classification of profile–profile alignments using support vector machines (SVMs). Using a number of both superfamily-specific and general features, SVMs were trained to identify true positive alignments of PD-(D/E)XK representatives. With this method we identified several PFAM families of uncharacterized proteins as putative new members of the PD-(D/E)XK superfamily. In addition, we assigned several unclassified restriction enzymes to the PD-(D/E)XK type. Results show that the new method is able to make confident assignments even for alignments that have statistically insignificant scores. We also implemented the method as a freely accessible web server at http://www.ibt.lt/bioinformatics/software/pdexk/.

Unusual target site disruption by the rare-cutting HNH restriction endonuclease PacI

The crystal structure of the rare-cutting HNH restriction endonuclease PacI in complex with its eight-base-pair target recognition sequence 5'-TTAATTAA-3' has been determined to 1.9 A resolution. The enzyme forms an extended homodimer, with each subunit containing two zinc-bound motifs surrounding a betabetaalpha-metal catalytic site. The latter is unusual in that a tyrosine residue likely initiates strand cleavage. PacI dramatically distorts its target sequence from Watson-Crick duplex DNA base pairing, with every base separated from its original partner. Two bases on each strand are unpaired, four are engaged in noncanonical A:A and T:T base pairs, and the remaining two bases are matched with new Watson-Crick partners. This represents a highly unusual DNA binding mechanism for a restriction endonuclease, and implies that initial recognition of the target site might involve significantly different contacts from those visualized in the DNA-bound cocrystal structures.

Folding, DNA recognition, and function of GIY-YIG endonucleases: crystal structures of R.Eco29kI

The GIY-YIG endonuclease family comprises hundreds of diverse proteins and a multitude of functions; none have been visualized bound to DNA. The structure of the GIY-YIG restriction endonuclease R.Eco29kI has been solved both alone and bound to its target site. The protein displays a domain-swapped homodimeric structure with several extended surface loops encircling the DNA. Only three side chains from each protein subunit contact DNA bases, two directly and one via a bridging solvent molecule. Both tyrosine residues within the GIY-YIG motif are positioned in the catalytic center near a putative nucleophilic water; the remainder of the active site resembles the HNH endonuclease family. The structure illustrates how the GIY-YIG scaffold has been adapted for the highly specific recognition of a DNA restriction site, in contrast to nonspecific DNA cleavage by GIY-YIG domains in homing endonucleases or structure-specific cleavage by DNA repair enzymes such as UvrC.

Deficiency of FANCD2-associated nuclease KIAA1018/FAN1 sensitizes cells to interstrand crosslinking agents

Cytotoxicity of cisplatin and mitomycin C (MMC) is ascribed largely to their ability to generate interstrand crosslinks (ICLs) in DNA, which block the progression of replication forks. The processing of ICLs requires the Fanconi anemia (FA) pathway, excision repair, and translesion DNA synthesis (TLS). It also requires homologous recombination (HR), which repairs double-strand breaks (DSBs) generated by cleavage of the blocked replication forks. Here we describe KIAA1018, an evolutionarily conserved protein that has an N-terminal ubiquitin-binding zinc finger (UBZ) and a C-terminal nuclease domain. KIAA1018 is a 5'-->3' exonuclease and a structure-specific endonuclease that preferentially incises 5' flaps. Like cells from FA patients, human cells depleted of KIAA1018 are sensitized to ICL-inducing agents and display chromosomal instability. The link of KIAA1018 to the FA pathway is further strengthened by its recruitment to DNA damage through interaction of its UBZ domain with monoubiquitylated FANCD2. We therefore propose to name KIAA1018 FANCD2-associated nuclease, FAN1.

Identification of KIAA1018/FAN1, a DNA repair nuclease recruited to DNA damage by monoubiquitinated FANCD2

DNA interstrand crosslinks (ICLs) are highly toxic because they block the progression of replisomes. The Fanconi Anemia (FA) proteins, encoded by genes that are mutated in FA, are important for repair of ICLs. The FA core complex catalyzes the monoubiquitination of FANCD2, and this event is essential for several steps of ICL repair. However, how monoubiquitination of FANCD2 promotes ICL repair at the molecular level is unknown. Here, we describe a highly conserved protein, KIAA1018/MTMR15/FAN1, that interacts with, and is recruited to sites of DNA damage by, the monoubiquitinated form of FANCD2. FAN1 exhibits endonuclease activity toward 5' flaps and has 5' exonuclease activity, and these activities are mediated by an ancient VRR_nuc domain. Depletion of FAN1 from human cells causes hypersensitivity to ICLs, defects in ICL repair, and genome instability. These data at least partly explain how ubiquitination of FANCD2 promotes DNA repair.