Premium
Simplified methods for pK a and acid pH‐dependent stability estimation in proteins: Removing dielectric and counterion boundaries
Author(s) -
Warwicker James
Publication year - 1999
Publication title -
protein science
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 3.353
H-Index - 175
eISSN - 1469-896X
pISSN - 0961-8368
DOI - 10.1110/ps.8.2.418
Subject(s) - counterion , chemistry , dielectric , solvation , computational chemistry , chemical physics , ionic bonding , solvent , thermodynamics , ion , organic chemistry , quantum mechanics , physics
Much computational research aimed at understanding ionizable group interactions in proteins has focused on numerical solutions of the Poisson–Boltzmann (PB) equation, incorporating protein exclusion zones for solvent and counterions in a continuum model. Poor agreement with measured pK a s and pH‐dependent stabilities for a (protein, solvent) relative dielectric boundary of (4,80) has lead to the adoption of an intermediate (20,80) boundary. It is now shown that a simple Debye‐Hückel (DH) calculation, removing both the low dielectric and counterion exclusion regions associated with protein, is equally effective in general pK a calculations. However, a broad‐based discrepancy to measured pH‐dependent stabilities is maintained in the absence of ionizable group interactions in the unfolded state. A simple model is introduced for these interactions, with a significantly improved match to experiment that suggests a potential utility in predicting and analyzing the acid pH‐dependence of protein stability. The methods are applied to the relative pH‐dependent stabilities of the pore‐forming domains of colicins A and N. The results relate generally to the well‐known preponderance of surface ionizable groups with solvent‐mediated interactions. Although numerical PB solutions do not currently have a significant advantage for overall pK a estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large Δ pK a s of a subset of ionizable groups with the overall robustness of the DH model.