Multi-alphabet consensus algorithm for identification of low specificity protein-DNA interactions

Analysis and development of the computer methods of nucleosome localization on DNA fragments with different AT-content

Trinucleotide parameter sets published previously were used for the development of the predictive method for the determining the nucleosome positions along the DNA. The choice of the type of parameter sets used depends upon AT-content of the fragment. Some limitations are imposed on these predictions due to the presence of A(n), T(n) tracts (in our case n>5 or =5) within the 145 bp fragment leading to the displacement or even the prohibition for the corresponding site to be occupied by nucleosomes. The predicted nucleosome positioning site with the large potential may influence on the choice of the proximal nucleosome positions with the weaker bending potentials as is revealed by the comparison with the micrococcal nuclease digestion map. Trinucleotide methods may be considered as advantageous in the comparison with the dinucleotide ones.

In silico approaches reveal the potential for DNA sequence-dependent histone octamer affinity to influence chromatin structure in vivo

Nucleosome positioning signals embedded within the DNA sequence have the potential to influence the detailed structure of the higher-order chromatin fibre. In two previous studies of long stretches of DNA, encompassing the chicken beta-globin and ovine beta-lactoglobulin genes, respectively, we mapped the relative affinity of every site for the core histone octamer. In both cases a periodic arrangement of the in vitro positioning sites suggests that they might influence the folding of a nucleosome chain into higher-order structure; this hypothesis was borne out in the case of the beta-lactoglobulin gene, where the distribution of the in vitro positioning sites is related to the positions nucleosomes actually occupy in sheep liver cells. Here, we have exploited the in vitro nucleosome positioning datasets to simulate nucleosomal organisation using in silico approaches. We use the high-resolution, quantitative positioning maps to define a one-dimensional positioning energy lattice, which can be populated with a defined number of nucleosomes. Monte Carlo techniques are employed to simulate the behaviour of the model at equilibrium to produce a set of configurations, which provide a probability-based occupancy map. Employing a variety of techniques we show that the occupancy maps are a sensitive function of the histone octamer density (nucleosome repeat length) and find that a minimal change in this property can produce dramatic localised changes in structure. Although simulations generally give rise to regular periodic nucleosomal arrangements, they often show octamer density-dependent discontinuities, which tend to co-localise with sequences that adopt distinctive chromatin structure in vivo. Furthermore, the overall organisation of simulated chromatin structures are more closely related to the situation in vivo than is the original in vitro positioning data, particularly at a nucleosome density corresponding to the in vivo state. Although our model is simplified, we argue that it provides a unique insight into the influence that DNA sequence can have in determining chromatin structure and could serve as a useful basis for the incorporation of other parameters.

Inferring regulatory elements from a whole genome. An analysis of Helicobacter pylori sigma(80) family of promoter signals

Helicobacter pylori is adapted to life in a unique niche, the gastric epithelium of primates. Its promoters may therefore be different from those of other bacteria. Here, we determine motifs possibly involved in the recognition of such promoter sequences by the RNA polymerase using a new motif identification method. An important feature of this method is that the motifs are sought with the least possible assumptions about what they may look like. The method starts by considering the whole genome of H. pylori and attempts to infer directly from it a description for a family of promoters. Thus, this approach differs from searching for such promoters with a previously established description. The two algorithms are based on the idea of inferring motifs by flexibly comparing words in the sequences with an external object, instead of between themselves. The first algorithm infers single motifs, the second a combination of two motifs separated from one another by strictly defined, sterically constrained distances. Besides independently finding motifs known to be present in other bacteria, such as the Shine-Dalgarno sequence and the TATA-box, this approach suggests the existence in H. pylori of a new, combined motif, TTAAGC, followed optimally 21 bp downstream by TATAAT. Between these two motifs, there is in some cases another, TTTTAA or, less frequently, a repetition of TTAAGC separated optimally from the TATA-box by 12 bp. The combined motif TTAAGCx(21+/-2)TATAAT is present with no errors immediately upstream from the only two copies of the ribosomal 23 S-5 S RNA genes in H. pylori, and with one error upstream from the only two copies of the ribosomal 16 S RNA genes. The operons of both ribosomal RNA molecules are strongly expressed, representing an encouraging sign of the pertinence of the motifs found by the algorithms. In 25 cases out of a possible 30, the combined motif is found with no more than three substitutions immediately upstream from ribosomal proteins, or operons containing a ribosomal protein. This is roughly the same frequency of occurrence as for TTGACAx(15-19)TATAAT (with the same maximum number of substitutions allowed) described as being the sigma(70 )promoter sequence consensus in Bacillus subtilis and Escherichia coli. The frequency of occurrence of the new motif obtained, TTAAGCx(19-23)TATAAT, remains high when all protein genes in H. pylori are considered, as is the case for the TTGACAx(15-19)TATAAT motif in B. subtilis but not in E. coli.

Promoter sequences and algorithmical methods for identifying them

This paper presents a survey of currently available mathematical models and algorithmical methods for trying to identify promoter sequences. The methods concern both searching in a genome for a previously defined consensus and extracting a consensus from a set of sequences. Such methods were often tailored for either eukaryotes or prokaryotes although this does not preclude use of the same method for both types of organisms. The survey therefore covers all methods; however, emphasis is placed on prokaryotic promoter sequence identification. Illustrative applications of the main extracting algorithms are given for three bacteria.

Recognition of regulatory sites by genomic comparison

Availability of complete bacterial genomes opens the way to the comparative approach to the recognition of transcription regulatory sites. Assumption of regulon conservation in conjunction with profile analysis provides two lines of independent evidence making it possible to make highly specific predictions. Recently this approach was used to analyze several regulons in eubacteria and archaebacteria. The present review covers recent advances in the comparative analysis of transcriptional regulation in prokaryotes and phylogenetic fingerprinting techniques in eukaryotes, and describes the emerging patterns of the evolution of regulatory systems.

Periodical distribution of transcription factor sites in promoter regions and connection with chromatin structure

Nucleosomes regulate transcriptional initiation when positioned in the promoter area. This may require the transcription factor (TF) sites to be correlated with the nucleosome positions and phased on the nucleosome surface. If this is the case, one would expect a periodical distribution of TF sites in the vicinity of promoters, with the nucleosomal period of 10.1-10.5 bp. We examined the distributions of putative binding sites of 323 different TFs along 1, 057 sequences of the Eukaryotic Promoter Database (release 50) [Cavin Perier, R., Junier, T. & Bucher, P. (1998) Nucleic Acids Res. 26, 353-357] and of 218 TFs on 673 sequences of the Lead Exon Database of human promoter sequences. We obtained a statistically significant overrepresentation of TF sites distributed with the main period of 10.1-10.5 bp in the region -50 to +120 around the transcription start site and in few locations nearby. Correlation of the positioning of the TF sites with the nucleosomes is further reinforced by sequence-directed mapping of the nucleosomes, a method previously developed.

An analysis of a class of DNA sequence reading molecules

Linked polyamides are a class of designed molecules that bind in the minor groove of double-stranded DNA in a partially sequence-specific manner but have limited sequence discriminatory abilities. This suggests a need for design alternatives to create molecules with enhanced sequence specificity. In this report we present formal proofs of the theoretical limits of the DNA sequence specificity of hypothetical sequence reading molecules as a function of their base recognition properties and sequence content and length of their target sequence. We prove that molecules containing nonspecific readers at critical positions within the molecule may have enhanced sequence specificity over molecules composed entirely of base specific reading elements. We also determine optimal patterns of base recognition for molecules in order to optimize their target sequence specificity. We also examine the effect of the length of a polyamide (i.e., the number of base pairs it binds) on its sequence discriminatory ability and determine necessary concentration dependent constraints on the binding free energies in order for longer polyamides to have greater sequence specificity than shorter ones. We show that unless the discriminatory ability of a ring for its preferred base is very strong, longer polyamides do not necessarily have greater sequence specificity over shorter ones when compared at the same molar concentration.

Chromosome translocation based on illegitimate recombination in human tumors

Recurrent chromosome translocations in nonhematological tumors are restricted to specific subtypes, and their mechanism is currently unknown. Analysis of the sequence data of 113 interchromosomal junctions derived from 77 Ewing's tumors carrying the characteristic t(11;22) translocation indicate that, in this tumor, translocations are initiated independently on each chromosome in regions that lack site specific recombination signal. Local sequence duplications, deletions, and, most importantly, inversions that are diagnostic of DNA hairpin formation indicate that, at the breakpoint, single-stranded DNA ends are processed individually before interchromosomal joining. Taken together, these observations suggest that chromosome translocations in Ewing's tumors are mediated through a genuine illegitimate recombination mechanism.

Searching for regulatory elements in human noncoding sequences

Important progress has been made in the past two years in the identification of Pol II promoters. For most other regulatory elements, however, current biological knowledge is still insufficient to allow the development of prediction tools. The phylogenetic-footprinting strategy, which is based on the comparative analysis of homologous sequences, is a very efficient approach to identify new unknown regulatory elements. The recent organization of large-scale sequencing projects for some model vertebrate organisms will be extremely valuable for the prediction of regulatory elements in the human genome.