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

Susmita Roy; Heiko Lammert; Ryan L. Hayes; Bin Chen; Regan LeBlanc; T. Kwaku Dayie; José N. Onuchic; Karissa Y. Sanbonmatsu

doi:10.1371/journal.pcbi.1005406

Abstract

Our ¹³C- and ¹H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.

Author summary

The presence of positively charged metal ions is essential to maintain the structural fold and function of RNA. Among different metal ions, magnesium is particularly important for the stability of RNA because it can efficiently support a close assembly of negatively charged phosphate groups in an RNA fold. The SAM-II riboswitch is an example of a classical pseudoknot fold, which binds S-adenosyl methionine, stabilizing an alternate folded form to inhibit gene expression. In our early ¹³C- and ¹H-chemical exchange saturation transfer (CEST) experiments, we found a conformational transition between a minor, partially closed and a major, open state conformation in the absence of ligand. Our CEST experiments at different magnesium concentrations now suggest that magnesium ions can induce a conformational pre-organization in the apo SAM-II riboswitch, which is expected to facilitate ligand binding. To understand the microscopic details of this magnesium-induced transition, we perform all-atom structure-based molecular simulations including electrostatics and explicit magnesium ions. Our free energy calculations reveal that the partially closed pre-organized state is further stabilized with increasing magnesium concentration. This is in excellent agreement with our ¹³C-CEST profile, SAXS, and size-exclusion chromatographic data, and with recent single molecule FRET experiments. Our results suggest that a sufficiently high concentration of magnesium is essential to pre-organize the apo SAM-II riboswitch.

Citation: Roy S, Lammert H, Hayes RL, Chen B, LeBlanc R, Dayie TK, et al. (2017) A magnesium-induced triplex pre-organizes the SAM-II riboswitch. PLoS Comput Biol 13(3): e1005406. https://doi.org/10.1371/journal.pcbi.1005406

Editor: Peter M. Kasson, University of Virginia, UNITED STATES

Received: November 11, 2016; Accepted: February 14, 2017; Published: March 1, 2017

This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This work was supported by National Institute of Health Grant (R01-GM110310 to JNO, R01- GM110310 to KYS, P50 GM103297 to TKD) and National Science Foundation (NSF) (DBI1040158 and NSF CHE1213668 to TKD). We acknowledge technical support and resources from the Center for Theoretical Biological Physics. Work at the Center for Theoretical Biological Physics was also supported by the NSF (Grant PHY-1427654 and NSF-CHE-1614101) and by the Cancer Prevention and Research Institute of Texas (CPRIT - grant R1110). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Non-coding RNAs are currently thought to account for over 75% of the human genome [1]. In bacteria, non-coding RNAs play important roles in gene regulation. One such class of RNAs, riboswitches, regulates metabolite production. Here, a single RNA sequence folds into one of two or more mutually exclusive folds depending on the metabolite concentration [2,3]. In some cases, such as the S-adenosylmethionine-I (SAM-I) riboswitch, the RNA contains a transcriptional terminator that forms when ligand is present, in effect silencing genes important for ligand production [4–6]. When the ligand is not present, the terminator does not form, allowing gene expression, and therefore ligand production, to continue efficiently. In other cases, such as the SAM-II riboswitch, ligand binding may lead to sequestration of the Shine/Dalgarno sequence, likely blocking ribosome binding and, as a consequence, protein synthesis [7]. While these examples of ligand-dependent secondary structure switches have been known for some time, a detailed thermodynamic understanding at the atomistic level, including the indispensable effect of the RNA’s ion-atmosphere, has not been achieved. In recent years, riboswitches have become canonical systems for studies of diverse RNA behaviors, as they possess quintessential characteristics of many RNA systems: ligand binding, Magnesium ion (Mg²⁺) sensitivity, conformational changes, secondary structure remodeling, and regulatory functions. Chemical footprinting, NMR, small-angle X-ray scattering (SAXS) and single molecule FRET techniques are being exploited to elucidate the folding kinetics, thermodynamics and the magnesium ion sensitivity in RNA systems such as the TPP riboswitch [8], glycine-dependent riboswitches [9], different variants of P4-P6 RNA [10–12] and P5abc subdomain of the Tetrahymena group I intron ribozyme [13]. Other work focuses on more complex functions, such as splicing and ligand recognition and their associations with proteins or different metabolites [14,15]. Success in understanding the structural, dynamical and functional aspects of riboswitch systems requires an integrated experimental and theoretical approach. Traditional crystallographic techniques produce static snapshots of the riboswitch. SmFRET, NMR, and SAXS methods obtain kinetic information and overall distributions of conformations. Molecular simulation allows one to integrate disparate experimental data into a single coherent picture, characterizing transitions in atomistic detail and the free energy landscape with fine resolution. A large number of riboswitches have been crystallized and have also been investigated via fluorescence and single molecule techniques [16–28].

Molecular simulations have also been used to study a number of riboswitches, including but not limited to the SAM-I, SAM-II, pre-Q, and adenine riboswitches [29–32]. Some of these studies replaced the essential atmosphere of divalent Mg²⁺ ions with other, monovalent ions. Others used a repulsive Debye-Hückel interaction between the phosphate groups of the RNA. Such implicit treatments of the ions, however, neglect important near-field effects that occur inside the core of the riboswitch, where the ions are strongly coupled to the RNA.

Although the role of Mg²⁺ in stabilizing the RNA tertiary structure has long been realized [33,34], the molecular basis of ion-RNA interactions, in terms of structure and function, is not well understood. In a pioneering study, Draper and co-workers distinguished three classes of ion environments: (i) the diffuse ions, which are not restrained to any particular region, (ii) water-surrounded ions separated from the RNA by a single hydration layer, which we call outer-sphere ions, and (iii) chelated ions in the inner sphere, which form direct contacts with at least two different phosphate groups of the RNA [34,35]. While the potential role of chelated ions has been emphasized in many studies, recent work using explicit solvent molecular dynamics instead highlights a dense layer of outer-sphere Mg²⁺ ions, which are primarily responsible for anchoring the RNA structure [36]. These outer-sphere Mg²⁺ ions are only transiently bound but nonetheless strongly coupled to the RNA dynamics. Such highly correlated Mg²⁺ ions may even reside in the core region of riboswitch RNAs [36,37]. This dynamic cloud of Mg²⁺ has also been simulated in all-atom reduced models that combine Manning theory with a background of monovalent ions, represented by Debye-Hückel interactions [37,38]. In addition to these native basin simulations, studies of metabolite recognition and specificity have also been initiated using conformational ensemble sampling, again in the absence of Mg²⁺ ions [39].

While most riboswitch studies have focused on riboswitches in the 5’-UTR of mRNA that control transcription, less attention has been paid to translational control by riboswitches through ligand-dependent sequestration of the ribosome binding site (i.e., the Shine/Dalgarno sequence). The SAM-II riboswitch is one such relatively small RNA element which regulates methionine and SAM biosynthesis. A single hairpin, classic H-type pseudoknot and triplex interaction near the ligand binding site make this RNA an interesting system to study RNA control over the translation initiation process [21,40–44].

A previous single molecule fluorescence resonance energy transfer (smFRET) [21] imaging study sheds light on the dynamic nature of ligand-free SAM-II riboswitch, which becomes conformationally restrained upon ligand binding. The flexibility of such highly transient conformations is tuned to ensure a viable time scale for conformational transitions in the absence of ligand. More rapid sensing could, however, be achieved if the riboswitch adopted a binding competent conformation in the ligand-free state. Mg²⁺ ions can act as effective anchors, aiding in the preservation of the structural integrity of the RNA. The emergence of two distinct FRET configurations in the presence of 2 mM Mg²⁺ in the ligand free system suggests that Mg²⁺ has the ability to compress the RNA structure, in such a way that it might pre-organize the RNA to form a binding competent conformation. In a series of small angle X-ray scattering (SAXS) experiments, we observed an analogous signature of Mg²⁺ induced structural collapse that can facilitate subsequent ligand binding [40]. Our studies provided direct insight into the global rearrangement induced by both Mg²⁺ and ligand. The compaction of RNA by Mg²⁺ was also studied by size-exclusion chromatography (SEC): changes in the measured elution volume suggested a decrease in the particles’ hydrodynamic radius [40].

Pre-organization by Mg²⁺ has also been observed in other riboswitches [21,28,45]. In the SAM-I system, our biochemical studies have shown that addition of Mg²⁺ yields the pre-organized partially folded state. In addition, we have shown that, in the absence of Mg²⁺, the fully folded state cannot be achieved, even at high ligand concentrations [23,28,46]. The presence of both, Mg²⁺ and ligand are required for the stabilization of the fully folded ligand-bound configuration.

While previous smFRET and SAXS data revealed that ligand free RNA undergoes substantial structural changes upon variation of Mg²⁺ concentration, these structural changes often remained undetected by traditional NMR and X-ray crystallography techniques because of the transient nature and low population levels of such intermediates [47]. The newly developed Chemical Exchange Saturation Transfer (CEST) measurements are now capable of probing these sparsely and transiently populated RNA conformations. Earlier we studied NMR dynamics of the SAM-II system with this new method [47]. The data indeed confirmed that SAM-II riboswitch can access a sparsely populated but bound-like pre-organized state even in the absence of ligand [40,47].

