Propeller Twist Research Articles

Discerning the Role of DNA Sequence, Shape, and Flexibility in Recognition by Drosophila Transcription Factors.

The precise spatial and temporal orchestration of gene expression is crucial for the ontogeny of an organism and is mainly governed by transcription factors (TFs). The mechanism of recognition of cognate sites amid millions of base pairs in the genome by TFs is still incompletely understood. In this study, we focus on DNA sequence composition, shape, and flexibility preferences of 28 quintessential TFs from Drosophila melanogaster that are critical to development and body patterning mechanisms. Our study finds that TFs exhibit distinct predilections for DNA shape, flexibility, and sequence compositions in the proximity of transcription factor binding sites (TFBSs). Notably, certain zinc finger proteins prefer GC-rich areas with less negative propeller twist, while homeodomains mainly seek AT-rich regions with a more negative propeller twist at their sites. Intriguingly, while numerous cofactors share similar binding site preferences and bind closer to each other in the genome, some cofactors that have different preferences bind farther apart. These findings shed light on TF DNA recognition and provide novel insights into possible cofactor binding and transcriptional regulation mechanisms.

Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques

BackgroundNaturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers.MethodsWe conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation.ResultsSimilar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs.ConclusionsThese data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.Graphical abstract

Triplet Excimer Formation in a DNA Duplex with Silver Ion-Mediated Base Pairs.

The dynamics of excited electronic states in self-assembled structures formed between silver(I) ions and cytosine-containing DNA strands or monomeric cytosine derivatives were investigated by time-resolved infrared (TRIR) spectroscopy and quantum mechanical calculations. The steady-state and time-resolved spectra depend sensitively on the underlying structures, which change with pH and the nucleobase and silver ion concentrations. At pH ∼ 4 and low dC20 strand concentration, an intramolecularly folded i-motif is observed, in which protons, and not silver ions, mediate C-C base pairing. However, at the higher strand concentrations used in the TRIR measurements, dC20 strands associate pairwise to yield duplex structures containing C-Ag+-C base pairs with a high degree of propeller twisting. UV excitation of the silver ion-mediated duplex produces a long-lived excited state, which we assign to a triplet excimer state localized on a pair of stacked cytosines. The computational results indicate that the propeller-twisted motifs induced by metal-ion binding are responsible for the enhanced intersystem crossing that populates the triplet state and not a generic heavy atom effect. Although triplet excimer states have been discussed frequently as intermediates in the formation of cyclobutane pyrimidine dimers, we find neither computational nor experimental evidence for cytosine-cytosine photoproduct formation in the systems studied. These findings provide a rare demonstration of a long-lived triplet excited state that is formed in a significant yield in a DNA duplex, demonstrating that supramolecular structural changes induced by metal ion binding profoundly affect DNA photophysics.

DNA structural properties of DNA binding sites for 21 transcription factors in the mycobacterial genome.

Mycobacterium tuberculosis, the causative agent of tuberculosis, has evolved over time into a multidrug resistance strain that poses a serious global pandemic health threat. The ability to survive and remain dormant within the host macrophage relies on multiple transcription factors contributing to virulence. To date, very limited structural insights from crystallographic and NMR studies are available for TFs and TF-DNA binding events. Understanding the role of DNA structure in TF binding is critical to deciphering MTB pathogenicity and has yet to be resolved at the genome scale. In this work, we analyzed the compositional and conformational preference of 21 mycobacterial TFs, evident at their DNA binding sites, in local and global scales. Results suggest that most TFs prefer binding to genomic regions characterized by unique DNA structural signatures, namely, high electrostatic potential, narrow minor grooves, high propeller twist, helical twist, intrinsic curvature, and DNA rigidity compared to the flanking sequences. Additionally, preference for specific trinucleotide motifs, with clear periodic signals of tetranucleotide motifs, are observed in the vicinity of the TF-DNA interactions. Altogether, our study reports nuanced DNA shape and structural preferences of 21 TFs.

Cover Feature: A Homoleptic Diphosphatetrahedrane Nickel(0) Complex (Z. Anorg. Allg. Chem. 19/2022)

A three-bladed windmill illustrates the molecular structure of the title complex tris(diphosphatetrahedrane)nickel(0). The three sails of the windmill represent the three diphosphatetrahedrane molecules, which coordinate via their P–P bonds to a single nickel(0) atom forming the windshaft. The three blades illustrate the threefold rotation axis of the complex and the propeller twist of the ligands, which makes the complex chiral.(DOI: 10.1002/zaac.202200124).

Symphony of the DNA flexibility and sequence environment orchestrates p53 binding to its responsive elements.

