The effect of base sequence on the stability of RNA and DNA single base bulges

Forty-eight RNA duplexes were constructed that contained all common single base bulges at six different locations. The stabilities of the RNAs were determined by temperature gradient gel electrophoresis (TGGE). The relative stability of a single base bulge was dependent on both base identity and the nearest neighbor context. The single base bulges were placed into two categories. A bulged base with no identical neighboring base was defined as a Group I base bulge. Group II-bulged bases had at least one neighboring base identical to it. Group II bulges were generally more stable than Group I bulges in the same nearest neighbor environments. This indicates that position degeneracy of an unpaired base enhances stability. Differences in the mobility transition temperatures between the RNA fragments with bulges and the completely base-paired reference RNAs were related to free energy differences. Simple models for estimating the free energy contribution of single base bulges were evaluated from the free energy difference data. The contribution of a Group I bulge 5'-(XNZ)-3'.5'-(Z'-X')-3' where N is the unpaired base and X.X' and Z.Z' the neighboring base pairs, could be well-represented (+/-0.34 kcal/mol) by the equation, DeltaG((X)(N)()(Z))(.)((Z)(')(-)(X)(')()) = 3.11 + 0. 40DeltaG(s)()((XZ))(.)((Z)(')(X)(')()). DeltaG(s)()((XZ))(. )((Z)(')(X)(')()) is the stacking energy of the closing base pair doublet. By adding a constant term, delta = -0.3 kcal/mol, to the right side of the above equation, free energies of Group II bulges could also be predicted with the same accuracy. The term delta represents the stabilizing effect due to position degeneracy. A similar equation/model was applied to previous data from 32 DNA fragments with single base bulges. It predicted the free energy differences with a similar standard deviation.     

Effect of environment, conformation, sequence and base substituents on the imino proton exchange rates in guanine and inosine-containing DNA, RNA, and DNA-RNA duplexes

High-resolution 1H nuclear magnetic resonance in H2O has been used to study the effect of sequence, conformation, environmental factors and base substituents on the exchange behavior of the hydrogen-bonded imino protons of guainine X cytosine and inosine X cytosine base-pairs in DNA, RNA, and DNA-RNA duplexes. The exchange rates were determined by measurement of the spin-lattice relaxation rates of the imino protons as a function of temperature. The exchange was not altered by the presence of high concentrations of salt, and the inability of phosphate to catalyze the exchange indicates that the exchange is limited by formation of a solvent-accessible "open" state. The exchange behavior depends on the duplex conformation and sequence. Exchange from the Z form polymers was orders of magnitude slower than the corresponding duplexes in the B conformation, and the A form RNA duplexes exchanged more slowly than the B form DNA polymers with the same sequence. The exchange behavior of the DNA-RNA hybrids was dependent on whether the purine or the pyrimidine strand contained the deoxyribose sugar. For both the guanine and inosine-containing duplexes, the homopolymer duplexes exchange more slowly than the more stable alternating copolymers. For the alternating duplexes, substitution of cytosine with 5-bromo- or 5-methylcytosine slowed the exchange and increased the activation energy for exchange. The inosine-containing duplexes exchanged more rapidly than the guanosine-containing duplexes, but both showed similar changes in exchange behavior in response to changes in sequence and base substituents. The activation energies for base-pair opening in B form DNA are correlated with the van der Waals contribution to the base-base interaction energy, suggesting that the purine base is partially unstacked in the open state. Using the relaxation measurements to set an upper limit on the exchange rate in poly(dG-dC) and the tritium exchange behavior at low temperature, we find that even though Z-DNA exchanges very slowly, the activation energy is similar to that observed in the A and B form duplexes, suggesting that exchange occurs from a similar open state.     

Differential hydration of dA.dT base pairing and dA and dT bulges in deoxyoligonucleotides

The role of water in the formation of stable duplexes of nucleic acids is being studied by determining the concurrent volume change, heats, and counterion uptake that accompany the duplexation process. The variability of the volume contraction that we have observed in the formation of a variety of homoduplexes suggests that sequence and conformation acutely affect the degree of hydration. We have used a combination of densimetric and calorimetric techniques to measure the change in volume and enthalpy resulting from the mixing of two complementary strands to form (a) fully paired duplexes with 10 or 11 base pairs and (b) bulged decameric duplexes with an extra dA or dT unmatched residue. We also monitored absorbance vs temperature profiles as a function of strand and salt concentration for all four duplexes. Relative to the decamer duplex, insertion of an extra dA.dT base pair to form an undecamer duplex results in a favorable enthalpy of -5.6 kcal/mol that is nearly compensated by an unfavorable entropy term of -5.1 kcal/mol. This enthalpy difference correlates with a differential uptake of water molecules, corresponding to an additional hydration of 16 mol of water molecules/mol of base pair. Relative to the fully paired duplexes, both bulged duplexes are 12-16 degrees C less stable and exhibit marginally larger counterion uptake on forming the duplex. The enthalpy change is slightly lower for the T-bulge duplex and less still for the A-bulge duplex. The volume change results indicate that an unmatched residue increases the amount of coulombic and/or structural hydration. The combined results strongly suggest that the destabilizing forces in bulged duplexes are partially compensated by an increase in hydration levels.     

2-Aminopurine fluorescence studies of base stacking interactions at abasic sites in DNA: metal-ion and base sequence effects.

Metal-ion and sequence dependent changes in the stacking interactions of bases surrounding abasic (AB) sites in 10 different DNA duplexes were examined by incorporating the fluorescent nucleotide probe 2-aminopurine (2-AP), opposite to the site (AB-APopp) or adjacent to the site (AB-APadj) on either strand. A detailed study of the fluorescence emission and excitation spectra of these AB duplexes and their corresponding parent duplexes indicates that AB-APoppis significantly less stacked than 2-AP in the corresponding normal duplex. In general, AB-APadjon the AB strand is stacked, but AB-APadjon the opposite strand shows destabilized stacking interactions. The results also indicate that divalent cation binding to the AB duplexes contributes to destabilizaton of the base stacking interactions of AB-APopp, but has little or no effect on the stacking interactions of AB-APadj. Consistent with these results, the fluorescence of AB-APoppis 18-30-fold more sensitive to an externally added quenching agent than the parent normal duplex. When uracil DNA glycosylase binds to AB-APoppin the presence of 2.5 mM MgCl2, a 3-fold decrease in fluorescence is observed ( K d = 400 +/- 90 nM) indicating that the unstacked 2-APoppbecomes more stacked upon binding. On the basis of these fluorescence studies a model for the local base stacking interactions at these AB sites is proposed.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Thermodynamic parameters for an expanded nearest-neighbor model for the formation of RNA duplexes with single nucleotide bulges

Thirty-four RNA duplexes containing single nucleotide bulges were optically melted, and the thermodynamic parameters deltaH degrees, deltaS degrees, deltaG degrees (37), and T(M) for each sequence were determined. Data from this study were combined with data from previous thermodynamic data [Longfellow, C. E., Kierzek, R., and Turner, D. H. (1990) Biochemistry 29, 278-85] to develop a model that will more accurately predict the free energy of an RNA duplex containing a single nucleotide bulge. Differences between purine and pyrimidine bulges as well as differences between Group I duplexes, those in which the bulge is not identical to either neighboring nucleotide, and Group II duplexes, those in which the bulge is identical to at least one neighboring nucleotide, were considered. The length of the duplex, non-nearest-neighbor effects, and bulge location were also examined. A model was developed which divides sequences into two groups: those with pyrimidine bulges and those with purine bulges. The proposed model for pyrimidine bulges predicts deltaG degrees (37,bulge) = 3.9 kcal/mol + 0.10deltaG degrees (37,nn) + beta, while the model for purine bulges predicts deltaG degrees (37,bulge) = 3.3 kcal/mol - 0.30deltaG degrees (37,nn) + beta, where beta has a value of 0.0 and -0.8 kcal/mol for Group I and Group II sequences, respectively, and deltaG degrees (37,nn) is the nearest-neighbor free energy of the base pairs surrounding the bulge. The conformation of bulge loops present in rRNA was examined. Three distinct families of structures were identified. The bulge loop was either extrahelical, intercalated, or in a "side-step" conformation.     

Comparison of the base pairing properties of a series of nitroazole nucleobase analogs in the oligodeoxyribonucleotide sequence 5'-d(CGCXAATTYGCG)-3'.

