Computational Kinetic Modeling of the Catalytic Cycle of Glutathione

Dec 7, 2017 - The catalytic cycle of a new derivative of ebselen, 1, was elucidated via three steps by the density functional theory and solvent-assis...
0 downloads 0 Views 2MB Size
Subscriber access provided by READING UNIV

Article

Computational Kinetic Modeling of the Catalytic Cycle of Glutathione Peroxidase Nanomimic: Effect of Nucleophilicity of Thiols on the Catalytic Activity Ramesh Kheirabadi, Mohammad Izadyar, and Mohammad Reza Housaindokht J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b09929 • Publication Date (Web): 07 Dec 2017 Downloaded from http://pubs.acs.org on December 9, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Computational Kinetic Modeling of the Catalytic Cycle of Glutathione Peroxidase Nanomimic: Effect of Nucleophilicity of Thiols on the Catalytic Activity Ramesh Kheirabadia, Mohammad Izadyarb*, Mohammad Reza Housaindokhtb a

Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, International Campus, Mashhad, Iran b

Department of Chemistry, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, Iran [email protected] Telefax: ++985138795457

ABSTRACT The catalytic cycle of a new derivative of ebselen, 1, was elucidated via three steps by the density functional theory and solvent-assisted proton exchange procedure involving indirect proton exchange through a hydrogen-bonded transfer network. Different behaviors of the aromatic and aliphatic thiols were investigated in the reduction of selenoxide (step 2→3) and selenurane (step 3→1) based on their nucleophilicity. The reduction of selenoxide in the presence of thiophenol (ΔG≠=15.9 kcal.mol-1) is faster than methanethiol (ΔG≠=29.3 kcal.mol-1), and methanethiol makes the reduction of selenoxide unspontaneous and kinetically unfavorable (ΔG=2.8 kcal.mol-1). The nucleophilic attack may be enhanced by using the thiophenol backbone at the selenium center to lower the energy barrier of the selenoxide reduction (ΔG≠=15.9 kcal.mol-1). On the 1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 44

basis of the turnover frequency calculations, during the catalytic cycle, the rate of the reaction was analyzed and discussed. Low values of the electron density and Laplacian at the transition states are the evidences of the covalent O-H and O-O bonds rupture the in the presence of methanethiol and thiophenol. The nature of the critical bond points was characterized, using the quantum theory of atoms in molecules, based on the electron location function and localized orbital locator values. Finally, the charge transfer process at the rate-determining step was investigated based on the natural bond orbital analysis. 1- INTRODUCTION In the last decades, the biological consistency of the characterized element, selenium, has been considerably designated in contrast with the toxic and beneficial effects.1 Although selenium plays the essential roles in the animal species and many bacteria,2 in 2011, selenazolidine and its relation with KashinBeck disease was studied in order to demonstrate that the selenium-containing proteins are composed of selenium as a beneficial trace element for humans and mammalian.3 In fact, selenium-enriched yeast and garlic were discovered as a crucial ingredient in our diet and inorganic selenium compounds, such as sodium 2 ACS Paragon Plus Environment

Page 3 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

selenite, restore the GPx activity.4 Many organoselenium compounds have been prepared to be used as antioxidant, anti-infective and antitumor agents5 and also enzyme inhibitors such as cytokine inducers.6 The first distinguished selenoprotein GPx1 is composed of five amino acids to protect their environment, biologically, against the oxidative stress through a catalytic process through the peroxides reduction. However, it mitigates the toxicity of the hydroperoxides through the reduction of hydrogen peroxide to water or similar compounds via a three-step catalytic mechanism.7 The commonly accepted catalytic cycle of the GPx is illustrated in Scheme 1.8 In the first step, the selenol part of the GPx is oxidized to selenenic acid and the peroxide is reduced to water. GSH is converted to the selenenylsulfide and H2O in the second step. In the third step, thioldisulfide, GSSG, is produced through the reduction of GSH, followed by the regeneration of selenol enzyme. Since GPx shows an original antioxidant activity in vivo, it cannot be substituted by other selenoproteins.9 The fundamental role of the mineral selenium and the selenocysteine incorporated in GPx during the translation has encouraged the researchers to synthesize numerous small molecules with GPx-like antioxidant

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 44

activity.10 These mimics obviously point out to the catalytic properties, substrate-specific and low toxicity.

Scheme 1. Catalytic mechanism of the GPx enzyme

The new generation of the GPx mimics are classified into the three main groups: cyclic selenenylamides (Ebselen: 2-phenyl-1,2-benzisoselenazol3(2H)-one), organodiselenides (DMBS: N,N-dimethylbenzylamine-2-selenol ) and aryl/alkyl monoselenides which are depicted in Scheme 2.11-14 The first proposed nontoxic scavenger of the reactive oxygen and nitrogen species is ebselen, which has been synthesized more than 30 years ago.15 In 2011, the mechanism of ebselen as a GPx-mimic in the catalytic cycle was studied, theoretically, by the proposed density functional theory-solvent assisted 4 ACS Paragon Plus Environment

Page 5 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

proton exchange (DFT-SAPE) model.16 On the basis of the obtained results, it was concluded that a realistic interpretation of the redox mechanism of ebselen is provided via the SAPE modeling.

Scheme 2. Chemical structures of the new generations of the GPx mimics

Explicit water molecules provide an indirect pathway for proton shuttle, showing the role of bulk water in the solution phase as the proton transfer procedure.17 One of these processes is the microsolvation technique as the solvent assisted proton exchange model which designs the explicit solvation effects.18 For example, in a research, the oxidation reaction of the ebselen was 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 44

considered by the SAPE-DFT modeling of the oxygen atom transfer to show an acceptable rate constant in comparison with the experimental results.19 The reduction mechanism of GPx mimic oxidation with thiols was studied by Heverly-Coulson and Boyd in 2012.12 According to their research, DMBS as a GPx mimic reduces the peroxides via a three-step catalytic mechanism in the presence of the thiol and thiolate. Thiol was considered in the proposed mechanism as the nucleophile and two molecules of water were applied as the solvent in the proton shuttle process. On the basis of the obtained results, neither methanethiol nor thiophenol has a significant effect on the energy barriers. Bayse and Ortwine, theoretically modeled the GPx-like selectivity of the seleninate by the DFT–SAPE approach in 2013.20 In this research, Seleninate which is produced by oxidation of the ω-hydroxyalkyl selenides, is more efficient than ebselen. According to the proposed models, Seleninate as a catalyst is regenerated by the selenenylsulfide.21-22 In order to analyze the mimetic aspects of the seleninate, its performance was examined via the ROS scavenging mechanism. For this purpose, several catalytic pathways were designed to determine the most probable intermediates and transition states

