Carbocation−π Interaction: Computational Study of Complexation of Methyl Cation with Benzene and Comparisons with Related Systems

Switchable Decarboxylation by Energy- or Electron-Transfer Photocatalysis

Kolbe dimerization and Hofer-Moest reactions are well-investigated carboxylic acid transformations, wherein new carbon-carbon and carbon-heteroatom bonds are constructed via electrochemical decarboxylation. These transformations can be switched by choosing an electrode that allows control of the reactive intermediate, such as carbon radical or carbocation. However, the requirement of a high current density diminishes the functional group compatibility with these electrochemical reactions. Here, we demonstrate the photocatalytic decarboxylative transformation of activated carboxylic acids in a switchable and functional group-compatible manner. We discovered that switching between Kolbe-type or Hofer-Moest-type reactions can be accomplished with suitable photocatalysts by controlling the reaction pathways: energy transfer (EnT) and single-electron transfer (SET). The EnT pathway promoted by an organo-photocatalyst yielded 1,2-diarylethane from arylacetic acids, whereas the ruthenium photoredox catalyst allows the construction of an ester scaffold with two arylmethyl moieties via the SET pathway. The resulting radical intermediates were coupled to olefins to realize multicomponent reactions. Consequently, four different products were selectively obtained from a simple carboxylic acid. This discovery offers new opportunities for selectively synthesizing multiple products via switchable reactions using identical substrates with minimal cost and effort.

Tryptophan Stabilization of a Biochemical Carbocation Evaluated by Analysis of π Complexes of 3-Ethylindole with the t-Butyl Cation

Understanding how the highly unstable carbocation intermediates in terpenoid biosynthesis are stabilized and protected during their transient existence in enzyme active sites is an intriguing challenge which has to be addressed computationally. Our efforts have focused on evaluating the stabilization afforded via carbocation-π complexation between a biochemical carbocation and an aromatic amino acid residue. This has involved making measurements on an X-ray structure of an enzyme active site that shows a π donor proximate to a putative carbocation site and using these to build models which are analyzed computationally to provide an estimated stabilization energy (SE). Previously, we reported estimated SEs for several such carbocation-π complexes involving phenylalanine. Herein, we report the first such estimate involving tryptophan as the π donor. Because there was almost no published information about indole as a π-complexation donor, we first located computationally equilibrium π and σ complexes of 3-ethylindole with the t-butyl cation as relevant background information. Then, measurements on the X-ray structure of the enzyme CotB2 complexed with geranylgeranyl thiodiphosphate (GGSPP), specifically on the geometric relationship of the putative carbocation at C15 of GGSPP to W186, were used to build a model that afforded a computed SE of -15.3 kcal/mol.

Carbocation chemistry confined in zeolites: spectroscopic and theoretical characterizations

Acid-catalyzed reactions inside zeolites are one type of broadly applied industrial reactions, where carbocations are the most common intermediates of these reaction processes, including methanol to olefins, alkene/aromatic alkylation, and hydrocarbon cracking/isomerization. The fundamental research on these acid-catalyzed reactions is focused on the stability, evolution, and lifetime of carbocations under the zeolite confinement effect, which greatly affects the efficiency, selectivity and deactivation of zeolite catalysts. Therefore, a profound understanding of the carbocations confined in zeolites is not only beneficial to explain the reaction mechanism but also drive the design of new zeolite catalysts with ideal acidity and cages/channels. In this review, we provide both an in-depth understanding of the stabilization of carbocations by the pore confinement effect and summary of the advanced characterization methods to capture carbocations in zeolites, including UV-vis spectroscopy, solid-state NMR, fluorescence microscopy, IR spectroscopy and Raman spectroscopy. Also, we clarify the relationship between the activity and stability of carbocations in zeolite-catalyzed reactions, and further highlight the role of carbocations in various hydrocarbon conversion reactions inside zeolites with diverse frameworks and varying acidic properties.

