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.

Determining RNA solution structure by segmental isotopic labeling and NMR: application to Caenorhabditis elegans spliced leader RNA 1.

Recent developments in multidimensional heteronuclear NMR spectroscopy and large-scale synthesis of uniformly 13C- and 15N-labeled oligonucleotides have greatly improved the prospects for determination of the solution structure of RNA. However, there are circumstances in which it may be advantageous to label only a segment of the entire RNA chain. For example, in a larger RNA molecule the structural question of interest may reside in a localized domain. Labeling only the corresponding nucleotides simplifies the spectrum and resonance assignments because one can filter proton spectra for coupling to 13C and 15N. Another example is in resolving alternative secondary structure models that are indistinguishable in imino proton connectivities. Here we report a general method for enzymatic synthesis of quantities of segmentally labeled RNA molecules required for NMR spectroscopy. We use the method to distinguish definitively two competing secondary structure models for the 5' half of Caenorhabditis elegans spliced leader RNA by comparison of the two-dimensional [15N] 1H heteronuclear multiple quantum correlation spectrum of the uniformly labeled sample with that of a segmentally labeled sample. The method requires relatively small samples; solutions in the 200-300 microM concentration range, with a total of 30 nmol or approximately 40 micrograms of RNA in approximately 150 microliters, give strong NMR signals in a short accumulation time. The method can be adapted to label an internal segment of a larger RNA chain for study of localized structural problems. This definitive approach provides an alternative to the more common enzymatic and chemical footprinting methods for determination of RNA secondary structure.

Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure

Single chicken erythrocyte chromatin fibers were stretched and released at room temperature with force-measuring laser tweezers. In low ionic strength, the stretch-release curves reveal a process of continuous deformation with little or no internucleosomal attraction. A persistence length of 30 nm and a stretch modulus of ≈5 pN is determined for the fibers. At forces of 20 pN and higher, the fibers are modified irreversibly, probably through the mechanical removal of the histone cores from native chromatin. In 40–150 mM NaCl, a distinctive condensation-decondensation transition appears between 5 and 6 pN, corresponding to an internucleosomal attraction energy of ≈2.0 kcal/mol per nucleosome. Thus, in physiological ionic strength the fibers possess a dynamic structure in which the fiber locally interconverting between “open” and “closed” states because of thermal fluctuations.

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.

Smoothness within ruggedness: the role of neutrality in adaptation.

RNA secondary structure folding algorithms predict the existence of connected networks of RNA sequences with identical structure. On such networks, evolving populations split into subpopulations, which diffuse independently in sequence space. This demands a distinction between two mutation thresholds: one at which genotypic information is lost and one at which phenotypic information is lost. In between, diffusion enables the search of vast areas in genotype space while still preserving the dominant phenotype. By this dynamic the success of phenotypic adaptation becomes much less sensitive to the initial conditions in genotype space.

Intron self-complementarity enforces exon inclusion in a yeast pre-mRNA

Skipping of internal exons during removal of introns from pre-mRNA must be avoided for proper expression of most eukaryotic genes. Despite significant understanding of the mechanics of intron removal, mechanisms that ensure inclusion of internal exons in multi-intron pre-mRNAs remain mysterious. Using a natural two-intron yeast gene, we have identified distinct RNA–RNA complementarities within each intron that prevent exon skipping and ensure inclusion of internal exons. We show that these complementarities are positioned to act as intron identity elements, bringing together only the appropriate 5′ splice sites and branchpoints. Destroying either intron self-complementarity allows exon skipping to occur, and restoring the complementarity using compensatory mutations rescues exon inclusion, indicating that the elements act through formation of RNA secondary structure. Introducing new pairing potential between regions near the 5′ splice site of intron 1 and the branchpoint of intron 2 dramatically enhances exon skipping. Similar elements identified in single intron yeast genes contribute to splicing efficiency. Our results illustrate how intron secondary structure serves to coordinate splice site pairing and enforce exon inclusion. We suggest that similar elements in vertebrate genes could assist in the splicing of very large introns and in the evolution of alternative splicing.

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

Homologous genetic recombination as an intrinsic dynamic property of a DNA structure induced by RecA/Rad51-family proteins: A possible advantage of DNA over RNA as genomic material