6 ACS Paragon Plus Environment

Page 7 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(TSs). On the basis of the DFT-SAPE results, its catalytic activity was confirmed.23 In 2015, the oxidation of the selenocysteine was investigated in the GPx catalytic cycle.7 The active site was theoretically modeled to validate the mechanism aspects and the efficiency of the GPx mimic. The destruction of the redox center is prohibited by the selenenylamide route due to selenium elimination or over-oxidation. In the previous work, we have employed the quantum chemistry methods as an important tool to inquire the catalytic activity of the new GPx mimic (selenenamide) which was synthesized from the ebselen and related compounds. 24

DFT-SAPE model, natural bond orbital (NBO) and quantum theory of atoms

in molecules (QTAIM) analyses were applied to demonstrate the indirect effect of the hydrogen transfer and the nature of bond formation and cleavage, respectively. Furthermore, the kinetic and mechanistic aspects of the proposed catalytic cycle of the GPx mimic and rate-determining step were investigated through a four-step mechanism. Considering the most recent researches on the enzyme mimics shows a great interest to the reactions of the amid-based organoselenium compounds.25 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 44

Although DMBS exhibited a poor antioxidant and low mimetic activities,12 it is possible to design new molecular structures containing both the Se−N covalent bond and the amine group. In this work, we consider the effects of different thiols as the substitute on the catalytic activity of new GPx nano mimic 1 (Figure 1). For this purpose, the molecular mechanisms were calculated by the DFT method from the energetic point of view to characterize the intermediates and TSs on the reaction pathways, relevant mechanism steps and rate-determining step (RDS).

2- COMPUTATIONAL DETAILS Geometry optimization and vibrational frequency calculations were performed by using the Gaussian 0926 at the M06-2X/6-31+G(d,p) level of theory which is an appropriate functional for organoselenium systems.27-29 The structures of the TSs were located by the Schlegel’s synchronous transit-guided quasi-Newton method (STQN) and their accuracy was checked via the intrinsic reaction coordinate (IRC) calculations.30,31 Thermodynamic parameters were obtained at the M06-2X/ 6-31+G (d,p) level within the zero-point vibrational energy, thermal and entropy corrections 8 ACS Paragon Plus Environment

Page 9 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

which are based on the harmonic oscillator approximation and solvation corrections at 298.15 K and 1 atm. Microsolvation technique as the solvent assisted proton exchange network was distinguished for designing the explicit solvation effects by two or three water molecules (dielectric constant of 78.36) which were applied in the conductor-like polarizable continuum model (CPCM).32-33 Additionally, there may be favorable free energies contribution due to an increased entropy which is modeled by the explicit solvent, because the entropy of the solvent is difficult to be modeled, implicitly.34. To include more accurate theoretical data, single point energies were obtained for all of the studied structures by using the mPW1PW91/6311++G(d,p) level of theory.35-36 The mPW1PW91 functional provides the best agreement with the several tested post-HF methods such as MP2 and this functional is accurate in the benchmarks of the thermochemical and kinetic parameters of the experimental data.37 On the basis of the relative energies of the catalytic cycle, state energies (E-representation) were generated by the computational study instead of the rate constants (k-representation) in an experimental study. The “state energy” is 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 44

defined by the turnover frequency (TOF). TOF is evaluated by the energetic span model which are used to demonstrate several catalytic cycles of the organometallic reactions. Accurate TOF (turnover frequency) was evaluated based on equation 1:38 TOF =

KB T h

e

−δE⁄ RT

(1)

Where kB and h are the Boltzmann and Planck constants, respectively. R is the gas constant (1.98 cal. mol-1. K-1) and T is the absolute temperature (298.15 K). δE is the energetic span as the energy difference between the TDI (TOFdetermining intermediate) and the TDTS (TOF-determining transition state) which is calculated by equation 2.39 TTDTS - ITDI δE =

if TDTS appears after TDI (2)

TTDTS - ITDI + ΔGt

if TDTS appears before TDI

TTDTS is the energy at the highest transition state, ITDI means the energy at the lowest intermediate state and ΔGt is the reaction Gibbs energy during the catalytic cycle. On the basis of the QTAIM analysis, the topological properties such as electron densities, kinetic and potential energy densities by AIM2000

10 ACS Paragon Plus Environment

Page 11 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

software package,40-41 electron location function (ELF)42-47 and localized orbital locator (LOL)48-49 were calculated by Multi WFN 3.1.50 Also, NBO method was applied to evaluate the electronic charge distribution into the atomic and molecular orbitals.51 Additionally, the secondorder perturbation energies (E2) of the donor-acceptor interactions at the RDS of the methanethiol and thiophenol were calculated to investigate the stability of the structures during the reaction proceeding. E(2) = qi

(Fij ) 2

(3)

εj − εi

Where qi is the orbital occupancy, εi, εj show the diagonal elements and Fi,j represents the off-diagonal NBO Fock matrix element.52 3- RESULTS AND DISCUSSION The catalytic activity of mimic 1 in the presence of various thiols is depicted in Figure 1. To investigate the influence of the aliphatic and aromatic substituent effects on the catalytic activity of the mimic 1, methanethiol and thiophenol were considered. According to Figure 1, mimic 1 is regenerated via two intermediates of selenoxide 2 and selenurane 3 which were trapped in vitro.53 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 44

Figure 1. Proposed catalytic cycle for the new GPx mimic 1, R = Ph, CH3.

The proposed mechanism, Figure 1, is composed of three main steps of the oxidation/reduction processes: the oxidation/reduction of 1 (step 1), the intermediate conversion (step 2), and the regeneration of the mimic 1 (step 3). Two and three water molecules (SAPE microsolvent model) are used to build a stable indirect proton exchange network as the hydrogen-bonded network in the proposed catalytic cycle. This technique makes a synchronous process of proton shuttle because water molecules play a key role in the limitation of heavy atom linkage/breakage. 12 ACS Paragon Plus Environment

