Free energy of imperfect nucleic acid helices. 3. Small internal loops resulting from mismatches

Physical studies of enzymioally synthesized oligoribonucleotides of defined sequence are used to evaluate quantitatively the destabilizing influence of mismatched bases in a double helix. The series (A-)₄G(-C)_n(-U)₄, N = 1 to 6, exist as imperfect dimer helices when N is equal to or less than 4, and as monomolecular hairpin helices when N is 5 and 6. Internal loops become progressively more destabilizing as their size increases from 2 to 4 to 6 nucleotides resulting from 1, 2 and 3 consecutive mismatched base pairs. However, the stability of a helix will generally be greater if a given number of mismatched pairs occur consecutively rather than in isolation from one another.These data may be used for improved calculations of stability of RNA secondary structure, to estimate the frequency of structural fluctuations in a double helix and to assess the stability of modified polynucleotide helices. An unmodified double helix of one million randomly arranged base pairs should contain on the time average approximately 10 G.C and 500 A.U pairs in non-hydrogen bonded, unstacked conformations at 25 °C. Our estimate of the effect of mismatching on T_m values of high polymers is less precise because of the long temperature extrapolation required. However, we estimate that DNA or RNA treated with mutagens which interrupt up to 20% of the nucleotide pairs will show a drop of about 1.2 deg. C in melting temperature with each unit per cent of modification.     

Determination of thermodynamic parameters for HIV DIS type loop-loop kissing complexes

The HIV-1 type dimerization initiation signal (DIS) loop was used as a starting point for the analysis of the stability of Watson–Crick (WC) base pairs in a tertiary structure context. We used ultraviolet melting to determine thermodynamic parameters for loop–loop tertiary interactions and compared them with regular secondary structure RNA helices of the same sequences. In 1 M Na⁺ the loop–loop interaction of a HIV-1 DIS type pairing is 4 kcal/mol more stable than its sequence in an equivalent regular and isolated RNA helix. This difference is constant and sequence independent, suggesting that the rules governing the stability of WC base pairs in the secondary structure context are also valid for WC base pairs in the tertiary structure context. Moreover, the effect of ion concentration on the stability of loop–loop tertiary interactions differs considerably from that of regular RNA helices. The stabilization by Na⁺ and Mg²⁺ is significantly greater if the base pairing occurs within the context of a loop–loop interaction. The dependence of the structural stability on salt concentration was defined via the slope of a T_m/log [ion] plot. The short base-paired helices are stabilized by 8°C/log [Mg²⁺] or 11°C/log [Na⁺], whereas base-paired helices forming tertiary loop–loop interactions are stabilized by 16°C/log [Mg²⁺] and 26°C/log [Na⁺]. The different dependence on ionic strength that is observed might reflect the contribution of specific divalent ion binding to the preformation of the hairpin loops poised for the tertiary kissing loop–loop contacts.     

The stability of RNA hairpin loops containing A-U-G: An-U-G-Um

RNA oligomers with the sequence A_n-U-G-U_m, n = 7–9, m = 5–10, have been synthesized and found to form hairpin loops in 21 mM or 101 mM Na⁺. The hairpin loops displayed melting temperatures 13°–29°C greater than that of the hairpin loop A₆-C₈-U₆ in the same solvent. The increased stability of these hairpin loops is attributed to the presence of the trinucleotide A-U-G in the loop. Circular dichroism (CD) spectra were taken of the hairpin loops A₇-U-G-U_6,7,8, A₈-U-G-U_6,7,8, and A₉-U-G-U₆ in 21 mM Na⁺, and compared with circular dichroism spectra of A₆-U₆ in 1 M Na⁺. Difference spectra were calculated between each A_n-U-G-U_m and 11.5 mM (nucleotide) A₆-U₆ at similar temperatures and identical singlestrand fractions to give the “experimental” CD spectrum of the unbase-paired nucleotides in the loop, assuming, five, six, or seven base pairs. CD spectra were calculated for each of the assumed unbase-paired sequences using the measured CD spectra of ApA, A-U-G-, A, and U, and compared with the experimental spectra. The best agreement was found for hairpin-loop models containing five base pairs and five to eight unbase-paired nucleotides in the loop.     

