On May 16, 2011, at 10:49 AM, Daniel Gurnon wrote:
As I understand it, if the script is to work with any pdb structure, I will need two segments of python code:
1) custom residue labeling to display only 1-letter code and number, and
2) something that creates a separate tube for each chain in a model, without connecting one tube to the next.
Hi Dan,
It's actually easier to combine the two, something like this:
from chimera import openModels, Molecule, runCommand
from chimera.resCode import res3to1
for m in openModels.list(modelTypes=[Molecule]):
for chain in m.sequences():
for r in chain.residues:
if not r:
# SEQRES sequence may include missing structure
continue
r.label = "%s %d" % (res3to1(r.type), r.id.position)
runCommand("shape tube %s:.%s@ca radius 0.35 bandLength 2" % (m.oslIdent(), chain.chainID))
--Eric