#!rasmol -script
# File: haem1.spt
# Creator: RasMol Version 2.6, modified by Eleanor Stainback
# Assumes pdb file already loaded in Rasmol Window
background [0,0,0]

set ambient 60
set specular off

reset
select all
slab off
rotate z 6
rotate y -12
rotate x 84
zoom 132

set bonds off
set axes off
set boundingbox off
set unitcell off
set bondmode and
dots off



#color molecule by chain, isolate basic residues surrounding BPG binding site
select all
wireframe
cpk off
ribbon off
color chain
select his2, his143, lys82, arg104
color red
wireframe 70
echo Haemoglobin
echo The molecule is coloured by chain with basic residues surrounding 
echo the BPG binding site in red.