Genomic variants influence phenotypes and organismal fitness, with their effects shaped by genomic context. In 3' untranslated regions, variants can alter phenotypes by influencing RNA-protein binding and subsequent post-transcriptional gene regulation. Here, we investigate how indel variants impact RNA-protein interactions from outside the binding site, through changes in RNA secondary structure. Our findings reveal that indels can significantly affect protein binding affinities over distances spanning tens of nucleotides, with longer indels exerting greater effects until saturation. Crucially, we find evidence that this effect is constrained by purifying selection. Naturally observed indel polymorphisms cause smaller changes in binding affinity than synthetic indels, and this signal of selection is more pronounced near HuR binding sites. We also find that the sequence context of HuR binding sites shows greater resilience to indel mutations compared to randomly selected sites. These results demonstrate that indel polymorphisms outside the binding site can modulate RNA-protein interactions through structural effects, and that purifying selection acts to filter out variants that disrupt critical interactions.

MOTIVATION: RNA-binding proteins are fundamental to many cellular processes. Double-stranded RNA-binding proteins (dsRBPs) in particular are crucial for RNA interference, mRNA elongation, A-to-I editing, host defense, splicing and a multitude of other important mechanisms. Since dsRBPs require double-stranded RNA to bind, their binding affinity depends on the competition among all possible secondary structures of the target RNA molecule. Here, we introduce a quantitative model that allows calculation of the effective affinity of dsRBPs to any RNA given a principal affinity and the sequence of the RNA, while fully taking into account the entire secondary structure ensemble of the RNA. RESULTS: We implement our model within the ViennaRNA folding package while maintaining its O(N3) time complexity. We validate our quantitative model by comparing with experimentally determined binding affinities and stoichiometries for transactivation response element RNA-binding protein (TRBP). We also find that the change in dsRBP binding affinity purely due to the presence of alternative RNA structures can be many orders of magnitude and that the predicted affinity of TRBP for pre-miRNA-like constructs correlates with experimentally measured processing rates. AVAILABILITY AND IMPLEMENTATION: Our modified version of the ViennaRNA package is available for download at http://bioserv.mps.ohio-state.edu/dsRBPBind, is free to use for research and educational purposes, and utilizes simple get/set methods for footprint size, concentration, cooperativity, principal dissociation constant and overlap.

Single nucleotide polymorphisms are widely associated with disease, but the ways in which they cause altered phenotypes are often unclear, especially when they appear in non-coding regions. One way in which non-coding polymorphisms could cause disease is by affecting crucial RNA-protein interactions. While it is clear that changing a protein binding motif will alter protein binding, it has been shown that single nucleotide polymorphisms can affect RNA secondary structure, and here we show that single nucleotide polymorphisms can affect RNA-protein interactions from outside binding motifs through altered RNA secondary structure. By using a modified version of the Vienna Package and PAR-CLIP data for HuR (ELAVL1) in humans we characterize the genome-wide effect of single nucleotide polymorphisms on HuR binding and show that they can have a many-fold effect on the affinity of HuR binding to RNA transcripts from tens of bases away. We also find some evidence that the effect of single nucleotide polymorphisms on protein binding might be under selection, with the non-reference alleles tending to make it harder for a protein to bind.