Issues when docking glycans

Hello,
I’m a webserver user of Haddock2.4 and I’m doing protein-glycan docking with glucans with mixed glycosidic linkages (beta1-3, beta1-4 and beta1-6), but, at the end of the run, beta1-6 bonds are missing. Is this a bug of the software or am I doing something wrong? here are attached:

  1. the input pdb files for the protein, for one of those glucans (6-O-β-D-glucopyranosyl-D-glucose) and for the docked structure.
  2. screenshots of the glucan before and after the run.

thank you for your support.

best

Francesco

INPUT_PROTEIN.pdb (344.7 KB)
6-O-β-D-glucopyranosyl-D-glucose.pdb (80.5 KB)
cluster1_1.pdb (433.4 KB)


Sorry to be annoying, but I need to fix this issue if possible

here is the parameter file of the job if it could be helpful.
job_params.json (440.6 KB)

Hi Francesco

It looks like that particular linkage was not implemented.
Should be soon fixed…

Ok thank you Alexandre for the feedback. Is it possible to receive a notification when the issue will be fixed?

I will let you know once properly tested and implemented on the server - probably tomorrow

Should now be fixed. Your original json file now results in a properly connected glycan:

https://bianca.science.uu.nl/haddock2.4/result/1111111111/279341-test-glycans-1-6

Yes, the glycan is properly connected. Thank you so much :pray:

I just need to ask for another help, always with glycans. I’ve tried to dock a molecule of 1,4-β-D-Galactobiose, but, at the end of the run, the beta1-4 bond is missing and a weird bond between the O3 and O4 is forming.

I’ve compared the .psf file of this glycan with the .psf file of 1,4-β-D-Glucobiose (that is a glycan with the same number of atoms and the same connectivity as the 1,4-β-D-Galactobiose) in order to see if there were differences during the topology generation and I’ve noticed that two different chemical types are assigned for the atom number 11 (O3): ‘OH1’ for 1,4-β-D-Glucobiose and ‘OAS’ for 1,4-β-D-Galactobiose, resulting in slightly different connectivity between atoms. Maybe could be this the explanation of the issue? Here are attached:

  1. the 1,4-β-D-Galactobiose input pdb file

  2. the output pdb file of the the complex

  3. the 1,4-β-D-Galactobiose.psf and the 1,4-β-D-Glucobiose.psf files (to see the difference that i was talking about)

  4. the parameter file of the job

  5. screenshots of 1,4-β-D-Galactobiose before and after the running

thank you again for the help

GA4GA_PSIMTI6.pdb (26.8 KB)
cluster1_1.pdb (433.1 KB)
1,4-β-D-Galactobiose.psf.txt (4.3 KB)
1,4-β-D-Glucobiose.psf.txt (4.3 KB)
job_params.json (395.5 KB)


Seems bad luck…

We use a distance cutoff to figure out the connectivity and in this case both O3 and O4 are within the distance cutoff (2.75A) of the C1 atom, which causes to patches to be applied.

We can shorten this cutoff a bit (might impact the detection of other linkages in case of bad geometries).
Let me run some tests first

Ok - please try again. Should be resolved now.

Yes it works! thank you :pray: