Elucidation of the Mechanism of Selenoprotein Glutathione Peroxidase (GPx)-Catalyzed Hydrogen Peroxide Reduction by Two Glutathione Molecules: A Density Functional Study

Structural basis for a distinct catalytic mechanism in Trypanosoma brucei tryparedoxin peroxidase

Trypanosoma brucei, the causative agent of African sleeping sickness, encodes three cysteine homologues (Px I-III) of classical selenocysteine-containing glutathione peroxidases. The enzymes obtain their reducing equivalents from the unique trypanothione (bis(glutathionyl)spermidine)/tryparedoxin system. During catalysis, these tryparedoxin peroxidases cycle between an oxidized form with an intramolecular disulfide bond between Cys(47) and Cys(95) and the reduced peroxidase with both residues in the thiol state. Here we report on the three-dimensional structures of oxidized T. brucei Px III at 1.4A resolution obtained by x-ray crystallography and of both the oxidized and the reduced protein determined by NMR spectroscopy. Px III is a monomeric protein unlike the homologous poplar thioredoxin peroxidase (TxP). The structures of oxidized and reduced Px III are essentially identical in contrast to what was recently found for TxP. In Px III, Cys(47), Gln(82), and Trp(137) do not form the catalytic triad observed in the selenoenzymes, and related proteins and the latter two residues are unaffected by the redox state of the protein. The mutational analysis of three conserved lysine residues in the vicinity of the catalytic cysteines revealed that exchange of Lys(107) against glutamate abrogates the reduction of hydrogen peroxide, whereas Lys(97) and Lys(99) play a crucial role in the interaction with tryparedoxin.

Selenium in chemistry and biochemistry in comparison to sulfur

What makes selenoenzymes--seen from a chemist's view--so special that they cannot be substituted by just more analogous or adapted sulfur proteins? This review compiles and compares physicochemical properties of selenium and sulfur, synthetic routes to selenocysteine (Sec) and its peptides, and comparative studies of relevant thiols and selenols and their (mixed) dichalcogens, required to understand the special role of selenium in selenoproteins on the atomic molecular level. The biochemically most relevant differences are the higher polarizability of Se- and the lower pKa of SeH. The latter has a strikingly different pH-dependence than thiols, with selenols being active at much lower pH. Finally, selected typical enzymatic mechanisms which involve selenocysteine are critically discussed, also in view of the authors' own results.

How similar are enzyme active site geometries derived from quantum mechanical theozymes to crystal structures of enzyme-inhibitor complexes? Implications for enzyme design

Quantum mechanical optimizations of theoretical enzymes (theozymes), which are predicted catalytic arrays of biological functionalities stabilizing a transition state, have been carried out for a set of nine diverse enzyme active sites. For each enzyme, the theozyme for the rate-determining transition state plus the catalytic groups modeled by side-chain mimics was optimized using B3LYP/6-31G(d) or, in one case, HF/3-21G(d) quantum mechanical calculations. To determine if the theozyme can reproduce the natural evolutionary catalytic geometry, the positions of optimized catalytic atoms, i.e., covalent, partial covalent, or stabilizing interactions with transition state atoms, are compared to the positions of the atoms in the X-ray crystal structure with a bound inhibitor. These structure comparisons are contrasted to computed substrate-active site structures surrounded by the same theozyme residues. The theozyme/transition structure is shown to predict geometries of active sites with an average RMSD of 0.64 A from the crystal structure, while the RMSD for the bound intermediate complexes are significantly higher at 1.42 A. The implications for computational enzyme design are discussed.

DFT study of the glutathione peroxidase-like activity of phenylselenol incorporating solvent-assisted proton exchange

Modeling of the glutathione peroxidase-like activity of phenylselenol has been accomplished using density-functional theory and solvent-assisted proton exchange (SAPE). SAPE is a modeling technique intended to mimic solvent participation in proton transfer associated with chemical reaction. Within this method, explicit water molecules incorporated into the gas-phase model allow relay of a proton through the water molecules from the site of protonation in the reactant to that in the product. The activation barriers obtained by SAPE for the three steps of the GPx-like mechanism of PhSeH fall within the limits expected for a catalytic system at physiological temperatures (DeltaG(1)++ = 19.1 kcal/mol; DeltaG(2)++= 6.6 kcal/mol; G(3)++ = 21.7 kcal/mol) and are significantly lower than studies which require direct proton transfer. The size of the SAPE network is also considered for the model of the reduction of the selenenic acid, step 2 of the GPx-like cycle. Use of a four-water network better accommodates the reaction pathway and reduces the activation barrier by 5 kcal/mol over the two-water model.

Crystal Structures of a Poplar Thioredoxin Peroxidase that Exhibits the Structure of Glutathione Peroxidases: Insights into Redox-driven Conformational Changes

Glutathione peroxidases (GPXs) are a group of enzymes that regulate the levels of reactive oxygen species in cells and tissues, and protect them against oxidative damage. Contrary to most of their counterparts in animal cells, the higher plant GPX homologues identified so far possess cysteine instead of selenocysteine in their active site. Interestingly, the plant GPXs are not dependent on glutathione but rather on thioredoxin as their in vitro electron donor. We have determined the crystal structures of the reduced and oxidized form of Populus trichocarpaxdeltoides GPX5 (PtGPX5), using a selenomethionine derivative. PtGPX5 exhibits an overall structure similar to that of the known animal GPXs. PtGPX5 crystallized in the assumed physiological dimeric form, displaying a pseudo ten-stranded beta sheet core. Comparison of both redox structures indicates that a drastic conformational change is necessary to bring the two distant cysteine residues together to form an intramolecular disulfide bond. In addition, a computer model of a complex of PtGPX5 and its in vitro recycling partner thioredoxin h1 is proposed on the basis of the crystal packing of the oxidized form enzyme. A possible role of PtGPX5 as a heavy-metal sink is also discussed.

Plant glutathione peroxidases are functional peroxiredoxins distributed in several subcellular compartments and regulated during biotic and abiotic stresses

We provide here an exhaustive overview of the glutathione (GSH) peroxidase (Gpx) family of poplar (Populus trichocarpa). Although these proteins were initially defined as GSH dependent, in fact they use only reduced thioredoxin (Trx) for their regeneration and do not react with GSH or glutaredoxin, constituting a fifth class of peroxiredoxins. The two chloroplastic Gpxs display a marked selectivity toward their electron donors, being exclusively specific for Trxs of the y type for their reduction. In contrast, poplar Gpxs are much less specific with regard to their electron-accepting substrates, reducing hydrogen peroxide and more complex hydroperoxides equally well. Site-directed mutagenesis indicates that the catalytic mechanism and the Trx-mediated recycling process involve only two (cysteine [Cys]-107 and Cys-155) of the three conserved Cys, which form a disulfide bridge with an oxidation-redox midpoint potential of -295 mV. The reduction/formation of this disulfide is detected both by a shift on sodium dodecyl sulfate-polyacrylamide gel electrophoresis or by measuring the intrinsic tryptophan fluorescence of the protein. The six genes identified coding for Gpxs are expressed in various poplar organs, and two of them are localized in the chloroplast, with one colocalizing in mitochondria, suggesting a broad distribution of Gpxs in plant cells. The abundance of some Gpxs is modified in plants subjected to environmental constraints, generally increasing during fungal infection, water deficit, and metal stress, and decreasing during photooxidative stress, showing that Gpx proteins are involved in the response to both biotic and abiotic stress conditions.

A DFT study of the mechanism of Ni superoxide dismutase (NiSOD): Role of the active site cysteine-6 residue in the oxidative half-reaction

In the present DFT study, the catalytic mechanism of H2O2 formation in the oxidative half-reaction of NiSOD, E-Ni(II) + O2- + 2H+ --> E-Ni(III) + H2O2, has been investigated. The main objective of this study is to investigate the source of two protons required in this half-reaction. The proposed mechanism consists of two steps: superoxide coordination and H2O2 formation. The effect of protonation of Cys6 and the proton donating roles of side chains (S) and backbones (B) of His1, Asp3, Cys6, and Tyr9 residues in these two steps have been studied in detail. For protonated Cys6, superoxide binding generates a Ni(III)-O2H species in a process that is exothermic by 17.4 kcal/mol (in protein environment using the continuum model). From the Ni(III)-O2H species, H2O2 formation occurs through a proton donation by His1 via Tyr9, which relative to the resting position of the enzyme is exothermic by 4.9 kcal/mol. In this pathway, a proton donating role of His1 residue is proposed. However, for unprotonated Cys6, a Ni(II)-O2- species is generated in a process that is exothermic by 11.3 kcal/mol. From the Ni(II)-O2- species, the only feasible pathway for H2O2 formation is through donation of protons by the Tyr9(S)-Asp3(S) pair. The results discussed in this study elucidate the role of the active site residues in the catalytic cycle and provide intricate details of the complex functioning of this enzyme.

Biomimetic Studies on Selenoenzymes: Modeling the Role of Proximal Histidines in Thioredoxin Reductases

The roles of built-in thiol cofactors and the basic histidine (His) residues in the active site of mammalian thioredoxin reductases (TrxRs) are described with the help of experimental and density functional theory calculations on small-molecule model compounds. The reduction of selenenyl sulfides by thiols in selenoenzymes such as glutathione peroxidase (GPx) and TrxR is crucial for the regeneration of the active site. Experimental as well as theoretical studies were carried out with model selenenyl sulfides to probe their reactivity toward incoming thiols. We have shown that the nucleophilic attack of thiols takes place at the selenium center in the selenenyl sulfides. These thiol exchange reactions would hamper the regeneration of the active species selenol. Therefore, the basic His residues are expected to play crucial roles in the selenenyl sulfide state of TrxR. Our model study with internal amino groups in the selenenyl sulfide state reveals that the basic His residues may play important roles by deprotonating the thiol moiety in the selenenic acid state and by interacting with the sulfur atom in the selenenyl sulfide state to facilitate the nucleophilic attack of thiol at sulfur rather than at selenium, thereby generating the catalytically active species selenol. This model study also suggests that the enzyme may use the internal cysteines as cofactors to overcome the thiol exchange reactions.

Is the Protein Surrounding the Active Site Critical for Hydrogen Peroxide Reduction by Selenoprotein Glutathione Peroxidase? An ONIOM Study

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.