Alkylation of zinc-bound thiolates occurs in both catalytic and structural zinc sites of enzymes. Recent biomimetic studies have led to a controversy as to which mechanism is operative in thiolate alkylation. Building on one of these biomimetic complexes, we have devised a series of models that allow for an appraisal of the roles of charge, ligand nature, and hydrogen bonding to sulfur on reactivity. The reactions of these complexes with methyl iodide, leading to thioethers and zinc iodide complexes, have been examined by density functional theory calculations, in the gas phase as well as in an aqueous solution. In all cases, a S(N)2 reaction is favored over sigma-bond metathesis. Both the net electronic charge and the hydrogen bond play a significant role in the nucleophilicity of the thiolate. We find that the mechanistic diversity observed experimentally can be explained by the difference in the net charge of the complexes. A dianionic complex follows a dissociative pathway, whereas an associative one is preferred for a neutral system.

Density Functional Theory and post-Hartree Fock calculations reveal an unusual energy profile for Zn-S and Zn-N bond dissociation reactions in several [Zn(SR)(4)](2-) and [Zn(Im)(SR)(3)](-) complexes. The Zn-S bond dissociation in tetrathiolate dianions, which is highly exothermic in the gas phase, proceeds through a late transition state which can be rationalized on the basis of an avoided crossing resulting from Coulomb repulsion between the anionic fragments and ligand-to-metal charge-transfer in the [Zn(SR)(4)](2-) complexes. When solvation models for water, DMSO, or acetonitrile are included, some complexes become stable while others are metastable, so this constitutes the first theoretical model which is in full agreement with the experimental data for various [Zn(SR)(4)](2-), [Zn(SR)(3)](-), and [Zn(Im)(SR)(3)](-) complexes. The analysis given here indicates that the Zn(Cys)(4) and Zn(His)(Cys)(3) cores of numerous proteins are metastable with respect to Zn-S and Zn-N bond dissociation, respectively. This is consistent with the kinetic lability at the zinc-centers and illustrates that in nature, thermodynamic stability is imparted upon the zinc cores by the protein environment.