Riboswitches: structures and mechanisms

SPRank─A Knowledge-Based Scoring Function for RNA-Ligand Pose Prediction and Virtual Screening

The growing interest in RNA-targeted drugs underscores the need for computational modeling of interactions between RNA molecules and small compounds. Having a reliable scoring function for RNA-ligand interactions is essential for effective computational drug screening. An ideal scoring function should not only predict the native pose for ligand binding but also rank the affinity of the binding for different ligands. However, existing scoring functions are primarily designed to predict the native binding modes for a given RNA-ligand pair and have not been thoroughly assessed for virtual screening purposes. In this paper, we introduce SPRank, a combination of machine-learning and knowledge-based scoring functions developed through a weighted iterative approach, specifically designed to tackle both binding mode prediction and virtual screening challenges. Our approach incorporates third-party docking software, such as rDock and AutoDock Vina, to sample flexible ligands against an ensemble of RNA structures, capturing the conformational flexibility of both the RNA and the ligand. Through rigorous testing, SPRank demonstrates improved performance compared to the tested scoring functions across four test sets comprising 122, 42, 55, and 71 nucleic acid-ligand complexes. Furthermore, SPRank exhibits improved performance in virtual screening tests targeting the HIV-1 TAR ensemble, which highlights its advantage in drug discovery. These results underscore the advantages of SPRank as a potentially promising tool for the RNA-targeted drug design. The source code of SPRank and the data sets are freely accessible at https://github.com/Vfold-RNA/SPRank.

Structure-based characterization and compound identification of the wild-type THF class-II riboswitch

Riboswitches are conserved regulatory RNA elements participating in various metabolic pathways. Recently, a novel RNA motif known as the folE RNA motif was discovered upstream of folE genes. It specifically senses tetrahydrofolate (THF) and is therefore termed THF-II riboswitch. To unravel the ligand recognition mechanism of this newly discovered riboswitch and decipher the underlying principles governing its tertiary folding, we determined both the free-form and bound-form THF-II riboswitch in the wild-type sequences. Combining structural information and isothermal titration calorimetry (ITC) binding assays on structure-based mutants, we successfully elucidated the significant long-range interactions governing the function of THF-II riboswitch and identified additional compounds, including alternative natural metabolites and potential lead compounds for drug discovery, that interact with THF-II riboswitch. Our structural research on the ligand recognition mechanism of the THF-II riboswitch not only paves the way for identification of compounds targeting riboswitches, but also facilitates the exploration of THF analogs in diverse biological contexts or for therapeutic applications.

START: A Versatile Platform for Bacterial Ligand Sensing with Programmable Performances

Recognition of signaling molecules for coordinated regulation of target genes is a fundamental process for biological systems. Cells often rely on transcription factors to accomplish these intricate tasks, yet the subtle conformational changes of protein structures, coupled with the complexity of intertwined protein interaction networks, pose challenges for repurposing these for bioengineering applications. This study introduces a novel platform for ligand-responsive gene regulation, termed START (Synthetic Trans-Acting Riboswitch with Triggering RNA). Inspired by the bacterial ligand sensing system, riboswitch, and the synthetic gene regulator, toehold switch, the START platform enables the implementation of synthetic biosensors for various ligands. Rational sequence design with targeted domain optimization yields high-performance STARTs with a dynamic range up to 67.29-fold and a tunable ligand sensitivity, providing a simple and intuitive strategy for sensor engineering. The START platform also exhibits modularity and composability to allow flexible genetic circuit construction, enabling seamless implementation of OR, AND, and NOT Boolean logic gates for multiple ligand inputs. The START design principle is capable of broadening the suite of synthetic biosensors for diverse chemical and protein ligands, providing a novel riboregulator chassis for synthetic biology and bioengineering applications.

Hfq mediates transcriptome-wide RNA structurome reprogramming under virulence-inducing conditions in a phytopathogen

Although RNA structures play important roles in regulating gene expression, the mechanism and function of mRNA folding in plant bacterial pathogens remain elusive. Therefore, we perform dimethyl sulfate sequencing (DMS-seq) on the Pseudomonas syringae under nutrition-rich and -deficient conditions, revealing that the mRNA structure changes substantially in the minimal medium (MM) that tunes global translation efficiency (TE), thereby inducing virulence. This process is led by the increased expression of hfq, which is directly activated by transcription regulators RpoS and CysB. The co-occurrence of Hfq and RpoS in diverse bacteria and the deep conservation of Hfq Y25 is critical for RNA-mediated regulation and implicates the wider biological importance of mRNA structure and feedback loops in the control of global gene expression.

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.

Therapeutic Applications of Aptamers

Affinity reagents, or target-binding molecules, are quite versatile and are major workhorses in molecular biology and medicine. Antibodies are the most famous and frequently used type and they have been used for a wide range of applications, including laboratory techniques, diagnostics, and therapeutics. However, antibodies are not the only available affinity reagents and they do have significant drawbacks, including laborious and costly production. Aptamers are one potential alternative that have a variety of unique advantages. They are single stranded DNA or RNA molecules that can be selected for binding to many targets including proteins, carbohydrates, and small molecules-for which antibodies typically have low affinity. There are also a variety of cost-effective methods for producing and modifying nucleic acids in vitro without cells, whereas antibodies typically require cells or even whole animals. While there are also significant drawbacks to using aptamers in therapeutic applications, including low in vivo stability, aptamers have had success in clinical trials for treating a variety of diseases and two aptamer-based drugs have gained FDA approval. Aptamer development is still ongoing, which could lead to additional applications of aptamer therapeutics, including antitoxins, and combinatorial approaches with nanoparticles and other nucleic acid therapeutics that could improve efficacy.

Probing RNA structure and dynamics using nanopore and next generation sequencing

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

Skin-associated Corynebacterium amycolatum shares cobamides

The underlying interactions that occur to maintain skin microbiome composition, function, and overall skin health are largely unknown. Often, these types of interactions are mediated by microbial metabolites. Cobamides, the vitamin B12 family of cofactors, are essential for metabolism in many bacteria, but are only synthesized by a small fraction of prokaryotes, including certain skin-associated species. Therefore, we hypothesize that cobamide sharing mediates skin community dynamics. Preliminary work predicts that several skin-associated Corynebacterium species encode de novo cobamide biosynthesis and that their abundance is associated with skin microbiome diversity. Here, we show that commensal Corynebacterium amycolatum produces cobamides and that this synthesis can be tuned by cobalt limitation. To demonstrate cobamide sharing by C. amycolatum, we employed a co-culture assay using an E. coli cobamide auxotroph and show that C. amycolatum produces sufficient cobamides to support E. coli growth, both in liquid co-culture and when separated spatially on solid medium. We also generated a C. amycolatum non-cobamide-producing strain (cob-) using UV mutagenesis that contains mutated cobamide biosynthesis genes cobK and cobO and confirm that disruption of cobamide biosynthesis abolishes support of E. coli growth through cobamide sharing. Our study provides a unique model to study metabolite sharing by microorganisms, which will be critical for understanding the fundamental interactions that occur within complex microbiomes and for developing approaches to target the human microbiota for health advances.

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.

Structured 3' UTRs destabilize mRNAs in plants

BACKGROUND

RNA secondary structure (RSS) can influence the regulation of transcription, RNA processing, and protein synthesis, among other processes. 3' untranslated regions (3' UTRs) of mRNA also hold the key for many aspects of gene regulation. However, there are often contradictory results regarding the roles of RSS in 3' UTRs in gene expression in different organisms and/or contexts.

RESULTS

Here, we incidentally observe that the primary substrate of miR159a (pri-miR159a), when embedded in a 3' UTR, could promote mRNA accumulation. The enhanced expression is attributed to the earlier polyadenylation of the transcript within the hybrid pri-miR159a-3' UTR and, resultantly, a poorly structured 3' UTR. RNA decay assays indicate that poorly structured 3' UTRs could promote mRNA stability, whereas highly structured 3' UTRs destabilize mRNA in vivo. Genome-wide DMS-MaPseq also reveals the prevailing inverse relationship between 3' UTRs' RSS and transcript accumulation in the transcriptomes of Arabidopsis, rice, and even human. Mechanistically, transcripts with highly structured 3' UTRs are preferentially degraded by 3'-5' exoribonuclease SOV and 5'-3' exoribonuclease XRN4, leading to decreased expression in Arabidopsis. Finally, we engineer different structured 3' UTRs to an endogenous FT gene and alter the FT-regulated flowering time in Arabidopsis.

CONCLUSIONS

We conclude that highly structured 3' UTRs typically cause reduced accumulation of the harbored transcripts in Arabidopsis. This pattern extends to rice and even mammals. Furthermore, our study provides a new strategy of engineering the 3' UTRs' RSS to modify plant traits in agricultural production and mRNA stability in biotechnology.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows either or both states are highly conserved. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

Harvesting phosphorus-containing moieties for their antibacterial effects

Clinically manifested resistance of bacteria to antibiotics has emerged as a global threat to society and there is an urgent need for the development of novel classes of antibacterial agents. Recently, the use of phosphorus in antibacterial agents has been explored in quite an unprecedent manner. In this comprehensive review, we summarize the use of phosphorus-containing moieties (phosphonates, phosphonamidates, phosphonopeptides, phosphates, phosphoramidates, phosphinates, phosphine oxides, and phosphoniums) in compounds with antibacterial effect, including their use as β-lactamase inhibitors and antibacterial disinfectants. We show that phosphorus-containing moieties can serve as novel pharmacophores, bioisosteres, and prodrugs to modify pharmacodynamic and pharmacokinetic properties. We further discuss the mechanisms of action, biological activities, clinical use and highlight possible future prospects.

Elimination of editing plasmid mediated by theophylline riboswitch in Zymomonas mobilis

Zymomonas mobilis is regarded as a potential chassis for the production of platform chemicals. Genome editing using the CRISPR-Cas system could meet the need for gene modification in metabolic engineering. However, the low curing efficiency of CRISPR editing plasmid is a common bottleneck in Z. mobilis. In this study, we utilized a theophylline-dependent riboswitch to regulate the expression of the replicase gene of the editing plasmid, thereby promoting the elimination of exogeneous plasmid. The riboswitch D (RSD) with rigorous regulatory ability was identified as the optimal candidate by comparing the transformation efficiency of four theophylline riboswitch-based backbone editing plasmids, and the optimal theophylline concentration for inducing RSD was determined to be 2 mM. A highly effective method for eliminating the editing plasmid, cells with RSD-based editing plasmid which were cultured in liquid and solid RM media in alternating passages at 37 °C without shaking, was established by testing the curing efficiency of backbone editing plasmids pMini and pMini-RSD in RM medium with or without theophylline at 30 °C or 37 °C. Finally, the RSD-based editing plasmid was applied to genome editing, resulting in an increase of more than 10% in plasmid elimination efficiency compared to that of pMini-based editing plasmid. KEY POINTS: • An effective strategy for curing CRISPR editing plasmid has been established in Z. mobilis. • Elimination efficiency of the CRISPR editing plasmid was enhanced by 10% to 20% under the regulation of theophylline-dependent riboswitch RSD.

AliNA - a deep learning program for RNA secondary structure prediction

Nowadays there are numerous discovered natural RNA variations participating in different cellular processes and artificial RNA, e. g., aptamers, riboswitches. One of the required tasks in the investigation of their functions and mechanism of influence on cells and interaction with targets is the prediction of RNA secondary structures. The classic thermodynamic-based prediction algorithms do not consider the specificity of biological folding and deep learning methods that were designed to resolve this issue suffer from homology-based methods problems. Herein, we present a method for RNA secondary structure prediction based on deep learning - AliNA (ALIgned Nucleic Acids). Our method successfully predicts secondary structures for non-homologous to train-data RNA families thanks to usage of the data augmentation techniques. Augmentation extends existing datasets with easily-accessible simulated data. The proposed method shows a high quality of prediction across different benchmarks including pseudoknots. The method is available on GitHub for free (https://github.com/Arty40m/AliNA).

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

A transient intermediate RNA structure underlies the regulatory function of the E. coli thiB TPP translational riboswitch

Riboswitches are cis-regulatory RNA elements that regulate gene expression in response to ligand binding through the coordinated action of a ligand-binding aptamer domain (AD) and a downstream expression platform (EP). Previous studies of transcriptional riboswitches have uncovered diverse examples that utilize structural intermediates that compete with the AD and EP folds to mediate the switching mechanism on the timescale of transcription. Here we investigate whether similar intermediates are important for riboswitches that control translation by studying the Escherichia coli thiB thiamin pyrophosphate (TPP) riboswitch. Using cellular gene expression assays, we first confirmed that the riboswitch acts at the level of translational regulation. Deletion mutagenesis showed the importance of the AD-EP linker sequence for riboswitch function. Sequence complementarity between the linker region and the AD P1 stem suggested the possibility of an intermediate nascent RNA structure called the antisequestering stem that could mediate the thiB switching mechanism. Experimentally informed secondary structure models of the thiB folding pathway generated from chemical probing of nascent thiB structures in stalled transcription elongation complexes confirmed the presence of the antisequestering stem, and showed it may form cotranscriptionally. Additional mutational analysis showed that mutations to the antisequestering stem break or bias thiB function according to whether the antisequestering stem or P1 is favored. This work provides an important example of intermediate structures that compete with AD and EP folds to implement riboswitch mechanisms.

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.

Riboswitches, from cognition to transformation

Riboswitches are functional RNA elements that regulate gene expression by directly detecting metabolites. Twenty years have passed since it was first discovered, researches on riboswitches are becoming increasingly standardized and refined, which could significantly promote people's cognition of RNA function as well. Here, we focus on some representative orphan riboswitches, enumerate the structural and functional transformation and artificial design of riboswitches including the coupling with ribozymes, hoping to attain a comprehensive understanding of riboswitch research.

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.

Visualizing a two-state conformational ensemble in stem-loop 3 of the transcriptional regulator 7SK RNA

Structural plasticity is integral to RNA function; however, there are currently few methods to quantitatively resolve RNAs that have multiple structural states. NMR spectroscopy is a powerful approach for resolving conformational ensembles but is size-limited. Chemical probing is well-suited for large RNAs but provides limited structural and no kinetics information. Here, we integrate the two approaches to visualize a two-state conformational ensemble for the central stem-loop 3 (SL3) of 7SK RNA, a critical element for 7SK RNA function in transcription regulation. We find that the SL3 distal end exchanges between two equally populated yet structurally distinct states in both isolated SL3 constructs and full-length 7SK RNA. We rationally designed constructs that lock SL3 into a single state and demonstrate that both chemical probing and NMR data fit to a linear combination of the two states. Comparison of vertebrate 7SK RNA sequences shows conservation of both states, suggesting functional importance. These results provide new insights into 7SK RNA structural dynamics and demonstrate the utility of integrating chemical probing with NMR spectroscopy to gain quantitative insights into RNA conformational ensembles.

RNA folding pathways from all-atom simulations with a variationally improved history-dependent bias

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. As a step forward in tackling these issues, we extend to RNA an enhanced path-sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force promoting the generation of productive folding trajectories in an all-atom force field with explicit solvent. A rigorous variational principle is then applied to minimize the effect of the bias. Here, we report on an application of this method to RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. By comparison with analog simulations performed on small proteins with similar size and architecture, we show that the RNA folding landscape is significantly more frustrated, even for relatively small chains with a simple topology. The predicted RNA folding mechanisms are found to be consistent with the available experiments and some of the existing coarse-grained models. Due to its computational performance, this scheme provides a promising platform to efficiently gather atomistic RNA folding trajectories, thus retain the information about the chemical composition of the sequence.

Translation regulation by a guanidine-II riboswitch is highly tunable in sensitivity, dynamic range, and apparent cooperativity

Riboswitches function as important translational regulators in bacteria. Comprehensive mutational analysis of transcriptional riboswitches has been used to probe the energetic intricacies of interplay between the aptamer and expression platform, but translational riboswitches have been inaccessible to massively parallel techniques. The guanidine-II (gdm-II) riboswitch is an exclusively translational class. We have integrated RelE cleavage with next-generation sequencing to quantify ligand-dependent changes in translation initiation for all single and double mutations of the Pseudomonas aeruginosa gdm-II riboswitch, a total of more than 23,000 variants. This extensive mutational analysis is consistent with the prominent features of the bioinformatic consensus. These data indicate, unexpectedly, that direct sequestration of the Shine-Dalgarno sequence is dispensable for riboswitch function. Additionally, this comprehensive data set reveals important positions not identified in previous computational and crystallographic studies. Mutations in the variable linker region stabilize alternate conformations. The double mutant data reveal the functional importance of the previously modeled P0b helix formed by the 5' and 3' tails that serves as the basis for translational control. Additional mutations to GU wobble base pairs in both P1 and P2 reveal how the apparent cooperativity of the system involves an intricate network of communication between the two binding sites. This comprehensive examination of a translational riboswitch's expression platform illuminates how the riboswitch is precisely tuned and tunable with regard to ligand sensitivity, the amplitude of expression between ON and OFF states, and the cooperativity of ligand binding.

Exploring the structure, function of thiamine pyrophosphate riboswitch, and designing small molecules for antibacterial activity

During the last decade, riboswitches emerged as new small-molecule sensing RNA in bacteria. Thiamine pyrophosphate (TPP) riboswitch is widely distributed and occurs in plants, bacteria, fungi, and archaea. Extensive biochemical, structural, and genetic studies have been carried out to elucidate the recognition mechanism of TPP riboswitches. However, a comprehensive report summarizing all information on recognition principles and newly designed ligands for TPP riboswitch is scarce in the literature. This review gives a comprehensive understanding of the TPP riboswitch's structure, mechanism, and methods applied to design ligands for the TPP riboswitch. The ligand-bound TPP riboswitch was studied with various experimental and theoretical techniques to elucidate the conformational dynamics. The mutation studies shed light on the significance of pyrimidine sensing helix for the binding of ligands. Further, the structure-activity relationship study and fragment-based approach lead to the development of ligands with Kd values at the sub-micromolar level. However, there is a need to design more potent inhibitors for TPP riboswitch for therapeutic applications. The recent advancements in ligand design highlight the TPP riboswitch as a promising target for developing new antibiotics. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Riboswitches Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K⁺, providing evidence for synergism between lysine/K⁺ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K⁺. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K⁺ binding, revealing that riboswitch folding (k_fold) and unfolding (k_unfold) rate constants increase and decrease dramatically with K⁺, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K⁺ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K⁺ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na⁺ and Cs⁺ studies to demonstrate that K⁺ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K⁺ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

A Short Review on Miniaturized Biosensors for the Detection of Nucleic Acid Biomarkers

Even today, most biomarker testing is executed in centralized, dedicated laboratories using bulky instruments, automated analyzers, and increased analysis time and expenses. The development of miniaturized, faster, low-cost microdevices is immensely anticipated for substituting for these conventional laboratory-oriented assays and transferring diagnostic results directly onto the patient's smartphone using a cloud server. Pioneering biosensor-based approaches might make it possible to test biomarkers with reliability in a decentralized setting, but there are still a number of issues and restrictions that must be resolved before the development and use of several biosensors for the proper understanding of the measured biomarkers of numerous bioanalytes such as DNA, RNA, urine, and blood. One of the most promising processes to address some of the issues relating to the growing demand for susceptible, quick, and affordable analysis techniques in medical diagnostics is the creation of biosensors. This article critically discusses a short review of biosensors used for detecting nucleic acid biomarkers, and their use in biomedical prognostics will be addressed while considering several essential characteristics.

Current Metabolic Engineering Strategies for Photosynthetic Bioproduction in Cyanobacteria

Cyanobacteria are photosynthetic microorganisms capable of using solar energy to convert CO₂ and H₂O into O₂ and energy-rich organic compounds, thus enabling sustainable production of a wide range of bio-products. More and more strains of cyanobacteria are identified that show great promise as cell platforms for the generation of bioproducts. However, strain development is still required to optimize their biosynthesis and increase titers for industrial applications. This review describes the most well-known, newest and most promising strains available to the community and gives an overview of current cyanobacterial biotechnology and the latest innovative strategies used for engineering cyanobacteria. We summarize advanced synthetic biology tools for modulating gene expression and their use in metabolic pathway engineering to increase the production of value-added compounds, such as terpenoids, fatty acids and sugars, to provide a go-to source for scientists starting research in cyanobacterial metabolic engineering.

Toehold-Exchange-Based Activation of Aptamer Switches Enables High Thermal Robustness and Programmability

Aptamer switches are attractive nature-inspired tools for developing smart materials and nanodevices. However, the thermal robustness and programmability of current aptamer switches are often limited by their activation processes that are coupled with high reaction enthalpy. Here, we present an enthalpy-independent activation approach that harnesses toehold-exchange as a general framework to design aptamer switches. We demonstrate mathematically and experimentally that this approach is highly effective in improving thermal robustness and thus leads to better analytical performances of aptamer switches. Enhanced programmability is also demonstrated through fine-grained and dynamic tuning of effective affinities and dynamic ranges, as well as the construction of a synthetic DNA network that resembled biological signaling cascades. Our study not only enriches the current toolbox for engineering and controlling synthetic molecular switches but also offers new insights into their thermodynamic basis, which is critical for diverse synthetic biological designs and applications.

RLBind: a deep learning method to predict RNA-ligand binding sites

Identification of RNA-small molecule binding sites plays an essential role in RNA-targeted drug discovery and development. These small molecules are expected to be leading compounds to guide the development of new types of RNA-targeted therapeutics compared with regular therapeutics targeting proteins. RNAs can provide many potential drug targets with diverse structures and functions. However, up to now, only a few methods have been proposed. Predicting RNA-small molecule binding sites still remains a big challenge. New computational model is required to better extract the features and predict RNA-small molecule binding sites more accurately. In this paper, a deep learning model, RLBind, was proposed to predict RNA-small molecule binding sites from sequence-dependent and structure-dependent properties by combining global RNA sequence channel and local neighbor nucleotides channel. To our best knowledge, this research was the first to develop a convolutional neural network for RNA-small molecule binding sites prediction. Furthermore, RLBind also can be used as a potential tool when the RNA experimental tertiary structure is not available. The experimental results show that RLBind outperforms other state-of-the-art methods in predicting binding sites. Therefore, our study demonstrates that the combination of global information for full-length sequences and local information for limited local neighbor nucleotides in RNAs can improve the model's predictive performance for binding sites prediction. All datasets and resource codes are available at https://github.com/KailiWang1/RLBind.

Structure-based investigations of the NAD+-II riboswitch

Riboswitches are conserved non-coding domains in bacterial mRNA with gene regulation function that are essential for maintaining enzyme co-factor metabolism. Recently, the pnuC RNA motif was reported to selectively bind nicotinamide adenine dinucleotide (NAD+), defining a novel class of NAD+ riboswitches (NAD+-II) according to phylogenetic analysis. To reveal the three-dimensional architecture and the ligand-binding mode of this riboswitch, we solved the crystal structure of NAD+-II riboswitch in complex with NAD+. Strikingly and in contrast to class-I riboswitches that form a tight recognition pocket for the adenosine diphosphate (ADP) moiety of NAD+, the class-II riboswitches form a binding pocket for the nicotinamide mononucleotide (NMN) portion of NAD+ and display only unspecific interactions with the adenosine. We support this finding by an additional structure of the class-II RNA in complex with NMN alone. The structures define a novel RNA tertiary fold that was further confirmed by mutational analysis in combination with isothermal titration calorimetry (ITC), and 2-aminopurine-based fluorescence spectroscopic folding studies. Furthermore, we truncated the pnuC RNA motif to a short RNA helical scaffold with binding affinity comparable to the wild-type motif to allude to the potential of engineering the NAD+-II motif for biotechnological applications.

Virus Evolution on Fitness Landscapes

The landscape paradigm is revisited in the light of evolution in simple systems. A brief overview of different classes of fitness landscapes is followed by a more detailed discussion of the RNA model, which is currently the only evolutionary model that allows for a comprehensive molecular analysis of a fitness landscape. Neutral networks of genotypes are indispensable for the success of evolution. Important insights into the evolutionary mechanism are gained by considering the topology of sequence and shape spaces. The dynamic concept of molecular quasispecies is viewed in the light of the landscape paradigm. The distribution of fitness values in state space is mirrored by the population structures of mutant distributions. Two classes of thresholds for replication error or mutations are important: (i) the-conventional-genotypic error threshold, which separates ordered replication from random drift on neutral networks, and (ii) a phenotypic error threshold above which the molecular phenotype is lost. Empirical landscapes are reviewed and finally, the implications of the landscape concept for virus evolution are discussed.

Riboswitches as therapeutic targets: promise of a new era of antibiotics

INTRODUCTION

The growth of antibiotic resistance among bacterial pathogens is an impending global threat that can only be averted through the development of novel antibacterial drugs. A promising answer could be the targeting of riboswitches, structured RNA elements found almost exclusively in bacteria.

AREAS COVERED

This review examines the potential of riboswitches as novel antibacterial drug targets. The limited mechanisms of action of currently available antibiotics are summarized, followed by a delineation of the functional mechanisms of riboswitches. We then discuss the potential for developing novel approaches that target paradigmatic riboswitches in the context of their bacterial gene expression machinery.

EXPERT OPINION

We highlight potential advantages of targeting riboswitches in their functional form, embedded within gene expression complexes critical for bacterial survival. We emphasize the benefits of this approach, including potentially higher species specificity and lower side effects.

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.

Rational selection of an ideal oncolytic virus to address current limitations in clinical translation

Oncolytic virus therapy (OVT) is a promising modality that leverages the propensity of natural or engineered viruses to selectively replicate in and kill cancer cells. Over the past decade, (pre)clinical studies have focused on the development and testing of adenovirus, herpes simplex virus, and vaccinia virus-based vectors. These studies have identified barriers to success confronting the field. Here, we propose a set of selection criteria or ideal properties of a successful oncolytic virus, which include lack of pathogenicity, low seroprevalence, selectivity (infection and replication), transgene carrying capacity, and genome stability. We use these requirements to analyze the oncolytic virus landscape, and then identify a potentially optimal species for platform development - vesicular stomatitis virus.

Development of DNA aptamer selection approach based on membrane ultrafiltration of aptamer/target complex

Background. Aptamers are small single-stranded DNA or RNA molecules that have an affinity for a specific target molecule. The main method of aptamers construction is the technology of systematic evolution of ligands with exponential enrichment (SELEX). However, the exact approach depends on the nature of target molecules, and is selected and optimized by each researcher independently. The article describes the technique of production of aptamers to the tick-borne encephalitis virus (TBEV) using membrane ultrafiltration with a molecular weight cut-off of 100 kDa. As a result, the pool of aptamers with observable affinity for TBEV is successfully selected and enriched.The aim. To develop the technique suitable for selection of specific DNA aptamers to a live, crude TBEV suspension directly in cell culture supernatant.Materials and methods. The selection of aptamers was carried out using a modified SELEX DNA aptamer technology in combination with semipermeable membrane ultrafiltration using Vivaspin 6 (Sartorius, Germany) concentrators of molecular weight cut-off of 100 kDa. Enrichment of a specific pool of aptamers was performed using real time polymerase chain reaction. Aptamers were sequenced with automated Sanger sequencing method. The direct virucidal effect of the aptamers was determined by the decrease in the titer of the infectious virus after incubation with the aptamer.Results. The pool of aptamers to TBEV was selected and enriched. This aptamer pool expressed affinity both to the infectious TBEV and to the TBEV antigen. Sixteen aptamers were sequenced from this pool and four of them were synthesized and tested for antiviral activity against TBEV. No antiviral activity was observed.Conclusions. The technique developed that can be successfully used to select aptamers to a live virus culture for the viruses comparable in size to TBEV or larger.

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.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.

Cobamide Sharing Is Predicted in the Human Skin Microbiome

The skin microbiome is a key player in human health, with diverse functions ranging from defense against pathogens to education of the immune system. While recent studies have begun to shed light on the valuable role that skin microorganisms have in maintaining the skin barrier, a detailed understanding of the complex interactions that shape healthy skin microbial communities is limited. Cobamides, the vitamin B₁₂ class of cofactor, are essential for organisms across the tree of life. Because this vitamin is only produced by a limited fraction of prokaryotes, cobamide sharing is predicted to mediate community dynamics within microbial communities. Here, we provide the first large-scale metagenomic assessment of cobamide biosynthesis and utilization in the skin microbiome. We show that while numerous and diverse taxa across the major bacterial phyla on the skin encode cobamide-dependent enzymes, relatively few species encode de novo cobamide biosynthesis. We show that cobamide producers and users are integrated into the network structure of microbial communities across the different microenvironments of the skin and that changes in microbiome community structure and diversity are associated with the abundance of cobamide producers in the Corynebacterium genus, for both healthy and diseased skin states. Finally, we find that de novo cobamide biosynthesis is enriched only in Corynebacterium species associated with hosts, including those prevalent on human skin. We confirm that the cofactor is produced in excess through quantification of cobamide production by human skin-associated species isolated in the laboratory. Taken together, our results reveal the potential for cobamide sharing within skin microbial communities, which we hypothesize mediates microbiome community dynamics and host interactions. IMPORTANCE The skin microbiome is essential for maintaining skin health and function. However, the microbial interactions that dictate microbiome structure, stability, and function are not well understood. Here, we investigate the biosynthesis and use of cobamides, a cofactor needed by many organisms but only produced by select prokaryotes, within the human skin microbiome. We found that while a large proportion of skin taxa encode cobamide-dependent enzymes, only a select few encode de novo cobamide biosynthesis. Further, the abundance of cobamide-producing Corynebacterium species is associated with skin microbiome diversity and structure, and within this genus, de novo biosynthesis is enriched in host-associated species compared to environment-associated species. These findings identify cobamides as a potential mediator of skin microbiome dynamics and skin health.

Targeting RNA structures with small molecules

RNA adopts 3D structures that confer varied functional roles in human biology and dysfunction in disease. Approaches to therapeutically target RNA structures with small molecules are being actively pursued, aided by key advances in the field including the development of computational tools that predict evolutionarily conserved RNA structures, as well as strategies that expand mode of action and facilitate interactions with cellular machinery. Existing RNA-targeted small molecules use a range of mechanisms including directing splicing - by acting as molecular glues with cellular proteins (such as branaplam and the FDA-approved risdiplam), inhibition of translation of undruggable proteins and deactivation of functional structures in noncoding RNAs. Here, we describe strategies to identify, validate and optimize small molecules that target the functional transcriptome, laying out a roadmap to advance these agents into the next decade.

Riboswitch-mediated regulation of riboflavin biosynthesis genes in prokaryotes

Prokaryotic organisms frequently use riboswitches to quantify intracellular metabolite concentration via high-affinity metabolite receptors. Riboswitches possess a metabolite-sensing system that controls gene regulation in a cis-acting fashion at the initiation of transcriptional/translational level by binding with a specific metabolite and controlling various biochemical pathways. Riboswitch binds with flavin mononucleotide (FMN), a phosphorylated form of riboflavin and controls gene expression involved in riboflavin biosynthesis and transport pathway. The first step of the riboflavin biosynthesis pathway is initiated by the conversion of guanine nucleotide triphosphate (GTP), which is an intermediate of the purine biosynthesis pathway. An alternative pentose phosphate pathway of riboflavin biosynthesis includes the enzymatic conversion of ribulose-5-phosphate into 3, 4 dihydroxy-2-butanone-4-phosphates by DHBP synthase. The product of ribAB interferes with both GTP cyclohydrolase II as well as DHBP synthase activities, which catalyze the cleavage of GTP and converts DHBP Ribu5P in the initial steps of both riboflavin biosynthesis branches. Riboswitches are located in the 5' untranslated region (5' UTR) of messenger RNAs and contain an aptamer domain (highly conserved in sequence) where metabolite binding leads to a conformational change in an aptamer domain, which modulate the regulation of gene expression located on bacterial mRNA. In this review, we focus on how riboswitch regulates the riboflavin biosynthesis pathway in Bacillus subtilis and Lactobacillus plantarum.

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.

