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

Deaminase-based RNA recording enables high throughput mutational profiling of protein-RNA interactions

Protein-RNA interactions govern nearly every aspect of RNA metabolism and are frequently dysregulated in disease. While individual protein residues and RNA nucleotides critical for these interactions have been characterized, scalable methods that jointly map protein- and RNA-level determinants remain limited. RNA deaminase fusions have emerged as a powerful strategy to identify transcriptome-wide targets of RNA-binding proteins by converting binding events into site-specific nucleotide edits. Here, we demonstrate that this 'RNA recording' approach enables high-throughput mutational scanning of protein-RNA interfaces. Using the λN-boxB system as a model, we show that editing by a fused TadA adenosine deaminase directly correlates with binding affinity between protein and RNA variants in vitro . Systematic variation of RNA sequence context reveals a strong bias for editing at UA dinucleotides by the engineered TadA8.20, mirroring wild-type TadA preferences. We further demonstrate that stepwise recruitment of the deaminase using nanobody and protein A/G fusions maintains both sequence and binding specificity. Stable expression of the TadA fusion in human cells reproduces in vitro editing patterns across a library of RNA variants. Finally, comprehensive single amino acid mutagenesis of λN in human cells reveals critical residues mediating RNA binding. Together, our results establish RNA recording as a versatile and scalable tool for dissecting protein-RNA interactions at nucleotide and residue resolution, both in vitro and in cells.

Unraveling the Bivalent and Rapid Interactions Between a Multivalent RNA Recognition Motif and RNA: A Kinetic Approach

The kinetics of the interaction between Musashi-1 (MSI1) and RNA have been characterized using surface plasmon resonance biosensor analysis. Truncated variants of human MSI1 encompassing the two homologous RNA recognition motifs (RRM1 and RRM2) in tandem (aa 1-200), and the two RRMs in isolation (aa 1-103 and aa 104-200, respectively) were produced. The proteins were injected over sensor surfaces with immobilized RNA, varying in sequence and length, and with one or two RRM binding motifs. The interactions of the individual RRMs with all RNA variants were well described by a 1:1 interaction model. The interaction between the MSI1 variant encompassing both RRM motifs was bivalent and rapid for all RNA variants. Due to difficulties in fitting this complex data using standard procedures, we devised a new method to quantify the interactions. It revealed that two RRMs in tandem resulted in a significantly longer residence time than a single RRM. It also showed that RNA with double UAG binding motifs and potential hairpin structures forms less stable bivalent complexes with MSI1 than the single UAG motif containing linear RNA. Substituting the UAG binding motif with a CAG sequence resulted in a reduction of the affinity of the individual RRMs, but for MSI1, this reduction was strongly enhanced, demonstrating the importance of bivalency for specificity. This study has provided new insights into the interaction between MSI1 and RNA and an understanding of how individual domains contribute to the overall interaction. It provides an explanation for why many RNA-binding proteins contain dual RRMs.

Causes, functions, and therapeutic possibilities of RNA secondary structure ensembles and alternative states

RNA secondary structure plays essential roles in encoding RNA regulatory fate and function. Most RNAs populate ensembles of alternatively paired states and are continually unfolded and refolded by cellular processes. Measuring these structural ensembles and their contributions to cellular function has traditionally posed major challenges, but new methods and conceptual frameworks are beginning to fill this void. In this review, we provide a mechanism- and function-centric compendium of the roles of RNA secondary structural ensembles and minority states in regulating the RNA life cycle, from transcription to degradation. We further explore how dysregulation of RNA structural ensembles contributes to human disease and discuss the potential of drugging alternative RNA states to therapeutically modulate RNA activity. The emerging paradigm of RNA structural ensembles as central to RNA function provides a foundation for a deeper understanding of RNA biology and new therapeutic possibilities.

Skipper analysis of eCLIP datasets enables sensitive detection of constrained translation factor binding sites

