DFT/B3LYP study of tocopherols and chromans antioxidant action energetics

Gas-phase reaction enthalpies related to the individual steps of three phenolic antioxidants action mechanisms – hydrogen atom transfer (HAT), single-electron transfer–proton transfer (SET-PT) and sequential proton loss electron transfer (SPLET) for four tocopherols and seven chromans – were calculated using DFT/B3LYP method. For α-tocopherol, one of the chromans and phenol, reaction enthalpies in water were computed. In comparison to gas phase, water causes severe changes in the energetics of studied compounds antioxidant action. From the thermodynamic point of view, entering SPLET mechanism represents the most probable process in water.     

DFT and QTAIM based investigation on the structure and antioxidant behavior of lichen substances Atranorin,Evernic acid and Diffractaic acid

In this study, the structural and antioxidant behavior of the three lichen-derived natural compounds such as atranorin (AT), evernic acid (EV) and diffractaic acid (DF) has been investigated in the gas and water phase using both B3LYP and M06-2X functional level of density functional theory (DFT) with two different basis sets 6-31+G (d, p) and 6-311++G (d, p). The intramolecular H–bonds (IHB) strength, aromaticity and noncovalent interactions (NCI) have been computed with the help of the quantum theory of atoms in molecules (QTAIM). This calculation gives major structural characteristics that indirectly influence the antioxidant behavior of the investigated compounds. The spin density (SD) delocalization of the unpaired electron is found to be the main stabilizing factor of neutral and cationic radical species. The main mechanisms, recommended in the literature, for the antioxidant action of polyphenols as radical scavengers such as hydrogen atom transfer (HAT), single electron transfer followed by proton transfer (SET-PT), and sequential proton loss electron transfer (SPLET), were examined. The result shows that the HAT and SPLET mechanism are the most conceivable one for the antioxidant action of this class of compounds in gas and water phase respectively. Preference of SPLET over HAT in water phase is due to the significantly lower value of proton affinity (PA) compared to the bond dissociation enthalpy (BDE) value. This study reveals that O₂-H₃, O₉-H₂₆ and O₄-H₄₅ respectively are the most favored site of AT, EV and DF for homolytic as well as heterolytic OH bond breaking.     

DFT/B3LYP study of the substituent effect on the reaction enthalpies of the individual steps of single electron transfer-proton transfer and sequential proton loss electron transfer mechanisms of phenols antioxidant action

The reaction enthalpies related to the individual steps of two phenolic antioxidants action mechanisms, single electron transfer-proton transfer (SET-PT) and sequential proton loss electron transfer (SPLET), for 30 meta and para-substituted phenols (ArOH) were calculated using DFT/B3LYP method. These mechanisms represent the alternative ways to the extensively studied hydrogen atom transfer (HAT) mechanism. Except the comparison of calculated reaction enthalpies with available experimental and/or theoretical values, obtained enthalpies were correlated with Hammett constants. We have found that electron-donating substituents induce the rise in the enthalpy of proton dissociation (PDE) from ArOH+* radical cation (second step in SET-PT) and in the proton affinities of phenoxide ions ArO- (reaction enthalpy of the first step in SPLET). Electron-withdrawing groups cause the increase in the reaction enthalpies of the processes where electron is abstracted, i.e., in the ionization potentials of ArOH (first step in SET-PT) and in the enthalpy of electron transfer from ArO- (second step in SPLET). Found results indicate that all dependences of reaction enthalpies on Hammett constants of the substituents are linear. The calculations of liquid-phase reaction enthalpies for several para-substituted phenols indicate that found trends hold also in water, although substituent effects are weaker. From the thermodynamic point of view, entering SPLET mechanism represents the most probable process in water.     

Novel 1,3,4-thiadiazole conjugates derived from protocatechuic acid: Synthesis,antioxidant activity,and computational and electrochemical studies

