Dear all, I want to align several proteins, but here is the problem - I want to use existing knowledge of protein domains in each protein that needs to be multiply aligned

It seems like existing sophisticated alignment software like POA, A-Bruijn and ProDA still are blind to domain annotations. Perhaps I should mention that each domain can be considered an independent evolutionary unit, just as an amino acid in proteins, or nucleotide in DNA would be.

Is anyone here aware of any tools in PERL or BioPERL that can take in a multifasta file, the Pfam domain data for these same sequence, and return an alignment that is a combination of global alignment along with local domain-constrained alignment?

The final complication is that several of these proteins can have rearranged order of linear domains, so is there a way around it like I think A-Bruijn attempts to do? But this problem is less important than just the 1st task of aligning proteins constraining using domain annotation coordinates...

I realize that this might not be the right forum to ask heavy computational biology theory questions, but hey I am a newbie, so I am taking a stab at seeking your help anyways...

In reply to Protein multiple sequence alignment by onlyIDleft

Title:
Use:  <p> text here (a paragraph) </p>
and:  <code> code here </code>
to format your post, it's "PerlMonks-approved HTML":


  • Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
  • Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
  • Read Where should I post X? if you're not absolutely sure you're posting in the right place.
  • Please read these before you post! —
  • Posts may use any of the Perl Monks Approved HTML tags:
    a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
  • You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
            For:     Use:
    & &amp;
    < &lt;
    > &gt;
    [ &#91;
    ] &#93;
  • Link using PerlMonks shortcuts! What shortcuts can I use for linking?
  • See Writeup Formatting Tips and other pages linked from there for more info.