A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35?V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35?V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.     

Prediction of thermostability from amino acid attributes by combination of clustering with attribute weighting: a new vista in engineering enzymes

The engineering of thermostable enzymes is receiving increased attention. The paper, detergent, and biofuel industries, in particular, seek to use environmentally friendly enzymes instead of toxic chlorine chemicals. Enzymes typically function at temperatures below 60°C and denature if exposed to higher temperatures. In contrast, a small portion of enzymes can withstand higher temperatures as a result of various structural adaptations. Understanding the protein attributes that are involved in this adaptation is the first step toward engineering thermostable enzymes. We employed various supervised and unsupervised machine learning algorithms as well as attribute weighting approaches to find amino acid composition attributes that contribute to enzyme thermostability. Specifically, we compared two groups of enzymes: mesostable and thermostable enzymes. Furthermore, a combination of attribute weighting with supervised and unsupervised clustering algorithms was used for prediction and modelling of protein thermostability from amino acid composition properties. Mining a large number of protein sequences (2090) through a variety of machine learning algorithms, which were based on the analysis of more than 800 amino acid attributes, increased the accuracy of this study. Moreover, these models were successful in predicting thermostability from the primary structure of proteins. The results showed that expectation maximization clustering in combination with uncertainly and correlation attribute weighting algorithms can effectively (100%) classify thermostable and mesostable proteins. Seventy per cent of the weighting methods selected Gln content and frequency of hydrophilic residues as the most important protein attributes. On the dipeptide level, the frequency of Asn-Glu was the key factor in distinguishing mesostable from thermostable enzymes. This study demonstrates the feasibility of predicting thermostability irrespective of sequence similarity and will serve as a basis for engineering thermostable enzymes in the laboratory.     

A rigid network of long-range contacts increases thermostability in a mutant endoglucanase

Thermodynamic stability of a protein at elevated temperatures is a key factor for thermostable enzymes to catalyze their specific reactions. Yet our understanding of biological determinants of thermostability is far from complete. Many different atomistic factors have been suggested as possible means for such proteins to preserve their activity at high temperatures. Among these factors are specific local interatomic interactions or enrichment of specific amino acid types. The case of glycosyl hydrolase family endoglucanase of Trichoderma reesei defies current hypotheses for thermostability because a single mutation far from the active site (A35?V) converts this mesostable protein into a thermostable protein without significant change in the protein structure. This substantial change in enzymatic activity cannot be explained on the basis of local intramolecular interactions alone. Here we present a more global view of the induced thermostability and show that the A35?V mutation affects the underlying structural rigidity of the whole protein via a number of long-range, non-local interactions. Our analysis of this structure reveals a precisely tuned, rigid network of atomic interactions. This cooperative, allosteric effect promotes the transformation of this mesostable protein into a thermostable one.     

Dynamics of endoglucanase catalytic domains: implications towards thermostability

Thermostable endoglucanases play a crucial role in the production of biofuels to breakdown plant cellulose. Analyzing their structure-dynamics relationship can inform about the origins of their thermostability. Although tertiary structures of many endoglucanase proteins are available, the relationship between thermostability, structure, and dynamics is not explored fully. We have generated elastic network models for thermostable and mesostable endoglucanases with the (αβ)? fold in substrate bound and unbound states. The comparative analyses shed light on the relation between protein dynamics, thermostability, and substrate binding. We observed specific differences in the dynamic behavior of catalytic residues in slow modes: while both the nucleophile and the acid/base donor residues show positively correlated motions in the thermophile, their dynamics is uncoupled in the mesophile. Our proof-of-concept comparison study suggests that global dynamics can be harnessed to further our understanding of thermostability.     

Discrimination of thermostable and thermophilic lipases using support vector machines

Discriminating thermophilic lipases from their similar thermostable counterparts is a challenging task and it would help to design stable proteins. In this study, the distributions of N (N=2, 3) neighboring amino acids and the non-adjacent di-residue coupling patterns in the sequences of 65 thermostable and 77 thermophilic lipases had been systematically analyzed. It was found that the hydrophobic residues Leu, Pro, Met, Phe, Trp, as well as the polar residue Tyr had higher occurrence in thermophilic lipases than thermostable ones. The occurrence frequencies of KC EE KE RE, VE, YI, EK, VK, EV, YV, EY, KY, VY and YY in thermophilic proteins were significantly higher, while the occurrence frequencies of QC, QH, QN, HQ, MQ, NQ, QQ, TQ, QS and QT were significantly lower. CXP or CPX showed significantly positive to lipase thermostability, while XXQ or QXX showed significantly negative to lipase thermostability. Non-adjacent di-residue coupling patterns of PR14, RY32, YR47, LE53, LE64, PP64, RP70 and PP101 were significantly different in thermophilic lipases and their thermostable counterparts. The composition of dipeptide, tripeptide and non-adjacent di-residue patterns contained more information than amino acid composition. A statistical method based on support vector machines (SVMs) was developed for discriminating thermophilic and thermostable lipases. The accuracy of this method for the training dataset was 97.17?. Furthermore, the highest accuracy of the method for testing datasets was 98.41?. The influence of some specific patterns on lipase thermostability was also discussed.     

Crystal structure of a thermostable lipase from Bacillus stearothermophilus P1

We describe the first lipase structure from a thermophilic organism. It shares less than 20% amino acid sequence identity with other lipases for which there are crystal structures, and shows significant insertions compared with the typical alpha/beta hydrolase canonical fold. The structure contains a zinc-binding site which is unique among all lipases with known structures, and which may play a role in enhancing thermal stability. Zinc binding is mediated by two histidine and two aspartic acid residues. These residues are present in comparable positions in the sequences of certain lipases for which there is as yet no crystal structural information, such as those from Staphylococcal species and Arabidopsis thaliana. The structure of Bacillus stearothermophilus P1 lipase provides a template for other thermostable lipases, and offers insight into mechanisms used to enhance thermal stability which may be of commercial value in engineering lipases for industrial uses.     

Directed evolution of Bacillus licheniformis lipase for improvement of thermostability

Lipases catalyze the hydrolysis of carboxylic acid esters and owing to their vast substrate specificity, they have many industrial applications. Due to the demand of thermostable lipases in industrial applications, we have enhanced the thermostability of lipase from Bacillus licheniformis RSP-09. The thermostable mutant lipases of Bacillus licheniformis RSP-09 were isolated following two rounds of directed evolution using error-prone PCR. The best mutant lipases obtained after first and second round of error-prone PCR were purified and characterized. The mutant lipases showed increased thermostability and retained catalytic function. The best mutant lipase (eP-231-51) showed 13.5-fold increase in percentage thermal stability (% remaining activity after incubation of purified enzyme at 60 °C for 1 h) than wild-type lipase. Also, this mutant lipase (ep-231-51) showed 30% improved catalytic efficiency compared with the wild-type which is due to significant decrease in K_m and marginal increase in k_cat. In addition, the thermostable mutant lipases have shown resistance to hydrophobic organic solvents. The role of mutations in the best mutant lipases of second round i.e. eP-231-51 (Asp72Gly, Asp61Gly, Tyr129His, and Thr101Pro) and eP-231-137 (Leu49Arg, Thr101Pro, Asp72Gly), that led to thermostability have been postulated after the comparison of molecular models of wild-type and mutated enzymes.     

Tiny TIM: a small, tetrameric, hyperthermostable triosephosphate isomerase

Comparative structural studies on proteins derived from organisms with growth optima ranging from 15 to 100 degrees C are beginning to shed light on the mechanisms of protein thermoadaptation. One means of sustaining hyperthermostability is for proteins to exist in higher oligomeric forms than their mesophilic homologues. Triosephosphate isomerase (TIM) is one of the most studied enzymes, whose fold represents one of nature's most common protein architectures. Most TIMs are dimers of approximately 250 amino acid residues per monomer. Here, we report the 2.7 A resolution crystal structure of the extremely thermostable TIM from Pyrococcus woesei, a hyperthermophilic archaeon growing optimally at 100 degrees C, representing the first archaeal TIM structure. P. woesei TIM exists as a tetramer comprising monomers of only 228 amino acid residues. Structural comparisons with other less thermostable TIMs show that although the central beta-barrel is largely conserved, severe pruning of several helices and truncation of some loops give rise to a much more compact monomer in the small hyperthermophilic TIM. The classical TIM dimer formation is conserved in P. woesei TIM. The extreme thermostability of PwTIM appears to be achieved by the creation of a compact tetramer where two classical TIM dimers interact via an extensive hydrophobic interface. The tetramer is formed through largely hydrophobic interactions between some of the pruned helical regions. The equivalent helical regions in less thermostable dimeric TIMs represent regions of high average temperature factor. The PwTIM seems to have removed these regions of potential instability in the formation of the tetramer. This study of PwTIM provides further support for the role of higher oligomerisation states in extreme thermal stabilisation.     

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.

     

