Thermodynamics of stabilization of RNA pseudoknots by cobalt(III) hexaammine

Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability

The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.

Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.

Molecular dynamics simulation of RNA pseudoknot unfolding pathway

Many biological functions of RNA molecules are related to their pseudoknot structures. It is significant for predicting the structure and function of RNA that learning about the stability and the process of RNA pseudoknot folding and unfolding. The structural features of mouse mammary tumor virus (MMTV) RNA pseudoknot in different ion concentration, the unfolding process of the RNA pseudoknot, and the two hairpin helices that constitute the RNA pseudoknot were studied with all atom molecule dynamics simulation method in this paper. We found that the higher cation concentration can cause structure of the RNA molecules more stable, and ions played an indispensable role in keeping the structure of RNA molecules stable; the unfolding process of hairpin structure was corresponding to the antiprocess of its folding process. The main pathway of pseudoknot unfolding was that the inner base pair opened first, and then, the two helices, which formed the RNA pseudoknot opened decussately, while the folding pathway of the RNA pseudoknot was a helix folding after formation of the other helix. Therefore, the unfolding process of RNA pseudoknot is different from the antiprocess of its folding process, and the unfolding process of each helix in the RNA pseudoknot is similar to the hairpin structure’s unfolding process, which means that both are the unzipping process.

The 3' splice site of influenza A segment 7 mRNA can exist in two conformations: a pseudoknot and a hairpin

The 3′ splice site of influenza A segment 7 is used to produce mRNA for the M2 ion-channel protein, which is critical to the formation of viable influenza virions. Native gel analysis, enzymatic/chemical structure probing, and oligonucleotide binding studies of a 63 nt fragment, containing the 3′ splice site, key residues of an SF2/ASF splicing factor binding site, and a polypyrimidine tract, provide evidence for an equilibrium between pseudoknot and hairpin structures. This equilibrium is sensitive to multivalent cations, and can be forced towards the pseudoknot by addition of 5 mM cobalt hexammine. In the two conformations, the splice site and other functional elements exist in very different structural environments. In particular, the splice site is sequestered in the middle of a double helix in the pseudoknot conformation, while in the hairpin it resides in a two-by-two nucleotide internal loop. The results suggest that segment 7 mRNA splicing can be controlled by a conformational switch that exposes or hides the splice site.

Atomistic analysis of pseudoknotted RNA unfolding

As important functional structures, RNA pseudoknots provide excellent models for studying the interplay between secondary and tertiary structures and the roles of triplexes, noncanonical interactions, and coaxial stacking in the folding/unfolding process. Here we report a first atomistic and statistical analysis of the unfolding of the pseudoknot within gene 32 mRNA of bacteriophage T2. Multiple unfolding pathways, diverse transition states, and various intermediate structures were observed. Water molecules were found to be coupled with the unfolding process via the expulsion or concurrent mechanism.

Salt dependence of nucleic acid hairpin stability

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na(+) and Mg(2+) in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na(+) and Mg(2+) concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.

Distinct reaction pathway promoted by non-divalent-metal cations in a tertiary stabilized hammerhead ribozyme

Divalent ion sensitivity of hammerhead ribozymes is significantly reduced when the RNA structure includes appropriate tertiary stabilization. Therefore, we investigated the activity of the tertiary stabilized "RzB" hammerhead ribozyme in several nondivalent ions. Ribozyme RzB is active in spermidine and Na(+) alone, although the cleavage rates are reduced by more than 1,000-fold relative to the rates observed in Mg(2+) and in transition metal ions. The trivalent cobalt hexammine (CoHex) ion is often used as an exchange-inert analog of hydrated magnesium ion. Trans-cleavage rates exceeded 8 min(-1) in 20 mM CoHex, which promoted cleavage through outersphere interactions. The stimulation of catalysis afforded by the tertiary structural interactions within RzB does not require Mg(2+), unlike other extended hammerhead ribozymes. Site-specific interaction with at least one Mg(2+) ion is suggested by CoHex competition experiments. In the presence of a constant, low concentration of Mg(2+), low concentrations of CoHex decreased the rate by two to three orders of magnitude relative to the rate in Mg(2+) alone. Cleavage rates increased as CoHex concentrations were raised further, but the final fraction cleaved was lower than what was observed in CoHex or Mg(2+) alone. These observations suggest that Mg(2+) and CoHex compete for binding and that they cause misfolded structures when they are together. The results of this study support the existence of an alternate catalytic mechanism used by nondivalent ions (especially CoHex) that is distinct from the one promoted by divalent metal ions, and they imply that divalent metals influence catalysis through a specific nonstructural role.

