RNA binding mediates the local cooperativity between the beta-sheet and the C-terminal tail of the human U1A RBD1 protein

A study on allosteric communication in U1A-snRNA binding interactions: network analysis combined with molecular dynamics data

The allosteric regulation during the binding interactions between small nuclear RNAs (snRNAs) and the associated protein factors is critical to the function of spliceosomes in alternative RNA splicing. Although network models combined with molecular dynamics simulations have shown to be powerful tools for the analysis of protein allostery, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems. In this work, we use a residual network model combined with a coarse-grained Gaussian network model (GNM) to investigate the binding interactions between the snRNA and the human U1A protein which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle, and to identify the residues that play an important role in the allosteric communication in U1A during this process. We also utilize the Girvan-Newman method to detect the structural organization in U1A-snRNA recognition and interactions. Our results reveal that: (Ι) not only the residues at the binding sites that are traditionally considered to play a major role in U1A-snRNA association, but those residues that are far away from the RNA binding interface participate in the U1A's allosteric signal transmission induced by the RNA binding; (Π) the structure of U1A protein is well organized with different communities acting different roles for its RNA binding and allosteric regulation. The study demonstrates that the combination of the residual network and elastic network models is an effective and efficient method which can be readily extended to the investigation of the allosteric communication for other macromolecular interaction systems.

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Comparative analyses of the thermodynamic RNA binding signatures of different types of RNA recognition motifs

RNA recognition motifs (RRMs) are structurally versatile domains important in regulation of alternative splicing. Structural mechanisms of sequence-specific recognition of single-stranded RNAs (ssRNAs) by RRMs are well understood. The thermodynamic strategies are however unclear. Therefore, we utilized microcalorimetry and semi-empirical analyses to comparatively analyze the cognate ssRNA binding thermodynamics of four different RRM domains, each with a different RNA binding mode. The different binding modes are: canonical binding to the β-sheet surface; canonical binding with involvement of N- and C-termini; binding to conserved loops; and binding to an α-helix. Our results identify enthalpy as the sole and general force driving association at physiological temperatures. Also, networks of weak interactions are a general feature regulating stability of the different RRM–ssRNA complexes. In agreement, non-polyelectrolyte effects contributed between ∼75 and 90% of the overall free energy of binding in the considered complexes. The various RNA binding modes also displayed enormous heat capacity differences, that upon dissection revealed large differential changes in hydration, conformations and dynamics upon binding RNA. Altogether, different modes employed by RRMs to bind cognate ssRNAs utilize various thermodynamics strategies during the association process.

Co-evolution of SNF spliceosomal proteins with their RNA targets in trans-splicing nematodes

Although the mechanism of pre-mRNA splicing has been well characterized, the evolution of spliceosomal proteins is poorly understood. The U1A/U2B″/SNF family (hereafter referred to as the SNF family) of RNA binding spliceosomal proteins participates in both the U1 and U2 small interacting nuclear ribonucleoproteins (snRNPs). The highly constrained nature of this system has inhibited an analysis of co-evolutionary trends between the proteins and their RNA binding targets. Here we report accelerated sequence evolution in the SNF protein family in Phylum Nematoda, which has allowed an analysis of protein:RNA co-evolution. In a comparison of SNF genes from ecdysozoan species, we found a correlation between trans-splicing species (nematodes) and increased phylogenetic branch lengths of the SNF protein family, with respect to their sister clade Arthropoda. In particular, we found that nematodes (~70-80 % of pre-mRNAs are trans-spliced) have experienced higher rates of SNF sequence evolution than arthropods (predominantly cis-spliced) at both the nucleotide and amino acid levels. Interestingly, this increased evolutionary rate correlates with the reliance on trans-splicing by nematodes, which would alter the role of the SNF family of spliceosomal proteins. We mapped amino acid substitutions to functionally important regions of the SNF protein, specifically to sites that are predicted to disrupt protein:RNA and protein:protein interactions. Finally, we investigated SNF's RNA targets: the U1 and U2 snRNAs. Both are more divergent in nematodes than arthropods, suggesting the RNAs have co-evolved with SNF in order to maintain the necessarily high affinity interaction that has been characterized in other species.

A compare-and-contrast NMR dynamics study of two related RRMs: U1A and SNF

The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/β sandwich topology, and these RRMs use their β-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the β-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 β-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.

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.

The role of the C-terminal helix of U1A protein in the interaction with U1hpII RNA

Previous kinetic investigations of the N-terminal RNA Recognition Motif (RRM) domain of spliceosomal A protein of the U1 small nuclear ribonucleoprotein particle (U1A) interacting with its RNA target U1 hairpin II (U1hpII) provided experimental evidence for a ‘lure and lock’ model of binding. The final step of locking has been proposed to involve conformational changes in an α-helix immediately C-terminal to the RRM domain (helix C), which occludes the RNA binding surface in the unbound protein. Helix C must shift its position to accommodate RNA binding in the RNA–protein complex. This results in a new hydrophobic core, an intraprotein hydrogen bond and a quadruple stacking interaction between U1A and U1hpII. Here, we used a surface plasmon resonance-based biosensor to gain mechanistic insight into the role of helix C in mediating the interaction with U1hpII. Truncation, removal or disruption of the helix exposes the RNA-binding surface, resulting in an increase in the association rate, while simultaneously reducing the ability of the complex to lock, reflected in a loss of complex stability. Disruption of the quadruple stacking interaction has minor kinetic effects when compared with removal of the intraprotein hydrogen bonds. These data provide new insights into the mechanism whereby sequences C-terminal to an RRM can influence RNA binding.

Human U2B″ protein binding to snRNA stemloops

The human U2B″ protein is one of the unique proteins that comprise the U2 snRNP, but it is also a representative of the U1A/U2B″ protein family. In the U2 snRNP, it is bound to Stem-Loop IV (SLIV) of the U2 snRNA. We find that in vitro it binds not only to human SLIV, but also to Stem-Loop II (SLII) from human U1 snRNA and to Drosophila U2 snRNA SLIV. The thermodynamics of these binding interactions show a striking similarity, leading to the conclusion that U2B″ has a relaxed specificity for its RNA targets. The binding properties of U2B″ are distinct from those of human U1A and of Drosophila SNF, despite its high homology to those proteins, and so provide important new information on how this protein family has modulated its target preferences.

Multistep kinetics of the U1A-SL2 RNA complex dissociation

The U1A-SL2 RNA complex is a model system for studying interactions between RNA and the RNA recognition motif (RRM), which is one of the most common RNA binding domains. We report here kinetic studies of dissociation of the U1A-SL2 RNA complex, using laser temperature jump and stopped-flow fluorescence methods with U1A proteins labeled with the intrinsic chromophore tryptophan. An analysis of the kinetic data suggests three phases of dissociation with time scales of ∼100 μs, ∼50 ms, and ∼2 s. We propose that the first step of dissociation is a fast rearrangement of the complex to form a loosely bound complex. The intermediate step is assigned to be the dissociation of the U1A-SL2 RNA complex, and the final step is assigned to a reorganization of the U1A protein structure into the conformation of the free protein. These assignments are consistent with previous proposals based on thermodynamic, NMR, and surface plasmon resonance experiments and molecular dynamics simulations. Together, these results begin to build a comprehensive model of the complex dynamic processes involved in the formation and dissociation of an RRM-RNA complex.

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.

Bruno protein contains an expanded RNA recognition motif

The RNA recognition motif (or RRM) is a ubiquitous RNA-binding module present in approximately 2% of the proteins encoded in the human genome. This work characterizes an expanded RRM, which is present in the Drosophila Bruno protein, and targets regulatory elements in the oskar mRNA through which Bruno controls translation. In this Bruno RRM, the deletion of 40 amino acids prior to the N-terminus of the canonical RRM resulted in a significantly decreased affinity of the protein for its RNA target. NMR spectroscopy showed that the expanded Bruno RRM contains the familiar RRM fold of four antiparallel beta-strands and two alpha-helices, preceded by a 10-residue loop that contacts helix alpha(1) and strand beta(2); additional amino acids at the N-terminus of the domain are relatively flexible in solution. NMR results also showed that a truncated form of the Bruno RRM, lacking the flexible N-terminal amino acids, forms a stable and complete canonical RRM, so that the loss of RNA binding activity cannot be attributed to disruption of the RRM fold. This expanded Bruno RRM provides a new example of the features that are important for RNA recognition by an RRM-containing protein.

Functional stabilization of an RNA recognition motif by a noncanonical N-terminal expansion

RNA recognition motifs (RRMs) constitute versatile macromolecular interaction platforms. They are found in many components of spliceosomes, in which they mediate RNA and protein interactions by diverse molecular strategies. The human U11/U12-65K protein of the minor spliceosome employs a C-terminal RRM to bind hairpin III of the U12 small nuclear RNA (snRNA). This interaction comprises one side of a molecular bridge between the U11 and U12 small nuclear ribonucleoprotein particles (snRNPs) and is reminiscent of the binding of the N-terminal RRMs in the major spliceosomal U1A and U2B'' proteins to hairpins in their cognate snRNAs. Here we show by mutagenesis and electrophoretic mobility shift assays that the beta-sheet surface and a neighboring loop of 65K C-terminal RRM are involved in RNA binding, as previously seen in canonical RRMs like the N-terminal RRMs of the U1A and U2B'' proteins. However, unlike U1A and U2B'', some 30 residues N-terminal of the 65K C-terminal RRM core are additionally required for stable U12 snRNA binding. The crystal structure of the expanded 65K C-terminal RRM revealed that the N-terminal tail adopts an alpha-helical conformation and wraps around the protein toward the face opposite the RNA-binding platform. Point mutations in this part of the protein had only minor effects on RNA affinity. Removal of the N-terminal extension significantly decreased the thermal stability of the 65K C-terminal RRM. These results demonstrate that the 65K C-terminal RRM is augmented by an N-terminal element that confers stability to the domain, and thereby facilitates stable RNA binding.

