How Mg2+ ions lower the SN2@P barrier in enzymatic triphosphate hydrolysis

Architecture, Function, Regulation, and Evolution of α-Glucans Metabolic Enzymes in Prokaryotes

Bacteria have acquired sophisticated mechanisms for assembling and disassembling polysaccharides of different chemistry. α-d-Glucose homopolysaccharides, so-called α-glucans, are the most widespread polymers in nature being key components of microorganisms. Glycogen functions as an intracellular energy storage while some bacteria also produce extracellular assorted α-glucans. The classical bacterial glycogen metabolic pathway comprises the action of ADP-glucose pyrophosphorylase and glycogen synthase, whereas extracellular α-glucans are mostly related to peripheral enzymes dependent on sucrose. An alternative pathway of glycogen biosynthesis, operating via a maltose 1-phosphate polymerizing enzyme, displays an essential wiring with the trehalose metabolism to interconvert disaccharides into polysaccharides. Furthermore, some bacteria show a connection of intracellular glycogen metabolism with the genesis of extracellular capsular α-glucans, revealing a relationship between the storage and structural function of these compounds. Altogether, the current picture shows that bacteria have evolved an intricate α-glucan metabolism that ultimately relies on the evolution of a specific enzymatic machinery. The structural landscape of these enzymes exposes a limited number of core catalytic folds handling many different chemical reactions. In this Review, we present a rationale to explain how the chemical diversity of α-glucans emerged from these systems, highlighting the underlying structural evolution of the enzymes driving α-glucan bacterial metabolism.

Computational insights into the inverse electron-demand Diels-Alder reaction of norbornenes with 1,2,4,5-tetrazines: norbornene substituents' effects on the reaction rate

The study of the reaction rates and mechanism of click chemistry reactions still remains an interesting challenge in organic chemistry. In this regard, the inverse electron demand Diels-Alder (IEDDA) reaction represents a promising metal-free alternative with enhanced reaction rates compared to other reactions of the click chemistry toolbox. Among the different types of dienophiles used in the IEDDA reactions, norbornenes have been widely used given their high stability and fast reaction rates. The inverse electron-demand Diels Alder reaction of 3,6-dipyridin-2-yl-1,2,4,5-tetrazine with a series of norbornene derivatives was studied with quantum mechanical calculations at the M06-2X/6-311+G(d,p) level of theory. The theoretical predictions were confirmed with the experimental data and analyzed with the use of the distortion/interaction model. The obtained results will help in obtaining a better understanding of the factors that affect the relative cycloaddition rates of norbornenes with tetrazines, which are crucial for selectively tuning their efficacy.

Chemical reactivity from an activation strain perspective

Chemical reactions are ubiquitous in the universe, they are at the core of life, and they are essential for industrial processes. The drive for a deep understanding of how something occurs, in this case, the mechanism of a chemical reaction and the factors controlling its reactivity, is intrinsically valuable and an innate quality of humans. The level of insight and degree of understanding afforded by computational chemistry cannot be understated. The activation strain model is one of the most powerful tools in our arsenal to obtain unparalleled insight into reactivity. The relative energy of interacting reactants is evaluated along a reaction energy profile and related to the rigidity of the reactants' molecular structure and the strength of the stabilizing interactions between the deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). Owing to the connectedness between the activation strain model and Kohn-Sham molecular orbital theory, one is able to obtain a causal relationship between both the sterics and electronics of the reactants and their mutual reactivity. Only when this is accomplished one can eclipse the phenomenological explanations that are commonplace in the literature and textbooks and begin to rationally tune and optimize chemical transformations. We showcase how the activation strain model is the ideal tool to elucidate fundamental organic reactions, the activation of small molecules by metallylenes, and the cycloaddition reactivity of cyclic diene- and dipolarophiles.

How Lewis Acids Catalyze Ring-Openings of Cyclohexene Oxide