Technology for crosslinking and immunoprecipitation (CLIP) followed by sequencing (CLIP-seq) has identified the transcriptomic targets of hundreds of RNA-binding proteins in cells. To increase the power of existing and future CLIP-seq datasets, we introduce Skipper, an end-to-end workflow that converts unprocessed reads into annotated binding sites using an improved statistical framework. Compared with existing methods, Skipper on average calls 210%-320% more transcriptomic binding sites and sometimes >1,000% more sites, providing deeper insight into post-transcriptional gene regulation. Skipper also calls binding to annotated repetitive elements and identifies bound elements for 99% of enhanced CLIP experiments. We perform nine translation factor enhanced CLIPs and apply Skipper to learn determinants of translation factor occupancy, including transcript region, sequence, and subcellular localization. Furthermore, we observe depletion of genetic variation in occupied sites and nominate transcripts subject to selective constraint because of translation factor occupancy. Skipper offers fast, easy, customizable, and state-of-the-art analysis of CLIP-seq data.

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

Cooperation and competition by RNA-binding proteins in cancer

Post-transcriptional regulation of gene expression plays a major role in determining the cellular proteome in health and disease. Post-transcriptional control mechanisms are disrupted in many cancers, contributing to multiple processes of tumorigenesis. RNA-binding proteins (RBPs), the main post-transcriptional regulators, often show altered expression and activity in cancer cells. Dysregulation of RBPs contributes to many cancer phenotypes, functioning in complex regulatory networks with other cellular players such as non-coding RNAs, signaling mediators and transcription factors to alter the expression of oncogenes and tumor suppressor genes. RBPs often function combinatorially, based on their binding to target sequences/structures on shared mRNA targets, to regulate the expression of cancer-related genes. This gives rise to cooperativity and competition between RBPs in mRNA binding and resultant functional outcomes in post-transcriptional processes such as mRNA splicing, stability, export and translation. Cooperation and competition is also observed in the case of interaction of RBPs and microRNAs with mRNA targets. RNA structural change is a common mechanism mediating the cooperative/competitive interplay between RBPs and between RBPs and microRNAs. RNA modifications, leading to changes in RNA structure, add a new dimension to cooperative/competitive binding of RBPs to mRNAs, further expanding the RBP regulatory landscape. Therefore, cooperative/competitive interplay between RBPs is a major determinant of the RBP interactome and post-transcriptional regulation of gene expression in cancer cells.

RNA-binding proteins La and HuR cooperatively modulate translation repression of PDCD4 mRNA

Posttranscriptional regulation of gene expression plays a critical role in controlling the inflammatory response. An uncontrolled inflammatory response results in chronic inflammation, often leading to tumorigenesis. Programmed cell death 4 (PDCD4) is a proinflammatory tumor-suppressor gene which helps to prevent the transition from chronic inflammation to cancer. PDCD4 mRNA translation is regulated by an interplay between the oncogenic microRNA miR-21 and the RNA-binding protein (RBP) human antigen R (HuR) in response to lipopolysaccharide stimulation, but the role of other regulatory factors remains unknown. Here, we report that the RBP lupus antigen (La) interacts with the 3'-untranslated region of PDCD4 mRNA and prevents miR-21-mediated translation repression. While lipopolysaccharide causes nuclear-cytoplasmic translocation of HuR, it enhances cellular La expression. Remarkably, La and HuR were found to bind cooperatively to the PDCD4 mRNA and mitigate miR-21-mediated translation repression. The cooperative action of La and HuR reduced cell proliferation and enhanced apoptosis, reversing the pro-oncogenic function of miR-21. Together, these observations demonstrate a cooperative interplay between two RBPs, triggered differentially by the same stimulus, which exerts a synergistic effect on PDCD4 expression and thereby helps maintain a balance between inflammation and tumorigenesis.

Transcriptome-wide analysis reveals spatial correlation between N6-methyladenosine and binding sites of microRNAs and RNA-binding proteins

