Question: Perturbation Of A Rna Secondary Structure
gravatar for Nicolas Rosewick
8.6 years ago by
Belgium, Brussels
Nicolas Rosewick9.0k wrote:


Is there a tool to pertub a secondary structure of a RNA by permutating its nucleotide sequence. In my case, I've a sequence encoding mutliple miRNAs (each pre-miR has a hairpin structure). I want to delete each of the hairpin structures by changing as little as possible the sequence.

The goal is to take this sequence encoding multiple miRNAs, to modify it to perturb the structure to delete the hairpins and to put it in a lentiviral vector and to see if it is processed

I posted a more biological question for the same topic on biology.stackexchange :

Thanks a lot,


mirna • 1.8k views
ADD COMMENTlink modified 8.6 years ago by Ib84100 • written 8.6 years ago by Nicolas Rosewick9.0k

nobody ? I'll try to pertub only the 3' part of the stemloop (so the miRNA star) to delete the complementarity between the miR-5p and miR-3p . Anyone has an another idea to do this ?


ADD REPLYlink written 8.6 years ago by Nicolas Rosewick9.0k
gravatar for Alastair Kerr
8.6 years ago by
Alastair Kerr5.2k
Manchester/UK/Cancer Biomarker Centre at CRUK-MI
Alastair Kerr5.2k wrote:

I am unaware of a package that will perform that task without some intervention of your own.

I suggest that the best approach at the moment is to use tools from the Vienna RNA Package along with some scripting (bioperl, biopython...) to systematically change that sequence.

My workflow would be:

  1. Create mutants using scripting language or manually
  2. Use RNAfold* to create the 2d structure
  3. Use RNAdistance* to compare secondary structure to original.
  4. Sort all mutations by (3) and judge yourself where you should make the cut-off

(*both from the Vienna RNA package)

Note that RNAfold uses Watson-Crick base pairing only. There are other folding programs in development see the wikiomics page but I have yet to try one better then RNAfold for miRNAs.

ADD COMMENTlink written 8.6 years ago by Alastair Kerr5.2k

Yes that's a good start. I'll add a step at the end: the scrambled hairpin sequence must have a score under 17 with sRNAloop. Thanks

ADD REPLYlink written 8.6 years ago by Nicolas Rosewick9.0k

It works well with RNAdistance ;)

ADD REPLYlink written 8.6 years ago by Nicolas Rosewick9.0k
gravatar for Ib84
8.0 years ago by
Ib84100 wrote:

I'm not sure if it's what you are after, but the GeneGA tool does permutate sequences using a genetic / evolutionary algorithm with an objective function which maximizes free energy. Hence it strives to decrease secondary structure, but not sure if this is suited well to explicitely eliminate hairpins. Also, afaik, it is restricted to synomymous mutations, i.e. resulting protein sequence remains identical. Its purpose is to improve heterologeous expression by optimizing both codon usage and RNA secondary structure, with respect to recent findings of the so-called "ramp" (Tuller 2010, Science). Recommend to use with rstudio.

ADD COMMENTlink modified 8.0 years ago • written 8.0 years ago by Ib84100
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1492 users visited in the last hour