Structure and function of preQ1 riboswitches

A (Scientific) Lifetime Affair With Nucleic Acids

I am Distinguished Professor in the Chemistry and Biochemistry Department at University of California, Los Angeles, where I was hired in 1985 as the first female assistant professor in the department. I received my PhD from University of California, San Diego, under the guidance of Professor David Kearns, where I used NMR spectroscopy to study drug binding to random sequence DNA and published the first two-dimensional NMR spectra of short synthetic DNA duplexes. From 1982 to 1985 I was a Damon Runyon-Walter Winchell Postdoctoral fellow in the Professor Alexander Rich laboratory, where I investigated structures of Z-DNA by NMR. At UCLA, my lab pioneered the application of macromolecular NMR spectroscopy to the study of DNA and RNA structure, folding, and interactions with cations, drugs, and proteins. We published the first NMR structures of DNA triplexes, quadruplexes, and aptamers, and our work has provided fundamental insights into DNA A-tract bending, cation interactions with DNA, Hoogsteen base pairs, and drug binding to DNA. My lab has made major contributions to understanding RNA folding, dynamics, and function, including pseudoknots, aptamers, ribozymes, and riboswitches, and recognition of RNA by proteins. Over the past 2 decades, the Feigon laboratory pioneered structure-function studies of telomerase, from solution NMR and X-ray crystal structures and dynamics studies of RNA and RNA-protein domains of human and Tetrahymena telomerase, to the first structure of a telomerase holoenzyme, by negative stain EM in 2013, and subsequent cryo-EM structures of telomerase and associated proteins. Recent work also includes structural biology of 7SK RNP.

Mechanistic analysis of Riboswitch Ligand interactions provides insights into pharmacological control over gene expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

MicroRNA (miRNA) as a biomarker for diagnosis, prognosis, and therapeutics molecules in neurodegenerative disease

Neurodegenerative diseases that include Alzheimer's disease (AD), amyotrophic lateral sclerosis (ALS), Parkinson's disease (PD), Huntington's disease (HD), and multiple sclerosis (MS) that arise due to numerous causes like protein accumulation and autoimmunity characterized by neurologic depletion which lead to incapacity in normal physiological function such as thinking and movement in these patients. Glial cells perform an important role in protective neuronal function; in the case of neuroinflammation, glial cell dysfunction can promote the development of neurodegenerative diseases. miRNA that participates in gene regulation and plays a vital role in many biological processes in the body; in the central nervous system (CNS), it can play an essential part in neural maturation and differentiation. In neurodegenerative diseases, miRNA dysregulation occurs, enhancing the development of these diseases. In this review, we discuss neurodegenerative disease (Alzheimer's disease (AD), Parkinson's disease (PD), amyotrophic lateral sclerosis (ALS), and multiple sclerosis (MS)) and how miRNA is preserved as a diagnostic biomarker or therapeutic agent in these disorders. Finally, we highlight miRNA as therapy.

Mechanistic Analysis of Riboswitch Ligand Interactions Provides Insights into Pharmacological Control over Gene Expression

Riboswitches are structured RNA elements that regulate gene expression upon binding to small molecule ligands. Understanding the mechanisms by which small molecules impact riboswitch activity is key to developing potent, selective ligands for these and other RNA targets. We report the structure-informed design of chemically diverse synthetic ligands for PreQ1 riboswitches. Multiple X-ray co-crystal structures of synthetic ligands with the Thermoanaerobacter tengcongensis (Tte)-PreQ1 riboswitch confirm a common binding site with the cognate ligand, despite considerable chemical differences among the ligands. Structure probing assays demonstrate that one ligand causes conformational changes similar to PreQ1 in six structurally and mechanistically diverse PreQ1 riboswitch aptamers. Single-molecule force spectroscopy is used to demonstrate differential modes of riboswitch stabilization by the ligands. Binding of the natural ligand brings about the formation of a persistent, folded pseudoknot structure, whereas a synthetic ligand decreases the rate of unfolding through a kinetic mechanism. Single round transcription termination assays show the biochemical activity of the ligands, while a GFP reporter system reveals compound activity in regulating gene expression in live cells without toxicity. Taken together, this study reveals that diverse small molecules can impact gene expression in live cells by altering conformational changes in RNA structures through distinct mechanisms.