A series of 15 novel 1,3,4-thiadiazole amide derivatives containing a protocatechuic acid moiety were synthesized and structurally characterized. In addition, the corresponding imino (4) and amino (5) analogues of a phenyl-substituted 1,3,4-thiadiazole amide derivative 3a were prepared to compare the effects of the structural changes on the radical-scavenging activity. The obtained compounds were examined for their antioxidative potential by 2,2-diphenyl-1-picrylhydrazyl (DPPH) and 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) assays. In addition, selected compounds were studied by density functional theory (DFT) and cyclic voltammetry experiments. The tested compounds showed high potential to scavenging DPPH radical and ABTS radical cation compared with the referent antioxidants ascorbic acid and nordihydroguaiaretic acid (NDGA). On the basis of the calculated thermodynamic parameters, it can be concluded that the sequential proton loss electron transfer (SPLET) mechanism represents the most probable reaction path in a polar solvent for DPPH radical–scavenging activity. On the other hand, the single electron transfer followed by proton transfer (SET-PT) can be a likely mechanistic pathway in the case of an ABTS radical cation.     

Scavenging of dpph* radicals by vitamin E is accelerated by its partial ionization: the role of sequential proton loss electron transfer

[reaction: see text] Rate constants for reaction of alpha-tocopherol, 2,2,5,7,8-pentamethyl-6-hydroxychroman, and 2,6-di-tert-butyl-4-methylphenol with 2,2-diphenyl-1-picrylhydrazyl radical were measured in solvents of different polarity and H-bond basicity. In ionization supporting solvents besides hydrogen atom transfer (HAT), the kinetics of the process is partially governed by sequential proton loss electron transfer (SPLET). Addition of acetic acid reduces the rate by eliminating SPLET to leave only HAT, while addition of water increases the rate by enhancing phenol deprotonation.     

On the role of ethylene bridge elongation in the antioxidant activity of polyhydroxylated stilbenes: A theoretical approach

The aim of this work is to investigate the elongation effect of the conjugated links of the 7–8 double bond of trans-resveratrol and its analogs on the antioxidant activity in vacuo and water using a quantum chemistry calculation by the Density Functional Theory (DFT) method. H atom transfer (HAT), single-electron transfer–proton transfer (SET–PT) and sequential proton loss electron transfer (SPLET) mechanisms were investigated. The highest occupied molecular orbital (HOMO), the lowest unoccupied molecular orbital (LUMO), and the spin density were calculated. The results reveal that the elongation of the conjugated links plays an important role in promoting the antioxidant properties of molecules because of its lowering effect on BDE, spin density, AIP, and PA values. The higher antioxidant activity of 3,4 dihydroxystilbene (A₄) and trans,trans-3,4-dihydroxybistyryl (B₄) may be from the abstraction of the hydrogen atoms of the ortho-position hydroxyls. This abstraction can occur continuously to form a semiquinone structure, or even a quinone structure. On the other hand, the compounds bearing the 4,4′-DHS skeleton exhibit strong antioxidant activity due to their para-quinone structure. The results indicate correspondences between the theoretical and the experimental results. Moreover, our calculations suggest that the HAT mechanism is the most important and dominant mechanism in vacuo, the SPLET mechanism is the most thermodynamically favourable pathway in water, while the SET–PT mechanism is not preferred in all the environments studied.     

Abnormal solvent effects on hydrogen atom abstraction. 2. Resolution of the curcumin antioxidant controversy. The role of sequential proton loss electron transfer

The rates of reaction of 1,1-diphenyl-2-picrylhydrazyl (dpph*) radicals with curcumin (CU, 1,7-bis(4-hydroxy-3-methoxyphenyl)-1,6-heptadiene-3,5-dione), dehydrozingerone (DHZ, "half-curcumin"), and isoeugenol (IE) have been measured in methanol and ethanol and in two non-hydroxylic solvents, dioxane and ethyl acetate, which have about the same hydrogen-bond-accepting abilities as the alcohols. The reactions of all three substrates are orders of magnitude faster in the alcohols, but these high rates can be suppressed to values essentially equal to those in the two non-hydroxylic solvents by the addition of acetic acid. The fast reactions in alcohols are attributed to the reaction of dpph* with the CU, DHZ, and IE anions (see J. Org. Chem. 2003, 68, 3433), a process which we herein name sequential proton loss electron transfer (SPLET). The most acidic group in CU is the central keto-enol moiety. Following CU's ionization to a monoanion, ET from the [-(O)CCHC(O)-](-) moiety to dpph* yields the neutral [-(O)CCHC(O)-]* radical moiety which will be strongly electron withdrawing. Consequently, a phenolic proton is quickly lost into the alcohol solvent. The phenoxide anion so formed undergoes charge migration to produce a neutral phenoxyl radical and the keto-enol anion, i.e., the same product as would be formed by a hydrogen atom transfer (HAT) from the phenolic group of the CU monoanion. The SPLET process cannot occur in a nonionizing solvent. The controversy as to whether the central keto-enol moiety or the peripheral phenolic hydroxyl groups of CU are involved in its radical trapping (antioxidant) activity is therefore resolved. In ionizing solvents, electron-deficient radicals will react with CU by a rapid SPLET process but in nonionizing solvents, or in the presence of acid, they will react by a slower HAT process involving one of the phenolic hydroxyl groups.     

Proton-coupled electron transfer versus hydrogen atom transfer in benzyl/toluene,methoxyl/methanol,and phenoxyl/phenol self-exchange reactions

Degenerate hydrogen atom exchange reactions have been studied using calculations, based on density functional theory (DFT), for (i) benzyl radical plus toluene, (ii) phenoxyl radical plus phenol, and (iii) methoxyl radical plus methanol. The first and third reactions occur via hydrogen atom transfer (HAT) mechanisms. The transition structure (TS) for benzyl/toluene hydrogen exchange has C(2)(h)() symmetry and corresponds to the approach of the 2p-pi orbital on the benzylic carbon of the radical to a benzylic hydrogen of toluene. In this TS, and in the similar C(2) TS for methoxyl/methanol hydrogen exchange, the SOMO has significant density in atomic orbitals that lie along the C-H vectors in the former reaction and nearly along the O-H vectors in the latter. In contrast, the SOMO at the phenoxyl/phenol TS is a pi symmetry orbital within each of the C(6)H(5)O units, involving 2p atomic orbitals on the oxygen atoms that are essentially orthogonal to the O.H.O vector. The transferring hydrogen in this reaction is a proton that is part of a typical hydrogen bond, involving a sigma lone pair on the oxygen of the phenoxyl radical and the O-H bond of phenol. Because the proton is transferred between oxygen sigma orbitals, and the electron is transferred between oxygen pi orbitals, this reaction should be described as a proton-coupled electron transfer (PCET). The PCET mechanism requires the formation of a hydrogen bond, and so is not available for benzyl/toluene exchange. The preference for phenoxyl/phenol to occur by PCET while methoxyl/methanol exchange occurs by HAT is traced to the greater pi donating ability of phenyl over methyl. This results in greater electron density on the oxygens in the PCET transition structure for phenoxyl/phenol, as compared to the PCET hilltop for methoxyl/methanol, and the greater electron density on the oxygens selectively stabilizes the phenoxyl/phenol TS by providing a larger binding energy of the transferring proton.     

Free radical scavenging and COX-2 inhibition by simple colon metabolites of polyphenols: A theoretical approach

Free radical scavenging and inhibitory potency against cyclooxygenase-2 (COX-2) by two abundant colon metabolites of polyphenols, i.e., 3-hydroxyphenylacetic acid (3-HPAA) and 4-hydroxyphenylpropionic acid (4-HPPA) were theoretically studied. Different free radical scavenging mechanisms are investigated in water and pentyl ethanoate as a solvent. By considering electronic properties of scavenged free radicals, hydrogen atom transfer (HAT) and sequential proton loss electron transfer (SPLET) mechanisms are found to be thermodynamically probable and competitive processes in both media. The Gibbs free energy change for reaction of inactivation of free radicals indicates 3-HPAA and 4-HPPA as potent scavengers. Their reactivity toward free radicals was predicted to decrease as follows: hydroxyl >> alkoxyls > phenoxyl ≈ peroxyls >> superoxide. Shown free radical scavenging potency of 3-HPAA and 4-HPPA along with their high μM concentration produced by microbial colon degradation of polyphenols could enable at least in situ inactivation of free radicals. Docking analysis with structural forms of 3-HPAA and 4-HPPA indicates dianionic ligands as potent inhibitors of COX-2, an inducible enzyme involved in colon carcinogenesis. Obtained results suggest that suppressing levels of free radicals and COX-2 could be achieved by 3-HPAA and 4-HPPA indicating that these compounds may contribute to reduced risk of colon cancer development.     

The reactivity of neurotransmitters and their metabolites towards various nitrogen-centered radicals: Experimental,theoretical, and biotoxicity evaluation

