Large favorable enthalpy changes drive specific RNA recognition by RNA recognition motif proteins

The RNA recognition motif (RRM) is a prevalent class of RNA binding domains. Although a number of RRM/RNA structures have been determined, thermodynamic analyses are relatively uncommon. Here, we use isothermal titration calorimetry to characterize single-stranded (ss)RNA binding by four representative RRM-containing proteins: (i) U2AF(65), (ii) SXL, (iii) TIA-1, and (iv) PAB. In all cases, ssRNA binding is accompanied by remarkably large favorable enthalpy changes (-30 to -60 kcal mol(-1)) and unfavorable entropy changes. Alterations of key RRM residues and binding sites indicate that under the nearly physiological conditions of these studies, large thermodynamic changes represent a signature of specific ssRNA recognition by RRMs.     

Surface shapes and surrounding environment analysis of single‐ and double‐stranded DNA‐binding proteins in protein‐DNA interface

Protein‐DNA bindings are critical to many biological processes. However, the structural mechanisms underlying these interactions are not fully understood. Here, we analyzed the residues shape (peak, flat, or valley) and the surrounding environment of double‐stranded DNA‐binding proteins (DSBs) and single‐stranded DNA‐binding proteins (SSBs) in protein‐DNA interfaces. In the results, we found that the interface shapes, hydrogen bonds, and the surrounding environment present significant differences between the two kinds of proteins. Built on the investigation results, we constructed a random forest (RF) classifier to distinguish DSBs and SSBs with satisfying performance. In conclusion, we present a novel methodology to characterize protein interfaces, which will deepen our understanding of the specificity of proteins binding to ssDNA (single‐stranded DNA) or dsDNA (double‐stranded DNA). Proteins 2016; 84:979–989. © 2016 Wiley Periodicals, Inc.     

Characterization of the RNA-binding domains in the replicase proteins of tomato bushy stunt virus

Tomato bushy stunt virus (TBSV), a tombusvirus with a nonsegmented, plus-stranded RNA genome, codes for two essential replicase proteins. The sequence of one of the replicase proteins, namely p33, overlaps with the N-terminal domain of p92, which contains the signature motifs of RNA-dependent RNA polymerases (RdRps) in its nonoverlapping C-terminal portion. In this work, we demonstrate that both replicase proteins bind to RNA in vitro based on gel mobility shift and surface plasmon resonance measurements. We also show evidence that the binding of p33 to single-stranded RNA (ssRNA) is stronger than binding to double-stranded RNA (dsRNA), ssDNA, or dsDNA in vitro. Competition experiments with ssRNA revealed that p33 binds to a TBSV-derived sequence with higher affinity than to other nonviral ssRNA sequences. Additional studies revealed that p33 could bind to RNA in a cooperative manner. Using deletion derivatives of the Escherichia coli-expressed recombinant proteins in gel mobility shift and Northwestern assays, we demonstrate that p33 and the overlapping domain of p92, based on its sequence identity with p33, contain an arginine- and proline-rich RNA-binding motif (termed RPR, which has the sequence RPRRRP). This motif is highly conserved among tombusviruses and related carmoviruses, and it is similar to the arginine-rich motif present in the Tat trans-activator protein of human immunodeficiency virus type 1. We also find that the nonoverlapping C-terminal domain of p92 contains additional RNA-binding regions. Interestingly, the location of one of the RNA-binding domains in p92 is similar to the RNA-binding domain of the NS5B RdRp protein of hepatitis C virus.     

The Fanconi anemia associated protein FAAP24 uses two substrate specific binding surfaces for DNA recognition

To maintain the integrity of the genome, multiple DNA repair systems exist to repair damaged DNA. Recognition of altered DNA, including bulky adducts, pyrimidine dimers and interstrand crosslinks (ICL), partially depends on proteins containing helix-hairpin-helix (HhH) domains. To understand how ICL is specifically recognized by the Fanconi anemia proteins FANCM and FAAP24, we determined the structure of the HhH domain of FAAP24. Although it resembles other HhH domains, the FAAP24 domain contains a canonical hairpin motif followed by distorted motif. The HhH domain can bind various DNA substrates; using nuclear magnetic resonance titration experiments, we demonstrate that the canonical HhH motif is required for double-stranded DNA (dsDNA) binding, whereas the unstructured N-terminus can interact with single-stranded DNA. Both DNA binding surfaces are used for binding to ICL-like single/double-strand junction-containing DNA substrates. A structural model for FAAP24 bound to dsDNA has been made based on homology with the translesion polymerase iota. Site-directed mutagenesis, sequence conservation and charge distribution support the dsDNA-binding model. Analogous to other HhH domain-containing proteins, we suggest that multiple FAAP24 regions together contribute to binding to single/double-strand junction, which could contribute to specificity in ICL DNA recognition.     

