A modified resonant recognition model to predict protein-protein interaction

Proteins are fundamental components of all living cells and the protein-protein interaction plays an important role in vital movement. This paper briefly introduced the original Resonant Recognition Model (RRM), and then modified it by using the wavelet transform to acquire the Modified Resonant Recognition Model (MRRM). The key characteristic of the new model is that it can predict directly the protein-protein interaction from the primary sequence, and the MRRM is more suitable than the RRM for this prediction. The results of numerical experiments show that the MRRM is effective for predicting the protein-protein interaction. Translated from Journal of Shanghai University (Natural Science), 2006, 12(1): 69–73 [译自: 上海大学学报(自然科学版)]     

Prediction of "hot spots" in interleukin-2 based on informational spectrum characteristics of growth-regulating factors. Comparison with experimental data

The Resonant Recognition Model (RRM) is a theoretical method for analysis of protein and nucleotide sequences, based on the Fourier transform of the numerical representation of sequences. The amplitude spectrum of this transform is designated Informational Spectrum (IS). There are certain common frequencies in IS of growth-regulating factors. These characteristic frequencies may correlate with their roles in cell proliferation and metabolism, and in antitumor activity. IS of IL-2 has prominent characteristics in the main frequency domain of growth factors, frequency domain of antitumor factors, and frequency domain characteristic for IL-2-alpha receptor. By means of the inverse method for these 3 domains, the amino acids in the sequence of human IL-2 that may be relevant to its biological function, the so-called "hot spots", were predicted. The most probable hot spots, obtained in this way, are in the potential binding site of IL-2 to its receptor, which agrees with experimental data.     

Bioactive peptide design using the Resonant Recognition Model

With a large number of DNA and protein sequences already known, the crucial question is to find out how the biological function of these macromolecules is "written" in the sequence of nucleotides or amino acids. Biological processes in any living organism are based on selective interactions between particular bio-molecules, mostly proteins. The rules governing the coding of a protein's biological function, i.e. its ability to selectively interact with other molecules, are still not elucidated. In addition, with the rapid accumulation of databases of protein primary structures, there is an urgent need for theoretical approaches that are capable of analysing protein structure-function relationships. The Resonant Recognition Model (RRM) [1, 2] is one attempt to identify the selectivity of protein interactions within the amino acid sequence. The RRM [1, 2] is a physico-mathematical approach that interprets protein sequence linear information using digital signal processing methods. In the RRM the protein primary structure is represented as a numerical series by assigning to each amino acid in the sequence a physical parameter value relevant to the protein's biological activity. The RRM concept is based on the finding that there is a significant correlation between spectra of the numerical presentation of amino acids and their biological activity. Once the characteristic frequency for a particular protein function/interaction is identified, it is possible then to utilize the RRM approach to predict the amino acids in the protein sequence, which predominantly contribute to this frequency and thus, to the observed function, as well as to design de novo peptides having the desired periodicities. As was shown in our previous studies of fibroblast growth factor (FGF) peptidic antagonists [2, 3] and human immunodeficiency virus (HIV) envelope agonists [2, 4], such de novo designed peptides express desired biological function. This study utilises the RRM computational approach to the analysis of oncogene and proto-oncogene proteins. The results obtained have shown that the RRM is capable of identifying the differences between the oncogenic and proto-oncogenic proteins with the possibility of identifying the "cancer-causing" features within their protein primary structure. In addition, the rational design of bioactive peptide analogues displaying oncogenic or proto-oncogenic-like activity is presented here.     

Biological effects of a de novo designed myxoma virus peptide analogue: evaluation of cytotoxicity on tumor cells

Background

The Resonant Recognition Model (RRM) is a physico-mathematical model that interprets protein sequence linear information using digital signal processing methods. In this study the RRM concept was employed for structure-function analysis of myxoma virus (MV) proteins and the design of a short bioactive therapeutic peptide with MV-like antitumor/cytotoxic activity.

Methodology/Principal Findings

