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

Live virus-free or die: coupling of antivirus immunity and programmed suicide or dormancy in prokaryotes

Background

The virus-host arms race is a major theater for evolutionary innovation. Archaea and bacteria have evolved diverse, elaborate antivirus defense systems that function on two general principles: i) immune systems that discriminate self DNA from nonself DNA and specifically destroy the foreign, in particular viral, genomes, whereas the host genome is protected, or ii) programmed cell suicide or dormancy induced by infection.

Presentation of the hypothesis

Almost all genomic loci encoding immunity systems such as CRISPR-Cas, restriction-modification and DNA phosphorothioation also encompass suicide genes, in particular those encoding known and predicted toxin nucleases, which do not appear to be directly involved in immunity. In contrast, the immunity systems do not appear to encode antitoxins found in typical toxin-antitoxin systems. This raises the possibility that components of the immunity system themselves act as reversible inhibitors of the associated toxin proteins or domains as has been demonstrated for the Escherichia coli anticodon nuclease PrrC that interacts with the PrrI restriction-modification system. We hypothesize that coupling of diverse immunity and suicide/dormancy systems in prokaryotes evolved under selective pressure to provide robustness to the antivirus response. We further propose that the involvement of suicide/dormancy systems in the coupled antivirus response could take two distinct forms:

1) induction of a dormancy-like state in the infected cell to ‘buy time’ for activation of adaptive immunity; 2) suicide or dormancy as the final recourse to prevent viral spread triggered by the failure of immunity.

Testing the hypothesis

This hypothesis entails many experimentally testable predictions. Specifically, we predict that Cas2 protein present in all cas operons is a mRNA-cleaving nuclease (interferase) that might be activated at an early stage of virus infection to enable incorporation of virus-specific spacers into the CRISPR locus or to trigger cell suicide when the immune function of CRISPR-Cas systems fails. Similarly, toxin-like activity is predicted for components of numerous other defense loci.

Implications of the hypothesis

The hypothesis implies that antivirus response in prokaryotes involves key decision-making steps at which the cell chooses the path to follow by sensing the course of virus infection.

Reviewers

This article was reviewed by Arcady Mushegian, Etienne Joly and Nick Grishin. For complete reviews, go to the Reviewers’ reports section.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

Rational engineering of sequence specificity in R.MwoI restriction endonuclease

R.MwoI is a Type II restriction endonucleases enzyme (REase), which specifically recognizes a palindromic interrupted DNA sequence 5′-GCNNNNNNNGC-3′ (where N indicates any nucleotide), and hydrolyzes the phosphodiester bond in the DNA between the 7th and 8th base in both strands. R.MwoI exhibits remote sequence similarity to R.BglI, a REase with known structure, which recognizes an interrupted palindromic target 5′-GCCNNNNNGGC-3′. A homology model of R.MwoI in complex with DNA was constructed and used to predict functionally important amino acid residues that were subsequently targeted by mutagenesis. The model, together with the supporting experimental data, revealed regions important for recognition of the common bases in DNA sequences recognized by R.BglI and R.MwoI. Based on the bioinformatics analysis, we designed substitutions of the S310 residue in R.MwoI to arginine or glutamic acid, which led to enzyme variants with altered sequence selectivity compared with the wild-type enzyme. The S310R variant of R.MwoI preferred the 5′-GCCNNNNNGGC-3′ sequence as a target, similarly to R.BglI, whereas the S310E variant preferentially cleaved a subset of the MwoI sites, depending on the identity of the 3rd and 9th nucleotide residues. Our results represent a case study of a REase sequence specificity alteration by a single amino acid substitution, based on a theoretical model in the absence of a crystal structure.

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.

Crystal structure and mechanism of action of the N6-methyladenine-dependent type IIM restriction endonuclease R.DpnI

DNA methylation-dependent restriction enzymes have many applications in genetic engineering and in the analysis of the epigenetic state of eukaryotic genomes. Nevertheless, high-resolution structures have not yet been reported, and therefore mechanisms of DNA methylation-dependent cleavage are not understood. Here, we present a biochemical analysis and high-resolution DNA co-crystal structure of the N⁶-methyladenine (m6A)-dependent restriction enzyme R.DpnI. Our data show that R.DpnI consists of an N-terminal catalytic PD-(D/E)XK domain and a C-terminal winged helix (wH) domain. Surprisingly, both domains bind DNA in a sequence- and methylation-sensitive manner. The crystal contains R.DpnI with fully methylated target DNA bound to the wH domain, but distant from the catalytic domain. Independent readout of DNA sequence and methylation by the two domains might contribute to R.DpnI specificity or could help the monomeric enzyme to cut the second strand after introducing a nick.

Structural, functional and evolutionary relationships between homing endonucleases and proteins from their host organisms

Homing endonucleases (HEs) are highly specific DNA-cleaving enzymes that are encoded by invasive DNA elements (usually mobile introns or inteins) within the genomes of phage, bacteria, archea, protista and eukaryotic organelles. Six unique structural HE families, that collectively span four distinct nuclease catalytic motifs, have been characterized to date. Members of each family display structural homology and functional relationships to a wide variety of proteins from various organisms. The biological functions of those proteins are highly disparate and include non-specific DNA-degradation enzymes, restriction endonucleases, DNA-repair enzymes, resolvases, intron splicing factors and transcription factors. These relationships suggest that modern day HEs share common ancestors with proteins involved in genome fidelity, maintenance and gene expression. This review summarizes the results of structural studies of HEs and corresponding proteins from host organisms that have illustrated the manner in which these factors are related.

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.

Plasmid pP62BP1 isolated from an Arctic Psychrobacter sp. strain carries two highly homologous type II restriction-modification systems and a putative organic sulfate metabolism operon