Page 13 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

- Oxidation of 1 (Step 1 → 2). Although the oxidation reaction of ebselen was considered via a two-step hydrogen transfer mechanism,54 the oxidation of selenol zwitterion was investigated by the protons via a single-step process.55 On the basis of the obtained data, a single-step hydrogen transfer via SN2 reaction was considered as the most probable path to complete the oxidation of 1. In the experimental studies, it was suggested that the selenocysteine residue in the GPx enzyme, directly enters into the catalytic process of hydroperoxide reduction .56-58 In the first step, the reaction of 1 with H2O2, the oxidative stress conditions, has been designed by Fisher and Dereu59 to acquire selenoxide 2. In hydrogen peroxide structure, the O-O bond is broken to alleviate the formation of the Se=O bond in the selenoxide. It is noteworthy that, in the first step, the peroxide bond is broken which is followed by the proton transfer, yielding the O....H bond. The calculated bond lengths of the Se=O are 3.36Å, 2.17Å, 1.77Å in the reactant, TS and product, respectively, confirming the formation of selenoxide 2. The Gibbs energy change of this step (Figure 2a, Figure 3a) shows an 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 44

exothermic nature, -44.1 kcal·mol-1, and the activation energy of the oxidation of 1 is 26.8 kcal.mol-1. Although the calculated activation energy is higher than the experimental data for the GPx (14.9 kcal.mol-1)60, it significantly overestimates the reaction barrier of the GPx1 more than 10.0 kcal mol-1, but the activation energy of the oxidation of 1 (26.8 kcal.mol-1) is in good agreement with the hybrid QM/MM calculations on the first redox step of the Bovine GPx (26.3 kcal.mol-1).27 - Reduction of selenoxide 2 in the presence of thiols (Step 2 → 3). It is well worth mentioning that selenylsulfide formation (E-Se-SG) through the seleno-sulfide linkage has been experimentally observed.61 Hence, on the basis of the experimental data, the reduction mechanism of selenoxide 2 was modeled in the presence of thiols. The H-transfer process from the thiol proceeds through a concerted four-center TS, in which three-hydrogen atom transfer is related to: 1) From the thiol to water molecule, 2) from the water molecule to the oxygen atom of the second water molecule, 3) from the second water molecule to the oxygen atom of selenoxide which is followed by the Se−S bond formation, yielding selenurane and also three water molecules. This model of proton shuttle 14 ACS Paragon Plus Environment

Page 15 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

is possible through the incorporation of the protic solvent that plays an important role in the thiolic hydrogen atom transfer to the oxygen atom of selenoxide. On the basis of the computational studies of Coulson and Boyd, two water molecules are necessary for hydrogen transfers.62 An effective process of the proton shuttle is not possible through a single water molecule, because a single water molecule is not sufficiently flexible to accept a hydrogen bond from the thiol. Accordingly, adding a second water molecule produces a chain in the reactant complex, which accepts a hydrogen bond from the thiol and donates it to the hydroxyl group. The optimized structures of the stationary points on the reaction paths are depicted in Figures 2 and 3 in the presence of methanethiol and thiophenol, respectively. In these figures, the reaction of the selenoxide 2 with the methanethiol and thiophenol has been considered in some parts of the Figures 2b and 3b, respectively. On the basis of the proposed mechanism, the hydrogen atom of the thiols is transferred to the oxygen atom of selenoxide 2 in the presence of the explicit water molecules to produce selenurane and three water molecules, which are necessary and sufficient for a normal chain of the proton transfer. 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

Page 16 of 44

Figure 2. Optimized stationary states on the reaction pathways within the imaginary vibrational frequencies (in parentheses) of the TSs in the presence of methanethiol. 16 ACS Paragon Plus Environment

Page 17 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

The Journal of Physical Chemistry

Figure 3. Optimized stationary states on the reaction pathways within the imaginary vibrational frequencies (in parentheses) of the TSs in the presence of thiophenol. 17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 44

This behavior can be explained based on the research reported by Wang and coworkers on the Se-N bond cleavage.63 They investigated the effects of the implicit water as the solvent on the Se-N bond splitting and concluded that the sulfur atom of the thiols is bonded to the selenium atom, without any hydrogen transfer. Since the hydrogen transfer is an important step in the proposed mechanism, using the explicit solvent molecules is not only necessary to describe comprehensively the H-transfer path, but also is suitable for completion of the catalytic cycle.

- Reduction of 3 (Selenurane) by thiols (Step 3 → 1). The final step of the catalytic cycle of 1 as the GPx mimic is considered by the reduction of thiols with selenurane in a concerted manner. Thiol acts as a nucleophile and hydrogen transfer occurs through a chain of three water molecules from the thiol to the hydroxyl group of selenium as a leaving group at the TS 31. The sulfur atom of the selenurane is linked to the sulfur atom of the thiol via the collinear interaction which is followed by the production of disulfide and water (see Figures 2c and 3c). The regeneration of 1 through the reduction of selenurane with methanethiol is 2.5 kcal.mol-1 more stable than 18 ACS Paragon Plus Environment

Page 19 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

thiophenol (Table 1). Furthermore, the nucleophilic attack to the sulfur atom in selenurane is an exothermic reaction. Calculated Gibbs energies of -29.3 kcal.mol-1 and -26.8 kcal.mol-1 by methanethiol and thiophenol, respectively, show that this step of total reaction (3 →1) is thermodynamically favorable. The comparison of the S-S bond length in disulfide formation in the presence of methanethiol is 1.4 Å less than thiophenol, confirming that step 3→1 is more probable in methanethiol.

- Structural and Energy Analysis. Calculated thermodynamic parameters for three steps of the catalytic cycle of 1 are represented in Table 1. According to Table 1, Gibbs energy for the reduction of hydrogen peroxide in the proposed mimic 1 is -44.1 kcal.mol-1 with an activation Gibbs energy of 26.8 kcal.mol-1. A positive value of ΔG (2→3) in the presence of methanethiol indicates that this step of the proposed mechanism is not spontaneous. Since this value is on the border line of the spontaneous and non-spontaneous reaction, it may be changed by an increment in the applied theoretical level of the calculation.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 44

