Structure and mechanism of a metal-sensing regulatory RNA

Calculating metalation in cells reveals CobW acquires CoII for vitamin B12 biosynthesis while related proteins prefer ZnII

Protein metal-occupancy (metalation) in vivo has been elusive. To address this challenge, the available free energies of metals have recently been determined from the responses of metal sensors. Here, we use these free energy values to develop a metalation-calculator which accounts for inter-metal competition and changing metal-availabilities inside cells. We use the calculator to understand the function and mechanism of GTPase CobW, a predicted Co^II-chaperone for vitamin B₁₂. Upon binding nucleotide (GTP) and Mg^II, CobW assembles a high-affinity site that can obtain Co^II or Zn^II from the intracellular milieu. In idealised cells with sensors at the mid-points of their responses, competition within the cytosol enables Co^II to outcompete Zn^II for binding CobW. Thus, Co^II is the cognate metal. However, after growth in different [Co^II], Co^II-occupancy ranges from 10 to 97% which matches CobW-dependent B₁₂ synthesis. The calculator also reveals that related GTPases with comparable Zn^II affinities to CobW, preferentially acquire Zn^II due to their relatively weaker Co^II affinities. The calculator is made available here for use with other proteins.

Current Understandings on Magnesium Deficiency and Future Outlooks for Sustainable Agriculture

The productivity of agricultural produce is fairly dependent on the availability of nutrients and efficient use. Magnesium (Mg2+) is an essential macronutrient of living cells and is the second most prevalent free divalent cation in plants. Mg2+ plays a role in several physiological processes that support plant growth and development. However, it has been largely forgotten in fertilization management strategies to increase crop production, which leads to severe reductions in plant growth and yield. In this review, we discuss how the Mg2+ shortage induces several responses in plants at different levels: morphological, physiological, biochemical and molecular. Additionally, the Mg2+ uptake and transport mechanisms in different cellular organelles and the role of Mg2+ transporters in regulating Mg2+ homeostasis are also discussed. Overall, in this review, we critically summarize the available information about the responses of Mg deficiency on plant growth and development, which would facilitate plant scientists to create Mg2+-deficiency-resilient crops through agronomic and genetic biofortification.

Biocompatible and functional inorganic magnesium ceramic particles for biomedical applications

Magnesium ceramics hold promise for numerous biological applications. This review covers the synthesis of magnesium ceramic particles with specific morphologies and potential modification techniques. Magnesium ceramic particles possess multiple characteristics directly applicable to human biology; they are anti-inflammatory, antibacterial, antiviral, and offer anti-cancer effects. Based on these advantages, magnesium hydroxide nanoparticles have been extensively utilized across biomedical fields. In a vascular stent, the incorporation of magnesium ceramic nanoparticles enhances re-endothelialization. Additionally, tissue regeneration for bone, cartilage, and kidney can be promoted by magnesium ceramics. This review enables researchers to identify the optimum synthetic conditions to prepare magnesium ceramics with specific morphologies and sizes and select the appropriate modification protocols. It is also intended to elucidate the desirable physicochemical properties and biological benefits of magnesium ceramics.

Cooperativity and Allostery in RNA Systems

Allostery is among the most basic biological principles employed by biological macromolecules to achieve a biologically active state in response to chemical cues. Although initially used to describe the impact of small molecules on the conformation and activity of protein enzymes, the definition of this term has been significantly broadened to describe long-range conformational change of macromolecules in response to small or large effectors. Such a broad definition could be applied to RNA molecules, which do not typically serve as protein-free cellular enzymes but fold and form macromolecular assemblies with the help of various ligand molecules, including ions and proteins. Ligand-induced allosteric changes in RNA molecules are often accompanied by cooperative interactions between RNA and its ligand, thus streamlining the folding and assembly pathways. This chapter provides an overview of the interplay between cooperativity and allostery in RNA systems and outlines methods to study these two biological principles.

RNA structure prediction using positive and negative evolutionary information

Knowing the structure of conserved structural RNAs is important to elucidate their function and mechanism of action. However, predicting a conserved RNA structure remains unreliable, even when using a combination of thermodynamic stability and evolutionary covariation information. Here we present a method to predict a conserved RNA structure that combines the following three features. First, it uses significant covariation due to RNA structure and removes spurious covariation due to phylogeny. Second, it uses negative evolutionary information: basepairs that have variation but no significant covariation are prevented from occurring. Lastly, it uses a battery of probabilistic folding algorithms that incorporate all positive covariation into one structure. The method, named CaCoFold (Cascade variation/covariation Constrained Folding algorithm), predicts a nested structure guided by a maximal subset of positive basepairs, and recursively incorporates all remaining positive basepairs into alternative helices. The alternative helices can be compatible with the nested structure such as pseudoknots, or overlapping such as competing structures, base triplets, or other 3D non-antiparallel interactions. We present evidence that CaCoFold predictions are consistent with structures modeled from crystallography.

The availability of deeper comparative sequence alignments and recent advances in statistical analysis of RNA sequence covariation have made it possible to identify a reliable set of conserved base pairs, as well as a reliable set of non-basepairs (positions that vary without covarying). Predicting an overall consensus secondary structure consistent with a set of individual inferred pairs and non-pairs remains a problem. Current RNA structure prediction algorithms that predict nested secondary structures cannot use the full set of inferred covarying pairs, because covariation analysis also identifies important non-nested pairing interactions such as pseudoknots, base triples, and alternative structures. Moreover, although algorithms for incorporating negative constraints exist, negative information from covariation analysis (inferred non-pairs) has not been systematically exploited. Here I introduce an efficient approximate RNA structure prediction algorithm that incorporates all inferred pairs and excludes all non-pairs. Using this, and an improved visualization tool, I show that the method correctly identifies many non-nested structures in agreement with known crystal structures, and improves many curated consensus secondary structure annotations in RNA sequence alignment databases.

Mg2+ Sensing by an RNA Fragment: Role of Mg2+-Coordinated Water Molecules

RNA molecules selectively bind to specific metal ions to populate their functional active states, making it important to understand their source of ion selectivity. In large RNA systems, metal ions interact with the RNA at multiple locations, making it difficult to decipher the precise role of ions in folding. To overcome this complexity, we studied the role of different metal ions (Mg²⁺, Ca²⁺, and K⁺) in the folding of a small RNA hairpin motif (5'-ucCAAAga-3') using unbiased all-atom molecular dynamics simulations. The advantage of studying this system is that it requires specific binding of a single metal ion to fold to its native state. We find that even for this small RNA, the folding free energy surface (FES) is multidimensional as different metal ions present in the solution can simultaneously facilitate folding. The FES shows that specific binding of a metal ion is indispensable for its folding. We further show that in addition to the negatively charged phosphate groups, the spatial organization of electronegative nucleobase atoms drives the site-specific binding of the metal ions. Even though the binding site cannot discriminate between different metal ions, RNA folds efficiently only in a Mg²⁺ solution. We show that the rigid network of Mg²⁺-coordinated water molecules facilitates the formation of important interactions in the transition state. The other metal ions such as K⁺ and Ca²⁺ cannot facilitate the formation of such interactions. These results allow us to hypothesize possible metal-sensing mechanisms in large metalloriboswitches and also provide useful insights into the design of appropriate collective variables for studying large RNA molecules using enhanced sampling methods.

19 F NMR-Based Fragment Screening for 14 Different Biologically Active RNAs and 10 DNA and Protein Counter-Screens

We report here the nuclear magnetic resonance ¹⁹F screening of 14 RNA targets with different secondary and tertiary structure to systematically assess the druggability of RNAs. Our RNA targets include representative bacterial riboswitches that naturally bind with nanomolar affinity and high specificity to cellular metabolites of low molecular weight. Based on counter‐screens against five DNAs and five proteins, we can show that RNA can be specifically targeted. To demonstrate the quality of the initial fragment library that has been designed for easy follow‐up chemistry, we further show how to increase binding affinity from an initial fragment hit by chemistry that links the identified fragment to the intercalator acridine. Thus, we achieve low‐micromolar binding affinity without losing binding specificity between two different terminator structures.

Crystal structure of an atypical cobalamin riboswitch reveals RNA structural adaptability as basis for promiscuous ligand binding

Cobalamin riboswitches encompass a structurally diverse group of cis-acting, gene regulatory elements found mostly in bacterial messenger RNA and are classified into subtypes based on secondary and tertiary characteristics. An unusual variant of the cobalamin riboswitch with predicted structural features was identified in Bacillus subtilis over a decade ago, but its structure and mechanisms of cobalamin selectivity and translational control have remained unsolved. We present the crystal structure of the aptamer domain of this atypical cobalamin riboswitch and a model for the complete riboswitch, including its expression platform domain. We demonstrate that this riboswitch binds to multiple cobalamin derivatives and correlate its promiscuous behavior to its structure and unique arrangement of peripheral elements. Comparative structural analyses between conventional cobalamin riboswitches and the B. subtilis cobalamin riboswitch reveal that the likely basis for this promiscuous ligand binding is intrinsic structural adaptability encoded in the RNA structure. It suggests that cobalamin selectivity might ultimately be viewed as existing on a spectrum of affinity for each derivative rather than as belonging to distinct types based on ligand specificities. Our work provides an interesting and notable example of functional coupling of ligand-sensing and adaptive folding by a structured RNA molecule.

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

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

Former orphan riboswitches reveal unexplored areas of bacterial metabolism, signaling, and gene control processes

Comparative sequence analyses have been used to discover numerous classes of structured noncoding RNAs, some of which are riboswitches that specifically recognize small-molecule or elemental ion ligands and influence expression of adjacent downstream genes. Determining the correct identity of the ligand for a riboswitch candidate typically is aided by an understanding of the genes under its regulatory control. Riboswitches whose ligands were straightforward to identify have largely been associated with well-characterized metabolic pathways, such as coenzyme or amino acid biosynthesis. Riboswitch candidates whose ligands resist identification, collectively known as orphan riboswitches, are often associated with genes coding for proteins of unknown function, or genes for various proteins with no established link to one another. The cognate ligands for 16 former orphan riboswitch motifs have been identified to date. The successful pursuit of the ligands for these classes has provided insight into areas of biology that are not yet fully explored, such as ion homeostasis, signaling networks, and other previously underappreciated biochemical or physiological processes. Herein we discuss the strategies and methods used to match ligands with orphan riboswitch classes, and overview the lessons learned to inform and motivate ongoing efforts to identify ligands for the many remaining candidates.

Discovery of ANTAR-RNAs and their Mechanism of Action in Mycobacteria

Non-coding RNAs play pivotal roles in bacterial signaling. However, RNAs from certain phyla (specially high-GC actinobacteria) still remain elusive. Here, by re-engineering the existing genome-wide search approach, we discover a family of structurally conserved RNAs that are present ubiquitously across actinobacteria, including mycobacteria. In vitro analysis shows that RNAs belonging to this family bind response-regulator proteins that contain the widely prevalent ANTAR domain. The Mycobacterium tuberculosis ANTAR protein gets phosphorylated by a histidine kinase and interacts with RNA only in its phosphorylated state. These newly identified RNAs reside only in certain transcripts and typically overlap with the ribosome-binding site, regulating translation of these transcripts. In this way, the RNAs directly link signaling pathways to translational control, thus expanding the mechanistic tool kit available for ANTAR-based control of gene expression. In mycobacteria, we find that RNAs targeted by ANTAR proteins majorly encode enzymes of lipid metabolism and associated redox pathways. This now allows us to identify the key genes that mediate ANTAR-dependent control of lipid metabolism. Our study establishes the identity and wide prevalence of ANTAR-target RNAs in mycobacteria, bringing RNA-mediated regulation in these bacteria to the center stage.

The czcD (NiCo) Riboswitch Responds to Iron(II)

