Reduction of α,β-Unsaturated Ketones by Old Yellow Enzymes: Mechanistic Insights from Quantum Mechanics/Molecular Mechanics Calculations

Unmasking the reverse catalytic activity of 'ene'-reductases for asymmetric carbonyl desaturation

Carbonyl desaturation is a fundamental reaction widely practised in organic synthesis. While numerous methods have been developed to expand the scope of this important transformation, most of them necessitate multi-step protocols or suffer from the use of high loadings of metal or strong oxidizing conditions. Moreover, approaches that can achieve precise stereochemical control of the desaturation process are extremely rare. Here we report a biocatalytic platform for desymmetrizing desaturation of cyclohexanones to generate diverse cyclohexenones bearing a remote quaternary stereogenic centre, by reengineering 'ene'-reductases to efficiently mediate dehydrogenation, the reverse process of their native activity. This 'ene'-reductase-based desaturation system operates under mild conditions with air as the terminal oxidant, tolerates oxidation-sensitive or metal-incompatible functional groups and, more importantly, exhibits unparalleled stereoselectivity compared with those achieved with small-molecule catalysts. Mechanistic investigations suggest that the reaction proceeded through α-deprotonation followed by a rate-determining β-hydride transfer.

Photobiocatalytic Platform for the Efficient Enantio-Divergent Synthesis of β-Fluoromethylated Ketones

β-Fluoromethyl (CH2F, CHF2, and CF3)-substituted chiral ketones are essential moieties and are vital building blocks in pharmaceutical and agrochemistry. However, general and convenient methods for enantio-diverse access to diverse β-fluoromethylated ketones are lacking, hindering the further development of these functional moieties. In this study, we developed an ene-reductase-based photobiocatalytic platform for efficient synthesis of enantio-divergent β-fluoromethylated chiral ketones. Our method highlights substrate-type diversity, excellent enantioselectivity, enzymatic enantio-divergent synthesis, as well as a dicyanopyrazine (DPZ)-type photosensitizer for biocompatible olefin E/Z isomerization in enzymatic stereoconvergent olefin asymmetric reduction, thereby providing a general photobiocatalytic solution to diverse β-fluoromethylated chiral ketones.

Applications of Ene-Reductases in the Synthesis of Flavors and Fragrances

Flavors and fragrances (F&F) are interesting organic compounds in chemistry. These compounds are widely used in the food, cosmetic, and medical industries. Enzymatic synthesis exhibits several advantages over natural extraction and chemical preparation, including a high yield, stable quality, mildness, and environmental friendliness. To date, many oxidoreductases and hydrolases have been used to biosynthesize F&F. Ene-reductases (ERs) are a class of biocatalysts that can catalyze the asymmetric reduction of α,β-unsaturated compounds and offer superior specificity and selectivity; therefore, ERs have been increasingly considered an ideal alternative to their chemical counterparts. This review summarizes the research progress on the use of ERs in F&F synthesis over the past 20 years, including the achievements of various scholars, the differences and similarities among the findings, and the discussions of future research trends related to ERs. We hope this review can inspire researchers to promote the development of biotechnology in the F&F industry.

Strategic enzymatic enantioselective desymmetrization of prochiral cyclohexa-2,5-dienones

Asymmetric desymmetrization through the selective reduction of one double bond of prochiral 2,5-cyclohexadienones is highly challenging. A novel method has been developed for synthesizing chiral cyclohexenones by employing an ene-reductase (Bacillus subtilis YqjM) enzyme that belongs to the OYE family. Our strategy demonstrates high substrate scope and enantioselectivity towards substrates containing all-carbon as well as heteroatom (O, N)-containing quaternary centers. The mechanistic studies (kH/D = ∼1.8) indicate that hydride transfer is probably the rate-limiting step. Mutation of several active site residues did not affect the stereochemical outcomes. This work provides a convenient way of synthesizing various enantioselective γ,γ-disubstituted cyclohexanones using enzymes.

Deciphering the Unconventional Reduction of C=N Bonds by Old Yellow Enzymes Using QM/MM

The reduction of C=X (X = N, O) bonds is a cornerstone in both synthetic organic chemistry and biocatalysis. Conventional reduction mechanisms usually involve a hydride ion targeting the less electronegative carbon atom. In a departure from this paradigm, our investigation into Old Yellow Enzymes (OYEs) reveals a mechanism involving transfer of hydride to the formally more electronegative nitrogen atom within a C=N bond. Beyond their known ability to reduce electronically activated C=C double bonds, e.g., in α, β-unsaturated ketones, these enzymes have recently been shown to reduce α-oximo-β-ketoesters to the corresponding amines. It has been proposed that this transformation involves two successive reduction steps and proceeds via imine intermediates formed by the reductive dehydration of the oxime moieties. We employ advanced quantum mechanics/molecular mechanics (QM/MM) simulations, enriched by a two-tiered approach incorporating QM/MM (UB3LYP-6-31G*/OPLS2005) geometry optimization, QM/MM (B3LYP-6-31G*/amberff19sb) steered molecular dynamics simulations, and detailed natural-bond-orbital analyses to decipher the unconventional hydride transfer to nitrogen in both reduction steps and to delineate the role of active site residues as well as of substituents present in the substrates. Our computational results confirm the proposed mechanism and agree well with experimental mutagenesis and enzyme kinetics data. According to our model, the catalysis of OYE involves hydride transfer from the flavin cofactor to the nitrogen atom in oximoketoesters as well as iminoketoesters followed by protonation at the adjacent oxygen or carbon atoms by conserved tyrosine residues and active site water molecules. Two histidine residues play a key role in the polarization and activation of the C=N bond, and conformational changes of the substrate observed along the reaction coordinate underline the crucial importance of dynamic electron delocalization for efficient catalysis.

Asymmetric α-benzylation of cyclic ketones enabled by concurrent chemical aldol condensation and biocatalytic reduction

Chemoenzymatic cascade catalysis has emerged as a revolutionary tool for streamlining traditional retrosynthetic disconnections, creating new possibilities for the asymmetric synthesis of valuable chiral compounds. Here we construct a one-pot concurrent chemoenzymatic cascade by integrating organobismuth-catalyzed aldol condensation with ene-reductase (ER)-catalyzed enantioselective reduction, enabling the formal asymmetric α-benzylation of cyclic ketones. To achieve this, we develop a pair of enantiocomplementary ERs capable of reducing α-arylidene cyclic ketones, lactams, and lactones. Our engineered mutants exhibit significantly higher activity, up to 37-fold, and broader substrate specificity compared to the parent enzyme. The key to success is due to the well-tuned hydride attack distance/angle and, more importantly, to the synergistic proton-delivery triade of Tyr28-Tyr69-Tyr169. Molecular docking and density functional theory (DFT) studies provide important insights into the bioreduction mechanisms. Furthermore, we demonstrate the synthetic utility of the best mutants in the asymmetric synthesis of several key chiral synthons.

Ene Reductase Enabled Intramolecular β-C-H Functionalization of Substituted Cyclohexanones for Efficient Synthesis of Bridged Bicyclic Nitrogen Scaffolds

Herein we report that ene reductases (EREDs) can facilitate an unprecedented intramolecular β-C-H functionalization reaction for the synthesis of bridged bicyclic nitrogen heterocycles containing the 6-azabicyclo[3.2.1]octane scaffold. To streamline the synthesis of these privileged motifs, we developed a gram-scale one-pot chemoenzymatic cascade by combining iridium photocatalysis with EREDs, using readily available N-phenylglycines and cyclohexenones that can be obtained from biomass. Further derivatization using enzymatic or chemical methods can convert 6-azabicyclo[3.2.1]octan-3-one into 6-azabicyclo[3.2.1]octan-3α-ols, which can be potentially utilized for the synthesis of azaprophen and its analogues for drug discovery. Mechanistic studies revealed the reaction requires oxygen, presumably to produce oxidized flavin, which can selectively dehydrogenate the 3-substituted cyclohexanone derivatives to form the α,β-unsaturated ketone, which subsequently undergoes spontaneous intramolecular aza-Michael addition under basic conditions.

Mechanistic Insights into the Ene-Reductase-Catalyzed Promiscuous Reduction of Oximes to Amines

The biocatalytic reduction of the oxime moiety to the corresponding amine group has only recently been found to be a promiscuous activity of ene-reductases transforming α-oximo β-keto esters. However, the reaction pathway of this two-step reduction remained elusive. By studying the crystal structures of enzyme oxime complexes, analyzing molecular dynamics simulations, and investigating biocatalytic cascades and possible intermediates, we obtained evidence that the reaction proceeds via an imine intermediate and not via the hydroxylamine intermediate. The imine is reduced further by the ene-reductase to the amine product. Remarkably, a non-canonical tyrosine residue was found to contribute to the catalytic activity of the ene-reductase OPR3, protonating the hydroxyl group of the oxime in the first reduction step.

Redesigning Enzymes for Biocatalysis: Exploiting Structural Understanding for Improved Selectivity

The discovery of new enzymes, alongside the push to make chemical processes more sustainable, has resulted in increased industrial interest in the use of biocatalytic processes to produce high-value and chiral precursor chemicals. Huge strides in protein engineering methodology and in silico tools have facilitated significant progress in the discovery and production of enzymes for biocatalytic processes. However, there are significant gaps in our knowledge of the relationship between enzyme structure and function. This has demonstrated the need for improved computational methods to model mechanisms and understand structure dynamics. Here, we explore efforts to rationally modify enzymes toward changing aspects of their catalyzed chemistry. We highlight examples of enzymes where links between enzyme function and structure have been made, thus enabling rational changes to the enzyme structure to give predictable chemical outcomes. We look at future directions the field could take and the technologies that will enable it.

Asymmetric Hydrogenation of C = C Bonds in a SpinChem Reactor by Immobilized Old Yellow Enzyme and Glucose Dehydrogenase

The application of immobilized enzymes in pharmaceutical and bulk chemical production has been shown to be economically viable. We demonstrate the exceptional performance of a method that immobilizes the old yellow enzyme YqjM and glucose dehydrogenase (GDH) on resin for the asymmetric hydrogenation (AH) of C = C bonds in a SpinChem reactor. When immobilized YqjM and GDH are reused 10 times, the conversion of 2-methylcyclopentenone could reach 78%. Which is because the rotor of the SpinChem reactor effectively reduces catalyst damage caused by shear force in the reaction system. When the substrate concentration is 175 mM, an 87% conversion of 2-methylcyclopentenone is obtained. The method is also observed to perform well for the AH of C = C bonds in other unsaturated carbonyl compounds with the SpinChem reactor. Thus, this method has great potential for application in the enzymatic production of chiral compounds.

Introduction of Chiral Centers to α- and/or β-Positions of Carbonyl Groups by Biocatalytic Asymmetric Reduction of α,β-Unsaturated Carbonyl Compounds

Biocatalytic asymmetric reductions of acyclic and cyclic α,β-unsaturated carbonyl compounds are favorable protocols for introduction of chiral centers to α- and/or β-positions of the carbonyl groups. Representative biocatalytic reductions of electron deficient olefins are compiled from a synthetic point of view according to compound types from the papers in 2012 to early 2022. Applications to syntheses of some enantiomericaly enriched perfumery ingredients are presented to show the feasibility of the biocatalytic reductions.

"A Study in Yellow": Investigations in the Stereoselectivity of Ene-Reductases

Ene‐reductases from the Old Yellow Enzyme (OYE) superfamily are a well‐known and efficient biocatalytic alternative for the asymmetric reduction of C=C bonds. Considering the broad variety of substituents that can be tolerated, and the excellent stereoselectivities achieved, it is apparent why these enzymes are so appealing for preparative and industrial applications. Different classes of C=C bonds activated by at least one electron‐withdrawing group have been shown to be accepted by these versatile biocatalysts in the last decades, affording a vast range of chiral intermediates employed in the synthesis of pharmaceuticals, agrochemicals, flavours, fragrances and fine chemicals. In order to access both enantiomers of reduced products, stereodivergent pairs of OYEs are desirable, but their natural occurrence is limited. The detailed knowledge of the stereochemical course of the reaction can uncover alternative strategies to orient the selectivity via mutagenesis, evolution, and substrate engineering. An overview of the ongoing studies on OYE‐mediated bioreductions will be provided, with particular focus on stereochemical investigations by deuterium labelling.

