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.

Mass spectrometric amino acid sequencing of a mixture of seed storage proteins (napin) from Brassica napus, products of a multigene family.

The amino acid sequences of a number of closely related proteins ("napin") isolated from Brassica napus were determined by mass spectrometry without prior separation into individual components. Some of these proteins correspond to those previously deduced (napA, BngNAP1, and gNa), chiefly from DNA sequences. Others were found to differ to a varying extent (BngNAP1', BngNAP1A, BngNAP1B, BngNAP1C, gNa', and gNaA). The short chains of gNa and gNa' and of BngNAP1 and BngNAP1' differ by the replacement of N-terminal proline by pyroglutamic acid; the long chains of gNaA and BngNAP1B contain a six amino acid stretch, MQGQQM, which is present in gNa (according to its DNA sequence) but absent from BngNAP1 and BngNAP1C. These alternations of sequences between napin isoforms are most likely due to homologous recombination of the genetic material, but some of the changes may also be due to RNA editing. The amino acids that follow the untruncated C termini of those napin chains for which the DNA sequences are known (napA, BngNAP1, and gNa) are aromatic amino acids. This suggests that the processing of the proprotein leading to the C termini of the two chains is due to the action of a protease that specifically cleaves a G/S-F/Y/W bond.

Small regulatory RNAs in mammals

Mammalian cells harbor numerous small non-protein-coding RNAs, including small nucleolar RNAs (snoRNAs), microRNAs (miRNAs), short interfering RNAs (siRNAs) and small double-stranded RNAs, which regulate gene expression at many levels including chromatin architecture, RNA editing, RNA stability, translation, and quite possibly transcription and splicing. These RNAs are processed by multistep pathways from the introns and exons of longer primary transcripts, including protein-coding transcripts. Most show distinctive temporal- and tissue-specific expression patterns in different tissues, including embryonal stem cells and the brain, and some are imprinted. Small RNAs control a wide range of developmental and physiological pathways in animals, including hematopoietic differentiation, adipocyte differentiation and insulin secretion in mammals, and have been shown to be perturbed in cancer and other diseases. The extent of transcription of non-coding sequences and the abundance of small RNAs suggests the existence of an extensive regulatory network on the basis of RNA signaling which may underpin the development and much of the phenotypic variation in mammals and other complex organisms and which may have different genetic signatures from sequences encoding proteins.

Non-coding RNAs in the nervous system

Increasing evidence suggests that the development and function of the nervous system is heavily dependent on RNA editing and the intricate spatiotemporal expression of a wide repertoire of non-coding RNAs, including micro RNAs, small nucleolar RNAs and longer non-coding RNAs. Non-coding RNAs may provide the key to understanding the multi-tiered links between neural development, nervous system function, and neurological diseases.

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

Sendai virus RNA polymerase scanning for mRNA start sites at gene junctions

Mini-genomes expressing two reporter genes and a variable gene junction were used to study Sendai virus RNA polymerase (RdRp) scanning for the mRNA start signal of the downstream gene (gs2). We found that RdRp could scan the template efficiently as long as the initiating uridylate of gs2 (3' UCCCnnUUUC) was preceded by the conserved intergenic region (3' GAA) and the last 3 uridylates of the upstream gene end signal (ge1; 3' AUUCUUUUU). The end of the leader sequence (3' CUAAAA, which precedes gs1) could also be used for gene2 expression, but this sequence was considerably less efficient. Increasing the distance between ge1 and gs2 (up to 200 nt) led to the progressive loss of gene2 expression, in which half of gene2 expression was lost for each 70 nucleotides of intervening sequence. Beyond 200 nt, gene2 expression was lost more slowly. Our results suggest that there may be two populations of RdRp that scan at gene junctions, which can be distinguished by the efficiency with which they can scan the genome template for gs.

Photochemically induced RNA and DNA abasic sites

Two phosphoramidite building blocks were synthesized that can easily be deprotected by UV light to reveal natural abasic sites in oligoribonucleotides as well as in oligodeoxyribonucleotides. Another building block which releases a 2 ′-O-methylated abasic site upon UV radiation is also described.

Messenger RNA levels and binding sites of muscarinic acetylcholine receptors in gastrointestinal muscle layers from healthy dairy cows

Acetylcholine interacts with muscarinic receptors (M) to mediate gastrointestinal (GI) smooth muscle contractions. We have compared mRNA levels and binding sites of M(1)to M(5) in muscle tissues from fundus abomasi, pylorus, ileum, cecum, proximal loop of the ascending colon (PLAC), and external loop of the spiral colon (ELSC) of healthy dairy cows. The mRNA levels were measured by quantitative RT-PCR. The inhibition of [(3)H]-QNB (1-quinuclidinyl-[phenyl-4-(3)H]-benzilate) binding by M antagonists [atropine (M(1 - 5)), pirenzepine (M(1)), methoctramine (M(2)), 4-DAMP (M(3)), and tropicamide (M(4))] was used to identify receptors at the functional level. Maximal binding (B(max)) was determined through saturation binding with atropine as a competitor. The mRNA levels of M(1), M(2), M(3), and M(5) represented 0.2, 48, 50, and 1.8%, respectively, of the total M population, whereas mRNA of M(4) was undetectable. The mRNA levels of M(2) and of M(3) in the ileum were lower (P < 0.05) than in other GI locations, which were similar among each other. Atropine, pirenzepine, methoctramine, and 4-DAMP inhibited [(3)H]-QNB binding according to an either low- or high-affinity receptor pattern, whereas tropicamide had no effect on [(3)H]-QNB binding. The [(3)H]-QNB binding was dose-dependent and saturable. B(max) in fundus, pylorus, and PLAC was lower (P < 0.05) than in the ELSC, and in the pylorus lower (P < 0.05) than in the ileum. B(max) and mRNA levels were negatively correlated (r = -0.3; P < 0.05). In conclusion, densities of M are different among GI locations, suggesting variable importance of M for digestive functions along the GI tract.

