Mechanistic Insight into the Reactivities of Aqueous-Phase Singlet Oxygen with Organic Compounds

Singlet oxygen (¹O₂) is a selective reactive oxygen species that plays a key role for the fate of various organic compounds in the aquatic environment under sunlight irradiation, engineered water oxidation systems, atmospheric water droplets, and biomedical systems. While the initial rate-determining charge-transfer reaction mechanisms and kinetics of ¹O₂ have been studied extensively, no comprehensive studies have been performed to elucidate the reaction mechanisms with organic compounds that have various functional groups. In this study, we use density functional theory calculations to determine elementary reaction mechanisms with a wide variety of organic compounds. The theoretically calculated aqueous-phase free energies of activation of single electron transfer and ¹O₂ addition reactions are compared to the experimentally determined rate constants in the literature to determine linear free-energy relationships. The theoretically calculated free energies of activation for the groups of phenolates and phenols show excellent correlations with the Hammett constants that accept electron densities by through-resonance. The dominant elementary reaction mechanism is discussed for each group of compounds. As a practical implication, we demonstrate the fate of environmentally relevant organic compounds induced by photochemically produced intermediate species at different pH and evaluate the impact of predicting rate constants to the half-life.

Hammett Relationship in Oxidase-Mimicking Metal-Organic Frameworks Revealed through a Protein-Engineering-Inspired Strategy

While the unique physicochemical properties of nanomaterials that enable regulation of nanozyme activities are demonstrated in many systems, quantitative relationships between the nanomaterials structure and their enzymatic activities remain poorly understood, due to the heterogeneity of compositions and active sites in these nanomaterials. Here, inspired by metalloenzymes with well-defined metal-ligand coordination, a set of substituted metal-organic frameworks (MOFs) with similar coordination is employed to investigate the relationship between structure and oxidase-mimicking activity. Both experimental results and density functional theory calculations reveal a Hammett-type structure-activity linear free energy relationship (H-SALR) of MIL-53(Fe) (MIL = Materials of Institute Lavoisier) nanozymes, in which increasing the Hammett σm value with electron-withdrawing ligands increases the oxidase-mimicking activity. As a result, MIL-53(Fe) NO2 with the strongest electron-withdrawing NO2 substituent shows a tenfold higher activity than the unsubstituted MIL-53(Fe). Furthermore, the generality of H-SALR is demonstrated for a range of substrates, one other metal (Cr), and even one other MOF type (MIL-101). Such biologically inspired quantitative studies demonstrate that it is possible to identify quantitative structure-activity relationships of nanozymes, and to provide detailed insight into the catalytic mechanisms as those in native enzymes, making it possible to use these relationships to develop high-performance nanomaterials.

Editing Domain Motions Preorganize the Synthetic Active Site of Prolyl-tRNA Synthetase

Prolyl-tRNA synthetases (ProRSs) catalyze the covalent attachment of proline onto cognate tRNAs, an indispensable step for protein synthesis in all living organisms. ProRSs are modular enzymes and the "prokaryotic-like" ProRSs are distinguished from "eukaryotic-like" ProRSs by the presence of an editing domain (INS) inserted between motifs 2 and 3 of the main catalytic domain. Earlier studies suggested the presence of coupled-domain dynamics could contribute to catalysis; however, the role that the distal, highly mobile INS domain plays in catalysis at the synthetic active site is not completely understood. In the present study, a combination of theoretical and experimental approaches has been used to elucidate the precise role of INS domain dynamics. Quantum mechanical/molecular mechanical simulations were carried out to model catalytic Pro-AMP formation by Enterococcus faecalis ProRS. The energetics of the adenylate formation by the wild-type enzyme was computed and contrasted with variants containing active site mutations, as well as a deletion mutant lacking the INS domain. The combined results revealed that two distinct types of dynamics contribute to the enzyme's catalytic power. One set of motions is intrinsic to the INS domain and leads to conformational preorganization that is essential for catalysis. A second type of motion, stemming from the electrostatic reorganization of active site residues, impacts the height and width of the energy profile and has a critical role in fine tuning the substrate orientation to facilitate reactive collisions. Thus, motions in a distal domain can preorganize the active site of an enzyme to optimize catalysis.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

Remarkable Reactivity Differences between Glucosides with Identical Leaving Groups

Two isomeric aryl 2-deoxy-2-fluoro-β-glucosides react with a β-glucosidase at rates differing by 10⁶-fold, despite the fact that they release the same aromatic aglycone. In contrast, the equivalent glucoside substrates react with essentially identical rate constants. Insight into the source of these surprising rate differences was obtained through a comprehensive study of the nonenzymatic (spontaneous) hydrolysis of these same substrates, wherein an approximate 10⁵-fold difference in rates was measured, clarifying that the differences were inherent rather than being due to specific interactions with the enzyme. The possibility that an alternate nucleophilic aryl substitution mechanism was responsible for the rapid reaction of the faster substrate was excluded through ¹⁸O-labeling studies. Further exploration of the origins of these rate differences involved analysis of X-ray crystal structures as well as quantum chemical calculations, which surprisingly revealed that ground state destabilization and transition state stabilizing effects contribute almost equally to the observed reactivity differences. These studies highlight the dangers of using simple reference equilibria such as pK_a values as measures of leaving group ability.

