RDKit | 子结构搜索和MCS算法
【摘要】 MCS:Maximum Common Substructure
即使查询分子与数据库分子不完全匹配,也想知道其中有多少个通用结构,即最大公共子结构(MCS)。可以使用与子图同构相同的方法来计算。也可以按原样使用公共键(边)的数量作为阈值,或将其转换为相似性指标,例如Jaccard / Tanimoto系数。
&nbs...
MCS:Maximum Common Substructure
即使查询分子与数据库分子不完全匹配,也想知道其中有多少个通用结构,即最大公共子结构(MCS)。可以使用与子图同构相同的方法来计算。也可以按原样使用公共键(边)的数量作为阈值,或将其转换为相似性指标,例如Jaccard / Tanimoto系数。
部分结构匹配的情况下,可以在结构匹配时(或确定它们不匹配时)中止搜索,但是在MCS的情况下,可以输出最优解,直到搜索到所有可能性为止。有必要设计诸如确定计算时间的上限,当公共边缘的数量超过阈值时中止搜索或者使用高速近似解算法的手段。
尽管即使使用VF2算法也可以计算MCS,但已经开发了许多更高效且针对特定应用的算法。作为开源的MCS实现,RDKit具有基于DFS的算法FMCS和VF2。
导入库
-
from rdkit import rdBase, Chem
-
from rdkit.Chem import AllChem, Draw
-
from rdkit.Chem import rdFMCS
-
from rdkit.Chem.Scaffolds import MurckoScaffold
-
import pubchempy as pcp
载入数据
文章来源: drugai.blog.csdn.net,作者:DrugAI,版权归原作者所有,如需转载,请联系作者。
原文链接:drugai.blog.csdn.net/article/details/102626236
【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱:
cloudbbs@huaweicloud.com
- 点赞
- 收藏
- 关注作者
评论(0)