Structure of an alternating-B DNA helix and its relationship to A-tract DNA

X-ray Structure Characterization of the Selective Recognition of AT Base Pair Sequences

The rational design of small molecules that target specific DNA sequences is a promising strategy to modulate gene expression. This report focuses on a diamidinobenzimidazole compound, whose selective binding to the minor groove of AT DNA sequences holds broad significance in the molecular recognition of AT-rich human promoter sequences. The objective of this study is to provide a more detailed and systematized understanding, at an atomic level, of the molecular recognition mechanism of different AT-specific sequences by a rationally designed minor groove binder. The specialized method of X-ray crystallography was utilized to investigate how the sequence-dependent recognition properties in general, A-tract, and alternating AT sequences affect the binding of diamidinobenzimidazole in the DNA minor groove. While general and A-tract AT sequences give a narrower minor groove, the alternating AT sequences intrinsically have a wider minor groove which typically constricts upon binding. A strong and direct hydrogen bond between the N-H of the benzimidazole and an H-bond acceptor atom in the minor groove is essential for DNA recognition in all sequences described. In addition, the diamidine compound specifically utilizes an interfacial water molecule for its DNA binding. DNA complexes of AATT and AAAAAA recognition sites show that the diamidine compound can bind in two possible orientations with a preference for water-assisted hydrogen bonding at either cationic end. The complex structures of AAATTT, ATAT, ATATAT, and AAAA are bound in a singular orientation. Analysis of the helical parameters shows a minor groove expansion of about 1 Å across all the nonalternating DNA complexes. The results from this systematic approach will convey a greater understanding of the specific recognition of a diverse array of AT-rich sequences by small molecules and more insight into the design of small molecules with enhanced specificity to AT and mixed DNA sequences.

Structural insights into the functional divergence of WhiB-like proteins in Mycobacterium tuberculosis

WhiB7 represents a distinct subclass of transcription factors in the WhiB-Like (Wbl) family, a unique group of iron-sulfur (4Fe-4S] cluster-containing proteins exclusive to the phylum of Actinobacteria. In Mycobacterium tuberculosis (Mtb), WhiB7 interacts with domain 4 of the primary sigma factor (σ^A₄) in the RNA polymerase holoenzyme and activates genes involved in multiple drug resistance and redox homeostasis. Here, we report crystal structures of the WhiB7:σ^A₄ complex alone and bound to its target promoter DNA at 1.55-Å and 2.6-Å resolution, respectively. These structures show how WhiB7 regulates gene expression by interacting with both σ^A₄ and the AT-rich sequence upstream of the -35 promoter DNA via its C-terminal DNA-binding motif, the AT-hook. By combining comparative structural analysis of the two high-resolution σ^A₄-bound Wbl structures with molecular and biochemical approaches, we identify the structural basis of the functional divergence between the two distinct subclasses of Wbl proteins in Mtb.

Natural and Synthetic Oligoarylamides: Privileged Structures for Medical Applications

The term "privileged structure" refers to a single molecular substructure or scaffold that can serve as a starting point for high-affinity ligands for more than one receptor type. In this report, a hitherto overlooked group of privileged substructures is addressed, namely aromatic oligoamides, for which there are natural models in the form of cystobactamids, albicidin, distamycin A, netropsin, and others. The aromatic and heteroaromatic core, together with a flexible selection of substituents, form conformationally well-defined scaffolds capable of specifically binding to conformationally well-defined regions of biomacromolecules such as helices in proteins or DNA often by acting as helices mimics themselves. As such, these aromatic oligoamides have already been employed to inhibit protein-protein and nucleic acid-protein interactions. This article is the first to bring together the scattered knowledge about aromatic oligoamides in connection with biomedical applications.

Overlapping but Distinct Sequences Play Roles in the Insulator and Promoter Activities of the Drosophila BEAF-Dependent scs' Insulator

Chromatin domain insulators are thought to help partition the genome into genetic units called topologically associating domains (TADs). In Drosophila, TADs are often separated by inter-TAD regions containing active housekeeping genes and associated insulator binding proteins. This raises the question of whether insulator binding proteins are involved primarily in chromosomal TAD architecture or gene activation, or if these two activities are linked. The Boundary Element-Associated Factor of 32 kDa (BEAF-32, or BEAF for short) is usually found in inter-TADs. BEAF was discovered based on binding to the scs' insulator, and is important for the insulator activity of scs' and other BEAF binding sites. There are divergent promoters in scs' with a BEAF binding site by each. Here, we dissect the scs' insulator to identify DNA sequences important for insulator and promoter activity, focusing on the half of scs' with a high affinity BEAF binding site. We find that the BEAF binding site is important for both insulator and promoter activity, as is another sequence we refer to as LS4. Aside from that, different sequences play roles in insulator and promoter activity. So while there is overlap and BEAF is important for both, insulator and promoter activity can be separated.

Sequence-Dependent Deviations of Constrained DNA from Canonical B-Form

Decades of crystallographic and NMR studies have produced canonical structural models of short DNA. However, no experimental method so far has been able to test these models in vivo, where DNA is long and constrained by interactions with membranes, proteins, and other molecules. Here, we employ high-resolution frequency-modulation AFM to image single long poly(dA)-poly(dT), poly(dG)-poly(dC), and lambda DNA molecules interacting with an underlying substrate that emulates the effect of biological constraints on molecular structure. We find systematic sequence-dependent variations in groove dimensions, indicating that the structure of DNA subject to realistic interactions may differ profoundly from canonical models. These findings highlight the value of AFM as a unique, single molecule characterization tool.

