Single-stranded RNA Binding Protein Research Articles

Transcriptome-wide identification of single-stranded RNA binding proteins.

RNA-protein interactions are precisely regulated by RNA secondary structures in various biological processes. Large-scale identification of proteins that interact with particular RNA structure is important to the RBPome. Herein, a kethoxal assisted single-stranded RNA interactome capture (KASRIC) strategy was developed to globally identify single-stranded RNA binding proteins (ssRBPs). This approach combines RNA secondary structure probing technology with the conventional method of RNA-binding proteins profiling, realizing the transcriptome-wide identification of ssRBPs. Applying KASRIC, we identified 3180 candidate RBPs and 244 candidate ssRBPs in HeLa cells. Importantly, the 244 candidate ssRBPs contained 55 previously reported ssRBPs and 189 novel ssRBPs. Function analysis of the candidate ssRBPs exhibited enrichment in cellular processes related to RNA splicing and RNA degradation. The KASRIC strategy will facilitate the investigation of RNA-protein interactions.

Identification of Igf2bp1 gene family and effect on chicken myoblast proliferation

ABSTRACT The function of Igf2bp1 regulated muscle development in Lueyang black-bone chickens is to be investigated. The insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) belongs to the family of single-stranded RNA-binding proteins (IGF2BP1-3). IGF2BP1 is abundantly expressed in a broad range of foetal tissues, but only in a limited number of normal adult tissues, and is involved in biological regulation processes, including cell proliferation, differentiation, metabolism and migration, especially in tumour cells. Recent study has reported Igf2bp1 was a key gene that regulated the body size of Pekin ducks. We found that the expression level of Igf2bp1 gradually decreased with age in the leg muscle of Lueyang black-bone chicken, but its function in chicken muscle development is unknown. Firstly, the evolutionary relationship of Igf2bp1 gene family was analysed in chicken. Additionally, Igf2bp1-overexpressed myoblasts of Lueyang black-bone chicken were constructed. The results showed that Igf2bp1 plays a negative role in cell proliferation by the CCK8 and 5-ethynyl-2′-deoxyuridine assay. Moreover, Igf2bp1 inhibited myogenesis by detecting the myoblast proliferation markers of pax7, the expression of myogenic differentiation markers, MyoD, MyoG and Myf5 in myoblasts. Thus, Igf2bp1 can inhibit chicken myoblast proliferation and differentiation by regulating key gene expression of myogenesis, which is not conducive to the growth and development of skeletal muscle in Lueyang black chickens.

Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins.

Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA–ssDBP and ssRNA–ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base–aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA–protein and ssRNA–protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA–ssRBP complexes compared with the ssDNA–ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (koff) for complexes involving the former.

NS38 is required for aquareovirus replication via interaction with viral core proteins and host eIF3A

Aquareoviruses contain an 11-segmented double-stranded RNA genome. Previous studies indicated that NS38, a virus-encoded putative single-stranded RNA binding protein, interacts with NS80 in viral inclusion bodies (VIBs). However, the role of NS38 in aquareovirus infection remained unclear. Here, we found that NS38 interacts with inner-capsid proteins (VP1–VP4 and VP6) and the NS80-RNA complex in both transfected and infected cells. Knockdown of NS38 by siRNAs-115/219 clearly reduced viral infection, with decreased mRNA and protein yields. Moreover, NS38 can interact with host cellular eukaryotic translation initiation factor 3 subunit A (eIF3A) in transfected cells, while no association was detected between eIF3A and NS80. This study is the first to define that the NS38 is essential to viral replication. Together, our findings indicate that NS38 might function as a mediator by interacting with viral and host cellular components in VIBs during replication.

Abstract 5181: KSRP suppresses cell invasion and metastasis through miR-23a-mediated EGR3 mRNA degradation in non-small cell lung cancer

Abstract KH-type splicing regulatory protein (KSRP) is a single-strand RNA binding protein which regulates mRNA stability either by binding to AU-rich elements (AREs) of mRNA 3′UTR or by facilitating miRNA biogenesis to target mRNA. Unlike its well-characterized function at the molecular level in maintaining RNA homeostasis, the role of KSRP in cancer progression remains largely unknown. Here we investigate the role of KSRP in non-small cell lung cancer (NSCLC). We first examined KSRP expression by immunohistochemistry in a cohort containing 196 NSCLC patients and observed a strong positive correlation between KSRP expression and survival of NSCLC patients. Multivariate analysis further identified KSRP as an independent prognostic factor. Manipulating KSRP expression significantly affected in vitro cell mobility and in vivo metastatic ability of NSCLC cells. Microarray analysis identified an ARE-containing gene, EGR3, as a downstream effector of KSRP in NSCLC. Interestingly, we found that KSRP decreased EGR3 mRNA stability in an ARE-independent manner. By screening KSRP-regulated miRNAs in NSCLC cells, we further found that miR-23a directly binds to EGR3 3′UTR, reducing EGR3 expression and thereby inhibiting NSCLC cell mobility. Our findings implicate a targetable KSRP/miR-23a/EGR3 signaling axis in advanced tumor phenotypes Citation Format: Ke-Fan Pan, Ming-Hsien Chien, Wei-Jiunn Lee, Michael Hsiao, Kuo-Tai Hua. KSRP suppresses cell invasion and metastasis through miR-23a-mediated EGR3 mRNA degradation in non-small cell lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5181.

β-arrestin-biased agonism of β-adrenergic receptor regulates Dicer-mediated microRNA maturation to promote cardioprotective signaling

RationaleMicroRNAs (miRs) are small, non-coding RNAs that function to post-transcriptionally regulate target genes. First transcribed as primary miR transcripts (pri-miRs), they are enzymatically processed by Drosha into premature miRs (pre-miRs) and further cleaved by Dicer into mature miRs. Initially discovered to desensitize β-adrenergic receptor (βAR) signaling, β-arrestins are now well-appreciated to modulate multiple pathways independent of G protein signaling, a concept known as biased signaling. Using the β-arrestin-biased βAR ligand carvedilol, we previously showed that β-arrestin1 (not β-arrestin2)-biased β1AR (not β2AR) cardioprotective signaling stimulates Drosha-mediated processing of six miRs by forming a multi-protein nuclear complex, which includes β-arrestin1, the Drosha microprocessor complex and a single-stranded RNA binding protein hnRNPA1. ObjectiveHere, we investigate whether β-arrestin-mediated βAR signaling induced by carvedilol could regulate Dicer-mediated miR maturation in the cytoplasm and whether this novel mechanism promotes cardioprotective signaling. Methods and resultsIn mouse hearts, carvedilol indeed upregulates three mature miRs, but not their pre-miRs and pri-miRs, in a β-arrestin 1- or 2-dependent manner. Interestingly, carvedilol-mediated activation of miR-466g or miR-532-5p, and miR-674 is dependent on β2ARs and β1ARs, respectively. Mechanistically, β-arrestin 1 or 2 regulates maturation of three newly identified βAR/β-arrestin-responsive miRs (β-miRs) by associating with the Dicer maturation RNase III enzyme on three pre-miRs of β-miRs. Myocardial cell approaches uncover that despite their distinct roles in different cell types, β-miRs act as gatekeepers of cardiac cell functions by repressing deleterious targets. ConclusionsOur findings indicate a novel role for βAR-mediated β-arrestin signaling activated by carvedilol in Dicer-mediated miR maturation, which may be linked to its protective mechanisms.

Engineering RNA-Binding Proteins by Modular Assembly of RanBP2-Type Zinc Fingers.

Deciphering the function of the nonprotein-coding portion of genomes represents one of the major challenges that molecular biology is facing today. Numerous classes of RNAs have been discovered over the last past decade and appear to play important regulatory roles in gene expression and disease. The ability to study and manipulate these RNAs relies on the development of programmable RNA-binding molecules such as RNA-binding proteins. Most RNA-binding proteins have modular architectures and combine different RNA-binding domains that provide binding affinity toward a specific RNA sequence and/or structure. Herein, we describe a general strategy to design single-stranded RNA-binding proteins using RanBP2-type zinc-finger (ZF) domains that can recognize a given RNA sequence with high affinity and specificity.

