Specific structural elements of the T-box riboswitch drive the two-step binding of the tRNA ligand

  1. Jiacheng Zhang
  2. Bhaskar Chetnani
  3. Eric D Cormack
  4. Dulce Alonso
  5. Wei Liu
  6. Alfonso Mondragón
  7. Jingyi Fei  Is a corresponding author
  1. University of Chicago, United States
  2. Northwestern University, United States
  3. University of Chicago College, United States

Abstract

T-box riboswitches are cis-regulatory RNA elements that regulate the expression of proteins involved in amino acid biosynthesis and transport by binding to specific tRNAs and sensing their aminoacylation state. While the T-box modular structural elements that recognize different parts of a tRNA have been identified, the kinetic trajectory describing how these interactions are established temporally remains unclear. Using smFRET, we demonstrate that tRNA binds to the riboswitch in two steps, first anticodon recognition followed by the sensing of the 3’ NCCA end, with the second step accompanied by a T-box riboswitch conformational change. Studies on site-specific mutants highlight that specific T-box structural elements drive the two-step binding process in a modular fashion. Our results set up a kinetic framework describing tRNA binding by T-box riboswitches, and suggest such binding mechanism is kinetically beneficial for efficient, co-transcriptional recognition of the cognate tRNA ligand.

https://doi.org/10.7554/eLife.39518.001

eLife digest

Living organisms depend upon a group of chemicals called amino acids to survive. Amino acids are the building blocks of proteins, and proteins have many important roles within and around cells. Bacteria regulate certain genes to ensure they have the right balance of different amino acids to survive. By controlling the availability of certain proteins that help them to make or collect certain amino acids, bacteria can control their overall amino acid balance.

Before a protein is made, a molecular machine called RNA polymerase must first copy the information in a gene to make a molecule called a messenger RNA (mRNA). The mRNA is then translated to make the protein from individual amino acids. In this process, each amino acid needs to be first attached to another molecule called a transfer RNA (tRNA). In many bacteria species, the mRNAs involved in making or transporting amino acids contain structures called T-boxes. These structures guide the RNA polymerase to make more of the mRNAs when the levels of the amino acid become too low. A T-box, however, does not sense the level of the amino acid directly. Instead it senses the number of tRNA molecules that do not carry an amino acid.

Zhang, Chetnani et al. examined a particular T-box interacting with tRNA using pairs of fluorescent dyes to detect distances between molecules. The T-box first recognizes a part of the tRNA called the anticodon to make sure it binds the correct type of tRNA. It then changes its shape to detect whether the tRNA is attached to an amino acid. This two-step process is driven by multiple structural elements within the T-box, and the flexibility of the T-box plays a critical role.

A cell’s survival depends on it keeping amino acid levels under control. Understanding how bacteria do this could lead to new antibiotic drugs that target the T-box to kill cells. This study also provides insights into the workings of mRNA components like T-boxes – a type of riboswitch – which is an unusual means of controlling gene activity.

https://doi.org/10.7554/eLife.39518.002

Introduction

