ONIOM meets : efficient, accurate, and robust multi-layer simulations across the periodic table

The computational treatment of large molecular structures is of increasing interest in fields of modern chemistry. Accordingly, efficient quantum chemical approaches are needed to perform sophisticated investigations on such systems. This engaged the development of the well-established "Our own N-layered integrated molecular orbital and molecular mechanics" (ONIOM) multi-layer scheme [L. W. Chung et al., Chem. Rev., 2015, 115, 5678-5796]. In this work, we present the specific implementation of the ONIOM scheme into the xtb semi-empirical extended tight-binding program package and its application to challenging transition-metal complexes. The efficient and broadly applicable GFNn-xTB and -FF methods are applied in the ONIOM framework to elucidate reaction energies, geometry optimizations, and explicit solvation effects for metal-organic systems with up to several hundreds of atoms. It is shown that an ONIOM-based combination of density functional theory, semi-empirical, and force-field methods can be used to drastically reduce the computational costs and thus enable the investigation of huge systems at almost no significant loss in accuracy.

Elucidating the Role of O2 Uncoupling in the Oxidative Biodegradation of Organic Contaminants by Rieske Non-heme Iron Dioxygenases

Oxygenations of aromatic soil and water contaminants with molecular O₂ catalyzed by Rieske dioxygenases are frequent initial steps of biodegradation in natural and engineered environments. Many of these non-heme ferrous iron enzymes are known to be involved in contaminant metabolism, but the understanding of enzyme–substrate interactions that lead to successful biodegradation is still elusive. Here, we studied the mechanisms of O₂ activation and substrate hydroxylation of two nitroarene dioxygenases to evaluate enzyme- and substrate-specific factors that determine the efficiency of oxygenated product formation. Experiments in enzyme assays of 2-nitrotoluene dioxygenase (2NTDO) and nitrobenzene dioxygenase (NBDO) with methyl-, fluoro-, chloro-, and hydroxy-substituted nitroaromatic substrates reveal that typically 20–100% of the enzyme’s activity involves unproductive paths of O₂ activation with generation of reactive oxygen species through so-called O₂ uncoupling. The ¹⁸O and ¹³C kinetic isotope effects of O₂ activation and nitroaromatic substrate hydroxylation, respectively, suggest that O₂ uncoupling occurs after generation of Fe^III-(hydro)peroxo species in the catalytic cycle. While 2NTDO hydroxylates ortho-substituted nitroaromatic substrates more efficiently, NBDO favors meta-substituted, presumably due to distinct active site residues of the two enzymes. Our data implies, however, that the O₂ uncoupling and hydroxylation activity cannot be assessed from simple structure–reactivity relationships. By quantifying O₂ uncoupling by Rieske dioxygenases, our work provides a mechanistic link between contaminant biodegradation, the generation of reactive oxygen species, and possible adaptation strategies of microorganisms to the exposure of new contaminants.

NRVS Studies of the Peroxide Shunt Intermediate in a Rieske Dioxygenase and Its Relation to the Native FeII O2 Reaction

The Rieske dioxygenases are a major subclass of mononuclear nonheme iron enzymes that play an important role in bioremediation. Recently, a high-spin FeIII-(hydro)peroxy intermediate (BZDOp) has been trapped in the peroxide shunt reaction of benzoate 1,2-dioxygenase. Defining the structure of this intermediate is essential to understanding the reactivity of these enzymes. Nuclear resonance vibrational spectroscopy (NRVS) is a recently developed synchrotron technique that is ideal for obtaining vibrational, and thus structural, information on Fe sites, as it gives complete information on all vibrational normal modes containing Fe displacement. In this study, we present NRVS data on BZDOp and assign its structure using these data coupled to experimentally calibrated density functional theory calculations. From this NRVS structure, we define the mechanism for the peroxide shunt reaction. The relevance of the peroxide shunt to the native FeII/O2 reaction is evaluated. For the native FeII/O2 reaction, an FeIII-superoxo intermediate is found to react directly with substrate. This process, while uphill thermodynamically, is found to be driven by the highly favorable thermodynamics of proton-coupled electron transfer with an electron provided by the Rieske [2Fe-2S] center at a later step in the reaction. These results offer important insight into the relative reactivities of FeIII-superoxo and FeIII-hydroperoxo species in nonheme Fe biochemistry.

