Question: Secondary Structure Alignment Method
2
gravatar for virpatel3
5.7 years ago by
virpatel320
United States
virpatel320 wrote:

I have many genes for many bacteriophage with precomputed secondary structures such as:

 

HHHHHHHHHHHHHCCCCCCHHHHHHHHHHHHHHHHHHHHHCCHHHHHHHHHCHHHHHHHH

I want to align each gene's secondary structure with the SSE (secondary structure elements) of the other bacteriophages. So for instance, I may also have

 

HHHHHHHHHHHHHCCCCCCHHHHHHHHHHHHHHHHHHHHHCCHHHHHHHHHCCHHHHHHH 

for another gene. I would like to be able to align these genes so that I can find the evolutionary relationships between the bacteriophages. 

Does anyone have a suggestion for the algorithms I should use for this project?

 

secondary structure alignment • 2.2k views
ADD COMMENTlink modified 14 months ago by The160 • written 5.7 years ago by virpatel320
1

I never faced the problem of aligning secondary structure so I'm probably overlooking a lot of complications, but by the way the question is asked the first and easiest thing that comes to my mind is to align one set of genes vs the other with BLAST.

ADD REPLYlink written 5.7 years ago by dariober11k
1

Blast is good but him ask for the alignment of the secondary structures and i think that blast is not good in this case. yes is true that the blast is a good on line tools for alignment but the problem is that blast compute the alignment on the basis of a database searching the possible releated sequences.

in your case virpatel3 i'm suggest to create an algorithm of alignment.

you can search on google examples about this.

sure i never heard about the alignment of SS elements.

ADD REPLYlink modified 5.7 years ago • written 5.7 years ago by a.polo88110

You could use Biopython's alignment module pairwise2 with a custom scoring function. Documentation

ADD REPLYlink written 14 months ago by cschu1812.5k
1
gravatar for a.polo88
5.7 years ago by
a.polo88110
Italy
a.polo88110 wrote:

ok i tried on the website clustalW

it seems to work 

http://www.ebi.ac.uk/Tools/services/web/toolresult.ebi?jobId=clustalw2-I20150222-140308-0817-61234194-pg

try it!

if you will use it remember that fro each of sequences the format is 

> name of protein or what do you want

sequence

example :

>  name

HCCCHHHC

>name1

HCCCHHHHHHHC

results of your examples

http://www.ebi.ac.uk/Tools/services/web/toolresult.ebi?jobId=clustalw2-I20150222-140838-0639-99240772-oy

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by a.polo88110
0
gravatar for The
14 months ago by
The160
United States
The160 wrote:

If I correctly remember , there was a paper in PNAS by George D Rose in late 90s which dealt with this. The goal was to find out structurally similar sequences based on predicted secondary structure. So far I remember they employed a simple scoring system (match/mismatch binary type): enter image description here

Edit: Here's the paper: "Seeking an ancient enzyme in Methanococcus jannaschii using orf, a program based on predicted secondary structure comparisons"

https://www.pnas.org/content/95/6/2818

ADD COMMENTlink modified 14 months ago • written 14 months ago by The160
Please log in to add an answer.

Help
Access

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