The complete nucleotide sequence of plasmid pP62BP1 (34,467 bp), isolated from Arctic Psychrobacter sp. DAB_AL62B, was determined and annotated. The conserved plasmid backbone is composed of several genetic modules, including a replication system (REP) with similarities to the REP region of the iteron-containing plasmid pPS10 of Pseudomonas syringae. The additional genetic load of pP62BP1 includes two highly related type II restriction-modification systems and a set of genes (slfRCHSL) encoding enzymes engaged in the metabolism of organic sulfates, plus a putative transcriptional regulator (SlfR) of the AraC family. The pP62BP1 slf locus has a compact and unique structure. It is predicted that the enzymes SlfC, SlfH, SlfS and SlfL carry out a chain of reactions leading to the transformation of alkyl sulfates into acyl-CoA, with dodecyl sulfate (SDS) as a possible starting substrate. Comparative analysis of the nucleotide sequences of pP62BP1 and other Psychrobacter spp. plasmids revealed their structural diversity. However, the presence of a few highly conserved DNA segments in pP62BP1, plasmid 1 of P. cryohalolentis K5 and pRWF-101 of Psychrobacter sp. PRwf-1 is indicative of recombinational shuffling of genetic information, and is evidence of lateral gene transfer in the Arctic environment.

Viral polymerases

Viral polymerases play a central role in viral genome replication and transcription. Based on the genome type and the specific needs of particular virus, RNA-dependent RNA polymerase, RNA-dependent DNA polymerase, DNA-dependent RNA polymerase, and DNA-dependent RNA polymerases are found in various viruses. Viral polymerases are generally active as a single protein capable of carrying out multiple functions related to viral genome synthesis. Specifically, viral polymerases use variety of mechanisms to recognize initial binding sites, ensure processive elongation, terminate replication at the end of the genome, and also coordinate the chemical steps of nucleic acid synthesis with other enzymatic activities. This review focuses on different viral genome replication and transcription strategies, and the polymerase interactions with various viral proteins that are necessary to complete genome synthesis.

Towards artificial metallonucleases for gene therapy: recent advances and new perspectives

The process of DNA targeting or repair of mutated genes within the cell, induced by specifically positioned double-strand cleavage of DNA near the mutated sequence, can be applied for gene therapy of monogenic diseases. For this purpose, highly specific artificial metallonucleases are developed. They are expected to be important future tools of modern genetics. The present state of art and strategies of research are summarized, including protein engineering and artificial ‘chemical’ nucleases. From the results, we learn about the basic role of the metal ions and the various ligands, and about the DNA binding and cleavage mechanism. The results collected provide useful guidance for engineering highly controlled enzymes for use in gene therapy.

Endonuclease active site plasticity allows DNA cleavage with diverse alkaline Earth and transition metal ions

A majority of enzymes show a high degree of specificity toward a particular metal ion in their catalytic reaction. However, Type II restriction endonuclease (REase) R.KpnI, which is the first member of the HNH superfamily of REases, exhibits extraordinary diversity in metal ion dependent DNA cleavage. Several alkaline earth and transition group metal ions induce high fidelity and promiscuous cleavage or inhibition depending upon their concentration. The metal ions having different ionic radii and co-ordination geometries readily replace each other from the enzyme's active site, revealing its plasticity. Ability of R.KpnI to cleave DNA with both alkaline earth and transition group metal ions having varied ionic radii could imply utilization of different catalytic site(s). However, mutation of the invariant His residue of the HNH motif caused abolition of the enzyme activity with all of the cofactors, indicating that the enzyme follows a single metal ion catalytic mechanism for DNA cleavage. Indispensability of His in nucleophile activation together with broad cofactor tolerance of the enzyme indicates electrostatic stabilization function of metal ions during catalysis. Nevertheless, a second metal ion is recruited at higher concentrations to either induce promiscuity or inhibit the DNA cleavage. Regulation of the endonuclease activity and fidelity by a second metal ion binding is a unique feature of R.KpnI among REases and HNH nucleases. The active site plasticity of R.KpnI opens up avenues for redesigning cofactor specificities and generation of mutants specific to a particular metal ion.

Defense islands in bacterial and archaeal genomes and prediction of novel defense systems

The arms race between cellular life forms and viruses is a major driving force of evolution. A substantial fraction of bacterial and archaeal genomes is dedicated to antivirus defense. We analyzed the distribution of defense genes and typical mobilome components (such as viral and transposon genes) in bacterial and archaeal genomes and demonstrated statistically significant clustering of antivirus defense systems and mobile genes and elements in genomic islands. The defense islands are enriched in putative operons and contain numerous overrepresented gene families. A detailed sequence analysis of the proteins encoded by genes in these families shows that many of them are diverged variants of known defense system components, whereas others show features, such as characteristic operonic organization, that are suggestive of novel defense systems. Thus, genomic islands provide abundant material for the experimental study of bacterial and archaeal antivirus defense. Except for the CRISPR-Cas systems, different classes of defense systems, in particular toxin-antitoxin and restriction-modification systems, show nonrandom clustering in defense islands. It remains unclear to what extent these associations reflect functional cooperation between different defense systems and to what extent the islands are genomic "sinks" that accumulate diverse nonessential genes, particularly those acquired via horizontal gene transfer. The characteristics of defense islands resemble those of mobilome islands. Defense and mobilome genes are nonrandomly associated in islands, suggesting nonadaptive evolution of the islands via a preferential attachment-like mechanism underpinned by the addictive properties of defense systems such as toxins-antitoxins and an important role of horizontal mobility in the evolution of these islands.

The phage-host arms race: shaping the evolution of microbes

