A Role for N6-Methyladenine in DNA Damage Repair.

The leading cause of mutation due to oxidative damage is 8-oxo-2'-deoxyguanosine (8-oxoG) mispairing with adenine (Ade), which can occur in two ways. First, guanine of a G:C DNA base pair can be oxidized. If not repaired in time, DNA polymerases can mispair Ade with 8-oxoG in the template. This 8-oxoG:A can be repaired by enzymes that remove Ade opposite to template 8-oxoG, or 8-oxoG opposite to Cyt. Second, free 8-oxo-dGTP can be misincorporated by DNA polymerases into DNA opposite template Ade. However, there is no known repair activity that removes 8-oxoG opposite to template Ade. We suggest that a major role of N6-methyladenine in mammalian DNA is minimizing incorporation of 8-oxoG opposite to Ade by DNA polymerases following adduct formation.

Structures of CTCF-DNA complexes including all 11 zinc fingers.

The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12-15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1-ZF7 and ZF3-ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1-ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9-ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1-ZF7 and ZF8-ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1-ZF7 and bulkier and polar/charged residues for ZF8-ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine-guanine and glutamine-adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9-ZF11. The differences between ZF1-ZF7 and ZF8-ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.

Allosteric autoregulation of DNA binding via a DNA-mimicking protein domain: a biophysical study of ZNF410-DNA interaction using small angle X-ray scattering.

ZNF410 is a highly-conserved transcription factor, remarkable in that it recognizes a 15-base pair DNA element but has just a single responsive target gene in mammalian erythroid cells. ZNF410 includes a tandem array of five zinc-fingers (ZFs), surrounded by uncharacterized N- and C-terminal regions. Unexpectedly, full-length ZNF410 has reduced DNA binding affinity, compared to that of the isolated DNA binding ZF array, both in vitro and in cells. AlphaFold predicts a partially-folded N-terminal subdomain that includes a 30-residue long helix, preceded by a hairpin loop rich in acidic (aspartate/glutamate) and serine/threonine residues. This hairpin loop is predicted by AlphaFold to lie against the DNA binding interface of the ZF array. In solution, ZNF410 is a monomer and binds to DNA with 1:1 stoichiometry. Surprisingly, the single best-fit model for the experimental small angle X-ray scattering profile, in the absence of DNA, is the original AlphaFold model with the N-terminal long-helix and the hairpin loop occupying the ZF DNA binding surface. For DNA binding, the hairpin loop presumably must be displaced. After combining biophysical, biochemical, bioinformatic and artificial intelligence-based AlphaFold analyses, we suggest that the hairpin loop mimics the structure and electrostatics of DNA, and provides an additional mechanism, supplementary to sequence specificity, of regulating ZNF410 DNA binding.

Biochemical and structural characterization of the first-discovered metazoan DNA cytosine-N4 methyltransferase from the bdelloid rotifer Adineta vaga.

Much is known about the generation, removal, and roles of 5-methylcytosine (5mC) in eukaryote DNA, and there is a growing body of evidence regarding N6-methyladenine, but very little is known about N4-methylcytosine (4mC) in the DNA of eukaryotes. The gene for the first metazoan DNA methyltransferase generating 4mC (N4CMT) was reported and characterized recently by others, in tiny freshwater invertebrates called bdelloid rotifers. Bdelloid rotifers are ancient, apparently asexual animals, and lack canonical 5mC DNA methyltransferases. Here, we characterize the kinetic properties and structural features of the catalytic domain of the N4CMT protein from the bdelloid rotifer Adineta vaga. We find that N4CMT generates high-level methylation at preferred sites, (a/c)CG(t/c/a), and low-level methylation at disfavored sites, exemplified by ACGG. Like the mammalian de novo 5mC DNA methyltransferase 3A/3B (DNMT3A/3B), N4CMT methylates CpG dinucleotides on both DNA strands, generating hemimethylated intermediates and eventually fully methylated CpG sites, particularly in the context of favored symmetric sites. In addition, like DNMT3A/3B, N4CMT methylates non-CpG sites, mainly CpA/TpG, though at a lower rate. Both N4CMT and DNMT3A/3B even prefer similar CpG-flanking sequences. Structurally, the catalytic domain of N4CMT closely resembles the Caulobacter crescentus cell cycle-regulated DNA methyltransferase. The symmetric methylation of CpG, and similarity to a cell cycle-regulated DNA methyltransferase, together suggest that N4CMT might also carry out DNA synthesis-dependent methylation following DNA replication.

