[AMBER] ATP parameters: tleap error

From: Valentina Romano <valentina.romano.unibas.ch>
Date: Wed, 16 Jul 2014 16:41:55 +0000

Dear Amber users

I I want to run a MD of a kinase in complex with ATP.
To build ATP parameters I used files .frcmod and .prep found on http://www.pharmacy.manchester.ac.uk/bryce/amber/ and I loaded them in tleap.

I used the following script:

source leaprc.ff99SB
loadamberparams ATPphos.frcmod
loadamberprpe ATP.prep
ATP = loadpdb ATP.pdb
saveoff ATP ATP.lib
saveamberparm ATP ATP.prmtop ATP.inpcrd
quit

I got the following error:

Loading parameters: ./ATPphos.frcmod
Reading force field modification type file (frcmod)
Reading title:
# Modifications to the AMBER94 force field for polyphosphates
ERROR: syntax error

Loading PDB file: ./ATP.pdb
Unknown residue: ATP number: 0 type: Terminal/last
..relaxing end constraints to try for a dbase match
  -no luck
Creating new UNIT for residue: ATP sequence: 1279
Created a new atom named: PG within residue: .R<ATP 1279>
Created a new atom named: O1G within residue: .R<ATP 1279>
Created a new atom named: O2G within residue: .R<ATP 1279>
Created a new atom named: O3 within residue: .R<ATP 1279>
Created a new atom named: PB within residue: .R<ATP 1279>
Created a new atom named: O1B within residue: .R<ATP 1279>
Created a new atom named: O2B within residue: .R<ATP 1279>
Created a new atom named: O1 within residue: .R<ATP 1279>
Created a new atom named: PA within residue: .R<ATP 1279>
Created a new atom named: O1A within residue: .R<ATP 1279>
Created a new atom named: O2A within residue: .R<ATP 1279>
Created a new atom named: O3A within residue: .R<ATP 1279>
Created a new atom named: O5' within residue: .R<ATP 1279>
Created a new atom named: C5' within residue: .R<ATP 1279>
Created a new atom named: C4' within residue: .R<ATP 1279>
Created a new atom named: O4' within residue: .R<ATP 1279>
Created a new atom named: C3' within residue: .R<ATP 1279>
Created a new atom named: O3' within residue: .R<ATP 1279>
Created a new atom named: C2' within residue: .R<ATP 1279>
Created a new atom named: O2' within residue: .R<ATP 1279>
Created a new atom named: C1' within residue: .R<ATP 1279>
Created a new atom named: N9 within residue: .R<ATP 1279>
Created a new atom named: C8 within residue: .R<ATP 1279>
Created a new atom named: N7 within residue: .R<ATP 1279>
Created a new atom named: C5 within residue: .R<ATP 1279>
Created a new atom named: C6 within residue: .R<ATP 1279>
Created a new atom named: N6 within residue: .R<ATP 1279>
Created a new atom named: N1 within residue: .R<ATP 1279>
Created a new atom named: C2 within residue: .R<ATP 1279>
Created a new atom named: N3 within residue: .R<ATP 1279>
Created a new atom named: C4 within residue: .R<ATP 1279>
Created a new atom named: H5'1 within residue: .R<ATP 1279>
Created a new atom named: H5'2 within residue: .R<ATP 1279>
Created a new atom named: H4' within residue: .R<ATP 1279>
Created a new atom named: H3' within residue: .R<ATP 1279>
Created a new atom named: HO3' within residue: .R<ATP 1279>
Created a new atom named: H2' within residue: .R<ATP 1279>
Created a new atom named: HO2' within residue: .R<ATP 1279>
Created a new atom named: H1' within residue: .R<ATP 1279>
Created a new atom named: H8 within residue: .R<ATP 1279>
Created a new atom named: H61 within residue: .R<ATP 1279>
Created a new atom named: H62 within residue: .R<ATP 1279>
Created a new atom named: H2 within residue: .R<ATP 1279>
  total atoms in file: 43
  The file contained 43 atoms not in residue templates
Building topology.
Building atom parameters.
Checking Unit.
FATAL: Atom .R<ATP 1279>.A<PG 1> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O1G 2> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O2G 3> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O3 4> does not have a type.
FATAL: Atom .R<ATP 1279>.A<PB 5> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O1B 6> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O2B 7> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O1 8> does not have a type.
FATAL: Atom .R<ATP 1279>.A<PA 9> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O1A 10> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O2A 11> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O3A 12> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O5' 13> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C5' 14> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C4' 15> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O4' 16> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C3' 17> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O3' 18> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C2' 19> does not have a type.
FATAL: Atom .R<ATP 1279>.A<O2' 20> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C1' 21> does not have a type.
FATAL: Atom .R<ATP 1279>.A<N9 22> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C8 23> does not have a type.
FATAL: Atom .R<ATP 1279>.A<N7 24> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C5 25> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C6 26> does not have a type.
FATAL: Atom .R<ATP 1279>.A<N6 27> does not have a type.
FATAL: Atom .R<ATP 1279>.A<N1 28> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C2 29> does not have a type.
FATAL: Atom .R<ATP 1279>.A<N3 30> does not have a type.
FATAL: Atom .R<ATP 1279>.A<C4 31> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H5'1 32> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H5'2 33> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H4' 34> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H3' 35> does not have a type.
FATAL: Atom .R<ATP 1279>.A<HO3' 36> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H2' 37> does not have a type.
FATAL: Atom .R<ATP 1279>.A<HO2' 38> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H1' 39> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H8 40> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H61 41> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H62 42> does not have a type.
FATAL: Atom .R<ATP 1279>.A<H2 43> does not have a type.
Failed to generate parameters
Parameter file was not saved.
        Quit

It looked like tleap did not recognize the ATP atom types.
I changed all atom names in the pdb file of ATP according to the atom types used in the file .prep and I got the same error.

Any suggestion to understand and solve that problem?

Thank you
Valentina




~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Valentina Romano | PhD Student | Biozentrum, University of Basel & SIB Swiss Institute of Bioinformatics
Klingelbergstrasse 61 | CH-4056 Basel |

Phone: +41 61 267 15 80

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jul 16 2014 - 10:00:02 PDT
Custom Search