Shape Readout Research Articles

Role of Shape Deformation of DNA-Binding Sites in Regulating the Efficiency and Specificity in Their Recognition by DNA-Binding Proteins.

Accurate transcription of genetic information is crucial, involving precise recognition of the binding motifs by DNA-binding proteins. While some proteins rely on short-range hydrophobic and hydrogen bonding interactions at binding sites, others employ a DNA shape readout mechanism for specific recognition. In this mechanism, variations in DNA shape at the binding motif resulted from either inherent flexibility or binding of proteins at adjacent sites are sensed and capitalized by the searching proteins to locate them specifically. Through extensive computer simulations, we investigate both scenarios to uncover the underlying mechanism and origin of specificity in the DNA shape readout mechanism. Our findings reveal that deformation in shape at the binding motif creates an entropy funnel, allowing information about altered shapes to manifest as fluctuations in minor groove widths. This signal enhances the efficiency of nonspecific search of nearby proteins by directing their movement toward the binding site, primarily driven by a gain in entropy. We propose this as a generic mechanism for DNA shape readout, where specificity arises from the alignment between the molecular frustration of the searching protein and the ruggedness of the entropic funnel governed by molecular features of the protein and arrangement of the DNA bases at the binding site, respectively.

Predicting DNA structure using a deep learning method

Understanding the mechanisms of protein-DNA binding is critical in comprehending gene regulation. Three-dimensional DNA structure, also described as DNA shape, plays a key role in these mechanisms. In this study, we present a deep learning-based method, Deep DNAshape, that fundamentally changes the current k-mer based high-throughput prediction of DNA shape features by accurately accounting for the influence of extended flanking regions, without the need for extensive molecular simulations or structural biology experiments. By using the Deep DNAshape method, DNA structural features can be predicted for any length and number of DNA sequences in a high-throughput manner, providing an understanding of the effects of flanking regions on DNA structure in a target region of a sequence. The Deep DNAshape method provides access to the influence of distant flanking regions on a region of interest. Our findings reveal that DNA shape readout mechanisms of a core target are quantitatively affected by flanking regions, including extended flanking regions, providing valuable insights into the detailed structural readout mechanisms of protein-DNA binding. Furthermore, when incorporated in machine learning models, the features generated by Deep DNAshape improve the model prediction accuracy. Collectively, Deep DNAshape can serve as versatile and powerful tool for diverse DNA structure-related studies.

Noncanonical binding of transcription factors: time to revisit specificity?

Transcription factors (TFs) are one of the most studied classes of DNA-binding proteins that have a direct functional impact on gene transcription and thus, on human physiology and disease. The mechanisms that TFs use for recognizing target DNA binding sites have been studied for nearly five decades, yet they remain poorly understood. It is classically assumed that a TF recognizes a specific sequence pattern, or motif, as its binding sites. However, recent studies are consistently finding examples of noncanonical binding, that is, TFs binding at sites that do not resemble their sequence motifs. Here we review the current literature on four major types of noncanonical TF binding, namely binding based on DNA shape readout, at Guanine-quadruplex structures, at repeat sequences, and bispecific binding. These examples point to a critical need for studies to unify our current observations, many of which are at odds with the "one TF, one motif" view, into a more comprehensive definition of the DNA-binding specificity of TFs.

Structural insights into human TFIIIC promoter recognition.

Transcription factor (TF) IIIC recruits RNA polymerase (Pol) III to most of its target genes. Recognition of intragenic A- and B-box motifs in transfer RNA (tRNA) genes by TFIIIC modules τA and τB is the first critical step for tRNA synthesis but is mechanistically poorly understood. Here, we report cryo-electron microscopy structures of the six-subunit human TFIIIC complex unbound and bound to a tRNA gene. The τB module recognizes the B-box via DNA shape and sequence readout through the assembly of multiple winged-helix domains. TFIIIC220 forms an integral part of both τA and τB connecting the two subcomplexes via a ~550-amino acid residue flexible linker. Our data provide a structural mechanism by which high-affinity B-box recognition anchors TFIIIC to promoter DNA and permits scanning for low-affinity A-boxes and TFIIIB for Pol III activation.

Structural insights into CodY activation and DNA recognition.

Virulence factors enable pathogenic bacteria to infect host cells, establish infection, and contribute to disease progressions. In Gram-positive pathogens such as Staphylococcus aureus (Sa) and Enterococcus faecalis (Ef), the pleiotropic transcription factor CodY plays a key role in integrating metabolism and virulence factor expression. However, to date, the structural mechanisms of CodY activation and DNA recognition are not understood. Here, we report the crystal structures of CodY from Sa and Ef in their ligand-free form and their ligand-bound form complexed with DNA. Binding of the ligands-branched chain amino acids and GTP-induces conformational changes in the form of helical shifts that propagate to the homodimer interface and reorient the linker helices and DNA binding domains. DNA binding is mediated by a non-canonical recognition mechanism dictated by DNA shape readout. Furthermore, two CodY dimers bind to two overlapping binding sites in a highly cooperative manner facilitated by cross-dimer interactions and minor groove deformation. Our structural and biochemical data explain how CodY can bind a wide range of substrates, a hallmark of many pleiotropic transcription factors. These data contribute to a better understanding of the mechanisms underlying virulence activation in important human pathogens.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)6GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKCC-type MTFs (MC-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. MC-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which MC-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of MC-MTFs). A special focus of our review is on the structural prerequisites of MC-MTFs that enable the specific recognition of target genes.

MELD-DNA: Accurate structural prediction of protein-DNA interactions.

Protein-DNA interaction is one of the most important aspects of biology. It is a main step of gene regulation in all organisms and is a potential target for therapeutic and biotech applications. Protein-DNA interactions are distinct from protein-protein interactions due to several factors: 1) DNA sequence and specific interactions (base readout); 2) DNA shape and its deformability (shape readout); 3) Exploration and searching on the DNA for alternative binding states. Such factors dictate affinity and specificity of transcription factors. Traditionally, these interactions were only predicted at the sequence level (e.g. using position weight matrices) and more recently using docking and machine learning models. Here we present MELD-DNA, a too that combines cheap information from experiment (which can be noisy, ambiguous and sparse) with physics-based molecular mechanics force fields using bayesian inference. MELD will accelerate sampling of the binding process by temperature and hamiltonian replica exchanging that is guided by experimental data. We have tested the tool on a diverse set of 15 systems from various TF families and we observed native-like bound structures most of them with room for more improvements. We also compare and benchmark with the state-of-the-art in docking (HADDOCK) and ML (RosettaFold2NA). In this study we discuss several issues in the available methods and techniques that need to be addressed for more accurate predictions. We believe MELD-DNA can be an essential tool in the community from drug design, regulation and fundamental research.

Brassinosteroid-induced gene repression requires specific and tight promoter binding of BIL1/BZR1 via DNA shape readout.

BRZ-INSENSITIVE-LONG 1 (BIL1)/BRASSINAZOLE-RESISTANT 1 (BZR1) and its homologues are plant-specific transcription factors that convert the signalling of the phytohormones brassinosteroids (BRs) to transcriptional responses, thus controlling various physiological processes in plants. Although BIL1/BZR1 upregulates some BR-responsive genes and downregulates others, the molecular mechanism underlying the dual roles of BIL1/BZR1 is still poorly understood. Here we show that BR-responsive transcriptional repression by BIL1/BZR1 requires the tight binding of BIL1/BZR1 alone to the 10 bp elements of DNA fragments containing the known 6 bp core-binding motifs at the centre. Furthermore, biochemical and structural evidence demonstrates that the selectivity for two nucleobases flanking the core motifs is realized by the DNA shape readout of BIL1/BZR1 without direct recognition of the nucleobases. These results elucidate the molecular and structural basis of transcriptional repression by BIL1/BZR1 and contribute to further understanding of the dual roles of BIL1/BZR1 in BR-responsive gene regulation.

Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation.

Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.

Defining the Critical DNA Features Targeted by RNA Polymerase I Core Factor