Biological Chemistry of Hydrogen Selenide

There are no two main-group elements that exhibit more similar physical and chemical properties than sulfur and selenium. Nonetheless, Nature has deemed both essential for life and has found a way to exploit the subtle unique properties of selenium to include it in biochemistry despite its congener sulfur being 10,000 times more abundant. Selenium is more easily oxidized and it is kinetically more labile, so all selenium compounds could be considered to be "Reactive Selenium Compounds" relative to their sulfur analogues. What is furthermore remarkable is that one of the most reactive forms of selenium, hydrogen selenide (HSe- at physiologic pH), is proposed to be the starting point for the biosynthesis of selenium-containing molecules. This review contrasts the chemical properties of sulfur and selenium and critically assesses the role of hydrogen selenide in biological chemistry.

Substrate and Transition State Binding in Alkaline Phosphatase Analyzed by Computation of Oxygen Isotope Effects

Enzymes are powerful catalysts, and a thorough understanding of the sources of their catalytic power will facilitate many medical and industrial applications. Here we have studied the catalytic mechanism of alkaline phosphatase (AP), which is one of the most catalytically proficient enzymes known. We have used quantum mechanics calculations and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to model a variety of isotope effects relevant to the reaction of AP. We have calculated equilibrium isotope effects (EIEs), binding isotope effects (BIEs), and kinetic isotope effects (KIEs) for a range of phosphate mono- and diester substrates. The results agree well with experimental values, but the model for the reaction's transition state (TS) differs from the original interpretation of those experiments. Our model indicates that isotope effects on binding make important contributions to measured KIEs on V/K, which complicated interpretation of the measured values. Our results provide a detailed interpretation of the measured isotope effects and make predictions that can test the proposed model. The model indicates that the substrate is deformed in the ground state (GS) of the reaction and partially resembles the TS. The highly preorganized active site preferentially binds conformations that resemble the TS and not the GS, which induces the substrate to adapt to the enzyme, rather than the other way around-as with classic "induced fit" models. The preferential stabilization of the TS over the GS is what lowers the barrier to the chemical step.

Enantiodivergent Fluorination of Allylic Alcohols: Data Set Design Reveals Structural Interplay between Achiral Directing Group and Chiral Anion

Enantioselectivity values represent relative rate measurements that are sensitive to the structural features of the substrates and catalysts interacting to produce them. Therefore, well-designed enantioselectivity data sets are information rich and can provide key insights regarding specific molecular interactions. However, if the mechanism for enantioselection varies throughout a data set, these values cannot be easily compared. This premise, which is the crux of free energy relationships, exposes a challenging issue of identifying mechanistic breaks within multivariate correlations. Herein, we describe an approach to addressing this problem in the context of a chiral phosphoric acid catalyzed fluorination of allylic alcohols using aryl boronic acids as transient directing groups. By designing a data set in which both the phosphoric and boronic acid structures were systematically varied, key enantioselectivity outliers were identified and analyzed. A mechanistic study was executed to reveal the structural origins of these outliers, which was consistent with the presence of several mechanistic regimes within the data set. While 2- and 4-substituted aryl boronic acids favored the (R)-enantiomer with most of the studied catalysts, meta-alkoxy substituted aryl boronic acids resulted in the (S)-enantiomer when used in combination with certain (R)-phosphoric acids. We propose that this selectivity reversal is the result of a lone pair-π interaction between the substrate ligated boronic acid and the phosphate. On the basis of this proposal, a catalyst system was identified, capable of producing either enantiomer in high enantioselectivity (77% (R)-2 to 92% (S)-2) using the same chiral catalyst by subtly changing the structure of the achiral boronic acid.

Development of linear free energy relationships for aqueous phase radical-involved chemical reactions

Aqueous phase advanced oxidation processes (AOPs) produce hydroxyl radicals (HO•) which can completely oxidize electron rich organic compounds. The proper design and operation of AOPs require that we predict the formation and fate of the byproducts and their associated toxicity. Accordingly, there is a need to develop a first-principles kinetic model that can predict the dominant reaction pathways that potentially produce toxic byproducts. We have published some of our efforts on predicting the elementary reaction pathways and the HO• rate constants. Here we develop linear free energy relationships (LFERs) that predict the rate constants for aqueous phase radical reactions. The LFERs relate experimentally obtained kinetic rate constants to quantum mechanically calculated aqueous phase free energies of activation. The LFERs have been applied to 101 reactions, including (1) HO• addition to 15 aromatic compounds; (2) addition of molecular oxygen to 65 carbon-centered aliphatic and cyclohexadienyl radicals; (3) disproportionation of 10 peroxyl radicals, and (4) unimolecular decay of nine peroxyl radicals. The LFERs correlations predict the rate constants within a factor of 2 from the experimental values for HO• reactions and molecular oxygen addition, and a factor of 5 for peroxyl radical reactions. The LFERs and the elementary reaction pathways will enable us to predict the formation and initial fate of the byproducts in AOPs. Furthermore, our methodology can be applied to other environmental processes in which aqueous phase radical-involved reactions occur.