Bacteria, the most abundant organisms on the planet, are outnumbered by a factor of 10 to 1 by phages that infect them. Faced with the rapid evolution and turnover of phage particles, bacteria have evolved various mechanisms to evade phage infection and killing, leading to an evolutionary arms race. The extensive co-evolution of both phage and host has resulted in considerable diversity on the part of both bacterial and phage defensive and offensive strategies. Here, we discuss the unique and common features of phage resistance mechanisms and their role in global biodiversity. The commonalities between defense mechanisms suggest avenues for the discovery of novel forms of these mechanisms based on their evolutionary traits.

Identification and characterization of CbeI, a novel thermostable restriction enzyme from Caldicellulosiruptor bescii DSM 6725 and a member of a new subfamily of HaeIII-like enzymes

Potent HaeIII-like DNA restriction activity was detected in cell-free extracts of Caldicellulosiruptor bescii DSM 6725 using plasmid DNA isolated from Escherichia coli as substrate. Incubation of the plasmid DNA in vitro with HaeIII methyltransferase protected it from cleavage by HaeIII nuclease as well as cell-free extracts of C. bescii. The gene encoding the putative restriction enzyme was cloned and expressed in E. coli with a His-tag at the C-terminus. The purified protein was 38 kDa as predicted by the 981-bp nucleic acid sequence, was optimally active at temperatures between 75°C and 85°C, and was stable for more than 1 week when stored at 35°C. The cleavage sequence was determined to be 5'-GG/CC-3', indicating that CbeI is an isoschizomer of HaeIII. A search of the C. bescii genome sequence revealed the presence of both a HaeIII-like restriction endonuclease (Athe 2438) and DNA methyltransferase (Athe 2437). Preliminary analysis of other Caldicellulosiruptor species suggested that this restriction/modification activity is widespread in this genus. A phylogenetic analysis based on sequence alignment and conserved motif searches identified features of CbeI distinct from other members of this group and classified CbeI as a member of a novel subfamily of HaeIII-like enzymes.

Domain movement within a gene: a novel evolutionary mechanism for protein diversification

A protein function is carried out by a specific domain localized at a specific position. In the present study, we report that, within a gene, a specific amino acid sequence can move between a certain position and another position. This was discovered when the sequences of restriction-modification systems within the bacterial species Helicobacter pylori were compared. In the specificity subunit of Type I restriction-modification systems, DNA sequence recognition is mediated by target recognition domain 1 (TRD1) and TRD2. To our surprise, several sequences are shared by TRD1 and TRD2 of genes (alleles) at the same locus (chromosomal location); these domains appear to have moved between the two positions. The gene/protein organization can be represented as x-(TRD1)-y-x-(TRD2)-y, where x and y represent repeat sequences. Movement probably occurs by recombination at these flanking DNA repeats. In accordance with this hypothesis, recombination at these repeats also appears to decrease two TRDs into one TRD or increase these two TRDs to three TRDs (TRD1-TRD2-TRD2) and to allow TRD movement between genes even at different loci. Similar movement of domains between TRD1 and TRD2 was observed for the specificity subunit of a Type IIG restriction enzyme. Similar movement of domain between TRD1 and TRD2 was observed for Type I restriction-modification enzyme specificity genes in two more eubacterial species, Streptococcus pyogenes and Mycoplasma agalactiae. Lateral domain movements within a protein, which we have designated DOMO (domain movement), represent novel routes for the diversification of proteins.

Conflicts targeting epigenetic systems and their resolution by cell death: novel concepts for methyl-specific and other restriction systems

Epigenetic modification of genomic DNA by methylation is important for defining the epigenome and the transcriptome in eukaryotes as well as in prokaryotes. In prokaryotes, the DNA methyltransferase genes often vary, are mobile, and are paired with the gene for a restriction enzyme. Decrease in a certain epigenetic methylation may lead to chromosome cleavage by the partner restriction enzyme, leading to eventual cell death. Thus, the pairing of a DNA methyltransferase and a restriction enzyme forces an epigenetic state to be maintained within the genome. Although restriction enzymes were originally discovered for their ability to attack invading DNAs, it may be understood because such DNAs show deviation from this epigenetic status. DNAs with epigenetic methylation, by a methyltransferase linked or unlinked with a restriction enzyme, can also be the target of DNases, such as McrBC of Escherichia coli, which was discovered because of its methyl-specific restriction. McrBC responds to specific genome methylation systems by killing the host bacterial cell through chromosome cleavage. Evolutionary and genomic analysis of McrBC homologues revealed their mobility and wide distribution in prokaryotes similar to restriction–modification systems. These findings support the hypothesis that this family of methyl-specific DNases evolved as mobile elements competing with specific genome methylation systems through host killing. These restriction systems clearly demonstrate the presence of conflicts between epigenetic systems.

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

Hpy188I-DNA pre- and post-cleavage complexes--snapshots of the GIY-YIG nuclease mediated catalysis

The GIY-YIG nuclease domain is present in all kingdoms of life and has diverse functions. It is found in the eukaryotic flap endonuclease and Holliday junction resolvase Slx1–Slx4, the prokaryotic nucleotide excision repair proteins UvrC and Cho, and in proteins of ‘selfish’ genetic elements. Here we present the structures of the ternary pre- and post-cleavage complexes of the type II GIY-YIG restriction endonuclease Hpy188I with DNA and a surrogate or catalytic metal ion, respectively. Our structures suggest that GIY-YIG nucleases catalyze DNA hydrolysis by a single substitution reaction. They are consistent with a previous proposal that a tyrosine residue (which we expect to occur in its phenolate form) acts as a general base for the attacking water molecule. In contrast to the earlier proposal, our data identify the general base with the GIY and not the YIG tyrosine. A conserved glutamate residue (Glu149 provided in trans in Hpy188I) anchors a single metal cation in the active site. This metal ion contacts the phosphate proS oxygen atom and the leaving group 3′-oxygen atom, presumably to facilitate its departure. Taken together, our data reveal striking analogy in the absence of homology between GIY-YIG and ββα-Me nucleases.

