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.     

Selection of an RNA domain that binds Zn2+.

We have selected an RNA that depends on zinc for affinity to a column, starting from a pool of ribooligonucleotides with 50 randomized positions. This RNA's chemical sensitivities, calculated folding thermodynamics, and activity when fragmented suggest that an ion binding site lies within a complex 21-nt hairpin loop, near the junction with an imperfect helical stem. This RNA site has an unselected selectivity among divalents, preferring nickel, cobalt, and cadmium to calcium, magnesium, and manganese, as expected for a simple site of chelation. A moderate zinc-dependent change in loop structure accompanies divalent binding and can be detected by chemical probing and zinc-dependent UV-induced crosslinking. The latter also demonstrates the apposition of loop sequences to make a structure that may be related to the E-loop motif found in a number of other RNA molecules; the E-loop motif, accordingly, may be a divalent site.     

Contribution of the C-terminal tail of U1A RBD1 to RNA recognition and protein stability.

The N-terminal RNA binding domain (RBD1) of the human U1A protein interacts specifically with a short RNA hairpin containing the U1 snRNA stem/loop II sequence. Previous RNA binding studies have suggested that the C-terminal tail of RBD1 contributes to RNA recognition in addition to interactions on the beta-sheet surface of the protein. To evaluate the contributions of these C-terminal residues in RBD1 to RNA binding affinity and specificity, as well as to study the thermodynamic stability of RBDs, a number of RBD1 mutants with truncated tails, with single amino acid substitutions, and with both a truncation and an amino acid substitution, have been constructed. The thermodynamic stabilities of these mutants have been measured and compared by GdnHCI unfolding experiments. The RNA binding affinity and specificity of these mutant proteins have been assessed by measuring the binding of each protein to the wild-type RNA hairpin and to selected RNA mutants with nucleotide substitutions in the RNA loop. The results demonstrate first that, although the C-terminal tail of RBD1 makes significant contributions to RNA binding affinity, it is not required for RNA binding, and second, its contributions to binding specificity are mediated only through selected nucleotides in the RNA loop, for in the absence of the tail, the protein continues to use other nucleotides to discriminate among RNAs. In these truncated proteins, the secondary structure intrinsic to the C-terminal tail is absent, yet their affinity and discrimination for RNAs are not lost. Thus, a structured tail is not required for RNA recognition.     

Screening for engineered neomycin riboswitches that control translation initiation

Riboswitches are genetic control elements that regulate gene expression in a small molecule-dependent way. We developed a two-stage strategy of in vitro selection followed by a genetic screen and identified several artificial small molecule-binding riboswitches that respond to the aminoglycoside neomycin. Structure-function relationships and structural probing revealed that they adopt the general neomycin-binding motif. They display no sequence similarities to in vitro selected neomycin aptamers but contain parts of the decoding site that is the binding site for neomycin on the ribosomal RNA. We propose a model of a composed binding pocket of an internal loop as primary docking site and a terminal flaplike loop structure fixing neomycin in a sandwich-like manner. Such binding pockets characterized by multiple contacts between ligand and RNA are described for both natural and engineered riboswitches. We anticipate that combination of in vitro selection and in vivo screening is a useful strategy to identify RNA molecules with a desired functionality.     

Selection of high affinity RNA ligands to the bacteriophage R17 coat protein.

RNA ligands with high affinity for the bacteriophage R17 coat protein were isolated from a pool of random RNA molecules using SELEX. Of the 38 ligands isolated, 36 were found to contain a hairpin very similar to the naturally occurring coat protein binding site in the R17 genome. The common features of these 36 sequences provide a consensus binding site and predict components of a hairpin that promote favorable interaction with the coat protein. These include a tetraloop of primary sequence AUCA and a variable-length stem with a bulged adenosine residue at a specific stem position. The predicted consensus agrees well with the highest-affinity RNA binding site of the R17 coat protein, identified through classical but laborious techniques. These results demonstrate the value of SELEX as a tool for isolating high affinity RNA ligands to a specific target protein, and the further value of those ligands to point the researcher toward natural sequences for that target protein.     

In vitro selection of RNA against kanamycin B.

