An experimental-theoretical approach for the elucidation of protein stability is proposed. The theoretical prediction of pH-dependent protein stability is based on the macroscopic electrostatic model for calculation of the pH-dependent electrostatic free energy of proteins. As a test of the method we have considered the pH-dependent stability of sperm whale metmyoglobin. Two theoretical methods for evaluation of the electrostatic free energy and p K values are applied: the finite-difference Poisson-Boltzmann method and the semiempirical approach based on the modified Tanford-Kirkwood theory. The theoretical results for electrostatic free energy of unfolding are compared with the experimental data for guanidine hydrochloride unfolding under equilibrium conditions over a wide pH range. Using the optical parameters of the Soret absorbance to monitor conformational equilibrium and Tanford's method to estimate the resulting data, it was found that the conformational free energy of unfolding of metmyoglobin is 16.3 kcal mol(-1) at neutral pH values. The total unfolding free energies were calculated on the basis of the theoretically predicted electrostatic unfolding free energies and the experimentally measured midpoints (pH(1/2)) of acidic and alkaline denaturation transitions. Experimental data for alkaline denaturation were used for the first time in theoretical analysis of the pH-dependent unfolding of myoglobin. The present results demonstrate that the simultaneous application of appropriate theoretical and experimental methods permits a more complete analysis of the pH-dependent and pH-independent properties and stability of globular proteins.