Recent advances in RNA-protein recognition

Nucleolar stress induces nucleolar stress body formation via the NOSR-1/NUMR-1 axis in Caenorhabditis elegans

Environmental stimuli not only alter gene expression profiles but also induce structural changes in cells. How distinct nuclear bodies respond to cellular stress is poorly understood. Here, we identify a subnuclear organelle named the nucleolar stress body (NoSB), the formation of which is induced by the inhibition of rRNA transcription or inactivation of rRNA processing and maturation in C. elegans. NoSB does not colocalize with other previously described subnuclear organelles. We conduct forward genetic screening and identify a bZIP transcription factor, named nucleolar stress response-1 (NOSR-1), that is required for NoSB formation. The inhibition of rRNA transcription or inactivation of rRNA processing and maturation increases nosr-1 expression. By using transcriptome analysis of wild-type animals subjected to different nucleolar stress conditions and nosr-1 mutants, we identify that the SR-like protein NUMR-1 (nuclear localized metal responsive) is the target of NOSR-1. Interestingly, NUMR-1 is a component of NoSB and itself per se is required for the formation of NoSB. We conclude that the NOSR-1/NUMR-1 axis likely responds to nucleolar stress and mediates downstream stress-responsive transcription programs and subnuclear morphology alterations in C. elegans.

RNA binding proteins: Diversity from microsurgeons to cowboys

The conformation and mechanism of proteins that degrade and bind RNA, which has provided key insights into post-transcriptional gene regulation, is explored here. During the twentieth century's last decades, the characterization of ribonucleases and RNA binding domains revealed the diversity of their reaction mechanisms and modes of RNA recognition, and the bases of protein folding, substrate specificity and binding affinity. More recent research showed how these domains combine through oligomerization or genetic recombination to create larger proteins with highly specific and readily programmable ribonucleolytic activity. In the last 15 years, the study of the capacity of proteins, usually disordered, to pool RNAs into discrete, non-aqueous microdroplets to facilitate their transport, modification and degradation - analogous to cowboys herding cattle - has advanced our comprehension of gene expression. Finally, the current uses of RNA binding proteins and the future applications of protein/RNA microdroplets are highlighted.

Coupling between conformational dynamics and catalytic function at the active site of the lead-dependent ribozyme

In common with other self-cleaving RNAs, the lead-dependent ribozyme (leadzyme) undergoes dynamic fluctuations to a chemically activated conformation. We explored the connection between conformational dynamics and self-cleavage function in the leadzyme using a combination of NMR spin-relaxation analysis of ribose groups and conformational restriction via chemical modification. The functional studies were performed with a North-methanocarbacytidine modification that prevents fluctuations to C2'-endo conformations while maintaining an intact 2'-hydroxyl nucleophile. Spin-relaxation data demonstrate that the active-site Cyt-6 undergoes conformational exchange attributed to sampling of a minor C2'-endo state with an exchange lifetime on the order of microseconds to tens of microseconds. A conformationally restricted species in which the fluctuations to the minor species are interrupted shows a drastic decrease in self-cleavage activity. Taken together, these data indicate that dynamic sampling of a minor species at the active site of this ribozyme, and likely of related naturally occurring motifs, is strongly coupled to catalytic function. The combination of NMR dynamics analysis with functional probing via conformational restriction is a general methodology for dissecting dynamics-function relationships in RNA.

The structure of transcription termination factor Nrd1 reveals an original mode for GUAA recognition

Transcription termination of non-coding RNAs is regulated in yeast by a complex of three RNA binding proteins: Nrd1, Nab3 and Sen1. Nrd1 is central in this process by interacting with Rbp1 of RNA polymerase II, Trf4 of TRAMP and GUAA/G terminator sequences. We lack structural data for the last of these binding events. We determined the structures of Nrd1 RNA binding domain and its complexes with three GUAA-containing RNAs, characterized RNA binding energetics and tested rationally designed mutants in vivo. The Nrd1 structure shows an RRM domain fused with a second α/β domain that we name split domain (SD), because it is formed by two non-consecutive segments at each side of the RRM. The GUAA interacts with both domains and with a pocket of water molecules, trapped between the two stacking adenines and the SD. Comprehensive binding studies demonstrate for the first time that Nrd1 has a slight preference for GUAA over GUAG and genetic and functional studies suggest that Nrd1 RNA binding domain might play further roles in non-coding RNAs transcription termination.

Non-canonical binding interactions of the RNA recognition motif (RRM) domains of P34 protein modulate binding within the 5S ribonucleoprotein particle (5S RNP)

RNA binding proteins are involved in many aspects of RNA metabolism. In Trypanosoma brucei, our laboratory has identified two trypanosome-specific RNA binding proteins P34 and P37 that are involved in the maturation of the 60S subunit during ribosome biogenesis. These proteins are part of the T. brucei 5S ribonucleoprotein particle (5S RNP) and P34 binds to 5S ribosomal RNA (rRNA) and ribosomal protein L5 through its N-terminus and its RNA recognition motif (RRM) domains. We generated truncated P34 proteins to determine these domains’ interactions with 5S rRNA and L5. Our analyses demonstrate that RRM1 of P34 mediates the majority of binding with 5S rRNA and the N-terminus together with RRM1 contribute the most to binding with L5. We determined that the consensus ribonucleoprotein (RNP) 1 and 2 sequences, characteristic of canonical RRM domains, are not fully conserved in the RRM domains of P34. However, the aromatic amino acids previously described to mediate base stacking interactions with their RNA target are conserved in both of the RRM domains of P34. Surprisingly, mutation of these aromatic residues did not disrupt but instead enhanced 5S rRNA binding. However, we identified four arginine residues located in RRM1 of P34 that strongly impact L5 binding. These mutational analyses of P34 suggest that the binding site for 5S rRNA and L5 are near each other and specific residues within P34 regulate the formation of the 5S RNP. These studies show the unique way that the domains of P34 mediate binding with the T. brucei 5S RNP.

Solution structure of the first RNA recognition motif domain of human spliceosomal protein SF3b49 and its mode of interaction with a SF3b145 fragment

The spliceosomal protein SF3b49, a component of the splicing factor 3b (SF3b) protein complex in the U2 small nuclear ribonucleoprotein, contains two RNA recognition motif (RRM) domains. In yeast, the first RRM domain (RRM1) of Hsh49 protein (yeast orthologue of human SF3b49) reportedly interacts with another component, Cus1 protein (orthologue of human SF3b145). Here, we solved the solution structure of the RRM1 of human SF3b49 and examined its mode of interaction with a fragment of human SF3b145 using NMR methods. Chemical shift mapping showed that the SF3b145 fragment spanning residues 598–631 interacts with SF3b49 RRM1, which adopts a canonical RRM fold with a topology of β1‐α1‐β2‐β3‐α2‐β4. Furthermore, a docking model based on NOESY measurements suggests that residues 607–616 of the SF3b145 fragment adopt a helical structure that binds to RRM1 predominantly via α1, consequently exhibiting a helix–helix interaction in almost antiparallel. This mode of interaction was confirmed by a mutational analysis using GST pull‐down assays. Comparison with structures of all RRM domains when complexed with a peptide found that this helix–helix interaction is unique to SF3b49 RRM1. Additionally, all amino acid residues involved in the interaction are well conserved among eukaryotes, suggesting evolutionary conservation of this interaction mode between SF3b49 RRM1 and SF3b145.

PDB Code(s): 5GVQ

A structural perspective of RNA recognition by intrinsically disordered proteins

Protein-RNA recognition is essential for gene expression and its regulation, which is indispensable for the survival of the living organism at one hand, on the other hand, misregulation of this recognition may lead to their extinction. Polymorphic conformation of both the interacting partners is a characteristic feature of such molecular recognition that promotes the assembly. Many RNA binding proteins (RBP) or regions in them are found to be intrinsically disordered, and this property helps them to play a central role in the regulatory processes. Sequence composition and the length of the flexible linkers between RNA binding domains in RBPs are crucial in making significant contacts with its partner RNA. Polymorphic conformations of RBPs can provide thermodynamic advantage to its binding partner while acting as a chaperone. Prolonged extensions of the disordered regions in RBPs also contribute to the stability of the large cellular machines including ribosome and viral assemblies. The involvement of these disordered regions in most of the significant cellular processes makes RBPs highly associated with various human diseases that arise due to their misregulation.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

A structural perspective of RNA recognition by intrinsically disordered proteins

Protein-RNA recognition is essential for gene expression and its regulation, which is indispensable for the survival of the living organism at one hand, on the other hand, misregulation of this recognition may lead to their extinction. Polymorphic conformation of both the interacting partners is a characteristic feature of such molecular recognition that promotes the assembly. Many RNA binding proteins (RBP) or regions in them are found to be intrinsically disordered, and this property helps them to play a central role in the regulatory processes. Sequence composition and the length of the flexible linkers between RNA binding domains in RBPs are crucial in making significant contacts with its partner RNA. Polymorphic conformations of RBPs can provide thermodynamic advantage to its binding partner while acting as a chaperone. Prolonged extensions of the disordered regions in RBPs also contribute to the stability of the large cellular machines including ribosome and viral assemblies. The involvement of these disordered regions in most of the significant cellular processes makes RBPs highly associated with various human diseases that arise due to their misregulation.

Attenuating HIV Tat/TAR-mediated protein expression by exploring the side chain length of positively charged residues

RNA is a drug target involved in diverse cellular functions and viral processes. Molecules that inhibit the HIV TAR RNA-Tat protein interaction may attenuate Tat/TAR-dependent protein expression and potentially serve as anti-HIV therapeutics. By incorporating positively charged residues with mixed side chain lengths, we designed peptides that bind TAR RNA with enhanced intracellular activity. Tat-derived peptides that were individually substituted with positively charged residues with varying side chain lengths were evaluated for TAR RNA binding. Positively charged residues with different side chain lengths were incorporated at each Arg and Lys position in the Tat-derived peptide to enhance TAR RNA binding. The resulting peptides showed enhanced TAR RNA binding affinity, cellular uptake, nuclear localization, proteolytic resistance, and inhibition of intracellular Tat/TAR-dependent protein expression compared to the parent Tat-derived peptide with no cytotoxicity. Apparently, the enhanced inhibition of protein expression by these peptides was not determined by RNA binding affinity, but by proteolytic resistance. Despite the high TAR binding affinity, a higher binding specificity would be necessary for practical purposes. Importantly, altering the positively charged residue side chain length should be a viable strategy to generate potentially useful RNA-targeting bioactive molecules.

Regulation of alternative splicing in Drosophila by 56 RNA binding proteins

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.

Computational challenges, tools, and resources for analyzing co- and post-transcriptional events in high throughput

Co- and post-transcriptional regulation of gene expression is complex and multifaceted, spanning the complete RNA lifecycle from genesis to decay. High-throughput profiling of the constituent events and processes is achieved through a range of technologies that continue to expand and evolve. Fully leveraging the resulting data is nontrivial, and requires the use of computational methods and tools carefully crafted for specific data sources and often intended to probe particular biological processes. Drawing upon databases of information pre-compiled by other researchers can further elevate analyses. Within this review, we describe the major co- and post-transcriptional events in the RNA lifecycle that are amenable to high-throughput profiling. We place specific emphasis on the analysis of the resulting data, in particular the computational tools and resources available, as well as looking toward future challenges that remain to be addressed.

Structural and biochemical analysis of the Hordeum vulgare L. HvGR-RBP1 protein, a glycine-rich RNA-binding protein involved in the regulation of barley plant development and stress response

The timing of whole-plant senescence influences important agricultural traits such as yield and grain protein content. Post-transcriptional regulation by plant RNA-binding proteins is essential for proper control of gene expression, development, and stress responses. Here, we report the three-dimensional solution NMR structure and nucleic acid-binding properties of the barley glycine-rich RNA-binding protein HvGR-RBP1, whose transcript has been identified as being >45-fold up-regulated in early-as compared to late-senescing near-isogenic barley germplasm. NMR analysis reveals that HvGR-RBP1 is a multidomain protein comprising a well-folded N-terminal RNA Recognition Motif (RRM) and a structurally disordered C-terminal glycine-rich domain. Chemical shift differences observed in 2D (1)H-(15)N correlation (HSQC) NMR spectra of full-length HvGR-RBP1 and N-HvGR-RBP1 (RRM domain only) suggest that the two domains can interact both in-trans and intramolecularly, similar to what is observed in the tobacco NtGR-RBP1 protein. Further, we show that the RRM domain of HvGR-RBP1 binds single-stranded DNA nucleotide fragments containing the consensus nucleotide sequence 5'-TTCTGX-3' with low micromolar affinity in vitro. We also demonstrate that the C-terminal glycine-rich (HvGR) domain of Hv-GR-RBP1 can interact nonspecifically with ssRNA in vitro. Structural similarities with other plant glycine-rich RNA-binding proteins suggest that HvGR-RBP1 may be multifunctional. Based on gene expression analysis following cold stress in barley and E. coli growth studies following cold shock treatment, we conclude that HvGR-RBP1 functions in a manner similar to cold-shock proteins and harbors RNA chaperone activity. HvGR-RBP1 is therefore not only involved in the regulation of barley development including senescence, but also functions in plant responses to environmental stress.

Computational Methods for CLIP-seq Data Processing

RNA-binding proteins (RBPs) are at the core of post-transcriptional regulation and thus of gene expression control at the RNA level. One of the principal challenges in the field of gene expression regulation is to understand RBPs mechanism of action. As a result of recent evolution of experimental techniques, it is now possible to obtain the RNA regions recognized by RBPs on a transcriptome-wide scale. In fact, CLIP-seq protocols use the joint action of CLIP, crosslinking immunoprecipitation, and high-throughput sequencing to recover the transcriptome-wide set of interaction regions for a particular protein. Nevertheless, computational methods are necessary to process CLIP-seq experimental data and are a key to advancement in the understanding of gene regulatory mechanisms. Considering the importance of computational methods in this area, we present a review of the current status of computational approaches used and proposed for CLIP-seq data.