Conformationally restricted nucleotides as a probe of structure-function relationships in RNA

We introduce the use of commercially available locked nucleic acids (LNAs) as a functional probe in RNA. LNA nucleotides contain a covalent linkage that restricts the pseudorotation phase of the ribose to C3'-endo (A-form). Introduction of an LNA at a single site thus allows the role of ribose structure and dynamics in RNA function to be assessed. We apply LNA probing at multiple sites to analyze self-cleavage in the lead-dependent ribozyme (leadzyme), thermodynamic stability in the UUCG tetraloop, and the kinetics of recognition of U1A protein by U1 snRNA hairpin II. In the leadzyme, locking a single guanosine residue into the C3'-endo pucker increases the catalytic rate by a factor of 20, despite the fact that X-ray crystallographic and NMR structures of the leadzyme ground state reported a C2'-endo conformation at this site. These results strongly suggest that a conformational change at this position is critical for catalytic function. Functional insights obtained in all three systems demonstrate the highly general applicability of LNA probing in analysis of the role of ribose orientation in RNA structure, dynamics, and function.

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.

Recognition of essential purines by the U1A protein

Background

The RNA recognition motif (RRM) is one of the largest families of RNA binding domains. The RRM is modulated so that individual proteins containing RRMs can specifically recognize RNA targets with diverse sequences and structures. Understanding the principles governing this specificity will be important for the rational modification and design of RRM-RNA complexes.

Results

In this paper we have investigated the origins of specificity of the N terminal RRM of the U1A protein for stem loop 2 (SL2) of U1 snRNA by substituting modified bases for essential purines in SL2 RNA. In one series of modified bases, hydrogen bond donors and acceptors were replaced by aliphatic groups to probe the importance of these functional groups to binding. In a second series of modified bases, hydrogen bond donors and acceptors were incorrectly placed on the purine bases to analyze the origins of discrimination between cognate and non-cognate RNA. The results of these experiments show that three different approaches are used by the U1A protein to gain specificity for purines. Specificity for the first base in the loop, A1, is based primarily on discrimination against RNA containing the incorrect base, specificity for the fourth base in the loop, G4, is based largely on recognition of the donors and acceptors of G4, while specificity for the sixth base in the loop, A6, results from a combination of direct recognition of the base and discrimination against incorrectly placed functional groups.

Conclusion

These investigations identify different roles that hydrogen bond donors and acceptors on bases in both cognate and non-cognate RNA play in the specific recognition of RNA by the U1A protein. Taken together with investigations of other RNA-RRM complexes, the results contribute to a general understanding of the origins of RNA-RRM specificity and highlight, in particular, the contribution of steric and electrostatic repulsion to binding specificity.

Do collective atomic fluctuations account for cooperative effects? Molecular dynamics studies of the U1A-RNA complex

A complete understanding of gene expression relies on a comprehensive understanding of the protein-RNA recognition process. However, the study of protein-RNA recognition is complicated by many factors that contribute to both binding affinity and specificity, including structure, energetics, dynamical motions, and cooperative interactions. Several recent studies have suggested that energetic coupling between residues contributes to formation of the complex between the U1A protein and stem loop 2 of U1 snRNA as a consequence of a cooperative network of interactions. We have performed molecular dynamics simulations on the U1A-RNA complex, including explicit water and counterions, and analyzed the results based on the calculated positional cross-correlations of atomic fluctuations. The results indicate that cross-correlations calculated on a per residue basis agree well with the observed inter-residue cooperativity and predict that the networks identified to date may also be coupled into an extensive hyper-network that reflects the intrinsic rigidity of the RNA recognition motif. In addition, we report a comparison of the MD calculated correlations with the results of a positional covariance analysis based on the sequences of 330 RNA recognition motifs, including U1A. The calculated inter-residue cross-correlations agree very well with the results of the sites exhibiting positional covariance. Collectively, these results strongly support the hypothesis that collective fluctuations contribute to cooperativity and the corresponding observed thermodynamic coupling. Predictions of additional sites in U1A that may be involved in cooperative networks are advanced.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

The role of RNA structure in the interaction of U1A protein with U1 hairpin II RNA

The N-terminal RNA Recognition Motif (RRM1) of the spliceosomal protein U1A interacting with its target U1 hairpin II (U1hpII) has been used as a paradigm for RRM-containing proteins interacting with their RNA targets. U1A binds to U1hpII via direct interactions with a 7-nucleotide (nt) consensus binding sequence at the 5' end of a 10-nt loop, and via hydrogen bonds with the closing C-G base pair at the top of the RNA stem. Using surface plasmon resonance (Biacore), we have examined the role of structural features of U1hpII in binding to U1A RRM1. Mutational analysis of the closing base pair suggests it plays a minor role in binding and mainly prevents "breathing" of the loop. Lengthening the stem and nontarget part of the loop suggests that the increased negative charge of the RNA might slightly aid association. However, this is offset by an increase in dissociation, which may be caused by attraction of the RRM to nontarget parts of the RNA. Studies of a single stranded target and RNAs with untethered loops indicate that structure is not very relevant for association but is important for complex stability. In particular, breaking the link between the stem and the 5' side of the loop greatly increases complex dissociation, presumably by hindering simultaneous contacts between the RRM and stem and loop nucleotides. While binding of U1A to a single stranded target is much weaker than to U1hpII, it occurs with nanomolar affinity, supporting recent evidence that binding of unstructured RNA by U1A has physiological significance.

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.

