Function of RNA secondary structures in transcriptional attenuation of the Bacillus subtilis pyr operon

The Bacillus subtilis pyr operon is regulated by exogenous pyrimidines by a transcriptional attenuation mechanism. Transcription in vitro from pyr DNA templates specifying attenuation regions yielded terminated and read-through transcripts of the expected lengths. Addition of the PyrR regulatory protein plus UMP led to greatly increased termination. Synthetic antisense deoxyoligonucleotides were used to probe possible secondary structures in the pyr mRNA that were proposed to play roles in controlling attenuation. Oligonucleotides predicted to disrupt terminator structures suppressed termination, whereas oligonucleotides predicted to disrupt the stem of antiterminator stem-loops strongly promoted termination at the usual termination site. Oligonucleotides that disrupt a previously unrecognized stem-loop structure, called the anti-antiterminator, the formation of which interferes with formation of the downstream antiterminator, suppressed termination. We propose that transcriptional attenuation of the pyr operon is governed by switching between alternative antiterminator versus anti-antiterminator plus terminator structures, and that PyrR acts by UMP-dependent binding to and stabilization of the anti-antiterminator.

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.

Comparative analyses of the secondary structures of synthetic and intracellular yeast MFA2 mRNAs

The overall folded (global) structure of mRNA may be critical to translation and turnover control mechanisms, but it has received little experimental attention. Presented here is a comparative analysis of the basic features of the global secondary structure of a synthetic mRNA and the same intracellular eukaryotic mRNA by dimethyl sulfate (DMS) structure probing. Synthetic MFA2 mRNA of Saccharomyces cerevisiae first was examined by using both enzymes and chemical reagents to determine single-stranded and hybridized regions; RNAs with and without a poly(A) tail were compared. A folding pattern was obtained with the aid of the mfold program package that identified the model that best satisfied the probing data. A long-range structural interaction involving the 5′ and 3′ untranslated regions and causing a juxtaposition of the ends of the RNA, was examined further by a useful technique involving oligo(dT)-cellulose chromatography and antisense oligonucleotides. DMS chemical probing of A and C nucleotides of intracellular MFA2 mRNA was then done. The modification data support a very similar intracellular structure. When low reactivity of A and C residues is found in the synthetic RNA, ≈70% of the same sites are relatively more resistant to DMS modification in vivo. A slightly higher sensitivity to DMS is found in vivo for some of the A and C nucleotides predicted to be hybridized from the synthetic structural model. With this small mRNA, the translation process and mRNA-binding proteins do not block DMS modifications, and all A and C nucleotides are modified the same or more strongly than with the synthetic RNA.

Preferential interaction of the his pause RNA hairpin with RNA polymerase β subunit residues 904–950 correlates with strong transcriptional pausing

RNA secondary structures (hairpins) that form as the nascent RNA emerges from RNA polymerase are important components of many signals that regulate transcription, including some pause sites, all ρ-independent terminators, and some antiterminators. At the his leader pause site, a 5-bp-stem, 8-nt-loop pause RNA hairpin forms 11 nt from the RNA 3′ end and stabilizes a transcription complex conformation slow to react with NTP substrate. This stabilization appears to depend at least in part on an interaction with RNA polymerase. We tested for RNA hairpin interaction with the paused polymerase by crosslinking 5-iodoUMP positioned specifically in the hairpin loop. In the paused conformation, strong and unusual crosslinking of the pause hairpin to β904–950 replaced crosslinking to β′ and to other parts of β that occurred in nonpaused complexes prior to hairpin formation. These changes in nascent RNA interactions may inhibit reactive alignment of the RNA 3′ end in the paused complex and be related to events at ρ-independent terminators.

Targeting nucleic acid secondary structures by antisense oligonucleotides designed through in vitro selection.

Using an in vitro selection approach, we have isolated oligonucleotides that can bind to a DNA hairpin structure. Complex formation of these oligonucleotides with the target hairpin involves some type of triple-stranded structure with noncanonical interaction, as indicated by bandshift assays and footprinting studies. The selected oligomers can block restriction endonuclease cleavage of the target hairpin in a sequence-specific manner. We demonstrate that in vitro selection can extend the antisense approach to functional targeting of secondary structure motifs. This could provide a basis for interfering with regulatory processes mediated by a variety of nucleic acid structures.

