Stability Of RNA Hairpins Research Articles

CCG•CGG interruptions in high‐penetrance SCA8 families increase RAN translation and protein toxicity

Spinocerebellar ataxia type 8 (SCA8), a dominantly inherited neurodegenerative disorder caused by a CTG•CAG expansion, is unusual because most individuals that carry the mutation do not develop ataxia. To understand the variable penetrance of SCA8, we studied the molecular differences between highly penetrant families and more common sporadic cases (82%) using a large cohort of SCA8 families (n = 77). We show that repeat expansion mutations from individuals with multiple affected family members have CCG•CGG interruptions at a higher frequency than sporadic SCA8 cases and that the number of CCG•CGG interruptions correlates with age at onset. At the molecular level, CCG•CGG interruptions increase RNA hairpin stability, and in cell culture experiments, increase p‐eIF2α and polyAla and polySer RAN protein levels. Additionally, CCG•CGG interruptions, which encode arginine interruptions in the polyGln frame, increase toxicity of the resulting proteins. In summary, SCA8 CCG•CGG interruptions increase polyAla and polySer RAN protein levels, polyGln protein toxicity, and disease penetrance and provide novel insight into the molecular differences between SCA8 families with high vs. low disease penetrance.

Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations.

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χOL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Investigating RNA Hairpin Stability and Protein Binding Using Optical Tweezers

More than 1.2 million people in the US are living with HIV. This RNA virus must be reverse‐transcribed into DNA before it can be integrated into the host genome and subsequently replicated. Retroviral DNA synthesis is generally understood, but questions remain about many identified critical elements. The transactivation response (TAR) RNA hairpin, specifically, inhibits the minus‐strand transfer step in reverse transcription and thereby regulates the life cycle of the HIV‐1 retrovirus. Before proceeding with reverse transcription, the TAR hairpin must be destabilized: chaperone assistance from the nucleocapsid (NC) protein can stimulate this rate‐limiting step 3000‐fold. Recent interest in the folding landscape of the TAR hairpin has revealed in‐vitro thermodynamic properties in the presence and absence of NC. These discoveries were enabled by optical tweezers force spectroscopy, the nature of which requires a complex TAR hairpin construct. Unfortunately, design and synthesis of the necessary TAR hairpin construct has limited the quantity and diversity of thermodynamic analyses regarding this essential facet of HIV‐1. We aim to develop efficient protocols to produce TAR RNA hairpin constructs optimized for optical tweezer manipulations.In order to manipulate individual hairpin molecules, the optical tweezers require that the RNA hairpin be ligated to two long double‐stranded DNA “handles.” One handle of the final construct specifically binds a bead in the optical tweezer dual‐laser trap, while the other handle binds a second bead that is manipulated via a micropipette. Tiny displacements of the second bead in the laser trap are used to calculate the force required to unzip the hairpin, and the thermodynamic parameters for force‐induced unfolding. The two ~3,000 base pair handles were constructed by amplification of bacterial plasmid DNA by PCR, restricted to create specific sticky ends, and purified by gel electrophoresis. RNA hairpins were transcribed from short DNA templates that were themselves synthesized by automated solid‐phase synthesis with specific sequences for RNA pol binding and the hairpin structure of interest. The handles were then individually ligated via several short DNA and RNA linker sequences to the RNA hairpin. Gel purification between ligation steps removed alternative ligation products, allowing size‐based verification and isolation of intermediates and the complete hairpin construct. Initial unzipping data support the correct length and composition of the final construct. We anticipate that subsequent data will contribute to a broader understanding of RNA hairpin stability as well as provide useful insight about a potential HIV‐1 drug target.Support or Funding InformationNIH‐NIGMS (Grant 2R01‐GM072462‐10), NSF‐MCB (Grant 1213883), Camille and Henry Dreyfus FoundationThis abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

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.

Effect of loop composition on the stability and folding kinetics of RNA hairpins with large loops.

RNA hairpins are ubiquitous structural elements in biological RNAs, where they have the potential to regulate RNA folding and interactions with other molecules. There are established methods for predicting the thermodynamic stability of an RNA hairpin, but there are still relatively few detailed examinations of the kinetics of folding. Nonetheless, several recent studies indicate that hairpin folding does not proceed via a simple two-state model. Here, we monitor fluorescence from hairpins constructed as molecular beacons in ensemble, fluorescence correlation spectroscopy, and stopped-flow experiments to describe the folding of RNA hairpins with long (15 nucleotide) loops. Our results show that folding of these hairpins occurs through more than two states and that the mechanism of folding includes a fast intermediate phase observed on the tens of microseconds time scale and a slow phase, attributed to formation of the native folded hairpin loop and stem, observed on the milliseconds time scale. The composition of the RNA loop determines the time scale of intermediate and native folded states. Hairpins with a polyuracil loop sequence exhibit slower relaxation of the intermediate state and faster relaxation of the native folded state when compared to that of hairpins with cytosine or adenine in the loop. We hypothesize this composition dependence could be attributed to nucleobase stacking in cytosine and adenine containing regions of the loop, which would be absent in hairpins containing polyuracil loops. Such base stacking could destabilize the intermediate folds, thereby speeding the relaxation of the intermediate relative to similar sized hairpins with no base stacking in the loop. Likewise, the lower intermediate stability could prolong the relaxation of the native folded state.

