Autonomous closed-loop mechanistic investigation of molecular electrochemistry via automation

Electrochemical research often requires stringent combinations of experimental parameters that are demanding to manually locate. Recent advances in automated instrumentation and machine-learning algorithms unlock the possibility for accelerated studies of electrochemical fundamentals via high-throughput, online decision-making. Here we report an autonomous electrochemical platform that implements an adaptive, closed-loop workflow for mechanistic investigation of molecular electrochemistry. As a proof-of-concept, this platform autonomously identifies and investigates an EC mechanism, an interfacial electron transfer (E step) followed by a solution reaction (C step), for cobalt tetraphenylporphyrin exposed to a library of organohalide electrophiles. The generally applicable workflow accurately discerns the EC mechanism's presence amid negative controls and outliers, adaptively designs desired experimental conditions, and quantitatively extracts kinetic information of the C step spanning over 7 orders of magnitude, from which mechanistic insights into oxidative addition pathways are gained. This work opens opportunities for autonomous mechanistic discoveries in self-driving electrochemistry laboratories without manual intervention.

Enantioselective Sulfonimidamide Acylation via a Cinchona Alkaloid-Catalyzed Desymmetrization: Scope, Data Science, and Mechanistic Investigation

Methods to access chiral sulfur(VI) pharmacophores are of interest in medicinal and synthetic chemistry. We report the desymmetrization of unprotected sulfonimidamides via asymmetric acylation with a cinchona-phosphinate catalyst. The desired products are formed in excellent yield and enantioselectivity with no observed bis-acylation. A data-science-driven approach to substrate scope evaluation was coupled to high throughput experimentation (HTE) to facilitate statistical modeling in order to inform mechanistic studies. Reaction kinetics, catalyst structural studies, and density functional theory (DFT) transition state analysis elucidated the turnover-limiting step to be the collapse of the tetrahedral intermediate and provided key insights into the catalyst-substrate structure-activity relationships responsible for the origin of the enantioselectivity. This study offers a reliable method for accessing enantioenriched sulfonimidamides to propel their application as pharmacophores and serves as an example of the mechanistic insight that can be gleaned from integrating data science and traditional physical organic techniques.

Electrocatalytic Asymmetric Nozaki-Hiyama-Kishi Decarboxylative Coupling: Scope, Applications, and Mechanism

The first general enantioselective alkyl-Nozaki-Hiyama-Kishi (NHK) coupling reactions are disclosed herein by employing a Cr-electrocatalytic decarboxylative approach. Using easily accessible aliphatic carboxylic acids (via redox-active esters) as alkyl nucleophile synthons, in combination with aldehydes and enabling additives, chiral secondary alcohols are produced in a good yield with high enantioselectivity under mild reductive electrolysis. This reaction, which cannot be mimicked using stoichiometric metal or organic reductants, tolerates a broad range of functional groups and is successfully applied to dramatically simplify the synthesis of multiple medicinally relevant structures and natural products. Mechanistic studies revealed that this asymmetric alkyl e-NHK reaction was enabled by using catalytic tetrakis(dimethylamino)ethylene, which acts as a key reductive mediator to mediate the electroreduction of the CrIII/chiral ligand complex.

14N to 15N Isotopic Exchange of Nitrogen Heteroaromatics through Skeletal Editing

The selective modification of nitrogen heteroaromatics enables the development of new chemical tools and accelerates drug discovery. While methods that focus on expanding or contracting the skeletal structures of heteroaromatics are emerging, methods for the direct exchange of single core atoms remain limited. Here, we present a method for 14N → 15N isotopic exchange for several aromatic nitrogen heterocycles. This nitrogen isotope transmutation occurs through activation of the heteroaromatic substrate by triflylation of a nitrogen atom, followed by a ring-opening/ring-closure sequence mediated by 15N-aspartate to effect the isotopic exchange of the nitrogen atom. Key to the success of this transformation is the formation of an isolable 15N-succinyl intermediate, which undergoes elimination to give the isotopically labeled heterocycle. These transformations occur under mild conditions in high chemical and isotopic yields.

Computational Methods Enable the Prediction of Improved Catalysts for Nickel-Catalyzed Cross-Electrophile Coupling

Cross-electrophile coupling has emerged as an attractive and efficient method for the synthesis of C(sp2)-C(sp3) bonds. These reactions are most often catalyzed by nickel complexes of nitrogenous ligands, especially 2,2'-bipyridines. Precise prediction, selection, and design of optimal ligands remains challenging, despite significant increases in reaction scope and mechanistic understanding. Molecular parameterization and statistical modeling provide a path to the development of improved bipyridine ligands that will enhance the selectivity of existing reactions and broaden the scope of electrophiles that can be coupled. Herein, we describe the generation of a computational ligand library, correlation of observed reaction outcomes with features of the ligands, and the in silico design of improved bipyridine ligands for Ni-catalyzed cross-electrophile coupling. The new nitrogen-substituted ligands display a 5-fold increase in selectivity for product formation versus homodimerization when compared to the current state of the art. This increase in selectivity and yield was general for several cross-electrophile couplings, including the challenging coupling of an aryl chloride with an N-alkylpyridinium salt.