Kinetic analysis of the role of the tyrosine 13, phenylalanine 56 and glutamine 54 network in the U1A/U1 hairpin II interaction

The A protein of the U1 small nuclear ribonucleoprotein particle, interacting with its stem-loop RNA target (U1hpII), is frequently used as a paradigm for RNA binding by recognition motif domains (RRMs). U1A/U1hpII complex formation has been proposed to consist of at least two steps: electrostatically mediated alignment of both molecules followed by locking into place, based on the establishment of close-range interactions. The sequence of events between alignment and locking remains obscure. Here we examine the roles of three critical residues, Tyr13, Phe56 and Gln54, in complex formation and stability using Biacore. Our mutational and kinetic data suggest that Tyr13 plays a more important role than Phe56 in complex formation. Mutational analysis of Gln54, combined with molecular dynamics studies, points to Arg52 as another key residue in association. Based on our data and previous structural and modeling studies, we propose that electrostatic alignment of the molecules is followed by hydrogen bond formation between the RNA and Arg52, and the sequential establishment of interactions with loop bases (including Tyr13). A quadruple stack, sandwiching two bases between Phe56 and Asp92, would occur last and coincide with the rearrangement of a C-terminal helix that partially occludes the RRM surface in the free protein.

Protein and RNA Dynamics Play Key Roles in Determining the Specific Recognition of GU-rich Polyadenylation Regulatory Elements by Human Cstf-64 Protein

The N-terminal domain of the 64 kDa subunit of the cleavage stimulation factor (CstF-64) recognizes GU-rich elements within the 3'-untranslated region of eukaryotic mRNAs. This interaction is essential for mRNA 3' end processing and transcription termination, and its strength affects the efficiency of utilization of different polyadenylation sites. The structure of the RNA-binding N-terminal domain of CstF-64 showed how the N-terminal RNA recognition motif of CstF-64 recognizes GU-rich RNAs. However, it is still perplexing how this protein can bind selectively to RNAs that are rich in G and U residues regardless of their detailed sequence composition, yet discriminate effectively against non-GU-RNAs. We investigated by NMR the dynamics of the CstF-64 RNA-binding domain, both free and bound to two GU-rich RNA sequences that represent polyadenylation regulatory elements. While the free protein displays the motional properties typical of a well-folded protein domain and is uniformly rigid, the protein-RNA interface acquires significant mobility on the micro- to millisecond time-scale once GU-rich RNAs binds to it. These motional features, we propose, are intrinsic to the functional requirement to bind all GU-rich sequences and yet to discriminate against non-GU-rich RNAs. This behavior may be a general mechanism by which some RNA-binding proteins are able to bind to classes of sequences, as opposed to a well-defined sequence or consensus.

A binding mechanism in protein-nucleotide interactions: implication for U1A RNA binding

We present a close electronic view of the protein-base interface for the N-terminal domain of the human protein U1A. Combining accurate mixed quantum mechanics/molecular mechanics techniques and protein structure prediction methods, we provide a detailed electronic structure description of the protein-RNA stacking interactions. Our analysis indicates the evolution of the protein structure optimizing the interaction between Asp-92 and the RNA bases. The results show a direct coupling of the C-terminal tail and Asp-92, providing a direct rationalization of the experimentally determined role of the C-terminal domain in RNA binding. Here, we propose a mechanism where a protein side chain, with a delocalized electronic pi system, assists in the nucleotide binding. The binding mechanism involves a short-range interaction of the side chain with the nucleotide base and an electronic long-range interaction through a sandwich-stacking motif. The structural motif of the binding mechanism is observed in similar protein-RNA interactions and in various protein-ATP-binding sites.

Altering the RNA-binding mode of the U1A RBD1 protein

The N-terminal RNA-binding domain (RBD1) of the human U1A protein is evolutionarily designed to bind its RNA targets with great affinity and specificity. The physical mechanisms that modulate the coupling (local cooperativity) among amino acid residues on the extensive binding surface of RBD1 are investigated here, using mutants that replace a highly conserved glycine residue. This glycine residue, at the strand/loop junction of beta3/loop3, is found in U1A RBD1, and in most RBD domains, suggesting it has a specific role in modulation of RNA binding. Here, two RBD1 proteins are constructed in which that residue (Gly53) is replaced by either alanine or valine. These new proteins are shown by NMR methods and molecular dynamics simulations to be very similar to the wild-type RBD1, both in structure and in their backbone dynamics. However, RNA-binding assays show that affinity for the U1 snRNA stem-loop II RNA target is reduced by nearly 200-fold for the RBD1-G53A protein, and by 1.6 x 10(4)-fold for RBD1-G53V. The mode of RNA binding by RBD1-G53A is similar to that of RBD1-WT, displaying its characteristic non-additive free energies of base recognition and its salt-dependence. The binding mode of RBD1-G53V is altered, having lost its salt-dependence and displaying site-independence of base recognition. The molecular basis for this alteration in RNA-binding properties is proposed to result from the inability of the RNA to induce a change in the structure of the free protein to produce a high-affinity complex.

