Re: AMBER: GLUTATHIONE

From: Ilyas Yildirim <yildirim.pas.rochester.edu>
Date: Wed, 11 Jan 2006 10:38:56 -0500 (EST)

Dear Joachim,

You dont need the double_a info in your .pdb file while using antechamber.
Extract the necessary ATOM lines to a file. The .pdb file u need to create
must have the right format. See more on .pdb format at the following
website:

http://www.rcsb.org/pdb/file_formats/pdb/pdbguide2.2/guide2.2_frame.html

What I would do is to use the following command first, and play with the
.pdb file to make the right .pdb file format.

grep "ATOM " aps.dat > deneme.pdb

check the deneme.pdb file and be sure that the file has the ATOM info.s in
it. Then use the awk script as follows:

more test.pdb | awk '{print $1" "$2" "$3" "$4" "$5" "$7"
"$8" "$9}' > deneme_format.pdb

Check out deneme_format.pdb file and play with it a little bit to bring it
to the right .pdb format. Now, u have a deneme_format.pdb which needs to
be protonated. There is a program in AMBER which does that (never used
that, but must not be hard to use). U can use that protonated file in
antechamber, and create your .prepi file.

Or, if u want to see what u are doing, use the deneme_format.pdb
(which is in right .pdb format) file in xleap. Load it, and connect the
atoms by hand (it is not a big structure, so must not be hard to do).

Good luck,

On Wed, 11 Jan 2006, Joachim Reichelt wrote:

> Dear all,
>
> I want to model a structure with GLUTATHIONE.
> I extracted the GTT-Lines from pdb1gsa.ent with grep an tried
> to run antechamber. But all I get is:
>
> Cannot successfully assign bond type for this molecule, please :
> (1) double check the structure (the connectivity) and/or
> (2) adjust atom valence penalty parameters in APS.DAT, and/or
> (3) increase MAXVASTATE in define.h and recompile bondtype.C
>
> I increase MAXVASTATE, but it did not help, so I set debug to 1.
> Now I see:
> ...
> double_a C1 0 0 0 0
> double_a O11 0 0 0 0
> double_a O12 0 0 0 0
> double_a CB1 0 0 0 0
> double_a CG1 0 0 0 0
> double_a CD1 0 0 0 0
> double_a OE1 0 0 0 0
> double_a N2 0 0 0 0
> double_a CA2 0 1 -4 -1
> double_a C2 1 2 0 3
> double_a O2 0 0 0 0
> double_a CB2 0 0 0 0
> double_a SG2 0 0 0 0
> double_a N3 1 1 0 0
> double_a CA3 0 0 0 0
> double_a C3 0 0 0 0
> double_a O31 0 0 0 0
> double_a O32 0 0 0 0
> double 11 CA2 C2 6
> Cannot assign bond types for the current valence state
>
> The current valence state (35) is listed as the following
> ...
>
> More see in ant.log
>
>

-- 
  Ilyas Yildirim
  ---------------------------------------------------------------
  - Department of Chemisty       -				-
  - University of Rochester      -				-
  - Hutchison Hall, # B10        -				-
  - Rochester, NY 14627-0216     - Ph.:(585) 275 67 66 (Office)	-
  - http://www.pas.rochester.edu/~yildirim/			-
  ---------------------------------------------------------------
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Jan 12 2006 - 06:10:06 PST
Custom Search