DNA bending by an adenine--thymine tract and its role in gene regulation

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

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

Discovering DNA shape motifs with multiple DNA shape features: generalization, methods, and validation

DNA motifs are crucial patterns in gene regulation. DNA-binding proteins (DBPs), including transcription factors, can bind to specific DNA motifs to regulate gene expression and other cellular activities. Past studies suggest that DNA shape features could be subtly involved in DNA-DBP interactions. Therefore, the shape motif annotations based on intrinsic DNA topology can deepen the understanding of DNA-DBP binding. Nevertheless, high-throughput tools for DNA shape motif discovery that incorporate multiple features altogether remain insufficient. To address it, we propose a series of methods to discover non-redundant DNA shape motifs with the generalization to multiple motifs in multiple shape features. Specifically, an existing Gibbs sampling method is generalized to multiple DNA motif discovery with multiple shape features. Meanwhile, an expectation-maximization (EM) method and a hybrid method coupling EM with Gibbs sampling are proposed and developed with promising performance, convergence capability, and efficiency. The discovered DNA shape motif instances reveal insights into low-signal ChIP-seq peak summits, complementing the existing sequence motif discovery works. Additionally, our modelling captures the potential interplays across multiple DNA shape features. We provide a valuable platform of tools for DNA shape motif discovery. An R package is built for open accessibility and long-lasting impact: https://zenodo.org/doi/10.5281/zenodo.10558980.

Probing the conformational dynamics of an Ago-RNA complex in water/methanol solution

Argonaute (Ago) proteins mediate target recognition guiding miRNA to bind complementary mRNA primarily in the seed region. However, additional pairing can occur beyond the seed, forming a supplementary duplex that can contribute to the guide-target affinity. In order to shed light on the connection, between protein-RNA interactions and miRNA-mRNA seed and supplementary duplex mobility, we carried out molecular dynamics simulations at the microsecond time-scale using a different approach compared to the ones normally used. Until now, theoretical investigations with classical MD on Ago-RNA complexes have been focused primarily on pure water solvent, which mimics the natural environment of biological molecules. Here, we explored the conformational space of a human Ago2 (hAgo2) bound to the seed + supplementary miRNA-mRNA duplex, using the solvent environment as a molecular probe. MD simulations have been performed in a mixture of water/MeOH at a molar ratio of 70 : 30 as well as in pure water for comparison. Our findings revealed that the mixed solvent promotes protein RNA association, principally enhancing salt-linkages between basic amino acid side-chains and acidic phosphates of the sugar-phosphate backbone. The primary effect registered was the restriction of supplementary duplex flexibility and the stabilization of the miRNA 3' terminus. Interestingly, we observed that the influence of the solvent appears to have almost no impact on the conformation of the seed duplex.

Molecular basis for the functional roles of the multimorphic T95R mutation of IRF4 causing human autosomal dominant combined immunodeficiency

Interferon regulatory factor 4 (IRF4) is a transcription factor that regulates the development and function of immune cells. Recently, a new multimorphic mutation T95R was identified in the IRF4 DNA-binding domain (DBD) in patients with autosomal dominant combined immune deficiency. Here, we characterized the interactions of the wild-type IRF4-DBD (IRF4-DBDWT) and T95R mutant (IRF4-DBDT95R) with a canonical DNA sequence and several noncanonical DNA sequences. We found that compared to IRF4-DBDWT, IRF4-DBDT95R exhibits higher binding affinities for both canonical and noncanonical DNAs, with the highest preference for the noncanonical GATA sequence. The crystal structures of IRF4-DBDWT in complex with the GATA sequence and IRF4-DBDT95R in complexes with both canonical and noncanonical DNAs were determined, showing that the T95R mutation enhances the interactions of IRF4-DBDT95R with the canonical and noncanonical DNAs to achieve higher affinity and specificity. Collectively, our data provide the molecular basis for the gain-of-function and new function of IRF4T95R.

Structural insights into the DNA recognition mechanism by the bacterial transcription factor PdxR

Specificity in protein-DNA recognition arises from the synergy of several factors that stem from the structural and chemical signatures encoded within the targeted DNA molecule. Here, we deciphered the nature of the interactions driving DNA recognition and binding by the bacterial transcription factor PdxR, a member of the MocR family responsible for the regulation of pyridoxal 5'-phosphate (PLP) biosynthesis. Single particle cryo-EM performed on the PLP-PdxR bound to its target DNA enabled the isolation of three conformers of the complex, which may be considered as snapshots of the binding process. Moreover, the resolution of an apo-PdxR crystallographic structure provided a detailed description of the transition of the effector domain to the holo-PdxR form triggered by the binding of the PLP effector molecule. Binding analyses of mutated DNA sequences using both wild type and PdxR variants revealed a central role of electrostatic interactions and of the intrinsic asymmetric bending of the DNA in allosterically guiding the holo-PdxR-DNA recognition process, from the first encounter through the fully bound state. Our results detail the structure and dynamics of the PdxR-DNA complex, clarifying the mechanism governing the DNA-binding mode of the holo-PdxR and the regulation features of the MocR family of transcription factors.

Z-flipon variants reveal the many roles of Z-DNA and Z-RNA in health and disease

Identifying roles for Z-DNA remains challenging given their dynamic nature. Here, we perform genome-wide interrogation with the DNABERT transformer algorithm trained on experimentally identified Z-DNA forming sequences (Z-flipons). The algorithm yields large performance enhancements (F1 = 0.83) over existing approaches and implements computational mutagenesis to assess the effects of base substitution on Z-DNA formation. We show Z-flipons are enriched in promoters and telomeres, overlapping quantitative trait loci for RNA expression, RNA editing, splicing, and disease-associated variants. We cross-validate across a number of orthogonal databases and define BZ junction motifs. Surprisingly, many effects we delineate are likely mediated through Z-RNA formation. A shared Z-RNA motif is identified in SCARF2, SMAD1, and CACNA1 transcripts, whereas other motifs are present in noncoding RNAs. We provide evidence for a Z-RNA fold that promotes adaptive immunity through alternative splicing of KRAB domain zinc finger proteins. An analysis of OMIM and presumptive gnomAD loss-of-function datasets reveals an overlap of Z-flipons with disease-causing variants in 8.6% and 2.9% of Mendelian disease genes, respectively, greatly extending the range of phenotypes mapped to Z-flipons.