GOBASE: the organelle genome database

GOBASE (http://megasun.bch.umontreal.ca/gobase/) is a network-accessible biological database, which is unique in bringing together diverse biological data on organelles with taxonomically broad coverage, and in furnishing data that have been exhaustively verified and completed by experts. So far, we have focused on mitochondrial data: GOBASE contains all published nucleotide and protein sequences encoded by mitochondrial genomes, selected RNA secondary structures of mitochondria-encoded molecules, genetic maps of completely sequenced genomes, taxonomic information for all species whose sequences are present in the database and organismal descriptions of key protistan eukaryotes. All of these data have been integrated and organized in a formal database structure to allow sophisticated biological queries using terms that are inherent in biological concepts. Most importantly, data have been validated, completed, corrected and standardized, a prerequisite of meaningful analysis. In addition, where critical data are lacking, such as genetic maps and RNA secondary structures, they are generated by the GOBASE team and collaborators, and added to the database. The database is implemented in a relational database management system, but features an object-oriented view of the biological data through a Web/Genera-generated World Wide Web interface. Finally, we have developed software for database curation (i.e. data updates, validation and correction), which will be described in some detail in this paper.

Conserved stem–loop structures in the HIV-1 RNA region containing the A3 3′ splice site and its cis-regulatory element: possible involvement in RNA splicing

The HIV-1 transcript is alternatively spliced to over 30 different mRNAs. Whether RNA secondary structure can influence HIV-1 RNA alternative splicing has not previously been examined. Here we have determined the secondary structure of the HIV-1/BRU RNA segment, containing the alternative A3, A4a, A4b, A4c and A5 3′ splice sites. Site A3, required for tat mRNA production, is contained in the terminal loop of a stem–loop structure (SLS2), which is highly conserved in HIV-1 and related SIVcpz strains. The exon splicing silencer (ESS2) acting on site A3 is located in a long irregular stem–loop structure (SLS3). Two SLS3 domains were protected by nuclear components under splicing condition assays. One contains the A4c branch points and a putative SR protein binding site. The other one is adjacent to ESS2. Unexpectedly, only the 3′ A residue of ESS2 was protected. The suboptimal A3 polypyrimidine tract (PPT) is base paired. Using site-directed mutagenesis and transfection of a mini-HIV-1 cDNA into HeLa cells, we found that, in a wild-type PPT context, a mutation of the A3 downstream sequence that reinforced SLS2 stability decreased site A3 utilization. This was not the case with an optimized PPT. Hence, sequence and secondary structure of the PPT may cooperate in limiting site A3 utilization.

RNA editing in Trypanosoma brucei: characterization of gRNA U-tail interactions with partially edited mRNA substrates

Guide RNAs (gRNAs), key components of the RNA editing reaction in Trypanosoma brucei, direct the insertion and deletion of uridylate (U) residues. Analyses of gRNAs reveal three functional elements. The 5′-end of the gRNA contains the anchor, which is responsible for selection and binding to the pre-edited mRNA. The second element (the guiding region) provides the information required for editing. At the 3′-end of the gRNA is a non-encoded U-tail, whose function remains unclear. However, the cleavage–ligation model for editing proposes that the U-tail binds to purine-rich regions upstream of editing sites, thereby strengthening the interaction and holding onto the 5′ cleavage product. Our previous studies demonstrated that the U-tail interacts with upstream sequences and may play roles in both stabilization and tethering. These studies also indicated that the U-tail interactions involved mRNA regions that were to be subsequently edited. This raised the question of what happens to the mRNA–U-tail interaction as editing proceeds in the 3′→5′ direction. We examined gCYb-558 and its U-tail interaction with 5′CYbUT and two partially edited 5′CYb substrates. Our results indicate that the 3′-end of the U-tail interacts with the same sequence in all three mRNAs. Predicted secondary structures using crosslinking data suggest that a similar structure is maintained as editing proceeds. These results indicate that the role of the U-tail may also involve maintenance of important secondary structure motifs.

Ribosomal protein S7 from Escherichia coli uses the same determinants to bind 16S ribosomal RNA and its messenger RNA

Ribosomal protein S7 from Escherichia coli binds to the lower half of the 3′ major domain of 16S rRNA and initiates its folding. It also binds to its own mRNA, the str mRNA, and represses its translation. Using filter binding assays, we show in this study that the same mutations that interfere with S7 binding to 16S rRNA also weaken its affinity for its mRNA. This suggests that the same protein regions are responsible for mRNA and rRNA binding affinities, and that S7 recognizes identical sequence elements within the two RNA targets, although they have dissimilar secondary structures. Overexpression of S7 is known to inhibit bacterial growth. This phenotypic growth defect was relieved in cells overexpressing S7 mutants that bind poorly the str mRNA, confirming that growth impairment is controlled by the binding of S7 to its mRNA. Interestingly, a mutant with a short deletion at the C-terminus of S7 was more detrimental to cell growth than wild-type S7. This suggests that the C-terminal portion of S7 plays an important role in ribosome function, which is perturbed by the deletion.

RNA–protein hybrid ribozymes that efficiently cleave any mRNA independently of the structure of the target RNA

Ribozyme activity in vivo depends on achieving high-level expression, intracellular stability, target colocalization, and cleavage site access. At present, target site selection is problematic because of unforeseeable secondary and tertiary RNA structures that prevent cleavage. To overcome this design obstacle, we wished to engineer a ribozyme that could access any chosen site. To create this ribozyme, the constitutive transport element (CTE), an RNA motif that has the ability to interact with intracellular RNA helicases, was attached to our ribozymes so that the helicase-bound, hybrid ribozymes would be produced in cells. This modification significantly enhanced ribozyme activity in vivo, permitting cleavage of sites previously found to be inaccessible. To confer cleavage enhancement, the CTE must retain helicase-binding activity. Binding experiments demonstrated the likely involvement of RNA helicase(s). We found that attachment of the RNA motif to our tRNA ribozymes leads to cleavage in vivo at the chosen target site regardless of the local RNA secondary or tertiary structure.

Cold shock induces a major ribosomal-associated protein that unwinds double-stranded RNA in Escherichia coli.

A 70-kDa protein was specifically induced in Escherichia coli when the culture temperature was shifted from 37 to 15 degrees C. The protein was identified to be the product of the deaD gene (reassigned csdA) encoding a DEAD-box protein. Furthermore, after the shift from 37 to 15 degrees C, CsdA was exclusively localized in the ribosomal fraction and became a major ribosomal-associated protein in cells grown at 15 degrees C. The csdA deletion significantly impaired cell growth and the synthesis of a number of proteins, specifically the derepression of heat-shock proteins, at low temperature. Purified CsdA was found to unwind double-stranded RNA in the absence of ATP. Therefore, the requirement for CsdA in derepression of heat-shock protein synthesis is a cold shock-induced function possibly mediated by destabilization of secondary structures previously identified in the rpoH mRNA.

Analysis of the relation between the sequence and secondary and three-dimensional structures of immunoglobulin molecules.

Methods of structural and statistical analysis of the relation between the sequence and secondary and three-dimensional structures are developed. About 5000 secondary structures of immunoglobulin molecules from the Kabat data base were predicted. Two statistical analyses of amino acids reveal 47 universal positions in strands and loops. Eight universally conservative positions out of the 47 are singled out because they contain the same amino acid in > 90% of all chains. The remaining 39 positions, which we term universally alternative positions, were divided into five groups: hydrophobic, charged and polar, aromatic, hydrophilic, and Gly-Ala, corresponding to the residues that occupied them in almost all chains. The analysis of residue-residue contacts shows that the 47 universal positions can be distinguished by the number and types of contacts. The calculations of contact maps in the 29 antibody structures revealed that residues in 24 of these 47 positions have contacts only with residues of antiparallel beta-strands in the same beta-sheet and residues in the remaining 23 positions always have far-away contacts with residues from other beta-sheets as well. In addition, residues in 6 of the 47 universal positions are also involved in interactions with residues of the other variable or constant domains.

Crystal structure of the ribosomal RNA domain essential for binding elongation factors

The structure of a 29-nucleotide RNA containing the sarcin/ricin loop (SRL) of rat 28 S rRNA has been determined at 2.1 Å resolution. Recognition of the SRL by elongation factors and by the ribotoxins, sarcin and ricin, requires a nearly universal dodecamer sequence that folds into a G-bulged cross-strand A stack and a GAGA tetraloop. The juxtaposition of these two motifs forms a distorted hairpin structure that allows direct recognition of bases in both grooves as well as recognition of nonhelical backbone geometry and two 5′-unstacked purines. Comparisons with other RNA crystal structures establish the cross-strand A stack and the GNRA tetraloop as defined and modular RNA structural elements. The conserved region at the top is connected to the base of the domain by a region presumed to be flexible because of the sparsity of stabilizing contacts. Although the conformation of the SRL RNA previously determined by NMR spectroscopy is similar to the structure determined by x-ray crystallography, significant differences are observed in the “flexible” region and to a lesser extent in the G-bulged cross-strand A stack.

The extended environment of mononuclear metal centers in protein structures

The objectives of this and the following paper are to identify commonalities and disparities of the extended environment of mononuclear metal sites centering on Cu, Fe, Mn, and Zn. The extended environment of a metal site within a protein embodies at least three layers: the metal core, the ligand group, and the second shell, which is defined here to consist of all residues distant less than 3.5 Å from some ligand of the metal core. The ligands and second-shell residues can be characterized in terms of polarity, hydrophobicity, secondary structures, solvent accessibility, hydrogen-bonding interactions, and membership in statistically significant residue clusters of different kinds. Findings include the following: (i) Both histidine ligands of type I copper ions exclusively attach the Nδ1 nitrogen of the histidine imidazole ring to the metal, whereas histidine ligands for all mononuclear iron ions and nearly all type II copper ions are ligated via the Nɛ2 nitrogen. By contrast, multinuclear copper centers are coordinated predominantly by histidine Nɛ2, whereas diiron histidine contacts are predominantly Nδ1. Explanations in terms of steric differences between Nδ1 and Nɛ2 are considered. (ii) Except for blue copper (type I), the second-shell composition favors polar residues. (iii) For blue copper, the second shell generally contains multiple methionine residues, which are elements of a statistically significant histidine–cysteine–methionine cluster. Almost half of the second shell of blue copper consists of solvent-accessible residues, putatively facilitating electron transfer. (iv) Mononuclear copper atoms are never found with acidic carboxylate ligands, whereas single Mn2+ ion ligands are predominantly acidic and the second shell tends to be mostly buried. (v) The extended environment of mononuclear Fe sites often is associated with histidine–tyrosine or histidine–acidic clusters.

SRPDB (Signal Recognition Particle Database)

Signal recognition particle (SRP) is a stable cytoplasmic ribonucleoprotein complex that serves to translocate secretory proteins across membranes during translation. The SRP Database (SRPDB) provides compilations of SRP components, ordered alphabetically and phylogenetically. Alignments emphasize phylogenetically-supported base pairs in SRP RNA and conserved residues in the proteins. Data are provided in various formats including a column arrangement for improved access and simplified computational usability. Included are motifs for identification of new sequences, SRP RNA secondary structure diagrams, 3-D models and links to high-resolution structures. This release includes 11 new SRP RNA sequences (total of 129), two protein SRP9 sequences (total of seven), two protein SRP14 sequences (total of 10), two protein SRP19 sequences (total of 16), 10 new SRP54 (ffh) sequences (total of 66), two protein SRP68 sequences (total of seven) and two protein SRP72 sequences (total of nine). Seven sequences of the SRP receptor α-subunit and its FtsY homolog (total of 51) are new. Also considered are β-subunit of SRP receptor, Flhf, Hbsu, CaM kinase II and cpSRP43. Access to SRPDB is at http://psyche.uthct.edu/dbs/SRPDB/SRPDB.html and the European mirror http://www.medkem.gu.se/dbs/SRPDB/SRPDB.html