RDKit:基于RECAP生成片段

#导入模块,载入分子
#!usr/bin/python3
from rdkit import Chem
from rdkit.Chem import Recap
from rdkit.Chem import Draw
from rdkit.Chem.Draw import IPythonConsole
mol = Chem.MolFromSmiles('O=c1nc([nH]c2NCC(Nc12)CN(c1ccc(cc1)C(=O)NC(CCC(=O)O)C(=O)O)C=O)N')
#输出分子
mol
RDKit:基于RECAP生成片段_第1张图片
#结果以树结构形式获得
node = Recap.RecapDecompose(mol)
print (node)
#输出叶节点片段
leaves = [leaf.mol for leaf in node.GetLeaves().values()]
img = Draw.MolsToGridImage(leaves)
img.save('leaves.png')
RDKit:基于RECAP生成片段_第2张图片
#输出所有片段
all_nodes = [node.mol for node in node.GetAllChildren().values()]
img = Draw.MolsToGridImage(all_nodes)
img.save('all_nodes.png')
RDKit:基于RECAP生成片段_第3张图片

RDKit:基于RECAP生成片段_第4张图片
DrugAI

https://blog.csdn.net/u012325865/article/details/81570965

你可能感兴趣的:(RDKit:基于RECAP生成片段)