A computational approach to understand the role of metals and axial ligands in artificial heme enzyme catalyzed C-H insertion

Elucidating the Enzymatic Mechanism of Dihydrocoumarin Degradation: Insight into the Functional Evolution of Methyl-Parathion Hydrolase from QM/MM and MM MD Simulations

Methyl-parathion hydrolase (MPH), which evolved from dihydrocoumarin hydrolase, offers one of the most efficient enzymes for the hydrolysis of methyl-parathion. Interestingly, the substrate preference of MPH shifts from the methyl-parathion to the lactone dihydrocoumarin (DHC) after its mutation of five specific residues (R72L, L273F, L258H, T271I, and S193Δ, m5-MPH). Here, extensive QM/MM calculations and MM MD simulations have been used to delve into the structure-function relationship of MPH enzymes and plausible mechanisms for the chemical and nonchemical steps, including the transportation and binding of the substrate DHC to the active site, the hydrolysis reaction, and the product release. The results reveal that the five mutations remodel the active pocket and reposition DHC within the active site, leading to stronger enzyme-substrate interactions. The MM/GBSA-estimated binding free energies are about -20.7 kcal/mol for m5-MPH and -17.1 kcal/mol for wild-type MPH. Furthermore, this conformational adjustment of the protein may facilitate the chemical step of DHC hydrolysis and the product release, although there is a certain influence on the substrate transport. The hydrolytic reaction begins with the nucleophilic attack of the bridging OH- with the energy barriers of 22.0 and 18.0 kcal/mol for the wild-type and m5-MPH enzymes, respectively, which is rate-determining for the entire process. Unraveling these mechanistic intricacies may help in the understanding of the natural evolution of enzymes for diverse substrates and establish the enzyme structure-function relationship.

Elucidation of the catalytic mechanism of a single-metal dependent homing endonuclease using QM and QM/MM approaches: the case study of I-PpoI

Homing endonucleases (HEs) are highly specific DNA cleaving enzymes, with I-PpoI having been suggested to use a single metal to accelerate phosphodiester bond cleavage. Although an I-PpoI mechanism has been proposed based on experimental structural data, no consensus has been reached regarding the roles of the metal or key active site amino acids. This study uses QM cluster and QM/MM calculations to provide atomic-level details of the I-PpoI catalytic mechanism. Minimal QM cluster and large-scale QM/MM models demonstrate that the experimentally-proposed pathway involving direct Mg2+ coordination to the substrate coupled with leaving group protonation through a metal-activated water is not feasible due to an inconducive I-PpoI active site alignment. Despite QM cluster models of varying size uncovering a pathway involving leaving group protonation by a metal-activated water, indirect (water-mediated) metal coordination to the substrate is required to afford this pathway, which renders this mechanism energetically infeasible. Instead, QM cluster models reveal that the preferred pathway involves direct Mg2+-O3' coordination to stabilize the charged substrate and assist leaving group departure, while H98 activates the water nucleophile. These calculations also underscore that both catalytic residues that directly interact with the substrate and secondary amino acids that position or stabilize these residues are required for efficient catalysis. QM/MM calculations on the solvated enzyme-DNA complex verify the preferred mechanism, which is fully consistent with experimental kinetic, structural, and mutational data. The fundamental understanding of the I-PpoI mechanism of action, gained from the present work can be used to further explore potential uses of this enzyme in biotechnology and medicine, and direct future computational investigations of other members of the understudied HE family.

Is Enol Always the Culprit? The Curious Case of High Enantioselectivity in a Chiral Rh(II) Complex Catalyzed Carbene Insertion Reaction

The mechanism of Rh2 (S-NTTL)4 catalyzed carbene insertion into C(3)-H of indole is investigated using DFT methods. Since the commonly accepted enol mechanism cannot account for enantioinduction, a concerted oxocarbenium pathway was proposed in an earlier work using a model catalyst. However, after considering the full catalytic system, this study finds that akin to other reactions, here, too, the enol pathway is of lower energy, which now naturally raises a conundrum regarding the mode of chiral induction. Herein, a new water promoted mechanistic pathway involving a metal-associated enol intermediate hydrogen bonding and stereochemical model are proposed to solve this puzzle. It is shown how the catalyst bowl-shaped structure along with substrate-catalyst binding is crucial for achieving high levels of enantioselectivity. A stereodetermining water-assisted proton transfer is proposed and confirmed through deuterium-labeling experiments. The water molecules are held together by H-bonding interactions with the carboxylate ligands that is reminiscent of enzyme catalysis. Although several previous studies have aimed at understanding the mechanism of metal catalyzed carbene insertion reactions, the origin of high stereoinduction especially with chiral metal complexes remains unclear, and till date there is no transition state model that can explain the high enantioselectivity with such chiral Rh complexes. The metal-associated enol pathway is currently underrepresented in catalytic cycles and may play a crucial role in catalyst design. Since the enol pathway is commonly adopted in other metal-catalyzed X-H insertion reactions involving a diazoester, the presented results are not specific to the current reaction. Therefore, this study could provide the direction for achieving high levels of enantioselectivity which is otherwise difficult to achieve with a single metal catalyst.

Computational Mechanistic Investigations of Biocatalytic Nitrenoid C-H Functionalizations via Engineered Heme Proteins

Engineered heme proteins were developed to possess numerous excellent biocatalytic nitrenoid C-H functionalizations. Computational approaches such as density functional theory (DFT), hybrid quantum mechanics/molecular mechanics (QM/MM), and molecular dynamics (MD) calculations were employed to help understand some important mechanistic aspects of these heme nitrene transfer reactions. This review summarizes advances of computational reaction pathway results of these biocatalytic intramolecular and intermolecular C-H aminations/amidations, focusing on mechanistic origins of reactivity, regioselectivity, enantioselectivity, diastereoselectivity as well as effects of substrate substituent, axial ligand, metal center, and protein environment. Some important common and distinctive mechanistic features of these reactions were also described with brief outlook of future development.

Role of mutations in a chemoenzymatic enantiodivergent C(sp3)-H insertion: exploring the mechanism and origin of stereoselectivity

New-to-nature enzymes have emerged as powerful catalysts in recent years for streamlining various stereoselective organic transformations. While synthetic strategies employing engineered enzymes have witnessed proliferating success, there is limited clarity on the mechanistic front and more so when considering molecular-level insights into the role of selected mutations, dramatically escalating catalytic competency and selectivity. We have investigated the mechanism and correlation between mutations and exquisite stereoselectivity of a lactone carbene insertion into the C(sp3)-H bond of substituted aniline, catalyzed by two mutants of a cytochrome P450 variant, "P411" (engineered through directed evolution) in which the axial cysteine has been mutated to serine, utilizing various computational tools. The pivotal role of S264 and L/R328 mutations in the active site has been delineated computationally using two cluster models, thus rationalizing the enantiodivergence. This report provides much-needed insights into the origin of enantiodivergence, furnishing a mechanistic framework for understanding the anchoring effects of H-bond donor residues with the lactone ring. This study is expected to have important implications in the rational design of stereodivergent enzymes and toward successful in silico enzyme designing.

Machine Learning and Quantum Calculation for Predicting Yield in Cu-Catalyzed P-H Reactions

The paper discussed the use of machine learning (ML) and quantum chemistry calculations to predict the transition state and yield of copper-catalyzed P-H insertion reactions. By analyzing a dataset of 120 experimental data points, the transition state was determined using density functional theory (DFT). ML algorithms were then applied to analyze 16 descriptors derived from the quantum chemical transition state to predict the product yield. Among the algorithms studied, the Support Vector Machine (SVM) achieved the highest prediction accuracy of 97%, with over 80% correlation in Leave-One-Out Cross-Validation (LOOCV). Sensitivity analysis was performed on each descriptor, and a comprehensive investigation of the reaction mechanism was conducted to better understand the transition state characteristics. Finally, the ML model was used to predict reaction plans for experimental design, demonstrating strong predictive performance in subsequent experimental validation.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

YfeX - A New Platform for Carbene Transferase Development with High Intrinsic Reactivity

Carbene transfer biocatalysis has evolved from basic science to an area with vast potential for the development of new industrial processes. In this study, we show that YfeX, naturally a peroxidase, has great potential for the development of new carbene transferases, due to its high intrinsic reactivity, especially for the N-H insertion reaction of aromatic and aliphatic primary and secondary amines. YfeX shows high stability against organic solvents (methanol and DMSO), greatly improving turnover of hydrophobic substrates. Interestingly, in styrene cyclopropanation, WT YfeX naturally shows high enantioselectivity, generating the trans product with 87 % selectivity for the (R,R) enantiomer. WT YfeX also catalyzes the Si-H insertion efficiently. Steric effects in the active site were further explored using the R232A variant. Quantum Mechanics/Molecular Mechanics (QM/MM) calculations reveal details on the mechanism of Si-H insertion. YfeX, and potentially other peroxidases, are exciting new targets for the development of improved carbene transferases.

Biocatalytic Intramolecular C-H aminations via Engineered Heme Proteins: Full Reaction Pathways and Axial Ligand Effects

Engineered heme protein biocatalysts provide an efficient and sustainable approach to develop amine-containing compounds through C-H amination. A quantum chemical study to reveal the complete heme catalyzed intramolecular C-H amination pathway and protein axial ligand effect was reported, using reactions of an experimentally used arylsulfonylazide with hemes containing L=none, SH^- , MeO^- , and MeOH to simulate no axial ligand, negatively charged Cys and Ser ligands, and a neutral ligand for comparison. Nitrene formation was found as the overall rate-determining step (RDS) and the catalyst with Ser ligand has the best reactivity, consistent with experimental reports. Both RDS and non-RDS (nitrene transfer) transition states follow the barrier trend of MeO^- <SH^- <MeOH<None due to the charge donation capability of the axial ligand to influence the key charge transfer process as the electronic driving forces. Results also provide new ideas for future biocatalyst design with enhanced reactivities.

Does an Enol Pathway Preclude High Stereoselectivity in Iron-Catalyzed Indole C-H Functionalization via Carbene Insertion?

C-H functionalization of indoles via Fe carbenoids presents an attractive strategy to obtain biologically important structural motifs. However, obtaining good stereoselectivity with Fe has been a significant challenge. It is unclear whether the low selectivity is due to a radical pathway or an ionic mechanism involving metal-free species. We therefore present a density functional theory (DFT) study of indole alkylation with diazoacetates catalyzed by Fe(ClO₄)TMEDA/spirobisoxazoline and myoglobin. We explore three mechanistic pathways: nucleophilic, radical, and oxocarbenium routes. The nucleophilic pathway is the most feasible with the formation of an enol species that tautomerizes to furnish the alkylated indole. While this mechanism is routinely proposed, the stereochemical model has been conspicuously absent until now. We show that the conventionally invoked enol pathway is not responsible for the low enantiomeric excess. The enol intermediate can stay coordinated to the catalyst via different binding sites placing the enol in proximity to the chiral environment and affecting the stereoselective proton transfer. Both the binding strength and the chiral environment are crucial for obtaining high selectivity. Our study provides the much needed insights for the modest-low selectivities of Fe systems and could help in expediting the discovery of an efficient catalytic system. These mechanistic underpinnings could also be applicable to other metal (Rh, Pd, Cu, etc.)-catalyzed X-H insertion reactions.

Recent advances in transition-metal-catalyzed carbene insertion to C-H bonds

C-H functionalization has been emerging as a powerful method to establish carbon-carbon and carbon-heteroatom bonds. Many efforts have been devoted to transition-metal-catalyzed direct transformations of C-H bonds. Metal carbenes generated in situ from transition-metal compounds and diazo or its equivalents are usually applied as the transient reactive intermediates to furnish a catalytic cycle for new C-C and C-X bond formation. Using this strategy compounds from unactivated simple alkanes to complex molecules can be further functionalized or transformed to multi-functionalized compounds. In this area, transition-metal-catalyzed carbene insertion to C-H bonds has been paid continuous attention. Diverse catalyst design strategies, synthetic methods, and potential applications have been developed. This critical review will summarize the advance in transition-metal-catalyzed carbene insertion to C-H bonds dated up to July 2021, by the categories of C-H bonds from aliphatic C(sp³)-H, aryl (aromatic) C(sp²)-H, heteroaryl (heteroaromatic) C(sp²)-H bonds, alkenyl C(sp²)-H, and alkynyl C(sp)-H, as well as asymmetric carbene insertion to C-H bonds, and more coverage will be given to the recent work. Due to the rapid development of the C-H functionalization area, future directions in this topic are also discussed. This review will give the authors an overview of carbene insertion chemistry in C-H functionalization with focus on the catalytic systems and synthetic applications in C-C bond formation.

DFT study on Ir-quinoid catalyzed C-H functionalization: new radical reactivity or direct carbene transfer?

DFT methods are used to probe the mechanism of a newly developed Ir-quinoid catalyzed C(sp³)-H functionalization of 1,4 dienes. The lowest energy pathway proceeds via an old-school concerted C-H insertion as opposed to a unique hydrogen atom transfer process proposed previously. The concertedness of the reaction shows an intriguing dependence on sterics of the diene leading to either inserted or dehydrogenated products. We use these new insights to tune the axial ligand, and design a more efficient catalyst.