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, 1016🔥, 0💬
Popular Posts:
How to display a 3-D molecule structure with 3Dmol.js? If you have a 3-D molecule structure stored i...
Molecule Summary: ID: FYI-1002995 Names: InChIKey: XDFDFUJTCOCLHM-UHFFFAOYS A-NSMILES: Cc1ccccc1c5cc...
Molecule Summary: ID: FYI-1002537 Names: InChIKey: WSLBJQQQZZTFBA-QQSPYBPKS A-NSMILES: CCCCC/C=C/C/C...
Molecule Summary: ID: FYI-1002057 Names: InChIKey: CQERVFFAOOUFEQ-UHFFFAOYS A-OSMILES: O=C(c1cncc(Br...
Molecule Summary: ID: FYI-1005761 Names: InChIKey: JACRWUWPXAESPB-UHFFFAOYS A-NSMILES: C1=CC=C(C=C1)...