A Benchmark Study of Quantum Mechanics and Quantum Mechanics-Molecular Mechanics Methods for Carbocation Chemistry

Carbocations play key roles in classical organic reactions and have also been implicated in several enzyme families. A hallmark of carbocation chemistry is multitudes of competing reaction pathways, and to be able to distinguish between pathways with quantum chemical calculations, it is necessary to approach chemical accuracy for relative energies between carbocations. Here, we present an extensive study of the performance of selected density functional theory (DFT) methods in describing the thermochemistry and kinetics of carbocations and their corresponding neutral alkenes both in the gas-phase and within a hybrid quantum mechanics-molecular mechanics (QM/MM) framework. The density functionals are benchmarked against accurate ab initio methods such as CBS-QB3 and DLPNO-CCSD(T). Based on the findings in the gas-phase calculations of carbocations and alkenes, the best functionals are chosen and tested further for non-covalent interactions in model systems using QM and QM/MM methods. We compute the interaction energies between a model carbocation/alkane and model π, dipole, and hydrophobic systems using DFT and QM(DFT)/MM and compare with DLPNO-CCSD(T). These latter model systems are representative of side chains of amino acids such as phenylalanine/tyrosine, tryptophan, asparagine/glutamine, serine/threonine, methionine, and other hydrophobic groups. The Lennard-Jones parameters of the QM atoms in QM(DFT)/MM calculations are modified to obtain an optimal fit with the QM energies. Finally, a selected carbocation reaction is studied in the gas phase and in implicit chloroform solvent using QM and in explicit chloroform solvent using QM/MM and umbrella sampling simulations. This study highlights the highest accuracy possible with selected density functionals and QM/MM methods but also some limitations in using QM/MM methods for carbocation systems.

A simpler method affords evaluation of π stabilization by phenylalanine of several biochemical carbocations

Carbocations are important intermediates in the biosynthesis of terpenes and steroids, and it is challenging to try to understand how these relatively unstable species survive even transiently during biochemical reactions. Carbocation-π interaction with aromatic amino acid residues is an important factor in helping to stabilize these positively charged species. However, the short lifetimes of these active site carbocations makes experimental evaluation of the stabilization afforded by such interaction impossible. Computational studies, however, have provided some insight into this phenomenon. Herein we report a simple, computationally efficient method to estimate such stabilization energies afforded by phenylalanine to biochemical carbocation intermediates. A model is constructed in which the biochemical carbocation is replaced by an appropriate carbocation mimic (t-butyl or dimethylallyl). This substitute carbocation is then aligned with an ethylbenzene serving as a surrogate for each proximate phenylalanine in a geometry that replicates as closely as possible the orientation of that phenylalanine using measurements made on an X-ray structure of an enzyme active site in which a carbocation surrogate is bound. Density functional theory computations on such models were then used to yield estimates of stabilization energies. Application of this method to the tertiary carbocation formed in the reaction catalyzed by geranyl diphosphate C-methyl transferase gave a stabilization energy (-12.3 kcal mol-1) that was essentially identical to that obtained previously by analysis of a much more computationally demanding model of the active site. As a check on the accuracy of the simpler method, it was applied with similar success to the farnesyl cation formed in the reaction catalyzed by aristolochene synthase that is stabilized by cation-π interaction with two phenylalanines. Application of this method is also described to estimate carbocation-π stabilization, by the same two phenylalanines, of the final carbocation intermediate leading to aristolochene through analysis of the X-ray structure of an inhibitor of that carbocation bound in the active site of aristolochene synthase. Finally, the stabilization, by either of two phenylalanines, of six carbocation intermediates in the oxidosqualene cyclase-catalyzed formation of lanosterol is estimated by comparable analysis of an X-ray structure of that reaction product bound in the enzyme active site.

A Protonated Quinone Methide Stabilized by a Combination of Partial Aromatization and π-Interaction: Spectroscopic and Crystallographic Analysis

