Riboswitch Aptamer Domain Research Articles

Opportunities for Riboswitch Inhibition by Targeting Co-Transcriptional RNA Folding Events

Antibiotic resistance is a critical global health concern, causing millions of prolonged bacterial infections every year and straining our healthcare systems. Novel antibiotic strategies are essential to combating this health crisis and bacterial non-coding RNAs are promising targets for new antibiotics. In particular, a class of bacterial non-coding RNAs called riboswitches has attracted significant interest as antibiotic targets. Riboswitches reside in the 5′-untranslated region of an mRNA transcript and tune gene expression levels in cis by binding to a small-molecule ligand. Riboswitches often control expression of essential genes for bacterial survival, making riboswitch inhibitors an exciting prospect for new antibacterials. Synthetic ligand mimics have predominated the search for new riboswitch inhibitors, which are designed based on static structures of a riboswitch’s ligand-sensing aptamer domain or identified by screening a small-molecule library. However, many small-molecule inhibitors that bind an isolated riboswitch aptamer domain with high affinity in vitro lack potency in vivo. Importantly, riboswitches fold and respond to the ligand during active transcription in vivo. This co-transcriptional folding is often not considered during inhibitor design, and may explain the discrepancy between a low Kd in vitro and poor inhibition in vivo. In this review, we cover advances in riboswitch co-transcriptional folding and illustrate how intermediate structures can be targeted by antisense oligonucleotides—an exciting new strategy for riboswitch inhibitor design.

Generative Modeling of RNA Sequence Families with Restricted Boltzmann Machines.

In this chapter, we discuss the potential application of Restricted Boltzmann machines (RBM) to model sequence families of structured RNA molecules. RBMs are a simple two-layer machine learning model able to capture intricate sequence dependencies induced by secondary and tertiary structure, as well as mechanisms of structural flexibility, resulting in a model that can be successfully used for the design of allosteric RNA such as riboswitches. They have recently been experimentally validated as generative models for the SAM-I riboswitch aptamer domain sequence family. We introduce RBM mathematically and practically, providing self-contained code examples to download the necessary training sequence data, train the RBM, and sample novel sequences. We present in detail the implementation of algorithms necessary to use RBMs, focusing on applications in biological sequence modeling.

Ionic Cooperativity between Lysine and Potassium in the Lysine Riboswitch: Single-Molecule Kinetic and Thermodynamic Studies.

Functionality in many biological systems, including proteins and nucleic acid structures, including protein and nucleic acid riboswitch structures, can depend on cooperative kinetic behavior between multiple small molecule ligands. In this work, single-molecule FRET data on the Bacillus subtilis lysine riboswitch reveals that affinity for the cognate lysine ligand increases significantly with K+, providing evidence for synergism between lysine/K+ binding to the aptamer and successful folding of the riboswitch. To describe/interpret this more complex kinetic scenario, we explore the conventional 4-state ("square") model for aptamer binding as a function of K+. Extension into this additional dimension generates a novel "cube" model for riboswitch folding dynamics with respect to lysine/K+ binding, revealing that riboswitch folding (kfold) and unfolding (kunfold) rate constants increase and decrease dramatically with K+, respectively. Furthermore, temperature-dependent single-molecule kinetic studies indicate that the presence of K+ entropically enhances the transition state barrier to folding but partially compensates for this by increasing the overall exothermicity for lysine binding. We rationalize this behavior as evidence that K+ facilitates hydrogen bonding between the negatively charged carboxyl group of lysine and the RNA, increasing structural rigidity and lowering entropy in the binding pocket. Finally, we explore the effects of cation size with Na+ and Cs+ studies to demonstrate that K+ is optimally suited for bridging interactions between lysine and the riboswitch aptamer domain. Regulation of lysine production and transport, dictated by the riboswitch's ability to recognize and bind lysine, is therefore intimately tied to the presence of K+ in the binding pocket and is strongly modulated by local cation conditions. The results suggest an increase in lysine riboswitch functionality by sensitivity to additional species in the cellular riboswitch environment.

Visualizing RNA conformational and architectural heterogeneity in solution

