[gmx-users] system has non zero net charge

Yogesh Sharma yogesh.rma13 at gmail.com
Wed Jan 8 12:01:55 CET 2020


Greetings,
I am performing membrane protein simulation in the presence of a
biomolecule. topology for the molecule was downloaded from ATB server.
After simulation run, I found unexpected behviour of the added
biomolecule.  This ligand named UINL was showing affinity to ASP and GLU
residues. ligand didnt seperate through out simulation once bound, which is
highly unexpected at neutral pH.

I was getting these two notes during solvation.

Generated 1062 of the 3486 non-bonded parameter combinations
Excluding 3 bonded neighbours molecule type 'Protein'
turning H bonds into constraints...
Excluding 3 bonded neighbours molecule type 'POPC'
turning H bonds into constraints...
Excluding 3 bonded neighbours molecule type 'U1NL'
turning H bonds into constraints...
Excluding 2 bonded neighbours molecule type 'SOL'
turning H bonds into constraints...

NOTE 1 [file topol2.top, line 18875]:
  System has non-zero total charge: -1.000002
  Total charge should normally be an integer.

Is it possible this -1 charge is producing artifact. I have addded 20
molecules of same ligand in my system and 6 of them are binding to negative
charged residues.
I checked total charge on ligand.itp file, its zero. then i checked line
18875 in topol2.top its SOL which is highly unlikely to have any charge.

U1NL     3
[ atoms ]
;  nr  type  resnr  resid  atom  cgnr  charge    mass
    1  HS14    1    U1NL     H4     1    0.480   1.0080
    2 OEOpt    1    U1NL     O4    2   -0.814  15.9994
    3    SI         1    U1NL     Si     3    1.336   28.0800
    4 OEOpt    1    U1NL     O1    4   -0.814  15.9994
    5  HS14     1    U1NL     H1     5    0.480   1.0080
    6 OEOpt    1    U1NL     O2    6   -0.814  15.9994
    7  HS14     1    U1NL     H3     7    0.480   1.0080
    8 OEOpt    1    U1NL     O3    8   -0.814  15.9994
    9  HS14    1    U1NL     H2     9    0.480   1.0080
; total charge of the molecule:  -0.000

How can i check what molecule is throwing this negative -1 charge to my
system?
Thankyou in advance for your valuable time.


More information about the gromacs.org_gmx-users mailing list