Developing methods to study conformational changes in RNA crystals using a photocaged ligand

Crystallographic observation of structural changes in real time requires that those changes be uniform both spatially and temporally. A primary challenge with time-resolved ligand-mixing diffraction experiments is asynchrony caused by variable factors, such as efficiency of mixing, rate of diffusion, crystal size, and subsequently, conformational heterogeneity. One method of minimizing such variability is use of a photolabile caged ligand, which can fully saturate the crystal environment (spatially), and whose photoactivation can rapidly (temporally) trigger the reaction in a controlled manner. Our recently published results on a ligand-mixing experiment using time-resolved X-ray crystallography (TRX) with an X-ray free electron laser (XFEL) demonstrated that large conformational changes upon ligand binding resulted in a solid-to-solid phase transition (SSPT), while maintaining Bragg diffraction. Here we investigate this SSPT by polarized video microscopy (PVM) after light-triggered release of a photo-caged adenine (pcADE). In general, the mean transition times and transition widths of the SSPT were less dependent on crystal size than what was observed in previous PVM studies with direct ADE mixing. Instead, the photo-induced transition appears to be heavily influenced by the equilibrium between caged and uncaged ADE due to relatively low sample exposure and uncaging efficiency. Nevertheless, we successfully demonstrate a method for the characterization of phase transitions in RNA crystals that are inducible with a photocaged ligand. The transition data for three crystals of different sizes were then applied to kinetic analysis by fitting to the known four-state model associated with ligand-induced conformational changes, revealing an apparent concentration of uncaged ADE in crystal of 0.43–0.46 mM. These results provide further insight into approaches to study time-resolved ligand-induced conformational changes in crystals, and in particular, highlight the feasibility of triggering phase transitions using a light-inducible system. Developing such approaches may be paramount for the rapidly emerging field of time-resolved crystallography.

NupR Responding to Multiple Signals Is a Nucleoside Permease Regulator in Bacillus thuringiensis BMB171

Nucleoside transport is essential for maintaining intracellular nucleoside and nucleobase homeostasis for living cells. Here, we identified an uncharacterized GntR/HutC family transcriptional regulator, NagR2, renamed NupR (nucleoside permease regulator), that mainly controls nucleoside transport in the Bacillus thuringiensis BMB171 strain. The deletion or overexpression of nupR affected the bacteria's utilization of guanosine, adenosine, uridine, and cytidine rather than thymidine. We further demonstrated that zinc ion is an effector for the NupR, dissociating NupR from its target DNA. Moreover, the expression of nupR is inhibited by NupR, ComK, and PurR, while it is promoted by CcpA. Also, a purine riboswitch located in its 5′ noncoding region influences the expression of nupR. Guanine is the ligand of the riboswitch, reducing the expression of nupR by terminating the transcription of nupR in advance. Hence, our results reveal an exquisite regulation mechanism enabling NupR to respond to multiple signals, control genes involved in nucleoside transport, and contribute to nucleoside substance utilization. Overall, this study provides essential clues for future studies exploring the function of the NupR homolog in other bacteria, such as Bacillus cereus, Bacillus anthracis, Klebsiella pneumoniae, and Streptococcus pneumoniae.

