Conformational and structural constraints in double-helical polynucleotides

Conformational analysis of antiparallel double-helical polynucleotides with Watson-Crick base pairing was reduced to a four-dimensional problem using original mathematical methods. In the four-dimensional conformational space the family of structures, characterized by the base-pair stacking with the most stable conformations in water solution as well as in the solid state, was localized. For the C′₂-endo sugar pucker, both right-handed and left-handed structures were found; right-handed structures only, however, seem to be allowed for the C′₃-endo pucker, the only possible one for ribonucleotides with base stacking.     

Mechanism for unwinding of double-helical polynucleotides by formaldehyde

We have formulated a mathematical model for the unwinding process of polynucleotides induced by formaldehyde by making use of the master equation approach. Our model assumes that the formaldehyde kinetics follows the helix–coil theory of Zimm and Bragg. The model incorporates experimental features such as the interplay between pH-dependent and -independent reactions and the dependence of the initial and final helix stability on the unwinding. That is, our model incorporates the existence of the critical point such that if the pH of the system is high enough (i.e., pH ? 7), the unwinding of polynucleotides occurs by means of two separate chemical reactions, either by the HCHO–imino reaction alone or by the HCHO–amino reaction alone. The critical point depends on the ionic strength, temperature, and the formaldehyde concentration satisfying the relation s = 1 + K_Uλ, where s is the helix stability parameter, K_U is the binding constant for the HCHO–imino reaction, and λ is the formaldehyde concentration. Thus above the critical point, i.e., s < (1 + K_Uλ), the unwinding is due to the reaction of imino groups with formaldehyde, and below the critical point, i.e., s > (1 + K_Uλ), the unwinding is due to the reaction of amino groups with formaldehyde. Although, below the critical point, the imino group does not participate directly in the rate-determining step, it participates indirectly in such a manner that it reduces the initial helix stability parameter s to s/(1 + K_Uλ) by forming a complex. The kinetic constants in the master equation have been determined by making use of the principle of detailed balance and the breathing mechanism proposed in the past. With this model we have quantitatively explained the anomalous temperature dependence observed in the kinetics of formaldehyde-induced poly[d(A-T)] and the salt dependence observed in poly(A·U).     

Effect of base-pair sequence on the conformations and thermally induced transitions in oligodeoxyribonucleotides containing only AT base pairs

Tm curves, CD spectra, and kinetics results of the self-complementary DNA dodecamers d(A6T6), d(A3T3A3T3), d(A2T2A2T2A2T2), d(ATATATATATAT), and d(T6A6) demonstrate that the thermal transitions of these oligomers at low salt concentration involve a hairpin intermediate. At high salt concentrations (greater than 0.1 M Na+) only a duplex to denatured-strand transition appears to occur. The temperature and salt-concentration regions of the transitions are very sequence dependent. Alternating-type AT sequences have a lower duplex stability and a greater tendency to form hairpins than sequences containing more nonalternating AT base pairs. Of the two nonalternating sequences, d(T6A6) is significantly less stable than d(A6T6). Both oligomers have CD curves that are very similar to the unusual CD spectrum of poly(dA).poly(dT). The Raman spectra of these two oligomers are also quite similar, but at low temperature, small intensity differences in two backbone modes and three nucleoside vibrations are obtained. The hairpin to duplex transition for the AT dodecamers was examined by salt-jump kinetics measurements. The transition is faster than transitions for palindromic-sequence oligomers containing terminal GC base pairs. Stopped-flow kinetics studies indicate that the transition is second order and has a relatively low activation energy. The reaction rate increases with increasing ionic strength. These results are consistent with a three-step mechanism for the hairpin to duplex reaction: (i) fraying of the hairpin oligomers' terminal base pairs, (ii) a rate-determining bimolecular step involving formation of a cruciform-type intermediate from two hairpin oligomers with open terminal base pairs, and (iii) base-pair migration and formation in the intermediate to give the duplex.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.     