RNA flexibility is reflected in its heterogeneous conformation. Through direct visualization using atomic force microscopy (AFM) and the adenosylcobalamin riboswitch aptamer domain as an example, we show that a single RNA sequence folds into conformationally and architecturally heterogeneous structures under near-physiological solution conditions. Recapitulated 3D topological structures from AFM molecular surfaces reveal that all conformers share the same secondary structural elements. Only a population-weighted cohort, not any single conformer, including the crystal structure, can account for the ensemble behaviors observed by small-angle X-ray scattering (SAXS). All conformers except one are functionally active in terms of ligand binding. Our findings provide direct visual evidence that the sequence-structure relationship of RNA under physiologically relevant solution conditions is more complex than the one-to-one relationship for well-structured proteins. The direct visualization of conformational and architectural ensembles at the single-molecule level in solution may suggest new approaches to RNA structural analyses.

Mechanism of ligand binding signal transduction in a transcriptional riboswitch.

Riboswitches are noncoding RNA molecules that regulate gene expression in bacteria. The riboswitch's aptamer domain (AD) senses the cognate ligand and induces a structural change in the expression platform (EP) to signal on/off to the gene expression machinery. Despite extensive studies, the question of how the ligand binding signal is transmitted from the AD to EP resulting in the structural transitions in the EP remains unanswered. We studied the ligand binding signal transduction from AD to EP in the preQ1-sensing transcriptional riboswitch. In the absence of preQ1, both AD and EP individually form stem-loops with no structural overlap between the two stem-loops. The segment between the AD and EP stem loops is rich in adenine. In the presence of preQ1, AD forms a pseudoknot to encapsulate preQ1. We performed coarse-grained molecular dynamics simulations of full-length preQ1 riboswitch in the presence of Mg2+ and preQ1 ligand. In the absence of preQ1 and variation in Mg2+ concentration ([Mg2+]), we observed that AD populates only stem-loop and does not form the pseudoknot required to capture the ligand. Only in the presence of preQ1 AD attains the pseudoknot conformation by wrapping the A-rich region around the pre-existing stem-loop to facilitate ligand binding. The pseudoknot formation at AD allows EP to form the terminator helix and sequester the poly-U through strand invasion. The wrapping of the A-rich region around the AD stem-loop forces the stem-loop at EP to rupture and slide by 10 nucleotides towards 3'-end and form the terminator helix. The propensity of pseudoknot and terminator formation through strand invasion in the presence of preQ1 is enhanced with increment in [Mg2+]. We show that the competition between base-pair formation and helix sliding rate is crucial for the ligand binding signal transduction.

Disentangling contact and ensemble epistasis in a riboswitch

Mutations introduced into macromolecules often exhibit epistasis, where the effect of one mutation alters the effect of another. Knowing the mechanisms that lead to epistasis is important for understanding how macromolecules work and evolve, as well as for effective macromolecular engineering. Here, we investigate the interplay between “contact epistasis” (epistasis arising from physical interactions between mutated residues) and “ensemble epistasis” (epistasis that occurs when a mutation redistributes the conformational ensemble of a macromolecule, thus changing the effect of the second mutation). We argue that the two mechanisms can be distinguished in allosteric macromolecules by measuring epistasis at differing allosteric effector concentrations. Contact epistasis manifests as nonadditivity in the microscopic equilibrium constants describing the conformational ensemble. This epistatic effect is independent of allosteric effector concentration. Ensemble epistasis manifests as nonadditivity in thermodynamic observables—such as ligand binding—that are determined by the distribution of ensemble conformations. This epistatic effect strongly depends on allosteric effector concentration. Using this framework, we experimentally investigated the origins of epistasis in three pairwise mutant cycles introduced into the adenine riboswitch aptamer domain by measuring ligand binding as a function of allosteric effector concentration. We found evidence for both contact and ensemble epistasis in all cycles. Furthermore, we found that the two mechanisms of epistasis could interact with each other. For example, in one mutant cycle we observed 6 kcal/mol of contact epistasis in a microscopic equilibrium constant. In that same cycle, the maximum epistasis in ligand binding was only 1.5 kcal/mol: shifts in the ensemble masked the contribution of contact epistasis. Finally, our work yields simple heuristics for identifying contact and ensemble epistasis based on measurements of a biochemical observable as a function of allosteric effector concentration.

RNA Chemical Labeling with Site-Specific, Relative Quantification by Mass Spectrometry for the Structural Study of a Neomycin-Sensing Riboswitch Aptamer Domain.

High-resolution mass spectrometry was used for the label-free, direct localization and relative quantification of CMC+ -modifications of a neomycin-sensing riboswitch aptamer domain in the absence and presence of the aminoglycoside ligands neomycin B, ribostamycin, and paromomycin. The chemical probing and MS data for the free riboswitch show high exposure to solvent of the uridine nucleobases U7, U8, U13, U14, U18 as part of the proposed internal and apical loops, but those of U10 and U21 as part of the proposed internal loop were found to be far less exposed than expected. Thus, our data are in better agreement with the proposed secondary structure of the riboswitch in complexes with aminoglycosides than with that of free RNA. For the riboswitch in complexes with neomycin B, ribostamycin, and paromomycin, we found highly similar CMC+ -modification patterns and excellent agreement with previous NMR studies. Differences between the chemical probing and MS data in the absence and presence of the aminoglycoside ligands were quantitative rather than qualitative (i. e., the same nucleobases were labeled, but to different extents) and can be rationalized by stabilization of both the proposed bulge and the apical loop by aminoglycoside binding. Our study shows that chemical probing and mass spectrometry can provide important structural information and complement other techniques such as NMR spectroscopy.

An uncommon [K+(Mg2+)2] metal ion triad imparts stability and selectivity to the Guanidine-I riboswitch.

The widespread ykkC-I riboswitch class exemplifies divergent riboswitch evolution. To analyze how natural selection has diversified its versatile RNA fold, we determined the X-ray crystal structure of the Burkholderia sp. TJI49 ykkC-I subtype-1 (Guanidine-I) riboswitch aptamer domain. Differing from the previously reported structures of orthologs from Dickeya dadantii and Sulfobacillus acidophilus, our Burkholderia structure reveals a chelated K+ ion adjacent to two Mg2+ ions in the guanidine-binding pocket. Thermal melting analysis shows that K+ chelation, which induces localized conformational changes in the binding pocket, improves guanidinium-RNA interactions. Analysis of ribosome structures suggests that the [K+(Mg2+)2] ion triad is uncommon. It is, however, reminiscent of metal ion clusters found in the active sites of ribozymes and DNA polymerases. Previous structural characterization of ykkC-I subtype-2 RNAs, which bind the effector ligands ppGpp and PRPP, indicate that in those paralogs, an adenine responsible for K+ chelation in the Burkholderia Guanidine-I riboswitch is replaced by a pyrimidine. This mutation results in a water molecule and Mg2+ ion binding in place of the K+ ion. Thus, our structural analysis demonstrates how ion and solvent chelation tune divergent ligand specificity and affinity among ykkC-I riboswitches.

Riboswitches for Controlled Expression of Therapeutic Transgenes Delivered by Adeno-Associated Viral Vectors.

Vectors developed from adeno-associated virus (AAV) are powerful tools for in vivo transgene delivery in both humans and animal models, and several AAV-delivered gene therapies are currently approved for clinical use. However, AAV-mediated gene therapy still faces several challenges, including limited vector packaging capacity and the need for a safe, effective method for controlling transgene expression during and after delivery. Riboswitches, RNA elements which control gene expression in response to ligand binding, are attractive candidates for regulating expression of AAV-delivered transgene therapeutics because of their small genomic footprints and non-immunogenicity compared to protein-based expression control systems. In addition, the ligand-sensing aptamer domains of many riboswitches can be exchanged in a modular fashion to allow regulation by a variety of small molecules, proteins, and oligonucleotides. Riboswitches have been used to regulate AAV-delivered transgene therapeutics in animal models, and recently developed screening and selection methods allow rapid isolation of riboswitches with novel ligands and improved performance in mammalian cells. This review discusses the advantages of riboswitches in the context of AAV-delivered gene therapy, the subsets of riboswitch mechanisms which have been shown to function in human cells and animal models, recent progress in riboswitch isolation and optimization, and several examples of AAV-delivered therapeutic systems which might be improved by riboswitch regulation.

RNA Ensembles from Solvent Accessibility Data: Application to the SAM-I Riboswitch Aptamer Domain

