Modelling Heme

Dear Modellers,

I am trying to model a protein containing a heme group. There is a residue type called heme in top_heav.lib, but I don`t know how to utilize it. I am aware of changing my alignment to include HETATMs but I haven`t succeeded yet.

Do I have to change the pdb-templates to include the heme residue? How can I make it clear that MODELLER 6v1 should use it?

Is it possible to receive the model with the heme group in it?

Thanks for any suggestions,

Christoph Nimptsch