In the present study, we performed molecular simulations to predict the effect of Mg²⁺ on the conformational landscape of the SAM-II riboswitch. We then tested these predictions with ¹³C-CEST data. Analysis of our simulations yields the free energy landscape of the SAM-II riboswitch, the effect of Mg²⁺ on this landscape and insight into the microscopic origins of these effects. More specifically, reappraisal of the ¹³C-CEST data for the ligand-free SAM-II riboswitch at different Mg²⁺ concentrations enabled us to probe the influence of Mg²⁺ on sparsely populated bound-like pre-organized states. We then revisited earlier smFRET, SAXS and SEC elution profiles and compared with our present equilibrium simulation results to integrate these data into a unified scenario of Mg²⁺-induced collapse. We calculated the free energy landscape of the SAM-II riboswitch using our recently developed all-atom structure-based model (SBM) that includes explicit Mg²⁺ ions and the effects of Manning condensation and Debye-Hückel Potassium and Chloride interactions. We specifically predict that, as Mg²⁺ concentration is increased from 0.25 mM to 2 mM, the SAM-II riboswitch collapses from an extended, partially unfolded state to a highly compact, pre-organized state, in agreement with the ¹³C-CEST studies, where we observe a shift in population towards a bound-like conformation. In addition, our simulations characterize this collapse transition in terms of the radius of gyration as a function of Mg²⁺ concentration, which is qualitatively similar to previous SAXS measurements. This agreement gives us confidence in the microscopic details of our simulations, showing that the triplex formation between helix P2b and loop L1 plays an important role in the collapse process.

Results

Simulations of the SAM-II riboswitch were performed with and without SAM at different Mg²⁺ concentrations. The simulation started from the crystal structure of SAM-II riboswitch (pdb accession code: 2QWY) [20]. A global view of this structure and details of the secondary structure are presented in Fig 1. As ligand-free conformations of SAM-II are still under investigation in experiments, we explore the full folding free energy landscape both in the presence of ligand (SAM) and in the absence of ligand to characterize the entire accessible conformational space. To cover the folding landscape, we employed umbrella sampling over the reaction coordinate, Q, which is the fraction of intra-molecular native contacts, present in both the free and bound states, formed by the riboswitch.

Download:

Fig 1. Secondary and tertiary structure of SAM-II riboswitch in ligand-bound state (pdb:2QWY).

(a) Sequence-aligned secondary structure of the SAM-II riboswitch where base pair and stacking interactions are indicated. There are three helices and two strands highlighted with different colors. (P1: Orange, P2a: Red, P2b: Blue, L1: Green, L3: Magenta). (b) Tertiary structure displays triple helix between helix P2b and loop L1.

https://doi.org/10.1371/journal.pcbi.1005406.g001

Detailed agreement between ¹³C-CEST, smFRET, SAXS, SEC elution profile and generalized Manning model corroborates pre-organization by Mg²⁺

As mentioned earlier, CEST experiments are able to capture transiently populated dynamic conformations [13]. This strategy was applied to the ligand-free SAM-II riboswitch in the presence of 0.25 mM and 2 mM Mg²⁺. The ¹³C-CEST profiles of the labeled ribose C1’ and base C6 carbons of C43 were recorded at three different B₁-fields (17.5, 27.9, 37.8 Hz) with a mixing time of 0.3 s at 298 K [47]. We compare the data for 0.25 mM and 2 mM Mg²⁺ concentrations at B₁-field of 17.5 Hz (Fig 2a). The data were fit with a two-state model where the low population of the partially closed state (peaks around 300 Hz spin-lock offset) appears to increase with addition of 2 mM Mg²⁺. Consistent results were obtained from CEST profiles for other B₁-fields of 27.9 and 37.8 Hz (Fig S1 in the S1 Text). Trajectory plots of the fraction of native contacts extracted from the generalized Manning equilibrium simulations of ligand-free SAM-II at these two concentrations clearly show the hopping between different conformations (Fig 2b). Furthermore, the dynamic transitions between the two major states (bound-like: Q≈0.9 and open: Q≈0.7) visit native-like conformations (Fig 2c) more frequently at 2 mM than at 0.25 mM Mg²⁺, as summarized in the corresponding contact histograms, P(Q) (Fig 2d). Both CEST experiments and simulation data indicate that the equilibrium shifts from the open conformations toward the native bound-like state as we increase Mg²⁺ concentration.

Download:

Fig 2. Detailed agreement between ¹³C CEST, smFRET, SAXS, and SEC elution profiles with the generalized Manning model corroborates pre-organization by Mg²⁺.

(a)-(d) Comparison between CEST-NMR experiments and Manning model simulations. (a) ¹³C CEST profiles for C43-C6 of SAM-II riboswitch in the ligand-free situation in the presence of two different concentrations of Mg²⁺ (black dots: 0.25 mM, green dots: 2.0 mM) at B₁ field of 17.5 Hz. The population distribution involving two major states was obtained by fitting each ¹³C CEST profile (solid lines) into a two-state model. (b) The fraction of native contacts (Q) dynamics extracted from simulations indicate the transition between two major states at 0.25 mM. (c) Same as (b) for 2.0 mM Mg²⁺. (d) The fraction of native contact population histograms at these two concentrations. (e)-(f) Steady state smFRET analysis from fluorescently labeled experiments and theoretical prediction from simulations. (e) Experimental smFRET occupancy histogram in the presence of 2.0 mM Mg²⁺ (adapted from ref. 22) [22]. (f) Theoretical prediction of FRET occupancy in the presence of 2.0 mM Mg²⁺. (g) Superposition of a set of different conformations at 2.0 mM Mg²⁺ depicts the prevalence of two distinct sets of ensembles. (h) Comparison of experimental scattering profiles of SAM-II collected from four different buffer conditions. Theoretical SAXS predictions are depicted by solid lines and the experimental data by dots (experimental data adapted from ref. 40). (i) Compaction of SAM-II as a function of [Mg²⁺] obtained from SEC elution profiles is consistent with a decrease in hydrodynamic radius (experimental data adapted from ref. 40), [40] in qualitative agreement with the observed folding transition in the average radius of gyration (Rg) obtained from equilibrium simulations at Mg_1/2≈6.0 mM.

https://doi.org/10.1371/journal.pcbi.1005406.g002

The signature of the existence of such Mg²⁺ induced bound-like states has also been reported in previous smFRET experiments (Fig 2e) [22]. To support our observations we have revisited some of these smFRET efficiency assessments [22] and compared them with theoretical FRET predictions obtained from our generalized Manning model simulations under similar buffer conditions. The equation used for theoretical FRET prediction is described in section S1 in the S1 Text. We tracked the dynamics of positions 14 and 52, where acceptor (cy5) and donor (cy3) fluorophore labels were placed in the smFRET experiments (Fig 2f). Both experimental and simulation FRET confirm the coexistence of two states at 2 mM Mg²⁺ (Fig 2g) [22].

Previous SAXS data corroborates well the existing smFRET observations [22,40]. The SAXS data also indicated both ligand and Mg²⁺ ions are required to effectively fold this riboswitch. To microscopically understand their mutual and stand-alone effects from the present simulations, we studied the conformational differences of this riboswitch in four extreme buffer conditions and compared our computational results with experimental SAXS data. For this comparison, we extracted multiple snapshots from several long trajectories and computed ensemble averaged SAXS profiles using the Debye formula for spherical scatterers parameterized in the FoXS web server [48,49] as described in section S2 in the S1 Text. The predicted SAXS curves here show qualitative agreement with experiments (Fig 2h) [40]. The Kratky representation of SAXS data presented in Fig s2 in the S1 Text shows a pronounced peak, indicating the emergence of more extended conformations with decreasing Mg²⁺ concentrations. We note that capturing the entire conformational heterogeneity of an extended state is computationally challenging. This mostly applies for the extreme case where neither ligand nor Mg²⁺ is present. In this case, the correlation between theoretical and experimental SAXS profiles leaves room for improvement. Values for chi-square reflect that and are shown in Table S1 in the S1 Text. These analyses indeed suggest the potential impact of both, ligand and Mg²⁺, stabilizing the closed conformations, which we characterize further below with contact data to describe the pre-organization and the ligand-organized closing.

A significant Mg²⁺ induced collapse transition, as indicated by the SAXS data, has been followed over a wide concentration range (up to 100 mM) of Mg²⁺ in SEC elution volume profile (Fig 2i). Here RNA elutes after longer retention times with increasing Mg²⁺ concentration ([Mg²⁺]) in the mobile phase [40]. Bigger elution volume signifies decreasing hydrodynamic radius of a monomeric RNA molecule. The folding transitions, both from experimental elution volume data and from average Rg measured from the equilibrium simulation analysis as functions of [Mg²⁺], follow sigmoid curves with transition midpoint, Mg_1/2 at 6 mM (Fig 2i).

Mg²⁺ remodels the free energy landscape, favoring pre-organized states

At this point, a range of experimental techniques and simulation data support the existence of pre-organized states. Here we aim to obtain a thermodynamic description of how Mg²⁺ governs the energy landscape of RNA from our model simulation study. In Fig 3a, we show the free energy landscape for the folding transition of SAM-II riboswitch in SAM-bound (in the presence of explicit ligand) and SAM-free (in the absence of ligand) conditions near the physiological concentration of Mg²⁺ ([Mg²⁺] = 2.0 mM). During this folding transition, each secondary structural segment folds sequentially illuminating the pathway of folding (Fig 3b). The free energy profile, in the presence of explicit SAM has a distinct bound-state-well, reflecting the ligand-induced stabilization of the closed conformations (designated as (i) in Fig 3c). In the apo-form of the riboswitch, the fully closed bound state does not correspond to a minimum in the landscape. At lower Q than this ligand-bound state, the free energy profile for apo-SAM-II riboswitch reveals three distinct minima. They involve: a ligand-free partially closed state (state (ii) in Fig 3c), which has a substantial overlap with the ligand-bound closed conformation. In this state, the nonlocal contacts (involving base-pairing contacts) including base-stacking contacts in P1, the P1-L3 pseudo-knot interaction, and major segments of P2b and the L1-P2b triplex interactions remain secured, while the contacts involved in Shine/Dalgarno sequence (AAAG50G51A523´), and in the part of L1-P2b are disrupted (state (ii) in Fig 3c).