The analogue RRM-MV was designed by RRM as a linear 18 aa 2.3 kDa peptide. The biological activity of this computationally designed peptide analogue against cancer and normal cell lines was investigated. The cellular cytotoxicity effects were confirmed by confocal immunofluorescence microscopy, by measuring the levels of cytoplasmic lactate dehydrogenase (LDH) and by Prestoblue cell viability assay for up to 72 hours in peptide treated and non-treated cell cultures. Our results revealed that RRM-MV induced a significant dose and time-dependent cytotoxic effect on murine and human cancer cell lines. Yet, when normal murine cell lines were similarly treated with RRM-MV, no cytotoxic effects were observed. Furthermore, the non-bioactive RRM designed peptide RRM-C produced negligible cytotoxic effects on these cancer and normal cell lines when used at similar concentrations. The presence/absence of phosphorylated Akt activity in B16F0 mouse melanoma cells was assessed to indicate the possible apoptosis signalling pathway that could be affected by the peptide treatment. So far, Akt activity did not seem to be significantly affected by RRM-MV as is the case for the original viral protein.

Conclusions/Significance

Our findings indicate the successful application of the RRM concept to design a bioactive peptide analogue (RRM-MV) with cytotoxic effects on tumor cells only. This 2.345 kDa peptide analogue to a 49 kDa viral protein may be suitable to be developed as a potential cancer therapeutic. These results also open a new direction to the rational design of therapeutic agents for future cancer treatment.     

Identification of domains in apobec-1 complementation factor required for RNA binding and apolipoprotein-B mRNA editing

The C-to-U editing of apolipoprotein-B (apo-B) mRNA is catalyzed by an enzyme complex that recognizes an 11-nt mooring sequence downstream of the editing site. A minimal holoenzyme that edits apo-B mRNA in vitro has been defined. This complex contains apobec-1, the catalytic subunit, and apobec-1 complementation factor (ACF), the RNA-binding subunit that binds to the mooring sequence. Here, we show that ACF binds with high affinity to single-stranded but not double-stranded apo-B mRNA. ACF contains three nonidentical RNA recognition motifs (RRM) and a unique C-terminal auxiliary domain. In many multi-RRM proteins, the RRMs mediate RNA binding and an auxiliary domain functions in protein-protein interactions. Here we show that ACF does not fit this simple model. Based on deletion mutagenesis, the RRMs in ACF are necessary but not sufficient for binding to apo-B mRNA. Amino acids in the pre-RRM region are required for complementing activity and RNA binding, but not for interaction with apobec-1. The C-terminal 196 amino acids are not absolutely essential for function. However, further deletion of an RG-rich region from the auxiliary domain abolished complementing activity, RNA binding, and apobec-1 interaction. The auxiliary domain alone did not bind apobec-1. Although all three RRMs are required for complementing activity and apobec-1 interaction, the individual motifs contribute differently to RNA binding. Point mutations in RRM1 or RRM2 decreased the Kd for apo-B mRNA by two orders of magnitude whereas mutations in RRM3 reduced binding affinity 13-fold. The pairwise expression of RRM1 with RRM2 or RRM3 resulted in moderate affinity binding.     

Dynamic property is a key determinant for protein-protein interactions

Dynamic property is highly correlated with the biological functions of macromolecules, such as the activity and specificity of enzymes and the allosteric regulation in the signal transduction process. Applications of the dynamic property to protein function researches have been discussed and encouraging progresses have been achieved, for example, in enzyme activity and protein-protein docking studies. However, how the global dynamic property contributes to protein-protein interaction was still unclear. We have studied the dynamic property in protein-protein interactions based on Gaussian Network Model and applied it to classify biological and nonbiological protein-protein complexes in crystal structures. The global motion correlation between residues from the two protomers was found to be remarkably different for biological and nonbiological complexes. This correlation has been used to discriminate biological and nonbiological complexes in crystal and gave a classification rate of 86.9% in the cross-validation test. The innovation of this feature is that it is a global dynamic property which does not rely directly on the interfacial properties of the complex. In addition, the correlation of the global motions was found to be weakly correlated with the dissociation rate constant of protein complexes. We suggest that the dynamic property is a key determinant for protein-protein interaction, which can be used to discriminate native and crystal complexes and potentially be applied in protein-protein dynamic rate constants estimations.     

Solution structure of the first RNA recognition motif domain of human spliceosomal protein SF3b49 and its mode of interaction with a SF3b145 fragment

