Respected Sir
i am trying to add my model one zinc HETATM, but my efforts are failed , *iam getting *ATOM in my file pdb instead HETATM in the template structure Experimentally solved Residues Region *191-609* , protein Total Length: *419* in that Five residue are Missing Residue
Here i mentioned
Alignment file
>P1;2DW0 structureX:2DW0.pdb: 196 :A:+414 :A:MOL_ID 1; MOLECULE CATROCOLLASTATIN; CHAIN A, B; FRAGMENT RESIDUES 191-609; SYNONYM VAP2, VASCULAR APOPTOSIS-INDUCING PROTEIN-2, CATROCOLLASTATIN/VAP2B:MOL_ID 1; ORGANISM_SCIENTIFIC CROTALUS ATROX; ORGANISM_COMMON WESTERN DIAMONDBACK RATTLESNAKE; ORGANISM_TAXID 8730: 2.15: 0.17 PFRFVELVLVVDKAMVTKNNGDLDKIKTRMYEIVNTVNEIYRYMYIHVALVGLEIWSNEDKITVKPEAGYTLNAF GEWRKTDLLTRKKHDNAQLLTAIDLD-RVIGLAYVGSMCHPKRSTGIIQDYSEINLVVAVIMAHEMGHNLGINHD SGYCSCGDYACIMRPEISPEPSTFFSNCSYFECWDFIMNHNPECILNEPLGTDIISPPVCGNELLEVGEECDCGT PENCQNECCDAATCKLKSGSQCGHGDCCEQCKFSKSGTECRASMSECDPAEHCTGQSSECPADVFHKNGQPCLDN YGYCYNGNCPIMYHQCYDLFGADVYEAEDSCFERNQKGNYYGYCRKENGNKIPCAPEDVKCGRLYCKDNSPGQNN PCKMFYSNEDEHKGMVLPGTKCADGKVCSNGHCVDVATAY.*
>P1;TvLDH sequence:TvLDH: : : : ::: 0.00: 0.00 P-RYIELVIVADRAMVTKNNGDLTAIRTWVHQIVNDMTVMYRDLNIHITLAAVVIWNKRDLITVTSSAEDTLNLF GTWRETKYLKHRKHDNAQLLTGIILDDDTIGLAYVGGMCDPKKSVGIIEDHSTEHLLVAATMAHEMGHNLGMNHD ANQCNCGANGCVMSATLTEQTSYQFSDCSKDEYQNYLTDHNPQCILNQPLRTD-T-P-VSRNELLQ-NSAHPC-- --------CDPVTCQPK--RFCVSGPCCENCQFVRAGTVCNPAKGD-WMDDYCTGISSDCPRNPL---------- -----NGFR------------------------------------------------------------------ ----------------------------------------.*
and
Model Ligand python script which i used as followed
# Homology modeling with ligand transfer from the template from modeller import * # Load standard Modeller classes from modeller.automodel import * # Load the automodel class
log.verbose() # request verbose output env = environ() # create a new MODELLER environment to build this model in
# directories for input atom files env.io.atom_files_directory = ['.', '../atom_files']
# Read in HETATM records from template PDBs env.io.hetatm = True
a = automodel(env, alnfile = 'TvLDH-2DW0A.ali', # alignment filename knowns = '2DW0', # codes of the templates sequence = 'TvLDH') # code of the target a.starting_model= 1 # index of the first model a.ending_model = 5 # index of the last model # (determines how many models to calculate) a.make() # do the actual homology modeling
kindly give provide your valuable suggestion regarding error