Comparative in silico study of congocidine congeners as potential inhibitors of African swine fever virus

African swine fever virus (ASFV) infection is fatal in domesticated pigs, with a mortality rate approaching 100%. This may result in economic losses and threats to food security. Currently, there are no approved vaccines or antiviral therapies for ASFV. Therefore, in this study, we evaluated congocidine congeners and a tris-benzimidazole as potential inhibitors of ASFV transcription using an in silico approach. We applied redocking of congocidine and docking of its congeners and a tris-benzimidazole to a receptor containing B-DNA with AT-motifs as a target to mimic conserved ASFV late gene promoters. Subsequently, the binding scores of DNA-ligand docked complexes were evaluated and their binding affinity was estimated. Molecular dynamics (MD) simulation was then used to assess ligand behavior within the minor groove. From our results, it is evident the less toxic congocidine congeners and tris-benzimidazole could dock to AT-rich regions significantly. Additionally, the predicted binding affinities had suitable values comparable to other experimentally determined minor groove binders, MD simulation of the docked DNA-ligand complexes and subsequent molecular trajectory visualization further showed that the ligands remained embedded in the minor groove during the time course of simulation, indicating that these ligands may have potential applications in abrogating ASFV transcription.

Adenine radicals generated in alternating AT duplexes by direct absorption of low-energy UV radiation

There is increasing evidence that the direct absorption of photons with energies that are lower than the ionization potential of nucleobases may result in oxidative damage to DNA. The present work, which combines nanosecond transient absorption spectroscopy and quantum mechanical calculations, studies this process in alternating adenine-thymine duplexes (AT)n. We show that the one-photon ionization quantum yield of (AT)10 at 266 nm (4.66 eV) is (1.5 ± 0.3) × 10-3. According to our PCM/TD-DFT calculations carried out on model duplexes composed of two base pairs, (AT)1 and (TA)1, simultaneous base pairing and stacking does not induce important changes in the absorption spectra of the adenine radical cation and deprotonated radical. The adenine radicals, thus identified in the time-resolved spectra, disappear with a lifetime of 2.5 ms, giving rise to a reaction product that absorbs at 350 nm. In parallel, the fingerprint of reaction intermediates other than radicals, formed directly from singlet excited states and assigned to AT/TA dimers, is detected at shorter wavelengths. PCM/TD-DFT calculations are carried out to map the pathways leading to such species and to characterize their absorption spectra; we find that, in addition to the path leading to the well-known TA* photoproduct, an AT photo-dimerization path may be operative in duplexes.

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

Motivation

The 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.

Results

We 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.

Non-programmed transcriptional frameshifting is common and highly RNA polymerase type-dependent

Background

The viral or host systems for a gene expression assume repeatability of the process and high quality of the protein product. Since level and fidelity of transcription primarily determines the overall efficiency, all factors contributing to their decrease should be identified and optimized. Among many observed processes, non-programmed insertion/deletion (indel) of nucleotide during transcription (slippage) occurring at homopolymeric A/T sequences within a gene can considerably impact its expression. To date, no comparative study of the most utilized Escherichia coli and T7 bacteriophage RNA polymerases (RNAP) propensity for this type of erroneous mRNA synthesis has been reported. To address this issue we evaluated the influence of shift-prone A/T sequences by assessing indel-dependent phenotypic changes. RNAP-specific expression profile was examined using two of the most potent promoters, P_araBAD of E. coli and φ10 of phage T7.

Results

Here we report on the first systematic study on requirements for efficient transcriptional slippage by T7 phage and cellular RNAPs considering three parameters: homopolymer length, template type, and frameshift directionality preferences. Using a series of out-of-frame gfp reporter genes fused to a variety of A/T homopolymeric sequences we show that T7 RNAP has an exceptional potential for generating frameshifts and is capable of slipping on as few as three adenine or four thymidine residues in a row, in a flanking sequence-dependent manner. In contrast, bacterial RNAP exhibits a relatively low ability to baypass indel mutations and requires a run of at least 7 tymidine and even more adenine residues. This difference comes from involvement of various intrinsic proofreading properties. Our studies demonstrate distinct preference towards a specific homopolymer in slippage induction. Whereas insertion slippage performed by T7 RNAP (but not deletion) occurs tendentiously on poly(A) rather than on poly(T) runs, strong bias towards poly(T) for the host RNAP is observed.

Conclusions

Intrinsic RNAP slippage properties involve trade-offs between accuracy, speed and processivity of transcription. Viral T7 RNAP manifests far greater inclinations to the transcriptional slippage than E. coli RNAP. This possibly plays an important role in driving bacteriophage adaptation and therefore could be considered as beneficial. However, from biotechnological and experimental viewpoint, this might create some problems, and strongly argues for employing bacterial expression systems, stocked with proofreading mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1034-4) contains supplementary material, which is available to authorized users.

Sequence-Dependent Persistence Length of Long DNA

Using a high-throughput genome-mapping approach, we obtained circa 50 million measurements of the extension of internal human DNA segments in a 41  nm×41  nm nanochannel. The underlying DNA sequences, obtained by mapping to the reference human genome, are 2.5-393 kilobase pairs long and contain percent GC contents between 32.5% and 60%. Using Odijk's theory for a channel-confined wormlike chain, these data reveal that the DNA persistence length increases by almost 20% as the percent GC content increases. The increased persistence length is rationalized by a model, containing no adjustable parameters, that treats the DNA as a statistical terpolymer with a sequence-dependent intrinsic persistence length and a sequence-independent electrostatic persistence length.

SelexGLM differentiates androgen and glucocorticoid receptor DNA-binding preference over an extended binding site

The DNA-binding interfaces of the androgen (AR) and glucocorticoid (GR) receptors are virtually identical, yet these transcription factors share only about a third of their genomic binding sites and regulate similarly distinct sets of target genes. To address this paradox, we determined the intrinsic specificities of the AR and GR DNA-binding domains using a refined version of SELEX-seq. We developed an algorithm, SelexGLM, that quantifies binding specificity over a large (31-bp) binding site by iteratively fitting a feature-based generalized linear model to SELEX probe counts. This analysis revealed that the DNA-binding preferences of AR and GR homodimers differ significantly, both within and outside the 15-bp core binding site. The relative preference between the two factors can be tuned over a wide range by changing the DNA sequence, with AR more sensitive to sequence changes than GR. The specificity of AR extends to the regions flanking the core 15-bp site, where isothermal calorimetry measurements reveal that affinity is augmented by enthalpy-driven readout of poly(A) sequences associated with narrowed minor groove width. We conclude that the increased specificity of AR is correlated with more enthalpy-driven binding than GR. The binding models help explain differences in AR and GR genomic binding and provide a biophysical rationale for how promiscuous binding by GR allows functional substitution for AR in some castration-resistant prostate cancers.

DNA's Chiral Spine of Hydration

The iconic helical structure of DNA is stabilized by the solvation environment, where a change in the hydration state can lead to dramatic changes to the DNA structure. X-ray diffraction experiments at cryogenic temperatures have shown crystallographic water molecules in the minor groove of DNA, which has led to the notion of a spine of hydration of DNA. Here, chiral nonlinear vibrational spectroscopy of two DNA sequences shows that not only do such structural water molecules exist in solution at ambient conditions but that they form a chiral superstructure: a chiral spine of hydration. This is the first observation of a chiral water superstructure templated by a biomolecule. While the biological relevance of a chiral spine of hydration is unknown, the method provides a direct way to interrogate the properties of the hydration environment of DNA and water in biological settings without the use of labels.

5-Formylcytosine does not change the global structure of DNA

The mechanism by which the recently identified DNA modification 5-formylcytosine (^fC) is recognized by epigenetic writer and reader proteins is not known. Recently, an unusual DNA structure, F-DNA, has been proposed as the basis for enzyme recognition of clusters of ^fC. We used NMR and X-ray crystallography to compare several modified DNA duplexes with unmodified analogs and found that in the crystal state the duplexes all belong to the A family, whereas in solution they are all members of the B family. We found that, contrary to previous findings, ^fC does not significantly affect the structure of DNA, although there are modest local differences at the modification sites. Hence, global conformation changes are unlikely to account for the recognition of this modified base, and our structural data favor a mechanism that operates at base-pair resolution for the recognition of ^fC by epigenome-modifying enzymes.

2D-IR Spectroscopy Shows that Optimized DNA Minor Groove Binding of Hoechst33258 Follows an Induced Fit Model

The induced fit binding model describes a conformational change occurring when a small molecule binds to its biomacromolecular target. The result is enhanced noncovalent interactions between the ligand and biomolecule. Induced fit is well-established for small molecule-protein interactions, but its relevance to small molecule-DNA binding is less clear. We investigate the molecular determinants of Hoechst33258 binding to its preferred A-tract sequence relative to a suboptimal alternating A-T sequence. Results from two-dimensional infrared spectroscopy, which is sensitive to H-bonding and molecular structure changes, show that Hoechst33258 binding results in loss of the minor groove spine of hydration in both sequences, but an additional perturbation of the base propeller twists occurs in the A-tract binding region. This induced fit maximizes favorable ligand-DNA enthalpic contributions in the optimal binding case and demonstrates that controlling the molecular details that induce subtle changes in DNA structure may hold the key to designing next-generation DNA-binding molecules.

H-NS, Its Family Members and Their Regulation of Virulence Genes in Shigella Species

The histone-like nucleoid structuring protein (H-NS) has played a key role in shaping the evolution of Shigella spp., and provides the backdrop to the regulatory cascade that controls virulence by silencing many genes found on the large virulence plasmid. H-NS and its paralogue StpA are present in all four Shigella spp., but a second H-NS paralogue, Sfh, is found in the Shigella flexneri type strain 2457T, which is routinely used in studies of Shigella pathogenesis. While StpA and Sfh have been proposed to serve as “molecular backups” for H-NS, the apparent redundancy of these proteins is questioned by in vitro studies and work done in Escherichia coli. In this review, we describe the current understanding of the regulatory activities of the H-NS family members, the challenges associated with studying these proteins and their role in the regulation of virulence genes in Shigella.

Context based computational analysis and characterization of ARS consensus sequences (ACS) of Saccharomyces cerevisiae genome

Genome-wide experimental studies in Saccharomyces cerevisiae reveal that autonomous replicating sequence (ARS) requires an essential consensus sequence (ACS) for replication activity. Computational studies identified thousands of ACS like patterns in the genome. However, only a few hundreds of these sites act as replicating sites and the rest are considered as dormant or evolving sites. In a bid to understand the sequence makeup of replication sites, a content and context-based analysis was performed on a set of replicating ACS sequences that binds to origin-recognition complex (ORC) denoted as ORC-ACS and non-replicating ACS sequences (nrACS), that are not bound by ORC. In this study, DNA properties such as base composition, correlation, sequence dependent thermodynamic and DNA structural profiles, and their positions have been considered for characterizing ORC-ACS and nrACS. Analysis reveals that ORC-ACS depict marked differences in nucleotide composition and context features in its vicinity compared to nrACS. Interestingly, an A-rich motif was also discovered in ORC-ACS sequences within its nucleosome-free region. Profound changes in the conformational features, such as DNA helical twist, inclination angle and stacking energy between ORC-ACS and nrACS were observed. Distribution of ACS motifs in the non-coding segments points to the locations of ORC-ACS which are found far away from the adjacent gene start position compared to nrACS thereby enabling an accessible environment for ORC-proteins. Our attempt is novel in considering the contextual view of ACS and its flanking region along with nucleosome positioning in the S. cerevisiae genome and may be useful for any computational prediction scheme.

Structure elucidation of the Pribnow box consensus promoter sequence by racemic DNA crystallography

It has previously been shown that the use of racemic mixtures of naturally chiral macromolecules such as protein and DNA can significantly aid the crystallogenesis process, thereby addressing one of the major bottlenecks to structure determination by X-ray crystallographic methods-that of crystal growth. Although previous studies have provided convincing evidence of the applicability of the racemic crystallization technique to DNA through the study of well-characterized DNA structures, we sought to apply this method to a historically challenging DNA sequence. For this purpose we chose a non-self-complementary DNA duplex containing the biologically-relevant Pribnow box consensus sequence 'TATAAT'. Four racemic crystal structures of this previously un-crystallizable DNA target are reported (with resolutions in the range of 1.65-2.3 Å), with further crystallographic studies and structural analysis providing insight into the racemic crystallization process as well as structural details of this highly pertinent DNA sequence.

High-Energy Long-Lived Mixed Frenkel-Charge-Transfer Excitons: From Double Stranded (AT)n to Natural DNA

The electronic excited states populated upon absorption of UV photons by DNA are extensively studied in relation to the UV-induced damage to the genetic code. Here, we report a new unexpected relaxation pathway in adenine-thymine double-stranded structures (AT)n . Fluorescence measurements on (AT)n hairpins (six and ten base pairs) and duplexes (20 and 2000 base pairs) reveal the existence of an emission band peaking at approximately 320 nm and decaying on the nanosecond time scale. Time-dependent (TD)-DFT calculations, performed for two base pairs and exploring various relaxation pathways, allow the assignment of this emission band to excited states resulting from mixing between Frenkel excitons and adenine-to-thymine charge-transfer states. Emission from such high-energy long-lived mixed (HELM) states is in agreement with their fluorescence anisotropy (0.03), which is lower than that expected for π-π* states (≥0.1). An increase in the size of the system quenches π-π* fluorescence while enhancing HELM fluorescence. The latter process varies linearly with the hypochromism of the absorption spectra, both depending on the coupling between π-π* and charge-transfer states. Subsequently, we identify the common features between the HELM states of (AT)n structures with those reported previously for alternating (GC)n : high emission energy, low fluorescence anisotropy, nanosecond lifetimes, and sensitivity to conformational disorder. These features are also detected for calf thymus DNA in which HELM states could evolve toward reactive π-π* states, giving rise to delayed fluorescence.

Drug-DNA Interaction Studies of Acridone-Based Derivatives

N10-alkylated 2-bromoacridones are a novel series of potent antitumor compounds. DNA binding studies of these compounds were carried out using spectrophotometric titrations, Circular dichroism (CD) measurements using Calf Thymus DNA (CT DNA). The binding constants were identified at a range of K=0.3 to 3.9×10(5) M(-1) and the percentage of hypochromism from the spectral titrations at 28-54%. This study has identified a compound 9 with the good binding affinity of K=0.39768×10(5) M(-1) with CT DNA. Molecular dynamics (MD) simulations have investigated the changes in structural and dynamic features of native DNA on binding to the active compound 9. All the synthesized compounds have increased the uptake of Vinblastine in MDR KBChR-8-5 cells to an extent of 1.25- to1.9-fold than standard modulator Verapamil of similar concentration. These findings allowed us to draw preliminary conclusions about the structural features of 2-bromoacridones and further chemical enhancement will improve the binding affinity of the acridone derivatives to CT-DNA for better drug-DNA interaction. The molecular modeling studies have shown mechanism of action and the binding modes of the acridones to DNA.

Functional interplay between the RK motif and linker segment dictates Oct4-DNA recognition

The POU family transcription factor Oct4 plays pivotal roles in regulating pluripotency and somatic cell reprogramming. Previous studies have indicated an important role for major groove contacts in Oct4–DNA recognition; however, the contributions of the RK motif in the POUh domain and the linker segment joining the two DNA-binding domains remain poorly understood. Here, by combining molecular modelling and functional assays, we find that the RK motif is essential for Oct4–DNA association by recognizing the narrowed DNA minor groove. Intriguingly, computational simulations reveal that the function of the RK motif may be finely tuned by H-bond interactions with the partially disordered linker segment and that breaking these interactions significantly enhances the DNA binding and reprogramming activities of Oct4. These findings uncover a self-regulatory mechanism for specific Oct4–DNA recognition and provide insights into the functional crosstalk at the molecular level that may illuminate mechanistic studies of the Oct protein family and possibly transcription factors in the POU family. Our gain-of-function Oct4 mutants might also be useful tools for use in reprogramming and regenerative medicine.

Specific minor groove solvation is a crucial determinant of DNA binding site recognition

The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein's ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.

Sequence-specific DNA interactions with calixarene-based langmuir monolayers

The interactions of an amphiphilic calixarene, namely p-guanidino-dodecyloxy-calix[4]arene, 1, self-assembled as Langmuir monolayers, with short double stranded DNA, were investigated by surface pressure-area (π-A) isotherms, surface ellipsometry and Brewster angle microscopy (BAM). Three DNA 30mers were used as models, poly(AT), poly(GC) and a random DNA sequence with 50% of G:C base pairs. The interactions of these model DNA duplexes with 1-based Langmuir monolayers were studied by measuring compression isotherms using increasing DNA concentrations (10(-6), 10(-5), 10(-4), and 5 × 10(-4) g L(-1)) in the aqueous subphase. The isotherms of 1 showed an expansion of the monolayer with, interestingly, significant differences depending on the duplex DNA sequence studied. Indeed, the interactions of 1-based monolayers with poly(AT) led to an expansion of the monolayer that was significantly more pronounced that for monolayers on subphases of poly(GC) and the random DNA sequence. The structure and thickness of 1-based Langmuir monolayers were investigated by BAM and surface ellipsometry that showed differences in thickness and structure between a monolayer formed on pure water or on a DNA subphase, with here again relevant dissimilarities depending on the DNA composition.

The presence of an RNA:DNA hybrid that is prone to slippage promotes termination by T7 RNA polymerase

Intrinsic termination signals for multisubunit bacterial RNA polymerases (RNAPs) encode a GC-rich stem-loop structure followed by a polyuridine [poly(U)] tract, and it has been proposed that steric clash of the stem-loop with the exit pore of the RNAP imposes a shearing force on the RNA in the downstream RNA:DNA hybrid, resulting in misalignment of the active site. The structurally unrelated T7 RNAP terminates at a similar type of signal (TΦ), suggesting a common mechanism for termination. In the absence of a hairpin (passive conditions), T7 RNAP slips efficiently in both homopolymeric A and U tracts, and we have found that replacement of the U tract in TΦ with a slippage-prone A tract still allows efficient termination. Under passive conditions, incorporation of a single G residue following a poly(U) tract (which is the situation during termination at TΦ) results in a "locked" complex that is unable to extend the transcript. Our results support a model in which transmission of the shearing force generated by steric clash of the hairpin with the exit pore is promoted by the presence of a slippery tracts downstream, resulting in alterations in the active site and the formation of a locked complex that represents an early step in the termination pathway.

Conformational modulation of DNA by polyamide binding: structural effects of f-Im-Py-Im based derivatives on 5'-ACGCGT-3'

The DNA sequence 5'-ACGCGT-3' is in the core site of the Mlu 1 cell-cycle box, a transcriptional element in the promoter region of human Dbf4 gene that is highly correlated with a large number of aggressive solid cancers. The polyamide formamido-imidazole-pyrrole-imidazole-amine(+) (f-Im-Py-Im-Am(+) ) can target the minor groove of 5'-ACGCGT-3' as an antiparallel stacked dimer and has shown good activity in inhibiting transcription factor binding. Recently, f-Im-Py-Im-Am(+) derivatives that involve different orthogonally positioned substituents were synthesized to target the same binding site, and some of them have displayed improved binding and pharmacological properties. In this study, the gel electrophoresis-ligation ladders assay was used to evaluate the conformational effects of f-Im-Py-Im-Am(+) and derivatives on the target DNA, an essential factor for establishing the molecular basis of polyamide-DNA complexes and their transcription factor inhibition. The results show that the ACGCGT site in DNA has a relatively wide minor groove and a B-form like overall structure. After binding with f-Im-Py-Im-Am(+) derivatives, the DNA conformation is changed as indicated by the different mobilities in the gel. These conformational effects on DNA will at least help to point to the mechanism for the observed Mlu 1 inhibition activity of these polyamides. Therefore, modulating DNA transcription by locking the DNA shape or altering the minor groove geometry to affect the binding affinity of certain transcription factors is an attractive possible therapeutic mechanism for polyamides. Some of the substituents are charged with electrostatic interactions with DNA phosphate groups, and their charge effects on DNA gel mobility have been observed.

Kinetics of DNA duplex formation: A-tracts versus AT-tracts

A-tracts (AAAA…:TTTT…) form much faster than AT-tracks (ATAT…:TATA…).

Molecular dynamics simulations of DNA-free and DNA-bound TAL effectors

TAL (transcriptional activator-like) effectors (TALEs) are DNA-binding proteins, containing a modular central domain that recognizes specific DNA sequences. Recently, the crystallographic studies of TALEs revealed the structure of DNA-recognition domain. In this article, molecular dynamics (MD) simulations are employed to study two crystal structures of an 11.5-repeat TALE, in the presence and absence of DNA, respectively. The simulated results indicate that the specific binding of RVDs (repeat-variable diresidues) with DNA leads to the markedly reduced fluctuations of tandem repeats, especially at the two ends. In the DNA-bound TALE system, the base-specific interaction is formed mainly by the residue at position 13 within a TAL repeat. Tandem repeats with weak RVDs are unfavorable for the TALE-DNA binding. These observations are consistent with experimental studies. By using principal component analysis (PCA), the dominant motions are open-close movements between the two ends of the superhelical structure in both DNA-free and DNA-bound TALE systems. The open-close movements are found to be critical for the recognition and binding of TALE-DNA based on the analysis of free energy landscape (FEL). The conformational analysis of DNA indicates that the 5′ end of DNA target sequence has more remarkable structural deformability than the other sites. Meanwhile, the conformational change of DNA is likely associated with the specific interaction of TALE-DNA. We further suggest that the arrangement of N-terminal repeats with strong RVDs may help in the design of efficient TALEs. This study provides some new insights into the understanding of the TALE-DNA recognition mechanism.

Role of minor groove width and hydration pattern on amsacrine interaction with DNA

