Riboswitch RNAs: using RNA to sense cellular metabolism

Light Control in Microbial Systems

Light is a key environmental component influencing many biological processes, particularly in prokaryotes such as archaea and bacteria. Light control techniques have revolutionized precise manipulation at molecular and cellular levels in recent years. Bacteria, with adaptability and genetic tractability, are promising candidates for light control studies. This review investigates the mechanisms underlying light activation in bacteria and discusses recent advancements focusing on light control methods and techniques for controlling bacteria. We delve into the mechanisms by which bacteria sense and transduce light signals, including engineered photoreceptors and light-sensitive actuators, and various strategies employed to modulate gene expression, protein function, and bacterial motility. Furthermore, we highlight recent developments in light-integrated methods of controlling microbial responses, such as upconversion nanoparticles and optical tweezers, which can enhance the spatial and temporal control of bacteria and open new horizons for biomedical applications.

Evolution of the RNA world: From signals to codes

The accumulated material in evolutionary biology, greatly enhanced by the achievements of modern synthetic biology, allows us to envision certain key hypothetical stages of prebiotic (chemical) evolution. This is often understood as the further evolution in the RNA World towards the RNA-protein World. It is a path towards the emergence of translation and the genetic code (I), signaling pathways with signaling molecules (II), and the appearance of RNA-based components of future gene regulatory networks (III). We believe that these evolutionary paths can be constructively viewed from the perspective of the concept of biological codes (Barbieri, 2003). Crucial evolutionary events in these directions would involve the emergence of RNA-based adaptors. Such adaptors connect two families of functionally and chemically distinct molecules into one functional entity. The emergence of primitive translation processes is undoubtedly the major milestone in the evolutionary path towards modern life. The key aspect here is the appearance of adaptors between amino acids and their cognate triplet codons. The initial steps are believed to involve the emergence of proto-transfer RNAs capable of self-aminoacylation. The second significant evolutionary breakthrough is the development of biochemical regulatory networks based on signaling molecules of the RNA World (ribonucleotides and their derivatives), as well as receptors and effectors (riboswitches) for these messengers. Some authors refer to this as the "lost language of the RNA World." The third evolutionary step is the emergence of signal sequences for ribozymes on the molecules of their RNA targets. This level of regulation in the RNA World is comparable to the gene regulatory networks of modern organisms. We believe that the signal sequences on target molecules have been rediscovered and developed by evolution into the gene regulatory networks of modern cells. In conclusion, the immense diversity of modern biological codes, in some of its key characteristics, can be traced back to the achievements of prebiotic evolution.

Mechanism of Ligand Discrimination by the NMT1 Riboswitch

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

Establishing artificial gene connections through RNA displacement-assembly-controlled CRISPR/Cas9 function

Construction of synthetic circuits that can reprogram genetic networks and signal pathways is a long-term goal for manipulation of biosystems. However, it is still highly challenging to build artificial genetic communications among endogenous RNA species due to their sequence independence and structural diversities. Here we report an RNA-based synthetic circuit that can establish regulatory linkages between expression of endogenous genes in both Escherichiacoli and mammalian cells. This design employs a displacement-assembly approach to modulate the activity of guide RNA for function control of CRISPR/Cas9. Our experiments demonstrate the great effectiveness of this RNA circuit for building artificial connections between expression of originally unrelated genes. Both exogenous and naturally occurring RNAs, including small/microRNAs and long mRNAs, are capable of controlling expression of another endogenous gene through this approach. Moreover, an artificial signal pathway inside mammalian cells is also successfully established to control cell apoptosis through our designed synthetic circuit. This study provides a general strategy for constructing synthetic RNA circuits, which can introduce artificial connections into the genetic networks of mammalian cells and alter the cellular phenotypes.

A Reduction of Transcriptional Regulation in Aquatic Oligotrophic Microorganisms Enhances Fitness in Nutrient-Poor Environments

In this review, we consider the regulatory strategies of aquatic oligotrophs, microbial cells that are adapted to thrive under low-nutrient concentrations in oceans, lakes, and other aquatic ecosystems. Many reports have concluded that oligotrophs use less transcriptional regulation than copiotrophic cells, which are adapted to high nutrient concentrations and are far more common subjects for laboratory investigations of regulation. It is theorized that oligotrophs have retained alternate mechanisms of regulation, such as riboswitches, that provide shorter response times and smaller amplitude responses and require fewer cellular resources. We examine the accumulated evidence for distinctive regulatory strategies in oligotrophs. We explore differences in the selective pressures copiotrophs and oligotrophs encounter and ask why, although evolutionary history gives copiotrophs and oligotrophs access to the same regulatory mechanisms, they might exhibit distinctly different patterns in how these mechanisms are used. We discuss the implications of these findings for understanding broad patterns in the evolution of microbial regulatory networks and their relationships to environmental niche and life history strategy. We ask whether these observations, which have emerged from a decade of increased investigation of the cell biology of oligotrophs, might be relevant to recent discoveries of many microbial cell lineages in nature that share with oligotrophs the property of reduced genome size.

A variant of guanidine-IV riboswitches exhibits evidence of a distinct ligand specificity

Riboswitches are regulatory RNAs that specifically bind a small molecule or ion. Like metabolite-binding proteins, riboswitches can evolve new ligand specificities, and some examples of this phenomenon have been validated. As part of work based on comparative genomics to discover novel riboswitches, we encountered a candidate riboswitch with striking similarities to the recently identified guanidine-IV riboswitch. This candidate riboswitch, the Gd4v motif, is predicted in four distinct bacterial phyla, thus almost as widespread as the guanidine-IV riboswitch. Bioinformatic and experimental analysis suggest that the Gd4v motif is a riboswitch that binds a ligand other than guanidine. It is found associated with gene classes that differ from genes regulated by confirmed guanidine riboswitches. In inline-probing assays, we showed that free guanidine binds only weakly to one of the tested sequences of the variant. Further tested compounds did not show binding, attenuation of transcription termination, or activation of a genetic reporter construct. We characterized an N-acetyltransferase frequently associated with the Gd4v motif and compared its substrate preference to an N-acetyltransferase that occurs under control of guanidine-IV riboswitches. The substrates of this Gd4v-motif-associated enzyme did not show activity for Gd4v RNA binding or transcription termination. Hence, the ligand of the candidate riboswitch motif remains unidentified. The variant RNA motif is predominantly found in gut metagenome sequences, hinting at a ligand that is highly relevant in this environment. This finding is a first step to determining the identity of this unknown ligand, and understanding how guanidine-IV-riboswitch-like structures can evolve to bind different ligands.

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

