Discovering the interaction propensities of amino acids and nucleotides from protein-RNA complexes

With the availability of many genome sequences, the mining of biological data is attracting much attention, most of it limited to the sequences of macromolecules. Sequence data are easy to analyze as they can be treated as strings of characters, whereas the structure of a macromolecule is much more complex. We developed a set of algorithms to analyze the structures of protein-RNA complexes at the atomic level and used them to analyze protein-RNA interactions using structural data on 51 protein-RNA complexes. The analysis revealed, among other things, that: (1) polar and charged amino acids have a strong tendency to interact with nucleotides, (2) arginine and asparagine tend to hydrogen bond with uracil, and (3) histidine favors uracil in water-mediated bonding with RNA. We analyzed a large set of structural data of protein-RNA complexes involving water-mediated hydrogen bonds as well as direct hydrogen bonds. The interaction patterns discovered from the analysis provide useful information for predicting the structure of RNA that binds proteins, and of proteins that bind RNA.     

Molecular architecture of protein-RNA recognition sites

The molecular architecture of protein-RNA interfaces are analyzed using a non-redundant dataset of 152 protein-RNA complexes. We find that an average protein-RNA interface is smaller than an average protein-DNA interface but larger than an average protein–protein interface. Among the different classes of protein-RNA complexes, interfaces with tRNA are the largest, while the interfaces with the single-stranded RNA are the smallest. Significantly, RNA contributes more to the interface area than its partner protein. Moreover, unlike protein–protein interfaces where the side chain contributes less to the interface area compared to the main chain, the main chain and side chain contributions flipped in protein-RNA interfaces. We find that the protein surface in contact with the RNA in protein-RNA complexes is better packed than that in contact with the DNA in protein-DNA complexes, but loosely packed than that in contact with the protein in protein–protein complexes. Shape complementarity and electrostatic potential are the two major factors that determine the specificity of the protein-RNA interaction. We find that the H-bond density at the protein-RNA interfaces is similar with that of protein-DNA interfaces but higher than the protein–protein interfaces. Unlike protein-DNA interfaces where the deoxyribose has little role in intermolecular H-bonds, due to the presence of an oxygen atom at the 2′ position, the ribose in RNA plays significant role in protein-RNA H-bonds. We find that besides H-bonds, salt bridges and stacking interactions also play significant role in stabilizing protein-nucleic acids interfaces; however, their contribution at the protein–protein interfaces is insignificant.     

PRI-Modeler: extracting RNA structural elements from PDB files of protein-RNA complexes

A complete understanding of protein and RNA structures and their interactions is important for determining the binding sites in protein-RNA complexes. Computational approaches exist for identifying secondary structural elements in proteins from atomic coordinates. However, similar methods have not been developed for RNA, due in part to the very limited structural data so far available. We have developed a set of algorithms for extracting and visualizing secondary and tertiary structures of RNA and for analyzing protein-RNA complexes. These algorithms have been implemented in a web-based program called PRI-Modeler (protein-RNA interaction modeler). Given one or more protein data bank files of protein-RNA complexes, PRI-Modeler analyzes the conformation of the RNA, calculates the hydrogen bond (H bond) and van der Waals interactions between amino acids and nucleotides, extracts secondary and tertiary RNA structure elements, and identifies the patterns of interactions between the proteins and RNAs. This paper presents PRI-Modeler and its application to the hydrogen bond and van der Waals interactions in the most representative set of protein-RNA complexes. The analysis reveals several interesting interaction patterns at various levels. The information provided by PRI-Modeler should prove useful for determining the binding sites in protein-RNA complexes. PRI-Modeler is accessible at http://wilab.inha.ac.kr/primodeler/, and supplementary materials are available in the analysis results section at http://wilab.inha.ac.kr/primodeler/.     

Non-random localization of ribonucleoprotein (RNP) structures within an adenovirus mRNA precursor.