Table 1. Calculated thermodynamic and activation parameters (kcal.mol-1) for the studied step of the proposed catalytic cycle of mimic 1. Reaction

Methanethiol

Thiophenol

pathway

ΔH≠

ΔH

ΔG≠

ΔG

ΔE≠

ΔE*

ΔH≠

ΔH

ΔG≠

ΔG

ΔE≠

ΔE*

1→2

25.5

-44.6

26.8

-44.1

18.4

-43.3

25.5

-44.6

26.8

-44.1

18.4

-43.3

2→3

26.6

1.8

29.3

2.8

29.6

-0.5

12.7

-1.2

15.9

-0.5

8.7

-4.6

3→1

4.5

-26.5

5.5

-29.3

9.3

-28.5

5.8

-21.6

8.9

-26.8

9.7

-26.1

*Single point energies were calculated by mPW1PW91/6-311++G(d,p). A lower barrier in the presence of thiophenol is likely due to its lower pKa, which enhances the reactivity. Also, the phenol ring of thiophenol stabilizes the TS23 through the resonance effects which is absent in the case of methanthiol. Therefore, the nucleophilic attack may be enhanced by using the thiophenol (phSH) backbone at the selenium center that decreases the lower energy barrier of the selenoxide reduction. To obtain a more accuracy in the theoretical data, high level of the computational theory and larger 6-311++G(d,p) basis set were chosen, based on the fact that all mechanisms in the catalytic cycle of the mimic 1 are spontaneous in the presence of the thiols, which is according to the experimental data. The optimized structures were reconsidered from the energy viewpoint at the 20 ACS Paragon Plus Environment

Page 21 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

mPW1PW91/6-311++G(d,p) level. According to Table 1, a value of -0.5 kcal.mol-1 as the reaction energy, at this higher level of theory, shows that this step of the proposed mechanism is thermodynamically favorable in the presence of methanethiol. Theoretical trends of the thermodynamic stability of the studied products in three steps are as follows: (1→2) ˃ (3→1) ˃ (2→3). Kinetic studies on the reaction paths show that the highest activation Gibbs energy in the presence of methanethiol is related to step 2→3, having a value of 29.3 kcal.mol-1, which shows that this step is the RDS, in agreement with the experimental data.65 There is a similar discussion for the reaction in the presence of thiophenol. In this case, step 1→2 is the RDS with an activation energy of 26.8 kcal.mol -1, which is not according to methanethiol. Kinetic studies shows that thiophenol, with a lower activation energy at the RDS, accelerates the catalytic reaction of 1 more than methanethiol. Moreover, a reported lower pKa for the thiophenol makes it proper leaving groups in these types of reactions.66 Relative Gibbs energy changes as a function of the reaction progressin the presence of the methanthiol and thiophenol is

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 44

depctide in Figure 4 as the potential energy diagram (PED). On the basis of the PED, it is confirmed that the most stable intermediate is selenoxide (P2). Since methanethiol and thiophenol are involved in step 2→3 of the catalytic mechanism, the obtained activation energies in step 1→2 are equal which is in contrast to steps 2→3 and 3→1. Moreover, TS23 and TS12 are the most important critical structures on the PED, showing the RDS of the reaction in two cases.

Figure 4. Potential energy diagram (PED) of the proposed mechanism; the reactants Gibbs energy in each step has been considered in the presence of thiols as the reference, for more clarity.

-Accurate TOF Calculation 22 ACS Paragon Plus Environment

Page 23 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Turnover frequency is used in the theoretical studies to make a relation between the theoretically calculated energy representation (E-representation) and the rate constant (k-representation) in the experimental studies.67 TOFdetermining transition state and the TOF-determining intermediate have been innovated in the catalytic cycles by equation 1. Meanwhile, δE is the energetic span to make the kinetic assessment possible through a specific computationally calculated catalytic cycle. In Figure 5, TDTS and TDI are illustrated by the energy profile of the modeled catalytic cycle of 1 in the presence of the methanthiol and thiolphenol.

Figure 5. Energy diagram of the proposed catalytic cycle to demonstrate the TOF transition states (TDTS) and TOF determining intermediate (TDI) in the presence of thiols.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 44

Figure 5 shows that in the presence of the methanethiol, TDTS appears after the TDI, which means that δE value (δE= TTDTS – ITDI) is 72.9 kcal.mol-1 and the corresponding TOF is equal to 3.7 × 1013 s-1. An opposite behavior is observed in the case of thiophenol, TDI appears after TDTS, which means that δE (δE= TTDTS – ITDI + ΔGt) is 35.6 kcal.mol-1 and TOF= 3.9 × 1013 s-1. Hawkes and coworkers indicated that a smaller energetic span (δE) showed a faster reaction.68,69 As a result, the catalytic cycle of the mimic 1 in the presence of thiophenol is faster than methanethiol, according to the previous analyses. - Topological Analysis The proton transfer process between the methanethiol and selenoxide in the reduction of selenoxide (2→3) and the oxidation of 1 (1→2) in the presence of thiophenol, which are introduced as the rate determining step, were studied in the framework of the QTAIM. Electron density, ρ(r), and Laplacian of the electron density, L(r), were evaluated at the bond critical points (BCPs). Additionally, the bond ellipticity at the BCPs, 𝜀 = (𝜆1/𝜆2) – 1, which estimated the anisotropy of the curvature of the electron density in the direction perpendicular to the bond, was calculated. This molecular descriptor is sensitive to the bond cleavage and formation.70,71 In the mathematical equation of the 24 ACS Paragon Plus Environment

Page 25 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ellipticity, 𝜆1 and 𝜆2 are the values of the negative curvature of the electron density. The dimensionless ratios of the potential to the kinetic energy density, ǀV(r)ǀ/G(r), were calculated in order to determine the nature of bond formation and cleavage at the corresponding TSs in the RDSs.

Table 3. Topological properties of the BCPs in the RDS steps in the presence of methanthiol and thiophenol. Methanethiol Stepa 2→3

ρ(r)

ǀV(r)ǀ / G(r)

L(r)

ɛ

R

Ts

R

Ts

R

Ts

R

Ts

Se1-S6

0.000

0.015

0.000

-0.098

0.000

0.897

1.010

0.159

S6-H7

0.215

0.097

0.149

-0.016

4.693

2.556

0.058

0.174

