Ss

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.

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(selection, "stored.pairs.append((resi, ss))")
    num, currentType = stored.pairs[0]

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

    for sse in sses:
        print sse

cmd.extend("ss", ss)