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

-1 vote147 views1 answer
Code error correction
Python
1 vote136 views1 answer
Python connectionreseterror: [winerror 10054] an existing connection was forcibly closed by the remote host
Python
14 votes97 views2 answers
Tkinter clear canvas
Python
20 votes148 views5 answers
Join lists python
Python
0 vote105 views1 answer
FunÃ§Ã£o anonima python
Python
0 vote145 views1 answer
Queryset to list python
Python
2 votes140 views2 answers
Python get time difference in milliseconds
Python
4 votes163 views1 answer
Dataframe to csv without ids
Python
0 vote215 views1 answer
Typecasting python
Python
0 vote133 views1 answer
Sending email in django
Python