Amsacrine is an anilinoacridine derivative anticancer drug, used to treat a wide variety of malignancies. In cells, amsacrine poisons topoisomerase 2 by stabilizing DNA-drug-enzyme ternary complex. Presence of amsacrine increases the steady-state concentration of these ternary complexes which in turn hampers DNA replication and results in subsequent cell death. Due to reversible binding and rapid slip-out of amsacrine from DNA duplex, structural data is not available on amsacrine-DNA complexes. In the present work, we designed five oligonucleotide duplexes, differing in their minor groove widths and hydration pattern, and examined their binding with amsacrine using attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy. Complexes of amsacrine with calf thymus DNA were also evaluated for a comparison. Our results demonstrate for the first time that amsacrine is not a simple intercalator; rather mixed type of DNA binding (intercalation and minor groove) takes place between amsacrine and DNA. Further, this binding is highly sensitive towards the geometries and hydration patterns of different minor grooves present in the DNA. This study shows that ligand binding to DNA could be very sensitive to DNA base composition and DNA groove structures. Results demonstrated here could have implication for understanding cytotoxic mechanism of aminoacridine based anticancer drugs and provide directions to modify these drugs for better efficacy and few side-effects.

DNA A-tracts are not curved in solutions containing high concentrations of monovalent cations

The intrinsic curvature of seven 98 bp DNA molecules containing up to four centrally located A6-tracts has been measured by gel and capillary electrophoresis as a function of the number and arrangement of the A-tracts. At low cation concentrations, the electrophoretic mobility observed in polyacrylamide gels and in free solution decreases progressively with the increasing number of phased A-tracts, as expected for DNA molecules with increasingly curved backbone structures. Anomalously slow electrophoretic mobilities are also observed for DNA molecules containing two pairs of phased A-tracts that are out of phase with each other, suggesting that out-of-phase distortions of the helix backbone do not cancel each other out. The mobility decreases observed for the A-tract samples are due to curvature, not cation binding in the A-tract minor groove, because identical free solution mobilities are observed for a molecule with four out-of-phase A-tracts and one with no A-tracts. Surprisingly, the curvature of DNA A-tracts is gradually lost when the monovalent cation concentration is increased to ∼200 mM, regardless of whether the cation is a hydrophilic ion like Na+, NH4+, or Tris+ or a hydrophobic ion like tetrabutylammonium. The decrease in A-tract curvature with increasing ionic strength, along with the known decrease in A-tract curvature with increasing temperature, suggests that DNA A-tracts are not significantly curved under physiological conditions.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

CYTOTOXICITY OF POLYAMINE ANALOGS IS DIRECTLY RELATED TO THEIR DNA AFFINITY AS DETERMINED BY POLYACRYLAMIDE GEL COELECTROPHORESIS (PACE) METHOD

Polyamines are essential for cell growth. Polyamine analogs that can replace intracellular polyamines inhibit tumor cell proliferation both in culture as well as in animal models. The positively charged polyamines interact with the negatively charged DNA backbone both in a nonspecific manner, as well as sequence specifically through direct or water mediated hydrogen bonds. Therefore, it is difficult to ascertain the exact interactions that regulate the biological functions of polyamines. Several attempts have been made to determine the thermodynamic parameters of polyamine-DNA interactions with conflicting results. Here, we report a simple method of determining the apparent association constants for polyamine-DNA interaction by using polyacrylamide gel coelectrophoresis (PACE). We have used several cytotoxic polyamine analogs of different conformations and chain lengths. We observed that polyamine analogs with higher charge density or with conformational restrictions, which are absent in the naturally occurring polyamines, interact with DNA more strongly than do natural polyamines. A comparison of the cytotoxicities of the polyamine analogs against human tumor cell lines with their DNA affinities revealed that the higher the DNA affinity the more the cytoxicity of the analogs. The direct correlation between DNA affinities and cytotoxities provides a novel method for a rational design of therapeutically effective cytotoxic polyamine analogs.

Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Structural basis for recognition of AT-rich DNA by unrelated xenogeneic silencing proteins

H-NS and Lsr2 are nucleoid-associated proteins from Gram-negative bacteria and Mycobacteria, respectively, that play an important role in the silencing of horizontally acquired foreign DNA that is more AT-rich than the resident genome. Despite the fact that Lsr2 and H-NS proteins are dissimilar in sequence and structure, they serve apparently similar functions and can functionally complement one another. The mechanism by which these xenogeneic silencers selectively target AT-rich DNA has been enigmatic. We performed high-resolution protein binding microarray analysis to simultaneously assess the binding preference of H-NS and Lsr2 for all possible 8-base sequences. Concurrently, we performed a detailed structure-function relationship analysis of their C-terminal DNA binding domains by NMR. Unexpectedly, we found that H-NS and Lsr2 use a common DNA binding mechanism where a short loop containing a "Q/RGR" motif selectively interacts with the DNA minor groove, where the highest affinity is for AT-rich sequences that lack A-tracts. Mutations of the Q/RGR motif abolished DNA binding activity. Netropsin, a DNA minor groove-binding molecule effectively outcompeted H-NS and Lsr2 for binding to AT-rich sequences. These results provide a unified molecular mechanism to explain findings related to xenogeneic silencing proteins, including their lack of apparent sequence specificity but preference for AT-rich sequences. Our findings also suggest that structural information contained within the DNA minor groove is deciphered by xenogeneic silencing proteins to distinguish genetic material that is self from nonself.

Understanding the sequence-dependence of DNA groove dimensions: implications for DNA interactions

Background

