def main():
## trjconv -f 2LZM_wt.trr -o 2LZM_wt_trjconv.pdb -s 2LZM_wt_MD.tpr -skip 10 -sep
fd = open('2lzm.pdb','r')
lines1 = fd.readlines()
fd.close()
lines2 = []
for line in lines1:
record = line[:6].strip()
if record in ['HELIX','TURN','SHEET',]:
lines2 += [line]
for i in range(50):
fd = open('%i.pdb' %(i),'r')
lines = fd.readlines()
fd.close()
lines = lines2+lines
fd = open('nma%i.pdb' %(i),'w')
fd.writelines(lines)
fd.close()
stop
import sys, os, numpy
sys.path.append('/home/people/tc/svn/tc_sandbox/quakes/')
import phipsi_comparison
sys.path.append('/home/people/tc/svn/tc_sandbox/misc/')
import phipsi_plot
sys.path.append('/home/people/tc/svn/Protool/')
import geometry
instance_geometry = geometry.geometry()
d_residues = {
'ALA':'A','CYS':'C','ASP':'D','GLU':'E','PHE':'F',
'GLY':'G','HIS':'H','ILE':'I','LYS':'K','LEU':'L',
'MET':'M','ASN':'N','PRO':'P','GLN':'Q','ARG':'R',
'SER':'S','THR':'T','VAL':'V','TRP':'W','TYR':'Y',
}
d_max = {}
d_ramachandran = {}
for res in d_residues.values():
d_ramachandran[res] = {}
for phi in range(-180,180):
d_ramachandran[res][phi] = {}
for psi in range(-180,180):
d_ramachandran[res][phi][psi] = 0
d_ramachandran = phipsi_comparison.parse_dihedrals(d_residues) ## slow
## d_ramachandran[res][phi][psi] = 1
for res in d_ramachandran.keys():
max = 0
for phi in d_ramachandran[res].keys():
for psi in d_ramachandran[res][phi].keys():
count = d_ramachandran[res][phi][psi]
if count > max:
max = count
d_max[res] = max
option = 'gromacs'
option = 'pdbs'
if option == 'frames':
l = range(50)
elif option == 'pdbs':
l = ['150l',]
chain = 'D'
elif option == 'gromacs':
l = []
for i in range(101):
l += ['2LZM_wt_trjconv_%i' %(i)]
chain = ' '
fd = open('2LZM_wt_trjconv_0.pdb','r')
lines = fd.readlines()
fd.close()
l_coords_ref = []
for line in lines:
record = line[:6].strip()
if record == 'ATOM':
atom_name = line[12:16].strip()
if atom_name == 'CA':
x = float(line[30:38])
y = float(line[38:46])
z = float(line[46:54])
coord = numpy.array([x,y,z,])
l_coords_ref += [coord]
d_transformations = {}
for pdb in l:
if pdb == '2LZM_wt_trjconv_0.pdb':
continue
fd = open('%s.pdb' %(pdb),'r')
lines = fd.readlines()
fd.close()
l_coords = []
for line in lines:
record = line[:6].strip()
if record == 'ATOM':
atom_name = line[12:16].strip()
if atom_name == 'CA' and line[21] == chain:
x = float(line[30:38])
y = float(line[38:46])
z = float(line[46:54])
coord = numpy.array([x,y,z,])
l_coords += [coord]
l_coords_ref = l_coords_ref[:len(l_coords)]
rmsd = instance_geometry.superpose(l_coords_ref,l_coords)
tv1 = instance_geometry.fitcenter
rm = instance_geometry.rotation
tv2 = instance_geometry.refcenter
d_transformations[pdb] = [tv1,rm,tv2,]
for pdb in l:
print pdb
if option == 'gromacs':
## i = int(pdb[pdb.rindex('_')+1:])
prefix = 'rasmol%s' %(i)
else:
prefix = 'rasmol%s' %(pdb)
lines_rotated = []
## write rasmol script
lines_rasmol = [
'rasmol -nodisplay %s_rotated.pdb << EOF\n' %(pdb),
'cartoon\n',
'wireframe 0\n',
]
d_phipsi, lines_gnuplot = phipsi_plot.calculate_phipsi(pdb, chain = chain)
fd = open('%s.pdb' %(pdb))
lines1 = fd.readlines()
fd.close()
for line1 in lines1:
record = line1[:6].strip()
if record in ['HELIX','TURN','SHEET',]:
lines_rotated += [line1]
if record != 'ATOM':
continue
if line1[21] != chain:
continue
atom_name = line1[12:16].strip()
## if atom_name != 'CA':
## continue
if atom_name not in ['N','CA','C','O',]:
continue
bfactor = float(line1[60:66])
res_no = int(line1[22:26])
res_name = line1[17:20].strip()
if res_no not in d_phipsi.keys():
bfactor = 50.
continue
else:
phi = int(d_phipsi[res_no][0])
psi = int(d_phipsi[res_no][1])
if phi == 180.:
phi = -180.
if psi == 180.:
psi = -180.
bfactor = (
100
*d_ramachandran[d_residues[res_name]][phi][psi]
/float(d_max[d_residues[res_name]])
)
x = float(line1[30:38])
y = float(line1[38:46])
z = float(line1[46:54])
coordinate = numpy.array([x,y,z,])
if pdb != '2LZM_wt_trjconv_0':
tv1 = d_transformations[pdb][0]
rm = d_transformations[pdb][1]
tv2 = d_transformations[pdb][2]
coordinate = numpy.dot(coordinate-tv1,rm)+tv2
lines_rotated += ['%s%8.3f%8.3f%8.3f%s%6.2f%s' %(
line1[:30],
coordinate[0],coordinate[1],coordinate[2],
line1[54:60],bfactor,line1[66:]
)]
h = 159.+80*bfactor/100.
s = 240.
l = 120.+120.*(50-abs(bfactor-50))/50.
r,g,b = hsl2rgb(h,s,l,)
lines_rasmol += [
'select %i\n' %(res_no),
'color [%i,%i,%i]\n' %(
r, g, b,
)
]
fd = open('%s_rotated.pdb' %(pdb),'w')
fd.writelines(lines_rotated)
fd.close()
lines_rasmol += [
## 'ribbons\n',
'rotate x 100\n',
'rotate z 30\n',
'rotate x 100\n',
'rotate z 90\n',
'rotate x 40\n',
'rotate y -20\n',
'write %s.ppm\n' %(prefix),
'exit\n',
]
## write rasmol script to file
fd = open('%srasmol.src' %(prefix),'w')
fd.writelines(lines_rasmol)
fd.close()
## execute rasmol script
os.system('source %srasmol.src > %srasmol.log' %(prefix,prefix))
## convert rasmol output
## os.system('convert %s.ppm -resize x80 %s.gif' %(prefix,prefix))
os.system('convert %s.ppm %s.gif' %(prefix,prefix))
## clean up
os.remove('%s.ppm' %(prefix))
os.remove('%srasmol.log' %(prefix))
os.remove('%srasmol.src' %(prefix))
## os.remove('frame%s.pdb' %(str(i).zfill(2)))
if option == 'gromacs':
line = 'convert '
for i in range(101): ## 0 to 100
line += 'rasmol%s.gif ' %(i)
for i in range(101-1-1,-1+1,-1): ## 99 to 1
line += 'rasmol%s.gif ' %(i)
line += '-loop 0 2lzm_MD.gif'
print line
os.system(line)
## for frame in range(50):
## os.remove('rasmol%s.gif' %(i))
return
def main():
## trjconv -f 2LZM_wt.trr -o 2LZM_wt_trjconv.pdb -s 2LZM_wt_MD.tpr -skip 10 -sep
fd = open('2lzm.pdb', 'r')
lines1 = fd.readlines()
fd.close()
lines2 = []
for line in lines1:
record = line[:6].strip()
if record in [
'HELIX',
'TURN',
'SHEET',
]:
lines2 += [line]
for i in range(50):
fd = open('%i.pdb' % (i), 'r')
lines = fd.readlines()
fd.close()
lines = lines2 + lines
fd = open('nma%i.pdb' % (i), 'w')
fd.writelines(lines)
fd.close()
stop
import sys, os, numpy
sys.path.append('/home/people/tc/svn/tc_sandbox/quakes/')
import phipsi_comparison
sys.path.append('/home/people/tc/svn/tc_sandbox/misc/')
import phipsi_plot
sys.path.append('/home/people/tc/svn/Protool/')
import geometry
instance_geometry = geometry.geometry()
d_residues = {
'ALA': 'A',
'CYS': 'C',
'ASP': 'D',
'GLU': 'E',
'PHE': 'F',
'GLY': 'G',
'HIS': 'H',
'ILE': 'I',
'LYS': 'K',
'LEU': 'L',
'MET': 'M',
'ASN': 'N',
'PRO': 'P',
'GLN': 'Q',
'ARG': 'R',
'SER': 'S',
'THR': 'T',
'VAL': 'V',
'TRP': 'W',
'TYR': 'Y',
}
d_max = {}
d_ramachandran = {}
for res in d_residues.values():
d_ramachandran[res] = {}
for phi in range(-180, 180):
d_ramachandran[res][phi] = {}
for psi in range(-180, 180):
d_ramachandran[res][phi][psi] = 0
d_ramachandran = phipsi_comparison.parse_dihedrals(d_residues) ## slow
## d_ramachandran[res][phi][psi] = 1
for res in d_ramachandran.keys():
max = 0
for phi in d_ramachandran[res].keys():
for psi in d_ramachandran[res][phi].keys():
count = d_ramachandran[res][phi][psi]
if count > max:
max = count
d_max[res] = max
option = 'gromacs'
option = 'pdbs'
if option == 'frames':
l = range(50)
elif option == 'pdbs':
l = [
'150l',
]
chain = 'D'
elif option == 'gromacs':
l = []
for i in range(101):
l += ['2LZM_wt_trjconv_%i' % (i)]
chain = ' '
fd = open('2LZM_wt_trjconv_0.pdb', 'r')
lines = fd.readlines()
fd.close()
l_coords_ref = []
for line in lines:
record = line[:6].strip()
if record == 'ATOM':
atom_name = line[12:16].strip()
if atom_name == 'CA':
x = float(line[30:38])
y = float(line[38:46])
z = float(line[46:54])
coord = numpy.array([
x,
y,
z,
])
l_coords_ref += [coord]
d_transformations = {}
for pdb in l:
if pdb == '2LZM_wt_trjconv_0.pdb':
continue
fd = open('%s.pdb' % (pdb), 'r')
lines = fd.readlines()
fd.close()
l_coords = []
for line in lines:
record = line[:6].strip()
if record == 'ATOM':
atom_name = line[12:16].strip()
if atom_name == 'CA' and line[21] == chain:
x = float(line[30:38])
y = float(line[38:46])
z = float(line[46:54])
coord = numpy.array([
x,
y,
z,
])
l_coords += [coord]
l_coords_ref = l_coords_ref[:len(l_coords)]
rmsd = instance_geometry.superpose(l_coords_ref, l_coords)
tv1 = instance_geometry.fitcenter
rm = instance_geometry.rotation
tv2 = instance_geometry.refcenter
d_transformations[pdb] = [
tv1,
rm,
tv2,
]
for pdb in l:
print pdb
if option == 'gromacs':
## i = int(pdb[pdb.rindex('_')+1:])
prefix = 'rasmol%s' % (i)
else:
prefix = 'rasmol%s' % (pdb)
lines_rotated = []
## write rasmol script
lines_rasmol = [
'rasmol -nodisplay %s_rotated.pdb << EOF\n' % (pdb),
'cartoon\n',
'wireframe 0\n',
]
d_phipsi, lines_gnuplot = phipsi_plot.calculate_phipsi(pdb,
chain=chain)
fd = open('%s.pdb' % (pdb))
lines1 = fd.readlines()
fd.close()
for line1 in lines1:
record = line1[:6].strip()
if record in [
'HELIX',
'TURN',
'SHEET',
]:
lines_rotated += [line1]
if record != 'ATOM':
continue
if line1[21] != chain:
continue
atom_name = line1[12:16].strip()
## if atom_name != 'CA':
## continue
if atom_name not in [
'N',
'CA',
'C',
'O',
]:
continue
bfactor = float(line1[60:66])
res_no = int(line1[22:26])
res_name = line1[17:20].strip()
if res_no not in d_phipsi.keys():
bfactor = 50.
continue
else:
phi = int(d_phipsi[res_no][0])
psi = int(d_phipsi[res_no][1])
if phi == 180.:
phi = -180.
if psi == 180.:
psi = -180.
bfactor = (100 *
d_ramachandran[d_residues[res_name]][phi][psi] /
float(d_max[d_residues[res_name]]))
x = float(line1[30:38])
y = float(line1[38:46])
z = float(line1[46:54])
coordinate = numpy.array([
x,
y,
z,
])
if pdb != '2LZM_wt_trjconv_0':
tv1 = d_transformations[pdb][0]
rm = d_transformations[pdb][1]
tv2 = d_transformations[pdb][2]
coordinate = numpy.dot(coordinate - tv1, rm) + tv2
lines_rotated += [
'%s%8.3f%8.3f%8.3f%s%6.2f%s' %
(line1[:30], coordinate[0], coordinate[1], coordinate[2],
line1[54:60], bfactor, line1[66:])
]
h = 159. + 80 * bfactor / 100.
s = 240.
l = 120. + 120. * (50 - abs(bfactor - 50)) / 50.
r, g, b = hsl2rgb(
h,
s,
l,
)
lines_rasmol += [
'select %i\n' % (res_no),
'color [%i,%i,%i]\n' % (
r,
g,
b,
)
]
fd = open('%s_rotated.pdb' % (pdb), 'w')
fd.writelines(lines_rotated)
fd.close()
lines_rasmol += [
## 'ribbons\n',
'rotate x 100\n',
'rotate z 30\n',
'rotate x 100\n',
'rotate z 90\n',
'rotate x 40\n',
'rotate y -20\n',
'write %s.ppm\n' % (prefix),
'exit\n',
]
## write rasmol script to file
fd = open('%srasmol.src' % (prefix), 'w')
fd.writelines(lines_rasmol)
fd.close()
## execute rasmol script
os.system('source %srasmol.src > %srasmol.log' % (prefix, prefix))
## convert rasmol output
## os.system('convert %s.ppm -resize x80 %s.gif' %(prefix,prefix))
os.system('convert %s.ppm %s.gif' % (prefix, prefix))
## clean up
os.remove('%s.ppm' % (prefix))
os.remove('%srasmol.log' % (prefix))
os.remove('%srasmol.src' % (prefix))
## os.remove('frame%s.pdb' %(str(i).zfill(2)))
if option == 'gromacs':
line = 'convert '
for i in range(101): ## 0 to 100
line += 'rasmol%s.gif ' % (i)
for i in range(101 - 1 - 1, -1 + 1, -1): ## 99 to 1
line += 'rasmol%s.gif ' % (i)
line += '-loop 0 2lzm_MD.gif'
print line
os.system(line)
## for frame in range(50):
## os.remove('rasmol%s.gif' %(i))
return
def main():
import os
lines2 = []
d_residues = {
'ALA': 'A',
'CYS': 'C',
'ASP': 'D',
'GLU': 'E',
'PHE': 'F',
'GLY': 'G',
'HIS': 'H',
'ILE': 'I',
'LYS': 'K',
'LEU': 'L',
'MET': 'M',
'ASN': 'N',
'PRO': 'P',
'GLN': 'Q',
'ARG': 'R',
'SER': 'S',
'THR': 'T',
'VAL': 'V',
'TRP': 'W',
'TYR': 'Y',
}
d_max = {}
import sys
sys.path.append('/home/people/tc/svn/tc_sandbox/quakes/')
import phipsi_comparison
sys.path.append('/home/people/tc/svn/tc_sandbox/misc/')
import phipsi_plot
d_ramachandran = {}
for res in d_residues.values():
d_ramachandran[res] = {}
for phi in range(-180, 180):
d_ramachandran[res][phi] = {}
for psi in range(-180, 180):
## d_ramachandran[res][phi][psi] = 0
##d_ramachandran = phipsi_comparison.parse_dihedrals(d_residues) ## slow
d_ramachandran[res][phi][psi] = 1
for res in d_ramachandran.keys():
max = 0
for phi in d_ramachandran[res].keys():
for psi in d_ramachandran[res][phi].keys():
count = d_ramachandran[res][phi][psi]
if count > max:
max = count
d_max[res] = max
option = 'pdbs'
if option == 'frames':
l = range(50)
elif option == 'pdbs':
l = [
'150l',
]
chain = 'D'
elif option == 'gromacs':
l = []
for i in range(101):
l += ['2LZM_wt_trjconv_%i' % (i)]
for pdb in l:
if option == 'gromacs':
i = int(pdb[pdb.rindex('_') + 1:])
lines2 = []
if option in [
'frames',
'gromacs',
]:
lines2 += [
'HEADER frame t=%2i.000 \n'
% (i)
]
lines2 += [
'MODEL %2i \n'
% (i)
]
prefix = 'rasmol%s' % (i)
## write rasmol script
lines_rasmol = [
'rasmol -nodisplay frame%s.pdb << EOF\n' % (str(i).zfill(2)),
'cartoon\n',
'wireframe 0\n',
## 'ribbons\n',
## 'color temperature\n',
]
d_phipsi, lines_gnuplot = phipsi_plot.calculate_phipsi(pdb,
chain=chain)
fd = open('%s.pdb' % (pdb))
lines1 = fd.readlines()
fd.close()
for line1 in lines1:
record = line1[:6].strip()
if record != 'ATOM':
continue
atom_name = line1[12:16].strip()
## if atom_name != 'CA':
## continue
if atom_name not in [
'N',
'CA',
'C',
'O',
]:
continue
bfactor = float(line1[60:66])
res_no = int(line1[22:26])
res_name = line1[17:20].strip()
if res_no not in d_phipsi.keys():
bfactor = 50.
continue
else:
phi = int(d_phipsi[res_no][0])
psi = int(d_phipsi[res_no][1])
if phi == 180.:
phi = -180.
if psi == 180.:
psi = -180.
bfactor = (100 *
d_ramachandran[d_residues[res_name]][phi][psi] /
float(d_max[d_residues[res_name]]))
lines2 += ['%s%6.2f%s' % (line1[:60], bfactor, line1[66:])]
if option in [
'frames',
'gromacs',
]:
lines2 += [
'TER \n'
]
lines2 += [
'ENDMDL \n'
]
if option == 'pdbs':
fd = open('%s_phipsi_colored.pdb' % (pdb), 'w')
fd.writelines(lines2)
fd.close()
if option in [
'frames',
'gromacs',
]:
h = 159. + 80 * bfactor / 100.
s = 240.
l = 120. + 120. * (50 - abs(bfactor - 50)) / 50.
r, g, b = hsl2rgb(
h,
s,
l,
)
lines_rasmol += [
'select %i\n' % (res_no),
'color [%i,%i,%i]\n' % (
r,
g,
b,
)
]
if option in [
'frames',
'gromacs',
]:
fd = open('%s.pdb' % (pdb), 'w')
fd.writelines(lines2)
fd.close()
lines_rasmol += [
## 'ribbons\n',
'rotate x 100\n',
'rotate z 30\n',
'rotate x 100\n',
'rotate z 90\n',
'rotate x 40\n',
'rotate y -20\n',
'write %s.ppm\n' % (prefix),
'exit\n',
]
## write rasmol script to file
fd = open('%srasmol.src' % (prefix), 'w')
fd.writelines(lines_rasmol)
fd.close()
## execute rasmol script
os.system('source %srasmol.src > %srasmol.log' % (prefix, prefix))
## convert rasmol output
## os.system('convert %s.ppm -resize x80 %s.gif' %(prefix,prefix))
os.system('convert %s.ppm %s.gif' % (prefix, prefix))
## clean up
os.remove('%s.ppm' % (prefix))
os.remove('%srasmol.log' % (prefix))
os.remove('%srasmol.src' % (prefix))
## os.remove('frame%s.pdb' %(str(i).zfill(2)))
if option in [
'frames',
'gromacs',
]:
line = 'convert '
for frame in range(50):
line += 'rasmol%s.gif ' % (i)
for frame in range(50 - 1 - 1, -1 + 1, -1):
line += 'rasmol%s.gif ' % (i)
line += '-loop 0 mode%s.gif' % (mode)
print line
os.system(line)
for frame in range(50):
os.remove('rasmol%s.gif' % (i))
return
def main():
import os
lines2 = []
d_residues = {
'ALA':'A','CYS':'C','ASP':'D','GLU':'E','PHE':'F',
'GLY':'G','HIS':'H','ILE':'I','LYS':'K','LEU':'L',
'MET':'M','ASN':'N','PRO':'P','GLN':'Q','ARG':'R',
'SER':'S','THR':'T','VAL':'V','TRP':'W','TYR':'Y',
}
d_max = {}
import sys
sys.path.append('/home/people/tc/svn/tc_sandbox/quakes/')
import phipsi_comparison
sys.path.append('/home/people/tc/svn/tc_sandbox/misc/')
import phipsi_plot
d_ramachandran = {}
for res in d_residues.values():
d_ramachandran[res] = {}
for phi in range(-180,180):
d_ramachandran[res][phi] = {}
for psi in range(-180,180):
## d_ramachandran[res][phi][psi] = 0
##d_ramachandran = phipsi_comparison.parse_dihedrals(d_residues) ## slow
d_ramachandran[res][phi][psi] = 1
for res in d_ramachandran.keys():
max = 0
for phi in d_ramachandran[res].keys():
for psi in d_ramachandran[res][phi].keys():
count = d_ramachandran[res][phi][psi]
if count > max:
max = count
d_max[res] = max
option = 'pdbs'
if option == 'frames':
l = range(50)
elif option == 'pdbs':
l = ['150l',]
chain = 'D'
elif option == 'gromacs':
l = []
for i in range(101):
l += ['2LZM_wt_trjconv_%i' %(i)]
for pdb in l:
if option == 'gromacs':
i = int(pdb[pdb.rindex('_')+1:])
lines2 = []
if option in ['frames','gromacs',]:
lines2 += ['HEADER frame t=%2i.000 \n' %(i)]
lines2 += ['MODEL %2i \n' %(i)]
prefix = 'rasmol%s' %(i)
## write rasmol script
lines_rasmol = [
'rasmol -nodisplay frame%s.pdb << EOF\n' %(str(i).zfill(2)),
'cartoon\n',
'wireframe 0\n',
## 'ribbons\n',
## 'color temperature\n',
]
d_phipsi, lines_gnuplot = phipsi_plot.calculate_phipsi(pdb, chain = chain)
fd = open('%s.pdb' %(pdb))
lines1 = fd.readlines()
fd.close()
for line1 in lines1:
record = line1[:6].strip()
if record != 'ATOM':
continue
atom_name = line1[12:16].strip()
## if atom_name != 'CA':
## continue
if atom_name not in ['N','CA','C','O',]:
continue
bfactor = float(line1[60:66])
res_no = int(line1[22:26])
res_name = line1[17:20].strip()
if res_no not in d_phipsi.keys():
bfactor = 50.
continue
else:
phi = int(d_phipsi[res_no][0])
psi = int(d_phipsi[res_no][1])
if phi == 180.:
phi = -180.
if psi == 180.:
psi = -180.
bfactor = (
100
*d_ramachandran[d_residues[res_name]][phi][psi]
/float(d_max[d_residues[res_name]])
)
lines2 += ['%s%6.2f%s' %(line1[:60],bfactor,line1[66:])]
if option in ['frames','gromacs',]:
lines2 += ['TER \n']
lines2 += ['ENDMDL \n']
if option == 'pdbs':
fd = open('%s_phipsi_colored.pdb' %(pdb),'w')
fd.writelines(lines2)
fd.close()
if option in ['frames','gromacs',]:
h = 159.+80*bfactor/100.
s = 240.
l = 120.+120.*(50-abs(bfactor-50))/50.
r,g,b = hsl2rgb(h,s,l,)
lines_rasmol += [
'select %i\n' %(res_no),
'color [%i,%i,%i]\n' %(
r, g, b,
)
]
if option in ['frames','gromacs',]:
fd = open('%s.pdb' %(pdb),'w')
fd.writelines(lines2)
fd.close()
lines_rasmol += [
## 'ribbons\n',
'rotate x 100\n',
'rotate z 30\n',
'rotate x 100\n',
'rotate z 90\n',
'rotate x 40\n',
'rotate y -20\n',
'write %s.ppm\n' %(prefix),
'exit\n',
]
## write rasmol script to file
fd = open('%srasmol.src' %(prefix),'w')
fd.writelines(lines_rasmol)
fd.close()
## execute rasmol script
os.system('source %srasmol.src > %srasmol.log' %(prefix,prefix))
## convert rasmol output
## os.system('convert %s.ppm -resize x80 %s.gif' %(prefix,prefix))
os.system('convert %s.ppm %s.gif' %(prefix,prefix))
## clean up
os.remove('%s.ppm' %(prefix))
os.remove('%srasmol.log' %(prefix))
os.remove('%srasmol.src' %(prefix))
## os.remove('frame%s.pdb' %(str(i).zfill(2)))
if option in ['frames','gromacs',]:
line = 'convert '
for frame in range(50):
line += 'rasmol%s.gif ' %(i)
for frame in range(50-1-1,-1+1,-1):
line += 'rasmol%s.gif ' %(i)
line += '-loop 0 mode%s.gif' %(mode)
print line
os.system(line)
for frame in range(50):
os.remove('rasmol%s.gif' %(i))
return
