Structure-guided design of a high-affinity ligand for a riboswitch

1H, 13C, 15N and 31P chemical shift assignment of the first stem-loop Guanidine-II riboswitch from Escherichia coli

A comprehensive understanding of RNA-based gene regulation is a fundamental aspect for the development of innovative therapeutic options in medicine and for a more targeted response to environmental problems. Within the different mechanisms of RNA-based gene regulation, riboswitches are particularly interesting as they change their structure in response to the interaction with a low molecular weight ligand, often a well-known metabolite. Four distinct classes of riboswitches recognize the very small guanidinium cation. We are focused on the Guanidine-II riboswitch with the mini-ykkC motif. We report here the assignment of the 1H, 13C, 15N and 31P chemical shifts of the 23 nucleotide-long sequence of the first stem-loop of the Guanidine-II riboswitch aptamer from Escherichia coli. Despite its small size, the assignment of the NMR signals of this RNA proved to be challenging as it has symmetrical base pairs and palindromic character.

Cooperative binding of bivalent ligands yields new insights into the guanidine-II riboswitch

Riboswitches are involved in regulating the gene expression in bacteria. They are located within the untranslated regions of bacterial messenger RNA and function as switches by adjusting their shape, depending on the presence or absence of specific ligands. To decipher the fundamental aspects of bacterial gene control, it is therefore important to understand the mechanisms that underlie these conformational switches. To this end, a combination of an experimental binding study, molecular simulations and machine learning has been employed to obtain insights into the conformational changes and structural dynamics of the guanidine-II riboswitch. By exploiting the design of a bivalent ligand, we were able to study ligand binding in the aptamer dimer at the molecular level. Spontaneous ligand-binding events, which are usually difficult to simulate, were observed and the contributing factors are described. These findings were further confirmed by in vivo experiments, where the cooperative binding effects of the bivalent ligands resulted in increased binding affinity compared to the native guanidinium ligand. Beyond ligand binding itself, the simulations revealed a novel, ligand-dependent base-stacking interaction outside of the binding pocket that stabilizes the riboswitch.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.

Tying the knot in the tetrahydrofolate (THF) riboswitch: A molecular basis for gene regulation

Effective gene regulation by the tetrahydrofolate riboswitch depends not only on ligand affinity but also on the kinetics of ligand association, which involves two cooperative binding sites. We have determined a 1.9-Å resolution crystal structure of the ligand-free THF riboswitch aptamer. The pseudoknot binding site 'unwinds' in the absence of ligand, whereby the adjacent helical domains (P1, P2, and P3) become disjointed, resulting in rotation and misalignment of the gene-regulatory P1 helix with respect to P3. In contrast, the second binding site at the three-way junction, which is the first to fold, is structurally conserved between apo and holo forms. This suggests a kinetic role for this site, in which binding of the first ligand molecule to the stably folded three-way junction promotes formation of the regulatory pseudoknot site and subsequent binding of the second molecule. As such, these findings provide a molecular basis for both conformational switching and kinetic control.

Biochemical Validation of a Fourth Guanidine Riboswitch Class in Bacteria

An intriguing consequence of ongoing riboswitch discovery efforts is the occasional identification of metabolic or toxicity response pathways for unusual ligands. Recently, we reported the experimental validation of three distinct bacterial riboswitch classes that regulate gene expression in response to the selective binding of a guanidinium ion. These riboswitch classes, called guanidine-I, -II, and -III, regulate numerous genes whose protein products include previously misannotated guanidine exporters and enzymes that degrade guanidine via an initial carboxylation reaction. Guanidine is now recognized as the primal substrate of many multidrug efflux pumps that are important for bacterial resistance to certain antibiotics. Guanidine carboxylase enzymes had long been annotated as urea carboxylase enzymes but are now understood to participate in guanidine degradation. Herein, we report the existence of a fourth riboswitch class for this ligand, called guanidine-IV. Members of this class use a novel aptamer to selectively bind guanidine and use an unusual expression platform arrangement that is predicted to activate gene expression when ligand is present. The wide distribution of this abundant riboswitch class, coupled with the striking diversity of other guanidine-sensing RNAs, demonstrates that many bacterial species maintain sophisticated sensory and genetic mechanisms to avoid guanidine toxicity. This finding further highlights the mystery regarding the natural source of this nitrogen-rich chemical moiety.

Do the P1 and P2 hairpins of the Guanidine-II riboswitch interact?

Riboswitches regulate genes by adopting different structures in responds to metabolite binding. The guanidine-II riboswitch is the smallest representative of the ykkC class with the mechanism of its function being centred on the idea that its two stem loops P1 and P2 form a kissing hairpin interaction upon binding of guanidinium (Gdm⁺). This mechanism is based on in-line probing experiments with the full-length riboswitch and crystal structures of the truncated stem loops P1 and P2. However, the crystal structures reveal only the formation of the homodimers P1 | P1 and P2 | P2 but not of the proposed heterodimer P1 | P2. Here, site-directed spin labeling (SDSL) in combination with Pulsed Electron–Electron Double Resonance (PELDOR or DEER) is used to study their structures in solution and how they change upon binding of Gdm⁺. It is found that both hairpins adopt different structures in solution and that binding of Gdm⁺ does indeed lead to the formation of the heterodimer but alongside the homodimers in a statistical 1:2:1 fashion. These results do thus support the proposed switching mechanism.

Structural Insights into RNA Dimerization: Motifs, Interfaces and Functions

In comparison with the pervasive use of protein dimers and multimers in all domains of life, functional RNA oligomers have so far rarely been observed in nature. Their diminished occurrence contrasts starkly with the robust intrinsic potential of RNA to multimerize through long-range base-pairing ("kissing") interactions, self-annealing of palindromic or complementary sequences, and stable tertiary contact motifs, such as the GNRA tetraloop-receptors. To explore the general mechanics of RNA dimerization, we performed a meta-analysis of a collection of exemplary RNA homodimer structures consisting of viral genomic elements, ribozymes, riboswitches, etc., encompassing both functional and fortuitous dimers. Globally, we found that domain-swapped dimers and antiparallel, head-to-tail arrangements are predominant architectural themes. Locally, we observed that the same structural motifs, interfaces and forces that enable tertiary RNA folding also drive their higher-order assemblies. These feature prominently long-range kissing loops, pseudoknots, reciprocal base intercalations and A-minor interactions. We postulate that the scarcity of functional RNA multimers and limited diversity in multimerization motifs may reflect evolutionary constraints imposed by host antiviral immune surveillance and stress sensing. A deepening mechanistic understanding of RNA multimerization is expected to facilitate investigations into RNA and RNP assemblies, condensates, and granules and enable their potential therapeutical targeting.

Structure and ligand binding of the glutamine-II riboswitch

We have determined the structure of the glutamine-II riboswitch ligand binding domain using X-ray crystallography. The structure was solved using a novel combination of homology modeling and molecular replacement. The structure comprises three coaxial helical domains, the central one of which is a pseudoknot with partial triplex character. The major groove of this helix provides the binding site for L-glutamine, which is extensively hydrogen bonded to the RNA. Atomic mutation of the RNA at the ligand binding site leads to loss of binding shown by isothermal titration calorimetry, explaining the specificity of the riboswitch. A metal ion also plays an important role in ligand binding. This is directly bonded to a glutamine carboxylate oxygen atom, and its remaining inner-sphere water molecules make hydrogen bonding interactions with the RNA.