RNA dynamics by design: biasing ensembles towards the ligand-bound state

Rapid and accurate determination of atomistic RNA dynamic ensemble models using NMR and structure prediction

Biomolecules form dynamic ensembles of many inter-converting conformations which are key for understanding how they fold and function. However, determining ensembles is challenging because the information required to specify atomic structures for thousands of conformations far exceeds that of experimental measurements. We addressed this data gap and dramatically simplified and accelerated RNA ensemble determination by using structure prediction tools that leverage the growing database of RNA structures to generate a conformation library. Refinement of this library with NMR residual dipolar couplings provided an atomistic ensemble model for HIV-1 TAR, and the model accuracy was independently supported by comparisons to quantum-mechanical calculations of NMR chemical shifts, comparison to a crystal structure of a substate, and through designed ensemble redistribution via atomic mutagenesis. Applications to TAR bulge variants and more complex tertiary RNAs support the generality of this approach and the potential to make the determination of atomic-resolution RNA ensembles routine.

Determining dynamic ensembles of biomolecules is still challenging. Here the authors present an approach for rapid RNA ensemble determination that combines RNA structure prediction tools and NMR residual dipolar coupling data and use it to determine atomistic ensemble models for a variety of RNAs.

Design of a small molecule that stimulates vascular endothelial growth factor A enabled by screening RNA fold-small molecule interactions

Vascular endothelial growth factor A (VEGFA) stimulates angiogenesis in human endothelial cells, and increasing its expression is a potential treatment for heart failure. Here, we report the design of a small molecule (TGP-377) that specifically and potently enhances VEGFA expression by the targeting of a non-coding microRNA that regulates its expression. A selection-based screen, named two-dimensional combinatorial screening, revealed preferences in small-molecule chemotypes that bind RNA and preferences in the RNA motifs that bind small molecules. The screening program increased the dataset of known RNA motif–small molecule binding partners by 20-fold. Analysis of this dataset against the RNA-mediated pathways that regulate VEGFA defined that the microRNA-377 precursor, which represses Vegfa messenger RNA translation, is druggable in a selective manner. We designed TGP-377 to potently and specifically upregulate VEGFA in human umbilical vein endothelial cells. These studies illustrate the power of two-dimensional combinatorial screening to define molecular recognition events between ‘undruggable’ biomolecules and small molecules, and the ability of sequence-based design to deliver efficacious structure-specific compounds.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

Demonstration that Small Molecules can Bind and Stabilize Low-abundance Short-lived RNA Excited Conformational States

Many promising RNA drug targets have functions that require the formation of RNA-protein complexes, but inhibiting RNA-protein interactions can prove difficult using small molecules. Regulatory RNAs have been shown to transiently form excited conformational states (ESs) that remodel local aspects of secondary structure. In some cases, the ES conformation has been shown to be inactive and to be poorly recognized by protein binding partners. In these cases, specifically targeting and stabilizing the RNA ES using a small molecule provides a rational structure-based strategy for inhibiting RNA activity. However, this requires that a small molecule discriminates between two conformations of the same RNA to preferentially bind and stabilize the short-lived low-abundance ES relative to the long-lived more abundant ground state (GS). Here, we tested the feasibility of this approach by designing a mutant that inverts the conformational equilibrium of the HIV-1 transactivation response element (TAR) RNA, such that the native GS conformation becomes a low-abundance ES. Using this mutant and NMR chemical shift mapping experiments, we show that argininamide, a ligand mimic of TAR's cognate protein binding partner Tat, is able to restore a native-like conformation by preferentially binding and stabilizing the transient and low-populated ES. A synthetic small molecule optimized to bind the TAR GS also partially stabilized the ES, whereas an aminoglycoside molecule that binds RNAs nonspecifically did not preferentially stabilize the ES to a similar extent. These results support the feasibility of inhibiting RNA activity using small molecules that preferentially bind and stabilize the ES.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Single-Molecule FRET Kinetics of the Mn2+ Riboswitch: Evidence for Allosteric Mg2+ Control of "Induced-Fit" vs "Conformational Selection" Folding Pathways

