A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins.

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

Blind tests of RNA-protein binding affinity prediction.

Interactions between RNA and proteins are pervasive in biology, driving fundamental processes such as protein translation and participating in the regulation of gene expression. Modeling the energies of RNA-protein interactions is therefore critical for understanding and repurposing living systems but has been hindered by complexities unique to RNA-protein binding. Here, we bring together several advances to complete a calculation framework for RNA-protein binding affinities, including a unified free energy function for bound complexes, automated Rosetta modeling of mutations, and use of secondary structure-based energetic calculations to model unbound RNA states. The resulting Rosetta-Vienna RNP-ΔΔG method achieves root-mean-squared errors (RMSEs) of 1.3 kcal/mol on high-throughput MS2 coat protein-RNA measurements and 1.5 kcal/mol on an independent test set involving the signal recognition particle, human U1A, PUM1, and FOX-1. As a stringent test, the method achieves RMSE accuracy of 1.4 kcal/mol in blind predictions of hundreds of human PUM2-RNA relative binding affinities. Overall, these RMSE accuracies are significantly better than those attained by prior structure-based approaches applied to the same systems. Importantly, Rosetta-Vienna RNP-ΔΔG establishes a framework for further improvements in modeling RNA-protein binding that can be tested by prospective high-throughput measurements on new systems.

Demonstration of protein cooperativity mediated by RNA structure using the human protein PUM2.

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.

Lso2 is a conserved ribosome-bound protein required for translational recovery in yeast.

Ribosome-binding proteins function broadly in protein synthesis, gene regulation, and cellular homeostasis, but the complete complement of functional ribosome-bound proteins remains unknown. Using quantitative mass spectrometry, we identified late-annotated short open reading frame 2 (Lso2) as a ribosome-associated protein that is broadly conserved in eukaryotes. Genome-wide crosslinking and immunoprecipitation of Lso2 and its human ortholog coiled-coil domain containing 124 (CCDC124) recovered 25S ribosomal RNA in a region near the A site that overlaps the GTPase activation center. Consistent with this location, Lso2 also crosslinked to most tRNAs. Ribosome profiling of yeast lacking LSO2 (lso2Δ) revealed global translation defects during recovery from stationary phase with translation of most genes reduced more than 4-fold. Ribosomes accumulated at start codons, were depleted from stop codons, and showed codon-specific changes in occupancy in lso2Δ. These defects, and the conservation of the specific ribosome-binding activity of Lso2/CCDC124, indicate broadly important functions in translation and physiology.

Pseudouridine and N6-methyladenosine modifications weaken PUF protein/RNA interactions.

RNA modifications are ubiquitous in biology, with over 100 distinct modifications. While the vast majority were identified and characterized on abundant noncoding RNA such as tRNA and rRNA, the advent of sensitive sequencing-based approaches has led to the discovery of extensive and regulated modification of eukaryotic messenger RNAs as well. The two most abundant mRNA modifications-pseudouridine (Ψ) and N6-methyladenosine (m6A)-affect diverse cellular processes including mRNA splicing, localization, translation, and decay and modulate RNA structure. Here, we test the hypothesis that RNA modifications directly affect interactions between RNA-binding proteins and target RNA. We show that Ψ and m6A weaken the binding of the human single-stranded RNA binding protein Pumilio 2 (hPUM2) to its consensus motif, with individual modifications having effects up to approximately threefold and multiple modifications giving larger effects. While there are likely to be some cases where RNA modifications essentially fully ablate protein binding, here we see modest responses that may be more common. Such modest effects could nevertheless profoundly alter the complex landscape of RNA:protein interactions, and the quantitative rather than qualitative nature of these effects underscores the need for quantitative, systems-level accounting of RNA:protein interactions to understand post-transcriptional regulation.

Protein kinase A regulates gene-specific translational adaptation in differentiating yeast.

Cellular differentiation is driven by coordinately regulated changes in gene expression. Recent discoveries suggest that translation contributes as much as transcription to regulating protein abundance, but the role of translational regulation in cellular differentiation is largely unexplored. Here we investigate translational reprogramming in yeast during cellular adaptation to the absence of glucose, a stimulus that induces invasive filamentous differentiation. Using ribosome footprint profiling and RNA sequencing to assay gene-specific translation activity genome-wide, we show that prolonged glucose withdrawal is accompanied by gene-specific changes in translational efficiency that significantly affect expression of the majority of genes. Notably, transcripts from a small minority (<5%) of genes make up the majority of translating mRNA in both rapidly dividing and starved differentiating cells, and the identities of these highly translated messages are almost nonoverlapping between conditions. Furthermore, these two groups of messages are subject to condition-dependent translational privilege. Thus the "housekeeping" process of translation does not stay constant during cellular differentiation but is highly adapted to different growth conditions. By comparing glucose starvation to growth-attenuating stresses that do not induce invasive filamentation, we distinguish a glucose-specific translational response mediated through signaling by protein kinase A (PKA). Together, these findings reveal a high degree of growth-state specialization of the translatome and identify PKA as an important regulator of gene-specific translation activity.

A comprehensive thermodynamic model for RNA binding by the Saccharomyces cerevisiae Pumilio protein PUF4.

Genomic methods have been valuable for identifying RNA-binding proteins (RBPs) and the genes, pathways, and processes they regulate. Nevertheless, standard motif descriptions cannot be used to predict all RNA targets or test quantitative models for cellular interactions and regulation. We present a complete thermodynamic model for RNA binding to the S. cerevisiae Pumilio protein PUF4 derived from direct binding data for 6180 RNAs measured using the RNA on a massively parallel array (RNA-MaP) platform. The PUF4 model is highly similar to that of the related RBPs, human PUM2 and PUM1, with one marked exception: a single favorable site of base flipping for PUF4, such that PUF4 preferentially binds to a non-contiguous series of residues. These results are foundational for developing and testing cellular models of RNA-RBP interactions and function, for engineering RBPs, for understanding the biophysical nature of RBP binding and the evolutionary landscape of RNAs and RBPs.

How to measure and evaluate binding affinities.

Quantitative measurements of biomolecule associations are central to biological understanding and are needed to build and test predictive and mechanistic models. Given the advances in high-throughput technologies and the projected increase in the availability of binding data, we found it especially timely to evaluate the current standards for performing and reporting binding measurements. A review of 100 studies revealed that in most cases essential controls for establishing the appropriate incubation time and concentration regime were not documented, making it impossible to determine measurement reliability. Moreover, several reported affinities could be concluded to be incorrect, thereby impacting biological interpretations. Given these challenges, we provide a framework for a broad range of researchers to evaluate, teach about, perform, and clearly document high-quality equilibrium binding measurements. We apply this framework and explain underlying fundamental concepts through experimental examples with the RNA-binding protein Puf4.

RluD, a highly conserved pseudouridine synthase, modifies 50S subunits more specifically and efficiently than free 23S rRNA.

Pseudouridine modifications in helix 69 (H69) of 23S ribosomal RNA are highly conserved among all organisms. H69 associates with helix 44 of 16S rRNA to form bridge B2a, which plays a vital role in bridging the two ribosomal subunits and stabilizing the ribosome. The three pseudouridines in H69 were shown earlier to play an important role in 50S subunit assembly and in its association with the 30S subunit. In Escherichia coli, these three modifications are made by the pseudouridine synthase, RluD. Previous work showed that RluD is required for normal ribosomal assembly and function, and that it is the only pseudouridine synthase required for normal growth in E. coli. Here, we show that RluD is far more efficient in modifying H69 in structured 50S subunits, compared to free or synthetic 23S rRNA. Based on this observation, we suggest that pseudouridine modifications in H69 are made late in the assembly of 23S rRNA into mature 50S subunits. This is the first reported observation of a pseudouridine synthase being able to modify a highly structured ribonucleoprotein particle, and it may be an important late step in the maturation of 50S ribosomal subunits.