Heteroduplex joints are general intermediates of homologous genetic recombination in DNA genomes. A heteroduplex joint is formed between a single-stranded region (or tail), derived from a cleaved parental double-stranded DNA, and homologous regions in another parental double-stranded DNA, in a reaction mediated by the RecA/Rad51-family of proteins. In this reaction, a RecA/Rad51-family protein first forms a filamentous complex with the single-stranded DNA, and then interacts with the double-stranded DNA in a search for homology. Studies of the three-dimensional structures of single-stranded DNA bound either to Escherichia coli RecA or Saccharomyces cerevisiae Rad51 have revealed a novel extended DNA structure. This structure contains a hydrophobic interaction between the 2′ methylene moiety of each deoxyribose and the aromatic ring of the following base, which allows bases to rotate horizontally through the interconversion of sugar puckers. This base rotation explains the mechanism of the homology search and base-pair switch between double-stranded and single-stranded DNA during the formation of heteroduplex joints. The pivotal role of the 2′ methylene-base interaction in the heteroduplex joint formation is supported by comparing the recombination of RNA genomes with that of DNA genomes. Some simple organisms with DNA genomes induce homologous recombination when they encounter conditions that are unfavorable for their survival. The extended DNA structure confers a dynamic property on the otherwise chemically and genetically stable double-stranded DNA, enabling gene segment rearrangements without disturbing the coding frame (i.e., protein-segment shuffling). These properties may give an extensive evolutionary advantage to DNA.

The solution structure of the Zα domain of the human RNA editing enzyme ADAR1 reveals a prepositioned binding surface for Z-DNA

Double-stranded RNA deaminase I (ADAR1) contains the Z-DNA binding domain Zα. Here we report the solution structure of free Zα and map the interaction surface with Z-DNA, confirming roles previously assigned to residues by mutagenesis. Comparison with the crystal structure of the (Zα)2/Z-DNA complex shows that most Z-DNA contacting residues in free Zα are prepositioned to bind Z-DNA, thus minimizing the entropic cost of binding. Comparison with homologous (α+β)helix–turn–helix/B-DNA complexes suggests that binding of Zα to B-DNA is disfavored by steric hindrance, but does not eliminate the possibility that related domains may bind to both B- and Z-DNA.

Altered 3′-terminal RNA structure in phage Qβ adapted to host factor-less Escherichia coli

The RNA phage Qβ requires for the replication of its genome an RNA binding protein called Qβ host factor or Hfq protein. Our previous results suggested that this protein mediates the access of replicase to the 3′-end of the Qβ plus strand RNA. Here we report the results of an evolutionary experiment in which phage Qβ was adapted to an Escherichia coli Q13 host strain with an inactivated host factor (hfq) gene. This strain initially produced phage at a titer ≈10,000-fold lower than the wild-type strain and with minute plaque morphology, but after 12 growth cycles, phage titer and plaque size had evolved to levels near those of the wild-type host. RNAs isolated from adapted Qβ mutants were efficient templates for replicase without host factor in vitro. Electron microscopy showed that mutant RNAs, in contrast to wild-type RNA, efficiently interacted with replicase at the 3′-end in the absence of host factor. The same set of four mutations in the 3′-terminal third of the genome was found in several independently evolved phage clones. One mutation disrupts the base pairing of the 3′-terminal CCCoh sequence, suggesting that the host factor stimulates activity of the wild-type RNA template by melting out its 3′-end.

Crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus

We report the crystal structure of the RNA-dependent RNA polymerase of hepatitis C virus, a major human pathogen, to 2.8-Å resolution. This enzyme is a key target for developing specific antiviral therapy. The structure of the catalytic domain contains 531 residues folded in the characteristic fingers, palm, and thumb subdomains. The fingers subdomain contains a region, the “fingertips,” that shares the same fold with reverse transcriptases. Superposition to the available structures of the latter shows that residues from the palm and fingertips are structurally equivalent. In addition, it shows that the hepatitis C virus polymerase was crystallized in a closed fingers conformation, similar to HIV-1 reverse transcriptase in ternary complex with DNA and dTTP [Huang H., Chopra, R., Verdine, G. L. & Harrison, S. C. (1998) Science 282, 1669–1675]. This superposition reveals the majority of the amino acid residues of the hepatitis C virus enzyme that are likely to be implicated in binding to the replicating RNA molecule and to the incoming NTP. It also suggests a rearrangement of the thumb domain as well as a possible concerted movement of thumb and fingertips during translocation of the RNA template-primer in successive polymerization rounds.

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 Role of the Schizosaccharomyces pombe gar2 Protein in Nucleolar Structure and Function Depends on the Concerted Action of its Highly Charged N Terminus and its RNA-binding Domains

