Density functional theory and solvent-assisted proton exchange methods have been applied for computational modeling of the catalytic cycle of selenol zwitterion anion from the kinetic and thermodynamic viewpoints. Selenol zwitterion anion has been represented as an effective glutathione peroxidase nanomimic. It reduces peroxides through a three-step pathway. In the first step, seleninic acid is produced through deprotonating of the selenol zwitterion anion in the presence of the hydrogen peroxide. Seleninic acid reacts with a thiol to form selenylsulfide in the second step. In the last step, selenylsulfide is reduced by the second thiol and regenerates selenolate anion through disulfide formation. Selenol zwitterion anion in comparison to more widely studied compounds such as ebselen has a good activity to react with hydrogen peroxide and producing seleninic acid. The energy barrier of this reaction is 11.7 kcal.mol-1 which is smaller than the reported enzyme mimics. Moreover, the reactions of seleninic acid and selenylsulfide with 2 methanethiol, which is used as a nucleophile, are exothermic by -18.4 or -57.0 kcal.mol-1, respectively. Based on the global electron density transfer value of - 0.507 e from the natural atomic charge analysis, an electronic charge depletion at the transition state (TS), electron-donor substitutions on the selenolate facilitates the reduction reaction, effectively. Finally, the nature of the bond formation/cleavage at the TS has been quantitatively described by using the topological analyses.

Keywords