Download:

Fig 3. Magnesium dependence of the SAM-II riboswitch free energy landscape demonstrates Mg²⁺ pre-organization of bound-like state.

(a) The free energy landscape as a function of the fraction of native contacts (Q) of SAM-II near physiological Mg²⁺ concentration ([Mg²⁺] = 2.0 mM). The system explores three distinct barrier-separated minima on the free energy landscape. (b) The order of secondary structure formation as a function of the fraction of all native contacts (Q) as measured by the fraction of non-local regional contacts (Q_Regional). The transition displays high cooperativity. (c) Representative structures. The representative structure corresponding to each region of the energy landscape is designated as follows: (i) ligand-bound closed (C), (ii) ligand-free partially closed (PC) including triplex interactions, (iii) ligand-free partially open (PO), (iv) ligand-free open (O), and (v) the unfolded (U) state. Green/mauve, bases in native conformations. The ligand-bound closed minimum has been characterized from the free energy profile of the folding transition of SAM-bound RNA which indicates the ligand-free partially closed conformation has a substantial resemblance with the ligand-bound closed state. (d) Mg²⁺ concentration dependence of the folding transition of apo-SAM-II riboswitch. Average free energy profiles of folding transition at different [Mg²⁺] show significant stability difference between the ligand-free partially closed (PC) and the partially open (PO) minima (ΔG_PC-PO); and the unfolded (U) and the extended open (O) state minima (ΔG_U-O). In the inset ΔG_PC-PO and ΔG_U-O are plotted together as a function of [Mg²⁺]. Both sigmoid curves follow the same transition midpoint, Mg_1/2 ≈6.0 mM, as found in the SEC profile.

https://doi.org/10.1371/journal.pcbi.1005406.g003

Recent fluorescence and NMR spectroscopic data also indicated that C16 in P2a helix remains mostly unpaired in the absence of SAM [22]. The data also suggested that formation of the pseudoknot in the absence of SAM is highly transient in nature. Intermediate states, (iii) and (iv) in Fig 3c, although marginally separated by a small barrier, effectively belong to a broad, flat basin which involves an ensemble of partially folded open configurations. A representative unfolded structure ((v) in Fig 3c) is shown to describe the unfolded minimum. As we increase the concentration of Mg²⁺ we find enhanced stabilization of the pre-organized partially closed conformations (state (ii) in Fig 3d) relative to the open conformations. Our latest ¹³C-CEST chemical exchange data anticipated that the emergence of Mg-induced pre-organization can have immense consequences for rapid ligand recognition [47]. In the context of the simulation, Mg²⁺ induced thermodynamic stabilization is reflected by the difference in stability, ΔG_PC-PO between the bound-like partially closed (PC) conformation and the partially open (PO) conformation and by ΔG_U-O between the unfolded (U) and the open conformation (O). These two stability differences, ΔG_PC-PO and ΔG_U-O vary upon increasing [Mg²⁺] until they reach their saturation limits. Both ΔG_PC-PO and ΔG_U-O plotted as functions of [Mg²⁺], are fitted well to sigmoid curves with a Mg²⁺_1/2 value around 6 mM (inset of Fig 3d) which again correlates well with the SEC elution volume data (Fig 2i) [40].

Mg²⁺-induced phosphate interactions facilitate pre-organization

To address the open question of how Mg²⁺ ions regulate structural collapse, we have determined the Mg²⁺ distribution in the ion-solvation layer of SAM-II, which accommodates increasing numbers of Mg²⁺ up to 8 mM Mg²⁺ content (Fig 4a). Subsequent additions of Mg²⁺ beyond 8 mM do not effectively add to the 1st layer of Mg²⁺ solvation. How we characterize the ion-solvation layer from our simulated trajectories is described in section S3 in the S1 Text (Fig s3 in S1 Text). We have further classified the outer sphere Mg²⁺ present in the ion-solvation layer into two categories based on their number of associated phosphate groups: (i) Single phosphate coordinating Mg²⁺ (Fig 4b), which efficiently neutralize the negative charge of the adjacent phosphate (Fig 4d), and (ii) multiple phosphate coordinating Mg²⁺ (Fig 4c). The key role in stabilizing the structure is played by such Mg²⁺ bridging multiple phosphates, which can act as glue in compact structures by holding a number of negatively charged phosphates together in close proximity.

Download:

Fig 4. Probability distribution and free energy landscape of the SAM-II riboswitch as a function of the number of Mg²⁺ions in the ion solvation layer and Mg²⁺ mediated phosphate contacts (PHOS_Cont).

(a) The distribution of Mg²⁺ in the ion solvation layer shows gradual shift accommodating more number of Mg²⁺ with increasing Mg²⁺ concentration. Beyond 8.0 mM the first ion-solvation layer appears saturated. Subsequent additions of Mg²⁺ beyond 8.0 mM do not effectively participate in the 1st layer of solvation by ionic interaction. (b) Snapshot extracted from the simulation at buffer condition [Mg²⁺] = 0.4 mM has single phosphate coordinated Mg²⁺. (c) Snapshot extracted from the simulation at buffer condition [Mg²⁺] = 8.0 mM has multiple phosphate coordinated Mg²⁺. (d) Single phosphate coordinated Mg²⁺ as a function of Mg concentration. (e) Multiple phosphate coordinated Mg²⁺ as a function of [Mg²⁺]. The population shift of such multiple coordinated Mg²⁺ serves to connect negatively charged phosphate groups. (f) Free energy landscape as a function of N_cont and PHOS_Cont for [Mg²⁺] = 0.4 mM. (g) Same as (f) for [Mg²⁺] = 2.0 mM. (h) Same as (f) for [Mg²⁺] = 8.0 mM. Increasing [Mg²⁺] reshapes the landscape, shifting the equilibrium from a basin at lower N_cont to a basin at higher N_cont.

https://doi.org/10.1371/journal.pcbi.1005406.g004

The population shift coincident with multiple coordinated Mg²⁺ ions with increasing Mg²⁺ concentration directly supports their role in stabilizing the structure (Fig 4e). We have also investigated the thermodynamic impact of Mg²⁺-mediated phosphate contacts (PHOS_Cont: total number of pair-wise phosphate-phosphate contacts) on the energy landscape as a function of overall folding progress, expressed by the number of native contacts (N_Cont), as shown in Fig 4f–4h. As we increase the Mg²⁺ concentration the broad minimum that appeared around N_Cont~800, involving partially folded open conformations, gradually becomes more stabilized. Concurrent enrichment of phosphate-phosphate contacts extends the contour of the minimum asymmetrically toward higher PHOS_Cont. Additionally, by 8 mM [Mg²⁺], the bound-like pre-organized state grows with substantial population, stabilized again by phosphate connections (Fig 4h).

Mg²⁺-induced triplex interactions between helix P2b and L1 loop help pre-organize the SAM-II riboswitch

We analyzed long equilibrium trajectories of the apo- and bound-forms of SAM-II slightly below the folding temperature in order to capture the essential characteristics of the pre-organized state, and also to compare this state with the fully folded ligand bound state. We have evaluated the distribution of native contact formation in each segment of secondary structure as a function of the total Q at different Mg²⁺ concentrations. Plots show contact formation in P2b (Fig 5a–5d) and the triplex interaction between helix P2b and loop L1 (Fig 5e–5h), which are most affected by Mg²⁺ concentration. Data for the nonlocal contacts of P1, L3-P1, P2a, which appear only marginally affected by Mg²⁺ concentration, are shown in Fig s4 in S1 Text. The two distinct basins visible at low [Mg²⁺], for the P2b helix and L1-P2b triplex contacts correspond to the pre-organized (at higher Q) and open states (at lower Q). At increasing [Mg²⁺], the populations gradually shift towards the pre-organized state.

Download:

Fig 5. Thermodynamics of triplex formation.

Free energy landscape of the SAM-II riboswitch as a function of total native contacts (x-axis) and native contacts for an individual structural element (y-axis). (a)-(c) Free energy landscapes for helix P2b formation for increasing values of [Mg²⁺], without SAM. (d) Free energy landscapes for helix P2b formation for the case of [Mg²⁺] = 2.0 mM, with SAM. (e)-(g) Free energy landscapes for triplex formation (L1-P2b contacts) for increasing values of [Mg²⁺], without SAM. (h) Free energy landscapes for triplex formation for the case of [Mg²⁺] = 2.0 mM, with SAM. (a)-(c) and (e)-(g) Collective integrity of triplex interaction (involving both P2b and P2b-L1) appears most sensitive RNA element to [Mg²⁺].

https://doi.org/10.1371/journal.pcbi.1005406.g005

Around 8 mM Mg²⁺, the dominant contribution arising from this pre-organized triplex to the conformational space is evident from Fig 5c and 5g. Ligand binding also strongly favors structure formation, even at moderate [Mg²⁺], as the ligand bridges the gap between L1 strand and P2b helix, producing the fully formed triplex.

Discussion

Motivated by our ¹³C-CEST profiles for SAM-II and their [Mg²⁺] dependence we have explored the free energy landscape of the SAM-II riboswitch using a recently developed all-atom SBM that includes explicit Mg²⁺ ions, Debye-Hückel treatment of implicit KCl interactions, and the effects of Manning condensation to accurately account for the ion atmosphere around the RNA. Our results support a mechanism involving Mg²⁺ induced pre-organization followed by conformational selection by the ligand, SAM, as we speculated in an early study [47]. The free energy analysis validates the observations of that pre-organization, providing an atomistic and thermodynamic basis for the enhanced population of a partially collapsed, pre-organized ensemble at sufficiently high Mg²⁺ concentration in the absence of ligand. We observe three distinct sets of conformations in the folding free energy landscape of ligand-free SAM-II riboswitch: (i) an ensemble of unfolded conformations, (ii) a broad ensemble of partially folded open conformations, and (iii) an ensemble of pre-organized bound-like conformations. As we increase magnesium concentration beyond 2–4 mM, the bound-like ensemble is further stabilized, shifting the equilibrium toward the pre-organized states. All the experimental results from our ¹³C-CEST profile, recent SAXS, single molecule FRET, and size-exclusion chromatographic studies are assembled and found to be in good agreement with the present simulation results (Fig 2). At higher concentrations, Mg²⁺ stabilizes compact structures by coordinating multiple charged phosphate groups of RNA in close proximity. The experimental results, together with free energy landscapes confirm that sufficient Mg²⁺ can indeed promote stable ligand binding in the SAM-II riboswitch, and is likely the structural basis for the switching control of protein translation.

While this structural pre-organization of SAM-II can assist in rapid ligand recognition, our study suggests that a sufficiently high concentration of Mg²⁺ is necessary to capture those pre-organized states. Only when the system achieves a well-organized ion solvation layer at high [Mg²⁺], the effect of additional Mg²⁺ seems limited. This layer involves a number of Mg²⁺ ions, each coordinating with multiple phosphate groups. Mg²⁺ ions thus serve as glue to the negatively charged phosphates and facilitate the structural compaction. We note that while chelated Mg²⁺ may play an important role in other riboswitch RNAs, no specific chelated ions have been reported so far in the SAM-II system.

Our molecular simulation trajectories also allow us to pinpoint the structural basis of the effect, revealing that triplex interaction between the helix P2b and its association with the L1 strand dominate the process of pre-organization as summarized in Fig 5a–5c and 5e–5g, showing the gain of structure with increasing [Mg²⁺]. In the final step, ligand binding firmly bridges the extended gap between L1 and P2b, which seems otherwise not achievable through the addition of small, dynamic Mg²⁺ alone. But although P2b and its connection with L1 can be secured by the ligand, its presence again alone cannot fully stabilize the overall structure without addition of significant amount of Mg²⁺ (Fig 2h). These findings suggest that a sufficiently high concentration of Mg²⁺ is necessary to stabilize the pre-organized triplex and then the presence ligand promotes the native triplex formation, as summarized in Fig 5d and 5h. We note that triplexes have recently emerged as important players in gene regulation by non-coding RNAs [50–53]. Base triples also play a role in RNase P and the Diels-Alder ribozyme [54]. Heroic calculations, as such recent microsecond explicit solvent simulations of riboswitches, will also shed light on these effects, especially regarding the role of solvation [55].

Nucleic acid-ion interactions make a substantial energetic contribution in the stabilization of the native state of RNAs, including complex formation with proteins and other macromolecules [56]. The dynamics of nucleic acids are also found to be strongly influenced by the motion of their ion atmospheres. Relative to other ionic species, Mg²⁺ can efficiently support a close assembly of negatively charged phosphates by mediating favorable interactions among them. Other earth alkali metals/divalent ions (e.g. Ca²⁺) and even monovalent ions are also able to induce similar transitions, albeit at higher concentration. Our early SEC elution profiles for SAM-II show that the transition midpoint in presence of Potassium (K⁺) alone occurs only at [K]_1/2 ≈ 25 mM. The midpoint for Calcium (Ca²⁺) is [Ca]_1/2 ≈ 8 mM, compared to 6 mM for the Mg²⁺ ion [40]. This is a direct result of the larger charge/radius ratio of magnesium [40,57].

Thus, having these special characteristics, Mg²⁺ efficiently helps pre-organize the system and enables access to the partially collapsed states that are further stabilized by ligand binding. The general importance of Mg²⁺ for the stability of compact RNA structures supports a possibly universal role of conformational selection in ligand-binding RNAs, such as riboswitches, aptamers, and possibly protein-binding RNAs. A detailed thermodynamic understanding of the underlying landscape will indeed enable greater control of riboswitch regulation, highly sought after by researchers in synthetic biology who are currently employing riboswitches as ligand-dependent ‘knobs’ to control desired gene expression [58].

Methods

RNA electrostatic model

Our all-atom structure-based model (SBM) has proven successful in describing the dynamics of numerous proteins and macromolecular complexes [59–63]. To elucidate RNA free energy landscapes under the influence of Mg²⁺, models capable of quantitatively describing the ion atmosphere are needed, including ionic condensation around the negatively charged phosphate groups of RNA. Early studies have simply included electrostatic effects in SBM of RNA via repulsive Debye-Hückel interactions, thus treating all ions implicitly [29,30]. Recently, our group developed a more detailed model of RNA electrostatics and applied it within all-atom structure-based molecular dynamics simulations. Our model treats Mg²⁺ ions explicitly to account for ion-ion correlations neglected by mean-field theories [38]. The KCl buffer, which completes the experimental setup, is treated implicitly by a generalized Manning counter ion condensation model [38,64], since mean-field theories correctly assess the charge densities of monovalent K⁺ and Cl^- ions. Classical Manning counter-ion condensation theory was originally developed for understanding the low concentration limiting behavior of polyelectrolyte chains modeled with an infinite line of charge. Folded RNA, however, is not a line of charge. To account for the compact and irregular structures of RNAs and the effects of varying ion concentrations, we improve the Manning counter ion condensation model to handle electrostatic heterogeneity, making the condensed charge density a dynamical function of each phosphate coordinate. KCl screening is characterized by a Debye-Hückel potential. Removal of the continuum screening ions from the inaccessible volume of RNA is a substantial extension to Manning counter-ion condensation. The model has been tested against experimental measurements of excess Mg²⁺ associated with RNA, characterizing the Mg²⁺-RNA interaction free energy. This hybrid SBM has opened up new possibilities to study various structural and functional processes of RNA that are essentially controlled by ions [38]. In the present study we used this recently developed all-atom hybrid SBM to understand the conformational transition of SAM-II and the corresponding Mg²⁺ sensitivity. The energy function used in this model is given below, (1) where, Φ_SBM is the all-atom SBM potential ensuring a global minimum in the landscape for the native state of RNA. The SBM potential is composed of two general types of interactions: (2) where, Φ_local characterizes the local interactions that encode covalent bonds and torsional angles, maintaining the correct local geometry and chirality. Φ_non-local comprises two non-local contributions: (i) an attractive term that is applied specifically to all tertiary interactions determined from the native structure, (ii) the general repulsive interactions, that describe the excluded volume by symmetric hard potentials (to avoid any unwanted chain crossing). Φ_Mg-Size adds the excluded volume interactions involving the explicit Mg²⁺ ions, regulating RNA-Mg²⁺ and Mg²⁺-Mg²⁺ interactions. Φ_ion-effect accounts for all interactions between charges in the system which consist of the fixed charge distribution of the RNA and the dynamic contribution from the ions. Mg²⁺ and phosphate charges interact via a Debye-Hückel potential with a screening term that depends, in turn, on the distribution of the monovalent ions. The monovalent ions, K⁺ and Cl^- from the added salt, fall into two categories: screening ions and Manning condensed ions. The screening ion density is obtained using Debye-Hückel electrostatics. The density of the Manning-condensed ions is modeled as the sum of two normalized Gaussian distributions where the center of each Gaussian is located on the position of the negatively charged phosphate group. All the condensation variables along with the explicit Mg²⁺ and RNA coordinates are evolved with Langevin dynamics [38]. The mathematical formulations of all the terms and the related parameterizations are discussed in depth in section S4 in the S1 Text.

Free energy simulations

The umbrella sampling method [65] was used to sample the conformational space of SAM-II riboswitch along the reaction coordinate, Q, which is the fraction of intra-molecular native contacts in the riboswitch. The Weighted Histogram Analysis Method [66] was then used to calculate the thermodynamic quantity, G(Q). The detail is described in section S5 in the S1 Text.

¹³C-CEST experiments

CEST data were collected using a pseudo-3D HSQC experiment with the B₁ field offsets (-600 to 600 Hz) incremented in an interleaved manner with 3 references (no CEST period) [47]. A total of 1024x16 complex points were recorded [40] with 32 transients with a recovery delay of 1.5 s for a total experimental time of approximately 12 hr for each spin-lock field. A CEST saturation period of 100 ms was used for the base and 200 ms for ribose. The pulse program used was an adaptation of a previously published one without the need for selective pulses [47].

Analysis of relaxation parameters

We used a two-state model to fit each of the three profiles of the selectively labeled carbon (ribose C1’ and base C6) and quantitatively extracted the carbon chemical shift (Δω), the exchange rate, and the population of the minor state based on the Bloch−McConnell 7x7 matrix [47]. The CEST data was plotted as I(t)/I(0) versus spin-lock offset (Hz) and was fit by numerically solving the matrix exponential for the CEST spin-lock period based on this 7x7 two-state Bloch-McConnell equation as described earlier [47,67].

FRET

