Contribution of the C-terminal tail of U1A RBD1 to RNA recognition and protein stability.

The N-terminal RNA binding domain (RBD1) of the human U1A protein interacts specifically with a short RNA hairpin containing the U1 snRNA stem/loop II sequence. Previous RNA binding studies have suggested that the C-terminal tail of RBD1 contributes to RNA recognition in addition to interactions on the beta-sheet surface of the protein. To evaluate the contributions of these C-terminal residues in RBD1 to RNA binding affinity and specificity, as well as to study the thermodynamic stability of RBDs, a number of RBD1 mutants with truncated tails, with single amino acid substitutions, and with both a truncation and an amino acid substitution, have been constructed. The thermodynamic stabilities of these mutants have been measured and compared by GdnHCI unfolding experiments. The RNA binding affinity and specificity of these mutant proteins have been assessed by measuring the binding of each protein to the wild-type RNA hairpin and to selected RNA mutants with nucleotide substitutions in the RNA loop. The results demonstrate first that, although the C-terminal tail of RBD1 makes significant contributions to RNA binding affinity, it is not required for RNA binding, and second, its contributions to binding specificity are mediated only through selected nucleotides in the RNA loop, for in the absence of the tail, the protein continues to use other nucleotides to discriminate among RNAs. In these truncated proteins, the secondary structure intrinsic to the C-terminal tail is absent, yet their affinity and discrimination for RNAs are not lost. Thus, a structured tail is not required for RNA recognition.     

Correlated motions in the U1 snRNA stem/loop 2:U1A RBD1 complex

The complex formed by U1A RBD1 and the U1 snRNA stem/loop II is noted for its high affinity and exquisite specificity. Here, that complex is investigated by 5 ns molecular dynamics simulations and analyzed by reorientational eigenmode dynamics to determine the dynamic properties of the RNA:protein interface that could contribute to the binding mechanism. The analysis shows that there is extensive correlation between motions of the RNA and protein, involving 7 of the 10 RNA loop nucleotides, the protein beta-sheet surface, two of its loops, and its C-terminal tripeptide sequence. Order parameters of these regions of the complex are uniformly high, indicating restricted motion. However, several regions of both RNA and protein retain local flexibility, notably three nucleotides of the RNA loop and one loop of RBD1 that does not contact RNA. The highly correlated motions involving both molecules reflect the intricate network of interactions that characterize this complex and could account in part for the thermodynamic coupling observed for complex formation.     

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Mutational definition of RNA-binding and protein-protein interaction domains of heterogeneous nuclear RNP C1

The heterogeneous nuclear ribonucleoprotein (hn- RNP) C proteins, among the most abundant pre-mRNA-binding proteins in the eukaryotic nucleus, have a single RNP motif RNA-binding domain. The RNA-binding domain (RBD) is comprised of approximately 80-100 amino acids, and its structure has been determined. However, relatively little is known about the role of specific amino acids of the RBD in the binding to RNA. We have devised a phage display-based screening method for the rapid identification of amino acids in hnRNP C1 that are essential for its binding to RNA. The identified mutants were further tested for binding to poly(U)-Sepharose, a substrate to which wild type hnRNP C1 binds with high affinity. We found both previously predicted, highly conserved residues as well as additional residues in the RBD to be essential for C1 RNA binding. We also identified three mutations in the leucine-rich C1-C1 interaction domain near the carboxyl terminus of the protein that both abolished C1 oligomerization and reduced RNA binding. These results demonstrate that although the RBD is the primary determinant of C1 RNA binding, residues in the C1-C1 interaction domain also influence the RNA binding activity of the protein. The experimental approach we described should be generally applicable for the screening and identification of amino acids that play a role in the binding of proteins to nucleic acid substrates.     

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.     

Solution structure of the complex formed by the two N-terminal RNA-binding domains of nucleolin and a pre-rRNA target

