Understanding the errors of SHAPE-directed RNA structure modeling

Single-nucleotide-resolution chemical mapping for structured RNA is being rapidly advanced by new chemistries, faster readouts, and coupling to computational algorithms. Recent tests have shown that selective 2'-hydroxyl acylation by primer extension (SHAPE) can give near-zero error rates (0-2%) in modeling the helices of RNA secondary structure. Here, we benchmark the method using six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from Escherichia coli, the P4-P6 domain of the Tetrahymena group I ribozyme, and ligand-bound domains from riboswitches for adenine, cyclic di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs gave an overall false negative rate (FNR) of 17% and a false discovery rate (FDR) of 21%, with at least one helix prediction error in five of the six cases. Extensive variations of data processing, normalization, and modeling parameters did not significantly mitigate modeling errors. Only one varation, filtering out data collected with deoxyinosine triphosphate during primer extension, gave a modest improvement (FNR = 12%, and FDR = 14%). The residual structure modeling errors are explained by the insufficient information content of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping analysis. Beyond these benchmark cases, bootstrapping suggests a low level of confidence (<50%) in the majority of helices in a previously proposed SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA modeling is not always unambiguous, and helix-by-helix confidence estimates, as described herein, may be critical for interpreting results from this powerful methodology.Comment: Biochemistry, Article ASAP (Aug. 15, 2011

Phylogenetic and Molecular Characterization of a 23S Ribosomal-Rna Gene Positions the Genus Campylobacter in the Epsilon-Subdivision of the Proteobacteria and Shows That the Presence of Transcribed Spacers Is Common in Campylobacter Spp

The nucleotide sequence of a 23S rRNA gene of Campylobacter coli VC167 was determined. The primary sequence of the C. coli 23S rRNA was deduced, and a secondary-structure model was constructed. Comparison with Escherichia coli 23S rRNA showed a major difference in the C. coli rRNA at approximately position 1170 (E. coli numbering) in the form of an extra sequence block approximately 147 bp long. PCR analysis of 31 other strains of C. coli and C. jejuni showed that 69% carried a transcribed spacer of either ca, 147 or ca. 37 bp. Comparison of all sequenced Campylobacter transcribed spacers showed that the Campylobacter inserts were related in sequence and percent G+C content. All Campylobacter strains carrying transcribed spacers in their 23S rRNA genes produced fragmented 23S rRNAs. Other strains which produced unfragmented 23S rRNAs did not appear to carry transcribed spacers at this position in their 23S rRNA genes. At the 1850 region (E. coli numbering), Campylobacter 23S rRNA displayed a base pairing signature most like that of the beta and gamma subdivisions of the class Proteobacteria, but in the 270 region, Campylobacter 23S rRNA displayed a helix signature which distinguished it from the alpha, beta, and gamma subdivisions. Phylogenetic analysis comparing C. coli VC167 23S rRNA and a C. jejuni TGH9011 (ATCC 43431) 23S rRNA with 53 other completely sequenced (eu)bacterial 23S rRNAs showed that the two campylobacters form a sister group to the alpha, beta, and gamma proteobacterial 23S rRNAs, a positioning consistent with the idea that the genus Campylobacter belongs to the epsilon subdivision of the class Proteobacteria

Mortality in women given diethylstilbestrol during pregnancy

We used Cox regression analyses to assess mortality outcomes in a combined cohort of 7675 women who received diethylstilbestrol (DES) through clinical trial participation or prenatal care. In the combined cohort, the RR for DES in relation to all-cause mortality was 1.06 (95% CI=0.98–1.16), and 1.11 (95% CI=1.02–1.21) after adjusting for covariates and omitting breast cancer deaths. The RR was 1.07 (95% CI=0.94–1.23) for overall cancer mortality, and remained similar after adjusting for covariates and omitting breast cancer deaths. The RR was 1.27 (95% CI=0.96–1.69) for DES and breast cancer, and 1.38 (95% CI=1.03–1.85) after covariate adjustment. The RR was 1.82 in trial participants and 1.12 in the prenatal care cohort, but the DES–cohort interaction was not significant (P=0.15). Diethylstilbestrol did not increase mortality from gynaecologic cancers. In summary, diethylstilbestrol was associated with a slight but significant increase in all-cause mortality, but was not significantly associated with overall cancer or gynaecological cancer mortality. The association with breast cancer mortality was more evident in trial participants, who received high DES doses

A Modern Mode of Activation for Nucleic Acid Enzymes

Through evolution, enzymes have developed subtle modes of activation in order to ensure the sufficiently high substrate specificity required by modern cellular metabolism. One of these modes is the use of a target-dependent module (i.e. a docking domain) such as those found in signalling kinases. Upon the binding of the target to a docking domain, the substrate is positioned within the catalytic site. The prodomain acts as a target-dependent module switching the kinase from an off state to an on state. As compared to the allosteric mode of activation, there is no need for the presence of a third partner. None of the ribozymes discovered to date have such a mode of activation, nor does any other known RNA. Starting from a specific on/off adaptor for the hepatitis delta virus ribozyme, that differs but has a mechanism reminiscent of this signalling kinase, we have adapted this mode of activation, using the techniques of molecular engineering, to both catalytic RNAs and DNAs exhibiting various activities. Specifically, we adapted three cleaving ribozymes (hepatitis delta virus, hammerhead and hairpin ribozymes), a cleaving 10-23 deoxyribozyme, a ligating hairpin ribozyme and an artificially selected capping ribozyme. In each case, there was a significant gain in terms of substrate specificity. Even if this mode of control is unreported for natural catalytic nucleic acids, its use needs not be limited to proteinous enzymes. We suggest that the complexity of the modern cellular metabolism might have been an important selective pressure in this evolutionary process

