cmd.read_pdbstr("""\ """, "chain1Y") cmd.hide("all") cmd.color('grey70', "chain1Y") cmd.show('ribbon', "chain1Y") cmd.select("e5vp21Y1", "c. Y & i. 1-107") cmd.center("e5vp21Y1", state=0, origin=1) cmd.zoom("e5vp21Y1", animate=-1) cmd.show_as('cartoon', "e5vp21Y1") cmd.spectrum('count', 'rainbow', "e5vp21Y1") cmd.disable("e5vp21Y1") cmd.show('spheres', 'c. A & i. 3163 | c. A & i. 4116 | c. Y & i. 501 | c. Y & i. 502 | c. Y & i. 503 | c. Y & i. 504') util.cbag('c. A & i. 3163 | c. A & i. 4116 | c. Y & i. 501 | c. Y & i. 502 | c. Y & i. 503 | c. Y & i. 504')