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

Affinity and Structural Analysis of the U1A RNA Recognition Motif with Engineered Methionines to Improve Experimental Phasing

RNA plays a central role in all organisms and can fold into complex structures to orchestrate function. Visualization of such structures often requires crystallization, which can be a bottleneck in the structure-determination process. To promote crystallization, an RNA-recognition motif (RRM) of the U1A spliceosomal protein has been co-opted as a crystallization module. Specifically, the U1-snRNA hairpin II (hpII) single-stranded loop recognized by U1A can be transplanted into an RNA target to promote crystal contacts and to attain phase information via molecular replacement or anomalous diffraction methods using selenomethionine. Herein, we produced the F37M/F77M mutant of U1A to augment the phasing capability of this powerful crystallization module. Selenomethionine-substituted U1A(F37M/F77M) retains high affinity for hpII (K _D of 59.7 ± 11.4 nM). The 2.20 Å resolution crystal structure reveals that the mutated sidechains make new S-π interactions in the hydrophobic core and are useful for single-wavelength anomalous diffraction. Crystals were also attained of U1A(F37M/F77M) in complex with a bacterial preQ₁-II riboswitch. The F34M/F37M/F77M mutant was introduced similarly into a lab-evolved U1A variant (TBP6.9) that recognizes the internal bulged loop of HIV-1 TAR RNA. We envision that this short RNA sequence can be placed into non-essential duplex regions to promote crystallization and phasing of target RNAs. We show that selenomethionine-substituted TBP6.9(F34M/F37M/F77M) binds a TAR variant wherein the apical loop was replaced with a GNRA tetraloop (K _D of 69.8 ± 2.9 nM), laying the groundwork for use of TBP6.9(F34M/F37M/F77M) as a crystallization module. These new tools are available to the research community.

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.

Role of Electrostatics in Protein-RNA Binding: The Global vs the Local Energy Landscape

U1A protein-stem loop 2 RNA association is a basic step in the assembly of the spliceosomal U1 small nuclear ribonucleoprotein. Long-range electrostatic interactions due to the positive charge of U1A are thought to provide high binding affinity for the negatively charged RNA. Short range interactions, such as hydrogen bonds and contacts between RNA bases and protein side chains, favor a specific binding site. Here, we propose that electrostatic interactions are as important as local contacts in biasing the protein-RNA energy landscape toward a specific binding site. We show by using molecular dynamics simulations that deletion of two long-range electrostatic interactions (K22Q and K50Q) leads to mutant-specific alternative RNA bound states. One of these states preserves short-range interactions with aromatic residues in the original binding site, while the other one does not. We test the computational prediction with experimental temperature-jump kinetics using a tryptophan probe in the U1A-RNA binding site. The two mutants show the distinct predicted kinetic behaviors. Thus, the stem loop 2 RNA has multiple binding sites on a rough RNA-protein binding landscape. We speculate that the rough protein-RNA binding landscape, when biased to different local minima by electrostatics, could be one way that protein-RNA interactions evolve toward new binding sites and novel function.

An Evolved RNA Recognition Motif That Suppresses HIV-1 Tat/TAR-Dependent Transcription

Potent and selective recognition and modulation of disease-relevant RNAs remain a daunting challenge. We previously examined the utility of the U1A N-terminal RNA recognition motif as a scaffold for tailoring new RNA hairpin recognition and showed that as few as one or two mutations can result in moderate affinity (low μM dissociation constant) for the human immunodeficiency virus (HIV) trans-activation response element (TAR) RNA, an RNA hairpin controlling transcription of the human immunodeficiency virus (HIV) genome. Here, we use yeast display and saturation mutagenesis of established RNA-binding regions in U1A to identify new synthetic proteins that potently and selectively bind TAR RNA. Our best candidate has truly altered, not simply broadened, RNA-binding selectivity; it binds TAR with subnanomolar affinity (apparent dissociation constant of ∼0.5 nM) but does not appreciably bind the original U1A RNA target (U1hpII). It specifically recognizes the TAR RNA hairpin in the context of the HIV-1 5'-untranslated region, inhibits the interaction between TAR RNA and an HIV trans-activator of transcription (Tat)-derived peptide, and suppresses Tat/TAR-dependent transcription. Proteins described in this work are among the tightest TAR RNA-binding reagents-small molecule, nucleic acid, or protein-reported to date and thus have potential utility as therapeutics and basic research tools. Moreover, our findings demonstrate how a naturally occurring RNA recognition motif can be dramatically resurfaced through mutation, leading to potent and selective recognition-and modulation-of disease-relevant RNA.

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.

