Computational Study of Engineered Cytochrome P450-Catalyzed C–H Amination: The Origin of the Regio- and Stereoselectivity

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.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

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- Collapse

Key Words

• amination
• biocatalysis
• density functional calculations
• heme proteins
• reaction mechanism

Collapse

MESH Headings

• Hemeproteins
• Ligands
• Amination
• Heme/chemistry
• Amines

Collapse

Grants

• CHE-1300912 and CHE-2054897 National Science Foundation
• National Science Foundation

Collapse

Affiliation(s)

Yang Wei Department of Chemistry and Chemical BiologyStevens Institute of Technology1 Castle Point on HudsonHobokenNJ 07030USA Department of Chemistry and BiochemistryLoyola University Chicago1032 W Sheridan RdChicagoIL 60660USA
Melissa Conklin Department of Chemistry and Chemical BiologyStevens Institute of Technology1 Castle Point on HudsonHobokenNJ 07030USA
Yong Zhang Department of Chemistry and Chemical BiologyStevens Institute of Technology1 Castle Point on HudsonHobokenNJ 07030USA

Collapse

MD simulations and QM/MM calculations reveal the key mechanistic elements which are responsible for the efficient C-H amination reaction performed by a bioengineered P450 enzyme

An enzyme which is capable of catalyzing C–H amination reactions is considered to be a dream tool for chemists due to its pharmaceutical potential and greener approach. Recently, the Arnold group achieved this feat using an engineered CYP411 enzyme, which further undergoes a random directed evolution which increases its efficiency and selectivity. The present study provides mechanistic insight and the root cause of the success of these mutations to enhance the reactivity and selectivity of the mutant enzyme. This is achieved by means of comprehensive MD simulations and hybrid QM/MM calculations. The study shows that the efficient C–H amination by the engineered CYP411 is a combined outcome of electronic and steric effects. The mutation of the axial cysteine ligand to serine relays electron density to the Fe ion in the heme, and thereby enhances the bonding capability of the heme-iron to the nitrogen atom of the tosyl azide. In comparison, the native cysteine-ligated P450 cannot bind the tosyl azide. Additionally, the A78V and A82L mutations in P411 provide ‘bulk’ to the active site which increases the enantioselectivity via a steric effect. At the same time, the QM/MM calculations elucidate the C–H amination by the iron nitrenoid, revealing a mechanism analogous to Compound I in the native C–H hydroxylation by P450.

Computer simulation method reveals the mechanism of C–H amination reaction due to a single site mutation.

Nitrene transfers mediated by natural and artificial iron enzymes

Amines are ubiquitous in biology and pharmacy. As a consequence, introducing N functionalities in organic molecules is attracting strong continuous interest. The past decade has witnessed the emergence of very efficient and selective catalytic systems achieving this goal thanks to engineered hemoproteins. In this review, we examine how these enzymes have been engineered focusing rather on the rationale behind it than the methodology employed. These studies are put in perspective with respect to in vitro and in vivo nitrene transfer processes performed by cytochromes P450. An emphasis is put on mechanistic aspects which are confronted to current molecular knowledge of these reactions. Forthcoming developments are delineated.

Theoretical studies on iron-catalyzed azaindoline formation: mechanism and site-selectivity

The mechanism and site-selectivity for Fe-catalyzed azaindoline formation from 1,2,3,4-tetrazole were examined computationally. The H-atom abstraction/radical rebound stepwise mechanism is proposed. The aliphatic H-atom abstraction (HAA) vs. electrophilic aromatic substitution (EAS) steps are responsible for the sp³vs. sp² C-H amination site-selectivity and a larger steric congestion disfavors sp² EAS, thus resulting in Fe-catalyzed site-selectivity toward sp³ C-H amination.

Computational insights into Ir(iii)-catalyzed allylic C-H amination of terminal alkenes: mechanism, regioselectivity, and catalytic activity

