Characterization of TESH
The optimum temperature and pH of TESH from P. koreensis BK-9 were 45°C and pH 9.0, respectively (Fig. 1A, 1B). However enzyme activity was almost lost when incubated at 45°C for 30 min, only 10% of residual activity remained when pH was lower than 5 or pH was higher than 10 (Fig. 1A, 1C), indicating poor temperature and pH stability.
Insert Fig. 1
Statistical analysis indicated that Mg2+, Ca2+ were the activators of the enzyme, increased TESH activity significantly (P < 0.001). Zn2+and Mn2+ had strongly inhibited the activity of TESH ( P< 0.001). The addition of EDTA reduced the activity of the enzyme by 46.63%. SDS was a strong inhibitor, decreasing 95.85% activity of the enzyme (Fig. 1D). The kinetic parameters were determined by Michaelis-Menten plotting method (GraphPad Prism 9.5). The vales of Km and kcat were calculated to be 3.50 mM and 94.75 s− 1, respectively (Fig. S1).
Sequence comparison of TESH and reported EHs
The sequence alignment results indicated that TESH shows no homology to other reported ESHs (ESH(L) belongs to HAD like superfamily, ESH(D) belongs to BKACE superfamily) that belongs to MhpC superfamily. We found several structurally characterized proteins with low but considerable sequence identity (between 19% and 46%) were shown in a Clustal Omega (Fig. 2). Several residues were highly conserved: nine nonpolar amino acids (G23, G25, G33, G61, G63, L48, L59, A49, P235), two arginine (D58 and D104), two histidine (H103 and H272) and one lysine (Y211), polar amino acids play an important role in enzyme catalysis. Besides, E35 and H104 form ion pairs in StEH1, so we also designed corresponding Q36 mutated amino acid in TESH.
Insert Fig. 2
Homology modeling of TESH
We selected 4NVR as the template structure to generate the 3D structure of the TESH from P. koreensis BK-9 with Swiss-Model, this template possesses the greatest amino acid sequence identity with TESH (46.52%). The root mean square deviation value of the modeled TESH was 0.068 Å (Fig. S2). The Ramachandran plot (Fig. S3) showed that the percentages of residues falling in disallowed regions, generously allowed regions, favorable regions and core regions were 0, 1.2, 9.3, and 89.5%, respectively. An overall plot showed that more than 95% of the residues were within the favorable region, which revealed that the quality of the model was good.
The modeled 3D structure of TESH has the typical eight parallel β-sheets (residues 6–10, 16–22, 27–31, 53–57, 98–103, 121–127, 236–241 and 263–268) surrounded by α-helices and is covered by a cap-domain (Fig. 3A). It is similar with the structure from MhpC superfamily (StEH1, AD1, CorEH, and FAcD).
Insert Fig. 3
The trans-epoxysuccinate substrate was docked in the core of the modeled TESH structure, results showed that there were thirteen amino acids (H34, D104, R105, R108, D128, I129, Y147, H149, W150, V175, R179, Y211, H272) near the binding pocket, which were labeled on the modeled structure by PyMOL (Fig. 3B).
Site-directed mutagenesis
Combing the results of multiple sequence alignment and molecular docking, we selected fourteen amino acid sites to replaced their polar amino and conserved charged acid residues with another residues as follows: H and D by N; Q by E; R by K; I, V and W by A; Y by F. Purified protein were showed in Fig. S4. The results of site-directed mutagenesis revealed that six mutants showed relative activities that were only 2% or less of the wild-type enzyme’s activity: R108K, D128N, Y147F, H149N, W150A and Y211F. Four mutants showed the lower kcat/Km value than the wild-type: H34N, D104N, R105K, Y211F (Table 1). The results revealed that the ten residues (H34, D104, R105, R108, D128, Y147, H149, W150, Y211 and H272) played important roles in the catalysis.
Table 1
Characterization of wild-type and mutant TESH enzyme.a
Enzyme | Relative activity(%) | Km(mM) | kcat(s− 1) | kcat/Km(mM− 1 s− 1) |
Wile-type | 100.00 ± 2.34 | 3.50 | 94.75 | 27.08 |
H34N | 5.75 ± 0.39 | 7.88 | 4.73 | 0.60 |
Q36E | 13.40 ± 0.67 | 2.80 | 50.41 | 18.02 |
D104N | 54.53 ± 0.96 | 83.75 | 11.41 | 0.14 |
R105K | 7.94 ± 0.26 | 17.53 | 9.06 | 0.52 |
R108K | 0.20 ± 0.06 | - | - | - |
D128N | 0.53 ± 0.04 | - | - | |
I129A | 8.64 ± 0.83 | 10.48 | 21.25 | 2.03 |
Y147F | 1.81 ± 0.09 | - | - | - |
H149N | 0.34 ± 0.00 | - | - | - |
W150A | 0.55 ± 0.12 | - | - | - |
V175A | 20.79 ± 0.40 | 6.93 | 16.88 | 2.44 |
R179K | 19.04 ± 0.60 | 13.20 | 22.09 | 1.67 |
Y211F | 0.15 ± 0.03 | 2.82 | 0.42 | 0.15 |
H272N | 11.21 ± 0.20 | 3.78 | 15.39 | 4.07 |
aThe specific activity of the wild-type was 124.38 ± 2.34 µmol min− 1 mg− 1. |
Insert Table 1
Single and multiple turnover reaction of wild-type TESH in H218O
The single and multiple turnover reactions of TESH from P. koreensis BK-9 in H218O showed that most of the tartrate was 16O-labeled in the single turnover reaction (Fig. 4A), while most of the tartrate was 18O-labeled in the multiple turnover reaction (Fig. 4B). These results suggested that the catalytic mechanism of the TESH from P. koreensis BK-9 same as most α/β hydrolases that proceeded by a two-step catalytic reaction involving the forming of an enzyme-substrate ester intermediate.
Insert Fig. 4
The ten important catalytic residues of the TESH we obtained were aligned with the dehalogenase from R. palustris CGA009 (FAcD), epoxide hydrolase from Corynebacterium sp. C12 (CorEH), epoxide hydrolase from A. radiobacter (AD1), epoxide hydrolase from S. tuberosum (StEH1) and CFTR inhibitory factor Cif from P. aeruginosa (Cif). We found that these ten residues were similar to the active sites of the FAcD, CorEH, AD1, StEH1 and Cif (Fig. 5A), which have a two-step catalytic mechanism characterized by an aspartate-histidine-aspartate catalytic triad (Nardini et al. 1999; Chan et al. 2011; Bahl et al. 2016; Bauer et al. 2016; Schuiten et al. 2021). Therefore, the TESH may consists of the catalytic nucleophile Asp104, the charge relay acid Asp128 and the histidine base His272 to forming an Asp104-His272-Asp128 catalytic triad, which is shown in Fig. 5B.
Insert Fig. 5
Like most of EHs, Tyr residue that near the substrate plays an important role in help opening the ring. Firstly, they have hydrogen bond with the epoxide oxygen, which can position the substrate in the active site for nucleophilic attack. Then, tyrosine provide a proton to the epoxide oxygen during the ring-opening step (Yamada et al. 2000). While epoxide hydrolases usually use Tyr-Tyr pair to assist ring-opening (Nardini et al. 1999; Loo et al. 2006; Bauer et al. 2016), some epoxide hydrolases use a His-Tyr pair instead (Chan et al. 2011; Bahl et al. 2012; Schuiten et al. 2021). We propose Tyr211 co-action with Tyr147-His149-Trp150 to open the ring of trans-epoxysuccinate (Fig. 5C). When nucleophilic Asp104 attack oxirane ring in the opposite side forming enzyme-substrate intermediate, Tyr211 co-action Tyr147-His149-Trp150 to form an oxyanion hole that stabilizes the oxyanion of the alkyl-enzyme intermediate formed by ring opening (Amrein et al. 2015).
In FAcD Arg111-Arg114 pair are the binding sites for carboxylic acids (Chan et al. 2011), in CorEH Trp100 is the halide-stabilizing residue (Schuiten et al. 2021), Trp38 and Phe108 have hydrogen bonds with Asp107 thus stabilize the intermediate (Nardini et al. 1999). Glu35-His104 ion pair in StEH1 can control solvent access to the active site (Amrein et al. 2015). Same as AD1, CFTR inhibitory factor Cif also have hydrogen bond between Phe63 and Ile130 which can stabilize intermediate (Bahl et al. 2016). According to homologous molecular docking results, Arg105 and Arg108 from TESH have hydrogen bond with the carboxylate binding site of trans-epoxysuccinate (Fig. 5D). H34 have polar contacts with D104 and Y211 (Fig. 5E). We propose like StEH1, H34 in TESH also play the same role as E35 is the part of a hydrogen-bond network that connects the active site with the solvent interface (Amrein et al. 2015). More experiments are needed to confirm this conjecture.
Based on the results and the above analysis, we speculated that the TESH adopted a two-step catalytic mechanism via Asp104-His272-Asp128 catalytic triad (Fig. 6). The catalytic mechanism of ESH(meso) is similar to that of ESH(L) from R.opacus and Klebsiella sp. BK-58, which are also proceeded through the two-step mechanism involving the formation of a covalent intermediate (Pan et al. 2011; Cheng et al. 2014a). But it is different from the ESH(D) from Bordetella sp. BK-52, which catalyzes by means of a Zn2+-dependent, one-step mechanism (Dong et al. 2018).
Insert Fig. 6