Determinants within an 18-amino-acid U1A autoregulatory domain that uncouple cooperative RNA binding, inhibition of polyadenylation, and homodimerization

The human U1 snRNP-specific U1A protein autoregulates its own production by binding to and inhibiting the polyadenylation of its own pre-mRNA. Previous work demonstrated that a short sequence of U1A protein is essential for autoregulation and contains three distinct activities, which are (i) cooperative binding of two U1A proteins to a 50-nucleotide region of U1A pre-mRNA called polyadenylation-inhibitory element RNA, (ii) formation of a novel homodimerization surface, and (iii) inhibition of polyadenylation by inhibition of poly(A) polymerase (PAP). In this study, we purified and analyzed 11 substitution mutant proteins, each having one or two residues in this region mutated. In 5 of the 11 mutant proteins, we found that particular amino acids associate with one activity but not another, indicating that they can be uncoupled. Surprisingly, in three mutant proteins, these activities were improved upon, suggesting that U1A autoregulation is selected for suboptimal inhibitory efficiency. The effects of these mutations on autoregulatory activity in vivo were also determined. Only U1A and U170K are known to regulate nuclear polyadenylation by PAP inhibition; thus, these results will aid in determining how widespread this type of regulation is. Our molecular dissection of the consequences of conformational changes within an RNP complex presents a powerful example to those studying more complicated pre-mRNA-regulatory systems.

Molecular dynamics simulation studies of induced fit and conformational capture in U1A-RNA binding: do molecular substates code for specificity?

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions.

Substitution of an essential adenine in the U1A-RNA complex with a non-polar isostere

The RNA recognition motif (RRM) binds to single-stranded RNA target sites of diverse sequences and structures. A conserved mode of base recognition by the RRM involves the simultaneous formation of a network of hydrogen bonds with the base functional groups and a stacking interaction between the base and a highly conserved aromatic amino acid. We have investigated the energetic contribution of the functional groups involved in the recognition of an essential adenine, A6, in stem-loop 2 of U1 snRNA by the N-terminal RRM of the U1A protein. Previously, we found that elimination of individual hydrogen bond donors and acceptors on A6 destabilized the complex by 0.8-1.9 kcal/mol, while mutation of the aromatic amino acid (Phe56) that stacks with A6 to Ala destabilized the complex by 5.5 kcal/mol. Here we continue to probe the contribution of A6 to complex stability through mutation of both the RNA and protein. We have removed two hydrogen-bonding functional groups by introducing a U1A mutation, Ser91Ala, and replacing A6 with tubercidin, purine, or 1-deazaadenine. We find that the complex is destabilized an additional 1.2-2.6 kcal/mol by the elimination of the second hydrogen bond donor or acceptor. Surprisingly, deletion of all of the functional groups involved in hydrogen bonds with the U1A protein by substituting adenine with 4-methylindole reduced the binding free energy by only 2.0 kcal/mol. Experiments with U1A proteins containing mutations of Phe56 suggested that improved stacking interactions due to the greater hydrophobicity of 4-methylindole than adenine may be partly responsible for the small destabilization of the complex upon substitution of 4-methylindole for A6. The data imply that hydrophobic interactions can compensate energetically for the disruption of the complex hydrogen-bonding network between nucleotide and protein.

Intramolecular interactions at protein surfaces and their impact on protein function

Protein surfaces play a key role in the biological function of proteins. Consequently, structural features of protein surfaces are the basis for predicting function from structure. A well-established principle of binding by proteins is that ligands must compete with water and other small molecules to form interactions with protein surfaces. A less obvious issue, and the emphasis of this article, is that ligands must also compete with interactions among residues at protein surfaces. Results from structural surveys, a variety of experimental studies and computations suggest that intramolecular interactions are present at protein surfaces and that the energetics of these interactions can change when proteins bind to other molecules.

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.

Complex role of the beta 2-beta 3 loop in the interaction of U1A with U1 hairpin II RNA

