Integrated analysis of RNA-binding protein complexes using in vitro selection and high-throughput sequencing and sequence specificity landscapes (SEQRS)

The RNA-Binding Protein HuR Is Integral to the Function of Nociceptors in Mice and Humans

HuR is an RNA-binding protein implicated in RNA processing, stability, and translation. Previously, we examined protein synthesis in dorsal root ganglion (DRG) neurons treated with inflammatory mediators using ribosome profiling. We found that the HuR consensus binding element was enriched in transcripts with elevated translation. HuR is expressed in the soma of nociceptors and their axons. Pharmacologic inhibition of HuR with the small molecule CMLD-2 reduced the activity of mouse and human sensory neurons. Peripheral administration of CMLD-2 in the paw or genetic elimination of HuR from sensory neurons diminished behavioral responses associated with NGF- and IL-6-induced allodynia in male and female mice. Genetic disruption of HuR altered the proximity of mRNA decay factors near a key neurotrophic factor (TrkA). Collectively, the data suggest that HuR is required for local control of mRNA stability and reveals a new biological function for a broadly conserved post-transcriptional regulatory factor.SIGNIFICANCE STATEMENT Nociceptors undergo long-lived changes in excitability, which may contribute to chronic pain. Noxious cues that promote pain lead to rapid induction of protein synthesis. The underlying mechanisms that confer specificity to mRNA control in nociceptors are unclear. Here, we identify a conserved RNA-binding protein called HuR as a key regulatory factor in sensory neurons. Using a combination of genetics and pharmacology, we demonstrate that HuR is required for signaling in nociceptors. In doing so, we report an important mechanism of mRNA control in sensory neurons that ensures appropriate nociceptive responses to inflammatory mediators.

PRIESSTESS: interpretable, high-performing models of the sequence and structure preferences of RNA-binding proteins

Modelling both primary sequence and secondary structure preferences for RNA binding proteins (RBPs) remains an ongoing challenge. Current models use varied RNA structure representations and can be difficult to interpret and evaluate. To address these issues, we present a universal RNA motif-finding/scanning strategy, termed PRIESSTESS (Predictive RBP-RNA InterpretablE Sequence-Structure moTif regrESSion), that can be applied to diverse RNA binding datasets. PRIESSTESS identifies dozens of enriched RNA sequence and/or structure motifs that are subsequently reduced to a set of core motifs by logistic regression with LASSO regularization. Importantly, these core motifs are easily visualized and interpreted, and provide a measure of RBP secondary structure specificity. We used PRIESSTESS to interrogate new HTR-SELEX data for 23 RBPs with diverse RNA binding modes and captured known primary sequence and secondary structure preferences for each. Moreover, when applying PRIESSTESS to 144 RBPs across 202 RNA binding datasets, 75% showed an RNA secondary structure preference but only 10% had a preference besides unpaired bases, suggesting that most RBPs simply recognize the accessibility of primary sequences.

Methods for the directed evolution of biomolecular interactions

Noncovalent interactions between biomolecules such as proteins and nucleic acids coordinate all cellular processes through changes in proximity. Tools that perturb these interactions are and will continue to be highly valuable for basic and translational scientific endeavors. By taking cues from natural systems, such as the adaptive immune system, we can design directed evolution platforms that can generate proteins that bind to biomolecules of interest. In recent years, the platforms used to direct the evolution of biomolecular binders have greatly expanded the range of types of interactions one can evolve. Herein, we review recent advances in methods to evolve protein-protein, protein-RNA, and protein-DNA interactions.

Expanding the binding specificity for RNA recognition by a PUF domain

The ability to design a protein to bind specifically to a target RNA enables numerous applications, with the modular architecture of the PUF domain lending itself to new RNA-binding specificities. For each repeat of the Pumilio-1 PUF domain, we generate a library that contains the 8,000 possible combinations of amino acid substitutions at residues critical for RNA contact. We carry out yeast three-hybrid selections with each library against the RNA recognition sequence for Pumilio-1, with any possible base present at the position recognized by the randomized repeat. We use sequencing to score the binding of each variant, identifying many variants with highly repeat-specific interactions. From these data, we generate an RNA binding code specific to each repeat and base. We use this code to design PUF domains against 16 RNAs, and find that some of these domains recognize RNAs with two, three or four changes from the wild type sequence.

Principles of mRNA control by human PUM proteins elucidated from multimodal experiments and integrative data analysis

