[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Fwd:




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 subunits

INCLUDE                             # 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       1



I 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