Phylogenetic analysis and comparative genomics of purine riboswitch distribution in prokaryotes

A journey into the regulatory secrets of the de novo purine nucleotide biosynthesis

De novo purine nucleotide biosynthesis (DNPNB) consists of sequential reactions that are majorly conserved in living organisms. Several regulation events take place to maintain physiological concentrations of adenylate and guanylate nucleotides in cells and to fine-tune the production of purine nucleotides in response to changing cellular demands. Recent years have seen a renewed interest in the DNPNB enzymes, with some being highlighted as promising targets for therapeutic molecules. Herein, a review of two newly revealed modes of regulation of the DNPNB pathway has been carried out: i) the unprecedent allosteric regulation of one of the limiting enzymes of the pathway named inosine 5'-monophosphate dehydrogenase (IMPDH), and ii) the supramolecular assembly of DNPNB enzymes. Moreover, recent advances that revealed the therapeutic potential of DNPNB enzymes in bacteria could open the road for the pharmacological development of novel antibiotics.

Phylogenomic and comparative analysis of the distribution and regulatory patterns of TPP riboswitches in fungi

Riboswitches are metabolite or ion sensing cis-regulatory elements that regulate the expression of the associated genes involved in biosynthesis or transport of the corresponding metabolite. Among the nearly 40 different classes of riboswitches discovered in bacteria so far, only the TPP riboswitch has also been found in algae, plants, and in fungi where their presence has been experimentally validated in a few instances. We analyzed all the available complete fungal and related genomes and identified TPP riboswitch-based regulation systems in 138 fungi and 15 oomycetes. We find that TPP riboswitches are most abundant in Ascomycota and Basidiomycota where they regulate TPP biosynthesis and/or transporter genes. Many of these transporter genes were found to contain conserved domains consistent with nucleoside, urea and amino acid transporter gene families. The genomic location of TPP riboswitches when correlated with the intron structure of the regulated genes enabled prediction of the precise regulation mechanism employed by each riboswitch. Our comprehensive analysis of TPP riboswitches in fungi provides insights about the phylogenomic distribution, regulatory patterns and functioning mechanisms of TPP riboswitches across diverse fungal species and provides a useful resource that will enhance the understanding of RNA-based gene regulation in eukaryotes.

Comparative genomics and phylogenomic analyses of lysine riboswitch distributions in bacteria

Riboswitches are cis-regulatory elements that regulate the expression of genes involved in biosynthesis or transport of a ligand that binds to them. Among the nearly 40 classes of riboswitches discovered so far, three are known to regulate the concentration of biologically encoded amino acids glycine, lysine, and glutamine. While some comparative genomics studies of riboswitches focusing on their gross distribution across different bacterial taxa have been carried out recently, systematic functional annotation and analysis of lysine riboswitches and the genes they regulate are still lacking. We analyzed 2785 complete bacterial genome sequences to systematically identify 468 lysine riboswitches (not counting hits from multiple strains of the same species) and obtain a detailed phylogenomic map of gene-specific lysine riboswitch distribution across diverse prokaryotic phyla. We find that lysine riboswitches are most abundant in Firmicutes and Gammaproteobacteria where they are found upstream to both biosynthesis and/or transporter genes. They are relatively rare in all other prokaryotic phyla where if present they are primarily found upstream to operons containing many lysine biosynthesis genes. The genome-wide study of the genetic organisation of the lysine riboswitches show considerable variation both within and across different Firmicute orders. Correlating the location of a riboswitch with its genomic context and its phylogenetic relationship with other evolutionarily related riboswitch carrying species, enables identification and annotation of many lysine biosynthesis, transporter and catabolic genes. It also reveals previously unknown patterns of lysine riboswitch distribution and gene/operon regulation and allows us to draw inferences about the possible point of origin of lysine riboswitches. Additionally, evidence of horizontal transfer of riboswitches was found between Firmicutes and Actinobacteria. Our analysis provides a useful resource that will lead to a better understanding of the evolution of these regulatory elements and prove to be beneficial for exploiting riboswitches for developing targeted therapies.

Development of a new sequential block finding strategy for detection of conserved sequences in riboswitches

Introduction: Some non-coding RNAs have an important role in the regulation of gene expression and consequently cellular function. Riboswitches are examples of these regulatory RNAs. Riboswitches are classified into various families according to sequential and structural similarities.

Methods: In this study, a block finder algorithm for identification of frequently appearing sequential blocks in five families of riboswitches from Rfam 12.0 database, without the use of alignment methods, was developed.

Results: The developed program identified 21 frequently appearing blocks in five families of riboswitches.

Conclusion: Comparison of the results of the proposed algorithm with those of sequential alignment methods revealed that our method can recognize most of the patterns present in conserved areas of individual riboswitch families and determine them as specific blocks, implying potential of the developed program as a platform for further studies and developments.

Gene expression control by Bacillus anthracis purine riboswitches

In all kingdoms of life, cellular replication relies on the presence of nucleosides and nucleotides, the building blocks of nucleic acids and the main source of energy. In bacteria, the availability of metabolites sometimes directly regulates the expression of enzymes and proteins involved in purine salvage, biosynthesis, and uptake through riboswitches. Riboswitches are located in bacterial mRNAs and can control gene expression by conformational changes in response to ligand binding. We have established an inverse reporter gene system in Bacillus subtilis that allows us to monitor riboswitch-controlled gene expression. We used it to investigate the activity of five potential purine riboswitches from Bacillus anthracis in response to different purines and pyrimidines. Furthermore, in vitro studies on the aptamer domains of the riboswitches reveal their variation in guanine binding affinity ranging from namomolar to micromolar. These data do not only provide insight into metabolite sensing but can also aid in engineering artificial cell regulatory systems.