Structural basis for transcription factor ZBTB7A recognition of DNA and effects of ZBTB7A somatic mutations that occur in human acute myeloid leukemia.

ZBTB7A belongs to a small family of transcription factors having three members in humans (7A, 7B, and 7C). They share a BTB/POZ protein interaction domain at the amino end and a zinc-finger DNA-binding domain at the carboxyl end. They control the transcription of a wide range of genes, having varied functions in hematopoiesis, oncogenesis, and metabolism (in particular glycolysis). ZBTB7A-binding profiles at gene promoters contain a consensus G(a/c)CCC motif, followed by a CCCC sequence in some instances. Structural and mutational investigations suggest that DNA-specific contacts with the four-finger tandem array of ZBTB7A are formed sequentially, initiated from ZF1-ZF2 binding to G(a/c)CCC before spreading to ZF3-ZF4, which bind the DNA backbone and the 3' CCCC sequence, respectively. Here, we studied some mutations found in t(8;21)-positive acute myeloid leukemia patients that occur within the ZBTB7A DNA-binding domain. We determined that these mutations generally impair ZBTB7A DNA binding, with the most severe disruptions resulting from mutations in ZF1 and ZF2, and the least from a frameshift mutation in ZF3 that results in partial mislocalization. Information provided here on ZBTB7A-DNA interactions is likely applicable to ZBTB7B/C, which have overlapping functions with ZBTB7A in controlling primary metabolism.

A growth-based assay using fluorescent protein emission to screen for S-adenosylmethionine synthetase inhibitors.

The use of cell growth-based assays to identify inhibitory compounds is straightforward and inexpensive, but is also inherently insensitive and somewhat nonspecific. To overcome these limitations and develop a sensitive, specific cell-based assay, two different approaches were combined. To address the sensitivity limitation, different fluorescent proteins have been introduced into a bacterial expression system to serve as growth reporters. To overcome the lack of specificity, these protein reporters have been incorporated into a plasmid in which they are paired with different orthologs of an essential target enzyme, in this case l-methionine S-adenosyltransferase (MAT, AdoMet synthetase). Screening compounds that serve as specific inhibitors will reduce the growth of only a subset of strains, because these strains are identical, except for which target ortholog they carry. Screening several such strains in parallel not only reveals potential inhibitors but the strains also serve as specificity controls for one another. The present study makes use of an existing Escherichia coli strain that carries a deletion of metK, the gene for MAT. Transformation with these plasmids leads to a complemented strain that no longer requires externally supplied S-adenosylmethionine for growth, but its growth is now dependent on the activity of the introduced MAT ortholog. The resulting fluorescent strains provide a platform to screen chemical compound libraries and identify species-selective inhibitors of AdoMet synthetases. A pilot study of several chemical libraries using this platform identified new lead compounds that are ortholog-selective inhibitors of this enzyme family, some of which target the protozoal human pathogen Cryptosporidium parvum.

Enzymatic characterization of mRNA cap adenosine-N6 methyltransferase PCIF1 activity on uncapped RNAs.

The phosphorylated RNA polymerase II CTD interacting factor 1 (PCIF1) is a methyltransferase that adds a methyl group to the N6-position of 2'O-methyladenosine (Am), generating N6, 2'O-dimethyladenosine (m6Am) when Am is the cap-proximal nucleotide. In addition, PCIF1 has ancillary methylation activities on internal adenosines (both A and Am), although with much lower catalytic efficiency relative to that of its preferred cap substrate. The PCIF1 preference for 2'O-methylated Am over unmodified A nucleosides is due mainly to increased binding affinity for Am. Importantly, it was recently reported that PCIF1 can methylate viral RNA. Although some viral RNA can be translated in the absence of a cap, it is unclear what roles PCIF1 modifications may play in the functionality of viral RNAs. Here we show, using in vitro assays of binding and methyltransfer, that PCIF1 binds an uncapped 5'-Am oligonucleotide with approximately the same affinity as that of a cap analog (KM = 0.4 versus 0.3 µM). In addition, PCIF1 methylates the uncapped 5'-Am with activity decreased by only fivefold to sixfold compared with its preferred capped substrate. We finally discuss the relationship between PCIF1-catalyzed RNA methylation, shown here to have broader substrate specificity than previously appreciated, and that of the RNA demethylase fat mass and obesity-associated protein (FTO), which demonstrates PCIF1-opposing activities on capped RNAs.

Human MettL3-MettL14 RNA adenine methyltransferase complex is active on double-stranded DNA containing lesions.

MettL3-MettL14 methyltransferase complex has been studied widely for its role in RNA adenine methylation. This complex is also recruited to UV- and X-ray exposed DNA damaged sites, and its methyltransfer activity is required for subsequent DNA repair, though in theory this could result from RNA methylation of short transcripts made at the site of damage. We report here that MettL3-MettL14 is active in vitro on double-stranded DNA containing a cyclopyrimidine dimer - a major lesion of UV radiation-induced products - or an abasic site or mismatches. Furthermore, N6-methyladenine (N6mA) decreases misincorporation of 8-oxo-guanine (8-oxoG) opposite to N6mA by repair DNA polymerases. When 8-oxoG is nevertheless incorporated opposite N6mA, the methylation inhibits N6mA excision from the template (correct) strand by the adenine DNA glycosylase (MYH), implying that the methylation decreases inappropriate misrepair. Finally, we observed that the N6mA reader domain of YTHDC1, which is also recruited to sites of DNA damage, binds N6mA that is located across from a single-base gap between two canonical DNA helices. This YTHDC1 complex with a gapped duplex is structurally similar to DNA complexes with FEN1 and GEN1 - two members of the nuclease family that act in nucleotide excision repair, mismatch repair and homologous recombination, and which incise distinct non-B DNA structures. Together, the parts of our study provide a plausible mechanism for N6mA writer and reader proteins acting directly on lesion-containing DNA, and suggest in vivo experiments to test the mechanisms involving methylation of adenine.

Preferential CEBP binding to T:G mismatches and increased C-to-T human somatic mutations.

DNA cytosine methylation in mammals modulates gene expression and chromatin accessibility. It also impacts mutation rates, via spontaneous oxidative deamination of 5-methylcytosine (5mC) to thymine. In most cases the resulting T:G mismatches are repaired, following T excision by one of the thymine DNA glycosylases, TDG or MBD4. We found that C-to-T mutations are enriched in the binding sites of CCAAT/enhancer binding proteins (CEBP). Within a CEBP site, the presence of a T:G mismatch increased CEBPß binding affinity by a factor of >60 relative to the normal C:G base pair. This enhanced binding to a mismatch inhibits its repair by both TDG and MBD4 in vitro. Furthermore, repair of the deamination product of unmethylated cytosine, which yields a U:G DNA mismatch that is normally repaired via uracil DNA glycosylase, is also inhibited by CEBPß binding. Passage of a replication fork over either a T:G or U:G mismatch, before repair can occur, results in a C-to-T mutation in one of the daughter duplexes. Our study thus provides a plausible mechanism for accumulation of C-to-T human somatic mutations.

Enzymatic Characterization of In Vitro Activity of RNA Methyltransferase PCIF1 on DNA.

PCIF1 and FTO are a pair of human mRNA cap-specific modification enzymes that have opposing activities. PCIF1 adds a methyl group to the N6-position of 2'O-methyladenosine (Am), generating N6, 2'O-dimethyladenosine (m6Am), when Am is the cap-proximal nucleotide. FTO removes the N6-methyl group from m6Am. In addition, FTO has a demethylase activity on a broad spectrum of various RNA substrates, as well as on DNA N6-methyldeoxyadenosine (m6dA). While the existence of m6dA in mammalian DNA remains controversial, we show here that PCIF1 has significant methylation activity on single stranded DNA deoxyadenosine, double stranded RNA/DNA hybrids, and double stranded DNA, though with lower catalytic efficiency than that on its preferred RNA substrate. PCIF1 has activities in the order ssRNA > RNA/DNA hybrid > ssDNA > dsDNA. We discuss the implications of PCIF1 generation, and FTO removal, of DNA adenine methylation.

Enzymatic characterization of three human RNA adenosine methyltransferases reveals diverse substrate affinities and reaction optima.

RNA methylations of varied RNA species (mRNA, tRNA, rRNA, non-coding RNA) generate a range of modified nucleotides, including N6-methyladenosine. Here we study the enzymology of three human RNA methyltransferases that methylate the adenosine amino group in diverse contexts, when it is: the first transcribed nucleotide after the mRNA cap (PCIF1), at position 1832 of 18S rRNA (MettL5-Trm112 complex), and within a hairpin in the 3' UTR of the S-adenosyl-l-methionine synthetase (MettL16). Among these three enzymes, the catalytic efficiency ranges from PCIF1, with the fastest turnover rate of >230 h-1 µM-1 on mRNA cap analog, down to MettL16, which has the lowest rate of ∼3 h-1 µM-1 acting on an RNA hairpin. Both PCIF1 and MettL5 have a binding affinity (Km) of ∼1 µM or less for both substrates of SAM and RNA, whereas MettL16 has significantly lower binding affinities for both (Km >0.4 mM for SAM and ∼10 µM for RNA). The three enzymes are active over a wide pH range (∼5.4-9.4) and have different preferences for ionic strength. Sodium chloride at 200 mM markedly diminished methylation activity of MettL5-Trm112 complex, whereas MettL16 had higher activity in the range of 200 to 500 mM NaCl. Zinc ion inhibited activities of all three enzymes. Together, these results illustrate the diversity of RNA adenosine methyltransferases in their enzymatic mechanisms and substrate specificities and underline the need for assay optimization in their study.

Beta class amino methyltransferases from bacteria to humans: evolution and structural consequences.

S-adenosyl-l-methionine dependent methyltransferases catalyze methyl transfers onto a wide variety of target molecules, including DNA and RNA. We discuss a family of methyltransferases, those that act on the amino groups of adenine or cytosine in DNA, have conserved motifs in a particular order in their amino acid sequence, and are referred to as class beta MTases. Members of this class include M.EcoGII and M.EcoP15I from Escherichia coli, Caulobacter crescentus cell cycle-regulated DNA methyltransferase (CcrM), the MTA1-MTA9 complex from the ciliate Oxytricha, and the mammalian MettL3-MettL14 complex. These methyltransferases all generate N6-methyladenine in DNA, with some members having activity on single-stranded DNA as well as RNA. The beta class of methyltransferases has a unique multimeric feature, forming either homo- or hetero-dimers, allowing the enzyme to use division of labor between two subunits in terms of substrate recognition and methylation. We suggest that M.EcoGII may represent an ancestral form of these enzymes, as its activity is independent of the nucleic acid type (RNA or DNA), its strandedness (single or double), and its sequence (aside from the target adenine).

Biochemical and structural basis for YTH domain of human YTHDC1 binding to methylated adenine in DNA.

The recently characterized mammalian writer (methyltransferase) and eraser (demethylase) of the DNA N6-methyladenine (N6mA) methyl mark act on single-stranded (ss) and transiently-unpaired DNA. As YTH domain-containing proteins bind N6mA-containing RNA in mammalian cells, we investigated whether mammalian YTH domains are also methyl mark readers of N6mA DNA. Here, we show that the YTH domain of YTHDC1 (known to localize in the nucleus) binds ssDNA containing N6mA, with a 10 nM dissociation constant. This binding is stronger by a factor of 5 than in an RNA context, tested under the same conditions. However, the YTH domains of YTHDF2 and YTHDF1 (predominantly cytoplasmic) exhibited the opposite effect with ∼1.5-2נstronger binding to ssRNA containing N6mA than to the corresponding DNA. We determined two structures of the YTH domain of YTHDC1 in complex with N6mA-containing ssDNA, which illustrated that YTHDC1 binds the methylated adenine in a single-stranded region flanked by duplexed DNA. We discuss the hypothesis that the writer-reader-eraser of N6mA-containining ssDNA is associated with maintaining genome stability. Structural comparison of YTH and SRA domains (the latter a DNA 5-methylcytosine reader) revealed them to be diverse members of a larger family of DNA/RNA modification readers, apparently having originated from bacterial modification-dependent restriction enzymes.

Structural basis for preferential binding of human TCF4 to DNA containing 5-carboxylcytosine.

The psychiatric risk-associated transcription factor 4 (TCF4) is linked to schizophrenia. Rare TCF4 coding variants are found in individuals with Pitt-Hopkins syndrome-an intellectual disability and autism spectrum disorder. TCF4 contains a C-terminal basic-helix-loop-helix (bHLH) DNA binding domain which recognizes the enhancer-box (E-box) element 5'-CANNTG-3' (where N = any nucleotide). A subset of the TCF4-occupancy sites have the expanded consensus binding specificity 5'-C(A/G)-CANNTG-3', with an added outer Cp(A/G) dinucleotide; for example in the promoter for CNIH3, a gene involved in opioid dependence. In mammalian genomes, particularly brain, the CpG and CpA dinucleotides can be methylated at the 5-position of cytosine (5mC), and then may undergo successive oxidations to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), and 5-carboxyl (5caC) forms. We find that, in the context of 5'-0CG-1CA-2CG-3TG-3'(where the numbers indicate successive dinucleotides), modification of the central E-box 2CG has very little effect on TCF4 binding, E-box 1CA modification has a negative influence on binding, while modification of the flanking 0CG, particularly carboxylation, has a strong positive impact on TCF4 binding to DNA. Crystallization of TCF4 in complex with unmodified or 5caC-modified oligonucleotides revealed that the basic region of bHLH domain adopts multiple conformations, including an extended loop going through the DNA minor groove, or the N-terminal portion of a long helix binding in the DNA major groove. The different protein conformations enable arginine 576 (R576) to interact, respectively, with a thymine in the minor groove, a phosphate group of DNA backbone, or 5caC in the major groove. The Pitt-Hopkins syndrome mutations affect five arginine residues in the basic region, two of them (R569 and R576) involved in 5caC recognition. Our analyses indicate, and suggest a structural basis for, the preferential recognition of 5caC by a transcription factor centrally important in brain development.

Structural basis for effects of CpA modifications on C/EBPß binding of DNA.

CCAAT/enhancer binding proteins (C/EBPs) regulate gene expression in a variety of cells/tissues/organs, during a range of developmental stages, under both physiological and pathological conditions. C/EBP-related transcription factors have a consensus binding specificity of 5'-TTG-CG-CAA-3', with a central CpG/CpG and two outer CpA/TpG dinucleotides. Methylation of the CpG and CpA sites generates a DNA element with every pyrimidine having a methyl group in the 5-carbon position (thymine or 5-methylcytosine (5mC)). To understand the effects of both CpG and CpA modification on a centrally-important transcription factor, we show that C/EBPß binds the methylated 8-bp element with modestly-increased (2.4-fold) binding affinity relative to the unmodified cognate sequence, while cytosine hydroxymethylation (particularly at the CpA sites) substantially decreased binding affinity (36-fold). The structure of C/EBPß DNA binding domain in complex with methylated DNA revealed that the methyl groups of the 5mCpA/TpG make van der Waals contacts with Val285 in C/EBPß. Arg289 recognizes the central 5mCpG by forming a methyl-Arg-G triad, and its conformation is constrained by Val285 and the 5mCpG methyl group. We substituted Val285 with Ala (V285A) in an Ala-Val dipeptide, to mimic the conserved Ala-Ala in many members of the basic leucine-zipper family of transcription factors, important in gene regulation, cell proliferation and oncogenesis. The V285A variant demonstrated a 90-fold binding preference for methylated DNA (particularly 5mCpA methylation) over the unmodified sequence. The smaller side chain of Ala285 permits Arg289 to adopt two alternative conformations, to interact in a similar fashion with either the central 5mCpG or the TpG of the opposite strand. Significantly, the best-studied cis-regulatory elements in RNA polymerase II promoters and enhancers have variable sequences corresponding to the central CpG or reduced to a single G:C base pair, but retain a conserved outer CpA sequence. Our analyses suggest an important modification-dependent CpA recognition by basic leucine-zipper transcription factors.

Distribution of RecBCD and AddAB recombination-associated genes among bacteria in 33 phyla.

Homologous recombination plays key roles in fundamental processes such as recovery from DNA damage and in bacterial horizontal gene transfer, yet there are still open questions about the distribution of recognized components of recombination machinery among bacteria and archaea. RecBCD helicase-nuclease plays a central role in recombination among Gammaproteobacteria like Escherichia coli; while bacteria in other phyla, like the Firmicute Bacillus subtilis, use the related AddAB complex. The activity of at least some of these complexes is controlled by short DNA sequences called crossover hotspot instigator (Chi) sites. When RecBCD or AddAB complexes encounter an autologous Chi site during unwinding, they introduce a nick such that ssDNA with a free end is available to invade another duplex. If homologous DNA is present, RecA-dependent homologous recombination is promoted; if not (or if no autologous Chi site is present) the RecBCD/AddAB complex eventually degrades the DNA. We examined the distribution of recBCD and addAB genes among bacteria, and sought ways to distinguish them unambiguously. We examined bacterial species among 33 phyla, finding some unexpected distribution patterns. RecBCD and addAB are less conserved than recA, with the orthologous recB and addA genes more conserved than the recC or addB genes. We were able to classify RecB vs. AddA and RecC vs. AddB in some bacteria where this had not previously been done. We used logo analysis to identify sequence segments that are conserved, but differ between the RecBC and AddAB proteins, to help future differentiation between members of these two families.

Role for first zinc finger of WT1 in DNA sequence specificity: Denys-Drash syndrome-associated WT1 mutant in ZF1 enhances affinity for a subset of WT1 binding sites.

Wilms tumor protein (WT1) is a Cys2-His2 zinc-finger transcription factor vital for embryonic development of the genitourinary system. The protein contains a C-terminal DNA binding domain with four tandem zinc-fingers (ZF1-4). An alternative splicing of Wt1 can add three additional amino acids-lysine (K), threonine (T) and serine (S)-between ZF3 and ZF4. In the -KTS isoform, ZF2-4 determine the sequence-specificity of DNA binding, whereas the function of ZF1 remains elusive. Three X-ray structures are described here for wild-type -KTS isoform ZF1-4 in complex with its cognate DNA sequence. We observed four unique ZF1 conformations. First, like ZF2-4, ZF1 can be positioned continuously in the DNA major groove forming a 'near-cognate' complex. Second, while ZF2-4 make base-specific interactions with one DNA molecule, ZF1 can interact with a second DNA molecule (or, presumably, two regions of the same DNA molecule). Third, ZF1 can intercalate at the joint of two tail-to-head DNA molecules. If such intercalation occurs on a continuous DNA molecule, it would kink the DNA at the ZF1 binding site. Fourth, two ZF1 units can dimerize. Furthermore, we examined a Denys-Drash syndrome-associated ZF1 mutation (methionine at position 342 is replaced by arginine). This mutation enhances WT1 affinity for a guanine base. X-ray crystallography of the mutant in complex with its preferred sequence revealed the interactions responsible for this affinity change. These results provide insight into the mechanisms of action of WT1, and clarify the fact that ZF1 plays a role in determining sequence specificity of this critical transcription factor.

Methyl-dependent and spatial-specific DNA recognition by the orthologous transcription factors human AP-1 and Epstein-Barr virus Zta.

T: Activator protein 1 (AP-1) is a transcription factor that recognizes two versions of a 7-base pair response element, either 5΄- GAG CA-3΄ or 5΄- GAG CA-3΄ (where M = 5-methylcytosine). These two elements share the feature that 5-methylcytosine and thymine both have a methyl group in the same position, 5-carbon of the pyrimidine, so each of them has two methyl groups at nucleotide positions 1 and 5 from the 5΄ end, resulting in four methyl groups symmetrically positioned in duplex DNA. Epstein-Barr Virus Zta is a key transcriptional regulator of the viral lytic cycle that is homologous to AP-1. Zta recognizes several methylated Zta-response elements, including meZRE1 (5΄- GAG C A-3΄) and meZRE2 (5΄- GAG G A-3΄), where a methylated cytosine occupies one of the inner thymine residues corresponding to the AP-1 element, resulting in the four spatially equivalent methyl groups. Here, we study how AP-1 and Zta recognize these methyl groups within their cognate response elements. These methyl groups are in van der Waals contact with a conserved di-alanine in AP-1 dimer (Ala265 and Ala266 in Jun), or with the corresponding Zta residues Ala185 and Ser186 (via its side chain carbon Cß atom). Furthermore, the two ZRE elements differ at base pair 6 (C:G versus G:C), forming a pseudo-symmetric sequence (meZRE1) or an asymmetric sequence (meZRE2). In vitro DNA binding assays suggest that Zta has high affinity for all four sequences examined, whereas AP-1 has considerably reduced affinity for the asymmetric sequence (meZRE2). We ascribe this difference to Zta Ser186 (a unique residue for Zta) whose side chain hydroxyl oxygen atom interacts with the two half sites differently, whereas the corresponding Ala266 of AP-1 Jun protein lacks such flexibility. Our analyses demonstrate a novel mechanism of 5mC/T recognition in a methylation-dependent, spatial and sequence-specific approach by basic leucine-zipper transcriptional factors.

Structural basis of human PR/SET domain 9 (PRDM9) allele C-specific recognition of its cognate DNA sequence.

PRDM9 is the only mammalian gene that has been associated with speciation. The PR/SET domain 9 (PRDM9) protein is a major determinant of meiotic recombination hot spots and acts through sequence-specific DNA binding via its C2H2 zinc finger (ZF) tandem array, which is highly polymorphic within and between species. The most common human variant, PRDM9 allele A (PRDM9a), contains 13 fingers (ZF1-13). Allele C (PRDM9c) is the second-most common among African populations and differs from PRDM9a by an arginine-to-serine change (R764S) in ZF9 and by replacement of ZF11 with two other fingers, yielding 14 fingers in PRDM9c. Here we co-crystallized the six-finger fragment ZF8-13 of PRDM9c, in complex with an oligonucleotide representing a known PRDM9c-specific hot spot sequence, and compared the structure with that of a characterized PRDM9a-specific complex. There are three major differences. First, Ser764 in ZF9 allows PRDM9c to accommodate a variable base, whereas PRDM9a Arg764 recognizes a conserved guanine. Second, the two-finger expansion of ZF11 allows PRDM9c to recognize three-base-pair-longer sequences. A tryptophan in the additional ZF interacts with a conserved thymine methyl group. Third, an Arg-Asp dipeptide immediately preceding the ZF helix, conserved in two PRDM9a fingers and three PRDM9c fingers, permits adaptability to variations from a C:G base pair (G-Arg interaction) to a G:C base pair (C-Asp interaction). This Arg-Asp conformational switch allows identical ZF modules to recognize different sequences. Our findings illuminate the molecular mechanisms for flexible and conserved binding of human PRDM9 alleles to their cognate DNA sequences.

Distinctive Klf4 mutants determine preference for DNA methylation status.

Reprogramming of mammalian genome methylation is critically important but poorly understood. Klf4, a transcription factor directing reprogramming, contains a DNA binding domain with three consecutive C2H2 zinc fingers. Klf4 recognizes CpG or TpG within a specific sequence. Mouse Klf4 DNA binding domain has roughly equal affinity for methylated CpG or TpG, and slightly lower affinity for unmodified CpG. The structural basis for this key preference is unclear, though the side chain of Glu446 is known to contact the methyl group of 5-methylcytosine (5mC) or thymine (5-methyluracil). We examined the role of Glu446 by mutagenesis. Substituting Glu446 with aspartate (E446D) resulted in preference for unmodified cytosine, due to decreased affinity for 5mC. In contrast, substituting Glu446 with proline (E446P) increased affinity for 5mC by two orders of magnitude. Structural analysis revealed hydrophobic interaction between the proline's aliphatic cyclic structure and the 5-methyl group of the pyrimidine (5mC or T). As in wild-type Klf4 (E446), the proline at position 446 does not interact directly with either the 5mC N4 nitrogen or the thymine O4 oxygen. In contrast, the unmethylated cytosine's exocyclic N4 amino group (NH2) and its ring carbon C5 atom hydrogen bond directly with the aspartate carboxylate of the E446D variant. Both of these interactions would provide a preference for cytosine over thymine, and the latter one could explain the E446D preference for unmethylated cytosine. Finally, we evaluated the ability of these Klf4 mutants to regulate transcription of methylated and unmethylated promoters in a luciferase reporter assay.