Structural insights into the γ-lactamase activity and substrate enantioselectivity of an isochorismatase-like hydrolase from Microbacterium hydrocarbonoxydans

(+)-γ-lactamase catalyzes the specific hydrolysis of (+)-γ-lactam out of the racemic γ-lactam (2-Azabicyclo[2.2.1]hept-5-en-3-one) to leave optically pure (−)-γ-lactam, which is the key building block of antiviral drugs such as carbovir and abacavir. However, no structural data has been reported on how the enzymes bind the γ-lactams and achieve their enantioselectivities. We previously identified an isochorismatase-like hydrolase (IHL, Mh33H4-5540) with (+)-γ-lactamase activity, which constitutes a novel family of γ-lactamase. Here, we first discovered that this enzyme actually hydrolyzed both (+)- and (−)-γ-lactam, but with apparently different specificities. We determined the crystal structures of the apo-form, (+)-γ-lactam bound, and (−)-γ-lactam bound forms of the enzyme. The structures showed that the binding sites of both (+) and (−)-γ-lactam resemble those of IHLs, but the “cover” loop conserved in IHLs is lacking in the enzyme, probably resulting in its incomplete enantioselectivity. Structural, biochemical, and molecular dynamics simulation studies demonstrated that the steric clash caused by the binding-site residues, especially the side-chain of Cys111 would reduce the binding affinity of (−)-γ-lactam and possibly the catalytic efficiency, which might explain the different catalytic specificities of the enantiomers of γ-lactam. Our results would facilitate the directed evolution and application of Mh33H4-5540 in antiviral drug synthesis.

. MhIHL is an isochorismatase-like superfamily member with (+)-γ-lactamase activity. (a) Kinetic resolution of (rac)-γ -lactam catalyzed by (+ )-γ -lactamase. (b) Structures of racemic γ -lactam and its reduced forms. (C) Sequence alignment of IHL superfamily members from different species. The sequences of other IHL family members are marked with their respective PDB codes. The secondary structural elements of MhIHL are indicated above its amino acid sequence. The catalytic triad "D13-K78-C111" are marked with stars. but with apparently different specificities (Fig. 2a). To investigate the mechanisms underlying the hydrolysis of γ -lactams by MhIHL and its enantioselectivity, we solved the crystal structure of native MhIHL at 2.05 Å, and the structures of the inactive C111A variant in complex with (+ )-γ -lactam and (− )-γ -lactam at 1.79 Å and 1.99 Å, respectively (Table 1). Structural and biochemical analysis confirmed the similar fold of MhIHL to other IHLs, and suggested a mechanism by which MhIHL facilitates γ -lactam hydrolysis. Further molecular dynamics (MD) simulation results of the enzyme-substrate complex structures clearly explained the different kinetic parameters and catalytic efficiencies of (+ )-and (− )-γ -lactam. Taken together, our study represented the first structure of a (+ )-γ -lactam/lactamase complex and also the first structure of a γ -lactam-bound IHL family member. Structural, biochemical, and MD simulation studies provided the basis of MhIHL-catalyzed hydrolysis of γ -lactam and the enantioselectivity of the enzyme. The study will also shed light on the further directed evolution and application of MhIHL and other IHL-family proteins in antiviral drug synthesis.
Overall Structure of native MhIHL. To better understand the different specificities of MhIHL for (+ )-and (− )-γ -lactam, we firstly solved the high resolution structure of MhIHL (Table 1). Like other IHLs, MhIHL folds into a typical α /β domain with a six-stranded parallel β -sheet (β 1-6) in the middle, flanked by three helices (α 2-4), and a single long helix (α 1) on both sides of the sheet, respectively (Fig. 2c). Relative to the sequence, the order of the strands in the sheet is 3-2-1-4-5-6. An additional 2-strand antiparallel β -sheet (residues 138-140 and 143-145, which corresponds to strands β 1′ and β 2′ in Fig. 2c) resides between β 5 and α 4, which is a unique feature shared by several IHL family proteins 20 . MhIHL crystallized in space group P3 1 21 with one polypeptide chain in the asymmetric unit. The PISA server 23 predicts that MhIHL forms a dimer in solution, consistent with the gel filtration result of MhIHL ( Supplementary Fig. S1). The dimer is formed mainly by the hydrophobic interactions between the 2-strand antiparallel β -sheet from one protomer and α 4 from the other protomer. The total contact area of the monomer is 2526.1 Å 2 , and the decrease in this area upon dimerization is 488.6 Å 2 in one protomer.
Comparison with other IHL proteins. A search in the DALI server 24 with the structure of MhIHL returns several results with Z scores ranging from 17.2 to 26.3, all of which are IHL superfamily members. SCOP 25 divides this superfamily into six subgroups according to their catalyzed reactions which were pyrazinamidase/ nicotinamidase E.C.3.5.1.19, nicotinamidase-related proteins, N-carbamoylsarcosine amidohydrolase (CSHase) E.C. 3.5.1.59, isochorismatase E.C. 3.3.2.1, ribonuclease, and a family of bacterial proteins of unknown function exemplified by the Escherichia coli ycac gene product. The overall "α -β -α sandwich" fold is highly conserved among MhIHL and other IHL proteins ( Supplementary Fig. S2). In addition, a cis peptide bond is observed between residues A106 and Q107, which is a well conserved feature among IHL family members 19 . Out of the IHL proteins, MhIHL shows the highest structural similarity to a putative IHL protein with no defined function from Oleispira antarctica (PDB code 3LQY), with a root-mean-squared displacement (RMSD) of 1.3 Å between the two proteins. However, they are different in the loop regions between β 1 and α 1, β 2 and β 3, and β 6 and α 4 ( Supplementary Fig. S2). PhzD protein from Pseudomonas aeruginosa (PDB code 1NF9) represents one of the best characterized IHL family proteins. PhzD, which belongs to the isochorismatase E.C. 3.3.2.1 subgroup,

Data collection
MhIHL MhIHL-(+)-γ-lactam MhIHL-(−)-γ-lactam  Active Site. Among the six subgroups of the IHL superfamily, three of the groups, namely the pyrazinamidase/nicotinamidase, ribonuclease and N-carbamoylsarcosine amidohydrolase proteins, share a D-K-C catalytic triad, while the others lack the nucleophilic cysteine (for example, D38-K122-G155 in PhzD). MhIHL also has this conserved triad (D13-K78-C111), and mutation of either residue to Ala caused sharp decrease or totally abolishment of the activity of this enzyme (Fig. 2d). Specifically, the D13A mutant showed 1/10 activity of wildtype enzyme (153.7 U/mg) and the K78A and C111A mutant enzymes almost showed no activities (Fig. 2d). This is consistent with the important roles the three residues played in the catalysis. Structures of the complex of MhIHL with either (+ )-γ -lactam or (− )-γ -lactam were obtained by soaking the crystals of the C111A variant with the corresponding substrate, separately (Fig. 3). The active site of MhIHL is located at the C-termini of the six-stranded β -sheet, similar to other known structures with substrates ( Fig. 3a). A clear electron density for one ligand was observed in the active site of each protomer ( Supplementary Fig. S4) in both (+ )-and (− )-γ -lactam-bound MhIHL structures. Surprisingly, (+ )-γ -lactam and (− )-γ -lactam bind to nearly the same position in MhIHL ( Fig. 3b-d). In this case, the backbone nitrogen atoms of these two residues formed the oxyanion hole which could stabilize the oxyanion transition state of the catalyzed reaction. In both structures, the carbonyl oxygen atom of the γ -lactam forms two hydrogen bonds with the main-chain nitrogen atoms of cis-Q107 and A111 (C111 in the native enzyme), respectively ( Fig. 3b-d). In this case, the backbone nitrogen atoms of these two residues formed the oxyanion hole, which could stabilize the oxyanion transition state of the catalyzed reaction. The nitrogen atom of the ligand is hydrogen bonded with the side-chain of D13. The residues in the binding site are mostly conserved among isochorismatases (Fig. 1c). Although the γ -lactam binds to MhIHL through different mechanisms compared to other substrates in respective IHL proteins, the location of (+ )-γ -lactam is typically the same as other substrates ( Supplementary Fig. S3). PhzD contains two aforementioned additional α helices (α 5: Thr83-Gly88 and α 6: Leu90-Gly95) above the substrate cavity, together with the loops flanking them, serving as a "cover", which could possibly slow down the release of the product from the cavity. This feature is conserved among the members whose active cavities are located totally in one protomer, and this "cover" could also be contributed by an adjacent protomer as seen in protozoan IHL proteins 19 . In MhIHL, however, this "cover" is lacking and could not be provided by the other protomer, either ( Supplementary Fig. S3). We suggest that this open confirmation of the active cavity in MhIHL would facilitate both the binding of the substrate and the release of the product molecule, but simultaneously reduce the specificity of the substrate, consistent with the comparable enzymatic kinetic parameters of both enantiomers catalyzed by the MhIHL.
Possible mechanism of γ-lactam hydrolysis. A possible mechanism of MhIHL-catalyzed hydrolysis of (+ )-γ -lactam or (− )-γ -lactam, is presented in Fig. 4, based on the structures of the C111A-(+ )γ -lactam and C111A-(− )-γ -lactam complexes and the biochemical properties and mutational studies of the enzyme. In the mechanism, the carbonyl oxygen atom of the γ -lactam binds to the oxyanion hole formed by the main-chain nitrogen atoms of Q107 and C111 (Fig. 3b-d). In stage I of the modified mechanism, proton abstraction by D13 would generate a C111 thiolate facilitating nucleophilic attack at the carbonyl group, and the resulting tetrahedral intermediate is stabilized by the oxyanion hole. And then, proton donation from D13 to the ring nitrogen atom of the substrate would promote C-N bond cleavage and the tetrahedral intermediate collapses in stage II to give an acyl intermediate. In stage III, water activation by D13, whereby the D13 residue is protonated, generates a nucleophilic hydroxy group to attack the acyl intermediate to produce the product and a thiolate. In stage IV, the product is released and the thiolate accepts a proton from D13 to regenerate a thiol. In the whole mechanism, the K78 keeps its protonated form and does not directly participate in the catalysis process. This catalytic mechanism resembles the mechanism proposed for other α /β hydrolase enzymes such as Acinetobacter baumanii PncA 26 .
In this mechanism, the thiol group of C111 acted as a nucleophilic group in the hydrolytic process. Mutation of C111 to other nucleophilic residues (C111T, C111S, and C111K) all caused sharp decrease of the activity (Fig. 2d), suggesting the indispensible role of thiol group in the active site. The slightly higher activity of C111K than that of C111T/S suggested a possible weak nucleophilic activity of K111 in the active site.  Fig. 5a and b). During the simulation, (+ )-γ -lactam could be well accommodated in the binding site of the WT enzyme and stayed stably in the equilibrium position ( Fig. 5a and g). However, the (− )-γ -lactam deviated from the starting position with a larger RMSD after the 50 ns simulation ( Fig. 5a and h). In addition, the three hydrogen bond interactions ( Fig. 3b and C), which are essential for the binding of the γ -lactam to the enzyme, were also more disturbed in WT-(− )-γ -lactam system than in the WT-(+ )-γ -lactam complex (Fig. 5c-e). The molecular mechanics Poisson-Boltzmann/surface area (MM-PB/SA) calculations also suggested the (+ )-γ -lactam exhibited a better binding energy than (− )-γ -lactam, probably due to the more stable hydrogen bond network ( Table 3). As the nucleophilic attack of the carbonyl group of the γ -lactam by the thiol group of C111 was suggested as an essential step in the reaction, the distance between the sulfur atom of C111 and the carbon atom of the γ -lactam was also investigated in both structures. In the WT-(+ )-γ -lactam complex, the distance was stabilized at around 3.5 Å, a reasonable distance for the nucleophilic reaction to happen, in contrast to the distance greater than 5 Å in the WT-(− )-γ -lactam complex (Fig. 5f).
Based on this result, we speculated that the bulkier sp3 bridged carbon atom C1 of (− )-γ -lactam might generate steric clash with the binding-site residues, especially the sidechain of C111, which resulted in the unstable binding of (− )-γ -lactam and its low catalytic efficiency, which might elucidate the enantioselectivity of MhIHL ( Fig. 5g and h). To test this hypothesis, we also carried out MD simulations for the complex structures of C111A mutant with (+ )-or (− )-γ -lactam ( Supplementary Fig. S5). The results showed that the average RMSD values of ligands and deviations of proteins were both similar between the two complexes ( Supplementary Fig. S5). Moreover, the differences of the occupancies of the three above mentioned hydrogen bonds between (+ )-and (− )-γ -lactam-bound structures were also decreased in the C111A complexes, compared to the WT complexes ( Supplementary Fig. S5), indicating that both (+ )-γ -lactam and (− )-γ -lactam could be well accommodated in the binding pocket of the C111A enzyme ( Supplementary Fig. S5). The gap of the calculated binding energies were also reduced between the systems of C111A mutant with (+ )-and (− )-γ -lactam (Table 3). Moreover, since the saturated C1 atom and unsaturated C3 and C4 atoms constitute the major determinants of enantioselectity, we also tested the hydrolysis of hydrogenation-reduced form of γ -lactam enantiomers (Fig. 1b) by MhIHL (Fig. 5i).
The results showed that MhIHL hydrolyzed both reduced enantiomers at similar rates (Fig. 5i). This suggested that both enantiomers of the reduced forms of γ -lactam might cause similar steric clashes with the enzyme, further supporting our proposal that the saturated C1 atom and unsaturated C3 and C4 atoms constitute the major enantioselectity determinants of γ -lactam by MhIHL.

