Kai
Wu
a,
Jinrong
Yan
bd,
Qinde
Liu
ac,
Xiaojing
Wang
a,
Piaoru
Wu
a,
Yiyang
Cao
a,
Xiuhong
Lu
a,
Yixin
Xu
a,
Junhai
Huang
*bd and
Lei
Shao
*ad
aSchool of Pharmacy, Shanghai University of Medicine & Health Sciences, 279 Zhouzhu Highway, Pudong New Area, Shanghai 201318, China. E-mail: shaolei00@gmail.com
bCollege of Chemistry and Chemical Engineering, Shanghai University of Engineering Science, 333 Longteng Road, Shanghai 201620, China. E-mail: huang_junhai@163.com
cShanghai University of Traditional Chinese Medicine, 1200 Cailun Road, Shanghai 201203, China
dState Key Laboratory of New Drug and Pharmaceutical Process, Shanghai Institute of Pharmaceutical Industry, 285 Gebaini Rd., Shanghai 200040, China
First published on 15th December 2023
Imine reductases (IREDs) are important biocatalysts in the asymmetric synthesis of chiral amines. However, a detailed understanding of the stereocontrol mechanism of IRED remains incomplete, making the design of IRED for producing the desired amine enantiomers challenging. In this study, we investigated the stereoselective catalytic mechanism and designed an (R)-stereoselective IRED from Paenibacillus mucilaginosus (PmIR) using pharmaceutically relevant 2-aryl-substituted pyrrolines as substrates. A putative mechanism for controlling stereoselectivity was proposed based on the crucial role of electrostatic interactions in controlling iminium cation orientation and employed to achieve complete inversion of stereoselectivity in PmIR using computational design. The variant PmIR-Re (Q138M/P140M/Y187E/Q190A/D250M/R251N) exhibited opposite (S)-stereoselectivity, with >96% enantiomeric excess (ee) towards tested 2-aryl-substituted pyrrolines. Computational tools were employed to identify stabilizing mutations at the interface between the two subunits. The variant PmIR-6P (P140A/Q190S/R251N/Q217E/A257R/T277M) showed a nearly 5-fold increase in activity and a 12 °C increase in melting temperature. The PmIR-6P successfully produced (R)-2-(2,5-difluorophenyl)-pyrrolidine, a key chiral pharmaceutical intermediate, at a concentration of 400 mM with an ee exceeding 99%. This study provides insight into the stereocontrol elements of IREDs and demonstrates the potential of computational design for tailored stereoselectivity and thermal stability.
The understanding of the specific roles of varying active-site residues in stereocontrol remains limited. To gain insight into the correlation between active site residues and stereoselectivity, multiple studies have compared IREDs exhibiting stereocomplementary activity towards 2-methyl-1-pyrroline and other cyclic imines.6–8 These investigations have led to the categorization of IREDs into two groups based on their opposite stereoselectivities and the putative catalytic residues involved. Y-type IREDs possess tyrosine as a potential proton-donating amino acid, whereas D-type IREDs feature aspartate at equivalent position. Although it seemed that the different putative catalytic residues, tyrosine and aspartate, led to opposite stereoselectivity of the IREDs, there remain exceptions. In 2015, Wetzl and coworkers reported IREDs that possessed phenylalanine or asparagine residues at equivalent positions.9 In 2016, Turner et al. reported an intriguing finding with AoIRED, which demonstrated higher activity and stereoselectivity in asymmetric imine reduction despite the absence of tyrosine or aspartate residues.10 These discoveries suggest that potentially different stereocontrol mechanisms are formed by distinct active sites in IREDs. The varying active-site residues and unidentified key stereocontrol elements make the rational design of an IRED for the desired stereoselectivity challenging. Additionally, stability has been emphasized as a crucial factor in the practical application of evolved IREDs.2,11,12 The conventional approach for mass screening in directed evolution is time-consuming. However, a rational design approach for improving stereoselectivity and stability remains underdeveloped.
In this study, we investigated the possible stereocontrol mechanism of IRED in Paenibacillus mucilaginosus (PmIR). Furthermore, by employing rational design with catalytic constraints based on a putative stereocontrol mechanism, we successfully reversed the stereoselectivity of PmIR towards 2-(2,5-difluorophenyl)pyrroline (2-DFPL) and other 2-aryl-substituted pyrrolines in the PmIR-Re variant. PmIR-Re represents the first instance of overcoming the boundary between Y-type- and D-type IREDs, provides valuable insights into the key stereocontrol elements of IREDs, and demonstrates an effective mechanism-guided approach for manipulating stereoselectivity. Furthermore, a combinational method including salt-bridge construction and ΔΔGfold prediction was applied to enable the rapid identification of interface mutations for stability enhancement. Finally, the engineered variant PmIR-6P was used to successfully produce (R)-2-(2,5-difluorophenyl)-pyrrolidine (2-DFPD), a key chiral intermediate of larotrectinib, demonstrating its potential for practical applications.
After screening IREDs stored in our laboratory, PmIR from Paenibacillus mucilaginosus emerged as one of the few enzymes capable of catalyzing the reduction of 2-DFPL to yield (R)-2-DFPD, which displayed the highest ee, exceeding 96%. Owing to the conserved binding position of NADPH in IREDs,5,7,13 the PmIR-NADPH complex model was constructed by superimposing the PmIR crystal structure (PDB ID: 8KFK) on the NADPH-bound BcSIRED crystal structure (sequence identity 75.8%). To investigate the key residues of PmIR controlling the asymmetric reduction, 18 residues around Y187 within 10 Å, located at the subunit interface forming the binding pocket, were targeted for alanine scanning (Fig. S1 and Table S2†). Next, key residues with potential catalytic roles that dramatically affected the activity or stereoselectivity were subjected to directed mutagenesis to further investigate their role in asymmetric reduction (Table 1). Alanine scanning results indicated that the M194A, F194A, F232A, and F234A mutations led to negligible activity. Both D250A and W195A mutations resulted in a complete loss of activity. The alanine mutation of the potential proton donor Y187, which is often considered catalytic and conserved in IREDs, displayed a low yet detectable conversion with an almost racemic ee of 4%, suggesting the possible role of a bulky aromatic side chain in activity and stereoselectivity. To further elucidate the contribution of the phenolic hydroxyl group of Y187 to imine reduction, the Y187F mutant was constructed. The Y187F mutant exhibited a conversion of 34%, slightly lower than that observed for the WT (46.5%), but stereoselectivity increased from 96% to >99%. Enzymatic kinetic analysis of the Y187F mutant revealed a 3.6-fold increase in KM but only a 24% reduction in kcat compared with those in the WT (Table 2). The predicted pKa of 2-DFPL is 6.60 ± 0.20 (CAS Scifinder), suggesting that at least 75% of the substrate exists in the form of the iminium cation in the reaction pH of 6.0. In summary, the catalytic activity of the Y187F mutant and its decreased substrate affinity suggest that Y187 in PmIR primarily contributes to substrate binding during catalysis, rather than serving as a proton donor, which differs from the previously reported role of tyrosine at equivalent positions in other IREDs. This indicates that Y187 is a modifiable residue that can be used for further enzyme design.
Variants | Conversiona (%) | ee (%) | Stereoselectivity |
---|---|---|---|
a The reaction mixture (500 μL) included 275 μL crude enzyme extract from 50 mg mL−1 resting cells, 50 μL of 2.5 mM NADP stock, 50 mM glucose, 10 mM substrate, 100 μL crude enzyme of recombinant glucose dehydrogenase, and potassium phosphate buffer (100 mM, pH 6.0). The reaction mixture was stirred at 30 °C for 40 min and analyzed using chiral HPLC. b Not detected. | |||
Wild type | 46.5 ± 3.6 | 96 | R |
Y187A | 4.3 ± 1.2 | 4 | R |
Y187F | 34.8 ± 0.9 | >99 | R |
Y187D | 2.3 ± 0.9 | 58 | R |
Y187E | 4.5 ± 0.8 | 33 | S |
D250A | n.d.b | n.d. | n.d. |
D250E | 38.1 ± 4.4 | 97 | R |
Variants | K M (mM) | V max (μmol min−1 mg−1) | k cat (s−1) | k cat/KM (s−1 mM−1) |
---|---|---|---|---|
WT | 0.84 | 0.78 | 0.44 | 0.52 |
Y187F | 3.09 | 0.67 | 0.38 | 0.12 |
D250E | 0.68 | 0.72 | 0.41 | 0.60 |
Aspartate is another residue often considered catalytically active in IREDs. A previous report of a fungal IRED catalyzing imine reduction and reductive amination indicated the important role of the Asp active site in binding an imine cation and amine substrate.4 A systematic evaluation of imine-reducing enzymes by Pleiss et al. suggested that most IREDs require a negative electrostatic potential for substrate-binding.8 Residue D250 of PmIR was positioned above Y187, located in the binding pocket, and close to NADPH. A complete loss of activity was observed in the D250A mutant with no detectable products. Saturation mutagenesis was performed to investigate the role of D250. Most mutations resulted in complete inactivation, except for D250E, which featured a carboxyl group with one additional methyl moiety compared with aspartate. The D250E mutant exhibited a comparable ee value, kcat, and slightly reduced KM values (Table 2). The optimal reaction pH for PmIR was determined to be 6.0, corresponding to the deprotonated state of the acidic residues. Therefore, the complete inactivation of the alanine mutation suggests that the loss of activity is primarily due to the absence of a negatively charged deprotonated carboxyl group. Furthermore, using the previously reported (R)-selective BcSIRED17 catalyzing 2-DFPL with 89% conversion and 88% ee, negligible conversion (1.7%) and relatively high conversion (19%) were observed with BcSIRED mutants D251A and Y188F, respectively. These results were consistent with the findings observed for PmIR. Multiple sequence alignments of eight IREDs displaying (R)-stereoselectivity towards 2-aryl-substituted pyrrolines also revealed that D250 is conserved in certain IRED subgroups (Fig. 1A). Electrostatic interactions between the negatively charged D250 and iminium cation during catalysis cannot be excluded in PmIR.
![]() | ||
Fig. 1 (A) Multisequence alignment of reported IREDs with (R)-stereoselectivity towards 2-aryl-substituted pyrrolines. The IREDs used for sequence alignment include PmIR (this study), ScIR (WP_003961038.1),15 IRED-GF3546 (BAM99301),16BcsIRED (WP_000841637),17SrIRED (WP_012890722),7PeIRED (WP_010497949.1),18SaIRED (WP_016644884.1),16 and NhIRED (WP_017538739.1).17 The figure was prepared using ESPript19 (numbering as PmIR). Y187 and D250 were conserved in the (R)-selective enzyme. (B) Representative snapshots showing the reorientation of the D250 side chain. (C) Docking conformation in pro-R pose and following conformation analysis with molecular dynamics (MD) simulations. The distance distributions indicate the high proportion of reactive conformation in the pro-R pose. |
Several studies speculated that conformational changes and domain motion of IREDs can affect stereoselectivity.9,10,17 The crystal structure alone may have limitations in that it cannot reflect important motions. Although the crystal structure indicated that the carboxyl group of D250 was not ideally positioned to interact with the iminium cation, MD simulations revealed that D250 was in a flexible binding loop (residues 242–251) (Fig. S3†), and conformational reorientation of the D250 side chain was observed during 100 ns simulations of substrate-free PmIR (Fig. 1B). This reorientation suggests that the carboxyl group could dynamically adjusts to accommodate the positive iminium cation, stabilizing it in a reactive conformation and promoting hydride transfer. Hence, D250 was selected as a flexible residue for molecular docking. The docking results indicate preferential binding of the substrate in the pro-R pose with electrostatic interaction to the D250 side chain (Fig. 1C), as a lower interaction energy of pro-R pose than pro-S pose (−31.9 vs. −27.2 kcal mol−1) was observed. The Y187, F232, F194, and W195 form a hydrophobic region that stabilizes 2-DFPL aromatic ring binding, which could explain why alanine mutations at these sites resulted in inactive or low activities. These interactions could position the specific CN face at an appropriate distance from NADPH.
Based on the importance of electrostatic interactions and hydride transfer during reduction, the distances between the imine nitrogen atom and the D250 carboxyl carbon atom (d1), and the distance between the imine carbon and NADPH-C4 (d2) are critical for catalysis. Electrostatic interactions can be strong, even at relatively long distance.20,21 The ternary complex of AoIRED has been previously reported, which revealed a binding site for the amine product that placed the chiral carbon within 4.0 Å of the C4 atom of NADPH.9 In addition, studies of enzymes such as short-chain dehydrogenase, which catalyzes bio-reduction using NADPH as a cofactor, reported that hydride transfer could happen when the distance between the NADPH-C4 atom and electrophile carbon is ≤ 4.5 Å,22–24 and a few of them show moderate imine-reducing activity.8 Therefore, conformations of d1 less than 5.0 Å and d2 less than 4.0 Å were considered reactive. Based on these two distance criteria, MD simulations of the docking complex demonstrated the high occurrence of a reactive conformation in the pro-R pose of the protonated 2-DFPL (Fig. 1C). These findings suggest that negatively charged acidic residues play a role in promoting the formation of reactive pro-R positions during catalysis. Therefore, based on the results of site-directed mutagenesis and MD simulations, we inferred that D250 is a crucial residue for stereocontrol in PmIR.
Based on the above findings, we propose a putative mechanism for stereocontrol, in which the electrostatic interaction between the negatively charged side chains of Asp or Glu and the iminium cation governs the orientation of the pyrroline face, thereby determining the stereoconfiguration of the product. To manipulate the stereoselectivity, we hypothesized that the orientation of the 1-pyrroline face could be reversed by simply introducing new electrostatic interactions in the opposite direction, considering that a proton transfer system or any flanked residues that lower the pKa value of the proton donor are unnecessary. Because the residues at positions 250 and 187 are adjacent to each other and close to the nicotinamide moiety of NADPH, we speculated that replacing Y187 with an acidic amino acid could create an electrostatic interaction in the opposite direction, thereby altering the orientation of the 1-pyrroline face. The side chains of aspartate and glutamate in the model compounds, with the titratable residue in the middle and two alanines on each side, are determined experimentally by previous study in water, yielding pKa values of 3.6 and 4.2, respectively.25 The protonation state of these residues was taken into consideration to ensure that most of the side chain was negatively charged in the protein. Using a pKa predicting module in Discovery Studio 4.0,26 the pKa value of Y187E was predicted as 5.4 with a 1.2 shift compared with glutamate in water, and Y187D was predicted as 4.7 with a 1.1 shift compared with aspartate in water, indicating a negative electrostatic potential at pH 6.0. As anticipated, the Y187D mutation resulted in a dramatically reduced ee to 58%, but a negligible conversion. The Y187E mutation successfully resulted in slightly reversed (S)-stereoselectivity with an ee of 33% (Table 1), providing further evidence for the crucial role of electrostatic interactions in stereoselective control. The docking results were in accordance with the proposed stereoselective control mechanism, as a pro-S pose was observed in the Y187E variant complex using the 2-DFPL iminium as the substrate, showing an opposite orientation to the 1-pyrroline face (Fig. 2).
Encouraged by the putative stereoselective control mechanism and 33% ee observed for the Y187E variant, we used a complex with the pro-S pose in the reactive conformation from MD simulations as the initial structure for the redesign to enhance (S)-stereoselectivity and activity simultaneously. To achieve the distances, angles and dihedrals between Y187E and substrate and NADPH for the catalytic geometric constraints, the cluster model with 2-DFPL in the pro-S pose was optimized to determine the transition state (TS) structure (Fig. S4 and Table S3†). We focused on redesigning the loop residues surrounding 2-DFPL, including V249, D250, and R251 in loops 242–251, and Q138 and P140 in loops 135–140. Additionally, we included Q190 in the design because Q190A displayed a higher conversion than the WT in previous alanine scanning. The hydrophobic region formed by F232, F194, and W195 that stabilizes 2-DFPL aromatic ring binding was not selected for design. A virtual library with 30000 mutant structures was generated using the Rosetta Design to explore the sequence and conformation space and optimize the catalytic interactions. These structures were then sort based on the Rosetta parameters of Total Score, Constraint Energy, and Interface Energy.
The top four non-redundant variants, 5352, 6645, 25658, and PmIR-Re, were further analyzed to determine the occurrence of reactive pro-S conformations using MD simulations. Experimental validations were conducted for 5352 and PmIR-Re, which showed a relatively high probability of adopting the pro-S pose in the reactive conformation. Both variants exhibited significant improvements in stereoselectivity, producing (S)-2-DFPD in 93% ee and >99% ee, respectively. Interestingly, the D250S and D250M mutations, which resulted in a complete loss of activity in the WT, were present in the two active variants, suggesting that E187 replaced D250 to anchor the iminium cation. Importantly, this observation further supports the hypothesis that the inactivation caused by D250S and D250M in the WT can be attributed to the absence of electrostatic interactions rather than destabilization of the protein structure.
The initial complex used for design demonstrated that the lack of tyrosine bulky side chain and electrostatic interaction with Y187E resulted in a lower binding position of the pro-S pose (Fig. S5†). In contrast, the side-chain repacking of Q190A and Q138M in PmIR-Re elevated 2-DFPL binding position, allowing for better interaction with the hydrophobic and aromatic regions formed by F232, F194, and W195. The P140M side chain could also push substrate into the aromatic region and the D250M could form a π–sulphur interaction with 2-DFPL from above. In addition, the PmIR-Re complex displayed a geometric feature that more closely resembled the TS structure than the initial design template, potentially enhancing the hydride transfer process (Table S4†). This observation indicates that the increased S-stereoselectivity and activity achieved by PmIR-Re could be due to the more reactive conformation for hydride transfer while retaining the electrostatic interaction with Y187E. These results further support the idea that the stereoselectivity of PmIR can be directly manipulated by creating electrostatic interactions in opposite directions, thereby reversing the orientation of the 1-pyrroline face.
To the best of our knowledge, only a few studies have investigated the stereoselective control mechanisms of imine reduction reactions. Several studies have compared IREDs that display stereo-complementary activity towards 2-methyl-1-pyrroline and have categorized them into Y-type and D-type IREDs based on their stereoselectivity and putative catalytic residues (tyrosine or aspartate).6,7 In a recent study by Pleiss et al., the stereoselectivity of a D-type IRED from Myxococcus stipitatus towards 2-methyl-1-pyrroline was reversed.29 They obtained a variant with the opposite stereoselectivity and 91% ee through site-saturation mutagenesis. In this case, the stereoselectivity switch was attributed to a difference in the conformation of the nicotinamide moiety of NADPH rather than the substrate-binding conformation. More recently, Himo emphasized the importance of steric clashes in the transition-state energy of the pro-R and pro-S complexes of AoIRED.30 The pro-S conformation, which lacked steric repulsion, resulted in a lower energy state and led to (S)-stereoselectivity towards bulky imine substrates. These results suggest that varying the active sites of IREDs offers multiple possibilities for fine-tuning the relative spatial positions of the imine and nicotinamide groups of the cofactors.
In this study, we crossed the boundary between Y-type and D-type IREDs by redesigning Y-type PmIR into a D/E-type enzyme with opposite stereoselectivity. PmIR-Re represents an engineered IRED with completely reversed stereoselectivity (>99% ee) achieved through mechanism-guided computational design. In contrast to the previously reported residues W210 in AspRedAm31 and W234/F260 in IRED-M5,32 which potentially affect stereoselectivity by steric hindrance of the bulky side chain, the stereoselectivity in PmIR could be controlled by electrostatic interactions between the iminium cation and the acidic amino acids adjacent to the nicotinamide moiety of NADPH. Therefore, the stereoselectivity of IREDs can be tailored through rational reconstruction of electrostatic interactions and reshaping of the binding pocket.
The PmIR structure is composed of reciprocal homodimers consisting of an N-terminal Rossman fold domain and a C-terminal helical domain connected by a long interdomain helix. The two subunits are held together by interacting helices, and it is noteworthy that the residues from both subunits contribute to the formation of the binding pocket (Fig. S1†). Natural salt bridges were also observed in IREDs on the surface between the two subunits, providing stability to the dimeric structure.38 Based on these observations, we hypothesized that enhanced subunit interactions would result in more stable variants. Herein, we employed a comprehensive approach to obtain a stable variant, PmIR-6P (P140A/Q190S/R251N/Q217E/A257R/T277M), which involved identifying interface mutations through alanine scanning, designing salt bridges between interfaces, and identifying mutations using ΔΔGfold calculations (Table 3).
Mutation | Origin | Location | T m (°C) | ΔTm (°C) | Relative activitya |
---|---|---|---|---|---|
a The activity of PmIR-6P was determined at 35 °C and others were determined at 30 °C. | |||||
WT | — | — | 52.1 | 0 | 100 ± 0.8 |
P140A/Q190S/R251N | Alanine scan | Interface | 57.2 | 5.1 | 450.2 ± 2.9 |
P140A/Q190S/R251N/A257R | Salt bridge | Interface | 59.0 | 1.8 | 533.0 ± 2.5 |
P140A/Q190S/R251N/Q217E | Salt bridge | Interface | 58.7 | 1.5 | 418.5 ± 1.8 |
P140A/Q190S/R251N/G143M | ΔΔGfold prediction | Surface | 58.5 | 1.3 | 283.6 ± 2.1 |
P140A/Q190S/R251N/T277M | ΔΔGfold prediction | Surface | 59.3 | 2.1 | 390.3 ± 3.8 |
PmIR-6P (P140A/Q190S/R251N/Q217E/A257R/T277M) | Combination of mutations | 64.1 | 12.0 | 453.9 ± 4.6 |
The previous alanine scanning revealed that the interface residues P140, Q190, and R251 demonstrated higher conversions than the WT (Table S2†). Saturation mutagenesis was conducted to identify the most beneficial mutations at each site. Additive effects of multiple mutations were also investigated (Table S5†). As a result, variant P140A/Q190S/R251N showed a 4.5-fold increase in catalytic efficiency with ee exceeding 99% compared to the WT. The purified variant also exhibited a Tm of 57.2 °C, which was 5.1 °C higher than that of the WT. These findings are consistent with previous reports on evolved IREDs, where mutations in the binding pocket not only increased activity, but also improved stability.29,31 These dual effects may be attributed to enhanced subunit interactions, resulting in a more compact active pocket that stabilizes substrate binding.
Based on the geometric criteria for salt bridges (Fig. S5†) and the crystal structure of PmIR, potential residue pairs that met the requirements for salt bridge formation were identified. Residues with B-factors above the average value were selected as candidates for designing salt bridges. A257 and Q217 were mutated to A257R and Q217E, respectively, which introduced a salt bridge between R257 and D87, and between E217 and R176. The two new salt bridges increased Tm by 1.8 °C and 1.5 °C, respectively.
In silico saturation mutagenesis was performed for each residue of PmIR, and the ΔΔGfold values were calculated using FoldX and Rosetta_ddg to generate point mutation structures. To narrow down the large number of variants for experimental validation, we performed MD simulations and visual inspection on variants below cut-off. The variants exhibited reduced ΔΔGfold values, lower local flexibility than the original residues (ΔZ-score <0), and no obvious structural undesirable mutations were selected for experimental validation (Table S6†). Among the 12 single mutations tested, four variants (G143M, G143W, T277M, and T277L) showed higher residual activity than the template after incubation at 60 °C. G143M and T277M, which had higher residual activity, were further purified to determine their Tm values (Table 3). This in silico workflow allowed for the rapid identification of mutations with a higher likelihood of enzyme stabilization and has been successful in several biocatalysts.39–41 However, there is a risk that certain mutations may increase stability at the expense of reduced activity. Although the G143M mutant exhibited a higher Tm, it also showed a significant reduction in activity. Therefore, it was not included in the final variant. Finally, mutations identified from alanine scanning, salt bridge design, and ΔΔGfold prediction were combined, resulting in the highly stable variant PmIR-6P, which exhibited a 12 °C increase in Tm, from 52.1 to 64.1 °C (Fig. S9†).
The highest specific activity of PmIR-6P was observed at an optimal reaction temperature of 35 °C (Fig. S10A†). At this temperature, the kinetic parameters KM and kcat were determined to be 0.92 mM and 2.01 s−1, respectively, showing a 4.7-fold increase in catalytic efficiency compared with the WT. The half-lives of both the WT and PmIR-6P variants were measured to assess their kinetic stability. At 35 °C, PmIR-6P exhibited an 8.4-fold improvement in half-life compared to the WT, increasing from 47.5 h to 396.0 h. Furthermore, PmIR-6P exhibited a half-life of 280.6 h at 40 °C (Fig. S10B†).
MD simulations revealed a significant difference in the radius of gyration and root mean square fluctuation (RMSF) between the WT and PmIR-6P at 308 K (35 °C). PmIR-6P exhibited a more compact structure, which was attributed to mutations at the subunit interface (Fig. 4A). The RMSF of WT and PmIR-6P was calculated and compared (Fig. S11†). It indicated that the most region of PmIR-6P with a less flexibility during simulation compared with WT. Interestingly, the 242–251 binding loop of PmIR -6P was observed with a higher mobility than WT, which suggested the loop motion might affect activity for PmIR.
The distance between the backbone Cα atoms of Q217-R176 (Fig. 4B) and D87-R257 (Fig. 4C) clearly indicate that the salt bridges stretched the two subunits more tightly than in the WT. These findings suggest that subunit dissociation plays a crucial role in the irreversible inactivation of IREDs. This characteristic has also been observed in ω-transaminases, where stabilizing mutations are predominantly located at the interface of the two subunits as the active pocket is formed there.41 Therefore, focusing on the subunit interfaces of IREDs may be effective for identifying stabilizing mutations.
Previously reported engineered IREDs were mostly stabilized through directed evolution involving extensive screening. For example, 256 of 296 amino acid positions of IR-46 from Saccharothrix espanaensis were subjected to single-site saturation mutagenesis and screening before the best variant was identified.36 Similarly, SpRedAm from Streptomyces purpureus was engineered to produce cis-cyclobutyl-N-methylamine, with 93 sites (out of 296) selected for single-site saturation mutagenesis and screening.34 Additionally, in the case of ScIR from Streptomyces clavuligerus, significantly increased activity and stabilization were achieved by screening 52 sites.35 These examples demonstrated the power of directed evolution in IRED engineering. However, a rational design approach that requires significantly less screening than the conventional directed evolution still needs to be developed for IREDs. In this study, we combined different computational tools to enhance the thermal stability of IREDs and significantly reduce the library size for experimental testing. Furthermore, our results suggested that the subunit interface, including the binding pocket, plays a crucial role in stabilizing IREDs.
Footnote |
† Electronic supplementary information (ESI) available: Materials and methods; supplementary tables and figures; HPLC spectra; 1H and 13C NMR spectra. See DOI: https://doi.org/10.1039/d3sc04636b |
This journal is © The Royal Society of Chemistry 2024 |