Based on the mean‐field approximation, we have applied a modified Ising model to describe general protein unfolding behavior at thermodynamic equilibrium with the free energy contributed by the subgroup units (amino acids or peptide bonds) of the protein. With the thermodynamic properties of the protein, this model can associate the stepwise change of an unfolding fraction ratio profile with the local and global conformation unfolding. Taking cytochrome c (cyt c) as a model protein, we have observed, using small‐angle X‐ray scattering and circular dichroism (CD), the global and local structure changes for the protein in three kinds of denaturant environments: acid, urea and guanidine hydrochloride. The small‐angle X‐ray scattering and CD results are mapped to the unfolding fractions as a function of the pH value or denaturant concentration, from which we have extracted local and global unfolding free energies of cyt c in different denaturant environments using a modified Ising model. Based on the characteristics of the thermodynamic properties deduced from the local and global protein folding–unfolding, we discuss the thermodynamic stabilities of the protein in the three denaturant environments, and the possible correlation between the global conformation change of the protein and the local unfolding activities of the S—Fe bond in the Met80‐heme and the α‐helices.