A magnesium-induced triplex pre-organizes the SAM-II riboswitch

Understanding the Thermodynamics of Magnesium Binding to RNA Structural Motifs

Divalent magnesium ions (Mg2+) serve a vital role in defining the structural and catalytic chemistry of a wide array of RNA molecules. The body of structural information on RNA motifs continues to expand and, in turn, the functional importance of Mg2+ is revealed. A combination of prior work on the structural characterization of magnesium binding ligands with inner- and outer-sphere coordination modes, with recorded experimental binding energies for inner- and outer-sphere contacts, demonstrates the relative affinity and thermodynamic hierarchy for these sites. In turn, these can be correlated with cellular concentrations of free available magnesium ions, allowing the prioritization of populating important functional sites and a correlation with physiological function. This paper summarizes some of the key results of that analysis and provides predictive rules for the affinity and role of newly identified Mg binding sites on complex RNA structures. The influence of crystal packing on magnesium binding to RNA motifs, relative to their solution form, is addressed and caveats made.

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods

RNA, a highly charged biopolymer composed of negatively charged phosphate groups, defies electrostatic repulsion to adopt well-defined, compact structures. Hence, the presence of positively charged metal ions is crucial not only for RNA's charge neutralization, but they also coherently decorate the ion atmosphere of RNA to stabilize its compact fold. This feature article elucidates various modes of close RNA-ion interactions, with a special emphasis on Mg2+ as an outer-sphere and inner-sphere ion. Through examples, we highlight how inner-sphere chelated Mg2+ stabilizes RNA pseudoknots, while outer-sphere ions can also exert long-range electrostatic interactions, inducing groove narrowing, coaxial helical stacking, and RNA ring formation. In addition to investigating the RNA's ion environment, we note that the RNA's hydration environment is relatively underexplored. Our study delves into its profound interplay with the structural dynamics of RNA, employing state-of-the-art atomistic simulation techniques. Through examples, we illustrate how specific ions and water molecules are associated with RNA functions, leveraging atomistic simulations to identify preferential ion binding and hydration sites. However, understanding their impact(s) on the RNA structure remains challenging due to the involvement of large length and long time scales associated with RNA's dynamic nature. Nevertheless, our contributions and recent advances in coarse-grained simulation techniques offer insights into large-scale structural changes dynamically linked to the RNA ion atmosphere. In this connection, we also review how different cutting-edge computational simulation methods provide a microscopic lens into the influence of ions and hydration on RNA structure and dynamics, elucidating distinct ion atmospheric components and specific hydration layers and their individual and collective impacts.

Dynamic Counterion Condensation Model Decodes Functional Dynamics of RNA Pseudoknot in SARS-CoV-2: Control of Ion-Mediated Pierced Lasso Topology

The programmed frameshifting stimulatory element, a promising drug target for COVID-19 treatment, involves a RNA pseudoknot (PK) structure. This RNA PK facilitates frameshifting, enabling RNA viruses to translate multiple proteins from a single mRNA, which is a key strategy for their rapid evolution. Overcoming the challenges of capturing large-scale structural changes of RNA under the influence of a dynamic counterion environment (K+ and Mg2+), the study extended the applications of a newly developed dynamic counterion condensation (DCC) model. DCC simulations reveal potential folding pathways of this RNA PK, supported by the experimental findings obtained using optical tweezers. The study elucidates the pivotal role of Mg2+ ions in crafting a lasso-like RNA topology, a novel RNA motif that governs dynamic transitions between the ring-opened and ring-closed states of the RNA. The pierced lasso component guided by Mg2+-mediated interactions orchestrates inward and outward motion fine-tuning tension on the slippery segment, a critical factor for optimizing frameshifting efficiency.

Structural and dynamic mechanisms for coupled folding and tRNA recognition of a translational T-box riboswitch

T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg2+-induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg2+, the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.

Magnesium Ion-Driven Folding and Conformational Switching Kinetics of Tetracycline Binding Aptamer: Implications for in vivo Riboswitch Engineering

