Modeling the Glutathione Peroxidase-Like Activity of a Cyclic Seleninate by DFT and Solvent-Assisted Proton Exchange

Regulation of Tyrosinase Enzyme Activity by Glutathione Peroxidase Mimics

Hydrogen peroxide plays a crucial role in the melanogenesis process by regulating the activity of the key melanin-forming enzyme tyrosinase, responsible for the browning of fruits, vegetables, and seafood. Therefore, a molecule with dual activities, both efficient tyrosinase inhibition and strong hydrogen peroxide degrading ability, may act as a promising antibrowning agent. Herein, we report highly efficient selone-based mushroom tyrosinase inhibitors 2 and 3 with remarkable glutathione peroxidase (GPx) enzyme-like activity. The presence of benzimidazole moiety enhances the tyrosinase inhibition efficiency of selone 2 (IC₅₀ = 0.4 μM) by almost 600 times higher than imidazole-based selone 1 (IC₅₀ = 238 μM). Interestingly, the addition of another aromatic ring to the benzimidazole moiety has led to the development of an efficient lipid-soluble tyrosinase inhibitor 3 (IC₅₀ = 2.4 μM). The selenium center and the -NH group of 2 and 3 are extremely crucial to exhibit high GPx-like activity and tyrosinase inhibition potency. The hydrophobic moiety of the inhibitors (2 and 3) further assists them in tightly binding at the active site of the enzyme and facilitates the C═Se group to strongly coordinate with the copper ions. Inhibitor 2 exhibited excellent antibrowning and polyphenol oxidase inhibition properties in banana and apple juice extracts.

Non-covalent interactions in the glutathione peroxidase active center and their influence on the enzyme activity

In this computational work, the structure of the active center of the enzyme glutathione peroxidase (in three forms -SeH, -SeOH and -Se(O)OH) and the non-covalent interactions in it were investigated using modern quantum chemistry methods. The non-covalent interactions are described in detail. The presence of σ-hole interactions (chalcogen, tetrel and pnictogen bonds) formed mostly by a selenium atom as an electrophile in the glutathione peroxidase active center is confirmed for the first time. It is shown that a number of non-covalent interactions stabilize intermediates along the catalytic cycle and that modelling of the whole enzyme active center is necessary for accurate predictions of thermodynamic parameters, in particular, activation barriers.

Selenium-Catalyzed Reduction of Hydroperoxides in Chemistry and Biology

Among the chalcogens, selenium is the key element for catalyzed H2O2 reduction. In organic synthesis, catalytic amounts of organo mono- and di-selenides are largely used in different classes of oxidations, in which H2O2 alone is poorly efficient. Biological hydroperoxide metabolism is dominated by peroxidases and thioredoxin reductases, which balance hydroperoxide challenge and contribute to redox regulation. When their selenocysteine is replaced by cysteine, the cellular antioxidant defense system is impaired. Finally, classes of organoselenides have been synthesized with the aim of mimicking the biological strategy of glutathione peroxidases, but their therapeutic application has so far been limited. Moreover, their therapeutic use may be doubted, because H2O2 is not only toxic but also serves as an important messenger. Therefore, over-optimization of H2O2 reduction may lead to unexpected disturbances of metabolic regulation. Common to all these systems is the nucleophilic attack of selenium to one oxygen of the peroxide bond promoting its disruption. In this contribution, we revisit selected examples from chemistry and biology, and, by using results from accurate quantum mechanical modelling, we provide an accurate unified picture of selenium's capacity of reducing hydroperoxides. There is clear evidence that the selenoenzymes remain superior in terms of catalytic efficiency.

Equilibrium between tri- and tetra-coordinate chalcogenuranes is critical for cysteine protease inhibition

There have been significant advances in the biological use of hypervalent selenium and tellurium compounds as cysteine protease inhibitors. However, the full understanding of their reaction mechanisms for and cysteine proteases inhibition is still elusive. Kinetic studies suggest an irreversible inhibition mechanism, which was explained by forming a covalent bond between the enzyme sulfhydryl group and the chalcogen atom at its hypervalent state (+4). In this work, we performed a theoretical investigation using density functional theory to propose the active inhibitor form in an aqueous solution. To this end, we investigated chloride ligand exchange reactions by oxygen and sulfur nucleophiles on hypervalent selenium and tellurium compounds. All tetra- and tri-coordinated chalcogen compounds and distinct protonation states of the nucleophiles were considered, totaling 34 unique species, 7 nucleophiles, and 155 free energies reactions. We discovered that chloride is easily replaced by a nonprotonated nucleophile (SH^- or OH^- ) in R₂ SeCl₂ . We also found that tri-coordinate species are more stable than their tetra-coordinate counterparts, with selenoxide (R₂ SeO) protonation being strongly exergonic in acid pH. The thermodynamic and kinetic results suggest that the protonated selenoxide (R₂ SeOH⁺ ) is the most probable active chemical species in biological media. The computed energetic profiles paint a possible picture for selenuranes activity, with successive exergonic steps leading to a covalent inhibition of thiol-dependent enzymes, like cysteine proteases. A second pathway has also been uncovered, with a direct reaction to chalcogenonium cation (R₂ SeCl⁺ ) as the inhibition step. Tellurium compounds showed similar trends but formed telluroxide in a pH-independent fashion.

Effect of Methylmercury Binding on the Peroxide-Reducing Potential of Cysteine and Selenocysteine

Methylmercury (CH3Hg+) binding to catalytically fundamental cysteine and selenocysteine of peroxide-reducing enzymes has long been postulated as the origin of its toxicological activity. Only very recently, CH3Hg+ binding to the selenocysteine of thioredoxin reductase has been directly observed [Pickering, I. J. Inorg. Chem., 2020, 59, 2711-2718], but the precise influence of the toxicant on the peroxide-reducing potential of such a residue has never been investigated. In this work, we employ state-of-the-art density functional theory calculations to study the reactivity of molecular models of the free and toxified enzymes. Trends in activation energies are discussed with attention to the biological consequences and are rationalized within the chemically intuitive framework provided by the activation strain model. With respect to the free, protonated amino acids, CH3Hg+ binding promotes oxidation of the S or Se nucleus, suggesting that chalcogenoxide formation might occur in the toxified enzyme, even if the actual rate of peroxide reduction is almost certainly lowered as suggested by comparison with fully deprotonated amino acids models.

New insight into the role of glutathione reductase in glutathione peroxidase-like activity determination by coupled reductase assay: Molecular Docking Study

Previously we have shown that among 15 substituted salicyloyl (2-hydroxybenzoyl) 5-seleninic acids (SSAs) 4 compounds with longer side chains or a cyclohexyl group exhibit no glutathione peroxidase (GPx)-like activity in the coupled reductase assay. Experimental inhibition of glutathione reductase (GR) by the selenenylsulfide (a main intermediate in the catalytic cycle for GPx-like activity determination) of one of the inactive compounds led us to assess the interactions between 15 selenenylsulfide compounds and the active site of GR by molecular docking. Docking results showed that S and Se atoms in selenenylsulfides of the compounds with no GPx-like activity were beyond 5 Å from S atom of Cys-58 or N atom of imidazole ring of His-467 (Root Mean Square Distances for general assessment of 3 major distances were over 4.8 Å) in the active site, so that they could not be catalyzed to be reduced by GR. Furthermore, their docking scores over 89 Kcal/mol meant that the selenenylsulfides were bound too strongly to the active site to leave it, leading eventually to inhibition of GR. We also applied the molecular docking to other GPx mimics such as ebselen, cyclic seleninate esters and di(propylaminomethylphenyl) diselenides to explain the differences in their GPx-like activity depending to the assays used. Our results suggest that the reduction of a selenenylsulfide by GR plays a positive role in GPx-like activity of GPx mimics in the coupled assay and recommended the prediction of possibility and strength of GPx-like activity by molecular docking before entering experimental research.

The role of methoxy substituents in regulating the activity of selenides that serve as spirodioxyselenurane precursors and glutathione peroxidase mimetics

A series of o-(hydroxymethyl)phenyl selenides containing single or multiple methoxy substituents was synthesized, and the rate at which each compound catalyzed the oxidation of benzyl thiol to its disulfide with excess hydrogen peroxide was measured. This assay provided the means for comparing the relative abilities of the selenides to mimic the antioxidant selenoenzyme glutathione peroxidase. The mechanism for catalytic activity involves oxidation of the selenides to their corresponding selenoxides with hydrogen peroxide, cyclization to spirodioxyselenuranes, followed by reduction with two equivalents of thiol to regenerate the original selenide with concomitant disulfide formation. A single p-methoxy group on each aryl moiety afforded the highest catalytic activity, while methoxy groups in the meta position had little effect compared to the unsubstituted selenide, and o-methoxy groups suppressed activity. The installation of multiple methoxy groups on each aryl moiety provided no improvement. These results can be rationalized on the basis of dominating mesomeric and steric effects of the p- and o-substituents, respectively.

Effect of Methoxy Substituents on the Activation Barriers of the Glutathione Peroxidase-Like Mechanism of an Aromatic Cyclic Seleninate

Density functional theory (DFT) models including explicit water molecules have been used to model the redox scavenging mechanism of aromatic cyclic seleninates. Experimental studies have shown that methoxy substitutions affect the rate of scavenging of reactive oxygen species differently depending upon the position. Activities are enhanced in the para position, unaffected in the meta, and decreased in the ortho. DFT calculations show that the activation barrier for the oxidation of the selenenyl sulfide, a proposed key intermediate, is higher for the ortho methoxy derivative than for other positions, consistent with the low experimental conversion rate.