Curved DNA

Analysis of nucleoid-associated protein-binding regions reveals DNA structural features influencing genome organization in Mycobacterium tuberculosis

Nucleoid-Associated Proteins (NAPs) maintain bacterial nucleoid configuration through their architectural properties of DNA bending, wrapping, and bridging. However, the contribution of DNA structural alterations to DNA-NAP recognition at the genomic scale remains unresolved. Present work dissects the DNA sequence, shape and altered structural preferences at a genomic scale for six NAPs in Mycobacterium tuberculosis. Results suggest narrower minor groove width and higher DNA rigidity are marked for the binding sites of EspR and Lsr2, while mIHF, MtHU and NapM have heterogeneous DNA structural predilections. In contrast, WhiB4-DNA binding sites were characterized by wider minor groove width, highly deformable and less curved DNA. This work provides systematic insight into NAP-mediated genome organization as a function of DNA structural features.

Whole genome nucleosome sequencing identifies novel types of forensic markers in degraded DNA samples

In the case of mass disasters, missing persons and forensic caseworks, highly degraded biological samples are often encountered. It can be a challenge to analyze and interpret the DNA profiles from these samples. Here we provide a new strategy to solve the problem by taking advantage of the intrinsic structural properties of DNA. We have assessed the in vivo positions of more than 35 million putative nucleosome cores in human leukocytes using high-throughput whole genome sequencing, and identified 2,462 single nucleotide variations (SNVs), 128 insertion-deletion polymorphisms (indels). After comparing the sequence reads with 44 STR loci commonly used in forensics, five STRs (TH01, TPOX, D18S51, DYS391, and D10S1248)were matched. We compared these “nucleosome protected STRs” (NPSTRs) with five other non-NPSTRs using mini-STR primer design, real-time PCR, and capillary gel electrophoresis on artificially degraded DNA. Moreover, genotyping performance of the five NPSTRs and five non-NPSTRs was also tested with real casework samples. All results show that loci located in nucleosomes are more likely to be successfully genotyped in degraded samples. In conclusion, after further strict validation, these markers could be incorporated into future forensic and paleontology identification kits, resulting in higher discriminatory power for certain degraded sample types.

Genomic instability in the PARK2 locus is associated with Parkinson's disease

Parkinson’s disease (PD) is a common neurodegenerative disorder affecting mostly elderly people, although there is a group of patients developing so-called early-onset PD (EOPD). Mutations in the PARK2 gene are a common cause of autosomal recessive EOPD. PARK2 belongs to the family of extremely large human genes which are often localised in genomic common fragile sites (CFSs) and exhibit gross instability. PARK2 is located in the centre of FRA6E, the third most mutation-susceptible CFS of the human genome. The gene encompasses a region of 1.3 Mbp and, among its mutations, large rearrangements of single or multiple exons account for around 50 %. We performed an analysis of the PARK2 gene in a group of 344 PD patients with EOPD and classical form of the disease. Copy number changes were first identified using multiplex ligation probe amplification (MLPA), with their ranges characterised by array comparative genomic hybridisation (aCGH). Exact breakpoints were mapped using direct sequencing. Rearrangements were found in eight subjects, including five deletions and three duplications. Rearrangements were mostly non-recurrent and no repetitive sequences or extended homologies were identified in the regions flanking breakpoint junctions. However, in most cases, 1–3 bp microhomologies were present, strongly suggesting that microhomology-mediated mechanisms, specifically non-homologous end joining (NHEJ) and fork stalling and template switching (FoSTeS)/microhomology-mediated break-induced replication (MMBIR), are predominantly involved in the rearrangement processes in this genomic region.

Structural waters in the minor and major grooves of DNA--a major factor governing structural adjustments of the A-T mini-helix

The role of microhydration in structural adjustments of the AT-tract in B-DNA was studied at the B97-D/def2-SV(P) level. The (dA:dT)5 complexes with 10 water molecules in minor and 15 water molecules in major grooves were studied. The obtained network of hydrogen bonds revealed the dependence between the groove width and the types of water patterns. In the minor groove, the following patterns were observed: interstrand one-water bridges similar to that of the Dickerson "water spine" and interstrand two-water bridges. The network of structural waters in the major groove is more diverse than that in the minor groove, which agrees with crystallographic data. As the major groove is wider, it is enriched by water molecules forming two- and three-water bridges. Results suggest the nucleobase-water interactions in both grooves prevent AT-tract twisting and its "collapse" along the minor groove. Whereby, a helix structure with narrow minor and wide major grooves is formed. The structural waters affect the polynucleotide conformation so that it becomes similar to poly(dA)·poly(dT) in fibers and acquires features of the A-tracts in DNA in solution. We suggest that formation of specific water patterns in both grooves is the factor responsible for stabilization of A-tracts with a narrowed minor groove, leading in turn to their strong intrinsic bending in DNA.

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

NAP1-assisted nucleosome assembly on DNA measured in real time by single-molecule magnetic tweezers