Modeling of hydration of incorrect nucleic acid base pairs by the Monte Carlo method

The hydration of water bridged base pairs of nucleic acids have been simulated via the Monte Carlo method. The simulation have shown that water molecules forming H-bonds with both bases preserve this H-bonding with large probability in the water surrounding. This fact supports the supposition about the important role of water molecules in wrong base pair formation and about the role of these base pairs in the structure and functioning of nucleic acids.     

Adsorption of single-stranded and double-helical polynucleotides on the mercury electrode

The adsorption of single-stranded polynucleotides and double-helical DNA on the dropping mercury electrode has been studied with the aid of Breyer's alternating current (a.c.) polarography. Our results indicate that all three constituents of polynucleotides (residues of bases, sugar, and phosphoric acid) are involved in the adsorption. At neutral pH their participation in adsorption depends on the ionic strength, the potential of the electrode, and the conformation of the polynucleotide in the solution. At an ionic strength of about 0.1, double-helical DNA is adsorbed electrostatically on a positively charged electrode surface by inadequately masked negative charges of the phosphate groups. At a higher ionic srength (about 0.5), this electrostatic adsorption is no longer detectable by using a.c. polarography; under these conditions it is probable that native DNA is adsorbed around the potential of the electrocapillary maximum with the aid of sugar residues and a few bases. Single-stranded polynucleotides, on the other hand, are primarily adsorbed by means of the bases. Desorption of double-helical DNA occurs around a potential of ?1.2 V against SCE. At this potential, the helical regions of single-stranded polynucleotides are also desorbed. Desorption of the disordered regions of single-stranded polynucleotides occurs at more negative potentials. Adsorption and desorption of a small number of bases released from double-helical DNA was evident in the a.c. polarograms only at elevated temperature, or at room temperature after degradation of DNA by sonication.     

Thermodynamics of DNA containing very stable chemically modified base pairs

DNA chemical modifications caused by the binding of some antitumor drugs give rise to a very strong local stabilization of the double helix. These sites melt at a temperature that is well above the melting temperatures of ordinary AT and GC base pairs. In this work we have examined the melting behavior of DNA containing very stable sites. Analytical expressions were derived and used to evaluate the thermodynamic properties of homopolymer DNA with several different distributions of stable sites. The results were extended to DNA with a heterogeneous sequence of AT and GC base pairs. The results were compared to the melting properties of DNA with ordinary covalent interstrand cross-links. It was found that, as with an ordinary interstrand cross-link, a single strongly stabilized site makes a DNA's melting temperature (T(m)) independent of strand concentration. However in contrast to a DNA with an interstrand cross-link, a strongly stabilized site makes the DNA's T(m) independent of DNA length and equal to T(infinity), the melting temperature of an infinite length DNA with the same GC-content and without a stabilized site. Moreover, at a temperature where more than 80% of base pairs are melted, the number of ordinary (non-modified) helical base pairs (n) is independent of both the DNA length and the location of the stabilized sites. For this condition, n(T) = (2 omega-a)S/(1-S) and S = exp[DeltaS(T(infinity)-T)/(RT)] where omega is the number of strongly stabilized sites in the DNA chain, a is the number of DNA ends that contain a stabilized site, and DeltaS, T, and R are the base pair entropy change, the temperature, and the universal gas constant per mole. The above expression is valid for a temperature interval that corresponds to n<0.2N for omega=1, and n<0.1N for omega>1, where N is the number of ordinary base pairs in the DNA chain.     

Dissociation kinetics of 19 base paired oligonucleotide-DNA duplexes containing different single mismatched base pairs.

The dissociation kinetics of 19 base paired oligonucleotide-DNA duplex containing a various single mismatched base pair are studied on dried agarose gels. The kinetics of the dissociation are first order under our experimental conditions. The incorporation of a single mismatched base pair destabilizes the DNA duplexes to some extent, the amount depending on the nature of the mismatched base pair. G-T and G-A mismatches slightly destabilize a duplex, while A-A, T-T, C-T and C-A mismatches significantly destabilize it. The activation energy for the overall dissociation processes for these oligonucleotide-DNA duplexes containing 19 base pairs is 52 +/- 2 Kcal mol-1 as determined from the slope of Arrhenius plot.     

