Tertiary contacts control switching of the SAM-I riboswitch

Hierarchical Hydration Dynamics of RNA with Nano-Water-Pool at Its Core

Many functional RNAs fold into a compact, roughly globular shape by minimizing the electrostatic repulsion between their negatively charged phosphodiester backbone. The fold of such close, compact RNA architecture is often so designed that its outer surface and complex core both are predominately populated by phosphate groups loosely sequestering bases in the intermediate layers. A number of helical junctions maintain the RNA core and its nano-water-pool. While the folding of RNA is manifested by its counterion environment composed of mixed mono- and divalent salts, the concerted role of ion and water in maintaining an RNA fold is yet to be explored. In this work, detailed atomistic simulations of SAM-I and Add Adenine riboswitch aptamers, and subgenomic flavivirus RNA (sfRNA) have been performed in a physiological mixed mono- and divalent salt environment. All three RNA systems have compact folds with a core diameter of range 1-1.7 nm. The spatiotemporal heterogeneity of RNA hydration was probed in a layer-wise manner by distinguishing the core, the intermediate, and the outer layers. The layer-wise decomposition of hydrogen bonds and collective single-particle reorientational dynamics reveal a nonmonotonic relaxation pattern with the slowest relaxation observed at the intermediate layers that involves functionally important tertiary motifs. The slowness of this intermediate layer is attributed to two types of long-resident water molecules: (i) water from ion-hydration layers and (ii) structurally trapped water (distant from ions). The relaxation kinetics of the core and the surface water essentially exposed to the phosphate groups show well-separated time scales from the intermediate layers. In the slow intermediate layers, site-specific ions and water control the functional dynamics of important RNA motifs like kink-turn, observed in different structure-probing experiments. Most interestingly, we find that as the size of the RNA core increases (SAM1 core < sfRNAcore < Add adenine core), its hydration tends to show faster relaxation. The hierarchical hydration and the layer-wise base-phosphate composition uniquely portray the globular RNA to act like a soft vesicle with a quasi-dynamic nano-water-pool at its core.

Targeting SAM-I Riboswitch Using Antisense Oligonucleotide Technology for Inhibiting the Growth of Staphylococcus aureus and Listeria monocytogenes

