No abstract

Chemically modified antisense oligonucleotides (ASOs) are widely used as a tool to functionalize microRNAs (miRNAs). Reduction of miRNA level after ASO inhibition is commonly reported to show efficacy. Whether this is the most relevant endpoint for measuring miRNA inhibition has not been adequately addressed in the field although it has important implications for evaluating miRNA targeting studies. Using a novel approach to quantitate miRNA levels in the presence of excess ASO, we have discovered that the outcome of miRNA inhibition can vary depending on the chemical modification of the ASO. Although some miRNA inhibitors cause a decrease in mature miRNA levels, we have identified a novel 2′-fluoro/2′-methoxyethyl modified ASO motif with dramatically improved in vivo potency which does not. These studies show there are multiple mechanisms of miRNA inhibition by ASOs and that evaluation of secondary endpoints is crucial for interpreting miRNA inhibition studies.

Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5′GG/3′UU motif were added to the single representation in the previous database. This revises the ΔG°37 for the 5′GG/3′UU motif from an unfavorable 0.5 kcal/mol to a favorable −0.2 kcal/mol. Similarly, the ΔG°37 for the 5′UG/3′GU motif changes from 0.3 to −0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°37, ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.

Optical melting was used to determine the stabilities of 11 small RNA oligomers of defined secondary structure as a function of magnesium ion concentration. The oligomers included helices composed of Watson-Crick base pairs, GA tandem base pairs, GU tandem base pairs, and loop E motifs (both eubacterial and eukaryotic). The effect of magnesium ion concentration on stability was interpreted in terms of two simple models. The first assumes an uptake of metal ion upon duplex formation. The second assumes nonspecific electrostatic attraction of metal ions to the RNA oligomer. For all oligomers, except the eubacterial loop E, the data could best be interpreted as nonspecific binding of metal ions to the RNAs. The effect of magnesium ions on the stability of the eubacterial loop E was distinct from that seen with the other oligomers in two ways. First, the extent of stabilization by magnesium ions (as measured by either change in melting temperature or free energy) was three times greater than that observed for the other helical oligomers. Second, the presence of magnesium ions produces a doubling of the enthalpy for the melting transition. These results indicate that magnesium ion stabilizes the eubacterial loop E sequence by chelating the RNA specifically. Further, these results on a rather small system shed light on the large enthalpy changes observed upon thermal unfolding of large RNAs like group I introns. It is suggested that parts of those large enthalpy changes observed in the folding of RNAs may be assigned to variations in the hydration states and types of coordinating atoms in some specifically bound magnesium ions and to an increase in the observed cooperativity of the folding transition due to the binding of those magnesium ions coupling the two stems together. Brownian dynamic simulations, carried out to visualize the metal ion binding sites, reveal rather delocalized ionic densities in all oligomers, except for the eubacterial loop E, in which precisely located ion densities were previously calculated.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGCXUAAUYGCC, where XY is the set of 10 possible mismatch base pairs. A nearest neighbor analysis of the data indicates the free energy for loop formation at 37 C varies from 2.9 to 4.5 kcal/mol. Thermodynamic parameters are also reported for hairpin formation by RNA sequences of the type GGXGUAAUAYCC (where XY are CG, GC, AU, UA, GU, and UG), with the common naturally occurring GA first mismatch (45% of small and large subunit rRNA loops of six). These results allow the development of a model to predict the stability of RNA hairpin loops. The model includes the size of the loop, the identity of the closing base pair, the free energy increment (delta G zero 37MM) for interaction of the closing base pair with the first mismatch, and an additional stabilization term for GA and UU first mismatches. delta G zero 37L(n) = delta G zero 37i(n) + delta G zero 37MM + 0.4 (if closed by AU or UA) -0.7 (if first mismatch is GA or UU). Here delta G zero 37i(n) is the free energy for initiating a loop of n nucleotides. delta G zero 37i(n) for n = 4-9 is 4.9, 4.4, 5.0, 5.0, 5.1, and 5.2 kcal/mol, respectively. The delta G zero 37MM is derived from measurements of model duplexes with terminal mismatches. The model gives good agreement when tested against four naturally occurring hairpin sequences.

Thermodynamic parameters are reported for duplex formation of 48 self-complementary RNA duplexes containing Watson–Crick terminal base pairs (GC, AU and UA) with all 16 possible 3′ double-nucleotide overhangs; mimicking the structures of short interfering RNAs (siRNA) and microRNAs (miRNA). Based on nearest-neighbor analysis, the addition of a second dangling nucleotide to a single 3′ dangling nucleotide increases stability of duplex formation up to 0.8 kcal/mol in a sequence dependent manner. Results from this study in conjunction with data from a previous study [A. S. O'Toole, S. Miller and M. J. Serra (2005) RNA, 11, 512.] allows for the development of a refined nearest-neighbor model to predict the influence of 3′ double-nucleotide overhangs on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against five oligomers with various core duplex sequences. Phylogenetic analysis of naturally occurring miRNAs was performed to support our results. Selection of the effector miR strand of the mature miRNA duplex appears to be dependent upon the identity of the 3′ double-nucleotide overhang. Thermodynamic parameters for 3′ single terminal overhangs adjacent to a UA pair are also presented.

The rates of duplex formation for two octamers of DNA (5' d-CACGGCTC/5' d-GAGCCGTG and 5' d-CACAGCAC/5' d-GTGCTGTG), the homologous RNA, and both sets of hybrids in 1 M NaCl buffer have been measured using stopped-flow spectroscopy. In addition, the thermodynamic parameters, ΔH° and ΔS°, have been determined for the same sequences under the same buffer conditions using optical melting techniques. These data reveal a linear free energy relationship between the free energy of activation for denaturation and the change in free energy for formation of the duplexes. This relationship indicates that these duplex formation reactions occur through a common unstructured transition state that is more similar to the single strands in solution than to the ensuing duplex. In addition, these data confirm that the greater stability of RNA duplexes relative to that of homologous DNA and hybrid duplexes is controlled by the denaturation rate and not the duplex formation rate.

ABSTRACT