Curved DNA without AA/TT dinucleotide step.

The evidence is accumulating that dinucleotide steps other than AA/TT affect DNA flexure of AnTm (m + n greater than = 4) containing fragments. However, it is not clear whether macroscopic DNA flexure without AA/TT steps might occur. In this paper we demonstrate the anomaly in electrophoretic mobility of non AA/TT repetitive DNA sequences which is a function of sequence phasing. Therefore, our results show that PyPu (TA) and AG/CT steps, angulary separated by close to 180 degrees from Pu/Py (GC) and GG/CC steps, bend DNA, even in the absence of AnTm tracts.     

MOSAIC: segmenting multiple aligned DNA sequences

MOSAIC is a set of tools for the segmentation of multiple aligned DNA sequences into homogeneous zones. The segmentation is based on the distribution of mutational events along the alignment. As an example, the analysis of one repeated sequence belonging to the subtelomeric regions of the yeast genome is presented. AVAILABILITY: Free access from ftp://ftp.biomath.jussieu.fr/pub/papers/MOSAIC     

The relationship between periodic dinucleotides and the nucleosomal DNA deformation revealed by normal mode analysis

Nucleosomes, which contain DNA and proteins, are the basic unit of eukaryotic chromatins. Polymers such as DNA and proteins are dynamic, and their conformational changes can lead to functional changes. Periodic dinucleotide patterns exist in nucleosomal DNA chains and play an important role in the nucleosome structure. In this paper, we use normal mode analysis to detect significant structural deformations of nucleosomal DNA and investigate the relationship between periodic dinucleotides and DNA motions. We have found that periodic dinucleotides are usually located at the peaks or valleys of DNA and protein motions, revealing that they dominate the nucleosome dynamics. Also, a specific dinucleotide pattern CA/TG appears most frequently.     

Correlation between the flexibility and periodic dinucleotide patterns in yeast nucleosomal DNA sequences

Nucleosome formation and positioning, which play important roles in a number of biological processes, are thought to be related to the distinctive periodic dinucleotide patterns observed in the DNA sequence wrapped around the protein octamer. Previous research shows that flexibility is a key structural property of a nucleosomal DNA sequence. However, the relationship between the flexibility and the periodic dinucleotide patterns has received little attention in research in the past. In this study, we propose the use of three different models to measure the flexibility of yeast DNA sequences. Although the three models involve different parameters, they deliver consistent results showing that yeast nucleosomal DNA sequences are more flexible than non-nucleosomal ones. In contrast to random flexibility values along non-nucleosomal DNA sequences, the flexibility of nucleosomal DNA sequences shows a clear periodicity of 10.14 base pairs, which is consistent with the periodicity of dinucleotide distributions. We also demonstrate that there is a strong relationship between the peak positions of the flexibility and the dinucleotide frequencies. Correlation between the flexibility and the dinucleotide patterns of CA/TG, CG, GC, GG/CC, AG/CT, AC/GT and GA/TC are positive with an average value of 0.5946. The highest correlation is shown by CA/TG with a value of 0.7438 and the lowest correlation is shown by AA/TT with a value of −0.7424. The source codes and data sets are available for downloading on http://www.hy8.com/bioinformatics.htm.     

Stability of DNA duplexes containing GG, CC, AA, and TT mismatches