KSRP suppresses cell invasion and metastasis through miR-23a-mediated EGR3 mRNA degradation in non-small cell lung cancer

KH-type splicing regulatory protein (KSRP) is a single-strand RNA binding protein which regulates mRNA stability either by binding to AU-rich elements (AREs) of mRNA 3′UTR or by facilitating miRNA biogenesis to target mRNA. Unlike its well-characterized function at the molecular level in maintaining RNA homeostasis, the role of KSRP in cancer progression remains largely unknown. Here we investigate the role of KSRP in non-small cell lung cancer (NSCLC). We first examined KSRP expression by immunohistochemistry in a cohort containing 196 NSCLC patients and observed a strong positive correlation between KSRP expression and survival of NSCLC patients. Multivariate analysis further identified KSRP as an independent prognostic factor. Manipulating KSRP expression significantly affected in vitro cell mobility and in vivo metastatic ability of NSCLC cells. Microarray analysis identified an ARE-containing gene, EGR3, as a downstream effector of KSRP in NSCLC. Interestingly, we found that KSRP decreased EGR3 mRNA stability in an ARE-independent manner. By screening KSRP-regulated miRNAs in NSCLC cells, we further found that miR-23a directly binds to EGR3 3′UTR, reducing EGR3 expression and thereby inhibiting NSCLC cell mobility. Our findings implicate a targetable KSRP/miR-23a/EGR3 signaling axis in advanced tumor phenotypes.

Pseudouridine and N6-methyladenosine modifications weaken PUF protein/RNA interactions.

RNA modifications are ubiquitous in biology, with over 100 distinct modifications. While the vast majority were identified and characterized on abundant noncoding RNA such as tRNA and rRNA, the advent of sensitive sequencing-based approaches has led to the discovery of extensive and regulated modification of eukaryotic messenger RNAs as well. The two most abundant mRNA modifications—pseudouridine (Ψ) and N6-methyladenosine (m6A)—affect diverse cellular processes including mRNA splicing, localization, translation, and decay and modulate RNA structure. Here, we test the hypothesis that RNA modifications directly affect interactions between RNA-binding proteins and target RNA. We show that Ψ and m6A weaken the binding of the human single-stranded RNA binding protein Pumilio 2 (hPUM2) to its consensus motif, with individual modifications having effects up to approximately threefold and multiple modifications giving larger effects. While there are likely to be some cases where RNA modifications essentially fully ablate protein binding, here we see modest responses that may be more common. Such modest effects could nevertheless profoundly alter the complex landscape of RNA:protein interactions, and the quantitative rather than qualitative nature of these effects underscores the need for quantitative, systems-level accounting of RNA:protein interactions to understand post-transcriptional regulation.

MiRNA-Mediated KHSRP Silencing Rewires Distinct Post-transcriptional Programs during TGF-β-Induced Epithelial-to-Mesenchymal Transition

Epithelial-to-mesenchymal transition (EMT) confers several traits to cancer cells that are required for malignant progression. Here, we report that miR-27b-3p-mediated silencing of the single-strand RNA binding protein KHSRP is required for transforming growth factor β (TGF-β)-induced EMT in mammary gland cells. Sustained KHSRP expression limits TGF-β-dependent induction of EMT factors and cell migration, whereas its knockdown inuntreated cellsmimics TGF-β-induced EMT. Genome-wide sequencing analyses revealed that KHSRP controls (1) levels of mature miR-192-5p, a microRNA that targets a group of EMT factors, and (2) alternative splicing of a cohort of pre-mRNAs related to cell adhesion and motility including Cd44 and Fgfr2. KHSRP belongs to a ribonucleoprotein complex that includes hnRNPA1, and the two proteins cooperate in promoting epithelial-type exon usage of select pre-mRNAs. Thus, TGF-β-induced KHSRP silencing is central in a pathway leading to gene-expression changes that contribute to the cellular changes linked toEMT.

Synthesis and inhibitory effects of novel pyrimido-pyrrolo-quinoxalinedione analogues targeting nucleoproteins of influenza A virus H1N1

