[gmx-users] Handling of charge field in PDB files by x2top

[David van der Spoel]

On Wed, 2005-11-16 at 11:38 +0100, Steven Kirk wrote:
> Hello,
> 
> Can anyone inform me of a way to get atom charge information into a 
> topology produced by x2top from a PDB-format input file ? My PDB file 
> contains ATOM records with fields laid out as defined in the PDB 
> specification, containing the charge in characters 79 and 80 in each 
> ATOM record in the form '2+','1-' etc.
Not without programming, but in x2top it shouldn't be too hard. The pdb
fields are read as B-factor / Occupancy. You want to write a normal
float number (e.g. -2, 1) otherwise reading won't work.

> 
> Failing this, does anyone know any other (automated) way to get charge 
> information into a topology, possibly using another input format ? For 
> example, I know that some of the Gromacs programs can GENERATE 
> 'PDB-like' (e.g. grasp, .pqr) files containing charge information hidden 
> in B-factor, etc. fields for use in other programs, but how can I READ 
> such files into Gromacs and process the charge information correctly ?
tpbconv can selectively set charges to zero (for a subset). It wouldn't
be hard to modify it such that it sets charges to something else.


> 
> Many thanks in advance,
> Steve Kirk
-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,          75124 Uppsala, Sweden
phone:  46 18 471 4205          fax: 46 18 511 755
spoel at xray.bmc.uu.se    spoel at gromacs.org   http://xray.bmc.uu.se/~spoel
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++