FireProt: Energy- and Evolution-Based Computational Design of Thermostable Multiple-Point Mutants

There is great interest in increasing proteins’ stability to enhance their utility as biocatalysts, therapeutics, diagnostics and nanomaterials. Directed evolution is a powerful, but experimentally strenuous approach. Computational methods offer attractive alternatives. However, due to the limited reliability of predictions and potentially antagonistic effects of substitutions, only single-point mutations are usually predicted in silico, experimentally verified and then recombined in multiple-point mutants. Thus, substantial screening is still required. Here we present FireProt, a robust computational strategy for predicting highly stable multiple-point mutants that combines energy- and evolution-based approaches with smart filtering to identify additive stabilizing mutations. FireProt’s reliability and applicability was demonstrated by validating its predictions against 656 mutations from the ProTherm database. We demonstrate that thermostability of the model enzymes haloalkane dehalogenase DhaA and γ-hexachlorocyclohexane dehydrochlorinase LinA can be substantially increased (ΔT _m = 24°C and 21°C) by constructing and characterizing only a handful of multiple-point mutants. FireProt can be applied to any protein for which a tertiary structure and homologous sequences are available, and will facilitate the rapid development of robust proteins for biomedical and biotechnological applications.     

Improvement in Thermostability of an Achaetomium sp. Strain Xz8 Endopolygalacturonase via the Optimization of Charge-Charge Interactions

Improving enzyme thermostability is of importance for widening the spectrum of application of enzymes. In this study, a structure-based rational design approach was used to improve the thermostability of a highly active, wide-pH-range-adaptable, and stable endopolygalacturonase (PG8fn) from Achaetomium sp. strain Xz8 via the optimization of charge-charge interactions. By using the enzyme thermal stability system (ETSS), two residues—D244 and D299—were inferred to be crucial contributors to thermostability. Single (D244A and D299R) and double (D244A/D299R) mutants were then generated and compared with the wild type. All mutants showed improved thermal properties, in the order D244A < D299R < D244A/D299R. In comparison with PG8fn, D244A/D299R showed the most pronounced shifts in temperature of maximum enzymatic activity (T_max), temperature at which 50% of the maximal activity of an enzyme is retained (T₅₀), and melting temperature (T_m), of about 10, 17, and 10.2°C upward, respectively, with the half-life (t_1/2) extended by 8.4 h at 50°C and 45 min at 55°C. Another distinguishing characteristic of the D244A/D299R mutant was its catalytic activity, which was comparable to that of the wild type (23,000 ± 130 U/mg versus 28,000 ± 293 U/mg); on the other hand, it showed more residual activity (8,400 ± 83 U/mg versus 1,400 ± 57 U/mg) after the feed pelleting process (80°C and 30 min). Molecular dynamics (MD) simulation studies indicated that mutations at sites D244 and D299 lowered the overall root mean square deviation (RMSD) and consequently increased the protein rigidity. This study reveals the importance of charge-charge interactions in protein conformation and provides a viable strategy for enhancing protein stability.     

Adopting Selected Hydrogen Bonding and Ionic Interactions from Aspergillus fumigatus Phytase Structure Improves the Thermostability of Aspergillus niger PhyA Phytase

Although it has been widely used as a feed supplement to reduce manure phosphorus pollution of swine and poultry, Aspergillus niger PhyA phytase is unable to withstand heat inactivation during feed pelleting. Crystal structure comparisons with its close homolog, the thermostable Aspergillus fumigatus phytase (Afp), suggest associations of thermostability with several key residues (E35, S42, R168, and R248) that form a hydrogen bond network in the E35-to-S42 region and ionic interactions between R168 and D161 and between R248 and D244. In this study, loss-of-function mutations (E35A, R168A, and R248A) were introduced singularly or in combination into seven mutants of Afp. All seven mutants displayed decreases in thermostability, with the highest loss (25% [P < 0.05]) in the triple mutant (E35A R168A R248A). Subsequently, a set of corresponding substitutions were introduced into nine mutants of PhyA to strengthen the hydrogen bonding and ionic interactions. While four mutants showed improved thermostability, the best response came from the quadruple mutant (A58E P65S Q191R T271R), which retained 20% greater (P < 0.05) activity after being heated at 80°C for 10 min and had a 7°C higher melting temperature than that of wild-type PhyA. This study demonstrates the functional importance of the hydrogen bond network and ionic interaction in supporting the high thermostability of Afp and the feasibility of adopting these structural units to improve the thermostability of a homologous PhyA phytase.     