In experiments, Fluorescence Resonance Energy Transfer (FRET) efficiency is the quantum yield of the energy transfer where a donor chromophore from its excited electronic state may transfer its energy to an acceptor chromophore through a non-radiative dipole-dipole coupling. The FRET efficiency varies with the separation between donor and acceptor fluorophores following the Fӧrster relation. For theoretical FRET predictions we use the Fӧrster relation where the value of Fӧrster radius is taken as 53Å [68]. We described it in detail in section S1 in the S1 Text.

SAXS analysis

In SAXS experiments, the scattering intensity is measured from the electron density difference between the purified sample and that of the solvent/buffer. FoXs is a method that uses the Debye formula by which a theoretical scattering profile of a structure can be computed [48,49]. The detail is discussed in section S2 in the S1 Text.

Supporting information

S1 Text. Supporting figures, table and text.

https://doi.org/10.1371/journal.pcbi.1005406.s001

(PDF)

Author Contributions

Conceptualization: JNO KYS SR.
Data curation: SR BC RL.
Formal analysis: SR JNO KYS.
Funding acquisition: JNO KYS TKD.
Investigation: SR BC RL.
Methodology: SR RLH HL.
Project administration: JNO KYS.
Resources: JNO KYS.
Supervision: JNO TKD KYS.
Validation: SR BC RL TKD.
Visualization: SR JNO KYS.
Writing – original draft: SR KYS TKD.
Writing – review & editing: SR HL JNO KYS TKD.

References

