Protein-RNA recognition

DeepA-RBPBS: A hybrid convolution and recurrent neural network combined with attention mechanism for predicting RBP binding site

It's important to infer the binding site of RNA-binding proteins (RBP) for understanding the interaction between RBP and its RNA targets and decipher the mechanisms of transcriptional regulation. However, experimental detection of RBP binding sites is still time-intensive and expensive. Algorithms based on machine learning can speed up detection of RBP binding sites. In this article, we propose a new deep learning method, DeepA-RBPBS, which can use RNA sequences and structural features to predict RBP binding site. DeepA-RBPBS uses CNN and BiGRU to extract sequences and structural features without long-term dependence issues. It also utilizes an attention mechanism to enhance the contribution of key features. The comparison shows that the performance of DeepA-RBPBS is better than that of the state-of-the-art predictors. In the testing on 31 datasets of CLIP-seq experiments over 19 proteins, MCC (AUC) is 8% (5%) higher than those of the latest method based on deep learning, iDeepS. We also apply DeepA-RBPBS to the target RNA data of RBPs related to diabetes (LIN28, RBFOX2, FTO, IGF2BP2, CELF1 and HuR). The results show that DeepA-RBPBS correctly predicted 41,693 samples, where iDeepS predicted 31,381 samples.Communicated by Ramaswamy H. Sarma.

NMR Studies of Retroviral Genome Packaging

Nearly all retroviruses selectively package two copies of their unspliced RNA genomes from a cellular milieu that contains a substantial excess of non-viral and spliced viral RNAs. Over the past four decades, combinations of genetic experiments, phylogenetic analyses, nucleotide accessibility mapping, in silico RNA structure predictions, and biophysical experiments were employed to understand how retroviral genomes are selected for packaging. Genetic studies provided early clues regarding the protein and RNA elements required for packaging, and nucleotide accessibility mapping experiments provided insights into the secondary structures of functionally important elements in the genome. Three-dimensional structural determinants of packaging were primarily derived by nuclear magnetic resonance (NMR) spectroscopy. A key advantage of NMR, relative to other methods for determining biomolecular structure (such as X-ray crystallography), is that it is well suited for studies of conformationally dynamic and heterogeneous systems—a hallmark of the retrovirus packaging machinery. Here, we review advances in understanding of the structures, dynamics, and interactions of the proteins and RNA elements involved in retroviral genome selection and packaging that are facilitated by NMR.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

Importance of amino acids Leu135 and Tyr236 for the interaction between EhCFIm25 and RNA: a molecular dynamics simulation study

The CFIm25 subunit of the heterotetrameric cleavage factor Im (CFIm) is a critical factor in the formation of the poly(A) tail at mRNA 3' end, regulating the recruitment of polyadenylation factors, poly(A) site selection, and cleavage/polyadenylation reactions. We previously reported the homologous protein (EhCFIm25) in Entamoeba histolytica, the protozoan causing human amoebiasis, and showed the relevance of conserved Leu135 and Tyr236 residues for RNA binding. We also identified the GUUG sequence as the recognition site of EhCFIm25. To understand the interactions network that allows the EhCFIm25 to maintain its three-dimensional structure and function, here we performed molecular dynamics simulations of wild-type (WT) and mutant proteins, alone or interacting with the GUUG molecule. Our results indicated that in the presence of the GUUG sequence, WT converged more quickly to lower RMSD values in comparison with mutant proteins. However, RMSF values showed that movements of amino acids of WT and EhCFIm25*L135 T were almost identical, interacting or not with the GUUG molecule. Interestingly, EhCFIm25*L135 T, which is the only mutant with a slight RNA binding activity experimentally, presents the same stabilization of bend structures and alpha helices as WT, notably in the C-terminus. Moreover, WT and EhCFIm25*L135 T presented almost the same number of contacts that mainly involve lysine residues interacting with the G4 nucleotide. Overall, our data proposed a clear description of the structural and mechanistic data that govern the RNA binding capacity of EhCFIm25.

Cross- and Co-Packaging of Retroviral RNAs and Their Consequences

Retroviruses belong to the family Retroviridae and are ribonucleoprotein (RNP) particles that contain a dimeric RNA genome. Retroviral particle assembly is a complex process, and how the virus is able to recognize and specifically capture the genomic RNA (gRNA) among millions of other cellular and spliced retroviral RNAs has been the subject of extensive investigation over the last two decades. The specificity towards RNA packaging requires higher order interactions of the retroviral gRNA with the structural Gag proteins. Moreover, several retroviruses have been shown to have the ability to cross-/co-package gRNA from other retroviruses, despite little sequence homology. This review will compare the determinants of gRNA encapsidation among different retroviruses, followed by an examination of our current understanding of the interaction between diverse viral genomes and heterologous proteins, leading to their cross-/co-packaging. Retroviruses are well-known serious animal and human pathogens, and such a cross-/co-packaging phenomenon could result in the generation of novel viral variants with unknown pathogenic potential. At the same time, however, an enhanced understanding of the molecular mechanisms involved in these specific interactions makes retroviruses an attractive target for anti-viral drugs, vaccines, and vectors for human gene therapy.

A game of tag: MAPS catches up on RNA interactomes

In the last few decades, small regulatory RNA (sRNA) molecules emerged as key regulators in every kingdom of life. Resolving the full targetome of sRNAs has however remained a challenge. To address this, we used an in vivo tagging MS2-affinity purification protocol coupled with RNA sequencing technology, namely MAPS, to assemble full bacterial small RNAs targetomes. The impressive potential of MAPS has been supported by a number of reports. Here, we concisely overview RNA-tagging history that preceded the development of the MAPS assay and expose the range of possible uses of this technology.

PRBP: Prediction of RNA-Binding Proteins Using a Random Forest Algorithm Combined with an RNA-Binding Residue Predictor