H7-O4

0.000

0.188

0.000

0.143

0.000

3.702

0.126

0.024

O4-H3

0.331

0.260

0.513

0.324

9.457

6.537

0.019

0.022

H3-O8

0.031

0.079

-0.025

-0.033

1

1.428

0.047

0.011

O8-H10

0.032

0.017

0.489

0.164

9.119

5.072

0.020

0.023

H10-O2

0.041

0.095

-0.031

-0.013

1.009

1.809

0.051

0.004

Stepb

Thiophenol

1→2 Se1-O2

0.000

0.085

0.000

-0.031

0.000

1.372

0.143

0.062

O2-O4

0.306

0.103

0.045

-0.086

2.213

1.075

0.029

0.042

O2-H3

0.348

0.043

0.516

-0.022

9.384

1.347

0.026

0.043

O4-H3

0.000

0.103

0.000

-0.086

0.000

1.076

0.033

0.026

For atom numbering see a) Figure 2b, b) Figure 3a.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 44

A low value of the electron density and positive sign of the Laplacian are evidences for bond rupture of the O4-H3, O8-H10 at the TS23. On the other hand, the change of the Laplacian sign from the positive values in the reactant to negative values in the TS23 indicates a decrease in the charge density confirming the S6-H7, O4-H3 and O8-H10 bonds cleavage. The electron density increment of the Se1-S6, H7-O4, H3-O8 and H10-O2 bonds shows the corresponding bond formation. Moreover, an increase in the sigma bond character of the Se1-S6 and H10O2, from the reactant to TS23, has been specified by a decrease in the ellipticity values. A small value of 0.011 au. for the ellipticity character of the H3-O8 bond shows an advanced sigma character within a little distortion. Theoretical ratios of the ǀV(r)ǀ/G(r) at the BCPs are different, which confirms that they are different in chemical nature. On the basis of the electron densities and Laplacian values of the critical bonds in Table 3, it is confirmed that the O2-O4 and O2-H3 bonds are broken, while Se1-O2 and O4-H3 bonds are formed at the TS12. These changes are according to the electron density depletion/increment and Laplacian increment/depletion at the TS12, respectively. 26 ACS Paragon Plus Environment

Page 27 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

A depletion in the ellipticity values of the Se1-O2 and O2–H3 bonds from the reactant to the TS12 indicates a less distortion of the electron density and strong sigma bond. In the reactant, O2-O4 and O2-H3 bonds are the strong covalent bonds, while this character decreases at the TS confirming the bond rupture. However, Se1-O2 and O2-H3 bonds represent a middle covalent character at the TS12 which is an evidence for the corresponding bond formation. Furthermore, the analysis of electron location function and localized orbital locator of the Se1-S6 and S6-H7 bonds at the TS23, in the presence of the methanethiol, and Se1-O2 and O2-H3 bonds at TS12, in the presence of thiophenol, were shown in Figures 6 and 7, respectively. According to Figure 6, the lowest values of the ELF and LOL in the reactant regions between the Se and S atoms confirm the electrostatic interactions at the BCP, whereas higher values at the TS display the bond formation. Similarly, in Figure 7, the lowest values of the ELF and LOL in the reactant regions, between the Se and O atoms, indicate the electrostatic interactions at the BCP, while the corresponding high values at the TS regions

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 44

confirm the bond formation. These results are in agreement with the previous analyses.

Figure 6. Electron location function and localized orbital locater of the Se···S-H bonds in presence of methanethiol at the TS23.

Figure 7. Electron location function and localized orbital locater of the Se···O-H bonds in presence of thiophenol at the TS12. 28 ACS Paragon Plus Environment

Page 29 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

- NBO Analysis Donor-acceptor interaction energies were calculated at the TSs in the presence of methanethiol and thiophenol. These calculations show the effects of different thiols on the stabilization energies of the corresponding TSs, using the second order perturbation energies, E(2). Moreover, natural charges were calculated to describe the charge transfer process at the TSs. Main frontier orbital interactions of the TS23 and TS12 in the presence of methanthiol and thiophenol are reported in Table 4, respectively. In Table 4, the largest stabilization energies are related to the lpO8 O4

and lpSe1

σ*O2-H3 in the reactant in the presence of methanthiol and

thiophenol, respectively. The stabilization energies of the σSe1-O2 and lpO4

σ*H3-

σ*O8-H10

σ*Se1-O2 are the maximum values at the corresponding TSs, which

are according to the bond formation/cleavage at the TSs. The calculated sum of the main stabilization energies for the TSs are 51.67 kcal.mol −1 and 18.04 kcal.mol−1 in the presence of methanthiol and thiophenol, respectively, while the corresponding values for the reactants are, 36.67 kcal.mol−1 and 15.44 kcal.mol−1. The comparison of these stabilization energies shows that a higher 29 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 44

nucleophilicity value of methanthiol affects the overall energy of the donor−acceptor interactions. Table 4. Significant natural bond orbital interactions at the TSs of the proposed mechanism and stabilization energies E(2) (kcal.mol-1) in presence of methanethiol at thiophenol. Donor

Acceptor

E(2), (kcal.mol-1)

Methanthiol

Reactant

TS23

lpSe1

σ*S7-H8

1.04

lpS6

σ*O2-Se1

8.81

lpO8

σ*H3-O4

16.82

σO4-H3

σ*H10-O8

3.04

lpO8

σ*O4-H3

3.55

σN-Se1

σ*Se1-S6

15.81

σSe1-O2

σ*O8-H10

16.60

lpS6

σ*H7-O4

16.59

lpS6

σ*Se1-N

6.08

Thiophenol

Reactant

TS12

lpSe1

σ*O2-H3

12.1

lpO2

σ*O4-H5

1.94

lpO4

σ*O2-H3

1.40

lpO2

σ*Se1-O2

6.18

lpO4

σ*Se1-O2

6.88

σO2-H3

lpO4

1.68

σO16-H17

σ*Se1-O14

3.3

30 ACS Paragon Plus Environment