Nucleases: diversity of structure, function and mechanism

Nucleases cleave the phosphodiester bonds of nucleic acids and may be endo or exo, DNase or RNase, topoisomerases, recombinases, ribozymes, or RNA splicing enzymes. In this review, I survey nuclease activities with known structures and catalytic machinery and classify them by reaction mechanism and metal-ion dependence and by their biological function ranging from DNA replication, recombination, repair, RNA maturation, processing, interference, to defense, nutrient regeneration or cell death. Several general principles emerge from this analysis. There is little correlation between catalytic mechanism and biological function. A single catalytic mechanism can be adapted in a variety of reactions and biological pathways. Conversely, a single biological process can often be accomplished by multiple tertiary and quaternary folds and by more than one catalytic mechanism. Two-metal-ion-dependent nucleases comprise the largest number of different tertiary folds and mediate the most diverse set of biological functions. Metal-ion-dependent cleavage is exclusively associated with exonucleases producing mononucleotides and endonucleases that cleave double- or single-stranded substrates in helical and base-stacked conformations. All metal-ion-independent RNases generate 2',3'-cyclic phosphate products, and all metal-ion-independent DNases form phospho-protein intermediates. I also find several previously unnoted relationships between different nucleases and shared catalytic configurations.

DNA intercalation without flipping in the specific ThaI-DNA complex

The PD-(D/E)XK type II restriction endonuclease ThaI cuts the target sequence CG/CG with blunt ends. Here, we report the 1.3 Å resolution structure of the enzyme in complex with substrate DNA and a sodium or calcium ion taking the place of a catalytic magnesium ion. The structure identifies Glu54, Asp82 and Lys93 as the active site residues. This agrees with earlier bioinformatic predictions and implies that the PD and (D/E)XK motifs in the sequence are incidental. DNA recognition is very unusual: the two Met47 residues of the ThaI dimer intercalate symmetrically into the CG steps of the target sequence. They approach the DNA from the minor groove side and penetrate the base stack entirely. The DNA accommodates the intercalating residues without nucleotide flipping by a doubling of the CG step rise to twice its usual value, which is accompanied by drastic unwinding. Displacement of the Met47 side chains from the base pair midlines toward the downstream CG steps leads to large and compensating tilts of the first and second CG steps. DNA intercalation by ThaI is unlike intercalation by HincII, HinP1I or proteins that bend or repair DNA.

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.

Sequence-specific cleavage of RNA by Type II restriction enzymes

The ability of 223 Type II restriction endonucleases to hydrolyze RNA–DNA heteroduplex oligonucleotide substrates was assessed. Despite the significant topological and sequence asymmetry introduced when one strand of a DNA duplex is substituted by RNA we find that six restriction enzymes (AvaII, AvrII, BanI, HaeIII, HinfI and TaqI), exclusively of the Type IIP class that recognize palindromic or interrupted-palindromic DNA sequences, catalyze robust and specific cleavage of both RNA and DNA strands of such a substrate. Time-course analyses indicate that some endonucleases hydrolyze phosphodiester bonds in both strands simultaneously whereas others appear to catalyze sequential reactions in which either the DNA or RNA product accumulates more rapidly. Such strand-specific variation in cleavage susceptibility is both significant (up to orders of magnitude difference) and somewhat sequence dependent, notably in relation to the presence or absence of uracil residues in the RNA strand. Hybridization to DNA oligonucleotides that contain endonuclease recognition sites can be used to achieve targeted hydrolysis of extended RNA substrates produced by in vitro transcription. The ability to ‘restrict’ an RNA–DNA hybrid, albeit with a limited number of restriction endonucleases, provides a method whereby individual RNA molecules can be targeted for site-specific cleavage in vitro.

BspRI restriction endonuclease: cloning, expression in Escherichia coli and sequential cleavage mechanism

The GGCC-specific restriction endonuclease BspRI is one of the few Type IIP restriction endonucleases, which were suggested to be a monomer. Amino acid sequence information obtained by Edman sequencing and mass spectrometry analysis was used to clone the gene encoding BspRI. The bspRIR gene is located adjacently to the gene of the cognate modification methyltransferase and encodes a 304 aa protein. Expression of the bspRIR gene in Escherichia coli was dependent on the replacement of the native TTG initiation codon with an ATG codon, explaining previous failures in cloning the gene using functional selection. A plasmid containing a single BspRI recognition site was used to analyze kinetically nicking and second-strand cleavage under steady-state conditions. Cleavage of the supercoiled plasmid went through a relaxed intermediate indicating sequential hydrolysis of the two strands. Results of the kinetic analysis of the first- and second-strand cleavage are consistent with cutting the double-stranded substrate site in two independent binding events. A database search identified eight putative restriction-modification systems in which the predicted endonucleases as well as the methyltransferases share high sequence similarity with the corresponding protein of the BspRI system. BspRI and the related putative restriction endonucleases belong to the PD-(D/E)XK nuclease superfamily.

DNA synapsis through transient tetramerization triggers cleavage by Ecl18kI restriction enzyme

To cut DNA at their target sites, restriction enzymes assemble into different oligomeric structures. The Ecl18kI endonuclease in the crystal is arranged as a tetramer made of two dimers each bound to a DNA copy. However, free in solution Ecl18kI is a dimer. To find out whether the Ecl18kI dimer or tetramer represents the functionally important assembly, we generated mutants aimed at disrupting the putative dimer–dimer interface and analysed the functional properties of Ecl18kI and mutant variants. We show by atomic force microscopy that on two-site DNA, Ecl18kI loops out an intervening DNA fragment and forms a tetramer. Using the tethered particle motion technique, we demonstrate that in solution DNA looping is highly dynamic and involves a transient interaction between the two DNA-bound dimers. Furthermore, we show that Ecl18kI cleaves DNA in the synaptic complex much faster than when acting on a single recognition site. Contrary to Ecl18kI, the tetramerization interface mutant R174A binds DNA as a dimer, shows no DNA looping and is virtually inactive. We conclude that Ecl18kI follows the association model for the synaptic complex assembly in which it binds to the target site as a dimer and then associates into a transient tetrameric form to accomplish the cleavage reaction.