Improvement of 2'-hydroxybiphenyl-2-sulfinate desulfinase, an enzyme involved in the dibenzothiophene desulfurization pathway, from Rhodococcus erythropolis KA2-5-1 by site-directed mutagenesis

In the microbial dibenzothiophene desulfurization pathway, 2'-hydroxybiphenyl-2-sulfinate is converted to 2-hydroxybiphenyl and sulfinate by desulfinase (DszB) at the last step, and this reaction is rate-limiting for the whole pathway. The catalytic activity and thermostability of DszB were enhanced by the two amino acid substitutions. Based on information on the 3-D structure of DszB and a comparison of amino acid sequences between DszB and reported thermophilic and thermostable homologs (TdsB and BdsB), two amino acid residues, Tyr63 and Gln65, were selected as targets to mutate and improve DszB. These two residues were replaced by several amino acids, and the promising mutant enzymes were purified and their properties were examined. Among the wild-type and mutant enzymes, Y63F had higher catalytic activity but similar thermostability, and Q65H showed higher thermostability but less catalytic activity and affinity for the substrate. To compensate for these drawbacks, the double mutant enzyme Y63F-Q65H was purified and its properties were investigated. This mutant enzyme showed higher thermostability without loss of catalytic activity or affinity for the substrate. These superior properties of the mutant enzyme have also been confirmed with resting cells harboring the mutant gene.     

A mutant phospholipase D with enhanced thermostability from Streptomyces sp.

To investigate the contribution of amino acid residues to the thermostability of phospholipase D (PLD), a chimeric form of two Streptomyces PLDs (thermolabile K1PLD and thermostable TH-2PLD) was constructed. K/T/KPLD, in which residues 329–441 of K1PLD were recombined with the homologous region of TH-2PLD, showed a thermostability midway between those of K1PLD and TH-2PLD. By comparing the primary structures of Streptomyces PLDs, the seven candidates of thermostability-related amino acid residues of K1PLD were identified. The K1E346DPLD mutant, in which Glu346 of K1PLD was substituted with Asp by site-directed mutagenesis, exhibited enhanced thermostability, which was almost the same as that of TH-2PLD.     

Relative tolerance of mesostable and thermostable protein homologs to extensive mutation

Evolvability, designability, and plasticity of a protein are properties that are important to protein engineers, but difficult to quantify. Here, we directly compare homologous AroQ chorismate mutases from the thermophile Methanococcus jannaschii and the mesophile Escherichia coli with respect to their capacity to accommodate extensive mutation. The N-terminal helix comprising about 40% of these proteins was randomized at the genetic level using a binary pattern of hydrophobic and hydrophilic residues based on the respective wild-type sequences. Catalytically active library members were identified by a survival-selection assay in a chorismate mutase-deficient E. coli strain. Functional variants were found approximately approximately 10-times more frequently with the thermostable protein compared to its mesostable counterpart. Moreover, detailed sequence analysis revealed that functional M. jannaschii enzyme variants contained a smaller number of conserved residues and tolerated greater variability at individual sequence positions. Our results thus highlight the greater robustness of the thermostable protein with respect to amino acid substitution, while identifying specific sites important for constructing active enzymes. Overall, they support the notion that redesign projects will benefit from using a thermostable starting structure, even at very high mutational loads.     

A mutant phospholipase D with enhanced thermostability from Streptomyces sp

To investigate the contribution of amino acid residues to the thermostability of phospholipase D (PLD), a chimeric form of two Streptomyces PLDs (thermolabile K1PLD and thermostable TH-2PLD) was constructed. K/T/KPLD, in which residues 329-441 of K1PLD were recombined with the homologous region of TH-2PLD, showed a thermostability midway between those of K1PLD and TH-2PLD. By comparing the primary structures of Streptomyces PLDs, the seven candidates of thermostability-related amino acid residues of K1PLD were identified. The K1E346DPLD mutant, in which Glu346 of K1PLD was substituted with Asp by site-directed mutagenesis, exhibited enhanced thermostability, which was almost the same as that of TH-2PLD.     

Exploring the sequence patterns in the alpha-helices of proteins

