On the stability of nucleic acid structures in solution: enthalpy-entropy compensations, internal rotations and reversibility.

The thermodynamics of self-association (stacking) of free bases and nucleotides, intramolecular stacking in dinucleotides, nearest-neighbour base pair stacking interactions in duplex DNA and RNA, and the formation of hairpin loops illustrate enthalpy/entropy compensations. Large stacking exothermicities are associated with large negative entropy changes that ensure that delta G is small, permitting readily reversible associations in solution. We rationalise enthalpy/entropy compensations with reference to residual motions and torsional vibrations which make a larger entropic contribution to binding when - delta H approximately kT (thermal energy at room temperature), than when - delta H >> kT. We present a factorisation of experimental free energies for helix formation in terms of approximate contributions from the restriction of rotations, hydrophobic interactions, electrostatic interactions due to base stacking, and contributions from hydrogen bonding, and estimate the adverse free energy cost per rotor (mainly entropy) of ordering the phosphate backbone as between 1.9 and 5.4 kJ mol-1 [averaged over 12 rotors per base pair for A-U on A-U stacking (lower limit), and G-C on C-G stacking (upper limit)]. The largest cost is associated with the most exothermic stacking interactions, while the range of values is consistent with earlier conclusions from data on the fusion of hydrocarbon chains (lower value), and with entropy changes in covalent isomerisations of small molecules involving severe restrictions (upper value).     

Single strand targeted triplex formation: targeting purine-pyrimidine mixed sequences using abasic linkers.

Foldback triplex-forming oligonucleotides (FTFOs) that contain an abasic linker, [2-(4-aminobutyr-1-yl)-1,3-propanediol] (APD linker), in the Hoogsteen domain against pyrimidine bases of a C:G and a T:A base pair were studied for their relative stability and sequence specificity of triplex formation. In general, the APD linker has less destabilizing effect against a C:G base pair than a T:A base pair. Incorporation of three APD linker moieties resulted in decreased binding to the target, which was comparable to results observed with three imperfectly matched natural base triplets. The APD linker incorporation did not result in the loss of sequence specificity of FTFOs, unlike in the case of normal triplex-forming oligonucleotides (TFOs). The introduction of a positively charged abasic linker, however, resulted in decreased stability of the triplex, because of loss of hydrogen bonding and stacking interactions in the major groove. The results of a molecular modeling study show that APD linker can be readily incorporated without any change in the conformation of the natural sugar-phosphate backbone conserving overall triple helix geometry. Further, the modeling study suggests a hydrogen bond formation between the amino group of linker and N4 of cytosine mediated by a solvent molecule (water) in the floor of the base triplet in addition to a contribution from the positive charge on the APD linker amino group. Either a direct or water-mediated hydrogen bond between the amino group of the APD linker and the O4 of thymine is unlikely when the linker is placed against a T:A base pair.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

Vibrational fluctuations of hydrogen bonds in a DNA double helix with nonuniform base pairs.

The Green's function technique is applied to a study of breathing modes in a DNA double helix which contains a region of different base pairs from the rest of the double helix. The calculation is performed on a G-C helix in the B conformation with four consecutive base pairs replaced by A-T. The average stretch in hydrogen bonds is found amplified around the A-T base pair region compared with that of poly(dG)-poly(dC). This is likely related to the A-T regions lower stability against hydrogen bond melting. The A-T region may be considered to be the initiation site for melting in such a helix.     

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.     

Structural features and hydration of a dodecamer duplex containing two C.A mispairs.

X-ray diffraction techniques have been used to characterise the crystal and molecular structure of the deoxyoligomer d(C-G-C-A-A-A-T-T-C-G-C-G) at 2.5A resolution. The final R factor is 0.19 with the location of 78 solvent molecules. The oligomer crystallises in a B-DNA type conformation with two strands coiled about each other to produce a duplex. This double helix consists of four A.T and six G.C Watson-Crick base pairs and two C.A mispairs. The mismatched base pairs adopt a "wobble" type structure with the cytosine displaced laterally into the major groove, the adenine into the minor groove. We have proposed that the two close contacts observed in the C.A pairing represent two hydrogen bonds one of which results from protonation of adenine. The mispairs are accommodated in the double helix with small adjustments in the conformation of the sugar-phosphate backbone. Details of the backbone conformation, base stacking interactions, thermal parameters and the hydration are now presented and compared with those of the native oligomer d(C-G-C-G-A-A-T-T-C-G-C-G) and with variations of this sequence containing G.T and G.A mispairs.     

Stability of XGCGCp, GCGCYp, and XGCGCYp helixes: an empirical estimate of the energetics of hydrogen bonds in nucleic acids

The stabilizing effects of dangling ends and terminal base pairs on the core helix GCGC are reported. Enthalpy and entropy changes of helix formation were measured spectrophotometrically for AGCGCU, UGCGCA, GGCGCCp, CGCGCGp, and the corresponding pentamers XGCGCp and GCGCYp containing the GCGC core plus a dangling end. Each 5' dangling end increases helix stability at 37 degrees C roughly 0.2 kcal/mol and each 3' end from 0.8 to 1.7 kcal/mol. The free energy increments for dangling ends on GCGC are similar to the corresponding increments reported for the GGCC core [Freier, S. M., Alkema, D., Sinclair, A., Neilson, T., & Turner, D. H. (1985) Biochemistry 24, 4533-4539], indicating a nearest-neighbor model is adequate for prediction of stabilization due to dangling ends. Nearest-neighbor parameters for prediction of the free energy effects of adding dangling ends and terminal base pairs next to G.C pairs are presented. Comparison of these free energy changes is used to partition the free energy of base pair formation into contributions of "stacking" and "pairing". If pairing contributions are due to hydrogen bonding, the results suggest stacking and hydrogen bonding make roughly comparable favorable contributions to the stability of a terminal base pair. The free energy increment associated with forming a hydrogen bond is estimated to be -1 kcal/mol of hydrogen bond.     

Sequence dependence for the energetics of dangling ends and terminal base pairs in ribonucleic acid

Stability increments of terminal unpaired nucleotides (dangling ends) and terminal base pairs on the core helixes AUGCAU and UGCGCA are reported. Enthalpy, entropy, and free energy changes of helix formation were measured spectrophotometrically for 18 oligoribonucleotides containing the core sequences. The results indicate 3' dangling purines add more stability than 3' dangling pyrimidines. In most cases, the additional stability from a 3' dangling end on an AU base pair is less than that on a GC base pair [Freier, S.M., Burger, B.J., Alkema, D., Neilson, T., & Turner, D.H. (1985) Biochemistry 22, 6198-6206]. The sequence dependence provides a test for the importance of dangling ends for various RNA interactions. Correlations are suggested with codon context effects and with the three-dimensional structure of yeast phenylalanine transfer RNA. In the latter case, all terminal unpaired nucleotides having stability increments more favorable than -1 kcal/mol are stacked on the adjacent base pair. All terminal unpaired nucleotides having stability increments less favorable than -0.3 kcal/mol are not stacked on the adjacent base pair. In several cases, this lack of stacking is associated with a turn in the sugar-phosphate backbone. This suggests stability increments measured on oligoribonucleotides may be useful for predicting tertiary structure in large RNA molecules. Comparison of the stability increments for terminal dangling ends and base pairs, and of terminal GC and AU base pairs, indicates the free energy increment associated with forming a hydrogen bond can be about -1 kcal/mol of hydrogen bond.     

Interaction between the Z-type DNA duplex and 1,3-propanediamine: crystal structure of d(CACGTG)2 at 1.2 A resolution