Heterogeneous nuclear protein complexes (hnRNP) containing the precursor RNA from the adenovirus early region 2 were analysed to determine the specificity of protein-RNA interaction. RNA precursor sequences were present in isolated hnRNP complexes and endogenous 30S particles. At least 20-40 bases long fragments were protected when RNase A was used to remove unprotected RNA sequences in hnRNA complexes. Similarly around 40 bases of RNA were protected in 30S particles. These sequences represent discrete regions of the adenovirus genome. Especially sequences complementary to the EcoRI-F fragment encoding the first leader and the major intron for the DNA binding protein (DBP) RNA precursor, were analysed in detail. Tentatively, sequences resistant to RNase A were located in the middle of the intron and at the splice-donor junction of the first leader of the DBP precursor RNA. The same sequences were identified irrespective whether hnRNP complexes or 30S particles were used suggesting that 30S particles originate from hnRNP complexes. A 38.000 dalton protein appears to be in direct contact with RNA sequences complementary to the EcoRI-F fragment.     

Dissecting the protein-RNA interface: the role of protein surface shapes and RNA secondary structures in protein-RNA recognition

Protein-RNA interactions are essential for many biological processes. However, the structural mechanisms underlying these interactions are not fully understood. Here, we analyzed the protein surface shape (dented, intermediate or protruded) and the RNA base pairing properties (paired or unpaired nucleotides) at the interfaces of 91 protein-RNA complexes derived from the Protein Data Bank. Dented protein surfaces prefer unpaired nucleotides to paired ones at the interface, and hydrogen bonds frequently occur between the protein backbone and RNA bases. In contrast, protruded protein surfaces do not show such a preference, rather, electrostatic interactions initiate the formation of hydrogen bonds between positively charged amino acids and RNA phosphate groups. Interestingly, in many protein-RNA complexes that interact via an RNA loop, an aspartic acid is favored at the interface. Moreover, in most of these complexes, nucleotide bases in the RNA loop are flipped out and form hydrogen bonds with the protein, which suggests that aspartic acid is important for RNA loop recognition through a base-flipping process. This study provides fundamental insights into the role of the shape of the protein surface and RNA secondary structures in mediating protein-RNA interactions.     

Computational methods for prediction of protein-RNA interactions

Understanding the molecular mechanism of protein-RNA recognition and complex formation is a major challenge in structural biology. Unfortunately, the experimental determination of protein-RNA complexes by X-ray crystallography and nuclear magnetic resonance spectroscopy (NMR) is tedious and difficult. Alternatively, protein-RNA interactions can be predicted by computational methods. Although less accurate than experimental observations, computational predictions can be sufficiently accurate to prompt functional hypotheses and guide experiments, e.g. to identify individual amino acid or nucleotide residues. In this article we review 10 methods for predicting protein-RNA interactions, seven of which predict RNA-binding sites from protein sequences and three from structures. We also developed a meta-predictor that uses the output of top three sequence-based primary predictors to calculate a consensus prediction, which outperforms all the primary predictors. In order to fully cover the software for predicting protein-RNA interactions, we also describe five methods for protein-RNA docking. The article highlights the strengths and shortcomings of existing methods for the prediction of protein-RNA interactions and provides suggestions for their further development.     

Recognition of GU-rich polyadenylation regulatory elements by human CstF-64 protein

Vertebrate polyadenylation sites are identified by the AAUAAA signal and by GU-rich sequences downstream of the cleavage site. These are recognized by a heterotrimeric protein complex (CstF) through its 64 kDa subunit (CstF-64); the strength of this interaction affects the efficiency of poly(A) site utilization. We present the structure of the RNA-binding domain of CstF-64 containing an RNA recognition motif (RRM) augmented by N- and C-terminal helices. The C-terminal helix unfolds upon RNA binding and extends into the hinge domain where interactions with factors responsible for assembly of the polyadenylation complex occur. We propose that this conformational change initiates assembly. Consecutive Us are required for a strong CstF-GU interaction and we show how UU dinucleotides are recognized. Contacts outside the UU pocket fine tune the protein-RNA interaction and provide different affinities for distinct GU-rich elements. The protein-RNA interface remains mobile, most likely a requirement to bind many GU-rich sequences and yet discriminate against other RNAs. The structural distinction between sequences that form stable and unstable complexes provides an operational distinction between weakly and strongly processed poly(A) sites.     

Leakage induced in Escherichia coli cells by A protein-RNA complexes from bacteriophage f2

