skip to main content

Title: Chemical Bonding in Homogenous Catalysis – Seen Through the Eyes of Vibrational Spectroscopy
We introduce two new tools for the analysis of bond forming/breaking processes taking place during catalytic reactions, the Uni!ed Reaction Valley Approach (URVA) and the Local Mode Analysis (LMA), both being based on vibrational spectroscopy. We discuss how URVA and LMA complement currently used computational approaches and provide valuable insights into catalytic processes, supporting current design efforts aiming at more ef!cient and environmentally friendly catalysts. Three examples are presented; Au- catalyzed [3,3]-sigmatropic rearrangement of allyl acetate, Re-catalyzed CO2 cycloaddition to epoxides, and a-ketoamide inhibitors for SARS-CoV-2 main protease. We hope that URVA and LMA will become routinely applied tools in computational catalysis and also enter the classroom.
Authors:
Award ID(s):
2102461
Publication Date:
NSF-PAR ID:
10342531
Journal Name:
eference Module in Chemistry, Molecular Sciences and Chemical Engineering - Comprehensive Computational Chemistry
Page Range or eLocation-ID:
1-22
Sponsoring Org:
National Science Foundation
More Like this
  1. In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationarymore »points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.« less
  2. In this study we investigate the Diels–Alder reaction between methyl acrylate and butadiene, which is catalyzed by BF3 Lewis acid in explicit water solution, using URVA and Local Mode Analysis as major tools complemented with NBO, electron density and ring puckering analyses. We considered four different starting orientations of methyl acrylate and butadiene, which led to 16 DA reactions in total. In order to isolate the catalytic effects of the BF3 catalyst and those of the water environment and exploring how these effects are synchronized, we systematically compared the non-catalyzed reaction in gas phase and aqueous solution with the catalyzed reaction in gas phase and aqueous solution. Gas phase studies were performed at the B3LYP/6-311+G(2d,p) level of theory and studies in aqueous solution were performed utilizing a QM/MM approach at the B3LYP/6-311+G(2d,p)/AMBER level of theory. The URVA results revealed reaction path curvature profiles with an overall similar pattern for all 16 reactions showing the same sequence of CC single bond formation for all of them. In contrast to the parent DA reaction with symmetric substrates causing a synchronous bond formation process, here, first the new CC single bond on the CH2 side of methyl acrylate is formed followed by themore »CC bond at the ester side. As for the parent DA reaction, both bond formation events occur after the TS, i.e., they do not contribute to the energy barrier. What determines the barrier is the preparation process for CC bond formation, including the approach diene and dienophile, CC bond length changes and, in particular, rehybridization of the carbon atoms involved in the formation of the cyclohexene ring. This process is modified by both the BF3 catalyst and the water environment, where both work in a hand-in-hand fashion leading to the lowest energy barrier of 9.06 kcal/mol found for the catalyzed reaction R1 in aqueous solution compared to the highest energy barrier of 20.68 kcal/mol found for the non-catalyzed reaction R1 in the gas phase. The major effect of the BF3 catalyst is the increased mutual polarization and the increased charge transfer between methyl acrylate and butadiene, facilitating the approach of diene and dienophile and the pyramidalization of the CC atoms involved in the ring formation, which leads to a lowering of the activation energy. The catalytic effect of water solution is threefold. The polar environment leads also to increased polarization and charge transfer between the reacting species, similar as in the case of the BF3 catalyst, although to a smaller extend. More important is the formation of hydrogen bonds with the reaction complex, which are stronger for the TS than for the reactant, thus stabilizing the TS which leads to a further reduction of the activation energy. As shown by the ring puckering analysis, the third effect of water is space confinement of the reacting partners, conserving the boat form of the six-member ring from the entrance to the exit reaction channel. In summary, URVA combined with LMA has led to a clearer picture on how both BF3 catalyst and aqueous environment in a synchronized effort lower the reaction barrier. These new insights will serve to further fine-tune the DA reaction of methyl acrylate and butadiene and DA reactions in general.« less
  3. Abstract

    Electrocatalytic nanocarbon (EN) is a class of material receiving intense interest as a potential replacement for expensive, metal-based electrocatalysts for energy conversion and chemical production applications. The further development of EN will require an intricate knowledge of its catalytic behaviors, however, the true nature of their electrocatalytic activity remains elusive. This review highlights work that contributed valuable knowledge in the elucidation of EN catalytic mechanisms. Experimental evidence from spectroscopic studies and well-defined molecular models, along with the survey of computational studies, is summarized to document our current mechanistic understanding of EN-catalyzed oxygen, carbon dioxide and nitrogen electrochemistry. We hope this review will inspire future development of synthetic methods and in situ spectroscopic tools to make and study well-defined EN structures.

  4. Bacterial tyrosine kinases (BY-kinases) comprise a family of protein tyrosine kinases that are structurally distinct from their functional counterparts in eukaryotes and are highly conserved across the bacterial kingdom. BY-kinases act in concert with their counteracting phosphatases to regulate a variety of cellular processes, most notably the synthesis and export of polysaccharides involved in biofilm and capsule biogenesis. Biochemical data suggest that BY-kinase function involves the cyclic assembly and disassembly of oligomeric states coupled to the overall phosphorylation levels of a C-terminal tyrosine cluster. This process is driven by the opposing effects of intermolecular autophosphorylation, and dephosphorylation catalyzed by tyrosine phosphatases. In the absence of structural insight into the interactions between a BY-kinase and its phosphatase partner in atomic detail, the precise mechanism of this regulatory process has remained poorly defined. To address this gap in knowledge, we have determined the structure of the transiently assembled complex between the catalytic core of the Escherichia coli (K-12) BY-kinase Wzc and its counteracting low–molecular weight protein tyrosine phosphatase (LMW-PTP) Wzb using solution NMR techniques. Unambiguous distance restraints from paramagnetic relaxation effects were supplemented with ambiguous interaction restraints from static spectral perturbations and transient chemical shift changes inferred from relaxation dispersion measurements andmore »used in a computational docking protocol for structure determination. This structure presents an atomic picture of the mode of interaction between an LMW-PTP and its BY-kinase substrate, and provides mechanistic insight into the phosphorylation-coupled assembly/disassembly process proposed to drive BY-kinase function.« less
  5. Gram-negative bacteria expressing class A β-lactamases pose a serious health threat due to their ability to inactivate all β-lactam antibiotics. The acyl–enzyme intermediate is a central milestone in the hydrolysis reaction catalyzed by these enzymes. However, the protonation states of the catalytic residues in this complex have never been fully analyzed experimentally due to inherent difficulties. To help unravel the ambiguity surrounding class A β-lactamase catalysis, we have used ultrahigh-resolution X-ray crystallography and the recently approved β-lactamase inhibitor avibactam to trap the acyl–enzyme complex of class A β-lactamase CTX-M-14 at varying pHs. A 0.83-Å-resolution CTX-M-14 complex structure at pH 7.9 revealed a neutral state for both Lys73 and Glu166. Furthermore, the avibactam hydroxylamine-O-sulfonate group conformation varied according to pH, and this conformational switch appeared to correspond to a change in the Lys73 protonation state at low pH. In conjunction with computational analyses, our structures suggest that Lys73 has a perturbed acid dissociation constant (pKa) compared with acyl–enzyme complexes with β-lactams, hindering its function to deprotonate Glu166 and the initiation of the deacylation reaction. Further NMR analysis demonstrated Lys73 pKato be ∼5.2 to 5.6. Together with previous ultrahigh-resolution crystal structures, these findings enable us to follow the proton transfer process ofmore »the entire acylation reaction and reveal the critical role of Lys73. They also shed light on the stability and reversibility of the avibactam carbamoyl acyl–enzyme complex, highlighting the effect of substrate functional groups in influencing the protonation states of catalytic residues and subsequently the progression of the reaction.

    « less