A Synergistic Bimetallic Ti/Co-Catalyzed Isomerization of Epoxides to Allylic Alcohols Enabled by Two-State Reactivity

Isomerization of epoxides into versatile allylic alcohols is an atom-economical synthetic method to afford vicinal bifunctional groups. Comprehensive density functional theory (DFT) calculations were carried out to elucidate the complex mechanism of a bimetallic Ti/Co-catalyzed selective isomerization of epoxides to allyl alcohols by examining several possible pathways. Our results suggest a possible mechanism involving (1) radical-type epoxide ring opening catalyzed by Cp2Ti(III)Cl leading to a Ti(IV)-bound β-alkyl radical, (2) hydrogen-atom transfer (HAT) catalyzed by the Co(II) catalyst to form the Ti(IV)-enolate and Co(III)-H intermediate, (3) protonation to give the alcohols, and (4) proton abstraction to form the Co(I) species followed by electron transfer to regenerate the active Co(II) and Ti(III) species. Moreover, bimetallic catalysis and two-state reactivity enable the key rate-determining HAT step. Furthermore, a subtle balance between dispersion-driven bimetallic processes and entropy-driven monometallic processes determines the most favorable pathway, among which the monometallic process is energetically more favorable in all steps except the vital hydrogen-atom transfer step. Our study should provide an in-depth mechanistic understanding of bimetallic catalysis.

Metal-Ligand Proton Tautomerism, Electron Transfer, and C(sp3)-H Activation by a 4-Pyridinyl-Pincer Iridium Hydride Complex

The para-N-pyridyl-based PCP pincer proligand 3,5-bis(di-tert-butylphosphinomethyl)-2,6-dimethylpyridine (pN-tBuPCP-H) was synthesized and metalated to give the iridium complex (pN-tBuPCP)IrHCl (2-H). In marked contrast with its phenyl-based congeners, e.g., (tBuPCP)IrHCl and derivatives, 2-H is highly air-sensitive and reacts with oxidants such as ferrocenium, trityl cation, and benzoquinone. These oxidations ultimately lead to intramolecular activation of a phosphino-t-butyl C(sp3)-H bond and cyclometalation. Considering the greater electronegativity of N than C, 2-H is expected to be less easily oxidized than simple PCP derivatives; cyclic voltammetry and DFT calculations support this expectation. However, 2-H is calculated to undergo metal-ligand-proton tautomerism (MLPT) to give an N-protonated complex that can be described with resonance forms representing a zwitterionic complex (with a negative charge on Ir) and a p-N-pyridylidene (a remote N-heterocyclic carbene) Ir(I) complex. One-electron oxidation of this tautomer is calculated to be dramatically more favorable than direct oxidation of 2-H (ΔΔG° = -31.3 kcal/mol). The resulting Ir(II) oxidation product is easily deprotonated to give metalloradical 2• which is observed by NMR spectroscopy. 2• can be further oxidized to give cationic Ir(III) complex, 2+, which can oxidatively add a phosphino-t-butyl C-H bond and undergo deprotonation to give the observed cyclometalated product. DFT calculations indicate that less sterically hindered analogues of 2+ would preferentially undergo intermolecular addition of C(sp3)-H bonds, for example, of n-alkanes. The resulting iridium alkyl complexes could undergo facile β-H elimination to afford olefin, thereby completing a catalytic cycle for alkane dehydrogenation driven by one-electron oxidation and deprotonation, enabled by MLPT.

Electron-rich benzofulvenes as effective dipolarophiles in copper(i)-catalyzed asymmetric 1,3-dipolar cycloaddition of azomethine ylides

A series of benzofulvenes without any electron-withdrawing substituents were employed as 2π-type dipolarophiles for the first time to participate in Cu(i)-catalyzed asymmetric 1,3-dipolar cycloaddition (1,3-DC) reactions of azomethine ylides. An intrinsic non-benzenoid aromatic characteristic from benzofulvenes serves as a key driving force for activation of the electron-rich benzofulvenes. Utilizing the current methodology, a wide range of multi-substituted chiral spiro-pyrrolidine derivatives containing two contiguous all-carbon quaternary centers were formed in good yield with exclusive chemo-/regioselectivity and high to excellent stereoselectivity. Computational mechanistic studies elucidate the origin of the stereochemical outcome and the chemoselectivity, in which the thermostability of these cycloaddition products is the major factor.

Irreversible inactivation of lactate racemase by sodium borohydride reveals reactivity of the nickel-pincer nucleotide cofactor

The nickel-pincer nucleotide (NPN) cofactor discovered in lactate racemase from Lactiplantibacillus plantarum (LarALp) is essential for the activities of racemases/epimerases in the highly diverse LarA superfamily. Prior mechanistic studies have established a proton-coupled hydride-transfer mechanism for LarALp, but direct evidence showing that hydride attacks the C4 atom in the pyridinium ring of NPN has been lacking. Here, we show that sodium borohydride (NaBH4) irreversibly inactivates LarALp accompanied by a rapid color change of the enzyme. The altered ultraviolet-visible spectra during NaBH4 titration supported hydride transfer to C4 of NPN, and the concomitant Ni loss unraveled by mass spectrometry experiments accounted for the irreversible inactivation. High resolution structures of LarALp revealed a substantially weakened C-Ni bond in the metastable sulfite-NPN adduct where the NPN cofactor is in the reduced state. These findings allowed us to propose a mechanism of LarALp inactivation by NaBH4 that provides key insights into the enzyme-catalyzed reaction and sheds light on the reactivity of small molecule NPN mimetics.

Unveiling the mechanisms and biosynthesis of a novel nickel-pincer enzyme

The nickel-pincer nucleotide (NPN) coenzyme, a substituted pyridinium mononucleotide that tri-coordinates nickel, was first identified covalently attached to a lysine residue in the LarA protein of lactate racemase. Starting from nicotinic acid adenine dinucleotide, LarB carboxylates C5 of the pyridinium ring and hydrolyzes the phosphoanhydride, LarE converts the C3 and C5 carboxylates to thiocarboxylates, and LarC incorporates nickel to form a C-Ni and two S-Ni bonds, during the biosynthesis of this cofactor. LarB uses a novel carboxylation mechanism involving the transient formation of a cysteinyl-pyridinium adduct. Depending on the source of the enzyme, LarEs either catalyze a sacrificial sulfur transfer from a cysteinyl side chain resulting in the formation of dehydroalanine or they utilize a [4Fe-4S] cluster bound by three cysteine residues to accept and transfer a non-core sulfide atom. LarC is a CTP-dependent enzyme that cytidinylylates its substrate, adds nickel, then hydrolyzes the product to release NPN and CMP. Homologs of the four lar genes are widely distributed in microorganisms, with some species containing multiple copies of larA whereas others lack this gene, consistent with the cofactor serving other functions. Several LarA-like proteins were shown to catalyze racemase or epimerase activities using 2-hydroxyacid substrates other than lactic acid. Thus, lactate racemase is the founding member of a large family of NPN-containing enzymes.

New Insights and Predictions into Complex Homogeneous Reactions Enabled by Computational Chemistry in Synergy with Experiments: Isotopes and Mechanisms

Homogeneous catalysis and biocatalysis have been widely applied in synthetic, medicinal, and energy chemistry as well as synthetic biology. Driven by developments of new computational chemistry methods and better computer hardware, computational chemistry has become an essentially indispensable mechanistic "instrument" to help understand structures and decipher reaction mechanisms in catalysis. In addition, synergy between computational and experimental chemistry deepens our mechanistic understanding, which further promotes the rational design of new catalysts. In this Account, we summarize new or deeper mechanistic insights (including isotope, dispersion, and dynamical effects) into several complex homogeneous reactions from our systematic computational studies along with subsequent experimental studies by different groups. Apart from uncovering new mechanisms in some reactions, a few computational predictions (such as excited-state heavy-atom tunneling, steric-controlled enantioswitching, and a new geminal addition mechanism) based on our mechanistic insights were further verified by ensuing experiments.The Zimmerman group developed a photoinduced triplet di-π-methane rearrangement to form cyclopropane derivatives. Recently, our computational study predicted the first excited-state heavy-atom (carbon) quantum tunneling in one triplet di-π-methane rearrangement, in which the reaction rates and ¹²C/¹³C kinetic isotope effects (KIEs) can be enhanced by quantum tunneling at low temperatures. This unprecedented excited-state heavy-atom tunneling in a photoinduced reaction has recently been verified by an experimental ¹²C/¹³C KIE study by the Singleton group. Such combined computational and experimental studies should open up opportunities to discover more rare excited-state heavy-atom tunneling in other photoinduced reactions. In addition, we found unexpectedly large secondary KIE values in the five-coordinate Fe(III)-catalyzed hetero-Diels-Alder pathway, even with substantial C-C bond formation, due to the non-negligible equilibrium isotope effect (EIE) derived from altered metal coordination. Therefore, these KIE values cannot reliably reflect transition-state structures for the five-coordinate metal pathway. Furthermore, our density functional theory (DFT) quasi-classical molecular dynamics (MD) simulations demonstrated that the coordination mode and/or spin state of the iron metal as well as an electric field can affect the dynamics of this reaction (e.g., the dynamically stepwise process, the entrance/exit reaction channels).Moreover, we unveiled a new reaction mechanism to account for the uncommon Ru(II)-catalyzed geminal-addition semihydrogenation and hydroboration of silyl alkynes. Our proposed key gem-Ru(II)-carbene intermediates derived from double migrations on the same alkyne carbon were verified by crossover experiments. Additionally, our DFT MD simulations suggested that the first hydrogen migration transition-state structures may directly and quickly form the key gem-Ru-carbene structures, thereby "bypassing" the second migration step. Furthermore, our extensive study revealed the origin of the enantioselectivity of the Cu(I)-catalyzed 1,3-dipolar cycloaddition of azomethine ylides with β-substituted alkenyl bicyclic heteroarenes enabled by dual coordination of both substrates. Such mechanistic insights promoted our computational predictions of the enantioselectivity reversal for the corresponding monocyclic heteroarene substrates and the regiospecific addition to the less reactive internal C═C bond of one diene substrate. These predictions were proven by our experimental collaborators. Finally, our mechanistic insights into a few other reactions are also presented. Overall, we hope that these interactive computational and experimental studies enrich our mechanistic understanding and aid in reaction development.

