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?
•
written
by
virpatel3 •
2 answers
OK I tried on the website clustalW
It seems to work
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
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):
Edit: Here's the paper: "Seeking an ancient enzyme in Methanococcus jannaschii using orf, a program based on predicted secondary structure comparisons"
Log in to answer this question.
-
How to represent a summary of disagreement between two secondary structure predictors given predict…
written by user366312 •
-
Clustal W multi sequence alignment with secondary structure constraints
written by JJP •
-
alignment of secondary structure of template and target
written by amara.jabeen86 •
-
Well characterized circular RNA structure
written by ccag •
-
Consensus secondary RNA stractures
written by blur
-
PSSPred Alignment Alignment
written by virpatel3 •
-
RaptorX Alignment Method
written by virpatel3 •
-
protein inverse folding
written by m.movahedi •
-
GenBank Batch Download
written by virpatel3 •
-
Where Is The Experimental Rna Structure Data?
written by Veewhy •
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.
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.
You could use Biopython's alignment module pairwise2 with a custom scoring function. Documentation