Search

It is difficult to exhaustively screen all possible DNA binding sequences for a given transcription factor (TF). Here, we developed the KaScape method, in which TFs bind to all possible DNA sequences in the same DNA pool where DNA sequences are prepared by randomized oligo synthesis and the random length can be adjusted to a length such as 4, 5, 6, or 7. After separating bound from unbound double-stranded DNAs (dsDNAs), their sequences are determined by next-generation sequencing. To demonstrate the relative binding affinities of all possible DNA sequences determined by KaScape, we developed three-dimensional KaScape viewing software based on a K-mer graph. We applied KaScape to 12 plant TF family AtWRKY proteins and found that all AtWRKY proteins bound to the core sequence GAC with similar profiles. KaScape can detect not only binding sequences consistent with the consensus W-box "TTGAC(C/T)" but also other sequences with weak affinity. KaScape provides a high-throughput, easy-to-operate, sensitive, and exhaustive method for quantitatively characterizing the relative binding strength of a TF with all possible binding sequences, allowing us to comprehensively characterize the specificity and affinity landscape of transcription factors, particularly for moderate- and low-affinity binding sites., (© 2023. Springer Nature Limited.)

Plant-specific TCP transcription factors are key regulators of diverse plant functions. TCP transcription factors have long been annotated as basic helix-loop-helix (bHLH) transcription factors according to remote sequence homology without experimental validation, and their consensus DNA-binding sequences and protein-DNA recognition mechanisms have remained elusive. Here, we report the crystal structures of the class I TCP domain from AtTCP15 and the class II TCP domain from AtTCP10 in complex with different double-stranded DNA (dsDNA). The complex structures reveal that the TCP domain is a distinct DNA-binding motif and the homodimeric TCP domains adopt a unique three-site recognition mode, binding to dsDNA mainly through a central pair of β-strands formed by the dimer interface and two basic flexible loops from each monomer. The consensus DNA-binding sequence for class I TCPs is a perfectly palindromic 11 bp (GTGGGNCCCAC), whereas that for class II TCPs is a near-palindromic 11 bp (GTGGTCCCCAC). The unique DNA binding mode allows the TCP domains to display broad specificity for a range of DNA sequences even shorter than 11 bp, adding further complexity to the regulatory network of plant TCP transcription factors., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein-DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

DNA methylation on CpG sites is the most common epigenetic modification. Recently, methylation in a non-CpG context was found to occur widely on genomic DNA. Moreover, methylation of non-CpG sites is a highly controlled process, and its level may vary during cellular development. To study non-CpG methylation effects on DNA/protein interactions, we have chosen three human transcription factors (TFs): glucocorticoid receptor (GR), brain and muscle ARNT-like 1 (BMAL1) - circadian locomotor output cycles kaput (CLOCK) and estrogen receptor (ER) with methylated or unmethylated DNA binding sequences, using single-molecule and isothermal titration calorimetry assays. The results demonstrated that these TFs interact with methylated DNA with different effects compared with their cognate DNA sequences. The effects of non-CpG methylation on transcriptional regulation were validated by cell-based luciferase assay at protein level. The mechanisms of non-CpG methylation influencing DNA-protein interactions were investigated by crystallographic analyses and molecular dynamics simulation. With BisChIP-seq assays in HEK-293T cells, we found that GR can recognize highly methylated sites within chromatin in cells. Therefore, we conclude that non-CpG methylation of DNA can provide a mechanism for regulating gene expression through directly affecting the binding of TFs.

Allostery is well documented for proteins but less recognized for DNA-protein interactions. Here, we report that specific binding of a protein on DNA is substantially stabilized or destabilized by another protein bound nearby. The ternary complex's free energy oscillates as a function of the separation between the two proteins with a periodicity of ~10 base pairs, the helical pitch of B-form DNA, and a decay length of ~15 base pairs. The binding affinity of a protein near a DNA hairpin is similarly dependent on their separation, which-together with molecular dynamics simulations-suggests that deformation of the double-helical structure is the origin of DNA allostery. The physiological relevance of this phenomenon is illustrated by its effect on gene expression in live bacteria and on a transcription factor's affinity near nucleosomes.

