Synergistic Substrate and Oxygen Activation in Salicylate Dioxygenase Revealed by QM/MM Simulations

Proton Transfer via Arginine with Suppressed pKa Mediates Catalysis by Gentisate and Salicylate Dioxygenase

Gentisate and salicylate 1,2-dioxygenases (GDO and SDO) facilitate aerobic degradation of aromatic rings by inserting both atoms of dioxygen into their substrates, thereby participating in global carbon cycling. The role of acid-base catalysts in the reaction cycles of these enzymes is debatable. We present evidence of the participation of a proton shuffler during catalysis by GDO and SDO. The pH dependence of Michaelis-Menten parameters demonstrates that a single proton transfer is mandatory for the catalysis. Measurements at variable temperatures and pHs were used to determine the standard enthalpy of ionization (ΔHion°) of 51 kJ/mol for the proton transfer event. Although the observed apparent pKa in the range of 6.0-7.0 for substrates of both enzymes is highly suggestive of a histidine residue, ΔHion° establishes an arginine residue as the likely proton source, providing phylogenetic relevance for this strictly conserved residue in the GDO family. We propose that the atypical 3-histidine ferrous binding scaffold of GDOs contributes to the suppression of arginine pKa and provides support for this argument by employing a 2-histidine-1-carboxylate variant of the enzyme that exhibits elevated pKa. A reaction mechanism considering the role of the proton source in stabilizing key reaction intermediates is proposed.

Oxygen Activation in Aromatic Ring Cleaving Salicylate Dioxygenase: Detection of Reaction Intermediates with a Nitro-substituted Substrate Analog

Cupin dioxygenases such as salicylate 1,2-dioxygense (SDO) perform aromatic C-C bond scission via a 3-His motif tethered iron cofactor. Here, transient kinetics measurements are used to monitor the catalytic cycle of SDO by using a nitro-substituted substrate analog, 3-nitrogentisate. Compared to the natural substrate, the nitro group reduces the enzymatic kcat by 500-fold, thereby facilitating the detection and kinetic characterization of reaction intermediates. Sums and products of reciprocal relaxation times derived from kinetic measurements were found to be linearly dependent on O2 concentration, suggesting reversible formation of two distinct intermediates. Dioxygen binding to the metal cofactor takes place with a forward rate of 5.9×103 M-1 s-1: two orders of magnitude slower than other comparable ring-cleaving dioxygenses. Optical chromophore of the first intermediate is distinct from the in situ generated SDO Fe(III)-O2⋅- complex but closer to the enzyme-substrate precursor.

Decoding Regioselective Reaction Mechanism of the Gentisic Acid Catalyzed by Gentisate 1,2-Dioxygenase Enzyme

Gentisate 1,2-dioxygenase (GDO), a ring-fission non-heme dioxygenase enzyme, displays a unique regioselective reaction of gentisic acid (GTQ) in the presence of molecular oxygen. GTQ is an important intermediate in the...

N-Nitrosation Mechanism Catalyzed by Non-heme Iron-Containing Enzyme SznF Involving Intramolecular Oxidative Rearrangement

The non-heme iron-dependent enzyme SznF catalyzes a critical N-nitrosation step during the N-nitrosourea pharmacophore biosynthesis in streptozotocin. The intramolecular oxidative rearrangement process is known to proceed at the Fe^II-containing active site in the cupin domain of SznF, but its mechanism has not been elucidated to date. In this study, based on the density functional theory calculations, a unique mechanism was proposed for the N-nitrosation reaction catalyzed by SznF in which a four-electron oxidation process is accomplished through a series of complicated electron transferring between the iron center and substrate to bypass the high-valent Fe^IV═O species. In the catalytic reaction pathway, the O₂ binds to the iron center and attacks on the substrate to form the peroxo bridge intermediate by obtaining two electrons from the substrate exclusively. Then, instead of cleaving the peroxo bridge, the C^ε-N^ω bond of the substrate is homolytically cleaved first to form a carbocation intermediate, which polarizes the peroxo bridge and promotes its heterolysis. After O-O bond cleavage, the following reaction steps proceed effortlessly so that the N-nitrosation is accomplished without NO exchange among reaction species.

