Addition-Elimination or Nucleophilic Substitution? Understanding the Energy Profiles for the Reaction of Chalcogenolates with Dichalcogenides

Computational analysis of the interactions between Ebselen and derivatives with the active site of the main protease from SARS-CoV-2

The main protease (Mpro) of the novel coronavirus SARS-CoV-2 is a key target for developing antiviral drugs. Ebselen (EbSe) is a selenium-containing compound that has been shown to inhibit Mpro in vitro by forming a covalent bond with the cysteine (Cys) residue in the active site of the enzyme. However, EbSe can also bind to other proteins, like albumin, and low molecular weight compounds that have free thiol groups, such as Cys and glutathione (GSH), which may affect its availability and activity. In this study, we analyzed the Mpro interaction with EbSe, its analogues, and its metabolites with Cys, GSH, and albumin by molecular docking. We also simulated the electronic structure of the generated molecules by density functional theory (DFT) and explored the stability of EbSe and one of its best derivatives, EbSe-2,5-MeClPh, in the catalytic pocket of Mpro through covalent docking and molecular dynamics. Our results show that EbSe and its analogues bound to GSH/albumin have larger distance between the selenium atom of the ligands and the sulfur atom of Cys145 of Mpro than the other compounds. This suggests that EbSe and its GSH/albumin-analogues may have less affinity for the active site of Mpro. EbSe-2,5-MeClPh was found one of the best molecules, and in molecular dynamics simulations, it showed to undergo more conformational changes in the active site of Mpro, in relation to EbSe, which remained stable in the catalytic pocket. Moreover, this study also reveals that all compounds have the potential to interact closely with the active site of Mpro, providing us with a concept of which derivatives may be promising for in vitro analysis in the future. We propose that these compounds are potential covalent inhibitors of Mpro and that organoselenium compounds are molecules that should be studied for their antiviral properties.

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.

Silver-catalyzed direct selanylation of indoles: synthesis and mechanistic insights

Herein we describe the Ag(i)-catalyzed direct selanylation of indoles with diorganoyl diselenides. The reaction gave 3-selanylindoles with high regioselectivity and also allowed direct access to 2-selanylindoles when the C3 position of the indole ring was blocked via a process similar to Plancher rearrangement. Experimental analyses and density functional theory calculations were carried out in order to picture the reaction mechanism. Among the pathways considered (via concerted metalation-deprotonation, Ag(iii), radical, and electrophilic aromatic substitution), our findings support a classic electrophilic aromatic substitution via Lewis adducts between Ag(i) and diorganoyl diselenides. The results also afforded new insights into the interactions between Ag(i) and diorganoyl diselenides.

Reactivity and binding mode of disulfiram, its metabolites, and derivatives in SARS-CoV-2 PLpro: insights from computational chemistry studies

The papain-like protease (PLpro) from SARS-CoV-2 is an important target for the development of antivirals against COVID-19. The safe drug disulfiram (DSF) presents antiviral activity inhibiting PLpro in vitro, and it is under clinical trial studies, indicating to be a promising anti-COVID-19 drug. In this work, we aimed to understand the mechanism of PLpro inhibition by DSF and verify if DSF metabolites and derivatives could be potential inhibitors too. Molecular docking, DFT, and ADMET techniques were applied. The carbamoylation of the active site cysteine residue by DSF metabolite (DETC-MeSO) is kinetically and thermodynamically favorable (ΔG‡ = 3.15 and ΔG =  - 12.10 kcal mol-1, respectively). Our results strongly suggest that the sulfoxide metabolites from DSF are promising covalent inhibitors of PLpro and should be tested in in vitro and in vivo assays to confirm their antiviral action.

In silico studies of Mpro and PLpro from SARS-CoV-2 and a new class of cephalosporin drugs containing 1,2,4-thiadiazole

The SARS-CoV-2 proteases Mpro and PLpro are important targets for the development of antivirals against COVID-19. The functional group 1,2,4-thiadiazole has been indicated to inhibit cysteinyl proteases, such as papain and cathepsins. Of note, the 1,2,4-thiadiazole moiety is found in a new class of cephalosporin FDA-approved antibiotics: ceftaroline fosamil, ceftobiprole, and ceftobiprole medocaril. Here we investigated the interaction of these new antibiotics and their main metabolites with the SARS-CoV-2 proteases by molecular docking, molecular dynamics (MD), and density functional theory (DFT) calculations. Our results indicated the PLpro enzyme as a better in silico target for the new antibacterial cephalosporins. The results with ceftaroline fosamil and the dephosphorylate metabolite compounds should be tested as potential inhibitor of PLpro, Mpro, and SARS-CoV-2 replication in vitro. In addition, the data here reported can help in the design of new potential drugs against COVID-19 by exploiting the S atom reactivity in the 1,2,4-thiadiazole moiety.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11224-022-02036-5.

In the Chalcogenoxide Elimination Panorama: Systematic Insight into a Key Reaction

The selenoxide elimination is a well-known reaction in organochalcogen chemistry, with wide synthetic, biological, and toxicological implications. In this work, we apply benchmarked density functional theory (DFT) calculations to investigate different aspects of the title reaction in three (bio)chemically relevant models, spanning minimal systems of theoretical interests as well as biological or synthetic organochalcogenides. The activation strain analysis (ASA) methodology is employed along a suitable reaction coordinate to obtain insight into the role of the chalcogen and of the oxidation state, to pinpoint the factors that tune the elimination reactivity of the investigated systems. Lastly, we computationally validate the hypothesis that telluroxides eliminate more slowly than selenoxides because of a detrimental hydration process that leads to unreactive hydrates.

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

We estimate the kinetic constants of a series of archetypal S_N2 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 S_N2 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.

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.

Diorganyl diselenides: a powerful tool for the construction of selenium containing scaffolds

Organoselenium compounds find versatile applications in organic synthesis, materials synthesis, and ligand chemistry. Organoselenium heterocycles are widely studied agents with diverse applications in various biological processes. This review highlights the recent progress in the synthesis of selenium heterocycles using diorganyl diselenides with keen attention on green synthetic approaches, scopes, C-H selanylation, the mechanisms of different reactions and insights into the formation of metal complexes. The C-H selanylation using diorganyl diselenides with different catalysts, bases, transition metals, iodine salts, NIS, hypervalent iodine, and other reagents is summarised. Finally, the diverse binding modes of bis(2/4-pyridyl)diselenide with different metal complexes are also summarised.

Mechanistic Insight into SARS-CoV-2 Mpro Inhibition by Organoselenides: The Ebselen Case Study

The main protease (Mpro) of SARS-CoV-2 is a current target for the inhibition of viral replication. Through a combined Docking and Density Functional Theory (DFT) approach, we investigated in-silico the molecular mechanism by which ebselen (IUPAC: 2-phenyl-1,2-benzoselenazol-3-one), the most famous and pharmacologically active organoselenide, inhibits Mpro. For the first time, we report on a mechanistic investigation in an enzyme for the formation of the covalent -S-Se- bond between ebselen and a key enzymatic cysteine. The results highlight the strengths and weaknesses of ebselen and provide hints for a rational drug design of bioorganic selenium-based inhibitors.

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.

The Se…S/N interactions as a possible mechanism of δ-aminolevulinic acid dehydratase enzyme inhibition by organoselenium compounds: A computational study

•

DPDS and PSA interacts with cysteine residues from AlaD active site.

•

The Se…S interactions could be involved in the δ-AlaD inhibition.

•

δ-AlaD from Cucumis sativus does not present cysteine residues in the active site.

•

Se…N interactions could be involved in the organoselenium action.

Organoselenium compounds present many pharmacological properties and are promising drugs. However, toxicological effects associated with inhibition of thiol-containing enzymes, such as the δ-aminolevulinic acid dehydratase (δ-AlaD), have been described. The molecular mechanism(s) by which they inhibit thiol-containing enzymes at the atomic level, is still not well known. The use of computational methods to understand the physical–chemical properties and biological activity of chemicals is essential to the rational design of new drugs. In this work, we propose an in silico study to understand the δ-AlaD inhibition mechanism by diphenyl diselenide (DPDS) and its putative metabolite, phenylseleninic acid (PSA), using δ-AlaD enzymes from Homo sapiens (Hsδ-AlaD), Drosophila melanogaster (Dmδ-AlaD) and Cucumis sativus (Csδ-AlaD). Protein modeling homology, molecular docking, and DFT calculations are combined in this study. According to the molecular docking, DPDS and PSA might bind in the Hsδ-AlaD and Dmδ-AlaD active sites interacting with the cysteine residues by Se^…S interactions. On the other hand, the DPDS does not access the active site of the Csδ-AlaD (a non-thiol protein), while the PSA interacts with the amino acids residues from the active site, such as the Lys291. These interactions might lead to the formation of a covalent bond, and consequently, to the enzyme inhibition. In fact, DFT calculations (mPW1PW91/def2TZVP) demonstrated that the selenylamide bond formation is energetically favored. The in silico data showed here are in accordance with previous experimental studies, and help us to understand the reactivity and biological activity of organoselenium compounds.

Silicon-Free SuFEx Reactions of Sulfonimidoyl Fluorides: Scope, Enantioselectivity, and Mechanism

SuFEx reactions, in which an S−F moiety reacts with a silyl‐protected phenol, have been developed as powerful click reactions. In the current paper we open up the potential of SuFEx reactions as enantioselective reactions, analyze the role of Si and outline the mechanism of this reaction. As a result, fast, high‐yielding, “Si‐free” and enantiospecific SuFEx reactions of sulfonimidoyl fluorides have been developed, and their mechanism shown, by both experimental and theoretical methods, to yield chiral products.

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.

Towards New Catalytic Antioxidants: A Simple and Mild Synthesis of Selenenylsulfides

