Allosteric mechanism of the V. vulnificus adenine riboswitch resolved by four-dimensional chemical mapping

Decoding the effect of temperatures on conformational stability and order of ligand unbound thermosensing adenine riboswitch using molecular dynamics simulation

The structure-function relationship of the riboswitch is governed mainly by two factors, ligand binding and temperature. Most of the experimental studies shed light on structural dynamics and gene regulation function of Adenine riboswitch from the aspect of ligand instead of temperature. Two unliganded Adenine riboswitch conformations (apoA and apoB) from the thermophile Vibrio vulnificus draw particular attention due to their diverse and polymorphic structures. Ligand-free apoB Adenine riboswitch conformation is not able to interact with the ligand whereas ligand-free apoA Adenine riboswitch conformation adopts ligand-receptive form. The interconversion between apoA and apoB conformation is temperature-dependent and thermodynamically controlled. Therefore Adenine riboswitch is called temperature sensing RNA. The molecular mechanism underlying the thermosensitivity of ligand free Adenine riboswitch is not well known. Hence an attempt is made to examine conformational stability and order of apoA with respect to apoB Adenine riboswitch aptamer computing RMSD, RMSF, RG, principal component analysis, hydrogen bonding interaction and conformational thermodynamics derived from all-atom molecular dynamics trajectories in the temperature range 283K-400K. The temperatures corresponding to the conformational stability and order of apoA adenine riboswitch with respect to apoB adenine riboswitch whole aptamer are shown in descending order 293K∼303K> 313K∼283K>373K>323K. Residue wise and domain wise changes in conformational free energy and entropy of conformational degrees of freedom like pseudo-torsion angle ƞ and θ reflect apoA exhibits pronounced conformational stability compared to apoB at temperatures 293K and 303K, whereas both forms reveal decreased stability at 323K and 400K and may be inactivated, highlighting their role as temperature sensors.

SEISMICgraph: a web-based tool for RNA structure data visualization

In recent years, RNA has been increasingly recognized for its essential roles in biology, functioning not only as a carrier of genetic information but also as a dynamic regulator of gene expression through its interactions with other RNAs, proteins, and itself. Advances in chemical probing techniques have significantly enhanced our ability to identify RNA secondary structures and understand their regulatory roles. These developments, alongside improvements in experimental design and data processing, have greatly increased the resolution and throughput of structural analyses. Here, we introduce SEISMICgraph, a web-based tool designed to support RNA structure research by offering data visualization and analysis capabilities for a variety of chemical probing modalities. SEISMICgraph enables simultaneous comparison of data across different sequences and experimental conditions through a user-friendly interface that requires no programming expertise. We demonstrate its utility by investigating known and putative riboswitches and exploring how RNA modifications influence their structure and binding. SEISMICgraph's ability to rapidly visualize adenine-dependent structural changes and assess the impact of pseudouridylation on these transitions provides novel insights and establishes a roadmap for numerous future applications.

Probing RNA structure and dynamics using nanopore and next generation sequencing

It has become increasingly evident that the structures RNAs adopt are conformationally dynamic; the various structured states that RNAs sample govern their interactions with other nucleic acids, proteins, and ligands to regulate a myriad of biological processes. Although several biophysical approaches have been developed and used to study the dynamic landscape of structured RNAs, technical limitations have limited their application to all classes of RNA due to variable size and flexibility. Recent advances combining chemical probing experiments with next-generation- and direct sequencing have emerged as an alternative approach to exploring the conformational dynamics of RNA. In this review, we provide a methodological overview of the sequencing-based techniques used to study RNA conformational dynamics. We discuss how different techniques have enabled us to better understand the propensity of RNAs from a variety of different classes to sample multiple conformational states. Finally, we present examples of the ways these techniques have reshaped how we think about RNA structure.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Structural insights into translation regulation by the THF-II riboswitch

In bacteria, expression of folate-related genes is controlled by the tetrahydrofolate (THF) riboswitch in response to specific binding of THF and its derivatives. Recently, a second class of THF riboswitches, named THF-II, was identified in Gram-negative bacteria, which exhibit distinct architecture from the previously characterized THF-I riboswitches found in Gram-positive bacteria. Here, we present the crystal structures of the ligand-bound THF-II riboswitch from Mesorhizobium loti. These structures exhibit a long rod-like fold stabilized by continuous base pair and base triplet stacking across two helices of P1 and P2 and their interconnecting ligand-bound binding pocket. The pterin moiety of the ligand docks into the binding pocket by forming hydrogen bonds with two highly conserved pyrimidines in J12 and J21, which resembles the hydrogen-bonding pattern at the ligand-binding site FAPK in the THF-I riboswitch. Using small-angle X-ray scattering and isothermal titration calorimetry, we further characterized the riboswitch in solution and reveal that Mg2+ is essential for pre-organization of the binding pocket for efficient ligand binding. RNase H cleavage assay indicates that ligand binding reduces accessibility of the ribosome binding site in the right arm of P1, thus down-regulating the expression of downstream genes. Together, these results provide mechanistic insights into translation regulation by the THF-II riboswitch.

Potential effects of metal ion induced two-state allostery on the regulatory mechanism of add adenine riboswitch

Riboswitches normally regulate gene expression through structural changes in response to the specific binding of cellular metabolites or metal ions. Taking add adenine riboswitch as an example, we explore the influences of metal ions (especially for K⁺ and Mg²⁺ ions) on the structure and dynamics of riboswitch aptamer (with and without ligand) by using molecular dynamic (MD) simulations. Our results show that a two-state transition marked by the structural deformation at the connection of J12 and P1 (C_J12-P1) is not only related to the binding of cognate ligands, but also strongly coupled with the change of metal ion environments. Moreover, the deformation of the structure at C_J12-P1 can be transmitted to P1 directly connected to the expression platform in multiple ways, which will affect the structure and stability of P1 to varying degrees, and finally change the regulation state of this riboswitch.

Molecular dynamic simulations are employed to assess the influence of metal ions on riboswitch structure and dynamics, suggesting a conformational control of riboswitch aptamers by metal ions before ligand binding.

Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

Incorporation of a FRET Pair into a Riboswitch RNA to Measure Mg2+ Concentration and RNA Conformational Change in Cell

Riboswitches are natural biosensors that can regulate gene expression by sensing small molecules. Knowledge of the structural dynamics of riboswitches is crucial to elucidate their regulatory mechanism and develop RNA biosensors. In this work, we incorporated the fluorophore, Cy3, and its quencher, TQ3, into a full-length adenine riboswitch RNA and its isolated aptamer domain to monitor the dynamics of the RNAs in vitro and in cell. The adenine riboswitch was sensitive to Mg²⁺ concentrations and could be used as a biosensor to measure cellular Mg²⁺ concentrations. Additionally, the TQ3/Cy3-labeled adenine riboswitch yielded a Mg²⁺ concentration that was similar to that measured using a commercial assay kit. Furthermore, the fluorescence response to the adenine of the TQ3/Cy3-labeled riboswitch RNA was applied to determine the proportions of multiple RNA conformational changes in cells. The strategy developed in this work can be used to probe the dynamics of other RNAs in cells and may facilitate the developments of RNA biosensors, drugs and engineering.

A transient conformation facilitates ligand binding to the adenine riboswitch

RNAs adopt various conformations to perform different functions in cells. Incapable of acquiring intermediates, the key initiations of ligand recognition in the adenine riboswitch have not been characterized. In this work, stopped-flow fluorescence was used to track structural switches in the full-length adenine riboswitch in real time. We used PLOR (position-selective labeling of RNA) to incorporate fluorophores into desired positions in the RNA. The switching sequence P1 responded to adenine more rapidly than helix P4 and the binding pocket, followed by stabilization of the binding pocket, P4, and annealing of P1. Moreover, a transient intermediate consisting of an unwound P1 was detected during adenine binding. These events were observed in both the WT riboswitch and a functional mutant. The findings provide insight into the conformational changes of the riboswitch RNA triggered by a ligand.

•

Real-time tracking of the adenine riboswitch at nucleotide resolution

•

A transient conformation with unwound P1 is identified in the adenine riboswitch

•

Helix P1 responds to ligand quicker than the binding pocket or expression platform