Thermotoga maritima oriC involves a DNA unwinding element with distinct modules and a DnaA-oligomerizing region with a novel directional binding mode

Initiation of chromosomal replication requires dynamic nucleoprotein complexes. In most eubacteria, the origin oriC contains multiple DnaA box sequences to which the ubiquitous DnaA initiators bind. In Escherichia coli oriC, DnaA boxes sustain construction of higher-order complexes via DnaA-DnaA interactions, promoting the unwinding of the DNA unwinding element (DUE) within oriC and concomitantly binding the single-stranded DUE to install replication machinery. Despite the significant sequence homologies among DnaA proteins, bacterial oriC sequences are highly diverse. The present study investigated the design of oriC (tma-oriC) from Thermotoga maritima, an evolutionarily ancient eubacterium. The minimal tma-oriC sequence includes a DUE and a flanking region containing five DnaA boxes recognized by the cognate DnaA initiator (tmaDnaA). This DUE was comprised of two distinct functional modules, an unwinding module and a tmaDnaA-binding module. Three direct repeats of the trinucleotide TAG within DUE were essential for both unwinding and single-stranded DUE binding by tmaDnaA complexes constructed on the DnaA boxes. Its surrounding AT-rich sequences stimulated only duplex unwinding. Moreover, head-to-tail oligomers of ATP-bound tmaDnaA were constructed within tma-oriC, irrespective of the directions of the DnaA boxes. This binding mode was considered to be induced by flexible swiveling of DnaA domains III and IV, which were responsible for DnaA-DnaA interactions and DnaA box binding, respectively. Phasing of specific tmaDnaA boxes in tma-oriC DNA was also responsible for unwinding. These findings indicate that a single-stranded DUE recruitment mechanism was responsible for unwinding, and would enhance understanding of the fundamental molecular nature of the origin sequences present in evolutionarily divergent bacteria.

Structural insights into DNA recognition by the BEN domain of the transcription factor BANP

The BEN domain-containing transcription factors regulate transcription by recruiting chromatin-modifying factors to specific chromatin regions via their DNA-binding BEN domains. The BEN domain of BANP has been shown to bind to a CGCG DNA sequence or an AAA-containing MARs (matrix attachment regions) DNA sequence. Consistent with these in vivo observations, we identified an optimal DNA binding sequence of AAATCTCG by PBM (protein binding microarray), which was also confirmed by our ITC (Isothermal Titration Calorimetry) and mutagenesis results to uncover additional mechanistic details about DNA binding by the BEN domain of BANP. We then determined crystal structures of the BANP BEN domain in apo form and in complex with a CGCG-containing DNA, respectively, which revealed that the BANP BEN domain mainly used the electrostatic interactions to bind DNA with some base-specific interactions with the TC motifs. Our ITC results also showed that BANP bound to unmethylated and methylated DNAs with comparable binding affinities. Our complex structure of BANP-mCGCG revealed that the BANP BEN domain bound to the unmethylated and methylated DNAs in a similar mode and cytosine methylation did not get involved in binding, which is also consistent with our observations from the complex structures of the BEND6 BEN domain with the CGCG or CGmCG DNAs. Taken together, our results further elucidate the elements important for DNA recognition and transcriptional regulation by the BANP BEN domain-containing transcription factor.

Structural determinants of DNA recognition by the NO sensor NsrR and related Rrf2-type [FeS]-transcription factors

Several transcription factors of the Rrf2 family use an iron-sulfur cluster to regulate DNA binding through effectors such as nitric oxide (NO), cellular redox status and iron levels. [4Fe-4S]-NsrR from Streptomyces coelicolor (ScNsrR) modulates expression of three different genes via reaction and complex formation with variable amounts of NO, which results in detoxification of this gas. Here, we report the crystal structure of ScNsrR complexed with an hmpA1 gene operator fragment and compare it with those previously reported for [2Fe-2S]-RsrR/rsrR and apo-IscR/hyA complexes. Important structural differences reside in the variation of the DNA minor and major groove widths. In addition, different DNA curvatures and different interactions with the protein sensors are observed. We also report studies of NsrR binding to four hmpA1 variants, which indicate that flexibility in the central region is not a key binding determinant. Our study explores the promotor binding specificities of three closely related transcriptional regulators.

The crystal structure of the iron-sulfur protein NsrR from Streptomyces coelicolor bound to a gene operator fragment is reported and compared with other structures, giving insight into the structural determinants of DNA recognition by the NO sensor.

Computational Investigation of Bending Properties of RNA AUUCU, CCUG, CAG, and CUG Repeat Expansions Associated With Neuromuscular Disorders

Expansions of RNA AUUCU, CCUG, CAG, and CUG repeats cause spinocerebellar ataxia type 10, myotonic dystrophy type 2, Huntington’s disease, and myotonic dystrophy type 1, respectively. By performing extensive molecular dynamic simulations, we investigated the bending propensities and conformational landscapes adopted by 3×3, 2×2, and 1×1 internal loops observed in RNA AUUCU, CCUG, CAG, and CUG repeat expansions using model systems having biologically relevant repeat sizes. We show that the conformational variability experienced by these loops is more complex than previous reports where a variety of unconventional hydrogen bonds are formed. At the global scale, strong bending propensity was observed in r(AUUCU)₁₀, r(CCUG)₁₅, r(CAG)_20, and r(CUG)_20, and, to a lesser extent, in r(AUUCU)₄, r(CCUG)₁₀, r(CAG)₁₀, and r(CUG)₁₀. Furthermore, RNA CAG repeats exhibit a tendency toward bent states with more than 50% of observed conformations having bending angles greater than 50°, while RNA CUG repeats display relatively linear-like conformations with extremely bent conformations accounting for less than 25% of the observed structures. Conformations experienced by RNA AUUCU repeats are a combination of strongly bent and kinked structures. The bent states in RNA CCUG repeats mostly fall into the moderately bent category with a marginal ensemble experiencing extreme bending. The general pattern observed in all the bent structures indicates the collapse of the major groove width as the mechanical trigger for bending, which is caused by alteration of base pair step parameters at multiple locations along the RNA due to local distortions at the loop sites. Overextension is also observed in all the RNA repeats that is attributed to widening of the major groove width as well as undertwisting phenomenon. This information and the rich structural repository could be applied for structure based small molecule design targeting disease-causing RNAs. The bending propensities of these constructs, at the global level, could also have implications on how expanded RNA repeats interact with proteins.

