Design principles for riboswitch function

Scientific and technological advances that enable the tuning of integrated regulatory components to match network and system requirements are critical to reliably control the function of biological systems. RNA provides a promising building block for the construction of tunable regulatory components based on its rich regulatory capacity and our current understanding of the sequence–function relationship. One prominent example of RNA-based regulatory components is riboswitches, genetic elements that mediate ligand control of gene expression through diverse regulatory mechanisms. While characterization of natural and synthetic riboswitches has revealed that riboswitch function can be modulated through sequence alteration, no quantitative frameworks exist to investigate or guide riboswitch tuning. Here, we combined mathematical modeling and experimental approaches to investigate the relationship between riboswitch function and performance. Model results demonstrated that the competition between reversible and irreversible rate constants dictates performance for different regulatory mechanisms. We also found that practical system restrictions, such as an upper limit on ligand concentration, can significantly alter the requirements for riboswitch performance, necessitating alternative tuning strategies. Previous experimental data for natural and synthetic riboswitches as well as experiments conducted in this work support model predictions. From our results, we developed a set of general design principles for synthetic riboswitches. Our results also provide a foundation from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands

Dynamic Energy Landscapes of Riboswitches Help Interpret Conformational Rearrangements and Function

Riboswitches are RNAs that modulate gene expression by ligand-induced conformational changes. However, the way in which sequence dictates alternative folding pathways of gene regulation remains unclear. In this study, we compute energy landscapes, which describe the accessible secondary structures for a range of sequence lengths, to analyze the transcriptional process as a given sequence elongates to full length. In line with experimental evidence, we find that most riboswitch landscapes can be characterized by three broad classes as a function of sequence length in terms of the distribution and barrier type of the conformational clusters: low-barrier landscape with an ensemble of different conformations in equilibrium before encountering a substrate; barrier-free landscape in which a direct, dominant “downhill” pathway to the minimum free energy structure is apparent; and a barrier-dominated landscape with two isolated conformational states, each associated with a different biological function. Sharing concepts with the “new view” of protein folding energy landscapes, we term the three sequence ranges above as the sensing, downhill folding, and functional windows, respectively. We find that these energy landscape patterns are conserved in various riboswitch classes, though the order of the windows may vary. In fact, the order of the three windows suggests either kinetic or thermodynamic control of ligand binding. These findings help understand riboswitch structure/function relationships and open new avenues to riboswitch design

Molecular insights into the ligand-controlled organization of SAM-I riboswitch

S-adenosylmethionine (SAM) riboswitches are widespread in bacteria, and up to five different SAM riboswitch families have been reported, highlighting the relevance of SAM regulation. On the basis of crystallographic and biochemical data, it has been postulated, but never demonstrated, that ligand recognition by SAM riboswitches involves key conformational changes in the RNA architecture. We show here that the aptamer follows a two-step hierarchical folding selectively induced by metal ions and ligand binding, each of them leading to the formation of one of the two helical stacks observed in the crystal structure. Moreover, we find that the anti-antiterminator P1 stem is rotated along its helical axis upon ligand binding, a mechanistic feature that could be common to other riboswitches. We also show that the nonconserved P4 helical domain is used as an auxiliary element to enhance the ligand-binding affinity. This work provides the first comprehensive characterization, to our knowledge, of a ligand-controlled riboswitch folding pathway