开源分子对接程序rDock使用方法(1)-Docking in 3 steps

news/2024/4/19 14:22:51/文章来源:https://blog.csdn.net/weixin_40192882/article/details/136175761

欢迎浏览我的CSND博客! Blockbuater_drug …点击进入


前言

rDock是一个快速、多功能的开源对接程序,可用于将小分子配体与蛋白质或核酸受体的对接;选用不同的对接模式可以完成考虑受体结合水的分子对接(Docking with explicit waters)以及药效团限制性对接(Docking with pharmacophore restraints),也可以用来做高通量虚拟筛选(HTVS)。
本文介绍 rDock用于受体-配体的标准对接(Docking in 3 steps),为研究其他模式下的分子对接做准备。


一、Docking in 3 steps 标准对接

rDock官网
rDock的介绍、Linux系统上本地安装请参考系列博文:开源分子对接程序rDock的安装及使用流程

rDock 的基本对接步骤及注意事项

包括3步:定义对接体系、产生对接位点和分子对接。
NOTES:

  1. 为受体(蛋白质或核酸)准备MOL2文件:需要考虑到原子类型问题,建议准备一个全原子MOL2文件,rDock会自动删除非极性氢。确保在$RBT ROOT/data/sf/RbtionicToms.rm中定义了任何非标准原子名称和子结构名称,以便正确分配分布式形式电荷。确保Tripos原子类型设置正确。rDock使用Tripos类型来推导其他关键的原子性质,如原子序数和杂化态。rDock MOL2解析器是为读取CCDC/Astex protein.MOL2文件而开发的,因此此验证集是事实上的标准参考。如果您怀疑某个特定的MOL2文件是否适合rDock,则应将其与CCDC/Astex MOL2文件的格式进行比较。关于Mol2格式,可以参考化学分子Mol2文件格式与使用注意事项。
  2. 准备定义系统的.prm。必须定义得受体参数包括SECTION MAPPER和SECTION CAVITY。如果您希望激活对接位点附近末端OH和NH3+基团的采样,请确保您定义了RECEPTOR_FLEX参数。
  3. 使用rbcavity生成对接站点(.as)文件;如果你想使用参考配体腔定义对接位点的方法,需要受体结合位点中的配体结构。
  4. 准备您想要对接的配体SD文件,需要注意SD文件解析的原子类型问题。特别是,要确保形式电荷和形式键级是一致的,这样文件中就不会有价键错误。rDock将报告任何感知到的价态误差,但仍然会完成。请注意,rDock在对接过程中从不采样键长、键角、环构象或不可旋转键,因此用户需要确保初始构象应该是合理的。
  5. 事先运行一个小的测试计算,以检查系统是否定义正确。例如,使用一个小的配体SD文件,仅使用打分选项(-p score.prm)和-T 2选项运行rbdock以生成详细输出用于检查。输出将包括受体原子特性、配体原子特性、灵活性参数、评分功能参数和对接协议参数。满足要求后,进行较大规模的计算。
    在这里插入图片描述

二、 三步对接案例

rDock对接案例输入文件来源:
人雌激素受体α配体结合结构域与拮抗剂配体4-D的复合物,RCSB下载 pdb id 1SJ0 。
receptor文件:下载1SJ0,加氢,加电荷,删除水分子,选中受体结构,保存为1sj0_rec.mol2
ligand文件:选中配体文件,保存为1sj0_ligand.sd

Step 1. 结构文件准备

mkdir 1sj0_workdir
cd 1sj0_workdir

通过prm文件定义对接体系。

prm 文件是 rDock所特有的文件格式,有以下作用:系统定义文件,评分函数定义文件,搜索协议定义文件

以下是ASTEX数据集的.prm文件示例:

RBT_PARAMETER_FILE_V1.00
TITLE 1sj0_DUDRECEPTOR_FILE 1sj0_rdock.mol2
RECEPTOR_FLEX 3.0SECTION MAPPERSITE_MAPPER RbtLigandSiteMapperREF_MOL 1sj0_ligand.sdRADIUS 6.0SMALL_SPHERE 1.0MIN_VOLUME 100MAX_CAVITIES 1VOL_INCR 0.0GRIDSTEP 0.5
END_SECTIONSECTION CAVITYSCORING_FUNCTION RbtCavityGridSFWEIGHT 1.0
END_SECTION

将以上内容保存为cdk2_rdock.prm,受体结构mol 2文件为cdk2_rdock.mol2,位于结合位点的已知配体结合pose的文件为xtal-lig.sd。
使用的时候我们只需要修改以上内容即可。关于.prm文件的注意事项可以参考如下:

  1. 文件的第一行必须是RBT PARAMETER FILE V1.00,前面不能有空格;
  2. 注释行:注释行应在第一列中以#字符开头,前面不能有空格;
  3. 对于关键字TITLE, SECTION或END_SECTION:
    · 关键字必须从第一列开始,前面不能有空格
    ·关键字TITLE 应该在文件中只出现一次,用于提供标题字符串,通过各种脚本显示,如运行rbscreen.pl,关键字后面应该跟一个空格字符,然后是标题字符串,其中可能包含空格。如果标题行出现的次数一次以上,使用最后一次的记录。
    ·关键字SECTION可以出现多次,并且应该始终与结束END_SECTION配对;关键字后面应该跟一个空格字符,然后是节名,它本身可能不包含空格。在.prm文件中,所有节名称都必须是唯一的。所有SECTION / END_SECTION节中的参数名/值对属于该节。
    ·在TITLE和SECTION关键字之后需要有一个空格字符,否则该部分的后续参数将被忽略。
  4. 参数名/值对:参数名称/值对被读取为自由格式的文本,并且可以有前缀、后缀,并由任意空格分隔。这意味着参数名称和值字符串本身不允许包含任何空格。值字符串被解释为适合该参数的数值、字符串或布尔值。布尔值应输入为TRUE或NULL大写字符串。
  5. prm文件不允许TAB出现。

Step 2. 产生对接位点

在这里插入图片描述
以上文件准备就绪,进入到以上文件的目录,用rbcavity命令生成对接空间:

rbcavity -was -d -r 1sj0_rdock.prm

在这里插入图片描述使用-d参数将生成网格“.grd”文件。该文件可以在pymol中查看:

pymol 1sj0_rdock.mol2 1sj0_ligand.sd 1sj0_rdock_cav1.grd

在这里插入图片描述
在pymol命令行输入以下:

isomesh cavity, 1sj0_rdock_cav1, 0.99

便于查看,调整了透明度。cavity基本覆盖了配体的空间,就是在这个区域进行对接。如果不合适,可以调整 .prm文件中的参数MIN_VOLUME,GRIDSTEP和MAX_CAVITIES。
在这里插入图片描述

Step 3. 运行分子对接

3.1 检查输入文件

将以上生成文件置于1sj0_workdir文件夹中, 文件内容如下:
在这里插入图片描述

3.2 测试-只进行打分

rbdock -i 1sj0_ligand.sd -o output-score -r 1sj0_rdock.prm -p score.prm -T 2

结果为output-score.sd 文件。输出部分如下:
在这里插入图片描述

3.3 运行rdock

运行配体1sj0_ligand.sd的重对接,可以使用以下命令,对每个配体运行50次:

rbdock -i 1sj0_ligand.sd -o output-rdock -r 1sj0_rdock.prm -p dock.prm -n 50

运行片刻,结果为output-rdock.sd 文件。最后显示“END OF RUN”输出部分如下:
在这里插入图片描述

三、 结果查看

在MOE中查看,可以通过SCORE排序。

moe 1sj0_rdock.mol2 output-dock.sd

rdock的结果:

在这里插入图片描述


总结

本文介绍了rDock的基本对接方法,包括结构文件准备、产生对接位点,运行分子对接三个部分,为研究其他模式下的分子对接打好基础。

参考资料

  1. https://bbdrug.blog.csdn.net/article/details/136050880
  2. https://rdock.github.io/

欢迎浏览我的CSND博客! Blockbuater_drug …点击进入

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.luyixian.cn/news_show_997593.aspx

如若内容造成侵权/违法违规/事实不符,请联系dt猫网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Finetuned Multimodal Language Models Are High-Quality Image-Text Data Filters

Finetuned Multimodal Language Models Are High-Quality Image-Text Data Filters 相关链接:arxiv 关键字:Multimodal Language Models、Image-Text Data Filtering、Fine-tuning、Quality Assessment Metrics、Data Quality 摘要: 我们提出…

算法Day05_707.设计链表

推荐阅读 算法day01_ 27. 移除元素、977.有序数组的平方 算法day02_209.长度最小的子数组 算法day03_ 59.螺旋矩阵II 算法Day04_203.移除链表元素 目录 推荐阅读707.设计链表题目思路解法单链表解法双链表解法 707.设计链表 题目 你可以选择使用单链表或者双链表,设…

【翻译】零信任架构准则(二)Know your architecture

了解你的业务架构 在零信任网络模型中,了解你的用户,设备,服务和数据比以往任何时候都更加重要。为了使零信任架构威力发挥到最大,你需要了解你的架构中的每一个组件。因此,正确了解你的资产是最关键的一步&#xff0…

RNN实现退位减法

文章目录 前言RNNRNN架构图前向传播公式反向传播算法 用RNN实现退位减法代码变量的对应关系 总结 前言 最近深入学习了一下RNN,即循环神经网络。RNN是一类比较基础的神经网络,本文使用的是最基础、最简单的循环神经网络的形式。LSTM也是一种常见的循环神…

Qt 简约美观的动画 摆钟风格 第十季

&#x1f60a; 今天给大家分享一个摆钟风格的加载动画 &#x1f60a; 效果如下: 最近工作忙起来了 , 后续再分享其他有趣的加载动画吧. 一共三个文件 , 可以直接编译运行 //main.cpp #include "LoadingAnimWidget.h" #include <QApplication> #include <Q…

TCP收发——计算机网络——day02

今天主要讲了TCP的收发 TCP发端步骤 ①socket ②connect ③send ④closeTCP收端步骤 ①socket ②bind ③listen ④accept ⑤recv ⑥clise其函数主要有 connect int connect(int sockfd, const struct sockaddr *addr,socklen_t addrlen);功能:发送链接请求参数:sockfd:套接…

【论文阅读】(2024.03.05-2024.03.15)论文阅读简单记录和汇总

(2024.03.05-2024.03.15)论文阅读简单记录和汇总 2024/03/05&#xff1a;随便简单写写&#xff0c;以后不会把太详细的记录在CSDN&#xff0c;有道的Markdown又感觉不好用。 目录 &#xff08;ICMM 2024&#xff09;Quality Scalable Video Coding Based on Neural Represent…

Mac电脑可以玩《幻兽帕鲁》吗? 幻兽帕鲁配置要求

Mac电脑可以玩《幻兽帕鲁》吗? 各位只有苹果电脑的玩家也不用担心&#xff0c;现在你也可以在Mac上玩《幻兽帕鲁》了!系统兼容神器CrossOver在第一时间就支持了《幻兽帕鲁》&#xff0c;不用虚拟机就能流畅运行。除了《幻兽帕鲁》外&#xff0c;CrossOver还支持《赛博朋克207…

HarmonyOS NEXT应用开发案例——阻塞事件冒泡

介绍 本示例主要介绍在点击事件中&#xff0c;子组件enabled属性设置为false的时候&#xff0c;如何解决点击子组件模块区域会触发父组件的点击事件问题&#xff1b;以及触摸事件中当子组件触发触摸事件的时候&#xff0c;父组件如果设置触摸事件的话&#xff0c;如何解决父组…

6个维度分析实时渲染和Webgl技术异同

在日常交流中&#xff0c;对Webgl技术熟悉的合作伙伴&#xff0c;在初次了解实时渲染技术时&#xff0c;都会问二者之间的异同。目前很多要求B/S架构的项目&#xff0c;很多在用webgl技术路线&#xff0c;而且这个方案在行业里比较普&#xff0c;业主方对这个也比较熟悉&#x…

重要通告 | 公司更名为“浙江实在智能科技有限公司”

更名公告 升级蜕变、砥砺前行 因业务快速发展和战略升级&#xff0c;经相关政府机构批准&#xff0c;自2024年3月1日起&#xff0c;原“杭州实在智能科技有限公司”正式更名为“浙江实在智能科技有限公司”。 更名后&#xff0c;公司统一社会信用代码不变&#xff0c;业务主体…

蓝桥杯备赛之二分专题

常用的算法二分模板 1. 在数组a[]中找大于等于x的第一个数的下标 //int ans lower_bound(a, a n, x) - a //相当于下方 int l 0, r n - 1; while(l < r) {int mid l r >> 1;if(a[mid] > x) r mid;else l mid 1; } cout << r;2. 在数组a[]中找大于…

nginx,php-fpm

一&#xff0c;Nginx是异步非阻塞多进程&#xff0c;io多路复用 1、master进程&#xff1a;管理进程 master进程主要用来管理worker进程&#xff0c;具体包括如下4个主要功能&#xff1a; &#xff08;1&#xff09;接收来自外界的信号。 &#xff08;2&#xff09;向各worker进…

php常见的45个漏洞及解决方案

[TOC]&#xff08;太多了&#xff0c;目录只列出最重要的几个&#xff0c;剩下的同学们自己翻&#xff09; PHP作为一种广泛应用的服务器端脚本语言&#xff0c;在历史上曾曝出过多种安全漏洞。以下是一些PHP漏洞的类别及其简要解释&#xff0c;以及如何解决这些问题&#xff1…

vue+nodejs超市购物商城在线选品系统wtk87

在此基础上&#xff0c;结合现有超市货品信息管理体系的特点&#xff0c;运用新技术&#xff0c;构建了以 vue为基础的超市货品信息管理信息化管理体系。首先&#xff0c;以需求为依据&#xff0c;根据需求分析结果进行了系统的设计&#xff0c;并将其划分为管理员和用户二种角…

c++ 11 新特性 不同数据类型之间转换函数之reinterpret_cast

一.不同数据类型之间转换函数reinterpret_cast介绍 reinterpret_cast是C中的一种类型转换操作符&#xff0c;用于执行低级别的位模式转换。具体来说&#xff0c;reinterpret_cast可以实现以下功能&#xff1a; 指针和整数之间的转换&#xff1a;这种转换通常用于在指针中存储额…

MATLAB | MATLAB版玫瑰祝伟大女性节日快乐!!

妇女节到了&#xff0c;这里祝全体伟大的女性&#xff0c;节日快乐&#xff0c;事业有成&#xff0c;万事胜意。 作为MATLAB爱好者&#xff0c;这里还是老传统画朵花叭&#xff0c;不过感觉大部分样式的花都画过了&#xff0c;这里将一段很古老的2012年的html玫瑰花代码转成MA…

[OpenCv]频域处理

目录 前言 一、频域变换 1.傅里叶变换 2.代码实现 二、频域中图像处理 1.理解数字图片的频谱 2.频域图像处理步骤 3.使用低通滤波器实现图像平滑 4.使用高通滤波器实现图像锐化 三、总结 前言 数字图像处理的方法有两大类&#xff1a;一种是空间域处理法&#xff0c;…

C# LINQ基础

LINQ基础 1. 入门2. 运算符流语法2.1 连续使用查询运算符2.2 使用Lambda表达式2.2.1 Lambda表达式及Func的方法签名2.2.2 Lambda表达式和元素类型2.2.3 自然排序2.2.4 其他查询运算符 3 查询表达式4 延迟执行4.1 重复执行4.2 捕获的变量4.3 延迟加载的工作原理4.4 查询语句的执…

Python学习笔记-Flask实现简单的抽奖程序

1.导入flask包和randint包 from flask import Flask,render_template from random import randint 2.初始化 Flask 应用: app Flask(__name__) 3. 定义英雄列表 hero [黑暗之女,狂战士,正义巨像,卡牌大师,德邦总管,无畏战车,诡术妖姬,猩红收割者,远古恐惧,正义天使,无极剑…