A quantitative study of the flexibility contributed to RNA structures by nicks and single-stranded gaps.

Disulfide crosslinking via thiol-disulfide interchange was applied to quantitate the relative flexibility contributed by nicks and single-stranded gaps in an RNA structure. An RNA duplex comprised of three strands was constructed containing the disulfide crosslink precursors 1 and 2 at opposite ends of the duplex on opposite strands. The third strand was of varying length to yield a nick or single-stranded gaps of 1, 2, or 3 nt. Crosslinking rates Indicated relative flexibilities of the resulting two-helix junctions. Crosslinking in the nicked duplex occurred two orders of magnitude slower than in a duplex containing a 3-nt gap. Rates of crosslinking in duplexes with 3-and 2-nt gaps showed only modest dependence on the gap sequence. Many natural RNAs, including ribozymes, contain two-helix junctions related to the model system described here. The data suggest that two-helix junctions containing a nick in one strand will retain substantial rigidity, whereas one or more single-stranded nucleotides at a two-helix junction allow significant flexibility.     

Solution structure of the helicase-interaction domain of the primase DnaG: a model for helicase activation

The helicase-primase interaction is a critical event in DNA replication and is mediated by a putative helicase-interaction domain within the primase. The solution structure of the helicase-interaction domain of DnaG reveals that it is made up of two independent subdomains: an N-terminal six-helix module and a C-terminal two-helix module that contains the residues of the primase previously identified as important in the interaction with the helicase. We show that the two-helix module alone is sufficient for strong binding between the primase and the helicase but fails to activate the helicase; both subdomains are required for helicase activation. The six-helix module of the primase has only one close structural homolog, the N-terminal domain of the corresponding helicase. This surprising structural relationship, coupled with the differences in surface properties of the two molecules, suggests how the helicase-interaction domain may perturb the structure of the helicase and lead to activation.     

Two helix DNA binding motif of CAP found in lac repressor and gal repressor.

Comparison of both the DNA and protein sequences of catabolite gene activator protein (CAP) with the sequences of lac and gal repressors shows significant homologies between a sequence that forms a two alpha-helix motif in CAP and sequences near the amino terminus of both repressors. This two-helix motif is thought to be involved in specific DNA sequence recognition by CAP. The region in lac repressor to which CAP is homologous contains many i-d mutations that are defective in DNA binding. Less significant sequence homologies between CAP and phage repressors and activators are also shown. The amino acid residues that are critical to the formation of the two-helix motif are conserved, while those residues expected to interact with DNA are variable. These observations suggest the lac and gal repressors also have a two alpha-helix structural motif which is involved in DNA binding and that this two helix motif may be generally found in many bacterial and phage repressors. We conclude that one major mechanism by which proteins can recognize specific base sequences in double stranded DNA is via the amino acid side chains of alpha-helices fitting into the major groove of B-DNA.     

Synthetic peptides mimicking the interleukin-6/gp130 interaction: a two-helix bundle system. Design and conformational studies.

The objective of our study was to mimic in a typical reductionist approach the molecular interactions observed at the interface between the gp130 receptor and interleukin-6 during formation of their complex. A peptide system obtained by reproducing some of the interleukin-6/gp130 molecular interactions into a two-helix bundle structure was investigated. The solution conformational features of this system were determined by CD and NMR techniques. The CD titration experiments demonstrated that the interaction between the designed peptides is specific and based on a well-defined stoichiometry. The NMR data confirmed some of the structural features of the binding mechanism as predicted by the rational design and indicated that under our conditions the recognition specificity and affinity can be explained by the formation of a two-helix bundle. Thus, the data reported herein represent a promising indication on how to develop new peptides able to interfere with formation of the interleukin-6/gp130 complex.     

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.     

Basic faced alpha-helices are widespread in the peptide extensions of the eukaryotic aminoacyl-tRNA synthetases

Three aminoacyl-tRNA synthetases from yeast, one from plants and one from mammals possess unusual structures at their N termini, namely alpha helices with basic residues distributed asymmetrically, on a single face of the helix. It is unknown if these 'basic faced' alpha helices (BFAHs) are unique to the aminoacyl-tRNA synthetases. Analysis of the amino acid sequences of these five aminoacyl-tRNA synthetases using the hydrophobic moment algorithm failed to accurately identify the BFAHs. A new algorithm was therefore developed, called the 'basic moment'. This is a Fourier analysis procedure that predicts the distribution of basic residues within protein secondary structure. The basic moment identifies with a high degree of accuracy the five known BFAHs and also identifies further potential BFAHs at evolutionarily conserved positions in the peptide extensions of aspartyl-, lysyl- and valyl- tRNA synthetases from a range of eukaryotic species. In addition, the algorithm identifies the two-helix pair tRNA binding domain of alanyl-tRNA synthetase, implying that the domain includes a BFAH. The functional and evolutionary aspects of these structural features are discussed.     

Different modes of interaction by TIAR and HuR with target RNA and DNA

TIAR and HuR are mRNA-binding proteins that play important roles in the regulation of translation. They both possess three RNA recognition motifs (RRMs) and bind to AU-rich elements (AREs), with seemingly overlapping specificity. Here we show using SPR that TIAR and HuR bind to both U-rich and AU-rich RNA in the nanomolar range, with higher overall affinity for U-rich RNA. However, the higher affinity for U-rich sequences is mainly due to faster association with U-rich RNA, which we propose is a reflection of the higher probability of association. Differences between TIAR and HuR are observed in their modes of binding to RNA. TIAR is able to bind deoxy-oligonucleotides with nanomolar affinity, whereas HuR affinity is reduced to a micromolar level. Studies with U-rich DNA reveal that TIAR binding depends less on the 2'-hydroxyl group of RNA than HuR binding. Finally we show that SAXS data, recorded for the first two domains of TIAR in complex with RNA, are more consistent with a flexible, elongated shape and not the compact shape that the first two domains of Hu proteins adopt upon binding to RNA. We thus propose that these triple-RRM proteins, which compete for the same binding sites in cells, interact with their targets in fundamentally different ways.     

RNA binding specificity of Drosophila muscleblind

Members of the muscleblind family of RNA binding proteins found in Drosophila and mammals are key players in both the human disease myotonic dystrophy and the regulation of alternative splicing. Recently, the mammalian muscleblind-like protein, MBNL1, has been shown to have interesting RNA binding properties with both endogenous and disease-related RNA targets. Here we report the characterization of RNA binding properties of the Drosophila muscleblind protein Mbl. Mutagenesis of double-stranded CUG repeats demonstrated that Mbl requires pyrimidine-pyrimidine mismatches for binding and that the identity and location of the C-G and G-C base pairs within the repeats are essential for Mbl binding. Systematic evolution of ligands by exponential enrichment (SELEX) was used to identify RNA sequences that bind Mbl with much higher affinity than CUG repeats. The RNA sequences identified by SELEX are structured and contain a five-nucleotide consensus sequence of 5'-AGUCU-3'. RNase footprinting of one of the SELEX RNA sequences with Mbl showed that Mbl binds both double-stranded and single-stranded regions of the RNA. Three guanosines show the strongest footprint in the presence of Mbl; mutation of any of these three guanosines eliminates Mbl binding. It was also found that Mbl specifically bound a human MBNL1 RNA target, demonstrating the conservation of the muscleblind proteins in recognizing RNA targets. Our results reveal that Mbl recognizes complex RNA secondary structures.     

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

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

Analysis of RNA binding by the dengue virus NS5 RNA capping enzyme

Flaviviruses are small, capped positive sense RNA viruses that replicate in the cytoplasm of infected cells. Dengue virus and other related flaviviruses have evolved RNA capping enzymes to form the viral RNA cap structure that protects the viral genome and directs efficient viral polyprotein translation. The N-terminal domain of NS5 possesses the methyltransferase and guanylyltransferase activities necessary for forming mature RNA cap structures. The mechanism for flavivirus guanylyltransferase activity is currently unknown, and how the capping enzyme binds its diphosphorylated RNA substrate is important for deciphering how the flavivirus guanylyltransferase functions. In this report we examine how flavivirus NS5 N-terminal capping enzymes bind to the 5' end of the viral RNA using a fluorescence polarization-based RNA binding assay. We observed that the K(D) for RNA binding is approximately 200 nM Dengue, Yellow Fever, and West Nile virus capping enzymes. Removal of one or both of the 5' phosphates reduces binding affinity, indicating that the terminal phosphates contribute significantly to binding. RNA binding affinity is negatively affected by the presence of GTP or ATP and positively affected by S-adensyl methoninine (SAM). Structural superpositioning of the dengue virus capping enzyme with the Vaccinia virus VP39 protein bound to RNA suggests how the flavivirus capping enzyme may bind RNA, and mutagenesis analysis of residues in the putative RNA binding site demonstrate that several basic residues are critical for RNA binding. Several mutants show differential binding to 5' di-, mono-, and un-phosphorylated RNAs. The mode of RNA binding appears similar to that found with other methyltransferase enzymes, and a discussion of diphosphorylated RNA binding is presented.     

Design and analysis of RNA structure-specific agents as potential antivirals

A number of pathogenic RNA viruses, such as HIV-1, have extensive folded RNA conformations with imperfect A-form duplexes that are essential for virus function, and could serve as targets for structure-specific antiviral drugs. A method for the discovery of such drugs involves evaluation of the interactions with RNA of a wide variety of compounds that are known to bind to nucleic acids by different mechanisms. This approach has been initiated by using corresponding sequence RNA and DNA polymers as initial test systems for analysis of RNA binding strength and selectivity. Compounds that bind exclusively in the minor groove in AT sequences of DNA do not have significant interactions with RNA. Polycations, however, can show significant RNA affinity and binding selectivity, probably through complex formation in the RNA major groove. Some intercalators and a group of diphenylfuran cations have strong interactions with RNA that are very dependent on compound structure. RNA hairpin model systems for the RRE binding site of HIV-1 Rev protein were constructed for more detailed investigations. The diphenylfuran cations bind strongly to RRE and selectively inhibit Rev binding. CD, NMR, and fluorescence binding studies indicate that the active compounds bind in the internal loop region of RRE (with binding constants >10⁷M⁻¹), and cause a conformational change in the RNA. None of the standard nucleic acid binding modes appears to fit the results for complexes of the active compounds with RRE, and it is proposed that the diphenylfuran system threads through the internal loop region of RRE. Such a model allows contacts of the furan cationic substituents with both grooves of RRE in addition to the intercalation interactions with the bases.     

Influence of RNA Binding on the Structure and Dynamics of the Lassa Virus Nucleoprotein

Lassa virus protects its viral genome through the formation of a ribonucleoprotein complex in which the nucleoprotein (NP) encapsidates the single-stranded RNA genome. Crystal structures provide evidence that a conformational change must occur to allow for RNA binding. In this study, the mechanism by which NP binds to RNA and how the conformational changes in NP are achieved was investigated with molecular-dynamics simulations. NP was structurally characterized in an open configuration when bound to RNA and in a closed form in the absence of RNA. Our results show that when NP is bound to RNA, the protein is highly dynamic and the system undergoes spontaneous deviations away from the open-state configuration. The equilibrium simulations are supported by free-energy calculations that quantify the influence of RNA on the free-energy surface, which governs NP dynamics. We predict that the globally stable states are qualitatively in agreement with the observed crystal structures, but that both open and closed conformations are thermally accessible in the presence of RNA. The free-energy calculations also provide a prediction of the location of the transition state for RNA binding and identify an intermediate metastable state that exhibits correlated motions that could promote RNA binding.     

RBPBind: Quantitative Prediction of Protein-RNA Interactions

There are hundreds of RNA binding proteins in the human genome alone and their interactions with messenger and other RNAs in a cell regulate every step in an RNA's life cycle. To understand this interplay of proteins and RNA it is important to be able to know which protein binds which RNA how strongly and where. Here, we introduce RBPBind, a web-based tool for the quantitative prediction of the interaction of single-stranded RNA binding proteins with target RNAs that fully takes into account the effect of RNA secondary structure on binding affinity. Given a user-specified RNA and a protein selected from a set of several RNA-binding proteins, RBPBind computes their binding curve and effective binding constant. The server also computes the probability that, at a given protein concentration, a protein molecule will bind to any particular nucleotide along the RNA. The sequence specificity of the protein-RNA interaction is parameterized from public RNAcompete experiments and integrated into the recursions of the Vienna RNA package to simultaneously take into account protein binding and RNA secondary structure. We validate our approach by comparison to experimentally determined binding affinities of the HuR protein for several RNAs of different sequence contexts from the literature, showing that integration of raw sequence affinities into RNA secondary structure prediction significantly improves the agreement between computationally predicted and experimentally measured binding affinities. Our resource thus provides a quick and easy way to obtain reliable predicted binding affinities and locations for single-stranded RNA binding proteins based on RNA sequence alone.     

Structural features of adenovirus 2 virus-associated RNA required for binding to the protein kinase DAI.

The double-stranded RNA activated protein kinase DAI contains an RNA binding domain consisting of two copies of a double-stranded RNA binding motif. We have investigated the role of RNA structure in the interaction between DAI and the structured single-stranded RNA, adenovirus VA RNAI, which inhibits DAI activation. Mutations in the apical stem, terminal stem, and central domain of the RNA were tested to assess the contribution of these elements to DAI binding in vitro. The data demonstrate that over half a turn of intact apical stem is required for the interaction and that there is a correlation between the binding of apical stem mutants and their ability to function both in vivo and in vitro. There was also evidence of preference for GC-rich sequence in the proximal region of the apical stem. In the central domain the correlation between binding and function of mutant RNAs was poor, suggesting that at least some of this region plays no direct role in binding to DAI, despite its functional importance. Exceptionally, central domain mutations that encroached on the phylogenetically conserved stem 4 of VA RNA disrupted binding, and complementary mutations in this sequence partially restored binding. Measurement of the binding of wild-type VA RNAI to DAI and p20, a truncated form of the protein containing the RNA binding domains alone, under various ionic conditions imply that the major interactions are electrostatic and occur via the protein's RNA binding domain. However, differences between full-length DAI and p20 in their binding to mutants in the conserved stem suggest that regions outside the RNA binding domain also participate in the binding. The additional interactions are likely to be non-ionic, and may be important for preventing DAI activation during virus infection.     

