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

ROS-responsive hydrogels: from design and additive manufacturing to biomedical applications

Hydrogels with intricate 3D networks and high hydrophilicity have qualities resembling those of biological tissues, making them ideal candidates for use as smart biomedical materials. Reactive oxygen species (ROS) responsive hydrogels are an innovative class of smart hydrogels, and are cross-linked by ROS-responsive modules through covalent interactions, coordination interactions, or supramolecular interactions. Due to the introduction of ROS response modules, this class of hydrogels exhibits a sensitive response to the oxidative stress microenvironment existing in organisms. Simultaneously, due to the modularity of the ROS-responsive structure, ROS-responsive hydrogels can be manufactured on a large scale through additive manufacturing. This review will delve into the design, fabrication, and applications of ROS-responsive hydrogels. The main goal is to clarify the chemical principles that govern the response mechanism of these hydrogels, further providing new perspectives and methods for designing responsive hydrogel materials.

Reactivity of Zinc Fingers in Oxidizing Environments: Insight from Molecular Models Through Activation Strain Analysis

The reactivity of Zn2+ tetrahedral complexes with H2O2 was investigated in silico, as a first step in their disruption process. The substrates were chosen to represent the cores of three different zinc finger protein motifs, i. e., a Zn2+ ion coordinated to four cysteines (CCCC), to three cysteines and one histidine (CCCH), and to two cysteines and two histidines (CCHH). The cysteine and histidine ligands were further simplified to methyl thiolate and imidazole, respectively. H2O2 was chosen as an oxidizing agent due to its biological role as a metabolic product and species involved in signaling processes. The mechanism of oxidation of a coordinated cysteinate to sulfenate-κS and the trends for the different substrates were rationalized through activation strain analysis and energy decomposition analysis in the framework of scalar relativistic Density Functional Theory (DFT) calculations at ZORA-M06/TZ2P ae // ZORA-BLYP-D3(BJ)/TZ2P. CCCC is oxidized most easily, an outcome explained considering both electrostatic and orbital interactions. The isomerization to sulfenate-κO was attempted to assess whether this step may affect the ligand dissociation; however, it was found to introduce a kinetic barrier without improving the energetics of the dissociation. Lastly, ligand exchange with free thiolates and selenolates was investigated as a trigger for ligand dissociation, possibly leading to metal ejection; molecular docking simulations also support this hypothesis.

Reactions of hypobromous acid with dimethyl selenide, dimethyl diselenide and other organic selenium compounds: kinetics and product formation

Selenium (Se) is an essential micronutrient for many living organisms particularly due to its unique redox properties. We recently found that the sulfur (S) analog for dimethyl selenide (DMSe), i.e. dimethyl sulfide (DMS), reacts fast with the marine oxidant hypobromous acid (HOBr) which likely serves as a sink of marine DMS. Here we investigated the reactivity of HOBr with dimethyl selenide and dimethyl diselenide (DMDSe), which are the main volatile Se compounds biogenically produced in marine waters. In addition, the reactivity of HOBr with further organic Se compounds was tested, i.e., SeMet (as N-acetylated-SeMet), and selenocystine (SeCys2 as N-acetylated-SeCys2), as well as the phenyl-analogs of DMSe and DMDSe, respectively, diphenyl selenide (DPSe) and diphenyl diselenide (DPDSe). Apparent second-order rate constants at pH 8 for the reactions of HOBr with the studied Se compounds were (7.1 ± 0.7) × 107 M-1 s-1 for DMSe, (4.3 ± 0.4) × 107 M-1 s-1 for DMDSe, (2.8 ± 0.3) × 108 M-1 s-1 for SeMet, (3.8 ± 0.2) × 107 M-1 s-1 for SeCys2, (3.5 ± 0.1) × 107 M-1 s-1 for DPSe, and (8.0 ± 0.4) × 106 M-1 s-1 for DPDSe, indicating a very high reactivity of all selected Se compounds with HOBr. The reactivity between HOBr and DMSe is lower than for DMS and therefore this reaction is likely not relevant for marine DMSe abatement. However, the high reactivity of SeMet with HOBr suggests that SeMet may act as a relevant quencher of HOBr.

Role of Group 12 Metals in the Reduction of H2O2 by Santi's Reagent: A Computational Mechanistic Investigation

