Classifying RNA-binding proteins based on electrostatic properties

Protein structure can provide new insight into the biological function of a protein and can enable the design of better experiments to learn its biological roles. Moreover, deciphering the interactions of a protein with other molecules can contribute to the understanding of the protein's function within cellular processes. In this study, we apply a machine learning approach for classifying RNA-binding proteins based on their three-dimensional structures. The method is based on characterizing unique properties of electrostatic patches on the protein surface. Using an ensemble of general protein features and specific properties extracted from the electrostatic patches, we have trained a support vector machine (SVM) to distinguish RNA-binding proteins from other positively charged proteins that do not bind nucleic acids. Specifically, the method was applied on proteins possessing the RNA recognition motif (RRM) and successfully classified RNA-binding proteins from RRM domains involved in protein-protein interactions. Overall the method achieves 88% accuracy in classifying RNA-binding proteins, yet it cannot distinguish RNA from DNA binding proteins. Nevertheless, by applying a multiclass SVM approach we were able to classify the RNA-binding proteins based on their RNA targets, specifically, whether they bind a ribosomal RNA (rRNA), a transfer RNA (tRNA), or messenger RNA (mRNA). Finally, we present here an innovative approach that does not rely on sequence or structural homology and could be applied to identify novel RNA-binding proteins with unique folds and/or binding motifs.     

RNase-mediated protein footprint sequencing reveals protein-binding sites throughout the human transcriptome

Background

Sequence specific RNA binding proteins are important regulators of gene expression. Several related crosslinking-based, high-throughput sequencing methods, including PAR-CLIP, have recently been developed to determine direct binding sites of global protein-RNA interactions. However, no studies have quantitatively addressed the contribution of background binding to datasets produced by these methods.

Results

We measured non-specific RNA background in PAR-CLIP data, demonstrating that covalently crosslinked background binding is common, reproducible and apparently universal among laboratories. We show that quantitative determination of background is essential for identifying targets of most RNA-binding proteins and can substantially improve motif analysis. We also demonstrate that by applying background correction to an RNA binding protein of unknown binding specificity, Caprin1, we can identify a previously unrecognized RNA recognition element not otherwise apparent in a PAR-CLIP study.

Conclusions

Empirical background measurements of global RNA-protein crosslinking are a necessary addendum to other experimental controls, such as performing replicates, because covalently crosslinked background signals are reproducible and otherwise unavoidable. Recognizing and quantifying the contribution of background extends the utility of PAR-CLIP and can improve mechanistic understanding of protein-RNA specificity, protein-RNA affinity and protein-RNA association dynamics.     

cAMP-independent protein kinase from amphibian lens: identification, organ distribution and substrates of phosphorylation

Eye lens extracts of the frog Rana temporaria contain a cAMP-independent protein kinase which is quantitatively adsorbed on immobilized RNA at physiological salt concentrations. The enzyme activity is maximal in the lenticular cortex, medium in the epithelium and minimal in the lens nuclei. Crude preparations of RNA-binding protein kinase from the epithelium, cortex and nuclei of the eye lens were prepared by affinity chromatography on poly(U)-Sepharose. It was found that these preparations contain no active forms of phosphatases, ATPases or proteases which may interfere with the results of phosphorylation experiments on exogenous and endogenous substrates. The protein kinase under study catalyzes the binding of phosphate groups to threonine and serine residues in casein molecules, does not phosphorylate histones and utilizes GTP alongside with ATP as phosphate donors. Heparin and RNA used at low concentrations inhibit the protein kinase activity. The data obtained allow the identification of lenticular RNA-binding protein kinase(s) as a casein kinase type II. It was shown that incubation of RNA-binding proteins from epithelium and lenticular cortex with [gamma-32P]ATP results in the label incorporation into six to seven polypeptide chains with Mr of 27-130 kDa. Poly(U) and heparin inhibit the self-phosphorylation reaction, cAMP has no stimulating effect on this process, while Ca2+ ions inhibit the self-phosphorylation of RNA-binding proteins.     

Protein profiling with cleavable isotope-coded affinity tag (cICAT) reagents: the yeast salinity stress response