Iron is essential for nearly every organism, and mismanagement of its intracellular concentrations (either deficiency or excess) contributes to diminished virulence in human pathogens, necessitating intricate metalloregulatory mechanisms. To date, although several metal-responsive riboswitches have been identified in bacteria, none has been shown to respond to Fe^II. The czcD riboswitch, present in numerous human gut microbiota and pathogens, was recently shown to respond to Ni^II and Co^II but thought not to respond to Fe^II, on the basis of aerobic, in vitro assays; its function in vivo is not well understood. We constructed a fluorescent sensor using this riboswitch fused to the RNA aptamer, Spinach2. When assayed anaerobically, the resulting sensor responds in vitro to Fe^II, as well as to Mn^II, Co^II, Ni^II, and Zn^II, but only in the cases of Fe^II and Mn^II do the apparent K_d values (0.4 and 11 μM, respectively) fall within the range of labile metal concentrations maintained by known metalloregulators. We also show that the sensor-which is, to the best of our knowledge, the first reversible genetically encoded fluorescent sensor for Fe^II-responds to iron in Escherichia coli cells. Finally, we demonstrate that the putative metal exporters directly downstream of two czcD riboswitches efficiently rescue iron toxicity in a heterologous expression system. Together, our results indicate that iron merits consideration as a plausible physiological ligand for czcD riboswitches, although a response to general metal stress cannot be ruled out at present.

Ligand Binding Mechanism and Its Relationship with Conformational Changes in Adenine Riboswitch

Riboswitches are naturally occurring RNA aptamers that control the expression of essential bacterial genes by binding to specific small molecules. The binding with both high affinity and specificity induces conformational changes. Thus, riboswitches were proposed as a possible molecular target for developing antibiotics and chemical tools. The adenine riboswitch can bind not only to purine analogues but also to pyrimidine analogues. Here, long molecular dynamics (MD) simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) computational methodologies were carried out to show the differences in the binding model and the conformational changes upon five ligands binding. The binding free energies of the guanine riboswitch aptamer with C74U mutation complexes were compared to the binding free energies of the adenine riboswitch (AR) aptamer complexes. The calculated results are in agreement with the experimental data. The differences for the same ligand binding to two different aptamers are related to the electrostatic contribution. Binding dynamical analysis suggests a flexible binding pocket for the pyrimidine ligand in comparison with the purine ligand. The 18 μs of MD simulations in total indicate that both ligand-unbound and ligand-bound aptamers transfer their conformation between open and closed states. The ligand binding obviously affects the conformational change. The conformational states of the aptamer are associated with the distance between the mass center of two key nucleotides (U51 and A52) and the mass center of the other two key nucleotides (C74 and C75). The results suggest that the dynamical character of the binding pocket would affect its biofunction. To design new ligands of the adenine riboswitch, it is recommended to consider the binding affinities of the ligand and the conformational change of the ligand binding pocket.

Mg2+ homeostasis and transport in cyanobacteria - at the crossroads of bacterial and chloroplast Mg2+ import

Magnesium cation (Mg2+) is the most abundant divalent cation in living cells, where it is required for various intracellular functions. In chloroplasts and cyanobacteria, established photosynthetic model systems, Mg2+ is the central ion in chlorophylls, and Mg2+ flux across the thylakoid membrane is required for counterbalancing the light-induced generation of a ΔpH across the thylakoid membrane. Yet, not much is known about Mg2+ homoeostasis, transport and distribution within cyanobacteria. However, Mg2+ transport across membranes has been studied in non-photosynthetic bacteria, and first observations and findings are reported for chloroplasts. Cyanobacterial cytoplasmic membranes appear to contain the well-characterized Mg2+ channels CorA and/or MgtE, which both facilitate transmembrane Mg2+ flux down the electrochemical gradient. Both Mg2+ channels are typical for non-photosynthetic bacteria. Furthermore, Mg2+ transporters of the MgtA/B family are also present in the cytoplasmic membrane to mediate active Mg2+ import into the bacterial cell. While the cytoplasmic membrane of cyanobacteria resembles a 'classical' bacterial membrane, essentially nothing is known about Mg2+ channels and/or transporters in thylakoid membranes of cyanobacteria or chloroplasts. As discussed here, at least one Mg2+ channelling protein must be localized within thylakoid membranes. Thus, either one of the 'typical' bacterial Mg2+ channels has a dual localization in the cytoplasmic plus the thylakoid membrane, or another, yet unidentified channel is present in cyanobacterial thylakoid membranes.

Rhomboid intramembrane protease YqgP licenses bacterial membrane protein quality control as adaptor of FtsH AAA protease

Magnesium homeostasis is essential for life and depends on magnesium transporters, whose activity and ion selectivity need to be tightly controlled. Rhomboid intramembrane proteases pervade the prokaryotic kingdom, but their functions are largely elusive. Using proteomics, we find that Bacillus subtilis rhomboid protease YqgP interacts with the membrane‐bound ATP‐dependent processive metalloprotease FtsH and cleaves MgtE, the major high‐affinity magnesium transporter in B. subtilis. MgtE cleavage by YqgP is potentiated in conditions of low magnesium and high manganese or zinc, thereby protecting B. subtilis from Mn²⁺/Zn²⁺ toxicity. The N‐terminal cytosolic domain of YqgP binds Mn²⁺ and Zn²⁺ ions and facilitates MgtE cleavage. Independently of its intrinsic protease activity, YqgP acts as a substrate adaptor for FtsH, a function that is necessary for degradation of MgtE. YqgP thus unites protease and pseudoprotease function, hinting at the evolutionary origin of rhomboid pseudoproteases such as Derlins that are intimately involved in eukaryotic ER‐associated degradation (ERAD). Conceptually, the YqgP‐FtsH system we describe here is analogous to a primordial form of “ERAD” in bacteria and exemplifies an ancestral function of rhomboid‐superfamily proteins.

A Mn-sensing riboswitch activates expression of a Mn2+/Ca2+ ATPase transporter in Streptococcus

