分子对接教程 | (8) PyMOL可视化对接结果

TCGA | GEO | 文献阅读 | 数据库 理论知识

R语言 | Bioconductor | 服务器与Linux



接前文:

分子对接教程 | (1) 软件安装准备

分子对接教程 | (2) 选择合适的蛋白受体

分子对接教程 | (3) 配体分子文件格式转换

分子对接教程 | (4) 蛋白受体文件的预处理

分子对接教程 | (5) 配体小分子的预处理

分子对接教程 | (6) AutoDock对接操作与对接结果解读

分子对接教程 | (7) AutoDock对接中易错问题


在介绍之前,先简单介绍一下这个软件,虽然前面我们简单的使用,但没有过多介绍,这里就简单介绍一下,具体细节上的东西,需要你不断使用,才能熟悉。下面是软件界面。

1.设定工作目录,打开文件

打开PyMOL软件,设置我们的工作目录。

或者通过命令cd F:\AutoDock来实现,这与window的dos命令行和Linux系统的cd(Change Directory)命令一样。

导入刚才下载的复合物PDB文件(1e8y.pdb)。

注意:有两种方法导入文件,一种是鼠标点击菜单栏File -> Open…,另一种是通过命令行。通过load命令加载:load 1e8y.pdb,因为这样可以同时把工作目录设置好,方便后续保存文件。

假设我们的工作目录为F:\AutoDock。注意:路径中不要有中文和空格,这是导致很多计算类软件(尤其是免费软件)出现问题的原因。

如果是从PDB网站上下载蛋白,如 1e8y.pdb, 也可以直接通过PyMOL 下载蛋白,点击菜单栏中的 File->Get PDB,或者通过fetch 命令下载,这在前面有介绍。

2.对象属性操作介绍

加载到蛋白后。

  • 1个是all, all 不是真实的object,它代表了所有的object

  • 1个是1e8y ,1e8y就是我们刚刚载入的蛋白

每个object 都有对应的A S H L C操作:

A:代表Action ,主要包含了对object的常用操作的集合,如 复制、删除object,对object加氢,展示Object等,

S:代表Show ,将object 渲染成cartoon 、line、stick lines sphere surface mesh dots ribbon 等模式

H:Hide 根据object的状态或者描述进行相应的掩藏

L:Label 显示object中残基、原子等名称或者属性

C:Color 对Object 进行着色

知道这些,里面的内容多试试几次就知道怎么回事啦。比如下面2种操作。

  • 点击 A->preset->simple 显示蛋白的简单形式

  • 点击 A->preset->ball and stick 显示球棍模型

再如下面的一些操作介绍;

L->residue 在α碳原子上标记其残基名字和编号 (常用)

L->residue name 在所有原子上标记残基名字 (不常用)

L->clear 删除该对象上所有的Label

L->element symbol 显示对象上所有原子的元素名字

L->vdw radius 看原子的范德华半径

3.可视化窗口操作

  • 平移,按住鼠标中键不放,然后上下左右移动,进行体会,蛋白会随着鼠标而移动

  • 旋转,按住鼠标左键不放,然后上下左右移动鼠标,蛋白会进行旋转

  • 缩放,按住鼠标右键不放,然后上下移动,蛋白会进行缩放

  • 切割 滚动鼠标中键, 建议将蛋白渲染成surface模式,然后滚动鼠标中键

多的不介绍了,B站有很多视频教程,这里我给大家找到了一个文本教程。http://pymol.chenzhaoqiang.com/ ,对作者的贡献表示感谢。

下面我们就简单处理一下我们前面对接的结果

4.对接结果简单处理演示

我们前面对接的结果文件,result.pdbqt,我们同样用OpenBabel这个软件转换成pdb格式。

我们用pymol打开result.pdb文件,在右下角处Selecting 处点击切换到Chains,然后选中蛋白质,在左上角sele处选择A,然后选择rename selection,然后在左上角的窗口处会弹出一行,在:后面输入名称。我这里是protein。或者通过命令set_name sele,protein实现。

我下面通过命令set_name protein,p又将protein改为p。

为了更好的区分受体和配体,我们通过C(color)来更改,根据自己喜好。

Selecting 处切换到残基。

然后选择小分子,右侧又多了一个sele,我们同样可以更改名字。

我这里把名称改成了ligand,我们选中小分子,按下图选择,让小分子显示氢键。

我们鼠标旋转,可以看见4个氢键。

接下来我们显示这4个氢键对接在氨基酸上的那几个残基上。首先,点击蛋白质(p)的S,点击sticks。

然后点击ligand的A,点击center,将配体小分子设置为中心。后面可以通过鼠标放大缩小旋转都以它为中心。

然后我们放大(长按鼠标右键拖动),旋转到合适角度,可以看见。

然后鼠标选中对接的残基

同样将残基改一个名字,我这里是rr。

在p中的H点击sticks,将蛋白的棍状结构隐藏起来。

然后是在rr中选择S里面的sticks,单独显示残基的棍状结构。

我们可以给残基换一个颜色(C)。

改变背景色,这里改为白色,同时我再把残基改成了蓝色。

可以保存文件编辑好的图片。

至于美化,这里远远不够,就得去自己好好研究这个软件怎么使用了,反正上面的教程也挺详细的,后续在介绍另外一个软件VMD。

参考:

http://pymol.chenzhaoqiang.com/ 

  • 15
    点赞
  • 132
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
Autodock是一种常用的分子对接软件,用于模拟分析分子间的结合和相互作用情况。分析Autodock分子对接结果一般可包括以下几个方面的内容: 1. 结合自由能:分析分子对接结果中的结合自由能,它是一个衡量分子结合稳定性的指标。较低的结合自由能值通常表示分子间的结合更稳定。可以观察不同配体(小分子)在相同受体上的结合自由能值,以判断它们的结合能力。 2. 互作模式:根据分子对接得到的复合物结构,分析受体和配体之间的相互作用。通过观察氢键、离子对、疏水作用等相互作用方式,可了解受体和配体之间的结合模式和作用力。 3. 残基相互作用:分析分子对接过程中重要氨基酸残基(如受体中的关键氨基酸)与配体之间的相互作用情况。观察这些残基在结合过程中的关键作用可以帮助了解受体位置、构象和结合特点。 4. 槽口分析:对于受体家族,可能存在特定的配体结合槽口。通过分析分子对接结果,可以评估配体是否能够正确地进入槽口,并与受体结合。 5. 结构优化:在得到对接结果后,可以对复合物结构进行进一步优化。根据评估指标选择较为合理的构象,进行分子力学模拟或量子力学计算,并进行结构优化,以进一步验证分子对接结果的可信性。 通过以上分析,可以更全面地理解Autodock分子对接结果,评估受体和配体的相互作用强度和可行性,并为后续的药物设计、分子改造等提供指导和方向。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值