The prediction of RNA-binding proteins is an incredibly challenging problem in computational biology. Although great progress has been made using various machine learning approaches with numerous features, the problem is still far from being solved. In this study, we attempt to predict RNA-binding proteins directly from amino acid sequences. A novel approach, PRBP predicts RNA-binding proteins using the information of predicted RNA-binding residues in conjunction with a random forest based method. For a given protein, we first predict its RNA-binding residues and then judge whether the protein binds RNA or not based on information from that prediction. If the protein cannot be identified by the information associated with its predicted RNA-binding residues, then a novel random forest predictor is used to determine if the query protein is a RNA-binding protein. We incorporated features of evolutionary information combined with physicochemical features (EIPP) and amino acid composition feature to establish the random forest predictor. Feature analysis showed that EIPP contributed the most to the prediction of RNA-binding proteins. The results also showed that the information from the RNA-binding residue prediction improved the overall performance of our RNA-binding protein prediction. It is anticipated that the PRBP method will become a useful tool for identifying RNA-binding proteins. A PRBP Web server implementation is freely available at http://www.cbi.seu.edu.cn/PRBP/.

Sequence-Based Prediction of RNA-Binding Proteins Using Random Forest with Minimum Redundancy Maximum Relevance Feature Selection

The prediction of RNA-binding proteins is one of the most challenging problems in computation biology. Although some studies have investigated this problem, the accuracy of prediction is still not sufficient. In this study, a highly accurate method was developed to predict RNA-binding proteins from amino acid sequences using random forests with the minimum redundancy maximum relevance (mRMR) method, followed by incremental feature selection (IFS). We incorporated features of conjoint triad features and three novel features: binding propensity (BP), nonbinding propensity (NBP), and evolutionary information combined with physicochemical properties (EIPP). The results showed that these novel features have important roles in improving the performance of the predictor. Using the mRMR-IFS method, our predictor achieved the best performance (86.62% accuracy and 0.737 Matthews correlation coefficient). High prediction accuracy and successful prediction performance suggested that our method can be a useful approach to identify RNA-binding proteins from sequence information.

Conserved determinants of lentiviral genome dimerization

BACKGROUND

Retroviruses selectively package two copies of their unspliced genomes by what appears to be a dimerization-dependent RNA packaging mechanism. Dimerization of human immunodeficiency virus Type-1 (HIV-1) genomes is initiated by "kissing" interactions between GC-rich palindromic loop residues of a conserved hairpin (DIS), and is indirectly promoted by long-range base pairing between residues overlapping the gag start codon (AUG) and an upstream Unique 5' element (U5). The DIS and U5:AUG structures are phylogenetically conserved among divergent retroviruses, suggesting conserved functions. However, some studies suggest that the DIS of HIV-2 does not participate in dimerization, and that U5:AUG pairing inhibits, rather than promotes, genome dimerization. We prepared RNAs corresponding to native and mutant forms of the 5' leaders of HIV-1 (NL4-3 strain), HIV-2 (ROD strain), and two divergent strains of simian immunodeficiency virus (SIV; cpz-TAN1 and -US strains), and probed for potential roles of the DIS and U5:AUG base pairing on intrinsic and NC-dependent dimerization by mutagenesis, gel electrophoresis, and NMR spectroscopy.

RESULTS

Dimeric forms of the native HIV-2 and SIV leaders were only detectable using running buffers that contained Mg(2+), indicating that these dimers are more labile than that of the HIV-1 leader. Mutations designed to promote U5:AUG base pairing promoted dimerization of the HIV-2 and SIV RNAs, whereas mutations that prevented U5:AUG pairing inhibited dimerization. Chimeric HIV-2 and SIV leader RNAs containing the dimer-promoting loop of HIV-1 (DIS) exhibited HIV-1 leader-like dimerization properties, whereas an HIV-1NL4-3 mutant containing the SIVcpzTAN1 DIS loop behaved like the SIVcpzTAN1 leader. The cognate NC proteins exhibited varying abilities to promote dimerization of the retroviral leader RNAs, but none were able to convert labile dimers to non-labile dimers.

CONCLUSIONS

The finding that U5:AUG formation promotes dimerization of the full-length HIV-1, HIV-2, SIVcpzUS, and SIVcpzTAN1 5' leaders suggests that these retroviruses utilize a common RNA structural switch mechanism to modulate function. Differences in native and NC-dependent dimerization propensity and lability are due to variations in the compositions of the DIS loop residues rather than other sequences within the leader RNAs. Although NC is a well-known RNA chaperone, its role in dimerization has the hallmarks of a classical riboswitch.

Protein-specific prediction of mRNA binding using RNA sequences, binding motifs and predicted secondary structures

Background

RNA-binding proteins interact with specific RNA molecules to regulate important cellular processes. It is therefore necessary to identify the RNA interaction partners in order to understand the precise functions of such proteins. Protein-RNA interactions are typically characterized using in vivo and in vitro experiments but these may not detect all binding partners. Therefore, computational methods that capture the protein-dependent nature of such binding interactions could help to predict potential binding partners in silico.

Results

We have developed three methods to predict whether an RNA can interact with a particular RNA-binding protein using support vector machines and different features based on the sequence (the Oli method), the motif score (the OliMo method) and the secondary structure (the OliMoSS method). We applied these approaches to different experimentally-derived datasets and compared the predictions with RNAcontext and RPISeq. Oli outperformed OliMoSS and RPISeq, confirming our protein-specific predictions and suggesting that tetranucleotide frequencies are appropriate discriminative features. Oli and RNAcontext were the most competitive methods in terms of the area under curve. A precision-recall curve analysis achieved higher precision values for Oli. On a second experimental dataset including real negative binding information, Oli outperformed RNAcontext with a precision of 0.73 vs. 0.59.

Conclusions

Our experiments showed that features based on primary sequence information are sufficiently discriminating to predict specific RNA-protein interactions. Sequence motifs and secondary structure information were not necessary to improve these predictions. Finally we confirmed that protein-specific experimental data concerning RNA-protein interactions are valuable sources of information that can be used for the efficient training of models for in silico predictions. The scripts are available upon request to the corresponding author.

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.

