Structure of the SAM-II riboswitch bound to S-adenosylmethionine

Targeting RNA with small molecules, from RNA structures to precision medicines: IUPHAR review: 40

RNA plays important roles in regulating both health and disease biology in all kingdoms of life. Notably, RNA can form intricate three-dimensional structures, and their biological functions are dependent on these structures. Targeting the structured regions of RNA with small molecules has gained increasing attention over the past decade, because it provides both chemical probes to study fundamental biology processes and lead medicines for diseases with unmet medical needs. Recent advances in RNA structure prediction and determination and RNA biology have accelerated the rational design and development of RNA-targeted small molecules to modulate disease pathology. However, challenges remain in advancing RNA-targeted small molecules towards clinical applications. This review summarizes strategies to study RNA structures, to identify small molecules recognizing these structures, and to augment the functionality of RNA-binding small molecules. We focus on recent advances in developing RNA-targeted small molecules as potential therapeutics in a variety of diseases, encompassing different modes of actions and targeting strategies. Furthermore, we present the current gaps between early-stage discovery of RNA-binding small molecules and their clinical applications, as well as a roadmap to overcome these challenges in the near future.

Linker-Mediated Inactivation of the SAM-II Domain in the Tandem SAM-II/SAM-V Riboswitch

Tandem SAM-II/SAM-V riboswitch belongs to a class of riboswitches found in the marine bacterium 'Candidatus Pelagibacter ubique'. Previous studies have demonstrated that these riboswitches have the potential for digital modulation of gene expression at both the transcriptional and translational levels. In this study, we investigate the conformational changes in the tandem SAM-II/SAM-V riboswitch binding to S-adenosylmethionine (SAM) using selective 2'-hydroxyl acylation analyzed by the primer extension (SHAPE) assay, small-angle X-ray scattering (SAXS), and oligos depressing probing. Our findings reveal that the linker between SAM-II/SAM-V aptamers blocks the SAM response of the SAM-II domain. This result proposes a new mechanism for gene expression regulation, where the ligand-binding functions of tandem riboswitches can be selectively masked or released through a linker.

Metastasis-Associated Lung Adenocarcinoma Transcript 1 (MALAT1) lncRNA Conformational Dynamics in Complex with RNA-Binding Protein with Serine-Rich Domain 1 (RNPS1) in the Pan-cancer Splicing and Gene Expression

Alternative splicing events increase the transcriptomic and proteomic complexity in cancers. Overexpression of metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), a highly conserved lncRNA, is widely known to promote cancer development, one mechanism for which may be through the regulation of alternative splicing and, thereby, gene expression. Its regulatory interactions with proteins have been a subject of much interest, yet little research has been carried out on the mechanisms adopted. It has been observed that MALAT1 binds to RNA-binding protein with serine-rich domain 1 (RNPS1), being colocalized in the nuclear speckles, and together, these two binding partners may regulate alternative splicing. Upregulated RNPS1 is predicted to play a key role in the pan-cancer development. Experimental tertiary structure of full-length MALAT1 is currently lacking despite the availability of the 3D structure of 3' expression and nuclear retention element. We hypothesize that the computationally modeled tertiary structures of the specific binding motifs in the M-region, E-region, and full-length structures of MALAT1 may adopt a modular structure and bind to the RNPS1 loop region of RS/P domain involved in exon skipping, interacting in a manner fully consistent with the biochemical experiments. Extensive observations using the powerful molecular dynamics (MD) simulations of MALAT1 regions bound to RNPS1 suggested that all three regions form interactive, yet stable complexes. The ranking of the MM-GBSA- and MM-PBSA-derived binding free energies between these complexes corroborated well in the MD simulations and experiments. Energy decomposition analyses suggested that arginines in the RNPS1 protein are among the major contributors toward the binding free energies as calculated by MM-GBSA present in the Amber package; while among the nucleotides, the major contributors were nucleotides with G and A nucleobases, with more contributory effect in comparison to arginines, across the bound M-region, E-region, and full-length MALAT1. This suggests that specific purines play a greater role in the complex formation, in a loop-specific manner, and the more proactive approach in complexation tilts toward MALAT1. To the best of our knowledge, our studies are the first studies taking a unique approach, utilizing the binding motifs to deduce a tertiary structure of MALAT1, toward our understanding of the lncRNA-protein interactions, stability, and binding on a structural basis. The therapeutic implications of targeting this complex formation to regulate splicing and hence, oncogenesis, is further envisaged.

Antisense RNA C9orf72 hexanucleotide repeat associated with amyotrophic lateral sclerosis and frontotemporal dementia forms a triplex-like structure and binds small synthetic ligand

The abnormal expansion of GGGGCC/GGCCCC hexanucleotide repeats (HR) in C9orf72 is associated with amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Structural polymorphisms of HR result in the multifactorial pathomechanism of ALS/FTD. Consequently, many ongoing studies are focused at developing therapies targeting pathogenic HR RNA. One of them involves small molecules blocking sequestration of important proteins, preventing formation of toxic nuclear foci. However, rational design of potential therapeutics is hindered by limited number of structural studies of RNA-ligand complexes. We determined the crystal structure of antisense HR RNA in complex with ANP77 ligand (1.1 Å resolution) and in the free form (0.92 and 1.5 Å resolution). HR RNA folds into a triplex structure composed of four RNA chains. ANP77 interacted with two neighboring single-stranded cytosines to form pseudo-canonical base pairs by adopting sandwich-like conformation and adjusting the position of its naphthyridine units to the helical twist of the RNA. In the unliganded structure, the cytosines formed a peculiar triplex i-motif, assembled by trans C•C+ pair and a third cytosine located at the Hoogsteen edge of the C•C+ pair. These results extend our knowledge of the structural polymorphisms of HR and can be used for rational design of small molecules targeting disease-related RNAs.