We have quantum chemically studied the Lewis acid-catalyzed epoxide ring-opening reaction of cyclohexene epoxide by MeZH (Z = O, S, and NH) using relativistic dispersion-corrected density functional theory. We found that the reaction barrier of the Lewis acid-catalyzed epoxide ring-opening reactions decreases upon ascending in group 1 along the series Cs+ > Rb+ > K+ > Na+ > Li+ > H+. Our activation strain and Kohn-Sham molecular orbital analyses reveal that the enhanced reactivity of the Lewis acid-catalyzed ring-opening reaction is caused by the reduced steric (Pauli) repulsion between the filled orbitals of the epoxide and the nucleophile, as the Lewis acid polarizes the filled orbitals of the epoxide more efficiently away from the incoming nucleophile. Furthermore, we established that the regioselectivity of these ring-opening reactions is, aside from the "classical" strain control, also dictated by a hitherto unknown mechanism, namely, the steric (Pauli) repulsion between the nucleophile and the substrate, which could be traced back to the asymmetric orbital density on the epoxide. In all, this work again demonstrates that the concept of Pauli-lowering catalysis is a general phenomenon.

Mechanism of Deoxyguanosine Diphosphate Insertion by Human DNA Polymerase β

DNA polymerases play vital roles in the maintenance and replication of genomic DNA by synthesizing new nucleotide polymers using nucleoside triphosphates as substrates. Deoxynucleoside triphosphates (dNTPs) are the canonical substrates for DNA polymerases; however, some bacterial polymerases have been demonstrated to insert deoxynucleoside diphosphates (dNDPs), which lack a third phosphate group, the γ-phosphate. Whether eukaryotic polymerases can efficiently incorporate dNDPs has not been investigated, and much about the chemical or structural role played by the γ-phosphate of dNTPs remains unknown. Using the model mammalian polymerase (Pol) β, we examine how Pol β incorporates a substrate lacking a γ-phosphate [deoxyguanosine diphosphate (dGDP)] utilizing kinetic and crystallographic approaches. Using single-turnover kinetics, we determined dGDP insertion across a templating dC by Pol β to be drastically impaired when compared to dGTP insertion. We found the most significant impairment in the apparent insertion rate (k_pol), which was reduced 32000-fold compared to that of dGTP insertion. X-ray crystal structures revealed similar enzyme-substrate contacts for both dGDP and dGTP. These findings suggest the insertion efficiency of dGDP is greatly decreased due to impairments in polymerase chemistry. This work is the first instance of a mammalian polymerase inserting a diphosphate nucleotide and provides insight into the nature of polymerase mechanisms by highlighting how these enzymes have evolved to use triphosphate nucleotide substrates.

A Unified Framework for Understanding Nucleophilicity and Protophilicity in the SN 2/E2 Competition

The concepts of nucleophilicity and protophilicity are fundamental and ubiquitous in chemistry. A case in point is bimolecular nucleophilic substitution (SN 2) and base-induced elimination (E2). A Lewis base acting as a strong nucleophile is needed for SN 2 reactions, whereas a Lewis base acting as a strong protophile (i.e., base) is required for E2 reactions. A complicating factor is, however, the fact that a good nucleophile is often a strong protophile. Nevertheless, a sound, physical model that explains, in a transparent manner, when an electron-rich Lewis base acts as a protophile or a nucleophile, which is not just phenomenological, is currently lacking in the literature. To address this fundamental question, the potential energy surfaces of the SN 2 and E2 reactions of X- +C2 H5 Y model systems with X, Y = F, Cl, Br, I, and At, are explored by using relativistic density functional theory at ZORA-OLYP/TZ2P. These explorations have yielded a consistent overview of reactivity trends over a wide range in reactivity and pathways. Activation strain analyses of these reactions reveal the factors that determine the shape of the potential energy surfaces and hence govern the propensity of the Lewis base to act as a nucleophile or protophile. The concepts of "characteristic distortivity" and "transition state acidity" of a reaction are introduced, which have the potential to enable chemists to better understand and design reactions for synthesis.

Ambident Nucleophilic Substitution: Understanding Non-HSAB Behavior through Activation Strain and Conceptual DFT Analyses

