Collections:
"babel -i fs ... -s query_file" - Substructure Search
Why am I getting no results with "babel -i fs ... -s query_file" command? It should perform a substructure search in a fastsearch index file.
✍: FYIcenter.com
According to the "babel -H fs" help document, you can perform a substructure search in a fastsearch index file with the query molecule stored in a file in any format.
However, Open Babel 2.3.1 seems to have an issue to support this. Here is an example session:
fyicenter$ babel -V Open Babel 2.3.1 -- Oct 13 2011 -- 15:14:47 fyicenter$ babel -i fs molecules-index.fs -o smiles result.smi -s 'c1ccccc1' 210 candidates from fingerprint search phase 208 molecules converted 7241 audit log messages fyicenter$ ls -l result.smi -rw-r--r-- 1 fyicenter staff 13889 May 25 14:18 result.smi fyicenter$ babel -i fs molecules-index.fs -o smiles result.smi -s benzene.sdf 210 candidates from fingerprint search phase 0 molecules converted 3836 audit log messages fyicenter$ ls -l result.smi -rw-r--r-- 1 fyicenter staff 0 May 25 14:17 result.smi
As you can see, if the query molecule, benzene ring, is specified in a SMILES string in the command line, "babel" found 210 candidates in the first step, and confirmed/converted 208 molecules as the substructure match result.
However, if the query molecule, benzene ring, is specified in in a SDF file, "babel" found 210 candidates in the first step, and failed to confirm/convert any molecules as the substructure match result.
This is a bug in Open Babel 2.3.1. There is no problem, if you are using Open Babel 2.4.1.
⇒ "babel -i fs ... -s query_file" - Exact Match Search
⇐ "babel -i fs ... -s SMILES" - Substructure Search with Index
2020-05-25, 739🔥, 0💬
Popular Posts:
Molecule Summary: ID: FYI-1001101 SMILES: CCOP(=O)(C)SCCN(C(C)C)C( C)CReceived at FYIcenter.com on: ...
Molecule Summary: ID: FYI-1003548 Names: InChIKey: RUVINXPYWBROJD-ONEGZZNKS A-NSMILES: C/C=C/c1ccc(O...
Molecule Summary: ID: FYI-1003902 Names: InChIKey: CDXLMBAUDBJELC-AATRIKPKS A-NSMILES: CCc4nc(C(N)=O...
Molecule Summary: ID: FYI-1000176 SMILES: COCC(=O)Cc1ccc(CC#N)cc1 Received at FYIcenter.com on: 2020...
Molecule Summary: ID: FYI-1003057 Names: InChIKey: JUQSUBXNHVMDLB-WCQYABFAS A-NSMILES: N#Cc3cc(c2cc(...