一、代码来源：

计算rmsd参考：GitHub - 0ut0fcontrol/isoRMSD: calculating RMSD between 2 conformers with different atom names

这个github代码是在命令行中输入命令运行的，命令如下：

python isoRMSD.py -r crystal_ligand.sdf -p crystal_out_1.sdf -o rmsd.txt

这样执行计算一对小分子还算方便，小分子格式支持mol2，sdf，mol，pdb，虽然支持这么多种格式，但是计算的时候有些能用mol2算却不能用sdf算，如果计算不出来可以转换一下小分子的格式没准就可以了。但是这个代码批量计算小分子rmsd时实现就很麻烦，所以把主要用的的代码拿出来，需要批量的部分，保存结果的部分可以根据需求自己写写。-r后面接的是参考分子。

二、主要代码：

根据需求写读取多个文件的循环，结果保存也可以根据需求自己写

import os
import sys
import math
from numpy.lib.shape_base import column_stack
import pandas as pd
from os.path import split
from oddt.toolkits import rdk as toolkit
from rdkit import Chem
from rdkit.Chem import AllChem
from rdkit.Chem.AllChem import AlignMolrfile = 'crystal_ligand.sdf'#参考小分子
pfile = 'crystal_out_0.sdf'
ref_fmt = rfile.split(".")[-1]
ref_oddt = next(toolkit.readfile(ref_fmt, rfile))
ref_rdk = Chem.RemoveHs(ref_oddt.Mol)probe_fmt = pfile.split(".")[-1]
probe_oddt_supp = toolkit.readfile(probe_fmt, pfile)data = []
for i, probe_oddt in enumerate(probe_oddt_supp):if probe_oddt is None:name = "NaN"rmsd_notalign = 10000.0rmsd_align = 10000.0else:probe_rdk = Chem.RemoveHs(probe_oddt.Mol)try:name = probe_rdk.GetProp("_Name")name = "_".join(name.split())except KeyError as e:name = "NaN"ref = AllChem.AssignBondOrdersFromTemplate(probe_rdk, ref_rdk)rmsd_notalign = GetBestRMSD(probe_rdk, ref, align=False)data.append((name,rmsd_notalign))print(data)
#下面是定义计算的函数，可以不动，上面加上批处理的一些代码就好
def GetBestRMSD(probe, ref, refConfId=-1, probeConfId=-1, maps=None, align=True):# When mapping the coordinate of probe will changed!!!ref.pos = orginXYZ(ref)probe.pos = orginXYZ(probe)try:name = probe.GetProp("_Name")except KeyError as e:name = "NaN"if not maps:matches = ref.GetSubstructMatches(probe, uniquify=False)if not matches:raise ValueError("mol %s does not match mol %s"% (ref.GetProp("_Name"), probe.GetProp("_Name")))if len(matches) > 1e6:print("{} matches detected for molecule {}, this may lead to a performance slowdown.".format(len(matches), name))maps = [list(enumerate(match)) for match in matches]bestRMSD = 10000.0for amap in maps:if align:rmsd = AlignMol(probe, ref, probeConfId, refConfId, atomMap=amap)else:rmsd = RMSD_NotAlign(probe, ref, amap)bestRMSD = min(bestRMSD, rmsd)return bestRMSDdef RMSD_NotAlign(probe, ref, amap):rmsd = 0.0atomNum = ref.GetNumAtoms() + 0.0for (pi, ri) in amap:posp = probe.pos[pi]posf = ref.pos[ri]rmsd += dist_2(posp, posf)rmsd = math.sqrt(rmsd / atomNum)return rmsddef orginXYZ(mol):mol_pos = {}for i in range(0, mol.GetNumAtoms()):pos = mol.GetConformer().GetAtomPosition(i)mol_pos[i] = posreturn mol_posdef dist_2(atoma_xyz, atomb_xyz):dis2 = 0.0for i, j in zip(atoma_xyz, atomb_xyz):dis2 += (i - j) ** 2return dis2

阅读代码细节的一些参考链接：

1）Python中if __name__ == ‘__main__‘：的作用和原理

(20条消息) Python中if __name__ == ‘__main__‘：的作用和原理_农村詹姆斯的博客-CSDN博客

2）split(“.“)[-1]：(20条消息) split(“.“)[-1]的怎么用_花落雨微扬的博客-CSDN博客_split[-1]

3）parser = argparse.ArgumentParser()：

(20条消息) 常用编程记录——parser = argparse.ArgumentParser()_umbrellalalalala的博客-CSDN博客

4）python 判断变量是否是 None 的三种写法：

(20条消息) python 判断变量是否是 None 的三种写法_Python热爱者的博客-CSDN博客

5）python中{}，[]，()，字典，列表，元组详解： Python基础括号()[]{}的详解_python_脚本之家 (jb51.net)

6）GetSubstructMatches（）：rdkit.Chem.rdchem module — The RDKit 2022.09.1 documentation

RDKit|分子子结构的删除、替换与切割 - 简书 (jianshu.io)

7）AllChem.AssignBondOrdersFromTemplate():rdkit.Chem.AllChem module — The RDKit 2022.09.1 documentation

8）关于rdkit包中的函数可以用下面这个网页索引：Index — The RDKit 2022.09.1 documentation

计算小分子之间的RMSD相关推荐

Python：计算两个蛋白或小分子之间的RMSD
Python脚本:计算两个蛋白或小分子之间的RMSD 用法: python rmsd.py protein1.pdb protein2.pdb rmsd.py # Root-mean-square d ...
计算两个蛋白或小分子之间的RMSD
计算两个蛋白或小分子之间的RMSD # Root-mean-square deviation (RMSD) for proteins and/or ligands # in PDB files. #c ...
RDKit：计算不同小分子构象之间的RMSD
计算两个小分子之间的RMSD,可以用来判断两个构象的接近程度. RDKit:计算不同分子或构象之间的RMSD点击打开链接第一步:安装Anaconda Win或者Linux系统下Anaconda安装, ...
RDKit：计算不同分子或构象之间的RMSD
RDKit:计算不同分子或构象之间的RMSD Linux(CentOS 7_x64位)系统下安装RDkit(修正)点击打开链接有时我们需要计算两个(或更多)分子之间的RMSD,可以用来判断两个构象的 ...
薛定谔 | 小分子叠合
所谓分子叠合就是将同一种分子不同构象移动到一起,使得二者整体保证原子的叠合.这类操作之后往往会有一个RMSD的指标用于衡量这两个分子之间的空间相似度.数值越小,二者差异越小.蛋白.核酸的分子叠合很多软 ...
分子模拟软件amber_使用Amber创建小分子与蛋白质复合蛋白的坐标和拓扑文件
复合蛋白amber坐标和拓扑文件的创建作者:朱宁来源:大科研小分享前言分子动力学(Molecular Dynamics, MD)是一门结合物理,数学和化学的综合技术.目前主流分子动力学软 ...
微信小程序高德地图知道两点求道路_微信小程序——计算2点之间的距离
关于计算2点之间的距离都依赖了腾讯地图,所以请先在下面具体讲计算2点之间距离的方法. 方法一: getPosition: function() {var that = this; wx.getLoc ...
小分子药物logp预测/计算
众所周知,有机小分子拥有很多内禀属性,如logP, logD,溶解度,和pKa等.可以帮助我们在实验过程中快速评估和预测小分子行为.如化合物的亲脂性常常由logP和logD表示.所谓logP是对P(分 ...
微信小程序——计算2点之间的距离
关于计算2点之间的距离都依赖了腾讯地图,所以请先在腾讯地图官网申请key.具体流程看下图: 下面具体讲计算2点之间距离的方法. 方法一: 1.通过 wx.getLocation(Object obje ...

计算小分子之间的RMSD

一、代码来源：

二、主要代码：

计算小分子之间的RMSD相关推荐

最新文章

热门文章