Complexes of f2 phage RNA and its A protein, or maturation protein, transfect Escherichia coli cells much better than does protein-free RNA. We used these complexes to introduce the bacteriophage f2 lysis gene into cells. The A protein-RNA complex was found to kill cells, probably by causing them to leak large macromolecules. Previously induced beta-galactosidase leaked from cells treated either with the A protein-RNA complex or with lethal but noninfectious complexes that had been treated with formaldehyde. This observation was consistent with an earlier finding that formaldehyde-treated f2 RNA stimulates the in vitro synthesis of a lysis protein. The complexes did not stimulate the rate of leakage of beta-galactosidase from a streptomycin-resistant mutant known to be lysis defective. On the other hand, the rate of leakage was increased in a double mutant resistant to both streptomycin and rifampin and which is lysed normally by f2 bacteriophage.     

Study of the assembly of vesicular stomatitis virus N protein: role of the P protein

To derive structural information about the vesicular stomatitis virus (VSV) nucleocapsid (N) protein, the N protein and the VSV phosphoprotein (P protein) were expressed together in Escherichia coli. The N and P proteins formed soluble protein complexes of various molar ratios when coexpressed. The major N/P protein complex was composed of 10 molecules of the N protein, 5 molecules of the P protein, and an RNA. A soluble N protein-RNA oligomer free of the P protein was isolated from the N/P protein-RNA complex using conditions of lowered pH. The molecular weight of the N protein-RNA oligomer, 513,879, as determined by analytical ultracentrifugation, showed that it was composed of 10 molecules of the N protein and an RNA of approximately 90 nucleotides. The N protein-RNA oligomer had the appearance of a disk with outer diameter, inner diameter, and thickness of 148 +/- 10 A, 78 +/- 9 A, and 83 +/- 8 A, respectively, as determined by electron microscopy. RNA in the complexes was protected from RNase digestion and was stable at pH 11. This verified that N/P protein complexes expressed in E. coli were competent for encapsidation. In addition to coexpression with the full-length P protein, the N protein was expressed with the C-terminal 72 amino acids of the P protein. This portion of the P protein was sufficient for binding to the N protein, maintaining it in a soluble state, and for assembly of N protein-RNA oligomers. With the results provided in this report, we propose a model for the assembly of an N/P protein-RNA oligomer.     

Retention of conformational flexibility in HIV-1 Rev-RNA complexes

Sequence-specific recognition between HIV-1 Rev and viral RNA mediates the nuclear export of the viral mRNA and thus is important for the viral life cycle. HIV Rev binds to its viral RNA target with high affinity and specificity and also binds to an in vitro selected RNA aptamer that has a significantly different sequence from the viral RNA target with a 6-fold higher affinity than its natural target. The high-resolution structures of HIV Rev Arg-rich motif (ARM) in complexes with the wild-type RNA and the RNA aptamer reveal that, despite the significantly different RNA sequences, the two complexes share similar structural features and the protein-RNA interactions are mediated mostly by the Arg side chains in Rev ARM. To gain further insight into the role of these Arg side chains in the sequence-specific protein-RNA recognition, we have characterized the flexibility of these Arg side chains at the interfaces of the two high-affinity complexes using (15)N R(1), R(2), nuclear Overhauser effect, and chemical-shift anisotropy dipolar cross-correlation relaxation measurements. The ARM peptide contains uniformly (13)C/(15)N-labeled Arg residues, and the RNA samples were unlabeled. Despite the apparently similar roles of Arg side chains in both complexes, most of them display a different dynamic behavior in the context of different RNA molecules, and extensive and highly diverse motions have been observed for all of these side chains that interact with RNA. Most of the differences in side-chain dynamics between the complexes cannot be inferred from the three-dimensional structures. Additionally, more than half of the residues have increased flexibility in the Rev-RNA aptamer complex that has a higher affinity. This study provides new insights into ARM-RNA recognition and indicates that retention of conformational flexibility is likely important in high-affinity ARM-RNA recognition.     

Computational analysis of hydrogen bonds in protein-RNA complexes for interaction patterns

