How ViennaRNA Package deal with character other than AUGC in RNA secondary structure prediction?
1
0
Entering edit mode
7.2 years ago
wyj.jackson ▴ 20

I have a RNA sequence as following, where "o" connects two RNA sequences.

UAGGGoooooooAGUU

Most of the RNA secondary structure prediction tools do not allow characters other than AUGC, but I find out the ViennaRNA allows "o" and can predict structure with it. But I am wondering how ViennaRNA (RNA2Dfold) deal with "o". I can not find any clues on ViennaRNA website. Please help me if you know the answer:) Thanks in advance.

RNA-Seq sequence • 1.7k views
ADD COMMENT
2
Entering edit mode
7.2 years ago
Martombo ★ 3.0k

non AUGC bases are considered part of the RNA molecule, but they cannot form any pair. they can be in a loop separating two matching stretches, like in the following example:

AAAAAAAAAAooooooooTTTTTTTTTT
((((((((((........))))))))))

which is not the same as if they were completely ignored:

AAAAAAAAAATTTTTTTTTT
((((((((....))))))))

can I ask you what do these o represent in your sequence?

ADD COMMENT
0
Entering edit mode

Thanks for your reply.

Do you know how those non AUGC bases influence the energy model of hairpin, interior and multi loops in the package? Also treat them as unpaired bases?

For "o"s, our collaborator said that they are physically very flexible and can link two RNA sequences.

ADD REPLY
0
Entering edit mode

according to how your collaborators described it, it looks like the ViennaRNA package deals with them as you want it to. They are in all respects like unpaired bases.

ADD REPLY

Login before adding your answer.

Traffic: 2509 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6