Riboswitches are cis-regulatory RNA elements that recognize and respond to defined external signals to affect transcription or translation of downstream messenger RNAs (mRNAs) (Breaker, 2012; Serganov and Nudler, 2013; Sherwood and Henkin, 2016). Riboswitches generally consist of two domains: a sensory or aptamer domain and a regulatory domain or expression platform. The expression platform can adopt different conformations in response to ligand binding to the aptamer, and in this way control gene expression outcome (Breaker, 2012). The aptamer of each riboswitch class contains conserved sequence motifs and unique secondary or tertiary structural elements that help distinguish and bind specific ligands (McCown et al., 2017). Bacterial T-box riboswitches represent a unique class of riboswitches that do not bind small molecule ligands, instead they recognize and bind tRNA molecules and sense directly their aminoacylation state (Zhang and Ferré-D'Amaré, 2015). T-box riboswitches serve as excellent paradigms to understand RNA-RNA interactions and RNA-based regulation.

T-box riboswitches are found in Gram-positive bacteria and are usually located in the region upstream of mRNA sequences encoding aminoacyl tRNA synthetases and proteins involved in amino acid biosynthesis and transport and hence participate directly in amino acid homeostasis (Zhang and Ferré-D'Amaré, 2015). In general, the aptamer domain of all T-box riboswitches contains a long stem, Stem I, responsible for specific tRNA binding (Rollins et al., 1997). The expression platform can adopt either a terminator or anti-terminator conformation, depending on whether the bound tRNA is charged or uncharged (Henkin, 2014; Zhang and Ferré-D'Amaré, 2015). In most T-box riboswitches, binding of a charged tRNA to the T-box leads to rho independent transcription termination whereas an uncharged tRNA stabilizes the anti-terminator conformation and leads to transcription read-through (Henkin, 2014; Zhang and Ferré-D'Amaré, 2015). Whereas Stem I and the anti-terminator domain are highly conserved among T-box riboswitches, the region connecting them can vary. The Bacillus subtilis glyQS T-box riboswitch, involved in glycine regulation, represents one of the simplest T-box riboswitches (Grundy et al., 2002b) in which only a short linker and a small stem, Stem III, connect Stem I and the anti-terminator domain (Figure 1).

Figure 1 with 2 supplements see all
Secondary and tertiary structures of B. subtilis glyQS T-box riboswitch and tRNAGly.

(A) Secondary structure diagrams of the B. subtilis glyQS T-box riboswitch and B. subtilis tRNAGly used in this study. Green and orange lines indicate interactions between the T-box specifier loop and the tRNA anticodon and between the T-box t-box sequence and the tRNA 3’ NCCA, respectively. For the glyQS T-box sequence, the nucleotides in red were added for surface immobilization. (B) Ribbon diagram of a model of a complex between the B. subtilis glyQS T-box riboswitch (blue) and B. subtilis tRNAGly (green) based on SAXS data (Chetnani and Mondragón, 2017). Distances between the 5’ and 3’ ends of the T-box and the 5’ end of the tRNAGly are shown (black dash lines). The NCCA sequence at the 3’ end of the tRNA is shown in light green and the t-box sequence in the T-box is shown in yellow.

https://doi.org/10.7554/eLife.39518.003

Recognition of tRNA by a T-box riboswitch involves three main structural elements of the tRNA: the anticodon region, the ‘elbow’ region formed by the conserved T- and D-loops, and the 3’ NCCA sequence (Figure 1). The anticodon and elbow regions of the tRNA interact with Stem I directly. Stem I contains several phylogenetically conserved structural motifs (Rollins et al., 1997), including a K-turn motif, a specifier loop, a distal bulge, and an apical loop (Rollins et al., 1997) (Figure 1). Bioinformatics and structural analyses have collectively revealed the interactions between Stem I and the tRNA (Grigg et al., 2013; Lehmann et al., 2013; Zhang and Ferré-D'Amaré, 2013). Specifically, the co-crystal structures of Stem I/tRNA complexes show that Stem I flexes to follow closely the tRNA anticodon stem and interacts directly with the anticodon loop and the elbow through its proximal and distal ends, respectively (Zhang and Ferré-D'Amaré, 2013). The distal bulge and the apical loop fold into a compact structural module of interdigitated T-loops (Chan et al., 2013; Krasilnikov and Mondragón, 2003), which interact directly with conserved unstacked nucleobases at the tRNA elbow (Grigg et al., 2013; Zhang and Ferré-D'Amaré, 2013). In addition, the structures revealed that Stem I turns sharply around two hinge regions using a conserved dinucleotide bulge and the K-turn motif (Grigg and Ke, 2013; Zhang and Ferré-D'Amaré, 2013). Sensing of the aminoacylation state involves direct binding of the tRNA 3’ end to a highly conserved bulge in the T-box, the t-box sequence (Grundy et al., 2002a) (Figure 1). A free NCCA end can base pair with the t-box sequence, enabling the anti-terminator conformation, whereas a charged NCCA end prevents the formation of the NCCA/t-box interactions, leading to the more stable terminator conformation (Henkin, 2014; Zhang and Ferré-D'Amaré, 2015). Importantly, discrimination between the charged and uncharged tRNA does not require any additional proteins, such as EF-Tu (Suddala et al., 2018; Zhang and Ferré-D'Amaré, 2014), and is driven solely by RNA/RNA interactions.

While Small Angle X-ray Scattering (SAXS)-derived models of the entire B. subtilis glyQS T-box riboswitch in complex with tRNA are available (Chetnani and Mondragón, 2017; Fang et al., 2017), atomic-level structural details on the interactions between tRNA and the anti-terminator region are still lacking. In addition, there is a dearth of information on the kinetics of the binding process. Whereas it is clear that tRNA recognition involves several specific interactions, their binding temporal sequence remains elusive. In addition, it is unclear whether sensing of the 3’ end of the tRNA involves any additional conformational changes in the T-box. Here, by introducing donor-acceptor fluorophore pairs at several locations in the tRNA and the T-box riboswitch, and using single-molecule fluorescence resonance energy transfer (smFRET), we demonstrate the temporal order of events in the trajectory of tRNA binding. Our results demonstrate that tRNA binds to the riboswitch in two steps, with its anticodon being recognized first, followed by NCCA binding accompanied with an inward motion of the 3’ region of the T-box riboswitch, including Stem III and the anti-terminator stem, relative to Stem I. In addition, by introducing mutations at different locations of the T-box, we further show that the two-step binding kinetics is regulated by the modular structural elements in the T-box riboswitch.

Results

Binding of cognate tRNA by the glyQS T-box results in two distinct FRET states

To observe directly the binding of tRNA to the T-box, we placed the donor dye (Cy3) on the 3’ end of a T-box fragment (T-box182), and the acceptor dye (Cy5) on the 5’ end of the tRNAGly, where the subscript ‘182’ denotes the length of the T-box construct (Figure 1A). In vitro transcribed and labeled T-box and tRNA were purified and refolded according to published procedures (Chetnani and Mondragón, 2017; Zhang and Ferré-D'Amaré, 2013) (Figure 1—figure supplement 1). Labeling of the tRNA at the 5’ end had a modest effect on the binding affinity (Figure 1—figure supplement 2). T-box182 spans Stem I, the linker sequence, Stem III and the anti-terminator, but does not contain the terminator sequence, thereby preventing the transition to the terminator conformation. A short RNA extension sequence was added to the 5’ end of the T-box for surface immobilization (Figure 2A, Supplementary file 1). Single-molecule fluorescence images were recorded under equilibrium condition in the presence of 30 nM tRNAGly-Cy5. Binding of tRNAGly-Cy5 results in a major distribution of FRET values around 0.7, with 79 ± 4% of the traces showing a stable signal at 0.7 and 9 ± 5% traces sampling from 0.7 to 0.4 (Figure 2B,C). The SAXS model (Chetnani and Mondragón, 2017) predicts a distance between the labeling positions at the 3’ end of the T-box182 and the 5’ end of the tRNA to be around 52 Å. (Figure 1B). Based on a Förster distance of 54–60 Å (Ha et al., 2002; Hohng et al., 2004), our measured FRET value is within the range of estimated FRET values (0.56–0.70). Therefore, we assign the 0.7 FRET state to be the fully bound state of the tRNAGly by the T-box.

Figure 2 with 4 supplements see all
Two-step binding of uncharged tRNA to the glyQS T-box riboswitch.

(A) FRET labeling scheme for the T-box and tRNA. Cy3 (green star) and Cy5 (red star) fluorophores are attached at the 3’ of the T-box (blue) and the 5’ of the tRNA (black), respectively. glyQS T-box riboswitch molecules are anchored on slides through a biotinylated DNA probe (purple) hybridized to a 5’ extension sequence on the T-box. (B) smFRET vs. time trajectories of T-box182-Cy3(3’) with tRNAGly-Cy5, T-box182-Cy3(3’) with tRNAΔNCCA-Cy5 and T-box149-Cy3(3’) with tRNAGly-Cy5. Cy3 and Cy5 fluorescence intensity traces (upper panel), and their corresponding smFRET traces calculated as ICy5 / (ICy3+ICy5) (lower panel). (C) One-dimensional FRET histograms. FRET peaks are fit with a Gaussian distribution (black curve) and the peak centers are shown in red. ‘N’ denotes the total number of traces in each histogram from three independent experiments. (D) Transition density plot (TDP). Contours are plotted from white (less than 15% of the maximum population) to red (more than 85% of the maximum population). TDPs are generated from all smFRET traces from three independent experiments. (E) Representative smFRET trajectories showing real-time binding of tRNAGly-Cy5 to T-box182-Cy3(3’) in a steady-state measurement. Traces showing transitions from the unbound state (0 FRET) to fully bound state (0.7 FRET) through the partially bound state (0.4 FRET) (left) and unbound state directly to fully bound state (right). (F) Surface contour plot of time-evolved FRET histogram of T-box182-Cy3(3’) with tRNAGly-Cy5 (top) and tRNAΔNCCA-Cy5 (bottom). Contours are plotted from blue (less than 5% of the maximum population) to red (more than 75% of the maximum population). ‘N’ denotes the total number of traces in each histogram from three independent experiments, which are a subset of traces showing real-time binding events in the steady-state measurements. Total numbers of traces in steady-state measurements are indicated in (C). Traces that reach the 0.7 FRET state (cutoff >0.55) are included in the plot for tRNAGly-Cy5 to reveal better the transition from the 0.4 to the 0.7 FRET state. Time-evolved FRET histograms of all traces are shown in Figure 2—figure supplement 4D for comparison.

https://doi.org/10.7554/eLife.39518.006
Figure 2—source data 1

Cy3 and Cy5 intensity traces from three repeats for T-box182-Cy3(3’) and tRNAGly-Cy5.

Each two columns show intensities of the Cy3/Cy5 pair from the same molecule.

https://doi.org/10.7554/eLife.39518.014
Figure 2—source data 2

Cy3 and Cy5 intensity traces from three repeats for T-box182-Cy3(3’) and tRNAΔNCCA-Cy5.

https://doi.org/10.7554/eLife.39518.015
Figure 2—source data 3

Cy3 and Cy5 intensity traces from three repeats for T-box149-Cy3(3’) and tRNAGly-Cy5.

https://doi.org/10.7554/eLife.39518.016

In order to assign the 0.4 FRET value to specific tRNA binding states, tRNATyr-Cy5 and tRNAΔNCCA-Cy5 (‘ΔNCCA’ denotes a tRNAGly with a deleted 3’ NCCA sequence) were flowed in the flow-chamber with pre-immobilized T-box182-Cy3(3’) (3’ denotes that the label was added at the 3’ end). We did not observe any binding of tRNATyr-Cy5 (Figure 2—figure supplement 1), confirming that recognition of the anticodon by the specifier region is required for tRNA binding. In the presence of tRNAΔNCCA-Cy5, we observed a fluctuating signal between 0.4 and 0 FRET (Figure 2B,D), with a mean lifetime of the 0.4 FRET state of 3.6 ± 0.6 s and a mean waiting time before binding of 31.3 ± 5.3 s (Figure 2—figure supplement 2B). Taken together with the results from the tRNAGly, tRNAΔNCCA and tRNATyr binding experiments, we assign the 0.4 FRET state to a partially bound state where only the anticodon interactions have been established.

To further confirm the assignment of the FRET states, we generated T-box149, where the anti-terminator sequence is truncated (Figure 1A, Supplementary file 1). Based on the structure model from the SAXS data (Chetnani and Mondragón, 2017) we predicted that a Cy3 dye placed either at the end of Stem III (T-box149) or at the end of the anti-terminator stem (T-box182) are localized in close proximity in three dimensions, further confirmed by the distance measurement using smFRET (Figure 2—figure supplement 3). Therefore, we expect that if tRNAGly-Cy5 can reach the same fully bound state in T-box149 as in T-box182, a high FRET state centered at 0.7 would be observed. However, using T-box149-Cy3(3’) in combination with tRNAGly-Cy5, we again observed transient binding of tRNAGly with a FRET value centered at ~0.4 with the same mean lifetime as observed with the T-box182-Cy3(3’) and tRNAΔNCCA-Cy5 combination (Figure 2B–D, Figure 2—figure supplement 2C). Therefore these two complexes (T-box182 + tRNAΔNCCA and T-box149 + tRNAGly) represent the same binding state of the tRNA, that is the state where binding of the anticodon to the specifier region has been established, but is unstable without the further interactions between the NCCA and the t-box region.

Collectively, our results suggest a two-step binding model involving the separate establishment of the interactions with the anticodon and the NCCA. The fact that tRNATyr, which has a mismatched anticodon, but contains an intact NCCA 3’ end, does not show any binding activity suggests that interactions with the anticodon precede the interactions with the NCCA end of the tRNA and are necessary for the establishment of the NCCA contacts. Without the interaction between the NCCA and the t-box sequence the binding of tRNAGly is not stable. From the binding kinetics of tRNAΔNCCA, we estimated the association rate constant (k1) and the disassociation rate constant (k-1) for the first binding step to be (5.0 ± 1.6) x 105 M−1s−1 and 0.28 ± 0.05 s−1, respectively (Figure 5; Figure 6E).

The transition from anticodon recognition to NCCA binding is rapid for uncharged tRNA

We classified smFRET traces for T-box182-Cy3(3’) in complex with tRNAGly-Cy5 into three types (Figure 2B): (I) traces stably sampling the 0.7 state (79 ± 4% of total traces), (II) traces transiently transitioning from the 0.7 state to the 0.4 state (9 ± 5%), and (III) traces only sampling the 0.4 state without reaching 0.7 state (12 ± 5%). The low percentage of Type III traces indicates that once the anticodon is recognized, the commitment to the next binding step, NCCA interactions, is high. The majority of the traces showed that the tRNAGly remained mostly in the fully bound state (Type I) until the fluorophore photobleached, with the actual lifetime limited by photobleaching (τ0.7 > 24 s, where τ0.7 denotes the lifetime of the 0.7 FRET state) (Figure 2—figure supplement 2A). The observation that tRNAGly is able to transit from the fully bound state back to the partially bound state (Type II) suggests that the NCCA/t-box interaction can break occasionally (Figure 2B). We estimated the lifetime of the transiently sampled partially bound state to be 0.35 ± 0.09 s in the presence of full-length tRNAGly (Figure 2—figure supplement 2A), ~10 fold shorter than the partially bound state in the presence of tRNAΔNCCA.

While the majority of the T-box molecules were already bound to tRNAGly before starting data acquisition, we could detect that some molecules show real-time binding during imaging acquisition. We observed only a few traces briefly sampling the 0.4 FRET state from the zero FRET (unbound) state before reaching the 0.7 FRET state, while most traces directly sampled the 0.7 FRET state without a detectable 0.4 FRET, likely due to our imaging time resolution (100 ms per frame). We post-synchronized the FRET traces at the transition point from the zero FRET state to the first sampled 0.4 FRET state, and plotted them in a time-evolved FRET histogram. From the time-evolved FRET histogram (Figure 2F), we estimated roughly that the upper limit of the lifetime spent at the 0.4 FRET state is ~100 ms, very rapidly followed by establishment of NCCA/t-box interactions. In contrast, tRNAΔNCCA could not pass the 0.4 FRET state. To capture better real-time binding, we performed a flow experiment, where tRNAGly-Cy5 was flowed into a chamber with immobilized T-box182-Cy3(3’) during imaging acquisition. The corresponding post-synchronized time-evolved FRET histogram again shows a fast transition into the fully bound state (Figure 2—figure supplement 4). In addition, the association rate constant of tRNAGly in the real-time flow experiment is (7.5 ± 0.7) x 105 M−1s−1, consistent with the k1 of tRNAΔNCCA and confirming that the NCCA end of the tRNA does not participate in the first binding step.

From the real-time binding kinetics of tRNAGly to T-box182, we estimated a transition rate constant from the partially bound state to the fully bound state (k2) to be on the order of 10 s−1 (Figure 2F, Figure 2—figure supplement 4). On the other hand, as transitions back to the partially bound state from the fully bound state were only observed in ~10% traces, we interpreted this to mean that the reverse transition rate constant (k-2) is very small, and the second binding step in the wild-type (WT) T-box with uncharged tRNAGly is close to irreversible (Figure 5; Figure 6E and see Discussion).

Establishment of the NCCA/t-box interaction is accompanied by conformational changes in the T-box riboswitch

We next investigated whether tRNA binding requires any conformational changes in the T-box itself. Using doubly labeled T-box182, with Cy3 at the 3’ end and Cy5 at the 5’ hybridization extension, we observed a high FRET state (centered at ~0.75) in the absence of tRNA (Figure 3—figure supplement 1). Based on the structural model (Chetnani and Mondragón, 2017), we estimated the distance between the 5’ and 3’ ends of the T-box182 to be ~36 Å (Figure 1B). Our measured FRET value is slightly less than the predicted FRET value (~0.90), likely due to the engineered 5’ extension sequence used to immobilize the T-box. No noticeable change was detected upon incubation with unlabeled tRNAGly (Figure 3—figure supplement 1), indicating that the 3’ portion (Stem III plus the anti-terminator stem) does not move away from the 5’ portion (Stem I). Given that the measured FRET efficiency of 0.75 is already located beyond the FRET sensitive region, it is unlikely that any inward motion of the 3’ portion relative to the 5’ could be detected. To overcome this limitation, we added extensions at both the 3’ and 5’ ends (Figure 3A, Supplementary file 1). ITC experiments suggest that addition of a 5’ and/or a 3’ extension sequences to the T-box does not affect tRNA binding (Figure 1—figure supplement 2). With this intra-T-box FRET scheme, we observed a FRET shift from ~0.5 to~0.65 when tRNAGly was added (Figure 3B), indicating that the 3’ half of the T-box moves closer to the 5’ half, potentially with the T-box becoming more compact due to the presence of the cognate tRNAGly. Adding non-cognate tRNAPhe or tRNAΔNCCA gave similar FRET values as the T-box alone (Figure 3B), suggesting that the conformational change is associated with binding of both the anticodon and NCCA, not with anticodon recognition alone (Figure 5).

Figure 3 with 1 supplement see all
Conformational changes following tRNA binding in the glyQS T-box riboswitch revealed by an intra-T-box FRET pair.

(A) Intra-T-box FRET scheme. Cy3 (green star) and Cy5 (red star) are attached at the 5’ and 3’ extensions of T-box (blue), respectively. (B) One-dimensional FRET histograms of T-box182 alone, with tRNAPhe, with tRNAΔNCCA, and with tRNAGly. ‘N’ denotes the total number of traces in each histogram from three independent experiments.

https://doi.org/10.7554/eLife.39518.017
Figure 3—source data 1

Cy3 and Cy5 intensity traces from three repeats for intra T-box182 FRET with 3’ extension without tRNA.

https://doi.org/10.7554/eLife.39518.020
Figure 3—source data 2

Cy3 and Cy5 intensity traces from three repeats for intra T-box182 FRET with 3’ extension in the presence of tRNAPhe.

https://doi.org/10.7554/eLife.39518.021
Figure 3—source data 3

Cy3 and Cy5 intensity traces from three repeats for intra T-box182 FRET with 3’ extension in the presence of tRNAΔNCCA.

https://doi.org/10.7554/eLife.39518.022
Figure 3—source data 4

Cy3 and Cy5 intensity traces from three repeats for intra T-box182 FRET with 3’ extension in the presence of tRNAGly.

https://doi.org/10.7554/eLife.39518.023

The NCCA end of the uncharged tRNA maintains its relative position to the K-turn region during the second binding step

Using the above two FRET pairs, we observed that the 3’ portion of the T-box moves towards the base of Stem I as well as the NCCA end of the tRNA during the second binding step. To ascertain whether the NCCA end of the tRNA also moves relative to the base of Stem I, we measured FRET between a Cy3 placed at the 5’ end of the T-box (T-box182-Cy3(5’)) and tRNAGly-Cy5 (Figure 4A). Using this FRET pair, binding of both tRNAGly and tRNAΔNCCA generated a similar FRET value centered at ~0.35 (Figure 4B,C). However, the FRET traces behaved differently for these two tRNA molecules. For tRNAΔNCCA, the signal fluctuated between zero and 0.35 (Figure 4D), with a lifetime of the 0.35 FRET state of 4.5 ± 1.0 s, reminiscent of the 0.4 FRET state using the tRNA/T-box182-Cy3(3’) FRET pair (Figure 4—figure supplement 1). For tRNAGly, the signal was more stably centered at 0.35 (Figure 4B). Since the tRNA-Cy5/T-box182–Cy3(5’) FRET pair cannot distinguish the partially bound from the fully bound state, we fit the lifetime with a double-exponential decay. The fast dissociation fraction has a lifetime of 3.9 ± 0.7 s (46 ± 21% of population), consistent with the lifetime for the partially bound state, and the low dissociation fraction has a lifetime of 15.7 ± 0.8 s (54 ± 21%), representing the stable fully bound state (Figure 4—figure supplement 1). Overall, the measurements with the tRNA-Cy5/T-box–Cy3(5’) FRET pair further validate the two-step binding model and reveal that the NCCA end of the uncharged tRNA maintains its relative position to the base of Stem I during the second binding step.

Figure 4 with 1 supplement see all
FRET between fluorophores at the 5’ end of the glyQS T-box riboswitch and 5’ end of tRNAGly is insensitive to the two binding states.

(A) Cy3 (green star) and Cy5 (red star) are attached at the 5’ extension of the T-box (blue) and the 5’ of the tRNA (black), respectively. (B) smFRET trajectories of T-box-Cy3(5’) with tRNAΔNCCA-Cy5 (left) and tRNAGly-Cy5 (right). ‘N’ denotes the total number of traces in each histogram from three independent experiments. (C) One-dimensional FRET histograms of T-box-Cy3(5’) with tRNAΔNCCA-Cy5 (left) and tRNAGly-Cy5 (right). (D) TDP of T-box-Cy3(5’) with tRNAΔNCCA-Cy5 (left) and tRNAGly-Cy5 (right). Contours are plotted in the same way as in Figure 2D. TDPs are generated from all smFRET traces from three independent experiments.

https://doi.org/10.7554/eLife.39518.024
Figure 4—source data 1

Cy3 and Cy5 intensity traces from three repeats for T-box182-Cy3(5’) with tRNAGly-Cy5.

https://doi.org/10.7554/eLife.39518.026
Figure 4—source data 2

Cy3 and Cy5 intensity traces from three repeats for T-box182-Cy3(5’) with tRNAΔNCCA-Cy5.

https://doi.org/10.7554/eLife.39518.027

A mutation in the T-loop region affects the first binding step but has minimal effect on the second binding step

The interdigitated T-loops structure formed by the interactions between the distal bulge and the apical loop at the distal end of Stem I has been shown to be important for tRNA binding (Grigg et al., 2013; Lehmann et al., 2013; Zhang and Ferré-D'Amaré, 2013). Specifically, C56 of the T-box stacks on a nucleobase in the D-loop of the tRNA, and a point mutation of C56 to U has been shown to reduce the tRNA binding affinity by ~40 fold (Zhang and Ferré-D'Amaré, 2013). We introduced the same mutation in the T-box182 backbone (T-boxC56U) (Figure 6A, Supplementary file 1). The smFRET trajectories for tRNAGly binding to T-boxC56U are overall similar to the trajectories for WT T-box182, with a majority of traces (73 ± 6% of total traces) showing stable binding at 0.7 FRET state, and 13 ± 4% of the traces showing transitions back to the 0.4 FRET state (Figure 6C,D). τ0.7 was estimated to be at least ~23 s (limited by the photobleaching of the fluorophore) (Figure 6E). Post-synchronized time-evolved histogram on the subset of traces that demonstrated real-time binding shows fast transition to the fully bound state (Figure 6—figure supplement 1). Comparison of tRNAGly binding to T-boxC56U and T-box182 suggest that the C56U mutation does not affect the second binding step. To investigate whether the mutation at the T-loop region affects the first binding step, we analyzed the binding and dissociation of tRNAΔNCCA-Cy5 to T-boxC56U-Cy3(3’). We found that the k1 of tRNA binding to T-boxC56U was roughly 20-fold slower compared to tRNA binding to T-box182, and the dissociation was roughly 2.5-fold faster compared to T-box182 (Figure 6E), leading to a ~ 50 fold higher dissociation constant for the first binding step. Our results suggest that the T-loop region of the T-box is critical during the first binding step, potentially aiding in anticodon recognition, but does not contribute significantly to the second binding step.

Kinetic model for the two-step binding of glyQS T-box riboswitch and uncharged tRNAGly.

Details of the model are described in the text. Rate constants are summarized in Figure 6E.

https://doi.org/10.7554/eLife.39518.028
Figure 6 with 5 supplements see all
Regulation of the tRNAGly binding kinetics by structural elements in the glyQS T-box riboswitch.

(A) Schematic diagram of three different mutations introduced to the T-box182 backbone (T-boxC56U, T-boxSIII-Δ4bp and T-boxΔKT). (B) Representative smFRET traces of T-boxΔKT-Cy3(3’) and tRNAGly-Cy5. (C) FRET histograms of the T-box mutants with tRNAGly-Cy5. ‘N’ denotes the total number of traces in each histogram from three independent experiments. (D) TDP of the T-box mutants with tRNAGly-Cy5. TDPs are generated from all smFRET traces from three independent experiments and are plotted in the same way as in Figure 2D. (E) Table of kinetic parameters for tRNAGly-Cy5 binding to different T-box constructs. k-1, k2, and k-2 of T-boxΔKT-Cy3(3’) are apparent rate constants estimated to allow comparison as described in Materials and methods. All rate constants are reported as mean ± standard deviation (S.D.) from three or four independent experiments.

https://doi.org/10.7554/eLife.39518.029
Figure 6—source data 1

Cy3 and Cy5 intensity traces from three repeats for T-boxC56U-Cy3(3’) with tRNAGly-Cy5.

https://doi.org/10.7554/eLife.39518.035
Figure 6—source data 2

Cy3 and Cy5 intensity traces from four repeats for T-boxSIII-Δ4bp-Cy3(3’) with tRNAGly-Cy5.

https://doi.org/10.7554/eLife.39518.036
Figure 6—source data 3

Cy3 and Cy5 intensity traces from three repeats for T-boxΔKT-Cy3(3’) with tRNAGly-Cy5.

https://doi.org/10.7554/eLife.39518.037
Figure 6—source data 4

Cy3 and Cy5 intensity traces from three repeats for T-boxC56U-Cy3(3’) with tRNAΔNCCA-Cy5.

https://doi.org/10.7554/eLife.39518.038
Figure 6—source data 5

Cy3 and Cy5 intensity traces from three repeats for T-boxSIII-Δ4bp-Cy3(3’) with tRNAΔNCCA-Cy5.

https://doi.org/10.7554/eLife.39518.039
Figure 6—source data 6

Cy3 and Cy5 intensity traces from three repeats for T-boxΔKT-Cy3(3’) with tRNAΔNCCA-Cy5.

https://doi.org/10.7554/eLife.39518.040

A truncation of stem III has a minor effect on tRNA binding

The functional role of Stem III is unclear. It has been speculated that Stem III might serve as a transcription stalling site to allow co-transcriptional folding and regulation of the T-box riboswitch (Grundy and Henkin, 2004; Zhang and Landick, 2016). In addition, a SAXS data-derived model suggested coaxial stacking of Stem III and the anti-terminator stem, leading to a plausible role of Stem III in stabilizing the anti-terminator conformation in the presence of uncharged tRNAGly (Chetnani and Mondragón, 2017). To investigate the latter hypothesis, we generated a T-box mutant (T-boxSIII-Δ4bp), in which four base pairs in Stem III are deleted to significantly shorten its length (Figure 6A, Supplementary file 1). smFRET studies using T-boxSIII-Δ4bp-Cy3(3’) with tRNAΔNCCA-Cy5 and tRNAGly-Cy5 revealed insignificant difference in overall kinetics in the first and second step bindings (Figure 6C–E, Figure 6—figure supplement 2). Noticeably, the τ0.7 was around 50% shorter than that for the T-box182 (Figure 6E), indicating that Stem III may contribute to the stabilization of the fully bound state, potentially through coaxial stacking with the anti-terminator stem, but the effect is minor.

A K-turn mutation affects both binding steps

We next investigated the role of the K-turn in regulating tRNA binding kinetics. We disrupted the K-turn (T-boxΔKT) by changing the six bulged nucleotides to three nucleotides (UCA) to replace the K-turn with a three base pair stem (Figure 6A, Supplementary file 1). In contrast to binding of tRNAGly to T-box182, binding to T-boxΔKT results in three FRET states centered on ~0.2, 0.4, and 0.7. (Figure 6B). While the exact boundary of each FRET state is difficult to determine accurately from the FRET histogram (Figure 6C), a transition density plot (TDP) clearly reveals interconversion between the 0.2, 0.4, and 0.7 states (Figure 6D), with transitions between the 0.2 and 0.4 FRET states, and between the 0.4 and 0.7 FRET states more populated. Binding of tRNAΔNCCA-Cy5 to T-boxΔKT, on the other hand, leads to the loss of population of the 0.7 state; however, both the 0.2 and 0.4 FRET states and fluctuations between these two states are frequently sampled (Figure 6—figure supplement 3A,B).

Comparing the tRNAGly and tRNAΔNCCA binding, we speculate that both the 0.2 and 0.4 FRET states observed in the case of the T-boxΔKT represent the partially bound state in which only the anticodon and elbow are recognized. In contrast to T-box182, T-boxΔKT, with the K-turn replaced by an extension of Stem I, could potentially favor a relaxed conformation of Stem I, as observed in the NMR structure of an isolated K-turn and specifier loop domain (Wang and Nikonowicz, 2011), generating a lower FRET value centered at 0.2. However, the sampling of the 0.4 FRET state in the T-boxΔKT construct suggests that the interactions between the specifier and anticodon of the tRNA may transiently force open the extended base pair region and bend the T-box to adopt a similar conformation to the one observed in T-box182. The lifetimes of the 0.4 state before transition back to the 0.2 state and before transition forward to the 0.7 state are 0.30 ± 0.03 s and 0.13 ± 0.05 s, respectively (Figure 6—figure supplement 3C), indicating that this forced bent state is energetically unfavorable. However, this 0.4 FRET state is very likely to be required for the NCCA/t-box interaction to occur, as in the presence of tRNAGly the transitions from the 0.2 FRET to the 0.7 FRET state often pass through the 0.4 FRET state (Figure 6B). Furthermore, we observed a small region in the TDP corresponding to direct transitions between the 0.2 and 0.7 states. Given the very short lifetime of the 0.4 state, which is close to the time resolution of our experiments, the 0.2 to 0.7 state transition is likely to represent populations whose 0.4 FRET state lifetime is even shorter than the time resolution of the experiment.

The 3-step kinetic scheme for tRNA binding to T-boxΔKT is presented in Figure 6—figure supplement 3C. The association rate constant k1 ((5.4 ± 2.9) x104 M−1s−1), estimated from binding of tRNAΔNCCA to the T-boxΔKT, is ~9 fold smaller than binding to T-box182, suggesting that disruption of the K-turn affects anticodon recognition. Considering both the 0.2 and 0.4 FRET state as the partially bound state in T-boxΔKT, the apparent dissociation rate from the partially bound state k-1_app of tRNAΔNCCA to T-boxΔKT is 0.21 ± 0.06 s−1, similar to that for T-box182, suggesting that disruption of the K-turn does not affect the stability of the partially bound state. The 0.7 FRET state observed for T-boxΔKT in the presence of tRNAGly is consistent with the FRET value for the fully bound state in the WT T-box, indicating that the NCCA/t-box interactions in T-boxΔKT can still be formed. However, tRNAGly bound to T-boxΔKT0.7 = 1.6±0.3 s) is at least 15-fold less stable compared to tRNAGly bound to T-box1820.7 >24 s). Furthermore, in contrast to the signal observed for tRNAGly binding to T-box182, in which fewer than 10% of the FRET traces show transitions back to 0.4 FRET, in T-boxΔKT the vast majority of the traces show backward transitions to the 0.4 and 0.2 FRET states, contributing largely to the instability of the fully bound state. Based on the transition rates between 0.2, 0.4, and 0.7, we estimated the apparent forward (k2_app) and reverse (k-2_app) transition rates between the partially bound and the fully bound state of the T-boxΔKT to be 1.3 ± 0.3 s−1 and 1.0 ± 0.3 s−1, respectively (Figure 6E, Figure 6—figure supplement 3C). The dramatically reduced k2_app and increased k-2_app in T-boxΔKT leads to a ~ 150 fold change in the equilibrium constant in the second binding step compared to T-box182, implying that inflexibility of the K-turn region largely inhibits the conformational change in the T-box required to form the NCCA/t-box interaction, and strongly destabilizes the fully bound state.

Discussion

T-box riboswitches represent a unique class of riboswitches as they recognize a macromolecule and require interactions at multiple spatially separated sites on the ligand, unlike other riboswitches that respond to the binding of small ligands. Previous structural studies suggested that tRNA recognition by a T-box riboswitch is a bipartite process (Grigg and Ke, 2013; Zhang and Ferré-D'Amaré, 2013) with Stem I largely responsible for discriminating non-cognate tRNAs while the t-box sequence in the expression platform senses the charged state of the tRNA. We used smFRET to elucidate the binding kinetics of tRNAGly by the glyQS T-box riboswitch. With three FRET pairs between different T-box riboswitch and tRNA ligand constructs, our data collectively reveal a two-step model of uncharged tRNAGly binding to the glyQS T-box riboswitch (Figure 5). The first binding step involves recognition of the anticodon of the tRNA by the specifier sequence located in Stem I of the T-box riboswitch, leading to a partially bound state. In the second step, the 3’ end of the T-box docks into the NCCA end of the tRNA through interactions with the t-box sequence, which leads to a fully bound state. Without the NCCA interaction, the binding of tRNA is unstable, with a mean lifetime of ~4 s, whereas with interactions both with the anticodon and the NCCA end, the binding of tRNA is very stable, with a mean lifetime >24 s. It should be noted that the latter lifetime is likely to be much longer for two reasons: first, the measurement is limited by fluorophore photobleaching, i.e. the loss of signal is more likely to be due to photobleaching rather than the actual tRNA dissociation; second, the Cy5 label placed at the 5’ end of the tRNA reduces the overall binding affinity by ~4 fold (Figure 1—figure supplement 2), likely because the fluorophore impairs the NCCA/t-box interactions to some extent. In addition, an intra-T-box FRET pair at the 5’ and 3’ ends of the T-box demonstrates that, while the T-box is largely pre-organized in a folded state before tRNA binding, it still exhibits conformational rearrangement in a tRNA-dependent manner. Specifically, the 3’ half of the T-box (including Stem III and the anti-terminator) moves inward relative to the 5’ half (Stem I) of the T-box to accommodate the interaction with the NCCA end in the second binding step.

While our manuscript was in preparation, Suddala et al. (Suddala et al., 2018) reported a related single-molecule study on tRNA binding to the glyQS T-box riboswitch and proposed a similar two-step binding model. By using colocalized signals from the bound tRNA and the immobilized T-box, Suddala et al. (Suddala et al., 2018) uncovered two binding states distinguished by different dissociation rates of the tRNA, aided by using a Stem I-only mutant that cannot interact with the NCCA end of the tRNA. Specifically, in the model of Suddala et al. (Suddala et al., 2018) binding of the anticodon of the uncharged tRNA results in a relatively unstable state (with a lifetime of ~4–5 s), consistent with the partially bound state in our model, while binding both the anticodon and the NCCA end of the tRNA results in a stable state, consistent with the fully bound state in our model. Although both studies propose consistent kinetic models, in the study by Suddala et al. (Suddala et al., 2018), the FRET pair dyes are attached at the variable loop of the tRNA and the 3’ or 5’ end of the glyQS T-box, and do not generate different FRET signals to distinguish the partially bound state from the fully bound state. Therefore neither transitions between the partially bound and the fully bound states, nor the order of events during tRNA binding can be resolved in the study (Suddala et al., 2018). In contrast, by employing a FRET pair located at the 5’ end of the tRNA and the 3’ end of the glyQS T-box, we observed directly two FRET states corresponding to the recognition of the anticodon (0.4 FRET) and the binding of the NCCA (0.7 FRET), therefore our study allows the discrimination between different states and generates a more complete kinetic framework describing the full trajectory of the tRNA binding process. Our data reveal that anticodon recognition precedes the NCCA end interactions, and that after anticodon recognition the commitment to further establishment of the NCCA/t-box interaction is high. The T-box/tRNAGly complex transits rapidly from the partially bound state to the fully bound state, with a rate constant (k2) on the order of 10 s−1. Interestingly, our data also reveal that in the fully bound state the NCCA/t-box interaction is not highly-stable or ultra-stable. Brief disruption of the NCCA/t-box interactions can occur, but transition back to the fully bound state is rapid,~10 fold faster than dissociation of the tRNA from the partially bound state. Therefore tRNA can remain bound during the breaking and reforming of the NCCA/t-box interaction. However, such transient breaking of NCCA/t-box interaction was only observed in ~10% of the total population, suggesting that the reverse transition rate constant (k-2) is overall very small.

Two-step binding mechanisms have been observed in a variety of protein or nucleic acid mediated biological process, including T-cell receptor (TCR) recognition of the major histocompatibility complex (MHC) presenting peptides, where TCRs scan the MHC scaffold first, followed by sensing of specific MHC-presenting peptides (Wu et al., 2002); interaction of the signal recognition particle (SRP) receptor with the membrane, where SRP receptors interact with the membrane in a dynamic mode followed by an SRP-induced conformational transition into a stable binding mode (Hwang Fu et al., 2017), DNA interrogation by CRISPR Cas9-crRNA, where Cas9-crRNA recognizes the protospacer adjacent motif (PAM) on the target DNA followed by sensing of the spacer sequence and triggering R-loop formation (Sternberg et al., 2014); and RNA-induced silencing complexes (RISCs) binding to their mRNA targets, where dynamic sampling of the ‘sub-seed’ region occurs before targeting stably across the full seed region (Chandradoss et al., 2015; Herzog and Ameres, 2015; Salomon et al., 2015; Salomon et al., 2016). In all of these cases, a two-step binding mechanism provides a good balance between sensitivity and specificity. T-box riboswitches fold and function co-transcriptionally in the cell, and hence we propose that a two-step kinetic model is kinetically beneficial during co-transcriptional folding of the T-box and sensing of the tRNA ligand. First, considering that the intracellular concentration of tRNA is on the order of µM in bacteria (Avcilar-Kucukgoze et al., 2016; Dong et al., 1996; Emilsson and Kurland, 1990), the binding of the tRNA ligand to the T-box riboswitch is rate limited by the first step. Pausing at Stem III in the presence of NusG in vitro was estimated to last 30–60 s (Grundy and Henkin, 2004). Assuming an in vivo concentration of tRNA ~1 μM, the apparent tRNA binding rate in the first step is ~0.5 s−1, which is sufficient for the uncharged tRNA to reach the partially bound state during the pausing. Second, bacterial RNA polymerase (RNAP) elongates at a rate of 40–50 nt/s (Mosteller and Yanofsky, 1970; Vogel and Jensen, 1994), and the length between terminator and anti-terminator sequences is ~40 nts; therefore the decision between termination and anti-termination needs to be made within ~0.5–1 s. On the one hand, the relatively slow dissociation (k-1) from the partially bound state can ensure the tRNA ligand stays bound until the completion of the anti-terminator sequence. On the other hand, very rapid transition into the fully bound state (k2) helps secure the interaction between the NCCA end of the tRNA and the t-box sequence and trap the T-box in the anti-terminator conformation before completion of transcription of the terminator sequence, which is immediately downstream. Finally, the fully bound state has a lifetime of at least 24 s, allowing RNAP to elongate more than 1000 nts before the tRNA dissociates. Therefore, the second binding step in the WT T-box with uncharged tRNAGly is close to irreversible in such biological setting. It is worth mentioning that Suddala et al. (Suddala et al., 2018) revealed that the uncharged T-box/tRNAGly complex exists in two populations: a stable complex, with a tRNA dissociation rate of ~0.03 s−1, and an ultra-stable complex, with a tRNA dissociation rate <2.4×10−4 s−1. Our estimation of the lifetime of the fully bound state was limited by the introduction of Cy5 to the 5’ end of the tRNA and fluorophore photobleaching, and therefore our experiments cannot distinguish between the stable and ultra-stable complexes. Nevertheless, both potential configurations in the fully bound state are stable enough to allow transcription read-through.

Based on the known structures of Stem I in complex with tRNAGly (Grigg and Ke, 2013; Zhang and Ferré-D'Amaré, 2013), it was hypothesized that the intra-T-box conformational change is likely to involve the K-turn region, which sits at the junction between the 5’ and 3’ portions of the T-box. Our observation that the NCCA end of the uncharged tRNA moves towards the anti-terminator stem (revealed by the T-box182-Cy3(3’) and tRNAGly-Cy5 FRET), but maintains its relative position to the K-turn region (revealed by the T-box182-Cy3(5’) and tRNAGly-Cy5 FRET) during the second binding step is consistent with the role of the K-turn region as the hinge of the intra-T-box conformational change. In addition, such conformational changes mediated by the K-turn region are Mg2+ dependent, as suggested by Suddala et al. (Suddala et al., 2018). Importantly, our results hint at the critical role of the K-turn region in promoting the fast coordination between anticodon sensing and locking of the NCCA end. Crystal structures of the Stem I/tRNA complex (Grigg and Ke, 2013; Zhang and Ferré-D'Amaré, 2013) show that Stem I flexes around the K-turn region and that this flexing seems to be important to establish the interactions between the anticodon and specifier sequence. Using a mutant where the K-turn region is removed, our data show that in the absence of the K-turn motif both the first and second binding steps are affected, but with a much more dramatic effect on the second step. In this specific K-turn mutant, where the K-turn is replaced by an extended Stem I helix, the conformational change in the T-box that brings closer the 3’ and 5’ portions becomes highly energetically unfavorable, leading to a ~ 10 fold reduction in k2, and ~15 fold destabilization of the fully bound state. Our kinetic measurement of the K-turn mutant explains the in vivo loss-of-function K-turn mutants (Winkler et al., 2001), and emphasizes the importance of the flexing around the K-turn region and the associated conformational changes in the T-box itself in the overall recognition and binding process.

Our results with various T-box mutants demonstrate that the T-box structural elements involved in tRNA recognition can drive the two-step binding process in a modular fashion. For example, in the first binding step, interactions with Stem I involve both anticodon recognition and the contacts between the interdigitated T-loops in Stem I and the elbow region (D- and T-loops) of the tRNA (Grigg and Ke, 2013; Zhang and Ferré-D'Amaré, 2013). When introducing a point mutation that impairs the T-loops/elbow region interaction (Zhang and Ferré-D'Amaré, 2013), we observed a dramatic decrease in the association rate constant and a moderate increase in the dissociation rate constant, leading to an overall ~50 fold reduction on the binding affinity for the first step, however the second step is unaffected. This observation suggests that establishment of the T-loops/elbow interactions is an important part of the Stem I/tRNA recognition process in the first binding step, but does not play any role in NCCA recognition in the second binding step. Similarly, truncation of Stem III has a minor effect on the stability of the fully bound state, but with no influence on the first binding step. Finally, the K-turn region (discussed above), which links the 5’ and 3’ portions of the T-box, plays a key role in coordinating the two binding steps by providing structural flexibility. It is worth mentioning that the impact on the tRNA binding kinetics by mutations in the T-loops, the K-turn, and the Stem III are consistent with the sequence conservation and the in vivo impact on amino acid-mediated transcription read-through using tyrS T-box riboswitches as a model system (Henkin, 2014; Rollins et al., 1997; Winkler et al., 2001). Mutations in the highly conserved T-loop and K-turn region have a more dramatic influence on the tRNA binding kinetics, translating into a larger in vivo impact. While the less conserved Stem III contributes to the stabilization of the anti-terminator conformation in vitro, deletion of it does not appear to significantly impair the tRNA binding process in vitro. Potentially its major function is to create a pause site to coordinate with the co-transcriptional folding of the T-box (Grundy and Henkin, 2004; Zhang and Landick, 2016).

In conclusion, our study provides a comprehensive kinetic framework for describing tRNA recognition by the T-box riboswitch. The two-step binding process is driven by the specific structural elements of the T-box, and is likely to be kinetically beneficial for efficient, co-transcriptional recognition of the cognate tRNA ligand. Specific T-box structural elements drive the two-step binding process in a modular fashion, that is the 5’ and 3’ portions of the T-box are responsible for the first and second binding steps, respectively, with the K-turn region coordinating the two binding steps by allowing structural flexibility, providing a guideline for synthetic biology design of RNA regulatory modules, as well as for the development of new antibiotics using critical T-box structural elements as potential targets. Finally, the glyQS T-box riboswitch represents one of the simplest members of this class of riboswitches. Other T-box riboswitches are larger and have additional structural elements and can even appear in tandem arrangements (Gutiérrez-Preciado et al., 2009). While the two-step binding kinetics may be common to all T-box riboswitches, it is likely that the process in other T-box riboswitches shows differences modulated by the additional structural elements.

Materials and methods

Key resources table
Reagent type (species)
or resource
DesignationSource or referenceIdentifiersAdditional information
Gene (Bacillus subtilis)glyQS T-box riboswitchPMID: 28531275
Gene (Bacillus subtilis)tRNA-GlyPMID: 9023104RRID:SCR_008636
Strain, strain
background (E. coli)
DH5αThermo Fisher ScientificCatalogue # 18265017
Recombinant
DNA reagent
pUC19 VectorNew England BiolabsCatalogue # N3041S
Peptide,
recombinant protein
T7 RNA polymerasePMID: 3684574
Peptide,
recombinant protein
BsaINew England BiolabsCatalgue #: M0535S
Peptide,
recombinant protein
RppHNew England BiolabsCatalgue #: M0356S
Chemical
compound, drug
Cyanine5 NHS esterLumiprobeCatalogue # 13020
Chemical
compound, drug
Cyanine3 hydrazideLumiprobeCatalogue # 13070
Software, algorithmOrigin 5.0Microcal OriginRRID:SCR_002815
Software, algorithmImageJImageJRRID:SCR_003070
Software, algorithmMATLABMATLABRRID:SCR_001622
Software, algorithmNIS-ElementsNikonRRID:SCR_014329

RNA purification and mutagenesis

Request a detailed protocol

RNA transcription was performed in vitro using His6-tagged T7 RNA polymerase using standard protocols (Milligan et al., 1987). Cloning, design of a bicistronic DNA template encoding the B. subtilis glyQS T-box riboswitch and its cognate tRNAGly, and conditions for in vitro transcription were described before (Chetnani and Mondragón, 2017). For the experiments, all the RNAs from the crude transcription reaction were purified on a 7.5% denaturing (8 M Urea) polyacrylamide gel. The RNAs of interest were located on the gel by UV shadowing, the bands were cut out, and the RNAs were eluted into 50 mM sodium acetate (pH 7.0) buffer containing 200 mM potassium chloride by overnight rocking at 4°C. The eluted RNAs were precipitated by adding 3 volumes of cold 100% ethanol and stored overnight at −20°C. The precipitated RNAs were pelleted by centrifugation for 15 min at 20,000 g. The RNA pellets were washed three times in cold 80% ethanol, dried in a Speedvac and re-suspended in water. The concentration was estimated by its absorbance at 260 nm and was kept frozen at −20°C for long term storage. All mutant constructs were made using a commercial site directed mutagenesis protocol (Kunkel, 1985) (Quikchange, Stratagene) and the RNAs were produced and purified by using the same protocol described above. The sequence of all mutant T-box and tRNA constructs were confirmed by sequencing and validated by MFOLD (Zuker, 2003) to ensure that the secondary structural elements were not affected by the mutations.

Fluorescent labelling of the RNA

Request a detailed protocol

For smFRET experiments, end labelling of RNA molecules was performed by modifying standard labelling protocols (Rinaldi et al., 2015). For 3’ end labelling, 50 µg of RNA in 50 µL of reaction volume was incubated with 0.1 M Na-periodate in 0.1 M Na-acetate buffer at pH 5.2 for 90 min in the dark. The reaction was quenched by adding 5 µL of 2.5 M KCl and incubating on ice for 10 min. The resultant insoluble KIO4 was removed by centrifugation at 20,000 g for 30 min and the supernatant was passed through a P6 column (Bio-Rad) to exchange the buffer to 0.1 M HEPES (pH 7.0), 40% DMSO. The RNA was incubated with Cy3 hydrazide (Lumiprobe) dye for 45 min with a final RNA/dye ratio of ~1:200. The RNA was then ethanol precipitated as described above. The precipitated RNA was pelleted by centrifugation for 15 min at 20,000 g, dried in a Speedvac and re-suspended in water. The final RNA solution was passed through a P6 column to remove any residual free dye.

5’ end labelling of RNA was performed by N-(3-Dimethylaminopropyl)-N′-ethylcarbodiimide hydrochloride (EDC) – N-hydroxysuccinimide (NHS) coupling through activation of the 5’ monophosphate of the RNA by EDC and imidazole (Rinaldi et al., 2015). To improve the overall labelling efficiency of this method, a modified approach was used in which the 5’ triphosphate of 100 µg of RNA was converted to 5’ monophosphate in 100 µL reaction volume by incubating it with 100 units of RNA 5' Pyrophosphohydrolase (NEB) at 37°C for 1 hr. The enzyme was removed by phenol chloroform extraction and the supernatant was passed through a P6 column to exchange the buffer to 10 mM HEPES (pH 7.0), 150 mM NaCl, 10 mM EDTA. This was followed by addition of 12.5 mg of EDC to the RNA solution along with 50 µL of ethylene diamine and 200 µL of 0.1 M imidazole buffer (pH 6.0). The reaction was incubated for 3 hr at 37°C and the RNA was then ethanol precipitated as described above. The resultant RNA pellet was re-suspended in 0.1 M sodium carbonate buffer (pH 8.7) and residual EDC was removed by passing the solution through a P6 column. The resultant RNA solution in 0.1 M sodium carbonate buffer (pH 8.7) was incubated with Cy5 NHS (Lumiprobe) dye for 45 min with a final RNA/dye ratio of ~1:200. The RNA was ethanol precipitated as described above and re-suspended in water. The final RNA solution was passed through a P6 column to remove any residual free dye. The overall labelling efficiency of RNA constructs used in this study varied from 75% to 95%.

Isothermal titration calorimetry

Request a detailed protocol

The thermodynamic parameters associated with binding of tRNAGly to glyQS T-box were determined at 25°C by using an ITC-200 Micro-Calorimeter (MicroCal). Prior to the experiment, the two interacting RNAs (in 50 mM HEPES pH 7.0, 100 mM KCl) were refolded separately by first heating for 3 min at 90o C followed by incubation on ice for 2 min. At this point, MgCl2 was added to attain a final concentration of either 1 mM or 10 mM. The RNA solution was then heated to 50o C for 10 min and 37o C for 30 min followed by cooling to 25°C. The refolded RNA was then concentrated using a 10 kDa cutoff Amicon filter and washed three times with the final ITC buffer (50 mM HEPES pH 7.0, 100 mM KCl and 10 mM MgCl2 or 1 mM MgCl2). For testing the binding of T-box with unlabeled and 5’ Cy5 labelled tRNA in 10 mM MgCl2 buffer, the sample cell was filled with 6.6 µM and 11 µM T-box respectively and the corresponding concentration of tRNAGly in the syringe was 111 µM and 152 µM. For the ITC experiment in 1 mM MgCl2 buffer, the cell was filled with T-box at a concentration of 6.7 µM and the syringe concentration of unlabeled tRNAGly was 72 µM. For each ITC experiment, the titration was carried out by stepwise (2 µL) injection of tRNAGly from a stirred syringe (1000 rev/min) into the sample cell. Successive injections were spaced by 150 s and values for the change in enthalpy (ΔHb), association constant (Kb), change in enthalpy (ΔHb) and stoichiometry (n) were determined by nonlinear least-squares fitting of the data using Origin 5.0 software (OriginLab).

Electrophoretic mobility shift assay

Request a detailed protocol

The T-box and tRNA were denatured and refolded using the same protocol as described for the ITC and smFRET experiments. For binding of tRNA, 5 µL of 4 µM of folded T-box and 5 µL of 2 µM folded tRNA samples were mixed together for 30 min in final buffer containing 50 mM HEPES pH 7.0, 100 mM KCl and 15 mM MgCl2. The folded T-box or T-box +tRNA mixture were loaded on to 6% native polyacrylamide gel containing 15 mM MgCl2 and the gel was run at room temperature for 2–2.5 hr. RNAs were stained by SYBRTM green RNA staining dye (Invitrogen). The gel was imaged using ChemiDoc (Bio-rad) in SYBR green channel (for unlabeled T-box constructs), Cy5 channel (for tRNA-Cy5), or Cy3 channel (for labeled T-box constructs). Images were processed and analyzed by ImageJ (Schneider et al., 2012).

Fluorophore conjugation of DNA oligos

Request a detailed protocol

DNA oligos that hybridize to 5’ extension and 3’ extension of the T-box construct were purchased from Integrated DNA Technologies with an amine modification at the 5’ end and 3’ end respectively. 13.5 µL of 100 µM DNA oligo was mixed with 1.5 µL of 1 M NaHCO3 (pH 8.6). 25 µg of NHS conjugated fluorophore (Cy3 or Cy5) was dissolved with 0.5 µL DMSO and mixed with the DNA oligo solution. The mixture was incubated at 37°C overnight. 1.67 µL of 3 M NaOAc and 50 µL of pure ethanol was added to the mixture to precipitate the conjugated DNA oligo overnight at −20°C. The precipitated DNA oligo was pelleted by centrifugation for 30 min at 21000 g and re-suspended with 40 µL water. The DNA solution was passed through a P6 column to remove any residue free dye and salt. The overall fluorophore labeling efficiency is ~60%.

smFRET measurements

Request a detailed protocol

Slides containing microfluidic channels were prepared as previously described (Blanchard et al., 2004). Slides and coverslips were coated with a mixture of poly-ethylene glycol (PEG, Mw = 500,000) and PEG-biotin (Mw = 500,000) according to previously published protocol (Blanchard et al., 2004). The T-box and tRNA were denatured and refolded using the same protocol described for the ITC experiments with the final buffer containing 50 mM HEPES pH 7.0, 100 mM KCl, and 15 mM MgCl2. T-box RNA was hybridized to the biotinylated DNA oligo at the extension during refolding, and immobilized via biotin-streptavidin interactions to the surface. tRNA was diluted in imaging buffer (50 mM HEPES pH 7.0, 100 mM KCl, 15 mM MgCl2, 5 mM protocatechuic acid (PCA) (Sigma), 160 nM protocatechuate-3,4-dioxygenase (PCD) (Sigma), and 2 mM Trolox (Sigma)) and flowed into the microfluidic channels. smFRET measurements were performed with an objective based total internal reflection fluorescence (TIRF) microscope based on a Nikon Ti-E with 100X NA 1.49 CFI HP TIRF objective (Nikon). A 561 nm laser (Coherent Obis at a power density of 4.07 × 105 W/cm2) was used for the FRET measurement. A 647 nm laser (Cobolt MLD at a power density of 5.88 × 105 W/cm2) was used for direct excitation of Cy5 to check the presence of the acceptor. Emissions from both donor and acceptor were passed through an emission splitter (OptoSplit III, Cairn), and collected at different locations on an EMCCD (iXon Ultra 888, Andor). 1500 frames of time-lapse images were taken with 100 ms exposure time. Each independent measurement in our study should be considered as a technical replicate. The biological samples (i.e. the in vitro transcribed tRNAs and T-boxes) were generated once. In each measurement, T-box and tRNA were folded, smFRET images were recorded, and analysis were performed independently.

Lifetime analysis

Request a detailed protocol

Individual spots were picked from maximum intensity projection of Cy5 emission channel (Figure 2—figure supplement 1) using NIS Elements software, that is only the pixels with the highest intensity value of the same XY coordinates in time-lapse images are displayed in the maximum intensity projection image. Fluorescent intensity trajectories of these spots were generated from Cy3 and Cy5 channels of time-lapse images, and corrected for baseline and bleed-through in MATLAB as previously described (Fei et al., 2008). FRET traces were generated by calculating ICy5 / (ICy5+ICy3) at each time point from the intensity trajectories. smFRET traces were idealized by fitting with a Hidden Markov Model using vbFRET (Bronson et al., 2009). Dwell time of each FRET state before transition to another FRET state of individual traces was extracted from the idealized traces and the dwell time histograms were fit with exponential decay by Origin 7.0 (OriginLab) (Figure 2—figure supplement 2). In most of the cases, dwell time histograms fit well with a single exponential decay (A.exp(-t/t0) + y0). In a few cases, in which mixed populations with different lifetimes were expected, dwell time histograms were fit with a double exponential decay (A1.exp(-t/t1) + A2.exp(-t/t2) + y0) and the population-weighted average lifetime was calculated by (A1.t1 +A2.t2) / (A1 + A2). The latter cases are explicitly mentioned in the text.

Determination of kinetic parameters

(I) k-1

Request a detailed protocol

For all T-box constructs, k-1 was estimated from the dwell time of the partially bound state (0.4 FRET state of T-box182, T-boxC56U and T-boxSIII-Δ4bp) of tRNAΔNCCA-Cy5 (Figure 2—figure supplement 2). The dwell time of the partially bound state was fit with a single exponential decay. The lifetime of 0.4 FRET state (τ0.4-0) was extracted and k-1 was calculated using the following equation:

k-1=1τ0.4-0

(II) k1

Request a detailed protocol

We noticed that the lifetime of the unbound state (τ0-0.4) of tRNAΔNCCA-Cy5 with Tbox182-Cy3(3’) was 31.3 ± 5.3 s, while the average smFRET trace length (τtrace) was 70 ± 7 s under 4.07 × 105 W/cm2 561 nm laser (Figure 2—figure supplement 2), limited by Cy3 photobleaching. Therefore, there was a possibility of underestimating the lifetime of the unbound state, or overestimating k1. We hence characterized the photobleaching effect on unbound state lifetime analysis using reduced 561 nm laser power (1.96 × 105 W/cm2). τ0-0.4 was calculated to be 46 ± 4 s and τtrace was 103.0 ± 0.6 s. To correct for the photobleaching effect, we estimated k1, using the following equation (Bartley et al., 2003; Fei et al., 2009):

k1=1τ0-0.4-1τtrace×1[tRNA]

where [tRNA] represents the concentration of tRNAΔNCCA-Cy5. k1 values were consistent using both high laser power ((5.8 ± 0.2) x 105 M−1s−1) and low laser power ((4.1 ± 0.6) x 105 M−1s−1). In Figure 6E, we reported the mean ±S.D. of k1 obtained with both high and low laser powers.

We also measured τoff in the flow experiment and it was determined by the dwell time between the injection time and detection of the first binding event (Figure 2—figure supplement 4C). k1 calculated for tRNAGly in the flow experiment was consistent with the value for tRNAΔNCCA in the steady-state experiment. We therefore used k1 of tRNAΔNCCA-Cy5 binding to the T-box to approximate the rate constant for the first binding step in all T-box constructs.

For some T-box mutants (specifically, T-boxC56U-Cy3 and T-boxΔKT-Cy3), binding of tRNA-Cy5 became less efficient compared to the T-box182-Cy3(3’), already a hint of a deficient first binding step. In this case, estimation of k1 from the selected FRET traces will highly bias for the molecules that bind relatively fast within the photobleaching time of Cy3, therefore significantly overestimating k1. We therefore applied a different approach to estimate k1 for the T-box mutants. We counted trace percentage of the T-box mutants (number of FRET traces divided by the number of Cy3 spots) and normalized it to the trace percentage of T-box182 with tRNAΔNCCA-Cy5 to get the normalized trace percentage (f) (Figure 6—figure supplement 4).

We performed Gillespie simulations (Gillespie, 1976; Gillespie, 1977; Gillespie, 2007) on the observed signal as a function of k1 and k-1 during the first binding step. The simulated signal mimicked the experiment and imaging processing. As our FRET traces were selected from the maximum intensity projection of the acceptor signal (Figure 2—figure supplement 1) regardless of the lifetime of the bound state, as long as there was a binding event during the imaging time window (70 s considering the photobleaching time of the donor dye), it would generate a signal (as a ‘footprint’) that showed in the output image of maximum projection (Figure 6—figure supplement 5A). For each parameter set, 2000 binding trajectories were simulated, and the fraction of the molecules that generated signal in the maximum intensity projection (referred to as ‘f’) within the 70 s time window was calculated. The simulation shows that f is insensitive to k-1, at least within the k-1 we observed for the various T-box constructs (Figure 6—figure supplement 5B), and that the dependence of f on k1 is best described as

f=f0-aebk1

where f0 is the binding fraction at saturation, and a and b are constants (Figure 6—figure supplement 5B). By fitting ln(f0-f) vs. k1 with a linear function (Figure 6—figure supplement 5C) we estimated k1 as

k1=0.45-ln1.557-f72.09

The standard deviations reported for k1 were calculated by an error propagation function considering the standard deviation of f.

(III) k2 and k-2

Request a detailed protocol

Even though we observed that smFRET traces of tRNAGly with T-box182, T-boxC56U and T-boxSIII-Δ4bp exhibit real-time binding, the majority of the traces show single-step binding from unbound to fully bound state while only a few traces sample the partially unbound state within 0.1 s (Figure 2F, Figure 2—figure supplement 4B,Figure 6—figure supplements 1A and 2A). For this reason, our calculation of k2 is limited by the imaging time resolution and we estimate the transition from 0.4 to 0.7 is rapid with a k2 larger than 10 s−1.

We can also detect transition from 0.7 FRET to 0.4 FRET, and calculate the lifetime of 0.7 FRET to 0.4 FRET (τ0.7-0.4) (Figure 2—figure supplement 2). However, the fraction of traces showing this transition is small. Since most of the traces demonstrate long time 0.7 FRET, we reason that k-2 should be very small and the second binding step for these three constructs is close to irreversible.

(IV) Kinetic model for T-boxΔKT

Request a detailed protocol

The kinetic model for tRNA binding to the K-turn mutant T-boxΔKT is illustrated in Figure 6—figure supplement 3C, with kinetic parameters marked. Specifically, the partially bound state contains a relaxed conformation and a bent conformation at the original K-turn region. The interconversion rates between 0.2, 0.4 and 0.7 FRET states were determined by the lifetime of the each FRET state before transitioning to the other FRET state. To compare better the kinetic parameters with other T-box constructs, we calculated the apparent k2 considering both pathways (0.2→0.7 and 0.2→0.4→0.7) (Figure 6—figure supplement 3C) as below:

k2_app=p1k1'k2'(k1'+k-1')+p2k2

Where p1 and p2 are the probabilities of the 0.2→0.4→0.7 and 0.2→0.7 pathways respectively. And we calculated the apparent k-2 considering both pathways (0.7→0.2 and 0.7→0.4) (Figure 6—figure supplement 3C) as shown below:

k-2_app=p1k-2'+p2k-2

Where p1 and p2 are the probabilities of the 0.7→0.4 and 0.7→0.2 pathways, respectively. To compare with the k-1 of other T-box constructs, the apparent k-1_app of T-boxΔKT was estimated by considering both 0.4 and 0.2 FRET states.

Data availability

All data generated or analysed during this study are included in the manuscript and supporting files. smFRET trajectories of each data sets included in the manuscript are available in Source Data.

References

Article and author information

Author details

  1. Jiacheng Zhang

    Institute for Biophysical Dynamics, University of Chicago, Chicago, United States
    Contribution
    Designed the research, Performed single molecule FRET experiments, Data analysis, Wrote and revised the manuscript
    Contributed equally with
    Bhaskar Chetnani
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8671-5888
  2. Bhaskar Chetnani

    Department of Molecular Biosciences, Northwestern University, Evanston, United States
    Present address
    Innovation and Technology Transfer Division, Rush University Medical Center, Chicago, United States
    Contribution
    Designed the research, Prepared samples and performed ITC experiments, Wrote and revised the manuscript
    Contributed equally with
    Jiacheng Zhang
    Competing interests
    No competing interests declared
  3. Eric D Cormack

    University of Chicago College, Chicago, United States
    Contribution
    Helped with smFRET imaging and data analysis
    Competing interests
    No competing interests declared
  4. Dulce Alonso

    Department of Molecular Biosciences, Northwestern University, Evanston, United States
    Contribution
    Sample preparation and ITC experiments
    Competing interests
    No competing interests declared
  5. Wei Liu

    Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, United States
    Contribution
    Helped with data analysis
    Competing interests
    No competing interests declared
  6. Alfonso Mondragón

    Department of Molecular Biosciences, Northwestern University, Evanston, United States
    Contribution
    Designed the research, Wrote and revised the manuscript
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0423-6323
  7. Jingyi Fei

    1. Institute for Biophysical Dynamics, University of Chicago, Chicago, United States
    2. Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, United States
    Contribution
    Designed the research, Helped with data analysis, Wrote and revised the manuscript
    For correspondence
    jingyifei@uchicago.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9775-3820

Funding

Chicago Community Trust

  • Alfonso Mondragón
  • Jingyi Fei

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was funded by a grant to JF and AM by the Chicago Biomedical Consortium with support from the Searle Funds at the Chicago Community Trust. We thank E Heideman for the preparation of slides for all imaging experiments and managing the Fei laboratory.

Version history

  1. Received: June 27, 2018
  2. Accepted: September 24, 2018
  3. Accepted Manuscript published: September 25, 2018 (version 1)
  4. Version of Record published: October 22, 2018 (version 2)

Copyright

© 2018, Zhang 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

  • 2,872
    views
  • 381
    downloads
  • 24
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Jiacheng Zhang
  2. Bhaskar Chetnani
  3. Eric D Cormack
  4. Dulce Alonso
  5. Wei Liu
  6. Alfonso Mondragón
  7. Jingyi Fei
(2018)
Specific structural elements of the T-box riboswitch drive the two-step binding of the tRNA ligand
eLife 7:e39518.
https://doi.org/10.7554/eLife.39518

Share this article

https://doi.org/10.7554/eLife.39518

Further reading

    1. Cell Biology
    2. Structural Biology and Molecular Biophysics
    Aaron JO Lewis, Frank Zhong ... Ramanujan S Hegde
    Research Article

    The protein translocon at the endoplasmic reticulum comprises the Sec61 translocation channel and numerous accessory factors that collectively facilitate the biogenesis of secretory and membrane proteins. Here, we leveraged recent advances in cryo-electron microscopy (cryo-EM) and structure prediction to derive insights into several novel configurations of the ribosome-translocon complex. We show how a transmembrane domain (TMD) in a looped configuration passes through the Sec61 lateral gate during membrane insertion; how a nascent chain can bind and constrain the conformation of ribosomal protein uL22; and how the translocon-associated protein (TRAP) complex can adjust its position during different stages of protein biogenesis. Most unexpectedly, we find that a large proportion of translocon complexes contains RAMP4 intercalated into Sec61’s lateral gate, widening Sec61’s central pore and contributing to its hydrophilic interior. These structures lead to mechanistic hypotheses for translocon function and highlight a remarkably plastic machinery whose conformations and composition adjust dynamically to its diverse range of substrates.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Roberto Efraín Díaz, Andrew K Ecker ... James S Fraser
    Research Article

    Chitin is an abundant biopolymer and pathogen-associated molecular pattern that stimulates a host innate immune response. Mammals express chitin-binding and chitin-degrading proteins to remove chitin from the body. One of these proteins, Acidic Mammalian Chitinase (AMCase), is an enzyme known for its ability to function under acidic conditions in the stomach but is also active in tissues with more neutral pHs, such as the lung. Here, we used a combination of biochemical, structural, and computational modeling approaches to examine how the mouse homolog (mAMCase) can act in both acidic and neutral environments. We measured kinetic properties of mAMCase activity across a broad pH range, quantifying its unusual dual activity optima at pH 2 and 7. We also solved high-resolution crystal structures of mAMCase in complex with oligomeric GlcNAcn, the building block of chitin, where we identified extensive conformational ligand heterogeneity. Leveraging these data, we conducted molecular dynamics simulations that suggest how a key catalytic residue could be protonated via distinct mechanisms in each of the two environmental pH ranges. These results integrate structural, biochemical, and computational approaches to deliver a more complete understanding of the catalytic mechanism governing mAMCase activity at different pH. Engineering proteins with tunable pH optima may provide new opportunities to develop improved enzyme variants, including AMCase, for therapeutic purposes in chitin degradation.