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.