With the discovery of antibiotics, a productive period of antibacterial drug innovation and application in healthcare systems and agriculture resulted in saving millions of lives. Unfortunately, the misusage of antibiotics led to the emergence of many resistant pathogenic strains. Some riboswitches have risen as promising targets for developing antibacterial drugs. Here, we describe the design and applications of the chimeric antisense oligonucleotide (ASO) as a novel antibacterial agent. The pVEC-ASO-1 consists of a cell-penetrating oligopeptide known as pVEC attached to an oligonucleotide part with modifications of the first and the second generations. This combination of modifications enables specific mRNA degradation under multiple turnover conditions via RNase H. The pVEC-ASO targets the S-adenosyl methionine (SAM)-I riboswitch found in the genome of many Gram-positive bacteria. The SAM-I riboswitch controls not only the biosynthesis but also the transport of SAM. We have established an antibiotic dosage of 700 nM (4.5 µg/mL) of pVEC-ASO that inhibits 80% of the growth of Staphylococcus aureus and Listeria monocytogenes. The pVEC-ASO-1 does not show any toxicity in the human cell line at MIC80's concentration. We have proven that the SAM-I riboswitch is a suitable target for antibacterial drug development based on ASO. The approach is rational and easily adapted to other bacterial RNA targets.

Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (Q_TC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Exploring the energy landscape of a SAM-I riboswitch

SAM-I riboswitches regulate gene expression through transcription termination upon binding a S-adenosyl-L-methionine (SAM) ligand. In previous work, we characterized the conformational energy landscape of the full-length Bacillus subtilis yitJ SAM-I riboswitch as a function of Mg²⁺ and SAM ligand concentrations. Here, we have extended this work with measurements on a structurally similar ligand, S-adenosyl-l-homocysteine (SAH), which has, however, a much lower binding affinity. Using single-molecule Förster resonance energy transfer (smFRET) microscopy and hidden Markov modeling (HMM) analysis, we identified major conformations and determined their fractional populations and dynamics. At high Mg²⁺ concentration, FRET analysis yielded four distinct conformations, which we assigned to two terminator and two antiterminator states. In the same solvent, but with SAM added at saturating concentrations, four states persisted, although their populations, lifetimes and interconversion dynamics changed. In the presence of SAH instead of SAM, HMM revealed again four well-populated states and, in addition, a weakly populated ‘hub’ state that appears to mediate conformational transitions between three of the other states. Our data show pronounced and specific effects of the SAM and SAH ligands on the RNA conformational energy landscape. Interestingly, both SAM and SAH shifted the fractional populations toward terminator folds, but only gradually, so the effect cannot explain the switching action. Instead, we propose that the noticeably accelerated dynamics of interconversion between terminator and antiterminator states upon SAM binding may be essential for control of transcription.

Chelated Magnesium Logic Gate Regulates Riboswitch Pseudoknot Formation

Magnesium plays a critical role in the structure, dynamics, and function of RNA. The precise microscopic effect of chelated magnesium on RNA structure is yet to be explored. Magnesium is known to act through its diffuse cloud around RNA, through the outer sphere (water-mediated), inner sphere, and often chelated ion-mediated interactions. A mechanism is proposed for the role of experimentally discovered site-specific chelated magnesium ions on the conformational dynamics of SAM-I riboswitch aptamers in bacteria. This mechanism is observed with atomistic simulations performed in a physiological mixed salt environment at a high temperature. The simulations were validated with phosphorothioate interference mapping experiments that help to identify crucial inner-sphere Mg²⁺ sites prescribing an appropriate initial distribution of inner- and outer-sphere magnesium ions to maintain a physiological ion concentration of monovalent and divalent salts. A concerted role of two chelated magnesium ions is newly discovered since the presence of both supports the formation of the pseudoknot. This constitutes a logical AND gate. The absence of any of these magnesium ions instigates the dissociation of long-range pseudoknot interaction exposing the inner core of the RNA. A base triple is the epicenter of the magnesium chelation effect. It allosterically controls RNA pseudoknot by bolstering the direct effect of magnesium chelation in protecting the functional fold of RNA to control ON and OFF transcription switching.

Zinc-finger protein CNBP alters the 3-D structure of lncRNA Braveheart in solution

Long non-coding RNAs (lncRNAs) constitute a significant fraction of the transcriptome, playing important roles in development and disease. However, our understanding of structure-function relationships for this emerging class of RNAs has been limited to secondary structures. Here, we report the 3-D atomistic structural study of epigenetic lncRNA, Braveheart (Bvht), and its complex with CNBP (Cellular Nucleic acid Binding Protein). Using small angle X-ray scattering (SAXS), we elucidate the ensemble of Bvht RNA conformations in solution, revealing that Bvht lncRNA has a well-defined, albeit flexible 3-D structure that is remodeled upon CNBP binding. Our study suggests that CNBP binding requires multiple domains of Bvht and the RHT/AGIL RNA motif. We show that RHT/AGIL, previously shown to interact with CNBP, contains a highly flexible loop surrounded by more ordered helices. As one of the largest RNA-only 3-D studies, the work lays the foundation for future structural studies of lncRNA-protein complexes.

Magnesium controls aptamer-expression platform switching in the SAM-I riboswitch

Investigations of most riboswitches remain confined to the ligand-binding aptamer domain. However, during the riboswitch mediated transcription regulation process, the aptamer domain and the expression platform compete for a shared strand. If the expression platform dominates, an anti-terminator helix is formed, and the transcription process is active (ON state). When the aptamer dominates, transcription is terminated (OFF state). Here, we use an expression platform switching experimental assay and structure-based electrostatic simulations to investigate this ON-OFF transition of the full length SAM-I riboswitch and its magnesium concentration dependence. Interestingly, we find the ratio of the OFF population to the ON population to vary non-monotonically as magnesium concentration increases. Upon addition of magnesium, the aptamer domain pre-organizes, populating the OFF state, but only up to an intermediate magnesium concentration level. Higher magnesium concentration preferentially stabilizes the anti-terminator helix, populating the ON state, relatively destabilizing the OFF state. Magnesium mediated aptamer-expression platform domain closure explains this relative destabilization of the OFF state at higher magnesium concentration. Our study reveals the functional potential of magnesium in controlling transcription of its downstream genes and underscores the importance of a narrow concentration regime near the physiological magnesium concentration ranges, striking a balance between the OFF and ON states in bacterial gene regulation.

Light-activated chemical probing of nucleobase solvent accessibility inside cells

The discovery of functional RNAs that are critical for normal and disease physiology continues to expand at a breakneck pace. Many RNA functions are controlled by the formation of specific structures, and an understanding of each structural component is necessary to elucidate its function. Measuring solvent accessibility intracellularly with experimental ease is an unmet need in the field. Here, we present a novel method for probing nucleobase solvent accessibility, Light Activated Structural Examination of RNA (LASER). LASER depends on light activation of a small molecule, nicotinoyl azide (NAz), to measure solvent accessibility of purine nucleobases. In vitro, this technique accurately monitors solvent accessibility and identifies rapid structural changes resulting from ligand binding in a metabolite-responsive RNA. LASER probing can further identify cellular RNA-protein interactions and unique intracellular RNA structures. Our photoactivation technique provides an adaptable framework to structurally characterize solvent accessibility of RNA in many environments.

Genetic regulation mechanism of the yjdF riboswitch

The yjdF riboswitch resides in potential 5' UTRs of homologues of protein-coding gene yjdF in Firmicutes. Unlike other 30 riboswitch classes previously validated, this riboswitch class, can sense and bind to a broad collection of azaaromatic ligands. Among these compounds, some do activate production of yjdF protein driven by the riboswitch, while others are out of riboswitch-mediated modulation possibly because of the toxicity at high ligand concentrations. By incorporating the structures with pseudoknots and ligand binding kinetics into the co-transcriptional folding theory, we theoretically studied the co-transcriptional folding behaviors of the yjdF riboswitch from Bacillus subtilis at different transcription conditions. Like most riboswitches, the yjdF riboswitch can quickly fold into the aptamer structure without any trapped states during the transcription process. After the aptamer structure is formed, the riboswitch shows two main co-transcriptional folding pathways: aptamer→ON state→OFF state and aptamer → the ligand bound aptamer → the ligand bound ON state. Our results suggested that this translational riboswitch is coupled with the transcription process to exert its biological function and it is kinetically controlled. The threshold concentration for the ligand to activate the riboswitch depends on the transcription rate and the association rate of the ligand binding.

Ligand recognition and helical stacking formation are intimately linked in the SAM-I riboswitch regulatory mechanism

Riboswitches are noncoding mRNA elements that control gene expression by altering their structure upon metabolite binding. Although riboswitch crystal structures provide detailed information about RNA-ligand interactions, little knowledge has been gathered to understand how riboswitches modulate gene expression. Here, we study the molecular recognition mechanism of the S-adenosylmethionine SAM-I riboswitch by characterizing the formation of a helical stacking interaction involving the ligand-binding process. We show that ligand binding is intimately linked to the formation of the helical stacking, which is dependent on the presence of three conserved purine residues that are flanked by stacked helices. We also find that these residues are important for the formation of a crucial long-range base pair formed upon SAM binding. Together, our results lend strong support to a critical role for helical stacking in the folding pathway and suggest a particularly important function in the formation of the long-range base pair.

Co-Transcriptional Folding and Regulation Mechanisms of Riboswitches

Riboswitches are genetic control elements within non-coding regions of mRNA. These self-regulatory elements have been found to sense a range of small metabolites, ions, and other physical signals to exert regulatory control of transcription, translation, and splicing. To date, more than a dozen riboswitch classes have been characterized that vary widely in size and secondary structure. Extensive experiments and theoretical studies have made great strides in understanding the general structures, genetic mechanisms, and regulatory activities of individual riboswitches. As the ligand-dependent co-transcriptional folding and unfolding dynamics of riboswitches are the key determinant of gene expression, it is important to investigate the thermodynamics and kinetics of riboswitches both in the presence and absence of metabolites under the transcription. This review will provide a brief summary of the studies about the regulation mechanisms of the pbuE, S_MK, yitJ, and metF riboswitches based on the ligand-dependent co-transcriptional folding of the riboswitches.

A magnesium-induced triplex pre-organizes the SAM-II riboswitch

Our 13C- and 1H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.

Linking aptamer-ligand binding and expression platform folding in riboswitches: prospects for mechanistic modeling and design

The power of riboswitches in regulation of bacterial metabolism derives from coupling of two characteristics: recognition and folding. Riboswitches contain aptamers, which function as biosensors. Upon detection of the signaling molecule, the riboswitch transduces the signal into a genetic decision. The genetic decision is coupled to refolding of the expression platform, which is distinct from, although overlapping with, the aptamer. Early biophysical studies of riboswitches focused on recognition of the ligand by the aptamer-an important consideration for drug design. A mechanistic understanding of ligand-induced riboswitch RNA folding can further enhance riboswitch ligand design, and inform efforts to tune and engineer riboswitches with novel properties. X-ray structures of aptamer/ligand complexes point to mechanisms through which the ligand brings together distal strand segments to form a P1 helix. Transcriptional riboswitches must detect the ligand and form this P1 helix within the timescale of transcription. Depending on the cell's metabolic state and cellular environmental conditions, the folding and genetic outcome may therefore be affected by kinetics of ligand binding, RNA folding, and transcriptional pausing, among other factors. Although some studies of isolated riboswitch aptamers found homogeneous, prefolded conformations, experimental, and theoretical studies point to functional and structural heterogeneity for nascent transcripts. Recently it has been shown that some riboswitch segments, containing the aptamer and partial expression platforms, can form binding-competent conformers that incorporate an incomplete aptamer secondary structure. Consideration of the free energy landscape for riboswitch RNA folding suggests models for how these conformers may act as transition states-facilitating rapid, ligand-mediated aptamer folding.

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.

Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT

Integration and calibration of molecular dynamics simulations with experimental data remain a challenging endeavor. We have developed a novel method to integrate chemical probing experiments with molecular simulations of RNA molecules by using a native structure-based model. Selective 2'-hydroxyl acylation by primer extension (SHAPE) characterizes the mobility of each residue in the RNA. Our method, SHAPE-FIT, automatically optimizes the potential parameters of the force field according to measured reactivities from SHAPE. The optimized parameter set allows simulations of dynamics highly consistent with SHAPE probing experiments. Such atomistic simulations, thoroughly grounded in experiment, can open a new window on RNA structure-function relations.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Fluorescence tools to investigate riboswitch structural dynamics

Riboswitches are novel regulatory elements that respond to cellular metabolites to control gene expression. They are constituted of highly conserved domains that have evolved to recognize specific metabolites. Such domains, so-called aptamers, are folded into intricate structures to enable metabolite recognition. Over the years, the development of ensemble and single-molecule fluorescence techniques has allowed to probe most of the mechanistic aspects of aptamer folding and ligand binding. In this review, we summarize the current fluorescence toolkit available to study riboswitch structural dynamics. We fist describe those methods based on fluorescent nucleotide analogues, mostly 2-aminopurine (2AP), to investigate short-range conformational changes, including some key steady-state and time-resolved examples that exemplify the versatility of fluorescent analogues as structural probes. The study of long-range structural changes by Förster resonance energy transfer (FRET) is mostly discussed in the context of single-molecule studies, including some recent developments based on the combination of single-molecule FRET techniques with controlled chemical denaturation methods. This article is part of a Special Issue entitled: Riboswitches.

Common themes and differences in SAM recognition among SAM riboswitches

The recent discovery of short cis-acting RNA elements termed riboswitches has caused a paradigm shift in our understanding of genetic regulatory mechanisms. The three distinct superfamilies of S-adenosyl-l-methionine (SAM) riboswitches are the most commonly found riboswitch classes in nature. These RNAs represent three independent evolutionary solutions to achieve specific SAM recognition. This review summarizes research on 1) modes of gene regulatory mechanisms, 2) common themes and differences in ligand recognition, and 3) ligand-induced conformational dynamics among SAM riboswitch families. The body of work on the SAM riboswitch families constitutes a useful primer to the topic of gene regulatory RNAs as a whole. This article is part of a Special Issue entitled: Riboswitches.

