Selenocysteine versus cysteine reactivity: a theoretical study of their oxidation by hydrogen peroxide

The cysteine and selenocysteine oxidation by H2O2 in vacuo and in aqueous solution was studied using the integrated molecular orbital + molecular orbital (IMOMO) method combining the quadratic configuration method QCISD(T) and the spin projection of second-order perturbation theory PMP2. It is shown that including in the model system of cysteine (selenocysteine) residue up to 20 atoms has significant consequences upon the calculated reaction energy barrier. On the other hand, it is demonstrated that free cysteine and selenocysteine have very similar reaction energy barriers, 77-79 kJ mol(-1) in aqueous solution. It is thus concluded that the high experimental reaction rate constant reported for the oxidation of the selenocysteine residue in the glutathione peroxidase (GPx) active center is due to an important interaction between selenocysteine and its molecular environment. The sensitivity of the calculated energy barrier to the dielectric constant of the molecular environment observed for both cysteine and selenocysteine as well as the catalytic effect of the NH group emphasized in the case of cysteine supports this hypothesis.