start pose and native pose don’t match in lengh

Member Site Forums Rosetta 3 Rosetta 3 – General start pose and native pose don’t match in lengh

Viewing 1 reply thread
  • Author
    Posts
    • #1118
      Anonymous

        Hi,

        I want to use a crystal structure with a shorter protein sequence as my native structure during loop re-modelling so the rmsd of generated structures will be aligned to it. But there is always an error:

        ERROR: Start pose and native pose don’t match in length
        ERROR:: Exit from: src/protocols/loops/LoopRelaxMover.cc line: 222

        How to force to align the same region of the two structures?

        Many thanks.

      • #6420
        Anonymous

          I don’t think it’s possible with stock code. Rosetta is not performing a sequence alignment before the structural alignment; it assumes it is comparing copies of the same structure. If the pose lengths don’t match, it doesn’t know what atoms to ignore in RMSD calculation.

          The simple option is to edit your native pose so that it will sequence-align with your result poses; whether or not you can do this depends on the problem.

          The complex solution will be to rewire the code to take an alignment file between your pose and your native pose and align based on that.

          I think one of the other executables may let you do RMSD calculation with explicit alignment files easily, so you could postcalculate RMSDs? It looks like super_aln has this functionality, but the RMSD part is commented out, so I don’t know if it works or not.

      Viewing 1 reply thread
      • You must be logged in to reply to this topic.