Ligand binding systems at equilibrium: specificity, heterogeneity, cross-reactivity, and site-site interactions

The simplest ligand binding model that can account for systems which consist of heterogeneous receptors that show site-site interactions and can react with several types of ligands was presented. This model was based on the grand canonical partition function whose properties and potentialities for obtaining binding functions are briefly illustrated. A computer simulation of this model was carried out in order to examine the contributions of site-site interactions on the shape of binding isotherms and on specificity curves. The shape of the binding isotherms was shown to be independent regardless of whether the path of binding was considered. However, the shape of the specificity curves was strongly dependent on site-site interactions and on the particular sequence of ligand-receptor configurations formed during the ligand binding process, leading to the conclusion that site-site interactions may influence the specificity of a system more than even very strong cross-reactions. A method based on the Ising model of statistical mechanics was also described and was used for obtaining binding functions applicable to infinite lattices which show site-site interactions and competitive binding. The results presented here point out possible errors that can arise if standard statistical methods are used to fit binding data in order to prove a given binding model.