Attachment style, assertive communication, and safer-sex behavior

This research tested the proposition that the effect of attachment security on safer-sex practice may be mediated by communication patterns. One hundred eighty-five undergraduate students completed questionnaire measures of attachment, assertiveness, and attitudes to communication about AIDS. Eight weeks later, they reported on their practice of safer sex in the period since the first testing session. Hierarchical regressions showed that at Step 1, anxiety about relationships (a measure of insecure attachment) was associated with less safer-sex practice, for all outcome measures. Attitudes to communication about AIDS added to the prediction of general reports of safer-sex practice: in line with the mediational model, anxiety about relationships became unimportant as a predictor when communication variables were included. Communication variables failed to add to the prediction of safer sex on the most recent encounter, and both anxiety about relationships and attitudes to communication about AIDS predicted condom use. Some gender differences in patterns of prediction were noted. The results are discussed in terms of attachment style and its links with the negotiation of sexual practice and relationship issues

Loss of Ribosomal Protein L11 Affects Zebrafish Embryonic Development through a p53-Dependent Apoptotic Response

Ribosome is responsible for protein synthesis in all organisms and ribosomal proteins (RPs) play important roles in the formation of a functional ribosome. L11 was recently shown to regulate p53 activity through a direct binding with MDM2 and abrogating the MDM2-induced p53 degradation in response to ribosomal stress. However, the studies were performed in cell lines and the significance of this tumor suppressor function of L11 has yet to be explored in animal models. To investigate the effects of the deletion of L11 and its physiological relevance to p53 activity, we knocked down the rpl11 gene in zebrafish and analyzed the p53 response. Contrary to the cell line-based results, our data indicate that an L11 deficiency in a model organism activates the p53 pathway. The L11-deficient embryos (morphants) displayed developmental abnormalities primarily in the brain, leading to embryonic lethality within 6–7 days post fertilization. Extensive apoptosis was observed in the head region of the morphants, thus correlating the morphological defects with apparent cell death. A decrease in total abundance of genes involved in neural patterning of the brain was observed in the morphants, suggesting a reduction in neural progenitor cells. Upregulation of the genes involved in the p53 pathway were observed in the morphants. Simultaneous knockdown of the p53 gene rescued the developmental defects and apoptosis in the morphants. These results suggest that ribosomal dysfunction due to the loss of L11 activates a p53-dependent checkpoint response to prevent improper embryonic development

Structural Constraints Identified with Covariation Analysis in Ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab’s new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab’s Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair

Comparisons between Chemical Mapping and Binding to Isoenergetic Oligonucleotide Microarrays Reveal Unexpected Patterns of Binding to the Bacillus subtilis RNase P RNA Specificity Domain†

ABSTRACT: Microarrays with isoenergetic pentamer and hexamer 20-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on theRNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1MNaþ and 10 mMMg2þ. Binding of RNA to RNA is important for many natural func-tions, includingproteinsynthesis (1,2), translationregulation (3,4), gene silencing (5, 6), metabolic regulation (7), RNAmodification (8, 9), etc. (10-13). Binding of oligonucleotides toRNAs is impor-tant for therapeutic approaches, such as siRNA, ribozymes, and antisense therapy (14, 15).Much remains to bediscovered, however, of the rules for predicting binding sites andpotential therapeutics

Binding of Ribosomal Proteins to RNA Covalently Coupled to Agarose

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/65486/1/j.1432-1033.1977.tb11554.x.pd

Arginine Cofactors on the Polymerase Ribozyme

The RNA world hypothesis states that the early evolution of life went through a stage in which RNA served both as genome and as catalyst. The central catalyst in an RNA world organism would have been a ribozyme that catalyzed RNA polymerization to facilitate self-replication. An RNA polymerase ribozyme was developed previously in the lab but it is not efficient enough for self-replication. The factor that limits its polymerization efficiency is its weak sequence-independent binding of the primer/template substrate. Here we tested whether RNA polymerization could be improved by a cationic arginine cofactor, to improve the interaction with the substrate. In an RNA world, amino acid-nucleic acid conjugates could have facilitated the emergence of the translation apparatus and the transition to an RNP world. We chose the amino acid arginine for our study because this is the amino acid most adept to interact with RNA. An arginine cofactor was positioned at ten different sites on the ribozyme, using conjugates of arginine with short DNA or RNA oligonucleotides. However, polymerization efficiency was not increased in any of the ten positions. In five of the ten positions the arginine reduced or modulated polymerization efficiency, which gives insight into the substrate-binding site on the ribozyme. These results suggest that the existing polymerase ribozyme is not well suited to using an arginine cofactor