Using Computational Chemistry To Reveal Nature’s Blueprints for Single-Site Catalysis of C–H Activation

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

C-H Bond Oxidation by MnIV-Oxo Complexes: Hydrogen-Atom Tunneling and Multistate Reactivity

The reactivity of six MnIV-oxo complexes in C-H bond oxidation has been examined using a combination of kinetic experiments and computational methods. Variable-temperature studies of the oxidation of 9,10-dihydroanthracene (DHA) and ethylbenzene by these MnIV-oxo complexes yielded activation parameters suitable for evaluating electronic structure computations. Complementary kinetic experiments of the oxidation of deuterated DHA provided evidence for hydrogen-atom tunneling in C-H bond oxidation for all MnIV-oxo complexes. These results are in accordance with the Bell model, where tunneling occurs near the top of the transition-state barrier. Density functional theory (DFT) and DLPNO-CCSD(T1) computations were performed for three of the six MnIV-oxo complexes to probe a previously predicted multistate reactivity model. The DFT computations predicted a thermal crossing from the 4B1 ground state to a 4E state along the C-H bond oxidation reaction coordinate. DLPNO-CCSD(T1) calculations further confirm that the 4E transition state offers a lower energy barrier, reinforcing the multistate reactivity model for these complexes. We discuss how this multistate model can be reconciled with recent computations that revealed that the kinetics of C-H bond oxidation by this set of MnIV-oxo complexes can be well-predicted on the basis of the thermodynamic driving force for these reactions.

Investigating Ligand Sphere Perturbations on MnIII-Alkylperoxo Complexes

Manganese catalysts that activate hydrogen peroxide carry out several different hydrocarbon oxidation reactions with high stereoselectivity. The commonly proposed mechanism for these reactions involves a key manganese(III)-hydroperoxo intermediate, which decays via O-O bond heterolysis to generate a Mn(V)-oxo species that institutes substrate oxidation. Due to the scarcity of characterized MnIII-hydroperoxo complexes, MnIII-alkylperoxo complexes are employed to understand factors that affect the mechanism of the O-O cleavage. Herein, we report a series of novel complexes, including two room-temperature-stable MnIII-alkylperoxo species, supported by a new amide-containing pentadentate ligand (6Medpaq5NO2). We use a combination of spectroscopic methods and density functional theory computations to probe the effects of the electronic changes in the ligand sphere trans to the hydroxo and alkylperoxo units to thermal stability and reactivity. The structural characterizations for both MnII(OTf)(6Medpaq5NO2) and [MnIII(OH)(6Medpaq5NO2)](OTf) were obtained via single-crystal X-ray crystallography. A perturbation to the ligand sphere allowed for a marked increase in reactivity towards an organic substrate, a modest change in the distribution of the O-O cleavage products from homolytic and heterolytic pathways, and little change in thermal stability.

MetalHawk: Enhanced Classification of Metal Coordination Geometries by Artificial Neural Networks