The mysterious RAMP proteins and their roles in small RNA-based immunity

A new class of prokaryotic RNA binding proteins called Repeat Associated Mysterious Proteins (RAMPs), has recently been identified. These proteins play key roles in a novel type immunity in which the DNA of the host organism (e.g. a prokaryote) has sequence segments corresponding to the sequences of potential viral invaders. The sequences embedded in the host DNA confer immunity by directing selective destruction of the nucleic acid of the virus using an RNA-based strategy. In this viral defense mechanism, RAMP proteins have multiple functional roles including endoribonucleotic cleavage and ribonucleoprotein particle assembly. RAMPs contain the classical RNA recognition motif (RRM), often in tandem, and a conserved glycine-rich segment (G-loop) near the carboxyl terminus. However, unlike RRMs that bind single-stranded RNA using their β-sheet surface, RAMPs make use of both sides of the RRM fold and interact with both single-stranded and structured RNA. The unique spatial arrangement of the two RRM folds, facilitated by a hallmark G-loop, is crucial to formation of a composite surface for recognition of specific RNA. Evidence for RNA-dependent oligomerization is also observed in some RAMP proteins that may serve as an important strategy to increase specificity.

Structural insight into RNA recognition motifs: versatile molecular Lego building blocks for biological systems

'RNA recognition motifs (RRMs)' are common domain-folds composed of 80-90 amino-acid residues in eukaryotes, and have been identified in many cellular proteins. At first they were known as RNA binding domains. Through discoveries over the past 20 years, however, the RRMs have been shown to exhibit versatile molecular recognition activities and to behave as molecular Lego building blocks to construct biological systems. Novel RNA/protein recognition modes by RRMs are being identified, and more information about the molecular recognition by RRMs is becoming available. These RNA/protein recognition modes are strongly correlated with their biological significance. In this review, we would like to survey the recent progress on these versatile molecular recognition modules.

Pub1p C-terminal RRM domain interacts with Tif4631p through a conserved region neighbouring the Pab1p binding site

Pub1p, a highly abundant poly(A)+ mRNA binding protein in Saccharomyces cerevisiae, influences the stability and translational control of many cellular transcripts, particularly under some types of environmental stresses. We have studied the structure, RNA and protein recognition modes of different Pub1p constructs by NMR spectroscopy. The structure of the C-terminal RRM domain (RRM3) shows a non-canonical N-terminal helix that packs against the canonical RRM fold in an original fashion. This structural trait is conserved in Pub1p metazoan homologues, the TIA-1 family, defining a new class of RRM-type domains that we propose to name TRRM (TIA-1 C-terminal domain-like RRM). Pub1p TRRM and the N-terminal RRM1-RRM2 tandem bind RNA with high selectivity for U-rich sequences, with TRRM showing additional preference for UA-rich ones. RNA-mediated chemical shift changes map to β-sheet and protein loops in the three RRMs. Additionally, NMR titration and biochemical in vitro cross-linking experiments determined that Pub1p TRRM interacts specifically with the N-terminal region (1-402) of yeast eIF4G1 (Tif4631p), very likely through the conserved Box1, a short sequence motif neighbouring the Pab1p binding site in Tif4631p. The interaction involves conserved residues of Pub1p TRRM, which define a protein interface that mirrors the Pab1p-Tif4631p binding mode. Neither protein nor RNA recognition involves the novel N-terminal helix, whose functional role remains unclear. By integrating these new results with the current knowledge about Pub1p, we proposed different mechanisms of Pub1p recruitment to the mRNPs and Pub1p-mediated mRNA stabilization in which the Pub1p/Tif4631p interaction would play an important role.

The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA--protein interaction sites in 211 RNA--protein chain pairs, the largest set of RNA--protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA--protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein--RNA interactions, indicating their mutual importance in protein recognition.

Association of OLE RNA with bacterial membranes via an RNA-protein interaction

Ornate, large, extremophilic (OLE) RNAs are large, non-coding transcripts characterized by their ornate secondary structure and presence predominantly in Gram-positive, extremophilic bacteria. A gene for an OLE-associated protein (OAP) is almost always located immediately downstream of the OLE gene. OAP has no extensive homology to other proteins and is predicted to form multiple transmembrane domains. We show that this protein forms a ribonucleoprotein complex with OLE RNA using at least 2:1 protein : RNA stoichiometry. A series of truncated OLE RNA constructs was used to establish that most of the RNA can be deleted without eliminating protein binding. Two primary binding sites are present within the RNA, although additional binding determinants exist and extensive structural stabilization is induced by OAP. RNA fluorescence in situ hybridization (FISH) was used in Escherichia coli to demonstrate that ribonucleoprotein complex formation localizes the RNA near cell membranes of this heterologous system. Therefore, the majority of the complex structure formed by OLE RNA may perform a biochemical function that requires membrane localization.

Structural basis for the dual RNA-recognition modes of human Tra2-β RRM

Human Transformer2-β (hTra2-β) is an important member of the serine/arginine-rich protein family, and contains one RNA recognition motif (RRM). It controls the alternative splicing of several pre-mRNAs, including those of the calcitonin/calcitonin gene-related peptide (CGRP), the survival motor neuron 1 (SMN1) protein and the tau protein. Accordingly, the RRM of hTra2-β specifically binds to two types of RNA sequences [the CAA and (GAA)(2) sequences]. We determined the solution structure of the hTra2-β RRM (spanning residues Asn110-Thr201), which not only has a canonical RRM fold, but also an unusual alignment of the aromatic amino acids on the β-sheet surface. We then solved the complex structure of the hTra2-β RRM with the (GAA)(2) sequence, and found that the AGAA tetra-nucleotide was specifically recognized through hydrogen-bond formation with several amino acids on the N- and C-terminal extensions, as well as stacking interactions mediated by the unusually aligned aromatic rings on the β-sheet surface. Further NMR experiments revealed that the hTra2-β RRM recognizes the CAA sequence when it is integrated in the stem-loop structure. This study indicates that the hTra2-β RRM recognizes two types of RNA sequences in different RNA binding modes.

Binding of the bacteriophage P22 N-peptide to the boxB RNA motif studied by molecular dynamics simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-alpha-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.

Prediction of functional class of proteins and peptides irrespective of sequence homology by support vector machines

Various computational methods have been used for the prediction of protein and peptide function based on their sequences. A particular challenge is to derive functional properties from sequences that show low or no homology to proteins of known function. Recently, a machine learning method, support vector machines (SVM), have been explored for predicting functional class of proteins and peptides from amino acid sequence derived properties independent of sequence similarity, which have shown promising potential for a wide spectrum of protein and peptide classes including some of the low- and non-homologous proteins. This method can thus be explored as a potential tool to complement alignment-based, clustering-based, and structure-based methods for predicting protein function. This article reviews the strategies, current progresses, and underlying difficulties in using SVM for predicting the functional class of proteins. The relevant software and web-servers are described. The reported prediction performances in the application of these methods are also presented.

Mrd1p is required for release of base-paired U3 snoRNA within the preribosomal complex

In eukaryotes, ribosomes are made from precursor rRNA (pre-rRNA) and ribosomal proteins in a maturation process that requires a large number of snoRNPs and processing factors. A fundamental problem is how the coordinated and productive folding of the pre-rRNA and assembly of successive pre-rRNA-protein complexes is achieved cotranscriptionally. The conserved protein Mrd1p, which contains five RNA binding domains (RBDs), is essential for processing events leading to small ribosomal subunit synthesis. We show that full function of Mrd1p requires all five RBDs and that the RBDs are functionally distinct and needed during different steps in processing. Mrd1p mutations trap U3 snoRNA in pre-rRNP complexes both in base-paired and non-base-paired interactions. A single essential RBD, RBD5, is involved in both types of interactions, but its conserved RNP1 motif is not needed for releasing the base-paired interactions. RBD5 is also required for the late pre-rRNP compaction preceding A(2) cleavage. Our results suggest that Mrd1p modulates successive conformational rearrangements within the pre-rRNP that influence snoRNA-pre-rRNA contacts and couple U3 snoRNA-pre-rRNA remodeling and late steps in pre-rRNP compaction that are essential for cleavage at A(0) to A(2). Mrd1p therefore coordinates key events in biosynthesis of small ribosome subunits.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

Structural basis for the sequence-specific RNA-recognition mechanism of human CUG-BP1 RRM3

The CUG-binding protein 1 (CUG-BP1) is a member of the CUG-BP1 and ETR-like factors (CELF) family or the Bruno-like family and is involved in the control of splicing, translation and mRNA degradation. Several target RNA sequences of CUG-BP1 have been predicted, such as the CUG triplet repeat, the GU-rich sequences and the AU-rich element of nuclear pre-mRNAs and/or cytoplasmic mRNA. CUG-BP1 has three RNA-recognition motifs (RRMs), among which the third RRM (RRM3) can bind to the target RNAs on its own. In this study, we solved the solution structure of the CUG-BP1 RRM3 by hetero-nuclear NMR spectroscopy. The CUG-BP1 RRM3 exhibited a noncanonical RRM fold, with the four-stranded β-sheet surface tightly associated with the N-terminal extension. Furthermore, we determined the solution structure of the CUG-BP1 RRM3 in the complex with (UG)₃ RNA, and discovered that the UGU trinucleotide is specifically recognized through extensive stacking interactions and hydrogen bonds within the pocket formed by the β-sheet surface and the N-terminal extension. This study revealed the unique mechanism that enables the CUG-BP1 RRM3 to discriminate the short RNA segment from other sequences, thus providing the molecular basis for the comprehension of the role of the RRM3s in the CELF/Bruno-like family.

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

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.

Structures, kinetics, thermodynamics, and biological functions of RNA hairpins

Most RNA comprises one strand and therefore can fold back on itself to form complex structures. At the heart of these structures is the hairpin, which is composed of a stem having Watson-Crick base pairing and a loop wherein the backbone changes directionality. First, we review the structure of hairpins including diversity in the stem, loop, and closing base pair. The function of RNA hairpins in biology is discussed next, including roles for isolated hairpins, as well as hairpins in the context of complex tertiary structures. We describe the kinetics and thermodynamics of hairpin folding including models for hairpin folding, folding transition states, and the cooperativity of folding. Lastly, we discuss some ways in which hairpins can influence the folding and function of tertiary structures, both directly and indirectly. RNA hairpins provide a simple means of controlling gene expression that can be understood in the language of physical chemistry.

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R(1) and R(2)), and relaxation dispersion (R(1rho)) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (approximately 30 mus) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.

Overview of protein structural and functional folds

This overview provides an illustrated, comprehensive survey of some commonly observed protein‐fold families and structural motifs, chosen for their functional significance. It opens with descriptions and definitions of the various elements of protein structure and associated terminology. Following is an introduction into web‐based structural bioinformatics that includes surveys of interactive web servers for protein fold or domain annotation, protein‐structure databases, protein‐structure‐classification databases, structural alignments of proteins, and molecular graphics programs available for personal computers. The rest of the overview describes selected families of protein folds in terms of their secondary, tertiary, and quaternary structural arrangements, including ribbon‐diagram examples, tables of representative structures with references, and brief explanations pointing out their respective biological and functional significance.

RNA structure determination by NMR

This chapter reviews the methodologies for RNA structure determination by liquid-state nuclear magnetic resonance (NMR). The routine production of milligram quantities of isotopically labeled RNA remains critical to the success of NMR-based structure studies. The standard method for the preparation of isotopically labeled RNA for structural studies in solution is in vitro transcription from DNA oligonucleotide templates using T7 RNA polymerase and unlabeled or isotopically labeled nucleotide triphosphates (NTPs). The purification of the desired RNA can be performed by either denaturing polyacrylamide gel electrophoresis (PAGE) or anion-exchange chromatography. Our basic strategy for studying RNA in solution by NMR is outlined. The topics covered include RNA resonance assignment, restraint collection, and the structure calculation process. Selected examples of NMR spectra are given for a correctly folded 30 nucleotide-containing RNA.

A knowledge-based potential function predicts the specificity and relative binding energy of RNA-binding proteins

RNA-protein interactions are fundamental to gene expression. Thus, the molecular basis for the sequence dependence of protein-RNA recognition has been extensively studied experimentally. However, there have been very few computational studies of this problem, and no sustained attempt has been made towards using computational methods to predict or alter the sequence-specificity of these proteins. In the present study, we provide a distance-dependent statistical potential function derived from our previous work on protein-DNA interactions. This potential function discriminates native structures from decoys, successfully predicts the native sequences recognized by sequence-specific RNA-binding proteins, and recapitulates experimentally determined relative changes in binding energy due to mutations of individual amino acids at protein-RNA interfaces. Thus, this work demonstrates that statistical models allow the quantitative analysis of protein-RNA recognition based on their structure and can be applied to modeling protein-RNA interfaces for prediction and design purposes.

Efficacy of different protein descriptors in predicting protein functional families

Background

Sequence-derived structural and physicochemical descriptors have frequently been used in machine learning prediction of protein functional families, thus there is a need to comparatively evaluate the effectiveness of these descriptor-sets by using the same method and parameter optimization algorithm, and to examine whether the combined use of these descriptor-sets help to improve predictive performance. Six individual descriptor-sets and four combination-sets were evaluated in support vector machines (SVM) prediction of six protein functional families.

Results

The performance of these descriptor-sets were ranked by Matthews correlation coefficient (MCC), and categorized into two groups based on their performance. While there is no overwhelmingly favourable choice of descriptor-sets, certain trends were found. The combination-sets tend to give slightly but consistently higher MCC values and thus overall best performance such that three out of four combination-sets show slightly better performance compared to one out of six individual descriptor-sets.

Conclusion

Our study suggests that currently used descriptor-sets are generally useful for classifying proteins and the prediction performance may be enhanced by exploring combinations of descriptors.

Muscleblind-like 1 interacts with RNA hairpins in splicing target and pathogenic RNAs

The MBNL and CELF proteins act antagonistically to control the alternative splicing of specific exons during mammalian postnatal development. This process is dysregulated in myotonic dystrophy because MBNL proteins are sequestered by (CUG)_n and (CCUG)_n RNAs expressed from mutant DMPK and ZNF9 genes, respectively. While these observations predict that MBNL proteins have a higher affinity for these pathogenic RNAs versus their normal splicing targets, we demonstrate that MBNL1 possesses comparably high affinities for (CUG)_n and (CAG)_n RNAs as well as a splicing target, Tnnt3. Mapping of a MBNL1-binding site upstream of the Tnnt3 fetal exon indicates that a preferred binding site for this protein is a GC-rich RNA hairpin containing a pyrimidine mismatch. To investigate how pathogenic RNAs sequester MBNL1 in DM1 cells, we used a combination of chemical/enzymatic structure probing and electron microscopy to determine that MBNL1 forms a ring-like structure which binds to the dsCUG helix. While the MBNL1 N-terminal region is required for RNA binding, the C-terminal region mediates homotypic interactions which may stabilize intra- and/or inter-ring interactions. Our results provide a mechanistic basis for dsCUG-induced MBNL1 sequestration and highlight a striking similarity in the binding sites for MBNL proteins on splicing precursor and pathogenic RNAs.

Trypanosoma cruzi: molecular characterization of an RNA binding protein differentially expressed in the parasite life cycle

Molecular studies have shown several peculiarities in the regulatory mechanisms of gene expression in trypanosomatids. Protein coding genes are organized in long polycistronic units that seem to be constitutively transcribed. Therefore, post-transcriptional regulation of gene expression is considered to be the main point for control of transcript abundance and functionality. Here we describe the characterization of a 17 kDa RNA-binding protein from Trypanosoma cruzi (TcRBP19) containing an RNA recognition motive (RRM). This protein is coded by a single copy gene located in a high molecular weight chromosome of T. cruzi. Orthologous genes are present in the TriTryp genomes. TcRBP19 shows target selectivity since among the different homoribopolymers it preferentially binds polyC. TcRBP19 is a low expression protein only barely detected at the amastigote stage localizing in a diffuse pattern in the cytoplasm.

Pseudouridine synthases

Pseudouridine synthases are the enzymes responsible for the most abundant posttranscriptional modification of cellular RNAs. These enzymes catalyze the site-specific isomerization of uridine residues that are already part of an RNA chain, and appear to employ both sequence and structural information to achieve site specificity. Crystallographic analyses have demonstrated that all pseudouridine synthases share a common core fold and active site structure and that this core is modified by peripheral domains, accessory proteins, and guide RNAs to give rise to remarkable substrate versatility.

Prediction of the functional class of metal-binding proteins from sequence derived physicochemical properties by support vector machine approach

Metal-binding proteins play important roles in structural stability, signaling, regulation, transport, immune response, metabolism control, and metal homeostasis. Because of their functional and sequence diversity, it is desirable to explore additional methods for predicting metal-binding proteins irrespective of sequence similarity. This work explores support vector machines (SVM) as such a method. SVM prediction systems were developed by using 53,333 metal-binding and 147,347 non-metal-binding proteins, and evaluated by an independent set of 31,448 metal-binding and 79,051 non-metal-binding proteins. The computed prediction accuracy is 86.3%, 81.6%, 83.5%, 94.0%, 81.2%, 85.4%, 77.6%, 90.4%, 90.9%, 74.9% and 78.1% for calcium-binding, cobalt-binding, copper-binding, iron-binding, magnesium-binding, manganese-binding, nickel-binding, potassium-binding, sodium-binding, zinc-binding, and all metal-binding proteins respectively. The accuracy for the non-member proteins of each class is 88.2%, 99.9%, 98.1%, 91.4%, 87.9%, 94.5%, 99.2%, 99.9%, 99.9%, 98.0%, and 88.0% respectively. Comparable accuracies were obtained by using a different SVM kernel function. Our method predicts 67% of the 87 metal-binding proteins non-homologous to any protein in the Swissprot database and 85.3% of the 333 proteins of known metal-binding domains as metal-binding. These suggest the usefulness of SVM for facilitating the prediction of metal-binding proteins. Our software can be accessed at the SVMProt server http://jing.cz3.nus.edu.sg/cgi-bin/svmprot.cgi.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