IMPORTANCE The transport of nucleosides and their homeostasis within the cell are essential for growth and proliferation. Here, we have identified a novel transcription factor, NupR, which, to our knowledge, is the first GntR family transcription factor primarily involved in the regulation of nucleoside transport. Moreover, responding to diverse intracellular signals, NupR regulates nucleoside transport. It is vital for utilizing extracellular nucleosides and maintaining intracellular nucleoside homeostasis. NupR may also be involved in other pathways such as pH homeostasis, molybdenum cofactor biosynthesis, nitrate metabolism, and transport. In addition, nucleosides have various applications, such as antiviral drugs. Thus, the elucidation of the transport mechanism of nucleosides could be helpful for the construction of engineered strains for nucleoside production.

Affinity-Based Profiling of the Flavin Mononucleotide Riboswitch

Riboswitches are structural RNA elements that control gene expression. These naturally occurring RNA sensors are of continued interest as antibiotic targets, molecular sensors, and functional elements of synthetic circuits. Here, we describe affinity-based profiling of the flavin mononucleotide (FMN) riboswitch to characterize ligand binding and structural folding. We designed and synthesized photoreactive ligands and used them for photoaffinity labeling. We showed selective labeling of the FMN riboswitch and used this covalent interaction to quantitatively measure ligand binding, which we demonstrate with the naturally occurring antibiotic roseoflavin. We measured conditional riboswitch folding as a function of temperature and cation concentration. Furthermore, combining photoaffinity labeling with reverse transcription revealed ligand binding sites within the aptamer domain with single-nucleotide resolution. The photoaffinity probe was applied to cellular extracts of Bacillus subtilis to demonstrate conditional folding of the endogenous low-abundant ribD FMN riboswitch in biologically derived samples using quantitative PCR. Lastly, binding of the riboswitch-targeting antibiotic roseoflavin to the FMN riboswitch was measured in live bacteria using the photoaffinity probe.