INTRODUCTION

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

AREAS COVERED

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

EXPERT OPINION

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

A Riboswitch-Driven Era of New Antibacterials

Riboswitches are structured non-coding RNAs found in the 5′ UTR of important genes for bacterial metabolism, virulence and survival. Upon the binding of specific ligands that can vary from simple ions to complex molecules such as nucleotides and tRNAs, riboswitches change their local and global mRNA conformations to affect downstream transcription or translation. Due to their dynamic nature and central regulatory role in bacterial metabolism, riboswitches have been exploited as novel RNA-based targets for the development of new generation antibacterials that can overcome drug-resistance problems. During recent years, several important riboswitch structures from many bacterial representatives, including several prominent human pathogens, have shown that riboswitches are ideal RNA targets for new compounds that can interfere with their structure and function, exhibiting much reduced resistance over time. Most interestingly, mainstream antibiotics that target the ribosome have been shown to effectively modulate the regulatory behavior and capacity of several riboswitches, both in vivo and in vitro, emphasizing the need for more in-depth studies and biological evaluation of new antibiotics. Herein, we summarize the currently known compounds that target several main riboswitches and discuss the role of mainstream antibiotics as modulators of T-box riboswitches, in the dawn of an era of novel inhibitors that target important bacterial regulatory RNAs.

Lineage-specific insertions in T-box riboswitches modulate antibiotic binding and action

T-box riboswitches (T-boxes) are essential RNA regulatory elements with a remarkable structural diversity, especially among bacterial pathogens. In staphylococci, all glyS T-boxes synchronize glycine supply during synthesis of nascent polypeptides and cell wall formation and are characterized by a conserved and unique insertion in their antiterminator/terminator domain, termed stem Sa. Interestingly, in Staphylococcus aureus the stem Sa can accommodate binding of specific antibiotics, which in turn induce robust and diverse effects on T-box-mediated transcription. In the present study, domain swap mutagenesis and probing analysis were performed to decipher the role of stem Sa. Deletion of stem Sa significantly reduces both the S. aureus glyS T-box-mediated transcription readthrough levels and the ability to discriminate among tRNAGly isoacceptors, both in vitro and in vivo. Moreover, the deletion inverted the previously reported stimulatory effects of specific antibiotics. Interestingly, stem Sa insertion in the terminator/antiterminator domain of Geobacillus kaustophilus glyS T-box, which lacks this domain, resulted in elevated transcription in the presence of tigecycline and facilitated discrimination among proteinogenic and nonproteinogenic tRNAGly isoacceptors. Overall, stem Sa represents a lineage-specific structural feature required for efficient staphylococcal glyS T-box-mediated transcription and it could serve as a species-selective druggable target through its ability to modulate antibiotic binding.

Adaptation by Type III CRISPR-Cas Systems: Breakthrough Findings and Open Questions

CRISPR-Cas systems acquire heritable defense memory against invading nucleic acids through adaptation. Type III CRISPR-Cas systems have unique and intriguing features of defense and are important in method development for Genetics research. We started to understand the common and unique properties of type III CRISPR-Cas adaptation in recent years. This review summarizes our knowledge regarding CRISPR-Cas adaptation with the emphasis on type III systems and discusses open questions for type III adaptation studies.

QRNAstruct: a method for extracting secondary structural features of RNA via regression with biological activity

Recent technological advances have enabled the generation of large amounts of data consisting of RNA sequences and their functional activity. Here, we propose a method for extracting secondary structure features that affect the functional activity of RNA from sequence–activity data. Given pairs of RNA sequences and their corresponding bioactivity values, our method calculates position-specific structural features of the input RNA sequences, considering every possible secondary structure of each RNA. A Ridge regression model is trained using the structural features as feature vectors and the bioactivity values as response variables. Optimized model parameters indicate how secondary structure features affect bioactivity. We used our method to extract intramolecular structural features of bacterial translation initiation sites and self-cleaving ribozymes, and the intermolecular features between rRNAs and Shine–Dalgarno sequences and between U1 RNAs and splicing sites. We not only identified known structural features but also revealed more detailed insights into structure–activity relationships than previously reported. Importantly, the datasets we analyzed here were obtained from different experimental systems and differed in size, sequence length and similarity, and number of RNA molecules involved, demonstrating that our method is applicable to various types of data consisting of RNA sequences and bioactivity values.

Imaging Approaches for the Study of Metabolism in Real Time Using Genetically Encoded Reporters

Metabolism comprises of two axes in order to serve homeostasis: anabolism and catabolism. Both axes are interbranched with the so-called bioenergetics aspect of metabolism. There is a plethora of analytical biochemical methods to monitor metabolites and reactions in lysates, yet there is a rising need to monitor, quantify and elucidate in real time the spatiotemporal orchestration of complex biochemical reactions in living systems and furthermore to analyze the metabolic effect of chemical compounds that are destined for the clinic. The ongoing technological burst in the field of imaging creates opportunities to establish new tools that will allow investigators to monitor dynamics of biochemical reactions and kinetics of metabolites at a resolution that ranges from subcellular organelle to whole system for some key metabolites. This article provides a mini review of available toolkits to achieve this goal but also presents a perspective on the open space that can be exploited to develop novel methodologies that will merge classic biochemistry of metabolism with advanced imaging. In other words, a perspective of "watching metabolism in real time."

On the predictibility of A-minor motifs from their local contexts