Thermal adaptation of structural dynamics and regulatory function of adenine riboswitch

Ligand binding and temperature play important roles in riboswitch RNAs' structures and functions. However, most studies focused on studying structural dynamics or gene-regulation function of riboswitches from the aspect of ligand, instead of temperature. Here we combined NMR, ITC, stopped-flow and in vivo assays to investigate the ligand-triggered switch of adenine riboswitch from 10 to 45°C. Our results demonstrated that at single-nucleotide resolution, structural regions sensed ligand and temperature diversely. Temperature had opposite effects on ligand-binding and gene-regulation of adenine riboswitch. Compared with higher temperature, the RNA bound with its cognate ligand obviously stronger, while its regulatory capacity was weakened at lower temperature. In addition, application of specific-labelled RNAs to the stopped-flow experiments identified the real-time folding of the specific positions upon ligand addition at different temperatures. The kissing loop and internal loop at the riboswitch responded to ligand and temperature differently. The distinct thermo-dynamics of adenine riboswitch exposed here may contribute to the fields of RNA sensors and drug design.

RNAProbe: a web server for normalization and analysis of RNA structure probing data

RNA molecules play key roles in all living cells. Knowledge of the structural characteristics of RNA molecules allows for a better understanding of the mechanisms of their action. RNA chemical probing allows us to study the susceptibility of nucleotides to chemical modification, and the information obtained can be used to guide secondary structure prediction. These experimental results can be analyzed using various computational tools, which, however, requires additional, tedious steps (e.g., further normalization of the reactivities and visualization of the results), for which there are no fully automated methods. Here, we introduce RNAProbe, a web server that facilitates normalization, analysis, and visualization of the low-pass SHAPE, DMS and CMCT probing results with the modification sites detected by capillary electrophoresis. RNAProbe automatically analyzes chemical probing output data and turns tedious manual work into a one-minute assignment. RNAProbe performs normalization based on a well-established protocol, utilizes recognized secondary structure prediction methods, and generates high-quality images with structure representations and reactivity heatmaps. It summarizes the results in the form of a spreadsheet, which can be used for comparative analyses between experiments. Results of predictions with normalized reactivities are also collected in text files, providing interoperability with bioinformatics workflows. RNAProbe is available at https://rnaprobe.genesilico.pl.

Accelerated cryo-EM-guided determination of three-dimensional RNA-only structures

The discovery and design of biologically important RNA molecules is outpacing three-dimensional structural characterization. Here, we demonstrate that cryo-electron microscopy can routinely resolve maps of RNA-only systems and that these maps enable subnanometer-resolution coordinate estimation when complemented with multidimensional chemical mapping and Rosetta DRRAFTER computational modeling. This hybrid 'Ribosolve' pipeline detects and falsifies homologies and conformational rearrangements in 11 previously unknown 119- to 338-nucleotide protein-free RNA structures: full-length Tetrahymena ribozyme, hc16 ligase with and without substrate, full-length Vibrio cholerae and Fusobacterium nucleatum glycine riboswitch aptamers with and without glycine, Mycobacterium SAM-IV riboswitch with and without S-adenosylmethionine, and the computer-designed ATP-TTR-3 aptamer with and without AMP. Simulation benchmarks, blind challenges, compensatory mutagenesis, cross-RNA homologies and internal controls demonstrate that Ribosolve can accurately resolve the global architectures of RNA molecules but does not resolve atomic details. These tests offer guidelines for making inferences in future RNA structural studies with similarly accelerated throughput.

Determination of RNA structural diversity and its role in HIV-1 RNA splicing

Human immunodeficiency virus-1 (HIV-1) is a retrovirus with a 10-kb single-stranded RNA genome. HIV-1 must express all of its gene products from the same primary transcript, which undergoes alternative splicing to produce diverse protein products, including structural proteins and regulatory factors^1,2. Despite the critical role of alternative splicing, the mechanisms driving splice-site choice are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements³. We use DMS-MaPseq to probe the structure of HIV-1 RNA in cells and develop an algorithm called Detection of RNA folding Ensembles using Expectation-Maximization (DREEM), which reveals alternative conformations assumed by the same RNA sequence. Contrary to previous models, which analyzed population averages⁴, our results reveal the widespread heterogeneous nature of HIV-1 RNA structure. In addition to confirming that in vitro characterized alternative structures for the HIV-1 Rev Responsive Element (RRE) exist in cells, we discover alternative conformations at critical splice sites that influence the ratio of transcript isoforms. Our simultaneous measurement of splicing and intracellular RNA structure provides evidence for the long-standing hypothesis^5–7 that RNA conformation heterogeneity regulates splice site usage and viral gene expression.

RNA base-pairing complexity in living cells visualized by correlated chemical probing

RNA structure and dynamics are critical to biological function. However, strategies for determining RNA structure in vivo are limited, with established chemical probing and newer duplex detection methods each having deficiencies. Here we convert the common reagent dimethyl sulfate into a useful probe of all 4 RNA nucleotides. Building on this advance, we introduce PAIR-MaP, which uses single-molecule correlated chemical probing to directly detect base-pairing interactions in cells. PAIR-MaP has superior resolution compared to alternative experiments, can resolve multiple sets of pairing interactions for structurally dynamic RNAs, and enables highly accurate structure modeling, including of RNAs containing multiple pseudoknots and extensively bound by proteins. Application of PAIR-MaP to human RNase MRP and 2 bacterial messenger RNA 5' untranslated regions reveals functionally important and complex structures undetected by prior analyses. PAIR-MaP is a powerful, experimentally concise, and broadly applicable strategy for directly visualizing RNA base pairs and dynamics in cells.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Conformational switch in the ribosomal protein S1 guides unfolding of structured RNAs for translation initiation

Initiation of bacterial translation requires that the ribosome-binding site in mRNAs adopts single-stranded conformations. In Gram-negative bacteria the ribosomal protein S1 (rS1) is a key player in resolving of structured elements in mRNAs. However, the exact mechanism of how rS1 unfolds persistent secondary structures in the translation initiation region (TIR) is still unknown. Here, we show by NMR spectroscopy that Vibrio vulnificus rS1 displays a unique architecture of its mRNA-binding domains, where domains D3 and D4 provide the mRNA-binding platform and cover the nucleotide binding length of the full-length rS1. D5 significantly increases rS1’s chaperone activity, although it displays structural heterogeneity both in isolation and in presence of the other domains, albeit to varying degrees. The heterogeneity is induced by the switch between the two equilibrium conformations and is triggered by an order-to-order transition of two mutually exclusive secondary structures (β-strand-to-α-helix) of the ‘AERERI’ sequence. The conformational switching is exploited for melting of structured 5′-UTR’s, as the conformational heterogeneity of D5 can compensate the entropic penalty of complex formation. Our data thus provides a detailed understanding of the intricate coupling of protein and RNA folding dynamics enabling translation initiation of structured mRNAs.

Evaluating riboswitch optimality

Riboswitches are RNA elements that recognize diverse chemical and biomolecular inputs, and transduce this recognition process to genetic, fluorescent, and other engineered outputs using RNA conformational changes. These systems are pervasive in cellular biology and are a promising biotechnology with applications in genetic regulation and biosensing. Here, we derive a simple expression bounding the activation ratio-the proportion of RNA in the active vs. inactive states-for both ON and OFF riboswitches that operate near thermodynamic equilibrium: 1+[I]/K_d^I, where [I] is the input ligand concentration and K_d^I is the intrinsic dissociation constant of the aptamer module toward the input ligand. A survey of published studies of natural and synthetic riboswitches confirms that the vast majority of empirically measured activation ratios have remained well below this thermodynamic limit. A few natural and synthetic riboswitches achieve activation ratios close to the limit, and these molecules highlight important principles for achieving high riboswitch performance. For several applications, including "light-up" fluorescent sensors and chemically-controlled CRISPR/Cas complexes, the thermodynamic limit has not yet been achieved, suggesting that current tools are operating at suboptimal efficiencies. Future riboswitch studies will benefit from comparing observed activation ratios to this simple expression for the optimal activation ratio. We present experimental and computational suggestions for how to make these quantitative comparisons and suggest new molecular mechanisms that may allow non-equilibrium riboswitches to surpass the derived limit.