The crystal structure of D212 from sulfolobus spindle-shaped virus ragged hills reveals a new member of the PD-(D/E)XK nuclease superfamily

Structural studies have made significant contributions to our understanding of Sulfolobus spindle-shaped viruses (Fuselloviridae), an important model system for archaeal viruses. Continuing these efforts, we report the structure of D212 from Sulfolobus spindle-shaped virus Ragged Hills. The overall fold and conservation of active site residues place D212 in the PD-(D/E)XK nuclease superfamily. The greatest structural similarity is found to the archaeal Holliday junction cleavage enzymes, strongly suggesting a role in DNA replication, repair, or recombination. Other roles associated with nuclease activity are also considered.

Cell-free protein synthesis for structure determination by X-ray crystallography

Structure determination has been difficult for those proteins that are toxic to the cells and cannot be prepared in a large amount in vivo. These proteins, even when biologically very interesting, tend to be left uncharacterized in the structural genomics projects. Their cell-free synthesis can bypass the toxicity problem. Among the various cell-free systems, the wheat-germ-based system is of special interest due to the following points: (1) Because the gene is placed under a plant translational signal, its toxic expression in a bacterial host is reduced. (2) It has only little codon preference and, especially, little discrimination between methionine and selenomethionine (SeMet), which allows easy preparation of selenomethionylated proteins for crystal structure determination by SAD and MAD methods. (3) Translation is uncoupled from transcription, so that the toxicity of the translation product on DNA and its transcription, if any, can be bypassed. We have shown that the wheat-germ-based cell-free protein synthesis is useful for X-ray crystallography of one of the 4-bp cutter restriction enzymes, which are expected to be very toxic to all forms of cells retaining the genome. Our report on its structure represents the first report of structure determination by X-ray crystallography using protein overexpressed with the wheat-germ-based cell-free protein expression system. This will be a method of choice for cytotoxic proteins when its cost is not a problem. Its use will become popular when the crystal structure determination technology has evolved to require only a tiny amount of protein.

Cofactor requirement of HpyAV restriction endonuclease

Background

Helicobacter pylori is the etiologic agent of common gastritis and a risk factor for gastric cancer. It is also one of the richest sources of Type II restriction-modification (R-M) systems in microorganisms.

Principal Findings

We have cloned, expressed and purified a new restriction endonuclease HpyAV from H. pylori strain 26695. We determined the HpyAV DNA recognition sequence and cleavage site as CCTTC 6/5. In addition, we found that HpyAV has a unique metal ion requirement: its cleavage activity is higher with transition metal ions than in Mg⁺⁺. The special metal ion requirement of HpyAV can be attributed to the presence of a HNH catalytic site similar to ColE9 nuclease instead of the canonical PD-X-D/EXK catalytic site found in many other REases. Site-directed mutagenesis was carried out to verify the catalytic residues of HpyAV. Mutation of the conserved metal-binding Asn311 and His320 to alanine eliminated cleavage activity. HpyAV variant H295A displayed approximately 1% of wt activity.

Conclusions/Significance

Some HNH-type endonucleases have unique metal ion cofactor requirement for optimal activities. Homology modeling and site-directed mutagenesis confirmed that HpyAV is a member of the HNH nuclease family. The identification of catalytic residues in HpyAV paved the way for further engineering of the metal binding site. A survey of sequenced microbial genomes uncovered 10 putative R-M systems that show high sequence similarity to the HpyAV system, suggesting lateral transfer of a prototypic HpyAV-like R-M system among these microorganisms.

An Mrr-family nuclease motif in the single polypeptide restriction-modification enzyme LlaGI

Bioinformatic analysis of the putative nuclease domain of the single polypeptide restriction–modification enzyme LlaGI reveals amino acid motifs characteristic of the Escherichia coli methylated DNA-specific Mrr endonuclease. Using mutagenesis, we examined the role of the conserved residues in both DNA translocation and cleavage. Mutations in those residues predicted to play a role in DNA hydrolysis produced enzymes that could translocate on DNA but were either unable to cleave the polynucleotide track or had reduced nuclease activity. Cleavage by LlaGI is not targeted to methylated DNA, suggesting that the conserved motifs in the Mrr domain are a conventional sub-family of the PD-(D/E)XK superfamily of DNA nucleases.

A putative mobile genetic element carrying a novel type IIF restriction-modification system (PluTI)

Genome comparison and genome context analysis were used to find a putative mobile element in the genome of Photorhabdus luminescens, an entomopathogenic bacterium. The element is composed of 16-bp direct repeats in the terminal regions, which are identical to a part of insertion sequences (ISs), a DNA methyltransferase gene homolog, two genes of unknown functions and an open reading frame (ORF) (plu0599) encoding a protein with no detectable sequence similarity to any known protein. The ORF (plu0599) product showed DNA endonuclease activity, when expressed in a cell-free expression system. Subsequently, the protein, named R.PluTI, was expressed in vivo, purified and found to be a novel type IIF restriction enzyme that recognizes 5′-GGCGC/C-3′ (/ indicates position of cleavage). R.PluTI cleaves a two-site supercoiled substrate at both the sites faster than a one-site supercoiled substrate. The modification enzyme homolog encoded by plu0600, named M.PluTI, was expressed in Escherichia coli and shown to protect DNA from R.PluTI cleavage in vitro, and to suppress the lethal effects of R.PluTI expression in vivo. These results suggested that they constitute a restriction–modification system, present on the putative mobile element. Our approach thus allowed detection of a previously uncharacterized family of DNA-interacting proteins.