Multidimensional Applications and Challenges of Riboswitches in Biosensing and Biotherapy

Riboswitches have received significant attention over the last two decades for their multiple functionalities and great potential for applications in various fields. This article highlights and reviews the recent advances in biosensing and biotherapy. These fields involve a wide range of applications, such as food safety detection, environmental monitoring, metabolic engineering, live cell imaging, wearable biosensors, antibacterial drug targets, and gene therapy. The discovery, origin, and optimization of riboswitches are summarized to help readers better understand their multidimensional applications. Finally, this review discusses the multidimensional challenges and development of riboswitches in order to further expand their potential for novel applications.

Thermal titration molecular dynamics (TTMD): shedding light on the stability of RNA-small molecule complexes

Ribonucleic acids are gradually becoming relevant players among putative drug targets, thanks to the increasing amount of structural data exploitable for the rational design of selective and potent binders that can modulate their activity. Mainly, this information allows employing different computational techniques for predicting how well would a ribonucleic-targeting agent fit within the active site of its target macromolecule. Due to some intrinsic peculiarities of complexes involving nucleic acids, such as structural plasticity, surface charge distribution, and solvent-mediated interactions, the application of routinely adopted methodologies like molecular docking is challenged by scoring inaccuracies, while more physically rigorous methods such as molecular dynamics require long simulation times which hamper their conformational sampling capabilities. In the present work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of unbinding kinetics, to characterize RNA-ligand complexes. In this article, we explored its applicability as a post-docking refinement tool on RNA in complex with small molecules, highlighting the capability of this method to identify the native binding mode among a set of decoys across various pharmaceutically relevant test cases.

Dynamical characterization and multiple unbinding paths of two PreQ1 ligands in one pocket

Riboswitches naturally regulate gene expression in bacteria by binding to specific small molecules. Class 1 preQ1 riboswitch aptamer is an important model not only for RNA folding but also as a target for designing small molecule antibiotics due to its well-known minimal aptamer domain. Here, we ran a total of 62.4 μs conventional and enhanced-sampling molecular dynamics (MD) simulations to characterize the determinants underlying the binding of the preQ1-II riboswitch aptamer to two preQ1 ligands in one binding pocket. Decomposition of binding free energy suggested that preQ1 ligands at α and β sites interact with four nucleotides (G5, C17, C18, and A30) and two nucleotides (A12 and C31), respectively. Mg2+ ions play a crucial role in both stabilizing the binding pocket and facilitating ligand binding. The flexible preQ1 ligand at the β site leads to the top of the binding pocket loosening and thus pre-organizes the riboswitch for ligand entry. Enhanced sampling simulations further revealed that the preQ1 ligand at the α site unbinds through two orthogonal pathways, which are dependent on whether or not a β site preQ1 ligand is present. One of the two preQ1 ligands has been identified in the binding pocket, which will aid to identify the second preQ1 Ligand. Our work provides new information for designing robust ligands.

Discovering riboswitches: the past and the future

Riboswitches are structured noncoding RNA domains used by many bacteria to monitor the concentrations of target ligands and regulate gene expression accordingly. In the past 20 years over 55 distinct classes of natural riboswitches have been discovered that selectively sense small molecules or elemental ions, and thousands more are predicted to exist. Evidence suggests that some riboswitches might be direct descendants of the RNA-based sensors and switches that were likely present in ancient organisms before the evolutionary emergence of proteins. We provide an overview of the current state of riboswitch research, focusing primarily on the discovery of riboswitches, and speculate on the major challenges facing researchers in the field.

Probing Transient Riboswitch Structures via Single Molecule Accessibility Analysis

Riboswitches are a class of RNA motifs in the untranslated regions of bacterial messenger RNAs (mRNAs) that can adopt different conformations to regulate gene expression. The binding of specific small molecule or ion ligands, or other RNAs, influences the conformation the riboswitch adopts. Single Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) offers an approach for probing this structural isomerization, or conformational switching, at the level of single mRNA molecules. SiM-KARTS utilizes fluorescently labeled, short, sequence-complementary DNA or RNA oligonucleotide probes that transiently access a specific RNA conformation over another. Binding and dissociation to a surface-immobilized target RNA of arbitrary length are monitored by Total Internal Reflection Fluorescence Microscopy (TIRFM) and quantitatively analyzed, via spike train and burst detection, to elucidate the rate constants of isomerization, revealing mechanistic insights into riboswitching.

Synthetic biology tools to promote the folding and function of RNA aptamers in mammalian cells

RNA aptamers are structured RNAs that can bind to diverse ligands, including proteins, metabolites, and other small molecules. RNA aptamers are widely used as in vitro affinity reagents. However, RNA aptamers have not been highly successful as bioactive intracellular molecules that can bind target molecules and influence cellular processes. We describe how poor RNA aptamer expression and especially poor RNA aptamer folding have limited the use of RNA aptamers in RNA synthetic biology applications. We discuss innovative new approaches that promote RNA aptamer folding in living cells and how these approaches have improved the function of aptamers in mammalian cells. These new approaches are making RNA aptamer-based synthetic biology and RNA aptamer therapeutic applications much more achievable.

Bioinformatics and Genomic Analyses of the Suitability of Eight Riboswitches for Antibacterial Drug Targets

Antibiotic resistance (AR) is an acute problem that results in prolonged and debilitating illnesses. AR mortality worldwide is growing and causes a pressing need to research novel mechanisms of action and untested target molecules. This article presents in silico analyses of eight bacterial riboswitches for their suitability for antibacterial drug targets. Most bacterial riboswitches are located in the 5′-untranslated region of messenger RNAs, act as allosteric cis-acting gene control elements, and have not been found in humans before. Sensing metabolites, the riboswitches regulate the synthesis of vital cellular metabolites in various pathogenic bacteria. The analyses performed in this article represent a complete and informative genome-wide bioinformatics analysis of the adequacy of eight riboswitches as antibacterial drug targets in different pathogenic bacteria based on four criteria. Due to the ability of the riboswitch to control biosynthetic pathways and transport proteins of essential metabolites and the presence/absence of alternative biosynthetic pathways, we classified them into four groups based on their suitability for use as antibacterial drug targets guided by our in silico analyses. We concluded that some of them are promising targets for antibacterial drug discovery, such as the PreQ1, MoCo RNA, cyclic-di-GMP I, and cyclic-di-GMP II riboswitches.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

Precise tuning of bacterial translation initiation by non-equilibrium 5'-UTR unfolding observed in single mRNAs

Noncoding, structured 5′-untranslated regions (5′-UTRs) of bacterial messenger RNAs (mRNAs) can control translation efficiency by forming structures that either recruit or repel the ribosome. Here we exploit a 5′-UTR embedded preQ₁-sensing, pseudoknotted translational riboswitch to probe how binding of a small ligand controls recruitment of the bacterial ribosome to the partially overlapping Shine-Dalgarno (SD) sequence. Combining single-molecule fluorescence microscopy with mutational analyses, we find that the stability of 30S ribosomal subunit binding is inversely correlated with the free energy needed to unfold the 5′-UTR during mRNA accommodation into the mRNA binding cleft. Ligand binding to the riboswitch stabilizes the structure to both antagonize 30S recruitment and accelerate 30S dissociation. Proximity of the 5′-UTR and stability of the SD:anti-SD interaction both play important roles in modulating the initial 30S-mRNA interaction. Finally, depletion of small ribosomal subunit protein S1, known to help resolve structured 5′-UTRs, further increases the energetic penalty for mRNA accommodation. The resulting model of rapid standby site exploration followed by gated non-equilibrium unfolding of the 5′-UTR during accommodation provides a mechanistic understanding of how translation efficiency is governed by riboswitches and other dynamic structure motifs embedded upstream of the translation initiation site of bacterial mRNAs.

A small RNA that cooperatively senses two stacked metabolites in one pocket for gene control

