Unprecedented tunability of riboswitch structure and regulatory function by sub-millimolar variations in physiological Mg2

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.

Bacillus subtilis YisK possesses oxaloacetate decarboxylase activity and exhibits Mbl-dependent localization

YisK is an uncharacterized protein in Bacillus subtilis previously shown to interact genetically with the elongasome protein Mbl. YisK overexpression leads to cell widening and lysis, phenotypes that are dependent on mbl and suppressed by mbl mutations. In the present work, we characterize YisK's localization, structure, and enzymatic activity. We show that YisK localizes as puncta that depend on Mbl. YisK belongs to the fumarylacetoacetate hydrolase (FAH) superfamily, and crystal structures revealed close structural similarity to two oxaloacetate (OAA) decarboxylases: human mitochondrial FAHD1 and Corynebacterium glutamicum Cg1458. We demonstrate that YisK can also catalyze the decarboxylation of OAA (K m = 134 µM, K cat = 31 min-1). A catalytic dead variant (YisK E148A, E150A) retains wild-type localization and still widens cells following overexpression, indicating these activities are not dependent on YisK catalysis. Conversely, a non-localizing variant (YisK E30A) retains wild-type enzymatic activity in vitro but localizes diffusely and no longer widens cells following overexpression. Together, these results suggest that YisK may be subject to spatial regulation that depends on the cell envelope synthesis machinery. IMPORTANCE The elongasome is a multiprotein complex that guides lengthwise growth in some bacteria. We previously showed that, in B. subtilis, overexpression of an uncharacterized putative enzyme (YisK) perturbed function of the actin-like elongasome protein Mbl. Here, we show that YisK exhibits Mbl-dependent localization. Through biochemical and structural characterization, we demonstrate that, like its mitochondrial homolog FAHD1, YisK can catalyze the decarboxylation of the oxaloacetate to pyruvate and CO2. YisK is the first example of an enzyme implicated in central carbon metabolism with subcellular localization that depends on Mbl.

Analysis of aptamer-target binding and molecular mechanisms by thermofluorimetric analysis and molecular dynamics simulation

Introduction: Aptamers are valuable for bioassays, but aptamer-target binding is susceptible to reaction conditions. In this study, we combined thermofluorimetric analysis (TFA) and molecular dynamics (MD) simulations to optimize aptamer-target binding, explore underlying mechanisms and select preferred aptamer. Methods: Alpha-fetoprotein (AFP) aptamer AP273 (as the model) was incubated with AFP under various experimental conditions, and melting curves were measured in a real-time PCR system to select the optimal binding conditions. The intermolecular interactions of AP273-AFP were analysed by MD simulations with these conditions to reveal the underlying mechanisms. A comparative study between AP273 and control aptamer AP-L3-4 was performed to validate the value of combined TFA and MD simulation in selecting preferred aptamers. Results: The optimal aptamer concentration and buffer system were easily determined from the dF/dT peak characteristics and the melting temperature (Tm) values on the melting curves of related TFA experiments, respectively. A high Tm value was found in TFA experiments performed in buffer systems with low metal ion strength. The molecular docking and MD simulation analyses revealed the underlying mechanisms of the TFA results, i.e., the binding force and stability of AP273 to AFP were affected by the number of binding sites, frequency and distance of hydrogen bonds, and binding free energies; these factors varied in different buffer and metal ion conditions. The comparative study showed that AP273 was superior to the homologous aptamer AP-L3-4. Conclusion: Combining TFA and MD simulation is efficient for optimizing the reaction conditions, exploring underlying mechanisms, and selecting aptamers in aptamer-target bioassays.

Observation of structural switch in nascent SAM-VI riboswitch during transcription at single-nucleotide and single-molecule resolution

Growing RNAs fold differently as they are transcribed, which modulates their finally adopted structures. Riboswitches regulate gene expression by structural change, which are sensitive to co-transcriptionally structural biology. Here we develop a strategy to track the structural change of RNAs during transcription at single-nucleotide and single-molecule resolution and use it to monitor individual transcripts of the SAM-VI riboswitch (riboSAM) as transcription proceeds, observing co-existence of five states in riboSAM. We report a bifurcated helix in one newly identified state from NMR and single-molecule FRET (smFRET) results, and its presence directs the translation inhibition in our cellular translation experiments. A model is proposed to illustrate the distinct switch patterns and gene-regulatory outcome of riboSAM when SAM is present or absent. Our strategy enables the precise mapping of RNAs' conformational landscape during transcription, and may combine with detection methods other than smFRET for structural studies of RNAs in general.

