cmd.read_pdbstr("""\
""", "chain1P")
cmd.hide("all")
cmd.color('grey70', "chain1P")
cmd.show('ribbon', "chain1P")
cmd.select("e4y4o1P1", "c. P & i. 1-149")
cmd.center("e4y4o1P1", state=0, origin=1)
cmd.zoom("e4y4o1P1", animate=-1)
cmd.show_as('cartoon', "e4y4o1P1")
cmd.spectrum('count', 'rainbow', "e4y4o1P1")
cmd.disable("e4y4o1P1")
cmd.show('spheres', 'c. A & i. 3018 | c. A & i. 3020 | c. A & i. 3111 | c. A & i. 3218 | c. A & i. 3232 | c. A & i. 3301 | c. A & i. 3474 | c. A & i. 3504 | c. A & i. 3623 | c. A & i. 3713 | c. A & i. 3755 | c. A & i. 3821 | c. A & i. 4012 | c. A & i. 4013 | c. A & i. 4014 | c. A & i. 4049 | c. A & i. 4050 | c. A & i. 4051 | c. P & i. 201')
util.cbag('c. A & i. 3018 | c. A & i. 3020 | c. A & i. 3111 | c. A & i. 3218 | c. A & i. 3232 | c. A & i. 3301 | c. A & i. 3474 | c. A & i. 3504 | c. A & i. 3623 | c. A & i. 3713 | c. A & i. 3755 | c. A & i. 3821 | c. A & i. 4012 | c. A & i. 4013 | c. A & i. 4014 | c. A & i. 4049 | c. A & i. 4050 | c. A & i. 4051 | c. P & i. 201')