The spliceosomal protein SF3b49, a component of the splicing factor 3b (SF3b) protein complex in the U2 small nuclear ribonucleoprotein, contains two RNA recognition motif (RRM) domains. In yeast, the first RRM domain (RRM1) of Hsh49 protein (yeast orthologue of human SF3b49) reportedly interacts with another component, Cus1 protein (orthologue of human SF3b145). Here, we solved the solution structure of the RRM1 of human SF3b49 and examined its mode of interaction with a fragment of human SF3b145 using NMR methods. Chemical shift mapping showed that the SF3b145 fragment spanning residues 598–631 interacts with SF3b49 RRM1, which adopts a canonical RRM fold with a topology of β1‐α1‐β2‐β3‐α2‐β4. Furthermore, a docking model based on NOESY measurements suggests that residues 607–616 of the SF3b145 fragment adopt a helical structure that binds to RRM1 predominantly via α1, consequently exhibiting a helix–helix interaction in almost antiparallel. This mode of interaction was confirmed by a mutational analysis using GST pull‐down assays. Comparison with structures of all RRM domains when complexed with a peptide found that this helix–helix interaction is unique to SF3b49 RRM1. Additionally, all amino acid residues involved in the interaction are well conserved among eukaryotes, suggesting evolutionary conservation of this interaction mode between SF3b49 RRM1 and SF3b145.     

Structural and functional characterization of the TgDRE multidomain protein, a DNA repair enzyme from Toxoplasma gondii

The parasite Toxoplasma gondii expresses a 55 kDa protein or TgDRE that belongs to a novel family of proteins characterized by the presence of three domains, a human splicing factor 45-like motif (SF), a glycine-rich motif (G-patch), and a RNA recognition motif (RRM). The two latter domains are mainly known as RNA-binding domains, and their presence in TgDRE, whose partial DNA repair function was demonstrated, suggests that the protein could also be involved in the RNA metabolism. In this work, we characterized the structure and function of the different domains by using single or multidomain proteins to define their putative role. The SF45-like domain has a helical conformation and is involved in the oligomerization of the protein. The G-patch domain, mainly unstructured on its own as well as in the presence of the SF upstream and RRM downstream domains, is able to bind small RNA oligonucleotides. We also report the structure determination of the RRM domain from the NMR data. It adopts a classical betaalphabetabetaalphabeta topology consisting of a four-stranded beta sheet packed against two alpha helices but does not present the key residues for the RNA interaction. In contrast, our analysis shows that the RRM of TgDRE is not only unable to bind small RNA oligonucleotides but it also shares the protein-protein interaction characteristics with two unusual RRMs of the U2AF heterodimeric splicing factor. The presence of both RNA- and protein-binding domains seems to indicate that TgDRE could also be involved in RNA metabolism.     

Escher: A new docking procedure applied to the reconstruction of protein tertiary structure

Evaluation of Surface Complementarity, Hydrogen bonding, and Electrostatic interaction in molecular Recognition (ESCHER) is a new docking procedure consisting of three modules that work in series. The first module evaluates the geometric complementarity and produces a set of rough solutions for the docking problem. The second module identifies molecular collisions within those solutions, and the third evaluates their electrostatic complementarity. We describe the algorithm and its application to the docking of cocrystallized protein domains and unbound components of protein-protein complexes. Furthermore, ESCHER has been applied to the reassociation of secondary and supersecondary structure elements. The possibility of applying a docking method to the problem of protein structure prediction is discussed. Proteins 28:556–567, 1997. © 1997 Wiley-Liss, Inc.     

The conserved RNA recognition motif 3 of U2 snRNA auxiliary factor (U2AF 65) is essential in vivo but dispensable for activity in vitro

The general splicing factor U2AF(65) recognizes the polypyrimidine tract (Py tract) that precedes 3' splice sites and has three RNA recognition motifs (RRMs). The C-terminal RRM (RRM3), which is highly conserved, has been proposed to contribute to Py-tract binding and establish protein-protein contacts with splicing factors mBBP/SF1 and SAP155. Unexpectedly, we find that the human RRM3 domain is dispensable for U2AF(65) activity in vitro. However, it has an essential function in Schizosaccharomyces pombe distinct from binding to the Py tract or to mBBP/SF1 and SAP155. First, deletion of RRM3 from the human protein has no effect on Py-tract binding. Second, RRM123 and RRM12 select similar sequences from a random pool of RNA. Third, deletion of RRM3 has no effect on the splicing activity of U2AF(65) in vitro. However, deletion of the RRM3 domain of S. pombe U2AF(59) abolishes U2AF function in vivo. In addition, certain amino acid substitutions on the four-stranded beta-sheet surface of RRM3 compromise U2AF function in vivo without affecting binding to mBBP/SF1 or SAP155 in vitro. We propose that RRM3 has an unrecognized function that is possibly relevant for the splicing of only a subset of cellular introns. We discuss the implications of these observations on previous models of U2AF function.     

