Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Protein structure used does not change when using autodock_vina_descriptor.build() #162

Open
OliverT1 opened this issue Aug 24, 2022 · 0 comments

Comments

@OliverT1
Copy link

OliverT1 commented Aug 24, 2022

When scoring multiple protein / ligand files using autodock_vina_descriptor the protein structure used is not changed, but the ligand file is changed. This means that all ligands are scored against the first protein used for autodock_vina_descriptor.build().

e.g.

proteins = ['protein1.pdb', 'protein2.pdb', 'protein3.pdb']
ligands = ['ligand1.sdf', 'ligand2.sdf', 'ligand3.sdf']
AV = autodock_vina_descriptor()
out = np.zeros((3, 6))
for i in range(len(proteins)):
    ligand_sdf = next(oddt.toolkit.readfile('sdf', ligands[i]))
    protein = next(oddt.toolkit.readfile('pdb', proteins[i]))
    protein.protein = True
    AVDs_all[i,:] = AV.build(protein = protein, ligands = ligand_sdf)

Expected behaviour: ligand1 scored against protein1, ligand2 scored against protein2, ligand3 scored against protein3
Actual behaviour: ligand1 scored against protein1, ligand2 scored against protein1, ligand3 scored against protein1
Temporary fix: Reinstantiate autodocok_vina_descriptor each time.

The analogous oddt_vina_descriptor behaves as expected.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant