Re: [AMBER] bond missing in phosphate group after antechamber run

From: David A Case <case.biomaps.rutgers.edu>
Date: Thu, 8 Nov 2012 09:56:09 -0500

On Wed, Nov 07, 2012, Krisztina Feher wrote:
>
> I am trying to derive GAFF parameters for a modified DNA residue with
> antechamber. I capped the residue with H's and ran through antechamber
> and parmcheck without error. However the bond between the oxigen and
> phosphorus is missing in the mol2 file as well as the respective
> parameters in the .frcmod file. Could someone have a look at teh input
> pdb file for antechamber to figure out the reason for this?

Your pdb file is badly formatted. First, all the atoms should be in the same
residue, but atoms 47-50 are not. When you put everything in one residue,
make sure the atom names are all unique (you can't have two atoms named "O1",
for example.)

More importantly, atom names in pdb file must be in residues 13-16, and hence
can be at most 4 characters. You have atom names that are in column 17 (atoms
42-46). This leads antechamber to change your atom names, since it thinks
you have duplicates.

As an aside: a terminal phosphate group will only be doubly protonated under
extremely acidic conditions. Unless you are sure you really have what is in
your pdb file, you should probably create a residue library for the -1 or -2
charge state of this molecule.

...good luck...dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Nov 08 2012 - 07:00:02 PST
Custom Search