The ability to understand and predict ambident reactivity is key to the rational design of organic syntheses. An approach to understand trends in ambident reactivity is the hard and soft acids and bases (HSAB) principle. The recent controversy over the general validity of this principle prompted us to investigate the competing gas-phase SN 2 reaction channels of archetypal ambident nucleophiles CN- , OCN- , and SCN- with CH3 Cl (SN 2@C) and SiH3 Cl (SN 2@Si), using DFT calculations. Our combined analyses highlight the inability of the HSAB principle to correctly predict the reactivity trends of these simple, model reactions. Instead, we have successfully traced reactivity trends to the canonical orbital-interaction mechanism and the resulting nucleophile-substrate interaction energy. The HOMO-LUMO orbital interactions set the trend in both SN 2@C and SN 2@Si reactions. We provide simple rules for predicting the ambident reactivity of nucleophiles based on our Kohn-Sham molecular orbital analysis.

PyFrag 2019-Automating the exploration and analysis of reaction mechanisms

We present a substantial update to the PyFrag 2008 program, which was originally designed to perform a fragment-based activation strain analysis along a provided potential energy surface. The original PyFrag 2008 workflow facilitated the characterization of reaction mechanisms in terms of the intrinsic properties, such as strain and interaction, of the reactants. The new PyFrag 2019 program has automated and reduced the time-consuming and laborious task of setting up, running, analyzing, and visualizing computational data from reaction mechanism studies to a single job. PyFrag 2019 resolves three main challenges associated with the automated computational exploration of reaction mechanisms: it (1) computes the reaction path by carrying out multiple parallel calculations using initial coordinates provided by the user; (2) monitors the entire workflow process; and (3) tabulates and visualizes the final data in a clear way. The activation strain and canonical energy decomposition results that are generated relate the characteristics of the reaction profile in terms of intrinsic properties (strain, interaction, orbital overlaps, orbital energies, populations) of the reactant species. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Dual Activation of Aromatic Diels-Alder Reactions

The unusually fast Diels–Alder reactions of [5]cyclophanes were analyzed by DFT at the BLYP‐D3(BJ)/TZ2P level of theory. The computations were guided by an integrated activation‐strain and Kohn–Sham molecular orbital analysis. It is revealed why both [5]metacyclophane and [5]paracyclophane exhibit a significant rate enhancement compared to their planar benzene analogue. The activation strain analyses revealed that the enhanced reactivity originates from 1) predistortion of the aromatic core resulting in a reduced activation strain of the aromatic diene, and/or 2) enhanced interaction with the dienophile through a distortion‐controlled lowering of the HOMO–LUMO gap within the diene. Both of these physical mechanisms and thus the rate of Diels–Alder cycloaddition can be tuned through different modes of geometrical distortion (meta versus para bridging) and by heteroatom substitution in the aromatic ring. Judicious choice of the bridge and heteroatom in the aromatic core enables effective tuning of the aromatic Diels–Alder reactivity to achieve activation barriers as low as 2 kcal mol⁻¹, which is an impressive 35 kcal mol⁻¹ lower than that of benzene.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Nucleophilic Substitution (SN 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent

The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.

Nucleophilic Substitution in Solution: Activation Strain Analysis of Weak and Strong Solvent Effects

We have quantum chemically studied the effect of various polar and apolar solvents on the shape of the potential energy surface (PES) of a diverse collection of archetypal nucleophilic substitution reactions at carbon, silicon, phosphorus, and arsenic by using density functional theory at the OLYP/TZ2P level. In the gas phase, all our model SN 2 reactions have single-well PESs, except for the nucleophilic substitution reaction at carbon (SN 2@C), which has a double-well energy profile. The presence of the solvent can have a significant effect on the shape of the PES and, thus, on the nature of the SN 2 process. Solvation energies, charges on the nucleophile or leaving group, and structural features are compared for the various SN 2 reactions in a spectrum of solvents. We demonstrate how solvation can change the shape of the PES, depending not only on the polarity of the solvent, but also on how the charge is distributed over the interacting molecular moieties during different stages of the reaction. In the case of a nucleophilic substitution at three-coordinate phosphorus, the reaction can be made to proceed through a single-well [no transition state (TS)], bimodal barrier (two TSs), and then through a unimodal transition state (one TS) simply by increasing the polarity of the solvent.