While many proteins are involved in the assembly and (re)positioning of nucleosomes, the dynamics of protein-assisted nucleosome formation are not well understood. We study NAP1 (nucleosome assembly protein 1) assisted nucleosome formation at the single-molecule level using magnetic tweezers. This method allows to apply a well-defined stretching force and supercoiling density to a single DNA molecule, and to study in real time the change in linking number, stiffness and length of the DNA during nucleosome formation. We observe a decrease in end-to-end length when NAP1 and core histones (CH) are added to the dsDNA. We characterize the formation of complete nucleosomes by measuring the change in linking number of DNA, which is induced by the NAP1-assisted nucleosome assembly, and which does not occur for non-nucleosomal bound histones H3 and H4. By rotating the magnets, the supercoils formed upon nucleosome assembly are removed and the number of assembled nucleosomes can be counted. We find that the compaction of DNA at low force is about 56 nm per assembled nucleosome. The number of compaction steps and associated change in linking number indicate that NAP1-assisted nucleosome assembly is a two-step process.

Persistence length of human cardiac α-tropomyosin measured by single molecule direct probe microscopy

α-Tropomyosin (αTm) is the predominant tropomyosin isoform in adult human heart and constitutes a major component in Ca²+-regulated systolic contraction of cardiac muscle. We present here the first direct probe images of WT human cardiac αTm by atomic force microscopy, and quantify its mechanical flexibility with three independent analysis methods. Single molecules of bacterially-expressed human cardiac αTm were imaged on poly-lysine coated mica and their contours were analyzed. Analysis of tangent-angle (θ(s)) correlation along molecular contours, second moment of tangent angles (<θ²(s)>), and end-to-end length (L(e-e)) distributions respectively yielded values of persistence length (L(p)) of 41-46 nm, 40-45 nm, and 42-52 nm, corresponding to 1-1.3 molecular contour lengths (L(c)). We also demonstrate that a sufficiently large population, with at least 100 molecules, is required for a reliable L(p) measurement of αTm in single molecule studies. Our estimate that L(p) for αTm is only slightly longer than L(c) is consistent with a previous study showing there is little spread of cooperative activation into near-neighbor regulatory units of cardiac thin filaments. The L(p) determined here for human cardiac αTm perhaps represents an evolutionarily tuned optimum between Ca²+ sensitivity and cooperativity in cardiac thin filaments and likely constitutes an essential parameter for normal function in the human heart.

Differences in DNA curvature-related sequence periodicity between prokaryotic chromosomes and phages, and relationship to chromosomal prophage content

Background

Periodic spacing of A-tracts (short runs of A or T) with the DNA helical period of ~10–11 bp is characteristic of intrinsically bent DNA. In eukaryotes, the DNA bending is related to chromatin structure and nucleosome positioning. However, the physiological role of strong sequence periodicity detected in many prokaryotic genomes is not clear.

Results

We developed measures of intensity and persistency of DNA curvature-related sequence periodicity and applied them to prokaryotic chromosomes and phages. The results indicate that strong periodic signals present in chromosomes are generally absent in phage genomes. Moreover, chromosomes containing prophages are less likely to possess a persistent periodic signal than chromosomes with no prophages.

Conclusions

Absence of DNA curvature-related sequence periodicity in phages could arise from constraints associated with DNA packaging in the viral capsid. Lack of prophages in chromosomes with persistent periodic signal suggests that the sequence periodicity and concomitant DNA curvature could play a role in protecting the chromosomes from integration of phage DNA.

Variety of genomic DNA patterns for nucleosome positioning

Precise positioning of nucleosomes along DNA is important for a variety of gene regulatory processes. Among the factors directing nucleosome positioning, the DNA sequence is highly important. Two main classes of nucleosome positioning sequence (NPS) patterns have previously been described. In the first class, AA, TT, and other WW dinucleotides (where W is A or T) tend to occur together (in-phase) in the major groove of DNA closest to the histone octamer surface, while SS dinucleotides (where S is G or C) are predominantly positioned in the major groove facing outward. In the second class, AA and TT are structurally separated (AA backbone near the histone octamer, and TT backbone further away), but grouped with other RR (where R is purine A or G) and YY (where Y is pyrimidine C or T) dinucleotides. As a result, the RR/YY pattern includes counter-phase AA/TT distributions. We describe here anti-NPS patterns, which are inverse to the conventional NPS patterns: WW runs inverse to SS, and RR inverse to YY. Evidence for the biological relevance of anti-NPS patterns is presented.

Comparative analysis of sequence periodicity among prokaryotic genomes points to differences in nucleoid structure and a relationship to gene expression

Regular spacing of short runs of A or T nucleotides in DNA sequences with a period close to the helical period of the DNA double helix has been associated with intrinsic DNA bending and nucleosome positioning in eukaryotes. Analogous periodic signals were also observed in prokaryotic genomes. While the exact role of this periodicity in prokaryotes is not known, it has been proposed to facilitate the DNA packaging in the prokaryotic nucleoid and/or to promote negative or positive supercoiling. We developed a methodology for assessments of intragenomic heterogeneity of these periodic patterns and applied it in analysis of 1,025 prokaryotic chromosomes. This technique allows more detailed analysis of sequence periodicity than previous methods where sequence periodicity was assessed in an integral form across the whole chromosome. We found that most genomes have the periodic signal confined to several chromosomal segments while most of the chromosome lacks a strong sequence periodicity. Moreover, there are significant differences among different prokaryotes in both the intensity and persistency of sequence periodicity related to DNA curvature. We proffer that the prokaryotic nucleoid consists of relatively rigid sections stabilized by short intrinsically bent DNA segments and characterized by locally strong periodic patterns alternating with regions featuring a weak periodic signal, which presumably permits higher structural flexibility. This model applies to most bacteria and archaea. In genomes with an exceptionally persistent periodic signal, highly expressed genes tend to concentrate in aperiodic sections, suggesting that structural heterogeneity of the nucleoid is related to local differences in transcriptional activity.

Sequence periodicity in nucleosomal DNA and intrinsic curvature

BACKGROUND

Most eukaryotic DNA contained in the nucleus is packaged by wrapping DNA around histone octamers. Histones are ubiquitous and bind most regions of chromosomal DNA. In order to achieve smooth wrapping of the DNA around the histone octamer, the DNA duplex should be able to deform and should possess intrinsic curvature. The deformability of DNA is a result of the non-parallelness of base pair stacks. The stacking interaction between base pairs is sequence dependent. The higher the stacking energy the more rigid the DNA helix, thus it is natural to expect that sequences that are involved in wrapping around the histone octamer should be unstacked and possess intrinsic curvature. Intrinsic curvature has been shown to be dictated by the periodic recurrence of certain dinucleotides. Several genome-wide studies directed towards mapping of nucleosome positions have revealed periodicity associated with certain stretches of sequences. In the current study, these sequences have been analyzed with a view to understand their sequence-dependent structures.

RESULTS

Higher order DNA structures and the distribution of molecular bend loci associated with 146 base nucleosome core DNA sequence from C. elegans and chicken have been analyzed using the theoretical model for DNA curvature. The curvature dispersion calculated by cyclically permuting the sequences revealed that the molecular bend loci were delocalized throughout the nucleosome core region and had varying degrees of intrinsic curvature.

CONCLUSIONS

The higher order structures associated with nucleosomes of C.elegans and chicken calculated from the sequences revealed heterogeneity with respect to the deviation of the DNA axis. The results points to the possibility of context dependent curvature of varying degrees to be associated with nucleosomal DNA.

Fall and rise of satellite repeats in allopolyploids of Nicotiana over c. 5 million years

Allopolyploids represent natural experiments in which DNA sequences from different species are combined into a single nucleus and then coevolve, enabling us to follow the parental genomes, their interactions and evolution over time. Here, we examine the fate of satellite DNA over 5 million yr of divergence in plant genus Nicotiana (family Solanaceae). We isolated subtelomeric, tandemly repeated satellite DNA from Nicotiana diploid and allopolyploid species and analysed patterns of inheritance and divergence by sequence analysis, Southern blot hybridization and fluorescent in situ hybridization (FISH). We observed that parental satellite sequences redistribute around the genome in allopolyploids of Nicotiana section Polydicliae, formed c. 1 million yr ago (Mya), and that new satellite repeats evolved and amplified in section Repandae, which was formed c. 5 Mya. In some cases that process involved the complete replacement of parental satellite sequences. The rate of satellite repeat replacement is faster than theoretical predictions assuming the mechanism involved is unequal recombination and crossing-over. Instead we propose that this mechanism occurs with the deletion of large chromatin blocks and reamplification, perhaps via rolling circle replication.

A statistical thermodynamic approach for predicting the sequence-dependent nucleosome positioning along genomes

Nucleosomes are the fundamental repeating unit of chromatin and constitute the structural building blocks of the eukaryotic genome. The distribution of nucleosomes along the genome is a significant aspect of chromatin structure and influences gene regulation through modulation of DNA accessibility. For this reason, an increasing interest is arising in models capable of predicting the nucleosome positioning along genomes. Toward this goal, we propose a theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence. The model, based on a statistical mechanical approach allows the calculation of the canonical ensemble free energy involved in nucleosome formation. The theoretical free energies were evaluated for about one hundred nucleosome DNA tracts and successfully compared with those obtained in different laboratories with nucleosome competitive reconstitution (correlation coefficient equal to 0.92). We extended these results to the nucleosome positioning along genomes. To test our model, the theoretical nucleosome distribution was compared with that of yeast genome experimentally determined. The results are comparable with those obtained by different authors adopting models based on identifying some recurrent sequence features obtained from the statistical analysis of a very large pool of nucleosomal DNA sequences provided by the positioning maps of genomes.

Upstream curved sequences in E. coli are related to the regulation of transcription initiation

The advancement in Escherichia coli genome research has made the information regarding transcription start sites of many genes available. A study relying on the availability of transcription start locations was performed. The first question addressed was what an average DNA curvature profile upstream of genes would look like when these genes are aligned by transcription start sites in comparison to alignment by translation start sites. Since it was hypothesized that curvature plays a role in transcription regulation, the expectation was that curvature measurements relative to transcription starts, rather than translation, should strengthen the signal. Our study justified this expectation. The second question aimed to clarify the relation between DNA curvature and promoter strength. Through clustering based on DNA curvature profiles along promoter regions, a strong positive correlation between the promoter strength and the curved DNA was found. The third question dealt with dinucleotide periodicity in E. coli to see whether a periodicity pattern specific to promoter regions exists. Such unknown pattern might shed new light on transcription regulation mechanisms in E. coli. A sequence periodicity of about 11 bp is characteristic to the whole E. coli genome, and is especially well-expressed in intergenic regions. Here it was shown that regions of the size of about 100-150 bp centered 70-100 bp upstream to transcription starts carry hidden periodicity with a period of about 10.3 bp.

