Comparison of Five Protein Engineering Strategies for Stabilizing an α/β-Hydrolase

Enhancing the Thermal Stability and Enzyme Activity of Ketopantoate Hydroxymethyltransferase through Interface Modification Engineering

Ketopantoate hydroxymethyltransferase (KPHMT) plays a pivotal role in d-pantothenic acid biosynthesis. Most KPHMTs are homodecamers with low thermal stability, posing challenges for protein engineering and limiting output enhancement. Previously, a high-enzyme activity KPHMT mutant (K25A/E189S) from Corynebacterium glutamicum was screened as mother strain (M0). Building upon this strain, our study focused on interface engineering modifications, employing a multifaceted approach including integrating folding-free energy calculation, B-factor analysis, and conserved site analysis. Preliminary screening led to the selection of five mutants in the interface─E106S, E98T, E98N, S247I, and S247D─showing improved thermal stability, culminating in the double-site mutant M8 (M0-E98N/S247D). M8 exhibited a T1/2 value of 288.79 min at 50 °C, showing a 3.29-fold increase compared to M0. Meanwhile, the Tm value of M8 was elevated from 53.2 to 59.6 °C. Investigations of structural and molecular dynamics simulations revealed alterations in surface electrostatic charge distribution and the formation of increased hydrogen bonds between subunits, contributing to enhanced thermal stability. This investigation corroborates the efficacy of interface engineering modifications in bolstering KPHMT stability while showing its potential for positively impacting industrial d-pantothenic acid synthesis.

Computational Enzyme Redesign Enhances Tolerance to Denaturants for Peptide C-Terminal Amidation

The escalating demand for biocatalysts in pharmaceutical and biochemical applications underscores the critical imperative to enhance enzyme activity and durability under high denaturant concentrations. Nevertheless, the development of a practical computational redesign protocol for improving enzyme tolerance to denaturants is challenging due to the limitations of relying solely on model-driven approaches to adequately capture denaturant-enzyme interactions. In this study, we introduce an enzyme redesign strategy termed GRAPE_DA, which integrates multiple data-driven and model-driven computational methods to mitigate the sampling biases inherent in a single approach and comprehensively predict beneficial mutations on both the protein surface and backbone. To illustrate the methodology's effectiveness, we applied it to engineer a peptidylamidoglycolate lyase, resulting in a variant exhibiting up to a 24-fold increase in peptide C-terminal amidation activity under 2.5 M guanidine hydrochloride. We anticipate that this integrated engineering strategy will facilitate the development of enzymatic peptide synthesis and functionalization under denaturing conditions and highlight the role of engineering surface residues in governing protein stability.

Enzymatic Degradation of Deoxynivalenol with the Engineered Detoxification Enzyme Fhb7

In the era of global climate change, the increasingly severe Fusarium head blight (FHB) and deoxynivalenol (DON) contamination have caused economic losses and brought food and feed safety concerns. Recently, an FHB resistance gene Fhb7 coding a glutathione-S transferase (GST) to degrade DON by opening the critical toxic epoxide moiety was identified and opened a new window for wheat breeding and DON detoxification. However, the poor stability of Fhb7 and the elusiveness of the catalytic mechanism hinder its practical application. Herein, we report the first structure of Fhb7 at 2.41 Å and reveal a unique catalytic mechanism of epoxide opening transformation in GST family proteins. Furthermore, variants V29P and M10 showed that 5.5-fold and 266.7-fold longer half-life time than wild-type, respectively, were identified. These variants offer broad substrate scope, and the engineered biosafe Bacillus subtilis overexpressing the variants shows excellent DON degradation performance, exhibiting potential at bacterium engineering to achieve DON detoxification in the feed and biomedicine industry. This work provides a profound mechanistic insight into the enzymatic activities of Fhb7 and paves the way for further utilizing Fhb7-related enzymes in crop breeding and DON detoxification by synthetic biology.

Thermal Stabilization of a Bacterial Zn(II)-Dependent Phospholipase C through Consensus Sequence Design