We have expanded the repertoire of cation-π interactions to include a carbocation-π system resulting from the protonation of a π-stacked para-quinone methide (p-QM). This unusual carbocation is stabilized by a combination of partial aromatization of the QM moiety and through-space interaction with the π-system of the adjacent aromatic ring. Single crystal X-ray analysis of the protonated form reveals a structure consisting of a hydrogen-bound complex involving two molecules of the precursor and one proton.

Navigating Past a Fork in the Road: Carbocation-π Interactions Can Manipulate Dynamic Behavior of Reactions Facing Post-Transition-State Bifurcations

Dynamics calculations are described for carbocation rearrangements involving product-forming pathways with post-transition-state bifurcations. We show that noncovalent interactions with associated benzene rings (a simple model of aromatic amino acid side chains) can switch inherent dynamical tendencies for competing modes of disrotation, establishing that meaningful changes in dynamically controlled product selectivity can be achieved with few weak noncovalent interactions.

Carbocation–π interaction: evaluation of the stabilization by phenylalanine of a biochemical carbocation intermediate

Computational analyses, using primarily density functional theory, have been used to determine the stabilization associated with the carbocation–π interaction of a biochemical carbocation intermediate binding to a phenylalanine residue in an enzyme active site.

Gas-phase electrophilic aromatic substitution mechanism with strong electrophiles explained by ab initio non-adiabatic dynamics

Ab initio non-adiabatic dynamics is used to monitor the attack of CH3(+) to benzene. The results show that in the gas phase the reaction is ultrafast and is governed by a single electron transfer producing a radical pair.

Inherent dynamical preferences in carbocation rearrangements leading to terpene natural products

An introduction to the application of quantum chemical dynamics calculations to mechanistic problems in the field of terpene biosynthesis is provided. A bare bones introduction to the fundamentals of chemical dynamics is followed by a brief account of previous applications to terpene-forming carbocation reactions, a discussion of questions in this field that dynamics calculations may help answer, and a description of current problems to which dynamics calculations are being applied.

Computational study of isopropylbenzenium ions

A theoretical investigation of the isopropylbenzenium ion system has been carried out with structures determined with B3LYP. Energies are calculated with the high accuracy composite methods, G3 (G3B3) and CBS (CBS-QB3). The main goal has been to resolve the following issue: Are there any stable ion π-electron complexes of the type C(6)H(6)/C(3)H(7)(+) or C(6)H(7)(+)/C(3)H(6) in this system? Two minimum points on the potential energy surface (PES) corresponding to benzenium ion/propene complexes were found. Due to free internal rotation, they represent only one species. The barrier for forming an isopropylbenzenium ion from the complex is low, so the lifetime will be short. Computation of the IR spectrum of the complex shows that there is a very intense absorption line due to C-H stretching, in an otherwise empty region, that may be used to identify the complex, if present. No stable C(6)H(6)/C(3)H(7)(+) complex was found, but a quasi-stable species structurally corresponding to the earlier described stable C(6)H(6)/C(4)H(9)(+) complex was observed. A simplistic explanation why the benzene/isopropyl cation, in contrast to benzene/tert-butyl cation, does not form a stable ion/π-electron complex is given.

The taxadiene-forming carbocation cascade

A complete pathway (structures and energies of intermediates and transition state structures connecting them) from geranylgeranyl diphosphate to taxadiene, obtained using quantum chemical calculations, is described. This pathway is fully consistent with previous labeling experiments, despite differing in several subtle ways (in terms of conformations of certain carbocation intermediates and in the concertedness and synchronicity of certain bond-forming events) from previous mechanistic proposals. Also, on the basis of the theoretical results, it is proposed that the 2-fluoro-geranylgeranyl diphosphate substrate analogue in the recently reported X-ray crystal structure of taxadiene synthase is bound in a nonproductive orientation.

Alkylation of Phenol: A Mechanistic View

