Reversible-Switch Mechanism of the SAM-III ... - ACS Publications

Nov 8, 2016 - The Enterococcus faecalis SAM-III (SMK) riboswitch regulates downstream gene expression through conformational change by sensing ...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/JPCB

Reversible-Switch Mechanism of the SAM-III Riboswitch Sha Gong,†,‡ Yujie Wang,† Zhen Wang,† Yanli Wang,† and Wenbing Zhang*,† †

Department of Physics, Wuhan University, Wuhan, Hubei 430072, P. R. China College of Mathematics and Physics, Huanggang Normal University, Huanggang, Hubei 438000, P. R. China



ABSTRACT: Riboswitches are self-regulatory elements located at the 5′ untranslated region of certain mRNAs. The Enterococcus faecalis SAM-III (SMK) riboswitch regulates downstream gene expression through conformational change by sensing S-adenosylmethionine (SAM) at the translation level. Using the recently developed systematic helix-based computational method, we studied the co-transcriptional folding behavior of the SMK riboswitch and its shortened construct lacking the first six nucleotides. We find that there are no obvious misfolded structures formed during the transcription and refolding processes for this riboswitch. The full-length riboswitch quickly folds into the ON-state in the absence of SAM, and the coupling between transcription and translation is not required for the riboswitch to function. The potential to form helix P0 is necessary for the riboswitch to function as a switch. For this thermodynamically controlled reversible riboswitch, the fast helix-exchanging transition pathway between the two functional structures guaranteed that this riboswitch can act as a reversible riboswitch.

1. INTRODUCTION Riboswitches, as self-regulatory elements, are involved in regulating downstream gene expression, a function that has been traditionally attributed to proteins.1−5 Most riboswitches consist of two distinct domains: a conserved aptamer and a variable expression platform. 6−10 The aptamer, as an independent folding unit, is responsible for recognizing and binding the target metabolite, and the expression platform is involved in the genetic regulation.11−15 Riboswitches control gene expression through conformational changes triggered by ligand binding at the levels of translation, transcription, and splicing in cells.11,12,16−19 The Enterococcus faecalis SAM-III (SMK) riboswitch responds to the intracellular concentration of SAM or its derivative S-adenosylhomocysteine SAH to regulate the translation of the metK gene,20,21 which encodes the synthetase of SAM, a cofactor in methylation reactions of proteins, nucleic acids, and other biomolecules.22,23 The crystal structure of the SMK riboswitch with its ligand binding has been reported recently.24 This riboswitch was found to be a reversible switch in response to the ligand concentration in cells.25 Unlike most riboswitches, this riboswitch only has a simple architecture that utilizes a single domain for both ligand binding and gene regulation.21,26 Its Shine−Dalgarno (SD) sequence directly takes part in binding SAM and is sequestered by base pairing with the anti-SD (ASD) sequence in the presence of the ligand (Figure 1). Helices P1, P2, P3, and P4 form the binding pocket for SAM, and ligand binding stabilizes the OFF-state by introducing tertiary interactions. The sequestration of the SD sequence in the OFF-state inhibits initiation of the translation process and represses downstream gene expression.25,26 In the absence of SAM, the riboswitch © 2016 American Chemical Society

adopts an alternative folding pattern, at which the SD sequence is unstructured and free to bind the ribosomal subunit.27,28 Thus, the riboswitch acts as a genetic ON switch without the binding of the ligand. Previous studies have investigated many features of the SMK riboswitch such as ligand specificity and its structural basis,21 structure transition,2,28 and rates of ligand binding and dissociation.25,27 However, these studies have mainly focused on the full-length transcriptional products. In cells, folding of riboswitches, which are synthesized during the transcription process, follows a sequential progression.29,30 RNA folding pathways and kinetic traps during the transcription have been suggested to be different from that of the renaturation folding, and the transcription process could dictate RNA folding pathways and kinetic traps.31−33 For most riboswitches, their biological functions were found to be tied to the cotranscriptional folding.11,34−37 However, the co-transcriptional folding behaviors and the reversible switching mechanism of the SMK riboswitch remain unclear, for example, if transcription and translation are coupled for the riboswitch to function, then what is the transition pathway between the two function structures, and is the transition fast enough for the riboswitch to play the reversible switch function? To fully investigate regulation mechanisms of the SMK riboswitch, in this study we used the recently developed helix-based kinetic folding method to predict the co-transcriptional folding and refolding behavior of the riboswitch as well as its shortened construct without the Received: September 25, 2016 Revised: November 5, 2016 Published: November 8, 2016 12305

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B

conformation space is calculated using the nearest-neighbor model.38,39 The population kinetics, pi(t), for each state i at time t is described by the master equation dpi(t)/dt = ∑j[kj→i pj(t) − ki→j pi(t)], where ∑j denotes the sum over all the conformations, and ki→j and kj→i are the rate constants for the respective transitions. The above equation has a matrix form: dp(t)/dt = M·p(t), where p(t) denotes the vector of the population distribution, and M is the rate matrix with elements Mij = ki→j(i ≠ j) and Mii = −∑i≠j ki→j. For t > 0, the master equation yields the population kinetics p(t ) =

∑ Cmn m e−λ t m

m=1

(1)

where nm and −λm are the mth eigenvector and eigenvalue of rate matrix M, respectively, and coefficient Cm is determined by the initial conditions. The rate constants for formation (k+) and disruption (k−) of a stack are written as k+ = k0 e−TΔSstack/kBT and k− = k0 e−ΔHstack/kBT, where k0 is equal to 6.6 × 1012 s−1 for the formation/disruption of an AU base pair and 6.6 × 1013 s−1 for the formation/ disruption of a GC base pair,35,40−43 kB is the Boltzmann constant, and T is the temperature. The rates for formation and = k0 disruption of a loop-closing stack (and the loop) are kloop + e−(ΔSstack+ΔSloop)/kB and kloop = k0 e−ΔHstack/kBT, where −ΔSstack and − −ΔHstack are the entropy and enthalpy changes upon formation of the stack, respectively, and −ΔSloop is the entropy change of the loop. Recently, a helix-based kinetic model was proposed to predict the folding kinetics of long RNA sequences and had been used for a few sequences.40,41,43 In this model, the conformation space is constructed by using helices as the building blocks; the kinetic move is the addition or deletion of a helix or an exchange between two helices. If two states differ in one helix, the transition between them will be the formation and disruption of the helix. The most probable pathway for the helix formation is the zipping pathway, and the formation of the first stack is rate-limiting in the helix formation (Figure 2a). The zipping rate, kf, along the 1 → 2 → 3 pathway in Figure 2b can be calculated as43

Figure 1. Secondary structure model of the E. faecalis SMK riboswitch. The nucleotides within the paired region of helices P1, P2, P3, and P4 are colored differently: P1 (fuchsia), P2 (green), P3 (blue), and P4 (red). The SMK 84 construct lacks the first six nucleotides at the 5′ terminus (gray box).

first six nucleotides from the 5′ terminus of the RNA.32 We found that the riboswitch folded into the ON-state without obvious misfolded structures formed during both the transcription and the refolding processes. Although the transition between ON-state and OFF-state required a great conformational rearrangement, the fast helix-exchanging transition pathway guaranteed that it could function as a reversible switch in response to the ligand concentration.

2. THEORETICAL METHODS 2.1. Free-Folding Kinetics. For a given RNA chain, the conformation space is sampled by all possible secondary structures. The free energy of secondary structures in the

Figure 2. (a) Schematic free-energy landscape for a pathway from the open chain to a helix. (b) Formation of a helix along the zipping pathways. The first closing stack is denoted by the red circle. (c) Free-energy landscape for the transition between helices A and B. Solid line: the tunneling pathway. Dotted line: complete unfolding of helix A followed by refolding to helix B. The 3′ end of the RNA chain is denoted by purple color. 12306

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B ∞

k f = k1 → 2K1(1 − K 2′K1′∑ (K 2′K1)n ) n=0

⎞ ⎛ 1 = k1 → 2K1⎜1 − K 2′K1′ ⎟ 1 − K 2′K1 ⎠ ⎝

(2)

where ki→j denotes the rate for the transition from state i to state j, and Ki and Ki′ are the forward and reverse probabilities, respectively k 2→3 k 2→1 , K1′ = , k 2→3 + k 2→1 + k 2→4 k 2→3 + k 2→1 + k 2→4 k3→5 + k3→6 k3→2 K2 = , K 2′ = k3→2 + k3→5 + k3→6 k3→2 + k3→5 + k3→6

K1 =

Figure 3. Four types of relationships between N-nt and (N + 1)-nt structures. The triangle and circle denote the last transcribed nucleotides at steps N and N + 1, respectively.

