[gmx-users] output pdb reformatted
Justin Lemkul
jalemkul at vt.edu
Tue Jul 11 22:40:39 CEST 2017
On 7/11/17 4:01 PM, Alex wrote:
> Hi all,
>
> I have a tiny acetonitrile molecule here. The input pdb prior to EM is this:
>
> TITLE Gromacs Runs On Most of All Computer Systems
> REMARK THIS IS A SIMULATION BOX
> CRYST1 30.000 30.000 30.000 90.00 90.00 90.00 P 1
> MODEL 1
> ATOM 1 CT ACT 1 14.870 14.744 14.594 1.00 0.00
> C
> ATOM 2 HC ACT 1 15.390 13.760 14.594 1.00 0.00
> H
> ATOM 3 HC ACT 1 13.769 14.581 14.594 1.00 0.00
> H
> ATOM 4 HC ACT 1 15.162 15.317 13.685 1.00 0.00
> H
> ATOM 5 CZ ACT 1 15.253 15.501 15.794 1.00 0.00
> C
> ATOM 6 NZ ACT 1 15.555 16.097 16.740 1.00 0.00
> N
> TER
> ENDMDL
>
> The output, however, is this:
>
> TITLE Great Red Owns Many ACres of Sand
> REMARK THIS IS A SIMULATION BOX
> CRYST1 30.000 30.000 30.000 90.00 90.00 90.00 P 1 1
> MODEL 1
> ATOM 1 1 CT 1 15.013 15.013 14.819 1.00 0.00
> C
> ATOM 2 1 HC 1 14.952 13.831 14.385 1.00 0.00
> H
> ATOM 3 1 HC 1 13.827 14.934 14.400 1.00 0.00
> H
> ATOM 4 1 HC 1 15.140 15.157 13.573 1.00 0.00
> H
> ATOM 5 2 CZ 1 15.821 15.819 16.308 1.00 0.00
> C
> ATOM 6 2 NZ 1 16.401 16.397 17.446 1.00 0.00
> N
> TER
> ENDMDL
>
> Why? I don't think I've seen this happen before, and I'd like to keep the
> initial format. Any suggestions? My itp is for an ACT residue, which is
> included in the simulation topology input.
>
mdrun gets the residue and atom names from the topology, so a dump of the .tpr
should show something useful. It seems like something in the topology has
gotten broken - atom names becoming numbers, residue name becoming the atom
names. Bizarre. grompp should have complained loudly if there were mismatches
between coordinates and topology.
-Justin
--
==================================================
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul at outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul
==================================================
More information about the gromacs.org_gmx-users
mailing list