Gene expression in bacteria is often regulated dynamically by conformational changes in a riboswitch upon ligand binding, a detailed understanding of which is very much in its infancy. For example, the manganese riboswitch is a widespread RNA motif that conformationally responds in regulating bacterial gene expression to micromolar levels of its eponymous ligand, Mn²⁺, but the mechanistic pathways are poorly understood. In this work, we quantitatively explore the dynamic folding behavior of the manganese riboswitch by single-molecule fluorescence resonance energy transfer spectroscopy as a function of cation/ligand conditions. From the detailed analysis of the kinetics, the Mn²⁺ is shown to fold the riboswitch by a "bind-then-fold" (i.e., "induced-fit", IF) mechanism, whereby the ligand binds first and then promotes folding. On the other hand, the data also clearly reveal the presence of a folded yet ligand-free structure predominating due to the addition of physiological Mg²⁺ to a nonselective metal ion binding site. Of particular kinetic interest, such a Mg²⁺ "prefolded" conformation of the riboswitch is shown to exhibit a significantly increased affinity for Mn²⁺ and further stabilization by subsequent binding of the ligand, thereby promoting efficient riboswitch folding by a "fold-then-bind" (i.e., "conformational selection", CS) mechanism. Our results not only demonstrate Mg²⁺-controlled switching between IF and CS riboswitch folding pathways but also suggest a novel heterotropic allosteric control in the manganese riboswitch activity co-regulated by Mg²⁺ binding.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

Approaches to Validate and Manipulate RNA Targets with Small Molecules in Cells

RNA has become an increasingly important target for therapeutic interventions and for chemical probes that dissect and manipulate its cellular function. Emerging targets include human RNAs that have been shown to directly cause cancer, metabolic disorders, and genetic disease. In this review, we describe various routes to obtain bioactive compounds that target RNA, with a particular emphasis on the development of small molecules. We use these cases to describe approaches that are being developed for target validation, which include target-directed cleavage, classic pull-down experiments, and covalent cross-linking. Thus, tools are available to design small molecules to target RNA and to identify the cellular RNAs that are their targets.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

A Highly Coupled Network of Tertiary Interactions in the SAM-I Riboswitch and Their Role in Regulatory Tuning

RNA folding in vivo is significantly influenced by transcription, which is not necessarily recapitulated by Mg(2+)-induced folding of the corresponding full-length RNA in vitro. Riboswitches that regulate gene expression at the transcriptional level are an ideal system for investigating this aspect of RNA folding as ligand-dependent termination is obligatorily co-transcriptional, providing a clear readout of the folding outcome. The folding of representative members of the SAM-I family of riboswitches has been extensively analyzed using approaches focusing almost exclusively upon Mg(2+) and/or S-adenosylmethionine (SAM)-induced folding of full-length transcripts of the ligand binding domain. To relate these findings to co-transcriptional regulatory activity, we have investigated a set of structure-guided mutations of conserved tertiary architectural elements of the ligand binding domain using an in vitro single-turnover transcriptional termination assay, complemented with phylogenetic analysis and isothermal titration calorimetry data. This analysis revealed a conserved internal loop adjacent to the SAM binding site that significantly affects ligand binding and regulatory activity. Conversely, most single point mutations throughout key conserved features in peripheral tertiary architecture supporting the SAM binding pocket have relatively little impact on riboswitch activity. Instead, a secondary structural element in the peripheral subdomain appears to be the key determinant in observed differences in regulatory properties across the SAM-I family. These data reveal a highly coupled network of tertiary interactions that promote high-fidelity co-transcriptional folding of the riboswitch but are only indirectly linked to regulatory tuning.

Unveiling Inherent Degeneracies in Determining Population-Weighted Ensembles of Interdomain Orientational Distributions Using NMR Residual Dipolar Couplings: Application to RNA Helix Junction Helix Motifs