The p53 tumor suppressor protein maintains the genome fidelity and integrity by modulating several cellular activities. It regulates these events by interacting with a heterogeneous set of response elements (REs) of regulatory genes in the background of chromatin configuration. At the p53-RE interface, both the base readout and torsional-flexibility of DNA account for high-affinity binding. However, DNA structure is an entanglement of a multitude of physicochemical features, both local and global structure should be considered for dealing with DNA-protein interactions. The goal of current research work is to conceptualize and abstract basic principles of p53-RE binding affinity as a function of structural alterations in DNA such as bending, twisting, and stretching flexibility and shape. For this purpose, we have exploited high throughput in-vitro relative affinity information of responsive elements and genome binding events of p53 from HT-Selex and ChIP-Seq experiments respectively. Our results confirm the role of torsional flexibility in p53 binding, and further, we reveal that DNA axial bending, stretching stiffness, propeller twist, and wedge angles are intimately linked to p53 binding affinity when compared to homeodomain, bZIP, and bHLH proteins. Besides, a similar DNA structural environment is observed in the distal sequences encompassing the actual binding sites of p53 cistrome genes. Additionally, we revealed that p53 cistrome target genes have unique promoter architecture, and the DNA flexibility of genomic sequences around REs in cancer and normal cell types display major differences. Altogether, our work provides a keynote on DNA structural features of REs that shape up the in-vitro and in-vivo high-affinity binding of the p53 transcription factor.

Phase Enabled Circular Dichroism Reversal in Twisted Bi‐Chiral Propeller Metamolecule Arrays

AbstractArtificial chiral structures play a crucial role in the realization of strong chiroptical responses and flexible light manipulation. Here, a unique bi‐chiral propeller metamolecule array is demonstrated in which the geometric chirality of the subunit with intracellular blades is opposite to that of the subunit with intercellular blades. Benefited from the special propeller twists, it is demonstrated that the propeller metamolecule arrays support the phase controlled chiral surface lattice resonances, whose amplitude can be tailored by tuning the coupling of chiral centers to change the interference phase. Specifically, the circular dichroism of the metamolecule array is readily reversed by rotating the propeller blades without changing the local geometry and the lattice parameters. The proposed bi‐chiral propeller metamolecule array and the exotic circular dichroism reversal are further demonstrated in experiments by adopting a nano‐kirigami fabrication method. This work provides a novel paradigm to investigate the profound chiroptical responses and manipulate the advanced spin states of light.

Translesion synthesis of the major nitrogen mustard-induced DNA lesion by human DNA polymerase η.

Nitrogen mustards are among the first modern anticancer chemotherapeutics that are still widely used as non-specific anticancer alkylating agents. While the mechanism of action of mustard drugs involves the generation of DNA interstrand cross-links, the predominant lesions produced by these drugs are nitrogen half-mustard-N7-dG (NHMG) adducts. The bulky major groove lesion NHMG, if left unrepaired, can be bypassed by translesion synthesis (TLS) DNA polymerases. However, studies of the TLS past NHMG have not been reported so far. Here, we present the first synthesis of an oligonucleotide containing a site-specific NHMG. We also report kinetic and structural characterization of human DNA polymerase η (polη) bypassing NHMG. The templating NHMG slows dCTP incorporation ∼130-fold, while it increases the misincorporation frequency ∼10-30-fold, highlighting the promutagenic nature of NHMG. A crystal structure of polη incorporating dCTP opposite NHMG shows a Watson-Crick NHMG:dCTP base pair with a large propeller twist angle. The nitrogen half-mustard moiety fits snugly into an open cleft created by the Arg61-Trp64 loop of polη, suggesting a role of the Arg61-Trp64 loop in accommodating bulky major groove adducts during lesion bypass. Overall, our results presented here to provide first insights into the TLS of the major DNA adduct formed by nitrogen mustard drugs.

Double-stranded RNA bending by AU-tract sequences.

Sequence-dependent structural deformations of the DNA double helix (dsDNA) have been extensively studied, where adenine tracts (A-tracts) provide a striking example for global bending in the molecule. However, in contrast to dsDNA, sequence-dependent structural features of dsRNA have received little attention. In this work, we demonstrate that the nucleotide sequence can induce a bend in a canonical Watson-Crick base-paired dsRNA helix. Using all-atom molecular dynamics simulations, we identified a sequence motif consisting of alternating adenines and uracils, or AU-tracts, that strongly bend the RNA double-helix. This finding was experimentally validated using atomic force microscopy imaging of dsRNA molecules designed to display macroscopic curvature via repetitions of phased AU-tract motifs. At the atomic level, this novel phenomenon originates from a localized compression of the dsRNA major groove and a large propeller twist at the position of the AU-tract. Moreover, the magnitude of the bending can be modulated by changing the length of the AU-tract. Altogether, our results demonstrate the possibility of modifying the dsRNA curvature by means of its nucleotide sequence, which may be exploited in the emerging field of RNA nanotechnology and might also constitute a natural mechanism for proteins to achieve recognition of specific dsRNA sequences.