The human PUF-family proteins, PUM1 and PUM2, posttranscriptionally regulate gene expression by binding to a PUM recognition element (PRE) in the 3'-UTR of target mRNAs. Hundreds of PUM1/2 targets have been identified from changes in steady-state RNA levels; however, prior studies could not differentiate between the contributions of changes in transcription and RNA decay rates. We applied metabolic labeling to measure changes in RNA turnover in response to depletion of PUM1/2, showing that human PUM proteins regulate expression almost exclusively by changing RNA stability. We also applied an in vitro selection workflow to precisely identify the binding preferences of PUM1 and PUM2. By integrating our results with prior knowledge, we developed a "rulebook" of key contextual features that differentiate functional versus nonfunctional PREs, allowing us to train machine learning models that accurately predict the functional regulation of RNA targets by the human PUM proteins.

High throughput approaches to study RNA-protein interactions in vitro

To understand the regulation of gene expression it is critical to determine how proteins interact with and discriminate between different RNAs. In this review, we discuss experimental techniques that utilize high throughput approaches to characterize the interactions of proteins with large numbers of RNAs in vitro. We describe the underlying principles for the main methods, briefly discuss their scope and limitations, and outline how insight from the techniques contributes to our understanding of specificity for RNA-protein interactions.

Molecular entrapment by RNA: an emerging tool for disrupting protein-RNA interactions in vivo

mRNA function is controlled by RNA-binding proteins. The specificity of RNA-binding factors for their targets is critical in that it enables all subsequent regulation. Despite widespread recognition of the pervasive role RNA-binding proteins play in development and disease, they remain challenging to target with small molecules. A renaissance in RNA therapeutics has led to the identification of modifications that substantially increase RNA stability. When combined with information regarding specificity, a new class of oligonucleotide mimics has emerged as a means to competitively disrupt the regulation of endogenous substrates. These decoys have been used to inhibit RNA-binding proteins in living animals. Decoys will likely provide new insights into the expansive roles of RNA-binding proteins in biology and disease. Here, we describe examples where they have been used and discuss how they could be applied to new targets.

Unique repression domains of Pumilio utilize deadenylation and decapping factors to accelerate destruction of target mRNAs

Pumilio is an RNA-binding protein that represses a network of mRNAs to control embryogenesis, stem cell fate, fertility and neurological functions in Drosophila. We sought to identify the mechanism of Pumilio-mediated repression and find that it accelerates degradation of target mRNAs, mediated by three N-terminal Repression Domains (RDs), which are unique to Pumilio orthologs. We show that the repressive activities of the Pumilio RDs depend on specific subunits of the Ccr4-Not (CNOT) deadenylase complex. Depletion of Pop2, Not1, Not2, or Not3 subunits alleviates Pumilio RD-mediated repression of protein expression and mRNA decay, whereas depletion of other CNOT components had little or no effect. Moreover, the catalytic activity of Pop2 deadenylase is important for Pumilio RD activity. Further, we show that the Pumilio RDs directly bind to the CNOT complex. We also report that the decapping enzyme, Dcp2, participates in repression by the N-terminus of Pumilio. These results support a model wherein Pumilio utilizes CNOT deadenylase and decapping complexes to accelerate destruction of target mRNAs. Because the N-terminal RDs are conserved in mammalian Pumilio orthologs, the results of this work broadly enhance our understanding of Pumilio function and roles in diseases including cancer, neurodegeneration and epilepsy.

Linking RNA Sequence, Structure, and Function on Massively Parallel High-Throughput Sequencers

High-throughput sequencing methods have revolutionized our ability to catalog the diversity of RNAs and RNA-protein interactions that can exist in our cells. However, the relationship between RNA sequence, structure, and function is enormously complex, demonstrating the need for methods that can provide quantitative thermodynamic and kinetic measurements of macromolecular interaction with RNA, at a scale commensurate with the sequence diversity of RNA. Here, we discuss a class of methods that extend the core functionality of DNA sequencers to enable high-throughput measurements of RNA folding and RNA-protein interactions. Topics discussed include a description of the method and multiple applications to RNA-binding proteins, riboswitch design and engineering, and RNA tertiary structure energetics.

The target specificity of the RNA binding protein Pumilio is determined by distinct co-factors

Puf family proteins are translational regulators essential to a wide range of biological processes, including cell fate specification, stem cell self-renewal, and neural function. Yet, despite being associated with hundreds of RNAs, the underlying mechanisms of Puf target specification remain to be fully elucidated. In Drosophila, Pumilio – a sole Puf family protein – is known to collaborate with cofactors Nanos (Nos) and Brain Tumor (Brat); however, their roles in target specification are not clearly defined. Here, we identify Bag-of-marbles (Bam) as a new Pum cofactor in repression of Mothers against dpp (mad) mRNAs, for which Nos is known to be dispensable. Notably, our data show that Nos (but not Bam) was required for Pum association with hunchback (hb) mRNAs, a well-known target of Pum and Nos. In contrast, Bam (but not Nos) was required for Pum association with mad mRNAs. These findings show for the first time that Pum target specificity is determined not independently but in collaboration with cofactors.