HHsvm: fast and accurate classification of profile-profile matches identified by HHsearch

MOTIVATION

Recently developed profile-profile methods rival structural comparisons in their ability to detect homology between distantly related proteins. Despite this tremendous progress, many genuine relationships between protein families cannot be recognized as comparisons of their profiles result in scores that are statistically insignificant.

RESULTS

Using known evolutionary relationships among protein superfamilies in SCOP database, support vector machines were trained on four sets of discriminatory features derived from the output of HHsearch. Upon validation, it was shown that the automatic classification of all profile-profile matches was superior to fixed threshold-based annotation in terms of sensitivity and specificity. The effectiveness of this approach was demonstrated by annotating several domains of unknown function from the Pfam database.

AVAILABILITY

Programs and scripts implementing the methods described in this manuscript are freely available from http://hhsvm.dlakiclab.org/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Engineering BspQI nicking enzymes and application of N.BspQI in DNA labeling and production of single-strand DNA

BspQI is a thermostable Type IIS restriction endonuclease (REase) with the recognition sequence 5'GCTCTTC N1/N4 3'. Here we report the cloning and expression of the bspQIR gene for the BspQI restriction enzyme in Escherichia coli. Alanine scanning of the BspQI charged residues identified a number of DNA nicking variants. After sampling combinations of different amino acid substitutions, an Nt.BspQI triple mutant (E172A/E248A/E255K) was constructed with predominantly top-strand DNA nicking activity. Furthermore, a triple mutant of BspQI (Nb.BspQI, N235A/K331A/R428A) was engineered to create a bottom-strand nicking enzyme. In addition, we demonstrated the application of Nt.BspQI in optical mapping of single DNA molecules. Nt or Nb.BspQI-nicked dsDNA can be further digested by E. coli exonuclease III to create ssDNA for downstream applications. BspQI contains two potential catalytic sites: a top-strand catalytic site (Ct) with a D-H-N-K motif found in the HNH endonuclease family and a bottom-strand catalytic site (Cb) with three scattered Glu residues. BlastP analysis of proteins in GenBank indicated a putative restriction enzyme with significant amino acid sequence identity to BspQI from the sequenced bacterial genome Croceibacter atlanticus HTCC2559. This restriction gene was amplified by PCR and cloned into a T7 expression vector. Restriction mapping and run-off DNA sequencing of digested products from the partially purified enzyme indicated that it is an EarI isoschizomer with 6-bp recognition, which we named CatHI (CTCTTC N1/N4).

Structural mechanisms for the 5'-CCWGG sequence recognition by the N- and C-terminal domains of EcoRII

EcoRII restriction endonuclease is specific for the 5′-CCWGG sequence (W stands for A or T); however, it shows no activity on a single recognition site. To activate cleavage it requires binding of an additional target site as an allosteric effector. EcoRII dimer consists of three structural units: a central catalytic core, made from two copies of the C-terminal domain (EcoRII-C), and two N-terminal effector DNA binding domains (EcoRII-N). Here, we report DNA-bound EcoRII-N and EcoRII-C structures, which show that EcoRII combines two radically different structural mechanisms to interact with the effector and substrate DNA. The catalytic EcoRII-C dimer flips out the central T:A base pair and makes symmetric interactions with the CC:GG half-sites. The EcoRII-N effector domain monomer binds to the target site asymmetrically in a single defined orientation which is determined by specific hydrogen bonding and van der Waals interactions with the central T:A pair in the major groove. The EcoRII-N mode of the target site recognition is shared by the large class of higher plant transcription factors of the B3 superfamily.

On the divalent metal ion dependence of DNA cleavage by restriction endonucleases of the EcoRI family

Restriction endonucleases of the PD...D/EXK family need Mg(2+) for DNA cleavage. Whereas Mg(2+) (or Mn(2+)) promotes catalysis, Ca(2+) (without Mg(2+)) only supports DNA binding. The role of Mg(2+) in DNA cleavage by restriction endonucleases has elicited many hypotheses, differing mainly in the number of Mg(2+) involved in catalysis. To address this problem, we measured the Mg(2+) and Mn(2+) concentration dependence of DNA cleavage by BamHI, BglII, Cfr10I, EcoRI, EcoRII (catalytic domain), MboI, NgoMIV, PspGI, and SsoII, which were reported in co-crystal structure analyses to bind one (BglII and EcoRI) or two (BamHI and NgoMIV) Me(2+) per active site. DNA cleavage experiments were carried out at various Mg(2+) and Mn(2+) concentrations at constant ionic strength. All enzymes show a qualitatively similar Mg(2+) and Mn(2+) concentration dependence. In general, the Mg(2+) concentration optimum (between approximately 1 and 10 mM) is higher than the Mn(2+) concentration optimum (between approximately 0.1 and 1 mM). At still higher Mg(2+) or Mn(2+) concentrations, the activities of all enzymes tested are reduced but can be reactivated by Ca(2+). Based on these results, we propose that one Mg(2+) or Mn(2+) is critical for restriction enzyme activation, and binding of a second Me(2+) plays a role in modulating the activity. Steady-state kinetics carried out with EcoRI and BamHI suggest that binding of a second Mg(2+) or Mn(2+) mainly leads to an increase in K(m), such that the inhibitory effect of excess Mg(2+) or Mn(2+) can be overcome by increasing the substrate concentration. Our conclusions are supported by molecular dynamics simulations and are consistent with the structural observations of both one and two Me(2+) binding to these enzymes.

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3' overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding-thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.

