Dear all, The following code snippet compares two resonance structures of formate anion:
import rdkit from rdkit import Chem mol1=Chem.MolFromSmiles('C([O-])=O') mol2=Chem.MolFromSmiles('C(=O)[O-]') mol1.GetSubstructMatches(mol2, uniquify = False) ((0, 2, 1),) mol1.GetSubstructMatches(mol1, uniquify = False) ((0, 1, 2),) I would rather like to get, in both cases, the following output: ((0, 1, 2),(0, 2, 1)) which would account for the carboxylate group symmetry due to resonance. The same applies to amidinium, guanidinium, etc. Is that currently feasible within the RDKit API? Thanks in advance, cheers Paolo ------------------------------------------------------------------------------ _______________________________________________ Rdkit-discuss mailing list Rdkit-discuss@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/rdkit-discuss