RNA recognition motifs (RRMs) are characterized by highly conserved regions located centrally on a beta-sheet, which forms the RNA binding surface. Variable flanking regions, such as the loop connecting beta-strands 2 and 3, are thought to be important in determining the RNA-binding specificities of individual RRMs. The N-terminal RRM of the spliceosomal U1A protein mediates binding to an RNA hairpin (U1hpII) in the U1 small nuclear RNA. In this complex, the beta(2)-beta(3) loop protrudes through the 10-nucleotide RNA loop. Shortening of the RNA loop strongly perturbs binding, suggesting that an optimal "fit" of the beta(2)-beta(3) loop into the RNA loop is an important factor in complexation. To understand this interaction further, we mutated or deleted loop residues Lys(50) and Met(51), which protrude centrally into the RNA loop but do not make any direct contacts to the bases. Using BIACORE, we analyzed the ability of these U1A mutants to bind to wild type RNAs, or RNAs with shortened loops. Alanine replacement mutations only modestly affected binding to wild type U1hpII. Interestingly, simultaneous replacement of Lys(50) and Met(51) with alanine appeared to alleviate the loss of binding caused by shortening of the RNA loop. Deletion of Lys(50) or Met(51) caused a dramatic loss in stability of the U1A.U1hpII complex. However, deletion of both residues simultaneously was much less deleterious. Simulated annealing molecular dynamics analyses suggest this is due to the ability of this mutant to rearrange flanking amino acids to substitute for the two deleted residues. The double deletion mutant also exhibited substantially reduced negative effects of RNA loop shortening, suggesting the rearranged loop is better able to accommodate a short RNA loop. Our results indicate that one of the roles of the beta(2)-beta(3) loop is to provide a steric fit into the RNA loop, thereby stabilizing the RNA.protein complex.

Anti-U1A monoclonal antibodies recognize unique epitope targets of U1A which are involved in the binding of U1 RNA

The U1A (or nRNP A) protein is known to play a critical role in eukaryotic pre-mRNA splicing and polyadenylation. Previous studies revealed that several mouse monoclonal antibodies (MAbs) recognized U1A as part of the U1snRNP, while MAb 12E12 was unique in that it recognized an epitope that is masked when U1A is bound to U1 RNA. In order to further characterize and understand the antigenic targets of these MAbs, we undertook fine specificity epitope mapping studies. Anti-U1A MAbs 12E12 and 10E3 each recognize unique peptides from the U1A protein. Interestingly, these MAbs recognize epitopes which have been shown to be antigenic in human autoimmune diseases. When superimposed on structures of U1A derived from crystal and NMR data, the major epitope recognized by 12E12 (amino acids 103-108) localizes to the surface of the U1A molecule. The 12E12 epitope is immediately adjacent to a helix which probably reacts to U1 RNA binding by undergoing a conformational change. This modification of structure effectively masks the 12E12 epitope, thus preventing binding of the monoclonal to U1A/U1 RNA complexes. These findings suggest that the structure of the U1A protein may be different when not part of the U1snRNP.

Pre-steady-state and stopped-flow fluorescence analysis of Escherichia coli ribonuclease III: insights into mechanism and conformational changes associated with binding and catalysis

To better understand substrate recognition and catalysis by RNase III, we examined steady-state and pre-steady-state reaction kinetics, and changes in intrinsic enzyme fluorescence. The multiple turnover cleavage of a model RNA substrate shows a pre-steady-state burst of product formation followed by a slower phase, indicating that the steady-state reaction rate is not limited by substrate cleavage. RNase III catalyzed hydrolysis is slower at low pH, permitting the use of pre-steady-state kinetics to measure the dissociation constant for formation of the enzyme-substrate complex (K(d)=5.4(+/-0.6) nM), and the rate constant for phosphodiester bond cleavage (k(c)=1.160(+/-0.001) min(-1), pH 5.4). Isotope incorporation analysis shows that a single solvent oxygen atom is incorporated into the 5' phosphate of the RNA product, which demonstrates that the cleavage step is irreversible. Analysis of the pH dependence of the single turnover rate constant, k(c), fits best to a model for two or more titratable groups with pK(a) of ca 5.6, suggesting a role for conserved acidic residues in catalysis. Additionally, we find that k(c) is dependent on the pK(a) value of the hydrated divalent metal ion included in the reaction, providing evidence for participation of a metal ion hydroxide in catalysis, potentially in developing the nucleophile for the hydrolysis reaction. In order to assess whether conformational changes also contribute to the enzyme mechanism, we monitored intrinsic tryptophan fluorescence. During a single round of binding and cleavage by the enzyme we detect a biphasic change in fluorescence. The rate of the initial increase in fluorescence was dependent on substrate concentration yielding a second-order rate constant of 1.0(+/-0.1)x10(8) M(-1) s(-1), while the rate constant of the second phase was concentration independent (6.4(+/-0.8) s(-1); pH 7.3). These data, together with the unique dependence of each phase on divalent metal ion identity and pH, support the hypothesis that the two fluorescence transitions, which we attribute to conformational changes, correlate with substrate binding and catalysis.

Investigation of a conserved stacking interaction in target site recognition by the U1A protein

