Computational Insights into the Regeneration of Ovothiol and Ergothioneine and Their Selenium Analogues by Glutathione.

Ovothiol and ergothioneine are powerful antioxidants that readily react with oxidants by forming their respective disulfides. In fact, ovothiol is widely considered one of the most powerful natural antioxidants. However, for these antioxidants to be again involved in reacting with oxidants, they must be regenerated via the reduction of the disulfide bonds. In the present work, the regeneration of the antioxidants ovothiol and ergothioneine and their selenium analogues, by the closed-shell nucleophilic attack of glutathione, was investigated using density functional theory. From the calculated thermodynamic data, the attack of glutathione on OSSO and EYYE (where Y = S and/or Se) will readily occur in solution. Moreover, in comparison to the reference reaction GSH + GSSG → GSSG + GSH, all reactions are expected to be faster. Overall, the results presented herein show that the key antioxidant GSH should readily recycle ovothiol, ovoselenol, ergothioneine, and ergoseloneine from OYYO and EYYE (where Y = S and/or Se).

Density Functional Theory Investigation of As(III) S-Adenosylmethionine Methyltransferase.

Arsenic is one of the most pervasive environmental toxins. It enters our water and food supply through many different routes, including the burning of fossil fuels, the application of arsenic-based herbicides, and natural sources. Using a density functional theory (DFT) cluster approach, the mechanism of arsenic (III) S-adenosylmethionine methyltransferases and various selenium-containing analogues was investigated. Notably, the methylation of arsenic by arsenic (III) S-adenosylmethionine is proposed to be a way to remove arsenic from contaminated water or soil. From the DFT cluster results, it was found that the selective substitution of the active-site Cys44, Cys72, and Cys174 residues with selenocysteines had a marginal effect on the barrier for CH3 transfer. Specifically, the average Gibbs activation energy was calculated to be only 4.2 kJ mol-1 lower than the Gibbs activation energy of 107.4 kJ mol-1 for the WT enzyme. However, importantly, it was found that with selective mutation, the methylation process becomes considerably more exergonic, where the methylation reaction can be made to be 26.4 kJ mol-1 more exergonic than the reaction catalyzed by the WT enzyme. Therefore, we propose that the selective substitution of the active-site Cys44, Cys72 and Cys174 residues with selenocysteines could make the process of methylation and volatilization more advantageous for bioremediation.

Examination of sulfonamide-based inhibitors of MMP3 using the conditioned media of invasive glioma cells.

Glioblastoma multiforme (GBM) is the deadliest and the most common primary malignant brain tumour. The median survival for patients with GBM is around one year due to the nature of glioma cells to diffusely invade that make the complete surgical resection of tumours difficult. Based upon the connexin43 (Cx43) model of glioma migration we have developed a computational framework to evaluate MMP inhibition in materials relevant to GBM. Using the ilomastat Leu-Trp backbone, we have synthesised novel sulphonamides and monitored the performance of these compounds in conditioned media expressing MMP3. From the results discussed herein we demonstrate the performance of sulfonamide based MMPIs included AP-3, AP-6, and AP-7.

Generation and Reactions of a Benzodehydrotropylium Ion-Co2(CO)6 Complex.

A series of 7-methylenedehydrobenzo[7]annulen-5-ol hexacarbonyldicobalt complexes were generated by Hosomi-Sakurai reactions of allylsilanes containing o-alkynylarylaldehyde-Co2(CO)6 complexes. One of the cyclization products was converted into its corresponding dihydrobenzo[7]annulen-7-ol hexacarbonyldicobalt complex, an immediate precursor to a benzodehydrotropylium-Co2(CO)6. The cation was generated in situ and reacted with four nucleophiles, and its aromatic stabilization was determined by computational methods.

Computational Investigation into the Ni(SeNHC2(CN)2)2 and Ni(SNHC2(CN)2)2 Complexes as Potential Catalysts for Hydrogen Production.

