DNA Shape Research Articles

Pou5f3, SoxB1, and Nanog remodel chromatin on high nucleosome affinity regions at zygotic genome activation.

The zebrafish embryo is transcriptionally mostly quiescent during the first 10 cell cycles, until the main wave of zygotic genome activation (ZGA) occurs, accompanied by fast chromatin remodeling. At ZGA, homologs of the mammalian stem cell transcription factors (TFs) Pou5f3, Nanog, and Sox19b bind to thousands of developmental enhancers to initiate transcription. So far, how these TFs influence chromatin dynamics at ZGA has remained unresolved. To address this question, we analyzed nucleosome positions in wild-type and maternal-zygotic (MZ) mutants for pou5f3 and nanog by MNase-seq. We show that Nanog, Sox19b, and Pou5f3 bind to the high nucleosome affinity regions (HNARs). HNARs are spanning over 600 bp, featuring high in vivo and predicted in vitro nucleosome occupancy and high predicted propeller twist DNA shape value. We suggest a two-step nucleosome destabilization-depletion model, in which the same intrinsic DNA properties of HNAR promote both high nucleosome occupancy and differential binding of TFs. In the first step, already before ZGA, Pou5f3 and Nanog destabilize nucleosomes at HNAR centers genome-wide. In the second step, post-ZGA, Nanog, Pou5f3, and SoxB1 maintain open chromatin state on the subset of HNARs, acting synergistically. Nanog binds to the HNAR center, whereas the Pou5f3 stabilizes the flanks. The HNAR model will provide a useful tool for genome regulatory studies in a variety of biological systems.

Clog and Release, and Reverse Motions of DNA in a Nanopore.

Motions of circular and linear DNA molecules of various lengths near a nanopore of 100 or 200 nm diameter were experimentally observed and investigated by fluorescence microscopy. The movement of DNA molecules through nanopores, known as translocation, is mainly driven by electric fields near and inside the pores. We found significant clogging of nanopores by DNA molecules, particularly by circular DNA and linear T4 DNA (165.65 kbp). Here, the probabilities of DNA clogging events, depending on the DNA length and shape—linear or circular—were determined. Furthermore, two distinct DNA motions were observed: clog and release by linear T4 DNA, and a reverse direction motion at the pore entrance by circular DNA, after which both molecules moved away from the pore. Finite element method-based numerical simulations were performed. The results indicated that DNA molecules with pores 100–200 nm in diameter were strongly influenced by opposing hydrodynamic streaming flow, which was further enhanced by bulky DNA configurations.

A De Novo Shape Motif Discovery Algorithm Reveals Preferences of Transcription Factors for DNA Shape Beyond Sequence Motifs.

DNA shape adds specificity to sequence motifs but has not been explored systematically outside this context. We hypothesized that DNA-binding proteins (DBPs) preferentially occupy DNA with specific structures ("shape motifs") regardless of whether or notthese correspond to high information content sequence motifs. We present ShapeMF, a Gibbs sampling algorithm that identifies de novo shape motifs. Using binding data from hundreds of invivo and invitro experiments, we show that most DBPs have shape motifs and can occupy these in the absence of sequence motifs. This "shape-only binding" is common for many DBPs and in regions co-bound by multiple DBPs. When shape and sequence motifs co-occur, they can be overlapping, flanking, or separated by consistent spacing. Finally, DBPs within the same protein family have different shape motifs, explaining their distinct genome-wide occupancy despite having similar sequence motifs. These results suggest that shape motifs not only complement sequence motifs but also facilitate recognition of DNA beyond conventionally defined sequence motifs.

ShapeGTB: the role of local DNA shape in prioritization of functional variants in human promoters with machine learning.