Orientation Dependence of DNA Blunt-End Stacking Studied by Free-Energy Simulations

DNA blunt ends can associate mediated by stacking interactions between the terminal base pairs that form blunt ends. The blunt end association plays a role in DNA repair and recombination processes and can also be of importance for the design of DNA-based nano-materials. Its function depends on the sequence and on the geometric arrangement that leads to stable interaction. For a stacked state, the relative orientation (twisting) of the base pairs is important. Molecular dynamics and advanced sampling simulations were used to calculate free energy change associated with twist changes of the stacked blunt-end base pairs. The calculations reproduce blunt stacking arrangements found in crystal structures of DNA oligonucleotides as free energy minima. To elucidate the physical origin of the stabilization of certain angular arrangements, the interactions between backbone atoms in the blunt-end stack were switched off in additional free energy calculations. It allows us to decipher the contributions to stacking stabilization due to the nucleobases and the backbone and to analyze the sequence dependence of the angular stacking preferences. Good qualitative agreement was also found for the comparison with quantum mechanical calculations. The results may help in the design of novel DNA-based materials.

Nonalternating purine pyrimidine sequences can form stable left-handed DNA duplex by strong topological constraint

In vivo, left-handed DNA duplex (usually refers to Z-DNA) is mainly formed in the region of DNA with alternating purine pyrimidine (APP) sequence and plays significant biological roles. It is well known that d(CG)_n sequence can form Z-DNA most easily under negative supercoil conditions, but its essence has not been well clarified. The study on sequence dependence of Z-DNA stability is very difficult without modification or inducers. Here, by the strong topological constraint caused by hybridization of two complementary short circular ssDNAs, left-handed duplex part was generated for various sequences, and their characteristics were investigated by using gel-shift after binding to specific proteins, CD and T_m analysis, and restriction enzyme cleavage. Under the strong topological constraint, non-APP sequences can also form left-handed DNA duplex as stable as that of APP sequences. As compared with non-APP sequences, the thermal stability difference for APP sequences between Z-form and B-form is smaller, which may be the reason that Z-DNA forms preferentially for APP ones. This result can help us to understand why nature selected APP sequences to regulate gene expression by transient Z-DNA formation, as well as why polymer with chirality can usually form both duplexes with left- or right-handed helix.

DNA sequence-dependent positioning of the linker histone in a nucleosome: A single-pair FRET study

Linker histones (LHs) bind to nucleosomes with their globular domain (gH) positioned in either an on- or an off-dyad binding mode. Here, we study the effect of the linker DNA (L-DNA) sequence on the binding of a full-length LH, Xenopus laevis H1.0b, to a Widom 601 nucleosome core particle (NCP) flanked by two 40 bp long L-DNA arms, by single-pair FRET spectroscopy. We varied the sequence of the 11 bp of L-DNA adjoining the NCP on either side, making the sequence either A-tract, purely GC, or mixed with 64% AT. The labeled gH consistently exhibited higher FRET efficiency with the labeled L-DNA containing the A-tract than that with the pure-GC stretch, even when the stretches were swapped. However, it did not exhibit higher FRET efficiency with the L-DNA containing 64% AT-rich mixed DNA when compared to the pure-GC stretch. We explain our observations with a model that shows that the gH binds on dyad and that two arginines mediate recognition of the A-tract via its characteristically narrow minor groove. To investigate whether this on-dyad minor groove-based recognition was distinct from previously identified off-dyad major groove-based recognition, a nucleosome was designed with A-tracts on both the L-DNA arms. One A-tract was complementary to thymine and the other to deoxyuridine. The major groove of the thymine-tract was lined with methyl groups that were absent from the major groove of the deoxyuridine tract. The gH exhibited similar FRET for both these A-tracts, suggesting that it does not interact with the thymine methyl groups exposed on the major groove. Our observations thus complement previous studies that suggest that different LH isoforms may employ different ways of recognizing AT-rich DNA and A-tracts. This adaptability may enable the LH to universally compact scaffold-associated regions and constitutive heterochromatin, which are rich in such 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.

Allostery through DNA drives phenotype switching

Allostery is a pervasive principle to regulate protein function. Growing evidence suggests that also DNA is capable of transmitting allosteric signals. Yet, whether and how DNA-mediated allostery plays a regulatory role in gene expression remained unclear. Here, we show that DNA indeed transmits allosteric signals over long distances to boost the binding cooperativity of transcription factors. Phenotype switching in Bacillus subtilis requires an all-or-none promoter binding of multiple ComK proteins. We use single-molecule FRET to demonstrate that ComK-binding at one promoter site increases affinity at a distant site. Cryo-EM structures of the complex between ComK and its promoter demonstrate that this coupling is due to mechanical forces that alter DNA curvature. Modifications of the spacer between sites tune cooperativity and show how to control allostery, which allows a fine-tuning of the dynamic properties of genetic circuits.

Most insights on DNA-mediated allostery upon transcription factor (TF) binding were either based on artificial promoters or found to be short-ranged. Here authors use single-molecule FRET and cryo-EM to show that Bacillus subtilis bacteria utilize long-range allostery in a stochastic and reversible phenotype switch.

Molecular Mechanism of the Intercalation of the SOX-4 Protein into DNA Inducing Bends and Kinks