To reduce our carbon footprint, we must look at alternative non-carbon-containing fuels to prevent continued global climate change. One environmentally friendly alternative fuel is molecular hydrogen. Herein the Ni(SeNHC2(CN)2)2 complex was studied using DFT to determine the thermodynamics associated with the electrocatalytic formation of H2(g). From the calculated thermodynamics, it appears that the Ni(SeNHC2(CN)2)2 complex is predicted to catalyze the production of H2 gas under mildly reducing conditions relative to the SHE. Notably, the thermodynamics are better than the values calculated for the analogous Ni(SNHC2(CN)2)2 complex which has been shown experimentally to catalyze the formation of H2 gas in aqueous solution. Regarding possible kinetic reactivity, the HOMO-LUMO gap energies were calculated. From the gap energies, it is expected that the Se-containing compounds would be more reactive to electron transfer in the third reduction step, meaning therefore that a smaller overpotential would be needed to drive the reduction of Red2-H2 relative to SRed2-H2 in agreement with past experimental work. Thus, the use of Se in such compounds may offer a means to improve the catalysts for H2 production.

Assessment of several DFT functionals in calculation of the reduction potentials for Ni-, Pd-, and Pt-bis-ethylene-1,2-dithiolene and -diselenolene complexes.

We performed an assessment of 10 common DFT functionals to determine their suitability for calculating the reduction potentials of the ([M(S2C2H2)2](0)/[M(S2C2H2)2](1-)), ([M(Se2C2H2)2](0)/[M(Se2C2H2)2](1-)), ([M(S2C2H2)2](1-)/[M(S2C2H2)2](2-)), and ([M(Se2C2H2)2](1-)/[M(Se2C2H2)2](2-)) redox couples (M = Ni, Pd, and Pt). Overall it was found that the M06 functional leads to the best agreement with the gold standard CCSD(T) method with an average difference of only +0.07 V and a RMS of 0.07 V in calculated reduction potentials. The variability in calculated reduction potentials between the various DFT functionals arise, in part, from the multireference character of these systems, which was determined by the T1 diagnostic values. Thus, the bisdiselenolene complexes show similar multireference character as the bisdithiolene complexes, which were previously shown to have such character. In particular, for the Ni-, Pd-, and Pt-bisdiselenolene complexes the average T1 values are 0.05, 0.03, and 0.02, respectively. For the CCSD(T) calculations the similarities in the reduction potentials between analogous bisdithiolene and bisdiselenolene redox couples, which appear to be independent of the metal, is a result of the noninnocence of the dithiolene and diselenolene ligands. Thus, the reduction potential is more dependent on the ligand than the metal.

The catalytic formation of leukotriene C4: a critical step in inflammatory processes.

Leukotrienes (LT) are a family of drug-like molecules involved in the pathobiology of bronchial asthma and are responsible for smooth muscle contraction. Leukotriene C4 synthase (LTC4S) is a nuclear-membrane enzyme responsible for the conjugation of leukotriene A4 (LTA4) to glutathione to form LTC4, a cysteinyl leukotriene. In this study, the mechanism of LTA4 binding by LTC4S has been computationally examined. More specifically, docking and molecular dynamics simulations were used to gain insight into the substrate-bound active site. These studies identified two possible orientations for bound LTA4: 'tail-to-head' and 'head-to-tail'. An ONIOM(QM/MM) approach was then used to elucidate the mechanism by which glutathione may add to LTA4. In particular, the thiolate of glutathione acts as a nucleophile attacking C6 of LTA4 forming a S-C6 bond. Concomitantly, a proton is transferred from the guanidinium of Arg31 to the epoxide ring oxygen. This results in opening of the epoxide ring and stabilization of the LTC4 product complex. Within the present computational methodology the 'tail-to-head' orientation appears to be the most likely substrate orientation.

The one-electron oxidation of a dithiolate molecule: the importance of chemical intuition.

A series of nine commonly used density functional methods were assessed to accurately predict the oxidation potential of the (C2H2S2(-2)/C2H2S2(•-)) redox couple. It was found that due to their greater tendency for charge delocalization the GGA functionals predict a structure where the radical electron is delocalized within the alkene backbone of C2H2S2(•-), whereas the hybrid functionals and the reference QCISD/cc-pVTZ predict that the radical electron remains localized on the sulfurs. However, chemical intuition suggests that the results obtained with the GGA functionals should be correct. Indeed, with the use of the geometries obtained at the HCTH/6-311++G(3df,3pd) level of theory both the QCISD and hybrid DFT methods yield a molecule with a delocalized electron. Notably, this new molecule lies at least 53 kJ mol(-1) lower in energy than the previously optimized one that had a localized radical. Using these new structures the calculated oxidation potential was found to be 2.71-2.97 V for the nine DFT functionals tested. The M06-L functional provided the best agreement with the QCISD/cc-pVTZ reference oxidation potential of 3.28 V.

