Structural basis of single-stranded RNA recognition

Predicting the involvement of polyQ- and polyA in protein-protein interactions by their amino acid context

Homorepeats, specifically polyglutamine (polyQ) and polyalanine (polyA), are often implicated in protein-protein interactions (PPIs). So far, a method to predict the participation of homorepeats in protein interactions is lacking. We propose a machine learning approach to identify PPI-involved polyQ and polyA regions within the human proteome based on known interacting regions. Using the dataset of human homorepeats, we identified 157 polyQ and 745 polyA regions potentially involved in PPIs. Machine learning models, trained on amino acid context and homorepeat length, demonstrated high precision (0.90-0.98) but variable recall (0.42-0.85). Random forest outperformed other models (AUC polyQ = 0.686, AUC polyA = 0.732) using the positions surrounding the homorepeat -10 to +10. Integrating paralog information marginally improved predictions but was excluded for model simplicity. Further optimization revealed that for polyQ, using amino acid surrounding positions from -6 to +6 increased AUC to 0.715. For polyA, no improvement was found. Incorporating coiled coil overlap information enhanced polyA predictions (AUC = 0.745) but not polyQ. Finally, we applied these models to predict PPI involvement across all polyQ and polyA regions, identifying potential interactions. Case studies illustrated the method's predictive capacity, highlighting known interacting regions with high scores and elucidating potential false negatives.

Repurposing the mammalian RNA-binding protein Musashi-1 as an allosteric translation repressor in bacteria

The RNA recognition motif (RRM) is the most common RNA-binding protein domain identified in nature. However, RRM-containing proteins are only prevalent in eukaryotic phyla, in which they play central regulatory roles. Here, we engineered an orthogonal post-transcriptional control system of gene expression in the bacterium Escherichia coli with the mammalian RNA-binding protein Musashi-1, which is a stem cell marker with neurodevelopmental role that contains two canonical RRMs. In the circuit, Musashi-1 is regulated transcriptionally and works as an allosteric translation repressor thanks to a specific interaction with the N-terminal coding region of a messenger RNA and its structural plasticity to respond to fatty acids. We fully characterized the genetic system at the population and single-cell levels showing a significant fold change in reporter expression, and the underlying molecular mechanism by assessing the in vitro binding kinetics and in vivo functionality of a series of RNA mutants. The dynamic response of the system was well recapitulated by a bottom-up mathematical model. Moreover, we applied the post-transcriptional mechanism engineered with Musashi-1 to specifically regulate a gene within an operon, implement combinatorial regulation, and reduce protein expression noise. This work illustrates how RRM-based regulation can be adapted to simple organisms, thereby adding a new regulatory layer in prokaryotes for translation control.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Transcriptome-wide identification of single-stranded RNA binding proteins

RNA-protein interactions are precisely regulated by RNA secondary structures in various biological processes. Large-scale identification of proteins that interact with particular RNA structure is important to the RBPome. Herein, a kethoxal assisted single-stranded RNA interactome capture (KASRIC) strategy was developed to globally identify single-stranded RNA binding proteins (ssRBPs). This approach combines RNA secondary structure probing technology with the conventional method of RNA-binding proteins profiling, realizing the transcriptome-wide identification of ssRBPs. Applying KASRIC, we identified 3180 candidate RBPs and 244 candidate ssRBPs in HeLa cells. Importantly, the 244 candidate ssRBPs contained 55 previously reported ssRBPs and 189 novel ssRBPs. Function analysis of the candidate ssRBPs exhibited enrichment in cellular processes related to RNA splicing and RNA degradation. The KASRIC strategy will facilitate the investigation of RNA-protein interactions.

Biocomputational Identification of sRNAs in Leptospira interrogans Serovar Lai

Regulatory small RNAs (sRNA) are RNA transcripts that are not translated into proteins but act as functional RNAs. Pathogenic Leptospira cause an epidemic spirochaetal zoonosis, Leptospirosis. It is speculated that Leptospiral sRNAs are involved in orchestrating their pathogenicity. In this study, biocomputational approach was adopted to identify Leptospiral sRNAs. In this study, two sRNA prediction programs, i.e., RNAz and nocoRNAc, were employed to screen the reference genome of Leptospira interrogans serovar Lai. Out of 126 predicted sRNAs, there are 96 cis-antisense sRNAs, 28 trans-encoded sRNAs and 2 sRNAs that partially overlap with protein-coding genes in a sense orientation. To determine whether these candidates are expressed in the pathogen, they were compared with the coverage files generated from our RNA-seq datasets. It was found out that 7 predicted sRNAs are expressed in mid-log phase, stationary phase, serum stress, temperature stress and iron stress while 2 sRNAs are expressed in mid-log phase, stationary phase, serum stress, and temperature stress. Besides, their expressions were also confirmed experimentally via RT-PCR. These experimentally validated candidates were also subjected to mRNA target prediction using TargetRNA2. Taken together, our study demonstrated that biocomputational strategy can serve as an alternative or as a complementary strategy to the laborious and expensive deep sequencing methods not only to uncover putative sRNAs but also to predict their targets in bacteria. In fact, this is the first study that integrates computational approach to predict putative sRNAs in L. interrogans serovar Lai.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01050-9.

Characteristic Features of Protein Interaction with Single- and Double-Stranded RNA

The review discusses differences between the specific protein interactions with single- and double-stranded RNA molecules using the data on the structure of RNA-protein complexes. Proteins interacting with the single-stranded RNAs form contacts with RNA bases, which ensures recognition of specific nucleotide sequences. Formation of such contacts with the double-stranded RNAs is hindered, so that the proteins recognize unique conformations of the RNA spatial structure and interact mainly with the RNA sugar-phosphate backbone.

Poly(ADP-Ribosyl) Code Functions

Poly(ADP-ribosyl)ation plays a key role in cellular metabolism. Covalent poly(ADP-ribosyl)ation affects the activity of the proteins engaged in DNA repair, chromatin structure regulation, gene expression, RNA processing, ribosome biogenesis, and protein translation. Non-covalent PAR-dependent interactions are involved in the various types of cellular response to stress and viral infection, such as inflammation, hormonal signaling, and the immune response. The review discusses how structurally different poly(ADP-ribose) (PAR) molecules composed of identical monomers can differentially participate in various cellular processes acting as the so-called "PAR code." The article describes the ability of PAR polymers to form functional biomolecular clusters through a phase-separation in response to various signals. This phase-separation contributes to rapid spatial segregation of biochemical processes and effective recruitment of the necessary components. The cellular PAR level is tightly controlled by a network of regulatory proteins: PAR code writers, readers, and erasers. Impaired PAR metabolism is associated with the development of pathological processes causing oncological, cardiovascular, and neurodegenerative diseases. Pharmacological correction of the PAR level may represent a new approach to the treatment of various diseases.

Weak binding to the A2RE RNA rigidifies hnRNPA2 RRMs and reduces liquid-liquid phase separation and aggregation

hnRNPA2 is a major component of mRNA transport granules in oligodendrocytes and neurons. However, the structural details of how hnRNPA2 binds the A2 recognition element (A2RE) and if this sequence stimulates granule formation by enhancing phase separation of hnRNPA2 has not yet been studied. Using solution NMR and biophysical studies, we find that each of the two individual RRMs retain the domain structure observed in complex with RNA but are not rigidly confined (i.e. they move independently) in solution in the absence of RNA. hnRNPA2 RRMs bind the minimal rA2RE11 weakly but at least, and most likely, two hnRNPA2 molecules are able to simultaneously bind the longer 21mer myelin basic protein A2RE. Upon binding of the RNA, NMR chemical shift deviations are observed in both RRMs, suggesting both play a role in binding the A2RE11. Interestingly, addition of short A2RE RNAs or longer RNAs containing this sequence completely prevents in vitro phase separation of full-length hnRNPA2 and aggregation of the disease-associated mutants. These findings suggest that RRM interactions with specific recognition sequences alone do not account for nucleating granule formation, consistent with models where multivalent protein:RNA and protein:protein contacts form across many sites in granule proteins and long RNA transcripts.

Phosphorothioate modified oligonucleotide-protein interactions

Antisense oligonucleotides (ASOs) interact with target RNAs via hybridization to modulate gene expression through different mechanisms. ASO therapeutics are chemically modified and include phosphorothioate (PS) backbone modifications and different ribose and base modifications to improve pharmacological properties. Modified PS ASOs display better binding affinity to the target RNAs and increased binding to proteins. Moreover, PS ASO protein interactions can affect many aspects of their performance, including distribution and tissue delivery, cellular uptake, intracellular trafficking, potency and toxicity. In this review, we summarize recent progress in understanding PS ASO protein interactions, highlighting the proteins with which PS ASOs interact, the influence of PS ASO protein interactions on ASO performance, and the structure activity relationships of PS ASO modification and protein interactions. A detailed understanding of these interactions can aid in the design of safer and more potent ASO drugs, as illustrated by recent findings that altering ASO chemical modifications dramatically improves therapeutic index.

Site-specific replacement of phosphorothioate with alkyl phosphonate linkages enhances the therapeutic profile of gapmer ASOs by modulating interactions with cellular proteins

Phosphorothioate-modified antisense oligonucleotides (PS-ASOs) interact with a host of plasma, cell-surface and intracellular proteins which govern their therapeutic properties. Given the importance of PS backbone for interaction with proteins, we systematically replaced anionic PS-linkages in toxic ASOs with charge-neutral alkylphosphonate linkages. Site-specific incorporation of alkyl phosphonates altered the RNaseH1 cleavage patterns but overall rates of cleavage and activity versus the on-target gene in cells and in mice were only minimally affected. However, replacing even one PS-linkage at position 2 or 3 from the 5'-side of the DNA-gap with alkylphosphonates reduced or eliminated toxicity of several hepatotoxic gapmer ASOs. The reduction in toxicity was accompanied by the absence of nucleolar mislocalization of paraspeckle protein P54nrb, ablation of P21 mRNA elevation and caspase activation in cells, and hepatotoxicity in mice. The generality of these observations was further demonstrated for several ASOs versus multiple gene targets. Our results add to the types of structural modifications that can be used in the gap-region to enhance ASO safety and provide insights into understanding the biochemistry of PS ASO protein interactions.

RNA virus evasion of nonsense-mediated decay

Nonsense-mediated decay (NMD) is a host RNA control pathway that removes aberrant transcripts with long 3' untranslated regions (UTRs) due to premature termination codons (PTCs) that arise through mutation or defective splicing. To maximize coding potential, RNA viruses often contain internally located stop codons that should also be prime targets for NMD. Using an agroinfiltration-based NMD assay in Nicotiana benthamiana, we identified two segments conferring NMD-resistance in the carmovirus Turnip crinkle virus (TCV) genome. The ribosome readthrough structure just downstream of the TCV p28 termination codon stabilized an NMD-sensitive reporter as did a frameshifting element from umbravirus Pea enation mosaic virus. In addition, a 51-nt unstructured region (USR) at the beginning of the TCV 3' UTR increased NMD-resistance 3-fold when inserted into an unrelated NMD-sensitive 3' UTR. Several additional carmovirus 3' UTRs also conferred varying levels of NMD resistance depending on the construct despite no sequence similarity in the analogous region. Instead, these regions displayed a marked lack of RNA structure immediately following the NMD-targeted stop codon. NMD-resistance was only slightly reduced by conversion of 19 pyrimidines in the USR to purines, but resistance was abolished when a 2-nt mutation was introduced downstream of the USR that substantially increased the secondary structure in the USR through formation of a stable hairpin. The same 2-nt mutation also enhanced the NMD susceptibility of a subgenomic RNA expressed independently of the genomic RNA. The conserved lack of RNA structure among most carmoviruses at the 5' end of their 3' UTR could serve to enhance subgenomic RNA stability, which would increase expression of the encoded capsid protein that also functions as the RNA silencing suppressor. These results demonstrate that the TCV genome has features that are inherently NMD-resistant and these strategies could be widespread among RNA viruses and NMD-resistant host mRNAs with long 3' UTRs.

Discrimination against RNA Backbones by a ssDNA Binding Protein