This study investigates the importance of the structural context in the formation of a type I/II A-minor motif. This very frequent structural motif has been shown to be important in the spatial folding of RNA molecules. We developed an automated method to classify A-minor motif occurrences according to their 3D context similarities, and we used a graph approach to represent both the structural A-minor motif occurrences and their classes at different scales. This approach leads us to uncover new subclasses of A-minor motif occurrences according to their local 3D similarities. The majority of classes are composed of homologous occurrences, but some of them are composed of non-homologous occurrences. The different classifications we obtain allow us to better understand the importance of the context in the formation of A-minor motifs. In a second step, we investigate how much knowledge of the context around an A-minor motif can help to infer its presence (and position). More specifically, we want to determine what kind of information, contained in the structural context, can be useful to characterize and predict A-minor motifs. We show that, for some A-minor motifs, the topology combined with a sequence signal is sufficient to predict the presence and the position of an A-minor motif occurrence. In most other cases, these signals are not sufficient for predicting the A-minor motif, however we show that they are good signals for this purpose. All the classification and prediction pipelines rely on automated processes, for which we describe the underlying algorithms and parameters.

Unique properties of spacer acquisition by the type III-A CRISPR-Cas system

Type III CRISPR-Cas systems have a unique mode of interference, involving crRNA-guided recognition of nascent RNA and leading to DNA and RNA degradation. How type III systems acquire new CRISPR spacers is currently not well understood. Here, we characterize CRISPR spacer uptake by a type III-A system within its native host, Streptococcus thermophilus. Adaptation by the type II-A system in the same host provided a basis for comparison. Cas1 and Cas2 proteins were critical for type III adaptation but deletion of genes responsible for crRNA biogenesis or interference did not detectably change spacer uptake patterns, except those related to host counter-selection. Unlike the type II-A system, type III spacers are acquired in a PAM- and orientation-independent manner. Interestingly, certain regions of plasmids and the host genome were particularly well-sampled during type III-A, but not type II-A, spacer uptake. These regions included the single-stranded origins of rolling-circle replicating plasmids, rRNA and tRNA encoding gene clusters, promoter regions of expressed genes and 5′ UTR regions involved in transcription attenuation. These features share the potential to form DNA secondary structures, suggesting a preferred substrate for type III adaptation. Lastly, the type III-A system adapted to and protected host cells from lytic phage infection.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Approaches in the photosynthetic production of sustainable fuels by cyanobacteria using tools of synthetic biology

Cyanobacteria, photosynthetic prokaryotic microorganisms having a simple genetic composition are the prospective photoautotrophic cell factories for the production of a wide range of biofuel molecules. The simple genetic composition of cyanobacteria allows effortless genetic manipulation which leads to increased research endeavors from the synthetic biology approach. Various unicellular model cyanobacterial strains like Synechocystis sp. PCC 6803 and Synechococcus elongatus PCC 7942 have been successfully engineered for biofuels generation. Improved development of synthetic biology tools, genetic modification methods and advancement in transformation techniques to construct a strain that can contain multiple foreign genes in a single operon have vastly expanded the functions that can be used for engineering photosynthetic cyanobacteria for the generation of various biofuel molecules. In this review, recent advancements and approaches in synthetic biology tools used for cyanobacterial genome editing have been discussed. Apart from this, cyanobacterial productions of various fuel molecules like isoprene, limonene, α-farnesene, squalene, alkanes, butanol, and fatty acids, which can be a substitute for petroleum and fossil fuels in the future, have been elaborated.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Modular tuning engineering and versatile applications of genetically encoded biosensors

Genetically encoded biosensors have a diverse range of detectable signals and potential applications in many fields, including metabolism control and high-throughput screening. Their ability to be used in situ with minimal interference to the bioprocess of interest could revolutionize synthetic biology and microbial cell factories. The performance and functions of these biosensors have been extensively studied and have been rapidly improved. We review here current biosensor tuning strategies and attempt to unravel how to obtain ideal biosensor functions through experimental adjustments. Strategies for expanding the biosensor input signals that increases the number of detectable compounds have also been summarized. Finally, different output signals and their practical requirements for biotechnology and biomedical applications and environmental safety concerns have been analyzed. This in-depth review of the responses and regulation mechanisms of genetically encoded biosensors will assist to improve their design and optimization in various application scenarios.

Long Tracts of Guanines Drive Aggregation of RNA G-Quadruplexes in the Presence of Spermine

G-Quadruplexes (GQs) are compact, stable structures in DNA and RNA comprised of two or more tiers of quartets whose G-rich motif of tracts of two or more G's occurs commonly within genomes and transcriptomes. While thermodynamically stable in vitro, these structures remain difficult to study in vivo. One approach to understanding GQ in vivo behavior is to test whether conditions and molecules found in cells facilitate their folding. Polyamines are biogenic polycations that interact with RNA. Among common polyamines, spermine contains the highest charge and is found in eukaryotes, making it a good candidate for association with high-charge density nucleic acid structures like GQs. Using a variety of techniques, including ultraviolet-detected thermal denaturation, circular dichroism, size exclusion chromatography, and confocal microscopy, on an array of quadruplex sequence variants, we find that eukaryotic biological concentrations of spermine induce microaggregation of three-tiered G-rich sequences, but not of purely two-tiered structures, although higher spermine concentrations induce aggregation of even these. The formation of microaggregates can also be induced by addition of as little as a single G to a two-tiered structure; moreover, they form at biological temperatures, are sensitive to salt, and can form in the presence of at least some flanking sequence. Notably, GQ aggregation is not observed under prokaryotic-like conditions of no spermine and higher NaCl concentrations. The sequence, polyamine, and salt specificity of microaggregation reported herein have implications for the formation and stability of G-rich nucleic acid aggregates in vivo and for functional roles for understudied GQ sequences with only two quadruplex tiers.

The RNA structurome in the asexual blood stages of malaria pathogen plasmodium falciparum

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite's asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5' UTR, CDS and 3' UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

Synthesis of 5-Cyanomethyluridine (cnm5 U) and 5-Cyanouridine (cn5 U) Phosphoramidites and Their Incorporation into RNA Oligonucleotides