Mutations in salt-bridging residues at the interface of the core and lid domains of epoxide hydrolase StEH1 affect regioselectivity, protein stability and hysteresis

Epoxide hydrolase, StEH1, shows hysteretic behavior in the catalyzed hydrolysis of trans-2-methylstyrene oxide (2-MeSO)¹. Linkage between protein structure dynamics and catalytic function was probed in mutant enzymes in which surface-located salt-bridging residues were substituted. Salt-bridges at the interface of the α/β-hydrolase fold core and lid domains, as well as between residues in the lid domain, between Lys¹⁷⁹-Asp²⁰², Glu²¹⁵-Arg⁴¹ and Arg²³⁶-Glu¹⁶⁵ were disrupted by mutations, K179Q, E215Q, R236K and R236Q. All mutants displayed enzyme activity with styrene oxide (SO) and 2-MeSO when assayed at 30 °C. Disruption of salt-bridges altered the rates for isomerization between distinct Michaelis complexes, with (1R,2R)-2-MeSO as substrate, presumably as a result of increased dynamics of involved protein segments. Another indication of increased flexibility was a lowered thermostability in all mutants. We propose that the alterations to regioselectivity in these mutants derive from an increased mobility in protein segments otherwise stabilized by salt bridging interactions.     

Thermostability Improvement of a Streptomyces Xylanase by Introducing Proline and Glutamic Acid Residues

Protein engineering is commonly used to improve the robustness of enzymes for activity and stability at high temperatures. In this study, we identified four residues expected to affect the thermostability of Streptomyces sp. strain S9 xylanase XynAS9 through multiple-sequence analysis (MSA) and molecular dynamic simulations (MDS). Site-directed mutagenesis was employed to construct five mutants by replacing these residues with proline or glutamic acid (V81P, G82E, V81P/G82E, D185P/S186E, and V81P/G82E/D185P/S186E), and the mutant and wild-type enzymes were expressed in Pichia pastoris. Compared to the wild-type XynAS9, all five mutant enzymes showed improved thermal properties. The activity and stability assays, including circular dichroism and differential scanning calorimetry, showed that the mutations at positions 81 and 82 increased the thermal performance more than the mutations at positions 185 and 186. The mutants with combined substitutions (V81P/G82E and V81P/G82E/D185P/S186E) showed the most pronounced shifts in temperature optima, about 17°C upward, and their half-lives for thermal inactivation at 70°C and melting temperatures were increased by >9 times and approximately 7.0°C, respectively. The mutation combination of V81P and G82E in adjacent positions more than doubled the effect of single mutations. Both mutation regions were at the end of long secondary-structure elements and probably rigidified the local structure. MDS indicated that a long loop region after positions 81 and 82 located at the end of the inner β-barrel was prone to unfold. The rigidified main chain and filling of a groove by the mutations on the bottom of the active site canyon may stabilize the mutants and thus improve their thermostability.     

Identification of a new subfamily of salt-tolerant esterases from a metagenomic library of tidal flat sediment

To search for novel lipolytic enzymes, a metagenomic library was constructed from the tidal flat sediment of Ganghwa Island in South Korea. By functional screening using tributyrin agar plates, 3 clones were selected from among the 80,050 clones of the fosmid library. The sequence analysis revealed that those clones contained different open reading frames, which showed 50–57% amino acid identity with putative lipolytic enzymes in the database. Based on the phylogenetic analysis, they were identified to encode novel members, which form a distinct and new subfamily in the family IV of bacterial lipolytic enzymes. The consensus sequence, GT(S)SA(G)G, encompassing the active site serine of the enzymes was different from the GDSAG motif, conserved in the other subfamily. The genes were expressed in Escherichia coli and recombinant proteins were purified as active soluble forms. The enzymes showed the highest activity toward p-nitrophenyl valerate (C5) and exhibited optimum activities at mesophilic temperature ranges and slightly alkaline pH. In particular, the enzymes displayed salt tolerance with over 50% of the maximum activity remained in the presence of 3 M NaCl (or KCl). In this study, we demonstrated that the metagenomic approach using marine tidal flat sediment as a DNA source expanded the diversity of lipolytic enzyme-encoding genes.     

Effect of surface electrostatic interactions on the stability and folding of formate dehydrogenase from Candida methylica

NAD⁺-dependent formate dehydrogenase (FDH-EC 1.2.1.2) is an important enzyme to regenerate valuable NADH required by NAD⁺-dependent oxidoreductases in enzyme catalysis. The limitation in the thermostability of FDH enzyme is a crucial problem for development of biotechnological and industrial processes, despite of its advantages. In this study, to investigate the contribution of surface electrostatic interaction to the thermostability of FDH from Candida methylica (cmFDH) N187E, H13E, Q105R, N300E, N147R N300E/N147R, N187E/Q105R, N187E/N147R,Y160R, Y302R, Y160E and Y302E mutants were designed using a homology model of cmFDH based on Candida boidinii (cb) by considering electrostatic interactions on the protein surface. The effects of site-specific engineering on the stability of this molecule was analyzed according to minimal model of folding and assembly reaction and deduced equilibrium properties of the native system with respect to its thermal and denaturant sensitivities. It was observed that mutations did not change the unfolding pattern of native cmFDH and increased numbers of electrostatic interactions can cause either stabilizing or destabilizing effect on the thermostability of this protein. The thermodynamic and kinetic results suggested that except relatively improved mutants, three out of the nine single mutations increased the melting temperature of cmFDH enzyme.     

In Silico Rational Design and Systems Engineering of Disulfide Bridges in the Catalytic Domain of an Alkaline α-Amylase from Alkalimonas amylolytica To Improve Thermostability

High thermostability is required for alkaline α-amylases to maintain high catalytic activity under the harsh conditions used in textile production. In this study, we attempted to improve the thermostability of an alkaline α-amylase from Alkalimonas amylolytica through in silico rational design and systems engineering of disulfide bridges in the catalytic domain. Specifically, 7 residue pairs (P35-G426, Q107-G167, G116-Q120, A147-W160, G233-V265, A332-G370, and R436-M480) were chosen as engineering targets for disulfide bridge formation, and the respective residues were replaced with cysteines. Three single disulfide bridge mutants—P35C-G426C, G116C-Q120C, and R436C-M480C—of the 7 showed significantly enhanced thermostability. Combinational mutations were subsequently assessed, and the triple mutant P35C-G426C/G116C-Q120C/R436C-M480C showed a 6-fold increase in half-life at 60°C and a 5.2°C increase in melting temperature compared with the wild-type enzyme. Interestingly, other biochemical properties of this mutant also improved: the optimum temperature increased from 50°C to 55°C, the optimum pH shifted from 9.5 to 10.0, the stable pH range extended from 7.0 to 11.0 to 6.0 to 12.0, and the catalytic efficiency (k_cat/K_m) increased from 1.8 × 10⁴ to 2.4 × 10⁴ liters/g · min. The possible mechanism responsible for these improvements was explored through comparative analysis of the model structures of wild-type and mutant enzymes. The disulfide bridge engineering strategy used in this work may be applied to improve the thermostability of other industrial enzymes.     

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.     

Improved thermostability of a metagenomic glucose-tolerant β-glycosidase based on its X-ray crystal structure

MeBglD2, a metagenomic β-glycosidase, is stimulated by various saccharides, including d-glucose, d-xylose, and maltose, and it promotes the enzymatic saccharification of plant biomass. To improve the thermostability of MeBglD2, its X-ray crystal structure was analyzed, and the amino acid residues responsible for its thermostability were identified using the structural information. Mutations in His8, Asn59, and Gly295 improved the thermostability of MeBglD2, and the combination of these mutations resulted in the highest thermostability. Compared with wild-type MeBglD2, thermostable MeBglD2 mutants promoted plant biomass saccharification using Trichoderma reesei cellulase. In addition to thermostability, the thermostable mutants exhibited higher tolerance to ethanol, dimethyl sulfoxide, and copper ions, indicating that the MeBglD2 mutants generated in this study were improved in their tolerance to not only high temperature but also to organic solvents and metal ions.

     