The one-electron reduction of dithiolate and diselenolate ligands.

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

A molecular dynamics examination on mutation-induced catalase activity in coral allene oxide synthase.

Coral allene oxide synthase (cAOS) catalyzes the formation of allene oxides from fatty acid hydroperoxides. Interestingly, its active site differs from that of catalase by only a single residue yet is incapable of catalase activity. That is, it is unable to catalyze the decomposition of hydrogen peroxide to molecular oxygen and water. However, the single active-site mutation T66V allows cAOS to exhibit catalase activity. We have performed a series of molecular dynamics (MD) simulations in order to gain insights into the differences in substrate (8R-hydroperoxyeicosatetraenoic) and H2O2 active site binding between wild-type cAOS and the T66V mutant cAOS. It is observed that in wild-type cAOS the active site Thr66 residue consistently forms a strong hydrogen-bonding interaction with H2O2 (catalase substrate) and, importantly, with the aid of His67 helps to pull H2O2 away from the heme Fe center. In contrast, in the T66V-cAOS mutant the H2O2 is much closer to the heme's Fe center and now forms a consistent Fe···O2H2 interaction. In addition, the His67···H2O2 distance shortens considerably, increasing the likelihood of a Cpd I intermediate and hence exhibiting catalase activity.

Insights into the catalytic mechanism of coral allene oxide synthase: a dispersion corrected density functional theory study.

In this present work the mechanism by which cAOS catalyzes the formation of allene oxide from its hydroperoxy substrate was computationally investigated by using a DFT-chemical cluster approach. In particular, the effects of dispersion interactions and DFT functional choice (M06, B3LYP, B3LYP*, and BP86), as well as the roles of multistate reactivity and the tyrosyl proximal ligand, were examined. It is observed that the computed relative free energies of stationary points along the overall pathway are sensitive to the choice of DFT functional, while the mechanism obtained is generally not. Large reductions in relative free energies for stationary points along the pathway (compared to the initial reactant complex) of on average 46.3 and 97.3 kJ mol(-1) for the doublet and quartet states, respectively, are observed upon going from the M06 to BP86 functional. From results obtained by using the B3LYP* method, well-tested previously on heme-containing systems, the mechanism of cAOS appears to occur with considerably higher Gibbs free energies than that for the analogous pathway in pAOS, possibly due to the presence of a ligating tyrosyl residue in cAOS. Furthermore, at the IEFPCM-B3LYP*/6-311+G(2df,p)//B3LYP/BS1 level of theory the inclusion of dispersion effects leads to the suggestion that the overall mechanism of cAOS could occur without the need for spin inversion.

A density functional theory investigation into the binding of the antioxidants ergothioneine and ovothiol to copper.

The ability of hybrid, nonhybrid and meta-GGA density functional theory (DFT) based methods (B3LYP, BP86, M06 and M06L) to provide reliable structures and thermochemical properties of biochemically important Cu(I)/(II)···ESH (ergothioneine) and ···OSH (ovothiol) has been assessed. For all functionals considered, convergence in the optimized structures and Cu(I)/(II)···S/N bond lengths is only obtained using the 6-311+G(2df,p) basis set or larger, with the nonhybrid DFT method BP86 appearing, in general, to provide the most reliable structures. The reduction potentials associated with the reduction of Cu(II) to Cu(I) when complexed with either OSH and ESH were also determined. The implications for their ability to thus help protect against Cu-mediated oxidative damage are discussed. Importantly, the binding of OSH and ESH with Cu ions disfavors Cu(I)/Cu(II) recycling by increasing the reduction potential for the Cu(II) to Cu(I) reduction and as a result, inhibits the potential oxidative damage caused by such Cu ions.

Gaining insight into the chemistry of lipoxygenases: a computational investigation into the catalytic mechanism of (8R)-lipoxygenase.