Structure-guided design of fluorescent S-adenosylmethionine analogs for a high-throughput screen to target SAM-I riboswitch RNAs

Many classes of S-adenosylmethionine (SAM)-binding RNAs and proteins are of interest as potential drug targets in diverse therapeutic areas, from infectious diseases to cancer. In the former case, the SAM-I riboswitch is an attractive target because this structured RNA element is found only in bacterial mRNAs and regulates multiple genes in several human pathogens. Here, we describe the synthesis of stable and fluorescent analogs of SAM in which the fluorophore is introduced through a functionalizable linker to the ribose. A Cy5-labeled SAM analog was shown to bind several SAM-I riboswitches via in-line probing and fluorescence polarization assays, including one from Staphylococcus aureus that controls the expression of SAM synthetase in this organism. A fluorescent ligand displacement assay was developed and validated for high-throughput screening of compounds to target the SAM-I riboswitch class.

Sizing up long non-coding RNAs: do lncRNAs have secondary and tertiary structure?

Long noncoding RNAs (lncRNAs) play a key role in many important areas of epigenetics, stem cell biology, cancer, signaling and brain function. This emerging class of RNAs constitutes a large fraction of the transcriptome, with thousands of new lncRNAs reported each year. The molecular mechanisms of these RNAs are not well understood. Currently, very little structural data exist. We review the available lncRNA sequence and secondary structure data. Since almost no tertiary information is available for lncRNAs, we review crystallographic structures for other RNA systems and discuss the possibilities for lncRNAs in the context of existing constraints.

The impact of a ligand binding on strand migration in the SAM-I riboswitch

Riboswitches sense cellular concentrations of small molecules and use this information to adjust synthesis rates of related metabolites. Riboswitches include an aptamer domain to detect the ligand and an expression platform to control gene expression. Previous structural studies of riboswitches largely focused on aptamers, truncating the expression domain to suppress conformational switching. To link ligand/aptamer binding to conformational switching, we constructed models of an S-adenosyl methionine (SAM)-I riboswitch RNA segment incorporating elements of the expression platform, allowing formation of an antiterminator (AT) helix. Using Anton, a computer specially developed for long timescale Molecular Dynamics (MD), we simulated an extended (three microseconds) MD trajectory with SAM bound to a modeled riboswitch RNA segment. Remarkably, we observed a strand migration, converting three base pairs from an antiterminator (AT) helix, characteristic of the transcription ON state, to a P1 helix, characteristic of the OFF state. This conformational switching towards the OFF state is observed only in the presence of SAM. Among seven extended trajectories with three starting structures, the presence of SAM enhances the trend towards the OFF state for two out of three starting structures tested. Our simulation provides a visual demonstration of how a small molecule (<500 MW) binding to a limited surface can trigger a large scale conformational rearrangement in a 40 kDa RNA by perturbing the Free Energy Landscape. Such a mechanism can explain minimal requirements for SAM binding and transcription termination for SAM-I riboswitches previously reported experimentally.

Fluorescence monitoring of riboswitch transcription regulation using a dual molecular beacon assay

Riboswitches are mRNA elements that specifically bind cellular metabolites and control gene expression by modifying their structure. As riboswitches often control essential genes in pathogenic bacteria, riboswitches have been proposed as new targets for antibiotics. High-throughput screening provides a powerful approach to identify riboswitch ligand analogs that could act as powerful antibacterial drugs. Biochemical assays have already been used to find riboswitch-binding analogs, but those methods do take into account the transcriptional context for riboswitch regulation. As the importance of co-transcriptional ligand binding has been shown for several riboswitches, it is vital to develop an assay that screens riboswitch-binding analogs during the transcriptional process. Here, we describe the development of a dual molecular beacon system monitoring the transcriptional regulation activity of the Bacillus subtilis pbuE adenine riboswitch. This system relies on two molecular beacons that enable the monitoring of transcription efficiency, as well as the regulatory activity of the riboswitch. Different analogs were tested using our system, and a good correlation was observed between riboswitch activity and reported metabolite affinities. This method is specific, reliable and could be applied at the high-throughput level for the identification of new potential antibiotics targeting any riboswitch-regulating gene expression at the mRNA level.

Rise of the RNA machines: exploring the structure of long non-coding RNAs

Novel, profound and unexpected roles of long non-coding RNAs (lncRNAs) are emerging in critical aspects of gene regulation. Thousands of lncRNAs have been recently discovered in a wide range of mammalian systems, related to development, epigenetics, cancer, brain function and hereditary disease. The structural biology of these lncRNAs presents a brave new RNA world, which may contain a diverse zoo of new architectures and mechanisms. While structural studies of lncRNAs are in their infancy, we describe existing structural data for lncRNAs, as well as crystallographic studies of other RNA machines and their implications for lncRNAs. We also discuss the importance of dynamics in RNA machine mechanism. Determining commonalities between lncRNA systems will help elucidate the evolution and mechanistic role of lncRNAs in disease, creating a structural framework necessary to pursue lncRNA-based therapeutics.

The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch

Riboswitch operation involves the complex interplay between the aptamer domain and the expression platform. During transcription, these two domains compete against each other for shared sequence. In this study, we explore the cooperative effects of ligand binding and Magnesium interactions in the SAM-I riboswitch in the context of aptamer collapse and anti-terminator formation. Overall, our studies show the apo-aptamer acts as (i) a pre-organized aptamer competent to bind ligand and undergo structural collapse and (ii) a conformation that is more accessible to anti-terminator formation. We show that both Mg(2+) ions and SAM are required for a collapse transition to occur. We then use competition between the aptamer and expression platform for shared sequence to characterize the stability of the collapsed aptamer. We find that SAM and Mg(2+) interactions in the aptamer are highly cooperative in maintaining switch polarity (i.e. aptamer 'off-state' versus anti-terminator 'on-state'). We further show that the aptamer off-state is preferentially stabilized by Mg(2+) and similar divalent ions. Furthermore, the functional switching assay was used to select for phosphorothioate interference, and identifies potential magnesium chelation sites while characterizing their coordinated role with SAM in aptamer stabilization. In addition, we find that Mg(2+) interactions with the apo-aptamer are required for the full formation of the anti-terminator structure, and that higher concentrations of Mg(2+) (>4 mM) shift the equilibrium toward the anti-terminator on-state even in the presence of SAM.

Magnesium fluctuations modulate RNA dynamics in the SAM-I riboswitch

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.

Biomolecular dynamics: order-disorder transitions and energy landscapes

While the energy landscape theory of protein folding is now a widely accepted view for understanding how relatively weak molecular interactions lead to rapid and cooperative protein folding, such a framework must be extended to describe the large-scale functional motions observed in molecular machines. In this review, we discuss (1) the development of the energy landscape theory of biomolecular folding, (2) recent advances toward establishing a consistent understanding of folding and function and (3) emerging themes in the functional motions of enzymes, biomolecular motors and other biomolecular machines. Recent theoretical, computational and experimental lines of investigation have provided a very dynamic picture of biomolecular motion. In contrast to earlier ideas, where molecular machines were thought to function similarly to macroscopic machines, with rigid components that move along a few degrees of freedom in a deterministic fashion, biomolecular complexes are only marginally stable. Since the stabilizing contribution of each atomic interaction is on the order of the thermal fluctuations in solution, the rigid body description of molecular function must be revisited. An emerging theme is that functional motions encompass order-disorder transitions and structural flexibility provides significant contributions to the free energy. In this review, we describe the biological importance of order-disorder transitions and discuss the statistical-mechanical foundation of theoretical approaches that can characterize such transitions.

Basis for ligand discrimination between ON and OFF state riboswitch conformations: the case of the SAM-I riboswitch

Riboswitches are RNA elements that bind to effector ligands and control gene expression. Most consist of two domains. S-Adenosyl Methionine (SAM) binds the aptamer domain of the SAM-I riboswitch and induces conformational changes in the expression domain to form an intrinsic terminator (transcription OFF state). Without SAM the riboswitch forms the transcription ON state, allowing read-through transcription. The mechanistic link between the SAM/aptamer recognition event and subsequent secondary structure rearrangement by the riboswitch is unclear. We probed for those structural features of the Bacillus subtilis yitJ SAM-I riboswitch responsible for discrimination between the ON and OFF states by SAM. We designed SAM-I riboswitch RNA segments forming "hybrid" structures of the ON and OFF states. The choice of segment constrains the formation of a partial P1 helix, characteristic of the OFF state, together with a partial antiterminator (AT) helix, characteristic of the ON state. For most choices of P1 vs. AT helix lengths, SAM binds with micromolar affinity according to equilibrium dialysis. Mutational analysis and in-line probing confirm that the mode of SAM binding by hybrid structures is similar to that of the aptamer. Altogether, binding measurements and in-line probing are consistent with the hypothesis that when SAM is present, stacking interactions with the AT helix stabilize a partially formed P1 helix in the hybrids. Molecular modeling indicates that continuous stacking between the P1 and the AT helices is plausible with SAM bound. Our findings raise the possibility that conformational intermediates may play a role in ligand-induced aptamer folding.

Molecular recognition and function of riboswitches

