Human hnRNP protein A1: A model polypeptide for a structural and genetic investigation of a broad family of RNA binding proteins

The hnRNP fiber is the substrate on which pre-mRNA processing occurs. The protein moiety of the fiber (hnRNP proteins) constitutes a broad family of RNA binding proteins that revealed, upon molecular analysis, a number of interesting features.Heterogeneous nuclear ribonucleoprotein A1 is a major component of the human hnRNP complex. In recent years this protein has attracted great attention because of several emerging evidences of its direct involvement in pre-mRNA processing and it has become one of the best characterized RNA binding proteins. Detailed knowledge of the structure of protein A1 has laid the basis for the understanding of its function, and for this reason A1 can be considered as a model polypeptide for the investigation of a large number of RNA binding proteins.In this work we report recent findings regarding the binding properties of protein A1 as well as new data on the gene structure of A1 and of its closely related hnRNP protein A2. Our results show that a single A1 molecule contains the determinants for simultaneous binding of two single-stranded nucleic acid molecules and we demonstrate that the glycine-rich domain of A1, isolated from the rest of the molecule, is capable of sustaining protein-protein interactions. These features probably account for the reannealing activity of the protein and for its capacity to modulate the binding of snRNPs to intron sequencesin vitro. Comparison of A1 and A2 gene sequences revealed a remarkable conservation of the overall structural organization, suggesting important functions for the different structural elements.     

Evolutionary applications of MIRs and SINEs

It is believed that short interspersed elements (SINEs) are irreversibly inserted into genomes. We use this concept to try to deduce the evolution of whales using sequence and hybridization studies. The observation that microsatellites are associated with SINEs lead us to screen sequences surrounding cetacean microsatellites for artiodactyl-derived SINEs. Two sequences that were thought to be cetacean SINEs and the bovine SINE were aligned for comparison to sequences flanking microsatellites from ungulates. The bovine SINE was observed only in ruminants while CetSINE1 and 2 were found in mammals. Hybridization studies using these three SINEs revealed that CetSINE1 was found in all ungulates and cetaceans with the strongest hybridization signal observed in the hippopotamus and beluga; CetSINE2 hybridized to all ungulate suborders, while the bovine SINE was only observed in Ruminantia. It is proposed that these putative SINEs are not 'generic' SINEs but mammalian-wide interspersed repeats (MIRs). Caution is urged: what initially appears to be a SINE may instead be a MIR and have reduced evolutionary resolving power.     

Prediction and validation of the unexplored RNA‐binding protein atlas of the human proteome

Detecting protein‐RNA interactions is challenging both experimentally and computationally because RNAs are large in number, diverse in cellular location and function, and flexible in structure. As a result, many RNA‐binding proteins (RBPs) remain to be identified. Here, a template‐based, function‐prediction technique SPOT‐Seq for RBPs is applied to human proteome and its result is validated by a recent proteomic experimental discovery of 860 mRNA‐binding proteins (mRBPs). The coverage (or sensitivity) is 42.6% for 1217 known RBPs annotated in the Gene Ontology and 43.6% for 860 newly discovered human mRBPs. Consistent sensitivity indicates the robust performance of SPOT‐Seq for predicting RBPs. More importantly, SPOT‐Seq detects 2418 novel RBPs in human proteome, 291 of which were validated by the newly discovered mRBP set. Among 291 validated novel RBPs, 61 are not homologous to any known RBPs. Successful validation of predicted novel RBPs permits us to further analysis of their phenotypic roles in disease pathways. The dataset of 2418 predicted novel RBPs along with confidence levels and complex structures is available at http://sparks-lab.org (in publications) for experimental confirmations and hypothesis generation. Proteins 2014; 82:640–647. © 2013 Wiley Periodicals, Inc.     

The function of proteins that interact with mRNA

Specific proteins are associated with mRNA in the cytoplasm of eukaryotic cells. The complement of associated proteins depends upon whether the mRNA is an integral component of the polysomal complex being translated, or, alternatively, whether it is part of the non-translated free mRNP fraction. By subjecting cells to ultraviolet irradiation in vivo to cross-link proteins to mRNA, mRNP proteins have been shown to be associated with specific regions of the mRNA molecule. Examination of mRNP complexes containing a unique mRNA has suggested that not all mRNA contain the same family of associated RNA binding proteins. The function of mRNA associated proteins may include a role in providing stability for mRNA, and/or in modulating translation. With the recent demonstrations that both free and polysomal mRNPs are associated with the cytoskeletal framework, specific mRNP proteins may play a role in determining the subcellular localization of specific mRNPs.     

A predominantly nuclear protein affecting cytoplasmic localization of beta-actin mRNA in fibroblasts and neurons.

