Good afternoon folks, I would imagine that if you remove the hydrogens, the resulting molecule would be a substructure of the original molecule. However, when I do the following to the attached MDL mol file, there is no matches.
>>> from rdkit import Chem
>>> mol = Chem.MolFromMolFile('temp.mol', removeHs=False)
>>> mhvy = Chem.RemoveHs(mol)
>>> matches = mol.GetSubstructMatches(mhvy)
matches turns out to be empty.
Is it something to do with the difference between Smarts
and Smiles? If so, how can I work around this to obtain
the atomic index relationship between the two molecules?
Thank you for your insight.
Ling
temp.mol
Description: Binary data
------------------------------------------------------------------------------ Introducing Performance Central, a new site from SourceForge and AppDynamics. Performance Central is your source for news, insights, analysis and resources for efficient Application Performance Management. Visit us today! http://pubads.g.doubleclick.net/gampad/clk?id=48897511&iu=/4140/ostg.clktrk
_______________________________________________ Rdkit-discuss mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