Predicting success in Cu-catalyzed C-N coupling reactions using data science

Data science is assuming a pivotal role in guiding reaction optimization and streamlining experimental workloads in the evolving landscape of synthetic chemistry. A discipline-wide goal is the development of workflows that integrate computational chemistry and data science tools with high-throughput experimentation as it provides experimentalists the ability to maximize success in expensive synthetic campaigns. Here, we report an end-to-end data-driven process to effectively predict how structural features of coupling partners and ligands affect Cu-catalyzed C-N coupling reactions. The established workflow underscores the limitations posed by substrates and ligands while also providing a systematic ligand prediction tool that uses probability to assess when a ligand will be successful. This platform is strategically designed to confront the intrinsic unpredictability frequently encountered in synthetic reaction deployment.

Deconvoluting Nonlinear Catalyst-Substrate Effects in the Intramolecular Dirhodium-Catalyzed C-H Insertion of Donor/Donor Carbenes Using Data Science Tools

Interactions between catalysts and substrates can be highly complex and dynamic, often complicating the development of models to either predict or understand such processes. A dirhodium(II)-catalyzed C-H insertion of donor/donor carbenes into 2-alkoxybenzophenone substrates to form benzodihydrofurans was selected as a model system to explore nonlinear methods to achieve a mechanistic understanding. We found that the application of traditional methods of multivariate linear regression (MLR) correlating DFT-derived descriptors of catalysts and substrates leads to poorly performing models. This inspired the introduction of nonlinear descriptor relationships into modeling by applying the sure independence screening and sparsifying operator (SISSO) algorithm. Based on SISSO-generated descriptors, a high-performing MLR model was identified that predicts external validation points well. Mechanistic interpretation was aided by the deconstruction of feature relationships using chemical space maps, decision trees, and linear descriptors. Substrates were found to have a strong dependence on steric effects for determining their innate cyclization selectivity preferences. Catalyst reactive site features can then be matched to product features to tune or override the resultant diastereoselectivity within the substrate-dictated ranges. This case study presents a method for understanding complex interactions often encountered in catalysis by using nonlinear modeling methods and linear deconvolution by pattern recognition.

Dataset Design for Building Models of Chemical Reactivity

Models can codify our understanding of chemical reactivity and serve a useful purpose in the development of new synthetic processes via, for example, evaluating hypothetical reaction conditions or in silico substrate tolerance. Perhaps the most determining factor is the composition of the training data and whether it is sufficient to train a model that can make accurate predictions over the full domain of interest. Here, we discuss the design of reaction datasets in ways that are conducive to data-driven modeling, emphasizing the idea that training set diversity and model generalizability rely on the choice of molecular or reaction representation. We additionally discuss the experimental constraints associated with generating common types of chemistry datasets and how these considerations should influence dataset design and model building.

Data science enabled discovery of a highly soluble 2,2'-bipyrimidine anolyte for application in a flow battery

Development of non-aqueous redox flow batteries as a viable energy storage solution relies upon the identification of soluble charge carriers capable of storing large amounts of energy over extended time periods. A combination of metrics including number of electrons stored per molecule, redox potential, stability, and solubility of the charge carrier impact performance. In this context, we recently reported a 2,2'-bipyrimidine charge carrier that stores two electrons per molecule with reduction near -2.0 V vs. Fc/Fc+ and high stability. However, these first-generation derivatives showed a modest solubility of 0.17 M (0.34 M e-). Seeking to improve solubility without sacrificing stability, we harnessed the synthetic modularity of this scaffold to design a library of sixteen candidates. Using computed molecular descriptors and a single node decision tree, we found that minimization of the solvent accessible surface area (SASA) can be used to predict derivatives with enhanced solubility. This parameter was used in combination with a heatmap describing stability to de-risk a virtual screen that ultimately identified a 2,2'-bipyrimidine with significantly increased solubility and good stability metrics in the reduced states. This molecule was paired with a cyclopropenium catholyte in a prototype all-organic redox flow battery, achieving a cell potential up to 3 V.

Identifying structure-function relationships to modulate crossover in nonaqueous redox flow batteries

Nonaqueous redox flow batteries (NARFBs) offer a promising solution for large-scale storage of renewable energy. However, crossover of redox active molecules between the two sides of the cell is a major factor limiting their development, as most selective separators are designed for deployment in water, rather than organic solvents. This report describes a systematic investigation of the crossover rates of redox active organic molecules through an anion exchange separator under RFB-relevant non-aqueous conditions (in acetonitrile/KPF6) using a combination of experimental and computational methods. A structurally diverse set of neutral and cationic molecules was selected, and their rates of crossover were determined experimentally with the organic solvent-compatible anion exchange separator Fumasep FAP-375-PP. The resulting data were then fit to various descriptors of molecular size, charge, and hydrophobicity (overall charge, solution diffusion coefficient, globularity, dynamic volume, dynamic surface area, clogP). This analysis resulted in multiple statistical models of crossover rates for this separator. These models were then used to predict tether groups that dramatically slow the crossover of small organic molecules in this system.