Rate kF for the helix formation is the sum of rates along all zipping pathways with different nucleation base pairs (Figure 2b). The rate for the helix deletion is calculated from the detailed balance condition, kU = kF e−ΔG/kBT, where ΔG is the folding free energy of the helix. If two helices A and B overlap with each other, they cannot coexist in the same structure. Compared with fully unfolding one helix followed by forming the other helix, the exchanging transition between the two helices occurs mostly through the tunneling pathway and has a much lower barrier (Figure 2c), wherein after one helix is partially disrupted, in each subsequent step, disruption of a stack in the helix is accompanied by the formation of a stack in the other helix. The rates between helices A and B for helix exchange can be estimated as43

beginning of step N + 1 can be summarized by the equations below32 P(N + 1)begin = P(N )end for a , b , and c ; P(N + 1)begin = 0 for d

For consecutive steps, the folding results of the previous step turn into the initial condition of the next step. Applying this method from the first step to the end of the transcription, we can compute the folding kinetics for the RNA chain during the transcription.

n

kA → B =

∏i ki n−1 j n ∑ j = 0 (∏i = 1 ki′∏m = j + 2 km)

kB → A = kA → B e−ΔGAB / kBT

3. RESULTS AND DISCUSSION 3.1. Co-transcriptional Folding and Refolding Behavior of the SMK Riboswitch. To investigate whether the action of the SMK riboswitch in cells depends on the transcription process, we first calculate the co-transcriptional folding behavior of the full-length SMK riboswitch without the ligand at a typical elongation rate of 50 nt/s. As the chain grows, the nascent RNA chain folds into a series of discrete intermediate states (Figure 4). Along the transcription, helix P5 (state C2) is formed as the first 19 nt are transcribed. Almost all riboswitches stay in this structure till the 36th nucleotide is transcribed, from which a new helix can be nucleated and the structure transits to C3. At the 44th step, the upper part of helix P3 is nucleated, and state C4, which consists of helices P5 and P3, quickly occupies most of the population. At about step 60, helix P0 becomes stable, and state C4 transits to the ON-state, which consists of helices P3, P5, and P0, by forming the nonlocal helix P0. At step 63, the ON-state is fully formed and occupies almost all of the population. From step 88, the repression stem P1 could be nucleated, and less than 2% of the population would transit from the ON-state to the OFF-state, as the ONstate with a free energy of −18.80 kcal/mol is more stable than the OFF-state (ΔGOFF = −17.40 kcal/mol). Therefore, in the absence of the ligand, the SMK riboswitch folds through a series of discrete intermediate states and finally populates at the ONstate with the SD sequence unstructured during the transcription process. This implies that it acts as a genetic ON switch to induce downstream gene expression without the ligand, in agreement with previous studies.25,26 The refolding kinetics of the full-length SAM-III riboswitch also shows that the thermodynamically favored ON-state is formed in less than 0.1 s, whereas the OFF-state only occupies a small proportion of the population (Figure 5). The results in

, (3)

where kn and k′n are the rate constants for the processes of formation (disruption) and disruption (formation) of a base stack in A (B), respectively. 2.2. Co-transcriptional Folding Kinetics. In the theory,32 the transcript of an L-nt RNA chain is divided into L transcriptional steps, each corresponding to the release of one nucleotide by RNA polymerase (RNAP) to freely form secondary structures. If the RNAP elongates the RNA chain with a speed of ν nucleotides/s, the (real) time window for each step will be 1/v s, that is, the polymerase spends 1/v s to synthesize a nucleotide. From time t when an N-nt chain is newly transcribed to time t + 1/v when the (N + 1)th nucleotide is (newly) transcribed, the N-nt chain samples the conformation space, and its population distribution is relaxed from [p1(N)begin, p2(N)begin, ..., pΩ(N)begin] to [p1(N)end, p2(N)end, ..., pΩ(N)end], where pi(N)begin and pi(N)end are the population of state i at the beginning and end of step N, respectively. At each such transcriptional step, the population kinetics is calculated in the same manner as the free-folding kinetics: first, the conformation space is generated, and the transition rates are calculated; then, the population relaxation within the folding time window is described by eq 1. It is important to note that the initial population at the current step is determined by the final population distribution at the previous step. Depending on possible structural changes upon elongation of the chain by one nucleotide, the structures are classified into four types (Figure 3). Then the population distribution at the 12307

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B

Figure 4. Population kinetics of main states (C1−C4) formed during the transcription of the SMK riboswitch without SAM at a transcription rate of 50 nt/s (a). (b) The transition pathways and structures formed during the transcription.

on the transcription process. Taken together, these results suggested, unlike the riboswitches operating at the transcription level,15,37 that the coupling between transcription and translation was not required for the SAM-III riboswitch to efficiently perform gene regulation. The ligand should bind the target to induce the function switch after the full-length chain is transcribed, and the flipping of the switch could be limited by ON-state and OFF-state transitions. 3.2. Helix P0 Is Crucial for the SMK Riboswitch To Perform Its Function. Previous study suggested that variable sequences outside the SAM-binding core critically influenced the conformational dynamics of the SAM-III riboswitch.28 The RNA construct lacking the first several nucleotides at the 5′ end of the RNA was found to favor a conformation resembling the genetic OFF-state even in the absence of SAM.27 To study the effect of these nucleotides on the function of the riboswitch, we predict the co-transcriptional folding behavior of an 84-nt SMK riboswitch construct lacking the first 6 nucleotides (Figure 6). It shows the same folding behavior as that of the full-length riboswitch before the 60th step (Figure 4), from which the fulllength RNA transits from structure C4 to the ON-state by forming helix P0. However, for the 84-nt construct, helix P0 is unable to form due to the lack of nucleotides. C4 occupies more than 98% of the population for a long transcription time, until the 72nd nucleotide is transcribed, when a new helix can

Figure 5. Refolding population kinetics of the 90-nt full-length SMK riboswitch.

Figures 4 and 5 showed that during both the free folding and the co-transcriptional folding processes, almost all the SAM-III riboswitches quickly folded into the thermodynamically favored ON-state. Hence, the regulation of its downstream gene translation by the structure of the riboswitch does not depend

Figure 6. Co-transcriptional folding kinetics of the 84-nt SMK riboswitch construct in the absence of SAM with 50 nt/s; population relaxation after the transcription is shown at the right side of the dotted line, where the end of step 84 is set as time 0 s. 12308

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B

transition from ON-state to OFF-state should be extremely slow through completely unfolding of helices P0 and P5, followed by refolding of helices P1, P2, and P4. The tunneling pathway, in which at first the helices of the ON-state are partially disrupted and, in each subsequent step, disruption of a base stack in the ON-state is accompanied by a formation of a base stack in the OFF-state, is much faster with a lower energy barrier.43 So the most probable transition pathway is the helixexchanging tunneling pathway (Figure 8a). According to eqs 2 and 3,43 the transition rate from ON-state to OFF-state and the reverse rate are 0.04 and 0.38 s−1, respectively, at 37 °C. The binding of the ligand to the full-length riboswitch can be described as a three-state transition model: ON-state ↔ OFFstate + SAM ↔ OFFb-state, where OFFb-state is the OFF-state with ligand bound. Ligand binding can stabilize the ligandbound OFF-state by ΔGbinding = kBT ln(kon[SAM]/koff) kcal/ mol, where [SAM] is the concentration of SAM, and kon and koff are the association and disassociation rates, respectively, which have been experimentally measured for the full-length riboswitch as kon = 0.11 μM−1 s−1 and koff = 0.089 s−1.25 The transition rate from OFF-state to OFFb-state was assumed to be the effective binding rate, keff = kon[SAM], and the transition rate from OFFb-state to OFF-state is the dissociation rate, koff. Assuming the initial state is the ON-state, most riboswitches could transit to the OFFb-state within 100 s in the presence of 500 μM SAM (Figure 8b), which is in agreement with the experimental measured value of 1 min.28 When the ligand concentration decreases to 1 μM, more than half of the riboswitch would transit from the genetic OFF-state (OFFbstate) to the genetic ON-state (ON-state) within 60 s (Figure 8c). These results suggest that the SMK riboswitch can switch back and forth to regulate translation of its host gene quickly in response to the ligand concentration. The major time scales that control the function of the SMK riboswitch are (1) the ligand-binding rate, kon = 0.11 μM−1 s−1, (2) the dissociation rate, koff = 0.089 s−1, and (3) the mRNA degradation rate, kmRNA = 3 min−1.47 A thermodynamic regime is that the time needed to attain an RNA−ligand equilibrium is short compared with the transcriptional time scale. The switching rate for the SMK riboswitch is faster than the degradation rate of the riboswitch, so the functional status of the SMK riboswitch at any given time is governed by thermodynamic principles.

