select all
color grey
cartoons off
wireframe 40
select water
cpk off


pause
define major_groove_p (*.N4 and DC) or (*.O4 and DT) or (*.N6, *.N7 and DA) or (*.O6, *.N7 and DG)  
define major_groove_n (*.C6, *.C5, *.C5 and DA) or (*.C6, *.C5, *.C8 and DG) or (*.C6, *.C5, *.C4 and DC) or (*.C6, *.C5, *.C4, *.C7 and DT) 
define minor_groove_p (*.N3, *.N9 and DA) or (*.N2, *.N3, *.N9 and DG) or (*.O2, *.N1 and DC) or (*.O2, *.N1 and DT)
define minor_groove_n (*.C2, *.C4 and DA) or (*.C2, *.C4 and DG) or (*.C2 and DC) or (*.C2, *.C6 and DT)
define drib_p *.O? and dna
define drib_n *.C? and dna
define phosphate_p *.OP? and dna
define phosphate_n *.P and dna

pause
select within (3.5, drib_p) and protein and (*.N? or *.O?)
cpk 200
color [112,46,176]

pause
select within (4.5, drib_n) and protein and (*.C? or  *.S? or *.P?)
cpk 200
color [0,201,0]

pause
select within (3.5, phosphate_p) and protein and (*.N? or *.O?)
cpk 200
color [31,255,199]

pause
select within (4.5, phosphate_n) and protein and (*.C? or  *.S? or *.P?)
cpk 200                                                                  
color [235,0,38]

pause
select within (3.5, major_groove_p) and protein and (*.N? or *.O?)                                                                 
cpk 200                                                     
color [255,255,48]

pause
select within (4.5, major_groove_n) and protein and (*.C? or  *.S? or *.P?)  
cpk 200                                                                
color [255,13,13]

pause
select within (3.5, minor_groove_p) and protein and (*.N? or *.O?)                                                             
cpk 200                                                       
color [255,128,0]

pause
select within (4.5, minor_groove_n) and protein and (*.C? or  *.S? or *.P?)                                                   
cpk 200                                                                   
color [230,230,230]