Structure/function relationships within the RNA recognition motif family applied to the hermes gene product

The RNA Recognition Motif (RRM) family of RNA-binding domains comprises distinct structural subclasses which can be equated to various types of cognate RNA(s) in relation to biological functions. By identifying structural templates within the appropriate RRM subclass, we have homology-modelled the three-dimensional structure of the hermes gene-encoded RRM. Our findings lead us to propose potential RNA targets for the corresponding protein and to predict possible functions in RNA metabolism during heart development.     

An unusual RNA recognition motif acts as a scaffold for multiple proteins in the pre-mRNA retention and splicing complex

The yeast pre-mRNA retention and splicing complex counteracts the escape of unspliced pre-mRNAs from the nucleus and activates splicing of a subset of Mer1p-dependent genes. A homologous complex is present in activated human spliceosomes. In many components of the spliceosome, RNA recognition motifs (RRMs) serve as versatile protein-RNA or protein-protein interaction platforms. Here, we show that in the retention and splicing complex, an atypical RRM of the Snu17p (small nuclear ribonucleoprotein-associated protein 17) subunit acts as a scaffold that organizes the other two constituents, Bud13p (bud site selection 13) and Pml1p (pre-mRNA leakage 1). GST pull-down experiments and size exclusion chromatography revealed that Snu17p constitutes the central platform of the complex, whereas Bud13p and Pml1p do not interact with each other. Fluorimetric structure probing showed the entire Bud13p and the N-terminal third of Pml1p to be natively disordered in isolation. Mutational analysis and tryptophan fluorescence confirmed that a conserved tryptophan-containing motif in the C terminus of Bud13p binds to the core RRM of Snu17p, whereas a different interaction surface encompassing a C-terminal extension of the Snu17p RRM is required to bind an N-terminal peptide of Pml1p. Isothermal titration calorimetry revealed 1:1 interaction stoichiometries, large negative binding entropies, and dissociation constants in the low nanomolar and micromolar ranges for the Snu17p-Bud13p and the Snu17p-Pml1p interactions, respectively. Our results demonstrate that the noncanonical Snu17p RRM concomitantly binds multiple ligand proteins via short, intrinsically unstructured peptide epitopes and thereby acts as a platform that displays functional modules of the ligands, such as a forkhead-associated domain of Pml1p and a conserved polylysine motif of Bud13p.     

A novel peptide recognition mode revealed by the X-ray structure of a core U2AF35/U2AF65 heterodimer

U2 auxiliary factor (U2AF) is an essential splicing factor that recognizes the 3' splice site and recruits the U2 snRNP to the branch point. The X-ray structure of the human core U2AF heterodimer, consisting of the U2AF35 central domain and a proline-rich region of U2AF65, has been determined at 2.2 A resolution. The structure reveals a novel protein-protein recognition strategy, in which an atypical RNA recognition motif (RRM) of U2AF35 and the U2AF65 polyproline segment interact via reciprocal "tongue-in-groove" tryptophan residues. Complementary biochemical experiments demonstrate that the core U2AF heterodimer binds RNA, and that the interacting tryptophan side chains are essential for U2AF dimerization. Atypical RRMs in other splicing factors may serve as protein-protein interaction motifs elsewhere during spliceosome assembly.     

Structural characterization of Set1 RNA recognition motifs and their role in histone H3 lysine 4 methylation

The yeast Set1 histone H3 lysine 4 (H3K4) methyltransferase contains, in addition to its catalytic SET domain, a conserved RNA recognition motif (RRM1). We present here the crystal structure and the secondary structure assignment in solution of the Set1 RRM1. Although RRM1 has the expected betaalphabetabetaalphabeta RRM-fold, it lacks the typical RNA-binding features of these modules. RRM1 is not able to bind RNA by itself in vitro, but a construct combining RRM1 with a newly identified downstream RRM2 specifically binds RNA. In vivo, H3K4 methylation is not affected by a point mutation in RRM2 that preserves Set1 stability but affects RNA binding in vitro. In contrast mutating RRM1 destabilizes Set1 and leads to an increase of dimethylation of H3K4 at the 5'-coding region of active genes at the expense of trimethylation, whereas both, dimethylation decreases at the 3'-coding region. Taken together, our results suggest that Set1 RRMs bind RNA, but Set1 RNA-binding activity is not linked to H3K4 methylation.     