Biological pincer complexes

At least two types of pincer complexes are known to exist in biology. A metal-pyrroloquinolone quinone (PQQ) cofactor was first identified in bacterial methanol dehydrogenase, and later also found in selected short-chain alcohol dehydrogenases of other microorganisms. The PQQ-associated metal can be calcium, magnesium, or a rare earth element depending on the enzyme sequence. Synthesis of this organic ligand requires a series of accessory proteins acting on a small peptide, PqqA. Binding of metal to PQQ yields an ONO-type pincer complex. More recently, a nickel-pincer nucleotide (NPN) cofactor was discovered in lactate racemase, LarA. This cofactor derives from nicotinic acid adenine dinucleotide via action of a carboxylase/hydrolase, sulfur transferase, and nickel insertase, resulting in an SCS-type pincer complex. The NPN cofactor likely occurs in selected other racemases and epimerases of bacteria, archaea, and a few eukaryotes.

Unusual KIE and dynamics effects in the Fe-catalyzed hetero-Diels-Alder reaction of unactivated aldehydes and dienes

Hetero-Diels-Alder (HDA) reaction is an important synthetic method for many natural products. An iron(III) catalyst was developed to catalyze the challenging HDA reaction of unactivated aldehydes and dienes with high selectivity. Here we report extensive density-functional theory (DFT) calculations and molecular dynamics simulations that show effects of iron (including its coordinate mode and/or spin state) on the dynamics of this reaction: considerably enhancing dynamically stepwise process, broadening entrance channel and narrowing exit channel from concerted asynchronous transition states. Also, our combined computational and experimental secondary KIE studies reveal unexpectedly large KIE values for the five-coordinate pathway even with considerable C-C bond forming, due to equilibrium isotope effect from the change in the metal coordination. Moreover, steric and electronic effects are computationally shown to dictate the C=O chemoselectivity for an α,β-unsaturated aldehyde, which is verified experimentally. Our mechanistic study may help design homogeneous, heterogeneous and biological catalysts for this challenging reaction.

Functional Models of the Nickel Pincer Nucleotide Cofactor of Lactate Racemase

A novel nickel pincer cofactor was recently discovered in lactate racemase. Reported here are three synthetic nickel pincer complexes that are both structural and functional models of the pincer cofactor in lactate racemase. DFT computations suggest the ipso-carbon atom of the pyridinium pincer ligands act as a hydride acceptor for lactate isomerization, whereas an organometallic pathway involving nickel-mediated β-hydride elimination is less favored.

Theoretical Studies of Nickel-Dependent Enzymes

The advancements of quantum chemical methods and computer power allow detailed mechanistic investigations of metalloenzymes. In particular, both quantum chemical cluster and combined QM/MM approaches have been used, which have been proven to successfully complement experimental studies. This review starts with a brief introduction of nickel-dependent enzymes and then summarizes theoretical studies on the reaction mechanisms of these enzymes, including NiFe hydrogenase, methyl-coenzyme M reductase, nickel CO dehydrogenase, acetyl CoA synthase, acireductone dioxygenase, quercetin 2,4-dioxygenase, urease, lactate racemase, and superoxide dismutase.

Nickel-pincer nucleotide cofactor

A novel organometallic cofactor, nickel pyridinium-3,5-dithiocarboxylic acid mononucleotide, was recently discovered in lactate racemase (LarA) of Lactobacillus plantarum. This review summarizes the substantial progress made in uncovering the function of this cofactor as a transient hydride acceptor in the LarA mechanism. The latest developments related to cofactor biosynthesis reveal insights into a pathway in which LarB serves as a nicotinic acid adenine dinucleotide hydrolase/carboxylase, LarE acts as a sacrificial sulfur transferase, and LarC functions as a nickel insertase, forming the nickel-pincer nucleotide cofactor that becomes covalently tethered to LarA in some bacteria. Bioinformatic studies reveal a widespread occurrence of larA, larB, larC, and larE orthologs in microorganisms, and additional roles for the cofactor are considered.

Lactate Racemase Nickel-Pincer Cofactor Operates by a Proton-Coupled Hydride Transfer Mechanism

Lactate racemase (LarA) of Lactobacillus plantarum contains a novel organometallic cofactor with nickel coordinated to a covalently tethered pincer ligand, pyridinium-3-thioamide-5-thiocarboxylic acid mononucleotide, but its function in the enzyme mechanism has not been elucidated. This study presents direct evidence that the nickel-pincer cofactor facilitates a proton-coupled hydride transfer (PCHT) mechanism during LarA-catalyzed lactate racemization. No signal was detected by electron paramagnetic resonance spectroscopy for LarA in the absence or presence of substrate, consistent with a +2 metal oxidation state and inconsistent with a previously proposed proton-coupled electron transfer mechanism. Pyruvate, the predicted intermediate for a PCHT mechanism, was observed in quenched solutions of LarA. A normal substrate kinetic isotope effect ( k_H/ k_D of 3.11 ± 0.17) was established using 2-α-²H-lactate, further supporting a PCHT mechanism. UV-visible spectroscopy revealed a lactate-induced perturbation of the cofactor spectrum, notably increasing the absorbance at 340 nm, and demonstrated an interaction of the cofactor with the inhibitor sulfite. A crystal structure of LarA provided greater resolution (2.4 Å) than previously reported and revealed sulfite binding to the pyridinium C4 atom of the reduced pincer cofactor, mimicking hydride reduction during a PCHT catalytic cycle. Finally, computational modeling supports hydride transfer to the cofactor at the C4 position or to the nickel atom, but with formation of a nickel-hydride species requiring dissociation of the His200 metal ligand. In aggregate, these studies provide compelling evidence that the nickel-pincer cofactor acts by a PCHT mechanism.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Structural insights into the catalytic mechanism of a sacrificial sulfur insertase of the N-type ATP pyrophosphatase family, LarE

The lar operon in Lactobacillus plantarum encodes five Lar proteins (LarA/B/C/D/E) that collaboratively synthesize and incorporate a niacin-derived Ni-containing cofactor into LarA, an Ni-dependent lactate racemase. Previous studies have established that two molecules of LarE catalyze successive thiolation reactions by donating the sulfur atom of their exclusive cysteine residues to the substrate. However, the catalytic mechanism of this very unusual sulfur-sacrificing reaction remains elusive. In this work, we present the crystal structures of LarE in ligand-free and several ligand-bound forms, demonstrating that LarE is a member of the N-type ATP pyrophosphatase (PPase) family with a conserved N-terminal ATP PPase domain and a unique C-terminal domain harboring the putative catalytic site. Structural analysis, combined with structure-guided mutagenesis, leads us to propose a catalytic mechanism that establishes LarE as a paradigm for sulfur transfer through sacrificing its catalytic cysteine residue.

From NAD+ to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase

Lactate racemase (LarA), a new nickel enzyme discovered recently, catalyzes the racemization between d- and l-lactates with a novel nickel pincer cofactor (Ni-PTTMN) derived from nicotinic acid. In this study, by using DFT and a 200-atom active-site model, LarA is revealed to employ a modified proton-coupled hydride-transfer mechanism in which a hydride is transferred to a cofactor pyridine carbon from the substrate α-carbon along with proton transfer from the substrate hydroxy group to a histidine, and then moved back from the opposite side. Tyr294 and Lys298 provide significant acceleration effects by orientating substrates and stabilizing the negative charge developing at the substrate hydroxy oxygen. The barrier was determined to be 12.0 kcal mol^-1 , which reveals enhanced racemase activity relative to the LarA reaction using NAD⁺ -like cofactors. Compared with NAD⁺ , Ni-PTTMN has a stronger hydride-addition reactivity in moderate and high environmental polarity and may fit perfectly the moderately polar active site of LarA.