We employed salt-dependent differential scanning calorimetric measurements to characterize the stability of six oligomeric DNA duplexes (5'-GCCGGAXTGCCGG-3'/5'-CCGGCAYTCCGGC-3') that contain in the central XY position the GC, AT, GG, CC, AA, or TT base pair. The heat-induced helix-to-coil transitions of all the duplexes are associated with positive changes in heat capacity, DeltaC(p), ranging from 0.43 to 0.53 kcal/mol. Positive values of DeltaC(p) result in strong temperature dependences of changes in enthalpy, DeltaH degrees, and entropy, DeltaS degrees , accompanying duplex melting and cause melting free energies, DeltaG degrees, to exhibit characteristically curved shapes. These observations suggest that DeltaC(p) needs to be carefully taken into account when the parameters of duplex stability are extrapolated to temperatures distant from the transition temperature, T(M). Comparison of the calorimetric and van't Hoff enthalpies revealed that none of the duplexes studied in this work exhibits two-state melting. Within the context of the central AXT/TYA triplet, the thermal and thermodynamic stabilities of the duplexes in question change in the following order: GC > AT > GG > AA approximately TT > CC. Our estimates revealed that the thermodynamic impact of the GG, AA, and TT mismatches is confined within the central triplet. In contrast, the thermodynamic impact of the CC mismatch propagates into the adjacent helix domains and may involve 7-9 bp. We discuss implications of our results for understanding the origins of initial recognition of mismatched DNA sites by enzymes of the DNA repair machinery.     

Methylation of nucleosomal and nuclease sensitive DNA.

The proportion of cytosines methylated in the DNA of nucleosome oligomers and of core particles appears indistinguishable from that of total nuclear DNA from CHO cells. However the DNA in nucleoprotein which is initially released from nuclei by treatment with very low levels of micrococcal nuclease and the first 10% of material rendered acid soluble by treatment of nuclei with DNase I are enriched 2 fold in their content of 5 methylcytosine. (Cessation of hydrolysis by nuclease occurs concomitantly with precipitation of nucleosomal core particles).     

Flexibility of nucleosomal DNA

In this study transient electric birefringence (TEB) has been used to investigate the molecular flexibility of short fragments of DNA. Nucleosomal DNA always exhibits negative birefringence and Kerr behavior was observed up to high field strengths (6 KV/cm). The value of the Kerr constant is 3.5 10^?2 e.s.u.. Birefringence decays were single exponentials and a field dependence of the molecular orientational relaxation time τ was found: it is explained by an inherent flexibility of the DNA molecule. A 20 % decrease in the calculated length was observed with fields applied as low as 2 KV/cm. The results obtained at very low fields establish TEB as a method well suited to calculate accurate values for the length of small fragments of DNA: the τ value of 4.3 μsec corresponds to a DNA length of 660 Å.     

Anisotropic flexibility of DNA and the nucleosomal structure.

Potential energy calculations of the DNA duplex dimeric subunit show that the double helix may be bent in the direction of minor and major grooves much more easily than in other directions. It is found that the total winding angle of DNA decreases upon such bending. A new model for DNA folding in the nucleosome is proposed on the basis of these findings according to which the DNA molecule is kinked each fifth base pair to the side of the minor and major grooves alternatively. The model explains the known contradiction between a C-like circular dichroism for the nucleosomal DNA and the nuclease digestion data, which testify to the B-form of DNA.     

Effect of flanking residues on CA and AA dinucleotides: some rationale

Effect of flanking base pairs on CA and AA dinluceotide step-geometry has been studied using molecular dynamics method. Sixteen dodecameric sequences are constructed for each doublet with all possible bases at their 5' and 3' positions along with their complementary sequences. Structural parameter, formation of Extra Watson Crick bifurcated hydrogen bonds along or across the strands and effect of sodium ions are studied for these sequences. It is found that geometry of CA doublet step is perturbed by the neighboring base pairs, which might be due to inherent flexibility of the step. Flexible character of CA step is reflected in its low bifurcated hydrogen bond formation capability and lower preference of sodium ions to enter in minor or major grooves. AA step on the other hand is quite rigid according to different structural parameters and respond much less to environmental changes due to formation of strong Extra Watson-Crick hydrogen bonds.     

Integrating ambiguously aligned regions of DNA sequences in phylogenetic analyses without violating positional homology

