Approximating the Set of Local Minima in Partial RNA Folding Landscapes

Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes pL induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004).

Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O(n2D?ln?), where n is the sequence length, ? is the number of local minima in the partial landscape pL under consideration and D is the maximum number of steepest descent steps in attraction basins associated with pL.

Fast RNA conjugations on solid phase by strain-promoted cycloadditions

Strain promoted cycloaddition is presented as a tool for RNA conjugation on the solid phase; RNA-cyclooctyne conjugates are prepared by cycloaddition to both azide (strain-promoted azide-alkyne cycloaddition, SPAAC) and nitrile oxide dipoles (strain-promoted nitrile oxide-alkyne cycloaddition, SPNOAC). The conjugation is compatible with 2'-OMe blocks and with 2'-O-TBDMS protection on the ribose moieties of the sugar. Nitrile oxide dipoles are found to be more reactive click partners than azides. The conjugation proceeds within 10 min in aqueous solvents, at room temperature without any metal catalyst and tolerates dipoles of varying steric bulk and electronic demands, including pyrenyl, coumarin and dabcyl derivatives. © 2012 The Royal Society of Chemistry.

Detecting sequence dependent transcriptional pauses from RNA and protein number time series

Background: Evidence suggests that in prokaryotes sequence-dependent transcriptional pauses a?ect the dynamics of transcription and translation, as well as of small genetic circuits. So far, a few pause-prone sequences have been identi?ed from in vitro measurements of transcription elongation kinetics.

Results: Using a stochastic model of gene expression at the nucleotide and codon levels with realistic parameter values, we investigate three di?erent but related questions and present statistical methods for their analysis. First, we show that information from in vivo RNA and protein temporal numbers is su?cient to discriminate between models with and without a pause site in their coding sequence. Second, we demonstrate that it is possible to separate a large variety of models from each other with pauses of various durations and locations in the template by means of a hierarchical clustering and a random forest classi?er. Third, we introduce an approximate likelihood function that allows to estimate the location of a pause site.

Conclusions: This method can aid in detecting unknown pause-prone sequences from temporal measurements of RNA and protein numbers at a genome-wide scale and thus elucidate possible roles that these sequences play in the dynamics of genetic networks and phenotype.

Polysaccharide co-polymerases:the enigmatic conductors of the O-antigen assembly orchestra

The O-antigen lipopolysaccharides on bacterial surface contain variable number of oligosaccharide repeat units with their length having a modal distribution specific to the bacterial strain. The polysaccharide length distribution is controlled by the proteins called polysaccharide co-polymerases (PCPs), which are embedded in the inner membrane in Gram-negative bacteria and form homo oligomers. The 3D structures of periplasmic domains of several PCPs have been determined and provided the first insights into the possible mechanism of polysaccharide length determination mechanism. Here we review the current knowledge of structure and function of these polysaccharide length regulators.

RNA interference targeting cathepsin B of the carcinogenic liver fluke, Opisthorchis viverrini

Functional genomics have not been reported for Opisthorchis viverrini or the related fish-borne fluke, Clonorchis sinensis. Here we describe the introduction by square wave electroporation of Cy3-labeled small RNA into adult O. viverrini worms. Adult flukes were subjected to square wave electroporation employing a single pulse for 20 ms of 125V in the presence of 50 µg/ml of Cy3-siRNA. The parasites tolerated this manipulation and, at 24 and 48 h after electroporation, fluorescence from the Cy3-siRNA was evident throughout the parenchyma of the worms, with strong fluorescence evident in the guts and reproductive organs of the adult worms. Second, other worms were treated using the same electroporation settings with double stranded RNA targeting an endogenous papain-like cysteine protease, cathepsin B. This manipulation resulted in a significant reduction in specific mRNA levels encoding cathepsin B, and a significant reduction in cathepsin B activity against the diagnostic peptide, Z-Arg-Arg-AMC. This appears to be the first report of introduction of reporter genes into O. viverrini and the first report of experimental RNA interference (RNAi) in this fluke. The findings indicated the presence of an intact RNAi pathway in these parasites which, in turn, provides an opportunity to probe gene functions in this neglected tropical disease pathogen.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

Topoisomerase IIα promotes activation of RNA polymerase I transcription by facilitating pre-initiation complex formation

Type II DNA topoisomerases catalyse DNA double-strand cleavage, passage and re-ligation to effect topological changes. There is considerable interest in elucidating topoisomerase II roles, particularly as these proteins are targets for anti-cancer drugs. Here we uncover a role for topoisomerase IIa in RNA polymerase I-directed ribosomal RNA gene transcription, which drives cell growth and proliferation and is upregulated in cancer cells. Our data suggest that topoisomerase IIa is a component of the initiation-competent RNA polymerase Iß complex and interacts directly with RNA polymerase I-associated transcription factor RRN3, which targets the polymerase to promoter-bound SL1 in pre-initiation complex formation. In cells, activation of rDNA transcription is reduced by inhibition or depletion of topoisomerase II, and this is accompanied by reduced transient double-strand DNA cleavage in the rDNA-promoter region and reduced pre-initiation complex formation. We propose that topoisomerase IIa functions in RNA polymerase I transcription to produce topological changes at the rDNA promoter that facilitate efficient de novo pre-initiation complex formation.

Discovery of a Potent Boronic Acid Derived Inhibitor of the HCV RNA-Dependent RNA Polymerase

A boronic acid moiety was found to be a critical pharmacophore for enhanced in vitro potency against wild type hepatitis C replicons and known clinical polymorphic and resistant HCV mutant replicons. The synthesis, optimization, and structure-activity relationships associated with inhibition of HCV replication in a sub-genomic replication system for a series of non-nucleoside boron-containing HCV RNA-Dependent RNA Polymerase (NS5B) inhibitors are described. A summary of the discovery of GSK5852 (3), a molecule which entered clinical trials in subjects infected with HCV in 2011, is included.

