James W. Gauld

A DFT Study of Nucleobase Dealkylation by the DNA Repair Enzyme AlkB

Oxidative dealkylation is a unique mechanistic pathway found in the alpha-ketoglutarate-Fe(II)-dependent AlkB family of enzymes to remove the alkylation damage to DNA bases and regenerate nucleobases to their native state. The B3LYP density functional combined with a self-consistent reaction field was used to explore the triplet, quintet, and septet spin-state potential energy surfaces of the multistep catalytic mechanism of AlkB. The mechanism was found to consist of four stages. First, binding of dioxygen to iron in the active-site complex occurs concerted with electron transfer, thereby yielding a ferric-superoxido species. Second, competing initiation for the activation of oxygen to generate the high-valent iron-oxygen intermediates (ferryl-oxo Fe(IV)O and ferric-oxyl Fe(III)O(*) species) was found to occur on the quintet and septet surfaces. Then, conformational reorientation of the activated iron-oxygen ligand was found to be nearly thermoneutral with a barrier of ca. 50 kJ mol(-1). The final stage is the oxidative dealkylation of the damaged nucleobase with the rate-controlling step being the abstraction of a hydrogen atom from the damaging methyl group by the ferryl-oxo ligand. For this step, the calculated barrier of 87.4 kJ mol(-1) is in good agreement with the experimental activation energy of ca. 83 kJ mol(-1) for the enzyme-catalyzed reaction.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is inactivated by S-sulfuration in vitro

Hydrogen sulfide (H2S) is produced enzymatically by cystathionine β-synthase (CBS) and cystathionine γ-lyase (CSE), as well as other enzymes in mammalian tissues. These discoveries have led to the crowning of H2S as yet another toxic gas that serves as a gasotransmitter like NO and CO. H2S is thought to exert its biological effects through its reaction with cysteine thiols in proteins, yielding sulfurated thiol (-SSH) derivatives. One of the first proteins shown to be modified by H2S was glyceraldehyde 3-phosphate dehydrogenase (GAPDH) [1] where the S-sulfuration of the active site cysteine (Cys 152) resulted in ~7-fold increase in the activity of the enzyme. In the present study we have attempted to reproduce this result with no success. GAPDH in its reduced, or hydrogen peroxide, or glutathione disulfide, or nitrosonium oxidized forms was reacted with sulfide or polysulfides. Sulfide had no effect on reduced GAPDH activity, while polysulfides inhibited GAPDH to ~42% of control. S-sulfuration of GAPDH occurred at Cys 247 after sulfide treatment, Cys 156 and Cys 247 after polysulfide treatment. No evidence of S-sulfuration at active site Cys 152 was discovered. Both sulfide and polysulfide was able to restore the activity of glutathione disulfide oxidized GAPDH, but not to control untreated levels. Treatment of glutathione disulfide oxidized GAPDH with polysulfide also produced S-sulfuration of Cys 156. Treatment of a C156S mutant of GAPDH with sulfide and polysulfide resulted in S-sulfuration of Cys 152, which also caused a decrease and not an increase in enzymatic activity. Computational chemistry shows S-sulfuration of Cys 156 may affect the position of catalytic Cys 152, raising its pKa by 0.5, which may affect the nucleophilicity of Cys 152. The current study raises significant questions about the reported ability of H2S to activate GAPDH by the sulfuration of its active site thiol, and indicates that polysulfide is a stronger protein S-sulfurating agent than sulfide.

Redox Mechanism of Glycosidic Bond Hydrolysis Catalyzed by 6-Phospho-α-glucosidase: A DFT Study