Pot1 is the shelterin component responsible for the protection of the single-stranded DNA (ssDNA) overhang at telomeres in nearly all eukaryotic organisms. The C-terminal domain of the DNA-binding domain, Pot1pC, exhibits non-specific ssDNA recognition, achieved through thermodynamically equivalent alternative binding conformations. Given this flexibility, it is unclear how specificity for ssDNA over RNA, an activity required for biological function, is achieved. Examination of the ribose-position specificity of Pot1pC shows that ssDNA specificity is additive but not uniformly distributed across the ligand. High-resolution structures of several Pot1pC complexes with RNA-DNA chimeric ligands reveal Pot1pC discriminates against RNA by utilizing non-compensatory binding modes that feature significant rearrangement of the binding interface. These alternative conformations, accessed through both ligand and protein flexibility, recover much, but not all, of the binding energy, leading to the observed reduction in affinities. These findings suggest that intermolecular interfaces are remarkably sophisticated in their tuning of specificity toward flexible ligands.

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.

Tying up the Ends: Plasticity in the Recognition of Single-Stranded DNA at Telomeres

Telomeres terminate nearly exclusively in single-stranded DNA (ssDNA) overhangs comprised of the G-rich 3' end. This overhang varies widely in length from species to species, ranging from just a few bases to several hundred nucleotides. These overhangs are not merely a remnant of DNA replication but rather are the result of complex further processing. Proper management of the telomeric overhang is required both to deter the action of the DNA damage machinery and to present the ends properly to the replicative enzyme telomerase. This Current Topic addresses the biochemical and structural features used by the proteins that manage these variable telomeric overhangs. The Pot1 protein tightly binds the single-stranded overhang, preventing DNA damage sensors from binding. Pot1 also orchestrates the access of telomerase to that same substrate. The remarkable plasticity of the binding interface exhibited by the Schizosaccharomyces pombe Pot1 provides mechanistic insight into how these roles may be accomplished, and disease-associated mutations clustered around the DNA-binding interface in the hPOT1 highlight the importance of this function. The budding yeast Cdc13-Stn1-Ten1, a telomeric RPA complex closely associated with telomere function, also interacts with ssDNA in a fashion that allows degenerate sequences to be recognized. A related human complex composed of hCTC1, hSTN1, and hTEN1 has recently emerged with links to both telomere maintenance and general DNA replication and also exhibits mutations associated with telomere pathologies. Overall, these sequence-specific ssDNA binders exhibit a range of recognition properties that allow them to perform their unique biological functions.

Fragment-based modelling of single stranded RNA bound to RNA recognition motif containing proteins

Protein-RNA complexes are important for many biological processes. However, structural modeling of such complexes is hampered by the high flexibility of RNA. Particularly challenging is the docking of single-stranded RNA (ssRNA). We have developed a fragment-based approach to model the structure of ssRNA bound to a protein, based on only the protein structure, the RNA sequence and conserved contacts. The conformational diversity of each RNA fragment is sampled by an exhaustive library of trinucleotides extracted from all known experimental protein–RNA complexes. The method was applied to ssRNA with up to 12 nucleotides which bind to dimers of the RNA recognition motifs (RRMs), a highly abundant eukaryotic RNA-binding domain. The fragment based docking allows a precise de novo atomic modeling of protein-bound ssRNA chains. On a benchmark of seven experimental ssRNA–RRM complexes, near-native models (with a mean heavy-atom deviation of <3 Å from experiment) were generated for six out of seven bound RNA chains, and even more precise models (deviation < 2 Å) were obtained for five out of seven cases, a significant improvement compared to the state of the art. The method is not restricted to RRMs but was also successfully applied to Pumilio RNA binding proteins.

The role of RNA conformation in RNA-protein recognition

Interactions between protein and RNA play a key role in many biological processes in the gene expression pathway. Those interactions are mediated through a variety of RNA-binding protein domains, among them the highly abundant RNA recognition motif (RRM). Here we studied protein-RNA complexes from different RNA binding domain families solved by NMR and x-ray crystallography. Characterizing the structural properties of the RNA at the binding interfaces revealed an unexpected number of nucleotides with unusual RNA conformations, specifically found in RNA-RRM complexes. Moreover, we observed that the RNA nucleotides that are directly involved in interactions with the RRM domains, via hydrogen bonds and hydrophobic contacts, are significantly enriched with unique RNA conformations. Further examination of the sequences binding the RRM domain showed a preference for G nucleotides in syn conformation to precede or to follow U nucleotides in the anti-conformation, and U nucleotides in C2' endo conformation to precede U and G nucleotides possessing the more common C3' endo conformation. These findings imply a possible mode of RNA recognition by the RRM domains which enables the recognition of a wide variety of different RNA sequences and shapes. Overall, this study suggests an additional way by which the RRM domain recognizes its RNA target, involving a conformational readout.

The sweet side of RNA regulation: glyceraldehyde-3-phosphate dehydrogenase as a noncanonical RNA-binding protein

The glycolytic protein, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), has a vast array of extraglycolytic cellular functions, including interactions with nucleic acids. GAPDH has been implicated in the translocation of transfer RNA (tRNA), the regulation of cellular messenger RNA (mRNA) stability and translation, as well as the regulation of replication and gene expression of many single-stranded RNA viruses. A growing body of evidence supports GAPDH-RNA interactions serving as part of a larger coordination between intermediary metabolism and RNA biogenesis. Despite the established role of GAPDH in nucleic acid regulation, it is still unclear how and where GAPDH binds to its RNA targets, highlighted by the absence of any conserved RNA-binding sequences. This review will summarize our current understanding of GAPDH-mediated regulation of RNA function. WIREs RNA 2016, 7:53-70. doi: 10.1002/wrna.1315 For further resources related to this article, please visit the WIREs website.

Hybridization Properties of RNA Containing 8-Methoxyguanosine and 8-Benzyloxyguanosine

