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 vote298 views2 answers
Python multilevel list comprehension
Python
3 votes333 views2 answers
Getting text input pygame python
Python
1 vote349 views1 answer
Counter in django template
Python Django
0 vote308 views1 answer
Customizable TKinter Buttons Python
Python
1 vote326 views2 answers
Is : and :: the same in python slice
Python
0 vote304 views1 answer
Python skip line
Python
1 vote303 views2 answers
Python pyqt5 table widget remove horizontal header
Python
0 vote294 views2 answers
Python to c converter tool
Python
0 vote391 views1 answer
How to fix invalid salt in python flask
Python Flask
13 votes281 views1 answer
Install python2.7 ubuntu
Python