Proteins' extraordinary performance in recognition and catalysis has led to their use in a range of applications. However, proteins obtained from natural sources are oftentimes not suitable for direct use in industrial or diagnostic setups. Natural proteins, evolved to optimally perform a task in physiological conditions, usually lack the stability required to be used in harsher conditions. Therefore, the alteration of the stability of proteins is commonly pursued in protein engineering studies. Here, we achieved a substantial thermal stabilization of a bacterial Zn(II)-dependent phospholipase C by consensus sequence design. We retrieved and analyzed sequenced homologues from different sources, selecting a subset of examples for expression and characterization. A non-natural consensus sequence showed the highest stability and activity among those tested. Comparison of the stability parameters of this stabilized mutant and other natural variants bearing similar mutations allows us to pinpoint the sites most likely to be responsible for the enhancement. Point mutations in these sites alter the unfolding process of the consensus sequence. We show that the stabilized version of the protein retains full activity even in harsh oil degumming conditions, making it suitable for industrial applications.

Improving the thermal stability and catalytic activity of ulvan lyase by the combination of FoldX and KnowVolution campaign

Thermal stability is one of the most important properties of ulvan lyases for their application in algae biomass degradation. The Knowledge gaining directed eVolution (KnowVolution) protein engineering strategy could be employed to improve thermostability of ulvan lyase with less screening effort. Herein, the unfolding free energies (ΔΔG) of the loop region were calculated using FoldX and four sites (D103, G104, T113, Q229) were selected for saturation mutagenesis, resulting in the identification of a favorable single-site mutant Q229M. Subsequently, iteration mutation was carried out with the mutant N57P (previously obtained by our group) to further enhance the performance of ulvan lyase. The results showed that the most beneficial variant N57P/Q229M exhibited a 1.67-fold and 2-fold increase in residual activity compared to the wild type after incubation at 40 °C and 50 °C for 1 h, respectively. In addition, the variant produced 1.06 mg/mL of reducing sugar in 2 h, which was almost four times as much as the wild type. Molecular dynamics simulations revealed that N57P/Q229M mutant enhanced the structural rigidity by augmenting intramolecular hydrogen bonds. Meanwhile, the shorter proton transmission distance between the general base of the enzyme and the substrate contributed to the glycosidic bond breakage. Our research showed that in silico saturation mutagenesis using position scan module in FoldX allowed for faster screening of mutants with improved thermal stability, and combining it with KnowVolution enabled a balanced effect of thermal stability and enzyme activity in protein engineering.

Improving the Activity and Stability of Serine Protease ThAPT3 by Alleviating Self-Cleavage and Its Application in Deproteinization of Shrimp Shells

The self-cleavage properties of proteases result in low activity and instability, which limit their industrial application. In this study, the serine protease ThAPT3 from Torrubiella hemipterigena was successfully expressed in Komagataella phaffii. We investigated the self-degradation mechanism of ThAPT3 and presented a rational strategy to alleviate self-cleavage. A major self-degradation site (Leu²³⁸-Met²³⁹) and a primary autolysis region were identified. The autolysis regions (loop₁₈, α₈-helix, and loop₁₉) were redesigned and optimized using loop transplantation, energy calculations, surface cavity optimization, and loop anchoring. A triple-superposition mutant, ThAPT3-M9 (M²³⁹GKDGAVAAGLC²⁵⁰ → M²³⁹TLNRTTAANAC²⁵⁰/A251E/A254Q/R259L/A267E/S280N), was obtained. Compared to the wild type, the autolysis of M9 was significantly alleviated, and its half-life at 60 °C was increased approximately 39-fold (from 1.6 to 62.4 min). The optimal temperature and specific activity of M9 increased by 5 °C (from 60 to 65 °C) and 62% (4985 vs 3078 U/mg), respectively. M9 showed significant advantages in shrimp shell deproteinization.

Turning thermostability of Aspergillus terreus (R)-selective transaminase At-ATA by synthetic shuffling

As versatile and green biocatalysts for the asymmetric amination of ketones, the insufficient thermostability of transaminases always limits its broad application in the pharmaceutical and fine chemical industries. Here, synthetic shuffling technology was used to enhance stability of (R)-selective transaminase from Aspergillus terreus. The results showed that 30 out of 5000 mutants had improved thermostability by color-based screening method, among which mutants with residual enzyme activity higher than 50% at 45 °C for 10 min were selected for further analysis. Especially, the half-inactivation temperature (T₅₀¹⁰), half-life (t_1/2), and melting temperature (T_m) of the best mutant M14 (M280C-H210N-M150C-F115L) were 13.7 °C, 165.8 min, and 13.9 °C higher than that of the wild type (WT), respectively. M14 also exhibited a significant biocatalytic efficiency toward acetophenone and 1-acetylnaphthalene, the yield of which were 265.6% and 117.5% higher than WT, respectively. Based on molecular dynamics simulation, improved catalytic efficiency of M14 could be attributed to its increased hydrogen bonds interaction around the mutation sites. Additionally, the introduction of disulfide bond combined with above mutations has a synergistic effect on the improved protein thermostability.

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.

Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design

The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.

Protein Stability: Enhancement and Measurement

This chapter defines protein stability, emphasizes its importance, and surveys the field of protein stabilization, with summary reference to a selection of 2014-2021 publications. One can enhance stability, particularly by protein engineering strategies but also by chemical modification and by other means. General protocols are set out on how to measure a given protein's (i) kinetic thermal stability and (ii) oxidative stability and (iii) how to undertake chemical modification of a protein in solution.

Enzymatic Conjugation of Modified RNA Fragments by Ancestral RNA Ligase AncT4_2

Oligonucleotide therapeutics have great potential as a next-generation approach to treating intractable diseases. Large quantities of modified DNA/RNA containing xenobiotic nucleic acids (XNAs) must be synthesized before clinical application. In this study, the ancestral RNA ligase AncT4_2 was designed by ancestral sequence reconstruction (ASR) to perform the conjugation reaction of modified RNA fragments. AncT4_2 had superior properties to native RNA ligase 2 from T4 phage (T4Rnl2), including high productivity, a >2.5-fold-higher turnover number, and >10°C higher thermostability. One remarkable point is the broad substrate selectivity of AncT4_2; the activity of AncT4_2 toward 17 of the modified RNA fragments was higher than that of T4Rnl2. The activity was estimated by measuring the conjugation reaction of two RNA strands, 3'-OH (12 bp) and 5'-PO₄ (12 bp), in which the terminal and penultimate positions of the 3'-OH fragment and the first and second positions of the 5'-PO₄ fragment were substituted by 2'-fluoro, 2'-O-methyl, 2'-O-methoxyethyl, and 2'-H, respectively. The enzymatic properties of AncT4_2 allowed the enzyme to conjugate large quantities of double-stranded RNA coding for patisiran (>400 μM level), which was formed by four RNA fragments containing 2'-OMe-substituted nucleic acids. Structural analysis of modeled AncT4_2 suggested that protein dynamics were changed by mutation to Gly or indel during ASR and that this may positively impact the conjugation of modified RNA fragments with the enzyme. AncT4_2 is expected to be a key biocatalyst in synthesizing RNA therapeutics by an enzymatic reaction. IMPORTANCE RNA therapeutics is one of the next-generation medicines for treating various diseases. Our designed ancestral RNA ligase AncT4_2 exhibited excellent enzymatic properties, such as high thermal stability, productivity, specific activity, and broad substrate selectivity compared to native enzymes. These advantages create the potential for AncT4_2 to be applied in conjugating the modified RNA fragments containing various xenobiotic nucleic acids. In addition, patisiran, a known polyneuropathy therapeutic, could be synthesized from four fragmented oligonucleotides at a preparative scale. Taken together, these findings indicate AncT4_2 could open the door to synthesizing RNA therapeutics by enzymatic reaction at large-scale production.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Interactions between Artificial Channel Protein, Water Molecules, and Ions Based on Theoretical Approaches

Contemporary techniques of molecular modeling allow for rational design of several specific classes of artificial proteins. Transmembrane channels are among these classes. A recent successful synthesis of self-assembling, highly symmetrical 12- or 16-helix channels by David Baker’s group prompted us to study interactions between one of these proteins, TMHC6, and low-molecular-weight components of the environment: water molecules and ions. To examine protein stability in a polar environment, molecular dynamics (MD) with classical force fields of the AMBER family was employed. Further characteristics of the chosen interactions were obtained using interaction energy calculations with usage of partially polarizable GFN-FF force field of Spicher and Grimme, symmetry-adapted perturbation theory (SAPT) and atoms in molecules (AIM) approaches for models of residues from the channel entry, crucial for interactions with water molecules and ions. The comparison of the interaction energy values between the gas phase and solvent reaction field gives the quantitative estimation of the strength of the interactions. The energy decomposition via the SAPT method showed that the electrostatics forces play a dominant role in the substructure stabilization. An application of the AIM theory enabled a description of the intermolecular hydrogen bonds and other noncovalent interactions.

Engineering and evaluation of thermostable Is PETase variants for PET degradation

