Search

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

Comprehensive understanding of the structure–function relationship of RNA both in real time and at atomic level will have a profound impact in advancing our understanding of RNA functions in biology. Here, we describe the first example of a multifunctional nucleoside probe, containing a conformation‐sensitive fluorophore and an anomalous X‐ray diffraction label (5‐selenophene uracil), which enables the correlation of RNA conformation and recognition under equilibrium and in 3D. The probe incorporated into the bacterial ribosomal RNA decoding site, fluorescently reports antibiotic binding and provides diffraction information in determining the structure without distorting native RNA fold. Further, by comparing solution binding data and crystal structure, we gained insight on how the probe senses ligand‐induced conformational change in RNA. Taken together, our nucleoside probe represents a new class of biophysical tool that would complement available tools for functional RNA investigations.

Three-dimensional RNA structures are notoriously difficult to determine, and the link between secondary structure and RNA conformation is only beginning to be understood. These challenges have hindered the identification of guiding principles for small molecule:RNA recognition. We herein demonstrate that the strong and differential binding ability of aminoglycosides to RNA structures can be used to classify five canonical RNA secondary structure motifs through principal component analysis (PCA). In these analyses, the aminoglycosides act as receptors, while RNA structures labeled with a benzofuranyluridine fluorophore act as analytes. Complete (100%) predictive ability for this RNA training set was achieved by incorporating two exhaustively guanidinylated aminoglycosides into the receptor library. The PCA was then externally validated using biologically relevant RNA constructs. In bulge-stem-loop constructs of HIV-1 transactivation response element (TAR) RNA, we achieved nucleotide-specific classification of two independent secondary structure motifs. Furthermore, examination of cheminformatic parameters and PCA loading factors revealed trends in aminoglycoside:RNA recognition, including the importance of shape-based discrimination, and suggested the potential for size and sequence discrimination within RNA structural motifs. These studies present a new approach to classifying RNA structure and provide direct evidence that RNA topology, in addition to sequence, is critical for the molecular recognition of RNA.

Interrogating non-coding RNA structures and functions with small molecules is an area of rapidly increasing interest among biochemists and chemical biologists. However, many biochemical approaches to monitoring RNA structures are time-consuming and low-throughput, and thereby are only of limited utility for RNA-small molecule studies. Fluorescence-based techniques are powerful tools for rapid investigation of RNA conformations, dynamics, and interactions with small molecules. Many fluorescence methods are amenable to high-throughput analysis, enabling library screening for small molecule binders. In this review, we summarize numerous fluorescence-based approaches for identifying and characterizing RNA-small molecule interactions. We describe in detail a high-information content dual-reporter FRET assay we developed to characterize small molecule-induced conformational and stability changes. Our assay is uniquely suited as a platform for both small molecule discovery and thorough characterization of RNA-small molecule binding mechanisms. Given the growing recognition of non-coding RNAs as attractive targets for therapeutic intervention, we anticipate our FRET assay and other fluorescence-based techniques will be indispensable for the development of potent and specific small molecule inhibitors targeting RNA.

Deciphering the conformations of RNAs in their cellular environment allows identification of RNA elements with potentially functional roles within biological contexts. Insight into the conformation of RNA in cells has been achieved using chemical probes that were developed to react specifically with flexible RNA nucleotides, or the Watson–Crick face of single-stranded nucleotides. The most widely used probes are either selective SHAPE (2′-hydroxyl acylation and primer extension) reagents that probe nucleotide flexibility, or dimethyl sulfate (DMS), which probes the base-pairing at adenine and cytosine but is unable to interrogate guanine or uracil. The constitutively charged carbodiimide N-cyclohexyl-N′-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate (CMC) is widely used for probing G and U nucleotides, but has not been established for probing RNA in cells. Here, we report the use of a smaller and conditionally charged reagent, 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC), as a chemical probe of RNA conformation, and the first reagent validated for structure probing of unpaired G and U nucleotides in intact cells. We showed that EDC demonstrates similar reactivity to CMC when probing transcripts in vitro. We found that EDC specifically reacted with accessible nucleotides in the 7SK noncoding RNA in intact cells. We probed structured regions within the Xist lncRNA with EDC and integrated these data with DMS probing data. Together, EDC and DMS allowed us to refine predicted structure models for the 3′ extension of repeat C within Xist. These results highlight how complementing DMS probing experiments with EDC allows the analysis of Watson–Crick base-pairing at all four nucleotides of RNAs in their cellular context.

RNA structure is essential to understand RNA function and regulation in cellular processes. RNA acts either in isolation or as part of a complex with proteins. Both isolated and protein-complexed RNA represent a challenge for structural biology, due to the complexity of its conformational space and intrinsic dynamics. NMR, with its capability to cope with dynamic structures, is the optimal technique to study RNA conformation, ideally in solution. However, RNA–protein complexes are often very large, which limits the application of solution-state NMR. In this chapter we describe the methodology that we developed to determine the structure of RNA by solid-state NMR. Solid-state NMR is not limited by large molecular weights and can be optimally applied to study both the RNA and the protein components of large RNA–protein complexes. We review the methods for resonance assignments, collection of RNA-specific distance restraints, as well as detection of protein–RNA interfaces.

3' downstream inhibitory stem plays a crucial role in locking rpoS mRNA 5' untranslated region in a self-inhibitory state. Here, we used optical tweezers to study the unfolding/refolding of rpoS inhibitory stem in the absence and presence of Mg2+. We found adding Mg2+ decreased the free energy of the RNA junction without re-arranging its secondary structure, through confirming that this RNA formed a canonical RNA three-way junction. We suspected increased free energy might change the relative orientation of different stems of rpoS and confirmed this by small angle X-ray scattering. Such changed conformation may improve Hfq-bridged annealing between sRNA and rpoS RNA inhibitory stem. We established a convenient route to analyze the changes of RNA conformation and folding dynamics by combining optical tweezers with X-ray scattering methods. This route can be easily applied in the studies of other RNA structure and ligand-RNA.

The fate of cellular RNAs is largely dependent on their structural conformation, which determines the assembly of ribonucleoprotein (RNP) complexes. Consequently, RNA-binding proteins (RBPs) play a pivotal role in the lifespan of RNAs. The advent of highly sensitive in cellulo approaches for studying RNPs reveals the presence of unprecedented RNA-binding domains (RBDs). Likewise, the diversity of the RNA targets associated with a given RBP increases the code of RNA-protein interactions. Increasing evidence highlights the biological relevance of RNA conformation for recognition by specific RBPs and how this mutual interaction affects translation control. In particular, noncanonical RBDs present in proteins such as Gemin5, Roquin-1, Staufen, and eIF3 eventually determine translation of selective targets. Collectively, recent studies on RBPs interacting with RNA in a structure-dependent manner unveil new pathways for gene expression regulation, reinforcing the pivotal role of RNP complexes in genome decoding.

In diverse bacterial species, the global regulator Hfq contributes to post-transcriptional networks that control expression of numerous genes. Hfq of the opportunistic pathogen Pseudomonas aeruginosa inhibits translation of target transcripts by forming a regulatory complex with the catabolite repression protein Crc. This repressive complex acts as part of an intricate mechanism of preferred nutrient utilisation. We describe high-resolution cryo-EM structures of the assembly of Hfq and Crc bound to the translation initiation site of a target mRNA. The core of the assembly is formed through interactions of two cognate RNAs, two Hfq hexamers and a Crc pair. Additional Crc protomers are recruited to the core to generate higher-order assemblies with demonstrated regulatory activity in vivo. This study reveals how Hfq cooperates with a partner protein to regulate translation, and provides a structural basis for an RNA code that guides global regulators to interact cooperatively and regulate different RNA targets., eLife digest Living things can adapt rapidly to changes in their surroundings by switching whole groups of genes on and off. These responses must be controlled carefully, and they are often coordinated by regulatory proteins working together. Within a cell, the coded information in genes is copied to create molecules called mRNAs, which are then translated to produce proteins. Stopping the cell from reading the information in mRNAs is one way of shutting down specific genes. Pseudomonas aeruginosa is a species of bacteria that can infect humans and can cause cases of sepsis and pneumonia. It changes the activity of its genes in response to its environment. For example, certain genes are only active inside a human host. This allows the bacteria to make the best use of available nutrients and energy. In these cells, two proteins named Hfq and Crc cooperate to silence groups of genes. They do this by stopping the cell from reading specific mRNA molecules, but how they do this is not fully understood. By using a technique called cryo-electron microscopy, Pei, Dendooven et al. studied Hfq and Crc attached to mRNAs. The results show that two groups of six Hfq molecules and two Crcs attach themselves to two mRNA sections to create a structure that stops an mRNA from being translated into a protein. Since the structure only forms with certain mRNAs, the effect is specific to certain genes. The structure needs both Hfq and Crc to work together, which means it only forms in specific situations – when the affected genes are not needed. P. aeruginosa is highly antibiotic resistant and new drugs are urgently needed to control infections. Understanding how this disease-causing bacterium controls its genes could lead to new treatments. The mechanisms of gene regulation are also common to many other forms of life so may also aid the wider understanding of how cells adapt to rapidly changing environments.

The realization that non protein-coding RNA (ncRNA) is implicated in an increasing number of cellular processes, many related to human disease, makes it imperative to understand and predict RNA folding. RNA secondary structure prediction is more tractable than tertiary structure or protein structure. Yet insights into RNA structure-function relationships are complicated by coupling between RNA folding and ligand binding. Here, we introduce a simple statistical mechanical formalism to calculate perturbations to equilibrium secondary structure conformational distributions for RNA, in the presence of bound cognate ligands. For the first time, this formalism incorporates a key factor in coupling ligand binding to RNA conformation: the differential affinity of the ligand for a range of RNA-folding intermediates. We apply the approach to the SAM-I riboswitch, for which binding data is available for analogs of intermediate secondary structure conformers. Calculations of equilibrium secondary structure distributions during the transcriptional “decision window” predict subtle shifts due to the ligand, rather than an on/off switch. The results suggest how ligand perturbation can release a kinetic block to the formation of a terminator hairpin in the full-length riboswitch. Such predictions identify aspects of folding that are most affected by ligand binding, and can readily be compared with experiment.

The structural flexibility of RNA underlies fundamental biological processes, but there are no methods for exploring the multiple conformations adopted by RNAs in vivo. We developed cross-linking of matched RNAs and deep sequencing (COMRADES) for in-depth RNA conformation capture, and a pipeline for the retrieval of RNA structural ensembles. Using COMRADES, we determined the architecture of the Zika virus RNA genome inside cells, and identified multiple site-specific interactions with human noncoding RNAs. In vivo probing of RNA structures with COMRADES yields insight into RNA folding of the ZIKA virus genome and its interaction with host RNAs.

Fused in sarcoma (FUS) is an RNA binding protein involved in regulating many aspects of RNA process- ing and linked to several neurodegenerative dis- eases. Transcriptomics studies indicate that FUS binds a large variety of RNA motifs, suggesting that FUS RNA binding might be quite complex.Here, we present solution structures of FUS zinc finger (ZnF) and RNA recognition motif (RRM) domains bound to RNA. These structures show a bipartite binding mode of FUS comprising of sequence-spe- cific recognition of a NGGU motif via the ZnF and an unusual shape recognition of a stem-loop RNA via the RRM. In addition, sequence-independent interac- tions via the RGG repeats significantly increase bind- ing affinity and promote destabilization of structured RNA conformation, enabling additional binding. We further show that disruption of the RRM and ZnF domains abolishes FUS function in splicing. Alto- gether, our results rationalize why deciphering the RNA binding mode of FUS has been so challenging.

RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA's recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.

Structural RNA domains are widely involved in the regulation of biological functions, such as gene expression, gene modification, and gene repair. Activity of these dynamic regions depends sensitively on the global fold of the RNA, in particular, on the binding affinity of individual conformations to effector molecules in solution. Consequently, both the 1) structure and 2) conformational dynamics of noncoding RNAs prove to be essential in understanding the coupling that results in biological function. Toward this end, we recently reported observation of three conformational states in the metal-induced folding pathway of the tRNA-like structure domain of Brome Mosaic Virus, via single-molecule fluorescence resonance energy transfer studies. We report herein selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE)-directed structure predictions as a function of metal ion concentrations ([Mn+]) to confirm the three-state folding model, as well as test 2° structure models from the literature. Specifically, SHAPE reactivity data mapped onto literature models agrees well with the secondary structures observed at 0–10 mM [Mg2+], with only minor discrepancies in the E hairpin domain at low [Mg2+]. SHAPE probing and SHAPE-directed structure predictions further confirm the stepwise unfolding pathway previously observed in our single-molecule studies. Of special relevance, this means that reduction in the metal-ion concentration unfolds the 3′ pseudoknot interaction before unfolding the long-range stem interaction. This work highlights the synergistic power of combining 1) single-molecule Forster resonance energy transfer and 2) SHAPE-directed structure-probing studies for detailed analysis of multiple RNA conformational states. In particular, single-molecule guided deconvolution of the SHAPE reactivities permits 2° structure predictions of isolated RNA conformations, thereby substantially improving on traditional limitations associated with current structure prediction algorithms.

