You are here

snugdock run error: Sequence for partner1 does not match first member of ensemble1

1 post / 0 new
snugdock run error: Sequence for partner1 does not match first member of ensemble1
#1

Hello everyone,

I am trying to use Rosetta for antibody homology modeling and docking.  
The method I used was from article:"Modeling and docking of antibody structures with Rosetta" (https://www.nature.com/articles/nprot.2016.180
The Rosetta version I used is:  3.13  Release date: Tuesday, June 1, 2021
 
I performed antibody homology modeling and docking preparation according to the method in the article.  When I run snugDock, the software keeps saying something wrong:  

[ ERROR ] UtilityExitException
ERROR: Sequence for partner1:
EVQLVESGGGLVQPGRSLRLSCAASGFTFDDYAMHWVRQAPGKGLEWVSAITWNSGHIDYADSVEGRFTISRDNAKNSLYLQMNSLRAEDTAVYYCAKVSYLSTASSLDYWGQGTLVTVSDIQMTQSPSSLSASVGDRVTITCRASQGIRNYLAWYQQKPGKAPKLLIYAASTLQSGVPSRFSGSGSGTDFTLTISSLQPEDVATYYCQRYNRAPYTFGQGTKVEIKRT
 does not match first member of ensemble1:
DIQMTQSPSSLSASVGDRVTITCRASQGIRNYLAWYQQKPGKAPKLLIYAASTLQSGVPSRFSGSGSGTDFTLTISSLQPEDVATYYCQRYNRAPYTFGQGTKVEIKRTEVQLVESGGGLVQPGRSLRLSCAASGFTFDDYAMHWVRQAPGKGLEWVSAITWNSGHIDYADSVEGRFTISRDNAKNSLYLQMNSLRAEDTAVYYCAKVSYLSTASSLDYWGQGTLVTVS


AN INTERNAL ERROR HAS OCCURED. PLEASE SEE THE CONTENTS OF ROSETTA_CRASH.log FOR DETAILS.

 


The code I am executing:

docking_prepack_protocol.mpi.linuxgccrelease  \
    -in:file:s   model-0.relaxed_0168.antibody_antigen_start_new.pdb   \
    -ex1    \
    -ex2    \
    -partners  LH_A  \
    -ensemble1  antibody_ensemble.1.list   \
    -ensemble2  antigen_ensemble.list    \
    -docking:dock_rtmin

Based on the log message, I thought there might be a problem with the chain order in my partner1 file or ensemble1. But after checking, I found that the antibody sequence in the ppk file of ensemble1 is fine, the order is: L chain->H chain. in partner1.ppk file, the order of the chain is fine, also L chain->H chain, but the identifier is changed to A and B, which may be handled by the software docking_prepack_protocol.
When I searched for related questions in the forum, I found that the same question was mentioned before by
https://www.rosettacommons.org/node/10396

So I tried to swap the order of L chain and H chain in the partner1.ppk file and changed the chain identifiers to H and L. The rest of the file remained unchanged.
I reran the above command and it works!

However I'm even more confused!

My questions:
1, why does snugdock think my sequence is not match?
2, Can I get normal results with the way I'm running it now? Why I set "-partners LH_A" in the parameter and the sequence order of partner1.ppk is H chain->L chain, but snugdock works well?
3, why the sequence identifier in partner1.ppk be changed to A and B after prepack run?

Category: 
Post Situation: 
Fri, 2022-03-25 01:21
pagumaaaaa