Themes in RNA-protein recognition

Exploring the role of cation–π interactions in glycoproteins lipid-binding proteins and RNA-binding proteins

We have analyzed and compared the influence of cation-pi interactions in glycoproteins (GPs), lipid-binding proteins (LBPs) and RNA-binding proteins (RBPs) in this study. We observed that all the proteins included in the study had profound cation-pi interactions. There is an average of one energetically significant cation-pi interaction for every 71 residues in GPs, for every 58 residues in LBPs and for every 64 residues in RBPs. Long-range contacts are predominant in all the three types of proteins studied. The pair-wise cation-pi interaction energy between the positively charged and aromatic residues shows that Arg-Trp pair energy was the strongest among all six possible pairs in all the three types of proteins studied. There were considerable differences in the preference of cation-pi interacting residues to different secondary structure elements and ASA and these might contribute to differences in biochemical functions of GPs, LBPs and RBPs. It was interesting to note that all the five residues involved in cation-pi interactions were found to have stabilization centers in GPs, LBPs and RBPs. Majority of the cation-pi interacting residues investigated in the present study had a conservation score of 6, the cutoff value used to identify the stabilizing residues. A small percentage of cation-pi interacting residues were also present as stabilizing residues. The cation-pi interaction-forming residues play an important role in the structural stability of in GPs, LBPs and RBPs. The results obtained in this study will be helpful in further understanding the stability, specificity and differences in the biochemical functions of GPs, LBPs and RBPs.

Staufen1 regulates diverse classes of mammalian transcripts

It is currently unknown how extensively the double-stranded RNA-binding protein Staufen (Stau)1 is utilized by mammalian cells to regulate gene expression. To date, Stau1 binding to the 3'-untranslated region (3'-UTR) of ADP ribosylation factor (ARF)1 mRNA has been shown to target ARF1 mRNA for Stau1-mediated mRNA decay (SMD). ARF1 SMD depends on translation and recruitment of the nonsense-mediated mRNA decay factor Upf1 to the ARF1 3'-UTR by Stau1. Here, we demonstrate that Stau1 binds to a complex structure within the ARF1 3'-UTR. We also use microarrays to show that 1.1 and 1.0% of the 11 569 HeLa-cell transcripts that were analyzed are upregulated and downregulated, respectively, at least two-fold upon Stau1 depletion in three independently performed experiments. We localize the Stau1 binding site to the 3'-UTR of four mRNAs that we define as natural SMD targets. Additionally, we provide evidence that the efficiency of SMD increases during the differentiation of C2C12 myoblasts to myotubes. We propose that Stau1 influences the expression of a wide variety of physiologic transcripts and metabolic pathways.

Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.

RNABindR: a server for analyzing and predicting RNA-binding sites in proteins

Understanding interactions between proteins and RNA is key to deciphering the mechanisms of many important biological processes. Here we describe RNABindR, a web-based server that identifies and displays RNA-binding residues in known protein–RNA complexes and predicts RNA-binding residues in proteins of unknown structure. RNABindR uses a distance cutoff to identify which amino acids contact RNA in solved complex structures (from the Protein Data Bank) and provides a labeled amino acid sequence and a Jmol graphical viewer in which RNA-binding residues are displayed in the context of the three-dimensional structure. Alternatively, RNABindR can use a Naive Bayes classifier trained on a non-redundant set of protein–RNA complexes from the PDB to predict which amino acids in a protein sequence of unknown structure are most likely to bind RNA. RNABindR automatically displays ‘high specificity’ and ‘high sensitivity’ predictions of RNA-binding residues. RNABindR is freely available at http://bindr.gdcb.iastate.edu/RNABindR.

Characterization of the RNA-binding regions in protein p36 of Heliothis armigera cypovirus 14

Some proteins of cypovirus (CPV) bind to RNA, probably contributing to the replication of viral genome. However, little is known about whether any protein from Heliothis armigera cypovirus (HaCPV) could bind to RNA. In this study, we cloned the ORF of segment 9 (S9) of HaCPV, serotype 14, into pMAL-c2X for the generation and purification of maltose binding protein (MBP) fused protein p36 (MBP-p36). The analysis of the RNA-binding properties of MBP-p36 revealed that p36, but not MBP alone, bound to ssRNA of CPV. Furthermore, the ssRNA-binding activities of p36 were significantly inhibited or completely eliminated by protein denaturants or unsuitable concentrations of NaCl. Importantly, the formation of ssRNA/p36 was only competitively inhibited by a heavy dose of competitive non-viral ssRNA or dsRNA, but not by ssDNA and dsDNA, suggesting that p36 bound to both ssRNA and dsRNA, but not DNA. Moreover, the characterization of different mutants of p36 revealed that the regions 1-26aa, 154-170aa, and 229-238aa, but not region 291-320aa, may be crucial for the ssRNA-binding ability of p36. Conceivably, the sensitivity of p36 to denaturants and the synergetic effect of different regions suggest that the RNA-binding ability of p36 may be conformation-dependent. Thus, our findings provide new insights into understanding the genomic function of HaCPV-14.

The mRNA-binding site of annexin A2 resides in helices C-D of its domain IV

Annexin A2 (AnxA2) is a Ca(2+)-binding and phospholipid-binding protein involved in different intracellular processes including exocytosis, endocytosis and membrane-cytoskeleton movements. We have previously identified AnxA2 as an mRNA-binding protein present in cytoskeleton-bound polysomes, that binds to a specific approximately 100 nucleotide region in the 3'-untranslated region of c-myc and its cognate mRNA. In the present study, we show by UV cross-linking assays and surface plasmon resonance analyses that the mRNA-binding site of AnxA2 resides in its domain IV. Furthermore, the interaction of full-length AnxA2 with the 3'-untranslated region of anxA2 mRNA is Ca(2+)-dependent. By contrast, the interaction is Ca(2+)-independent for the isolated domain IV of AnxA2, suggesting that the mRNA-binding site is masked in Apo-AnxA2 and gains exposure through a Ca(2+)-induced conformational change of AnxA2 generating a favourable mRNA-binding site. The AnxA2-mRNA interaction is specific and involves helices C and D in domain IV of AnxA2, since point mutagenesis of several charged and polar exposed residues of these helices in the full-length protein strongly reduce RNA binding. The interaction appears to be sequential involving an initial phase of recognition dominated by electrostatic interactions, most likely between lysine residues and the phosphate backbone of RNA, followed by a second phase contributing to the specificity of the interaction.