Magnesium Modulates Bacillus subtilis Cell Division Frequency

By chance, we discovered a window of extracellular magnesium (Mg²⁺) availability that modulates the division frequency of Bacillus subtilis without affecting its growth rate. In this window, cells grown with excess Mg²⁺ produce shorter cells than do those grown in unsupplemented medium. The Mg²⁺-responsive adjustment in cell length occurs in both rich and minimal media as well as in domesticated and undomesticated strains. Of other divalent cations tested, manganese (Mn²⁺) and zinc (Zn²⁺) also resulted in cell shortening, but this occurred only at concentrations that affected growth. Cell length decreased proportionally with increasing Mg²⁺ from 0.2 mM to 4.0 mM, with little or no detectable change being observed in labile, intracellular Mg²⁺, based on a riboswitch reporter. Cells grown in excess Mg²⁺ had fewer nucleoids and possessed more FtsZ-rings per unit cell length, consistent with the increased division frequency. Remarkably, when shifting cells from unsupplemented to supplemented medium, more than half of the cell length decrease occurred in the first 10 min, consistent with rapid division onset. Relative to unsupplemented cells, cells growing at steady-state with excess Mg²⁺ showed an enhanced expression of a large number of SigB-regulated genes and the activation of the Fur, MntR, and Zur regulons. Thus, by manipulating the availability of one nutrient, we were able to uncouple the growth rate from the division frequency and identify transcriptional changes that suggest that cell division is accompanied by the general stress response and an enhanced demand to sequester and/or increase the uptake of iron, Mn²⁺, and Zn²⁺. IMPORTANCE The signals that cells use to trigger cell division are unknown. Although division is often considered intrinsic to the cell cycle, microorganisms can continue to grow and repeat rounds of DNA replication without dividing, indicating that cycles of division can be skipped. Here, we show that by manipulating a single nutrient, namely, Mg²⁺, cell division can be uncoupled from the growth rate. This finding can be applied to investigate the nature of the cell division signal(s).

Sarcin/Ricin Domain RNA Retains Its Structure Better Than A-RNA in Adaptively Biased Molecular Dynamics Simulations

Less than one in thirty of the RNA sequences transcribed in humans are translated into protein. The noncoding RNA (ncRNA) functions in catalysis, structure, regulation, and more. However, for the most part, these functions are poorly characterized. RNA is modular and described by motifs that include helical A-RNA with canonical Watson-Crick base-pairing as well as structures with only noncanonical base pairs. Understanding the structure and dynamics of motifs will aid in deciphering functions of specific ncRNAs. We present computational studies on a standard sarcin/ricin domain (SRD), citrus bark cracking viroid SRD, as well as A-RNA. We have applied enhanced molecular dynamics techniques that construct an inverse free-energy surface (iFES) determined by collective variables that monitor base-pairing and backbone conformation. Each SRD RNA is flanked on each side by A-RNA, allowing comparison of the behavior of these motifs in the same molecule. The RNA iFESs have single peaks, indicating that the combined motifs should denature as a single cohesive unit, rather than by regional melting. Local root-mean-square deviation (RMSD) analysis and communication propensity (CProp, variance in distances between residue pairs) reveal distinct motif properties. Our analysis indicates that the standard SRD is more stable than the viroid SRD, which is more stable than A-RNA. Base pairs at SRD to A-RNA transitions have limited flexibility. Application of CProp reveals extraordinary stiffness of the SRD, allowing residues on opposite sides of the motif to sense each other's motions.

Synergism in the Molecular Crowding of Ligand-Induced Riboswitch Folding: Kinetic/Thermodynamic Insights from Single-Molecule Spectroscopy

Conformational dynamics in riboswitches involves ligand binding and folding of RNA, each of which can be influenced by excluded volume effects under "crowded" in vivo cellular conditions and thus incompletely characterized by in vitro studies under dilute buffer conditions. In this work, temperature-dependent single-molecule fluorescence resonance energy transfer (FRET) spectroscopy is used to characterize the thermodynamics of (i) cognate ligand and (ii) molecular crowders (PEG, polyethylene glycol) on folding of the B. subtilis LysC lysine riboswitch. With the help of detailed kinetic analysis, we isolate and study the effects of PEG on lysine binding and riboswitch folding steps individually, from which we find that PEG crowding facilitates riboswitch folding primarily via a surprising increase in affinity for the cognate ligand. This is furthermore confirmed by temperature-dependent studies, which reveal that PEG crowding is not purely entropic and instead significantly impacts both enthalpic and entropic contributions to the free energy landscape for folding. The results indicate that PEG molecular crowding/stabilization of the lysine riboswitch is more mechanistically complex and requires extension beyond the conventional picture of purely repulsive solvent-solute steric interactions arising from excluded volume and entropy. Instead, the current experimental FRET data support an alternative multistep mechanism, whereby PEG first entropically crowds the unfolded riboswitch into a "pre-folded" conformation, which in turn greatly increases the ligand binding affinity and thereby enhances the overall equilibrium for riboswitch folding.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

Cellular Adaptations to Cytoplasmic Mg2+ Limitation

Mg²⁺ is the most abundant divalent cation in living cells. It is essential for charge neutralization, macromolecule stabilization, and the assembly and activity of ribosomes and as a cofactor for enzymatic reactions. When experiencing low cytoplasmic Mg²⁺, bacteria adopt two main strategies: They increase the abundance and activity of Mg²⁺ importers and decrease the abundance of Mg²⁺-chelating ATP and rRNA. These changes reduce regulated proteolysis by ATP-dependent proteases and protein synthesis in a systemic fashion. In many bacterial species, the transcriptional regulator PhoP controls expression of proteins mediating these changes. The 5' leader region of some mRNAs responds to low cytoplasmic Mg²⁺ or to disruptions in translation of open reading frames in the leader regions by furthering expression of the associated coding regions, which specify proteins mediating survival when the cytoplasmic Mg²⁺ concentration is low. Microbial species often utilize similar adaptation strategies to cope with low cytoplasmic Mg²⁺ despite relying on different genes to do so.

Riboswitch Mechanisms: New Tricks for an Old Dog

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

Dual-mode aptasensor for simultaneous detection of multiple food-borne pathogenic bacteria based on colorimetry and microfluidic chip using stir bar sorptive extraction

A dual-mode aptasensor using colorimetry and microfluidic chip (MC) together with stir bar sorptive extraction (SBSE) has been developed for firstly qualifying samples contaminated with Vibrio parahaemolyticus (V.P) and Salmonella typhimurium (S.T), then precisely determine both of them in positive samples. For this purpose, the aptamer-streptavidin encoded probes (Apt-SAEs) corresponding to different bacteria were prepared in advance. Then, a stir bar modified with 4-mercaptophenylboronic acid (MPBA) was made to extract bacteria together with Apt-SAE probes. The binding event of aptamer and target triggered the formation of two sandwich structures containing Apt-SAE, V.P or S.T. The concentration of bacteria could be enriched by 1000 times within 15 min to avoid long-time enrichment process. Finally, the stir bar was immersed in the 3,3',5,5'-Tetramethylbenzidine (TMB)-H₂O₂ solution for color development. The color could be observed by naked eyes to judge whether the analytes were present. The colorless samples were judged to be negative. For the positive samples, the adsorbed encoded probes corresponding to different bacteria would be eluted from the stir bar and rapidly analyzed by the MC. Under the optimized conditions, 100 CFU/mL of V.P or S.T or both of them could be observed by colorimetry and 35 CFU/mL of them could be detected (S/N = 3) by the MC. The assay has significant application value for on-site screening and multiple detection of food-borne pathogenic bacteria.

Microsecond Folding of preQ1 Riboswitch and Its Biological Significance Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

Riboswitches are regulatory elements of bacterial mRNA which function with conformational switching upon binding of specific cellular metabolites. In particular, transcriptional riboswitches regulate gene expression kinetically through the conformational change of the aptamer domain. In this study, we investigate the conformational dynamics and ligand binding mechanisms of the aptamer domain of a transcriptional prequeuosine (preQ₁) riboswitch from Bacillus subtilis using two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) with microsecond time resolution. The obtained time-resolved single-molecule data indicate that the aptamer domain undergoes folding/unfolding including three forms, which are attributed to hairpin (O), pseudoknot-like (pF), and H-type pseudoknot (fF) structures. It is found that a cofactor, Mg²⁺, binds only to the fF form with the conformational selection mechanism. In contrast, it is indicated that the ligand, preQ₁, binds to the O form with the induced-fit mechanism and significantly accelerates the microsecond O → pF folding process. It is also shown that the binding with preQ₁ substantially stabilizes the fF form that is generated from the pF form with a long time constant (>10 ms). Combining these results with the results of a former smFRET study on the slower time scale, we obtain an overall picture of the folding/unfolding dynamics of the aptamer domain as well as its energy landscape. On the basis of the picture obtained, we discuss the significance of the microsecond folding/unfolding of the aptamer domain for biological function of the riboswitch and propose the molecular mechanism of the gene expression controlled by the structural dynamics of the aptamer domain.

Thermal adaptation of structural dynamics and regulatory function of adenine riboswitch

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

Transcriptional Riboswitches Integrate Timescales for Bacterial Gene Expression Control

Transcriptional riboswitches involve RNA aptamers that are typically found in the 5' untranslated regions (UTRs) of bacterial mRNAs and form alternative secondary structures upon binding to cognate ligands. Alteration of the riboswitch's secondary structure results in perturbations of an adjacent expression platform that controls transcription elongation and termination, thus turning downstream gene expression "on" or "off." Riboswitch ligands are typically small metabolites, divalent cations, anions, signaling molecules, or other RNAs, and can be part of larger signaling cascades. The interconnectedness of ligand binding, RNA folding, RNA transcription, and gene expression empowers riboswitches to integrate cellular processes and environmental conditions across multiple timescales. For a successful response to an environmental cue that may determine a bacterium's chance of survival, a coordinated coupling of timescales from microseconds to minutes must be achieved. This review focuses on recent advances in our understanding of how riboswitches affect such critical gene expression control across time.

An NMR-based approach reveals the core structure of the functional domain of SINEUP lncRNAs

Long non-coding RNAs (lncRNAs) are attracting widespread attention for their emerging regulatory, transcriptional, epigenetic, structural and various other functions. Comprehensive transcriptome analysis has revealed that retrotransposon elements (REs) are transcribed and enriched in lncRNA sequences. However, the functions of lncRNAs and the molecular roles of the embedded REs are largely unknown. The secondary and tertiary structures of lncRNAs and their embedded REs are likely to have essential functional roles, but experimental determination and reliable computational prediction of large RNA structures have been extremely challenging. We report here the nuclear magnetic resonance (NMR)-based secondary structure determination of the 167-nt inverted short interspersed nuclear element (SINE) B2, which is embedded in antisense Uchl1 lncRNA and upregulates the translation of sense Uchl1 mRNAs. By using NMR 'fingerprints' as a sensitive probe in the domain survey, we successfully divided the full-length inverted SINE B2 into minimal units made of two discrete structured domains and one dynamic domain without altering their original structures after careful boundary adjustments. This approach allowed us to identify a structured domain in nucleotides 31-119 of the inverted SINE B2. This approach will be applicable to determining the structures of other regulatory lncRNAs.

High pressure single-molecule FRET studies of the lysine riboswitch: cationic and osmolytic effects on pressure induced denaturation

Deep sea biology is known to thrive at pressures up to ≈1 kbar, which motivates fundamental biophysical studies of biomolecules under such extreme environments. In this work, the conformational equilibrium of the lysine riboswitch has been systematically investigated by single molecule FRET (smFRET) microscopy at pressures up to 1500 bar. The lysine riboswitch preferentially unfolds with increasing pressure, which signals an increase in free volume (ΔV0 > 0) upon folding of the biopolymer. Indeed, the effective lysine binding constant increases quasi-exponentially with pressure rise, which implies a significant weakening of the riboswitch-ligand interaction in a high-pressure environment. The effects of monovalent/divalent cations and osmolytes on folding are also explored to acquire additional insights into cellular mechanisms for adapting to high pressures. For example, we find that although Mg2+ greatly stabilizes folding of the lysine riboswitch (ΔΔG0 < 0), there is negligible impact on changes in free volume (ΔΔV0 ≈ 0) and thus any pressure induced denaturation effects. Conversely, osmolytes (commonly at high concentrations in deep sea marine species) such as the trimethylamine N-oxide (TMAO) significantly reduce free volumes (ΔΔV0 < 0) and thereby diminish pressure-induced denaturation. We speculate that, besides stabilizing RNA structure, enhanced levels of TMAO in cells might increase the dynamic range for competent riboswitch folding by suppressing the pressure-induced denaturation response. This in turn could offer biological advantage for vertical migration of deep-sea species, with impacts on food searching in a resource limited environment.