The nucleoside analogs 1-(2'-deoxy-beta-D-ribofuranosyl)- 3-nitropyrrole (9), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitropyrazole (10), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitroimidazole (11) and 1-(2'-deoxy-beta-D-ribofuranosyl)-5-nitroindole (21) were incorporated into the oligonucleotide 5'-d(CGCXAATTYGCG)-3'in the fourth position from the 5'-end. Procedures for synthesis of two of the nitroazole nucleosides, 10 and 11, were developed for this study. Each of the nitroazoles was converted into a 3'-phosphoramidite for oligonucleotide synthesis by conventional automated protocols. Four oligonucleotides were synthesized for each modified nucleoside in order to obtain duplexes in which each of the four natural bases was placed opposite (position 9) the nitroazole. In order to assess the role of the nitro group on base stacking interaction, sequences were also synthesized in which the fourth base was 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole. Corresponding sequences containing an abasic site, as well as sequences containing inosine, were synthesized for comparison. Thermal melting studies yielded T m values and thermodynamic parameters. Each nucleoside analog displayed a unique pattern of base pairing preferences. The least discriminating analog was 3-nitropyrrole, for which T m values differed by 5 degrees C and Delta G 25 degrees C ranged from -6.1 to -6.5 kcal/mol. 5-Nitroindole gave duplexes with significantly higher thermal stability, with Tm values varying from 35.0 to 46.5 degrees C and -Delta G 25 degrees C ranging from 7.7 to 8.5 kcal/mol. Deoxyinosine (22), a natural analog which has found extensive use as a universal nucleoside, is far less non-discriminating than any of the nitroazole derivatives. Tm values ranged from 35.4 degrees C when paired with G to 62.3 degrees C when paired with C. The significance of the nitro substituent was determined by comparison of the base pairing properties of a simple azole nucleoside, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole (12). The pyrazole-containing sequences melt at 10-20 degrees C lower than the corresponding nitropyrazole-containing sequences. On average, the pyrazole-containing sequences were equivalent in stability (average Delta G = -4.8 kcal/mol) to the sequences containing an abasic site (average Delta G = -4.7 kcal/mol).     

Influence of an exocyclic guanine adduct on the thermal stability, conformation, and melting thermodynamics of a DNA duplex.

As part of an overall program to characterize the impact of mutagenic lesions on the physiochemical properties of DNA, we report here the results of a comparative spectroscopic study on pairs of DNA duplexes both with and without an exocyclic guanine lesion. Specifically, we have studied a family of four 13-mer duplexes of the form d(CGCATGYGTACGC).d(GCGTACZCATGCG) in which Y is either the normal deoxyguanosine residue (G) or the exocyclic guanine adduct 1,N2-propanodeoxyguanosine (X), while Z is either deoxycytosine (C) or deoxyadenosine (A). Thus, the four duplexes studied, which can be designated by the identity of their central Y.Z base pair, are a Watson-Crick duplex (GC), a duplex with a central mismatch (GA), and two duplexes with exocyclic guanine lesions (X), that differ only by the base opposite the lesion (XC and XA). The data derived from our spectroscopic measurements on these four duplexes have allowed us to evaluate the influence of the exocyclic guanine lesion, as well as the base opposite the lesion, on the conformation, thermal stability, and melting energetics of the host DNA duplex. To be specific, our circular dichroism (CD) spectra show that the exocyclic guanine lesion induces alterations in the duplex structure, while our temperature-dependent optical measurements reveal that these lesion-induced structural alterations reduce the thermal stability, the transition enthalpy, and the transition free energy of the duplex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Allosteric conversion of Z DNA to an intercalated right-handed conformation by daunomycin

Absorbance and fluorescence methods were used to measure the binding of the anticancer drug daunomycin to poly (dGdC) under ionic conditions that initially favor the left-handed Z conformation of the polymer. Drug binding was cooperative under these conditions and may be fully accounted for by an allosteric model in which the drug binds preferentially (but not exclusively) to the right-handed B conformation and shifts the polymer from the Z to an intercalated right-handed conformation. Quantitative analysis of binding isotherms in terms of the allosteric model allowed for estimation of the equilibrium constants for the conversion of a base pair at a B-Z interface from the Z to the B conformation and for the formation of a base pair in the B conformation within a stretch of helix in the Z conformation. The free energy of the Z to B conversion of a base pair was calculated from this data and ranges from +0.03 to +0.3 kcal/mol over the NaCl range of 2.4-3.5 M. The free energy for the formation of a B-Z junction was nearly constant at +4.0 kcal/mol over the same range of NaCl concentrations. The salt dependence of the free energy of the Z to B transition indicates preferential Na+ binding to the Z form and that there is a net release of Na+ upon conversion of a base pair from the Z to the B conformation. The energetically unfavorable Z to B transition was found by this analysis to be driven by coupling to the energetically favorable interaction of daunomycin with B form DNA. In 3.5 M NaCl, for example, the free energy change for the overall reaction (Z DNA base pairs) + (daunomycin) in equilibrium with (right-handed complex) is -7.0 kcal/mol, nearly all of which is contributed by the binding of drug to B DNA. Analysis using the allosteric model also shows that the number of base pairs converted from the Z to the B conformation per bound drug molecule is salt dependent and provides evidence that drug molecules partition into regions of the polymer in the right-handed conformation.     

Effect of base stacking on the relative thermodynamic stability of oligonucleotide complexes: a spectroscopic study

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be -1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Base stacking and unstacking as determined from a DNA decamer containing a fluorescent base

Time-resolved fluorescence decay of a single-stranded DNA decamer d(CTGAAT5CAG), where d5 is the fluorescent base 1-(beta-D-2'-deoxyribosyl)-5-methyl-2-pyrimidinone, was measured and analyzed at several temperatures. The d5 base in the decamer is resolved into three states according to their fluorescence decay lifetime characteristics and temperature dependence of their associated amplitudes: fully extended and completely unstacked state, loosely associated state, and fully stacked state. These states are in slow exchange compared to their fluorescence decay rates. The population of the fully extended and completely unstacked state is small and decreases further with increasing temperature. The loosely associated state, whose fluorescence can still be efficiently quenched by other DNA bases, occupies a large portion of the conventionally defined unstacked state. Stacking enthalpy and entropy for the d5 base with thymine or cytosine bases in the DNA decamer are calculated to be -6.6 kcal/mol and -22 cal/mol.K, respectively. This work shows that fluorescent bases in DNA can be useful to the study of local conformations of bases.     

Thermodynamic analysis of 5' and 3' single- and 3' double-nucleotide overhangs neighboring wobble terminal base pairs

Thermodynamic parameters are reported for duplex formation of 40 self-complementary RNA duplexes containing wobble terminal base pairs with all possible 3′ single and double-nucleotide overhangs, mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest neighbor analysis, the addition of a single 3′ dangling nucleotide increases the stability of duplex formation up to 1 kcal/mol in a sequence-dependent manner. The addition of a second dangling nucleotide increases the stability of duplexes closed with wobble base pairs in an idiosyncratic manner. The results allow for the development of a nearest neighbor model, which improves the predication of free energy and melting temperature for duplexes closed by wobble base pairs with 3′ single or double-nucleotide overhangs. Phylogenetic analysis of naturally occurring miRNAs was performed. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent on the orientation of the GU closing base pair rather than the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for the 5′ single terminal overhangs adjacent to wobble closing base pairs are also presented.     

DNA polymerase insertion fidelity. Gel assay for site-specific kinetics

A quantitative assay based on gel electrophoresis is described to measure nucleotide insertion kinetics at an arbitrary DNA template site. The assay is used to investigate kinetic mechanisms governing the fidelity of DNA synthesis using highly purified Drosophila DNA polymerase alpha holoenzyme complex and M13 primer-template DNA. Km and Vmax values are reported for correct insertion of A and misinsertion of G, C, and T opposite a single template T site. The misinsertion frequencies are 2 X 10(-4) for G-T and 5 X 10(-5) for both C-T and T-T relative to normal A-T base pairs. The dissociation constant of the polymerase-DNA-dNTP complex, as measured by Km, plays a dominant role in determining the rates of forming right and wrong base pairs. Compared with Km for insertion of A opposite T (3.7 +/- 0.7 microM), the Km value is 1100-fold greater for misinsertion of G opposite T (4.2 +/- 0.4 mM), and 2600-fold greater for misinsertion of either C or T opposite T (9.8 +/- 4.2 mM). These Km differences indicate that in the enzyme binding site the stability of A-T base pairs is 4.3 kcal/mol greater than G-T pairs and 4.9 kcal/mol greater than C-T or T-T pairs. In contrast to the large differences in Km, differences in Vmax are relatively small. There is only a 4-fold reduction in Vmax for insertion of G opposite T and an 8-fold reduction for C or T opposite T, compared with the correct insertion of A. For the specific template T site investigated, the nucleotide insertion fidelity for Drosophila polymerase alpha seems to be governed primarily by a Km discrimination mechanism.     

Thermodynamic contribution and nearest-neighbor parameters of pseudouridine-adenosine base pairs in oligoribonucleotides

Pseudouridine (Ψ) is the most common noncanonical nucleotide present in naturally occurring RNA and serves a variety of roles in the cell, typically appearing where structural stability is crucial to function. Ψ residues are isomerized from native uridine residues by a class of highly conserved enzymes known as pseudouridine synthases. In order to quantify the thermodynamic impact of pseudouridylation on U-A base pairs, 24 oligoribonucleotides, 16 internal and eight terminal Ψ-A oligoribonucleotides, were thermodynamically characterized via optical melting experiments. The thermodynamic parameters derived from two-state fits were used to generate linearly independent parameters for use in secondary structure prediction algorithms using the nearest-neighbor model. On average, internally pseudouridylated duplexes were 1.7 kcal/mol more stable than their U-A counterparts, and terminally pseudouridylated duplexes were 1.0 kcal/mol more stable than their U-A equivalents. Due to the fact that Ψ-A pairs maintain the same Watson-Crick hydrogen bonding capabilities as the parent U-A pair in A-form RNA, the difference in stability due to pseudouridylation was attributed to two possible sources: the novel hydrogen bonding capabilities of the newly relocated imino group as well as the novel stacking interactions afforded by the electronic configuration of the Ψ residue. The newly derived nearest-neighbor parameters for Ψ-A base pairs may be used in conjunction with other nearest-neighbor parameters for accurately predicting the most likely secondary structure of A-form RNA containing Ψ-A base pairs.     

Solution structure and dynamics of DNA duplexes containing the universal base analogues 5-nitroindole and 5-nitroindole 3-carboxamide

Universal bases hybridize with all other natural DNA or RNA bases, and have applications in PCR and sequencing. We have analysed by nuclear magnetic resonance spectroscopy the structure and dynamics of three DNA oligonucleotides containing the universal base analogues 5-nitroindole and 5-nitroindole-3-carboxamide. In all systems studied, both the 5-nitroindole nucleotide and the opposing nucleotide adopt a standard anti conformation and are fully stacked within the DNA duplex. The 5-nitroindole bases do not base pair with the nucleotide opposite them, but intercalate between this base and an adjacent Watson–Crick pair. In spite of their smooth accommodation within the DNA double-helix, the 5-nitroindole-containing duplexes exist as a dynamic mixture of two different stacking configurations exchanging fast on the chemical shift timescale. These configurations depend on the relative intercalating positions of the universal base and the opposing base, and their exchange implies nucleotide opening motions on the millisecond time range. The structure of these nitroindole-containing duplexes explains the mechanism by which these artificial moieties behave as universal bases.     

A search for base analogs to enhance third-strand binding to 'inverted' target base pairs of triplexes in the pyrimidine/parallel motif.

Eight base analogs were tested as third strand residues in otherwise homopyrimidine strands opposite each of the 'direct' (A.T and G.C) and 'inverted' (T.A and C.G) Watson-Crick base pairs, using UV melting profiles to assess triplex stability. The target duplexes contained 20 A.T base pairs and a central test base pair X.Y, while the third strand contained 20 T residues and a central Z test base. Z included 5-bromo-uracil, 5-propynyluracil, 5-propynylcytosine, 5-methyl-cytosine, 5-bromocytosine, hypoxanthine, 2-amino-purine and 2,6-diaminopurine. Some of the base analogs enhanced third strand binding to the target duplex with one or other 'inverted' central base pair relative to the binding afforded by any of the canonical bases. Other analogs did the same for the duplexes with the 'direct' target pairs. The increasing order of triplex stabilization by these base analogs is: opposite the 'inverted' base pairs, for T.A, A < C < 5-pC < 5-pU < T < 5-BrC < 5-meC < 5-BrU < 2-AP < 2,6-DAP < Hy < G, for C.G, 2-AP < A < Hy < G < 5-pC < 5-BrC < 5-meC < C < 2,6-DAP < T < 5-BrU < 5-pU; opposite the 'direct' base pairs, for A.T, 2-AP < A < 5-meC < C < G < Hy < 2,6-DAP < 5-pU < T = 5-BrU < 5-BrC < 5-pC, for G.C, G < 2,6-DAP < 2-AP < A < Hy < T < 5-BrU < 5-pU < 5-pC < 5-BrC < C < 5-meC.     

The thermodynamics of 3'-terminal pyrene and guanosine for the design of isoenergetic 2'-O-methyl-RNA-LNA chimeric oligonucleotide probes of RNA structure

To facilitate design of short isoenergetic hybridization probes for RNA, we report the influence of adding 5'- or 3'-terminal 2'-O-methylguanosine (GM), LNA-guanosine (GL), or 3'-terminal pyrene pseudo-nucleotide (PPN) on the thermodynamic stability of 2'-O-methyl-RNA/RNA (2'-O-Me-RNA/RNA) duplexes with sequences 5'CMGMGMCMAM/3'AAXGCCGUXAA, where X is A, C, G, or U. A 3'-terminal GM or GL added to the 2'-O-Me-RNA strand to form a G-A, G-G or G-U mismatch enhances thermodynamic stability (DeltaDeltaG degrees 37) of the 2'-O-Me-RNA/RNA duplexes on average by 0.7 and 1.5 kcal/mol, respectively. A 3'-terminal GM or GL in a GM-C or GL-C pair stabilizes the 2'-O-Me-RNA/RNA duplex by 2.6 and 3.4 kcal/mol, respectively. A 5'-terminal GM or GL in a G-A or G-G mismatch provided less stabilization in comparison with a 3'-terminal G-A or G-G mismatch, but more stabilization in a G-C or G-U pair. In contrast to guanosine derivatives, pyrene residue (P) as PPN at the 3'-terminal position enhances thermodynamic stability of the 2'-O-Me-RNA/RNA duplexes on average by 2.3 +/- 0.1 kcal/mol, relatively independent of the type of ribonucleotide placed in the opposite strand. The thermodynamic data can be applied to design 2'-O-Me-RNA/RNA duplexes with enhanced thermodynamic stability that is also sequence independent. This is useful for design of hybridization probes to interrogate RNA structure and/or expression by microarray and other methods.     

Thermodynamics of unpaired terminal nucleotides on short RNA helixes correlates with stacking at helix termini in larger RNAs.

Free energies for stacking of unpaired nucleotides (dangling ends) at the termini of oligoribonucleotide Watson-Crick helixes (DeltaG(0)37,stack) depend on sequence for 3' ends but are always small for 5' ends. Here, these free energies are correlated with stacking at helix termini in a database of 34 RNA structures determined by X-ray crystallography and NMR spectroscopy. Stacking involving GA pairs is considered separately. A base is categorized as stacked by its distance from (相似文献   

Intrinsic conformational properties of deoxyribonucleosides: implicated role for cytosine in the equilibrium among the A, B, and Z forms of DNA.

Structural properties of biomolecules are dictated by their intrinsic conformational energetics in combination with environmental contributions. Calculations using high-level ab initio methods on the deoxyribonucleosides have been performed to investigate the influence of base on the intrinsic conformational energetics of nucleosides. Energy minima in the north and south ranges of the deoxyribose pseudorotation surfaces have been located, allowing characterization of the influence of base on the structures and energy differences between those minima. With all bases, chi values associated with the south energy minimum are lower than in canonical B-DNA, while chi values associated with the north energy minimum are close to those in canonical A-DNA. In deoxycytidine, chi adopts an A-DNA conformation in both the north and south energy minima. Energy differences between the A and B conformations of the nucleosides are <0.5 kcal/mol in the present calculations, except with deoxycytidine, where the A form is favored by 2.3 kcal/mol, leading the intrinsic conformational energetics of GC basepairs to favor the A form of DNA by 1.5 kcal/mol as compared with AT pairs. This indicates that the intrinsic conformational properties of cytosine at the nucleoside level contribute to the A form of DNA containing predominately GC-rich sequences. In the context of a B versus Z DNA equilibrium, deoxycytidine favors the Z form over the B form by 1.6 kcal/mol as compared with deoxythymidine, suggesting that the intrinsic conformational properties of cytosine also contribute to GC-rich sequences occurring in Z DNA with a higher frequency than AT-rich sequences. Results show that the east pseudorotation energy barrier involves a decrease in the furanose amplitude and is systematically lower than the inversion barrier, with the energy differences influenced by the base. Energy barriers going from the south (B form) sugar pucker to the east pseudorotation barrier are lower in pyrimidines as compared with purines, indicating that the intrinsic conformational properties associated with base may also influence the sugar pseudorotational population distribution seen in DNA crystal structures and the kinetics of B to A transitions. The present work provides evidence that base composition, in addition to base sequence, can influence DNA conformation.