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