Telomeric DNA oligonucleotides form novel intramolecular structures containing guanine-guanine base pairs

Structural properties of DNA oligonucleotides corresponding to the single-stranded molecular terminus of telomeres from several organisms were analyzed. Based on physical studies including nondenaturing polyacrylamide gel electrophoresis, absorbance thermal denaturation analysis, and 1H and 31P nuclear magnetic resonance spectroscopy, we conclude that these molecules can self-associate by forming non-Watson-Crick, guanine.guanine based-paired, intramolecular structures. These structures form below 40 degrees C at moderate ionic strength and neutral pH and behave like hairpin duplexes in nondenaturing polyacrylamide gels. Detailed analysis of the hairpin structure formed by the telomeric sequence from Tetrahymena, (T2G4)4, shows that it is a unique structure stabilized by hydrogen bonds and contains G residues in the syn conformation. We propose that this novel form of DNA is important for telomere function and sets a precedent for the biological relevance of non-Watson-Crick base-paired DNA structures.     

Possible bioactive conformations of alpha-melanotropin

Rat brain microtubules were prepared at the adult stage and from immature (i.e., 4-day-old) animals. At an early stage of development, the composition of microtubule-associated proteins is qualitatively different from that found at the adult stage [(1982) Eur. J. Biochem. 129, 465-471]. The influence of calmodulin on the time course of assembly of second cycle microtubules was compared at both stages of brain development (i.e., microtubules originating from 4-day-old and adult animals). In the presence of Ca2+ the inhibition of microtubule assembly was more pronounced at a young stage of brain development than at the adult stage. Cross-linking studies with 125I-labeled calmodulin further established that the two major microtubule-associated proteins, MAP2 and TAU were able to bind to calmodulin at both stages of brain development but with different intensities. The labeling with 125I-labeled calmodulin was Ca2+-dependent, specific, displaced by unlabeled calmodulin and trifluoperazine.     

Parallel-stranded duplex DNA containing dA · dU base pairs

DNA oligonucleotides with dA and dU residues can form duplexes with trans d(A · U) base pairing and the sugar-phosphate backbone in a parallel-stranded orientation, as previously established for oligonucleotides with d(A · T) base pairs. The properties of such parallel-stranded DNA (ps-DNA) 25-mer duplexes have been characterized by absorption (uv), CD, ir, and fluorescence spectroscopy, as well as by nuclease sensitivity. Comparisons were made with duplex molecules containing (a) dT in both strands, (b) dU in one strand and dT in the second, and (c) the same base combinations in reference antiparallel-stranded (aps) structures. Thermodynamic analysis revealed that total replacement of deoxythymine by deoxyuridine was accompanied by destabilization of the ps-helix (reduction in T_m by −13°C in 2 mM MgGl₂, 10 mM Na-cacodylate). The U-containing ps-helix (U1 · U2) also melted 14°C lower than the corresponding aps-helix under the same ionic conditions; this difference was very close to that observed between ps and aps duplexes with d(A · T) base pairs. Force field minimized structures of the various ps and aps duplexes with either d(A · T) or d(A · U) base pairs ps/aps and dT/dU combinations are presented. The energy-minimized helical parameters did not differ significantly between the DNAs containing dT and dU. © 1996 John Wiley & Sons, Inc.     

Additional possibilities of the formation of incorrect nucleotide pairs

Calculations of intermolecular interaction energies for the systems consisting of two nitrogen bases and one or two molecules of water have been performed by the atom--atom potential function method. In some energy minima corresponding to the arrangement of the bases with one of the bases being in the syn-orientation with respect to the sugar two bases are linked through one H-bond and one or two water bridges. Such pairs provide an additional possible pathway for the appearance of errors during template biosynthesis.     

