Functional roles of a tetraloop/receptor interacting module in a cyclic di-GMP riboswitch

Riboswitches, from cognition to transformation

Riboswitches are functional RNA elements that regulate gene expression by directly detecting metabolites. Twenty years have passed since it was first discovered, researches on riboswitches are becoming increasingly standardized and refined, which could significantly promote people's cognition of RNA function as well. Here, we focus on some representative orphan riboswitches, enumerate the structural and functional transformation and artificial design of riboswitches including the coupling with ribozymes, hoping to attain a comprehensive understanding of riboswitch research.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.

The Theophylline Aptamer: 25 Years as an Important Tool in Cellular Engineering Research

The theophylline aptamer was isolated from an oligonucleotide library in 1994. Since that time, the aptamer has found wide utility, particularly in synthetic biology, cellular engineering, and diagnostic applications. The primary application of the theophylline aptamer is in the construction and characterization of synthetic riboswitches for regulation of gene expression. These riboswitches have been used to control cellular motility, regulate carbon metabolism, construct logic gates, screen for mutant enzymes, and control apoptosis. Other applications of the theophylline aptamer in cellular engineering include regulation of RNA interference and genome editing through CRISPR systems. Here we describe the uses of the theophylline aptamer for cellular engineering over the past 25 years. In so doing, we also highlight important synthetic biology applications to control gene expression in a ligand-dependent manner.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

Recognition of cyclic-di-GMP by a riboswitch conducts translational repression through masking the ribosome-binding site distant from the aptamer domain

The riboswitch is a class of RNA-based gene regulatory machinery that is dependent on recognition of its target ligand by RNA tertiary structures. Ligand recognition is achieved by the aptamer domain, and ligand-dependent structural changes of the expression platform then usually mediate termination of transcription or translational initiation. Ligand-dependent structural changes of the aptamer domain and expression platform have been reported for several riboswitches with short (<40 nucleotides) expression platforms. In this study, we characterized structural changes of the Vc2 c-di-GMP riboswitch that represses translation of downstream open reading frames in a ligand-dependent manner. The Vc2 riboswitch has a long (97 nucleotides) expression platform, but its structure and function are largely unknown. Through mutational analysis and chemical probing, we identified its secondary structures that are possibly responsible for switch-OFF and switch-ON states of translational initiation.

Cyclic di-GMP Regulates TfoY in Vibrio cholerae To Control Motility by both Transcriptional and Posttranscriptional Mechanisms

3',5'-Cyclic diguanylic acid (c-di-GMP) is a bacterial second messenger molecule that is a key global regulator in Vibrio cholerae, but the molecular mechanisms by which this molecule regulates downstream phenotypes have not been fully characterized. One such regulatory factor that may respond to c-di-GMP is the Vc2 c-di-GMP-binding riboswitch that is hypothesized to control the expression of the downstream putative transcription factor TfoY. Although much is known about the physical and structural properties of the Vc2 riboswitch aptamer, the nature of its expression and function in V. cholerae has not been investigated. Here, we show that Vc2 functions as an off switch to inhibit TfoY production at intermediate and high concentrations of c-di-GMP. At low c-di-GMP concentrations, TfoY production is induced to stimulate dispersive motility. We also observed increased transcription of tfoY at high intracellular concentrations of c-di-GMP, but this induction is independent of the Vc2 riboswitch and occurs via transcriptional control of promoters upstream of tfoY by the previously identified c-di-GMP dependent transcription factor VpsR. Our results show that TfoY is induced by c-di-GMP at both low and high intracellular concentrations of c-di-GMP via posttranscriptional and transcriptional mechanisms, respectively. This regulation contributes to the formation of three distinct c-di-GMP signaling states in V. choleraeIMPORTANCE The bacterial pathogen Vibrio cholerae must transition between life in aquatic environmental reservoirs and life in the gastrointestinal tract. Biofilm formation and bacterial motility, and their control by the second messenger molecule c-di-GMP, play integral roles in this adaptation. Here, we define the third major mechanism by which c-di-GMP controls bacterial motility. This pathway utilizes a noncoding RNA element known as a riboswitch that, when bound to c-di-GMP, inhibits the expression of the transcription factor TfoY. TfoY production switches V. cholerae motility from a dense to a dispersive state. Our results suggest that the c-di-GMP signaling network of V. cholerae can exist in at least three distinct states to regulate biofilm formation and motility.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Insights into Stability and Folding of GNRA and UNCG Tetraloops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15-μs-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the GL4 nucleobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that GL4-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.

Mutational analysis of structural elements in a class-I cyclic di-GMP riboswitch to elucidate its regulatory mechanism

The Vc2 riboswitch possesses an aptamer domain belonging to the class-I c-di-GMP riboswitch family. This domain has been analysed and the molecular mechanism by which it recognizes the c-di-GMP ligand has been elucidated. On the other hand, the regulatory mechanism of the full-length Vc2 riboswitch to control its downstream open reading frame (ORF) remains largely unknown. In this study, we performed in vivo reporter assays and in vitro biochemical analyses of the full-length riboswitch and its aptamer domain. We evaluated the results of in vivo and in vitro analyses to elucidate the regulatory mechanism of the Vc2 riboswitch. The present results suggest that recognition of c-di-GMP ligand by the Vc2 riboswitch aptamer domain downregulates expression of its downstream ORF primarily at the translational level.

Characterization of an RNA receptor motif that recognizes a GCGA tetraloop

Tertiary interactions between a new RNA motif and RNA tetraloops were analyzed to determine whether this new motif shows preference for a GCGA tetraloop. In the structural context of a ligase ribozyme, this motif discriminated GCGA loop from 3 other tetraloops. The affinity between the GCGA loop and its receptor is strong enough to carry out the ribozyme activity.

Optimization of RNA-based c-di-GMP fluorescent sensors through tuning their structural modules

Cyclic diguanylate (c-di-GMP) is a second messenger of bacteria and its detection is an important issue in basic and applied microbiology. As c-di-GMP riboswitch ligand-binding domains (aptamer domains) capture c-di-GMP with high affinity and selectivity, they are promising platforms for the development of RNA-based c-di-GMP sensors. We analyzed two previously reported c-di-GMP sensor RNAs derived from the Vc2 riboswitch. We also designed and tested their variants, some of which showed improved properties as RNA-based c-di-GMP sensors.

Antimicrobial aptamers for detection and inhibition of microbial pathogen growth

Discovery of alternative sources of antimicrobial agents are essential in the ongoing battle against microbial pathogens. Legislative and scientific challenges considerably hinder the discovery and use of new antimicrobial drugs, and new approaches are in urgent demand. On the other hand, rapid, specific and sensitive detection of airborne pathogens is becoming increasingly critical for public health. In this respect affinity oligonucleotides, aptamers, provide unique opportunities for the development of nanotechnological solutions for such medical applications. In recent years, aptamers specifically recognizing microbial cells and viruses showed great potential in a range of analytical and therapeutic applications. This article describes the significant advances in the development of aptamers targeting specific pathogens. Therapeutic application of aptamers as neutralizing agents demonstrates great potential as a future source of antimicrobial agent.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Identification of a candidate single-nucleotide polymorphism related to chemotherapeutic response through a combination of knowledge-based algorithm and hypothesis-free genomic data

Inter-individual variations in drug responses among patients are known to cause serious problems in medicine. Genome-wide association study (GWAS) is powerful for examining single-nucleotide polymorphisms (SNPs) and their relationships with drug response variations. However, no significant SNP has been identified using GWAS due to multiple testing problems. Therefore, we propose a combination method consisting of knowledge-based algorithm, two stages of screening, and permutation test for identifying SNPs in the present study. We applied this method to a genome-wide pharmacogenomics study for which 109,365 SNPs had been genotyped using Illumina Human-1 BeadChip for 119 gastric cancer patients treated with fluoropyrimidine. We identified rs2293347 in epidermal growth factor receptor (EGFR) is as a candidate SNP related to chemotherapeutic response. The p value for the rs2293347 was 2.19 × 10(-5) for Fisher's exact test, and the p value was 0.00360 for the permutation test (multiple testing problems are corrected). Additionally, rs2293347 was clearly superior to clinical parameters and showed a sensitivity value of 55.0% and specificity value of 94.4% in the evaluation by using multiple regression models. Recent studies have shown that combination chemotherapy of fluoropyrimidine and EGFR-targeting agents is effective for gastric cancer patients highly expressing EGFR. These results suggest that rs2293347 is a potential predictive factor for selecting chemotherapies, such as fluoropyrimidine alone or combination chemotherapies.