Comprehensive homing endonuclease target site specificity profiling reveals evolutionary constraints and enables genome engineering applications

Homing endonucleases (HEs) promote the evolutionary persistence of selfish DNA elements by catalyzing element lateral transfer into new host organisms. The high site specificity of this lateral transfer reaction, termed homing, reflects both the length (14–40 bp) and the limited tolerance of target or homing sites for base pair changes. In order to better understand molecular determinants of homing, we systematically determined the binding and cleavage properties of all single base pair variant target sites of the canonical LAGLIDADG homing endonucleases I-CreI and I-MsoI. These Chlorophyta algal HEs have very similar three-dimensional folds and recognize nearly identical 22 bp target sites, but use substantially different sets of DNA-protein contacts to mediate site-specific recognition and cleavage. The site specificity differences between I-CreI and I-MsoI suggest different evolutionary strategies for HE persistence. These differences also provide practical guidance in target site finding, and in the generation of HE variants with high site specificity and cleavage activity, to enable genome engineering applications.     

Reprogramming homing endonuclease specificity through computational design and directed evolution

Homing endonucleases (HEs) can be used to induce targeted genome modification to reduce the fitness of pathogen vectors such as the malaria-transmitting Anopheles gambiae and to correct deleterious mutations in genetic diseases. We describe the creation of an extensive set of HE variants with novel DNA cleavage specificities using an integrated experimental and computational approach. Using computational modeling and an improved selection strategy, which optimizes specificity in addition to activity, we engineered an endonuclease to cleave in a gene associated with Anopheles sterility and another to cleave near a mutation that causes pyruvate kinase deficiency. In the course of this work we observed unanticipated context-dependence between bases which will need to be mechanistically understood for reprogramming of specificity to succeed more generally.     

Context dependence between subdomains in the DNA binding interface of the I-CreI homing endonuclease

Homing endonucleases (HE) have emerged as precise tools for achieving gene targeting events. Redesigned HEs with tailored specificities can be used to cleave new sequences, thereby considerably expanding the number of targetable genes and loci. With HEs, as well as with other protein scaffolds, context dependence of DNA/protein interaction patterns remains one of the major limitations for rational engineering of new DNA binders. Previous studies have shown strong crosstalk between different residues and regions of the DNA binding interface. To investigate this phenomenon, we systematically combined mutations from three groups of amino acids in the DNA binding regions of the I-CreI HE. Our results confirm that important crosstalk occurs throughout this interface in I-CreI. Detailed analysis of success rates identified a nearest-neighbour effect, with a more pronounced level of dependence between adjacent regions. Taken together, these data suggest that combinatorial engineering does not necessarily require the identification of separable functional or structural regions, and that groups of amino acids provide acceptable building blocks that can be assembled, overcoming the context dependency of the DNA binding interface. Furthermore, the present work describes a sequential method to engineer tailored HEs, wherein three contiguous regions are individually mutated and assembled to create HEs with engineered specificity.     

Homing endonuclease I-CreI derivatives with novel DNA target specificities

Homing endonucleases are highly specific enzymes, capable of recognizing and cleaving unique DNA sequences in complex genomes. Since such DNA cleavage events can result in targeted allele-inactivation and/or allele-replacement in vivo, the ability to engineer homing endonucleases matched to specific DNA sequences of interest would enable powerful and precise genome manipulations. We have taken a step-wise genetic approach in analyzing individual homing endonuclease I-CreI protein/DNA contacts, and describe here novel interactions at four distinct target site positions. Crystal structures of two mutant endonucleases reveal the molecular interactions responsible for their altered DNA target specificities. We also combine novel contacts to create an endonuclease with the predicted target specificity. These studies provide important insights into engineering homing endonucleases with novel target specificities, as well as into the evolution of DNA recognition by this fascinating family of proteins.     

Computer design of obligate heterodimer meganucleases allows efficient cutting of custom DNA sequences