This article contains detailed synthetic protocols for preparation of 5-cyanomethyluridine (cnm⁵ U) and 5-cyanouridine (cn⁵ U) phosphoramidites. The synthesis of the cnm⁵ U phosphoramidite building block starts with commercially available 5-methyluridine (m⁵ C), followed by bromination of the 5-methyl group to install the cyano moiety using TMSCN/TBAF. The cn⁵ U phosphoramidite is obtained by regular Vorbrüggen glycosylation of the protected ribofuranose with silylated 5-cyanouracil. These two modified phosphoramidites are suitable for synthesis of RNA oligonucleotides on solid phase using conventional amidite chemistry. Our protocol provides access to two novel building blocks for constructing RNA-based therapeutics. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Preparation of cnm⁵ U and cn⁵ U phosphoramidites Basic Protocol 2: Synthesis, purification, and characterization of cnm⁵ U- and cn⁵ U-modified RNA oligonucleotides.

Regulation of Glycine Cleavage and Detoxification by a Highly Conserved Glycine Riboswitch in Burkholderia spp

The glycine riboswitch is a known regulatory element that is unique in having two aptamers that are joined by a linker region. In this study, we investigated a glycine riboswitch located in the 5' untranslated region of a glycine cleavage system homolog (gcvTHP) in Burkholderia spp. Structure prediction using the sequence generated a model with a glycine binding pocket composed of base-triple interactions (G62-A64-A86 and G65-U84-C85) that are supported by A/G minor interactions (A17-C60-G88 and G16-C61-G87, respectively) and two ribose-zipper motifs (C11-G12 interacting with A248-A247 and C153-U154 interacting with A79-A78) which had not been previously reported. The capacity of the riboswitch to bind to glycine was experimentally validated by native gel assays and the crucial role of interactions that make up the glycine binding pocket were proven by mutations of A17U and G16C which resulted in conformational differences that may lead to dysfunction. Using glycine supplemented minimal media, we were able to prove that the expression of the gcvTHP genes found downstream of the riboswitch responded to the glycine concentrations introduced thus confirming the role of this highly conserved Burkholderia riboswitch and its associated genes as a putative glycine detoxification system in Burkholderia spp.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures

RNA structures play a fundamental role in nearly every aspect of cellular physiology and pathology. Gaining insights into the functions of RNA molecules requires accurate predictions of RNA secondary structures. However, the existing thermodynamic folding models remain less accurate than desired, even when chemical probing data, such as selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) reactivities, are used as restraints. Unlike most SHAPE-directed algorithms that only consider SHAPE restraints for base pairing, we extract two-dimensional structural features encoded in SHAPE data and establish robust relationships between characteristic SHAPE patterns and loop motifs of various types (hairpin, internal, and bulge) and lengths (2-11 nucleotides). Such characteristic SHAPE patterns are closely related to the sugar pucker conformations of loop residues. Based on these patterns, we propose a computational method, SHAPELoop, which refines the predicted results of the existing methods, thereby further improving their prediction accuracy. In addition, SHAPELoop can provide information about local or global structural rearrangements (including pseudoknots) and help researchers to easily test their hypothesized secondary structures.

De Novo Cobalamin Biosynthesis, Transport, and Assimilation and Cobalamin-Mediated Regulation of Methionine Biosynthesis in Mycobacterium smegmatis

Alterations in cobalamin-dependent metabolism have marked the evolution of Mycobacterium tuberculosis into a human pathogen. However, the role(s) of cobalamin in mycobacterial physiology remain poorly understood.

Cobalamin is an essential cofactor in all domains of life, yet its biosynthesis is restricted to some bacteria and archaea. Mycobacterium smegmatis, an environmental saprophyte frequently used as surrogate for the obligate human pathogen M. tuberculosis, carries approximately 30 genes predicted to be involved in de novo cobalamin biosynthesis. M. smegmatis also encodes multiple cobalamin-dependent enzymes, including MetH, a methionine synthase that catalyzes the final reaction in methionine biosynthesis. In addition to metH, M. smegmatis possesses a cobalamin-independent methionine synthase, metE, suggesting that enzyme use—MetH versus MetE—is regulated by cobalamin availability. Consistent with this notion, we previously described a cobalamin-sensing riboswitch controlling metE expression in M. tuberculosis. Here, we apply a targeted mass spectrometry-based approach to confirm de novo cobalamin biosynthesis in M. smegmatis during aerobic growth in vitro. We also demonstrate that M. smegmatis can transport and assimilate exogenous cyanocobalamin (CNCbl; also known as vitamin B₁₂) and its precursor, dicyanocobinamide ([CN]₂Cbi). However, the uptake of CNCbl and (CN)₂Cbi in this organism is restricted and seems dependent on the conditional essentiality of the cobalamin-dependent methionine synthase. Using gene and protein expression analyses combined with single-cell growth kinetics and live-cell time-lapse microscopy, we show that transcription and translation of metE are strongly attenuated by endogenous cobalamin. These results support the inference that metH essentiality in M. smegmatis results from riboswitch-mediated repression of MetE expression. Moreover, differences observed in cobalamin-dependent metabolism between M. smegmatis and M. tuberculosis provide some insight into the selective pressures which might have shaped mycobacterial metabolism for pathogenicity.

IMPORTANCE Alterations in cobalamin-dependent metabolism have marked the evolution of Mycobacterium tuberculosis into a human pathogen. However, the role(s) of cobalamin in mycobacterial physiology remains poorly understood. Using the nonpathogenic saprophyte M. smegmatis, we investigated the production of cobalamin, transport and assimilation of cobalamin precursors, and the role of cobalamin in regulating methionine biosynthesis. We confirm constitutive de novo cobalamin biosynthesis in M. smegmatis, in contrast with M. tuberculosis, which appears to lack de novo cobalamin biosynthetic capacity. We also show that uptake of cyanocobalamin (vitamin B₁₂) and its precursors is restricted in M. smegmatis, apparently depending on the cofactor requirements of the cobalamin-dependent methionine synthase. These observations establish M. smegmatis as an informative foil to elucidate key metabolic adaptations enabling mycobacterial pathogenicity.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

RNA-controlled regulation in Caulobacter crescentus

In the past decades, Caulobacter crescentus has been extensively studied, mostly regarding its dimorphic, asymmetric life cycle. Its distinct mode of reproduction and the need to optimally adapt to ever-changing environmental conditions require tight coordination of gene regulation. Post-transcriptional regulation through non-coding RNAs and RNA-binding proteins constitutes an important layer of expression control in bacteria, but its principles and mechanisms in Caulobacter have only recently been explored. RNA-binding proteins including the RNA chaperone Hfq and ribonuclease RNase E contribute to the activity of regulatory RNAs. Riboswitches and RNA thermometers govern expression of downstream open reading frames, while the small regulatory RNAs CrfA, ChvR and GsrN instead control targets encoded in trans by direct base-pairing interactions.