The increasing awareness of RNA’s central role in biology calls for a new understanding of how RNAs, like proteins, recognize biological partners. Because RNA is inherently flexible, it assumes a variety of conformations. This conformational flexibility can be a critical aspect of how RNA attracts and binds molecular partners. Structurally, RNA consists of rigid basepaired duplexes, separated by flexible non-basepaired regions. Here, using an RNA system consisting of two short helices, connected by a single-stranded (non-basepaired) junction, we explore the role of helix length and junction sequence in determining the range of conformations available to a model RNA. Single-molecule Förster resonance energy transfer reports on the RNA conformation as a function of either mono- or divalent ion concentration. Electrostatic repulsion between helices dominates at low salt concentration, whereas junction sequence effects determine the conformations at high salt concentration. Near physiological salt concentrations, RNA conformation is sensitive to both helix length and junction sequence, suggesting a means for sensitively tuning RNA conformations.

An increasing number of known RNA 3D structures contributes to the recognition of various RNA families and identification of their features. These tasks are based on an analysis of RNA conformations conducted at different levels of detail. On the other hand, the knowledge of native nucleotide conformations is crucial for structure prediction and understanding of RNA folding. However, this knowledge is stored in structural databases in a rather distributed form. Therefore, only automated methods for sampling the space of RNA structures can reveal plausible conformational representatives useful for further analysis. Here, we present a machine learning-based approach to inspect the dataset of RNA three-dimensional structures and to create a library of nucleotide conformers. A median neural gas algorithm is applied to cluster nucleotide structures upon their trigonometric description. The clustering procedure is two-stage: (i) backbone- and (ii) ribose-driven. We show the resulting library that contains RNA nucleotide representatives over the entire data, and we evaluate its quality by computing normal distribution measures and average RMSD between data points as well as the prototype within each cluster.

Interactions between protein and RNA play a key role in many biological processes in the gene expression pathway. Those interactions are mediated through a variety of RNA-binding protein domains, among them the highly abundant RNA recognition motif (RRM). Here we studied protein-RNA complexes from different RNA binding domain families solved by NMR and x-ray crystallography. Characterizing the structural properties of the RNA at the binding interfaces revealed an unexpected number of nucleotides with unusual RNA conformations, specifically found in RNA-RRM complexes. Moreover, we observed that the RNA nucleotides that are directly involved in interactions with the RRM domains, via hydrogen bonds and hydrophobic contacts, are significantly enriched with unique RNA conformations. Further examination of the sequences binding the RRM domain showed a preference for G nucleotides in syn conformation to precede or to follow U nucleotides in the anti-conformation, and U nucleotides in C2' endo conformation to precede U and G nucleotides possessing the more common C3' endo conformation. These findings imply a possible mode of RNA recognition by the RRM domains which enables the recognition of a wide variety of different RNA sequences and shapes. Overall, this study suggests an additional way by which the RRM domain recognizes its RNA target, involving a conformational readout.

We have developed an advanced microintaglio printing method for fabricating fine and high-density micropatterns and applied it to the microarraying of RNA molecules. The microintaglio printing of RNA reported here is based on the hybridization of RNA with immobilized complementary DNA probes. The hybridization was controlled by switching the RNA conformation via the temperature, and an RNA microarray with a diameter of 1.5 µm and a density of 40,000 spots/mm(2) with high contrast was successfully fabricated. Specifically, no size effects were observed in the uniformity of patterned signals over a range of microarray feature sizes spanning one order of magnitude. Additionally, we have developed a microintaglio printing method for transcribed RNA microarrays on demand using DNA-immobilized magnetic beads. The beads were arrayed on wells fabricated on a printing mold and the wells were filled with in vitro transcription reagent and sealed with a DNA-immobilized glass substrate. Subsequently, RNA was in situ synthesized using the bead-immobilized DNA as a template and printed onto the substrate via hybridization. Since the microintaglio printing of RNA using DNA-immobilized beads enables the fabrication of a microarray of spots composed of multiple RNA sequences, it will be possible to screen or analyze RNA functions using an RNA microarray fabricated by temperature-controlled microintaglio printing (TC-µIP).

Riboswitches are short mRNA sequences that can change their structural conformation to regulate the expression of adjacent genes. They regulate various biological processes via utilizing their secondary structure and consist of two distinct regions: (i) an evolutionarily conserved ligand-binding aptamer region and (ii) a variable expression platform regulating the gene expression. Many ligands bind to riboswitches, and its accurate and selective recognition requires a specific architecture that completely matches a given molecule. In general, the ligand-binding site can be found within the adjacent junction or regions; however, certain ligands may interact with the distant riboswitch regions. Gene expressions can be regulated by switching between two alternative RNA conformations; one of these conformations is favored in the presence of bound metabolite, while the other is favored in its absence. Riboswitches can regulate genes via metabolic pathways that are involved in the biosynthesis of vitamins, amino acids, and purines. In the present chapter, we aim to explain the structure, functions, and biological significance of these molecules.

A comparative study on the interaction of the natural plant alkaloid chelerythrine with triple helical poly(U)·poly(A)*poly(U), double helical poly(A)·poly(U) and single stranded poly(U) (the dot and star representing the Watson-Crick and Hoogsteen base pairing) has been performed using various biophysical and thermodynamic techniques. Chelerythrine binds to the duplex and triplexes in a cooperative manner with affinity of the order of 106 M−1. A weaker binding (∼105 M−1) in a non-cooperative mode occurred with poly(U). Chelerythrine is more selective towards RNA triplex than its parent duplex. The triplex was stabilized specifically without affecting the stability of the duplex. Fluorescence quenching, fluorescence polarization and energy transfer from the nucleotides to the alkaloid, and viscosity results gave convincing evidence for a true intercalative binding of chelerythrine to the triplex and the duplex structures, and partial base stacking with poly(U). The conformations of both double and triple helices were perturbed on binding but no effect occurred to the single strand structure. The binding of the alkaloid to all three RNA helices was found to be exothermic; to the triplex it was entropy driven with favorable enthalpy change, to the duplex enthalpy driven and to the single strand it was enthalpy driven. These results provide new knowledge on the mode, mechanism and specificity, and energetics of binding of the natural alkaloid and putative anticancer agent chelerythrine to different RNA conformations.

Summary Ribonucleoprotein enzymes require dynamic conformations of their RNA constituents for regulated catalysis. Human telomerase employs a non-coding RNA (hTR) with a bipartite arrangement of domains—a template-containing core and a distal three-way junction (CR4/5) that stimulates catalysis through unknown means. Here, we show that telomerase activity unexpectedly depends upon the holoenzyme protein TCAB1, which in turn controls conformation of CR4/5. Cells lacking TCAB1 exhibit a marked reduction in telomerase catalysis without affecting enzyme assembly. Instead, TCAB1 inactivation causes unfolding of CR4/5 helices that are required for catalysis and for association with the telomerase reverse-transcriptase (TERT). CR4/5 mutations derived from patients with telomere biology disorders provoke defects in catalysis and TERT binding similar to TCAB1 inactivation. These findings reveal a conformational "activity switch" in human telomerase RNA controlling catalysis and TERT engagement. The identification of two discrete catalytic states for telomerase suggests an intramolecular means for controlling telomerase in cancers and progenitor cells.

Potato spindle tuber viroid (PSTVd) is a circular, single-stranded, noncoding RNA plant pathogen that is a useful model for studying the processing of noncoding RNA in eukaryotes. Infective PSTVd circles are replicated via an asymmetric rolling circle mechanism to form linear multimeric RNAs. An endonuclease cleaves these into monomers, and a ligase seals these into mature circles. All eukaryotes may have such enzymes for processing noncoding RNA. As a test, we investigated the processing of three PSTVd RNA constructs in the yeast Saccharomyces cerevisiae . Of these, only one form, a construct that adopts a previously described tetraloop-containing conformation (TL), produces circles. TL has 16 nucleotides of the 3′ end duplicated at the 5′ end and a 3′ end produced by self-cleavage of a delta ribozyme. The other two constructs, an exact monomer flanked by ribozymes and a trihelix-forming RNA with requisite 5′ and 3′ duplications, do not produce circles. The TL circles contain nonnative nucleotides resulting from the 3′ end created by the ribozyme and the 5′ end created from an endolytic cleavage by yeast at a site distinct from where potato enzymes cut these RNAs. RNAs from all three transcripts are cleaved in places not on path for circle formation, likely representing RNA decay. We propose that these constructs fold into distinct RNA structures that interact differently with host cell RNA metabolism enzymes, resulting in various susceptibilities to degradation versus processing. We conclude that PSTVd RNA is opportunistic and may use different processing pathways in different hosts. IMPORTANCE In higher eukaryotes, the majority of transcribed RNAs do not encode proteins. These noncoding RNAs are responsible for mRNA regulation, control of the expression of regulatory microRNAs, sensing of changes in the environment by use of riboswitches (RNAs that change shape in response to environmental signals), catalysis, and more roles that are still being uncovered. Some of these functions may be remnants from the RNA world and, as such, would be part of the evolutionary past of all forms of modern life. Viroids are noncoding RNAs that can cause disease in plants. Since they encode no proteins, they depend on their own RNA and on host proteins for replication and pathogenicity. It is likely that viroids hijack critical host RNA pathways for processing the host's own noncoding RNA. These pathways are still unknown. Elucidating these pathways should reveal new biological functions of noncoding RNA.

Comparison of kinetic and thermodynamic properties of IRP1 (iron regulatory protein1) binding to FRT (ferritin) and ACO2 (aconitase2) IRE-RNAs, with or without Mn2+, revealed differences specific to each IRE-RNA. Conserved among animal mRNAs, IRE-RNA structures are noncoding and bind Fe2+ to regulate biosynthesis rates of the encoded, iron homeostatic proteins. IRP1 protein binds IRE-RNA, inhibiting mRNA activity; Fe2+ decreases IRE-mRNA/IRP1 binding, increasing encoded protein synthesis. Here, we observed heat, 5 °C to 30 °C, increased IRP1 binding to IRE-RNA 4-fold (FRT IRE-RNA) or 3-fold (ACO2 IRE-RNA), which was enthalpy driven and entropy favorable. Mn2+ (50 µM, 25 °C) increased IRE-RNA/IRP1 binding (Kd) 12-fold (FRT IRE-RNA) or 6-fold (ACO2 IRE-RNA); enthalpic contributions decreased ~61% (FRT) or ~32% (ACO2), and entropic contributions increased ~39% (FRT) or ~68% (ACO2). IRE-RNA/IRP1 binding changed activation energies: FRT IRE-RNA 47.0 ± 2.5 kJ/mol, ACO2 IRE-RNA 35.0 ± 2.0 kJ/mol. Mn2+ (50 µM) decreased the activation energy of RNA-IRP1 binding for both IRE-RNAs. The observations suggest decreased RNA hydrogen bonding and changed RNA conformation upon IRP1 binding and illustrate how small, conserved, sequence differences among IRE-mRNAs selectively influence thermodynamic and kinetic selectivity of the protein/RNA interactions.

Several protein-targeted RNA aptamers have been identified for a variety of applications and although the affinities of numerous protein-aptamer complexes have been determined, the structural details of these complexes have not been widely explored. We examined the structural accommodation of an RNA aptamer that binds bacterial r-protein S8. The core of the primary binding site for S8 on helix 21 of 16S rRNA contains a pair of conserved base triples that mold the sugar-phosphate backbone to S8. The aptamer, which does not contain the conserved sequence motif, is specific for the rRNA binding site of S8. The protein-free RNA aptamer adopts a helical structure with multiple non-canonical base pairs. Surprisingly, binding of S8 leads to a dramatic change in the RNA conformation that restores the signature S8 recognition fold through a novel combination of nucleobase interactions. Nucleotides within the non-canonical core rearrange to create a G-(G-C) triple and a U-(A-U)-U quartet. Although native-like S8-RNA interactions are present in the aptamer-S8 complex, the topology of the aptamer RNA differs from that of the helix 21-S8 complex. This is the first example of an RNA aptamer that adopts substantially different secondary structures in the free and protein-bound states and highlights the remarkable plasticity of RNA secondary structure.

