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.

Human growth hormone (hGH) binds to its receptor (hGHr) in a three-body interaction: one molecule of the hormone and two identical monomers of the receptor form a trimer. Curiously, the hormone-receptor interactions in the trimer are not equivalent and the formation of the complex occurs in a specific kinetic order (Cunningham BC, Ultsch M, De Vos AM, Mulkerrin MG, Clauser KR, Wells JA, 1991, Science 254:821–825). In this paper, we model the recognition of hGH to the hGHr using shape complementarity of the three-dimensional structures and macromolecular docking to explore possible binding modes between the receptor and hormone. The method, reported previously (Hendrix DK, Kuntz ID, 1998, Pacific symposium on biocomputing 1998, pp 1234–1244), is based upon matching complementary-shaped strategic sites on the molecular surface. We modify the procedure to examine three-body systems. We find that the order of binding seen experimentally is also essential to our model. We explore the use of mutational data available for hGH to guide our model. In addition to docking hGH to the hGHr, we further test our methodology by successfully reproducing 16 macromolecular complexes from X-ray crystal structures, including enzyme-inhibitor, antibody-antigen, protein dimer, and protein-DNA complexes.