An Asymmetric Aromatic Finkelstein Reaction: A Platform for Remote Diarylmethane Desymmetrization

A first-of-its-kind enantioselective aromatic Finkelstein reaction is disclosed for the remote desymmetrization of diarylmethanes. The reaction operates through a copper-catalyzed C-I bond-forming event, and high levels of enantioselectivity are achieved through the deployment of a tailored guanidinylated peptide ligand. Strategic use of transition-metal-mediated reactions enables the chemoselective modification of the aryl iodide products; thus, the synthesis of a diverse set of otherwise difficult-to-access diarylmethanes with excellent levels of selectivity is realized from a common intermediate. A mixed experimental/computational analysis of steric parameters and substrate conformations identifies the importance of remote conformational effects as a key to achieving high enantioselectivity in this desymmetrization reaction.

Data Science-Enabled Palladium-Catalyzed Enantioselective Aryl-Carbonylation of Sulfonimidamides

New methods for the general asymmetric synthesis of sulfonimidamides are of great interest due to their applications in medicinal chemistry, agrochemical discovery, and academic research. We report a palladium-catalyzed cross-coupling method for the enantioselective aryl-carbonylation of sulfonimidamides. Using data science techniques, a virtual library of calculated bisphosphine ligand descriptors was used to guide reaction optimization by effectively sampling the catalyst chemical space. The optimized conditions identified using this approach provided the desired product in excellent yield and enantioselectivity. As the next step, a data science-driven strategy was also used to explore a diverse set of aryl and heteroaryl iodides, providing key information about the scope and limitations of the method. Furthermore, we tested a range of racemic sulfonimidamides for compatibility of this coupling partner. The developed method offers a general and efficient strategy for accessing enantioenriched sulfonimidamides, which should facilitate their application in industrial and academic settings.

Using Data Science for Mechanistic Insights and Selectivity Predictions in a Non-Natural Biocatalytic Reaction

The study of non-natural biocatalytic transformations relies heavily on empirical methods, such as directed evolution, for identifying improved variants. Although exceptionally effective, this approach provides limited insight into the molecular mechanisms behind the transformations and necessitates multiple protein engineering campaigns for new reactants. To address this limitation, we disclose a strategy to explore the biocatalytic reaction space and garner insight into the molecular mechanisms driving enzymatic transformations. Specifically, we explored the selectivity of an "ene"-reductase, GluER-T36A, to create a data-driven toolset that explores reaction space and rationalizes the observed and predicted selectivities of substrate/mutant combinations. The resultant statistical models related structural features of the enzyme and substrate to selectivity and were used to effectively predict selectivity in reactions with out-of-sample substrates and mutants. Our approach provided a deeper understanding of enantioinduction by GluER-T36A and holds the potential to enhance the virtual screening of enzyme mutants.

Data Science Enables the Development of a New Class of Chiral Phosphoric Acid Catalysts

The widespread success of BINOL-chiral phosphoric acids (CPAs) has led to the development of several high molecular weight, sterically encumbered variants. Herein, we disclose an alternative, minimalistic chiral phosphoric acid backbone incorporating only a single instance of point chirality. Data science techniques were used to select a diverse training set of catalysts, which were benchmarked against the transfer hydrogenation of an 8-aminoquinoline. Using a univariate classification algorithm and multivariate linear regression, key catalyst features necessary for high levels of selectivity were deconvoluted, revealing a simple catalyst model capable of predicting selectivity for out-of-set catalysts. This workflow enabled extrapolation to a catalyst providing higher selectivity than both reported peptide-type and BINOL-type catalysts (up to 95:5 er). These techniques were then successfully applied towards two additional transforms. Taken together, these examples illustrate the power of combining rational design with data science (ab initio) to efficiently explore reactivity during catalyst development.

Catalyst-Controlled Enantioselective and Regiodivergent Addition of Aryl Boron Nucleophiles to N-Alkyl Nicotinate Salts

Dihydropyridines are versatile building blocks for the synthesis of pyridines, tetrahydropyridines, and piperidines. Addition of nucleophiles to activated pyridinium salts allows synthesis of 1,2-, 1,4-, or 1,6-dihydropyridines; however, this process often leads to a mixture of constitutional isomers. Catalyst-controlled regioselective addition of nucleophiles to pyridiniums has the potential to solve this problem. Herein, we report that the regioselective addition of boron-based nucleophiles to pyridinium salts can be accomplished by the choice of a Rh catalyst.

Generality-oriented optimization of enantioselective aminoxyl radical catalysis

Catalytic enantioselective methods that are generally applicable to a broad range of substrates are rare. We report a strategy for the oxidative desymmetrization of meso-diols predicated on a nontraditional catalyst optimization protocol by using a panel of screening substrates rather than a singular model substrate. Critical to this approach was rational modulation of a peptide sequence in the catalyst incorporating a distinct aminoxyl-based active residue. A general catalyst emerged, providing high selectivity in the delivery of enantioenriched lactones across a broad range of diols, while also achieving up to ~100,000 turnovers.

A Physical Organic Approach towards Statistical Modeling of Tetrazole and Azide Decomposition

Nitrogen atom-rich heterocycles and organic azides have found extensive use in many sectors of modern chemistry from drug discovery to energetic materials. The prediction and understanding of their energetic properties are thus key to the safe and effective application of these compounds. In this work, we disclose the use of multivariate linear regression modeling for the prediction of the decomposition temperature and impact sensitivity of structurally diverse tetrazoles and organic azides. We report a data-driven approach for property prediction featuring a collection of quantum mechanical parameters and computational workflows. The statistical models reported herein carry predictive accuracy as well as chemical interpretability. Model validation was successfully accomplished via tetrazole test sets with parameters generated exclusively in silico. Mechanistic analysis of the statistical models indicated distinct divergent pathways of thermal and impact-initiated decomposition.

Interrogating the Mechanistic Features of Ni(I)-Mediated Aryl Iodide Oxidative Addition Using Electroanalytical and Statistical Modeling Techniques

While the oxidative addition of Ni(I) to aryl iodides has been commonly proposed in catalytic methods, an in-depth mechanistic understanding of this fundamental process is still lacking. Herein, we describe a detailed mechanistic study of the oxidative addition process using electroanalytical and statistical modeling techniques. Electroanalytical techniques allowed rapid measurement of the oxidative addition rates for a diverse set of aryl iodide substrates and four classes of catalytically relevant complexes (Ni(MeBPy), Ni(MePhen), Ni(Terpy), and Ni(BPP)). With >200 experimental rate measurements, we were able to identify essential electronic and steric factors impacting the rate of oxidative addition through multivariate linear regression models. This has led to a classification of oxidative addition mechanisms, either through a three-center concerted or halogen-atom abstraction pathway based on the ligand type. A global heat map of predicted oxidative addition rates was created and shown applicable to a better understanding of the reaction outcome in a case study of a Ni-catalyzed coupling reaction.

Data-Driven Multi-Objective Optimization Tactics for Catalytic Asymmetric Reactions Using Bisphosphine Ligands

Optimization of the catalyst structure to simultaneously improve multiple reaction objectives (e.g., yield, enantioselectivity, and regioselectivity) remains a formidable challenge. Herein, we describe a machine learning workflow for the multi-objective optimization of catalytic reactions that employ chiral bisphosphine ligands. This was demonstrated through the optimization of two sequential reactions required in the asymmetric synthesis of an active pharmaceutical ingredient. To accomplish this, a density functional theory-derived database of >550 bisphosphine ligands was constructed, and a designer chemical space mapping technique was established. The protocol used classification methods to identify active catalysts, followed by linear regression to model reaction selectivity. This led to the prediction and validation of significantly improved ligands for all reaction outputs, suggesting a general strategy that can be readily implemented for reaction optimizations where performance is controlled by bisphosphine ligands.

Exploring Structure-Function Relationships of Aryl Pyrrolidine-Based Hydrogen-Bond Donors in Asymmetric Catalysis Using Data-Driven Techniques

Hydrogen bond-based organocatalysts rely on networks of attractive noncovalent interactions (NCIs) to impart enantioselectivity. As a specific example, aryl pyrrolidine substituted urea, thiourea, and squaramide organocatalysts function cooperatively through hydrogen bonding and difficult-to-predict NCIs as a function of the reaction partners. To uncover the synergistic effect of the structural components of this catalyst class, we applied data science tools to study various model reactions using a derivatized, aryl pyrrolidine-based, hydrogen-bond donor (HBD) catalyst library. Through a combination of experimentally collected data and data mined from previous reports, statistical models were constructed, illuminating the general features necessary for high enantioselectivity. A distinct dependence on the identity of the electrophilic reaction partner and HBD catalyst is observed, suggesting that a general interaction is conserved throughout the reactions analyzed. The resulting models also demonstrate predictive capability by the successful improvement of a previously reported reaction using out-of-sample reaction components. Overall, this study highlights the power of data science in exploring mechanistic hypotheses in asymmetric HBD catalysis and provides a prediction platform applicable in future reaction optimization.

Catalytic asymmetric C-H insertion reactions of vinyl carbocations

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³)-H functionalization platform.

Data Science-Driven Analysis of Substrate-Permissive Diketopiperazine Reverse Prenyltransferase NotF: Applications in Protein Engineering and Cascade Biocatalytic Synthesis of (-)-Eurotiumin A