Riboswitches are regulatory ribonucleic acid (RNA) elements that act as ligand-dependent conformational switches that recognize their cognate ligand via a binding pocket located in their aptamer domain. In the apo form, the aptamer domain is dynamic, requiring an ensemble representation of its structure. Here, as a proof-of-concept, we used solvent accessibility information to construct a pair of dynamical ensembles of the aptamer domain of the well-studied S-adenosylmethionine (SAM) class-I riboswitch in the absence (-SAM) and presence (+SAM) of SAM. To achieve this, we first generated a large conformational library and then reweighted conformers in the library using solvent-accessible surface area (SASA) data derived from recently reported light-activated structural examination of RNA (LASER) reactivities measured in the -SAM and +SAM states of the riboswitch. The differences in the resulting -SAM and +SAM ensembles are consistent with a SAM-dependent reshaping of the free-energy landscape of the aptamer domain. Within our -SAM ensemble, we identified a "transient" state that is missing a critical long-range contact, leading us to speculate that it may be representative of a folding intermediate. Further structural analysis also revealed that the transient state harbors a hidden binding pocket that is distinct from the SAM-binding pocket and is predicted by docking calculations to selectively bind small-molecule ligands. The SASA-based method we applied to the SAM-I riboswitch aptamer domain is general and could be used to construct dynamical ensembles for other riboswitch aptamer domains and, more broadly, other classes of structured RNAs.

Tertiary RNA Folding-Targeted Drug Screening Strategy Using a Protein Nanopore.

Bacterial riboswitch RNAs are attractive targets for novel antibiotics against antibiotic-resistant superbacteria. Their binding to cognate metabolites is essential for the regulation of bacterial gene expression. Despite the importance of RNAs as therapeutic targets, the development of RNA-targeted, small-molecule drugs is limited by current biophysical methods. Here, we monitored the specific interaction between the adenine-sensing riboswitch aptamer domain (ARS) and adenine at the single-molecule level using α-hemolysin (αHL) nanopores. During adenine-induced tertiary folding, adenine-bound ARS intermediates exhibited characteristic nanopore events, including a two-level ionic current blockade and a ∼ 5.6-fold longer dwell time than that of free RNA. In a proof-of-concept experiment, tertiary RNA folding-targeted drug screening was performed using a protein nanopore, which resulted in the discovery of three new ARS-targeting hit compounds from a natural compound library. Taken together, these results reveal that αHL nanopores are a valuable platform for ultrasensitive, label-free, and single-molecule-based drug screening against therapeutic RNA targets.

Aromatic 19 F-13 C TROSY-[19 F, 13 C]-Pyrimidine Labeling for NMR Spectroscopy of RNA.

We present the access to [5‐19F, 5‐13C]‐uridine and ‐cytidine phosphoramidites for the production of site‐specifically modified RNAs up to 65 nucleotides (nts). The amidites were used to introduce [5‐19F, 5‐13C]‐pyrimidine labels into five RNAs—the 30 nt human immunodeficiency virus trans activation response (HIV TAR) 2 RNA, the 61 nt human hepatitis B virus ϵ (hHBV ϵ) RNA, the 49 nt SAM VI riboswitch aptamer domain from B. angulatum, the 29 nt apical stem loop of the pre‐microRNA (miRNA) 21 and the 59 nt full length pre‐miRNA 21. The main stimulus to introduce the aromatic 19F–13C‐spin topology into RNA comes from a work of Boeszoermenyi et al., in which the dipole‐dipole interaction and the chemical shift anisotropy relaxation mechanisms cancel each other leading to advantageous TROSY properties shown for aromatic protein sidechains. This aromatic 13C–19F labeling scheme is now transferred to RNA. We provide a protocol for the resonance assignment by solid phase synthesis based on diluted [5‐19F, 5‐13C]/[5‐19F] pyrimidine labeling. For the 61 nt hHBV ϵ we find a beneficial 19F–13C TROSY enhancement, which should be even more pronounced in larger RNAs and will facilitate the NMR studies of larger RNAs. The [19F, 13C]‐labeling of the SAM VI aptamer domain and the pre‐miRNA 21 further opens the possibility to use the biorthogonal stable isotope reporter nuclei in in vivo NMR to observe ligand binding and microRNA processing in a biological relevant setting.

