Tanimoto coefficient rdkit

from rdkit import Chem
from rdkit import DataStructs
from rdkit.Chem.Fingerprints import FingerprintMols
import pandas as pd

# read and Conconate the csv's
df_1 = pd.read_csv('first.csv')
df_2 = pd.read_csv('second.csv')
df_3 = pd.concat([df_1, df_2])

# proof and make a list of SMILES
df_smiles = df_3['smiles']
c_smiles = []
for ds in df_smiles:
    try:
        cs = Chem.CanonSmiles(ds)
        c_smiles.append(cs)
    except:
        print('Invalid SMILES:', ds)
print()

# make a list of mols
ms = [Chem.MolFromSmiles(x) for x in c_smiles]

# make a list of fingerprints (fp)
fps = [FingerprintMols.FingerprintMol(x) for x in ms]

# the list for the dataframe
qu, ta, sim = [], [], []

# compare all fp pairwise without duplicates
for n in range(len(fps)-1): # -1 so the last fp will not be used
    s = DataStructs.BulkTanimotoSimilarity(fps[n], fps[n+1:]) # +1 compare with the next to the last fp
    print(c_smiles[n], c_smiles[n+1:]) # witch mol is compared with what group
    # collect the SMILES and values
    for m in range(len(s)):
        qu.append(c_smiles[n])
        ta.append(c_smiles[n+1:][m])
        sim.append(s[m])
print()

# build the dataframe and sort it
d = {'query':qu, 'target':ta, 'Similarity':sim}
df_final = pd.DataFrame(data=d)
df_final = df_final.sort_values('Similarity', ascending=False)
print(df_final)

# save as csv
df_final.to_csv('third.csv', index=False, sep=',')

Source: stackoverflow.com

Similar codes

0 vote326 views1 answer
Percentage of null values for every variable in dataframe
Python
0 vote237 views2 answers
Python fraction simplify
Python
5 votes348 views2 answers
How to set icon in tkinter
Python
0 vote286 views1 answer
Call a function with prameters inm tkinter buttion
Python
10 votes325 views3 answers
Python tostring
Python
0 vote297 views1 answer
C# script for download music from telegram channel
Python
5 votes311 views2 answers
Read matlab file in python
Python
-1 vote287 views1 answer
Save a file as a pickle
Python
8 votes302 views3 answers
Len python
Python
0 vote292 views1 answer
Lists and for loops
Python