Computational studies on Ir(iii)-catalyzed intermolecular branch-selective allylic C–H amination of terminal olefins with methyl dioxazolone have been carried out to investigate the mechanism, including the origins of regioselectivity and catalytic activity difference. The result suggests that the reaction proceeds through generation of active species, alkene coordination, allylic C–H activation, decarboxylation, migratory insertion, and protodemetalation. The presence of AgNTf₂ could thermodynamically promote the formation of catalytically active species [Cp*Ir(OAc)]⁺. Both the weaker Ir–C(internal) bond and the closer interatomic distance of N⋯C(internal) in the key allyl-Ir(v)-nitrenoid intermediate make the migratory insertion into Ir–C(internal) bond easier than into the Ir–C(terminal) bond, leading to branch-selective allylic C–H amidation. The high energy barrier for allylic C–H activation in the Co system could account for the observed sluggishness, which is mainly ascribed to the weaker coordination capacity of alkenes to the triplet Cp*Co(OAc)⁺ and the deficient metal⋯H interaction to assist hydrogen transfer.

DFT studies on Ir(iii)-catalyzed branch-selective allylic C–H amination of terminal olefins with methyl dioxazolone have been carried out to investigate the mechanism, including the origins of regioselectivity and catalytic activity difference.

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

Engineered heme enzymes such as myoglobin and cytochrome P450s metalloproteins are gaining widespread importance due to their efficiency in catalyzing non-natural reactions. In a recent strategy, the naturally occurring Fe metal in the heme unit was replaced with non-native metals such as Ir, Rh, Co, Cu, etc., and axial ligands to generate artificial metalloenzymes. Determining the best metal-ligand for a chemical transformation is not a trivial task. Here we demonstrate how computational approaches can be used in deciding the best metal-ligand combination which would be highly beneficial in designing new enzymes as well as small molecule catalysts. We have used Density Functional Theory (DFT) to shed light on the enhanced reactivity of an Ir system with varying axial ligands. We look at the insertion of a carbene group generated from diazo precursors via N2 extrusion into a C-H bond. For both Ir(Me) and Fe systems, the first step, i.e., N2 extrusion is the rate determining step. Strikingly, neither the better ligand overlap with 5d orbitals on Ir nor the electrophilicity on the carbene centre play a significant role. A comparison of Fe and Ir systems reveals that a lower distortion in the Ir(Me)-porphyrin on moving from the reactant to the transition state renders it catalytically more active. We notice that for both metal porphyrins, the free energy barriers are affected by axial ligand substitution. Further, for Fe porphyrin, the axial ligand also changes the preferred spin state. We show that for the carbene insertion into the C-H bond, Fe porphyrin systems undergo a stepwise HAT (hydrogen atom transfer) instead of a concerted hydride transfer process. Importantly, we find that the substitution of the axial Me ligand on Ir to imidazole or chloride, or without an axial substitution changes the rate determining step of the reaction. Therefore, an optimum ligand that can balance the barriers for both steps of the catalytic cycle is essential. We subsequently used the QM cluster approach to delineate the protein environment's role and mutations in improving the catalytic activity of the Ir(Me) system.

Unravelling nitrene chemistry from acyclic precursors: recent advances and challenges

Recent advances in nitrene chemistry from acyclic precursors are reviewed in this paper.

The molecular mechanism of P450-catalyzed amination of the pyrrolidine derivative of lidocaine: insights from multiscale simulations

Nitrogen heterocycles are key and prevalent motifs in drugs. Evolved variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium employ high-valent oxo-iron(iv) species to catalyze the synthesis of imidazolidine-4-ones via an intramolecular C–H amination. Herein, we use multi-scale simulations, including classical molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical (QM/MM) calculations and QM calculations, to reveal the molecular mechanism of the intramolecular C–H amination of the pyrrolidine derivative of lidocaine bearing cyclic amino moieties catalyzed by the variant RP/FV/EV of P450_BM3, which bears five mutations compared to wild type. Our calculations show that overall catalysis includes both the enzymatic transformation in P450 and non-enzymatic transformation in water solution. The enzymatic transformation involves the exclusive hydroxylation of the C–H bond of the pyrrolidine derivative of lidocaine, leading to the hydroxylated intermediate, during which the substrate radical would be bypassed. The following dehydration and C–N coupling reactions are found to be much favored in aqueous situation compared to that in the non-polar protein environment. The present findings expand our understanding of the P450-catalyzed C(sp₃)–H amination reaction.

Nitrogen heterocycles are key and prevalent motifs in drugs.

Mechanism-Guided Design and Discovery of Efficient Cytochrome P450-Derived C-H Amination Biocatalysts

Cytochromes P450 have been recently identified as a promising class of biocatalysts for mediating C-H aminations via nitrene transfer, a valuable transformation for forging new C-N bonds. The catalytic efficiency of P450s in these non-native transformations is however significantly inferior to that exhibited by these enzymes in their native monooxygenase function. Using a mechanism-guided strategy, we report here the rational design of a series of P450_BM3-based variants with dramatically enhanced C-H amination activity acquired through disruption of the native proton relay network and other highly conserved structural elements within this class of enzymes. This approach further guided the identification of XplA and BezE, two "atypical" natural P450s implicated in the degradation of a man-made explosive and in benzastatins biosynthesis, respectively, as very efficient C-H aminases. Both XplA and BezE could be engineered to further improve their C-H amination reactivity, which demonstrates their evolvability for abiological reactions. These engineered and natural P450 catalysts can promote the intramolecular C-H amination of arylsulfonyl azides with over 10 000-14 000 catalytic turnovers, ranking among the most efficient nitrene transfer biocatalysts reported to date. Mechanistic and structure-reactivity studies provide insights into the origin of the C-H amination reactivity enhancement and highlight the divergent structural requirements inherent to supporting C-H amination versus C-H monooxygenation reactivity within this class of enzymes. Overall, this work provides new promising scaffolds for the development of nitrene transferases and demonstrates the value of mechanism-driven rational design as a strategy for improving the catalytic efficiency of metalloenzymes in the context of abiological transformations.

Cobalt-catalysed unactivated C(sp3)–H amination: two-state reactivity and multi-reference electronic character

A remarkable two-state reactivity scenario and an unusual multi-reference character have been computationally found in Co-catalysed C(sp³)–H amination. In addition, the investigation on the additive, aminating reagent, metal center, and auxiliary ligand provides implications for development of new catalytic C–H functionalization systems.

Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2

Vitamin K₁ (VK1) plays an important role in the modulation of bleeding disorders. It has been reported that ω-hydroxylation on the VK1 aliphatic chain is catalyzed by cytochrome P450 4F2 (CYP4F2), an enzyme responsible for the metabolism of eicosanoids. However, the mechanism of VK1 ω-hydroxylation by CYP4F2 has not been disclosed. In this study, we employed a combination of quantum mechanism (QM) calculations, homology modeling, molecular docking, molecular dynamics (MD) simulations, and combined quantum mechanism/molecular mechanism (QM/MM) calculations to investigate the metabolism profile of VK1 ω-hydroxylation. QM calculations based on the truncated VK1 model show that the energy barrier for ω-hydroxylation is about 6-25 kJ/mol higher than those at other potential sites of metabolism. However, results from the MD simulations indicate that hydroxylation at the ω-site is more favorable than at the other potential sites, which is in accordance with the experimental observation. The evaluation of MD simulations was further endorsed by the QM/MM calculation results. Our studies thus suggest that the active site residues of CYP4F2 play a determinant role in the ω-hydroxylation. Our results provide structural insights into the mechanism of VK1 ω-hydroxylation by CYP4F2 at the atomistic level and are helpful not only for characterizing the CYP4F2 functions but also for looking into the ω-hydroxylation mediated by other CYP4 enzymes.