A New Promising Anti-Infective Agent Inhibits Biofilm Growth by Targeting Simultaneously a Conserved RNA Function That Controls Multiple Genes

Combating single and multi-drug-resistant infections in the form of biofilms is an immediate challenge. The challenge is to discover innovative targets and develop novel chemistries that combat biofilms and drug-resistant organisms, and thwart emergence of future resistant strains. An ideal novel target would control multiple genes, and can be inhibited by a single compound. We previously demonstrated success against Staphylococcus aureus biofilms by targeting the tRNA-dependent regulated T-box genes, not present in the human host. Present in Gram-positive bacteria, T-box genes attenuate transcription with a riboswitch-like element that regulates the expression of aminoacyl-tRNA synthetases and amino acid metabolism genes required for cell viability. PKZ18, the parent of a family of compounds selected in silico from 305,000 molecules, inhibits the function of the conserved T-box regulatory element and thus blocks growth of antibiotic-resistant S. aureus in biofilms. The PKZ18 analog PKZ18-22 was 10-fold more potent than vancomycin in inhibiting growth of S. aureus in biofilms. In addition, PKZ18-22 has a synergistic effect with existing antibiotics, e.g., gentamicin and rifampin. PKZ18-22 inhibits the T-box regulatory mechanism, halts the transcription of vital genes, and results in cell death. These effects are independent of the growth state, planktonic or biofilm, of the bacteria, and could inhibit emergent strains.

Riboswitch-Mediated Detection of Metabolite Fluctuations During Live Cell Imaging of Bacteria

Riboswitches are a class of noncoding RNAs that regulate gene expression in response to changes in intracellular metabolite concentrations. When riboswitches are placed upstream of genetic reporters, the degree of reporter activity reflects the relative abundance of the metabolite that is sensed by the riboswitch. This method describes how reporters for live cell imaging, such as yellow fluorescent protein (YFP), can be placed under genetic control by metabolite-sensing riboswitches in the bacterium Bacillus subtilis. Specifically, a protocol for generating a fluorescent YFP reporter, based on a c-di-GMP responsive riboswitch, is outlined below.

From structure to function: Route to understanding lncRNA mechanism

RNAs have emerged as a major target for diagnostics and therapeutics approaches. Regulatory nonprotein-coding RNAs (ncRNAs) in particular display remarkable versatility. They can fold into complex structures and interact with proteins, DNA, and other RNAs, thus modulating activity, localization, or interactome of multi-protein complexes. Thus, ncRNAs confer regulatory plasticity and represent a new layer of regulatory control. Interestingly, long noncoding RNAs (lncRNAs) tend to acquire complex secondary and tertiary structures and their function-in many cases-is dependent on structural conservation rather than primary sequence conservation. Whereas for many proteins, structure and its associated function are closely connected, for lncRNAs, the structural domains that determine functionality and its interactome are still not well understood. Numerous approaches for analyzing the structural configuration of lncRNAs have been developed recently. Here, will provide an overview of major experimental approaches used in the field, and discuss the potential benefit of using combinatorial strategies to analyze lncRNA modes of action based on structural information.

How Pathogens Feel and Overcome Magnesium Limitation When in Host Tissues

Host organisms utilize nutritional immunity to limit the availability of nutrients essential to an invading pathogen. Nutrients may include amino acids, nucleotide bases, and transition metals, the essentiality of which varies among pathogens. The mammalian macrophage protein Slc11a1 (previously Nramp1) mediates resistance to several intracellular pathogens. Slc11a1 is proposed to restrict growth of Salmonella enterica serovar Typhimurium in host tissues by causing magnesium deprivation. This is intriguing because magnesium is the most abundant divalent cation in all living cells. A pathogen's response to factors such as Slc11a1 that promote nutritional immunity may therefore reflect what the pathogen 'feels' in its cytoplasm, rather than the nutrient concentration in host cell compartments.

First Come, First Served: Sui Generis Features of the First Intron

Most of the transcribed genes in eukaryotic cells are interrupted by intervening sequences called introns that are co-transcriptionally removed from nascent messenger RNA through the process of splicing. In Arabidopsis, 79% of genes contain introns and more than 60% of intron-containing genes undergo alternative splicing (AS), which ostensibly is considered to increase protein diversity as one of the intrinsic mechanisms for fitness to the varying environment or the internal developmental program. In addition, recent findings have prevailed in terms of overlooked intron functions. Here, we review recent progress in the underlying mechanisms of intron function, in particular by focusing on unique features of the first intron that is located in close proximity to the transcription start site. The distinct deposition of epigenetic marks and nucleosome density on the first intronic DNA sequence, the impact of the first intron on determining the transcription start site and elongation of its own expression (called intron-mediated enhancement, IME), translation control in 5′-UTR, and the new mechanism of the trans-acting function of the first intron in regulating gene expression at the post-transcriptional level are summarized.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

A riboswitch gives rise to multi-generational phenotypic heterogeneity in an auxotrophic bacterium

Auxotrophy, the inability to produce an organic compound essential for growth, is widespread among bacteria. Auxotrophic bacteria rely on transporters to acquire these compounds from their environment. Here, we study the expression of both low- and high-affinity transporters of the costly amino acid methionine in an auxotrophic lactic acid bacterium, Lactococcus lactis. We show that the high-affinity transporter (Met-transporter) is heterogeneously expressed at low methionine concentrations, resulting in two isogenic subpopulations that sequester methionine in different ways: one subpopulation primarily relies on the high-affinity transporter (high expression of the Met-transporter) and the other subpopulation primarily relies on the low-affinity transporter (low expression of the Met-transporter). The phenotypic heterogeneity is remarkably stable, inherited for tens of generations, and apparent at the colony level. This heterogeneity results from a T-box riboswitch in the promoter region of the met operon encoding the high-affinity Met-transporter. We hypothesize that T-box riboswitches, which are commonly found in the Lactobacillales, may play as-yet unexplored roles in the predominantly auxotrophic lifestyle of these bacteria.

