We discuss the rationale for using the reverse Boltzmann methodology to convert atomatom distance distributions from a knowledge base of protein-ligand complexes into energy-like functions. We also generate an updated version of the BLEEP statistical potential, using a dataset of 196 complexes. This performs similarly to the existing BLEEP. An algorithm is implemented to allow the automatic calculation of bond orders, and hence of the appropriate numbers of hydrogen atoms present. An attempt is made to generate a potential specific to strongly bound complexes; however, we find no evidence that this improves the prediction of binding affinities. We also discuss the range of binding energies available as a function of ligand molecular weight and derive some simple functions describing this behaviour.