Contributions of the RNA-binding and linker domains and RNA structure to the specificity and affinity of the nucleolin RBD12/NRE interaction

Nucleolin is a multidomain phosphoprotein involved in ribosome biogenesis. In vitro selection and binding studies with pre-rRNA fragments have shown that the first two RNA-binding domains (RBDs) in nucleolin (RBD12) recognize the consensus sequence (U/G)CCCG(A/G) in the context of a stem-loop structure (nucleolin-recognition element = NRE). Structural studies of nucleolin RBD12 in complex with an in vitro selected NRE (sNRE) and a natural pre-rRNA NRE (b2NRE) have revealed that sequence-specific binding of the consensus NRE is achieved in a similar manner in both complexes using residues in both RBDs as well as the linker connecting them. Using fluorescence anisotropy (FA) and nuclear magnetic resonance (NMR), we demonstrate the importance of the linker for NRE affinity by showing that only the individual RBDs with the linker attached retain the ability to specifically bind, albeit weakly, to sNRE and b2NRE. Binding of RBD1 and RBD2 to the NREs in trans is not detected even when one of the RBDs has the linker attached, which suggests that the linker also contributes to the affinity by tethering the two RBDs. To determine if binding of nucleolin RBD12 to natural NREs is dependent on a specific RNA stem-loop structure, as was the case for the sNRE, we conducted FA and NMR binding assays with nucleolin RBD12 and a single-stranded NRE. The results show that nucleolin RBD12 sequence-specifically binds a single-stranded NRE with an affinity similar to that for b2NRE, indicating that a stem-loop structure is not required for the nucleolin RBD12/pre-rRNA NRE interaction.     

Solution structures of stem-loop RNAs that bind to the two N-terminal RNA-binding domains of nucleolin

Nucleolin, a multi-domain protein involved in ribosome biogenesis, has been shown to bind the consensus sequence (U/G)CCCG(A/G) in the context of a hairpin loop structure (nucleolin recognition element; NRE). Previous studies have shown that the first two RNA-binding domains in nucleolin (RBD12) are responsible for the interaction with the in vitro selected NRE (sNRE). We have previously reported the structures of nucleolin RBD12, sNRE and nucleolin RBD12–sNRE complex. A comparison of free and bound sNRE shows that the NRE loop becomes structured upon binding. From this observation, we hypothesized that the disordered hairpin loop of sNRE facilitates conformational rearrangements when the protein binds. Here, we show that nucleolin RBD12 is also sufficient for sequence- specific binding of two NRE sequences found in pre-rRNA, b1NRE and b2NRE. Structural investigations of the free NREs using NMR spectroscopy show that the b1NRE loop is conformationally heterogeneous, while the b2NRE loop is structured. The b2NRE forms a hairpin capped by a YNMG-like tetraloop. Comparison of the chemical shifts of sNRE and b2NRE in complex with nucleolin RBD12 suggests that the NRE consensus nucleotides adopt a similar conformation. These results show that a disordered NRE consensus sequence is not a prerequisite for nucleolin RBD12 binding.     

Molecular basis of sequence-specific recognition of pre-ribosomal RNA by nucleolin

The structure of the 28 kDa complex of the first two RNA binding domains (RBDs) of nucleolin (RBD12) with an RNA stem-loop that includes the nucleolin recognition element UCCCGA in the loop was determined by NMR spectroscopy. The structure of nucleolin RBD12 with the nucleolin recognition element (NRE) reveals that the two RBDs bind on opposite sides of the RNA loop, forming a molecular clamp that brings the 5' and 3' ends of the recognition sequence close together and stabilizing the stem-loop. The specific interactions observed in the structure explain the sequence specificity for the NRE sequence. Binding studies of mutant proteins and analysis of conserved residues support the proposed interactions. The mode of interaction of the protein with the RNA and the location of the putative NRE sites suggest that nucleolin may function as an RNA chaperone to prevent improper folding of the nascent pre-rRNA.     

RNA recognition by the joint action of two nucleolin RNA-binding domains: genetic analysis and structural modeling.

The interaction of nucleolin with a short stem-loop structure (NRE) requires two contiguous RNA-binding domains (RBD 1+2). The structural basis for RNA recognition by these RBDs was studied using a genetic system in Escherichia coli. Within each of the two domains, we identified several mutations that severely impair interaction with the RNA target. Mutations that alter RNA-binding specificity were also isolated, suggesting the identity of specific contacts between RBD 1+2 amino acids and nucleotides within the NRE stem-loop. Our data indicate that both RBDs participate in a joint interaction with the NRE and that each domain uses a different surface to contact the RNA. The constraints provided by these genetic data and previous mutational studies have enabled us to propose a three-dimensional model of nucleolin RBD 1+2 bound to the NRE stem-loop.     

Recognition of pre-formed and flexible elements of an RNA stem-loop by nucleolin

Nucleolin is an abundant nucleolar protein which is essential for ribosome biogenesis. The first two of its four tandem RNA-binding domains (RBD12) specifically recognize a stem-loop structure containing a conserved UCCCGA sequence in the loop called the nucleolin-recognition element (NRE). We have determined the structure of the consensus SELEX NRE (sNRE) by NMR spectroscopy. In both the free and bound RNA the top part of the stem forms a loop E (or S-turn) motif. In the absence of protein, the structure of the hairpin loop is not well defined due to conformational heterogeneity, and appears to be in equilibrium between two families of conformations. Titrations of RBD1, RBD2, and RBD12 with the sNRE show that specific binding requires RBD12. In complex with RBD12, the hairpin loop interacts specifically with the protein and adopts a well-defined structure which shares some of the features of the free form. The loop E motif also has specific interactions with the protein. Implications of these findings for the mechanism of recognition of RNA structures by modular proteins are discussed.     

Two different combinations of RNA-binding domains determine the RNA binding specificity of nucleolin

Nucleolin is an abundant nucleolar protein involved in several steps of ribosome biogenesis. The protein is highly conserved through evolution and possesses four RNA-binding domains (RBD), which are likely to determine its RNA binding specificity. Previous studies have shown that nucleolin interacts with two different RNA targets. The first is a small stem-loop structure, the nucleolin recognition element (NRE), found all along the pre-ribosomal RNA. The second is a short single-stranded RNA sequence, the evolutionary conserved motif (ECM), located five nucleotides downstream of the first processing site in the pre-ribosomal RNA 5' external transcribed spacer. Biochemical, genetic, and structural studies have shown that the first two RBD of nucleolin are necessary and sufficient for the specific interaction of nucleolin with the NRE motif. In this work, we have studied the interaction of nucleolin with the ECM sequence. Deletion and mutational analyses showed that all four RBDs of hamster nucleolin were required for the interaction with the ECM sequence. This RNA binding specificity is conserved between hamster and Xenopus laevis, whereas the Xenopus protein does not interact with the NRE. Nucleolin is the first example of a protein that requires four RBDs for its interaction with an RNA target, demonstrating that a single protein can use different combinations of RBD to interact specifically with several RNA sequences.     

Solution structure of the two N-terminal RNA-binding domains of nucleolin and NMR study of the interaction with its RNA target

Nucleolin is an abundant 70 kDa nucleolar protein involved in many aspects of ribosomal RNA biogenesis. The central region of nucleolin contains four tandem consensus RNA-binding domains (RBD). The two most N-terminal domains (RBD12) bind with nanomolar affinity to an RNA stem-loop containing the consensus sequence UCCCGA in the loop. We have determined the solution structure of nucleolin RBD12 in its free form and have studied its interaction with a 22 nt RNA stem-loop using multidimensional NMR spectroscopy. The two RBDs adopt the expected beta alpha beta beta alpha beta fold, but the position of the beta 2 strand in both domains differs from what was predicted from sequence alignments. RBD1 and RBD2 are significantly different from each others and this is likely important in their sequence specific recognition of the RNA. RBD1 has a longer alpha-helix 1 and a shorter beta 2-beta 3 loop than RBD2, and differs from most other RBDs in these respects. The two RBDs are separated by a 12 amino acid flexible linker and do not interact with one another in the free protein. This linker becomes ordered when RBD12 binds to the RNA. Analysis of the observed NOEs between the protein and the RNA indicates that both RBDs interact with the RNA loop via their beta-sheet. Each domain binds residues on one side of the loop; specifically, RBD2 contacts the 5' side and RBD1 contacts the 3'.     

Evolutionary Conservation of the Ribosomal Biogenesis Factor Rbm19/Mrd1: Implications for Function

Ribosome biogenesis in eukaryotes requires coordinated folding and assembly of a pre-rRNA into sequential pre-rRNA-protein complexes in which chemical modifications and RNA cleavages occur. These processes require many small nucleolar RNAs (snoRNAs) and proteins. Rbm19/Mrd1 is one such protein that is built from multiple RNA-binding domains (RBDs). We find that Rbm19/Mrd1 with five RBDs is present in all branches of the eukaryotic phylogenetic tree, except in animals and Choanoflagellates, that instead have a version with six RBDs and Microsporidia which have a minimal Rbm19/Mrd1 protein with four RBDs. Rbm19/Mrd1 therefore evolved as a multi-RBD protein very early in eukaryotes. The linkers between the RBDs have conserved properties; they are disordered, except for linker 3, and position the RBDs at conserved relative distances from each other. All but one of the RBDs have conserved properties for RNA-binding and each RBD has a specific consensus sequence and a conserved position in the protein, suggesting a functionally important modular design. The patterns of evolutionary conservation provide information for experimental analyses of the function of Rbm19/Mrd1. In vivo mutational analysis confirmed that a highly conserved loop 5-β4-strand in RBD6 is essential for function.     

Multiple RNA interactions position Mrd1 at the site of the small subunit pseudoknot within the 90S pre-ribosome

Ribosomal subunit biogenesis in eukaryotes is a complex multistep process. Mrd1 is an essential and conserved small (40S) ribosomal subunit synthesis factor that is required for early cleavages in the 35S pre-ribosomal RNA (rRNA). Yeast Mrd1 contains five RNA-binding domains (RBDs), all of which are necessary for optimal function of the protein. Proteomic data showed that Mrd1 is part of the early pre-ribosomal complexes, and deletion of individual RBDs perturbs the pre-ribosomal structure. In vivo ultraviolet cross-linking showed that Mrd1 binds to the pre-rRNA at two sites within the 18S region, in helix 27 (h27) and helix 28. The major binding site lies in h27, and mutational analyses shows that this interaction requires the RBD1-3 region of Mrd1. RBD2 plays the dominant role in h27 binding, but other RBDs also contribute directly. h27 and helix 28 are located close to the sequences that form the central pseudoknot, a key structural feature of the mature 40S subunit. We speculate that the modular structure of Mrd1 coordinates pseudoknot formation with pre-rRNA processing and subunit assembly.     

Structure and interactions with RNA of the N-terminal UUAG-specific RNA-binding domain of hnRNP D0

Heterogeneous nuclear ribonucleoprotein (hnRNP) D0 has two ribonucleoprotein (RNP)-type RNA-binding domains (RBDs), each of which can bind solely to the UUAG sequence specifically. The structure of the N-terminal RBD (RBD1) determined by NMR is presented here. It folds into a compact alphabeta structure comprising a four-stranded antiparallel beta-sheet packed against two alpha-helices, which is characteristic of the RNP-type RBDs. Special structural features of RBD1 include N-capping boxes for both alpha-helices, a beta-bulge in the second beta-strand, and an additional short antiparallel beta-sheet coupled with a beta-turn-like structure in a loop. Two hydrogen bonds which restrict the positions of loops were identified. Backbone resonance assignments for RBD1 complexed with r(UUAGGG) revealed that the overall folding is maintained in the complex. The candidate residues involved in the interactions with RNA were identified by chemical shift perturbation analysis. They are located in the central and peripheral regions of the RNA-binding surface composed of the four-stranded beta-sheet, loops, and the C-terminal region. It is suggested that non-specific interactions with RNA are performed by the residues in the central region of the RNA-binding surface, while specific interactions are performed by those in the peripheral regions. It was also found that RBD1 has the ability to inhibit the formation of the quadruplex structure.     

Structure, backbone dynamics and interactions with RNA of the C-terminal RNA-binding domain of a mouse neural RNA-binding protein, Musashi1

Musashi1 is an RNA-binding protein abundantly expressed in the developing mouse central nervous system. Its restricted expression in neural precursor cells suggests that it is involved in the regulation of asymmetric cell division. Musashi1 contains two ribonucleoprotein (RNP)-type RNA-binding domains (RBDs), RBD1 and RBD2. Our previous studies showed that RBD1 alone binds to RNA, while the binding of RBD2 is not detected under the same conditions. Joining of RBD2 to RBD1, however, increases the affinity to greater than that of RBD1 alone, indicating that RBD2 contributes to RNA-binding. We have determined the three-dimensional solution structure of the C-terminal RBD (RBD2) of Musashi1 by NMR. It folds into a compact alpha beta structure comprising a four-stranded antiparallel beta-sheet packed against two alpha-helices, which is characteristic of RNP-type RBDs. Special structural features of RBD2 include a beta-bulge in beta2 and a shallow twist of the beta-sheet. The smaller 1H-15N nuclear Overhauser enhancement values for the residues of loop 3 between beta2 and beta3 suggest that this loop is flexible in the time-scale of nano- to picosecond order. The smaller 15N T2 values for the residues around the border between alpha2 and the following loop (loop 5) suggest this region undergoes conformational exchange in the milli- to microsecond time-scale. Chemical shift perturbation analysis indicated that RBD2 binds to an RNA oligomer obtained by in vitro selection under the conditions for NMR measurements, and thus the nature of the weak RNA-binding of RBD2 was successfully characterized by NMR, which is otherwise difficult to assess. Mainly the residues of the surface composed of the four-stranded beta-sheet, loops and C-terminal region are involved in the interaction. The appearance of side-chain NH proton resonances of arginine residues of loop 3 and imino proton resonances of RNA bases upon complex formation suggests the formation of intermolecular hydrogen bonds. The structural arrangement of the rings of the conserved aromatic residues of beta2 and beta3 is suitable for stacking interaction with RNA bases, known to be one of the major protein-RNA interactions, but a survey of the perturbation data suggested that the stacking interaction is not ideally achieved in the complex, which may be related to the weaker RNA-binding of RBD2.     

Structure of the C-terminal RNA-binding domain of hnRNP D0 (AUF1), its interactions with RNA and DNA, and change in backbone dynamics upon complex formation with DNA

Heterogeneous nuclear ribonucleoprotein (hnRNP) D0 has two ribonucleoprotein (RNP) -type RNA-binding domains (RBDs), each of which can specifically bind to the UUAG-sequence. hnRNP D0 also binds specifically to single-stranded d(TTAGGG)(n), the human telomeric DNA repeat. We have already reported the structure and interactions with RNA of the N-terminal RBD (RBD1). Here, the structure of the C-terminal RBD (RBD2) determined by NMR is presented. It folds into a compact alpha beta structure comprising an antiparallel beta-sheet packed against two alpha-helices, which is characteristic of RNP-type RBDs. In addition to the four beta-strands commonly found in RNP-type RBDs, an extra beta-strand, termed beta 4(-), was found just before the fourth beta-strand, yielding a five-stranded beta-sheet. Candidate residues of RBD2 involved in the interactions with RNA were identified by chemical shift perturbation analysis. Perturbation was detected on the beta-sheet side, not on the opposite alpha-helix side, as observed for RBD1. It is notable that the beta 4(-) to beta 4 region of RBD2 is involved in the interactions in contrast to the case of RBD1. The chemical shift perturbation analysis also showed that RBD2 interacts with DNA in essentially the same way as with RNA. Changes in the backbone dynamics upon complex formation with DNA were examined by means of model free analysis of relaxation data. In free RBD2, the beta 4(-) to beta 4 region exhibits slow conformational exchange on the milli- to microsecond time scale. The exchange is quenched upon complex formation. The flexibility of free RBD2 may be utilized in the recognition process by allowing different conformational states to be accessed and facilitating induced fit. Additionally, faster flexibility on the nano- to picosecond time scale was observed for loop 3 located between beta 2 and beta 3 in free RBD2, which is retained by the complex as well.     

Structure of Musashi1 in a complex with target RNA: the role of aromatic stacking interactions

