Protein stability is crucial in enzymatic catalysis. To improve the efficiency in the searching for thermostablizing mutations, we applied a sequence consensus approach focusing on dimeric interface residues of ketoreductase ChKRED20. The strategy returned a success rate of 43%, revealing 9 beneficial mutations from 21 candidates with improved kinetic or thermodynamic stability. Several combinatorial mutants were then constructed, and mutant M8K displayed the highest thermostability, with a melting temperature (Tm) of 89 °C and a half-inactivation temperature (T50) of 93.4 °C, both of over 35 °C increase compared to the wild-type. M8K could remain stable for at least 7 days at its optimal reaction temperature of 55 °C. Its inactivation half-life (t1/2) was 110 min at 90 °C, while the wild-type was 18.6 min at 60 °C. The results were interpreted in the context of structural and molecular dynamic simulation analysis, which revealed the addition of intramolecular interactions, decreased conformational flexibility and increased compactness, all in agreement with the observed effect.

t-Butyl 6-cyano-(3R,5R)-dihydroxyhexanoate ((3R,5R)-2) is an advanced chiral diol intermediate of the cholesterol-lowering drug atorvastatin. KmAKRM5 (W297H/Y296W/K29H/Y28A/T63M) constructed in our previous work, displayed good biocatalytic performance on (3R,5R)-2. In the present work, stepwise evolution was applied to further enhance the thermostability and activity of KmAKRM5. For thermostability enhancement, N109 and S196 located far from the active site were picked out by structure-guided consensus engineering, and mutated by site-directed mutagenesis (SDM). For catalytic efficiency improvement, the residues A30 and T302 adjacent to the substrate-binding pocket were subjected to site-saturation mutagenesis (SSM). As a result, the \"best\" mutant KmAKRM9 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C) was developed, of which T5015 and Tm were 5.0 °C and 8.2 °C higher than those of KmAKRM5. Moreover, compared to KmAKRM5, KmAKRM9 displayed a 1.9-fold (846 vs 2436 min) and 6.7-fold (126 vs 972 min) longer half-lives at 40 and 50 °C, respectively. Structural analysis suggested that beneficial mutations introduced additional hydrophobic interactions and hydrogen bonds, contributing rigidification of the flexible loops and the increase of internal forces, hence increasing the thermostability and activity. 5 g DCW (dry cell weight) L-1KmAKRM9 completely reduced 350 g L-1t-butyl 6-cyano-(5R)-hydroxy-3-oxo-hexanoate ((5R)-1), within 3.7 h at 40 °C, yielding optically pure (3R,5R)-2 (d.e.p > 99.5%) with a space-time yield (STY) of 1.82 kg L-1 d-1. Hence, KmAKRM9 is a robust biocatalyst for the synthesis of (3R,5R)-2.

BACKGROUND: Enzymatic quantification of creatinine has become an essential method for clinical evaluation of renal function. Although creatinase (CR) is frequently used for this purpose, its poor thermostability severely limits industrial applications. Herein, we report a novel creatinase from Alcaligenes faecalis (afCR) with higher catalytic activity and lower KM value, than currently used creatinases. Furthermore, we developed a non-biased phylogenetic consensus method to improve the thermostability of afCR.
RESULTS: We applied a non-biased phylogenetic consensus method to identify 59 candidate consensus residues from 24 creatinase family homologs for screening afCR mutants with improved thermostability. Twenty-one amino acids of afCR were selected to mutagenesis and 11 of them exhibited improved thermostability compared to the parent enzyme (afCR-M0). Combination of single-site mutations in sequential screens resulted in a quadruple mutant D17V/T199S/L6P/T251C (M4-2) which showed ~ 1700-fold enhanced half-life at 57 °C and a 4.2 °C higher T5015 than that of afCR-M0. The mutant retained catalytic activity equivalent to afCR-M0, and thus showed strong promise for application in creatinine detection. Structural homology modeling revealed a wide range of potential molecular interactions associated with individual mutations that contributed to improving afCR thermostability.
CONCLUSIONS: Results of this study clearly demonstrated that the non-biased-phylogenetic consensus design for improvement of thermostability in afCR is effective and promising in improving the thermostability of more enzymes.