Cloning and analysis of a bifunctional methyltransferase/restriction endonuclease TspGWI, the prototype of a Thermus sp. enzyme family

Background

Restriction-modification systems are a diverse class of enzymes. They are classified into four major types: I, II, III and IV. We have previously proposed the existence of a Thermus sp. enzyme family, which belongs to type II restriction endonucleases (REases), however, it features also some characteristics of types I and III. Members include related thermophilic endonucleases: TspGWI, TaqII, TspDTI, and Tth111II.

Results

Here we describe cloning, mutagenesis and analysis of the prototype TspGWI enzyme that recognises the 5'-ACGGA-3' site and cleaves 11/9 nt downstream. We cloned, expressed, and mutagenised the tspgwi gene and investigated the properties of its product, the bifunctional TspGWI restriction/modification enzyme. Since TspGWI does not cleave DNA completely, a cloning method was devised, based on amino acid sequencing of internal proteolytic fragments. The deduced amino acid sequence of the enzyme shares significant sequence similarity with another representative of the Thermus sp. family – TaqII. Interestingly, these enzymes recognise similar, yet different sequences in the DNA. Both enzymes cleave DNA at the same distance, but differ in their ability to cleave single sites and in the requirement of S-adenosylmethionine as an allosteric activator for cleavage. Both the restriction endonuclease (REase) and methyltransferase (MTase) activities of wild type (wt) TspGWI (either recombinant or isolated from Thermus sp.) are dependent on the presence of divalent cations.

Conclusion

TspGWI is a bifunctional protein comprising a tandem arrangement of Type I-like domains; particularly noticeable is the central HsdM-like module comprising a helical domain and a highly conserved S-adenosylmethionine-binding/catalytic MTase domain, containing DPAVGTG and NPPY motifs. TspGWI also possesses an N-terminal PD-(D/E)XK nuclease domain related to the corresponding domains in HsdR subunits, but lacks the ATP-dependent translocase module of the HsdR subunit and the additional domains that are involved in subunit-subunit interactions in Type I systems. The MTase and REase activities of TspGWI are autonomous and can be uncoupled. Structurally and functionally, the TspGWI protomer appears to be a streamlined 'half' of a Type I enzyme.

Contribution of RecFOR machinery of homologous recombination to cell survival after loss of a restriction-modification gene complex

Loss of a type II restriction-modification (RM) gene complex, such as EcoRI, from a bacterial cell leads to death of its descendent cells through attack by residual restriction enzymes on undermethylated target sites of newly synthesized chromosomes. Through such post-segregational host killing, these gene complexes impose their maintenance on their host cells. This finding led to the rediscovery of type II RM systems as selfish mobile elements. The host prokaryote cells were found to cope with such attacks through a variety of means. The RecBCD pathway of homologous recombination in Escherichia coli repairs the lethal lesions on the chromosome, whilst it destroys restricted non-self DNA. recBCD homologues, however, appear very limited in distribution among bacterial genomes, whereas homologues of the RecFOR proteins, responsible for another pathway, are widespread in eubacteria, just like the RM systems. In the present work, therefore, we examined the possible contribution of the RecFOR pathway to cell survival after loss of an RM gene complex. A recF mutation reduced survival in an otherwise rec-positive background and, more severely, in a recBC sbcBC background. We also found that its effect is prominent in the presence of specific non-null mutant forms of the RecBCD enzyme: the resistance to killing seen with recC1002, recC1004, recC2145 and recB2154 is severely reduced to the level of a null recBC allele when combined with a recF, recO or recR mutant allele. Such resistance was also dependent on RecJ and RecQ functions. UV resistance of these non-null recBCD mutants is also reduced by recF, recJ or recQ mutation. These results demonstrate that the RecFOR pathway of recombination can contribute greatly to resistance to RM-mediated host killing, depending on the genetic background.

Crystal structure of the beta beta alpha-Me type II restriction endonuclease Hpy99I with target DNA

The ββα-Me restriction endonuclease (REase) Hpy99I recognizes the CGWCG target sequence and cleaves it with unusual stagger (five nucleotide 5′-recessed ends). Here we present the crystal structure of the specific complex of the dimeric enzyme with DNA. The Hpy99I protomer consists of an antiparallel β-barrel and two β4α2 repeats. Each repeat coordinates a structural zinc ion with four cysteine thiolates in two CXXC motifs. The ββα-Me region of the second β4α2 repeat holds the catalytic metal ion (or its sodium surrogate) via Asp148 and Asn165 and activates a water molecule with the general base His149. In the specific complex, Hpy99I forms a ring-like structure around the DNA that contacts DNA bases on the major and minor groove sides via the first and second β4α2 repeats, respectively. Hpy99I interacts with the central base pair of the recognition sequence only on the minor groove side, where A:T resembles T:A and G:C is similar to C:G. The Hpy99I–DNA co-crystal structure provides the first detailed illustration of the ββα-Me site in REases and complements structural information on the use of this active site motif in other groups of endonucleases such as homing endonucleases (e.g. I-PpoI) and Holliday junction resolvases (e.g. T4 endonuclease VII).

Functional analysis of MmeI from methanol utilizer Methylophilus methylotrophus, a subtype IIC restriction-modification enzyme related to type I enzymes

