Bell: 5 revisions

2014-03-28T03:36:28Z

5 revisions

Newacct at 08:14, 27 November 2009

2009-11-27T08:14:58Z

New page

== Overview ==
This script will write the secondary structural elements for each PDB file in some specified directory, for each alpha carbon in the protein, to an output text file. Residues without secondary structure definition will get a "." in the string.

This is untested code and has little chance of working for anyone but me. Although, it does work for me.

== Instructions ==
# Copy the code to your machine
# Find '''FIXME''' in the code below and change "files" to either a glob like the one already there (for a whole directory) or a list with one element (for just one file).
# change the output filename. By default it's '''PDBCODE.ss'''.

== The Code ==
<source lang="python">
import os
import os.path
import glob

from pymol import cmd
from pymol import stored
from pymol import selector

files = glob.glob("/tmp/thy_model/*.pdb")

for file in files:
pdbName = os.path.basename(file).split(".")[0]
cmd.load(file, pdbName)
outFile = open(pdbName + '.ss', 'wb')
stored.ss = ""
cmd.iterate( '(n. CA)', 'stored.ss = stored.ss + ("%1s"%ss)')
for c in stored.ss:
if c == " ":
outFile.write('.')
else:
outFile.write(c)
cmd.delete(pdbName)
outFile.close()
</source>

== Example ==
<source lang="python">
python ss.pym # my directory has 1d7p.pdb in it
</source>
output for 1D7P.pdb:

...............HHH.......SSS..SSSHHHHH..................................................................SS..SS..................................................

[[Category:Script_Library|Write SS]]
[[Category:Structural_Biology_Scripts]]

← Older revision	Revision as of 03:36, 28 March 2014
(No difference)

WriteSS - Revision history

Bell: 5 revisions

Newacct at 08:14, 27 November 2009