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