MotivationThe identification of functional sequence variations in regulatory DNA regions is one of the major challenges of modern genetics. Here, we report results of a combined multifactor analysis of properties characterizing functional sequence variants located in promoter regions of genes.ResultsWe demonstrate that GC-content of the local sequence fragments and local DNA shape features play significant role in prioritization of functional variants and outscore features related to histone modifications, transcription factors binding sites, or evolutionary conservation descriptors. Those observations allowed us to build specialized machine learning classifier identifying functional single nucleotide polymorphisms within promoter regions—ShapeGTB. We compared our method with more general tools predicting pathogenicity of all non-coding variants. ShapeGTB outperformed them by a wide margin (average precision 0.93 vs. 0.47–0.55). On the external validation set based on ClinVar database it displayed worse performance but was still competitive with other methods (average precision 0.47 vs. 0.23–0.42). Such results suggest unique characteristics of mutations located within promoter regions and are a promising signal for the development of more accurate variant prioritization tools in the future.

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR’s core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

Active site residue identity regulates cleavage preference of LAGLIDADG homing endonucleases.

LAGLIDADG homing endonucleases (meganucleases) are site-specific mobile endonucleases that can be adapted for genome-editing applications. However, one problem when reprogramming meganucleases on non-native substrates is indirect readout of DNA shape and flexibility at the central 4 bases where cleavage occurs. To understand how the meganuclease active site regulates DNA cleavage, we used functional selections and deep sequencing to profile the fitness landscape of 1600 I-LtrI and I-OnuI active site variants individually challenged with 67 substrates with central 4 base substitutions. The wild-type active site was not optimal for cleavage on many substrates, including the native I-LtrI and I-OnuI targets. Novel combinations of active site residues not observed in known meganucleases supported activity on substrates poorly cleaved by the wild-type enzymes. Strikingly, combinations of E or D substitutions in the two metal-binding residues greatly influenced cleavage activity, and E184D variants had a broadened cleavage profile. Analyses of I-LtrI E184D and the wild-type proteins co-crystallized with the non-cognate AACC central 4 sequence revealed structural differences that correlated with kinetic constants for cleavage of individual DNA strands. Optimizing meganuclease active sites to enhance cleavage of non-native central 4 target sites is a straightforward addition to engineering workflows that will expand genome-editing applications.

Extending the Code of Sequence Readout by Gene Regulatory Proteins: The Role of Hoogsteen Base Pairing in p53-DNA Recognition

The tumor suppressor p53 regulates the transcription of a myriad of target genes. In this issue of Structure, Golovenko etal. (2018) elucidate the role of non-canonical Hoogsteen A/T base pairs in p53-DNA recognition, highlighting DNA shape as an additional layer of sequence readout.

Bend-Induced Twist Waves and the Structure of Nucleosomal DNA.

Recent work indicates that twist-bend coupling plays an important role in DNA micromechanics. Here we investigate its effect on bent DNA. We provide an analytical solution of the minimum-energy shape of circular DNA, showing that twist-bend coupling induces sinusoidal twist waves. This solution is in excellent agreement with both coarse-grained simulations of minicircles and nucleosomal DNA data, which is bent and wrapped around histone proteins in a superhelical conformation. Our analysis shows that the observed twist oscillation in nucleosomal DNA, so far attributed to the interaction with the histone proteins, is an intrinsic feature of free bent DNA, and should be observable in other protein-DNA complexes.

Intrinsic DNA Shape Accounts for Affinity Differences between Hox-Cofactor Binding Sites

Transcription factors bind to their binding sites over a wide range of affinities, yet how differences in affinity are encoded in DNA sequences is not well understood. Here, we report X-ray crystal structures of four heterodimers of the Hox protein AbdominalB bound with its cofactor Extradenticle to four target DNA molecules that differ in affinity by up to ∼20-fold. Remarkably, despite large differences in affinity, the overall structures are very similar in all four complexes. In contrast, the predicted shapes of the DNA binding sites (i.e., the intrinsic DNA shape) in the absence of bound protein are strikingly different from each other and correlate with affinity: binding sites that must change conformations upon protein binding have lower affinities than binding sites that have more optimal conformations prior to binding. Together, these observations suggest that intrinsic differences in DNA shape provide a robust mechanism for modulating affinity without affecting other protein-DNA interactions.

Shared nucleotide flanks confer transcriptional competency to bZip core motifs.

Sequence-specific DNA binding recruits transcription factors (TFs) to the genome to regulate gene expression. Here, we perform high resolution mapping of CEBP proteins to determine how sequence dictates genomic occupancy. We demonstrate a fundamental difference between the sequence repertoire utilized by CEBPs in vivo versus the palindromic sequence preference reported by classical in vitro models, by identifying a palindromic motif at <1% of the genomic binding sites. On the native genome, CEBPs bind a diversity of related 10 bp sequences resulting from the fusion of degenerate and canonical half-sites. Altered DNA specificity of CEBPs in cells occurs through heterodimerization with other bZip TFs, and approximately 40% of CEBP-binding sites in primary human cells harbor motifs characteristic of CEBP heterodimers. In addition, we uncover an important role for sequence bias at core-motif-flanking bases for CEBPs and demonstrate that flanking bases regulate motif function across mammalian bZip TFs. Favorable flanking bases confer efficient TF occupancy and transcriptional activity, and DNA shape may explain how the flanks alter TF binding. Importantly, motif optimization within the 10-mer is strongly correlated with cell-type-independent recruitment of CEBPβ, providing key insight into how sequence sub-optimization affects genomic occupancy of widely expressed CEBPs across cell types.

New Insights into the Role of DNA Shape on Its Recognition by p53 Proteins

The tumor suppressor p53 acts as a transcription factor recognizing diverse DNA response elements (REs). Previous structural studies of p53-DNA complexes revealed non-canonical Hoogsteen geometry of A/T base pairs at conserved CATG motifs leading to changes in DNA shape and its interface with p53. To study the effects of DNA shape on binding characteristics, we designed REs with modified base pairs "locked" into either Hoogsteen or Watson-Crick form. Here we present crystal structures of these complexes and their thermodynamic and kinetic parameters, demonstrating that complexes with Hoogsteen base pairs are stabilized relative to those with all-Watson-Crick base pairs. CATG motifs are abundant in p53REs such as GADD45 and p53R2 related to cell-cycle arrest and DNA repair. The high-resolution structures of these complexes validate their propensity to adopt the unique Hoogsteen-induced structure, thus providing insights into the functional role of DNA shape and broadening the mechanisms that contribute to DNA recognition by proteins.

Crystal Structure of Apo MEF2B Reveals New Insights in DNA Binding and Cofactor Interaction.

The myocyte enhancer factor 2 (MEF2) family of transcription factors plays important roles in developmental processes and adaptive responses. Although MEF2 proteins are known to bind DNA in the nucleus to regulate specific gene expression, there are reports that show that MEF2 also functions in the cytoplasm. Previous structural studies of MEF2 focused exclusively on DNA-bound MEF2 with and without various corepressors or coactivators. While these studies have established a comprehensive structural model of DNA recognition and cofactor recruitment by MEF2, the structure of MEF2 not bound to DNA, which include cytoplasmic MEF2 and free MEF2 in the nucleus, is unknown. Here we determined the structure of the MADS-box/MEF2 domain of MEF2B without DNA nor cofactor. The Apo structure of MEF2B reveals a largely preformed DNA binding interface that may be important for recognizing the shape of DNA from the minor groove side. In addition, our structure also reveals that the C-terminal helix of the MEF2-specific domain could flip up to bind to the hydrophobic groove that serves as the binding sites of MEF2 transcription cofactors. These observations shed new insights into DNA binding and cofactor interaction by MEF2 proteins.

DNA sequence and shape are predictive for meiotic crossovers throughout the plant kingdom.

A better understanding of genomic features influencing the location of meiotic crossovers (COs) in plant species is both of fundamental importance and of practical relevance for plant breeding. Using CO positions with sufficiently high resolution from four plant species [Arabidopsis thaliana, Solanum lycopersicum (tomato), Zea mays (maize) and Oryza sativa (rice)] we have trained machine-learning models to predict the susceptibility to CO formation. Our results show that CO occurrence within various plant genomes can be predicted by DNA sequence and shape features. Several features related to genome content and to genomic accessibility were consistently either positively or negatively related to COs in all four species. Other features were found as predictive only in specific species. Gene annotation-related features were especially predictive for maize, whereas in tomato and Arabidopsis propeller twist and helical twist (DNA shape features) and AT/TA dinucleotides were found to be the most important. In rice, high roll (another DNA shape feature) and low CA dinucleotide frequency in particular were found to be associated with CO occurrence. The accuracy of our models was sufficient for Arabidopsis and rice (area under receiver operating characteristic curve, AUROC>0.5), and was high for tomato and maize (AUROC≫0.5), demonstrating that DNA sequence and shape are predictive for meiotic COs throughout the plant kingdom.

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.