Nonribosomal nucleolar protein gar2 is required for 18S rRNA and 40S ribosomal subunit production in Schizosaccharomyces pombe. We have investigated the consequences of the absence of each structural domain of gar2 on cell growth, 18S rRNA production, and nucleolar structure. Deletion of gar2 RNA-binding domains (RBDs) causes stronger inhibition of growth and 18S rRNA accumulation than the absence of the whole protein, suggesting that other factors may be titrated by its remaining N-terminal basic/acidic serine-rich domain. These drastic functional defects correlate with striking nucleolar hypertrophy. Point mutations in the conserved RNP1 motifs of gar2 RBDs supposed to inhibit RNA–protein interactions are sufficient to induce severe nucleolar modifications but only in the presence of the N-terminal domain of the protein. Gar2 and its mutants also distribute differently in glycerol gradients: gar2 lacking its RBDs is found either free or assembled into significantly larger complexes than the wild-type protein. We propose that gar2 helps the assembly on rRNA of factors necessary for 40S subunit synthesis by providing a physical link between them. These factors may be recruited by the N-terminal domain of gar2 and may not be released if interaction of gar2 with rRNA is impaired.

Crystal structure of yeast initiation factor 4A, a DEAD-box RNA helicase

The eukaryotic translation initiation factor 4A (eIF4A) is a member of the DEA(D/H)-box RNA helicase family, a diverse group of proteins that couples an ATPase activity to RNA binding and unwinding. Previous work has provided the structure of the amino-terminal, ATP-binding domain of eIF4A. Extending those results, we have solved the structure of the carboxyl-terminal domain of eIF4A with data to 1.75 Å resolution; it has a parallel α-β topology that superimposes, with minor variations, on the structures and conserved motifs of the equivalent domain in other, distantly related helicases. Using data to 2.8 Å resolution and molecular replacement with the refined model of the carboxyl-terminal domain, we have completed the structure of full-length eIF4A; it is a “dumbbell” structure consisting of two compact domains connected by an extended linker. By using the structures of other helicases as a template, compact structures can be modeled for eIF4A that suggest (i) helicase motif IV binds RNA; (ii) Arg-298, which is conserved in the DEA(D/H)-box RNA helicase family but is absent from many other helicases, also binds RNA; and (iii) motifs V and VI “link” the carboxyl-terminal domain to the amino-terminal domain through interactions with ATP and the DEA(D/H) motif, providing a mechanism for coupling ATP binding and hydrolysis with conformational changes that modulate RNA binding.

Solution structure of the A loop of 23S ribosomal RNA

The A loop is an essential RNA component of the ribosome peptidyltransferase center that directly interacts with aminoacyl (A)-site tRNA. The A loop is highly conserved and contains a ubiquitous 2′-O-methyl ribose modification at position U2552. Here, we present the solution structure of a modified and unmodified A-loop RNA to define both the A-loop fold and the structural impact of the U2552 modification. Solution data reveal that the A-loop RNA has a compact structure that includes a noncanonical base pair between C2556 and U2552. NMR evidence is presented that the N3 position of C2556 has a shifted pKa and that protonation at C2556-N3 changes the C-U pair geometry. Our data indicate that U2552 methylation modifies the A-loop fold, in particular the dynamics and position of residues C2556 and U2555. We compare our structural data with the structure of the A loop observed in a recent 50S crystal structure [Ban, N., Nissen, P., Hansen, J., Moore, P. B. & Steitz, T. A. (2000) Science 289, 905–920; Nissen, P., Hansen, J., Ban, N., Moore, P. B. & Steitz, T. A. (2000) Science 289, 920–930]. The solution and crystal structures of the A loop are dramatically different, suggesting that a structural rearrangement of the A loop must occur on docking into the peptidyltransferase center. Possible roles of this docking event, the shifted pKa of C2556 and the U2552 2′-O-methylation in the mechanism of translation, are discussed.