Ultrafast excited state dynamics of silver ion-mediated cytosine-cytosine base pairs in metallo-DNA.

To better understand the nexus between structure and photophysics in metallo-DNA assemblies, the parallel-stranded duplex formed by the all-cytosine oligonucleotide, dC20, and silver nitrate was studied by circular dichroism (CD), femtosecond transient absorption spectroscopy, and time-dependent-density functional theory calculations. Silver(I) ions mediate Cytosine-Cytosine (CC) base pairs by coordinating to the N3 atoms of two cytosines. Although these silver(I) mediated CC base pairs resemble the proton-mediated CC base pairs found in i-motif DNA at first glance, a comparison of experimental and calculated CD spectra reveals that silver ion-mediated i-motif structures do not form. Instead, the parallel-stranded duplex formed between dC20 and silver ions is proposed to contain consecutive silver-mediated base pairs with high propeller twist-like ones seen in a recent crystal structure of an emissive, DNA-templated silver cluster. Femtosecond transient absorption measurements with broadband probing from the near UV to the near IR reveal an unusually long-lived (>10 ns) excited state in the dC20 silver ion complex that is not seen in dC20 in single-stranded or i-motif forms. This state is also absent in a concentrated solution of cytosine-silver ion complexes that are thought to assemble into planar ribbons or sheets that lack stacked silver(I) mediated CC base pairs. The large propeller twist angle present in metal-mediated base pairs may promote the formation of long-lived charged separated or triplet states in this metallo-DNA.

Delving into Eukaryotic Origins of Replication Using DNA Structural Features.

DNA replication in eukaryotes is an intricate process, which is precisely synchronized by a set of regulatory proteins, and the replication fork emanates from discrete sites on chromatin called origins of replication (Oris). These spots are considered as the gateway to chromosomal replication and are stereotyped by sequence motifs. The cognate sequences are noticeable in a small group of entire origin regions or totally absent across different metazoans. Alternatively, the use of DNA secondary structural features can provide additional information compared to the primary sequence. In this article, we report the trends in DNA sequence-based structural properties of origin sequences in nine eukaryotic systems representing different families of life. Biologically relevant DNA secondary structural properties, namely, stability, propeller twist, flexibility, and minor groove shape were studied in the sequences flanking replication start sites. Results indicate that Oris in yeasts show lower stability, more rigidity, and narrow minor groove preferences compared to genomic sequences surrounding them. Yeast Oris also show preference for A-tracts and the promoter element TATA box in the vicinity of replication start sites. On the contrary, Drosophila melanogaster, humans, and Arabidopsis thaliana do not have such features in their Oris, and instead, they show high preponderance of G-rich sequence motifs such as putative G-quadruplexes or i-motifs and CpG islands. Our extensive study applies the DNA structural feature computation to delve into origins of replication across organisms ranging from yeasts to mammals and including a plant. Insights from this study would be significant in understanding origin architecture and help in designing new algorithms for predicting DNA trans-acting factor recognition events.

Design of High Altitude Propeller Using Multilevel Optimization

An improved D’Angelo optimization framework based on the surrogate model and optimization algorithm is proposed to design and optimize a new high altitude propeller which is applied to the propulsion system of the stratospheric aircraft. The framework adopts a multilevel optimization strategy to determine a high efficiency and light weight propeller. The aerodynamic characteristics of S1223 airfoil are calculated by Computational Fluid Dynamics (CFD) method. The aerodynamic performance of the D’Angelo propeller which is obtained by the first-level optimization is predicted by Blade element momentum (BEM) theory. Then the D’Angelo propeller chord and twist angle, which are regarded as the initial conditions of the second-level optimization, are optimized to achieve the maximum efficiency by Multi-island genetic algorithm (MIGA). In addition, the non-dominated sorting genetic algorithm II (NSGA-II) is applied to maximize the propeller efficiency and minimize the propeller weight. And the Pareto frontier solutions about the efficiency and the blade area are obtained by NSGA-II. What is more, the aerodynamic characteristics of the D’Angelo propeller, the optimization propeller by MIGA, and the optimization propeller by NSGA-II are calculated by CFD simulation and compared with BEM results. It is shown that the CFD results are in fair agreement with BEM results and the aerodynamic performance of the NSGA-II propeller is prior to the MIGA propeller and is close to the D’Angelo propeller. Besides, the NSGA-II propeller is the lightest among them and can satisfy the cruise requirements of the high altitude propeller.

Deciphering the Gene Regulatory Landscape Encoded in DNA Biophysical Features.

