Ss: Difference between revisions

From PyMOLWiki
Jump to navigation Jump to search
m (changed a loop variable to avoid name clashes.)
mNo edit summary
Line 2: Line 2:


<source lang="python">
<source lang="python">
from pymol import cmd
from pymol import stored
def ss(selection):
def ss(selection):


Line 14: Line 17:


     stored.pairs = []
     stored.pairs = []
     cmd.iterate(selection, "stored.pairs.append((resi, ss))")
     cmd.iterate("%s and n. ca" % selection, "stored.pairs.append((resi, ss))")
     num, currentType = stored.pairs[0]
     num, currentType = stored.pairs[0]



Revision as of 09:53, 24 October 2007

A command to list a summary of the secondary structure for a selection. Use like "ss my_protein" where "my_protein" is the name of the chain or structure in view.

from pymol import cmd
from pymol import stored

def ss(selection):

    class SSE(object):

        def __init__(self, start, typ):
            self.start, self.typ = start, typ
            self.end = -1

        def __repr__(self):
            return "%s-%s %s" % (self.start, self.end, self.typ)

    stored.pairs = []
    cmd.iterate("%s and n. ca" % selection, "stored.pairs.append((resi, ss))")
    num, currentType = stored.pairs[0]

    sses = [SSE(num, currentType)]
    currentSSE = sses[0]
    for resi, ssType in stored.pairs:
        if ssType == currentType:
            currentSSE.end = resi
        else:
            sses.append(SSE(resi, ssType))
            currentSSE = sses[-1]
            currentType = ssType

    for sse in sses:
        print sse

cmd.extend("ss", ss)