Ranking of Prokaryotic Genomes Based on Maximization of Sortedness of Gene Lengths.

How variations of gene lengths (some genes become longer than their predecessors, while other genes become shorter and the sizes of these factions are randomly different from organism to organism) depend on organismal evolution and adaptation is still an open question. We propose to rank the genomes according to lengths of their genes, and then find association between the genome rank and variousproperties, such as growth temperature, nucleotide composition, and pathogenicity. This approach reveals evolutionary driving factors. The main purpose of this study is to test effectiveness and robustness of several ranking methods. The selected method of evaluation is measuring of overall sortedness of the data. We have demonstrated that all considered methods give consistent results and Bubble Sort and Simulated Annealing achieve the highest sortedness. Also, Bubble Sort is considerably faster than the Simulated Annealing method.

Activity of Serratia plymuthica IC1270 gene chiA promoter region in Escherichia coli mutants deficient in global regulators of transcription.

To study the regulation of expression of the Serratia plymuthica gene chiA encoding a 58-kDa endochitinase, its 586-bp-long upstream regulatory region was cloned, sequenced and fused to a promoterless lac operon in phage lambdaRS45 to obtain a single-copy transcriptional fusion (P F1chiA )-lac in lysogens of Escherichia coli wild-type strains or their mutants deficient in various global regulators of transcription. The level of P F1chiA -lac expression increased about 20- and 90-fold, respectively, in E. coli K12 Deltahns and double Deltahns stpA mutants deficient in H-NS, and in both H-NS and StpA DNA-binding histone-like proteins, as compared to levels in the wild-type strain. In a Deltalrp mutant deficient in the leucine-responsive transcriptional regulator Lrp, the level of P F1chiA -lac expression increased only up to threefold, whereas even smaller differences relative to the wild-type strain were observed in rpoS and Deltacrp mutants deficient in the sigmaS subunit of RNA polymerase and catabolite-repression protein (CRP), respectively. Deletion of the inverted-repeat sequences and curved DNA regions located in the upstream region of chiA essentially did not influence strain IC1270's chiA promoter activity in E. coli .

New elements of the termination of transcription in prokaryotes.

We report here that phased runs of adenines and thymines are very frequent in the neighborhood of 3' of the coding regions of Escherichia coli and Bacillus subtilis. These findings suggest that the DNA curvature could affect transcription termination either directly, through contacts with RNA polymerase, or indirectly, via contacts with some regulatory proteins.

Ecologic genomics of DNA: upstream bending in prokaryotic promoters.

After our analysis of the distribution of predicted intrinsic curvature along all available complete prokaryotic genomes, the genomes were divided into two groups. Curvature distribution in all prokaryotes of the first group indicated a substantial fraction of promoters characterized by intrinsic DNA curvature located within or upstream of the promoter region. We did not find this peculiar DNA curvature distribution in prokaryotes in the second group. Remarkably, all bacteria of the first group were mesophilic, whereas many prokaryotes of the second group were hyperthermophilic. We hypothesize that DNA curvature plays a biologic role in gene regulation in mesophilic as opposed to hyperthermophilic prokaryotes, i.e., DNA curvature presumably has a functional adaptive significance determined by temperature selection.

Sequence complexity and DNA curvature.

A linguistic complexity measure was applied to the complete genomes of HIV-1, Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Mycoplasma genitalium, and to long human and yeast genomic fragments. Complexity values averaged over entire genomic sequences were compared, as were predicted average values of intrinsic DNA curvature. We found that both the most curved and the least complex fragments are located preferentially in non-coding parts of the genome. Analysis of location of the most curved and the simplest regions in bacteria showed that the low-complexity segments are preferentially located in close proximity to the highly curved sequences, which are, in turn, placed from 100 to 200 bases upstream to the start of the nearest coding sequence. We conclude that the parallel analysis of sequence complexity and DNA curvature might provide important information about sequence-structure-function relationship in genomes.

Environmental influences on DNA curvature.