Prenyltransfer is an early-stage carbon-hydrogen bond (C-H) functionalization prevalent in the biosynthesis of a diverse array of biologically active bacterial, fungal, plant, and metazoan diketopiperazine (DKP) alkaloids. Toward the development of a unified strategy for biocatalytic construction of prenylated DKP indole alkaloids, we sought to identify and characterize a substrate-permissive C2 reverse prenyltransferase (PT). As the first tailoring event within the biosynthesis of cytotoxic notoamide metabolites, PT NotF catalyzes C2 reverse prenyltransfer of brevianamide F. Solving a crystal structure of NotF (in complex with native substrate and prenyl donor mimic dimethylallyl S-thiolodiphosphate (DMSPP)) revealed a large, solvent-exposed active site, intimating NotF may possess a significantly broad substrate scope. To assess the substrate selectivity of NotF, we synthesized a panel of 30 sterically and electronically differentiated tryptophanyl DKPs, the majority of which were selectively prenylated by NotF in synthetically useful conversions (2 to >99%). Quantitative representation of this substrate library and development of a descriptive statistical model provided insight into the molecular origins of NotF's substrate promiscuity. This approach enabled the identification of key substrate descriptors (electrophilicity, size, and flexibility) that govern the rate of NotF-catalyzed prenyltransfer, and the development of an "induced fit docking (IFD)-guided" engineering strategy for improved turnover of our largest substrates. We further demonstrated the utility of NotF in tandem with oxidative cyclization using flavin monooxygenase, BvnB. This one-pot, in vitro biocatalytic cascade enabled the first chemoenzymatic synthesis of the marine fungal natural product, (-)-eurotiumin A, in three steps and 60% overall yield.

Investigating Oxidative Addition Mechanisms of Allylic Electrophiles with Low-Valent Ni/Co Catalysts Using Electroanalytical and Data Science Techniques

The catalysis by a π-allyl-Co/Ni complex has drawn significant attention recently due to its distinct reactivity in reductive Co/Ni-catalyzed allylation reactions. Despite significant success in reaction development, the critical oxidative addition mechanism to form the π-allyl-Co/Ni complex remains unclear. Herein, we present a study to investigate this process with four catalysis-relevant complexes: Co(^MeBPy)Br₂, Co(^MePhen)Br₂, Ni(^MeBPy)Br₂, and Ni(^MePhen)Br₂. Enabled by an electroanalytical platform, Co(I)/Ni(I) species were found responsible for the oxidative addition of allyl acetate. Kinetic features of different substrates were characterized through linear free-energy relationship (Hammett-type) studies, statistical modeling, and a DFT computational study. In this process, a coordination-ionization-type transition state was proposed, sharing a similar feature with Pd(0)-mediated oxidative addition in Tsuji-Trost reactions. Computational and ligand structural analysis studies support this mechanism, which should provide key information for next-generation catalyst development.

Noncovalent Stabilization of Radical Intermediates in the Enantioselective Hydroamination of Alkenes with Sulfonamides

Noncovalent interactions (NCIs) are critical elements of molecular recognition in a wide variety of chemical contexts. While NCIs have been studied extensively for closed-shell molecules and ions, very little is understood about the structures and properties of NCIs involving free radical intermediates. In this report, we describe a detailed mechanistic study of the enantioselective radical hydroamination of alkenes with sulfonamides and present evidence suggesting that the basis for asymmetric induction in this process arises from attractive NCIs between a neutral sulfonamidyl radical intermediate and a chiral phosphoric acid (CPA). We describe experimental, computational, and data science-based evidence that identifies the specific radical NCIs that form the basis for the enantioselectivity. Kinetic studies support that C-N bond formation determines the enantioselectivity. Density functional theory investigations revealed the importance of both strong H-bonding between the CPA and the N-centered radical and a network of aryl-based NCIs that serve to stabilize the favored diastereomeric transition state. The contributions of these specific aryl-based NCIs to the selectivity were further confirmed through multivariate linear regression analysis by comparing the measured enantioselectivity to computed descriptors. These results highlight the power of NCIs to enable high levels of enantioselectivity in reactions involving uncharged open-shell intermediates and expand our understanding of radical-molecule interactions.

Leveraging Regio- and Stereoselective C(sp3)-H Functionalization of Silyl Ethers to Train a Logistic Regression Classification Model for Predicting Site-Selectivity Bias

The C-H functionalization of silyl ethers via carbene-induced C-H insertion represents an efficient synthetic disconnection strategy. In this work, site- and stereoselective C(sp³)-H functionalization at α, γ, δ, and even more distal positions to the siloxy group has been achieved using donor/acceptor carbene intermediates. By exploiting the predilections of Rh₂(R-TCPTAD)₄ and Rh₂(S-2-Cl-5-BrTPCP)₄ catalysts to target either more electronically activated or more spatially accessible C-H sites, respectively, divergent desired products can be formed with good diastereocontrol and enantiocontrol. Notably, the reaction can also be extended to enable desymmetrization of meso silyl ethers. Leveraging the broad substrate scope examined in this study, we have trained a machine learning classification model using logistic regression to predict the major C-H functionalization site based on intrinsic substrate reactivity and catalyst propensity for overriding it. This model enables prediction of the major product when applying these C-H functionalization methods to a new substrate of interest. Applying this model broadly, we have demonstrated its utility for guiding late-stage functionalization in complex settings and developed an intuitive visualization tool to assist synthetic chemists in such endeavors.

