[gmx-users] Top file
Justin A. Lemkul
jalemkul at vt.edu
Tue Jun 26 12:10:46 CEST 2012
On 6/26/12 1:52 AM, Shima Arasteh wrote:
> Dear gmx users,
> I have a problem about the topology output of pdb2gmx.
> There is a formyl residue as N-terminal in my .pdb file. I know that the pdb file is correct in agreement with RCSB .pdb files.
> force field I have to use is gmx ( I need to regenerate
> some results ), so I defined the formyl residue as a new residue in gmx
> force field. But when I get the output of pdb2gmx the topology shows
> that the N atom of the next residue has 3 H atoms ,I think It means that
> the formyl residue is not connected to the next residue ( here: Valine
> ). As a result the total charge derived by pdb2gmx is +2e and it is not
> correct regards to
> the protein structure. I don't know what is the problem? Is it related
> to the force field files?
Did you follow all of the instructions here?
It also sounds like you are not setting termini properly using pdb2gmx -ter.
Please provide your complete pdb2gmx command any any relevant screen output.
Justin A. Lemkul, Ph.D.
Department of Biochemistry
jalemkul[at]vt.edu | (540) 231-9080
More information about the gromacs.org_gmx-users