Search

Posttranslational gene regulation requires a complex network of RNA–protein interactions. Cooperativity, which tunes response sensitivities, originates from protein–protein interactions in many systems. For RNA-binding proteins, cooperativity can also be mediated through RNA structure. RNA structural cooperativity (RSC) arises when binding of one protein induces a redistribution of RNA conformational states that enhance access (positive cooperativity) or block access (negative cooperativity) to additional binding sites. As RSC does not require direct protein–protein interactions, it allows cooperativity to be tuned for individual RNAs, via alterations in sequence that alter structural stability. Given the potential importance of this mechanism of control and our desire to quantitatively dissect features that underlie physiological regulation, we developed a statistical mechanical framework for RSC and tested this model by performing equilibrium binding measurements of the human PUF family protein PUM2. Using 68 RNAs that contain two to five PUM2-binding sites and RNA structures of varying stabilities, we observed a range of structure-dependent cooperative behaviors. To test our ability to account for this cooperativity with known physical constants, we used PUM2 affinity and nearest-neighbor RNA secondary structure predictions. Our model gave qualitative agreement for our disparate set of 68 RNAs across two temperatures, but quantitative deviations arise from overestimation of RNA structural stability. Our results demonstrate cooperativity mediated by RNA structure and underscore the power of quantitative stepwise experimental evaluation of mechanisms and computational tools.

Nearest-neighbor (NN) rules provide a simple and powerful quantitative framework for RNA structure prediction that is strongly supported for canonical Watson-Crick duplexes from a plethora of thermodynamic measurements. Predictions of RNA secondary structure based on nearest-neighbor (NN) rules are routinely used to understand biological function and to engineer and control new functions in biotechnology. However, NN applications to RNA structural features such as internal and terminal loops rely on approximations and assumptions, with sparse experimental coverage of the vast number of possible sequence and structural features. To test to what extent NN rules accurately predict thermodynamic stabilities across RNAs with non-WC features, we tested their predictions using a quantitative high-throughput assay platform, RNA-MaP. Using a thermodynamic assay with coupled protein binding, we carried out equilibrium measurements for over 1000 RNAs with a range of predicted secondary structure stabilities. Our results revealed substantial scatter and systematic deviations between NN predictions and observed stabilities. Solution salt effects and incorrect or omitted loop parameters contribute to these observed deviations. Our results demonstrate the need to independently and quantitatively test NN computational algorithms to identify their capabilities and limitations. RNA-MaP and related approaches can be used to test computational predictions and can be adapted to obtain experimental data to improve RNA secondary structure and other prediction algorithms.Significance statementRNA secondary structure prediction algorithms are routinely used to understand, predict and design functional RNA structures in biology and biotechnology. Given the vast number of RNA sequence and structural features, these predictions rely on a series of approximations, and independent tests are needed to quantitatively evaluate the accuracy of predicted RNA structural stabilities. Here we measure the stabilities of over 1000 RNA constructs by using a coupled protein binding assay. Our results reveal substantial deviations from the RNA stabilities predicted by popular algorithms, and identify factors contributing to the observed deviations. We demonstrate the importance of quantitative, experimental tests of computational RNA structure predictions and present an approach that can be used to routinely test and improve the prediction accuracy.

RNA-guided nucleases (RGNs) provide sequence-specific gene regulation through base-pairing interactions between a small RNA guide and target RNA or DNA. RGN systems, which include CRISPR-Cas9 and RNA interference (RNAi), hold tremendous promise as programmable tools for engineering and therapeutic purposes. However, pervasive targeting of sequences that closely resemble the intended target has remained a major challenge, limiting the reliability and interpretation of RGN activity and the range of possible applications. Efforts to reduce off-target activity and enhance RGN specificity have led to a collection of empirically derived rules, which often paradoxically include decreased binding affinity of the RNA-guided nuclease to its target. Here we demonstrate that simple kinetic considerations of the targeting reaction can explain these and other literature observations. The kinetic models described provide a foundation for understanding RGN systems and a necessary physical and functional framework for their rational engineering.