I have a profile-profile alignment and I want to get the per column sequence similarity. That is, for every aligned column pair I need to get a score.
Do you know any program or programming library that can take already aligned profiles and calculate the per column scores (using any reasonable scoring function)?
•
written
by
Jan Kosinski
From the discussion in the comments, you want to calculate a similarity score between two columns which are aligned to each other.
A column can be represented as a vector of 20 values, i.e. the number of each amino acid in the column. So you can see it as calculate the distance between two vectors. It can be easily implemented as a first score.
Yes, and I think I could it quite easily with PyCogent. I may try. Although it would be extremely simple score, for example it does not take into account amino acid similarities.
Yes yes, of course, the point is that I feel like re-inventing the wheel, scores for comparing columns from two profiles have been improving for more than a decade and now, and they are quite accurate now. The equations are there, but they are quite complicated, so that's why I look for ready to use solution ;-)
The Matlab bioinformatics toolbox can do this easily. you can read in an alignment in virtually any format using multialignread. Then you can use seqconsensus, seqlogo or seqprofile to do almost anything you might need.
No no, I don't want to get scores of columns of a single multiple sequence alignment. I have TWO multiple sequence alignments aligned to each other, like from typical profile/HMM-profile/HMM comparison (e.g. HHSearch). Now I need to calculate the SIMILARITY SCORES OF ALIGNING EVERY ALIGNED PAIR OF COLUMNS.
Such scoring functions are listed for instance here (Section Scoring functions of Materials and Methods): http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2279992/
I don't think so. Lets say I have 1000 sequences in alignment A, and 50 sequences in alignment B. Lets say that a column from A with invariant Lys residue has been aligned with a column from B with invariant Val. The score should be low, but the score for the "combined" column would be high.
Log in to answer this question.
Do you have the alignment in fasta format ? If so, you can use MstatX
isn't COMPASS downloadable as a standalone? http://www.soton.ac.uk/~re1u06/software/compass/index.html
I can convert it to any format and MstaX looks cool, but again, it calculates scores for ONE multiple sequence alignment. I have TWO, aligned as profile to profile (definition of profile-profile alignment is e.g. here http://phylogenomics.berkeley.edu/profilealignment/).
Many programs like COMPASS for profile-profile or HMM-HMM comparison are downloadable, but I cannot find any that 1) takes existing alignment of two profiles, instead of aligning them by itself, and 2) returns the local scores
Jeremy, why it has been tagged with "sequence-logo". I don't look for sequence logos for profiles, and there is a tool for that http://www.sanger.ac.uk/cgi-bin/software/analysis/logomat-p.cgi
And what does the "PWM" tag mean?