Internal ribosome entry site (IRES) elements are high-order RNA structures that promote internal initiation of translation to allow protein synthesis under situations that compromise the general cap-dependent translation mechanism. Picornavirus IRES elements are highly efficient elements with a modular RNA structure organization. Here we investigated the effect of Mg(2+) concentration on the local flexibility and solvent accessibility of the foot-and-mouth disease virus (FMDV) IRES element measured on the basis of selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivity and hydroxyl radical cleavage. We have found that Mg(2+) concentration affects the organization of discrete IRES regions, mainly the apical region of domain 3, the 10 nt loop of domain 4, and the pyrimidine tract of domain 5. In support of the effect of RNA structure on IRES activity, substitution or deletion mutants of the 10 nt loop of domain 4 impair internal initiation. In addition, divalent cations affect the binding of eIF4G, a eukaryotic initiation factor that is essential for IRES-dependent translation that interacts with domain 4. Binding of eIF4G is favored by the local RNA flexibility adopted at low Mg(2+) concentration, while eIF4B interacts with the IRES independently of the compactness of the RNA structure. Our study shows that the IRES element adopts a near-native structure in the absence of proteins, shedding light on the influence of Mg(2+) ions on the local flexibility and binding of eIF4G in a model IRES element.

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

A, U, C, and G are boring. Truly, four nucleosides sometimes repeated in seemingly endless sequence has a numbing effect even in a computational and structural context. However, there is nothing boring about the chemistries and structures of post-transcriptional modifications of RNA as indicated by the rising interest and rapidly increasing numbers of publications. A plethora of papers demonstrate that modified nucleosides of RNA are functionally indispensable for RNA processing and translation. Over the last 20 years, efficiency, accuracy, enhanced activity and responsiveness are terms that accompany studies of modifications. Modifications are found in all types of RNAs, in every structural motif and in unstructured regions. The large numbers of modifications in tRNAs are quite familiar. The presence, chemistries and RNA sequence locations of modified nucleosides are only now beginning to be appreciated as both temporal—in response to environmental (internal and external) signals—and spatial, in their effect on RNA conformation, dynamics, and functions. Modifications can be recognition determinants for proteins as well as incredible determiners of RNA structure and function. Unfortunately, we tend to neglect the modified nucleosides. Using amplification and transcription in vitro that disregard the native modified nucleosides, we study the interesting low abundance, non-coding RNAs for which one molecule per cell is sometimes enough to affect gene expression. On the other hand, we lack technologies sensitive enough to explore modifications and their functions in these same low abundance RNAs. Our lack of understanding of RNA modification regrettably limits innovation in the use of modified nucleosides for synthetic biology and medicine. However, the past two decades have led to some of the most interesting of findings.

Folding messenger RNA into specific structures is a common regulatory mechanism involved in translation. In Escherichia coli, the operator of the rpsO gene transcript folds into a pseudoknot or double-hairpin conformation. S15, the gene product, binds only to the pseudoknot, thereby repressing its own synthesis when it is present in excess in the cell. The two RNA conformations have been proposed to exist in equilibrium. However, it remained unclear how structural changes can be achieved between these two topologically distinct conformations. We used optical tweezers to study the structural dynamics and rearrangements of the rpsO operator RNA at the single-molecule level. We discovered that the two RNA structures can be interchanged spontaneously and the pseudoknot can exist in conformations that exhibit various levels of stability. Conversion from the double hairpin to a pseudoknot through potential hairpin–hairpin interactions favoured the high-stability conformation. By contrast, mutations that blocked the formation of a hairpin typically resulted in alternative low-stability pseudoknots. These results demonstrate that specific tertiary interactions of RNA can be established and modulated based on the interactions and rearrangements between secondary structural components. Our findings provide new insight into the RNA folding pathway that leads to a regulatory conformation for target protein binding.

It has long been appreciated that Mg(2+) is essential for the stabilization of RNA tertiary structure. However, the problem of quantitative prediction for the ion effect in tertiary structure folding remains. By using the virtual bond RNA folding model (Vfold) to generate RNA conformations and the newly improved tightly bound ion model (TBI) to treat ion-RNA interactions, we investigate Mg(2+)-facilitated tetraloop-receptor docking. For the specific construct of the tetraloop-receptor system, the theoretical analysis shows that the Mg(2+)-induced stabilizing force for the docked state is predominantly entropic and the major contribution comes from the entropy of the diffusive ions. Furthermore, our results show that Mg(2+) ions promote tetraloop-receptor docking mainly through the entropy of the diffusive ions. The theoretical prediction agrees with experimental analysis. The method developed in this paper, which combines the theory for the (Mg(2+)) ion effects in RNA folding and RNA conformational sampling, may provide a useful framework for studying the ion effect in the folding of more complex RNA structures.

Neomycin dimers synthesized using "click chemistry" with varying functionality and length in the linker region have been shown to be effective in targeting the HIV-1 TAR RNA region of the HIV virus. TAR (Transactivation Response) RNA region, a 59 base pair stem loop structure located at the 5'-end of all nascent viral transcripts interacts with its target, a key regulatory protein, Tat, and necessitates the replication of HIV-1 virus. Ethidium bromide displacement and FRET competition assays have revealed nanomolar binding affinity between neomycin dimers and wildtype TAR RNA while in case of neomycin, only a weak binding was detected. Here, NMR and FID-based comparisons reveal an extended binding interface for neomycin dimers involving the upper stem of the TAR RNA thereby offering an explanation for increased affinities. To further explore the potential of these modified aminosugars we have extended binding studies to include four TAR RNA mutants that display conformational differences with minimal sequence variation. The differences in binding between neomycin and neomycin dimers is characterized with TAR RNA mutants that include mutations to the bulge region, hairpin region, and both the bulge and hairpin regions. Our results demonstrate the effect of these mutations on neomycin binding and our results show that linker functionalities between dimeric units of neomycin can distinguish between the conformational differences of mutant TAR RNA structures.

RNA sequences having up to three consecutive internal amide linkages were synthesized and studied using UV and NMR spectroscopy. The amide modifications did not interfere with normal base-pairing and A-type RNA conformation. Three consecutive amides were well tolerated in the passenger strand of siRNA and caused little change in RNAi activity.

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