Contribution of individual structural components to stabilization of the secondary structure of double-helical polynucleotides

Contributions of individual structural components of the double-helical polynucleotide to the stabilization of its secondary structure have been studied. The energy of intramolecular interactions was calculated by the method of atom-atom potentials. Sections of the energy function were constructed according to the parameters determining mutual location of base pairs with the optimal conformation of ribose-phosphate backbone in the points close to A- and B-forms of DNA. A complicated nature of the contribution of different structural components of the polynucleotide to its stabilization was revealed by means of different parameters. Relationships between the change of conformational parameters corresponding to optimal values of energy for A- and B-families of nucleic acids were found.     

Crystal structure of d(GCGAAAGCT) containing a parallel-stranded duplex with homo base pairs and an anti-parallel duplex with Watson-Crick base pairs

A DNA fragment d(GCGAAAGCT), known to adopt a stable mini-hairpin structure in solution, has been crystallized in the space group I4₁22 with the unit-cell dimensions a = b = 53.4 Å and c = 54.0 Å, and the crystal structure has been determined at 2.5 Å resolution. The four nucleotide residues CGAA of the first half of the oligomer form a parallel duplex with another half through the homo base pairs, C₂:C₂⁺ (singly-protonated between the Watson– Crick sites), G₃:G₃ (between the minor groove sites), A₄:A₄ (between the major groove sites) and A₅:A₅ (between the Watson–Crick sites). The two strands remaining in the half of the parallel duplex are split away in different directions, and they pair in an anti-parallel B-form duplex with the second half extending from a neighboring parallel duplex, so that an infinite column is formed in a head-to-tail fashion along the c-axis. It seems that a hexa-ammine cobalt cation supports such a branched and bent conformation of the oligomer. One end of the parallel duplex is stacked on the corresponding end of the adjacent parallel duplex; between them, the guanine base of the first residue is stacked on the fourth ribose of another duplex.     

Preparation and properties of a nucleosomal particle containing 180 base pairs of DNA

If chromatin from chicken erythrocytes is enzymatically degraded in the presence of histone H5, nucleosomal DNA usually exceeds the size of 140 base pairs. Conditions are derived allowing the isolation of a 180 base pair particle which is subsequently characterized by histone binding and thermal denaturation studies. Association of H5 to such a particle is cooperative and occurs with a larger affinity than binding to specimens in the range of 140– 170 base pairs. Thermal denaduration studies show that some of the extra DNA participates in the main transition at 74°C (1 mm Na-cacodylate) indicating a tight binding to the histone core. Another part of the extra segment is loosely associated with the core but also distinct from free DNA.     

Thermodynamics and kinetics of the helix-coil transition of oligomers containing GC base pairs

Absorbance-temperature profiles have been determined for the following self-complementary oligonucleotides or equimolar paris of complementary oligonucleotides containing GC base pairs: A₂GCU₂, A₃GCU₃, A₄GCU₄, A₆CG + CGU₆, A₈CG + CGU₈, A₄G₂ + C₂U₄, A₅G₂ + C₂U₅, A₄G₃ + C₃U₄, and A₅G₃ + C₃U₅. In all cases cooperative melting transitions indicate double-helix formation. As was found previously, the stability of GC containing oligomer helices is much higher than that of AU helices of corresponding length. Moreover, helices with the same length and base composition but different sequences also have quite different stabilites. The melting curves were andlyzed using a zipper model and the thermodynamic parameters for the AU pairs determined previously. The effect of single-strand stacking was considered separately. According to this model, the formation of a GC pair from unstacked single strands is associated with an ethalpy change of ?15 kcal/mole. Due to the high degree of single-strand stacking at room temperature the enthalpy change for the formation of GC pairs from unstacked single strands is only ?5 to ?6 kcal/mole. (The corresponding parameters for AU pairs are ?10.7 kcal/mole and ?5 to ?6 kcal/mole.) The sequence dependence of helix stability seems to be primarily entropic since no differences in ΔH were seen among the sequence isomers. The kinetics of helix formation was investigated for the same molecules using the temperature jump technique. Recombination of strands is second order with rate constants in the range of 10⁵ to 10⁷M^?1 sec^?1 depending on the chain length and the nucleotide sequence. Within a series of oligomers of a given type, the rates of recombination decrease with increasing chain length. Oligomers with the sequence A_nGCU_n recombine six to eight times slower than the other oligomers of corresponding chain length. The experimental enthalpies of activation of 6 to 9 kcal/mole suggest a nucleation length of one or two GC base pairs. The helix dissociation process has rate constants between 0.5 and 500 sec^?1 and enthalpies of activation of 25 to 50 kcal/mole. An increase of chain length within a given nucleotide series leads to decreased rates of dissociation and increased enthalpies of activation. An investigation of the effect of ionic strength on A_nGCU_n helix formation showed that the rates of recombination increase considerably with increased ionic strength.     