This paper reports an extensive sequence analysis of the alpha-helices of proteins. alpha-Helices were extracted from the Protein Data Bank (PDB) and were divided into groups according to their sizes. It was found that some amino acids had differential propensity values for adopting helical conformation in short, medium and long alpha-helices. Pro and Trp had a significantly higher propensity for helical conformation in short helices than in medium and long helices. Trp was the strongest helix conformer in short helices. Sequence patterns favoring helical conformation were derived from a neighbor-dependent sequence analysis of proteins, which calculated the effect of neighboring amino acid type on the propensity of residues for adopting a particular secondary structure in proteins. This method produced an enhanced statistical significance scale that allowed us to explore the positional preference of amino acids for alpha-helical conformations. It was shown that the amino acid pair preference for alpha-helix had a unique pattern and this pattern was not always predictable by assuming proportional contributions from the individual propensity values of the amino acids. Our analysis also yielded a series of amino acid dyads that showed preference for alpha-helix conformation. The data presented in this study, along with our previous study on loop sequences of proteins, should prove useful for developing potential 'codes' for recognizing sequence patterns that are favorable for specific secondary structural elements in proteins.     

Thermophilic nucleoside phosphorylases: Their properties,characteristics and applications

Nucleoside phosphorylases catalyze the reversible phosphorolysis of pyrimidine and purine nucleosides in the presence of phosphate. They are valuable catalysts in the synthesis of nucleosides and their analogues, which are often used as pharmaceuticals or their precursors. Thermostable nucleoside phosphorylases are promising biocatalysts, as they withstand harsh reaction conditions such as high pH or the addition of organic solvents. In this review, the characteristics and properties of thermostable nucleoside phosphorylases are described. Differences in amino acid content and protein structure were compared to their mesophilic homologues to identify features involved in thermostability. Substrate spectra of thermostable nucleoside phosphorylases were analyzed, and it is shown that thermostable nucleoside phosphorylases have a wider substrate spectrum than their mesophilic counterparts. Thus, thermostable nucleoside phosphorylases are interesting biocatalysts for industrial applications.     

The effect of proline insertions on the thermostability of a barley alpha-glucosidase

The thermal stability of alpha-glucosidase is important because the conversion of starch to fermentable sugars during industrial production of ethanol (e.g. brewing, fuel ethanol production) typically takes place at temperatures of 65-73 degrees C. In this study we investigate the thermostability of alpha-glucosidases from four plant species, compare their deduced amino acid sequences, and test the effect of substituting a proline for the residue present in the wild-type enzyme on the thermostability of alpha-glucosidase. The alpha-glucosidase from barley (Hordeum vulgare) was significantly less thermostable than the other three alpha-glucosidases. A comparison of the published deduced amino acid sequences of these four alpha-glucosidases revealed conserved proline residues in the three most thermostable alpha-glucosidases that were not found in the barley enzyme. Site-directed mutagenesis was done on recombinant barley alpha-glucosidase to create proteins with prolines at these conserved positions. The thermostability (T(50)) of one of these mutant enzymes, T340P, was 10 degrees C higher than the non-mutated enzyme.     

Crystal structure of D-Hydantoinase from Burkholderia pickettii at a resolution of 2.7 Angstroms: insights into the molecular basis of enzyme thermostability

D-Hydantoinase (D-HYD) is an industrial enzyme that is widely used in the production of D-amino acids which are precursors for semisynthesis of antibiotics, peptides, and pesticides. This report describes the crystal structure of D-hydantoinase from Burkholderia pickettii (HYD(Bp)) at a 2.7-A resolution. The structure of HYD(Bp) consists of a core (alpha/beta)(8) triose phosphate isomerase barrel fold and a beta-sheet domain, and the catalytic active site consists of two metal ions and six highly conserved amino acid residues. Although HYD(Bp) shares only moderate sequence similarity with D-HYDs from Thermus sp. (HYD(Tsp)) and Bacillus stearothermophilus (HYD(Bst)), whose structures have recently been solved, the overall structure and the structure of the catalytic active site are strikingly similar. Nevertheless, the amino acids that compose the substrate-binding site are less conserved and have different properties, which might dictate the substrate specificity. Structural comparison has revealed insights into the molecular basis of the differential thermostability of D-HYDs. The more thermostable HYD(Tsp) contains more aromatic residues in the interior of the structure than HYD(Bp) and HYD(Bst). Changes of large aromatic residues in HYD(Tsp) to smaller residues in HYD(Bp) or HYD(Bst) decrease the hydrophobicity and create cavities inside the structure. HYD(Tsp) has more salt bridges and hydrogen-bonding interactions and less oxidation susceptible Met and Cys residues on the protein surface than HYD(Bp) and HYD(Bst). Besides, HYD(Tsp) also contains more rigid Pro residues. These factors are likely to make major contributions to the varying thermostability of these enzymes. This information could be exploited in helping to engineer more thermostable mesophilic enzymes.     

Analysis of Comparative Sequence and Genomic Data to Verify Phylogenetic Relationship and Explore a New Subfamily of Bacterial Lipases

Thermostable and organic solvent-tolerant enzymes have significant potential in a wide range of synthetic reactions in industry due to their inherent stability at high temperatures and their ability to endure harsh organic solvents. In this study, a novel gene encoding a true lipase was isolated by construction of a genomic DNA library of thermophilic Aneurinibacillus thermoaerophilus strain HZ into Escherichia coli plasmid vector. Sequence analysis revealed that HZ lipase had 62% identity to putative lipase from Bacillus pseudomycoides. The closely characterized lipases to the HZ lipase gene are from thermostable Bacillus and Geobacillus lipases belonging to the subfamily I.5 with ≤ 57% identity. The amino acid sequence analysis of HZ lipase determined a conserved pentapeptide containing the active serine, GHSMG and a Ca²⁺-binding motif, GCYGSD in the enzyme. Protein structure modeling showed that HZ lipase consisted of an α/β hydrolase fold and a lid domain. Protein sequence alignment, conserved regions analysis, clustal distance matrix and amino acid composition illustrated differences between HZ lipase and other thermostable lipases. Phylogenetic analysis revealed that this lipase represented a new subfamily of family I of bacterial true lipases, classified as family I.9. The HZ lipase was expressed under promoter P_lac using IPTG and was characterized. The recombinant enzyme showed optimal activity at 65°C and retained ≥ 97% activity after incubation at 50°C for 1h. The HZ lipase was stable in various polar and non-polar organic solvents.     

Increasing in archaeal GlcNAc-1-P uridyltransferase activity by targeted mutagenesis while retaining its extreme thermostability

UDP-GlcNAc, an activated and essential form of GlcNAc which is an important component in the polysaccharide structure of most organisms, is synthesized from GlcNAc-1-P and UTP by GlcNAc-1-P UTase. We previously reported the identification of the extremely thermostable ST0452 protein, which has dual sugar-1-P NTase activities (Glc-1-P TTase and GlcNAc-1-P UTase activities) from an acidothermophilic archaeon, Sulfolobus tokodaii strain 7. Detailed analyses of the protein indicated that the activity is slightly lower than that of bacteria. For industrial applications, activity needs to be increased without decreasing thermostability. Therefore, to enhance this activity, we introduced mutations into the amino acid residues located within the predicted reaction centre by targeted mutagenesis. All 12 mutant ST0452 proteins showed no decrease in thermostability. Among them, six mutant proteins were found to have increased GlcNAc-1-P UTase activity under optimal reaction conditions with sufficient substrates or an appropriate metal ion. Our results indicate that targeted mutagenesis is a powerful technique for in vitro production of a thermostable enzyme with enhanced activity. The results of this study also indicate that the space for the metal ion is important for selecting the type of metal ion and also affects the rate of the reaction.     

Improving the thermostability of <Emphasis Type="Italic">N</Emphasis>-carbamyl-<Emphasis Type="SmallCaps">d</Emphasis>-amino acid amidohydrolase by error-prone PCR

To facilitate the easier production of d-amino acids using N-carbamyl-d-amino acid amidohydrolase (DCase) in an immobilized form, we improved the enzymatic thermostability of highly soluble DCase-M3 of Ralstonia pickettii using directed mutagenesis. Six novel mutation sites were identified in this study, apart from several thermostability-related amino acid sites reported previously. The most thermostable mutant, in which the 12th amino acid had been changed from glutamine to leucine, showed a 7 °C increase in thermostability. Comparative characterization of the parental and mutant DCases showed that although there was a slight reduction in the oxidative stability of the mutants, their kinetic properties and high solubility were not affected. The mutated enzymes are expected to be applied to the development of a fully enzymatic process for the industrial production of d-amino acids.