Residues Asn214, Gln211, Glu219 and Gln221 contained in the subfamily 3 catalytic signature of the isocitrate lyase from <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis> are involved in its catalytic and thermal properties

Isocitrate lyase, encoded by the aceA gene, plays an important role in the ability of Pseudomonas aeruginosa to grow on fatty acids, acetate, acyclic terpenes, and amino acids. Phylogenetic analysis indicated that the ICL superfamily is divided in two families: the ICL family, which includes five subfamilies, and the 2-methylisocitrate lyase (MICL) family. ICL from P. aeruginosa (ICL-Pa) was identified in a different ICL node (subfamily 3) than other Pseudomonas ICL enzymes (grouped in subfamily 1). Analysis also showed that psychrophilic bacteria are mainly grouped in ICL subfamily 3, whose ICL proteins contain the highly conserved catalytic pattern Q_IEN_QV_SDE_KQ_CGHQD. We performed site-directed mutagenesis, enzymatic activity, and structure modeling of conserved residues in mutated ICLs by using ICL-Pa as a model. Our results indicated that the N214 residue is essential for catalytic function, while mutating the Q211, E219, and Q221 residues impairs its catalytic and thermostability properties. Our findings suggest that conserved residues in the subfamily 3 signature of ICL-Pa play important roles in catalysis and thermostability and are likely associated with the catalytic loop structural conformation.     

Phosphorylation of myosin A regulates gliding motility and is essential for Plasmodium transmission

Malaria‐causing parasites rely on an actin–myosin‐based motor for the invasion of different host cells and tissue traversal in mosquitoes and vertebrates. The unusual myosin A of Plasmodium spp. has a unique N‐terminal extension, which is important for red blood cell invasion by P. falciparum merozoites in vitro and harbors a phosphorylation site at serine 19. Here, using the rodent‐infecting P. berghei we show that phosphorylation of serine 19 increases ookinete but not sporozoite motility and is essential for efficient transmission of Plasmodium by mosquitoes as S19A mutants show defects in mosquito salivary gland entry. S19A along with E6R mutations slow ookinetes and salivary gland sporozoites in both 2D and 3D environments. In contrast to data from purified proteins, both E6R and S19D mutations lower force generation by sporozoites. Our data show that the phosphorylation cycle of S19 influences parasite migration and force generation and is critical for optimal migration of parasites during transmission from and to the mosquito.     

Mutations in PIEZO2 Cause Gordon Syndrome,Marden-Walker Syndrome,and Distal Arthrogryposis Type 5

Gordon syndrome (GS), or distal arthrogryposis type 3, is a rare, autosomal-dominant disorder characterized by cleft palate and congenital contractures of the hands and feet. Exome sequencing of five GS-affected families identified mutations in piezo-type mechanosensitive ion channel component 2 (PIEZO2) in each family. Sanger sequencing revealed PIEZO2 mutations in five of seven additional families studied (for a total of 10/12 [83%] individuals), and nine families had an identical c.8057G>A (p.Arg2686His) mutation. The phenotype of GS overlaps with distal arthrogryposis type 5 (DA5) and Marden-Walker syndrome (MWS). Using molecular inversion probes for targeted sequencing to screen PIEZO2, we found mutations in 24/29 (82%) DA5-affected families and one of two MWS-affected families. The presence of cleft palate was significantly associated with c.8057G>A (Fisher’s exact test, adjusted p value < 0.0001). Collectively, although GS, DA5, and MWS have traditionally been considered separate disorders, our findings indicate that they are etiologically related and perhaps represent variable expressivity of the same condition.     

Simultaneous Enhancement of Thermostability and Catalytic Activity of Phospholipase A1 by Evolutionary Molecular Engineering