Hydrolysis of polynucleotides and the characterization of their secondary structure. A theoretical study

1. Single-stranded RNA may be regarded as an assembly of L hairpin loops each stabilized by N base pairs and each containing b unpaired residues; one loop is connected to another by c residues. 2. A theory based on the statistics of the random degradation of linear polymers was developed to relate N, b and c with the probability, p, of hydrolysing a diesterified phosphate bond. 3. The number of residues per hairpin loop, which is 2N+b, is related to the fraction, f, of the original loops remaining intact by the equation: 2N+b=logf/log(1–p). 4. The theory was extended to show that the number of residues per loop may be evaluated by fractionating the RNA after hydrolysis and examining the secondary structure of each fraction. Fragments that are shorter than the hairpin loop cannot reproduce the original secondary structure. The probability that a fragment will form an intact loop increases most rapidly for fragments of between 2N+b and 2(2N+b)+c residues. 5. The probability of producing a fragment capable of forming one, and only one, hairpin loop was related to N, b and c.     

Loop dependence of the stability and dynamics of nucleic acid hairpins

Hairpin loops are critical to the formation of nucleic acid secondary structure, and to their function. Previous studies revealed a steep dependence of single-stranded DNA (ssDNA) hairpin stability with length of the loop (L) as ~L^{8.5 ± 0.5}, in 100 mM NaCl, which was attributed to intraloop stacking interactions. In this article, the loop-size dependence of RNA hairpin stabilities and their folding/unfolding kinetics were monitored with laser temperature-jump spectroscopy. Our results suggest that similar mechanisms stabilize small ssDNA and RNA loops, and show that salt contributes significantly to the dependence of hairpin stability on loop size. In 2.5 mM MgCl₂, the stabilities of both ssDNA and RNA hairpins scale as ~L^{4 ± 0.5}, indicating that the intraloop interactions are weaker in the presence of Mg²⁺. Interestingly, the folding times for ssDNA hairpins (in 100 mM NaCl) and RNA hairpins (in 2.5 mM MgCl₂) are similar despite differences in the salt conditions and the stem sequence, and increase similarly with loop size, ~L^{2.2 ± 0.5} and ~L^{2.6 ± 0.5}, respectively. These results suggest that hairpins with small loops may be specifically stabilized by interactions of the Na⁺ ions with the loops. The results also reinforce the idea that folding times are dominated by an entropic search for the correct nucleating conformation.     

Stabilities of HIV-1 DIS type RNA loop-loop interactions in vitro and in vivo

RNA loop–loop interactions are a prevalent motif in the formation of tertiary structure and are well suited to trigger molecular recognition between RNA molecules. We determined the stabilities of several loop–loop interactions with a constant 6 bp core sequence and varying unpaired flanking nucleotides and found that the flanking bases have a strong influence on the stability and ion dependence of the kissing complex. In general, the stabilities determined in 1 M Na⁺ are equivalent to those in the presence of near physiological Mg²⁺ concentrations. Therefore we further tested whether the stabilities determined in vitro and within yeast cells correlate, using a recently developed yeast RNA-hybrid system. For the majority of the loop types analyzed here, the melting temperatures determined in vitro are in good agreement with the relative β-galactosidase activity in yeast cells, showing that data derived from in vitro measurements reflect in vivo properties. The most stable interactions are the naturally occurring HIV-1 DIS MAL and LAI derived loops with the motif (5′ A^A/_GN₆A 3′), emphasizing the crucial role of stable kissing complexes in HIV genome dimerization.     

Determination of the secondary structure of group II bulge loops using the fluorescent probe 2-aminopurine