Parallel-stranded duplex DNA containing blocks of trans purine-purine and purine-pyrimidine base pairs.

A 30 base pair parallel-stranded (ps) duplex ps-L1.L2 composed of two adjoined purine-purine and purine-pyrimidine sequence blocks has been characterized thermodynamically and spectroscopically. The 5'-terminal 15 residues in both strands ('left-half') consisted of the alternating d(GA)7G sequence that forms a ps homoduplex secondary structure stabilized by d(G.G) and d(A.A) base pairs. The 3'-terminal 15 positions of the sequence ('right-half') were combinations of A and T with complementary reverse Watson-Crick d(A.T) base pairing between the two strands. The characteristics of the full length duplex were compared to those of the constituent left and right halves in order to determine the compatibility of the two ps helical forms. The thermal denaturation curves and hyperchromicity profiles of all three duplexes determined by UV absorption spectroscopy were characteristic of ps-DNA, in accordance with previous studies. The thermodynamic properties of the 30 bp duplex corresponded within experimental error to the linear combination of the two 15-mers. Thus, the Tm and delta HvH of ps-L1.L2 in 10 mM MgCl2, derived from analyses according to a statistical mechanical formulation for the helix-coil transition, were 43 degrees C and 569 kJ mol-1, compared to 21 degrees C, 315 kJ mol-1 (ps-F5.F6) and 22 degrees C, 236 kJ mol-1 (ps-GA15). The UV absorption and CD spectra of ps-L1.L2 and the individual 15-mer ps motifs were also compared quantitatively. The sums of the two constituent native spectra (left+right halves) accurately matched that of the 30 bp duplex, with only small deviations in the 195-215 nm (CD) and 220-240 nm (absorption) regions. Based on analysis by native gel electrophoresis, the sequences studied formed duplex structures exclusively; there were no indications of higher order species. Chemical modification with diethyl pyrocarbonate showed no hyperreactivity of the junctional bases, indicating a smooth transition between the two parallel-stranded conformations. We conclude that under given salt conditions, oligonucleotides with normal primary chemical structures can readily form a parallel-stranded double helix based on blocks of very disparate non-canonical purine-purine and purine-pyrimidine base pairs and without perceptible destabilization at the junction. There are biological implications of these findings in relation to genetic structure and expression.     

Crystal structure of a Z-DNA fragment containing thymine/2-aminoadenine base pairs

The Z-DNA structure has been shown to form in two crystals made from self-complementary DNA hexamers d(CGTDCG) and d(CDCGTG) which contain thymine/2-aminoadenine (TD) base pairs. The latter structure has been solved and refined to 1.3 A resolution and it shows only small conformational changes due to the introduction of the TD base pairs in comparison with the structure of d(CG)3. Spectroscopic studies with these compounds demonstrate that DNA molecules containing 2-aminoadenine residues form Z-DNA slightly more easily than do those containing adenine nucleotides, but not as readily as the parent sequence containing only guanine-cytosine base pairs.