Bifidobacterium bifidum SAM-VI Riboswitch Conformation Change Requires Peripheral Helix Formation

The Bifidobacterium bifidum SAM-VI riboswitch undergoes dynamic conformational changes that modulate downstream gene expression. Traditional structural methods such as crystallography capture the bound conformation at high resolution, and additional efforts would reveal details from the dynamic transition. Here, we revealed a transcription-dependent conformation model for Bifidobacterium bifidum SAM-VI riboswitch. In this study, we combine small-angle X-ray scattering, chemical probing, and isothermal titration calorimetry to unveil the ligand-binding properties and conformational changes of the Bifidobacterium bifidum SAM-VI riboswitch and its variants. Our results suggest that the SAM-VI riboswitch contains a pre-organized ligand-binding pocket and stabilizes into the bound conformation upon binding to SAM. Whether the P1 stem formed and variations in length critically influence the conformational dynamics of the SAM-VI riboswitch. Our study provides the basis for artificially engineering the riboswitch by manipulating its peripheral sequences without modifying the SAM-binding core.

Biophysical Investigation of RNA ⋅ DNA : DNA Triple Helix and RNA : DNA Heteroduplex Formation by the lncRNAs MEG3 and Fendrr

Long non-coding RNAs (lncRNAs) are important regulators of gene expression and can associate with DNA as RNA : DNA heteroduplexes or RNA ⋅ DNA : DNA triple helix structures. Here, we review in vitro biochemical and biophysical experiments including electromobility shift assays (EMSA), circular dichroism (CD) spectroscopy, thermal melting analysis, microscale thermophoresis (MST), single-molecule Förster resonance energy transfer (smFRET) and nuclear magnetic resonance (NMR) spectroscopy to investigate RNA ⋅ DNA : DNA triple helix and RNA : DNA heteroduplex formation. We present the investigations of the antiparallel triplex-forming lncRNA MEG3 targeting the gene TGFB2 and the parallel triplex-forming lncRNA Fendrr with its target gene Emp2. The thermodynamic properties of these oligonucleotides lead to concentration-dependent heterogeneous mixtures, where a DNA duplex, an RNA : DNA heteroduplex and an RNA ⋅ DNA : DNA triplex coexist and their relative populations are modulated in a temperature-dependent manner. The in vitro data provide a reliable readout of triplex structures, as RNA ⋅ DNA : DNA triplexes show distinct features compared to DNA duplexes and RNA : DNA heteroduplexes. Our experimental results can be used to validate computationally predicted triple helix formation between novel disease-relevant lncRNAs and their DNA target genes.

S-adenosyl-l-methionine interaction signatures in methyltransferases

The switching on or off of methylation, a change from a normal methylation to hyper or hypo methylation is implicated in many diseases that include cancers, infectious, neurodegenerative diseases and others. Methyltransferases are one of the most sought targets that have diversified for the methylation of a variety of substrates. However, without S-adenosyl-l-methionine (SAM), the universal methyl donor, the majority of the methyltransferases remain functionally inactive. In this article, we did a comprehensive analysis of all available SAM-receptor crystal structures at atom, moiety and structure levels to gain deeper insights into the structure and function of SAM. SAM demonstrated flexibility in binding to a variety of receptors irrespective of the size of the binding pockets. Further analysis of the binding pockets resulted in all SAM conformations clustering into four natural shapes. The conserved interaction analysis provides an unambiguous orientation of SAM binding to receptors which has been elusive till now. SAM peptide moiety (SPM) and SAM nucleobase moiety (SNM) show up to 89% interactions with receptors whereas only 11% interactions with SAM ribose moiety (SRM). It is found that SPM and SNM terminal atoms anchor to the highly conserved receptor subsites creating a workbench for catalysis. It is seen that every interacting atom and its position is crucial in the methyl transfer phenomenon. A very unique observation is that the methyl group of SAM does not have even one interaction with the receptor. The deep insights gained help in the design and development of novel drugs against the methyltransferases.Communicated by Ramaswamy H. Sarma.

Flipping the script: Understanding riboswitches from an alternative perspective

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

Characterization and designing of an SAM riboswitch to establish a high-throughput screening platform for SAM overproduction in Saccharomyces cerevisiae

S-adenosyl- l-methionine (SAM) is a high-value compound widely used in the treatment of various diseases. SAM can be produced through fermentation, but further enhancing the microbial production of SAM requires novel high-throughput screening methods for rapid detection and screening of mutant libraries. In this work, an SAM-OFF riboswitch capable of responding to the SAM concentration was obtained and a high-throughput platform for screening SAM overproducers was established. SAM synthase was engineered by semirational design and directed evolution, which resulted in the SAM2S203F,W164R,T251S,Y285F,S365R mutant with almost twice higher catalytic activity than the parental enzyme. The best mutant was then introduced into Saccharomyces cerevisiae BY4741, and the resulting strain BSM8 produced a sevenfold higher SAM titer in shake-flask fermentation, reaching 1.25 g L-1 . This work provides a reference for designing biosensors to dynamically detect metabolite concentrations for high-throughput screening and the construction of effective microbial cell factories.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibility and the role of co-transcriptional folding in these processes.

Mechanism of Ligand Discrimination by the NMT1 Riboswitch

Riboswitches are conserved functional domains in mRNA that almost exclusively exist in bacteria. They regulate the biosynthesis and transport of amino acids and essential metabolites such as coenzymes, nucleobases, and their derivatives by specifically binding small molecules. Due to their ability to precisely discriminate between different cognate molecules as well as their common existence in bacteria, riboswitches have become potential antibacterial drug targets that could deliver urgently needed antibiotics with novel mechanisms of action. In this work, we report the recognition mechanisms of four oxidization products (XAN, AZA, UAC, and HPA) generated during purine degradation by an RNA motif termed the NMT1 riboswitch. Specifically, we investigated the physical interactions between the riboswitch and the oxidized metabolites by computing the changes in the free energy on mutating key nucleobases in the ligand binding pocket of the riboswitch. We discovered that the electrostatic interactions are central to ligand discrimination by this riboswitch. The relative binding free energies of the mutations further indicated that some of the mutations can also strengthen the binding affinities of the ligands (AZA, UAC, and HPA). These mechanistic details are also potentially relevant in the design of novel compounds targeting riboswitches.

Computational design and experimental verification of pseudoknotted ribozymes

The design of new RNA sequences that retain the function of a model RNA structure is a challenge in bioinformatics because of the structural complexity of these molecules. RNA can fold into its secondary and tertiary structures by forming stem-loops and pseudoknots. A pseudoknot is a set of base pairs between a region within a stem-loop and nucleotides outside of this stem-loop; this motif is very important for numerous functional structures. It is important for any computational design algorithm to take into account these interactions to give a reliable result for any structures that include pseudoknots. In our study, we experimentally validated synthetic ribozymes designed by Enzymer, which implements algorithms allowing for the design of pseudoknots. Enzymer is a program that uses an inverse folding approach to design pseudoknotted RNAs; we used it in this study to design two types of ribozymes. The ribozymes tested were the hammerhead and the glmS, which have a self-cleaving activity that allows them to liberate the new RNA genome copy during rolling-circle replication or to control the expression of the downstream genes, respectively. We demonstrated the efficiency of Enzymer by showing that the pseudoknotted hammerhead and glmS ribozymes sequences it designed were extensively modified compared to wild-type sequences and were still active.

Targeting Riboswitches with Beta-Axial-Substituted Cobalamins

RNA-targeting small-molecule therapeutics is an emerging field hindered by an incomplete understanding of the basic principles governing RNA-ligand interactions. One way to advance our knowledge in this area is to study model systems where these interactions are better understood, such as riboswitches. Riboswitches bind a wide array of small molecules with high affinity and selectivity, providing a wealth of information on how RNA recognizes ligands through diverse structures. The cobalamin-sensing riboswitch is a particularly useful model system, as similar sequences show highly specialized binding preferences for different biological forms of cobalamin. This riboswitch is also widely dispersed across bacteria and therefore holds strong potential as an antibiotic target. Many synthetic cobalamin forms have been developed for various purposes including therapeutics, but their interaction with cobalamin riboswitches is yet to be explored. In this study, we characterize the interactions of 11 cobalamin derivatives with three representative cobalamin riboswitches using in vitro binding experiments (both chemical footprinting and a fluorescence-based assay) and a cell-based reporter assay. The derivatives show productive interactions with two of the three riboswitches, demonstrating simultaneous plasticity and selectivity within these RNAs. The observed plasticity is partially achieved through a novel structural rearrangement within the ligand binding pocket, providing insight into how similar RNA structures can be targeted. As the derivatives also show in vivo functionality, they serve as several potential lead compounds for further drug development.

General Strategies for RNA X-ray Crystallography

An extremely small proportion of the X-ray crystal structures deposited in the Protein Data Bank are of RNA or RNA-protein complexes. This is due to three main obstacles to the successful determination of RNA structure: (1) low yields of pure, properly folded RNA; (2) difficulty creating crystal contacts due to low sequence diversity; and (3) limited methods for phasing. Various approaches have been developed to address these obstacles, such as native RNA purification, engineered crystallization modules, and incorporation of proteins to assist in phasing. In this review, we will discuss these strategies and provide examples of how they are used in practice.

Structural and Energetic Features of Base-Base Stacking Contacts in RNA

Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.

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.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

RNA as a Major-Groove Ligand: RNA-RNA and RNA-DNA Triplexes Formed by GAA and UUC or TTC Sequences

Friedreich's ataxia is associated with noncanonical nucleic acid structures that emerge when GAA:TTC repeats in the first intron of the FXN gene expand beyond a critical number of repeats. Specifically, the noncanonical repeats are associated with both triplexes and R-loops. Here, we present an in silico investigation of all possible triplexes that form by attaching a third RNA strand to an RNA:RNA or DNA:DNA duplex, complementing previous DNA-based triplex studies. For both new triplexes results are similar. For a pyridimine UUC⁺ third strand, the parallel orientation is stable while its antiparallel counterpart is unstable. For a neutral GAA third strand, the parallel conformation is stable. A protonated GA⁺A third strand is stable in both parallel and antiparallel orientations. We have also investigated Na⁺ and Mg²⁺ ion distributions around the triplexes. The presence of Mg²⁺ ions helps stabilize neutral, antiparallel GAA triplexes. These results (along with previous DNA-based studies) allow for the emergence of a complete picture of the stability and structural characteristics of triplexes based on the GAA and TTC/UUC sequences, thereby contributing to the field of trinucleotide repeats and the associated unusual structures that trigger expansion.

Molecular dynamics simulation on the Thermosinus carboxydivorans pfl ZTP riboswitch by ligand binding

Riboswitches are RNA molecules that can regulate gene expression which is affected by ligand-binding during cotranscriptional folding process. However, the role of ligand during the folding is still unclear. In this study, the pfl domain of Thermosinus carboxydivorans ZTP riboswitch was discussed. The ligand is molecule ZMP. We mainly analyzed the change of ZMP-free and ZMP-bound aptamer domain by the dynamics simulation method. Structural features by calculating their RMSD, RMSF, etc. are analyzed. The results demonstrate that the binding domain require the presence of ZMP to maintain a stable fold. It also suggested that ZMP specificly binding to ZTP can generate more hydrogen bonds in the binding domain. Through the calculation of binding free energy decomposition of each nucleotide, molecule ZMP was found to promote the recognition and binding process of ligands by controlling some special nucleotides in the process of ligand binding. At last, the dynamical correlation and components of conformational motions were both applied to explore the effect of molecule ZMP to ZTP riboswitch. In general, ZMP can effectively affect the motions of the pfl riboswitch and facilitate the folding process of the ZTP riboswitch.These results may provide some new ideas for structural changes in riboswitches and their cotranscriptional folding process.