Mammalian Musashi1 (Msi1) is an RNA-binding protein that regulates the translation of target mRNAs, and participates in the maintenance of cell 'stemness' and tumorigenesis. Msi1 reportedly binds to the 3'-untranslated region of mRNA of Numb, which encodes Notch inhibitor, and impedes initiation of its translation by competing with eIF4G for PABP binding, resulting in triggering of Notch signaling. Here, the mechanism by which Msi1 recognizes the target RNA sequence using its Ribonucleoprotein (RNP)-type RNA-binding domains (RBDs), RBD1 and RBD2 has been revealed on identification of the minimal binding RNA for each RBD and determination of the three-dimensional structure of the RBD1:RNA complex. Unique interactions were found for the recognition of the target sequence by Msi1 RBD1: adenine is sandwiched by two phenylalanines and guanine is stacked on the tryptophan in the loop between β1 and α1. The minimal recognition sequences that we have defined for Msi1 RBD1 and RBD2 have actually been found in many Msi1 target mRNAs reported to date. The present study provides molecular clues for understanding the biology involving Musashi family proteins.     

The roles of S1 RNA-binding domains in Rrp5's interactions with pre-rRNA

RNA-binding proteins mediate the function of all RNAs. Since few distinct RNA-binding domains (RBDs) exist, with most RBDs contacting only a few nucleotides, RNA-binding proteins often combine multiple RNA-binding motifs to achieve a higher affinity and selectivity for their targets. Rrp5, a ribosome assembly factor essential for both 40S and 60S ribosome maturation, is an extreme example as it contains 12 tandem S1 RNA-binding domains. In this study, we use a combination of RNA binding and DMS probing experiments to probe interactions of Rrp5 with pre-rRNA mimics. Our data localize Rrp5's binding site to three distinct regions within internal transcribed spacer 1 (ITS1), the sequence between 18S and 5.8S rRNAs. One of these regions is directly adjacent to a recently uncovered helical structure, which prevents premature cleavage at the 3'-end of 18S rRNA. This finding, together with previous results, suggests a role for Rrp5 in regulating the above-mentioned helical element. Furthermore, we have produced two truncated forms of the protein, Rrp5N and Rrp5C, which together encompass the entire protein and fully restore growth. Quantitative analysis of the RNA affinity of these Rrp5 fragments indicates that the first nine S1 motifs contribute much of Rrp5's RNA affinity, while the last three domains alone provide its specificity for the pre-rRNA. This surprising division of labor is unique, as it suggests that S1 domains can bind RNA both specifically as well as nonspecifically with high affinity; this has important implications for the molecular details of the Rrp5?pre-rRNA complex.     

NMR studies on functional structures of the AU-rich element-binding domains of Hu antigen C

Hu antigen C (HuC) has three RNA-binding domains (RBDs). The N-terminal two, RBD1 and RBD2, are linked in tandem and bind to the AU-rich elements (AREs) in the 3′-untranslated region of particular mRNAs. The solution structures of HuC RBD1 and RBD2 were determined by NMR methods. The HuC RBD1 and RBD2 structures are quite similar to those of Sxl RBD1 and RBD2, respectively. The individual RBDs of HuC, RBD1 and RBD2 in isolation can interact rather weakly with the minimal ARE motif, AUUUA, while the didomain fragment, RBD1–RBD2, of HuC binds more tightly to a longer ARE RNA, UAUUUAUUUU. Chemical shift perturbations by the longer RNA on HuC RBD1–RBD2 were mapped on and around the two β-sheets and on the C-terminal region of RBD1. The HuC RBD1–RBD2 residues that exhibited significant chemical shift perturbations coincide with those conserved in Sxl RBD1–RBD2. These data indicate that the RNA-binding characteristics of the HuC and Sxl didomain fragments are similar, even though the target RNAs and the biological functions of the proteins are different.     

Interaction of nucleolin with an evolutionarily conserved pre-ribosomal RNA sequence is required for the assembly of the primary processing complex

The first processing event of the precursor ribosomal RNA (pre-rRNA) takes place within the 5' external transcribed spacer. This primary processing requires conserved cis-acting RNA sequence downstream from the cleavage site and several nucleic acids (small nucleolar RNAs) and proteins trans-acting factors including nucleolin, a major nucleolar protein. The specific interaction of nucleolin with the pre-rRNA is required for processing in vitro. Xenopus laevis and hamster nucleolin interact with the same pre-rRNA site and stimulate the processing activity of a mouse cell extract. A highly conserved 11-nucleotide sequence located 5-6 nucleotides after the processing site is required for the interaction of nucleolin and processing. In vitro selection experiments with nucleolin have identified an RNA sequence that contains the UCGA motif present in the 11-nucleotide conserved sequence. The interaction of nucleolin with pre-rRNA is required for the formation of an active processing complex. Our findings demonstrate that nucleolin is a key factor for the assembly and maturation of pre-ribosomal ribonucleoparticles.     