Aminoglycosides are well-known antibiotics that function by interacting with ribosomal RNA in bacteria. In order to understand the molecular details between RNA and the drug, RNA aptamer was selected against kanamycin B. After 12 cycles of selection, RNA was cloned and sequenced. Among 9 clones, sequences of three clones were identical, suggesting the selected RNA was enriched. Among the cloned RNA molecules, the triplicated RNA was the maximum binding RNA. It showed a 180 nM affinity (KD) to the cognate aminoglycoside, as measured by a surface plasmon resonance, and a competition assay using a fluorescence anisotropy technique. The affinity of the maximum binding RNA to a similar aminoglycoside, tobramycin, was much stronger than 12 nM of KD. The binding site of the aminoglycoside in the maximum binding RNA was a stem loop located at the end of the 5' region. A stem loop structural motif, found in this study, was similar to those previously reported, even though the sequences of the RNA were totally different from the known sequences of the aminoglycoside binding site of other aptamers. The present study suggests that the aminoglycoside-binding region in RNA does not have a sequence specificity, but has a shape-specific bulged stem loop, even though it has a nanomolar affinity.     

Specificity in the binding of aminoglycosides to HIV-RRE RNA.

Quantitative studies of the binding of neomycin B to RRE constructs are carried out to determine the relationship between non-Watson Crick base-paired elements in the RNA and aminoglycoside binding. The RRE region contains two unpaired domains containing a single base bulge and a bubble structure, respectively. Deletion of the single base bulge has no effect on neomycin binding as the site of aminoglycoside binding is localized to the bubble region. Converting the bubble region into an A-form duplex gradually abolishes neomycin B binding in 3-5-fold steps in affinity over a 75-fold range. Thus, the binding of aminoglycoside is favored at domains in RNA that are nonduplex in nature, but aminoglycoside binding is only graded-specific in that affinities are enhanced gradually as the structure further deviates from a duplex form. It is likely that high-affinity aminoglycoside binding does not occur in duplex RNA because the major groove is too narrow to allow for aminoglycoside access and that structural perturbations that allow widening of the groove facilitate access. However, these interactions are only graded-specific with respect to both aminoglycoside structure and RNA domain structure.     

Discovering ligands for a microRNA precursor with peptoid microarrays

We have screened peptoid microarrays to identify specific ligands for the RNA hairpin precursor of miR-21, a microRNA involved in cancer and heart disease. Microarrays were printed by spotting a library of 7680 N-substituted oligoglycines (peptoids) onto glass slides. Two compounds on the array specifically bind RNA having the sequence and predicted secondary structure of the miR-21 precursor hairpin and have specific affinity for the target in solution. Their binding induces a conformational change around the hairpin loop, and the most specific compound recognizes the loop sequence and a bulged uridine in the proximal duplex. Functional groups contributing affinity and specificity were identified, and by varying a critical methylpyridine group, a compound with a dissociation constant of 1.9 μM for the miR-21 precursor hairpin and a 20-fold discrimination against a closely-related hairpin was created. This work describes a systematic approach to discovery of ligands for specific pre-defined novel RNA structures. It demonstrates discovery of new ligands for an RNA for which no specific lead compounds were previously known by screening a microarray of small molecules.     

A more complex isoleucine aptamer with a cognate triplet

The simplest RNA that can meet a column affinity selection for isoleucine was previously defined using selection amplification with decreasing numbers of randomized nucleotides. This simplest UAUU motif was a small asymmetric internal loop. Conserved positions of the loop include isoleucine codon and anticodon triplets (Lozupone C., Changayil, S., Majerfeld, I., and Yarus, M. (2003) RNA (N. Y.) 9, 1315-1322). Using new primer sequences, we now select a somewhat more complex isoleucine binding RNA, requiring 4.7 more bits of information to describe. The newly selected structure is a terminal or hairpin loop of 20 nucleotides, 15 being invariant. An information profile shows that the new binding site contains five short functional loop regions joined by less significant single nucleotide positions. Among the important nucleotides is a conserved isoleucine anticodon, supporting the escaped triplet theory, which posits a stereochemical genetic code originating in RNA amino acid binding sites.     

RNA recognition by the DNA end-binding Ku heterodimer

Most nucleic acid-binding proteins selectively bind either DNA or RNA, but not both nucleic acids. The Saccharomyces cerevisiae Ku heterodimer is unusual in that it has two very different biologically relevant binding modes: (1) Ku is a sequence-nonspecific double-stranded DNA end-binding protein with prominent roles in nonhomologous end-joining and telomeric capping, and (2) Ku associates with a specific stem–loop of TLC1, the RNA subunit of budding yeast telomerase, and is necessary for proper nuclear localization of this ribonucleoprotein enzyme. TLC1 RNA-binding and dsDNA-binding are mutually exclusive, so they may be mediated by the same site on Ku. Although dsDNA binding by Ku is well studied, much less is known about what features of an RNA hairpin enable specific recognition by Ku. To address this question, we localized the Ku-binding site of the TLC1 hairpin with single-nucleotide resolution using phosphorothioate footprinting, used chemical modification to identify an unpredicted motif within the hairpin secondary structure, and carried out mutagenesis of the stem–loop to ascertain the critical elements within the RNA that permit Ku binding. Finally, we provide evidence that the Ku-binding site is present in additional budding yeast telomerase RNAs and discuss the possibility that RNA binding is a conserved function of the Ku heterodimer.     

Modular construction for function of a ribonucleoprotein enzyme: the catalytic domain of Bacillus subtilis RNase P complexed with B. subtilis RNase P protein

The bacterial RNase P holoenzyme catalyzes the formation of the mature 5′-end of tRNAs and is composed of an RNA and a protein subunit. Among the two folding domains of the RNase P RNA, the catalytic domain (C-domain) contains the active site of this ribozyme. We investigated specific binding of the Bacillus subtilis C-domain with the B.subtilis RNase P protein and examined the catalytic activity of this C-domain–P protein complex. The C-domain forms a specific complex with the P protein with a binding constant of ~0.1 µM. The C-domain–P protein complex and the holoenzyme are equally efficient in cleaving single-stranded RNA (~0.9 min^–1 at pH 7.8) and substrates with a hairpin–loop 3′ to the cleavage site (~40 min^–1). The holoenzyme reaction is much more efficient with a pre-tRNA substrate, binding at least 100-fold better and cleaving 10–500 times more efficiently. These results demonstrate that the RNase P holoenzyme is functionally constructed in three parts. The catalytic domain alone contains the active site, but has little specificity and affinity for most substrates. The specificity and affinity for the substrate is generated by either the specificity domain of RNase P RNA binding to a T stem–loop-like hairpin or RNase P protein binding to a single-stranded RNA. This modular construction may be exploited to obtain RNase P-based ribonucleoprotein complexes with altered substrate specificity.     

Nucleotide-sequence-specific and non-specific interactions of T4 DNA polymerase with its own mRNA

The DNA-binding DNA polymerase (gp43) of phage T4 is also an RNA-binding protein that represses translation of its own mRNA. Previous studies implicated two segments of the untranslated 5′-leader of the mRNA in repressor binding, an RNA hairpin structure and the adjacent RNA to the 3′ side, which contains the Shine–Dalgarno sequence. Here, we show by in vitro gp43–RNA binding assays that both translated and untranslated segments of the mRNA contribute to the high affinity of gp43 to its mRNA target (translational operator), but that a Shine–Dalgarno sequence is not required for specificity. Nucleotide sequence specificity appears to reside solely in the operator’s hairpin structure, which lies outside the putative ribosome-binding site of the mRNA. In the operator region external to the hairpin, RNA length rather than sequence is the important determinant of the high binding affinity to the protein. Two aspects of the RNA hairpin determine specificity, restricted arrangement of purine relative to pyrimidine residues and an invariant 5′-AC-3′ in the unpaired (loop) segment of the RNA structure. We propose a generalized structure for the hairpin that encompasses these features and discuss possible relationships between RNA binding determinants of gp43 and DNA binding by this replication enzyme.     

The Affinity of an Oligodeoxynucleotide-Peptide Conjugate for an RNA Hairpin Loop Depends on Stereochemistry at the DNA-Peptide Junction

Abstract

Molecular models of an oligodeoxynucleotide-peptide conjugate complexed to an RNA hairpin loop were constructed to assess the effect of stereoisomerism at the point of attachment of the peptide to the oligodeoxynucleotide on the affinity of the conjugate for an RNA target. The peptide portion of the oligodeoxynucleotide-peptide conjugate, (L-lysine)₈, was covalently attached to the N-allyl group of (D)- or (L)-aspartic alcohol that was incorporated into the interior of an antisense oligodeoxynucleotide. The stereocenter in the oligodeoxynucleotide interior originates from either (D)- or (L)-aspartic alcohol. The oligodeoxynucleotide portion of the oligodeoxynucleotide-peptide conjugate forms Watson-Crick base pairs with the single-stranded RNA that flanks the RNA hairpin loop. The positively charged peptide makes specific electrostatic contacts with the negatively charged phosphate backbone of the RNA hairpin loop when attached to the N-allyl of (D)-aspartic alcohol but does not have the proper orientation to make these electrostatic contacts when attached to the N-allyl of (L)-aspartic alcohol. This modelling study emphasizes the importance of stereocontrol at the point of branching in synthesizing oligodeoxynucleotide-peptide conjugates for binding of RNA hairpin loops.     

Footprinting,circular dichroism and UV melting studies on neomycin B binding to the packaging region of human immunodeficiency virus type-1 RNA

We have studied the binding of neomycin to a 171mer RNA (ψ-RNA) from the packaging region of the LAI strain of human immunodeficiency virus type-1, HIV-1 (LAI). The RNase I footprinting studies reveal that the primary binding site for the drug is in stem–loop 1, which contains the dimer initiation site of HIV-1. Loading this site with neomycin causes a structural change in the RNA, allowing nucleotides in the neighboring stem–loop 2 to participate in the drug site. Drug binding to secondary sites induces structural changes in other stem–loops of the RNA. Footprinting plots, showing cutting at a site as a function of drug concentration, were analyzed using a two-state model to obtain relative site-specific binding constants. Circular dichroism measurements show that neomycin binding to ψ-RNA changes the intensity of the strong negative CD band at 208 nm, confirming that neomycin induces structural changes. Melting studies of the RNA showed melting transitions in the absence of drug at 28.2, 37.2, 47.4, 55.5 and 60.8°C. Only the first two were affected by drug binding, the reason for this being explained by our analysis.     

Specific RNA binding by Q beta coat protein

The interaction between the bacteriophage Q beta coat protein and its specific binding site on Q beta genomic RNA was characterized by using a nitrocellulose filter binding assay. Q beta coat protein bound to a synthetic 29-nucleotide RNA hairpin with an association constant of 400 microM-1 at 4 degrees C, 0.2 M ionic strength, pH 6.0. Complex formation had a broad pH optimum centered around pH 6.0 and was favored by both enthalpy and entropy. The salt dependence of Ka revealed that four to five ion pairs may be formed in the complex although approximately 80% of the free energy of complex formation is contributed by nonelectrostatic interactions. Truncation experiments revealed that coat protein binding required only the presence of a hairpin with an eight base pair stem and a three-base loop. Analysis of the binding properties of hairpin variants showed that the sequence of the stem was not important for coat protein recognition and only one of the three loop residues was essential. A bulged adenosine present in the coat protein binding site was not required for coat protein binding. Q beta coat protein binding specific is therefore primarily achieved by the structure and not by the sequence of the operator.     

Coupling of drug protonation to the specific binding of aminoglycosides to the A site of 16 S rRNA: elucidation of the number of drug amino groups involved and their identities