Modified nucleobase analogues can serve as powerful tools for changing physicochemical and biological properties of DNA or RNA. Guanosine derivatives containing bulky substituents at 8 position are known to adopt syn conformation of N-glycoside bond. On the contrary, in RNA the anti conformation is predominant in Watson-Crick base pairing. In this paper two 8-substituted guanosine derivatives, 8-methoxyguanosine and 8-benzyloxyguanosine, were synthesized and incorporated into oligoribonucleotides to investigate their influence on the thermodynamic stability of RNA duplexes. The methoxy and benzyloxy substituents are electron-donating groups, decreasing the rate of depurination in the monomers, as confirmed by N-glycoside bond stability assessments. Thermodynamic stability studies indicated that substitution of guanosine by 8-methoxy- or 8-benzyloxyguanosine significantly decreased the thermodynamic stability of RNA duplexes. Moreover, the presence of 8-substituted guanosine derivatives decreased mismatch discrimination. Circular dichroism spectra of modified RNA duplexes exhibited patterns typical for A-RNA geometry.

The ins and outs of lncRNA structure: How, why and what comes next?

The field of structural biology has the unique advantage of being able to provide a comprehensive picture of biological mechanisms at the molecular and atomic level. Long noncoding RNAs (lncRNAs) represent the new frontier in the molecular biology of complex organisms yet remain the least characterised of all the classes of RNA. Thousands of new lncRNAs are being reported each year yet very little structural data exists for this rapidly expanding field. The length of lncRNAs ranges from 200 nt to over 100 kb in length and they generally exhibit low cellular abundance. Therefore, obtaining sufficient quantities of lncRNA to use for structural analysis is challenging. However, as technologies develop structures of lncRNAs are starting to emerge providing important information regarding their mechanism of action. Here we review the current methods used to determine the structure of lncRNA and lncRNA:protein complexes and describe the significant contribution structural biology has and will make to the field of lncRNA research. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

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.

C. elegans RNA-binding protein GLD-1 recognizes its multiple targets using sequence, context, and structural information to repress translation

Caenorhabditis elegans GLD-1, a maxi-KH motif containing RNA-binding protein, has various functions mainly during female germ cell development, suggesting that it likely controls the expression of a selective group of maternal mRNAs. To gain an insight into how GLD-1 specifically recognizes these mRNA targets, we identified 38 biochemically proven GLD-1 binding regions from multiple mRNA targets that are among over 100 putative targets co-immunoprecipitated with GLD-1. The sequence information of these regions revealed three over-represented and phylogenetically conserved sequence motifs. We found that two of the motifs, one of which is novel, are important for GLD-1 binding in several GLD-1 binding regions but not in other regions. Further analyses indicate that the importance of one of the sequence motifs is dependent on two aspects: (1) surrounding sequence information, likely acting as an accessory feature for GLD-1 to efficiently select the sequence motif and (2) RNA secondary structural environment where the sequence motif resides, which likely provides "binding-site accessibility" for GLD-1 to effectively recognize its targets. Our data suggest some mRNAs recruit GLD-1 by a distinct mechanism, which involves more than one sequence motif that needs to be embedded in the correct context and structural environment.

Finding the target sites of RNA-binding proteins

RNA–protein interactions differ from DNA–protein interactions because of the central role of RNA secondary structure. Some RNA-binding domains (RBDs) recognize their target sites mainly by their shape and geometry and others are sequence-specific but are sensitive to secondary structure context. A number of small- and large-scale experimental approaches have been developed to measure RNAs associated in vitro and in vivo with RNA-binding proteins (RBPs). Generalizing outside of the experimental conditions tested by these assays requires computational motif finding. Often RBP motif finding is done by adapting DNA motif finding methods; but modeling secondary structure context leads to better recovery of RBP-binding preferences. Genome-wide assessment of mRNA secondary structure has recently become possible, but these data must be combined with computational predictions of secondary structure before they add value in predicting in vivo binding. There are two main approaches to incorporating structural information into motif models: supplementing primary sequence motif models with preferred secondary structure contexts (e.g., MEMERIS and RNAcontext) and directly modeling secondary structure recognized by the RBP using stochastic context-free grammars (e.g., CMfinder and RNApromo). The former better reconstruct known binding preferences for sequence-specific RBPs but are not suitable for modeling RBPs that recognize shape and geometry of RNAs. Future work in RBP motif finding should incorporate interactions between multiple RBDs and multiple RBPs in binding to RNA. WIREs RNA 2014, 5:111–130. doi: 10.1002/wrna.1201

Nonspecific recognition is achieved in Pot1pC through the use of multiple binding modes

Pot1 is the protein responsible for binding to and protecting the 3' single-stranded DNA (ssDNA) overhang at most eukaryotic telomeres. Here, we present the crystal structure of one of the two oligonucleotide/oligosaccharide-binding folds (Pot1pC) that make up the ssDNA-binding domain in S. pombe Pot1. Comparison with the homologous human domain reveals unexpected structural divergence in the mode of ligand binding that explains the differing ligand requirements between species. Despite the presence of apparently base-specific hydrogen bonds, Pot1pC is able to bind a wide range of ssDNA sequences with thermodynamic equivalence. To address how Pot1pC binds ssDNA with little to no specificity, multiple structures of Pot1pC bound to noncognate ssDNA ligands were solved. These structures reveal that this promiscuity is implemented through new binding modes that thermodynamically compensate for base-substitutions through alternate stacking interactions and new H-bonding networks.

Inhibition of pre-mRNA splicing by a synthetic Blom7α-interacting small RNA

Originally the novel protein Blom7α was identified as novel pre-mRNA splicing factor that interacts with SNEV^Prp19/Pso4, an essential protein involved in extension of human endothelial cell life span, DNA damage repair, the ubiquitin-proteasome system, and pre-mRNA splicing. Blom7α belongs to the heteronuclear ribonucleoprotein K homology (KH) protein family, displaying 2 KH domains, a well conserved and widespread RNA-binding motif. In order to identify specific sequence binding motifs, we here used Systematic Evolution of Ligands by Exponential Enrichment (SELEX) with a synthetic RNA library. Besides sequence motifs like (U/A)_1–4 C_2–6 (U/A)_1–5, we identified an AC-rich RNA-aptamer that we termed AK48 (Aptamer KH-binding 48), binding to Blom7α with high affinity. Addition of AK48 to pre-mRNA splicing reactions in vitro inhibited the formation of mature spliced mRNA and led to a slight accumulation of the H complex of the spliceosome. These results suggest that the RNA binding activity of Blom7α might be required for pre-mRNA splicing catalysis. The inhibition of in-vitro splicing by the small RNA AK48 indicates the potential use of small RNA molecules in targeting the spliceosome complex as a novel target for drug development.