Riboswitches are structured non-coding RNAs often located upstream of essential genes in bacterial messenger RNAs. Such RNAs regulate expression of downstream genes by recognizing a specific cellular effector. Although nearly 50 riboswitch classes are known, only a handful recognize multiple effectors. Here, we report the 2.60-Å resolution co-crystal structure of a class I type I preQ₁-sensing riboswitch that reveals two effectors stacked atop one another in a single binding pocket. These effectors bind with positive cooperativity in vitro and both molecules are necessary for gene regulation in bacterial cells. Stacked effector recognition appears to be a hallmark of the largest subgroup of preQ₁ riboswitches, including those from pathogens such as Neisseria gonorrhoeae. We postulate that binding to stacked effectors arose in the RNA World to closely position two substrates for RNA-mediated catalysis. These findings expand known effector recognition capabilities of riboswitches and have implications for antimicrobial development.

Riboswitches contain an aptamer domain that recognizes a metabolite and an expression platform that regulates gene expression. Here the authors report the crystal structure of a preQ₁-sensing riboswitch from Carnobacterium antarcticus that shows two metabolites in a single binding pocket.

A chemical probe based on the PreQ1 metabolite enables transcriptome-wide mapping of binding sites

The role of metabolite-responsive riboswitches in regulating gene expression in bacteria is well known and makes them useful systems for the study of RNA-small molecule interactions. Here, we study the PreQ₁ riboswitch system, assessing sixteen diverse PreQ₁-derived probes for their ability to selectively modify the class-I PreQ₁ riboswitch aptamer covalently. For the most active probe (11), a diazirine-based photocrosslinking analog of PreQ₁, X-ray crystallography and gel-based competition assays demonstrated the mode of binding of the ligand to the aptamer, and functional assays demonstrated that the probe retains activity against the full riboswitch. Transcriptome-wide mapping using Chem-CLIP revealed a highly selective interaction between the bacterial aptamer and the probe. In addition, a small number of RNA targets in endogenous human transcripts were found to bind specifically to 11, providing evidence for candidate PreQ₁ aptamers in human RNA. This work demonstrates a stark influence of linker chemistry and structure on the ability of molecules to crosslink RNA, reveals that the PreQ₁ aptamer/ligand pair are broadly useful for chemical biology applications, and provides insights into how PreQ₁, which is similar in structure to guanine, interacts with human RNAs.

Siblings or doppelgängers? Deciphering the evolution of structured cis-regulatory RNAs beyond homology

Structured cis-regulatory RNAs have evolved across all domains of life, highlighting the utility and plasticity of RNA as a regulatory molecule. Homologous RNA sequences and structures often have similar functions, but homology may also be deceiving. The challenges that derive from trying to assign function to structure and vice versa are not trivial. Bacterial riboswitches, viral and eukaryotic IRESes, CITEs, and 3′ UTR elements employ an array of mechanisms to exert their effects. Bioinformatic searches coupled with biochemical and functional validation have elucidated some shared and many unique ways cis-regulators are employed in mRNA transcripts. As cis-regulatory RNAs are resolved in greater detail, it is increasingly apparent that shared homology can mask the full spectrum of mRNA cis-regulator functional diversity. Furthermore, similar functions may be obscured by lack of obvious sequence similarity. Thus looking beyond homology is crucial for furthering our understanding of RNA-based regulation.

Comparison of clinicopathological parameters, prognosis, micro-ecological environment and metabolic function of Gastric Cancer with or without Fusobacterium sp. Infection

Background: Fusobacterium sp. plays a crucial role in the tumorigenesis and development of gastrointestinal tumors. Our research group previously disclosed that Fusobacterium sp. was more abundant in gastric cancer (GC) tissues than adjacent non-cancerous (NC) tissues. However, Fusobacterium sp. did not exist in all GC tissues and the differentiated features of GC with or without Fusobacterium sp. infection is not clear.

Methods: The expression data of 61 GC tissues came from 16S rRNA gene sequencing. Comparison groups were defined based on sOTU at the genus level of Fusobacterium sp., which was performed by the Qiime2 microbiome bioinformatics platform. We used Chi-square and Fisher's exact test to compare clinicopathological parameters, and used Kaplan-Meier analysis, Cox univariate and multivariate analysis to compare prognosis. Micro-ecological environment comparison was characterized by 16S rRNA gene sequencing, and the metabolic function prediction was applied by PICRUSt2. Results of microbial diversity, differential enrichment genus and metabolic function in GC with or without Fusobacterium sp. infection was validated with 229 GC tissues downloaded from an independent cohort in ENA database (PRJNA428883).

