PotSim2: Simple package to segment and compare protein potential grids
Using gridDataFormats to read the grids and MDAnalysis to make pdb atom selections
Installation
pip install potsim2
PIPSA like scores
Usage: Open two grids and calculateThe protein electrostatic potential grids can generated with APBS
from potsim2 import PotGrid
# read grids
grid1 = PotGrid('A__02_01.pdb', 'A__02_01.pkl')
grid2 = PotGrid('A__11_01.pdb', 'A__11_01.pkl')
# calculate skin for grid1
skin_mask1 = grid1.get_skin_mask()
grid1.apply_mask(skin_mask1)
# calculate skin for grid2
skin_mask2 = grid2.get_skin_mask()
grid2.apply_mask(skin_mask2)
# calc the Hodgkin similarity index and PIPSA like distance
hsi, dis = grid1.score(grid2)
# export the skins in uncompressed OpenDX format to be visualized in PyMol/ChimeraX
grid1.export('A__02_01.dx')
grid2.export('A__11_01.dx')