Three highly conserved aromatic residues in RNA recognition motifs (RRM) participate in stacking interactions with RNA bases upon binding RNA. We have investigated the contribution of one of these aromatic residues, Phe56, to the complex formed between the N-terminal RRM of the spliceosomal protein U1A and stem-loop 2 of U1 snRNA. Previous work showed that the aromatic group is important for high affinity binding. Here we probe how mutation of Phe56 affects the kinetics of complex dissociation, the strength of the hydrogen bonds formed between U1A and the base that stacks with Phe56 (A6) and specific target site recognition. Substitution of Phe56 with Trp or Tyr increased the rate of dissociation of the complex, consistent with previously reported results. However, substitution of Phe56 with His decreased the rate of complex association, implying a change in the initial formation of the complex. Simultaneous modification of residue 56 and A6 revealed energetic coupling between the aromatic group and the functional groups of A6 that hydrogen bond to U1A. Finally, mutation of Phe56 to Leu reduced the ability of U1A to recognize stem-loop 2 correctly. Taken together, these experiments suggest that Phe56 contributes to binding affinity by stacking with A6 and participating in networks of energetically coupled interactions that enable this conserved aromatic amino acid to play a complex role in target site recognition.

Calculations of free-energy contributions to protein-RNA complex stabilization

The problem of calculating binding affinities of protein-RNA complexes is addressed by analyzing a computational strategy of modeling electrostatic free energies based on a nonlinear Poisson-Boltzmann (NLPB) model and linear response approximation (LRA). The underlying idea is to treat binding as a two-step process. Solutions to the NLPB equation calculate free energies arising from electronic polarizability and the LRA is constructed from molecular dynamics simulations to model reorganization free energies due to conformational transitions. By implementing a consistency condition of requiring the NLPB model to reproduce the solute-solvent free-energy transitions determined by the LRA, a "macromolecule dielectric constant" (epsilon(m)) for treating reorganization is obtained. The applicability of this hybrid approach was evaluated by calculating the absolute free energy of binding and free-energy changes for amino acid substitutions in the complex between the U1A spliceosomal protein and its cognate RNA hairpin. Depending on the residue substitution, epsilon(m) varied from 3 to 18, and reflected dipolar reorientation not included in the polarization modeled by epsilon(m) = 2. Although the changes in binding affinities from substitutions modeled strictly at the implicit level by the NLPB equation with epsilon(m) = 4 reproduced the experimental values with good overall agreement, substitutions problematic to this simple treatment showed significant improvement when solved by the NLPB-LRA approach.

Two functionally distinct steps mediate high affinity binding of U1A protein to U1 hairpin II RNA

Binding of the U1A protein to its RNA target U1 hairpin II has been extensively studied as a model for a high affinity RNA/protein interaction. However, the mechanism and kinetics by which this complex is formed remain largely unknown. Here we use real-time biomolecular interaction analysis to dissect the roles various protein and RNA structural elements play in the formation of the U1A.U1 hairpin II complex. We show that neutralization of positive charges on the protein or increasing the salt concentration slows the association rate, suggesting that electrostatic interactions play an important role in bringing RNA and protein together. In contrast, removal of hydrogen bonding or stacking interactions within the RNA/protein interface, or reducing the size of the RNA loop, dramatically destabilizes the complex, as seen by a strong increase in the dissociation rate. Our data support a binding mechanism consisting of a rapid initial association based on electrostatic interactions and a subsequent locking step based on close-range interactions that occur during the induced fit of RNA and protein. Remarkably, these two steps can be clearly distinguished using U1A mutants containing single amino acid substitutions. Our observations explain the extraordinary affinity of U1A for its target and may suggest a general mechanism for high affinity RNA/protein interactions.

Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex

Molecular dynamics (MD) simulations and free energy component analysis have been performed to evaluate the molecular origins of the 5.5 kcal/mol destabilization of the complex formed between the N-terminal RNP domain of U1A and stem loop 2 of U1 snRNA upon mutation of a conserved aromatic residue, Phe56, to Ala. MD simulations, including counterions and water, have been carried out on the wild type and Phe56Ala peptide-stem loop 2 RNA complexes, the free wild type and Phe56Ala peptides, and the free stem loop 2 RNA. The MD structure of the Phe56Ala-stem loop 2 complex is similar to that of the wild type complex except the stacking interaction between Phe56 and A6 of stem loop 2 is absent and loop 3 of the peptide is more dynamic. However, the MD simulations predict large changes in the structure and dynamics of helix C and increased dynamic range of loop 3 for the free Phe56Ala peptide compared to the wild type peptide. Since helix C and loop 3 are highly variable regions of RNP domains, this indicates that a significant contribution to the reduced affinity of the Phe56Ala peptide for RNA results from cooperation between highly conserved and highly variable regions of the RNP domain of U1A. Surprisingly, these structural effects, which are manifested as cooperative free energy changes, occur in the free peptide, rather than in the complex, and are revealed only by study of both the initial and final states of the complexation process. Free energy component analysis correctly accounts for the destabilization of the Phe56Ala-stem loop 2 complex, and indicates that approximately 80% of the destabilization is due to the loss of the stacking interaction and approximately 20% is due to differences in U1A adaptation.

Recent advances in RNA-protein recognition