RNA Polymerase I(Pol I) is one of three essential DNA dependent RNA polymerases in eukaryotes and is responsible for synthesizing ribosomal RNA. A critical and essential Pol I transcription factor in yeast is Core Factor (CF) which binds to a ~24 bp region in the rDNA promoter called the Core Element (CE). CF plays fundamental roles in the Pol I transcription process helping to recruit Pol I and open Pol I promoter DNA before initiation. Previously, it was not yet known how CF precisely recognized the CE. When interacting with DNA, proteins use two main mechanisms: i) base‐pair readout, and ii) shape/structural readout. Base‐readout is the most common and is dominated by hydrogen bonding between the amino acid residues and base‐pair hydrogen bond donors and acceptors. The second mechanism centers around DNA shape readout where proteins target specific DNA features such as curvature, bendability, and groove width. Previous studies from our lab have shown that CF and its human orthologue, Selectivity Factor 1 (SL1), use an evolutionarily conserved mechanism to target DNA which is governed by interactions with the GC minor groove, a unique surface rarely targeted by DNA‐binding proteins. To further understand the extent of structural recognition as well as specific DNA features CF may be using to interact with CE, we have employed a variety of selection‐based methods to resolve the structural rules governing CF’s interaction with DNA. Our findings are consistent with a model that CF‐CE interaction is governed primarily by DNA shape‐based structural features rather than sequence.

The Role of Charge Density Coupled DNA Bending in Transcription Factor Sequence Binding Specificity: A Generic Mechanism for Indirect Readout.

The accurate reading of genetic information during transcription is essential for the expression of genes. Sequence binding specificity very often is attributed to short-range, usually specific interactions between amino acid residues and individual nucleotide bases through hydrogen bonding or hydrophobic contacts: "base readout" (direct readout). In contrast, many proteins recognize DNA sequences in an alternative fashion via "shape readout" (indirect readout), where many elements of the DNA sequence cooperate to localize the transcription factor. In this study, we use a coarse-grained protein-DNA model to investigate the origin of the sequence specificity of the protein PU.1 binding to its binding sites for a series of DNA sequences. We find that the binding specificity of PU.1 is achieved primarily via a nonspecific electrostatically driven DNA mechanism involving the change in the elastic properties of the DNA. To understand the underlying mechanism, we analyze how the mechanical properties of DNA change in relation to the location of the PU.1 bound along DNA. The simulations first show that electrostatic interactions between PU.1 and DNA can cause complex DNA conformational/dynamics changes. Using a semiflexible polymer theory, we find that PU.1 influences the DNA dynamics through a second-order mechanical effect. When PU.1 binds nonspecifically to the DNA via electrostatics, the DNA becomes stiffer and the protein slides along DNA in a search mode. In contrast, once the protein finds its specific binding site, the DNA becomes softer there. PU.1 thus locks into place through configurational entropy effects, which we suggest is a generic mechanism for indirect readout.

Cell-free genomics reveal intrinsic, cooperative and competitive determinants of chromatin interactions.

Metazoan transcription factors distinguish their response elements from a large excess of similar sequences. We explored underlying principles of DNA shape read-out and factor cooperativity in chromatin using a unique experimental system. We reconstituted chromatin on Drosophila genomes in extracts of preblastoderm embryos, mimicking the naïve state of the zygotic genome prior to developmental transcription activation. We then compared the intrinsic binding specificities of three recombinant transcription factors, alone and in combination, with GA-rich recognition sequences genome-wide. For MSL2, all functional elements reside on the X chromosome, allowing to distinguish physiological elements from non-functional ‘decoy’ sites. The physiological binding profile of MSL2 is approximated through interaction with other factors: cooperativity with CLAMP and competition with GAF, which sculpts the profile by occluding non-functional sites. An extended DNA shape signature is differentially read out in chromatin. Our results reveal novel aspects of target selection in a complex chromatin environment.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq

Key messageWe studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode.MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

Transcription Factor Binding Affinities and DNA Shape Readout.