Adaptable molecular interactions guide phosphorylation of the SR protein ASF/SF2 by SRPK1

The SR (arginine-serine rich) protein ASF/SF2 (also called human alternative splicing factor), an essential splicing factor, contains two functional modules consisting of tandem RNA recognition motifs (RRMs; RRM1-RRM2) and a C-terminal arginine-serine repeat region (RS domain, a domain rich in arginine-serine repeats). The SR-specific protein kinase (SRPK) 1 phosphorylates the RS domain at multiple serines using a directional (C-terminal-to-N-terminal) and processive mechanism—a process that directs the SR protein to the nucleus and influences protein-protein interactions associated with splicing function. To investigate how SRPK1 accomplishes this feat, the enzyme-substrate complex was analyzed using single-turnover and multiturnover kinetic methods. Deletion studies revealed that while recognition of the RS domain by a docking groove on SRPK1 is sufficient to initiate the processive and directional mechanism, continued processive phosphorylation in the presence of building repulsive charge relies on the fine-tuning of contacts with the RRM1-RRM2 module. An electropositive pocket in SRPK1 that stabilizes newly phosphorylated serines enhanced processive phosphorylation of later serines. These data indicate that SRPK1 uses stable, yet highly flexible protein-protein interactions to facilitate both early and late phases of the processive phosphorylation of SR proteins.     

The role of the C-terminal helix of U1A protein in the interaction with U1hpII RNA

Previous kinetic investigations of the N-terminal RNA Recognition Motif (RRM) domain of spliceosomal A protein of the U1 small nuclear ribonucleoprotein particle (U1A) interacting with its RNA target U1 hairpin II (U1hpII) provided experimental evidence for a ‘lure and lock’ model of binding. The final step of locking has been proposed to involve conformational changes in an α-helix immediately C-terminal to the RRM domain (helix C), which occludes the RNA binding surface in the unbound protein. Helix C must shift its position to accommodate RNA binding in the RNA–protein complex. This results in a new hydrophobic core, an intraprotein hydrogen bond and a quadruple stacking interaction between U1A and U1hpII. Here, we used a surface plasmon resonance-based biosensor to gain mechanistic insight into the role of helix C in mediating the interaction with U1hpII. Truncation, removal or disruption of the helix exposes the RNA-binding surface, resulting in an increase in the association rate, while simultaneously reducing the ability of the complex to lock, reflected in a loss of complex stability. Disruption of the quadruple stacking interaction has minor kinetic effects when compared with removal of the intraprotein hydrogen bonds. These data provide new insights into the mechanism whereby sequences C-terminal to an RRM can influence RNA binding.     

eIF4G, eIFiso4G, and eIF4B bind the poly(A)-binding protein through overlapping sites within the RNA recognition motif domains

The poly(A)-binding protein (PABP), a protein that contains four conserved RNA recognition motifs (RRM1-4) and a C-terminal domain, is expressed throughout the eukaryotic kingdom and promotes translation through physical and functional interactions with eukaryotic initiation factor (eIF) 4G and eIF4B. Two highly divergent isoforms of eIF4G, known as eIF4G and eIFiso4G, are expressed in plants. As little is known about how PABP can interact with RNA and three distinct translation initiation factors in plants, the RNA binding specificity and organization of the protein interaction domains in wheat PABP was investigated. Wheat PABP differs from animal PABP in that its RRM1 does not bind RNA as an individual domain and that RRM 2, 3, and 4 exhibit different RNA binding specificities to non-poly(A) sequences. The PABP interaction domains for eIF4G and eIFiso4G were distinct despite the functional similarity between the eIF4G proteins. A single interaction domain for eIF4G is present in the RRM1 of PABP, whereas eIFiso4G interacts at two sites, i.e. one within RRM1-2 and the second within RRM3-4. The eIFiso4G binding site in RRM1-2 mapped to a 36-amino acid region encompassing the C-terminal end of RRM1, the linker region, and the N-terminal end of RRM2, whereas the second site in RRM3-4 was more complex. A single interaction domain for eIF4B is present within a 32-amino acid region representing the C-terminal end of RRM1 of PABP that overlaps with the N-proximal eIFiso4G interaction domain. eIF4B and eIFiso4G exhibited competitive binding to PABP, supporting the overlapping nature of their interaction domains. These results support the notion that eIF4G, eIFiso4G, and eIF4B interact with distinct molecules of PABP to increase the stability of the interaction between the termini of an mRNA.