I tried again to dock zinc by changing the ID. But, I do not see any improvement. The zinc is >10 nm away from the C-C-H-H binding site. I am checking if I am doing anything wrong. I have modified my systems and trying again. The details are as follows:
Step-1: Inputing the PDB structure, Chain A
(NO ERROR)
Step-2: Inputing the zinc structure obtained from the NMR structure (PDB: 1SP1). Chain D
(NO ERROR) Kind of molecule selected : LIGAND
REMARK BIOVIA PDB file
REMARK Created: 2023-06-14T09:23:10Z
***CRYST1 1.000 1.000 1.000 90.00 90.00 90.00 P 1 1 ***
HETATM 1 ZN+2 ZN2 D 30 -0.998 4.266 2.578 1.00 0.81 ZN2+
***TER 2 ZN2 D 30 ***
END
Step-3: Active/Passive residue section (Default, I didn’t define the active residues for both Molecule 1, protein, and Molecule 2, Zinc).
Step-4: You can supply a HADDOCK restraints TBL file with restraints that will always be enforced (unambiguous restraints). FILE NAME: restraints.tbl as follows
assign (segid A and resid 38 and name SG)
(segid D and resid 1 and name “ZN+2”) 3.0 0.5 0.5
assign (segid A and resid 42 and name SG)
(segid D and resid 1 and name “ZN+2”) 3.0 0.5 0.5
assign (segid A and resid 53 and name NE2)
(segid D and resid 1 and name “ZN+2”) 3.0 0.5 0.5
assign (segid A and resid 61 and name NE2)
(segid D and resid 1 and name “ZN+2”) 3.0 0.5 0.5
All other options were default.
Using these steps, the cluster I am generating has no zinc on the active site (>10 nm far).
I am trying to run the docking again by modifying step-3. This time, I am defining the active residues for molecule 1 (38, 42, 53, and 61) and molecule 2 (30).
I would appreciate it if you let me know if I am doing anything wrong and why the zinc is not parsing to the active site at all.
-Sj