Effects of Trimethylamine N-Oxide (TMAO) and Crowding Agents on the Stability of RNA Hairpins

We study the effect of the osmolyte, Trimethylamine N-Oxide (TMAO), which accumulates in cells in response to osmotic stress, on the stability of RNA hairpins. All atom molecular dynamics (MD) simulations of a nucleotide and the 22-nucleotide RNA hairpin P5GA in an aqueous TMAO solution show that TMAO preferentially interacts with the base through the formation of a single hydrogen bond. To circumvent the difficulties of adequately sampling the conformational space of polynucleotides, we used coarse-grained models (including one that is inspired by the results of all-atom MD simulations of a single nucleotide) to probe the effects of osmoyltes on the stability of P5GA. If, as revealed by our MD simulations, the cosolute specifically interacts with only one base at a time, then we find practically no change in hairpin stability as measured by Delta T m = T m(Phi) - T m, where T m(Phi) and T m are the melting temperatures at volume fraction Phi of the osmolyte and Phi = 0, respectively. This finding is in qualitative agreement with recent experiments. If the interactions between the RNA and osmolytes are repulsive, which is appropriate for mimicking the effects of crowding, Delta T m can vary from 5 to 15 K depending on the size of the osmolyte and the nature of RNA-osmolyte interactions. Cosolutes that interact favorably with multiple bases simultaneously can stabilize the hairpin more than a crowding agent of the same size. The implications of our predictions for experiments are briefly outlined.

Measurement of the Effect of Monovalent Cations on RNA Hairpin Stability

Using optical tweezers, we have measured the effect of monovalent cation concentration and species on the folding free energy of five large (49-124 nt) RNA hairpins, including HIV-1 TAR and molecules approximating A.U and G.C homopolymers. RNA secondary structure thermodynamics are accurately described by a model consisting of nearest-neighbor interactions and additive loop and bulge terms. Melting of small (<15 bp) duplexes and hairpins in 1 M NaCl has been used to determine the parameters of this model, which is now used extensively to predict structure and folding dynamics. Few systematic measurements have been made in other ionic conditions or for larger structures. By applying mechanical force, we measured the work required to fold and unfold single hairpins at room temperature over a range of cation concentrations from 50 to 1000 mM. Free energies were then determined using the Crooks fluctuation theorem. We observed the following: (1) In most cases, the nearest-neighbor model accurately predicted the free energy of folding at 1 M NaCl. (2) Free energy was proportional to the logarithm of salt concentration. (3) Substituting potassium ions for sodium slightly decreased hairpin stability. The TAR hairpin also misfolded nearly twice as often in KCl, indicating a differential kinetic response. (4) Monovalent cation concentration affects RNA stability in a sequence-dependent manner. G.C helices were unaffected by changing salt concentration, A.U helices were modestly affected, and the hairpin loop was very sensitive. Surprisingly, the U.C.U bulge of TAR was found to be equally stable in all conditions tested. We also report a new estimate for the elastic parameters of single-stranded RNA.

Stability of RNA hairpins closed by wobble base pairs.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXANmAYCC, where XY is the wobble base pair, GU or UG, and the underlined loop sequences are three to eight nucleotides. A nearest-neighbor analysis indicates the free energy of loop formation is dependent upon loop size and closing base pair. Hairpin loops closed by UG base pairs are on average 1.3 kcal/mol less stable than hairpins closed by GU base pairs. The hairpin loops closed by UG have approximately the same stability as hairpin loops closed by AU/UA base pairs, while the loops closed by GU are approximately 0.7 kcal/mol more stable than hairpins loops closed by GC/CG base pairs. These results, combined with the model previously developed [Serra et al. (1997) Biochemistry 36, 4844] to predict the stability for hairpin loops closed by Watson-Crick base pairs, allow for the following model to predict the stability of hairpin loops: delta G degree 37L(n) = delta G degree 37iL(n) + delta G degree 37mm + 0.6 (if closed by AU, UA, or UB) - 0.7 (if closed by GU) - 0.7 (if first mismatch is GA or UU except for loops closed by GU). Here, delta G degree 37iL(n) is the free energy increment for initiating a loop of n nucleotides with a CG or GC pair, and delta G degree 37mm is the free energy for the interaction of the first mismatch with the closing base pair. For hairpin loops of n = 4-9, delta G037iL(n) is 4.9, 5.0, 5.0, 5.0, 4.9, and 5.5 kcal/mol, respectively. For hairpin loops of n = 3, delta G degree 37L(3) = +4.8 + 0.6 (if closed by AU, UA, or UG) kcal/mol. Thermodynamic parameters for hairpin formation in 1 M NaCl for 13 naturally occurring RNA hairpin sequences closed by wobble base pairs are reported. The model provides good agreement for both TM and delta G degree 37 for most hairpins studied. Thermodynamic values for five terminal mismatches adjacent to wobble base pairs are also reported.

