I am facing the error during protein peptide docking in haddock V3. My protein contains two domains defined as chainA and chainB. Both chains are reported to have functional role in ligand binding and activity. PDB file has been cleaned and prepared properly before docking using pdb-tools. Also similar cleaning and pdb file preparation was done for peptide using tools provided by haddock.
However, I keep ending up with error of “100.00% of output was not generated for this module and tolerance was set to 5.00%.”
I am new to haddock and finding it difficult to resolve the error.
Steps I have performed to solve the issue so far.
- Checked the pdb format many times.
- verified the *.cfg file and also tried adding mol_fix_origin = false
- added unambigious restraint file to define the segments of protein generated using command
haddock3-restraints restrain_bodies pdbs/fixed_apo_rep100_c0_clean.pdb > restraints/domain_restraints.tbl
- installed CNS library separately again cns_v1.3_r9, thinking this could be the source of error.
Error still persist with all the efforts.
CNS should not be the issue.
To diagnose the problem, set in the general parameter header of the config file: debug=true
The output files will then be kept and you can check those in the rigidbody directory for error messages.
Are all files present and ok in the topoaa module?
Did you check / try to run some of the provided examples?
Thank you for your prompt response.
I checked all the files and they are present in the specified folder. They are correct, to the best of my understanding based on the available tutorials.
I did run the protein-peptide tutorial examples. They are running fine with no error.
You still have most likely an issue either with your PDB files or the restraints.
Only way to help you is if you provide us all input data and your workflow.
Or set debug=true and rerun, then inspect the out files in the rigidbody directory for errors.
rigidbody_1.inp (25.8 KB)
domain_restraints.tbl (558 Bytes)
Pep_0001_ambig.tbl (2.3 KB)
I am attaching the files for your reference
run_dir = “runs/protein_Pep_0001”
mode = “local”
ncores = 20
debug = true
molecules = [
“pdbs/fixed_apo_rep100_c0_clean-fixed.pdb”,
“pdbs/Pep_0001_fix_clean.pdb”
]
disables the automatic determination of histidine protonation state is set as false
Pep_0001_fix_clean.pdb (6.4 KB)
fixed_apo_rep100_c0_clean-fixed.pdb (1.1 MB)
[topoaa]
autohis = true
[rigidbody]
ambig_fname = “restraints/Pep_0001_ambig.tbl”
Add inter-domain restraints here to keep Chains A and B together
unambig_fname = “restraints/domain_restraints.tbl”
sampling = 1000
[flexref]
tolerance = 5
ambig_fname = “restraints/Pep_0001_ambig.tbl”
Add inter-domain restraints here to keep Chains A and B together
unambig_fname = “restraints/domain_restraints.tbl”
fle_sta_2 = 1
fle_end_2 = 5
fle_seg_2 = “C”
ssdihed = “alphabeta”
mdsteps_cool1 = 2000
mdsteps_cool2 = 4000
mdsteps_cool3 = 4000
mdsteps_rigid = 2000
[emref]
ambig_fname = “restraints/Pep_0001_ambig.tbl”
Maintain restraints during final minimization
unambig_fname = “restraints/domain_restraints.tbl”
fle_sta_2 = 1
fle_end_2 = 5
fle_seg_2 = “C”
ssdihed = “alphabeta”
[seletop]
select = 200
[rmsdmatrix]
[clustrmsd]
clust_cutoff= 5
plot_matrix = true
[seletopclusts]
top_models = 4
Your receptor PDB file contains two chains.
For docking purposes an input PDB file should only contain a single chain.
So assign a single chainID to your protein, making sure there is not overlap in the numbering of residues.
Thank you for your response.
I will modify the pdb and submit the job accordingly.
However, I have question regarding haddock module. Currently it is not possible to dock multi chain protein in haddock?
Since in my case both chains are important for the binding and modulation of activity.
You simply have to treat it as one chain - as is done for example for antibodies - check our online tutorials.
And you can dock more than two molecules, but the rule remains: one molecule - one chainID