Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

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.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

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.

Examining the Effects of Netropsin on the Curvature of DNA A-Tracts Using Electrophoresis

A-tracts are sequences of repeated adenine bases that, under the proper conditions, are capable of mediating DNA curvature. A-tracts occur naturally in the regulatory regions of many organisms, yet their biological functions are not fully understood. Orienting multiple A-tracts together constructively or destructively in a phase has the potential to create different shapes in the DNA helix axis. One means of detecting these molecular shape differences is from altered DNA mobilities measured using electrophoresis. The small molecule netropsin binds the minor groove of DNA, particularly at AT-rich sequences including A-tracts. Here, we systematically test the hypothesis that netropsin binding eliminates the curvature of A-tracts by measuring the electrophoretic mobilities of seven 98-base pair DNA samples containing different numbers and arrangements of centrally located A-tracts under varying conditions with netropsin. We find that netropsin binding eliminates the mobility difference between the DNA fragments with different A-tract arrangements in a concentration-dependent manner. This work provides evidence for the straightening of A-tracts upon netropsin binding and illustrates an artificial approach to re-sculpt DNA shape.

Epigenetic competition reveals density-dependent regulation and target site plasticity of phosphorothioate epigenetics in bacteria

Phosphorothioate (PT) DNA modifications-in which a nonbonding phosphate oxygen is replaced with sulfur-represent a widespread, horizontally transferred epigenetic system in prokaryotes and have a highly unusual property of occupying only a small fraction of available consensus sequences in a genome. Using Salmonella enterica as a model, we asked a question of fundamental importance: How do the PT-modifying DndA-E proteins select their G_PSAAC/G_PSTTC targets? Here, we applied innovative analytical, sequencing, and computational tools to discover a novel behavior for DNA-binding proteins: The Dnd proteins are "parked" at the G^6mATC Dam methyltransferase consensus sequence instead of the expected GAAC/GTTC motif, with removal of the ^6mA permitting extensive PT modification of GATC sites. This shift in modification sites further revealed a surprising constancy in the density of PT modifications across the genome. Computational analysis showed that GAAC, GTTC, and GATC share common features of DNA shape, which suggests that PT epigenetics are regulated in a density-dependent manner partly by DNA shape-driven target selection in the genome.

Modulation of the helical properties of DNA: next-to-nearest neighbour effects and beyond

We used extensive molecular dynamics simulations to study the structural and dynamic properties of the central d(TpA) step in the highly polymorphic d(CpTpApG) tetranucleotide. Contrary to the assumption of the dinucleotide-model and its nearest neighbours (tetranucleotide-model), the properties of the central d(TpA) step change quite significantly dependent on the next-to-nearest (hexanucleotide) sequence context and in a few cases are modulated by even remote neighbours (beyond next-to-nearest from the central TpA). Our results highlight the existence of previously undescribed dynamical mechanisms for the transmission of structural information into the DNA and demonstrate the existence of certain sequences with special physical properties that can impact on the global DNA structure and dynamics.

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.

Ml proteins from Mesorhizobium loti and MucR from Brucella abortus: an AT-rich core DNA-target site and oligomerization ability

Mesorhizobium loti contains ten genes coding for proteins sharing high amino acid sequence identity with members of the Ros/MucR transcription factor family. Five of these Ros/MucR family members from Mesorhizobium loti (Ml proteins) have been recently structurally and functionally characterized demonstrating that Ml proteins are DNA-binding proteins. However, the DNA-binding studies were performed using the Ros DNA-binding site with the Ml proteins. Currently, there is no evidence as to when the Ml proteins are expressed during the Mesorhizobium lo ti life cycle as well as no information concerning their natural DNA-binding site. In this study, we examine the ml genes expression profile in Mesorhizobium loti and show that ml1, ml2, ml3 and ml5 are expressed during planktonic growth and in biofilms. DNA-binding experiments show that the Ml proteins studied bind a conserved AT-rich site in the promoter region of the exoY gene from Mesorhizobium loti and that the proteins make important contacts with the minor groove of DNA. Moreover, we demonstrate that the Ml proteins studied form higher-order oligomers through their N-terminal region and that the same AT-rich site is recognized by MucR from Brucella abortus using a similar mechanism involving contacts with the minor groove of DNA and oligomerization.

