[gmx-users] NVT Equilibration error

Justin Lemkul jalemkul at vt.edu
Tue May 5 18:26:32 CEST 2015



On 5/5/15 12:18 PM, Agnivo Gosai wrote:
> Hi,
>
> This is the warning message I received while preparing the topology. It
> occurred for one chain of the protein.
>
> Making bonds...
> Warning: Long Bond (2429-2431 = 0.80585 nm)
> Number of bonds was 4137, now 4136
> Generating angles, dihedrals and pairs...
> Before cleaning: 10853 pairs
> Before cleaning: 10993 dihedrals
> Keeping all generated dihedrals
> Making cmap torsions...There are 10993 dihedrals,  841 impropers, 7477
> angles
>            10769 pairs,     4136 bonds and     0 virtual sites
> Total mass 29166.760 a.m.u.
> Total charge 5.000 e
>
> I checked the two atoms mentioned in the warning part and found that they
> are in two residues of the protein chain in question. I am not sure if in
> physiological pH they bond with each other or not.
>
> This is the raw data from the pdb file :
>
> ATOM   2429  OG1 THR H 229      24.393  41.044  35.727  1.00  2.93
>    O
> ATOM   2430  CG2 THR H 229      23.066  40.670  37.642  1.00  6.21
>    C
> ATOM   2431  N      HIS H 230      25.280  38.177  39.029  1.00 10.19
>      N
>
> Should I stop this bond making by somehow ? Below the warning it is written
> : Number of bonds was 4137, now 4136 , so did pdb2gmx remove this bond ?
>

This comes from merging termini and .rtp entries; it's not relevant to your problem.

An 8-Angstrom separation is massive.  Either there are missing residues 
(checking for obvious MISSING lines in the .pdb file) or the assigned geometry 
is really bad.  There's probably some residue missing.  Never overlook long and 
short bond warnings.  They usually indicate trouble.

-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