Two riboswitch classes that share a common ligand-binding fold show major differences in the ability to accommodate mutations.

Riboswitches are structured RNAs that sense small molecules to control expression. Prequeuosine1 (preQ1)-sensing riboswitches comprise three classes (I, II and III) that adopt distinct folds. Despite this difference, class II and III riboswitches each use 10 identical nucleotides to bind the preQ1 metabolite. Previous class II studies showed high sensitivity to binding-pocket mutations, which reduced preQ1 affinity and impaired function. Here, we introduced four equivalent mutations into a class III riboswitch, which maintained remarkably tight preQ1 binding. Co-crystal structures of each class III mutant showed compensatory interactions that preserve the fold. Chemical modification analysis revealed localized RNA flexibility changes for each mutant, but molecular dynamics (MD) simulations suggested that each mutation was not overtly destabilizing. Although impaired, class III mutants retained tangible gene-regulatory activity in bacteria compared to equivalent preQ1-II variants; mutations in the preQ1-pocket floor were tolerated better than wall mutations. Principal component analysis of MD trajectories suggested that the most functionally deleterious wall mutation samples different motions compared to wildtype. Overall, the results reveal that formation of compensatory interactions depends on the context of mutations within the overall fold and that functionally deleterious mutations can alter long-range correlated motions that link the riboswitch binding pocket with distal gene-regulatory sequences.

A riboswitch separated from its ribosome-binding site still regulates translation.

Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10-G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch-akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

Structure and function analysis of a type III preQ1-I riboswitch from Escherichia coli reveals direct metabolite sensing by the Shine-Dalgarno sequence.

Riboswitches are small noncoding RNAs found primarily in the 5' leader regions of bacterial messenger RNAs where they regulate expression of downstream genes in response to binding one or more cellular metabolites. Such noncoding RNAs are often regulated at the translation level, which is thought to be mediated by the accessibility of the Shine-Dalgarno sequence (SDS) ribosome-binding site. Three classes (I-III) of prequeuosine1 (preQ1)-sensing riboswitches are known that control translation. Class I is divided into three subtypes (types I-III) that have diverse mechanisms of sensing preQ1, which is involved in queuosine biosynthesis. To provide insight into translation control, we determined a 2.30 Å-resolution cocrystal structure of a class I type III preQ1-sensing riboswitch identified in Escherichia coli (Eco) by bioinformatic searches. The Eco riboswitch structure differs from previous preQ1 riboswitch structures because it has the smallest naturally occurring aptamer and the SDS directly contacts the preQ1 metabolite. We validated structural observations using surface plasmon resonance and in vivo gene-expression assays, which showed strong switching in live E. coli. Our results demonstrate that the Eco riboswitch is relatively sensitive to mutations that disrupt noncanonical interactions that form the pseudoknot. In contrast to type II preQ1 riboswitches, a kinetic analysis showed that the type III Eco riboswitch strongly prefers preQ1 over the chemically similar metabolic precursor preQ0. Our results reveal the importance of noncanonical interactions in riboswitch-driven gene regulation and the versatility of the class I preQ1 riboswitch pseudoknot as a metabolite-sensing platform that supports SDS sequestration.

Full-Length NAD+-I Riboswitches Bind a Single Cofactor but Cannot Discriminate against Adenosine Triphosphate.

Bacterial riboswitches are structured RNAs that bind small metabolites to control downstream gene expression. Two riboswitch classes have been reported to sense nicotinamide adenine dinucleotide (NAD+), which plays a key redox role in cellular metabolism. The NAD+-I (class I) riboswitch stands out because it comprises two homologous, tandemly arranged domains. However, previous studies examined the isolated domains rather than the full-length riboswitch. Crystallography and ligand binding analyses led to the hypothesis that each domain senses NAD+ but with disparate equilibrium binding constants (KD) of 127 µM (domain I) and 3.4 mM (domain II). Here, we analyzed individual domains and the full-length riboswitch by isothermal titration calorimetry to quantify the cofactor affinity and specificity. Domain I senses NAD+ with a KD of 24.6 ± 8.4 µM but with a reduced ligand-to-receptor stoichiometry, consistent with nonproductive domain self-association observed by gel-filtration chromatography; domain II revealed no detectable binding. By contrast, the full-length riboswitch binds a single NAD+ with a KD of 31.5 ± 1.5 µM; dinucleotides NADH and AP2-ribavirin also bind with one-to-one stoichiometry. Unexpectedly, the full-length riboswitch also binds a single ATP equivalent (KD = 11.0 ± 3.5 µM). The affinity trend of the full-length riboswitch is ADP = ATP > NAD+ = AP2-ribavirin > NADH. Although our results support riboswitch sensing of a single NAD+ at concentrations significantly below the intracellular levels of this cofactor, our findings do not support the level of specificity expected for a riboswitch that exclusively senses NAD+. Gene regulatory implications and future challenges are discussed.

RNA target highlights in CASP15: Evaluation of predicted models by structure providers.

The first RNA category of the Critical Assessment of Techniques for Structure Prediction competition was only made possible because of the scientists who provided experimental structures to challenge the predictors. In this article, these scientists offer a unique and valuable analysis of both the successes and areas for improvement in the predicted models. All 10 RNA-only targets yielded predictions topologically similar to experimentally determined structures. For one target, experimentalists were able to phase their x-ray diffraction data by molecular replacement, showing a potential application of structure predictions for RNA structural biologists. Recommended areas for improvement include: enhancing the accuracy in local interaction predictions and increased consideration of the experimental conditions such as multimerization, structure determination method, and time along folding pathways. The prediction of RNA-protein complexes remains the most significant challenge. Finally, given the intrinsic flexibility of many RNAs, we propose the consideration of ensemble models.

Cyclic peptides with a distinct arginine-fork motif recognize the HIV trans-activation response RNA in vitro and in cells.

RNA represents a potential target for new antiviral therapies, which are urgently needed to address public health threats such as the human immunodeficiency virus (HIV). We showed previously that the interaction between the viral Tat protein and the HIV-1 trans-activation response (TAR) RNA was blocked by TB-CP-6.9a. This cyclic peptide was derived from a TAR-binding loop that emerged during lab evolution of a TAR-binding protein (TBP) family. Here we synthesized and characterized a next-generation, cyclic-peptide library based on the TBP scaffold. We sought to identify conserved RNA-binding interactions and the influence of cyclization linkers on RNA binding and antiviral activity. A diverse group of cyclization linkers, encompassing disulfide bonds to bicyclic aromatic staples, was used to restrain the cyclic peptide geometry. Thermodynamic profiling revealed specific arginine-rich sequences with low to submicromolar affinity driven by enthalpic and entropic contributions. The best compounds exhibited no appreciable off-target binding to related molecules, such as BIV TAR and human 7SK RNAs. A specific arginine-to-lysine change in the highest affinity cyclic peptide reduced TAR binding by tenfold, suggesting that TBP-derived cyclic peptides use an arginine-fork motif to recognize the TAR major groove while differentiating the mode of binding from other TAR-targeting molecules. Finally, we showed that HIV infectivity in cell culture was reduced in the presence of cyclic peptides constrained by methylene or naphthalene-based linkers. Our findings provide insight into the molecular determinants required for HIV-1 TAR recognition and antiviral activity. These findings are broadly relevant to the development of antivirals that target RNA molecules.

Analysis of a preQ1-I riboswitch in effector-free and bound states reveals a metabolite-programmed nucleobase-stacking spine that controls gene regulation.

Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1-exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1-sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor.

Nucleobase mutants of a bacterial preQ1-II riboswitch that uncouple metabolite sensing from gene regulation.

Riboswitches are a class of nonprotein-coding RNAs that directly sense cellular metabolites to regulate gene expression. They are model systems for analyzing RNA-ligand interactions and are established targets for antibacterial agents. Many studies have analyzed the ligand-binding properties of riboswitches, but this work has outpaced our understanding of the underlying chemical pathways that govern riboswitch-controlled gene expression. To address this knowledge gap, we prepared 15 mutants of the preQ1-II riboswitch-a structurally and biochemically well-characterized HLout pseudoknot that recognizes the metabolite prequeuosine1 (preQ1). The mutants span the preQ1-binding pocket through the adjoining Shine-Dalgarno sequence (SDS) and include A-minor motifs, pseudoknot-insertion helix P4, U·A-U base triples, and canonical G-C pairs in the anti-SDS. As predicted-and confirmed by in vitro isothermal titration calorimetry measurements-specific mutations ablated preQ1 binding, but most aberrant binding effects were corrected by compensatory mutations. In contrast, functional analysis in live bacteria using a riboswitch-controlled GFPuv-reporter assay revealed that each mutant had a deleterious effect on gene regulation, even when compensatory changes were included. Our results indicate that effector binding can be uncoupled from gene regulation. We attribute loss of function to defects in a chemical interaction network that links effector binding to distal regions of the fold that support the gene-off RNA conformation. Our findings differentiate effector binding from biological function, which has ramifications for riboswitch characterization. Our results are considered in the context of synthetic ligands and drugs that bind tightly to riboswitches without eliciting a biological response.

Co-crystal structures of HIV TAR RNA bound to lab-evolved proteins show key roles for arginine relevant to the design of cyclic peptide TAR inhibitors.

RNA-protein interfaces control key replication events during the HIV-1 life cycle. The viral trans-activator of transcription (Tat) protein uses an archetypal arginine-rich motif (ARM) to recruit the host positive transcription elongation factor b (pTEFb) complex onto the viral trans-activation response (TAR) RNA, leading to activation of HIV transcription. Efforts to block this interaction have stimulated production of biologics designed to disrupt this essential RNA-protein interface. Here, we present four co-crystal structures of lab-evolved TAR-binding proteins (TBPs) in complex with HIV-1 TAR. Our results reveal that high-affinity binding requires a distinct sequence and spacing of arginines within a specific ß2-ß3 hairpin loop that arose during selection. Although loops with as many as five arginines were analyzed, only three arginines could bind simultaneously with major-groove guanines. Amino acids that promote backbone interactions within the ß2-ß3 loop were also observed to be important for high-affinity interactions. Based on structural and affinity analyses, we designed two cyclic peptide mimics of the TAR-binding ß2-ß3 loop sequences present in two high-affinity TBPs (KD values of 4.2 ± 0.3 and 3.0 ± 0.3 nm). Our efforts yielded low-molecular weight compounds that bind TAR with low micromolar affinity (KD values ranging from 3.6 to 22 µm). Significantly, one cyclic compound within this series blocked binding of the Tat-ARM peptide to TAR in solution assays, whereas its linear counterpart did not. Overall, this work provides insight into protein-mediated TAR recognition and lays the ground for the development of cyclic peptide inhibitors of a vital HIV-1 RNA-protein interaction.

Face-time with TAR: Portraits of an HIV-1 RNA with diverse modes of effector recognition relevant for drug discovery.

Small molecules and short peptides that potently and selectively bind RNA are rare, making the molecular structures of these complexes highly exceptional. Accordingly, several recent investigations have provided unprecedented structural insights into how peptides and proteins recognize the HIV-1 transactivation response (TAR) element, a 59-nucleotide-long, noncoding RNA segment in the 5' long terminal repeat region of viral transcripts. Here, we offer an integrated perspective on these advances by describing earlier progress on TAR binding to small molecules, and by drawing parallels to recent successes in the identification of compounds that target the hepatitis C virus internal ribosome entry site (IRES) and the flavin-mononucleotide riboswitch. We relate this work to recent progress that pinpoints specific determinants of TAR recognition by: (i) viral Tat proteins, (ii) an innovative lab-evolved TAR-binding protein, and (iii) an ultrahigh-affinity cyclic peptide. New structural details are used to model the TAR-Tat-super-elongation complex (SEC) that is essential for efficient viral transcription and represents a focal point for antiviral drug design. A key prediction is that the Tat transactivation domain makes modest contacts with the TAR apical loop, whereas its arginine-rich motif spans the entire length of the TAR major groove. This expansive interface has significant implications for drug discovery and design, and it further suggests that future lab-evolved proteins could be deployed to discover steric restriction points that block Tat-mediated recruitment of the host SEC to HIV-1 TAR.

Arginine Forks Are a Widespread Motif to Recognize Phosphate Backbones and Guanine Nucleobases in the RNA Major Groove.

RNA recognition by proteins is central to biology. Here we demonstrate the existence of a recurrent structural motif, the "arginine fork", that codifies arginine readout of cognate backbone and guanine nucleobase interactions in a variety of protein-RNA complexes derived from viruses, metabolic enzymes, and ribosomes. Nearly 30 years ago, a theoretical arginine fork model was posited to account for the specificity between the HIV-1 Tat protein and TAR RNA. This model predicted that a single arginine should form four complementary contacts with nearby phosphates, yielding a two-pronged backbone readout. Recent high-resolution structures of TAR-protein complexes have unveiled new details, including (i) arginine interactions with the phosphate backbone and the major-groove edge of guanine and (ii) simultaneous cation-π contacts between the guanidinium group and flanking nucleobases. These findings prompted us to search for arginine forks within experimental protein-RNA structures retrieved from the Protein Data Bank. The results revealed four distinct classes of arginine forks that we have defined using a rigorous but flexible nomenclature. Examples are presented in the context of ribosomal and nonribosomal interfaces with analysis of arginine dihedral angles and structural (suite) classification of RNA targets. When arginine fork chemical recognition principles were applied to existing structures with unusual arginine-guanine recognition, we found that the arginine fork geometry was more consistent with the experimental data, suggesting the utility of fork classifications to improve structural models. Software to analyze arginine-RNA interactions has been made available to the community.

Structure of HIV TAR in complex with a Lab-Evolved RRM provides insight into duplex RNA recognition and synthesis of a constrained peptide that impairs transcription.

Natural and lab-evolved proteins often recognize their RNA partners with exquisite affinity. Structural analysis of such complexes can offer valuable insight into sequence-selective recognition that can be exploited to alter biological function. Here, we describe the structure of a lab-evolved RNA recognition motif (RRM) bound to the HIV-1 trans-activation response (TAR) RNA element at 1.80 Å-resolution. The complex reveals a trio of arginines in an evolved ß2-ß3 loop penetrating deeply into the major groove to read conserved guanines while simultaneously forming cation-π and salt-bridge contacts. The observation that the evolved RRM engages TAR within a double-stranded stem is atypical compared to most RRMs. Mutagenesis, thermodynamic analysis and molecular dynamics validate the atypical binding mode and quantify molecular contributions that support the exceptionally tight binding of the TAR-protein complex (KD,App of 2.5 ± 0.1 nM). These findings led to the hypothesis that the ß2-ß3 loop can function as a standalone TAR-recognition module. Indeed, short constrained peptides comprising the ß2-ß3 loop still bind TAR (KD,App of 1.8 ± 0.5 µM) and significantly weaken TAR-dependent transcription. Our results provide a detailed understanding of TAR molecular recognition and reveal that a lab-evolved protein can be reduced to a minimal RNA-binding peptide.

Observation of preQ1-II riboswitch dynamics using single-molecule FRET.

PreQ1 riboswitches regulate the synthesis of the hypermodified tRNA base queuosine by sensing the pyrrolopyrimidine metabolite preQ1. Here, we use single-molecule FRET to interrogate the structural dynamics of apo and preQ1-bound states of the preQ1-II riboswitch from Lactobacillus rhamnosus. We find that the apo-form of the riboswitch spontaneously samples multiple conformations. Magnesium ions and preQ1 stabilize conformations that sequester the ribosome-binding site of the mRNA within the pseudoknotted structure, thus inhibiting translation initiation. Our results reveal that folding of the preQ1-II riboswitch is complex and provide evidence favoring a conformational selection model of effector binding by riboswitches of this class.

Coupling Green Fluorescent Protein Expression with Chemical Modification to Probe Functionally Relevant Riboswitch Conformations in Live Bacteria.

Noncoding RNAs engage in numerous biological activities including gene regulation. To fully understand RNA function it is necessary to probe biologically relevant conformations in living cells. To address this challenge, we coupled RNA-mediated regulation of the green fluorescent protein (GFP)uv-reporter gene to icSHAPE (in cell Selective 2'-Hydroxyl Acylation analyzed by Primer Extension). Our transcript-specific approach provides sensitive, fluorescence-based readout of the regulatory-RNA status as a means to coordinate chemical modification experiments. We chose a plasmid-based reporter compatible with Escherichia coli to allow use of knockout strains that eliminate endogenous effector biosynthesis. The approach was piloted using the Lactobacillus rhamnosus ( Lrh) preQ1-II riboswitch, which senses the pyrrolopyrimidine metabolite preQ1. Using an E. coli Δ queF strain incapable of preQ1 anabolism, the Lrh riboswitch yielded nearly one log unit of GFPuv-gene repression resulting from exogenously added preQ1. We then subjected cells in gene "on" and "off" states to icSHAPE. The resulting differential analysis indicated reduction in Lrh riboswitch flexibility in the P3 helix of the pseudoknot, which comprises the ribosome-binding site (RBS) paired with the anti-RBS. Such expression platform modulation was not observed by in vitro chemical probing and demonstrates that the crowded cellular environment does not preclude detection of compact and loose RNA-regulatory conformations. Here we describe the design, methods, interpretation, and caveats of Reporter Coupled (ReCo) icSHAPE. We also describe mapping of the differential ReCo-icSHAPE results onto the Lrh riboswitch-preQ1 cocrystal structure. The approach should be readily applicable to functional RNAs triggered by effectors or environmental variations.

Metalloriboswitches: RNA-based inorganic ion sensors that regulate genes.

Divalent ions fulfill essential cellular roles and are required for virulence by certain bacteria. Free intracellular Mg2+ can approach 5 mm, but at this level Mn2+, Ni2+, or Co2+ can be growth-inhibitory, and magnesium fluoride is toxic. To maintain ion homeostasis, many bacteria have evolved ion sensors embedded in the 5'-leader sequences of mRNAs encoding ion uptake or efflux channels. Here, we review current insights into these "metalloriboswitches," emphasizing ion-specific binding by structured RNA aptamers and associated conformational changes in downstream signal sequences. This riboswitch-effector interplay produces a layer of gene regulatory feedback that has elicited interest as an antibacterial target.

Structural analysis of a class III preQ1 riboswitch reveals an aptamer distant from a ribosome-binding site regulated by fast dynamics.

PreQ1-III riboswitches are newly identified RNA elements that control bacterial genes in response to preQ1 (7-aminomethyl-7-deazaguanine), a precursor to the essential hypermodified tRNA base queuosine. Although numerous riboswitches fold as H-type or HLout-type pseudoknots that integrate ligand-binding and regulatory sequences within a single folded domain, the preQ1-III riboswitch aptamer forms a HLout-type pseudoknot that does not appear to incorporate its ribosome-binding site (RBS). To understand how this unusual organization confers function, we determined the crystal structure of the class III preQ1 riboswitch from Faecalibacterium prausnitzii at 2.75 Å resolution. PreQ1 binds tightly (KD,app 6.5 ± 0.5 nM) between helices P1 and P2 of a three-way helical junction wherein the third helix, P4, projects orthogonally from the ligand-binding pocket, exposing its stem-loop to base pair with the 3' RBS. Biochemical analysis, computational modeling, and single-molecule FRET imaging demonstrated that preQ1 enhances P4 reorientation toward P1-P2, promoting a partially nested, H-type pseudoknot in which the RBS undergoes rapid docking (kdock ∼ 0.6 s(-1)) and undocking (kundock ∼ 1.1 s(-1)). Discovery of such dynamic conformational switching provides insight into how a riboswitch with bipartite architecture uses dynamics to modulate expression platform accessibility, thus expanding the known repertoire of gene control strategies used by regulatory RNAs.

Molecular mechanism for preQ1-II riboswitch function revealed by molecular dynamics.

Riboswitches are RNA molecules that regulate gene expression using conformational change, affected by binding of small molecule ligands. A crystal structure of a ligand-bound class II preQ1 riboswitch has been determined in a previous structural study. To gain insight into the dynamics of this riboswitch in solution, eight total molecular dynamic simulations, four with and four without ligand, were performed using the Amber force field. In the presence of ligand, all four of the simulations demonstrated rearranged base pairs at the 3' end, consistent with expected base-pairing from comparative sequence analysis in a prior bioinformatic analysis; this suggests the pairing in this region was altered by crystallization. Additionally, in the absence of ligand, three of the simulations demonstrated similar changes in base-pairing at the ligand binding site. Significantly, although most of the riboswitch architecture remained intact in the respective trajectories, the P3 stem was destabilized in the ligand-free simulations in a way that exposed the Shine-Dalgarno sequence. This work illustrates how destabilization of two major groove base triples can influence a nearby H-type pseudoknot and provides a mechanism for control of gene expression by a fold that is frequently found in bacterial riboswitches.

Structure of a class II preQ1 riboswitch reveals ligand recognition by a new fold.

PreQ1 riboswitches regulate genes by binding the pyrrolopyrimidine intermediate preQ1 during the biosynthesis of the essential tRNA base queuosine. We report what is to our knowledge the first preQ1-II riboswitch structure at 2.3-Å resolution, which uses a previously uncharacterized fold to achieve effector recognition at the confluence of a three-way helical junction flanking a pseudoknotted ribosome-binding site. The results account for translational control mediated by the preQ1-II riboswitch class and expand the known repertoire of ligand-binding modes used by regulatory RNAs.

Single transcriptional and translational preQ1 riboswitches adopt similar pre-folded ensembles that follow distinct folding pathways into the same ligand-bound structure.

Riboswitches are structural elements in the 5' untranslated regions of many bacterial messenger RNAs that regulate gene expression in response to changing metabolite concentrations by inhibition of either transcription or translation initiation. The preQ1 (7-aminomethyl-7-deazaguanine) riboswitch family comprises some of the smallest metabolite sensing RNAs found in nature. Once ligand-bound, the transcriptional Bacillus subtilis and translational Thermoanaerobacter tengcongensis preQ1 riboswitch aptamers are structurally similar RNA pseudoknots; yet, prior structural studies have characterized their ligand-free conformations as largely unfolded and folded, respectively. In contrast, through single molecule observation, we now show that, at near-physiological Mg(2+) concentration and pH, both ligand-free aptamers adopt similar pre-folded state ensembles that differ in their ligand-mediated folding. Structure-based Go-model simulations of the two aptamers suggest that the ligand binds late (Bacillus subtilis) and early (Thermoanaerobacter tengcongensis) relative to pseudoknot folding, leading to the proposal that the principal distinction between the two riboswitches lies in their relative tendencies to fold via mechanisms of conformational selection and induced fit, respectively. These mechanistic insights are put to the test by rationally designing a single nucleotide swap distal from the ligand binding pocket that we find to predictably control the aptamers' pre-folded states and their ligand binding affinities.

Evidence for the role of active site residues in the hairpin ribozyme from molecular simulations along the reaction path.

The hairpin ribozyme accelerates a phosphoryl transfer reaction without catalytic participation of divalent metal ions. Residues A38 and G8 have been implicated as playing roles in general acid and base catalysis, respectively. Here we explore the structure and dynamics of key active site residues using more than 1 µs of molecular dynamics simulations of the hairpin ribozyme at different stages along the catalytic pathway. Analysis of results indicates hydrogen bond interactions between the nucleophile and proR nonbridging oxygen are correlated with active inline attack conformations. Further, the simulation results suggest a possible alternative role for G8 to promote inline fitness and facilitate activation of the nucleophile by hydrogen bonding, although this does not necessarily exclude an additional role as a general base. Finally, we suggest that substitution of G8 with N7- or N3-deazaguanosine which have elevated pKa values, both with and without thio modifications at the 5' leaving group position, would provide valuable insight into the specific role of G8 in catalysis.