Hijacking of the host-cell response and translational control during influenza virus infection

Influenza virus is a major public health problem with annual deaths in the US of 36,000 with pandemic outbreaks, such as in 1918, resulting in deaths exceeding 20 million worldwide. Recently, there is much concern over the introduction of highly pathogenic avian influenza H5N1 viruses into the human population. Influenza virus has evolved complex translational control strategies that utilize cap-dependent translation initiation mechanisms and involve the recruitment of both viral and host-cell proteins to preferentially synthesize viral proteins and prevent activation of antiviral responses. Influenza virus is a member of the Orthomyxoviridae family of negative-stranded, segmented RNA viruses and represents a particularly attractive model system as viral replication strategies are closely intertwined with normal cellular processes including the host defense and stress pathways. In this chapter, we review the parallels between translational control in influenza virus infected cells and in stressed cells with a focus on selective translation of viral mRNAs and the antagonism of the dsRNA and host antiviral responses. Moreover, we will discuss how the use of genomic technologies such as DNA microarrays and high through-put proteomics can be used to gain new insights into the control of protein synthesis during viral infection and provide a near comprehensive view of virus-host interactions.

Grabbing the message: structural basis of mRNA 3'UTR recognition by Hrp1

The recognition of specific signals encoded within the 3'-untranslated region of the newly transcribed mRNA triggers the assembly of a multiprotein machine that modifies its 3'-end. Hrp1 recognises one of such signals, the so-called polyadenylation enhancement element (PEE), promoting the recruitment of other polyadenylation factors in yeast. The molecular bases of this interaction are revealed here by the solution structure of a complex between Hrp1 and an oligonucleotide mimicking the PEE. Six consecutive bases (AUAUAU) are specifically recognised by two RNA-binding domains arranged in tandem. Both protein and RNA undergo significant conformational changes upon complex formation with a concomitant large surface burial of RNA bases. Key aspects of RNA specificity can be explained by the presence of intermolecular aromatic-aromatic contacts and hydrogen bonds. Altogether, the Hrp1-PEE structure represents one of the first steps towards understanding of the assembly of the cleavage and polyadenylation machinery at the atomic level.

Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure

We analyzed structural features of 11,038 direct atomic contacts (either electrostatic, H-bonds, hydrophobic, or other van der Waals interactions) extracted from 139 protein-DNA and 49 protein-RNA nonhomologous complexes from the Protein Data Bank (PDB). Globally, H-bonds are the most frequent interactions (approximately 50%), followed by van der Waals, hydrophobic, and electrostatic interactions. From the protein viewpoint, hydrophilic amino acids are over-represented in the interaction databases: Positively charged amino acids mainly contact nucleic acid phosphate groups but can also interact with base edges. From the nucleotide point of view, DNA and RNA behave differently: Most protein-DNA interactions involve phosphate atoms, while protein-RNA interactions involve more frequently base edge and ribose atoms. The increased participation of DNA phosphate involves H-bonds rather than salt bridges. A statistical analysis was performed to find the occurrence of amino acid-nucleotide pairs most different from chance. These pairs were analyzed individually. Finally, we studied the conformation of DNA in the interaction sites. Despite the prevalence of B-DNA in the database, our results suggest that A-DNA is favored in the interaction sites.

Drosophila Sex-lethal protein mediates polyadenylation switching in the female germline

The Drosophila master sex-switch protein Sex-lethal (SXL) regulates the splicing and/or translation of three known targets to mediate somatic sexual differentiation. Genetic studies suggest that additional target(s) of SXL exist, particularly in the female germline. Surprisingly, our detailed molecular characterization of a new potential target of SXL, enhancer of rudimentary (e(r)), reveals that SXL regulates e(r) by a novel mechanism--polyadenylation switching--specifically in the female germline. SXL binds to multiple SXL-binding sites, which include the GU-rich poly(A) enhancer, and competes for the binding of CstF64 in vitro. The SXL-binding sites are able to confer sex-specific poly(A) switching onto an otherwise nonresponsive polyadenylation signal in vivo. The sex-specific poly(A) switching of e(r) provides a means for translational regulation in germ cells. We present a model for the SXL-dependent poly(A) site choice in the female germline.

Mass spectrometry of RNA

A complex population of non-coding RNAs is present in higher organisms. These RNAs have a multitude of functions and execute control over gene expression through various, often poorly understood, mechanisms. At present, the identification and analysis of functional regulatory RNAs and disparate ribonucleoprotein complexes remain an experimental challenge for biologists. They require specially designed approaches and techniques in genomics and RNA biochemistry. Developments in technologies based on mass spectrometry could offer sensitive and efficient solutions to analysis of the sequence, structure, modification and composition of RNA.