In silico prediction and expression profile analysis of small non-coding RNAs in Herbaspirillum seropedicae SmR1

BACKGROUND

Herbaspirillum seropedicae is a diazotrophic bacterium from the β-proteobacteria class that colonizes endophytically important gramineous species, promotes their growth through phytohormone-dependent stimulation and can express nif genes and fix nitrogen inside plant tissues. Due to these properties this bacterium has great potential as a commercial inoculant for agriculture. The H. seropedicae SmR1 genome is completely sequenced and annotated but despite the availability of diverse structural and functional analysis of this genome, studies involving small non-coding RNAs (sRNAs) has not yet been done. We have conducted computational prediction and RNA-seq analysis to select and confirm the expression of sRNA genes in the H. seropedicae SmR1 genome, in the presence of two nitrogen independent sources and in presence of naringenin, a flavonoid secreted by some plants.

RESULTS

This approach resulted in a set of 117 sRNAs distributed in riboswitch, cis-encoded and trans-encoded categories and among them 20 have Rfam homologs. The housekeeping sRNAs tmRNA, ssrS and 4.5S were found and we observed that a large number of sRNAs are more expressed in the nitrate condition rather than the control condition and in the presence of naringenin. Some sRNAs expression were confirmed in vitro and this work contributes to better understand the post transcriptional regulation in this bacterium.

CONCLUSIONS

H. seropedicae SmR1 express sRNAs in the presence of two nitrogen sources and/or in the presence of naringenin. The functions of most of these sRNAs remains unknown but their existence in this bacterium confirms the evidence that sRNAs are involved in many different cellular activities to adapt to nutritional and environmental changes.

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

The 5' untranslated regions (5' UTR) of mRNAs play an important role in the eukaryotic translation initiation process. Additional levels of translational regulation may be mediated through interactions between structured mRNAs that can adopt interchangeable secondary or tertiary structures and the regulatory protein/RNA factors or components of the translational apparatus. Here we report a regulatory function of the 5' UTR mRNA of the spe2 gene (SAM decarboxylase) in polyamine metabolism of the fission yeast Schizosaccharomyces pombe Reporter assays, biochemical experiments, and mutational analysis demonstrate that this 5' UTR mRNA of spe2 can bind to spermidine to regulate translation. A tertiary structure transition in the 5' UTR RNA upon spermidine binding is essential for translation regulation. This study provides biochemical evidence for spermidine binding to regulate translation of the spe2 gene through interactions with the 5' UTR mRNA. The identification of such a regulatory RNA that is directly associated with an essential eukaryotic metabolic process suggests that other ligand-binding RNAs may also contribute to eukaryotic gene regulation.

In vivo evolutionary engineering of riboswitch with high-threshold for N-acetylneuraminic acid production

Riboswitches with desired properties, such as sensitivity, threshold, dynamic range, is important for its application. However, the property change of a natural riboswitch is difficult due to the lack of the understanding of aptamer ligand binding properties and a proper screening method for both rational and irrational design. In this study, an effective method to change the threshold of riboswitch was established in vivo based on growth coupled screening by combining both positive and negative selections. The feasibility of the method was verified by the model library. Using this method, an N-acetylneuraminic acid (NeuAc) riboswitch was evolved and modified riboswitches with high threshold and large dynamic range were obtained. Then, using a new NeuAc riboswitch, both ribosome binding sites and key gene in NeuAc biosynthesis pathway were optimized. The highest NeuAc production of 14.32 g/l that has been reported using glucose as sole carbon source was obtained.

Hidden cases of tRNA gene duplication and remolding in mitochondrial genomes of amphipods

The evolution of tRNA genes in mitochondrial (mt) genomes is a complex process that includes duplications, degenerations, and transpositions, as well as a specific process of identity change through mutations in the anticodon (tRNA gene remolding or tRNA gene recruitment). Using amphipod-specific tRNA models for annotation, we show that tRNA duplications are more common in the mt genomes of amphipods than what was revealed by previous annotations. Seventeen cases of tRNA gene duplications were detected in the mt genomes of amphipods, and ten of them were tRNA genes that underwent remolding. The additional tRNA gene findings were verified using phylogenetic analysis and genetic distance analysis. The majority of remolded tRNA genes (seven out of ten cases) were found in the mt genomes of endemic amphipod species from Lake Baikal. All additional mt tRNA genes arose independently in the Baikalian amphipods, indicating the unusual plasticity of tRNA gene evolution in these species assemblages. The possible reasons for the unusual abundance of additional tRNA genes in the mt genomes of Baikalian amphipods are discussed. The amphipod-specific tRNA models developed for MiTFi refine existing predictions of tRNA genes in amphipods and reveal additional cases of duplicated tRNA genes overlooked by using less specific Metazoa-wide models. The application of these models for mt tRNA gene prediction will be useful for the correct annotation of mt genomes of amphipods and probably other crustaceans.

A Modified In Vitro Transcription Approach to Improve RNA Synthesis and Ribozyme Cleavage Efficiency

RNA elements such as catalytic RNA, riboswitch, microRNA, and long non-coding RNA perform a major role in cellular processes. A complete understanding of cellular processes is impossible without knowing the structure-function relationship of participating RNA molecules that ultimately requires large quantities of pure RNAs. Thus, structural/functional analyses of emerging RNAs necessitate revised protocols for improved RNA quantity and quality. Here we present a modified in vitro transcription protocol to enhance ribozyme cleaving efficiency and RNA yield by working on two variables, i.e., incubation temperature and limiting GTPs. Following an improved RNA synthesis, the target RNA is purified from transcription mixture components through denaturing size-exclusion chromatography. The protocol confirms that cyclic elevated incubation temperatures during transcription and increased concentrations of GTPs improve the production rate of RNA. Our modified in vitro transcription method improves the ribozyme cleaving efficiency and targets RNA yield by four- to fivefold that can benefit almost any RNA-related study from protein-RNA interaction analysis to crystallography.

Mammalian RNA switches: Molecular rheostats in gene regulation, disease, and medicine

