You are here

flexpep docking problem

6 posts / 0 new
Last post
flexpep docking problem

Dear All:
I am very new to the peptide flexible docking. When I use the locally installed rosetta 3.2.1 to perform a peptide docking. I got the error as :

ERROR: total_residue() != 0
ERROR:: Exit from: src/core/pose/ line: 1369

The command line is :
FlexPepDocking.linuxgccrelease -database /usr/local/src/rosetta_database/ -s renum.pdb -flexpep_prepack > flex.log

I also attached the pdb file and log file.

Hope you can help me out! thank you so much!!

Best wishes,

flex.log_.txt5.8 KB
renum.pdb_.txt99.14 KB
Post Situation: 
Mon, 2011-04-18 17:15

I am having the same problem. The line refrred to in the error talks about the occupancy entries in the pdb file, but it doesn't help when I set them all to 1.00.
command: /home/rosetta-3.2/rosetta_source/bin/FlexPepDocking.linuxgccrelease -s input.pdb -database /home/rosetta-3.2/rosetta_database/ -flexpep_prepack -ex1 -ex2aro
An idea of what the problem is would be appreciated.

EDIT: OK it works so long as the native structure is included in the command

Tue, 2011-04-19 03:36

Dear pardave:
Could you kindly show me the command line and the input file(pdbs, native structure, flags etc.). So I can find the problems in my case. Thank you so much!!

Tue, 2011-04-19 18:32

Dear All:
by following the comment of pardave, I add the -native renum.pdb in my command. It goes well and finished the calculation successfully. So in summary:
1, use some script to renumber your protein, best way is start from residue number 1. ( I don't know this is the real issue. But If I use the discontinued residue numbers, the flexpep_docking program crashed with the segment fault error.)
2, use the command including the -native *.pdb ( can be same as the -s *.pdb)

Hoping this is useful to you.

Best wishes,

Tue, 2011-04-19 18:41

Dear Michael and Dave,

The problem is definitely on the account of a missing "-native" structure. As repeatedly noted in the documentation the -native flag is required by FlexPepDock. It is used for RMSD calculations, as michael mentioned you can input the starting structure (same one you use for '-s') or a different reference structure (like the actual solution, or a different model, etc...)

Regrading discontinuous numbering, this should pose no problem, Michael, if you can reproduce the scenario please send us the cmd-line and PDB and we'll check it out.


Wed, 2011-04-20 02:21

Dear Nir:
Yes, you are right. I check the structure I used before and found some occupancy column is 0.00. So rosetta recognize as no atoms in residue. Then it printed out many missing atoms error, and finally it crashed. So in general, discontinuous numbers of residue DID NOT have problems.
Thanks all the discussions.

Best wishes,


Wed, 2011-04-20 22:06