Architectures and complex functions of tandem riboswitches

Riboswitch architectures that involve the binding of a single ligand to a single RNA aptamer domain result in ordinary dose-response curves that require approximately a 100-fold change in ligand concentration to cover nearly the full dynamic range for gene regulation. However, by using multiple riboswitches or aptamer domains in tandem, these ligand-sensing structures can produce additional, complex gene control outcomes. In the current study, we have computationally searched for tandem riboswitch architectures in bacteria to provide a more complete understanding of the diverse biological and biochemical functions of gene control elements that are made exclusively of RNA. Numerous different arrangements of tandem homologous riboswitch architectures are exploited by bacteria to create more 'digital' gene control devices, which operate over a narrower ligand concentration range. Also, two heterologous riboswitch aptamers are sometimes employed to create two-input Boolean logic gates with various types of genetic outputs. These findings illustrate the sophisticated genetic decisions that can be made by using molecular sensors and switches based only on RNA.

Incoherent dual regulation by a SAM-II riboswitch controlling translation at a distance

In Sinorhizobium meliloti, the methionine biosynthesis genes metA and metZ are preceded by S-adenosyl-L-methionine (SAM) riboswitches of the SAM-II class. Upon SAM binding, structural changes in the metZ riboswitch were predicted to cause transcriptional termination, generating the sRNA RZ. By contrast, the metA riboswitch was predicted to regulate translation from an AUG1 codon. However, downstream of the metA riboswitch, we found a putative Rho-independent terminator and an in-frame AUG2 codon, which may contribute to metA regulation. We validated the terminator between AUG1 and AUG2, which generates the sRNA RA1 that is processed to RA2. Under high SAM conditions, the activities of the metA and metZ promoters and the steady-state levels of the read-through metA and metZ mRNAs were decreased, while the levels of the RZ and RA2 sRNAs were increased. Under these conditions, the sRNAs and the mRNAs were stabilized. Reporter fusion experiments revealed that the Shine–Dalgarno (SD) sequence in the metA riboswitch is required for translation, which, however, starts 74 nucleotides downstream at AUG2, suggesting a novel translation initiation mechanism. Further, the reporter fusion data supported the following model of RNA-based regulation: Upon SAM binding by the riboswitch, the SD sequence is sequestered to downregulate metA translation, while the mRNA is stabilized. Thus, the SAM-II riboswitches fulfil incoherent, dual regulation, which probably serves to ensure basal metA and metZ mRNA levels under high SAM conditions. This probably helps to adapt to changing conditions and maintain SAM homoeostasis.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Charge Assisted S/Se Chalcogen Bonds in SAM Riboswitches: A Combined PDB and ab Initio Study

In this study, we provide experimental (Protein Data Bank (PDB) inspection) and theoretical (RI-MP2/def2-TZVP level of theory) evidence of the involvement of charge assisted chalcogen bonding (ChB) interactions in the recognition and folding mechanisms of S-adenosylmethionine (SAM) riboswitches. Concretely, an initial PDB search revealed several examples where ChBs between S-adenosyl methionine (SAM)/adenosyl selenomethionine (EEM) molecules and uracil (U) bases belonging to RNA take place. While these interactions are usually described as a merely Coulombic attraction between the positively charged S/Se group and RNA, theoretical calculations indicated that the σ holes of S and Se are involved. Moreover, computational models shed light on the strength and directionality properties of the interaction, which was also further characterized from a charge-density perspective using Bader’s “Atoms in Molecules” (AIM) theory, Non-Covalent Interaction plot (NCIplot) visual index, and Natural Bonding Orbital (NBO) analyses. As far as our knowledge extends, this is the first time that ChBs in SAM–RNA complexes have been systematically analyzed, and we believe the results might be useful for scientists working in the field of RNA engineering and chemical biology as well as to increase the visibility of the interaction among the biological community.

Riboswitch Mechanisms: New Tricks for an Old Dog

Discovered almost twenty years ago, riboswitches turned out to be one of the most common regulatory systems in bacteria, with representatives found in eukaryotes and archaea. Unlike many other regulatory elements, riboswitches are entirely composed of RNA and capable of modulating expression of genes by direct binding of small cellular molecules. While bacterial riboswitches had been initially thought to control production of enzymes and transporters associated with small organic molecules via feedback regulatory circuits, later findings identified riboswitches directing expression of a wide range of genes and responding to various classes of molecules, including ions, signaling molecules, and others. The 5'-untranslated mRNA regions host a vast majority of riboswitches, which modulate transcription or translation of downstream genes through conformational rearrangements in the ligand-sensing domains and adjacent expression-controlling platforms. Over years, the repertoire of regulatory mechanisms employed by riboswitches has greatly expanded; most recent studies have highlighted the importance of alternative mechanisms, such as RNA degradation, for the riboswitch-mediated genetic circuits. This review discusses the plethora of bacterial riboswitch mechanisms and illustrates how riboswitches utilize different features and approaches to elicit various regulatory responses.

Supersized Ribosomal RNA Expansion Segments in Asgard Archaea

The ribosome’s common core, comprised of ribosomal RNA (rRNA) and universal ribosomal proteins, connects all life back to a common ancestor and serves as a window to relationships among organisms. The rRNA of the common core is similar to rRNA of extant bacteria. In eukaryotes, the rRNA of the common core is decorated by expansion segments (ESs) that vastly increase its size. Supersized ESs have not been observed previously in Archaea, and the origin of eukaryotic ESs remains enigmatic. We discovered that the large ribosomal subunit (LSU) rRNA of two Asgard phyla, Lokiarchaeota and Heimdallarchaeota, considered to be the closest modern archaeal cell lineages to Eukarya, bridge the gap in size between prokaryotic and eukaryotic LSU rRNAs. The elongated LSU rRNAs in Lokiarchaeota and Heimdallarchaeota stem from two supersized ESs, called ES9 and ES39. We applied chemical footprinting experiments to study the structure of Lokiarchaeota ES39. Furthermore, we used covariation and sequence analysis to study the evolution of Asgard ES39s and ES9s. By defining the common eukaryotic ES39 signature fold, we found that Asgard ES39s have more and longer helices than eukaryotic ES39s. Although Asgard ES39s have sequences and structures distinct from eukaryotic ES39s, we found overall conservation of a three-way junction across the Asgard species that matches eukaryotic ES39 topology, a result consistent with the accretion model of ribosomal evolution.

The potential versatility of RNA catalysis

It is commonly thought that in the early development of life on this planet RNA would have acted both as a store of genetic information and as a catalyst. While a number of RNA enzymes are known in contemporary cells, they are largely confined to phosphoryl transfer reactions, whereas an RNA based metabolism would have required a much greater chemical diversity of catalysis. Here we discuss how RNA might catalyze a wider variety of chemistries, and particularly how information gleaned from riboswitches could suggest how ribozymes might recruit coenzymes to expand their chemical range. We ask how we might seek such activities in modern biology. This article is categorized under: RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions Regulatory RNAs/RNAi/Riboswitches > Riboswitches RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Transcription of Cystathionine β-Lyase (MetC) Is Repressed by HeuR in Campylobacter jejuni, and Methionine Biosynthesis Facilitates Colonocyte Invasion

A previously identified transcriptional regulator in Campylobacter jejuni, termed HeuR, was found to positively regulate heme utilization. Additionally, transcriptomic work demonstrated that the putative operons CJJ81176_1390 to CJJ81176_1394 (CJJ81176_1390-1394) and CJJ81176_1214-1217 were upregulated in a HeuR mutant, suggesting that HeuR negatively regulates expression of these genes. Because genes within these clusters include a cystathionine β-lyase (metC) and a methionine synthase (metE), it appeared HeuR negatively regulates C. jejuni methionine biosynthesis. To address this, we confirmed mutation of HeuR reproducibly results in metC overexpression under nutrient-replete conditions but did not affect expression of metE, while metC expression in the wild type increased to heuR mutant levels during iron limitation. We subsequently determined that both gene clusters are operonic and demonstrated the direct interaction of HeuR with the predicted promoter regions of these operons. Using DNase footprinting assays, we were able to show that HeuR specifically binds within the predicted -35 region of the CJJ81176_1390-1394 operon. As predicted based on transcriptional results, the HeuR mutant was able to grow and remain viable in a defined medium with and without methionine, but we identified significant impacts on growth and viability in metC and metE mutants. Additionally, we observed decreased adherence, invasion, and persistence of metC and metE mutants when incubated with human colonocytes, while the heuR mutant exhibited increased invasion. Taken together, these results suggest that HeuR regulates methionine biosynthesis in an iron-responsive manner and that the ability to produce methionine is an important factor for adhering to and invading the gastrointestinal tract of a susceptible host. IMPORTANCE As the leading cause of bacterium-derived gastroenteritis worldwide, Campylobacter jejuni has a significant impact on human health. Investigating colonization factors that allow C. jejuni to successfully infect a host furthers our understanding of genes and regulatory elements necessary for virulence. In this study, we have begun to characterize the role of the transcriptional regulatory protein, HeuR, on methionine biosynthesis in C. jejuni. When the ability to synthesize methionine is impaired, detrimental impacts on growth and viability are observed during growth in limited media lacking methionine and/or iron. Additionally, mutations in the methionine biosynthetic pathway result in decreased adhesion, invasion, and intracellular survival of C. jejuni when incubated with human colonocytes, indicating the importance of regulating methionine biosynthesis.

Fundamental photophysics of isomorphic and expanded fluorescent nucleoside analogues

Fluorescent nucleoside analogues (FNAs) are structurally diverse mimics of the natural essentially non-fluorescent nucleosides which have found numerous applications in probing the structure and dynamics of nucleic acids as well as their interactions with various biomolecules. In order to minimize disturbance in the labelled nucleic acid sequences, the FNA chromophoric groups should resemble the natural nucleobases in size and hydrogen-bonding patterns. Isomorphic and expanded FNAs are the two groups that best meet the criteria of non-perturbing fluorescent labels for DNA and RNA. Significant progress has been made over the past decades in understanding the fundamental photophysics that governs the spectroscopic and environmentally sensitive properties of these FNAs. Herein, we review recent advances in the spectroscopic and computational studies of selected isomorphic and expanded FNAs. We also show how this information can be used as a rational basis to design new FNAs, select appropriate sequences for optimal spectroscopic response and interpret fluorescence data in FNA applications.

Recent Insights into the Structure and Function of Mycobacterial Membrane Proteins Facilitated by Cryo-EM

Mycobacterium tuberculosis (Mtb) is one of the deadliest pathogens encountered by humanity. Over the decades, its characteristic membrane organization and composition have been understood. However, there is still limited structural information and mechanistic understanding of the constituent membrane proteins critical for drug discovery pipelines. Recent advances in single-particle cryo-electron microscopy and cryo-electron tomography have provided the much-needed impetus towards structure determination of several vital Mtb membrane proteins whose structures were inaccessible via X-ray crystallography and NMR. Important insights into membrane composition and organization have been gained via a combination of electron tomography and biochemical and biophysical assays. In addition, till the time of writing this review, 75 new structures of various Mtb proteins have been reported via single-particle cryo-EM. The information obtained from these structures has improved our understanding of the mechanisms of action of these proteins and the physiological pathways they are associated with. These structures have opened avenues for structure-based drug design and vaccine discovery programs that might help achieve global-TB control. This review describes the structural features of selected membrane proteins (type VII secretion systems, Rv1819c, Arabinosyltransferase, Fatty Acid Synthase, F-type ATP synthase, respiratory supercomplex, ClpP1P2 protease, ClpB disaggregase and SAM riboswitch), their involvement in physiological pathways, and possible use as a drug target. Tuberculosis is a deadly disease caused by Mycobacterium tuberculosis. The Cryo-EM and tomography have simplified the understanding of the mycobacterial membrane organization. Some proteins are located in the plasma membrane; some span the entire envelope, while some, like MspA, are located in the mycomembrane. Cryo-EM has made the study of such membrane proteins feasible.

RNA stabilization by a poly(A) tail 3'-end binding pocket and other modes of poly(A)-RNA interaction

Polyadenylate [poly(A)] tail addition to the 3' end of a wide range of RNAs is a highly conserved modification that plays a central role in cellular RNA function. Elements for nuclear expression (ENEs) are cis-acting RNA elements that stabilize poly(A) tails by sequestering them in RNA triplex structures. A crystal structure of a double ENE from a rice hAT transposon messenger RNA complexed with poly(A)₂₈ at a resolution of 2.89 angstroms reveals multiple modes of interaction with poly(A), including major-groove triple helices, extended minor-groove interactions with RNA double helices, a quintuple-base motif that transitions poly(A) from minor-groove associations to major-groove triple helices, and a poly(A) 3'-end binding pocket. Our findings both expand the repertoire of motifs involved in long-range RNA interactions and provide insights into how polyadenylation can protect an RNA's extreme 3' end.

Tertiary Base Triple Formation in the SRV-1 Frameshifting Pseudoknot Stabilizes Secondary Structure Components

Minor-groove base triples formed between stem 1 and loop 2 of the simian retrovirus type 1 (SRV-1) mRNA frameshifting pseudoknot are essential in stimulating -1 ribosomal frameshifting. How tertiary base triple formation affects the local stabilities of secondary structures (stem 1 and stem 2) and thus ribosomal frameshifting efficiency is not well understood. We made a short peptide nucleic acid (PNA) that is expected to invade stem 1 of the SRV-1 pseudoknot by PNA-RNA duplex formation to mimic the stem 1 unwinding process by a translating ribosome. In addition, we used a PNA for invading stem 2 in the SRV-1 pseudoknot. Our nondenaturing polyacrylamide gel electrophoresis data for the binding of PNA to the SRV-1 pseudoknot and mutants reveal that mutations in loop 2 disrupting base triple formation between loop 2 and stem 1 in the SRV-1 pseudoknot result in enhanced invasion by both PNAs. Our data suggest that tertiary stem 1-loop 2 base triple interactions in the SRV-1 pseudoknot can stabilize both of the secondary structural components, stem 1 and stem 2. Stem 2 stability is thus coupled to the structural stability of stem 1-loop 2 base triples, mediated through a long-range effect. The apparent dissociation constants of both PNAs are positively correlated with the pseudoknot mechanical stabilities and frameshifting efficiencies. The relatively simple PNA local invasion experiment may be used to characterize the energetic contribution of tertiary interactions and ligand binding in many other RNA and DNA structures.

Substitution-Inert Polynuclear Platinum Complexes Inhibit Reverse Transcription Preferentially in RNA Triplex-Forming Templates

RNA triplexes are significant tertiary structure motifs that are found in many functional RNAs. Hence, small molecules capable of recognition, binding, and stabilization of the triple-helical RNA structures are emerging as attractive potential molecular biology tools and therapeutic agents. Here, we utilize methods of molecular biology and biophysics to study the interactions of a series of antitumor substitution-inert polynuclear platinum complexes (SI-PPCs) with triple-helical RNA structures. We show that SI-PPCs recognize and stabilize RNA triplexes and inhibit reverse transcription preferentially in the RNA template prone to the triplex formation. These so far unexplored properties of SI-PPCs suggest that the targeting of triple-stranded regions in RNA might contribute to the biological effects of SI-PPCs.

Imaginary Ribozymes

Biocatalysis is dominated by protein enzymes, and only a few classes of ribozymes are known to contribute to the task of promoting biochemical transformations. The RNA World theory encompasses the notion that earlier forms of life made use of a much greater diversity of ribozymes and other functional RNAs to guide complex metabolic states long before proteins had emerged in evolution. In recent years, the discoveries of various classes of ribozymes, riboswitches, and other noncoding RNAs in bacteria have provided additional support for the hypothesis that RNA molecules indeed have the catalytic competence to promote diverse chemical reactions without the aid of protein enzymes. Herein, some of the most striking observations made from examinations of natural riboswitches that bind small ligands are highlighted and used as a basis to imagine the characteristics and functions of long-extinct ribozymes from the RNA World.

Crystal structure and ligand-induced folding of the SAM/SAH riboswitch

While most SAM riboswitches strongly discriminate between SAM and SAH, the SAM/SAH riboswitch responds to both ligands with similar apparent affinities. We have determined crystal structures of the SAM/SAH riboswitch bound to SAH, SAM and other variant ligands at high resolution. The riboswitch forms an H-type pseudoknot structure with coaxial alignment of the stem–loop helix (P1) and the pseudoknot helix (PK). An additional three base pairs form at the non-open end of P1, and the ligand is bound at the interface between the P1 extension and the PK helix. The adenine nucleobase is stacked into the helix and forms a trans Hoogsteen–Watson–Crick base pair with a uridine, thus becoming an integral part of the helical structure. The majority of the specific interactions are formed with the adenosine. The methionine or homocysteine chain lies in the groove making a single hydrogen bond, and there is no discrimination between the sulfonium of SAM or the thioether of SAH. Single-molecule FRET analysis reveals that the riboswitch exists in two distinct conformations, and that addition of SAM or SAH shifts the population into a stable state that likely corresponds to the form observed in the crystal. A model for translational regulation is presented whereby in the absence of ligand the riboswitch is largely unfolded, lacking the PK helix so that translation can be initiated at the ribosome binding site. But the presence of ligand stabilizes the folded conformation that includes the PK helix, so occluding the ribosome binding site and thus preventing the initiation of translation.

A novel riboswitch classification based on imbalanced sequences achieved by machine learning

Riboswitch, a part of regulatory mRNA (50-250nt in length), has two main classes: aptamer and expression platform. One of the main challenges raised during the classification of riboswitch is imbalanced data. That is a circumstance in which the records of a sequences of one group are very small compared to the others. Such circumstances lead classifier to ignore minority group and emphasize on majority ones, which results in a skewed classification. We considered sixteen riboswitch families, to be in accord with recent riboswitch classification work, that contain imbalanced sequences. The sequences were split into training and test set using a newly developed pipeline. From 5460 k-mers (k value 1 to 6) produced, 156 features were calculated based on CfsSubsetEval and BestFirst function found in WEKA 3.8. Statistically tested result was significantly difference between balanced and imbalanced sequences (p < 0.05). Besides, each algorithm also showed a significant difference in sensitivity, specificity, accuracy, and macro F-score when used in both groups (p < 0.05). Several k-mers clustered from heat map were discovered to have biological functions and motifs at the different positions like interior loops, terminal loops and helices. They were validated to have a biological function and some are riboswitch motifs. The analysis has discovered the importance of solving the challenges of majority bias analysis and overfitting. Presented results were generalized evaluation of both balanced and imbalanced models, which implies their ability of classifying, to classify novel riboswitches. The Python source code is available at https://github.com/Seasonsling/riboswitch.

Molecular structure of a U•A-U-rich RNA triple helix with 11 consecutive base triples

Three-dimensional structures have been solved for several naturally occurring RNA triple helices, although all are limited to six or fewer consecutive base triples, hindering accurate estimation of global and local structural parameters. We present an X-ray crystal structure of a right-handed, U•A-U-rich RNA triple helix with 11 continuous base triples. Due to helical unwinding, the RNA triple helix spans an average of 12 base triples per turn. The double helix portion of the RNA triple helix is more similar to both the helical and base step structural parameters of A′-RNA rather than A-RNA. Its most striking features are its wide and deep major groove, a smaller inclination angle and all three strands favoring a C3′-endo sugar pucker. Despite the presence of a third strand, the diameter of an RNA triple helix remains nearly identical to those of DNA and RNA double helices. Contrary to our previous modeling predictions, this structure demonstrates that an RNA triple helix is not limited in length to six consecutive base triples and that longer RNA triple helices may exist in nature. Our structure provides a starting point to establish structural parameters of the so-called ‘ideal’ RNA triple helix, analogous to A-RNA and B-DNA double helices.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

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

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

Computational study of switching mechanism in add A-riboswitch

Riboswitch can bind small molecules to regulate gene expression. Unlike other RNAs, riboswitch relies on its conformational switching for regulation. However, the understanding of the switching mechanism is still limited. Here, we focussed on the add A-riboswitch to illustrate the dynamical switching mechanism as an example. We performed molecular dynamics simulation, conservation and co-evolution calculations to infer the dynamical motions and evolutionary base pairings. The results suggest that the binding domain is stable for molecule recognition and binding, whereas the switching base pairings are co-evolutionary for translation. The understanding of the add A-riboswitch switching mechanism provides a potential solution for riboswitch drug design.

Interactions between SAM and the 5' UTR mRNA of the sam1 gene regulate translation in S. pombe

The 5' untranslated region (5' UTR) of eukaryotic mRNA plays an important role in translation. Here we report the function of the 5' UTR mRNA of S-adenosylmethionine synthetase (sam1) in translational modulation in the presence of SAM in fission yeast Schizosaccharomyces pombe Reporter assays, binding and chemical probing experiments, and mutational analysis show that the 5' UTR mRNA of sam1 binds to SAM to effect translation. Translational modulation is dependent on a tertiary structure transition in the RNA upon SAM binding. The characterization of such an RNA that is directly associated with an essential metabolic process in eukaryotes provides additional evidence that ligand binding by RNAs plays an important role in eukaryotic gene regulation.

Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations

Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.

SAM-VI riboswitch structure and signature for ligand discrimination

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

Structure and ligand binding of the glutamine-II riboswitch

We have determined the structure of the glutamine-II riboswitch ligand binding domain using X-ray crystallography. The structure was solved using a novel combination of homology modeling and molecular replacement. The structure comprises three coaxial helical domains, the central one of which is a pseudoknot with partial triplex character. The major groove of this helix provides the binding site for L-glutamine, which is extensively hydrogen bonded to the RNA. Atomic mutation of the RNA at the ligand binding site leads to loss of binding shown by isothermal titration calorimetry, explaining the specificity of the riboswitch. A metal ion also plays an important role in ligand binding. This is directly bonded to a glutamine carboxylate oxygen atom, and its remaining inner-sphere water molecules make hydrogen bonding interactions with the RNA.

Cryo-EM structure of a 40 kDa SAM-IV riboswitch RNA at 3.7 Å resolution

Specimens below 50 kDa have generally been considered too small to be analyzed by single-particle cryo-electron microscopy (cryo-EM). The high flexibility of pure RNAs makes it difficult to obtain high-resolution structures by cryo-EM. In bacteria, riboswitches regulate sulfur metabolism through binding to the S-adenosylmethionine (SAM) ligand and offer compelling targets for new antibiotics. SAM-I, SAM-I/IV, and SAM-IV are the three most commonly found SAM riboswitches, but the structure of SAM-IV is still unknown. Here, we report the structures of apo and SAM-bound SAM-IV riboswitches (119-nt, ~40 kDa) to 3.7 Å and 4.1 Å resolution, respectively, using cryo-EM. The structures illustrate homologies in the ligand-binding core but distinct peripheral tertiary contacts in SAM-IV compared to SAM-I and SAM-I/IV. Our results demonstrate the feasibility of resolving small RNAs with enough detail to enable detection of their ligand-binding pockets and suggest that cryo-EM could play a role in structure-assisted drug design for RNA.

High-affinity recognition of specific tRNAs by an mRNA anticodon-binding groove

T-box riboswitches are modular bacterial noncoding RNAs that sense and regulate amino acid availability through direct interactions with tRNAs. Between the 5' anticodon-binding stem I domain and the 3' amino acid sensing domains of most T-boxes lies the stem II domain of unknown structure and function. Here, we report a 2.8-Å cocrystal structure of the Nocardia farcinica ileS T-box in complex with its cognate tRNA^Ile. The structure reveals a perpendicularly arranged ultrashort stem I containing a K-turn and an elongated stem II bearing an S-turn. Both stems rest against a compact pseudoknot, dock via an extended ribose zipper and jointly create a binding groove specific to the anticodon of its cognate tRNA. Contrary to proposed distal contacts to the tRNA elbow region, stem II locally reinforces the codon-anticodon interactions between stem I and tRNA, achieving low-nanomolar affinity. This study illustrates how mRNA junctions can create specific binding sites for interacting RNAs of prescribed sequence and structure.

Stability of an RNA•DNA-DNA triple helix depends on base triplet composition and length of the RNA third strand

Recent studies suggest noncoding RNAs interact with genomic DNA, forming an RNA•DNA–DNA triple helix that regulates gene expression. However, base triplet composition of pyrimidine motif RNA•DNA–DNA triple helices is not well understood beyond the canonical U•A–T and C•G–C base triplets. Using native gel-shift assays, the relative stability of 16 different base triplets at a single position, Z•X–Y (where Z = C, U, A, G and X–Y = A–T, G–C, T–A, C–G), in an RNA•DNA–DNA triple helix was determined. The canonical U•A–T and C•G–C base triplets were the most stable, while three non-canonical base triplets completely disrupted triple-helix formation. We further show that our RNA•DNA–DNA triple helix can tolerate up to two consecutive non-canonical A•G–C base triplets. Additionally, the RNA third strand must be at least 19 nucleotides to form an RNA•DNA–DNA triple helix but increasing the length to 27 nucleotides does not increase stability. The relative stability of 16 different base triplets in DNA•DNA–DNA and RNA•RNA–RNA triple helices was distinctly different from those in RNA•DNA–DNA triple helices, showing that base triplet stability depends on strand composition being DNA and/or RNA. Multiple factors influence the stability of triple helices, emphasizing the importance of experimentally validating formation of computationally predicted triple helices.

The structure of the SAM/SAH-binding riboswitch

S-adenosylmethionine (SAM) is a central metabolite since it is used as a methyl group donor in many different biochemical reactions. Many bacteria control intracellular SAM concentrations using riboswitch-based mechanisms. A number of structurally different riboswitch families specifically bind to SAM and mainly regulate the transcription or the translation of SAM-biosynthetic enzymes. In addition, a highly specific riboswitch class recognizes S-adenosylhomocysteine (SAH)—the product of SAM-dependent methyl group transfer reactions—and regulates enzymes responsible for SAH hydrolysis. High-resolution structures are available for many of these riboswitch classes and illustrate how they discriminate between the two structurally similar ligands SAM and SAH. The so-called SAM/SAH riboswitch class binds both ligands with similar affinities and is structurally not yet characterized. Here, we present a high-resolution nuclear magnetic resonance structure of a member of the SAM/SAH-riboswitch class in complex with SAH. Ligand binding induces pseudoknot formation and sequestration of the ribosome binding site. Thus, the SAM/SAH-riboswitches are translational ‘OFF’-switches. Our results establish a structural basis for the unusual bispecificity of this riboswitch class. In conjunction with genomic data our structure suggests that the SAM/SAH-riboswitches might be an evolutionary late invention and not a remnant of a primordial RNA-world as suggested for other riboswitches.

Finely tuned conformational dynamics regulate the protective function of the lncRNA MALAT1 triple helix

Nucleic acid triplexes may regulate many important biological processes. Persistent accumulation of the oncogenic 7-kb long noncoding RNA MALAT1 is dependent on an unusually long intramolecular triple helix. This triplex structure is positioned within a conserved ENE (element for nuclear expression) motif at the lncRNA 3′ terminus and protects the entire transcript from degradation in a polyA-independent manner. A requisite 3′ maturation step leads to triplex formation though the precise mechanism of triplex folding remains unclear. Furthermore, the contributions of several peripheral structural elements to triplex formation and protective function have not been determined. We evaluated the stability, conformational fluctuations, and function of this MALAT1 ENE triple helix (M1^TH) protective element using in vitro mutational analyses coupled with biochemical and biophysical characterizations. Using fluorescence and UV melts, FRET, and an exonucleolytic decay assay we define a concerted mechanism for triplex formation and uncover a metastable, dynamic triplex population under near-physiological conditions. Structural elements surrounding the triplex regulate the dynamic M1^TH conformational variability, but increased triplex dynamics lead to M1^TH degradation. Taken together, we suggest that finely tuned dynamics may be a general mechanism regulating triplex-mediated functions.