#This script is for display of Zif268-DNA by Pavletich and Pabo
#If molecule isn't loaded...
#zap
#load 1aay.pdb
#Script starts here:
define g_major (G and (*.N7 or *.O6 or *.C5 or *.C6 or *.C8 or *.H8))
define a_major (A and (*.N6 or *.N7 or *.C5 or *.C6 or *.2H6 or *.C8 or *.H8)) 
define t_major (T and (*.C5 or *.O4 or *.C4 or *.C5M or *.1H5M or *.2H5M or *.3H5M))
define c_major (C and (*.C5 or *.N4 or *.C4 or *.1H4 or *.2H4 or *.H5))
define g_minor (G and (*.N3 or *.C4 or *.C2 or *.N2 or *.1H2 or *.2H2))
define a_minor (A and (*.N3 or *.C4 or *.C2 or *.H2))
define c_minor (C and (*.C2 or *.O2))
define t_minor (T and (*.C2 or *.O2))
define major_groove (g_major or a_major or c_major or t_major)
define minor_groove (g_minor or a_minor or c_minor or t_minor)

#--------------------------------------------------------------------------
#
#--------------------------------------------------------------------------
# Now highlight the Arg's (and one His) which make specific major groove
# contacts to G's.  Note that all of the recognized G's are on one strand.
#
reset
rotate z 163
rotate y -69
rotate x -102
zoom 131
select zn
color [100,100,100]
select dna and not *.P
color white
wireframe 120
backbone 200
select protein
color [143,255,143]
select g2,g4,g6,g7,g8,g10
color cyan
select (g.o6 or g.n7) and (g2,g4,g6,g7,g8,g10)
color red
spacefill 300
select arg180,arg174,arg146,arg124,arg118
color yellow
wireframe 150
select (arg.NH1 or arg.NH2) and (arg180,arg174,arg146,arg124,arg118)
color blue
spacefill 300
select his149
color violet
wireframe 150
select his149.NE2
color green
spacefill 300
set specular on
set shadow on
#