Search

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The modification of DNA bases is a classic hallmark of epigenetics. Four forms of modified cytosine-5-methylcytosine, 5-hydroxymethylcytosine, 5-formylcytosine, and 5-carboxylcytosine-have been discovered in eukaryotic DNA. In addition to cytosine carbon-5 modifications, cytosine and adenine methylated in the exocyclic amine-N4-methylcytosine and N6-methyladenine-are other modified DNA bases discovered even earlier. Each modified base can be considered a distinct epigenetic signal with broader biological implications beyond simple chemical changes. Since 1994, several crystal structures of proteins and enzymes involved in writing, reading, and erasing modified bases have become available. Here, we present a structural synopsis of writers, readers, and erasers of the modified bases from prokaryotes and eukaryotes. Despite significant differences in structures and functions, they are remarkably similar regarding their engagement in flipping a target base/nucleotide within DNA for specific recognitions and/or reactions. We thus highlight base flipping as a common structural framework broadly applied by distinct classes of proteins and enzymes across phyla for epigenetic regulations of DNA., (© 2022. The Author(s), under exclusive license to Springer Nature Switzerland AG.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Metazoan transcription factors typically regulate large numbers of genes. Here we identify via a CRISPR-Cas9 genetic screen ZNF410, a pentadactyl DNA-binding protein that in human erythroid cells directly activates only a single gene, the NuRD component CHD4. Specificity is conveyed by two highly evolutionarily conserved clusters of ZNF410 binding sites near the CHD4 gene with no counterparts elsewhere in the genome. Loss of ZNF410 in adult-type human erythroid cell culture systems and xenotransplantation settings diminishes CHD4 levels and derepresses the fetal hemoglobin genes. While previously known to be silenced by CHD4, the fetal globin genes are exposed here as among the most sensitive to reduced CHD4 levels.. In vitro DNA binding assays and crystallographic studies reveal the ZNF410-DNA binding mode. ZNF410 is a remarkably selective transcriptional activator in erythroid cells, and its perturbation might offer new opportunities for treatment of hemoglobinopathies., Competing Interests: Declaration of Interests X.J., J.S., and G.A.B. are contributors to a patent filed on behalf of the Children’s Hospital of Philadelphia., (Copyright © 2020 Elsevier Inc. All rights reserved.)

HhaI, a Type II restriction endonuclease, recognizes the symmetric sequence 5'-GCG↓C-3' in duplex DNA and cleaves ('↓') to produce fragments with 2-base, 3'-overhangs. We determined the structure of HhaI in complex with cognate DNA at an ultra-high atomic resolution of 1.0 Å. Most restriction enzymes act as dimers with two catalytic sites, and cleave the two strands of duplex DNA simultaneously, in a single binding event. HhaI, in contrast, acts as a monomer with only one catalytic site, and cleaves the DNA strands sequentially, one after the other. HhaI comprises three domains, each consisting of a mixed five-stranded β sheet with a defined function. The first domain contains the catalytic-site; the second contains residues for sequence recognition; and the third contributes to non-specific DNA binding. The active-site belongs to the 'PD-D/EXK' superfamily of nucleases and contains the motif SD-X11-EAK. The first two domains are similar in structure to two other monomeric restriction enzymes, HinP1I (G↓CGC) and MspI (C↓CGG), which produce fragments with 5'-overhangs. The third domain, present only in HhaI, shifts the positions of the recognition residues relative to the catalytic site enabling this enzyme to cleave the recognition sequence at a different position. The structure of M.HhaI, the biological methyltransferase partner of HhaI, was determined earlier. Together, these two structures represent the first natural pair of restriction-modification enzymes to be characterized in atomic detail., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2018. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Tandem zinc finger (ZF) proteins are the largest and most rapidly diverging family of DNA-binding transcription regulators in mammals. ZFP568 represses a transcript of placental-specific insulin like growth factor 2 (Igf2-P0) in mice. ZFP568 binds a 24-base pair sequence-specific element upstream of Igf2-P0 via the eleven-ZF array. Both DNA and protein conformations deviate from the conventional one finger-three bases recognition, with individual ZFs contacting 2, 3, or 4 bases and recognizing thymine on the opposite strand. These interactions arise from a shortened minor groove caused by an AT-rich stretch, suggesting adaptability of ZF arrays to sequence variations. Despite conservation in mammals, mutations at Igf2 and ZFP568 reduce their binding affinity in chimpanzee and humans. Our studies provide important insights into the evolutionary and structural dynamics of ZF-DNA interactions that play a key role in mammalian development and evolution., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The multidomain CCCTC-binding factor (CTCF), containing a tandem array of 11 zinc fingers (ZFs), modulates the three-dimensional organization of chromatin. We crystallized the human CTCF DNA-binding domain in complex with a known CTCF-binding site. While ZF2 does not make sequence-specific contacts, each finger of ZF3-7 contacts three bases of the 15-bp consensus sequence. Each conserved nucleotide makes base-specific hydrogen bonds with a particular residue. Most of the variable base pairs within the core sequence also engage in interactions with the protein. These interactions compensate for deviations from the consensus sequence, allowing CTCF to adapt to sequence variations. CTCF is sensitive to cytosine methylation at position 2, but insensitive at position 12 of the 15-bp core sequence. These differences can be rationalized structurally. Although included in crystallizations, ZF10 and ZF11 are not visible, while ZF8 and ZF9 span the backbone of the DNA duplex, conferring no sequence specificity but adding to overall binding stability., (Copyright © 2017 Elsevier Inc. All rights reserved.)

DNA binding specificities of transcription factors (TFs) are a key component of gene regulatory processes. Underlying mechanisms that explain the highly specific binding of TFs to their genomic target sites are poorly understood. A better understanding of TF-DNA binding requires the ability to quantitatively model TF binding to accessible DNA as its basic step, before additional in vivo components can be considered. Traditionally, these models were built based on nucleotide sequence. Here, we integrated 3D DNA shape information derived with a high-throughput approach into the modeling of TF binding specificities. Using support vector regression, we trained quantitative models of TF binding specificity based on protein binding microarray (PBM) data for 68 mammalian TFs. The evaluation of our models included cross-validation on specific PBM array designs, testing across different PBM array designs, and using PBM-trained models to predict relative binding affinities derived from in vitro selection combined with deep sequencing (SELEX-seq). Our results showed that shape-augmented models compared favorably to sequence-based models. Although both k-mer and DNA shape features can encode interdependencies between nucleotide positions of the binding site, using DNA shape features reduced the dimensionality of the feature space. In addition, analyzing the feature weights of DNA shape-augmented models uncovered TF family-specific structural readout mechanisms that were not revealed by the DNA sequence. As such, this work combines knowledge from structural biology and genomics, and suggests a new path toward understanding TF binding and genome function.

Until now, it has been reasonably assumed that specific base-pair recognition is the only mechanism controlling the specificity of transcription factor (TF)-DNA binding. Contrary to this assumption, here we show that nonspecific DNA sequences possessing certain repeat symmetries, when present outside of specific TF binding sites (TFBSs), statistically control TF-DNA binding preferences. We used high-throughput protein-DNA binding assays to measure the binding levels and free energies of binding for several human TFs to tens of thousands of short DNA sequences with varying repeat symmetries. Based on statistical mechanics modeling, we identify a new protein-DNA binding mechanism induced by DNA sequence symmetry in the absence of specific base-pair recognition, and experimentally demonstrate that this mechanism indeed governs protein-DNA binding preferences.

MspJI belongs to a family of restriction enzymes that cleave DNA containing 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC). MspJI is specific for the sequence 5(h)mC-N-N-G or A and cleaves with some variability 9/13 nucleotides downstream. Earlier, we reported the crystal structure of MspJI without DNA and proposed how it might recognize this sequence and catalyze cleavage. Here we report its co-crystal structure with a 27-base pair oligonucleotide containing 5mC. This structure confirms that MspJI acts as a homotetramer and that the modified cytosine is flipped from the DNA helix into an SRA-like-binding pocket. We expected the structure to reveal two DNA molecules bound specifically to the tetramer and engaged with the enzyme's two DNA-cleavage sites. A coincidence of crystal packing precluded this organization, however. We found that each DNA molecule interacted with two adjacent tetramers, binding one specifically and the other non-specifically. The latter interaction, which prevented cleavage-site engagement, also involved base flipping and might represent the sequence-interrogation phase that precedes specific recognition. MspJI is unusual in that DNA molecules are recognized and cleaved by different subunits. Such interchange of function might explain how other complex multimeric restriction enzymes act., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

