I am trying to add a ligand as new residue type. My ligand is acetylcholine. I found its topology in CHARMM /toppar/top_all27_port_lipid.rtf as ACHO. This is how I did:
I copy the part from CHARMM file for ACHO to the end of $MODLIB/top_heav.lib. Then add a line at the end of $MODLIB/restyp.lib 142 | ACH | # | ACHO | acetylcholine
I use # in the alignment file
I have the following in my top file: SET TOPLIB = '$(LIB)/top_heav.lib' SET PARLIB = '$(LIB)/par.lib' SET HETATM_IO = 'ON'
However, mod7v7 still indicate that I have unknown residue type. Where I was wrong?