In Vitro Directed Evolution of a Restriction Endonuclease with More Stringent Specificity

Restriction endonuclease (REase) specificity engineering is extremely difficult. Here we describe a multistep protocol that helps to produce REase ...

International Institute of Molecular and Cell Biology, Warsaw (IIMCB)

Specific peptidases exist for nearly every amide linkage in peptidoglycan. In several cases, families of peptidoglycan hydrolases with different specificities turned out to be related. Here we show that lysostaphin-type peptidases and D-Ala-D-Ala metallopeptidases have similar active sites and share a core folding motif in otherwise highly divergent folds. The central Zn(2+) is tetrahedrally coordinated by two histidines, an aspartate, and a water molecule. The Zn(2+) chelating residues occur in the order histidine, aspartate, histidine in all sequences and contact the metal via the Nepsilon, the Odelta, and the Ndelta, respectively. The identity of the other active-site residues varies, but in all enzymes of known structure except for VanX, a conserved histidine is present two residues upstream of the second histidine ligand to the Zn(2+). As the same arrangement of active-site residues is also found in the N-terminal, cryptic peptidase domain of sonic hedgehog, we propose that this arrangement of active-site residues be called the "LAS" arrangement, because it is present in lysostaphin-type enzymes, D-Ala-D-Ala metallopeptidases, and in the cryptic peptidase in the N-domain of sonic hedgehog.

Similar active sites in lysostaphins and D-Ala-D-Ala metallopeptidases.

Restricion endonuclease Ecl18kI is specific for the sequence /CCNGG and cleaves it before the outer C to generate 5 nt 5'-overhangs. It has been suggested that Ecl18kI is evolutionarily related to NgoMIV, a 6-bp cutter that cleaves the sequence G/CCGGC and leaves 4 nt 5'-overhangs. Here, we report the crystal structure of the Ecl18kI-DNA complex at 1.7 A resolution and compare it with the known structure of the NgoMIV-DNA complex. We find that Ecl18kI flips both central nucleotides within the CCNGG sequence and buries the extruded bases in pockets within the protein. Nucleotide flipping disrupts Watson-Crick base pairing, induces a kink in the DNA and shifts the DNA register by 1 bp, making the distances between scissile phosphates in the Ecl18kI and NgoMIV cocrystal structures nearly identical. Therefore, the two enzymes can use a conserved DNA recognition module, yet recognize different sequences, and form superimposable dimers, yet generate different cleavage patterns. Hence, Ecl18kI is the first example of a restriction endonuclease that flips nucleotides to achieve specificity for its recognition site.

Nucleotide flips determine the specificity of the Ecl18kI restriction endonuclease.

The Gumbel-Fisher-Tippett (GFT) extreme-value analysis is applied to evaluate the distribution, expectation value and standard deviation of the intensity J of the strongest reflection in a given resolution shell in the X-ray diffraction pattern of a crystal with many scattering atoms in the unit cell. For convenience, intensities are measured in units of the average reflection intensity in the resolution shell and, for simplicity, centric and acentric reflections are treated separately. For acentric reflections, the expectation value mu and standard deviation sigma of J are mu = ln n + gamma and sigma = pi/6(1/2), where n is the number of crystallographically independent reflections in the resolution shell and gamma approximately 0.577 is the Euler-Mascheroni constant. For centric reflections with expectation value 1 for the intensity, the corresponding expressions are mu = 2(ln n + gamma) - ln(pi ln n) and sigma = 2pi/6(1/2) - pi/(6(1/2) ln n). Extensive numerical simulations show that these formulas are excellent approximations for random atom configurations at all resolutions, and good approximations for real protein crystal structures in the resolution range between 2.5 and 1.0 A.

The highest reflection intensity in a resolution shell.