Protein expression profiles in yeast cells, in response to salinity stress, were determined using the cleavable isotope-coded affinity tag (cICAT) labeling strategy. The analysis included separation of the mixed protein samples by SDS-PAGE, followed by excision of the entire gel lane, and division of the lane into 14 gel regions. Regions were subjected to in-gel digestion, biotin affinity chromatography, and analysis by nano-scale microcapillary liquid chromatography coupled to tandem mass spectrometry. The novel (13)C-labeled ICAT reagents have identical elution profiles for labeled peptide pairs and broadly spread the distribution of labeled peptides during reversed-phase chromatography. A total of 560 proteins were identified and quantified, with 51 displaying more than 2-fold expression differences. In addition to some known proteins involved in salt stress, four RNA-binding proteins were found to be up-regulated by high salinity, suggesting that selective RNA export from the nucleus is important for the salt-stress response. Some proteins involved in amino acid synthesis, which have been observed to be up-regulated by amino acid starvation, were also found to increase their abundance on salt stress. These results indicate that salt stress and amino acid starvation cause overlapping cellular responses and are likely to be physiologically linked.     

Epitranscriptomic technologies and analyses

RNA can interact with RNA-binding proteins(RBPs), mRNA, or other non-coding RNAs(ncRNAs) to form complex regulatory networks. High-throughput CLIP-seq, degradome-seq, and RNA-RNA interactome sequencing methods represent powerful approaches to identify biologically relevant ncRNA-target and protein-ncRNA interactions. However, assigning ncRNAs to their regulatory target genes or interacting RNA-binding proteins(RBPs) remains technically challenging. Chemical modifications to mRNA also play important roles in regulating gene expression. Investigation of the functional roles of these modifications relies highly on the detection methods used. RNA structure is also critical at nearly every step of the RNA life cycle. In this review, we summarize recent advances and limitations in CLIP technologies and discuss the computational challenges of and bioinformatics tools used for decoding the functions and regulatory networks of ncRNAs. We also summarize methods used to detect RNA modifications and to probe RNA structure.     

Visinin-like protein (VILIP) is a neuron-specific calcium-dependent double-stranded RNA-binding protein.

Double-stranded RNA-binding proteins function in regulating the stability, translation, and localization of specific mRNAs. In this study, we have demonstrated that the neuron-specific, calcium-binding protein, visinin-like protein (VILIP) contains one double-stranded RNA-binding domain, a protein motif conserved among many double-stranded RNA-binding proteins. We showed that VILIP can specifically bind double-stranded RNA, and this interaction specifically requires the presence of calcium. Mobility shift studies indicated that VILIP binds double-stranded RNA as a single protein-RNA complex with an apparent equilibrium dissociation constant of 9.0 x 10(-6) M. To our knowledge, VILIP is the first double-stranded RNA-binding protein shown to be calcium-dependent. Furthermore, VILIP specifically binds the 3'-untranslated region of the neurotrophin receptor, trkB, an mRNA localized to hippocampal dendrites in an activity-dependent manner. Given that VILIP is also expressed in the hippocampus, these data suggest that VILIP may employ a novel, calcium-dependent mechanism to regulate its binding to important localized mRNAs in the central nervous system.     

Characterization of RNA-binding properties of three types of RNA-binding proteins in Anabaena sp. PPC 7120.

The Rbp proteins in cyanobacteria are RNA-binding proteins with a single RNA recognition motif or RRM. A comprehensive assembly of genomic data suggests that there are two major classes of Rbp proteins (classes I and II) that diverged before the diversification of cyanobacteria. Class I proteins are further classified into two types with or without a C-terminal glycine-rich domain. The results of selection from a random RNA pool suggest that RbpA1 (class I) has affinity to C-rich and G-rich sequences. In vitro RNA binding assay with homopolymers indicated that class II protein has low affinity to poly(G) in contrast with class I proteins. Site-specific mutagenesis analysis of the RRM in RbpA1 showed that the aromatic residues Tyr4 or Phe46 are important in RNA binding as well as maintenance of secondary structure. We also tested various truncated proteins lacking the C-terminal domain as well as point mutants. Most of these proteins exhibited decreased affinity to RNA. Circular dichroism analysis as well as chromatographic analysis showed that Tyr4 and Phe46 are also important in maintaining the structure of RbpA1 protein. The C-terminal glycine-rich domain itself does not contribute much to the RNA-binding, but Arg83 which is located close to the C-terminal end of RRM is important in the RNA-binding.     

Identification and characterization of proteins binding A + U-rich elements

