apbs psize.py

When using psize.py with my PQR file (psize.py myfile.pqr) it gives:

No ATOM entires in file!

However, when I use the Pymol plugin with the same PQR and the .mol to match, it works just fine.

And when APBS using the same PQR file (after the cglen & fglen have been set in the input file manually), it works fine as well.

It appears the error is because the PQR is for a small molecule (not a protein) and, as such, only has HETATM declarations. Nonetheless, I'm not sure how to proceed.

Also, do I correctly understand that the input file values, cglen and fglen are usually calculated by psize.py and that you do need these values in your input file prior to running APBS. In other words, APBS does not call psize.py during run time -- correct?

thanks!

Scott

Mar 22 '21 23:03 sbembenek18

Hi --

Much of the psize.py functionality has been integrated into pdb2pqr; we are planning to deprecate the psize.py tool provided with APBS. Have you tried PDB2PQR?

Thanks,

Nathan

Mar 26 '21 01:03 sobolevnrm

Will do -- thanks!

Mar 27 '21 06:03 sbembenek18

Unfortunately, I am unable to generate a PQR for the ligand. I'm using the command:

pdb2pqr30 --ff=PARSE --ligand=frag.mol2 frag.pdb out.pqr

I actually have a pqr file that works with APBS -- I'm really looking for a way to generate a (new) input for APBS and get parameter suggestion therein. However -- as noted -- psize isn't working with my current PQR, so I am not sure how to get the info in an automated fashion otherwise (other than going through the Pymol Plugin).

I've attached the files.

thanks for the help!

frag.mol2.txt frag.pdb.txt

Mar 27 '21 08:03 sbembenek18