Probing the catalytic mechanism of bovine CD38/NAD+ glycohydrolase by site directed mutagenesis of key active site residues

Bovine CD38/NAD(+) glycohydrolase catalyzes the hydrolysis of NAD(+) to nicotinamide and ADP-ribose and the formation of cyclic ADP-ribose via a stepwise reaction mechanism. Our recent crystallographic study of its Michaelis complex and covalently-trapped intermediates provided insights into the modalities of substrate binding and the molecular mechanism of bCD38. The aim of the present work was to determine the precise role of key conserved active site residues (Trp118, Glu138, Asp147, Trp181 and Glu218) by focusing mainly on the cleavage of the nicotinamide-ribosyl bond. We analyzed the kinetic parameters of mutants of these residues which reside within the bCD38 subdomain in the vicinity of the scissile bond of bound NAD(+). To address the reaction mechanism we also performed chemical rescue experiments with neutral (methanol) and ionic (azide, formate) nucleophiles. The crucial role of Glu218, which orients the substrate for cleavage by interacting with the N-ribosyl 2'-OH group of NAD(+), was highlighted. This contribution to catalysis accounts for almost half of the reaction energy barrier. Other contributions can be ascribed notably to Glu138 and Asp147 via ground-state destabilization and desolvation in the vicinity of the scissile bond. Key interactions with Trp118 and Trp181 were also proven to stabilize the ribooxocarbenium ion-like transition state. Altogether we propose that, as an alternative to a covalent acylal reaction intermediate with Glu218, catalysis by bCD38 proceeds through the formation of a discrete and transient ribooxocarbenium intermediate which is stabilized within the active site mostly by electrostatic interactions.

Efficient, crosswise catalytic promiscuity among enzymes that catalyze phosphoryl transfer

The observation that one enzyme can accelerate several chemically distinct reactions was at one time surprising because the enormous efficiency of catalysis was often seen as inextricably linked to specialization for one reaction. Originally underreported, and considered a quirk rather than a fundamental property, enzyme promiscuity is now understood to be important as a springboard for adaptive evolution. Owing to the large number of promiscuous enzymes that have been identified over the last decade, and the increased appreciation for promiscuity's evolutionary importance, the focus of research has shifted to developing a better understanding of the mechanistic basis for promiscuity and the origins of tolerant or restrictive specificity. We review the evidence for widespread crosswise promiscuity amongst enzymes that catalyze phosphoryl transfer, including several members of the alkaline phosphatase superfamily, where large rate accelerations between 10(6) and 10(17) are observed for both native and multiple promiscuous reactions. This article is part of a Special Issue entitled: Chemistry and mechanism of phosphatases, diesterases and triesterases.

Linear free energy relationships between aqueous phase hydroxyl radical reaction rate constants and free energy of activation

The hydroxyl radical (HO(•)) is a strong oxidant that reacts with electron-rich sites on organic compounds and initiates complex radical chain reactions in aqueous phase advanced oxidation processes (AOPs). Computer based kinetic modeling requires a reaction pathway generator and predictions of associated reaction rate constants. Previously, we reported a reaction pathway generator that can enumerate the most important elementary reactions for aliphatic compounds. For the reaction rate constant predictor, we develop linear free energy relationships (LFERs) between aqueous phase literature-reported HO(•) reaction rate constants and theoretically calculated free energies of activation for H-atom abstraction from a C-H bond and HO(•) addition to alkenes. The theoretical method uses ab initio quantum mechanical calculations, Gaussian 1-3, for gas phase reactions and a solvation method, COSMO-RS theory, to estimate the impact of water. Theoretically calculated free energies of activation are found to be within approximately ±3 kcal/mol of experimental values. Considering errors that arise from quantum mechanical calculations and experiments, this should be within the acceptable errors. The established LFERs are used to predict the HO(•) reaction rate constants within a factor of 5 from the experimental values. This approach may be applied to other reaction mechanisms to establish a library of rate constant predictions for kinetic modeling of AOPs.

Glucosamine-6-sulfamate analogues of heparan sulfate as inhibitors of endosulfatases

Keeping Sulfate. The extracellular endosulfatases, which modulate signalling pathways by removing sulfate groups from heparan, can be inhibited by replacing the 6-sulfate destined for cleavage with an inhibitory sulfamate motif, as demonstrated by simple glucosamine-6-sulfamate analogs of heparan sulfate.