Search

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5'GGUC3' paired to 3'CUGG5' separately from other 5'GU3'/3'UG5' stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The influenza A virus (IAV) is a continuous health threat to humans as well as animals due to its recurring epidemics and pandemics. The IAV genome is segmented and the eight negative-sense viral RNAs (vRNAs) are transcribed into positive sense complementary RNAs (cRNAs) and viral messenger RNAs (mRNAs) inside infected host cells. A role for the secondary structure of IAV mRNAs has been hypothesized and debated for many years, but knowledge on the structure mRNAs adopt in vivo is currently missing. Here we solve, for the first time, the in vivo secondary structure of IAV mRNAs in living infected cells. We demonstrate that, compared to the in vitro refolded structure, in vivo IAV mRNAs are less structured but exhibit specific locally stable elements. Moreover, we show that the targeted disruption of these high-confidence structured domains results in an extraordinary attenuation of IAV replicative capacity. Collectively, our data provide the first comprehensive map of the in vivo structural landscape of IAV mRNAs, hence providing the means for the development of new RNA-targeted antivirals., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA folding free energy change nearest neighbor parameters are widely used to predict folding stabilities of secondary structures. They were determined by linear regression to datasets of optical melting experiments on small model systems. Traditionally, the optical melting experiments are analyzed assuming a two-state model, i.e. a structure is either complete or denatured. Experimental evidence, however, shows that structures exist in an ensemble of conformations. Partition functions calculated with existing nearest neighbor parameters predict that secondary structures can be partially denatured, which also directly conflicts with the two-state model. Here, a new approach for determining RNA nearest neighbor parameters is presented. Available optical melting data for 34 Watson-Crick helices were fit directly to a partition function model that allows an ensemble of conformations. Fitting parameters were the enthalpy and entropy changes for helix initiation, terminal AU pairs, stacks of Watson-Crick pairs and disordered internal loops. The resulting set of nearest neighbor parameters shows a 38.5% improvement in the sum of residuals in fitting the experimental melting curves compared to the current literature set.