Phytopathogen transcription activator-like effectors (TALEs) bind DNA in a sequence specific manner in order to manipulate host transcription. TALE specificity correlates with repeat variable diresidues in otherwise highly stereotypical 34-35mer repeats. Recently, the crystal structures of two TALE DNA-binding domains have illustrated the molecular basis of the TALE cipher. The structures show that the TALE repeats form a right-handed superhelix that is wound around largely undistorted B-DNA to match its helical parameters. Surprisingly, repeat variable residue 1 is not in contact with the bases. Instead, it is involved in hydrogen bonding interactions that stabilize the overall structure of the protein. Repeat variable residue 2 contacts the top strand base and forms sequence-specific hydrogen bonds and/or van der Waals contacts. Very unexpectedly, bottom strand bases are exposed to solvent and do not make any direct contacts with the protein. This review contains a summary of TALE biology and applications and a detailed description of the recent breakthroughs that have provided insights into the molecular basis of the TALE code.

Structural basis of the TAL effector-DNA interaction.

PvuRts1I is a prototype for a larger family of restriction endonucleases that cleave DNA containing 5-hydroxymethylcytosine (5hmC) or 5-glucosylhydroxymethylcytosine (5ghmC), but not 5-methylcytosine (5mC) or cytosine. Here, we report a crystal structure of the enzyme at 2.35 Å resolution. Although the protein has been crystallized in the absence of DNA, the structure is very informative. It shows that PvuRts1I consists of an N-terminal, atypical PD-(D/E)XK catalytic domain and a C-terminal SRA domain that might accommodate a flipped 5hmC or 5ghmC base. Changes to predicted catalytic residues of the PD-(D/E)XK domain or to the putative pocket for a flipped base abolish catalytic activity. Surprisingly, fluorescence changes indicative of base flipping are not observed when PvuRts1I is added to DNA substrates containing pyrrolocytosine in place of 5hmC (5ghmC). Despite this caveat, the structure suggests a model for PvuRts1I activity and presents opportunities for protein engineering to alter the enzyme properties for biotechnological applications. 

Crystal structure of the 5hmC specific endonuclease PvuRts1I.

R.DpnI consists of N-terminal catalytic and C-terminal winged helix domains that are separately specific for the Gm6ATC sequences in Dam-methylated DNA. Here we present a crystal structure of R.DpnI with oligoduplexes bound to the catalytic and winged helix domains and identify the catalytic domain residues that are involved in interactions with the substrate methyl groups. We show that these methyl groups in the Gm6ATC target sequence are positioned very close to each other. We further show that the presence of the two methyl groups requires a deviation from B-DNA conformation to avoid steric conflict. The methylation compatible DNA conformation is complementary with binding sites of both R.DpnI domains. This indirect readout of methylation adds to the specificity mediated by direct favorable interactions with the methyl groups and solvation/desolvation effects. We also present hydrogen/deuterium exchange data that support 'crosstalk' between the two domains in the identification of methylated DNA, which should further enhance R.DpnI methylation specificity. 

Structural basis of the methylation specificity of R.DpnI.

In mammals, a family of TET enzymes producing oxidized forms of 5-methylcytosine (5mC) plays an important role in modulating DNA demethylation dynamics. In contrast, nothing is known about the function of a single TET orthologue present in invertebrates. Here, we show that the honeybee TET (AmTET) catalytic domain has dioxygenase activity and converts 5mC to 5-hydroxymethylcytosine (5hmC) in a HEK293T cell assay. In vivo, the levels of 5hmC are condition-dependent and relatively low, but in testes and ovaries 5hmC is present at approximately 7-10% of the total level of 5mC, which is comparable to that reported for certain mammalian cells types. AmTET is alternatively spliced and highly expressed throughout development and in adult tissues with the highest expression found in adult brains. Our findings reveal an additional level of flexible genomic modifications in the honeybee that may be important for the selection of multiple pathways controlling contrasting phenotypic outcomes in this species. In a broader context, our study extends the current, mammalian-centred attention to TET-driven DNA hydroxymethylation to an easily manageable organism with attractive and unique biology. 

Insights into DNA hydroxymethylation in the honeybee from in-depth analyses of TET dioxygenase.

In this issue of Structure, Lambert et al. (2016) describe extensive structural and functional work on meganucleases, the group of homing endonucleases most commonly adapted to genome engineering applications. The data are of interest to structural biologists, evolutionary biologists, protein designers, and genome engineers.

Indirect DNA Sequence Readout by LAGLIDADG Homing Endonucleases.