DNA curvature plays an important role in many biological processes. To study environmental influences on DNA curvature we compared the anomalous migration on polyacrylamide gels of ligation ladders of 11 specifically-designed oligonucleotides. At low temperatures (25 degrees C and below) most of the sequences exhibited a degree of anomalous migration. Increased temperature had a significant effect on the anomalous migration (curvature) of some sequences but limited effects on others; at 50 degrees C only 1 sequence migrated anomalously. Mg2+ had a strong influence on the migration of certain sequences, whilst spermine enhanced the anomalous migration of a different set of sequences. Sequences with a GGC motif exhibited greater curvature than predicted by the presently-used angles for the nearest-neighbour wedge model and are especially sensitive to Mg2+. The data have implications for models for DNA curvature and for environmentally-sensitive DNA conformations in the regulation of gene expression.

Curved DNA in promoter sequences.

Sequence-dependent DNA curvature is known to play an important role in initiation of transcription of many genes. We compared the distribution of predicted intrinsic curvature of Escherichia coli and human promoter sequences with the distribution of curvature of randomly selected coding and non-coding fragments from these organisms. Different methods of curvature calculation were found to yield mostly similar overall tendencies of DNA curvature in all groups of sequences. According to all methods of calculation, E. coli promoters were found to be more curved than coding sequences from the same genome and random sequences with the same nucleotide composition. By contrast, the average curvature of human promoter sequences was only marginally greater than the average curvature of human coding sequences. Non-coding intron sequences were found to be the most curved of the human sequences examined. Based on these observations, we hypothesize about the role of DNA curvature in promoter sequences.

Enhancement of the nucleosomal pattern in sequences of lower complexity.

Intuitively, the complexity of a given DNA sequence is related to the number of various superimposed biological messages it contains. Here we assess the expectation that in nucleosome DNA sequences of lower linguistic complexity, the nucleosome DNA positioning pattern would be more pronounced than in those of higher linguistic complexity. The nucleosome DNA positioning pattern is one of the weakest (highly degenerate) sequence patterns. It has been extracted recently by specially designed multiple alignment procedures. We applied the most sensitive of these procedures to nearly equal subsets of a nucleosome database separated according to linguistic complexity. The pattern extracted from the subset of the simpler nucleosome sequences not only possesses all major attributes of the known nucleosomal pattern, but is substantially stronger with respect to amplitude in comparison with the total database. This result constitutes the first demonstration that a weak pattern can be significantly enhanced by selective treatment of a lower complexity subset of the sequence ensemble under consideration.

Applicability of the multiple alignment algorithm for detection of weak patterns: periodically distributed DNA pattern as a study case.

MOTIVATION: A nucleosome DNA positioning pattern is known to be one of the weakest (highly degenerated) patterns. The alignment procedure that has been developed recently for the extraction of such a pattern is based on a statistical matching of the sequences, and its success depends on the pattern/background ratio in the individual sequences and in the generated pattern. The heuristic nature of the method and distinctive properties of the pattern bring up the question of efficiency and sensitivity in the procedure. This paper presents a method of verification for this multiple sequence alignment algorithm. RESULTS: To verify the applicability of the multiple alignment approach, we constructed a set of sequences carrying the hidden pattern. The pattern was presented by weak ('signal') oscillations of occurrences of AA and TT dinucleotides along otherwise random sequences. Only a few dinucleotides of any given 145 base long sequence would correspond to the signal, appearing in about the same phase within the simulated periodic pattern. The novelty of our simulation approach is that we simulated a database as a whole, as opposed to simulating each sequence separately. The correlation between the hidden pattern and a sequence from the database is negligible on average, but our statistical multicycle alignment procedure produced the pattern with attributes very close to the simulated ones. The accuracy of the procedure was tested and calibrated. The presence in a typical sequence of as little as three dinucleotides corresponding to the signal is sufficient to generate (detect) the pattern hidden in a collection of 204 sequences.

Nucleosome DNA sequence pattern revealed by multiple alignment of experimentally mapped sequences.

Five different algorithms have been applied for detecting DNA sequence pattern hidden in 204 DNA sequences collected from the literature which are experimentally found to be involved in nucleosome formation. Each algorithm was used to perform a multiple alignment of the nucleosome DNA sequences within the window 145 nt, the size of a nucleosome core DNA. From these alignments five pairs of AA and TT dinucleotide positional frequency distributions have been computed. The frequency profiles calculated by different algorithms are rather different due to substantial noise. They, however, share several important features. Both AA and TT dinucleotide positional frequencies display periodicity with the period of 10.3(+/- 0.2) bases. TT dinucleotides appear to be distributed symmetrically relative to AA dinucleotides of the same DNA strand, with the center of symmetry at the midpoint of the nucleosome core DNA. The phase shift between the AA and TT patterns is about 6 bp. Superposition of the five pairs of the AA (TT) positional frequency profiles has produced the refined pattern, with the above features well pronounced. An interesting novel feature of the pattern is an absence of central peaks in the periodical AA and TT distributions. This may indicate that the central section of nucleosome DNA, 15 bp around the dyad axis of the nucleosome, is not bent. Positional distributions of other dinucleotides were not found in this study to be as informative as the ones for AA and TT.

CURVATURE: software for the analysis of curved DNA.

Software is presented to plot the sequence-dependent spatial trajectory of the DNA double helix and/or distribution of curvature along the DNA molecule. The nearest-neighbor wedge model is implemented to calculate overall DNA path using local helix parameters: helix twist angle, wedge (deflection) angle and direction (of deflection) angle. The procedures described proved to be very convenient as tools for investigation of a relationship between overall DNA curvature and its gel electrophoretic mobility. All parameters of the model had been estimated from experimental data. Using these wedge parameters the program takes, as input, any DNA sequence and calculates the likely degree of curvature at each point along the molecule. This information is displayed both graphically and in the form of simplified representations of curved double helices. The Software, CURVATURE, can thus be used to investigate possible roles of curvature in modulation of gene expression and for location of curved portions of DNA, which may play an important role in sequence-specific protein--DNA interactions.

Preferred positions of AA and TT dinucleotides in aligned nucleosomal DNA sequences.

Multiple alignment of 118 nucleosomal DNA sequences by maximizing simultaneously match of AA dinucleotides and match of TT dinucleotides results in a pattern of the dinucleotide distributions which is characteristic of the nucleosomal DNA sequences. The AA dinucleotides are found to be distributed symmetrically relative to the TT dinucleotide distribution, around the middle point of the nucleosomal DNA sequence. The distances between major peaks of the distributions are multiples of about 10.4 bases. The peaks of the TT distribution are shifted by 6 bases downstream from the peaks of the AA distribution.

Curved DNA without A-A: experimental estimation of all 16 DNA wedge angles.

The principal sequence feature responsible for intrinsic DNA curvature is generally assumed to be runs of adenines. However, according to the wedge model of DNA curvature, each dinucleotide step is associated with a characteristic deflection of the local helix axis. Thus, an important test of a more general view of sequence-dependent DNA curvature is whether sequence elements other than A-A cause the DNA axis to deflect. To address this question, we have applied the wedge model to a large body of experimental data. The axial path of DNA can be described at each step by three Eulerian angles: the helical twist, the deflection angle (wedge angle), and the direction of the deflection. Circularization and gel electrophoretic mobility data on 54 synthetic DNA fragments, both from other laboratories and from our own, were used to compare the theoretical predictions of the wedge model with experiment. By minimizing misfit between calculated and observed DNA curvature, we have found that the stacks AG/CT, CG/CG, GA/TC, and GC/GC, in addition to AA/TT, have large wedge values. We have also synthesized seven sequences without AA/TT elements but with these other wedges correctly phased to cause appreciable predicted curvature. All appear curved as demonstrated by anomalous gel mobilities. The full set of 16 roll and tilt wedge angles is estimated and, together with the known 10 helical twists, these allow prediction of the general sequence-dependent trajectory of the DNA axis.

Sequence-dependent kinks induced in curved DNA.

In certain curved DNA fragments without AA dinucleotides, the gel retardation anomaly associated with curvature passes through a maximum with fragment length, indicating length (and electric field) dependent structural transitions in the DNA. We suggest that thermally induced stereochemical kinks in DNA are stabilized in the gel, thus relieving the effects of curvature. These kinks are shown to occur specifically at CA/TG and TA/TA stacks. Other physical and biological evidence points to frequent structural dislocations at CA and TA steps. These reversible sequence dependent kinks may therefore represent a novel class of structural protein-DNA recognition elements.