be formed and the population of structure C4 begins to transit to structure C5, which does not form for the full-length riboswitch (Figure 4). As the 77th nucleotide is free to form structure, structure C5 begins to transit to the OFF-state. At the end of the transcription, structure C5 and the OFF-state occupy about 52.5 and 45.1% of the population, respectively. The free energy of structure C5 and the OFF-state are −16.32 and −17.40 kcal/mol, respectively, so they do not reach equilibrium at the end of transcription. Another 0.5 s are required to significantly reach equilibrium, and then the OFFstate occupies about 80% of the population. Although the ribosome-binding time is not well known, the time required to transcribe the mRNA footprint for the ribosome can provide a lower bound.44 Crystal structures and footprinting studies suggest that this involves another 15-nt downstream of the 3′ end of the riboswitch expression platform,45,46 requiring about 0.3 s at typical transcription rates. During the transcription of the 84-nt RNA construct, the ON-state is not formed due to the lack of the nucleotides at the left shoulder of helix P0. The refolding process of the 84-nt RNA also shows that structure C5 as an intermediate can be formed, but it is quickly replaced by the OFF-state (Figure 7), and the thermodynami-

Figure 7. Refolding population kinetics of the 84-nt SMK riboswitch.

4. CONCLUSIONS Riboswitches, as typical functional mRNAs, which are synthesized during the transcription process, adopt alternative folds to control downstream gene expression.48−50 It has been found that for many riboswitches, their functions are sensitive to co-transcriptional folding events.15,29,35,37 To understand the functions of riboswitches, it is important and necessary to investigate their co-transcriptional folding behaviors. However, sequential folding during the transcription of mRNAs poses a serious challenge in detecting these structures as they are formed within a short time.30,51 Using the recently developed systematic helix-based computational method to predict the co-transcriptional folding kinetics,32 we theoretically studied the co-transcriptional folding behavior and refolding kinetics of the E. faecalis SAM-III riboswitch. The riboswitch, with a simple architecture to bind the ligand and regulate gene expression, displays many unique features in its regulation. Both the refolding and the cotranscriptional folding kinetics of the full-length riboswitch showed that most RNAs predominantly folded into the ON-

cally favored OFF-state can be formed in less than 1 s. This implies that most riboswitches fold into the OFF-state before the ribosome can bind to the riboswitch. Thus, the 84-nt SMK riboswitch shortened construct functions as a “constitutiveOFF” switch to repress downstream gene expression even in the absence of SAM in cells, consistent with experimental observation.28 In other words, the potential to form helix P0 is crucial for the SMK riboswitch to perform its switch function. 3.3. ON−OFF-State Transition Is through a HelixExchanging Pathway. The free-folding kinetics and the cotranscriptional folding kinetics indicated that without SAM, the SMK riboswitch would predominantly fold to the thermodynamically favored ON-state. The OFF-state was suggested to be a SAM-binding-competent structure.28 But the transition from ON-state to OFF-state needs to undergo enormous conformational rearrangements, as it needs to break helices P0 and P5 and form helices P1, P2, and P4. What is the transition pathway, and is the transition rate fast enough for this reversible riboswitch to function? Helices P0 and P5 of the ON-state overlap with those of P1, P2, P4 of the OFF-state. The 12309

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B

Figure 8. (a) Exchanging transition pathway between the ON-state and the OFF-state and the schematic free-energy landscape of the pathway (unit kcal/mol). (b) Transition from ON-state to OFFb-state through the pathway: ON ↔ OFF + SAM ↔ OFFb-state in the presence of 500 μM SAM. (c) Transition from OFFb-state to ON-state when the ligand concentration decreases to 1 μM.



ACKNOWLEDGMENTS This work was partly supported by the National Natural Science Foundation of China under grant nos. 31270761 and 1157324 (W.Z.).

