autodataming的个人博客分享 http://blog.sciencenet.cn/u/autodataming

博文

pdb文件的分割脚本splitpdbs.pl

已有 3664 次阅读 2013-10-10 17:03 |个人分类:DrugDesign|系统分类:科研笔记| PDB, 文件分割

splitpdbs.pl

背景:

vina结果构想文件是多个pdbqt,之前把它转成pdb格式的时候没有进行分隔,

而之前的统计氢键信息的脚本不能处理 多状态的pdb文件

所以写了一个分割pdb脚本

==============================================================

#!/usr/bin/perl -w
use strict;

#author: chen zhaoqiang
#usage: perl splitpdbs.pl -f xxx.pdb
#contact: 744891290@qq.com

#the remark must be {MODEL id} AND ENDMDL

if(@ARGV<2)
{
     print "#usage: perl splitpdbs.pl -f xxx.pdb";
     print "notice:  #the remark must be {MODEL id} AND ENDMDLn";
     #$ARGV[1]="ar3_9.pdb";
     exit 0;
}


open FH,$ARGV[1];

my $basename=$ARGV[1];
$basename=~s/\.pdb//;

while(<FH>)
{
     if($_=~/MODEL\s+(\d+)/)
     {
           my $newfile=$basename.'_'.$1.'.pdb';
           open FF,">$newfile";
     }
     print FF $_ ;          
           
}

 



https://blog.sciencenet.cn/blog-950202-731771.html

上一篇:vina中的构象保存为pdb格式
下一篇:批量统计配体和受体氢键的信息
收藏 IP: 202.127.19.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-3-29 01:30

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部