Eleven RNA hairpins containing 2-aminopurine (2-AP) in either base-paired or single nucleotide bulge loop positions were optically melted in 1 M NaCl; and, the thermodynamic parameters ΔH°, ΔS°, ΔG°₃₇, and T_M for each hairpin were determined. Substitution of 2-AP for an A (adenosine) at a bulge position (where either the 2-AP or A is the bulge) in the stem of a hairpin, does not affect the stability of the hairpin. For group II bulge loops such as AA/U, where there is ambiguity as to which of the A residues is paired with the U, hairpins with 2-AP substituted for either the 5′ or 3′ position in the hairpin stem have similar stability. Fluorescent melts were performed to monitor the environment of the 2-AP. When the 2-AP was located distal to the hairpin loop on either the 5′ or 3′ side of the hairpin stem, the change in fluorescent intensity upon heating was indicative of an unpaired nucleotide. A database of phylogenetically determined RNA secondary structures was examined to explore the presence of naturally occurring bulge loops embedded within a hairpin stem. The distribution of bulge loops is discussed and related to the stability of hairpin structures.     

Thermodynamics of melting of the circular dumbbell d〈pCGC-TT-GCG-TT〉

The conformational behavior of DNA minihairpin loops is sensitive to the directionality of the base pair that closes the loop. Especially tailored circular dumbbells, consisting of a stem of three Watson–Crick base pairs capped on each side with a minihairpin loop, serve as excellent model compounds by means of which deeper insight is gained into the relative stability and melting properties of hairpin loops that differ only in directionality of the closing pair: C-G vs G-C. For this reason the thermodynamic properties of the circular DNA decamers 5′-d〈pCGC-TT-GCG-TT〉-3′( I ) and reference compounds 5′-d〈pGGC-TT-GCC-TT≤-3′( II ) and 5′-d(GCG-TC-CGC)-3′( III ) are studied by means of nmr spectroscopy. Molecules I and II adopt dumbbell structures closed on both sides by a two-membered hairpin hop. At low temperature I consists of a mixture of two slowly exchanging forms, denoted L2L2 and L2L4 . The low-temperature L2L2 form is the fully intact minihairpin structure with three Watson–Crick C-G base pairs. The high-temperature form, L2L4 ,contains a partially disrupted closing G-C base pair in the 5′-GTTC-3′ loop, with the cytosine base placed in a syn orientation. The opposite 5′-CTTG-3′ loop remains stable. A study of the noncircular hairpin structure III shows similar conformational behavior for the 5′-GTTC-3′ loop as found in I a syn orientation for C(6) and two slowly exchanging imino proton signals for G(3). The melting point T_m of II was estimated to lie above 365 K. The T_m value of the duplex stem and the 5′-CTTG-3′ loop of the L2L4 form ofIis 352 ± 2 K. The ΔH° is calculated as ?89 ± 10 kJ/mol. The T_m value determined for the individual residues of the 5′-GTTC-3′ loop lies 4°–11° lower. The enthalpy ΔH° of melting the thymine residues in the 5′-GTTC-3′ loop is calculated to be -61± 7 kJ/mol. Thermodynamic data of the equilibrium between the slowly exchanging two- and four-membered loop conformers of I reveal an upper limit for ΔH° of +30 kJ/mol in going from a two-memberedto a four-membered loop, in agreement with the enthalpy difference of +28 k.j/mol between the two loops at the T_m midpoint. For hairpin III the upper limit for ΔH° going from a two-membered to a four-membered loop amounts to ±21 kJ/mol. The mutual exchange rate between the L2 and L4 form in III is estimated as 13.6 s^?1. Our results clearly suggest that small four-way DNA junctions(model for immobilized Holliday junctions) can be designed that consist of a single DNA strandthat features -CTTG-caps on three of the four arms of the junction. © 1995 John Wiley & Sons, Inc.     

Thermodynamic and kinetic parameters of an oligonucleotide hairpin helix

The thermodynamics of the hairpin helix-single strand transition of A₆C₆U₆ has been analyzed by a staggering zipper model with consideration of single strand stacking. This analysis yields an enthalpy change of +11 kcal/mole for the formation of a first, isolated base pair. The stability constant of a first (intramolecular) base pair in A₆C₆U₆ is around 2 × 1O^?5 at 25°C, whereas a first (intermoleciilar) base pair in an A₆ · U₆ helix is characterised by a stability constant of about 4 × 10^?3M^?1 (25°C, extrapolated from A_n · V_n oligomer measurements). These data indicate a destabilizing effect of the C₆ loop.The rate constant of hairpin helix formation is 2 to 3 × 10⁴ sec^?1 associated with an activation enthalpy of +2.5 kcal/mote. The rate of helix dissociation of the A₆C₆U₆ hairpin is in the range of 10³ to lO⁵ sec^?1 with an activation enthalpy of 21 $\text{kcal}\text{mole}$. A comparison with the kinetic parameters obtained for A · U oligomer helices shows a specific influence of the C₆ loop due to the stacking tendency of the cytosine residues. This intluence is preferentially reflected in the relatively low value of the rate constant of helix formation.     

