Search

De novo enzyme design has sought to introduce active sites and substrate-binding pockets that are predicted to catalyse a reaction of interest into geometrically compatible native scaffolds1,2, but has been limited by a lack of suitable protein structures and the complexity of native protein sequence–structure relationships. Here we describe a deep-learning-based ‘family-wide hallucination’ approach that generates large numbers of idealized protein structures containing diverse pocket shapes and designed sequences that encode them. We use these scaffolds to design artificial luciferases that selectively catalyse the oxidative chemiluminescence of the synthetic luciferin substrates diphenylterazine3 and 2-deoxycoelenterazine. The designed active sites position an arginine guanidinium group adjacent to an anion that develops during the reaction in a binding pocket with high shape complementarity. For both luciferin substrates, we obtain designed luciferases with high selectivity; the most active of these is a small (13.9 kDa) and thermostable (with a melting temperature higher than 95 °C) enzyme that has a catalytic efficiency on diphenylterazine (kcat/Km = 106 M−1 s−1) comparable to that of native luciferases, but a much higher substrate specificity. The creation of highly active and specific biocatalysts from scratch with broad applications in biomedicine is a key milestone for computational enzyme design, and our approach should enable generation of a wide range of luciferases and other enzymes.

The Diels-Alder cycloaddition is one of the most powerful approaches in organic synthesis and is often used in the synthesis of important pharmaceuticals. Yet, strictly controlling the stereoselectivity of the Diels-Alder reactions is challenging, and great efforts are needed to construct complex molecules with desired chirality via organocatalysis or transition metal strategies. Nature has evolved different types of enzymes to exquisitely control cyclization stereochemistry, however, most of the reported Diels-Alderases have been shown to only facilitate the energetically favorable diastereoselective cycloadditions. Here we report the discovery and characterization of CtdP, as a member of a new class of bifunctional oxidoreductase/Diels-Alderase, which was previously annotated as an NmrA-like transcriptional regulator. We demonstrate that CtdP catalyses the inherently unfavored cycloaddition to form the bicyclo[2,2,2]diazaoctane scaffold with a strict α-anti-selectivity. Guided by computational studies, we reveal a NADP(+)/NADPH-dependent redox mechanism for the CtdP-catalysed inverse electron demand Diels-Alder cycloaddition, which serves as the first example of a bifunctional Diels-Alderase that utilizes this mechanism.

An efficient alkoxyl radical-triggered ring expansion/cross-coupling cascade was developed under cheap metal catalysis. Through the metal-catalyzed radical relay strategy, a wide range of medium-sized lactones (9-11 membered) and macrolactones (12, 13, 15, 18, and 19-membered) were constructed in moderate to good yields, along with diverse functional groups including CN, N3, SCN, and X groups installed concurrently. Density functional theory (DFT) calculations revealed that reductive elimination of the cycloalkyl-Cu(iii) species is a more favorable reaction pathway for the cross-coupling step. Based on the results of experiments and DFT, a Cu(i)/Cu(ii)/Cu(iii) catalytic cycle is proposed for this tandem reaction.

The importance of fluoro and trifluoromethyl substituents in drug effectiveness prompted the computational exploration of fluorine-containing substituents in valuable synthetic cycloadditions. Diels-Alder or 1,3-dipolar cycloaddition reactions of typical reactants, cyclopentadiene

Keteniminium ions, the nitrogen analogues of ketenes, exhibit high reactivity toward olefins and π-systems. Previous results from the Maulide group demonstrated an unexpected propensity for an alternative intramolecular Belluš-Claisen-type rearrangement rather than an expected intramolecular (2 + 2) cycloaddition. We have conducted a cooperative density functional theory/experimental investigation of this process, seeking insights into the competition between the observed Claisen-type reaction and the historically expected (2 + 2) cyclization. Our calculations revealed a surprisingly small difference in the free energy barrier between these two intramolecular reactions. Further theoretical and experimental investigations probe the electronics of the substrate, rationalize a competing deallylation side reaction, and demonstrate the proof-of-concept for an enantioselective (2 + 2) variant.