The influenza nucleoprotein (NP) is a single-strand RNA-binding protein and the core of the influenza ribonucleoprotein (RNP) particle that serves many critical functions for influenza replication. NP has been considered as a promising anti-influenza target. A new class of anti-influenza compounds, nucleozin and analogues were reported recently in several laboratories to inhibit the synthesis of influenza macromolecules and prevent the cytoplasmic trafficking of the influenza RNP. In this study, pyrimido-pyrrolo-quinoxalinedione (PPQ) analogues as a new class of novel anti-influenza agents are reported. Compound PPQ-581 was identified as a potential anti-influenza lead with EC50 value of 1 μM for preventing virus-induced cytopathic effects. PPQ produces similar anti-influenza effects as nucleozin does in influenza-infected cells. Treatment with PPQ at the beginning of H1N1 infection inhibited viral protein synthesis, while treatment at later times blocked the RNP nuclear export and the appearance of cytoplasmic RNP aggregation. PPQ resistant H1N1 (WSN) viruses were isolated and found to have a NPS377G mutation. Recombinant WSN carrying the S377G NP is resistant to PPQ in anti-influenza and RNA polymerase assays. The WSN virus with the NPS377G mutation also is devoid of the PPQ-mediated RNP nuclear retention and cytoplasmic aggregation. The NPS377G expressing WSN virus is not resistant to the reported NP inhibitors nucleozin. Similarly, the nucleozin resistant WSN viruses are not resistant to PPQ, suggesting that PPQ targets a different site from the nucleozin-binding site. Our results also suggest that NP can be targeted through various binding sites to interrupt the crucial RNP trafficking, resulting in influenza replication inhibition.

The design and structural characterization of a synthetic pentatricopeptide repeat protein.

Proteins of the pentatricopeptide repeat (PPR) superfamily are characterized by tandem arrays of a degenerate 35-amino-acid α-hairpin motif. PPR proteins are typically single-stranded RNA-binding proteins with essential roles in organelle biogenesis, RNA editing and mRNA maturation. A modular, predictable code for sequence-specific binding of RNA by PPR proteins has recently been revealed, which opens the door to the de novo design of bespoke proteins with specific RNA targets, with widespread biotechnological potential. Here, the design and production of a synthetic PPR protein based on a consensus sequence and the determination of its crystal structure to 2.2 Å resolution are described. The crystal structure displays helical disorder, resulting in electron density representing an infinite superhelical PPR protein. A structural comparison with related tetratricopeptide repeat (TPR) proteins, and with native PPR proteins, reveals key roles for conserved residues in directing the structure and function of PPR proteins. The designed proteins have high solubility and thermal stability, and can form long tracts of PPR repeats. Thus, consensus-sequence synthetic PPR proteins could provide a suitable backbone for the design of bespoke RNA-binding proteins with the potential for high specificity.

RNA structure generates natural cooperativity between single-stranded RNA binding proteins targeting 5' and 3'UTRs.

In post-transcriptional regulation, an mRNA molecule is bound by many proteins and/or miRNAs to modulate its function. To enable combinatorial gene regulation, these binding partners of an RNA must communicate with each other, exhibiting cooperativity. Even in the absence of direct physical interactions between the binding partners, such cooperativity can be mediated through RNA secondary structures, since they affect the accessibility of the binding sites. Here we propose a quantitative measure of this structure-mediated cooperativity that can be numerically calculated for an arbitrary RNA sequence. Focusing on an RNA with two binding sites, we derive a characteristic difference of free energy differences, i.e. ΔΔG, as a measure of the effect of the occupancy of one binding site on the binding strength of another. We apply this measure to a large number of human and Caenorhabditis elegans mRNAs, and find that structure-mediated cooperativity is a generic feature. Interestingly, this cooperativity not only affects binding sites in close proximity along the sequence but also configurations in which one binding site is located in the 5′UTR and the other is located in the 3′UTR of the mRNA. Furthermore, we find that this end-to-end cooperativity is determined by the UTR sequences while the sequences of the coding regions are irrelevant.

Virus-mediated compartmentalization of the host translational machinery.