Identification of a minimal region of the HIV-1 5'-leader required for RNA dimerization, NC binding, and packaging

Assembly of human immunodeficiency virus type 1 (HIV-1) particles is initiated in the cytoplasm by the formation of a ribonucleoprotein complex comprising the dimeric RNA genome and a small number of viral Gag polyproteins. Genomes are recognized by the nucleocapsid (NC) domains of Gag, which interact with packaging elements believed to be located primarily within the 5'-leader (5'-L) of the viral RNA. Recent studies revealed that the native 5'-L exists as an equilibrium of two conformers, one in which dimer-promoting residues and NC binding sites are sequestered and packaging is attenuated, and one in which these sites are exposed and packaging is promoted. To identify the elements within the dimeric 5'-L that are important for packaging, we generated HIV-1 5'-L RNAs containing mutations and deletions designed to eliminate substructures without perturbing the overall structure of the leader and examined effects of the mutations on RNA dimerization, NC binding, and packaging. Our findings identify a 159-residue RNA packaging signal that possesses dimerization and NC binding properties similar to those of the intact 5'-L and contains elements required for efficient RNA packaging.

Structural determinants and mechanism of HIV-1 genome packaging

Like all retroviruses, the human immunodeficiency virus selectively packages two copies of its unspliced RNA genome, both of which are utilized for strand-transfer-mediated recombination during reverse transcription-a process that enables rapid evolution under environmental and chemotherapeutic pressures. The viral RNA appears to be selected for packaging as a dimer, and there is evidence that dimerization and packaging are mechanistically coupled. Both processes are mediated by interactions between the nucleocapsid domains of a small number of assembling viral Gag polyproteins and RNA elements within the 5'-untranslated region of the genome. A number of secondary structures have been predicted for regions of the genome that are responsible for packaging, and high-resolution structures have been determined for a few small RNA fragments and protein-RNA complexes. However, major questions regarding the RNA structures (and potentially the structural changes) that are responsible for dimeric genome selection remain unanswered. Here, we review efforts that have been made to identify the molecular determinants and mechanism of human immunodeficiency virus type 1 genome packaging.

Analysis of electric moments of RNA-binding proteins: implications for mechanism and prediction

Background

Protein-RNA interactions play important role in many biological processes such as gene regulation, replication, protein synthesis and virus assembly. Although many structures of various types of protein-RNA complexes have been determined, the mechanism of protein-RNA recognition remains elusive. We have earlier shown that the simplest electrostatic properties viz. charge, dipole and quadrupole moments, calculated from backbone atomic coordinates of proteins are biased relative to other proteins, and these quantities can be used to identify DNA-binding proteins. Closely related, RNA-binding proteins are investigated in this study. In particular, discrimination between various types of RNA-binding proteins, evolutionary conservation of these bulk electrostatic features and effect of conformational changes by complex formation are investigated. Basic binding mechanism of a putative RNA-binding protein (HI1333 from Haemophilus influenza) is suggested as a potential application of this study.

Results

We found that similar to DNA-binding proteins (DBPs), RNA-binding proteins (RBPs) also show significantly higher values of electric moments. However, higher moments in RBPs are found to strongly depend on their functional class: proteins binding to ribosomal RNA (rRNA) constitute the only class with all three of the properties (charge, dipole and quadrupole moments) being higher than control proteins. Neural networks were trained using leave-one-out cross-validation to predict RBPs from control data as well as pair-wise classification capacity between proteins binding to various RNA types. RBPs and control proteins reached up to 78% accuracy measured by the area under the ROC curve. Proteins binding to rRNA are found to be best distinguished (AUC = 79%). Changes in dipole and quadrupole moments between unbound and bound structures were small and these properties are found to be robust under complex formation.

Conclusions

Bulk electric moments of proteins considered here provide insights into target recognition by RNA-binding proteins, as well as ability to recognize one type of RBP from others. These results help in understanding the mechanism of protein-RNA recognition, and identifying RNA-binding proteins.

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

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

Classifying RNA-binding proteins based on electrostatic properties

Protein structure can provide new insight into the biological function of a protein and can enable the design of better experiments to learn its biological roles. Moreover, deciphering the interactions of a protein with other molecules can contribute to the understanding of the protein's function within cellular processes. In this study, we apply a machine learning approach for classifying RNA-binding proteins based on their three-dimensional structures. The method is based on characterizing unique properties of electrostatic patches on the protein surface. Using an ensemble of general protein features and specific properties extracted from the electrostatic patches, we have trained a support vector machine (SVM) to distinguish RNA-binding proteins from other positively charged proteins that do not bind nucleic acids. Specifically, the method was applied on proteins possessing the RNA recognition motif (RRM) and successfully classified RNA-binding proteins from RRM domains involved in protein–protein interactions. Overall the method achieves 88% accuracy in classifying RNA-binding proteins, yet it cannot distinguish RNA from DNA binding proteins. Nevertheless, by applying a multiclass SVM approach we were able to classify the RNA-binding proteins based on their RNA targets, specifically, whether they bind a ribosomal RNA (rRNA), a transfer RNA (tRNA), or messenger RNA (mRNA). Finally, we present here an innovative approach that does not rely on sequence or structural homology and could be applied to identify novel RNA-binding proteins with unique folds and/or binding motifs.

Gene expression in all living organisms is regulated by a complex set of events at both transcriptional and posttranscriptional levels. RNA-binding proteins play a key role in posttranscriptional events including splicing, stability, transport, and translation. Nowadays, there is increasing evidence that many other cellular processes may be mediated by RNA. Identifying new proteins involved in interaction with RNA is thus essential to unraveling the cellular processes in which these interactions are involved. In the current study we present a successful computational approach for classifying RNA-binding proteins and distinguishing them from other proteins based on structural and electrostatic properties. We test the method on a unique protein domain, the RNA recognition motif (RRM), which mediates both RNA and protein interactions. We show that we can discriminate RNA-binding RRMs from protein-binding RRMs. Further, we demonstrate that we can classify known RNA-binding proteins based on their RNA target (mRNA, rRNA, or tRNA). Our method does not rely on any kind of evolutionary information and thus can be applied to identify RNA-binding proteins with novel modes of RNA recognition.

