An Evaluation of Various Computational Methods for the Treatment of Organoselenium Compounds

Electrochemical Modification of Polypeptides at Selenocysteine

Mild strategies for the selective modification of peptides and proteins are in demand for applications in therapeutic peptide and protein discovery, and in the study of fundamental biomolecular processes. Herein, we describe the development of an electrochemical selenoetherification (e-SE) platform for the efficient site-selective functionalization of polypeptides. This methodology utilizes the unique reactivity of the 21st amino acid, selenocysteine, to effect formation of valuable bioconjugates through stable selenoether linkages under mild electrochemical conditions. The power of e-SE is highlighted through late-stage C-terminal modification of the FDA-approved cancer drug leuprolide and assembly of a library of anti-HER2 affibody conjugates bearing complex cargoes. Following assembly by e-SE, the utility of functionalized affibodies for in vitro imaging and targeting of HER2 positive breast and lung cancer cell lines is also demonstrated.

Site-selective photocatalytic functionalization of peptides and proteins at selenocysteine

The importance of modified peptides and proteins for applications in drug discovery, and for illuminating biological processes at the molecular level, is fueling a demand for efficient methods that facilitate the precise modification of these biomolecules. Herein, we describe the development of a photocatalytic method for the rapid and efficient dimerization and site-specific functionalization of peptide and protein diselenides. This methodology, dubbed the photocatalytic diselenide contraction, involves irradiation at 450 nm in the presence of an iridium photocatalyst and a phosphine and results in rapid and clean conversion of diselenides to reductively stable selenoethers. A mechanism for this photocatalytic transformation is proposed, which is supported by photoluminescence spectroscopy and density functional theory calculations. The utility of the photocatalytic diselenide contraction transformation is highlighted through the dimerization of selenopeptides, and by the generation of two families of protein conjugates via the site-selective modification of calmodulin containing the 21^st amino acid selenocysteine, and the C-terminal modification of a ubiquitin diselenide.

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.

The one-electron reduction of dithiolate and diselenolate ligands

Herein we present an assessment to determine which of nine well-established DFT functionals best describes the reduction of C2H2Se2(-)˙. In addition, we have also studied the effects of changing the substituents bound to the alkene functional group of dithiolene and diselenolene ligands. Such ligands are important due to their unique electrochemical and physical properties when ligated to metals. The M06-L functional shows best agreement with the QCISD/cc-pVTZ value of -2.45 V for the reduction potential of the (C2H2Se2˙(-)/C2H2Se2(-2)) redox couple. At the M06-L/6-311+G(d,p) level of theory the calculated reduction potential for the (C2H2Se2˙(-)/C2H2Se2(-2)) redox couple is only 0.09 V in error. However, as a result of the nature of the oxidized species for the respective ligands the absolute reduction potential of the (C2H2Se2˙(-)/C2H2Se2(-2)) redox couple is 0.57 V more oxidizing than the (C2H2S2˙(-)/C2H2S2(-2)) redox couple. This is due to the radical electron in C2H2S2˙(-) being delocalized within the alkene backbone, whereas in C2H2Se2˙(-) the electron is largely localized on the Se atoms. The relative reducing power of the S- and Se-containing redox couples is shown to vary depending on the choice of substituents. In particular the reduction potential of the various S-containing redox couples range from being 0.34 V more reducing to 0.28 V more oxidizing than the analogous Se-containing redox couples. This difference in the relative reducing power appears to be a result of the nature of the oxidized ligand. Thus, depending on the choice of moiety very different chemistry is seen between the analogous dithiolate and diselenolate ligands.

Homoselenocysteine — An oxygen or selenium acid in the gas phase?