The past few years have witnessed remarkable progress in knowledge of the structure and function of RNA-binding proteins and their RNA complexes. X-ray crystallography and NMR spectroscopy have provided structures for all major classes of RNA-binding proteins, both alone and complexed with RNA. New computational and experimental tools have provided unprecedented insight into the molecular basis of RNA recognition.

Structure and thermodynamics of RNA-protein binding: using molecular dynamics and free energy analyses to calculate the free energies of binding and conformational change

An adaptive binding mechanism, requiring large conformational rearrangements, occurs commonly with many RNA-protein associations. To explore this process of reorganization, we have investigated the conformational change upon spliceosomal U1A-RNA binding with molecular dynamics (MD) simulations and free energy analyses. We computed the energetic cost of conformational change in U1A-hairpin and U1A-internal loop binding using a hybrid of molecular mechanics and continuum solvent methods. Encouragingly, in all four free energy comparisons (two slightly different proteins, two different RNAs), the free macromolecule was more stable than the bound form by the physically reasonable value of approximately 10 kcal/mol. We calculated the absolute binding free energies for both complexes to be in the same range as that found experimentally.

Investigating the binding specificity of U1A-RNA by computational mutagenesis

The mammalian spliceosomal protein U1A binds a hairpin RNA with picomolar affinity. To examine the origin of this binding specificity, we carried out computational mutagenesis on protein and RNA residues in the U1A-RNA binding interface. Our computational mutagenesis methods calculate the relative binding affinity between mutant and wild-type as the sum of molecular mechanical energies and solvation free energies estimated with a continuum solvent model. We obtained good agreement with experimental studies and we verified mutations that abolish and improve binding. Therefore, we offer these methods as computationally inexpensive tools for investigating and predicting the effects of site-specific mutagenesis.

Calorimetric examination of high-affinity Src SH2 domain-tyrosyl phosphopeptide binding: dissection of the phosphopeptide sequence specificity and coupling energetics

SH2 domains are protein modules which interact with specific tyrosine phosphorylated sequences in target proteins. The SH2 domain of the Src kinase binds with high affinity to a tyrosine phosphorylated peptide containing the amino acids Glu, Glu, and Ile (EEI) at the positions +1, +2, and +3 C-terminal to the phosphotyrosine, respectively. To investigate the degree of selectivity of the Src SH2 domain for each amino acid of the EEI motif, the binding thermodynamics of a panel of substitutions at the +1 (Gln, Asp, Ala, Gly), +2 (Gln, Asp, Ala, Gly), and +3 (Leu, Val, Ala, Gly) positions were examined using titration microcalorimetry. It was revealed that the Src SH2 domain is insensitive (DeltaDeltaG degrees </= 0.6 kcal/mol) to conservative substitutions at all three peptide positions. However, mutation to Ala resulted in moderate reductions in DeltaG degrees, with the substitution at the +3 position showing the largest loss in affinity (DeltaDeltaG degrees = 1.4 kcal/mol), followed by the +2 (DeltaDeltaG degrees = 1.0 kcal/mol) and +1 (DeltaDeltaG degrees = 0.5 kcal/mol) positions. This hierarchy of binding was not reflected in the values of the heat capacity change, since only the peptide substituted to Ala at the +3 position showed a DeltaCp degrees that was reduced in magnitude compared to wild-type. To assess the degree of cooperation upon binding (or coupling) between the amino acids of the EEI sequence, the binding of a series of singly, doubly, and triply Ala substituted phosphopeptides was examined and analyzed using double mutant cycles. It was revealed that the effects of the Ala substitutions on DeltaG degrees were additive. However, nonadditive binding enthalpies were observed between the +1 Glu and +3 Ile, as well as the +2 Glu and +3 Ile, suggesting that communication occurs between residues of the EEI motif upon binding.

RNA recognition by the human U1A protein is mediated by a network of local cooperative interactions that create the optimal binding surface

One of the most common structural motifs in RNA-binding proteins is the RNA-binding domain (RBD). These domains share a common alpha/beta sandwich tertiary fold, and are highly conserved, though they bind diverse RNA targets with a wide range of binding affinities. The N-terminal RNA-binding domain (RBD1) of the human U1A protein binds specifically to stem/loop II of the U1 snRNA with sub-nanomolar affinity. Solvent-exposed aromatic residues on the beta-sheet surface are highly conserved among RBD domains; in RBD1, these are Tyr13 and Phe56, with a unique Gln at position 54. Effects of substitutions at these positions were examined using energetic pairwise coupling to describe the communication between these residues in both the free and RNA-bound states of the protein. 15N NMR experiments were used to determine effects of the beta-sheet substitutions on the structural and dynamic properties of this domain. The combination of thermodynamic pairwise coupling and 15N-backbone dynamics provides direct evidence for local cooperative interactions among Y13, Q54, and F56, and a non-conserved loop that directly affect RNA-binding. The results describe how conserved and non-conserved regions of an RBD can communicate with each other to mediate recognition of the RNA.