Discussion
Specific hydrolysis of (+ )-γ -lactam out of racemic γ -lactam to leave optically pure (− )-γ -lactam is of great importance, since (− )-γ -lactam is the key building block of antiviral drugs such as carbovir and abacavir. To date, several (+ )-γ -lactamases have been identified, all of which resemble formamidases and acetamidases, but no structure has been described for the complex between (+ )-γ -lactam and a (+ )-γ -lactamase. The IHLs superfamily is a well-characterized hydrolase family, which is best known in pathogenic organisms 19 . Recently, several structures of isochorismatases or putative isochorismatases have been elucidated. However, most of the isochorismatase homologs were poorly characterized biochemically and enzymatically. Moreover, no γ -lactamase activity has been reported for the members in this family. This study is the first report of the structures and mechanisms of an IHL family member with γ -lactamase activity. Although the quaternary structures observed for members of this family are quite diverse, their tertiary α /β hydrolase fold is highly homologous, and key features of the active site are also conserved. First, most of the family numbers contain a D-K-C catalytic triad (D13-K78-C111 in MhIHL), out of which the thiol group of the cysteine residue acts as a nucleophile. Second, all structures determined to date from this superfamily exhibit a rare nonproline cis-peptide bond (A106-Q107 in MhIHL) in the active site, which has been proposed to help position the substrate-binding residues appropriately 27 . Inspired by this, we propose that members in the IHL family harboring these two features should also have the γ -lactamase activity. Consistently, another three IHL family proteins (designated isochorismatases A, B, and C, respectively) all displayed (+ )-γ -lactamase activities (Supplementary Table 1), although the activities of these enzymes are lower than that of MhIHL. We also tested the activity of an isochorismatases-like protein Yecd (isochorismatase D) which lacks the suggested nucleophilic cysteine (D26-K112-G145 in Yecd). Yecd displayed no γ -lactamase activity, further supporting the indispensible role of a nucleophilic residue, such as cystein, serine, threonine, et al. (Supplementary Table 1).
Interestingly, Hamano et al. reported that the enzyme SttH from Streptomyces albulus could catalyze the hydrolysis of the amide bond of streptothricins (STs), which are broad-spectrum antibiotics. Streptothricins (STs), produced by Streptomyces strains, consist of a carbamoylated D-gulosamine to which the β -lysine homopolymer (1 to 7 residues) and the amide form of the streptolidine (streptolidine lactam) are attached 28 . The authors found that upon hydrolysis of the amide bond of streptolidine lactam, the toxicity of ST-D (possessing three-β -lysine moiety) was altered from broad-spectrum to bacteria-specific. They showed that hydrolyzed form of ST-D (ST-D-acid) exhibits potent antibacterial activity but its toxicity against eukaryotic cells is reduced. SttH is also a 276-amino-acid IHL family member, with 28% identity to MhIHL among 178 residues. A catalytic D-K-C triad (D70-K143-C176) also exists in SttH, and it will be interesting to investigate whether SttH catalyzes the hydrolysis of the amide bond of streptolidine lactam through the same mechanism as MhIHL, although the binding mechanism must be different due to the more complex structure of ST-D.
MhIHL exhibited different hydrolysis activities towards (+ )-and (− )-γ -lactam, with a better specificity for (+ )-γ -lactam. Therefore, in the beginning of the hydrolysis reaction of racemic γ -lactam, (+ )-γ -lactam might act both as a substrate of MhIHL, and also an "inhibitor" for (− )-γ -lactam. That explains why it seems that MhIHL first hydrolyzed the (+ )-γ -lactam, and then the (− )-γ -lactam (Fig. 2a). To investigate the mechanism underlying the different specificities towards (+ )-γ -lactam and (− )-γ -lactam, we solved the structures of apo-form MhIHL, and both substrate-bound complex with the C111A mutant enzyme, and performed MD simulation studies. The structures revealed that both γ -lactam enantiomers bind to almost the same position and their nitrogen and oxygen atoms are stabilized by similar hydrogen bonds. However, MD simulation results suggested that the bulkier sp3 bridged carbon atom C1 of (− )-γ -lactam might generate steric clash with the binding-site residues of MhIHL, especially the side chain of C111, which would decrease the binding specificity of the (− )-γ -lactam in MhIHL and possibly the catalytic efficiency as well. In the (+ )-γ -lactam-bound complex structure, however, the thiol group of C111 is facing the unsaturated carbon atoms C3 and C4, leading to less steric clash. Importantly, although the differences of the occupancies of the three hydrogen bonds between (+ )-and (− )-γ -lactam-bound structures were decreased in the C111A complexes, (− )-γ -lactam still could not be accommodated well in the binding cavity ( Supplementary Fig. S5), suggesting that other residues in the binding site could also contribute to the enantioselectivity of MhIHL. In addition, the similar hydrolysis rate of MhIHL towards the reduced forms of (+ )-and (− )-γ -lactam (Fig. 5i) further supported that the saturated C1 atom and the unsaturated C3 and C4 atoms might be the determinants of the enantioselectivity towards racemic γ -lactam of MhIHL.
In all, we have presented three structures of MhIHL and its complex with both (+ ) and (− )-γ -lactam substrate and conducted biochemical and mutational studies into the catalytic mechanism and enantioselectivity of MhIHL. MhIHL from the isochorismatase-like superfamily, represents a novel class of (+ )-γ -lactamases, thus further expanding the (+ )-γ -lactamase superfamily. MhIHL displayed different hydrolysis activities towards (+ )-and (− )-γ -lactam, and our structural, biochemical, and MD simulation studies demonstrated the underlying mechanisms for the enantioselectivity of MhIHL. Due to the wide applications of (− )-γ -lactam in antiviral drugs, this study provided much needed insight into the catalytic mechanism of a new type of (+ )-γ -lactamase and would shed light on the molecular manipulation of (+ )-γ -lactamases for higher efficiency and enantioselectivity. Expression and purification of proteins. E. coli Rosetta (DE3) cells harboring recombinant plasmids were grown in LB medium containing ampicillin (50 μ g/mL). The cells were induced with isopropyl β -D-thiogalactopyranoside (IPTG) at a final concentration of 0.1 mM when the absorbance reached 0.8-1.0 at 600 nm. The incubation was then continued at 30 °C for another 5 h. The cells were harvested by centrifugation at 3,500 × g for 10 min, resuspended in 50 mM Tris-HCl (pH 7.8), and disrupted by sonication. The cell debris was removed by centrifugation at 12,000 × g for 50 min. The cell-free extract obtained was applied onto a nickel nitrilotriacetic acid (Ni-NTA) agarose (GE Healthcare) column equilibrated with 20 mM imidazole in 50 mM Tris-HCl buffer (pH 8.0). Next, the enzyme was further purified on a Superdex 200 column (GE Healthcare) with a 50 mM Tris-HCl buffer (pH 7.8) containing 300 mM NaCl and 2 mM β -mercaptoethanol using AKTA FPLC system (GE Healthcare). The protein concentration was determined spectrophotometrically at 280 nm using the Thermo NanoDrop 2000 and further calibrated according to the theoretical extinction coefficient based on the amino acid sequences. The MhIHL mutants and other isochorismatases were expressed and purified in a manner identical to that used for the WT MhIHL enzyme.

Measurement of enzymatic activity and kinetic constants.
The catalytic activities of MhIHL and mutants were examined by the addition of 0.2 μ g pure MhIHL to 1000 μ L of the 46 mM substrate solution (50 mM Tris-HCl, pH 7.5). The catalytic activities of isochorismatase A, B, C, and D were detected by the addition of 80 μ g pure enzyme to 1000 μ L of the 46 mM substrate solution (50 mM Tris-HCl, pH 7.5). The reactions were incubated at 25 °C for 15 min and then extracted with 1000 μ L ethyl acetate. The (+ )-γ -lactamase activity was detected as described previously 10 . One unit (U) of enzyme activity was defined as the amount of enzyme that catalyzed the conversion of one nanomole of substrate per minute. Measurements of kinetic parameters for bioresolution reactions were performed at 25 °C. The hydrolyzing reactions were performed in 50 mM Tris-HCl (pH 7.5) with Table 3. Calculated binding energy (Unit: kcal/mol) by molecular mechanics Poisson-Boltzmann/surface area (MM-PB/SA).
(− )-γ -lactam or (+ )-γ -lactam as the substrates. The concentration of the γ -lactam was varied from 5 to 200 mM. All the experiments were conducted with three replicates. Km and V max values were calculated by nonlinear fitting. Conversion and enantiomeric excess (ee) were calculated as previously described 18 . Crystallization. All the crystals were obtained at 20 °C using the sitting drop vapor diffusion method by mixing 1 μ L of protein solution (15 mg/mL) with an equal volume of reservoir solution containing 0.2 M sodium malonate (pH 6.0) and 20% (w/v) polyethylene glycol 3,350. Crystals grew to a full size in two or three days. Crystals were cryoprotected in the reservoir solutions supplemented with 25% (w/v) glycerol, flash cooled in the liquid nitrogen, and stored at cryogenic temperatures for data collection. Both co-crystallization and soaking methods were conducted to obtain the enzyme-substrate structure. The final crystals of MhIHL-(− )-γ -lactam complex were obtained by soaking the crystal of the substrate-free C111A mutant enzyme in a cryoprotectant containing 250 mM (− )-γ -lactam for 1.5 h. The final MhIHL-(+ )-γ -lactam complex crystals were obtained by soaking the crystal of substrate-free C111A mutant enzyme in cryoprotectant containing 500 mM (+ )-γ -lactam for 10 min.
Data collection, and structure determination. Diffraction data were collected at the Shanghai Synchrotron Radiation Facility (SSRF), integrated and scaled using HKL2000 29 . The initial phase calculation, phase improvement, and automated model building were performed using the molecular replacement (MR) method with PHASER in the CCP4 program suite 30 and PHENIX suite 31 . The structure of an isochorismatase (Pp1826) from Pseudomonas putida Kt2440 (PDB code 4H17) was used as the search model. The final model rebuilding was performed using COOT 32 and the protein structure was refined with PHENIX 33 against the data using stereochemistry information. For the complex of C111A mutant with either (+ )-γ -lactam or (− )-γ -lactam, the molecular replacement method was also used with the structure of the native enzyme as the search model. The subsequent model building and refinement was performed as mentioned for the native enzyme. Data-collection and refinement statistics are presented in Table 1. All structural model figures were created with PyMOL 34 .

