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
[email protected]
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss
