The 1.3 A crystal structure of a biotin-binding pseudoknot and the basis for RNA molecular recognition

A pseudoknot-containing aptamer isolated from a pool of random sequence molecules has been shown previously to represent an optimal RNA solution to the problem of binding biotin. The affinity of this RNA molecule is nonetheless orders of magnitude weaker than that of its highly evolved protein analogs, avidin and streptavidin. To understand the structural basis for biotin binding and to compare directly strategies for ligand recognition available to proteins and RNA molecules, we have determined the 1.3 A crystal structure of the aptamer complexed with its ligand. Biotin is bound at the interface between the pseudoknot's stacked helices in a pocket defined almost entirely by base-paired nucleotides. In comparison to the protein avidin, the aptamer packs more tightly around the biotin headgroup and makes fewer contacts with its fatty acid tail. Whereas biotin is deeply buried within the hydrophobic core in the avidin complex, the aptamer relies on a combination of hydrated magnesium ions and immobilized water molecules to surround its ligand. In addition to demonstrating fundamentally different approaches to molecular recognition by proteins and RNA, the structure provides general insight into the mechanisms by which RNA function is mediated by divalent metals.     

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA''s secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.     

Evolution of aptamers with a new specificity and new secondary structures from an ATP aptamer

Small changes in target specificity can sometimes be achieved, without changing aptamer structure, through mutation of a few bases. Larger changes in target geometry or chemistry may require more radical changes in an aptamer. In the latter case, it is unknown whether structural and functional solutions can still be found in the region of sequence space close to the original aptamer. To investigate these questions, we designed an in vitro selection experiment aimed at evolving specificity of an ATP aptamer. The ATP aptamer makes contacts with both the nucleobase and the sugar. We used an affinity matrix in which GTP was immobilized through the sugar, thus requiring extensive changes in or loss of sugar contact, as well as changes in recognition of the nucleobase. After just five rounds of selection, the pool was dominated by new aptamers falling into three major classes, each with secondary structures distinct from that of the ATP aptamer. The average sequence identity between the original aptamer and new aptamers is 76%. Most of the mutations appear to play roles either in disrupting the original secondary structure or in forming the new secondary structure or the new recognition loops. Our results show that there are novel structures that recognize a significantly different ligand in the region of sequence space close to the ATP aptamer. These examples of the emergence of novel functions and structures from an RNA molecule with a defined specificity and fold provide a new perspective on the evolutionary flexibility and adaptability of RNA.     

Structural probing and damage selection of citrulline- and arginine-specific RNA aptamers identify base positions required for binding.

In a recent study, an RNA aptamer for the specific recognition of the amino acid L-arginine was evolved from an in vitro selected L-citrulline binding parent sequence [M. Famulok (1994) J. Am. Chem. Soc. 116, 1698-1706]. We have now carried out a structural analysis of these aptamers by using chemical modification experiments. Footprinting experiments and a damage selection approach were performed to identify those positions protected from modification in the presence of the amino acids and modifications that interfere with the binding of the ligand. It is shown that of the two bulged regions present in both aptamers one can be modified without loss of binding activity whereas in the other bulge nearly every position is shown to be involved in the recognition of the ligands. This might be indicative for non-canonical base pairing to occur within the non-Watson-Crick paired regions which might be stabilized by the complexed amino acid. Binding to the cognate amino acid significantly enhances the conformational stability of the RNA. We also tested the sensitivity of both aptamers towards lead (II) ion induced cleavage and identified a hypersensitive cleavage site within the invariant bulged region. Lead cleavage is inhibited by the complexed amino acid, indicating a conformational change of the aptamer upon ligand binding. NMR titration data obtained with both aptamers and their cognate ligands confirm the proposed conformational changes and indicate the formation of a 1:1 complex of RNA:amino acid.     

Biophysical characterization of DNA and RNA aptamer interactions with hen egg lysozyme

This work characterized the binding of an RNA aptamer recognizing hen egg white lysozyme, as well as a literature-reported single-stranded DNA analog of sequence identical to the original RNA aptamer, using fluorescence anisotropy, isothermal titration calorimetry (ITC) and analytical ultracentrifugation. The polyanionic DNA aptamer analog is selective for lysozyme even over cationic cytochrome c and has been reported to be successfully used in biosensing applications. The association however, is predominantly of electrostatic character, strongly salt-sensitive and entropically-driven, in contrast to previously described enthalpically-driven antibody-lysozyme and DNA aptamer-VEGF interactions. With a moderate selectivity for their target, high salt-sensitivity along with fast association and dissociation behavior, these molecules might serve as pseudo-affinity ligands for biomolecular separations.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Characterisation of aptamer–target interactions by branched selection and high-throughput sequencing of SELEX pools

Nucleic acid aptamer selection by systematic evolution of ligands by exponential enrichment (SELEX) has shown great promise for use in the development of research tools, therapeutics and diagnostics. Typically, aptamers are identified from libraries containing up to 10¹⁶ different RNA or DNA sequences by 5–10 rounds of affinity selection towards a target of interest. Such library screenings can result in complex pools of many target-binding aptamers. New high-throughput sequencing techniques may potentially revolutionise aptamer selection by allowing quantitative assessment of the dynamic changes in the pool composition during the SELEX process and by facilitating large-scale post-SELEX characterisation. In the present study, we demonstrate how high-throughput sequencing of SELEX pools, before and after a single round of branched selection for binding to different target variants, can provide detailed information about aptamer binding sites, preferences for specific target conformations, and functional effects of the aptamers. The procedure was applied on a diverse pool of 2′-fluoropyrimidine-modified RNA enriched for aptamers specific for the serpin plasminogen activator inhibitor-1 (PAI-1) through five rounds of standard selection. The results demonstrate that it is possible to perform large-scale detailed characterisation of aptamer sequences directly in the complex pools obtained from library selection methods, thus without the need to produce individual aptamers.     

Role of a heterogeneous free state in the formation of a specific RNA-theophylline complex

The helical regions of RNA are generally very stable, but the single-stranded and loop regions often exist as an ensemble of conformations in solution. The theophylline-binding RNA aptamer forms a very stable structure when bound to the bronchodilator theophylline, but the theophylline binding site is not stably formed in the absence of ligand. The kinetics for theophylline binding were measured here by stopped-flow fluorescence spectroscopy to probe the mechanism for theophylline binding in this RNA aptamer. The kinetic studies showed that formation of the RNA-theophylline complex is over 1000 times slower than a diffusion-controlled rate, and the high affinity of the RNA-theophylline complex arises primarily from a slow dissociation rate for the complex. A theophylline-independent rate was observed for formation of the theophylline-RNA complex at high theophylline concentration, indicating that a conformational change in the RNA is the rate-limiting step in complex formation under these conditions. The RNA-theophylline complex requires divalent metal ions, such as Mg2+, to form a high-affinity complex, and there is a greater than 10000-fold reduction in affinity for theophylline in the absence of Mg2+. This decrease in binding affinity in the absence of Mg2+ results primarily from an increased dissociation rate for the complex. The implications of an ensemble of conformations in the free state of this theophylline-binding RNA are discussed and compared with mechanisms for formation of protein-ligand complexes.     

Simple fluorescent sensors engineered with catalytic DNA 'MgZ' based on a non-classic allosteric design

Most NAE (nucleic acid enzyme) sensors are composed of an RNA-cleaving catalytic motif and an aptameric receptor. They operate by activating or repressing the catalytic activity of a relevant NAE through the conformational change in the aptamer upon target binding. To transduce a molecular recognition event to a fluorescence signal, a fluorophore-quencher pair is attached to opposite ends of the RNA substrate such that when the NAE cleaves the substrate, an increased level of fluorescence can be generated. However, almost all NAE sensors to date harbor either NAEs that cannot accommodate a fluorophore-quencher pair near the cleavage site or those that can accept such a modification but require divalent transition metal ions for catalysis. Therefore, the signaling magnitude and the versatility of current NAE sensors might not suffice for analytical and biological applications. Here we report an RNA-cleaving DNA enzyme, termed 'MgZ', which depends on Mg(2+) for its activity and can accommodate bulky dye moieties next to the cleavage site. MgZ was created by in vitro selection. The selection scheme entailed acidic buffering and ethanol-based reaction stoppage to remove selfish DNAs. Characterization of MgZ revealed a three-way junction structure, a cleavage rate of 1 min(-1), and 26-fold fluorescence enhancement. Two ligand-responsive NAE sensors were rationally designed by linking an aptamer sequence to the substrate of MgZ. In the absence of the target, the aptamer-linked substrate is locked into a conformation that prohibits MgZ from accessing the substrate. In the presence of the target, the aptamer releases the substrate, which induces MgZ-mediated RNA cleavage. The discovery of MgZ and the introduction of the above NAE sensor design strategy should facilitate future efforts in sensor engineering.     

Molecular recognition of amino acids by RNA aptamers: the evolution into an L-tyrosine binder of a dopamine-binding RNA motif

We report the evolution of an RNA aptamer to change its binding specificity. RNA aptamers that bind the free amino acid tyrosine were in vitro selected from a degenerate pool derived from a previously selected dopamine aptamer. Three independent sequences bind tyrosine in solution, the winner of the selection binding with a dissociation constant of 35 microM. Competitive affinity chromatography with tyrosine-related ligands indicated that the selected aptamers are highly L-stereo selective and also recognize L-tryptophan and L-dopa with similar affinity. The binding site was localized by sequence comparison, analysis of minimal boundaries, and structural probing upon ligand binding. Tyrosine-binding sites are characterized by the presence of both tyrosine (UAU and UAC) and termination (UAG and UAA) triplets.     

Dual-colour imaging of RNAs using quencher- and fluorophore-binding aptamers

In order to gain deeper insight into the functions and dynamics of RNA in cells, the development of methods for imaging multiple RNAs simultaneously is of paramount importance. Here, we describe a modular approach to image RNA in living cells using an RNA aptamer that binds to dinitroaniline, an efficient general contact quencher. Dinitroaniline quenches the fluorescence of different fluorophores when directly conjugated to them via ethylene glycol linkers by forming a non-fluorescent intramolecular complex. Since the binding of the RNA aptamer to the quencher destroys the fluorophore-quencher complex, fluorescence increases dramatically upon binding. Using this principle, a series of fluorophores were turned into fluorescent turn-on probes by conjugating them to dinitroaniline. These probes ranged from fluorescein-dinitroaniline (green) to TexasRed-dinitroaniline (red) spanning across the visible spectrum. The dinitroaniline-binding aptamer (DNB) was generated by in vitro selection, and was found to bind all probes, leading to fluorescence increase in vitro and in living cells. When expressed in E. coli, the DNB aptamer could be labelled and visualized with different-coloured fluorophores and therefore it can be used as a genetically encoded tag to image target RNAs. Furthermore, combining contact-quenched fluorogenic probes with orthogonal DNB (the quencher-binding RNA aptamer) and SRB-2 aptamers (a fluorophore-binding RNA aptamer) allowed dual-colour imaging of two different fluorescence-enhancing RNA tags in living cells, opening new avenues for studying RNA co-localization and trafficking.     

Molecular interactions and metal binding in the theophylline-binding core of an RNA aptamer