The current work utilizes the ab initio density functional theory (DFT) to develop a molecular level of the mechanistic understanding on the phenol alkylation in the presence of a cation-exchange resin catalyst, Amberlyst-15. The catalyst is modeled with the benzene sulfonic acid, and the effect of this acid on olefins such as isopropene (i-Pr) and tributene (t-Bu) in a phenol solution mimics the experimental condition. A neutral-pathway mechanism is established to account for early-stage high concentration of the phenolic ether observed in experiments. The mechanism involves an exothermic reaction between olefin and the benzene sulfonic acid to form ester followed by three reaction pathways leading to direct O-alkylation, o-C-alkylation, and p-C-alkylation. Our calculations conclude that O-alkylation to form the phenolic ether is the most energetically favorable in the neutral condition. An ionic rearrangement mechanism describes intramolecular migrations of the alkyl group from the phenolic ether to form C-alkylphenols, while the positively charged protonation significantly lowers transition barriers for these migrations. The ionic rearrangement mechanism accounts for high yields of o-C-alkylphenol and p-C-alkylphenol. Competition between the H atom and the alkyl R group at the substitutive site of the protonated ortho configuration is found to be the determining factor to the ortho/para ratio of C-alkylation products.

Cation−π Interactions with a Model for the Side Chain of Tryptophan: Structures and Absolute Binding Energies of Alkali Metal Cation−Indole Complexes†

Threshold collision-induced dissociation techniques are employed to determine bond dissociation energies (BDEs) of mono- and bis-complexes of alkali metal cations, Li+, Na+, K+, Rb+, and Cs+, with indole, C8H7N. The primary and lowest energy dissociation pathway in all cases is endothermic loss of an intact indole ligand. Sequential loss of a second indole ligand is observed at elevated energies for the bis-complexes. Density functional theory calculations at the B3LYP/6-31G level of theory are used to determine the structures, vibrational frequencies, and rotational constants of these complexes. Theoretical BDEs are determined from single point energy calculations at the MP2(full)/6-311+G(2d,2p) level using the B3LYP/6-31G* geometries. The agreement between theory and experiment is very good for all complexes except Li+ (C8H7N), where theory underestimates the strength of the binding. The trends in the BDEs of these alkali metal cation-indole complexes are compared with the analogous benzene and naphthalene complexes to examine the influence of the extended pi network and heteroatom on the strength of the cation-pi interaction. The Na+ and K+ binding affinities of benzene, phenol, and indole are also compared to those of the aromatic amino acids, phenylalanine, tyrosine, and tryptophan to elucidate the factors that contribute to the binding in complexes to the aromatic amino acids. The nature of the binding and trends in the BDEs of cation-pi complexes between alkali metal cations and benzene, phenol, and indole are examined to help understand nature's preference for engaging tryptophan over phenylalanine and tyrosine in cation-pi interactions in biological systems.

Insight into steroid scaffold formation from the structure of human oxidosqualene cyclase

In higher organisms the formation of the steroid scaffold is catalysed exclusively by the membrane-bound oxidosqualene cyclase (OSC; lanosterol synthase). In a highly selective cyclization reaction OSC forms lanosterol with seven chiral centres starting from the linear substrate 2,3-oxidosqualene. Valuable data on the mechanism of the complex cyclization cascade have been collected during the past 50 years using suicide inhibitors, mutagenesis studies and homology modelling. Nevertheless it is still not fully understood how the enzyme catalyses the reaction. Because of the decisive role of OSC in cholesterol biosynthesis it represents a target for the discovery of novel anticholesteraemic drugs that could complement the widely used statins. Here we present two crystal structures of the human membrane protein OSC: the target protein with an inhibitor that showed cholesterol lowering in vivo opens the way for the structure-based design of new OSC inhibitors. The complex with the reaction product lanosterol gives a clear picture of the way in which the enzyme achieves product specificity in this highly exothermic cyclization reaction.