Results: The infection rate of Fusobacterium sp. in 61 GC tissues was 52.46% and elderly GC patients were more prone to Fusobacterium sp. infection. GC patients infected with Fusobacterium sp. were more likely to have tumor-infiltrating lymphocytes and p53 expression. The microbial diversity and microbial structure showed significant differences between two GC tissue groups with 42 differential enrichment genera. The metabolic function of Fusobacterium sp.-positive GC tissues was related to the biosynthesis of lysine, peptidoglycan, and tRNA. The differences in microbial structure, the existence of some differential enrichment genera and the metabolic function of Fusobacterium sp.-positive GC tissues, were then validated by 229 GC tissues of an independent cohort.

Conclusions: Fusobacterium sp. infection can affect the phenotypic characteristics, micro-ecological environment, and metabolic functions of GC, which may provide a basis for further exploring the relationship between Fusobacterium sp. infection and carcinogenesis of GC.

Widespread Protection of RNA Cleavage Sites by a Riboswitch Aptamer that Folds as a Compact Obstacle to Scanning by RNase E

Riboswitches are thought generally to function by modulating transcription elongation or translation initiation. In rare instances, ligand binding to a riboswitch has been found to alter the rate of RNA degradation by directly stimulating or inhibiting nearby cleavage. Here, we show that guanidine-induced pseudoknot formation by the aptamer domain of a guanidine III riboswitch from Legionella pneumophila has a different effect, stabilizing mRNA by protecting distal cleavage sites en masse from ribonuclease attack. It does so by creating a coaxially base-paired obstacle that impedes scanning from a monophosphorylated 5' end to those sites by the regulatory endonuclease RNase E. Ligand binding by other riboswitch aptamers peripheral to the path traveled by RNase E does not inhibit distal cleavage. These findings reveal that a riboswitch aptamer can function independently of any overlapping expression platform to regulate gene expression by acting directly to prolong mRNA longevity in response to ligand binding.

Protein unties the pseudoknot: S1-mediated unfolding of RNA higher order structure

Ribosomal protein S1 plays important roles in the translation initiation step of many Escherichia coli mRNAs, particularly those with weak Shine-Dalgarno sequences or structured 5′ UTRs, in addition to a variety of cellular processes beyond the ribosome. In all cases, the RNA-binding activity of S1 is a central feature of its function. While sequence determinants of S1 affinity and many elements of the interactions of S1 with simple secondary structures are known, mechanistic details of the protein's interactions with RNAs of more complex secondary and tertiary structure are less understood. Here, we investigate the interaction of S1 with the well-characterized H-type pseudoknot of a class-I translational preQ₁ riboswitch as a highly structured RNA model whose conformation and structural dynamics can be tuned by the addition of ligands of varying binding affinity, particularly preQ₁, guanine, and 2,6-diaminopurine. Combining biochemical and single molecule fluorescence approaches, we show that S1 preferentially interacts with the less folded form of the pseudoknot and promotes a dynamic, partially unfolded conformation. The ability of S1 to unfold the RNA is inversely correlated with the structural stability of the pseudoknot. These mechanistic insights delineate the scope and limitations of S1-chaperoned unfolding of structured RNAs.

Circulating RNAs as Potential Biomarkers in Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a complex multi-system neurodegenerative disorder with currently limited diagnostic and no therapeutic options. Despite the intense efforts no clinically applicable biomarkers for ALS are yet established. Most current research is thus focused, in particular, in identifying potential non-invasive circulating biomarkers for more rapid and accurate diagnosis and monitoring of the disease. In this review, we have focused on messenger RNA (mRNA), non-coding RNAs (lncRNAs), micro RNAs (miRNAs) and circular RNA (circRNAs) as potential biomarkers for ALS in peripheral blood serum, plasma and cells. The most promising miRNAs include miR-206, miR-133b, miR-27a, mi-338-3p, miR-183, miR-451, let-7 and miR-125b. To test clinical potential of this miRNA panel, a useful approach may be to perform such analysis on larger multi-center scale using similar experimental design. However, other types of RNAs (lncRNAs, circRNAs and mRNAs) that, together with miRNAs, represent RNA networks, have not been yet extensively studied in blood samples of patients with ALS. Additional research has to be done in order to find robust circulating biomarkers and therapeutic targets that will distinguish key RNA interactions in specific ALS-types to facilitate diagnosis, predict progression and design therapy.