MmeI from Methylophilus methylotrophus belongs to the type II restriction-modification enzymes. It recognizes an asymmetric DNA sequence, 5'-TCCRAC-3' (R indicates G or A), and cuts both strands at fixed positions downstream of the specific site. This particular feature has been exploited in transcript profiling of complex genomes (using serial analysis of gene expression technology). We have shown previously that the endonucleolytic activity of MmeI is strongly dependent on the presence of S-adenosyl-l-methionine (J. Nakonieczna, J. W. Zmijewski, B. Banecki, and A. J. Podhajska, Mol. Biotechnol. 37:127-135, 2007), which puts MmeI in subtype IIG. The same cofactor is used by MmeI as a methyl group donor for modification of an adenine in the upper strand of the recognition site to N(6)-methyladenine. Both enzymatic activities reside in a single polypeptide (919 amino acids [aa]), which puts MmeI also in subtype IIC of the restriction-modification systems. Based on a molecular model, generated with the use of bioinformatic tools and validated by site-directed mutagenesis, we were able to localize three functional domains in the structure of the MmeI enzyme: (i) the N-terminal portion containing the endonucleolytic domain with the catalytic Mg2+-binding motif D(70)-X(9)-EXK(82), characteristic for the PD-(D/E)XK superfamily of nucleases; (ii) a central portion (aa 310 to 610) containing nine sequence motifs conserved among N(6)-adenine gamma-class DNA methyltransferases; (iii) the C-terminal portion (aa 610 to 919) containing a putative target recognition domain. Interestingly, all three domains showed highest similarity to the corresponding elements of type I enzymes rather than to classical type II enzymes. We have found that MmeI variants deficient in restriction activity (D70A, E80A, and K82A) can bind and methylate specific nucleotide sequence. This suggests that domains of MmeI responsible for DNA restriction and modification can act independently. Moreover, we have shown that a single amino acid residue substitution within the putative target recognition domain (S807A) resulted in a MmeI variant with a higher endonucleolytic activity than the wild-type enzyme.

The Helicobacter pylori HpyAXII restriction-modification system limits exogenous DNA uptake by targeting GTAC sites but shows asymmetric conservation of the DNA methyltransferase and restriction endonuclease components

The naturally competent organism Helicobacter pylori encodes a large number of restriction-modification (R-M) systems that consist of a restriction endonuclease and a DNA methyltransferase. R-M systems are not only believed to limit DNA exchange among bacteria but may also have other cellular functions. We report a previously uncharacterized H. pylori type II R-M system, M.HpyAXII/R.HpyAXII. We show that this system targets GTAC sites, which are rare in the H. pylori chromosome but numerous in ribosomal RNA genes. As predicted, this type II R-M system showed attributes of a selfish element. Deletion of the methyltransferase M.HpyAXII is lethal when associated with an active endonuclease R.HpyAXII unless compensated by adaptive mutation or gene amplification. R.HpyAXII effectively restricted both unmethylated plasmid and chromosomal DNA during natural transformation and was predicted to belong to the novel 'half pipe' structural family of endonucleases. Analysis of a panel of clinical isolates revealed that R.HpyAXII was functional in a small number of H. pylori strains (18.9%, n = 37), whereas the activity of M.HpyAXII was highly conserved (92%, n = 50), suggesting that GTAC methylation confers a selective advantage to H. pylori. However, M.HpyAXII activity did not enhance H. pylori fitness during stomach colonization of a mouse infection model.

Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily, but exhibits an unusual active site

BACKGROUND

Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds.

RESULTS

R.Hpy188I is a Type II REase with unknown structure. PSI-BLAST searches of the non-redundant protein sequence database reveal only 1 homolog (R.HpyF17I, with nearly identical amino acid sequence and the same DNA sequence specificity). Standard application of state-of-the-art protein fold-recognition methods failed to predict the relationship of R.Hpy188I to proteins with known structure or to other protein families. In order to increase the amount of evolutionary information in the multiple sequence alignment, we have expanded our sequence database searches to include sequences from metagenomics projects. This search resulted in identification of 23 further members of R.Hpy188I family, both from metagenomics and the non-redundant database. Moreover, fold-recognition analysis of the extended R.Hpy188I family revealed its relationship to the GIY-YIG domain and allowed for computational modeling of the R.Hpy188I structure. Analysis of the R.Hpy188I model in the light of sequence conservation among its homologs revealed an unusual variant of the active site, in which the typical Tyr residue of the YIG half-motif had been substituted by a Lys residue. Moreover, some of its homologs have the otherwise invariant Arg residue in a non-homologous position in sequence that nonetheless allows for spatial conservation of the guanidino group potentially involved in phosphate binding.

CONCLUSION

The present study eliminates a significant "white spot" on the structural map of REases. It also provides important insight into sequence-structure-function relationships in the GIY-YIG nuclease superfamily. Our results reveal that in the case of proteins with no or few detectable homologs in the standard "non-redundant" database, it is useful to expand this database by adding the metagenomic sequences, which may provide evolutionary linkage to detect more remote homologs.

Activity and specificity of the bacterial PD-(D/E)XK homing endonuclease I-Ssp6803I

The restriction endonuclease fold [a three-layer alpha-beta sandwich containing variations of the PD-(D/E)XK nuclease motif] has been greatly diversified during evolution, facilitating its use for many biological functions. Here we characterize DNA binding and cleavage by the PD-(D/E)XK homing endonuclease I-Ssp6803I. Unlike most restriction endonucleases harboring the same core fold, the specificity profile of this enzyme extends over a long (17 bp) target site. The DNA binding and cleavage specificity profiles of this enzyme were independently determined and found to be highly correlated. However, the DNA target sequence contains several positions where binding and cleavage activities are not tightly coupled: individual DNA base-pair substitutions at those positions that significantly decrease cleavage activity have minor effects on binding affinity. These changes in the DNA target sequence appear to correspond to substitutions that uniquely increase the free energy change between the ground state and the transition state, rather than simply decreasing the overall DNA binding affinity. The specificity of the enzyme reflects constraints on its host gene and limitations imposed by the enzyme's quaternary structure and illustrate the highly diverse repertoire of DNA recognition specificities that can be adopted by the related folds surrounding the PD-(D/E)XK nuclease motif.