Glutathione peroxidase revisited--simulation of the catalytic cycle by computer-assisted molecular modelling

Glutathione peroxidase, the first example of selenoproteins identified in mammals, was subjected to force field calculations and molecular dynamics in order to enable a clearer comprehension of enzymatic selenium catalysis. Starting from the established X-ray structure of bovine GPX, all kinetically defined intermediates and enzyme substrate complexes were modelled. The models thus obtained support the hypothesis that the essential steps of the catalysis are three distinct redox changes of the active site selenium which, in the ground state, presents itself at the surface of selenoperoxidases as the center of a characteristic triad built by selenocysteine, glutamine and tryptophan. In GPX, four arginine residues and a lysine residue provide an electrostatic architecture which, in each reductive step, directs the donor substrate GSH towards the catalytic center in such a way that its sulfhydryl group must react with the selenium moiety. To this end, different equally efficient modes of substrate binding appear possible. The models are consistent with substrate specificity data, kinetic pattern and other functional characteristics of the enzyme. Comparison of molecular models of GPX with those of other members of the GPX superfamily reveals that the cosubstrate binding mechanisms are unique for the classical type of cytosolic glutathione peroxidases but cannot operate e. g. in plasma GPX and phospholipid hydroperoxide GPX. The structural differences between the selenoperoxidases, shown to be relevant to their specificities, are discussed in terms of functional diversification within the GPX superfamily.