Monte Carlo simulations are used to describe the charging behavior of metal oxide nanoparticles thus enabling a novel and original approach to predict nanoparticle reactivity and possible interactions with biological and environmental molecules. The charging behavior of spherical nanoparticles is investigated by adjusting the pH of the media and the influence of surface site distribution, density and dielectric constant as well as acid/base properties of the surface sites and values (difference between two successive deprotonation constants) is systematically studied using a grand canonical Monte Carlo method. A primitive Coulomb model is applied to describe the interaction energies between explicit discrete sites. Homogeneous, heterogeneous surfaces and patches with homogeneous and heterogeneous distributions are considered in order to reproduce possible site distributions of metal oxide nanoparticles. Two models are used. In the model (one deprotonation step) results indicate that the deprotonation process is controlled by inter-site distances which are defined by site distributions and densities. It is shown that the homogeneous surface is the most efficient site distribution to obtain high ionization degrees. In the model (two deprotonation steps), the value is found to control the surface charge properties with regards to pH changes. By considering the variation of the total NPs surface charge as a function of pH our results help in the distinction between the zero charge and isoelectric point and interpretation of experimental NPs titration curves.