Engineering of restriction endonucleases: using methylation activity of the bifunctional endonuclease Eco57I to select the mutant with a novel sequence specificity

Type II restriction endonucleases (REs) are widely used tools in molecular biology, biotechnology and diagnostics. Efforts to generate new specificities by structure-guided design and random mutagenesis have been unsuccessful so far. We have developed a new procedure called the methylation activity-based selection (MABS) for generating REs with a new specificity. MABS uses a unique property of bifunctional type II REs to methylate DNA targets they recognize. The procedure includes three steps: (1) conversion of a bifunctional RE into a monofunctional DNA-modifying enzyme by cleavage center disruption; (2) mutagenesis and selection of mutants with altered DNA modification specificity based on their ability to protect predetermined DNA targets; (3) reconstitution of the cleavage center's wild-type structure. The efficiency of the MABS technique was demonstrated by altering the sequence specificity of the bifunctional RE Eco57I from 5'-CTGAAG to 5'-CTGRAG, and thus generating the mutant restriction endonuclease (and DNA methyltransferase) of a specificity not known before. This study provides evidence that MABS is a promising technique for generation of REs with new specificities.     

The generation of new protein functions by the combination of domains

During evolution, many new proteins have been formed by the process of gene duplication and combination. The genes involved in this process usually code for whole domains. Small proteins contain one domain; medium and large proteins contain two or more domains. We have compared homologous domains that occur in both one-domain proteins and multidomain proteins. We have determined (1) how the functions of the individual domains in the multidomain proteins combine to produce their overall functions and (2) the extent to which these functions are similar to those in the one-domain homologs. We describe how domain combinations increase the specificity of enzymes; act as links between domains that have functional roles; regulate activity; combine within one chain functions that can act either independently, in concert or in new contexts; and provide the structural framework for the evolution of entirely new functions.     

Grouping together highly diverged PD-(D/E)XK nucleases and identification of novel superfamily members using structure-guided alignment of sequence profiles

The PD-(D/E)XK nuclease domains, initially identified in type II restriction enzymes, serve as models for studying aspects of protein-DNA interactions, mechanisms of phosphodiester hydrolysis, and provide indispensable tools for techniques in genetic engineering and molecular medicine. However, the low degree of amino acid conservation hampers the possibility of identification of PD-(D/E)XK superfamily members based solely on sequence comparisons. In several proteins implicated in DNA recombination and repair the restriction enzyme-like nuclease domain has been found only after the corresponding structures were determined experimentally. Here, we identified highly diverged variants of the PD-(D/E)XK domain in many proteins and open reading frames using iterative database searches and progressive, structure-guided alignment of sequence profiles. We predicted the possible cellular function for many hypothetical proteins based on their relative similarity to characterized nucleases or observed presence of additional domains. We also identified the nuclease domain in genuine recombinases and restriction enzymes, whose homology to other PD-(D/E)XK enzymes has not been demonstrated previously. The first superfamily-wide comparative analysis, not limited to nucleases of known structure, will guide cloning and characterization of novel enzymes and planning new experiments to better understand those already studied.     

SURVEY AND SUMMARY: holliday junction resolvases and related nucleases: identification of new families, phyletic distribution and evolutionary trajectories

Holliday junction resolvases (HJRs) are key enzymes of DNA recombination. A detailed computer analysis of the structural and evolutionary relationships of HJRs and related nucleases suggests that the HJR function has evolved independently from at least four distinct structural folds, namely RNase H, endonuclease, endonuclease VII-colicin E and RusA. The endonuclease fold, whose structural prototypes are the phage lambda exonuclease, the very short patch repair nuclease (Vsr) and type II restriction enzymes, is shown to encompass by far a greater diversity of nucleases than previously suspected. This fold unifies archaeal HJRs, repair nucleases such as RecB and Vsr, restriction enzymes and a variety of predicted nucleases whose specific activities remain to be determined. Within the RNase H fold a new family of predicted HJRs, which is nearly ubiquitous in bacteria, was discovered, in addition to the previously characterized RuvC family. The proteins of this family, typified by Escherichia coli YqgF, are likely to function as an alternative to RuvC in most bacteria, but could be the principal HJRs in low-GC Gram-positive bacteria and AQUIFEX: Endonuclease VII of phage T4 is shown to serve as a structural template for many nucleases, including MCR:A and other type II restriction enzymes. Together with colicin E7, endonuclease VII defines a distinct metal-dependent nuclease fold. As a result of this analysis, the principal HJRs are now known or confidently predicted for all bacteria and archaea whose genomes have been completely sequenced, with many species encoding multiple potential HJRs. Horizontal gene transfer, lineage-specific gene loss and gene family expansion, and non-orthologous gene displacement seem to have been major forces in the evolution of HJRs and related nucleases. A remarkable case of displacement is seen in the Lyme disease spirochete Borrelia burgdorferi, which does not possess any of the typical HJRs, but instead encodes, in its chromosome and each of the linear plasmids, members of the lambda exonuclease family predicted to function as HJRs. The diversity of HJRs and related nucleases in bacteria and archaea contrasts with their near absence in eukaryotes. The few detected eukaryotic representatives of the endonuclease fold and the RNase H fold have probably been acquired from bacteria via horizontal gene transfer. The identity of the principal HJR(s) involved in recombination in eukaryotes remains uncertain; this function could be performed by topoisomerase IB or by a novel, so far undetected, class of enzymes. Likely HJRs and related nucleases were identified in the genomes of numerous bacterial and eukaryotic DNA viruses. Gene flow between viral and cellular genomes has probably played a major role in the evolution of this class of enzymes. This analysis resulted in the prediction of numerous previously unnoticed nucleases, some of which are likely to be new restriction enzymes.     

Understanding the evolution of restriction-modification systems: clues from sequence and structure comparisons

Restriction-modification (RM) systems comprise two opposing enzymatic activities: a restriction endonuclease, that targets specific DNA sequences and performs endonucleolytic cleavage, and a modification methyltransferase that renders these sequences resistant to cleavage. Studies on molecular genetics and biochemistry of RM systems have been carried out over the past four decades, laying foundations for modern molecular biology and providing important models for mechanisms of highly specific protein-DNA interactions. Although the number of known, relevant sequences 3D structures of RM proteins is growing steadily, we do not fully understand their functional diversities from an evolutionary perspective and we are not yet able to engineer new sequence specificities based on rational approaches. Recent findings on the evolution of RM systems and on their structures and mechanisms of action have led to a picture in which conserved modules with defined function are shared between different RM proteins and other enzymes involved in nucleic acid biochemistry. On the other hand, it has been realized that some of the modules have been replaced in the evolution by unrelated domains exerting similar function. The aim of this review is to give a survey on the recent progress in the field of structural phylogeny of RM enzymes with special emphasis on studies of sequence-structure-function relationships and emerging potential applications in biotechnology.     

Treble clef finger--a functionally diverse zinc-binding structural motif

Detection of similarity is particularly difficult for small proteins and thus connections between many of them remain unnoticed. Structure and sequence analysis of several metal-binding proteins reveals unexpected similarities in structural domains classified as different protein folds in SCOP and suggests unification of seven folds that belong to two protein classes. The common motif, termed treble clef finger in this study, forms the protein structural core and is 25-45 residues long. The treble clef motif is assembled around the central zinc ion and consists of a zinc knuckle, loop, beta-hairpin and an alpha-helix. The knuckle and the first turn of the helix each incorporate two zinc ligands. Treble clef domains constitute the core of many structures such as ribosomal proteins L24E and S14, RING fingers, protein kinase cysteine-rich domains, nuclear receptor-like fingers, LIM domains, phosphatidylinositol-3-phosphate-binding domains and His-Me finger endonucleases. The treble clef finger is a uniquely versatile motif adaptable for various functions. This small domain with a 25 residue structural core can accommodate eight different metal-binding sites and can have many types of functions from binding of nucleic acids, proteins and small molecules, to catalysis of phosphodiester bond hydrolysis. Treble clef motifs are frequently incorporated in larger structures or occur in doublets. Present analysis suggests that the treble clef motif defines a distinct structural fold found in proteins with diverse functional properties and forms one of the major zinc finger groups.     

Toprim--a conserved catalytic domain in type IA and II topoisomerases, DnaG-type primases, OLD family nucleases and RecR proteins.

Iterative profile searches and structural modeling show that bacterial DnaG-type primases, small primase-like proteins from bacteria and archaea, type IA and type II topoisomerases, bacterial and archaeal nucleases of the OLD family and bacterial DNA repair proteins of the RecR/M family contain a common domain, designated Toprim (topoisomerase-primase) domain. The domain consists of approximately 100 amino acids and has two conserved motifs, one of which centers at a conserved glutamate and the other one at two conserved aspartates (DxD). Examination of the structure of Topo IA and Topo II and modeling of the Toprim domains of the primases reveal a compact beta/alpha fold, with the conserved negatively charged residues juxtaposed, and inserts seen in Topo IA and Topo II. The conserved glutamate may act as a general base in nucleotide polymerization by primases and in strand rejoining by topoisomerases and as a general acid in strand cleavage by topoisomerases and nucleases. The role of this glutamate in catalysis is supported by site-directed mutagenesis data on primases and Topo IA. The DxD motif may coordinate Mg2+that is required for the activity of all Toprim-containing enzymes. The common ancestor of all life forms could encode a prototype Toprim enzyme that might have had both nucleotidyl transferase and polynucleotide cleaving activity.     