Maintaining manganese (Mn) homeostasis is important for the virulence of numerous bacteria. In the human respiratory pathogen Streptococcus pneumoniae, the Mn-specific importer PsaBCA, exporter MntE, and transcriptional regulator PsaR establish Mn homeostasis. In other bacteria, Mn homeostasis is controlled by yybP-ykoY family riboswitches. Here, we characterize a yybP-ykoY family riboswitch upstream of the mgtA gene encoding a P_II-type ATPase in S. pneumoniae, suggested previously to function in Ca²⁺ efflux. We show that the mgtA riboswitch aptamer domain adopts a canonical yybP-ykoY structure containing a three-way junction that is compacted in the presence of Ca²⁺ or Mn²⁺ at a physiological Mg²⁺ concentration. Although Ca²⁺ binds to the RNA aptamer with higher affinity than Mn²⁺, in vitro activation of transcription read-through of mgtA by Mn²⁺ is much greater than by Ca²⁺. Consistent with this result, mgtA mRNA and protein levels increase ≈5-fold during cellular Mn stress, but only in genetic backgrounds of S. pneumoniae and Bacillus subtilis that exhibit Mn²⁺ sensitivity, revealing that this riboswitch functions as a failsafe ‘on’ signal to prevent Mn²⁺ toxicity in the presence of high cellular Mn²⁺. In addition, our results suggest that the S. pneumoniae yybP-ykoY riboswitch functions to regulate Ca²⁺ efflux under these conditions.

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

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

Alternate RNA Structures

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