Substitutions of Thr30 provide mechanistic insight into tryptophan-mediated activation of TRAP binding to RNA

TRAP is an 11 subunit RNA binding protein that regulates expression of genes involved in tryptophan biosynthesis and transport in Bacillus subtilis. TRAP is activated to bind RNA by binding up to 11 molecules of l-tryptophan in pockets formed by adjacent subunits. The precise mechanism by which tryptophan binding activates TRAP is not known. Thr30 is in the tryptophan binding pocket. A TRAP mutant in which Thr30 is substituted with Val (T30V) does not bind tryptophan but binds RNA constitutively, suggesting that Thr30 plays a key role in the activation mechanism. We have examined the effects of other substitutions of Thr30. TRAP proteins with small beta-branched aliphatic side chains at residue 30 bind RNA constitutively, whereas those with a small polar side chain show tryptophan-dependent RNA binding. Several mutant proteins exhibited constitutive RNA binding that was enhanced by tryptophan. Although the tryptophan and RNA binding sites on TRAP are distinct and are separated by approximately 7.5 A, several substitutions of residues that interact with the bound RNA restored tryptophan binding to T30V TRAP. These observations support the hypothesis that conformational changes in TRAP relay information between the tryptophan and RNA binding sites of the protein.     

Specificity of RNA Binding by CPEB: Requirement for RNA Recognition Motifs and a Novel Zinc Finger

CPEB is an RNA binding protein that interacts with the maturation-type cytoplasmic polyadenylation element (CPE) (consensus UUUUUAU) to promote polyadenylation and translational activation of maternal mRNAs in Xenopus laevis. CPEB, which is conserved from mammals to invertebrates, is composed of three regions: an amino-terminal portion with no obvious functional motif, two RNA recognition motifs (RRMs), and a cysteine-histidine region that is reminiscent of a zinc finger. In this study, we investigated the physical properties of CPEB required for RNA binding. CPEB can interact with RNA as a monomer, and phosphorylation, which modifies the protein during oocyte maturation, has little effect on RNA binding. Deletion mutations of CPEB have been overexpressed in Escherichia coli and used in a series of RNA gel shift experiments. Although a full-length and a truncated CPEB that lacks 139 amino-terminal amino acids bind CPE-containing RNA avidly, proteins that have had either RRM deleted bind RNA much less efficiently. CPEB that has had the cysteine-histidine region deleted has no detectable capacity to bind RNA. Single alanine substitutions of specific cysteine or histidine residues within this region also abolish RNA binding, pointing to the importance of this highly conserved domain of the protein. Chelation of metal ions by 1,10-phenanthroline inhibits the ability of CPEB to bind RNA; however, RNA binding is restored if the reaction is supplemented with zinc. CPEB also binds other metals such as cobalt and cadmium, but these destroy RNA binding. These data indicate that the RRMs and a zinc finger region of CPEB are essential for RNA binding.     

Utilization of a mammalian cell-based RNA binding assay to characterize the RNA binding properties of picornavirus 3C proteinases.

Using an assay capable of detecting sequence-specific RNA/protein interactions in mammalian cells, we demonstrate that the poliovirus and rhinovirus 3C proteinases are able to bind structured target RNA sequences derived from their respective 5' noncoding regions in vivo. Specific RNA binding by poliovirus 3C was found to be dependent on the integrity of stem-loop d of the RNA cloverleaf structure located at the 5' end of poliovirus genomic RNA. In contrast, mutation of stem-loop b did not prevent this in vivo interaction. However, mutation of stem-loop b, which serves as the RNA binding site for a cellular co-factor important for efficient poliovirus replication, did significantly attenuate the efficiency of 3C RNA binding in vivo and 3CD RNA binding in vitro. This in vivo protein:RNA binding assay was also used to identify several residues in 3C that are critical for RNA binding, but dispensable for 3C proteinase activity. The mammalian cell-based RNA binding assay described in this study may have considerable potential utility in the future detection or analysis of in vivo RNA/protein interactions unrelated to the 3C/RNA interaction described here.