An RNA aptamer containing a 15-nt binding site shows high affinity and specificity for the bronchodilator theophylline. A variety of base modifications or 2' deoxyribose substitutions in binding-site residues were tested for theophyllinebinding affinity and the results were compared with the previously determined three-dimensional structure of the RNA-theophylline complex. The RNA-theophylline complex contains a U6-A28-U23 base triple, and disruption of this A28-U23 Hoogsteen-pair by a 7-deaza, 2'-deoxy A28 mutant reduces theophylline binding >45-fold at 25 degrees C. U24 is part of a U-turn in the core of the RNA, and disruption of this U-turn motif by a 2'-deoxy substitution of U24 also reduces theophylline binding by >90-fold. Several mutations outside the "conserved core" of the RNA aptamer showed reduced binding affinity, and these effects could be rationalized by comparison with the three-dimensional structure of the complex. Divalent ions are absolutely required for high-affinity theophylline binding. High-affinity binding was observed with 5 mM Mg2+, Mn2+, or Co2+ ions, whereas little or no significant binding was observed for other divalent or lanthanide ions. A metal-binding site in the core of the complex was revealed by paramagnetic Mn2+-induced broadening of specific RNA resonances in the NMR spectra. When caffeine is added to the aptamer in tenfold excess, the NMR spectra show no evidence for binding in the conserved core and instead the drug stacks on the terminal helix. The lack of interaction between caffeine and the theophylline-binding site emphasizes the extreme molecular discrimination of this RNA aptamer.     

Solution structure of an informationally complex high-affinity RNA aptamer to GTP

Higher-affinity RNA aptamers to GTP are more informationally complex than lower-affinity aptamers. Analog binding studies have shown that the additional information needed to improve affinity does not specify more interactions with the ligand. In light of those observations, we would like to understand the structural characteristics that enable complex aptamers to bind their ligands with higher affinity. Here we present the solution structure of the 41-nt Class I GTP aptamer (K(d) = 75 nM) as determined by NMR. The backbone of the aptamer forms a reverse-S that shapes the binding pocket. The ligand nucleobase stacks between purine platforms and makes hydrogen bonds with the edge of another base. Interestingly, the local modes of interaction for the Class I aptamer and an RNA aptamer that binds ATP with a K(d) of 6 microM are very much alike. The aptamers exhibit nearly identical levels of binding specificity and fraction of ligand sequestered from the solvent (81%-85%). However, the GTP aptamer is more informationally complex (approximately 45 vs. 35 bits) and has a larger recognition bulge (15 vs. 12 nucleotides) with many more stabilizing base-base interactions. Because the aptamers have similar modes of ligand binding, we conclude that the stabilizing structural elements in the Class I aptamer are responsible for much of the difference in K(d). These results are consistent with the hypothesis that increasing the number of intra-RNA interactions, rather than adding specific contacts to the ligand, is the simplest way to improve binding affinity.     

Recognition of a cognate RNA aptamer by neomycin B: quantitative evaluation of hydrogen bonding and electrostatic interactions

Aminoglycosides are an important class of antibiotic that selectively target RNA structural motifs. Recently we have demonstrated copper derivatives of amino-glycosides to be efficient cleavage agents for cognate RNA motifs. To fully develop their potential as pharmaceutical agents it is necessary to understand both the structural mechanisms used by aminoglycosides to target RNA, and the relative contributions of hydrogen bonding and electrostatic interactions to recognition selectivity. Herein we report results from a calorimetric analysis of a stem-loop 23mer RNA aptamer complexed to the aminoglycoside neomycin B. Key thermodynamic parameters for complex formation have been determined by isothermal titration calorimetry, and from the metal-ion dependence of these binding parameters the relative contributions of electrostatics and hydrogen bonding toward binding affinity have been assessed. The principal mechanism for recognition and binding of neomycin B to the RNA major groove is mediated by hydrogen bonding.     

Investigating the role of metal ions in the catalytic mechanism of the yeast RNA triphosphatase

The Saccharomyces cerevisiae RNA triphosphatase (Cet1) requires the presence of metal ion cofactors to catalyze its phosphohydrolase activity, the first step in the formation of the 5'-terminal cap structure of mRNAs. We have used endogenous tryptophan fluorescence studies to elucidate both the nature and the role(s) of the metal ions in the Cet1-mediated phosphohydrolase reaction. The association of Mg2+, Mn2+, and Co2+ ions with the enzyme resulted in a decrease in the intensity of the tryptophan emission spectrum. This decrease was then used to determine the apparent dissociation constants for these ions. Subsequent dual ligand titration experiments demonstrated that the metal ions bind to a common site, for which they compete. The kinetics of real-time metal ion binding to the Cet1 protein were also investigated, and the effects on RNA and nucleotide binding were evaluated. To provide additional insight into the relationship between Cet1 structure and metal ion binding, we correlated the effect of ion binding on protein structure using both circular dichroism and guanidium hydrochloride-induced denaturation as structural indicators. Our data indicate that binding of RNA, nucleotides, and metal ion cofactors does not lead to significant structural modifications of the Cet1 architecture. This suggests a model in which Cet1 possesses a preformed active site, and where major domain rearrangements are not required to form an active catalytic site. Finally, denaturation studies demonstrate that the metal ion cofactors can act by stabilizing the ground state binding of the phosphohydrolase substrate.     

Molecular dynamics simulation on the allosteric analysis of the c‐di‐GMP class I riboswitch induced by ligand binding

Riboswitches are RNA molecules that regulate gene expression using conformation change, affected by binding of small molecule ligands. Although a number of ligand‐bound aptamer complex structures have been solved, it is important to know ligand‐free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, we use dynamics simulations on a series of models to characterize the ligand‐free and ligand‐bound aptamer domain of the c‐di‐GMP class I (GEMM‐I) riboswitch. The results revealed that the ligand‐free aptamer has a stable state with a folded P2 and P3 helix, an unfolded P1 helix and open binding pocket. The first Mg ions binding to the aptamer is structurally favorable for the successive c‐di‐GMP binding. The P1 helix forms when c‐di‐GMP is successive bound. Three key junctions J1/2, J2/3 and J1/3 in the GEMM‐I riboswitch contributing to the formation of P1 helix have been found. The binding of the c‐di‐GMP ligand to the GEMM‐I riboswitch induces the riboswitch's regulation through the direct allosteric communication network in GEMM‐I riboswitch from the c‐di‐GMP binding sites in the J1/2 and J1/3 junctions to the P1 helix, the indirect ones from those in the J2/3 and P2 communicating to P1 helix via the J1/2 and J1/3 media.     

Mutational analysis of the purine riboswitch aptamer domain

The purine riboswitch is one of a number of mRNA elements commonly found in the 5'-untranslated region capable of controlling expression in a cis-fashion via its ability to directly bind small-molecule metabolites. Extensive biochemical and structural analysis of the nucleobase-binding domain of the riboswitch, referred to as the aptamer domain, has revealed that the mRNA recognizes its cognate ligand using an intricately folded three-way junction motif that completely encapsulates the ligand. High-affinity binding of the purine nucleobase is facilitated by a distal loop-loop interaction that is conserved between both the adenine and guanine riboswitches. To understand the contribution of conserved nucleotides in both the three-way junction and the loop-loop interaction of this RNA, we performed a detailed mutagenic survey of these elements in the context of an adenine-responsive variant of the xpt-pbuX guanine riboswitch from Bacillus subtilis. The varying ability of these mutants to bind ligand as measured by isothermal titration calorimetry uncovered the conserved nucleotides whose identity is required for purine binding. Crystallographic analysis of the bound form of five mutants and chemical probing of their free state demonstrate that the identity of several universally conserved nucleotides is not essential for formation of the RNA-ligand complex but rather for maintaining a binding-competent form of the free RNA. These data show that conservation patterns in riboswitches arise from a combination of formation of the ligand-bound complex, promoting an open form of the free RNA, and participating in the secondary structural switch with the expression platform.