PhSeZnCl, which is also known as Santi's reagent, can catalyze the reduction of hydrogen peroxide by thiols with a GPx-like mechanism. In this work, the first step of this catalytic cycle, i.e., the reduction of H2O2 by PhSeZnCl, is investigated in silico using state-of-the-art density functional theory calculations. Then, the role of the metal is evaluated by replacing Zn with its group 12 siblings (Cd and Hg). The thermodynamic and kinetic factors favoring Zn are elucidated. Furthermore, the role of the halogen is considered by replacing Cl with Br in all three metal compounds, and this turns out to be negligible. Finally, the overall GPx-like mechanism of PhSeZnCl and PhSeZnBr is discussed by evaluating the energetics of the mechanistic path leading to the disulfide product.

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 (M^pro) and the papain-like protease (PL^pro) 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)₂ which is an archetypal model of diselenides and a renowned potential therapeutic agent. The in vitro inhibitory concentration of (PhSe)₂ 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)₂ + M^pro product formation involving the catalytic residue C145, respectively. The in vitro CC₅₀ (24.61 μM) and EC₅₀ (2.39 μM) data indicate that (PhSe)₂ 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)₂; in particular, the results of the covalent inhibition here discussed for M^pro, whose thermodynamics is approximatively isoergonic, prompt further investigation in the design of antiviral organodiselenides.

Voyage of selenium from environment to life: Beneficial or toxic?

Selenium (Se), a naturally occurring metalloid, is an essential micronutrient for life as it is incorporated as selenocysteine in proteins. Although beneficial at low doses, Se is hazardous at high concentrations and poses a serious threat to various ecosystems. Due to this contrasting 'dual' nature, Se has garnered the attention of researchers wishing to unravel its puzzling properties. In this review, we describe the impact of selenium's journey from environment to diverse biological systems, with an emphasis on its chemical advantage. We describe the uneven distribution of Se and how this affects the bioavailability of this element, which, in turn, profoundly affects the habitat of a region. Once taken up, the subsequent incorporation of Se into proteins as selenocysteine and its antioxidant functions are detailed here. The causes of improved protein function due to the incorporation of redox-active Se atom (instead of S) are examined. Subsequently, the reasons for the deleterious effects of Se, which depend on its chemical form (organo-selenium or the inorganic forms) in different organisms are elaborated. Although Se is vital for the function of many antioxidant enzymes, how the pro-oxidant nature of Se can be potentially exploited in different therapies is highlighted. Furthermore, we succinctly explain how the presence of Se in biological systems offsets the toxic effects of heavy metal mercury. Finally, the different avenues of research that are fundamental to expand our understanding of selenium biology are suggested.

Radical Scavenging Potential of the Phenothiazine Scaffold: A Computational Analysis

The reactivity of phenothiazine (PS), phenoselenazine (PSE), and phenotellurazine (PTE) with different reactive oxygen species (ROS) has been studied using density functional theory (DFT) in combination with the QM‐ORSA (Quantum Mechanics‐based Test for Overall Free Radical Scavenging Activity) protocol for an accurate kinetic rate calculation. Four radical scavenging mechanisms have been screened, namely hydrogen atom transfer (HAT), radical adduct formation (RAF), single electron transfer (SET), and the direct oxidation of the chalcogen atom. The chosen ROS are HO^., HOO^., and CH₃OO^.. PS, PSE, and PTE exhibit an excellent antioxidant activity in water regardless of the ROS due to their characteristic diffusion‐controlled regime processes. For the HO^. radical, the primary active reaction mechanism is, for all antioxidants, RAF. But, for HOO^. and CH₃OO^., the dominant mechanism strongly depends on the antioxidant: HAT for PS and PSE, and SET for PTE. The scavenging efficiency decreases dramatically in lipid environment and remains only significant (via RAF) for the most reactive radical (HO^.). Therefore, PS, PSE, and PTE are excellent antioxidant molecules, especially in aqueous, physiological environments where they are active against a broad spectrum of harmful radicals. There is no advantage or significant difference in the scavenging efficiency when changing the chalcogen since the reactivity mainly derives from the amino hydrogen and the aromatic sites.

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.

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).

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.