Principles of mRNA transport in yeast

mRNA localization and localized translation is a common mechanism by which cellular asymmetry is achieved. In higher eukaryotes the mRNA transport machinery is required for such diverse processes as stem cell division and neuronal plasticity. Because mRNA localization in metazoans is highly complex, studies at the molecular level have proven to be cumbersome. However, active mRNA transport has also been reported in fungi including Saccharomyces cerevisiae, Ustilago maydis and Candida albicans, in which these events are less difficult to study. Amongst them, budding yeast S. cerevisiae has yielded mechanistic insights that exceed our understanding of other mRNA localization events to date. In contrast to most reviews, we refrain here from summarizing mRNA localization events from different organisms. Instead we give an in-depth account of ASH1 mRNA localization in budding yeast. This approach is particularly suited to providing a more holistic view of the interconnection between the individual steps of mRNA localization, from transcriptional events to cytoplasmic mRNA transport and localized translation. Because of our advanced mechanistic understanding of mRNA localization in yeast, the present review may also be informative for scientists working, for example, on mRNA localization in embryogenesis or in neurons.

Protein-RNA and protein-protein recognition by dual KH1/2 domains of the neuronal splicing factor Nova-1

Nova onconeural antigens are neuron-specific RNA-binding proteins implicated in paraneoplastic opsoclonus-myoclonus-ataxia (POMA) syndrome. Nova harbors three K-homology (KH) motifs implicated in alternate splicing regulation of genes involved in inhibitory synaptic transmission. We report the crystal structure of the first two KH domains (KH1/2) of Nova-1 bound to an in vitro selected RNA hairpin, containing a UCAG-UCAC high-affinity binding site. Sequence-specific intermolecular contacts in the complex involve KH1 and the second UCAC repeat, with the RNA scaffold buttressed by interactions between repeats. Whereas the canonical RNA-binding surface of KH2 in the above complex engages in protein-protein interactions in the crystalline state, the individual KH2 domain can sequence-specifically target the UCAC RNA element in solution. The observed antiparallel alignment of KH1 and KH2 domains in the crystal structure of the complex generates a scaffold that could facilitate target pre-mRNA looping on Nova binding, thereby potentially explaining Nova's functional role in splicing regulation.

Solvent effects on structural and thermochemical properties of p53 tumor-suppressor gene: a molecular modeling approach in drug design

The p53 tumor-suppressor protein is a cellular phosphoprotein and a negative regulator of cell growth. Most p53 mutations occur in exons 5-8 within the DNA-binding domain. Therefore, p53 can potentially be targeted with novel drugs designed to bind to a mutation and restore its stability or wild-type conformation. For the current study, Hartree-Fock calculations were used to investigate the solvent-induced effects of five different solvent media (acetone, ethanol, methanol, dimethyl sulfoxide, and water) on the thermochemical parameters and relative energies, and on the multinuclear nuclear magnetic resonance shielding tensors of oxygen, nitrogen, and phosphorus nuclei, of GAT. To understand how the solvent affects the mutation region (the "hot spot") of p53, the relative energies of GAT in selected solvent media were determined. Some biological evidence suggested the structural stabilities of hot spots of GAT have the optimum temperature and solvent type for mutation. All the authors' findings are in accordance with common biological phenomena. Another important objective of this study was to compare the hydration Gibbs free energies of CUA and GAT in water using two different approaches where the solvent was treated as a continuum of the constant at different levels of Hartree-Fock theory. The Gibbs hydration energy values obtained in water with the polarized continuum model directly applied on the isolated CUA and GAT sequences were compared with those determined from the hydrated models with four, six, and eight water molecule clusters around the hot spots uracil and adenine. The clustered structures of water molecules around the hot spots of GAT (in DNA level) and CUA (in transcriptional level) were found to be energetically favored. The results of this study provide a reliable insight into the nature of mutation processes, which is of utmost importance for the study of biochemical structures, and provide a basis for drug design.

Investigation of solvent effect and NMR shielding tensors of p53 tumor-suppressor gene in drug design

The p53 tumor-suppressor gene encodes a nuclear phosphoprotein with cancer- inhibiting properties. The most probable cancerous mutations occur as point mutations in exons 5 up to 8 of p53, as a base pair substitution that encompasses CUA and GAT sequences. As DNA drug design represents a direct genetic treatment of cancer, in the research reported computational drug design was carried out to explore, at the Hartree-Fock level, effects of solvents on the thermochemical properties and nuclear magnetic resonance (NMR) shielding tensors of some atoms of CUA involved in the hydrogen-bonding network. The observed NMR shielding variations of the solutes caused by solvent change seemed significant and were attributed to solvent polarity, and solute-solvent and solvent-solute hydrogen-bonding interactions. The results provide a reliable insight into the nature of mutation processes. However, to improve our knowledge of the hydration pattern more rigorous computations of the hydrated complexes are needed.

Bioinformatics for RNomics

Rapid improvements in high-throughput experimental technologies make it nowadays possible to study the expression, as well as changes in expression, of whole transcriptomes under different environmental conditions in a detailed view. We describe current approaches to identify genome-wide functional RNA transcripts (experimentally as well as computationally), and focus on computational methods that may be utilized to disclose their function. While genome databases offer a wealth of information about known and putative functions for protein-coding genes, functional information for novel non-coding RNA genes is almost nonexistent. This is mainly explained by the lack of established software tools to efficiently reveal the function and evolutionary origin of non-coding RNA genes. Here, we describe in detail computational approaches one may follow to annotate and classify an RNA transcript.

A portable RNA sequence whose recognition by a synthetic antibody facilitates structural determination