Cobalt-electrocatalytic HAT for functionalization of unsaturated C-C bonds

The study and application of transition metal hydrides (TMHs) has been an active area of chemical research since the early 1960s¹, for energy storage, through the reduction of protons to generate hydrogen^2,3, and for organic synthesis, for the functionalization of unsaturated C-C, C-O and C-N bonds^4,5. In the former instance, electrochemical means for driving such reactivity has been common place since the 1950s⁶ but the use of stoichiometric exogenous organic- and metal-based reductants to harness the power of TMHs in synthetic chemistry remains the norm. In particular, cobalt-based TMHs have found widespread use for the derivatization of olefins and alkynes in complex molecule construction, often by a net hydrogen atom transfer (HAT)⁷. Here we show how an electrocatalytic approach inspired by decades of energy storage research can be made use of in the context of modern organic synthesis. This strategy not only offers benefits in terms of sustainability and efficiency but also enables enhanced chemoselectivity and distinct, tunable reactivity. Ten different reaction manifolds across dozens of substrates are exemplified, along with detailed mechanistic insights into this scalable electrochemical entry into Co-H generation that takes place through a low-valent intermediate.

Predicting relative efficiency of amide bond formation using multivariate linear regression

Amides are ubiquitous in biologically active natural products and commercial drugs. The most common strategy for introducing this functional group is the coupling of a carboxylic acid with an amine, which requires the use of a coupling reagent to facilitate elimination of water. However, the optimal reaction conditions often appear rather arbitrary to the specific reaction. Herein, we report the development of statistical models correlating measured rates to physical organic descriptors to enable the prediction of reaction rates for untested carboxylic acid/amine pairs. The key to the success of this endeavor was the development of an end-to-end data science–based workflow to select a set of coupling partners that are appropriately distributed in chemical space to facilitate statistical model development. By using a parameterization, dimensionality reduction, and clustering protocol, a training set was identified. Reaction rates for a range of carboxylic acid and primary alkyl amine couplings utilizing carbonyldiimidazole (CDI) as the coupling reagent were measured. The collected rates span five orders of magnitude, confirming that the designed training set encompasses a wide range of chemical space necessary for effective model development. Regressing these rates with high-level density functional theory (DFT) descriptors allowed for identification of a statistical model wherein the molecular features of the carboxylic acid are primarily responsible for the observed rates. Finally, out-of-sample amide couplings are used to determine the limitations and effectiveness of the model.

A Comprehensive Discovery Platform for Organophosphorus Ligands for Catalysis

The design of molecular catalysts typically involves reconciling multiple conflicting property requirements, largely relying on human intuition and local structural searches. However, the vast number of potential catalysts requires pruning of the candidate space by efficient property prediction with quantitative structure-property relationships. Data-driven workflows embedded in a library of potential catalysts can be used to build predictive models for catalyst performance and serve as a blueprint for novel catalyst designs. Herein we introduce kraken, a discovery platform covering monodentate organophosphorus(III) ligands providing comprehensive physicochemical descriptors based on representative conformer ensembles. Using quantum-mechanical methods, we calculated descriptors for 1558 ligands, including commercially available examples, and trained machine learning models to predict properties of over 300000 new ligands. We demonstrate the application of kraken to systematically explore the property space of organophosphorus ligands and how existing data sets in catalysis can be used to accelerate ligand selection during reaction optimization.

Stereoconvergent and -divergent Synthesis of Tetrasubstituted Alkenes by Nickel-Catalyzed Cross-Couplings

We report the development of a method to diastereoselectively access tetrasubstituted alkenes via nickel-catalyzed Suzuki-Miyaura cross-couplings of enol tosylates and boronic acid esters. Either diastereomeric product was selectively accessed from a mixture of enol tosylate starting material diastereomers in a convergent reaction by judicious choice of the ligand and reaction conditions. A similar protocol also enabled a divergent synthesis of each product isomer from diastereomerically pure enol tosylates. Notably, high-throughput optimization of the monophosphine ligands was guided by chemical space analysis of the kraken library to ensure a diverse selection of ligands was examined. Stereoelectronic analysis of the results provided insight into the requirements for reactive and selective ligands in this transformation. The synthetic utility of the optimized catalytic system was then probed in the stereoselective synthesis of various tetrasubstituted alkenes, with yields up to 94% and diastereomeric ratios up to 99:1 Z/E and 93:7 E/Z observed. Moreover, a detailed computational analysis and experimental mechanistic studies provided key insights into the nature of the underlying isomerization process impacting selectivity in the cross-coupling.

A Data-Driven Approach to the Development and Understanding of Chiroptical Sensors for Alcohols with Remote γ-Stereocenters

