Design principles for riboswitch function

Scientific and technological advances that enable the tuning of integrated regulatory components to match network and system requirements are critical to reliably control the function of biological systems. RNA provides a promising building block for the construction of tunable regulatory components based on its rich regulatory capacity and our current understanding of the sequence–function relationship. One prominent example of RNA-based regulatory components is riboswitches, genetic elements that mediate ligand control of gene expression through diverse regulatory mechanisms. While characterization of natural and synthetic riboswitches has revealed that riboswitch function can be modulated through sequence alteration, no quantitative frameworks exist to investigate or guide riboswitch tuning. Here, we combined mathematical modeling and experimental approaches to investigate the relationship between riboswitch function and performance. Model results demonstrated that the competition between reversible and irreversible rate constants dictates performance for different regulatory mechanisms. We also found that practical system restrictions, such as an upper limit on ligand concentration, can significantly alter the requirements for riboswitch performance, necessitating alternative tuning strategies. Previous experimental data for natural and synthetic riboswitches as well as experiments conducted in this work support model predictions. From our results, we developed a set of general design principles for synthetic riboswitches. Our results also provide a foundation from which to investigate how natural riboswitches are tuned to meet systems-level regulatory demands

Transforaminal endoscopic surgery for symptomatic lumbar disc herniations: a systematic review of the literature

The study design includes a systematic literature review. The objective of the study was to evaluate the effectiveness of transforaminal endoscopic surgery and to compare this with open microdiscectomy in patients with symptomatic lumbar disc herniations. Transforaminal endoscopic techniques for patients with symptomatic lumbar disc herniations have become increasingly popular. The literature has not yet been systematically reviewed. A comprehensive systematic literature search of the MEDLINE and EMBASE databases was performed up to May 2008. Two reviewers independently checked all retrieved titles and abstracts and relevant full text articles for inclusion criteria. Included articles were assessed for quality and outcomes were extracted by the two reviewers independently. One randomized controlled trial, 7 non-randomized controlled trials and 31 observational studies were identified. Studies were heterogeneous regarding patient selection, indications, operation techniques, follow-up period and outcome measures and the methodological quality of these studies was poor. The eight trials did not find any statistically significant differences in leg pain reduction between the transforaminal endoscopic surgery group (89%) and the open microdiscectomy group (87%); overall improvement (84 vs. 78%), re-operation rate (6.8 vs. 4.7%) and complication rate (1.5 vs. 1%), respectively. In conclusion, current evidence on the effectiveness of transforaminal endoscopic surgery is poor and does not provide valid information to either support or refute using this type of surgery in patients with symptomatic lumbar disc herniations. High-quality randomized controlled trials with sufficiently large sample sizes are direly needed to evaluate if transforaminal endoscopic surgery is more effective than open microdiscectomy

Is Thermosensing Property of RNA Thermometers Unique?

A large number of studies have been dedicated to identify the structural and sequence based features of RNA thermometers, mRNAs that regulate their translation initiation rate with temperature. It has been shown that the melting of the ribosome-binding site (RBS) plays a prominent role in this thermosensing process. However, little is known as to how widespread this melting phenomenon is as earlier studies on the subject have worked with a small sample of known RNA thermometers. We have developed a novel method of studying the melting of RNAs with temperature by computationally sampling the distribution of the RNA structures at various temperatures using the RNA folding software Vienna. In this study, we compared the thermosensing property of 100 randomly selected mRNAs and three well known thermometers - rpoH, ibpA and agsA sequences from E. coli. We also compared the rpoH sequences from 81 mesophilic proteobacteria. Although both rpoH and ibpA show a higher rate of melting at their RBS compared with the mean of non-thermometers, contrary to our expectations these higher rates are not significant. Surprisingly, we also do not find any significant differences between rpoH thermometers from other -proteobacteria and E. coli non-thermometers

Efficient Algorithms for Probing the RNA Mutation Landscape

The diversity and importance of the role played by RNAs in the regulation and development of the cell are now well-known and well-documented. This broad range of functions is achieved through specific structures that have been (presumably) optimized through evolution. State-of-the-art methods, such as McCaskill's algorithm, use a statistical mechanics framework based on the computation of the partition function over the canonical ensemble of all possible secondary structures on a given sequence. Although secondary structure predictions from thermodynamics-based algorithms are not as accurate as methods employing comparative genomics, the former methods are the only available tools to investigate novel RNAs, such as the many RNAs of unknown function recently reported by the ENCODE consortium. In this paper, we generalize the McCaskill partition function algorithm to sum over the grand canonical ensemble of all secondary structures of all mutants of the given sequence. Specifically, our new program, RNAmutants, simultaneously computes for each integer k the minimum free energy structure MFE(k) and the partition function Z(k) over all secondary structures of all k-point mutants, even allowing the user to specify certain positions required not to mutate and certain positions required to base-pair or remain unpaired. This technically important extension allows us to study the resilience of an RNA molecule to pointwise mutations. By computing the mutation profile of a sequence, a novel graphical representation of the mutational tendency of nucleotide positions, we analyze the deleterious nature of mutating specific nucleotide positions or groups of positions. We have successfully applied RNAmutants to investigate deleterious mutations (mutations that radically modify the secondary structure) in the Hepatitis C virus cis-acting replication element and to evaluate the evolutionary pressure applied on different regions of the HIV trans-activation response element. In particular, we show qualitative agreement between published Hepatitis C and HIV experimental mutagenesis studies and our analysis of deleterious mutations using RNAmutants. Our work also predicts other deleterious mutations, which could be verified experimentally. Finally, we provide evidence that the 3′ UTR of the GB RNA virus C has been optimized to preserve evolutionarily conserved stem regions from a deleterious effect of pointwise mutations. We hope that there will be long-term potential applications of RNAmutants in de novo RNA design and drug design against RNA viruses. This work also suggests potential applications for large-scale exploration of the RNA sequence-structure network. Binary distributions are available at http://RNAmutants.csail.mit.edu/

Overexpression of miR-128 specifically inhibits the truncated isoform of NTRK3 and upregulates BCL2 in SH-SY5Y neuroblastoma cells

<p>Abstract</p> <p>Background</p> <p>Neurotrophins and their receptors are key molecules in the regulation of neuronal differentiation and survival. They mediate the survival of neurons during development and adulthood and are implicated in synaptic plasticity. The human neurotrophin-3 receptor gene <it>NTRK3 </it>yields two major isoforms, a full-length kinase-active form and a truncated non-catalytic form, which activates a specific pathway affecting membrane remodeling and cytoskeletal reorganization. The two variants present non-overlapping 3'UTRs, indicating that they might be differentially regulated at the post-transcriptional level. Here, we provide evidence that the two isoforms of <it>NTRK3 </it>are targeted by different sets of microRNAs, small non-coding RNAs that play an important regulatory role in the nervous system.</p> <p>Results</p> <p>We identify one microRNA (miR-151-3p) that represses the full-length isoform of <it>NTRK3 </it>and four microRNAs (miR-128, miR-485-3p, miR-765 and miR-768-5p) that repress the truncated isoform. In particular, we show that the overexpression of miR-128 - a brain enriched miRNA - causes morphological changes in SH-SY5Y neuroblastoma cells similar to those observed using an siRNA specifically directed against truncated <it>NTRK3</it>, as well as a significant increase in cell number. Accordingly, transcriptome analysis of cells transfected with miR-128 revealed an alteration of the expression of genes implicated in cytoskeletal organization as well as genes involved in apoptosis, cell survival and proliferation, including the anti-apoptotic factor <it>BCL2</it>.</p> <p>Conclusions</p> <p>Our results show that the regulation of <it>NTRK3 </it>by microRNAs is isoform-specific and suggest that neurotrophin-mediated processes are strongly linked to microRNA-dependent mechanisms. In addition, these findings open new perspectives for the study of the physiological role of miR-128 and its possible involvement in cell death/survival processes.</p

Maximum expected accuracy structural neighbors of an RNA secondary structure

International audienceBACKGROUND: Since RNA molecules regulate genes and control alternative splicing by allostery, it is important to develop algorithms to predict RNA conformational switches. Some tools, such as paRNAss, RNAshapes and RNAbor, can be used to predict potential conformational switches; nevertheless, no existent tool can detect general (i.e., not family specific) entire riboswitches (both aptamer and expression platform) with accuracy. Thus, the development of additional algorithms to detect conformational switches seems important, especially since the difference in free energy between the two metastable secondary structures may be as large as 15-20 kcal/mol. It has recently emerged that RNA secondary structure can be more accurately predicted by computing the maximum expected accuracy (MEA) structure, rather than the minimum free energy (MFE) structure. RESULTS: Given an arbitrary RNA secondary structure S₀ for an RNA nucleotide sequence a = a₁,..., a(n), we say that another secondary structure S of a is a k-neighbor of S₀, if the base pair distance between S₀ and S is k. In this paper, we prove that the Boltzmann probability of all k-neighbors of the minimum free energy structure S₀ can be approximated with accuracy ε and confidence 1 - p, simultaneously for all 0 ≤ k N(ε,p,K)=Φ⁻¹(p/2K)²/4ε², where Φ(z) is the cumulative distribution function (CDF) for the standard normal distribution. We go on to describe the algorithm RNAborMEA, which for an arbitrary initial structure S₀ and for all values 0 ≤ k < K, computes the secondary structure MEA(k), having maximum expected accuracy over all k-neighbors of S₀. Computation time is O(n³ * K²), and memory requirements are O(n² * K). We analyze a sample TPP riboswitch, and apply our algorithm to the class of purine riboswitches. CONCLUSIONS: The approximation of RNAbor by sampling, with rigorous bound on accuracy, together with the computation of maximum expected accuracy k-neighbors by RNAborMEA, provide additional tools toward conformational switch detection. Results from RNAborMEA are quite distinct from other tools, such as RNAbor, RNAshapes and paRNAss, hence may provide orthogonal information when looking for suboptimal structures or conformational switches. Source code for RNAborMEA can be downloaded from http://sourceforge.net/projects/rnabormea/ or http://bioinformatics.bc.edu/clotelab/RNAborMEA/