RNA molecules have highly versatile structures that can fold into myriad conformations, providing many potential pockets for binding small molecules. The increasing number of available RNA structures, in complex with proteins, small ligands and in free form, enables the design of new therapeutically useful RNA-binding ligands. Here we studied RNA ligand complexes from 10 RNA groups extracted from the protein data bank (PDB), including adaptive and non-adaptive complexes. We analyzed the chemical, physical, structural and conformational properties of binding pockets around the ligand. Comparing the properties of ligand-binding pockets to the properties of computed pockets extracted from all available RNA structures and RNA-protein interfaces, revealed that ligand-binding pockets, mainly the adaptive pockets, are characterized by unique properties, specifically enriched in rare conformations of the nucleobase and the sugar pucker. Further, we demonstrate that nucleotides possessing the rare conformations are preferentially involved in direct interactions with the ligand. Overall, based on our comprehensive analysis of RNA-ligand complexes, we suggest that the unique conformations adopted by RNA nucleotides play an important role in RNA recognition by small ligands. We term the recognition of a binding site by a ligand via the unique RNA conformations “RNA conformational readout.” We propose that “conformational readout” is a general way by which RNA binding pockets are recognized and selected from an ensemble of different RNA states.

Picornavirus RNAs initiate translation using a 5′ end-independent mechanism based on internal ribosome entry site (IRES) elements. Despite performing similar functions, IRES elements present in genetically distant RNAs differ in primary sequence, RNA secondary structure and trans-acting factors requirement. The lack of conserved features amongst IRESs represents obstacles for the understanding of the internal initiation process. RNA structure is tightly linked to picornavirus IRES activity, consistent with the conservation of RNA motifs. This study extends the functional relevance of evolutionary conserved motifs of foot-and-mouth disease virus (FMDV) IRES. SHAPE structural analysis of mutant IRESs revealed local changes in RNA flexibility indicating the existence of an interactive structure constrained by lateral bulges that maintain the RNA conformation necessary for IRES-mediated translation. © 2013 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved., BFU2011-25437; Fundación Ramón Areces

Genuine La and La-related proteins group 7 (LARP7) bind to the non-coding RNAs transcribed by RNA polymerase III (RNAPIII), which end in UUU-3'OH. The La motif and RRM1 of these proteins (the La module) cooperate to bind the UUU-3'OH, protecting the RNA from degradation, while other domains may be important for RNA folding or other functions. Among the RNAPIII transcripts is ciliate telomerase RNA (TER). p65, a member of the LARP7 family, is an integral Tetrahymena thermophila telomerase holoenzyme protein required for TER biogenesis and telomerase RNP assembly. p65, together with TER and telomerase reverse transcriptase (TERT), form the Tetrahymena telomerase RNP catalytic core. p65 has an N-terminal domain followed by a La module and a C-terminal domain, which binds to the TER stem 4. We recently showed that the p65 C-terminal domain harbors a cryptic, atypical RRM, which uses a unique mode of single- and double-strand RNA binding and is required for telomerase RNP catalytic core assembly. This domain, which we named xRRM, appears to be present in and unique to genuine La and LARP7 proteins. Here we review the structure of the xRRM, discuss how this domain could recognize diverse substrates of La and LARP7 proteins and discuss the functional implications of the xRRM as an RNP chaperone.

One of the multitasking proteins, transactive response DNA-binding protein 43 (tdp43) plays a key role in RNA regulation and the two pathogenic mutations such as D169G and K263E, located at the RNA Recognition Motif (RRM) of tdp43, are reported to cause neurological disorders such as Amyotrophic Lateral Sclerosis and FrontoTemporal Lobar Degeneration. As the exploration of the proteinopathy demands both structural and functional characterizations of mutants, a comparative analysis on the wild type and mutant tdp43 (D169G and K263E) and their complexes with RNA has been performed using computational approaches. Molecular dynamics simulations revealed comparatively stable mutant structures compared to wild type tdp43. Both mutants show lesser binding affinity toward RNA molecule when compared to the wild type tdp43. Some of the observed features, including the increased solvent-accessible surface area, conformational flexibility as well as unfolding of tdp43, and the altered RNA conformation in tp43-RNA complex, reveal the susceptibility of these mutants to induce conformational changes in tdp43 for a possible aggregation in the cytoplasm. Particularly, the enhanced aggregation propensity of both mutants also evidences the higher probability of cytoplasmic aggregation of tdp43 mutants. Hence, the present analysis highlighting the structural and functional aspects of wild and mutant tdp43 will form the basis to gain insight into the proteinopathy of tdp43 and the related structure-based drug discovery. Thus, tdp43 can be used as target to develop novel therapeutic approaches or drug designing.

Antibodies that bind RNA represent an unrealized source of reagents for synthetic biology and for characterizing cellular transcriptomes. However, facile access to RNA-binding antibodies requires the engineering of effective Fab libraries guided by the knowledge of the principles that govern RNA recognition. Here, we describe a Fab identified from a minimalist synthetic library during phage display against a branched RNA target. The Fab (BRG) binds with 20 nM dissociation constant to a single-stranded RNA (ssRNA) sequence adjacent to the branch site and can block the action of debranchase enzyme. We report the crystal structure in complex with RNA target at 2.38 A. The Fab traps the RNA in a hairpin conformation that contains a 2-bp duplex capped by a tetraloop. The paratope surface consists of residues located in four complementarity-determining regions including a major contribution from H3, which adopts a helical structure that projects into a deep, wide groove formed by the RNA. The amino acid composition of the paratope reflects the library diversity, consisting mostly of tyrosine and serine residues and a small but significant contribution from a single arginine residue. This structure, involving the recognition of ssRNA via a stem–loop conformation, together with our two previous structures involving the recognition of an RNA hairpin loop and an RNA tertiary structure, reveals the capacity of minimalist libraries biased with tyrosine, serine, glycine, and arginine to form binding surfaces for specific RNA conformations and distinct levels of RNA structural hierarchy.