The smu.134 gene encodes a putative transcriptional regulator of 217 residues in Streptococcus mutans, a major pathogen for human dental caries. The gene was cloned into expression vector pET28alpha and expressed in soluble form in E. coli strain BL21 (DE3) with a His tag at its N-terminus. The recombinant protein SMU.134 was purified to homogeneity in a two step procedure of Ni(2+ ) chelating and size exclusion chromatography. Crystals suitable for X-ray diffraction were obtained by hanging-drop vapor diffusion method and diffracted to 2.6 A;. The crystal belonged to space group P2(1)2(1)2(1), with unit-cell parameters a=55.03 A, b=80.84 A, c=107.96 A.

WRKY proteins, defined by the conserved WRKYGQK sequence, are comprised of a large superfamily of transcription factors identified specifically from the plant kingdom. This superfamily plays important roles in plant disease resistance, abiotic stress, senescence as well as in some developmental processes. In this study, the Arabidopsis WRKY1 was shown to be involved in the salicylic acid signaling pathway and partially dependent on NPR1; a C-terminal domain of WRKY1, AtWRKY1-C, was constructed for structural studies. Previous investigations showed that DNA binding of the WRKY proteins was localized at the WRKY domains and these domains may define novel zinc-binding motifs. The crystal structure of the AtWRKY1-C determined at 1.6 A resolution has revealed that this domain is composed of a globular structure with five beta strands, forming an antiparallel beta-sheet. A novel zinc-binding site is situated at one end of the beta-sheet, between strands beta4 and beta5. Based on this high-resolution crystal structure and site-directed mutagenesis, we have defined and confirmed that the DNA-binding residues of AtWRKY1-C are located at beta2 and beta3 strands. These results provided us with structural information to understand the mechanism of transcriptional control and signal transduction events of the WRKY proteins.

The C-terminal WRKY domain of Arabidopsis thaliana WRKY1 protein, a transcription factor, was cloned and expressed. The expressed protein was then purified and crystallized. The preliminary X-ray analysis was undertaken. The crystal diffracted to 2.50 A resolution in-house and belongs to space group P2(1) with unit-cell parameters a=64.10 A, b=34.88 A, c=114.72 A, beta=90.49 degrees .

CLOCK (circadian locomotor output cycles kaput) and BMAL1 (brain and muscle ARNT-like 1) are both transcription factors of the circadian core loop in mammals. Recently published mouse CLOCK-BMAL1 bHLH (basic helix-loop-helix)-PAS (period-ARNT-single-minded) complex structure sheds light on the mechanism for heterodimer formation, but the structural details of the protein-DNA recognition mechanisms remain elusive. Here we have elucidated the crystal structure of human CLOCK-BMAL1 bHLH domains bound to a canonical E-box DNA. We demonstrate that CLOCK and BMAL1 bHLH domains can be mutually selected, and that hydrogen-bonding networks mediate their E-box recognition. We identified a hydrophobic contact between BMAL1 Ile80 and a flanking thymine nucleotide, suggesting that CLOCK-BMAL1 actually reads 7-bp DNA and not the previously believed 6-bp DNA. To find potential non-canonical E-boxes that could be recognized by CLOCK-BMAL1, we constructed systematic single-nucleotide mutations on the E-box and measured their relevant affinities. We defined two non-canonical E-box patterns with high affinities, AACGTGA and CATGTGA, in which the flanking A7-T7′ base pair is indispensable for recognition. These results will help us to identify functional CLOCK-BMAL1-binding sites in vivo and to search for clock-controlled genes. Furthermore, we assessed the inhibitory role of potential phosphorylation sites in bHLH regions. We found that the phospho-mimicking mutation on BMAL1 Ser78 could efficiently block DNA binding as well as abolish normal circadian oscillation in cells. We propose that BMAL1 Ser78 should be a key residue mediating input signal-regulated transcriptional inhibition for external cues to entrain the circadian clock by kinase cascade. [ABSTRACT FROM AUTHOR]