接前文:
分子对接教程 | (1) 软件安装准备
分子对接教程 | (2) 选择合适的蛋白受体
分子对接教程 | (3) 配体分子文件格式转换
对于蛋白质,在对接之前我们需要预处理,比如,去水,加全氢,导出为PDBQT文件设为受体,下面我们来介绍一下怎么做。
我们这里处理蛋白文件需要借助pyMOL这个软件,我在前面提到过,我也上传到了 盘,有2个版本,一个2.2,一个2.4,前面说安装时需要安装2.0序列的Python,这是针对2.2版本的pyMOL,如果你安装2.4版本,需要安装Python3.7。我用的是2.4的版本,有点喜新厌旧啦。
首先我们打开pyMOL这个软件
或者通过File里面选择get PDB…,弹出窗口输入信息后点击下载。
然后我们可以通过序列,左右滑动,删除我们不需要的部分,比如水分子,把0全部选中,右键点击remove,就删除了。如果有其他金属离子,那就需要根据需要了,如果该离子是该蛋白的组成部分,那就不能去除。
如果我们去掉的组分比较多,比如有多条链,去掉了其中的一些链,把不需要的水分子,离子,溶剂分子去掉,只保留对对接需要的部分,对剩下的结构我们需要进行一个修复。这里需要一个 页工具,地址:https://swift.cmbi.umcn.nl/servers/html/model.html. 做法是找到 页最下面的Prepare PDB file for docking programs,点进去,上传自己的蛋白结构文件,然后点击send,稍等一下可以直接下载处理过的蛋白结构文件。
关于什么颜色代表什么原子,如下:
接下来就是加氢,因为从pdb数据库中下载蛋白质晶体结构是没有氢原子的(除了很少分辨率小于1A的蛋白质有H),这是一个技术问题。所以我们需要把氢原子加上,这一步是必须的。官方解释也是:在选择一个分子作为配体或受体之前,必须把所有的氢都加到这个分子上。注意是全氢。按照下面操作加氢。
我们就可以看见多了好多白色,白色就是氢原子。
上面这一步操作方式选择大分子会自动执行以下一系列初始化步骤,ADT检查分子是否带有电荷。如果没有,它会给每个原子加上Gasteiger电荷。记住,所有的氢必须先加到大分子上,然后才被选中。如果分子已经带电荷,ADT会问你是否想保留输入电荷而不是增加Gasteiger电荷。还有,该步骤,ADT会合并非极性氢,如果你不需要,需要手动设置。下面是官方文档的解释。
好了,保存文件后,我们可以删除分子,后续处理小分子。
文章知识点与官方知识档案匹配,可进一步学习相关知识Python入门技能树首页概览210166 人正在系统学习中
声明:本站部分文章及图片源自用户投稿,如本站任何资料有侵权请您尽早请联系jinwei@zod.com.cn进行处理,非常感谢!