#This script is for display of GCN4, a yeast transcription 
#factor containing a basic region leucine zipper (bZIP) DNA-binding motif.
#From the structure determined by Ellenberger et al., 1992, Cell 71, pp.
# 1223-1237.  Rasmol commands by Hannah Chen, Jmol by Jason Kahn.
#If molecule isn't loaded...
#zap
#load GCN4-1.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)


#--------------------------------------------------------------------------
#View H-bonds between 
#	Arg243d		G32(chain b)
#	NH1(145)  --  	 O6(652)
#	NH2(1456) -- 	 N7(649)
#Due to asymmetry of site, Arg243c can not form H-bonds with corresponding G.
restrict all
restrict arg243d, g32, 31, 33
color white
wireframe 100
select arg243d, g32
spacefill 300
reset
rotate z -174
rotate y -4
rotate x 156
translate x 36
translate y -8
center atomno=1456
zoom 400

color cpk
select atomno=1455, atomno=652
color green
select atomno=1456, atomno=649
color blue
set measure 2
measure 652 1455
measure 649 1456
set echo top
echo Central Arg-Guanine interaction

set fontsize 20
select G32.O5*
label G32
select arg243d.NE
label Arg243

#--------------------------------------------------------------------------