A growing number of studies employ time-averaged experimental data to determine dynamic ensembles of biomolecules. While it is well-known that different ensembles can satisfy experimental data to within error, the extent and nature of these degeneracies, and their impact on the accuracy of the ensemble determination remains poorly understood. Here, we use simulations and a recently introduced metric for assessing ensemble similarity to explore degeneracies in determining ensembles using NMR residual dipolar couplings (RDCs) with specific application to A-form helices in RNA. Various target ensembles were constructed representing different domain-domain orientational distributions that are confined to a topologically restricted (<10%) conformational space. Five independent sets of ensemble averaged RDCs were then computed for each target ensemble and a "sample and select" scheme used to identify degenerate ensembles that satisfy RDCs to within experimental uncertainty. We find that ensembles with different ensemble sizes and that can differ significantly from the target ensemble (by as much as ∑Ω ∼ 0.4 where ∑Ω varies between 0 and 1 for maximum and minimum ensemble similarity, respectively) can satisfy the ensemble averaged RDCs. These deviations increase with the number of unique conformers and breadth of the target distribution, and result in significant uncertainty in determining conformational entropy (as large as 5 kcal/mol at T = 298 K). Nevertheless, the RDC-degenerate ensembles are biased toward populated regions of the target ensemble, and capture other essential features of the distribution, including the shape. Our results identify ensemble size as a major source of uncertainty in determining ensembles and suggest that NMR interactions such as RDCs and spin relaxation, on their own, do not carry the necessary information needed to determine conformational entropy at a useful level of precision. The framework introduced here provides a general approach for exploring degeneracies in ensemble determination for different types of experimental data.

Slowdown of Interhelical Motions Induces a Glass Transition in RNA

RNA function depends crucially on the details of its dynamics. The simplest RNA dynamical unit is a two-way interhelical junction. Here, for such a unit--the transactivation response RNA element--we present evidence from molecular dynamics simulations, supported by nuclear magnetic resonance relaxation experiments, for a dynamical transition near 230 K. This glass transition arises from the freezing out of collective interhelical motional modes. The motions, resolved with site-specificity, are dynamically heterogeneous and exhibit non-Arrhenius relaxation. The microscopic origin of the glass transition is a low-dimensional, slow manifold consisting largely of the Euler angles describing interhelical reorientation. Principal component analysis over a range of temperatures covering the glass transition shows that the abrupt slowdown of motion finds its explanation in a localization transition that traps probability density into several disconnected conformational pools over the low-dimensional energy landscape. Upon temperature increase, the probability density pools then flood a larger basin, akin to a lakes-to-sea transition. Simulations on transactivation response RNA are also used to backcalculate inelastic neutron scattering data that match previous inelastic neutron scattering measurements on larger and more complex RNA structures and which, upon normalization, give temperature-dependent fluctuation profiles that overlap onto a glass transition curve that is quasi-universal over a range of systems and techniques.

Small-angle X-ray scattering: a bridge between RNA secondary structures and three-dimensional topological structures

Whereas the structures of small to medium-sized well folded RNA molecules often can be determined by either X-ray crystallography or NMR spectroscopy, obtaining structural information for large RNAs using experimental, computational, or combined approaches remains a major interest and challenge. RNA is very sensitive to small-angle X-ray scattering (SAXS) due to high electron density along phosphate-sugar backbones, whose scattering contribution dominates SAXS intensity. For this reason, SAXS is particularly useful in obtaining global RNA structural information that outlines backbone topologies and, therefore, molecular envelopes. Such information is extremely valuable in bridging the gap between the secondary structures and three-dimensional topological structures of RNA molecules, particularly those that have proven difficult to study using other structure-determination methods. Here we review published results of RNA topological structures derived from SAXS data or in combination with other experimental data, as well as details on RNA sample preparation for SAXS experiments.

Dynamic motions of the HIV-1 frameshift site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDOint) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDOint = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.

The design of RNA binders: targeting the HIV replication cycle as a case study

The human immunodeficiency virus 1 (HIV-1) replication cycle is finely tuned with many important steps involving RNA-RNA or protein-RNA interactions, all of them being potential targets for the development of new antiviral compounds. This cycle can also be considered as a good benchmark for the evaluation of early-stage strategies aiming at designing drugs that bind to RNA, with the possibility to correlate in vitro activities with antiviral properties. In this review, we highlight different approaches developed to interfere with four important steps of the HIV-1 replication cycle: the early stage of reverse transcription, the transactivation of viral transcription, the nuclear export of partially spliced transcripts and the dimerization step.

Hierarchy of RNA functional dynamics

RNA dynamics play a fundamental role in many cellular functions. However, there is no general framework to describe these complex processes, which typically consist of many structural maneuvers that occur over timescales ranging from picoseconds to seconds. Here, we classify RNA dynamics into distinct modes representing transitions between basins on a hierarchical free-energy landscape. These transitions include large-scale secondary-structural transitions at >0.1-s timescales, base-pair/tertiary dynamics at microsecond-to-millisecond timescales, stacking dynamics at timescales ranging from nanoseconds to microseconds, and other "jittering" motions at timescales ranging from picoseconds to nanoseconds. We review various modes within these three different tiers, the different mechanisms by which they are used to regulate function, and how they can be coupled together to achieve greater functional complexity.

Coarse grained models reveal essential contributions of topological constraints to the conformational free energy of RNA bulges

Recent studies have shown that simple stereochemical constraints encoded at the RNA secondary structure level significantly restrict the orientation of RNA helices across two-way junctions and yield physically reasonable distributions of RNA 3D conformations. Here we develop a new coarse-grain model, TOPRNA, that is optimized for exploring detailed aspects of these topological constraints in complex RNA systems. Unlike prior models, TOPRNA effectively treats RNAs as collections of semirigid helices linked by freely rotatable single strands, allowing us to isolate the effects of secondary structure connectivity and sterics on 3D structure. Simulations of bulge junctions show that TOPRNA captures new aspects of topological constraints, including variations arising from deviations in local A-form structure, translational displacements of the helices, and stereochemical constraints imposed by bulge-linker nucleotides. Notably, these aspects of topological constraints define free energy landscapes that coincide with the distribution of bulge conformations in the PDB. Our simulations also quantitatively reproduce NMR RDC measurements made on HIV-1 TAR at low salt concentrations, although not for different TAR mutants or at high salt concentrations. Our results confirm that topological constraints are an important determinant of bulge conformation and dynamics and demonstrate the utility of TOPRNA for studying the topological constraints of complex RNAs.

Methods to enable the design of bioactive small molecules targeting RNA

RNA is an immensely important target for small molecule therapeutics or chemical probes of function. However, methods that identify, annotate, and optimize RNA-small molecule interactions that could enable the design of compounds that modulate RNA function are in their infancies. This review describes recent approaches that have been developed to understand and optimize RNA motif-small molecule interactions, including structure-activity relationships through sequencing (StARTS), quantitative structure-activity relationships (QSAR), chemical similarity searching, structure-based design and docking, and molecular dynamics (MD) simulations. Case studies described include the design of small molecules targeting RNA expansions, the bacterial A-site, viral RNAs, and telomerase RNA. These approaches can be combined to afford a synergistic method to exploit the myriad of RNA targets in the transcriptome.

Distinct conformational transition patterns of noncoding 7SK snRNA and HIV TAR RNAs upon Tat binding

Noncoding 7SK snRNA is believed to play an important role in the recruitment of P-TEFb by viral protein Tat to stimulate HIV processive transcription. Because HIV-2 TAR RNA and 7SK both evolved to feature a dinucleotide bulge region, compared to the trinucleotide bulge for HIV-1 TAR, ultrafast time-resolved fluorescence spectroscopy has been used to probe the conformational landscape of HIV-2 TAR and 7SK-SL4 RNA to monitor the conformational changes upon Tat binding. Our studies demonstrate that both HIV-1/2 TAR and 7SK-SL4 sample heterogeneous ensembles in the free state and undergo distinct conformational transitions upon Tat binding. These findings provide exquisite knowledge on the conformational complexity and intricate mechanism of molecular recognition and pave the way for drug design and discovery that incorporate dynamics information.

Molecular mechanism of GTPase activation at the signal recognition particle (SRP) RNA distal end

The signal recognition particle (SRP) RNA is a universally conserved and essential component of the SRP that mediates the co-translational targeting of proteins to the correct cellular membrane. During the targeting reaction, two functional ends in the SRP RNA mediate distinct functions. Whereas the RNA tetraloop facilitates initial assembly of two GTPases between the SRP and SRP receptor, this GTPase complex subsequently relocalizes ∼100 Å to the 5',3'-distal end of the RNA, a conformation crucial for GTPase activation and cargo handover. Here we combined biochemical, single molecule, and NMR studies to investigate the molecular mechanism of this large scale conformational change. We show that two independent sites contribute to the interaction of the GTPase complex with the SRP RNA distal end. Loop E plays a crucial role in the precise positioning of the GTPase complex on these two sites by inducing a defined bend in the RNA helix and thus generating a preorganized recognition surface. GTPase docking can be uncoupled from its subsequent activation, which is mediated by conserved bases in the next internal loop. These results, combined with recent structural work, elucidate how the SRP RNA induces GTPase relocalization and activation at the end of the protein targeting reaction.

Influence of dimethylsulfoxide on RNA structure and ligand binding

Dimethyl sulfoxide (DMSO) is widely used as a cosolvent to solubilize hydrophobic compounds in RNA-ligand binding assays. Although it is known that high concentrations of DMSO (>75%) can significantly affect RNA structure and folding energetics, a thorough analysis of how lower concentrations (<10%) of DMSO typically used in binding assays affects RNA structure and ligand binding has not been undertaken. Here, we use NMR and 2-aminopurine fluorescence spectroscopy to examine how DMSO affects the structure, dynamics, and ligand binding properties of two flexible hairpin RNAs: the transactivation response element from HIV-1 and bacterial ribosomal A-site. In both cases, 5-10% DMSO decreased stacking interactions and increased local disorder in noncanonical residues within bulges and loops and resulted in 0.3-4-fold reduction in the measured binding affinities for different small molecules, with the greatest reduction observed for an intercalating compound that binds RNA nonspecifically. Our results suggest that, by competing for hydrophobic interactions, DMSO can have a small but significant effect on RNA structure and ligand binding. These effects should be considered when developing ligand binding assays and high throughput screens.

A general method for constructing atomic-resolution RNA ensembles using NMR residual dipolar couplings: the basis for interhelical motions revealed

The ability to modulate alignment and measure multiple independent sets of NMR residual dipolar couplings (RDCs) has made it possible to characterize internal motions in proteins at atomic resolution and with time scale sensitivity ranging from picoseconds up to milliseconds. The application of such methods to the study of RNA dynamics, however, remains fundamentally limited by the inability to modulate alignment and by strong couplings between internal and overall motions that complicate the quantitative interpretation of RDCs. Here, we address this problem by showing that RNA alignment can be generally modulated, in a controlled manner, by variable elongation of A-form helices and that the information contained within the measured RDCs can be extracted even in the presence of strong couplings between motions and overall alignment via structure-based prediction of alignment. Using this approach, four RDC data sets, and a broad conformational pool obtained from a 8.2 μs molecular dynamics simulation, we successfully construct and validate an atomic resolution ensemble of human immunodeficiency virus type I transactivation response element RNA. This ensemble reveals local motions in and around the bulge involving changes in stacking and hydrogen-bonding interactions, which are undetectable by traditional spin relaxation and drive global changes in interhelical orientation. This new approach broadens the scope of using RDCs in characterizing the dynamics of nucleic acids.

Nucleotides adjacent to the ligand-binding pocket are linked to activity tuning in the purine riboswitch

Direct sensing of intracellular metabolite concentrations by riboswitch RNAs provides an economical and rapid means to maintain metabolic homeostasis. Since many organisms employ the same class of riboswitch to control different genes or transcription units, it is likely that functional variation exists in riboswitches such that activity is tuned to meet cellular needs. Using a bioinformatic approach, we have identified a region of the purine riboswitch aptamer domain that displays conservation patterns linked to riboswitch activity. Aptamer domain compositions within this region can be divided into nine classes that display a spectrum of activities. Naturally occurring compositions in this region favor rapid association rate constants and slow dissociation rate constants for ligand binding. Using X-ray crystallography and chemical probing, we demonstrate that both the free and bound states are influenced by the composition of this region and that modest sequence alterations have a dramatic impact on activity. The introduction of non-natural compositions result in the inability to regulate gene expression in vivo, suggesting that aptamer domain activity is highly plastic and thus readily tunable to meet cellular needs.

Dominant role of the 5' TAR bulge in dimerization of HIV-1 genomic RNA, but no evidence of TAR-TAR kissing during in vivo virus assembly

