autodock 结果pdb的生成
Is there a way to save a protein-ligand complex as a PDB file in AutoDock?
Extracting Dockings from DLG Files
AutoDock 4 writes out the coordinates of the atoms in the ligand (and any moving parts of the receptor, if you are doing a flexible sidechain docking). It does not write out the coordinates of the fixed part of the receptor. Each docked conformation is written in PDBQT format in the DLG (docking log file). (Note that AutoDock 3 writes in PDBQ format). So if you did 10 dockings, there should be 10 different docked conformations in the DLG.
Each line of the PDBQT-formatted docked conformation is prefixed by the string "DOCKED: ", so it is possible to extract these lines from the DLG using a couple of simple UNIX commands. You need to go to the UNIX/Linux/Cygwin/Mac OS X Terminal and change directory (cd
) to the directory that contains your DLG, then type the following line at the command line and press < Return > (substitute my_docking.dlg
with the name of your DLG):
grep '^DOCKED' my_docking.dlg | cut -c9- > my_docking.pdbqt
The grep
command is a UNIX command that prints lines that match a pattern. Here, the pattern is ^DOCKED
, and the ^
or caret symbol (Shift-6 on most keyboards) means "at the start of the line", so this pattern matches all lines that begin with the prefix "DOCKED". The |
character (Shift-\ or Shift-backslash) is called a pipe, and it takes the output of the command on its left and feeds it into the input of the command on the right. The cut
command selects portions of each line, and the flag -c9-
means "cut out all the characters after column 9', which has the effect of removing the "DOCKED: " prefix from the line. The last part of the command, > my_docking.pdbqt
, uses the >
redirect command (Shift-. = greater-than symbol) to save the output into a new file called my_docking.pdbqt
.
Converting from PDBQT to PDB
To convert from PDBQT format to PDB format, the simplest thing to do is to remove the charge (Q
) and atom type (T
) columns; this can be achieved using a simple UNIX command. Make sure you are in the same directory where you created my_docking.pdbqt
, and type:
cut -c-66 my_docking.pdbqt > my_docking.pdb
This will create a PDB file, containing all of the docking results. Each docking will appear as a single "MODEL", which is the PDB record usually used to denote an NMR model. Each "MODEL" will contain the ligand and any moving parts of the receptor. If you would like to view the models in this PDB file, you can go load the multi-model PDB file in a program like "PyMol"http://pymol.sourceforge.net/ and then click on the "Play" button to play through all the docked conformations. Click the "Stop" button to halt the play-back, and click on the ">" and "<" buttons to step through the conformations one-at-a-time. It is possible to load the PDB file of the receptor, too, to see how the ligand interacts.
Splitting a Multi-Model PDB File into Separate PDB Files
If you want to split the PDB file that contains all the docked conformations, my_docking.pdb
, into separate PDB files each containing just one docking, then use these commands:
set a=`grep ENDMDL my_docking.pdb | wc -l`
set b=`expr $a - 2`
csplit -k -s -n 3 -f my_docking. mydocking.pdb '/^ENDMDL/+1' '{'$b'}'
foreach f ( mydocking.[0-9][0-9][0-9] )
mv $f $f.pdb
end
For example, if there were 50 ENDMDL records in the file my_docking.pdb, these commands would create 50 separate PDB files numbered from 000 to 049, and they would be named my_docking.000.pdb
, my_docking.001.pdb
, my_docking.002.pdb
and so on, all the way up to my_docking.049.pdb
.
Creating a PDB File of the Complex
To create a single PDB file that contains a complex of both the the receptor and all the models of the docked ligand, you can use the following command to combine the PDB file of the receptor (receptor.pdb
) and all the docked conformations of the ligand stored in 'my_docking.pdb':
cat receptor.pdb my_docking.pdb | grep -v '^END ' | grep -v '^END$' > complex.pdb
This uses the UNIX cat
command which concatenates files together. This command will create a new PDB file called complex.pdb
.
To create a PDB file that contains a complex of both the receptor and a single ligand, then use the commands for splitting the multi-model docked PDB file, select the docked conformation of the ligand, and then use the following command to combine the PDB file of the receptor (receptor.pdb
) and the docked conformation of the ligand; say we chose the ligand conformation in 'my_docking.042.pdb':
cat receptor.pdb my_docking.042.pdb | grep -v '^END ' | grep -v '^END$' > complex.pdb
This uses the UNIX cat
command which concatenates files together, and the grep
command with the -v
flag which extracts all the lines except lines containing the END
record. This command then creates a new PDB file called complex.pdb
, that contains the coordinates of the receptor, followed by all the models of the docked ligand stored in the PDB file my_docking.pdb
.
autodock 结果pdb的生成的更多相关文章
- C#输出到Release VS中Release模式下生成去掉生成pdb文件
Release 与 Debug 的区别就不多说了, 简单来说 Release 优化过, 性能高一些. Debug 为方便调试. 默认情况下是 Debug, 那如何改成 Release 呢? 项目上右键 ...
- pdb文件 小结
.pdb文件,是VS生成的用于调试的符号文件(program database),保存着调试的信息.在VS的工程属性,C/C++,调试信息格式,设置/Zi,那么VS就会在构建项目时创建PDB文件. 在 ...
- PDB文件详解
原文地址:http://blog.csdn.net/feihe0755/article/details/54233714 PDB文件的介绍 PDB(Program Data Base),意即程序的基本 ...
- VS PDB文件
转 PDB文件 PDB(Program Data Base),意即程序的基本数据,是VS编译链接时生成的文件.DPB文件主要存储了VS调试程序时所需要的基本信息,主要包括源文件名.变量名.函数名.FP ...
- VS生成Map文件
一.右键项目属性->链接器->启用增量链接:关闭,选择 否 (/INCREMENTAL:NO) 二.右键项目属性->链接器->调试-> 生成调试信息:是 (/DEBUG ...
- .net开发过程中Bin目录下面几种文件格式的解释
在.NET开发中,我们经常会在bin目录下面看到这些类型的文件: .pdb..xsd..vshost.exe..exe..exe.config..vshost.exe.config 项目发布的时候,往 ...
- VS2010环境下C++工程相关问题汇总
1.链接其他库调试时产生告警: warning LNK4099: 未找到 PDB“vc100.pdb” 解决方案:属性 -> C/C++ -> 输出文件 -> 程序数据库文件名 -& ...
- 如何用VS调试不属于解决方案的EXE和DLL程序
如果你手里有一个现成的EXE, 以及EXE相关联PDB文件, 还有相关联的CPP文件和H文件. 你如何用VS调试? (当然你可以选择WinDbg.不过这里就讨论VS) 你或许想问我干嘛不从一开始就用V ...
- windows下捕获dump之Google breakpad_client的理解
breakpad是Google开源的一套跨平台工具,用于dump的处理.很全的一套东西,我这里只简单涉及breakpad客户端,不涉及纯文本符号生成,不涉及dump解析. 一.使用 最简单的是使用进程 ...
随机推荐
- R语言绘图002-页面布局
par().layout().split.screen()函数 1. par()函数的参数详解 函数par()可以用来设置或者获取图形参数,par()本身(括号中不写任何参数)返回当前的图形参数设置( ...
- 断点续传 (HTTP) 归纳
由于最近项目中要上传较大的文件,基于公司原有的底层框架的局限性,对于大文件的传输都束手无策,基于文件传输的安全性,考虑用断点续传(HTTP)以及FTP上传两种方式实现下面归纳下HTTP续传和FTP上传 ...
- Vue2.0 vue-source.js jsonp demo vue跨域请求
以调用百度的输入提示接口为例 ===================================================================================== ...
- 【BZOJ】【4011】【HNOI2015】落忆枫音
拓扑排序+DP 题解:http://blog.csdn.net/PoPoQQQ/article/details/45194103 http://www.cnblogs.com/mmlz/p/44487 ...
- Model-View-Controller - 杂谈
经典的B/S结构+MVC模式: 注意:MVC层次均位于server中. 下面通过一个简单并且具体的例子来体会下这种MVC模式: --- 功能:查询所有图书:涉及到数据库,JSP,Servlet等 S ...
- MSDN:Code First 迁移
来自:MSDN:Code First 迁移,完全照搬! 本演练将提供对实体框架中 Code First 迁移的概述.您可以完成整个演练,也可以跳至自己感兴趣的主题.主题如下: 启用迁移 生成并运行迁移 ...
- jquery ajax/post/get 传参数给 mvc的action
jquery ajax/post/get 传参数给 mvc的action1.ActionResult Test1 2.View Test1.aspx3.ajax page4.MetaObjec ...
- Spring.net Could not load type from string value
最近有点懒了啊,都没有按时上来博客园更新下,个人觉得遇到难题的时候在这里留下脚印也亦造福他人,进来 晓镜水月 被项目围的团团转,asp.net MVC项目来的,但是我还是不务正业啊,在弄网络爬虫,这个 ...
- PE文件结构详解
(注:最左边是文件头的偏移量.) IMAGE_DOS_HEADER STRUCT { +0h WORD e_magic // Magic DOS signature MZ(4Dh 5Ah) DOS可执 ...
- Gitlab仓库规范实践建议
记录一下Gitlab仓库实践信息: 仓库是指一个可以git clone的地址,用于存储某个服务,模块,子系统或某类辅助代码的地方 仓库的visibility level一般设置为Private(访问需 ...