A new role for HIV nucleocapsid protein in modulating the specificity of plus strand priming

The current study indicates a new role for HIV nucleocapsid protein (NC) in modulating the specificity of plus strand priming. RNase H cleavage by reverse transcriptase (RT) during minus strand synthesis gives rise to RNA fragments that could potentially be used as primers for synthesis of the plus strand, leading to the initiation of priming from multiple points as has been observed for other retroviruses. For HIV, the central and 3' polypurine tracts (PPTs) are the major sites of plus strand initiation. Using reconstituted in vitro assays, results showed that NC greatly reduced the efficiency of extension of non-PPT RNA primers, but not PPT. Experiments mimicking HIV replication showed that RT generated and used both PPT and non-PPT RNAs to initiate "plus strand" synthesis, but non-PPT usage was strongly inhibited by NC. The results support a role for NC in specifying primer usage during plus strand synthesis.

Nucleocapsid protein function in early infection processes

The role of nucleocapsid protein (NC) in the early steps of retroviral replication appears largely that of a facilitator for reverse transcription and integration. Using a wide variety of cell-free assay systems, the properties of mature NC proteins (e.g. HIV-1 p7(NC) or MLV p10(NC)) as nucleic acid chaperones have been extensively investigated. The effect of NC on tRNA annealing, reverse transcription initiation, minus-strand-transfer, processivity of reverse transcription, plus-strand-transfer, strand-displacement synthesis, 3' processing of viral DNA by integrase, and integrase-mediated strand-transfer has been determined by a large number of laboratories. Interestingly, these reactions can all be accomplished to varying degrees in the absence of NC; some are facilitated by both viral and non-viral proteins and peptides that may or may not be involved in vivo. What is one to conclude from the observation that NC is not strictly required for these necessary reactions to occur? NC likely enhances the efficiency of each of these steps, thereby vastly improving the productivity of infection. In other words, one of the major roles of NC is to enhance the effectiveness of early infection, thereby increasing the probability of productive replication and ultimately of retrovirus survival.

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.

Complexes of tRNA and maturation enzymes: shaping up for translation

Several significant structures of transfer ribonucleic acid (tRNA) maturation enzymes complexed with precursor tRNA or fragments thereof have been published recently, providing detailed knowledge of enzyme-tRNA recognition and catalytic strategies. In addition to reinforcing the general principles of RNA-protein interaction, the new structures highlight both the features of composite RNA recognition by multiple enzyme subunits and the pronounced RNA structural flexibility in or near the active site in all cases. These structural principles provide plausible explanations for the exquisite specificity and catalytic power of these enzymes and, in the case of evolutionary adaptation, for the ability of some enzymes to develop novel specificities.

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.

Structure/function mapping of amino acids in the N-terminal zinc finger of the human immunodeficiency virus type 1 nucleocapsid protein: residues responsible for nucleic acid helix destabilizing activity

The nucleocapsid protein (NC) of HIV-1 is 55 amino acids in length and possesses two CCHC-type zinc fingers. Finger one (N-terminal) contributes significantly more to helix destabilizing activity than finger two (C-terminal). Five amino acids differ between the two zinc fingers. To determine at the amino acid level the reason for the apparent distinction between the fingers, each different residue in finger one was incrementally replaced by the one at the corresponding location in finger two. Mutants were analyzed in annealing assays with unstructured and structured substrates. Three groupings emerged: (1) those similar to wild-type levels (N17K, A25M), (2) those with diminished activity (I24Q, N27D), and (3) mutant F16W, which had substantially greater helix destabilizing activity than that of the wild type. Unlike I24Q and the other mutants, N27D was defective in DNA binding. Only I24Q and N27D showed reduced strand transfer in in vitro assays. Double and triple mutants F16W/I24Q, F16W/N27D, and F16W/I24Q/N27D all showed defects in DNA binding, strand transfer, and helix destabilization, suggesting that the I24Q and N27D mutations have a dominant negative effect and abolish the positive influence of F16W. Results show that amino acid differences at positions 24 and 27 contribute significantly to finger one's helix destabilizing activity.

Different modes and potencies of translational repression by sequence-specific RNA-protein interaction at the 5'-UTR

To determine whether sequence-specific RNA–protein interaction at the 5′-untranslated region (5′-UTR) can potently repress translation in mammalian cells, a bicistronic translational repression assay was developed to permit direct assessment of RNA–protein interaction and translational repression in transiently transfected living mammalian cells. Changes in cap-dependent yellow fluorescent protein (YFP) and internal ribosome entry sequence (IRES)-dependent cyan fluorescent protein (CFP) translation were monitored by fluorescence microscopy. Selective repression of YFP or coordinate repression of both YFP and CFP translation occurred, indicating two distinct modes by which RNA-binding proteins repress translation through the 5′-UTR. Interestingly, a single-stranded RNA-binding protein from Bacillus subtilis, tryptophan RNA-binding attenuation protein (TRAP), showed potent translational repression, dependent on the level of TRAP expression and position of its cognate binding site within the bicistronic reporter transcript. As the first of its class to be examined in mammalian cells, its potency in repression of translation through the 5′-UTR may be a general feature for this class of single-stranded RNA-binding proteins. Finally, a one-hybrid screen based on translational repression through the 5′-UTR identified linkers supporting full-translational repression as well as a range of partial repression by TRAP within the context of a fusion protein.

Dissecting the protein-RNA and RNA-RNA interactions in the nucleocapsid-mediated dimerization and isomerization of HIV-1 stemloop 1