SummaryAn essential event in gene regulation is the binding of a transcription factor (TF) to its target DNA. Models considering the interactions between the TF and the DNA geometry proved to be successful approaches to describe this binding event, while conserving data interpretability. However, a direct characterization of the DNA shape contribution to binding is still missing due to the lack of accurate and large-scale binding affinity data. Here, we use a binding assay we recently established to measure with high sensitivity the binding specificities of 13 Drosophila TFs, including dinucleotide dependencies to capture non-independent amino acid-base interactions. Correlating the binding affinities with all DNA shape features, we find that shape readout is widely used by these factors. A shape readout/TF-DNA complex structure analysis validates our approach while providing biological insights such as positively charged or highly polar amino acids often contact nucleotides that exhibit strong shape readout.

Landscape of DNA binding signatures of myocyte enhancer factor-2B reveals a unique interplay of base and shape readout.

Myocyte enhancer factor-2B (MEF2B) has the unique capability of binding to its DNA target sites with a degenerate motif, while still functioning as a gene-specific transcriptional regulator. Identifying its DNA targets is crucial given regulatory roles exerted by members of the MEF2 family and MEF2B’s involvement in B-cell lymphoma. Analyzing structural data and SELEX-seq experimental results, we deduced the DNA sequence and shape determinants of MEF2B target sites on a high-throughput basis in vitro for wild-type and mutant proteins. Quantitative modeling of MEF2B binding affinities and computational simulations exposed the DNA readout mechanisms of MEF2B. The resulting binding signature of MEF2B revealed distinct intricacies of DNA recognition compared to other transcription factors. MEF2B uses base readout at its half-sites combined with shape readout at the center of its degenerate motif, where A-tract polarity dictates nuances of binding. The predominant role of shape readout at the center of the core motif, with most contacts formed in the minor groove, differs from previously observed protein–DNA readout modes. MEF2B, therefore, represents a unique protein for studies of the role of DNA shape in achieving binding specificity. MEF2B–DNA recognition mechanisms are likely representative for other members of the MEF2 family.

Atomistic insight into sequence-directed DNA bending and minicircle formation propensity in the absence and presence of phased A-tracts.

Bending of double-stranded (ds) DNA plays a crucial role in many important biological processes and is relevant for nanotechnological applications. Among all the elements that have been studied in relation to dsDNA bending, A-tracts stand out as one of the most controversial. The "ApA wedge" theory was disproved when a series of linear polynucleotides containing phased 5'-A4T4-3' or 5'-T4A4-3' runs were shown to be bent or straight, respectively, and crystallographic evidence revealed that A-tracts are unbent. Furthermore, some of the smallest dsDNA minicircles described to date (~ 100bp in size) lack A-tracts and are subjected to varying levels of torsional stress. Representative DNA sequences from this experimental background were modeled in atomic detail and their dynamic behavior was simulated over hundreds of nanoseconds using the AMBER force field ParmBSC1. Subsequent analysis of the resulting trajectories allowed us to (i) unambiguously establish the location of the bends in all cases; (ii) identify the structural elements that are directly responsible for the macroscopically detected curvature; and (iii) reveal the importance not only of coherently summing the effects of the bending elements when they are in synchrony with the natural repeat of the helix (i.e. separated by an integral number of helical turns) but also when alternated with a half-integral separation of opposite effects. We conclude that the major determinant of the macroscopically observed bending is the proper grouping and phasing of the positive roll imposed by pyrimidine-purine (YR) steps and the negative or null roll characteristic of RY steps and A-tracts, respectively. This conclusion is in very good agreement with the structural parameters experimentally derived for much smaller DNA molecules either alone or as found in DNA-protein complexes. We expect that this work will pave the way for future studies on drug-induced DNA bending, DNA shape readout by transcription factors, structure of circular extrachromosomal DNA, and custom design of curved DNA origami scaffolds.

How B-DNA Dynamics Decipher Sequence-Selective Protein Recognition