On loop folding in nucleic acid hairpin-type structures

In a series of studies, combining NMR, optical melting and T-jump experiments, it was found that DNA hairpins display a maximum stability when the loop part of the molecule comprises four or five nucleotide residues. This is in contrast with the current notion based on RNA hairpin studies, from which it had been established that a maximum hairpin stability is obtained for six or seven residues in the loop. Here we present a structural model to rationalize these observations. This model is based on the notion that to a major extent base stacking interactions determine the stability of nucleic acid conformations. The model predicts that loop folding in RNA is characterized by an extension of the base stacking at the 5'-side of the double helix by five or six bases; the remaining gap can then easily be closed by two nucleotides. Conversely, loop folding in DNA is characterized by extending base stacking at the 3'-side of the double helical stem by two or three residues; again bridging of the remaining gap can then be achieved by one or two nucleotides. As an example of loop folding in RNA the anticodon loop of yeast tRNAPhe is discussed. For the DNA hairpin formed by d(ATCCTAT4TAGGAT) it is shown that the loop structure obtained from molecular mechanics calculations obeys the above worded loop folding principles.     

Solution conformation of an RNA hairpin loop.

The hairpin conformation adopted by the RNA sequence 5'GCGAUUUCUGACCGCC3' has been studied by one- and two-dimensional NMR spectroscopy. Exchangeable imino spectra in 60 mM Na+ indicate that the hairpin has a stem of six base pairs (indicated by boldface type) and a loop of three nucleotides. NOESY spectra of nonexchangeable protons confirm the formation of the stem region. The duplex has an A-conformation and contains an A.C apposition; a G.U base pair closes the loop region. The stem nucleotides have C3'-endo sugar conformations, as expected of an A-form duplex, whereas the three loop nucleotides adopt C2'-endo sugar puckers. Stacking within the loop, C8 upon the sugar of U7, stabilizes the structure. The pH dependence of both the exchangeable and nonexchangeable NMR spectra is consistent with the formation of an A+.C base pair, protonated at the N1 position of adenine. The stability of the hairpin was probed by using absorbance melting curves. The hairpin structure with the A+.C base pair is about +2 kcal/mol less stable in free energy at 37 degrees C than the hairpin formed with an A.U pair replacing the A+.C pair.     

The Conformation of a Conserved Stem-Loop Structure in Ribosomal RNA

Abstract

The RNA of small ribosomal subunits contains a conserved stem-loop structure near the 3′ end. Characteristics for the hairpins are: (a) a nine-basepairs stem; (b) a conserved ^A-U _U-Gjunction in the stem; (c) a conserved sequence Gm⁶ ₂Am⁶ ₂A sequence in the loop (except yeast mitochondria and mutants from bacteria). We are using UV-optics, micro-calorimetry and 500 MHz-NMR to investigate fragments of about 50 nucleotides cleaved from the 3′ ends of small ribosomal subunit RNA's by bacteriocins. Our preliminary conclusions are: (1) Dimethylation of the adenines in the loop destabilizes the hairpin because of an increased stacking; (2) melting of the hairpin starts at the ends as well as in the middle at the 5-H junction; (3) basepair substitutions have an unexpectedly large effect on thermal stability.     

Preference for Ribose Over Deoxyribose in Loop-Closing Base Pairs of Extra Stable Nucleic Acid Hairpins