The specific binding of HIV-1 nucleocapsid protein (NC) to the different forms assumed in vitro by the stemloop 1 (Lai variant) of the genome's packaging signal has been investigated using electrospray ionization-Fourier transform mass spectrometry (ESI-FTMS). The simultaneous observation of protein-RNA and RNA-RNA interactions in solution has provided direct information about the role of NC in the two-step model of RNA dimerization and isomerization. In particular, two distinct binding sites have been identified on the monomeric stemloop structure, corresponding to the apical loop and stem-bulge motifs. These sites share similar binding affinities that are intermediate between those of stemloop 3 (SL3) and the putative stemloop 4 (SL4) of the packaging signal. Binding to the apical loop, which contains the dimerization initiation site (DIS), competes directly with the annealing of self-complementary sequences to form a metastable kissing-loop (KL) dimer. In contrast, binding to the stem-bulge affects indirectly the monomer-dimer equilibrium by promoting the rearrangement of KL into the more stable extended duplex (ED) conformer. This process is mediated by the duplex-melting activity of NC, which destabilizes the intramolecular base-pairs surrounding the KL stem-bulges and enables their exchange to form the inter-strand pairs that define the ED structure. In this conformer, high-affinity binding takes place at stem-bulge sites that are identical to those present in the monomeric and KL forms. In this case, however, the NC-induced "breathing" does not result in dissociation of the double-stranded structure because of the large number of intermolecular base-pairs. The different binding modes manifested by conformer-specific mutants have shown that NC can also provide low affinity interactions with the bulged-out adenine bases flanking the DIS region of the ED conformer, thus supporting the hypothesis that these exposed nucleotides may constitute "base-grips" for protein contacts during the late stages of the viral lifecycle.

Comprehensive alanine-scanning mutagenesis of Escherichia coli CsrA defines two subdomains of critical functional importance

The RNA-binding protein CsrA (carbon storage regulator) of Escherichia coli is a global regulator of gene expression and is representative of the CsrA/RsmA family of bacterial proteins. These proteins act by regulating mRNA translation and stability and are antagonized by binding to small noncoding RNAs. Although the RNA target sequence and structure for CsrA binding have been well defined, little information exists concerning the protein requirements for RNA recognition. The three-dimensional structures of three CsrA/RsmA proteins were recently solved, revealing a novel protein fold consisting of two interdigitated monomers. Here, we performed comprehensive alanine-scanning mutagenesis on csrA of E. coli and tested the 58 resulting mutants for regulation of glycogen accumulation, motility, and biofilm formation. Quantitative effects of these mutations on expression of glgCA'-'lacZ, flhDC'-'lacZ, and pgaA'-'lacZ translational fusions were also examined, and eight of the mutant proteins were purified and tested for RNA binding. These studies identified two regions of the amino acid sequence that were critical for regulation and RNA binding, located within the first (beta1, residues 2-7) and containing the last (beta5, residues 40-47) beta-strands of CsrA. The beta1 and beta5 strands of opposite monomers lie adjacent and parallel to each other in the three-dimensional structure of this protein. Given the symmetry of the CsrA dimer, these findings imply that two distinct RNA binding surfaces or functional subdomains lie on opposite sides of the protein.

Real-time RNA profiling within a single bacterium

Characterizing the dynamics of specific RNA levels requires real-time RNA profiling in a single cell. We show that the combination of a synthetic modular genetic system with fluorescence correlation spectroscopy allows us to directly measure in real time the activity of any specific promoter in prokaryotes. Using a simple inducible gene expression system, we found that induced RNA levels within a single bacterium of Escherichia coli exhibited a pulsating profile in response to a steady input of inducer. The genetic deletion of an efflux pump system, a key determinant of antibiotic resistance, altered the pulsating transcriptional dynamics and caused overexpression of induced RNA. In contrast with population measurements, real-time RNA profiling permits identifying relationships between genotypes and transcriptional dynamics that are accessible only at the level of the single cell.

Structural basis of single-stranded RNA recognition

RNA is an ancient and highly versatile molecule that plays fundamental roles in all living organisms. Its molecular functions range from being a mediator of genetic information to the regulation of essential cellular processes. These functions are often accomplished in close association with RNA binding proteins. Over the past few years, a considerable number of high-resolution three-dimensional structures of important protein-RNA complexes have been determined. Here, we wish to discuss recent examples and highlight principles and distinct features of single-stranded RNA recognition by conserved RNA binding domains.

An evolutionarily conserved role for SRm160 in 3'-end processing that functions independently of exon junction complex formation

SRm160 (the SR-related nuclear matrix protein of 160 kDa) functions as a splicing coactivator and 3'-end cleavage-stimulatory factor. It is also a component of the splicing-dependent exon-junction complex (EJC), which has been implicated in coupling of pre-mRNA splicing with mRNA turnover and mRNA export. We have investigated whether the association of SRm160 with the EJC is important for efficient 3'-end cleavage. The EJC components RNPS1, REF, UAP56, and Y14 interact with SRm160. However, when these factors were tethered to transcripts, only SRm160 and RNPS1 stimulated 3'-end cleavage. Whereas SRm160 stimulated cleavage to a similar extent in the presence or absence of an active intron, stimulation of 3'-end cleavage by tethered RNPS1 is dependent on an active intron. Assembly of an EJC adjacent to the cleavage and polyadenylation signal in vitro did not significantly affect cleavage efficiency. These results suggest that SRm160 stimulates cleavage independently of its association with EJC components and that the cleavage-stimulatory activity of RNPS1 may be an indirect consequence of its ability to stimulate splicing. Using RNA interference (RNAi) in Caenorhabditis elegans, we determined whether interactions between SRm160 and the cleavage machinery are important in a whole organism context. Simultaneous RNAi of SRm160 and the cleavage factor CstF-50 (Cleavage stimulation factor 50-kDa subunit) resulted in late embryonic developmental arrest. In contrast, RNAi of CstF-50 in combination with RNPS1 or REFs did not result in an apparent phenotype. Our combined results provide evidence for an evolutionarily conserved interaction between SRm160 and the 3'-end cleavage machinery that functions independently of EJC formation.

