APBS

From PyMOLWiki
Revision as of 15:05, 26 July 2018 by Speleo3 (talk | contribs) (update lots of outdated information)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
APBS-generated electrostatic surface displayed in PyMOL

APBS, the Adaptive Poisson-Boltzmann Solver, is a freely available macromolecular electrostatics calculation program released under a BSD license. PyMOL can display the results of the calculations as an electrostatic potential molecular surface.

APBS Plugins for PyMOL

Both plugins make it possible to run APBS from within PyMOL, and then display the results as a color-coded electrostatic surface (units ) in the molecular display window (as with the image to the right).

Required Dependencies

The plugins require apbs and pdb2pqr.

apt-get install apbs pdb2pqr

After all components are installed, open the plugin and browse for apbs and pdb2pqr on the "Program Locations" tab (if they haven't been found automatically).

Troubleshooting

  • (this might be outdated information): If the B-factor is then APBS doesn't properly read in the PDB file and thus outputs garbage (or dies). To fix this, set all b factors to be less than 100.
    alter all, b=min(b,99.9)
    
    The problem stems from how to parse a PDB file. The PDB file originally was written when most people used FORTRAN programs, and so the file format was specified by columns, not by the more modern comma separated value format we tend to prefer today. For the latest on the PDB format see the new PDB format docs.
  • APBS has problems, sometimes, in reading atoms with alternate conformations. You can remove the alternate locations with a simple script removeAlt.
  • For pdb2pqr, RNA resdiue names must be RA, RC, RG, and RU.
    alter polymer & resn A+C+G+U, resn = "R" + resn
    
  • Incomplete Residues: Some truncated PDB files include a single backbone atom of the next residue, e.g. 2xwu chain B residue 954 atom N. pdb2pqr reports: Error encountered: Too few atoms present to reconstruct or cap residue LEU B 954 in structure!. The easiest solution is to remove that atom:
    remove /2xwu//B/954
    
PyMOL visualizing two maps at once

Using APBS

See APBS Electrostatics Plugin and APBS Tools2.1.