分子对接|使用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. 结果分析

相关推荐
吴声子夜歌15 小时前
Java数据结构与算法——基本数学问题
java·开发语言·windows
这儿有一堆花18 小时前
用原生脚本编写无害恶作剧
windows
因我你好久不见18 小时前
Windows部署springboot jar支持开机自启动
windows·spring boot·jar
夜流冰19 小时前
Excel - MS Support for Excel: 2 Collaborate
数据库·windows·excel
林瞅瞅19 小时前
PowerShell 启动卡顿?内存飙升?原来是 800MB 的历史记录在作祟!
windows
Shepherd061920 小时前
【Windows Server 实战】WAC 反向代理配置
windows
云小逸20 小时前
【windows系统编程】第一章 Windows 系统核心架构与基础概念
windows·架构
怣疯knight1 天前
Docker Desktop 4.55.0版本安装成功教程
windows·docker
liulilittle1 天前
VEthernet 框架实现 tun2socks 的技术原理
网络·windows·c#·信息与通信·通信
独钓寒江雨1 天前
win11在安全模式下删除360tray.exe
windows·电脑