We have investigated the effect of switching ribose to deoxyribose at the closing base-pair of an extra-stable RNA hairpin. Specifically, we studied the sequence 5′-GGAC(UUCG)GUCC, a dodecanucleotide that folds into a well-characterized, “extra stable” RNA hairpin structure. Recently, we showed that hairpins containing a 2′,5′-linked (UUCG) loop instead of the native 3′,5′-linked loop also exhibit extra-stability (Hannoush and Damha, J. Am. Chem. Soc., 2001, 123, 12368–12374). In this article, we show that the ribose units located at the loop-closing positions (i.e., rC ₄ and rG ₉ ) contribute significantly to the stabilization of RNA hairpins, particularly those containing the 3′,5′-UUCG loop. Interestingly, the requirement of rC4 and rG9 is more relaxed for DNA hairpins containing the 2′,5′-UUCG loop and, in fact, they may be replaced altogether (ribose → deoxyribose) without affecting stability. The results broaden our understanding of the behavior of highly stable (UUCG) hairpin loops and how they respond to structural perturbation of the loop-closing base pairs.     

How long is too long? Effects of loop size on G-quadruplex stability

We compared here 80 different sequences containing four tracts of three guanines with loops of variable length (between 1 and 15 bases for unmodified sequences, up to 30 for fluorescently labeled oligonucleotides). All sequences were capable of forming stable quadruplexes, with T_m above physiological temperature in most cases. Unsurprisingly, the melting temperature was systematically lower in sodium than in potassium but the difference between both ionic conditions varied between 1 and >39°C (average difference: 18.3°C). Depending on the sequence context, and especially for G4 sequences involving two very short loops, the third one may be very long without compromising the stability of the quadruplex. A strong inverse correlation between total loop length and T_m was found in K⁺: each added base leads to a 2°C drop in T_m or ∼0.3 kcal/mol loss in ΔG°. The trend was less clear in Na⁺, with a longer than expected optimal loop length (up to 5 nt). This study will therefore extend the sequence repertoire of quadruplex-prone sequences, arguing for a modification of the widely used consensus (maximal loop size of 7 bases).     

A Short Sequence Motif in the 5′ Leader of the HIV-1 Genome Modulates Extended RNA Dimer Formation and Virus Replication

The 5′ leader of the HIV-1 RNA genome encodes signals that control various steps in the replication cycle, including the dimerization initiation signal (DIS) that triggers RNA dimerization. The DIS folds a hairpin structure with a palindromic sequence in the loop that allows RNA dimerization via intermolecular kissing loop (KL) base pairing. The KL dimer can be stabilized by including the DIS stem nucleotides in the intermolecular base pairing, forming an extended dimer (ED). The role of the ED RNA dimer in HIV-1 replication has hardly been addressed because of technical challenges. We analyzed a set of leader mutants with a stabilized DIS hairpin for in vitro RNA dimerization and virus replication in T cells. In agreement with previous observations, DIS hairpin stability modulated KL and ED dimerization. An unexpected previous finding was that mutation of three nucleotides immediately upstream of the DIS hairpin significantly reduced in vitro ED formation. In this study, we tested such mutants in vivo for the importance of the ED in HIV-1 biology. Mutants with a stabilized DIS hairpin replicated less efficiently than WT HIV-1. This defect was most severe when the upstream sequence motif was altered. Virus evolution experiments with the defective mutants yielded fast replicating HIV-1 variants with second site mutations that (partially) restored the WT hairpin stability. Characterization of the mutant and revertant RNA molecules and the corresponding viruses confirmed the correlation between in vitro ED RNA dimer formation and efficient virus replication, thus indicating that the ED structure is important for HIV-1 replication.     

Stability of RNA hairpin loops: A 6 -C m -U 6

The thermodynamics and circular dichroism of a series of A₆-C_m-U₆ (m = 4, 5, 6 or 8) oligoribonucleotides have been studied. These molecules form intramolecular hairpin loops at low temperatures and therefore are useful models for similar structures which occur in larger, natural RNA molecules. The stability of the helix forming the stem of these loops was found to be considerably greater than an intermolecular helix with the same length and composition. The most stable loop is m = 6. The enthalpy for initiation of the loop is unfavorable; it ranges from + 24 kcal, for m = 4 to + 21 kcal, for m = 6. The maximum in stability for the C₆ loop and the large positive enthalpy for loop initiation are in disagreement with expectations from simple theories assuming a Gaussian distribution of end-to-end distances. Loop strain for m = 4 and m = 5 and the unstacking of the cytosines on loop formation are likely physical explanations for these thermodynamic data. The circular dichroism spectrum of cytosine residues in the C₆ and C₈ loops is very similar to the spectrum of single-stranded oligoribocytidylate. However, the cytosine residues in the C₅ loop have a very different circular dichroism spectrum from the corresponding oligo(C₅) spectrum. In accordance with the thermodynamic data, we conclude from the circular dichroism data that the C₅ loop has an altered conformation from the C₅ and C₈ loops.     

