2011-10-07 70 views
0

我想从PDB记录中去除各种配体。 它只是足以删除HET,HETNAM,HETATM ....,即。那些在哪里用3字母代码标识化合物,还是有必要清理其他一些领域?pdb蛋白库格式 - 配体去除

是否有任何python | perl脚本已经为此写作?

回答

1
open(FILE,"file.pdb"); 
@file=<FILE>; 
foreach (@file){ 
if (/^HETATM/){ 
print $_,"\n"; 
}} 

这分离出配体。去除配体,保持不等于正则表达式的前面。