Structural determinants of human APOBEC3A enzymatic and nucleic acid binding properties

Human APOBEC3A (A3A) is a single-domain cytidine deaminase that converts deoxycytidine residues to deoxyuridine in single-stranded DNA (ssDNA). It inhibits a wide range of viruses and endogenous retroelements such as LINE-1, but it can also edit genomic DNA, which may play a role in carcinogenesis. Here, we extend our recent findings on the NMR structure of A3A and report structural, biochemical and cell-based mutagenesis studies to further characterize A3A’s deaminase and nucleic acid binding activities. We find that A3A binds ssRNA, but the RNA and DNA binding interfaces differ and no deamination of ssRNA is detected. Surprisingly, with only one exception (G105A), alanine substitution mutants with changes in residues affected by specific ssDNA binding retain deaminase activity. Furthermore, A3A binds and deaminates ssDNA in a length-dependent manner. Using catalytically active and inactive A3A mutants, we show that the determinants of A3A deaminase activity and anti-LINE-1 activity are not the same. Finally, we demonstrate A3A’s potential to mutate genomic DNA during transient strand separation and show that this process could be counteracted by ssDNA binding proteins. Taken together, our studies provide new insights into the molecular properties of A3A and its role in multiple cellular and antiviral functions.     

Characterization of a family of RanBP2-type zinc fingers that can recognize single-stranded RNA

The recognition of single-stranded RNA (ssRNA) is an important aspect of gene regulation, and a number of different classes of protein domains that recognize ssRNA in a sequence-specific manner have been identified. Recently, we demonstrated that the RanBP2-type zinc finger (ZnF) domains from the human splicing factor ZnF Ran binding domain-containing protein 2 (ZRANB2) can bind to a sequence containing the consensus AGGUAA. Six other human proteins, namely, Ewing's sarcoma (EWS), translocated in liposarcoma (TLS)/FUS, RNA-binding protein 56 (RBP56), RNA-binding motif 5 (RBM5), RNA-binding motif 10 (RBM10) and testis-expressed sequence 13A (TEX13A), each contains a single ZnF with homology to the ZRANB2 ZnFs, and several of these proteins have been implicated in the regulation of mRNA processing. Here, we show that all of these ZnFs are able to bind with micromolar affinities to ssRNA containing a GGU motif. NMR titration data reveal that binding is mediated by the corresponding surfaces on each ZnF, and we also show that sequence selectivity is largely limited to the GGU core motif and that substitution of the three flanking adenines that were selected in our original selection experiment has a minimal effect on binding affinity. These data establish a subset of RanBP2-type ZnFs as a new family of ssRNA-binding motifs.     

Predicting nucleic acid binding interfaces from structural models of proteins

The function of DNA‐ and RNA‐binding proteins can be inferred from the characterization and accurate prediction of their binding interfaces. However, the main pitfall of various structure‐based methods for predicting nucleic acid binding function is that they are all limited to a relatively small number of proteins for which high‐resolution three‐dimensional structures are available. In this study, we developed a pipeline for extracting functional electrostatic patches from surfaces of protein structural models, obtained using the I‐TASSER protein structure predictor. The largest positive patches are extracted from the protein surface using the patchfinder algorithm. We show that functional electrostatic patches extracted from an ensemble of structural models highly overlap the patches extracted from high‐resolution structures. Furthermore, by testing our pipeline on a set of 55 known nucleic acid binding proteins for which I‐TASSER produces high‐quality models, we show that the method accurately identifies the nucleic acids binding interface on structural models of proteins. Employing a combined patch approach we show that patches extracted from an ensemble of models better predicts the real nucleic acid binding interfaces compared with patches extracted from independent models. Overall, these results suggest that combining information from a collection of low‐resolution structural models could be a valuable approach for functional annotation. We suggest that our method will be further applicable for predicting other functional surfaces of proteins with unknown structure. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