Nucleolin is a 70 kDa multidomain protein involved in several steps of eukaryotic ribosome biogenesis. In vitro selection in combination with mutagenesis and structural analysis identified binding sites in pre-rRNA with the consensus (U/G)CCCG(A/G) in the context of a hairpin structure, the nucleolin recognition element (NRE). The central region of the protein contains four tandem RNA-binding domains (RBDs), of which the first two are responsible for the RNA-binding specificity and affinity for NREs. Here, we present the solution structure of the 28 kDa complex formed by the two N-terminal RNA-binding domains of nucleolin (RBD12) and a natural pre-rRNA target, b2NRE. The structure demonstrates that the sequence-specific recognition of the pre-rRNA NRE is achieved by intermolecular hydrogen bonds and stacking interactions involving mainly the beta-sheet surfaces of the two RBDs and the linker residues. A comparison with our previously determined NMR structure of RBD12 in complex with an in vitro selected RNA target, sNRE, shows that although the sequence-specific recognition of the loop consensus nucleotides is the same in the two complexes, they differ in several aspects. While the protein makes numerous specific contacts to the non-consensus nucleotides in the loop E motif (S-turn) in the upper part of the sNRE stem, nucleolin RBD12 contacts only consensus nucleotides in b2NRE. The absence of these upper stem contacts from the RBD12/b2NRE complex results in a much less stable complex, as demonstrated by kinetic analyses. The role of the loop E motif in high-affinity binding is supported by gel-shift analyses with a series of sNRE mutants. The less stable interaction of RBD12 with the natural RNA target is consistent with the proposed role of nucleolin as a chaperone that interacts transiently with pre-rRNA to prevent misfolding.     

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.     

Contribution of the tyrosines to the structure and function of the human U1A N-terminal RNA binding domain.

RNA binding domains (RBDs) are members of a large family of proteins that share minimal sequence conservation but adopt an alpha beta sandwich global fold. Defining the contributions of specific amino acids to RBD structure and RNA binding is critical to understanding the functions of these proteins. In these experiments with the human U1A N-terminal RNA binding domain (RBD1), the contributions from each of its four tyrosines to protein structure, stability, and RNA binding were measured. Each tyrosine was substituted with phenylalanine and one other selected residue, and the resulting proteins were characterized by chemical denaturation to measure their unfolding free energy, by binding free energies to the wild-type RNA hairpin, and by 19F NMR to probe for structural changes. Features of the protein identified in these experiments include a possible tyrosine/lysine contact in an alpha-helix, which may be an example of an energetically favorable aromatic/amino side chain interaction. One long loop of the protein, which shows unusual 15N backbone and tyrosine side-chain dynamics, is implicated in protein:protein association. The diverse interactions of the four tyrosine residues in the organization of RBD1 illustrate how each member of this family of proteins will have unique molecular details that contribute to function.     

Defining the orientation of the human U1A RBD1 on its UTR by tethered-EDTA(Fe) cleavage.

The N-terminal RNA binding domain of the human U1A protein (RBD1) specifically binds an RNA hairpin of U1 snRNA as well as two internal loops in the 3' UTR of its own mRNA. Here, a single cysteine has been introduced into Loop 1 of RBD1, which is subsequently used to attach (EDTA-2-aminoethyl) 2-pyridyl disulfide-Fe3+ (EPD-Fe). This EDTA-Fe derivative is used to generate hydroxyl radicals to cleave the proximal RNA sugar-phosphate backbone in the RNA-RBD complexes. RBD1(K20C)-EPD-Fe cleaves the 5' strand of the RNA hairpin stem, centered four base pairs away from the base of the loop, and cleaves the UTR in two places, again centered on the 5' side of the fourth base pair from each internal loop. These data, extrapolated to the position of Lys 20 in RBD1, orient the two proteins bound to the UTR, and provide direct biochemical evidence for the proposed model of the RBD1:UTR complex.     

Crosslinking of an iodo-uridine-RNA hairpin to a single site on the human U1A N-terminal RNA binding domain.

The N-terminal RNA binding domain (RBD) of the human U1A snRNP protein binds tightly and specifically to an RNA hairpin that contains a 10-nucleotide loop. The protein is one of a class of RNA binding proteins that adopts a beta alpha beta beta alpha beta global fold, which in turn forms a four-stranded antiparallel beta-sheet. This sheet forms the primary binding surface for the RNA, as shown by the crosslinking results described here, and in more detail by a recently described co-crystal of this RBD with an RNA hairpin (Oubridge C, et al., 1994, Nature 372:432-438). The RNA hairpin sequence used in the crosslinking experiments, containing 5-iodo-uridine, is a variant of the normal U1 snRNA sequence which is able to form a crosslink with the protein, in contrast to the wild-type sequence, which does not. This single uridine substitution in the 10-nucleotide loop is the site of cross-linking to one tyrosine (Tyr 13) in the beta 1 strand of the U1A N-terminal RBD. This same uridine is also crosslinked to a mutant Tyr 13 Phe RBD, at this Phe 13 substitution.     

