Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein-inhibitor association: application to plasmepsin, cathepsin D and endothiapepsin-pepstatin complexes.

The protein-inhibitor binding energies of enzymes are often pH dependent, and binding induces either proton uptake or proton release. The proton uptake/release and the binding energy for three complexes with available experimental data were numerically studied: pepstatin-cathepsin D, pepstatin-plasmepsin II and pepstatin-endothiapepsin. Very good agreement with the experimental data was achieved when conformational changes were taken into account. The role of the desolvation energy and the conformational changes was revealed by modeling the complex, the separated molecules in the complex conformation and the free molecules. It was shown that the conformational changes induced by the complex formation are as important for the proton transfer as the loss of solvation energy caused by the burial of interface residues. The residues responsible for the proton transfer were identified and their contribution to the proton uptake/release calculated. These residues were found to be scattered along the whole protein rather than being localized only at the active site. In the case of cathepsin D, these residues were found to be highly conserved among the cathepsin D sequences of other species. It was shown that conformation and ionization changes induced by the complex formation are critical for the correct calculation of the binding energy. Taking into account the electrostatics and the van der Waals (vdW) energies within the Boltzmann distribution of energies and allowing ionization and conformation changes to occur makes the calculated binding energy more realistic and closer to the experimental value. The interplay between electrostatic and vdW forces makes the pH dependence of the binding energy smoother, because the vdW force acts in reaction to the changes of the electrostatic energy. It was found that a small fraction of the ionizable groups remain uncharged in both the free and complexed molecules. The sequence and structural position of these groups aligns well within the three proteases, suggesting that these may have specific role.