The genetic code is nearly optimal for allowing additional information within protein-coding sequences

DNA sequences that code for proteins need to convey, in addition to the protein-coding information, several different signals at the same time. These "parallel codes" include binding sequences for regulatory and structural proteins, signals for splicing, and RNA secondary structure. Here, we show that the universal genetic code can efficiently carry arbitrary parallel codes much better than the vast majority of other possible genetic codes. This property is related to the identity of the stop codons. We find that the ability to support parallel codes is strongly tied to another useful property of the genetic code--minimization of the effects of frame-shift translation errors. Whereas many of the known regulatory codes reside in nontranslated regions of the genome, the present findings suggest that protein-coding regions can readily carry abundant additional information.

Cucumber mosaic virus as a presentation system for a double hepatitis C virus-derived epitope

Chimeric plant viruses are emerging as promising vectors for use in innovative vaccination strategies. In this context, cucumber mosaic virus (CMV) has proven to be a suitable carrier of the hepatitis C virus (HCV)-derived R9 mimotope. In the present work, a new chimeric CMV, expressing on its surface the HCV-derived R10 mimotope, was produced but lost the insert after the first passage on tobacco. A comparative analysis between R10- and R9-CMV properties indicated that R9-CMV stability was related to structural features typical of the foreign insert. Thus, in order to combine high virus viability with strong immuno-stimulating activity, we doubled R9 copies on each of the 180 coat protein (CP) subunits of CMV. One of the chimeras produced by this approach (2R9-CMV) was shown to systemically infect the host, stably maintaining both inserts. Notably, it was strongly recognized by sera of HCV-infected patients and, as compared with R9-CMV, displayed an enhanced ability to stimulate lymphocyte IFN-gamma production. The high immunogen levels achievable in plants or fruits infected with 2R9-CMV suggest that this chimeric form of CMV may be useful in the development of oral vaccines against HCV.

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 of dual function iron regulatory protein 1 complexed with ferritin IRE-RNA

Iron regulatory protein 1 (IRP1) binds iron-responsive elements (IREs) in messenger RNAs (mRNAs), to repress translation or degradation, or binds an iron-sulfur cluster, to become a cytosolic aconitase enzyme. The 2.8 angstrom resolution crystal structure of the IRP1:ferritin H IRE complex shows an open protein conformation compared with that of cytosolic aconitase. The extended, L-shaped IRP1 molecule embraces the IRE stem-loop through interactions at two sites separated by approximately 30 angstroms, each involving about a dozen protein:RNA bonds. Extensive conformational changes related to binding the IRE or an iron-sulfur cluster explain the alternate functions of IRP1 as an mRNA regulator or enzyme.

Amino acid residue doublet propensity in the protein-RNA interface and its application to RNA interface prediction

Protein-RNA interactions play essential roles in a number of regulatory mechanisms for gene expression such as RNA splicing, transport, translation and post-transcriptional control. As the number of available protein-RNA complex 3D structures has increased, it is now possible to statistically examine protein-RNA interactions based on 3D structures. We performed computational analyses of 86 representative protein-RNA complexes retrieved from the Protein Data Bank. Interface residue propensity, a measure of the relative importance of different amino acid residues in the RNA interface, was calculated for each amino acid residue type (residue singlet interface propensity). In addition to the residue singlet propensity, we introduce a new residue-based propensity, which gives a measure of residue pairing preferences in the RNA interface of a protein (residue doublet interface propensity). The residue doublet interface propensity contains much more information than the sum of two singlet propensities alone. The prediction of the RNA interface using the two types of propensities plus a position-specific multiple sequence profile can achieve a specificity of about 80%. The prediction method was then applied to the 3D structure of two mRNA export factors, TAP (Mex67) and UAP56 (Sub2). The prediction enables us to point out candidate RNA interfaces, part of which are consistent with previous experimental studies and may contribute to elucidation of atomic mechanisms of mRNA export.

Structure and dynamics of DNA-dendrimer complexation: role of counterions, water, and base pair sequence

We study sequence-dependent complexation between oligonucleotides (single-strand DNA) and various generation ethylene diamine (EDA) cored poly amido amide (PAMAM) dendrimers through atomistic molecular dynamics simulations accompanied by free energy calculations and inherent structure determination. Simulations reveal formation of a stable complex and provide a detailed molecular level understanding of the structure and dynamics of such a complexation. The reaction free energy surface in the initial stage is found to be funnel-like, with a significant barrier arising in the late stage due to the occurrence of misfolded states of DNA. Complexation shows surprisingly strong sensitivity to the ssDNA sequence, which is found to arise from a competition between enthalpic versus entropic rigidity of ssDNA.

Stability of the 'L12 stalk' in ribosomes from mesophilic and (hyper)thermophilic Archaea and Bacteria

The ribosomal stalk complex, consisting of one molecule of L10 and four or six molecules of L12, is attached to 23S rRNA via protein L10. This complex forms the so-called ‘L12 stalk’ on the 50S ribosomal subunit. Ribosomal protein L11 binds to the same region of 23S rRNA and is located at the base of the ‘L12 stalk’. The ‘L12 stalk’ plays a key role in the interaction of the ribosome with translation factors. In this study stalk complexes from mesophilic and (hyper)thermophilic species of the archaeal genus Methanococcus and from the Archaeon Sulfolobus solfataricus, as well as from the Bacteria Escherichia coli, Geobacillus stearothermophilus and Thermus thermophilus, were overproduced in E.coli and purified under non-denaturing conditions. Using filter-binding assays the affinities of the archaeal and bacterial complexes to their specific 23S rRNA target site were analyzed at different pH, ionic strength and temperature. Affinities of both archaeal and bacterial complexes for 23S rRNA vary by more than two orders of magnitude, correlating very well with the growth temperatures of the organisms. A cooperative effect of binding to 23S rRNA of protein L11 and the L10/L12₄ complex from mesophilic and thermophilic Archaea was shown to be temperature-dependent.

MIMOX: a web tool for phage display based epitope mapping

Background

Phage display is widely used in basic research such as the exploration of protein-protein interaction sites and networks, and applied research such as the development of new drugs, vaccines, and diagnostics. It has also become a promising method for epitope mapping. Research on new algorithms that assist and automate phage display based epitope mapping has attracted many groups. Most of the existing tools have not been implemented as an online service until now however, making it less convenient for the community to access, utilize, and evaluate them.

Results