state. This implies that the function of the riboswitch is not tied to the transcription process and that the ligand binding occurs after the transcription. Without ligand, the ON-state is the most stable state. Upon SAM binding, the equilibrium is shifted toward the thermodynamically favorable ligand-bound OFFstate with the repression stem formed. The SMK riboswitch exerts its function by thermodynamic control, and its switch efficiency is linked to the stability of the two function structures instead of the transcription context. Although the ON−OFF transition requires enormous conformation rearrangements, the fast helix-exchanging transition pathway between the two function structures allows it to perform as a reversible switch in response to the ligand concentration. The shortened or mutated RNA constructs, which disrupt helix P0 of the ON-state, are not capable of switching and behave as a constitutive-OFF switch. These results demonstrated that the riboswitch is highly evolved and has selected sequences. However, the current method has limitations. The tertiary interactions and ion effects,52,53 which may alter the folding pathways, are neglected. Future development of the model would incorporate these effects into the model.





REFERENCES

(1) Gong, Z.; Zhao, Y.; Chen, C.; Xiao, Y. Role of Ligand Binding in Structural Organization of Add A-Riboswitch Aptamer: A Molecular Dynamics Simulation. J. Biomol. Struct. Dyn. 2011, 29, 403−416. (2) Lin, J. C.; Thirumalai, D. Kinetics of Allosteric Transitions in SAdenosylmethionine Riboswitch Are Accurately Predicted from the Folding Landscape. J. Am. Chem. Soc. 2013, 135, 16641−16650. (3) Mellin, J. R.; Koutero, M.; Dar, D.; Nahori, M. a.; Sorek, R.; Cossart, P. Riboswitches. Sequestration of a Two-Component Response Regulator by a Riboswitch-Regulated Noncoding RNA. Science 2014, 345, 940−943. (4) Reining, A.; Nozinovic, S.; Schlepckow, K.; Buhr, F.; Fürtig, B.; Schwalbe, H. Three-State Mechanism Couples Ligand and Temperature Sensing in Riboswitches. Nature 2013, 499, 355−359. (5) Roth, A.; Breaker, R. R. The Structural and Functional Diversity of Metabolite-Binding Riboswitches. Annu. Rev. Biochem. 2009, 78, 305−334. (6) Serganov, A.; Nudler, E. A Decade of Riboswitches. Cell 2013, 152, 17−24. (7) Winkler, W. C. Riboswitches and the Role of Noncoding RNAs in Bacterial Metabolic Control. Curr. Opin. Chem. Biol. 2005, 9, 594− 602. (8) Breaker, R. R. Riboswitches and the RNA World. Cold Spring Harbor Perspect. Biol. 2012, 4, a003566. (9) Mandal, M.; Breaker, R. R. Gene Regulation by Riboswitches. Nat. Rev. Mol. Cell Biol. 2004, 5, 451−463. (10) Winkler, W. C.; Breaker, R. R. Genetic Control by MetaboliteBinding Riboswitches. ChemBioChem 2003, 4, 1024−1032.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Wenbing Zhang: 0000-0001-6880-7087 Notes

The authors declare no competing financial interest. 12310

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311

Article

The Journal of Physical Chemistry B