Using dissolved H₂O in rhyolitic glasses to estimate palaeo-ice thickness during a subglacial eruption at Bláhnúkur(Torfajökull, Iceland)

The last decade has seen the refinement of a technique for reconstructing palaeo-ice thicknesses based on using the retained H2O and CO2 content in glassy eruptive deposits to infer quenching pressures and therefore ice thicknesses. The method is here applied to Bláhnúkur, a subglacially erupted rhyolitic edifice in Iceland. A decrease in water content from ~0.7 wt.% at the base to ~0.3 wt.% at the top of the edifice suggests that the ice was 400 m thick at the time of the eruption. As Bláhnúkur rises 350 m above the surrounding terrain, this implies that the eruption occurred entirely within ice, which corroborates evidence obtained from earlier lithofacies studies. This paper presents the largest data set (40 samples) so far obtained for the retained volatile contents of deposits from a subglacial eruption. An important consequence is that it enables subtle but significant variations in water content to become evident. In particular, there are anomalous samples which are either water-rich (up to 1 wt.%) or water-poor (~0.2 wt.%), with the former being interpreted as forming intrusively within hyaloclastite and the latter representing batches of magma that were volatile-poor prior to eruption. The large data set also provides further insights into the strengths and weaknesses of using volatiles to infer palaeo-ice thicknesses and highlights many of the uncertainties involved. By using examples from Bláhnúkur, the quantitative use of this technique is evaluated. However, the relative pressure conditions which have shed light on Bláhnúkur’s eruption mechanisms and syn-eruptive glacier response show that, despite uncertainties in absolute values, the volatile approach can provide useful insight into the mechanisms of subglacial rhyolitic eruptions, which have never been observed

An accurate and interpretable model for siRNA efficacy prediction

BACKGROUND: The use of exogenous small interfering RNAs (siRNAs) for gene silencing has quickly become a widespread molecular tool providing a powerful means for gene functional study and new drug target identification. Although considerable progress has been made recently in understanding how the RNAi pathway mediates gene silencing, the design of potent siRNAs remains challenging. RESULTS: We propose a simple linear model combining basic features of siRNA sequences for siRNA efficacy prediction. Trained and tested on a large dataset of siRNA sequences made recently available, it performs as well as more complex state-of-the-art models in terms of potency prediction accuracy, with the advantage of being directly interpretable. The analysis of this linear model allows us to detect and quantify the effect of nucleotide preferences at particular positions, including previously known and new observations. We also detect and quantify a strong propensity of potent siRNAs to contain short asymmetric motifs in their sequence, and show that, surprisingly, these motifs alone contain at least as much relevant information for potency prediction as the nucleotide preferences for particular positions. CONCLUSION: The model proposed for prediction of siRNA potency is as accurate as a state-of-the-art nonlinear model and is easily interpretable in terms of biological features. It is freely available on the web a

Identification of Continuous Human B-Cell Epitopes in the Envelope Glycoprotein of Dengue Virus Type 3 (DENV-3)

BACKGROUND:Dengue virus infection is a growing global public health concern in tropical and subtropical regions of the world. Dengue vaccine development has been hampered by concerns that cross-reactive immunological memory elicited by a candidate vaccine could increase the risk of development of more severe clinical forms. One possible strategy to reduce risks associated with a dengue vaccine is the development of a vaccine composed of selected critical epitopes of each of the serotypes. METHODOLOGY/PRINCIPAL FINDINGS:Synthetic peptides were used to identify B-cell epitopes in the envelope (E) glycoprotein of dengue virus type 3 (DENV-3). Eleven linear, immunodominant epitopes distributed in five regions at amino acid (aa) positions: 51-65, 71-90, 131-170, 196-210 and 246-260 were identified by employing an enzyme- linked immunosorbent assay (ELISA), using a pool of human sera from dengue type 3 infected individuals. Peptides 11 (aa51-65), 27 and 28 (aa131-150) also reacted with dengue 1 (DENV-1) and dengue 2 (DENV-2) patient sera as analyzed through the ROC curves generated for each peptide by ELISA and might have serotype specific diagnostic potential. Mice immunized against each one of the five immunogenic regions showed epitopes 51-65, 131-170, 196-210 and 246-260 elicited the highest antibody response and epitopes131-170, 196-210 and 246-260, elicited IFN-gamma production and T CD4+ cell response, as evaluated by ELISA and ELISPOT assays respectively. CONCLUSIONS/SIGNIFICANCE:Our study identified several useful immunodominant IgG-specific epitopes on the envelope of DENV-3. They are important tools for understanding the mechanisms involved in antibody dependent enhancement and immunity. If proven protective and safe, in conjunction with others well-documented epitopes, they might be included into a candidate epitope-based vaccine