Modeling the Oxidation of Ebselen and Other Organoselenium Compounds Using Explicit Solvent Networks

How solvent determines the molecular reactive conformation and the selectivity: Solvation spheres and energy

In solution, the solvent determines the molecular conformation and the chemical reaction viability and selectivity. When solvent-solute and solvent-solvent interactions present similar strengths, explicit salvation is the best way to describe a system. The problem to solve is how big the explicit shell should be. In this paper, we want to answer one of the fundamental questions in the implementation of explicit solvation, exactly how many solvent molecules should be added and where they should be placed. Here we determine the first solvent sphere around a molecule and describe how it controls the conformation and selectivity of a selected reaction. NMR experiments were carried out to identify the number of solvent molecules around the solute that constitutes the first solvent sphere, and the interaction between this solvent sphere and the solute was detected using DFT and QTAIM calculations. A new approach to the solvation energy is presented. Finally, we established the role of solvent molecules in the conformation of the solute and in the transition states that produce the two possible products of the reaction.

Topological and quantum molecular descriptors as effective tools for analyzing cytotoxic activity achieved by a series of (diselanediyldibenzene-4,1-diylnide)biscarbamate derivatives

A molecular modeling study has been carried out on a previously reported series of (diselanediyldibenzene-4,1-diylnide)biscarbamate derivatives that show cytotoxic and antiproliferative in vitro activity against MCF-7 human cell line; radical scavenging properties were also confirmed when these compounds were tested for their ability to scavenge DPPH and ABTS radicals. The data obtained allowed us to classify the compounds into two different groups: (a) aliphatic carbamates for which the activity could be related with a first nucleophilic attack (mediated by H₂O, for example) on the selenium atoms of the central scaffold, followed by the release of the alkyl N-(4-selanylphenyl) and N-(4-selenenophenyl)carbamate moieties. Then, a second nucleophilic attack on the carbamate moiety, to yield 4-aminobenzeneselenol and 4-selenenoaniline respectively, which can ultimately be responsible for the activity of the compounds; (b) aromatic carbamates, for which we propose a preferred nucleophilic attack on the carbamate moiety, yielding 4-[(4-aminophenyl)diselanyl]aniline, the common structural fragment for this series, for which we have previously demonstrated its cytotoxic profile. Then, selenium atoms of the central fragment may later undergo a new nucleophilic attack, to yield 4-selenenoaniline and 4-aminobenzeneselenol. The phenolic moieties released in this process may also have a synergistic cytotoxic and redox activity. The data that support this connection include the conformational behavior and the molecular topography of the derivatives which can influence the accessibility of the hydrolysis points, and some quantum descriptors (bond order, atomic charges, total valences, ionization potential, electron affinity, HOMO 0 and LUMO 0 location, etc.) that have been related to the biological activity of the compounds.

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.

Computational kinetic modeling of the selenol catalytic activity as the glutathione peroxidase nanomimic

Density functional theory and solvent-assisted proton exchange methods have been applied for computational modeling of the catalytic cycle of selenol zwitterion anion from the kinetic and thermodynamic viewpoints. Selenol zwitterion anion has been represented as an effective glutathione peroxidase nanomimic. It reduces peroxides through a three-step pathway. In the first step, seleninic acid is produced through deprotonating of the selenol zwitterion anion in the presence of the hydrogen peroxide. Seleninic acid reacts with a thiol to form selenylsulfide in the second step. In the last step, selenylsulfide is reduced by the second thiol and regenerates selenolate anion through disulfide formation. Selenol zwitterion anion in comparison to more widely studied compounds such as ebselen has a good activity to react with hydrogen peroxide and producing seleninic acid. The energy barrier of this reaction is 11.7kcalmol^-1 which is smaller than the reported enzyme mimics. Moreover, the reactions of seleninic acid and selenylsulfide with methanethiol, which is used as a nucleophile, are exothermic by -18.4 or -57.0kcalmol^-1, respectively. Based on the global electron density transfer value of -0.507 e from the natural atomic charge analysis, an electronic charge depletion at the transition state (TS), electron-donor substitutions on the selenolate facilitates the reduction reaction, effectively. Finally, the nature of the bond formation/cleavage at the TS has been quantitatively described by using the topological analyses.

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

We have quantum chemically explored the mechanism of the substitution reaction between CH3X(-) and the homo- and heterodichalcogenides CH3X'X″CH3 (X, X', X″ = S, Se, Te) using relativistic density functional theory at ZORA-OLYP/TZ2P and COSMO for simulating the effect of aqueous solvation. In the gas phase, all substitution reactions proceed via a triple-well addition-elimination mechanism that involves a stable three-center intermediate. Aqueous solvation, in some cases, switches the character of the mechanism to double-well SN2 in which the stable three-center intermediate has become a labile transition state. We rationalize reactivity trends and some puzzling aspects of these elementary reactions, in particular, vanishing activation energies and ghost three-center intermediates, using the activation strain model (ASM).

Organochalcogen peroxidase mimetics as potential drugs: a long story of a promise still unfulfilled

Organochalcogen compounds have attracted the interest of a multitude of studies to design potential therapeutic agents mimicking the peroxidase activity of selenium-based glutathione peroxidases (GPx's). Starting from the pioneering ebselen, various compounds have been synthesized over the years, which may be traced in three major classes of molecules: cyclic selenenyl amides, diaryl diselenides, and aromatic or aliphatic monoselenides. These compounds share common features and determinants needed to exert an efficient GPx-like activity, such as polarizing groups in close proximity to selenium and steric effects. Nonetheless, the reactivity of selenium, and tellurium as well, poses serious problems for the predictability of the biological effects of these compounds in vivo when used as potential drugs. These molecules, indeed, interfere with thiols of redox-regulated proteins and enzymes, leading to unexpected biological effects. The various chemical aspects of the reaction mechanism of peroxidase mimetics are surveyed here, focusing on experimental evidence and quantum mechanics calculations of organochalcogen representatives of the various classes.

Transition states for cysteine redox processes modeled by DFT and solvent-assisted proton exchange

Oxidation and disulfide coupling of cysteine, processes central to oxidative stress and biochemical signaling, are modeled using DFT and solvent-assisted proton exchange, a method of microsolvation. Calculated barriers are consistent with experimental kinetics and observed product ratios and suggest a dependence on the polarity of the surrounding medium.