Meganucleases cut long (>12 bp) unique sequences in genomes and can be used to induce targeted genome engineering by homologous recombination in the vicinity of their cleavage site. However, the use of natural meganucleases is limited by the repertoire of their target sequences, and considerable efforts have been made to engineer redesigned meganucleases cleaving chosen targets. Homodimeric meganucleases such as I-CreI have provided a scaffold, but can only be modified to recognize new quasi-palindromic DNA sequences, limiting their general applicability. Other groups have used dimer-interface redesign and peptide linkage to control heterodimerization between related meganucleases such as I-DmoI and I-CreI, but until now there has been no application of this aimed specifically at the scaffolds from existing combinatorial libraries of I-CreI. Here, we show that engineering meganucleases to form obligate heterodimers results in functional endonucleases that cut non-palindromic sequences. The protein design algorithm (FoldX v2.7) was used to design specific heterodimer interfaces between two meganuclease monomers, which were themselves engineered to recognize different DNA sequences. The new monomers favour functional heterodimer formation and prevent homodimer site recognition. This design massively increases the potential repertoire of DNA sequences that can be specifically targeted by designed I-CreI meganucleases and opens the way to safer targeted genome engineering.     

Generation of single-chain LAGLIDADG homing endonucleases from native homodimeric precursor proteins

Homing endonucleases (HEs) cut long DNA target sites with high specificity to initiate and target the lateral transfer of mobile introns or inteins. This high site specificity of HEs makes them attractive reagents for gene targeting to promote DNA modification or repair. We have generated several hundred catalytically active, monomerized versions of the well-characterized homodimeric I-CreI and I-MsoI LAGLIDADG family homing endonuclease (LHE) proteins. Representative monomerized I-CreI and I-MsoI proteins (collectively termed mCreIs or mMsoIs) were characterized in detail by using a combination of biochemical, biophysical and structural approaches. We also demonstrated that both mCreI and mMsoI proteins can promote cleavage-dependent recombination in human cells. The use of single chain LHEs should simplify gene modification and targeting by requiring the expression of a single small protein in cells, rather than the coordinate expression of two separate protein coding genes as is required when using engineered heterodimeric zinc finger or homing endonuclease proteins.     

Alteration of the C-Terminal Ligand Specificity of the Erbin PDZ Domain by Allosteric Mutational Effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.     

Rational engineering of type II restriction endonuclease DNA binding and cleavage specificity

The type II restriction endonucleases are indispensible tools for molecular biology. Although enzymes recognizing nearly 300 unique sequences are known, the ability to engineer enzymes to recognize any sequence of choice would be valuable. However, previous attempts to engineer new recognition specificity have met limited success. Here we report the rational engineering of multiple new type II specificities. We recently identified a family of MmeI-like type II endonucleases that have highly similar protein sequences but different recognition specificity. We identified the amino-acid positions within these enzymes that determine position specific DNA base recognition at three positions within their recognition sequences through correlations between their aligned amino-acid residues and aligned recognition sequences. We then altered the amino acids at the identified positions to those correlated with recognition of a desired new base to create enzymes that recognize and cut at predictable new DNA sequences. The enzymes so altered have similar levels of endonuclease activity compared to the wild-type enzymes. Using simple and predictable mutagenesis in this family it is now possible to create hundreds of unique new type II restriction endonuclease specificities. The findings suggest a simple mechanism for the evolution of new DNA specificity in Nature.     

Recombination of hybrid target sites by binary combinations of Flp variants: mutations that foster interprotomer collaboration and enlarge substrate tolerance

Strategies of directed evolution and combinatorial mutagenesis applied to the Flp site-specific recombinase have yielded recombination systems that utilize bi-specific hybrid target sites. A hybrid site is assembled from two half-sites, each harboring a distinct binding specificity. Satisfying the two specificities by a binary combination of Flp variants, while necessary, may not be sufficient to elicit recombination. We have identified amino acid substitutions that foster interprotomer collaboration between partner Flp variants to potentiate strand exchange in hybrid sites. One such substitution, A35T, acts specifically in cis with one of the two partners of a variant pair, Flp(K82M) and Flp(A35T, R281V). The same A35T mutation is also present within a group of mutations that rescue a Flp variant, Flp(Y60S), that is defective in establishing monomer-monomer interactions on the native Flp target site. Strikingly, these mutations are localized to peptide regions involved in interdomain and interprotomer interactions within the recombination complex. The same group of mutations, when transferred to the context of wild-type Flp, can relax its specificity to include non-native target sites. The hybrid Flp systems described here mimic the naturally occurring XerC/XerD recombination system that utilizes two recombinases with distinct DNA binding specificities. The ability to overcome the constraints of binding site symmetry in Flp recombination has important implications in the targeted manipulations of genomes.     