Dynamic covalent chemistry-based sensors have recently emerged as powerful tools to rapidly determine the enantiomeric excess of organic small molecules. While a bevy of sensors have been developed, those for flexible molecules with stereocenters remote to the functional group that binds the chiroptical sensor remain scarce. In this study, we develop an iterative, data-driven workflow to design and analyze a chiroptical sensor capable of assessing challenging acyclic γ-stereogenic alcohols. Following sensor optimization, the mechanism of sensing was probed with a combination of computational parametrization of the sensor molecules, statistical modeling, and high-level density functional theory (DFT) calculations. These were used to elucidate the mechanism of stereochemical recognition and revealed that competing attractive noncovalent interactions (NCIs) determine the overall performance of the sensor. It is anticipated that the data-driven workflows developed herein will be generally applicable to the development and understanding of dynamic covalent and supramolecular sensors.

The Evolution of Data-Driven Modeling in Organic Chemistry

Organic chemistry is replete with complex relationships: for example, how a reactant's structure relates to the resulting product formed; how reaction conditions relate to yield; how a catalyst's structure relates to enantioselectivity. Questions like these are at the foundation of understanding reactivity and developing novel and improved reactions. An approach to probing these questions that is both longstanding and contemporary is data-driven modeling. Here, we provide a synopsis of the history of data-driven modeling in organic chemistry and the terms used to describe these endeavors. We include a timeline of the steps that led to its current state. The case studies included highlight how, as a community, we have advanced physical organic chemistry tools with the aid of computers and data to augment the intuition of expert chemists and to facilitate the prediction of structure-activity and structure-property relationships.

Development of High Energy Density Diaminocyclopropenium-Phenothiazine Hybrid Catholytes for Non-Aqueous Redox Flow Batteries

This report describes the design of diaminocyclopropenium-phenothiazine hybrid catholytes for non-aqueous redox flow batteries. The molecules are synthesized in a rapid and modular fashion by appending a diaminocyclopropenium (DAC) substituent to the nitrogen of the phenothiazine. Combining a versatile C-N coupling protocol (which provides access to diverse derivatives) with computation and structure-property analysis enabled the identification of a catholyte that displays stable two-electron cycling at potentials of 0.64 and 1.00 V vs. Fc/Fc⁺ as well as high solubility in all oxidation states (≥0.45 M in TBAPF₆ /MeCN). This catholyte was deployed in a high energy density two-electron RFB, exhibiting >90 % capacity retention over 266 hours of flow cell cycling at >0.5 M electron concentration.

Univariate classification of phosphine ligation state and reactivity in cross-coupling catalysis

[Figure: see text].

Simultaneously Enhancing the Redox Potential and Stability of Multi-Redox Organic Catholytes by Incorporating Cyclopropenium Substituents

High redox potential, two-electron organic catholytes for nonaqueous redox flow batteries were developed by appending diaminocyclopropenium (DAC) substituents to phenazine and phenothiazine cores. The parent heterocycles exhibit two partially reversible oxidations at moderate potentials [both at lower than 0.7 V vs ferrocene/ferrocenium (Fc/Fc⁺)]. The incorporation of DAC substituents has a dual effect on these systems. The DAC groups increase the redox potential of both couples by ∼300 mV while simultaneously rendering the second oxidation (which occurs at 1.20 V vs Fc/Fc⁺ in the phenothiazine derivative) reversible. The electron-withdrawing nature of the DAC unit is responsible for the increase in redox potential, while the DAC substituents stabilize oxidized forms of the molecules through resonance delocalization of charge and unpaired spin density. These new catholytes were deployed in two-electron redox flow batteries that exhibit voltages of up to 2.0 V and no detectable crossover over 250 cycles.

Data Science Meets Physical Organic Chemistry

ConspectusAt the heart of synthetic chemistry is the holy grail of predictable catalyst design. In particular, researchers involved in reaction development in asymmetric catalysis have pursued a variety of strategies toward this goal. This is driven by both the pragmatic need to achieve high selectivities and the inability to readily identify why a certain catalyst is effective for a given reaction. While empiricism and intuition have dominated the field of asymmetric catalysis since its inception, enantioselectivity offers a mechanistically rich platform to interrogate catalyst-structure response patterns that explain the performance of a particular catalyst or substrate.In the early stages of an asymmetric reaction development campaign, the overarching mechanism of the reaction, catalyst speciation, the turnover limiting step, and many other details are unknown or posited based on related reactions. Considering the unclear details leading to a successful reaction, initial enantioselectivity data are often used to intuitively guide the ultimate direction of optimization. However, if the conditions of the Curtin-Hammett principle are satisfied, then measured enantioselectivity can be directly connected to the ensemble of diastereomeric transition states (TSs) that lead to the enantiomeric products, and the associated free energy difference between competing TSs (ΔΔG⧧ = -RT ln[(S)/(R)], where (S) and (R) represent the concentrations of the enantiomeric products). We, and others, speculated that this important piece of information can be leveraged to guide reaction optimization in a quantitative way.Although traditional linear free energy relationships (LFERs), such as Hammett plots, have been used to illuminate important mechanistic features, we sought to develop data science derived tools to expand the power of LFERs in order to describe complex reactions frequently encountered in modern asymmetric catalysis. Specifically, we investigated whether enantioselectivity data from a reaction can be quantitatively connected to the attributes of reaction components, such as catalyst and substrate structural features, to harness data for asymmetric catalyst design.In this context, we developed a workflow to relate computationally derived features of reaction components to enantioselectivity using data science tools. The mathematical representation of molecules can incorporate many aspects of a transformation, such as molecular features from substrate, product, catalyst, and proposed transition states. Statistical models relating these features to reaction outputs can be used for various tasks, such as performance prediction of untested molecules. Perhaps most importantly, statistical models can guide the generation of mechanistic hypotheses that are embedded within complex patterns of reaction responses. Overall, merging traditional physical organic experiments with statistical modeling techniques creates a feedback loop that enables both evaluation of multiple mechanistic hypotheses and future catalyst design. In this Account, we highlight the evolution and application of this approach in the context of a collaborative program based on chiral phosphoric acid catalysts (CPAs) in asymmetric catalysis.

