Cytosine analogues as DNA methyltransferase substrates.

DNA methyltransferases are drug targets for myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML), acute myelogenous leukemia (AML) and possibly ß-hemoglobinopathies. We characterize the interaction of nucleoside analogues in DNA with a prokaryotic CpG-specific DNA methyltransferase (M.MpeI) as a model for mammalian DNMT1 methyltransferases. We tested DNA containing 5-hydroxymethylcytosine (5hmC), 5-hydroxycytosine (5OHC), 5-methyl-2-pyrimidinone (in the ribosylated form known as 5-methylzebularine, 5mZ), 5,6-dihydro-5-azacytosine (dhaC), 5-fluorocytosine (5FC), 5-chlorocytosine (5ClC), 5-bromocytosine (5BrC) and 5-iodocytosine (5IC). Covalent complex formation was by far most efficient for 5FC. Non-covalent complexes were most abundant for dhaC and 5mZ. Surprisingly, we observed methylation of 5IC and 5BrC, and to a lesser extent 5ClC and 5FC, in the presence, but not the absence of small molecule thiol nucleophiles. For 5IC and 5BrC, we demonstrated by mass spectrometry that the reactions were due to methyltransferase driven dehalogenation, followed by methylation. Crystal structures of M.MpeI-DNA complexes capture the 'in' conformation of the active site loop for analogues with small or rotatable (5mZ) 5-substituents and its 'out' form for bulky 5-substituents. Since very similar 'in' and 'out' loop conformations were also observed for DNMT1, it is likely that our conclusions generalize to other DNA methyltransferases.

Crystal structures of the EVE-HNH endonuclease VcaM4I in the presence and absence of DNA.

Many modification-dependent restriction endonucleases (MDREs) are fusions of a PUA superfamily modification sensor domain and a nuclease catalytic domain. EVE domains belong to the PUA superfamily, and are present in MDREs in combination with HNH nuclease domains. Here, we present a biochemical characterization of the EVE-HNH endonuclease VcaM4I and crystal structures of the protein alone, with EVE domain bound to either 5mC modified dsDNA or to 5mC/5hmC containing ssDNA. The EVE domain is moderately specific for 5mC/5hmC containing DNA according to EMSA experiments. It flips the modified nucleotide, to accommodate it in a hydrophobic pocket of the enzyme, primarily formed by P24, W82 and Y130 residues. In the crystallized conformation, the EVE domain and linker helix between the two domains block DNA binding to the catalytic domain. Removal of the EVE domain and inter-domain linker, but not of the EVE domain alone converts VcaM4I into a non-specific toxic nuclease. The role of the key residues in the EVE and HNH domains of VcaM4I is confirmed by digestion and restriction assays with the enzyme variants that differ from the wild-type by changes to the base binding pocket or to the catalytic residues.

Restriction endonucleases that cleave RNA/DNA heteroduplexes bind dsDNA in A-like conformation.

Restriction endonucleases naturally target DNA duplexes. Systematic screening has identified a small minority of these enzymes that can also cleave RNA/DNA heteroduplexes and that may therefore be useful as tools for RNA biochemistry. We have chosen AvaII (G↓GWCC, where W stands for A or T) as a representative of this group of restriction endonucleases for detailed characterization. Here, we report crystal structures of AvaII alone, in specific complex with partially cleaved dsDNA, and in scanning complex with an RNA/DNA hybrid. The specific complex reveals a novel form of semi-specific dsDNA readout by a hexa-coordinated metal cation, most likely Ca2+ or Mg2+. Substitutions of residues anchoring this non-catalytic metal ion severely impair DNA binding and cleavage. The dsDNA in the AvaII complex is in the A-like form. This creates space for 2'-OH groups to be accommodated without intra-nucleic acid steric conflicts. PD-(D/E)XK restriction endonucleases of known structure that bind their dsDNA targets in the A-like form cluster into structurally similar groups. Most such enzymes, including some not previously studied in this respect, cleave RNA/DNA heteroduplexes. We conclude that A-form dsDNA binding is a good predictor for RNA/DNA cleavage activity.

The Nϵ-Rule for Serine, but Not Cysteine Catalytic Triads.