(32) Zhao, P.; Zhang, W.; Chen, S.-J. Cotranscriptional Folding Kinetics of Ribonucleic Acid Secondary Structures. J. Chem. Phys. 2011, 135, No. 245101. (33) Kramer, F. R.; Mills, D. R. Secondary Structure Formation during RNA Synthesis. Nucleic Acids Res. 1981, 9, 5109−5124. (34) Gong, S.; Wang, Y.; Zhang, W. The Regulation Mechanism of yitJ and metF Riboswitches. J. Chem. Phys. 2015, 143, No. 045103. (35) Gong, S.; Wang, Y.; Zhang, W. Kinetic Regulation Mechanism of pbuE Riboswitch. J. Chem. Phys. 2015, 142, No. 015103. (36) Tomšič, J.; McDaniel, B. A.; Grundy, F. J.; Henkin, T. M. Natural Variability in S-Adenosylmethionine (SAM)-Dependent Riboswitches: S-Box Elements in Bacillus subtilis Exhibit Differential Sensitivity to SAM in Vivo and in Vitro. J. Bacteriol. 2008, 190, 823− 833. (37) Wickiser, J. K.; Winkler, W. C.; Breaker, R. R.; Crothers, D. M. The Speed of RNA Transcription and Metabolite Binding Kinetics Operate an FMN Riboswitch. Mol. Cell 2005, 18, 49−60. (38) Xia, T.; Lucia, J. S., Jr.; Burkard, M. E.; Kierzek, R.; Schroeder, S. J.; Jiao, X.; Cox, C.; Turner, D. H. Thermodynamic Parameters for an Expanded Nearest-Neighbor Model for Formation of RNA Duplexes with Watson-Crick Base Pairs. Biochemistry 1998, 37, 14719−14735. (39) Mathews, D. H.; Sabina, J.; Zuker, M.; Turner, D. H. Expanded Sequence Dependence of Thermodynamic Parameters Improves Prediction of RNA Secondary Structure. J. Mol. Biol. 1999, 288, 911−940. (40) Chen, J.; Zhang, W. Kinetic Analysis of the Effects of Target Structure on siRNA Efficiency. J. Chem. Phys. 2012, 137, No. 225102. (41) Chen, J.; Gong, S.; Wang, Y.; Zhang, W. Kinetic Partitioning Mechanism of HDV Ribozyme Folding. J. Chem. Phys. 2014, 140, No. 025102. (42) Zhang, W.; Chen, S.-J. Exploring the Complex Folding Kinetics of RNA Hairpins: I. General Folding Kinetics Analysis. Biophys. J. 2006, 90, 765−777. (43) Zhao, P.; Zhang, W. B.; Chen, S. J. Predicting Secondary Structural Folding Kinetics for Nucleic Acids. Biophys. J. 2010, 98, 1617−1625. (44) Neupane, K.; Yu, H.; Foster, D. A. N.; Wang, F.; Woodside, M. T. Single-Molecule Force Spectroscopy of the Add Adenine Riboswitch Relates Folding to Regulatory Mechanism. Nucleic Acids Res. 2011, 39, 7677−7687. (45) Yusupova, G. Z.; Yusupov, M. M.; Cate, J. H. D.; Noller, H. F. The Path of Messenger RNA through the Ribosome. Cell 2001, 106, 233−241. (46) Hüttenhofer, A.; Noller, H. F. Footprinting mRNA-Ribosome Complexes with Chemical Probes. EMBO J. 1994, 13, 3892−3901. (47) Lin, J. C.; Yoon, J.; Hyeon, C.; Thirumalai, D. Using Simulations and Kinetic Network Models to Reveal the Dynamics and Functions of Riboswitches. Methods Enzymol. 2015, 553, 235−258. (48) Sudarsan, N.; Lee, E. R.; Weinberg, Z.; Moy, R. H.; Kim, J. N.; Link, K. H.; Breaker, R. R. Riboswitches in Eubacteria Sense the Second Messenger Cyclic Di-GMP. Science 2008, 321, 411−413. (49) Chen, J.; Gottesman, S. Riboswitch Regulates RNA. Science 2014, 345, 876−877. (50) Perdrizet, G. A., II; Artsimovitch, I.; Furman, R.; Sosnick, T. R.; Pan, T. Transcriptional Pausing Coordinates Folding of the Aptamer Domain and the Expression Platform of a Riboswitch. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 3323−3328. (51) Mahen, E. M.; Watson, P. Y.; Cottrell, J. W.; Fedor, M. J. mRNA Secondary Structures Fold Sequentially but Exchange Rapidly in Vivo. PLoS Biol. 2010, 8, No. e1000307. (52) Tan, Z.-J.; Chen, S.-J. Salt Contribution to RNA Tertiary Structure Folding Stability. Biophys. J. 2011, 101, 176−187. (53) Shi, Y. Z.; Wang, F. H.; Wu, Y. Y.; Tan, Z. J. A Coarse-Grained Model with Implicit Salt for RNAs: Predicting 3D Structure, Stability and Salt Effect. J. Chem. Phys. 2014, 141, No. 105102.