Carbon Atom Insertion into Pyrroles and Indoles Promoted by Chlorodiazirines

Herein, we report a reaction that selectively generates 3-arylpyridine and quinoline motifs by inserting aryl carbynyl cation equivalents into pyrrole and indole cores, respectively. By employing α-chlorodiazirines as thermal precursors to the corresponding chlorocarbenes, the traditional haloform-based protocol central to the parent Ciamician-Dennstedt rearrangement can be modified to directly afford 3-(hetero)arylpyridines and quinolines. Chlorodiazirines are conveniently prepared in a single step by oxidation of commercially available amidinium salts. Selectivity as a function of pyrrole substitution pattern was examined, and a predictive model based on steric effects is put forward, with DFT calculations supporting a selectivity-determining cyclopropanation step. Computations surprisingly indicate that the stereochemistry of cyclopropanation is of little consequence to the subsequent electrocyclic ring opening that forges the pyridine core, due to a compensatory homoaromatic stabilization that counterbalances orbital-controlled torquoselectivity effects. The utility of this skeletal transform is further demonstrated through the preparation of quinolinophanes and the skeletal editing of pharmaceutically relevant pyrroles.

Data-science driven autonomous process optimization

Autonomous process optimization involves the human intervention-free exploration of a range process parameters to improve responses such as product yield and selectivity. Utilizing off-the-shelf components, we develop a closed-loop system for carrying out parallel autonomous process optimization experiments in batch. Upon implementation of our system in the optimization of a stereoselective Suzuki-Miyaura coupling, we find that the definition of a set of meaningful, broad, and unbiased process parameters is the most critical aspect of successful optimization. Importantly, we discern that phosphine ligand, a categorical parameter, is vital to determination of the reaction outcome. To date, categorical parameter selection has relied on chemical intuition, potentially introducing bias into the experimental design. In seeking a systematic method for selecting a diverse set of phosphine ligands, we develop a strategy that leverages computed molecular feature clustering. The resulting optimization uncovers conditions to selectively access the desired product isomer in high yield.

N-Ammonium Ylide Mediators for Electrochemical C-H Oxidation

The site-specific oxidation of strong C(sp³)-H bonds is of uncontested utility in organic synthesis. From simplifying access to metabolites and late-stage diversification of lead compounds to truncating retrosynthetic plans, there is a growing need for new reagents and methods for achieving such a transformation in both academic and industrial circles. One main drawback of current chemical reagents is the lack of diversity with regard to structure and reactivity that prevents a combinatorial approach for rapid screening to be employed. In that regard, directed evolution still holds the greatest promise for achieving complex C-H oxidations in a variety of complex settings. Herein we present a rationally designed platform that provides a step toward this challenge using N-ammonium ylides as electrochemically driven oxidants for site-specific, chemoselective C(sp³)-H oxidation. By taking a first-principles approach guided by computation, these new mediators were identified and rapidly expanded into a library using ubiquitous building blocks and trivial synthesis techniques. The ylide-based approach to C-H oxidation exhibits tunable selectivity that is often exclusive to this class of oxidants and can be applied to real-world problems in the agricultural and pharmaceutical sectors.

Nickel-catalyzed asymmetric reductive cross-coupling of α-chloroesters with (hetero)aryl iodides

An asymmetric reductive cross-coupling of α-chloroesters and (hetero)aryl iodides is reported. This nickel-catalyzed reaction proceeds with a chiral BiOX ligand under mild conditions, affording α-arylesters in good yields and enantioselectivities. The reaction is tolerant of a variety of functional groups, and the resulting products can be converted to pharmaceutically-relevant chiral building blocks. A multivariate linear regression model was developed to quantitatively relate the influence of the α-chloroester substrate and ligand on enantioselectivity.

A Ni-catalyzed enantioselective reductive cross-coupling of α-chloroesters and (hetero)aryl iodides is reported. A MLR model was developed to quantitatively relate the influence of the α-chloroester substrate and ligand on enantioselectivity.