Stepwise constants can be used to describe competitive proton and metal binding to macromolecules with a large number of sites. With the aim of accessing information on the microscopic binding model, we report an expression that connects the stepwise constants to the site-specific metal constants. This expression holds for a very general complexation model including heterogeneity, interactions, and chelate complexation. Assuming bidentate binding of the Cd ions to adjacent carboxylate groups in poly(acrylic acid), stepwise and intrinsic stability constants for proton and cadmium binding were estimated from the experimental data. Intrinsic values were split into specific and electrostatic contributions (by means of the Poisson - Boltzmann equation under cylindrical geometry). Free of the electrostatic contribution, the remaining Cd binding energy showed almost no dependence on the coverage and ionic strength, and the corresponding average values allowed for a reasonable reproduction of raw binding data. Small systematic discrepancies from the homogeneous behavior are critically discussed.