SummaryGene regulation in higher organisms involves a sophisticated interplay between genetic and epigenetic mechanisms. Despite advances, the logic in selective usage of certain genomic regions as regulatory elements remains unclear. Here we show that the inherent biophysical properties of the DNA encode epigenetic state and the underlying regulatory potential. We find that the propeller twist (ProT) level is indicative of genomic location of the regulatory elements, their strength, the affinity landscape of transcription factors, and distribution in the nuclear 3D space. We experimentally show that ProT levels confer increased DNA flexibility and surface accessibility, and thus potentially primes usage of high ProT regions as regulatory elements. ProT levels also correlate with occurrence and phenotypic consequences of mutations. Interestingly, cell-fate switches involve a transient usage of low ProT regulatory elements. Altogether, our work provides unprecedented insights into the gene regulatory landscape encoded in the DNA biophysical features.

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes.

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2′-OH group. The weaker Watson–Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in 3 E-envelope C3' -endo, and B-type in 2 E-envelope C2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

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.

Ligand-induced DNA conformational changes in proflavine minor groove-bound complexes studied by molecular dynamics simulation

Background: Minor groove binding is a rate-limiting step in proflavine-DNA intercalation reaction. This step is believed also to be responsible for the sequence-dependent kinetics of proflavine binding to DNA. At the same time, most studies are focused on the final stage of the reaction – the intercalation complex, and there is a lack of data concerning the structure and stability of proflavine-DNA minor groove-bound complexes. Objectives: The objective of this study was to investigate the stability of proflavine minor groove-bound complexes with DNA oligonucleotides of different sequence by molecular dynamics simulation and to analyze the DNA conformational changes caused by the proflavine binding. Materials and methods: The molecular dynamics simulations of proflavine minor groove-bound complexes with poly(dA)·poly(dT) and poly(dCG)·poly(dCG) oligonucleotides of 30 bp length were done in program package AMBER12 with explicit water (SPC/E) and ions (NaCl 0.15 M) using force fields FF14SB for DNA and GAFF for ligand. The starting configurations of complexes were obtained by docking method in AutoDock 3.05. After multi-stage equilibration protocol, each system was simulated at T=300 K and p=1 bar for a 50 ns production phase. Then trajectories were post-processed in AMBERTools17 and VMD-1.9.3 packages. Results: Our simulations confirm that proflavine-DNA minor groove-bound complexes are stable in the 50 ns time range but there are some structural rearrangements in them with respect to the initial structures. The narrowing of the DNA minor groove is observed in the proflavine binding site. In proflavine-poly(dCG)·poly(dCG) complex it is more pronounced and is accompanied by the BI/BII transitions in DNA and the reorientation of ligand. In proflavine-poly(dA)·poly(dT) complex the specific intermolecular hydrogen bonds are formed, which are optimized by the changes in opening and propeller twisting of involved AT-base pairs. Complexes are stabilized by the van der Waals and hydrophobic interactions, which are more favorable in the proflavine-poly(dA)·poly(dT) complex. Conclusions: Our results show that the binding of proflavine to a minor groove of DNA induces the conformational changes in the DNA that are important for the resulting complex stability.

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.

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.

Systematic prediction of DNA shape changes due to CpG methylation explains epigenetic effects on protein\u2013DNA binding

BackgroundDNA shape analysis has demonstrated the potential to reveal structure-based mechanisms of protein–DNA binding. However, information about the influence of chemical modification of DNA is limited. Cytosine methylation, the most frequent modification, represents the addition of a methyl group at the major groove edge of the cytosine base. In mammalian genomes, cytosine methylation most frequently occurs at CpG dinucleotides. In addition to changing the chemical signature of C/G base pairs, cytosine methylation can affect DNA structure. Since the original discovery of DNA methylation, major efforts have been made to understand its effect from a sequence perspective. Compared to unmethylated DNA, however, little structural information is available for methylated DNA, due to the limited number of experimentally determined structures. To achieve a better mechanistic understanding of the effect of CpG methylation on local DNA structure, we developed a high-throughput method, methyl-DNAshape, for predicting the effect of cytosine methylation on DNA shape.ResultsUsing our new method, we found that CpG methylation significantly altered local DNA shape. Four DNA shape features—helix twist, minor groove width, propeller twist, and roll—were considered in this analysis. Distinct distributions of effect size were observed for different features. Roll and propeller twist were the DNA shape features most strongly affected by CpG methylation with an effect size depending on the local sequence context. Methylation-induced changes in DNA shape were predictive of the measured rate of cleavage by DNase I and suggest a possible mechanism for some of the methylation sensitivities that were recently observed for human Pbx-Hox complexes.ConclusionsCpG methylation is an important epigenetic mark in the mammalian genome. Understanding its role in protein–DNA recognition can further our knowledge of gene regulation. Our high-throughput methyl-DNAshape method can be used to predict the effect of cytosine methylation on DNA shape and its subsequent influence on protein–DNA interactions. This approach overcomes the limited availability of experimental DNA structures that contain 5-methylcytosine.