The thermal stability and catalytic activity of phospholipase A₁ from Serratia sp. strain MK1 were improved by evolutionary molecular engineering. Two thermostable mutants were isolated after sequential rounds of error-prone PCR performed to introduce random mutations and filter-based screening of the resultant mutant library; we determined that these mutants had six (mutant TA3) and seven (mutant TA13) amino acid substitutions. Different types of substitutions were found in the two mutants, and these substitutions resulted in an increase in nonploar residues (mutant TA3) or in differences between side chains for polar or charged residues (mutant TA13). The wild-type and mutant enzymes were purified, and the effect of temperature on the stability and catalytic activity of the enzymes was investigated. The melting temperatures of the TA3 and TA13 enzymes were increased by 7 and 11°C, respectively, compared with the melting temperature of the wild-type enzyme. Thus, we found that evolutionary molecular engineering was an effective and efficient approach for increasing thermostability without compromising enzyme activity.     

High-throughput sequencing reveals suppressors of Vibrio cholerae rpoE mutations: one fewer porin is enough

Analyses of suppressor mutations have been extremely valuable in understanding gene function. However, techniques for mapping suppressor mutations are not available for most bacterial species. Here, we used high-throughput sequencing technology to identify spontaneously arising suppressor mutations that enabled disruption of rpoE (which encodes σ^E) in Vibrio cholerae, the agent of cholera. The alternative sigma factor σ^E, which is activated by envelope stress, promotes expression of factors that help preserve and/or restore cell envelope integrity. In Escherichia coli, rpoE is an essential gene that can only be disrupted in the presence of additional suppressor mutations. Among a panel of independent V. cholerae rpoE mutants, more than 75% contain suppressor mutations that reduce production of OmpU, V. cholerae’s principal outer membrane porin. OmpU appears to be a key determinant of V. cholerae’s requirement for and production of σ^E. Such dependence upon a single factor contrasts markedly with regulation of σ^E in E. coli, in which numerous factors contribute to its activation and none is dominant. We also identified a suppressor mutation that differs from all previously described suppressors in that it elevates, rather than reduces, σ^E’s activity. Finally, analyses of a panel of rpoE mutants shed light on the mechanisms by which suppressor mutations may arise in V. cholerae.     

Stress Tolerance in Doughs of Saccharomyces cerevisiae Trehalase Mutants Derived from Commercial Baker’s Yeast

Accumulation of trehalose is widely believed to be a critical determinant in improving the stress tolerance of the yeast Saccharomyces cerevisiae, which is commonly used in commercial bread dough. To retain the accumulation of trehalose in yeast cells, we constructed, for the first time, diploid homozygous neutral trehalase mutants (Δnth1), acid trehalase mutants (Δath1), and double mutants (Δnth1 ath1) by using commercial baker’s yeast strains as the parent strains and the gene disruption method. During fermentation in a liquid fermentation medium, degradation of intracellular trehalose was inhibited with all of the trehalase mutants. The gassing power of frozen doughs made with these mutants was greater than the gassing power of doughs made with the parent strains. The Δnth1 and Δath1 strains also exhibited higher levels of tolerance of dry conditions than the parent strains exhibited; however, the Δnth1 ath1 strain exhibited lower tolerance of dry conditions than the parent strain exhibited. The improved freeze tolerance exhibited by all of the trehalase mutants may make these strains useful in frozen dough.     

Improving the thermostability of methyl parathion hydrolase from Ochrobactrum sp. M231 using a computationally aided method

Good protein thermostability is very important for the protein application. In this report, we propose a strategy which contained a prediction method to select residues related to protein thermal stability, but not related to protein function, and an experiment method to screen the mutants with enhanced thermostability. The prediction strategy was based on the calculated site evolutionary entropy and unfolding free energy difference between the mutant and wild-type (WT) methyl parathion hydrolase enzyme from Ochrobactrum sp. M231 [Ochr-methyl parathion hydrolase (MPH)]. As a result, seven amino acid sites within Ochr-MPH were selected and used to construct seven saturation mutagenesis libraries. The results of screening these libraries indicated that six sites could result in mutated enzymes exhibiting better thermal stability than the WT enzyme. A stepwise evolutionary approach was designed to combine these selected mutants and a mutant with four point mutations (S274Q/T183E/K197L/S192M) was selected. The T _m and T ₅₀ of the mutant enzyme were 11.7 and 10.2 °C higher, respectively, than that of the WT enzyme. The success of this design methodology for Ochr-MPH suggests that it was an efficient strategy for enhancing protein thermostability and suitable for protein engineering.