The 5' untranslated region of HIV-1 genomic RNA (gRNA) contains two stem-loop structures that appear to be equally important for gRNA dimerization: the 57-nucleotide 5' TAR, at the very 5' end, and the 35-nucleotide SL1 (nucleotides 243-277). SL1 is well-known for containing the dimerization initiation site (DIS) in its apical loop. The DIS is a six-nucleotide palindrome. Here, we investigated the mechanism of TAR-directed gRNA dimerization. We found that the trinucleotide bulge (UCU24) of the 5' TAR has dominant impacts on both formation of HIV-1 RNA dimers and maturation of the formed dimers. The ΔUCU trinucleotide deletion strongly inhibited the first process and blocked the other, thus impairing gRNA dimerization as severely as deletion of the entire 5' TAR, and more severely than deletion of the DIS, inactivation of the viral protease, or most severe mutations in the nucleocapsid protein. The apical loop of TAR contains a 10-nucleotide palindrome that has been postulated to stimulate gRNA dimerization by a TAR-TAR kissing mechanism analogous to the one used by SL1 to stimulate dimerization. Using mutations that strongly destabilize formation of the TAR palindrome duplex, as well as compensatory mutations that restore duplex formation to a wild-type-like level, we found no evidence of TAR-TAR kissing, even though mutations nullifying the kissing potential of the TAR palindrome could impair dimerization by a mechanism other than hindering of SL1. However, nullifying the kissing potential of TAR had much less severe effects than ΔUCU. By not uncovering a dimerization mechanism intrinsic to TAR, our data suggest that TAR mutations exert their effect 3' of TAR, yet not on SL1, because TAR and SL1 mutations have synergistic effects on gRNA dimerization.

Functional complexity and regulation through RNA dynamics

Changes to the conformation of coding and non-coding RNAs form the basis of elements of genetic regulation and provide an important source of complexity, which drives many of the fundamental processes of life. Although the structure of RNA is highly flexible, the underlying dynamics of RNA are robust and are limited to transitions between the few conformations that preserve favourable base-pairing and stacking interactions. The mechanisms by which cellular processes harness the intrinsic dynamic behaviour of RNA and use it within functionally productive pathways are complex. The versatile functions and ease by which it is integrated into a wide variety of genetic circuits and biochemical pathways suggests there is a general and fundamental role for RNA dynamics in cellular processes.

The dynamic nature of RNA as key to understanding riboswitch mechanisms

Riboswitches are gene regulation elements within RNA that recognize specific metabolites. They predominantly occur in the untranslated leader regions of bacterial messenger RNA (mRNA). Upon metabolite binding to the aptamer domain, a structural change in the adjoining downstream expression platform signals "on" or "off" for gene expression. Researchers have achieved much progress in characterizing ligand-bound riboswitch states at the molecular level; an impressive number of high-resolution structures of aptamer-ligand complexes is now available. These structures have significantly contributed toward our understanding of how riboswitches interact with their natural ligands and with structurally related analogues. In contrast, relatively little is known about the nature of the unbound (apo) form of riboswitches. Moreover, the details of how changes in the aptamer domain are transduced into conformational changes in the decision-making expression platform remain murky. In this Account, we report on recent efforts aimed at the characterization of free states, ligand recognition, and ligand-induced folding in riboswitches. Riboswitch action is best approached as a cotranscriptional process, which implies sequential folding and release of the aptamer prior to the signaling of the expression platform. Thus, a complex interplay of several factors has to be taken into account, such as speed of transcription, transcriptional pausing, kinetics and thermodynamics of RNA structure formation, and kinetics and thermodynamics of ligand binding. The response mechanism appears to be best described as a process in which ligand recognition critically dictates the folding pathway of the nascent mRNA during its expression; the resulting structures determine the interactions with the transcriptional or translational apparatus. We discuss experimental methods that offer insight into the dynamics of the free riboswitch state. These include probing experiments, such as in-line and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) techniques, small-angle X-ray scattering (SAXS) analysis, NMR spectroscopy, and fluorescence spectroscopy, including single-molecule fluorescence resonance energy transfer (smFRET) imaging. One of our research contributions is an approach, termed 2ApFold, that incorporates noninvasive 2-aminopurine modifications in riboswitches. The fluorescence response of these moieties is used to delineate the order of secondary-tertiary structure formation and rearrangements taking place during ligand-induced folding. This information can be used to explore the kinetics of ligand recognition and to analyze the degree of structure preorganization of the free riboswitch state. Furthermore, we discuss a recent smFRET study on the SAM-II riboswitch; this report underscores the importance of choosing strategic labeling patterns that leave maximal conformational freedom to the regulatory interaction. Finally, we comment on how riboswitch ligand recognition appeals to the concepts of conformational selection and induced fit, and on the question of whether riboswitches act under thermodynamic or kinetic control. This Account highlights the fact that a thorough understanding of RNA dynamics in vitro is required to shed light on cellular riboswitch mechanisms. Elucidating these mechanisms will contribute not only to ongoing efforts to target riboswitches with antibiotics but also to attempts to engineer artificial cell regulation systems.

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.

