OEGet2dPSA

OEGet2dPSA(mol: Union[OEGraphMol,OEMol,OEQMol],
           atomPSA: OEFloatArray = None, SandP: bool = False) -> bool

Returns the topological polar surface area for a given molecule as described in the Polar Surface Area section. The SandP parameter controls whether sulfur and phosphorus should be counted towards the total surface area. See example in Figure: Example of depicting the atom contributions of the polar surface area.

Note

TPSA values are mildly sensitive to the protonation state of a molecule.

The atomPSA parameter can be used to retrieve the contribution of each atom to the total polar surface area as shown in Listing 1.

Listing 1: Example of retrieving individual atom contributions to PSA.

atomPSA = oechem.OEFloatArray(mol.GetMaxAtomIdx())
psa = oemolprop.OEGet2dPSA(mol, atomPSA)

print("PSA =", psa)
for atom in mol.GetAtoms():
    idx = atom.GetIdx()
    print(idx, atomPSA[idx])

Full listing.