Balancing the Stability-Activity Trade-Off by Fine-Tuning Dehalogenase Access Tunnels

Reconstructing dynamics correlation network to simultaneously improve activity and stability of 2,3-butanediol dehydrogenase by design of distal interchain disulfide bonds

The complete enzyme catalytic cycle includes substrate binding, chemical reaction and product release, in which different dynamic conformations are adopted. Due to the complex relationship among enzyme activity, stability and dynamics, the directed evolution of enzymes for improved activity or stability commonly leads to a trade-off in stability or activity. It hence remains a challenge to engineer an enzyme to have both enhanced activity and stability. Here, we have attempted to reconstruct the dynamics correlation network involved with active center to improve both activity and stability of a 2,3-butanediol dehydrogenase (2,3-BDH) by introducing inter-chain disulfide bonds. A computational strategy was first applied to evaluate the effect of introducing inter-chain disulfide bond on activity and stability of three 2,3-BDHs, and the N258C mutation of 2,3-BDH from Corynebacterium glutamicum (CgBDH) was proved to be effective in improving both activity and stability. In the results, CgBDH-N258C showed a different unfolding curve from the wild type, with two melting temperatures (Tm) of 68.3 °C and 50.8 °C, 19.7 °C and 2 °C higher than 48.6 °C of the wild type. Its half-life was also improved by 14.8-fold compared to the wild type. Catalytic efficiency (kcat/Km) of the mutant was increased by 7.9-fold toward native substrate diacetyl and 8.8-fold toward non-native substrate 2,5-hexanedione compared to the wild type. Molecular dynamics simulations revealed that an interaction network formed by Cys258, Arg162, Ala144 and the catalytic residues was reconstructed in the mutant and the dynamics change caused by the disulfide bond could be propagated through the interactions network. This improved the enzyme stability and activity by decreasing the flexibility and locking more "reactive" pose, respectively. Further construction of mutations including A144G showing a 44-fold improvement in catalytic efficiency toward meso-2,3-BD confirmed the role of modifying dynamics correlation network in tunning enzyme activity and selectivity. This study provided important insights into the relationship among dynamics, enzyme catalysis and stability, and will be useful in the designing new enzymes with co-evolution of stability, activity and selectivity.

Modulation of the catalytic performance of OYE3 by engineering key residues at the entrance of the catalytic pocket

The amino acid residues at the entrance of the catalytic pocket may impose steric hindrance on the substrate to enter the active center of the enzyme. Based on the analysis of the three-dimensional structure of the Saccharomyces cerevisiae old yellow enzyme 3 (OYE3), four bulky residues were chosen and mutated to small amino acids. The results showed that mutation of the W116 residue had interesting impacts on the catalytic performance. All four variants became inactive for the reduction of (R)-carvone and (S)-carvone, but inverted the stereoselectivity for the reduction of (E/Z)-citral. The mutation of the F250 residue had a more positive effect on the activity and stereoselectivity. Two variants, F250A and F250S, showed excellent diastereoselectivity and activity for the reduction of (R)-carvone (de > 99%, c > 99%) and increased diastereoselectivity and activity for the reduction of (S)-carvone (de > 96%, c > 80%). One variant of the P295 residue, P295G, displayed excellent diastereoselectivity and activity only for the reduction of (R)-carvone (de > 99%, c > 99%). Mutation of the Y375 residue had a negative impact on the activity of the enzyme. These findings provide some solutions for rational enzyme engineering of OYE3.

Counteraction of stability-activity trade-off of Nattokinase through flexible region shifting

The widespread trade-off between stability and activity severely limits enzyme evolution. Although some progresses have been made to overcome this limitation, the counteraction mechanism for enzyme stability-activity trade-off remains obscure. Here, we clarified the counteraction mechanism of the Nattokinase stability-activity trade-off. A combinatorial mutant M4 was obtained by multi-strategy engineering, exhibiting a 20.7-fold improved half-life; meanwhile, the catalytic efficiency was doubled. Molecular dynamics simulation revealed that an obvious flexible region shifting in the structure of mutant M4 was occurred. The flexible region shifting which contributed to maintain the global structural flexibility, was considered to be the key factor for counteracting the stability-activity trade-off. Further analysis illustrated that the flexible region shifting was driven by region dynamical networks reshaping. This work provided deep insight into the counteraction mechanism of enzyme stability-activity trade-off, suggesting that flexible region shifting would be an effective strategy for enzyme evolution through computational protein engineering.

Coevolving stability and activity of LsCR by a single point mutation and constructing neat substrate bioreaction system

Carbonyl reductase (CR)-catalyzed bioreduction in the organic phase and the neat substrate reaction system is a lasting challenge, placing higher requirements on the performance of enzymes. Protein engineering is an effective method to enhance the properties of enzymes for industrial applications. In the present work, a single point mutation E145A on our previously constructed CR mutant LsCR_M3 , coevolved thermostability, and activity. Compared with LsCR_M3 , the catalytic efficiency k_cat /K_M of LsCR_M3 -E145A (LsCR_M4 ) was increased from 6.6 to 21.9 s^-1  mM^-1 . Moreover, E145A prolonged the half-life t_1/2 at 40°C from 4.1 to 117 h, T m ${T}_{m}$ was increased by 5°C, T 50 30 ${T}_{50}^{30}$ was increased by 14.6°C, and T_opt was increased by 15°C. Only 1 g/L of lyophilized Escherichia coli cells expressing LsCR_M4 completely reduced up to 600 g/L 2-chloro-1-(3,4-difluorophenyl)ethanone (CFPO) within 13 h at 45°C, yielding the corresponding (1S)-2-chloro-1-(3,4-difluorophenyl)ethanol ((S)-CFPL) in 99.5% ee_P , with a space-time yield of 1.0 kg/L d, the substrate to catalyst ratios (S/C) of 600 g/g. Compared with LsCR_M3 , the substrate loading was increased by 50%, with the S/C increased by 14 times. Compared with LsCR_WT , the substrate loading was increased by 6.5 times. In contrast, LsCR_M4 completely converted 600 g/L CFPO within 12 h in the neat substrate bioreaction system.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Engineering diacetylchitobiose deacetylase from Pyrococcus horikoshii towards an efficient glucosamine production

In this study, semi-rational design based on site-directed saturation mutagenesis and surface charge modification was used to improve the catalytic efficiency of the diacetylchitobiose deacetylase derived from Pyrococcus horikoshii (PhDac). PhDac mutant M14, which was screened by site-directed saturation mutagenesis, showed a ~ 2.21 -fold enhanced catalytic efficiency (kcat/Km) and the specific activity was improved by 70.02%. To keep the stability of glucosamine (GlcN), we reduced the optimal pH of M14 by modifying the surface charge from -35 to -59 to obtain mutant M20, whose specific activity reached 2 -fold of the wild-type. The conversion rate of N-acetylglucosamine (GlcNAc) to GlcN catalyzed by M20 reached 94.3%. Moreover, the decline of GlcN production was slowed down by the reduction of pH when temperature was higher than 50 ℃. Our results would accelerate the process of industrial production of GlcN by biocatalysis.

Stabilizing AqdC, a Pseudomonas Quinolone Signal-Cleaving Dioxygenase from Mycobacteria, by FRESCO-Based Protein Engineering

The mycobacterial PQS dioxygenase AqdC, a cofactor-less protein with an α/β-hydrolase fold, inactivates the virulence-associated quorum-sensing signal molecule 2-heptyl-3-hydroxy-4(1H)-quinolone (PQS) produced by the opportunistic pathogen Pseudomonas aeruginosa and is therefore a potential anti-virulence tool. We have used computational library design to predict stabilizing amino acid replacements in AqdC. While 57 out of 91 tested single substitutions throughout the protein led to stabilization, as judged by increases in T app m of >2 °C, they all impaired catalytic activity. Combining substitutions, the proteins AqdC-G40K-A134L-G220D-Y238W and AqdC-G40K-G220D-Y238W showed extended half-lives and the best trade-off between stability and activity, with increases in T app m of 11.8 and 6.1 °C and relative activities of 22 and 72 %, respectively, compared to AqdC. Molecular dynamics simulations and principal component analysis suggested that stabilized proteins are less flexible than AqdC, and the loss of catalytic activity likely correlates with an inability to effectively open the entrance to the active site.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

Exploiting the activity-stability trade-off of glucose oxidase from Aspergillus niger using a simple approach to calculate thermostability of mutants

Glucose oxidase (Gox) is a biocatalyst that is widely applied in the food industry, as well as other biotechnological industries. However, the industrial application of Gox is hampered by its low thermostability and activity. Here, we aimed to improve the thermostability of GoxM4 from Aspergillus niger without reducing its activity due to the activity-stability trade-off. A simple and effective approach combining enzyme activity and structure stability was adopted to evaluate the thermostability of GoxM4 and its mutants. After four rounds of computer-aided rational design, the best mutant, GoxM8, was obtained. The melting temperature (T_m) of GoxM8 was increased by 9 °C compared with GoxM4. The catalytic efficiency of GoxM8 was similar to GoxM4, suggesting that the enzyme activity-stability trade-off was counteracted. To explore its mechanism, we performed molecular dynamics simulations of GoxM4 and its mutants. Our findings provided a typical example for researching the enzyme activity-stability trade-off.

Systematic Evaluation of Imine-Reducing Enzymes: Common Principles in Imine Reductases, β-Hydroxy Acid Dehydrogenases, and Short-Chain Dehydrogenases/ Reductases

The enzymatic, asymmetric reduction of imines is catalyzed by imine reductases (IREDs), members of the short-chain dehydrogenase/reductase (SDR) family, and β-hydroxy acid dehydrogenase (βHAD) variants. Systematic evaluation of the structures and substrate-binding sites of the three enzyme families has revealed four common principles for imine reduction: structurally conserved cofactor-binding domains; tyrosine, aspartate, or glutamate as proton donor; at least four characteristic flanking residues that adapt the donor's pKa and polarize the substrate; and a negative electrostatic potential in the substrate-binding site to stabilize the transition state. As additional catalytically relevant positions, we propose alternative proton donors in IREDs and βHADs as well as proton relays in IREDs, βHADs, and SDRs. The functional role of flanking residues was experimentally confirmed by alanine scanning of the imine-reducing SDR from Zephyranthes treatiae. Mutating the "gatekeeping" phenylalanine at standard position 200 resulted in a tenfold increase in imine-reducing activity.

Decoding the intricate network of molecular interactions of a hyperstable engineered biocatalyst

Computational design of protein catalysts with enhanced stabilities for use in research and enzyme technologies is a challenging task. Using force-field calculations and phylogenetic analysis, we previously designed the haloalkane dehalogenase DhaA115 which contains 11 mutations that confer upon it outstanding thermostability (T m = 73.5 °C; ΔT m > 23 °C). An understanding of the structural basis of this hyperstabilization is required in order to develop computer algorithms and predictive tools. Here, we report X-ray structures of DhaA115 at 1.55 Å and 1.6 Å resolutions and their molecular dynamics trajectories, which unravel the intricate network of interactions that reinforce the αβα-sandwich architecture. Unexpectedly, mutations toward bulky aromatic amino acids at the protein surface triggered long-distance (∼27 Å) backbone changes due to cooperative effects. These cooperative interactions produced an unprecedented double-lock system that: (i) induced backbone changes, (ii) closed the molecular gates to the active site, (iii) reduced the volumes of the main and slot access tunnels, and (iv) occluded the active site. Despite these spatial restrictions, experimental tracing of the access tunnels using krypton derivative crystals demonstrates that transport of ligands is still effective. Our findings highlight key thermostabilization effects and provide a structural basis for designing new thermostable protein catalysts.

Caver Web 1.0: identification of tunnels and channels in proteins and analysis of ligand transport

Caver Web 1.0 is a web server for comprehensive analysis of protein tunnels and channels, and study of the ligands’ transport through these transport pathways. Caver Web is the first interactive tool allowing both the analyses within a single graphical user interface. The server is built on top of the abundantly used tunnel detection tool Caver 3.02 and CaverDock 1.0 enabling the study of the ligand transport. The program is easy-to-use as the only required inputs are a protein structure for a tunnel identification and a list of ligands for the transport analysis. The automated guidance procedures assist the users to set up the calculation in a way to obtain biologically relevant results. The identified tunnels, their properties, energy profiles and trajectories for ligands’ passages can be calculated and visualized. The tool is very fast (2–20 min per job) and is applicable even for virtual screening purposes. Its simple setup and comprehensive graphical user interface make the tool accessible for a broad scientific community. The server is freely available at https://loschmidt.chemi.muni.cz/caverweb.

Directed evolution methods for overcoming trade-offs between protein activity and stability

Engineered proteins are being widely developed and employed in applications ranging from enzyme catalysts to therapeutic antibodies. Directed evolution, an iterative experimental process composed of mutagenesis and library screening, is a powerful technique for enhancing existing protein activities and generating entirely new ones not observed in nature. However, the process of accumulating mutations for enhanced protein activity requires chemical and structural changes that are often destabilizing, and low protein stability is a significant barrier to achieving large enhancements in activity during multiple rounds of directed evolution. Here we highlight advances in understanding the origins of protein activity/stability trade-offs for two important classes of proteins (enzymes and antibodies) as well as innovative experimental and computational methods for overcoming such trade-offs. These advances hold great potential for improving the generation of highly active and stable proteins that are needed to address key challenges related to human health, energy and the environment.

Tunnel engineering to accelerate product release for better biomass-degrading abilities in lignocellulolytic enzymes

BACKGROUND

For enzymes with buried active sites, transporting substrates/products ligands between active sites and bulk solvent via access tunnels is a key step in the catalytic cycle of these enzymes. Thus, tunnel engineering is becoming a powerful strategy to refine the catalytic properties of these enzymes. The tunnel-like structures have been described in enzymes catalyzing bulky substrates like glycosyl hydrolases, while it is still uncertain whether these structures involved in ligands exchange. Till so far, no studies have been reported on the application of tunnel engineering strategy for optimizing properties of enzymes catalyzing biopolymers.

RESULTS

In this study, xylanase S7-xyl (PDB: 2UWF) with a deep active cleft was chosen as a study model to evaluate the functionalities of tunnel-like structures on the properties of biopolymer-degrading enzymes. Three tunnel-like structures in S7-xyl were identified and simultaneously reshaped through multi-sites saturated mutagenesis; the most advantageous mutant 254RL1 (V207N/Q238S/W241R) exhibited 340% increase in specific activity compared to S7-xyl. Deconvolution analysis revealed that all three mutations contributed synergistically to the improved activity of 254RL1. Enzymatic characterization showed that larger end products were released in 254RL1, while substrate binding and structural stability were not changed. Dissection of the structural alterations revealed that both the tun_1 and tun_2 in 254RL1 have larger bottleneck radius and shorter length than those of S7-xyl, suggesting that these tunnel-like structures may function as products transportation pathways. Attributed to the improved catalytic efficiency, 254RL1 represents a superior accessory enzyme to enhance the hydrolysis efficiency of cellulase towards different pretreated lignocellulose materials. In addition, tunnel engineering strategy was also successfully applied to improve the catalytic activities of three other xylanases including xylanase NG27-xyl from Bacillus sp. strain NG-27, TSAA1-xyl from Geobacillus sp. TSAA1 and N165-xyl from Bacillus sp. N16-5, with 80%, 20% and 170% increase in specific activity, respectively.

CONCLUSIONS

This study represents a pilot study of engineering and functional verification of tunnel-like structures in enzymes catalyzing biopolymer. The specific activities of four xylanases with buried active sites were successfully improved by tunnel engineering. It is highly likely that tunnel reshaping can be used to engineer better biomass-degrading abilities in other lignocellulolytic enzymes with buried active sites.

Deciphering the Structural Basis of High Thermostability of Dehalogenase from Psychrophilic Bacterium Marinobacter sp. ELB17

Haloalkane dehalogenases are enzymes with a broad application potential in biocatalysis, bioremediation, biosensing and cell imaging. The new haloalkane dehalogenase DmxA originating from the psychrophilic bacterium Marinobacter sp. ELB17 surprisingly possesses the highest thermal stability (apparent melting temperature Tm,app = 65.9 °C) of all biochemically characterized wild type haloalkane dehalogenases belonging to subfamily II. The enzyme was successfully expressed and its crystal structure was solved at 1.45 Å resolution. DmxA structure contains several features distinct from known members of haloalkane dehalogenase family: (i) a unique composition of catalytic residues; (ii) a dimeric state mediated by a disulfide bridge; and (iii) narrow tunnels connecting the enzyme active site with the surrounding solvent. The importance of narrow tunnels in such paradoxically high stability of DmxA enzyme was confirmed by computational protein design and mutagenesis experiments.

Engineering enzyme access tunnels

Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.

Exploiting correlated molecular-dynamics networks to counteract enzyme activity-stability trade-off

Rigidifying flexible sites is a powerful method to improve enzyme stability. However, if the highly flexible regions form the active site, modifying them risks losing activity due to the activity–stability trade-off. We hypothesized here that regions outside the active site whose dynamics were highly correlated to flexible active sites, would provide good targets for stabilizing mutations. To test this hypothesis, six variants were constructed in the 3M variant of Escherichia coli transketolase. The best variant had a 10.8-fold improved half-life at 55 °C, and increased the T_m and T_agg by 3 °C and 4.3 °C, respectively. The variants even increased the activity, by up to threefold. This study highlights how protein engineering strategies could be potentially improved by considering long-range dynamics.

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.

Filling the Void: Introducing Aromatic Interactions into Solvent Tunnels To Enhance Lipase Stability in Methanol

An enhanced stability of enzymes in organic solvents is desirable under industrial conditions. The potential of lipases as biocatalysts is mainly limited by their denaturation in polar alcohols. In this study, we focused on selected solvent tunnels in lipase from Geobacillus stearothermophilus T6 to improve its stability in methanol during biodiesel synthesis. Using rational mutagenesis, bulky aromatic residues were incorporated to occupy solvent channels and induce aromatic interactions leading to a better inner core packing. The chemical and structural characteristics of each solvent tunnel were systematically analyzed. Selected residues were replaced with Phe, Tyr, or Trp. Overall, 16 mutants were generated and screened in 60% methanol, from which 3 variants showed an enhanced stability up to 81-fold compared with that of the wild type. All stabilizing mutations were found in the longest tunnel detected in the "closed-lid" X-ray structure. The combination of Phe substitutions in an A187F/L360F double mutant resulted in an increase in unfolding temperature (T_m ) of 7°C in methanol and a 3-fold increase in biodiesel synthesis yield from waste chicken oil. A kinetic analysis with p-nitrophenyl laurate revealed that all mutants displayed lower hydrolysis rates (k _cat), though their stability properties mostly determined the transesterification capability. Seven crystal structures of different variants were solved, disclosing new π-π or CH/π intramolecular interactions and emphasizing the significance of aromatic interactions for improved solvent stability. This rational approach could be implemented for the stabilization of other enzymes in organic solvents.IMPORTANCE Enzymatic synthesis in organic solvents holds increasing industrial opportunities in many fields; however, one major obstacle is the limited stability of biocatalysts in such a denaturing environment. Aromatic interactions play a major role in protein folding and stability, and we were inspired by this to redesign enzyme voids. The rational protein engineering of solvent tunnels of lipase from Geobacillus stearothermophilus is presented here, offering a promising approach to introduce new aromatic interactions within the enzyme core. We discovered that longer tunnels leading from the surface to the enzyme active site were more beneficial targets for mutagenesis for improving lipase stability in methanol during biodiesel biosynthesis. A structural analysis of the variants confirmed the generation of new interactions involving aromatic residues. This work provides insights into stability-driven enzyme design by targeting the solvent channel void.

Enhanced trypsin thermostability in Pichia pastoris through truncating the flexible region

Background

High thermostability is required for trypsin to have wider industrial applications. Target mutagenesis at flexible regions has been proved to be an efficient protein engineering method to enhance the protein thermostability.

Results

The flexible regions in porcine trypsin were predicted using the methods including molecular dynamic simulation, FlexPred, and FoldUnfold. The amino acids 78–90 was predicted to be the highly flexible region simultaneously by the three methods and hence selected to be the mutation target. We constructed five variants (D3, D5, D7, D9, and D11) by truncating the region. And the variant D9 showed higher thermostability, with a 5 °C increase in T_opt, 5.8 °C rise in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{50}^{10}$$\end{document}T5010, and a 4.5 °C rise in T_m, compared to the wild-type. Moreover, the half-life value of the variant D9 was also found to be dramatically improved by 46 min. Circular dichroism and intrinsic fluorescence indicated that the structures had no significant change between the variant D9 and the wild-type. The surface hydrophobicity of D9 was measured to be lower than that of wild-type, indicating the increased hydrophobic interaction, which could have contributed to the improved thermostability of D9.

Conclusions

These results showed that the thermostability of variant D9 was increased. The variant D9 could be expected to be a promising tool enzyme for its wider industrial applications. The method of truncating the flexible region used in our study has the potential to be used for enhancing the thermostability of other proteins.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1012-x) contains supplementary material, which is available to authorized users.

Rational and Semirational Protein Design

This mini review gives an overview over different design approaches and methodologies applied in rational and semirational enzyme engineering. The underlying principles for engineering novel activities, enantioselectivity, substrate specificity, stability, and pH optimum are summarized.

Computational Analysis of Protein Tunnels and Channels

Protein tunnels connecting the functional buried cavities with bulk solvent and protein channels, enabling the transport through biological membranes, represent the structural features that govern the exchange rates of ligands, ions, and water solvent. Tunnels and channels are present in a vast number of known proteins and provide control over their function. Modification of these structural features by protein engineering frequently provides proteins with improved properties. Here we present a detailed computational protocol employing the CAVER software that is applicable for: (1) the analysis of tunnels and channels in protein structures, and (2) the selection of hot-spot residues in tunnels or channels that can be mutagenized for improved activity, specificity, enantioselectivity, or stability.

Ancestral Haloalkane Dehalogenases Show Robustness and Unique Substrate Specificity

Ancestral sequence reconstruction (ASR) represents a powerful approach for empirical testing structure-function relationships of diverse proteins. We employed ASR to predict sequences of five ancestral haloalkane dehalogenases (HLDs) from the HLD-II subfamily. Genes encoding the inferred ancestral sequences were synthesized and expressed in Escherichia coli, and the resurrected ancestral enzymes (AncHLD1-5) were experimentally characterized. Strikingly, the ancestral HLDs exhibited significantly enhanced thermodynamic stability compared to extant enzymes (ΔT_m up to 24 °C), as well as higher specific activities with preference for short multi-substituted halogenated substrates. Moreover, multivariate statistical analysis revealed a shift in the substrate specificity profiles of AncHLD1 and AncHLD2. This is extremely difficult to achieve by rational protein engineering. The study highlights that ASR is an efficient approach for the development of novel biocatalysts and robust templates for directed evolution.

Simultaneous engineering of an enzyme's entrance tunnel and active site: the case of monoamine oxidase MAO-N

A new directed evolution approach is presented to enhance the activity of an enzyme and to manipulate stereoselectivity by focusing iterative saturation mutagenesis (ISM) simultaneously on residues lining the entrance tunnel and the binding pocket. This combined mutagenesis strategy was applied successfully to the monoamine oxidase from Aspergillus niger (MAO-N) in the reaction of sterically demanding substrates which are of interest in the synthesis of chiral pharmaceuticals based on the benzo-piperidine scaffold. Reversal of enantioselectivity of Turner-type deracemization was achieved in the synthesis of (S)-1,2,3,4-tetrahydro-1-methyl-isoquinoline, (S)-1,2,3,4-tetrahydro-1-ethylisoquinoline and (S)-1,2,3,4-tetrahydro-1-isopropylisoquinoline. Extensive molecular dynamics simulations indicate that the altered catalytic profile is due to increased hydrophobicity of the entrance tunnel acting in concert with the altered shape of the binding pocket.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

AnimoAminoMiner: Exploration of Protein Tunnels and their Properties in Molecular Dynamics

In this paper we propose a novel method for the interactive exploration of protein tunnels. The basic principle of our approach is that we entirely abstract from the 3D/4D space the simulated phenomenon is embedded in. A complex 3D structure and its curvature information is represented only by a straightened tunnel centerline and its width profile. This representation focuses on a key aspect of the studied geometry and frees up graphical estate to key chemical and physical properties represented by surrounding amino acids. The method shows the detailed tunnel profile and its temporal aggregation. The profile is interactively linked with a visual overview of all amino acids which are lining the tunnel over time. In this overview, each amino acid is represented by a set of colored lines depicting the spatial and temporal impact of the amino acid on the corresponding tunnel. This representation clearly shows the importance of amino acids with respect to selected criteria. It helps the biochemists to select the candidate amino acids for mutation which changes the protein function in a desired way. The AnimoAminoMiner was designed in close cooperation with domain experts. Its usefulness is documented by their feedback and a case study, which are included.