Biocatalytic Reduction Reactions from a Chemist's Perspective

Reductions play a key role in organic synthesis, producing chiral products with new functionalities. Enzymes can catalyse such reactions with exquisite stereo-, regio- and chemoselectivity, leading the way to alternative shorter classical synthetic routes towards not only high-added-value compounds but also bulk chemicals. In this review we describe the synthetic state-of-the-art and potential of enzymes that catalyse reductions, ranging from carbonyl, enone and aromatic reductions to reductive aminations.

Finding Reactive Configurations: A Machine Learning Approach for Estimating Energy Barriers Applied to Sirtuin 5

Sirtuin 5 is a class III histone deacetylase that, unlike its classification, mainly catalyzes desuccinylation and demanoylation reactions. It is an interesting drug target that we use here to test new ideas for calculating reaction pathways of large molecular systems such as enzymes. A major issue with most schemes (e.g., adiabatic mapping) is that the resulting activation barrier height heavily depends on the chosen educt conformation. This makes the selection of the initial structure decisive for the success of the characterization. Here, we apply machine learning to a large number of molecular dynamics frames and potential energy barriers obtained by quantum mechanics/molecular mechanics calculations in order to identify (1) suitable start-conformations for reaction path calculations and (2) structural features relevant for the first step of the desuccinylation reaction catalyzed by Sirtuin 5. The latter generally aids the understanding of reaction mechanisms and important interactions in active centers. Using our novel approach, we found eleven key features that govern the reactivity. We were able to estimate reaction barriers with a mean absolute error of 3.6 kcal/mol and identified reactive configurations.

Linear enzyme cascade for the production of (-)-iso-isopulegol

Biocatalysis has developed enormously in the last decade and now offers solutions for the sustainable production of chiral and highly functionalised asset molecules. Products generated by enzymatic transformations are already being used in the food, feed, chemical, pharmaceutical and cosmetic industry, and the accessible compound panoply is expected to expand even further. In particular, the combination of stereo-selective enzymes in linear cascade reactions is an elegant strategy toward enantiomeric pure compounds, as it reduces the number of isolation and purification steps and avoids accumulation of potentially unstable intermediates. Here, we present the set-up of an enzyme cascade to selectively convert citral to (-)-iso-isopulegol by combining an ene reductase and a squalene hopene cyclase. In the initial reaction step, the ene reductase YqjM from Bacillus subtilis selectively transforms citral to (S)-citronellal, which is subsequently cyclised exclusively to (-)-iso-isopulegol by a mutant of the squalene hopene cyclase from Alicyclobacillus acidocaldarius (AacSHC). With this approach, we can convert citral to an enantiopure precursor for isomenthol derivatives.

Pd-catalyzed synthesis of α,β-unsaturated ketones by carbonylation of vinyl triflates and nonaflates

A general and highly chemoselective Pd-catalyzed protocol for the synthesis of α,β-unsaturated ketones by carbonylation of vinyl triflates and nonaflates is presented. Applying the specific monophosphine ligand cataCXium® A, the synthesis of various vinyl ketones as well as carbonylated natural product derivatives proceeds in good yields.

Asymmetric Reductive Carbocyclization Using Engineered Ene Reductases

Ene reductases from the Old Yellow Enzyme (OYE) family reduce the C=C double bond in α,β-unsaturated compounds bearing an electron-withdrawing group, for example, a carbonyl group. This asymmetric reduction has been exploited for biocatalysis. Going beyond its canonical function, we show that members of this enzyme family can also catalyze the formation of C-C bonds. α,β-Unsaturated aldehydes and ketones containing an additional electrophilic group undergo reductive cyclization. Mechanistically, the two-electron-reduced enzyme cofactor FMN delivers a hydride to generate an enolate intermediate, which reacts with the internal electrophile. Single-site replacement of a crucial Tyr residue with a non-protic Phe or Trp favored the cyclization over the natural reduction reaction. The new transformation enabled the enantioselective synthesis of chiral cyclopropanes in up to >99 % ee.

In silico and in vitro studies of the reduction of unsaturated α,β bonds of trans-2-hexenedioic acid and 6-amino-trans-2-hexenoic acid - Important steps towards biobased production of adipic acid