A new methodology for the synthesis of small molecules containing the S-Se bond is reported. Aryl- and alkyl-selenols react smoothly with N-thiophthalimides to afford the corresponding selenenylsulfides through a clean SN2 path occurring at the sulfur atom. The reaction proceeds under very mild conditions in DMF in absence of catalysts for most of the substrates. The scope of the reaction was found to be broad, allowing a wide series of selenols and N-thiophtalimides to be efficiently employed in this procedure. Owing to the instability of the S-Se bond, selenenylsulfides exhibited a remarkable tendency to disproportionate to the corresponding symmetric diselenides and disulfides. Preliminary evaluation of the catalytic antioxidant properties of novel selenenylsulfides showed their behaviour as GPx mimics.

Psychiatric Disorders and Oxidative Injury: Antioxidant Effects of Zolpidem Therapy disclosed In Silico

Zolpidem (N,N-Dimethyl-2-[6-methyl-2-(4-methylphenyl)imidazo[1,2-a]pyridin-3-yl]acetamide) is a well-known drug for the treatment of sleeping disorders. Recent literature reports on positive effects of zolpidem therapy on improving renal damage after cisplatin and on reducing akinesia without sleep induction. This has been ascribed to the antioxidant and neuroprotective capacity of this molecule, and tentatively explained according to a generic structural similarity between zolpidem and melatonin. In this work, we investigate in silico the antioxidant potential of zolpidem as scavenger of five ROSs, acting via hydrogen atom transfer (HAT) mechanism; computational methodologies based on density functional theory are employed. For completeness, the analysis is extended to six metabolites. Thermodynamic and kinetic results disclose that indeed zolpidem is an efficient radical scavenger, similarly to melatonin and Trolox, supporting the biomedical evidence that the antioxidant potential of zolpidem therapy may have a beneficial effect against oxidative injury, which is emerging as an important etiopathogenesis in numerous severe diseases, including psychiatric disorders.

Nucleophilic Substitution (SN 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent

The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.

Antioxidant activity of selenenamide-based mimic as a function of the aromatic thiols nucleophilicity, a DFT-SAPE model

The mechanism of action of the selenenamide 1 as a mimic of the glutathione peroxidase (GPx) was investigated by the density functional theory. The solvent-assisted proton exchange procedure was applied to model the catalytic behavior and antioxidant activity of this mimic. To have an insight into the charge transfer effect, different aromatic thiols, including electron donating substituents on the phenyl ring were considered. The catalytic behavior of the selenenamide was modeled in a four-step mechanism, described by the oxidation of the mimic, the reduction of the obtained product, selenoxide, the reduction of the selenenylsulfide and dehydration of selenenic acid. On the basis of the activation parameters, the final step of the proposed mechanism is the rate determining states of the catalytic cycle. Turnover frequency (TOF) analysis showed that the electron donating groups at the para-position of the phenyl ring of the PhSH do not affect the catalytic activity of the selenenamide in contrast to p-methyl thiophenol which indicates the highest nucleophilicity. The evaluation of the electronic contribution of the various donating groups on the phenyl ring of the aromatic thiols shows that the antioxidant activity of the selenenamide sufficiently increases in the presence of the electron-donating substitutions. Finally, the charge transfer process at the rate-determining state was investigated based on the natural bond orbital analysis.

Mechanistic Insight into the Oxidation of Organic Phenylselenides by H2 O2

The oxidation of organic phenylselenides by H₂ O₂ 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 (¹ H and ⁷⁷ 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.

Computational Modeling of the Catalytic Cycle of Glutathione Peroxidase Nanomimic

To elucidate the role of a derivative of ebselen as a mimic of the antioxidant selenoenzyme glutathione peroxidase, density functional theory and solvent-assisted proton exchange (SAPE) were applied to model the reaction mechanism in a catalytic cycle. This mimic plays the role of glutathione peroxidase through a four-step catalytic cycle. The first step is described as the oxidation of 1 in the presence of hydrogen peroxide, while selenoxide is reduced by methanthiol at the second step. In the third step of the reaction, the reduction of selenenylsulfide occurs by methanthiol, and the selenenic acid is dehydrated at the final step. Based on the kinetic parameters, step 4 is the rate-determining step (RDS) of the reaction. The bond strength of the atoms involved in the RDS is discussed with the quantum theory of atoms in molecules (QTAIM). Low value of electron density, ρ(r), and positive Laplacian values are the evidence for the covalent nature of the hydrogen bonds rupture (O₃₀-H₃₁, O₃₃-H₃₄). A change in the sign of the Laplacian, L(r), from the positive value in the reactant to a negative character at the transition state indicates the depletion of the charge density, confirming the N₅-H₁₀ and O₁₁-Se₁ bond breaking. The analysis of electron location function (ELF) and localized orbital locator (LOL) of the Se₁-N₅ and Se₁-O₁₁ bonds have been done by multi-WFN program. High values of ELF and LOL at the transition state regions between the Se, N, and O atoms display the bond formation. Finally, the main donor-acceptor interaction energies were analyzed using the natural bond orbital analysis for investigation of their stabilization effects on the critical bonds at the RDS.