Connectivity Mapping for Candidate Therapeutics Identification Using Next Generation Sequencing RNA-Seq Data

The advent of next generation sequencing technologies (NGS) has expanded the area of genomic research, offering high coverage and increased sensitivity over older microarray platforms. Although the current cost of next generation sequencing is still exceeding that of microarray approaches, the rapid advances in NGS will likely make it the platform of choice for future research in differential gene expression. Connectivity mapping is a procedure for examining the connections among diseases, genes and drugs by differential gene expression initially based on microarray technology, with which a large collection of compound-induced reference gene expression profiles have been accumulated. In this work, we aim to test the feasibility of incorporating NGS RNA-Seq data into the current connectivity mapping framework by utilizing the microarray based reference profiles and the construction of a differentially expressed gene signature from a NGS dataset. This would allow for the establishment of connections between the NGS gene signature and those microarray reference profiles, alleviating the associated incurring cost of re-creating drug profiles with NGS technology. We examined the connectivity mapping approach on a publicly available NGS dataset with androgen stimulation of LNCaP cells in order to extract candidate compounds that could inhibit the proliferative phenotype of LNCaP cells and to elucidate their potential in a laboratory setting. In addition, we also analyzed an independent microarray dataset of similar experimental settings. We found a high level of concordance between the top compounds identified using the gene signatures from the two datasets. The nicotine derivative cotinine was returned as the top candidate among the overlapping compounds with potential to suppress this proliferative phenotype. Subsequent lab experiments validated this connectivity mapping hit, showing that cotinine inhibits cell proliferation in an androgen dependent manner. Thus the results in this study suggest a promising prospect of integrating NGS data with connectivity mapping. © 2013 McArt et al.

The Arabidopsis RNA-binding protein FCA requires a lysine-specific demethylase 1 homolog to downregulate FLC.

A repressor of the transition to flowering in Arabidopsis is the MADS box protein FLOWERING LOCUS C (FLC). FCA, an RNA-binding protein, and FY, a homolog of the yeast RNA 3' processing factor Pfs2p, downregulate FLC expression and therefore promote flowering. FCA/FY physically interact and alter polyadenylation/3' processing to negatively autoregulate FCA. Here, we show that FCA requires FLOWERING LOCUS D (FLD), a homolog of the human lysine-specific demethylase 1 (LSD1) for FLC downregulation. FCA also partially depends on DICER-LIKE 3, involved in chromatin silencing. fca mutations increased levels of unspliced sense FLC transcript, altered processing of antisense FLC transcripts, and increased H3K4 dimethylation in the central region of FLC. These data support a close association of FCA and FLD in mediating H3K4 demethylation and thus transcriptional silencing of FLC and reveal roles for antisense RNA processing and DCL3 function in this regulation.

Small RNA-mediated chromatin silencing directed to the 3' region of the Arabidopsis gene encoding the developmental regulator, FLC.

Small RNA-mediated chromatin silencing is well characterized for repeated sequences and transposons, but its role in regulating single-copy endogenous genes is unclear. We have identified two small RNAs (30 and 24 nucleotides) corresponding to the reverse strand 3' to the canonical poly(A) site of FLOWERING LOCUS C (FLC), an Arabidopsis gene encoding a repressor of flowering. Genome searches suggest that these RNAs originate from the FLC locus in a genomic region lacking repeats. The 24-nt small RNA, which is most abundant in developing fruits, is absent in mutants defective in RNA polymerase IVa, RNA-DEPENDENT RNA POLYMERASE 2, and DICER-LIKE 3, components required for RNAi-mediated chromatin silencing. The corresponding genomic region shows histone 3 lysine 9 dimethylation, which was reduced in a dcl2,3,4 triple mutant. Investigations into the origins of the small RNAs revealed a polymerase IVa-dependent spliced, antisense transcript covering the 3' FLC region. Mutation of this genomic region by T-DNA insertion led to FLC misexpression and delayed flowering, suggesting that RNAi-mediated chromatin modification is an important component of endogenous pathways that function to suppress FLC expression.

Investigation of the mechanism of repression of rRNA synthesis by an anticancer drug

Ribosome biogenesis is a fundamental cellular process which is tightly regulated in normal cells. A number of tumour suppressors and oncogenes could affect the production of ribosomes at different levels and an upregulation could lead to increased protein biosynthesis which is one of the characteristic features of all cancer cells. Ribosome biogenesis is a very complex process which requires coordinated transcription by all three nucleolar polymerases and the first event in this process is synthesis of ribosomal RNA (rRNA) by RNA Polymerase I (Pol I). Importantly, recent data has pictured rRNA transcription as a key regulator of whole ribosome biogenesis and therefore makes it a valid and very attractive target for anticancer therapy, as well as a perspective biomarker. However, at the moment there is only one known specific inhibitor of Pol I transcription (at stage one of clinical trials) and this makes it very difficult for the development of drugs which would target rRNA transcription and consequently ribosome biogenesis. We have recently discovered that antitumor alkaloid ellipticine (isolated in 1959 from the plant species Ochrosia) is a potent inhibitor of Pol I transcription (both in vitro and in vivo). Ellipticine and its derivatives are known as efficient topoisomerase II inhibitors and inhibitors of some kinases, however we have shown that these inhibitory activities and the ability of ellipticine to repress Pol I activity are unrelated. Moreover, our preliminary data suggests that ellipticine specifically targets Pol I transcription and it has no effect on transcription by Pol II and Pol III at the same time scale. The possible mechanisms of inhibition of Pol I transcription by ellipticines will be discussed.