Restriction landmark genomic scanning (RLGS) spot identification by second generation virtual RLGS in multiple genomes with multiple enzyme combinations

Background

Restriction landmark genomic scanning (RLGS) is one of the most successfully applied methods for the identification of aberrant CpG island hypermethylation in cancer, as well as the identification of tissue specific methylation of CpG islands. However, a limitation to the utility of this method has been the ability to assign specific genomic sequences to RLGS spots, a process commonly referred to as "RLGS spot cloning."

Results

We report the development of a virtual RLGS method (vRLGS) that allows for RLGS spot identification in any sequenced genome and with any enzyme combination. We report significant improvements in predicting DNA fragment migration patterns by incorporating sequence information into the migration models, and demonstrate a median Euclidian distance between actual and predicted spot migration of 0.18 centimeters for the most complex human RLGS pattern. We report the confirmed identification of 795 human and 530 mouse RLGS spots for the most commonly used enzyme combinations. We also developed a method to filter the virtual spots to reduce the number of extra spots seen on a virtual profile for both the mouse and human genomes. We demonstrate use of this filter to simplify spot cloning and to assist in the identification of spots exhibiting tissue-specific methylation.

Conclusion

The new vRLGS system reported here is highly robust for the identification of novel RLGS spots. The migration models developed are not specific to the genome being studied or the enzyme combination being used, making this tool broadly applicable. The identification of hundreds of mouse and human RLGS spot loci confirms the strong bias of RLGS studies to focus on CpG islands and provides a valuable resource to rapidly study their methylation.

Unique Physical Signature of DNA Curvature and Its Implications for Structure and Dynamics

A particularly sensitive birefringence technique is used to analyze a curved DNA fragment with 118 bp and a standard DNA with 119 bp. At salt concentrations from 0.5 to 10 mM, both fragments show the usual negative stationary birefringence and monotonic transients - differences are relatively small. At 100 mM salt the curved DNA shows a positive stationary birefringence and non-monotonic transients with processes having amplitudes of opposite sign, whereas signals of the standard DNA remain as usual. Transients induced by reversal of the field vector indicate the existence of a permanent dipole for the curved DNA. 2-MHz-ac pulses induce a negative stationary birefringence in both DNAs. These results are consistent with calculations on models for curved DNA predicting a quasi-permanent dipole and a positive dichroism/birefringence. The quasi-permanent dipole results from the loss of symmetry in the charge distribution of the curved polyelectrolyte. The appearance of the unique signature of curvature at high salt is mainly due to a strong decrease of the polarizability by about 2 orders of magnitude. The special mode of orientation resulting from the quasi-permanent dipole is expected to contribute to the gel migration anomaly. The time constants of birefringence decay for the curved fragment are shorter than those of the 119 bp fragment by a factor of approximately 1.10 at 0.6 mM salt, whereas this factor is approximately 1.20 at 100 mM Na+. If both fragments were normal DNA with 3.4 A rise per base pair, the factor would be approximately 1.02. At high salt and high electric field strengths the factor increases up to 1.37. The implications for the bending dynamics and the potential to distinguish static from dynamic persistence by field reversal experiments are discussed. The dependence of the curvature on the salt concentration indicated by the time constants is consistent with a clear decrease of the electrophoretic anomaly at decreasing salt concentration.

Sequence-directed mapping of nucleosome positions

A nucleosome DNA sequence probe is designed that combines recently derived RR/YY counter-phase and AA/TT in-phase periodical patterns. A simple nucleosome mapping procedure is introduced for prediction of the nucleosome positions in the sequence of interest, to serve as a guide for experimental studies of the chromatin structure.

Phylogenetic reconstruction of Aegilops section Sitopsis and the evolution of tandem repeats in the diploids and derived wheat polyploids

The evolution of 2 tandemly repeated sequences Spelt1 and Spelt52 was studied in Triticum species representing 2 evolutionary lineages of wheat and in Aegilops sect. Sitopsis, putative donors of their B/G genomes. Using fluorescence in situ hybridization we observed considerable polymorphisms in the hybridization patterns of Spelt1 and Spelt52 repeats between and within Triticum and Aegilops species. Between 2 and 28 subtelomeric sites of Spelt1 probe were detected in Ae. speltoidies, depending on accession. From 8 to 12 Spelt1 subtelomeric sites were observed in species of Timopheevi group (GAt genome), whereas the number of signals in emmer/aestivum accessions was significantly less (from 0 to 6). Hybridization patterns of Spelt52 in Ae. speltoides, Ae. longissima, and Ae. sharonensis were species specific. Subtelomeric sites of Spelt52 repeat were detected only in T. araraticum (T. timopheevii), and their number and chromosomal location varied between accessions. Superimposing copy number data onto our phylogenetic scheme constructed from RAPD data suggests 2 major independent amplifications of Spelt52 and 1 of Spelt1 repeats in Aegilops divergence. It is likely that the Spelt1 amplification took place in the ancient Ae. speltoides before the divergence of polyploid wheats. The Spelt52 repeat was probably amplified in the lineage of Ae. speltoides prior to divergence of the allopolyploid T. timopheevii but after the divergence of T. durum. In a separate amplification event, Spelt52 copy number expanded in the common ancestor of Ae. longissima and Ae. sharonensis.

Sequence Structure of Human Nucleosome DNA

Positional distributions of various dinucleotides in experimentally derived human nucleosome DNA sequences are analyzed. Nucleosome positioning in this species is found to depend largely on GG and CC dinucleotides periodically distributed along the nucleosome DNA sequence, with the period of 10.4 bases. The GG and CC dinucleotides oscillate counterphase, i.e., their respective preferred positions are shifted about a half-period from one another, as it was observed earlier for AA and TT dinucleotides. Other purine-purine and pyrimidine-pyrimidine dinucleotides (RR and YY) display the same periodical and counterphase pattern. The dominance of oscillating GG and CC dinucleotides in human nucleosomes and the contribution of AG(CT), GA(TC), and AA(TT) suggest a general nucleosome DNA sequence pattern - counterphase oscillation of RR and YY dinucleotides. AA and TT dinucleotides, commonly accepted as major players, are only weak contributors in the case of human nucleosomes.

Plant 5S rDNA has multiple alternative nucleosome positions

In plants, 5S ribosomal DNA (5S rDNA) is typically found in hundreds of copies of tandemly arranged units. Nucleotide database searches revealed that the majority of 5S genes (>90%) have repeat lengths that are not simple multiples of a plant nucleosomal unit, ranging in plants from 175-185 bp. To get insight into the chromatin structure, we have determined positions of nucleosomes in the Nicotiana sylvestris and Nicotiana tomentosiformis 5S rDNA units with repeat lengths of about 430 and 645 bp, respectively. Mapping experiments carried out on isolated nucleo somal DNA revealed many (>50) micrococcal nuclease cleavage sites in each class of repeats. Permutation analysis and theoretical computer prediction showed multiple DNA bend sites, mostly located in the nontranscribed spacer region. The distance between bend sites, however, did not correspond to the average spacing of nucleosomes in 5S chromatin (approximately 180 bp). These data indicate that 5S rDNA does not have fixed nucleosomal positioning sites and that units can be wrapped in a number of alternative nucleosome frames. Consequently, accessibility of transcription factors to cognate motifs might vary across the tandem array, potentially influencing gene expression.

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816-3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed 'quasi-regularly' throughout the genome, including both the coding and intergenic sequences. Second, there is a 10-12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in approximately 100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the approximately 100 bp long clusters of the phased A-tracts constitute the 'structural code' for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.

Gene splice sites correlate with nucleosome positions

Gene sequences in the vicinity of splice sites are found to possess dinucleotide periodicities, especially RR and YY, with the period close to the pitch of nucleosome DNA. This confirms previously reported findings about preferential positioning of splice junctions within the nucleosomes. The RR and YY dinucleotides oscillate counter-phase, i.e., their respective preferred positions are shifted about half-period from one another, as it was observed earlier for AA and TT dinucleotides. Species specificity of nucleosome positioning DNA pattern is indicated by the predominant use of the periodical GG(CC) dinucleotides in human and mouse genes, as opposed to predominant AA(TT) dinucleotides in Arabidopsis and C. elegans.

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.

Centromere parC of plasmid R1 is curved

The centromere sequence parC of Escherichia coli low-copy-number plasmid R1 consists of two sets of 11 bp iterated sequences. Here we analysed the intrinsic sequence-directed curvature of parC by its migration anomaly in polyacrylamide gels. The 159 bp long parC is strongly curved with anomaly values (k-factors) close to 2. The properties of the parC curvature agree with those of other curved DNA sequences. parC contains two regions of 5-fold repeated iterons separated by 39 bp. We modified 4 bp within this intermediate sequence so that we could analyse the two 5-fold repeated regions independently. The analysis shows that the two repeat regions are not independently curved parts of parC but that the overall curvature is a property of the whole fragment. Since the centromere sequence of an E.coli plasmid as well as eukaryotic centromere sequences show DNA curvature, we speculate that curvature might be a general property of centromeres.

Sequence periodicity of Escherichia coli is concentrated in intergenic regions

BACKGROUND

Sequence periodicity with a period close to the DNA helical repeat is a very basic genomic property. This genomic feature was demonstrated for many prokaryotic genomes. The Escherichia coli sequences display the period close to 11 base pairs.

RESULTS

Here we demonstrate that practically only ApA/TpT dinucleotides contribute to overall dinucleotide periodicity in Escherichia coli. The noncoding sequences reveal this periodicity much more prominently compared to protein-coding sequences. The sequence periodicity of ApC/GpT, ApT and GpC dinucleotides along the Escherichia coli K-12 is found to be located as well mainly within the intergenic regions.

CONCLUSIONS

The observed concentration of the dinucleotide sequence periodicity in the intergenic regions of E. coli suggests that the periodicity is a typical property of prokaryotic intergenic regions. We suppose that this preferential distribution of dinucleotide periodicity serves many biological functions; first of all, the regulation of transcription.

Structure of DNA Hexamer Sequence d-CGATCG by Two-dimensional Nuclear Magnetic Resonance Spectroscopy and Restrained Molecular Dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.

Structure-specific DNA-binding proteins as the foundation for three-dimensional chromatin organization

