#!rasmol -script
# File: haem4.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
label off
ribbon off
cpk off
wireframe
slab off
rotate z 43
rotate y 11
rotate x -52
translate x -16
translate y 6
zoom 331

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



#show FE bonds
connect true

#color nonpolar residues vs others
select polar
color green
select hydrophobic
color white
cpk 300

#highlight heme group
select hem
color red
wireframe 80
cpk 120
select iron
color orange

#highlight important histidine residues
select his
color blue
wireframe 80
cpk 120

# Bonds
select all
wireframe 80

# Ribbons
ribbons off

# Backbone
backbone off

# Labels
labels off
colour labels none
set fontsize 8
select atomno=439
label "%n%r"
select atomno=648
label "%n%r"

# Monitors
monitors off

select all
ssbonds off
hbonds off


#echo statements
echo Heme group of Haemoglobin coordinated
echo with a proximal histidine residue.
echo nonpolar residues in white, polar in green,
echo histidines in blue, heme in red, iron in orange