示例#1
0
import build
import rotate
import analyze

m=14
n=9
l=13

path_in='/data/liaot/toolbox/python/protein/patchdock_sample/pdb-for-link/after-model/'
path_out='/data/liaot/toolbox/python/protein/patchdock_sample/rassf-bended-ready-for-solvate/'

for k in range(1,2):
	pdb_filename='%d-%d-%d-%d.pdb' % (m, n, l, k+1)
	protein_obj=build.protein()
	protein_obj.read_residue_info(path_in+pdb_filename)

##### move protein to center of mass
	protein_obj.move_to_com(protein_obj.com)

##### rotate protein to fit insdie a smallest box
##### print vortex to check the box size
	print rotate.find_vortex(protein_obj)
	protein_obj=rotate.minbox(protein_obj)
	print rotate.find_vortex(protein_obj)

	for natom in range(len(protein_obj.resid)):
        	protein_obj.resid[natom]-=204

	pdb_out='rassf-%d-%d.crd' % (m, k+1)
	protein_obj.write_crd(path_out+pdb_out)
	del protein_obj
示例#2
0
文件: main.py 项目: tjliao/vault
# pdb_filename='sarah.pdb'
# protein_obj=build.protein()
# protein_obj.read_residue_info(pdb_filename)

##### move protein to center of mass
# protein_obj.move_to_com(protein_obj.com)

##### rotate protein to fit insdie a smallest box
##### print vortex to check the box size
# print rotate.find_vortex(protein_obj)
# protein_obj=rotate.minbox(protein_obj)
# print rotate.find_vortex(protein_obj)


# superimpose two proteins
target_filename = "mra.pdb"
protein_target = build.protein()
protein_target.read_residue_info(target_filename)

arrow_filename = "ra.pdb"
protein_arrow = build.protein()
protein_arrow.read_residue_info(arrow_filename)

select_arrow = ["P", 205, 366]
protein_arrow = rotate.superimpose(protein_arrow, select_arrow, protein_target)
# protein_arrow=rotate.superimpose(protein_arrow, select_arrow, protein_target, select_target)

# protein_obj.write_pdb('mra-chain.pdb')
# protein_obj.write_crd('test.crd')