Variants of the guanine riboswitch class exhibit altered ligand specificities for xanthine, guanine, or 2'-deoxyguanosine

The aptamer portions of previously reported riboswitch classes that sense guanine, adenine, or 2′-deoxyguanosine are formed by a highly similar three-stem junction with distinct nucleotide sequences in the regions joining the stems. The nucleotides in these joining regions form the major features of the selective ligand-binding pocket for each aptamer. Previously, we reported the existence of additional, rare variants of the predominant guanine-sensing riboswitch class that carry nucleotide differences in the ligand-binding pocket, suggesting that these RNAs have further diversified their structures and functions. Herein, we report the discovery and analysis of three naturally occurring variants of guanine riboswitches that are narrowly distributed across Firmicutes. These RNAs were identified using comparative sequence analysis methods, which also revealed that some of the gene associations for these variants are atypical for guanine riboswitches or their previously known natural variants. Binding assays demonstrate that the newfound variant riboswitch representatives recognize xanthine, guanine, or 2′-deoxyguanosine, with the guanine class exhibiting greater discrimination against related purines than the more common guanine riboswitch class reported previously. These three additional variant classes, together with the four previously discovered riboswitch classes that employ the same three-stem junction architecture, reveal how a simple structural framework can be diversified to expand the range of purine-based ligands sensed by RNA.

The Second Class of Tetrahydrofolate (THF-II) Riboswitches Recognizes the Tetrahydrofolic Acid Ligand via Local Conformation Changes

Riboswitches are regulatory noncoding RNAs found in bacteria, fungi and plants, that modulate gene expressions through structural changes in response to ligand binding. Understanding how ligands interact with riboswitches in solution can shed light on the molecular mechanisms of this ancient regulators. Previous studies showed that riboswitches undergo global conformation changes in response to ligand binding to relay information. Here, we report conformation switching models of the recently discovered tetrahydrofolic acid-responsive second class of tetrahydrofolate (THF-II) riboswitches in response to ligand binding. Using a combination of selective 2′-hydroxyl acylation, analyzed by primer extension (SHAPE) assay, 3D modeling and small-angle X-ray scattering (SAXS), we found that the ligand specifically recognizes and reshapes the THF-II riboswitch loop regions, but does not affect the stability of the P3 helix. Our results show that the THF-II riboswitch undergoes only local conformation changes in response to ligand binding, rearranging the Loop1-P3-Loop2 region and rotating Loop1 from a ~120° angle to a ~75° angle. This distinct conformation changes suggest a unique regulatory mechanism of the THF-II riboswitch, previously unseen in other riboswitches. Our findings may contribute to the fields of RNA sensors and drug design.

tRNA Fusion to Streamline RNA Structure Determination: Case Studies in Probing Aminoacyl-tRNA Sensing Mechanisms by the T-Box Riboswitch

RNAs are prone to misfolding and are often more challenging to crystallize and phase than proteins. Here, we demonstrate that tRNA fusion can streamline the crystallization and structure determination of target RNA molecules. This strategy was applied to the T-box riboswitch system to capture a dynamic interaction between the tRNA 3′-UCCA tail and the T-box antiterminator, which senses aminoacylation. We fused the T-box antiterminator domain to the tRNA anticodon arm to capture the intended interaction through crystal packing. This approach drastically improved the probability of crystallization and successful phasing. Multiple structure snapshots captured the antiterminator loop in an open conformation with some resemblance to that observed in the recent co-crystal structures of the full-length T box riboswitch–tRNA complex, which contrasts the resting, closed conformation antiterminator observed in an earlier NMR study. The anticipated tRNA acceptor–antiterminator interaction was captured in a low-resolution crystal structure. These structures combined with our previous success using prohead RNA–tRNA fusions demonstrates tRNA fusion is a powerful method in RNA structure determination.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.