The thermal dimerization of cycloheptatriene is predicted to occur by a concerted [6 + 4] cycloaddition via an ambimodal [6 + 4]/[4 + 6] transition state (TS) and a competing stepwise diradical (6 + 2) cycloaddition; both dimers subsequently undergo intramolecular [4 + 2] cycloadditions to afford thermally stable tetracyclic products. The ambimodal TS is the 10π-electron version of the prototype bispericyclic dimerization of cyclopentadiene discovered by Caramella et al. in 2002. Quantum mechanical studies using several common DFT functionals and post-HF methods, ωB97X-D, M06-2X, DLPNO-CCSD(T), NEVPT2, and PWPB95-D3(BJ), and quasiclassical molecular dynamics simulations provide details of bond timing and bifurcation pathways. By comparing the ambimodal [6 + 4]/[4 + 6] TS for cycloheptatriene dimerization with the ambimodal [4 + 2]/[2 + 4] TS of cyclopentadiene dimerization, we found that the high distortion energy in cycloheptatriene dimerization is the key to its relatively high energy barrier. The computational investigations were coupled with experimental studies of the cycloheptatriene dimerization, which resulted in the isolation of the two tetracyclic dimers. At lower temperature, the product from the predicted exo-[6 + 4]/[4 + 6] cycloaddition, followed by a subsequent intramolecular [4 + 2] cycloaddition, predominantly forms, while at higher temperature, the diradical (6 + 2) cycloadduct is the major product.

ConspectusThis Account describes our quest to understand and predict organic reactivity, a principal goal of physical and theoretical organic chemistry. The focus is on the development and testing of models for the prediction of cycloaddition reactivities and selectivities. We describe the involvement of the Houk group, and other groups, in the evolution of theoretical models that can achieve ever greater accuracy as well as provide practical heuristic models for understanding and prediction.Is the venerable frontier molecular orbital (FMO) model, the basis of Kenichi Fukui's 1981 Nobel Prize, still useful, or must it be replaced with more advanced models? In particular, models such as Conceptual Density Functional, the Pauli Exclusion Model, and the recent popularity of Electrostatic Potential Plots and Dispersion Energies have not only added to our understanding, but they have also created uncertainty about whether the simple FMO heuristic model has a place in 21st century discussions. This Account addresses this issue and asserts the value of the FMO model.Beginning with brief descriptions of selected models for cycloaddition reactivity starting with early donor-acceptor (nucleophile-electrophile) charge-transfer concepts, this Account reviews Fukui's frontier molecular orbital model, Salem and Klopman's orbital, electrostatic and Pauli repulsion model, the conceptual DFT model by Parr and later by Domingo and others, the recent Houk and Bickelhaupt Distortion/Interaction Activation Strain model, and the Bickelhaupt-Hamlin's Pauli-repulsion lowering model.Computations and analyses of four well-studied Diels-Alder cycloadditions, both normal and inverse electron-demand types, are presented. Most were studied earlier in our published work but are presented here with new insights from calculations with modern methods. Depending on the types of substrates (cycloaddends), the dominant factors controlling reactivity can be orbital interactions, electrostatics and polarization, or Pauli repulsion and dispersion effects, or a combination of all of these.By comparing orbital interactions, especially the frontier molecular orbital interactions, with the other factors that influence reactivity, we show why the FMO model is such a powerful─and theoretically meaningful─heuristic for understanding and predicting reactivity. We also present a method to understand Pauli repulsion effects on activation barriers, ρ(1.1). The use of a new reaction coordinate, the extent of Pauli repulsion along the reaction path, is advocated to emphasize the role of repulsive occupied orbital interactions on reactivity.Fukui's frontier molecular orbital model is effective because FMO interactions parallel all the quantities that influence reactivity. The FMO model continues to provide a practical model to understand and guide experiments.

1,2,3-Triazines and 1,2,3,5-tetrazines react rapidly, efficiently, and selectively with amidines to form pyrimidines/1,3,5-triazines, exhibiting an orthogonal reactivity with 1,2,4,5-tetrazine-based conjugation chemistry. Whereas the mechanism of the reaction of the isomeric 1,2,4-triazines and 1,2,4,5-tetrazines with alkenes is well understood, the mechanism of the 1,2,3-triazine/1,2,3,5-tetrazine-amidine reaction as well as its intrinsic reactivity remains underexplored. By using 15N-labeling, kinetic investigations, and kinetic isotope effect studies, complemented by extensive computational investigations, we show that this reaction proceeds through an addition/N2 elimination/cyclization pathway, rather than the generally expected concerted or stepwise Diels-Alder/retro Diels-Alder sequence. The rate-limiting step in this transformation is the initial nucleophilic attack of an amidine on azine C4, with a subsequent energetically favored N2 elimination step compared with a disfavored stepwise formation of a Diels-Alder cycloadduct. The proposed reaction mechanism is in agreement with experimental and computational results, which explains the observed reactivity of 1,2,3-triazines and 1,2,3,5-tetrazines with amidines.

