Antisense RNA control of gene expression in bacteriophage P22. II. Kinetic mechanism and cation specificity of the pairing reaction.

The bacteriophage P22 sar RNA-ant mRNA pairing reaction was characterized kinetically. The pairing reaction proceeds by a three-step pathway. First, reversible base pairs form between complementary hairpin loops in sar RNA and ant RNA (Kd = 270 nM). Next, stable duplex formation initiates between single-stranded nucleotides in sar RNA and ant RNA; the ant RNA nucleotides are at the bottom of a hairpin stem that is partially accessible. Concomitant unwinding of one sar RNA hairpin and the complementary ant RNA hairpin then occurs, to form a partially paired intermediate (k2 = 12 min(-1). Finally, a complete duplex forms after unwinding of the other sar RNA hairpin and the complementary ant RNA hairpin (k3 = 7 min(-1). Experiments with sar RNA sequence and length variants demonstrate that the precise structures of both sar RNA hairpins affect the kinetic parameters. The pairing reaction is Mg2+-dependent, and shows high specificity for the required cation. Maximal pairing rates are achieved when more than one Mg2+ ion is bound. The cation-dependence and specificity indicate a requirement for Mg2+-dependent tertiary structure.     

Perturbation of DNA hairpins containing the EcoRI recognition site by hairpin loops of varying size and composition: physical (NMR and UV) and enzymatic (EcoRI) studies.

We have investigated loop-induced structural perturbation of the stem structure in hairpins d(GAATTCXnGAATTC) (X = A, T and n = 3, 4, 5 and 6) that contain an EcoRI restriction site in close proximity to the hairpin loop. Oligonucleotides containing either a T3 or a A3 loop were not hydrolyzed by the restriction enzyme and also showed only weak binding to EcoRI in the absence of the cofactor Mg2+. In contrast, hairpins with larger loops are hydrolyzed by the enzyme at the scission site next to the loop although the substrate with a A4 loop is significantly more resistant than the oligonucleotide containing a T4 loop. The hairpin structures with 3 loop residues were found to be thermally most stable while larger hairpin loops resulted in structures with lower melting temperatures. The T-loop hairpins are thermally more stable than the hairpins containing the same number of A residues in the loop. As judged from proton NMR spectroscopy and the thermodynamic data, the base pair closest to the hairpin loop did form in all cases studied. The hairpin loops did, however, affect the conformation of the stem structure of the hairpins. From 31P and 1H NMR spectroscopy we conclude that the perturbation of the stem structure is stronger for smaller hairpin loops and that the extent of the perturbation is limited to 2-3 base pairs for hairpins with T3 or A4 loops. Our results demonstrate that hairpin loops modulate the conformation of the stem residues close to the loop and that this in turn reduces the substrate activity for DNA sequence specific proteins.     

Syntheses of alternating oligo-2'-O-methylribonucleoside methylphosphonates and their interactions with HIV TAR RNA

Oligonucleotide analogues 15-20 nucleotides in length have been prepared, whose sequences are complementary to nucleotides in the upper hairpin of HIV TAR RNA. These alternating oligonucleoside methylphosphonates, mr-AOMPs, contain 2'-O-methylribonucleosides and alternating methylphosphonate and phosphodiester internucleotide linkages. The methylphosphonate and phosphodiester linkages of these oligomers are highly resistant to hydrolysis by exonuclease activity found in mammalian serum and to endonucleases, such as S1 nuclease. The oligomers were prepared using automated phosphoramidite chemistry and terminate with a 5'-phosphate group, which provides an affinity handle for purification by strong anion exchange HPLC. A 15-mer mr-AOMP, 1676, that is complementary to the 5'-side of the TAR RNA hairpin, including the 3-base bulge and 6-base loop region, forms a 1:1 duplex with a complementary RNA 18-mer, mini-TAR RNA. The T(m) of this duplex is 71 degrees C, which is similar to that of the duplex formed by the corresponding all phosphodiester 15-mer. Introduction of two mismatched bases reduces the T(m) by 17 degrees C. The apparent dissociation constant, K(d), for the 1676/mini-TAR RNA duplex as determined by an electrophoretic mobility shift assay at 37 degrees C is 0.3 nM. Oligomer 1676 also binds tightly to the full length TAR RNA target under physiological conditions (K(d) = 20 nM), whereas no binding was observed by the mismatched oligomer. A 19-mer that is complementary to the entire upper hairpin also binds to TAR RNA with a K(d) that is similar to that of 1676, a result that suggests only part of the oligomer binds. When two of the methylphosphonate linkages in the region complementary to the 6-base loop are replaced with phosphodiester linkages, the K(d) is reduced by approximately a factor of 10. This result suggests that interactions between TAR RNA and the oligomer occur initially with nucleotides in the 6-base loop, and that these interactions are sensitive to presence and possibly the chirality of the methylphosphonate linkages in the oligomer. The high affinities of mr-AOMPs for TAR RNA and their resistance to nuclease hydrolysis suggests their potential utility as antisense agents in cell culture.     

Molecular recognition of 6'-N-5-hexynoate kanamycin A and RNA 1x1 internal loops containing CA mismatches

In our previous study to identify the RNA internal loops that bind an aminoglycoside derivative, we determined that 6'-N-5-hexynoate kanamycin A prefers to bind 1x1 nucleotide internal loops containing C·A mismatches. In this present study, the molecular recognition between a variety of RNAs that are mutated around the C·A loop and the ligand was investigated. Studies show that both loop nucleotides and loop closing pairs affect binding affinity. Most interestingly, it was shown that there is a correlation between the thermodynamic stability of the C·A internal loops and ligand affinity. Specifically, C·A loops that had relatively high or low stability bound the ligand most weakly whereas loops with intermediate stability bound the ligand most tightly. In contrast, there is no correlation between the likelihood that a loop forms a C-A(+) pair at lower pH and ligand affinity. It was also found that a 1x1 nucleotide C·A loop that bound to the ligand with the highest affinity is identical to the consensus site in RNAs that are edited by adenosine deaminases acting on RNA type 2 (ADAR2). These studies provide a detailed investigation of factors affecting small molecule recognition of internal loops containing C·A mismatches, which are present in a variety of RNAs that cause disease.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

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.     

Two-dimensional combinatorial screening and the RNA Privileged Space Predictor program efficiently identify aminoglycoside–RNA hairpin loop interactions

Herein, we report the identification of RNA hairpin loops that bind derivatives of kanamycin A, tobramycin, neamine, and neomycin B via two-dimensional combinatorial screening, a method that screens chemical and RNA spaces simultaneously. An arrayed aminoglycoside library was probed for binding to a 6-nucleotide RNA hairpin loop library (4096 members). Members of the loop library that bound each aminoglycoside were excised from the array, amplified and sequenced. Sequences were analyzed with our newly developed RNA Privileged Space Predictor (RNA-PSP) program, which analyzes selected sequences to identify statistically significant trends. RNA-PSP identified the following unique trends: 5′UNNNC3′ loops for the kanamycin A derivative (where N is any nucleotide); 5′UNNC3′ loops for the tobramycin derivative; 5′UNC3′ loops for the neamine derivative; and 5′UNNG3′ loops for the neomycin B derivative. The affinities and selectivities of a subset of the ligand–hairpin loop interactions were determined. The selected interactions have K_d values ranging from 10 nM to 605 nM. Selectivities ranged from 0.4 to >200-fold. Interestingly, the results from RNA-PSP are able to qualitatively predict specificity based on overlap between the RNA sequences selected for the ligands. These studies expand the information available on small molecule–RNA motif interactions, which could be useful to design ligands targeting RNA.     

Direct mass spectrometric determination of the stoichiometry and binding affinity of the complexes between nucleocapsid protein and RNA stem-loop hairpins of the HIV-1 Psi-recognition element

The formation of noncovalent complexes between the HIV-1 nucleocapsid protein p7 (NC) and RNA hairpins SL2-SL4 of the Psi-recognition element was investigated by direct infusion electrospray ionization-Fourier transform mass spectrometry (ESI-FTMS). The high resolution afforded by this method provided the unambiguous characterization of the stoichiometry and composition of complexes formed by multiple equilibria in solution. For each hairpin, the formation of a 1:1 complex was found to be the primary binding mode in solutions of intermediate salt content (150 mM ammonium acetate). Binding of multiple units of NC was observed with lower affinity and a maximum stoichiometry matching the limit calculated from the number of nucleotides in the construct and the size of the footprint of NC onto single-stranded nucleic acids, thus implying the defolding of the hairpin three-dimensional (3D) structure. Dissociation constants of 62 +/- 22 nM, 178 +/- 64 nM, and 1.3 +/- 0.5 microM were determined for SL2, SL3-2, and SL4, respectively, which are similar to values obtained by spectroscopic and calorimetric methods with the additional confidence offered by a direct, rather than inferred, knowledge of the binding stoichiometry. Competitive binding experiments carried out in solutions of intermediate ionic strength, which has the effect of weakening the electrostatic interactions in solution, provided a direct way of evaluating the stabilizing contributions of H-bonding and hydrophobic interactions that are more sensitive to the sequence and structural context of the different hairpins. The relative scale of binding affinity obtained in this environment reflects the combination of contributions provided by the different structures of both the tetraloop and the double-stranded stem. The importance of the stem 3D structure in modulating the binding activity was tested by a competitive binding experiment that included the SL3-2 RNA construct, a DNA analogue of SL3 (SL3(DNA)), and a DNA analogue in which all four loop bases were replaced with abasic nucleotides (SL3(abasic)). NC was found to bind the A-type double-stranded stem of SL3-2 RNA at least 30 times more tightly than the B-type helical structure of SL3(DNA). Eliminating the stabilization provided by the interactions with the tetraloop bases made the binding of SL3(abasic) approximately 50 times weaker than that of SL3(DNA).     

An RNA aptamer to the xanthine/guanine base with a distinctive mode of purine recognition.

RNAs that bind to xanthine (2,6-dioxypurine) were isolated from a population of 10(12) random sequences by in vitro selection. These xanthine-binding RNAs were found to have a 10 nt consensus sequence at an internal loop in the most probable secondary structure. By trimming one of the xanthine-binding RNAs, a representative xanthine-binding RNA (designated as XBA) of 32 nt residues was prepared. The dissociation constant of this RNA for xanthine was determined to be 3.3 microM by equilibrium filtration experiments. The XBA RNA can bind to guanine as well, whereas it hardly accommodates adenine, cytosine or uracil. The K d values for various xanthine/guanine analogues were determined, and revealed that the N1H, N7 and O6 moieties of the ligand are involved in the binding with the XBA RNA. The ribonuclease sensitivities of some internal-loop residues changed upon the addition of xanthine, suggesting that the internal loop of the XBA RNA is involved in the ligand binding. Interestingly, the consensus sequence of the xanthine/guanine-binding RNAs is the same as a sequence in one of the internal loops of the hairpin ribozyme, except for a substitution that is neutral with respect to xanthine/guanine binding.     

Measuring dissociation constants of RNA and aminoglycoside antibiotics by electrospray ionization mass spectrometry

Electrospray ionization mass spectrometry (ESI-MS) has been used to determine the dissociation constants (K(D)s) and binding stoichiometry for tobramycin and paromomycin with a 27-nucleotide RNA construct representing the A-site of the 16S ribosomal RNA. K(D) values determined by holding the ligand concentration fixed are compared with K(D) values derived by holding the RNA target concentration fixed. Additionally, the effect of solution conditions such as the amount of organic solvent present and the amount of salt present in the solution on the K(D) measurement is investigated. It is shown that the preferred method for determining dissociation constants using ESI-MS is holding the RNA target concentration fixed below the expected K(D) and titrating the ligand. K(D) measurements should also be carried out at as high as possible salt concentration to minimize nonspecific binding due primarily to electrostatic interactions. For tobramycin, two nonequivalent binding sites were found with K(D1) = 352 nM and K(D2) = 9 microM. For paromomycin, there is only one binding site with K(D) = 52 nM.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Human UP1 as a model for understanding purine recognition in the family of proteins containing the RNA recognition motif (RRM)

Heterogeneous ribonucleoprotein A1 (hnRNP A1) is a prototype for the family of eukaryotic RNA processing proteins containing the common RNA recognition motif (RRM). The region consisting of residues 1-195 of hnRNP A1 is referred to as UP1. This region has two RRMs and has a high affinity for both single-stranded RNA and the human telomeric repeat sequence d(TTAGGG)(n). We have used UP1's novel DNA binding to investigate how RRMs bind nucleic acid bases through their highly conserved RNP consensus sequences. Nine complexes of UP1 bound to modified telomeric repeats were investigated using equilibrium fluorescence binding and X-ray crystallography. In two of the complexes, alteration of a guanine to either 2-aminopurine or nebularine resulted in an increase in K(d) from 88nM to 209nM and 316nM, respectively. The loss of these orienting interactions between UP1 and the substituted base allows it to flip between syn and anti conformations. Substitution of the same base with 7-deaza-guanine preserves the O6/N1 contacts but still increases the K(d) to 296nM and suggests that it is not simply the loss of affinity that gives rise to the base mobility, but also the stereochemistry of the specific contact to O6. Although these studies provide details of UP1 interactions to nucleic acids, three general observations about RRMs are also evident: (1) as suggested by informatic studies, main-chain to base hydrogen bonding makes up an important aspect of ligand recognition (2) steric clashes generated by modification of a hydrogen bond donor-acceptor pair to a donor-donor pair are poorly tolerated and (3) a conserved lysine position proximal to RNP-2 (K(106)-IFVGGI) orients the purine to allow stereochemical discrimination between adenine and guanine based on the 6-position. This single interaction is well-conserved in known RRM structures and appears to be a broad indicator for purine preference in the larger family of RRM proteins.     

A double hairpin structure is necessary for the efficient encapsidation of spleen necrosis virus retroviral RNA.

We conducted a mutational analysis within the previously defined encapsidation sequence (E) for spleen necrosis virus (SNV), an avian retrovirus. We found that two regions are necessary for efficient SNV replication. The first region is a double hairpin structure as proposed by Konings et al. (1992, J. Virol., 66, 632-640); the second region is located downstream of the hairpins. We showed further that the double hairpin structure is required for efficient SNV RNA encapsidation. Our work is the first to demonstrate, via linker-scanning and site-directed mutagenesis, that a specific RNA secondary structure is required for the encapsidation of retroviral RNA. Analysis of a series of mutations within the E region indicates (i) that preserving the secondary structure of the two hairpins is important for efficient encapsidation and (ii) that the stem regions of the hairpins contain specific sequences critical for encapsidation. Within the hairpins, the presence of at least one of the two conserved GACG four-residue loops, but not the moderately conserved bulge sequence of the first hairpin, is crucial for function. The function of the hairpins is independent of the relative order of the two hairpins. However, the two hairpins are not redundant and are not functionally identical. Replacement of SNV double hairpin sequences with those of Moloney murine leukemia virus (M-MLV) has no detectable effect on the replication of SNV-based retrovirus vectors with reticuloendotheliosis virus strain A (REV-A) helper virus. Furthermore, replacement of the entire E sequence of SNV with that of Moloney murine sarcoma virus (M-MSV) and M-MLV results in retroviral vectors that replicate as well as SNV vectors with wild type SNV E. This result indicates that the encapsidation sequences of M-MSV/M-MLV and SNV are not virus specific and that, during packaging of SNV and MLV RNA with viral proteins from REV-A, the encapsidation sequences are recognized largely by their secondary or tertiary structures.     

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.