通过smina对接产生的小分子mol2文件中,每一个pose的默认名字都为"*****",用起来很不方便,按照ranking的排序将每一个小分子文件重命名为“hit000“+ranking顺序:
def mol2_rename(mol2):
f = open(mol2, 'r+')
lines = f.readlines()
f1 = open(mol2, 'w')
count = 0
for lines in lines:
if "*****" in lines:
count +=1
counts = str(count).zfill(4) #
f1.write(lines.replace("*****", f'hit{counts}'))
else:
f1.write(lines)
f.close()
f1.close()