Blake's Space

生物🧬+信息🖥️!

Latest

GroMacs 跑起来!

开始之前感谢一下Justin A. Lemkul, Ph.D.教授的GROMACS Tutorials。现在已经几乎成为新手必学行业标准了。强烈推荐可以去上他的课。这篇blog是我根据他的教程总结出来,并加上了自动化的部分,让一切更加方便! 我们会采用3NPO pdb file 和一个autodock vina对接后的pdbqt文件开始。你可以从下面下载他们,但也可以参考我之前的教程一步一步拿到自己的pdbqt flie: 3NPO3NPO.pdb143 KBdownload-circleKaempferol_dockingKaempferol_docking.pdbqt53 KBdownload-circle 一、准备pdb file 你会发现autodock vina dock之后的pdbqt里面有多个构象,他们都被MODEL分割。如下图: 所以我们处理的第一步就是分割,把这个含有多个model的变成单独的一个model。利用下面几行sh就可以轻松拿到第一个构象。当然,也可以直接用手复制粘贴出来。不过可能使用vina的vina_split会更可靠一些 vina_spl

By Blake Jia
Autodock vina对接!

Autodock vina对接!

相信你已经准备好了两个pdbqt,如果没有可以去看我准备pdbqt的那个blog。如果你想要直接下载,请下载: KaempferolKaempferol.pdbqt2 KBdownload-circle3NPO3NPO.pdbqt122 KBdownload-circle 这个blog采用3NPO作为模拟蛋白,可以在这里看到!上面的pdbqt是处理好的ph=7.0的,Kaempferol就是ligand,3NPO是receptor config file准备 pdbqt是最难准备的,config相对就容易很多。autodock 4的config file准备也很看水平,但是vina把这个过程简化了,就只要最基础的几个事情 # 受体和配体设置 # 应该你已经准备好了两个pdbqt receptor = receptor.pdbqt ligand = ligand.pdbqt out = output.pdbqt # 输出格式也是pdbqt!是包含多个对接结果的pdbqt! # 网格中心与大小 (Grid Box) # 这一部分要看蛋白质的状态! ce

By Blake Jia
Autodock vina の pdbqt准备

Autodock vina の pdbqt准备

这篇blog中假设你毫无autodock经验但是有最最基础linux使用经验,本blog着重描述vina的一些文件和操作流程,以避免以后自己操作的时候犯错。 我们会先从配体准备的软件选择和安装开始,然后会分别讲解ligand和receptor的准备过程。 pdbqt生成-软件的选择和安装 这里的软件指的是从sdf或者pdb文件到pdbqt文件的过程。本质上来说他们都是笛卡尔坐标系xyz的信息储存文件,都可以用文本文档直接打开。 head your_pdbqt.pdbqt 如果你打开看了,就会发现其实pdbqt比pdb多出了两列:也就是q和t。Partial Charge (Q)和Atom Type (T)。Q用于计算静电相互作用,T用于计算van derr waal力和氢键,他们正常pdb是不需要的,但是预测docking的时候非常有用。 一直以来Autodock tools或者MGL tools是比较被推荐的,也是官方标准工具,也是图形化工具。另外也有PyRx等别的第三方pdbqt生成器的库,但是兼容性相对差一点点(生成QT就说白了是很简单的活儿,但是既然有官方的为啥不

By Blake Jia
Autodock vina自动化运行

Autodock vina自动化运行

autodock vina 现在还是比较常用了,虽说config file之类的还是比较直观易懂的,但是一般来说我们做一个东西都要同时dock很多个ligand。下面就是我常用的autodock vina自动化运行脚本。 * 自动化批量对接(Batch Docking): 自动扫描指定目录下的所有配体文件(.pdbqt),并逐个与受体(Receptor)进行对接。 * 智能参数记忆(State Persistence): 首次输入参数后,脚本会自动将其缓存到隐藏文件 .dockenv 中。下次运行时直接回车即可跳过输入,极大提升复用效率。 * 动态生成配置与日志(Auto-Logging): 自动为每个配体生成专属的 config.txt 配置文件,并实时保存带有时间戳的详细运行日志(.out)。 * 重名防覆盖保护(File Protection): 如果输出目录已存在同名对接结果,脚本会自动给新文件加上序号(如 _1, _2),防止实验数据被意外覆盖。 用户输入内容: 序号允许用户输入的内容 (Input)默认值 (Default)说明 (Notes)1

By Blake Jia