Evolutionary transitions to new DNA methyltransferases through target site expansion and shrinkage

DNA-binding and modifying proteins show high specificity but also exhibit a certain level of promiscuity. Such latent promiscuous activities comprise the starting points for new protein functions, but this hypothesis presents a paradox: a new activity can only evolve if it already exists. How then, do novel activities evolve? DNA methyltransferases, for example, are highly divergent in their target sites, but how transitions toward novel sites occur remains unknown. We performed laboratory evolution of the DNA methyltransferase M.HaeIII. We found that new target sites emerged primarily through expansion of the original site, GGCC, and the subsequent shrinkage of evolved expanded sites. Variants evolved for sites that are promiscuously methylated by M.HaeIII [GG(^A/_T)CC and GGCGCC] carried mutations in ‘gate-keeper’ residues. They could thereby methylate novel target sites such as GCGC and GGATCC that were neither selected for nor present in M.HaeIII. These ‘generalist’ intermediates were further evolved to obtain variants with novel target specificities. Our results demonstrate the ease by which new DNA-binding and modifying specificities evolve and the mechanism by which they occur at both the protein and DNA levels.     

Computational design of protein antigens that interact with the CDR H3 loop of HIV broadly neutralizing antibody 2F5

Rational design of proteins with novel binding specificities and increased affinity is one of the major goals of computational protein design. Epitope‐scaffolds are a new class of antigens engineered by transplanting viral epitopes of predefined structure to protein scaffolds, or by building protein scaffolds around such epitopes. Epitope‐scaffolds are of interest as vaccine components to attempt to elicit neutralizing antibodies targeting the specified epitope. In this study we developed a new computational protocol, MultiGraft Interface, that transplants epitopes but also designs additional scaffold features outside the epitope to enhance antibody‐binding specificity and potentially influence the specificity of elicited antibodies. We employed MultiGraft Interface to engineer novel epitope‐scaffolds that display the known epitope of human immunodeficiency virus 1 (HIV‐1) neutralizing antibody 2F5 and that also interact with the functionally important CDR H3 antibody loop. MultiGraft Interface generated an epitope‐scaffold that bound 2F5 with subnanomolar affinity (K_D = 400 pM) and that interacted with the antibody CDR H3 loop through computationally designed contacts. Substantial structural modifications were necessary to engineer this antigen, with the 2F5 epitope replacing a helix in the native scaffold and with 15% of the native scaffold sequence being modified in the design stage. This epitope‐scaffold represents a successful example of rational protein backbone engineering and protein–protein interface design and could prove useful in the field of HIV vaccine design. MultiGraft Interface can be generally applied to engineer novel binding partners with altered specificity and optimized affinity. Proteins 2014; 82:2770–2782. © 2014 Wiley Periodicals, Inc.     

Target DNA bending is an important specificity determinant in target site selection in Tn10 transposition

The bacterial transposon Tn10 inserts preferentially into specific DNA sequences. DNA footprinting and interference studies have revealed that the Tn10-encoded transposase protein contacts a large stretch of target DNA ( approximately 24 bp) and that the target DNA structure is deformed upon incorporation into the transpososome. Target DNA deformation might contribute significantly to target site selection and thus it is of interest to further define the nature of this deformation. Circular permutation analysis was used to demonstrate that the target DNA is bent upon its incorporation into the transpososome. Two lines of evidence are presented that target DNA bending is an important event in target site selection. First, we demonstrate a correlation between increased target site usage and an increased level of target DNA bending. Second, transposase mutants with relaxed target specificity are shown to cause increased target DNA bending relative to wild-type transposase. This latter observation provides new insight into how relaxed specificity may be achieved. We also show that Ca(2+) facilitates target capture by stabilizing transposase interactions with sequences immediately flanking the insertion site. Ca(2+) could, in theory, exert this effect by stabilizing bends in the target DNA.     

I-Ssp6803I: the first homing endonuclease from the PD-(D/E)XK superfamily exhibits an unusual mode of DNA recognition