Solution structure of a luteoviral P1-P2 frameshifting mRNA pseudoknot

A hairpin-type messenger RNA pseudoknot from pea enation mosaic virus RNA1 (PEMV-1) regulates the efficiency of programmed -1 ribosomal frameshifting. The solution structure and 15N relaxation rates reveal that the PEMV-1 pseudoknot is a compact-folded structure composed almost entirely of RNA triple helix. A three nucleotide reverse turn in loop 1 positions a protonated cytidine, C(10), in the correct orientation to form an A((n-1)).C(+).G-C(n) major groove base quadruple, like that found in the beet western yellows virus pseudoknot and the hepatitis delta virus ribozyme, despite distinct structural contexts. A novel loop 2-loop 1 A.U Hoogsteen base-pair stacks on the C(10)(+).G(28) base-pair of the A(12).C(10)(+).G(28)-C(13) quadruple and forms a wedge between the pseudoknot stems stabilizing a bent and over-rotated global conformation. Substitution of key nucleotides that stabilize the unique conformation of the PEMV-1 pseudoknot greatly reduces ribosomal frameshifting efficacy.

Effects of cobalt hexammine on folding and self-cleavage of the Neurospora VS ribozyme

We have investigated the effects of Co(NH3)6(3+), an analog of hexahydrated Mg2+, on folding and catalysis of the Neurospora VS ribozyme. Most of the metal ion-induced changes detected by chemical modification structure probing in either metal ion are similar, but occur at approximately 33-fold lower concentrations of Co(NH3)6(3+) than Mg2+. However, Co(NH3)6(3+) is not as effective at inducing two functionally important structural changes: stabilizing the pseudoknot interaction between loops I and V, and rearranging the secondary structure of helix Ib. Comparison of the folding of the precursor and the downstream cleavage product, which lacks helix Ia, shows that helix Ia inhibits stable pseudoknot formation and rearrangement of helix Ib. The VS ribozyme does not self-cleave with Co(NH3)6(3+) as the sole polyvalent cation; however, mixed-metal kinetic experiments show that Co(NH3)6(3+) does not inhibit Mg2+-induced self-cleavage. In contrast, at sub-saturating concentrations of Mg2+, Co(NH3)6(3+) increases the rate of Mg2+-induced self-cleavage, indicating that Co(NH3)6(3+) contributes to the functionally relevant folding of the VS ribozyme.

Structure, stability and function of RNA pseudoknots involved in stimulating ribosomal frameshifting

Programmed -1 ribosomal frameshifting has become the subject of increasing interest over the last several years, due in part to the ubiquitous nature of this translational recoding mechanism in pathogenic animal and plant viruses. All cis-acting frameshift signals encoded in mRNAs are minimally composed of two functional elements: a heptanucleotide "slippery sequence" conforming to the general form X XXY YYZ, followed by an RNA structural element, usually an H-type RNA pseudoknot, positioned an optimal number of nucleotides (5 to 9) downstream. The slippery sequence itself promotes a low level ( approximately 1 %) of frameshifting; however, downstream pseudoknots stimulate this process significantly, in some cases up to 30 to 50 %. Although the precise molecular mechanism of stimulation of frameshifting remains poorly understood, significant advances have been made in our knowledge of the three-dimensional structures, thermodynamics of folding, and functional determinants of stimulatory RNA pseudoknots derived from the study of several well-characterized frameshift signals. These studies are summarized here and provide new insights into the structural requirements and mechanism of programmed -1 ribosomal frameshifting.

