Flipping off the riboswitch: RNA structures that control gene expression

Pharmacological perturbation of thiamine metabolism sensitizes Pseudomonas aeruginosa to multiple antibacterial agents

New therapeutic concepts are critically needed for carbapenem-resistant Pseudomonas aeruginosa, an opportunistic pathogen particularly recalcitrant to antibiotics. The screening of around 230,000 small molecules yielded a very low hit rate of 0.002% after triaging for known antibiotics. The only novel hit that stood out was the antimetabolite oxythiamine. Oxythiamine is a known transketolase inhibitor in eukaryotic cells, but its antibacterial potency has not been reported. Metabolic and transcriptomic analyses indicated that oxythiamine is intracellularly converted to oxythiamine pyrophosphate and subsequently inhibits several vitamin-B1-dependent enzymes, sensitizing the bacteria to several antibiotic and non-antibiotic drugs such as tetracyclines, 5-fluorouracil, and auranofin. The positive interaction between 5-fluorouracil and oxythiamine was confirmed in a murine ocular infection model, indicating relevance during infection. Together, this study revealed a system-level significance of thiamine metabolism perturbation that sensitizes P. aeruginosa to multiple small molecules, a property that could inform on the development of a rational drug combination.

Efficient Method to Identify Synthetic Riboswitches Using RNA-Based Capture-SELEX Combined with In Vivo Screening

Synthetic riboswitches are a promising tool for conditional gene expression. In vitro selected aptamers used as binding domains for the design of RNA-based switches have to exhibit excellent binding affinity as well as ligand binding-induced structural changes. Selection via Capture-SELEX favors the enrichment of aptamers which exhibit both characteristics. For the Capture-SELEX, an RNA pool is used that gets immobilized onto a capture oligonucleotide by hybridization. Addition of the ligand frees the aptamers by their binding to the ligand, resulting in the release from the capture oligonucleotide through structural changes. These sequences get reverse transcribed, PCR amplified, and used for the following selection rounds. In this publication, we present a detailed protocol for Capture-SELEX, followed by screening in yeast to identify aptamers suitable for the design of synthetic riboswitches.

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.

Towards an understanding of RNA structural modalities: a riboswitch case study

A riboswitch is a type of RNA molecule that regulates important biological functions by changing structure, typically under ligand-binding. We assess the extent that these ligand-bound structural alternatives are present in the Boltzmann sample, a standard RNA secondary structure prediction method, for three riboswitch test cases. We use the cluster analysis tool RNAStructProfiling to characterize the different modalities present among the suboptimal structures sampled. We compare these modalities to the putative base pairing models obtained from independent experiments using NMR or fluorescence spectroscopy. We find, somewhat unexpectedly, that profiling the Boltzmann sample captures evidence of ligand-bound conformations for two of three riboswitches studied. Moreover, this agreement between predicted modalities and experimental models is consistent with the classification of riboswitches into thermodynamic versus kinetic regulatory mechanisms. Our results support cluster analysis of Boltzmann samples by RNAStructProfiling as a possible basis for de novo identification of thermodynamic riboswitches, while highlighting the challenges for kinetic ones.

A 5'-cap for DNA probes binding RNA target strands

Detecting short RNA strands with high fidelity at any of the bases of their sequence, including the termini, can be challenging, since fraying, wobbling, and refolding all compete with canonical base pairing. We performed a search for 5'-substituents of oligodeoxynucleotides that increase base pairing fidelity at the terminus of duplexes with RNA target strands. From a total of over 70 caps, differing in stacking moiety and linker, a phosphodiester-linked sequence of the residues of L-prolinol, glycine, and oxolinic acid, dubbed ogOA, was identified as a 5'-cap that stabilizes any of the four canonical base pairs, with ΔT(m) values of up to +13.1 °C for an octamer. At the same time, the cap increases discrimination against any of the 12 possible terminal mismatches, including mismatches that are more stable than their perfectly matched counterparts in the control duplex, such as A:A. A probe with the cap also showed increased selectivity in the detection of two closely related microRNAs, let7c and let7a, with a ΔT(m) value of 9.2 °C. Melting curves also yielded thermodynamic data that shed light on the uniformity of molecular recognition in the sequence space of DNA:DNA and DNA:RNA duplexes. Hybridization probes with fidelity-enhancing caps should find applications in the individual and parallel detection of biologically active RNA species.

SAM recognition and conformational switching mechanism in the Bacillus subtilis yitJ S box/SAM-I riboswitch

S-box (SAM-I) riboswitches are a widespread class of riboswitches involved in the regulation of sulfur metabolism in Gram-positive bacteria. We report here the 3.0-Å crystal structure of the aptamer domain of the Bacillus subtilis yitJ S-box (SAM-I) riboswitch bound to S-adenosyl-L-methionine (SAM). The RNA folds into two sets of helical stacks spatially arranged by tertiary interactions including a K-turn and a pseudoknot at a four-way junction. The tertiary structure is further stabilized by metal coordination, extensive ribose zipper interactions, and SAM-mediated tertiary interactions. Despite structural differences in the peripheral regions, the SAM-binding core of the B. subtilis yitJ riboswitch is virtually superimposable with the previously determined Thermoanaerobacter tengcongensis yitJ riboswitch structure, suggesting that a highly conserved ligand-recognition mechanism is utilized by all S-box riboswitches. SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemical probing analysis further revealed that the alternative base-pairing element in the expression platform controls the conformational switching process. In the absence of SAM, the apo yitJ aptamer domain folds predominantly into a pre-binding conformation that resembles, but is not identical with, the SAM-bound state. We propose that SAM enters the ligand-binding site through the "J1/2-J3/4" gate and "locks" down the SAM-bound conformation through an induced-fit mechanism. Temperature-dependent SHAPE revealed that the tertiary interaction-stabilized SAM-binding core is extremely stable, likely due to the cooperative RNA folding behavior. Mutational studies revealed that certain modifications in the SAM-binding region result in loss of SAM binding and constitutive termination, which suggests that these mutations lock the RNA into a form that resembles the SAM-bound form in the absence of SAM.

High fidelity base pairing at the 3'-terminus

Binding target strands with single base selectivity at a terminal position is difficult with natural DNA or RNA hybridization probes. Nature uses a degenerate genetic code that is based on RNA:RNA codon:anticodon duplexes tolerating wobble base pairs at the terminus. The importance of short RNA strands in regulatory processes in the cell make it desirable to develop receptor-like approaches for high fidelity binding, even at the very 3'-terminus of a probe. Here, we report the three-dimensional structure of a DNA duplex with a 3'-terminal 2'-anthraquinoylamido-2'-deoxyuridine (Uaq) residue that was solved by NMR and restrained molecular dynamics. The Uaq residue binds the 5'-terminus of the target strand through a combination of pi-stacking, hydrogen bonding, and interactions in the minor groove. The acylated aminonucleoside is the best molecular cap for 3'-termini reported to date. The Uaq motif assists binding of DNA strands, but is particularly effective in enhancing the affinity for RNA target strands, with a DeltaT(m) in the UV melting point of up to +18.2 degrees C per residue. Increased base pairing selectivity is induced for all sequence motifs tested, even in cases where unmodified duplexes show no preference for the canonical base pair at all. A single mismatched nucleobase facing the 3'-terminus gives DeltaDeltaT(m) values as large as -23.9 degrees C (RNA) or -29.5 degrees C (DNA). The 5'-phosphoramidite of the Uaq cap reported here allows for routine incorporation during automated syntheses.

Riboswitches: ancient and promising genetic regulators

BAIT AND SWITCH: Metabolite-sensing riboswitches make use of RNA structural modulation to regulate gene expression, as illustrated in the scheme, in response to subtle changes in metabolite concentrations. This review describes the current knowledge about naturally occurring riboswitches and their growing potential as antibacterial cellular targets and as molecular biosensors. Newly discovered metabolite-sensing riboswitches have revealed that cellular processes extensively make use of RNA structural modulation to regulate gene expression in response to subtle changes in metabolite concentrations. Riboswitches are involved at various regulation levels of gene expression, such as transcription attenuation, translation initiation, mRNA splicing and mRNA processing. Riboswitches are found in the three kingdoms of life, and in various cases, are involved in the regulation of essential genes, which makes their regulation an essential part of cell survival. Because riboswitches operate without the assistance of accessory proteins, they are believed to be remnants of an ancient time, when gene regulation was strictly based on RNA, from which are left numerous "living molecular fossils", as exemplified by ribozymes, and more spectacularly, by the ribosome. Due to their nature, riboswitches hold high expectations for the manipulation of gene expression and the detection of small metabolites, and also offer an unprecedented potential for the discovery of novel classes of antimicrobial agents.

Riboswitches: emerging themes in RNA structure and function

Riboswitches are RNAs capable of binding cellular metabolites using a diverse array of secondary and tertiary structures to modulate gene expression. The recent determination of the three-dimensional structures of parts of six different riboswitches illuminates common features that allow riboswitches to be grouped into one of two types. Type I riboswitches, as exemplified by the purine riboswitch, are characterized by a single, localized binding pocket supported by a largely pre-established global fold. This arrangement limits ligand-induced conformational changes in the RNA to a small region. In contrast, Type II riboswitches, such as the thiamine pyrophosphate riboswitch, contain binding pockets split into at least two spatially distinct sites. As a result, binding induces both local changes to the binding pocket and global architecture. Similar organizational themes are found in other noncoding RNAs, making it possible to begin to build a hierarchical classification of RNA structure based on the spatial organization of their active sites and associated secondary structural elements.

Identification of metabolites with anticancer properties by computational metabolomics

Background

Certain endogenous metabolites can influence the rate of cancer cell growth. For example, diacylglycerol, ceramides and sphingosine, NAD⁺ and arginine exert this effect by acting as signaling molecules, while carrying out other important cellular functions. Metabolites can also be involved in the control of cell proliferation by directly regulating gene expression in ways that are signaling pathway-independent, e.g. by direct activation of transcription factors or by inducing epigenetic processes. The fact that metabolites can affect the cancer process on so many levels suggests that the change in concentration of some metabolites that occurs in cancer cells could have an active role in the progress of the disease.

Results

CoMet, a fully automated Computational Metabolomics method to predict changes in metabolite levels in cancer cells compared to normal references has been developed and applied to Jurkat T leukemia cells with the goal of testing the following hypothesis: Up or down regulation in cancer cells of the expression of genes encoding for metabolic enzymes leads to changes in intracellular metabolite concentrations that contribute to disease progression. All nine metabolites predicted to be lowered in Jurkat cells with respect to lymphoblasts that were examined (riboflavin, tryptamine, 3-sulfino-L-alanine, menaquinone, dehydroepiandrosterone, α-hydroxystearic acid, hydroxyacetone, seleno-L-methionine and 5,6-dimethylbenzimidazole), exhibited antiproliferative activity that has not been reported before, while only two (bilirubin and androsterone) of the eleven tested metabolites predicted to be increased or unchanged in Jurkat cells displayed significant antiproliferative activity.

Conclusion

These results: a) demonstrate that CoMet is a valuable method to identify potential compounds for experimental validation, b) indicate that cancer cell metabolism may be regulated to reduce the intracellular concentration of certain antiproliferative metabolites, leading to uninhibited cellular growth and c) suggest that many other endogenous metabolites with important roles in carcinogenesis are awaiting discovery.

RNA switches regulate initiation of translation in bacteria

A large variety of RNA-based mechanisms have been uncovered in all living organisms to regulate gene expression in response to internal and external changes, and to rapidly adapt cell growth in response to these signals. In bacteria, structural elements in the 5′ leader regions of mRNAs have direct effects on translation initiation of the downstream coding sequences. The docking and unfolding of these mRNAs on the 30S subunit are critical steps in the initiation process directly modulating and timing translation. Structural elements can also undergo conformational changes in response to environmental cues (i.e., temperature sensors) or upon binding of a variety oftrans-acting factors, such as metabolites, non-coding RNAs or regulatory proteins. These RNA switches can temporally regulate translation, leading either to repression or to activation of protein synthesis.

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.

Ligand-induced folding of the thiM TPP riboswitch investigated by a structure-based fluorescence spectroscopic approach

Riboswitches are genetic control elements within non-coding regions of mRNA. They consist of a metabolite-sensitive aptamer and an adjoining expression platform. Here, we describe ligand-induced folding of a thiamine pyrophosphate (TPP) responsive riboswitch from Escherichia coli thiM mRNA, using chemically labeled variants. Referring to a recent structure determination of the TPP/aptamer complex, each variant was synthesized with a single 2-aminopurine (AP) nucleobase replacement that was selected to monitor formation of tertiary interactions of a particular region during ligand binding in real time by fluorescence experiments. We have determined the rate constants for conformational adjustment of the individual AP sensors. From the 7-fold differentiation of these constants, it can be deduced that tertiary contacts between the two parallel helical domains (P2/J3-2/P3/L3 and P4/P5/L5) that grip the ligand's ends in two separate pockets, form significantly faster than the function-critical three-way junction with stem P1 fully developed. Based on these data, we characterize the process of ligand binding by an induced fit of the RNA and propose a folding model of the TPP riboswitch aptamer. For the full-length riboswitch domain and for shorter constructs that represent transcriptional intermediates, we have additionally evaluated ligand-induced folding via AP-modified variants and provide insights into the sequential folding pathway that involves a finely balanced equilibrium of secondary structures.

Ligand-induced folding of the adenosine deaminase A-riboswitch and implications on riboswitch translational control

By using a structure-based fluorescence spectroscopic approach, we have examined the folding of an adenine-responsive riboswitch that regulates translation initiation. We observed adaptive recognition of the ligand for the aptamer domain of adenosine deaminase (add) mRNA from Vibrio vulnificus, and revealed pronounced conformational changes even in the preorganized loop-loop region that is distant from the binding site. Importantly, the full-length riboswitch domain, which has a potential translational repressor stem is able to form a binary complex with adenine, and does not act as a folding trap to inhibit binding. The aptamer that is extended by the expression platform therefore remains fully responsive to its ligand; this is in contrast to the previously investigated pbuE A-riboswitch, which becomes trapped in a nonresponsive terminator fold. Consequently, the latter must employ complex response mechanisms, such as operating in kinetic-control mode or using transcriptional pausing, to provide time for the aptamer portion to fold and to bind. The different behavior of the riboswitches can be rationalized by their distinct sequence interface between the aptamer and expression platform. For the add A-riboswitch, our data suggest a thermodynamically driven response mechanism.