Catalytic triads, composed of a serine or cysteine nucleophile, a histidine, and a third triad residue (typically Asp/Glu/Asn), are common in enzyme active sites and catalyze a wide variety of chemical reactions. Two types of triads can be distinguished: We refer to them as Nδ- or Nϵ-configured, depending on whether the histidine imidazole Nδ or Nϵ atom is close to the nucleophile Oγ/Sγ. In this study, we have analyzed triad configuration. In structural triads, the more stable Nδ-configuration predominates. For catalytic triads, the configuration depends on the nucleophile. When it is a cysteine residue, both configuration types occur, depending on the family. However, when the nucleophile is a serine residue, the less stable Nϵ-configuration is almost exclusively found. We posit that the energetically less favored conformation is selected for in serine triads to facilitate the otherwise difficult proton transfer from the nucleophile to the histidine residue.

Crystal structure of the EcoKMcrA N-terminal domain (NEco): recognition of modified cytosine bases without flipping.

EcoKMcrA from Escherichia coli restricts CpG methylated or hydroxymethylated DNA, and may act as a barrier against host DNA. The enzyme consists of a novel N-terminal specificity domain that we term NEco, and a C-terminal catalytic HNH domain. Here, we report that NEco and full-length EcoKMcrA specificities are consistent. NEco affinity to DNA increases more from hemi- to full-methylation than from non- to hemi-methylation, indicating cooperative binding of the methyl groups. We determined the crystal structures of NEco in complex with fully modified DNA containing three variants of the Y5mCGR EcoKMcrA target sequence: C5mCGG, T5mCGA and T5hmCGA. The structures explain the specificity for the two central base pairs and one of the flanking pairs. As predicted based on earlier biochemical experiments, NEco does not flip any DNA bases. The proximal and distal methyl groups are accommodated in separate pockets. Changes to either pocket reduce DNA binding by NEco and restriction by EcoKMcrA, confirming the relevance of the crystallographically observed binding mode in solution.

A protein architecture guided screen for modification dependent restriction endonucleases.

Modification dependent restriction endonucleases (MDREs) often have separate catalytic and modification dependent domains. We systematically looked for previously uncharacterized fusion proteins featuring a PUA or DUF3427 domain and HNH or PD-(D/E)XK catalytic domain. The enzymes were clustered by similarity of their putative modification sensing domains into several groups. The TspA15I (VcaM4I, CmeDI), ScoA3IV (MsiJI, VcaCI) and YenY4I groups, all featuring a PUA superfamily domain, preferentially cleaved DNA containing 5-methylcytosine or 5-hydroxymethylcytosine. ScoA3V, also featuring a PUA superfamily domain, but of a different clade, exhibited 6-methyladenine stimulated nicking activity. With few exceptions, ORFs for PUA-superfamily domain containing endonucleases were not close to DNA methyltransferase ORFs, strongly supporting modification dependent activity of the endonucleases. DUF3427 domain containing fusion proteins had very little or no endonuclease activity, despite the presence of a putative PD-(D/E)XK catalytic domain. However, their expression potently restricted phage T4gt in Escherichia coli cells. In contrast to the ORFs for PUA domain containing endonucleases, the ORFs for DUF3427 fusion proteins were frequently found in defense islands, often also featuring DNA methyltransferases.

Crystal structure of the modification-dependent SRA-HNH endonuclease TagI.

TagI belongs to the recently characterized SRA-HNH family of modification-dependent restriction endonucleases (REases) that also includes ScoA3IV (Sco5333) and TbiR51I (Tbis1). Here, we present a crystal structure of dimeric TagI, which exhibits a DNA binding site formed jointly by the nuclease domains, and separate binding sites for modified DNA bases in the two protomers. The nuclease domains have characteristic features of HNH/ßßα-Me REases, and catalyze nicks or double strand breaks, with preference for /RY and RYN/RY sites, respectively. The SRA domains have the canonical fold. Their pockets for the flipped bases are spacious enough to accommodate 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC), but not glucosyl-5-hydroxymethylcytosine (g5hmC). Such preference is in agreement with the biochemical determination of the TagI modification dependence and the results of phage restriction assays. The ability of TagI to digest plasmids methylated by Dcm (C5mCWGG), M.Fnu4HI (G5mCNGC) or M.HpyCH4IV (A5mCGT) suggests that the SRA domains of the enzyme are tolerant to different sequence contexts of the modified base.