Chemical and Biological Sensing Using Hybridization Chain Reaction.

Since the advent of its theoretical discovery more than 30 years ago, DNA nanotechnology has been used in a plethora of diverse applications in both the fundamental and applied sciences. The recent prominence of DNA-based technologies in the scientific community is largely due to the programmable features stored in its nucleobase composition and sequence, which allow it to assemble into highly advanced structures. DNA nanoassemblies are also highly controllable due to the precision of natural and artificial base-pairing, which can be manipulated by pH, temperature, metal ions, and solvent types. This programmability and molecular-level control have allowed scientists to create and utilize DNA nanostructures in one, two, and three dimensions (1D, 2D, and 3D). Initially, these 2D and 3D DNA lattices and shapes attracted a broad scientific audience because they are fundamentally captivating and structurally elegant; however, transforming these conceptual architectural blueprints into functional materials is essential for further advancements in the DNA nanotechnology field. Herein, the chemical and biological sensing applications of a 1D DNA self-assembly process known as hybridization chain reaction (HCR) are reviewed. HCR is a one-dimensional (1D) double stranded (ds) DNA assembly process initiated only in the presence of a specific short ssDNA (initiator) and two kinetically trapped DNA hairpin structures. HCR is considered an enzyme-free isothermal amplification process, which shows substantial promise and offers a wide range of applications for in situ chemical and biological sensing. Due to its modular nature, HCR can be programmed to activate only in the presence of highly specific biological and/or chemical stimuli. HCR can also be combined with different types of molecular reporters and detection approaches for various analytical readouts. While the long dsDNA HCR product may not be as structurally attractive as the 2D and 3D DNA networks, HCR is highly instrumental for applied biological, chemical, and environmental sciences, and has therefore been studied to foster a variety of objectives. In this review, we have focused on nucleic acid, protein, metabolite, and heavy metal ion detection using this 1D DNA nanotechnology via fluorescence, electrochemical, and nanoparticle-based methodologies.

Structural Insights into the CRTC2–CREB Complex Assembly on CRE

The cAMP response element (CRE) binding protein (CREB) is central in the transcription regulation by cAMP, and the CREB-regulated transcriptional coactivators (CRTCs) play critical roles in CREB-mediated transcription activation. Upon stimulation, CRTCs translocate into the nucleus and complex with CREB on CRE promoters to activate target gene transcription. Their physiological importance is underscored by their function in energy balance, long-term memory, longevity and other processes. The CREB binding domain on CRTCs has been mapped, which interacts with the CREB basic leucine zipper domain that also mediates interaction with CRE-containing DNA. We report here crystal structures of a complex containing the CRTC2 CREB binding domain, the CREB basic leucine zipper domain and a CRE-containing DNA. The structures revealed that CRTC and CREB form a 2:2 complex on CRE-containing DNA, and CRTC interacts with both CREB and DNA through highly conserved residues. Structure-guided functional studies revealed that both interactions are crucial for the complex assembly and CREB stabilization on DNA. Interestingly, we found that the CRTC–DNA interaction confers selectivity toward the intrinsic DNA shape, which may play a role in selective transcription activation of the CRE genes.

Direct Application of Vitis Vinifera (Grape) Leaves Extract on Cotton Fabric: A Potential to Prevent UV Induced Skin Problems

UV radiation causes various skin problems starting from reddening of skin to the change in the shape of DNA which further leads to the skin cancer. Considering harmful effect of UV radiation present study was designed to develop an environmental friendly UV protective finish for cotton fabric using Vitis vinifera leaves. Process parameters of the study were optimized using response surface methodology for getting maximum UV protection while maintaining other physical properties of the fabric. It was found that 9.038 % concentration of extract with 60 minutes of exhaustion time at 40.909oC exhaustion temperature were selected as optimum conditions for application of Vitis vinifera extract on cotton fabric. It can be concluded that application of Vitis vinifera extract on cotton fabric incorporates UV protection properties on cotton fabric which can be helpful to prevent various UV induced skin problems.

