Modeling the reduction of hydrogen peroxide by glutathione peroxidase mimics

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

Effect of substituents on the GPx-like activity of ebselen: steric versus electronic

The direct oxidation of ebselen and several derivatives by hydrogen peroxide is investigated using the B3LYP/6-31G(d,p) method to elucidate the effects of substituents on GPx-like activity. While previous studies have attributed the differences in GPx activity of substituted ebselen compounds to the electronic nature of the substituents, the influence of functional groups is poorly understood. The effects of various solvents are incorporated by employing the CPCM method. It is shown that a substituent in the ortho position to the selenium atom sterically hinders attack of a nucleophile at selenium and thus increases the barrier to reaction. The observed increase in GPx-like activity of an ebselen derivative with an ortho substituent is explained by the fact that the steric hindrance prevents thiol exchange reactions.     

取代基对依布硒啉电子结构及光谱性质影响的理论研究

采用密度泛函理论B3LYP方法在6-311++G**水平对依布硒啉及其6个衍生物进行结构优化,并分析了取代基对分子几何构型、电荷分布、前线轨道能级分布的影响.在相同计算水平上采用含时密度泛函理论方法进行了电子光谱研究,讨论了取代基对电子光谱的影响.计算结果表明,标题化合物分子为非平面结构,苯并异硒唑酮环与苯环之间存在一...     

Synthesis and antioxidant activity of peptide-based ebselen analogues

A series of di- and tripeptide-based ebselen analogues has been synthesized. The compounds were characterized by (1)H, (13)C, and (77)Se NMR spectroscopy and mass spectral techniques. The glutathione peroxidase (GPx)-like antioxidant activity has been studied by using H(2)O(2) , tert-butyl hydroperoxide (tBuOOH), and cumene hydroperoxide (Cum-OOH) as substrates, and glutathione (GSH) as a cosubstrate. Although all the peptide-based compounds have a selenazole ring similar to that of ebselen, the GPx activity of these compounds highly depends on the nature of the peptide moiety attached to the nitrogen atom of the selenazole ring. It was observed that the introduction of a phenylalanine (Phe) amino acid residue in the N-terminal reduces the activity in all three peroxide systems. On the other hand, the introduction of aliphatic amino acid residues such as valine (Val) significantly enhances the GPx activity of the ebselen analogues. The difference in the catalytic activity of dipeptide-based ebselen derivatives can be ascribed mainly to the change in the reactivity of these compounds toward GSH and peroxide. Although the presence of the Val-Ala-CO(2) Me moiety facilitates the formation of a catalytically active selenol species, the reaction of ebselen analogues that has a Phe-Ile-CO(2) Me residue with GSH does not generate the corresponding selenol. To understand the antioxidant activity of the peptide-based ebselen analogues in the absence of GSH, these compounds were studied for their ability to inhibit peroxynitrite (PN)-mediated nitration of bovine serum albumin (BSA) and oxidation of dihydrorhodamine 123. In contrast to the GPx activity, the PN-scavenging activity of the Phe-based peptide analogues was found to be comparable to that of the Val-based compounds. However, the introduction of an additional Phe residue to the ebselen analogue that had a Val-Ala dipeptide significantly reduced the potency of the parent compound in PN-mediated nitration.     

Synthesis, characterization, and antioxidant activity of some ebselen analogues

Simple synthetic routes for several analogues of the anti-inflammatory organoselenium drug, ebselen, are described. The compounds are characterized by (1)H, (13)C, and (77)Se NMR spectroscopy and mass spectral techniques and, in some cases, by single-crystal X-ray diffraction studies. The glutathione peroxidase (GPx)-like antioxidant activity has been studied by using H(2)O(2), tBuOOH, and Cum-OOH as substrates, and thiophenol (PhSH, 4-Me-C(6)H(4)SH) and glutathione (GSH) as cosubstrates. Density functional theory (DFT) calculations have been performed on these systems to understand the effects of various substituents on the (77)Se NMR chemical shifts; these results have been compared with the experimental data. The experimental and theoretical results suggest that the presence of a phenyl substituent on the nitrogen atom is important for the antioxidant activity of ebselen. While ebselen and its analogues are poor catalysts in aromatic thiol assays, these compounds exhibit high GPx activity when GSH is used as the cosubstrate. The poor catalytic activity of ebselen analogues in the presence of aromatic thiols such as PhSH and 4-Me-C(6)H(4)SH can be ascribed to the undesired thiol exchange reaction that takes place at the selenium center due to SeO nonbonding interactions. To understand the effects of different peroxides on the catalytic activities, we have determined the initial rates at various concentrations of GSH and peroxides. These data suggest that the nature of peroxide has little effect on the catalytic efficiencies, although the initial reaction rates observed with hydrogen peroxide were found to be higher than that with tBuOOH and Cum-OOH. In contrast to the effect of peroxides, the nature of thiols appears to have a dramatic effect on the catalytic activity of ebselen and its related derivatives.     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Pyridoxine‐Derived Organoselenium Compounds with Glutathione Peroxidase‐Like and Chain‐Breaking Antioxidant Activity