The role of coevolutionary signatures in protein interaction dynamics, complex inference, molecular recognition, and mutational landscapes

Evolution imposes constraints at the interface of interacting biomolecules in order to preserve function or maintain fitness. This pressure may have a direct effect on the sequence composition of interacting biomolecules. As a result, statistical patterns of amino acid or nucleotide covariance that encode for physical and functional interactions are observed in sequences of extant organisms. In recent years, global pairwise models of amino acid and nucleotide coevolution from multiple sequence alignments have been developed and utilized to study molecular interactions in structural biology. In proteins, for which the energy landscape is funneled and minimally frustrated, a direct connection between the physical and sequence space landscapes can be established. Estimating coevolutionary information from sequences of interacting molecules has a broad impact in molecular biology. Applications include the accurate determination of 3D structures of molecular complexes, inference of protein interaction partners, models of protein-protein interaction specificity, the elucidation, and design of protein-nucleic acid recognition as well as the discovery of genome-wide epistatic effects. The current state of the art of coevolutionary analysis includes biomedical applications ranging from mutational landscapes and drug-design to vaccine development.

A crystal structure of a collaborative RNA regulatory complex reveals mechanisms to refine target specificity

In the Caenorhabditis elegans germline, fem-3 Binding Factor (FBF) partners with LST-1 to maintain stem cells. A crystal structure of an FBF-2/LST-1/RNA complex revealed that FBF-2 recognizes a short RNA motif different from the characteristic 9-nt FBF binding element, and compact motif recognition coincided with curvature changes in the FBF-2 scaffold. Previously, we engineered FBF-2 to favor recognition of shorter RNA motifs without curvature change (Bhat et al., 2019). In vitro selection of RNAs bound by FBF-2 suggested sequence specificity in the central region of the compact element. This bias, reflected in the crystal structure, was validated in RNA-binding assays. FBF-2 has the intrinsic ability to bind to this shorter motif. LST-1 weakens FBF-2 binding affinity for short and long motifs, which may increase target selectivity. Our findings highlight the role of FBF scaffold flexibility in RNA recognition and suggest a new mechanism by which protein partners refine target site selection.

Preparation of cooperative RNA recognition complexes for crystallographic structural studies

It is essential that mRNA-binding proteins recognize specific motifs in target mRNAs to control their processing, localization, and expression. Although mRNAs are typically targets of many different regulatory factors, our understanding of how they work together is limited. In some cases, RNA-binding proteins work cooperatively to regulate an mRNA target. A classic example is Drosophila melanogaster Pumilio (Pum) and Nanos (Nos). Pum is a sequence-specific RNA-binding protein. Nos also binds RNA, but interaction with some targets requires Pum to bind first. We recently determined crystal structures of complexes of Pum and Nos with two different target RNA sequences. A crystal structure in complex with the hunchback mRNA element showed how Pum and Nos together can recognize an extended RNA sequence with Nos binding to an A/U-rich sequence 5' of the Pum sequence element. Nos also enables recognition of elements that contain an A/U-rich 5' sequence, but imperfectly match the Pum sequence element. We determined a crystal structure of Pum and Nos in complex with the Cyclin B mRNA element, which demonstrated how Nos clamps the Pum-RNA complex and enables recognition of the imperfect element. Here, we describe methods for expression and purification of stable Pum-Nos-RNA complexes for crystallization, details of the crystallization and structure determination, and guidance on how to analyze protein-RNA structures and evaluate structure-driven hypotheses. We aim to provide tips and guidance that can be applied to other protein-RNA complexes. With hundreds of mRNA-binding proteins identified, combinatorial control is likely to be common, and much work remains to understand them structurally.

Global pairwise RNA interaction landscapes reveal core features of protein recognition

RNA–protein interactions permeate biology. Transcription, translation, and splicing all hinge on the recognition of structured RNA elements by RNA-binding proteins. Models of RNA–protein interactions are generally limited to short linear motifs and structures because of the vast sequence sampling required to access longer elements. Here, we develop an integrated approach that calculates global pairwise interaction scores from in vitro selection and high-throughput sequencing. We examine four RNA-binding proteins of phage, viral, and human origin. Our approach reveals regulatory motifs, discriminates between regulated and non-regulated RNAs within their native genomic context, and correctly predicts the consequence of mutational events on binding activity. We design binding elements that improve binding activity in cells and infer mutational pathways that reveal permissive versus disruptive evolutionary trajectories between regulated motifs. These coupling landscapes are broadly applicable for the discovery and characterization of protein–RNA recognition at single nucleotide resolution.

RNA–protein interactions often depend on the recognition of extended RNA elements but the identification of these motifs is challenging. Here, the authors present a global integrated approach to analyze RNA–protein binding landscapes, mapping extended RNA interaction motifs for four RNA-binding proteins.

Inhibition of Poly(A)-binding protein with a synthetic RNA mimic reduces pain sensitization in mice

Nociceptors rely on cap-dependent translation to rapidly induce protein synthesis in response to pro-inflammatory signals. Comparatively little is known regarding the role of the regulatory factors bound to the 3′ end of mRNA in nociceptor sensitization. Poly(A)-binding protein (PABP) stimulates translation initiation by bridging the Poly(A) tail to the eukaryotic initiation factor 4F complex associated with the mRNA cap. Here, we use unbiased assessment of PABP binding specificity to generate a chemically modified RNA-based competitive inhibitor of PABP. The resulting RNA mimic, which we designated as the Poly(A) SPOT-ON, is more stable than unmodified RNA and binds PABP with high affinity and selectivity in vitro. We show that injection of the Poly(A) SPOT-ON at the site of an injury can attenuate behavioral response to pain. Collectively, these results suggest that PABP is integral for nociceptive plasticity. The general strategy described here provides a broad new source of mechanism-based inhibitors for RNA-binding proteins and is applicable for in vivo studies.

Poly(A)-binding protein (PABP) is an RNA binding protein with translation function. Here, Barragán-Iglesias and colleagues devise an RNA mimic that inhibits PABP activity, and show that inhibitors can reduce animal’s pain response in vivo when injected locally.

The contribution of the C5 protein subunit of Escherichia coli ribonuclease P to specificity for precursor tRNA is modulated by proximal 5' leader sequences

Recognition of RNA by RNA processing enzymes and RNA binding proteins often involves cooperation between multiple subunits. However, the interdependent contributions of RNA and protein subunits to molecular recognition by ribonucleoproteins are relatively unexplored. RNase P is an endonuclease that removes 5' leaders from precursor tRNAs and functions in bacteria as a dimer formed by a catalytic RNA subunit (P RNA) and a protein subunit (C5 in E. coli). The P RNA subunit contacts the tRNA body and proximal 5' leader sequences [N(-1) and N(-2)] while C5 binds distal 5' leader sequences [N(-3) to N(-6)]. To determine whether the contacts formed by P RNA and C5 contribute independently to specificity or exhibit cooperativity or anti-cooperativity, we compared the relative k_cat/K_m values for all possible combinations of the six proximal 5' leader nucleotides (n = 4096) for processing by the E. coli P RNA subunit alone and by the RNase P holoenzyme. We observed that while the P RNA subunit shows specificity for 5' leader nucleotides N(-2) and N(-1), the presence of the C5 protein reduces the contribution of P RNA to specificity, but changes specificity at N(-2) and N(-3). The results reveal that the contribution of C5 protein to RNase P processing is controlled by the identity of N(-2) in the pre-tRNA 5' leader. The data also clearly show that pairing of the 5' leader with the 3' ACCA of tRNA acts as an anti-determinant for RNase P cleavage. Comparative analysis of genomically encoded E. coli tRNAs reveals that both anti-determinants are subject to negative selection in vivo.

Combinatorial control of messenger RNAs by Pumilio, Nanos and Brain Tumor Proteins

Eukaryotes possess a vast array of RNA-binding proteins (RBPs) that affect mRNAs in diverse ways to control protein expression. Combinatorial regulation of mRNAs by RBPs is emerging as the rule. No example illustrates this as vividly as the partnership of 3 Drosophila RBPs, Pumilio, Nanos and Brain Tumor, which have overlapping functions in development, stem cell maintenance and differentiation, fertility and neurologic processes. Here we synthesize 30 y of research with new insights into their molecular functions and mechanisms of action. First, we provide an overview of the key properties of each RBP. Next, we present a detailed analysis of their collaborative regulatory mechanism using a classic example of the developmental morphogen, hunchback, which is spatially and temporally regulated by the trio during embryogenesis. New biochemical, structural and functional analyses provide insights into RNA recognition, cooperativity, and regulatory mechanisms. We integrate these data into a model of combinatorial RNA binding and regulation of translation and mRNA decay. We then use this information, transcriptome wide analyses and bioinformatics predictions to assess the global impact of Pumilio, Nanos and Brain Tumor on gene regulation. Together, the results support pervasive, dynamic post-transcriptional control.