Ion binding to polyelectrolytes: Monte Carlo simulations versus classical mean field theories