RNA crystallization and phasing represent major bottlenecks in RNA structure determination. Seeking to exploit antibody fragments as RNA crystallization chaperones, we have used an arginine-enriched synthetic Fab library displayed on phage to obtain Fabs against the class I ligase ribozyme. We solved the structure of a Fab:ligase complex at 3.1Å using molecular replacement with Fab coordinates, confirming the ribozyme architecture and revealing the chaperone’s role in RNA recognition and crystal contacts. The epitope resides in the GAAACAC sequence that caps the P5 helix and retains high-affinity Fab binding within the context of other structured RNAs. This portable epitope provides a new RNA crystallization chaperone system that easily can be screened in parallel to the U1A RNA-binding protein, with the advantages of the smaller size of the loop and high molecular weight, large surface area, and phasing power provided by Fabs.

Structural characterization of naturally occurring RNA single mismatches

RNA is known to be involved in several cellular processes; however, it is only active when it is folded into its correct 3D conformation. The folding, bending and twisting of an RNA molecule is dependent upon the multitude of canonical and non-canonical secondary structure motifs. These motifs contribute to the structural complexity of RNA but also serve important integral biological functions, such as serving as recognition and binding sites for other biomolecules or small ligands. One of the most prevalent types of RNA secondary structure motifs are single mismatches, which occur when two canonical pairs are separated by a single non-canonical pair. To determine sequence–structure relationships and to identify structural patterns, we have systematically located, annotated and compared all available occurrences of the 30 most frequently occurring single mismatch-nearest neighbor sequence combinations found in experimentally determined 3D structures of RNA-containing molecules deposited into the Protein Data Bank. Hydrogen bonding, stacking and interaction of nucleotide edges for the mismatched and nearest neighbor base pairs are described and compared, allowing for the identification of several structural patterns. Such a database and comparison will allow researchers to gain insight into the structural features of unstudied sequences and to quickly look-up studied sequences.

The role of G-quadruplex in RNA metabolism: involvement of FMRP and FMR2P

Regulation of post-transcriptional gene expression is a cellular process that is accomplished through the activity of multiple mRNP (messenger RiboNucleoProtein) complexes which are composed of mRNA-binding proteins and RNA molecules interacting with those proteins. The specificity of these interactions is mediated by the ability of the RNA-binding proteins to precisely recognize and bind RNA sequences or structures. Alterations of their function may have some dramatic consequences, resulting in different pathologies. An increasing body of data is emerging showing the impact of a G-quadruplex forming structure in the maturation and expression of some RNA molecules. We review here the role of the G-quadruplex RNA structure in the regulation of translation and splicing, when it interacts with two RNA-binding proteins: FMRP (Fragile X Mental Retardation Protein) and FMR2P (Fragile X Mental Retardation 2 protein). Impaired expression of these proteins causes two forms of intellectual disability: the Fragile X Mental Retardation syndrome (FXS) and the FRAXE-associated mental retardation (FRAXE), respectively. FMRP is involved in different steps of RNA metabolism and, in particular, in translational regulation. FMR2P has been initially described as a transcription factor and we recently showed also its role in regulation of alternative splicing. By the study of the functional significance of the interaction of both FMRP and FMR2P with a G-quadruplex forming RNA we were able to show an impact of this structure in translational regulation and also in splicing, behaving as an Exonic Splicing Enhancer.

A novel RNA-binding peptide regulates the establishment of the Medicago truncatula-Sinorhizobium meliloti nitrogen-fixing symbiosis

Plants use a variety of small peptides for cell to cell communication during growth and development. Leguminous plants are characterized by their ability to develop nitrogen-fixing nodules via an interaction with symbiotic bacteria. During nodule organogenesis, several so-called nodulin genes are induced, including large families that encode small peptides. Using a three-hybrid approach in yeast cells, we identified two new small nodulins, MtSNARP1 and MtSNARP2 (for small nodulin acidic RNA-binding protein), which interact with the RNA of MtENOD40, an early induced nodulin gene showing conserved RNA secondary structures. The SNARPs are acidic peptides showing single-stranded RNA-binding activity in vitro and are encoded by a small gene family in Medicago truncatula. These peptides exhibit two new conserved motifs and a putative signal peptide that redirects a GFP fusion to the endoplasmic reticulum both in protoplasts and during symbiosis, suggesting they are secreted. MtSNARP2 is expressed in the differentiating region of the nodule together with several early nodulin genes. MtSNARP2 RNA interference (RNAi) transgenic roots showed aberrant early senescent nodules where differentiated bacteroids degenerate rapidly. Hence, a functional symbiotic interaction may be regulated by secreted RNA-binding peptides.

Bruno protein contains an expanded RNA recognition motif

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

Crystal structure of Escherichia coli polynucleotide phosphorylase core bound to RNase E, RNA and manganese: implications for catalytic mechanism and RNA degradosome assembly

Polynucleotide phosphorylase (PNPase) is a processive exoribonuclease that contributes to messenger RNA turnover and quality control of ribosomal RNA precursors in many bacterial species. In Escherichia coli, a proportion of the PNPase is recruited into a multi-enzyme assembly, known as the RNA degradosome, through an interaction with the scaffolding domain of the endoribonuclease RNase E. Here, we report crystal structures of E. coli PNPase complexed with the recognition site from RNase E and with manganese in the presence or in the absence of modified RNA. The homotrimeric PNPase engages RNase E on the periphery of its ring-like architecture through a pseudo-continuous anti-parallel beta-sheet. A similar interaction pattern occurs in the structurally homologous human exosome between the Rrp45 and Rrp46 subunits. At the centre of the PNPase ring is a tapered channel with an adjustable aperture where RNA bases stack on phenylalanine side chains and trigger structural changes that propagate to the active sites. Manganese can substitute for magnesium as an essential co-factor for PNPase catalysis, and our crystal structure of the enzyme in complex with manganese suggests how the metal is positioned to stabilise the transition state. We discuss the implications of these structural observations for the catalytic mechanism of PNPase, its processive mode of action, and its assembly into the RNA degradosome.

The zinc fingers of the SR-like protein ZRANB2 are single-stranded RNA-binding domains that recognize 5' splice site-like sequences

The alternative splicing of mRNA is a critical process in higher eukaryotes that generates substantial proteomic diversity. Many of the proteins that are essential to this process contain arginine/serine-rich (RS) domains. ZRANB2 is a widely-expressed and highly-conserved RS-domain protein that can regulate alternative splicing but lacks canonical RNA-binding domains. Instead, it contains 2 RanBP2-type zinc finger (ZnF) domains. We demonstrate that these ZnFs recognize ssRNA with high affinity and specificity. Each ZnF binds to a single AGGUAA motif and the 2 domains combine to recognize AGGUAA(N(x))AGGUAA double sites, suggesting that ZRANB2 regulates alternative splicing via a direct interaction with pre-mRNA at sites that resemble the consensus 5' splice site. We show using X-ray crystallography that recognition of an AGGUAA motif by a single ZnF is dominated by side-chain hydrogen bonds to the bases and formation of a guanine-tryptophan-guanine "ladder." A number of other human proteins that function in RNA processing also contain RanBP2 ZnFs in which the RNA-binding residues of ZRANB2 are conserved. The ZnFs of ZRANB2 therefore define another class of RNA-binding domain, advancing our understanding of RNA recognition and emphasizing the versatility of ZnF domains in molecular recognition.

Heterogeneous nuclear ribonuclear protein K interacts with the enterovirus 71 5' untranslated region and participates in virus replication

Enterovirus 71 (EV71) is a picornavirus that can cause severe neurological complications in children. Like other picornaviruses, the genomic RNA of EV71 contains a long 5' untranslated region (UTR). Cellular proteins interact with the EV71 5' UTR, and these interactions are important for virus replication. Using an RNA pull-down assay and proteomics approaches, this study identified the heterogeneous nuclear ribonucleoprotein K (hnRNP K) as one of the EV71 5' UTR-associated proteins. The interaction between hnRNP K and the 5' UTR was further confirmed by mapping the interaction regions to stem-loops I-II and IV in the 5' UTR. During EV71 infection, hnRNP K was enriched in the cytoplasm where virus replication occurs, whereas hnRNP K was localized in the nucleus in mock-infected cells. Viral yields were found to be significantly lower in hnRNP K knockdown cells and viral RNA synthesis was delayed in hnRNP K knockdown cells in comparison with negative-control cells treated with small interfering RNA. These results suggest that hnRNP K interacts with the EV71 5' UTR and participates in virus replication.

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.

Sequence and transcriptional study of HNRPK pseudogenes, and expression and molecular modeling analysis of hnRNP K isoforms

The heterogeneous nuclear ribonucleoproteins (hnRNPs) comprise a large family of proteins that play important roles in telomere biogenesis, DNA repair, cellular signaling, and the regulation of expression at both the transcriptional and translational levels. One of the most extensively studied hnRNP family members, hnRNP K, has been implicated in a variety of processes, including chromatin remodeling, transcription, splicing, and translation events. In this study, we analyzed processed HNRPK pseudogenes (HNRPK psi1-psi4) and coding sequences. HNRPK pseudogenes are apparently nonfunctional, and psi1 might correspond to transcripts from an ancestral gene. Phylogenetic and sequence analyses suggest that HNRP genes arose by duplication, and that new structural and sequence features expanded the functions of hnRNPs. The expression analysis of hnRNP K isoforms showed that isoform a is expressed in normal testis and in non-small cell lung cancer (NCI-H1155 NSCLC cell line), although the shorter isoform (isoform b) is expressed in different tumor cell lines (IM9 B-lymphoblastoid, Hs578T human breast cancer epithelial, T98G human glioma cell lines). Using molecular modeling, we obtained KH1 and KH3 models, which pointed to important residues for DNA-protein binding and no structural differences between isoforms a and b. To our knowledge, this is the first phylogenetic study including vertebrate HNRP genes and HNRPK pseudogenes, and the first report comparing the KH1 and KH3 domains of isoforms a and b of the hnRNP K protein. New investigations in tumor samples must be done to validate the differential expression observed here. The results shown are important because the hnRNP K protein might represent a new target for pharmacologic intervention in virus replication and cancer.

Structural insight into a molecular switch in tandem winged-helix motifs from elongation factor SelB

Elongation factor SelB is responsible for co-translational incorporation of selenocysteine (Sec) into proteins. The UGA stop codon is recoded as a Sec codon in the presence of a downstream mRNA hairpin. In prokaryotes, in addition to the EF-Tu-like N-terminal domains, a C-terminal extension containing four tandem winged-helix motifs (WH1-4) recognizes the mRNA hairpin. The 2.3-A resolution crystal structure of the Escherichia coli WH3/4 domains bound to mRNA with mutagenesis data reveal that the two WH motifs use the same structural elements to bind RNA. The structure together with the 2.6-A resolution structure of the WH1-4 domains from Moorella thermoacetica bound to RNA revealed that a salt bridge connecting WH2 to WH3 modules is disrupted upon mRNA binding. The results provide a structural basis for the molecular switch that may allow communication between tRNA and mRNA binding sites and illustrate how RNA acts as an activator of the switch. The structures show that tandem WH motifs not only provide an excellent scaffold for RNA binding but can also have an active role in the function of protein-RNA complexes.

Kinetic analysis of the entire RNA amplification process by Qbeta replicase

The kinetics of the RNA replication reaction by Qbeta replicase were investigated. Qbeta replicase is an RNA-dependent RNA polymerase responsible for replicating the RNA genome of coliphage Qbeta and plays a key role in the life cycle of the Qbeta phage. Although the RNA replication reaction using this enzyme has long been studied, a kinetic model that can describe the entire RNA amplification process has yet to be determined. In this study, we propose a kinetic model that is able to account for the entire RNA amplification process. The key to our proposed kinetic model is the consideration of nonproductive binding (i.e. binding of an enzyme to the RNA where the enzyme cannot initiate the reaction). By considering nonproductive binding and the notable enzyme inactivation we observed, the previous observations that remained unresolved could also be explained. Moreover, based on the kinetic model and the experimental results, we determined rate and equilibrium constants using template RNAs of various lengths. The proposed model and the obtained constants provide important information both for understanding the basis of Qbeta phage amplification and the applications using Qbeta replicase.