Thymine DNA glycosylase (TDG) and Nei-like 1 (NEIL1) have both been implicated in the base excision repair step of active DNA demethylation. The robust glycosylase activity of TDG on DNA substrates containing 5-formylcytosine (5fC) or 5-carboxylcytosine (5caC) is universally accepted, but the mode of action of NEIL1 is still debated. Based on genetic experiments, it has been suggested that NEIL1 acts redundantly with TDG and excises 5fC and 5caC directly. However, this result has been disputed, and it was suggested instead that NEIL1 is recruited by the monofunctional TDG for the 2'-deoxyribose excision step. Using purified human NEIL1 and its catalytically impaired P2T and E3Q variants as controls, we detect NEIL1 activity on 5caC, but not a 5fC containing dsDNA substrate. We confirm direct NEIL1 TDG binding and NEIL1 mediated 2'-deoxyribose excision downstream of TDG glycosylase activity. NEIL1 acts not only downstream of TDG, but also enhances TDG activity on 5fC or 5caC containing DNA. NEIL1 mediated enhancement of the TDG glycosylase activity is substrate specific and does not occur for dsDNA with a T/G mismatch.

Nei-like 1 (NEIL1) excises 5-carboxylcytosine directly and stimulates TDG-mediated 5-formyl and 5-carboxylcytosine excision.

Ancestral β-subunit (Anbu) is homologous to HslV and 20S proteasomes. Based on its phylogenetic distribution and sequence clustering, Anbu has been proposed as the "ancestral" form of proteasomes. Here, we report biochemical data, small-angle X-ray scattering results, negative-stain electron microscopy micrographs and a crystal structure of the Anbu particle from Yersinia bercovieri (YbAnbu). All data are consistent with YbAnbu forming defined 12-14 subunit multimers that differ in shape from both HslV and 20S proteasomes. The crystal structure reveals that YbAnbu subunits form tight dimers, held together in part by the Anbu specific C-terminal helices. These dimers ("protomers") further assemble into a low-rise left-handed staircase. The lock-washer shape of YbAnbu is consistent with the presence of defined multimers, X-ray diffraction data in solution and negative-stain electron microscopy images. The presented structure suggests a possible evolutionary pathway from helical filaments to highly symmetric or pseudosymmetric multimer structures. YbAnbu subunits have the Ntn-hydrolase fold, a putative S pocket and conserved candidate catalytic residues Thr1, Asp17 and Lys32(33). Nevertheless, we did not detect any YbAnbu peptidase or amidase activity. However, we could document orthophosphate production from ATP catalyzed by the ATP-grasp protein encoded in the Y. bercovieri Anbu operon.

The Y. bercovieri Anbu crystal structure sheds light on the evolution of highly (pseudo)symmetric multimers.

CXXC domains have traditionally been considered as CpG specific DNA binding domains that are repelled by cytosine modifications. This view has recently been challenged by the demonstration that CXXC domain of TET3 has relaxed sequence specificity and binds with the highest affinity to symmetric DNA duplex containing 5caCpG. Here, we present a comparative analysis of the MLL1-CXXC and TET3-CXXC sequence specificity and tolerance to cytosine modifications (5-methyl, 5-hydroxymethyl, 5-formyl, 5-carboxyl) in CpG and non-CpG context. For the first time, we take into consideration possible interference from cytosine bases elsewhere in the sequence. We show that despite similar overall structure, MLL1-CXXC has greater sequence and modification specificity than TET3-CXXC. MLL1-CXXC is specific only for CpG and does not tolerate any cytosine modifications. In contrast, TET3-CXXC does not require the CpG context of cytosine bases. Methyl-, formyl- and carboxyl-modifications are tolerated by TET3-CXXC, but only preceding G. Based on our and other data we propose a parsimonious model of MLL1-CXXC and TET3-CXXC DNA binding. This model explains why the binding of modified DNA duplexes by TET3-CXXC requires in some cases a register shift and is therefore context-dependent.

Matthias Bochtler

International Institute of Molecular and Cell Biology, Warsaw (IIMCB)

In Vitro Directed Evolution of a Restriction Endonuclease with More Stringent Specificity