1. Huarte M (2015) The emerging role of lncRNAs in cancer. Nat Med 21: 1253–1261. pmid:26540387
- View Article
- PubMed/NCBI
- Google Scholar
2. Montange RK, Batey RT (2008) Riboswitches: Emerging Themes in RNA Structure and Function. Annu Rev Biophys 37: 117–133. pmid:18573075
- View Article
- PubMed/NCBI
- Google Scholar
3. Roth A, Breaker RR (2009) The Structural and Functional Diversity of Metabolite-Binding Riboswitches. Annu Rev Biochem 78: 305–334. pmid:19298181
- View Article
- PubMed/NCBI
- Google Scholar
4. Lin JC, Yoon J, Hyeon C, Thirumalai D (2015) Using simulations and kinetic network models to reveal the dynamics and functions of riboswitches. Methods Enzymol 553: 235–258. pmid:25726468
- View Article
- PubMed/NCBI
- Google Scholar
5. Serganov A, Patel DJ (2007) Ribozymes, riboswitches and beyond: regulation of gene expression without proteins. Nature Reviews Genetics 8: 776–790. pmid:17846637
- View Article
- PubMed/NCBI
- Google Scholar
6. Winkler WC, Nahvi A, Sudarsan N, Barrick JE, Breaker RR (2003) An mRNA structure that controls gene expression by binding S-adenosylmethionine. Nature Structural Biology 10: 701–707. pmid:12910260
- View Article
- PubMed/NCBI
- Google Scholar
7. Corbino KA, Barrick JE, Lim J, Welz R, Tucker BJ, et al. (2005) Genome Biology 6: R70. pmid:16086852
- View Article
- PubMed/NCBI
- Google Scholar
8. Yamauchi T, Miyoshi D, Kubodera T, Nishimura A, Nakai S, et al. (2005) Roles of Mg2+ in TPP-dependent riboswitch. FEBS Lett 579: 2583–2588. pmid:15862294
- View Article
- PubMed/NCBI
- Google Scholar
9. Lipfert J, Das R, Chu VB, Kudaravalli M, Boyd N, et al. (2007) Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae. J Mol Biol 365: 1393–1406. pmid:17118400
- View Article
- PubMed/NCBI
- Google Scholar
10. Bisaria N, Greenfeld M, Limouse C, Pavlichin DS, Mabuchi H, et al. (2016) Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway. Proc Natl Acad Sci U S A 113: E4956–4965. pmid:27493222
- View Article
- PubMed/NCBI
- Google Scholar
11. Bisaria N, Herschlag D (2015) Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET. Biochem Soc Trans 43: 172–178. pmid:25849913
- View Article
- PubMed/NCBI
- Google Scholar
12. Takamoto K, Das R, He Q, Doniach S, Brenowitz M, et al. (2004) Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations. J Mol Biol 343: 1195–1206. pmid:15491606
- View Article
- PubMed/NCBI
- Google Scholar
13. Xue Y, Gracia B, Herschlag D, Russell R, Al-Hashimi HM (2016) Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch. Nat Commun 7: ncomms11768. pmid:27292179
- View Article
- PubMed/NCBI
- Google Scholar
14. Bocobza SE, Aharoni A (2014) Small molecules that interact with RNA: riboswitch-based gene control and its involvement in metabolic regulation in plants and algae. Plant J 79: 693–703. pmid:24773387
- View Article
- PubMed/NCBI
- Google Scholar
15. Russell R (2015) Unwinding the mechanisms of a DEAD-box RNA helicase in cancer. J Mol Biol 427: 1797–1800. pmid:25836982
- View Article
- PubMed/NCBI
- Google Scholar
16. Batey RT, Gilbert SD, Montange RK (2004) Structure of a natural guanine-responsive riboswitch complexed with the metabolite hypoxanthine. Nature 432: 411–415. pmid:15549109
- View Article
- PubMed/NCBI
- Google Scholar
17. Edwards AL, Reyes FE, Heroux A, Batey RT (2010) Structural basis for recognition of S-adenosylhomocysteine by riboswitches. RNA 16: 2144–2155. pmid:20864509
- View Article
- PubMed/NCBI
- Google Scholar
18. Garst AD, Heroux A, Rambo RP, Batey RT (2008) Crystal structure of the lysine riboswitch regulatory mRNA element. J Biol Chem 283: 22347–22351. pmid:18593706
- View Article
- PubMed/NCBI
- Google Scholar
19. Gilbert SD, Montange RK, Stoddard CD, Batey RT (2006) Structural Studies of the Purine and SAM Binding Riboswitches. Cold Spring Harbor Symposia on Quantitative Biology 71: 259–268. pmid:17381305
- View Article
- PubMed/NCBI
- Google Scholar
20. Gilbert SD, Rambo RP, Van Tyne D, Batey RT (2008) Structure of the SAM-II riboswitch bound to S-adenosylmethionine. Nat Struct Mol Biol 15: 177–182. pmid:18204466
- View Article
- PubMed/NCBI
- Google Scholar
21. Haller A, Altman RB, Souliere MF, Blanchard SC, Micura R (2013) Folding and ligand recognition of the TPP riboswitch aptamer at single-molecule resolution. Proceedings of the National Academy of Sciences 110: 4188–4193.
- View Article
- Google Scholar
22. Haller A, Rieder U, Aigner M, Blanchard SC, Micura R (2011) Conformational capture of the SAM-II riboswitch. Nature Chemical Biology 7: 393–400. pmid:21532598
- View Article
- PubMed/NCBI
- Google Scholar
23. Hennelly SP, Sanbonmatsu KY (2010) Tertiary contacts control switching of the SAM-I riboswitch. Nucleic Acids Research 39: 2416–2431. pmid:21097777
- View Article
- PubMed/NCBI
- Google Scholar
24. Heppell B, Mulhbacher J, Penedo JC, Lafontaine DA (2009) Application of Fluorescent Measurements for Characterization of Riboswitch–Ligand Interactions. Methods in Molecular Biology: Springer Science + Business Media. pp. 25–37.
25. Lemay J-F, Penedo JC, Tremblay R, Lilley DMJ, Lafontaine Daniel A (2006) Folding of the Adenine Riboswitch. Chemistry & Biology 13: 857–868.
- View Article
- Google Scholar
26. Montange RK, Batey RT (2006) Structure of the S-adenosylmethionine riboswitch regulatory mRNA element. Nature 441: 1172–1175. pmid:16810258
- View Article
- PubMed/NCBI
- Google Scholar
27. Stoddard CD, Montange RK, Hennelly SP, Rambo RP, Sanbonmatsu KY, et al. (2010) Free State Conformational Sampling of the SAM-I Riboswitch Aptamer Domain. Structure 18: 787–797. pmid:20637415
- View Article
- PubMed/NCBI
- Google Scholar
28. Wacker A, Buck J, Mathieu D, Richter C, Wohnert J, et al. (2011) Structure and dynamics of the deoxyguanosine-sensing riboswitch studied by NMR-spectroscopy. Nucleic Acids Research 39: 6802–6812. pmid:21576236
- View Article
- PubMed/NCBI
- Google Scholar
29. Cho SS, Pincus DL, Thirumalai D (2009) Assembly mechanisms of RNA pseudoknots are determined by the stabilities of constituent secondary structures. Proceedings of the National Academy of Sciences 106: 17349–17354.
- View Article
- Google Scholar
30. Hyeon C, Thirumalai D (2005) Mechanical unfolding of RNA hairpins. Proceedings of the National Academy of Sciences 102: 6789–6794.
- View Article
- Google Scholar
31. Whitford PC, Schug A, Saunders J, Hennelly SP, Onuchic JN, et al. (2009) Nonlocal Helix Formation Is Key to Understanding S-Adenosylmethionine-1 Riboswitch Function. Biophysical Journal 96: L7–L9. pmid:19167285
- View Article
- PubMed/NCBI
- Google Scholar
32. Xue X, Yongjun W, Zhihong L (2015) Folding of SAM-II riboswitch explored by replica-exchange molecular dynamics simulation. Journal of Theoretical Biology 365: 265–269. pmid:25451761
- View Article
- PubMed/NCBI
- Google Scholar
33. Di Palma F, Colizzi F, Bussi G (2014) Role of Magnesium Ions and Ligand Stacking in the Adenine Riboswitch Folding. Biophysical Journal 106: 284a–284a.
- View Article
- Google Scholar
34. Draper DE (2004) A guide to ions and RNA structure. RNA 10: 335–343. pmid:14970378
- View Article
- PubMed/NCBI
- Google Scholar
35. Misra VK, Draper DE (2001) A thermodynamic framework for Mg2+ binding to RNA. Proceedings of the National Academy of Sciences 98: 12456–12461.
- View Article
- Google Scholar
36. Hayes RL, Noel JK, Mohanty U, Whitford PC, Hennelly SP, et al. (2012) Magnesium Fluctuations Modulate RNA Dynamics in the SAM-I Riboswitch. J Am Chem Soc 134: 12043–12053. pmid:22612276
- View Article
- PubMed/NCBI
- Google Scholar
37. Hayes Ryan L, Noel Jeffrey K, Whitford Paul C, Mohanty U, Sanbonmatsu Karissa Y, et al. (2014) Reduced Model Captures Mg2+-RNA Interaction Free Energy of Riboswitches. Biophysical Journal 106: 1508–1519. pmid:24703312
- View Article
- PubMed/NCBI
- Google Scholar
38. Hayes RL, Noel JK, Mandic A, Whitford PC, Sanbonmatsu KY, et al. (2015) Generalized Manning Condensation Model Captures the RNA Ion Atmosphere. Phys Rev Lett 114.
- View Article
- Google Scholar
39. Quarta G, Kim N, Izzo JA, Schlick T (2009) Analysis of Riboswitch Structure and Function by an Energy Landscape Framework. Journal of Molecular Biology 393: 993–1003. pmid:19733179
- View Article
- PubMed/NCBI
- Google Scholar
40. Chen B, Zuo X, Wang YX, Dayie TK (2011) Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy. Nucleic Acids Research 40: 3117–3130. pmid:22139931
- View Article
- PubMed/NCBI
- Google Scholar
41. Doshi U, Kelley JM, Hamelberg D (2011) Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch. RNA 18: 300–307. pmid:22194311
- View Article
- PubMed/NCBI
- Google Scholar
42. Fuchs RT, Grundy FJ, Henkin TM (2007) S-adenosylmethionine directly inhibits binding of 30S ribosomal subunits to the SMK box translational riboswitch RNA. Proceedings of the National Academy of Sciences 104: 4876–4880.
- View Article
- Google Scholar
43. Lu C, Smith AM, Ding F, Chowdhury A, Henkin TM, et al. (2011) Variable Sequences outside the Sam-Binding Core Critically Influence the Conformational Dynamics of the SAM-III/SMK Box Riboswitch. Journal of Molecular Biology 409: 786–799. pmid:21549712
- View Article
- PubMed/NCBI
- Google Scholar
44. Smith AM, Fuchs RT, Grundy FJ, Henkin TM (2010) The SAM-responsive SMK box is a reversible riboswitch. Molecular Microbiology 78: 1393–1402. pmid:21143313
- View Article
- PubMed/NCBI
- Google Scholar
45. Suddala KC, Wang J, Hou Q, Walter NG (2015) Mg 2+ Shifts Ligand-Mediated Folding of a Riboswitch from Induced-Fit to Conformational Selection. J Am Chem Soc 137: 14075–14083. pmid:26471732
- View Article
- PubMed/NCBI
- Google Scholar
46. Hennelly SP, Novikova IV, Sanbonmatsu KY (2012) The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch. Nucleic Acids Research 41: 1922–1935. pmid:23258703
- View Article
- PubMed/NCBI
- Google Scholar
47. Chen B, LeBlanc R, Dayie TK (2016) SAM-II Riboswitch Samples at least Two Conformations in Solution in the Absence of Ligand: Implications for Recognition. Angewandte Chemie International Edition 55: 2724–2727. pmid:26800479
- View Article
- PubMed/NCBI
- Google Scholar
48. Pelikan M, Hura GL, Hammel M (2009) Structure and flexibility within proteins as identified through small angle X-ray scattering. General Physiology and Biophysics 28: 174–189.
- View Article
- Google Scholar
49. Schneidman-Duhovny D, Hammel M, Sali A (2010) FoXS: a web server for rapid computation and fitting of SAXS profiles. Nucleic Acids Research 38: W540–W544. pmid:20507903
- View Article
- PubMed/NCBI
- Google Scholar
50. Bergmann JH, Spector DL (2014) Long non-coding RNAs: modulators of nuclear structure and function. Current Opinion in Cell Biology 26: 10–18. pmid:24529241
- View Article
- PubMed/NCBI
- Google Scholar
51. Brown JA, Valenstein ML, Yario TA, Tycowski KT, Steitz JA (2012) Formation of triple-helical structures by the 3'-end sequences of MALAT1 and MENbeta noncoding RNAs. Proc Natl Acad Sci U S A 109: 19202–19207. pmid:23129630
- View Article
- PubMed/NCBI
- Google Scholar
52. Conrad NK (2014) The emerging role of triple helices in RNA biology. Wiley Interdiscip Rev RNA 5: 15–29. pmid:24115594
- View Article
- PubMed/NCBI
- Google Scholar
53. Wilusz JE, JnBaptiste CK, Lu LY, Kuhn CD, Joshua-Tor L, et al. (2012) A triple helix stabilizes the 3 ' ends of long noncoding RNAs that lack poly(A) tails. Genes & Development 26: 2392–2407.
- View Article
- Google Scholar
54. Lescoute A, Westhof E (2006) The interaction networks of structured RNAs. Nucleic Acids Research 34: 6587–6604. pmid:17135184
- View Article
- PubMed/NCBI
- Google Scholar
55. Villa A, Wohnert J, Stock G (2009) Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch. Nucleic Acids Research 37: 4774–4786. pmid:19515936
- View Article
- PubMed/NCBI
- Google Scholar
56. Lipfert J, Doniach S, Das R, Herschlag D (2014) Understanding Nucleic Acid-Ion Interactions. Annual Review of Biochemistry, Vol 83 83: 813–841.
- View Article
- Google Scholar
57. Ahmadi MS, Fattahi A (2011) On the binding of Mg2+, Ca2+, Zn2+ and Cu+ metal cations to 2 '-deoxyguanosine: Changes on sugar puckering and strength of the N-glycosidic bond. Scientia Iranica 18: 1343–1352.
- View Article
- Google Scholar
58. Espah Borujeni A, Mishler DM, Wang J, Huso W, Salis HM (2016) Automated physics-based design of synthetic riboswitches from diverse RNA aptamers. Nucleic Acids Research 44: 1–13. pmid:26621913
- View Article
- PubMed/NCBI
- Google Scholar
59. Jana B, Morcos F, Onuchic JN (2014) From structure to function: the convergence of structure based models and co-evolutionary information. Physical Chemistry Chemical Physics 16: 6496–6507. pmid:24603809
- View Article
- PubMed/NCBI
- Google Scholar
60. Noel JK, Whitford PC, Onuchic JN (2012) The shadow map: a general contact definition for capturing the dynamics of biomolecular folding and function. Journal of Physical Chemistry B 116: 8692–8702.
- View Article
- Google Scholar
61. Noel JK, Whitford PC, Sanbonmatsu KY, Onuchic JN (2010) SMOG@ctbp: simplified deployment of structure-based models in GROMACS. Nucleic Acids Research 38: W657–661. pmid:20525782
- View Article
- PubMed/NCBI
- Google Scholar
62. Onuchic JN, Luthey-Schulten Z, Wolynes PG (1997) Theory of protein folding: the energy landscape perspective. Annu Rev Phys Chem 48: 545–600. pmid:9348663
- View Article
- PubMed/NCBI
- Google Scholar
63. Wolynes PG (2015) Evolution, energy landscapes and the paradoxes of protein folding. Biochimie 119: 218–230. pmid:25530262
- View Article
- PubMed/NCBI
- Google Scholar
64. Manning GS (1978) The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides. Q Rev Biophys 11: 179–246. pmid:353876
- View Article
- PubMed/NCBI
- Google Scholar
65. Torrie GM, Valleau JP (1974) Monte-Carlo Free-Energy Estimates Using Non-Boltzmann Sampling—Application to Subcritical Lennard-Jones Fluid. Chemical Physics Letters 28: 578–581.
- View Article
- Google Scholar
66. Kumar S, Bouzida D, Swendsen RH, Kollman PA, Rosenberg JM (1992) The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules .1. The Method. Journal of Computational Chemistry 13: 1011–1021.
- View Article
- Google Scholar
67. Garcia-Mayoral MF, Trevino MA, Perez-Pinar T, Caballero ML, Knaute T, et al. (2014) Relationships between IgE/IgG4 epitopes, structure and function in Anisakis simplex Ani s 5, a member of the SXP/RAL-2 protein family. PLoS Negl Trop Dis 8: e2735. pmid:24603892
- View Article
- PubMed/NCBI
- Google Scholar
68. Kobitski AY, Nierth A, Helm M, Jaschke A, Nienhaus GU (2007) Mg2+-dependent folding of a Diels-Alderase ribozyme probed by single-molecule FRET analysis. Nucleic Acids Research 35: 2047–2059. pmid:17344321
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Huarte M (2015) The emerging role of lncRNAs in cancer. Nat Med 21: 1253–1261. pmid:26540387
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Montange RK, Batey RT (2008) Riboswitches: Emerging Themes in RNA Structure and Function. Annu Rev Biophys 37: 117–133. pmid:18573075
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Roth A, Breaker RR (2009) The Structural and Functional Diversity of Metabolite-Binding Riboswitches. Annu Rev Biochem 78: 305–334. pmid:19298181
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Lin JC, Yoon J, Hyeon C, Thirumalai D (2015) Using simulations and kinetic network models to reveal the dynamics and functions of riboswitches. Methods Enzymol 553: 235–258. pmid:25726468
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Serganov A, Patel DJ (2007) Ribozymes, riboswitches and beyond: regulation of gene expression without proteins. Nature Reviews Genetics 8: 776–790. pmid:17846637
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Winkler WC, Nahvi A, Sudarsan N, Barrick JE, Breaker RR (2003) An mRNA structure that controls gene expression by binding S-adenosylmethionine. Nature Structural Biology 10: 701–707. pmid:12910260
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Corbino KA, Barrick JE, Lim J, Welz R, Tucker BJ, et al. (2005) Genome Biology 6: R70. pmid:16086852
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Yamauchi T, Miyoshi D, Kubodera T, Nishimura A, Nakai S, et al. (2005) Roles of Mg2+ in TPP-dependent riboswitch. FEBS Lett 579: 2583–2588. pmid:15862294
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Lipfert J, Das R, Chu VB, Kudaravalli M, Boyd N, et al. (2007) Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae. J Mol Biol 365: 1393–1406. pmid:17118400
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Bisaria N, Greenfeld M, Limouse C, Pavlichin DS, Mabuchi H, et al. (2016) Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway. Proc Natl Acad Sci U S A 113: E4956–4965. pmid:27493222
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Bisaria N, Herschlag D (2015) Probing the kinetic and thermodynamic consequences of the tetraloop/tetraloop receptor monovalent ion-binding site in P4-P6 RNA by smFRET. Biochem Soc Trans 43: 172–178. pmid:25849913
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Takamoto K, Das R, He Q, Doniach S, Brenowitz M, et al. (2004) Principles of RNA compaction: insights from the equilibrium folding pathway of the P4-P6 RNA domain in monovalent cations. J Mol Biol 343: 1195–1206. pmid:15491606
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Xue Y, Gracia B, Herschlag D, Russell R, Al-Hashimi HM (2016) Visualizing the formation of an RNA folding intermediate through a fast highly modular secondary structure switch. Nat Commun 7: ncomms11768. pmid:27292179
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Bocobza SE, Aharoni A (2014) Small molecules that interact with RNA: riboswitch-based gene control and its involvement in metabolic regulation in plants and algae. Plant J 79: 693–703. pmid:24773387
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Russell R (2015) Unwinding the mechanisms of a DEAD-box RNA helicase in cancer. J Mol Biol 427: 1797–1800. pmid:25836982
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Batey RT, Gilbert SD, Montange RK (2004) Structure of a natural guanine-responsive riboswitch complexed with the metabolite hypoxanthine. Nature 432: 411–415. pmid:15549109
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Edwards AL, Reyes FE, Heroux A, Batey RT (2010) Structural basis for recognition of S-adenosylhomocysteine by riboswitches. RNA 16: 2144–2155. pmid:20864509
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref18] 18. Garst AD, Heroux A, Rambo RP, Batey RT (2008) Crystal structure of the lysine riboswitch regulatory mRNA element. J Biol Chem 283: 22347–22351. pmid:18593706
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref19] 19. Gilbert SD, Montange RK, Stoddard CD, Batey RT (2006) Structural Studies of the Purine and SAM Binding Riboswitches. Cold Spring Harbor Symposia on Quantitative Biology 71: 259–268. pmid:17381305
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref20] 20. Gilbert SD, Rambo RP, Van Tyne D, Batey RT (2008) Structure of the SAM-II riboswitch bound to S-adenosylmethionine. Nat Struct Mol Biol 15: 177–182. pmid:18204466
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref21] 21. Haller A, Altman RB, Souliere MF, Blanchard SC, Micura R (2013) Folding and ligand recognition of the TPP riboswitch aptamer at single-molecule resolution. Proceedings of the National Academy of Sciences 110: 4188–4193.
View Article
Google Scholar

