Sequence-specific dynamics of DNA response elements and their flanking sites regulate the recognition by AP-1 transcription factors

Updated understanding of the protein-DNA recognition code used by C2H2 zinc finger proteins

C2H2 zinc-finger (ZF) proteins form the largest family of DNA-binding transcription factors coded by mammalian genomes. In a typical DNA-binding ZF module, there are twelve residues (numbered from -1 to -12) between the last zinc-coordinating cysteine and the first zinc-coordinating histidine. The established C2H2-ZF "recognition code" suggests that residues at positions -1, -4, and -7 recognize the 5', central, and 3' bases of a DNA base-pair triplet, respectively. Structural studies have highlighted that additional residues at positions -5 and -8 also play roles in specific DNA recognition. The presence of bulky and either charged or polar residues at these five positions determines specificity for given DNA bases: guanine is recognized by arginine, lysine, or histidine; adenine by asparagine or glutamine; thymine or 5-methylcytosine by glutamate; and unmodified cytosine by aspartate. This review discusses recent structural characterizations of C2H2-ZFs that add to our understanding of the principles underlying the C2H2-ZF recognition code.

C2H2 Zinc Finger Transcription Factors Associated with Hemoglobinopathies

In humans, specific aberrations in β-globin results in sickle cell disease and β-thalassemia, symptoms of which can be ameliorated by increased expression of fetal globin (HbF). Two recent CRISPR-Cas9 screens, centered on ∼1500 annotated sequence-specific DNA binding proteins and performed in a human erythroid cell line that expresses adult hemoglobin, uncovered four groups of candidate regulators of HbF gene expression. They are (1) members of the nucleosome remodeling and deacetylase (NuRD) complex proteins that are already known for HbF control; (2) seven C2H2 zinc finger (ZF) proteins, including some (ZBTB7A and BCL11A) already known for directly silencing the fetal γ-globin genes in adult human erythroid cells; (3) a few other transcription factors of different structural classes that might indirectly influence HbF gene expression; and (4) DNA methyltransferase 1 (DNMT1) that maintains the DNA methylation marks that attract the MBD2-associated NuRD complex to DNA as well as associated histone H3 lysine 9 methylation. Here we briefly discuss the effects of these regulators, particularly C2H2 ZFs, in inducing HbF expression for treating β-hemoglobin disorders, together with recent advances in developing safe and effective small-molecule therapeutics for the regulation of this well-conserved hemoglobin switch.

Structural and dynamical aspect of DNA motif sequence specific binding of AP-1 transcription factor

Activator protein-1 (AP-1) comprises one of the largest and most evolutionary conserved families of ubiquitous eukaryotic transcription factors that act as a pioneer factor. Diversity in DNA binding interaction of AP-1 through a conserved basic-zipper (bZIP) domain directs in-depth understanding of how AP-1 achieves its DNA binding selectivity and consequently gene regulation specificity. Here, we address the structural and dynamical aspects of the DNA target recognition process of AP-1 using microsecond-long atomistic simulations based on the structure of the human AP-1 FosB/JunD bZIP-DNA complex. Our results show the unique role of DNA shape features in selective base specific interactions, characteristic ion population, and solvation properties of DNA grooves to form the motif sequence specific AP-1-DNA complex. The TpG step at the two terminals of the AP-1 site plays an important role in the structural adjustment of DNA by modifying the helical twist in the AP-1 bound state. We addressed the role of intrinsic motion of the bZIP domain in terms of opening and closing gripper motions of DNA binding helices, in target site recognition and binding of AP-1 factors. Our observations suggest that binding to the cognate motif in DNA is mainly accompanied with the precise adjustment of closing gripper motion of DNA binding helices of the bZIP domain.

Nonspecific vs. specific DNA binding free energetics of a transcription factor domain protein

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling.

Decoding the dual recognition mechanism of the glucocorticoid receptor for DNA and RNA: sequence versus shape

Transcription factors (TFs) regulate eukaryotic transcription through selective DNA-binding, can also specifically interact with RNA, which may present another layer of transcriptional control. The mechanisms of the TFs-DNA recognition are often well-characterised, while the details of TFs-RNA complexation are less understood. Here we investigate the dual recognition mechanism of the glucocorticoid receptor (GR), which interacts with similar affinities with consensus DNA and diverse RNA hairpin motifs but discriminates against uniform dsRNA. Using atomic molecular dynamics simulations, we demonstrate that the GR binding to nucleic acids requires a wide and shallow groove pocket. The protein effectively moulds its binding site within DNA major groove, which enables base-specific interactions. Contrary, the GR binding has little effect on the grooves geometry of RNA systems, most notably in uniform dsRNA. Instead, a hairpin motif in RNA yields a wide and shallow major groove pocket, allowing the protein to anchor itself through nonspecific electrostatic contacts with RNA backbone. Addition of a bulge increases RNA hairpin flexibility, which leads to a greater number of GR-RNA contacts and, thus, higher affinity. Thus, the combination of structural motifs defines the GR-RNA selective binding: a recognition mechanism, which may be shared by other zinc finger TFs.

How acidic amino acid residues facilitate DNA target site selection

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine-purine dinucleotides to adopt the BII conformation.

Mechanical determinants of chromatin topology and gene expression

The compaction of linear DNA into micrometer-sized nuclear boundaries involves the establishment of specific three-dimensional (3D) DNA structures complexed with histone proteins that form chromatin. The resulting structures modulate essential nuclear processes such as transcription, replication, and repair to facilitate or impede their multi-step progression and these contribute to dynamic modification of the 3D-genome organization. It is generally accepted that protein–protein and protein–DNA interactions form the basis of 3D-genome organization. However, the constant generation of mechanical forces, torques, and other stresses produced by various proteins translocating along DNA could be playing a larger role in genome organization than currently appreciated. Clearly, a thorough understanding of the mechanical determinants imposed by DNA transactions on the 3D organization of the genome is required. We provide here an overview of our current knowledge and highlight the importance of DNA and chromatin mechanics in gene expression.

Homologous basic helix–loop–helix transcription factors induce distinct deformations of torsionally-stressed DNA: a potential transcription regulation mechanism

Changing torsional restraints on DNA is essential for the regulation of transcription. Torsional stress, introduced by RNA polymerase, can propagate along chromatin facilitating topological transitions and modulating the specific binding of transcription factors (TFs) to DNA. Despite the importance, the mechanistic details on how torsional stress impacts the TFs-DNA complexation remain scarce. Herein, we address the impact of torsional stress on DNA complexation with homologous human basic helix–loop–helix (BHLH) hetero- and homodimers: MycMax, MadMax and MaxMax. The three TF dimers exhibit specificity towards the same DNA consensus sequence, the E-box response element, while regulating different transcriptional pathways. Using microseconds-long atomistic molecular dynamics simulations together with the torsional restraint that controls DNA total helical twist, we gradually over- and underwind naked and complexed DNA to a maximum of ± 5°/bp step. We observe that the binding of the BHLH dimers results in a similar increase in DNA torsional rigidity. However, under torsional stress the BHLH dimers induce distinct DNA deformations, characterised by changes in DNA grooves geometry and a significant asymmetric DNA bending. Supported by bioinformatics analyses, our data suggest that torsional stress may contribute to the execution of differential transcriptional programs of the homologous TFs by modulating their collaborative interactions.