In the past few years, there has been a certain interest in nitrogen-centered radicals, biologically important radicals that play a vital role in various processes and constitute many important biological molecules. In this paper, there was an attempt to bridge a gap in the literature that concerns the antiradical potency of monoamine neurotransmitters (dopamine, epinephrine, and norepinephrine) and their metabolites towards these radicals. The most probable radical quenching mechanism was determined for each radical out of three common mechanisms, namely Hydrogen Atom Transfer (HAT), Single Electron Transfer followed by the Proton Transfer (SET-PT), and Sequential Proton Loss Electron Transfer (SPLET). Marcus’ theory was then used to determine the reaction rates for the electron transfer process. SPLET was the most probable mechanism for both reactions with the aminyl and hydrazyl radicals, while HAT and SPLET were plausible mechanisms for reactions with the imidazolyl radical. Special emphasis was put on the investigation of the substituent effect on the preferred mechanism. The necessity of both thermodynamic and kinetic parameters for the comparison of the antiradical potency of compounds was discussed. The same methodology was applied for the theoretical investigation of the reactivity towards DPPH^⦁, a member of the hydrazyl radicals. An ecotoxicity analysis was performed to assess the impact the investigated radicals have on the ecosystem. Except for histidine, every other neutral form was either toxic or highly toxic to some of the analyzed marine organisms.     

Are pterins able to modulate oxidative stress?

Pterins (also known as pteridines) are common animal colorants that constitute heterocyclic compounds and have the highest nitrogen content of any pigment analyzed from animals. It has been reported that pterins modulate oxidative stress as these molecules are able to scavenge free radicals. Previous reports suggest three possible mechanisms that are responsible for scavenging free radicals; these are electron transfer (ET) reaction, hydrogen atom transfer (HAT) and radical addition. In this paper, the facility to scavenge free radicals (antiradical power) of pterins is analyzed, using density functional theory calculations and considering two possible mechanisms: ET and HAT. For the electron transfer process, considering the electron donor facility of the free radical scavenger molecules, vertical ionization energy of pterins indicates that the antiradical power of those pterins is lower than the antiradical power of any carotenoids (except for tetrahydrobiopterin). In terms of the HAT mechanism, the bond dissociation energy involved in the removal of one hydrogen atom from pterins is higher than for carotenoids (except for sepiapterin and 7,8-dihydrobiopterin). It can be expected that the most reactive molecules are those that have the smallest dissociation energy since the dissociation of the hydrogen atom is the first step of the reaction. This could indicate that some pterins are depicted as poorer antiradicals than carotenoids in terms of the HAT mechanism. Further studies focusing on the third mechanism (radical addition) and the kinetics of the reactions are necessary in order to fully understand the antiradical power of these substances. For this reason, work continues in order to clarify these aspects.     

A DFT study on reaction of eupatilin with hydroxyl radical in solution

Antioxidants scavenge reactive oxygen species and, therefore, are vitally important in the living cells. The antioxidant properties of eupatilin have recently been reported. In this article, the reactions of eupatilin with the hydroxyl radical (OH^?) in solution are studied using density functional theory calculations and the polarizable continuum model. Three mechanisms are considered including: sequential electron proton transfer (SEPT), sequential proton loss electron transfer (SPLET), and hydrogen abstraction (HA). Three solvents with different polarities, that is, benzene, methanol, and water, are used to investigate the effect of the environment on the mechanisms. The relative Gibbs free energies and enthalpies corresponding to different mechanisms are calculated. Our results show that SEPT is thermodynamically favored in aqueous solution. Once the eupatilin anion is produced, the second step in SPLET mechanism is thermodynamically favored in methanol and water. The HA mechanism is thermodynamically favored in gas, benzene, methanol, and water. This mechanism is more energetically favorable to occur in a more polar solvent. The natural bond orbital charges and spin densities as well as the singly occupied molecular orbital are then analyzed. It is concluded that the HA process is governed by proton coupled electron transfer mechanism. The attack of the radical takes place preferentially at position 7 of eupatilin. © 2012 Wiley Periodicals, Inc.     

Proton-coupled electron transfer in solution, proteins, and electrochemistry