The biobased production of adipic acid, a precursor in the production of nylon, is of great interest in order to replace the current petrochemical production route. Glucose-rich lignocellulosic raw materials have high potential to replace the petrochemical raw material. A number of metabolic pathways have been proposed for the microbial conversion of glucose to adipic acid, but achieved yields and titers remain to be improved before industrial applications are feasible. One proposed pathway starts with lysine, an essential metabolite industrially produced from glucose by microorganisms. However, the drawback of this pathway is that several reactions are involved where there is no known efficient enzyme. By changing the order of the enzymatic reactions, we were able to identify an alternative pathway with one unknown enzyme less compared to the original pathway. One of the reactions lacking known enzymes is the reduction of the unsaturated α,β bond of 6-amino-trans-2-hexenoic acid and trans-2-hexenedioic acid. To identify the necessary enzymes, we selected N-ethylmaleimide reductase from Escherichia coli and Old Yellow Enzyme 1 from Saccharomyces pastorianus. Despite successful in silico docking studies, where both target substrates could fit in the enzyme pockets, and hydrogen bonds with catalytic residues of both enzymes were predicted, no in vitro activity was observed. We hypothesize that the lack of activity is due to a difference in electron withdrawing potential between the naturally reduced aldehyde and the carboxylate groups of our target substrates. Suggestions for protein engineering to induce the reactions are discussed, as well as the advantages and disadvantages of the two metabolic pathways from lysine. We have highlighted bottlenecks associated with the lysine pathways, and proposed ways of addressing them.

Engineering the biocatalytic selectivity of iridoid production in Saccharomyces cerevisiae

Monoterpene indole alkaloids (MIAs) represent a structurally diverse, medicinally essential class of plant derived natural products. The universal MIA building block strictosidine was recently produced in the yeast Saccharomyces cerevisiae, setting the stage for optimization of microbial production. However, the irreversible reduction of pathway intermediates by yeast enzymes results in a non-recoverable loss of carbon, which has a strong negative impact on metabolic flux. In this study, we identified and engineered the determinants of biocatalytic selectivity which control flux towards the iridoid scaffold from which all MIAs are derived. Development of a bioconversion based production platform enabled analysis of the metabolic flux and interference around two critical steps in generating the iridoid scaffold: oxidation of 8-hydroxygeraniol to the dialdehyde 8-oxogeranial followed by reductive cyclization to form nepetalactol. In vitro reconstitution of previously uncharacterized shunt pathways enabled the identification of two distinct routes to a reduced shunt product including endogenous 'ene'-reduction and non-productive reduction by iridoid synthase when interfaced with endogenous alcohol dehydrogenases. Deletion of five genes involved in α,β-unsaturated carbonyl metabolism resulted in a 5.2-fold increase in biocatalytic selectivity of the desired iridoid over reduced shunt product. We anticipate that our engineering strategies will play an important role in the development of S. cerevisiae for sustainable production of iridoids and MIAs.

Structural investigation into the C-terminal extension of the ene-reductase from Ralstonia (Cupriavidus) metallidurans

Ene-reductases (ERs), or Old Yellow Enzymes, catalyze the asymmetric reduction of various activated alkenes. This class of biocatalysts is considered an attractive alternative to current chemical technologies for hydrogenation due to their high selectivity and specificity. Here the X-ray crystal structure of RmER, a "thermophilic"-like ER from Ralstonia (Cupriavidus) metallidurans, is reported. Unlike other members of this class of ERs, RmER is monomeric in solution which we previously related to its atypical elongated C-terminus. A typical dimer interface was however observed in our crystal structure, with the conserved Arg-"finger" forming part of the adjacent monomer's active site and the elongated C-terminus extending into the active site through contacting the "capping" domain. This dimerization also resulted in the loss of one FMN cofactor from each dimer pair. This potential transient dimerization and dissociation of FMN could conceivably explain the rapid rates previously observed when an FMN light-driven cofactor regeneration system was used during catalysis with RmER.

Redox-dependent substrate-cofactor interactions in the Michaelis-complex of a flavin-dependent oxidoreductase

How an enzyme activates its substrate for turnover is fundamental for catalysis but incompletely understood on a structural level. With redox enzymes one typically analyses structures of enzyme–substrate complexes in the unreactive oxidation state of the cofactor, assuming that the interaction between enzyme and substrate is independent of the cofactors oxidation state. Here, we investigate the Michaelis complex of the flavoenzyme xenobiotic reductase A with the reactive reduced cofactor bound to its substrates by X-ray crystallography and resonance Raman spectroscopy and compare it to the non-reactive oxidized Michaelis complex mimics. We find that substrates bind in different orientations to the oxidized and reduced flavin, in both cases flattening its structure. But only authentic Michaelis complexes display an unexpected rich vibrational band pattern uncovering a strong donor–acceptor complex between reduced flavin and substrate. This interaction likely activates the catalytic ground state of the reduced flavin, accelerating the reaction within a compressed cofactor–substrate complex.

Due to their transient nature, enzyme-substrate complexes are difficult to characterize structurally. Here, the authors capture the reactive reduced form of xenobiotic reductase A bound to its substrate and show that the oxidation state of the flavin cofactor affects the interaction of the substrate with the enzyme.

Cofactor-Free, Direct Photoactivation of Enoate Reductases for the Asymmetric Reduction of C=C Bonds

Enoate reductases from the family of old yellow enzymes (OYEs) can catalyze stereoselective trans-hydrogenation of activated C=C bonds. Their application is limited by the necessity for a continuous supply of redox equivalents such as nicotinamide cofactors [NAD(P)H]. Visible light-driven activation of OYEs through NAD(P)H-free, direct transfer of photoexcited electrons from xanthene dyes to the prosthetic flavin moiety is reported. Spectroscopic and electrochemical analyses verified spontaneous association of rose bengal and its derivatives with OYEs. Illumination of a white light-emitting-diode triggered photoreduction of OYEs by xanthene dyes, which facilitated the enantioselective reduction of C=C bonds in the absence of NADH. The photoenzymatic conversion of 2-methylcyclohexenone resulted in enantiopure (ee>99 %) (R)-2-methylcyclohexanone with conversion yields as high as 80-90 %. The turnover frequency was significantly affected by the substitution of halogen atoms in xanthene dyes.

Plant tissue cultures as sources of new ene- and ketoreductase activities

While many redox enzymes are nowadays available for synthetic applications, the toolbox of ene-reductases is still limited. Consequently, the screening for these enzymes from diverse sources in the search of new biocatalyst suitable for green chemistry approaches is needed. Among 13 plant tissue cultures, Medicago sativa and Tessaria absinthioides calli, as well as Capsicum annuum hairy roots, were selected due to their ability to hydrogenate the CC double bond of the model substrate 2-cyclohexene-1-one. The three axenic plant cultures showed more preference toward highly activated molecules such as nitrostyrene and maleimide rather than the classical substrates of the well-known Old Yellow Enzymes, resembling the skills of the NAD(P)H-dependent flavin-independent enzymes. When the three biocatalytic systems were applied in the reduction of chalcones, T. absinthioides showed high chemoselectivity toward the CC double bond whereas the other two demonstrated abilities to biohydrogenate the CC double bounds and the carbonyl groups in a sequential fashion.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

How Many Conformations of Enzymes Should Be Sampled for DFT/MM Calculations? A Case Study of Fluoroacetate Dehalogenase

The quantum mechanics/molecular mechanics (QM/MM) method (e.g., density functional theory (DFT)/MM) is important in elucidating enzymatic mechanisms. It is indispensable to study "multiple" conformations of enzymes to get unbiased energetic and structural results. One challenging problem, however, is to determine the minimum number of conformations for DFT/MM calculations. Here, we propose two convergence criteria, namely the Boltzmann-weighted average barrier and the disproportionate effect, to tentatively address this issue. The criteria were tested by defluorination reaction catalyzed by fluoroacetate dehalogenase. The results suggest that at least 20 conformations of enzymatic residues are required for convergence using DFT/MM calculations. We also tested the correlation of energy barriers between small QM regions and big QM regions. A roughly positive correlation was found. This kind of correlation has not been reported in the literature. The correlation inspires us to propose a protocol for more efficient sampling. This saves 50% of the computational cost in our current case.