AbaSI, a member of the PvuRts1I-family of modification-dependent restriction endonucleases, cleaves deoxyribonucleic acid (DNA) containing 5-hydroxymethylctosine (5hmC) and glucosylated 5hmC (g5hmC), but not DNA containing unmodified cytosine. AbaSI has been used as a tool for mapping the genomic locations of 5hmC, an important epigenetic modification in the DNA of higher organisms. Here we report the crystal structures of AbaSI in the presence and absence of DNA. These structures provide considerable, although incomplete, insight into how this enzyme acts. AbaSI appears to be mainly a homodimer in solution, but interacts with DNA in our structures as a homotetramer. Each AbaSI subunit comprises an N-terminal, Vsr-like, cleavage domain containing a single catalytic site, and a C-terminal, SRA-like, 5hmC-binding domain. Two N-terminal helices mediate most of the homodimer interface. Dimerization brings together the two catalytic sites required for double-strand cleavage, and separates the 5hmC binding-domains by ∼70 Å, consistent with the known activity of AbaSI which cleaves DNA optimally between symmetrically modified cytosines ∼22 bp apart. The eukaryotic SET and RING-associated (SRA) domains bind to DNA containing 5-methylcytosine (5mC) in the hemi-methylated CpG sequence. They make contacts in both the major and minor DNA grooves, and flip the modified cytosine out of the helix into a conserved binding pocket. In contrast, the SRA-like domain of AbaSI, which has no sequence specificity, contacts only the minor DNA groove, and in our current structures the 5hmC remains intra-helical. A conserved, binding pocket is nevertheless present in this domain, suitable for accommodating 5hmC and g5hmC. We consider it likely, therefore, that base-flipping is part of the recognition and cleavage mechanism of AbaSI, but that our structures represent an earlier, pre-flipped stage, prior to actual recognition., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Motivation: The DNA binding specificity of a transcription factor (TF) is typically represented using a position weight matrix model, which implicitly assumes that individual bases in a TF binding site contribute independently to the binding affinity, an assumption that does not always hold. For this reason, more complex models of binding specificity have been developed. However, these models have their own caveats: they typically have a large number of parameters, which makes them hard to learn and interpret., Results: We propose novel regression-based models of TF-DNA binding specificity, trained using high resolution in vitro data from custom protein-binding microarray (PBM) experiments. Our PBMs are specifically designed to cover a large number of putative DNA binding sites for the TFs of interest (yeast TFs Cbf1 and Tye7, and human TFs c-Myc, Max and Mad2) in their native genomic context. These high-throughput quantitative data are well suited for training complex models that take into account not only independent contributions from individual bases, but also contributions from di- and trinucleotides at various positions within or near the binding sites. To ensure that our models remain interpretable, we use feature selection to identify a small number of sequence features that accurately predict TF-DNA binding specificity. To further illustrate the accuracy of our regression models, we show that even in the case of paralogous TF with highly similar position weight matrices, our new models can distinguish the specificities of individual factors. Thus, our work represents an important step toward better sequence-based models of individual TF-DNA binding specificity., Availability: Our code is available at http://genome.duke.edu/labs/gordan/ISMB2013. The PBM data used in this article are available in the Gene Expression Omnibus under accession number GSE47026.