2-Deoxystreptamine (2-DOS) aminoglycoside antibiotics bind specifically to the central region of the 16S rRNA A site and interfere with protein synthesis. Recently, we have shown that the binding of 2-DOS aminoglycosides to an A site model RNA oligonucleotide is linked to the protonation of drug amino groups. Here, we extend these studies to define the number of amino groups involved as well as their identities. Specifically, we use pH-dependent 15N NMR spectroscopy to determine the pK(a) values of the amino groups in neomycin B, paromomycin I, and lividomycin A sulfate, with the resulting pK(a) values ranging from 6.92 to 9.51. For each drug, the 3-amino group was associated with the lowest pK(a), with this value being 6.92 in neomycin B, 7.07 in paromomycin I, and 7.24 in lividomycin A. In addition, we use buffer-dependent isothermal titration calorimetry (ITC) to determine the number of protons linked to the complexation of the three drugs with the A site model RNA oligomer at pH 5.5, 8.8, or 9.0. At pH 5.5, the binding of the three drugs to the host RNA is independent of drug protonation effects. By contrast, at pH 9.0, the RNA binding of paromomycin I and neomycin B is coupled to the uptake of 3.25 and 3.80 protons, respectively, with the RNA binding of lividomycin A at pH 8.8 being coupled to the uptake of 3.25 protons. A comparison of these values with the protonation states of the drugs predicted by our NMR-derived pK(a) values allows us to identify the specific drug amino groups whose protonation is linked to complexation with the host RNA. These determinations reveal that the binding of lividomycin A to the host RNA is coupled to the protonation of all five of its amino groups, with the RNA binding of paromomycin I and neomycin B being linked to the protonation of four and at least five amino groups, respectively. For paromomycin I, the protonation reactions involve the 1-, 3-, 2'-, and 2"'-amino groups, while, for neomycin B, the binding-linked protonation reactions involve at least the 1-, 3-, 2', 6'-, and 2"'-amino groups. Our results clearly identify drug protonation reactions as important thermodynamic participants in the specific binding of 2-DOS aminoglycosides to the A site of 16S rRNA.     

2-aminopurine fluorescence: discrimination between specific and unspecific ligand binding to the kissing-loop dimer of the HIV-1 RNA

The fluorescent 2-aminopurine probe (2-AP) incorporated into the loop of 23-mer RNA hairpin of HIV-1 genome dimerization initiation site (DIS) was used for discrimination of specific and unspecific binding of paromomycin and spermine to the kissing loop dimer (KD) formed in solution. While both ligands stabilized the KD RNA structure, only paromomycin binding resulted in significant increase of 2-AP fluorescence. These observations suggest that the 2-AP fluorescent RNA construct might be useful for selecting ligands specifically binding the HIV-1 kissing loop RNA dimer.     

RNA binding properties of the coat protein from bacteriophage GA.

The coat protein of bacteriophage GA, a group II RNA phage, binds to a small RNA hairpin corresponding to its replicase operator. Binding is specific, with a Ka of 71 microM -1. This interaction differs kinetically from the analogous coat protein-RNA hairpin interactions of other RNA phage and also deviates somewhat in its pH and salt dependence. Despite 46 of 129 amino acid differences between the GA and group I phage R17 coat proteins, the binding sites are fairly similar. The essential features of the GA coat protein binding site are a based-paired stem with an unpaired purine and a four nucleotide loop having an A at position -4 and a purine at -7. Unlike the group I phage proteins, the GA coat protein does not distinguish between two alternate positions for the unpaired purine and does not show high specificity for a pyrimidine at position -5 of the loop.     

Ni2+-binding RNA motifs with an asymmetric purine-rich internal loop and a G-A base pair.

RNA molecules with high affinity for immobilized Ni2+ were isolated from an RNA pool with 50 randomized positions by in vitro selection-amplification. The selected RNAs preferentially bind Ni2+ and Co2+ over other cations from first series transition metals. Conserved structure motifs, comprising about 15 nt, were identified that are likely to represent the Ni2+ binding sites. Two conserved motifs contain an asymmetric purine-rich internal loop and probably a mismatch G-A base pair. The structure of one of these motifs was studied with proton NMR spectroscopy and formation of the G-A pair at the junction of helix and internal loop was demonstrated. Using Ni2+ as a paramagnetic probe, a divalent metal ion binding site near this G-A base pair was identified. Ni2+ ions bound to this motif exert a specific stabilization effect. We propose that small asymmetric purine-rich loops that contain a G-A interaction may represent a divalent metal ion binding site in RNA.