Synthetic RNA recognition motifs that selectively recognize HIV-1 trans-activation response element hairpin RNA

A multitude of RNA hairpins are directly implicated in human disease. Many of these RNAs are potentially valuable targets for drug discovery and basic research. However, very little is known about the molecular requirements for achieving sequence-selective recognition of a particular RNA sequence and structure. Although a relatively modest number of synthetic small to medium-sized RNA-binding molecules have been reported, rapid identification of sequence-selective RNA-binding molecules remains a daunting challenge. RNA recognition motif (RRM) domains may represent unique privileged scaffolds for the generation of synthetic proteins that selectively recognize structured disease-relevant RNAs, including RNA hairpins. As a demonstration of this potential, we mutated putative RNA-binding regions within the U1A RRM and a variant thereof and screened these synthetic proteins for affinity to HIV-1 trans-activation response (TAR) element hairpin RNA. Some of these U1A-derived proteins bind TAR with single-digit micromolar dissociation constants, and they do so preferentially over the native protein's original target RNA (U1hpII) and a DNA TAR variant. Binding affinity is not appreciably diminished by addition of 10 molar equivalents of cellular tRNAs from Escherichia coli. Taken together, our findings represent the first synthetic RRMs that selectively bind a disease-relevant RNA hairpin and may represent a general approach for achieving sequence-selective recognition of RNA hairpins, which are the focus of therapeutic discovery and basic research.

Analysis of protein-RNA complexes involving a RNA recognition motif engineered to bind hairpins with seven- and eight-nucleotide loops

U1A binds U1hpII, a hairpin RNA with a 10-nucleotide loop. A U1A mutant (ΔK50ΔM51) binds U1hpII-derived hairpins with shorter loops, making it an interesting scaffold for engineering or evolving proteins that bind similarly sized disease-related hairpin RNAs. However, a more detailed understanding of complexes involving ΔK50ΔM51 is likely a prerequisite to generating such proteins. Toward this end, we measured mutational effects for complexes involving U1A ΔK50ΔM51 and U1hpII-derived hairpin RNAs with seven- or eight-nucleotide loops and identified contacts that are critical to the stabilization of these complexes. Our data provide valuable insight into sequence-selective recognition of seven- or eight-nucleotide loop hairpins by an engineered RNA binding protein.

Structure of the heterodimer of human NONO and paraspeckle protein component 1 and analysis of its role in subnuclear body formation

Proteins of the Drosophila behavior/human splicing (DBHS) family include mammalian SFPQ (PSF), NONO (p54nrb), PSPC1, and invertebrate NONA and Hrp65. DBHS proteins are predominately nuclear, and are involved in transcriptional and posttranscriptional gene regulatory functions as well as DNA repair. DBHS proteins influence a wide gamut of biological processes, including the regulation of circadian rhythm, carcinogenesis, and progression of cancer. Additionally, mammalian DBHS proteins associate with the architectural long noncoding RNA NEAT1 (Menε/β) to form paraspeckles, subnuclear bodies that alter gene expression via the nuclear retention of RNA. Here we describe the crystal structure of the heterodimer of the multidomain conserved region of the DBHS proteins, PSPC1 and NONO. These proteins form an extensively intertwined dimer, consistent with the observation that the different DBHS proteins are typically copurified from mammalian cells, and suggesting that they act as obligate heterodimers. The PSPC1/NONO heterodimer has a right-handed antiparallel coiled-coil that positions two of four RNA recognition motif domains in an unprecedented arrangement on either side of a 20-Å channel. This configuration is supported by a protein:protein interaction involving the NONA/paraspeckle domain, which is characteristic of the DBHS family. By examining various mutants and truncations in cell culture, we find that DBHS proteins require an additional antiparallel coiled-coil emanating from either end of the dimer for paraspeckle subnuclear body formation. These results suggest that paraspeckles may potentially form through self-association of DBHS dimers into higher-order structures.

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.