The B-DNA major and minor groove dimensions are crucial for DNA-protein interactions. It has long been thought that the groove dimensions depend on the DNA sequence, however this relationship has remained elusive. Here, our aim is to elucidate how the DNA sequence intrinsically shapes the grooves.

Methodology/Principal Findings

The present study is based on the analysis of datasets of free and protein-bound DNA crystal structures, and from a compilation of NMR ³¹P chemical shifts measured on free DNA in solution on a broad range of representative sequences. The ³¹P chemical shifts can be interpreted in terms of the BI↔BII backbone conformations and dynamics. The grooves width and depth of free and protein-bound DNA are found to be clearly related to the BI/BII backbone conformational states. The DNA propensity to undergo BI↔BII backbone transitions is highly sequence-dependent and can be quantified at the dinucleotide level. This dual relationship, between DNA sequence and backbone behavior on one hand, and backbone behavior and groove dimensions on the other hand, allows to decipher the link between DNA sequence and groove dimensions. It also firmly establishes that proteins take advantage of the intrinsic DNA groove properties.

Conclusions/Significance

The study provides a general framework explaining how the DNA sequence shapes the groove dimensions in free and protein-bound DNA, with far-reaching implications for DNA-protein indirect readout in both specific and non specific interactions.

Vacuum-ultraviolet circular dichroism spectroscopy of DNA: a valuable tool to elucidate topology and electronic coupling in DNA

Circular dichroism (CD) is a powerful technique to obtain information on electronic transitions and has been used extensively for studies on DNA. Most experiments are done in the UV region but new information is often revealed from extending the wavelength region down into the vacuum ultraviolet (VUV) region. Such experiments are most easily carried out with synchrotron radiation (SR) light sources that provide large photon fluxes. Here we provide a summary of the SRCD data taken on different DNA strands with emphasis on results from our own laboratory within the last five years.(1-3) Signal intensities in the VUV are often significantly larger than those in the UV, and the electronic coupling between bases may increase with excitation energy. CD spectroscopy is particularly useful for investigating the extent of electronic coupling within a strand, i.e., the degree of delocalisation of the excited-state electronic wavefunction. The spatial extent of the wavefunction may be limited to just one base or it extends over two or more bases in a stack or between bases on different strands.(4,5) The actual character of the electronically excited state is linked to base composition and sequence as well as DNA folding motif (A-, B-, Z-DNA, triplexes, quadruplexes, etc.). The latter depends on experimental conditions such as solution acidity, temperature, ionic strength, and solvent.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

Fingerprints of bonding motifs in DNA duplexes of adenine and thymine revealed from circular dichroism: synchrotron radiation experiments and TDDFT calculations

Synchrotron radiation circular dichroism (SRCD) spectra were recorded for a family of 12 DNA duplexes that all contain nine adenines (A) and nine thymines (T) in each strand but in different combinations. The total number of AT Watson-Crick (WC) base pairs is constant (18), but the number of cross-strand (CS) hydrogen bonds between A and T varies between 0 and 16, the maximum possible. Eleven of the duplexes have one or more A tracts, and one duplex has T tracts. The signals due to hybridization were found from subtraction of spectra of single strands from spectra of the duplexes. The residual spectrum of the T-tract duplex T(9)A(9):A(9)T(9) (5'-3':3'-5') significantly differs from that of the A-tract duplex A(9)T(9):T(9)A(9), but only below 210 nm, which suggests that the signal in this region depends on the superhelicity of the duplex. A principal component analysis of all residual spectra reveals that spectra of A-tract duplexes can be obtained to a good approximation as a linear combination of just two basis spectra. The first component is assigned to the spectrum of 18 WC and 8 CS pairs, whereas the second component is that of 8 CS pairs. This interpretation is supported by separate experiments on duplexes of varying lengths but with similar arrangements of the A and T's and by experiments on two other duplex families of 14 and 30 base pairs. The best correlation is obtained by the assumption that cross-strand interactions occur as long as there are two adenine neighbors in a strand. Our data indicate that a circular dichroism spectrum of a duplex containing only A and T can simply be inferred from the number of WC base pairs and the number of CS interactions, and we provide reference spectra for these two interactions. Finally, time dependent density functional theory calculations of the circular dichroism spectra for an isolated WC base pair and two different CS base pairs (between adenine N-6 amine and thymine O-4 or between adenine C-2-H and thymine O-2) were performed to provide some additional support for the interpretation of the experimental spectra. We find large differences between the two calculated CS spectra. However, there is a reasonable qualitative agreement between the calculated WC and the C-2-H...O-2 CS spectra and those deduced from the experimental data.

The unique structure of A-tracts and intrinsic DNA bending

Short runs of adenines are a ubiquitous DNA element in regulatory regions of many organisms. When runs of 4–6 adenine base pairs (‘A-tracts’) are repeated with the helical periodicity, they give rise to global curvature of the DNA double helix, which can be macroscopically characterized by anomalously slow migration on polyacrylamide gels. The molecular structure of these DNA tracts is unusual and distinct from that of canonical B-DNA. We review here our current knowledge about the molecular details of A-tract structure and its interaction with sequences flanking them of either side and with the environment. Various molecular models were proposed to describe A-tract structure and how it causes global deflection of the DNA helical axis. We review old and recent findings that enable us to amalgamate the various findings to one model that conforms to the experimental data. Sequences containing phased repeats of A-tracts have from the very beginning been synonymous with global intrinsic DNA bending. In this review, we show that very often it is the unique structure of A-tracts that is at the basis of their widespread occurrence in regulatory regions of many organisms. Thus, the biological importance of A-tracts may often be residing in their distinct structure rather than in the global curvature that they induce on sequences containing them.