Phylogenetic analyses of non-protein-coding nucleotide sequences such as ribosomal RNA genes, internal transcribed spacers, and introns are often impeded by regions of the alignments that are ambiguously aligned. These regions are characterized by the presence of gaps and their uncertain positions, no matter which optimization criteria are used. This problem is particularly acute in large-scale phylogenetic studies and when aligning highly diverged sequences. Accommodating these regions, where positional homology is likely to be violated, in phylogenetic analyses has been dealt with very differently by molecular systematists and evolutionists, ranging from the total exclusion of these regions to the inclusion of every position regardless of ambiguity in the alignment. We present a new method that allows the inclusion of ambiguously aligned regions without violating homology. In this three-step procedure, first homologous regions of the alignment containing ambiguously aligned sequences are delimited. Second, each ambiguously aligned region is unequivocally coded as a new character, replacing its respective ambiguous region. Third, each of the coded characters is subjected to a specific step matrix to account for the differential number of changes (summing substitutions and indels) needed to transform one sequence to another. The optimal number of steps included in the step matrix is the one derived from the pairwise alignment with the greatest similarity and the least number of steps. In addition to potentially enhancing phylogenetic resolution and support, by integrating previously nonaccessible characters without violating positional homology, this new approach can improve branch length estimations when using parsimony.     

Preferred crossover sites on polyomavirus DNA.

RmI is a hybrid replicon consisting of polyomavirus (Py) and mouse sequences that yields unit-length polyomavirus DNA via recombination between two directly repeated viral sequences of 182 base pairs (S repeats). To define the contribution of the S repeats in this intramolecular recombination, we derived from RmI a series of replicons containing the original S repeats as well as additional direct viral repeats which were 1 to 2 kilobases in length (L repeats). After mouse 3T6 cells were transfected with these constructs, recombination products that displayed the physical properties of homologous recombinants were detected. The structures of these recombinants indicated that whereas repeat length influences the likelihood of recombination, crossover occurs preferentially near the S repeats, provided that one of them is proximal to the viral origin of replication. This finding suggests that recombination near the S repeats depends on a process initiated near the viral origin of replication.     

Rapid spontaneous accessibility of nucleosomal DNA

DNA wrapped in nucleosomes is sterically occluded, creating obstacles for proteins that must bind it. How proteins gain access to DNA buried inside nucleosomes is not known. Here we report measurements of the rates of spontaneous nucleosome conformational changes in which a stretch of DNA transiently unwraps off the histone surface, starting from one end of the nucleosome, and then rewraps. The rates are rapid. Nucleosomal DNA remains fully wrapped for only approximately 250 ms before spontaneously unwrapping; unwrapped DNA rewraps within approximately 10-50 ms. Spontaneous unwrapping of nucleosomal DNA allows any protein rapid access even to buried stretches of the DNA. Our results explain how remodeling factors can be recruited to particular nucleosomes on a biologically relevant timescale, and they imply that the major impediment to entry of RNA polymerase into a nucleosome is rewrapping of nucleosomal DNA, not unwrapping.     

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.

     

Geometry of the nucleosomal DNA superhelix

Nucleosome stability is largely an indirect measure of DNA sequence based on the material properties of DNA and the ability of a sequence to assume the required left-handed superhelical conformation. Here we focus attention only on the geometry of the superhelix and present two distinct mathematical expressions that rely on the DNA helical parameters (Shift, Slide, Rise, Tilt, Roll, Twist). One representation requires torsion for superhelix formation; the other requires shear. To compare these mathematical expressions to experimental data we develop a strategy for Fourier-filtering the helical parameters that identifies necessary and sufficient conditions to achieve a high-resolution model of the nucleosome superhelix. We apply this filtering strategy to 24 high-resolution structures of the nucleosome and demonstrate that all structures have a highly conserved distribution of Roll, Slide and Twist that involves two length scales. One length scale spans the entire length of nucleosomal DNA. The other is associated with the helix repeat. Our strategy also enables us to identify ground state or simple nucleosomes and altered nucleosome structures. These results form a basis for characterizing structural variations in the emerging family of nucleosome structures and a method for further developing structure-based models of nucleosome stability.