This paper presents a Monte Carlo simulation (MCS) method for estimating the parameters that characterize ligand-receptor binding directly from experimentally derived binding isotherms. Binding parameters are estimated by incorporating an MCS algorithm for ligand binding to a two-dimensional receptor array into a nonlinear regression program. The MCS method was tested by analyzing experimental isotherms of avidin binding to biotinylated lipid in Langmuir-Blodgett (LB) monolayers. The MCS-derived cooperativity coefficients and intrinsic association constants for avidin-biotin binding to LB films are correlated strongly (R2 > 0.93) with the binding parameters determined from the same experimental data by a thermodynamic equilibrium binding model (Zhao et al. 1993. Langmuir. 9:3166-3173). This result shows MCS to be an accurate and potentially more versatile method for characterizing biomolecular interactions at surfaces.