分子对接|使用Schrödinger中的Glide模块进行分子对接

1. 蛋白准备

1.1 为什么进行蛋白准备?

从PDB晶体结构数据库中下载下来的晶体结构不适合直接用于分子对接。PDB数据库中的晶体结构仅含有重原子,并且可能含有共结晶配体、水分子、金属离子和辅因子等。多聚体结构则可能需要取出单体结构。另外,PDB数据库中的结构可能缺少连接信息,因此必须重新分配连接信息、键序和形式电荷。Schrödinger的Protein Preparation Wizard模块便是用于准备蛋白结构以解决上述问题。

1.2 蛋白准备流程

下面列出了蛋白准备的基本流程,在该流程中假设初始的蛋白结构是PDB格式,里面含有共结晶的配体分子,不含氢原子。在其他的蛋白结构准备过程中需要针对具体问题具体分析,在一些案例中,并不需要把下面所有步骤都走一遍。

(1) 将PDB晶体结构库中的蛋白配体复合物结构导入Maestro中。

(2) 定位需要保留的水分子,删掉其他水分子。通常和金属原/离子有相互作用的以及桥连蛋白和配体的水分子需要保留,其他的水分子删掉。

(3) 简化多聚复合物。如果PDB数据库中蛋白结构是多聚体,查看是否该多聚体有重复的结合位点和重复的链,如果有的话仅保留一个,删掉其他冗余的结构。

(4) 调整蛋白、金属离/原子和辅因子。

修补缺失的残基,如果缺失的残基距离结合位点比较远则一般不影响对接结果,如果缺失的残基在结合位点附近则一定要补全;检查蛋白结构中的金属原/离子;如果有键连着金属离子,则要将键删掉,并调整与金属离子相连原子和金属离子的形式电荷;设置金属原子的电荷并且纠正原子类型;设置辅因子的键序和形式电荷;修正摆向错误基团的方向。

(5) 调整配体的键序和形式电荷。

(6) 如果需要的话,调整蛋白和配体的离子化状态和互变异构状态。

(7) 精炼结构。能量最小化,调整基团摆向。

(8) 检查准备好的结构。

检查整体结构的形式电荷和质子化状态,并根据需要进行调整。检查水分子以及其他基团的取向,如酰胺、羟基的取向等。

1.3 Schrödinger中Protein Preparation Wizard进行蛋白准备 (实际操作)

简介:Protein Preparation Wizard模块可以帮助用户自动进行蛋白准备。经过该模块准备后的蛋白结构可用于后续的其他计算任务,如Glide、Prime、QSite和MacroModel。

Protein Preparation Wizard模块的打开方式如下图:

Protein Preparation Wizard模块的主界面简介如下所示:

(1) 蛋白结构输入

(2) 蛋白结构预处理

(3) 修复蛋白结构的一些问题

问题(Problem)对话框列出了出现原子缺失(通常是侧链)、原子类型错配、原子重叠或原子位置交替等问题的原子和残基。如果检测到任何这些问题,它会在预处理结束时自动打开。

(4) 检查结构并删除不想要的基团

(5) 优化蛋白氢键网络

优化氢键网络主要是重新优化羟基、硫醇、水分子、天冬酰胺和谷氨酰胺中氨基、组氨酸的咪唑环的摆向,预测组氨酸、天冬氨酸和谷氨酸的质子化状态以及组氨酸的互变异构态。这一步的优化是非常有必要的,因为X晶体衍射中并不能确定羟基、硫醇、水分子、天冬酰胺和谷氨酰胺中氨基、组氨酸的咪唑环的摆向。

(6) 氢键网络的交互优化

(7) 能量最小化

(8) 具体操作总结:

2. 对接盒子生成

2.1 设置用于生成盒子的受体参数

2.2 选择对接盒子中心并设置盒子尺寸

2.3 为对接盒子生成设置约束条件

设置约束使 Glide 能够在评估对接适用性的早期即可筛除不符合这些约束的配体、构象或姿势。

2.4 为对接盒子生成设置可旋转基团

2.5 为对接盒子生成设置排除体积

在某些情况下,您可能希望防止配体占据某些空间区域。例如,如果您在已知配体不结合的活性位点附近有一个口袋,您可能希望阻止配体占据该口袋。另一种情况是蛋白质的某些部分缺失,而您想要防止配体占据该区域。

2.6 生成盒子

run即可

2.7 具体操作总结

需要具体问题具体分析,一般只需要进行盒子中心位点选择和盒子尺寸设置即可。其他根据项目的具体情况根据上述参数说明进行参数设置即可。最终可生成glide-grid_1.zip的盒子信息文件。

3. 配体准备

LigPrep是薛定谔中用于准备配体高质量、全原子的三维结构的软件集。LigPrep最简单的用法是去产生单一的低能量并且手性正确的三维结构。LigPrep还可以用于产生各种电离状态、拓扑异构状态、立体化学和环构象的分子。并且还可以i指定条件来过滤分子。

3.1 配体准备LigPrep流程

(1) 转换分子格式:将输入分子的格式 (sdf, mol, mol2) 转换为mae格式。

(2) 加氢

(3) 移除不想要的原子:去除盐

(4) 中和带电基团

(5) 产生质子化状态

(6) 产生互变异构结构

(7) 设置条件过滤结构

(8) 产生正确的手性结构

(9) 产生三维结构并且能量最小化

(10) 转换输出文件格式

3.2 LigPrep实际操作

配体准备完成之后可生成ligprep_1文件夹,文件夹下即配体准备完成之后的文件。

4. 分子对接

薛定谔的Ligand Docking模块用于进行分子对接。

4.1 加载对接盒子

4.2 加载配体

4.3 对接设置

4.4 限制配体对接到指定的核心上

4.5 设置约束条件

4.6 为对接设置扭转约束

4.7 设置输出

4.8 对接具体操作总结

一般需要导入对接盒子和配体,然后设置对接精度,设置输出选项,其他默认即可,具体问题需要具体分析。最后如下设置核数和子任务进行对接即可。

5. 结果分析

相关推荐
hairenjing11236 小时前
使用 Mac 数据恢复从 iPhoto 图库中恢复照片
windows·stm32·嵌入式硬件·macos·word
九鼎科技-Leo8 小时前
了解 .NET 运行时与 .NET 框架:基础概念与相互关系
windows·c#·.net
九鼎科技-Leo11 小时前
什么是 ASP.NET Core?与 ASP.NET MVC 有什么区别?
windows·后端·c#·asp.net·mvc·.net
黎明晓月15 小时前
Java之字符串分割转换List
java·windows·list
九鼎科技-Leo15 小时前
在 C# 中,ICollection 和 IList 接口有什么区别?
windows·c#·.net
顾辰呀16 小时前
实现uniapp-微信小程序 搜索框+上拉加载+下拉刷新
前端·windows
Bunny Chen18 小时前
如何缩小PPT演示文稿的大小?
windows·microsoft·powerpoint
如光照18 小时前
Linux与Windows中的流量抓取工具:wireshark与tcpdump
linux·windows·测试工具·网络安全
wwc_boke19 小时前
Linux查看端口占用及Windows查看端口占用
linux·运维·windows
WangMing_X20 小时前
C# 一个工具类让winform自动根据窗体大小缩放所有控件
开发语言·windows·c#