[82] View Article

[83] Google Scholar

[ref22] 22. Haller A, Rieder U, Aigner M, Blanchard SC, Micura R (2011) Conformational capture of the SAM-II riboswitch. Nature Chemical Biology 7: 393–400. pmid:21532598
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref23] 23. Hennelly SP, Sanbonmatsu KY (2010) Tertiary contacts control switching of the SAM-I riboswitch. Nucleic Acids Research 39: 2416–2431. pmid:21097777
View Article
PubMed/NCBI
Google Scholar

[89] View Article

[90] PubMed/NCBI

[91] Google Scholar

[ref24] 24. Heppell B, Mulhbacher J, Penedo JC, Lafontaine DA (2009) Application of Fluorescent Measurements for Characterization of Riboswitch–Ligand Interactions. Methods in Molecular Biology: Springer Science + Business Media. pp. 25–37.

[ref25] 25. Lemay J-F, Penedo JC, Tremblay R, Lilley DMJ, Lafontaine Daniel A (2006) Folding of the Adenine Riboswitch. Chemistry & Biology 13: 857–868.
View Article
Google Scholar

[94] View Article

[95] Google Scholar

[ref26] 26. Montange RK, Batey RT (2006) Structure of the S-adenosylmethionine riboswitch regulatory mRNA element. Nature 441: 1172–1175. pmid:16810258
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref27] 27. Stoddard CD, Montange RK, Hennelly SP, Rambo RP, Sanbonmatsu KY, et al. (2010) Free State Conformational Sampling of the SAM-I Riboswitch Aptamer Domain. Structure 18: 787–797. pmid:20637415
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref28] 28. Wacker A, Buck J, Mathieu D, Richter C, Wohnert J, et al. (2011) Structure and dynamics of the deoxyguanosine-sensing riboswitch studied by NMR-spectroscopy. Nucleic Acids Research 39: 6802–6812. pmid:21576236
View Article
PubMed/NCBI
Google Scholar

[105] View Article

[106] PubMed/NCBI

[107] Google Scholar

[ref29] 29. Cho SS, Pincus DL, Thirumalai D (2009) Assembly mechanisms of RNA pseudoknots are determined by the stabilities of constituent secondary structures. Proceedings of the National Academy of Sciences 106: 17349–17354.
View Article
Google Scholar

[109] View Article

[110] Google Scholar

[ref30] 30. Hyeon C, Thirumalai D (2005) Mechanical unfolding of RNA hairpins. Proceedings of the National Academy of Sciences 102: 6789–6794.
View Article
Google Scholar

[112] View Article

[113] Google Scholar

[ref31] 31. Whitford PC, Schug A, Saunders J, Hennelly SP, Onuchic JN, et al. (2009) Nonlocal Helix Formation Is Key to Understanding S-Adenosylmethionine-1 Riboswitch Function. Biophysical Journal 96: L7–L9. pmid:19167285
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref32] 32. Xue X, Yongjun W, Zhihong L (2015) Folding of SAM-II riboswitch explored by replica-exchange molecular dynamics simulation. Journal of Theoretical Biology 365: 265–269. pmid:25451761
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref33] 33. Di Palma F, Colizzi F, Bussi G (2014) Role of Magnesium Ions and Ligand Stacking in the Adenine Riboswitch Folding. Biophysical Journal 106: 284a–284a.
View Article
Google Scholar

[123] View Article

[124] Google Scholar

[ref34] 34. Draper DE (2004) A guide to ions and RNA structure. RNA 10: 335–343. pmid:14970378
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref35] 35. Misra VK, Draper DE (2001) A thermodynamic framework for Mg2+ binding to RNA. Proceedings of the National Academy of Sciences 98: 12456–12461.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref36] 36. Hayes RL, Noel JK, Mohanty U, Whitford PC, Hennelly SP, et al. (2012) Magnesium Fluctuations Modulate RNA Dynamics in the SAM-I Riboswitch. J Am Chem Soc 134: 12043–12053. pmid:22612276
View Article
PubMed/NCBI
Google Scholar

[133] View Article

[134] PubMed/NCBI

[135] Google Scholar

[ref37] 37. Hayes Ryan L, Noel Jeffrey K, Whitford Paul C, Mohanty U, Sanbonmatsu Karissa Y, et al. (2014) Reduced Model Captures Mg2+-RNA Interaction Free Energy of Riboswitches. Biophysical Journal 106: 1508–1519. pmid:24703312
View Article
PubMed/NCBI
Google Scholar

[137] View Article

[138] PubMed/NCBI

[139] Google Scholar

[ref38] 38. Hayes RL, Noel JK, Mandic A, Whitford PC, Sanbonmatsu KY, et al. (2015) Generalized Manning Condensation Model Captures the RNA Ion Atmosphere. Phys Rev Lett 114.
View Article
Google Scholar

[141] View Article

[142] Google Scholar

[ref39] 39. Quarta G, Kim N, Izzo JA, Schlick T (2009) Analysis of Riboswitch Structure and Function by an Energy Landscape Framework. Journal of Molecular Biology 393: 993–1003. pmid:19733179
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref40] 40. Chen B, Zuo X, Wang YX, Dayie TK (2011) Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy. Nucleic Acids Research 40: 3117–3130. pmid:22139931
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref41] 41. Doshi U, Kelley JM, Hamelberg D (2011) Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch. RNA 18: 300–307. pmid:22194311
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref42] 42. Fuchs RT, Grundy FJ, Henkin TM (2007) S-adenosylmethionine directly inhibits binding of 30S ribosomal subunits to the SMK box translational riboswitch RNA. Proceedings of the National Academy of Sciences 104: 4876–4880.
View Article
Google Scholar

[156] View Article

[157] Google Scholar

