50 Years of Organoselenium Chemistry, Biochemistry and Reactivity: Mechanistic Understanding, Successful and Controversial Stories.

In 1973, two major discoveries changed the face of selenium chemistry: the identification of the first mammal selenoenzyme, glutathione peroxidase 1, and the discovery of the synthetic utility of the so-called selenoxide elimination.  While the chemical mechanism behind the catalytic activity of glutathione peroxidases appears to be mostly unveiled, little is known about the mechanisms of other selenoproteins and, for some of them, even the function lies in the dark. In chemistry, the capacity of organoselenides of catalyzing hydrogen peroxide activation for the practical manipulation of organic functional groups has been largely explored, and some mechanistic details have been clearly elucidated. As a paradox, despite the long-standing experience in the field, the nature of the active oxidant in various reactions still remains matter of debate. While many successes characterize these fields, the pharmacological use of organoselenides still lacks any true application, and while some organoselenides were found to be non-toxic and safe to use, to date no therapeutically approved use was granted. In this review, some fundamental and chronologically aligned topics spanning organoselenium biochemistry, chemistry and pharmacology are discussed, focusing on the current mechanistic picture describing their activity as either bioactive compounds or catalysts.

Five Hypotheses on the Origins of Temperature Dependence of 77Se NMR Shifts in Diselenides.

Notable thermal shifts in diselenides have been documented in 77Se NMR for more than 50 years, but no satisfactory explanation has been found. Here, five hypotheses are considered as possible explanations for the large temperature dependence of the 77Se chemical shifts of diaryl and dialkyl diselenides compared to monoselenides and selenols. Density functional theory calculations are provided to bolster hypotheses and better understand the effects of barrier height and dipole energies. It is proposed that the temperature dependence of diselenide 77Se NMR chemical shifts is due to rotation around the Se-Se bond and sampling of twisted conformers at higher temperatures. The molecular twisting is solvent dependent; here, DMSO-d6 and toluene-d8 were evaluated. No correlation was established between para-substituents on diaryl diselenides and the magnitude of the change in the 77Se NMR shift (Δδ) with temperature.

Diphenyl Diselenide and SARS-CoV-2: in silico Exploration of the Mechanisms of Inhibition of Main Protease (Mpro) and Papain-like Protease (PLpro).

The SARS-CoV-2 pandemic has prompted global efforts to develop therapeutics. The main protease of SARS-CoV-2 (Mpro) and the papain-like protease (PLpro) are essential for viral replication and are key targets for therapeutic development. In this work, we investigate the mechanisms of SARS-CoV-2 inhibition by diphenyl diselenide (PhSe)2 which is an archetypal model of diselenides and a renowned potential therapeutic agent. The in vitro inhibitory concentration of (PhSe)2 against SARS-CoV-2 in Vero E6 cells falls in the low micromolar range. Molecular dynamics (MD) simulations and density functional theory (DFT) calculations [level of theory: SMD-B3LYP-D3(BJ)/6-311G(d,p), cc-pVTZ] are used to inspect non-covalent inhibition modes of both proteases via π-stacking and the mechanism of covalent (PhSe)2 + Mpro product formation involving the catalytic residue C145, respectively. The in vitro CC50 (24.61 µM) and EC50 (2.39 µM) data indicate that (PhSe)2 is a good inhibitor of the SARS-CoV-2 virus replication in a cell culture model. The in silico findings indicate potential mechanisms of proteases' inhibition by (PhSe)2; in particular, the results of the covalent inhibition here discussed for Mpro, whose thermodynamics is approximatively isoergonic, prompt further investigation in the design of antiviral organodiselenides.

Antioxidant Chimeric Molecules: Are Chemical Motifs Additive? The Case of a Selenium-Based Ligand.