Interaction of N-terminal domain of U1A protein with an RNA stem/loop.

The U1A protein is a sequence-specific RNA binding protein found in the U1 snRNP particle where it binds to stem/loop II of U1 snRNA. U1A contains two 'RNP' or 'RRM' (RNA Recognition Motif) domains, which are common to many RNA-binding proteins. The N-terminal RRM has been shown to bind specifically to the U1 RNA stem/loop, while the RNA target of the C-terminal domain is unknown. Here, we describe experiments using a 102 amino acid N-terminal RRM of U1A (102A) and a 25-nucleotide RNA stem/loop to measure the binding constants and thermodynamic parameters of this RNA:protein complex. Using nitrocellulose filter binding, we measure a dissociation constant KD = 2 x 10(-11) M in 250 mM NaCl, 2 mM MgC2, and 10 mM sodium cacodylate, pH 6 at room temperature, and a half-life for the complex of 5 minutes. The free energy of association (delta G degrees) of this complex is about -14 kcal/mol in these conditions. Determination of the salt dependence of the binding suggests that at least 8 ion-pairs are formed upon complex formation. A mutation in the RNA loop sequence reduces the affinity 10 x, or about 10% of the total free energy.     

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'.     

Climbing the vertebrate branch of U1A/U2B″ protein evolution

In the vertebrate lineage of the U1A/U2B″/SNF protein family, the U1A and U2B″ proteins bind to RNA stem–loops in the U1 or U2 snRNPs, respectively. However, their specialization is fairly recent, as they evolved from a single ancestral protein. The progress of their specialization (subfunctionalization) can be monitored by the amino acid sequence changes that give rise to their modern RNA-binding specificity. Using ancestral sequence reconstruction to predict the intermediates on the evolutionary branch, a probable path of sequential changes is defined for U1A and U2B″. The RNA-binding affinity for U1A/U2B″ protein ancestors was measured using modern U1 and U2 snRNA stem–loops and RNA stem–loop variants to understand how the proteins’ RNA specificities evolved.     

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.     

Characterization of yeast U1 snRNP A protein: identification of the N-terminal RNA binding domain (RBD) binding site and evidence that the C-terminal RBD functions in splicing.

The yeast U1A protein is a U1 snRNP-specific protein. Like its human counterpart (hU1A), it has two conserved RNA binding domains (RBDs). The N-terminal RBD is quite different from the human protein, and a binding site on yeast U1 snRNA is not readily apparent. The C-terminal RBD is of unknown function. Using in vivo dimethyl sulfate (DMS) protection of mutant strains, we defined a region in yeast U1 snRNA as the likely U1A N-terminal RBD binding site. This was confirmed by direct in vitro binding assays. The site is very different from its vertebrate counterpart, but its location within yeast U1 snRNA suggests a conserved structural relationship to other U1 snRNP components. Genetic studies and sensitive in vivo splicing measurements indicate that the yeast U1A C-terminal RBD also functions in pre-mRNA splicing. We propose that the N-terminal RBD serves to tether the splicing-relevant C-terminal RBD to the snRNP.     

Polyadenylation proteins CstF-64 and tauCstF-64 exhibit differential binding affinities for RNA polymers

CstF-64 (cleavage stimulation factor-64), a major regulatory protein of polyadenylation, is absent during male meiosis. Therefore a paralogous variant, tauCstF-64 is expressed in male germ cells to maintain normal spermatogenesis. Based on sequence differences between tauCstF-64 and CstF-64, and on the high incidence of alternative polyadenylation in testes, we hypothesized that the RBDs (RNA-binding domains) of tauCstF-64 and CstF-64 have different affinities for RNA elements. We quantified K(d) values of CstF-64 and tauCstF-64 RBDs for various ribopolymers using an RNA cross-linking assay. The two RBDs had similar affinities for poly(G)18, poly(A)18 or poly(C)18, with affinity for poly(C)18 being the lowest. However, CstF-64 had a higher affinity for poly(U)18 than tauCstF-64, whereas it had a lower affinity for poly(GU)9. Changing Pro-41 to a serine residue in the CstF-64 RBD did not affect its affinity for poly(U)18, but changes in amino acids downstream of the C-terminal alpha-helical region decreased affinity towards poly(U)18. Thus we show that the two CstF-64 paralogues differ in their affinities for specific RNA sequences, and that the region C-terminal to the RBD is mportant in RNA sequence recognition. This supports the hypothesis that tauCstF-64 promotes germ-cell-specific patterns of polyadenylation by binding to different downstream sequence elements.     