The rules governing sequence-specific DNA–protein recognition are under a long-standing debate regarding the prevalence of base versus shape readout mechanisms to explain sequence specificity and of the conformational selection versus induced fit binding paradigms to explain binding-related conformational changes in DNA. Using a combination of atomistic simulations on a subset of representative sequences and mesoscopic simulations at the protein–DNA interactome level, we demonstrate the prevalence of the shape readout model in determining sequence-specificity and of the conformational selection paradigm in defining the general mechanism for binding-related conformational changes in DNA. Our results suggest that the DNA uses a double mechanism to adapt its structure to the protein: it moves along the easiest deformation modes to approach the bioactive conformation, while final adjustments require localized rearrangements at the base-pair step and backbone level. Our study highlights the large impact of B-DNA dynamics in modulating DNA–protein binding.

Crystal Structures of REF6 and Its Complex with DNA Reveal Diverse Recognition Mechanisms

RELATIVE OF EARLY FLOWERING 6 (REF6) is a DNA-sequence-specific H3K27me3/2 demethylase containing four zinc finger (ZnF) domains, targeting several thousand genes in Arabidopsis thaliana. The ZnF domains are essential for binding target genes, but the structural basis and mechanisms remain unclear. Here, we determined crystal structures of the ZnF domains and its dsDNA complex. It is composed of a unique REF6-family-specific half-cross-braced ZnF (RCZ) domain and two C2H2-type ZnFs. DNA binding induces a profound conformational change of REF6 at the hinge region. Each REF6 recognizes six bases and DNA methylation reduces the binding affinity. Both the acidic region and basic region are crucial for the self-association of REF6. The binding ability and specificity are enhanced by combining base and shape readout, and cooperativity. The resulting conformational plasticity enables REF6 to function as a global transcriptional regulator that directly binds many diverse genes, revealing the structural basis of recognizing epigenetic modification.

Transcriptional activation by MafR, a global regulator of Enterococcus faecalis

Proteins that act as global transcriptional regulators play key roles in bacterial adaptation to new niches. These proteins recognize multiple DNA sites across the bacterial genome by different mechanisms. Enterococcus faecalis is able to survive in various niches of the human host, either as a commensal or as a leading cause of serious infections. Nonetheless, the regulatory pathways involved in its adaptive responses remain poorly understood. We reported previously that the MafR protein of E. faecalis causes genome-wide changes in the transcriptome. Here we demonstrate that MafR functions as a transcription activator. In vivo, MafR increased the activity of the P12294 and P11486 promoters and also the transcription levels of the two genes controlled by those promoters. These genes are predicted to encode a calcium-transporting P-type ATPase and a QueT transporter family protein, respectively. Thus, MafR could have a regulatory role in calcium homeostasis and queuosine synthesis. Furthermore, MafR recognized in vitro specific DNA sites that overlap the −35 element of each target promoter. The MafR binding sites exhibit a low sequence identity, suggesting that MafR uses a shape readout mechanism to achieve DNA-binding specificity.

The floral homeotic protein SEPALLATA3 recognizes target DNA sequences by shape readout involving a conserved arginine residue in the MADS-domain.

SEPALLATA3 of Arabidopsis thaliana is a MADS-domain transcription factor (TF) and a key regulator of flower development. MADS-domain proteins bind to sequences termed 'CArG-boxes' [consensus 5'-CC(A/T)6 GG-3']. Because only a fraction of the CArG-boxes in the Arabidopsis genome are bound by SEPALLATA3, more elaborate principles have to be discovered to better understand which features turn CArG-boxes into genuine recognition sites. Here, we investigate to what extent the shape of the DNA is involved in a 'shape readout' that contributes to the binding of SEPALLATA3. We determined invitro binding affinities of SEPALLATA3 to DNA probes that all contain the CArG-box motif, but differ in their predicted DNA shape. We found that binding affinity correlates well with a narrow minor groove of the DNA. Substitution of canonical bases with non-standard bases supports the hypothesis of minor groove shape readout by SEPALLATA3. Analysis of mutant SEPALLATA3 proteins further revealed that a highly conserved arginine residue, which is expected to contact the DNA minor groove, contributes significantly to the shape readout. Our studies show that the specific recognition of cis-regulatory elements by a plant MADS-domain TF, and by inference probably also of other TFs of this type, heavily depends on shape readout mechanisms.