TPP riboswitch characterization in Alishewanella tabrizica and Alishewanella aestuarii and comparison with other TPP riboswitches

Riboswitches are located in non-coding areas of mRNAs and act as sensors of cellular small molecules, regulating gene expression in response to ligand binding. The TPP riboswitch is the most widespread riboswitch occurring in all three domains of life. However, it has been rarely characterized in environmental bacteria other than Escherichia coli and Bacillus subtilis. In this study, TPP riboswitches located in the 5' UTR of thiC operon from Alishewanella tabrizica and Alishewanella aestuarii were identified and characterized. Moreover, affinity analysis of TPP binding to the TPP aptamer domains originated from A. tabrizica, A. aestuarii, E.coli, and B. subtilis were studied and compared using In-line probing and Surface Plasmon Resonance (SPR). TPP binding to the studied RNAs from A. tabrizica and A. aestuarii caused distinctive changes of the In-line cleavage pattern, demonstrating them as functional TPP riboswitches. With dissociation constant of 2-4nM (depending on the method utilized), the affinity of TPP binding was highest in A. tabrizica, followed by the motifs sourced from A. aestuarii, E. coli, and B. subtilis. The observed variation in their TPP-binding affinity might be associated with adaptation to the different environments of the studied bacteria.

De Novo Guanine Biosynthesis but Not the Riboswitch-Regulated Purine Salvage Pathway Is Required for Staphylococcus aureus Infection In Vivo

De novo guanine biosynthesis is an evolutionarily conserved pathway that creates sufficient nucleotides to support DNA replication, transcription, and translation. Bacteria can also salvage nutrients from the environment to supplement the de novo pathway, but the relative importance of either pathway during Staphylococcus aureus infection is not known. In S. aureus, genes important for both de novo and salvage pathways are regulated by a guanine riboswitch. Bacterial riboswitches have attracted attention as a novel class of antibacterial drug targets because they have high affinity for small molecules, are absent in humans, and regulate the expression of multiple genes, including those essential for cell viability. Genetic and biophysical methods confirm the existence of a bona fide guanine riboswitch upstream of an operon encoding xanthine phosphoribosyltransferase (xpt), xanthine permease (pbuX), inosine-5'-monophosphate dehydrogenase (guaB), and GMP synthetase (guaA) that represses the expression of these genes in response to guanine. We found that S. aureus guaB and guaA are also transcribed independently of riboswitch control by alternative promoter elements. Deletion of xpt-pbuX-guaB-guaA genes resulted in guanine auxotrophy, failure to grow in human serum, profound abnormalities in cell morphology, and avirulence in mouse infection models, whereas deletion of the purine salvage genes xpt-pbuX had none of these effects. Disruption of guaB or guaA recapitulates the xpt-pbuX-guaB-guaA deletion in vivo In total, the data demonstrate that targeting the guanine riboswitch alone is insufficient to treat S. aureus infections but that inhibition of guaA or guaB could have therapeutic utility.

IMPORTANCE

De novo guanine biosynthesis and purine salvage genes were reported to be regulated by a guanine riboswitch in Staphylococcus aureus We demonstrate here that this is not true, because alternative promoter elements that uncouple the de novo pathway from riboswitch regulation were identified. We found that in animal models of infection, the purine salvage pathway is insufficient for S. aureus survival in the absence of de novo guanine biosynthesis. These data suggest targeting the de novo guanine biosynthesis pathway may have therapeutic utility in the treatment of S. aureus infections.

An Efficient Minimum Free Energy Structure-Based Search Method for Riboswitch Identification Based on Inverse RNA Folding

Riboswitches are RNA genetic control elements that were originally discovered in bacteria and provide a unique mechanism of gene regulation. They work without the participation of proteins and are believed to represent ancient regulatory systems in the evolutionary timescale. One of the biggest challenges in riboswitch research is to find additional eukaryotic riboswitches since more than 20 riboswitch classes have been found in prokaryotes but only one class has been found in eukaryotes. Moreover, this single known class of eukaryotic riboswitch, namely the TPP riboswitch class, has been found in bacteria, archaea, fungi and plants but not in animals. The few examples of eukaryotic riboswitches were identified using sequence-based bioinformatics search methods such as a combination of BLAST and pattern matching techniques that incorporate base-pairing considerations. None of these approaches perform energy minimization structure predictions. There is a clear motivation to develop new bioinformatics methods, aside of the ongoing advances in covariance models, that will sample the sequence search space more flexibly using structural guidance while retaining the computational efficiency of sequence-based methods. We present a new energy minimization approach that transforms structure-based search into a sequence-based search, thereby enabling the utilization of well established sequence-based search utilities such as BLAST and FASTA. The transformation to sequence space is obtained by using an extended inverse RNA folding problem solver with sequence and structure constraints, available within RNAfbinv. Examples in applying the new method are presented for the purine and preQ1 riboswitches. The method is described in detail along with its findings in prokaryotes. Potential uses in finding novel eukaryotic riboswitches and optimizing pre-designed synthetic riboswitches based on ligand simulations are discussed. The method components are freely available for use.