Structural analysis of protein-RNA complexes is labor-intensive, yet provides insight into the interaction patterns between a protein and RNA. As the number of protein-RNA complex structures reported has increased substantially in the last few years, a systematic method is required for automatically identifying interaction patterns. This paper presents a computational analysis of the hydrogen bonds in the most representative set of protein-RNA complexes. The analysis revealed several interesting interaction patterns. (1) While residues in the beta-sheets favored unpaired nucleotides, residues in the helices showed no preference and residues in turns favored paired nucleotides. (2) The backbone hydrogen bonds were more dominant than the base hydrogen bonds in the paired nucleotides, but the reverse was observed in the unpaired nucleotides. (3) The protein-RNA complexes contained more paired nucleotides than unpaired nucleotides, but the unpaired nucleotides were observed more frequently interacting with the proteins. And (4) Arg-U, Thr-A, Lys-A, and Asn-U were the most frequently observed pairs. The interaction patterns discovered from the analysis will provide us with useful information in predicting the structure of the RNA binding protein and the structure of the protein binding RNA.     

Use of biotin-labeled nucleic acids for protein purification and agarose-based chemiluminescent electromobility shift assays

We have employed biotin-labeled RNA to serve two functions. In one, the biotin tethers the RNA to streptavidin-agarose beads, creating an affinity resin for protein purification. In the other, the biotin functions as a label for use in a modified chemiluminescent electromobility shift assay (EMSA), a technique used to detect the formation of protein-RNA complexes. The EMSA that we describe avoids the use not only of radioactivity but also of neurotoxic acrylamide by using agarose as the gel matrix in which the free nucleic acid is separated from protein-nucleic acid complexes. After separation of free from complexed RNA in agarose, the RNA is electroblotted to positively charged nylon. The biotin-labeled RNA is readily bound by a streptavidin-alkaline phosphatase conjugate, allowing for very sensitive chemiluminescent detection ( approximately 0.1-1.0 fmol limit). Using our system, we were able to purify both known iron-responsive proteins (IRPs) from rat liver and assess their binding affinity to RNA containing the iron-responsive element (IRE) using the same batch of biotinylated RNA. We show data indicating that agarose is especially useful for cases when large complexes are formed, although smaller complexes are even better resolved.     

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.     

A protein-RNA docking benchmark (I): nonredundant cases

We have developed a nonredundant protein-RNA docking benchmark dataset, which is derived from the available bound and unbound structures in the Protein Data Bank involving polypeptide and nucleic acid chains. It consists of nine unbound-unbound cases where both the protein and the RNA are available in the free form. The other 36 cases are of unbound-bound type where only the protein is available in the free form. The conformational change upon complex formation is calculated by a distance matrix alignment method, and based on that, complexes are classified into rigid, semi-flexible, and full flexible. Although in the rigid body category, no significant conformational change accompanies complex formation, the fully flexible test cases show large domain movements, RNA base flips, etc. The benchmark covers four major groups of RNA, namely, t-RNA, ribosomal RNA, duplex RNA, and single-stranded RNA. We find that RNA is generally more flexible than the protein in the complexes, and the interface region is as flexible as the molecule as a whole. The structural diversity of the complexes in the benchmark set should provide a common ground for the development and comparison of the protein-RNA docking methods. The benchmark can be freely downloaded from the internet.     

Replication of poliovirus RNA with complete internal ribosome entry site deletions

cis-acting RNA sequences and structures in the 5' and 3' nontranslated regions of poliovirus RNA interact with host translation machinery and viral replication proteins to coordinately regulate the sequential translation and replication of poliovirus RNA. The poliovirus internal ribosome entry site (IRES) in the 5' nontranslated region (NTR) has been implicated as a cis-active RNA required for both viral mRNA translation and viral RNA replication. To evaluate the role of the IRES in poliovirus RNA replication, we exploited the advantages of cell-free translation-replication reactions and preinitiation RNA replication complexes. Genetic complementation with helper mRNAs allowed us to create preinitiation RNA replication complexes containing RNA templates with defined deletions in the viral open reading frame and the IRES. A series of deletions revealed that no RNA elements of either the viral open reading frame or the IRES were required in cis for negative-strand RNA synthesis. The IRES was dispensable for both negative- and positive-strand RNA syntheses. Intriguingly, although small viral RNAs lacking the IRES replicated efficiently, the replication of genome length viral RNAs was stimulated by the presence of the IRES. These results suggest that RNA replication is not directly dependent on a template RNA first functioning as an mRNA. These results further suggest that poliovirus RNA replication is not absolutely dependent on any protein-RNA interactions involving the IRES.     

Template recognition mechanisms by replicase proteins differ between bipartite positive-strand genomic RNAs of a plant virus

Recognition of RNA templates by viral replicase proteins is one of the key steps in the replication process of all RNA viruses. However, the mechanisms underlying this phenomenon, including primary RNA elements that are recognized by the viral replicase proteins, are not well understood. Here, we used aptamer pulldown assays with membrane fractionation and protein-RNA coimmunoprecipitation in a cell-free viral translation/replication system to investigate how viral replicase proteins recognize the bipartite genomic RNAs of the Red clover necrotic mosaic virus (RCNMV). RCNMV replicase proteins bound specifically to a Y-shaped RNA element (YRE) located in the 3' untranslated region (UTR) of RNA2, which also interacted with the 480-kDa replicase complexes that contain viral and host proteins. The replicase-YRE interaction recruited RNA2 to the membrane fraction. Conversely, RNA1 fragments failed to interact with the replicase proteins supplied in trans. The results of protein-RNA coimmunoprecipitation assays suggest that RNA1 interacts with the replicase proteins coupled with their translation. Thus, the initial template recognition mechanisms employed by the replicase differ between RCNMV bipartite genomic RNAs and RNA elements are primary determinants of the differential replication mechanism.     

A basic region neighboring the lysine-rich C-terminus of protein SRP19 is required for binding to signal recognition particle RNA.

To identify some of the determinants in the 19-kilodalton protein of signal recognition particle (SRP19) for binding to signal recognition particle RNA, two mutant derivatives of the SRP19 were constructed, lacking 14 and 24 C-terminal amino acids. Polypeptides were transcribed and translated in vitro and tested for their ability to bind to signal recognition particle RNA by retention of protein-RNA complexes on DEAE-Sepharose. Both mutant polypeptides form complexes with the RNA, demonstrating that the 24 C-terminal amino acids, which include a lysine-rich sequence at positions 136-144, are dispensable. A third mutant polypeptide, in which eight additional amino acids were removed by oligonucleotide-directed digestion of the mRNA, was unable to bind. The amino acids in the sequence PKLKTRTQ correspond to positions 113-120; they are suggested to be involved in interaction with signal recognition particle RNA.     

Template activity of complexes formed between bacteriophage f2 RNA and coat protein.

Formation of complexes between f2 RNA polymerase cistron was partially inhibited, some RNA and coat protein was studied using salt conditions which are optimum for phage protein synthesis. In this ionic environment, coat protein precipitation can be prevented by sulfhydryl group-protecting agents. Complexes formed at different protein-RNA input molar ratios were isolated and tested for template activity in an in vitro protein synthesizing system. Simultaneously, the number of protein molecules bound per RNA strand in such complexes was measured by the membrane (Millipore) filtration technique. Under conditions in which translation of the RNA strands were complexed with six molecules of coat protein, whereas some remained unbound. Strong inhibition of the translation of the RNA polymerase cistron was observed when each of the RNA strands present in the mixture was associated with six molecules of coat protein.     

Structure-based analysis of protein-RNA interactions using the program ENTANGLE

Until recently, drawing general conclusions about RNA recognition by proteins has been hindered by the paucity of high-resolution structures. We have analyzed 45 PDB entries of protein-RNA complexes to explore the underlying chemical principles governing both specific and non-sequence specific binding. To facilitate the analysis, we have constructed a database of interactions using ENTANGLE, a JAVA-based program that uses available structural models in their PDB format and searches for appropriate hydrogen bonding, stacking, electrostatic, hydrophobic and van der Waals interactions. The resulting database of interactions reveals correlations that suggest the basis for the discrimination of RNA from DNA and for base-specific recognition. The data illustrate both major and minor interaction strategies employed by families of proteins such as tRNA synthetases, ribosomal proteins, or RNA recognition motifs with their RNA targets. Perhaps most surprisingly, specific RNA recognition appears to be mediated largely by interactions of amide and carbonyl groups in the protein backbone with the edge of the RNA base. In cases where a base accepts a proton, the dominant amino acid donor is arginine, whereas in cases where the base donates a proton, the predominant acceptor is the backbone carbonyl group, not a side-chain group. This is in marked contrast to DNA-protein interactions, which are governed predominantly by amino acid side-chain interactions with functional groups that are presented in the accessible major groove. RNA recognition often proceeds through loops, bulges, kinks and other irregular structures that permit use of all the RNA functional groups and this is seen throughout the protein-RNA interaction database.