Engineering in vitro selected RNA aptamers into in vivo functional riboswitches represents a long-standing challenge in molecular biology. The highly specific aptamer domain of the riboswitch undergoes a conformational adjustment in response to ligand sensing, which in turn exerts the regulatory function. Besides essential factors like structural complexity and ligand binding kinetics, the active role of magnesium ions in stabilizing RNA tertiary structures and assisting in ligand binding can be a vital criterion. We present spectroscopic studies on the magnesium ion-driven folding of the Tetracycline binding aptamer. Using fluorescent labels, the aptamer pre-folding and subsequent ligand binding is monitored by magnesium titration experiments and time-resolved stopped-flow measurements. A minimum concentration of 0.5 mM magnesium is required to fold into a magnesium ion-stabilized binding-competent state with a preformed binding pocket. Tetracycline binding causes a pronounced conformational change that results in the establishment of the triple helix core motif, and that further propagates towards the closing stem. By a dynamic acquisition of magnesium ions, a kink motif is formed at the intersection of the triple helix and closing stem regions. This ultimately entails a stabilization of the closing stem which is discussed as a key element in the regulatory function of the Tetracycline aptamer.

Hierarchical Hydration Dynamics of RNA with Nano-Water-Pool at Its Core

Many functional RNAs fold into a compact, roughly globular shape by minimizing the electrostatic repulsion between their negatively charged phosphodiester backbone. The fold of such close, compact RNA architecture is often so designed that its outer surface and complex core both are predominately populated by phosphate groups loosely sequestering bases in the intermediate layers. A number of helical junctions maintain the RNA core and its nano-water-pool. While the folding of RNA is manifested by its counterion environment composed of mixed mono- and divalent salts, the concerted role of ion and water in maintaining an RNA fold is yet to be explored. In this work, detailed atomistic simulations of SAM-I and Add Adenine riboswitch aptamers, and subgenomic flavivirus RNA (sfRNA) have been performed in a physiological mixed mono- and divalent salt environment. All three RNA systems have compact folds with a core diameter of range 1-1.7 nm. The spatiotemporal heterogeneity of RNA hydration was probed in a layer-wise manner by distinguishing the core, the intermediate, and the outer layers. The layer-wise decomposition of hydrogen bonds and collective single-particle reorientational dynamics reveal a nonmonotonic relaxation pattern with the slowest relaxation observed at the intermediate layers that involves functionally important tertiary motifs. The slowness of this intermediate layer is attributed to two types of long-resident water molecules: (i) water from ion-hydration layers and (ii) structurally trapped water (distant from ions). The relaxation kinetics of the core and the surface water essentially exposed to the phosphate groups show well-separated time scales from the intermediate layers. In the slow intermediate layers, site-specific ions and water control the functional dynamics of important RNA motifs like kink-turn, observed in different structure-probing experiments. Most interestingly, we find that as the size of the RNA core increases (SAM1 core < sfRNAcore < Add adenine core), its hydration tends to show faster relaxation. The hierarchical hydration and the layer-wise base-phosphate composition uniquely portray the globular RNA to act like a soft vesicle with a quasi-dynamic nano-water-pool at its core.

Structural insights into translation regulation by the THF-II riboswitch

In bacteria, expression of folate-related genes is controlled by the tetrahydrofolate (THF) riboswitch in response to specific binding of THF and its derivatives. Recently, a second class of THF riboswitches, named THF-II, was identified in Gram-negative bacteria, which exhibit distinct architecture from the previously characterized THF-I riboswitches found in Gram-positive bacteria. Here, we present the crystal structures of the ligand-bound THF-II riboswitch from Mesorhizobium loti. These structures exhibit a long rod-like fold stabilized by continuous base pair and base triplet stacking across two helices of P1 and P2 and their interconnecting ligand-bound binding pocket. The pterin moiety of the ligand docks into the binding pocket by forming hydrogen bonds with two highly conserved pyrimidines in J12 and J21, which resembles the hydrogen-bonding pattern at the ligand-binding site FAPK in the THF-I riboswitch. Using small-angle X-ray scattering and isothermal titration calorimetry, we further characterized the riboswitch in solution and reveal that Mg2+ is essential for pre-organization of the binding pocket for efficient ligand binding. RNase H cleavage assay indicates that ligand binding reduces accessibility of the ribosome binding site in the right arm of P1, thus down-regulating the expression of downstream genes. Together, these results provide mechanistic insights into translation regulation by the THF-II riboswitch.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.