Page 31 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In order to have an insight into the charge transfer during the reaction, the change in the atomic charges of the atoms involved at the center of the TSs was evaluated and reported in Table S1. According to Table S1 and on the basis of the charge evolution or depletion on the Se1, S6, O4, O2 and H3 atoms at the TS23 and TS12, the corresponding chemical bond formation and cleavage are confirmed. Additionally, the global electron density transfer (GEDT) values of 0.03e and -0.06e in the presence of methanethiol and thiophenol were obtained, respectively. These values are according to the electronic charge increment at the TS23 in the presence of methanethiol and electronic charge depletion at the TS12 in the presence of thiophenol. 4- CONCLUSION In this research, a kinetic thermodynamic study on the designed glutathione peroxidase nano mimic 1 was done by the DFT-SAPE model. The mechanism of action of the proposed GPx-mimic 1 was analyzed by a three-step catalytic cycle in which the effects of the methanethiol and thiophenol were considered as the nucleophiles. A higher nucleophilicity of the thiophenol than methanethiol reduces the energy barrier in the catalytic reaction. The 31 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 44

significance of the dispersion interactions in the thermodynamic and kinetic studies of the proposed mechanism was considered via the mPW1PW91 functional, confirming a favorable nature of the catalytic reaction, thermodynamically. To describe the kinetic aspects of the reaction, comprehensively, energetic span model (δE) was used to evaluate the turnover frequency through which TDTS and TDI parameters were analyzed and discussed. On the basis of the location of the TDTS and TDI on the PED, the catalytic reaction of mimic 1 in the presence of thiophenol is faster than methanethiol. These analyzes show that the phenyl ring of the thiophenol is able to enrich the antioxidant activity of the new GPx-mimic 1. Moreover, QTAIM and NBO analyses were applied for investigation of the possible effects of the different thiols on the topological parameters and donor-acceptor interactions, respectively, at the TSs in the rate determining steps. Finally, the comparison of the theoretical results on the GPx mimic 1 and the experimental thermodynamic results on the glutathione peroxidase shows that mimic 1 is good enough to simulate the enzymatic behavior. Also, the

32 ACS Paragon Plus Environment

Page 33 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

comparison of the various GPx mimics which investigated theoretically, indicates that the proposed GPx mimic 1 is an efficient catalyst. Supporting Information Natural atomic charges calculated by the NBO method for the reactant and TSs at the RDS in the presence of methanethiol and thiophenol at the M06-2X/6 31+G(d,p) level of the theory; Cartesian coordinates of the studied compounds in the catalytic cycle of the GPX mimic 1 in water.

ACKNOWLEDGEMENT Research Council of Ferdowsi University of Mashhad, international Campus is acknowledged for financial supports (Grant No. 3/38560). REFERENCES: 1- Bodnar, M.; Konieczka, P.; Namiesnik, J. Polycyclic Aromatic Hydrocarbons and Digestive Tract Cancers: A Perspective. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2012, 30, 225-252. 2- Nogueira, C.W.; Rocha, J. B. Toxicology and Pharmacology of Selenium: Emphasis on Synthetic Organoselenium Compounds. Arch Toxicol. 2011, 85, 1313–1359. 33 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 44

3- Yao, Y.; Pei, F.; Kang, P. Selenium, Iodine, and the Relation with KashinBeck Disease. J. Nutr. 2011, 27, 1095−1100. 4- Steinbrenner, H.; Sies, H. Protection against Reactive Oxygen Species by Selenoproteins. Biochimica et Biophysica Acta. 2009, 1790, 1478–1485. 5- Lu, X.; Mestres, G.; Singh, V. P.; Effati, P.; Poon, J. F.; Engman. L; Ott, M. K. Selenium- and Tellurium-Based Antioxidants for Modulating Inflammation and Effects on Osteoblastic Activity. Antioxidants. 2017, 6, 13- 25. 6- Yu, F.; Li, P.; Li, G.; Zhao, G.; Zhao, G.; Chu, T.; Han. K. A Near-IR Reversible Fluorescent Probe Modulated by Selenium for Monitoring Peroxynitrite and Imaging in Living Cells. J. Am. Chem. Soc. 2011, 133, 11030– 11033. 7- Orian, L.; Mauri, P.; Roveri, A.; Toppo, S.; Benazzi,. L.; Bosello-Travain, V.; Palma, A. D.; Maiorino, M.; Miotto, G.; Zaccarin, M.; Polimeno, A.; Flohé, L.; Ursini, F. Selenocysteine Oxidation in Glutathione Peroxidase Catalysis: An MS-Supported Quantum Mechanics Study. Free Radical Biol. Med. 2015, 87, 1−14. 8- Toppo, S.; Flohé, L.; Ursini, F.; Vanin, S.; Maiorino, M. Catalytic Mechanisms and Specificities of Glutathione Peroxidases: Variations of a Basic Scheme. Biochim. Biophys. Acta. 2009, 1790, 1486–1500. 9- Lubos, E.; Loscalzo, J.; Handy, D. E. Glutathione Peroxidase-1 in Health and Disease: from Molecular Mechanisms to Therapeutic Opportunities. Antioxid Redox Signal. 2011, 15, 1957-1997. 34 ACS Paragon Plus Environment

Page 35 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

10- Tideia, C.; Piroddib, M.; Gallib, F.; Santi, C. Oxidation of Thiols Promoted by Phsezncl. Tetrahedron Lett. 2012, 53, 232–234. 11- Antony, S.; Bayse, C. A. Density Functional Theory Study of the Attack of Ebselen on a Zinc-Finger Model. Inorg. Chem. 2013, 52, 13803–13805. 12- Heverly-Coulson, G. S.; Boyd, R. J. Mechanism of the Reduction of an Oxidized Glutathione Peroxidase Mimic with Thiols. J. Chem. Theory Comput. 2012, 8, 5052−5057. 13- Arai, K.; Kumakura, F.; Takahira, M.; Sekiyama, N.; Kuroda, N.; Suzuki, T.; Iwaoka, M. Effects of Ring Size and Polar Functional Groups on the Glutathione Peroxidase-Like Antioxidant Activity of Water-Soluble Cyclic Selenides. J. Org. Chem. 2015, 80, 5633–5642. 14- Bhabak, K. P.; Mugesh, G. A Simple and Efficient Strategy to Enhance the Antioxidant Activities of Amino-Substituted Glutathione Peroxidase Mimics. Chem. Eur. J. 2008, 14, 8640–8651. 15- Schewe, T. Molecular Actions of Ebselen—an Antiinflammatory Antioxidant. Gen. Pharmacol. 1995, 26, 1153–169. 16- Antony, S.; Bayse, C. A. Modeling the Mechanism of the Glutathione Peroxidase Mimic Ebselen. Inorg. Chem. 2011, 50, 12075–12084. 17- Wang, B.; Cao, Z. Mechanism of Acid-Catalyzed Hydrolysis of Formamide from Cluster-Continuum Model Calculations: Concerted versus Stepwise Pathway. J. Phys. Chem. A. 2010, 114, 12918–12927.

