I have a HETATM entry at the end of one of my protein chains. I have included the "env.io.hetatm = True" statement in my code. I have also included "." residues in my alignment file for each HETATM residue and formed a separate chain out of them using the "/" separator. I did this in both the template sequence and in the sequence to be modelled. However, I still get the following error:
-- read_te_290E> Number of residues in the alignment and pdb files are different: 1566 1560 For alignment entry: 1 1GWP.pdb x (mismatch at alignment position 1561)
Alignment KLFHIRTGNMDTYYWHRYT.. PDB KLFHIRTGNMDTYYWHRYT Match ******************** -- It appears that Modeller cannot read the HETATM entries in the PDB file. I wonder if anyone has a suggestion.
Thanks.