Molecular dynamics (MD) simulation.
To further explore the molecular mechanism of enantioselectivity of the enzyme, MD simulations were carried out for the complex structures of the WT/C111A enzyme with (+ )-γ -lactam and (− )-γ -lactam respectively. The initial structures of the two WT complexes were modeled based on the C111A mutant-(+ )/(− )-γ -lactam crystal structures using Protein Local Optimization Program (PLOP) [35][36][37] . The MD simulations were carried out using AMBER10.0 suite 38 . The AMBER99SB force filed 39 was applied for protein and the general Amber force field (GAFF) 40 for the ligands. ANTECHARMBER was used for calculation of the AM1-BCC charges of ligands. All system setups were performed in TLEAP module. A 12 Å pad of TIP3P waters was added to solvate to each system. Neutralizing counter ions were added to each system. Simulations were carried out using the SANDER module. Each system was minimized using 5000 steps of steepest descent algorithm with gradually reduced restraint. Following minimization, the system was gradually heated to 300 K in the canonical NPT ensemble using a Berendsen thermostat 41 , with a 50 kcal/mol/Å 2 harmonic position restraint applied to all heavy atoms of the protein. The restraints were gradually removed and the production run was performed in MTK-NPT (1 bar, 300 K) ensemble for 50 ns (WT enzyme) and 10 ns (C111A enzyme), respectively. The SHAKE algorithm 42 was applied to constrain all bonds involving hydrogen atoms with a time step of 2 fs.