Polyethylene terephthalate (PET) is a mass‐produced petroleum‐based synthetic polymer. Enzymatic PET degradation using, for example, Ideonella sakaiensis PETase (IsPETase) can be a more environmentally friendly and energy‐saving alternative to the chemical recycling of PET. However, IsPETase is a mesophilic enzyme with an optimal reaction temperature lower than the glass transition temperature (T_g) of PET, where the amorphous polymers can be readily accessed for enzymatic breakdown. In this study, we used error‐prone PCR to generate a mutant library based on a thermostable triple mutant (TM) of IsPETase. The library was screened against the commercially available polyester‐polyurethane Impranil DLN W 50 for more thermostable IsPETase variants, yielding four variants with higher melting points. The most promising IsPETaseTM^K95N/F201I variant had a 5.0°C higher melting point than IsPETaseTM. Although this variant showed a slightly lower activity on PET at lower incubation temperatures, its increased thermostability makes it a more active PET hydrolase at higher reaction temperatures up to 60°C. Several other variants were compared and combined with selected previously published IsPETase mutants in terms of thermostability and hydrolytic activity against PET nanoparticles and amorphous PET films. Our findings indicate that thermostability is one of the most important characteristics of an effective PET hydrolase.

Enhanced Thermostability of Pseudomonas nitroreducens Isoeugenol Monooxygenase by the Combinatorial Strategy of Surface Residue Replacement and Consensus Mutagenesis

Vanillin has many applications in industries. Isoeugenol monooxygenase (IEM) can catalyze the oxidation of isoeugenol to vanillin in the presence of oxygen under mild conditions. However, the low thermal stability of IEM limits its practical application in the biosynthesis of natural vanillin. Herein, two rational strategies were combined to improve the thermostability of IEM from Pseudomonas nitroreducens Jin1. Two variants (K83R and K95R) with better thermostability and one mutant (G398A) with higher activity were identified from twenty candidates based on the Surface Residue Replacement method. According to the Consensus Mutagenesis method, one mutant (I352R) with better thermostability and another mutant (L273F) with higher activity were also identified from nine candidates. After combinatorial mutation, a triple mutant K83R/K95R/L273F with the best thermostability and catalytic efficiency was generated. Compared with the wild-type IEM, the thermal inactivation half-lives (t1/2) of K83R/K95R/L273F at 25 °C, 30 °C, and 35 °C increased 2.9-fold, 11.9-fold, and 24.7-fold, respectively. Simultaneously, it also exhibited a 4.8-fold increase in kcat, leading to a 1.2-fold increase in catalytic efficiency (kcat/Km). When the whole cell of K83R/K95R/L273F was applied to the biotransformation of isoeugenol on preparative scale, the vanillin concentration reached 240.1 mM with space-time yield of 109.6 g/L/d, and vanillin was achieved in 77.6% isolated yield and >99% purity.

Partial Consensus Design and Enhancement of Protein Function by Secondary-Structure-Guided Consensus Mutations

Consensus design (CD) is a representative sequence-based protein design method that enables the design of highly functional proteins by analyzing vast amounts of protein sequence data. This study proposes a partial consensus design (PCD) of a protein as a derivative approach of CD. The method replaces the target protein sequence with a consensus sequence in a secondary-structure-dependent manner (i.e., regionally dependent and divided into α-helix, β-sheet, and loop regions). In this study, we generated several artificial partial consensus l-threonine 3-dehydrogenases (PcTDHs) by PCD using the TDH from Cupriavidus necator (CnTDH) as a target protein. Structural and functional analysis of PcTDHs suggested that thermostability would be independently improved when consensus mutations are introduced into the loop region of TDHs. On the other hand, enzyme kinetic parameters (k_cat/K_m) and average productivity would be synergistically enhanced by changing the combination of the mutations-replacement of one region of CnTDH with a consensus sequence provided only negative effects, but the negative effects were nullified when the two regions were replaced simultaneously. Taken together, we propose the hypothesis that there are protein regions that encode individual protein properties, such as thermostability and activity, and that the introduction of consensus mutations into these regions could additively or synergistically modify their functions.

Assessment of the PETase conformational changes induced by poly(ethylene terephthalate) binding

Recently, a bacterium strain of Ideonella sakaiensis was identified with the uncommon ability to degrade the poly(ethylene terephthalate) (PET). The PETase from I. sakaiensis strain 201-F6 (IsPETase) catalyzes the hydrolysis of PET converting it to mono(2-hydroxyethyl) terephthalic acid (MHET), bis(2-hydroxyethyl)-TPA (BHET), and terephthalic acid (TPA). Despite the potential of this enzyme for mitigation or elimination of environmental contaminants, one of the limitations of the use of IsPETase for PET degradation is the fact that it acts only at moderate temperature due to its low thermal stability. Besides, molecular details of the main interactions of PET in the active site of IsPETase remain unclear. Herein, molecular docking and molecular dynamics (MD) simulations were applied to analyze structural changes of IsPETase induced by PET binding. Results from the essential dynamics revealed that the β1-β2 connecting loop is very flexible. This loop is located far from the active site of IsPETase and we suggest that it can be considered for mutagenesis to increase the thermal stability of IsPETase. The free energy landscape (FEL) demonstrates that the main change in the transition between the unbound to the bound state is associated with the β7-α5 connecting loop, where the catalytic residue Asp206 is located. Overall, the present study provides insights into the molecular binding mechanism of PET into the IsPETase structure and a computational strategy for mapping flexible regions of this enzyme, which can be useful for the engineering of more efficient enzymes for recycling plastic polymers using biological systems.

Revolutionizing enzyme engineering through artificial intelligence and machine learning

The combinatorial space of an enzyme sequence has astronomical possibilities and exploring it with contemporary experimental techniques is arduous and often ineffective. Multi-target objectives such as concomitantly achieving improved selectivity, solubility and activity of an enzyme have narrow plausibility under approaches of restricted mutagenesis and combinatorial search. Traditional enzyme engineering approaches have a limited scope for complex optimization due to the requirement of a priori knowledge or experimental burden of screening huge protein libraries. The recent surge in high-throughput experimental methods including Next Generation Sequencing and automated screening has flooded the field of molecular biology with big-data, which requires us to re-think our concurrent approaches towards enzyme engineering. Artificial Intelligence (AI) and Machine Learning (ML) have great potential to revolutionize smart enzyme engineering without the explicit need for a complete understanding of the underlying molecular system. Here, we portray the role and position of AI techniques in the field of enzyme engineering along with their scope and limitations. In addition, we explain how the traditional approaches of directed evolution and rational design can be extended through AI tools. Recent successful examples of AI-assisted enzyme engineering projects and their deviation from traditional approaches are highlighted. A comprehensive picture of current challenges and future avenues for AI in enzyme engineering are also discussed.

Improving the Thermostability of Xylanase A from Bacillus subtilis by Combining Bioinformatics and Electrostatic Interactions Optimization

The rational improvement of the enzyme catalytic activity is one of the most significant challenges in biotechnology. Most conventional strategies used to engineer enzymes involve selecting mutations to increase their thermostability. Determining good criteria for choosing these substitutions continues to be a challenge. In this work, we combine bioinformatics, electrostatic analysis, and molecular dynamics to predict beneficial mutations that may improve the thermostability of XynA from Bacillus subtilis. First, the Tanford-Kirkwood surface accessibility method is used to characterize each ionizable residue contribution to the protein native state stability. Residues identified to be destabilizing were mutated with the corresponding residues determined by the consensus or ancestral sequences at the same locations. Five mutants (K99T/N151D, K99T, S31R, N151D, and K154A) were investigated and compared with 12 control mutants derived from experimental approaches from the literature. Molecular dynamics results show that the mutants exhibited folding temperatures in the order K99T > K99T/N151D > S31R > N151D > WT > K154A. The combined approaches employed provide an effective strategy for low-cost enzyme optimization needed for large-scale biotechnological and medical applications.