RNA-Guided Genome Editing

Ciliates have evolved highly complex and intricately controlled pathways to ensure the precise and complete removal of all genomic sequences not required for vegetative growth. At the same time, they retain a reference copy of all their genetic information for future generations. This chapter describes how different ciliates use RNA-mediated DNA comparison processes to form new somatic nuclei from germline nuclei. While these processes vary in their precise mechanisms, they all use RNA to target genomic DNA sequences—either for retention or elimination. They also all consist of more than one individual pathway acting cooperatively—the two subsets of small RNAs in Paramecium and the guide RNAs and Piwi-interacting RNAs in Oxytricha—to ensure a strong belt-and-braces approach to consistent and precise somatic nucleus development. Nonetheless, this genome comparison approach to somatic nucleus development provides an elegant method for trans-generational environmental adaptation. Conceptually, it is easy to imagine how somatic changes that occur during vegetative growth could be transferred to meiotic offspring, while an unaltered germline genome is retained. Further research in this area will have far-reaching implications for the trans-generational adaptation of more distantly related eukaryotes, such as humans.

Organization of Highly Acetylated Chromatin around Sites of Heterogeneous Nuclear RNA Accumulation

Histones found within transcriptionally competent and active regions of the genome are highly acetylated. Moreover, these highly acetylated histones have very short half-lives. Thus, both histone acetyltransferases and histone deacetylases must enrich within or near these euchromatic regions of the interphase chromatids. Using an antibody specific for highly acetylated histone H3, we have investigated the organization of transcriptionally active and competent chromatin as well as nuclear histone acetyltransferase and deacetylase activities. We observe an exclusion of highly acetylated chromatin around the periphery of the nucleus and an enrichment near interchromatin granule clusters (IGCs). The highly acetylated chromatin is found in foci that may reflect the organization of highly acetylated chromatin into “chromonema” fibers. Transmission electron microscopy of Indian muntjac fibroblast cell nuclei indicates that the chromatin associated with the periphery of IGCs remains relatively condensed, most commonly found in domains containing chromatin folded beyond 30 nm. Using electron spectroscopic imaging, we demonstrate that IGCs are clusters of ribonucleoprotein particles. The individual granules comprise RNA-rich fibrils or globular regions that fold into individual granules. Quantitative analysis of individual granules indicates that they contain variable amounts of RNA estimated between 1.5 and >10 kb. We propose that interchromatin granules are heterogeneous nuclear RNA-containing particles, some of which may be pre-mRNA generated by nearby transcribed chromatin. An intermediary zone between the IGC and surrounding chromatin is described that contains factors with the potential to provide specificity to the localization of sequences near IGCs.

Mapping the σ70 subunit contact sites on Escherichia coli RNA polymerase with a σ70-conjugated chemical protease

The core enzyme of Escherichia coli RNA polymerase acquires essential promoter recognition and transcription initiation activities by binding one of several σ subunits. To characterize the proximity between σ70, the major σ for transcription of the growth-related genes, and the core enzyme subunits (α2ββ′), we analyzed the protein-cutting patterns produced by a set of covalently tethered FeEDTA probes [FeBABE: Fe (S)-1-(p-bromoacetamidobenzyl)EDTA]. The probes were positioned in or near conserved regions of σ70 by using seven mutants, each carrying a single cysteine residue at position 132, 376, 396, 422, 496, 517, or 581. Each FeBABE-conjugated σ70 was bound to the core enzyme, which led to cleavage of nearby sites on the β and β′ subunits (but not α). Unlike the results of random cleavage [Greiner, D. P., Hughes, K. A., Gunasekera, A. H. & Meares, C. F. (1996) Proc. Natl. Acad. Sci. USA 93, 71–75], the cut sites from different probe-modified σ70 proteins are clustered in distinct regions of the subunits. On the β subunit, cleavage is observed in two regions, one between residues 383 and 554, including the conserved C and Rif regions; and the other between 854 and 1022, including conserved region G, regions of ppGpp sensitivity, and one of the segments forming the catalytic center of RNA polymerase. On the β′ subunit, the cleavage was identified within the sequence 228–461, including β′ conserved regions C and D (which comprise part of the catalytic center).

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Different cleavage sites are aligned differently in the active site of M1 RNA, the catalytic subunit of Escherichia coli RNase P.

We have studied RNase P RNA (M1 RNA) cleavage of model tRNA precursors that are cleaved at two independent positions. Here we present data demonstrating that cleavage at both sites depends on the 2'-OH immediately 5' of the respective cleavage site. However, we show that the 2-amino group of a guanosine at the cleavage site plays a significant role in cleavage at one of these sites but not at the other. These data suggest that these two cleavage sites are handled differently by the ribozyme. This theory is supported by our finding that the cross-linking pattern between Ml RNA and tRNA precursors carrying 4-thioU showed distinct differences, depending on the location of the 4-thioU relative to the respective cleavage site. These findings lead us to suggest that different cleavage sites are aligned differently in the active site, possibly as a result of different binding modes of a substrate to M1 RNA. We discuss a model in which the interaction between the 3'-terminal "RCCA" motif (first three residues interact) of a tRNA precursor and M1 RNA plays a significant role in this process.