python - 使用 RDKit 计算 sdf 文件和结构 SMILE 之间的 Tanimoto 相似性？

Question

我正在使用带有 Python 3.7 的 RDKIt 来计算 sdf 中的数据库（每个结构的微笑）与一个分子的相似性，我对此微笑。我找到了一种仅使用以下代码在两个 SMILES 之间计算 Tanimoto 指数的方法：

import numpy as np 
import scipy
import matplotlib
import matplotlib.pyplot as plt
import rdkit as  rd
from rdkit import Chem

ref = Chem.MolFromSmiles('Nc1nc2nc(N)nc(N)c2nc1-c1cccc(Cl)c1')
mol1 = Chem.MolFromSmiles('structure smiles')
fp1 = Chem.RDKFingerprint(ref)
fp2 = Chem.RDKFingerprint(mol1)

Tan  =DataStructs.TanimotoSimilarity(fp1,fp2)

print (Tan)

有没有办法用 sdf 文件替换 mol1 ？

Answer 1

您可以使用SDMolSupplier.

from rdkit import Chem, DataStructs

ref = Chem.MolFromSmiles('Nc1nc2nc(N)nc(N)c2nc1-c1cccc(Cl)c1')
fp1 = Chem.RDKFingerprint(ref)

suppl = Chem.SDMolSupplier('yourSDF.sdf')
for mol in suppl:
    fp2 = Chem.RDKFingerprint(mol)
    Tan = DataStructs.TanimotoSimilarity(fp1,fp2)
    print(Tan)