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

相关推荐
web守墓人1 小时前
【goed/ui】自定义组件设计思想篇
linux·windows·ui·golang
2601_962364975 小时前
红外活体模组,Windows 设备安全身份入口解析
windows·安全·电脑
易点云徐源5 小时前
(已解决)Photoshop 2020 点击“选择主体”闪退:不用重装,1 个配置切回旧算法解决
windows·photoshop
云计算练习生6 小时前
什么是系统调用?为什么程序访问硬件必须经过它
linux·windows·操作系统·系统调用·操作系统原理
其实防守也摸鱼13 小时前
DeepSeek Harness 开源贡献手记:从入门到合入主线
服务器·数据库·windows·https·ssl
qqww15514 小时前
Win11更新后Edge卡顿掉帧?隐私窗口流畅解决办法
windows·edge·优化·卡顿
其实防守也摸鱼15 小时前
堆叠注入(Stacked Injection)详解
服务器·数据库·windows·https·ssl
郝学胜-神的一滴15 小时前
C++20模板元编程 05:吃透变参模板,解锁编译期万能参数能力
服务器·开发语言·c++·windows·vscode
是店小二呀17 小时前
鸿蒙PC_Dart-Sass-ohos适配全记录
linux·运维·windows
2601_962364971 天前
不止指纹和人脸:Windows Hello如何重塑PC生物识别安全?
windows·安全·电脑