Conscious uncoupling of riboswitch functions

Riboswitches alter gene expression in response to ligand binding, coupling sensing and regulatory functions to help bacteria respond to their environment. The structural determinants of ligand binding in the prequeuosine (7-aminomethyl-7-deazaguanine, preQ₁) bacterial riboswitches have been studied, but the functional consequences of structural perturbations are less known. A new article combining biophysical and cell-based readouts of 15 mutants of the preQ₁-II riboswitch from Lactobacillus rhamnosus demonstrates that ligand binding does not ensure successful gene regulation, providing new insights into these shapeshifting sequences.

Exploring the RNA-Recognition Mechanism Using Supervised Molecular Dynamics (SuMD) Simulations: Toward a Rational Design for Ribonucleic-Targeting Molecules?

Although proteins have represented the molecular target of choice in the development of new drug candidates, the pharmaceutical importance of ribonucleic acids has gradually been growing. The increasing availability of structural information has brought to light the existence of peculiar three-dimensional RNA arrangements, which can, contrary to initial expectations, be recognized and selectively modulated through small chemical entities or peptides. The application of classical computational methodologies, such as molecular docking, for the rational development of RNA-binding candidates is, however, complicated by the peculiarities characterizing these macromolecules, such as the marked conformational flexibility, the singular charges distribution, and the relevant role of solvent molecules. In this work, we have thus validated and extended the applicability domain of SuMD, an all-atoms molecular dynamics protocol that allows to accelerate the sampling of molecular recognition events on a nanosecond timescale, to ribonucleotide targets of pharmaceutical interest. In particular, we have proven the methodological ability by reproducing the binding mode of viral or prokaryotic ribonucleic complexes, as well as that of artificially engineered aptamers, with an impressive degree of accuracy.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

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.

MicroRNA in Brain pathology: Neurodegeneration the Other Side of the Brain Cancer

The mammalian brain is made up of billions of neurons and supporting cells (glial cells), intricately connected. Molecular perturbations often lead to neurodegeneration by progressive loss of structure and malfunction of neurons, including their death. On the other side, a combination of genetic and cellular factors in glial cells, and less frequently in neurons, drive oncogenic transformation. In both situations, microenvironmental niches influence the progression of diseases and therapeutic responses. Dynamic changes that occur in cellular transcriptomes during the progression of developmental lineages and pathogenesis are controlled through a variety of regulatory networks. These include epigenetic modifications, signaling pathways, and transcriptional and post-transcriptional mechanisms. One prominent component of the latter is small non-coding RNAs, including microRNAs, that control the vast majority of these networks including genes regulating neural stemness, differentiation, apoptosis, projection fates, migration and many others. These cellular processes are also profoundly dependent on the microenvironment, stemness niche, hypoxic microenvironment, and interactions with associated cells including endothelial and immune cells. Significantly, the brain of all other mammalian organs expresses the highest number of microRNAs, with an additional gain in expression in the early stage of neurodegeneration and loss in expression in oncogenesis. However, a mechanistic explanation of the concept of an apparent inverse correlation between the odds of cancer and neurodegenerative diseases is only weakly developed. In this review, we thus will discuss widespread de-regulation of microRNAome observed in these two major groups of brain pathologies. The deciphering of these intricacies is of importance, as therapeutic restoration of pre-pathological microRNA landscape in neurodegeneration must not lead to oncogenesis and vice versa. We thus focus on microRNAs engaged in cellular processes that are inversely regulated in these diseases. We also aim to define the difference in microRNA networks between pro-survival and pro-apoptotic signaling in the brain.

Classification of riboswitch sequences using k-mer frequencies