QM/MM Study of the Reaction Catalyzed by Alkyladenine DNA Glycosylase: Examination of the Substrate Specificity of a DNA Repair Enzyme

Human alkyladenine DNA glycosylase (AAG) functions as part of the base excision repair pathway to excise structurally diverse oxidized and alkylated DNA purines. Specifically, AAG uses a water molecule activated by a general base and a nonspecific active site lined with aromatic residues to cleave the N-glycosidic bond. Despite broad substrate specificity, AAG does not target the natural purines (adenine (A) and guanine (G)). Using the ONIOM(QM:MM) methodology, we provide fundamental atomic level details of AAG bound to DNA-containing a neutral substrate (hypoxanthine (Hx)), a nonsubstrate (G), or a cationic substrate (7-methylguanine (7MeG)) and probe changes in the reaction pathway that occur when AAG targets different nucleotides. We reveal that subtle differences in protein-DNA contacts upon binding different substrates within the flexible AAG active site can significantly affect the deglycosylation reaction. Notably, we predict that AAG excises Hx in a concerted mechanism that is facilitated through correct alignment of the (E125) general base due to hydrogen bonding with a neighboring aromatic amino acid (Y127). Hx departure is further stabilized by π-π interactions with aromatic amino acids and hydrogen bonds with active site water. Despite possessing a similar structure to Hx, G is not excised since the additional exocyclic amino group leads to misalignment of the general base due to disruption of the key E125-Y127 hydrogen bond, the catalytically unfavorable placement of water within the active site, and weakened π-contacts between aromatic amino acids and the nucleobase. In contrast, cationic 7MeG does not occupy the same position within the AAG active site as G due to steric clashes with the additional N7 methyl group, which results in the correct alignment of the general base and permits nucleobase excision as observed for neutral Hx. Overall, our structural data rationalizes the observed substrate specificity of AAG and contributes to our fundamental understanding of enzymes with flexible active sites and broad substrate specificities.

Resolving Discrepancy between Theory and Experiment in 4-Nitrotoluene Oxidation

We have performed calculations of possible oxidation pathways of 4-nitrotoluene (4NT) by permanganate anion and evaluated relative contributions of oxidation of the methyl group and aromatic ring. Although a few theory levels matched the experimental results obtained by compound specific isotope analysis (CSIA) for 4NT, they failed in reproducing results for other nitrotoluene derivatives studied previously [Wijker, R.S.; Adamczyk, P.; Bolotin, J.; Paneth, P.; Hofstetter, T.B. Environ. Sci. Technol., 2013, 47, 13459-13468]. This discrepancy prompted us to reevaluate the experimental isotopic fractionation of carbon and hydrogen for 4NT on which the relative contributions of the oxidation channels has been based. Using position specific isotope analysis (PSIA) for hydrogen isotopic fractionation we have found that the previously determined value was incorrect. Reexamination of theory levels that are in agreement with these new findings indicated that while a better agreement for this particular case can be reached, overall, the previously used B3LYP functional expressed in the 6-31+G(d,p) basis set with inclusion of the polarized continuum model (PCM) of aqueous solution remains the theoretical level of choice in modeling oxidation of nitrotoluene derivatives by permanganate.

Insights into the Reaction Mechanism of Aromatic Ring Cleavage by Homogentisate Dioxygenase: A Quantum Mechanical/Molecular Mechanical Study

To elucidate the reaction mechanism of the ring cleavage of homogentisate by homogentisate dioxygenase, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out by using two systems in different protonation states of the substrate C2 hydroxyl group. When the substrate C2 hydroxyl group is ionized (the ionized pathway), the superoxo attack on the substrate is the rate-limiting step in the catalytic cycle, with a barrier of 15.9 kcal/mol. Glu396 was found to play an important role in stabilizing the bridge species and its O-O cleavage product by donating a proton via a hydrogen-bonded water molecule. When the substrate C2 hydroxyl group is not ionized (the nonionized pathway), the O-O bond cleavage of the bridge species is the rate-limiting step, with a barrier of 15.3 kcal/mol. The QM/MM-optimized geometries for the dioxygen and alkylperoxo complexes using the nonionized model (for the C2 hydroxyl group) are in agreement with the experimental crystal structures, suggesting that the C2 hydroxyl group is more likely to be nonionized.