One of the vitamin B6 vitamers, pyridoxine, was modified to incorporate selenium in various oxidation states in place of the methyl group in position 2. Such compounds were conveniently accessed by treatment of bis‐4,5‐(carboethoxy)‐2‐iodo‐3‐pyridinol with disodium diselenide and LiAlH₄‐reduction. After work‐up, selone 7 was isolated in good yield as an air‐stable crystalline material. Hydrogen bonding to the neighboring hydroxyl group, as revealed by the short intramolecular Se ??? H distance in the crystal structure is likely to provide extra stabilization to the compound. Computational studies showed that selone 7 is more stable than the corresponding selenol tautomer by 12.2 kcal mol^?1. Hydrogen peroxide oxidation of the selone 7 afforded diselenide 12 , and, on further oxidation, seleninic acid 13 . Treatment of the seleninic acid with thiophenol provided an isolable selenosulfide 14 . The glutathione peroxidase‐like properties of the pyridoxine‐derived compounds were assessed by using the coupled reductase method. Seleninic acid 13 was found to be twofold more active than ebselen. The chain‐breaking capacity of the pyridoxine compounds were studied in a water/chlorobenzene membrane model containing linoleic acid as an oxidizable substrate and N‐acetylcysteine as a thiol reducing agent. Diselenide 15 could match α‐tocopherol when it comes to reactivity towards peroxyl radicals and inhibition time.     

含氮ZSM-5分子筛骨架中氮取代位置的密度泛函计算

运用Gaussian 98程序包, 采用密度泛函理论B3LYP方法, 基于ZSM-5分子筛的8T模型, 分别通过6-31G, 6-31G(d)和6-311G(d,p)基组计算了ZSM-5分子筛中氮原子取代前后各O原子和各N原子的能量, 从而得到各O原子与各N原子在骨架中的稳定性及其对氮化取代反应的影响. 计算结果表明, N原子在骨架中的稳定性对氮取代反应的影响较大. ZSM-5分子筛晶体结构中与B酸位处于同一个四面体的O11位置, 为氮原子的最佳取代位置, 因此氮化后分子筛表面的B酸强度得到较大程度的减弱.     

Importance of selenium in antioxidant behavior of ebselen: A theoretical study

The electronic structure, geometrical parameters, proton affinities, lone pair occupancies, and orbital interactions have been explored to understand the antioxidant behavior of ebselen at B3LYP/6-31+G1 level. The role of substituents in the nucleophilicity of selenium in eight model molecules containing SeN bond has been analyzed in terms of relative proton affinities. The free energy changes associated with probable catalytic reaction steps have also been evaluated.     

Kinetic Study of Di‐Tert‐Butyl Peroxide: Thermal Decomposition and Product Reaction Pathways

The major reaction path for oxidation of di‐tert‐butyl peroxide (DTBP) is generally considered to occur via fission of the weak peroxide RO? OR bond at temperatures above 393 K. The initial stable intermediates in the thermal decomposition or combustion of DTBP are acetone and ethane, and the overall reaction is accompanied by an important heat release which when mixed with air (oxygen) may exceed the self‐ignition temperatures. A kinetic study on plausible DTBP reaction paths was initiated in this work, and a detailed study of the thermochemistry of new intermediates, transition state structures, and products is reported. The density functional theory (DFT; B3LYP/6‐311g(d,p)), which is practical for large compounds along with the composite ab initio G3MP2B3 and G3 calculations, (when possible), are used. Computational chemistry results from DFT and ab initio calculations are coupled with isodesmic reaction analysis which, as demonstrated in previous studies, results in good accuracy. Over 10 unimolecular decomposition pathways are identified and reported.     

Spectroscopic and quantum chemical study of the novel compound cyclopropylmethylselenol

An investigation into the properties of the novel compound cyclopropylmethylselenol has been undertaken by use of Stark-modulation microwave spectroscopy and high-level quantum chemical calculations. Ground-state spectra belonging to six isotopomers of a single conformer of the molecule were recorded and assigned. This conformer, predicted to be the lowest in energy by a series of quantum chemical calculations, possesses a synclinal arrangement of the H-C-C-Se atoms. In addition to the assignment of these ground-state spectra, transitions attributable to vibrationally excited states of the 78Se- and 80Se-containing isotopomers were identified. A tentative assignment of these excited-state spectra to specific vibrational modes has been made with the assistance of a density functional theory calculation at the B3LYP/6-311++G(3df,2pd) level of theory. Close agreement was found between experimental ground-state rotational constants and ab initio equilibrium values calculated at the MP2/aug-cc-pVTZ level of theory. Good agreement was also noted between certain r(s) principal axis coordinates of atoms in the molecule and the corresponding ab initio r(e) values. Limited evidence in favor of the formation of a weak intramolecular hydrogen bond between the H atom of the selenol group and electron density associated with the cyclopropyl ring is discussed.     