Riboswitches are non-coding RNAs that regulate gene expression by altering the structural conformation of mRNA transcripts. Their regulation mechanism might be exploited for interesting biomedical applications such as drug targets and biosensors. A major challenge consists in accurately identifying metabolite-binding RNA switches which are structurally complex and diverse. In this regard, we investigated the classification of 16 riboswitch families using supervised learning algorithms trained solely with sequence-based features. We generated a reduced feature set and proposed a visual representation to explore its components. We induced Support Vector Machine, Random Forest, Naive Bayes, J48, and HyperPipes classifiers with our proposed feature set and tested their performance over independent data. Our best multi-class classifier achieved F-measure values of 0.996 and 0.966 in the training and test phases, respectively, outperforming those of a previous approach. When compared against BLAST, our best classifiers yielded competitive results. This work shows that the classifiers trained with our sequence-based feature set accurately discriminate riboswitches.

Biophysical properties, thermal stability and functional impact of 8-oxo-7,8-dihydroguanine on oligonucleotides of RNA-a study of duplex, hairpins and the aptamer for preQ1 as models

A better understanding of the effects that oxidative lesions have on RNA is of importance to understand their role in the development/progression of disease. 8-oxo-7,8-dihydroguanine was incorporated into RNA to understand its structural and functional impact on RNA:RNA and RNA:DNA duplexes, hairpins and pseudoknots. One to three modifications were incorporated into dodecamers of RNA [AAGAGGGAUGAC] resulting in thermal destabilization (ΔT_m – 10°C per lesion). Hairpins with tetraloops c-UUCG*-g* (8-10), a-ACCG-g* (11-12), c-UUG*G*-g* (13-16) and c-ACG*G*-g* (17-20) were modified and used to determine thermal stabilities, concluding that: (i) modifying the stem leads to destabilization unless adenosine is the opposing basepair of 8-oxoGua; (ii) modification at the loop is position- and sequence-dependent and varies from slight stabilization to large destabilization, in some cases leading to formation of other secondary structures (hairpin→duplex). Functional effects were established using the aptamer for preQ₁ as model. Modification at G5 disrupted the stem P1 and inhibited recognition of the target molecule 7-methylamino-7-deazaguanine (preQ₁). Modifying G11 results in increased thermal stability, albeit with a K_d 4-fold larger than its canonical analog. These studies show the capability of 8-oxoG to affect structure and function of RNA, resulting in distinct outcomes as a function of number and position of the lesion.

NMR monitoring of the SELEX process to confirm enrichment of structured RNA

RNA aptamers are RNA molecules that bind to a target molecule with high affinity and specificity using uniquely-folded tertiary structures. RNA aptamers are selected from an RNA pool typically comprising up to 1015 different sequences generated by iterative steps of selection and amplification known as Systematic Evolution of Ligands by EXponential enrichment (SELEX). Over several rounds of SELEX, the diversity of the RNA pool decreases and the aptamers are enriched. Hence, monitoring of the enrichment of these RNA pools is critical for the successful selection of aptamers, and several methods for monitoring them have been developed. In this study, we measured one-dimensional imino proton NMR spectra of RNA pools during SELEX. The spectrum of the initial RNA pool indicates that the RNAs adopt tertiary structures. The structural diversity of the RNA pools was shown to depend highly on the design of the primer-binding sequence. Furthermore, we demonstrate that enrichment of RNA aptamers can be monitored using NMR. The RNA pools can be recovered from the NMR tube after measurement of NMR spectra. We also can monitor target binding in the NMR tubes. Thus, we propose using NMR to monitor the enrichment of structured aptamers during the SELEX process.

Characterization of Engineered PreQ1 Riboswitches for Inducible Gene Regulation in Mycobacteria