The chemical properties of metal complexes are strongly dependent on the number and geometrical arrangement of ligands coordinated to the metal center. Existing methods for determining either coordination number or geometry rely on a trade-off between accuracy and computational costs, which hinders their application to the study of large structure data sets. Here, we propose MetalHawk (https://github.com/vrettasm/MetalHawk), a machine learning-based approach to perform simultaneous classification of metal site coordination number and geometry through artificial neural networks (ANNs), which were trained using the Cambridge Structural Database (CSD) and Metal Protein Data Bank (MetalPDB). We demonstrate that the CSD-trained model can be used to classify sites belonging to the most common coordination numbers and geometry classes with balanced accuracy equal to 96.51% for CSD-deposited metal sites. The CSD-trained model was also found to be capable of classifying bioinorganic metal sites from the MetalPDB database, with balanced accuracy equal to 84.29% on the whole PDB data set and to 91.66% on manually reviewed sites in the PDB validation set. Moreover, we report evidence that the output vectors of the CSD-trained model can be considered as a proxy indicator of metal-site distortions, showing that these can be interpreted as a low-dimensional representation of subtle geometrical features present in metal site structures.

Why Nonheme Iron Halogenases Do Not Fluorinate C-H Bonds: A Computational Investigation

Selective halogenation is necessary for a range of fine chemical applications, including the development of therapeutic drugs. While synthetic processes to achieve C-H halogenation require harsh conditions, enzymes such as nonheme iron halogenases carry out some types of C-H halogenation, i.e., chlorination or bromination, with ease, while others, i.e., fluorination, have never been observed in natural or engineered nonheme iron enzymes. Using density functional theory and correlated wave function theory, we investigate the differences in structural and energetic preferences of the smaller fluoride and the larger chloride or bromide intermediates throughout the catalytic cycle. Although we find that the energetics of rate-limiting hydrogen atom transfer are not strongly impacted by fluoride substitution, the higher barriers observed during the radical rebound reaction for fluoride relative to chloride and bromide contribute to the difficulty of C-H fluorination. We also investigate the possibility of isomerization playing a role in differences in reaction selectivity, and our calculations reveal crucial differences in terms of isomer energetics of the key ferryl intermediate between fluoride and chloride/bromide intermediates. While formation of monodentate isomers believed to be involved in selective catalysis is shown for chloride and bromide intermediates, we find that formation of the fluoride monodentate intermediate is not possible in our calculations, which lack additional stabilizing interactions with the greater protein environment. Furthermore, the shorter Fe-F bonds are found to increase isomerization reaction barriers, suggesting that incorporation of residues that form a halogen bond with F and elongate Fe-F bonds could make selective C-H fluorination possible in nonheme iron halogenases. Our work highlights the differences between the fluoride and chloride/bromide intermediates and suggests potential steps toward engineering nonheme iron halogenases to enable selective C-H fluorination.

Heterogenization of molecular catalysts within porous solids: the case of Ni-catalyzed ethylene oligomerization from zeolites to metal-organic frameworks

The last decade has seen a tremendous expansion of the field of heterogenized molecular catalysis, especially with the growing interest in metal-organic frameworks and related porous hybrid solids. With successful achievements in the transfer from molecular homogeneous catalysis to heterogenized processes come the necessary discussions on methodologies used and a critical assessment on the advantages of heterogenizing molecular catalysis. Here we use the example of nickel-catalyzed ethylene oligomerization, a reaction of both fundamental and applied interest, to review heterogenization methodologies of well-defined molecular catalysts within porous solids while addressing the biases in the comparison between original molecular systems and heterogenized counterparts.

Electronic structure contributions to O-O bond cleavage reactions for MnIII-alkylperoxo complexes

Synthetic manganese catalysts that activate hydrogen peroxide perform a variety of hydrocarbon oxidation reactions. The most commonly proposed mechanism for these catalysts involves the generation of a manganese(III)-hydroperoxo intermediate that decays via heterolytic O-O bond cleavage to generate a Mn(V)-oxo species that initiates substrate oxidation. Due to the paucity of well-defined MnIII-hydroperoxo complexes, MnIII-alkylperoxo complexes are often employed to understand the factors that affect the O-O cleavage reaction. Herein, we examine the decay pathways of the MnIII-alkylperoxo complexes [MnIII(OOtBu)(6Medpaq)]+ and [MnIII(OOtBu)(N4S)]+, which have distinct coordination environments (N5- and N4S-, respectively). Through the use of density functional theory (DFT) calculations and comparisons with published experimental data, we are able to rationalize the differences in the decay pathways of these complexes. For the [MnIII(OOtBu)(N4S)]+ system, O-O homolysis proceeds via a two-state mechanism that involves a crossing from the quintet reactant to a triplet state. A high energy singlet state discourages O-O heterolysis for this complex. In contrast, while quintet-triplet crossing is unfavorable for [MnIII(OOtBu)(6Medpaq)]+, a relatively low-energy single state accounts for the observation of both O-O homolysis and heterolysis products for this complex. The origins of these differences in decay pathways are linked to variations in the electronic structures of the MnIII-alkylperoxo complexes.

Computational Discovery of Stable Metal-Organic Frameworks for Methane-to-Methanol Catalysis

The challenge of direct partial oxidation of methane to methanol has motivated the targeted search of metal-organic frameworks (MOFs) as a promising class of materials for this transformation because of their site-isolated metals with tunable ligand environments. Thousands of MOFs have been synthesized, yet relatively few have been screened for their promise in methane conversion. We developed a high-throughput virtual screening workflow that identifies MOFs from a diverse space of experimental MOFs that have not been studied for catalysis, yet are thermally stable, synthesizable, and have promising unsaturated metal sites for C-H activation via a terminal metal-oxo species. We carried out density functional theory calculations of the radical rebound mechanism for methane-to-methanol conversion on models of the secondary building units (SBUs) from 87 selected MOFs. While we showed that oxo formation favorability decreases with increasing 3d filling, consistent with prior work, previously observed scaling relations between oxo formation and hydrogen atom transfer (HAT) are disrupted by the greater diversity in our MOF set. Accordingly, we focused on Mn MOFs, which favor oxo intermediates without disfavoring HAT or leading to high methanol release energies─a key feature for methane hydroxylation activity. We identified three Mn MOFs comprising unsaturated Mn centers bound to weak-field carboxylate ligands in planar or bent geometries with promising methane-to-methanol kinetics and thermodynamics. The energetic spans of these MOFs are indicative of promising turnover frequencies for methane to methanol that warrant further experimental catalytic studies.

Transition metal oxide complexes as molecular catalysts for selective methane to methanol transformation: any prospects or time to retire?

Transition metal oxides have been extensively used in the literature for the conversion of methane to methanol. Despite the progress made over the past decades, no method with satisfactory performance or economic viability has been detected. The main bottleneck is that the produced methanol oxidizes further due to its weaker C-H bond than that of methane. Every improvement in the efficiency of a catalyst to activate methane leads to reduction of the selectivity towards methanol. Is it therefore prudent to keep studying (both theoretically and experimentally) metal oxides as catalysts for the quantitative conversion of methane to methanol? This perspective focuses on molecular metal oxide complexes and suggests strategies to bypass the current bottlenecks with higher weight on the computational chemistry side. We first discuss the electronic structure of metal oxides, followed by assessing the role of the ligands in the reactivity of the catalysts. For better selectivity, we propose that metal oxide anionic complexes should be explored further, while hydrophylic cavities in the vicinity of the metal oxide can perturb the transition-state structure for methanol increasing appreciably the activation barrier for methanol. We also emphasize that computational studies should target the activation reaction of methanol (and not only methane), the study of complete catalytic cycles (including the recombination and oxidation steps), and the use of molecular oxygen as an oxidant. The titled chemical conversion is an excellent challenge for theory and we believe that computational studies should lead the field in the future. It is finally shown that bottom-up approaches offer a systematic way for exploration of the chemical space and should still be applied in parallel with the recently popular machine learning techniques. To answer the question of the title, we believe that metal oxides should still be considered provided that we change our focus and perform more systematic investigations on the activation of methanol.

Quantifiable polarity match effect on C-H bond cleavage reactivity and its limits in reaction design

When oxidants favour cleaving a strong C-H bond at the expense of weaker ones, which are otherwise inherently preferred due to their favourable reaction energy, reactivity factors such as the polarity match effect are often invoked. Polarity match follows the intuition of electrophilic (nucleophilic) oxidants reacting faster with nucleophilic (electrophilic) C-H bonds. Nevertheless, this concept is purely qualitative and is best suited for a posteriori rationalization of experimental observations. Here, we propose and inspect two methods to quantify polar effects in C-H cleavage reactions, one by computation via the difference of atomic charges (Δq) of reacting atoms, and one amenable to experimental measurement through asynchronicity factors, η. By their application to three case studies, we observe that both Δq and η faithfully capture the notion of polarity match. The polarity match model, however, proves insufficient as a predictor of H-atom abstraction reactivity and we discourage its use as a standalone variable in reaction design. Besides this caveat, η and Δq (through its mapping on η) allow the implementation of polarity match into a Marcus-type model of reactivity, alleviating its shortcomings and making reaction planning feasible.