Ligand binding distributions in nucleic acids.

We illustrate a new method for the determination of the complete binding polynomial for nucleic acids based on experimental titration data with respect to ligand concentration. From the binding polynomial, one can then calculate the distribution function for the number of ligands bound at any ligand concentration. The method is based on the use of a finite set of moments of the binding distribution function, which are obtained from the titration curve. Using the maximum-entropy method, the moments are then used to construct good approximations to the binding distribution function. Given the distribution functions at different ligand concentrations, one can calculate all of the coefficients in the binding polynomial no matter how many binding sites a molecule has. Knowledge of the complete binding polynomial in turn yields the thermodynamics of binding. This method gives all of the information that can be obtained from binding isotherms without the assumption of any specific molecular model for the nature of the binding. Examples are given for the binding of Mn(2+) and Mg(2+) to t-RNA and for the binding of Mg(2+) and I(6) to poly-C using literature data.