拆分 multiple mol2文件

import os

def write_new_mol2(path, context):
    s = ''
    for con in context:
        s += con
    with open(path,'w') as f:
        f.write(s)
        f.close()

path = 'actives_final.mol2'

with open(path, 'r') as f:
    lines = f.readlines()

#full_mol = False
ith_m = 0
mol2_list = ['@<TRIPOS>MOLECULE\n']
for line in lines[1:]:
    if '@<TRIPOS>MOLECULE' not in line:
        mol2_list.append(line)
    else:
        #full_mol = False
        new_file_name = '{0}_{1}_ligand.mol2'.format(ith_m, 'actives')
        write_new_mol2(new_file_name, mol2_list)
        
        mol2_list = ['@<TRIPOS>MOLECULE\n']
        ith_m += 1
        
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值