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