Alteration of RNA structure by environmental signals is a fundamental mechanism of gene regulation. For example, the riboswitch is a noncoding RNA regulatory element that binds a small molecule and causes a structural change in the RNA, thereby regulating transcription, splicing, or translation of an mRNA. The role of riboswitches in metabolite sensing and gene regulation in bacteria and other lower species was reported almost two decades ago, but riboswitches have not yet been discovered in mammals. An analog of the riboswitch, the protein-directed RNA switch (PDRS), has been identified as an important regulatory mechanism of gene expression in mammalian cells. RNA-binding proteins and microRNAs are two major executors of PDRS via their interaction with target transcripts in mammals. These protein-RNA interactions influence cellular functions by integrating environmental signals and intracellular pathways from disparate stimuli to modulate stability or translation of specific mRNAs. The discovery of a riboswitch in eukaryotes that is composed of a single class of thiamine pyrophosphate (TPP) suggests that additional ligand-sensing RNAs may be present to control eukaryotic or mammalian gene expression. In this review, we focus on protein-directed RNA switch mechanisms in mammals. We offer perspectives on the potential discovery of mammalian protein-directed and compound-dependent RNA switches that are related to human disease and medicine.

On the nature and origin of biological information: The curious case of RNA

Biological information is most commonly thought of in terms of biology's Central Dogma where DNA is viewed as a linearized code used to synthesize proteins. Using DNA's chemical cousin, RNA, as a case study we consider how biological information operates outside the linear arrangement of its polymeric subunits. Much like individual pieces of a jigsaw puzzle, particular structures enable biomolecules to undergo precise molecular interactions with one another based on their respective shapes. By exploring the relationship between sequence and structure in RNA we argue that biological information finds its ultimate functional fulfillment in the three-dimensional structural arrangement of its atoms. We show how recurrent structural RNA motifs-operating at the tertiary level of a molecule-provide robust building blocks for the formation of new structural configurations and thereby convey the information required for emergent biological functions. We posit that these same RNA structures, guided by their respective thermodynamic stabilities, experience selective pressure to maintain particular three-dimensional architectures over and above pressures to maintain a particular sequence of nucleotides. Ultimately, this framework for understanding the nature of biological information provides a useful paradigm for understanding its origins and how biological information can result from chaotic prebiotic conditions.

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.

Synthetic 5' UTRs Can Either Up- or Downregulate Expression upon RNA-Binding Protein Binding

The construction of complex gene-regulatory networks requires both inhibitory and upregulatory modules. However, the vast majority of RNA-based regulatory "parts" are inhibitory. Using a synthetic biology approach combined with SHAPE-seq, we explored the regulatory effect of RNA-binding protein (RBP)-RNA interactions in bacterial 5' UTRs. By positioning a library of RNA hairpins upstream of a reporter gene and co-expressing them with the matching RBP, we observed a set of regulatory responses, including translational stimulation, translational repression, and cooperative behavior. Our combined approach revealed three distinct states in vivo: in the absence of RBPs, the RNA molecules can be found in either a molten state that is amenable to translation or a structured phase that inhibits translation. In the presence of RBPs, the RNA molecules are in a semi-structured phase with partial translational capacity. Our work provides new insight into RBP-based regulation and a blueprint for designing complete gene-regulatory circuits at the post-transcriptional level.

Redefining the Small Regulatory RNA Transcriptome in Streptococcus pneumoniae Serotype 2 Strain D39

Streptococcus pneumoniae (pneumococcus) is a major human respiratory pathogen and a leading cause of bacterial pneumonia worldwide. Small regulatory RNAs (sRNAs), which often act by posttranscriptionally regulating gene expression, have been shown to be crucial for the virulence of S. pneumoniae and other bacterial pathogens. Over 170 putative sRNAs have been identified in the S. pneumoniae TIGR4 strain (serotype 4) through transcriptomic studies, and a subset of these sRNAs has been further implicated in regulating pneumococcal pathogenesis. However, there is little overlap in the sRNAs identified among these studies, which indicates that the approaches used for sRNA identification were not sufficiently sensitive and robust and that there are likely many more undiscovered sRNAs encoded in the S. pneumoniae genome. Here, we sought to comprehensively identify sRNAs in Avery's virulent S. pneumoniae strain D39 using two independent RNA sequencing (RNA-seq)-based approaches. We developed an unbiased method for identifying novel sRNAs from bacterial RNA-seq data and have further tested the specificity of our analysis program toward identifying sRNAs encoded by both strains D39 and TIGR4. Interestingly, the genes for 15% of the putative sRNAs identified in strain TIGR4, including ones previously implicated in virulence, are not present in the strain D39 genome, suggesting that the differences in sRNA repertoires between these two serotypes may contribute to their strain-specific virulence properties. Finally, this study has identified 66 new sRNA candidates in strain D39, 30 of which have been further validated, raising the total number of sRNAs that have been identified in strain D39 to 112.IMPORTANCE Recent work has shown that sRNAs play crucial roles in S. pneumoniae pathogenesis, as inactivation of nearly one-third of the putative sRNA genes identified in one study led to reduced fitness or virulence in a murine model. Yet our understanding of sRNA-mediated gene regulation in S. pneumoniae has been hindered by limited knowledge about these regulatory RNAs, including which sRNAs are synthesized by different S. pneumoniae strains. We sought to address this problem by developing a sensitive sRNA detection technique to identify sRNAs in S. pneumoniae D39. A comparison of our data set reported here to those of other RNA-seq studies for S. pneumoniae strain D39 and TIGR4 has provided new insights into the S. pneumoniae sRNA transcriptome.

An evolving tale of two interacting RNAs-themes and variations of the T-box riboswitch mechanism