The localization of beta-actin mRNA to the leading lamellae of chicken fibroblasts and neurite growth cones of developing neurons requires a 54-nt localization signal (the zipcode) within the 3' untranslated region. In this study we have identified and isolated five proteins binding to the zipcode. One of these we previously identified as zipcode binding protein (ZBP)1, a 4-KH domain protein. A second is now investigated in detail: a 92-kD protein, ZBP2, that is especially abundant in extracts from embryonic brain. We show that ZBP2 is a homologue of the human hnRNP protein, KSRP, that appears to mediate pre-mRNA splicing. However, ZBP2 has a 47-amino acid (aa) sequence not present in KSRP. Various portions of ZBP2 fused to GFP indicate that the protein most likely shuttles between the nucleus and the cytoplasm, and that the 47-aa insert promotes the nuclear localization. Expression of a truncated ZBP2 inhibits the localization of beta-actin mRNA in both fibroblast and neurons. These data suggest that ZBP2, although predominantly a nuclear protein, has a role in the cytoplasmic localization of beta-actin mRNA.     

Evolutionary modes of emergence of short interspersed nuclear element (SINE) families in grasses

Short interspersed nuclear elements (SINEs) are non‐autonomous transposable elements which are propagated by retrotransposition and constitute an inherent part of the genome of most eukaryotic species. Knowledge of heterogeneous and highly abundant SINEs is crucial for de novo (or improvement of) annotation of whole genome sequences. We scanned Poaceae genome sequences of six important cereals (Oryza sativa, Triticum aestivum, Hordeum vulgare, Panicum virgatum, Sorghum bicolor, Zea mays) and Brachypodium distachyon to examine the diversity and evolution of SINE populations. We comparatively analyzed the structural features, distribution, evolutionary relation and abundance of 32 SINE families and subfamilies within grasses, comprising 11 052 individual copies. The investigation of activity profiles within the Poaceae provides insights into their species‐specific diversification and amplification. We found that Poaceae SINEs (PoaS) fall into two length categories: simple SINEs of up to 180 bp and dimeric SINEs larger than 240 bp. Detailed analysis at the nucleotide level revealed that multimerization of related and unrelated SINE copies is an important evolutionary mechanism of SINE formation. We conclude that PoaS families diversify by massive reshuffling between SINE families, likely caused by insertion of truncated copies, and provide a model for this evolutionary scenario. Twenty‐eight of 32 PoaS families and subfamilies show significant conservation, in particular either in the 5′ or 3′ regions, across Poaceae species and share large sequence stretches with one or more other PoaS families.     

Additional box B of RNA polymerase III promoter in SINE B1 can be functional