Glycosidic bonds are remarkably resistant to cleavage by chemical hydrolysis. Glycoside hydrolases catalyze their selective hydrolysis in oligosaccharides, polysaccharides, and glycoconjugates by following nonredox catalytic pathways or a net redox-neutral catalytic pathway using NAD(+) and divalent metal ions as cofactors. GlvA (6-phospho-alpha-glucosidase) is a glycosidase belonging to family GH4 and follows a regioselective redox-neutral mechanism of glycosidic-bond hydrolysis that favors alpha- over beta-glycosides. Its proposed catalytic mechanism can be divided into two half-reactions: the first one activates the glucopyranose ring by successively forming intermediates that are oxidized at the 3-, 2-, and 1-positions of the ring, which ultimately facilitate the heterolytic deglycosylation. The second half-reaction is essentially the reverse of the first half-reaction, beginning with the pyranose ring hydroxylation at the anomeric carbon, and it is followed by 3-reduction and regeneration of the active forms of the catalytic site and its cofactors. We investigated the NAD(+)-dependent redox mechanism of glycosidic bond hydrolysis as catalyzed by GlvA through the combined application of density functional theory and a self-consistent reaction field to a large active-site model obtained from the crystallographic structure of the enzyme, then we applied natural bond orbital and second-order perturbation analyses to monitor the electron flow and change in oxidation state on each atomic center along the reaction coordinate to rationalize the energetics and regioselectivity of this catalytic mechanism. We find that in GlvA, the redox catalytic mechanism of hydrolysis is driven by the gradual strengthening of the axial endo-anomeric component within the hexose ring along the reaction coordinate to facilitate the heterolytic dissociation of the axial C1-O bond. In addition, the combined influence of specific components of the generalized anomeric effect fully explains the regioselectivity observed in the catalytic activity of GlvA.

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.

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.

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 histidines 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.

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 hemes 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.

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.

A molecular dynamics and quantum mechanics/molecular mechanics study of the catalytic reductase mechanism of methionine sulfoxide reductase A: formation and reduction of a sulfenic acid.

The catalytic mechanism of MsrA in Mycobacterium tuberculosis, in which S-methionine sulfoxide (Met-O) is reduced to methionine (Met), has been investigated using docking, molecular dynamics (MD) simulations, and ONIOM (quantum mechanics/molecular mechanics) methods. In addition, the roles of specific active site residues, including an aspartyl (Asp87) near the recycling cysteine, tyrosyls (Tyr44 and Tyr92), and glutamyl (Glu52), have been examined, as well as the general effects of the protein and active site on the nature and properties of mechanistic intermediates. The mechanism is initiated by the transfer of a proton from the catalytic cysteines thiol (Cys13SH) via a bridging water to the R group carboxylate of Glu52. The now anionic sulfur of Cys13 nucleophilically attacks the substrates sulfur with concomitant transfer of a proton from Glu52 to the sulfoxide oxygen, generating a sulfurane. The active site enhances the proton affinity of the sulfurane oxygen, which can readily accept a proton from the phenolic hydroxyls of Tyr44 or Tyr92 to give a sulfonium cation. Subsequently, Asp87 and the recycling cysteine (Cys154) can facilitate nucleophilic attack of a solvent water at the Cys13S center of the sulfonium to give a sulfenic acid (Cys13SOH) and Met. For the subsequent reduction of Cys13SOH with intramolecular disulfide bond formation, Asp87 can help facilitate nucleophilic attack of Cys154S at the sulfur of Cys13SOH by deprotonating its thiol. This reduction is found likely to occur readily upon suitable positioning of the active site hydrogen bond network and the sulfur centers of both Cys13 and Cys154. The calculated rate-limiting barrier is in good agreement with experiment.

Computational Insights into the Mechanism of Porphobilinogen Synthase

Porphobilinogen synthase (PBGS) is a key enzyme in heme biosynthesis that catalyzes the formation of porphobilinogen (PBG) from two 5-aminolevulinic acid (5-ALA) molecules via formation of intersubstrate C-N and C-C bonds. The active site consists of several invariant residues, including two lysyl residues (Lys210 and Lys263; yeast numbering) that bind the two substrate moieties as Schiff bases. Based on experimental studies, various reaction mechanisms have been proposed for this enzyme that generally can be classified according to whether the intersubstrate C-C or C-N bond is formed first. However, the detailed catalytic mechanism of PBGS remains unclear. In the present study, we have employed density functional theory methods in combination with chemical models of the two key lysyl residues and two substrate moieties in order to investigate various proposed reaction steps and gain insight into the mechanism of PBGS. Importantly, it is found that mechanisms in which the intersubstrate C-N bond is formed first have a rate-limiting barrier (17.5 kcal/mol) that is lower than those in which the intersubstrate C-C bond is formed first (22.8 kcal/mol).