the minimizedAffinity value of output.sdf results in one of:
-10.98858
-10.48566
-9.52101
Using --seed doesn't give deterministic results either.
I attached another pair of examples "8ACL_LQL" that shows the same problem.
As I tested, most of the inputs give deterministic results. The above problem occurs about ~4% of cases.
Please let me know if I can help finding out the cause.
Smina version: Smina Oct 15 2019. Based on AutoDock Vina 1.1.2.
OS: Red Hat Enterprise Linux release 9.3 (Plow)
Attached files are from: https://github.com/maabuu/posebusters
The ligand is charged which means OpenBabel wants to add a hydrogen to the charged nitrogen, but it doesn't really know how and uses a random vector resulting in non-deterministic output. If the charge is removed or if the ligand is first converted to a different format (e.g. pdbqt) you get deterministic output. You also get deterministic output with gnina, since gnina explicitly initializes the random number generator that OpenBabel uses with the user provided seed. I've backported this behavior to smina here: https://sourceforge.net/p/smina/code/ci/02b978cc88e71be64ed3e40036969b435714c665/
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
There are some cases when repeated run with options
--local_only
and/or--minimize
gives varying results.For example, when I repeatedly run the following command with the files I attached:
the
minimizedAffinity
value ofoutput.sdf
results in one of:-10.98858
-10.48566
-9.52101
Using
--seed
doesn't give deterministic results either.I attached another pair of examples "8ACL_LQL" that shows the same problem.
As I tested, most of the inputs give deterministic results. The above problem occurs about ~4% of cases.
Please let me know if I can help finding out the cause.
Smina version: Smina Oct 15 2019. Based on AutoDock Vina 1.1.2.
OS: Red Hat Enterprise Linux release 9.3 (Plow)
Attached files are from: https://github.com/maabuu/posebusters
Last edit: Sang-Yeon Hwang 2024-07-09
The ligand is charged which means OpenBabel wants to add a hydrogen to the charged nitrogen, but it doesn't really know how and uses a random vector resulting in non-deterministic output. If the charge is removed or if the ligand is first converted to a different format (e.g. pdbqt) you get deterministic output. You also get deterministic output with gnina, since gnina explicitly initializes the random number generator that OpenBabel uses with the user provided seed. I've backported this behavior to smina here:
https://sourceforge.net/p/smina/code/ci/02b978cc88e71be64ed3e40036969b435714c665/