Chalcogen-Nitrogen Bond: Insights into a Key Chemical Motif

Chalcogen-nitrogen chemistry deals with systems in which sulfur, selenium, or tellurium is linked to a nitrogen nucleus. This chemical motif is a key component of different functional structures, ranging from inorganic materials and polymers, to rationally designed catalysts, to bioinspired molecules and enzymes. The formation of a selenium–nitrogen bond, typically occurring upon condensation of an amine and the unstable selenenic acid, often leading to intramolecular cyclizations, and its disruption, mainly promoted by thiols, are rather common events in organic Se-catalyzed processes. In this work, focusing on examples taken from selenium organic chemistry and biochemistry, the selenium–nitrogen bond is described, and its strength and reactivity are quantified using accurate computational methods applied to model molecular systems. The intermediate strength of the Se–N bond, which can be tuned to necessity, gives rise to significant trends when comparing it to the stronger S– and weaker Te–N bonds, reaffirming also in this context the peculiar and valuable role of selenium in chemistry and life.

Proton Transfer and SN 2 Reactions as Steps of Fast Selenol and Thiol Oxidation in Proteins: A Model Molecular Study Based on GPx

The so-called peroxidatic cysteines and selenocysteines in proteins reduce hydroperoxides through a dual attack to the peroxide bond in a two-step mechanism. First, a proton dislocation from the thiol/selenol to a close residue of the enzymatic pocket occurs. Then, a nucleophilic attack of the anionic cysteine/selenocysteine to one O atom takes place, while the proton is shuttled back to the second O atom, promoting the formation of a water molecule. In this computational study, we use a molecular model of GPx to demonstrate that the enzymatic environment significantly lowers the barrier of the latter S_N 2 step. Particularly, in our Se-based model the energy barriers for the two steps are 29.82 and 2.83 kcal mol^-1 , both higher than the corresponding barriers computed in the enzymatic cluster, i. e., 21.60 and null, respectively. Our results, obtained at SMD-B3LYP-D3(BJ)/6-311+G(d,p), cc-pVTZ//B3LYP-D3(BJ)/6-311G(d,p), cc-pVTZ level of theory, show that the mechanistic details can be well reproduced using an oversimplified model, but the energetics is definitively more favorable in the GPx active site. In addition, we pinpoint the role of the chalcogen in the peroxide reduction process, rooting the advantages of the presence of selenium in its acidic and nucleophilic properties.

Low-Dose X-ray-Responsive Diselenide Nanocarriers for Effective Delivery of Anticancer Agents

X-ray-responsive nanocarriers for anticancer drug delivery have shown great promise for enhancing the efficacy of chemoradiotherapy. A critical challenge remains for development of such radiation-controlled drug delivery systems (DDSs), which is to minimize the required X-ray dose for triggering the cargo release. Herein, we design and fabricate an effective DDS based on diselenide block copolymers (as nanocarrier), which can be triggered to release their cargo with a reduced radiation dose of 2 Gy due to their sensitivity to both X-ray and the high level of reactive oxygen species (ROS) in the microenvironment of cancer cells. The underlying molecular mechanism is further illustrated by proton nuclear magnetic resonance (¹H NMR) experiments and density functional theory (DFT) calculations. In vivo experiments on tumor-bearing mice validated that the loaded drugs are effectively delivered to the tumor site and exert remarkable antitumor effects (minimum tumor volume/weight) along with X-ray. Furthermore, the diselenide nanocarriers exhibit no noticeable cytotoxicity. These findings provide new insights for the de novo design of radiation-controlled DDSs for cancer chemoradiotherapy.

Antioxidant Potential of Psychotropic Drugs: From Clinical Evidence to In Vitro and In Vivo Assessment and toward a New Challenge for in Silico Molecular Design

Due to high oxygen consumption, the brain is particularly vulnerable to oxidative stress, which is considered an important element in the etiopathogenesis of several mental disorders, including schizophrenia, depression and dependencies. Despite the fact that it is not established yet whether oxidative stress is a cause or a consequence of clinic manifestations, the intake of antioxidant supplements in combination with the psychotropic therapy constitutes a valuable solution in patients' treatment. Anyway, some drugs possess antioxidant capacity themselves and this aspect is discussed in this review, focusing on antipsychotics and antidepressants. In the context of a collection of clinical observations, in vitro and in vivo results are critically reported, often highlighting controversial aspects. Finally, a new challenge is discussed, i.e., the possibility of assessing in silico the antioxidant potential of these drugs, exploiting computational chemistry methodologies and machine learning. Despite the physiological environment being incredibly complex and the detection of meaningful oxidative stress biomarkers being all but an easy task, a rigorous and systematic analysis of the structural and reactivity properties of antioxidant drugs seems to be a promising route to better interpret therapeutic outcomes and provide elements for the rational design of novel drugs.

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].

Fluoxetine scaffold to design tandem molecular antioxidants and green catalysts

Fluoxetine finds application in the treatment of depression and mood disorders. This selective serotonin-reuptake inhibitor (SSRI) also contrasts oxidative stress by direct ROS scavenging, modulation of the endogenous antioxidant defense system, and/or enhancement of the serotonin antioxidant capacity. We synthesised some fluoxetine analogues incorporating a selenium nucleus, thus expanding its antioxidant potential by enabling a hydroperoxides-inactivating, glutathione peroxidase (GPx)-like activity. Radical scavenging and peroxidatic activity were combined in a water-soluble, drug-like, tandem antioxidant molecule. Selenofluoxetine derivatives were reacted with H₂O₂ in water, and the mechanistic details of the reaction were unravelled combining nuclear magnetic resonance (NMR), electrospray ionisation-mass spectrometry (ESI-MS) and quantum chemistry calculations. The observed oxidation–elimination process led to the formation of seleninic acid and cinnamylamine in a trans-selective manner. This mechanism is likely to be extended to other substrates for the preparation of unsaturated cinnamylamines.

We modified fluoxetine by incorporating a selenium nucleus enabling a hydroperoxide-inactivating, glutathione peroxidase (GPx)-like activity and paving the way for its use as green catalyst.

Sequential oxidations of phenylchalcogenides by H2O2: insights into the redox behavior of selenium via DFT analysis

The biological activity of sulfur and selenium, despite their similarity, shows some remarkable differences that have been recognized in many different scenarios.

Major Depressive Disorder and Oxidative Stress: In Silico Investigation of Fluoxetine Activity against ROS

Major depressive disorder is a psychiatric disease having approximately a 20% lifetime prevalence in adults in the United States (U.S.), as reported by Hasin et al. in JAMA Psichiatry 2018 75, 336–346. Symptoms include low mood, anhedonia, decreased energy, alteration in appetite and weight, irritability, sleep disturbances, and cognitive deficits. Comorbidity is frequent, and patients show decreased social functioning and a high mortality rate. Environmental and genetic factors favor the development of depression, but the mechanisms by which stress negatively impacts on the brain are still not fully understood. Several recent works, mainly published during the last five years, aim at investigating the correlation between treatment with fluoxetine, a non-tricyclic antidepressant drug, and the amelioration of oxidative stress. In this work, the antioxidant activity of fluoxetine was investigated using a computational protocol based on the density functional theory approach. Particularly, the scavenging of five radicals (HO•, HOO•, CH3OO•, CH2=CHOO•, and CH3O•) was considered, focusing on hydrogen atom transfer (HAT) and radical adduct formation (RAF) mechanisms. Thermodynamic as well as kinetic aspects are discussed, and, for completeness, two metabolites of fluoxetine and serotonin, whose extracellular concentration is enhanced by fluoxetine, are included in our analysis. Indeed, fluoxetine may act as a radical scavenger, and exhibits selectivity for HO• and CH3O•, but is inefficient toward peroxyl radicals. In contrast, the radical scavenging efficiency of serotonin, which has been demonstrated in vitro, is significant, and this supports the idea of an indirect antioxidant efficiency of fluoxetine.

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

The interest in diphenyl ditelluride (Ph₂Te₂) is related to its strict analogy to diphenyl diselenide (Ph₂Se₂), 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 Ph₂Se₂ and Ph₂Te₂ 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 Ph₂Te₂ 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 ¹²⁵Te chemical shift.

Nature and strength of chalcogen–π bonds

We have analyzed the chalcogen–π bonding mechanism in a systematic series of model systems through Kohn–Sham molecular orbital theory and a quantitative energy decomposition scheme.