RDKit：运用RDKit计算USRCAT（形状相似性）

2024-04-22 12:35:16

USRCAT

USRCAT是基于形状的方法，它的工作速度非常快。代码是免费提供的，如果要使用代码，用户需要安装它。

形状相似性是分子的3D形状的比对，其基于如果分子的形状重叠良好如果可以相似的原理。基于形状筛选有效地用于先导化合物优化，并通过叠加相似化合物了解构效关系。此外，许多研究表明其在虚拟筛选实验中的成功应用。基于形状相似性得分作为选择的评价标准，并对命中化合物进行排序。

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3505738/

RDKit代码：

import os

import seaborn as sns

import pandas as pd

from rdkit import Chem

from rdkit.Chem import rdBase

from rdkit.Chem import RDConfig

from rdkit.Chem import AllChem

from rdkit.Chem.rdMolDescriptors import GetUSRScore, GetUSRCAT

from rdkit.Chem import DataStructs

print( rdBase.rdkitVersion )

mols = [ mol for mol in Chem.SDMolSupplier( "cdk2.sdf" ) ]

for mol in mols:

AllChem.EmbedMolecule( mol,

useExpTorsionAnglePrefs = True,

useBasicKnowledge = True )

usrcats = [ GetUSRCAT( mol ) for mol in mols ]

fps = [ AllChem.GetMorganFingerprintAsBitVect( mol, 2 ) for mol in mols ]

data = { "tanimoto":[], "usrscore":[] }

for i in range( len( usrcats )):

for j in range( i ):

tc = DataStructs.TanimotoSimilarity( fps[ i ], fps[ j ] )

score = GetUSRScore( usrcats[ i ], usrcats[ j ] )

data["tanimoto"].append( tc )

data["usrscore"].append( score )

print( score, tc )

df = pd.DataFrame( data )

fig = sns.pairplot( df )

fig.savefig( 'plot.png' )

（阅读原文查看完整代码，公众号排版不好）

jupyter notebook效果

RDKit：运用RDKit计算USRCAT（形状相似性）相关推荐

RDKit | RDKit（2019.09）新增相似性图函数
RDKit(2019.09)新增相似性图函数导入库 from rdkit import Chem from rdkit.Chem import Draw from rdkit.Chem.Draw i ...
化合物相似性搜索_分子形状相似性及其应用
1.分子形状(molecular shape) 分子形状是药物化学中的基本概念,基于形状的虚拟筛选已经成功识别出新的抑制剂.形状用于评估候选分子与一组已知活性物质的相似性,并评估分子形状与靶受体上结合 ...
提出用于“L”形建筑物识别的形状相似性计算方法。
为了识别"L"形建筑物,可以使用形状相似性计算方法.具体步骤如下: 将建筑物的轮廓提取出来,可以使用边缘检测等方法. 对于每个边缘点,计算其到建筑物轮廓的距离和方向.可以使用极坐标 ...
RDKit | 基于RDKit(≥2020.09.1)的相似图绘制新方法
导入库 from rdkit import Chem from rdkit.Chem import Draw from rdkit.Chem.Draw import SimilarityMaps fr ...
RDKit | 基于RDKit绘制带原子和键的索引、注释和立体化学
基于RDKit绘制带原子和键的索引.注释和立体化学导入库 from rdkit import Chem from rdkit.Chem import Draw from rdkit.Chem.Dra ...
RDKit | 基于RDKit可视化脂溶性在溶解度贡献中的应用
尝试了logP,仅使用原子的权重即可可视化. 导入库 import base64 from io import BytesIO import warnings warnings.filterwarni ...
RDKit | 基于RDKit和Cytoscape绘制分子相似图
化学信息学中有许多网络结构化数据.例如分子,分子相似图和MMP等. Cytoscape是一款图形化显示网络并进行分析和编辑的软件基于RDKit和Cytoscape绘制分子的相似图. py2cytos ...
ITK:计算PCA形状模型
ITK:计算PCA形状模型内容提要 C++实现代码内容提要根据训练样本计算PCA形状模型. C++实现代码 #include "itkBinaryThresholdImageFilte ...
opencv学习——Moments（）函数，计算物体形状方向
opencv中的矩主要包括以下几种:空间矩,中心矩和中心归一化矩. class Moments { public: - // 空间矩 double m00, m10, m01, m20, m11, m ...

最新文章

热门文章