Any functions of tandem repetitive sequences need proteins that specifically bind to them. Telomere-binding TRF2/MTBP attaches telomeres to the nuclear envelope in interphase due to its rod-domain-like motif. Interphase nuclei organized as a number of sponge-like ruffly round chromosome territories that could be rotated from outside. SAF-A/hnRNP-U and p68-helicase are proteins suitable to do that. Their location in the interchromosome territory space, ATPase domains, and the ability to be bound by satellite DNAs (satDNA) make them part of the wires used to help chromosome territory rotates. In case of active transcription p68-helicase can be involved in the formation of local "gene expression matrices" and due to its satDNA-binding specificity cause the rearrangement of the local chromosome territory. The marks of chromatin rearrangement, which have to be heritable, could be provided by SAF-A/hnRNP-U. During telophase unfolding the proper chromatin arrangement is restored according to these marks. The structural specificity of both proteins to the satDNAs provides a regulative but relatively stable mode of binding. The structural specificity of protein binding could help to find the "magic" centromeric sequence. With future investigations of proteins with the structural specificity of binding during early embryogenesis, when heterochromatin formation goes on, the molecular mechanisms of the "gene gating" hypothesis (Blobel, 1985) will be confirmed.

Sequence-dependent pausing of single lambda exonuclease molecules

Lambda exonuclease processively degrades one strand of duplex DNA, moving 5'-to-3' in an ATP-independent fashion. When examined at the single-molecule level, the speeds of digestion were nearly constant at 4 nanometers per second (12 nucleotides per second), interspersed with pauses of variable duration. Long pauses, occurring at stereotypical locations, were strand-specific and sequence-dependent. Pause duration and probability varied widely. The strongest pause, GGCGAT TCT, was identified by gel electrophoresis. Correlating single-molecule dwell positions with sequence independently identified the motif GGCGA. This sequence is found in the left lambda cohesive end, where exonuclease inhibition may contribute to the reduced recombination efficiency at that end.

Analysis of DNA curvature distribution in mycobacterial promoters using theoretical models

In this paper, 125 different mycobacterial promoters are analyzed for their DNA curvature distribution using several di- and tri-nucleotide dependent models of DNA curvature. Different models give similar behavior and therefore qualitative validation of the results. Mycobacterial promoters resembling the E. coli sigma(70) type have almost 81% (85%) sequences having medium and high curvature profiles using dinucleotide-dependent models. Non-E. coli sigma(70) type mycobacterial promoters have comparatively higher percent of low curvature profiles. Very few extended -10 promoters have low curvature profiles. Mycobacterial promoters having A(n)T(m) (n+m > or =3) tract in the upstream region of -35 box and repeated in phase with each other have high curvature profiles. M. smegmatis promoters have high curvature profiles compared to M. tuberculosis promoters.

Time-resolved fluorescence resonance energy transfer studies of DNA bending in double-stranded oligonucleotides and in DNA-protein complexes

Time-resolved Förster resonance energy transfer (trFRET) has been used to obtain interdye distance distributions. These distributions give the most probable distance as well as a parameter, sigma, that characterize the width of the distribution. This latter parameter contains information not only on the flexibility of the dyes tethered to macromolecules, but on the flexibility of the macromolecules. Both the most probable interdye distance as well as sigma provide insight into DNA static bending and DNA flexibility. Time-resolved fluorescence anisotropy and static anisotropy measurements can be combined to provide a measure of the cone angle within which the tethered dyes appear to wobble. When this motion is an order of magnitude faster than the average lifetime that characterizes transfer, an average value of the dipolar orientational parameter kappa2 can be calculated for various mutual dye orientations. The resulting kappa2 distribution is very much narrower than the limiting values of 0 and 4, allowing more precise distances and distance changes to be determined. Static and time-resolved fluorescence data can be combined to constrain the analyses of DNA-protein kinetics to provide thermodynamic parameters for binding and for conformational changes along a reaction coordinate. The parameter sigma can be used to model multiple DNA-protein complexes with varying DNA bend angles in a global fitting of trFRET data. Such a global fitting approach has shown how the range of bends in single base DNA variants, when bound by the TATA binding protein (TBP), can be understood in terms of two limiting forms. Time-resolved FRET, combined with steady-state FRET, can be used to show not only how osmolytes affect the binding of DNA to proteins, but also how DNA bending depends on osmolyte concentration in the DNA-protein complexes.

Contribution of the intrinsic curvature to measured DNA persistence length

The persistence length of DNA, a, depends both on the intrinsic curvature of the double helix and on the thermal fluctuations of the angles between adjacent base-pairs. We have evaluated two contributions to the value of a by comparing measured values of a for DNA containing a generic sequence and for an "intrinsically straight" DNA. In each 10 bp segment of the intrinsically straight DNA an initial sequence of five bases is repeated in the sequence of the second five bases, so any bends in the first half of the segment are compensated by bends in the opposite direction in the second half. The value of a for the latter DNA depends, to a good approximation, on thermal fluctuations only; there is no intrinsic curvature. The values of a were obtained from measurements of the cyclization efficiency for short DNA fragments, about 200 bp in length. This method determines the persistence length of DNA with exceptional accuracy, due to the very strong dependence of the cyclization efficiency of short fragments on the value of a. We find that the values of a for the two types of DNA fragment are very close and conclude that the contribution of the intrinsic curvature to a is at least 20 times smaller than the contribution of thermal fluctuations. The relationship between this result and the angles between adjacent base-pairs, which specify the intrinsic curvature, is analyzed.