A + U-Rich elements (AREs) have been extensively investigated as cis-acting determinants of rapid mRNA turnover. Recently, a number of RNA-binding proteins interacting with AREs have been described. This article presents strategies and techniques used by our laboratory to identify and characterize a family of ARE-binding proteins collectively termed AUF1. However, these techniques may be applied to the study of any protein displaying sequence-specific RNA binding activity. The techniques described here include the purification of native AUF1 from cultured cells as well as the preparation of recombinant AUF1 proteins using a bacterial expression system. Analyses of RNA-protein interactions are also described, including the use of gel mobility shift assays with synthetic RNA probes to monitor specific RNA binding activity in cell extracts or with recombinant proteins. Variations of this technique are also described to evaluate the RNA binding affinity of recombinant proteins and the use of specific RNA competitors to assess RNA determinants of protein binding specificity. Other techniques presented include the identification of specific proteins in RNA:protein complexes using antibody supershifts and the estimation of molecular weights of RNA-binding proteins by UV crosslinking. Results of individual experiments are presented as examples of some techniques. Throughout the article, suggestions are included to avoid commonly encountered problems and to assist in the optimization of these techniques for the study of other RNA-binding proteins.     

Partner-specific prediction of RNA-binding residues in proteins: A critical assessment

RNA-protein interactions play essential roles in regulating gene expression. While some RNA-protein interactions are “specific”, that is, the RNA-binding proteins preferentially bind to particular RNA sequence or structural motifs, others are “non-RNA specific.” Deciphering the protein-RNA recognition code is essential for comprehending the functional implications of these interactions and for developing new therapies for many diseases. Because of the high cost of experimental determination of protein-RNA interfaces, there is a need for computational methods to identify RNA-binding residues in proteins. While most of the existing computational methods for predicting RNA-binding residues in RNA-binding proteins are oblivious to the characteristics of the partner RNA, there is growing interest in methods for partner-specific prediction of RNA binding sites in proteins. In this work, we assess the performance of two recently published partner-specific protein-RNA interface prediction tools, PS-PRIP, and PRIdictor, along with our own new tools. Specifically, we introduce a novel metric, RNA-specificity metric (RSM), for quantifying the RNA-specificity of the RNA binding residues predicted by such tools. Our results show that the RNA-binding residues predicted by previously published methods are oblivious to the characteristics of the putative RNA binding partner. Moreover, when evaluated using partner-agnostic metrics, RNA partner-specific methods are outperformed by the state-of-the-art partner-agnostic methods. We conjecture that either (a) the protein-RNA complexes in PDB are not representative of the protein-RNA interactions in nature, or (b) the current methods for partner-specific prediction of RNA-binding residues in proteins fail to account for the differences in RNA partner-specific versus partner-agnostic protein-RNA interactions, or both.     

Probing the mechanisms of DEAD-box proteins as general RNA chaperones: the C-terminal domain of CYT-19 mediates general recognition of RNA

The DEAD-box protein CYT-19 functions in the folding of several group I introns in vivo and a diverse set of group I and group II RNAs in vitro. Recent work using the Tetrahymena group I ribozyme demonstrated that CYT-19 possesses a second RNA-binding site, distinct from the unwinding active site, which enhances unwinding activity by binding nonspecifically to the adjacent RNA structure. Here, we probe the region of CYT-19 responsible for that binding by constructing a C-terminal truncation variant that lacks 49 amino acids and terminates at a domain boundary, as defined by limited proteolysis. This truncated protein unwinds a six-base-pair duplex, formed between the oligonucleotide substrate of the Tetrahymena ribozyme and an oligonucleotide corresponding to the internal guide sequence of the ribozyme, with near-wild-type efficiency. However, the truncated protein is activated much less than the wild-type protein when the duplex is covalently linked to the ribozyme or single-stranded or double-stranded extensions. Thus, the active site for RNA unwinding remains functional in the truncated CYT-19, but the site that binds the adjacent RNA structure has been compromised. Equilibrium binding experiments confirmed that the truncated protein binds RNA less tightly than the wild-type protein. RNA binding by the compromised site is important for chaperone activity, because the truncated protein is less active in facilitating the folding of a group I intron that requires CYT-19 in vivo. The deleted region contains arginine-rich sequences, as found in other RNA-binding proteins, and may function by tethering CYT-19 to structured RNAs, so that it can efficiently disrupt exposed, non-native structural elements, allowing them to refold. Many other DExD/H-box proteins also contain arginine-rich ancillary domains, and some of these domains may function similarly as nonspecific RNA-binding elements that enhance general RNA chaperone activity.     

Identification of Sam68 arginine glycine-rich sequences capable of conferring nonspecific RNA binding to the GSG domain