Myotonic dystrophy type 1 RNA crystal structures reveal heterogeneous 1 × 1 nucleotide UU internal loop conformations

RNA internal loops often display a variety of conformations in solution. Herein, we visualize conformational heterogeneity in the context of the 5'CUG/3'GUC repeat motif present in the RNA that causes myotonic dystrophy type 1 (DM1). Specifically, two crystal structures of a model DM1 triplet repeating construct, 5'r[UUGGGC(CUG)(3)GUCC](2), refined to 2.20 and 1.52 Å resolution are disclosed. Here, differences in the orientation of the 5' dangling UU end between the two structures induce changes in the backbone groove width, which reveals that noncanonical 1 × 1 nucleotide UU internal loops can display an ensemble of pairing conformations. In the 2.20 Å structure, CUGa, the 5' UU forms a one hydrogen-bonded pair with a 5' UU of a neighboring helix in the unit cell to form a pseudoinfinite helix. The central 1 × 1 nucleotide UU internal loop has no hydrogen bonds, while the terminal 1 × 1 nucleotide UU internal loops each form a one-hydrogen bond pair. In the 1.52 Å structure, CUGb, the 5' UU dangling end is tucked into the major groove of the duplex. While the canonically paired bases show no change in base pairing, in CUGb the terminal 1 × 1 nucleotide UU internal loops now form two hydrogen-bonded pairs. Thus, the shift in the major groove induced by the 5' UU dangling end alters noncanonical base patterns. Collectively, these structures indicate that 1 × 1 nucleotide UU internal loops in DM1 may sample multiple conformations in vivo. This observation has implications for the recognition of this RNA, and other repeating transcripts, by protein and small molecule ligands.

Discovery of selective bioactive small molecules by targeting an RNA dynamic ensemble

Current approaches used to identify protein-binding small molecules are not suited for identifying small molecules that can bind emerging RNA drug targets. By docking small molecules onto an RNA dynamic ensemble constructed by combining NMR spectroscopy and computational molecular dynamics, we virtually screened small molecules that target the entire structure landscape of the transactivation response element (TAR) from HIV type 1 (HIV-1). We quantitatively predict binding energies for small molecules that bind different RNA conformations and report the de novo discovery of six compounds that bind TAR with high affinity and inhibit its interaction with a Tat peptide in vitro (K(i) values of 710 nM-169 μM). One compound binds HIV-1 TAR with marked selectivity and inhibits Tat-mediated activation of the HIV-1 long terminal repeat by 81% in T-cell lines and HIV replication in an HIV-1 indicator cell line (IC(50) ∼23.1 μM).

Use of RNA structure flexibility data in nanostructure modeling

In the emerging field of RNA-based nanotechnology there is a need for automation of the structure design process. Our goal is to develop computer methods for aiding in this process. Towards that end, we created the RNA junction database, which is a repository of RNA junctions, i.e. internal, multi-branch and kissing loops with emanating stem stubs, extracted from the larger RNA structures stored in the PDB database. These junctions can be used as building blocks for nanostructures. Two programs developed in our laboratory, NanoTiler and RNA2D3D, can combine such building blocks with idealized fragments of A-form helices to produce desired 3D nanostructures. Initially, the building blocks are treated as rigid objects and the resulting geometry is tested against the design objectives. Experimental data, however, shows that RNA accommodates its shape to the constraints of larger structural contexts. Therefore we are adding analysis of the flexibility of our building blocks to the full design process. Here we present an example of RNA-based nanostructure design, putting emphasis on the need to characterize the structural flexibility of the building blocks to induce ring closure in the automated exploration. We focus on the use of kissing loops (KL) in nanostructure design, since they have been shown to play an important role in RNA self-assembly. By using an experimentally proven system, the RNA tectosquare, we show that considering the flexibility of the KLs as well as distortions of helical regions may be necessary to achieve a realistic design.

Dynamic ensemble view of the conformational landscape of HIV-1 TAR RNA and allosteric recognition

RNA conformational dynamics and the resulting structural heterogeneity play an important role in RNA functions, e.g., recognition. Recognition of HIV-1 TAR RNA has been proposed to occur via a conformational capture mechanism. Here, using ultrafast time-resolved fluorescence spectroscopy, we have probed the complexity of the conformational landscape of HIV-1 TAR RNA and monitored the position-dependent changes in the landscape upon binding of a Tat protein-derived peptide and neomycin B. In the ligand-free state, the TAR RNA samples multiple families of conformations with various degrees of base stacking around the three-nucleotide bulge region. Some subpopulations partially resemble those ligand-bound states, but the coaxially stacked state is below the detection limit. When Tat or neomycin B binds, the bulge region as an ensemble undergoes a conformational transition in a position-dependent manner. Tat and neomycin B induce mutually exclusive changes in the TAR RNA underlying the mechanism of allosteric inhibition at an ensemble level with residue-specific details. Time-resolved anisotropy decay measurements revealed picosecond motions of bases in both ligand-free and ligand-bound states. Mutation of a base pair at the bulge--stem junction has differential effects on the conformational distributions of the bulge bases. A dynamic model of the ensemble view of the conformational landscape for HIV-1 TAR RNA is proposed, and the implication of the general mechanism of RNA recognition and its impact on RNA-based therapeutics are discussed.

Conformational capture of the SAM-II riboswitch

Riboswitches are gene regulation elements in mRNA that function by specifically responding to metabolites. Although the metabolite-bound states of riboswitches have proven amenable to structure determination efforts, knowledge of the structural features of riboswitches in their ligand-free forms and their ligand-response mechanisms giving rise to regulatory control is lacking. Here we explore the ligand-induced folding process of the S-adenosylmethionine type II (SAM-II) riboswitch using chemical and biophysical methods, including NMR and fluorescence spectroscopy, and single-molecule fluorescence imaging. The data reveal that the unliganded SAM-II riboswitch is dynamic in nature, in that its stem-loop element becomes engaged in a pseudoknot fold through base-pairing with nucleosides in the 3' overhang containing the Shine-Dalgarno sequence. Although the pseudoknot structure is highly transient in the absence of its ligand, S-adenosylmethionine (SAM), it becomes conformationally restrained upon ligand recognition, through a conformational capture mechanism. These insights provide a molecular understanding of riboswitch dynamics that shed new light on the mechanism of riboswitch-mediated translational regulation.

Topological constraints: using RNA secondary structure to model 3D conformation, folding pathways, and dynamic adaptation

Accompanying recent advances in determining RNA secondary structure is the growing appreciation for the importance of relatively simple topological constraints, encoded at the secondary structure level, in defining the overall architecture, folding pathways, and dynamic adaptability of RNA. A new view is emerging in which tertiary interactions do not define RNA 3D structure, but rather, help select specific conformers from an already narrow, topologically pre-defined conformational distribution. Studies are providing fundamental insights into the nature of these topological constraints, how they are encoded by the RNA secondary structure, and how they interplay with other interactions, breathing new meaning to RNA secondary structure. New approaches have been developed that take advantage of topological constraints in determining RNA backbone conformation based on secondary structure, and a limited set of other, easily accessible constraints. Topological constraints are also providing a much-needed framework for rationalizing and describing RNA dynamics and structural adaptation. Finally, studies suggest that topological constraints may play important roles in steering RNA folding pathways. Here, we review recent advances in our understanding of topological constraints encoded by the RNA secondary structure.

Recognition of S-adenosylmethionine by riboswitches

Riboswitches are regulatory elements commonly found in the 5' leader sequences of bacterial mRNAs that bind cellular metabolites to direct expression at either the transcriptional or translational level. The effectors of these RNAs are chemically diverse, including nucleobases and nucleosides, amino acids, cofactors, and second messenger molecules. Over the last few years, a number of structures have revealed the architectural means by which RNA creates binding pockets of high affinity and specificity for these compounds. For most effectors, there is a single class of associated riboswitches. However, eight individual classes of S-adenosylmethionine (SAM) and/or S-adenosylhomocysteine (SAH) responsive riboswitches that control various aspects of sulfur metabolism have been validated, revealing a diverse set of solutions to the recognition of these ubiquitous metabolites. This review focuses upon the structures of RNAs that bind SAM and SAH and how they discriminate between these compounds.