Solving the Structure of PTB in Complex with Pyrimidine Tracts: An NMR Study of Protein–RNA Complexes of Weak Affinities†

NMR spectroscopy has proven to be a powerful tool for the structure determination of protein/RNA complexes. However, the quality of these structures depends critically on the number of unambiguous intermolecular and intra-RNA nuclear Overhauser effect (NOE) constraints that can be derived. This number is often limited due to exchange phenomena that can cause signal line broadening and the fact that unambiguous NOE assignments are challenging in systems that exchange between different conformations in the intermediate to fast exchange limit. These exchange processes can include exchange between free and bound form, as well as exchange of the ligand between different binding sites on the protein. Furthermore, for the large class of RNA metabolizing proteins that bind repetitive low-complexity RNA sequences in multiple register, exchange of the protein between these overlapping binding sites introduces additional exchange pathways. Here, we describe the strategy we used to overcome these exchange processes and to reduce significantly the line width of the RNA resonances in complexes of the RNA recognition motifs (RRMs) of the polypyrimidine tract-binding protein (PTB) in complex with pyrimidine tracts and hence allowed a highly precise structure determination. This method could be employed to derive structures of other protein/single-stranded nucleic acid complexes by NMR spectroscopy. Furthermore, we have determined the affinities of the individual RRMs of PTB for pyrimidine tracts of different length and sequence. These measurements show that PTB binds preferentially to long pyrimidine tracts that contain cytosine and hence confirm the structure of PTB in complex with RNA. Furthermore, they provide quantitative insight into the question of which pyrimidine sequences within alternatively spliced pre-mRNAs will be preferentially bound by PTB.

Substrate determinants for RNA editing and editing complex interactions at a site for full-round U insertion

Multisubunit RNA editing complexes catalyze uridylate insertion/deletion RNA editing directed by complementary guide RNAs (gRNAs). Editing in trypanosome mitochondria is transcript-specific and developmentally controlled, but the molecular mechanisms of substrate specificity remain unknown. Here we used a minimal A6 pre-mRNA/gRNA substrate to define functional determinants for full-round insertion and editing complex interactions at the editing site 2 (ES2). Editing begins with pre-mRNA cleavage within an internal loop flanked by upstream and downstream duplexes with gRNA. We found that substrate recognition around the internal loop is sequence-independent and that completely artificial duplexes spanning a single helical turn are functional. Furthermore, after our report of cross-linking interactions at the deletion ES1 (35), we show for the first time editing complex contacts at an insertion ES. Our studies using site-specific ribose 2' substitutions defined 2'-hydroxyls within the (a) gRNA loop region and (b) flanking helixes that markedly stimulate both pre-mRNA cleavage and editing complex interactions at ES2. Modification of the downstream helix affected scissile bond specificity. Notably, a single 2'-hydroxyl at ES2 is essential for cleavage but dispensable for editing complex cross-linking. This study provides new insights on substrate recognition during full-round editing, including the relevance of secondary structure and the first functional association of specific (pre-mRNA and gRNA) riboses with both endonuclease cleavage and cross-linking activities of editing complexes at an ES. Importantly, most observed cross-linking interactions are both conserved and relatively stable at ES2 and ES1 in hybrid substrates. However, they were also detected as transient low-stability contacts in a non-edited transcript.

Common mode of DNA binding to cold shock domains. Crystal structure of hexathymidine bound to the domain-swapped form of a major cold shock protein from Bacillus caldolyticus

Bacterial cold shock proteins (CSPs) regulate cellular adaptation to cold stress. Functions ascribed to CSP include roles as RNA chaperones and in transcription antitermination. We present the crystal structure of the Bacillus caldolyticus CSP (Bc-Csp) in complex with hexathymidine (dT(6)) at a resolution of 1.29 A. Bound to dT(6), crystalline Bc-Csp forms a domain-swapped dimer in which beta strands 1-3 associate with strands 4 and 5 from the other subunit to form a closed beta barrel and vice versa. The globular units of dimeric Bc-Csp closely resemble the well-known structure of monomeric CSP. Structural reorganization from the monomer to the domain-swapped dimer involves a strictly localized change in the peptide bond linking Glu36 and Gly37 of Bc-Csp. Similar structural reorganizations have not been found in any other CSP or oligonucleotide/oligosaccharide-binding fold structures. Each dT(6) ligand is bound to one globular unit of Bc-Csp via an amphipathic protein surface. Individual binding subsites interact with the DNA bases through stacking and hydrogen bonding. The sugar-phosphate backbone remains solvent exposed. Based on crystallographic and biochemical studies of deoxyoligonucleotide binding to CSP, we suggest a common mode of binding of single-stranded heptanucleotide motifs to proteins containing cold shock domains, including the eukaryotic Y-box factors.

Specificity of Phage Display Selected Peptides for Modified Anticodon Stem and Loop Domains of tRNA

Protein recognition of RNA has been studied using Peptide Phage Display Libraries, but in the absence of RNA modifications. Peptides from two libraries, selected for binding the modified anticodon stem and loop (ASL) of human tRNA(LyS3) having 2-thiouridine (s(2)U34) and pseudouridine (psi39), bound the modified human ASL(Lys3)(s(2)U34;psi39) preferentially and had significant homology with RNA binding proteins. Selected peptides were narrowed to a manageable number using a less sensitive, but inexpensive assay before conducting intensive characterization. The affinity and specificity of the best binding peptide (with an N-terminal fluorescein) were characterized by fluorescence spectrophotometry. The peptide exhibited the highest binding affinity for ASL(LYS3)(s(2)U34; psi39), followed by the hypermodified ASL(Lys3) (mcm(5)s(2) U34; ms(2)t(6)A37) and the unmodified ASL(Lys3), but bound poorly to singly modified ASL(Lys3) constructs (psi39, ms(2)t(6)A37, s(2)34), ASL(Lys1,2) (t(6)A37) and Escherichia coli ASL(Glu) (s(2)U34). Thus, RNA modifications are potentially important recognition elements for proteins and can be targets for selective recognition by peptides.