ABSTRACTViruses require the host translational apparatus to synthesize viral proteins. Host stress response mechanisms that suppress translation, therefore, represent a significant obstacle that viruses must overcome. Here, we report a strategy whereby the mammalian orthoreoviruses compartmentalize the translational machinery within virus-induced inclusions known as viral factories (VF). VF are the sites of reovirus replication and assembly but were thought not to contain ribosomes. It was assumed viral mRNAs exited the VF to undergo translation by the cellular machinery, and proteins reentered the factory to participate in assembly. Here, we used ribopuromycylation to visualize active translation in infected cells. These studies revealed that active translation occurs within VF and that ribosomal subunits and proteins required for translation initiation, elongation, termination, and recycling localize to the factory. Interestingly, we observed components of the 43S preinitiation complex (PIC) concentrating primarily at factory margins, suggesting a spatial and/or dynamic organization of translation within the VF. Similarly, the viral single-stranded RNA binding protein σNS localized to the factory margins and had a tubulovesicular staining pattern that extended a short distance from the margins of the factories and colocalized with endoplasmic reticulum (ER) markers. Consistent with these colocalization studies, σNS was found to associate with both eukaryotic translation initiation factor 3 subunit A (eIF3A) and the ribosomal subunit pS6R. Together, these findings indicate that σNS functions to recruit 43S PIC machinery to the primary site of viral translation within the viral factory. Pathogen-mediated compartmentalization of the translational apparatus provides a novel mechanism by which viruses might avoid host translational suppression.

Crystal structure of the single-stranded RNA binding protein HutP from Geobacillus thermodenitrificans

RNA binding proteins control gene expression by the attenuation/antitermination mechanism. HutP is an RNA binding antitermination protein. It regulates the expression of hut operon when it binds with RNA by modulating the secondary structure of single-stranded hut mRNA. HutP necessitates the presence of l-histidine and divalent metal ion to bind with RNA. Herein, we report the crystal structures of ternary complex (HutP–l-histidine–Mg2+) and EDTA (0.5M) treated ternary complex (HutP–l-histidine–Mg2+), solved at 1.9Å and 2.5Å resolutions, respectively, from Geobacillus thermodenitrificans. The addition of 0.5M EDTA does not affect the overall metal-ion mediated ternary complex structure and however, the metal ions at the non-specific binding sites are chelated, as evidenced from the results of structural features.

Identification of proteins that form specific complexes with the highly conserved protein Translin in Schizosaccharomyces pombe

Translin is a single-stranded DNA and RNA binding protein that has a high affinity for G-rich sequences. TRAX is a Translin paralog that associates with Translin. Both Translin and TRAX were highly conserved in eukaryotes. The nucleic acid binding form of Translin is a barrel-shaped homo-octamer. A Translin–TRAX hetero-octamer having a similar structure also binds nucleic acids. Previous reports suggested that Translin may be involved in chromosomal translocations, telomere metabolism and the control of mRNA transport and translation. More recent studies have indicated that Translin–TRAX hetero-octamers are involved in RNA silencing. To gain a further insight into the functions of Translin, we have undertaken to systematically search for proteins with which it forms specific complexes in living cells. Here we report the results of such a search conducted in the fission yeast Schizosaccharomyces pombe, a suitable model system. This search was carried out by affinity purification and immuno-precipitation techniques, combined with differential labeling of the intracellular proteins with the stable isotopes 15N and 14N. We identified for the first time two proteins containing an RNA Recognition Motif (RRM), which are specifically associated with the yeast Translin: (1) the pre-mRNA-splicing factor srp1 that belongs to the highly conserved SR family of proteins and (2) vip1, a protein conserved in fungi. Our data also support the presence of RNA in these intracellular complexes. Our experimental approach should be generally applicable to studies of weak intracellular protein–protein interactions and provides a clear distinction between false positive vs. truly interacting proteins.

EndoU is a novel regulator of AICD during peripheral B cell selection