Mechanistic insights into a non-heme 2-oxoglutarate-dependent ethylene-forming enzyme: selectivity of ethylene-formation versusl-Arg hydroxylation

The ethylene-forming enzyme (EFE) is a unique member of the Fe(ii)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenases. It converts 2OG into ethylene plus three CO2 molecules (ethylene-forming reaction) and also catalyzes the C5 hydroxylation of l-arginine coupled to the oxidative decarboxylation of 2OG (l-Arg hydroxylation reaction). To uncover the mechanisms of the dual transformations by EFE, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out. Based on the results, a branched mechanism was proposed. An FeII-peroxysuccinate complex with a dissociated CO2 generated through the nucleophilic attack of the superoxo moiety of the Fe-O2 species on the keto carbon of 2OG is the key common intermediate in both reactions. A competition between the subsequent CO2 insertion (a key step in the ethylene-forming pathway) and the O-O bond cleavage (leading to the formation of succinate) governs the product selectivity. The calculated reaction barriers suggested that the CO2 insertion is favored over the O-O bond cleavage. This is consistent with the product preference observed in experiments. By comparison with the results of AsqJ (an Fe/2OG oxygenase that leads to substrate oxidation exclusively), the protein environment was found to be crucial for the selectivity. Further calculations demonstrated that the local electric field of the protein environment in EFE promotes ethylene formation by acting as a charge template, exemplifying the importance of the electrostatic interaction in enzyme catalysis. These findings offer mechanistic insights into the EFE catalysis and provide important clues for better understanding the unique ethylene-forming capability of EFE compared with other Fe/2OG oxygenases.

Substrate promiscuity and active site differences in gentisate 1,2-dioxygenases: electron paramagnetic resonance study

Gentisate 1,2-dioxygenases (GDOs) are non-heme iron enzymes that catalyze the oxidation of dihydroxylated aromatic substrate, gentisate (2,5-dihydroxybenzoate). Salicylate 1,2-dioxygenase (SDO), a member of the GDO family, performs the ring scission of monohydroxylated substrates such as salicylate, thereby oxidizing a broader range of substrates compared to GDOs. Although the two types of enzymes share a high degree of sequence similarity, the origin of substrate specificity between SDO and GDOs is not understood. We present electron paramagnetic resonance (EPR) investigation of ferrous-nitrosyl complexes of SDO and a GDO from the bacterium Corynebacterium glutamicum (GDO_Cg). The EPR spectra of these complexes, which mimic the Fe-substrate-O₂ intermediates in the catalytic cycle, show unexpected differences in the substrate binding mode and the coordination geometry of the metal cofactor in the two enzymes. Binding of substrate to the ferrous center increases the symmetry of the Fe(II)-NO complex in SDO, while a reverse trend is observed in GDO_Cg where substrate ligation reduces the symmetry of the nitrosyl complex. Identical EPR spectra were obtained for the NO derivatives of a variant of GDO_Cg(A112G), which can oxidize salicylate, and wild-type GDO_Cg revealing that the A112G mutation does not alter the nature of the Fe-substrate-O₂ ternary complex.

Mechanistic Insights into a Stibene Cleavage Oxygenase NOV1 from Quantum Mechanical/Molecular Mechanical Calculations