Recent advances in the theoretical treatment of proton-coupled electron transfer (PCET) reactions are reviewed. These reactions play an important role in a wide range of biological processes, as well as in fuel cells, solar cells, chemical sensors, and electrochemical devices. A unified theoretical framework has been developed to describe both sequential and concerted PCET, as well as hydrogen atom transfer (HAT). A quantitative diagnostic has been proposed to differentiate between HAT and PCET in terms of the degree of electronic nonadiabaticity, where HAT corresponds to electronically adiabatic proton transfer and PCET corresponds to electronically nonadiabatic proton transfer. In both cases, the overall reaction is typically vibronically nonadiabatic. A series of rate constant expressions have been derived in various limits by describing the PCET reactions in terms of nonadiabatic transitions between electron-proton vibronic states. These expressions account for the solvent response to both electron and proton transfer and the effects of the proton donor-acceptor vibrational motion. The solvent and protein environment can be represented by a dielectric continuum or described with explicit molecular dynamics. These theoretical treatments have been applied to numerous PCET reactions in solution and proteins. Expressions for heterogeneous rate constants and current densities for electrochemical PCET have also been derived and applied to model systems.     

A DFT study of proton transfers for the reaction of phenol and hydroxyl radical leading to dihydroxybenzene and H2O in the water cluster

Reactions of phenol and hydroxyl radical were studied under the aqueous environment to investigate the antioxidant characters of phenolic compounds. M06‐2X/6‐311 + G(d,p) calculations were carried out, where proton transfers via water molecules were examined carefully. Stepwise paths from phenol + OH^• + (H₂O)_n (n = 3, 7, and 12) to the phenoxyl radical (Ph O^•) via dihydroxycyclohexadienyl radicals (ipso, ortho, meta, and para OH‐adducts) were obtained. In those paths, the water dimer was computed to participate in the bond interchange along hydrogen bonds. The concerted path corresponding to the hydrogen atom transfer (HAT, apparently Ph OH + OH^• → Ph O^• + H₂O) was found. In the path, the hydroxyl radical located on the ipso carbon undergoes the charge transfer to prompt the proton (not hydrogen) transfer. While the present new mechanism is similar to the sequential proton loss electron transfer (SPLET) one, the former is of the concerted character. Tautomerization reactions of ortho or para (OH)C₆H₅=O + (H₂O)_n → C₆H₄(OH)₂(H₂O)_n were traced with n = 2, 3, 4, and 14. The n = 3 (and n = 14) model of ortho and para was calculated to be most likely by the strain‐less hydrogen‐bond circuit.     

Antioxidant activity of trans-resveratrol toward hydroxyl and hydroperoxyl radicals: a quantum chemical and computational kinetics study

In this work, we have carried out a systematic study of the antioxidant activity of trans-resveratrol toward hydroxyl ((?)OH) and hydroperoxyl ((?)OOH) radicals in aqueous simulated media using density functional quantum chemistry and computational kinetics methods. All possible mechanisms have been considered: hydrogen atom transfer (HAT), proton-coupled electron transfer (PCET), sequential electron proton transfer (SEPT), and radical adduct formation (RAF). Rate constants have been calculated using conventional transition state theory in conjunction with the Collins-Kimball theory. Branching ratios for the different paths contributing to the overall reaction, at 298 K, are reported. For the global reactivity of trans-resveratrol toward (?)OH radicals, in water at physiological pH, the main mechanism of reaction is proposed to be the sequential electron proton transfer (SEPT). However, we show that trans-resveratrol always reacts with (?)OH radicals at a rate that is diffusion-controlled, independent of the reaction pathway. This explains why trans-resveratrol is an excellent but very unselective (?)OH radical scavenger that provides antioxidant protection to the cell. Reaction between trans-resveratrol and the hydroperoxyl radical occurs only by phenolic hydrogen abstraction. The total rate coefficient is predicted to be 1.42 × 10(5) M(-1) s(-1), which is much smaller than the ones for reactions of trans-resveratrol with (?)OH radicals, but still important. Since the (?)OOH half-life time is several orders larger than the one of the (?)OH radical, it should contribute significantly to trans-resveratrol oxidation in aqueous biological media. Thus, trans-resveratrol may act as an efficient (?)OOH, and also presumably (?)OOR, radical scavenger.     

On the direct scavenging activity of melatonin towards hydroxyl and a series of peroxyl radicals