Concerted dynamics of metallo-base pairs in an A/B-form helical transition

Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-Hg^II-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-Hg^II-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-Hg^II-(N4)C base pairs and greater A-form characteristics. Ten-fold different ¹J coupling constants (¹⁵N,¹⁹⁹Hg) are observed for (N3)C-Hg^II (114 Hz) versus (N4)C-Hg^II (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-Hg^II-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding Hg^II, demonstrating that C-Hg^II-T has unique switching potential in DNA-based materials and devices.

Local-to-global signal transduction at the core of a Mn2+ sensing riboswitch

The widespread Mn2+-sensing yybP-ykoY riboswitch controls the expression of bacterial Mn2+ homeostasis genes. Here, we first determine the crystal structure of the ligand-bound yybP-ykoY riboswitch aptamer from Xanthomonas oryzae at 2.96 Å resolution, revealing two conformations with docked four-way junction (4WJ) and incompletely coordinated metal ions. In >100 µs of MD simulations, we observe that loss of divalents from the core triggers local structural perturbations in the adjacent docking interface, laying the foundation for signal transduction to the regulatory switch helix. Using single-molecule FRET, we unveil a previously unobserved extended 4WJ conformation that samples transient docked states in the presence of Mg2+. Only upon adding sub-millimolar Mn2+, however, can the 4WJ dock stably, a feature lost upon mutation of an adenosine contacting Mn2+ in the core. These observations illuminate how subtly differing ligand preferences of competing metal ions become amplified by the coupling of local with global RNA dynamics.

Genetic screens reveal novel major and minor players in magnesium homeostasis of Staphylococcus aureus

Magnesium is one of the most abundant metal ions in living cells. Very specific and devoted transporters have evolved for transporting Mg2+ ions across the membrane and maintain magnesium homeostasis. Using genetic screens, we were able to identify the main players in magnesium homeostasis in the opportunistic pathogen Staphylococcus aureus. Here, we show that import of magnesium relies on the redundant activity of either CorA2 or MgtE since in absence of these two importers, bacteria require increased amounts of magnesium in the medium. A third CorA-like importer seems to play a minor role, at least under laboratory conditions. For export of magnesium, we identified two proteins, MpfA and MpfB. MpfA, is the main actor since it is essential for growth in high magnesium concentrations. We show that gain of function mutations or overexpression of the minor factor, MpfB, which is part of a sigmaB controlled stress response regulon, can compensate for the absence of MpfA.

Single-Cell Microscopy Reveals That Levels of Cyclic di-GMP Vary among Bacillus subtilis Subpopulations

The synthesis of signaling molecules is one strategy bacteria employ to sense alterations in their environment and rapidly adjust to those changes. In Gram-negative bacteria, bis-(3'-5')-cyclic dimeric GMP (c-di-GMP) regulates the transition from a unicellular motile state to a multicellular sessile state. However, c-di-GMP signaling has been less intensively studied in Gram-positive organisms. To that end, we constructed a fluorescent yfp reporter based on a c-di-GMP-responsive riboswitch to visualize the relative abundance of c-di-GMP for single cells of the Gram-positive model organism Bacillus subtilis Coupled with cell-type-specific fluorescent reporters, this riboswitch reporter revealed that c-di-GMP levels are markedly different among B. subtilis cellular subpopulations. For example, cells that have made the decision to become matrix producers maintain higher intracellular c-di-GMP concentrations than motile cells. Similarly, we find that c-di-GMP levels differ between sporulating and competent cell types. These results suggest that biochemical measurements of c-di-GMP abundance are likely to be inaccurate for a bulk ensemble of B. subtilis cells, as such measurements will average c-di-GMP levels across the population. Moreover, the significant variation in c-di-GMP levels between cell types hints that c-di-GMP might play an important role during B. subtilis biofilm formation. This study therefore emphasizes the importance of using single-cell approaches for analyzing metabolic trends within ensemble bacterial populations.IMPORTANCE Many bacteria have been shown to differentiate into genetically identical yet morphologically distinct cell types. Such population heterogeneity is especially prevalent among biofilms, where multicellular communities are primed for unexpected environmental conditions and can efficiently distribute metabolic responsibilities. Bacillus subtilis is a model system for studying population heterogeneity; however, a role for c-di-GMP in these processes has not been thoroughly investigated. Herein, we introduce a fluorescent reporter, based on a c-di-GMP-responsive riboswitch, to visualize the relative abundance of c-di-GMP for single B. subtilis cells. Our analysis shows that c-di-GMP levels are conspicuously different among B. subtilis cellular subtypes, suggesting a role for c-di-GMP during biofilm formation. These data highlight the utility of riboswitches as tools for imaging metabolic changes within individual bacterial cells. Analyses such as these offer new insight into c-di-GMP-regulated phenotypes, especially given that other biofilms also consist of multicellular communities.

Structural Basis of Transcription: RNA Polymerase Backtracking and Its Reactivation

Regulatory sequences or erroneous incorporations during DNA transcription cause RNA polymerase backtracking and inactivation in all kingdoms of life. Reactivation requires RNA transcript cleavage. Essential transcription factors (GreA and GreB, or TFIIS) accelerate this reaction. We report four cryo-EM reconstructions of Escherichia coli RNA polymerase representing the entire reaction pathway: (1) a backtracked complex; a backtracked complex with GreB (2) before and (3) after RNA cleavage; and (4) a reactivated, substrate-bound complex with GreB before RNA extension. Compared with eukaryotes, the backtracked RNA adopts a different conformation. RNA polymerase conformational changes cause distinct GreB states: a fully engaged GreB before cleavage; a disengaged GreB after cleavage; and a dislodged, loosely bound GreB removed from the active site to allow RNA extension. These reconstructions provide insight into the catalytic mechanism and dynamics of RNA cleavage and extension and suggest how GreB targets backtracked complexes without interfering with canonical transcription.

Riboswitches in Archaea

Background:

Riboswitches are cis-acting, non-coding RNA elements found in the 5’UTR of bacterial mRNA and 3’ UTR of eukaryotic mRNA, that fold in a complex manner to act as receptors for specific metabolites hence altering their conformation in response to the change in concentrations of a ligand or metabolite. Riboswitches function as gene regulators in numerous bacteria, archaea, fungi, algae and plants.

Aim and Objective:

This study identifies different classes of riboswitches in the Archaeal domain of life. Previous studies have suggested that riboswitches carry a conserved aptameric domain in different domains of life. Since Archaea are considered to be the most idiosyncratic organisms it was interesting to look for the conservation pattern of riboswitches in these obviously strange microorganisms.

Materials and Methods:

Completely sequenced Archaeal Genomes present in the NCBI repository were used for studying riboswitches and other ncRNAs. The sequence files in FASTA format were downloaded from NCBI Genome database and information related to these genomes was retrieved from GenBank. Three bioinformatics approaches were used namely, ab initio, consensus structure prediction and statistical model-based prediction for identifying riboswitches.

Results:

Archaeal genomes have a sporadic distribution of putative riboswitches like the TPP, FMN, Guanidine, Lysine and c-di-AMP riboswitches, which are known to occur in bacteria. Also, a class of riboswitch sensing c-di-GMP, a second messenger, has been identified in a few Archaeal organisms.

Conclusion:

This study clearly reveals that bioinformatics methods are likely to play a major role in identifying conserved riboswitches and in establishing how widespread these classes are in all domains of life, even though the final confirmation may come from wet lab methods.

Evaluating riboswitch optimality

Riboswitches are RNA elements that recognize diverse chemical and biomolecular inputs, and transduce this recognition process to genetic, fluorescent, and other engineered outputs using RNA conformational changes. These systems are pervasive in cellular biology and are a promising biotechnology with applications in genetic regulation and biosensing. Here, we derive a simple expression bounding the activation ratio-the proportion of RNA in the active vs. inactive states-for both ON and OFF riboswitches that operate near thermodynamic equilibrium: 1+[I]/K_d^I, where [I] is the input ligand concentration and K_d^I is the intrinsic dissociation constant of the aptamer module toward the input ligand. A survey of published studies of natural and synthetic riboswitches confirms that the vast majority of empirically measured activation ratios have remained well below this thermodynamic limit. A few natural and synthetic riboswitches achieve activation ratios close to the limit, and these molecules highlight important principles for achieving high riboswitch performance. For several applications, including "light-up" fluorescent sensors and chemically-controlled CRISPR/Cas complexes, the thermodynamic limit has not yet been achieved, suggesting that current tools are operating at suboptimal efficiencies. Future riboswitch studies will benefit from comparing observed activation ratios to this simple expression for the optimal activation ratio. We present experimental and computational suggestions for how to make these quantitative comparisons and suggest new molecular mechanisms that may allow non-equilibrium riboswitches to surpass the derived limit.

Structural prediction of RNA switches using conditional base-pair probabilities

An RNA switch triggers biological functions by toggling between two conformations. RNA switches include bacterial riboswitches, where ligand binding can stabilize a bound structure. For RNAs with only one stable structure, structural prediction usually just requires a straightforward free energy minimization, but for an RNA switch, the prediction of a less stable alternative structure is often computationally costly and even problematic. The current sampling-clustering method predicts stable and alternative structures by partitioning structures sampled from the energy landscape into two clusters, but it is very time-consuming. Instead, we predict the alternative structure of an RNA switch from conditional probability calculations within the energy landscape. First, our method excludes base pairs related to the most stable structure in the energy landscape. Then, it detects stable stems (“seeds”) in the remaining landscape. Finally, it folds an alternative structure prediction around a seed. While having comparable riboswitch classification performance, the conditional-probability computations had fewer adjustable parameters, offered greater predictive flexibility, and were more than one thousand times faster than the sampling step alone in sampling-clustering predictions, the competing standard. Overall, the described approach helps traverse thermodynamically improbable energy landscapes to find biologically significant substructures and structures rapidly and effectively.

Bacillus subtilis FolE is sustained by the ZagA zinc metallochaperone and the alarmone ZTP under conditions of zinc deficiency

Bacteria tightly regulate intracellular zinc levels to ensure sufficient zinc to support essential functions, while preventing toxicity. The bacterial response to zinc limitation includes the expression of putative zinc metallochaperones belonging to subfamily 1 of the COG0523 family of G3E GTPases. However, the client proteins and the metabolic processes served by these chaperones are unclear. Here, we demonstrate that the Bacillus subtilis YciC zinc metallochaperone (here renamed ZagA for ZTP activated GTPase A) supports de novo folate biosynthesis under conditions of zinc limitation, and interacts directly with the zinc-dependent GTP cyclohydrolase IA, FolE (GCYH-IA). Furthermore, we identify a role for the alarmone ZTP, a modified purine biosynthesis intermediate, in the response to zinc limitation. ZTP, a signal of 10-formyl-tetrahydrofolate (10f-THF) deficiency in bacteria, transiently accumulates as FolE begins to fail, stimulates the interaction between ZagA and FolE, and thereby helps to sustain folate synthesis despite declining zinc availability.

Riboswitch distribution, structure, and function in bacteria

Riboswitches are gene control elements that directly bind to specific ligands to regulate gene expression without the need for proteins. They are found in all three domains of life, including Bacteria, Archaea, and Eukaryota. Riboswitches are mostly spread in bacteria and archaea. In this paper, we discuss the general distribution, structure, and function of 28 different riboswitch classes as we focus our attention on riboswitches in bacteria. Bacterial riboswitches regulate gene expression by four distinct mechanisms. They regulate the expression of a limited number of genes. However, most of these genes are responsible for the synthesis of essential metabolites without which the cell cannot function. Therefore, riboswitch distribution is also important for antibacterial drug development.

Going beyond base-pairs: topology-based characterization of base-multiplets in RNA

Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.

Single-Molecule FRET Kinetics of the Mn2+ Riboswitch: Evidence for Allosteric Mg2+ Control of "Induced-Fit" vs "Conformational Selection" Folding Pathways

Gene expression in bacteria is often regulated dynamically by conformational changes in a riboswitch upon ligand binding, a detailed understanding of which is very much in its infancy. For example, the manganese riboswitch is a widespread RNA motif that conformationally responds in regulating bacterial gene expression to micromolar levels of its eponymous ligand, Mn²⁺, but the mechanistic pathways are poorly understood. In this work, we quantitatively explore the dynamic folding behavior of the manganese riboswitch by single-molecule fluorescence resonance energy transfer spectroscopy as a function of cation/ligand conditions. From the detailed analysis of the kinetics, the Mn²⁺ is shown to fold the riboswitch by a "bind-then-fold" (i.e., "induced-fit", IF) mechanism, whereby the ligand binds first and then promotes folding. On the other hand, the data also clearly reveal the presence of a folded yet ligand-free structure predominating due to the addition of physiological Mg²⁺ to a nonselective metal ion binding site. Of particular kinetic interest, such a Mg²⁺ "prefolded" conformation of the riboswitch is shown to exhibit a significantly increased affinity for Mn²⁺ and further stabilization by subsequent binding of the ligand, thereby promoting efficient riboswitch folding by a "fold-then-bind" (i.e., "conformational selection", CS) mechanism. Our results not only demonstrate Mg²⁺-controlled switching between IF and CS riboswitch folding pathways but also suggest a novel heterotropic allosteric control in the manganese riboswitch activity co-regulated by Mg²⁺ binding.

Small-Molecule-Binding Riboswitches

RNA is a versatile biomolecule capable of transferring information, taking on distinct three-dimensional shapes, and reacting to ambient conditions. RNA molecules utilize a wide range of mechanisms to control gene expression. An example of such regulation is riboswitches. Consisting exclusively of RNA, they are able to control important metabolic processes, thus providing an elegant and efficient RNA-only regulation system. Existing across all domains of life, riboswitches appear to represent one of the most highly conserved mechanisms for the regulation of a broad range of biochemical pathways. Through binding of a wide range of small-molecule ligands to their so-called aptamer domain, riboswitches undergo a conformational change in their downstream "expression platform." In consequence, the pattern of gene expression changes, which in turn results in increased or decreased protein production. Riboswitches unite the sensing and transduction of a signal that can directly be coupled to the metabolism of the cell; thus they constitute a very potent regulatory mechanism for many organisms. Highly specific RNA-binding domains not only occur in vivo but can also be evolved by means of the SELEX (systematic evolution of ligands by exponential enrichment) method, which allows in vitro selection of aptamers against almost any ligand. Coupling of these aptamers with an expression platform has led to the development of synthetic riboswitches, a highly active research field of great relevance and immense potential. The aim of this review is to summarize developments in the riboswitch field over the last decade and address key questions of recent research.

Predicting RNA-Metal Ion Binding with Ion Dehydration Effects

Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg²⁺ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.

Genome-wide Analysis of the Distribution of Riboswitches and Function Analyses of the Corresponding Downstream Genes in Prokaryotes

Background:

Riboswitches are structured elements that usually reside in the noncoding regions of mRNAs, with which various ligands bind to control a wide variety of downstream gene expressions. To date, more than twenty different classes of riboswitches have been characterized to sense various metabolites, including purines and their derivatives, coenzymes, amino acids, and metal ions, etc. </P><P> Objective: This study aims to study the genome-wide analysis of the distribution of riboswitches and function analyses of the corresponding downstream genes in prokaryotes.

Results:

In this study, we have completed a genome context analysis of 27 riboswitches to elucidate their metabolic capacities of riboswitch-mediated gene regulation from the completely-sequenced 3,079 prokaryotic genomes. Furthermore, Cluster of Orthologous Groups of proteins (COG) annotation was applied to predict and classify the possible functions of corresponding downstream genes of these riboswitches. We found that they could all be successfully annotated and grouped into 20 different COG functional categories, in which the two main clusters &quot;coenzyme metabolism [H]&quot; and &quot;amino acid transport and metabolism [E]&quot; were the most significantly enriched.

Conclusion:

Riboswitches are found to be widespread in bacteria, among which three main classes of TPP-, cobalamin- and SAM-riboswitch were the most widely distributed. We found a wide variety of functions were associated with the corresponding downstream genes, suggesting that a wide extend of regulatory roles were mediated by these riboswitches in prokaryotes.

Multiple prebiotic metals mediate translation

Today, Mg2+ is an essential cofactor with diverse structural and functional roles in life's oldest macromolecular machine, the translation system. We tested whether ancient Earth conditions (low O2, high Fe2+, and high Mn2+) can revert the ribosome to a functional ancestral state. First, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) was used to compare the effect of Mg2+, Fe2+, and Mn2+ on the tertiary structure of rRNA. Then, we used in vitro translation reactions to test whether Fe2+ or Mn2+ could mediate protein production, and quantified ribosomal metal content. We found that (i) Mg2+, Fe2+, and Mn2+ had strikingly similar effects on rRNA folding; (ii) Fe2+ and Mn2+ can replace Mg2+ as the dominant divalent cation during translation of mRNA to functional protein; and (iii) Fe and Mn associate extensively with the ribosome. Given that the translation system originated and matured when Fe2+ and Mn2+ were abundant, these findings suggest that Fe2+ and Mn2+ played a role in early ribosomal evolution.

Structural Studies of the 3',3'-cGAMP Riboswitch Induced by Cognate and Noncognate Ligands Using Molecular Dynamics Simulation

Riboswtich RNAs can control gene expression through the structural change induced by the corresponding small-molecule ligands. Molecular dynamics simulations and free energy calculations on the aptamer domain of the 3′,3′-cGAMP riboswitch in the ligand-free, cognate-bound and noncognate-bound states were performed to investigate the structural features of the 3′,3′-cGAMP riboswitch induced by the 3′,3′-cGAMP ligand and the specificity of ligand recognition. The results revealed that the aptamer of the 3′,3′-cGAMP riboswitch in the ligand-free state has a smaller binding pocket and a relatively compact structure versus that in the 3′,3′-cGAMP-bound state. The binding of the 3′,3′-cGAMP molecule to the 3′,3′-cGAMP riboswitch induces the rotation of P1 helix through the allosteric communication from the binding sites pocket containing the J1/2, J1/3 and J2/3 junction to the P1 helix. Simultaneously, these simulations also revealed that the preferential binding of the 3′,3′-cGAMP riboswitch to its cognate ligand, 3′,3′-cGAMP, over its noncognate ligand, c-di-GMP and c-di-AMP. The J1/2 junction in the 3′,3′-cGAMP riboswitch contributing to the specificity of ligand recognition have also been found.

Classification of riboswitch sequences using k-mer frequencies

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

Molecular dynamics simulation on the allosteric analysis of the c-di-GMP class I riboswitch induced by ligand binding

Riboswitches are RNA molecules that regulate gene expression using conformation change, affected by binding of small molecule ligands. Although a number of ligand-bound aptamer complex structures have been solved, it is important to know ligand-free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, we use dynamics simulations on a series of models to characterize the ligand-free and ligand-bound aptamer domain of the c-di-GMP class I (GEMM-I) riboswitch. The results revealed that the ligand-free aptamer has a stable state with a folded P2 and P3 helix, an unfolded P1 helix and open binding pocket. The first Mg ions binding to the aptamer is structurally favorable for the successive c-di-GMP binding. The P1 helix forms when c-di-GMP is successive bound. Three key junctions J1/2, J2/3 and J1/3 in the GEMM-I riboswitch contributing to the formation of P1 helix have been found. The binding of the c-di-GMP ligand to the GEMM-I riboswitch induces the riboswitch's regulation through the direct allosteric communication network in GEMM-I riboswitch from the c-di-GMP binding sites in the J1/2 and J1/3 junctions to the P1 helix, the indirect ones from those in the J2/3 and P2 communicating to P1 helix via the J1/2 and J1/3 media.

Convergent Use of Heptacoordination for Cation Selectivity by RNA and Protein Metalloregulators

The large yybP-ykoY family of bacterial riboswitches is broadly distributed phylogenetically. Previously, these gene-regulatory RNAs were proposed to respond to Mn²⁺. X-ray crystallography revealed a binuclear cation-binding pocket. This comprises one hexacoordinate site, with six oxygen ligands, which preorganizes the second, with five oxygen and one nitrogen ligands. The relatively soft nitrogen ligand was proposed to confer affinity for Mn²⁺, but how this excludes other soft cations remained enigmatic. By subjecting representative yybP-ykoY riboswitches to diverse cations in vitro, we now find that these RNAs exhibit limited transition metal ion selectivity. Among the cations tested, Cd²⁺ and Mn²⁺ bind most tightly, and comparison of three new Cd²⁺-bound crystal structures suggests that these riboswitches achieve selectivity by enforcing heptacoordination (favored by high-spin Cd²⁺ and Mn²⁺, but otherwise uncommon) in the softer site. Remarkably, the Cd²⁺- and Mn²⁺-selective bacterial transcription factor MntR also uses heptacoordination within a binuclear site to achieve selectivity.

Blind prediction of noncanonical RNA structure at atomic accuracy

Prediction of RNA structure from nucleotide sequence remains an unsolved grand challenge of biochemistry and requires distinct concepts from protein structure prediction. Despite extensive algorithmic development in recent years, modeling of noncanonical base pairs of new RNA structural motifs has not been achieved in blind challenges. We report a stepwise Monte Carlo (SWM) method with a unique add-and-delete move set that enables predictions of noncanonical base pairs of complex RNA structures. A benchmark of 82 diverse motifs establishes the method's general ability to recover noncanonical pairs ab initio, including multistrand motifs that have been refractory to prior approaches. In a blind challenge, SWM models predicted nucleotide-resolution chemical mapping and compensatory mutagenesis experiments for three in vitro selected tetraloop/receptors with previously unsolved structures (C7.2, C7.10, and R1). As a final test, SWM blindly and correctly predicted all noncanonical pairs of a Zika virus double pseudoknot during a recent community-wide RNA-Puzzle. Stepwise structure formation, as encoded in the SWM method, enables modeling of noncanonical RNA structure in a variety of previously intractable problems.

Functional roles of Mg2+ binding sites in ion-dependent gating of a Mg2+ channel, MgtE, revealed by solution NMR

Magnesium ions (Mg²⁺) are divalent cations essential for various cellular functions. Mg²⁺ homeostasis is maintained through Mg²⁺ channels such as MgtE, a prokaryotic Mg²⁺ channel whose gating is regulated by intracellular Mg²⁺ levels. Our previous crystal structure of MgtE in the Mg²⁺-bound, closed state revealed the existence of seven crystallographically-independent Mg²⁺-binding sites, Mg1-Mg7. The role of Mg²⁺-binding to each site in channel closure remains unknown. Here, we investigated Mg²⁺-dependent changes in the structure and dynamics of MgtE using nuclear magnetic resonance spectroscopy. Mg²⁺-titration experiments, using wild-type and mutant forms of MgtE, revealed that the Mg²⁺ binding sites Mg1, Mg2, Mg3, and Mg6, exhibited cooperativity and a higher affinity for Mg²⁺, enabling the remaining Mg²⁺ binding sites, Mg4, Mg5, and Mg7, to play important roles in channel closure. This study revealed the role of each Mg²⁺-binding site in MgtE gating, underlying the mechanism of cellular Mg²⁺ homeostasis.

Comparative genomic analysis of fungal TPP-riboswitches

Riboswitches are conserved RNA structures located in non-coding regions of mRNA and able to bind small molecules (e.g. metabolites) changing conformation upon binding. This feature enables them to function as regulators of gene expression. The thiamin pyrophosphate (TPP) riboswitch is the only type of riboswitches found not only in bacteria, but also in eukaryotes - in plants, green algae, protists, and fungi. Two main mechanisms of fungal TPP riboswitch action, involving alternative splicing, have been established so far. Here, we report a large-scale bioinformatic study of riboswitch structural features, action mechanisms, and distribution along the fungal taxonomy groups. For each putatively regulated gene, we reconstruct the riboswitch structure, identify other components of the regulation machinery, and establish mechanisms of riboswitch-mediated regulation. In addition to three genes known to be regulated by TPP riboswitches, thiazole synthase THI4, hydroxymethilpyrimidine-syntase NMT1, and putative transporter NCU01977, we identify two new genes, a putative thiamin transporter THI9 and a transporter of unknown specificity. While the riboswitch sequence and structure remain highly conserved in all species and genes, the mode of riboswitch-mediated regulation varies between regulated genes. The riboswitch usage varies strongly between fungal taxa, with the largest number of riboswitch-regulated genes found in Pezizomycotina and no riboswitch-mediated regulation established in Saccaromycotina.

A Glycine Riboswitch in Streptococcus pyogenes Controls Expression of a Sodium:Alanine Symporter Family Protein Gene

Regulatory RNAs play important roles in the control of bacterial gene expression. In this study, we investigated gene expression regulation by a putative glycine riboswitch located in the 5'-untranslated region of a sodium:alanine symporter family (SAF) protein gene in the group A Streptococcus pyogenes serotype M49 strain 591. Glycine-dependent gene expression mediated by riboswitch activity was studied using a luciferase reporter gene system. Maximal reporter gene expression was observed in the absence of glycine and in the presence of low glycine concentrations. Differences in glycine-dependent gene expression were not based on differential promoter activity. Expression of the SAF protein gene and the downstream putative cation efflux protein gene was investigated in wild-type bacteria by RT-qPCR transcript analyses. During growth in the presence of glycine (≥1 mM), expression of the genes were downregulated. Northern blot analyses revealed premature transcription termination in the presence of high glycine concentrations. Growth in the presence of 0.1 mM glycine led to the production of a full-length transcript. Furthermore, stability of the SAF protein gene transcript was drastically reduced in the presence of glycine. We conclude that the putative glycine riboswitch in S. pyogenes serotype M49 strain 591 represses expression of the SAF protein gene and the downstream putative cation efflux protein gene in the presence of high glycine concentrations. Sequence and secondary structure comparisons indicated that the streptococcal riboswitch belongs to the class of tandem aptamer glycine riboswitches.

LncRNA Structural Characteristics in Epigenetic Regulation

The rapid development of new generation sequencing technology has deepened the understanding of genomes and functional products. RNA-sequencing studies in mammals show that approximately 85% of the DNA sequences have RNA products, for which the length greater than 200 nucleotides (nt) is called long non-coding RNAs (lncRNA). LncRNAs now have been shown to play important epigenetic regulatory roles in key molecular processes, such as gene expression, genetic imprinting, histone modification, chromatin dynamics, and other activities by forming specific structures and interacting with all kinds of molecules. This paper mainly discusses the correlation between the structure and function of lncRNAs with the recent progress in epigenetic regulation, which is important to the understanding of the mechanism of lncRNAs in physiological and pathological processes.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Design, Optimization and Application of Small Molecule Biosensor in Metabolic Engineering

The development of synthetic biology and metabolic engineering has painted a great future for the bio-based economy, including fuels, chemicals, and drugs produced from renewable feedstocks. With the rapid advance of genome-scale modeling, pathway assembling and genome engineering/editing, our ability to design and generate microbial cell factories with various phenotype becomes almost limitless. However, our lack of ability to measure and exert precise control over metabolite concentration related phenotypes becomes a bottleneck in metabolic engineering. Genetically encoded small molecule biosensors, which provide the means to couple metabolite concentration to measurable or actionable outputs, are highly promising solutions to the bottleneck. Here we review recent advances in the design, optimization and application of small molecule biosensor in metabolic engineering, with particular focus on optimization strategies for transcription factor (TF) based biosensors.

Crystal structure of a novel RNA motif that allows for precise positioning of a single metal ion

We have determined a crystal structure of an RNA duplex containing a novel metal-binding motif. The motif is composed of two sheared G○A base pairs, two unpaired A residues and four phosphate groups in close proximity. Four A residues make an A-A-A-A stacking column at the minor groove side and two G bases are highly inclined, thereby forming the pocket-shaped motif at the major groove side. In the present structure, a hydrated Sr²⁺ ion exists in the pocket and binds to the O6 and N7 atoms of the two G bases and four phosphate groups. According to the previously-reported metal-binding properties to RNA molecules, many of divalent cations, such as Mg²⁺, Mn²⁺, Co²⁺, Zn²⁺, Ba²⁺, Pb²⁺ and Cd²⁺, may bind to the motif. This metal-binding motif can be used as a modular building block that allows for precise positioning of a single metal ion in functional nucleic acid molecules.