We present MIMOX, a free web tool that helps to map the native epitope of an antibody based on one or more user supplied mimotopes and the antigen structure. MIMOX was coded in Perl using modules from the Bioperl project. It has two sections. In the first section, MIMOX provides a simple interface for ClustalW to align a set of mimotopes. It also provides a simple statistical method to derive the consensus sequence and embeds JalView as a Java applet to view and manage the alignment. In the second section, MIMOX can map a single mimotope or a consensus sequence of a set of mimotopes, on to the corresponding antigen structure and search for all of the clusters of residues that could represent the native epitope. NACCESS is used to evaluate the surface accessibility of the candidate clusters; and Jmol is embedded to view them interactively in their 3D context. Initial case studies show that MIMOX can reproduce mappings from existing tools such as FINDMAP and 3DEX, as well as providing novel, rational results.

Conclusion

A web-based tool called MIMOX has been developed for phage display based epitope mapping. As a publicly available online service in this area, it is convenient for the community to access, utilize, and evaluate, complementing other existing programs. MIMOX is freely available at .

Protein-RNA interactions: exploring binding patterns with a three-dimensional superposition analysis of high resolution structures

MOTIVATION

The recognition of specific RNA sequences and structures by proteins is critical to our understanding of RNA processing, gene expression and viral replication. The diversity of RNA structures suggests that RNA recognition is substantially different than that of DNA.

RESULTS

The atomic coordinates of 41 protein-RNA complexes have been used to probe composite nucleoside binding pockets that form the structural and chemical underpinnings of base recognition. Composite nucleoside binding pockets were constructed using three-dimensional superpositions of each RNA nucleoside. Unlike protein-DNA interactions which are dominated by accessibility, RNA recognition frequently occurs in non-canonical and single-strand-like structures that allow interactions to occur from a much wider set of geometries and make fuller use of unique base shapes and hydrogen-bonding ability. By constructing composites that include all van der Waals, hydrogen-bonding, stacking and general non-polar interactions made to a particular nucleoside, the strategies employed are made readily visible. Protein-RNA interactions can result in the formation of a glove-like tight binding pocket around RNA bases, but the size, shape and non-polar binding patterns differ between specific RNA bases. We show that adenine can be distinguished from guanine based on the size and shape of the binding pocket and steric exclusion of the guanine N2 exocyclic amino group. The unique shape and hydrogen-bonding pattern for each RNA base allow proteins to make specific interactions through a very small number of contacts, as few as two in some cases.

AVAILABILITY

The program ENTANGLE is available from http://www.bioc.rice.edu/~shamoo

BindN: a web-based tool for efficient prediction of DNA and RNA binding sites in amino acid sequences