MOTIVATION: Restriction endonucleases (REases) and homing endonucleases (HEases) are biotechnologically important enzymes. Nearly all structurally characterized REases belong to the PD-(D/E)XK superfamily of nucleases, while most HEases belong to an unrelated LAGLIDADG superfamily. These two protein folds are typically associated with very different modes of protein-DNA recognition, consistent with the different mechanisms of action required to achieve high specificity. REases recognize short DNA sequences using multiple contacts per base pair, while HEases recognize very long sites using a few contacts per base pair, thereby allowing for partial degeneracy of the target sequence. Thus far, neither REases with the LAGLIDADG fold, nor HEases with the PD-(D/E)XK fold, have been found. RESULTS: Using protein fold recognition, we have identified the first member of the PD-(D/E)XK superfamily among homing endonucleases, a cyanobacterial enzyme I-Ssp6803I. We present a model of the I-Ssp6803I-DNA complex based on the structure of Type II restriction endonuclease R.BglI and predict the active site and residues involved in specific DNA sequence recognition by I-Ssp6803I. Our finding reveals a new unexpected evolutionary link between HEases and REases and suggests how PD-(D/E)XK nucleases may develop a 'HEase-like' way of interacting with the extended DNA sequence. This in turn may be exploited to study the evolution of DNA sequence specificity and to engineer nucleases with new substrate specificities.     

The two homeodomains of the ZmHox2a gene from maize originated as an internal gene duplication and have evolved different target site specificities.

The maize ZmHox2a gene encodes two homeodomains which originated by a 699 bp duplication within an ancestral precursor. The sequences of the two ZmHox2a homeodomains are highly diverged in the N-terminal arm, while residues in the helical part have mostly been conserved. We show here that both ZmHox2a homeodomains are functional DNA-binding motifs but exhibit different target site specificities. CASTing experiments reveal a TCCT motif recognized by HD1 but a GATC tetranucleotide as the recognition sequence of HD2. Mutation of the central nucleotides in both tetranucleotide core motifs abolishes DNA binding. A domain swap experiment indicates that target site specificity is achieved in a combinatorial manner by the contributions of the diverged N-terminal arms together with the slightly different recognition helices. Computer modelling suggests that K47 and H54 in the recognition helices preferentially contact the bases at the 3'-terminus of the tetranucleotide target sequences.     

Combinatorial redesign of the DNA binding specificity of a prokaryotic helix-turn-helix repressor

Redesign of the bacteriophage 434 Cro repressor was accomplished by using an in vivo genetic screening system to identify new variants that specifically bound previously unrecognized DNA sequences. Site-directed, combinatorial mutagenesis of the 434 Cro helix-turn-helix (HTH) motif generated libraries of new variants which were screened for binding to new target sequences. Multiple mutations of 434 Cro that functionally converted wild-type (wt) 434 Cro DNA binding-sequence specificity to that of a lambda bacteriophage-specific repressor were identified. The libraries contained variations within the HTH sequence at only three positions. In vivo and in vitro analysis of several of the identified 434 Cro variants showed that the relatively few changes in the recognition helix of the HTH motif of 434 Cro resulted in specific and tight binding of the target DNA sequences. For the best 434 Cro variant identified, an apparent K(d) for lambda O(R)3 of 1 nM was observed. In competition experiments, this Cro variant was observed to be highly selective. We conclude that functional 434 Cro repressor variants with new DNA binding specificities can be generated from wt 434 Cro by mutating just the recognition helix. Important characteristics of the screening system responsible for the successful identifications are discussed. Application of the techniques presented here may allow the identification of DNA binding protein variants that functionally affect DNA regulatory sequences important in disease and industrial and biotechnological processes.     

Mining endonuclease cleavage determinants in genomic sequence data

Homing endonucleases have great potential as tools for targeted gene therapy and gene correction, but identifying variants of these enzymes capable of cleaving specific DNA targets of interest is necessary before the widespread use of such technologies is possible. We identified homologues of the LAGLIDADG homing endonuclease I-AniI and their putative target insertion sites by BLAST searches followed by examination of the sequences of the flanking genomic regions. Amino acid substitutions in these homologues that were located close to the target site DNA, and thus potentially conferring differences in target specificity, were grafted onto the I-AniI scaffold. Many of these grafts exhibited novel and unexpected specificities. These findings show that the information present in genomic data can be exploited for endonuclease specificity redesign.