35 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 44

18- Bayse, C. A. Transition States for Cysteine Redox Processes Modeled by DFT and Solvent-Assisted Proton Exchange. Org. Biomol. Chem. 2011, 9, 4748–4751. 19- Bayse, C. A.; Antony, S. Modeling the Oxidation of Ebselen and Other Organoselenium Compounds Using Explicit Solvent Networks. J. Phys. Chem. A. 2009, 113, 5780–5785. 20- Bayse, C. A.; Ortwine, K. N. Modeling the Glutathione Peroxidase-Like Activity of a Cyclic Seleninate by DFT and Solvent-Assisted Proton Exchange. Eur. J. Inorg. Chem. 2013, 21, 3680–3688. 21- Mercier, E. A.; Smith, C. D.; Parvez, M.; Back, T. G. Cyclic Seleninate Esters as Catalysts for the Oxidation of Sulfides to Sulfoxides, Epoxidation of Alkenes, and Conversion of Enamines to α-Hydroxyketones. J. Org. Chem. 2012, 77, 3508–3517. 22- Press, D. J.; McNeil, N. M. R.; Hambrook, M.; Back, T. G. Effects of Methoxy Substituents on the Glutathione Peroxidase-like Activity of Cyclic Seleninate Esters. J. Org. Chem. 2014, 79, 9394–9401. 23- Lando, P. W.; Orian, L. Peroxidase Activity of Organic Selenides: Mechanistic Insights from Quantum Chemistry. Curr. Org. Chem. 2016, 20, 189-197. 24- Kheirabadi, R.; Izadyar, M. Computational Modeling of the Catalytic Cycle of Glutathione Peroxidase Nanomimic. J. Phys. Chem. A. 2016, 120, 10108−10115. 36 ACS Paragon Plus Environment

Page 37 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

25- Bhabak, K. P.; Bhowmik, D.; Mugesh, G. Synthetic Glutathione Peroxidase Mimics: Effect of Nucleophilicity of the Aryl Thiol Cofactor on the Antioxidant Activity. Indian J. Chem. 2013, 52, 1019-1025. 26- Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A., et al. Gaussian 09, Revision A.1; Gaussian, Inc.: Wallingford, CT, 2009. 27- Kona, J.; Fabian. W. M. F. Hybrid QM/MM Calculations on the First Redox Step of the Catalytic Cycle of Bovine Glutathione Peroxidase GPx1. J. Chem. Theory Comput. 2011, 7, 2610–2616. 28- Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 other Functionals. Theor Chem Account. 2008, 120, 215–241. 29- Becke, A. D. Becke's Three Parameter Hybrid Method Using the LYP Correlation Functional. J. Chem. Phys. 1993, 98, 5648−5652. 30- Peng, C.; Schlegel, H. B.; Frisch, M. J. Using Redundant Internal Coordinates to Optimize Equilibrium Geometries and Transition States. J. Comput. Chem. 1996, 17, 49-56. 31- Niu, S.; Hall, M. B. Theoretical Studies on Reactions of Transition-Metal Complexes. Chem. Rev. 2000, 100, 353-406.

37 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 44

32- Heverly-Coulson, G. S.; Boyd, R. J. Theoretical Investigations on the Reaction of Monosubstituted Tertiary-Benzylamine Selenols with Hydrogen Peroxide.J. Phys. Chem. A. 2010, 114, 1996−2000. 33- Scalmani, G.; Frisch, M. J. Continuous Surface Charge Polarizable Continuum Models of Solvation. I. General Formalism. J. Chem. Phys. 2010, 132, 114110−114118. 34- Zhang, J.; Zhang, H.; Wu, T.; Wang, Q.; Spoel, D. V. D. Comparison of Implicit and Explicit Solvent Models for the Calculation of Solvation Free Energy in Organic Solvents. J. Chem. Theory Comput. 2017, 13, 1034−1043. 35- Adamo, C.; Barone, V. Exchange Functionals with Improved Long-Range Behavior and Adiabatic Connection Methods without Adjustable Parameters: The mpw and mpw1pw Models. J. Chem. Phys. 1998, 108, 664-670. 36- Perdew, J. P.; Wang, Y. Accurate and Simple Analytic Representation of the Electron-Gas Correlation Energy. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 45, 13244-13252. 37- Bayse, C. A.; Antony, S. Molecular Modeling of Bioactive Selenium Compounds. Main Group Chem. 2007, 6, 185 – 200. 38- Astruc, D. Organometallic Chemistry and Catalysis; Springer press: Berlin, Heidelberg, 2007. 39- Flegeau, E.; Bruneau, C.; Dixneuf, P. H.; Jutand, A. Autocatalysis for C–H Bond Activation by Ruthenium (II) Complexes in Catalytic Arylation of Functional Arenes. J. Am. Chem. Soc. 2011, 133,10161–10170. 38 ACS Paragon Plus Environment

Page 39 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