Evolution of protein function, from a structural perspective.

The recent growth in structural data, and ensuing analyses, have revealed the structural and functional versatility of protein families. With respect to enzymes, local active-site mutations, variations in surface loops and recruitment of additional domains accommodate the diverse substrate specificities and catalytic activities observed within several superfamilies. Conversely, some functions have more than one structural solution, having evolved independently several times during evolution. Combined with the existence of multi-functional genes, which have arisen by gene recruitment, these phenomena must be considered in the process of genome annotation.     

Theoretical model of restriction endonuclease HpaI in complex with DNA, predicted by fold recognition and validated by site-directed mutagenesis

Type II restriction enzymes are commercially important deoxyribonucleases and very attractive targets for protein engineering of new specificities. At the same time they are a very challenging test bed for protein structure prediction methods. Typically, enzymes that recognize different sequences show little or no amino acid sequence similarity to each other and to other proteins. Based on crystallographic analyses that revealed the same PD-(D/E)XK fold for more than a dozen case studies, they were nevertheless considered to be related until the combination of bioinformatics and mutational analyses has demonstrated that some of these proteins belong to other, unrelated folds PLD, HNH, and GIY-YIG. As a part of a large-scale project aiming at identification of a three-dimensional fold for all type II REases with known sequences (currently approximately 1000 proteins), we carried out preliminary structure prediction and selected candidates for experimental validation. Here, we present the analysis of HpaI REase, an ORFan with no detectable homologs, for which we detected a structural template by protein fold recognition, constructed a model using the FRankenstein monster approach and identified a number of residues important for the DNA binding and catalysis. These predictions were confirmed by site-directed mutagenesis and in vitro analysis of the mutant proteins. The experimentally validated model of HpaI will serve as a low-resolution structural platform for evolutionary considerations in the subgroup of blunt-cutting REases with different specificities. The research protocol developed in the course of this work represents a streamlined version of the previously used techniques and can be used in a high-throughput fashion to build and validate models for other enzymes, especially ORFans that exhibit no sequence similarity to any other protein in the database.     

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

Linear adenovirus DNA is organized into supercoiled domains in virus particles.

Electron microscopic analysis of bis-psoralen crosslinked adenovirus type 5 virion DNA revealed supercoiled domains in an otherwise linear DNA. The existence of supercoiled arrangement in all the virion DNA was demonstrated by the sensitivity of Ad5 DNA in pentonless virus particles to the supercoiling-dependent endonucleolytic activity of Bal31 and S1 nucleases. These nucleases were found to cleave Ad5 virion DNA at specific sites. The observation of stable cleavage sites in the limit digestion of virion DNA by Bal31 suggests that cleavage sites represent boundaries of core proteins which impede the exonuclease activity of Bal31. These data suggest that specific arrangement of core proteins on Ad5 virion DNA. Based on this analysis we determined positions of core proteins in viral genome using indirect end labeling technique. The size of supercoiled domains of virion DNA was estimated by electron microscopy and also by boundaries of mutually exclusive Bal31 cleavage sites at limit digestion condition. Our data suggest each supercoiled domain is equal to about 12% of Ad5 genome length and about 8 loops can be accommodated in Ad5 virion. However sequences at two extreme ends of the viral genome were found to be outside of supercoiled domains. An interesting correlation between supercoiled domains and gene domains of Ad5 genome was noticed.     

Structure model of core proteins in photosystem I inferred from the comparison with those in photosystem II and bacteria; an application of principal component analysis to detect the similar regions between distantly related families of proteins.

A principal component analysis based on the physico-chemical properties of amino acid residues is developed to assign similar regions between distantly related families of proteins, taking account of the species diversities in respective families. The most important advantage of this analysis should be that it reflects different physico-chemical properties and thus can predict more detailed structural properties, including the transmembrane helices, than the hydropathy analysis. Its first application reconfirms the similarity between the core proteins of photosynthetic reaction center in purple bacteria and those of photosystem II, indicating that the low percentage of identical amino acid residues estimated previously between them is due to much allowance for amino acid substitutions in purple bacteria. The application of this analysis to the core proteins of photosystem I reveals that any of these proteins includes two domains, each showing high similarity to the amino acid sequences of core proteins in photosystem II and purple bacteria. A core structure model of A1 and A2 proteins folded into four layers of sheets of transmembrane helices is proposed to provide a molecular basis for the electron pathway suggested by spectroscopic experiments as well as for the interaction sites with plastocyanin, 9 kDa protein and LHC proteins.     

Resolving protein structure‐function‐binding site relationships from a binding site similarity network perspective

