0    0    2

Topology file for haem for the transition state for the aliphatic hydrogen-abstraction 
in cytochromes P450; P. Rydberg, et al., J. Chem. Theor. Comput, 2007.

HES   INT     0
CORRECT  OMIT DU   BEG
   0.50000
   1   DUMM  DU    M      0   0   0       0.0000      0.0000      0.0000      0.000000
   2   DUMM  DU    M      1   0   0       1.4490      0.0000      0.0000      0.000000
   3   DUMM  DU    M      2   1   0       1.5220    111.1000      0.0000      0.000000
  19   FE    FE    M     16  14  12       2.1000    124.0000    180.0000     -0.158237
  20   NA    NP    S     19  16  14       2.0800     98.0000     90.0000      0.152116
  21   C1A   Cc    S     20  19  16       1.3800    125.4000     90.0000      0.238449
  22   C2A   Cb    B     21  20  19       1.4100    109.0000    180.0000     -0.232754
  23   CAA   c3    3     22  21  20       1.5100    124.0000    180.0000     -0.629069
  24   HAA1  hc    E     23  22  21       1.0900    109.5000     60.0000      0.277025
  25   HAA2  hc    E     23  22  21       1.0900    109.5000    300.0000      0.277025
  26   CBA   c3    3     23  22  21       1.5400    111.0000    180.0000     -0.333974
  27   HBA1  hc    E     26  23  22       1.0900    109.5000     60.0000      0.162019
  28   HBA2  hc    E     26  23  22       1.0900    109.5000    300.0000      0.162019
  29   CGA   c     B     26  23  22       1.5270    109.4000    180.0000      1.088328
  30   O1A   o     E     29  26  23       1.2600    117.2000     90.0000     -0.755717
  31   O2A   o     E     29  26  23       1.2600    117.2000    270.0000     -0.755717
  32   C3A   Cb    B     22  21  20       1.4100    107.0000      0.0000     -0.190936
  33   CMA   c3    3     32  22  21       1.5100    125.0000    180.0000     -0.233673
  34   HMA1  hc    E     33  32  22       1.0900    109.5000     60.0000      0.092911
  35   HMA2  hc    E     33  32  22       1.0900    109.5000    180.0000      0.092911
  36   HMA3  hc    E     33  32  22       1.0900    109.5000    300.0000      0.092911
  37   C4A   Cc    S     32  22  21       1.4100    107.0000      0.0000      0.136015
  38   CHB   Cd    B     37  32  22       1.3700    127.0000    180.0000     -0.283571
  39   HB    HA    E     38  37  32       1.0800    120.0000      0.0000      0.148888
  40   C1B   Cc    B     38  37  32       1.3700    127.0000    180.0000      0.147768
  41   NB    NO    E     40  38  37       1.3800    124.0000      0.0000     -0.095150
  42   C2B   Cb    B     40  38  37       1.4100    127.0000    180.0000     -0.144803
  43   CMB   c3    3     42  40  38       1.5100    125.0000      0.0000     -0.235651
  44   HMB1  hc    E     43  42  40       1.0900    109.5000     60.0000      0.093601
  45   HMB2  hc    E     43  42  40       1.0900    109.5000    180.0000      0.093601
  46   HMB3  hc    E     43  42  40       1.0900    109.5000    300.0000      0.093601
  47   C3B   Cb    B     42  40  38       1.4100    107.0000    180.0000     -0.266283
  48   CAB   cf    B     47  42  40       1.5100    126.0000    180.0000     -0.042694
  49   HAB   ha    E     48  47  42       1.0800    120.0000      0.0000      0.080372
  50   CBB   c2    B     48  47  42       1.3300    120.0000    180.0000     -0.174415
  51   HBB1  ha    E     50  48  47       1.0800    120.0000      0.0000      0.059614
  52   HBB2  ha    E     50  48  47       1.0800    120.0000    180.0000      0.059614
  53   C4B   Cc    S     47  42  40       1.4100    107.0000      0.0000      0.313680
  54   CHC   Cd    B     53  47  42       1.3700    127.0000    180.0000     -0.337061
  55   HC    HA    E     54  53  47       1.0800    120.0000      0.0000      0.144143
  56   C1C   Cc    B     54  53  47       1.3700    130.0000    180.0000      0.249094
  57   NC    NP    E     56  54  53       1.3800    124.0000      0.0000     -0.038372
  58   C2C   Cb    B     56  54  53       1.4100    127.0000    180.0000     -0.214761
  59   CMC   c3    3     58  56  54       1.5100    125.0000      0.0000     -0.253534
  60   HMC1  hc    E     59  58  56       1.0900    109.5000     60.0000      0.100347
  61   HMC2  hc    E     59  58  56       1.0900    109.5000    180.0000      0.100347
  62   HMC3  hc    E     59  58  56       1.0900    109.5000    300.0000      0.100347
  63   C3C   Cb    B     58  56  54       1.4100    107.0000    180.0000     -0.211987
  64   CAC   cf    B     63  58  56       1.5100    126.0000    180.0000      0.036650
  65   HAC   ha    E     64  63  58       1.0800    120.0000      0.0000      0.067469
  66   CBC   c2    B     64  63  58       1.3300    120.0000    180.0000     -0.234416
  67   HBC1  ha    E     66  64  63       1.0800    120.0000      0.0000      0.068238
  68   HBC2  ha    E     66  64  63       1.0800    120.0000    180.0000      0.068238
  69   C4C   Cc    S     63  58  56       1.4100    107.0000      0.0000      0.236376
  70   CHD   Cd    B     69  63  58       1.3700    127.0000    180.0000     -0.324695
  71   HD    HA    E     70  69  63       1.0800    120.0000      0.0000      0.139859
  72   C1D   Cc    B     70  69  63       1.3700    130.0000    180.0000      0.285203
  73   ND    NO    E     72  70  69       1.3800    124.0000      0.0000     -0.018192
  74   C2D   Cb    B     72  70  69       1.4100    127.0000    180.0000     -0.258314
  75   CMD   c3    3     74  72  70       1.5100    125.0000      0.0000     -0.115607
  76   HMD1  hc    E     75  74  72       1.0900    109.5000     60.0000      0.058633
  77   HMD2  hc    E     75  74  72       1.0900    109.5000    180.0000      0.058633
  78   HMD3  hc    E     75  74  72       1.0900    109.5000    300.0000      0.058633
  79   C3D   Cb    B     74  72  70       1.4100    107.0000    180.0000     -0.170110
  80   C4D   Cc    S     79  74  72       1.4100    107.0000      0.0000      0.161164
  81   CHA   Cd    S     80  79  74       1.3700    127.0000    180.0000     -0.321104
  82   HA    HA    E     81  80  79       1.0800    120.0000      0.0000      0.143117
  83   CAD   c3    3     79  74  72       1.5100    124.0000    180.0000     -0.087199
  84   HAD1  hc    E     83  79  74       1.0900    109.5000     60.0000      0.115702
  85   HAD2  hc    E     83  79  74       1.0900    109.5000    300.0000      0.115702
  86   CBD   c3    3     83  79  74       1.5400    111.0000    180.0000     -0.332992
  87   HBD1  hc    E     86  83  79       1.0900    109.5000     60.0000      0.099706
  88   HBD2  hc    E     86  83  79       1.0900    109.5000    300.0000      0.099706
  89   CGD   c     B     86  83  79       1.5300    109.4000    180.0000      0.952366
  90   O1D   o     E     89  86  83       1.2600    117.2000     90.0000     -0.762515
  91   O2D   o     E     89  86  83       1.2600    117.2000    270.0000     -0.762515

LOOP  EXPLICIT
 NA   C4A
 FE   NB
 FE   NC
 FE   ND
 NB   C4B
 NC   C4C
 ND   C4D
 C1A  CHA
 
DONE
oxygen bindning to iron

OFE   INT     0
CORRECT  OMIT DU   BEG
   0.50000
   1   DUMM  DU    M      0   0   0       0.0000      0.0000      0.0000      .0000
   2   DUMM  DU    M      1   0   0       1.4490      0.0000      0.0000      .0000
   3   DUMM  DU    M      2   1   0       1.5220    111.1000      0.0000      .0000
   4   OF    OQ    E      3   2   1       1.540     111.208     180.000     -0.529965

DONE