The potential energy surface of l-homoselenocysteine (HSEC) has been explored through the use of B3LYP/6-311+G(d,p) calculations. In this survey, seventy-seven different conformers have been located. These local minima can be classified in four groups, A–D. Structures A, B, and D are stabilized by intramolecular hydrogen bonds (IMHBs) with the amino group acting as the hydrogen bond (HB) donor and the carbonyl group (structures A and D) or the hydroxyl group (structure B) as HB acceptors. The structures in set C present an IMHB with the amino group acting as the HB acceptor and the hydroxyl group as the HB donor. The stability order decreases in the following order: A > B > C > D. From their relative stabilities it can be concluded that only three of these conformers, namely A1, A4, and A5, would exist in the gas phase at room temperature. The most stable deprotonated form corresponds to a Se-deprotated species stabilized by a strong IMHB between the hydroxyl group and the Se atom. However, a direct deprotonation of the most stable neutrals lead to O-deprotonated species, which eventually isomerize to yield the global minimum. Hence, we can conclude that, quite unexpectedly, HSEC behaves as a Se acid in the gas phase, its intrinsic acidity being 1374 kJ mol–1at the B3LYP/6-311++G(3df,2p) level of theory. The most stable protonated forms are systematically the N-protonated ones, the global minimum being a structure stabilized through an IMHB involving the protonated amino group as the HB donor and the SeH group as the HB acceptor. The calculated gas-phase proton affinity (PA) at the B3LYP/6-311++G(3df,2p) level of theory is 930 kJ mol–1.

Conformationally Locked Bridged Bicyclic Diselenides: Synthesis, Structure, Se···O Interaction, and Theoretical Studies

A general synthetic methodology has been developed for the synthesis of a conformationally locked, bridged diselena-bicyclo[3.2.1]octane skeleton by regio- and stereospecific tandem nucleophilic ring opening of cis-1,4-aziridino-epoxides with tetraethylammonium tetraselenotungstate [Et4N]2WSe4, 1, in a one-pot synthesis. Some correlations have been made on the physicochemical characteristics of the diselenides with a change in the dihedral angles.

Density Functional Theory Study of the Reaction Mechanism and Energetics of the Reduction of Hydrogen Peroxide by Ebselen, Ebselen Diselenide, and Ebselen Selenol

Density functional theory calculations at the B3LYP/6-311++G(3df,3pd)//B3LYP/6-31G(d,p) level have been performed to elucidate the mechanism and reaction energetics for the reduction of hydrogen peroxide by ebselen, ebselen diselenide, ebselen selenol, and their sulfur analogues. The effects of solvation have been included with the CPCM model, and in the case of the selenol anion reaction, diffuse functions were used on heavy atoms for the geometry optimizations and thermochemical calculations. The topology of the electron density in each system was investigated using the quantum theory of atoms in molecules, and a detailed interpretation of the electronic charge and population data as well as the atomic energies is presented. Reaction free energy barriers for the oxidation of ebselen, ebselen diselenide, and ebselen selenol are 36.8, 38.4, and 32.5 kcal/mol, respectively, in good qualitative agreement with experiment. It is demonstrated that the oxidized selenium atom is significantly destabilized in all cases and that the exothermicity of the reactions is attributed to the peroxide oxygen atoms via reduction. The lower barrier to oxidation exhibited by the selenol is largely due to entropic effects in the reactant complex.

Modeling the Reduction of Hydrogen Peroxide by Glutathione Peroxidase Mimics

Theoretical calculations have been performed on three model reactions representing the reduction of hydrogen peroxide by ebselen, ebselen selenol, and ebselen diselenide. The reaction surfaces have been investigated at the B3PW91/6-311G(2df,p) level, and single-point energies were calculated using the 6-311++G(3df,3pd) basis set. Solvent effects were included implicitly with the conductor-like polarizable continuum model and in one case with explicit inclusion of three water molecules. Mechanistic information is gained from investigating the critical points using the quantum theory of atoms in molecules. The barriers for the reduction of hydrogen peroxide with the ebselen, ebselen selenol, and ebselen diselenide models are 56.7, 53.4, and 35.3 kcal/mol, respectively, suggesting that ebselen diselenide may be the most active antioxidant in the ebselen GPx redox pathway. Results are also compared to that of the sulfur analogues of the model compounds.