DNA-protein interactions regulate several biophysical functions, yet the mechanism of only a few is investigated in molecular detail. An important example is the intercalation of transcription factor proteins into DNA that produce bent and kinked DNA. Here, we have studied the molecular mechanism of the intercalation of a transcription factor SOX4 into DNA with a goal to understand the sequence of molecular events that precede the bending and kinking of the DNA. Our long well-tempered metadynamics and molecular dynamics (MD) simulations show that the protein primarily binds to the backbone of DNA and rotates around it to form an intercalative native state. We show that although there are multiple pathways for intercalation, the deintercalation pathway matches with the most probable intercalation pathway. In both cases, bending and kinking happen simultaneously, driven by the onset of the intercalation of the amino acid.

Beyond the double helix: DNA structural diversity and the PDB

The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Novel organization of mitochondrial minicircles and guide RNAs in the zoonotic pathogen Trypanosoma lewisi

Kinetoplastid flagellates are known for several unusual features, one of which is their complex mitochondrial genome, known as kinetoplast (k) DNA, composed of mutually catenated maxi- and minicircles. Trypanosoma lewisi is a member of the Stercorarian group of trypanosomes which is, based on human infections and experimental data, now considered a zoonotic pathogen. By assembling a total of 58 minicircle classes, which fall into two distinct categories, we describe a novel type of kDNA organization in T. lewisi. RNA-seq approaches allowed us to map the details of uridine insertion and deletion editing events upon the kDNA transcriptome. Moreover, sequencing of small RNA molecules enabled the identification of 169 unique guide (g) RNA genes, with two differently organized minicircle categories both encoding essential gRNAs. The unprecedented organization of minicircles and gRNAs in T. lewisi broadens our knowledge of the structure and expression of the mitochondrial genomes of these human and animal pathogens. Finally, a scenario describing the evolution of minicircles is presented.

Induced DNA bending by unique dimerization of HigA antitoxin

The bacterial toxin-antitoxin (TA) system regulates cell growth under various environmental stresses. Mycobacterium tuberculosis, the causative pathogen of tuberculosis (TB), has three HigBA type II TA systems with reverse gene organization, consisting of the toxin protein HigB and labile antitoxin protein HigA. Most type II TA modules are transcriptionally autoregulated by the antitoxin itself. In this report, we first present the crystal structure of the M. tuberculosis HigA3 antitoxin (MtHigA3) and MtHigA3 bound to its operator DNA complex. We also investigated the interaction between MtHigA3 and DNA using NMR spectroscopy. The MtHigA3 antitoxin structure is a homodimer that contains a structurally well conserved DNA-binding domain at the N-terminus and a dimerization domain at the C-terminus. Upon comparing the HigA homologue structures, a distinct difference was found in the C-terminal region that possesses the β-lid, and diverse orientations of two helix-turn-helix (HTH) motifs from HigA homologue dimers were observed. The structure of MtHigA3 bound to DNA reveals that the promoter DNA is bound to two HTH motifs of the MtHigA3 dimer presenting 46.5° bending, and the distance between the two HTH motifs of each MtHigA3 monomer was increased in MtHigA3 bound to DNA. The β-lid, which is found only in the tertiary structure of MtHigA3 among the HigA homologues, causes the formation of a tight dimerization network and leads to a unique arrangement for dimer formation that is related to the curvature of the bound DNA. This work could contribute to the understanding of the HigBA system of M. tuberculosis at the atomic level and may contribute to the development of new antibiotics for TB treatment.

H-NS Silences Gene Expression of LeuO, the Master Regulator of the Cyclic(Phe-Pro)-dependent Signal Pathway, in Vibrio vulnificus

The histone-like nucleoid structuring protein (H-NS) is an abundant global regulator of environmentally controlled gene expression. Herein, we demonstrate that H-NS represses the expression of LeuO, the master regulator of the cyclic(Phe-Pro)-dependent signaling pathway, by directly binding to the upstream region of the gene. H-NS binds to a long stretched region (more than 160-bp long), which overlaps with binding sites for ToxR and LeuO. A high quantity of H-NS outcompetes ToxR for binding to the cis-acting element of leuO. However, our footprinting analyses suggests that the binding of H-NS is relatively weaker than LeuO or ToxR at the same molarity. Considering that the DNA nucleotide sequences of the upstream regions of leuO genes are highly conserved among various Vibrio, such patterns as those found in V. vulnificus would be a common feature in the regulation of leuO gene expression in Vibrionaceae. Taken together, these results suggest that, in species belonging to Vibrionaceae, H-NS regulates the expression of leuO as a basal stopper when cFP-ToxR mediated signaling is absent.

Nearest-neighbor parameters for predicting DNA duplex stability in diverse molecular crowding conditions

The intracellular environment is crowded and heterogeneous. Although the thermodynamic stability of nucleic acid duplexes is predictable in dilute solutions, methods of predicting such stability under specific intracellular conditions are not yet available. We recently showed that the nearest-neighbor model for self-complementary DNA is valid under molecular crowding condition of 40% polyethylene glycol with an average molecular weight of 200 (PEG 200) in 100 mM NaCl. Here, we determined nearest-neighbor parameters for DNA duplex formation under the same crowding condition to predict the thermodynamics of DNA duplexes in the intracellular environment. Preferential hydration of the nucleotides was found to be the key factor for nearest-neighbor parameters in the crowding condition. The determined parameters were shown to predict the thermodynamic parameters (∆H°, ∆S°, and ∆G°₃₇) and melting temperatures (T _m) of the DNA duplexes in the crowding condition with significant accuracy. Moreover, we proposed a general method for predicting the stability of short DNA duplexes in different cosolutes based on the relationship between duplex stability and the water activity of the cosolute solution. The method described herein would be valuable for investigating biological processes that occur under specific intracellular crowded conditions and for the application of DNA-based biotechnologies in crowded environments.

Multi-Factor Regulation of the Master Modulator LeuO for the Cyclic-(Phe-Pro) Signaling Pathway in Vibrio vulnificus