RNA splicing is a critical step in the human immunodeficiency virus type 1 (HIV-1) replication cycle because it controls the expression of the complex viral proteome. The major 5′ splice site (5′ss) that is positioned in the untranslated leader of the HIV-1 RNA transcript is of particular interest because it is used for the production of the more than 40 differentially spliced subgenomic mRNAs. HIV-1 splicing needs to be balanced tightly to ensure the proper levels of all viral proteins, including the Gag-Pol proteins that are translated from the unspliced RNA. We previously presented evidence that the major 5′ss is regulated by a repressive local RNA structure, the splice donor (SD) hairpin, that masks the 11 nucleotides (nts) of the 5′ss signal for recognition by U1 small nuclear RNA (snRNA) of the spliceosome machinery. A strikingly different multiple-hairpin RNA conformation was recently proposed for this part of the HIV-1 leader RNA. We therefore inspected the sequence of natural HIV-1 isolates in search for support, in the form of base pair (bp) co-variations, for the different RNA conformations.

The structure of RNA, which is considered to be a second layer of information alongside the genetic code, provides fundamental insights into the cellular function of both coding and non-coding RNAs. Several high-throughput technologies have been developed to profile transcriptome-wide RNA structures, i.e., the structurome. However, it is challenging to interpret the profiling data because the observed data represent an average over different RNA conformations and isoforms with different abundance. To address this challenge, we developed an RNA structurome quantification method (RSQ) to statistically model the distribution of reads over both isoforms and RNA conformations, and thus provide accurate quantification of the isoform-specific structurome. The quantified RNA structurome enables the comparison of isoform-specific conformations between different conditions, the exploration of RNA conformation variation affected by single nucleotide polymorphism (SNP),and the measurement of RNA accessibility for binding of either small RNAs in RNAi-based assays or RNA binding protein in transcriptional regulation. The model used in our method sheds new light on the potential impact of the RNA structurome on gene regulation.

Single-nucleotide polymorphisms (SNPs) are often linked to critical phenotypes such as diseases or responses to vaccines, medications and environmental factors. However, the specific molecular mechanisms by which a causal SNP acts is usually not obvious. Changes in RNA secondary structure emerge as a possible explanation necessitating the development of methods to measure the impact of single-nucleotide variation on RNA structure. Despite the recognition of the importance of considering the changes in Boltzmann ensemble of RNA conformers in this context, a formal method to perform directly such comparison was lacking. Here, we solved this problem and designed an efficient method to compute the relative entropy between the Boltzmann ensembles of the native and a mutant structure. On the basis of this theoretical progress, we developed a software tool, remuRNA, and investigated examples of its application. Comparing the impact of common SNPs naturally occurring in populations with the impact of random point mutations, we found that structural changes introduced by common SNPs are smaller than those introduced by random point mutations. This suggests a natural selection against mutations that significantly change RNA structure and demonstrates, surprisingly, that randomly inserted point mutations provide inadequate estimation of random mutations effects. Subsequently, we applied remuRNA to determine which of the disease-associated non-coding SNPs are potentially related to RNA structural changes.

Genome packaging is an essential step in virus replication and a potential drug target. Single-stranded RNA viruses have been thought to encapsidate their genomes by gradual co-assembly with capsid subunits. In contrast, using a single molecule fluorescence assay to monitor RNA conformation and virus assembly in real time, with two viruses from differing structural families, we have discovered that packaging is a two-stage process. Initially, the genomic RNAs undergo rapid and dramatic (approximately 20–30%) collapse of their solution conformations upon addition of cognate coat proteins. The collapse occurs with a substoichiometric ratio of coat protein subunits and is followed by a gradual increase in particle size, consistent with the recruitment of additional subunits to complete a growing capsid. Equivalently sized nonviral RNAs, including high copy potential in vivo competitor mRNAs, do not collapse. They do support particle assembly, however, but yield many aberrant structures in contrast to viral RNAs that make only capsids of the correct size. The collapse is specific to viral RNA fragments, implying that it depends on a series of specific RNA–protein interactions. For bacteriophage MS2, we have shown that collapse is driven by subsequent protein–protein interactions, consistent with the RNA–protein contacts occurring in defined spatial locations. Conformational collapse appears to be a distinct feature of viral RNA that has evolved to facilitate assembly. Aspects of this process mimic those seen in ribosome assembly.

Experiments demonstrate that Mg(2+) is crucial for structure and function of RNA systems, yet the detailed molecular mechanism of Mg(2+) action on RNA is not well understood. We investigate the interplay between RNA and Mg(2+) at atomic resolution through ten 2-μs explicit solvent molecular dynamics simulations of the SAM-I riboswitch with varying ion concentrations. The structure, including three stemloops, is very stable on this time scale. Simulations reveal that outer-sphere coordinated Mg(2+) ions fluctuate on the same time scale as the RNA, and that their dynamics couple. Locally, Mg(2+) association affects RNA conformation through tertiary bridging interactions; globally, increasing Mg(2+) concentration slows RNA fluctuations. Outer-sphere Mg(2+) ions responsible for these effects account for 80% of Mg(2+) in our simulations. These ions are transiently bound to the RNA, maintaining interactions, but shuttled from site to site. Outer-sphere Mg(2+) are separated from the RNA by a single hydration shell, occupying a thin layer 3-5 Å from the RNA. Distribution functions reveal that outer-sphere Mg(2+) are positioned by electronegative atoms, hydration layers, and a preference for the major groove. Diffusion analysis suggests transient outer-sphere Mg(2+) dynamics are glassy. Since outer-sphere Mg(2+) ions account for most of the Mg(2+) in our simulations, these ions may change the paradigm of Mg(2+)-RNA interactions. Rather than a few inner-sphere ions anchoring the RNA structure surrounded by a continuum of diffuse ions, we observe a layer of outer-sphere coordinated Mg(2+) that is transiently bound but strongly coupled to the RNA.

The importance of RNA tertiary structure is evident from the growing number of published high resolution NMR and X-ray crystallographic structures of RNA molecules. These structures provide insights into function and create a knowledge base that is leveraged by programs such as Assemble, ModeRNA, RNABuilder, NAST, FARNA, Mc-Sym, RNA2D3D, and iFoldRNA for tertiary structure prediction and design. While these methods sample native-like RNA structures during simulations, all struggle to capture the native RNA conformation after scoring. We propose RSIM, an improved RNA fragment assembly method that preserves RNA global secondary structure while sampling conformations. This approach enhances the quality of predicted RNA tertiary structure, provides insights into the native state dynamics, and generates a powerful visualization of the RNA conformational space. RSIM is available for download from http://www.github.com/jpbida/rsim.