Energetics of a strongly pH dependent RNA tertiary structure in a frameshifting pseudoknot

Retroviruses employ -1 translational frameshifting to regulate the relative concentrations of structural and non-structural proteins critical to the viral life cycle. The 1.6 A crystal structure of the -1 frameshifting pseudoknot from beet western yellows virus reveals, in addition to Watson-Crick base-pairing, many loop-stem RNA tertiary structural interactions and a bound Na(+). Investigation of the thermodynamics of unfolding of the beet western yellows virus pseudoknot reveals strongly pH-dependent loop-stem tertiary structural interactions which stabilize the molecule, contributing a net of DeltaH approximately -30 kcal mol(-1) and DeltaG degrees (37) of -3.3 kcal mol(-1) to a total DeltaH and DeltaG degrees (37) of -121 and -16 kcal mol(-1), respectively, at pH 6.0, 0.5 M K(+) by DSC. Characterization of mutant RNAs supports the presence of a C8(+).G12-C26 loop 1-stem 2 base-triple (pK(a)=6.8), protonation of which contributes nearly -3.5 kcal mol(-1) in net stability in the presence of a wild-type loop 2. Substitution of the nucleotides in loop 2 with uridine bases, which would eliminate the minor groove triplex, destroys pseudoknot formation. An examination of the dependence of the monovalent ion and type on melting profiles suggests that tertiary structure unfolding occurs in a manner quantitatively consistent with previous studies on the stabilizing effects of K(+), NH(4)(+) and Na(+) on other simple duplex and pseudoknotted RNAs.

Equilibrium unfolding pathway of an H-type RNA pseudoknot which promotes programmed -1 ribosomal frameshifting

The equilibrium unfolding pathway of a 41-nucleotide frameshifting RNA pseudoknot from the gag-pro junction of mouse intracisternal A-type particles (mIAP), an endogenous retrovirus, has been determined through analysis of dual optical wavelength, equilibrium thermal melting profiles and differential scanning calorimetry. The mIAP pseudoknot is an H-type pseudoknot proposed to have structural features in common with the gag-pro frameshifting pseudoknots from simian retrovirus-1 (SRV-1) and mouse mammary tumor virus (MMTV). In particular, the mIAP pseudoknot is proposed to contain an unpaired adenosine base at the junction of the two helical stems (A15), as well as one in the middle of stem 2 (A35). A mutational analysis of stem 1 hairpins and compensatory base-pair substitutions incorporated into helical stem 2 was used to assign optical melting transitions to molecular unfolding events. The optical melting profile of the wild-type RNA is most simply described by four sequential two-state unfolding transitions. Stem 2 melts first in two closely coupled low-enthalpy transitions at low tmin which the stem 3' to A35, unfolds first, followed by unfolding of the remainder of the helical stem. The third unfolding transition is associated with some type of stacking interactions in the stem 1 hairpin loop not present in the pseudoknot. The fourth transition is assigned to unfolding of stem 1. In all RNAs investigated, DeltaHvH approximately DeltaHcal, suggesting that DeltaCpfor unfolding is small. A35 has the thermodynamic properties expected for an extrahelical, unpaired nucleotide. Deletion of A15 destabilizes the stem 2 unfolding transition in the context of both the wild-type and DeltaA35 mutant RNAs only slightly, by DeltaDeltaG degrees approximately 1 kcal mol-1(at 37 degrees C). The DeltaA15 RNA is considerably more susceptible to thermal denaturation in the presence of moderate urea concentrations than is the wild-type RNA, further evidence of a detectable global destabilization of the molecule. Interestingly, substitution of the nine loop 2 nucleotides with uridine residues induces a more pronounced destabilization of the molecule (DeltaDeltaG degrees approximately 2.0 kcal mol-1), a long-range, non-nearest neighbor effect. These findings provide the thermodynamic basis with which to further refine the relationship between efficient ribosomal frameshifting and pseudoknot structure and stability.