[gmx-users] g_sas values of zero and Warning

Nash, Anthony a.nash at ucl.ac.uk
Tue Jan 13 07:43:06 CET 2015


HI All,

I'm trying to get a measure of the solvent accessible surface area of a protein's catalytic site. It is unknown precisely how the substrate actually fits in the site given that in the crystal structure the site it too enclosed for the bulky substrate. I am using a progressive measure of g_sas -probe to get an idea of how accessible the catalytic site is during the simulation.

So I've tried six different -probe sizes, 0.1 (nm) 0.8, default , 1.6, 2.4, 3.2. Only in the first two instances do I actually get any read out. It appears as though I am unable to access the active site residues as the probe radius is near to a vdw radius. I could asume that the catalytic site is very enclosed, which could physiologically make sense. Or I have done something wrong, after all I am getting the warning:


WARNING: could not find a Van der Waals radius for 2 atoms
205 out of 205 atoms were classified as hydrophobic

Generated values per probe radius below:


-probe 0.1 (nm)
         0     8.30466           0     8.30466
       100     8.04708           0     8.04708
       200     7.82871           0     7.82871
.....

DEFAULT: -probe       real   0.14    Radius of the solvent probe (nm)
         0     4.30913           0     4.30913
       100     4.24779           0     4.24779
       200     3.81957           0     3.81957

-probe 0.8 (nm)
         0           0           0           0
       100           0           0           0
       200           0           0           0
......

-probe 1.6 (nm)
         0           0           0           0
       100           0           0           0
       200           0           0           0
......

-probe 2.4
         0           0           0           0
       100           0           0           0
       200           0           0           0
....

-probe 3.2
         0           0           0           0
       100           0           0           0
       200           0           0           0

Many thanks
Anthony


More information about the gromacs.org_gmx-users mailing list