Quantum mechanical treatment of As3+-thiol model compounds: implication for the core structure of As(III)-metallothionein

Exposure to inorganic arsenic (As) is one of the major health concerns in several regions around the world. Binding of As(III) with thiols is central to the mechanisms related to its toxicity, detoxification, and therapeutic effects. Due to its high thiol content, metallothionein (MT) is presumed to play an important role in case of arsenic toxicity. Consequences of these As-thiol interactions are not yet clear due to various difficulties in the characterization of arsenic bound proteins by spectroscopic techniques. Computational modeling can be a reliable approach in predicting the molecular structures of such complexes. This paper presents the results of a systematic study on different As(III)-thiol model compounds conducted by both ab initio and DFT methods with different Gaussian type basis sets. Proficiency of these theoretical methods has been evaluated in terms of bond lengths, bond angles, free energy, partial atomic charges, computational cost, and comparison with the experimental data. It has been demonstrated that the DFT-B3LYP/6-311+G(3df) functional offers better accuracy in predicting the structure and the UV absorption spectra of As(III)-thiol complexes. The results of the present study also helps in defining the boundaries for the core of arsenic bound MT so that quantum mechanical/molecular mechanical (QM/MM) methods can be employed to predict the structural and functional aspects of the protein.