The reactions of melatonin (MLT) with hydroxyl and several peroxyl radicals have been studied using the Density Functional Theory, specifically the M05-2X functional. Five mechanisms of reaction have been considered: radical adduct formation (RAF), Hydrogen atom transfer (HAT), single electron transfer (SET), sequential electron proton transfer (SEPT) and proton coupled electron transfer (PCET). It has been found that MLT reacts with OH radicals in a diffusion-limited way, regardless of the polarity of the environment, which indicates that MLT is an excellent OH radical scavenger. The calculated values of the overall rate coefficient of MLT + ˙OH reaction in benzene and water solutions are 2.23 × 10(10) and 1.85 × 10(10) M(-1) s(-1), respectively. MLT is also predicted to be a very good ˙OOCCl(3) scavenger but rather ineffective for scavenging less reactive peroxyl radicals, such as alkenyl peroxyl radicals and ˙OOH. Therefore it is concluded that the protective effect of MLT against lipid peroxidation does not take place by directly trapping peroxyl radicals, but rather by scavenging more reactive species, such as ˙OH, which can initiate the degradation process. Branching ratios for the different channels of reaction are reported for the first time. In aqueous solutions SEPT was found to be the main mechanism for the MLT + ˙OH reaction, accounting for about 44.1% of the overall reactivity of MLT towards this radical. The good agreement between the calculated and the available experimental data, on the studied processes, supports the reliability of the results presented in this work.     

MPW1K Performs Much Better than B3LYP in DFT Calculations on Reactions that Proceed by Proton-Coupled Electron Transfer (PCET)

DFT calculations have been performed with the B3LYP and MPW1K functional on the hydrogen atom abstraction reactions of ethenoxyl with ethenol and of phenoxyl with both phenol and alpha-naphthol. Comparison with the results of G3 calculations shows that B3LYP seriously underestimates the barrier heights for the reaction of ethenoxyl with ethenol by both proton-coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms. The MPW1K functional also underestimates the barrier heights, but by much less than B3LYP. Similarly, comparison with the results of experiments on the reaction of phenoxyl radical with alpha-naphthol indicates that the barrier height for the preferred PCET mechanism is calculated more accurately by MPW1K than by B3LYP. These findings indicate that the MPW1K functional is much better suited than B3LYP for calculations on hydrogen abstraction reactions by both HAT and PCET mechanisms.     

Exploring the transfer of hydrogen atom from kaempferol-based compounds to hydroxyl radical at ground state using PCM-DFT approach

Thermodynamic and kinetic studies of the hydrogen atom transfer (HAT) from hydroxyl (OH) groups of four kaempferol-based compounds, namely kaempferol, morin, morin-5^*-sulfonate and morin-7-O-sulfate to hydroxyl radical (·OH), have been carried out using density functional theory (DFT) methods at the CAM-B3LYP/6–311++G(d,p) level equipped with polarizable continuum model (PCM) of solvation. All HAT reactions in aqueous solution are exothermic and spontaneous. For most compounds, the most preferable OH group for HAT is situated at position C3 (O3-H3) on the pyrone ring. The reaction potential of such a reactive group is found to be highest in morin-7-O-sulfate. The rate constants for the HAT reactions at different OH groups of each compound have been determined based on the transition state theory. The presence of substituents leads to the variation in either the characteristic interactions at the reactive site or the charge distribution on transition-state geometries, hence significantly affecting the kinetics of HAT. The highest rate of HAT is resulted for the OH group at position C4^* (O4^*-H4^*) on the phenyl ring (ring B) of morin-5^*-sulfonate because a hydrogen bond between ·OH and the sulfonate group favors the formation of transition state. However, for most compounds under study, the HAT reaction at O3-H3 initiated by ·OH is highly favorable both thermodynamically and kinetically.

     

A combined experimental and theoretical approach for radical-scavenging activity of edaravone and its related derivatives

The experimental and theoretical study for evaluation of scavenging activity of edaravone (S1) and related derivatives, such as antipyrine (S2), dipyrone (S3), and phenylbutazone (S4), was carried out against DPPH and ABTS radicals. Structure–activity relationship study was performed using quantum chemical calculations at the DFT/B3LYP level of theory along with the 6-31G* basis sets. S1 and S4 are more effective scavengers against DPPH and ABTS. We observed little effects of S2 and S3 at several concentrations against these two free radicals. The calculations of HOMO, ionization potential, and bond dissociation energy confirmed that a hydrogen transfer is more preferential than an electron transfer. The radical stability of these compounds is related with spin densities. In accordance with experimental and theoretical results, edaravone is more active than phenylbutazone as scavenging drug.