Improving the Thermostability and Activity of Transaminase From Aspergillus terreus by Charge-Charge Interaction

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 (t _1/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 t _1/2 at 40°C, and 6.08°C increase in T 50 10 . 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.

Membrane Protein Stabilization Strategies for Structural and Functional Studies

Accounting for nearly two-thirds of known druggable targets, membrane proteins are highly relevant for cell physiology and pharmacology. In this regard, the structural determination of pharmacologically relevant targets would facilitate the intelligent design of new drugs. The structural biology of membrane proteins is a field experiencing significant growth as a result of the development of new strategies for structure determination. However, membrane protein preparation for structural studies continues to be a limiting step in many cases due to the inherent instability of these molecules in non-native membrane environments. This review describes the approaches that have been developed to improve membrane protein stability. Membrane protein mutagenesis, detergent selection, lipid membrane mimics, antibodies, and ligands are described in this review as approaches to facilitate the production of purified and stable membrane proteins of interest for structural and functional studies.

Engineering a thermostable chondroitinase for production of specifically distributed low-molecular-weight chondroitin sulfate

Chondroitinase ABC I (csABC I) has attracted intensive attention because of its great potential in heparin refining and the enzymatic preparation of low-molecular-weight chondroitin sulfate (LMW-CS). However, low thermal resistance (<30℃) restricts its applications. Herein, structure-guided and sequence-assisted combinatorial engineering approaches were applied to improve the thermal resistance of Proteus vulgaris csABC I. By integrating the deletion of the flexible fragment R166-L170 at the N-terminal domain and the mutation of E694P at the C-terminal domain, variant NΔ5/E694P exhibited 247-fold improvement of its half-life at 37℃ and a 2.3-fold increase in the specific activity. Through batch fermentation in a 3-L fermenter, the expression of variant NΔ5/E694P in an Escherichia coli host reached 1.7 g L^-1 with the activity of 1.0 × 10⁵ U L^-1 . Finally, the enzymatic approach for the preparation of LMW-CS was established. By modulating enzyme concentration and controlling depolymerization time, specifically distributed LMW-CS (7000, 3400, and 1900 Da) with low polydispersity was produced, demonstrating the applicability of these processes for the industrial production of LMW-CS in a more environmentally friendly way.

Development of a structure-analysis pipeline using multiple-solvent crystal structures of barrier-to-autointegration factor

The multiple-solvent crystal structure (MSCS) approach uses high concentrations of organic solvents to characterize the interactions and effects of solvents on proteins. Here, the method has been further developed and an MSCS data-handling pipeline is presented that uses the Detection of Related Solvent Positions (DRoP) program to improve data quality. DRoP is used to selectively model conserved water molecules, so that an advanced stage of structural refinement is reached quickly. This allows the placement of organic molecules more accurately and convergence on high-quality maps and structures. This pipeline was applied to the chromatin-associated protein barrier-to-autointegration factor (BAF), resulting in structural models with better than average statistics. DRoP and Phenix Structure Comparison were used to characterize the data sets and to identify a binding site that overlaps with the interaction site of BAF with emerin. The conserved water-mediated networks identified by DRoP suggested a mechanism by which water molecules are used to drive the binding of DNA. Normalized and differential B-factor analysis is shown to be a valuable tool to characterize the effects of specific solvents on defined regions of BAF. Specific solvents are identified that cause stabilization of functionally important regions of the protein. This work presents tools and a standardized approach for the analysis and comprehension of MSCS data sets.

Protein Sequence Selection Method That Enables Full Consensus Design of Artificial l-Threonine 3-Dehydrogenases with Unique Enzymatic Properties

Exponentially increasing protein sequence data enables artificial enzyme design using sequence-based protein design methods, including full-consensus protein design (FCD). The success of artificial enzyme design is strongly dependent on the nature of the sequences used. Hence, sequences must be selected from databases and curated libraries prepared to enable a successful design by FCD. In this study, we proposed a selection approach regarding several key residues as sequence motifs. We used l-threonine 3-dehydrogenase (TDH) as a model to test the validity of this approach. In the classification, four residues (143, 174, 188, and 214) were used as key residues. We classified thousands of TDH homologous sequences into five groups containing hundreds of sequences. Utilizing sequences in the libraries, we designed five artificial TDHs by FCD. Among the five, we successfully expressed four in soluble form. Biochemical analysis of artificial TDHs indicated that their enzymatic properties vary; half of the maximum measured enzyme activity (t_1/2) and activation energies were distributed from 53 to 65 °C and from 38 to 125 kJ/mol, respectively. The artificial TDHs had unique kinetic parameters, distinct from one another. Structural analysis indicates that consensus mutations are mainly introduced in the secondary or outer shell. The functional diversity of the artificial TDHs is due to the accumulation of mutations that affect their physicochemical properties. Taken together, our findings indicate that our proposed approach can help generate artificial enzymes with unique enzymatic properties.

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.

Improving the catalytic thermostability of Bacillus altitudinis W3 ω-transaminase by proline substitutions

As a green biocatalyst, transaminase with high thermostability can be better employed to synthesize many pharmaceutical intermediates in industry. To improve the thermostability of (R)-selective amine transaminase from Bacillus altitudinis W3, related mutation sites were determined by multiple amino acid sequence alignment between wild-type ω-transaminase and four potential thermophilic ω-transaminases, followed by replacement of the related amino acid residues with proline by site-directed mutagenesis. Three stabilized mutants (D192P, T237P, and D192P/T237P) showing the highest stability were obtained and used for further analysis. Comparison with the wild-type enzyme revealed that the double mutant D192P/T237P exhibited the largest shift in thermostability, with a 2.5-fold improvement of t _1/2 at 40 °C, and a 6.3 °C increase in T ₅₀ ¹⁵, and a 5 °C higher optimal catalytic temperature. Additionally, this mutant exhibited an increase in catalytic efficiency (k _cat/K _m) relative to the wild-type enzyme. Modeling analysis indicated that the improved thermostability of the mutants could be associated with newly formed hydrophobic interactions and hydrogen bonds. This study shown that proline substitutions guided by sequence alignment to improve the thermostability of (R)-selective amine transaminase was effective and this method can also be used to engineering other enzymes.

Larger active site in an ancestral hydroxynitrile lyase increases catalytically promiscuous esterase activity

Hydroxynitrile lyases (HNL's) belonging to the α/β-hydrolase-fold superfamily evolved from esterases approximately 100 million years ago. Reconstruction of an ancestral hydroxynitrile lyase in the α/β-hydrolase fold superfamily yielded a catalytically active hydroxynitrile lyase, HNL1. Several properties of HNL1 differ from the modern HNL from rubber tree (HbHNL). HNL1 favors larger substrates as compared to HbHNL, is two-fold more catalytically promiscuous for ester hydrolysis (p-nitrophenyl acetate) as compared to mandelonitrile cleavage, and resists irreversible heat inactivation to 35 °C higher than for HbHNL. We hypothesized that the x-ray crystal structure of HNL1 may reveal the molecular basis for the differences in these properties. The x-ray crystal structure solved to 1.96-Å resolution shows the expected α/β-hydrolase fold, but a 60% larger active site as compared to HbHNL. This larger active site echoes its evolution from esterases since related esterase SABP2 from tobacco also has a 38% larger active site than HbHNL. The larger active site in HNL1 likely accounts for its ability to accept larger hydroxynitrile substrates. Site-directed mutagenesis of HbHNL to expand the active site increased its promiscuous esterase activity 50-fold, consistent with the larger active site in HNL1 being the primary cause of its promiscuous esterase activity. Urea-induced unfolding of HNL1 indicates that it unfolds less completely than HbHNL (m-value = 0.63 for HNL1 vs 0.93 kcal/mol·M for HbHNL), which may account for the ability of HNL1 to better resist irreversible inactivation upon heating. The structure of HNL1 shows changes in hydrogen bond networks that may stabilize regions of the folded structure.

Consensus Design of an Evolved High-Redox Potential Laccase

Among the broad repertory of protein engineering methods that set out to improve stability, consensus design has proved to be a powerful strategy to stabilize enzymes without compromising their catalytic activity. Here, we have applied an in-house consensus method to stabilize a laboratory evolved high-redox potential laccase. Multiple sequence alignments were carried out and computationally refined by applying relative entropy and mutual information thresholds. Through this approach, an ensemble of 20 consensus mutations were identified, 18 of which were consensus/ancestral mutations. The set of consensus variants was produced in Saccharomyces cerevisiae and analyzed individually, while site directed recombination of the best mutations did not produce positive epistasis. The best single variant carried the consensus-ancestral A240G mutation in the neighborhood of the T2/T3 copper cluster, which dramatically improved thermostability, kinetic parameters and secretion.

Improving the thermostability and activity of Paenibacillus pasadenensis chitinase through semi-rational design

Chitinase is a promising biocatalyst for chitin biotransformation in the field of recalcitrant biomass degradation. Excellent catalytic performance is conducive to its commercial utilization. In this work, sequence- and structure-based semi-rational design was performed to evolve the thermostability and activity of a previously identified chitinase PpChi1 from Paenibacillus pasadenensis CS0611. After combinational mutagenesis, the mutant S244C-I319C/T259P with disulfide bond introduction and proline substitution exhibited higher specific activity at higher temperature, 26.3-fold in half-life value at 50 °C, and a 7.9 °C rise in half-inactivation temperature T_1/2^15min compared to the wild-type enzyme. The optimal reaction temperature of the mutant was shifted from 45 °C to 52.5 °C. Molecular dynamic simulation and structure analysis confirmed that these improvements of the mutant were attributed to its stabilized folding form, possibly caused by the decreased entropy of unfolding. This work gives an initial insight into the effect of conserved proline residues in thermostable chitinases and proposes a feasible approach for improving chitinase thermostability to facilitate its application in chitin hydrolysis to valuable oligosaccharides.

High-Throughput, Lysis-Free Screening for Sulfatase Activity Using Escherichia coli Autodisplay in Microdroplets

Directed evolution of enzymes toward improved catalytic performance has become a powerful tool in protein engineering. To be effective, a directed evolution campaign requires the use of high-throughput screening. In this study we describe the development of an ultra high-throughput lysis-free procedure to screen for improved sulfatase activity by combining microdroplet-based single-variant activity sorting with E. coli autodisplay. For the first step in a 4-step screening procedure, we quantitatively screened >10⁵ variants of the homodimeric arylsulfatase from Silicibacter pomeroyi (SpAS1), displayed on the E. coli cell surface, for improved sulfatase activity using fluorescence activated droplet sorting. Compartmentalization of the fluorescent reaction product with living E. coli cells autodisplaying the sulfatase variants ensured the continuous linkage of genotype and phenotype during droplet sorting and allowed for direct recovery by simple regrowth of the sorted cells. The use of autodisplay on living cells simplified and reduced the degree of liquid handling during all steps in the screening procedure to the single event of simply mixing substrate and cells. The percentage of apparent improved variants was enriched >10-fold as a result of droplet sorting. We ultimately identified 25 SpAS1 variants with improved performance toward 4-nitrophenyl sulfate (up to 6.2-fold) and/or fluorescein disulfate (up to 30-fold). In SpAS1 variants with improved performance toward the bulky fluorescein disulfate, many of the beneficial mutations occur in residues that form hydrogen bonds between α-helices in the C-terminal oligomerization region, suggesting a previously unknown role for the dimer interface in shaping the substrate binding site of SpAS1.

Construction of stabilized (R)-selective amine transaminase from Aspergillus terreus by consensus mutagenesis

Amine transaminases are a class of efficient and industrially-desired biocatalysts for the production of chiral amines. In this study, stabilized variants of the (R)-selective amine transaminase from Aspergillus terreus (AT-ATA) were constructed by consensus mutagenesis. Using Consensus Finder (http://cbs-kazlab.oit.umn.edu/), six positions with the most prevalent amino acid (over 60% threshold) among the homologous family members were identified. Subsequently, these six residues were individually mutated to match the consensus sequence (I77 L, Q97E, H210N, N245D, G292D, and I295 V) using site-directed mutagenesis. Compared to that of the wild-type, the thermostability of all six single variants was improved. The H210N variant displayed the largest shift in thermostability, with a 3.3-fold increase in half-life (t_1/2) at 40 °C, and a 4.6 °C increase in T₅₀¹⁰ among the single variants. In addition, the double mutant H210N/I77L displayed an even larger shift with 6.1-fold improvement of t_1/2 at 40 °C, and a 6.6 °C increase in T₅₀¹⁰. Furtherly, the H210N/I77L mutation was introduced into the previously engineered thermostable AT-ATA by the introduction of disulfide bonds, employing B-factor and folding free energy (ΔΔG^fold) calculations. Our results showed that the combined variant H210N/I77L/M150C-M280C had the largest shift in thermostability, with a 16.6-fold improvement of t_1/2 and a 11.8 °C higher T₅₀¹⁰.

"Site and Mutation"-Specific Predictions Enable Minimal Directed Evolution Libraries

Directed evolution experiments designed to improve the activity of a biocatalyst have increased in sophistication from the early days of completely random mutagenesis. Sequence-based and structure-based methods have been developed to identify "hotspot" positions that when randomized provide a higher frequency of beneficial mutations that improve activity. These focused mutagenesis methods reduce library sizes and therefore reduce screening burden, accelerating the rate of finding improved enzymes. Looking for further acceleration in finding improved enzymes, we investigated whether two existing methods, one sequence-based (Protein GPS) and one structure-based (using Bioluminate and MOE), were sufficiently predictive to provide not just the hotspot position, but also the amino acid substitution that improved activity at that position. By limiting the libraries to variants that contained only specific amino acid substitutions, library sizes were kept to less than 100 variants. For an initial round of ATA-117 R-selective transaminase evolution, we found that the methods used produced libraries where 9% and 18% of the amino acid substitutions chosen were amino acids that improved reaction performance in lysates. The ability to create combinations of mutations as part of the initial design was confounded by the relatively large number of predicted mutations that were inactivating (30% and 45% for the sequence-based and structure-based methods, respectively). Despite this, combining several mutations identified within a given method produced variant lysates 7- and 9-fold more active than the wild-type lysate, highlighting the capability of mutations chosen this way to generate large advances in activity in addition to the reductions in screening.