Alternative RNA Conformations: Companion or Combatant

RNA molecules, in one form or another, are involved in almost all aspects of cell physiology, as well as in disease development. The diversity of the functional roles of RNA comes from its intrinsic ability to adopt complex secondary and tertiary structures, rivaling the diversity of proteins. The RNA molecules form dynamic ensembles of many interconverting conformations at a timescale of seconds, which is a key for understanding how they execute their cellular functions. Given the crucial role of RNAs in various cellular processes, we need to understand the RNA molecules from a structural perspective. Central to this review are studies aimed at revealing the regulatory role of conformational equilibria in RNA in humans to understand genetic diseases such as cancer and neurodegenerative diseases, as well as in pathogens such as bacteria and viruses so as to understand the progression of infectious diseases. Furthermore, we also summarize the prior studies on the use of RNA structures as platforms for the rational design of small molecules for therapeutic applications.

Dynamic effects of the spine of hydrated magnesium on viral RNA pseudoknot structure

In the cellular environment, a viral RNA Pseudoknot (PK) structure is responsive to its prevailing ion atmosphere created by a mixture of monovalent and divalent cations. We investigate the influence of such a mixed-salt environment on RNA-PK structure at an atomic resolution through three sets of 1.5 μs explicit solvent molecular dynamics (MD) simulations and also by building a dynamic counterion-condensation (DCC) model at varying divalent Mg²⁺ concentrations. The DCC model includes explicit interaction of the ligand and adjacent chelated Mg²⁺ by extending the recently developed generalized Manning condensation model. Its implementation within an all-atom structure-based molecular dynamics framework bolsters its opportunity to explore large-length scale and long-timescale phenomena associated with complex RNA systems immersed in its dynamic ion environment. In the present case of RNA-PK, both explicit MD and DCC simulations reveal a spine of hydrated-Mg²⁺ to induce stem-I and stem-II closure where the minor groove between these stems is akin to breathing. Mg²⁺ mediated minor-groove narrowing is coupled with local base-flipping dynamics of a base triple and quadruple, changing the stem structure of such RNA. Contrary to the conversational view of the indispensable need for Mg²⁺ for the tertiary structure of RNA, the study warns about the plausible detrimental effect of specific Mg²⁺-phosphate interactions on the RNA-PK structure beyond a certain concentration of Mg²⁺.

Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (Q_TC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Incoherent dual regulation by a SAM-II riboswitch controlling translation at a distance

In Sinorhizobium meliloti, the methionine biosynthesis genes metA and metZ are preceded by S-adenosyl-L-methionine (SAM) riboswitches of the SAM-II class. Upon SAM binding, structural changes in the metZ riboswitch were predicted to cause transcriptional termination, generating the sRNA RZ. By contrast, the metA riboswitch was predicted to regulate translation from an AUG1 codon. However, downstream of the metA riboswitch, we found a putative Rho-independent terminator and an in-frame AUG2 codon, which may contribute to metA regulation. We validated the terminator between AUG1 and AUG2, which generates the sRNA RA1 that is processed to RA2. Under high SAM conditions, the activities of the metA and metZ promoters and the steady-state levels of the read-through metA and metZ mRNAs were decreased, while the levels of the RZ and RA2 sRNAs were increased. Under these conditions, the sRNAs and the mRNAs were stabilized. Reporter fusion experiments revealed that the Shine–Dalgarno (SD) sequence in the metA riboswitch is required for translation, which, however, starts 74 nucleotides downstream at AUG2, suggesting a novel translation initiation mechanism. Further, the reporter fusion data supported the following model of RNA-based regulation: Upon SAM binding by the riboswitch, the SD sequence is sequestered to downregulate metA translation, while the mRNA is stabilized. Thus, the SAM-II riboswitches fulfil incoherent, dual regulation, which probably serves to ensure basal metA and metZ mRNA levels under high SAM conditions. This probably helps to adapt to changing conditions and maintain SAM homoeostasis.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

Chelated Magnesium Logic Gate Regulates Riboswitch Pseudoknot Formation

Magnesium plays a critical role in the structure, dynamics, and function of RNA. The precise microscopic effect of chelated magnesium on RNA structure is yet to be explored. Magnesium is known to act through its diffuse cloud around RNA, through the outer sphere (water-mediated), inner sphere, and often chelated ion-mediated interactions. A mechanism is proposed for the role of experimentally discovered site-specific chelated magnesium ions on the conformational dynamics of SAM-I riboswitch aptamers in bacteria. This mechanism is observed with atomistic simulations performed in a physiological mixed salt environment at a high temperature. The simulations were validated with phosphorothioate interference mapping experiments that help to identify crucial inner-sphere Mg²⁺ sites prescribing an appropriate initial distribution of inner- and outer-sphere magnesium ions to maintain a physiological ion concentration of monovalent and divalent salts. A concerted role of two chelated magnesium ions is newly discovered since the presence of both supports the formation of the pseudoknot. This constitutes a logical AND gate. The absence of any of these magnesium ions instigates the dissociation of long-range pseudoknot interaction exposing the inner core of the RNA. A base triple is the epicenter of the magnesium chelation effect. It allosterically controls RNA pseudoknot by bolstering the direct effect of magnesium chelation in protecting the functional fold of RNA to control ON and OFF transcription switching.

Mutations of N1 Riboswitch Affect its Dynamics and Recognition by Neomycin Through Conformational Selection

Short, structured fragments of non-coding mRNA may act as molecular switches upon binding specific ligands, regulating the translation of proteins encoded downstream this mRNA sequence. One switch, called riboswitch N1, is regulated by aminoglycosides such as neomycin. Nucleobase mutations in the apical loop, although distant from the binding pocket, significantly affect neomycin affinity and riboswitch regulatory efficiency. To explain this influence, we conducted molecular dynamics simulations using generalized replica exchange with solute tempering (gREST). Translation assay of a reporter protein in a yeast system shows that mutating A17 to G in the riboswitch apical loop reduces 6-fold the translation regulation efficiency of the mutant. Indeed, simulations of the unbound riboswitch show that G17 frequently stacks with base 7, while base 8 is stabilized towards the binding site in a way that it may interfere with the conformational selection mechanism and decrease riboswitch regulatory activity. In the riboswitch complexes, this single-point A to G mutation disrupts a strong hydrogen bond between nucleotides 5 and 17 and, instead, a new hydrogen bond between residue 17 and neomycin is created. This change forces neomycin to occupy a slightly shifted position in the binding pocket, which increases neomycin flexibility. Our simulations of the U14C mutation suggest that the riboswitch complex with neomycin is more stable if cytosine 14 is protonated. A hydrogen bond between the RNA phosphate and protonated cytosine appears as the stabilizing factor. Also, based on the cell-free translation assay and isothermal titration calorimetry experiments, mutations of nucleotides 14 and 15 affect only slightly the riboswitch ability to bind the ligand and its activity. Indeed, the simulation of the unbound U15A mutant suggests conformations preformed for ligand binding, which may explain slightly higher regulatory activity of this mutant. Overall, our results corroborate the in vivo and in vitro experiments on the N1 riboswitch-neomycin system, detail the relationship between nucleobase mutations and RNA dynamics, and reveal the conformations playing the major role in the conformational selection mechanism.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

Unraveling the structure and biological functions of RNA triple helices

It has been nearly 63 years since the first characterization of an RNA triple helix in vitro by Gary Felsenfeld, David Davies, and Alexander Rich. An RNA triple helix consists of three strands: A Watson–Crick RNA double helix whose major‐groove establishes hydrogen bonds with the so‐called “third strand”. In the past 15 years, it has been recognized that these major‐groove RNA triple helices, like single‐stranded and double‐stranded RNA, also mediate prominent biological roles inside cells. Thus far, these triple helices are known to mediate catalysis during telomere synthesis and RNA splicing, bind to ligands and ions so that metabolite‐sensing riboswitches can regulate gene expression, and provide a clever strategy to protect the 3′ end of RNA from degradation. Because RNA triple helices play important roles in biology, there is a renewed interest in better understanding the fundamental properties of RNA triple helices and developing methods for their high‐throughput discovery. This review provides an overview of the fundamental biochemical and structural properties of major‐groove RNA triple helices, summarizes the structure and function of naturally occurring RNA triple helices, and describes prospective strategies to isolate RNA triple helices as a means to establish the “triplexome”.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure and Dynamics

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

Magnesium controls aptamer-expression platform switching in the SAM-I riboswitch

Investigations of most riboswitches remain confined to the ligand-binding aptamer domain. However, during the riboswitch mediated transcription regulation process, the aptamer domain and the expression platform compete for a shared strand. If the expression platform dominates, an anti-terminator helix is formed, and the transcription process is active (ON state). When the aptamer dominates, transcription is terminated (OFF state). Here, we use an expression platform switching experimental assay and structure-based electrostatic simulations to investigate this ON-OFF transition of the full length SAM-I riboswitch and its magnesium concentration dependence. Interestingly, we find the ratio of the OFF population to the ON population to vary non-monotonically as magnesium concentration increases. Upon addition of magnesium, the aptamer domain pre-organizes, populating the OFF state, but only up to an intermediate magnesium concentration level. Higher magnesium concentration preferentially stabilizes the anti-terminator helix, populating the ON state, relatively destabilizing the OFF state. Magnesium mediated aptamer-expression platform domain closure explains this relative destabilization of the OFF state at higher magnesium concentration. Our study reveals the functional potential of magnesium in controlling transcription of its downstream genes and underscores the importance of a narrow concentration regime near the physiological magnesium concentration ranges, striking a balance between the OFF and ON states in bacterial gene regulation.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg²⁺ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

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.

Structure and ligand binding of the SAM-V riboswitch

SAM-V is one of the class of riboswitches that bind S-adenosylmethione, regulating gene expression by controlling translation. We have solved the crystal structure of the metY SAM-V riboswitch bound to its SAM ligand at 2.5 Å resolution. The RNA folds as an H-type pseudoknot, with a major-groove triple helix in which resides the SAM ligand binding site. The bound SAM adopts an elongated conformation aligned with the axis of the triple helix, and is held at either end by hydrogen bonding to the adenine and the amino acid moieties. The central sulfonium cation makes electrostatic interactions with an U:A.U base triple, so conferring specificity. We propose a model in which SAM binding leads to association of the triplex third strand that stabilizes a short helix and occludes the ribosome binding site. Thus the new structure explains both ligand specificity and the mechanism of genetic control.

Cooperation between Magnesium and Metabolite Controls Collapse of the SAM-I Riboswitch

The S-adenosylmethionine (SAM)-I riboswitch is a noncoding RNA that regulates the transcription termination process in response to metabolite (SAM) binding. The aptamer portion of the riboswitch may adopt an open or closed state depending on the presence of metabolite. Although the transition between the open and closed states is critical for the switching process, its atomistic details are not well understood. Using atomistic simulations, we calculate the effect of SAM and magnesium ions on the folding free energy landscape of the SAM-I riboswitch. These molecular simulation results are consistent with our previous wetlab experiments and aid in interpreting the SHAPE probing measurements. Here, molecular dynamics simulations explicitly identify target RNA motifs sensitive to magnesium ions and SAM. In the simulations, we observe that, whereas the metabolite mostly stabilizes the P1 and P3 helices, magnesium serves an important role in stabilizing a pseudoknot interaction between the P2 and P4 helices, even at high metabolite concentrations. The pseudoknot stabilization by magnesium, in combination with P1 stabilization by SAM, explains the requirement of both SAM and magnesium to form the fully collapsed metabolite-bound closed state of the SAM-I riboswitch. In the absence of SAM, frequent open-to-closed conformational transitions of the pseudoknot occur, akin to breathing. These pseudoknot fluctuations disrupt the binding site by facilitating fluctuations in the 5'-end of helix P1. Magnesium biases the landscape toward a collapsed state (preorganization) by coordinating pseudoknot and 5'-P1 fluctuations. The cooperation between SAM and magnesium in stabilizing important tertiary interactions elucidates their functional significance in transcription regulation.