We have solved at 1. 07: Å resolution the X-ray crystal structure of a polyriboadenylic acid (poly(rA)) parallel and continuous double helix. Fifty-nine years ago, double helices of poly(rA) were first proposed to form at acidic pH. Here, we show that 7-mer oligo(rA), i.e. rA7, hybridizes and overlaps in all registers at pH 3.5 to form stacked double helices that span the crystal. Under these conditions, rA7 forms well-ordered crystals, whereas rA6 forms fragile crystalline-like structures, and rA5, rA8 and rA11 fail to crystallize. Our findings support studies from ∼50 years ago: one showed using spectroscopic methods that duplex formation at pH 4.5 largely starts with rA7 and begins to plateau with rA8; another proposed a so-called 'staggered zipper' model in which oligo(rA) strands overlap in multiple registers to extend the helical duplex. While never shown, protonation of adenines at position N1 has been hypothesized to be critical for helix formation. Bond angles in our structure suggest that N1 is protonated on the adenines of every other rAMP-rAMP helix base pair. Our data offer new insights into poly(rA) duplex formation that may be useful in developing a pH sensor., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Oligonucleotide microarrays are widely used in various biological studies. In this review, application of oligonucleotide microarrays for identifying binding sites and probing structure of RNAs is described. Deep sequencing allows fast determination of DNA and RNA sequence. High-throughput methods for determination of secondary structures of RNAs have also been developed. Those methods, however, do not reveal binding sites for oligonucleotides. In contrast, microarrays directly determine binding sites while also providing structural insights. Microarray mapping can be used over a wide range of experimental conditions, including temperature, pH, various cations at different concentrations and the presence of other molecules. Moreover, it is possible to make universal microarrays suitable for investigations of many different RNAs, and readout of results is rapid. Thus, microarrays are used to provide insight into oligonucleotide sequences potentially able to interfere with biological function. Better understanding of structure-function relationships of RNA can be facilitated by using microarrays to find RNA regions capable to bind oligonucleotides. That information is extremely important to design optimal sequences for antisense oligonucleotides and siRNA because both bind to single-stranded regions of target RNAs., (© The Author(s) 2014. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Thermodynamic data are reported revealing that pseudouridine (Ψ) can stabilize RNA duplexes when replacing U and forming Ψ-A, Ψ-G, Ψ-U and Ψ-C pairs. Stabilization is dependent on type of base pair, position of Ψ within the RNA duplex, and type and orientation of adjacent Watson-Crick pairs. NMR spectra demonstrate that for internal Ψ-A, Ψ-G and Ψ-U pairs, the N3 imino proton is hydrogen bonded to the opposite strand nucleotide and the N1 imino proton may also be hydrogen bonded. CD spectra show that general A-helix structure is preserved, but there is some shifting of peaks and changing of intensities. Ψ has two hydrogen donors (N1 and N3 imino protons) and two hydrogen bond acceptors because the glycosidic bond is C-C rather than C-N as in uridine. This greater structural potential may allow Ψ to behave as a kind of structurally driven universal base because it can enhance stability relative to U when paired with A, G, U or C inside a double helix. These structural and thermodynamic properties may contribute to the biological functions of Ψ.

The Nearest Neighbor Database (NNDB, http://rna.urmc.rochester.edu/NNDB) is a web-based resource for disseminating parameter sets for predicting nucleic acid secondary structure stabilities. For each set of parameters, the database includes the set of rules with descriptive text, sequence-dependent parameters in plain text and html, literature references to experiments and usage tutorials. The initial release covers parameters for predicting RNA folding free energy and enthalpy changes.

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2' oxygen and 4' carbon of ribose) and 2,6-diaminopurine were incorporated into 2'-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5' coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5' RNA). This R2Bm 5' RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5' RNA.

Modified nucleotides are useful tools to study the structures, biological functions and chemical and thermodynamic stabilities of nucleic acids. Derivatives of 2,6-diaminopurine riboside (D) are one type of modified nucleotide. The presence of an additional amino group at position 2 relative to adenine results in formation of a third hydrogen bond when interacting with uridine. New method for chemical synthesis of protected 3'-O-phosphoramidite of LNA-2,6-diaminopurine riboside is described. The derivatives of 2'-O-methyl-2,6-diaminopurine and LNA-2,6-diaminopurine ribosides were used to prepare complete 2'-O-methyl RNA and LNA-2'-O-methyl RNA chimeric oligonucleotides to pair with RNA oligonucleotides. Thermodynamic stabilities of these duplexes demonstrated that replacement of a single internal 2'-O-methyladenosine with 2'-O-methyl-2,6-diaminopurine riboside (D(M)) or LNA-2,6-diaminopurine riboside (D(L)) increases the thermodynamic stability (DeltaDeltaG degrees 37) on average by 0.9 and 2.3 kcal/mol, respectively. Moreover, the results fit a nearest neighbor model for predicting duplex stability at 37 degrees C. D-A and D-G but not D-C mismatches formed by D(M) or D(L) generally destabilize 2'-O-methyl RNA/RNA and LNA-2'-O-methyl RNA/RNA duplexes relative to the same type of mismatches formed by 2'-O-methyladenosine and LNA-adenosine, respectively. The enhanced thermodynamic stability of fully complementary duplexes and decreased thermodynamic stability of some mismatched duplexes are useful for many RNA studies, including those involving microarrays.

Results from optical melting studies of Watson-Crick complementary heteroduplexes formed between 2'-O-methyl RNA and RNA oligonucleotides are used to determine nearest neighbor thermodynamic parameters for predicting the stabilities of such duplexes. The results are consistent with the physical model assumed by the individual nearest neighbor-hydrogen bonding model, which contains terms for helix initiation, base pair stacking and base pair composition. The sequence dependence is similar to that for Watson-Crick complementary RNA/RNA duplexes, which suggests that the sequence dependence may also be similar to that for other backbones that favor A-form RNA conformations.

A complete set of nearest neighbor parameters to predict the enthalpy change of RNA secondary structure formation was derived. These parameters can be used with available free energy nearest neighbor parameters to extend the secondary structure prediction of RNA sequences to temperatures other than 37 degrees C. The parameters were tested by predicting the secondary structures of sequences with known secondary structure that are from organisms with known optimal growth temperatures. Compared with the previous set of enthalpy nearest neighbor parameters, the sensitivity of base pair prediction improved from 65.2 to 68.9% at optimal growth temperatures ranging from 10 to 60 degrees C. Base pair probabilities were predicted with a partition function and the positive predictive value of structure prediction is 90.4% when considering the base pairs in the lowest free energy structure with pairing probability of 0.99 or above. Moreover, a strong correlation is found between the predicted melting temperatures of RNA sequences and the optimal growth temperatures of the host organism. This indicates that organisms that live at higher temperatures have evolved RNA sequences with higher melting temperatures.

The influence of locked nucleic acid (LNA) residues on the thermodynamic properties of 2'-O-methyl RNA/RNA heteroduplexes is reported. Optical melting studies indicate that LNA incorporated into an otherwise 2'-O-methyl RNA oligonucleotide usually, but not always, enhances the stabilities of complementary duplexes formed with RNA. Several trends are apparent, including: (i) a 3' terminal U LNA and 5' terminal LNAs are less stabilizing than interior and other 3' terminal LNAs; (ii) most of the stability enhancement is achieved when LNA nucleotides are separated by at least one 2'-O-methyl nucleotide; and (iii) the effects of LNA substitutions are approximately additive when the LNA nucleotides are separated by at least one 2'-O-methyl nucleotide. An equation is proposed to approximate the stabilities of complementary duplexes formed with RNA when at least one 2'-O-methyl nucleotide separates LNA nucleotides. The sequence dependence of 2'-O-methyl RNA/RNA duplexes appears to be similar to that of RNA/RNA duplexes, and preliminary nearest-neighbor free energy increments at 37 degrees C are presented for 2'-O-methyl RNA/RNA duplexes. Internal mismatches with LNA nucleotides significantly destabilize duplexes with RNA.

The X-ray crystallographic structure of the RNA duplex [r(CGCAIGCG)]2 has been refined to 2.5 A. It shows a symmetric internal loop of two non-Watson-Crick base pairs which form in the middle of the duplex. The tandem A-I/I-A pairs are related by a crystallographic two-fold axis. Both A(anti)-I(anti) mismatches are in a head-to-head conformation forming hydrogen bonds using the Watson-Crick positions. The octamer duplexes stack above one another in the cell forming a pseudo-infinite helix throughout the crystal. A hydrated calcium ion bridges between the 3'-terminal of one molecule and the backbone of another. The tandem A-I mismatches are incorporated with only minor distortion to the backbone. This is in contrast to the large helical perturbations often produced by sheared G-A pairs in RNA oligonucleotides.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXAUAAUAYCC, where X and Y are CG, GC, AU, UA, GU, or UG. A nearest neighbor analysis of the data indicates the free energy change for loop formation at 37 degrees C, delta degrees Gl,37, averages 3.4 kcal/mol for hairpin loops closed with C.G, G.C, and G.U pairs. In contrast, delta G degree l,37 averages 4.6 kcal/mol for loops closed with A.U, U.A, or U.G pairs. Thus the stability of an RNA hairpin depends on the closing base pair. The hairpin with a GA mismatch that is formed by GGCGUAAUAGCC is more stable than the corresponding hairpin with an AA mismatch. Thus hairpin stability also depends on loop sequence. These effects are not included in current algorithms for prediction of RNA structure from sequence.

An effect of 5'-phosphorylation on the stability of triple helical DNA containing pyrimidine:purine:pyrimidine strands has been demonstrated by both gel electrophoresis and UV melting. A 5'-phosphate on the purine-rich middle strand of a triple helix lowers the stability of triple helix formation by approximately 1 kcal/mol at 25 degrees C. The middle strand is involved in both Watson-Crick and Hoogsteen base pairing. In contrast, a 5'-phosphate on the pyrimidine-rich strands, which are involved in either Watson-Crick or Hoogsteen base pairing, has a smaller effect on the stability of triple helix. The order of stability is: no phosphate on either strand > phosphate on both pyrimidine strands > phosphate on purine strand > phosphate on all three strands. Differential stability of triple helix species is postulated to stem from an increase in rigidity due to steric hindrance from the 5'-phosphate. This result indicates that labelling with 32P affect equilibrium in triplex formation.

Oligoribonucleotides with 2'-5' linkages have been synthesized on solid support. UV melting and CD experiments indicate complementary strands associate to give complexes with melting temperatures 30 to 40 degrees C lower than for duplexes formed by 3'-5' oligoribonucleotides with the same sequence. UV melting and imino proton NMR spectra and NOEs for (2'-5') CGGCGCCG are consistent with formation of an antiparallel duplex. The results suggest greater duplex stability was one factor favoring 3'-5' over 2'-5' linkages in evolution.

This article describes the latest version of an RNA folding algorithm that predicts both optimal and suboptimal solutions based on free energy minimization. A number of RNA's with known structures deduced from comparative sequence analysis are folded to test program performance. The group of solutions obtained for each molecule is analysed to determine how many of the known helixes occur in the optimal solution and in the best suboptimal solution. In most cases, a structure about 80% correct is found with a free energy within 2% of the predicted lowest free energy structure.

The self-splicing intervening sequence from the precursor rRNA of Tetrahymena thermophila cyclizes to form a covalently closed circle. This circle can be reopened by reaction with oligonucleotides or water. The kinetics of circle opening as a function of substrate and Mg2+ concentrations have been measured for dCrU, rCdU, dCdT, and H2O addition. Comparisons with previous results for rCrU suggest: (1) the 2' OH of the 5' sugar of a dinucleoside phosphate is involved in substrate binding, and (2) the 2' OH of the 3' sugar of a dimer substrate is involved in Mg2+ binding. Evidently, the binding site for a required Mg2+ ion is dependent on both the ribozyme and the dimer substrate. The apparent activation energy and entropy for circle opening by hydrolysis are 31 kcal/mol and 50 eu, respectively. The large, positive activation entropy suggests a partial unfolding of the ribozyme is required for reaction.

The Mg2+ contribution to the reaction of circular intervening sequence (CIVS) from rRNA precursor of Tetrahymena thermophila with a dinucleotide CU has been investigated. The results indicated that the circle opening of CIVS may involve binding of a weakly held Mg2+ ion.

The first photochemical crosslinking of a protein to a nucleic acid using laser excitation is reported. A single, 120 mJ, 20 ns pulse at 248 nm crosslinks about 10% of bound E. coli RNA polymerase to T7 DNA under the conditions studied. The crosslinking yield depends on mercaptoethanol concentration, and is a linear function of laser intensity. The protein subunits crosslinked to DNA are beta, beta' and sigma.

The kinetics of U-U-C-A binding to the dodecanucleotide (A-Cm-U-Gm-A-A-Y-A-psi-m5C-U-Gp) isolated from the anticodon region of yeast tRNA-Phe are similar to the kinetics of binding of U-U-C-A to intact tRNA-Phe. A large enhancement in binding constant over that predicted for U-U-C-A-U-G-A-A is observed for both the complexes of dodecanucleotide and tRNA-Phe with U-U-C-A. This strongly suggests that both the anticodon loop in tRNA-Phe and the dodecanucleotide can form four base pairs with U-U-C-A. Furthermore, the enhanced stability cannot be attributed to a special conformation of the anticodon loop, but instead the anticodon loop is probably flexible. A likely explanation for the increased binding is the effect of non-base-paired ends. This increased thermodynamic stability comes from a larger entropy gain rather than a larger enthalpy decrease.