BindN (http://bioinformatics.ksu.edu/bindn/) takes an amino acid sequence as input and predicts potential DNA or RNA-binding residues with support vector machines (SVMs). Protein datasets with known DNA or RNA-binding residues were selected from the Protein Data Bank (PDB), and SVM models were constructed using data instances encoded with three sequence features, including the side chain pK(a) value, hydrophobicity index and molecular mass of an amino acid. The results suggest that DNA-binding residues can be predicted at 69.40% sensitivity and 70.47% specificity, while prediction of RNA-binding residues achieves 66.28% sensitivity and 69.84% specificity. When compared with previous studies, the SVM models appear to be more accurate and more efficient for online predictions. BindN provides a useful tool for understanding the function of DNA and RNA-binding proteins based on primary sequence data.

Prediction of RNA binding sites in proteins from amino acid sequence

RNA-protein interactions are vitally important in a wide range of biological processes, including regulation of gene expression, protein synthesis, and replication and assembly of many viruses. We have developed a computational tool for predicting which amino acids of an RNA binding protein participate in RNA-protein interactions, using only the protein sequence as input. RNABindR was developed using machine learning on a validated nonredundant data set of interfaces from known RNA-protein complexes in the Protein Data Bank. It generates a classifier that captures primary sequence signals sufficient for predicting which amino acids in a given protein are located in the RNA-protein interface. In leave-one-out cross-validation experiments, RNABindR identifies interface residues with >85% overall accuracy. It can be calibrated by the user to obtain either high specificity or high sensitivity for interface residues. RNABindR, implementing a Naive Bayes classifier, performs as well as a more complex neural network classifier (to our knowledge, the only previously published sequence-based method for RNA binding site prediction) and offers the advantages of speed, simplicity and interpretability of results. RNABindR predictions on the human telomerase protein hTERT are in good agreement with experimental data. The availability of computational tools for predicting which residues in an RNA binding protein are likely to contact RNA should facilitate design of experiments to directly test RNA binding function and contribute to our understanding of the diversity, mechanisms, and regulation of RNA-protein complexes in biological systems. (RNABindR is available as a Web tool from http://bindr.gdcb.iastate.edu.).

The role of RNA structure in the interaction of U1A protein with U1 hairpin II RNA

The N-terminal RNA Recognition Motif (RRM1) of the spliceosomal protein U1A interacting with its target U1 hairpin II (U1hpII) has been used as a paradigm for RRM-containing proteins interacting with their RNA targets. U1A binds to U1hpII via direct interactions with a 7-nucleotide (nt) consensus binding sequence at the 5' end of a 10-nt loop, and via hydrogen bonds with the closing C-G base pair at the top of the RNA stem. Using surface plasmon resonance (Biacore), we have examined the role of structural features of U1hpII in binding to U1A RRM1. Mutational analysis of the closing base pair suggests it plays a minor role in binding and mainly prevents "breathing" of the loop. Lengthening the stem and nontarget part of the loop suggests that the increased negative charge of the RNA might slightly aid association. However, this is offset by an increase in dissociation, which may be caused by attraction of the RRM to nontarget parts of the RNA. Studies of a single stranded target and RNAs with untethered loops indicate that structure is not very relevant for association but is important for complex stability. In particular, breaking the link between the stem and the 5' side of the loop greatly increases complex dissociation, presumably by hindering simultaneous contacts between the RRM and stem and loop nucleotides. While binding of U1A to a single stranded target is much weaker than to U1hpII, it occurs with nanomolar affinity, supporting recent evidence that binding of unstructured RNA by U1A has physiological significance.

PROFEAT: a web server for computing structural and physicochemical features of proteins and peptides from amino acid sequence

Sequence-derived structural and physicochemical features have frequently been used in the development of statistical learning models for predicting proteins and peptides of different structural, functional and interaction profiles. PROFEAT (Protein Features) is a web server for computing commonly-used structural and physicochemical features of proteins and peptides from amino acid sequence. It computes six feature groups composed of ten features that include 51 descriptors and 1447 descriptor values. The computed features include amino acid composition, dipeptide composition, normalized Moreau-Broto autocorrelation, Moran autocorrelation, Geary autocorrelation, sequence-order-coupling number, quasi-sequence-order descriptors and the composition, transition and distribution of various structural and physicochemical properties. In addition, it can also compute previous autocorrelations descriptors based on user-defined properties. Our computational algorithms were extensively tested and the computed protein features have been used in a number of published works for predicting proteins of functional classes, protein-protein interactions and MHC-binding peptides. PROFEAT is accessible at http://jing.cz3.nus.edu.sg/cgi-bin/prof/prof.cgi.

Escherichia coli Hfq binds A18 and DsrA domain II with similar 2:1 Hfq6/RNA stoichiometry using different surface sites

Hfq is a RNA-binding protein in Escherichia coli that plays an essential role in post-transcriptional regulation of mRNAs by facilitating pairing of noncoding RNAs (ncRNAs) to mRNA target sites. Recent work has provided evidence that E. coli Hfq has two distinct RNA-binding surfaces. In this study, a comparative sequence-structure analysis of hfq genes in bacterial genomes was employed to identify conserved residues that may be involved in binding RNA. A covariance of residue properties at neighboring positions 12 and 39 and conserved surface residues with high propensities at binding sites of RNA-binding proteins suggested several sites for Hfq-RNA interactions. On the basis of these predictions, eight mutant Hfq proteins were produced and their interactions were examined with the 38 nucleotide (nt) domain II of DsrA ncRNA (DsrA(DII)) and A(18) by a gel-mobility shift assay, fluorescence anisotropy, and fluorescence quenching. Mutations on the proximal surface of Hfq had a small affect on Hfq binding to A(18) (<or=2-fold), while the mutations Y25A and K31A on the distal surface decreased affinity to A(18) by 100-fold in solution. Mutations F39A and R16A on the proximal surface reduced affinity to DsrA(DII) by 6-8-fold, while other mutations on the distal or proximal surfaces affected affinity to DsrA(DII) by <or=2-fold using the gel-mobility shift assay. The F39A/L12F double mutation partially regained the affinity for DsrA(DII) lost by the F39A mutation. The latter observation is consistent with the implied importance of an aromatic residue at position 12 or 39 suggested by the sequence covariance. Titration experiments indicate a 2:1 Hfq(6)/RNA stoichiometry for the strong binding complexes of Hfq with either A(18) or DsrA(DII) and suggests that RNA-induced dimer formation of Hfq(6) is a common feature of Hfq-RNA interactions.

Development of small molecules designed to modulate protein-protein interactions

Protein-protein interactions are ubiquitous, essential to almost all known biological processes, and offer attractive opportunities for therapeutic intervention. Developing small molecules that modulate protein-protein interactions is challenging, owing to the large size of protein-complex interface, the lack of well-defined binding pockets, etc. We describe a general approach based on the "privileged-structure hypothesis" [Che, Ph.D. Thesis, Washington University, 2003] - that any organic templates capable of mimicking surfaces of protein-recognition motifs are potential privileged scaffolds as protein-complex antagonists--to address the challenges inherent in the discovery of small-molecule inhibitors of protein-protein interactions.

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

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

Determination of the amino acids in yeast ribosomal protein YS11 essential for the recognition of nucleotides in 18 S ribosomal RNA

The nucleotides in domain I of 18 S rRNA that are important for the binding of the essential yeast ribosomal protein YS11 are mainly in a kink-turn motif and the terminal loop of helix 11 (H11). In the atomic structure of the Thermus thermophilus 30 S subunit, 16 amino acids in S17, the homolog of YS11, are within hydrogen bonding distance of nucleotides in 16 S rRNA. The homologous or analogous 16 amino acids in YS11 were replaced with alanine; nine of the substitutions slowed the growth of yeast cells. The most severe effects were caused by mutations R103A, N106A, K133A, T134A, and K151A. The T. thermophilus analogs of Arg103, Asn106, Thr134, and Lys151 contact nucleotides in the kink-turn motif of 16 S rRNA, whereas Lys133 contacts nucleotides in the terminal loop of H11. These contacts are predominantly with backbone phosphate and sugar oxygens in regions that deviate from A-form geometry, suggesting that YS11 recognizes the shape of its rRNA-binding site rather than reading the sequence of nucleotides. The effect of the mutations on the binding of YS11 to a domain I fragment of 18 S rRNA accorded, in general, with their effect on growth. Mutations of seven YS11 amino acids (Ser77, Met80, Arg88, Tyr97, Pro130, Ser132, and Arg136) whose homologs or analogs in S17 are within hydrogen bonding distance of nucleotides in 16 S rRNA did not affect binding. Apparently, proximities alone do not define either the amino acids or the nucleotides that are important for recognition.

ATP- and ADP-dependent modulation of RNA unwinding and strand annealing activities by the DEAD-box protein DED1

DEAD-box RNA helicases, which are involved in virtually all aspects of RNA metabolism, are generally viewed as enzymes that unwind RNA duplexes or disrupt RNA-protein interactions in an ATP-dependent manner. Here, we show in vitro that the DEAD-box protein DED1 from Saccharomyces cerevisiae promotes not only RNA unwinding but also strand annealing, the latter in such a profound fashion that the physical limit for a bimolecular association rate constant is approached. We further demonstrate that DED1 establishes an ATP-dependent steady state between unwinding and annealing, which enables the enzyme to modulate the balance between the two opposing activities through ATP and ADP concentrations. The ratio between unwinding and annealing and the degree to which both activities are ATP- and ADP-modulated are strongly influenced by structured as well as unstructured regions in the RNA substrate. Collectively, these findings expand the known functional repertoire of DEAD-box proteins and reveal the capacity of DED1 to remodel RNA in response to ADP and ATP concentrations by facilitating not only disruption but also formation of RNA duplexes.

Prediction of the functional class of lipid binding proteins from sequence-derived properties irrespective of sequence similarity

Lipid binding proteins play important roles in signaling, regulation, membrane trafficking, immune response, lipid metabolism, and transport. Because of their functional and sequence diversity, it is desirable to explore additional methods for predicting lipid binding proteins irrespective of sequence similarity. This work explores the use of support vector machines (SVMs) as such a method. SVM prediction systems are developed using 14,776 lipid binding and 133,441 nonlipid binding proteins and are evaluated by an independent set of 6,768 lipid binding and 64,761 nonlipid binding proteins. The computed prediction accuracy is 78.9, 79.5, 82.2, 79.5, 84.4, 76.6, 90.6, 79.0, and 89.9% for lipid degradation, lipid metabolism, lipid synthesis, lipid transport, lipid binding, lipopolysaccharide biosynthesis, lipoprotein, lipoyl, and all lipid binding proteins, respectively. The accuracy for the nonmember proteins of each class is 99.9, 99.2, 99.6, 99.8, 99.9, 99.8, 98.5, 99.9, and 97.0%, respectively. Comparable accuracies are obtained when homologous proteins are considered as one, or by using a different SVM kernel function. Our method predicts 86.8% of the 76 lipid binding proteins nonhomologous to any protein in the Swiss-Prot database and 89.0% of the 73 known lipid binding domains as lipid binding. These findings suggest the usefulness of SVMs for facilitating the prediction of lipid binding proteins. Our software can be accessed at the SVMProt server (http://jing.cz3.nus.edu.sg/cgi-bin/svmprot.cgi).

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

The NMR and X-ray structures of the Saccharomyces cerevisiae Vts1 SAM domain define a surface for the recognition of RNA hairpins

The SAM domain of the Saccharomyces cerevisiae post-transcriptional regulator Vts1 has a high affinity towards RNA hairpins containing a CUGGC pentaloop. We present the 1.6 Angstroms X-ray crystal structure of the Vts1 SAM domain in its unliganded state, and the NMR solution structure of this domain in its RNA-bound state. Both structures reveal a canonical five helix SAM domain flanked by additional secondary structural elements at the N and C termini. The two structures are essentially identical, implying that no major structural rearrangements occur upon RNA binding. Amide chemical shift changes map the RNA-binding site to a shallow, basic patch at the junction of helix alpha5 and the loop connecting helices alpha1 and alpha2.

Transcripts and transcript-binding proteins in mitochondria of Neurospora crassa

We analyzed expression elements of three disparate groups of mitochondrial genes in Neurospora crassa, apocytochrome b (COB), cytochrome c oxidase 1 (COX1), and the clustered ATP8-ATP6-mtATP9-COX2. To identify promoter sequences we employed the published N. crassa consensus sequence for COB and rRNA genes, and we found closely related sequences within the 5'-regions of both COX1 and the ATP8-COX2 transcriptional units. We determined that the mature COX1 RNA includes two flanking unassigned reading frame (URF) sequences, but the 3'-flanking ND1 is not included in the COX1 mRNA. The ATP8-ATP6-mtATP9-COX2 polycistronic transcript does not include an adjacent 5'-URF sequence. Primer extension analysis showed one likely 5'-end for the COX1 transcript, which is 73 nucleotides downstream of the consensus promoter sequence and is the first nucleotide 3' of the sequence for the tRNA(cys). Primer extension analysis and S1 nuclease mapping of the ATP8-COX2 RNA showed that the 5'-end for this transcript is the first nucleotide 3' of the consensus promoter sequence. We performed gel-shift experiments to detect proteins in mitochondria that bind to transcripts as possible regulatory proteins. The 5'-untranslated region (UTR) RNAs of COB, COX1, and ATP8-COX2 appear to bind both unique proteins and an overlapping group of two to four proteins of approximately 155-45 M(r). We successively deleted regions of the RNA 5'-UTRs to identify sequences that bound these proteins. Similar predicted stem-loop secondary structures were detected in the protein-binding regions of all three UTRs.

A simple motif for protein recognition in DNA secondary structures

DNA in a single-stranded form (ssDNA) exists transiently within the cell and comprises the telomeres of linear chromosomes and the genomes of some DNA viruses. As with RNA, in the single-stranded state, some DNA sequences are able to fold into complex secondary and tertiary structures that may be recognized by proteins and participate in gene regulation. To better understand how such DNA elements might fold and interact with proteins, and to compare recognition features to those of a structured RNA, we used in vitro selection to identify ssDNAs that bind an RNA-binding peptide from the HIV Rev protein with high affinity and specificity. The large majority of selected binders contain a non-Watson-Crick G.T base-pair and an adjacent C:G base-pair and both are essential for binding. This GT motif can be presented in different DNA contexts, including a nearly perfect duplex and a branched three-helix structure, and appears to be recognized in large part by arginine residues separated by one turn of an alpha-helix. Interestingly, a very similar GT motif is necessary also for protein binding and function of a well-characterized model ssDNA regulatory element from the proenkephalin promoter.

Common and specific amino acid residues in the prokaryotic polypeptide release factors RF1 and RF2: possible functional implications

Termination of protein synthesis is promoted in ribosomes by proper stop codon discrimination by class 1 polypeptide release factors (RFs). A large set of prokaryotic RFs differing in stop codon specificity, RF1 for UAG and UAA, and RF2 for UGA and UAA, was analyzed by means of a recently developed computational method allowing identification of the specificity-determining positions (SDPs) in families composed of proteins with similar but not identical function. Fifteen SDPs were identified within the RF1/2 superdomain II/IV known to be implicated in stop codon decoding. Three of these SDPs had particularly high scores. Five residues invariant for RF1 and RF2 [invariant amino acid residues (IRs)] were spatially clustered with the highest-scoring SDPs that in turn were located in two zones within the SDP/IR area. Zone 1 (domain II) included PxT and SPF motifs identified earlier by others as ‘discriminator tripeptides’. We suggest that IRs in this zone take part in the recognition of U, the first base of all stop codons. Zone 2 (domain IV) possessed two SDPs with the highest scores not identified earlier. Presumably, they also take part in stop codon binding and discrimination. Elucidation of potential functional role(s) of the newly identified SDP/IR zones requires further experiments.

Evolution from DNA to RNA recognition by the bI3 LAGLIDADG maturase

LAGLIDADG endonucleases bind across adjacent major grooves via a saddle-shaped surface and catalyze DNA cleavage. Some LAGLIDADG proteins, called maturases, facilitate splicing by group I introns, raising the issue of how a DNA-binding protein and an RNA have evolved to function together. In this report, crystallographic analysis shows that the global architecture of the bI3 maturase is unchanged from its DNA-binding homologs; in contrast, the endonuclease active site, dispensable for splicing facilitation, is efficiently compromised by a lysine residue replacing essential catalytic groups. Biochemical experiments show that the maturase binds a peripheral RNA domain 50 A from the splicing active site, exemplifying long-distance structural communication in a ribonucleoprotein complex. The bI3 maturase nucleic acid recognition saddle interacts at the RNA minor groove; thus, evolution from DNA to RNA function has been mediated by a switch from major to minor groove interaction.

The social life of ribosomal proteins

Ribosomal proteins hold a unique position in biology because their function is so closely tied to the large rRNAs of the ribosomes in all kingdoms of life. Following the determination of the complete crystal structures of both the large and small ribosomal subunits from bacteria, the functional role of the proteins has often been overlooked when focusing on rRNAs as the catalysts of translation. In this review we highlight some of the many known and important functions of ribosomal proteins, both during translation on the ribosome and in a wider context.

Is Hairpin Formation in Single-Stranded Polynucleotide Diffusion-Controlled?

An intriguing puzzle in biopolymer science is the observation that single-stranded DNA and RNA oligomers form hairpin structures on time scales of tens of microseconds, considerably slower than the estimated time for loop formation for a semiflexible polymer of similar length. To address the origin of the slow kinetics and to determine whether hairpin dynamics are diffusion-controlled, the effect of solvent viscosity (eta) on hairpin kinetics was investigated using laser temperature-jump techniques. The viscosity was varied by addition of glycerol, which significantly destabilizes hairpins. A previous study on the viscosity dependence of hairpin dynamics, in which all the changes in the measured rates were attributed to a change in solvent viscosity, reported an apparent scaling of relaxation times (tau(r)) on eta as tau(r) approximately eta(0.8). In this study, we demonstrate that if the effect of viscosity on the measured rates is not deconvoluted from the inevitable effect of change in stability, then separation of tau(r) into opening (tau(o)) and closing (tau(c)) times yields erroneous behavior, with different values (and opposite signs) of the apparent scaling exponents, tau(o) approximately eta(-0.4) and tau(c) approximately eta(1.5). Under isostability conditions, obtained by varying the temperature to compensate for the destabilizing effect of glycerol, both tau(o) and tau(c) scale as approximately eta(1.1+/-0.1). Thus, hairpin dynamics are strongly coupled to solvent viscosity, indicating that diffusion of the polynucleotide chain through the solvent is involved in the rate-determining step.

Kinetic analysis of the role of the tyrosine 13, phenylalanine 56 and glutamine 54 network in the U1A/U1 hairpin II interaction

The A protein of the U1 small nuclear ribonucleoprotein particle, interacting with its stem-loop RNA target (U1hpII), is frequently used as a paradigm for RNA binding by recognition motif domains (RRMs). U1A/U1hpII complex formation has been proposed to consist of at least two steps: electrostatically mediated alignment of both molecules followed by locking into place, based on the establishment of close-range interactions. The sequence of events between alignment and locking remains obscure. Here we examine the roles of three critical residues, Tyr13, Phe56 and Gln54, in complex formation and stability using Biacore. Our mutational and kinetic data suggest that Tyr13 plays a more important role than Phe56 in complex formation. Mutational analysis of Gln54, combined with molecular dynamics studies, points to Arg52 as another key residue in association. Based on our data and previous structural and modeling studies, we propose that electrostatic alignment of the molecules is followed by hydrogen bond formation between the RNA and Arg52, and the sequential establishment of interactions with loop bases (including Tyr13). A quadruple stack, sandwiching two bases between Phe56 and Asp92, would occur last and coincide with the rearrangement of a C-terminal helix that partially occludes the RRM surface in the free protein.

The p92 polymerase coding region contains an internal RNA element required at an early step in Tombusvirus genome replication

The replication of positive-strand RNA viral genomes involves various cis-acting RNA sequences. Generally, regulatory RNA sequences are present at or near genomic termini; however, internal replication elements (IREs) also exist. Here we report the structural and functional characterization of an IRE present in the readthrough portion of the p92 polymerase gene of Tomato bushy stunt virus. Analysis of this element in the context of a noncoding defective interfering RNA revealed a functional core structure composed of two noncontiguous segments of sequence that interact with each other to form an extended helical conformation. IRE activity required maintenance of several base-paired sections as well as two distinct structural features: (i) a short, highly conserved segment that can potentially form two different and mutually exclusive structures and (ii) an internal loop that contains a critical CC mismatch. The IRE was also shown to play an essential role within the context of the viral genome. In vivo analysis with novel RNA-based temperature-sensitive genomic mutants and translationally active subgenomic viral replicons revealed the following about the IRE: (i) it is active in the positive strand, (ii) it is dispensable late in the viral RNA replication process, and (iii) it is functionally inhibited by active translation over its sequence. Together, these results suggest that IRE activity is required in the cytosol at an early step in the viral replication process, such as template recruitment and/or replicase complex assembly.

Specific binding of tombusvirus replication protein p33 to an internal replication element in the viral RNA is essential for replication

The mechanism of template selection for genome replication in plus-strand RNA viruses is poorly understood. Using the prototypical tombusvirus, Tomato bushy stunt virus (TBSV), we show that recombinant p33 replicase protein binds specifically to an internal replication element (IRE) located within the p92 RNA-dependent RNA polymerase coding region of the viral genome. Specific binding of p33 to the IRE in vitro depends on the presence of a C.C mismatch within a conserved RNA helix. Interestingly, the absence of the p33:p33/p92 interaction domain in p33 prevented specific but allowed nonspecific RNA binding, suggesting that a multimeric form of this protein is involved in the IRE-specific interaction. Further support for the selectivity of p33 binding in vitro was provided by the inability of the replicase proteins of the closely related Turnip crinkle virus and distantly related Hepatitis C virus to specifically recognize the TBSV IRE. Importantly, there was also a strong correlation between p33:IRE complex formation in vitro and viral replication in vivo, where mutations in the IRE that disrupted selective p33 binding in vitro also abolished TBSV RNA replication both in plant and in Saccharomyces cerevisiae cells. Based on these findings and the other known properties of p33 and the IRE, it is proposed that the p33:IRE interaction provides a mechanism to selectively recruit viral RNAs into cognate viral replicase complexes. Since all genera in Tombusviridae encode comparable replicase proteins, these results may be relevant to other members of this large virus family.

The recognition of nucleotides with model beta-hairpin receptors: investigation of critical contacts and nucleotide selectivity

[reaction: see text] We have investigated the factors that contribute to binding of ATP by a designed 12-residue beta-hairpin peptide, WKWK, and have determined its selectivity for binding to the naturally occurring nucleotide triphosphates. We have previously shown that WKWK creates an ATP binding pocket on one face of the beta-hairpin consisting of two Trp and two Lys residues. Mutation of the two Lys residues on the binding face of the beta-hairpin resulted in a lower affinity, indicating that each is involved in ATP binding and that each residue contributes approximately -1.5 kcal/mol to the energy of complexation. Replacement of either Trp residue of the ATP binding pocket with Phe or Leu destabilizes the complex formed with ATP by approximately 1 kcal/mol, indicating that both Trp residues participate in interactions with ATP. For binding to the nucleotide triphosphates, the order of binding affinity was shown to follow dTTP > GTP > ATP > CTP, with differences in binding energies spanning as much as 1.6 kcal/mol. NMR analysis demonstrates that both aromatic interactions with the Trp side chains and CH-pi interactions between the ribose protons and the Trp residues may contribute significantly to binding. The results from our model system provide useful thermodynamic information regarding protein-nucleic acid interactions that occur at the surface of a beta-sheet.

An efficient synthesis of a probe for protein function: 2,3-diaminopropionic acid with orthogonal protecting groups

[reaction: see text] An efficient and cost-effective synthesis of N(alpha)-Boc2-N(beta)-Cbz-2,3-diaminopropionic acid is reported. The synthesis starts from commercially available N(alpha)-Boc-Asp(OBn)-OH and employs a Curtius rearrangement to establish the beta-nitrogen. Proper protection of the alpha-nitrogen is essential for the success of the Curtius rearrangement.

De novo design of peptides with L-alpha-nucleobase amino acids and their binding properties to the P22 boxB RNA and its mutants

A method to design novel molecules that specifically recognize a structured RNA would be a promising tool for the development of drugs or probes targeting RNA. In this study, the de novo design of the alpha-helical peptides having L-alpha-amino acids with nucleobases (nucleobase amino acids, NBAs) was carried out. Binding affinities of the peptides for a hairpin RNA derived from P22 phage were dependent on the types and positions of the NBA units they have. Some NBA peptides bound to the wild-type RNA or its mutant with high affinity and high specificity compared with the native P22 N peptide. These results indicate that the NBA units on the peptides interact with the RNA bases in a specific manner. It is demonstrated that the de novo design of peptides with the NBA units is an effective way to construct novel RNA-binding molecules.

A determination of the identity elements in yeast 18 S ribosomal RNA for the recognition of ribosomal protein YS11: the role of the kink-turn motif in helix 11

A description of the site of interaction of YS11, the yeast homolog of eubacterial S17, with 18 S rRNA was obtained by assessing the binding of the ribosomal protein, in a filter retention assay, to oligoribonucleotides that reproduce regions of 18 S rRNA. YS11 binds predominantly to domain I; the Kd value is 113 nM. The dimensions of the YS11 binding site were refined, guided by chemical protection data and by the atomic structure of the Thermus thermophilus 30 S subunit, which has the S17 recognition site in 16 S rRNA. An oligoribonucleotide that mimics helix 11, a phylogenetically conserved region in domain I, binds YS11 with a Kd value of 230 nM; a second oligoribonucleotide that contains only the kink-turn motif in helix 11 binds YS11 with a Kd value of 528 nM. Thus, helix 11 has most of the nucleotides required for the recognition of YS11. To identify those nucleotides a set of 27 transversion mutations in H11 was constructed and their contribution to the binding of YS11 determined. Mutations of nine nucleotides (U313, C314, A316, G337, C338, G347, U348, U350, and C351) increased the Kd value for YS11 binding by at least eightfold; G325U and U349A mutations increased the Kd value fivefold. Eight of the 11 mutations are in the kink-turn in H11, confirming the critical importance of the motif for YS11 recognition. The other three nucleotides are in the lower stem and the terminal loop of H11, which makes a lesser, but still important, contribution to YS11 binding. The identity elements for YS11 recognition are: A316, G325, G337, G347, U348, U349, U350, and C351. The effect of the other nucleotides that decrease binding is probably indirect, presumably they affect the conformation of the binding site but do not have contacts to YS11 amino acid residues. The eight identity element nucleotides are in regions of H11 that deviate from A-form geometry and the contacts are predominantly, if not exclusively, to backbone phosphate and sugar oxygen atoms, indicating that YS11 recognizes the shape of the rRNA binding site rather than reading the sequence of nucleotides.

The double-stranded-RNA-binding motif: interference and much more

RNA duplexes have been catapulted into the spotlight by the discovery of RNA interference and related phenomena. But double-stranded and highly structured RNAs have long been recognized as key players in cell processes ranging from RNA maturation and localization to the antiviral response in higher organisms. Penetrating insights into the metabolism and functions of such RNAs have come from the identification and study of proteins that contain the double-stranded-RNA-binding motif.

Modelling protein–RNA interactions: an electron density study of the guanidinium and formate complexes with RNA bases

The complexes formed by the double interaction established between RNA bases and guanidinium and formate ions, as a model for the interacting groups of arginine and glutamic or aspartic amino acid side chains, have been theoretically studied. A density functional theory method (B3LYP/6-31 + G**) has been used for this study. The range of interaction energies obtained allowed for a distinction between bidentate and bifurcate hydrogen bond interactions. The analysis of the electron density and the natural bond orbital analysis shows that these complexes are bound by double hydrogen bonds established between the donor and acceptor groups of guanidinium and formate respectively and those of the RNA bases. Comparisons are made with the results obtained in some previous theoretical and experimental studies.

Solution structure of a consensus stem-loop D RNA domain that plays important roles in regulating translation and replication in enteroviruses and rhinoviruses

Stem-loop D from the cloverleaf RNA is a highly conserved domain within the 5'-UTR of enteroviruses and rhinoviruses. Interaction between the stem-loop D RNA and the viral 3C or 3CD proteins constitutes an essential feature of a ribonucleoprotein complex that plays a critical role in regulating viral translation and replication. Here we report the solution NMR structure of a 38-nucleotide RNA with a sequence that encompasses the entire stem-loop D domain and corresponds to the consensus sequence found in enteroviruses and rhinoviruses. Sequence variants corresponding to Poliovirus type 1 and Coxsackievirus B3 have virtually the same structure, based on small differences in chemical shifts. A substantial number (136) of (1)H-(13)C one-bond residual dipolar coupling (RDC) values were used in the structure determination in addition to conventional distance and torsion angle restraints. Inclusion of the RDC restraints was essential for achieving well-defined structures, both globally and locally. The structure of the consensus stem-loop D is an elongated A-type helical stem capped by a UACG tetraloop with a wobble UG closing base pair. Three consecutive pyrimidine base pairs (two UU and one CU pair) are present in the middle of the helical stem, creating distinctive local structural features such as a dramatically widened major groove. A dinucleotide bulge is located near the base of the stem. The bulge itself is flexible and not as well defined as the other parts of the molecule, but the flanking base pairs are intact. The peculiar spatial arrangement of the distinctive structural elements implies that they may work synergistically to achieve optimal binding affinity and specificity toward the viral 3C or 3CD proteins.

Secondary structure of the 3' UTR of bicoid mRNA

Formation of the Bicoid morphogen gradient in early Drosophila embryos requires the pre-localization of bicoid mRNA to the anterior pole of the egg. The program of bcd mRNA localization involves multiples steps and proceeds from oogenesis until early embryogenesis. This process requires cis-elements in the 3' UTR of bcd mRNA and successive and/or concomitant critical protein interactions. Furthermore, numerous RNA elements and binding proteins contribute to regulate bcd expression. In the present paper, we investigated the secondary structure of the full length 3' UTR of the bcd mRNA, using a variety of chemical and enzymatic structural probes. This RNA probing analysis allowed us to give a detailed description of the 3' UTR of the bcd mRNA and its organization into five well-defined and independent domains (I-V). One prominent result that emerges from our data is the unexpected high degree of flexibility of the different domains relative to each others. This plasticity relies upon the open conformation of the central hinge region interconnecting domains II, III, and IV + V. Otherwise, dimerization of the 3' UTR, which participates to anchoring bcd mRNA at the anterior pole of the embryo, only results in discrete and local change in domain III. Domain I that contains sites for trans-acting factors exhibiting single stranded RNA binding specificity is mainly unstructured. By contrast, each core domains (II-V) is highly organized and folds into helices interrupted by bulges and interior loops and closed by very exposed apical loops. These elements mostly built specific determinants for trans-acting factors. Besides, these findings provide a valuable database for structure/function studies.

Shift in nucleotide conformational equilibrium contributes to increased rate of catalysis of GpAp versus GpA in barnase

The microbial ribonuclease barnase exhibits low catalytic activity toward GpN dinucleotides, where G is guanosine, p is phosphate and N represents any nucleoside. When a phosphate is added to the 3'-end, as in GpNp, substrate affinity is enhanced by one order of magnitude, and the catalytic rate by two. In order to gain insight into this phenomenon, we analyzed the nucleotide conformations and protein-nucleotide interactions of 4 ns molecular dynamics (MD) trajectories of complexes of barnase with guanylyl(3'-5') adenosine (GpA) and guanylyl(3'-5') adenosine 3'-monophosphate (GpAp), respectively, in the presence of solvent and counter ions. We found that, in a majority of the bound GpA conformations, the guanine base was firmly bound to the recognition site. The phosphate and adenosine moieties pointed into the solvent, and interactions with key catalytic residues were absent. In contrast, the bound GpAp adopted conformations in which all of the nucleotide portions remained tightly bound to the enzyme and interactions with key catalytic residues were maintained. These observations indicate that, for GpA, a significant proportion of the bound nucleotide adopts non-productive conformations and that adding the terminal phosphate as in GpAp shifts the equilibrium of the bound conformations towards structures capable of undergoing catalysis. Incorporating this property into the kinetic equations yields an increase in both the apparent rate constant (kcat) and the apparent dissociation constant (K(M)) for GpAp versus GpA. The increase in K(M), caused by the presence of additional non-productive binding modes for GpA, should however be counterbalanced by the propensity of free GpA to adopt folded conformations in solution, which are unable to bind the enzyme and by the tighter binding of GpAp (Giraldo J, Wodak SJ, Van Belle D. Conformational analysis of GpA and GpAp in aqueous solution by molecular dynamics and statistical methods. J Mol Biol 1998; 283:863-882). Addition of the terminal phosphate is shown to significantly influence the collective motion of the enzyme in a manner that fosters interactions with key catalytic residues, representing a further likely contribution to the catalytic rate enhancement.

Structural elements required for association of the Saccharomyces cerevisiae telomerase RNA with the Est2 reverse transcriptase

Telomere synthesis in most organisms depends on the action of the telomerase enzyme, which contains an RNA subunit that is stably associated with the reverse transcriptase subunit as well as additional telomerase proteins. In the budding yeast Saccharomyces cerevisiae, several structural domains that are responsible for mediating protein interactions with the telomerase RNA TLC1 have been identified. We report here the identification and characterization of a TLC1 stem-loop that is required for its interaction with the Est2 reverse transcriptase protein. This hairpin, which does not contain any bulges in the duplex stem that commonly mediate protein-RNA interaction, appears to be a part of a larger structure, as nucleotides immediately to either side of this stem-loop contribute to the interaction of TLC1 with the Est2 protein. Surprisingly, replacement of a 95-nucleotide region of the yeast telomerase RNA that is required for Est2 interaction with a 39-nucleotide pseudoknot from a distantly related telomerase RNA results in a functional telomerase enzyme. These findings suggest that the ability of the budding yeast reverse transcriptase to associate with the telomerase RNA depends on a highly structured region rather than specific sequence elements.

Electrostatic contribution of serine phosphorylation to the Drosophila SLBP--histone mRNA complex

Unlike all other metazoan mRNAs, mRNAs encoding the replication-dependent histones are not polyadenylated but end in a unique 26 nucleotide stem-loop structure. The protein that binds the 3' end of histone mRNA, the stem-loop binding protein (SLBP), is essential for histone pre-mRNA processing, mRNA translation, and mRNA degradation. Using biochemical, biophysical, and nuclear magnetic resonance (NMR) experiments, we report the first structural insight into the mechanism of SLBP-RNA recognition. In the absence of RNA, phosphorylated and unphosphorylated forms of the RNA binding and processing domain (RPD) of Drosophila SLBP (dSLBP) possess helical secondary structure but no well-defined tertiary fold. Drosophila SLBP is phosphorylated at four out of five potential serine or threonine sites in the sequence DTAKDSNSDSDSD at the extreme C-terminus, and phosphorylation at these sites is necessary for histone pre-mRNA processing. Here, we provide NMR evidence for serine phosphorylation of the C-terminus using (31)P direct-detect experiments and show that both serine phosphorylation and RNA binding are necessary for proper folding of the RPD. The electrostatic effect of protein phosphorylation can be partially mimicked by a mutant form of SLBP wherein four C-terminal serines are replaced with glutamic acids. Hence, both RNA binding and protein phosphorylation are necessary for stabilization of the SLBP RPD.