Experimental data shows that the binding of human prolactin (hPRL) to human prolactin receptor (hPRLr-ECD) is strongly pH-dependent, while the binding of the same receptor to human growth hormone (hGH) is pH-independent. Here we carry in silico analysis of the molecular effects causing such a difference and reveal the role of individual amino acids. It is shown that the computational modeling correctly predicts experimentally determined pKa’s of histidine residues in an unbound state in the majority of the cases and the pH-dependence of the binding free energy. Structural analysis carried in conjunction with calculated pH-dependence of the binding revealed that the main reason for pH-dependence of the binding of hPRL-hPRLr-ECD is a number of salt-bridges across the interface of the complex, while no salt-bridges are formed in the hGH-hPRlr-ECD. Specifically, most of the salt-bridges involve histidine residues and this is the reason for the pH-dependence across a physiological range of pH. The analysis not only revealed the molecular mechanism of the pH-dependence of the hPRL-hPRLr-ECD, but also provided critical insight into the underlying physic-chemical mechanism.

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 pKas 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 pKa 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 pKa estimations, development based on consideration of microscopic solvation energetics in tandem with the continuum model could combine the large ΔpKas of a subset of ionizable groups with the overall robustness of the DH model.