Spin-center shift (SCS) is a radical process involving 1,2-radical translocation along with a two-electron ionic movement, such as elimination of an adjacent leaving group. Such a process was initially observed in some important biochemical transformations, and the unique property has also attracted considerable interest in synthetic chemistry. Experimental, kinetic, as well as computational studies have been performed, and a series of useful radical transformations have been developed and applied in organic synthesis based on SCS processes in the last 20 years. This Perspective is an overview of radical transformations involving the SCS mechanism.

Density functional theory calculations (ωB97X-D) are reported for the reactions of methoxy

In nature and synthetic chemistry, stereoselective [2+1] cyclopropanation is the most prevalent strategy for the synthesis of chiral cyclopropanes, a class of key pharmacophores in pharmaceuticals and bioactive natural products. One of the most extensively studied reactions in the organic chemist’s arsenal, stereoselective [2+1] cyclopropanation, largely relies on the use of stereodefined olefins, which require elaborate laboratory synthesis or tedious separation to ensure high stereoselectivity. Here we report engineered hemoproteins derived from a bacterial cytochrome P450 that catalyze the synthesis of chiral 1,2,3-polysubstituted cyclopropanes, regardless of the stereopurity of the olefin substrates used. Cytochrome P450BM3 variant IC-G3 exclusively converts (Z)-enol acetates to enantio- and diastereoenriched cyclopropanes and in our model reaction delivers a leftover (E)-enol acetate with 98% stereopurity, using whole Escherichia coli cells. IC-G3 was further engineered with a single mutation to enable the biotransformation of (E)-enol acetates to α-branched ketones with high levels of enantioselectivity while simultaneously catalyzing the cyclopropanation of (Z)-enol acetates with excellent activities and selectivities. We conducted docking studies and molecular dynamics simulations to understand how active-site residues distinguish between the substrate isomers and enable the enzyme to perform these distinct transformations with such high selectivities. Computational studies suggest the observed enantio- and diastereoselectivities are achieved through a stepwise pathway. These biotransformations streamline the synthesis of chiral 1,2,3-polysubstituted cyclopropanes from readily available mixtures of (Z/E)-olefins, adding a new dimension to classical cyclopropanation methods.

Comprehensive computational studies were carried out to explore the mechanisms of enantioselective Cu/Pd and stereodivergent Cu/Ir dual-catalytic syntheses of α,α-disubstituted α-amino acids (α-AAs). A chiral copper azomethine ylide undergoes facile α-allylation with racemic π-allylpalladium species or stereopure π-allyliridium complex to stereoconvergently or stereodivergently furnish single/double stereocenters, respectively. Stereoselectivity at the α-center is controlled by the facial selectivity of α-allylation with respect to the prochiral nucleophile. Despite apparently similar transition-state assemblies, computational models and distortion/interaction analyses disclose versatile modes of stereoinduction wherein the copper azomethine ylide species can face-selectively intercept metal-π-allyl intermediates utilizing attractive dispersion interactions and/or sterically caused distortions. Generation of the β-stereocenter in the Cu/Ir system relies on a stereospecifically generated allyliridium complex and electronically controlled branched-to-linear selectivity, while the dual Cu/Pd system yields a linear monochiral product due to steric factors and π-π stacking interactions. The studies demonstrate on a molecular level how ligand-encoded chiral information is transferred to the α-/β-sites of the resulting α-AAs and how the mode of regio-/stereoselection is altered by differences in transition-metal-stabilized coupling partners. To facilitate studies of stereoselective catalysis, a suite of analytical tools to extract controlling factors for asymmetric induction is demonstrated.

The catalysis of the Hajos–Parrish reaction bycis- andtrans-4,5-ethano-proline was explored experimentally and computationally with DFT (ωB97X-D and MN15) and DLPNO-CCSD(T).

