HI Xiaobo, There's an explicit hydrogen in the SMARTS that shouldn't be there. I also wouldn't include the single bonds around the ring closures.
'[#8]=[#6]-3-c1c2c(ccc1)cccc2-[#6](-[#7]-3-*[#1]*)=[#8]') from rdkit import Chem from rdkit.Chem import Draw smi = "O=C(C1=C2C(C=CC=C23)=CC=C1)N([H])C3=O" mol = Chem.MolFromSmiles(smi) mol_list = [mol] core = Chem.MolFromSmarts("[#8]=[#6]3-c1c2c(ccc1)cccc2-[#6](-[#7H]3)=[#8]") Draw.MolsToGridImage(mol_list,highlightAtomLists=[x.GetSubstructMatch(core) for x in mol_list]) [image: image.png]
_______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss