Front. Chem. Frontiers in Chemistry Front. Chem. 2296-2646 Frontiers Media S.A. 10.3389/fchem.2021.664156 Chemistry Original Research Improving the Thermostability and Activity of Transaminase From Aspergillus terreus by Charge-Charge Interaction Cao Jia-Ren 1 † Fan Fang-Fang 1 † Lv Chang-Jiang 1 † Wang Hong-Peng 1 Li Ye 1 Hu Sheng 2 Zhao Wei-Rui 2 Chen Hai-Bin 3 Huang Jun 1 * Mei Le-He 2 4 5 * 1School of Biological and Chemical Engineering, Zhejiang University of Science and Technology, Hangzhou, China 2School of Biotechnology and Chemical Engineering, NingboTech University, Ningbo, China 3Enzymaster (Ningbo) Bio-Engineering Co., Ltd., Ningbo, China 4Jinhua Advanced Research Institute, Jinhua, China 5Department of Chemical and Biological Engineering, Zhejiang University, Hangzhou, China

Edited by: Gao-Wei Zheng, East China University of Science and Technology, China

Reviewed by: Aitao Li, Hubei University, China; Feng Cheng, Zhejiang University of Technology, China

 *Correspondence: Jun Huang huangjun@zust.edu.cn Le-He Mei meilh@zju.edu.cn

This article was submitted to Chemical Biology, a section of the journal Frontiers in Chemistry

†These authors have contributed equally to this work

 14 04 2021 2021 9 664156 05 02 2021 23 03 2021 Copyright © 2021 Cao, Fan, Lv, Wang, Li, Hu, Zhao, Chen, Huang and Mei. 2021 Cao, Fan, Lv, Wang, Li, Hu, Zhao, Chen, Huang and Mei

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Transaminases that promote the amination of ketones into amines are an emerging class of biocatalysts for preparing a series of drugs and their intermediates. One of the main limitations of (R)-selective amine transaminase from Aspergillus terreus (At-ATA) is its weak thermostability, with a half-life (t1/2) of only 6.9 min at 40°C. To improve its thermostability, four important residue sites (E133, D224, E253, and E262) located on the surface of At-ATA were identified using the enzyme thermal stability system (ETSS). Subsequently, 13 mutants (E133A, E133H, E133K, E133R, E133Q, D224A, D224H, D224K, D224R, E253A, E253H, E253K, and E262A) were constructed by site-directed mutagenesis according to the principle of turning the residues into opposite charged ones. Among them, three substitutions, E133Q, D224K, and E253A, displayed higher thermal stability than the wild-type enzyme. Molecular dynamics simulations indicated that these three mutations limited the random vibration amplitude in the two α-helix regions of 130–135 and 148–158, thereby increasing the rigidity of the protein. Compared to the wild-type, the best mutant, D224K, showed improved thermostability with a 4.23-fold increase in t1/2 at 40°C, and 6.08°C increase in T5010. Exploring the three-dimensional structure of D224K at the atomic level, three strong hydrogen bonds were added to form a special “claw structure” of the α-helix 8, and the residues located at 151–156 also stabilized the α-helix 9 by interacting with each other alternately.

 amine transaminase thermostability enzyme thermal stability system site-directed mutagenesis molecular dynamics simulations Introduction

Chiral amines are important components of many significant bioactive compounds, pharmaceutical intermediates and agrochemical industry products (Bornscheuer et al., 2012; Mathew and Yun, 2012a; Ghislieri and Turner, 2014; Park et al., 2014; Fuchs et al., 2015; Ferrandi and Monti, 2017; Dawood et al., 2018; Cai et al., 2020). In addition to optical rotation, enantiomers of chiral drugs have the same physical properties, but they are absorbed, activated or degraded by the metabolic system of the human body in different ways during the pharmacological action, resulting in different efficacies and toxicities (Burke and Henderson, 2002). Because of their broad-spectrum biological activities and high purity in the synthesis of enantiomeric amines, they can act as chiral building blocks for the synthesis of more complex structural drugs (Breuer et al., 2004; Nugent, 2010). For instance, aromatic chiral amine derivatives are intermediates of the highly potent KCNQ2 opener (Wu et al., 2004), rivastigmine agents to treat Alzheimer's disease (Hua et al., 2018), and anti-arthritic drugs (Dyckman et al., 2011). Unfortunately, synthetic routes for these compounds are still challenging.

Many enzymes have been employed for the synthesis of chiral amines, including transaminases, imine reductases, amine dehydrogenases, and reductive aminases (Turner and Truppo, 2010; Łyskowski et al., 2014; Godwin et al., 2017; Gomm and O'Reilly, 2018; Jiang et al., 2020). Among them, the coenzyme PLP of transaminase can be recycled, unlike many expensive coenzymes, which gradually deplete the amino products synthesized in situ (Paul et al., 2014). In particular, amine transaminases (ATAs) have many industrial advantages and chemical properties over conventional chemical synthesis of optically pure chiral amines, including excellent stereoselectivity (Svedendahl et al., 2010) and broad substrate spectrum (Guo and Berglund, 2016). More importantly, ATAs can achieve continuous flow biotransformations under mild conditions (Andrade et al., 2014) as an alternative technology to replace toxic, non-recyclable chemical catalysts and reduce the use of high- temperature and high-pressure conditions in the chemical production process.

Transaminases have been extensively studied during the past few years for the synthesis of chemically pure chiral amines (Hhne and Bornscheuer, 2009; Zhu and Hua, 2009; Hailes et al., 2010; Nugent, 2010; Tufvesson et al., 2011; Malik et al., 2012; Mathew and Yun, 2012b; Kroutil et al., 2013). However, transaminases usually require better reaction rates, higher temperature adaptability in industrial production, and reduce risk of microbial contamination. Protein engineering plays a vital role in enhancing the thermal stability of (R)-selective At-ATA to expand its applicability in industrial processes (Liu et al., 2019). To date, the rational design of protein engineering involves many factors, such as surface electrostatic interactions, hydrophobic interactions, B-factor values, consensus mutagenesis, disulphide bridges, coevolution networks, and hydrogen bonding interactions (Pace et al., 2011; Wang et al., 2014; Zhang et al., 2014, 2020; Huang et al., 2017; Xie et al., 2018, 2019; Moon et al., 2019; Zhu et al., 2019; Cao et al., 2020). All of these have been employed to develop stable proteins.

Recently, the enzyme thermal stability system (ETSS), a suite of computational programs based on TK-SA model calculation and surface charge-charge interaction analysis was released (Zhang et al., 2014). The TK model was constructed by Tanford and Kirkwood in 1957 (Tanford and Kirkwood, 1957), based on the charge-charge interaction to describe the electrostatic properties of the whole protein (Matthew et al., 1985; Matthew and Gurd, 1986a,b). After solvent accessibility (SA) was introduced to refine the TK-SA model, Bashford and Karplus majorized the TK-SA model with the effects of partition function (Sanchez-Ruiz et al., 1999) and Gibbs free energy (Matthew et al., 1979; Richmond, 1984; Bashford and Karplus, 1991), which has been successfully applied in protein modification engineering (Elcock, 2001; Sanchez-Ruiz and Makhatadze, 2001; Ibarra-Molero and Sanchez-Ruiz, 2002; Makhatadze et al., 2003, 2004; Strickler et al., 2006; Gribenko et al., 2009; Schweiker and Makhatadze, 2009).

In this work, the potential mutation sites were replaced with electrically neutral amino acids on the protein surface based on the TK-SA model, which was used to construct the mutant with enhanced thermostability. Subsequently, the prospective stabilizing effects of these mutations were verified by thermal inactivation experiments, and basic-amino-acid scanning was more accurate in finding the most thermal-stability-improved mutation at each site. In addition, these variants at each site were analyzed by molecular dynamics (MD) simulation (Purmonen et al., 2007), with the aim of exploring the improved thermostability and catalytic activity of the mutants at the atomic level.

 Materials and Methods Materials

The At-ATA cDNA from Aspergillus terreus sequence, including the NcoI and XhoI restriction sites, was synthesized by General Biosystems (Chuzhou, China), and the plasmid pET-28a(+) was used for gene cloning and DNA sequencing. All PCR primers were synthesized by Qingke Biology Co., Ltd. (Hangzhou, China). PrimeSTAR® Max DNA polymerase was obtained from Takara Biotechnology (Dalian, China) for the polymerase chain reaction (PCR). Dpn I, Yeast extract and tryptone were obtained from Thermo Fisher Scientific (Shanghai, China). Dimethyl sulfoxide (DMSO), 1-(R)-PEA and pyruvate were obtained from Aladdin Biochemical Technology Co., Ltd. (Shanghai, China). NaCl, NaH2PO4, Na2HPO4, NaOH, DNA ladder, protein marker, protein loading buffer, kanamycin sulfate, isopropyl-β-d-thiogalactoside (IPTG), Ni-NTA Sefinose (TM) Resin (Settled Resin) kit, SDS-PAGE gel kit, and Modified Bradford Protein Assay Kit were obtained from Sangon (Shanghai, China). E. coli BL21(DE3) Chemically Competent Cell, EasyPure® HiPure Plasmid MaxiPrep Kit, EasyPure® Quick Gel Extraction Kit and EasyPure® PCR Purification Kit were purchased from TransGen Biotech (Beijing, China).

 Location of the Mutant Sites

Based on the crystal structure of At-ATA (PDB ID: 4CE5) obtained from the Protein Data Bank (http://www.rcsb.org), ETSS was used to calculate the total interaction energy between charged amino acids at points i and j (Eij) from wild-type At-ATA. Based on these results, we selected the modification of residues with positive Eij values, which signified an unfavorable interaction, to enhance the At-ATA thermostability. In addition, three principals were used to redesign At-ATA mutants: (i) residues with high positive Eij values were priorities, (ii) residues predicted to form hydrogen bonds with the nearby conserved in At-ATA, and (iii) residues far from the catalytic and binding pocket were selected to retain the activity. When negatively charged amino acids are involved in the mutation residues, it is a good strategy to turn the target to be amidated to approach electric neutrality. When both negatively and positively charged residues located at the same position were failed to depress the Eij values or turned to the opposite ones, the charged residues will attempt to be alanine. All three-dimensional (3D) structures of At-ATA were visualized using PyMOL 2.0.7 software (http://pymol.org).

 Site-Directed Mutagenesis

The At-ATA gene from A. terreus was cloned into pET-28a(+), using Nco I and Xho I, and transformed into E. coli DH5a. All primers used for site-directed mutagenesis are listed in Supplementary Table 1. For the mutagenic PCR, first stage: 98°C for 1 min, one cycle; second stage: 98°C, 15 s/55°C, 15 s/72°C, 2 min, 30 cycles; and third stage: 72°C for 10 min, one cycle; 2 × PrimeSTAR® Max DNA polymerase (1 ×, 25 μL), forward and reverse primers (10 μM, 1 μL each) together with template (pET-28a(+)-At-ATA; 0.1 ng/μL, 2 μL), and diluted with ddH2O to 50 μL. The PCR products were purified using an EasyPure® PCR Purification Kit. Following depuration, Dpn I (20 U) was added to the buffer, and the mixture was incubated for 3 h at 37°C and transformed into E. coli BL21 (DE3) bacteria by heat shock. After screening on a Luria-Bertani (LB) plate containing 50 μg/mL kanamycin and incubated at 37°C overnight, the sequence of the mutated plasmid DNA was aligned and checked.

 Protein Expression and Purification

The colony was designed by transferring the recombinant plasmid into E. coli BL21 (DE3) using the EasyPure® HiPure Plasmid MaxiPrep Kit. Protein was expressed by adding a colony to LB-Kan medium (5 mL; 50 μg/mL), and the strain was grown at 37°C and 200 rpm for 6 h, as the OD600 reached ~0.6. The culture was then added to LB-Kan medium (200 mL; 50 μg/mL) and allowed to continue growing at 37°C and 180 rpm for another 3 h. Thereafter, the protein expression was induced by adding IPTG to a final concentration of 1 mM and returned to 25°C, 150 rpm, for 20 h. Subsequently, the cells were centrifuged at 4°C (6,000 × g, 6 min), then washed twice with buffer A (50 mM sodium phosphate buffer, and pH 8.0). The cells were later dissolved in 55 mL of buffer B (300 mM NaCl, 50 mM sodium phosphate buffer, 20 mM imidazole, and pH 8.0) and disrupted using a high-pressure homogenizer (ATS, Jiangsu, People's Republic of China) for 1 min in 140 mPa until the mixture became clear. After the protein solution was centrifuged at 8,000 rpm for 55 min at 4°C, the purified protein, containing an N-terminal His6-tag, was attained using a Ni-NTA Sefinose column and eluted with buffer C (300 mM NaCl, 50 mM sodium phosphate buffer, 250 mM imidazole, and pH 8.0). Furthermore, SDS-PAGE (12% separating and 5% stacking gels) and a Modified Bradford Protein Assay Kit (Sangon Biotech Co., Ltd. Shanghai, China) were used to analyze purified proteins. Protein concentrations were determined by the Bradford method using BSA as a standard.

 Thermostability of the <italic>At</italic>-ATA and Variants

Both the half-life (t1/2) and half-life temperature (T5010) values characterize the thermostability of the enzyme. t1/2 is defined as the time when the residual activity of At-ATA and its mutants was reduced to 50% of its original activity at 40°C. Similarly, T5010 refers to the temperature at which the enzyme activity was reduced to half of the original activity after heat treatment at a continuous temperature for 10 min.

The purified At-ATA and its mutants were incubated for 0–30 min at 40°C, and then cooled on ice for 10 min. Enzyme activity test was performed at 25°C for 3 min. An exponential function model: Exp2PMod1 [formula: y = exp(-kd·t)] by nonlinear regression was used to fit our data using Origin 8.0, then the first-order rate constants (kd) and 50% of relative enzyme activity were determined. In addition, the enzyme solution was incubated for 10 min at temperatures of 4, 25, 35, 37, 40, 42, 45, 50, and 55°C, and then cooled on ice for 10 min. The data were fitted to a four-parameter Boltzmann sigmoidal function reformed with the Levenberg-Marquardt iterative algorithm. The formula is presented in Equation (1).

R=A+B-A1+e(Tm-T)C

Where R is the percentage of residual activity at temperature T, A, and B are the pretransitional and posttransitional percentages of residual activity, respectively, and C is the slope factor. In the absence of treatment, Tm = T0.

Differential scanning fluorimetry (DSF) is a rapid and highly efficient method for identifying protein thermostability (Niesen et al., 2007). The protein unfolding temperature was gauged by increased fluorescence of the dye molecules. The dye molecules have affinity for hydrophobic portions of the protein, exposed as the protein unfolds. The mixture consisted of 1 mg/mL pure enzyme, 1× SYPRO Orange dye (dissolved in DMSO) and diluted with buffer C (150 mM NaCl, 50 mM sodium phosphate buffer, and pH 8.0) to 50 μL. The sample with buffer C instead of pure enzyme was used as the negative control. The measurements were performed on a StepOne Real-Time PCR System (Applied Biosystems, USA). The temperature from 25 to 70°C was scanned in 0.7°C increments, with each temperature maintained for 30 s. The excitation and emission wavelengths were 490 and 605 nm, respectively. The melting temperature Tm was calculated using the formula shown in Equation (2):

y=UF+ NF-UF1+eTm-x  α

Where UF and NF are the minimum and maximum emission fluorescence intensities, respectively, and α is the slope of the curve within Tm.

 Enzyme Activity Assay and Kinetic Parameters

The activities of the wild-type and mutant strains were measured as described by Schätzle et al. (2009). One unit of activity was defined as the amount of enzyme required to release 1 μmol of acetophenone per min under the assay conditions. The substrate pre-mixture was prepared containing 2.5 mM 1-(R)-PEA, and 2.5 mM pyruvate, 0.1 mM PLP and 0.25% (w/v) DMSO in 180 μL of buffer A, and 20 μL of enzyme concentrated with buffer A. The reaction was measured in UV 96-well microtiter plates at 25°C and 245 nm for 3 min, and the production of acetophenone was monitored using an MD 190 photometer (Molecular Devices, Sunnyvale, CA, USA). The kinetic parameters of At-ATA for 1-(R)-PEA and pyruvate were determined by measuring the activities at different substrate concentrations when either R-MBA or pyruvate was a fixed concentration of 2.5 mM, and the other substrate was changed at 0-3.0 mM until substrate inhibition was observed. The kinetic results were fitted to the Michaelis-Menten equation in Origin 8.0.

 Molecular Dynamics Simulation

We redesigned the 3D structure of At-ATA combined with PLP based on the crystal structure of At-ATA (PDB ID: 4CE5) from the Protein Data Bank (http://www.rcsb.org). In view of the pretreatment of enzyme without substrate, we had to delete the free Lys180 and PLP-amino donor compound (PDG), and combined PLP with Lys180. According to the experimental conditions, the protonation states of all ionizable residues were assigned based on pKa values from the PROPKA software in combination with visual inspection of local hydrogen bonding networks using Discovery Studio 2018 at pH 8.0, including the protonation status of ionizable residues and lysine affected by the pocket environment.

MD simulation was performed at a constant temperature (313 K) for 50 ns using the Amber 14 force field of YASARA (version 16.4.6) software (http://www.yasara.org). The 3D structures were filled with water with a density of 0.998 mg/L and inserted into a cube with edge lengths of 10 Å. Sodium and chloride ions (0.9%) were added as counter ions to form an electrically neutral system, and the ionizable groups were protonated according to their pKa values at pH 8.0 in the medium. The systems were optimized by three-step energy minimization at the molecular mechanics (MM) level to adjust the poor interatomic interactions. First, the water molecules were minimized while keeping the protein and substrate constrained. Then, the side chains were allowed to relax while the main protein chains were restrained. Finally, the entire system was completely relaxed without any restrictions. After each energy minimization was completed, a 1000-step conjugate gradient iteration loop was performed. Subsequently, the optimized system was gradually heated from 0 to 313 K in a constant volume environment for 150 ps, and then balanced for 150 ps in a constant pressure environment with the density of the system gradually becoming 0.997 g/cm3. Finally, a 50 ns MD simulation with a step length of 2.5 fs was completed under constant pressure conditions, and the trajectory was collected every 25 ps. The cutoff value of the van der Waals force and electrostatic interaction during the simulation was handled at 8.0 Å. Analyses of protein structures including root mean square deviation (RMSD) of backbone atom positions, and root mean square fluctuation (RMSF) for individual residues were performed using YASARA (Purmonen et al., 2007; Dong et al., 2018).

 Results Selection of Mutants With Increased Thermal Stability <italic>in silico</italic>

To predict the residue with an optimized value and analyze the flexibility of the surrounding amino acid residues, ETSS was used to evaluate the interaction of charged amino acids (Zhang et al., 2014). The program was used to calculate the interaction parameters of At-ATA, and the total Eij shown in Figure 1 reveals that there are 95 charged amino acids in the monomer of At-ATA. Thus, we selected amino acids with positive values and far from the active center to mutate to amino acids with opposite or neutral charges. Therefore, four high-value residues located in different loop regions from the surface of the protein were selected (E133, D224, E253, and E262 shown in Supplementary Figure 1).

Total energy of each chargeable residue. Positive value indicates the overall contribution of repulsive force to the thermal stability of protein structure, whereas negative value represents the overall contributions from gravitation, which is beneficial to the thermostability of the protein structure.

 Thermostability Analysis of the Mutant <italic>At</italic>-ATA by Alanine Scanning

Wild-type was scanned for alanine by site-directed mutagenesis and successfully expressed in E. coli BL21(DE3). As shown in Supplementary Figure 2, the four purified variants (E133A, D224A, E253A, and E262A) showed a single band with an apparent molecular mass of 35-45 kDa as the wild-type enzyme (36.1 kDa). Except for the mutant E262A, other mutants showed an increased thermostability and activity comparable to that of the wild-type enzyme (data shown in Figure 2A), and its t1/2 and T5010 values were measured. Among them, the D224A mutant showed a significant 3.05-fold increase in t1/2 at 40°C and withstood a higher temperature value (T5010) of 43.1°C, compared to the wild-type enzyme with a t1/2 of 6.9 min and T5010 of 38.5°C (Figures 2B,C and Supplementary Table 2). The results indicate that eliminating the charge of acidic residues on the At-ATA surface may improve the thermostability.

The thermal stability of wild-type At-ATA and its mutants by the alanine scanning: (A) the enzymatic activity of wild-type At-ATA and the four variants, (B) the thermal deactivation of wild-type At-ATA and the four variants at various temperatures for 10 min (T5010), (C) the thermal deactivation half-life (t1/2) of wild-type At-ATA and the three variants at 40°C. Black, Wild-type; red, E133A; blue, D224A; and purple, E253A.

 Thermostability Analysis of the Mutant <italic>At</italic>-ATAs by Basic-Amino-Acid Scanning

Suitable basic amino acids were selected as targets, and nine mutants (E133H, E133K, E133R, E133Q, D224H, D224K, D224R, E253H, and E253K) were obtained to analyze the change rule of thermostability affected by the electrostatic interaction combined with the results of alanine scanning. The new mutant consistently showed higher activity at 40°C. Their stabilities were then compared at the preference temperatures. As shown in Figure 3 and Table 1, the T5010 values for the E133Q, D224K, and E253A mutants, which is the best one in each mutation group, were 41.12, 44.59, and 40.69°C, respectively, compared with that of the wild-type (38.5°C). The three mutants had consistently higher residual activity at 40°C, which may keep half of the activity 1.71-fold, 2.17-fold, and 4.24-fold higher than the wild-type. The melting profiles of the four enzymes were determined by monitoring the fluorescence of SYPRO Orange dye using DSF over the temperature range of 25–70°C, as a method for enzyme thermostability (Zhu et al., 2019). A melting temperature (Tm) of 47.6 ± 0.2°C was obtained with the fitting data for D224K (Table 1 and Supplementary Figure 3), with R2 > 0.99 in both cases. The data indicate the consistency of thermodynamic stability.

The thermal stability of wild-type At-ATA and its mutants in each site. (A–C) the thermal deactivation of wild-type At-ATA and the 12 variants at various temperatures for 10 min (T5010), (D–F) the thermal deactivation half-life (t1/2) of wild-type At-ATA and the 12 variants at 40°C.

The stability of the wild-type and stabilized mutant At-ATAs.

Name T5010 (°C) t1/2 (min) Tm (°C)
WT 38.5 ± 0.5 6.9 ± 0.6 41.3 ± 0.2
E133Q 41.1 ± 0.2 15.0 ± 0.5 42.6 ± 0.2
D224K 44.5 ± 0.2 29.2 ± 0.2 47.6 ± 0.1
E253A 40.7 ± 0.5 11.8 ± 0.6 42.4 ± 0.2
Kinetics of <italic>At</italic>-ATA and Its Mutants

The At-ATA fit the single-substrate Michaelis-Menten kinetics at 25°C when pyruvate or 1-(R)-PEA was used as the single variable. Kinetic parameters were measured at 25°C by carrying out 3 min reactions over a range of pyruvate or 1-(R)-PEA from 0.125 to 3.0 mM giving the kinetic values shown in Table 2. The kcatpyruvate values for E133Q, D224K, and E253A were 2.72-fold, 2.60-fold, and 2.17-fold higher than the wild-type, and for 1-(R)-PEA, the kcat1-(R)-PEA values for the E133Q, D224K, and E253A were 1.82-fold, 1.01-fold, and 2.20-fold higher than that of the WT, respectively, showing lower reaction energy requirements for both substrates. The Kmpyruvate for three mutants were 0.50, 0.35, and 0.69 mM, which were 2.17-fold, 1.52-fold, and 3.00-fold higher than that of the wild-type. In contrast, the Km1-(R)-PEA values for the three mutants were decreased to 0.16, 0.14, and 0.17 mM compared with that of the wild-type (0.23 mM). Thus, the kcat/Km values of the E133Q, D224K, and E253A mutants were 2.74, 3.76, and 2.51 L/(s·mmol) for 1-(R)-PEA, and 7.46, 4.64, and 8.38 L/(s·mmol) for pyruvate, respectively.

Steady-state kinetic constants of wild-type and stabilized mutant At-ATAs.

Name kcatpyruvate Kmpyruvate kcat/Kmpyruvate kcat1−(R)−PEA Km1−(R)−PEA kcat/Km1−(R)−PEA
(s−1) (mM) (L/(s·mmol)) (s−1) (mM) (L/(s·mmol))
WT 0.50 ± 0.01 0.23 ± 0.02 2.22 0.64 ± 0.01 0.23 ± 0.03 2.82
E133Q 1.36 ± 0.03 0.50 ± 0.01 2.74 1.17 ± 0.01 0.16 ± 0.01 7.46
D224K 1.30 ± 0.02 0.35 ± 0.02 3.76 0.65 ± 0.02 0.14 ± 0.01 4.64
E253A 1.73 ± 0.02 0.69 ± 0.01 2.51 1.41 ± 0.01 0.17 ± 0.02 8.38
Conformation and Energy Change Revealed by Molecular Dynamics Simulation

The stability structure of At-ATA and its mutants at 313 K was established by MD simulation and further analyzed. Analysis of the residual level fluctuations demonstrated that the wild-type enzyme is more flexible than mutants at the 130-135 site (belonging to α-helix 8) and 148-158 site (belonging to α-helix 9) (Figures 4a,b and Supplementary Figures 4A,B). In particular, the largest fluctuations at residue Arg131 in the D224K model was 2.81 nm which was 1.09 less than that of the wild-type as simulated using the RMSF. The reason for this increased stability is likely to be the powerful hydrogen bond formed by Ile135 and Arg131, Pro132, and Glu13 from the α-helix 8, which are 2.6, 2.8, and 3.4 Å in length, respectively (Figure 4c). More importantly, these strong hydrogen bonds reduce the distance between α-helix 8 and α-helix 6 as shown in Figure 5, which shorten from 12.7 to 6.5 Å in the D224K model. At the same time, three new hydrogen bond interactions were also added at α-helix 9, as shown in Figure 4d (Ala 276 and Val 149; Pro 152 and Met 154; Pro152 and Gln155 shown in red). In addition, the distance of two hydrogen bonds (the NH in Met154 and the O in Glu151, the NH in Arg156 and the O in Pro152, shown in blue) were reduced to 1.9 and 2.3 Å, respectively, compared with that of the wild-type. These strong hydrogen bonds reduce the flexibility of α-helix 8 and α-helix 9 and greatly stabilize the structure of the protein in the D224K mutant. The RMSF shown in Figures 4a,b reveals that both E133Q and E253A mutants are less flexible than the wild-type, although not as stable as D224K, especially the region at α-helix 8 and α-helix 9. Supplementary Figures 4C–F show four newly added (red) and four shortened (blue) hydrogen bond interactions acting on α-helix 8 and α-helix 9 in the E133Q mutant. For the E253A mutant, five newly added (red) and two shortened (blue) hydrogen bond interactions directly affected the activity of α-helix 8 and α-helix 9, reducing the radius of the helix (Supplementary Figures 4G,H).

MD analysis of At-ATA and D224K using YASARA at 313 K in the last 20 ns. (a) The RMSF values of At-ATA and D224K; (b) the detail RMSF values of α-helix 8 and α-helix 9; (c,d) 3D view of two α-helix structures. The hydrogen bonds were displayed by dotted line with red and blue.

3D view of the surrounding structure of α-helix 8 surrounding on the wild type (blue) and D224K (green).

 Discussion

Transaminase is a biocatalyst which can efficiently catalyze the synthesis of chiral amines. It has become an important research topic to improve the thermal stability of transaminase by site-directed mutagenesis to expand its application in industry. At present, different protein engineering strategies have been used to improve the thermal stability of transaminases. Among them, some good protein engineering strategies can reliably predict mutants with improved stability, and are easy to implement (Jones et al., 2017). In this work, the results showed that the combination of computer-aided rational design and amino acid scanning proved to be promising for finding mutation candidates with enhanced thermal stability.

In this study, four key residues of At-ATA, E133, D244, E253, and E262, were identified to redesign by site-directed mutagenesis combined with the open-source software ETSS. Among them, the alanine scanning results showed that three mutants (E133A, D224A, and E253A) could enhance the thermostability of transaminase. Based on the above results, we screened the most suitable charged amino acids for each site. Finally, E133Q, D224K, and E253A were confirmed to have positive effects on the stability of transaminase. Viewed in PyMOL 2.0.7, Glu133 is located at helix 5, whereas Asp224 belongs to the loop area connected to β-sheet 12 and β-sheet 13, and Glu253 is located at a loop area between helix 8 and β-sheet 12.

To further study the results of the experiment, a virtual mutation model was built by Discovery Studio 2018, and recalculated by YASARA and ETSS. Intermolecular interactions, such as hydrogen bonds, hydrophobic interactions, van der Waals forces, ionic bonds and disulphide bonds play a decisive role in stabilizing the tertiary structure of proteins (Fan et al., 2018; Xie et al., 2019). MD simulation results also support the key role of substitution of positively charged amino acids (glutamine, lysine, histidine, and arginine) or neutral amino acids (alanine) in thermal stability. As shown in Figure 4c, a special “claw structure” was found in D224K mutant, in which Ile135 combined with Arg131, Pro132, and Glu133 three residues. These three hydrogen bonds were formed compared with the wild-type, resulting in the great stability of α-helix 8. The “claw” was connected between α-helix 8 and α-helix 6, thus it would be more inward and closer to the center of the mutant than the wild-type (Figure 5). In addition, the O atom in the main chain of Pro152 formed a hydrogen bond with the N atom in the main chain of Met154, Gln155, and Arg156, which led to the Pro152 firmly catching hold of the center of α-helix 9. Further, a hydrogen bond is formed between Val149 and Ala276 with 2.3 Å, thus approaching the β-sheet 15 and 16 with Ala276 to avoid large fluctuation (Figure 4d). For E133Q, the interactions between Thr146 and Glu155 shortened the distance between α-helix 9 and β-sheet 7 and 16, respectively (Supplementary Figures 4C–F). Similar to the structure of D224K, the strategy of alternating hydrogen bonding reduces the fluctuation range of α-helix 9 in E253A (Supplementary Figures 4G,H). In addition, based on the ETSS analysis, when Asp224 was mutated to Lys224, the Eij value was decreased from 22.16 to 17.37 kJ/mol in the D224K mutants. For E133Q and E253A mutant, the Eij values were decreased by 2.17 and 2.28 kJ/mol, respectively, compared with that of the wild-type. This indicates that the entire energy contribution is transferred from an unfavorable to favorable state, which improves the thermal stability significantly. Furthermore, it was found that the effect of oppositely charged residues on the charged amino acids at different positions was different.

 Conclusion

In this study, we obtained three stable mutants, E133Q, D224K, and E253A, according to the ETSS strategy based on the TK-SA algorithm. The mutant D224K not only improved the thermal stability, but also increased the enzyme activity to a certain extent. At the same time, we analyzed hydrogen bond interactions in the two helix regions (α-helix 8 and α-helix 9) by MD simulation, and these were significantly different from the wild type. This study explored the important influence of charge interaction on the structure of transaminase, and provided a feasible strategy for improving the thermal stability of transaminase and promoting its application in industrial production.

 Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

 Author Contributions

JH and L-HM designed research and assigned personnel. J-RC and C-JL conducted experiments. F-FF performed the molecular dynamics simulation. H-PW analyzed the data. H-BC provided the technical help. J-RC, C-JL, and F-FF prepared the original draft of the manuscript. YL, SH, and W-RZ reviewed the manuscript and gave suggestions. All authors read and approved the manuscript.

 Conflict of Interest

H-BC was employed by the company Enzymaster (Ningbo) Bio-Engineering. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We would like to thank Associate Professor Lujia Zhang (East China Normal University, Shanghai, China) for providing ETSS and helpful technical assistance.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.664156/full#supplementary-material

The 3D structure of At-ATA with putative important residues introduced by ETSS and the four residues with strong charge-charge interactions were visualized in the crystal structure of At-ATA (PDB ID: 4CE5) using the program PyMOL (http://pymol.org).

The SDS-PAGE analysis of wild-type and At-ATA mutants. Lanes: 1, Protein marker; 2, Wild-type; 3, E133A; 4, D224A; 5, E253A; 6, E262A.

Thermal unfolding of wild-type At-ATA and its mutants were monitored by DSF.

MD analysis of At-ATA and E133Q and E253A using YASARA at 313 K in the last 20 ns. (A) The RMSF values of At-ATA and its mutants. (B) The detail RMSF values of α-helix 8 and α-helix. (C–F) The hydrogen bond interaction in two α-helix for the mutant E133Q; (G,H) the hydrogen bond interaction in two α-helix for the mutant E253A; The added bonds were displayed by red, and the blue one showed a shortened change compared with the wild type.

Primers used for site-directed mutagenesis.

The stability of the wild-type and stabilized mutant At-ATAs.

 References Andrade L. H. Kroutil W. Jamison T. F. (2014). Continuous flow synthesis of chiral amines in organic solvents: immobilization of E. coli cells containing both ω-transaminase and PLP. Org. Lett. 16, 6092–6095. 10.1021/ol502712v Bashford D. Karplus M. (1991). Multiple-site titrationcurves of proteins: an analysis of exact and approximate methods for their calculation. J. Phys. Chem. 95, 9556–9561. 10.1021/j100176a093 Bornscheuer U. T. Huisman G. W. Kazlauskas R. J. Lutz S. Moore J. C. Robins K. (2012). Engineering the third wave of biocatalysis. Nature 485, 185–194. 10.1038/nature11117 Breuer M. Ditrich K. Habicher T. Hauer B.;, Kesseler, M. Strmer R. Zelinski T. (2004). Industrial methods for the production of optically active intermediates. Angew. Chem. 43, 788–824. 10.1002/anie.20030059914767950 Burke D. Henderson D. J. (2002). Chirality: a blueprint for the future. Br. J. Anaesth. 88, 563–576. 10.1093/bja/88.4.56312066734 Cai R.-F. Liu L. Chen F.-F.;, Li, A.-T. Xu J.-H. Zheng G.-W. (2020). Reductive amination of biobased levulinic acid to unnatural chiral γ-amino acid using an engineered amine dehydrogenase. ACS Sustain. Chem. Eng. 8, 17054–17061. 10.1021/acssuschemeng.0c04647 Cao J.-R. Li F.-W. Li Y.-Y.;, Chen, H.-P. Liao X.-J. Zhang Y. (2020). Hydrophobic interaction driving the binding of soybean protein isolate and chlorophyll: improvements to the thermal stability of chlorophyll. Food Hydrocoll. 113:106465. 10.1016/j.foodhyd.2020.106465 Dawood A. W. H. deSouza R. O. M. A. Bornscheuer U. T. (2018). Asymmetric synthesis of chiral halogenated amines using amine transaminases. ChemCatChem. 10, 951–955. 10.1002/cctc.201701962 Dong Y.-W. Liao M.-L. Meng X.-L. Somero G. N. (2018). Structural flexibility and protein adaptation to temperature: molecular dynamics analysis of malate dehydrogenases of marine molluscs. Proc. Natl. Acad. Sci. U. S. A. 115, 1274–1279. 10.1073/pnas.171891011529358381 Dyckman A. J. Langevine C. M. Quesnelle C. Kempson J. Guo J. Gill P. . (2011). Imidazo[4,5-d]thiazolo[5,4-b]pyridine based inhibitors of IKK2: synthesis, SAR, PK/PD and activity in a preclinical model of rheumatoid arthritis. Bioorg. Med. Chem. Lett. 21, 383–386. 10.1016/j.bmcl.2010.10.133 Elcock A. H. (2001). Prediction of functionally important residues based solely on the computed energetics of protein structure. J. Mol. Biol. 312, 885–896. 10.1006/jmbi.2001.500911575940 Fan L.-Q. Li M.-W. Qiu Y.-J. Chen Q.-M. Jiang S.-J. Shang Y.-J. . (2018). Increasing thermal stability of glutamate decarboxylase from Escherichia coli, by site-directed saturation mutagenesis and its application in GABA production. J. Biotechnol. 278, 1–9. 10.1016/j.jbiotec.2018.04.00929660473 Ferrandi E. E. Monti D. (2017). Amine transaminases in chiral amines synthesis: recent advances and challenges. World J. Microbiol. Biotechnol. 34:13. 10.1007/s11274-017-2395-229255954 Fuchs M. Farnberger J. E. Kroutil W. (2015). The industrial age of biocatalytic transamination. Eur. J. Org. Chem. 32, 6965–6982. 10.1002/ejoc.201500852 Ghislieri D. Turner N. J. (2014). Biocatalytic approaches to the synthesis of enantiomerically pure chiral amines. Top Catal. 57, 284–300. 10.1007/s11244-013-0184-1 Godwin A. France S. P. Man H. Mangas-Sanchez J. Montgomery S. L. Sharma M. . (2017). A reductive aminase from Aspergillus oryzae. Nat. Chem. 9, 961–969. 10.1038/nchem.278228937665 Gomm A. O'Reilly E. (2018). Transaminases for chiral amine synthesis. Curr. Opin. Chem. Biol. 43, 106–112. 10.1016/j.cbpa.2017.12.00729278779 Gribenko A. V. Patel M. M. Liu J. McCallum S. A. Wang C.-Y. Makhatadze G. I. (2009). Rational stabilization of enzymes by computational redesign of surface charge-charge interactions. Proc. Natl. Acad. Sci. U. S. A. 106, 2601–2606. 10.1073/pnas.080822010619196981 Guo F. Berglund P. (2016). Transaminase biocatalysis: optimization and application. Green Chem. 19, 333–360. 10.1039/C6GC02328B Hailes H. C. Dalby P. A. Lye G. J. Baganz F. Micheletti M. Szita N. . (2010). α, α '-dihydroxy ketones and 2-amino-1,3-diols: synthetic and process strategies using biocatalysts. Curr. Org. Chem. 14:1883. 10.2174/138527210792927555 Hhne M. Bornscheuer U. T. (2009). Biocatalytic routes to optically active amines. ChemCatChem. 1, 42–51. 10.1002/cctc.200900110 Hua Y.-G. Han L.-P. Yang Q.-Q. Wang M.-J. Zhang E. Liu H.-M. (2018). A practical and efficient stereoselective synthesis of (S)-rivastigmine and (R)-rivastigmine. ChemistrySelect. 3, 1385–1387. 10.1002/slct.201703032 Huang J. Xie D.-F. Feng Y. (2017). Engineering thermostable (R)-selective amine transaminase from Aspergillus terreus through in silico design employing B-factor and folding free energy calculations. Biochem. Biophys. Res. Comm. 483, 397–402. 10.1016/j.bbrc.2016.12.13128017723 Ibarra-Molero B. Sanchez-Ruiz J. M. (2002). Genetic algo-rithm to design stabilizing surface-charge distributions in proteins. J. Phys. Chem. B. 106, 6609–6613. 10.1021/jp020483o Jiang H.-W. Chen Q. Pan J. Zheng G.-W. Xu J.-H. (2020). Rational engineering of formate dehydrogenase substrate/cofactor affinity for better performance in NADPH regeneration. Appl. Biochem. Biotechnol. 192, 530–543. 10.1007/s12010-020-03317-732405732 Jones B. J. Lim H. Y. Huang J. Kazlauskas R. J. (2017). Comparison of five protein engineering strategies for stabilizing an alpha/beta-Hydrolase. Biochemistry 56, 6521–6532. 10.1021/acs.biochem.7b0057129087185 Kroutil W. Fischereder E. M. Fuchs C. S. Lechner H. Mutti F. G. Pressnitz D. . (2013). Asymmetric preparation of prim-, sec-, and tert-amines employing selected biocatalysts. Org. Process Res. Dev. 17, 751–759. 10.1021/op400023723794796 Liu B. Qu G. Li J. Fan W. Ma J. Xu Y. . (2019). Conformational dynamics-guided loop engineering of an alcohol dehydrogenase: capture, turnover and enantioselective transformation of difficult to reduce ketones. Adv. Synth. Catal. 361, 3182–3190. 10.1002/adsc.201900249 Łyskowski A. Gruber C. Steinkellner G. Schu rmann M. Schwab H. Gruber K. . (2014). Crystal structure of an (R)-selective ω-transaminase from Aspergillus terreus. PLoS ONE 9:87350. 10.1371/journal.pone.008735024498081 Makhatadze G. I. Loladze V. V. Ermolenko D. N. Chen X.-F. Thomas S. T. (2003). Contribution of surface salt bridges to protein stability: guidelines for protein engineering. J. Mol. Biol. 327, 1135–1148. 10.1016/S0022-2836(03)00233-X12662936 Makhatadze G. I. Loladze V. V. Gribenko A. V. Lopez M. M. (2004). Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions. J. Mol. Biol. 336, 929–942. 10.1016/j.jmb.2003.12.05815095870 Malik M. S. Park E. S. Shin J. S. (2012). Features and technical applications of ω-transaminases. Appl. Microbiol. Biotechnol. 94, 1163–1171. 10.1007/s00253-012-4103-322555915 Mathew S. Yun H. (2012a). Deracemization of unnatural amino acid: homoalanine using D-amino acid oxidase and ω-transaminase. Org. Biomol. Chem. 10, 2482–2485. 10.1039/c2ob07161d22344532 Mathew S. Yun H. (2012b). ω-Transaminases for the production of optically pure amines and unnatural amino acids. ACS Catal. 2, 993–1001. 10.1021/cs300116n Matthew J. Hanania G. Gurd F. (1979). Electrostatic effects in hemoglobin: hydrogen ion equilibria in human deoxy-and oxyhemoglobin A. Biochemistry 18:1919. 10.1021/bi00577a011435457 Matthew J. B. Gurd F. (1986a). Calculation of electrostatic interactions in proteins. Methods Enzymol. 130, 413–436. 10.1016/0076-6879(86)30019-3 Matthew J. B. Gurd F. (1986b). Stabilization and destabili-zation of protein structure by charge interactions. Methods Enzymol. 130, 437–453. 10.1016/0076-6879(86)30020-X3773743 Matthew J. B. Gurd F. R. Garcia-Moreno B. E. Flanagan M. A. March K. L. Shire S. J. (1985). pH-dependent processes in protein. Crit. Rev. Biochem. Mol. Biol. 18, 91–197. 10.3109/10409238509085133 Moon S. Kim J. Koo J. Bae E. (2019). Structural and mutational analyses of psychrophilic and mesophilic adenylate kinases highlight the role of hydrophobic interactions in protein thermal stability. Struct. Dyn. 6:024702. 10.1063/1.508970731111079 Niesen F. H. Berglund H. Vedadi M. (2007). The use of differential scanning fluorimetry to detect ligand interactions that promote protein stability. Nat. Protoc. 2, 2212–2221. 10.1038/nprot.2007.32117853878 Nugent T. C. (2010). Chiral Amine Synthesis: Methods, Developments and Applications. Weinheim: Wiley-VCH. 10.1002/978352762954121321974 Pace C. N. Fu H. Fryar K. L. Landua J. Trevino S. R. Shirley B. A. . (2011). Contribution of hydrophobic interactions to protein stability. J. Mol. Biol. 408, 514–528. 10.1016/j.jmb.2011.02.05321377472 Park E. S. Dong J. Y. Shin J. S. (2014). Active site model of (R)-selective ω-transaminase and its application to the production of D-amino acids. Appl. Microbiol. Biotechnol. 98, 651–660. 10.1007/s00253-013-4846-523576035 Paul C. E. Rodríguez-Mata M. Busto E. Lavandera I. Gotor-Fernández V. Gotor V. . (2014). Transaminases applied to the synthesis of high added-value enantiopure amines. Org. Process Res. Dev. 18, 788–792. 10.1021/op4003104 Purmonen M. Valjakka J. T akkinen K. Laitinen T. Rouvinen J. (2007). Molecular dynamics studies on the thermostability of family II xylanases. Protein Eng. Des. Sel. 20, 551–559. 10.1093/protein/gzm056 Richmond T. J. (1984). Solvent accessible surface area and excluded volume in proteins: analytical equations for overlapping spheres and implications for the hydro- phobic effect. J. Mol. Biol. 178, 63–89. 10.1016/0022-2836(84)90231-6 Sanchez-Ruiz J. M. Ibarra-Molero B. Loladze V. V. Makhatadze G. I. (1999). Thermal versus guanidine induced unfolding of ubiquitin. An analysis in terms of the contributions from charge-charge interactions to protein stability. Biochemistry 38, 8138–8149. 10.1021/bi990581910387059 Sanchez-Ruiz J. M. Makhatadze G. I. (2001). To charge or not to charge? Trends Biotechnol. 19, 132–135. 10.1016/S0167-7799(00)01548-1 Schätzle S. Höhne M. Redestad E. Robins K. Bornscheuer U. T. (2009). Rapid and sensitive kinetic assay for characterization of ω-transaminases. Anal. Chem. 81, 8244–8248. 10.1021/ac901640q19739593 Schweiker K. L. Makhatadze G. I. (2009). Protein stabilization by the rational design of surface charge-charge interactions. Methods Mol. Biol. 490, 261–283. 10.1007/978-1-59745-367-7_1119157087 Strickler S. S. Gribenko A. V. Keiffer T. R. Tomlinson J. Reihle T. Loladze V. V. . (2006). Protein stability and surface electrostatics: a charged relationship. Biochemistry 45, 2761–2766. 10.1021/bi060014316503630 Svedendahl M. Branneby C. Lindbergm L. Berglund P. (2010). Reversed enantiopreference of an ω-transaminase by a single-point mutation. Chemcatchem 2, 976–980. 10.1002/cctc.201000107 Tanford C. Kirkwood J. G. (1957). Theory of protein titration curves. I. General equations for impenetrable spheres. J. Am. Chem. Soc. 79, 5333–5339. 10.1021/ja01577a001 Tufvesson P. Lima-Ramos J. Jensen J. S. Al-Haque N. Neto W. Woodley J. M. (2011). Process considerations for the asymmetric synthesis of chiral amines using transaminases. Biotechnol. Bioeng. 108, 1479–1493. 10.1002/bit.2315421455931 Turner N. J. Truppo M. D. (2010). Biocatalytic Routes to Nonracemic Chiral Amines. Weinheim: Wiley-VCH 431–459. 10.1002/9783527629541.ch14 Wang J. Tan Z. Wu M. Li J. Wu J. (2014). Improving the thermostability of a mesophilic family 10 xylanase, AuXyn10A, from Aspergillus usamii by in silico design. J. Ind. Microbiol. Biotechnol. 41, 217–225. 10.1007/s10295-014-1463-y24879480 Wu Y.-J. He H. Sun L.-Q. L'Heureux A. Chen J. Dextraze P. . (2004). Synthesis and structure-activity relationship of acrylamides as KCNQ2 potassium channel openers. J. Med. Chem. 47, 2887–2896. 10.1021/jm030582615139767 Xie D.-F. Fang H. Mei J.-Q. Gong J.-Y. Wang H.-P. Shen X.-Y. . (2018). Improving thermostability of (R)-selective amine transaminase from Aspergillus terreus through introduction of disulfide bonds. Biotechnol. Appl. Biochem. 65, 255–262. 10.1002/bab.157228639260 Xie D.-F. Yang J.-X. Lv C.-J. Mei J.-Q. Wang H.-P. Hu S. . (2019). Construction of stabilized (R)-selective amine transaminase from Aspergillus terreus by consensusmutagenesis. J. Biotechnol. 293, 8–16. 10.1016/j.jbiotec.2019.01.00730703468 Zhang L.-J. Tang X.-M. Cui D.-B. Yao Z.-Q. Gao B. (2014). A method to rationally increase protein stability based on the charge–charge interaction, with application to lipase LipK107. Protein Sci. 23, 110–116. 10.1002/pro.238824353171 Zhang Y.-X. Yang J.-W. Yu X.-Q. Hu X.-Q. Zhang H.-B. (2020). Engineering leuconostoc mesenteroides dextransucrase by inserting disulfide bridges for enhanced thermotolerance. Enzyme Microb. Technol. 139:09603. 10.1016/j.enzmictec.2020.10960332732025 Zhu D. Hua L. (2009). Biocatalytic asymmetric amination of carbonyl functional groups - a synthetic biology approach to organic chemistry. Biotechnol. J. 4, 1420–1431. 10.1002/biot.20090011019757497 Zhu W.-L. Hu S. Lv C.-J. Zhao W.-R. Wang H.-P. Mei J.-Q. . (2019). A single mutation increases the thermostability and activity of Aspergillus terreus amine transaminase. Molecules 2019, 24:1194. 10.3390/molecules2407119430934681

Funding. This research was funded by the National Natural Science Foundation of China (Nos. 32071268, 31971372, and 31670804), Zhejiang Natural Science Foundation (Nos. LQ18B060002, LZ13B060002, and LY16B060008), and Ningbo Science and Technology Innovation 2025 Major Project (No. 2020YZCXYD180178).