Recognition of two distinct elements in the RNA substrate by the RNA-binding domain of the T. thermophilus DEAD box helicase Hera

DEAD box helicases catalyze the ATP-dependent destabilization of RNA duplexes. Whereas duplex separation is mediated by the helicase core shared by all members of the family, flanking domains often contribute to binding of the RNA substrate. The Thermus thermophilus DEAD-box helicase Hera (for “heat-resistant RNA-binding ATPase”) contains a C-terminal RNA-binding domain (RBD). We have analyzed RNA binding to the Hera RBD by a combination of mutational analyses, nuclear magnetic resonance and X-ray crystallography, and identify residues on helix α1 and the C-terminus as the main determinants for high-affinity RNA binding. A crystal structure of the RBD in complex with a single-stranded RNA resolves the RNA–protein interactions in the RBD core region around helix α1. Differences in RNA binding to the Hera RBD and to the structurally similar RBD of the Bacillus subtilis DEAD box helicase YxiN illustrate the versatility of RNA recognition motifs as RNA-binding platforms. Comparison of chemical shift perturbation patterns elicited by different RNAs, and the effect of sequence changes in the RNA on binding and unwinding show that the RBD binds a single-stranded RNA region at the core and simultaneously contacts double-stranded RNA through its C-terminal tail. The helicase core then unwinds an adjacent RNA duplex. Overall, the mode of RNA binding by Hera is consistent with a possible function as a general RNA chaperone.     

Target discrimination by RNA-binding proteins: role of the ancillary protein U2A' and a critical leucine residue in differentiating the RNA-binding specificity of spliceosomal proteins U1A and U2B".

The spliceosomal proteins U1A and U2B" each use a homologous RRM domain to bind specifically to their respective snRNA targets, U1hpll and U2hpIV, two stem-loops that are similar yet distinct in sequence. Previous studies have shown that binding of U2B" to U2hpIV is facilitated by the ancillary protein U2A', whereas specific binding of U1A to U1hpll requires no cofactor. Here we report that U2A' enables U2B" to distinguish the loop sequence of U2hpIV from that of U1hpll but plays no role in stem sequence discrimination. Although U2A' can also promote heterospecific binding of U1A to U2hpIV, a much higher concentration of the ancillary protein is required due to the approximately 500-fold greater affinity of U2A' for U2B". Additional experiments have identified a single leucine residue in U1A(Leu-44) that is critical for the intrinsic specificity of this protein for the loop sequence of U1 hpll in preference to that of U2hpIV. Our data suggest that most of the difference in RNA-binding specificity between U1A and U2B" can be accounted for by this leucine residue and by the contribution of the ancillary protein U2A' to the specificity of U2B".     

Anti-(Raf-1) RNA aptamers that inhibit Ras-induced Raf-1 activation.

RNA aptamers with affinity for the Ras-binding domain (RBD) of Raf-1 were isolated from a degenerate pool by in vitro selection. These aptamers efficiently inhibited the Ras interaction with the Raf-1 RBD, and also inhibited Ras-induced Raf-1 activation in a cell-free system. The RNA aptamer with the most potent inhibitory effect specifically inhibited the Ras-Raf-1 interaction and had no affinity for the RBD of the RGL protein, a homolog of the Ral GDP dissociation stimulator. Although the aptamer was capable of binding to the B-Raf RBD, the RNA did not inhibit the interaction between Ras and the B-Raf RBD. Enzymatic and chemical probing experiments indicated that the aptamer was folded into a pseudoknot structure, and some loop regions of the pseudoknot were located at the binding interface for the Raf-1 RBD.