Beena G.
Singh
*ab,
Pavitra
Kumar
ab,
P.
Phadnis
c,
Michio
Iwaoka
d and
K. Indira
Priyadarsini
*bc
aRadiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India. E-mail: beenam@barc.gov.in
bHomi Bhabha National Institute, Anushaktinagar, Mumbai-400092, India. E-mail: kindira@barc.gov.in
cChemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
dDepartment of Chemistry, School of Science, Tokai University, Kitakaname, Hiratsuka-shi, Kanagawa 259-1292, Japan
First published on 26th July 2019
With an aim to understand the effect of structural modifications on antioxidant activity, two structurally isomeric selenium compounds, linear, bis(ethan-2-ol)selenide (SeEOH), and cyclic, DL-trans-3,4-dihydroxy-1-selenolane (DHS), were compared. The antioxidant activity was assigned through their ability to scavenge reactive oxygen species (ROS), where both compounds indicated formation of selenoxides. The ROS identified for the present study are peroxynitrite, hydroxyl (˙OH) and carbonate (CO3˙−) radicals. The rate constant for the scavenging of both peroxynitrite and the CO3˙− radical was higher for DHS as compared to that for SeEOH. The one-electron oxidation during these reactions indicated the formation of a selenium centred dimer radical cation (Se∴Se
)+, as the intermediate, which undergoes the disproportionation reaction to form the corresponding selenoxide. Quantitative estimation of selenoxides indicated that cyclic compound DHS produced a higher yield of selenoxides than the linear compound SeEOH due to higher stability of (
Se∴Se
)+. The selenoxides undergo reduction by thiols, thereby regenerating back to the parent compound. The results thereby suggest that formation of selenoxides improves the antioxidant activity of organoselenium compounds and should be considered as an important step in the design of new selenium based antioxidants.
In search of water soluble and low molecular weight antioxidants/radioprotectors, our group has studied the hydroxyl (˙OH) radical reaction with different classes of organoselenium compounds by the pulse radiolysis technique.11 The studies have shown that the initial attack of the ˙OH radical takes place at the selenium centre to form a hydroxyselenouranyl radical adduct (Se∴OH) and the fate of this transient depends on the nature of the functional group (hetero atom) attached to the selenium centre. Of the several compounds screened in our laboratory, trans-3,4-dihydroxyselenolane (DHS), a water-soluble organoselenium compound, has been found to reduce ˙OH radicals, hydroperoxides and peroxynitrite catalytically through the formation of a stable selenoxide, which is reverted back to DHS by thiols.12 Furthermore, this antioxidant activity was also verified using in vitro and in vivo models.13 In contrast, its linear analogue bis(ethan-2-ol)selenide (SeEOH) showed low GPx like activity and was toxic to cells.14 Therefore, to understand the contrasting behaviour, two structural analogues have been investigated for their reactivity with ROS and the yields of the products formed were estimated. To further complement the observations, quantum chemical calculations have been carried out. The structures of DHS, SeEOH and their selenoxides are given in Scheme 1.
![]() | ||
Scheme 1 Structure of 3,4-dihydroxy selenolane (DHS), bis(ethan-2-ol)selenide (SeEOH) and their corresponding selenoxide DHSox and SeEOHox. |
The reaction between DHS/SeEOH and peroxynitrite was studied by employing competition kinetics using dihydrorhodamine 123 (DHR123) as a reference solute.15
![]() | (1) |
![]() | (2) |
The bimolecular rate constant (k2) was estimated by plotting (F0/F) − 1 as a function of ((SeEOH or DHS)/DHR123) according to eqn (3):
![]() | (3) |
The rate constant for the reaction of DHS with the CO3˙− radical was reported to be (1.2 ± 0.2) × 109 M−1 s−1.12 In the present study, SeEOH was tested for its ability to scavenge CO3˙− and NO2˙ radicals using the pulse radiolysis technique. The kinetics between SeEOH and the CO3˙− radical was estimated by following the characteristic decay of CO3˙− at ∼600 nm (ε600nm = 1680 M−1 cm−1) in the presence of varying concentration of SeEOH. Fig. S2 (ESI†) shows the absorption–time plot of the CO3˙− radical in the absence and presence of 10 μM SeEOH. In the absence of SeEOH, the CO3˙− radical decayed by second order kinetics with a 2k/εl value of (5.7 ± 0.1) × 107 s−1 at pH 7.4. On addition of 10 μM SeEOH, the CO3˙− radical decayed very fast and the reaction was found to follow pseudo-first order kinetics (kobs). The bimolecular rate constant was estimated from the slope of the linear plot of kobs at different concentrations of SeEOH (10–100 μM), and was found to be (6.5 ± 0.3) × 108 M−1 s−1. Both DHS and SeEOH did not show any reaction with the NO2˙ radical, which could be due to the unfavourable redox potential of these compounds.
To further understand the mechanism of electron transfer in these compounds, the nature of the transients produced during these ROS reactions was investigated as discussed below.
The relative yield of the radical cations, estimated by the ABTS˙−18,19 reaction, was found to be 1.8 times higher for DHS than SeEOH. As reported earlier, the (
Se∴Se
)+ radical undergoes radical–radical disproportionation to form selenoxides, as estimated below.
To explain the differential yield of selenoxides obtained during reaction with ˙OH radicals and their ROS scavenging activity, quantum chemical calculations were performed and the results are discussed below.
To explain the nature and stability of the resultant (Se∴Se
)+ radicals of SeEOH and DHS, their geometries were optimized followed by correction for BSSE (basis-set superposition error) using the counter-poise method. Fig. 2 represents the optimized ground state structure of the (
Se∴Se
)+ radical at the B3LYP/6.31+G(d,p) level. The bond distance between two selenium atoms is less (3.035 Å) in the (
Se∴Se
)+ radical of SeEOH which enhances the inter electronic repulsion and destabilizes the bond. On the other hand in the (
Se∴Se
)+ radical of DHS, the Se–Se bond distance (3.071 Å) is slightly more leading to a decrease in interelectronic repulsion and stabilisation of the two centred three electron bonds in the (
Se∴Se
)+ radical. Additionally, (
Se∴Se
)+ of DHS is stabilized by two non-bonded Se⋯O interactions from the hydroxyl functional group, which is not seen in SeEOH. Therefore, the (
Se∴Se
)+ radical of DHS is more stable than that of SeEOH.
![]() | ||
Fig. 2 Optimized geometries of the (![]() ![]() |
Furthermore the energetics for the initial reaction of the ˙OH radical with these compounds and consequent formation of the (Se∴Se
)+ radicals was calculated. The results showed that the initial attack of the ˙OH radical at the Se-centre to form the (
Se∴OH) radical is more preferred in DHS (ΔE = −17.29 kcal mol−1) than SeEOH (ΔE = +3.49 kcal mol−1). This (
Se∴OH) generally decays by proton catalyzed dehydration to form a selenium centered radical cation (
Se˙+) which on reaction with another parent molecule forms (
Se∴Se
)+. The dimerization process was found to be more exothermic in DHS (−19.3 kcal mol−1) as compared to that in SeEOH (−14.9 kcal mol−1). These results suggest that the formation of the (
Se∴Se
)+ radical is more favourable in the case of DHS than SeEOH and support the experimental observations. The preferential formation of the (
Se∴Se
)+ radical can also be explained on the basis of HOMO energy values of different transients involved in its formation. Table 1 lists the HOMO, LUMO and HOMO–LUMO energy gaps in different transients formed for SeEOH and DHS. The HOMO energy values of DHS and the DHS˙+ radical are comparable and therefore an efficient overlap is expected.
Compounds | HOMO (eV) | LUMO (eV) | HOMO–LUMO energy gap (eV) |
---|---|---|---|
SeEOH | −5.7822 | −0.2829 | −5.4992 |
(SeEOH)˙+ | −7.8693 | −5.2380 | −2.6312 |
(SeEOH)2˙+ | −5.7414 | −3.5918 | −2.1496 |
DHS | −5.6625 | −0.1006 | −5.5618 |
(DHS)˙+ | −6.8081 | −1.3496 | −5.4584 |
(DHS)2˙+ | −6.0489 | −0.9088 | −5.1401 |
The monomer radical cation of organoselenium compounds tends to undergo irreversible deprotonation to form a carbon centred radical (α-reducing radical) on the carbon adjacent to the chalcogen centre. In our earlier work, we have shown that the deprotonation is favoured by the effective overlap between the σ* of the selenium centre and the σ-C–H bond at the α-position. The energy of such non-bonding interactions can theoretically be calculated by employing natural bond orbital (NBO) analysis.22 The non-bonding interaction energy (Enb) indicates the extent of the interactions taking place between two orbitals. The optimized structure of (Se˙+) derived from DHS/SeEOH was obtained and the Enb value of DHS and SeEOH was calculated to be 0.36 kcal mol−1 and 1.61 kcal mol−1, respectively. From this it can be inferred that the (
Se˙+) derived from DHS is less prone to undergo irreversible deprotonation as compared to SeEOH. Furthermore, the energetics for conversion of (
Se∴Se
)+ to form selenoxides indicated that this radical in DHS disproportionates in an energetically favored path (ΔE = +31.1 kcal mol−1), as compared to SeEOH (ΔE = +32.5 kcal mol−1).
On reaction with ˙OH radical, both the compounds were converted to their corresponding selenoxide, where DHS produced a significantly higher amount of selenoxide than SeEOH.
The experimental results were complemented by calculating the energetics determined by quantum chemical calculations. The HOMO energy of DHS is higher than that of SeEOH indicating the easier oxidation of DHS. The reaction of initial addition of the ˙OH radical on the selenium atom to form a (Se∴OH) adduct is endothermic in SeEOH (ΔE = +3.48 kcal mol−1) but exothermic in DHS (ΔE = −17.29 kcal mol−1). This (
Se∴OH) adduct generally decays to form (
Se∴Se
)+. The dimerization process was found to be more exothermic in DHS as compared to SeEOH. Additionally, the stability of the dimer radical cation of DHS is due to the direct interaction between selenium and oxygen of DHS and this prevents deprotonation to form a carbon centred radical. Thus a higher HOMO value coupled with higher stability of the intermediate (
Se∴Se
)+ results in a higher yield of selenoxide formation and thereby resulting in better antioxidant activity of DHS as compared to SeEOH. The selenoxide being reversible in the presence of thiols imparts higher antioxidant activity in DHS than in SeEOH against ROS.
Peroxynitrite was prepared by the ozonolysis of 0.1 M sodium azide solution containing 0.1 M NaOH.23 For all experiments stock peroxynitrite solutions were diluted using 10 mM NaOH and during the experiments, the pH was maintained at 7.4 by using 70 mM phosphate buffer. The concentration of peroxynitrite was estimated by measuring the absorbance at 302 nm (ε302nm = 1705 M−1 cm−1) on a UV-Visible spectrophotometer (Jasco V-639).
Pulse radiolysis studies were carried using a 7 MeV electron beam with 100 ns pulse width and the details of the instrument are reported. An average dose of 9 Gy as estimated by using a thiocyanate dosimeter (aerated aqueous solution of 10 mM KSCN, Gε475nm = 2.59 × 10−4 m2 J−1) was used for all experiments.24 Solutions were saturated with N2O to monitor the ˙OH radical reaction (G˙OH = 0.6 μmole J−1). The carbonate (CO3˙−) radical was generated on pulse radiolyzing N2O saturated solution containing sodium carbonate (0.1 M) at pH 7. The product formed on reaction of the ˙OH radical with the organoselenium compound was quantified by the HPLC technique. For this, N2O saturated aqueous solutions of 10 mM SeEOH were radiolysed using the 60Co γ-source with a dose rate of 40 Gy min−1. The absorbed dose was set in a way that not more than 10% of the compound undergoes radiolysis. Selenoxide estimation was performed as per the reported method. In brief, 180 μl of radiolysed samples were treated with 20 μl DTT (10 mM) and incubated for 5 minutes. 20 μl of the mixture was injected into the C-18 reverse phase column. The different compounds were eluted using an acetonitrile:water (5
:
95) mixture as the mobile phase and absorption peaks were detected at 240 nm. The yield of formaldehyde formed during radiolysis of the organoselenium compound was estimated by derivatizing the aldehyde with DNPH. The radiated solutions were mixed with 1 mM DNPH dissolved in 10 mM HCl and stirred for 5 minutes. 20 μl of the resulting solutions were injected into the C-18 reverse phase column. Acetonitrile:water (60
:
40) was used as the mobile phase and the detector was set at 345 nm.
Reaction of the organoselenium compounds with peroxynitrite was studied by competition kinetics using dihydrorhodamine (DHR)123 as standard. The fluorescence measurements were done on a Hitachi F-4500 fluorescence spectrophotometer with excitation and emission wavelengths of 510 nm and 536 nm respectively. Peroxynitrite solution (5 μM) was added to DHR123 (10 μM) containing 0.1 mM DTPA in 70 mM phosphate buffer (pH 7.5) in the absence and presence of the organoselenium compound (10–100 μM).
The structures of the transients were optimized in vacuo by regressive variation in the starting geometry at the B3LYP/6-31++G(d,p) level (the Becke nonlocal model and Lee–Yang–Parr nonlocal correlation functionals). The geometries obtained were checked by frequency calculations. The global minimum structures were then further optimized in water at the B3LYP/6-31+G(d,p) level in water using the PCM-SMD model. Geometry optimization and frequency calculations were performed by adopting the GAMESS suite of programs on a PC-based LINUX cluster platform.25 The absorption maximum wavelength was calculated in water by using the UCIS model at the B3LYP;6-31-G(d,p) level using Gaussian 09. The transient's geometries and molecular orbitals were visualized using chemissian V4.38 software.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9nj02227a |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2019 |