LeuO plays the role of a master regulator in the cyclic-L-phenylalanine-L-proline (cFP)-dependent signaling pathway in Vibrio vulnificus. cFP, as shown through isothermal titration calorimetry analysis, binds specifically to the periplasmic domain of ToxR. Binding of cFP triggers a change in the cytoplasmic domain of ToxR, which then activates transcription of leuO encoding a LysR-type regulator. LeuO binds to the region upstream of its own coding sequence, inhibiting its own transcription and maintaining a controlled level of expression. A five-bp deletion in this region abolished expression of LeuO, but a ten-bp deletion did not, suggesting that a DNA bending mechanism is involved in the regulation. Furthermore, binding of RNA polymerase was significantly lower both in the deletion of the ToxR binding site and in the five-bp deletion, but not in the ten-bp deletion, as shown in pull-down assays using an antibody against RNA polymerase subunit α. In summary, multiple factors are involved in control of the expression of LeuO, a master regulator that orchestrates downstream regulators to modulate factors required for survival and pathogenicity of the pathogen.

Specificity of protein-DNA interactions in hypersaline environment: structural studies on complexes of Halobacterium salinarum oxidative stress-dependent protein hsRosR

Interactions between proteins and DNA are crucial for all biological systems. Many studies have shown the dependence of protein–DNA interactions on the surrounding salt concentration. How these interactions are maintained in the hypersaline environments that halophiles inhabit remains puzzling. Towards solving this enigma, we identified the DNA motif recognized by the Halobactrium salinarum ROS-dependent transcription factor (hsRosR), determined the structure of several hsRosR–DNA complexes and investigated the DNA-binding process under extreme high-salt conditions. The picture that emerges from this work contributes to our understanding of the principles underlying the interplay between electrostatic interactions and salt-mediated protein–DNA interactions in an ionic environment characterized by molar salt concentrations.

Novel structural features drive DNA binding properties of Cmr, a CRP family protein in TB complex mycobacteria

Mycobacterium tuberculosis (Mtb) encodes two CRP/FNR family transcription factors (TF) that contribute to virulence, Cmr (Rv1675c) and CRP_Mt (Rv3676). Prior studies identified distinct chromosomal binding profiles for each TF despite their recognizing overlapping DNA motifs. The present study shows that Cmr binding specificity is determined by discriminator nucleotides at motif positions 4 and 13. X-ray crystallography and targeted mutational analyses identified an arginine-rich loop that expands Cmr’s DNA interactions beyond the classical helix-turn-helix contacts common to all CRP/FNR family members and facilitates binding to imperfect DNA sequences. Cmr binding to DNA results in a pronounced asymmetric bending of the DNA and its high level of cooperativity is consistent with DNA-facilitated dimerization. A unique N-terminal extension inserts between the DNA binding and dimerization domains, partially occluding the site where the canonical cAMP binding pocket is found. However, an unstructured region of this N-terminus may help modulate Cmr activity in response to cellular signals. Cmr’s multiple levels of DNA interaction likely enhance its ability to integrate diverse gene regulatory signals, while its novel structural features establish Cmr as an atypical CRP/FNR family member.

A duplex connection can further illuminate G-quadruplex/crystal violet complex

Duplex-connected G-quadruplex (dsG4) has a more evident promotion of the fluorescent emission of crystal violet than pure G-quadruplex (G4), which is firstly found and investigated systematically using fluorescence spectra, circular dichroism, and density functional theory. This new phenomenon also shows potential in application for target nucleic acid detection.

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

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

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Genes are sections of DNA that carry the instructions needed to build other molecules including all the proteins that the cell needs to fulfill its role. The information in the DNA is stored as a code consisting of four chemical bases, often referred to simply as “A”, “C”, “G” and “T”. The order or sequence of these bases determines the role of a protein. Many organisms – including humans – are built of many different types of cells that perform unique roles. Almost all cells carry the same genetic information, but proteins called transcription factors can regulate the activity of genes so that only a relevant subset of genes is switched on at a particular time.

Transcription factors glide along DNA and bind to short DNA sequences by attaching to the DNA bases directly or through bridges made up of water molecules. Two physical concepts known as enthalpy and entropy determine the strength of the connection. Enthalpy relates to how strong the chemical bonds that form between the transcription factors and the DNA bases are, compared to a situation where the transcription factor and DNA do not form a complex and bind to water molecules around them. Entropy measures the disorder of the system – the more disordered the solvent and protein-DNA complex are compared to solvent-containing free DNA and protein, the stronger the binding. A water molecule that bridges a DNA base with an amino-acid of a protein contributes to enthalpy, but results in loss of entropy, because the system becomes more ordered since the water molecule can no longer move freely.

Most transcription factors can only bind to DNA sequences that are very similar to each other, but some transcription factors can recognize several different kinds of sequences, and until now it was not clear how they could do this. Morgunova et al. studied four different human transcription factors that can each bind to two distinct DNA sequences. The results showed that the transcription factors bound to both DNA sequences with similar strength, but via different mechanisms. For one DNA sequence, an enthalpy-based mechanism essentially ‘froze’ the transcription factor to the DNA through rigid water bridges. The other DNA sequence was bound equally strongly but through moving water molecules, because this increased the entropy of the system. It is possible that these mechanisms could also apply to many other molecules that interact with each other through water-molecule bridges.

A better knowledge of the chemical bonds between transcription factors and DNA bases may in future help efforts to develop new treatments that depend on molecules being able to bind to other molecules. In addition, these findings may one day help scientists to predict how strongly two molecules will interact simply by knowing the structures of the molecules involved.

A unified approach for quantifying and interpreting DNA shape readout by transcription factors

Transcription factors (TFs) interpret DNA sequence by probing the chemical and structural properties of the nucleotide polymer. DNA shape is thought to enable a parsimonious representation of dependencies between nucleotide positions. Here, we propose a unified mathematical representation of the DNA sequence dependence of shape and TF binding, respectively, which simplifies and enhances analysis of shape readout. First, we demonstrate that linear models based on mononucleotide features alone account for 60-70% of the variance in minor groove width, roll, helix twist, and propeller twist. This explains why simple scoring matrices that ignore all dependencies between nucleotide positions can partially account for DNA shape readout by a TF Adding dinucleotide features as sequence-to-shape predictors to our model, we can almost perfectly explain the shape parameters. Building on this observation, we developed a post hoc analysis method that can be used to analyze any mechanism-agnostic protein-DNA binding model in terms of shape readout. Our insights provide an alternative strategy for using DNA shape information to enhance our understanding of how cis-regulatory codes are interpreted by the cellular machinery.

Probing the role of intercalating protein sidechains for kink formation in DNA

Protein binding can induce DNA kinks, which are for example important to enhance the specificity of the interaction and to facilitate the assembly of multi protein complexes. The respective proteins frequently exhibit amino acid sidechains that intercalate between the DNA base steps at the site of the kink. However, on a molecular level there is only little information available about the role of individual sidechains for kink formation. To unravel structural principles of protein-induced DNA kinking we have performed molecular dynamics (MD) simulations of five complexes that varied in their architecture, function, and identity of intercalated residues. Simulations were performed for the DNA complexes of wildtype proteins (Sac7d, Sox-4, CcpA, TFAM, TBP) and for mutants, in which the intercalating residues were individually or combined replaced by alanine. The work revealed that for systems with multiple intercalated residues, not all of them are necessarily required for kink formation. In some complexes (Sox-4, TBP), one of the residues proved to be essential for kink formation, whereas the second residue has only a very small effect on the magnitude of the kink. In other systems (e.g. Sac7d) each of the intercalated residues proved to be individually capable of conferring a strong kink suggesting a partially redundant role of the intercalating residues. Mutation of the key residues responsible for kinking either resulted in stable complexes with reduced kink angles or caused conformational instability as evidenced by a shift of the kink to an adjacent base step. Thus, MD simulations can help to identify the role of individual inserted residues for kinking, which is not readily apparent from an inspection of the static structures. This information might be helpful for understanding protein-DNA interactions in more detail and for designing proteins with altered DNA binding properties in the future.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Large-scale chromosome folding versus genomic DNA sequences: A discrete double Fourier transform technique

Using state-of-the-art techniques combining imaging methods and high-throughput genomic mapping tools leaded to the significant progress in detailing chromosome architecture of various organisms. However, a gap still remains between the rapidly growing structural data on the chromosome folding and the large-scale genome organization. Could a part of information on the chromosome folding be obtained directly from underlying genomic DNA sequences abundantly stored in the databanks? To answer this question, we developed an original discrete double Fourier transform (DDFT). DDFT serves for the detection of large-scale genome regularities associated with domains/units at the different levels of hierarchical chromosome folding. The method is versatile and can be applied to both genomic DNA sequences and corresponding physico-chemical parameters such as base-pairing free energy. The latter characteristic is closely related to the replication and transcription and can also be used for the assessment of temperature or supercoiling effects on the chromosome folding. We tested the method on the genome of E. coli K-12 and found good correspondence with the annotated domains/units established experimentally. As a brief illustration of further abilities of DDFT, the study of large-scale genome organization for bacteriophage PHIX174 and bacterium Caulobacter crescentus was also added. The combined experimental, modeling, and bioinformatic DDFT analysis should yield more complete knowledge on the chromosome architecture and genome organization.

Polarizable Force Field for DNA Based on the Classical Drude Oscillator: II. Microsecond Molecular Dynamics Simulations of Duplex DNA

The structure and dynamics of DNA are governed by a sensitive balance between base stacking and pairing, hydration, and interactions with ions. Force-field models that include explicit representations of electronic polarization are capable of more accurately modeling the subtle details of these interactions versus commonly used additive force fields. In this work, we validate our recently refined polarizable force field for DNA based on the classical Drude oscillator model, in which electronic degrees of freedom are represented as negatively charged particles attached to their parent atoms via harmonic springs. The previous version of the force field, called Drude-2013, produced stable A- and B-DNA trajectories on the order of hundreds of nanoseconds, but deficiencies were identified that included weak base stacking ultimately leading to distortion of B-DNA duplexes and unstable Z-DNA. As a result of extensive refinement of base nonbonded terms and bonded parameters in the deoxyribofuranose sugar and phosphodiester backbone, we demonstrate that the new version of the Drude DNA force field is capable of simulating A- and B-forms of DNA on the microsecond time scale and the resulting conformational ensembles agree well with a broad set of experimental properties, including solution X-ray scattering profiles. In addition, simulations of Z-form duplex DNA in its crystal environment are stable on the order of 100 ns. The revised force field is to be called Drude-2017.

Differential regulation of ssb genes in the nitrogen-fixing cyanobacterium, Anabaena sp. strain PCC71201

Anabaena sp. PCC7120 possesses three genes coding for single-stranded DNA-binding (SSB) protein, of which ssb1 was a single gene, and ssb2 and ssb3 are the first genes of their corresponding operons. Regulation of the truncated ssb genes, ssb1 (alr0088) and ssb2 (alr7559), was unaffected by N-status of growth. They were negatively regulated by the SOS-response regulatory protein LexA, as indicated by the (i) binding of Anabaena LexA to the LexA box of regulatory regions of ssb1 and ssb2, and (ii) decreased expression of the downstream gfp reporter gene in Escherichia coli upon co-expression of LexA. However, the full-length ssb gene, ssb3 (all4779), was regulated by the availability of Fe²⁺ and combined nitrogen, as indicated by (i) increase in the levels of SSB3 protein on Fe²⁺ -depletion and decrease under Fe²⁺ -excess conditions, and (ii) 1.5- to 1.6-fold decrease in activity under nitrogen-fixing conditions compared to nitrogen-supplemented conditions. The requirement of Fe²⁺ as a co-factor for repression by FurA and the increase in levels of FurA under nitrogen-deficient conditions in Anabaena (Lopez-Gomollon et al. 2007) indicated a possible regulation of ssb3 by FurA. This was substantiated by (i) the binding of FurA to the regulatory region of ssb3, (ii) repression of the expression of the downstream gfp reporter gene in E. coli upon co-expression of FurA, and (iii) negative regulation of ssb3 promoter activity by the upstream AT-rich region in Anabaena. This is the first report on possible role of FurA, an important protein for iron homeostasis, in DNA repair of cyanobacteria.

Study of DNA binding and bending by Bacillus subtilis GabR, a PLP-dependent transcription factor

BACKGROUND

GabR is a transcriptional regulator belonging to the MocR/GabR family, characterized by a N-terminal wHTH DNA-binding domain and a C-terminal effector binding and/or oligomerization domain, structurally homologous to aminotransferases (ATs). In the presence of γ-aminobutyrate (GABA) and pyridoxal 5'-phosphate (PLP), GabR activates the transcription of gabT and gabD genes involved in GABA metabolism.

METHODS

Here we report a biochemical and atomic force microscopy characterization of Bacillus subtilis GabR in complex with DNA. Complexes were assembled in vitro to study their stoichiometry, stability and conformation.

RESULTS

The fractional occupancy of the GabR cognate site suggests that GabR binds as a dimer with Kd of 10nM. Upon binding GabR bends the DNA by 80° as measured by anomalous electrophoretic mobility. With GABA we observed a decrease in affinity and conformational rearrangements compatible with a less compact nucleo-protein complex but no changes of the DNA bending angle. By employing promoter and GabR mutants we found that basic residues of the positively charged groove on the surface of the AT domain affect DNA affinity.

CONCLUSIONS

The present data extend current understanding of the GabR-DNA interaction and the effect of GABA and PLP. A model for the GabR-DNA complex, corroborated by a docking simulation, is proposed.

GENERAL SIGNIFICANCE

Characterization of the GabR DNA binding mode highlights the key role of DNA bending and interactions with bases outside the canonical direct repeats, and might be of general relevance for the action mechanism of MocR transcription factors.

Predictive Models of Recombination Rate Variation across the Drosophila melanogaster Genome

In all eukaryotic species examined, meiotic recombination, and crossovers in particular, occur non‐randomly along chromosomes. The cause for this non-random distribution remains poorly understood but some specific DNA sequence motifs have been shown to be enriched near crossover hotspots in a number of species. We present analyses using machine learning algorithms to investigate whether DNA motif distribution across the genome can be used to predict crossover variation in Drosophila melanogaster, a species without hotspots. Our study exposes a combinatorial non-linear influence of motif presence able to account for a significant fraction of the genome-wide variation in crossover rates at all genomic scales investigated, from 20% at 5-kb to almost 70% at 2,500-kb scale. The models are particularly predictive for regions with the highest and lowest crossover rates and remain highly informative after removing sub-telomeric and -centromeric regions known to have strongly reduced crossover rates. Transcriptional activity during early meiosis and differences in motif use between autosomes and the X chromosome add to the predictive power of the models. Moreover, we show that population-specific differences in crossover rates can be partly explained by differences in motif presence. Our results suggest that crossover distribution in Drosophila is influenced by both meiosis-specific chromatin dynamics and very local constitutive open chromatin associated with DNA motifs that prevent nucleosome stabilization. These findings provide new information on the genetic factors influencing variation in recombination rates and a baseline to study epigenetic mechanisms responsible for plastic recombination as response to different biotic and abiotic conditions and stresses.

Structures of an all-α protein running along the DNA major groove

Despite over 3300 protein–DNA complex structures have been reported in the past decades, there remain some unknown recognition patterns between protein and target DNA. The silkgland-specific transcription factor FMBP-1 from the silkworm Bombyx mori contains a unique DNA-binding domain of four tandem STPRs, namely the score and three amino acid peptide repeats. Here we report three structures of this STPR domain (termed BmSTPR) in complex with DNA of various lengths. In the presence of target DNA, BmSTPR adopts a zig-zag structure of three or four tandem α-helices that run along the major groove of DNA. Structural analyses combined with binding assays indicate BmSTPR prefers the AT-rich sequences, with each α-helix covering a DNA sequence of 4 bp. The successive AT-rich DNAs adopt a wider major groove, which is in complementary in shape and size to the tandem α-helices of BmSTPR. Substitutions of DNA sequences and affinity comparison further prove that BmSTPR recognizes the major groove mainly via shape readout. Multiple-sequence alignment suggests this unique DNA-binding pattern should be highly conserved for the STPR domain containing proteins which are widespread in animals. Together, our findings provide structural insights into the specific interactions between a novel DNA-binding protein and a unique deformed B-DNA.

Flanking A·T basepairs destabilize the B(∗) conformation of DNA A-tracts

Capillary electrophoresis has been used to characterize the interaction of monovalent cations with 26-basepair DNA oligomers containing A-tracts embedded in flanking sequences with different basepair compositions. A 26-basepair random-sequence oligomer was used as the reference; lithium and tetrabutylammonium (TBA(+)) ions were used as the probe ions. The free solution mobilities of the A-tract and random-sequence oligomers were identical in solutions containing <∼ 100 mM cation. At higher cation concentrations, the A-tract oligomers migrated faster than the reference oligomer in TBA(+) and slower than the reference in Li(+). Hence, cations of different sizes can interact very differently with DNA A-tracts. The increased mobilities observed in TBA(+) suggest that the large hydrophobic TBA(+) ions are preferentially excluded from the vicinity of the A-tract minor groove, increasing the effective net charge of the A-tract oligomers and increasing the mobility. By contrast, Li(+) ions decrease the mobility of A-tract oligomers because of the preferential localization of Li(+) ions in the narrow A-tract minor groove. Embedding the A-tracts in AT-rich flanking sequences markedly alters preferential interactions of monovalent cations with the B(∗) conformation. Hence, A-tracts embedded in genomic DNA may or may not interact preferentially with monovalent cations, depending on the relative number of A · T basepairs in the flanking sequences.

Binding of transcription factor GabR to DNA requires recognition of DNA shape at a location distinct from its cognate binding site

Mechanisms for transcription factor recognition of specific DNA base sequences are well characterized and recent studies demonstrate that the shape of these cognate binding sites is also important. Here, we uncover a new mechanism where the transcription factor GabR simultaneously recognizes two cognate binding sites and the shape of a 29 bp DNA sequence that bridges these sites. Small-angle X-ray scattering and multi-angle laser light scattering are consistent with a model where the DNA undergoes a conformational change to bend around GabR during binding. In silico predictions suggest that the bridging DNA sequence is likely to be bendable in one direction and kinetic analysis of mutant DNA sequences with biolayer interferometry, allowed the independent quantification of the relative contribution of DNA base and shape recognition in the GabR-DNA interaction. These indicate that the two cognate binding sites as well as the bendability of the DNA sequence in between these sites are required to form a stable complex. The mechanism of GabR-DNA interaction provides an example where the correct shape of DNA, at a clearly distinct location from the cognate binding site, is required for transcription factor binding and has implications for bioinformatics searches for novel binding sites.