Many genes of small RNAs and short interspersed elements (SINEs) are transcribed by RNA polymerase III due to an internal promoter that is composed of two boxes (A and B) spaced by 30-45 bp. Rodent SINE B1 originated from 7SL RNA, and a 29-bp tandem duplication took place in B1 at an early stage of its evolution. As a result of this duplication, an additional box B (named B') located at a distance of 79-82 bp from box A arose in SINE B1. Here we have shown that despite the unusually large distance between boxes A and B', they can form an active promoter. In chinchillas, guinea pigs, and other rodents belonging to clade Ctenohystrica, structure of the B' box was well preserved and closely resembles the canonical B box. One may suggest therefore, that box B' can functionally replace box B in those copies of B1 where the latter has lost activity due to mutations.     

Divergence of 3′ ends as a driver of short interspersed nuclear element (SINE) evolution in the Salicaceae

Short interspersed nuclear elements (SINEs) are small, non‐autonomous and heterogeneous retrotransposons that are widespread in plants. To explore the amplification dynamics and evolutionary history of SINE populations in representative deciduous tree species, we analyzed the genomes of the six following Salicaceae species: Populus deltoides, Populus euphratica, Populus tremula, Populus tremuloides, Populus trichocarpa, and Salix purpurea. We identified 11 Salicaceae SINE families (SaliS‐I to SaliS‐XI), comprising 27 077 full‐length copies. Most of these families harbor segmental similarities, providing evidence for SINE emergence by reshuffling or heterodimerization. We observed two SINE groups, differing in phylogenetic distribution pattern, similarity and 3′ end structure. These groups probably emerged during the ‘salicoid duplication’ (~65 million years ago) in the Salix–Populus progenitor and during the separation of the genus Salix (45–65 million years ago), respectively. In contrast to conserved 5′ start motifs across species and SINE families, the 3′ ends are highly variable in sequence and length. This extraordinary 3′‐end variability results from mutations in the poly(A) tail, which were fixed by subsequent amplificational bursts. We show that the dissemination of newly evolved 3′ ends is accomplished by a displacement of older motifs, leading to various 3′‐end subpopulations within the SaliS families.     

Evaluating conformational changes in protein structures binding RNA

Many protein-RNA recognition events are known to exhibit conformational changes from qualitative observations of individual complexes. However, a quantitative estimation of conformational changes is required if protein-RNA docking and template-based methods for RNA binding site prediction are to be developed. This study presents the first quantitative evaluation of conformational changes that occur when proteins bind RNA. The analysis of twelve RNA-binding proteins in the bound and unbound states using error-scaled difference distance matrices is presented. The binding site residues are mapped to each structure, and the conformational changes that affect these residues are evaluated. Of the twelve proteins four exhibit greater movements in nonbinding site residues, and a further four show the greatest movements in binding site residues. The remaining four proteins display no significant conformational change. When interface residues are found to be in conformationally variable regions of the protein they are typically seen to move less than 2 A between the bound and unbound conformations. The current data indicate that conformational changes in the binding site residues of RNA binding proteins may not be as significant as previously suggested, but a larger data set is required before wider conclusions may be drawn. The implications of the observed conformational changes for protein function prediction are discussed.     

Single nucleotide polymorphisms (SNPs) within Bov-A2 SINE in the second intron of bovine and buffalo k-casein (CSN3) gene

The genetic polymorphism of an entire Bov-A2 element located in the second intron of the buffalo and bovine k-casein (CSN3) gene was investigated by amplification and sequencing of PCR products. Single nucleotide polymorphisms were detected. A PCR-RFLP method was developed to detect an A or G mutation at position 605 of bovine Bov-A2 element which creates a BfaI polymorphic site. The frequencies of the B allele, with the BfaI site, were for 0.275, 0.775, 0.750, 0.975, respectively, for Italian Holstein Friesian, Grey Alpine, Friuli Red Pied and Reggio bovine breeds. The mutation rate (substitutions and deletions/insertions per nucleotide site per year) was 2.5 x 10(-9) for Bov-A2 sequences in the second intron of CSN3. The comparison with other Bov-A2 elements suggests that this retroelement might be an important source of single nucleotide polymorphism for analysis of Bovidae genomes.     

A new size-independent score for pairwise protein structure alignment and its application to structure classification and nucleic-acid binding prediction

A structure alignment program aligns two structures by optimizing a scoring function that measures structural similarity. It is highly desirable that such scoring function is independent of the sizes of proteins in comparison so that the significance of alignment across different sizes of the protein regions aligned is comparable. Here, we developed a new score called SP‐score that fixes the cutoff distance at 4 Å and removed the size dependence using a normalization prefactor. We further built a program called SPalign that optimizes SP‐score for structure alignment. SPalign was applied to recognize proteins within the same structure fold and having the same function of DNA or RNA binding. For fold discrimination, SPalign improves sensitivity over TMalign for the chain‐level comparison by 12% and over DALI for the domain‐level comparison by 13% at the same specificity of 99.6%. The difference between TMalign and SPalign at the chain level is due to the inability of TMalign to detect single domain similarity between multidomain proteins. For recognizing nucleic acid binding proteins, SPalign consistently improves over TMalign by 12% and DALI by 31% in average value of Mathews correlation coefficients for four datasets. SPalign with default setting is 14% faster than TMalign. SPalign is expected to be useful for function prediction and comparing structures with or without domains defined. The source code for SPalign and the server are available at http://sparks.informatics.iupui.edu . Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Identification of FUSE-binding proteins as interacting partners of TIA proteins

TIA-1 and TIAR are closely related RNA-binding proteins involved in several mechanisms of RNA metabolism, including alternative hnRNA splicing and mRNA translation regulation. In particular, TIA-1 represses tumor necrosis factor (TNF) mRNA translation by binding to the AU-rich element (ARE) present in the mRNA 3' untranslated region. Here, we demonstrate that TIA proteins interact with FUSE-binding proteins (FBPs) and that fbp genes are co-expressed with tia genes during Xenopus embryogenesis. FBPs participate in various steps of RNA processing and degradation. In Cos cells, FBPs co-localize with TIA proteins in the nucleus and migrate into TIA-enriched cytoplasmic granules upon oxidative stress. Overexpression of FBP2-KH3 RNA-binding domain fused to EGFP induces the specific sequestration of TIA proteins in cytoplasmic foci, thereby precluding their nuclear accumulation. In cytosolic RAW 264.7 macrophage extracts, FBPs are found associated in EMSA to the TIA-1/TNF-ARE complex. Together, our results indicate that TIA and FBP proteins may thus be relevant biological involved in common events of RNA metabolism occurring both in the nucleus and the cytoplasm.     

Divergent RNA binding specificity of yeast Puf2p

PUF proteins bind mRNAs and regulate their translation, stability, and localization. Each PUF protein binds a selective group of mRNAs, enabling their coordinate control. We focus here on the specificity of Puf2p and Puf1p of Saccharomyces cerevisiae, which copurify with overlapping groups of mRNAs. We applied an RNA-adapted version of the DRIM algorithm to identify putative binding sequences for both proteins. We first identified a novel motif in the 3' UTRs of mRNAs previously shown to associate with Puf2p. This motif consisted of two UAAU tetranucleotides separated by a 3-nt linker sequence, which we refer to as the dual UAAU motif. The dual UAAU motif was necessary for binding to Puf2p, as judged by gel shift, yeast three-hybrid, and coimmunoprecipitation from yeast lysates. The UAAU tetranucleotides are required for optimal binding, while the identity and length of the linker sequences are less critical. Puf1p also binds the dual UAAU sequence, consistent with the prior observation that it associates with similar populations of mRNAs. In contrast, three other canonical yeast PUF proteins fail to bind the Puf2p recognition site. The dual UAAU motif is distinct from previously known PUF protein binding sites, which invariably possess a UGU trinucleotide. This study expands the repertoire of cis elements bound by PUF proteins and suggests new modes by which PUF proteins recognize their mRNA targets.