Dear modellers,I've succesfully modelled a dimer with MODELLER 6. Now, I would like to put the substrates in my model and optimize the interactions, but I might be doing something wrong. I've initially placed the substrates (one per subunit) as they are located at the template, and I'm using the following alignment and top files:ALINGMENT: The template pdb corresponds to a dimer with two AMPPNP molecules (one per subunit). The target sequence has been written twice to align it with the two subunits of the template.>P1;ag_anp structureX:ag_anp:1:A:518:A MMNPLIIKLGGVLLDSEEALERLFSALVNYRESHQRPLVIVHGGGCVVDELMKGLNLPV--- KKKNGLRVTPADQIDIIT-GALAGTANKTLLAWAKKHQIAAVGLFLGDGDSVKVTQL DEELGHVGLAQPGSPKLINSLLENG----YLPVVSSIGVTDE------GQLMNV NADQAATALAATLGA-DLILLSDVSGILDGKGQRIAEMTAAKAEQLIEQGIITDGMIVK VNAALDAARTLGR-----PVDIASWRHAEQLPALFNGMPMGTRILA------. / MMNPLIIKLGGVLLDSEEALERLFSALVNYRESHQRPLVIVHGGGCVVDELMKGLNLPV--- KKKNGLRVTPADQIDIIT-GALAGTANKTLLAWAKKHQIAAVGLFLGDGDSVKVTQL DEELGHVGLAQPGSPKLINSLLENG----YLPVVSSIGVTDE------GQLMNV NADQAATALAATLGA-DLILLSDVSGILDGKGQRIAEMTAAKAEQLIEQGIITDGMIVK VNAALDAARTLGR-----PVDIASWRHAEQLPALFNGMPMGTRILA------.* >P1;akIII sequence:akIII:1::502 MSEIVVSKFGGTSVADFDAMNRSADIVLSDANVR---LVVLSASAGITNLLVALAEGLEPGE /-------SPALTDELVSH-GELMSTLLFVEILRER-DVQAQW-F----DVRKVMRT NDRFGRAE---PDIAALAELAALQLLPRLNEGLVITQGFIGSENKGRTTTLGRG GSDYTAALLAEALHASRVDIWTDVPGIY--TTDPRVV-SAAKRIDEIAFA-EAAEMATF GAKVLHPATLLPAVRSDIPVFVGSSKD----PRAG-----GTLVCNKTENPP. / MSEIVVSKFGGTSVADFDAMNRSADIVLSDANVR---LVVLSASAGITNLLVALAEGLEPGE /-------SPALTDELVSH-GELMSTLLFVEILRER-DVQAQW-F----DVRKVMRT NDRFGRAE---PDIAALAELAALQLLPRLNEGLVITQGFIGSENKGRTTTLGRG GSDYTAALLAEALHASRVDIWTDVPGIY--TTDPRVV-SAAKRIDEIAFA-EAAEMATF GAKVLHPATLLPAVRSDIPVFVGSSKD----PRAG-----GTLVCNKTENPP.*TOP FILE: an initial model named ter_anp.pdb is given as the starting model. A subroutine is included to restrain the same changes in both subunitsINCLUDE # Include the predefined TOP routines SET OUTPUT_CONTROL = 1 1 1 1 1 SET ALNFILE = 'dimer_anp.ali' # alignment filename SET KNOWNS = 'ag_anp' # codes of the templates SET SEQUENCE = 'akIII' # code of the target SET ATOM_FILES_DIRECTORY = '.' # directories for input atom files SET MODEL = 'ter_anp.pdb' , GENERATE_METHOD = 'read_xyz' SET STARTING_MODEL= 1 # index of the first model SET ENDING_MODEL = 1 # index of the last model# (determines how many models to calculate)SET HETATM_IO = on CALL ROUTINE = 'model' # do homology modelling SUBROUTINE ROUTINE = 'special_restraints' SET RES_TYPES = 'ALL' SET ATOM_TYPES = 'ALL' SET SELECTION_STATUS = 'INITIALIZE' SET SELECTION_SEARCH = 'SEGMENT' SET SYMMETRY_WEIGHT = 0.5 PICK_ATOMS PICK_ATOMS_SET = 2, SELECTION_SEGMENT = '1:''251:' PICK_ATOMS PICK_ATOMS_SET = 3, SELECTION_SEGMENT = '252:''502:' DEFINE_SYMMETRY ADD_SYMMETRY = on off RETURN END_SUBROUTINE The log file gives the following error: preppdf_203E> CHARMM atom type is out of range; Probably because GENERATE_TOPOLOGY was not called. Atom index, atom type: 0 recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1I would be very thankful if anyone could tell me what I am doing wrong and if there is a better strategy to do it.Many thanks in advance
Santiago