Preferential counterion binding to A-tract DNA oligomers

The free solution mobility of four 20 bp DNA oligomers, with and without A-tracts, has been measured by capillary electrophoresis in Tris-acetate buffer, to test the hypothesis that site-specific binding of monovalent counterions can occur in the narrow minor groove of A-tract DNAs. Preferential counterion binding has been proposed to cause A-tract bending because of asymmetric charge neutralization and collapse of the helix backbone toward the minor groove. Preferential counterion binding in A-tract DNAs should be manifested by a decrease in the electrophoretic mobility observed in free solution, compared to that of non-A-tract DNAs of the same size. Of the four sequences studied here, the slowest absolute mobility, indicative of the greatest counterion binding, was observed for a 20 bp oligomer containing two runs of A3T3 in phase with the helix repeat. A 20-mer containing phased CACA sequences migrated with the fastest mobility; 20-mers containing phased A5 tracts or phased runs of T3A3 migrated with intermediate mobilities. Very similar mobility differences were observed when 1-20 mM NaCl was added to the buffer. The results suggest that preferential counterion binding occurs in A-tract DNAs, especially those containing the AnTn sequence motif.

A refined prediction method for gel retardation of DNA oligonucleotides from dinucleotide step parameters: reconciliation of DNA bending models with crystal structure data

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinulcleotide step parameters are described. The method is formulated using the Babcock-Olson equations for base pair step geometry (1) and employs Monte Carlo simulated annealing for parameter optimization against experimental data. The refined base pair step parameters define a stuctural construct which, when the width of observed parameter distributions is taken into account, is consistent with the results of DNA oligonucleotide crystal structures. The predictive power of the method is demonstrated and tested via comparisons with DNA bending data on sets of sequences not included in the training set, including A-tracts with and without periodic helix phasing, phased A4T4 and T4A4 motifs, a sequence with a phased GGGCCC motif, some "unconventional" helix phasing sequences, and three short fragments of kinetoplast DNA from Crithidia fasiculata that exhibit significantly different behavior on non-denaturing polyacrylamide gels. The nature of the structural construct produced by the methodology is discussed with respect to static and dynamic models of structure and representations of bending and bendability. An independent theoretical account of sequence dependent chemical footprinting results is provided. Detailed analysis of sequences with A-tract induced axis bending forms the basis for a critical discussion of the applicability of wedge models,junction models and non A-tract, general sequence models for understanding the origin of DNA curvature at the molecular level.

Global structure and mechanical properties of a 10-bp nucleosome positioning motif

The method of DNA cyclization kinetics reveals special properties of the TATAAACGCC sequence motif found in DNA sequences that have high affinity for core histones. Replacement of 30 bp of generic DNA by three 10-bp repeats of the motif in small cyclization constructs increases cyclization rates by two orders of magnitude. We document a 13 degrees bend in the motif and characterize the direction of curvature. The bending force constant is smaller by nearly 2-fold and there is a 35% decrease in the twist modulus, relative to generic DNA. These features are the likely source of the high affinity for bending around core histones to form nucleosomes. Our results establish a protocol for determination of the ensemble-averaged global solution structure and mechanical properties of any approximately 10-bp DNA sequence element of interest, providing information complementary to that from NMR and crystallographic structural studies.

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.

Effects of local transcription and H-NS on inversion of the fim switch of Escherichia coli

The fim switch of Escherichia coli is responsible for phase-variable expression of type 1 fimbriae. Switching in the ON-to-OFF and OFF-to-ON directions is promoted by the FimB recombinase, while the FimE recombinase directs switching predominantly in the ON-to-OFF direction. The effects of local promoter activity and the H-NS nucleoid-associated protein on inversion of the switch were assessed. In contrast to FimB-mediated inversion, inversion of the switch by the FimE recombinase was unaffected by the H-NS status of the cell. Transcription towards the switch from within a translationally inactivated fimE gene was found to bias the switch strongly in the OFF direction, creating a FimE+-like phenotype in the absence of the FimE protein. This biasing was H-NS dependent and was also contingent on transcription from within the switch. These data show that local transcription and a nucleoid-associated protein both contribute to the modulation of a site-specific recombination event on the bacterial chromosome.

Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure

Single chicken erythrocyte chromatin fibers were stretched and released at room temperature with force-measuring laser tweezers. In low ionic strength, the stretch-release curves reveal a process of continuous deformation with little or no internucleosomal attraction. A persistence length of 30 nm and a stretch modulus of approximately 5 pN is determined for the fibers. At forces of 20 pN and higher, the fibers are modified irreversibly, probably through the mechanical removal of the histone cores from native chromatin. In 40-150 mM NaCl, a distinctive condensation-decondensation transition appears between 5 and 6 pN, corresponding to an internucleosomal attraction energy of approximately 2.0 kcal/mol per nucleosome. Thus, in physiological ionic strength the fibers possess a dynamic structure in which the fiber locally interconverting between "open" and "closed" states because of thermal fluctuations.

Transcriptional activation by the myb proteins requires a specific local promoter structure