DNA sequence is a major determinant of the binding specificity of transcription factors (TFs) for their genomic targets. However, eukaryotic cells often express, at the same time, TFs with highly similar DNA binding motifs but distinct in vivo targets. Currently, it is not well understood how TFs with seemingly identical DNA motifs achieve unique specificities in vivo. Here, we used custom protein-binding microarrays to analyze TF specificity for putative binding sites in their genomic sequence context. Using yeast TFs Cbf1 and Tye7 as our case studies, we found that binding sites of these bHLH TFs (i.e., E-boxes) are bound differently in vitro and in vivo, depending on their genomic context. Computational analyses suggest that nucleotides outside E-box binding sites contribute to specificity by influencing the three-dimensional structure of DNA binding sites. Thus, the local shape of target sites might play a widespread role in achieving regulatory specificity within TF families., (Copyright © 2013 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: We hypothesized that single-nucleotide polymorphisms (SNPs) associated with heightened in vitro platelet function during aspirin exposure (which we define as "laboratory aspirin resistance") would be associated with greater risk for death, myocardial infarction (MI) or stroke among patients with coronary artery disease regularly using aspirin., Methods: Duke Databank for Cardiovascular Disease patients with (n = 3,449, CATHeterization GENetics cohort) or without (n = 11,754, nongenetic cohort) banked DNA with ≥1 coronary stenosis >75% were followed up at 6 months, then annually for death, MI, or stroke occurring during periods of reported aspirin use. We evaluated associations of candidate SNPs from GNB3, PEAR1, ITGB3, VAV3, ITGA2, GPVI, PTGS1, F2R, THBS1, A2AR, and GP1BA with events during follow-up using Cox proportional hazards modeling adjusted for clinical characteristics associated with outcomes in the nongenetic cohort., Results: Over a median of 3.5 years, 2,762 (24%) nongenetic cohort patients and 648 (19%) CATHeterization GENetics cohort patients had the composite outcome during reported aspirin use. No candidate SNPs were significantly associated with death, MI, or stroke in either univariable or multivariable analyses. A prospective analysis demonstrated 80% to 88% power to detect a hazard ratio of ≥1.3 for minor allele carriers., Conclusion: Patients with angiographically significant coronary artery disease regularly using aspirin and carrying SNPs associated with laboratory aspirin resistance were not at higher risk for death, MI, or stroke. Using these SNPs to guide more aggressive antiplatelet therapy is not justified by these results. Direct extrapolation from in vitro findings to the clinical setting should be avoided., (Copyright © 2011 Mosby, Inc. All rights reserved.)

Maintenance methylation of hemimethylated CpG dinucleotides at DNA replication forks is the key to faithful mitotic inheritance of genomic methylation patterns. UHRF1 (ubiquitin-like, containing PHD and RING finger domains 1) is required for maintenance methylation by interacting with DNA nucleotide methyltransferase 1 (DNMT1), the maintenance methyltransferase, and with hemimethylated CpG, the substrate for DNMT1 (refs 1 and 2). Here we present the crystal structure of the SET and RING-associated (SRA) domain of mouse UHRF1 in complex with DNA containing a hemimethylated CpG site. The DNA is contacted in both the major and minor grooves by two loops that penetrate into the middle of the DNA helix. The 5-methylcytosine has flipped completely out of the DNA helix and is positioned in a binding pocket with planar stacking contacts, Watson-Crick polar hydrogen bonds and van der Waals interactions specific for 5-methylcytosine. Hence, UHRF1 contains a previously unknown DNA-binding module and is the first example of a non-enzymatic, sequence-specific DNA-binding protein domain to use the base flipping mechanism to interact with DNA.

HinP1I recognizes and cleaves the palindromic tetranucleotide sequence G downward arrowCGC in DNA. We report three structures of HinP1I-DNA complexes: in the presence of Ca(2+) (pre-reactive complex), in the absence of metal ion (binary complex) and in the presence of Mg(2+) (post-reactive complex). HinP1I forms a back-to-back dimer with two active sites and two DNA duplexes bound on the outer surfaces of the dimer facing away from each other. The 10 bp DNA duplexes undergo protein-induced distortions exhibiting features of A-, B- and Z-conformations: bending on one side (by intercalation of a phenylalanine side chain into the major groove), base flipping on the other side of the recognition site (by expanding the step rise distance of the local base pair to Z-form) and a local A-form conformation between the two central C:G base pairs of the recognition site (by binding of the N-terminal helix in the minor groove). In the pre- and post-reactive complexes, two metals (Ca(2+) or Mg(2+)) are found in the active site. The enzyme appears to cleave DNA sequentially, hydrolyzing first one DNA strand, as seen in the post-reactive complex in the crystalline state, and then the other, as supported by the observation that, in solution, a nicked DNA intermediate accumulates before linearization.

DNA methyltransferases methylate target bases within specific nucleotide sequences. Three structures are described for bacteriophage T4 DNA-adenine methyltransferase (T4Dam) in ternary complexes with partially and fully specific DNA and a methyl-donor analog. We also report the effects of substitutions in the related Escherichia coli DNA methyltransferase (EcoDam), altering residues corresponding to those involved in specific interaction with the canonical GATC target sequence in T4Dam. We have identified two types of protein-DNA interactions: discriminatory contacts, which stabilize the transition state and accelerate methylation of the cognate site, and antidiscriminatory contacts, which do not significantly affect methylation of the cognate site but disfavor activity at noncognate sites. These structures illustrate the transition in enzyme-DNA interaction from nonspecific to specific interaction, suggesting that there is a temporal order for formation of specific contacts.

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket ('base flipping') is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI-DNA-AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.

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 CpGflanking 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. [ABSTRACT FROM AUTHOR]

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 30 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. [ABSTRACT FROM AUTHOR]

Epigenetically targeted therapeutic development, particularly for SAM-dependent methylations of DNA, mRNA and histones has been proceeding rapidly for cancer treatments over the past few years. However, this approach has barely begun to be exploited for developing new antibiotics, despite an overwhelming global need to counter antimicrobial resistance. Here, we explore whether SAM analogues, some of which are in (pre)clinical studies as inhibitors of human epigenetic enzymes, can also inhibit Clostridioides difficile-specific DNA adenine methyltransferase (CamA), a sporulation regulator present in all C. difficile genomes sequenced to date, but found in almost no other bacteria. We found that SGC0946 (an inhibitor of DOT1L), JNJ-64619178 (an inhibitor of PRMT5) and SGC8158 (an inhibitor of PRMT7) inhibit CamA enzymatic activity in vitro at low micromolar concentrations. Structural investigation of the ternary complexes of CamA-DNA in the presence of SGC0946 or SGC8158 revealed conformational rearrangements of the N-terminal arm, with no apparent disturbance of the active site. This N-terminal arm and its modulation of exchanges between SAM (the methyl donor) and SAH (the reaction product) during catalysis of methyl transfer are, to date, unique to CamA. Our work presents a substantial first step in generating potent and selective inhibitors of CamA that would serve in the near term as chemical probes to investigate the cellular mechanism(s) of CamA in controlling spore formation and colonization, and eventually as therapeutic antivirulence agents useful in treating C. difficile infection. [ABSTRACT FROM AUTHOR]

DNA binding specificities of transcription factors (TFs) are a key component of gene regulatory processes. Underlying mechanisms that explain the highly specific binding of TFs to their genomic target sites are poorly understood. A better understanding of TF-DNA binding requires the ability to quantitatively model TF binding to accessible DNA as its basic step, before additional in vivo components can be considered. Traditionally, these models were built based on nucleotide sequence. Here, we integrated 3D DNA shape information derived with a high-throughput approach into the modeling of TF binding specificities. Using support vector regression, we trained quantitative models of TF binding specificity based on protein binding microarray (PBM) data for 68 mammalian TFs. The evaluation of our models included cross-validation on specific PBM array designs, testing across different PBM array designs, and using PBM-trained models to predict relative binding affinities derived from in vitro selection combined with deep sequencing (SELEX-seq). Our results showed that shape-augmented models compared favorably to sequence-based models. Although both k-mer and DNA shape features can encode interdependencies between nucleotide positions of the binding site, using DNA shape features reduced the dimensionality of the feature space. In addition, analyzing the feature weights of DNA shape-augmented models uncovered TF family-specific structural readout mechanisms that were not revealed by the DNA sequence. As such, this work combines knowledge from structural biology and genomics, and suggests a new path toward understanding TF binding and genome function. [ABSTRACT FROM AUTHOR]

Abstract: M-phase phosphoprotein 8 (MPP8) harbors an N-terminal chromodomain and a C-terminal ankyrin repeat domain. MPP8, via its chromodomain, binds histone H3 peptide tri- or di-methylated at lysine 9 (H3K9me3/H3K9me2) in submicromolar affinity. We determined the crystal structure of MPP8 chromodomain in complex with H3K9me3 peptide. MPP8 interacts with at least six histone H3 residues from glutamine 5 to serine 10, enabling its ability to distinguish lysine-9-containing peptide (QTARKS) from that of lysine 27 (KAARKS), both sharing the ARKS sequence. A partial hydrophobic cage with three aromatic residues (Phe59, Trp80 and Tyr83) and one aspartate (Asp87) encloses the methylated lysine 9. MPP8 has been reported to be phosphorylated in vivo, including the cage residue Tyr83 and the succeeding Thr84 and Ser85. Modeling a phosphate group onto the side-chain hydroxyl oxygen of Tyr83 suggests that the negatively charged phosphate group could enhance the binding of positively charged methyl-lysine or create a regulatory signal by allowing or inhibiting binding of other protein(s). [Copyright &y& Elsevier]

The crystal structure of the Escherichia coli DNA adenine methyltransferase (EcoDam) in a binary complex with the cofactor product S-adenosyl-L-homocysteine (AdoHcy) unexpectedly showed the bound AdoHcy in two alternative conformations, extended or folded. The extended conformation represents the catalytically competent conformation, identical to that of EcoDam-DNA-AdoHcy ternary complex. The folded conformation prevents catalysis, because the homocysteine moiety occupies the target Ade binding pocket. The largest difference between the binary and ternary structures is in the conformation of the N-terminal hexapeptide (9KWAGGK14), Cofactor binding leads to a strong change in the fluorescence of Trp10, whose indole ring approaches the cofactor by 3.3 Å. Stopped-flow kinetics and AdoMet cross-linking studies indicate that the cofactor prefers binding to the enzyme after preincubation with DNA. In the presence of DNA, AdoMet binding is ~2-fold stronger than AdoHcy binding. In the binary complex the side chain of Lys14 is disordered, whereas Lys14 stabilizes the active site in the ternary complex. Fluorescence stopped-flow experiments indicate that Lys14 is important for EcoDam binding of the extrahelical target base into the active site pocket. This suggests that the hexapeptide couples specific DNA binding (Lys9), AdoMet binding (Trp10), and insertion of the flipped target base into the active site pocket (Lys14). [ABSTRACT FROM AUTHOR]

The structure of the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase in complex with cognate DNA was determined at 1.89Å resolution in the presence of S-adenosyl-l-homocysteine. DNA recognition and the dynamics of base-flipping were studied by site-directed mutagenesis, DNA methylation kinetics and fluorescence stopped-flow experiments. Our data illustrate the mechanism of coupling of DNA recognition and base-flipping. Contacts to the non-target strand in the second (3′) half of the GATC site are established by R124 to the fourth base-pair, and by L122 and P134 to the third base-pair. The aromatic ring of Y119 intercalates into the DNA between the second and third base-pairs, which is essential for base-flipping to occur. Compared to previous published structures of bacteriophage T4 Dam, three major new observations are made in E.coli Dam. (1) The first Gua is recognized by K9, removal of which abrogates the first base-pair recognition. (2) The flipped target Ade binds to the surface of EcoDam in the absence of S-adenosyl-l-methionine, which illustrates a possible intermediate in the base-flipping pathway. (3) The orphaned Thy residue displays structural flexibility by adopting an extrahelical or intrahelical position where it is in contact to N120. [Copyright &y& Elsevier]

DNA-adenine methylation at certain GATC sites plays a pivotal role in bacterial and phage gene expression as well as bacterial virulence. We report here the crystal structures of the bacteriophage T4Dam DNA adenine methyltransferase (MTase) in a binary complex with the methyl-donor product S-adenosyl-L-homocysteine (AdoHcy) and in a ternary complex with a synthetic 12-bp DNA duplex and AdoHcy. T4Dam contains two domains: a seven-stranded catalytic domain that harbors the binding site for AdoHcy and a DNA binding domain consisting of a five-helix bundle and a ß-hairpin that is conserved in the family of GATC-related MTase orthologs. Unexpectedly, the sequence-specific T4Dam bound to DNA in a nonspecific mode that contained two Dam monomers per synthetic duplex, even though the DNA contains a single GATC site. The ternary structure provides a rare snapshot of an enzyme poised for linear diffusion along the DNA. [ABSTRACT FROM AUTHOR]

Three structures have been determined for complexes between HhaI methyltransferase (M.HhaI) and oligonucleotides containing a G:A, G:U or G:AP (AP = abasic or apurinic/apyrimidinic) mismatch at the target base pair. The mismatched adenine, uracil and abasic site are all flipped out of the DNA helix and located in the enzyme's active-site pocket, adopting the same conformation as in the flipped-out normal substrate. These results, particularly the flipped-out abasic deoxyribose sugar, provide insight into the mechanism of base flipping. If the process involves the protein pushing the base out of the helix, then the push must take place not on the base, but rather on the sugar-phosphate backbone. Thus rotation of the DNA backbone is probably the key to base flipping. [ABSTRACT FROM AUTHOR]

The Caulobacter crescentus cell cycle-regulated DNA methyltransferase (CcrM) methylates the adenine of hemimethylated GANTC after replication. Here we present the structure of CcrM in complex with double-stranded DNA containing the recognition sequence. CcrM contains an N-terminal methyltransferase domain and a C-terminal nonspecific DNA-binding domain. CcrM is a dimer, with each monomer contacting primarily one DNA strand: the methyltransferase domain of one molecule binds the target strand, recognizes the target sequence, and catalyzes methyl transfer, while the C-terminal domain of the second molecule binds the non-target strand. The DNA contacts at the 5-base pair recognition site results in dramatic DNA distortions including bending, unwinding and base flipping. The two DNA strands are pulled apart, creating a bubble comprising four recognized base pairs. The five bases of the target strand are recognized meticulously by stacking contacts, van der Waals interactions and specific Watson–Crick polar hydrogen bonds to ensure high enzymatic specificity. CcrM is a cell cycle-regulated DNA methyltransferase that methylates an adenine within a specific sequence following replication in the gram negative bacterium Caulobacter crescentus. Here the authors present a crystal structure of DNA-bound CcrM that reveals the molecular mechanism leading to sequence-specific methylation. [ABSTRACT FROM AUTHOR]

DIM-5 is a SUV39-type histone H3 Lys9 methyltransferase that is essential for DNA methylation in N. crassa. We report the structure of a ternary complex including DIM-5, S-adenosyl-L-homocysteine, and a substrate H3 peptide. The histone tail inserts as a parallel strand between two DIM-5 strands, completing a hybrid sheet. Three post-SET cysteines coordinate a zinc atom together with Cys242 from the SET signature motif (NHXCXPN) near the active site. Consequently, a narrow channel is formed to accommodate the target Lys9 side chain. The sulfur atom of S-adenosyl-L-homocysteine, where the transferable methyl group is to be attached in S-adenosyl-L-methionine, lies at the opposite end of the channel, ∼4 A˚ away from the target Lys9 nitrogen. Structural comparison of the active sites of DIM-5, an H3 Lys9 trimethyltransferase, and SET7/9, an H3 Lys4 monomethyltransferase, allowed us to design substitutions in both enzymes that profoundly alter their product specificities without affecting their catalytic activities. [Copyright &y& Elsevier]