Functional annotation is seldom straightforward with complexities arising due to functional divergence in protein families or functional convergence between non‐homologous protein families, leading to mis‐annotations. An enzyme may contain multiple domains and not all domains may be involved in a given function, adding to the complexity in function annotation. To address this, we use binding site information from bound cognate ligands and catalytic residues, since it can help in resolving fold‐function relationships at a finer level and with higher confidence. A comprehensive database of 2,020 fold‐function‐binding site relationships has been systematically generated. A network‐based approach is employed to capture the complexity in these relationships, from which different types of associations are deciphered, that identify versatile protein folds performing diverse functions, same function associated with multiple folds and one‐to‐one relationships. Binding site similarity networks integrated with fold, function, and ligand similarity information are generated to understand the depth of these relationships. Apart from the observed continuity in the functional site space, network properties of these revealed versatile families with topologically different or dissimilar binding sites and structural families that perform very similar functions. As a case study, subtle changes in the active site of a set of evolutionarily related superfamilies are studied using these networks. Tracing of such similarities in evolutionarily related proteins provide clues into the transition and evolution of protein functions. Insights from this study will be helpful in accurate and reliable functional annotations of uncharacterized proteins, poly‐pharmacology, and designing enzymes with new functional capabilities. Proteins 2017; 85:1319–1335. © 2017 Wiley Periodicals, Inc.     

Thematic minireview series on enzyme evolution in the post-genomic era

The emergence of genomics; ongoing computational advances; and the development of large-scale sequence, structural, and functional databases have created important new interdisciplinary linkages between molecular evolution, molecular biology, and enzymology. The five minireviews in this series survey advances and challenges in this burgeoning field from complementary perspectives. The series has three major themes. The first is the evolution of enzyme superfamilies, in which members exhibit increasing sequence, structural, and functional divergence with increasing time of divergence from a common ancestor. The second is the evolutionary role of promiscuous enzymes, which, in addition to their primary function, have adventitious secondary activities that frequently provide the starting point for the evolution of new enzymes. The third is the importance of in silico approaches to the daunting challenge of assigning and predicting the functions of the many uncharacterized proteins in the large-scale sequence and structural databases that are now available. A recent computational advance, the use of protein similarity networks that map functional data onto proteins clustered by similarity, is presented as an approach that can improve functional insight and inference. The three themes are illustrated with several examples of enzyme superfamilies, including the amidohydrolase, metallo-β-lactamase, and enolase superfamilies.     

Stimulation of homologous recombination through targeted cleavage by chimeric nucleases

Chimeric nucleases that are hybrids between a nonspecific DNA cleavage domain and a zinc finger DNA recognition domain were tested for their ability to find and cleave their target sites in living cells. Both engineered DNA substrates and the nucleases were injected into Xenopus laevis oocyte nuclei, in which DNA cleavage and subsequent homologous recombination were observed. Specific cleavage required two inverted copies of the zinc finger recognition site in close proximity, reflecting the need for dimerization of the cleavage domain. Cleaved DNA molecules were activated for homologous recombination; in optimum conditions, essentially 100% of the substrate recombined, even though the DNA was assembled into chromatin. The original nuclease has an 18-amino-acid linker between the zinc finger and cleavage domains, and this enzyme cleaved in oocytes at paired sites separated by spacers in the range of 6 to 18 bp, with a rather sharp optimum at 8 bp. By shortening the linker, we found that the range of effective site separations could be narrowed significantly. With no intentional linker between the binding and cleavage domains, only binding sites exactly 6 bp apart supported efficient cleavage in oocytes. We also showed that two chimeric enzymes with different binding specificities could collaborate to stimulate recombination when their individual sites were appropriately placed. Because the recognition specificity of zinc fingers can be altered experimentally, this approach holds great promise for inducing targeted recombination in a variety of organisms.     

The prokaryotic enzyme DsbB may share key structural features with eukaryotic disulfide bond forming oxidoreductases

Three different classes of thiol-oxidoreductases that facilitate the formation of protein disulfide bonds have been identified. They are the Ero1 and SOX/ALR family members in eukaryotic cells, and the DsbB family members in prokaryotic cells. These enzymes transfer oxidizing potential to the proteins PDI or DsbA, which are responsible for directly introducing disulfide bonds into substrate proteins during oxidative protein folding in eukaryotes and prokaryotes, respectively. A comparison of the recent X-ray crystal structure of Ero1 with the previously solved structure of the SOX/ALR family member Erv2 reveals that, despite a lack of primary sequence homology between Ero1 and Erv2, the core catalytic domains of these two proteins share a remarkable structural similarity. Our search of the DsbB protein sequence for features found in the Ero1 and Erv2 structures leads us to propose that, in a fascinating example of structural convergence, the catalytic core of this integral membrane protein may resemble the soluble catalytic domain of Ero1 and Erv2. Our analysis of DsbB also identified two new groups of DsbB proteins that, based on sequence homology, may also possess a catalytic core similar in structure to the catalytic domains of Ero1 and Erv2.