Activity and structure of EcoKMcrA.

Escherichia coli McrA (EcoKMcrA) acts as a methylcytosine and hydroxymethylcytosine dependent restriction endonuclease. We present a biochemical characterization of EcoKMcrA that includes the first demonstration of its endonuclease activity, small angle X-ray scattering (SAXS) data, and a crystal structure of the enzyme in the absence of DNA. Our data indicate that EcoKMcrA dimerizes via the anticipated C-terminal HNH domains, which together form a single DNA binding site. The N-terminal domains are not homologous to SRA domains, do not interact with each other, and have separate DNA binding sites. Electrophoretic mobility shift assay (EMSA) and footprinting experiments suggest that the N-terminal domains can sense the presence and sequence context of modified cytosines. Pyrrolocytosine fluorescence data indicate no base flipping. In vitro, EcoKMcrA DNA endonuclease activity requires Mn2+ ions, is not strictly methyl dependent, and is not observed when active site variants of the enzyme are used. In cells, EcoKMcrA specifically restricts DNA that is modified in the correct sequence context. This activity is impaired by mutations of the nuclease active site, unless the enzyme is highly overexpressed.

On the role of steric clashes in methylation control of restriction endonuclease activity.

Restriction-modification systems digest non-methylated invading DNA, while protecting host DNA against the endonuclease activity by methylation. It is widely believed that the methylated DNA would not 'fit' into the binding site of the endonuclease in the productive orientation, and thus steric clashes should account for most of the protection. We test this concept statistically by grafting methyl groups in silico onto non-methylated DNA in co-crystal structures with restriction endonucleases. Clash scores are significantly higher for protective than non-protective methylation (P < 0.05% according to the Wilcoxon rank sum test). Structural data alone are sufficient to distinguish between protective and non-protective DNA methylation with 90% confidence and decision thresholds of 1.1 Å and 48 Å(3) for the most severe distance-based and cumulative volume-based clash with the protein, respectively (0.1 Å was deducted from each interatomic distance to allow for coordinate errors). The most severe clashes are more pronounced for protective methyl groups attached to the nitrogen atoms (N6-methyladenines and N4-methylcytosines) than for C5-methyl groups on cytosines. Cumulative clashes are comparable for all three types of protective methylation.

Crystal structure of the 5hmC specific endonuclease PvuRts1I.

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.

Structural basis of the methylation specificity of R.DpnI.

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.

CpG underrepresentation and the bacterial CpG-specific DNA methyltransferase M.MpeI.

Cytosine methylation promotes deamination. In eukaryotes, CpG methylation is thought to account for CpG underrepresentation. Whether scarcity of CpGs in prokaryotic genomes is diagnostic for methylation is not clear. Here, we report that Mycoplasms tend to be CpG depleted and to harbor a family of constitutively expressed or phase variable CpG-specific DNA methyltransferases. The very CpG poor Mycoplasma penetrans and its constitutively active CpG-specific methyltransferase M.MpeI were chosen for further characterization. Genome-wide sequencing of bisulfite-converted DNA indicated that M.MpeI methylated CpG target sites both in vivo and in vitro in a locus-nonselective manner. A crystal structure of M.MpeI with DNA at 2.15-Å resolution showed that the substrate base was flipped and that its place in the DNA stack was taken by a glutamine residue. A phenylalanine residue was intercalated into the "weak" CpG step of the nonsubstrate strand, indicating mechanistic similarities in the recognition of the short CpG target sequence by prokaryotic and eukaryotic DNA methyltransferases.

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

Importance of single molecular determinants in the fidelity of expanded genetic codes.

The site-selective encoding of noncanonical amino acids (NAAs) is a powerful technique for the installation of novel chemical functional groups in proteins. This is often achieved by recoding a stop codon and requires two additional components: an evolved aminoacyl tRNA synthetase (AARS) and a cognate tRNA. Analysis of the most successful AARSs reveals common characteristics. The highest fidelity NAA systems derived from the Methanocaldococcus jannaschii tyrosyl AARS feature specific mutations to two residues reported to interact with the hydroxyl group of the substrate tyrosine. We demonstrate that the restoration of just one of these determinants for amino acid specificity results in the loss of fidelity as the evolved AARSs become noticeably promiscuous. These results offer a partial explanation of a recently retracted strategy for the synthesis of glycoproteins. Similarly, we reinvestigated a tryptophanyl AARS reported to allow the site-selective incorporation of 5-hydroxy tryptophan within mammalian cells. In multiple experiments, the enzyme displayed elements of promiscuity despite its previous characterization as a high fidelity enzyme. Given the many similarities of the TyrRSs and TrpRSs reevaluated here, our findings can be largely combined, and in doing so they reinforce the long-established central dogma regarding the molecular basis by which these enzymes contribute to the fidelity of translation. Thus, our view is that the central claims of fidelity reported in several NAA systems remain unproven and unprecedented.