Properties and Anti-HIV Activity of Nicked Dumbbell Oligonucleotides

Abstract

We have designed a new type of oligodeoxyribonucleotide. These oligodeoxyribonucleotides form two hairpin loop structures with base pairs (sense and antisense) in the double helical stem at the 3′ and 5′-ends (nicked dumbbell oligonucleotides). The nicked dumbbell oligonucleotides are molecules with free ends that are more resistant to exonuclease attack. Furthermore, the nicked dumbbell oligonucleotide containing phosphorothioate (P=S) bonds in the hairpin loops has increased nuclease resistance, as compared to the unmodified nicked oligonucleotide. The binding of the nicked dumbbell oligonucleotide to RNA is lower than that of a single-stranded DNA. We also describe the anti-HIV activity of nicked dumbbell oligonucleotides.

     

Stability of single-nucleotide bulge loops embedded in a GAAA RNA hairpin stem

Forty-six RNA hairpins containing combinations of 3' or 5' bulge loops and a 3' or 5' fluorescein label were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each hairpin were determined. The bulge loops were of the group I variety, in which the identity of the bulge is known, and the group II variety, in which the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The fluorescein label at either the 3' end or 5' end of the hairpin did not significantly influence the stability of the hairpin. As observed with bulge loops inserted into a duplex motif, the insertion of a bulge loop into the stem of a hairpin loop was destabilizing. The model developed to predict the influence of bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability. Specifically, the influence of the bulge is related to the stability of the hairpin stem distal from the hairpin loop.     

Loopstructures in synthetic oligonucleotides. Hairpin stability and structure studied as a function of loop elongation

The formation of hairpin structures in the homologous, (partly) self-complementary DNA fragments d(ATCCTAT_nTAGGAT), n = 0–7, was studied by means of nuclear magnetic resonance, T-jump and ultra-violet techniques. It is shown that all compounds in the series may adopt hairpin-like conformations, albeit for n < 3 this only occurs to a significant amount at relatively low concentrations (∼ 10μM). For the present series of oligonucleotides, hairpin formation is accompanied by an apparent loop enthalpy significantly different from zero. The stability of the DNA hairpins turns out to be at its maximum for loop lengths of four or five residues, whereas earlier experiments (Tinocoet al., 1973) indicated that loop lengths of six to seven residues are most favourable for RNA hairpins. This is explained by considering the difference in geometry of A-RNA and B-DNA helices.     

A Test and Refinement of Folding Free Energy Nearest Neighbor Parameters for RNA Including N6-Methyladenosine

RNA folding free energy change parameters are widely used to predict RNA secondary structure and to design RNA sequences. These parameters include terms for the folding free energies of helices and loops. Although the full set of parameters has only been traditionally available for the four common bases and backbone, it is well known that covalent modifications of nucleotides are widespread in natural RNAs. Covalent modifications are also widely used in engineered sequences. We recently derived a full set of nearest neighbor terms for RNA that includes N⁶-methyladenosine (m⁶A). In this work, we test the model using 98 optical melting experiments, matching duplexes with or without N⁶-methylation of A. Most experiments place RRACH, the consensus site of N⁶-methylation, in a variety of contexts, including helices, bulge loops, internal loops, dangling ends, and terminal mismatches. For matched sets of experiments that include either A or m⁶A in the same context, we find that the parameters for m⁶A are as accurate as those for A. Across all experiments, the root mean squared deviation between estimated and experimental free energy changes is 0.67 kcal/mol. We used the new experimental data to refine the set of nearest neighbor parameter terms for m⁶A. These parameters enable prediction of RNA secondary structures including m⁶A, which can be used to model how N⁶-methylation of A affects RNA structure.