Regulatory mRNAs elements termed riboswitches respond to elevated concentrations of cellular metabolites by modulating expression of associated genes. Riboswitches attain their high metabolite selectivity by capitalizing on the intrinsic tertiary structures of their sensor domains. Over the years, riboswitch structure and folding have been amongst the most researched topics in the RNA field. Most recently, novel structures of single-ligand and cooperative double-ligand sensors have broadened our knowledge of architectural and molecular recognition principles exploited by riboswitches. The structural information has been complemented by extensive folding studies, which have provided several important clues on the formation of ligand-competent conformations and mechanisms of ligand discrimination. These studies have greatly improved our understanding of molecular events in riboswitch-mediated gene expression control and provided the molecular basis for intervention into riboswitch-controlled genetic circuits.

Conformational heterogeneity of the SAM-I riboswitch transcriptional ON state: a chaperone-like role for S-adenosyl methionine

Riboswitches are promising targets for the design of novel antibiotics and engineering of portable genetic regulatory elements. There is evidence that variability in riboswitch properties allows tuning of expression for genes involved in different stages of biosynthetic pathways by mechanisms that are not currently understood. Here, we explore the mechanism for tuning of S-adenosyl methionine (SAM)-I riboswitch folding. Most SAM-I riboswitches function at the transcriptional level by sensing the cognate ligand SAM. SAM-I riboswitches orchestrate the biosynthetic pathways of cysteine, methionine, SAM, and so forth. We use base-pair probability predictions to examine the secondary-structure folding landscape of several SAM-I riboswitch sequences. We predict different folding behaviors for different SAM-I riboswitch sequences. We identify several "decoy" base-pairing interactions involving 5' riboswitch residues that can compete with the formation of a P1 helix, a component of the ligand-bound "transcription OFF" state, in the absence of SAM. We hypothesize that blockage of these interactions through SAM contacts contributes to stabilization of the OFF state in the presence of ligand. We also probe folding patterns for a SAM-I riboswitch RNA using constructs with different 3' truncation points experimentally. Folding was monitored through fluorescence, susceptibility to base-catalyzed cleavage, nuclear magnetic resonance, and indirectly through SAM binding. We identify key decision windows at which SAM can affect the folding pathway towards the OFF state. The presence of decoy conformations and differential sensitivities to SAM at different transcript lengths is crucial for SAM-I riboswitches to modulate gene expression in the context of global cellular metabolism.

Dynamic energy landscapes of riboswitches help interpret conformational rearrangements and function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way that sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we mimic transcription by computing energy landscapes which describe accessible secondary structures for a range of sequence lengths. Consistent with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: Low-barrier landscape with an ensemble of conformations in equilibrium before encountering a substrate; barrier-free landscape with a dominant “downhill” pathway to the minimum free energy structure; and barrier-dominated landscape with two isolated conformational states with different functions. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved between riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design.

RNA tertiary interactions in a riboswitch stabilize the structure of a kink turn

The kink turn is a widespread RNA motif that introduces an acute kink into the axis of duplex RNA, typically comprising a bulge followed by a G⋅A and A⋅G pairs. The kinked conformation is stabilized by metal ions, or the binding of proteins including L7Ae. We now demonstrate a third mechanism for the stabilization of k-turn structure, involving tertiary interactions within a larger RNA structure. The SAM-I riboswitch contains an essential standard k-turn sequence that kinks a helix so that its terminal loop can make a long-range interaction. We find that some sequence variations in the k-turn within the riboswitch do not prevent SAM binding, despite preventing the folding of the k-turn in isolation. Furthermore, two crystal structures show that the sequence-variant k-turns are conventionally folded within the riboswitch. This study shows that the folded structure of the k-turn can be stabilized by tertiary interactions within a larger RNA structure.

Molecular insights into the ligand-controlled organization of the SAM-I riboswitch

S-adenosylmethionine (SAM) riboswitches are widespread in bacteria, and up to five different SAM riboswitch families have been reported, highlighting the relevance of SAM regulation. On the basis of crystallographic and biochemical data, it has been postulated, but never demonstrated, that ligand recognition by SAM riboswitches involves key conformational changes in the RNA architecture. We show here that the aptamer follows a two-step hierarchical folding selectively induced by metal ions and ligand binding, each of them leading to the formation of one of the two helical stacks observed in the crystal structure. Moreover, we find that the anti-antiterminator P1 stem is rotated along its helical axis upon ligand binding, a mechanistic feature that could be common to other riboswitches. We also show that the nonconserved P4 helical domain is used as an auxiliary element to enhance the ligand-binding affinity. This work provides the first comprehensive characterization, to our knowledge, of a ligand-controlled riboswitch folding pathway.