Improved parameters for the prediction of RNA hairpin stability.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXANmAYCC, where XY is the set of four Watson-Crick base pairs and the underlined loop sequences are three to nine nucleotides. A nearest neighbor analysis of the data indicates the free energy of loop formation at 37 degrees C is dependent upon loop size and closing base pair. The model previously developed to predict the stability for RNA hairpin loops (n > 3) includes contributions from the size of the loop, the identity of the closing base pair, the free energy increment (deltaGo(37mm)) for the interaction of the closing base pair with the first mismatch and an additional stabilization term for GA and UU first mismatches [Serra, M. J., Axenson, T. J., & Turner, D. H. (1994) Biochemistry 33, 14289]. The results presented here allow improvements in the parameters used to predict RNA hairpin stability. For hairpin loops of n = 4-9, deltaGo(37iL)(n) is 4.9, 5.0, 5.0, 5.0, 4.9, and 5.5 kcal/mol, respectively, and the penalty for hairpin closure by AU or UA is +0.6 kcal/mol. deltaGo(37iL)(n) is the free energy for initiating a loop of n nucleotides. The model for predicting hairpin loop stability for loops larger than three becomes deltaGo(37L)(n) = deltaGo(37iL)(n) + deltaGo(37mm) + 0.6(if closed by AU or UA) - 0.7(if first mismatch is GA or UU). Hairpin loops of three are modeled as independent of loop sequence with deltaGo(37iL)(3) = 4.8 and the penalty for AU closure of +0.6 kcal/mol. Thermodynamic parameters for hairpin formation in 1 M NaCl for 11 naturally occurring RNA hairpin sequences are reported. The model provides good agreement with the measured values for both T(M) (within 10 degrees C of the measured value) and deltaGo(37) (within 0.8 kcal/mol of the measured value) for hairpin formation. In general, the nearest neighbor model allows prediction of RNA hairpin stability to within 5-10% of the experimentally measured values.

Structural Features that Give Rise to the Unusual Stability of RNA Hairpins Containing GNRA Loops

The most frequently occurring RNA hairpins in 16S and 23S ribosomal RNA contain a tetranucleotide loop that has a GNRA consensus sequence. The solution structures of the GCAA and GAAA hairpins have been determined by nuclear magnetic resonance spectroscopy. Both loops contain an unusual G-A base pair between the first and last residue in the loop, a hydrogen bond between a G base and a phosphate, extensive base stacking, and a hydrogen bond between a sugar 2'-end OH and a base. These interactions explain the high stability of these hairpins and the sequence requirements for the variant and invariant nucleotides in the GNRA tetranucleotide loop family.

Prediction of rho-independent Escherichia coli transcription terminators : A statistical analysis of their RNA stem-loop structures

Escherichia coli rho-independent transcription terminators are characterized by an RNA structure having a G+C-rich stem-loop followed by a series of uridine residues, but they can be only partially predicted by the stability of this structure or by its primary sequence. A large number of such terminators have been identified or proposed in the literature, and we have constituted a list of them (148 found in 1021 x 10(3) base-pairs of E. coli DNA sequences) in order to analyze statistically the corresponding RNA hairpins. We show that the size of the loops presents a narrow distribution, that their sequences are not random, and that most loops are closed by a C.G base-pair. In particular, 55% of the loops are tetranucleotides and the most abundant loop sequences are UUCG and GAAA. These loops are abundant in prokaryotic and eukaryotic RNAs, and are known to enhance the stability of RNA hairpins. We propose that these tetraloops play an important role in the nucleation of the nascent RNA structures, as does also the presence of a C.G base-pair closing a hairpin loop. This analysis allows us to propose a model of formation of an RNA hairpin during the termination process and to construct an algorithm of prediction of the terminators in a given DNA sequence. For the E. coli sequences, it clearly distinguishes inter- from intracistronic terminator-like structures, and selects 141 of the 148 rho-independent terminators given in the literature, with a very low background. It also predicts with reasonable accuracy the in vitro termination efficiency of known rho-independent terminators, as well as predicting the existence of 35 as yet uncharacterized terminators.

Thermal stability of RNA hairpins containing a four-membered loop and a bulge nucleotide

Fourteen RNA hairpins containing a four-membered loop and a bulge nucleotide were synthesized and their thermal stabilities determined. The combined contribution of a four-membered loop and bulge A to the free energy of a hairpin is calculated to be 9.3 kcal/mol at 37 degrees C and successfully predicts the stability of an independent RNA hairpin. The introduction of a bulge nucleotide to the helical stem of an RNA hairpin destabilizes the molecule in a sequence-dependent manner. The individual thermodynamic contributions of a four-membered loop and bulge A, G, and U residues to the stability of an RNA hairpin loop are presented.