Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

Thermodynamic characterization of deoxyribooligonucleotide duplexes containing bulges

Ultraviolet absorption techniques were used to study the thermodynamics of duplex formation for a DNA decamer, d(GCGAAAAGCG).d(CGCTTTTCGC), and a series of related duplexes, each of which contains a bulged base centered in the A.T tract. Thermodynamic parameters were obtained from nonlinear least-squares fits of the melting curves and the concentration dependences of the melting temperatures. Duplexes containing a localized single-base bulge were found to be 3.5-4.6 kcal/mol less stable than the decamer at 37 degrees C. These results indicate that both the identity of the bulged base and the strand in which it is located may influence the amount by which the duplex is destabilized. Bulged bases located in the T-strand, d(CGCTTYTTCGC), in position Y, were observed to be slightly more destabilizing than those located in the A-strand, d(GCGAAXAAGCG), in position X. Bulged purines may be more destabilizing than bulged pyrimidines.     

On the use of phasing experiments to measure helical repeat and bulge loop-associated twist in RNA.

In a phasing experiment, two bends are introduced into a long duplex RNA or DNA and the number of base pairs between them varied. When electrophoresed in a gel, the set of molecules may show a periodic variation in mobility that contains information about the twist associated with the bends and the intervening helix. We show how a set of three phasing experiments can be used to extract this information, and apply it to an RNA helix bend at the bulge sequence A2. The bulge introduces a negative (left-handed) twist of approximately 30 degrees; at low temperatures, it is mostly confined to the 5' side of the bulge. The apparent helical repeat of random sequence RNA measured in these experiments was 10.2 +/- 0.1 base pairs, an unexpectedly low value. It is likely that moderate curvative of the RNA helix axis (30-40 degrees over 80 bp) has affected the measurement.     

Scanning tunneling microscopy of mercapto-hexyl-oligonucleotides attached to gold.