Balanced transmembrane signals maintain a competent peripheral B cell pool limited in self-reactive B cells that may produce pathogenic autoantibodies. To identify molecules regulating peripheral B cell survival and tolerance to self-antigens (Ags), a gene modifier screen was performed with B cells from CD22-deficient C57BL/6 (CD22(-/-[B6])) mice that undergo activation-induced cell death (AICD) and fail to up-regulate c-Myc expression after B cell Ag receptor ligation. Likewise, lysozyme auto-Ag-specific B cells in Ig(Tg) hen egg lysozyme (HEL) transgenic mice inhabit the spleen but undergo AICD after auto-Ag encounter. This gene modifier screen identified EndoU, a single-stranded RNA-binding protein of ancient origin, as a major regulator of B cell survival in both models. EndoU gene disruption prevents AICD and normalizes c-Myc expression. These findings reveal that EndoU is a critical regulator of an unexpected and novel RNA-dependent pathway controlling peripheral B cell survival and Ag responsiveness that may contribute to peripheral B cell tolerance.

Bone Morphogenetic Protein/SMAD Signaling Orients Cell Fate Decision by Impairing KSRP-Dependent MicroRNA Maturation

MicroRNAs (miRNAs) are essential regulators of development, physiology, and evolution, and their biogenesis is strictly controlled at multiple levels. Regulatory proteins, such as KSRP, modulate rates and timing of enzymatic reactions responsible for maturation of select miRNAs from their primary transcripts in response to specific stimuli. Here, we show that KSRP silencing in mesenchymal C2C12 cells produces a change in the transcriptome largely overlapping that induced by bone morphogenetic protein 2 (BMP2) signaling activation. This induces osteoblastic differentiation while preventing myogenic differentiation. KSRP silencing- and BMP2-dependent myogenic miRNA (myomiR) maturation blockade is required for osteoblastic differentiation of C2C12 cells. Our results demonstrate that phosphorylated R-SMAD proteins, the transducers of BMP2 signal, associate with phosphorylated KSRP and block its interaction with primary myomiRs. This abrogates KSRP-dependent myomiR maturation, with SMAD4, SMAD5, and SMAD9 silencing being able to rescue KSRP function. Thus, SMAD-induced blockade of KSRP-dependent myomiR maturation is critical for orienting C2C12 cell differentiation toward osteoblastic lineage.

Synthetic oligonucleotides recruit ILF2/3 to RNA transcripts to modulate splicing

We describe a new technology for recruiting specific proteins to RNA through selective recognition of heteroduplexes formed with chemically modified antisense oligonucleotides (ASOs). Typically, ASOs function by hybridizing to their RNA targets and blocking the binding of single-stranded RNA-binding proteins. Unexpectedly, we found that ASOs with 2'-deoxy-2'-fluoro (2'-F) nucleotides, but not with other 2' chemical modifications, have an additional property: they form heteroduplexes with RNA that are specifically recognized by the interleukin enhancer-binding factor 2 and 3 complex (ILF2/3). 2'-F ASO-directed recruitment of ILF2/3 to RNA can be harnessed to control gene expression by modulating alternative splicing of target transcripts. ILF2/3 recruitment to precursor mRNA near an exon results in omission of the exon from the mature mRNA, both in cell culture and in mice. We discuss the possibility of using chemically engineered ASOs that recruit specific proteins to modulate gene expression for therapeutic intervention.

From face to interface recognition: a differential geometric approach to distinguish DNA from RNA binding surfaces

Protein nucleic acid interactions play a critical role in all steps of the gene expression pathway. Nucleic acid (NA) binding proteins interact with their partners, DNA or RNA, via distinct regions on their surface that are characterized by an ensemble of chemical, physical and geometrical properties. In this study, we introduce a novel methodology based on differential geometry, commonly used in face recognition, to characterize and predict NA binding surfaces on proteins. Applying the method on experimentally solved three-dimensional structures of proteins we successfully classify double-stranded DNA (dsDNA) from single-stranded RNA (ssRNA) binding proteins, with 83% accuracy. We show that the method is insensitive to conformational changes that occur upon binding and can be applicable for de novo protein-function prediction. Remarkably, when concentrating on the zinc finger motif, we distinguish successfully between RNA and DNA binding interfaces possessing the same binding motif even within the same protein, as demonstrated for the RNA polymerase transcription-factor, TFIIIA. In conclusion, we present a novel methodology to characterize protein surfaces, which can accurately tell apart dsDNA from an ssRNA binding interfaces. The strength of our method in recognizing fine-tuned differences on NA binding interfaces make it applicable for many other molecular recognition problems, with potential implications for drug design.