N-Glycosylation influences the structure and self-association abilities of recombinant nucleolin

Nucleolin is a major nucleolar protein involved in fundamental processes of ribosome biogenesis, regulation of cell proliferation and growth. Nucleolin is known to shuttle between nucleus, cytoplasm and cell surface. We have previously found that nucleolin undergoes complex N- and O-glycosylations in extra-nuclear isoforms. We found that surface nucleolin is exclusively glycosylated and that N-glycosylation is required for its expression on the cells. Interestingly, the two N-glycans are located in the RNA-binding domains (RBDs) which participate in the self-association properties of nucleolin. We hypothesized that the occupancy of RBDs by N-glycans plays a role in these self-association properties. Here, owing to the inability to quantitatively produce full-size nucleolin, we expressed four N-glycosylation nucleolin variants lacking the N-terminal acidic domain in a baculovirus/insect cell system. As assessed by heptafluorobutyrate derivatization and mass spectrometry, this strategy allowed the production of proteins bearing or not paucimannosidic-type glycans on either one or two of the potential N-glycosylation sites. Their structure was investigated by circular dichroism and fluorimetry, and their ability to self-interact was analyzed by electrophoresis and surface plasmon resonance. Our results demonstrate that all nucleolin-derived variants are able to self-interact and that N-glycosylation on both RBD1 and RBD3, or RBD3 alone, but not RBD1 alone, modifies the structure of the N-terminally truncated nucleolin and enhances its self-association properties. In contrast, N-glycosylation does not modify interaction with lactoferrin, a ligand of cell surface nucleolin. Our results suggest that the occupancy of the N-glycosylation sites may contribute to expression and functions of surface nucleolin.     

Origin of higher affinity to RNA of the N-terminal RNA-binding domain than that of the C-terminal one of a mouse neural protein,musashi1, as revealed by comparison of their structures,modes of interaction,surface electrostatic potentials,and backbone dynamics

Musashi1 is an RNA-binding protein abundantly expressed in the developing mouse central nervous system. Its restricted expression in neural precursor cells suggests that it is involved in maintenance of the character of progenitor cells. Musashi1 contains two ribonucleoprotein-type RNA-binding domains (RBDs), RBD1 and RBD2, the affinity to RNA of RBD1 being much higher than that of RBD2. We previously reported the structure and mode of interaction with RNA of RBD2. Here, we have determined the structure and mode of interaction with RNA of RBD1. We have also analyzed the surface electrostatic potential and backbone dynamics of both RBDs. The two RBDs exhibit the same ribo-nucleoprotein-type fold and commonly make contact with RNA on the beta-sheet side. On the other hand, there is a remarkable difference in surface electrostatic potential, the beta-sheet of RBD1 being positively charged, which is favorable for binding negatively charged RNA, but that of RBD2 being almost neutral. There is also a difference in backbone dynamics, the central portion of the beta-sheet of RBD1 being flexible, but that of RBD2 not being flexible. The flexibility of RBD1 may be utilized in the recognition process to facilitate an induced fit. Thus, comparative studies have revealed the origin of the higher affinity of RBD1 than that of RBD2 and indicated that the affinity of an RBD to RNA is not governed by its fold alone but is also determined by its surface electrostatic potential and/or backbone dynamics. The biological role of RBD2 with lower affinity is also discussed.     

Covariance analysis of RNA recognition motifs identifies functionally linked amino acids.

The RNA recognition motif (RRM) is one of the most common eukaryotic protein motifs. RRM sequences form a conserved globular structure known as the RNA-binding domain (RBD) or the ribonucleoprotein domain. Many proteins that contain RRM sequences bind RNA in a sequence-specific manner. To investigate the basis for the RNA-binding specificity of RRMs, we subjected 330 aligned RRM sequences to covariance analysis. The analysis revealed a single network of covariant amino acid pairs comprising the buried core of the RBD and a surface patch. Structural studies have implicated a subset of these residues in RNA binding. The covariance linkages identify a larger set of amino acid residues, including some not directly in contact with bound RNA, that may influence RNA-binding specificity.