[ref43] 43. Lu C, Smith AM, Ding F, Chowdhury A, Henkin TM, et al. (2011) Variable Sequences outside the Sam-Binding Core Critically Influence the Conformational Dynamics of the SAM-III/SMK Box Riboswitch. Journal of Molecular Biology 409: 786–799. pmid:21549712
View Article
PubMed/NCBI
Google Scholar

[159] View Article

[160] PubMed/NCBI

[161] Google Scholar

[ref44] 44. Smith AM, Fuchs RT, Grundy FJ, Henkin TM (2010) The SAM-responsive SMK box is a reversible riboswitch. Molecular Microbiology 78: 1393–1402. pmid:21143313
View Article
PubMed/NCBI
Google Scholar

[163] View Article

[164] PubMed/NCBI

[165] Google Scholar

[ref45] 45. Suddala KC, Wang J, Hou Q, Walter NG (2015) Mg 2+ Shifts Ligand-Mediated Folding of a Riboswitch from Induced-Fit to Conformational Selection. J Am Chem Soc 137: 14075–14083. pmid:26471732
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref46] 46. Hennelly SP, Novikova IV, Sanbonmatsu KY (2012) The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch. Nucleic Acids Research 41: 1922–1935. pmid:23258703
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref47] 47. Chen B, LeBlanc R, Dayie TK (2016) SAM-II Riboswitch Samples at least Two Conformations in Solution in the Absence of Ligand: Implications for Recognition. Angewandte Chemie International Edition 55: 2724–2727. pmid:26800479
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref48] 48. Pelikan M, Hura GL, Hammel M (2009) Structure and flexibility within proteins as identified through small angle X-ray scattering. General Physiology and Biophysics 28: 174–189.
View Article
Google Scholar

[179] View Article

[180] Google Scholar

[ref49] 49. Schneidman-Duhovny D, Hammel M, Sali A (2010) FoXS: a web server for rapid computation and fitting of SAXS profiles. Nucleic Acids Research 38: W540–W544. pmid:20507903
View Article
PubMed/NCBI
Google Scholar

[182] View Article

[183] PubMed/NCBI

[184] Google Scholar

[ref50] 50. Bergmann JH, Spector DL (2014) Long non-coding RNAs: modulators of nuclear structure and function. Current Opinion in Cell Biology 26: 10–18. pmid:24529241
View Article
PubMed/NCBI
Google Scholar

[186] View Article

[187] PubMed/NCBI

[188] Google Scholar

[ref51] 51. Brown JA, Valenstein ML, Yario TA, Tycowski KT, Steitz JA (2012) Formation of triple-helical structures by the 3'-end sequences of MALAT1 and MENbeta noncoding RNAs. Proc Natl Acad Sci U S A 109: 19202–19207. pmid:23129630
View Article
PubMed/NCBI
Google Scholar

[190] View Article

[191] PubMed/NCBI

[192] Google Scholar

[ref52] 52. Conrad NK (2014) The emerging role of triple helices in RNA biology. Wiley Interdiscip Rev RNA 5: 15–29. pmid:24115594
View Article
PubMed/NCBI
Google Scholar

[194] View Article

[195] PubMed/NCBI

[196] Google Scholar

[ref53] 53. Wilusz JE, JnBaptiste CK, Lu LY, Kuhn CD, Joshua-Tor L, et al. (2012) A triple helix stabilizes the 3 ' ends of long noncoding RNAs that lack poly(A) tails. Genes & Development 26: 2392–2407.
View Article
Google Scholar

[198] View Article

[199] Google Scholar

[ref54] 54. Lescoute A, Westhof E (2006) The interaction networks of structured RNAs. Nucleic Acids Research 34: 6587–6604. pmid:17135184
View Article
PubMed/NCBI
Google Scholar

[201] View Article

[202] PubMed/NCBI

[203] Google Scholar

[ref55] 55. Villa A, Wohnert J, Stock G (2009) Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch. Nucleic Acids Research 37: 4774–4786. pmid:19515936
View Article
PubMed/NCBI
Google Scholar

[205] View Article

[206] PubMed/NCBI

[207] Google Scholar

[ref56] 56. Lipfert J, Doniach S, Das R, Herschlag D (2014) Understanding Nucleic Acid-Ion Interactions. Annual Review of Biochemistry, Vol 83 83: 813–841.
View Article
Google Scholar

[209] View Article

[210] Google Scholar

[ref57] 57. Ahmadi MS, Fattahi A (2011) On the binding of Mg2+, Ca2+, Zn2+ and Cu+ metal cations to 2 '-deoxyguanosine: Changes on sugar puckering and strength of the N-glycosidic bond. Scientia Iranica 18: 1343–1352.
View Article
Google Scholar

[212] View Article

[213] Google Scholar

[ref58] 58. Espah Borujeni A, Mishler DM, Wang J, Huso W, Salis HM (2016) Automated physics-based design of synthetic riboswitches from diverse RNA aptamers. Nucleic Acids Research 44: 1–13. pmid:26621913
View Article
PubMed/NCBI
Google Scholar

[215] View Article

[216] PubMed/NCBI

[217] Google Scholar

[ref59] 59. Jana B, Morcos F, Onuchic JN (2014) From structure to function: the convergence of structure based models and co-evolutionary information. Physical Chemistry Chemical Physics 16: 6496–6507. pmid:24603809
View Article
PubMed/NCBI
Google Scholar

[219] View Article

[220] PubMed/NCBI

[221] Google Scholar

[ref60] 60. Noel JK, Whitford PC, Onuchic JN (2012) The shadow map: a general contact definition for capturing the dynamics of biomolecular folding and function. Journal of Physical Chemistry B 116: 8692–8702.
View Article
Google Scholar

[223] View Article

[224] Google Scholar

[ref61] 61. Noel JK, Whitford PC, Sanbonmatsu KY, Onuchic JN (2010) SMOG@ctbp: simplified deployment of structure-based models in GROMACS. Nucleic Acids Research 38: W657–661. pmid:20525782
View Article
PubMed/NCBI
Google Scholar

[226] View Article

[227] PubMed/NCBI

[228] Google Scholar

[ref62] 62. Onuchic JN, Luthey-Schulten Z, Wolynes PG (1997) Theory of protein folding: the energy landscape perspective. Annu Rev Phys Chem 48: 545–600. pmid:9348663
View Article
PubMed/NCBI
Google Scholar

[230] View Article

[231] PubMed/NCBI

[232] Google Scholar

[ref63] 63. Wolynes PG (2015) Evolution, energy landscapes and the paradoxes of protein folding. Biochimie 119: 218–230. pmid:25530262
View Article
PubMed/NCBI
Google Scholar

[234] View Article

[235] PubMed/NCBI

[236] Google Scholar

[ref64] 64. Manning GS (1978) The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides. Q Rev Biophys 11: 179–246. pmid:353876
View Article
PubMed/NCBI
Google Scholar

[238] View Article

[239] PubMed/NCBI

[240] Google Scholar

[ref65] 65. Torrie GM, Valleau JP (1974) Monte-Carlo Free-Energy Estimates Using Non-Boltzmann Sampling—Application to Subcritical Lennard-Jones Fluid. Chemical Physics Letters 28: 578–581.
View Article
Google Scholar

[242] View Article

[243] Google Scholar

[ref66] 66. Kumar S, Bouzida D, Swendsen RH, Kollman PA, Rosenberg JM (1992) The Weighted Histogram Analysis Method for Free-Energy Calculations on Biomolecules .1. The Method. Journal of Computational Chemistry 13: 1011–1021.
View Article
Google Scholar

[245] View Article

[246] Google Scholar

[ref67] 67. Garcia-Mayoral MF, Trevino MA, Perez-Pinar T, Caballero ML, Knaute T, et al. (2014) Relationships between IgE/IgG4 epitopes, structure and function in Anisakis simplex Ani s 5, a member of the SXP/RAL-2 protein family. PLoS Negl Trop Dis 8: e2735. pmid:24603892
View Article
PubMed/NCBI
Google Scholar

[248] View Article

[249] PubMed/NCBI

[250] Google Scholar

[ref68] 68. Kobitski AY, Nierth A, Helm M, Jaschke A, Nienhaus GU (2007) Mg2+-dependent folding of a Diels-Alderase ribozyme probed by single-molecule FRET analysis. Nucleic Acids Research 35: 2047–2059. pmid:17344321
View Article
PubMed/NCBI
Google Scholar

[252] View Article

[253] PubMed/NCBI

[254] Google Scholar

Figures

Abstract

Author summary

Introduction

Results

Detailed agreement between 13C-CEST, smFRET, SAXS, SEC elution profile and generalized Manning model corroborates pre-organization by Mg2+

Mg2+ remodels the free energy landscape, favoring pre-organized states

Mg2+-induced phosphate interactions facilitate pre-organization

Mg2+-induced triplex interactions between helix P2b and L1 loop help pre-organize the SAM-II riboswitch

Discussion

Methods

RNA electrostatic model

Free energy simulations

13C-CEST experiments

Analysis of relaxation parameters

FRET

SAXS analysis

Supporting information

S1 Text. Supporting figures, table and text.

Author Contributions

References

Detailed agreement between ¹³C-CEST, smFRET, SAXS, SEC elution profile and generalized Manning model corroborates pre-organization by Mg²⁺

Mg²⁺ remodels the free energy landscape, favoring pre-organized states

Mg²⁺-induced phosphate interactions facilitate pre-organization

Mg²⁺-induced triplex interactions between helix P2b and L1 loop help pre-organize the SAM-II riboswitch

¹³C-CEST experiments