Hi all,

I am trying to model an ion, water and ligand from a template and below is my alignment file.


>P1;SAM_docked
structureX:SAM_docked.pdb:3:A:+1048:B:::-1.00:-1.00
-S-QVEHPAGGYKKLFETVEELSSPLTAHVTGRIPLWLTGSLLRCGPGLFEVGSEPFYHLFDGQALLHKFDFKEG
HVTYHRRFIRTDAYVRAMTEKRIVITEFGTCAFPDPCKNIFSRFFSYFRGVEVTDNALVNIYPVGEDYYACTETN
FITKVNPETLETIKQVDLCNYVSVNGATAHPHIENDGTVYNIGNCF-----IAYNIVKIPPLQADKEDPISKSEI
VVQFPCSDRFKPSYVHSFGLTPNYIVFVETPVKINLFKFLSS----GANYMDCFESNETMGVWLHIADKKRKKYI
NNKYRTSPFNLFHHINTYEDHEFLIVDLCCWKGFEFVYNYLYLANLRENWEEVKKNARKAPQPEVRRYVLPLNID
KADTGKNLVTLPNTTATAILCSDETIWLEPEVLFSGPRQAFEFPQINYQKYGGKPYTYAYGLGLNHFVPDRLCKL
NVKTKETWVWQEPDSYPSEPIFVSHPDALEEDDGVVLSVVVSPGAGQKPAYLLILNAKDLSEVARAEVEINIPVT
FHGLFKKS.w./S--QVEHPAGGYKKLFETVEELSSPLTAHVTGRIPLWLTGSLLRCGPGLFEVGSEPFYHLFDG
QALLHKFDFKEGHVTYHRRFIRTDAYVRAMTEKRIVITEFGTCAFPDPCKNIFSRFFSYFRGVEVTDNALVNIYP
VGEDYYACTETNFITKVNPETLETIKQVDLCNYVSVNGATAHPHIENDGTVYNIGNCF-----IAYNIVKIPPLQ
ADKEDPISKSEIVVQFPCSDRFKPSYVHSFGLTPNYIVFVETPVKINLFKFLG------ANYMDCFESNETMGVW
LHIADKKRKKYINNKYRTSPFNLFHHINTYEDHEFLIVDLCCWKGFEFVYNYLYLANLRENWEEVKKNARKAPQP
EVRRYVLPLNIDKADTGKNLVTLPNTTATAILCSDETIWLEPEVLFSGPRQAFEFPQINYQKYGGKPYTYAYGLG
LNHFVPDRLCKLNVKTKETWVWQEPDSYPSEPIFVSHPDALEEDDGVVLSVVVSPGAGQKPAYLLILNAKDLSEV
ARAEVEINIPVTFHGLFKKS.w.*

>P1;EGR
sequence:EGR:: :: ::: 0.00: 0.00
MSIQVEHPAGGYKKLFETVEELSSPLTAHVTGRIPLWLTGSLLRCGPGLFEVGSEPFYHLFDGQALLHKFDFKEG
HVTYHRRFIRTDAYVRAMTEKRIVITEFGTCAFPDPCKNIFSRFFSYFRGVEVTDNALVNVYPVGEDYYACTETN
FITKINPETLETIKQVDLCNYVSVNGATAHPHIENDGTVYNIGNCFGKNFSIAYNIVKIPPLQADKEDPISKSEI
VVQFPCSDRFKPSYVHSFGLTPNYIVFVETPVKINLFKFLSSWSLWGANYMDCFESNETMGVWLHIADKKRKKYL
NNKYRTSPFNLFHHINTYEDNGFLIVDLCCWKGFEFVYNYLYLANLRENWEEVKKNARKAPQPEVRRYVLPLNID
KADTGKNLVTLPNTTATAILCSDETIWLEPEVLFSGPRQAFEFPQINYQKYCGKPYTYAYGLGLNHFVPDRLCKL
NVKTKETWVWQEPDSYPSEPIFVSHPDALEEDDGVVLSVVVSPGAGQKPAYLLILNAKDLSEVARAEVEINIPVT
FHGLFKKS.w./MSIQVEHPAGGYKKLFETVEELSSPLTAHVTGRIPLWLTGSLLRCGPGLFEVGSEPFYHLFDG
QALLHKFDFKEGHVTYHRRFIRTDAYVRAMTEKRIVITEFGTCAFPDPCKNIFSRFFSYFRGVEVTDNALVNVYP
VGEDYYACTETNFITKINPETLETIKQVDLCNYVSVNGATAHPHIENDGTVYNIGNCFGKNFSIAYNIVKIPPLQ
ADKEDPISKSEIVVQFPCSDRFKPSYVHSFGLTPNYIVFVETPVKINLFKFLSSWSLWGANYMDCFESNETMGVW
LHIADKKRKKYLNNKYRTSPFNLFHHINTYEDNGFLIVDLCCWKGFEFVYNYLYLANLRENWEEVKKNARKAPQP
EVRRYVLPLNIDKADTGKNLVTLPNTTATAILCSDETIWLEPEVLFSGPRQAFEFPQINYQKYCGKPYTYAYGLG
LNHFVPDRLCKLNVKTKETWVWQEPDSYPSEPIFVSHPDALEEDDGVVLSVVVSPGAGQKPAYLLILNAKDLSEV
ARAEVEINIPVTFHGLFKKS.w.*


and I tried a model using the below specified script.

from modeller import *
from modeller.automodel import *
#from modeller import soap_protein_od

env = Environ()
env.io.hetatm = True
env.io.water = True
class MyModel(AutoModel):
    def special_restraints(self, aln):
        rsr = self.restraints
        at = self.atoms
#       Add some restraints from a file:
        rsr.add(secondary_structure.Alpha(self.residue_range('111:A', '123:A')))
        rsr.add(secondary_structure.Alpha(self.residue_range('646:B', '658:B')))

a = MyModel(env, alnfile='alignment.ali',
              knowns='SAM_docked', sequence='EGR',
              assess_methods=(assess.DOPE,
                              #soap_protein_od.Scorer(),
                              assess.GA341))


a.starting_model = 1
a.ending_model = 1
a.make()

The modelled structure generated without any error containing ion and water. But only two atoms were observed in the modelled structure but the ligand in the template contains 38 atoms. The atoms in ligand were not retained in the model from the template.

Please provide me with suggestions to solve this error?

Samdani