40- The Quantum Theory of Atoms in Molecules; Matta, C. F.; Boyd, R. J., Eds.; Wiley-VCH, 2007. 41- Spackman, M. A.; Byrom, P. G. A Novel Definition of a Molecule in a Crystal. Chem. Phys. Lett. 1997, 267, 215−220. 42- Scherer, W.; Sirsch, P.; Shorokhov, D.; Tafipolsky, M.; McGrady, G. S.; Gullo, E. Valence Charge Concentrations, Electron Delocalization and Β‐ Agostic Bonding in D0 Metal Alkyl Complexes. Chem. - Eur. J. 2003, 9, 6057−6070. 43- Shurki, A.; Hiberty, P.C.; Shaik, S. Charge-Shift Bonding in Group IVB Halides: A Valence Bond Study of MH3− Cl (M= C, Si, Ge, Sn, Pb) Molecules. J. Am. Chem. Soc. 1999, 121, 822−834. 44- Becke, A. D.; Edgecombe, K. E. A Simple Measure of Electron Localization in Atomic and Molecular Systems. J. Chem. Phys. 1990, 92, 5397−5403. 45- Savin, A.; Nesper, R.; Wengert, S.; Fassler, T. F. ELF: The Electron Localization Function. Chem., Int. Ed. Engl. 1997, 36, 1808−1832. 46- Burdett, J. K.; McCormick, T. A. Electron Localization in Molecules and Solids: The Meaning of ELF. J. Phys. Chem. A. 1998, 102, 6366−6372. 47- Poater, J.; Duran, M.; Solà, M.; Silvi, B. Theoretical Evaluation of Electron Delocalization in Aromatic Molecules by Means of Atoms in Molecules (AIM) and Electron Localization Function (ELF) Topological Approaches. Chem. Rev., 2005, 105, 3911–3947.

39 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 44

48- Tsirelson, V.; Stash, A. Winxpro: A Program for Calculating Crystal and Molecular Properties Using Multipole Parameters of the Electron Density. Chem. Phys. Lett. 2002, 351, 142−148. 49- Jacobsen, H. Localized-Orbital Locator (LOL) Profiles of Chemical Bonding. Can. J. Chem. 2008, 86, 695−702. 50- Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580−592. 51- Pyykko, P.; Atsumi, M. Molecular Single-Bond Covalent Radii for Elements 1-118. Chem. Eur. J. 2008, 15, 186−197. 52- Gangadharana, R. P.; Krishnanb, S. S. Natural Bond Orbital (NBO) Population Analysis of 1-Azanapthalene-8-ol. Acta Physica Polonica A. 2014, 125, 18-22. 53- Parnham, M. J.; Sies, H. The Early Research and Development of Ebselen. - NCBI. Biochem. Pharmacol. 2013, 86, 1248−1253. 54- Pearson, J. K.; Boyd, R. J. Effect of Substituents on the GPx-Like Activity of Ebselen: Steric Versus Electronic. J. Phys. Chem. A. 2008, 112, 1013−1017. 55- Kheirabadi, R.; Izadyar, M.; Housiandokht, M. R. Computational Kinetic Modeling of the Selenol Catalytic Activity as the Glutathione Peroxidase Nanomimic. J. Theor. Biol. 2016, 409, 108-114. 56- Flohe. L; The Selenoprotein Glutathione Peroxidase; Dolphin, D., Avramovic, O., Poulson. R., Eds.; John Wiley & Sons: New York, 1989, pp. 644–731 40 ACS Paragon Plus Environment

Page 41 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

57- Epp, O.; Ladenstein, R.; Wendel, A. The Refined Structure of the Selenoenzyme Glutathione Peroxidase at 0.2‐nm Resolution. FEBS J. 1983, 133, 51-69. 58- Flohe, L.; Loschen, G.; Gunzler, W.A.; Eichele, E. Glutathione Peroxidase: A Selenoenzyme. Physiol. Chem. 1973, 353, 987-999. 59- Fischer, H.; Dereu, N. Mechanism of the Catalytic Reduction of Hydroperoxides by Ebselen: A Selenium-77 NMR Study. Bull. Soc. Chim. Belg. 1987, 96, 757−68 60- Roy, G.; Nethaji, M.; Mugesh, G. Biomimetic Studies on Anti-Thyroid Drugs and Thyroid Hormone Synthesis. J. Am. Chem. Soc. 2004,126, 27122713. 61- Carmagnol. F.; Sinet. P.M.; Jerome, H. Selenium-Dependent and NonSelenium-Dependent Glutathione Peroxidases in Human Tissue Extracts. Biochim. Biophys. Acta. 1983, 759, 49-57. 62- Coulson, G. S.; Boyd, R. J. Systematic Study of the Performance of Density Functional Theory Methods for Prediction of Energies and Geometries of Organoselenium Compounds. J. Phys. Chem. A. 2011, 115, 4827–4831. 63- Wang, R.; Chen, L.; Liu, P.; Zhang, Q.; Wang, Y. Sensitive Near‐Infrared Fluorescent Probes for Thiols Based on Se-N Bond Cleavage: Imaging in Living Cells and Tissues. Chem. Eur. J. 2012, 18, 11343−11349.

41 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 44

64- Sarma, B. K.; Mugesh G. Antioxidant Activity of the Anti-Inflammatory Compound Ebselen: A Reversible Cyclization Pathway via Selenenic and Seleninic Acid Intermediates. Chem Eur J. 2008, 14, 10603-10614. 65- Zhang, C.; Lv, X.; Lu. G.; Wang, Z. X. Metal-Free Homolytic Hydrogen Activation: A Quest Through Density Functional Theory Computations. New J. Chem., 2016,40, 8141-8148. 66- Thapa, P. Zhang, R. Y. Native Chemical Ligation: A Boon to Peptide Chemistry. Molecules. 2014, 19, 14461-14483. 67- Solis, B.; Hammes-Schiffer, S. Computational Study of Anomalous Reduction Potentials for Hydrogen Evolution Catalyzed by Cobalt Dithiolene Complexes. J. Am. Chem. Soc. 2012, 134, 15253–15256. 68- Hawkes, K. J.; Cavell, K. J.; Yates, B. F. Rhodium-Catalyzed C−C Coupling Reactions: Mechanistic Considerations. Organometallics. 2008, 27, 4758–4771. 69- Guiducci, A. E.; Boyd, C. L.; Clot, E. Reactions of CyclopentadienylAmidinate Titanium Imido Compounds with CO2: Cycloaddition-Extrusion vs. Cycloaddition-Insertion. Dalton Trans. 2009, 30, 5960–5979. 70- Vahedpour, M.; Baghary, R.; Khalili. F. Prediction of Mechanism and Thermochemical Properties of O𝟑+H𝟐S Atmospheric Reaction. J. Chem. 2013, 2013, 1-12. 71- Priya, A. M.; Senthilkumar. L. Degradation of Methyl Salicylate through Cl Initiated Atmospheric Oxidation – a Theoretical Study. RSC Adv. 2014, 4, 23464-23475. 42 ACS Paragon Plus Environment

Page 43 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

"TOC Graphic"

43 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

84x65mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 44 of 44