γ-Aminobutyrate (GABA) is an important bioactive compound synthesized through decarboxylation of L-glutamate by the glutamate decarboxylase (GAD). In this study, stabilized variants of the GAD from Lactobacillus brevis were constructed by consensus mutagenesis. Using Consensus Finder ( http://cbs-kazlab.oit.umn.edu/ ), eight positions with the most prevalent amino acid (over 60% threshold) among the homologous family members were identified. Subsequently, these eight residues were individually mutated to match the consensus sequence using site-directed mutagenesis. Compared to the wild-type, T383K variant displayed the largest shift in thermostability among the single variants, with a 3.0 °C increase in semi-inactivation temperature (T5015), a 1.7-fold improvement of half-life (t1/2) at 55 °C, and a 1.2-fold improvement of t1/2 at 37 °C, respectively, while its catalytic efficiency (kcat/Km) was reduced. To obtain the mutant with improvement in both thermostability and catalytic activity, we performed a site-saturation mutation at T383. Notably, mutants T383V and T383G exhibited an increasement in thermostability and kcat/Km than that of wild-type. This study not only emphasizes the value of consensus mutagenesis for improving the thermostability of GAD but also sheds a powerful guidance to study the thermal stability of other enzymes.

Membrane proteins are generally challenging to work with because of their notorious instability. Protein engineering has been used increasingly to thermostabilize labile membrane proteins such as G-protein-coupled receptors for structural and functional studies in recent years. Two major strategies exist. Scanning mutagenesis systematically eliminates destabilizing residues, whereas the consensus approach assembles mutants with the most frequent residues among selected homologs, bridging sequence conservation with stability. Here, we applied the consensus concept to stabilize a fungal homolog of the human sterol Δ8-7 isomerase, a 26.4 kDa protein with five transmembrane helices. The isomerase is also called emopamil-binding protein (EBP), as it binds this anti-ischemic drug with high affinity. The wild-type had an apparent melting temperature (Tm) of 35.9 °C as measured by the fluorescence-detection size-exclusion chromatography-based thermostability assay. A total of 87 consensus mutations sourced from 22 homologs gained expression level and thermostability, increasing the apparent Tm to 69.9 °C at the cost of partial function loss. Assessing the stability and activity of several systematic chimeric constructs identified a construct with an apparent Tm of 79.8 °C and two regions for function rescue. Further back-mutations of the chimeric construct in the two target regions yielded the final construct with similar apparent activity to the wild-type and an elevated Tm of 88.8 °C, totaling an increase of 52.9 °C. The consensus approach is effective and efficient because it involves fewer constructs compared with scanning mutagenesis. Our results should encourage more use of the consensus strategy for membrane protein thermostabilization.

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 (t1/2) at 40 °C, and a 4.6 °C increase in T5010 among the single variants. In addition, the double mutant H210N/I77L displayed an even larger shift with 6.1-fold improvement of t1/2 at 40 °C, and a 6.6 °C increase in T5010. 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 (ΔΔGfold) 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 t1/2 and a 11.8 °C higher T5010.

Lipoxygenase (LOX) from Anabaena sp. PCC 7120 (Ana-rLOX) offers important applications in the food industry, especially for improving aroma and dough rheological properties. However, industrial applications of LOXs have been limited by their poor thermostability. Herein, we report a bioinformatics-based consensus concept approach for the engineering of thermostable Ana-rLOX.
A series of mutations (N130D, G260A, S437T, N130D/G260Q, N130D/S437Y) showed higher thermostability and activity than the wild-type enzyme. Thus, N130D/G260Q exhibited a 6.6-fold increase in half-life and 2.45 °C increase in unfolding temperature; N130D/S437Y showed a 10 °C increase in optimal temperature. The secondary structure did not change much that contributed to improved thermostability were investigated in detail using circular dichroism. Homology modeling suggested that enhanced thermostability and specific activity may result from favorable hydrophobic interactions.
A series of mutations were achieved, showing higher thermostability and activity than the wild-type enzyme by semi-rational mutagenesis with limited structure information. Our findings provide important new insights into molecular modifications aimed at improving Ana-rLOX thermostability and activity.

Aldehyde-deformylating oxygenase (ADO) is an essential enzyme for production of long-chain alkanes as drop-in biofuels, which are compatible with existing fuel systems. The most active ADOs are present in mesophilic cyanobacteria, especially Nostoc punctiforme Given the potential applications of thermostable enzymes in biorefineries, here we generated a thermostable (Cts)-ADO based on a consensus of ADO sequences from several thermophilic cyanobacterial strains. Using an in silico design pipeline and a metagenome library containing 41 hot-spring microbial communities, we created Cts-ADO. Cts-ADO displayed a 3.8-fold increase in pentadecane production on raising the temperature from 30 to 42 °C, whereas ADO from N. punctiforme (Np-ADO) exhibited a 1.7-fold decline. 3D structure modeling and molecular dynamics simulations of Cts- and Np-ADO at different temperatures revealed differences between the two enzymes in residues clustered on exposed loops of these variants, which affected the conformation of helices involved in forming the ADO catalytic core. In Cts-ADO, this conformational change promoted ligand binding to its preferred iron, Fe2, in the di-iron cluster at higher temperature, but the reverse was observed in Np-ADO. Detailed mapping of residues conferring Cts-ADO thermostability identified four amino acids, which we substituted individually and together in Np-ADO. Among these substitution variants, A161E was remarkably similar to Cts-ADO in terms of activity optima, kinetic parameters, and structure at higher temperature. A161E was located in loop L6, which connects helices H5 and H6, and supported ligand binding to Fe2 at higher temperatures, thereby promoting optimal activity at these temperatures and explaining the increased thermostability of Cts-ADO.

The multidomain, catalytically self-sufficient cytochrome P450 BM-3 from Bacillus megaterium (P450BM3 ) constitutes a versatile enzyme for the oxyfunctionalization of organic molecules and natural products. However, the limited stability of the diflavin reductase domain limits the utility of this enzyme for synthetic applications. In this work, a consensus-guided mutagenesis approach was applied to enhance the thermal stability of the reductase domain of P450BM3 . Upon phylogenetic analysis of a set of distantly related P450s (>38 % identity), a total of 14 amino acid substitutions were identified and evaluated in terms of their stabilizing effects relative to the wild-type reductase domain. Recombination of the six most stabilizing mutations generated two thermostable variants featuring up to tenfold longer half-lives at 50 °C and increased catalytic performance at elevated temperatures. Further characterization of the engineered P450BM3 variants indicated that the introduced mutations increased the thermal stability of the FAD-binding domain and that the optimal temperature (Topt ) of the enzyme had shifted from 25 to 40 °C. This work demonstrates the effectiveness of consensus mutagenesis for enhancing the stability of the reductase component of a multidomain P450. The stabilized P450BM3 variants developed here could potentially provide more robust scaffolds for the engineering of oxidation biocatalysts.

The UPF0016 family is a recently identified group of poorly characterized membrane proteins whose function is conserved through evolution and that are defined by the presence of 1 or 2 copies of the E-φ-G-D-[KR]-[TS] consensus motif in their transmembrane domain. We showed that 2 members of this family, the human TMEM165 and the budding yeast Gdt1p, are functionally related and are likely to form a new group of Ca2+ transporters. Mutations in TMEM165 have been demonstrated to cause a new type of rare human genetic diseases denominated as Congenital Disorders of Glycosylation. Using site-directed mutagenesis, we generated 17 mutations in the yeast Golgi-localized Ca2+ transporter Gdt1p. Single alanine substitutions were targeted to the highly conserved consensus motifs, 4 acidic residues localized in the central cytosolic loop, and the arginine at position 71. The mutants were screened in a yeast strain devoid of both the endogenous Gdt1p exchanger and Pmr1p, the Ca2+ -ATPase of the Golgi apparatus. We show here that acidic and polar uncharged residues of the consensus motifs play a crucial role in calcium tolerance and calcium transport activity and are therefore likely to be architectural components of the cation binding site of Gdt1p. Importantly, we confirm the essential role of the E53 residue whose mutation in humans triggers congenital disorders of glycosylation.