1,2-Azaborines represent a unique class of benzene isosteres that holds great potential for various applications. However, it remains a long-standing challenge to prepare monocyclic 1,2-azaborines in an efficient and modular manner. Here we report a straightforward method to directly access diverse multi-substituted 1,2-azaborines from readily available cyclopropyl imines/ketones and dibromoboranes under relatively mild conditions. The reaction is scalable, shows a broad substrate scope, and tolerates a range of functional groups. The utility of this method is demonstrated in the concise syntheses of BN isosteres of a PD-1/PD-L1 inhibitor and pyrethroid insecticide bifenthrin. Combined experimental and computational mechanistic studies suggest that the reaction pathway involves boron-mediated cyclopropane ring-opening and base-mediated elimination, followed by an unusual low-barrier 6π-electrocyclization accelerated by the BN/CC isomerism.

Lignin valorization is being intensely pursued via tandem catalytic depolymerization and biological funneling to produce single products. In many lignin depolymerization processes, aromatic dimers and oligomers linked by carbon–carbon bonds remain intact, necessitating the development of enzymes capable of cleaving these compounds to monomers. Recently, the catabolism of erythro -1,2-diguaiacylpropane-1,3-diol ( erythro -DGPD), a ring-opened lignin-derived β-1 dimer, was reported in Novosphingobium aromaticivorans . The first enzyme in this pathway, LdpA (formerly LsdE), is a member of the nuclear transport factor 2 (NTF-2)-like structural superfamily that converts erythro -DGPD to lignostilbene through a heretofore unknown mechanism. In this study, we performed biochemical, structural, and mechanistic characterization of the N. aromaticivorans LdpA and another homolog identified in Sphingobium sp. SYK-6, for which activity was confirmed in vivo. For both enzymes, we first demonstrated that formaldehyde is the C 1 reaction product, and we further demonstrated that both enantiomers of erythro -DGPD were transformed simultaneously, suggesting that LdpA, while diastereomerically specific, lacks enantioselectivity. We also show that LdpA is subject to a severe competitive product inhibition by lignostilbene. Three-dimensional structures of LdpA were determined using X-ray crystallography, including substrate-bound complexes, revealing several residues that were shown to be catalytically essential. We used density functional theory to validate a proposed mechanism that proceeds via dehydroxylation and formation of a quinone methide intermediate that serves as an electron sink for the ensuing deformylation. Overall, this study expands the range of chemistry catalyzed by the NTF-2-like protein family to a prevalent lignin dimer through a cofactorless deformylation reaction.

Semi-transparent organic photovoltaics (OPVs) are an emerging solar-energy-harvesting technology with promising applications, such as rooftop energy supplies for environmentally friendly greenhouses. However, the poor operational stability of OPVs poses challenges to their feasibility as incessantly serving facilities. Here we report a reductive interlayer structure for semi-transparent OPVs that improves the operational stability of OPVs under continuous solar radiation. The interlayer effectively suppresses the generation of radicals from the electron transport layer under sunlight and prevents the structural decomposition of the organic photoactive layer during operation. The defects that serve as the charge carrier recombination sites are nullified by the electron-donating functional groups of the reduced molecules, which improves photovoltaic performance. The semi-transparent OPVs demonstrate a power conversion efficiency of 13.5% and an average visible transmittance of 21.5%, with remarkable operational stability (84.8% retention after 1,008 h) under continuous illumination. Greenhouse results show that the semi-transparent OPV roof benefits the survival rate and growth of the crops, indicating the importance of our approach in addressing food and energy challenges.

International audience; We report a Rh(III)-catalyzed ortho-C-H bond functionalization of nitroarenes with 1,2-diarylalkynes and carboxylic anhydrides. The reaction unpredictably affords 3,3-disubstituted oxindoles with the formal reduction of the nitro group under redox-neutral conditions. Besides good functional group tolerance, this transformation allows the preparation of oxindoles with a quaternary carbon stereocenter using nonsymmetrical 1,2-diarylalkynes. This protocol is facilitated by the use of a functionalized cyclopentadienyl (Cp(TMP)*)Rh(III) [Cp(TMP)* = 1-(3,4,5-trimethoxyphenyl)-2,3,4,5-tetramethylcyclopentadienyl] catalyst we developed, which combines an electron-rich character with an elliptical shape. Mechanistic investigations, including the isolation of three rhodacyle intermediates and extensive density functional theory calculations, indicate that the reaction proceeds through nitrosoarene intermediates via a cascade of C-H bond activation─O-atom transfer─[1,2]-aryl shift─deoxygenation─N-acylation.