MuMoD: a Bayesian approach to detect multiple modes of protein–DNA binding from genome-wide ChIP data

High-throughput chromatin immunoprecipitation has become the method of choice for identifying genomic regions bound by a protein. Such regions are then investigated for overrepresented sequence motifs, the assumption being that they must correspond to the binding specificity of the profiled protein. However this approach often fails: many bound regions do not contain the ‘expected’ motif. This is because binding DNA directly at its recognition site is not the only way the protein can cause the region to immunoprecipitate. Its binding specificity can change through association with different co-factors, it can bind DNA indirectly, through intermediaries, or even enforce its function through long-range chromosomal interactions. Conventional motif discovery methods, though largely capable of identifying overrepresented motifs from bound regions, lack the ability to characterize such diverse modes of protein–DNA binding and binding specificities. We present a novel Bayesian method that identifies distinct protein–DNA binding mechanisms without relying on any motif database. The method successfully identifies co-factors of proteins that do not bind DNA directly, such as mediator and p300. It also predicts literature-supported enhancer–promoter interactions. Even for well-studied direct-binding proteins, this method provides compelling evidence for previously uncharacterized dependencies within positions of binding sites, long-range chromosomal interactions and dimerization.     

Binding Site Identification and Flexible Docking of Single Stranded RNA to Proteins Using a Fragment-Based Approach

Protein-RNA docking is hampered by the high flexibility of RNA, and particularly single-stranded RNA (ssRNA). Yet, ssRNA regions typically carry the specificity of protein recognition. The lack of methodology for modeling such regions limits the accuracy of current protein-RNA docking methods. We developed a fragment-based approach to model protein-bound ssRNA, based on the structure of the protein and the sequence of the RNA, without any prior knowledge of the RNA binding site or the RNA structure. The conformational diversity of each fragment is sampled by an exhaustive RNA fragment library that was created from all the existing experimental structures of protein-ssRNA complexes. A systematic and detailed analysis of fragment-based ssRNA docking was performed which constitutes a proof-of-principle for the fragment-based approach. The method was tested on two 8-homo-nucleotide ssRNA-protein complexes and was able to identify the binding site on the protein within 10 Å. Moreover, a structure of each bound ssRNA could be generated in close agreement with the crystal structure with a mean deviation of ~1.5 Å except for a terminal nucleotide. This is the first time a bound ssRNA could be modeled from sequence with high precision.     

Fragment-based modelling of single stranded RNA bound to RNA recognition motif containing proteins

Protein-RNA complexes are important for many biological processes. However, structural modeling of such complexes is hampered by the high flexibility of RNA. Particularly challenging is the docking of single-stranded RNA (ssRNA). We have developed a fragment-based approach to model the structure of ssRNA bound to a protein, based on only the protein structure, the RNA sequence and conserved contacts. The conformational diversity of each RNA fragment is sampled by an exhaustive library of trinucleotides extracted from all known experimental protein–RNA complexes. The method was applied to ssRNA with up to 12 nucleotides which bind to dimers of the RNA recognition motifs (RRMs), a highly abundant eukaryotic RNA-binding domain. The fragment based docking allows a precise de novo atomic modeling of protein-bound ssRNA chains. On a benchmark of seven experimental ssRNA–RRM complexes, near-native models (with a mean heavy-atom deviation of <3 Å from experiment) were generated for six out of seven bound RNA chains, and even more precise models (deviation < 2 Å) were obtained for five out of seven cases, a significant improvement compared to the state of the art. The method is not restricted to RRMs but was also successfully applied to Pumilio RNA binding proteins.     

Molecular basis of double-stranded RNA-protein interactions: structure of a dsRNA-binding domain complexed with dsRNA.