Differing roles of the N- and C-terminal zinc fingers in human immunodeficiency virus nucleocapsid protein-enhanced nucleic acid annealing

The replication process of human immunodeficiency virus requires a number of nucleic acid annealing steps facilitated by the hybridization and helix-destabilizing activities of human immunodeficiency virus nucleocapsid (NC) protein. NC contains two CCHC zinc finger motifs numbered 1 and 2 from the N terminus. The amino acids surrounding the CCHC residues differ between the two zinc fingers. Assays were preformed to investigate the activities of the fingers by determining the effect of mutant and wild-type proteins on annealing of 42-nucleotide RNA and DNA complements. The mutants 1.1 NC and 2.2 NC had duplications of the N- and C-terminal zinc fingers in positions 1 and 2. The mutant 2.1 NC had the native zinc fingers with their positions switched. Annealing assays were completed with unstructured and highly structured oligonucleotide complements. 2.2 NC had a near wild-type level of annealing of unstructured nucleic acids, whereas it was completely unable to stimulate annealing of highly structured nucleic acids. In contrast, 1.1 NC was able to stimulate annealing of both unstructured and structured substrates, but to a lesser degree than the wild-type protein. Results suggest that finger 1 has a greater role in unfolding of strong secondary structures, whereas finger 2 serves an accessory role that leads to a further increase in the rate of annealing.

Sequence context dependence of tandem guanine:adenine mismatch conformations in RNA: a continuum solvent analysis

Guanine:adenine (G:A) mismatches and in particular tandem G:A (tG:A) mismatches are frequently observed in biological RNA molecules and can serve as sites for tertiary interaction, metal binding and protein recognition. Depending on the surrounding sequence tG:A mismatches can adopt different basepairing topologies. In the sequence context (5'-) GGAC (tandem G:A in bold) a face-to-face (imino or Watson-Crick-like) pairing is preferred whereas in the CGAG context, G and A adopt a sheared arrangement. Systematic conformational searches with a generalized Born continuum model and molecular dynamics simulations including explicit water molecules and ions have been used to generate face-to-face and sheared tG:A mismatches in both CGAG and GGAC sequence contexts. Conformations from both approaches were evaluated using the same force field and a Poisson-Boltzmann continuum solvent model. Although the substate analysis predicted the sheared arrangement to be energetically preferred in both sequence contexts, a significantly greater preference of the sheared form was found for the CGAG context. In agreement with the experimental observation, the analysis of molecular dynamics trajectories indicated a preference of the sheared form in the case of the CGAG-context and a favorization of the face-to-face form in the case of the GGAC context. The computational studies allowed to identify energetic contributions that stabilize or destabilize the face-to-face and sheared tandem mismatch topologies. The calculated nonpolar solvation and Lennard-Jones packing interaction were found to stabilize the sheared topology independent of the sequence context. Electrostatic contributions are predicted to make the most significant contribution to the sequence context dependence on the structural preference of tG:A mismatches.

Crystal structure of the human CCA-adding enzyme: insights into template-independent polymerization

All tRNA molecules carry the invariant sequence CCA at their 3'-terminus for amino acid attachment. The post-transcriptional addition of CCA is carried out by ATP(CTP):tRNA nucleotidyltransferase, also called CCase. This enzyme catalyses a unique template-independent but sequence-specific nucleotide polymerization reaction. In order to reveal the molecular mechanism of this activity, we solved the crystal structure of human CCase by single isomorphous replacement. The structure reveals a four domain architecture with a cluster of conserved residues forming a positively charged cleft between the first two domains. Structural homology of the N-terminal CCase domain to other nucleotidyltransferases could be exploited for modeling a tRNA-substrate complex. The model places the tRNA 3'-end into the N-terminal nucleotidyltransferase site, close to a patch of conserved residues that provide the binding sites for CTP and ATP. Based on our results, we introduce a corkscrew model for CCA addition that includes a fixed active site and a traveling tRNA-binding region formed by flexible parts of the protein.

Structure and function of the PWI motif: a novel nucleic acid-binding domain that facilitates pre-mRNA processing

The PWI motif is a highly conserved domain of unknown function in the SRm160 splicing and 3'-end cleavage-stimulatory factor, as well as in several other known or putative pre-mRNA processing components. We show here that the PWI motif is a new type of RNA/DNA-binding domain that has an equal preference for single- and double-stranded nucleic acids. Deletion of the motif prevents SRm160 from binding RNA and stimulating 3'-end cleavage, and its substitution with a heterologous RNA-binding domain restores these functions. The NMR solution structure of the SRm160-PWI motif reveals a novel, four-helix bundle and represents the first example of an alpha-helical fold that can bind single-stranded (ss)RNA. Structure-guided mutagenesis indicates that the same surface is involved in RNA and DNA binding and requires the cooperative action of a highly conserved, adjacent basic region. Thus, the PWI motif is a novel type of nucleic acid-binding domain that likely has multiple important functions in pre-mRNA processing, including SRm160-dependent stimulation of 3'-end formation.

Protein-RNA interactions and virus stability as probed by the dynamics of tryptophan side chains

The correlation between dynamics and stability of icosahedral viruses was studied by steady-state and time-resolved fluorescence approaches. We compared the environment and dynamics of tryptophan side chains of empty capsids and ribonucleoprotein particles of two icosahedral viruses from the comovirus group: cowpea mosaic virus (CPMV) and bean pod mottle virus (BPMV). We found a great difference between tryptophan fluorescence emission spectra of the ribonucleoprotein particles and the empty capsids of BPMV. For CPMV, time-resolved fluorescence revealed differences in the tryptophan environments of the capsid protein. The excited-state lifetimes of tryptophan residues were significantly modified by the presence of RNA in the capsid. More than half of the emission of the tryptophans in the ribonucleoprotein particles of CPMV originates from a single exponential decay that can be explained by a similar, nonpolar environment in the local structure of most of the tryptophans, even though they are physically located in different regions of the x-ray structure. CPMV particles without RNA lost this discrete component of emission. Anisotropy decay measurements demonstrated that tryptophans rotate faster in empty particles when compared with the ribonucleoprotein particles. The increased structural breathing facilitates the denaturation of the empty particles. Our studies bring new insights into the intricate interactions between protein and RNA where part of the missing structural information on the nucleic acid molecule is compensated for by the dynamics.

Automated selection of aptamers against protein targets translated in vitro: from gene to aptamer

Reagents for proteome research must of necessity be generated by high throughput methods. Aptamers are potentially useful as reagents to identify and quantitate individual proteins, yet are currently produced for the most part by manual selection procedures. We have developed automated selection methods, but must still individually purify protein targets. Therefore, we have attempted to select aptamers against protein targets generated by in vitro transcription and translation of individual genes. In order to specifically immobilize the protein targets for selection, they are also biotinylated in vitro. As a proof of this method, we have selected aptamers against translated human U1A, a component of the nuclear spliceosome. Selected sequences demonstrated exquisite mimicry of natural binding sequences and structures. These results not only reveal a potential path to the high throughput generation of aptamers, but also yield insights into the incredible specificity of the U1A protein for its natural RNA ligands.

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

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

HIV-1 nucleocapsid protein as a nucleic acid chaperone: spectroscopic study of its helix-destabilizing properties, structural binding specificity, and annealing activity

Assembly of infectious retroviral particles involves recognition of specific sequences on the viral RNA by the nucleocapsid (NC) domain of the Gag polyprotein, and subsequent stoichiometric binding of the processed NC protein along the entire length of the RNA. NC proteins also act as nucleic acid chaperones. They accelerate nucleic acid hybridization and strand exchange, which may be critical during the initial stages of reverse transcription. In order to better understand these properties, we have studied the nucleic acid helix-destabilizing t(m)-depressing) and binding activities of HIV-1 NCp7 protein with a variety of substrates, and the real-time kinetics of NC-induced strand exchange. At low ionic strength (0.01 M Na phosphate, pH 7.0) and saturating levels of protein, NCp7 displays moderate helix-destabilizing activity on double-stranded DNA. Saturating levels of NCp7 lowered the t(m) of a synthetic 28 base-pair 28(+)/28(-) oligonucleotide duplex by about 10 deg. C (51 to 41 degrees C). The presence of single-stranded calf thymus DNA (equimolar with duplex) eliminated the t(m) depression, whereas double-stranded calf thymus DNA only altered the t(m) of the 28-mer duplex by about 2 deg. C. Similar effects were seen with duplexes with single-stranded overhangs or internal single-stranded gaps. Binding experiments utilizing intrinsic tryptophan quenching indicated significant affinity (K(d) about 0.1 microM) for both single-stranded and double-stranded forms of the 28-mer in 0.01 M sodium phosphate at 25 degrees C, although long-chain (calf thymus double-stranded) DNA displayed a much lower affinity. The effects of NCp7 on the kinetics of nucleic acid annealing, strand exchange, and strand displacement were determined by use of oligonucleotides with end-labeled fluorophores serving as donor-acceptor pairs. NCp7 accelerated all these reactions. In the strand exchange reaction, an imperfect duplex, 28(+)/21(-), was reacted with a perfect complement, 28(-). The kinetics of 28(+)/28(-) annealing in this reaction did not conform to a simple bimolecular model, but could be well fit to the sum of two exponential decays. Addition of stoichiometric levels of NCp7 increased the rate constants of both components, and significantly increased the fraction of exchange associated with the rapid process. Increasing levels of 28(-) also increased the rapid fraction, as well as the rapid rate constant. This concentration dependence indicates that, although the kinetic decays appear biexponential, at least one of the steps is bimolecular. Simple annealing reactions, 28(+) with 28(-), could be fit to single-exponential decays, and their magnitudes in the presence of NCp7 were comparable to the rapid step of annealing observed for exchange reactions, suggesting that this step is connected with annealing. Strand dissociation during exchange was monitored by placing the fluorescent acceptor on the 21(-) strand. The results, though complex, suggest that the slow step of exchange is largely associated with the dissociation of the shorter oligonucleotide. Analogous experiments were performed with variants of these oligonucleotides, and the results are in line with the 28(+)/21(-)/28(-) experiments. On the basis of an analysis of the effect of increasing levels of 28(-) on the formation of the perfect 28 bp duplex from the imperfect duplex, we propose that NCp7 forms a ternary complex intermediate with imperfect duplex and 28(-), and suggest several ways by which such an intermediate would facilitate strand exchange.

Modelling ion binding to AA platform motifs in RNA: a continuum solvent study including conformational adaptation

Binding of monovalent and divalent cations to two adenine-adenine platform structures from the Tetrahymena group I intron ribozyme has been studied using continuum solvent models based on the generalised Born and the finite-difference Poisson-Boltzmann approaches. The adenine-adenine platform RNA motif forms an experimentally characterised monovalent ion binding site important for ribozyme folding and function. Qualitative agreement between calculated and experimental ion placements and binding selectivity was obtained. The inclusion of solvation effects turned out to be important to obtain low energy structures and ion binding placements in agreement with the experiment. The calculations indicate that differences in solvation of the isolated ions contribute to the calculated ion binding preference. However, Coulomb attraction and van der Waals interactions due to ion size differences and RNA conformational adaptation also influence the calculated ion binding affinity. The calculated alkali ion binding selectivity for both platforms followed the order K(+) > Na(+) > Rb(+) > Cs(+) > Li(+) (Eisenman series VI) in the case of allowing RNA conformational relaxation during docking. With rigid RNA an Eisenman series V was obtained (K(+) > Rb(+) > Na(+) > Cs(+) > Li(+)). Systematic energy minimisation docking simulations starting from several hundred initial placements of potassium ions on the surface of platform containing RNA fragments identified a coordination geometry in agreement with the experiment as the lowest energy binding site. The approach could be helpful to identify putative ion binding sites in nucleic acid structures determined at low resolution or with experimental methods that do not allow identification of ion binding sites.

Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and protein-nucleic acid noncovalent interactions

Computer modeling has been developed and widely applied in studying molecules of biological interest. The force field is the cornerstone of computer simulations, and many force fields have been developed and successfully applied in these simulations. Two interesting areas are (a) studying enzyme catalytic mechanisms using a combination of quantum mechanics and molecular mechanics, and (b) studying macromolecular dynamics and interactions using molecular dynamics (MD) and free energy (FE) calculation methods. Enzyme catalysis involves forming and breaking of covalent bonds and requires the use of quantum mechanics. Noncovalent interactions appear ubiquitously in biology, but here we confine ourselves to review only noncovalent interactions between protein and protein, protein and ligand, and protein and nucleic acids.

Crystal structure of a heptameric Sm-like protein complex from archaea: implications for the structure and evolution of snRNPs

The Sm/Lsm proteins associate with small nuclear RNA to form the core of small nuclear ribonucleoproteins, required for processes as diverse as pre-mRNA splicing, mRNA degradation and telomere formation. The Lsm proteins from archaea are likely to represent the ancestral Sm/Lsm domain. Here, we present the crystal structure of the Lsm alpha protein from the thermophilic archaeon Methanobacterium thermoautotrophicum at 2.0 A resolution. The Lsm alpha protein crystallizes as a heptameric ring comprised of seven identical subunits interacting via beta-strand pairing and hydrophobic interactions. The heptamer can be viewed as a propeller-like structure in which each blade consists of a seven-stranded antiparallel beta-sheet formed from neighbouring subunits. There are seven slots on the inner surface of the heptamer ring, each of which is lined by Asp, Asn and Arg residues that are highly conserved in the Sm/Lsm sequences. These conserved slots are likely to form the RNA-binding site. In archaea, the gene encoding Lsm alpha is located next to the L37e ribosomal protein gene in a putative operon, suggesting a role for the Lsm alpha complex in ribosome function or biogenesis.

A comprehensive characterization of a group IB intron and its encoded maturase reveals that protein-assisted splicing requires an almost intact intron RNA

The group I intron (AnCOB) of the mitochondrial apocytochrome b gene from Aspergillus nidulans encodes a bi-functional maturase protein that is also a DNA endonuclease. Although the AnCOB intron self-splices, the encoded maturase protein greatly facilitates splicing, in part, by stabilizing RNA tertiary structure. To determine their role in self-splicing and in protein-assisted splicing, several peripheral RNA sub-domains in the 313 nucleotide intron were deleted (P2, P9, P9.1) or truncated (P5ab, P6a). The sequence in two helices (P2 and P9) was also inverted. Except for P9, the deleted regions are not highly conserved among group I introns and are often dispensable for catalytic activity. Nevertheless, despite the very tight binding of AnCOB RNA to the maturase and the high activity of the bimolecular complex (the rate of 5' splice-site cleavage was >20 min(-1) with guanosine as the cofactor), the intron was surprisingly sensitive to these modifications. Several mutations inactivated splicing completely and virtually all impaired splicing to varying degrees. Mutants containing comparatively small deletions in various regions of the intron significantly decreased binding affinity (generally >10(4)-fold), indicating that none of the domains that remained constitutes the primary recognition site of the maturase. The data argue that tight binding requires tertiary interactions that can be maintained by only a relatively intact intron RNA, and that the binding mechanism of the maturase differs from those of two other well-characterized group I intron splicing factors, CYT-18 and Cpb2. A model is proposed in which the protein promotes widespread cooperative folding of an RNA lacking extensive initial tertiary structure.

Protein-RNA interactions: a structural analysis

A detailed computational analysis of 32 protein-RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein-double-stranded DNA and protein-single-stranded DNA complexes. The interface properties of the protein-RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein-RNA and protein-DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein-RNA complexes, backbone contacts were more dominant in the protein-DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Binding of a cyclic BIV beta-Tat peptide with its TAR RNA construct

The ability of RNA structures to adopt diverse yet complex tertiary structures has resulted in numerous fascinating RNA-protein recognition events. It was recently reported that a close relative of the HIV Rev peptide, namely a 17 residue Tat peptide from bovine immuno-deficiency virus (BIV), is able to bind to the 28 nucleotide BIV TAR RNA construct. Here we report that by simply converting the 17 residue beta-ribbon peptide structure to a 19 residue cyclopeptide, the binding affinity (Kd) of the resulting cyclopeptide to the TAR RNA target, observed by fluorescence binding study, was enhanced approximately 5-fold.

Transfer RNA recognition by aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases are an ancient group of enzymes that catalyze the covalent attachment of an amino acid to its cognate transfer RNA. The question of specificity, that is, how each synthetase selects the correct individual or isoacceptor set of tRNAs for each amino acid, has been referred to as the second genetic code. A wealth of structural, biochemical, and genetic data on this subject has accumulated over the past 40 years. Although there are now crystal structures of sixteen of the twenty synthetases from various species, there are only a few high resolution structures of synthetases complexed with cognate tRNAs. Here we review briefly the structural information available for synthetases, and focus on the structural features of tRNA that may be used for recognition. Finally, we explore in detail the insights into specific recognition gained from classical and atomic group mutagenesis experiments performed with tRNAs, tRNA fragments, and small RNAs mimicking portions of tRNAs.

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

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

Stitching together RNA tertiary architectures

The powerful explanatory paradigm of molecular biology requiring form to co-evolve with function has again been proven successful when, over the recent two decades, a wealth of biological functions have been uncovered for RNA. Previously considered as a mere mediator of the genetic code, RNA is now acknowledged as a key player in a wide variety of cellular processes. Along with the discovery of novel biological functions of RNA molecules, a number of RNA three-dimensional structures have been solved which beautifully demonstrate the molecular adaptability which allows RNA to participate as a key player in these functions. A distinct repertoire of molecular motifs provides a basis for the assembly of complex RNA tertiary architectures.