We report here the behavior of naturally occurring and rationally engineered preQ1 riboswitches and their application to inducible gene regulation in mycobacteria. Because mycobacteria lack preQ1 biosynthetic genes, we hypothesized that preQ1 could be used as an exogenous nonmetabolite ligand to control riboswitches in mycobacteria. Selected naturally occurring preQ1 riboswitches were assayed and successfully drove preQ1-dependent repression of a green fluorescent protein reporter in Mycobacterium smegmatis Using structure-based design, we engineered three preQ1 riboswitches from Thermoanaerobacter tencongensis, Bacillus subtilis, and Lactobacillus rhamnosus toward achieving higher response ratios and increased repression. Assuming a steady-state model, variants of the T. tencongensis riboswitch most closely followed the predicted trends. Unexpectedly, the preQ1 dose response was best described by a model with a second, independent preQ1 binding site. This behavior was general to the preQ1 riboswitch family, since the wild type and rationally designed mutants of riboswitches from all three bacteria behaved analogously. Across all variants, the response ratios, which describe expression in the absence versus the presence of preQ1, ranged from <2 to ∼10, but repression in all cases was incomplete up to 1 mM preQ1. By reducing the transcript expression level, we obtained a preQ1 riboswitch variant appropriate for inducible knockdown applications. We further showed that the preQ1 response is reversible, is titratable, and can be used to control protein expression in mycobacteria within infected macrophages. By engineering naturally occurring preQ1 riboswitches, we have not only extended the tools available for inducible gene regulation in mycobacteria but also uncovered new behavior of these riboswitches.IMPORTANCE Riboswitches are elements found in noncoding regions of mRNA that regulate gene expression, typically in response to an endogenous metabolite. Riboswitches have emerged as important tools for inducible gene expression in diverse organisms. We noted that mycobacteria lack the biosynthesis genes for preQ1, a ligand for riboswitches from diverse bacteria. Predicting that preQ1 is not present in mycobacteria, we showed that it controls optimized riboswitches appropriate for gene knockdown applications. Further, the riboswitch response is subject to a second independent preQ1 binding event that has not been previously documented. By engineering naturally occurring riboswitches, we have uncovered a new behavior, with implications for riboswitch function in its native context, and extended the tools available for inducible gene regulation in mycobacteria.

The Shine-Dalgarno sequence of riboswitch-regulated single mRNAs shows ligand-dependent accessibility bursts

In response to intracellular signals in Gram-negative bacteria, translational riboswitches—commonly embedded in messenger RNAs (mRNAs)—regulate gene expression through inhibition of translation initiation. It is generally thought that this regulation originates from occlusion of the Shine-Dalgarno (SD) sequence upon ligand binding; however, little direct evidence exists. Here we develop Single Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) to investigate the ligand-dependent accessibility of the SD sequence of an mRNA hosting the 7-aminomethyl-7-deazaguanine (preQ₁)-sensing riboswitch. Spike train analysis reveals that individual mRNA molecules alternate between two conformational states, distinguished by ‘bursts' of probe binding associated with increased SD sequence accessibility. Addition of preQ₁ decreases the lifetime of the SD's high-accessibility (bursting) state and prolongs the time between bursts. In addition, ligand-jump experiments reveal imperfect riboswitching of single mRNA molecules. Such complex ligand sensing by individual mRNA molecules rationalizes the nuanced ligand response observed during bulk mRNA translation.

In response to intracellular signals, bacterial translational riboswitches embedded in mRNAs can regulate gene expression through inhibition of translation initiation. Here, the authors describe SiM-KARTS, a novel approach for detecting changes in the structure of single RNA molecules in response to a ligand.

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.

Mg(2+) shifts ligand-mediated folding of a riboswitch from induced-fit to conformational selection

Bacterial riboswitches couple small-molecule ligand binding to RNA conformational changes that widely regulate gene expression, rendering them potential targets for antibiotic intervention. Despite structural insights, the ligand-mediated folding mechanisms of riboswitches are still poorly understood. Using single-molecule fluorescence resonance energy transfer (smFRET), we have investigated the folding mechanism of an H-type pseudoknotted preQ1 riboswitch in dependence of Mg(2+) and three ligands of distinct affinities. We show that, in the absence of Mg(2+), both weakly and strongly bound ligands promote pseudoknot docking through an induced-fit mechanism. By contrast, addition of as low as 10 μM Mg(2+) generally shifts docking toward conformational selection by stabilizing a folded-like conformation prior to ligand binding. Supporting evidence from transition-state analysis further highlights the particular importance of stacking interactions during induced-fit and of specific hydrogen bonds during conformational selection. Our mechanistic dissection provides unprecedented insights into the intricate synergy between ligand- and Mg(2+)-mediated RNA folding.

Complete Genome Sequence of Enterotoxigenic Escherichia coli Siphophage Seurat

Enterotoxigenic Escherichia coli (ETEC) is one of the leading causes of diarrhea in developing countries. Bacteriophage therapy has the potential to aid in the prevention and treatment of ETEC-related illness. To that end, we present here the complete genome of ETEC siphophage Seurat and describe its major features.