How optimal are the binding energetics of barnase and barstar?
The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.