Multiple Electronic and Structural Factors Control Cyclobutane Pyrimidine Dimer and 6-4 Thymine-Thymine Photodimerization in a DNA Duplex

The T-T photodimerization paths leading to the formation of cyclobutane pyrimidine dimer (CPD) and 6-4 pyrimidine pyrimidone (64-PP), the two main DNA photolesions, have been resolved for a T-T step in a DNA duplex by two complementary state-of-the-art quantum mechanical approaches: QM(CASPT2//CASSCF)/MM and TD-DFT/PCM. Based on the analysis of several different representative structures, we define a new-ensemble of cooperating geometrical and electronic factors (besides the distance between the reacting bonds) ruling T-T photodimerization in DNA. CPD is formed by a barrierless path on an exciton state delocalized over the two bases. Large interbase stacking and shift values, together with a small pseudorotation phase angle for T at the 3'-end, favor this reaction. The oxetane intermediate, leading to a 64-PP adduct, is formed on a singlet T→T charge-transfer state and is favored by a large interbase angle and slide values. A small energy barrier (<0.3 eV) is associated to this path, likely contributing to the smaller quantum yield observed for this process. Eventually, a clear directionality is always shown by the electronic excitation characterizing the singlet photoactive state driving the photodimerization process: an exciton that is more localized on T³ and a 5'-T→3'-T charge transfer for CPD and oxetane formation, respectively, thus calling for specific electronic constraints.

Structural dynamics of wild type and mutated forms of human L1 endonuclease and insights into its sequence specific nucleic acid binding mechanism: A molecular dynamics study

Biomolecular recognition of proteins and nucleic acids is mainly mediated by their structural features and the molecular dynamics simulations approach has been used to explore this recognition processes at the atomic level. L1-Endonuclease, an enzyme involved in L1 retrotransposition, cleaves the TA junction DNA (5'-TTTT/AA-3') and expresses high specificity for target site recognition. The present study highlights the structural features of L1-endonuclease as well as DNA responsible for such specific recognition. Especially, the importance of βB6-B5 hairpin loop in DNA recognition has been elucidated by analyzing the dynamics of Thr192 mutated L1-endonuclease. In addition, simulations of the endonuclease complexed with DNA substrates (sequences having TA and CG junctions) revealed the specificity of L1 endonuclease towards TA junction. Molecular dynamics simulations revealed that the βB6-B5 hairpin loop protrudes well into the minor groove of DNA having TA junction and induces DNA bending such that the width of minor groove is increased. Such endonuclease induced bending of TA junction DNA sequence positions the scissile phosphodiester bond of DNA for cleavage. The innate property of minor groove widening in TA junction than in CG junction is utilized by the βB6-βB5 hairpin loop of endonuclease while recognizing the DNA sequences. The present study also highlights the role of Mg²⁺ cation in catalysis and attempts to explore the possible target site DNA cleavage mechanism.

3D-NuS: A Web Server for Automated Modeling and Visualization of Non-Canonical 3-Dimensional Nucleic Acid Structures

The inherent conformational flexibility of nucleic acids facilitates the formation of a range of conformations such as duplex, triplex, quadruplex, etc. that play crucial roles in biological processes. Elucidation of the influence of non-canonical base pair mismatches on DNA/RNA structures at different sequence contexts to understand the mismatch repair, misregulation of alternative splicing mechanisms and the sequence-dependent effect of RNA-DNA hybrid in relevance to antisense strategy demand their three-dimensional structural information. Furthermore, structural insights about nucleic acid triplexes, which are generally not tractable to structure determination by X-ray crystallography or NMR techniques, are essential to establish their biological function(s). A web server, namely 3D-NuS (http://iith.ac.in/3dnus/), has been developed to generate energy-minimized models of 80 different types of triplexes, 64 types of G-quadruplexes, left-handed Z-DNA/RNA duplexes, and RNA-DNA hybrid duplex along with inter- and intramolecular DNA or RNA duplexes comprising a variety of mismatches and their chimeric forms for any user-defined sequence and length. It also generates an ensemble of conformations corresponding to the modeled structure. These structures may serve as good starting models for docking proteins and small molecules with nucleic acids, NMR structure determination, cryo-electron microscope modeling, DNA/RNA nanotechnology applications and molecular dynamics simulation studies.

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.

μABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

All-atom crystal simulations of DNA and RNA duplexes

BACKGROUND

Molecular dynamics simulations can complement experimental measures of structure and dynamics of biomolecules. The quality of such simulations can be tested by comparisons to models refined against experimental crystallographic data.

METHODS

We report simulations of DNA and RNA duplexes in their crystalline environment. The calculations mimic the conditions for PDB entries 1D23 [d(CGATCGATCG)2] and 1RNA [(UUAUAUAUAUAUAA)2], and contain 8 unit cells, each with 4 copies of the Watson-Crick duplex; this yields in aggregate 64μs of duplex sampling for DNA and 16μs for RNA.

RESULTS

The duplex structures conform much more closely to the average structure seen in the crystal than do structures extracted from a solution simulation with the same force field. Sequence-dependent variations in helical parameters, and in groove widths, are largely maintained in the crystal structure, but are smoothed out in solution. However, the integrity of the crystal lattice is slowly degraded in both simulations, with the result that the interfaces between chains become heterogeneous. This problem is more severe for the DNA crystal, which has fewer inter-chain hydrogen bond contacts than does the RNA crystal.

CONCLUSIONS

Crystal simulations using current force fields reproduce many features of observed crystal structures, but suffer from a gradual degradation of the integrity of the crystal lattice.

GENERAL SIGNIFICANCE

The results offer insights into force-field simulations that test their ability to preserve weak interactions between chains, which will be of importance also in non-crystalline applications that involve binding and recognition. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

NMR studies of DNA support the role of pre-existing minor groove variations in nucleosome indirect readout

We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI ↔ BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.

Local dynamics coupled to hydration water determines DNA-sequence-dependent deformability

Molecular dynamics (MD) simulations and quasielastic neutron scattering (QENS) experiments were conducted on two hydrated DNA dodecamers with distinct deformability: 5'CGCG[under AATT]̲CGCG3' and 5'CGCG[under TTAA]̲CGCG3'. The former is known to be rigid and the latter to be flexible. The mean-square displacements of DNA dodecamers exhibit so-called dynamical transition around 200-240 K for both sequences. To investigate the DNA-sequence-dependent dynamics, the dynamics of DNA and hydration water above the transition temperature were examined using both MD simulations and QENS experiments. The fluctuation amplitude of the AATT central tetramer is smaller, and its relaxation time is longer, than that observed in TTAA, suggesting that the AT step is kinetically more stable than TA. The sequence-dependent local base pair step dynamics correlates with the kinetics of breaking the hydrogen bond between DNA and hydration water. The sequence-dependent DNA base pair step fluctuations appear above the dynamical transition temperature. Together with these results, we conclude that DNA deformability is related to the local dynamics of the base pair steps, themselves coupled to hydration water in the minor groove.

Making the bend: DNA tertiary structure and protein-DNA interactions

DNA structure functions as an overlapping code to the DNA sequence. Rapid progress in understanding the role of DNA structure in gene regulation, DNA damage recognition and genome stability has been made. The three dimensional structure of both proteins and DNA plays a crucial role for their specific interaction, and proteins can recognise the chemical signature of DNA sequence ("base readout") as well as the intrinsic DNA structure ("shape recognition"). These recognition mechanisms do not exist in isolation but, depending on the individual interaction partners, are combined to various extents. Driving force for the interaction between protein and DNA remain the unique thermodynamics of each individual DNA-protein pair. In this review we focus on the structures and conformations adopted by DNA, both influenced by and influencing the specific interaction with the corresponding protein binding partner, as well as their underlying thermodynamics.

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.

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.

Indirect readout of DNA sequence by p22 repressor: roles of DNA and protein functional groups in modulating DNA conformation

The repressor of bacteriophage P22 (P22R) discriminates between its various DNA binding sites by sensing the identity of non-contacted base pairs at the center of its binding site. The "indirect readout" of these non-contacted bases is apparently based on DNA's sequence-dependent conformational preferences. The structures of P22R-DNA complexes indicate that the non-contacted base pairs at the center of the binding site are in the B' state. This finding suggests that indirect readout and therefore binding site discrimination depend on P22R's ability to either sense and/or impose the B' state on the non-contacted bases of its binding sites. We show here that the affinity of binding sites for P22R depends on the tendency of the central bases to assume the B'-DNA state. Furthermore, we identify functional groups in the minor groove of the non-contacted bases as the essential modulators of indirect readout by P22R. In P22R-DNA complexes, the negatively charged E44 and E48 residues are provocatively positioned near the negatively charged DNA phosphates of the non-contacted nucleotides. The close proximity of the negatively charged groups on protein and DNA suggests that electrostatics may play a key role in the indirect readout process. Changing either of two negatively charged residues to uncharged residues eliminates the ability of P22R to impose structural changes on DNA and to recognize non-contacted base sequence. These findings suggest that these negatively charged amino acids function to force the P22R-bound DNA into the B' state and therefore play a key role in indirect readout by P22R.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

Role of sequence encoded κB DNA geometry in gene regulation by Dorsal

Many proteins of the Rel family can act as both transcriptional activators and repressors. However, mechanism that discerns the ‘activator/repressor’ functions of Rel-proteins such as Dorsal (Drosophila homologue of mammalian NFκB) is not understood. Using genomic, biophysical and biochemical approaches, we demonstrate that the underlying principle of this functional specificity lies in the ‘sequence-encoded structure’ of the κB-DNA. We show that Dorsal-binding motifs exist in distinct activator and repressor conformations. Molecular dynamics of DNA-Dorsal complexes revealed that repressor κB-motifs typically have A-tract and flexible conformation that facilitates interaction with co-repressors. Deformable structure of repressor motifs, is due to changes in the hydrogen bonding in A:T pair in the ‘A-tract’ core. The sixth nucleotide in the nonameric κB-motif, ‘A’ (A₆) in the repressor motifs and ‘T’ (T₆) in the activator motifs, is critical to confer this functional specificity as A₆ → T₆ mutation transformed flexible repressor conformation into a rigid activator conformation. These results highlight that ‘sequence encoded κB DNA-geometry’ regulates gene expression by exerting allosteric effect on binding of Rel proteins which in turn regulates interaction with co-regulators. Further, we identified and characterized putative repressor motifs in Dl-target genes, which can potentially aid in functional annotation of Dorsal gene regulatory network.

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes -2.3 kcal mol(-1) to -0.9 kcal mol(-1) (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of +1.3 kcal mol(-1). Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.

Pyrimidine-purine and pyrimidine heterodinucleosides synthesis containing a triazole linkage

This article describes a synthetic route to generate two purine-pyrimidine and pyrimidine heterodinucleosides. Both microwave activated regioselective alkylation using hydride and copper-catalyzed-azide-alkyne-cycloaddition (CuAAC) were used in order to perform the synthesis.

Structural insight into the sequence dependence of nucleosome positioning

Nucleosome positioning displays sequence dependency and contributes to genomic regulation in a site-specific manner. We solved the structures of nucleosome core particle composed of strong positioning TTTAA elements flanking the nucleosome center. The positioning strength of the super flexible TA dinucleotide is consistent with its observed central location within minor groove inward regions, where it can contribute maximally to energetically challenging minor groove bending, kinking and compression. The marked preference for TTTAA and positioning power of the site 1.5 double helix turns from the nucleosome center relates to a unique histone protein motif at this location, which enforces a sustained, extremely narrow minor groove via a hydrophobic "sugar clamp." Our analysis sheds light on the basis of nucleosome positioning and indicates that the histone octamer has evolved not to fully minimize sequence discrimination in DNA binding.

Base-base recognition of nonionic dinucleotide analogues in an apolar environment studied by low-temperature NMR spectroscopy

Two self-complementary dinucleotide analogues T(Si)A and A(Si)T with a nonionic diisopropylsilyl-modified backbone were synthesized, and their association in a nonaqueous aprotic environment was studied by NMR spectroscopy. Using a CDClF(2)/CDF(3) solvent mixture, measurements at temperatures as low as 113 K allowed the observation and structural characterization of individual complexes in the slow exchange regime. The A(Si)T analogue associates to exclusively form a dinucleotide antiparallel duplex with regular Watson-Crick base pairing, but both A and T nucleosides exhibit a predominant C3'-endo sugar pucker reminiscent of an A-type conformation. In contrast to A(Si)T, the T(Si)A dinucleotide is found to exhibit significant variability and flexibility. Thus, different secondary structures with weaker hydrogen bonds for all T(Si)A structures are observed at low temperatures. Although a B-like Watson-Crick antiparallel dinucleotide duplex with a preferred C2'-endo sugar pucker largely predominates at temperatures above 153 K, two additional species, namely a dinucleotide Hoogsteen duplex with a syn glycosidic torsion angle of the adenosine nucleoside and a presumably intramolecularly folded structure, are increasingly populated upon further cooling. By adding typical DNA intercalators like anthracene or benz[c]acridine derivatives to the A(Si)T dinucleotide duplex in the aprotic solvent environment, no binding of the polycyclic aromatic molecules can be detected even at lower temperatures. Obviously, van der Waals and stacking interactions are insufficient to compensate for the other unfavorable contributions to the overall free energy of binding, and only in the presence of additional hydrophobic effects in an aqueous environment does binding occur.

Origins of specificity in protein-DNA recognition

Specific interactions between proteins and DNA are fundamental to many biological processes. In this review, we provide a revised view of protein-DNA interactions that emphasizes the importance of the three-dimensional structures of both macromolecules. We divide protein-DNA interactions into two categories: those when the protein recognizes the unique chemical signatures of the DNA bases (base readout) and those when the protein recognizes a sequence-dependent DNA shape (shape readout). We further divide base readout into those interactions that occur in the major groove from those that occur in the minor groove. Analogously, the readout of the DNA shape is subdivided into global shape recognition (for example, when the DNA helix exhibits an overall bend) and local shape recognition (for example, when a base pair step is kinked or a region of the minor groove is narrow). Based on the >1500 structures of protein-DNA complexes now available in the Protein Data Bank, we argue that individual DNA-binding proteins combine multiple readout mechanisms to achieve DNA-binding specificity. Specificity that distinguishes between families frequently involves base readout in the major groove, whereas shape readout is often exploited for higher resolution specificity, to distinguish between members within the same DNA-binding protein family.

A measure of bending in nucleic acids structures applied to A-tract DNA

A method is proposed to measure global bending in DNA and RNA structures. It relies on a properly defined averaging of base-fixed coordinate frames, computes mean frames of suitably chosen groups of bases and uses these mean frames to evaluate bending. The method is applied to DNA A-tracts, known to induce considerable bend to the double helix. We performed atomistic molecular dynamics simulations of sequences containing the A(4)T(4) and T(4)A(4) tracts, in a single copy and in two copies phased with the helical repeat. Various temperature and salt conditions were investigated. Our simulations indicate bending by roughly 10 degrees per A(4)T(4) tract into the minor groove, and an essentially straight structure containing T(4)A(4), in agreement with electrophoretic mobility data. In contrast, we show that the published NMR structures of analogous sequences containing A(4)T(4) and T(4)A(4) tracts are significantly bent into the minor groove for both sequences, although bending is less pronounced for the T(4)A(4) containing sequence. The bending magnitudes obtained by frame averaging are confirmed by the analysis of superhelices composed of repeated tract monomers.

Structure-based analysis of DNA sequence patterns guiding nucleosome positioning in vitro

Recent studies of genome-wide nucleosomal organization suggest that the DNA sequence is one of the major determinants of nucleosome positioning. Although the search for underlying patterns encoded in nucleosomal DNA has been going on for about 30 years, our knowledge of these patterns still remains limited. Based on our evaluations of DNA deformation energy, we developed new scoring functions to predict nucleosome positioning. There are three principal differences between our approach and earlier studies: (i) we assume that the length of nucleosomal DNA varies from 146 to 147 bp; (ii) we consider the anisotropic flexibility of pyrimidine-purine (YR) dimeric steps in the context of their neighbors (e.g., YYRR versus RYRY); (iii) we postulate that alternating AT-rich and GC-rich motifs reflect sequence-dependent interactions between histone arginines and DNA in the minor groove. Using these functions, we analyzed 20 nucleosome positions mapped in vitro at single nucleotide resolution (including clones 601, 603, 605, the pGUB plasmid, chicken beta-globin and three 5S rDNA genes). We predicted 15 of the 20 positions with 1-bp precision, and two positions with 2-bp precision. The predicted position of the '601' nucleosome (i.e., the optimum of the computed score) deviates from the experimentally determined unique position by no more than 1 bp - an accuracy exceeding that of earlier predictions. Our analysis reveals a clear heterogeneity of the nucleosomal sequences which can be divided into two groups based on the positioning 'rules' they follow. The sequences of one group are enriched by highly deformable YR/YYRR motifs at the minor-groove bending sites SHL+/- 3.5 and +/- 5.5, which is similar to the alpha-satellite sequence used in most crystallized nucleosomes. Apparently, the positioning of these nucleosomes is determined by the interactions between histones H2A/H2B and the terminal parts of nucleosomal DNA. In the other group (that includes the '601' clone) the same YR/YYRR motifs occur predominantly at the sites SHL +/- 1.5. The interaction between the H3/H4 tetramer and the central part of the nucleosomal DNA is likely to be responsible for the positioning of nucleosomes of this group, and the DNA trajectory in these nucleosomes may differ in detail from the published structures. Thus, from the stereochemical perspective, the in vitro nucleosomes studied here follow either an X-ray-like pattern (with strong deformations in the terminal parts of nucleosomal DNA), or an alternative pattern (with the deformations occurring predominantly in the central part of the nucleosomal DNA). The results presented here may be useful for genome-wide classification of nucleosomes, linking together structural and thermodynamic characteristics of nucleosomes with the underlying DNA sequence patterns guiding their positions.

Direct measurements of the nucleosome-forming preferences of periodic DNA motifs challenge established models

Several periodic motifs have been implicated in facilitating the bending of DNA around the histone core of the nucleosome. For example, di-nucleotides AA/TT/TA and GC at ∼10-bp periods, but offset by 5 bp, are found with higher-than-expected occurrences in aligned nucleosomal DNAs in vitro and in vivo. Additionally, regularly oscillating period-10 trinucleotide motifs non-T, A/T, G and their complements have been implicated in the formation of regular nucleosome arrays. The effects of these periodic motifs on nucleosome formation have not been systematically tested directly by competitive reconstitution assays. We show that, in general, none of these period-10 motifs, except TA, in certain sequence contexts, facilitates nucleosome formation. The influence of periodic TAs on nucleosome formation is appreciable; with some of the 200-bp DNAs out-competing bulk nucleosomal DNA by more than 400-fold. Only the nucleotides immediately flanking TA influence its nucleosome-forming ability. Period-10 TA, when flanked by a pair of permissive nucleotides, facilitates DNA bending through compression of the minor groove. The free energy change for nucleosome formation decreases linearly with the number of consecutive TAs, up to eight. We suggest how these data can be reconciled with previous findings.

The shape of the DNA minor groove directs binding by the DNA-bending protein Fis

The bacterial nucleoid-associated protein Fis regulates diverse reactions by bending DNA and through DNA-dependent interactions with other control proteins and enzymes. In addition to dynamic nonspecific binding to DNA, Fis forms stable complexes with DNA segments that share little sequence conservation. Here we report the first crystal structures of Fis bound to high- and low-affinity 27-base-pair DNA sites. These 11 structures reveal that Fis selects targets primarily through indirect recognition mechanisms involving the shape of the minor groove and sequence-dependent induced fits over adjacent major groove interfaces. The DNA shows an overall curvature of approximately 65 degrees , and the unprecedented close spacing between helix-turn-helix motifs present in the apodimer is accommodated by severe compression of the central minor groove. In silico DNA structure models show that only the roll, twist, and slide parameters are sufficient to reproduce the changes in minor groove widths and recreate the curved Fis-bound DNA structure. Models based on naked DNA structures suggest that Fis initially selects DNA targets with intrinsically narrow minor grooves using the separation between helix-turn-helix motifs in the Fis dimer as a ruler. Then Fis further compresses the minor groove and bends the DNA to generate the bound structure.

Sequence dependencies of DNA deformability and hydration in the minor groove

DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

The promoter spacer influences transcription initiation via sigma70 region 1.1 of Escherichia coli RNA polymerase

Transcription initiation is a dynamic process in which RNA polymerase (RNAP) and promoter DNA act as partners, changing in response to one another, to produce a polymerase/promoter open complex (RPo) competent for transcription. In Escherichia coli RNAP, region 1.1, the N-terminal 100 residues of sigma(70), is thought to occupy the channel that will hold the DNA downstream of the transcription start site; thus, region 1.1 must move from this channel as RPo is formed. Previous work has also shown that region 1.1 can modulate RPo formation depending on the promoter. For some promoters region 1.1 stimulates the formation of open complexes; at the P(minor) promoter, region 1.1 inhibits this formation. We demonstrate here that the AT-rich P(minor) spacer sequence, rather than promoter recognition elements or downstream DNA, determines the effect of region 1.1 on promoter activity. Using a P(minor) derivative that contains good sigma(70)-dependent DNA elements, we find that the presence of a more GC-rich spacer or a spacer with the complement of the P(minor) sequence results in a promoter that is no longer inhibited by region 1.1. Furthermore, the presence of the P(minor) spacer, the GC-rich spacer, or the complement spacer results in different mobilities of promoter DNA during gel electrophoresis, suggesting that the spacer regions impart differing conformations or curvatures to the DNA. We speculate that the spacer can influence the trajectory or flexibility of DNA as it enters the RNAP channel and that region 1.1 acts as a "gatekeeper" to monitor channel entry.

CSI-Tree: a regression tree approach for modeling binding properties of DNA-binding molecules based on cognate site identification (CSI) data

The identification and characterization of binding sites of DNA-binding molecules, including transcription factors (TFs), is a critical problem at the interface of chemistry, biology and molecular medicine. The Cognate Site Identification (CSI) array is a high-throughput microarray platform for measuring comprehensive recognition profiles of DNA-binding molecules. This technique produces datasets that are useful not only for identifying binding sites of previously uncharacterized TFs but also for elucidating dependencies, both local and nonlocal, between the nucleotides at different positions of the recognition sites. We have developed a regression tree technique, CSI-Tree, for exploring the spectrum of binding sites of DNA-binding molecules. Our approach constructs regression trees utilizing the CSI data of unaligned sequences. The resulting model partitions the binding spectrum into homogeneous regions of position specific nucleotide effects. Each homogeneous partition is then summarized by a position weight matrix (PWM). Hence, the final outcome is a binding intensity rank-ordered collection of PWMs each of which spans a different region in the binding spectrum. Nodes of the regression tree depict the critical position/nucleotide combinations. We analyze the CSI data of the eukaryotic TF Nkx-2.5 and two engineered small molecule DNA ligands and obtain unique insights into their binding properties. The CSI tree for Nkx-2.5 reveals an interaction between two positions of the binding profile and elucidates how different nucleotide combinations at these two positions lead to different binding affinities. The CSI trees for the engineered DNA ligands exhibit a common preference for the dinucleotide AA in the first two positions, which is consistent with preference for a narrow and relatively flat minor groove. We carry out a reanalysis of these data with a mixture of PWMs approach. This approach is an advancement over the simple PWM model and accommodates position dependencies based on only sequence data. Our analysis indicates that the dependencies revealed by the CSI-Tree are challenging to discover without the actual binding intensities. Moreover, such a mixture model is highly sensitive to the number and length of the sequences analyzed. In contrast, CSI-Tree provides interpretable and concise summaries of the complete recognition profiles of DNA-binding molecules by utilizing binding affinities.

Magnitude and direction of DNA bending induced by screw-axis orientation: influence of sequence, mismatches and abasic sites

DNA-bending flexibility is central for its many biological functions. A new bending restraining method for use in molecular mechanics calculations and molecular dynamics simulations was developed. It is based on an average screw rotation axis definition for DNA segments and allows inducing continuous and smooth bending deformations of a DNA oligonucleotide. In addition to controlling the magnitude of induced bending it is also possible to control the bending direction so that the calculation of a complete (2-dimensional) directional DNA-bending map is now possible. The method was applied to several DNA oligonucleotides including A(adenine)-tract containing sequences known to form stable bent structures and to DNA containing mismatches or an abasic site. In case of G:A and C:C mismatches a greater variety of conformations bent in various directions compared to regular B-DNA was found. For comparison, a molecular dynamics implementation of the approach was also applied to calculate the free energy change associated with bending of A-tract containing DNA, including deformations significantly beyond the optimal curvature. Good agreement with available experimental data was obtained offering an atomic level explanation for stable bending of A-tract containing DNA molecules. The DNA-bending persistence length estimated from the explicit solvent simulations is also in good agreement with experiment whereas the adiabatic mapping calculations with a GB solvent model predict a bending rigidity roughly two times larger.

Nucleosome Positioning Determinants

A previous report demonstrated that one site in a nucleosome assembled onto a synthetic positioning sequence known as Fragment 67 is hypersensitive to permanganate. The site is required for positioning activity and is located 1.5 turns from the dyad, which is a region of high DNA curvature in the nucleosome. Here, the permanganate sensitivity of the nucleosome positioning Fragment 601 was examined in order to expand the dataset of nucleosome sequences containing KMnO(4) hypersensitive sites. The hyperreactive T residue in the six sites detected as well as the one in Fragment 67 and three in the 5 S rDNA positioning sequence were contained within a TA step. Seven of the ten sequences were of the form CTAGPuG or the related sequence TTAAPu. These motifs were also found in the binding sites of several transcriptional regulatory proteins that kink DNA. In order to assess the significance of these sites, the 10 bp positioning determinant in Fragment 67 was removed and replaced by the nine sequences from the 5 S rDNA and Fragment 601. The results demonstrated that these derivative fragments promoted high nucleosome stability and positioning as compared to a control sequence that contained an AT step in place of the TA step. The importance of the TA step was further tested by making single base-pair substitutions in Fragment 67 and the results revealed that stability and positioning activity followed the order: TA>TG>TT>/=TC approximately GG approximately GA approximately AT. Sequences flanking the TA step were also shown to be critical for nucleosome stability and positioning. Nucleosome positioning was restored to near wild-type levels with (CTG)(3), which can form slipped stranded structures and with one base bulges that kink DNA. The results of this study suggest that local DNA structures are important for positioning and that single base-pair changes at these sites could have profound effects on those genomic functions that depend on ordered nucleosomes.

Hydration pattern of A4T4 and T4A4 DNA: A molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.

Shaping the Borrelia burgdorferi genome: crystal structure and binding properties of the DNA-bending protein Hbb

The genome of the Lyme disease-causing spirochete Borrelia burgdorferi encodes only a single polypeptide from the integration host factor (IHF)/HU or 'DNABII' family of nucleoid-associated proteins - Hbb. DNABII proteins induce large bends in DNA and serve as architectural factors in a variety of prokaryotic cellular processes. We have solved the crystal structure of an Hbb-DNA complex in which the DNA is bent by over 180 degrees . We find that like IHF, Hbb relies exclusively on indirect readout to recognize its cognate site. Additional binding studies show that the sequence preferences of Hbb are related to, yet distinct from those of IHF. Defining these binding characteristics may help to uncover additional roles for Hbb in Borrelia DNA metabolism as well as further our understanding of the mechanism of indirect readout.