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

Re: HETATM_IO and strange log



Hi

yes, you need to use a . for each HETATM  you include in your ali. The
same for water. They treated separately. 
so a PDB file with 2 ion ligands and 4 waters may look like this:

>P1;1abc
JSDHFLSKJAGHSALKJ/..wwww*

these residues are not selected on the basis of making or not making
contact,
but simply following the order as they appear in the PDB file.


"/pipe1/databases/pdb/uncompressed_files/pdb1hci.ent" was our local
setup here, I do not see why it is popping up in your setup. 

Andras


Frank J A B wrote:
> 
> Dear All
> 
> I'm facing two problems with modeller.
> In the FAQ's question #16 there is an explanation on how to use hetero
> atoms in modelling adding the line HETATM_IO = on in the top script.
> 
> There is also a comment (not clear to me) about the
> use of Ws for water and .s for other atoms in the alignment that still
> remains unsolved to me.
> In my case I'm trying to use FAD during the modelling.
> I can't wonder how to edit the alignment to include this Ws and .s
> 
> I understand that the backslash is a sequence break, but I deduce (not
> sure) that you must to repeat for each residue a symbol like w, . or
> whatever in the same .pir file after the backslash. I didn't find any
> explanation along the manual regard this.
> 
> Check the example given.
> 
> INCLUDE
> SET ALIGNMENT_FILE = 'align1.ali', KNOWNS='templ1', SEQUENCE='targ1'
> SET HETATM_IO = on, WATER_IO = on
> CALL ROUTINE = 'model'
> 
> The alignment file:
> 
> >P1;templ1
> structureX:templ1:1::10::
> FAYVI/.wwww*
> 
> >P1;targ1
> sequence:targ1:1::8::
> -GWIV/.ww-w*
> 
> Need I to repeat the entire alignment (template/target) inside the same
> .pir file (but now using only signs like - , ., w , #) where
> the .s/Ws positions correpond to those aa that make contact with the
> HETATM in the template ?
> 
> The second question is about my log file.
> I'm habituated to use regularly Modeller since version 4 but suddenly the
> job was finished abruptly with this :
> 
>                          MODELLER 6v2, 17 Feb 2002
>      PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS
> 
> assgn___544E> Variable name not recognized:
> /pipe1/databases/pdb/uncompressed_files/pdb1hci.ent
>               Command line:  INCLUDE
> assgn___525E> Error in TOP variable assignment.
> recover____E> ERROR_STATUS >= STOP_ON_ERROR:  1   1
> 
> Is this some kind of bug or a user error ?
> I'm sure that the files are typed correctly (case-sensitive) and the
> folder where the job is started has all the files in it (the tops, ali
> and pdb files).
> I'm not using any file located in the path shown above :
> /pipe1/databases/pdb/uncompressed_files/pdb1hci.ent
> 
> Any help will be welcome
> 
> Thanks in advance
> 
> Frank Alarcon
> LNCC - Brasil
> www.lncc.br

-- 
    ,
Andras Fiser, PhD            # phone: (212) 327 7216
The Rockefeller University   # fax:   (212) 327 7540 
Box 270, 1230 York Avenue    # e-mail:
New York, NY 10021-6399, USA # http://salilab.org/~andras