Lipoxygenases (LOXs) are ubiquitous in nature and catalyze a range of life-essential reactions within organisms. In particular they are critical to the formation of eicosanoids, which are critical for normal cell function. However, a number of important questions about the reactivity and mechanism of these enzymes still remain. Specifically, although the initial step in the mechanism of LOXs has been well studied, little is known of subsequent steps. Thus, with use of a quantum mechanical/molecular mechanical approach, the complete catalytic mechanism of (8R)-LOX was investigated. The results have provided a better understanding of the general chemistry of LOXs as a whole. In particular, from comparisons with soybean LOX-1, it appears that the initial proton-coupled electron transfer may be very similar among all LOXs. Furthermore, LOXs appear to undergo multistate reactivity where potential spin inversion of an electron may occur either in the attack of O(2) or in the regeneration of the active site. Lastly, it is shown that with the explicit modeling of the environment, the regeneration of the active center likely occurs via the rotation of the intermediate followed by an outer-sphere [Formula: see text] transfer as opposed to the formation of a "purple intermediate" complex.

An assessment of pure, hybrid, meta, and hybrid-meta GGA density functional theory methods for open-shell systems: the case of the nonheme iron enzyme 8R-LOX.

The performance of a range density functional theory functionals combined in a quantum mechanical (QM)/molecular mechanical (MM) approach was investigated in their ability to reliably provide geometries, electronic distributions, and relative energies of a multicentered open-shell mechanistic intermediate in the mechanism 8R-Lipoxygenase. With the use of large QM/MM active site chemical models, the smallest average differences in geometries between the catalytically relevant quartet and sextet complexes were obtained with the B3LYP(*) functional. Moreover, in the case of the relative energies between (4) II and (6) II, the use of the B3LYP(*) functional provided a difference of 0.0 kcal mol(-1). However, B3LYP(±) and B3LYP also predicted differences in energies of less than 1 kcal mol(-1). In the case of describing the electronic distribution (i.e., spin density), the B3LYP(*), B3LYP, or M06-L functionals appeared to be the most suitable. Overall, the results obtained suggest that for systems with multiple centers having unpaired electrons, the B3LYP(*) appears most well rounded to provide reliable geometries, electronic structures, and relative energies.

A Molecular Dynamics (MD) and Quantum Mechanics/Molecular Mechanics (QM/MM) study on Ornithine Cyclodeaminase (OCD): a tale of two iminiums.

Ornithine cyclodeaminase (OCD) is an NAD+-dependent deaminase that is found in bacterial species such as Pseudomonas putida. Importantly, it catalyzes the direct conversion of the amino acid L-ornithine to L-proline. Using molecular dynamics (MD) and a hybrid quantum mechanics/molecular mechanics (QM/MM) method in the ONIOM formalism, the catalytic mechanism of OCD has been examined. The rate limiting step is calculated to be the initial step in the overall mechanism: hydride transfer from the L-ornithine's C(α)-H group to the NAD+ cofactor with concomitant formation of a C(α)=NH(2)+ Schiff base with a barrier of 90.6 kJ mol-1. Importantly, no water is observed within the active site during the MD simulations suitably positioned to hydrolyze the C(α)=NH(2)+ intermediate to form the corresponding carbonyl. Instead, the reaction proceeds via a non-hydrolytic mechanism involving direct nucleophilic attack of the δ-amine at the C(α)-position. This is then followed by cleavage and loss of the α-NH(2) group to give the Δ1-pyrroline-2-carboxylate that is subsequently reduced to L-proline.

Model iron-oxo species and the oxidation of imidazole: insights into the mechanism of OvoA and EgtB?

A density functional theory cluster and first-principles quantum and statistical mechanics approach have been used to investigate the ability of iron-oxygen intermediates to oxidize a histidine cosubstrate, which may then allow for the possible formation of 2- and 5-histidylcysteine sulfoxide, respectively. Namely, the ability of ferric superoxo (Fe(III)O(2)(•-)), Fe(IV)═O, and ferrous peroxysulfur (Fe(III)OOS) complexes to oxidize the imidazole of histidine via an electron transfer (ET) or a proton-coupled electron transfer (PCET) was considered. While the high-valent mononuclear Fe(IV)═O species is generally considered the ultimate biooxidant, the free energies for its reduction (via ET or PCET) suggest that it is unable to directly oxidize histidine's imidazole. Instead, only the ferrous peroxysulfur complexes are sufficiently powerful enough oxidants to generate a histidyl-derived radical via a PCET process. Furthermore, while this process preferably forms a HisN(δ)(-H)(•) radical, several such oxidants are also suggested to be capable of generating the higher-energy HisC(δ)(-H)(•) and HisC(ε)(-H)(•) radicals. Importantly, the present results suggest that formation of the sulfoxide-containing products (seen in both OvoA and EgtB) is a consequence of the reduction of a powerful Fe(III)OOS oxidant via a PCET.

