Resolution of steroid binding heterogeneity by Fourier-derived affinity spectrum analysis (FASA)

A new mathematical method of analyzing radioreceptor assay data is presented. When there are many binding classes with different affinities, the probability-density function B(p) is described by the equation B(p) = (integral negative infinity to infinity) q(k)f(p-k)dk, where q(k) is the affinity spectrum (density of a particular binding class as a function of affinity) and f(p-k) is a probability function (probability that dissociation constants will fall between k and p-k, where p is the free ligand concentration). This equation is solved for q(k) and evaluated explicitly by Fourier transformation, namely, q(w) = b(w)/f(w), where w is frequency. Since division by f(w) can amplify and high frequency noise present in the experimental data, a Gaussian smoothing function is introduced thus: qs(w) = q(w)e(-w/W0)2, where W0 is a constant. This produces an affinity spectrum defined as a plot of the number of binding sites, qs(k), versus their respective dissociation constants, k. Using a FORTRAN computer program, we verify this algorithm using simulated data. We also apply the procedure to resolve heterogeneous populations of estrogen binders in human endometrium using [3H]estradiol as ligand. Two estrogen binder classes are revealed with dissociation constants approximately 2.5 natural logarithmic units apart. We identify one high-affinity (Kd = 0.18 nM)-low density (70 pM [or 72 fmol/mg protein]) subpopulation and one low affinity (Kd = 2.5 nM)-high density (101 pM [or 102 fmol/mg protein]) subpopulation of estradiol binders. The management of experimental error, sampling limitations, and nonspecific binding are discussed. This method directly transforms experimental data into an easily interpretable representation without mathematical modeling or statistical procedures.