NOV1, a stilbene cleavage oxygenase, catalyzes the cleavage of the central double bond of stilbenes to two phenolic aldehydes, using a 4-His Fe(II) center and dioxygen. Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the reaction mechanism of the central double bond cleavage of phytoalexin resveratrol by NOV1. Our results showed that the oxygen molecule prefers to bind to the iron center in a side-on fashion, as suggested from the experiment. The quintet Fe-O2 complex with the side-on superoxo antiferromagnetic coupled to the resveratrol radical is identified as the reactive oxygen species. The QM/MM results support the dioxygenase mechanism involving a dioxetane intermediate with a rate-limiting barrier of 10.0 kcal mol-1. The alternative pathway through an epoxide intermediate is ruled out due to a larger rate-limiting barrier (26.8 kcal mol-1). These findings provide important insight into the catalytic mechanism of carotenoid cleavage oxygenases and also the dioxygen activation of non-heme enzymes.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Density Functional Theory Study on the Demethylation Reaction between Methylamine, Dimethylamine, Trimethylamine, and Tamoxifen Catalyzed by a Fe(IV)-Oxo Porphyrin Complex

In this work, we studied computationally the N-demethylation reaction of methylamine, dimethylamine, and trimethylamine as archetypal examples of primary, secondary, and tertiary amines catalyzed by high-field low-spin Fe-containing enzymes such as cytochromes P450. Using DFT calculations, we found that the expected C-H hydroxylation process was achieved for trimethylamine. When dimethylamine and methylamine were studied, two different reaction mechanisms (C-H hydroxylation and a double hydrogen atom transfer) were computed to be energetically accessible and both are equally preferred. Both processes led to the formation of formaldehyde and the N-demethylated substrate. Finally, as an illustrative example, the relative contribution of the three primary oxidation routes of tamoxifen was rationalized through energetic barriers obtained from density functional calculations and docking experiments involving CYP3A4 and CYP2D6 isoforms. We found that the N-demethylation process was the intrinsically favored one, whereas other oxidation reactions required most likely preorganization imposed by the residues close to the active sites.

Insights into the dioxygen activation and catalytic mechanism of the nickel-containing quercetinase

The catalytic mechanism of Ni-QueD^FLA was elucidated by QM/MM calculations, and the different reactivities of nickel and iron were illuminated.

Theoretical Study of the Mechanism of the Nonheme Iron Enzyme EgtB

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

Cleavage mechanism of the aliphatic C–C bond catalyzed by 2,4′-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP: a QM/MM study

Calculations suggest that the reactant complex may firstly undergo a triplet–quintet crossing to initiate the reaction and then the subsequent chemistry occurs on the multiple-states surfaces. The key C–C bond cleavage is accompanied by an insertion reaction of oxygen radical.

Expansion of the substrate range of the gentisate 1,2-dioxygenase from Corynebacterium glutamicum for the conversion of monohydroxylated benzoates

The gentisate 1,2-dioxygenases (GDOs) from Corynebacterium glutamicum and various other organisms oxidatively cleave the aromatic nucleus of gentisate (2,5-dihydroxybenzoate), but are not able to convert salicylate (2-hydroxybenzoate). In contrast, the α-proteobacterium Pseudaminobacter salicylatoxidans synthesises an enzyme (‘salicylate dioxygenase’, SDO) which cleaves gentisate, but also (substituted) salicylate(s). Sequence comparisons showed that the SDO belongs to a group of GDOs mainly originating from Gram-positive bacteria which also include the GDO from C. glutamicum ATCC 13032. The combination of sequence comparisons with previously performed structural and mutational analyses of the SDO allowed to identify an amino acid residue (Ala112) which might prevent the oxidation of (substituted) salicylate(s) by the GDO from C. glutamicum. Therefore, the relevant mutation (Ala→Gly) was introduced into the GDO from C. glutamicum. The GDO variant obtained gained the ability to oxidise salicylate and several other monohydroxylated substrates. In order to screen a broader range of enzyme variants a chromogenic assay was developed which allowed the detection of bacterial colonies converting salicylate. The applicability of this test system was proven by screening a set of GDO variants obtained by saturation mutagenesis at different positions. This demonstrated that also GDO variants carrying the mutations Ala112→Ser, Ala112→Ile and Ala112→Asp converted salicylate.