CADD之分子对接二：简单autodockvina对接流程——redocking
1、需要软件：autodocktools，vina，pymol
pymol使用教程(包含安装)：基础教程 — PyMOL中文教程 2022.09 文档
autodocktools及vina安装包：
「autodocktools」https://www.aliyundrive.com/s/88pQFmFBNKK 点击链接保存，或者复制本段内容，打开「阿里云盘」APP ，无需下载极速在线查看，视频原画倍速播放。
2、需要蛋白质文件及小分子文件
需要文件为蛋白质文件：.pdb(PDB上下载，或直接在pymol上处理：fetch ####)
小分子配体文件：.mol2格式或.pdbqt格式(利用openbabel进行格式转换)
简单介绍一下openbabel：
openbabel是化学信息学中常用的一款格式转换的工具，用openbabel可实现smi，mol2，pdbqt等多种格式之间的相互转换。
Open Babel由C/C++编写，并提供C ++， Perl， Python等的API接口方便开发。
想深入了解openbabel的，请参考下面的文献：
Open Babel: An open chemical toolbox | Journal of Cheminformatics | Full TextA frequent problem in computational modeling is the interconversion of chemical structures between different formats. While standard interchange formats exist (for example, Chemical Markup Language) and de facto standards have arisen (for example, SMILES format), the need to interconvert formats is a continuing problem due to the multitude of different application areas for chemistry data, differences in the data stored by different formats (0D versus 3D, for example), and competition between software along with a lack of vendor-neutral formats. We discuss, for the first time, Open Babel, an open-source chemical toolbox that speaks the many languages of chemical data. Open Babel version 2.3 interconverts over 110 formats. The need to represent such a wide variety of chemical and molecular data requires a library that implements a wide range of cheminformatics algorithms, from partial charge assignment and aromaticity detection, to bond order perception and canonicalization. We detail the implementation of Open Babel, describe key advances in the 2.3 release, and outline a variety of uses both in terms of software products and scientific research, including applications far beyond simple format interconversion. Open Babel presents a solution to the proliferation of multiple chemical file formats. In addition, it provides a variety of useful utilities from conformer searching and 2D depiction, to filtering, batch conversion, and substructure and similarity searching. For developers, it can be used as a programming library to handle chemical data in areas such as organic chemistry, drug design, materials science, and computational chemistry. It is freely available under an open-source license from http://openbabel.org .
https://jcheminf.biomedcentral.com/articles/10.1186/1758-2946-3-33openbabel生成分子3D构象时，对于有机及类药小分子一般采用MMFF94力场(目前公认最精确的小分子力场)，对于包含任何元素周期表或复杂几何形状的分子(既不受MMFF94支持的)，可以使用UFF(通用力场)力场来代替。GAFF力场也可考虑使用。
下面演示openbabel的安装及简单的使用demo：
# anaconda中openbabel的安装(windows)
conda install openbabel -c conda-forge
# 验证安装成功——显示版本号
obabel -V'''openbabel可用于文件格式转换，过滤(按smarts，标题，描述符值或属性字段), 3D
或2D结构生成，氢从隐式转换显式，并删除小片段或重复结构'''# 最常见使用场景将smiles格式转化为sdf格式
obabel -ismi imatinaib.smi -osdf imatinaib.sdf
# 在转换过程中去掉所有氢
obabel -isdf imatinaib.sdf -osmi iamtinaib.smi -d
# 在转换过程中加上所有氢 ———— 一般不建议在openbabel过程中加氢
obabel -isdf imatinaib.smi -osmi imatinaib.smi -h
# 在转换时添加合适pH=7.4的氢
obabel -isdf imatinaib.sdf -osmi imatinaib.smi -p
# 添加所有氢但不更改文件类型
obabel -isdf imatinaib.sdf -osdf imatinaib.sdf -h
# 一些官能团，例如硝基或砜，可以表示为[N+]([O-])=O或N(=O)=O,以将全部转化为配位键形式
obabel -isdf imatinaib.sdf -osmi imatinaib.smi -b
# 转换分子的子集，可以使用-f和-I定义它们，demo:要转换文件mymols.sdf类型的分子2-4
obabel mymols.sdf -f 2 -l 4 -osdf outputfile.sdf
# 筛选与SMARTS模式匹配的子集，demo：选择出所有含溴苯的分子
obabel mymols.sdf -osdf 'selected.sdf' -s 'c1ccccc1Br'
# 筛选与SMARTS模式不匹配的子集,demo: 选择出不含溴苯的分子
obabel mymols.sdf -osdf 'selected.sdf' -v 'c1ccccc1Br'产生三维结构的时候顺带通过MMFF94高精度有机小分子力场进行优化，demo: 运行如下命令，得到的pdbqt文件里的结构就都是优化过的了
obabel -ismi imatinaib.smi -opdb imatinaib.pdb --gen3d -m --minimize -ff MMFF943、处理大分子蛋白及有机小分子demo
大分子蛋白处理过程：
在处理蛋白前，先做准备工作，在PDB及PDBSum中搜索蛋白结构，下面以有蛋白配体复合物晶体结构的一个样本为例：1SM2与星型孢菌素(STU)结合物
在PDBSum中可以得到一些关于1SM2的信息：
prot-prot中可得到chain A 与chain B是等效的，两条链通过氢键连接到一起，我们在对接时选择其中的一条处理即可。
Ligand中可以得到蛋白与STU作用的方式及位点，可以看到主要相互作用为氢键和疏水作用，STU与Arg(486),Glu(436),Met(438)之间分别有较强的氢键相互作用(键长3埃米以下)。
首先，处理蛋白质分子，打开pymol
# 打开要处理的蛋白结构
fetch 1SM2 ####(蛋白的四位编码，如：1SM2)
# 除去所有非键合水
remove resn HOH
# 删除所有与对接无关的配体分子，金属离子等
# 1SM2两条链是等效的，我们处理其中一条
create 1SM2A, Chain A
# 方便操作，删掉不需要的对象
delete 1SM2
# 将复合物的配体复制为一份单独的对象
create ligand, RESN STU
# 将复合物的蛋白质复制为一份单独的对象
create protein, 1SM2A and not RESN STU
# 保存蛋白为本地PDB文件
save d:/untitled/protein.pdb protein
# 保存配体为本地PDB文件
save d:/untitled/ligand.pdb ligand
根据对蛋白配体结合物的分析，氢键长度范围在2.5-3.0埃米之间，将配体STU周围5埃米残基定义为结合位点并创建对象
# active是选择的命名，sele要选的对象
select active, byres all within 5 of ligand这一步可以对结合位点进行更好的渲染，非必须。
在autodocktools中打开protein.pdb以及ligand.pdb两个文件，进行加氢—合并非极性氢—电荷优化等操作(在电荷优化中，autodock工具中默认的一些金属离子，如Mg，Zn等默认带电荷为零，这里需要手动在.pdbqt文件中进行更改电荷参数)。
加氢：Edit—Hydrogens—Add—OK
合并非极性氢：Edit—Hydrogens—Merge Non—Polar
进行电荷优化：Edit—Charges—Computer Gasteiger (Gasteiger是电荷优化计算的一种，轨道电负性部分均衡(PEOE)电荷，这种方法运用了电负性均衡概念，即电负性不同的原子成键时，电负性较小的原子附近电子密度会流向电负性较大的原子。在这个过程中原先电负性小的原子电负性会增大，当所有原子间电负性相等时，电子密度非分布就是平衡状态分布。实际上PEOE电荷并没有达到完全的电负性均衡，而只达到了部分均衡)
处理蛋白后将其保存成pdbqt格式，方面后面使用。
小分子的处理过程：
redocking过程：什么是redocking，如果你的受体蛋白在PDB数据库里有和其它小分子复合的结构，你可以先做这个小分子对接。然后，把对接的结果和原来的真实结构重叠比较看一下这个对接软件在你的受体蛋白上表现怎么样，是不是可以准确预测结合方式。
对于redocking的过程，只需要处理单个小分子，在ZINC或TCMSP官网下载小分子mol2格式，用pymol打开保存为.pdb格式(或者也可以用从PDB蛋白小分子配合物中提取出来的.pdb文件，都是可以的)。
预处理过程：使用autodocktools预处理小分子。首先，读入小分子：file—read molecule(导入小分子pdb文件)；加氢：Edit—hydrogens—Add；计算电荷：Edit—charge—compute gasteiger；指定为对接小分子：ligand—input—choose—选中小分子—select molecule for autodock4；Ligand—torsion tree—detect root-torsion tree—choose torsions-done；输出：ligand—output—save as pdbqt。
Vina对接过程：
打开处理好的大分子pdbqt结构：Grid—Macromolecule—open(打开大分子.pdbqt文件)。变成二级结构看起来会比较美观：display—Secondary Structure—display二级结构的形式展示。
打开处理好的小分子配体pdbqt结构：Grid—Set Map Types—open ligand(打开你的小分子.pdbqt文件)。
建立对接盒子：Grid—Grid Box，选择盒子中心：Center—Center on ligand—File close saving current，保存对接盒子的参数：Grid—output—Save GPF
写config文件：建立一个config.txt文件，内容参考如下：
receptor = protein.pdbqt
ligand = molecule.pdbqt
# 注释：center部分的内容填写参考GPF
center_x = 24.273
center_y = -4.659
center_z = 28.987
# 注释：size部分的内容填写也参考GPF
size_x = 40
size_y = 40
size_z = 40
energy_range = 4
out = molecule_out.pdbqt
log = molecule_out.pdbqt
# 注释：exhaustiveness是一个精度相关的参数，数值越大，时间开销越大，精准度也越高
exhaustiveness = 24
命令行操作：
vina --config config.txt --out result.pdbqt运行后等待查看结果即可。
对于redocking的过程，主要看对接的均方根偏差RMSD(Root Mean Square Deviation)值与单晶结构中的配体接近程度，越接近，对接能量越低则越好。当对接构象与晶体构象间的重原子RMSD值<=2埃米时，认为对接成功，反之则为失败(所有RMSD值可采用Schrodinger程序重的rmsd.py脚本计算)。
下一篇将介绍批量小分子对单一蛋白靶点的对接。
CADD之分子对接二：简单autodockvina对接流程——redocking相关推荐
- CADD之分子对接三：利用autodockvina套件进行批量对接——virtual drug screening
目录 1.准备需要的对接文件 靶蛋白的pdbqt文件 小分子库(.pdbqt)格式 对接盒子及信息文件config.txt 2.需要用到的工具 linux系统(服务器,虚拟机或windows上的WLS ...
- CADD之分子对接一：背景介绍
在药物筛选过程中,当靶标的结构已知时,可以用计算工具虚拟筛选找到能够调控蛋白质功能的小分子,这种方法叫做基于结构的虚拟筛选(structure-based virtual screening, SBV ...
- 分子对接教程 | (6) AutoDock对接操作与对接结果解读
TCGA | GEO | 文献阅读 | 数据库 | 理论知识 R语言 | Bioconductor | 服务器与Linux 接前文: 分子对接教程 | (1) 软件安装准备 分子对接教程 | (2) ...
- 支付宝支付 第九集：产品数据和支付二维码对接
支付宝支付 第九集:产品数据和支付二维码对接 一.问题 javascript中单引号和双引号的区别是什么? vue中的:作用 v-bind指令 二.代码 目录结构 更新AlipayServiceImp ...
- 【SAP PO】SAP PO 接口配置完整教程之二REST服务对接
SAP PO 接口配置完整教程之二REST服务对接 1.了解服务协议 1.1.服务通讯协议 1.2.具体接口协议 1.3.接口服务测试 2.PO端接口配置 2.1.PO端ESR配置 2.2.PO端IB ...
- 浅谈人行二代征信对接申请及接入验收流程
2018年6月,中国人民银行征信中心发布<中国人民银行征信中心关于抓紧开展与二代征信系统对接有关技术准备的通知>(银征信中心[2018]41号文),通知各分中心组织辖内接入机构在试运行期间 ...
- 泛微E9二次开发,对接金蝶云星空，数据同步，表单同步。
文章目录 泛微E9二次开发,对接金蝶云星空 一.搭建开发环境,引入相关依赖 一.创建项目 二.配置javaSDK 三.配置项目依赖 四.编写测试代码 五.配置编译 六.resin远程debug配置 二 ...
- 柔性对接和半柔性对接的区别_柔性屏风
柔性对接和半柔性对接的区别 I have grown up watching science fiction movies fantasizing about having a watch whose ...
- 移远BC28指令对接电信IOT平台基本流程
目录 概述 一.Coap协议基本流程 1.AT //串口波特率自行匹配 1.AT+NRB //软重启模组 2.AT+ ...
最新文章
- StoryBoard布局注意事项
- SDUT 2401 最大矩形面积
- apply和call
- 听说你做IT的，帮我做个网站呗！
- android 自定义edittext方框样式,Android之EditText自定义边框和边框颜色（转载）
- 游戏试玩站打码平台系统可运营源码
- 2万字硬核剖析网页自定义字体解析（css样式表解析、字体点阵图绘制与本地图像识别等）...
- 为什么低通滤波器的作用相当于积分器
- 备用计算机机房管理制度,计算机机房管理制度
- O2O电子商务盈利模式是什么 O2O电子商务模式遇到的困难有哪些？
- ASTC on Android
- linux免密码互信,Linux建立SSH免密码互信方法
- 黑灰白箱测试+Ubuntu wireshark wifibluetooth
- *python中的字符串
- 网站容器化升级---各模块分别运行一个容器
- 2023年开网店还能赚钱吗？去哪里找货源?
- CGAL例程：地理信息系统----点云数据生成DSM、DTM、等高线和数据分类
- as 运行java_解决AndroidStudio无法运行java中的mian方法问题
- CTSC2016时空旅行
- oslo_messaging使用及测试