#!rasmol -script
# File: haem5.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 on

reset
select all
label off
cpk off
ribbon off
wireframe
slab off
rotate z -1
rotate y 2
rotate x -5
translate x 14
translate y -65
zoom 466

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


# Avoid Colour Problems!
select all
colour bonds none
colour backbone none
colour hbonds none
colour ssbonds none
colour ribbons none
colour white


#color by modified cpk
select carbon
color greenblue
select nitrogen
color blue
select oxygen
color red
select iron
color orange

#show all connections, including iron in heme group
select all
connect true

#isolate heme group and important histidines
restrict (hem, his87, his58) and *a
select (hem, his87, his58) and *a
wireframe 40
cpk 120

#centre rotation of molecule
centre atomno=1071

#echo statements

echo Oxygen bound to heme group of haemoglobin
echo Carbon in green, nitrogen in blue, oxygen in red
echo iron in orange
echo Note important histidine residues proximal to the heme group