from pymol import cmd
import time
cmd.fetch("1ozj")
cmd.show_as("cartoon","all")
time.sleep(3)
cmd.show_as("lines","all")
cmd.remove("(all)-((/1ozj/A/C)+(/1ozj/B/D))")
time.sleep(3)
cmd.show_as("spheres","/1ozj/*/*/*/O4'")
cmd.show_as("spheres","/1ozj/*/*/*/O3'")
time.sleep(3)
cmd.show_as("lines","all")
cmd.show_as("spheres","/1ozj/*/*/*/OP*")
time.sleep(3)
cmd.show_as("lines","all")
cmd.show_as("spheres","/1ozj/*/*/*/N*")