6-mercapto hexyl-oligonucleotides bind to a gold surface strongly enough to permit imaging by a scanning tunneling microscope (STM). STM images showed worm-like chains that were approximately 12-(A-wide for single-stranded DNA and 20-(A-wide for double-stranded DNA. The chain lengths corresponded to 3.4 +/- 0.4 A per basepair for double-stranded DNA and 2.2 +/- 0.4 A per base for single-stranded DNA. This unexpectedly short length for single-stranded DNA was confirmed using oligomers with both single- and double-stranded regions. When the attachment of the samples was weakened (by imaging in water or scraping with the STM tip) the images changed to pairs of "blobs," apparently reflecting the attachment points of the molecules to the gold surface. Given this interpretation, images of DNA containing a five-base bulge imply that the bulge bends the oligomer by 90 degrees +/- 20 degrees.     

Solution structure of a trinucleotide A-T-A bulge loop within a DNA duplex.

We have synthesized an oligodeoxynucleotide duplex, d(G-C-A-T-C-G-A-T-A-G-C-T-A-C-G).d(C-G-T-A-G-C-C-G-A-T-C-G), with a three-base bulge loop (A-T-A) at a central site in the first strand. Nuclear Overhauser experiments (NOESY) in H2O indicate that the GC base pairs flanking the bulge loop are intact between 0 and 25 degrees C. Nuclear Overhauser effects in both H2O and D2O indicate that all bases within the bulge loop are stacked into the helix. These unpaired bases retain an anti conformation about their glycosidic bonds as they stack within the duplex. The absence of normal sequential connectivities between the two cytosine residues flanking the bulge site on the opposite strand indicates a disruption in the geometry of this base step upon insertion of the bulged bases into the helix. This conformational perturbation is more akin to a shearing apart of the bases, which laterally separates the two halves of the molecule, rather than the "wedge" model often invoked for single-base bulges. Using molecular dynamics calculations, with both NOE-derived proton-proton distances and relaxation matrix-calculated NOESY cross peak volumes as restraints, we have determined the solution structure of an A-T-A bulge loop within a DNA duplex. The bulged bases are stacked among themselves and with the guanine bases on either side of the loop. All three of the bulged bases are displaced by 2-3 A into the major groove, increasing the solvent accessibility of these residues. The ATA-bulge duplex is significantly kinked at the site of the lesion, in agreement with previously reported electron microscopy and gel retardation studies on bulge-containing duplexes [Hsieh, C.-H., & Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A 86, 4833-4837; Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840]. Bending occurs in a direction away from the bulge-containing strand, and we find a significant twist difference of 84 degrees between the two base pairs flanking the bulge loop site. This value represents 58% of the twist difference for base pairs four steps apart in B-DNA. These results suggest a structural mechanism for the bending of DNA induced by unpaired bases, as well as accounting for the effect bulge loops may have on the secondary and tertiary structures of nucleic acids.     

Positional effect of single bulge nucleotide on PNA(peptide nucleic acid)/DNA hybrid stability.

We report positional effect of bulge nucleotide on PNA/DNA hybrid stability. CD spectra showed that PNA/DNA hybrids required at least seven base pairings at a stem region to form a bulged structure. On the other hand, DNA/DNA could form bulged structure when there are only four base pairings adjacent to the bulge nucleotide. We discuss why PNA requests such a many base pairings to form bulged structure from a nearest neighbor standpoint.     

Visualization of diagnostic heteroduplex DNAs from cystic fibrosis deletion heterozygotes provides an estimate of the kinking of DNA by bulged bases.

Previous studies (Hsieh, C.-H., and Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837) of DNAs containing extra or bulged bases on one of the two strands of a duplex showed that they exhibit slower electrophoretic mobility than non-bulged DNAs, indicating that bulges create stiff kinks in the DNA. Here we paired a 97-base single-stranded DNA from the intact cystic fibrosis (CF) gene with a complementary 94-base strand containing a central 3-base deletion (delta F508), common to many CF patients. This produced a 94-base pair DNA with a central 3-base bulge. Visualization of these DNAs by electron microscopy showed that twice as many bulge-containing DNAs had a central kink as compared with the non-bulged controls. Examination of the distribution of kinking angles showed that the bulged population contained 5-7-fold more molecules with a central kink of 80 +/- 10 degrees than did the control molecules. When the 3-base bulge was replaced by a 3-base gap, the resulting duplex DNA showed central kinks with a somewhat lower frequency but greater range of kinking angles.     

Bulge loops used to measure the helical twist of RNA in solution.

Bulge loops are commonly found in helical segments of cellular RNAs. When incorporated into long double-stranded RNAs, they may introduce points of flexibility or permanent bend that can be detected by the altered electrophoretic gel mobility of the RNA. We find that a single An or Un bulge loop near the middle of a long RNA helix significantly retards the RNA during polyacrylamide gel electrophoresis if n greater than or equal to 2. The mobility of an RNA containing two A2 bulges various periodically with the number of base pairs between the bulges. We interpret this to mean that A2 bulges varies periodically with the number of base pairs between the bulges. We interpret this to mean that Z2 bulges form torsionally stiff bends in the helix; the gel mobility reaches a minimum when the total helical twist between the bulges rotates the arms of the molecule into a cis conformation. The gel mobilities are proportional to the predicted end-to-end distance of the RNA if the average RNA helical repeat is 11.8 +/- 0.2 bp/turn and there is no helical twist (3 +/- 9 degrees) associated with the bulge (data obtained in 0.15 M Na+). Other sizes and sequences of bulges have very different effects on RNA helix conformation and flexibility. U2 bulges bend the helix to a much smaller degree than A2 bulges, while longer A or U bulge sequences probably allow bends of 90 degrees or more; all of these may be fairly flexible joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Stereochemical control of small molecule binding to bulged DNA: comparison of structures of spirocyclic enantiomer-bulged DNA complexes

The solution structure of the complex formed between an oligonucleotide containing a two-base bulge (5'-CACGCAGTTCGGAC.5'-GTCCGATGCGTG) and ent-DDI, a designed synthetic agent, has been elucidated using high-resolution NMR spectroscopy and restrained molecular dynamic simulation. Ent-DDI is a left-handed wedge-shaped spirocyclic molecule whose aglycone portion is an enantiomer of DDI, which mimics the spirocyclic geometry of the natural product, NCSi-gb, formed by base-catalyzed activation of the enediyne antibiotic neocarzinostatin. The benzindanone moiety of ent-DDI intercalates between the A6.T21 and the T9.A20 base pairs, overlapping with portions of the purine bases; the dihydronaphthalenone moiety is positioned in the minor groove along the G7-T8-T9 bulge sequence; and the aminoglycoside is in the middle of the minor groove, approaching A20 of the nonbulged strand. This alignment of ent-DDI along the DNA helical duplex is in the reverse direction to that of DDI. The aminoglycoside moiety of ent-DDI is positioned in the 3' direction from the bulge region, whereas that of the DDI is positioned in the 5' direction from the same site. This reverse binding orientation within the bulge site is the natural consequence of the opposite handedness imposed by the spirocyclic ring junction and permits the aromatic ring systems of the two spirocyclic enantiomers access to the bulge region. NMR and CD data indicate that the DNA in the DDI-bulged DNA complex undergoes a larger conformational change upon complex formation in comparison to the ent-DDI-bulged DNA, explaining the different binding affinities of the two drugs to the bulged DNA. In addition, there are different placements of the bulge bases in the helical duplex in the two complexes. One bulge base (G7) stacks inside the helix, and the other one (T8) is extrahelical in the DDI-bulged DNA complex, whereas both bulge bases in the ent-DDI-bulged DNA complex prefer extrahelical positions for drug binding. Elucidation of the detailed binding characteristics of the synthetic spirocyclic enantiomers provides a rational basis for the design of stereochemically controlled drugs for bulge binding sites.     

NMR evidence of the stabilisation by the carcinogen N-2-acetylaminofluorene of a frameshift mutagenesis intermediate.

Two heteroduplexes d(C1A2C3T4C5G6C7A8C9A10C11)-d (G12T13G14T15G16G17A18G19T20G21) containing a bulged guanine either unmodified or modified with the carcinogen N-2-acetylaminofluorene (AAF) have been studied by nuclear magnetic resonance (NMR) as models of slipped mutagenic intermediates (SMI). Conformational equilibria are observed in both the unmodified and the AAF-modified heteroduplexes. The major conformation of the unmodified duplex is one where the extra guanine is stacked in the helix and the major conformation of the AAF-modified heteroduplex is one where the AAF is external to the helix. Unusual sugar proton chemical shifts of C5- and G6-AAF indicate that the AAF ring is pointing out in the 5' direction. A strong increase in the modified heteroduplex melting temperature (+15 degrees C) is observed. Moreover, in contrast to the unmodified heteroduplex, which shows extensive melting in the vicinity of the bulged guanine, the base pairs around the bulge in the AAF-modified heteroduplex remain paired at temperatures up to 30 degrees C. This exceptional stability of the site around the bulged modified guanine is suggested to be responsible for the high rate of -1 mutation induced by AAF at repetitive sequences.     

DNA bending induced by high mobility group proteins studied by fluorescence resonance energy transfer.

The HMG domains of the chromosomal high mobility group proteins homologous to the vertebrate HMG1 and HMG2 proteins preferentially recognize distorted DNA structures. DNA binding also induces a substantial bend. Using fluorescence resonance energy transfer (FRET), we have determined the changes in the end-to-end distance consequent on the binding of selected insect counterparts of HMG1 to two DNA fragments, one of 18 bp containing a single dA(2) bulge and a second of 27 bp with two dA(2) bulges. The observed changes are consistent with overall bend angles for the complex of the single HMG domain with one bulge and of two domains with two bulges of approximately 90-100 degrees and approximately 180-200 degrees, respectively. The former value contrasts with an inferred value of 150 degrees reported by Heyduk et al. (1) for the bend induced by a single domain. We also observe that the induced bend angle is unaffected by the presence of the C-terminal acidic region. The DNA bend of approximately 95 degrees observed in the HMG domain complexes is similar in magnitude to that induced by the TATA-binding protein (80 degrees), each monomeric unit of the integration host factor (80 degrees), and the LEF-1 HMG domain (107 degrees). We suggest this value may represent a steric limitation on the extent of DNA bending induced by a single DNA-binding motif.     

Extent of helix perturbation associated with DNA modification by the o-acetyl derivative of the carcinogen 4-hydroxyaminoquinoline-1-oxide

Duplex unwinding associated with DNA modification by 4-acetoxyaminoquinoline-1-oxide, a model ultimate carcinogen of 4-nitroquinoline-1-oxide, has been determined by the agarose gel electrophoresis band-shift method. An average unwinding angle per stable adduct of -15.1 degrees +/- 1.5 degrees for negatively supercoiled topoisomers and -6.5 degrees +/- 1.4 degrees for positively supercoiled topoisomers was obtained. Because of the different proportion of stable adducts (dGuo-N2-AQO, dGuo-C8-AQO, dAdo-N6-AQO) between negatively (8:1.5:0.5) and positively (5:2.5:1) supercoiled topoisomers, the difference in unwinding angles is suggestive of a diverse contribution of the various adducts to the overall conformational change. Since the largest unwinding angle was coupled with the highest proportion of dGuo-N2-AQO adduct, it is likely that this adduct is the most distortive lesion. A contribution of sites of base loss to DNA unwinding was also observed.     

Comparative NMR study of A(n)-bulge loops in DNA duplexes: intrahelical stacking of A, A-A, and A-A-A bulge loops.

We have prepared a series of deoxyoligonucleotide duplexes of the sequence d(G-C-A-T-C-G-X-G-C-T-A-C-G).d(C-G-T-A-G-C-C-G-A-T-G-C), in which X represents either one (A), two (A-A), or three (A-A-A) unpaired adenine basis. Using two-dimensional proton and phosphorus NMR spectroscopy, we have characterized conformational features of these bulge-loop duplexes in solution. We find that Watson-Crick hydrogen bonding is intact for all 12 base pairs, including the GC bases that flank the bulge loop. Observation of NOE connectivities in both H2O and D2O allows us to unambiguously localize all of the bulged adenine residues to intrahelical positions within the duplex. This is in contrast to an earlier model for multiple-base bulge loops in DNA [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840], in which all but the most 5' bulged base are looped out into solution. We find that insertion of two or three bases into the duplex results in the disruption of specific sequential NOEs for the base step across from the bulge loop site on the opposite strand. This disruption is characterized by a partial shearing apart of these bases, such that certain sequential NOEs for this base step are preserved. We observe a downfield-shifted phosphorus resonance, which we assign in the A-A-A bulge duplex to the 3' side of the last bulged adenine residue. Proton and phosphorus chemical shift trends within the An-bulge duplex series indicate that there is an additive effect on the structural perturbations caused by additional unpaired bases within the bulge loop. This finding parallels previous observations [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840; Hsieh, C.-H., & Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837] on the magnitude of the induced bending of DNA duplexes by multiple-base bulge loops.     

Solution structure of a wedge-shaped synthetic molecule at a two-base bulge site in DNA

The solution structure of the complex formed between an oligonucleotide containing a two-base bulge (5'-CACGCAGTTCGGAC.5'-GTCCGATGCGTG) and DDI, a designed synthetic agent, has been elucidated using high-resolution NMR spectroscopy and restrained molecular dynamic simulation. DDI, which has been found to modulate DNA strand slippage synthesis by DNA polymerase I [Kappen, L. S., Xi, Z., Jones, G. B., and Goldberg, I. H. (2003) Biochemistry 42, 2166-2173], is a wedge-shaped spirocyclic molecule whose aglycone structure closely resembles that of the natural product, NCSi-gb, which strongly binds to an oligonucleotide containing a two-base bulge. Changes in chemical shifts of the DNA upon complex formation and intermolecular NOEs between DDI and the bulged DNA duplex indicate that agent specifically binds to the bulge site of DNA. The benzindanone moiety of DDI intercalates via the minor groove into the G7-T8-T9.A20 pocket, which consists of a helical base pair and two unpaired bulge bases, stacking with the G7 and A20 bases. On the other hand, the dihydronaphthalenone and aminoglycoside moieties are positioned in the minor groove. The aminoglycoside, which is attached to spirocyclic ring, aligns along the A20T21G22 sequence of the nonbulged strand, while the dihydronaphthalenone, which is restrained by the spirocyclic structure, is positioned near the G7-T8-T9 bulge site. The aminoglycoside is closely aligned with the dihydronaphthalenone, preventing its intercalation into the bulge site. In the complex, the unpaired purine (G7) is intrahelical and stacks with the intercalating moiety of DDI, whereas the unpaired pyrimidine (T8) is extrahelical. The structure of the complex formed by binding of the synthetic agent to the two-base bulged DNA reveals a binding mode that differs in important details from that of the natural product, explaining the different binding specificity for the bulge sites of DNA. The structure of the DDI-bulged DNA complex provides insight into the structure-binding affinity relationship, providing a rational basis for the design of specific, high-affinity probes of the role of bulged nucleic acid structures in various biological processes.     

Effects of bulge composition and flanking sequence on the kinking of DNA by bulged bases

We recently showed that bulged bases kink duplex DNA, with the degree of kinking increasing in roughly equal increments as the number of bases in the bulge increases from one to four [Hsieh, C.-H., & Griffith, J.D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837]. Here we have examined the kinking of DNA by single A, C, G, or T bulges with different neighboring base pairs. Synthetic 30 base pair (bp) duplex DNAs containing 2 single-base bulges spaced by 10 bp were ligated head to tail, and their electrophoretic behavior in highly cross-linked gels was examined. All bulge-containing DNAs showed marked electrophoretic retardations as compared to non-bulge-containing DNA. Regardless of the sequence of the flanking base pairs, purine bulges produced greater retardations than pyrimidine bulges. Furthermore, C and T bulges produced the same retardations as did G and A bulges. Bulged DNA containing different flanking base pairs showed marked differences in electrophoretic mobility. For C-bulged DNA, the greatest retardations were observed with G.C neighbors, the least with T.A neighbors, and an intermediate amount with a mixture of neighboring base pairs. For A-bulged DNA, the retardations were greatest with G.C neighbors, less with T.A neighbors, even less with a mixture of neighboring base pairs, and finally least with C.G neighbors. Thus flanking base pairs affect C-bulged DNA and A-bulged DNA differently, and G.C and C.G flanking base pairs were seen to have very different effects. These results imply an important role of base stacking in determining how neighboring base pairs influence the kinking of DNA by a single-base bulge.     

Thermodynamic parameters based on a nearest-neighbor model for DNA sequences with a single-bulge loop

All 64 possible thermodynamic parameters for a single-bulge loop in the middle of a sequence were derived from optical melting studies. The relative stability of a single bulge depended on both the type of bulged base and its flanking base pairs. The contribution of the single bulge to helix stability ranged from 3.69 kcal/mol for a TAT bulge to -1.05 kcal/mol for an ACC bulge. Thermodynamics for 10 sequences with a GTG bulge were determined to test the applicability of the nearest-neighbor model to a single-bulge loop. Thermodynamic parameters for the GTG bulge and Watson-Crick base pairs predict, DeltaH degrees, DeltaS degrees, and T(M)(50 microM) values with average deviations of 3.0%, 4.3%, 4.7%, and 0.9 degrees C, respectively. The prediction accuracy was within the limits of what can be expected for a nearest-neighbor model. This certified that the thermodynamics for single-bulge loops can be estimated adequately using a nearest-neighbor model.     

Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.