load=1rio

spacefill off
wireframe off
ribbon 0.4
pause

define set1 *.C?' and dna
define set2 *.P, *.OP? and dna
define polar nitrogen, oxygen
define nonpolar carbon, phosphorus, sulfur

define minor ((da and (*.C2 or *.N3 or *.C4 or *.N9)) or (dt and (*.N1 or *.C2 or *.O2)) or (dg and (*.C2 or *.N2 or *.N3 or *.C4 or *.N9)) or (dc and (*.N1 or *.C2 or *.O2))) and dna

define major ((da and (*.N1 or *.C6 or *.N6 or *.C5 or *.N7 or *.C8)) or (dt and (*.N3 or *.C4 or *.O4 or *.C5 or *.C7 or *.C6)) or (dg and (*.C5 or *.C6 or *.O6 or *.N1 or *.N7 or *.C8)) or (dc and (*.N3 or *.C4 or *.N4 or *.C5 or *.C6))) and dna

select protein and within(3.5, set1) and polar
color [197,159,201]
spacefill 0.7
echo "polar with desoxyribose"

pause

color cpk
spacefill off

select protein and within(4.5, set1) and nonpolar
color [242,156,163]
spacefill 0.7
echo "nonpolar with desoxyribose"

pause

color cpk
spacefill off

select protein and within(3.5, set2) and polar
color [227,210,111]
spacefill 0.7
echo "polar with phosphate"

pause

color cpk
spacefill off

select protein and within(4.5, set2) and nonpolar
color [107,241,120]
spacefill 0.7
echo "nonpolar with phosphate"

pause


color cpk
spacefill off

select protein and within(3.5, minor) and polar
color [106,213,203]
spacefill 0.7
echo "polar with minor groove"

pause

color cpk
spacefill off

select protein and within(4.5, minor) and nonpolar
color [126,224,129]
spacefill 0.7
echo "nonpolar with minor groove"

pause

color cpk
spacefill off

select protein and within(3.5, major) and polar
color [247,154,211]
spacefill 0.7
echo "polar with major groove"

pause

color cpk
spacefill off

select protein and within(4.5, major) and nonpolar
color [231,187,65]
spacefill 0.7
echo "nonpolar with major groove"

pause