The crystal structure of a hexamer duplex d(CACGTG)(2) has been determined and refined to an R-factor of 18.3% using X-ray data up to 1.2 A resolution. The sequence crystallizes as a left-handed Z-form double helix with Watson-Crick base pairing. There is one hexamer duplex, a spermine molecule, 71 water molecules, and an unexpected diamine (Z-5, 1,3-propanediamine, C(3)H(10)N(2)) in the asymmetric unit. This is the high-resolution non-disordered structure of a Z-DNA hexamer containing two AT base pairs in the interior of a duplex with no modifications such as bromination or methylation on cytosine bases. This structure does not possess multivalent cations such as cobalt hexaammine that are known to stabilize Z-DNA. The overall duplex structure and its crystal interactions are similar to those of the pure-spermine form of the d(CGCGCG)(2) structure. The spine of hydration in the minor groove is intact except in the vicinity of the T5A8 base pair. The binding of the Z-5 molecule in the minor grove of the d(CACGTG)(2) duplex appears to have a profound effect in conferring stability to a Z-DNA conformation via electrostatic complementarity and hydrogen bonding interactions. The successive base stacking geometry in d(CACGTG)(2) is similar to the corresponding steps in d(CG)(3). These results suggest that specific polyamines such as Z-5 could serve as powerful inducers of Z-type conformation in unmodified DNA sequences with AT base pairs. This structure provides a molecular basis for stabilizing AT base pairs incorporated into an alternating d(CG) sequence.     

Degradation of nucleic acid in aqueous solution by ionizing radiation: III. The correlation of radiation damage with change in melting transition - model experiments

The melting behavior of polydeoxynucleotide double helices of known structure is analyzed in terms of the thermodynamics of helix stability, taking into account separately those contributions to the transition free energy that are proportional to the numbers of polymer molecules and those that are proportional to the numbers of base pairs formed. From the analysis of the melting transitions of helices having an alternating (d-) A-T, G-C base-pair sequence and containing either single-strand nicks or both nicks and damaged thymine bases, the effects of these structural lesions are assessed; it is concluded that, in a moderately long helix of this sequence (400 base pairs), the initial introduction of one mid-chain doublestrand break or single-strand break produces respectively some 3.5 or 4 times as much depression in the transition temperature (T_m) as does the destruction of a single internal A-T base pair.     

Contributions of dangling end stacking and terminal base-pair formation to the stabilities of XGGCCp, XCCGGp, XGGCCYp, and XCCGGYp helixes

The role of stacking in terminal base-pair formation was studied by comparison of the stability increments for dangling ends to those for fully formed base pairs. Thermodynamic parameters were measured spectrophotometrically for helix formation of the hexanucleotides AGGCCUp, UGGCCAp, CGGCCGp, GCCGGCp, and UCCGGAp and for the corresponding pentanucleotides containing a 5'-dangling end on the GGCCp or CCGGp core helix. In 1 M NaCl at 1 X 10(-4) M strands, a 5'-dangling nucleotide in this series increases the duplex melting temperature (Tm) only 0-4 degrees C, about the same as adding a 5'-phosphate. In contrast, a 3'-dangling nucleotide increases the Tm at 1 X 10(-4) M strands 7-23 degrees C, depending on the sequence [Freier, S. M., Burger, B. J., Alkema, D., Neilson, T., & Turner, D. H. (1983) Biochemistry 22, 6198-6206]. These results are consistent with stacking patterns observed in A-form RNA. The stability increments from terminal A.U, C.G, or U.A base pairs on GGCC or a terminal U.A pair on CCGG are nearly equal to the sums of the stability increments from the corresponding dangling ends. This suggests stacking plays a large role in nucleic acid stability. The stability increment from the terminal base pairs in GCCGGCp, however, is about 5 times the sum of the corresponding dangling ends, suggesting hydrogen bonding can also make important contributions.     

G.T wobble base-pairing in Z-DNA at 1.0 A atomic resolution: the crystal structure of d(CGCGTG).