T-box riboswitches are a widespread class of structured noncoding RNAs in Gram-positive bacteria that regulate the expression of amino acid-related genes. They form negative feedback loops to maintain steady supplies of aminoacyl-transfer RNAs (tRNAs) to the translating ribosomes. T-box riboswitches are located in the 5' leader regions of mRNAs that they regulate and directly bind to their cognate tRNA ligands. T-boxes further sense the aminoacylation state of the bound tRNAs and, based on this readout, regulate gene expression at the level of transcription or translation. T-box riboswitches consist of two conserved domains-a 5' Stem I domain that is involved in specific tRNA recognition and a 3' antiterminator/antisequestrator (or discriminator) domain that senses the amino acid on the 3' end of the bound tRNA. Interaction of the 3' end of an uncharged but not charged tRNA with a thermodynamically weak discriminator domain stabilizes it to promote transcription readthrough or translation initiation. Recent biochemical, biophysical, and structural studies have provided high-resolution insights into the mechanism of tRNA recognition by Stem I, several structural models of full-length T-box-tRNA complexes, mechanism of amino acid sensing by the antiterminator domain, as well as kinetic details of tRNA binding to the T-box riboswitches. In addition, translation-regulating T-box riboswitches have been recently characterized, which presented key differences from the canonical transcriptional T-boxes. Here, we review the recent developments in understanding the T-box riboswitch mechanism that have employed various complementary approaches. Further, the regulation of multiple essential genes by T-boxes makes them very attractive drug targets to combat drug resistance. The recent progress in understanding the biochemical, structural, and dynamic aspects of the T-box riboswitch mechanism will enable more precise and effective targeting with small molecules. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1167-1180, 2019.

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.

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.

Comparative Transcriptomic Profiling of Yersinia enterocolitica O:3 and O:8 Reveals Major Expression Differences of Fitness- and Virulence-Relevant Genes Indicating Ecological Separation

Yersinia enterocolitica is a major diarrheal pathogen and is associated with a large range of gut-associated diseases. Members of this species have evolved into different phylogroups with genotypic variations. We performed the first characterization of the Y. enterocolitica transcriptional landscape and tracked the consequences of the genomic variations between two different pathogenic phylogroups by comparing their RNA repertoire, promoter usage, and expression profiles under four different virulence-relevant conditions. Our analysis revealed major differences in the transcriptional outputs of the closely related strains, pointing to an ecological separation in which one is more adapted to an environmental lifestyle and the other to a mostly mammal-associated lifestyle. Moreover, a variety of pathoadaptive alterations, including alterations in acid resistance genes, colonization factors, and toxins, were identified which affect virulence and host specificity. This illustrates that comparative transcriptomics is an excellent approach to discover differences in the functional output from closely related genomes affecting niche adaptation and virulence, which cannot be directly inferred from DNA sequences.

Yersinia enterocolitica is a zoonotic pathogen and an important cause of bacterial gastrointestinal infections in humans. Large-scale population genomic analyses revealed genetic and phenotypic diversity of this bacterial species, but little is known about the differences in the transcriptome organization, small RNA (sRNA) repertoire, and transcriptional output. Here, we present the first comparative high-resolution transcriptome analysis of Y. enterocolitica strains representing highly pathogenic phylogroup 2 (serotype O:8) and moderately pathogenic phylogroup 3 (serotype O:3) grown under four infection-relevant conditions. Our transcriptome sequencing (RNA-seq) approach revealed 1,299 and 1,076 transcriptional start sites and identified strain-specific sRNAs that could contribute to differential regulation among the phylogroups. Comparative transcriptomics further uncovered major gene expression differences, in particular, in the temperature-responsive regulon. Multiple virulence-relevant genes are differentially regulated between the two strains, supporting an ecological separation of phylogroups with certain niche-adapted properties. Strong upregulation of the ystA enterotoxin gene in combination with constitutive high expression of cell invasion factor InvA further showed that the toxicity of recent outbreak O:3 strains has increased. Overall, our report provides new insights into the specific transcriptome organization of phylogroups 2 and 3 and reveals gene expression differences contributing to the substantial phenotypic differences that exist between the lineages.

IMPORTANCEYersinia enterocolitica is a major diarrheal pathogen and is associated with a large range of gut-associated diseases. Members of this species have evolved into different phylogroups with genotypic variations. We performed the first characterization of the Y. enterocolitica transcriptional landscape and tracked the consequences of the genomic variations between two different pathogenic phylogroups by comparing their RNA repertoire, promoter usage, and expression profiles under four different virulence-relevant conditions. Our analysis revealed major differences in the transcriptional outputs of the closely related strains, pointing to an ecological separation in which one is more adapted to an environmental lifestyle and the other to a mostly mammal-associated lifestyle. Moreover, a variety of pathoadaptive alterations, including alterations in acid resistance genes, colonization factors, and toxins, were identified which affect virulence and host specificity. This illustrates that comparative transcriptomics is an excellent approach to discover differences in the functional output from closely related genomes affecting niche adaptation and virulence, which cannot be directly inferred from DNA sequences.

Discovery of Small-Molecule Antibiotics against a Unique tRNA-Mediated Regulation of Transcription in Gram-Positive Bacteria

The emergence of multidrug-resistant bacteria necessitates the identification of unique targets of intervention and compounds that inhibit their function. Gram-positive bacteria use a well-conserved tRNA-responsive transcriptional regulatory element in mRNAs, known as the T-box, to regulate the transcription of multiple operons that control amino acid metabolism. T-box regulatory elements are found only in the 5'-untranslated region (UTR) of mRNAs of Gram-positive bacteria, not Gram-negative bacteria or the human host. Using the structure of the 5'UTR sequence of the Bacillus subtilis tyrosyl-tRNA synthetase mRNA T-box as a model, in silico docking of 305 000 small compounds initially yielded 700 as potential binders that could inhibit the binding of the tRNA ligand. A single family of compounds inhibited the growth of Gram-positive bacteria, but not Gram-negative bacteria, including drug-resistant clinical isolates at minimum inhibitory concentrations (MIC 16-64 μg mL^-1 ). Resistance developed at an extremely low mutational frequency (1.21×10^-10 ). At 4 μg mL^-1 , the parent compound PKZ18 significantly inhibited in vivo transcription of glycyl-tRNA synthetase mRNA. PKZ18 also inhibited in vivo translation of the S. aureus threonyl-tRNA synthetase protein. PKZ18 bound to the Specifier Loop in vitro (K_d ≈24 μm). Its core chemistry necessary for antibacterial activity has been identified. These findings support the T-box regulatory mechanism as a new target for antibiotic discovery that may impede the emergence of resistance.