#########################################################################
#                                                                       #
# PyMol Script to render chorismate mutase mutase from E. coli          #
# Few active site residues surrounding the bound inhibitor              #
#                                                                       #
# Esther Zhuang and Kalju Kahn, UCSB, 2004                              #
#                                                                       #
#########################################################################

hide 

set stick_radius, 0.1
set sphere_scale, 0.17

select tsa, resn tsa
cmd.disable('tsa')

color gray
color marine, elem n
color salmon, elem o
color white, elem h
color wheat, elem c and tsa

show sticks
hide sticks, resi 52
show sphere, tsa
show sphere, resn hoh
hide sticks, elem h and neighbor elem c and (!tsa)
hide sphere, elem h and neighbor elem c and tsa
hide sticks, elem h and neighbor elem c and tsa

# Show hydrogen bonds
cmd.select('lb',"(1ECM_chor_mut_asite_EC`15)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`48)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`13)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`51)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`14)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`157)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`157)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`15)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`78)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`14)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`76)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`11)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`13)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`111)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`16)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`156)")
dist
cmd.select('lb',"(1ECM_chor_mut_asite_EC`119)")
cmd.select('rb',"(1ECM_chor_mut_asite_EC`16)")
dist
hide label
cmd.delete("lb")
cmd.delete("rb")

# Hide unimportant atoms
cmd.select('lb',"(1ECM_chor_mut_asite_EC`121)")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`133))) and not (((1ECM_chor_mut_asite_EC`133)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`136))) and not (((1ECM_chor_mut_asite_EC`136)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`130))) and not (((1ECM_chor_mut_asite_EC`130)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`129))) and not (((1ECM_chor_mut_asite_EC`129)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`128))) and not (((1ECM_chor_mut_asite_EC`128)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`138))) and not (((1ECM_chor_mut_asite_EC`138)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`153))) and not (((1ECM_chor_mut_asite_EC`153)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`150))) and not (((1ECM_chor_mut_asite_EC`150)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`147))) and not (((1ECM_chor_mut_asite_EC`147)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`145))) and not (((1ECM_chor_mut_asite_EC`145)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`146))) and not (((1ECM_chor_mut_asite_EC`146)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`105))) and not (((1ECM_chor_mut_asite_EC`105)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`104))) and not (((1ECM_chor_mut_asite_EC`104)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`97))) and not (((1ECM_chor_mut_asite_EC`97)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`106))) and not (((1ECM_chor_mut_asite_EC`106)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`112))) and not (((1ECM_chor_mut_asite_EC`112)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`118))) and not (((1ECM_chor_mut_asite_EC`118)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`87))) and not (((1ECM_chor_mut_asite_EC`87)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`88))) and not (((1ECM_chor_mut_asite_EC`88)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`86))) and not (((1ECM_chor_mut_asite_EC`86)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`85))) and not (((1ECM_chor_mut_asite_EC`85)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`53))) and not (((1ECM_chor_mut_asite_EC`53)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`64))) and not (((1ECM_chor_mut_asite_EC`64)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`68))) and not (((1ECM_chor_mut_asite_EC`68)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`72))) and not (((1ECM_chor_mut_asite_EC`72)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`63))) and not (((1ECM_chor_mut_asite_EC`63)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`62))) and not (((1ECM_chor_mut_asite_EC`62)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`37))) and not (((1ECM_chor_mut_asite_EC`37)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`36))) and not (((1ECM_chor_mut_asite_EC`36)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`27))) and not (((1ECM_chor_mut_asite_EC`27)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`46))) and not (((1ECM_chor_mut_asite_EC`46)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`42))) and not (((1ECM_chor_mut_asite_EC`42)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`38))) and not (((1ECM_chor_mut_asite_EC`38)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`54))) and not (((1ECM_chor_mut_asite_EC`54)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`55))) and not (((1ECM_chor_mut_asite_EC`55)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`139))) and not (((1ECM_chor_mut_asite_EC`139)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`140))) and not (((1ECM_chor_mut_asite_EC`140)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`98))) and not (((1ECM_chor_mut_asite_EC`98)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`99))) and not (((1ECM_chor_mut_asite_EC`99)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`28))) and not (((1ECM_chor_mut_asite_EC`28)) in (lb)))")
cmd.select('lb',"( ((lb) or ((1ECM_chor_mut_asite_EC`29))) and not (((1ECM_chor_mut_asite_EC`29)) in (lb)))")
select extra, lb
hide sticks, extra
cmd.delete("lb")
cmd.disable('extra')

# Zoom into the active site
set_view (\
    0.558672965,    0.473125309,    0.681182563,\
    0.803632975,   -0.511842549,   -0.303590119,\
    0.205030248,    0.717038035,   -0.666183352,\
   -0.000440143,   -0.000899255,  -35.670017242,\
   27.960599899,   14.778103828,   49.769565582,\
   23.474666595,   51.443298340,    0.000000000 )