OxyR-dependent formation of DNA methylation patterns in OpvABOFF and OpvABON cell lineages of Salmonella enterica

Phase variation of the Salmonella enterica opvAB operon generates a bacterial lineage with standard lipopolysaccharide structure (OpvAB(OFF)) and a lineage with shorter O-antigen chains (OpvAB(ON)). Regulation of OpvAB lineage formation is transcriptional, and is controlled by the LysR-type factor OxyR and by DNA adenine methylation. The opvAB regulatory region contains four sites for OxyR binding (OBSA-D), and four methylatable GATC motifs (GATC1-4). OpvAB(OFF) and OpvAB(ON) cell lineages display opposite DNA methylation patterns in the opvAB regulatory region: (i) in the OpvAB(OFF) state, GATC1 and GATC3 are non-methylated, whereas GATC2 and GATC4 are methylated; (ii) in the OpvAB(ON) state, GATC2 and GATC4 are non-methylated, whereas GATC1 and GATC3 are methylated. We provide evidence that such DNA methylation patterns are generated by OxyR binding. The higher stability of the OpvAB(OFF) lineage may be caused by binding of OxyR to sites that are identical to the consensus (OBSA and OBSc), while the sites bound by OxyR in OpvAB(ON) cells (OBSB and OBSD) are not. In support of this view, amelioration of either OBSB or OBSD locks the system in the ON state. We also show that the GATC-binding protein SeqA and the nucleoid protein HU are ancillary factors in opvAB control.

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.

Single-Walled Carbon Nanotubes Modulate the B- to A-DNA Transition

We study the conformational equilibrium between B-to-A forms of ds-DNA adsorbed onto a single-walled carbon nanotube (SWNT) using free energy profile calculations based on all-atom molecular dynamics simulations. The potential of mean force (PMF) of the B-to-A transition of ds-DNA in the presence of an uncharged (10,0) carbon nanotube for two dodecamers with poly-AT or poly-GC sequences is calculated as a function of a root-mean-square-distance (ΔRMSD) difference metric for the B-to-A transition. The calculations reveal that in the presence of a SWNT DNA favors B-form DNA significantly in both poly-GC and poly-AT sequences. Furthermore, the poly-AT DNA:SWNT complex shows a higher energy penalty for adopting an A-like conformation than poly-GC DNA:SWNT by several kcal/mol. The presence of a SWNT on either poly-AT or poly-GC DNA affects the PMF of the transition such that the B form is favored by as much as 10 kcal/mol. In agreement with published data, we find a potential energy minimum between A and B-form DNA at ΔRMSD ≈ -1.5 Å and that the presence of the SWNT moves this minimum by as much as ΔRMSD = 3 Å.

Mechanical Model of DNA Allostery

The importance of allosteric effects in DNA is becoming increasingly appreciated, but the underlying mechanisms remain poorly understood. In this work, we propose a general modeling framework to study DNA allostery. We describe DNA in a coarse-grained manner by intra-base pair and base pair step coordinates, complemented by groove widths. Quadratic deformation energy is assumed, yielding linear relations between the constraints and their effect. Model parameters are inferred from standard unrestrained, explicit-solvent molecular dynamics simulations of naked DNA. We applied the approach to study minor groove binding of diamidines and pyrrole-imidazole polyamides. The predicted DNA bending is in quantitative agreement with experiment and suggests that diamidine binding to the alternating TA sequence brings the DNA closer to the A-tract conformation, with potentially important functional consequences. The approach can be readily applied to other allosteric effects in DNA and generalized to model allostery in various molecular systems.

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

An A-tract at the AtzR binding site assists DNA binding, inducer-dependent repositioning and transcriptional activation of the PatzDEF promoter

The LysR-type regulator AtzR activates the Pseudomonas sp. ADP atzDEF operon in response to nitrogen limitation and cyanuric acid. Activation involves repositioning of the AtzR tetramer on the PatzDEF promoter and relaxation of an AtzR-induced DNA bend. Here we examine the in vivo and in vitro contribution of an A5 -tract present at the PatzDEF promoter region to AtzR binding and transcriptional activation. Substitution of the A-tract for the sequence ACTCA prevented PatzDEF activation and high-affinity AtzR binding, impaired AtzR contacts with the activator binding site and shifted the position of the AtzR-induced DNA bend. Analysis of a collection of mutants bearing different alterations in the A-tract sequence showed that the extent of AtzR-dependent activation does not correlate with the magnitude or orientation of the spontaneous DNA bend generated at this site. Our results support the notion that indirect readout of the A-tract-associated narrow minor groove is essential for the AtzR-DNA complex to achieve a conformation competent for activation of the PatzDEF promoter. Conservation of this motif in several binding sites of LysR-type regulators suggests that this mechanism may be shared by other proteins in this family.

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.

Structural ensemble and microscopic elasticity of freely diffusing DNA by direct measurement of fluctuations

Precisely measuring the ensemble of conformers that a macromolecule populates in solution is highly challenging. Thus, it has been difficult to confirm or falsify the predictions of nanometer-scale dynamical modeling. Here, we apply an X-ray interferometry technique to probe the solution structure and fluctuations of B-form DNA on a length scale comparable to a protein-binding site. We determine an extensive set of intrahelix distance distributions between pairs of probes placed at distinct points on the surface of the DNA duplex. The distributions of measured distances reveal the nature and extent of the thermally driven mechanical deformations of the helix. We describe these deformations in terms of elastic constants, as is common for DNA and other polymers. The average solution structure and microscopic elasticity measured by X-ray interferometry are in striking agreement with values derived from DNA-protein crystal structures and measured by force spectroscopy, with one exception. The observed microscopic torsional rigidity of DNA is much lower than is measured by single-molecule twisting experiments, suggesting that torsional rigidity increases when DNA is stretched. Looking forward, molecular-level interferometry can provide a general tool for characterizing solution-phase structural ensembles.