Structure-guided steric hindrance engineering of Bacillus badius phenylalanine dehydrogenase for efficient l-homophenylalanine synthesis
Biotechnology for Biofuels volume 14, Article number: 207 (2021)
Direct reductive amination of prochiral 2-oxo-4-phenylbutyric acid (2-OPBA) catalyzed by phenylalanine dehydrogenase (PheDH) is highly attractive in the synthesis of the pharmaceutical chiral building block l-homophenylalanine (l-HPA) given that its sole expense is ammonia and that water is the only byproduct. Current issues in this field include a poor catalytic efficiency and a low substrate loading.
In this study, we report a structure-guided steric hindrance engineering of PheDH from Bacillus badius to create an enhanced biocatalyst for efficient l-HPA synthesis. Mutagenesis libraries based on molecular docking, double-proximity filtering, and a degenerate codon significantly increased catalytic efficiency. Seven superior mutants were acquired, and the optimal triple-site mutant, V309G/L306V/V144G, showed a 12.7-fold higher kcat value, and accordingly a 12.9-fold higher kcat/Km value, than that of the wild type. A paired reaction system comprising V309G/L306V/V144G and glucose dehydrogenase converted 1.08 M 2-OPBA to l-HPA in 210 min, and the specific space–time conversion was 30.9 mmol g−1 L−1 h−1. The substrate loading and specific space–time conversion are the highest values to date. Docking simulation revealed increases in substrate-binding volume and additional degrees of freedom of the substrate 2-OPBA in the pocket. Tunnel analysis suggested the formation of new enzyme tunnels and the expansion of existing ones.
Overall, the results show that the mutant V309G/L306V/V144G has the potential for the industrial synthesis of l-HPA. The modified steric hindrance engineering approach can be a valuable addition to the current enzyme engineering toolbox.
l-Homophenylalanine (l-HPA) is an unnatural aromatic amino acid and contains a single-carbon extended side chain compared with l-phenylalanine . l-HPA can serve as a versatile chiral building block for the synthesis of several blockbuster pharmaceuticals, such as angiotensin-converting enzyme inhibitors [2, 3], proteasome inhibitors [4, 5], acetylcholinesterase inhibitors , and β-lactam antibiotics . Chemical methods of l-HPA synthesis are constrained by lengthy steps, high cost, environmental pollution, and poor enantioselectivity, rendering industrial processes ultimately unsustainable [7, 8]. By contrast, biocatalysts are environmentally friendly, mild reaction conditions, and have excellent selectivity [9, 10]. Thus, various enzymatic routes including kinetic resolution catalyzed by lipase , in situ racemization catalyzed by hydantoinase [12, 13], and asymmetric transamination catalyzed by transaminase [14, 15] were recently developed and applied to l-HPA synthesis. Nevertheless, the foregoing reactions also have limitations, such as a theoretical maximum 50% lipase yield, poor applicability because of the strict d-enantioselectivity of hydantoinase, and undesirable equilibrium shifts in the reversible transaminase transamination .
Direct asymmetric reductive amination of prochiral 2-oxo-4-phenylbutyric acid (2-OPBA) catalyzed by phenylalanine dehydrogenase (PheDH) is a promising approach in l-HPA synthesis. It consumes only free ammonia and generates only water as a byproduct (Fig. 1). Asano et al.  reported the use of PheDH from Bacillus sphaericus (BsPheDH) to prepare l-HPA. Bradshaw et al.  investigated the application of recombinant E. coli cells harboring PheDH from Rhodococcus sp. M4 (RsPheDH) to produce l-HPA. They achieved 63% conversion using 60 mM substrate. Ahmad et al.  enhanced the applicability of RsPheDH in l-HPA production via cell immobilization. They achieved > 80% conversion using 67 mM substrate. Notably, Zhang et al.  stated that the combination of substrate fed-batch with continuous product removal improved substrate loading in l-HPA production catalyzed by PheDH from Thermoactinomyces intermedius (TiPheDH). A final substrate concentration of 510 mM was achieved after eight cycles of a fed-batch operation. An important constraint in PheDH-catalyzed l-HPA production is the poor catalytic efficiency toward non-native, bulky 2-OPBA compared with native phenylpyruvic acid (PPA). Consequently, the substrate loading is low and fails to meet the requirements of industrially and commercially viable biocatalysts and biocatalysis processes .
Directed evolution techniques are powerful tools for generating prolific sources of biocatalysts [22, 23]. Steric hindrance of the substrate-binding pocket plays a significant role in modulating the catalytic properties of an enzyme, including enantioselectivities [24, 25], substrate specificities [26,27,28], and catalytic activities [25, 28,29,30]. Recent enzyme engineering studies focused on increasing or conferring enzyme activity towards non-native bulky substrates by attenuating the steric hindrance of the substrate-binding pocket. The mutant constructed by Merck and Codexis  using transaminase ATA-117 as a mutagenic template is a remarkable example of such engineering, as the mutant successfully showed detectable reductive amination activity toward the bulky prositagliptin ketone that lacked any activity before, which was successfully used to catalyze the synthesis of the antidiabetic compound sitagliptin. Chen et al.  used rational design based on steric hindrance attenuation to induce asymmetric reductive amination activity of engineered amine dehydrogenases towards several previously unreactive long-chain α-aliphatic ketones. The hydroxynitrile lyase mutants engineered by Zheng et al.  showed remarkably increased asymmetric hydrocyanation activity toward the rigid pharmaco aldehydes, and successfully applied them to achieve the synthesis of the β2-adrenoreceptor agonist salmeterol. These successful advances demonstrated the feasibility of steric hindrance engineering of substrate-binding pockets to modulate the catalytic activity of enzymes towards non-native, bulky substrates.
In the present study, we report the structure-guided steric hindrance engineering of PheDH from Bacillus badius (BbPheDH) to improve its catalytic efficiency in l-HPA synthesis. The steric hindrance of the substrate-binding pocket of BbPheDH was subjected to crucial residue identification and evaluation of the mutagenesis libraries constructed based on degenerate codons. Seven superior mutants were generated after four rounds of steric hindrance engineering. The optimal triple-site mutants significantly enhanced catalytic efficiency, substrate loading, and specific space–time conversion. Structure analysis elucidated the factors contributing to the enhanced catalytic efficiency observed in the enzyme-kinetics and transformation experiments.
Results and discussion
Preliminary activity studies and molecular modeling
Initially, two PheDHs from B. badius (BbPheDH) and G. kaustophilus (GkPheDH) were selected as candidate enzymes in the development of an enhanced biocatalyst for l-HPA synthesis. As 2-OPBA had a larger side chain than the native substrate PPA, its catalysis by BbPheDH and GkPheDH was characterized by lower turnover frequency (kcat) and less favorable catalytic efficiency (kcat/Km) (Table 1), which motivated us to consider enhancing catalytic efficiency toward 2-OPBA by the steric hindrance engineering. Biotransformation experiments were subsequently performed using BbPheDH and GkPheDH as biocatalysis, giving 67.5% and 53.6% conversions at 0.2 M substrate concentration within 8 h (Table 1 and Additional file 2: Figure S1). Due to its relatively higher conversion, BbPheDH was selected as the starting mutagenesis template for steric hindrance engineering.
To assess the feasibility of developing an enhanced biocatalyst for efficient l-HPA synthesis, we first generated a structural homology model of BbPheDH (Fig. 2a) based on the reported crystal structure of Rhodococcus sp. M4 PheDH  to develop hypotheses for mutagenesis library designs. The model revealed that the residues surrounding the catalytic active center of BbPheDH could be classified into two categories on the basis of their function in the substrate-binding process. Group one encompasses the catalytic residues that are appointed in the binding of the α-carboxylic and α-amino groups of the substrate, namely, K78, K90, D125, and N276, respectively, which were highly conserved among the PheDH family (Additional file 2: Figure S2). These residues create a hydrophilic environment that combines the hydrophilic α-amino acid moiety of the substrate (Fig. 2b). On the opposite side, a panel of hydrophobic residues creates a hindered hydrophobic environment forcing the substrate side chain into its ideal binding pose (Fig. 2c). Furthermore, limited structure-guided engineering was previously conducted on PheDHs and showed that most of the beneficial mutations introduced into the substrate-binding pocket were hydrophobic amino acid residues [35,36,37]. Based on the foregoing information, we selected amino acids with strong hydrophobicity and weak steric hindrance as target substitutes in the subsequent mutagenesis experiments.
Mutagenesis library construction and evaluation based on attenuated steric hindrance
To identify mutagenesis candidates for initial library construction, 2-OPBA was docked into the catalytic active center of BbPheDH to investigate the unfavorable hindrance interactions. Based on the acquired binding pose, all amino acid residues except the smallest glycine and the conserved catalytic residues within 6 Å of the side-chain phenyl ring of 2-OPBA were selected for site-directed mutagenesis in the first round of steric hindrance engineering (Fig. 3a). Twenty-three single-site mutants were constructed at nine amino acid sites based on the designed degenerate codon (Additional file 1: Table S1) by whole-plasmid PCR . All target mutants were isolated from the mixed libraries and confirmed by DNA sequencing. Three superior mutants L50V (M1–1), V309A (M1–2), and V309G (M1–3) were identified by initial activity screening of cell-free extracts (Additional file 2: Figure S3). Subsequent enzyme-kinetics experiments with the purified proteins demonstrated that the Km value of these mutants toward 2-OPBA did not change significantly (Fig. 3h), whereas the kcat values were 25, 61, and 89 s−1, which were 2.3-, 5.7-, and 8.3-fold higher than that of the wild type (11 s−1), respectively (Fig. 3g).
As mutant M1–3 had the highest catalytic efficiency (74 s−1 mM−1) towards 2-OPBA, it was selected as the mutagenesis template in the second round of steric hindrance engineering to investigate the potential synergistic effects. The binding pose of M1–3 with 2-OPBA indicated eight amino acid residues within 6 Å of the side-chain phenyl ring of 2-OPBA. For the rapid detection of any synergy of a proximal steric hindrance with residue G309, the residues were simplified by double-proximity filtering. We identified four mutagenesis candidates within 6 Å of both the phenyl ring of 2-OPBA and G309 (Fig. 3b). Thus, a mutagenesis library consisting of ten double-site mutants was constructed, and the superior mutants V309G/L50V (M2–1) and V309G/L306V (M2–2) were identified by the initial activity screening of cell-free extracts (Additional file 2: Figure S4). The kcat/Km values of the purified mutant proteins, M2–1 and M2–2, toward 2-OPBA were determined to be 143 and 179 s−1 mM−1, respectively, which were 8.4- and 10.5-fold higher than that of the wild type, respectively (Fig. 3f). Notably, this change was also caused by the enhanced kcat value of the mutants (Fig. 3g).
With this result, the mutant M2–2 was selected as the template in the third round of steric hindrance engineering. Here, we investigated proximal steric hindrance synergy with V306. Three mutagenesis candidates were identified by double-proximity filtering using V306 as the key residue (Fig. 3c). The initial activity screening of the cell-free extracts of the ten triple-site mutants revealed two superior mutants V309G/L306V/V144A (M3–1) and V309G/L306V/V144G (M3–2) (Additional file 2: Figure S5). The kcat values of the purified mutant proteins, M3–1 and M3–2, toward 2-OPBA reached 131 and 136 s−1, respectively, which were 12.2- and 12.7-fold higher than that of the wild type, respectively (Fig. 3g); this change was accompanied by a 11.3- and 12.9-fold increase in the kcat/Km values, respectively (Fig. 3f).
Inspired by the above results, we explored the possibility of increasing catalytic activity towards 2-OPBA using M3–2 as the template in the fourth round of steric hindrance engineering. Two mutagenesis candidates were identified (Fig. 3d) and four corresponding quadruple-site mutants were constructed. However, none of them exhibited higher activity towards 2-OPBA than M3–2 (Additional file 2: Figure S6).
Together, seven superior mutants were successfully screened from the 47 elements in the site-directed mutagenesis library (Fig. 3e). The optimal triple-site mutant, M3–2 (V309G/L306V/V144G), showed a 12.7-fold increase in the kcat value (Fig. 3g) compared with that of the wild-type BbPheDH, whereas no obvious change was observed in the Km value (Fig. 3h), which resulted in a 12.9-fold increase in the kcat/Km value (Fig. 3f). This result was similar to that of other studies [25, 32]. Subsequently, the kinetic parameters of BbPheDH and obtained mutants with superior activity toward NADH were assessed. The results showed that the catalytic efficiency of the mutants toward NADH was improved to varying degrees, which is beneficial for the practical application of these mutants to l-HPA synthesis. Moreover, the thermostabilities of BbPheDH and its superior mutants were assessed by measuring the T5030 value (the temperature at which 50% of enzyme activity is lost following a heat treatment for 30 min). It is not unexpected that some kind of stability-activity tradeoff occurred in all mutants, although this tradeoff is of no practical importance, since the T5030 value of all mutants still exceeded 50 ℃ (Additional file 1: Table S2).
Biocatalytic optimization analysis on the best mutant M3–2
The dependence of temperature and pH value were investigated for the reductive amination of 2-OPBA catalyzed by M3–2 paired with glucose dehydrogenase (GluDH) at 0.2 M substrate concentration. Figure 4a shows that 2-OPBA conversion was > 99% after 180 min when the reactions proceeded at 30–40 °C. Catalytic efficiency decreased slightly at 25 °C, with the reaction requiring 240 min to attain 99% conversion. At 50 °C, however, quantitative conversion was < 90% after 240 min. We then evaluated the effects of NH4OH/HCOONH4 buffer (1 M) at various pH values and 30 °C. Figure 4b shows that the fastest conversions occurred at pH 8.5 and the quantitative conversion was > 99% within 90 min. We subsequently performed the reductive amination of 2-OPBA in various concentrations of NH4OH/HCOONH4 buffer (pH 8.5) and NAD+ to determine the influences of NH4+ and NAD+ concentration. The 2-OPBA conversion was > 99% at 1–3 M NH4+. By contrast, 4 M and 5 M NH4+ lowered the conversion rate (Fig. 4c). An activity assay indicated that low GluDH activity might explain the observed decrease in conversion rate in the presence of high NH4+ concentrations (Additional file 2: Figure S7). Moreover, 2-OPBA could be fully converted when the NAD+ concentration was 0.3 mM, 0.5 mM, or 1 mM; however, the conversion dropped to < 50% when the NAD+ concentration was < 0.3 mM (Fig. 4d).
Preparative scale synthesis of l-HPA employing mutant M3–2
Using the optimized reaction system (30 °C; pH 8.5; 3 M NH4+; 0.3 mM NAD+), the preparative scale synthesis of l-HPA by reductive amination of 2-OPBA was carried out in 100 mL reaction volume. The pH value was maintained at 8.5 by titrating concentrated NH3·H2O during the reaction course. Figure 4e shows that 0.2 M and 0.3 M 2-OPBA were converted by > 99% conversion within 20 min and 30 min, respectively. However, substrate inhibition was observed at high concentrations. The quantitative conversion rates were < 70% and < 30% after 240 min in the presence of 0.4 M and 0.5 M substrate, respectively (Additional file 1: Table S3). Similar results were reported for TiPheDH from T. intermedius . We speculate that high substrate concentration may have an irreversible adverse effect on the enzyme .
Continuous substrate fed-batch l-HPA synthesis
Substrate concentration is a crucial factor in industrially feasible biocatalytic processes. Therefore, a combination of substrate fed-batch with product removal strategy was applied for continuous l-HPA synthesis. During the reaction course, 2-OPBA was repeatedly added to the initial NH4OH/HCOONH4 reaction medium (pH 8.5) containing 0.3 mM NAD+, 3 M NH4+, 10 g L−1 cell-free extracts of M3–2, and 12 g L−1 cell-free extracts of GluDH. The substrate concentration was fixed at ~ 0.3 M to maintain high reaction efficiency. Meanwhile, l-HPA product was continuously precipitated from the reaction solution and isolated after filtration and drying. Consequently, 1.08 M 2-OPBA with a quantitative conversion of 90.2% was effectively transformed after four times fed-batch over 210 min (Fig. 4f), and the specific space–time conversion reached 30.9 mmol g−1 L−1 h−1. Comparisons with previously reported biocatalysis processes revealed that our M3–2 had the highest substrate loading and specific space–time conversion in l-HPA production (Additional file 1: Table S4). Furthermore, M3–2 displayed perfect stereoselectivity and yielded l-configuration HPA with up to 99% enantiomeric excess (ee). These results underscore the strong potential of M3–2 as a biocatalyst for continuous l-HPA synthesis.
Docking simulation and tunnel analysis of wild type and derived mutants
To better understand the enhancement in the turnover frequency (kcat) observed in enzyme-kinetics experiments at the atomic and molecular levels, we performed docking simulations using the available structural information. The substrate 2-OPBA was docked into the catalytic active center of the wild-type BbPheDH and of the derived mutants M1–3, M2–2, and M3–2 (Fig. 5). The docking simulations revealed poor adaptability of 2-OPBA to the substrate-binding pocket of wild-type BbPheDH, compared with the binding pose of the native substrate PPA (Additional file 2: Figure S8). The relatively bulky hydrophobic phenyl group side chain of 2-OPBA could not fit into the substrate-binding pocket of wild-type BbPheDH in a relatively stretched configuration (Fig. 5a, b), which may prevent the hydride transfer of the substrate in proper orientation. The low degree of freedom of the phenyl group side chain may be the reason for the low kcat value, and thereby, for the low kcat/Km value of the wild-type BbPheDH toward 2-OPBA. This observation was similar to that in a previous study . In contrast, the incorporation of three crucial mutations, V309G, L306V, and V144G, gradually enlarged the volume of the substrate-binding pocket compared with that of BbPheDH. With this change, the bulky phenyl group side chain of 2-OPBA could fit into the enlarged substrate-binding pocket with a relatively stretched configuration at a high degree of freedom (Fig. 5c–h). Hence, the mutant M3–2 showed an enhanced kcat value toward 2-OPBA compared with that of the wild-type BbPheDH. Moreover, the observed enlargement of the substrate-binding pocket can be further supported by calculating the volume of the substrate-binding pocket, which revealed that the volume of the eventual mutant M3–2 was about 196.9 Å3 larger than that of the wild-type BbPheDH (Additional file 1: Table S5). These results demonstrated that the enlargement of the substrate-binding pocket in a reasonable range may be beneficial to enhance the kcat value of enzymes toward the substrates with bulky hydrophobic side chain .
Tunnels represent potential transport pathways for small molecules, water molecules, and ions, and play a significant role in the functioning of a large variety of proteins [42, 43]. The accessibility of a tunnel depends largely on its shape, size, and amino acid composition. To a certain extent, it can be modified by protein engineering [44, 45]. Therefore, we investigated whether steric hindrance engineering alters BbPheDH tunnels. Tunnels analyses based on the structure homology model of M3–2 revealed that, compared with BbPheDH, the incorporation of V144G, L306V, and V309G formed a new tunnel (shown in magenta) in the vicinity of the catalytic active center (Fig. 6), suggesting a new transport pathway for substrate access and product egress. Furthermore, V144G, L306V, and V309G mutations markedly altered the shapes of the existing tunnels (Fig. 6). The widened tunnels (shown in red and orange) could further accelerate the transport rates of substrate and product. Both pieces of information lend support to the enhancement of the kcat value toward 2-OPBA catalyzed by M3–2 observed in the enzyme-kinetics experiments. This observation demonstrated that the formation of new tunnels and widening of existing ones can modulate the catalytic activity of enzymes. It also resembled the findings reported for prior studies [46, 47].
Here, we employed a modified steric hindrance engineering approach to create an enhanced biocatalyst with BbPheDH as a mutagenesis template for efficient l-HPA synthesis. Seven superior mutants were screened from the mutagenesis libraries constructed based on molecular docking, double-proximity filtering, and a degenerate codon. The kcat and kcat/Km values of the mutant M3–2 (V309G/L306V/V144G) were 12.7- and 12.9-fold higher than those of the wild type, respectively, it converted 1.08 M 2-OPBA to l-HPA in 210 min, and its specific space–time conversion was 30.9 mmol g−1 L−1 h−1. The substrate loading and specific space–time conversion are the highest values reported so far, demonstrating that mutant M3–2 is a useful biocatalyst with potential for the industrial production of l-HPA. Structure analysis demonstrated that mutations altered the substrate-binding pocket and tunnel, and variations in steric hindrance modulate the catalytic activity of BbPheDH towards 2-OPBA. Moreover, our study indicates the combination of molecular docking with a double-proximity filter as a promising workflow for capturing the proximal steric hindrance synergistic effects. The modified steric hindrance engineering approach can be a valuable addition and extension of the current enzyme engineering toolbox.
Strains, plasmids, and chemicals
PheDHs from Bacillus badius (BbPheDH) and Geobacillus kaustophilus (GkPheDH) were generated in the laboratory. Glucose dehydrogenase from Bacillus amyloliquefaciens (BaGluDH) was purchased from Sigma-Aldrich Corp. (Beijing, China). Escherichia coli BL21 (DE3) and plasmid pET-28a (+) were procured from Novagen (Nanjing, China) and used as the gene expression host and vector, respectively. The 2-oxo-4-phenylbutyric acid (2-OPBA) and l-homophenylalanine (l-HPA) were acquired from J&K Chemical (Beijing, China). Isopropyl-β-d-thiogalactopyranoside (IPTG), kanamycin, NADH, and NAD+ were obtained from TCI (Shanghai, China). All other chemicals were purchased from commercial sources, at least reagent grade, and used without further purification.
Sequence alignment and site-directed mutagenesis
Amino acid sequence alignment was performed using the MUSCLE server (https://www.ebi.ac.uk/Tools/msa/muscle/)  and displayed using the Esprit server (http://espript.ibcp.fr) . Site-directed mutagenesis was introduced by whole-plasmid PCR  using primers containing the required degenerate codons and the pET-28a (+) plasmid inserts the pdh gene as the template. The PCR products were analyzed by agarose gel electrophoresis and transformed into E. coli BL21 (DE3) cells. All target mutants were isolated from the constructed mixed libraries and confirmed by DNA sequencing (Talen-Bio, Wuxi, China).
Protein expression and purification
Escherichia coli BL21 (DE3) cells containing recombinant pET-28a (+) plasmids were cultured in Luria–Bertani (LB) medium with kanamycin (50 mg L−1) at 37 °C and with shaking at 200 rpm for 2 h. Enzyme expression was induced with IPTG (final concentration = 0.2 mM) when the bacteria reached OD600 = 0.6–0.8. The cells were incubated at 20 °C and with shaking at 200 rpm for 12 h. His-tagged enzymes were purified with an ÄKTA purifier system (GE Healthcare, Little Chalfont, UK). The cells were resuspended in buffer A (100 mM potassium phosphate, 150 mM NaCl, and 20 mM imidazole; pH 7.5) and disrupted with an ultrasonic cell disruption system (SCIENTZ-IID; Ningbo Scientz Biotechnology, Ningbo, Zhejiang, China) followed by centrifugation at 12,000 rpm and 4 °C for 30 min to remove cell debris and obtain cell-free extracts. The latter were passed through a 0.22-μm filter and loaded into a standard Ni–NTA affinity column (Thermo Fisher Scientific, Waltham, MA, USA) pre-equilibrated with buffer A. The column was then gradient-eluted with buffer B (100 mM potassium phosphate, 150 mM NaCl, and 500 mM imidazole; pH 7.5). Eluted fraction purity was determined by SDS-PAGE (Additional file 2: Figure S9).
Activity and kinetic parameters analysis
Reductive amination activity was determined at a 200-μL scale in 96-well microtiter plates at 30 °C. The rate of initial decrease of the absorbance was measured at 340 nm (indicating NADH consumption) using a MultiSkan GO UV-spectrometer (Thermo Fisher Scientific). The reaction mixture consisted of NH4Cl/NH4OH buffer (2 M; pH 9.5), 10 mM 2-OPBA, 0.2 mM NADH, and a certain amount of cell-free extract. One unit of enzyme activity (1 U) was defined as the amount of enzyme oxidizing 1 μM NADH per minute. Kinetic parameters with respect to 2-OPBA (PPA) were evaluated in NH4Cl/NH4OH buffer (2 M; pH 9.5) at 30 °C in the presence of 1–20 mM 2-OPBA (PPA) and 0.5 mM NADH with the purified enzymes. Kinetic parameters with respect to NADH were evaluated in NH4Cl/NH4OH buffer (2 M; pH 9.5) at 30 °C in the presence of 0.1–0.5 mM NADH and 20 mM 2-OPBA with the purified enzymes (Additional file 1: Table S6). Km and Vmax were calculated in GraphPad Prism 8 (GraphPad Software, La Jolla, CA, USA) by nonlinear curve fitting of the initial rate versus substrate concentration data to the Michaelis–Menten equation.
Thermostability assays of wild-type BbPheDH and its mutant were conducted by measuring T5030, which is the temperature at which 50% of the enzyme activity is lost after 30 min heat treatment. The enzymes were incubated at various temperatures (30–70 °C) for 30 min, and their residual activity was measured. Specific activity before incubation was normalized as 100%. Enzyme activity was measured in NH4Cl/NH4OH buffer (2 M; pH 9.5) containing 0.2 mM NADH and 10 mM 2-OPBA.
Asymmetric reductive amination of 2-OPBA to l-HPA
PheDH was coupled with GluDH for the asymmetric reductive amination of 2-OPBA. The total volume of the biocatalytic reaction mixture was 5 mL. It contained cell-free extracts of E. coli wet cells (10 g L−1 PheDH and 12 g L−1 GluDH), 2-OPBA (0.1–0.5 M), glucose (0.12–0.6 M), NAD+ (0.05–1.0 mM), and NH4OH/HCOONH4 buffer (1–5 M; pH 7.5–9.5). The reaction mixture was incubated at 25–50 °C with shaking at 200 rpm. Samples (200 μL) were drawn at intervals and alkalized with 200 μL NaOH (10 M) to terminate the reaction. The solutions were then passed through a 0.22 μm filter and subjected to High-Performance Liquid Chromatography (HPLC).
Fed-batch synthesis of l-HPA
The initial reaction mixture (100 mL) consisted of cell-free extract (10 g L−1 PheDH and 12 g L−1 GluDH), 2-OPBA (0.3 M), glucose (0.36 M), NAD+ (0.3 mM), and NH4OH/HCOONH4 buffer (3 M; pH 8.5). The reaction was initiated by adding PheDH and GluDH and incubated at 30 °C with shaking at 200 rpm. The pH value was maintained at 8.5 by titrating concentrated NH3·H2O during the reaction. Over the 30–75 min reaction period, the precipitated l-HPA was separated by filtration. The filtrate contained enzymes and cofactors and was used as the reaction medium for the next reaction period. The latter was initiated by adding 0.3 M 2-OPBA and 0.36 M glucose. Four fed-batch cycles were performed, and the 2-OPBA conversion and l-HPA enantiomeric excess (ee) values were determined by HPLC.
The conversion of 2-OPBA and the ee value of l-HPA were determined with an Agilent 1260 HPLC (Agilent Technologies, Santa Clara, CA, USA). The absolute l-HPA configuration was established by comparing it with authentic reference material after derivatization. The conditions for determining the 2-OPBA conversion rate were as follows: Diamonsil C18(2) column (5 μm; 250 mm × 4.6 mm) (Phenomenex, Torrance, CA, USA); injection volume, 10 μL; mobile phase A [55% (v/v) methanol plus 0.1% (v/v) trifluoroacetic acid], 20 min at 30 °C; flow rate, 0.8 mL min−1; and UV detection at 230 nm. The 2-OPBA retention time was 9.04 min. The ee value of l-HPA was determined using a Develosil® ODS-UG-5 column (5 μm; 150 mm × 4.6 mm) (Phenomenex) after derivatization with 1-fluoro-2,4-dinitrophenyl-5-l-alanineamide (FDAA). A 5-μL reaction sample, 4 μL of 1 M NaHCO3, and 20 μL of 1% (w/v) FDAA in acetone were mixed and heated at 40 °C for 60 min. Then 4 μL of 1 M HCl and 467 μL of 40% (v/v) aqueous acetonitrile were added to the mixture. The latter was then passed through a 0.22-μm filter and subjected to HPLC. The conditions for determining the ee value of l-HPA were as follows: injection volume, 20 μL; mobile phase B (5% (v/v) acetonitrile, 0.05% (v/v) trifluoroacetic acid, and 1% (v/v) methanol); mobile phase C (60% (v/v) acetonitrile, 0.05% (v/v) trifluoroacetic acid, and 1% (v/v) methanol); linear gradient from 0 to 100%; mobile phase C, 45 min at 30 °C; flow rate, 1.0 mL min−1; and UV detection at 340 nm. The l-HPA retention time was 33.65 min.
Homology modeling and molecular docking
Structural homology models of BbPheDH and its mutants were constructed using the SWISS-MODEL server (https://swissmodel.expasy.org/)  based on the crystal structure of the PheDH from Rhodococcus sp. M4 (PDB: 1C1D) . The 3D structure of the ligand was obtained from the PubChem website (https://pubchem.ncbi.nlm.nih.gov/) . Ligand energy was minimized with the OPTIMIZE plugin of PyMOL . AutoDock v. 4.2  was used for docking simulations. Prior to docking, each protein was protonated at pH 9.5 using the H++ server (http://biophysics.cs.vt.edu/H++)  to mimic the experimental condition. All water molecules were removed, and nonpolar hydrogen atoms were added with MGLTools v. 1.5.4. The components of the catalytic triad (K78, K90, and N276) were defined as flexible residues. A grid box of 50 × 50 × 50 Å with 0.375 Å spacing encompassed the active cavities of BbPheDH and its mutants and was set as the search space for suitable substrate-binding poses. Energies of the interactions between the ligands and the receptors were calculated with MGLTools v. 1.5.4. The lowest energy docking pose was selected for the subsequent molecular dynamics simulations.
The tunnels of BbPheDH and its mutant M3–2 were calculated with the CAVER v. 3.0 plugin of PyMOL . The probe radius was set to 1.2 Å, the default settings were applied, and the center of the catalytic triad (K78, K90, D125, and N276) was the starting point. All water molecules were removed, nonpolar hydrogen atoms were added with MGLTools v. 1.5.4, and the volumes of the substrate-binding pockets of BbPheDH and its mutants were calculated with the ProteinsPlus tool (https://proteins.plus/) . The PDB format files were uploaded for the calculations.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
High-Performance Liquid Chromatography
Koketsu K, Mitsuhashi S, Tabata K. Identification of homophenylalanine biosynthetic genes from the cyanobacterium Nostoc punctiforme PCC73102 and application to its microbial production by Escherichia coli. Appl Environ Microbiol. 2013;79:2201–8.
Liu Z, Lei D, Qiao B, Li S, Qiao J, Zhao G-R. Integrative biosynthetic gene cluster mining to optimize a metabolic pathway to efficiently produce l-homophenylalanine in Escherichia coli. ACS Synth Biol. 2020;9:2943–54.
Regulski M, Regulska K, Stanisz BJ, Murias M, Gieremek P, Wzgarda A, Niznik B. Chemistry and pharmacology of angiotensin-converting enzyme inhibitors. Curr Pharm Des. 2015;21:1764–75.
Kuhn DJ, Chen Q, Voorhees PM, Strader JS, Shenk KD, Sun CM, Derno SD, Bennett MK, Van Leeuwen FWB, Chanan-Khan AA, Orlowski RZ. Potent activity of carfilzomib, a novel, irreversible inhibitor of the ubiquitin-proteasome pathway, against preclinical models of multiple myelorna. Blood. 2007;110:3281–90.
Alanazi F, Kwa FAA, Burchall G, Jackson DE. New generation drugs for treatment of multiple myeloma. Drug Discov Today. 2020;25:367–79.
Liese A, Seelbach K, Wandrey C. Industrial biotransformations. 2nd ed. Weinheim: Wiley-VCH Verlag GmbH; 2007. p. 94–396.
Yu LT, Huang JL, Chang CY, Yang TK. Formal synthesis of the ACE inhibitor benazepril·HCl via an asymmetric Aza-Michael reaction. Molecules. 2006;11:641–8.
Xie Y, Lou R, Zhi L, Mi A, Jiang Y. DPAMPP in catalytic asymmetric reactions: enantioselective synthesis of l-homophenylalanine. Tetrahedron Asymmetry. 2000;11:1487–94.
Sheldon RA, Woodley JM. Role of biocatalysis in sustainable chemistry. Chem Rev. 2018;118:801–38.
Bornscheuer UT, Huisman GW, Kazlauskas RJ, Lutz S, Moore JC, Robins K. Engineering the third wave of biocatalysis. Nature. 2012;485:185–94.
Larissegger-Schnell B, Glueck S, Kroutil W, Faber K. Enantio-complementary deracemization of (±)-2-hydroxy-4-phenylbutanoic acid and (±)-3-phenyllactic acid using lipase-catalyzed kinetic resolution combined with biocatalytic racemization. Tetrahedron. 2006;62:2912–6.
Rodríguez-Alonso MJ, Clemente-Jiménez JM, Rodrígueš Vico F, Vázquez F. Rational re-design of the “double-racemase hydantoinase process” for optically pure production of natural and non-natural l-amino acids. Biochem Eng J. 2015;101:68–76.
Lo HH, Kao CH, Lee DS, Yang TK, Hsu WH. Enantioselective synthesis of (S)-2-amino-4-phenylbutanoic acid by the hydantoinase method. Chirality. 2010;15:699–702.
Joon-Young H, Jihyang P, Joo-Hyun S, Minho C, Byung-Kwan C, Juhan K, Byung-Gee K. Simultaneous synthesis of 2-phenylethanol and l-homophenylalanine using aromatic transaminase with yeast Ehrlich pathway. Biotechnol Bioeng. 2009;102:1323–9.
Luo C, Lin Q, Lin S, Meng C, Wang H. Cosynthesis of l-homophenylalanine and 2-phenylethanol by recombinant Saccharomyces cerevisiae expressing aspartate aminotransferase from Escherichia coli BL21(DE3). J Biosci Bioeng. 2017;123:1–7.
Ahmad AL, Oh PC, Shukor SRA. Sustainable biocatalytic synthesis of l-homophenylalanine as pharmaceutical drug precursor. Biotechnol Adv. 2009;27:286–96.
Asano Y, Yamada A, Kato Y, Yamaguchi K, Hibino Y, Hirai K, Kondo K. Enantioselective synthesis of (S)-amino acids by phenylalanine dehydrogenase from Bacillus sphaericus: use of aatural and recombinant enzymes. J Org Chem. 1990;55:5567–71.
Bradshaw CW, Wong CH, Hummel W, Kula MR. Enzyme-catalyzed asymmetric synthesis of (S)-2-amino-4-phenylbutanoic acid and (R)-2-hydroxy-4-phenylbutanoic acid. Bioorg Chem. 1991;19:29–39.
Ahmad AL, Low EM, Abd Shukor SR. Immobilization of phenylalanine dehydrogenase onto Eupergit CM for the synthesis of (S)-2-amino-4-phenylbutyric acid. J Mol Catal B Enzym. 2013;88:26–31.
Zhang J, Zhu T, Wu X, Chen Y. Enhancement of biocatalytic efficiency by increasing substrate loading: enzymatic preparation of l-homophenylalanine. Appl Microbiol Biotechnol. 2013;97:8487–94.
Luetz S, Giver L, Lalonde J. Engineered enzymes for chemical production. Biotechnol Bioeng. 2008;101:647–53.
Currin A, Swainston N, Day PJ, Kell DB. Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently. Chem Soc Rev. 2015;44:1172–239.
Reetz MT. Laboratory evolution of stereoselective enzymes: a prolific source of catalysts for asymmetric reactions. Angew Chem Int Ed. 2011;50:138–74.
Nie Y, Wang S, Xu Y, Luo S, Zhao YL, Xiao R, Montelione G, Hunt JF, Szyperski T. Enzyme engineering based on X-ray structures and kinetic profiling of substrate libraries: alcohol dehydrogenases for stereospecific synthesis of a broad range of chiral alcohols. ACS Catal. 2018;8:5145–52.
Wu K, Yang Z, Meng X, Chen R, Shao L. Engineering an alcohol dehydrogenase with enhanced activity and stereoselectivity toward diaryl Ketones: reduction of steric hindrance and change of stereocontrol element. Catal Sci Technol. 2020;10:1050–60.
Yamamoto K, Miyake H, Kusunoki M, Osaki S. Steric hindrance by 2 amino acid residues determines the substrate specificity of isomaltase from Saccharomyces cerevisiae. J Biosci Bioeng. 2011;112:545–50.
Tian K, Li Z. A simple biosystem for the high-yielding cascade conversion of racemic alcohols to enantiopure amines. Angew Chem Int Ed. 2020;59:21745–51.
Franklin RD, Mount CJ, Bommarius BR, Bommarius AS. Separate sets of mutations enhance activity and substrate scope of amine dehydrogenase. ChemCatChem. 2020;12:2436–9.
Lewis JC, Mantovani SM, Fu Y, Snow CD, Komor RS, Wong CH, Arnold FH. Combinatorial alanine substitution enables rapid optimization of cytochrome P450BM3 for selective hydroxylation of large substrates. ChemBioChem. 2010;11:2502–5.
Ye LJ, Toh HH, Yang Y, Adams JP, Snajdrova R, Li Z. Engineering of amine dehydrogenase for asymmetric reductive amination of ketone by evolving Rhodococcus phenylalanine dehydrogenase. ACS Catal. 2015;5:1119–22.
Savile CK, Janey JM, Mundorff EC, Moore JC, Tarn S, Jarvis WR, Colbeck JC, Krebber A, Fleitz FJ, Brands J. Biocatalytic asymmetric synthesis of chiral amines from ketones applied to sitagliptin manufacture. Science. 2010;329:305–9.
Chen FF, Zheng GW, Liu L, Li H, Chen Q, Li FL, Li CX, Xu JH. Reshaping the active pocket of amine dehydrogenases for asymmetric synthesis of bulky aliphatic amines. ACS Catal. 2018;8:2622–8.
Zheng Y, Li F, Lin Z, Lin G, Hong R, Yu H, Xu J. Structure-guided tuning of a hydroxynitrile lyase to accept rigid pharmaco aldehydes. ACS Catal. 2020;10:5757–63.
Brunhuber NMW, Thoden JB, Blanchard JS, Vanhooke JL. Rhodococcus l-phenylalanine dehydrogenase: kinetics, mechanism, and structural basis for catalytic specifity. Biochemistry. 2000;39:9174–87.
Yousefi F, Ataei F, Arab SS, Hosseinkhani S. Increase of Bacillus badius phenylalanine dehydrogenase specificity towards phenylalanine substrate by site-directed mutagenesis. Arch Biochem Biophys. 2017;635:44–51.
Paradisi F, Conway PA, Maguire AR, Engel PC. Engineered dehydrogenase biocatalysts for non-natural amino acids: efficient isolation of the d-enantiomer from racemic mixtures. Org Biomol Chem. 2008;6:3611–5.
Busca P, Paradisi F, Moynihan E, Maguire AR, Engel PC. Enantioselective synthesis of non-natural amino acids using phenylalanine dehydrogenases modified by site-directed mutagenesis. Org Biomol Chem. 2004;2:2684–91.
Miyazaki K, Takenouchi M. Creating random mutagenesis libraries using megaprimer PCR of whole plasmid. Biotechniques. 2002;33:1033–8.
Zhang J, Tao S, Zhang B, Wu X, Chen Y. Microparticle-based strategy for controlled release of substrate for the biocatalytic preparation of l-homophenylalanine. ACS Catal. 2014;4:1584–7.
Zhou F, Mu X, Nie Y, Xu Y. Enhanced catalytic efficiency and coenzyme affinity of leucine dehydrogenase by comprehensive screening strategy for l-tert-leucine synthesis. Appl Microbiol Biotechnol. 2021;105:3625–34.
Yin XJ, Liu YY, Meng LJ, Zhou HS, Wu JP, Yang LR. Rational molecular engineering of glutamate dehydrogenases for enhancing asymmetric reductive amination of bulky α-keto acids. Adv Synth Catal. 2019;361:803–12.
Gouaux E, Mackinnon R. Principles of selective Ion transport in channels and pumps. Science. 2005;310:1461–5.
Jiang Y, Lee A, Chen J, Cadene M, Chait BT, Mackinnon R. Crystal structure and mechanism of a calcium-gated potassium channel. Nature. 2002;417:515–22.
Klvana M, Pavlova M, Koudelakova T, Chaloupkova R, Dvorak P, Prokop Z, Stsiapanava A, Kuty M, Kuta-Smatanova I, Dohnalek J. Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations. J Mol Biol. 2009;392:1339–56.
Pavlova M, Klvana M, Prokop Z, Chaloupkova R, Banas P, Otyepka M, Wade RC, Tsuda M, Nagata Y, Damborsky J. Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate. Nat Chem Biol. 2009;5:727–33.
Kre N, Halder JM, Rapp LR, Hauer B. Unlocked potential of dynamic elements in protein structures: channels and loops. Curr Opin Chem Biol. 2018;47:109–16.
Heinemann PM, Armbruster D, Hauer B. Active-site loop variations adjust activity and selectivity of the cumene dioxygenase. Nat Commun. 2021;12:1–12.
Madeira F, Park YM, Lee J, Buso N, Gur T, Madhusoodanan N, Basutkar P, Tivey ARN, Potter SC, Finn RD, Lopez R. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 2019;47:W636–41.
Robert X, Gouet P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014;42:W320–4.
Andrew W, Martino B, Stefan B, Gabriel S, Gerardo T, Rafal G, Heer FT, Christine R, Lorenza B. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296–303.
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S, Li Q, Shoemaker BA, Thiessen PA, Yu B, Zaslavsky L, Zhang J, Bolton EE. PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res. 2021;49:D1388–95.
O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open Babel: an open chemical toolbox. J Cheminform. 2011;3:1–14.
Morris GM, Huey R, Lindstrom W, Sanner MF, Olson AJ. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J Comput Chem. 2009;30:2785–91.
Gordon JC, Myers JB, Folta T, Shoja V, Heath LS, Onufriev A. H++: a server for estimating PK(a)s and adding missing hydrogens to macromolecules. Nucleic Acids Res. 2005;33:W368–71.
Chovancova E, Pavelka A, Benes P, Strnad O, Brezovsky J, Kozlikova B, Gora A, Sustr V, Klvana M, Medek P, Biedermannova L, Sochor J, Damborsky J. CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures. PLoS Comput Biol. 2012;8:23–30.
Schoening-Stierand K, Diedrich K, Faehrrolfes R, Flachsenberg F, Meyder A, Nittinger E, Steinegger R, Rarey M. ProteinsPlus: interactive analysis of protein-ligand binding interfaces. Nucleic Acids Res. 2020;48:W48–53.
This work was supported by the National Key Research and Development Program of China (Grant Numbers 2021YFC2100100), the National Natural Science Foundation of China (NSFC) (Grant Numbers 21336009 and 21176103), the 111 Project (Grant Number 111-2-06), the High-end Foreign Experts Recruitment Program (Grant Number G20190010083), and the National First-Class Discipline Program of Light Industry Technology and Engineering (Grant Number LITE2018-09).
Ethics approval and consent to participate
Consent for publication
All authors approved to publish to article.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers used for site-directed mutagenesis. Table S2. Enzyme thermostabilities of BbPheDH and its superior mutants. Table S3. Effects of substrate concentration on the asymmetric reductive amination of 2-OPBA catalyzed by mutant M3–2. Table S4. Comparison between M3–2 and other reported PheDHs for the synthesis of l-HPA from 2-OPBA. Table S5. Calculated substrate-binding pocket volume and area of BbPheDH and its superior mutants. Table S6. Kinetic parameters of BbPheDH and its superior mutants toward NADH.
Time course of asymmetric reductive amination of 2-OPBA catalyzed by BbPheDH and GkPheDH. Figure S2. Amino acid sequence alignment of the PheDHs from different sources. Figure S3. Relative activity of the single-site mutants constructed in the first round of steric hindrance engineering. Figure S4. Relative activity of the double-site mutants constructed in the second round of steric hindrance engineering. Figure S5. Relative activity of the triple-site mutants constructed in the third round of steric hindrance engineering. Figure S6. Relative activity of the quadruple-site mutants constructed in the fourth round of steric hindrance engineering. Figure S7. Relative activity of BbPheDH and GluDH in different concentrations of NH4OH/HCOONH4 buffer (pH 8.5). Figure S8. Binding poses of native substrate PPA (a) and bulky substrate 2-OPBA (b) with BbPheDH. Figure S9. SDS-PAGE analysis of BbPheDH and its superior mutants.
About this article
Cite this article
Wu, T., Mu, X., Xue, Y. et al. Structure-guided steric hindrance engineering of Bacillus badius phenylalanine dehydrogenase for efficient l-homophenylalanine synthesis. Biotechnol Biofuels 14, 207 (2021). https://doi.org/10.1186/s13068-021-02055-0
- Phenylalanine dehydrogenase
- Steric hindrance
- Enzyme engineering
- Catalytic efficiency
- Reductive amination