Theoretical calculations of base-base interactions in nucleic acids: II. Stacking interactions in polynucleotides

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.

The binding affinity of Ff gene 5 protein depends on the nearest-neighbor composition of the ssDNA substrate

The Ff gene 5 protein (g5p) is considered to be a nonspecific single-stranded DNA binding protein, because it binds cooperatively to and saturates the Ff bacteriophage single-stranded DNA genome and other single-stranded polynucleotides. However, the binding affinity Komega (the intrinsic binding constant times a cooperativity factor) differs by over an order of magnitude for binding to single-stranded polynucleotides such as poly[d(A)] and poly[d(C)]. A polynucleotide that is more stacked, like poly[d(A)], binds more weakly than one that is less stacked, like poly[d(C)]. To test the hypothesis that DNA base stacking, a nearest-neighbor property, is involved in the binding affinity of the Ff g5p for different DNA sequences, Komega values were determined as a function of NaCl concentration for binding to six synthetic sequences 48 nucleotides in length: dA48, dC48, d(AAC)16, d(ACC)16, d(AACC)12, and d(AAACC)9A3. The binding affinities of the protein for these sequences were indeed found to be related to the nearest-neighbor compositions of the sequences, rather than to simple base compositions. That is, the g5p binding site, which is spanned by four nucleotides, discriminates among these sequences on the basis of the relative numbers of nearest neighbors (AA, CC, and AC plus CA) in the sequence. The results support the hypothesis that the extent of base stacking/unstacking of the free, nonbound ssDNA plays an important role in the binding affinity of the Ff gene 5 protein.

Energetics of left and right handed models of DNA

It has been shown by model building studies that various right handed and left handed models are compatible with X-ray data of B-DNA and C-DNA. These models are also found to be in good agreement with infrared dichroism data. Detailed potential energy calculations have now been carried out for these models, viz., right and left handed B-DNA and right and left handed C-DNA. It is found that base sugar stacking and interactions involving the phosphate groups are the dominant forces for stabilizing a particular structure. For some sequences, viz., A-A, T-A and C-A, left handed stacking is quite favourable in both B and C structures. But intranucleotide interactions make the left B-DNA unfavourable while the left C-DNA structure is more stable, for all the sequences, than the right C-DNA structure, proposed from fibre data. For the hexanucleoside pentaphosphate fragments the same trend is observed, with the right handed B-DNA being the most stable of the four models studied. However, the left C-DNA structure is only marginally higher in energy, particularly if the shielding effect of the counter ions, on the phosphate group is taken into consideration.

Conformational flexibility of DNA: polymorphism and handedness

It is shown that left-handed duplexes are possible for A, B, and D forms of DNA. These duplexes are stereochemically satisfactory and are consistent with the observed x-ray intensity data. On scrutiny the refined right-handed models of B and D DNA by Arnott and coworkers are found to be stereochemically unacceptable. It was possible to formulate a stereochemical guideline for molecular model building based on theory and analysis of single-crystal structure data of dinucleoside monophosphate and higher oligomers. This led to both right- and left-handed DNA duplexes. The right-handed B and D DNA duplexes so obtained are stereochemically superior to earlier models and agree well with the observed x-ray intensity data. The observation that DNA can exist in either handedness for all the polymorphous forms of DNA at once explained A in equilibrium B and B in equilibrium D transitions. Hence it is confirmed that polymorphism of DNA is a reflection on the conformational flexibility inherent in DNA, the same cause that ultimately allows DNA in either handedness. The possibility of various types of right- and left-handed duplexes generated by using dinucleoside monophosphate and trinucleoside diphosphate as repeating units resulted in a variety of models, called RL models. All these models have alternating right and left helical segments and inverted stacking at the bend region as suggested by us earlier. It turns out that the B-Z DNA model of Wang et al. is only an example of RL models.

Conformational transitions in closed circular DNA molecules. II. Biological implications

A model of regulation of gene action based on the theory of conformational transitions in closed circular DNA molecules is proposed and discussed in connection with the mechanisms of cellular differentiation. The model predicts two main types of regulation of gene action (1) the change in the topological linking number of the DNA loops leading to the change in the amount of the DNA segments present in the transcriptionally active A-form and (2) the change of some nucleotide sequences in closed superhelical DNA loops resulting in conformational transitions of some of the other sequences in the same loop. The first type of regulation may explain the mechanism of terminal differentiation of the stem cells and the changes accompanying the malignant transformation. The second one may explain the variegated position effect of the gene and determination of the stem cells during ontogenesis.

Unlinked strands as a topological constraint on chromosomal DNA, plasmid integration, and DNA repair

It is proposed that circular chromosomal DNA must be constrained such that the two strands are topologically unlinked. This structure can be replicated without strand breakage (nicking) by locally acting enzymes. The recently discovered form V DNA satisfies this topological constraint. The structure is likely to consist of left-handed and right-handed segments of double helix. The constraint of zero linkage has to be preserved by DNA repair, plasmid insertion and by crossing over. The argument presented in this paper is a topological one, following W.F. Pohl, that linkage is a global property that cannot be measured by locally acting enzymes. In contrast to Pohl no argument in favor of a side-by-side structure is presented. A zero linkage constraint would be hereditary and compatible with a multitude of local structures.