2022-01-10 利用RDKit基于分子指纹相似性对sdf小分子数据进行筛选

作者&投稿:支备 (若有异议请与网页底部的电邮联系)
~ import rdkit

from rdkit import Chem

from rdkit.Chem import Draw

from rdkit.Chem import AllChem

from rdkit.Chem import Draw

from rdkit import DataStructs

from rdkit.Chem import Descriptors

#导入相应的包和函数

suppl = Chem.SDMolSupplier('C:/Users/zhangguhua/Desktop/pycharm_job_path/Fragment_Library_4063.sdf')

#导入带筛选数据

mols = [mol for mol in suppl if mol is not None]

#读取为分子对象

fo = open('C:/Users/zhangguhua/Desktop/pycharm_job_path/Sulcotrione.smi').read()

ref_mol = Chem.MolFromSmiles(fo)

#读入参考分子 为分子对象

fo_fp = AllChem.GetMorganFingerprint(ref_mol,2)

mols_fps = [(mol,AllChem.GetMorganFingerprint(mol,2)) for mol in mols]

#计算分子指纹

fo_mols = [(m,DataStructs.TanimotoSimilarity(fo_fp,mols_fp)) for m,mols_fp in mols_fps]

sorted_fo_mols = sorted(fo_mols,reverse=True,key=lambda x: x[1])

#计算分子指纹相似性并根据大小进行排序

result = sorted_fo_mols[:200]

test_20 = open('C:/Users/zhangguhua/Desktop/pycharm_job_path/test20.sdf','w')

writer = Chem.SDWriter(test_200)

for x,y in result:

  writer.write(x)

writer.close

创建sdf文件,并将筛选得到的结果写入到sdf文件中用于后续分析


三国群英传11的作弊密码
还是楼主强,三国群英传11都玩过?偶只记得最高版本是7吧!拿2分,走人!补:到网上下个修改器,愿意怎么修改就怎么修改!

五星体育预告节目表
01:25,022极速终点线4 01:55,实况录像:2022北京冬奥会:男冰壶循环赛:意大利-中国(9-12)55'02:55,实况录像:2022北京冬奥会:男冰壶循环赛:中国-瑞士(6-5)55'03:50,实况录像:2022中国女足超级联赛第七轮:江苏无锡女足-上海农商银行55'04:50,实况录像:2022伦敦马拉松(10.2)50'06...

名侦探柯南共多少集
到现在为止TV是154集,10月20日出155集 与KID有关 很期待 最新的第十二集剧场版十月份就有中文字幕了 TV版 1996年 001 云霄飞车杀人事件 002 董事长千金诱拐事件 003 偶像密室杀人事件 004 大都会暗号地图藏宝事件 005 新干线大爆炸事件 006 情人节杀人事件 007 每月一份礼物威胁事件 008 美术馆...

022Cr22Ni5Mo3N是什么材料
022Cr22Ni5Mo3N是中国新国标双相不锈钢,对应牌号1.4462(X12CrNiMoN225)双相不锈钢,中国旧国标牌号为00Cr22Ni5Mo3N ,美国牌号为UNS S31803\/F51\/S32205\/F60\/2205(13524351833)S31803(F51)S32205(F60)双相不锈钢 标准号:ASTMA240\/A240M--01 双相不锈钢2205合金是由21%铬,2.5%钼及4.5%...

柯南全集目录
022 豪华客轮连续杀人事件(上集)023 豪华客轮连续杀人事件(下集)024 丧失记忆的美少女事件 025 真假人质绑票事件 026 爱犬约翰杀人事件 027 小五郎同学会杀人事件(上集)028 小五郎同学会杀人事件(下集)029 计算机杀人事件 030 不在场证明杀人事件 031 电视台杀人事件 032 咖啡厅杀人事件 033 ...

