cmd.read_pdbstr("""\
""", "chain16")
cmd.hide("all")
cmd.color('grey70', "chain16")
cmd.show('ribbon', "chain16")
cmd.select("e4z8c161", "c. 6 & i. 2-54")
cmd.center("e4z8c161", state=0, origin=1)
cmd.zoom("e4z8c161", animate=-1)
cmd.show_as('cartoon', "e4z8c161")
cmd.spectrum('count', 'rainbow', "e4z8c161")
cmd.disable("e4z8c161")
cmd.show('spheres', 'c. A & i. 3459 | c. A & i. 4009 | c. A & i. 4021')
util.cbag('c. A & i. 3459 | c. A & i. 4009 | c. A & i. 4021')