DNA‐binding properties of MafR, a global regulator of Enterococcus faecalis

Global transcriptional regulators play key roles during bacterial adaptation to environmental fluctuations. Protein MafR from Enterococcus faecalis was shown to activate the transcription of many genes on a genome‐wide scale. We proposed that MafR is a global regulator of the Mga/AtxA family. Here, we purified an untagged form of the MafR protein and found that it binds to linear double‐stranded DNAs in a nonsequence‐specific manner. Moreover, multiple MafR units (likely dimers) bind sequentially to the DNA molecule generating multimeric complexes. On DNAs that contain the promoter of the mafR gene, MafR recognizes a potentially curved DNA region. We discuss that a characteristic of the Mga/AtxA regulators might be their ability to recognize particular DNA shapes across the bacterial genomes.

Genome-wide determinants of sequence-specific DNA binding of general regulatory factors.

General regulatory factors (GRFs), such as Reb1, Abf1, Rap1, Mcm1, and Cbf1, positionally organize yeast chromatin through interactions with a core consensus DNA sequence. It is assumed that sequence recognition via direct base readout suffices for specificity and that spurious nonfunctional sites are rendered inaccessible by chromatin. We tested these assumptions through genome-wide mapping of GRFs in vivo and in purified biochemical systems at near–base pair (bp) resolution using several ChIP-exo–based assays. We find that computationally predicted DNA shape features (e.g., minor groove width, helix twist, base roll, and propeller twist) that are not defined by a unique consensus sequence are embedded in the nonunique portions of GRF motifs and contribute critically to sequence-specific binding. This dual source specificity occurs at GRF sites in promoter regions where chromatin organization starts. Outside of promoter regions, strong consensus sites lack the shape component and consequently lack an intrinsic ability to bind cognate GRFs, without regard to influences from chromatin. However, sites having a weak consensus and low intrinsic affinity do exist in these regions but are rendered inaccessible in a chromatin environment. Thus, GRF site-specificity is achieved through integration of favorable DNA sequence and shape readouts in promoter regions and by chromatin-based exclusion from fortuitous weak sites within gene bodies. This study further revealed a severe G/C nucleotide cross-linking selectivity inherent in all formaldehyde-based ChIP assays, which includes ChIP-seq. However, for most tested proteins, G/C selectivity did not appreciably affect binding site detection, although it does place limits on the quantitativeness of occupancy levels.

Comparison of discriminative motif optimization using matrix and DNA shape-based models

BackgroundTranscription factor (TF) binding site specificity is commonly represented by some form of matrix model in which the positions in the binding site are assumed to contribute independently to the site’s activity. The independence assumption is known to be an approximation, often a good one but sometimes poor. Alternative approaches have been developed that use k-mers (DNA “words” of length k) to account for the non-independence, and more recently DNA structural parameters have been incorporated into the models. ChIP-seq data are often used to assess the discriminatory power of motifs and to compare different models. However, to measure the improvement due to using more complex models, one must compare to optimized matrix models.ResultsWe describe a program “Discriminative Additive Model Optimization” (DAMO) that uses positive and negative examples, as in ChIP-seq data, and finds the additive position weight matrix (PWM) that maximizes the Area Under the Receiver Operating Characteristic Curve (AUROC). We compare to a recent study where structural parameters, serving as features in a gradient boosting classifier algorithm, are shown to improve the AUROC over JASPAR position frequency matrices (PFMs). In agreement with the previous results, we find that adding structural parameters gives the largest improvement, but most of the gain can be obtained by an optimized PWM and nearly all of the gain can be obtained with a di-nucleotide extension to the PWM.ConclusionTo appropriately compare different models for TF bind sites, optimized models must be used. PWMs and their extensions are good representations of binding specificity for most TFs, and more complex models, including the incorporation of DNA shape features and gradient boosting classifiers, provide only moderate improvements for a few TFs.