口袋妖怪红金手指代码大全
01e 力量粉(恢复体力50,很苦粉,减底怪兽亲密度) 01f 力量根(恢复体力200,很苦根,减底怪兽亲密度) 020 万能粉(全恢复状态,非常苦的粉,大大减低与怪兽的亲密度) 021 复活草(死亡复活,很苦的草) 022 pp艾依(1种技能值恢复10) 023 pp力卡(1种技能值全恢复) 024 ???(???) 025 pp最大(1只怪兽的全部...

名侦探柯南各集标题?
022 豪华客轮连续杀人事件(前篇) 023 豪华客轮连续杀人事件(后篇) 024 丧失记忆的美女事件 025 真假人质绑票事件 026 小五郎同学会杀人事件(前篇) 027 小五...201 第10个乘客(前篇) 202 第10个乘客(后篇) 203 伊卡路斯的黑色之翼(前篇) 204 伊卡路斯的黑色之翼(后篇) 205 本厅刑事恋爱物语3(前篇) 206 ...

柯南到现在每集的名字...现在到500集了!
022 豪华客轮连续杀人事件(上集)023 豪华客轮连续杀人事件(下集)024 丧失记忆的美少女事件 025 真假人质绑票事件 026 爱犬约翰杀人事件 027 小五郎同学会杀人事件(上集)028 小五郎同学会杀人事件(下集)029 计算机杀人事件 030 不在场证明杀人事件 031 电视台杀人事件 032 咖啡厅杀人事件 033 ...

conan 全集列表
KK019 电梯杀人事件 (エレベ?タ??人事件) 1996\/06\/10 动画原创 KK020 鬼屋杀人事件 (幽?屋敷?人事件) 1996\/06\/17 KK021 电视剧外景队杀人事件 (TVドラマロケ?人事件) 1996\/06\/24 动画原创 KK022 豪华客轮连续杀人事件(前编) (豪?客船???人事件(前?)) 1996\/07\/01 KK023 豪华客轮连续杀人事件(...

辐射4代码大全-控制台代码大全
001831B6 X-01五型躯干内衬 001831AC X-01五型头盔 001831B0 X-01五型腿甲 001D0266 X-01校准冲击 00203744 X-01液压腕甲 001D030C X-01医疗用帮浦 00203743 X-01移动辅助伺服机 00183566 X-01隐形小子 00183584 X-01侦察感应器 001EC0E8 X-01紫色头灯 001CFE46 X-01最佳化护腕 001D022F X-01最佳化...

息县13072736870: 2022年春节是几月几日(法定假日安排)
逯亨牛黄: 一、2022年春节是几月几日: 2022年春节是2022年02月01日星期二.二、这一年有大年三十吗: 2022年为什么没有大年三十 天文专家解释说,我国农历中的月是按月亮...

息县13072736870: 私家车年检新规定2022年10月
逯亨牛黄: 随着2022年新交通法规出台,关于车辆的的年检,也推出了很多新的规定.下面一起来看看2022年10月1日起私家车年检新规定有哪些变化吧.据公安部发布的新规显示,...

息县13072736870: 69亿8千万韩元等于多少人民币 -
逯亨牛黄: 6980000000韩元≈37104055.0799人民币元,由于外汇汇率是实时变动的,此数据仅供参考,具体汇率以您交易时实际交割的成交价为准,货币换算更新时间:2022-01-10 08:25.如要查询平安银行汇率,您可以登陆平安口袋银行-首页-更多-跨境金融—外汇行情,在这里您可以看到实时汇率.温馨提示: 1、以上内容仅供参考,不作任何建议. 2、入市有风险,投资需谨慎.您在做任何投资之前,应确保自己完全明白该产品的投资性质和所涉及的风险,详细了解和谨慎评估产品后,再自身判断是否参与交易. 应答时间:2022-01-11,最新业务变化请以平安银行官网公布为准.

息县13072736870: 从惠州到上海 坐什么交通工具最省钱 -
逯亨牛黄: 正解:坐火车最省钱.惠州到上海列车共3趟 T102惠州 - 上海南17:10 - 10:2617小时16分硬座 203 硬卧下 367 T212惠州 - 上海南14:33 - 07:5317小时20分硬座 203 硬卧下 367 K2056惠州 - 上海南15:00 - 14:0923小时9分硬座 203 硬卧下 367

息县13072736870: 2022虎年1月领证黄道吉日 - 结婚登记佳期
逯亨牛黄: 选择吉祥的日子对我们来说是一种美好的祝愿和期望.不仅如此,吉日会让我们感受... 那么对于日子的采用也不许大概,要择吉.2022年一月领证备案黄道好日子2022年1...

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网