Antioxidant activity of the anti-inflammatory compound ebselen: a reversible cyclization pathway via selenenic and seleninic acid intermediates

A revised mechanism that accounts for the glutathione peroxidase (GPx)-like catalytic activity of the organoselenium compound ebselen is described. It is shown that the reaction of ebselen with H(2)O(2) yields seleninic acid as the only oxidized product. The X-ray crystal structure of the seleninic acid shows that the selenium atom is involved in a noncovalent interaction with the carbonyl oxygen atom. In the presence of excess thiol, the Se--N bond in ebselen is readily cleaved by the thiol to produce the corresponding selenenyl sulfide. The selenenyl sulfide thus produced undergoes a disproportionation in the presence of H(2)O(2) to produce the diselenide, which upon reaction with H(2)O(2), produces a mixture of selenenic and seleninic acids. The addition of thiol to the mixture containing selenenic and seleninic acids leads to the formation of the selenenyl sulfide. When the concentration of the thiol is relatively low in the reaction mixture, the selenenic acid undergoes a rapid cyclization to produce ebselen. The seleninic acid, on the other hand, reacts with the diselenide to produce ebselen as the final product. DFT calculations show that the cyclization of selenenic acids to the corresponding selenenyl amides is more favored than that of sulfenic acids to the corresponding sulfenyl amides. This indicates that the regeneration of ebselen under a variety of conditions protects the selenium moiety from irreversible inactivation, which may be responsible for the biological activities of ebselen.     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

肝醇脱氢酶催化乙醇氧化生成乙醛反应机理的理论研究

 采用 B3LYP 方法研究了肝醇脱氢酶催化烟酰胺腺嘌呤二核苷酸氧化乙醇生成乙醛的反应机理. 优化得到了反应物、过渡态、中间体和产物的几何构型, 并计算了在蛋白质或水环境下有或没有肝醇脱氢酶时的反应势垒. 结果表明, 没有催化剂时, 乙醇负离子的形成及其被氧化生成乙醛的反应势垒都很高, 常温下反应难以进行; 当肝醇脱氢酶存在时, 乙醇负离子可以与肝醇脱氢酶中的 Zn2+配位形成络合物, 从而极大地降低了这两步的反应势垒, 使得反应在常温下容易进行.     

Computational study of the aminolysis of anhydrides: effect of the catalysis to the reaction of succinic anhydride with methylamine in gas phase and nonpolar solution

Different possible pathways of the aminolysis reaction of succinic anhydride were investigated by applying high level electronic structure theory, examining the general base catalysis by amine and the general acid catalysis by acetic acid, and studying the effect of solvent. The density functional theory at the B3LYP/6-31G(d) and B3LYP/6-311++G(d,p) levels was employed to investigate the reaction pathways for the aminolysis reaction between succinic anhydride and methylamine. The single point ab initio calculations were based on the second-order M?ller-Plesset perturbation theory (MP2) with 6-31G(d) and 6-311++G(d,p) basis sets and CCSD(T)/6-31G(d) level calculations for geometries optimized at the B3LYP/6-311++G(d,p) level of theory. A detailed analysis of the atomic movements during the process of concerted aminolysis was further obtained by intrinsic reaction coordinate calculations. Solvent effects were assessed by the polarized continuum model method. The results show that the concerted mechanism of noncatalyzed aminolysis has distinctly lower activation energy compared with the addition/elimination stepwise mechanism. In the case of the process catalyzed by a second methylamine molecule, asynchronous proton transfer takes place, while the transition vectors of the acid-catalyzed transition states correspond to the simultaneous motion of protons. The most favorable pathway of the reaction was found through the bifunctional acid catalyzed stepwise mechanism that involves formation of eight-membered rings in the transition state structures. The difference between the activation barriers for the two mechanisms averages 2 kcal/mol at various levels of theory.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Isomerization of B_6,B_6~- and B_6~+ clusters

The interconversions between isomers with the same spin multiplicity of neutral B6 and charged B6-and B6+ clusters have been investigated at the B3LYP/6-311+G level of theory,including determination of the minimum energy pathways with transition states connecting the corresponding reactants and products.In dynamic calculations,26 isomers were optimized,including 11 novel isomers.In order to further refine the energies,single-point B3LYP/6-311+G(3df) calculations were carried out on the corresponding B3LYP/6...     

O（~1D）与CF_2HCl反应的理论研究

用量子化学密度泛函理论（DFT）和G3B3方法，对为（~1D）与CF_2HCl的反应 进行了研究，在B3LYP/6-311+G(d)，B3LYP/6-311+G(2df,2pd)和G3B3计算水平上， 优化了反应热能面上各驻点的几何结构，通过内禀反应坐标（IRC）计算和振动分 析，对反应过渡态进行了确认，并确定了反应机理。     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.