Changes in RNA conformation can alter gene expression. The guanine quadruplex sequence (GQS) is an RNA motif that folds in the presence of K(+) ions. Changes in the conformation of this motif could be especially important in regulating gene expression in plants because intracellular K(+) concentrations often increase during drought stress. Little is known about the folding thermodynamics of RNA GQS. We show here that RNA GQS with tracts containing three G's [e.g., (GGGxx)(4)] have a modest dependence on the K(+) concentration, folding with no or even negative cooperativity (Hill coefficients ≤1), and are associated with populated folding intermediates. In contrast, GQS with tracts containing just two G's [e.g., (GGxx)(4)] have a steep dependence on the K(+) concentration and fold with positive cooperativity (Hill coefficients of 1.7-2.7) without significantly populating intermediate states. We postulate that in plants, the more stable G3 sequences are largely folded even under unstressed conditions, while the less stable G2 sequences fold only at the higher K(+) concentrations associated with cellular stress, wherein they respond sharply to changing K(+) concentrations. Given the binary nature of their folding, G2 sequences may find application in computation with DNA and in engineering of genetic circuits.

The local conformation of RNA molecules is an important factor in determining their catalytic and binding properties. The analysis of such conformations is particularly difficult due to the large number of degrees of freedom, such as the measured torsion angles per residue and the interatomic distances among interacting residues. In this work, we use a nearest-neighbor search method based on the statistical mechanical Potts model to find clusters in the RNA conformational space. The proposed technique is mostly automatic and may be applied to problems, where there is no prior knowledge on the structure of the data space in contrast to many other clustering techniques. Results are reported for both single residue conformations, where the parameter set of the data space includes four to seven torsional angles, and base pair geometries, where the data space is reduced to two dimensions. Moreover, new results are reported for base stacking geometries. For the first two cases, i.e., single residue conformations and base pair geometries, we get a very good match between the results of the proposed clustering method and the known classifications with only few exceptions. For the case of base stacking geometries, we validate our classification with respect to geometrical constraints and describe the content, and the geometry of the new clusters.

Bamboo mosaic virus (BaMV) satellite RNA (satBaMV) depends on BaMV for its replication and encapsidation. SatBaMV-encoded P20 protein is an RNA-binding protein that facilitates satBaMV systemic movement in co-infected plants. Here, we examined phosphorylation of P20 and its regulatory functions. Recombinant P20 (rP20) was phosphorylated by host cellular kinase(s) in vitro, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and mutational analyses revealed Ser-11 as the phosphorylation site. The phosphor-mimic rP20 protein interactions with satBaMV-translated mutant P20 were affected. In overlay assay, the Asp mutation at S11 (S11D) completely abolished the self-interaction of rP20 and significantly inhibited the interaction with both the WT and S11A rP20. In chemical cross-linking assays, S11D failed to oligomerize. Electrophoretic mobility shift assay and subsequent Hill transformation analysis revealed a low affinity of the phospho-mimicking rP20 for satBaMV RNA. Substantial modulation of satBaMV RNA conformation upon interaction with nonphospho-mimic rP20 in circular dichroism analysis indicated formation of stable satBaMV ribonucleoprotein complexes. The dissimilar satBaMV translation regulation of the nonphospho- and phospho-mimic rP20 suggests that phosphorylation of P20 in the ribonucleoprotein complex converts the translation-incompetent satBaMV RNA to messenger RNA. The phospho-deficient or phospho-mimicking P20 mutant of satBaMV delayed the systemic spread of satBaMV in co-infected Nicotiana benthamiana with BaMV. Thus, satBaMV likely regulates the formation of satBaMV RNP complex during co-infection in planta.

Riboswitch RNAs fold into complex tertiary structures upon binding to their cognate ligand. Ligand recognition is accomplished by key residues in the binding pocket. In addition, it often crucially depends on the stability of peripheral structural elements. The ligand-bound complex of the guanine-sensing riboswitch from Bacillus subtilis, for example, is stabilized by extensive interactions between apical loop regions of the aptamer domain. Previously, we have shown that destabilization of this tertiary loop-loop interaction abrogates ligand binding of the G37A/C61U-mutant aptamer domain (Gsw(loop)) in the absence of Mg(2+). However, if Mg(2+) is available, ligand-binding capability is restored by a population shift of the ground-state RNA ensemble toward RNA conformations with pre-formed loop-loop interactions. Here, we characterize the striking influence of long-range tertiary structure on RNA folding kinetics and on ligand-bound complex structure, both by X-ray crystallography and time-resolved NMR. The X-ray structure of the ligand-bound complex reveals that the global architecture is almost identical to the wild-type aptamer domain. The population of ligand-binding competent conformations in the ground-state ensemble of Gsw(loop) is tunable through variation of the Mg(2+) concentration. We quantitatively describe the influence of distinct Mg(2+) concentrations on ligand-induced folding trajectories both by equilibrium and time-resolved NMR spectroscopy at single-residue resolution.

Riboswitches are functional mRNA that control gene expression. Thiamine pyrophosphate (TPP) binds to thi-box riboswitch RNA and allosterically inhibits genes that code for proteins involved in the biosynthesis and transport of thiamine. Thiamine binding to the pyrimidine sensor helix and pyrophosphate binding to the pyrophosphate sensor helix cause changes in RNA conformation that regulate gene expression. Here we examine the thermodynamic properties of the internal loop of the pyrophosphate binding domain by comparing the wild-type construct (RNA WT) with six modified 2 × 2 bulged RNA and one 2 × 2 bulged DNA. The wild-type construct retains five conserved bases of the pyrophosphate sensor domain, two of which are in the 2 × 2 bulge (C65 and G66). The RNA WT construct was among the most stable (ΔG°37 = −7.7 kcal/mol) in 1 M KCl at pH 7.5. Breaking the A•G mismatch of the bulge decreases the stability of the construct ∼0.5–1 kcal/mol, but does not affect magnesium binding to the RNA WT. Guanine at position 48 is important for RNA–Mg2+ interactions of the TPP-binding riboswitch at pH 7.5. In the presence of 9.5 mM magnesium at pH 5.5, the bulged RNA constructs gained an average of 1.1 kcal/mol relative to 1 M salt. Formation of a single A+•C mismatch base pair contributes about 0.5 kcal/mol at pH 5.5, whereas two tandem A+•C mismatch base pairs together contribute about 2 kcal/mol.