Pcetk

Pcetk calculates proton binding energetics within proteins to model protonation states, electrostatic energies, and titration behavior for studies of protein electrostatics.

Key Features:

• Object-oriented implementation: Provided as an object-oriented module of the pDynamo software library implemented in Python with performance-critical components in C and Cython.
• Protonation-state calculation: Computes protonation states for protonatable groups in proteins.
• Electrostatic energy calculation: Calculates electrostatic energies relevant to proton binding and protein electrostatics.
• Titration curve computation: Generates titration curves and related protonation-dependent properties.
• Poisson–Boltzmann solver interface: Interfaces with the extended-MEAD Poisson-Boltzmann solver for electrostatic modeling.
• Aqueous pKa estimation: Implements a computational method to estimate aqueous pKa (pKa aq) values with reported accuracy of ±0.5 pKa units or better.

Scientific Applications:

• pKa determination: Estimating pKa values of protonatable groups within proteins.
• Ligand protonation analysis: Predicting protonation states of protonatable ligands bound to proteins when experimental data are limited.
• Protein electrostatics studies: Modeling electrostatic contributions to protein structure, function, and energetics.

Methodology:

Pcetk is implemented as a pDynamo module in Python with C and Cython components, interfaces with the extended-MEAD Poisson-Boltzmann solver, and applies a computational method for estimating aqueous pKa (pKa aq) values with ~±0.5 pKa-unit accuracy.