The pH-induced conformational changes of proteins are systematically studied in the framework of a hydrophobic-polar (HP) model, in which proteins are dramatically simplified as chains of hydrophobic (H) and polar (P) beads on a lattice. We express the electrostatic interaction, the principal driving force of pH-induced unfolding that is not included in the conventional HP model, as the repulsive energy term between P monomers. As a result of the exact enumeration of all of the 14- to 18-mers, it is found that lowest-energy states in many sequences change from single "native" conformations to multiple sets of "denatured" conformations with an increase in the electrostatic repulsion. The switching of the lowest-energy states occurs in quite a similar way to real proteins: it is almost always between two states, while in a small fraction of ≥16-mers it is between three states. We also calculate the structural fluctuations for all of the denatured states and find that the denatured states contain a broad range of incompletely unfolded conformations, similar to "molten globule" states referred to in acid or alkaline denatured real proteins. These results show that the proposed model provides a simple physical picture of pH-induced protein denaturation.