N6-methyladenosine (m⁶A), the most prevalent epitranscriptomic modification in eukaryotes, is enriched in 3'-untranslated regions (3'UTRs) of mRNAs. As 3'UTRs are major binding sites of RNA-binding proteins (RBPs) and microRNAs (miRNAs), m⁶A-dependent local RNA structure change may alter the accessibility of RBPs and miRNAs to their target sites and regulate mRNA function. Using a human transcriptome-wide computational analysis to investigate the relation between m⁶A, RBPs and miRNAs, we find a strong positive correlation between number of m⁶A sites, miRNAs and RBPs binding to mRNAs, suggesting m⁶A-modified mRNAs are more targeted by miRNAs and RBPs. Moreover, m⁶A sites are located proximally to miRNA target sites and binding sites of multiple RBPs. Further, miRNA target sites and RBP-binding sites located close to each other are also located proximally to m⁶A. This study indicates three-way interplay between m⁶A, microRNA and RBP binding, suggesting the influence of mRNA modifications on the miRNA and RBP interactomes.

An evolutionarily conserved RNA structure in the functional core of the lincRNA Cyrano

The wide prevalence and regulated expression of long noncoding RNAs (lncRNAs) highlight their functional roles, but the molecular basis for their activities and structure-function relationships remains to be investigated, with few exceptions. Among the relatively few lncRNAs conserved over significant evolutionary distances is the long intergenic noncoding RNA (lincRNA) Cyrano (orthologous to human OIP5-AS1), which contains a region of 300 highly conserved nucleotides within tetrapods, which in turn contains a functional stretch of 26 nt of deep conservation. This region binds to and facilitates the degradation of the microRNA miR-7, a short ncRNA with multiple cellular functions, including modulation of oncogenic expression. We probed the secondary structure of Cyrano in vitro and in cells using chemical and enzymatic probing, and validated the results using comparative sequence analysis. At the center of the functional core of Cyrano is a cloverleaf structure maintained over the >400 million years of divergent evolution that separates fish and primates. This strikingly conserved motif provides interaction sites for several RNA-binding proteins and masks a conserved recognition site for miR-7. Conservation in this region strongly suggests that the function of Cyrano depends on the formation of this RNA structure, which could modulate the rate and efficiency of degradation of miR-7.

Competitive endogenous network of lncRNA, miRNA, and mRNA in the chemoresistance of gastrointestinal tract adenocarcinomas

Chemotherapy is one of the main therapeutic strategies used for gastrointestinal tract adenocarcinomas (GTAs), but resistance to anticancer drugs is a substantial obstacle in successful chemotherapy. Accumulating evidence shows that non-coding RNAs, especially long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), can affect the drug resistance of tumor cells by forming a ceRNA regulatory network with mRNAs. The efficiency of the competing endogenous RNAs (ceRNAs) network can be affected by the number and integrality of miRNA recognition elements (MREs). Dynamic factors such as RNA editing, alternative splicing, single nucleotide polymorphism (SNP), RNA-binding proteins and RNA secondary structure can influence the MRE activity, which may in turn be involved in the regulation of chemoresistance-associated ceRNA network by prospective approaches. Besides activities in a single tumor cell, the components of the tumor micoenvironment (TME) also affect the ceRNA network by regulating the expression of non-coding RNA directly or indirectly. The alternation of the ceRNA network often has an impact on the malignant phenotype of tumor including chemoresistance. In this review, we focused on how MRE-associated dynamic factors and components of TME affected the ceRNA network and speculated the potential association of ceRNA network with chemoresistance. We also summarized the ceRNA network of lncRNAs, miRNAs, and mRNAs which efficiently triggers chemoresistance in the specific types of GTAs and analyzed the role of each RNA as a "promoter" or "suppressor" of chemoresistance.

High throughput sequencing of in vitro selections of mRNA-displayed peptides: data analysis and applications

In vitro selection using mRNA display is currently a widely used method to isolate functional peptides with desired properties. The analysis of high throughput sequencing (HTS) data from in vitro evolution experiments has proven to be a powerful technique but only recently has it been applied to mRNA display selections. In this Perspective, we introduce aspects of mRNA display and HTS that may be of interest to physical chemists. We highlight the potential of HTS to analyze in vitro selections of peptides and review recent advances in the application of HTS analysis to mRNA display experiments. We discuss some possible issues involved with HTS analysis and summarize some strategies to alleviate them. Finally, the potential for future impact of advancing HTS analysis on mRNA display experiments is discussed.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

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.

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.