DESIGN: computerized optimization of experimental design for estimating Kd and Bmax in ligand binding experiments. I. Homologous and heterologous binding to one or two classes of sites.

We have developed a versatile computer program for optimization of ligand binding experiments (e.g., radioreceptor assay system for hormones, drugs, etc.). This optimization algorithm is based on an overall measure of precision of the parameter estimates (D-optimality). The program DESIGN uses an exact mathematical model of the equilibrium ligand binding system with up to two ligands binding to any number of classes of binding sites. The program produces a minimal list of the optimal ligand concentrations for use in the binding experiment. This potentially reduces the time and cost necessary to perform a binding experiment. The program allows comparison of any proposed experimental design with the D-optimal design or with assay protocols in current use. The level of nonspecific binding is regarded as an unknown parameter of the system, along with the affinity constant (Kd) and binding capacity (Bmax). Selected parameters can be fixed at constant values and thereby excluded from the optimization algorithm. Emphasis may be placed on improving the precision of a single parameter or on improving the precision of all the parameters simultaneously. We present optimal designs for several of the more commonly used assay protocols (saturation binding with a single labeled ligand, competition or displacement curve, one or two classes of binding sites), and evaluate the robustness of these designs to changes in parameter values of the underlying models. We also derive the theoretical D-optimal design for the saturation binding experiment with a homogeneous receptor class.