Protein interactions with double-stranded RNA (dsRNA) are critical for many cell processes; however, in contrast to protein-dsDNA interactions, surprisingly little is known about the molecular basis of protein-dsRNA interactions. A large and diverse class of proteins that bind dsRNA do so by utilizing an approximately 70 amino acid motif referred to as the dsRNA-binding domain (dsRBD). We have determined a 1.9 A resolution crystal structure of the second dsRBD of Xenopus laevis RNA-binding protein A complexed with dsRNA. The structure shows that the protein spans 16 bp of dsRNA, interacting with two successive minor grooves and across the intervening major groove on one face of a primarily A-form RNA helix. The nature of these interactions explains dsRBD specificity for dsRNA (over ssRNA or dsDNA) and the apparent lack of sequence specificity. Interestingly, the dsRBD fold resembles a portion of the conserved core structure of a family of polynucleotidyl transferases that includes RuvC, MuA transposase, retroviral integrase and RNase H. Structural comparisons of the dsRBD-dsRNA complex and models proposed for polynucleotidyl transferase-nucleic acid complexes suggest that similarities in nucleic acid binding also exist between these families of proteins.     

RNA recognition by the DNA end-binding Ku heterodimer

Most nucleic acid-binding proteins selectively bind either DNA or RNA, but not both nucleic acids. The Saccharomyces cerevisiae Ku heterodimer is unusual in that it has two very different biologically relevant binding modes: (1) Ku is a sequence-nonspecific double-stranded DNA end-binding protein with prominent roles in nonhomologous end-joining and telomeric capping, and (2) Ku associates with a specific stem–loop of TLC1, the RNA subunit of budding yeast telomerase, and is necessary for proper nuclear localization of this ribonucleoprotein enzyme. TLC1 RNA-binding and dsDNA-binding are mutually exclusive, so they may be mediated by the same site on Ku. Although dsDNA binding by Ku is well studied, much less is known about what features of an RNA hairpin enable specific recognition by Ku. To address this question, we localized the Ku-binding site of the TLC1 hairpin with single-nucleotide resolution using phosphorothioate footprinting, used chemical modification to identify an unpredicted motif within the hairpin secondary structure, and carried out mutagenesis of the stem–loop to ascertain the critical elements within the RNA that permit Ku binding. Finally, we provide evidence that the Ku-binding site is present in additional budding yeast telomerase RNAs and discuss the possibility that RNA binding is a conserved function of the Ku heterodimer.     

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.     

PARP activation regulates the RNA-binding protein NONO in the DNA damage response to DNA double-strand breaks

After the generation of DNA double-strand breaks (DSBs), poly(ADP-ribose) polymerase-1 (PARP-1) is one of the first proteins to be recruited and activated through its binding to the free DNA ends. Upon activation, PARP-1 uses NAD⁺ to generate large amounts of poly(ADP-ribose) (PAR), which facilitates the recruitment of DNA repair factors. Here, we identify the RNA-binding protein NONO, a partner protein of SFPQ, as a novel PAR-binding protein. The protein motif being primarily responsible for PAR-binding is the RNA recognition motif 1 (RRM1), which is also crucial for RNA-binding, highlighting a competition between RNA and PAR as they share the same binding site. Strikingly, the in vivo recruitment of NONO to DNA damage sites completely depends on PAR, generated by activated PARP-1. Furthermore, we show that upon PAR-dependent recruitment, NONO stimulates nonhomologous end joining (NHEJ) and represses homologous recombination (HR) in vivo. Our results therefore place NONO after PARP activation in the context of DNA DSB repair pathway decision. Understanding the mechanism of action of proteins that act in the same pathway as PARP-1 is crucial to shed more light onto the effect of interference on PAR-mediated pathways with PARP inhibitors, which have already reached phase III clinical trials but are until date poorly understood.     

Both helix topology and counterion distribution contribute to the more effective charge screening in dsRNA compared with dsDNA

The recent discovery of the RNA interference mechanism emphasizes the biological importance of short, isolated, double-stranded (ds) RNA helices and calls for a complete understanding of the biophysical properties of dsRNA. However, most previous studies of the electrostatics of nucleic acid duplexes have focused on DNA. Here, we present a comparative investigation of electrostatic effects in RNA and DNA. Using resonant (anomalous) and non-resonant small-angle X-ray scattering, we characterized the charge screening efficiency and counterion distribution around short (25 bp) dsDNA and RNA molecules of comparable sequence. Consistent with theoretical predictions, we find counterion mediated screening to be more efficient for dsRNA than dsDNA. Furthermore, the topology of the RNA A-form helix alters the spatial distribution of counterions relative to B-form DNA. The experimental results reported here agree well with ion-size-corrected non-linear Poisson–Boltzmann calculations. We propose that differences in electrostatic properties aid in selective recognition of different types of short nucleic acid helices by target binding partners.