RrA, an enzyme from Rhodospirillum rubrum, is a prototype of a new family of short-chain L-asparaginases.

L-Asparaginases (ASNases) catalyze the hydrolysis of L-Asn to L-Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300-400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum, reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N-terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L-Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short-chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild-type RrA, both with and without bound L-Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand-alone subfamily of ASNases that can catalyze the hydrolysis of L-Asn to L-Asp despite the lack of the C-terminal domain that is present in all ASNases described structurally to date.

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.

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.

Structural basis for efficient chromophore communication and energy transfer in a constructed didomain protein scaffold.

The construction of useful functional biomolecular components not currently part of the natural repertoire is central to synthetic biology. A new light-capturing ultra-high-efficiency energy transfer protein scaffold has been constructed by coupling the chromophore centers of two normally unrelated proteins: the autofluorescent protein enhanced green fluorescent protein (EGFP) and the heme-binding electron transfer protein cytochrome b(562) (cyt b(562)). Using a combinatorial domain insertion strategy, a variant was isolated in which resonance energy transfer from the donor EGFP to the acceptor cyt b(562) was close to 100% as evident by virtually full fluorescence quenching on heme binding. The fluorescence signal of the variant was also sensitive to the reactive oxygen species H(2)O(2), with high signal gain observed due to the release of heme. The structure of oxidized holoprotein, determined to 2.75 Å resolution, revealed that the two domains were arranged side-by-side in a V-shape conformation, generating an interchromophore distance of ~17 Å (14 Å edge-to-edge). Critical to domain arrangement is the formation of a molecular pivot point between the two domains as a result of different linker sequence lengths at each domain junction and formation of a predominantly polar interdomain interaction surface. The retrospective structural analysis has provided an explanation for the basis of the observed highly efficient energy transfer through chromophore arrangement in the directly evolved protein scaffold and provides an insight into the molecular principles by which to design new proteins with coupled functions.

Competition between electrostatic interactions and halogen bonding in the protein-ligand system: structural and thermodynamic studies of 5,6-dibromobenzotriazole-hCK2α complexes.

CK2 is a member of the CMGC group of eukaryotic protein kinases and a cancer drug target. It can be efficiently inhibited by halogenated benzotriazoles and benzimidazoles. Depending on the scaffold, substitution pattern, and pH, these compounds are either neutral or anionic. Their binding poses are dictated by a hydrophobic effect (desolvation) and a tug of war between a salt bridge/hydrogen bond (to K68) and halogen bonding (to E114 and V116 backbone oxygens). Here, we test the idea that binding poses might be controllable by pH for ligands with near-neutral pKa, using the conditionally anionic 5,6-DBBt and constitutively anionic TBBt as our models. We characterize the binding by low-volume Differential Scanning Fluorimetry (nanoDSF), Isothermal Calorimetry (ITC), Hydrogen/Deuterium eXchange (HDX), and X-ray crystallography (MX). The data indicate that the ligand pose away from the hinge dominates for the entire tested pH range (5.5-8.5). The insensitivity of the binding mode to pH is attributed to the perturbation of ligand pKa upon binding that keeps it anionic in the ligand binding pocket at all tested pH values. However, a minor population of the ligand, detectable only by HDX, shifts towards the hinge in acidic conditions. Our findings demonstrate that electrostatic (ionic) interactions predominate over halogen bonding.

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

The beta beta alpha-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 beta-barrel and two beta 4 alpha 2 repeats. Each repeat coordinates a structural zinc ion with four cysteine thiolates in two CXXC motifs. The beta beta alpha-Me region of the second beta 4 alpha 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 beta 4 alpha 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 beta beta alpha-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).