Structure and novel functional mechanism of Drosophila SNF in sex-lethal splicing

Sans-fille (SNF) is the Drosophila homologue of mammalian general splicing factors U1A and U2B'', and it is essential in Drosophila sex determination. We found that, besides its ability to bind U1 snRNA, SNF can also bind polyuridine RNA tracts flanking the male-specific exon of the master switch gene Sex-lethal (Sxl) pre-mRNA specifically, similar to Sex-lethal protein (SXL). The polyuridine RNA binding enables SNF directly inhibit Sxl exon 3 splicing, as the dominant negative mutant SNF(1621) binds U1 snRNA but not polyuridine RNA. Unlike U1A, both RNA recognition motifs (RRMs) of SNF can recognize polyuridine RNA tracts independently, even though SNF and U1A share very high sequence identity and overall structure similarity. As SNF RRM1 tends to self-associate on the opposite side of the RNA binding surface, it is possible for SNF to bridge the formation of super-complexes between two introns flanking Sxl exon 3 or between a intron and U1 snRNP, which serves the molecular basis for SNF to directly regulate Sxl splicing. Taken together, a new functional model for SNF in Drosophila sex determination is proposed. The key of the new model is that SXL and SNF function similarly in promoting Sxl male-specific exon skipping with SNF being an auxiliary or backup to SXL, and it is the combined dose of SXL and SNF governs Drosophila sex determination.

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.

Apoptosis-linked changes in the phosphorylation status and subcellular localization of the spliceosomal autoantigen U1-70K

Apoptosis consists of highly regulated pathways involving post-translational modifications and cleavage of proteins leading to sequential inactivation of the main cellular processes. Here, we focused on the apoptotic processing of one of the essential components of the mRNA splicing machinery, the U1-70K snRNP protein. We found that at an early stage of apoptosis, before the cleavage of the C-terminal part of the protein by caspase-3, the basal phosphorylation of the Ser140 residue located within the RNA recognition motif, increases very significantly. A caspase-dependent, PP1-mediated dephosphorylation of other serine residues takes place in a subset of U1-70K proteins. The U1-70K protein phosphorylated at Ser140 is clustered in heterogeneous ectopic RNP-derived structures, which are finally extruded in apoptotic bodies. The elaborate processing of the spliceosomal U1-70K protein we identified might play an important role in the regulated breakdown of the mRNA splicing machinery during early apoptosis. In addition, these specific changes in the phosphorylation/dephosphorylation balance and the subcellular localization of the U1-70K protein might explain why the region encompassing the Ser140 residue becomes a central autoantigen during the autoimmune disease systemic lupus erythematosus.

Biosensor-surface plasmon resonance methods for quantitative analysis of biomolecular interactions

The surface plasmon resonance (SPR) biosensor method has emerged as a very flexible and powerful approach for detecting a wide diversity of biomolecular interactions. SPR monitors molecular interactions in real time and provides significant advantages over optical or calorimetric methods for systems with strong binding and low spectroscopic signals or reaction heats. The SPR method simultaneously provides kinetic and equilibrium characterization of the interactions of biomolecules. Such information is essential for development of a full understanding of molecular recognition as well as for areas such as the design of receptor-targeted therapeutics. This article presents basic, practical procedures for conducting SPR experiments. Initial preparation of the SPR instrument, sensor chips, and samples are described. This is followed by suggestions for experimental design, data analysis, and presentation. Steady-state and kinetic studies of some small molecule-DNA complexes are used to illustrate the capability of this technique. Examples of the agreement between biosensor-SPR and solution studies are presented.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

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.

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.

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.

A survey of the year 2002 commercial optical biosensor literature

We have compiled 819 articles published in the year 2002 that involved commercial optical biosensor technology. The literature demonstrates that the technology's application continues to increase as biosensors are contributing to diverse scientific fields and are used to examine interactions ranging in size from small molecules to whole cells. Also, the variety of available commercial biosensor platforms is increasing and the expertise of users is improving. In this review, we use the literature to focus on the basic types of biosensor experiments, including kinetics, equilibrium analysis, solution competition, active concentration determination and screening. In addition, using examples of particularly well-performed analyses, we illustrate the high information content available in the primary response data and emphasize the impact of including figures in publications to support the results of biosensor analyses.