Human metallothionein 1a, a protein with two cysteine-rich metal-binding domains (α with 11 Cys and β with 9), was analyzed in its metal-free form by selective, covalent Cys modification coupled with ESI-MS. The modification profiles of the isolated β- and α-fragments reacted with p-benzoquinone (Bq), N-ethylmalemide (NEM) and iodoacetamide (IAM) were compared with the full length protein using ESI-mass spectral data to follow the reaction pathway. Under denaturing conditions at low pH, the reaction profile with each modifier followed pathways that resulted in stochastic, Normal distributions of species whose maxima was equal to the mol. eq. of modifier added. Our interpretation of modification at this pH is that reaction with the cysteines is unimpeded when the full protein or those of its isolated domains are denatured. At neutral pH, where the protein is expected to be folded in a more compact structure, there is a difference in the larger Bq and NEM modification, whose reaction profiles indicate a cooperative pattern. The reaction profile with IAM under native conditions follows a similar stochastic distribution as at low pH, suggesting that this modifier is small enough to access the cysteines unimpeded by the compact structure. The data emphasize the utility of residue modification coupled with electrospray ionization mass spectrometry for the study of protein structure.