Sam68 is an RNA-binding protein that contains a heterogeneous nuclear ribonucleoprotein K homology domain embedded in a larger RNA binding domain called the GSG (GRP33, Sam68, GLD-1) domain. This family of proteins is often referred to as the STAR (signal transduction and activators of RNA metabolism) proteins. It is not known whether Sam68 is a general nonspecific RNA-binding protein or whether it recognizes specific response elements in mRNAs with high affinity. Sam68 has been shown to bind homopolymeric RNA and a synthetic RNA sequence called G8-5 that has a core UAAA motif. Here we performed a structure function analysis of Sam68 and identified two arginine glycine (RG)-rich regions that confer nonspecific RNA binding to the Sam68 GSG domain. In addition, by using chimeric proteins between Sam68 and QKI-7, we demonstrated that one of the Sam68 RG-rich sequences of 26 amino acids was sufficient to confer homopolymeric RNA binding to the GSG domain of QKI-7, another STAR protein. Furthermore, that minimal sequence can also give QKI-7 the ability (as Sam68) to functionally substitute for HIV-1 REV to facilitate the nuclear export of RNAs. Our studies suggest that neighboring RG-rich sequences may impose nonspecific RNA binding to GSG domains. Because the Sam68 RNA binding activity is negatively regulated by tyrosine phosphorylation, our data lead us to propose that Sam68 might be a specific RNA-binding protein when tyrosine phosphorylated.     

RIP-Chip: the isolation and identification of mRNAs, microRNAs and protein components of ribonucleoprotein complexes from cell extracts

RNA targets of multitargeted RNA-binding proteins (RBPs) can be studied by various methods including mobility shift assays, iterative in vitro selection techniques and computational approaches. These techniques, however, cannot be used to identify the cellular context within which mRNAs associate, nor can they be used to elucidate the dynamic composition of RNAs in ribonucleoprotein (RNP) complexes in response to physiological stimuli. But by combining biochemical and genomics procedures to isolate and identify RNAs associated with RNA-binding proteins, information regarding RNA-protein and RNA-RNA interactions can be examined more directly within a cellular context. Several protocols--including the yeast three-hybrid system and immunoprecipitations that use physical or chemical cross-linking--have been developed to address this issue. Cross-linking procedures in general, however, are limited by inefficiency and sequence biases. The approach outlined here, termed RNP immunoprecipitation-microarray (RIP-Chip), allows the identification of discrete subsets of RNAs associated with multi-targeted RNA-binding proteins and provides information regarding changes in the intracellular composition of mRNPs in response to physical, chemical or developmental inducements of living systems. Thus, RIP-Chip can be used to identify subsets of RNAs that have related functions and are potentially co-regulated, as well as proteins that are associated with them in RNP complexes. Using RIP-Chip, the identification and/or quantification of RNAs in RNP complexes can be accomplished within a few hours or days depending on the RNA detection method used.     

Intrinsically disordered electronegative clusters improve stability and binding specificity of RNA-binding proteins

RNA-binding proteins play crucial roles in various cellular functions and contain abundant disordered protein regions. The disordered regions in RNA-binding proteins are rich in repetitive sequences, such as poly-K/R, poly-N/Q, poly-A, and poly-G residues. Our bioinformatic analysis identified a largely neglected repetitive sequence family we define as electronegative clusters (ENCs) that contain acidic residues and/or phosphorylation sites. The abundance and length of ENCs exceed other known repetitive sequences. Despite their abundance, the functions of ENCs in RNA-binding proteins are still elusive. To investigate the impacts of ENCs on protein stability, RNA-binding affinity, and specificity, we selected one RNA-binding protein, the ribosomal biogenesis factor 15 (Nop15), as a model. We found that the Nop15 ENC increases protein stability and inhibits nonspecific RNA binding, but minimally interferes with specific RNA binding. To investigate the effect of ENCs on sequence specificity of RNA binding, we grafted an ENC to another RNA-binding protein, Ser/Arg-rich splicing factor 3. Using RNA Bind-n-Seq, we found that the engineered ENC inhibits disparate RNA motifs differently, instead of weakening all RNA motifs to the same extent. The motif site directly involved in electrostatic interaction is more susceptible to the ENC inhibition. These results suggest that one of functions of ENCs is to regulate RNA binding via electrostatic interaction. This is consistent with our finding that ENCs are also overrepresented in DNA-binding proteins, whereas underrepresented in halophiles, in which nonspecific nucleic acid binding is inhibited by high concentrations of salts.