Zinc endopeptidase thermolysin can be inhibited by a series of phosphorus-containing peptide analogues, CbzGly-$(PO2)-X-Leu-Y-R (ZGP(X)L(Y)R), where X = NH, 0, or CH2; Y = NH or 0; and R = Leu, Ala, Gly, Phe, H , or CH3. The affinity correlation as well as an X-ray crystallography study suggest that these inhibitors bind to thermolysin in an identical mode. In this work, we calculate the electrostatic binding free energies for a series of 13 phosphorus-containing inhibitors with modifications at X, Y, and R moieties using finite difference solution to the Poisson-Boltzmann equation. A method has been developed to include the solvation entropy changes due to binding different ligands to a macromolecule. We demonstrate that the electrostatic energy and empirically derived solvation entropy can account for most of the binding energy differences in this series. By analyzing the binding contribution from individual residues, we show that the energy of a hydrogen bond is not confined to the donor and acceptor. In particular, the positive charges on Zn and Arg 203, which are not the acceptors, contribute significantly to the hydrogen bonds between two amides of ZGPLL and the thermolysin. Morgan et al., 1991). It has been shown that modifications on X or Y yield uniform binding constant changes for those inhibitors with identical R. In particular, the average binding energy change between phosphonamidates and phosphonates is 4.0 * 0.1 kcal/mol. These correlations, as well as an X-ray crystallography study Tronrud et al., 1987;Matthews, 1988), suggest that these inhibitors bind to thermolysin in an identical mode.Previous theoretical studies focused on the variation of the X moiety of these inhibitors. Bash et al. (1987) and Merz and Kollman (1989) carried out free energy perturbation (FEP) simulations for ZGP(X)LL with thermolysin. The good agreement between the calculations and experiments indicates that the solvation of the inhibitors and the interaction between the X moiety and a hydrogen bond (H-bond) acceptor of the enzyme are two primary determinants of the binding affinity difference. Through a qualitative analysis, Bartlett and Marlowe (1987) suggested that the intrinsic H-bond energy for phosphonamidates Reprint requests to: Jian Shen at his present address: Marion Merre11 Dow, lnc., 2110 East Galbraith Road, P.O. Box 156300, Cincinnati, is significantly higher than the observed binding energy change of 4 kcal/mol from the PI' NH to 0 substitution. However, the calculated solvation energies as well as the interpretation were challenged by partition coefficient and pK, measurements. Grohelny et al. (1989) proposed that the basicities of the inhibitors and a much smaller H-bond energy of 1.5 kcal/mol are responsible for the loss of binding affinity in going from phosphonamidates to phosphonates. Although Morgan et al. (1991) showed that the basicities of another series of phosphorus-containing inhibitors are not critical to the binding, a final conclusion has not been reached due to the lack of relevant ...