We set up an in silico experiment and designed a chimeric compound integrating molecular features from different efficient ROS (Reactive Oxygen Species) scavengers, with the purpose of investigating potential relationships between molecular structure and antioxidant activity. Furthermore, a selenium centre was inserted due to its known capacity to reduce hydroperoxides, acting as a molecular mimic of glutathione peroxidase; finally, since this organoselenide is a precursor of a N-heterocyclic carbene ligand, its Au(I) carbene complex was designed and examined. A validated protocol based on DFT (Density Functional Theory) was employed to investigate the radical scavenging activity of available sites on the organoselenide precursor ((SMD)-M06-2X/6-311+G(d,p)//M06-2X/6-31G(d)), as well as on the organometallic complex ((SMD)-M06-2X/SDD (Au), 6-311+G(d,p)//ZORA-BLYP-D3(BJ)/TZ2P), considering HAT (Hydrogen Atom Transfer) and RAF (Radical Adduct Formation) regarding five different radicals. The results of this case study suggest that the antioxidant potential of chemical motifs should not be considered as an additive property when designing a chimeric compound, but rather that the relevance of a molecular topology is derived from a chemical motif combined with an opportune chemical space of the molecule. Thus, the direct contributions of single functional groups which are generally thought of as antioxidants per se do not guarantee the efficient radical scavenging potential of a molecular species.

Parameter free evaluation of SN2 reaction rates for halide substitution in halomethane.

We estimate the kinetic constants of a series of archetypal SN2 reactions, i.e., the nucleophilic substitutions of halides in halomethane. A parameter free, multiscale approach recently developed [Campeggio et al., Phys. Chem. Chem. Phys., 2020, 22, 3455] is employed. The protocol relies on quantum mechanical calculations for the description of the energy profile along the intrinsic reaction coordinate, which is then mapped onto a reaction coordinate conveniently built for the reactive process. A Kramers-Klein equation is used to describe the stochastic time evolution of the reaction coordinate and its velocity; friction is parameterized using a hydrodynamic model and Kramers theory is used to derive the rate constant of the reaction. The method is here applied to six SN2 reactions in water at 295.15 K, which differ in the nucleophile and the leaving group. The computed reaction rates are in good agreement with the experimental data and correlate well with the trends observed for the activation energies.

Methylmercury Can Facilitate the Formation of Dehydroalanine in Selenoenzymes: Insight from DFT Molecular Modeling.

Experimental studies have indicated that electrophilic mercury forms (e.g., methylmercury, MeHg+) can accelerate the breakage of selenocysteine in vitro. Particularly, in 2009, Khan et al. (Environ. Toxicol. Chem. 2009, 28, 1567-1577) proposed a mechanism for the degradation of a free methylmercury selenocysteinate complex that was theoretically supported by Asaduzzaman et al. (Inorg. Chem. 2010, 50, 2366-2372). However, little is known about the fate of methylmercury selenocysteinate complexes embedded in an enzyme, especially in conditions of oxidative stress in which methylmercury target enzymes operate. Here, an accurate computational study on molecular models (level of theory: COSMO-ZORA-BLYP-D3(BJ)/TZ2P) was carried out to investigate the formation of dehydroalanine (Dha) in selenoenzymes, which irreversibly impairs their function. Methylselenocysteine as well as methylcysteine and methyltellurocysteine were included to gain insight on the peculiar behavior of selenium. Dha forms in a two-step process, i.e., the oxidation of the chalcogen nucleus followed by a syn-elimination leading to the alkene and the chalcogenic acid. The effect of an excess of hydrogen peroxide, which may lead to the formation of chalcogenones before the elimination, and of MeHg+, a severe toxicant targeting selenoproteins, which leads to the formation of methylmercury selenocysteinate, are also studied with the aim of assessing whether these pathological conditions facilitate the formation of Dha. Indeed, elimination occurs after chalcogen oxidation and MeHg+ facilitates the process. These results indicate a possible mechanism of toxicity of MeHg+ in selenoproteins.

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.

Drug Design: Where We Are and Future Prospects.

Medicinal chemistry is facing new challenges in approaching precision medicine. Several powerful new tools or improvements of already used tools are now available to medicinal chemists to help in the process of drug discovery, from a hit molecule to a clinically used drug. Among the new tools, the possibility of considering folding intermediates or the catalytic process of a protein as a target for discovering new hits has emerged. In addition, machine learning is a new valuable approach helping medicinal chemists to discover new hits. Other abilities, ranging from the better understanding of the time evolution of biochemical processes to the comprehension of the biological meaning of the data originated from genetic analyses, are on their way to progress further in the drug discovery field toward improved patient care. In this sense, the new approaches to the delivery of drugs targeted to the central nervous system, together with the advancements in understanding the metabolic pathways for a growing number of drugs and relating them to the genetic characteristics of patients, constitute important progress in the field.

Selenoxide Elimination Triggers Enamine Hydrolysis to Primary and Secondary Amines: A Combined Experimental and Theoretical Investigation.

We discuss a novel selenium-based reaction mechanism consisting in a selenoxide elimination-triggered enamine hydrolysis. This one-pot model reaction was studied for a set of substrates. Under oxidative conditions, we observed and characterized the formation of primary and secondary amines as elimination products of such compounds, paving the way for a novel strategy to selectively release bioactive molecules. The underlying mechanism was investigated using NMR, mass spectrometry and density functional theory (DFT).

Chalcogen-mercury bond formation and disruption in model Rabenstein's reactions: A computational analysis.

Methylmercury is a highly toxic compound and human exposure is mainly related to consumption of polluted fish and seafood. The inactivation of thiol-based enzymes, promoted by the strong affinity binding of electrophilic mercuric ions to thiol and selenol groups of proteins, is likely an important factor explaining its toxicity. A key role is played by the chemistry and reactivity of the mercury-chalcogens bond, particularly HgS and HgSe, which is the focus of this computational work (level of theory: (COSMO)-ZORA-BLYP-D3(BJ)/TZ2P). We analyze nine ligand-exchange model reactions (the so-called Rabenstein's reactions) involving an entering ligand (methylchalcogenolate) and a substrate (methylchalcogenolatemethylmercury). Trends in reaction and activation energies are discussed and a change in mechanism is reported for all cases when going from gas phase to water, that is, from a single-well potential energy surface (PES) to a canonical SN 2-like mechanism. The reasons accounting for the biochemically challenging and desired displacement of methylmercury from a seleno/thiol protein can be found already in these model reactions, as can be seen from the similarities of the ligand exchange reactions in solution in thermodynamics and kinetics.

Human elastin-like polypeptides as a versatile platform for exploitation of ultrasensitive bilirubin detection by UnaG.

A new, bifunctional recombinant protein was expressed as the fusion product of human elastin-like polypeptide (HELP) and the bilirubin-binding protein UnaG. The engineered product displays both the HELP-specific property of forming a functional hydrogel matrix and the UnaG-specific capacity of emitting green fluorescence upon ligand binding. The new fusion protein has been proven to be effective at detecting bilirubin in complex environments with high background noise. A cell culture model of the stress response, consisting of bilirubin released in the cell culture medium, was set up to assess the bilirubin-sensing properties of the functional matrix obtained by cross-linking the HELP moiety. Our engineered protein allowed us to monitor cell induction by the release of bilirubin in the culture medium on a nanomolar scale. This study shows that elastin-like protein fusion represents a versatile platform for the development of novel and commercially viable analytical and biosensing devices.

Multiscale modeling of reaction rates: application to archetypal SN2 nucleophilic substitutions.

We propose an approach to the evaluation of kinetic rates of elementary chemical reactions within Kramers' theory based on the definition of the reaction coordinate as a linear combination of natural, pseudo Z-matrix, internal coordinates of the system. The element of novelty is the possibility to evaluate the friction along the reaction coordinate, within a hydrodynamic framework developed recently [J. Campeggio et al., J. Comput. Chem. 2019, 40, 679-705]. This, in turn, allows to keep into account barrier recrossing, i.e. the transmission coefficient that is employed in correcting transition state theory evaluations. To test the capabilities and the flaws of the approach we use as case studies two archetypal SN2 reactions. First, we consider to the standard substitution of chloride ion to bromomethane. The rate constant at 295.15 K is evaluated to k/c⊖ = 2.7 × 10-6 s-1 (with c⊖ = 1 M), which compares well to the experimental value of 3.3 × 10-6 s-1 [R. H. Bathgate and E. A. Melwyn-Hughes, J. Chem. Soc 1959, 2642-2648]. Then, the method is applied to the SN2 reaction of methylthiolate to dimethyl disulfide in water. In biology, such an interconversion of thiols and disulfides is an important metabolic topic still not entirely rationalized. The predicted rate constant is k/c⊖ = 7.7 × 103 s-1. No experimental data is available for such a reaction, but it is in accord with the fact that the alkyl thiolates to dialkyl disulfides substitutions in water have been found to be fast reactions [S. M. Bachrach, J. M. Hayes, T. Dao and J. L. Mynar, Theor. Chem. Acc. 2002, 107, 266-271].

Radical scavenging activity of natural antioxidants and drugs: Development of a combined machine learning and quantum chemistry protocol.

Many natural substances and drugs are radical scavengers that prevent the oxidative damage to fundamental cell components. This process may occur via different mechanisms, among which, one of the most important, is hydrogen atom transfer. The feasibility of this process can be assessed in silico using quantum mechanics to compute ΔGHAT ○. This approach is accurate, but time consuming. The use of machine learning (ML) allows us to reduce tremendously the computational cost of the assessment of the scavenging properties of a potential antioxidant, almost without affecting the quality of the results. However, in many ML implementations, the description of the relevant features of a molecule in a machine-friendly language is still the most challenging aspect. In this work, we present a newly developed machine-readable molecular representation aimed at the application of automatized ML algorithms. In particular, we show an application on the calculation of ΔGHAT ○.

The 125Te Chemical Shift of Diphenyl Ditelluride: Chasing Conformers over a Flat Energy Surface.

The interest in diphenyl ditelluride (Ph2Te2) is related to its strict analogy to diphenyl diselenide (Ph2Se2), whose capacity to reduce organic peroxides is largely exploited in catalysis and green chemistry. Since the latter is also a promising candidate as an antioxidant drug and mimic of the ubiquitous enzyme glutathione peroxidase (GPx), the use of organotellurides in medicinal chemistry is gaining importance, despite the fact that tellurium has no recognized biological role and its toxicity must be cautiously pondered. Both Ph2Se2 and Ph2Te2 exhibit significant conformational freedom due to the softness of the inter-chalcogen and carbon⁻chalcogen bonds, preventing the existence of a unique structure in solution. Therefore, the accurate calculation of the NMR chemical shifts of these flexible molecules is not trivial. In this study, a detailed structural analysis of Ph2Te2 is carried out using a computational approach combining classical molecular dynamics and relativistic density functional theory methods. The goal is to establish how structural changes affect the electronic structure of diphenyl ditelluride, particularly the 125Te chemical shift.

Nature and strength of chalcogen-π bonds.

Chalcogen-π interactions occur between a covalently bound chalcogen atom that enters into a non-covalent interaction with an unsaturated moiety, a bonding motif found in various structures, such as, proteins. In this work, we have systematically explored and analyzed chalcogen-π interactions in model systems X2DA (with D = O, S, Se, Te; X = halogen; A = acetylene, ethylene and 2-butyne), using relativistic density functional theory (DFT). The nature and trends in stability of the chalcogen-π bonds are analyzed and interpreted in terms of quantitative MO theory in combination with a matching canonical energy decomposition analysis (EDA) scheme. We find that chalcogen-π bonds increase in strength as the X-D electronegativity difference becomes greater. Moreover, 2-butyne was found to participate in the strongest non-covalent interaction due to enhanced orbital interactions.

Oxidation of organic diselenides and ditellurides by H2O2 for bioinspired catalyst design.

The reactivity of diselenides and ditellurides of general formula (RX)2 (X = Se, Te; R = H, CH3, Ph) toward hydrogen peroxide was studied through a computational approach based on accurate Density Functional Theory (DFT) calculations. The aliphatic and aromatic dichalcogenides have been chosen in light of their activity in glutathione peroxidase (GPx)-like catalytic cycles and their promising features as efficient antioxidant compounds. The reaction products, the energetics and the mechanistic details of these oxidations are discussed. Analogous disulfides are included in our analysis for completeness. We find that the barrier for oxidation of dichalcogenides decreases from disulfides to diselenides to ditellurides. On the other hand, variation of the substituents at the chalcogen nucleus has relatively little effect on the reactivity.

Mechanistic Insight into the Oxidation of Organic Phenylselenides by H2 O2.

The oxidation of organic phenylselenides by H2 O2 is investigated in model compounds, namely, n-butyl phenyl selenide (PhSe(nBu)), bis(phenylselanyl)methane (PhSeMeSePh), diphenyl diselenide (PhSeSePh), and 1,2-bis(phenylselanyl)ethane (PhSeEtSePh). Through a combined experimental (1 H and 77 Se NMR) and computational approach, we characterize the direct oxidation of monoselenide to selenoxide, the stepwise double oxidation of PhSeMeSePh that leads to different diastereomeric diselenoxides, the complete oxidation of the diphenyldiselenide that leads to selenium-selenium bond cleavage, and the subsequent formation of the phenylseleninic product. The oxidation of PhSeEtSePh also results in the formation of phenylseleninic acid along with 1-(vinylseleninyl)benzene, which is derived from a side elimination reaction. The evidence of a direct mechanism, in addition to an autocatalytic mechanism that emerges from kinetic studies, is discussed. By considering our observations of diselenides with chalcogen atoms that are separated by alkyl spacers of different length, a rationale for the advantage of diselenide versus monoselenide catalysts is presented.

Role of the Chalcogen (S, Se, Te) in the Oxidation Mechanism of the Glutathione Peroxidase Active Site.

The oxidation by H2 O2 of the human phospholipid hydroperoxide glutathione peroxidase (GPx4), used as a model peroxidase selenoenzyme, as well as that of its cysteine (Cys) and tellurocysteine (Tec) mutants, was investigated in silico through a combined classic and quantum mechanics approach to assess the role of the different chalcogens. To perform this analysis, new parameters for selenocysteine (Sec) and tellurocysteine (Tec) were accurately derived for the AMBER ff14SB force field. The oxidation represents the initial step of the antioxidant activity of GPx, which catalyzes the reduction of H2 O2 and organic hydroperoxides by glutathione (GSH). A mechanism involving a charge-separation intermediate is feasible for the Cys and Sec enzymes, leading from the initial thiol/selenol form to sulfenic/selenenic acid, whereas for the Tec mutant a direct oxidation pathway is proposed. Activation strain analyses, performed for Cys-GPx and Sec-GPx, provided insight into the rate-accelerating effect of selenium as compared to sulfur and the role of specific amino acids other than Cys/Sec that are typically conserved in the catalytic pocket.

Antioxidant activity of Zuccagnia-type propolis: A combined approach based on LC-HRMS analysis of bioanalytical-guided fractions and computational investigation.

This study reports a combined approach to assess the antioxidant activity of Zuccagnia-type propolis. Fractions exhibiting the highest antioxidant activities evidenced by DPPH, a ß-carotene bleaching and superoxide radical scavenging activity-non-enzymatic assays, were processed by LC-HRMS/MS to characterize the relevant chemical compounds. A computational protocol based on the DFT calculations was used to rationalize the main outcomes. Among the 28 identified flavonoids, caffeic acids derivatives were in the fraction exhibiting the highest antioxidant activity, with 1-methyl-3-(4'-hydroxyphenyl)-propyl caffeic acid ester and 1-methyl-3-(3',4'-dihydroxyphenyl)-propyl caffeic acid ester as major components. Results clearly showed roles of specific chemical motifs, which can be supported by the computational analysis. This is the first report ascribing the antioxidant ability of Zuccagnia-type propolis to its content in specific caffeic acid derivatives, a potential source of radical scavenging phytochemicals. The proposed protocol can be extended to the study of other plant-products to address the most interesting bioactive compounds.

Radical Scavenging Potential of Ginkgolides and Bilobalide: Insight from Molecular Modeling.

The reactive oxygen species (ROS) scavenging capacities of ginkgolides and bilobalide, which are the peculiar constituents of the extract of Ginkgo biloba, are investigated in silico (level of theory: (SMD)-M06-2X/6-311+G(d,p)//M06-2X/6-31G(d)). Unlike other popular antioxidant natural substances, the carbon backbones of these compounds are entirely aliphatic and exclusively single C-C bonds are present. The selectivity for alkoxyl radicals via hydrogen-atom transfer (HAT) is assessed; importantly, the scavenging of peroxyl radicals is also possible from a peculiar site, here labeled C10 both for ginkgolides and bilobalide. The energetics are described in detail, and the analysis discloses that the studied compounds are powerful scavengers, with thermodynamic and kinetic properties similar to those of Trolox and melatonin, and that, in addition, they display selectivity for peroxyl radicals. These are all chemical-reactivity features contributing to the therapeutic action of the extract of G. biloba.