薛定谔 | 分子对接及基于受体的虚拟筛选

本文详细介绍了蛋白质受体的准备过程,包括使用ProteinPreparation模块处理可变残基,以及通过ReceptorGridGeneration模块生成受体分子格点文件。此外,还涉及到了LigPrep模块下的配体准备和Ligand-Docking模块的对接操作。最后,讨论了如何评估对接结果,并强调了分数负值与对接效果的关系。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

1.准备受体

最好直接从File >> Get PDB输入PDB ID获取。因为原始pdb文件保留了蛋白完整序列信息,这样在准备蛋白时薛定谔可以自动补全缺失残基。

选择Protein Preparation模块。按照下图操作:

在正常准备蛋白的步骤中弹出了蛋白质中一些位置可变的残基(如下图所示),选中任意一种残基,点击next position会展现出这个残基可能的其他构象,点击commit即选定这个位置,表格中的数字表示这种位置的残基可能的概率大小,当所有残基都commit后,再点击update,将选定的这一组氨基酸的定位给固定下来(当然每个残基要选哪一个残基要靠自己的判断进行适当的取舍,对于不在活性口袋的可变残基可忽略其可变的位置。)


2. 生成受体分子格点文件

选择receptor grid generation模块。鼠标点击蛋白-配体复合物中的配体分子即可生成下图所示的格点盒子。

  • 选小配体分子时,需要在显示视图中选中ligand,要让其出现格子才算选上。
  • 在site项中设定盒子中心,在高级选项中可以调整小盒子的大小。
  • 点击对话框页面下方的文件名右边的小齿轮图标,可以设置生成格点文件的存放路径。
  • Constraintstab栏项:可以添加一些原子或残基,在对接过程中他们的位置不会变。
  • Rotatable Groupstab栏项:系统会列出可以旋转的残基以及他们可旋转的键,感觉作用不大,基本不会用上。
  • Exclusive Volumestab栏项:可以选中一些原子或残基,系统会认定这些选中的基团周围的球状体积都是不能被配体分子占据的。

3. 准备配体

选择LigPrep模块


4. Docking

选择Ligand-docking模块


5.查看结果

在右上角table中查看打分结果,分数越负对接效果越好。

### Schrödinger Molecular Docking 和 Electrostatic Potential Map 教程 #### 使用Schrödinger软件包进行分子对接 为了执行分子对接研究,通常采用Glide模块作为核心工具之一。该模块能够模拟配体与受体之间的相互作用过程,并评估其结合模式和亲和力。 1. **准备结构文件** - 需要准备好蛋白质(PDB格式)以及小分子(SDF或MOL2格式)的三维坐标文件。 - 利用PrepWizard功能处理蛋白结构,优化氢键网络并修复可能存在的缺陷[^1]。 2. **定义活性位点** - 可通过指定已知底物位置来限定搜索范围;也可以基于残基列表手动划定区域。 3. **设置参数选项** - 对于初步筛选可选用HTVS(high-throughput virtual screening),后续验证阶段推荐SP(standard precision) 或 XP(extra precision)。 4. **运行对接计算** - 提交任务至Maestro界面中的Jobs面板等待完成即可获得预测构象集合及其打分情况。 5. **分析结果输出** - 查看最佳姿态对应的能量值、根均方偏差(RMSD)等指标判断合理性; - 运用LigPlot+绘制二维交互图辅助理解关键接触点分布特征。 ```bash # 执行命令行形式的Glide对接流程 (假设环境变量配置完毕) $ schr_run glide -i input_file.maegz -o output_folder/ ``` #### 创建静电势能映射图像 对于可视化展示生物大分子表面特性而言,EPik/Prime QM Polarized Force Field联合方案提供了高效途径: - 启动Maestro图形化操作平台加载目标对象后,在属性栏里选取“Electrostatics”标签页下的相应按钮启动计算进程。 - 用户可以自定义网格分辨率、截断距离等细节设定以满足特定需求。 - 完成之后将以彩色编码方式呈现正负电荷聚集区,便于直观识别潜在药物靶向部位。 ```python from pymol import cmd, stored def show_ep_map(selection='all'): """显示选定原子的选择性静电势""" cmd.set('surface_quality', 0) cmd.show('surf', selection=selection) cmd.spectrum('q', 'rainbow', selection) cmd.extend('show_ep_map', show_ep_map) ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

药研猿

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值