Configurational statistics of polysaccharide chains. Part II. Cellulose

The characteristic ratio of unperturbed cellulose chain was computed as a function of the angle τ at, the bridge oxygen atom and the degree of polymerization. Very high values of the order of 40 or more, depending on the angle at the bridge oxygen atom, have been obtained for this ratio, indicating that cellulose chains are highly extended. The average dimensions of cellulose chains are found to be sensitive even for small changes in the angle at the bridge oxygen, and these chains attain the character of a random coil in very high molecular weight range (degree of polymerization greater than 2000). The large differences in the unperturbed dimensions of cellulosic chains observed in different solvents have been attributed to the possible small changes in the angle τ caused by specific solvent, interactions.     

Configurational statistics of 1,6-linked glycans

We have calculated the unperturbed dimensions of some 1,6-linked glycans by regarding the homopolymer as a ‘copolymer’ in which the different ‘comonomers’ are conformers characterized by the dihedral angle about the C5C6 bond.We find that the characteristic ratio is small (e.g. for the α-glucan it is about 3·3 to 3·4). From the results of some model calculations in which the dihedral angle about the C5C6 bond is fixed, we argue that these low values arise from bonding geometry effects, which are at least as important as the additional conformational freedom from rotation about the C5C6 bond.     

Base sequence effects in double helical DNA. I. Potential energy estimates of local base morphology

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.     

Hydration of the phosphate group in double-helical DNA.

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.     

Effects of long-range correlations in DNA on sequence alignment score statistics.

Long-range correlations in genomic base composition are a ubiquitous statistical feature among many eukaryotic genomes. In this article, these correlations are shown to substantially influence the statistics of sequence alignment scores. Using a Gaussian approximation to model the correlated score landscape, we calculate the corrections to the scale parameter lambda of the extreme value distribution of alignment scores. Our approximate analytic results are supported by a detailed numerical study based on a simple algorithm to efficiently generate long-range correlated random sequences. We find both, mean and exponential tail of the score distribution for long-range correlated sequences to be substantially shifted compared to random sequences with independent nucleotides. The significance of measured alignment scores will therefore change upon incorporation of the correlations in the null model. We discuss the magnitude of this effect in a biological context.     

Complete DNA sequence of yeast chromosome II.

In the framework of the EU genome-sequencing programmes, the complete DNA sequence of the yeast Saccharomyces cerevisiae chromosome II (807 188 bp) has been determined. At present, this is the largest eukaryotic chromosome entirely sequenced. A total of 410 open reading frames (ORFs) were identified, covering 72% of the sequence. Similarity searches revealed that 124 ORFs (30%) correspond to genes of known function, 51 ORFs (12.5%) appear to be homologues of genes whose functions are known, 52 others (12.5%) have homologues the functions of which are not well defined and another 33 of the novel putative genes (8%) exhibit a degree of similarity which is insufficient to confidently assign function. Of the genes on chromosome II, 37-45% are thus of unpredicted function. Among the novel putative genes, we found several that are related to genes that perform differentiated functions in multicellular organisms of are involved in malignancy. In addition to a compact arrangement of potential protein coding sequences, the analysis of this chromosome confirmed general chromosome patterns but also revealed particular novel features of chromosomal organization. Alternating regional variations in average base composition correlate with variations in local gene density along chromosome II, as observed in chromosomes XI and III. We propose that functional ARS elements are preferably located in the AT-rich regions that have a spacing of approximately 110 kb. Similarly, the 13 tRNA genes and the three Ty elements of chromosome II are found in AT-rich regions. In chromosome II, the distribution of coding sequences between the two strands is biased, with a ratio of 1.3:1. An interesting aspect regarding the evolution of the eukaryotic genome is the finding that chromosome II has a high degree of internal genetic redundancy, amounting to 16% of the coding capacity.     

Hydrodynamic properties of a double-helical model for DNA.

The translational and rotational diffusion coefficients of very short DNA fragments have been calculated using a double-helical bead model in which each nucleotide is represented by one bead. The radius of the helix is regarded as an adjustable parameter. The translational coefficient and the perpendicular rotation coefficient agree very well with experimental values for oligonuclotides with 8, 12, and 20 base pairs, for a single value of the helical radius of about 10 A. We have also calculated a nuclear magnetic resonance relaxation time in which the coefficient for rotation about the main axis is involved. As found previously with cylindrical models, the results deviate from experimental values, indicating that the internal motion of the bases has a remarkable amplitude. An attempt to quantify the extent of internal motions is presented.     

A structural similarity analysis of double-helical DNA

A database of the structural properties of all 32,896 unique DNA octamer sequences has been calculated, including information on stability, the minimum energy conformation and flexibility. The contents of the database have been analysed using a variety of Euclidean distance similarity measures. A global comparison of sequence similarity with structural similarity shows that the structural properties of DNA are much less diverse than the sequences, and that DNA sequence space is larger and more diverse than DNA structure space. Thus, there are many very different sequences that have very similar structural properties, and this may be useful for identifying DNA motifs that have similar functional properties that are not apparent from the sequences. On the other hand, there are also small numbers of almost identical sequences that have very different structural properties, and these could give rise to false-positives in methods used to identify function based on sequence alignment. A simple validation test demonstrates that structural similarity can differentiate between promoter and non-promoter DNA. Combining structural and sequence similarity improves promoter recall beyond that possible using either similarity measure alone, demonstrating that there is indeed information available in the structure of double-helical DNA that is not readily apparent from the sequence.     

Nucleotide sequence analysis of DNA. II. Complete nucleotide sequence of the cohesive ends of bacteriophage lambda DNA

At the ends of bacteriophage λ DNA, the 5′-terminated strands are 12 nucleotides longer than the 3′-terminated strands. The complete sequence of deoxynucleotides in both the protruding 5′-terminated single strands of λ DNA has been determined by partial repair and by complete repair followed by sequencing of isolated oligonucleotides. Starting from the 5′-end of the left-hand cohesive end, the 12 nucleotides are in the sequence dpGpGpGpCpGpGpCpGpApCpCpT. The sequence from the right-hand cohesive end is exactly complementary to that from the left-hand end.     

Hexaamminecobalt(III) binding environments on double-helical DNA.

Previous cation nmr evidence suggests that univalent cations such as Na+ bind to DNA in a diffuse, nonspecific manner, whereas di- and trivalent cations show distinct binding heterogeneity. Here are reported 59Co- and 23Na-nmr measurements of the %GC dependence of the DNA binding behavior of the trivalent hexaamminecobalt(III) cation. When Co(NH3)6Cl3 titrations are performed on one mammalian and three bacterial DNAs, evidence is found for at least three distinct classes of bound Co(NH3)6(3+). A comparison of titration curves for all four DNAs demonstrates that an increase in GC content correlates with an increase in the fraction of specific Co(NH3)6(3+). binding sites. For M. lysodeikticus DNA (72% GC), a slowly exchanging class of bound 59Co(NH3)6(3+) is apparent. This class of sites is saturated at very low binding densities (between 0.02 and 0.03 cobalt cations per DNA phosphate). At higher binding densities (greater than 0.03), the signal due to slowly exchanging 59Co(NH3)6(3+) disappears into the noise, and a single 59Co(NH3)6(3+) signal is observed. Within the sensitivity limitations of these measurements, no evidence for slowly exchanging bound 59Co(NH3)6(3+) could be found for any of the other DNAs, for which a single, rapidly exchanging 59Co(NH3)6(3+) signal is observed at all binding densities. For this rapidly exchanging signal, for all four DNAs, the measured 59Co(NH3)6(3+) nmr parameters depend significantly on (a) binding density and (b) GC content of the DNA.(ABSTRACT TRUNCATED AT 250 WORDS)     

Base pair stacking in nucleosome DNA and bendability sequence pattern

DNA deformation in the nucleosome involves partial unstacking between bases and base pairs. By adjusting orientations of different base-pair stacks relative to the histone octamer surface, the optimal set of stacks and their positions is derived, resulting in a sequence pattern, theoretically best suitable for nucleosome DNA. The sequence is very much consistent with available experimental data, thus, suggesting a common eukaryotic nucleosome DNA bendability sequence pattern based exclusively on the very basics of DNA.     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be −0.5° per base pair (bp). This is equivalent to a change of −0.15±0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55±0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38±0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained.     

Effect of 1-methyladenine on double-helical DNA structures

Methylation at the N1 site of adenine leads to the formation of cytotoxic 1-methyladenine (m1A). Since the N1 site of adenine is involved in the hydrogen bonding of T·A and A·T Watson–Crick base pairs, it is expected that the pairing interactions will be disrupted upon 1-methylation. In this study, high-resolution NMR investigations were performed to determine the effect of m1A on double-helical DNA structures. Interestingly, instead of disrupting hydrogen bonding, we found that 1-methylation altered the T·A Watson–Crick base pair to T(anti)·m1A(syn) Hoogsteen base pair, providing insights into the observed differences in AlkB-repair efficiency between dsDNA and ssDNA.     

Superhelicity of nucleosomal DNA changes its double-helical repeat

Winding DNA in a superhelix can be considered a process consisting of two smooth deformations: bending and twisting. The extra twist angle introduced by winding DNA into the nucleosomal superhelix is calculated by means of the Crick formula to be -0.5 degrees per base pair (bp). This is equivalent to a change of -0.15 +/- 0.015 bp in the DNA double-helical repeat. Free DNA in solution is known to have a helical repeat of 10.55 +/- 0.1 bp. On the other hand, a weighted average of various estimates of the DNA repeat in the nucleosome is 10.38 +/- 0.02. The difference happens to be perfectly accounted for by the superhelicity of the nucleosomal DNA. This implies that the latter is essentially nonconstrained .     

On the hydration of DNA. II. Base composition dependence of the net hydration of DNA

The dependence of the net hydration of DNA on its base composition has been measured by density gradient ultracentrifugation of three DNA's in a series of cesium and lithium salt solutions of different water activities. Extrapolation to zero water activity showed the dependence of the partial specific volume on base composition to be very small for CsDNA and aero for LiDNA. At least 99% of the dependence of buoyant density on base composition can be accounted for on the basis of a differential hydration, with a mole of adenine–thymine pairs binding about 2 moles more water than a mole of guanine–cytosine pairs in CsCl.