Hello!

I am trying to calculate the shape Tanimoto distance between two molecules (m1 
and m2) but I am finding that I get different values depending on the order in 
which I input m1 and m2 into the ShapeTanimotoDist function. Should they not be 
the same? Here is my code and I attach the two ligand sdfs:


In [1]: from rdkit import Chem

In [2]: m1 = Chem.MolFromMolFile("5qgn_aligned_lig1.sdf")

In [3]: m2 = Chem.MolFromMolFile("5qgi_aligned_lig.sdf")

In [4]: from rdkit.Chem import rdShapeHelpers

In [5]: rdShapeHelpers.ShapeTanimotoDist(m1, m2)
Out[5]: 0.8132543103448275

In [6]: rdShapeHelpers.ShapeTanimotoDist(m2, m1)
Out[6]: 0.8145196036191297

Thanks,

Susan

Attachment: 5qgn_aligned_lig1.sdf
Description: 5qgn_aligned_lig1.sdf

Attachment: 5qgi_aligned_lig.sdf
Description: 5qgi_aligned_lig.sdf

_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to