Effective factors in thermostability of thermophilic proteins

Bringing the heat: Thermostable analogs of Bst polymerase allow high-temperature LAMP

BACKGROUND

Loop-mediated isothermal amplification (LAMP) is a nucleic acid amplification method that gained prominence during the early months of the COVID-19 pandemic due to its simplicity, sensitivity and robustness. However, this technique is susceptible to non-specific amplifications, raising concerns about false-positive results and reduced diagnostic accuracy. A primary contributor to false-positive testing is primer dimerization, which can theoretically be mitigated by performing reactions at higher temperatures. Unfortunately, the strand-displacing DNA polymerases typically used in LAMP, such as Bst, exhibit reduced efficiency at elevated temperatures. To address this limitation, we hypothesised that naturally occurring thermophilic analogs of Bst may be capable of supporting LAMP at higher temperatures, thereby improving reaction specificity.

METHODS

Bioinformatics and recombinant enzyme production allowed the identification and synthesis of several Bst analogs. These were tested in real-time LAMP assays to detect diverse targets, in a wide range of reaction temperatures (63°C-75°C) and in the presence of typical qPCR inhibitors.

RESULTS

Three polymerases-Bst_7, Bst_8 and Bst_15-demonstrated exceptional activity and robust stability at higher temperature conditions (up to 72.5°C), while displaying considerable resistance to common qPCR inhibitors.

CONCLUSIONS

The identified thermophilic Bst analogs represent a potential solution for the mitigation of non-specific amplification in LAMP, further boosting the application of this technique in molecular diagnostic settings.

Discovery and In Silico Characterization of Anatolian Water Buffalo Rumen-Derived Bacterial Thermostable Xylanases: A Sequence-Based Metagenomic Approach

This study involved shotgun sequencing of rumen metagenomes from three Anatolian water buffalos, an exploration of the relationship between microbial flora and xylanases, and in silico analyses of thermostable xylanases, focusing on their sequence, structure, and dynamic properties. For this purpose, the rumen metagenome of three Anatolian water buffalos was sequenced and bioinformatically analyzed to determine microbial diversity and full-length xylanases. Analyses of BLAST, biophysicochemical characteristics, phylogenetic tree, and multiple sequence alignment were performed with Blastp, ProtParam, MEGA11 software, and Clustal Omega, respectively. Three-dimensional homology models of three xylanases (AWBRMetXyn5, AWBRMetXyn10, and AWBRMetXyn19) were constructed by SWISS-MODEL and validated by ProSA, ProCheck, and Verify3D. Also, their 3D models were structurally analyzed by PyMOL, BANΔIT, thermostability predictor, What If, and Protein Interaction Calculator (PIC) software. Protein-ligand interactions were examined by docking and MD simulation. Shotgun sequence and Blastp analyses showed that Clostridium (Clostridiales bacterial order), Ruminococcus (Oscillospiraceae bacterial family), Prevotella (Bacteroidales bacterial order), and Butyrivibrio (Lachnospiraceae bacterial family) were found as dominant potential xylanase-producer genera in three rumen samples. Furthermore, the biophysicochemical analysis indicated that three xylanases exhibited an aliphatic index above 80, an instability index below 40, and melting temperatures (T m) surpassing 65 °C. Phylogenetic analysis placed three xylanases within the GH10 family, clustering them with thermophilic xylanases, while homology modeling identified the optimal template as a xylanase from a thermophilic bacterium. The structural analysis indicated that three xylanases possessed the number of salt bridges, hydrophobic interactions, and T m score higher than 50, 165, and 70 °C, respectively; however, the reference thermophilic XynAS9 had 43, 145, and 54.41 °C, respectively. BANΔIT analysis revealed that three xylanases exhibited lower B'-factor values in the β3-α1 loop/short-helix at the N-terminal site compared to the reference thermophilic XynAS9. In contrast, six residues (G79, M123, D150, T199, A329, and G377) possessed higher B'-factor values in AWBRMetXyn5 and their aligned positions in AWBRMetXyn10 and AWBRMetXyn19, relative to XynAS9 including Gln, Glu, Ile, Lys, Ser, and Val at these positions, respectively. MD simulation results showed that the β9-η5 loop including catalytic nucleophile glutamic acid in the RMSF plot of three xylanases had a higher fluctuation than the aligned region in XynAS9. The distance analysis from the MD simulation showed that the nucleophile residue in AWBRMetXyn5 and AWBRMetXyn10 remained closer to the ligand throughout the simulation compared with XynAS9 and AWBRMetXyn19. The most notable difference between AWBRMetXyn5 and AWBRMetXyn10 was the increased amino acid fluctuations in two specific regions, the η3 short-helix and the η3-α3 loop, despite a minimal sequence difference of only 1.24%, which included three key amino acid variations (N345, N396, and T397 in AWBRMetXyn5; D345, K396, and A397 in AWBRMetXyn10). Thus, this study provided computational insights into xylanase function and thermostability, which could inform future protein engineering efforts. Additionally, three xylanases, especially AWBRMetXyn5, are promising candidates for various high-temperature industrial applications. In a forthcoming study, three xylanases will be experimentally characterized and considered for potential industrial applications. In addition, the amino acid substitutions (G79Q, M123E, D150I, T199K, A329S, and G377V) and the residues in the β3-α1 loop will be targeted for thermostability improvement of AWBRMetXyn5. The amino acids (N345, N396, and T397) and the residues on the β9-η5 loop, η3 short-helix, and η3-α3 loop will also be focused on development of the catalytic efficiency.

ProCeSa: Contrast-Enhanced Structure-Aware Network for Thermostability Prediction with Protein Language Models

Proteins play a fundamental role in biology, and their thermostability is essential for their proper functionality. The precise measurement of thermostability is crucial, traditionally relying on resource-intensive experiments. Recent advances in deep learning, particularly in protein language models (PLMs), have significantly accelerated the progress in protein thermostability prediction. These models utilize various biological characteristics or deep representations generated by PLMs to represent the protein sequences. However, effectively incorporating structural information, based on the PLM embeddings, while not considering atomic protein structures, remains an open and formidable challenge. Here, we propose a novel Protein Contrast-enhanced Structure-Aware (ProCeSa) model that seamlessly integrates both sequence and structural information extracted from PLMs to enhance thermostability prediction. Our model employs a contrastive learning scheme guided by the categories of amino acid residues, allowing it to discern intricate patterns within protein sequences. Rigorous experiments conducted on publicly available data sets establish the superiority of our method over state-of-the-art approaches, excelling in both classification and regression tasks. Our results demonstrate that ProCeSa addresses the complex challenge of predicting protein thermostability by utilizing PLM-derived sequence embeddings, without requiring access to atomic structural data.

Enhanced thermostability of nattokinase by rational design of disulfide bond

Nattokinase, the thrombolytically active substance in the health food natto, nevertheless, its lower thermostability restricts its use in food and pharmaceutical applications. In this study, two heat-resistant variants of nattokinase, designated 50-109 (M1) and 15-271 (M2), were successfully obtained by introducing a disulfide bonding strategy. Their half-lives at 55℃ were found to be 2.50-fold and 5.17-fold higher, respectively, than that of the wild type. Furthermore, the specific enzyme activities of the variants, M1 and M2, were also increased by 2.37 and 1.66-fold, respectively. Meanwhile, the combination of two mutants increased the thermostability of nattokinase by 8.0-fold. Bioinformatics analyses indicated that the enhanced thermostability of the M1 and M2 variants was due to the increased rigidity and structural contraction of the overall structure. Finally, the fermentation process of mutant M1 was optimized to increase the expression of nattokinase. Study provides substantial molecular and theoretical support for the industrial production and application of nattokinase.

Dual-domain superoxide dismutase: In silico prediction directed combinatorial mutation for enhanced robustness and catalytic efficiency

The robustness and catalytic activity of superoxide dismutase (SOD) are still the main factors limiting their application in industrial fields. This study aims to further improve the properties of a natural thermophilic iron/manganese dual-domain SOD (Fe/Mn-SODA fused with N-terminal polypeptide) from Geobacillus thermodenitrificans NG80-2 (GtSOD) by modifying its each domain using in-depth in silico prediction analysis as well as protein engineering. First, computational analysis of the N-terminal domain and GtSODA domain was respectively performed by using homologous sequence alignment and virtual mutagenesis. Seven proposed mutation sites favoring increased robustness were screened out for single-point mutants (SPMs) construction. Enzymatic characterization of these SPMs identified the most favorable mutation sites E107 and S265 located in two different domains. Subsequently, the dual-domain site combinatorial mutant (DDSCM) E107L/S265K showed significant superposition effects and additional improvement in catalytic efficiency, with a Kcat/Km value of 145.45 %, 33.66 %, and 60.33 % higher than the wild type (WT), the SPMs E107L and S265K, respectively. Molecular dynamics simulations, structural and surface charge analysis revealed the possible mechanism by which combinatorial mutations improve the robustness and catalytic activity of GtSOD. Furthermore, DDSCM showed more significant resistance to ultraviolet B and various stress than WT, indicating its highly competitive industrial application prospects.

Construction and characterization of a novel fusion alginate lyase with endolytic and exolytic cleavage activity for industrial preparation of alginate oligosaccharides

Alginate lyases with high activity and good thermostability are lacking for the preparation of alginate oligosaccharides (AOS) with various biological activities. We constructed a fusion alginate lyase with both endo-and exo-activities. AlyRm6A-Zu7 was successfully constructed by connecting the highly thermostable AlyRm6A to a new exotype lyase, AlyZu7. The fusion enzyme exhibited high catalytic activity and thermostability. It transformed sodium alginate into oligosaccharides with degrees of polymerization (DP) of 2-4 while producing 4-deoxy-L-erythro-5-hexoseulose uronic acid (DEH). The maximum reducing sugar, AOS, and DP1 + DEH yields were 75 %, 45 %, and 40 %, respectively. Molecular docking confirmed the formation of a stable complex between the substrate and AlyRm6A-Zu7. Protein interactions increased the thermostability of AlyZu7. This work provides new insights into the industrial formation of AOS and monosaccharide DEH using thermally stable fusion enzymes, which has a positive effect in the fields of functional oligosaccharide production and biofuel formation.

TemBERTure: advancing protein thermostability prediction with deep learning and attention mechanisms

Motivation

Understanding protein thermostability is essential for numerous biotechnological applications, but traditional experimental methods are time-consuming, expensive, and error-prone. Recently, deep learning (DL) techniques from natural language processing (NLP) was extended to the field of biology, since the primary sequence of proteins can be viewed as a string of amino acids that follow a physicochemical grammar.

Results

In this study, we developed TemBERTure, a DL framework that predicts thermostability class and melting temperature from protein sequences. Our findings emphasize the importance of data diversity for training robust models, especially by including sequences from a wider range of organisms. Additionally, we suggest using attention scores from Deep Learning models to gain deeper insights into protein thermostability. Analyzing these scores in conjunction with the 3D protein structure can enhance understanding of the complex interactions among amino acid properties, their positioning, and the surrounding microenvironment. By addressing the limitations of current prediction methods and introducing new exploration avenues, this research paves the way for more accurate and informative protein thermostability predictions, ultimately accelerating advancements in protein engineering.

Availability and implementation

TemBERTure model and the data are available at: https://github.com/ibmm-unibe-ch/TemBERTure.

Thermophilic β-mannanases from bacteria: production, resources, structural features and bioengineering strategies

β-mannanases are pivotal enzymes that cleave the mannan backbone to release short chain mannooligosaccharides, which have tremendous biotechnological applications including food/feed, prebiotics and biofuel production. Due to the high temperature conditions in many industrial applications, thermophilic mannanases seem to have great potential to overcome the thermal impediments. Thus, structural analysis of thermostable β-mannanases is extremely important, as it could open up new avenues for genetic engineering, and protein engineering of these enzymes with enhanced properties and catalytic efficiencies. Under this scope, the present review provides a state-of-the-art discussion on the thermophilic β-mannanases from bacterial origin, their production, engineering and structural characterization. It covers broad insights into various molecular biology techniques such as gene mutagenesis, heterologous gene expression, and protein engineering, that are employed to improve the catalytic efficiency and thermostability of bacterial mannanases for potential industrial applications. Further, the bottlenecks associated with mannanase production and process optimization are also discussed. Finally, future research related to bioengineering of mannanases with novel protein expression systems for commercial applications are also elaborated.

Molecular dynamics simulations at high temperatures of the Aeropyrum pernix L7Ae thermostable protein: Insight into the unfolding pathway

Most life forms on earth live at temperatures below 50 °C. Within these organisms are proteins that form the three-dimensional structures essential to their biological activity and function. However, some thermophilic life forms can resist higher temperatures and have corresponding adaptations to preserve protein function at these high temperatures. Among the structural factors responsible for this resistance of thermophilic proteins to high temperatures is the presence of additional hydrogen bonds in the thermophilic proteins, which means that the structure of the protein is more resistant to unfolding. Similarly, thermostable proteins are rich in structure-stabilizing salt bridges and/or disulfide bridges. In this context, we perform multiple replica molecular dynamics simulations at different temperatures on the Aeropyrum pernix (L7Ae) protein (from the crenarchaeal species A. pernix), known for its high melting temperature, and this in the aim to elucidate the structural factors responsible for its high thermostability. The results reveal that between the most sensitive regions of the protein to the increase of temperature are the loops L1, and L5, which surround the hydrophobic core region of the protein, besides the loop L9, and the C-terminal α5 region. This latter is the longer alpha helix of the protein secondary structure motifs and it is the first to be denaturated at 450 K, while the rest of the protein secondary structure motifs at this temperature were intact. The mechanism of unfolding that follows this protein at 550 K is similar to other thermophile proteins found in literature, with the opening of the loops that surround the hydrophobic core of the protein. So, the latter is completely exposed to the solvent, and partially denatured. The total denaturation process of the protein takes an average time of 40 ns to be achieved. Our investigation also shows that all the calculated salt bridges, with distances less than or equal to 6 A°, are on the periphery part of the protein, exposed to the solvent. However, the hydrophobic core of the protein is not involved in the formation of salt bridges, but rather with formation of some important hydrogen bondings that still persist even at 450 K. So, optimizing hydrogen bonding, near or within the core region, at high temperatures is a strategy that follows this thermostable protein to protect its hydrophobic core from denaturation, and ensure the thermal stability of the protein.

DeepPPThermo: A Deep Learning Framework for Predicting Protein Thermostability Combining Protein-Level and Amino Acid-Level Features

Using wet experimental methods to discover new thermophilic proteins or improve protein thermostability is time-consuming and expensive. Machine learning methods have shown powerful performance in the study of protein thermostability in recent years. However, how to make full use of multiview sequence information to predict thermostability effectively is still a challenge. In this study, we proposed a deep learning-based classifier named DeepPPThermo that fuses features of classical sequence features and deep learning representation features for classifying thermophilic and mesophilic proteins. In this model, deep neural network (DNN) and bi-long short-term memory (Bi-LSTM) are used to mine hidden features. Furthermore, local attention and global attention mechanisms give different importance to multiview features. The fused features are fed to a fully connected network classifier to distinguish thermophilic and mesophilic proteins. Our model is comprehensively compared with advanced machine learning algorithms and deep learning algorithms, proving that our model performs better. We further compare the effects of removing different features on the classification results, demonstrating the importance of each feature and the robustness of the model. Our DeepPPThermo model can be further used to explore protein diversity, identify new thermophilic proteins, and guide directed mutations of mesophilic proteins.

Prediction of thermophilic protein using 2-D general series correlation pseudo amino acid features

The demand for thermophilic protein has been increasing in protein engineering recently. Many machine-learning methods for identifying thermophilic proteins have emerged during this period. However, most machine learning-based thermophilic protein identification studies have only focused on accuracy. The relationship between the features' meaning and the proteins' physicochemical properties has yet to be studied in depth. In this article, we focused on the relationship between the features and the thermal stability of thermophilic proteins. This method used 2-D general series correlation pseudo amino acid (SC-PseAAC-General) features and realized accuracy of 82.76% using the J48 classifier. In addition, this research found the presence of higher frequencies of glutamic acid in thermophilic proteins, which help thermophilic proteins maintain their thermal stability by forming hydrogen bonds and salt bridges that prevent denaturation at high temperatures.

Thermostable enzyme research advances: a bibliometric analysis

Thermostable enzymes are enzymes that can withstand elevated temperatures as high as 50 °C without altering their structure or distinctive features. The potential of thermostable enzymes to increase the conversion rate at high temperature has been identified as a key factor in enhancing the efficiency of industrial operations. Performing procedures at higher temperatures with thermostable enzymes minimises the risk of microbial contamination, which is one of the most significant benefits. In addition, it helps reduce substrate viscosity, improve transfer speeds, and increase solubility during reaction operations. Thermostable enzymes offer enormous industrial potential as biocatalysts, especially cellulase and xylanase, which have garnered considerable amount of interest for biodegradation and biofuel applications. As the usage of enzymes becomes more common, a range of performance-enhancing applications are being explored. This article offers a bibliometric evaluation of thermostable enzymes. Scopus databases were searched for scientific articles. The findings indicated that thermostable enzymes are widely employed in biodegradation as well as in biofuel and biomass production. Japan, the United States, China, and India, as along with the institutions affiliated with these nations, stand out as the academically most productive in the field of thermostable enzymes. This study's analysis exposed a vast number of published papers that demonstrate the industrial potential of thermostable enzymes. These results highlight the significance of thermostable enzyme research for a variety of applications.

Whole-genome sequencing analysis of two heat-evolved Escherichia coli strains

BACKGROUND

High temperatures cause a suite of problems for cells, including protein unfolding and aggregation; increased membrane fluidity; and changes in DNA supercoiling, RNA stability, transcription and translation. Consequently, enhanced thermotolerance can evolve through an unknown number of genetic mechanisms even in the simple model bacterium Escherichia coli. To date, each E. coli study exploring this question resulted in a different set of mutations. To understand the changes that can arise when an organism evolves to grow at higher temperatures, we sequenced and analyzed two previously described E. coli strains, BM28 and BM28 ΔlysU, that have been laboratory adapted to the highest E. coli growth temperature reported to date.

RESULTS

We found three large deletions in the BM28 and BM28 ΔlysU strains of 123, 15 and 8.5 kb in length and an expansion of IS10 elements. We found that BM28 and BM28 ΔlysU have considerably different genomes, suggesting that the BM28 culture that gave rise to BM28 and BM28 ΔlysU was a mixed population of genetically different cells. Consistent with published findings of high GroESL expression in BM28, we found that BM28 inexplicitly carries the groESL bearing plasmid pOF39 that was maintained simply by high-temperature selection pressure. We identified over 200 smaller insertions, deletions, single nucleotide polymorphisms and other mutations, including changes in master regulators such as the RNA polymerase and the transcriptional termination factor Rho. Importantly, this genome analysis demonstrates that the commonly cited findings that LysU plays a crucial role in thermotolerance and that GroESL hyper-expression is brought about by chromosomal mutations are based on a previous misinterpretation of the genotype of BM28.

CONCLUSIONS

This whole-genome sequencing study describes genetically distinct mechanisms of thermotolerance evolution from those found in other heat-evolved E. coli strains. Studying adaptive laboratory evolution to heat in simple model organisms is important in the context of climate change. It is important to better understand genetic mechanisms of enhancing thermotolerance in bacteria and other organisms, both in terms of optimizing laboratory evolution methods for various organisms and in terms of potential genetic engineering of organisms most at risk or most important to our societies and ecosystems.

Structural features, temperature adaptation and industrial applications of microbial lipases from psychrophilic, mesophilic and thermophilic origins

Microbial lipases are very prominent biocatalysts because of their ability to catalyze a wide variety of reactions in aqueous and non-aqueous media. Here microbial lipases from different origins (psychrophiles, mesophiles, and thermophiles) have been reviewed. This review emphasizes an update of structural diversity in temperature adaptation and industrial applications, of psychrophilic, mesophilic, and thermophilic lipases. The microbial origins of lipases are logically dynamic, proficient, and also have an extensive range of industrial uses with the manufacturing of altered molecules. It is therefore of interest to understand the molecular mechanisms of adaptation to temperature in occurring lipases. However, lipases from extremophiles (psychrophiles, and thermophiles) are widely used to design biotransformation reactions with higher yields, fewer byproducts, or useful side products and have been predicted to catalyze those reactions also, which otherwise are not possible with the mesophilic lipases. Lipases as a multipurpose biological catalyst have given a favorable vision in meeting the needs of several industries such as biodiesel, foods, and drinks, leather, textile, detergents, pharmaceuticals, and medicals.

Enhancement of the polyethylene terephthalate and mono-(2-hydroxyethyl) terephthalate degradation activity of Ideonella sakaiensis PETase by an electrostatic interaction-based strategy

The serious environmental pollution that came up with the continuously growing demand for polyethylene terephthalate (PET) has attracted global concern. The IsPETase which has shown the highest PET degradation activity under ambient temperature is a promising enzyme for PET biodegradation, while poor thermostability limited its practical application. Herein, an electrostatic interaction-based strategy was applied for rational design of IsPETase towards enhanced thermostability. The IsPETaseI139R variant displayed the highest Tm value of 56.4 °C and 3.6-times higher PET degradation activity. Molecular simulations demonstrated that the introduction of salt bridges stabilized the local structures, resulting in robust thermostability. Meanwhile, the IsPETaseS92K/D157E/R251A not only exhibited higher thermostability but also showed a 1.74-fold kcat increase towards mono-(2-hydroxyethyl) terephthalate, which ultimately achieved PET depolymerization to complete monomer TPA. Collectively, the electrostatic interaction-based strategy, together with the derived IsPETase variants, could help promote the bio-recycle of PET, reducing the severe global burden of PET waste.

SAPPHIRE: A stacking-based ensemble learning framework for accurate prediction of thermophilic proteins

Thermophilic proteins (TPPs) are important in the field of protein biochemistry and development of new enzymes. Thus, computational methods must be urgently developed to accurately and rapidly identify TPPs. To date, several computational methods have been developed for TPP identification; however, few limitations in terms of performance and utility remain. In this study, we present a novel computational method, SAPPHIRE, to achieve more accurate identification of TPPs using only sequence information without any need for structural information. We combined twelve different feature encodings representing different perspectives and six popular machine learning algorithms to train 72 baseline models and extract the key information of TPPs. Subsequently, the informative predicted probabilities from the baseline models were mined and selected using a genetic algorithm in conjunction with a self-assessment-report approach. Finally, the final meta-predictor, SAPPHIRE, was built and optimized by applying an optimal feature set. The performance of SAPPHIRE in the 10-fold cross-validation test showed that a superior predictive performance compared with several baseline models could be achieved. Moreover, SAPPHIRE yielded an accuracy of 0.942 and Matthew's coefficient correlation of 0.884, which were 7.68 and 5.12% higher than those of the current existing methods, respectively, as indicated by the independent test. The proposed computational approach is anticipated to facilitate large-scale identification of TPPs and accelerate their applications in the food industry. The codes and datasets are available at https://github.com/plenoi/SAPPHIRE.

Glycosyl hydrolase 11 (xynA) gene with xylanase activity from thermophilic bacteria isolated from thermal springs

Background

Hemicellulose is one of the copious polymer in lignocellulosic biomass (LCB). It is primarily composed of xylan linked by β-1,4 glycosidic bonds. Xylanase preferentially cleaves the β-1,4-glycosidic bonds in the xylan backbone resulting in complete hydrolysis of the biomass. Thermostable variants of glycoside hydrolases act as robust catalysts, not only in degradation but also during processing, to obtain specific carbohydrate-containing chemicals and materials (Ramasamy et al. in Madras Agric J 107(special):1. 10.29321/MAJ.2020.000382, 2020).

Results

The xylanase production by two thermophilic bacteria isolated from thermal springs was evaluated. In addition, the gene encoding this industrially vital enzyme was isolated and characterized, and its protein structure was analyzed. The thermophilic bacteria producing xylanases were isolated from augmented sawdust and banana fiber biomass from hot springs of Himachal Pradesh and identified as Bacillus subtilis VSDB5 and Bacillus licheniformis KBFB4 using 16S rRNA gene sequencing. The persistent xylanase activity revealed that the enzyme is secreted extracellularly with the maximum activity of 0.76 IU mL−1 and 1.0 IU mL−1 at 6 h and 12 h of growth by KBFB4 and VSDB5, respectively, under submerged fermentation. Both the strains exhibited the maximum activity at pH 6 and a temperature of 50 °C. The xylanases of KBFB4 and VSDB5 were thermostable and retained 40% of their activity at 60 °C after incubation for 30 min. Xylanase of VSDB5 had wide thermotolerance and retained 20% of its activity from 60 to 80 °C, whereas xylanase of KBFB4 showed wide alkali tolerance and retained 80% of its activity until pH 10. The xylanase (xynA)-encoding gene (650 bp) cloned from both the strains using specific primers showed 98 to 99% homology to β-1,4-endoxylanase gene. Further in silico analysis predicted that the xylanase protein, with a molecular weight of 23 kDa, had a high pI (9.44–9.65), which explained the alkaline nature of the enzyme and greater aliphatic index (56.29). This finding suggested that the protein is thermostable. Multiple sequence alignment and homology modeling of the protein sequence revealed that the gene product belonged to the GH11 family, indicating its possible application in bioconversion.

Conclusion

The strains B. subtilis VSDB5 and B. licheniformis KBFB4 obtained from hot springs of Himachal Pradesh produced potent and alkali-tolerant thermostable xylanases, which belong to the GH11 family. The enzyme can be supplemented in industrial applications for biomass conversion at high temperatures and pH (or in processes involving alkali treatment).

Graphical Abstract

iThermo: A Sequence-Based Model for Identifying Thermophilic Proteins Using a Multi-Feature Fusion Strategy

Thermophilic proteins have important application value in biotechnology and industrial processes. The correct identification of thermophilic proteins provides important information for the application of these proteins in engineering. The identification method of thermophilic proteins based on biochemistry is laborious, time-consuming, and high cost. Therefore, there is an urgent need for a fast and accurate method to identify thermophilic proteins. Considering this urgency, we constructed a reliable benchmark dataset containing 1,368 thermophilic and 1,443 non-thermophilic proteins. A multi-layer perceptron (MLP) model based on a multi-feature fusion strategy was proposed to discriminate thermophilic proteins from non-thermophilic proteins. On independent data set, the proposed model could achieve an accuracy of 96.26%, which demonstrates that the model has a good application prospect. In order to use the model conveniently, a user-friendly software package called iThermo was established and can be freely accessed at http://lin-group.cn/server/iThermo/index.html. The high accuracy of the model and the practicability of the developed software package indicate that this study can accelerate the discovery and engineering application of thermally stable proteins.

Novel Proteome and N-Glycoproteome of the Thermophilic Fungus Chaetomium thermophilum in Response to High Temperature

Thermophilic fungi are eukaryotic species that grow at high temperatures, but little is known about the underlying basis of thermophily at cell and molecular levels. Here the proteome and N-glycoproteome of Chaetomium thermophilum at varying culture temperatures (30, 50, and 55°C) were studied using hydrophilic interaction liquid chromatography enrichment and high-resolution liquid chromatography–tandem mass spectroscopy analysis. With respect to the proteome, the numbers of differentially expressed proteins were 1,274, 1,374, and 1,063 in T50/T30, T55/T30, and T55/T50, respectively. The upregulated proteins were involved in biological processes, such as protein folding and carbohydrate metabolism. Most downregulated proteins were involved in molecular functions, including structural constituents of the ribosome and other protein complexes. For the N-glycoproteome, the numbers of differentially expressed N-glycoproteins were 160, 176, and 128 in T50/T30, T55/T30, and T55/T50, respectively. The differential glycoproteins were mainly involved in various types of N-glycan biosynthesis, mRNA surveillance pathway, and protein processing in the endoplasmic reticulum. These results indicated that an efficient protein homeostasis pathway plays an essential role in the thermophily of C. thermophilum, and N-glycosylation is involved by affecting related proteins. This is the novel study to reveal thermophilic fungi’s physiological response to high-temperature adaptation using omics analysis, facilitating the exploration of the thermophily mechanism of thermophilic fungi.

Dynamic cross correlation analysis of Thermus thermophilus alkaline phosphatase and determinants of thermostability

BACKGROUND

Understanding the determinants of protein thermostability is very important both from the theoretical and applied perspective. One emerging view in thermostable enzymes seems to indicate that a salt bridge/charged residue network plays a fundamental role in their thermostability.

METHODS

The structure of alkaline phosphatase (AP) from Thermus thermophilus HB8 was solved by X-ray crystallography at 2.1 Å resolution. The obtained structure was further analyzed by molecular dynamics studies at different temperatures (303 K, 333 K and 363 K) and compared to homologous proteins from the cold-adapted organisms Shewanella sp. and Vibrio strain G15-21. To analyze differences in measures of dynamic variation, several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis were used. Using hierarchical clustering, the obtained results were combined to determine residues showing high degree dynamical variations due to temperature jumps. Furthermore, dynamic cross correlation (DCC) analysis was carried out to characterize networks of charged residues.

RESULTS

Top clustered residues showed a higher propensity for thermostabilizing mutations, indicating evolutionary pressure acting on thermophilic organisms. The description of rotamer distributions by Gini coefficients and Kullback-Leibler (KL) divergence both revealed significant correlations with temperature. DCC analysis revealed a significant trend to de-correlation of the movement of charged residues at higher temperatures.

SIGNIFICANCE

The de-correlation of charged residues detected in Thermus thermophilus AP, highlights the importance of dynamic electrostatic network interactions for the thermostability of this enzyme.

Crystal structure of a Cu,Zn superoxide dismutase from the thermophilic fungus Chaetomium thermophilum

BACKGROUND

Thermophilic fungi have recently emerged as a promising source of thermostable enzymes. Superoxide dismutases are key antioxidant metalloenzymes with promising therapeutic effects in various diseases, both acute and chronic. However, structural heterogeneity and low thermostability limit their therapeutic efficacy.

OBJECTIVE

Although several studies from hypethermophilic superoxide dismutases (SODs) have been reported, information about Cu,Zn-SODs from thermophilic fungi is scarce. Chaetomium thermophilum is a thermophilic fungus that could provide proteins with thermophilic properties.

METHOD

The enzyme was expressed in Pichia pastoris cells and crystallized using the vapor-diffusion method. X-ray data were collected, and the structure was determined and refined to 1.56 Å resolution. Structural analysis and comparisons were carried out.

RESULTS

The presence of 8 molecules (A through H) in the asymmetric unit resulted in four different interfaces. Molecules A and F form the typical homodimer which is also found in other Cu,Zn-SODs. Zinc was present in all subunits of the structure while copper was found in only four subunits with reduced occupancy (C, D, E and F).

CONCLUSION

The ability of the enzyme to form oligomers and the elevated Thr:Ser ratio may be contributing factors to its thermal stability. Two hydrophobic residues that participate in interface formation and are not present in other CuZn-SODs may play a role in the formation of new interfaces and the oligomerization process. The CtSOD crystal structure reported here is the first Cu,Zn-SOD structure from a thermophilic fungus.

An In-Silico Comparative Study of Lipases from the Antarctic Psychrophilic Ciliate Euplotes focardii and the Mesophilic Congeneric Species Euplotes crassus: Insight into Molecular Cold-Adaptation

Cold-adapted enzymes produced by psychrophilic organisms have elevated catalytic activities at low temperatures compared to their mesophilic counterparts. This is largely due to amino acids changes in the protein sequence that often confer increased molecular flexibility in the cold. Comparison of structural changes between psychrophilic and mesophilic enzymes often reveal molecular cold adaptation. In the present study, we performed an in-silico comparative analysis of 104 hydrolytic enzymes belonging to the family of lipases from two evolutionary close marine ciliate species: The Antarctic psychrophilic Euplotes focardii and the mesophilic Euplotes crassus. By applying bioinformatics approaches, we compared amino acid composition and predicted secondary and tertiary structures of these lipases to extract relevant information relative to cold adaptation. Our results not only confirm the importance of several previous recognized amino acid substitutions for cold adaptation, as the preference for small amino acid, but also identify some new factors correlated with the secondary structure possibly responsible for enhanced enzyme activity at low temperatures. This study emphasizes the subtle sequence and structural modifications that may help to transform mesophilic into psychrophilic enzymes for industrial applications by protein engineering.

Evolution of Protein Structure and Stability in Global Warming

This review focuses on the molecular signatures of protein structures in relation to evolution and survival in global warming. It is based on the premise that the power of evolutionary selection may lead to thermotolerant organisms that will repopulate the planet and continue life in general, but perhaps with different kinds of flora and fauna. Our focus is on molecular mechanisms, whereby known examples of thermoresistance and their physicochemical characteristics were noted. A comparison of interactions of diverse residues in proteins from thermophilic and mesophilic organisms, as well as reverse genetic studies, revealed a set of imprecise molecular signatures that pointed to major roles of hydrophobicity, solvent accessibility, disulfide bonds, hydrogen bonds, ionic and π-electron interactions, and an overall condensed packing of the higher-order structure, especially in the hydrophobic regions. Regardless of mutations, specialized protein chaperones may play a cardinal role. In evolutionary terms, thermoresistance to global warming will likely occur in stepwise mutational changes, conforming to the molecular signatures, such that each "intermediate" fits a temporary niche through punctuated equilibrium, while maintaining protein functionality. Finally, the population response of different species to global warming may vary substantially, and, as such, some may evolve while others will undergo catastrophic mass extinction.

Discrimination of Thermophilic Proteins and Non-thermophilic Proteins Using Feature Dimension Reduction

Thermophilicity is a very important property of proteins, as it sometimes determines denaturation and cell death. Thus, methods for predicting thermophilic proteins and non-thermophilic proteins are of interest and can contribute to the design and engineering of proteins. In this article, we describe the use of feature dimension reduction technology and LIBSVM to identify thermophilic proteins. The highest accuracy obtained by cross-validation was 96.02% with 119 parameters. When using only 16 features, we obtained an accuracy of 93.33%. We discuss the importance of the different characteristics in identification and report a comparison of the performance of support vector machine to that of other methods.

Thermophilic Chitinases: Structural, Functional and Engineering Attributes for Industrial Applications

Chitin is the second most widely found natural polymer next to cellulose. Chitinases degrade the insoluble chitin to bioactive chitooligomers and monomers for various industrial applications. Based on their function, these enzymes act as biocontrol agents against pathogenic fungi and invasive pests compared with conventional chemical fungicides and insecticides. They have other functional roles in shellfish waste management, fungal protoplast generation, and Single-Cell Protein production. Among the chitinases, thermophilic and thermostable chitinases are gaining popularity in recent years, as they can withstand high temperatures and maintain the enzyme stability for longer periods. Not all chitinases are thermostable; hence, tailor-made thermophilic chitinases are designed to enhance their thermostability by direct evolution, genetic engineering involving mutagenesis, and proteomics approach. Although research has been done extensively on cloning and expression of thermophilic chitinase genes, there are only few papers discussing on the mechanism of chitin degradation using thermophiles. The current review discusses the sources of thermophilic chitinases, improvement of protein stability by gene manipulation, metagenomics approaches, chitin degradation mechanism in thermophiles, and their prospective applications for industrial, agricultural, and pharmaceutical purposes.

Insights on protein thermal stability: a graph representation of molecular interactions

Motivation

Understanding the molecular mechanisms of thermal stability is a challenge in protein biology. Indeed, knowing the temperature at which proteins are stable has important theoretical implications, which are intimately linked with properties of the native fold, and a wide range of potential applications from drug design to the optimization of enzyme activity.

Results

Here, we present a novel graph-theoretical framework to assess thermal stability based on the structure without any a priori information. In this approach we describe proteins as energy-weighted graphs and compare them using ensembles of interaction networks. Investigating the position of specific interactions within the 3D native structure, we developed a parameter-free network descriptor that permits to distinguish thermostable and mesostable proteins with an accuracy of 76% and area under the receiver operating characteristic curve of 78%.

Availability and implementation

Code is available upon request to edoardo.milanetti@uniroma1.it

Supplementary information

Supplementary data are available at Bioinformatics online.

Crystal structure of Thermus thermophilus methylenetetrahydrofolate dehydrogenase and determinants of thermostability

The elucidation of mechanisms behind the thermostability of proteins is extremely important both from the theoretical and applied perspective. Here we report the crystal structure of methylenetetrahydrofolate dehydrogenase (MTHFD) from Thermus thermophilus HB8, a thermophilic model organism. Molecular dynamics trajectory analysis of this protein at different temperatures (303 K, 333 K and 363 K) was compared with homologous proteins from the less temperature resistant organism Thermoplasma acidophilum and the mesophilic organism Acinetobacter baumannii using several data reduction techniques like principal component analysis (PCA), residue interaction network (RIN) analysis and rotamer analysis. These methods enabled the determination of important residues for the thermostability of this enzyme. The description of rotamer distributions by Gini coefficients and Kullback–Leibler (KL) divergence both revealed significant correlations with temperature. The emerging view seems to indicate that a static salt bridge/charged residue network plays a fundamental role in the temperature resistance of Thermus thermophilus MTHFD by enhancing both electrostatic interactions and entropic energy dispersion. Furthermore, this analysis uncovered a relationship between residue mutations and evolutionary pressure acting on thermophilic organisms and thus could be of use for the design of future thermostable enzymes.

Structural insights on binding mechanism of CAD complexes (CPSase, ATCase and DHOase)

Pyrimidine biosynthetic pathway enzymes constitute an important target for the development of antitumor drugs. To understand the role of binding mechanisms underlying the inborn errors of pyrimidine biosynthetic pathway, structure and function of enzymes have been analyzed. Pyrimidine biosynthetic pathway is initiated by CAD enzymes that harbor the first three enzymatic activities facilitated by Carbamoyl Phosphate Synthetase (CPSase), Aspartate Transcarbamoylase (ATCase) and Dihydroorotase (DHOase). While being an attractive therapeutic target, the lack of data driven us to study the CPSase (CarA and CarB) and its mode of binding to ATCase and DHOase which are the major limitation for its structural optimization. Understanding the binding mode of CPSase, ATCase and DHOase could help to identify the potential interface hotspot residues that favor the mechanism behind it. The mechanistic insight into the CAD complexes were achieved through Molecular modeling, Protein-Protein docking, Alanine scanning and Molecular dynamics (MD) Studies. The hotspot residues present in the CarB region of carboxy phosphate and carbamoyl phosphate synthetic domains are responsible for the assembly of CAD (CPSase-ATCase-DHOase) complexes. Overall analysis suggests that the identified hotspot residues were confirmed by alanine scanning and important for the regulation of pyrimidine biosynthesis. MD simulations analysis provided the prolonged stability of the interacting complexes. The present study reveals the novel hotspot residues such as Glu134, Glu147, Glu154, Asp266, Lys269, Glu274, Asp333, Trp459, Asp526, Asp528, Glu533, Glu544, Glu546, Glu800, Val855, Asp877, Tyr884 and Gln919 which could be targeted for structure-based inhibitor design to potentiate the CAD mediated regulation of aggressive tumors.Communicated by Ramaswamy H. Sarma.

Lysosomal protein thermal stability does not correlate with cellular half-life: global observations and a case study of tripeptidyl-peptidase 1

Late-infantile neuronal ceroid lipofuscinosis (LINCL) is a neurodegenerative lysosomal storage disorder caused by mutations in the gene encoding the protease tripeptidyl-peptidase 1 (TPP1). Progression of LINCL can be slowed or halted by enzyme replacement therapy, where recombinant human TPP1 is administered to patients. In this study, we utilized protein engineering techniques to increase the stability of recombinant TPP1 with the rationale that this may lengthen its lysosomal half-life, potentially increasing the potency of the therapeutic protein. Utilizing multiple structure-based methods that have been shown to increase the stability of other proteins, we have generated and evaluated over 70 TPP1 variants. The most effective mutation, R465G, increased the melting temperature of TPP1 from 55.6°C to 64.4°C and increased its enzymatic half-life at 60°C from 5.4 min to 21.9 min. However, the intracellular half-life of R465G and all other variants tested in cultured LINCL patient-derived lymphoblasts was similar to that of WT TPP1. These results provide structure/function insights into TPP1 and indicate that improving in vitro thermal stability alone is insufficient to generate TPP1 variants with improved physiological stability. This conclusion is supported by a proteome-wide analysis that indicates that lysosomal proteins have higher melting temperatures but also higher turnover rates than proteins of other organelles. These results have implications for similar efforts where protein engineering approaches, which are frequently evaluated in vitro, may be considered for improving the physiological properties of proteins, particularly those that function in the lysosomal environment.

Evaluation of Luciferase Thermal Stability by Arginine Saturation in the Flexible Loops

Background:

The firefly luciferase enzyme is widely used in protein engineering and diverse areas of biotechnology, but the main problem with this enzyme is low-temperature stability. Previous reports indicated that surface areas of thermostable proteins are rich in arginine, which increased their thermal stability. In this study, this aspect of thermophilic proteins evaluated by mutations of surface residues to Arg. Here, we report the construction, purification, and studying of these mutated luciferases.

Methods:

For mutagenesis, the QuikChange site-directed mutagenesis was used and the I108R, T156R, and N177R mutant luciferases were created. In the following, the expression and purification of wild-type and mutant luciferases were conducted and their kinetic and structural properties were analyzed. To analyze the role of these Arg in these loops, the 3D models of these mutants’ enzymes were constructed in the I-TASSER server and the exact situation of these mutants was studied by the SPDBV and PyMOL software.

Results:

Overall, the optimum temperature of these mutated enzymes was not changed. However, after 30 min incubation of these mutated enzymes at 30°C, the I108R, T156R, N177R, and wild-type kept the 80%, 50%, 20%, and 20% of their original activity, respectively. It should be noted that substitution of these residues by Arg preserved the specific activity of firefly luciferase.

Conclusion:

Based on these results, it can be concluded that T156R and N177R mutants by compacting local protein structure, increased the thermostability of luciferase. However, insertion of positively charged residues in these positions create the new hydrogen bonds that associated with a series of structural changes and confirmed by intrinsic and extrinsic fluorescence spectroscopy and homology modeling studies.

Confirmation of the formation of salt bridges in the denatured state of CutA1 protein using molecular dynamics simulations

It remains unclear how the abundant charged residues in proteins from hyperthermophiles contribute to the stabilization of proteins. Previously, based on molecular dynamics (MD) simulations, we proposed that these charged residues decrease the entropic effect by forming salt bridges in the denatured state under physiological conditions (Yutani et al., Sci. Rep. 8, 7613 (2018)). Because the quality of MD results is strongly dependent on the force fields used, in this study we performed the MD simulations using a different force field (AMBER99SB) along with the one we used before (Gromos43a1), at the same temperatures examined previously as well as at higher temperatures. In these experiments, we used the same ionic mutant (Ec0VV6) of CutA1 from Escherichia coli as in the previous study. In MD simulations at 300 K, Lys87 and Arg88 in the loop region of Ec0VV6 formed salt bridges with different favorable pairs in different force fields. Furthermore, the helical content and radius of gyration differed slightly between two force fields. However, at a higher temperature (600 K), the average numbers of salt bridges for the six substituted residues of Ec0VV6 were 0.87 per residue for Gromos43a1 and 0.88 for AMBER99SB in 400-ns MD simulation, indicating that the values were similar despite the use of different force fields. These observations suggest that the charged residues in Ec0VV6 can form a considerable number of salt bridges, even in the denatured state with drastic fluctuation at 600 K. These results corroborate our previous proposal.

Crystal Structure of a GH3 β-Glucosidase from the Thermophilic Fungus Chaetomium thermophilum

Beta-glucosidases (β-glucosidases) have attracted considerable attention in recent years for use in various biotechnological applications. They are also essential enzymes for lignocellulose degradation in biofuel production. However, cost-effective biomass conversion requires the use of highly efficient enzymes. Thus, the search for new enzymes as better alternatives of the currently available enzyme preparations is highly important. Thermophilic fungi are nowadays considered as a promising source of enzymes with improved stability. Here, the crystal structure of a family GH3 β-glucosidase from the thermophilic fungus Chaetomium thermophilum (CtBGL) was determined at a resolution of 2.99 Å. The structure showed the three-domain architecture found in other β-glucosidases with variations in loops and linker regions. The active site catalytic residues in CtBGL were identified as Asp287 (nucleophile) and Glu517 (acid/base). Structural comparison of CtBGL with Protein Data Bank (PDB)-deposited structures revealed variations among glycosylated Asn residues. The enzyme displayed moderate glycosylation compared to other GH3 family β-glucosidases with similar structure. A new glycosylation site at position Asn504 was identified in CtBGL. Moreover, comparison with respect to several thermostability parameters suggested that glycosylation and charged residues involved in electrostatic interactions may contribute to the stability of the enzyme at elevated temperatures. The reported CtBGL structure provides additional insights into the family GH3 enzymes and could offer new ideas for further improvements in β-glucosidases for more efficient use in biotechnological applications regarding cellulose degradation.

Selection and analyses of variants of a designed protein suggest importance of hydrophobicity of partially buried sidechains for protein stability at high temperatures

Computationally designed proteins of high stability provide specimen in addition to natural proteins for the study of sequence-structure stability relationships at the very high end of protein stability spectrum. The melting temperature of E_1r26, a protein we previously designed using the A Backbone-based Amino aCid Usage Survey (ABACUS) sequence design program, is above 110 °C, more than 50 °C higher than that of the natural thioredoxin protein whose backbone (PDB ID 1R26) has been used as the design target. Using an experimental selection approach, we obtained variants of E_1r26 that remain folded but are of reduced stability, including one whose unfolding temperature and denaturing guanidine concentration are similar to those of 1r26. The mutant unfolds with a certain degree of cooperativity. Its structure solved by X-ray crystallography agrees with that of 1r26 by a root mean square deviation of 1.3 Å, adding supports to the accuracy of the ABACUS method. Analyses of intermediate mutants indicate that the substitution of two partially buried hydrophobic residues (isoleucine and leucine) by polar residues (threonine and serine, respectively) are responsible for the dramatic change in the unfolding temperature. It is suggested that the effects of mutations located in rigid secondary structure regions, but not those in loops, may be well predicted through ABACUS mutation energy analysis. The results also suggest that hydrophobic effects involving intermediately buried sidechains can be critically important for protein stability at high temperatures.

Biochemical profiles of two thermostable and organic solvent-tolerant esterases derived from a compost metagenome

Owing to the functional versatility and potential applications in industry, interest in lipolytic enzymes tolerant to organic solvents is increasing. In this study, functional screening of a compost soil metagenome resulted in identification of two lipolytic genes, est1 and est2, encoding 270 and 389 amino acids, respectively. The two genes were heterologously expressed and characterized. Est1 and Est2 are thermostable enzymes with optimal enzyme activities at 80 and 70 °C, respectively. A second-order rotatable design, which allows establishing the relationship between multiple variables with the obtained responses, was used to explore the combined effects of temperature and pH on esterase stability. The response curve indicated that Est1, and particularly Est2, retained high stability within a broad range of temperature and pH values. Furthermore, the effects of organic solvents on Est1 and Est2 activities and stabilities were assessed. Notably, Est2 activity was significantly enhanced (two- to tenfold) in the presence of ethanol, methanol, isopropanol, and 1-propanol over a concentration range between 6 and 30% (v/v). For the short-term stability (2 h of incubation), Est2 exhibited high tolerance against 60% (v/v) of ethanol, methanol, isopropanol, DMSO, and acetone, while Est1 activity resisted these solvents only at lower concentrations (below 30%, v/v). Est2 also displayed high stability towards some water-immiscible organic solvents, such as ethyl acetate, diethyl ether, and toluene. With respect to long-term stability, Est2 retained most of its activity after 26 days of incubation in the presence of 30% (v/v) ethanol, methanol, isopropanol, DMSO, or acetone. All of these features indicate that Est1 and Est2 possess application potential.

Taxonomic Landscape of the Dark Proteomes: Whole-Proteome Scale Interplay Between Structural Darkness, Intrinsic Disorder, and Crystallization Propensity

Growth rate of the protein sequence universe dramatically exceeds the speed of expansion for the protein structure universe, generating an immense dark proteome that includes proteins with unknown structure. A whole-proteome scale analysis of 5.4 million proteins from 987 proteomes in the three domains of life and viruses to systematically dissect an interplay between structural coverage, degree of putative intrinsic disorder, and predicted propensity for structure determination is performed. It has been found that Archaean and Bacterial proteomes have relatively high structural coverage and low amounts of disorder, whereas Eukaryotic and Viral proteomes are characterized by a broad spread of structural coverage and higher disorder levels. The analysis reveals that dark proteomes (i.e., proteomes containing high fractions of proteins with unknown structure) have significantly elevated amounts of intrinsic disorder and are predicted to be difficult to solve structurally. Although the majority of dark proteomes are of viral origin, many dark viral proteomes have at least modest crystallization propensity and only a handful of them are enriched in the intrinsic disorder. The disorder, structural coverage, and propensity are mapped for structural determination onto a novel proteome-level sequence similarity network to analyze the interplay of these characteristics in the taxonomic landscape.

RankProt: A multi criteria-ranking platform to attain protein thermostabilizing mutations and its in vitro applications - Attribute based prediction method on the principles of Analytical Hierarchical Process

Attaining recombinant thermostable proteins is still a challenge for protein engineering. The complexity is the length of time and enormous efforts required to achieve the desired results. Present work proposes a novel and economic strategy of attaining protein thermostability by predicting site-specific mutations at the shortest possible time. The success of the approach can be attributed to Analytical Hierarchical Process and the outcome was a rationalized thermostable mutation(s) prediction tool- RankProt. Briefly the method involved ranking of 17 biophysical protein features as class predictors, derived from 127 pairs of thermostable and mesostable proteins. Among the 17 predictors, ionic interactions and main-chain to main-chain hydrogen bonds were the highest ranked features with eigen value of 0.091. The success of the tool was judged by multi-fold in silico validation tests and it achieved the prediction accuracy of 91% with AUC 0.927. Further, in vitro validation was carried out by predicting thermostabilizing mutations for mesostable Bacillus subtilis lipase and performing the predicted mutations by multi-site directed mutagenesis. The rationalized method was successful to render the lipase thermostable with optimum temperature stability and T_m increase by 20°C and 7°C respectively. Conclusively it can be said that it was the minimum number of mutations in comparison to the number of mutations incorporated to render Bacillus subtilis lipase thermostable, by directed evolution techniques. The present work shows that protein stabilizing mutations can be rationally designed by balancing the biophysical pleiotropy of proteins, in accordance to the selection pressure.

Ion-ion interactions in the denatured state contribute to the stabilization of CutA1 proteins

In order to elucidate features of the denatured state ensembles that exist in equilibrium with the native state under physiological conditions, we performed 1.4-μs molecular dynamics (MD) simulations at 400 K and 450 K using the monomer subunits of three CutA1 mutants from Escherichia coli: an SH-free mutant (Ec0SH) with denaturation temperature (T_d) = 85.6 °C, a hydrophobic mutant (Ec0VV) with T_d = 113.3 °C, and an ionic mutant (Ec0VV_6) with T_d = 136.8 °C. The occupancy of salt bridges by the six substituted charged residues in Ec0VV_6 was 140.1% at 300 K and 89.5% at 450 K, indicating that even in the denatured state, salt bridge occupancy was high, approximately 60% of that at 300 K. From these results, we can infer that proteins from hyperthermophiles with a high ratio of charged residues are stabilized by a decrease in conformational entropy due to ion–ion interactions in the denatured state. The mechanism must be comparable to the stabilization conferred by disulfide bonds within a protein. This suggests that introduction of charged residues, to promote formation of salt bridges in the denatured state, would be a simple way to rationally design stability-enhanced mutants.

Strategy for Stabilization of CutA1 Proteins Due to Ion-Ion Interactions at Temperatures of over 100 °C

In order to elucidate the contribution of charged residues to protein stabilization at temperatures of over 100 °C, we constructed many mutants of the CutA1 protein ( EcCutA1) from Escherichia coli. The goal was to see if one can achieve the same stability as for a CutA1 from hyperthermophile Pyrococcus horikoshii that has the denaturation temperature near 150 °C. The hydrophobic mutant of EcCutA1 ( Ec0VV) with denaturation temperature ( T_d) of 113.2 °C was used as a template for mutations. The highest T_d of Ec0VV mutants substituted by a single charged residue was 118.4 °C. Multiple ion mutants were also constructed by combination of single mutants and found to have an increased thermostability. The highest stability of multiple mutants was a mutant substituted by nine charged residues that had a T_d of 142.2 °C. To evaluate the energy of ion-ion interactions of mutant proteins, we used the structural ensemble obtained by a molecular dynamics simulation at 300 K. The T_d of ionic mutants linearly increases with the increments of the computed energy of ion-ion interactions for ionic mutant proteins even up to the temperatures near 140 °C, suggesting that ion-ion interactions cumulatively contribute to the stabilization of a protein at high temperatures.

Structural and Dynamics Comparison of Thermostability in Ancient, Modern, and Consensus Elongation Factor Tus

Rationally engineering thermostability in proteins would create enzymes and receptors that function under harsh industrial applications. Several sequence-based approaches can generate thermostable variants of mesophilic proteins. To gain insight into the mechanisms by which proteins become more stable, we use structural and dynamic analyses to compare two popular approaches, ancestral sequence reconstruction (ASR) and the consensus method, used to generate thermostable variants of Elongation Factor Thermo-unstable (EF-Tu). We present crystal structures of ancestral and consensus EF-Tus, accompanied by molecular dynamics simulations aimed at probing the strategies employed to enhance thermostability. All proteins adopt crystal structures similar to extant EF-Tus, revealing no difference in average structure between the methods. Molecular dynamics reveals that ASR-generated sequences retain dynamic properties similar to extant, thermostable EF-Tu from Thermus aquaticus, while consensus EF-Tu dynamics differ from evolution-based sequences. This work highlights the advantage of ASR for engineering thermostability while preserving natural motions in multidomain proteins.

Density discriminates between thermophilic and mesophilic proteins

Despite an intense interest and a remarkable number of studies on the subject, the relationships between thermostability and (primary, secondary and tertiary) structure of proteins are still not fully understood. Here, comparing the protein density - defined by the ratio between the residue number and protein excluded volume - for a set of thermophilic/mesophilic pairs, we provide evidence that this property is connected to the optimal growth temperature. In particular, our results indicate that thermophilic proteins have - in general - a lower density with respect to the mesophilic counterparts, being such a correlation more pronounced for optimal growth temperature differences greater than 40°C. The effect of the protein thermostability changes on the molecular shape is also presented.

Role of an N-terminal extension in stability and catalytic activity of a hyperthermostable α/β hydrolase fold esterase

The carbohydrate esterase family 7 (CE7) enzymes catalyze the deacetylation of acetyl esters of a broad range of alcohols and is unique in its activity towards cephalosporin C. The CE7 fold contains a conserved N-terminal extension that distinguishes it from the canonical α/β hydrolase fold. The hexameric quaternary structure indicates that the N-terminus may affect activity and specificity by controlling access of substrates to the buried active sites via an entrance tunnel. In this context, we characterized the catalytic parameters, conformation and thermal stability of two truncation variants lacking four and ten residues of the N-terminal region of the hyperthermostable Thermotoga maritima CE7 acetyl esterase (TmAcE). The truncations did not affect the secondary structure or the fold but modulated the oligomerization dynamics. A modest increase was observed in substrate specificity for acetylated xylose compared with acetylated glucose. A drastic reduction of ~30-40°C in the optimum temperature for activity of the variants indicated lower thermal stability. The loss of hyperthermostability appears to be an indirect effect associated with an increase in the conformational flexibility of an otherwise rigid neighboring loop containing a catalytic triad residue. The results suggest that the N-terminal extension was evolutionarily selected to preserve the stability of the enzyme.

Structural analysis and interaction studies of acyl-carrier protein (acpP) of Staphylococcus aureus, an extraordinarily thermally stable protein

Acyl-carrier-protein (acpP) is an essential protein in fatty acid biosynthesis of Staphylococcus aureus [Cronan, J.E. and Thomas, J. (2009). Complex enzymes in microbial natural product biosynthesis, part B: polyketides, aminocoumarins and carbohydrates.

METHOD

Enzymol. 459, 395-433; Halavaty, A.S., Kim, Y., Minasov, G., Shuvalova, L., Dubrovska, I., Winsor, J., Zhou, M., Onopriyenko, O., Skarina, T., Papazisi, L., et al. (2012). Structural characterization and comparison of three acyl-carrier-protein synthases from pathogenic bacteria. Acta Crystallogr. Sect. D Biol. Crystallogr. 68, 1359-1370]. The inactive apo-form is converted to the active holo-enzyme by acyl-carrier protein synthase (acpS) through addition of a 4'-phosphopantetheine group from coenzyme A to a conserved serine residue of acpP [Flugel, R.S., Hwangbo, Y., Lambalot, R.H., Cronan, J.E., and Walsh, C.T. (2000). Holo-(acyl-carrier protein) synthase and phosphopantetheinyl transfer in Escherichia coli. J. Biol. Chem. 275, 959-968; Lambalot, R.H. and Walsh, C.T. (1995). Cloning, overproduction, and characterization of the Escherichia coli holo-acyl-carrier protein synthase. J. Biol. Chem. 270, 24658-24661]. Once activated, acpP acts as an anchor for the growing fatty acid chain. Structural data from X-ray crystallographic analysis reveals that, despite its small size (8 kDa), acpP adopts a distinct, mostly α-helical structure when complexed with acpS [Halavaty, A.S., Kim, Y., Minasov, G., Shuvalova, L., Dubrovska, I., Winsor, J., Zhou, M., Onopriyenko, O., Skarina, T., Papazisi, L., et al. (2012). Structural characterization and comparison of three acyl-carrier-protein synthases from pathogenic bacteria. Acta Crystallogr. Sect. D Biol. Crystallogr. 68, 1359-1370; Byers, D.M. and Gong, H. (2007). Acyl carrier protein: structure-function relationships in a conserved multifunctional protein family. Biochem. Cell Biol. 85, 649-662]. We expressed and purified recombinant, active S. aureus acpP from Escherichia coli and mimicked the beginning of fatty acid biosynthesis by employing an [14C]-acp loading assay. Surprisingly, acpP remained functional even after heat treatment at 95°C for up to 10 min. NMR data from 2D-HSQC experiments as well as interaction studies with acpS confirmed that acpP is structured and active both before and after heat treatment, with no significant differences between the two. Thus, our data suggest that S. aureus acpP is a highly stable protein capable of maintaining its structure at high temperatures.

Extremozymes from metagenome: Potential applications in food processing

The long-established use of enzymes for food processing and product formulation has resulted in an increased enzyme market compounding to 7.0% annual growth rate. Advancements in molecular biology and recognition that enzymes with specific properties have application for industrial production of infant, baby and functional foods boosted research toward sourcing the genes of microorganisms for enzymes with distinctive properties. In this regard, functional metagenomics for extremozymes has gained attention on the premise that such enzymes can catalyze specific reactions. Hence, metagenomics that can isolate functional genes of unculturable extremophilic microorganisms has expanded attention as a promising tool. Developments in this field of research in relation to food sector are reviewed.

Exploring Covalent Allosteric Inhibition of Antigen 85C from Mycobacterium tuberculosis by Ebselen Derivatives

Previous studies identified ebselen as a potent in vitro and in vivo inhibitor of the Mycobacterium tuberculosis (Mtb) antigen 85 (Ag85) complex, comprising three homologous enzymes required for the biosynthesis of the mycobacterial cell wall. In this study, the Mtb Ag85C enzyme was cocrystallized with azido and adamantyl ebselen derivatives, resulting in two crystallographic structures of 2.01 and 1.30 Å resolution, respectively. Both structures displayed the anticipated covalent modification of the solvent accessible, noncatalytic Cys209 residue forming a selenenylsulfide bond. Continuous difference density for both thiol modifiers allowed for the assessment of interactions that influence ebselen binding and inhibitor orientation that were unobserved in previous Ag85C ebselen structures. The k_inact/K_I values for ebselen, adamantyl ebselen, and azido ebselen support the importance of observed constructive chemical interactions with Arg239 for increased in vitro efficacy toward Ag85C. To better understand the in vitro kinetic properties of these ebselen derivatives, the energetics of specific protein-inhibitor interactions and relative reaction free energies were calculated for ebselen and both derivatives using density functional theory. These studies further support the different in vitro properties of ebselen and two select ebselen derivatives from our previously published ebselen library with respect to kinetics and protein-inhibitor interactions. In both structures, the α9 helix was displaced farther from the enzyme active site than the previous Ag85C ebselen structure, resulting in the restructuring of a connecting loop and imparting a conformational change to residues believed to play a role in substrate binding specific to Ag85C. These notable structural changes directly affect protein stability, reducing the overall melting temperature by up to 14.5 °C, resulting in the unfolding of protein at physiological temperatures. Additionally, this structural rearrangement due to covalent allosteric modification creates a sizable solvent network that encompasses the active site and extends to the modified Cys209 residue. In all, this study outlines factors that influence enzyme inhibition by ebselen and its derivatives while further highlighting the effects of the covalent modification of Cys209 by said inhibitors on the structure and stability of Ag85C. Furthermore, the results suggest a strategy for developing new classes of Ag85 inhibitors with increased specificity and potency.

The Missing Electrostatic Interactions Between DNA Substrate and Sulfolobus solfataricus DNA Photolyase: What is the Role of Charged Amino Acids in Thermophilic DNA Binding Proteins?

DNA photolyase can be used to study how a protein with its required cofactor has adapted over a large temperature range. The enzymatic activity and thermodynamics of substrate binding for protein from Sulfolobus solfataricus were directly compared to protein from Escherichia coli. Turnover numbers and catalytic activity were virtually identical, but organic cosolvents may be necessary to maintain activity of the thermophilic protein at higher temperatures. UV-damaged DNA binding to the thermophilic protein is less favorable by ∼2 kJ/mol. The enthalpy of binding is ∼10 kJ/mol less exothermic for the thermophile, but the amount and type of surface area buried upon DNA binding appears to be somewhat similar. The most important finding was observed when ionic strength studies were used to separate binding interactions into electrostatic and nonelectrostatic contributions; DNA binding to the thermophilic protein appears to lack the electrostatic contributions observed with the mesophilic protein.