plot
  window 10;
  
  read mol "had12.pdb";
  map omit "had_12_fofc.omap" format brix around in residue A401 plus 5;
  
  transform all
    by centre position require atom C* and in residue A401
    by rotation x 180.0
    by rotation z 146.100
    by rotation y -7.22053
    by rotation x 46.4024
    
by rotation
  -0.747261 -0.522499 0.410604
  0.490607 -0.0170095 0.871215
  -0.448225 0.85247 0.269053
by rotation
  0.958724 -0.279512 -0.0521659
  0.279898 0.96003 9.50947e-05
  0.0500542 -0.0146923 0.998638  ;
  
colour ball-and-stick by atom;
set maplinewidth 0.4;
set mapcolour blue;
set map2colour blue;
set coverradius 10;

contour in-out omit at 3.0  sigma covering in residue A300;   
colour fractional bonds
          from black via black at 0.04 via white at 0.04 to white at 0.05
        by atom; 


    set atomcolour atom C* green;
    set atomcolour atom N* blue;
    set atomcolour atom O* red;
    !set atomcolour atom FE* magenta;   
    set atomcolour atom CL* yellow;  
    
set atomradius atom * 0.30000;
cpk require atom FE and in residue A300;
set sticksegments 2;
set stickradius 0.15;
set atomradius atom * 0.150000;


cpk in residue A401;
ball-and-stick in residue A401;
cpk in residue A310;
ball-and-stick in residue A310;
end_plot