Ligand response of guanidine-IV riboswitch at single-molecule level
eLife assessment
This study presents valuable findings on the ligand- and ion-dependent structural dynamics of a transcriptional riboswitch. The single-molecule data presented are solid and prompts intriguing hypotheses and models, which will undoubtedly stimulate future structural analyses. These findings are of considerable interest to biochemists and biophysicists engaged in the study of RNA structure and riboswitch mechanisms.
https://doi.org/10.7554/eLife.94706.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Solid: Methods, data and analyses broadly support the claims with only minor weaknesses
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.
Introduction
Riboswitches, located at the 5′-untranslated region of mRNA, are capable of regulating gene expression through structural changes after binding to their specific ligands (Jones and Ferré-D’Amaré, 2017; Mandal and Breaker, 2004; Serganov and Nudler, 2013; Sherlock and Breaker, 2020). A typical riboswitch consists of an aptamer domain and an expression platform. The aptamer domain binds to a specific ligand, leading to structural changes and regulation of gene expression (Garst et al., 2011; Kavita and Breaker, 2023; Winkler and Breaker, 2003). To date, more than 55 classes of riboswitches have been identified that can sense a diverse range of ligands, including small metabolites, anions, cations, amino acids, and nucleotides (Breaker, 2022; McCown et al., 2017). Furthermore, riboswitches responding to the same ligands have been found to exhibit distinct consequences and structures. For instance, there are four types of guanidine riboswitches, namely guanidine-I, guanidine-II, guanidine-III, and guanidine-IV riboswitch, which specifically bind to guanidine. Notably, the guanidine-IV riboswitch differs significantly from other guanidine riboswitches (Battaglia and Ke, 2018; Huang et al., 2017a; Huang et al., 2017b; Lenkeit et al., 2020; Nelson et al., 2017; Salvail et al., 2020). Specifically, the guanidine-IV riboswitch is able to enhance the transcription of mepA by forming an anti-terminator structure upon guanidine binding. This structure plays a crucial role in exporting guanidine and reducing its toxicity in cells, making the guanidine-IV riboswitch a potential target for antibiotics (Kermani et al., 2018). Although the atomic-resolution structures of the guanidine-IV riboswitch have not been solved, it is proposed that a kissing loop (KL) forms in the riboswitch in the presence of guanidine. This KL facilitates the dissociation of the terminator and promotes transcriptional read-through (Lenkeit et al., 2020; Salvail et al., 2020). KL or pseudoknot formation upon ligand binding is a common occurrence in various riboswitches, including adenine, guanine, ZTP, pre-Q1, fluoride, and c-di-AMP riboswitches (Jones and Ferré-D’Amaré, 2017; Lenkeit et al., 2020; Salvail et al., 2020). However, the mechanism by which KL affects gene expression in a transcriptional riboswitch, particularly during transcription, remains elusive.
The folding of RNA during transcription plays a critical role in its function, especially in riboswitches involved in gene regulation. Previous studies have investigated the co-transcriptional folding of various riboswitches using different techniques. For instance, Frieda and Block, 2012 employed the optical-trapping assay to examine the co-transcriptional folding of an adenine riboswitch at the single-molecule resolution. Helmling et al., 2017 and Binas et al., 2020, on the other hand, used NMR to map the transcriptional intermediates of the 2′-deoxyguanosine (2′-dG) and ZMP riboswitches at the single-nucleotide resolution. Landgraf et al., 2022 utilized NMR and computational modeling to identify the ligand-sensing transcriptional window for the c-di-GMP and c-GAMP riboswitches. Hua et al. investigated the co-transcriptional folding of the twister ribozyme and the ZTP riboswitch by mimicking the folding process using helicases to dissociate an RNA-DNA hybrid (Hua et al., 2020; Hua et al., 2018). Lou and Woodson, 2024, Uhm et al., 2018, Widom et al., 2018, and Yadav et al., 2022 employed the single-molecule fluorescence resonance energy transfer (smFRET) strategy to study the co-transcriptional folding of the glmS ribozyme, TPP, pre-Q1, and fluoride riboswitches by assembling transcriptional complexes. Xue et al., 2023 monitored the folding of a growing SAM-VI riboswitch at the single-molecule resolution by FRET. Despite these previous investigations, there is lack of information regarding the conformational changes and ligand-sensing behavior of the guanidine-IV riboswitch during transcription. In this study, we used position-specific labeling of RNA (PLOR) to introduce a FRET pair, Cy3 and Cy5, at sites 78 and 35, respectively (Liu et al., 2018; Liu et al., 2015). Subsequently, we examined the ligand and Mg2+ responses of the aptamer domain, terminator/anti-terminator, and full-length guanidine-IV riboswitch using smFRET. Our findings revealed that both guanidine and Mg2+ facilitated the formation of a KL in the aptamer domain and significantly increased the KL-formed conformation in that domain. Interestingly, unlike the aptamer domain, the addition of guanidine and Mg2+ only had a minimal effect on the formation of the KL-formed anti-terminator structure in the terminator/anti-terminator or full-length RNA. Additionally, we mimicked the conformational changes and ligand-sensing behavior from the aptamer domain to the terminator sequence of the guanidine-IV riboswitch during transcription using smFRET and PLOR. The formation of the KL increased as the riboswitch was transcribed up to 88 nucleotides (corresponding to the aptamer), but decreased as the riboswitch was transcribed from 88 to 105 nucleotides (corresponding to the terminator). Based on our data, it became evident that the transcriptional fate did not determined by the terminator sequence, but rather by the aptamer region, which exhibited distinct ligand-sensing abilities.
Results
The construction and characteristics of the guanidine-IV riboswitch
The Clostridium botulinum guanidine-IV riboswitch was investigated in this study. The secondary structures of the guanidine-IV riboswitch are depicted in Figure 1A and B. In the absence of the guanidinium cation (Gua+), the terminator consists of the aptamer domain (in black) and poly-U tail (in red). This terminator is proposed to terminate transcription by dissociating RNA polymerase from the transcription system (Figure 1A). However, in the presence of Gua+, the formation of a KL (in green, Figure 1B) hinders the formation of the complete terminator stem. This leads to the formation of the anti-terminator conformation, and transcription of the full-length riboswitch comprises the extended nucleotides (in blue), in addition to the terminator/anti-terminator region (in red, Figure 1B and C). Consequently, the anti-terminator conformation activates the expression of the downstream genes (Figure 1C and D). In our study, we noted that the transcriptional read-through of the guanidine-IV riboswitch during the single-round PLOR reaction was sensitive to Gua+, exhibiting an apparent EC50 value of 68.7±7.3 μM (Figure 1D; Chien et al., 2023). This is comparable to the reported KD of approximately 64 μM that was determined by in-line probing (Salvail et al., 2020). Higher read-through efficiencies were detected at higher Gua+ concentrations, consistent with previous findings (Lenkeit et al., 2020). In order to probe the structural changes of different domains in the guanidine-IV riboswitch under various conditions, we applied PLOR to introduce Cy3 and Cy5 at sites 78 and 35 in the aptamer (riboG-apt), terminator/anti-terminator (riboG-term), and full-length guanidine-IV riboswitch (riboG). Sites 78 and 35 are close to the KL, making them suitable for monitoring the KL formation at different Mg2+ and Gua+ concentrations.
smFRET revealed the significant effects of Mg2+ and Gua+ on riboG-apt
smFRET is a powerful method that can detect multiple conformations by measuring the energy transfer efficiency (EFRET) of a pair of fluorophores at the single molecule level. This technique has been extensively employed to investigate the conformational dynamics of RNA at the single-molecule level (Manz et al., 2017; Uhm et al., 2018). When the riboG-apt labeled with dual fluorophore-labeled (Figure 2A, Figure 2—figure supplement 1) was examined, a single peak (EFRET ~0.2) was observed in the presence of 0 mM Mg2+ and 0 mM Gua+ (Figure 2B). This peak is believed to correspond to the unfolded structure of riboG-apt, where the KL is not formed and thus, the labeled sites are distant from each other. Furthermore, in the presence of 2.0 mM Mg2+, three peaks (EFRET ~0.2, 0.5, and 0.8) were observed, as revealed by hidden Markov modeling (Figure 2C). This suggests the coexistence of at least three structures when 2.0 mM Mg2+ is present, with the peak at an EFRET of approximately 0.8 likely representing the KL formation. The peak characterized by an EFRET ~0.5 emerged at varying collection rates and was indicative of a more compact conformation compared to the unfolded structure, which we referred to as the pre-folded structure (Figure 2D, Figure 2—figure supplements 2 and 3). Notably, at 2.0 mM Mg2+, these three states exhibited high dynamics, with transitions occurring particularly between the pre-folded and the unfolded states (Figure 2B,C, Figure 2—figure supplements 2 and 3). Our results demonstrate that the presence of Mg2+ induces a conformational change in the unfolded riboG-apt, leading to the adoption of more compact pre-folded and folded conformations. The proportion of the folded conformation significantly increased to approximately 90% at 20.0 mM Mg2+ (Figure 2E, Figure 2—figure supplement 4G). At high concentrations of Mg2+, riboG-apt remained flexible and dynamic in its transition between states (Figure 2—figure supplement 4). However, at a higher concentration of 50.0 mM Mg2+, the proportion of the pre-folded and unfolded conformations were more prevalent at 50.0 mM Mg2+ than at 20.0 mM Mg2+. This suggests that an excess of Mg2+ may promote the pre-folded and even unfolded conformations.
Gua+ also had the ability to promote the folded conformation (EFRET ~0.8) of riboG-apt. However, Gua+ exhibited lower efficacy compared to Mg2+, as evidenced by the proportion of the folded conformation being approximately 25% in the presence of 20.0 mM Gua+. This value is significantly lower than the approximately 90% observed in the presence of 20.0 mM Mg2+ (Figure 2E, Figure 2—figure supplements 4G and 5F). Our observations also revealed that riboG-apt exhibited pronounced dynamics, frequently transiting among the three states at 10.0–100.0 mM Gua+ (Figure 2—figure supplement 5E–H). The presence of Mg2+ was critical in facilitating the induction of the folded state by Gua+. Notably, riboG-apt showed increased susceptibility to Gua+ in the presence of Mg2+ (Figure 2F, Figure 2—figure supplement 6). Upon the addition of 0.5 mM Gua+ and 0.5 mM Mg2+, the KL-folded conformation predominated in riboG-apt. The proportion of the KL-folded conformation increased to over 40%, which was comparable to the effect observed with a 20-fold increase in 0.5 mM Gua+ in the absence of Mg2+ (Figure 2—figure supplements 5B and 6E). Interestingly, at 0.5 mM Mg2+, Gua+ bound to riboG-apt with a Kd of 286.0±18.1 μM (Figure 2—figure supplement 6H), higher than the Kd of 64 μM that was determined by in-line probing in the presence of 20 mM Mg2+ (Salvail et al., 2020). This difference is likely due to the low concentration of Mg2+, which has a weak ability to help guanidine induced the KL conformation. In addition, we observed that riboG-apt became more stable as the concentration of Gua+ increased, particularly in the presence of high Mg2+ (Figure 2—figure supplement 7). Surprisingly, the proportion of the folded conformation hardly increased when Mg2+ exceeded 2.0 mM in the presence of 1.0 mM Gua+. Taking into account that Gua+ was present in a positively charged state in the buffer employed for our FRET studies, we conducted additional experiments to investigate the impact of commonly used cations, including Na+, K+, and urea (a structural analog of guanidine). Interestingly, we observed minimal changes in the folded conformation of riboG-apt upon the addition of Na+, K+, or urea. These findings suggest that the effects of Gua+ on the folded conformation are distinct from those of the tested cations. Furthermore, our results indicate that Gua+ specifically binds to riboG-apt and induces structural changes as a ligand, rather than simply acting as a general cation or chaotrope (Figure 2G).
The highly conserved nucleotides surrounding the KL are crucial for its formation (Lenkeit et al., 2020). To test our hypothesis that the state with EFRET ~0.8 corresponds to the conformation with the KL, we preformed smFRET analysis on several mutations at these crucial nucleotides (Figure 2—figure supplements 8–10). Consistent with our expectations, the peaks with EFRET ~0.8 was significantly diminished in the riboG-G71C mutant, which features a single nucleotide mutation at site 71 (with 97% nucleotide conservation) in the KL (Figure 2—figure supplement 8A and B). It is worth noting that the C30G and G71C mutant, which were initially expected to restore a base pair in the KL, did not successfully bring about the anticipated peak of EFRET ~0.8 (Figure 2—figure supplement 8C and D). On the other hand, the riboG-U72C mutant exhibited a lower proportion at the state with EFRET ~0.8 than riboG-apt. However, the A29G and U72C mutations restored a base pair in the KL, as well as the formation of the KL (Figure 2—figure supplement 9). Furthermore, our investigation revealed that the G77C mutant, involving a single nucleotide mutation at a highly conversed site, 77 (with 97% nucleotide conservation), also hindered the formation of the KL (Figure 2—figure supplement 10). This finding aligns with previous research (Lenkeit et al., 2020) and the predicted second structure of G77C mutation by Mfold (Zuker, 2003).
smFRET analysis revealed the slight effects of Mg2+ and Gua+ on riboG-term and full-length riboG
Unlike riboG-apt, riboG-term labeled with dual fluorophores exhibited the coexistence of three conformations, characterized by EFRET values of approximately 0.2, 0.5, and 0.8 (Figure 3A, Figure 3—figure supplements 1–5). This coexistence was observed even in the absence of Mg2+ and Gua+ (Figure 3B, Figure 3—figure supplement 2A and 2B). This indicates that monovalent ions in the buffer can facilitate the formation of stable guanidine-IV riboswitch. The sensitivity of RiboG-term to Gua+ was significantly lower compared to riboG-apt. This is evidenced by the fact that, in the presence of 2.0 mM Mg2+, the addition of 1.0 mM Gua+ only resulted in an increase of approximately 4% the folded conformation for riboG-term, while it was approximately 60% for riboG-apt (Figures 2B and 3B). With the increase of either Mg2+ from 0 to 50.0 mM or Gua+ from 0 to 100.0 mM, less than 12% of the molecules transited dynamically among the three states of riboG-term (Figure 3—figure supplements 2 and 3). This finding suggests that the unfolded structure of riboG-term without the KL (EFRET ~0.2) is stable, and high concentrations of Mg2+ or Gua+ only slightly affect the KL formation (Figure 3C). In the presence of 0.5 mM Mg2+, the structural change of riboG-term in response to increasing Gua+ was negligible (Figure 3D, Figure 3—figure supplement 4). Surprisingly, in the presence of 1.0 mM Gua+, the conformation dynamics of riboG-term became more evident as Mg2+ increased. Moreover, at high Mg2+ concentrations, such as 10.0 mM, the abundance of the folded conformation (EFRET ~0.8) significantly elevated with increasing Gua+ (Figure 3E, Figure 3—figure supplement 5). This suggests that Gua+ enhances the flexibility of riboG-term, potentially triggering its structural changes and facilitating the regulatory function of riboG in the presence of Mg2+. In contrast to riboG-term, both its G71C and C30G-G71C mutants displayed a reduced proportion of the state with EFRET ~0.8. Remarkably, the fractions of EFRET ~0.8 remained unaffected by the addition of 1.0 mM Gua+ in these mutants. Distinct from riboG-term, no structural transitions between states were observed in the two mutants (Figure 3—figure supplement 6). Regarding the U72C mutant of riboG-term, the mutation at the site 72 had a reduced impact on the KL conformation in the presence of 1.0 mM Gua+ and 2.0 mM Mg2+. However, the increased proportion of EFRET ~0.8 in the A29G-U72C mutant of riboG-term suggests that these mutations can restore the base-pairing between sites 29 and 72, as well as facilitate the formation of the KL (Figure 3—figure supplement 7).
In the presence of Mg2+ and/or Gua+, the FRET histogram of full-length riboG (Figure 4A, Figure 4—figure supplement 1) displayed a greater number of peaks compared to riboG-apt and riboG-term (Figure 4B, Figure 4—figure supplement 2). Additionally, the smFRET traces revealed structural transitions among the states, as shown in Figure 4—figure supplements 2–5. The presence of more structures in the full-length riboswitch compared to the aptamer domain is consistent with findings in other RNA molecules, such as the 2′-dG riboswitch (Helmling et al., 2017). Notably, upon the addition of 1.0 mM Gua+ in the presence of 2.0 mM Mg2+, the proportion of the new peak (EFRET ~0.4) increased from 19% to 29%, and then decreased as the Mg2+ concentrations increased from 2.0 to 10.0 mM (Figure 4—figure supplements 2D and 4B–D). In the absence of Gua+, the addition of 0–50.0 mM Mg2+ only slightly altered the proportion of the folded conformation (EFRET ~0.8) of the full-length riboG (Figure 4—figure supplement 2). Even in the presence of 1.0 mM Gua+, an increase in Mg2+ concentration from 0 to 10.0 mM had a slight effect on the proportion of the folded conformation (Figure 4—figure supplement 4B–D). However, the effect of Gua+ on full-length riboG in the presence of 10.0 mM Mg2+ was significantly more pronounced compared to its effect in the absence of Mg2+ (Figure 4C, Figure 4—figure supplement 5B-F). The impact of Mg2+ and Gua+ on the folded conformation of full-length riboG showed similarities to their effects on riboG-term but had distinct differences compared to their effects on riboG-apt (Figure 4D and E). Overall, the effects of Mg2+ and Gua+ on the structures of full-length riboG and riboG-term were weaker compared to their effects on riboG-apt. Upon comparing the G71C and C30G-G71C mutants of the full-length riboG with their wild-type counterpart, it was observed that the wild-type adopted higher proportions of the state with EFRET ~0.8 (Figure 4—figure supplement 6). Regarding the U72C and A29G-U72C mutants of the full-length riboG, their behaviors with regards to the peak with EFRET ~0.8 were similar to that of their counterparts in riboG-term (Figure 4—figure supplement 7).
Kinetic analysis of smFRET data to identify the folding pathway of the guanidine-IV riboswitch
The guanidine-IV riboswitch exhibits dynamic behavior, with its structural transitions in response to Mg2+ and Gua+ closely linked to transcription termination. The relative free energy (ΔΔG) was calculated based on the state populations in the FRET histogram (Figure 5A and B). For riboG-apt, both the pre-folded and folded states displayed higher ΔΔG values compared to the unfolded state, indicating that the former states were less stable in the presence of 2.0 mM Mg2+. However, the addition of 1.0 mM Gua+ successfully reduced the ΔΔG of the folded state to a level lower than that of the unfolded state (Figure 5A). Meanwhile, only slight decreases in ΔΔG were observed for riboG-term and riboG upon the addition of 1.0 mM Gua+ in the presence of 2.0 mM Mg2+ (Figure 5B). The dwell time distributions were obtained from the single-molecule trajectories of FRET data for riboG-apt, riboG-term, and riboG (Figure 5C–E, Figure 5—figure supplements 1–4). The exponential decays of the dwell time were fitted to a first-order kinetic equation, allowing for the determination of dwell time and rate constants for the individual states of riboG-apt (Figure 5C–E). Based on these results, we proposed folding pathways to illustrate the transcriptional fate of the guanidine-IV riboswitch under different conditions (Figure 5F, Figure 5—figure supplement 4). The analysis of the dwell time data from single-molecule trajectories indicated that the folding process of riboG-apt was influenced by Mg2+ and/or Gua+. In the presence of 2.0 mM Mg2+, there was inefficient structural switching from the unfolded riboG-apt (τ~0.69 s) to the pre-folded state (τ~0.11 s) or folded riboG-apt (τ~0.40 s), resulting in the production of the terminator during transcription (red arrows, Figure 5F). Conversely, in the presence of 1.0 mM Gua+, the primary switching of riboG-apt occurred between the unfolded (τ~0.60 s) and folded states (τ~1.04 s), which can subsequently be transcribed into the anti-terminator (blue arrows, Figure 5F).
smFRET analysis revealed the Gua+-sensitive transcription window of riboG
The sensitivity of Gua+ varied among the post-transcriptional riboG-apt, riboG-term, and full-length riboG sequences. The folding and interactions of RNA during transcription play a critical role in its control functions for transcriptional riboswitch. To mimic how nascent riboG responds to Gua+ and its regulatory effects during transcription, we utilized PLOR to control its transcription at a single-nucleotide pace. Simultaneously, we collected smFRET data from active transcription elongation complexes (ECs) containing nascent RNA of different lengths (Figure 6A,B, Figure 6—figure supplement 1). PLOR is a liquid–solid hybrid phase transcription process that allows for pausing transcription at specific positions by omitting certain nucleotide triphosphates (NTPs). This paused transcription can then be restarted by adding the necessary NTPs. This unique ‘pause–restart’ transcriptional mode enables the evaluation of instant changes in the environment, such as the addition or removal of Gua+ at specific steps, facilitating the monitoring of nascent RNA under varied conditions. To investigate the range of ligand-responsiveness in riboG, we employed PLOR to generate various transcriptional complexes with different RNA length, this allowed us to assess the impact of Gua+ on the guanidine-IV riboswitch throughout the transcription process. We monitored the behavior of active transcription complexes containing guanidine-IV riboswitches with varying lengths of 87–105 nucleotides (nt) using smFRET (Figure 6C–I, Figure 6—figure supplements 2–5). To distinguish between the different complexes, we named them based on their nucleotide length, for example, the EC-87 for the 87 nt complex. We collected smFRET histograms for EC-87, EC-88, EC-89, EC-91, EC-94, EC-96, and EC-105 in the presence of 0.5 mM Mg2+ and either 0 or 10.0 mM Gua+ (Figure 6C–I, Figure 6—figure supplements 2 and 3). Three distinct peaks (EFRET ~0.2, 0.5, and 0.8) were observed, with the high-FRET peak suggesting the presence of a newly folded RNA structure, potentially containing a KL, as the transcript size exceeded 87 nt. In the presence of 0.5 mM Mg2+, the proportion of the folded structure slightly decreased during transcription from 87 to 88 nt (from approximately 15% to 13%), but then sharply decreased to 3% at 89 nt. The addition of Gua+ significantly increased the proportion of the folded state for EC-88 compared to EC-87. For longer ECs, such as EC-89, EC-91, and EC-94, the proportion of the folded state remained relatively stable (Figure 6E–I, Figure 6—figure supplements 2 and 3). At 0.5 and 2.0 mM Mg2+, the proportions of the peaks changed differently, but EC-88 exhibited more dramatic changes than other ECs upon the addition of Gua+ (Figure 6J, Figure 6—figure supplements 2–5). These results suggest that the optimal formation of KL occurs during transcription of the 88 nt nascent RNA, and a drastic structural switch takes place between 88 and 89 nt in the presence of Mg2+. Interestingly, the addition of 1.0 mM Gua+ induced more pronounced changes in the majority of ECs at higher Mg2+ concentrations (Figure 6J, Figure 6—figure supplements 4 and 5). These findings suggest that the ECs exhibit a greater propensity for structural transitions upon the addition of Gua+ at high Mg2+ concentrations. Furthermore, under conditions of 0.5 mM Mg2+ and 10.0 mM Gua+, the proportions of folded structures (EFRET ~0.8) in EC-88 and EC-105 were approximately 36% and 14% lower, respectively, compared to those observed in the isolated post-transcriptional riboG-apt and riboG-term (Figure 6D,I, Figure 2—figure supplement 6 and Figure 3—figure supplement 4F). The T7 RNA polymerase (RNAP) sequestered about 8 nt of the nascent RNA, preventing the EC-88 construct from forming the P2 stem (Durniak et al., 2008; Huang and Sousa, 2000; Lubkowska et al., 2011; Tahirov et al., 2002; Wang et al., 2022; Yin and Steitz, 2002). Consequently, a pseudoknot structure potentially formed instead of the expected KL. This distinction may account for the observed heterogeneity between EC-88 and riboG-apt.
Gua+ modulated the transcription of the guanidine-IV riboswitch at specific sequences
To investigate the regulatory functions of the nascent riboswitch at different lengths during transcription, we performed a transcriptional termination assay of the guanidine-IV riboswitch using PLOR (Figure 7A, Figure 7—figure supplement 1). Specifically, we carefully designed eight PLOR experiments to temporarily pause the transcription at nucleotides 69, 77, 78, 83, 88, 91, 94, and 105, respectively, before transcribing the remaining RNA fragment in the last step, following strategies 1–8 (Figure 7B, Figure 7—figure supplement 1B and 1C). In the absence of Gua+, there was minimal change in transcriptional termination among strategies 1–7, but a sharp decrease of approximately 20% was observed for strategy 8 (Figure 7C and D). However, with the addition of Gua+ in the last step, there was little variation in transcriptional read-through of the guanidine-IV riboswitch among strategies 1–5, but a dramatic decrease in anti-termination from strategy 6 to strategy 8 (Figure 7C and D). The detailed procedures of strategies 1–8 were shown in Figure 7—figure supplement 1. Furthermore, we observed a significant increase in read-through upon the addition of Gua+ in our single-round transcription reactions when the RNA length reached 88 nt (Figure 7—figure supplement 1D). Additionally, a slight increase in read-through was observed at the pausing site 105, which is consistent with the smFRET data presented in Figure 6J. This indicates that Gua+ can still bind and induce slight structural changes in the guanidine-IV riboswitch, even if its full terminator sequence has been transcribed. Our findings from in vitro transcriptional termination assays and smFRET results, indicating that EC-91 exhibits a more stable fold with a longer stem, making it less sensitive to environmental changes compared to EC-88 (Figure 6F, Figure 7—figure supplement 1D).
-
Figure 7—source data 1
- https://cdn.elifesciences.org/articles/94706/elife-94706-fig7-data1-v1.zip
-
Figure 7—source data 2
- https://cdn.elifesciences.org/articles/94706/elife-94706-fig7-data2-v1.zip
-
Figure 7—source data 3
- https://cdn.elifesciences.org/articles/94706/elife-94706-fig7-data3-v1.xlsx
Discussion
The guanidine-IV riboswitch exhibits similarities to the guanidine-I riboswitch in gene regulatory mechanism, functioning as a transcriptional riboswitch. Structurally, it resembles the guanidine-II riboswitch through the formation of loop-loop interactions upon binding to guanidine (Battaglia and Ke, 2018; Huang et al., 2017a; Huang et al., 2017b; Lenkeit et al., 2020; Nelson et al., 2017; Reiss and Strobel, 2017; Salvail et al., 2020). In this study, we have demonstrated that both Mg2+ and Gua+ play critical roles in promoting structural folding and transcriptional anti-termination in the guanidine-IV riboswitch. Based on our findings, we propose a model illustrating the termination/anti-termination process of the guanidine-IV riboswitch as transcription progresses (Figure 8). In the presence of Mg2+, the proportion of the anti-terminator induced by Gua+ increases as transcription progresses up to EC-88. However, after EC-88, a significant decrease is observed up to EC-105. The elongation of transcription from EC-88 to EC-105 leads to impaired pseudoknot or KL formation and an increase in terminator abundance. This incremental switch from the anti-terminator to the terminator from EC-88 to EC-105 is influenced by the elongating P2 structure (Figure 8). The competition between the P2 and KL structures determines the proportions of the terminator and the anti-terminator in the guanidine-IV riboswitch. The stability of the terminator is enhanced by the presence of P2, and the increase in base pairs within P2 strengthens the terminator as transcription proceeds. On the other hand, the stability of the anti-terminator is dependent on the formation of pseudoknot or KL. The elongation of P2 shortens the flexible linker between P1 and P2, which creates steric hindrance for KL after EC-88. Furthermore, the addition of Gua+ facilitates the functional switch from the terminator to the anti-terminator for constructs ranging from EC-87 to EC-105. Here, we propose a narrow Gua+-sensitive transcription window for the guanidine-IV riboswitch, with a dramatic switch observed during the extension from EC-87 to EC-88 upon Gua+ addition (shown in a green box, Figure 8). It is worth noting that the ligand-sensitive sequence for the guanidine-IV riboswitch is located within the aptamer domain, spanning 88 nucleotides.
Our findings demonstrate that the guanidine-IV riboswitch possesses a distinct binding window that can be modulated by guanidine. The mechanism of anti-termination supports a model in which the formation of a long-range KL prevents the assembly of a rigid terminator. This distinguishes the guanidine-IV riboswitch from typical transcriptional riboswitches that rely on overlapping sequences and base-pair rearrangements for functional change (Strobel et al., 2019; Yadav et al., 2022; Zhao et al., 2017). Moreover, our results reveal the existence of at least three distinct states in the guanidine-IV riboswitch, including a folded state characterized by the presence of a KL that plays a critical role in ligand binding and transcriptional regulation. Importantly, the formation of the long-range KL in the aptamer domain is specifically responsive to guanidine, as other cations such as K+, Na+, and urea do not produce the same effect. While guanidine contributes to structural switching in the aptamer domain, its impact on the terminator or full-length guanidine-IV riboswitch is comparatively weak. Additionally, our data indicate that guanidine induces less significant changes once the RNA polymerase has transcribed the aptamer domain. It is worth noting that the RNAP utilized in our study is T7 RNAP, which exhibits distinct characteristics compared to bacterial RNAP in terms of transcriptional speed, dynamics, and interactions. However, Xue et al. have reported similarities between T7 and E. coli RNAP in the folding of nascent RNA. Additionally, Lou and Woodson have provided valuable insights into the co-transcriptional folding of the glmS ribozyme using T7 RNAP (Xue et al., 2023; Lou and Woodson, 2024).
In the transcriptional mimicking assay, the nascent RNA displayed a weaker response to guanidine compared to the post-transcriptional counterparts. The guanidine sensitivity of EC-87 indicates that guanidine regulation for the guanidine-IV riboswitch occurs before transcription of the aptamer domain is completed. Additionally, the sensitivity to guanidine gradually increases in the RNA within ECs until EC-88, after which it decreases. In other words, there exhibits a limited transcriptional window for the guanidine-IV riboswitch to dynamically respond to guanidine, and the inclination to form the terminator increases after transcribing 88 nucleotides, resulting in a decrease in the abundance of the KL-folded conformation. This aligns with previous findings that the transcription decision of a transcriptional riboswitch is made prior to terminator formation (Sherwood and Henkin, 2016; Wickiser et al., 2005). Nevertheless, the ligand-sensitive windows of riboswitches during transcription vary. In a study conducted by Helmling et al., 2017 using NMR spectroscopy, they proposed a broad transcriptional window for deoxyguanosine-sensing riboswitches, whereby the ligand binding capability gradually diminishes over several nucleotide lengths. However, more recent research by Binas et al., 2020 and Landgraf et al., 2022 on riboswitches sensing ZMP, c-di-GMP, and c-GAMP revealed a narrow window with a sharp transition in binding capability, even with transcript lengths differing by only one or three nucleotides. In line with the findings for the c-GAMP-sensing riboswitch, our study on the guanidine-IV riboswitch also demonstrated a sharp transition in binding capability with just a single nucleotide extension. Our data also suggest that guanidine can induce slight structural changes in the guanidine-IV riboswitch even after the terminator domain has been transcribed and stabilized. Similar observations have been reported for other riboswitches, although the transcriptional window sensitive to the ligand differs among riboswitches (Wickiser et al., 2005).
Materials and methods
DNA template preparation
Request a detailed protocolThe biotin labeled DNA templates were prepared by PCR as described earlier (Liu et al., 2018; Liu et al., 2015). The sequences of the DNA templates are listed in Supplementary file 1, table S1. The primers used in PCR were ordered from Sangon Biotech (Shanghai, China). The DNA products generated from PCR were purified by 12% denaturing PAGE (polyacrylamide electrophoresis) before immobilized on streptavidin-coated agarose beads (Smart-Life Science, Changzhou, China) in the buffer (40 mM Tris-HCl, 100 mM K2SO4, 6 mM MgSO4, pH 8.0) as described earlier (Liu et al., 2018). The solid-phase DNA templates were stored at 4 °C before usage.
Cy3Cy5-labeled riboG-apt, riboG-term and full-length riboG synthesis
Request a detailed protocolThe Cy3Cy5-labeled RNA used in smFRET were synthesized by PLOR as described earlier (Liu et al., 2018; Liu et al., 2015), and RNA sequences were depicted in Supplementary file 1, tables S1, S2. The schemes and reagents of synthesis of Cy3Cy5-riboG-apt, Cy3Cy5-riboG-term and Cy3Cy5-full-length riboG are listed in Supplementary file 1,tables S3, S4 and S5. For Cy3Cy5-riboG-apt synthesis, 10 μM DNA-beads gently shook with 10 μM T7 RNAP, 0.96 mM ATP, 0.96 mM GTP and 96 μM UTP in the buffer (40 mM Tris-HCl, 100 mM K2SO4, 6 mM MgSO4, 10 mM DTT, pH 8.0) at 37 °C for 15 min in the first step. The reaction mixture was filtered by solid-phase extraction and rinsed for 3~4 times in the buffer (40 mM Tris-HCl, 6 mM MgSO4, pH 8.0). And the bead-rinsing procedure was performed before the addition of NTPs in each step except noted. The residual steps were proceeded in the elongation buffer (40 mM Tris-HCl, 6 mM MgSO4, 10 mM DTT, pH 8.0) at 25 °C for 10 min with the addition of different NTP mixture at individual step as listed in Supplementary file 1, table S3. The Cy3Cy5-labeled riboG-apt was collected at the last step, purified by 12% denaturing PAGE and reversed phase HPLC loaded with C8 column (4.6*250 mm, Cat. No. EXL-122–2546 U, Phenomenex Luna, USA) as described earlier. The Cy3Cy5-riboG-term and Cy3Cy5-full-length riboG were synthesized as Cy3Cy5-riboG-apt except different NTPs were added. The purified RNA was stored at –80 °C before usage. T7 RNAP was utilized in the PLOR and in vitro transcription reactions except noted.
Preparation of transcription complexes containing RNA at different lengths
Request a detailed protocolThe RNA sequences in the transcription complexes (ECs) are listed in Supplementary file 1, table S6, and the ECs were collected at different steps of PLOR as listed in Supplementary file 1,table S7. For example, EC-87 was obtained as the following procedures: took 10 μL 10 μM DNA-beads from 200 μL of the PLOR reaction at step 14 and incubated with 1 μL 1 mM biotin at 25 °C for 20 min. The addition of biotin was to displace the biotin-labeled DNA template on the streptavidin-labeled beads, filtered by solid-phase exchange and collected the elution containing EC-87. EC-88, 89, 91, 94, 96, and 105 were obtained at steps15, 16, 17, 18, 20, and 21 similarly as EC-87 (Supplementary file 1, table S7).
smFRET experiments and data analysis
Request a detailed protocol0.5 μM Cy3Cy5-riboG-apt, Cy3Cy5-riboG-term or Cy3Cy5-full-length riboG were hybridized with 0.75 μM biotinylated 13 nt-DNA by incubating at 90 °C for 5 min and then cooled to room temperature in the T50 buffer (10 mM Tris-HCl, 50 mM NaCl, pH 8.0). The sequence of the biotinylated DNA is listed in Supplementary file 1, table S1. The hybridized Cy3Cy5-labeled RNA was then diluted to 10~100 pM for the subsequent smFRET measurements. Polymer-coated quartz slides were immersed at the 0.1 mg/μL biotin-PEG and 0.3 mg/μL m-PEG for 2 h (Roy et al., 2008), then coated with 0.2 mg/mL streptavidin for 5 min. 10~100 pM Cy3Cy5-labeled RNA or ECs in the buffer (10 mM Tris-HCl, 50 mM NaCl, 0~100 mM guanidine, 0~100 mM Mg2+) was injected and then rinsed 3 times to remove the free RNA or ECs. 3 mM Trolox (MCE, USA), 5 mM 3,4-protocatechuic acid (PCA, Shanghai Aladdin Bio-Chem Technology Co., China) and 5 nM protocatechuate dioxygenase (PCD) was added to reduce the photobleaching (Aitken et al., 2008). smFRET experiments were performed by using an objective-type total internal reflection fluorescence (TRIF) microscopy and an inverted microscope (Eclipse Ti, Nikon, Japan) at 20 °C, with FRET marker dyes Cy3 (donor) and Cy5 (acceptor) was excited by 532 nm and 640 nm laser, respectively. Open-sourced software iSMS (Preus et al., 2015) was used to process the single-molecule videos to extract the time-dependent signals. Time resolution of 100ms was used for each single-molecule video collection excepted noted. Single-molecule trajectories were identified by Deep FRET (Thomsen et al., 2020). The FRET efficiency, EFRET was calculated using the equation:
where IA and ID are acceptor and donor fluorescence intensity, respectively. We used empirical Bayesian method of the hidden Markov modelling (HMM) package vbFRET to estimate the FRET states and the time points of transitions (Bronson et al., 2009). The transition density plots (TDP) were plotted by Python module matplotlib, and the lifetimes were plotted in histogram and fitted with a simple exponential decay function by Origin 8.5 to obtain the rate constants. The relative free energy, ΔΔG was measured using the equation ΔΔGab = −RT ln(Pb/Pa), where R, T, Pa and Pb are the gas constant, absolute temperature, population of the observed state b and population of the control state a, respectively (Manz et al., 2017).
In vitro transcription termination assay
Request a detailed protocolThe single-round transcription termination assay was carried out by PLOR (Chien et al., 2023). The detailed procedure and reagents are listed in Supplementary file 1,tables S8–16. 0~ 20mM guanidine and 6 mM Mg2+ were added at step 8 in Supplementary file 1, table S8 to measure the effect of guanidine on transcription read-through of the guanidine-IV riboswitch. The products collected at the last step were loaded to 12% denaturing PAGE to compare the transcription read-through efficiencies at different conditions.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Guanidine-sensing riboswitches: How do they work and what do they regulate?Wiley Interdisciplinary Reviews RNA 9:e1482.https://doi.org/10.1002/wrna.1482
-
The biochemical landscape of riboswitch ligandsBiochemistry 61:137–149.https://doi.org/10.1021/acs.biochem.1c00765
-
Quantitative analysis of transcription termination via Position-Selective Labeling of RNA (PLOR) methodInternational Journal of Molecular Sciences 24:4934.https://doi.org/10.3390/ijms24054934
-
Riboswitches: structures and mechanismsCold Spring Harbor Perspectives in Biology 3:a003533.https://doi.org/10.1101/cshperspect.a003533
-
NMR structural profiling of transcriptional intermediates reveals riboswitch regulation by metastable RNA conformationsJournal of the American Chemical Society 139:2647–2656.https://doi.org/10.1021/jacs.6b10429
-
Mimicking co-transcriptional RNA bolding using a superhelicaseJournal of the American Chemical Society 140:10067–10070.https://doi.org/10.1021/jacs.8b03784
-
T7 RNA polymerase elongation complex structure and movementJournal of Molecular Biology 303:347–358.https://doi.org/10.1006/jmbi.2000.4150
-
The structure of the guanidine-II riboswitchCell Chemical Biology 24:695–702.https://doi.org/10.1016/j.chembiol.2017.05.014
-
Structure of the guanidine III riboswitchCell Chemical Biology 24:1407–1415.https://doi.org/10.1016/j.chembiol.2017.08.021
-
Long-range interactions in riboswitch control of gene expressionAnnual Review of Biophysics 46:455–481.https://doi.org/10.1146/annurev-biophys-070816-034042
-
Discovering riboswitches: the past and the futureTrends in Biochemical Sciences 48:119–141.https://doi.org/10.1016/j.tibs.2022.08.009
-
Discovery and characterization of a fourth class of guanidine riboswitchesNucleic Acids Research 48:12889–12899.https://doi.org/10.1093/nar/gkaa1102
-
Co-transcriptional folding of the glmS ribozyme enables a rapid response to metaboliteNucleic Acids Research 52:872–884.https://doi.org/10.1093/nar/gkad1120
-
RNA folding in transcription elongation complexJournal of Biological Chemistry 286:31576–31585.https://doi.org/10.1074/jbc.M111.249359
-
Gene regulation by riboswitchesNature Reviews Molecular Cell Biology 5:451–463.https://doi.org/10.1038/nrm1403
-
Single-molecule FRET reveals the energy landscape of the full-length SAM-I riboswitchNature Chemical Biology 13:1172–1178.https://doi.org/10.1038/nchembio.2476
-
iSMS: single-molecule FRET microscopy softwareNature Methods 12:593–594.https://doi.org/10.1038/nmeth.3435
-
A practical guide to single-molecule FRETNature Methods 5:507–516.https://doi.org/10.1038/nmeth.1208
-
Riboswitch-mediated gene regulation: novel RNA architectures dictate gene expression responsesAnnual Review of Microbiology 70:361–374.https://doi.org/10.1146/annurev-micro-091014-104306
-
A ligand-gated strand displacement mechanism for ZTP riboswitch transcription controlNature Chemical Biology 15:1067–1076.https://doi.org/10.1038/s41589-019-0382-7
-
Short oligonucleotides facilitate co-transcriptional labeling of RNA at specific positionsJournal of the American Chemical Society 144:5494–5502.https://doi.org/10.1021/jacs.2c00020
-
Genetic control by metabolite-binding riboswitchesChembiochem 4:1024–1032.https://doi.org/10.1002/cbic.200300685
-
An excited state underlies gene regulation of a transcriptional riboswitchNature Chemical Biology 13:968–974.https://doi.org/10.1038/nchembio.2427
-
Mfold web server for nucleic acid folding and hybridization predictionNucleic Acids Research 31:3406–3415.https://doi.org/10.1093/nar/gkg595
Article and author information
Author details
Funding
National Natural Science Foundation of China (32071300)
- Yu Liu
National Natural Science Foundation of China (31872628)
- Yu Liu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The work was supported by the National Natural Science Foundation of China [grant no. 32071300 and 31872628].
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.94706. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Gao et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 203
- views
-
- 20
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins, combined with the exposure of their residues, accounts for this sensitivity. One context in which IDPs play important roles that are concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat-soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet the mechanisms underlying this synergy differ between IDP families.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Dynamic conformational and structural changes in proteins and protein complexes play a central and ubiquitous role in the regulation of protein function, yet it is very challenging to study these changes, especially for large protein complexes, under physiological conditions. Here, we introduce a novel isobaric crosslinker, Qlinker, for studying conformational and structural changes in proteins and protein complexes using quantitative crosslinking mass spectrometry. Qlinkers are small and simple, amine-reactive molecules with an optimal extended distance of ~10 Å, which use MS2 reporter ions for relative quantification of Qlinker-modified peptides derived from different samples. We synthesized the 2-plex Q2linker and showed that the Q2linker can provide quantitative crosslinking data that pinpoints key conformational and structural changes in biosensors, binary and ternary complexes composed of the general transcription factors TBP, TFIIA, and TFIIB, and RNA polymerase II complexes.