Protein ion charge state distributions in electrospray ionization mass spectra have a potential to provide a wealth of information on protein dynamics, because they contain contributions from all protein conformers present in solution. Such ionic contributions often overlap, limiting the amount of useful information that can be extracted from the spectra. This difficulty is overcome in the present work by using a chemometric approach, which allows spectral deconvolution to be carried out and information on individual protein conformers to be extracted. Experiments are carried out by acquiring a series of spectra over a range of near-native and denaturing conditions to ensure significant changes in the conformers' populations. A total number of protein conformers sampled in the course of the experiments is determined by subjecting the set of collected spectra to singular value decomposition. Ionic contributions of each conformer to the total signal are then determined using a supervised optimization routine. Validation of the method has been carried out by monitoring acid- and alcohol-induced equilibrium states of well-characterized model proteins, chymotrypsin inhibitor 2 (two states), ubiquitin (three states) and apo-myoglobin (four states). For each of the model proteins, a new chemometric procedure yielded a picture of protein dynamics that was in excellent agreement with their documented behavior (as studied with other biophysical tools). The new method appears to be well-suited for monitoring protein dynamics in highly heterogeneous systems consisting of multiple proteins sampling a range of conformational states.