Crystallographic analysis of a sex-specific enhancer element: sequence-dependent DNA structure, hydration, and dynamics

The crystal structure of a sex-specific enhancer element is described at a resolution of 1.6 A. This 16-bp site, designated Dsx(A), functions in the regulation of a genetic switch between male and female patterns of gene expression in Drosophila melanogaster. Related sites are broadly conserved in metazoans, including in the human genome. This enhancer element is unusually rich in general regulatory sequences related to DNA recognition by multiple classes of eukaryotic transcription factors, including the DM motifs, homeodomain, and high mobility group box. Whereas free DNA is often crystallized as an A-form double helix, Dsx(A) was crystallized as B-DNA and thus provides a model for the prebound conformation of diverse regulatory DNA complexes. Sequence-dependent conformational properties that extend features of shorter B-DNA fragments with respect to double helical parameters, groove widths, hydration, and binding of divalent metal ions are observed. The structure also exhibits a sequence-dependent pattern of isotropic thermal B-factors, suggesting possible variation in the local flexibility of the DNA backbone. Such fluctuations are in accord with structural variability observed in prior B-DNA structures. We speculate that sites of intrinsic flexibility within a DNA control element provide hinges for its protein-directed reorganization in a transcriptional preinitiation complex.

Probing the interactions of the solvated electron with DNA by molecular dynamics simulations: bromodeoxyuridine substituted DNA

Solvated electrons ((e-)(aq)) are produced during water radiolysis and can interact with biological substrates, including DNA. To augment DNA damage, radiosensitizers such as bromo-deoxyuridine (BUdR), often referred to as an "electron affinic radiosensitizer", are incorporated in place of isosteric thymidine. However, little is known about the primary interactions of (e-)(aq) with DNA. In the present study we addressed this problem by applying molecular modeling and molecular dynamics (MD) simulations to a system of normal (BUdR.A)-DNA and a hydrated electron, where the excess electron was modeled as a localized (e-)(H2O)6 anionic cluster. Our goals were to evaluate the suitability of the MD simulations for this application; to characterize the motion of (e-)(aq) around DNA (e.g., diffusion coefficients); to identify and describe configurational states of close (e-)(aq) localization to DNA; and to evaluate the structural dynamics of DNA in the presence of (e-)(aq). The results indicate that (e-)(aq) has distinct space-preferences for forming close contacts with DNA and is more likely to interact directly with nucleotides other than BUdR. Several classes of DNA - (e-)(aq) contact sites, all within the major groove, were distinguished depending on the structure of the intermediate water layer H-bonding pattern (or its absence, i.e., a direct H-bonding of (e-)(aq) with DNA bases). Large-scale structural perturbations were identified during and after the (e-)(aq) approached the DNA from the major groove side, coupled with deeper penetration of sodium counterions in the minor groove.

Complex microsatellite dynamics in the myostatin gene within ruminants

A microsatellite has previously been identified in myostatin in cattle. Sequencing of this region from other artiodactyls coupled with phylogenetic analysis has been used to uncover the potential origins of the microsatellite event, which appears either to have been born twice or to have been gained and lost within ruminants. While caprids and ovids share the ancestral state with pigs and other mammals, microsatellite activity (length polymorphism) is uncovered in both deer and bovids. The dynamic process of microsatellite evolution, including birth, is discussed here in light of several models. Finally, these models are evaluated in the context of patterns of microsatellite conservation between closely related mammalian genomes.

DNA A-tracts bending: polarization effects on electrostatic interactions across their minor groove

Bending by the DNA A-tracts constitutes a contentious issue, suggesting deficiencies in the physics employed so far. Here, we inquire as to the importance in this bending of many-body polarization effects on the electrostatic interactions across their narrow minor groove. We have done this on the basis of the findings of Jarque and Buckingham who developed a procedure based on a Monte Carlo simulation for two charges of the same sign embedded in a polarizable medium. Remarkably, the present analysis reveals that for compact DNA conformations, which result from dynamic effects, an overall attractive interaction operates between the phosphate charges; this interaction is especially strong for the narrow minor groove of the A-tracts, suggesting a tendency for DNA to bend toward this groove. This tendency is in agreement with the conclusions of electrophoretic and NMR solution studies. The present analysis is also consistent with the experimental observations that the minor groove is much more easily compressible than the major groove and the bending propensity of the A-tracts is greatly reduced at "premelting" temperatures. By contrast, the dielectric screening model predicts a repulsion between the phosphate charges and is not consistent with the aforementioned bending tendency or experimental observations.

Principles of protein-DNA recognition revealed in the structural analysis of Ndt80-MSE DNA complexes

The Saccharomyces cerevisiae transcription factor Ndt80 selectively binds a DNA consensus sequence (the middle sporulation element [MSE]) to activate gene expression after the successful completion of meiotic recombination. Here we report the X-ray crystal structures of Ndt80 bound to ten distinct MSE variants. Comparison of these structures with the structure of Ndt80 bound to a consensus MSE reveals structural principles that determine the DNA binding specificity of this transcription factor. The 5' GC-rich end of the MSE contains distinct 5'-YpG-3' steps that are recognized by arginine side chains through a combination of hydrogen bonding and cation-pi interactions. The 3' AT-rich region is recognized via minor groove contacts that sterically exclude the N2 atom of GC base pairs. The conformation of the AT-rich region is fixed by interactions with the protein that favor recognition of poly(A)-poly(T) versus mixed AT sequences through an avoidance of major groove steric clashes at 5'-ApT-3' steps.