(11) Miranda-Ríos, J.; Navarro, M.; Soberon, M. A Conserved RNA Structure (Thi Box) Is Involved in Regulation of Thiamin Biosynthetic Gene. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 9736−9741. (12) Huang, W.; Kim, J.; Jha, S.; Aboul-ela, F. A Mechanism for SAdenosyl Methionine Assisted Formation of a Riboswitch Conformation: A Small Molecule with a Strong Arm. Nucleic Acids Res. 2009, 37, 6528−6539. (13) Winkler, W. C.; Cohen-Chalamish, S.; Breaker, R. R. An mRNA Structure That Controls Gene Expression by Binding FMN. Proc. Natl. Acad. Sci. U.S.A. 2002, 99, 15908−15913. (14) Winkler, W. C.; Nahvi, A.; Sudarsan, N.; Barrick, J. E.; Breaker, R. R. An mRNA Structure That Controls Gene Expression by Binding S-Adenosylmethionine. Nat. Struct. Biol. 2003, 10, 701−707. (15) Lemay, J.-F.; Desnoyers, G.; Blouin, S.; Heppell, B.; Bastet, L.; St-Pierre, P.; Massé, E.; Lafontaine, D. A. Comparative Study between Transcriptionally- and Translationally-Acting Adenine Riboswitches Reveals Key Differences in Riboswitch Regulatory Mechanisms. PLoS Genet. 2011, 7, No. e1001278. (16) Mandal, M.; Lee, M.; Barrick, J. E.; Weinberg, Z.; Emilsson, G. M.; Ruzzo, W. L.; Breaker, R. R. A Glycine-Dependent Riboswitch That Uses Cooperative Binding to Control Gene Expression. Science 2004, 306, 275−279. (17) McDaniel, B. A.; Grundy, F. J.; Henkin, T. M. A Tertiary Structural Element in S Box Leader RNAs Is Required for SAdenosylmethionine-Directed Transcription Termination. Mol. Microbiol. 2005, 57, 1008−1021. (18) Mandal, M.; Breaker, R. R. Adenine Riboswitches and Gene Activation by Disruption of a Transcription Terminator. Nat. Struct. Mol. Biol. 2004, 11, 29−35. (19) Tucker, B. J.; Breaker, R. R. Riboswitches as Versatile Gene Control Elements. Curr. Opin. Struct. Biol. 2005, 15, 342−348. (20) Doshi, U.; Kelley, J. M.; Hamelberg, D. Atomic-Level Insights into Metabolite Recognition and Specificity of the SAM-II Riboswitch. RNA 2012, 18, 300−307. (21) Priyakumar, U. D. Atomistic Details of the Ligand Discrimination Mechanism of S MK/SAM-III Riboswitch. J. Phys. Chem. B 2010, 114, 9920−9925. (22) Batey, R. T. Recognition of S-Adenosylmethionine by Riboswitches. Wiley Interdiscip. Rev.: RNA 2011, 2, 299−311. (23) Fuchs, R. T.; Grundy, F. J.; Henkin, T. M. S-Adenosylmethionine Directly Inhibits Binding of 30 S Ribosomal Subunits to the SMK Box Translational Riboswitch RNA. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 4876−4880. (24) Lu, C.; Smith, A. M.; Fuchs, R. T.; Ding, F.; Rajashankar, K.; Henkin, T. M.; Ke, A. Crystal Structures of the SAM-III/S(MK) Riboswitch Reveal the SAM-Dependent Translation Inhibition Mechanism. Nat. Struct. Mol. Biol. 2008, 15, 1076−1083. (25) Smith, A. M.; Fuchs, R. T.; Grundy, F. J.; Henkin, T. M. The SAM-Responsive S(MK) Box Is a Reversible Riboswitch. Mol. Microbiol. 2010, 78, 1393−1402. (26) Fuchs, R. T.; Grundy, F. J.; Henkin, T. M. The S(MK) Box Is a New SAM-Binding RNA for Translational Regulation of SAM Synthetase. Nat. Struct. Mol. Biol. 2006, 13, 226−233. (27) Wilson, R. C.; Smith, A. M.; Fuchs, R. T.; Kleckner, I. R.; Henkin, T. M.; Foster, M. P. Tuning Riboswitch Regulation through Conformational Selection. J. Mol. Biol. 2011, 405, 926−938. (28) Lu, C.; Smith, A. M.; Ding, F.; Chowdhury, A.; Henkin, T. M.; Ke, A. Variable Sequences Outside the SAM-Binding Core Critically Influence the Conformational Dynamics of the SAM-III/SMK Box Riboswitch. J. Mol. Biol. 2011, 409, 786−799. (29) Lutz, B.; Faber, M.; Verma, A.; Klumpp, S.; Schug, A. Differences between Cotranscriptional and Free Riboswitch Folding. Nucleic Acids Res. 2014, 42, 2687−2696. (30) Frieda, K. L.; Block, S. M. Direct Observation of Cotranscriptional Folding in an Adenine Riboswitch. Science 2012, 338, 397−400. (31) Wong, T. N.; Sosnick, T. R.; Pan, T. Folding of Noncoding RNAs during Transcription Facilitated by Pausing-Induced Nonnative Structures. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 17995−18000. 12311

DOI: 10.1021/acs.jpcb.6b09698 J. Phys. Chem. B 2016, 120, 12305−12311