A 300-fold enhancement of imino nucleic acid resonances by hyperpolarized water provides a new window for probing RNA refolding by 1D and 2D NMR

NMR sensitivity-enhancement methods involving hyperpolarized water could be of importance for solution-state biophysical investigations. Hyperpolarized water (HyperW) can enhance the 1H NMR signals of exchangeable sites by orders of magnitude over their thermal counterparts, while providing insight into chemical exchange and solvent accessibility at a site-resolved level. As HyperW's enhancements are achieved by exploiting fast solvent exchanges associated with minimal interscan delays, possibilities for the rapid monitoring of chemical reactions and biomolecular (re)folding are opened. HyperW NMR can also accommodate heteronuclear transfers, facilitating the rapid acquisition of 2-dimensional (2D) 15N-1H NMR correlations, and thereby combining an enhanced spectral resolution with speed and sensitivity. This work demonstrates how these qualities can come together for the study of nucleic acids. HyperW injections were used to target the guanine-sensing riboswitch aptamer domain (GSRapt) of the xpt-pbuX operon in Bacillus subtilis Unlike what had been observed in proteins, where residues benefited of HyperW NMR only if/when sufficiently exposed to water, these enhancements applied to every imino resonance throughout the RNA. The >300-fold enhancements observed in the resulting 1H NMR spectra allowed us to monitor in real time the changes that GSRapt undergoes upon binding hypoxanthine, a high-affinity interaction leading to conformational refolding on a ∼1-s timescale at 36 °C. Structural responses could be identified for several nucleotides by 1-dimensional (1D) imino 1H NMR as well as by 2D HyperW NMR spectra acquired upon simultaneous injection of hyperpolarized water and hypoxanthine. The folding landscape revealed by this HyperW strategy for GSRapt, is briefly discussed.

Conformational flexibility of adenine riboswitch aptamer in apo and bound states using NMR and an X-ray free electron laser

Riboswitches are structured cis-regulators mainly found in the untranslated regions of messenger RNA. The aptamer domain of a riboswitch serves as a sensor for its ligand, the binding of which triggers conformational changes that regulate the behavior of its expression platform. As a model system for understanding riboswitch structures and functions, the add adenine riboswitch has been studied extensively. However, there is a need for further investigation of the conformational dynamics of the aptamer in light of the recent real-time crystallographic study at room temperature (RT) using an X-ray free electron laser (XFEL) and femtosecond X-ray crystallography (SFX). Herein, we investigate the conformational motions of the add adenine riboswitch aptamer domain, in the presence or absence of adenine, using nuclear magnetic resonance relaxation measurements and analysis of RT atomic displacement factors (B-factors). In the absence of ligand, the P1 duplex undergoes a fast exchange where the overall molecule exhibits a motion at kex ~ 319 s−1, based on imino signals. In the presence of ligand, the P1 duplex adopts a highly ordered conformation, with kex~ 83 s−1, similar to the global motion of the molecule, excluding the loops and binding pocket, at 84 s−1. The µs–ms motions in both the apo and bound states are consistent with RT B-factors. Reduced spatial atomic fluctuation, ~ 50%, in P1 upon ligand binding coincides with significantly attenuated temporal dynamic exchanges. The binding pocket is structured in the absence or presence of ligand, as evidenced by relatively low and similar RT B-factors. Therefore, despite the dramatic rearrangement of the binding pocket, those residues exhibit similar spatial thermal fluctuation before and after binding.

Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: crystal structure analysis, classification, optimal structures, and accurate interaction energies.