Molecular dynamics investigation into substrate binding and identity of the catalytic base in the mechanism of Threonyl-tRNA synthetase.

The structure and nature of the fully bound active site of Threonyl-tRNA Synthetase (ThrRS) for the second half-reaction has been investigated using molecular dynamics simulations. More specifically, we examined the ThrRS active site with both the substrate Threonyl-AMP and the cosubstrate cognate Threonyl-tRNA bound. Furthermore, we also considered the cases in which an active-site histidyl residue (His309) is either neutral or protonated. Moreover, we considered the role a water molecule may play in formation of a viable Michaelis complex. From the results it is found that the most likely role of His309 is in binding and properly orientating the ribose of the Ado76 nucleotidyl residue of the threonyl-tRNA via formation of a direct His309···Ado76 hydrogen bond, i.e., without involvement of a water. In addition, the imidazole of the His309 residue is likely neutral. It was found that upon protonation the positioning of the Ado76-3'-OH was perturbed, leading to a reduced chance for nucleophilic attack of the threonyl's C1 center.

A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase.

Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the α-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and α-aminoadipate-δ-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, α-aminoadipate-δ-semialdehyde and NADPH in a manner and polar environment conducive to reaction.

The α-amino group of the threonine substrate as the general base during tRNA aminoacylation: a new version of substrate-assisted catalysis predicted by hybrid DFT.

Density functional theory-based methods in combination with large chemical models have been used to investigate the mechanism of the second half-reaction catalyzed by Thr-tRNA synthetase: aminoacyl transfer from Thr-AMP onto the (A76)3'OH of the cognate tRNA. In particular, we have examined pathways in which an active site His309 residue is either protonated or neutral (i.e., potentially able to act as a base). In the protonated His309-assisted mechanism, the rate-limiting step is formation of the tetrahedral intermediate. The barrier for this step is 155.0 kJ mol(-1), and thus, such a pathway is concluded to not be enzymatically feasible. For the neutral His309-assisted mechanism, two models were used with the difference being whether Lys465 was included. For either model, the barrier of the rate-limiting step is below the upper thermodynamic enzymatic limit of ~125 kJ mol(-1). Specifically, without Lys465, the rate-limiting barrier is 122.1 kJ mol(-1) and corresponds to a rotation about the tetrahedral intermediate C(carb)-OH bond. For the model with Lys465, the rate-limiting barrier is slightly lower and corresponds to the formation of the tetrahedral intermediate. Importantly, for both "neutral His309" models, the neutral amino group of the threonyl substrate directly acts as the proton acceptor; in the formation of the tetrahedral intermediate, the (A76)3'OH proton is directly transferred onto the Thr-NH(2). Therefore, the overall mechanism follows a general substrate-assisted catalytic mechanism.

The first branching point in porphyrin biosynthesis: a systematic docking, molecular dynamics and quantum mechanical/molecular mechanical study of substrate binding and mechanism of uroporphyrinogen-III decarboxylase.

In humans, uroporphyrinogen decarboxylase is intimately involved in the synthesis of heme, where the decarboxylation of the uroporphyrinogen-III occurs in a single catalytic site. Several variants of the mechanistic proposal exist; however, the exact mechanism is still debated. Thus, using an ONIOM quantum mechanical/molecular mechanical approach, the mechanism by which uroporphyrinogen decarboxylase decarboxylates ring D of uroporphyrinogen-III has been investigated. From the study performed, it was found that both Arg37 and Arg50 are essential in the decarboxylation of ring D, where experimentally both have been shown to be critical to the catalytic behavior of the enzyme. Overall, the reaction was found to have a barrier of 10.3 kcal mol(-1) at 298.15 K. The rate-limiting step was found to be the initial proton transfer from Arg37 to the substrate before the decarboxylation. In addition, it has been found that several key interactions exist between the substrate carboxylate groups and backbone amides of various active site residues as well as several other functional groups.