The first examples of cross-metathesis between two different allenes is disclosed. First- and second-generation Ru complexes were found to be ineffective, at most affording only oligomeric products. The exception was a first-generation complex bearing a bidentate phenyl isopropoxy ligand (i.e., PCy3 is not released upon initiation), reactions with which afforded a 1,3-disubstituted allenyl boronate in 22% yield. On the basis of mechanistic studies designed to gain deeper understanding of the reasons for the ineffectiveness of different Ru catalysts, it was discovered that phosphine-free Ru-CAAC complexes have the steric and electronic attributes to be highly effective. The results of these investigations pave the way for development of additional olefin metathesis reactions that generate allenes.

We report a computational investigation of two reaction cascades occurring following the Claisen rearrangements of aryl propargyl ethers to the alternate ortho positions in unsymmetrical reactants. Our computations explain how substituents influence reactivity and regioselectivity. Rearrangement to the substituted ortho carbon leads to a tricyclo[3.2.1.0]octane core, while rearrangement to an unsubstituted ortho carbon leads to a benzopyran. Density functional theory with ωB97X-D indicates that these reactions involve rate-determining Claisen rearrangements followed by subsequent reaction cascades of the Claisen rearrangement products depending on the presence or absence of a substituent at the ortho carbon.

Hypervalent chloranes are a class of rare and poorly explored reagents. Their unique electronic properties confer reactivity that is complementary to that of the common iodanes and emerging bromanes. Highly chemo- and regioselective, metal-free, and mild C-C and C-O couplings are reported here. Experimental and computational mechanistic studies elucidate the unprecedented reactivities and selectivities of these systems and the intermediacy of aryne intermediates. The synthetic potential of these transformations is further demonstrated

From the preparation of pharmaceuticals to enzymatic construction of natural products, carbocations are central to molecular synthesis. Although these reactive intermediates are engaged in stereoselective processes in nature, exerting enantiocontrol over carbocations with synthetic catalysts remains challenging. Many resonance-stabilized tricoordinated carbocations, such as iminium and oxocarbenium ions, have been applied in catalytic enantioselective reactions. However, their dicoordinated counterparts (aryl and vinyl carbocations) have not, despite their emerging utility in chemical synthesis. We report the discovery of a highly enantioselective vinyl carbocation carbon–hydrogen (C–H) insertion reaction enabled by imidodiphosphorimidate organocatalysts. Active site confinement featured in this catalyst class not only enables effective enantiocontrol but also expands the scope of vinyl cation C–H insertion chemistry, which broadens the utility of this transition metal–free C(sp 3 )–H functionalization platform.

Pericyclic reactions, which involve cyclic concerted transition states without ionic or radical intermediates, have been extensively studied since their definition in the 1960s, and the famous Woodward-Hoffmann rules predict their stereoselectivity and chemoselectivity. Here, we describe the application of a fully automated reaction-path search method, that is, the artificial force induced reaction (AFIR), to trace an input compound back to reasonable starting materials through thermally allowed pericyclic reactions via product-based quantum-chemistry-aided retrosynthetic analysis (QCaRA) without using any a priori experimental knowledge. All categories of pericyclic reactions, including cycloadditions, ene reactions, group-transfer, cheletropic, electrocyclic, and sigmatropic reactions, were successfully traced back via concerted reaction pathways, and starting materials were computationally obtained with the correct stereochemistry. Furthermore, AFIR was used to predict whether the identified reaction pathway can be expected to occur in good yield relative to other possible reactions of the identified starting material. In order to showcase its practical utility, this state-of-the-art technology was also applied to the retrosynthetic analysis of a natural product with a relatively high number of atoms (52 atoms: endiandric acid C methyl ester), which was first synthesized by Nicolaou in 1982 and provided the corresponding starting polyenes with the correct stereospecificity via three pericyclic reaction cascades (one Diels-Alder reaction as well as 6π and 8π electrocyclic reactions). Moreover, not only systems that obey the Woodward-Hoffmann rules but also systems that violate these rules, such as those recently calculated by Houk, can be retrosynthesized accurately.