Monte Carlo simulations are used to calculate the binding constant of two Ca2+ ions to the protein bovine calbindin D9k. The change in binding constant with respect to mutation of charged amino acids, presence of various electrolytes, protein concentration, solution pH, and competitive binding of monovalent ions is investigated. Each of these factors may have a large influence on the binding constant. The simulations are performed in a dielectric continuum model, the so-called primitive model of electrolyte theory, with a fixed protein structure and a uniform dielectric permittivity. The calculated binding constants are in excellent agreement with experimental data and describe changes in the binding constant over six orders of magnitude.