The biological effects of the cellular c-Myb and the viral v-Myb proteins are strikingly different. While c-Myb is indispensable for normal hematopoiesis, v-Myb induces acute leukemia. The v-Myb DNA-binding domain (DBD) differs from that of c-Myb mainly by deletion of the first of three repeats which correlates with efficient oncogenic transformation and a decrease in DNA-binding activity. To investigate the difference in DNA-binding and transcriptional activation, oligonucleotide selection and electrophoretic mobility shift assays were employed. The v-Myb DBD (R2R3) shows an intrinsic DNA-binding specificity for an AT-rich downstream extension of the Myb recognition element (MRE) PyAAC(T)/(G)G for efficient binding to this site, whereas R1 within the c-Myb DBD allows for more flexibility for this downstream extension. Therefore, due to the presence of repeat R1, c-Myb can bind to a greater number of target sites. The intrinsic DNA-binding specificity of R2R3 is further supported with the results from in vivo transcriptional activation experiments which demonstrated that both the v-Myb and c-Myb DBDs require an extension of the MRE (motif #1) by a downstream T-stretch (motif #2) for full activity. Surprisingly, the T-stretch improves binding when present on either strand, but is required on a specific strand for transcriptional activation.

DNA multimode interaction with berenil and pentamidine; double helix stiffening, unbending and bending

The antitrypanosomal drugs berenil (Ber) and pentamidine (Pm) preferentially bind to DNA in the minor groove of A.T-rich domains. The properties of A.T clusters are essential for sequence-mediated helix bending. Groove binding drugs locally stiffen the DNA helix but may also change intrinsic helix bends or may bend straight DNA. Ligand binding to randomly distributed sites alters the apparent DNA persistence length, a. Criteria permit the distinction of the underlying mechanism(s). Helix bends, if phased with the helix screw, however, generate solenoidal superhelix components mediating an apparent change of the hydrodynamically effective DNA contour length, L. The measurement of relative changes of both, a and L, as induced by Ber or Pm is performed by titration rotational viscometry. The determination of the two quantities requires two independent measurements: the relative change of DNA intrinsic viscosity, deltay, for short (tending to rod-like) DNA molecules and for comparably long (almost coil-like) ones as a function of r, the bound drug molecules per DNA-P, and this under conditions effectively excluding intramolecular DNA-DNA crosslinking effects. At least at r< or =0.05 and < or =0.03, respectively, the two drugs virtually bind completely to a eukaryotic DNA. r ranges of different drug binding strength and, concomitantly, of different specific conformational response, could be resolved. They represent (sub)modes of different DNA sequences... Whereas the mode-specific elongation effects are fairly similar for both systems, there are pronounced quantitative differences in the relative change of DNA persistence length. The sites of highest Ber-binding strength are correlated to unbent alternating helical A.T segments followed by bent and by less bent or unbent dAn.dTn tracts straightened on Ber-binding. For Pm-DNA interaction the ligand bends the sites of highest Pm affinity. Generally, ligand induced and sequence mediated local DNA-bend removal or DNA bending, as observed for several modes of interaction with A.T rich DNA, are considered to be of gene regulatory relevance.

Protein-DNA interactions: A structural analysis

A detailed analysis of the DNA-binding sites of 26 proteins is presented using data from the Nucleic Acid Database (NDB) and the Protein Data Bank (PDB). Chemical and physical properties of the protein-DNA interface, such as polarity, size, shape, and packing, were analysed. The DNA-binding sites shared common features, comprising many discontinuous sequence segments forming hydrophilic surfaces capable of direct and water-mediated hydrogen bonds. These interface sites were compared to those of protein-protein binding sites, revealing them to be more polar, with many more intermolecular hydrogen bonds and buried water molecules than the protein-protein interface sites. By looking at the number and positioning of protein residue-DNA base interactions in a series of interaction footprints, three modes of DNA binding were identified (single-headed, double-headed and enveloping). Six of the eight enzymes in the data set bound in the enveloping mode, with the protein presenting a large interface area effectively wrapped around the DNA.A comparison of structural parameters of the DNA revealed that some values for the bound DNA (including twist, slide and roll) were intermediate of those observed for the unbound B-DNA and A-DNA. The distortion of bound DNA was evaluated by calculating a root-mean-square deviation on fitting to a canonical B-DNA structure. Major distortions were commonly caused by specific kinks in the DNA sequence, some resulting in the overall bending of the helix. The helix bending affected the dimensions of the grooves in the DNA, allowing the binding of protein elements that would otherwise be unable to make contact. From this structural analysis a preliminary set of rules that govern the bending of the DNA in protein-DNA complexes, are proposed.

DNA bendability--a novel feature in E. coli promoter recognition

The distribution of deformable base-pair steps in the structure of bacterial promoters is analyzed with respect to their possible structural and functional role. A regular positioning of TA and TG stacks is detected with the best fit period 5.6 bp. This value is interpreted as a half of the sequence period 11.2 bp, somewhat higher than the structural helical repeat of B-DNA (10.55 bp). The difference, +0.65 bp, suggests a sequence-dependent helical writhe of the promoter DNA--a right-handed superhelix. Apparently, to favour rotational setting of DNA on the surface of RNA polymerase the flexible steps deformable largely towards the grooves, follow the half-period spacing. Such rotational setting is consistent with the DNase I footprinting data. Periodical distribution of deformable base-pair stacks shows negative correlation with the presence of -35 canonical hexamer, suggesting the functional significance of this novel element for promoter recognition. The RNA polymerase--DNA recognition is discussed as interaction of distributional type that involves many elements of different nature which are in partially compensatory relations.