The DNA oligomer d(CGCGTG) crystallizes as a Z-DNA double helix containing two guanine-thymine base pair mismatches of the wobble type. The crystal diffracts to 1 A resolution and the structure has been solved and refined. At this resolution, a large amount of information is revealed about the organization of the water molecules in the lattice generally and more specifically around the wobble base pairs. By comparing this structure with the analogous high resolution structure of d(CGCGCG) we can visualize the structural changes as well as the reorganization of the solvent molecules associated with wobble base pairing. There is only a small distortion of the Z-DNA backbone resulting from introduction of the GT mismatched base pairs. The water molecules cluster around the wobble base pair taking up all of the hydrogen bonding capabilities of the bases due to wobble pairing. These bridging water molecules serve to stabilize the base-base interaction and, thus, may be generally important for base mispairing either in DNA or in RNA molecules.     

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.     

Electrochemical detection of lesions in DNA

Electrochemical DNA-based sensors that exploit the inherent sensitivity of DNA-mediated charge transport (CT) to base pair stacking perturbations are capable of detecting base pair mismatches and some common base damage products. Here, using DNA-modified gold electrodes, monitoring the electrocatalytic reduction of DNA-bound methylene blue, we examine a wide range of base analogues and DNA damage products. Among those detected are base damage products O4-methyl-thymine, O6-methyl-guanine, 8-oxo-guanine, and 5-hydroxy-cytosine, as well as a therapeutic base, nebularine. The efficiency of DNA-mediated CT is found not to depend on the thermodynamic stability of the helix. However, general trends in how base modifications affect CT efficiency are apparent. Modifications to the hydrogen bonding interface in Watson-Crick base pairs yields a substantial loss in CT efficiency, as does added steric bulk. Base structure modifications that may induce base conformational changes also appear to attenuate CT in DNA as do those that bury hydrophilic groups within the DNA helix. Addition and subtraction of methyl groups that do not disrupt hydrogen bonding interactions do not have a large effect on CT efficiency. This sensitive detection methodology based upon DNA-mediated CT may have utility in diagnostic applications and implicates DNA-mediated CT as a possible damage detection mechanism for DNA repair enzymes.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Solution structure of duplex DNA containing the mutagenic lesion 1,N(2)-etheno-2'-deoxyguanine

We have used high-resolution NMR spectroscopy and molecular dynamics simulations to determine the solution structure of DNA containing the genotoxic lesion 1, N (2)-etheno-2'-deoxyguanosine (epsilonG), paired to dC. The NMR data suggest the presence of a major, minimally perturbed structure at neutral pH. NOESY spectra indicate the presence of a right-handed helix with all nucleotides in anti, 2'-deoxyribose conformations within the C2'-endo/C1'-exo range and proper Watson-Crick base pair alignments outside the lesion site. The epsilonG residue remains deeply embedded inside the helix and stacks between the flanking base pairs. The lesion partner dC is extrahelical and is located in the minor groove of the duplex, where it is highly exposed to solvent. Upon acidification of the sample, a second conformation at the lesion site of the duplex emerges, with protonation of the lesion partner dC and possible formation of a Hoogsteen base pair. Restrained molecular dynamics simulations of the neutral-pH structure generated a set of three-dimensional models that show epsilonG inside the helix, where the lesion is stabilized by stacking interactions with flanking bases but without participating in hydrogen bonding. The lesion counterbase dC is displaced in the minor groove of the duplex where it can form a hydrogen bond with the sugar O4' atom of a residue 2 bp away.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

The energetics of small internal loops in RNA

The energetics of small internal loops are important for prediction of RNA secondary and tertiary structure, selection of drug target sites, and understanding RNA structure-function relationships. Hydrogen bonding, base stacking, electrostatic interactions, backbone distortion, and base-pair size compatibility all contribute to the energetics of small internal loops. Thus, the sequence dependence of these energetics are idiosyncratic. Current approximations for predicting the free energies of internal loops consider size, asymmetry, closing base pairs, and the potential to form GA, GG, or UU pairs. The database of known three-dimensional structures allows for comparison with the models used for predicting stability from sequence.