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

Set demerit threshold #19

Open
dence opened this issue Jun 21, 2022 · 1 comment
Open

Set demerit threshold #19

dence opened this issue Jun 21, 2022 · 1 comment

Comments

@dence
Copy link

dence commented Jun 21, 2022

Hi, I'd like to get the demerit score for all the input smiles strings. Basically, I'd like to remove the cutoff threshold for output.

In the README.md file, it says, "For details adjusting the settings, see the program's documentation". I can't find any commandline option in documentation.md or in the Lilly_Medchem_Rules.rb file where the threshold can be adjusted.

Is that something that I can set on the command-line?

@nbehrnd
Copy link
Contributor

nbehrnd commented Jun 22, 2022

@dence From limited code archeology, my speculation is that the rules around the critical threshold reside in the C/C++ files. In file /Molecule/iwdemerit.cc for example contains dozens of functions around demerit and the structure related demerits ("butyl, pentyl, hexyl, heptyl, and cyclohexyl group in a molecule equate 10, 25, 50, 100, and 170" README.md
mentions likely compose from file Molecule/substructure_demerits.cc.

Indeed, running the program in a pattern like

ruby Lilly_Medchem_Rules.rb -v input.smi > okmedchem.smi

to enable the verbose mode does not automatically yield a report about all molecules submitted .and. their assigned demerits either (even if, with the set of rules applied, these would be considered as good enough to enter okmedchem.smi).

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

2 participants