Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • How do matrix weights for amino acid 'X' affect protein alignments in MUSCLE?

    Sorry, this is a cross-post from biostars, but I haven't had any luck there in figuring this out...

    I'm developing a protein weight matrix from scratch to use with MUSCLE. I started with the 20 amino acids and the stop character, but using just those, I get a warning from MUSCLE that makes no sense^. But if I add a row & column for the amino acid code 'X' (with all 0s - everything else ranges from 1-25), the warning goes away. However, even though the sequences I'm aligning have no occurrences of 'X', the alignments produced using the matrix without 'X' and the matrix with 'X' are different. I don't understand why that is.

    How do alignment algorithms treat 'X' in the weight matrix when there are no X's in the sequences being aligned? Why would an alignment produced using a matrix with 'X' versus an alifgnment produced using a matrix without 'X' be different? Are there weights I can insert for 'X' in the matrix that would not affect the alignment?

    Side question - I don't have 'B' or 'Z' in my matrix either (nor do they occur in my sequences). Do they need to be in the matrix for the alignment software to yield good results?

    ^ [I]The warning from muscle is "*** WARNING *** Matrix is not symmetrical, ?->?=5, ?->?=0". This doesn't make sense for 3 reasons: 1. the matrix is symmetrical. 2. There were no 0s in my matrix when the warning occurred. 3. Simply adding a row/column for 'X' silenced the warning. I looked briefly at the code in muscle where the warning is generated and it's simply checking that
    Code:
    matrix[i][j] == matrix[j][i][/i]
    . The fact that there's an item that doesn't match is unsettling. Perhaps this is a hidden requirement that 'X' must be present in the matrix - but still, why would the alignment change with/without it?

  • #2
    The answer is, the presence or absence of 'X' in a protein weight matrix should not affect the alignment of sequences not containing 'X', however there is a bug in muscle that is fixed by the following patch. Supplying muscle a custom weight matrix without any of the ambiguous amino acids causes the alignment to differ because muscle was not setting the weights in the last column and that was causing the warnings.

    https://github.com/Bioconductor-mirror/muscle/compare/master...hepcat72atch-2

    Originally posted by hepcat72 View Post
    Sorry, this is a cross-post from biostars, but I haven't had any luck there in figuring this out...

    I'm developing a protein weight matrix from scratch to use with MUSCLE. I started with the 20 amino acids and the stop character, but using just those, I get a warning from MUSCLE that makes no sense^. But if I add a row & column for the amino acid code 'X' (with all 0s - everything else ranges from 1-25), the warning goes away. However, even though the sequences I'm aligning have no occurrences of 'X', the alignments produced using the matrix without 'X' and the matrix with 'X' are different. I don't understand why that is.

    How do alignment algorithms treat 'X' in the weight matrix when there are no X's in the sequences being aligned? Why would an alignment produced using a matrix with 'X' versus an alifgnment produced using a matrix without 'X' be different? Are there weights I can insert for 'X' in the matrix that would not affect the alignment?

    Side question - I don't have 'B' or 'Z' in my matrix either (nor do they occur in my sequences). Do they need to be in the matrix for the alignment software to yield good results?

    ^ [I]The warning from muscle is "*** WARNING *** Matrix is not symmetrical, ?->?=5, ?->?=0". This doesn't make sense for 3 reasons: 1. the matrix is symmetrical. 2. There were no 0s in my matrix when the warning occurred. 3. Simply adding a row/column for 'X' silenced the warning. I looked briefly at the code in muscle where the warning is generated and it's simply checking that
    Code:
    matrix[i][j] == matrix[j][i][/i]
    . The fact that there's an item that doesn't match is unsettling. Perhaps this is a hidden requirement that 'X' must be present in the matrix - but still, why would the alignment change with/without it?

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    9 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    49 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    67 views
    0 likes
    Last Post seqadmin  
    Working...
    X