In the present work, 67 crystal structures of the aptamer domains of RNA riboswitches are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of 80 unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson–Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate, and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that 67 out of 80 pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal a rich hydrogen bonding landscape, ranging from weak interactions (base:other, −3 kcal/mol) to strong (base:phosphate, −48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in the context of RNA:ligand recognition. This can, in turn, contribute toward current efforts to develop antimicrobials that target RNAs.

A Mn-sensing riboswitch activates expression of a Mn2+/Ca2+ ATPase transporter in Streptococcus.

Maintaining manganese (Mn) homeostasis is important for the virulence of numerous bacteria. In the human respiratory pathogen Streptococcus pneumoniae, the Mn-specific importer PsaBCA, exporter MntE, and transcriptional regulator PsaR establish Mn homeostasis. In other bacteria, Mn homeostasis is controlled by yybP-ykoY family riboswitches. Here, we characterize a yybP-ykoY family riboswitch upstream of the mgtA gene encoding a PII-type ATPase in S. pneumoniae, suggested previously to function in Ca2+ efflux. We show that the mgtA riboswitch aptamer domain adopts a canonical yybP-ykoY structure containing a three-way junction that is compacted in the presence of Ca2+ or Mn2+ at a physiological Mg2+ concentration. Although Ca2+ binds to the RNA aptamer with higher affinity than Mn2+, in vitro activation of transcription read-through of mgtA by Mn2+ is much greater than by Ca2+. Consistent with this result, mgtA mRNA and protein levels increase ≈5-fold during cellular Mn stress, but only in genetic backgrounds of S. pneumoniae and Bacillus subtilis that exhibit Mn2+ sensitivity, revealing that this riboswitch functions as a failsafe ‘on’ signal to prevent Mn2+ toxicity in the presence of high cellular Mn2+. In addition, our results suggest that the S. pneumoniae yybP-ykoY riboswitch functions to regulate Ca2+ efflux under these conditions.

Dynamics of metal ions around an RNA molecule.

Noncoding RNA molecules take part in many biological processes, while metal ions play crucial roles in helping RNAs to perform their functions. However, the statics and dynamics of these metal ions around RNA molecules are still not well understood. In this work, we report a detailed molecular dynamics study of the type-I preQ_{1}-bound riboswitch aptamer domain (PRAD) at different ionic conditions (K^{+}, Na^{+}, and Mg^{2+}). The results show that the structural properties and flexibility of the PRAD molecule greatly influence the distributions and dynamics of metal ions around it. Simultaneously, Na^{+} ions show a stronger competitiveness with Mg^{2+} ions than K^{+} ions, and the three types of metal ions have different modes of interaction with the RNA molecule. Furthermore, we have also investigated specific binding sites of metal ions on the PRAD molecule and found that the dynamics and hydration structures of metal ions located at the ion-binding sites were obviously affected by the RNA structure near these ion-binding sites. These results may be useful to understand the role of the metal ions in noncoding RNA functions.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches.

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

Crowding Shifts the FMN Recognition Mechanism of Riboswitch Aptamer from Conformational Selection to Induced Fit.

In bacteria, the binding between the riboswitch aptamer domain and ligand is regulated by environmental cues, such as low Mg2+ in macrophages during pathogenesis to ensure spatiotemporal expression of virulence genes. Binding was investigated between the flavin mononucleotide (FMN) riboswitch aptamer and its anionic ligand in the presence of molecular crowding agent without Mg2+ ion, which mimics pathogenic conditions. Structural, kinetic, and thermodynamic analyses under the crowding revealed more dynamic conformational rearrangements of the FMN riboswitch aptamer compared to dilute Mg2+ -containing solution. It is hypothesized that under crowding conditions FMN binds through an induced fit mechanism in contrast to the conformational selection mechanism previously demonstrated in dilute Mg2+ solution. Since these two mechanisms involve different conformational intermediates and rate constants, these findings have practical significance in areas such as drug design and RNA engineering.

Crowding Shifts the FMN Recognition Mechanism of Riboswitch Aptamer from Conformational Selection to Induced Fit

AbstractIn bacteria, the binding between the riboswitch aptamer domain and ligand is regulated by environmental cues, such as low Mg2+ in macrophages during pathogenesis to ensure spatiotemporal expression of virulence genes. Binding was investigated between the flavin mononucleotide (FMN) riboswitch aptamer and its anionic ligand in the presence of molecular crowding agent without Mg2+ ion, which mimics pathogenic conditions. Structural, kinetic, and thermodynamic analyses under the crowding revealed more dynamic conformational rearrangements of the FMN riboswitch aptamer compared to dilute Mg2+‐containing solution. It is hypothesized that under crowding conditions FMN binds through an induced fit mechanism in contrast to the conformational selection mechanism previously demonstrated in dilute Mg2+solution. Since these two mechanisms involve different conformational intermediates and rate constants, these findings have practical significance in areas such as drug design and RNA engineering.